Feminising Wolbachia disrupt Armadillidium vulgare insulin-like signalling pathway.

The endosymbiont Wolbachia feminises male isopods by making them refractory to the insulin-like masculinising hormone, which shunts the autocrine development of the androgenic glands. It was, therefore, proposed that Wolbachia silences the IR receptors, either by preventing their expression or by inactivating them. We describe here the two IR paralogs of Armadillidium vulgare. They displayed a conventional structure and belonged to a family widespread among isopods. Av-IR1 displayed an ubiquist expression, whereas the expression of Av-IR2 was restricted to the gonads. Both were constitutively expressed in males and females and throughout development. However, upon silencing, altered gland physiology and gene expression therein suggested antagonistic roles for Av-IR1 (androinhibiting) and Av-IR2 (androstimulating). They may function in tandem with regulating neurohormones, as a conditional platform that conveys insulin signalling. Wolbachia infection did not alter their expression patterns: leaving the IRs unscathed, the bacteria would suppress the secretion of the neurohormones, thus inducing body-wide IR deactivation and feminisation. Adult males injected with Wolbachia acquired an intersexed physiology. Their phenotypes and gene expressions mirrored the silencing of Av-IR1 only, suggesting that imperfect feminisation stems from a flawed invasion of the androstimulating centre, whereas in fully feminised males invasion would be complete in early juveniles. TAKE AWAY: Two antagonistic Insulin Receptors were characterised in Armadillidium vulgare. The IRs were involved in androstimulating and androinhibiting functions. Wolbachia-induced feminisation did not prevent the expression of the IRs. Imperfectly feminised intersexes phenocopied the silencing of Av-IR1 only. Wolbachia would deactivate the IRs by suppressing neurosecretory co-factors.

Lignocellulose degradation in isopods: new insights into the adaptation to terrestrial life.

BACKGROUND: Isopods constitute a particular group of crustaceans that has successfully colonized all environments including marine, freshwater and terrestrial habitats. Their ability to use various food sources, especially plant biomass, might be one of the reasons of their successful spread. All isopods, which feed on plants and their by-products, must be capable of lignocellulose degradation. This complex composite is the main component of plants and is therefore an important nutrient source for many living organisms. Its degradation requires a large repertoire of highly specialized Carbohydrate-Active enZymes (called CAZymes) which are produced by the organism itself and in some cases, by its associated microbiota. The acquisition of highly diversified CAZymes could have helped isopods to adapt to their diet and to their environment, especially during land colonization. RESULTS: To test this hypothesis, isopod host CAZomes (i.e. the entire CAZyme repertoire) were characterized in marine, freshwater and terrestrial species through a transcriptomic approach. Many CAZymes were identified in 64 isopod transcriptomes, comprising 27 de novo datasets. Our results show that marine, freshwater and terrestrial isopods exhibit different CAZomes, illustrating different strategies for lignocellulose degradation. The analysis of variations of the size of CAZy families shows these are expanded in terrestrial isopods while they are contracted in aquatic isopods; this pattern is probably resulting from the evolution of the host CAZomes during the terrestrial adaptation of isopods. We show that CAZyme gene duplications and horizontal transfers can be involved in adaptive divergence between isopod CAZomes. CONCLUSIONS: Our characterization of the CAZomes in 64 isopods species provides new insights into the evolutionary processes that enabled isopods to conquer various environments, especially terrestrial ones.

IGFBP-rP1, a strongly conserved member of the androgenic hormone signalling pathway in Isopoda.

