Automixis in Artemia: solving a century-old controversy

A novel strategy to study apomixis, automixis, and autogamy in plants

KEY MESSAGE

The combination of a flow cytometric seed screen and genotyping of each single seed offers a cost-effective approach to detecting complex reproductive pathways in flowering plants. Reproduction may be seen as one of the driving forces of evolution. Flow cytometric seed screen and genotyping of parents and progeny are commonly employed techniques to discern various modes of reproduction in flowering plants. Nevertheless, both methods possess limitations constraining their individual capacity to investigate reproductive modes thoroughly. We implemented both methods in a novel manner to analyse reproduction pathways using a carefully selected material of parental individuals and their seed progeny. The significant advantage of this approach lies in its ability to apply both methods to a single seed. The introduced methodology provides valuable insights into discerning the levels of apomixis, sexuality, and selfing in complex Rubus taxa. The results may be explained by the occurrence of automixis in Rubus, which warrants further investigation. The approach showcased its effectiveness in a different apomictic system, specifically in Taraxacum. Our study presents a comprehensive methodological approach for determining the mode of reproduction where flow cytometry loses its potential. It provides a reliable and cost-effective method with significant potential in biosystematics, population genetics, and crop breeding.

Evidence for selfing in a vertebrate from whole-genome sequencing

A growing number of recent genomic studies report asexual parthenogenetic reproduction in a wide range of taxa, including vertebrate species from the reptile, bird, and fish lineages. Yet, self-fertilization (selfing) has been recorded only in a single vertebrate, the mangrove killifish Kryptolebias marmoratus In cichlid fishes, sex determination is notably diverse and can be influenced by the environment, and sequential hermaphroditism has been reported for some species. Here, we present evidence for a case of facultative selfing in the cichlid fish Benitochromis nigrodorsalis, which is otherwise known as biparentally reproducing ovophilic mouthbrooder from Western Africa. Our laboratory observations revealed that a wild-caught individual produced repeatedly viable offspring in absence of a mating partner. By analyzing genome-wide single-nucleotide polymorphism (SNP) data, we compare that individual and two of its offspring to shed light on its reproductive mode. First, our results confirm uniparental reproduction. Second, overall heterozygosity is reduced in the offspring compared with outbred individuals. Retained maternal heterozygosity in the offspring is ∼51%, which is close to the theoretically expected value of a heterozygosity reduction of 50% by selfing. Heterozygosity patterns along individual chromosomes do not point to alternative parthenogenetic reproductive mechanisms like automixis by terminal or central fusion. Facultative selfing may represent an adaptive strategy ensuring reproduction when mating partners are absent and, hence, contribute to the cichlids' enormous evolutionary success.

A cytological revisit on parthenogenetic Artemia and the deficiency of a meiosis-specific recombinase DMC1 in the possible transition from bisexuality to parthenogenesis

Although parthenogenesis is widespread in nature and known to have close relationships with bisexuality, the transitional mechanism is poorly understood. Artemia is an ideal model to address this issue because bisexuality and "contagious" obligate parthenogenesis independently exist in its congeneric members. In the present study, we first performed chromosome spreading and immunofluorescence to compare meiotic processes of Artemia adopting two distinct reproductive ways. The results showed that, unlike conventional meiosis in bisexual Artemia, meiosis II in parthenogenic Artemia is entirely absent and anaphase I is followed by a single mitosis-like equational division. Interspecific comparative transcriptomics showed that two central molecules in homologous recombination (HR), Dmc1 and Rad51, exhibited significantly higher expression in bisexual versus parthenogenetic Artemia. qRT-PCR indicated that the expression of both genes peaked at the early oogenesis and gradually decreased afterward. Knocking-down by RNAi of Dmc1 in unfertilized females of bisexual Artemia resulted in a severe deficiency of homologous chromosome pairing and produced univalents at the middle oogenesis stage, which was similar to that of parthenogenic Artemia, while in contrast, silencing Rad51 led to no significant chromosome morphological change. Our results indicated that Dmc1 is vital for HR in bisexual Artemia, and the deficiency of Dmc1 may be correlated with or even possibly one of core factors in the transition from bisexuality to parthenogenesis.

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.

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 Artemia 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 Artemia sp. Kazakhstan and several asexual lineages of Artemia 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.

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.

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.

