The Combined Effect of Temperature and Host Clonal Line on the Microbiota of a Planktonic Crustacean.

Host-associated microbiota vary across host individuals and environmental conditions, but the relative importance of their genetic background versus their environment is difficult to disentangle. We sought to experimentally determine the factors shaping the microbiota of the planktonic Crustacean, Daphnia magna. We used clonal lines from a wide geographic distribution, which had been kept under standardized conditions for over 75 generations. Replicate populations were kept for three generations at 20 and 28 °C. The interaction of the host clonal line and environment (i.e., temperature) influenced microbiota community characteristics, including structure, the relative abundance of common microbial species, and the microbial richness and phylogenetic diversity. We did not find any correlation between host-associated microbiota and the geographic origin of the clones or their temperature tolerance. Our results highlight the prominent effects that host clonal lineage and its interaction with the environment has on host-associated microbiota composition.

Microbial ecosystems are dominated by specialist taxa.

Abundance and specificity are two key characteristics of species distribution and biodiversity. Theories of species assembly aim to reproduce the empirical joint patterns of specificity and abundance, with the goal to explain patterns of biodiversity across habitats. The specialist-generalist paradigm predicts that specialists should have a local advantage over generalists and thus be more abundant. We developed a specificity index to analyse abundance-specificity relationships in microbial ecosystems. By analysing microbiota spanning 23 habitats from three very different data sets covering a wide range of sequencing depths and environmental conditions, we find that habitats are consistently dominated by specialist taxa, resulting in a strong, positive correlation between abundance and specificity. This finding is consistent over several levels of taxonomic aggregation and robust to errors in abundance measures. The relationship explains why shallow sequencing captures similar ß-diversity as deep sequencing, and can be sufficient to capture the habitat-specific functions of microbial communities.

Estimating bacterial diversity for ecological studies: methods, metrics, and assumptions.

Methods to estimate microbial diversity have developed rapidly in an effort to understand the distribution and diversity of microorganisms in natural environments. For bacterial communities, the 16S rRNA gene is the phylogenetic marker gene of choice, but most studies select only a specific region of the 16S rRNA to estimate bacterial diversity. Whereas biases derived from from DNA extraction, primer choice and PCR amplification are well documented, we here address how the choice of variable region can influence a wide range of standard ecological metrics, such as species richness, phylogenetic diversity, ß-diversity and rank-abundance distributions. We have used Illumina paired-end sequencing to estimate the bacterial diversity of 20 natural lakes across Switzerland derived from three trimmed variable 16S rRNA regions (V3, V4, V5). Species richness, phylogenetic diversity, community composition, ß-diversity, and rank-abundance distributions differed significantly between 16S rRNA regions. Overall, patterns of diversity quantified by the V3 and V5 regions were more similar to one another than those assessed by the V4 region. Similar results were obtained when analyzing the datasets with different sequence similarity thresholds used during sequences clustering and when the same analysis was used on a reference dataset of sequences from the Greengenes database. In addition we also measured species richness from the same lake samples using ARISA Fingerprinting, but did not find a strong relationship between species richness estimated by Illumina and ARISA. We conclude that the selection of 16S rRNA region significantly influences the estimation of bacterial diversity and species distributions and that caution is warranted when comparing data from different variable regions as well as when using different sequencing techniques.

De novo transcriptome hybrid assembly and validation in the European earwig (Dermaptera, Forficula auricularia).

BACKGROUND: The European earwig (Forficula auricularia) is an established system for studies of sexual selection, social interactions and the evolution of parental care. Despite its scientific interest, little knowledge exists about the species at the genomic level, limiting the scope of molecular studies and expression analyses of genes of interest. To overcome these limitations, we sequenced and validated the transcriptome of the European earwig. METHODOLOGY AND PRINCIPAL FINDINGS: To obtain a comprehensive transcriptome, we sequenced mRNA from various tissues and developmental stages of female and male earwigs using Roche 454 pyrosequencing and Illumina HiSeq. The reads were de novo assembled independently and screened for possible microbial contamination and repeated elements. The remaining contigs were combined into a hybrid assembly and clustered to reduce redundancy. A comparison with the eukaryotic core gene dataset indicates that we sequenced a substantial part of the earwig transcriptome with a low level of fragmentation. In addition, a comparative analysis revealed that more than 8,800 contigs of the hybrid assembly show significant similarity to insect-specific proteins and those were assigned for Gene Ontology terms. Finally, we established a quantitative PCR test for expression stability using commonly used housekeeping genes and applied the method to five homologs of known sex-biased genes of the honeybee. The qPCR pilot study confirmed sex specific expression and also revealed significant expression differences between the brain and antenna tissue samples. CONCLUSIONS: By employing two different sequencing approaches and including samples obtained from different tissues, developmental stages, and sexes, we were able to assemble a comprehensive transcriptome of F. auricularia. The transcriptome presented here offers new opportunities to study the molecular bases and evolution of parental care and sociality in arthropods.

