Genome-Wide Patterns of Bracovirus Chromosomal Integration into Multiple Host Tissues during Parasitism.

Bracoviruses are domesticated viruses found in parasitic wasp genomes. They are composed of genes of nudiviral origin that are involved in particle production and proviral segments containing virulence genes that are necessary for parasitism success. During particle production, proviral segments are amplified and individually packaged as DNA circles in nucleocapsids. These particles are injected by parasitic wasps into host larvae together with their eggs. Bracovirus circles of two wasp species were reported to undergo chromosomal integration in parasitized host hemocytes, through a conserved sequence named the host integration motif (HIM). Here, we used bulk Illumina sequencing to survey integrations of Cotesia typhae bracovirus circles in the DNA of its host, the maize corn borer (Sesamia nonagrioides), 7 days after parasitism. First, assembly and annotation of a high-quality genome for C. typhae enabled us to characterize 27 proviral segments clustered in proviral loci. Using these data, we characterized large numbers of chromosomal integrations (from 12 to 85 events per host haploid genome) for all 16 bracovirus circles containing a HIM. Integrations were found in four S. nonagrioides tissues and in the body of a caterpillar in which parasitism had failed. The 12 remaining circles do not integrate but are maintained at high levels in host tissues. Surprisingly, we found that HIM-mediated chromosomal integration in the wasp germ line has occurred accidentally at least six times during evolution. Overall, our study furthers our understanding of wasp-host genome interactions and supports HIM-mediated chromosomal integration as a possible mechanism of horizontal transfer from wasps to their hosts. IMPORTANCE Bracoviruses are endogenous domesticated viruses of parasitoid wasps that are injected together with wasp eggs into wasp host larvae during parasitism. Several studies have shown that some DNA circles packaged into bracovirus particles become integrated into host somatic genomes during parasitism, but the phenomenon has never been studied using nontargeted approaches. Here, we use bulk Illumina sequencing to systematically characterize and quantify bracovirus circle integrations that occur in four tissues of the Mediterranean corn borer (Sesamia nonagrioides) during parasitism by the Cotesia typhae wasp. Our analysis reveals that all circles containing a HIM integrate at substantial levels (from 12 to 85 integrations per host cell, in total) in all tissues, while other circles do not integrate. In addition to shedding new light on wasp-bracovirus-host interactions, our study supports HIM-mediated chromosomal integration of bracovirus as a possible source of wasp-to-host horizontal transfer, with long-term evolutionary consequences.

Characterization of a Sex-Determining Region and Its Genomic Context via Statistical Estimates of Haplotype Frequencies in Daughters and Sons Sequenced in Pools.

Sex chromosomes are generally derived from a pair of autosomes that have acquired a locus controlling sex. Sex chromosomes may evolve reduced recombination around this locus and undergo a long process of molecular divergence. At that point, the original loci controlling sex may be difficult to pinpoint. This difficulty has affected many model species from mammals to birds to flies, which present highly diverged sex chromosomes. Identifying sex-controlling loci is easier in species with molecularly similar sex chromosomes. Here we aimed at pinpointing the sex-determining region (SDR) of Armadillidium vulgare, a terrestrial isopod with female heterogamety (ZW females and ZZ males) and whose sex chromosomes appear to show low genetic divergence. To locate the SDR, we assessed single-nucleotide polymorphism (SNP) allele frequencies in F1 daughters and sons sequenced in pools (pool-seq) in several families. We developed a Bayesian method that uses the SNP genotypes of individually sequenced parents and pool-seq data from F1 siblings to estimate the genetic distance between a given genomic region (contig) and the SDR. This allowed us to assign more than 43 Mb of contigs to sex chromosomes, and to demonstrate extensive recombination and very low divergence between these chromosomes. By taking advantage of multiple F1 families, we delineated a very short genomic region (∼65 kb) that presented no evidence of recombination with the SDR. In this short genomic region, the comparison of sequencing depths between sexes highlighted female-specific genes that have undergone recent duplication, and which may be involved in sex determination in A. vulgare.

