Intra-lineage microevolution of Wolbachia leads to the emergence of new cytoplasmic incompatibility patterns.

Mosquitoes of the Culex pipiens complex are worldwide vectors of arbovirus, filarial nematodes, and avian malaria agents. In these hosts, the endosymbiotic bacteria Wolbachia induce cytoplasmic incompatibility (CI), i.e., reduced embryo viability in so-called incompatible crosses. Wolbachia infecting Culex pipiens (wPip) cause CI patterns of unparalleled complexity, associated with the amplification and diversification of cidA and cidB genes, with up to 6 different gene copies described in a single wPip genome. In wPip, CI is thought to function as a toxin-antidote (TA) system where compatibility relies on having the right antidotes (CidA) in the female to bind and neutralize the male's toxins (CidB). By repeating crosses between Culex isofemale lines over a 17 years period, we documented the emergence of a new compatibility type in real time and linked it to a change in cid genes genotype. We showed that loss of specific cidA gene copies in some wPip genomes results in a loss of compatibility. More precisely, we found that this lost antidote had an original sequence at its binding interface, corresponding to the original sequence at the toxin's binding interface. We showed that these original cid variants are recombinant, supporting a role for recombination rather than point mutations in rapid CI evolution. These results strongly support the TA model in natura, adding to all previous data acquired with transgenes expression.

Emergence and clonal expansion of Vibrio aestuarianus lineages pathogenic for oysters in Europe.

Crassostrea gigas oysters represent a significant global food source with 4.7 million tons harvested per year. In 2001, the bacterium Vibrio aestuarianus subsp. francensis emerged as a pathogen that causes adult oyster mortality in France and Ireland. Its impact on oyster aquaculture has increased in Europe since its re-emergence in 2012. To better understand the evolutionary mechanisms leading to the emergence and persistence over time of this pathogen, we conducted a survey of mollusc diseases through national reference laboratories across Europe. We analysed 54 new genomes of Vibrio aestuarianus (Va) isolated from multiple environmental compartments since 2001, in areas with and without bivalve mortalities. We used a combination of comparative genomics and population genetics approaches and show that Va has a classical epidemic population structure from which the pathogenic Va francensis subspecies emerged and clonally expanded. Furthermore, we identified a specific cus-cop-containing island conferring copper resistance to Va francensis whose acquisition may have favoured the emergence of pathogenic lineages adapted and specialized to oysters.

The Toxin-Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications.

Wolbachia bacteria inhabit the cells of about half of all arthropod species, an unparalleled success stemming in large part from selfish invasive strategies. Cytoplasmic incompatibility (CI), whereby the symbiont makes itself essential to embryo viability, is the most common of these and constitutes a promising weapon against vector-borne diseases. After decades of theoretical and experimental struggle, major recent advances have been made toward a molecular understanding of this phenomenon. As pieces of the puzzle come together, from yeast and Drosophila fly transgenesis to CI diversity patterns in natural mosquito populations, it becomes clearer than ever that the CI induction and rescue stem from a toxin-antidote (TA) system. Further, the tight association of the CI genes with prophages provides clues to the possible evolutionary origin of this phenomenon and the levels of selection at play.

RNA interference identifies domesticated viral genes involved in assembly and trafficking of virus-derived particles in ichneumonid wasps.

There are many documented examples of viral genes retained in the genomes of multicellular organisms that may in some cases bring new beneficial functions to the receivers. The ability of certain ichneumonid parasitic wasps to produce virus-derived particles, the so-called ichnoviruses (IVs), not only results from the capture and domestication of single viral genes but of almost entire ancestral virus genome(s). Indeed, following integration into wasp chromosomal DNA, the putative and still undetermined IV ancestor(s) evolved into encoding a 'virulence gene delivery vehicle' that is now required for successful infestation of wasp hosts. Several putative viral genes, which are clustered in distinct regions of wasp genomes referred to as IVSPERs (Ichnovirus Structural Protein Encoding Regions), have been assumed to be involved in virus-derived particles morphogenesis, but this question has not been previously functionally addressed. In the present study, we have successfully combined RNA interference and transmission electron microscopy to specifically identify IVSPER genes that are responsible for the morphogenesis and trafficking of the virus-derived particles in ovarian cells of the ichneumonid wasp Hyposoter didymator. We suggest that ancestral viral genes retained within the genomes of certain ichneumonid parasitoids possess conserved functions which were domesticated for the purpose of assembling viral vectors for the delivery of virulence genes to parasitized host animals.

