Signs of neutralization in a redundant gene involved in homologous recombination in Wolbachia endosymbionts

Temporal stability of sex ratio distorter prevalence in natural populations of the isopod Armadillidium vulgare

In the terrestrial isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios due to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and a nuclear non-mendelian locus called the f element. To investigate the potential impact of these SRD on the evolution of host sex determination, we analyzed their temporal distribution in six A. vulgare populations sampled between 2003 and 2017, for a total of 29 time points. SRD distribution was heterogeneous among populations despite their close geographic locations, so that when one SRD was frequent in a population, the other SRD was rare. In contrast with spatial heterogeneity, our results overall did not reveal substantial temporal variability in SRD prevalence within populations, suggesting equilibria in SRD evolutionary dynamics may have been reached or nearly so. Temporal stability was also generally reflected in mitochondrial and nuclear variation. Nevertheless, in a population, a Wolbachia strain replacement coincided with changes in mitochondrial composition but no change in nuclear composition, thus constituting a typical example of mitochondrial sweep caused by endosymbiont rise in frequency. Rare incongruence between Wolbachia strains and mitochondrial haplotypes suggested the occurrence of intraspecific horizontal transmission, making it a biologically relevant parameter for Wolbachia evolutionary dynamics in A. vulgare. Overall, our results provide an empirical basis for future studies on SRD evolutionary dynamics in the context of multiple sex determination factors co-existing within a single species, to ultimately evaluate the impact of SRD on the evolution of host sex determination mechanisms and sex chromosomes.

Three feminizing Wolbachia strains in a single host species: comparative genomics paves the way for identifying sex reversal factors

Introduction

Endosymbiotic bacteria in the genus Wolbachia have evolved numerous strategies for manipulating host reproduction in order to promote their own transmission. This includes the feminization of males into functional females, a well-studied phenotype in the isopod Armadillidium vulgare. Despite an early description of this phenotype in isopods and the development of an evolutionary model of host sex determination in the presence of Wolbachia, the underlying genetic mechanisms remain elusive.

Methods

Here we present the first complete genomes of the three feminizing Wolbachia (wVulC, wVulP, and wVulM) known to date in A. vulgare. These genomes, belonging to Wolbachia B supergroup, contain a large number of mobile elements such as WO prophages with eukaryotic association modules. Taking advantage of these data and those of another Wolbachia-derived feminizing factor integrated into the host genome (f element), we used a comparative genomics approach to identify putative feminizing factors.

Results

This strategy has enabled us to identify three prophage-associated genes secreted by the Type IV Secretion System: one ankyrin repeat domain-containing protein, one helix-turn-helix transcriptional regulator and one hypothetical protein. In addition, a latrotoxin-related protein, associated with phage relic genes, was shared by all three genomes and the f element.

Conclusion

These putative feminization-inducing proteins shared canonical interaction features with eukaryotic proteins. These results pave the way for further research into the underlying functional interactions.

Heterogeneous distribution of sex ratio distorters in natural populations of the isopod Armadillidium vulgare

In the isopod Armadillidium vulgare, many females produce progenies with female-biased sex ratios, owing to two feminizing sex ratio distorters (SRD): Wolbachia endosymbionts and the f element. We investigated the distribution and population dynamics of these SRD and mitochondrial DNA variation in 16 populations from Europe and Japan. Confirming and extending results from the 1990s, we found that the SRD are present at variable frequencies in populations and that the f element is overall more frequent than Wolbachia. The two SRD never co-occur at high frequency in any population, suggesting an apparent mutual exclusion. We also detected Wolbachia or the f element in some males, which probably reflects insufficient titer to induce feminization or presence of masculinizing alleles. Our results are consistent with a single integration event of a Wolbachia genome in the A. vulgare genome at the origin of the f element, which contradicts an earlier hypothesis of frequent losses and gains. We identified strong linkage between Wolbachia strains and mitochondrial haplotypes, but no association between the f element and mitochondrial background. Our results open new perspectives on SRD evolutionary dynamics in A. vulgare, the evolution of genetic conflicts and their impact on the variability of sex determination systems.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

