Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

The parasitic nematode Strongyloides ratti exists predominantly as populations of long-lived asexual lineages

Nematodes are important parasites of people and animals, and in natural ecosystems they are a major ecological force. Strongyloides ratti is a common parasitic nematode of wild rats and we have investigated its population genetics using single-worm, whole-genome sequencing. We find that S. ratti populations in the UK consist of mixtures of mainly asexual lineages that are widely dispersed across a host population. These parasite lineages are likely very old and may have originated in Asia from where rats originated. Genes that underly the parasitic phase of the parasite's life cycle are hyperdiverse compared with the rest of the genome, and this may allow the parasites to maximise their fitness in a diverse host population. These patterns of parasitic nematode population genetics have not been found before and may also apply to Strongyloides spp. that infect people, which will affect how we should approach their control.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Hermaphroditic origins of anisogamy

Anisogamy-the size dimorphism of gametes-is the defining difference between the male and female sexual strategies. Game-theoretic thinking led to the first convincing explanation for the evolutionary origins of anisogamy in the 1970s. Since then, formal game-theoretic models have continued to refine our understanding of when and why anisogamy should evolve. Such models typically presume that the earliest anisogamous organisms had separate sexes. However, in most taxa, there is no empirical evidence to support this assumption. Here, we present a model of the coevolution of gamete size and sex allocation, which allows for anisogamy to emerge alongside either hermaphroditism or separate sexes. We show that hermaphroditic anisogamy can evolve directly from isogamous ancestors when the average size of spawning groups is small and fertilization is relatively efficient. Sex allocation under hermaphroditism becomes increasingly female-biased as group size decreases and the degree of anisogamy increases. When spawning groups are very small, our model also predicts the existence of complex isogamous organisms in which individuals allocate resources equally to two large gamete types. We discuss common, but potentially unwarranted, assumptions in the literature that could be relaxed in future models. This article is part of the theme issue 'Half a century of evolutionary games: a synthesis of theory, application and future directions'.

Conflicting effects of recombination on the evolvability and robustness in neutrally evolving populations

Understanding the benefits and costs of recombination under different scenarios of evolutionary adaptation remains an open problem for theoretical and experimental research. In this study, we focus on finite populations evolving on neutral networks comprising viable and unfit genotypes. We provide a comprehensive overview of the effects of recombination by jointly considering different measures of evolvability and mutational robustness over a broad parameter range, such that many evolutionary regimes are covered. We find that several of these measures vary non-monotonically with the rates of mutation and recombination. Moreover, the presence of unfit genotypes that introduce inhomogeneities in the network of viable states qualitatively alters the effects of recombination. We conclude that conflicting trends induced by recombination can be explained by an emerging trade-off between evolvability on the one hand, and mutational robustness on the other. Finally, we discuss how different implementations of the recombination scheme in theoretical models can affect the observed dependence on recombination rate through a coupling between recombination and genetic drift.

Genetic variation and clonal diversity in floating aquatic plants: Comparative genomic analysis of water hyacinth species in their native range

Many eukaryotic organisms reproduce by sexual and asexual reproduction. Genetic diversity in populations can be strongly dependent on the relative importance of these two reproductive modes. Here, we compare the amounts and patterns of genetic diversity in related water hyacinths that differ in their propensity for clonal propagation - highly clonal Eichhornia crassipes and moderately clonal E. azurea (Pontederiaceae). Our comparisons involved genotype-by-sequencing (GBS) of 137 E. crassipes ramets from 60 locations (193,495 nucleotide sites) and 118 E. azurea ramets from 53 locations (198,343 nucleotide sites) among six hydrological basins in central South America, the native range of both species. We predicted that because of more prolific clonal propagation, E. crassipes would exhibit lower clonal diversity than E. azurea. This prediction was supported by all measures of clonal diversity that we examined. Eichhornia crassipes also had a larger excess of heterozygotes at variant sites, another signature of clonality. However, genome-wide heterozygosity was not significantly different between the species. Eichhornia crassipes had weaker spatial genetic structure and lower levels of differentiation among hydrological basins than E. azurea, probably because of higher clonality and more extensive dispersal of its free-floating life form. Our findings for E. crassipes contrast with earlier studies from the invasive range which have reported very low levels of clonal diversity and extensive geographic areas of genetic uniformity.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Genome sequencing and haplotype assembly of two cyprinid teleosts, a sexual tetraploid and an unisexual hexaploid, reveal insights into the evolutionary mechanisms underpinning the reproductive success of unisexual polyploid vertebrates.

Genetic Variation in Reproductive Investment Across an Ephemerality Gradient in Daphnia pulex

