Haplotype divergence supports long-term asexuality in the oribatid mite Oppiella nova.

Sex strongly impacts genome evolution via recombination and segregation. In the absence of these processes, haplotypes within lineages of diploid organisms are predicted to accumulate mutations independently of each other and diverge over time. This so-called "Meselson effect" is regarded as a strong indicator of the long-term evolution under obligate asexuality. Here, we present genomic and transcriptomic data of three populations of the asexual oribatid mite species Oppiella nova and its sexual relative Oppiella subpectinata We document strikingly different patterns of haplotype divergence between the two species, strongly supporting Meselson effect-like evolution and long-term asexuality in O. nova: I) variation within individuals exceeds variation between populations in O. nova but vice versa in O. subpectinata; II) two O. nova sublineages feature a high proportion of lineage-specific heterozygous single-nucleotide polymorphisms (SNPs), indicating that haplotypes continued to diverge after lineage separation; III) the deepest split in gene trees generally separates the two haplotypes in O. nova, but populations in O. subpectinata; and IV) the topologies of the two haplotype trees match each other. Our findings provide positive evidence for the absence of canonical sex over evolutionary time in O. nova and suggest that asexual oribatid mites can escape the dead-end fate usually associated with asexual lineages.

What Ecological Factors Favor Parthenogenesis over Sexual Reproduction? A Study on the Facultatively Parthenogenetic Mayfly Alainites muticus in Natural Populations.

AbstractDifferent reproductive modes are characterized by costs and benefits that depend on ecological contexts. For example, sex can provide benefits under complex biotic interactions, while its costs increase under mate limitation. Furthermore, ecological contexts often vary along abiotic gradients. Here, we study how these factors simultaneously influence the frequency of sex in the facultatively parthenogenetic mayfly Alainites muticus. We first verified that parthenogenesis translates into female-biased population sex ratios. We then measured the density of individuals (a proxy for mate limitation) and community diversity (biotic interaction complexity) for 159 A. muticus populations covering a broad altitudinal gradient and used structural equation modeling to investigate their direct and indirect influences on sex ratios. We found no effect of community diversity or altitude on sex ratios. Furthermore, even when females can reproduce parthenogenetically, they generally reproduce sexually, indicating that the benefits of sex exceed its costs in most situations. Sex ratios become female-biased only under low population densities, as expected if mate limitation was the main factor selecting for parthenogenesis. Mate limitation might be widespread in mayflies because of their short adult life span and limited dispersal, which can generate strong selection for reproductive assurance and may provide a stepping stone toward obligate parthenogenesis.

Evidence for cryptic sex in parthenogenetic stick insects of the genus Timema.

Obligately parthenogenetic species are expected to be short lived since the lack of sex and recombination should translate into a slower adaptation rate and increased accumulation of deleterious alleles. Some, however, are thought to have been reproducing without males for millions of years. It is not clear how these old parthenogens can escape the predicted long-term costs of parthenogenesis, but an obvious explanation is cryptic sex. In this study, we screen for signatures of cryptic sex in eight populations of four parthenogenetic species of Timema stick insects, some estimated to be older than 1 Myr. Low genotype diversity, homozygosity of individuals and high linkage disequilibrium (LD) unaffected by marker distances support exclusively parthenogenetic reproduction in six populations. However, in two populations (namely, of the species Timema douglasi and T. monikensis) we find strong evidence for cryptic sex, most likely mediated by rare males. These populations had comparatively high genotype diversities, lower LD, and a clear LD decay with genetic distance. Rare sex in species that are otherwise largely parthenogenetic could help explain the unusual success of parthenogenesis in the Timema genus and raises the question whether episodes of rare sex are in fact the simplest explanation for the persistence of many old parthenogens in nature.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species.

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Dynamics of sex-biased gene expression during development in the stick insect Timema californicum.