Narrow Genetic Diversity of Wolbachia Symbionts in Acrididae Grasshopper Hosts (Insecta, Orthoptera)

Bacteria of the Wolbachia genus are maternally inherited symbionts of Nematoda and numerous Arthropoda hosts. There are approximately 20 lineages of Wolbachia, which are called supergroups, and they are designated alphabetically. Wolbachia strains of the supergroups A and B are predominant in arthropods, especially in insects, and supergroup F seems to rank third. Host taxa have been studied very unevenly for Wolbachia symbionts, and here, we turn to one of largely unexplored insect families: Acrididae. On the basis of five genes subject to multilocus sequence typing, we investigated the incidence and genetic diversity of Wolbachia in 41 species belonging three subfamilies (Gomphocerinae, Oedipodinae, and Podisminae) collected in Turkey, Kazakhstan, Tajikistan, Russia, and Japan, making 501 specimens in total. Our results revealed a high incidence and very narrow genetic diversity of Wolbachia. Although only the strains belonging to supergroups A and B are commonly present in present, the Acrididae hosts here proved to be infected with supergroups B and F without A-supergroup variants. The only trace of an A-supergroup lineage was noted in one case of an inter-supergroup recombinant haplotype, where the ftsZ gene came from supergroup A, and the others from supergroup B. Variation in the Wolbachia haplotypes in Acrididae hosts within supergroups B and F was extremely low. A comprehensive genetic analysis of Wolbachia diversity confirmed specific features of the Wolbachia allelic set in Acrididae hosts. This result can help to elucidate the crucial issue of Wolbachia biology: the route(s) and mechanism(s) of Wolbachia horizontal transmission.

DNA recombination and repair in Wolbachia: RecA and related proteins

Wolbachia is an obligate intracellular bacterium that has undergone extensive genomic streamlining in its arthropod and nematode hosts. Because the gene encoding the bacterial DNA recombination/repair protein RecA is not essential in Escherichia coli, abundant expression of this protein in a mosquito cell line persistently infected with Wolbachia strain wStri was unexpected. However, RecA's role in the lytic cycle of bacteriophage lambda provides an explanation for retention of recA in strains known to encode lambda-like WO prophages. To examine DNA recombination/repair capacities in Wolbachia, a systematic examination of RecA and related proteins in complete or nearly complete Wolbachia genomes from supergroups A, B, C, D, E, F, J and S was undertaken. Genes encoding proteins including RecA, RecF, RecO, RecR, RecG and Holliday junction resolvases RuvA, RuvB and RuvC are uniformly absent from Wolbachia in supergroup C and have reduced representation in supergroups D and J, suggesting that recombination and repair activities are compromised in nematode-associated Wolbachia, relative to strains that infect arthropods. An exception is filarial Wolbachia strain wMhie, assigned to supergroup F, which occurs in a nematode host from a poikilothermic lizard. Genes encoding LexA and error-prone polymerases are absent from all Wolbachia genomes, suggesting that the SOS functions induced by RecA-mediated activation of LexA do not occur, despite retention of genes encoding a few proteins that respond to LexA induction in E. coli. Three independent E. coli accessions converge on a single Wolbachia UvrD helicase, which interacts with mismatch repair proteins MutS and MutL, encoded in nearly all Wolbachia genomes. With the exception of MutL, which has been mapped to a eukaryotic association module in Phage WO, proteins involved in recombination/repair are uniformly represented by single protein annotations. Putative phage-encoded MutL proteins are restricted to Wolbachia supergroups A and B and show higher amino acid identity than chromosomally encoded MutL orthologs. This analysis underscores differences between nematode and arthropod-associated Wolbachia and describes aspects of DNA metabolism that potentially impact development of procedures for transformation and genetic manipulation of Wolbachia.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

Molecular diversity of Wolbachia in Lepidoptera: Prevalent allelic content and high recombination of MLST genes

