Not so clonal asexuals: Unraveling the secret sex life of Artemia parthenogenetica

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.

Recombination in bdelloid rotifer genomes: asexuality, transfer and stress

Bdelloid rotifers constitute a class of microscopic animals living in freshwater habitats worldwide. Several strange features of bdelloids have drawn attention: their ability to tolerate desiccation and other stresses, a lack of reported males across the clade despite centuries of study, and unusually high numbers of horizontally acquired, non-metazoan genes. Genome sequencing is transforming our understanding of their lifestyle and its consequences, while in turn providing wider insights about recombination and genome organisation in animals. Many questions remain, not least how to reconcile apparent genomic signatures of sex with the continued absence of reported males, why bdelloids have so many horizontally acquired genes, and how their remarkable ability to survive stress interacts with recombination and other genomic processes.

Spatio-temporal patterns of ovarian development and VgR gene silencing reduced fecundity in parthenogenetic Artemia

The halophilic zooplankton brine shrimp Artemia has been used as an experimental animal in multidisciplinary studies. However, the reproductive patterns and its regulatory mechanisms in Artemia remain unclear. In this study, the ovarian development process of parthenogenetic Artemia (A. parthenogenetica) was divided into five stages, and oogenesis or egg formation was identified in six phases. The oogenesis mode was assumed to be polytrophic. We also traced the dynamic translocation of candidate germline stem cells (cGSCs) using EdU labelling and elucidated several key cytological events in oogenesis through haematoxylin and eosin staining and fluorescence imaging. Distinguished from the ovary structure of insects and crustaceans, Artemia germarium originated from ovariole buds and are located at the base of the ovarioles. RNA-seq based on five stages of ovarian development identified 2657 upregulated genes related to reproduction by pair-to-pair comparison. Gbb, Dpp, piwi, vasa, nanos, VgA and VgR genes associated with cGSCs recognition and reproductive development were screened and verified using qPCR. Silencing of the VgR gene in A. parthenogenetica (Ap-VgR) at ovarian development Stage II led to a low level of gene expression (less than 10%) within 5 days, which resulted in variations in oogenesis-related gene expression and significantly inhibited vitellogenesis, impeded oocyte maturation, and eventually decreased the number of offspring. In conclusion, we have illustrated the patterns of ovarian development, outlined the key spatio-temporal features of oogenesis and identified the negative impacts of VgR gene knockdown on oogenesis using A. parthenogenetica as an experimental animal. The findings of this study also lay a foundation for the further study of reproductive biology of invertebrates.

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

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

No recombination suppression in asexually produced males of Daphnia pulex

Obligate parthenogenesis (OP) is often thought to evolve by disruption of reductional meiosis and suppression of crossover recombination. In the crustacean Daphnia pulex, OP lineages, which have evolved from cyclical parthenogenetic (CP) ancestors, occasionally produce males that are capable of reductional meiosis. Here, by constructing high-density linkage maps, we find that these males show only slightly and nonsignificantly reduced recombination rates compared to CP males and females. Both meiosis disruption and recombination suppression are therefore sex-limited (or partly so), which speaks against the evolution of OP by disruption of a gene that is essential for meiosis or recombination in both sexes. The findings may be explained by female-limited action of genes that suppress recombination, but previously identified candidate genes are known to be expressed in both sexes. Alternatively, and equally consistent with the data, OP might have evolved through a reuse of the parthenogenesis pathways already present in CP and through their extension to all events of oogenesis. The causal mutations for the CP to OP transition may therefore include mutations in genes involved in oogenesis regulation and may not necessarily be restricted to genes of the "meiosis toolkit." More generally, our study emphasizes that there are many ways to achieve asexuality, and elucidating the possible mechanisms is key to ultimately identify the genes and traits involved.

Evolutionary Potential of Parthenogenesis-Bisexual Lineages within Triploid Apomictic Thelytoky in Cacopsylla ledi (Flor, 1861) (Hemiptera, Psylloidea) in Fennoscandia

A widely accepted hypothesis is that parthenogenesis is an evolutionary dead end since it is selectively advantageous in the short term only but results in lowered diversification rates. Triploid apomictic parthenogenesis might represent an exception, as in favorable environments, triploid females are able to produce rare males and diploid females. The aim of the present study was to analyze the modes of reproduction and their evolutionary implications in the parthenogenetic psyllid Cacopsylla ledi (Flor, 1861) from Fennoscandia. The cytogenetic assessment of ploidy levels and the analysis of the COI haplotype revealed two geographically separated bisexual lineages implying genuine bisexual populations. The southern lineage occurring south of latitude 65° N in Finland showed a COI haplotype different from that of parthenogenetic triploids in the same population but identical to the haplotype of specimens in a genuine bisexual population in the Czech Republic. This allows us to suggest that bisexuals in southern Fennoscandia represent the original bisexual C. ledi. By contrast, in the northern bisexual lineage north of latitude 65° N, rare males and diploid females carried the same haplotype as triploids in the same population, having been produced by the triploids. In the Kola Peninsula, a genuine bisexual population of presumably rare male/diploid female origin was discovered. As this population is geographically isolated from populations of the ancestral bisexual C. ledi, it can develop into a new bisexual species through peripatric speciation during evolution. Our findings demonstrate that apomictic triploid parthenogenesis is not necessarily an evolutionary dead end but is able to lead to the emergence of a new bisexual species of parthenogenetic origin.

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.