The cellular phenotype of cytoplasmic incompatibility in Culex pipiens in the light of cidB diversity.

Wolbachia are maternally inherited endosymbiotic bacteria, widespread among arthropods thanks to host reproductive manipulations that increase their prevalence into host populations. The most commonly observed manipulation is cytoplasmic incompatibility (CI). CI leads to embryonic death in crosses between i) infected males and uninfected females and ii) individuals infected with incompatible Wolbachia strains. CI can be conceptualized as a toxin-antidote system where a toxin deposited by Wolbachia in the sperm would induce embryonic death unless countered by an antidote produced by Wolbachia present in the eggs. In Drosophila melanogaster, transgenic expression of Wolbachia effector cidB revealed its function of CI-inducing toxin. Moreover in Culex pipiens, the diversity of cidB variants present in wPip strains accounts for the diversity in crossing-types. We conducted cytological analyses to determine the CI mechanisms that lead to embryonic death in C. pipiens, and assess whether diversity in crossing-types could be based on variations in these mechanisms. We revealed that paternal chromatin condensation and segregation defects during the first embryonic division are always responsible for embryonic death. The strongest observed defects lead to an exclusion of the paternal chromatin from the first zygotic division, resulting in haploid embryos unable to hatch. The proportion of unhatched haploid embryos, developing with only maternal chromatin, which reflects the frequency of strong defects can be considered as a proxy of CI intensity at the cellular level. We thus studied the putative effect of variations in crossing types and cidB diversification on CI defects intensity. Incompatible crosses involving distinct wPip strains revealed that CI defects intensity depends on the Wolbachia strains hosted by the males and is linked to the diversity of cidB genes harbored in their genomes. These results support that, additionally to its implication in C. pipiens crossing type variability, cidB diversification also influences the strength of CI embryonic defects.

Wolbachia modulates prevalence and viral load of Culex pipiens densoviruses in natural populations.

The inadequacy of standard mosquito control strategies calls for ecologically safe novel approaches, for example the use of biological agents such as the endosymbiotic α-proteobacteria Wolbachia or insect-specific viruses (ISVs). Understanding the ecological interactions between these "biocontrol endosymbionts" is thus a fundamental step. Wolbachia are transmitted vertically from mother to offspring and modify their hosts' phenotypes, including reproduction (e.g., cytoplasmic incompatibility) and survival (e.g., viral interference). In nature, Culex pipiens (sensu lato) mosquitoes are always found infected with genetically diverse Wolbachia called wPip that belong to five phylogenetic groups. In recent years, ISVs have also been discovered in these mosquito species, although their interactions with Wolbachia in nature are unknown. Here, we studied the interactions between a widely prevalent ISV, the Culex pipiens densovirus (CpDV, Densovirinae), and Wolbachia in northern Tunisian C. pipiens populations. We showed an influence of different Wolbachia groups on CpDV prevalence and a general positive correlation between Wolbachia and CpDV loads. By investigating the putative relationship between CpDV diversification and wPip groups in the different sites, we detected a signal linked to wPip groups in CpDV phylogeny in sites where all larvae were infected by the same wPip group. However, no such signal was detected where the wPip groups coexisted, suggesting CpDV horizontal transfer between hosts. Overall, our results provide good evidence for an ecological influence of Wolbachia on an ISV, CpDV, in natural populations and highlight the importance of integrating Wolbachia in our understanding of ISV ecology in nature.

Variation in Wolbachia cidB gene, but not cidA, is associated with cytoplasmic incompatibility mod phenotype diversity in Culex pipiens.