Wolbachia are common endosymbiotic bacteria of Arthropoda and Nematoda that are ordinarily transmitted vertically in host lineages through the egg cytoplasm. Despite the great interest in the Wolbachia symbiont, many issues of its biology remain unclear, including its evolutionary history, routes of transfer among species, and the molecular mechanisms underlying the symbiont's effect on its host. In this report, we present data relating to Wolbachia infection in 120 species of 13 Lepidoptera families, mostly butterflies, from West Siberian localities based on Multilocus sequence typing (MLST) and the wsp locus and perform a comprehensive survey of the distribution of Wolbachia and its genetic diversity in Lepidoptera worldwide. We observed a high infection incidence in the studied region; this finding is probably also true for other temperate latitude regions because many studied species have broad Palearctic and even Holarctic distribution. Although 40 new MLST alleles and 31 new STs were described, there was no noticeable difference in the MLST allele content in butterflies and probably also in moths worldwide. A genetic analysis of Wolbachia strains revealed the MLST allele core in lepidopteran hosts worldwide, viz. the ST-41 allele content. The key finding of our study was the detection of rampant recombination among MLST haplotypes. High rates of homologous recombination between Wolbachia strains indicate a substantial contribution of genetic exchanges to the generation of new STs. This finding should be considered when discussing issues related to the reconstruction of Wolbachia evolution, divergence time, and the routes of Wolbachia transmission across arthropod hosts.

Supergroup C Wolbachia, mutualist symbionts of filarial nematodes, have a distinct genome structure

Wolbachia pipientis is possibly the most widespread endosymbiont of arthropods and nematodes. While all Wolbachia strains have historically been defined as a single species, 16 monophyletic clusters of diversity (called supergroups) have been described. Different supergroups have distinct host ranges and symbiotic relationships, ranging from mutualism to reproductive manipulation. In filarial nematodes, which include parasites responsible for major diseases of humans (such as Onchocerca volvulus, agent of river blindness) and companion animals (Dirofilaria immitis, the dog heartworm), Wolbachia has an obligate mutualist role and is the target of new treatment regimens. Here, we compare the genomes of eight Wolbachia strains, spanning the diversity of the major supergroups (A–F), analysing synteny, transposable element content, GC skew and gene loss or gain. We detected genomic features that differ between Wolbachia supergroups, most notably in the C and D clades from filarial nematodes. In particular, strains from supergroup C (symbionts of O. volvulus and D. immitis) present a pattern of GC skew, conserved synteny and lack of transposable elements, unique in the Wolbachia genus. These features could be the consequence of a distinct symbiotic relationship between C Wolbachia strains and their hosts, highlighting underappreciated differences between the mutualistic supergroups found within filarial nematodes.

Feminization of the Isopod Cylisticus convexus after Transinfection of the wVulC Wolbachia Strain of Armadillidium vulgare

Reproductive parasites such as Wolbachia are able to manipulate the reproduction of their hosts by inducing parthenogenesis, male-killing, cytoplasmic incompatibility or feminization of genetic males. Despite extensive studies, no underlying molecular mechanism has been described to date. The goal of this study was to establish a system with a single Wolbachia strain that feminizes two different isopod species to enable comparative analyses aimed at elucidating the genetic basis of feminization. It was previously suggested that Wolbachia wVulC, which naturally induces feminization in Armadillidium vulgare, induces the development of female secondary sexual characters in transinfected Cylisticus convexus adult males. However, this does not demonstrate that wVulC induces feminization in C. convexus since feminization is the conversion of genetic males into functional females that occurs during development. Nevertheless, it suggests that C. convexus may represent a feminization model suitable for further development. Knowledge about C. convexus sexual differentiation is also essential for comparative analyses, as feminization is thought to take place just before or during sexual differentiation. Consequently, we first described gonad morphological differentiation of C. convexus and compared it with that of A. vulgare. Then, wVulC was injected into male and female C. convexus adult individuals. The feminizing effect was demonstrated by the combined appearance of female secondary sexual characters in transinfected adult males, as well as the presence of intersexes and female biases in progenies in which wVulC was vertically transmitted from transinfected mothers. The establishment of a new model of feminization of a Wolbachia strain in a heterologous host constitutes a useful tool towards the understanding of the molecular mechanism of feminization.