The expression of one ankyrin pk2 allele of the WO prophage is correlated with the Wolbachia feminizing effect in isopods.

BACKGROUND: The maternally inherited α-Proteobacteria Wolbachia pipientis is an obligate endosymbiont of nematodes and arthropods, in which they induce a variety of reproductive alterations, including Cytoplasmic Incompatibility (CI) and feminization. The genome of the feminizing wVulC Wolbachia strain harboured by the isopod Armadillidium vulgare has been sequenced and is now at the final assembly step. It contains an unusually high number of ankyrin motif-containing genes, two of which are homologous to the phage-related pk1 and pk2 genes thought to contribute to the CI phenotype in Culex pipiens. These genes encode putative bacterial effectors mediating Wolbachia-host protein-protein interactions via their ankyrin motifs. RESULTS: To test whether these Wolbachia homologs are potentially involved in altering terrestrial isopod reproduction, we determined the distribution and expression of both pk1 and pk2 genes in the 3 Wolbachia strains that induce CI and in 5 inducing feminization of their isopod hosts. Aside from the genes being highly conserved, we found a substantial copy number variation among strains, and that is linked to prophage diversity. Transcriptional analyses revealed expression of one pk2 allele (pk2b2) only in the feminizing Wolbachia strains of isopods. CONCLUSIONS: These results reveal the need to investigate the functions of Wolbachia ankyrin gene products, in particular those of Pk2, and their host targets with respect to host sex manipulation.

Widespread Wolbachia infection in terrestrial isopods and other crustaceans.

Wolbachia bacteria are obligate intracellular alpha-Proteobacteria of arthropods and nematodes. Although widespread among isopod crustaceans, they have seldom been found in non-isopod crustacean species. Here, we report Wolbachia infection in fourteen new crustacean species. Our results extend the range of Wolbachia infections in terrestrial isopods and amphipods (class Malacostraca). We report the occurrence of two different Wolbachia strains in two host species (a terrestrial isopod and an amphipod). Moreover, the discovery of Wolbachia in the goose barnacle Lepas anatifera (subclass Thecostraca) establishes Wolbachia infection in class Maxillopoda. The new bacterial strains are closely related to B-supergroup Wolbachia strains previously reported from crustacean hosts. Our results suggest that Wolbachia infection may be much more widespread in crustaceans than previously thought. The presence of related Wolbachia strains in highly divergent crustacean hosts suggests that Wolbachia endosymbionts can naturally adapt to a wide range of crustacean hosts. Given the ability of isopod Wolbachia strains to induce feminization of genetic males or cytoplasmic incompatibility, we speculate that manipulation of crustacean-borne Wolbachia bacteria might represent potential tools for controlling crustacean species of commercial interest and crustacean or insect disease vectors.

Conservation of the Type IV secretion system throughout Wolbachia evolution.

The Type IV Secretion System (T4SS) is an efficient pathway with which bacteria can mediate the transfer of DNA and/or proteins to eukaryotic cells. In Wolbachia pipientis, a maternally inherited obligate endosymbiont of arthropods and nematodes, two operons of vir genes, virB3-B6 and virB8-D4, encoding a T4SS were previously identified and characterized at two separate genomic loci. Using the largest data set of Wolbachia strains studied so far, we show that vir gene sequence and organization are strictly conserved among 37 Wolbachia strains inducing various phenotypes such as cytoplasmic incompatibility, feminization, or oogenesis in their arthropod hosts. In sharp contrast, extensive variation of genomic sequences flanking the virB8-D4 operon suggested its distinct location among Wolbachia genomes. Long term conservation of the T4SS may imply maintenance of a functional effector translocation system in Wolbachia, thereby suggesting the importance for the T4SS in Wolbachia biology and survival inside host cells.