Endosymbiotic Wolbachia bacteria are, to date, considered the most widespread symbionts in arthropods and are the cornerstone of major biological control strategies. Such a high prevalence is based on the ability of Wolbachia to manipulate their hosts' reproduction. One manipulation called cytoplasmic incompatibility (CI) is based on the death of the embryos generated by crosses between infected males and uninfected females or between individuals infected with incompatible Wolbachia strains. CI can be seen as a modification-rescue system (or mod-resc) in which paternal Wolbachia produce mod factors, inducing embryonic defects, unless the maternal Wolbachia produce compatible resc factors. Transgenic experiments in Drosophila melanogaster and Saccharomyces cerevisiae converged towards a model where the cidB Wolbachia gene is involved in the mod function while cidA is involved in the resc function. However, as cidA expression in Drosophila males was required to observe CI, it has been proposed that cidA could be involved in both resc and mod functions. A recent correlative study in natural Culex pipiens mosquito populations has revealed an association between specific cidA and cidB variations and changes in mod phenotype, also suggesting a role for both these genes in mod diversity. Here, by studying cidA and cidB genomic repertoires of individuals from newly sampled natural C. pipiens populations harbouring wPipIV strains from North Italy, we reinforce the link between cidB variation and mod phenotype variation fostering the involvement of cidB in the mod phenotype diversity. However, no association between any cidA variants or combination of cidA variants and mod phenotype variation was observed. Taken together our results in natural C. pipiens populations do not support the involvement of cidA in mod phenotype variation.

Sharing cells with Wolbachia: the transovarian vertical transmission of Culex pipiens densovirus.

Culex pipiens densovirus (CpDV), a single stranded DNA virus, has been isolated from Culex pipiens mosquitoes but differs from other mosquito densoviruses in terms of genome structure and sequence identity. Its transmission from host to host, the nature of its interactions with both its host and host's endosymbiotic bacteria Wolbachia are not known. Here, we report the presence of CpDV in the ovaries and eggs of Cx. pipiens mosquitoes in close encounters with Wolbachia. In the ovaries, CpDV amount significantly differed between mosquito lines harbouring different strains of Wolbachia and these differences were not linked to variations in Wolbachia densities. CpDV was vertically transmitted in all laboratory lines to 17%-20% of the offspring. For some females, however, the vertical transmission reached 90%. Antibiotic treatment that cured the host from Wolbachia significantly decreased both CpDV quantity and vertical transmission suggesting an impact of host microbiota, including Wolbachia, on CpDV transmission. Overall our results show that CpDV is transmitted vertically via transovarian path along with Wolbachia with which it shares the same cells. Our results are primordial to understand the dynamics of densovirus infection, their persistence and spread in populations considering their potential use in the regulation of mosquito vector populations.

Phenotypic shift in Wolbachia virulence towards its native host across serial horizontal passages.

Vertical transmission mode is predicted to decrease the virulence of symbionts. However, Wolbachia, a widespread vertically transmitted endosymbiont, exhibits both negative and beneficial effects on arthropod fitness. This 'Jekyll and Hyde' behaviour, as well as its ability to live transiently outside host cells and to establish new infections via horizontal transmission, may reflect the capacity of Wolbachia to exhibit various phenotypes depending on the prevailing environmental constraints. To study the ability of Wolbachia to readily cope with new constraints, we forced this endosymbiont to spread only via horizontal transmission. To achieve this, we performed serial horizontal transfers of haemolymph from Wolbachia-infected to naive individuals of the isopod Armadillidium vulgare. Across passages, we observed phenotypic changes in the symbiotic relationship: (i) The Wolbachia titre increased in both haemolymph and nerve cord but remained stable in ovaries; (ii) Wolbachia infection was benign at the beginning of the experiment, but highly virulent, killing most hosts after only a few passages. Such a phenotypic shift after recurrent horizontal passages demonstrates that Wolbachia can rapidly change its virulence when facing new environmental constraints. We thoroughly discuss the potential mechanism(s) underlying this phenotypic change, which are likely to be crucial for the ongoing radiation of Wolbachia in arthropods.

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.

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.

High virulence of Wolbachia after host switching: when autophagy hurts.