The first protein which has been described to interact with the malacostracan Androgenic Gland Hormone (AGH) is a binding protein called IGFBP-rP1. It has been identified and studied in several species of decapods, in which its interaction with the masculinizing hormone and its expression patterns have been established in several ways. However, this protein remains uncharacterised to date in the other malacostracan orders, like Amphipoda and Isopoda, although they were historically the first ones in which the androgenic gland and the corresponding hormone were respectively described. In this article, we identified the IGFBP-rP1 of isopods and established its implication in the pathway of the AGH with a silencing approach in the model species Armadillidium vulgare. We also showed that this gene is expressed in all the tissues of males and females, with a similar pattern in animals infected with Wolbachia, a feminizing endosymbiont of several isopod species. The expression pattern did not differ during the development of uninfected and infected animals either. We finally studied the evolution of the IGFBP-rP1 in 68 isopod species, looking for conserved motifs and evidence of natural selection. Altogether, our results showed that this gene is constitutively expressed and strongly conserved in isopods, in which it likely constitutes a key element of the insulin/IGF signalling pathway. However, we also illustrated that IGFBP-rP1 is not sufficient on its own to explain the different developmental paths taken by the males and the females or feminized genetic males.

Birth of a W sex chromosome by horizontal transfer of Wolbachia bacterial symbiont genome.

Sex determination is a fundamental developmental pathway governing male and female differentiation, with profound implications for morphology, reproductive strategies, and behavior. In animals, sex differences between males and females are generally determined by genetic factors carried by sex chromosomes. Sex chromosomes are remarkably variable in origin and can differ even between closely related species, indicating that transitions occur frequently and independently in different groups of organisms. The evolutionary causes underlying sex chromosome turnover are poorly understood, however. Here we provide evidence indicating that Wolbachia bacterial endosymbionts triggered the evolution of new sex chromosomes in the common pillbug Armadillidium vulgare We identified a 3-Mb insert of a feminizing Wolbachia genome that was recently transferred into the pillbug nuclear genome. The Wolbachia insert shows perfect linkage to the female sex, occurs in a male genetic background (i.e., lacking the ancestral W female sex chromosome), and is hemizygous. Our results support the conclusion that the Wolbachia insert is now acting as a female sex-determining region in pillbugs, and that the chromosome carrying the insert is a new W sex chromosome. Thus, bacteria-to-animal horizontal genome transfer represents a remarkable mechanism underpinning the birth of sex chromosomes. We conclude that sex ratio distorters, such as Wolbachia endosymbionts, can be powerful agents of evolutionary transitions in sex determination systems in animals.

Divergent evolution and clade-specific duplications of the Insulin-like Receptor in malacostracan crustaceans.

The Insulin-like Receptors (IRs) are an important protein family, represented by three members in vertebrates, two of which are well-known for their implication in metabolism (Insulin Receptor) and growth (IGF Receptor). In contrast, little is known about these receptors in invertebrates, in which a single gene generally exists except for a part of insects and other occasional species-specific duplications. In this study, we used publicly available sequences as well as de novo assembled transcriptomes to investigate the IR evolution in malacostracan crustaceans, animals in which the Insulin/IGF pathway is known to be implicated in sexual development through the androgenic gland hormone. We described the evolutionary divergence of malacostracan IRs compared to all the other metazoan sequences, including other pancrustaceans. We also demonstrated two well conserved duplications of IRs: one specific to the whole malacostracan class, another one specific to the decapod order. The potential implications for malacostracan biology are discussed.

Isolation, characterization and PCR multiplexing of microsatellite loci for two sub-species of terrestrial isopod Porcellio dilatatus (Crustacea, Oniscidea).

Several microsatellite markers have already been developed for different terrestrial isopod species such as Armadillidium vulgare, A. nasatum and Porcellionides pruinosus. In all these species, the endosymbiont Wolbachia has a feminizing effect that generates a female bias in sex ratio and reduces the number of reproductive males. Thus this can potentially decrease the genetic diversity of host populations. However, in some other isopod species, Wolbachia induces cytoplasmic incompatibility (CI); the most commonly described effect of Wolbachia in arthropods. The CI by rendering some crossings incompatible can reduce the gene flow and strengthen genetic differentiation between isopod populations. To date, the influence of Wolbachia inducing CI on population structure of terrestrial isopods has never been investigated. In this study, we developed 10 polymorphic microsatellite markers shared by two sub-species of Porcellio dilatatus. Crossings between the two sub-species are partially incompatible due to two CI-inducing Wolbachia strains. These new microsatellite markers will allow us to investigate the effect of CI on host genetic differentiation in this species complex.

