Consequences of Asexuality in Natural Populations: Insights from Stick Insects

Parthenogenetic Stick Insects Exhibit Signatures of Preservation in the Molecular Architecture of Male Reproduction

After the loss of a trait, theory predicts that the molecular machinery underlying its phenotypic expression should decay. Yet, empirical evidence is contrasting. Here, we test the hypotheses that (i) the molecular ground plan of a lost trait could persist due to pleiotropic effects on other traits and (ii) that gene co-expression network architecture could constrain individual gene expression. Our testing ground has been the Bacillus stick insect species complex, which contains close relatives that are either bisexual or parthenogenetic. After the identification of genes expressed in male reproductive tissues in a bisexual species, we investigated their gene co-expression network structure in two parthenogenetic species. We found that gene co-expression within the male gonads was partially preserved in parthenogens. Furthermore, parthenogens did not show relaxed selection on genes upregulated in male gonads in the bisexual species. As these genes were mostly expressed in female gonads, this preservation could be driven by pleiotropic interactions and an ongoing role in female reproduction. Connectivity within the network also played a key role, with highly connected-and more pleiotropic-genes within male gonad also having a gonad-biased expression in parthenogens. Our findings provide novel insight into the mechanisms which could underlie the production of rare males in parthenogenetic lineages; more generally, they provide an example of the cryptic persistence of a lost trait molecular architecture, driven by gene pleiotropy on other traits and within their co-expression network.

Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

Genetic and Phenotypic Consequences of Local Transitions between Sexual and Parthenogenetic Reproduction in the Wild

AbstractTransitions from sexual to asexual reproduction have occurred in numerous lineages, but it remains unclear why asexual populations rarely persist. In facultatively parthenogenetic animals, all-female populations can arise when males are absent or become extinct, and such populations could help to understand the genetic and phenotypic changes that occur in the initial stages of transitions to asexuality. We investigated a naturally occurring spatial mosaic of mixed-sex and all-female populations of the facultatively parthenogenetic Australian phasmid Megacrania batesii. Analysis of single-nucleotide polymorphisms indicated multiple independent transitions between reproductive modes. All-female populations had much lower heterozygosity and allelic diversity than mixed-sex populations, but we found few consistent differences in fitness-related traits between population types. All-female populations exhibited more frequent and severe deformities in their (flight-incapable) wings but did not show higher rates of appendage loss. All-female populations also harbored more ectoparasites in swamp (but not beach) habitats. Reproductive mode explained little variation in female body size, fecundity, or egg hatch rate. Our results suggest that transitions to parthenogenetic reproduction can lead to dramatic genetic changes with little immediate effect on performance. All-female M. batesii populations appear to consist of high-fitness genotypes that might be able to thrive for many generations in relatively constant and benign environments but could be vulnerable to environmental challenges, such as increased parasite abundance.

AlexandrusPS: A User-Friendly Pipeline for the Automated Detection of Orthologous Gene Clusters and Subsequent Positive Selection Analysis

