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 re-evaluate 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 towards 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.

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.

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.

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.

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.

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.

Habitat heterogeneity favors asexual reproduction in natural populations of grassthrips

Explaining the overwhelming success of sex among eukaryotes is difficult given the obvious costs of sex relative to asexuality. Different studies have shown that sex can provide benefits in spatially heterogeneous environments under specific conditions, but whether spatial heterogeneity commonly contributes to the maintenance of sex in natural populations remains unknown. We experimentally manipulated habitat heterogeneity for sexual and asexual thrips lineages in natural populations and under seminatural mesocosm conditions by varying the number of hostplants available to these herbivorous insects. Asexual lineages rapidly replaced the sexual ones, independently of the level of habitat heterogeneity in mesocosms. In natural populations, the success of sexual thrips decreased with increasing habitat heterogeneity, with sexual thrips apparently only persisting in certain types of hostplant communities. Our results illustrate how genetic diversity-based mechanisms can favor asexuality instead of sex when sexual lineages co-occur with genetically variable asexual lineages.

Extreme genetic diversity in asexual grass thrips populations

The continuous generation of genetic variation has been proposed as one of the main factors explaining the maintenance of sexual reproduction in nature. However, populations of asexual individuals may attain high levels of genetic diversity through within-lineage diversification, replicate transitions to asexuality from sexual ancestors and migration. How these mechanisms affect genetic variation in populations of closely related sexual and asexual taxa can therefore provide insights into the role of genetic diversity for the maintenance of sexual reproduction. Here, we evaluate patterns of intra- and interpopulation genetic diversity in sexual and asexual populations of Aptinothrips rufus grass thrips. Asexual A. rufus populations are found throughout the world, whereas sexual populations appear to be confined to few locations in the Mediterranean region. We found that asexual A. rufus populations are characterized by extremely high levels of genetic diversity, both in comparison with their sexual relatives and in comparison with other asexual species. Migration is extensive among asexual populations over large geographic distances, whereas close sexual populations are strongly isolated from each other. The combination of extensive migration with replicate evolution of asexual lineages, and a past demographic expansion in at least one of them, generated high local clone diversities in A. rufus. These high clone diversities in asexual populations may mimic certain benefits conferred by sex via genetic diversity and could help explain the extreme success of asexual A. rufus populations.

Co-evolutionary dynamics between public good producers and cheats in the bacterium Pseudomonas aeruginosa

The production of beneficial public goods is common in the microbial world, and so is cheating--the exploitation of public goods by nonproducing mutants. Here, we examine co-evolutionary dynamics between cooperators and cheats and ask whether cooperators can evolve strategies to reduce the burden of exploitation, and whether cheats in turn can improve their exploitation abilities. We evolved cooperators of the bacterium Pseudomonas aeruginosa, producing the shareable iron-scavenging siderophore pyoverdine, together with cheats, defective in pyoverdine production but proficient in uptake. We found that cooperators managed to co-exist with cheats in 56% of all replicates over approximately 150 generations of experimental evolution. Growth and competition assays revealed that co-existence was fostered by a combination of general adaptions to the media and specific adaptions to the co-evolving opponent. Phenotypic screening and whole-genome resequencing of evolved clones confirmed this pattern, and suggest that cooperators became less exploitable by cheats because they significantly reduced their pyoverdine investment. Cheats, meanwhile, improved exploitation efficiency through mutations blocking the costly pyoverdine-signalling pathway. Moreover, cooperators and cheats evolved reduced motility, a pattern that likely represents adaptation to laboratory conditions, but at the same time also affects social interactions by reducing strain mixing and pyoverdine sharing. Overall, we observed parallel evolution, where co-existence of cooperators and cheats was enabled by a combination of adaptations to the abiotic and social environment and their interactions.

Switching between apparently redundant iron-uptake mechanisms benefits bacteria in changeable environments

Bacteria often possess multiple siderophore-based iron uptake systems for scavenging this vital resource from their environment. However, some siderophores seem redundant, because they have limited iron-binding efficiency and are seldom expressed under iron limitation. Here, we investigate the conundrum of why selection does not eliminate this apparent redundancy. We focus on Pseudomonas aeruginosa, a bacterium that can produce two siderophores-the highly efficient but metabolically expensive pyoverdine, and the inefficient but metabolically cheap pyochelin. We found that the bacteria possess molecular mechanisms to phenotypically switch from mainly producing pyoverdine under severe iron limitation to mainly producing pyochelin when iron is only moderately limited. We further show that strains exclusively producing pyochelin grew significantly better than strains exclusively producing pyoverdine under moderate iron limitation, whereas the inverse was seen under severe iron limitation. This suggests that pyochelin is not redundant, but that switching between siderophore strategies might be beneficial to trade off efficiencies versus costs of siderophores. Indeed, simulations parameterized from our data confirmed that strains retaining the capacity to switch between siderophores significantly outcompeted strains defective for one or the other siderophore under fluctuating iron availabilities. Finally, we discuss how siderophore switching can be viewed as a form of collective decision-making, whereby a coordinated shift in behaviour at the group level emerges as a result of positive and negative feedback loops operating among individuals at the local scale.