The impact of endosymbionts on the evolution of host sex-determination mechanisms.

The past years have revealed that inherited bacterial endosymbionts are important sources of evolutionary novelty for their eukaryotic hosts. In this review we discuss a fundamental biological process of eukaryotes influenced by bacterial endosymbionts: the mechanisms of sex determination. Because they are maternally inherited, several endosymbionts of arthropods, known as reproductive parasites, have developed strategies to convert non-transmitting male hosts into transmitting females through feminization of genetic males and parthenogenesis induction. Recent investigations have also highlighted that endosymbionts can impact upon host sex determination more subtly through genetic conflicts, resulting in selection of host nuclear genes resisting endosymbiont effects. Paradoxically, it is because of their selfish nature that reproductive parasites are such powerful agents of evolutionary change in their host sex-determination mechanisms. They might therefore represent excellent models for studying transitions between sex-determining systems and, more generally, the evolution of sex-determination mechanisms in eukaryotes.

A host as an ecosystem: Wolbachia coping with environmental constraints.

The Wolbachia are intracellular endosymbionts widely distributed among invertebrates. These primarily vertically transmitted α-proteobacteria have been intensively studied during the last decades because of their intriguing interactions with hosts, ranging from reproductive manipulations to mutualism. To optimize their vertical transmission from mother to offspring, the Wolbachia have developed fine-tuned strategies. However, the Wolbachia are not restricted to the female gonads and frequently exhibit wide intra-host distributions. This extensive colonization of somatic organs might be necessary for Wolbachia to develop their diverse extended phenotypes. From an endosymbiont's perspective, the within-host environment potentially presents different environmental constraints. Hence, the Wolbachia have to face different intracellular habitats, their host's immune system as well as other microorganisms co-occurring in the same host individual and sometimes even in the same cell. A means for the Wolbachia to protect themselves from these environmental constraints may be to live 'hidden' in vacuoles within host cells. In this review, we summarize the current knowledge regarding the extent of the Wolbachia pandemic and discuss the various environmental constraints these bacteria may have to face within their 'host ecosystem'. Finally, we identify new avenues for future research to better understand the complexity of Wolbachia's interactions with their intracellular environment.

Bidirectional cytoplasmic incompatibility caused by Wolbachia in the terrestrial isopod Porcellio dilatatus.

In the terrestrial isopod species Porcellio dilatatus, unidirectional Cytoplasmic Incompatibility (CI) between two morphs (P. d. dilatatus and P. d. petiti) caused by a Wolbachia strain (wPet) infecting the morph P. d. petiti has been previously described by experiments initiated four decades ago. Here, we studied another Wolbachia that has been recently detected in a population of the morph P. d. dilatatus. The MLST markers reveal that this Wolbachia is a new strain called wDil distinct from wPet also belonging to the isopod clade of Wolbachia. Quantifications of both Wolbachia strains in the gonads of the two P. dilatatus morphs revealed that all males exhibit similar Wolbachia titers while the titers in females depend on the Wolbachia strain they host. Crossing experiments showed that both wDil and wPet induced partial unidirectional CI with different intensities. Moreover, these two strains induced bidirectional CI when individuals were both infected with one of the two different Wolbachia strains. This way, we demonstrated that P. dilatatus can be infected by two closely related Wolbachia strains (wDil and wPet), that seem to have different modification-rescue systems.

The terrestrial isopod symbiont 'Candidatus Hepatincola porcellionum' is a potential nutrient scavenger related to Holosporales symbionts of protists.