ZW sex-chromosome evolution and contagious parthenogenesis in Artemia brine shrimp

Eurasian brine shrimp (genus Artemia) have closely related sexual and asexual lineages of parthenogenetic females, which produce rare males at low frequencies. Although they are known to have ZW chromosomes, these are not well characterized, and it is unclear whether they are shared across the clade. Furthermore, the underlying genetic architecture of the transmission of asexuality, which can occur when rare males mate with closely related sexual females, is not well understood. We produced a chromosome-level assembly for the sexual Eurasian species A. sinica and characterized in detail the pair of sex chromosomes of this species. We combined this new assembly with short-read genomic data for the sexual species A. sp. Kazakhstan and several asexual lineages of A. parthenogenetica, allowing us to perform an in-depth characterization of sex-chromosome evolution across the genus. We identified a small differentiated region of the ZW pair that is shared by all sexual and asexual lineages, supporting the shared ancestry of the sex chromosomes. We also inferred that recombination suppression has spread to larger sections of the chromosome independently in the American and Eurasian lineages. Finally, we took advantage of a rare male, which we backcrossed to sexual females, to explore the genetic basis of asexuality. Our results suggest that parthenogenesis is likely partly controlled by a locus on the Z chromosome, highlighting the interplay between sex determination and asexuality.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Y recombination arrest and degeneration in the absence of sexual dimorphism

Current theory proposes that degenerated sex chromosomes-such as the mammalian Y-evolve through three steps: (i) recombination arrest, linking male-beneficial alleles to the Y chromosome; (ii) Y degeneration, resulting from the inefficacy of selection in the absence of recombination; and (iii) dosage compensation, correcting the resulting low expression of X-linked genes in males. We investigate a model of sex chromosome evolution that incorporates the coevolution of cis and trans regulators of gene expression. We show that the early emergence of dosage compensation favors the maintenance of Y-linked inversions by creating sex-antagonistic regulatory effects. This is followed by degeneration of these nonrecombining inversions caused by regulatory divergence between the X and Y chromosomes. In contrast to current theory, the whole process occurs without any selective pressure related to sexual dimorphism.

Chromosome-level genome assembly reveals homologous chromosomes and recombination in asexual rotifer Adineta vaga

Bdelloid rotifers are notorious as a speciose ancient clade comprising only asexual lineages. Thanks to their ability to repair highly fragmented DNA, most bdelloid species also withstand complete desiccation and ionizing radiation. Producing a well-assembled reference genome is a critical step to developing an understanding of the effects of long-term asexuality and DNA breakage on genome evolution. To this end, we present the first high-quality chromosome-level genome assemblies for the bdelloid Adineta vaga, composed of six pairs of homologous (diploid) chromosomes with a footprint of paleotetraploidy. The observed large-scale losses of heterozygosity are signatures of recombination between homologous chromosomes, either during mitotic DNA double-strand break repair or when resolving programmed DNA breaks during a modified meiosis. Dynamic subtelomeric regions harbor more structural diversity (e.g., chromosome rearrangements, transposable elements, and haplotypic divergence). Our results trigger the reappraisal of potential meiotic processes in bdelloid rotifers and help unravel the factors underlying their long-term asexual evolutionary success.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Diversity of Modes of Reproduction and Sex Determination Systems in Invertebrates, and the Putative Contribution of Genetic Conflict

About eight million animal species are estimated to live on Earth, and all except those belonging to one subphylum are invertebrates. Invertebrates are incredibly diverse in their morphologies, life histories, and in the range of the ecological niches that they occupy. A great variety of modes of reproduction and sex determination systems is also observed among them, and their mosaic-distribution across the phylogeny shows that transitions between them occur frequently and rapidly. Genetic conflict in its various forms is a long-standing theory to explain what drives those evolutionary transitions. Here, we review (1) the different modes of reproduction among invertebrate species, highlighting sexual reproduction as the probable ancestral state; (2) the paradoxical diversity of sex determination systems; (3) the different types of genetic conflicts that could drive the evolution of such different systems.

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.