Wolbachia are widespread endosymbionts found in a large variety of arthropods. While these bacteria are generally transmitted vertically and exhibit weak virulence in their native hosts, a growing number of studies suggests that horizontal transfers of Wolbachia to new host species also occur frequently in nature. In transfer situations, virulence variations can be predicted since hosts and symbionts are not adapted to each other. Here, we describe a situation where a Wolbachia strain (wVulC) becomes a pathogen when transfected from its native terrestrial isopod host species (Armadillidium vulgare) to another species (Porcellio d. dilatatus). Such transfer of wVulC kills all recipient animals within 75 days. Before death, animals suffer symptoms such as growth slowdown and nervous system disorders. Neither those symptoms nor mortalities were observed after injection of wVulC into its native host A. vulgare. Analyses of wVulC's densities in main organs including Central Nervous System (CNS) of both naturally infected A. vulgare and transfected P. d. dilatatus and A. vulgare individuals revealed a similar pattern of host colonization suggesting an overall similar resistance of both host species towards this bacterium. However, for only P. d. dilatatus, we observed drastic accumulations of autophagic vesicles and vacuoles in the nerve cells and adipocytes of the CNS from individuals infected by wVulC. The symptoms and mortalities could therefore be explained by this huge autophagic response against wVulC in P. d. dilatatus cells that is not triggered in A. vulgare. Our results show that Wolbachia (wVulC) can lead to a pathogenic interaction when transferred horizontally into species that are phylogenetically close to their native hosts. This change in virulence likely results from the autophagic response of the host, strongly altering its tolerance to the symbiont and turning it into a deadly pathogen.

Strength of the pathogenicity caused by feminizing Wolbachia after transfer in a new host: strain or dose effect?

The alphaproteobacteria Wolbachia pipientis are among the most common and widespread symbionts in the animal world. Their vertical transmission mode is predicted to favour genotypes with low virulence. On the contrary, horizontal transfers of Wolbachia from one host to another have been shown to possibly increase the symbiont virulence. This situation has been previously described when two feminizing Wolbachia strains, wVulC and wVulM, from the ovaries of the woodlouse Armadillidium vulgare were introduced into another woodlouse named Porcellio dilatatus. These two Wolbachia strains induced severe symptoms and eventually caused the death of the recipient host. However, symptoms and death appeared sooner with wVulC than with wVulM. To know whether this difference was due to variation in the dose of infection or a difference in virulence between the two Wolbachia strains, we performed controlled and gradual doses of injection with wVulC and wVulM in P. dilatatus. We showed that the two strains differed intrinsically in their virulence against P. dilatatus and that their virulence is related to the injection dose. Moreover, we showed that wVulC reached higher concentrations in the recipient host than wVulM suggesting a potential link between the bacterial titers and the levels of virulence. We also addressed the impact of the tissue source of the Wolbachia used for the transinfection and demonstrated that Wolbachia transinfected via hemolymph colonized the body of the recipient more quickly and caused accelerated symptoms compared to Wolbachia introduced via a crushed ovaries suspension.

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.

Horizontal transfers of feminizing versus non-feminizing Wolbachia strains: from harmless passengers to pathogens.

The endosymbiont Wolbachia pipientis infects various hosts in which it navigates vertically from mothers to offspring. However, horizontal transfers of Wolbachia can occur between hosts. The virulence of the horizontally acquired Wolbachia can change in the new host as it has been illustrated by the case of the feminizing strain wVulC from the woodlouse Armadillidium vulgare that turns to a pathogen when introduced into Porcellio dilatatus dilatatus. In the present study, we aim to show whether symbiotic traits, such as (i) host sex manipulation and (ii) colonization patterns, which differ between eight isopod Wolbachia strains, are connected to their virulence towards the recipient host P. d. dilatatus. Among the transferred Wolbachia, some feminizing strains gradually differing in feminizing intensity in their native hosts induced different levels of pathogenicity to P. d. dilatatus. Not a single feminizing strain passed vertically with high titres to the next generation. The non-feminizing Wolbachia strains, even if they reached high densities in the host, did not impact host life-history traits and some vertically passed with high titres to the offspring. These results suggest that a potential link between the manners Wolbachia manipulates its native host reproduction, its virulence and its ability to vertically infect the offspring.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction.

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