The order Holosporales (Alphaproteobacteria) encompasses obligate intracellular bacterial symbionts of diverse Eukaryotes. These bacteria have highly streamlined genomes and can have negative fitness effects on the host. Herein, we present a comparative analysis of the first genome sequences of 'Ca. Hepatincola porcellionum', a facultative symbiont occurring extracellularly in the midgut glands of terrestrial isopods. Using a combination of long-read and short-read sequencing, we obtained the complete circular genomes of two Hepatincola strains and an additional metagenome-assembled draft genome. Phylogenomic analysis validated its phylogenetic position as an early-branching family-level clade relative to all other established Holosporales families associated with protists. A 16S rRNA gene survey revealed that this new family encompasses diverse bacteria associated with both marine and terrestrial host species, which expands the host range of Holosporales bacteria from protists to several phyla of the Ecdysozoa (Arthropoda and Priapulida). Hepatincola has a highly streamlined genome with reduced metabolic and biosynthetic capacities as well as a large repertoire of transmembrane transporters. This suggests that this symbiont is rather a nutrient scavenger than a nutrient provider for the host, likely benefitting from a nutrient-rich environment to import all necessary metabolites and precursors. Hepatincola further possesses a different set of bacterial secretion systems compared to protist-associated Holosporales, suggesting different host-symbiont interactions depending on the host organism.

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.

Feminizing Wolbachia: a transcriptomics approach with insights on the immune response genes in Armadillidium vulgare.

BACKGROUND: Wolbachia are vertically transmitted bacteria known to be the most widespread endosymbiont in arthropods. They induce various alterations of the reproduction of their host, including feminization of genetic males in isopod crustaceans. In the pill bug Armadillidium vulgare, the presence of Wolbachia is also associated with detrimental effects on host fertility and lifespan. Deleterious effects have been demonstrated on hemocyte density, phenoloxidase activity, and natural hemolymph septicemia, suggesting that infected individuals could have defective immune capacities. Since nothing is known about the molecular mechanisms involved in Wolbachia-A. vulgare interactions and its secondary immunocompetence modulation, we developed a transcriptomics strategy and compared A. vulgare gene expression between Wolbachia-infected animals (i.e., "symbiotic" animals) and uninfected ones (i.e., "asymbiotic" animals) as well as between animals challenged or not challenged by a pathogenic bacteria. RESULTS: Since very little genetic data is available on A. vulgare, we produced several EST libraries and generated a total of 28 606 ESTs. Analyses of these ESTs revealed that immune processes were over-represented in most experimental conditions (responses to a symbiont and to a pathogen). Considering canonical crustacean immune pathways, these genes encode antimicrobial peptides or are involved in pathogen recognition, detoxification, and autophagy. By RT-qPCR, we demonstrated a general trend towards gene under-expression in symbiotic whole animals and ovaries whereas the same gene set tends to be over-expressed in symbiotic immune tissues. CONCLUSION: This study allowed us to generate the first reference transcriptome ever obtained in the Isopoda group and to identify genes involved in the major known crustacean immune pathways encompassing cellular and humoral responses. Expression of immune-related genes revealed a modulation of host immunity when females are infected by Wolbachia, including in ovaries, the crucial tissue for the Wolbachia route of transmission.

Targeted genome enrichment for efficient purification of endosymbiont DNA from host DNA.

