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.

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.

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.

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.

Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization.

Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction.

Convergent consequences of parthenogenesis on stick insect genomes.

The shift from sexual reproduction to parthenogenesis has occurred repeatedly in animals, but how the loss of sex affects genome evolution remains poorly understood. We generated reference genomes for five independently evolved parthenogenetic species in the stick insect genus Timema and their closest sexual relatives. Using these references and population genomic data, we show that parthenogenesis results in an extreme reduction of heterozygosity and often leads to genetically uniform populations. We also find evidence for less effective positive selection in parthenogenetic species, suggesting that sex is ubiquitous in natural populations because it facilitates fast rates of adaptation. Parthenogenetic species did not show increased transposable element (TE) accumulation, likely because there is little TE activity in the genus. By using replicated sexual-parthenogenetic comparisons, our study reveals how the absence of sex affects genome evolution in natural populations, providing empirical support for the negative consequences of parthenogenesis as predicted by theory.

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.

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.

First annotated draft genomes of nonmarine ostracods (Ostracoda, Crustacea) with different reproductive modes.

Ostracods are one of the oldest crustacean groups with an excellent fossil record and high importance for phylogenetic analyses but genome resources for this class are still lacking. We have successfully assembled and annotated the first reference genomes for three species of nonmarine ostracods; two with obligate sexual reproduction (Cyprideis torosa and Notodromas monacha) and the putative ancient asexual Darwinula stevensoni. This kind of genomic research has so far been impeded by the small size of most ostracods and the absence of genetic resources such as linkage maps or BAC libraries that were available for other crustaceans. For genome assembly, we used an Illumina-based sequencing technology, resulting in assemblies of similar sizes for the three species (335-382 Mb) and with scaffold numbers and their N50 (19-56 kb) in the same orders of magnitude. Gene annotations were guided by transcriptome data from each species. The three assemblies are relatively complete with BUSCO scores of 92-96. The number of predicted genes (13,771-17,776) is in the same range as Branchiopoda genomes but lower than in most malacostracan genomes. These three reference genomes from nonmarine ostracods provide the urgently needed basis to further develop ostracods as models for evolutionary and ecological research.

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.

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.

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.

Sex-biased gene expression is repeatedly masculinized in asexual females.

Males and females feature strikingly different phenotypes, despite sharing most of their genome. A resolution of this apparent paradox is through differential gene expression, whereby genes are expressed at different levels in each sex. This resolution, however, is likely to be incomplete, leading to conflict between males and females over the optimal expression of genes. Here we test the hypothesis that gene expression in females is constrained from evolving to its optimum level due to sexually antagonistic selection on males, by examining changes in sex-biased gene expression in five obligate asexual species of stick insect, which do not produce males. We predicted that the transcriptome of asexual females would be feminized as asexual females do not experience any sexual conflict. Contrary to our prediction we find that asexual females feature masculinized gene expression, and hypothesise that this is due to shifts in female optimal gene expression levels following the suppression of sex.

Mapping of Multiple Complementary Sex Determination Loci in a Parasitoid Wasp.

Sex determination has evolved in a variety of ways and can depend on environmental and genetic signals. A widespread form of genetic sex determination is haplodiploidy, where unfertilized, haploid eggs develop into males and fertilized diploid eggs into females. One of the molecular mechanisms underlying haplodiploidy in Hymenoptera, the large insect order comprising ants, bees, and wasps, is complementary sex determination (CSD). In species with CSD, heterozygosity at one or several loci induces female development. Here, we identify the genomic regions putatively underlying multilocus CSD in the parasitoid wasp Lysiphlebus fabarum using restriction-site associated DNA sequencing. By analyzing segregation patterns at polymorphic sites among 331 diploid males and females, we identify up to four CSD candidate regions, all on different chromosomes. None of the candidate regions feature evidence for homology with the csd gene from the honey bee, the only species in which CSD has been characterized, suggesting that CSD in L. fabarum is regulated via a novel molecular mechanism. Moreover, no homology is shared between the candidate loci, in contrast to the idea that multilocus CSD should emerge from duplications of an ancestral single-locus system. Taken together, our results suggest that the molecular mechanisms underlying CSD in Hymenoptera are not conserved between species, raising the question as to whether CSD may have evolved multiple times independently in the group.

Asexual reproduction reduces transposable element load in experimental yeast populations.

Theory predicts that sexual reproduction can either facilitate or restrain transposable element (TE) accumulation by providing TEs with a means of spreading to all individuals in a population, versus facilitating TE load reduction via purifying selection. By quantifying genomic TE loads over time in experimental sexual and asexual Saccharomyces cerevisiae populations, we provide direct evidence that TE loads decrease rapidly under asexual reproduction. We show, using simulations, that this reduction may occur via evolution of TE activity, most likely via increased excision rates. Thus, sex is a major driver of genomic TE loads and at the root of the success of TEs.

Hybridogenesis.

Hybridogenesis is an unusual form of reproduction that is found in hybrids between different species. It involves the selective transmission of one of the parental genomes, while the other one is renewed by mating with the corresponding species. It is seen as a form of sexual parasitism, in which the hybridogenetic genome gains a twofold transmission advantage and exploits the reproductive effort of another species.

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.

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.