Comparative Genomics of Strictly Vertically Transmitted, Feminizing Microsporidia Endosymbionts of Amphipod Crustaceans.

Microsporidia are obligate intracellular eukaryotic parasites of vertebrates and invertebrates. Microsporidia are usually pathogenic and undergo horizontal transmission or a mix of horizontal and vertical transmission. However, cases of nonpathogenic microsporidia, strictly vertically transmitted from mother to offspring, have been reported in amphipod crustaceans. Some of them further evolved the ability to feminize their nontransmitting male hosts into transmitting females. However, our understanding of the evolution of feminization in microsporidia is hindered by a lack of genomic resources. We report the sequencing and analysis of three strictly vertically transmitted microsporidia species for which feminization induction has been demonstrated (Nosema granulosis) or is strongly suspected (Dictyocoela muelleri and Dictyocoela roeselum), along with a draft genome assembly of their host Gammarus roeselii. Contrary to horizontally transmitted microsporidia that form environmental spores that can be purified, feminizing microsporidia cannot be easily isolated from their host cells. Therefore, we cosequenced symbiont and host genomic DNA and devised a computational strategy to obtain genome assemblies for the different partners. Genomic comparison with feminizing Wolbachia bacterial endosymbionts of isopod crustaceans indicated independent evolution of feminization in microsporidia and Wolbachia at the molecular genetic level. Feminization thus represents a remarkable evolutionary convergence of eukaryotic and prokaryotic microorganisms. Furthermore, a comparative genomics analysis of microsporidia allowed us to identify several candidate genes for feminization, involving functions such as DNA binding and membrane fusion. The genomic resources we generated contribute to establish Gammarus roeselii and its microsporidia symbionts as a new model to study the evolution of symbiont-mediated feminization.

Sex chromosomes control vertical transmission of feminizing Wolbachia symbionts in an isopod.

Microbial endosymbiosis is widespread in animals, with major ecological and evolutionary implications. Successful symbiosis relies on efficient vertical transmission through host generations. However, when symbionts negatively affect host fitness, hosts are expected to evolve suppression of symbiont effects or transmission. Here, we show that sex chromosomes control vertical transmission of feminizing Wolbachia endosymbionts in the isopod Armadillidium nasatum. Theory predicts that the invasion of an XY/XX species by cytoplasmic sex ratio distorters is unlikely because it leads to fixation of the unusual (and often lethal or infertile) YY genotype. We demonstrate that A. nasatum X and Y sex chromosomes are genetically highly similar and that YY individuals are viable and fertile, thereby enabling Wolbachia spread in this XY-XX species. Nevertheless, we show that Wolbachia cannot drive fixation of YY individuals, because infected YY females do not transmit Wolbachia to their offspring, unlike XX and XY females. The genetic basis fits the model of a Y-linked recessive allele (associated with an X-linked dominant allele), in which the homozygous state suppresses Wolbachia transmission. Moreover, production of all-male progenies by infected YY females restores a balanced sex ratio at the host population level. This suggests that blocking of Wolbachia transmission by YY females may have evolved to suppress feminization, thereby offering a whole new perspective on the evolutionary interplay between microbial symbionts and host sex chromosomes.

The Genome of Armadillidium vulgare (Crustacea, Isopoda) Provides Insights into Sex Chromosome Evolution in the Context of Cytoplasmic Sex Determination.

The terrestrial isopod Armadillidium vulgare is an original model to study the evolution of sex determination and symbiosis in animals. Its sex can be determined by ZW sex chromosomes, or by feminizing Wolbachia bacterial endosymbionts. Here, we report the sequence and analysis of the ZW female genome of A. vulgare. A distinguishing feature of the 1.72 gigabase assembly is the abundance of repeats (68% of the genome). We show that the Z and W sex chromosomes are essentially undifferentiated at the molecular level and the W-specific region is extremely small (at most several hundreds of kilobases). Our results suggest that recombination suppression has not spread very far from the sex-determining locus, if at all. This is consistent with A. vulgare possessing evolutionarily young sex chromosomes. We characterized multiple Wolbachia nuclear inserts in the A. vulgare genome, none of which is associated with the W-specific region. We also identified several candidate genes that may be involved in the sex determination or sexual differentiation pathways. The A. vulgare genome serves as a resource for studying the biology and evolution of crustaceans, one of the most speciose and emblematic metazoan groups.