Viral cystatin evolution and three-dimensional structure modelling: a case of directional selection acting on a viral protein involved in a host-parasitoid interaction.

BACKGROUND: In pathogens, certain genes encoding proteins that directly interact with host defences coevolve with their host and are subject to positive selection. In the lepidopteran host-wasp parasitoid system, one of the most original strategies developed by the wasps to defeat host defences is the injection of a symbiotic polydnavirus at the same time as the wasp eggs. The virus is essential for wasp parasitism success since viral gene expression alters the immune system and development of the host. As a wasp mutualist symbiont, the virus is expected to exhibit a reduction in genome complexity and evolve under wasp phyletic constraints. However, as a lepidopteran host pathogenic symbiont, the virus is likely undergoing strong selective pressures for the acquisition of new functions by gene acquisition or duplication. To understand the constraints imposed by this particular system on virus evolution, we studied a polydnavirus gene family encoding cyteine protease inhibitors of the cystatin superfamily. RESULTS: We show that cystatins are the first bracovirus genes proven to be subject to strong positive selection within a host-parasitoid system. A generated three-dimensional model of Cotesia congregata bracovirus cystatin 1 provides a powerful framework to position positively selected residues and reveal that they are concentrated in the vicinity of actives sites which interact with cysteine proteases directly. In addition, phylogenetic analyses reveal two different cystatin forms which evolved under different selective constraints and are characterized by independent adaptive duplication events. CONCLUSION: Positive selection acts to maintain cystatin gene duplications and induces directional divergence presumably to ensure the presence of efficient and adapted cystatin forms. Directional selection has acted on key cystatin active sites, suggesting that cystatins coevolve with their host target. We can strongly suggest that cystatins constitute major virulence factors, as was already proposed in previous functional studies.

Characterization and transcriptional analysis of two gene clusters for type IV secretion machinery in Wolbachia of Armadillidium vulgare.

Wolbachia are maternally inherited alpha-proteobacteria that induce feminization of genetic males in most terrestrial crustacean isopods. Two clusters of vir genes for a type IV secretion machinery have been identified at two separate loci and characterized for the first time in a feminizing Wolbachia. Furthermore, we demonstrated that these operons are transcriptionally active in ovaries and in all other tissues tested, suggesting that T4SS has a significant role in Wolbachia biology. These observations and the identification of homologous vir genes in Wolbachia strains infecting insects or nematodes show that vir genes are conserved among Wolbachia strains whatever the phenotype induced by the bacteria.

Intense transpositional activity of insertion sequences in an ancient obligate endosymbiont.

The streamlined genomes of ancient obligate endosymbionts generally lack transposable elements, such as insertion sequences (IS). Yet, the genome of Wolbachia, one of the most abundant bacterial endosymbionts on Earth, is littered with IS. Such a paradox raises the question as to why there are so many ISs in the genome of this ancient endosymbiont. To address this question, we investigated IS transpositional activity in the unculturable Wolbachia by tracking the evolutionary dynamics and history of ISWpi1 elements. We show that 1) ISWpi1 is widespread in Wolbachia, being present in at least 55% of the 40 sampled strains, 2) ISWpi1 copies exhibit virtually identical nucleotide sequences both within and among Wolbachia genomes and possess an intact transposase gene, 3) individual ISWpi1 copies are differentially inserted among Wolbachia genomes, and 4) ISWpi1 occurs at variable copy numbers among Wolbachia genomes. Collectively, our results provide compelling evidence for intense ISWpi1 transpositional activity and frequent ISWpi1 horizontal transmission among strains during recent Wolbachia evolution. Thus, the genomes of ancient obligate endosymbionts can carry high loads of functional and transpositionally active transposable elements. Our results also indicate that Wolbachia genomes have experienced multiple and temporally distinct ISWpi1 invasions during their evolutionary history. Such recurrent exposition to new IS invasions may explain, at least partly, the unusually high density of transposable elements found in the genomes of Wolbachia endosymbionts.