Species across the tree of life can switch between asexual and sexual reproduction. In facultatively sexual species, the ability to switch between reproductive modes is often environmentally dependent and subject to local adaptation. However, the ecological and evolutionary factors that influence the maintenance and turnover of polymorphism associated with facultative sex remain unclear. We studied the ecological and evolutionary dynamics of reproductive investment in the facultatively sexual model species, Daphnia pulex. We found that patterns of clonal diversity, but not genetic diversity varied among ponds consistent with the predicted relationship between ephemerality and clonal structure. Reconstruction of a multi-year pedigree demonstrated the coexistence of clones that differ in their investment into male production. Mapping of quantitative variation in male production using lab-generated and field-collected individuals identified multiple putative quantitative trait loci (QTL) underlying this trait, and we identified a plausible candidate gene. The evolutionary history of these QTL suggests that they are relatively young, and male limitation in this system is a rapidly evolving trait. Our work highlights the dynamic nature of the genetic structure and composition of facultative sex across space and time and suggests that quantitative genetic variation in reproductive strategy can undergo rapid evolutionary turnover.

Invasions of an obligate asexual daphnid species support the nearly neutral theory

To verify the "nearly neutral theory (NNT)," the ratio of nonsynonymous to synonymous substitutions (dN/dS) was compared among populations of different species. To determine the validity of NNT, however, populations that are genetically isolated from each other but share the same selection agents and differ in size should be compared. Genetically different lineages of obligate asexual Daphnia pulex invading Japan from North America are an ideal example as they satisfy these prerequisites. Therefore, we analyzed the whole-genome sequences of 18 genotypes, including those of the two independently invaded D. pulex lineages (JPN1 and JPN2) and compared the dN/dS ratio between the lineages. The base substitution rate of each genotype demonstrated that the JPN1 lineage having a larger distribution range diverged earlier and thus was older than the JPN2 lineage. Comparisons of the genotypes within lineages revealed that changes in dN/dS occurred after the divergence and were larger in the younger lineage, JPN2. These results imply that the JPN1 lineage has been more effectively subjected to purification selections, while slightly deteriorating mutations are less purged in JPN2 with smaller population size. Altogether, the lineage-specific difference in the dN/dS ratio for the obligate asexual D. pulex was well explained by the NNT.

Evidence from automixis with inverted meiosis for the maintenance of sex by loss of complementation

The adaptive value of sexual reproduction is still debated. A short-term disadvantage of asexual reproduction is loss of heterozygosity, which leads to the unmasking of recessive deleterious mutations. The cost of this loss of complementation is predicted to be higher than the twofold cost of meiosis for most types of asexual reproduction. Automixis with terminal fusion of sister nuclei is especially vulnerable to the effect of loss of complementation. It is found, however, in some taxa including oribatid mites, the most prominent group of ancient asexuals. Here, I show that automixis with terminal fusion is stable if it is associated with inverted meiosis and that this appears to be the case in nature, notably in oribatid mites. The existence of automixis with terminal fusion, and its co-occurrence with inverted meiosis, therefore, is consistent with the hypothesis that loss of complementation is important in the evolution of sexual reproduction.

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene's functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.

Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

Environmental fluctuations dampen the effects of clonal reproduction on evolutionary rescue

Evolutionary rescue occurs when genetic change allows a population to persist in response to an environmental change that would otherwise have led to extinction. Most studies of evolutionary rescue assume that species have either fully clonal or fully sexual reproduction; however, many species have partially clonal reproductive strategies in which they reproduce both clonally and sexually. Furthermore, the few evolutionary rescue studies that have evaluated partially clonal reproduction did not consider fluctuations in the environment, which are nearly ubiquitous in nature. Here, we use individual-based simulations to investigate how environmental fluctuations (either uncorrelated or positively autocorrelated) influence the effect of clonality on evolutionary rescue. We show that, for moderate magnitudes of environmental fluctuations, as was found in the absence of fluctuations, increasing the degree of clonality increases the probability of population persistence in response to an abrupt environmental change, but decreases persistence in response to a continuous, directional environmental change. However, with large magnitudes of fluctuations, both the benefits of clonality following a step change and the detrimental effects of clonality following a continuous, directional change are generally reduced; in fact, in the latter scenario, increasing clonality can even become beneficial if environmental fluctuations are autocorrelated. We also show that increased generational overlap dampens the effects of environmental fluctuations. Overall, we demonstrate that understanding the evolutionary rescue of partially clonal organisms requires not only knowledge of the species life history and the type of environmental change, but also an understanding of the magnitude and autocorrelation of environmental fluctuations.