Sexually dimorphic phenotypes are thought to arise primarily from sex-biased gene expression during development. Major changes in developmental strategies, such as the shift from hemimetabolous to holometabolous development, are therefore expected to have profound consequences for the dynamics of sex-biased gene expression. However, no studies have previously examined sex-biased gene expression during development in hemimetabolous insects, precluding comparisons between developmental strategies. Here we characterized sex-biased gene expression at three developmental stages in a hemimetabolous stick insect (Timema californicum): hatchlings, juveniles, and adults. As expected, the proportion of sex-biased genes gradually increased during development, mirroring the gradual increase of phenotypic sexual dimorphism. Sex-biased genes identified at early developmental stages were generally consistently male- or female-biased at later stages, suggesting their importance in sexual differentiation. Additionally, we compared the dynamics of sex-biased gene expression during development in T. californicum to those of the holometabolous fly Drosophila melanogaster by reanalyzing publicly available RNA-seq data from third instar larval, pupal and adult stages. In D. melanogaster, 84% of genes were sex-biased at the adult stage (compared to only 20% in T. californicum), and sex-biased gene expression increased abruptly at the adult stage when morphological sexual dimorphism is manifested. Our findings are consistent with the prediction that the dynamics of sex-biased gene expression during development differ extensively between holometabolous and hemimetabolous insect species.

A single QTL with large effect is associated with female functional virginity in an asexual parasitoid wasp.

During the transition from sexual to asexual reproduction, a suite of reproduction-related sexual traits become superfluous, and may be selected against if costly. Female functional virginity refers to asexual females resisting to mate or not fertilizing eggs after mating. These traits appear to be among the first that evolve during transitions from sexual to asexual reproduction. The genetic basis of female functional virginity remains elusive. Previously, we reported that female functional virginity segregates as expected for a single recessive locus in the asexual parasitoid wasp Asobara japonica. Here, we investigate the genetic basis of this trait by quantitative trait loci (QTL) mapping and candidate gene analyses. Consistent with the segregation of phenotypes, we found a single QTL of large effect, spanning over 4.23 Mb and comprising at least 131 protein-coding genes, of which 15 featured sex-biased expression in the related sexual species Asobara tabida. Two of the 15 sex-biased genes were previously identified to differ between related sexual and asexual population/species: CD151 antigen and nuclear pore complex protein Nup50. A third gene, hormone receptor 4, is involved in steroid hormone mediated mating behaviour. Overall, our results are consistent with a single locus, or a cluster of closely linked loci, underlying rapid evolution of female functional virginity in the transition to asexuality. Once this variant, causing rejection to mate, has swept through a population, the flanking region does not get smaller owing to lack of recombination in asexuals.

Extremely Widespread Parthenogenesis and a Trade-Off Between Alternative Forms of Reproduction in Mayflies (Ephemeroptera).

Studying alternative forms of reproduction in natural populations is of fundamental importance for understanding the costs and benefits of sex. Mayflies are one of the few animal groups where sexual reproduction co-occurs with different types of parthenogenesis, providing ideal conditions for identifying benefits of sex in natural populations. Here, we establish a catalog of all known mayfly species capable of reproducing by parthenogenesis, as well as species unable to do so. Overall, 1.8% of the described species reproduce parthenogenetically, which is an order of magnitude higher than reported in other animal groups. This frequency even reaches 47.8% if estimates are based on the number of studied rather than described mayfly species, as reproductive modes have thus far been studied in only 17 out of 42 families. We find that sex is a more successful strategy than parthenogenesis (associated with a higher hatching success of eggs), with a trade-off between the hatching success of parthenogenetic and sexual eggs. This means that improving the capacity for parthenogenesis may come at a cost for sexual reproduction. Such a trade-off can help explain why facultative parthenogenesis is extremely rare among animals despite its potential to combine the benefits of sexual and parthenogenetic reproduction. We argue that parthenogenesis is frequently selected in mayflies in spite of this probable trade-off because their typically low dispersal ability and short and fragile adult life may frequently generate situations of mate limitation in females. Mayflies are currently clearly underappreciated for understanding the benefits of sex under natural conditions.

Genomic Features of Parthenogenetic Animals.

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum.

BACKGROUND: Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts. RESULTS: We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes. CONCLUSIONS: These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Repeated Evolution of Asexuality Involves Convergent Gene Expression Changes.