Cost of cooperation rules selection for cheats in bacterial metapopulations

Bacteria secrete a large variety of beneficial metabolites into the environment, which can be shared as public goods among producing bacteria, but also be exploited by nonproducing cheats. Here, we focus on cooperative production of iron-chelating molecules (siderophores) in the bacterium Pseudomonas aeruginosa to study how relevant ecological factors influence selection for cheating. We designed patch-structured metapopulations that allowed us introducing among-patch ecological variation. We found that cheating readily evolved in uniform iron-limited environments. This finding is explained by severe iron limitation demanding high siderophore-production efforts, which results in high metabolic costs accruing to cooperators, and thereby facilitates the spread of cheats. In contrast, we observed a significant reduction or even negation of selection for cheating in metapopulations where we introduced patches with increased iron availability and/or opportunities to recycle siderophores. These findings are compatible with the view that cheats are less likely to invade in environments that allow bacteria to reduce siderophore-production efforts, as this lowers the overall metabolic costs accruing to cooperators. Because we increased iron availability and siderophore recycling opportunities moderately, and only in some patches, our findings demonstrate that already-small local variations in ecological conditions as occurring in nature can significantly affect selection for public-goods secretion in microbes. In addition, we found that most (84.6%) of the evolved cheats were partially deficient for siderophore production and not loss-of-function mutants. Genetic considerations indicate that mutations leading to partial deficiency occur more frequent than mutations leading to loss of function, but also suggest that partially deficient mutants might often be the more competitive cheats.

Ever-young sex chromosomes in European tree frogs

Non-recombining sex chromosomes are expected to undergo evolutionary decay, ending up genetically degenerated, as has happened in birds and mammals. Why are then sex chromosomes so often homomorphic in cold-blooded vertebrates? One possible explanation is a high rate of turnover events, replacing master sex-determining genes by new ones on other chromosomes. An alternative is that X-Y similarity is maintained by occasional recombination events, occurring in sex-reversed XY females. Based on mitochondrial and nuclear gene sequences, we estimated the divergence times between European tree frogs (Hyla arborea, H. intermedia, and H. molleri) to the upper Miocene, about 5.4–7.1 million years ago. Sibship analyses of microsatellite polymorphisms revealed that all three species have the same pair of sex chromosomes, with complete absence of X-Y recombination in males. Despite this, sequences of sex-linked loci show no divergence between the X and Y chromosomes. In the phylogeny, the X and Y alleles cluster according to species, not in groups of gametologs. We conclude that sex-chromosome homomorphy in these tree frogs does not result from a recent turnover but is maintained over evolutionary timescales by occasional X-Y recombination. Seemingly young sex chromosomes may thus carry old-established sex-determining genes, a result at odds with the view that sex chromosomes necessarily decay until they are replaced. This raises intriguing perspectives regarding the evolutionary dynamics of sexually antagonistic genes and the mechanisms that control X-Y recombination.

Non-recombining sex chromosomes, such as the Y chromosome, are expected to degenerate over evolutionary times because they accumulate deleterious mutations that cannot be corrected by recombination with a pristine copy. In most cold-blooded vertebrates, such as frogs, however, sex chromosomes are undifferentiated. Why is that so? On the one hand, the “high-turnover” hypothesis holds that these sex chromosomes are regularly replaced before they had time to decay. On the other hand, the “fountain-of-youth” hypothesis posits that they are regularly rejuvenated by X-Y recombination in sex-reversed XY females. Here, we show that three species of tree frogs that diverged more than 5.4 million years ago share the same pair of undifferentiated sex chromosomes. Although male recombination stopped before species divergence, X and Y alleles show no differentiation, and cluster by species, not gametologs. We conclude that their sex chromosome homomorphy is not due to a recent turnover but is maintained over long evolutionary times by occasional recombination. Such rare episodes of X-Y recombination are expected to have long-lasting consequences on the evolution of sex chromosomes and sex antagonistic genes.