Genomic signatures of clonality in the deep water kelp Laminaria rodriguezii

The development of population genomic approaches in non-model species allows for renewed studies of the impact of reproductive systems and genetic drift on population diversity. Here, we investigate the genomic signatures of partial clonality in the deep water kelp Laminaria rodriguezii, known to reproduce by both sexual and asexual means. We compared these results with the species Laminaria digitata, a closely related species that differs by different traits, in particular its reproductive mode (no clonal reproduction). We analysed genome-wide variation with dd-RAD sequencing using 4,077 SNPs in L. rodriguezii and 7,364 SNPs in L. digitata. As predicted for partially clonal populations, we show that the distribution of FIS within populations of L. rodriguezii is shifted toward negative values, with a high number of loci showing heterozygote excess. This finding is the opposite of what we observed within sexual populations of L. digitata, characterized by a generalized deficit in heterozygotes. Furthermore, we observed distinct distributions of FIS among populations of L. rodriguezii, which is congruent with the predictions of theoretical models for different levels of clonality and genetic drift. These findings highlight that the empirical distribution of FIS is a promising feature for the genomic study of asexuality in natural populations. Our results also show that the populations of L. rodriguezii analysed here are genetically differentiated and probably isolated. Our study provides a conceptual framework to investigate partial clonality on the basis of RAD-sequencing SNPs. These results could be obtained without any reference genome, and are therefore of interest for various non-model species.

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

The maintenance of sex is paradoxical as sexual species pay the “twofold cost of males” and should thus quickly be replaced by asexual mutants reproducing clonally. However, asexuals may not be strictly clonal and engage in “cryptic sex,” challenging this simple scenario. We study the cryptic sex life of the brine shrimp Artemia parthenogenetica, which has once been termed an “ancient asexual” and where no genetic differences have ever been observed between parents and offspring. This asexual species rarely produces males, which can hybridize with sexual females of closely related species and transmit asexuality to their offspring. Using such hybrids, we show that recombination occurs in asexual lineages, causing loss‐of‐heterozygosity and parent‐offspring differences. These differences cannot generally be observed in field‐sampled asexuals because once heterozygosity is lost, subsequent recombination leaves no footprint. Furthermore, using extensive paternity tests, we show that hybrid females can reproduce both sexually and asexually, and transmit asexuality to both sexually and asexually produced offspring in a dominant fashion. Finally, we show that, contrary to previous reports, field‐sampled asexual females also rarely reproduce sexually (rate ∼2‰). Overall, most previously known facts about Artemia asexuality turned out to be erroneous. More generally, our findings suggest that the evidence for strictly clonal reproduction of asexual species needs to be reconsidered, and that rare sex and consequences of nonclonal asexuality, such as gene flow within asexuals, need to be more widely taken into account in more realistic models for the maintenance of sex and the persistence of asexual lineages.

The discernible and hidden effects of clonality on the genotypic and genetic states of populations: Improving our estimation of clonal rates

Partial clonality is widespread across the tree of life, but most population genetic models are designed for exclusively clonal or sexual organisms. This gap hampers our understanding of the influence of clonality on evolutionary trajectories and the interpretation of population genetic data. We performed forward simulations of diploid populations at increasing rates of clonality (c), analysed their relationships with genotypic (clonal richness, R, and distribution of clonal sizes, Pareto β) and genetic (F_IS and linkage disequilibrium) indices, and tested predictions of c from population genetic data through supervised machine learning. Two complementary behaviours emerged from the probability distributions of genotypic and genetic indices with increasing c. While the impact of c on R and Pareto β was easily described by simple mathematical equations, its effects on genetic indices were noticeable only at the highest levels (c > 0.95). Consequently, genotypic indices allowed reliable estimates of c, while genetic descriptors led to poorer performances when c < 0.95. These results provide clear baseline expectations for genotypic and genetic diversity and dynamics under partial clonality. Worryingly, however, the use of realistic sample sizes to acquire empirical data systematically led to gross underestimates (often of one to two orders of magnitude) of c, suggesting that many interpretations hitherto proposed in the literature, mostly based on genotypic richness, should be reappraised. We propose future avenues to derive realistic confidence intervals for c and show that, although still approximate, a supervised learning method would greatly improve the estimation of c from population genetic data.