How and When Does Outcrossing Occur in the Predominantly Selfing Species Medicago truncatula?

Empirical studies on natural populations of Medicago truncatula revealed selfing rates higher than 80%, but never up to 100%. Similarly, several studies of predominantly selfing species show variability in the level of residual outcrossing between populations and also between temporal samples of the same population. However, these studies measure global selfing rates at the scale of the population and we do not know whether there is intra-population variation and how outcrossing events are distributed, between genotypes, plants, flowers, or seeds. Theoretical studies predict the maintenance of residual outcrossing in highly selfing species due to environmental (e.g., pollen biology) and/or genetic determinants and decompositions of the variation in outcrossing rate using experimental data can be very informative to test these hypotheses. Here, we focus on one natural population of M. truncatula in order to describe precisely its mating system. In particular, we investigated the determinants of the selfing rate by testing for seasonal variations (environmental determinism) and variations between genotypes (genetic determinism). We measured selfing rates in maternal progenies from plants collected widely across a natural population. For each plant, we collected pods from flowers produced at the beginning and at the end of the flowering season to test for a seasonal variation in the outcrossing rate. For each collected offspring, we also estimated the likelihood that it was issued from a self-fertilization event and assessed the genetic component of variation of this mating system measure. We found a significant, albeit small, increase in outcrossing rate in progenies collected at the end [t _m = 0.137 (SD = 0.025)] compared to those collected at the beginning [t _m = 0.083 (0.016)] of the flowering season. A significant between genotypes variation in selfing rate was also detected, resulting in a heritability of 9% for the rate of residual outcrossing. Altogether, our work shows that despite a predominantly selfing reproductive mode, M. truncatula displays variation in residual outcrossing rate, and that this trait is likely under a complex determinism combining environmental and genetic factors. We discuss the evolutionary implications of our results for the population.

Approximating the Coalescent Under Facultative Sex

Genome studies of facultative sexual species, which can either reproduce sexually or asexually, are providing insight into the evolutionary consequences of mixed reproductive modes. It is currently unclear to what extent the evolutionary history of facultative sexuals' genomes can be approximated by the standard coalescent, and if a coalescent effective population size Ne exists. Here, I determine if and when these approximations can be made. When sex is frequent (occurring at a frequency much greater than 1/N per reproduction per generation, for N the actual population size), the underlying genealogy can be approximated by the standard coalescent, with a coalescent Ne≈N. When sex is very rare (at frequency much lower than 1/N), approximations for the pairwise coalescent time can be obtained, which is strongly influenced by the frequencies of sex and mitotic gene conversion, rather than N. However, these terms do not translate into a coalescent Ne. These results are used to discuss the best sampling strategies for investigating the evolutionary history of facultative sexual species.

Sex and Asex: A Clonal Lexicon

Organisms across the tree of life have complex life cycles that include both sexual and asexual reproduction or that are obligately asexual. These organisms include ecologically dominant species that structure many terrestrial and marine ecosystems, as well as many pathogens, pests, and invasive species. We must consider both the evolution and maintenance of these various reproductive modes and how these modes shape the genetic diversity, adaptive evolution, and ability to persist in the species that exhibit them. Thus, having a common framework is a key aspect of understanding the biodiversity that shapes our planet. In the 2019 AGA President's Symposium, Sex and Asex: The genetics of complex life cycles, researchers investigating a wide range of taxonomic models and using a variety of modes of investigation coalesced around a common theme-understanding not only how such complex life cycles may evolve, but how they are shaped by the evolutionary and ecological forces around them. In this introduction to the Special Issue from the symposium, we give an overview of some of the key ideas and areas of investigation (a common clonal lexicon, we might say) and introduce the breadth of work submitted by symposium participants.

Does Gene Expression Noise Play a Functional Role in Plants?

Gene expression in individual cells can be surprisingly noisy. In unicellular organisms this noise can be functional; for example, by allowing a subfraction of the population to prepare for environmental stress. The role of gene expression noise in multicellular organisms has, however, remained unclear. In this review, we discuss how new techniques are revealing an unexpected level of variability in gene expression between and within genetically identical plants. We describe recent progress as well as speculate on the function of transcriptional noise as a mechanism for generating functional phenotypic diversity in plants.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Genomic and proteomic biases inform metabolic engineering strategies for anaerobic fungi