The detection of adaptive selection in a system approach considering all protein-coding genes allows for the identification of mechanisms and pathways that enabled adaptation to different environments. Currently, available programs for the estimation of positive selection signals can be divided into two groups. They are either easy to apply but can analyze only one gene family at a time, restricting system analysis; or they can handle larger cohorts of gene families, but require considerable prerequisite data such as orthology associations, codon alignments, phylogenetic trees, and proper configuration files. All these steps require extensive computational expertise, restricting this endeavor to specialists. Here, we introduce AlexandrusPS, a high-throughput pipeline that overcomes technical challenges when conducting transcriptome-wide positive selection analyses on large sets of nucleotide and protein sequences. The pipeline streamlines 1) the execution of an accurate orthology prediction as a precondition for positive selection analysis, 2) preparing and organizing configuration files for CodeML, 3) performing positive selection analysis using CodeML, and 4) generating an output that is easy to interpret, including all maximum likelihood and log-likelihood test results. The only input needed from the user is the CDS and peptide FASTA files of proteins of interest. The pipeline is provided in a Docker image, requiring no program or module installation, enabling the application of the pipeline in any computing environment. AlexandrusPS and its documentation are available via GitHub (https://github.com/alejocn5/AlexandrusPS).

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Gene transcriptional profiles in gonads of Bacillus taxa (Phasmida) with different cytological mechanisms of automictic parthenogenesis

The evolution of automixis - i.e., meiotic parthenogenesis - requires several features, including ploidy restoration after meiosis and maintenance of fertility. Characterizing the relative contribution of novel versus pre-existing genes and the similarities in their expression and sequence evolution is fundamental to understand the evolution of reproductive novelties. Here we identify gonads-biased genes in two Bacillus automictic stick-insects and compare their expression profile and sequence evolution with a bisexual congeneric species. The two parthenogens restore ploidy through different cytological mechanisms: in Bacillus atticus, nuclei derived from the first meiotic division fuse to restore a diploid egg nucleus, while in Bacillus rossius, diploidization occurs in some cells of the haploid blastula through anaphase restitution. Parthenogens' gonads transcriptional program is found to be largely assembled from genes that were already present before the establishment of automixis. The three species transcriptional profiles largely reflect their phyletic relationships, yet we identify a shared core of genes with gonad-biased patterns of expression in parthenogens which are either male gonads-biased in the sexual species or are not differentially expressed there. At the sequence level, just a handful of gonads-biased genes were inferred to have undergone instances of positive selection exclusively in the parthenogen species. This work is the first to explore the molecular underpinnings of automixis in a comparative framework: it delineates how reproductive novelties can be sustained by genes whose origin precedes the establishment of the novelty itself and shows that different meiotic mechanisms of reproduction can be associated with a shared molecular ground plan.

Parthenogenesis doubles the rate of amino acid substitution in whiptail mitochondria

Sexual reproduction is ubiquitous in the natural world, suggesting that sex must have extensive benefits to overcome the cost of males compared to asexual reproduction. One hypothesized advantage of sex with strong theoretical support is that sex plays a role in removing deleterious mutations from the genome. Theory predicts that transitions to asexuality should lead to the suppression of recombination and segregation and, in turn, weakened natural selection, allowing for the accumulation of slightly deleterious mutations. We tested this prediction by estimating the dN/dS ratios in asexual vertebrate lineages in the genus Aspidoscelis using whole mitochondrial genomes from seven asexual and five sexual species. We found higher dN/dS ratios in asexual Aspidoscelis species, indicating that asexual whiptails accumulate nonsynonymous substitutions due to weaker purifying selection. Additionally, we estimated nucleotide diversity and found that asexuals harbor significantly less diversity. Thus, despite their recent origins, slightly deleterious mutations accumulated rapidly enough in asexual lineages to be detected. We provide empirical evidence to corroborate the connection between asexuality and increased amino acid substitutions in asexual vertebrate lineages.

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.

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.

Genome Fractionation and Loss of Heterozygosity in Hybrids and Polyploids: Mechanisms, Consequences for Selection, and Link to Gene Function

Hybridization and genome duplication have played crucial roles in the evolution of many animal and plant taxa. The subgenomes of parental species undergo considerable changes in hybrids and polyploids, which often selectively eliminate segments of one subgenome. However, the mechanisms underlying these changes are not well understood, particularly when the hybridization is linked with asexual reproduction that opens up unexpected evolutionary pathways. To elucidate this problem, we compared published cytogenetic and RNAseq data with exome sequences of asexual diploid and polyploid hybrids between three fish species; Cobitis elongatoides, C. taenia, and C. tanaitica. Clonal genomes remained generally static at chromosome-scale levels but their heterozygosity gradually deteriorated at the level of individual genes owing to allelic deletions and conversions. Interestingly, the impact of both processes varies among animals and genomic regions depending on ploidy level and the properties of affected genes. Namely, polyploids were more tolerant to deletions than diploid asexuals where conversions prevailed, and genomic restructuring events accumulated preferentially in genes characterized by high transcription levels and GC-content, strong purifying selection and specific functions like interacting with intracellular membranes. Although hybrids were phenotypically more similar to C. taenia, we found that they preferentially retained C. elongatoides alleles. This demonstrates that favored subgenome is not necessarily the transcriptionally dominant one. This study demonstrated that subgenomes in asexual hybrids and polyploids evolve under a complex interplay of selection and several molecular mechanisms whose efficiency depends on the organism's ploidy level, as well as functional properties and parental ancestry of the genomic region.

Effects of Selection at Linked Sites on Patterns of Genetic Variability

Patterns of variation and evolution at a given site in a genome can be strongly influenced by the effects of selection at genetically linked sites. In particular, the recombination rates of genomic regions correlate with their amount of within-population genetic variability, the degree to which the frequency distributions of DNA sequence variants differ from their neutral expectations, and the levels of adaptation of their functional components. We review the major population genetic processes that are thought to lead to these patterns, focusing on their effects on patterns of variability: selective sweeps, background selection, associative overdominance, and Hill–Robertson interference among deleterious mutations. We emphasize the difficulties in distinguishing among the footprints of these processes and disentangling them from the effects of purely demographic factors such as population size changes. We also discuss how interactions between selective and demographic processes can significantly affect patterns of variability within genomes.

Genetic diversity and the origins of parthenogenesis in the teiid lizard Aspidoscelis laredoensis

Unisexual vertebrates typically form through hybridization events between sexual species in which reproductive mode transitions occur in the hybrid offspring. This evolutionary history is thought to have important consequences for the ecology of unisexual lineages and their interactions with congeners in natural communities. However, these consequences have proven challenging to study owing to uncertainty about patterns of population genetic diversity in unisexual lineages. Of particular interest is resolving the contribution of historical hybridization events versus post formational mutation to patterns of genetic diversity in nature. Here we use restriction site associated DNA genotyping to evaluate genetic diversity and demographic history in Aspidoscelis laredoensis, a diploid unisexual lizard species from the vicinity of the Rio Grande River in southern Texas and northern Mexico. The sexual progenitor species from which one or more lineages are derived also occur in the Rio Grande Valley region, although patterns of distribution across individual sites are quite variable. Results from population genetic and phylogenetic analyses resolved the major axes of genetic variation in this species and highlight how these match predictions based on historical patterns of hybridization. We also found discordance between results of demographic modelling using different statistical approaches with the genomic data. We discuss these insights within the context of the ecological and evolutionary mechanisms that generate and maintain lineage diversity in unisexual species. As one of the most dynamic, intriguing, and geographically well investigated groups of whiptail lizards, these species hold substantial promise for future studies on the constraints of diversification in unisexual vertebrates.

BASE: A novel workflow to integrate nonubiquitous genes in comparative genomics analyses for selection

Inferring the selective forces that orthologous genes underwent across different lineages can help us understand the evolutionary processes that have shaped their extant diversity and the phenotypes they underlie. The most widespread metric to estimate the selection regimes of coding genes-across sites and phylogenies-is the ratio of nonsynonymous to synonymous substitutions (dN/dS, also known as ω). Nowadays, modern sequencing technologies and the large amount of already available sequence data allow the retrieval of thousands of orthologous genes across large numbers of species. Nonetheless, the tools available to explore selection regimes are not designed to automatically process all genes, and their practical usage is often restricted to the single-copy ones which are found across all species considered (i.e., ubiquitous genes). This approach limits the scale of the analysis to a fraction of single-copy genes, which can be as low as an order of magnitude in respect to those which are not consistently found in all species considered (i.e., nonubiquitous genes). Here, we present a workflow named BASE that-leveraging the CodeML framework-eases the inference and interpretation of gene selection regimes in the context of comparative genomics. Although a number of bioinformatics tools have already been developed to facilitate this kind of analyses, BASE is the first to be specifically designed to allow the integration of nonubiquitous genes in a straightforward and reproducible manner. The workflow-along with all relevant documentation-is available at github.com/for-giobbe/BASE.

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.

Development of microsatellite markers for the geographically parthenogenetic stick insect Phraortes elongatus (Insecta: Phasmatodea)

Many plant and animal species exhibit geographic parthenogenesis, wherein unisexual (= parthenogenetic) lineages are more common in their marginal habitats such as high latitude or altitudes than their closely related bisexual counterparts. The Japanese stick insect, Phraortes elongatus (Thunberg) (Insecta: Phasmatodea), is known as a geographically parthenogenetic species due to the existence of both bisexual and unisexual populations. Here, we developed microsatellite markers to infer the genetic variation among populations of P. elongatus. Totally, 13 primer pairs were developed for the species, and they were tested on 47 samples collected from both a bisexual population and a unisexual population. All 13 loci were polymorphic in the bisexual population, whereas no loci were polymorphic in the unisexual population. The loss of variation in the unisexual population implies automixis with terminal fusion or gamete duplication as the mode of parthenogenesis. The markers developed in this study will be helpful for further comprehensive analysis of the genetic diversity and gene flow between bisexual and parthenogenetic lineages of P. elongatus.

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits-mode of genomic inheritance-influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S-5.8S-28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA-histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction-with potential functional consequences related to the loss of sexual reproduction.

Population Size, Sex and Purifying Selection: Comparative Genomics of Two Sister Taxa of the Wild Yeast Saccharomyces paradoxus

This study uses population genomic data to estimate demographic and selection parameters in two sister lineages of the wild yeast Saccharomyces paradoxus and compare their evolution. We first estimate nucleotide and recombinational diversities in each of the two lineages to infer their population size and frequency of sex and then analyze the rate of mutation accumulation since divergence from their inferred common ancestor to estimate the generation time and efficacy of selection. We find that one of the lineages has significantly higher silent nucleotide diversity and lower linkage disequilibrium, indicating a larger population with more frequent sexual generations. The same lineage also shows shorter generation time and higher efficacy of purifying selection, the latter consistent with the finding of larger population size and more frequent sex. Similar analyses are also performed on the ancestries of individual strains within lineages and we find significant differences between strains implying variation in rates of mitotic cell divisions. Our sample includes some strains originating in the Chernobyl nuclear-accident exclusion zone, which has been subjected to high levels of radiation for nearly 30 years now. We find no evidence, however, for increased rates of mutation. Finally, there is a positive correlation between rates of mutation accumulation and length of growing period, as measured by latitude of the place of origin of strains. Our study illustrates the power of genomic analyses in estimating population and life history parameters and testing predictions based on population genetic theory.

A New Diploid Parthenogenetic Whiptail Lizard from Sonora, Mexico, Is the "Missing Link" in the Evolutionary Transition to Polyploidy

AbstractTransitions between sexual and unisexual reproductive modes have significant consequences for the evolutionary trajectories of species. These transitions have occurred numerous times in vertebrates and are frequently mediated by hybridization events. Triploid unisexual vertebrates are thought to arise through hybridization between individuals of a diploid unisexual lineage and a sexual species, although additional evidence that confirms this mechanism is needed in numerous groups. North American whiptail lizards (Aspidoscelis) are notable for being one of the largest radiations of unisexual vertebrates, and the most diverse group of Aspidoscelis includes numerous triploid lineages that have no known diploid unisexual ancestors. This pattern of "missing" ancestors may result from the short evolutionary life span of unisexual lineages or the selective advantages of polyploidy, or it could suggest that alternative mechanisms of triploid formation are operating in nature. We leverage genomic, morphological, and karyotypic data to describe a new diploid unisexual whiptail and show that it is likely the unisexual progenitor of an extant triploid lineage, A. opatae. We also resolve patterns of polyploidization within the A. sexlineatus species group and test predictions about the phenotypic outcomes of hybridization.

Charting the genomic landscape of seed-free plants

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

Evolutionary Dynamics and Consequences of Parthenogenesis in Vertebrates

Parthenogenesis is asexual reproduction without any required participation from males and, as such, is a null model for sexual reproduction. In a comparative context, we can expand our understanding of the evolution and ecology of sex by investigating the consequences of parthenogenesis. In this review, we examine the theoretical predictions of and empirical results on the evolution of asexual reproduction in vertebrates, focusing on recent studies addressing the origins and geographic spread of parthenogenetic lineages and the genomic consequences of an asexual life history. With advances in computational methods and genome technologies, researchers are poised to make rapid and significant progress in studying the origin and evolution of parthenogenesis in vertebrates, thus providing an important perspective on understanding biodiversity patterns of both asexual and sexual populations.

Multigenerational Mitigating Effects of Ocean Acidification on In Vivo Endpoints, Antioxidant Defense, DNA Damage Response, and Epigenetic Modification in an Asexual Monogonont Rotifer

Ocean acidification (OA) is caused by changes in ocean carbon chemistry due to increased atmospheric pCO₂ and is predicted to have deleterious effects on marine ecosystems. While the potential impacts of OA on many marine species have been studied, the multigenerational effects on asexual organisms remain unknown. We found that low seawater pH induced oxidative stress and DNA damage, decreasing growth rates, fecundity, and lifespans in the parental generation, whereas deleterious effects on in vivo endpoints in F1 and F2 offspring were less evident. The findings suggest that multigenerational adaptive effects play a role in antioxidant abilities and other defense mechanisms. OA-induced DNA damage, including double-strand breaks (DSBs), was fully repaired in F1 offspring of parents exposed to OA for 7 days, indicating that an adaptation mechanism may be the major driving force behind multigenerational adaptive effects. Analysis of epigenetic modification in response to OA involved examination of histone modification of DNA repair genes and a chromatin immunoprecipitation assay, as Bombus koreanus has no methylation pattern for CpG in its genome. We conclude that DSBs, DNA repair, and histone modification play important roles in multigenerational plasticity in response to OA in an asexual monogonont rotifer.

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.

Patterns of Genetic Variability in Genomic Regions with Low Rates of Recombination

The amount of DNA sequence variability in a genomic region is often positively correlated with its rate of crossing over (CO) [1-3]. This pattern is caused by selection acting on linked sites, which reduces genetic variability and biases the frequency distribution of segregating variants toward more rare variants than are expected without selection (skew). These effects may involve the spread of beneficial mutations (selective sweeps [SSWs]), the elimination of deleterious mutations (background selection [BGS]), or both, and are expected to be stronger with lower CO rates [1-3]. However, in a recent study of human populations, the skew was reduced in the lowest CO regions compared with regions with somewhat higher CO rates [4]. A low skew in very low CO regions, compared with theoretical predictions, is seen in the population genomic studies of Drosophila simulans described here and in other Drosophila species. Here, we propose an explanation for lower than expected skew in low CO regions, and validate it using computer simulations; explanations for higher skew with higher CO rates, as in D. simulans, will be explored elsewhere. Partially recessive, linked deleterious mutations can increase neutral variability when the product of the effective population size (N_e) and the selection coefficient against homozygous carriers of mutations (s) is ≤1, i.e., there is associative overdominance (AOD) rather than BGS [5]. AOD can operate in low CO regions, producing a lower skew than in its absence. This opens up a new perspective on how selection affects patterns of variability at linked sites.

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.

Loss and gain of sexual reproduction in the same stick insect

The outcome of competition between different reproductive strategies within a single species can be used to infer selective advantage of the winning strategy. Where multiple populations have independently lost or gained sexual reproduction it is possible to investigate whether the advantage is contingent on local conditions. In the New Zealand stick insect Clitarchus hookeri, three populations are distinguished by recent change in reproductive strategy and we determine their likely origins. One parthenogenetic population has established in the United Kingdom and we provide evidence that sexual reproduction has been lost in this population. We identify the sexual population from which the parthenogenetic population was derived, but show that the UK females have a post‐mating barrier to fertilisation. We also demonstrate that two sexual populations have recently arisen in New Zealand within the natural range of the mtDNA lineage that otherwise characterizes parthenogenesis in this species. We infer independent origins of males at these two locations using microsatellite genotypes. In one population, a mixture of local and nonlocal alleles suggested males were the result of invasion. Males in another population were most probably the result of loss of an X chromosome that produced a male phenotype in situ. Two successful switches in reproductive strategy suggest local competitive advantage for outcrossing over parthenogenetic reproduction. Clitarchus hookeri provides remarkable evidence of repeated and rapid changes in reproductive strategy, with competitive outcomes dependent on local conditions.

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.

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.

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.

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

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