Convergent recombination cessation between mating-type genes and centromeres in selfing anther-smut fungi

The degree of selfing has major impacts on adaptability and is often controlled by molecular mechanisms determining mating compatibility. Changes in compatibility systems are therefore important evolutionary events, but their underlying genomic mechanisms are often poorly understood. Fungi display frequent shifts in compatibility systems, and their small genomes facilitate elucidation of the mechanisms involved. In particular, linkage between the pre- and postmating compatibility loci has evolved repeatedly, increasing the odds of gamete compatibility under selfing. Here, we studied the mating-type chromosomes of two anther-smut fungi with unlinked mating-type loci despite a self-fertilization mating system. Segregation analyses and comparisons of high-quality genome assemblies revealed that these two species displayed linkage between mating-type loci and their respective centromeres. This arrangement renders the same improved odds of gamete compatibility as direct linkage of the two mating-type loci under the automictic mating (intratetrad selfing) of anther-smut fungi. Recombination cessation was found associated with a large inversion in only one of the four linkage events. The lack of trans-specific polymorphism at genes located in nonrecombining regions and linkage date estimates indicated that the events of recombination cessation occurred independently in the two sister species. Our study shows that natural selection can repeatedly lead to similar genomic patterns and phenotypes, and that different evolutionary paths can lead to distinct yet equally beneficial responses to selection. Our study further highlights that automixis and gene linkage to centromeres have important genetic and evolutionary consequences, while being poorly recognized despite being present in a broad range of taxa.

Long-term prevalence data reveals spillover dynamics in a multi-host (Artemia), multi-parasite (Microsporidia) community

In the study of multi-host parasites, it is often found that host species contribute asymmetrically to parasite transmission. Yet in natural populations, identifying which hosts contribute to parasite transmission and maintenance is a recurring challenge. Here, we approach this issue by taking advantage of natural variation in the composition of a host community. We studied the brine shrimps Artemia franciscana and Artemia parthenogenetica and their microsporidian parasites Anostracospora rigaudi and Enterocytospora artemiae. Previous laboratory experiments had shown that each host can transmit both parasites, but could not predict their actual contributions to the parasites' maintenance in the field. To resolve this, we gathered long-term prevalence data from a metacommunity of these species. Metacommunity patches could contain either or both of the Artemia host species, so that the presence of the hosts could be linked directly to the persistence of the parasites. First, we show that the microsporidian A. rigaudi is a spillover parasite: it was unable to persist in the absence of its maintenance host A. parthenogenetica. This result was particularly striking, as A. rigaudi displayed both high prevalence (in the field) and high infectivity (when tested in the laboratory) in both hosts. Moreover, the seasonal presence of A. parthenogenetica imposed seasonality on the rate of spillover, causing cyclical pseudo-endemics in the spillover host A. franciscana. Second, while our prevalence data was sufficient to identify E. artemiae as either a spillover or a facultative multi-host parasite, we could not distinguish between the two possibilities. This study supports the importance of studying the community context of multi-host parasites, and demonstrates that in appropriate multi-host systems, sampling across a range of conditions and host communities can lead to clear conclusions about the drivers of parasite persistence.

The repeated emergence of asexuality, the hidden genomes and the role of parthenogenetic rare males in the brine shrimp Artemia

The backbone of this endeavour consists of three major components as they appear in the title. My intention is to summarise, as explicitly as possible, both existing and novel data on the occurrence of parthenogenetic rare males assessing their role in conveying sets of genetic information between asexual strains and sexual Artemia species to and fro. Additionally, an assemblage of strong indications and evidence is quoted aiming to unravel possible scenarios of the repeated emergence of asexuality in the brine shrimp and its significance in evolutionary processes involved in speciation.