Asexual reproduction has evolved repeatedly from sexual ancestors across a wide range of taxa. Whereas the costs and benefits associated with asexuality have received considerable attention, the molecular changes underpinning the evolution of asexual reproduction remain relatively unexplored. In particular, it is completely unknown whether the repeated evolution of asexual phenotypes involves similar molecular changes, as previous studies have focused on changes occurring in single lineages. Here, we investigate the extent of convergent gene expression changes across five independent transitions to asexuality in stick insects. We compared gene expression of asexual females to females of close sexual relatives in whole-bodies, reproductive tracts, and legs. We identified a striking amount of convergent gene expression change (up to 8% of genes), greatly exceeding that expected by chance. Convergent changes were also tissue-specific, and most likely driven by selection for functional changes. Genes showing convergent changes in the reproductive tract were associated with meiotic spindle formation and centrosome organization. These genes are particularly interesting as they can influence the production of unreduced eggs, a key barrier to asexual reproduction. Changes in legs and whole-bodies were likely involved in female sexual trait decay, with enrichment in terms such as sperm-storage and pigmentation. By identifying changes occurring across multiple independent transitions to asexuality, our results provide a rare insight into the molecular basis of asexual phenotypes and suggest that the evolutionary path to asexuality is highly constrained, requiring repeated changes to the same key genes.

Evolutionary dynamics of specialisation in herbivorous stick insects.

Understanding the evolutionary dynamics underlying herbivorous insect mega-diversity requires investigating the ability of insects to shift and adapt to different host plants. Feeding experiments with nine related stick insect species revealed that insects retain the ability to use ancestral host plants after shifting to novel hosts, with host plant shifts generating fundamental feeding niche expansions. These expansions were, however, not accompanied by expansions of the realised feeding niches, as species on novel hosts are generally ecologically specialised. For shifts from angiosperm to chemically challenging conifer hosts, generalist fundamental feeding niches even evolved jointly with strong host plant specialisation, indicating that host plant specialisation is not driven by constraints imposed by plant chemistry. By coupling analyses of plant chemical compounds, fundamental and ecological feeding niches in multiple insect species, we provide novel insights into the evolutionary dynamics of host range expansion and contraction in herbivorous insects.

Consequences of Asexuality in Natural Populations: Insights from Stick Insects.

Recombination is a fundamental process with significant impacts on genome evolution. Predicted consequences of the loss of recombination include a reduced effectiveness of selection, changes in the amount of neutral polymorphisms segregating in populations, and an arrest of GC-biased gene conversion. Although these consequences are empirically well documented for nonrecombining genome portions, it remains largely unknown if they extend to the whole genome scale in asexual organisms. We identify the consequences of asexuality using de novo transcriptomes of five independently derived, obligately asexual lineages of stick insects, and their sexual sister-species. We find strong evidence for higher rates of deleterious mutation accumulation, lower levels of segregating polymorphisms and arrested GC-biased gene conversion in asexuals as compared with sexuals. Taken together, our study conclusively shows that predicted consequences of genome evolution under asexuality can indeed be found in natural populations.

Niche differentiation among clones in asexual grass thrips.

Many asexual animal populations comprise a mixture of genetically different lineages, but to what degree this genetic diversity leads to ecological differences remains often unknown. Here, we test whether genetically different clonal lineages of Aptinothrips grass thrips differ in performance on a range of plants used as hosts in natural populations. We find a clear clone-by-plant species interactive effect on reproductive output, meaning that clonal lineages perform differently on different plant species and thus are characterized by disparate ecological niches. This implies that local clonal diversities can be driven and maintained by frequency-dependent selection and that resource heterogeneity can generate diverse clone assemblies.

Fundamental and realized feeding niche breadths of sexual and asexual stick insects.

The factors contributing to the maintenance of sex over asexuality in natural populations remain unclear. Ecological divergences between sexual and asexual lineages could help to maintain reproductive polymorphisms, at least transiently, but the consequences of asexuality for the evolution of ecological niches are unknown. Here, we investigated how niche breadths change in transitions from sexual reproduction to asexuality. We used host plant ranges as a proxy to compare the realized feeding niche breadths of five independently derived asexual Timema stick insect species and their sexual relatives at both the species and population levels. Asexual species had systematically narrower realized niches than sexual species, though this pattern was not apparent at the population level. To investigate how the narrower realized niches of asexual species arise, we performed feeding experiments to estimate fundamental niche breadths but found no systematic differences between reproductive modes. The narrow realized niches found in asexual species are therefore probably a consequence of biotic interactions such as predation or competition, that constrain realized niche size in asexuals more strongly than in sexuals.