Wolbachia endosymbionts are widespread in arthropods and are generally considered reproductive parasites, inducing various phenotypes including cytoplasmic incompatibility, parthenogenesis, feminization and male killing, which serve to promote their spread through populations. In contrast, Wolbachia infecting filarial nematodes that cause human diseases, including elephantiasis and river blindness, are obligate mutualists. DNA purification methods for efficient genomic sequencing of these unculturable bacteria have proven difficult using a variety of techniques. To efficiently capture endosymbiont DNA for studies that examine the biology of symbiosis, we devised a parallel strategy to an earlier array-based method by creating a set of SureSelect™ (Agilent) 120-mer target enrichment RNA oligonucleotides ("baits") for solution hybrid selection. These were designed from Wolbachia complete and partial genome sequences in GenBank and were tiled across each genomic sequence with 60 bp overlap. Baits were filtered for homology against host genomes containing Wolbachia using BLAT and sequences with significant host homology were removed from the bait pool. Filarial parasite Brugia malayi DNA was used as a test case, as the complete sequence of both Wolbachia and its host are known. DNA eluted from capture was size selected and sequencing samples were prepared using the NEBNext® Sample Preparation Kit. One-third of a 50 nt paired-end sequencing lane on the HiSeq™ 2000 (Illumina) yielded 53 million reads and the entirety of the Wolbachia genome was captured. We then used the baits to isolate more than 97.1 % of the genome of a distantly related Wolbachia strain from the crustacean Armadillidium vulgare, demonstrating that the method can be used to enrich target DNA from unculturable microbes over large evolutionary distances.

Legionella spp. All Ears? The Broad Occurrence of Quorum Sensing Elements outside Legionella pneumophila.

Legionella spp. are ubiquitous bacteria principally found in water networks and ∼20 species are implicated in Legionnaire's disease. Among them, Legionella pneumophila is an intracellular pathogen of environmental protozoa, responsible for ∼90% of cases in the world. Legionella pneumophila regulates in part its virulence by a quorum sensing system named "Legionella quorum sensing," composed of a signal synthase LqsA, two histidine kinase membrane receptors LqsS and LqsT and a cytoplasmic receptor LqsR. To date, this communication system was only found in L. pneumophila. Here, we investigated 58 Legionella genomes to determine the presence of a lqs cluster or homologous receptors using TBlastN. This analysis revealed three categories of species: 19 harbored a complete lqs cluster, 20 did not possess lqsA but maintained the receptor lqsR and/or lqsS, and 19 did not have any of the lqs genes. No correlation was observed between pathogenicity and the presence of a quorum sensing system. We determined by RT-qPCR that the lqsA gene was expressed at least in four strains among different species available in our laboratory. Furthermore, we showed that the lqs genomic region was conserved even in species possessing only the receptors of the quorum sensing system, indicating an ancestral acquisition and various loss dynamics during evolution. This system could therefore function in interspecific communication as well.

Experimental evidence of Wolbachia introgressive acquisition between terrestrial isopod subspecies.

Wolbachia are the most widespread endosymbiotic bacteria in animals. In many arthropod host species, they manipulate reproduction via several mechanisms that favor their maternal transmission to offspring. Among them, cytoplasmic incompatibility (CI) promotes the spread of the symbiont by specifically decreasing the fertility of crosses involving infected males and uninfected females, via embryo mortality. These differences in reproductive efficiency may select for the avoidance of incompatible mating, a process called reinforcement, and thus contribute to population divergence. In the terrestrial isopod Porcellio dilatatus, the Wolbachia wPet strain infecting the subspecies P. d. petiti induces unidirectional CI with uninfected individuals of the subspecies P. d. dilatatus. To study the consequences of CI on P. d. dilatatus and P. d. petiti hybridization, mitochondrial haplotypes and Wolbachia infection dynamics, we used population cages seeded with different proportions of the 2 subspecies in which we monitored these genetic parameters 5 and 7 years after the initial setup. Analysis of microsatellite markers allowed evaluating the degree of hybridization between individuals of the 2 subspecies. These markers revealed an increase in P. d. dilatatus nuclear genetic signature in all mixed cages, reflecting an asymmetry in hybridization. Hybridization led to the introgressive acquisition of Wolbachia and mitochondrial haplotype from P. d. petiti into nuclear genomes dominated by alleles of P. d. dilatatus. We discuss these results with regards to Wolbachia effects on their host (CI and putative fitness cost), and to a possible reinforcement that may have led to assortative mating, as possible factors contributing to the observed results.

Variations of immune parameters in terrestrial isopods: a matter of gender, aging and Wolbachia.