Untangling Heteroplasmy, Structure, and Evolution of an Atypical Mitochondrial Genome by PacBio Sequencing.

The highly compact mitochondrial (mt) genome of terrestrial isopods (Oniscidae) presents two unusual features. First, several loci can individually encode two tRNAs, thanks to single nucleotide polymorphisms at anticodon sites. Within-individual variation (heteroplasmy) at these loci is thought to have been maintained for millions of years because individuals that do not carry all tRNA genes die, resulting in strong balancing selection. Second, the oniscid mtDNA genome comes in two conformations: a ∼14 kb linear monomer and a ∼28 kb circular dimer comprising two monomer units fused in palindrome. We hypothesized that heteroplasmy actually results from two genome units of the same dimeric molecule carrying different tRNA genes at mirrored loci. This hypothesis, however, contradicts the earlier proposition that dimeric molecules result from the replication of linear monomers-a process that should yield totally identical genome units within a dimer. To solve this contradiction, we used the SMRT (PacBio) technology to sequence mirrored tRNA loci in single dimeric molecules. We show that dimers do present different tRNA genes at mirrored loci; thus covalent linkage, rather than balancing selection, maintains vital variation at anticodons. We also leveraged unique features of the SMRT technology to detect linear monomers closed by hairpins and carrying noncomplementary bases at anticodons. These molecules contain the necessary information to encode two tRNAs at the same locus, and suggest new mechanisms of transition between linear and circular mtDNA. Overall, our analyses clarify the evolution of an atypical mt genome where dimerization counterintuitively enabled further mtDNA compaction.

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.

Comparative paleovirological analysis of crustaceans identifies multiple widespread viral groups.

BACKGROUND: The discovery of many fragments of viral genomes integrated in the genome of their eukaryotic host (endogenous viral elements; EVEs) has recently opened new avenues to further our understanding of viral evolution and of host-virus interactions. Here, we report the results of a comprehensive screen for EVEs in crustaceans. Following up on the recent discovery of EVEs in the terrestrial isopod, Armadillidium vulgare, we scanned the genomes of six crustacean species: a terrestrial isopod (Armadillidium nasatum), two water fleas (Daphnia pulex and D. pulicaria), two copepods (the salmon louse, Lepeophtheirus salmonis and Eurytemora affinis), and a freshwater amphipod (Hyalella azteca). RESULTS: In total, we found 210 EVEs representing 14 different lineages belonging to five different viral groups that are present in two to five species: Bunyaviridae (-ssRNA), Circoviridae (ssDNA), Mononegavirales (-ssRNA), Parvoviridae (ssDNA) and Totiviridae (dsRNA). The identification of shared orthologous insertions between A. nasatum and A. vulgare indicates that EVEs have been maintained over several millions of years, although we did not find any evidence supporting exaptation. Overall, the different degrees of EVE degradation (from none to >10 nonsense mutations) suggest that endogenization has been recurrent during the evolution of the various crustacean taxa. Our study is the first to report EVEs in D. pulicaria, E. affinis and H. azteca, many of which are likely to result from recent endogenization of currently circulating viruses. CONCLUSIONS: In conclusion, we have unearthed a large diversity of EVEs from crustacean genomes, and shown that four of the five viral groups we uncovered (Bunyaviridae, Circoviridae, Mononegavirales, Parvoviridae) were and may still be present in three to four highly divergent crustacean taxa. In addition, the discovery of recent EVEs offers an interesting opportunity to characterize new exogenous viruses currently circulating in economically or ecologically important copepod species.