Low recombination rates in sexual species and sex-asex transitions

In most sexual, diploid eukaryotes, at least one crossover occurs between each pair of homologous chromosomes during meiosis, presumably in order to ensure proper segregation. Well-known exceptions to this rule are species in which one sex does not recombine and specific chromosomes lacking crossover. We review other possible exceptions, including species with chromosome maps of less than 50 cM in one or both sexes. We discuss the idea that low recombination rates may favour sex-asex transitions, or, alternatively may be a consequence of it. We then show that a yet undescribed species of brine shrimp Artemia from Kazakhstan (A sp. Kazakhstan), the closest known relative of the asexual Artemia parthenogenetica, has one of the shortest genetic linkage maps known. Based on a family of 42 individuals and 589 RAD markers, we find that many linkage groups are considerably shorter than 50 cM, suggesting either no obligate crossover or crossovers concentrated at terminal positions with little effect on recombination. We contrast these findings with the published map of the more distantly related sexual congener, A. franciscana, and conclude that the study of recombination in non-model systems is important to understand the evolutionary causes and consequences of recombination.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Evolutionary mysteries in meiosis

Meiosis is a key event of sexual life cycles in eukaryotes. Its mechanistic details have been uncovered in several model organisms, and most of its essential features have received various and often contradictory evolutionary interpretations. In this perspective, we present an overview of these often 'weird' features. We discuss the origin of meiosis (origin of ploidy reduction and recombination, two-step meiosis), its secondary modifications (in polyploids or asexuals, inverted meiosis), its importance in punctuating life cycles (meiotic arrests, epigenetic resetting, meiotic asymmetry, meiotic fairness) and features associated with recombination (disjunction constraints, heterochiasmy, crossover interference and hotspots). We present the various evolutionary scenarios and selective pressures that have been proposed to account for these features, and we highlight that their evolutionary significance often remains largely mysterious. Resolving these mysteries will likely provide decisive steps towards understanding why sex and recombination are found in the majority of eukaryotes.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Resurrection ecology in Artemia

Resurrection ecology (RE) is a very powerful approach to address a wide range of question in ecology and evolution. This approach rests on using appropriate model systems, and only few are known to be available. In this study, we show that Artemia has multiple attractive features (short generation time, cyst bank and collections, well‐documented phylogeography, and ecology) for a good RE model. We show in detail with a case study how cysts can be recovered from sediments to document the history and dynamics of a biological invasion. We finally discuss with precise examples the many RE possibilities with this model system: adaptation to climate change, to pollution, to parasites, to invaders and evolution of reproductive systems.

Asexual but Not Clonal: Evolutionary Processes in Automictic Populations

Many parthenogenetically reproducing animals produce offspring not clonally but through different mechanisms collectively referred to as automixis. Here, meiosis proceeds normally but is followed by a fusion of meiotic products that restores diploidy. This mechanism typically leads to a reduction in heterozygosity among the offspring compared to the mother. Following a derivation of the rate at which heterozygosity is lost at one and two loci, depending on the number of crossovers between loci and centromere, a number of models are developed to gain a better understanding of basic evolutionary processes in automictic populations. Analytical results are obtained for the expected neutral genetic variation, effective population size, mutation-selection balance, selection with overdominance, the spread of beneficial mutations, and selection on crossover rates. These results are complemented by numerical investigations elucidating how associative overdominance (two off-phase deleterious mutations at linked loci behaving like an overdominant locus) can in some cases maintain heterozygosity for prolonged times, and how clonal interference affects adaptation in automictic populations. These results suggest that although automictic populations are expected to suffer from the lack of gene shuffling with other individuals, they are nevertheless, in some respects, superior to both clonal and outbreeding sexual populations in the way they respond to beneficial and deleterious mutations. Implications for related genetic systems such as intratetrad mating, clonal reproduction, selfing, as well as different forms of mixed sexual and automictic reproduction are discussed.

The ontogeny of tolerance curves: habitat quality vs. acclimation in a stressful environment