No Accumulation of Transposable Elements in Asexual Arthropods.

Transposable elements (TEs) and other repetitive DNA can accumulate in the absence of recombination, a process contributing to the degeneration of Y-chromosomes and other nonrecombining genome portions. A similar accumulation of repetitive DNA is expected for asexually reproducing species, given their entire genome is effectively nonrecombining. We tested this expectation by comparing the whole-genome TE loads of five asexual arthropod lineages and their sexual relatives, including asexual and sexual lineages of crustaceans (Daphnia water fleas), insects (Leptopilina wasps), and mites (Oribatida). Surprisingly, there was no evidence for increased TE load in genomes of asexual as compared to sexual lineages, neither for all classes of repetitive elements combined nor for specific TE families. Our study therefore suggests that nonrecombining genomes do not accumulate TEs like nonrecombining genomic regions of sexual lineages. Even if a slight but undetected increase of TEs were caused by asexual reproduction, it appears to be negligible compared to variance between species caused by processes unrelated to reproductive mode. It remains to be determined if molecular mechanisms underlying genome regulation in asexuals hamper TE activity. Alternatively, the differences in TE dynamics between nonrecombining genomes in asexual lineages versus nonrecombining genome portions in sexual species might stem from selection for benign TEs in asexual lineages because of the lack of genetic conflict between TEs and their hosts and/or because asexual lineages may only arise from sexual ancestors with particularly low TE loads.

Diploid males support a two-step mechanism of endosymbiont-induced thelytoky in a parasitoid wasp.

BACKGROUND: Haplodiploidy, where females develop from diploid, fertilized eggs and males from haploid, unfertilized eggs, is abundant in some insect lineages. Some species in these lineages reproduce by thelytoky that is caused by infection with endosymbionts: infected females lay haploid eggs that undergo diploidization and develop into females, while males are very rare or absent. It is generally assumed that in thelytokous wasps, endosymbionts merely diploidize the unfertilized eggs, which would then trigger female development. RESULTS: We found that females in the parasitoid wasp Asobara japonica infected with thelytoky-inducing Wolbachia produce 0.7-1.2% male offspring. Seven to 39% of these males are diploid, indicating that diploidization and female development can be uncoupled in A. japonica. Wolbachia titer in adults was correlated with their ploidy and sex: diploids carried much higher Wolbachia titers than haploids, and diploid females carried more Wolbachia than diploid males. Data from introgression lines indicated that the development of diploid individuals into males instead of females is not caused by malfunction-mutations in the host genome but that diploid males are most likely produced when the endosymbiont fails to activate the female sex determination pathway. Our data therefore support a two-step mechanism by which endosymbionts induce thelytoky in A. japonica: diploidization of the unfertilized egg is followed by feminization, whereby each step correlates with a threshold of endosymbiont titer during wasp development. CONCLUSIONS: Our new model of endosymbiont-induced thelytoky overthrows the view that certain sex determination mechanisms constrain the evolution of endosymbiont-induced thelytoky in hymenopteran insects. Endosymbionts can cause parthenogenesis through feminization, even in groups in which endosymbiont-diploidized eggs would develop into males following the hosts' sex determination mechanism. In addition, our model broadens our understanding of the mechanisms by which endosymbionts induce thelytoky to enhance their transmission to the next generation. Importantly, it also provides a novel window to study the yet-poorly known haplodiploid sex determination mechanisms in haplodiploid insects.

Female Stick Insects Mate Multiply to Find Compatible Mates.

Why females of many species mate multiply in the absence of direct benefits remains an open question in evolutionary ecology. Interacting and mating with multiple males can be costly to females in terms of time, resources, predation risk, and disease transmission. A number of indirect genetic benefits have been proposed to explain such behaviors, but the relative importance of these mechanisms in natural systems remains unclear. We tested for several direct and indirect benefits of polyandry in the walking stick Timema cristinae. We found no evidence of direct benefits with respect to longevity or fecundity. However, male × female genotypic interactions affected egg-hatching success and offspring production independent of relatedness, suggesting that mating with certain males benefits females and that the best male may differ for each female. Furthermore, multiply mated females biased paternity toward one or few males, and the extent of this bias was positively correlated to egg-hatching success. Our data, therefore, provide evidence for indirect benefits through compatibility effects in this species. By mating multiply, females may improve their chances of mating with a compatible male if compatibility cannot be assessed before mating. Such compatibility effects can explain the evolution and maintenance of polyandry in Timema and many other species.