From Wolbachia genomics to phenotype: molecular models of cytoplasmic incompatibility must account for the multiplicity of compatibility types.

Wolbachia endosymbionts commonly induce cytoplasmic incompatibility, making infected males' sperm lethal to the embryos unless these are rescued by the same bacterium, inherited from their mother. Causal genes were recently identified but two families of mechanistic models are still opposed. In the toxin-antidote model, interaction between the toxin and the antidote is required for rescuing the embryos. In host modification models, a host factor is misregulated in sperm and rescue occurs through compensation or withdrawal of this modification. While these models have been thoroughly discussed, the multiplicity of compatibility types, that is, the existence of many mutually incompatible strains, as seen in Culex mosquitoes, has not received sufficient attention. To explain such a fact, host modification models must posit that the same embryonic defects can be induced and rescued through a large variety of host targets. Conversely, the toxin-antidote model simply accommodates this pattern in a lock-key fashion, through variations in the toxin-antidote interaction sites.

Paternal transmission of the Wolbachia CidB toxin underlies cytoplasmic incompatibility.

Wolbachia are widespread endosymbiotic bacteria that manipulate the reproduction of arthropods through a diversity of cellular mechanisms. In cytoplasmic incompatibility (CI), a sterility syndrome originally discovered in the mosquito Culex pipiens, uninfected eggs fertilized by sperm from infected males are selectively killed during embryo development following the abortive segregation of paternal chromosomes in the zygote. Despite the recent discovery of Wolbachia CI factor (cif) genes, the mechanism by which they control the fate of paternal chromosomes at fertilization remains unknown. Here, we have analyzed the cytological distribution and cellular impact of CidA and CidB, a pair of Cif proteins from the Culex-infecting Wolbachia strain wPip. We show that expression of CidB in Drosophila S2R+ cells induces apoptosis unless CidA is co-expressed and associated with its partner. In transgenic Drosophila testes, both effectors colocalize in germ cells until the histone-to-protamine transition in which only CidB is retained in maturing spermatid nuclei. We further show that CidB is similarly targeted to maturing sperm of naturally infected Culex mosquitoes. At fertilization, CidB associates with paternal DNA regions exhibiting DNA replication stress, as a likely cause of incomplete replication of paternal chromosomes at the onset of the first mitosis. Importantly, we demonstrate that inactivation of the deubiquitylase activity of CidB does not abolish its cell toxicity or its ability to induce CI in Drosophila. Our study thus demonstrates that CI functions as a transgenerational toxin-antidote system and suggests that CidB acts by poisoning paternal DNA replication in incompatible crosses.

The mosquito microbiome includes habitat-specific but rare symbionts.

Microbial communities are known to influence mosquito lifestyles by modifying essential metabolic and behavioral processes that affect reproduction, development, immunity, digestion, egg survival, and the ability to transmit pathogens. Many studies have used 16S rRNA gene amplicons to characterize mosquito microbiota and investigate factors that influence host-microbiota dynamics. However, a relatively low taxonomic resolution due to clustering methods based on arbitrary threshold and the overall dominance of Wolbachia or Asaia symbionts obscured the investigation of rare members of mosquito microbiota in previous studies. Here, we used high resolution Shannon entropy-based oligotyping approach to analyze the microbiota of Culex pipiens, Culex quinquefasciatus and Aedes individuals from continental Southern France and overseas Guadeloupe as well as from laboratories with or without antibiotics treatment. Our experimental design that resulted in a series of mosquito samples with a gradient of Wolbachia density and relative abundance along with high-resolution analyses of amplicon sequences enabled the recovery of a robust signal from typically less accessible bacterial taxa. Our data confirm species-specific mosquito-bacteria associations with geography as a primary factor that influences bacterial community structure. But interestingly, they also reveal co-occurring symbiotic bacterial variants within single individuals for both Elizabethkingia and Erwinia genera, distinct and specific Asaia and Chryseobacterium in continental and overseas territories, and a putative rare Wolbachia variant. Overall, our study reveals the presence of previously overlooked microdiversity and multiple closely related symbiotic strains within mosquito individuals with a remarkable habitat-specificity.