Anaerobic fungi (Neocallimastigomycota) are emerging non-model hosts for biotechnology due to their wealth of biomass-degrading enzymes, yet tools to engineer these fungi have not yet been established. Here, we show that the anaerobic gut fungi have the most GC depleted genomes among 443 sequenced organisms in the fungal kingdom, which has ramifications for heterologous expression of genes as well as for emerging CRISPR-based genome engineering approaches. Comparative genomic analyses suggest that anaerobic fungi may contain cellular machinery to aid in sexual reproduction, yet a complete mating pathway was not identified. Predicted proteomes of the anaerobic fungi also contain an unusually large fraction of proteins with homopolymeric amino acid runs consisting of five or more identical consecutive amino acids. In particular, threonine runs are especially enriched in anaerobic fungal carbohydrate active enzymes (CAZymes) and this, together with a high abundance of predicted N-glycosylation motifs, suggests that gut fungal CAZymes are heavily glycosylated, which may impact heterologous production of these biotechnologically useful enzymes. Finally, we present a codon optimization strategy to aid in the development of genetic engineering tools tailored to these early-branching anaerobic fungi.

Population Genomics of Transitions to Selfing in Brassicaceae Model Systems

Many plants harbor complex mechanisms that promote outcrossing and efficient pollen transfer. These include floral adaptations as well as genetic mechanisms, such as molecular self-incompatibility (SI) systems. The maintenance of such systems over long evolutionary timescales suggests that outcrossing is favorable over a broad range of conditions. Conversely, SI has repeatedly been lost, often in association with transitions to self-fertilization (selfing). This transition is favored when the short-term advantages of selfing outweigh the costs, primarily inbreeding depression. The transition to selfing is expected to have major effects on population genetic variation and adaptive potential, as well as on genome evolution. In the Brassicaceae, many studies on the population genetic, gene regulatory, and genomic effects of selfing have centered on the model plant Arabidopsis thaliana and the crucifer genus Capsella. The accumulation of population genomics datasets have allowed detailed investigation of where, when and how the transition to selfing occurred. Future studies will take advantage of the development of population genetics theory on the impact of selfing, especially regarding positive selection. Furthermore, investigation of systems including recent transitions to selfing, mixed mating populations and/or multiple independent replicates of the same transition will facilitate dissecting the effects of mating system variation from processes driven by demography.

When Synchrony Makes the Best of Both Worlds Even Better: How Well Do We Really Understand Facultative Sex?

Biological diversity abounds in potential study topics. Studies of model systems have their advantages, but reliance on a few well-understood cases may create false impressions of what biological phenomena are the norm. Here I focus on facultative sex, which is often hailed as offering the best of both worlds, in that rare sex offers benefits almost equal to obligate sex and avoids paying most of the demographic costs. How well do we understand when and why this form of sexual reproduction is expected to prevail? I show several gaps in the theoretical literature and, by contrasting asynchronous with synchronous sex, highlight the need to link sex theories to the theoretical underpinnings of bet hedging, on the one hand, and to mate limitation considerations, on the other. Condition-dependent sex and links between sex with dispersal or dormancy appear understudied. While simplifications are justifiable as a simple assumption structure enhances analytical tractability, much remains to be done to incorporate key features of real sex to the main theoretical edifice.

Population genomics of the facultatively asexual duckweed Spirodela polyrhiza

Clonal propagation allows some plant species to achieve massive population sizes quickly but also reduces the evolutionary independence of different sites in the genome. We examine genome-wide genetic diversity in Spirodela polyrhiza, a duckweed that reproduces primarily asexually. We find that this geographically widespread and numerically abundant species has very low levels of genetic diversity. Diversity at nonsynonymous sites relative to synonymous sites is high, suggesting that purifying selection is weak. A potential explanation for this observation is that a very low frequency of sex renders selection ineffective. However, there is a pronounced decay in linkage disequilibrium over 40 kb, suggesting that though sex may be rare at the individual level it is not too infrequent at the population level. In addition, neutral diversity is affected by the physical proximity of selected sites, which would be unexpected if sex was exceedingly rare at the population level. The amount of genetic mixing as assessed by the decay in linkage disequilibrium is not dissimilar from selfing species such as Arabidopsis thaliana, yet selection appears to be much less effective in duckweed. We discuss alternative explanations for the signature of weak purifying selection.

Hybridization and a loss of sex shape genome-wide diversity and the origin of species in the evening primroses (Oenothera, Onagraceae)

Hybridization is thought to promote speciation in at least two ways - by fixation of heterozygosity from diploid progenitors in allopolyploids, and by generation of transgressive phenotypes and shifting fitness optima during homoploid hybrid speciation. While recent studies support a hybrid origin for a growing number of species, the extent to which hybrid origins shape patterns of diversity in asexual species remains underexplored. Here we employed transcriptome sequencing and population genomic analysis to describe patterns of genomic variation in the 13 species belonging to Oenothera subsection Oenothera. Eight of these species are functionally asexual and arose by hybrid speciation from parents spanning a range of phylogenetic divergence. We showed that genomic divergence between parents has been retained as heterozygosity in functionally asexual species, and that genome-wide levels of heterozygosity in these asexuals largely reflects the divergence of parental haplotypes coupled with a breakdown in recombination and segregation across the genome. These results show that divergence between parental species and loss of sex in hybrids shape patterns of whole-genome diversity and the origin of asexual species.