Ecological factors modulate animal immunocompetence and potentially shape the evolution of their immune systems. Not only environmental parameters impact on immunocompetence: Aging is one major cause of variability of immunocompetence between individuals, and sex-specific levels of immunocompetence have also been frequently described. Moreover, a growing core of data put in light that vertically transmitted symbionts can dramatically modulate the immunocompetence of their hosts. In this study, we addressed the influence of gender, age and the feminising endosymbiont Wolbachia (wVulC) on variations in haemocyte density, total PO activity and bacterial load in the haemolymph of the terrestrial isopod Armadillidium vulgare. This host-symbiont system is of particular interest to address this question since: (1) wVulC was previously shown as immunosuppressive in middle-aged females and (2) wVulC influences sex determination. We show that age, gender and Wolbachia modulate together immune parameters in A. vulgare. However, wVulC, which interacts with aging, appears to be the prominent factor interfering with both PO activity and haemocyte density. This interference with immune parameters is not the only aspect of wVulC virulence on its host, as reproduction and survival are also altered.

Isopod holobionts as promising models for lignocellulose degradation.

BACKGROUND: Isopods have colonized all environments, partly thanks to their ability to decompose the organic matter. Their enzymatic repertoire, as well as the one of their associated microbiota, has contributed to their colonization success. Together, these holobionts have evolved several interesting life history traits to degrade the plant cell walls, mainly composed of lignocellulose. It has been shown that terrestrial isopods achieve lignocellulose degradation thanks to numerous and diverse CAZymes provided by both the host and its microbiota. Nevertheless, the strategies for lignocellulose degradation seem more diversified in isopods, in particular in aquatic species which are the least studied. Isopods could be an interesting source of valuable enzymes for biotechnological industries of biomass conversion. RESULTS: To provide new features on the lignocellulose degradation in isopod holobionts, shotgun sequencing of 36 metagenomes of digestive and non-digestive tissues was performed from several populations of four aquatic and terrestrial isopod species. Combined to the 15 metagenomes of an additional species from our previous study, as well as the host transcriptomes, this large dataset allowed us to identify the CAZymes in both the host and the associated microbial communities. Analyses revealed the dominance of Bacteroidetes and Proteobacteria in the five species, covering 36% and 56% of the total bacterial community, respectively. The identification of CAZymes and new enzymatic systems for lignocellulose degradation, such as PULs, cellulosomes and LPMOs, highlights the richness of the strategies used by the isopods and their associated microbiota. CONCLUSIONS: Altogether, our results show that the isopod holobionts are promising models to study lignocellulose degradation. These models can provide new enzymes and relevant lignocellulose-degrading bacteria strains for the biotechnological industries of biomass conversion.

The shutting down of the insulin pathway: a developmental window for Wolbachia load and feminization.

Using the isopod Armadillidium vulgare as a case study, we review the significance of the "bacterial dosage model", which connects the expression of the extended phenotype to the rise of the Wolbachia load. In isopods, the Insulin-like Androgenic Gland hormone (IAG) induces male differentiation: Wolbachia feminizes males through insulin resistance, presumably through defunct insulin receptors. This should prevent an autocrine development of the androgenic glands so that females differentiate instead: feminization should translate as IAG silencing and increased Wolbachia load in the same developmental window. In line with the autocrine model, uninfected males expressed IAG from the first larval stage on, long before the androgenic gland primordia begin to differentiate, and exponentially throughout development. In contrast in infected males, expression fully stopped at stage 4 (juvenile), when male differentiation begins. This co-occurred with the only significant rise in the Wolbachia load throughout the life-stages. Concurrently, the raw expression of the bacterial Secretion Systems co-increased, but they were not over-expressed relative to the number of bacteria. The isopod model leads to formulate the "bacterial dosage model" throughout extended phenotypes as the conjunction between bacterial load as the mode of action, timing of multiplication (pre/post-zygotic), and site of action (soma vs. germen).

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.

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.