Evolution of alternative reproductive systems in Bacillus stick insects.

Reproduction is a key feature of all organisms, yet the way in which it is achieved varies greatly across the tree of life. One striking example of this variation is the stick insect genus Bacillus, in which five different reproductive modes have been described: sex, facultative and obligate parthenogenesis, and two highly unusual reproductive modes: hybridogenesis and androgenesis. Under hybridogenesis, the entire genome from the paternal species is eliminated and replaced each generation by mating with the corresponding species. Under androgenesis, an egg is fertilized, but the developing diploid offspring bear two paternal genomes and no maternal genome, as a consequence of unknown mechanisms. Here, we reevaluate the previous descriptions of Bacillus lineages and the proposed F1 hybrid ancestries of the hybridogenetic and obligately parthenogenetic lineages (based on allozymes and karyotypes) from Sicily, where all these reproductive modes are found. We generate a chromosome-level genome assembly for a facultative parthenogenetic species (B. rossius) and combine extensive field sampling with RADseq and mtDNA data. We identify and genetically corroborate all previously described species and confirm the ancestry of hybrid lineages. All hybrid lineages have fully retained their F1 hybrid constitution throughout the genome, indicating that the elimination of the paternal genome in hybridogens is always complete and that obligate parthenogenesis in Bacillus hybrid species is not associated with an erosion of heterozygosity as known in other hybrid asexuals. Our results provide a stepping stone toward understanding the transitions between reproductive modes and the proximate mechanisms of genome elimination.

Hydrocarbon divergence and reproductive isolation in Timema stick insects.

BACKGROUND: Individuals commonly prefer certain trait values over others when choosing their mates. If such preferences diverge between populations, they can generate behavioral reproductive isolation and thereby contribute to speciation. Reproductive isolation in insects often involves chemical communication, and cuticular hydrocarbons, in particular, serve as mate recognition signals in many species. We combined data on female cuticular hydrocarbons, interspecific mating propensity, and phylogenetics to evaluate the role of cuticular hydrocarbons in diversification of Timema walking-sticks. RESULTS: Hydrocarbon profiles differed substantially among the nine analyzed species, as well as between partially reproductively-isolated T. cristinae populations adapted to different host plants. In no-choice trials, mating was more likely between species with similar than divergent hydrocarbon profiles, even after correcting for genetic divergences. The macroevolution of hydrocarbon profiles, along a Timema species phylogeny, fits best with a punctuated model of phenotypic change concentrated around speciation events, consistent with change driven by selection during the evolution of reproductive isolation. CONCLUSION: Altogether, our data indicate that cuticular hydrocarbon profiles vary among Timema species and populations, and that most evolutionary change in hydrocarbon profiles occurs in association with speciation events. Similarities in hydrocarbon profiles between species are correlated with interspecific mating propensities, suggesting a role for cuticular hydrocarbon profiles in mate choice and speciation in the genus Timema.

Deleterious mutation accumulation in asexual Timema stick insects.

Sexual reproduction is extremely widespread in spite of its presumed costs relative to asexual reproduction, indicating that it must provide significant advantages. One postulated benefit of sex and recombination is that they facilitate the purging of mildly deleterious mutations, which would accumulate in asexual lineages and contribute to their short evolutionary life span. To test this prediction, we estimated the accumulation rate of coding (nonsynonymous) mutations, which are expected to be deleterious, in parts of one mitochondrial (COI) and two nuclear (Actin and Hsp70) genes in six independently derived asexual lineages and related sexual species of Timema stick insects. We found signatures of increased coding mutation accumulation in all six asexual Timema and for each of the three analyzed genes, with 3.6- to 13.4-fold higher rates in the asexuals as compared with the sexuals. In addition, because coding mutations in the asexuals often resulted in considerable hydrophobicity changes at the concerned amino acid positions, coding mutations in the asexuals are likely associated with more strongly deleterious effects than in the sexuals. Our results demonstrate that deleterious mutation accumulation can differentially affect sexual and asexual lineages and support the idea that deleterious mutation accumulation plays an important role in limiting the long-term persistence of all-female lineages.