Stressful environments affect life-history components of fitness through (i) instantaneous detrimental effects, (ii) historical (carry-over) effects and (iii) history-by-environment interactions, including acclimation effects. The relative contributions of these different responses to environmental stress are likely to change along life, but such ontogenic perspective is often overlooked in studies of tolerance curves, precluding a better understanding of the causes of costs of acclimation, and more generally of fitness in temporally fine-grained environments. We performed an experiment in the brine shrimp Artemia to disentangle these different contributions to environmental tolerance, and investigate how they unfold along life. We placed individuals from three clones of A. parthenogenetica over a range of salinities during a week, before transferring them to a (possibly) different salinity for the rest of their lives. We monitored individual survival at repeated intervals throughout life, instead of measuring survival or performance at a given point in time, as commonly done in acclimation experiments. We then designed a modified survival analysis model to estimate phase-specific hazard rates, accounting for the fact that individuals may share the same treatment for only part of their lives. Our approach allowed us to distinguish effects of salinity on (i) instantaneous mortality in each phase (habitat quality effects), (ii) mortality later in life (history effects) and (iii) their interaction. We showed clear effects of early salinity on late survival and interactions between effects of past and current environments on survival. Importantly, analysis of the ontogenetic dynamics of the tolerance curve reveals that acclimation affects different parts of the curve at different ages. Adopting a dynamical view of the ontogeny of tolerance curve should prove useful for understanding niche limits in temporally changing environments, where the full sequence of environments experienced by an individual determines its overall environmental tolerance, and how it changes throughout life.

What does the geography of parthenogenesis teach us about sex?

Theory predicts that sexual reproduction is difficult to maintain if asexuality is an option, yet sex is very common. To understand why, it is important to pay attention to repeatably occurring conditions that favour transitions to, or persistence of, asexuality. Geographic parthenogenesis is a term that has been applied to describe a large variety of patterns where sexual and related asexual forms differ in their geographic distribution. Often asexuality is stated to occur in a habitat that is, in some sense, marginal, but the interpretation differs across studies: parthenogens might not only predominate near the margin of the sexuals' distribution, but might also extend far beyond the sexual range; they may be disproportionately found in newly colonizable areas (e.g. areas previously glaciated), or in habitats where abiotic selection pressures are relatively stronger than biotic ones (e.g. cold, dry). Here, we review the various patterns proposed in the literature, the hypotheses put forward to explain them, and the assumptions they rely on. Surprisingly, few mathematical models consider geographic parthenogenesis as their focal question, but all models for the evolution of sex could be evaluated in this framework if the (often ecological) causal factors vary predictably with geography. We also recommend broadening the taxa studied beyond the traditional favourites.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Uncovering Cryptic Asexuality in Daphnia magna by RAD Sequencing

The breeding systems of many organisms are cryptic and difficult to investigate with observational data, yet they have profound effects on a species' ecology, evolution, and genome organization. Genomic approaches offer a novel, indirect way to investigate breeding systems, specifically by studying the transmission of genetic information from parents to offspring. Here we exemplify this method through an assessment of self-fertilization vs. automictic parthenogenesis in Daphnia magna. Self-fertilization reduces heterozygosity by 50% compared to the parents, but under automixis, whereby two haploid products from a single meiosis fuse, the expected heterozygosity reduction depends on whether the two meiotic products are separated during meiosis I or II (i.e., central vs. terminal fusion). Reviewing the existing literature and incorporating recombination interference, we derive an interchromosomal and an intrachromosomal prediction of how to distinguish various forms of automixis from self-fertilization using offspring heterozygosity data. We then test these predictions using RAD-sequencing data on presumed automictic diapause offspring of so-called nonmale producing strains and compare them with "self-fertilized" offspring produced by within-clone mating. The results unequivocally show that these offspring were produced by automixis, mostly, but not exclusively, through terminal fusion. However, the results also show that this conclusion was only possible owing to genome-wide heterozygosity data, with phenotypic data as well as data from microsatellite markers yielding inconclusive or even misleading results. Our study thus demonstrates how to use the power of genomic approaches for elucidating breeding systems, and it provides the first demonstration of automictic parthenogenesis in Daphnia.