Outstanding intraindividual genetic diversity in fissiparous planarians (Dugesia, Platyhelminthes) with facultative sex

BACKGROUND

Predicted genetic consequences of asexuality include high intraindividual genetic diversity (i.e., the Meselson effect) and accumulation of deleterious mutations (i.e., Muller's Ratchet), among others. These consequences have been largely studied in parthenogenetic organisms, but studies on fissiparous species are scarce. Differing from parthenogens, fissiparous organisms inherit part of the soma of the progenitor, including somatic mutations. Thus, in the long term, fissiparous reproduction may also result in genetic mosaicism, besides the presence of the Meselson effect and Muller's Ratchet. Dugesiidae planarians show outstanding regeneration capabilities, allowing them to naturally reproduce by fission, either strictly or combined with sex (facultative). Therefore, they are an ideal model to analyze the genetic footprint of fissiparous reproduction, both when it is alternated with sex and when it is the only mode of reproduction.

RESULTS

In the present study, we generate and analyze intraindividual cloned data of a nuclear and a mitochondrial gene of sexual, fissiparous and facultative wild populations of the species Dugesia subtentaculata. We find that most individuals, independently of their reproductive strategy, are mosaics. However, the intraindividual haplotype and nucleotide diversity of fissiparous and facultative individuals is significantly higher than in sexual individuals, with no signs of Muller's Ratchet. Finally, we also find that this high intraindividual genetic diversity of fissiparous and facultative individuals is composed by different combinations of ancestral and derived haplotypes of the species.

CONCLUSIONS

The intraindividual analyses of genetic diversity point out that fissiparous reproduction leaves a very special genetic footprint in individuals, characterized by mosaicism combined with the Meselson effect (named in the present study as the mosaic Meselson effect). Interestingly, the different intraindividual combinations of ancestral and derivate genetic diversity indicate that haplotypes generated during periods of fissiparous reproduction can be also transmitted to the progeny through sexual events, resulting in offspring showing a wide range of genetic diversity and putatively allowing purifying selection to act at both intraindividual and individual level. Further investigations, using Dugesia planarians as model organisms, would be of great value to delve into this new model of genetic evolution by the combination of fission and sex.

How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species

Most theoretical works predict that selfing should reduce the level of additive genetic variance available for quantitative traits within natural populations. Despite a growing number of quantitative genetic studies undertaken during the last two decades, this prediction is still not well supported empirically. To resolve this issue and confirm or reject theoretical predictions, we reviewed quantitative trait heritability estimates from natural plant populations with different rates of self-fertilization and carried out a meta-analysis. In accordance with models of polygenic traits under stabilizing selection, we found that the fraction of additive genetic variance is negatively correlated with the selfing rate. Although the mating system explains a moderate fraction of the variance, the mean reduction of narrow-sense heritability values between strictly allogamous and predominantly selfing populations is strong, around 60%. Because some nonadditive components of genetic variance become selectable under inbreeding, we determine whether self-fertilization affects the relative contribution of these components to genetic variance by comparing narrow-sense heritability estimates from outcrossing populations with broad-sense heritability estimated in autogamous populations. Results suggest that these nonadditive components of variance may restore some genetic variance in predominantly selfing populations; it remains, however, uncertain how these nonadditive components will contribute to adaptation.

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.

How clonal are clones? A quest for loss of heterozygosity during asexual reproduction in Daphnia magna

Due to the lack of recombination, asexual organisms are predicted to accumulate mutations and show high levels of within‐individual allelic divergence (heterozygosity); however, empirical evidence for this prediction is largely missing. Instead, evidence of genome homogenization during asexual reproduction is accumulating. Ameiotic crossover recombination is a mechanism that could lead to long genomic stretches of loss of heterozygosity (LOH) and unmasking of mutations that have little or no effect in heterozygous state. Therefore, LOH might be an important force for inducing variation among asexual offspring and may contribute to the limited longevity of asexual lineages. To investigate the genetic consequences of asexuality, here we used high‐throughput sequencing of Daphnia magna for assessing the rate of LOH over a single generation of asexual reproduction. Comparing parthenogenetic daughters with their mothers at several thousand genetic markers generated by restriction site‐associated DNA (RAD) sequencing resulted in high LOH rate estimation that largely overlapped with our estimates for the error rate. To distinguish these two, we Sanger re‐sequenced the top 17 candidate RAD‐loci for LOH, and all of them proved to be false positives. Hence, even though we cannot exclude the possibility that short stretches of LOH occur in genomic regions not covered by our markers, we conclude that LOH does not occur frequently during asexual reproduction in D. magna and ameiotic crossovers are very rare or absent. This finding suggests that clonal lineages of D. magna will remain genetically homogeneous at least over time periods typically relevant for experimental work.

No signal of deleterious mutation accumulation in conserved gene sequences of extant asexual hexapods

Loss of sex and recombination is generally assumed to impede the effectiveness of purifying selection and to result in the accumulation of slightly deleterious mutations. Empirical evidence for this has come from several studies investigating mutational load in a small number of individual genes. However, recent whole transcriptome based studies have yielded inconsistent results, hence questioning the validity of the assumption of mutational meltdown in asexual populations. Here, we study the effectiveness of purifying selection in eight asexual hexapod lineages and their sexual relatives, as present in the 1 K Insect Transcriptome Evolution (1KITE) project, covering eight hexapod groups. We analyse the accumulation of slightly deleterious nonsynonymous and synonymous point mutations in 99 single copy orthologue protein-coding loci shared among the investigated taxa. While accumulation rates of nonsynonymous mutations differed between genes and hexapod groups, we found no effect of reproductive mode on the effectiveness of purifying selection acting at nonsynonymous and synonymous sites. Although the setup of this study does not fully rule out nondetection of subtle effects, our data does not support the established consensus of asexual lineages undergoing ‘mutational meltdown’.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

Coalescence and Linkage Disequilibrium in Facultatively Sexual Diploids

Under neutrality, linkage disequilibrium results from physically linked sites having nonindependent coalescent histories. In obligately sexual organisms, meiotic recombination is the dominant force separating linked variants from one another, and thus in determining the decay of linkage disequilibrium with physical distance. In facultatively sexual diploid organisms that principally reproduce clonally, mechanisms of mitotic exchange are expected to become relatively more important in shaping linkage disequilibrium. Here we outline mathematical and computational models of a facultative-sex coalescent process that includes meiotic and mitotic recombination, via both crossovers and gene conversion, to determine how linkage disequilibrium is affected with facultative sex. We demonstrate that the degree to which linkage disequilibrium is broken down by meiotic recombination simply scales with the probability of sex if it is sufficiently high (much greater than [Formula: see text] for population size N). However, with very rare sex (occurring with frequency on the order of [Formula: see text]), mitotic gene conversion plays a particularly important and complicated role because it both breaks down associations between sites and removes within-individual diversity. Strong population structure under rare sex leads to lower average linkage disequilibrium values than in panmictic populations, due to the influence of low-frequency polymorphisms created by allelic sequence divergence acting in individual subpopulations. These analyses provide information on how to interpret observed linkage disequilibrium patterns in facultative sexuals and to determine what genomic forces are likely to shape them.

Low genetic diversity contrasts with high phenotypic variability in heptaploid Spartina densiflora populations invading the Pacific coast of North America

Species can respond to environmental pressures through genetic and epigenetic changes and through phenotypic plasticity, but few studies have evaluated the relationships between genetic differentiation and phenotypic plasticity of plant species along changing environmental conditions throughout wide latitudinal ranges. We studied inter‐ and intrapopulation genetic diversity (using simple sequence repeats and chloroplast DNA sequencing) and inter‐ and intrapopulation phenotypic variability of 33 plant traits (using field and common‐garden measurements) for five populations of the invasive cordgrass Spartina densiflora Brongn. along the Pacific coast of North America from San Francisco Bay to Vancouver Island. Studied populations showed very low genetic diversity, high levels of phenotypic variability when growing in contrasted environments and high intrapopulation phenotypic variability for many plant traits. This intrapopulation phenotypic variability was especially high, irrespective of environmental conditions, for those traits showing also high phenotypic plasticity. Within‐population variation represented 84% of the total genetic variation coinciding with certain individual plants keeping consistent responses for three plant traits (chlorophyll b and carotenoid contents, and dead shoot biomass) in the field and in common‐garden conditions. These populations have most likely undergone genetic bottleneck since their introduction from South America; multiple introductions are unknown but possible as the population from Vancouver Island was the most recent and one of the most genetically diverse. S. densiflora appears as a species that would not be very affected itself by climate change and sea‐level rise as it can disperse, establish, and acclimate to contrasted environments along wide latitudinal ranges.

Demography and mating system shape the genome-wide impact of purifying selection in Arabis alpina

Plant mating systems have profound effects on levels and structuring of genetic variation and can affect the impact of natural selection. Although theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population-genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer Arabis alpina We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian A. alpina show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles, whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, mating system and demography shape the impact of purifying selection on genomic variation in A. alpina These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.

Sexual conflict and the evolution of asexuality at low population densities

Theories for the evolution of sex rarely include facultatively sexual reproduction. Sexual harassment by males is an underappreciated factor: it should at first sight increase the relative advantage of asexual reproduction by increasing the cost of sex. However, if the same females can perform either sexual or asexual life cycles, then females trying to reproduce asexually may not escape harassment. If resisting male harassment is costly, it might be beneficial for a female to accept a mating and undertake a sexual life cycle rather than 'insist' on an asexual one. We investigate the effects of sexual harassment on the maintenance of sex under different population densities. Our model shows that resisting matings pays off at low population densities, which leads to the complete extinction of males, and thus to the evolution of completely asexual populations. Facultative sex persists in a narrow range of slightly higher densities. At high densities, selection favours giving up resisting male mating attempts and thus sexual reproduction takes over. These interactions between the outcomes of sexual conflict and population density suggest an explanation for the rarity of facultative sex and also patterns of geographical parthenogenesis, where marginal environments with potentially low densities are associated with asexuality.

The frequency of sex in fungi

Fungi are a diverse group of organisms with a huge variation in reproductive strategy. While almost all species can reproduce sexually, many reproduce asexually most of the time. When sexual reproduction does occur, large variation exists in the amount of in- and out-breeding. While budding yeast is expected to outcross only once every 10 000 generations, other fungi are obligate outcrossers with well-mixed panmictic populations. In this review, we give an overview of the costs and benefits of sexual and asexual reproduction in fungi, and the mechanisms that evolved in fungi to reduce the costs of either mode. The proximate molecular mechanisms potentiating outcrossing and meiosis appear to be present in nearly all fungi, making them of little use for predicting outcrossing rates, but also suggesting the absence of true ancient asexual lineages. We review how population genetic methods can be used to estimate the frequency of sex in fungi and provide empirical data that support a mixed mode of reproduction in many species with rare to frequent sex in between rounds of mitotic reproduction. Finally, we highlight how these estimates might be affected by the fungus-specific mechanisms that evolved to reduce the costs of sexual and asexual reproduction.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

On the origin of asexual species by means of hybridization and drift

Several species of asexuals appear to have existed for millions of years. This is despite the prevalent view that natural selection is weakened without gene exchange, which should cause these organisms to rapidly go extinct. In theory, one can identify evolutionary long-lived asexuals from their allelic sequence divergence, also known as the 'Meselson effect', which leads to elevated within-individual diversity. Yet evidence that this phenomenon exists is mixed. Furthermore, several confounding factors can lead to similar outcomes, including the formation of asexual species by hybridization. Disentangling these factors has proved to be tricky, but Ament-Velásquez et al. (2016) have provided an elegant solution in this issue of Molecular Ecology. They studied transcriptomes and mitochondrial DNA from the Lineus genus of nemertean worms, which contains both sexual and asexual types, to first show that the asexual L. pseudolactus is a hybrid between a sexual and an asexual species. After isolating out diversity arising from this hybridization, they find subsequent evidence for the Meselson effect. This study sets a new standard for differentiating between the complex causes and consequences of asexuality.

Effective purifying selection in ancient asexual oribatid mites

Sex is beneficial in the long term because it can prevent mutational meltdown through increased effectiveness of selection. This idea is supported by empirical evidence of deleterious mutation accumulation in species with a recent transition to asexuality. Here, we study the effectiveness of purifying selection in oribatid mites which have lost sex millions of years ago and diversified into different families and species while reproducing asexually. We compare the accumulation of deleterious nonsynonymous and synonymous mutations between three asexual and three sexual lineages using transcriptome data. Contrasting studies of young asexual lineages, we find evidence for strong purifying selection that is more effective in asexual as compared to sexual oribatid mite lineages. Our results suggest that large populations likely sustain effective purifying selection and facilitate the escape of mutational meltdown in the absence of sex. Thus, sex per se is not a prerequisite for the long-term persistence of animal lineages.

Asexual reproduction is thought to be an evolutionary dead end in eukaryotes because deleterious mutations will not be purged effectively. Here, Brandt and colleagues show that anciently asexual oribatid mites in fact have reduced accumulation of deleterious mutations compared to their sexual relatives.

Chemosensory and hyperoxia circuits in C. elegans males influence sperm navigational capacity

The sperm’s crucial function is to locate and fuse with a mature oocyte. Under laboratory conditions, Caenorhabditis elegans sperm are very efficient at navigating the hermaphrodite reproductive tract and locating oocytes. Here, we identify chemosensory and oxygen-sensing circuits that affect the sperm’s navigational capacity. Multiple Serpentine Receptor B (SRB) chemosensory receptors regulate Gα pathways in gustatory sensory neurons that extend cilia through the male nose. SRB signaling is necessary and sufficient in these sensory neurons to influence sperm motility parameters. The neuropeptide Y pathway acts together with SRB-13 to antagonize negative effects of the GCY-35 hyperoxia sensor on spermatogenesis. SRB chemoreceptors are not essential for sperm navigation under low oxygen conditions that C. elegans prefers. In ambient oxygen environments, SRB-13 signaling impacts gene expression during spermatogenesis and the sperm’s mitochondria, thereby increasing migration velocity and inhibiting reversals within the hermaphrodite uterus. The SRB-13 transcriptome is highly enriched in genes implicated in pathogen defense, many of which are expressed in diverse tissues. We show that the critical time period for SRB-13 signaling is prior to spermatocyte differentiation. Our results support the model that young C. elegans males sense external environment and oxygen tension, triggering long-lasting downstream signaling events with effects on the sperm’s mitochondria and navigational capacity. Environmental exposures early in male life may alter sperm function and fertility.

Habitat loss, disease, climate change, and pollution are thought to negatively affect animal fertility. Sperm are a potential target, but the molecular mechanisms are not understood. The nematode C. elegans is a powerful genetic model to investigate the relationship between environment and male fertility. The hermaphrodite’s transparent epidermis permits the direct visualization of migrating male sperm and fertilization. In this study, we identified multiple serpentine receptor B (SRB) chemosensory receptors that are expressed in amphid sensory neurons, which extend cilia through the male nose. These SRB chemoreceptors are necessary to produce sperm that are efficient at navigating the hermaphrodite reproductive tract to the fertilization site. We show that SRB-13 signaling counteracts the negative effect of GCY-35 O₂ sensor activity, thereby maintaining sperm mitochondrial function and navigational capacity in hyperoxic conditions. Of particular interest, SRB-13 acts in early larval stage males prior to testis maturation. We propose that young males respond to specific stressful environments by altering SRB neural circuits, which in turn impact sperm mitochondrial function and motility. This chemosensory mechanism may be part of a systemic response in C. elegans males to external environment and oxygen levels.

Aging asexual lineages and the evolutionary maintenance of sex

Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within-population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages "senesce" in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.

The Evolutionary Interplay between Adaptation and Self-Fertilization

Genome-wide surveys of nucleotide polymorphisms, obtained from next-generation sequencing, have uncovered numerous examples of adaptation in self-fertilizing organisms, especially regarding changes to climate, geography, and reproductive systems. Yet existing models for inferring attributes of adaptive mutations often assume idealized outcrossing populations, which risks mischaracterizing properties of these variants. Recent theoretical work is emphasizing how various aspects of self-fertilization affects adaptation, yet empirical data on these properties are lacking. We review theoretical and empirical studies demonstrating how self-fertilization alters the process of adaptation, illustrated using examples from current sequencing projects. We propose ideas for how future research can more accurately quantify aspects of adaptation in self-fertilizers, including incorporating the effects of standing variation, demographic history, and polygenic adaptation.

Analysis of large-scale next-generation sequencing datasets are finding more examples of adaptive evolution at the genomic level.

Advances in theoretical work has demonstrated how self-fertilisation affects different aspects of adaptation in these organisms, compared to outcrossers.

Current software and statistical methods do not take different mating systems into account, which risks mischaracterising the presence or strength of adaptive mutations from genome scans.

Development of new mathematical and statistical methods that explicitly consider self-fertilization and associated demographic effects will enable researchers to more accurately quantify adaptation in these organisms.

Determinants of genetic diversity

Genetic polymorphism varies among species and within genomes, and has important implications for the evolution and conservation of species. The determinants of this variation have been poorly understood, but population genomic data from a wide range of organisms now make it possible to delineate the underlying evolutionary processes, notably how variation in the effective population size (Ne) governs genetic diversity. Comparative population genomics is on its way to providing a solution to 'Lewontin's paradox' - the discrepancy between the many orders of magnitude of variation in population size and the much narrower distribution of diversity levels. It seems that linked selection plays an important part both in the overall genetic diversity of a species and in the variation in diversity within the genome. Genetic diversity also seems to be predictable from the life history of a species.