Abortive meiosis in the oogenesis of parthenogenetic Daphnia pulex

Maternal inheritance of functional centrioles in two parthenogenetic nematodes

Centrioles are the core constituent of centrosomes, microtubule-organizing centers involved in directing mitotic spindle assembly and chromosome segregation in animal cells. In sexually reproducing species, centrioles degenerate during oogenesis and female meiosis is usually acentrosomal. Centrioles are retained during male meiosis and, in most species, are reintroduced with the sperm during fertilization, restoring centriole numbers in embryos. In contrast, the presence, origin, and function of centrioles in parthenogenetic species is unknown. We found that centrioles are maternally inherited in two species of asexual parthenogenetic nematodes and identified two different strategies for maternal inheritance evolved in the two species. In Rhabditophanes diutinus, centrioles organize the poles of the meiotic spindle and are inherited by both the polar body and embryo. In Disploscapter pachys, the two pairs of centrioles remain close together and are inherited by the embryo only. Our results suggest that maternally-inherited centrioles organize the embryonic spindle poles and act as a symmetry-breaking cue to induce embryo polarization. Thus, in these parthenogenetic nematodes, centrioles are maternally-inherited and functionally replace their sperm-inherited counterparts in sexually reproducing species.

Efficient CRISPR genome editing and integrative genomic analyses reveal the mosaicism of Cas-induced mutations and pleiotropic effects of scarlet gene in an emerging model system

Despite the revolutionary impacts of CRISPR-Cas gene editing systems, the effective and widespread use of CRISPR technologies in emerging model organisms still faces significant challenges. These include the inefficiency in generating heritable mutations at the organismal level, limited knowledge about the genomic consequences of gene editing, and an inadequate understanding of the inheritance patterns of CRISPR-Cas-induced mutations. This study addresses these issues by 1) developing an efficient microinjection delivery method for CRISPR editing in the microcrustacean Daphnia pulex; 2) assessing the editing efficiency of Cas9 and Cas12a nucleases, examining mutation inheritance patterns, and analyzing the local and global mutation spectrum in the scarlet mutants; and 3) investigating the transcriptomes of scarlet mutants to understand the pleiotropic effects of scarlet underlying their swimming behavior changes. Our reengineered CRISPR microinjection method results in efficient biallelic editing with both nucleases. While indels are dominant in Cas-induced mutations, a few on-site large deletions (>1kb) are observed, most likely caused by microhomology-mediated end joining repair. Knock-in of a stop codon cassette to the scarlet locus was successful, despite complex induced mutations surrounding the target site. Moreover, extensive germline mosaicism exists in some mutants, which unexpectedly produce different phenotypes/genotypes in their asexual progenies. Lastly, our transcriptomic analyses unveil significant gene expression changes associated with scarlet knock-out and altered swimming behavior in mutants, including several genes (e.g., NMDA1, ABAT, CNTNAP2) involved in human neurodegenerative diseases. This study expands our understanding of the dynamics of gene editing in the tractable model organism Daphnia and highlights its promising potential as a neurological disease model.

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.

Transcriptomics and the origin of obligate parthenogenesis

Despite the presence of obligately parthenogenetic (OP) lineages derived from sexual ancestors in diverse phylogenetic groups, the genetic mechanisms giving rise to the OP lineages remain poorly understood. The freshwater microcrustacean Daphnia pulex typically reproduces via cyclical parthenogenesis. However, some populations of OP D. pulex have emerged due to ancestral hybridization and introgression events between two cyclically parthenogenetic (CP) species D. pulex and D. pulicaria. These OP hybrids produce both subitaneous and resting eggs parthenogenetically, deviating from CP isolates where resting eggs are produced via conventional meiosis and mating. This study examines the genome-wide expression and alternative splicing patterns of early subitaneous versus early resting egg production in OP D. pulex isolates to gain insight into the genes and mechanisms underlying this transition to obligate parthenogenesis. Our differential expression and functional enrichment analyses revealed a downregulation of meiosis and cell cycle genes during early resting egg production, as well as divergent expression patterns of metabolism, biosynthesis, and signaling pathways between the two reproductive modes. These results provide important gene candidates for future experimental verification, including the CDC20 gene that activates the anaphase-promoting complex in meiosis.

The Transcriptomic Signature of Cyclical Parthenogenesis

Cyclical parthenogenesis, where females can engage in sexual or asexual reproduction depending on environmental conditions, represents a novel reproductive phenotype that emerged during eukaryotic evolution. The fact that environmental conditions can trigger cyclical parthenogens to engage in distinct reproductive modes strongly suggests that gene expression plays a key role in the origin of cyclical parthenogenesis. However, the genetic basis underlying cyclical parthenogenesis remains understudied. In this study, we characterize the female transcriptomic signature of sexual versus asexual reproduction in the cyclically parthenogenetic microcrustacean Daphnia pulex and Daphnia pulicaria. Our analyses of differentially expressed genes (DEGs), pathway enrichment, and gene ontology (GO) term enrichment clearly show that compared with sexual reproduction, the asexual reproductive stage is characterized by both the underregulation of meiosis and cell cycle genes and the upregulation of metabolic genes. The consensus set of DEGs that this study identifies within the meiotic, cell cycle, and metabolic pathways serves as candidate genes for future studies investigating how the two reproductive cycles in cyclical parthenogenesis are mediated at a molecular level. Furthermore, our analyses identify some cases of divergent expression among gene family members (e.g., doublesex and NOTCH2) associated with asexual or sexual reproductive stage, suggesting potential functional divergence among gene family members.

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.

The maternal effects of dietary restriction on Dnmt expression and reproduction in two clones of Daphnia pulex

The inheritance of epigenetic marks induced by environmental variation in a previous generation is broadly accepted as a mediator of phenotypic plasticity. Transgenerational effects linking maternal experiences to changes in morphology, gene expression, and life history of successive generations are known across many taxa. While the number of studies linking epigenetic variation to ecological maternal effects is increasing rapidly, few if any attempts have been made to investigate molecular mechanisms governing epigenetic functions in the context of ecologically relevant maternal effects. Daphnia make an ideal model for investigating molecular epigenetic mechanisms and ecological maternal effects because they will reproduce asexually in the lab. Daphnia are also known to have strong maternal effects, involving a variety of traits and environmental variables. Using two clones of Daphnia pulex, we investigated the plasticity of life history and DNA methyltransferase (Dnmt) gene expression with respect to food limitation within and across generations. We found strong evidence of genotypic variation of responses of life history and Dnmt expression to low food diets, both within and across generations. In general, effects of offspring diet were larger than either the direct maternal effect or offspring-maternal environment interactions, but the direction of the maternal effect was usually in the opposite direction of the within-generation effect. For both life history and Dnmt expression, we also found that when offspring had low food, effects of the maternal environment were stronger than when offspring had high food.

DNA repair during nonreductional meiosis in the asexual rotifer Adineta vaga

Rotifers of the class Bdelloidea are microscopic animals notorious for their long-term persistence in the apparent absence of sexual reproduction and meiotic recombination. This evolutionary paradox is often counterbalanced by invoking their ability to repair environmentally induced genome breakage. By studying the dynamics of DNA damage response in the bdelloid species Adineta vaga, we found that it occurs rapidly in the soma, producing a partially reassembled genome. By contrast, germline DNA repair is delayed to a specific time window of oogenesis during which homologous chromosomes adopt a meiotic-like juxtaposed configuration, resulting in accurate reconstitution of the genome in the offspring. Our finding that a noncanonical meiosis is the mechanism of germline DNA repair in bdelloid rotifers gives previously unidentified insights on their enigmatic long-term evolution.

Metabolomic analysis predicted changes in growth rate in Daphnia magna exposed to acetaminophen

As urbanization and the global population increases, pollutants associated with municipal wastewater such as pharmaceuticals are becoming more prevalent in aquatic environments. Acetaminophen (paracetamol) is a widely used drug worldwide and one of the most frequently detected pharmaceuticals in freshwater ecosystems. This study investigated the impact of acetaminophen on the metabolite profile of Daphnia magna at two life stages; and used these metabolomic findings to hypothesize a potential impact at a higher organismal level which was subsequently tested experimentally. Targeted polar metabolite analysis by liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to measure changes in the concentration of 51 metabolites in the neonate (> 24 h old) and adult (8 day-old) daphnids following a 48-h exposure to sub-lethal concentrations of acetaminophen. The impact of acetaminophen on the metabolic profile of neonates was widely different from adults. Also, acetaminophen exposure perturbed the abundance of nucleotides more extensively than other metabolites. The acute metabolomic experimental results led to the hypotheses that exposure to sub-lethal concentrations of acetaminophen upregulates protein synthesis in D. magna and subsequently increases growth during early life stages and has an opposite impact on adults. Accordingly, a 10 day growth rate experiment indicated that exposure to acetaminophen elevated biomass production in neonates but not in adults. These novel findings demonstrate that a targeted analysis and interpretation of the changes in the polar metabolic profile of organisms in response to environmental stressors could be used as a tool to predict changes at higher biological levels. As such, this study further emphasizes the incorporation of molecular-level platforms as critical and robust tools in environmental assessment frameworks and biomonitoring of aquatic ecosystems.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

The transcriptomic signature of obligate parthenogenesis

Investigating the origin of parthenogenesis through interspecific hybridization can provide insight into how meiosis may be altered by genetic incompatibilities, which is fundamental for our understanding of the formation of reproductive barriers. Yet the genetic mechanisms giving rise to obligate parthenogenesis in eukaryotes remain understudied. In the microcrustacean Daphnia pulex species complex, obligately parthenogenetic (OP) isolates emerged as backcrosses of two cyclically parthenogenetic (CP) parental species, D. pulex and D. pulicaria, two closely related but ecologically distinct species. We examine the genome-wide expression in OP females at the early resting egg production stage, a life-history stage distinguishing OP and CP reproductive strategies, in comparison to CP females of the same stage from the two parental species. Our analyses of the expression data reveal that underdominant and overdominant genes are abundant in OP isolates, suggesting widespread regulatory incompatibilities between the parental species. More importantly, underdominant genes (i.e., genes with expression lower than both parentals) in the OP isolates are enriched in meiosis and cell-cycle pathways, indicating an important role of underdominance in the origin of obligate parthenogenesis. Furthermore, metabolic and biosynthesis pathways enriched with overdominant genes (i.e., expression higher than both parentals) are another genomic signature of OP isolates.

Genome-Wide Allele-Specific Expression in Obligately Asexual Daphnia pulex and the Implications for the Genetic Basis of Asexuality

Although obligately asexual lineages are thought to experience selective disadvantages associated with reduced efficiency of fixing beneficial mutations and purging deleterious mutations, such lineages are phylogenetically and geographically widespread. However, despite several genome-wide association studies, little is known about the genetic elements underlying the origin of obligate asexuality and how they spread. Because many obligately asexual lineages have hybrid origins, it has been suggested that asexuality is caused by the unbalanced expression of alleles from the hybridizing species. Here, we investigate this idea by identifying genes with allele-specific expression (ASE) in a Daphnia pulex population, in which obligate parthenogens (OP) and cyclical parthenogens (CP) coexist, with the OP clones having been originally derived from hybridization between CP D. pulex and its sister species, Daphnia pulicaria. OP D. pulex have significantly more ASE genes (ASEGs) than do CP D. pulex. Whole-genomic comparison of OP and CP clones revealed ∼15,000 OP-specific markers and 42 consistent ASEGs enriched in marker-defined regions. Ten of the 42 ASEGs have alleles coding for different protein sequences, suggesting functional differences between the products of the two parental alleles. At least three of these ten genes appear to be directly involved in meiosis-related processes, for example, RanBP2 can cause abnormal chromosome segregation in anaphase I, and the presence of Wee1 in immature oocytes leads to failure to enter meiosis II. These results provide a guide for future molecular resolution of the genetic basis of the transition to ameiotic parthenogenesis.

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.

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.

Reproduction stage specific dysregulation of Daphnia magna metabolites as an early indicator of reproductive endocrine disruption

Rapid biomolecular observation in model indicator organisms has been considered as a potential predictor of water pollution from chronic and trace toxicants. This study evaluated the use of Daphnia magna metabolomic measurements as indicators for exposure to reproductive endocrine disruptors by using the model juvenile hormone analogue fenoxycarb. Because D. magna reproduction controls metabolic regulation, the reproduction stage was also carefully considered in metabolic observations and data analysis to examine differences. Comparisons of metabolite abundance regulation between 1 and 12 days of fenoxycarb exposure were performed to investigate the predictability of the sub-chronic (12 days) adverse impacts on reproduction and metabolic regulation based on acute (1 day) metabolic observations. ANOVA-simultaneous component analysis (ASCA) detected reversed patterns in direction of time-course metabolite abundance regulation with fenoxycarb exposure. For example, decreases in the abundances of leucine, asparagine, methionine, and isoleucine which then changed to increases were observed with time during fenoxycarb exposures. The reversed regulation pattern was observed at the last reproduction stage (stage 3), exclusively. Pearson correlation analysis showed that correlations of pairwise metabolites were disrupted with fenoxycarb exposure. Similar to ASCA, data normalization based on the reproduction stage improved the detectability of significant correlations. The disruption on ambient metabolite regulation patterns and pairwise metabolite correlations was consistently observed with both 1 and 12 days of fenoxycarb exposures for sets of select metabolites. The observed regulatory disruptions to these specific metabolites suggest altered oogenesis as the affected metabolites and the specific reproduction stage are related to successful oogenesis. This study demonstrates that D. magna metabolic dysregulation is a predictor of water contamination by endocrine disrupting compounds. The high predictability of sub-chronic (12 days) endocrine disruption was confirmed based on acute (1 day) metabolic observations. Furthermore, integration of the reproduction cycle information in D. magna metabolomics was validated by observing a reproduction stage specific dysregulation in metabolite abundance regulation, which was not observable from the broader data analysis. Consequently, this study confirms the potential for establishing a quantitative relationship between water quality and indicator species metabolic observations. Additionally, it was found that constraining variables relevant to toxicity mechanisms of interest, such as the reproduction stage, is a key consideration for extraction of ecologically meaningful information in environmental metabolomics.

Methyl farnesoate regulatory mechanisms underlying photoperiod-dependent sex determination in the freshwater crustacean Daphnia magna

Freshwater zooplankton Daphnia magna has been widely used in ecotoxicology studies. During the last 20 years, it has been demonstrated that the topical application of juvenile hormone (JH) or JH analogs to mother daphnids induce male offspring production. Based on this finding, an in vivo screening validation method for chemicals with JH agonistic effect has developed. Although this screening system successfully identified a number of JH-like chemicals, molecular mechanisms underlying the male sex-determining process remain largely unknown. To address this issue, we established a reliable male- or female-producing system using Daphnia pulex WTN6 strain by changing the rearing photoperiod. Taking advantage of this rearing system, we successfully found several factors involving male sex determination such as ionotropic glutamate receptors, protein kinase C and pantothenate. Here, we used two D. magna strains that can also control the production of female or male offspring by photoperiod differences as model species for ecotoxicology studies. We demonstrated that either treatment of antagonist of ionotropic glutamate receptors or inhibitor of protein kinase C strongly suppressed male offspring production even under male-producing conditions. Moreover, we revealed that male sex-determining processes are likely diverged between D. magna and D. pulex based on the current experiment. This study provides a fine experimental method for in vivo screening not only JH agonists but also JH antagonists. Moreover, using daphnids with photoperiod-dependent sex determination manner will hugely contribute to understanding the mode-of-action of JH in daphnids.

Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data

Several methods based on the Sequential Markovian coalescence (SMC) have been developed that make use of genome sequence data to uncover population demographic history, which is of interest in its own right and is a key requirement to generate a null model for selection tests. While these methods can be applied to all possible kind of species, the underlying assumptions are sexual reproduction in each generation and non-overlapping generations. However, in many plants, invertebrates, fungi and other taxa, those assumptions are often violated due to different ecological and life history traits, such as self-fertilization or long term dormant structures (seed or egg-banking). We develop a novel SMC-based method to infer 1) the rates/parameters of dormancy and of self-fertilization, and 2) the populations' past demographic history. Using simulated data sets, we demonstrate the accuracy of our method for a wide range of demographic scenarios and for sequence lengths from one to 30 Mb using four sampled genomes. Finally, we apply our method to a Swedish and a German population of Arabidopsis thaliana demonstrating a selfing rate of ca. 0.87 and the absence of any detectable seed-bank. In contrast, we show that the water flea Daphnia pulex exhibits a long lived egg-bank of three to 18 generations. In conclusion, we here present a novel method to infer accurate demographies and life-history traits for species with selfing and/or seed/egg-banks. Finally, we provide recommendations for the use of SMC-based methods for non-model organisms, highlighting the importance of the per site and the effective ratios of recombination over mutation.

Hybridization promotes asexual reproduction in Caenorhabditis nematodes

Although most unicellular organisms reproduce asexually, most multicellular eukaryotes are obligately sexual. This implies that there are strong barriers that prevent the origin or maintenance of asexuality arising from an obligately sexual ancestor. By studying rare asexual animal species we can gain a better understanding of the circumstances that facilitate their evolution from a sexual ancestor. Of the known asexual animal species, many originated by hybridization between two ancestral sexual species. The balance hypothesis predicts that genetic incompatibilities between the divergent genomes in hybrids can modify meiosis and facilitate asexual reproduction, but there are few instances where this has been shown. Here we report that hybridizing two sexual Caenorhabditis nematode species (C. nouraguensis females and C. becei males) alters the normal inheritance of the maternal and paternal genomes during the formation of hybrid zygotes. Most offspring of this interspecies cross die during embryogenesis, exhibiting inheritance of a diploid C. nouraguensis maternal genome and incomplete inheritance of C. becei paternal DNA. However, a small fraction of offspring develop into viable adults that can be either fertile or sterile. Fertile offspring are produced asexually by sperm-dependent parthenogenesis (also called gynogenesis or pseudogamy); these progeny inherit a diploid maternal genome but fail to inherit a paternal genome. Sterile offspring are hybrids that inherit both a diploid maternal genome and a haploid paternal genome. Whole-genome sequencing of individual viable worms shows that diploid maternal inheritance in both fertile and sterile offspring results from an altered meiosis in C. nouraguensis oocytes and the inheritance of two randomly selected homologous chromatids. We hypothesize that hybrid incompatibility between C. nouraguensis and C. becei modifies maternal and paternal genome inheritance and indirectly induces gynogenetic reproduction. This system can be used to dissect the molecular mechanisms by which hybrid incompatibilities can facilitate the emergence of asexual reproduction.

Live observation of the oviposition process in Daphnia magna

In favorable conditions, Daphnia magna undergoes parthenogenesis to increase progeny production in a short time. However, in unfavorable conditions, Daphnia undergoes sexual reproduction instead and produces resting eggs. Here, we report live observations of the oviposition process in Daphnia magna. We observed that the cellular contents flowed irregularly through the narrow egg canal during oviposition. Amorphous ovarian eggs developed an oval shape immediately after oviposition and, eventually, a round shape. Oviposition of resting eggs occurred in a similar way. Based on the observations, we propose that, unlike Drosophila eggs, Daphnia eggs cannot maintain cytoplasmic integrity during oviposition. We also determined that the parthenogenetic eggs were activated within 20 min, as demonstrated by vitelline envelope formation. Therefore, it is plausible that the eggs of Daphnia magna may be activated by squeezing pressure during oviposition.

Parthenogenesis and developmental constraints

The absence of a paternal contribution in an unfertilized ovum presents two developmental constraints against the evolution of parthenogenesis. We discuss the constraint caused by the absence of a centrosome and the one caused by the missing set of chromosomes and how they have been broken in specific taxa. They are examples of only a few well-underpinned examples of developmental constraints acting at macro-evolutionary scales in animals. Breaking of the constraint of the missing chromosomes is the best understood and generally involves rare occasions of drastic changes of meiosis. These drastic changes can be best explained by having been induced, or at least facilitated, by sudden cytological events (e.g., repeated rounds of hybridization, endosymbiont infections, and contagious infections). Once the genetic and developmental machinery is in place for regular or obligate parthenogenesis, shifts to other types of parthenogenesis can apparently rather easily evolve, for example, from facultative to obligate parthenogenesis, or from pseudoarrhenotoky to haplodiploidy. We argue that the combination of the two developmental constraints forms a near-absolute barrier against the gradual evolution from sporadic to obligate or regular facultative parthenogenesis, which can probably explain why the occurrence of the highly advantageous mode of regular facultative parthenogenesis is so rare and entirely absent in vertebrates.

Evolution of Gene Expression during a Transition from Environmental to Genetic Sex Determination

Genetic sex determination (GSD) can evolve from environmental sex determination (ESD) via an intermediate state in which both coexist in the same population. Such mixed populations are found in the crustacean Daphnia magna, where non-male producers (NMP, genetically determined females) coexist with male producers (MP), in which male production is environmentally inducible and can also artificially be triggered by exposure to juvenile hormone. This makes Daphnia magna a rare model species for the study of evolutionary transitions from ESD to GSD. Although the chromosomal location of the NMP-determining mutation has been mapped, the actual genes and pathways involved in the evolution of GSD from ESD remain unknown. Here, we present a transcriptomic analysis of MP and NMP females under control (female producing) and under hormone exposure conditions. We found ∼100 differentially expressed genes between MP and NMP under control conditions. Genes in the NMP-determining chromosome region were especially likely to show such constitutive expression differences. Hormone exposure led to expression changes of an additional ∼100 (MP) to ∼600 (NMP) genes, with an almost systematic upregulation of those genes in NMP. These observations suggest that the NMP phenotype is not determined by a simple "loss-of-function" mutation. Rather, homeostasis of female offspring production under hormone exposure appears to be an active state, tightly regulated by complex mechanisms involving many genes. In a broader view, this illustrates that the evolution of GSD, while potentially initiated by a single mutation, likely leads to secondary integration involving many genes and pathways.

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

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

Use of 5-azacytidine in a proof-of-concept study to evaluate the impact of pre-natal and post-natal exposures, as well as within generation persistent DNA methylation changes in Daphnia

Short-term exposures at critical stages of development can lead to delayed adverse effects long after the initial stressor has been removed, a concept referred to as developmental origin of adult disease. This indicates that organisms' phenotypes may epigenetically reflect their past exposure history as well as reflecting chemicals currently present in their environment. This concept has significant implications for environmental monitoring. However, there is as yet little or no implementation of epigenetics in environmental risk assessment. In a proof-of-principle study we exposed Daphnia magna to 5-azacytidine, a known DNA de-methylating agent. Exposures covered combinations of prenatal and postnatal exposures as well as different exposure durations and recovery stages. Growth, the transcription of genes and levels of metabolites involved in regulating DNA methylation, and methylation levels of several genes were measured. Our data shows that prenatal exposures caused significant changes in the methylome of target genes, indicating that prenatal stages of Daphnia are also susceptible to same level of change as post-natal stages of Daphnia. While the combination of pre- and postnatal exposures caused the most extreme reduction in DNA methylation compared to the control group. Furthermore, some of the changes in the methylation patterns were persistent even after the initial stressor was removed. Our results suggest that epigenetic biomarkers have the potential to be used as indicators of past chemical exposure history of organisms and provide strong support for implementing changes to the current regimes for chemical risk assessment to mimic realistic environmental scenarios.

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'.

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.

Population Genomics of Daphnia pulex

Using data from 83 isolates from a single population, the population genomics of the microcrustacean Daphnia pulex are described and compared to current knowledge for the only other well-studied invertebrate, Drosophila melanogaster These two species are quite similar with respect to effective population sizes and mutation rates, although some features of recombination appear to be different, with linkage disequilibrium being elevated at short ([Formula: see text] bp) distances in D. melanogaster and at long distances in D. pulex The study population adheres closely to the expectations under Hardy-Weinberg equilibrium, and reflects a past population history of no more than a twofold range of variation in effective population size. Fourfold redundant silent sites and a restricted region of intronic sites appear to evolve in a nearly neutral fashion, providing a powerful tool for population genetic analyses. Amino acid replacement sites are predominantly under strong purifying selection, as are a large fraction of sites in UTRs and intergenic regions, but the majority of SNPs at such sites that rise to frequencies [Formula: see text] appear to evolve in a nearly neutral fashion. All forms of genomic sites (including replacement sites within codons, and intergenic and UTR regions) appear to be experiencing an [Formula: see text] higher level of selection scaled to the power of drift in D. melanogaster, but this may in part be a consequence of recent demographic changes. These results establish D. pulex as an excellent system for future work on the evolutionary genomics of natural populations.

Promoter Architecture and Sex-Specific Gene Expression in Daphnia pulex

Large-scale transcription start site (TSS) profiling produces a high-resolution, quantitative picture of transcription initiation and core promoter locations within a genome. However, application of TSS profiling to date has largely been restricted to a small set of prominent model systems. We sought to characterize the cis-regulatory landscape of the water flea Daphnia pulex, an emerging model arthropod that reproduces both asexually (via parthenogenesis) and sexually (via meiosis). We performed Cap Analysis of Gene Expression (CAGE) with RNA isolated from D. pulex within three developmental states: sexual females, asexual females, and males. Identified TSSs were utilized to generate a "Daphnia Promoter Atlas," i.e., a catalog of active promoters across the surveyed states. Analysis of the distribution of promoters revealed evidence for widespread alternative promoter usage in D. pulex, in addition to a prominent fraction of compactly-arranged promoters in divergent orientations. We carried out de novo motif discovery using CAGE-defined TSSs and identified eight candidate core promoter motifs; this collection includes canonical promoter elements (e.g., TATA and Initiator) in addition to others lacking obvious orthologs. A comparison of promoter activities found evidence for considerable state-specific differential gene expression between states. Our work represents the first global definition of transcription initiation and promoter architecture in crustaceans. The Daphnia Promoter Atlas presented here provides a valuable resource for comparative study of cis-regulatory regions in metazoans, as well as for investigations into the circuitries that underpin meiosis and parthenogenesis.

High mutational rates of large-scale duplication and deletion in Daphnia pulex

Knowledge of the genome-wide rate and spectrum of mutations is necessary to understand the origin of disease and the genetic variation driving all evolutionary processes. Here, we provide a genome-wide analysis of the rate and spectrum of mutations obtained in two Daphnia pulex genotypes via separate mutation-accumulation (MA) experiments. Unlike most MA studies that utilize haploid, homozygous, or self-fertilizing lines, D. pulex can be propagated ameiotically while maintaining a naturally heterozygous, diploid genome, allowing the capture of the full spectrum of genomic changes that arise in a heterozygous state. While base-substitution mutation rates are similar to those in other multicellular eukaryotes (about 4 × 10(-9) per site per generation), we find that the rates of large-scale (>100 kb) de novo copy-number variants (CNVs) are significantly elevated relative to those seen in previous MA studies. The heterozygosity maintained in this experiment allowed for estimates of gene-conversion processes. While most of the conversion tract lengths we report are similar to those generated by meiotic processes, we also find larger tract lengths that are indicative of mitotic processes. Comparison of MA lines to natural isolates reveals that a majority of large-scale CNVs in natural populations are removed by purifying selection. The mutations observed here share similarities with disease-causing, complex, large-scale CNVs, thereby demonstrating that MA studies in D. pulex serve as a system for studying the processes leading to such alterations.

Male meiosis in Crustacea: synapsis, recombination, epigenetics and fertility in Daphnia magna

We present the first detailed cytological study of male meiosis in Daphnia (Crustacea: Branchiopoda: Cladocera)—an aquatic microcrustacean with a cyclical parthenogenetic life cycle. Using immunostaining of the testes in Daphnia magna for baseline knowledge, we characterized the different stages of meiotic division and spermiogenesis in relation to the distribution of proteins involved in synapsis, early recombination events and sister chromatid cohesion. We also studied post-translational histone modifications in male spermatocytes, in relation to the dynamic chromatin progression of meiosis. Finally, we applied a DNA fragmentation test to measure sperm quality of D. magna, with respect to levels of inbreeding. As a proxy for fertility, this technique may be used to assess the reproductive health of a sentinel species of aquatic ecosystems. Daphnia proves to be a model species for comparative studies of meiosis that is poised to improve our understanding of the cytological basis of sexual and asexual reproduction.

Hybridization and the Origin of Contagious Asexuality in Daphnia pulex

Hybridization plays a potentially important role in the origin of obligate parthenogenesis (OP) in many organisms. However, it remains controversial whether hybridization directly triggers the transition from sexual reproduction to obligate asexuality or a hybrid genetic background enables asexual species to persist. Furthermore, we know little about the specific genetic elements from the divergent, yet still hybridizing lineages responsible for this transition and how these elements are further spread to create other OP lineages. In this study, we address these questions in Daphnia pulex, where cyclically parthenogenetic (CP) and OP lineages coexist. Ancestry estimates and whole-genome association mapping using 32 OP isolates suggest that a complex hybridization history between the parental species D. pulex and D. pulicaria is responsible for the introgression of a set of 647 D. pulicaria single nucleotide polymorphism alleles that show perfect association with OP. Crossing experiments using males of OP lineages and females of CP lineages strongly support a polygenic basis for OP. Single-sperm analyses show that although normal meiotic recombination occurs in the production of haploid sperm by males of OP lineages, a significant proportion of such sperm are polyploid, suggesting that the spread of asexual elements through these males (i.e., contagious asexuality) is much less efficient than previously envisioned. Although the current Daphnia genome annotation does not provide mechanistic insight into the nature of the asexuality-associated alleles, these alleles should be considered as candidates for future investigations on the genetic underpinnings of OP.

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.

Population genetic dynamics of an invasion reconstructed from the sediment egg bank

Biological invasions are a global issue with far-reaching consequences for single species, communities and whole ecosystems. Our understanding of modes and mechanisms of biological invasions requires knowledge of the genetic processes associated with successful invasions. In many instances, this information is particularly difficult to obtain as the initial phases of the invasion process often pass unnoticed and we rely on inferences from contemporary population genetic data. Here, we combined historic information with the genetic analysis of resting eggs to reconstruct the invasion of Daphnia pulicaria into Lower Lake Constance (LLC) in the 1970s from the resting egg bank in the sediments. We identified the invader as 'European D. pulicaria' originating from meso- and eutrophic lowland lakes and ponds in Central Europe. The founding population was characterized by extremely low genetic variation in the resting egg bank that increased considerably over time. Furthermore, strong evidence for selfing and/or biparental inbreeding was found during the initial phase of the invasion, followed by a drop of selfing rate to low levels in subsequent decades. Moreover, the increase in genetic variation was most pronounced during early stages of the invasion, suggesting additional introductions during this period. Our study highlights that genetic data covering the entire invasion process from its beginning can be crucial to accurately reconstruct the invasion history of a species. We show that propagule banks can preserve such information enabling the study of population genetic dynamics and sources of genetic variation in successful invasive populations.

Evolutionary dynamics of R2 retroelement and insertion inheritance in the genome of bisexual and parthenogenetic Bacillus rossius populations (Insecta Phasmida)

Theoretical and empirical studies have shown differential management of transposable elements in organisms with different reproductive strategies. To investigate this issue, we analysed the R2 retroelement structure and variability in parthenogenetic and bisexual populations of Bacillus rossius stick insects, as well as insertions inheritance in the offspring of parthenogenetic isolates and of crosses. The B. rossius genome hosts a functional (R2Br(fun) ) and a degenerate (R2Br(deg) ) element, their presence correlating with neither reproductive strategies nor population distribution. The median-joining network method indicated that R2Br(fun) duplicates through a multiple source model, while R2Br(deg) is apparently still duplicating via a master gene model. Offspring analyses showed that unisexual and bisexual offspring have a similar number of R2Br-occupied sites. Multiple or recent shifts from gonochoric to parthenogenetic reproduction may explain the observed data. Moreover, insertion frequency spectra show that higher-frequency insertions in unisexual offspring significantly outnumber those in bisexual offspring. This suggests that unisexual offspring eliminate insertions with lower efficiency. A comparison with simulated insertion frequencies shows that inherited insertions in unisexual and bisexual offspring are significantly different from the expectation. On the whole, different mechanisms of R2 elimination in unisexual vs bisexual offspring and a complex interplay between recombination effectiveness, natural selection and time can explain the observed data.

Targeted gene disruption by use of transcription activator-like effector nuclease (TALEN) in the water flea Daphnia pulex

BACKGROUND

The cosmopolitan microcrustacean Daphnia pulex provides a model system for both human health research and monitoring ecosystem integrity. It is the first crustacean to have its complete genome sequenced, an unprecedented ca. 36% of which has no known homologs with any other species. Moreover, D. pulex is ideally suited for experimental manipulation because of its short reproductive cycle, large numbers of offspring, synchronization of oocyte maturation, and other life history characteristics. However, existing gene manipulation techniques are insufficient to accurately define gene functions. Although our previous investigations developed an RNA interference (RNAi) system in D. pulex, the possible time period of functional analysis was limited because the effectiveness of RNAi is transient. Thus, in this study, we developed a genome editing system for D. pulex by first microinjecting transcription activator-like effector nuclease (TALEN) mRNAs into early embryos and then evaluating TALEN activity and mutation phenotypes.

RESULTS

We assembled a TALEN construct specific to the Distal-less gene (Dll), which is a homeobox transcription factor essential for distal limb development in invertebrates and vertebrates, and evaluated its activity in vitro by single-strand annealing assay. Then, we injected TALEN mRNAs into eggs within 1 hour post-ovulation. Injected embryos presented with defects in the second antenna and altered appendage development, and indel mutations were detected in Dll loci, indicating that this technique successfully knocked out the target gene.

CONCLUSIONS

We succeeded, for the first time in D. pulex, in targeted mutagenesis by use of Platinum TALENs. This genome editing technique makes it possible to conduct reverse genetic analysis in D. pulex, making this species an even more appropriate model organism for environmental, evolutionary, and developmental genomics.

Formation and structure of the ephippium (resting egg case) in relation to molting and egg laying in the water flea Daphnia pulex De Geer (Cladocera: Daphniidae)

Resting eggs produced by daphnid species in response to environmental deterioration play an important role in colonizing new habitats or in re-establishing extinct populations. Females lay resting eggs into the space within the dorsal part of their carapace and form an egg case called the ephippium to protect them. Previous studies mainly reported the morphology of the completely formed ephippium and/or the forming ephippium of an uncertain stage. To understand ephippium formation and to clarify key transitions in the formation of resting eggs, we examined the structure and formation of the ephippium in the water flea Daphnia pulex De Geer (Cladocera: Daphniidae) by stereomicroscopy, histology, and scanning electron microscopy. The females used in this study produced resting eggs by obligate parthenogenesis. We divided ephippium formation into four stages based on two molts and a single ovulation, as follows. Stage I begins 13 min after molting in adult females that do not ovulate. In Stage II, immediately after the first molt, a protuberance appears beneath the neck region and the carapace begins to thicken. In Stage III, the resting eggs ovulate and the carapace in the area of the forming ephippium becomes much thicker than the normal carapace and accumulates dark pigmentation. In Stage IV, following the second molt, the female sheds the ephippium with the enclosed resting eggs and forms a new carapace. These stages will provide a useful reference for future studies on resting egg formation.

Development of a microinjection system for RNA interference in the water flea Daphnia pulex

Background

The ubiquitous, freshwater microcrustacean Daphnia pulex provides a model system for both human health research and monitoring ecosystem integrity. It is the first crustacean to have a well annotated, reference genome assembly that revealed an unusually high gene count highlighted by a large gene orphanage,-i.e., previously uncharacterized genes. Daphnia are capable of either clonal or sexual reproduction, making them ideally suited for genetic manipulation, but the establishment of gene manipulation techniques is needed to accurately define gene functions. Although previous investigations developed an RNA interference (RNAi) system for one congener D. magna, these methods are not appropriate for D. pulex because of the smaller size of their early embryos. In these studies, we develop RNAi techniques for D. pulex by first determining the optimum culture conditions of their isolated embryos and then applying these conditions to the development of microinjection techniques and proof-of-principle RNAi experiments.

Results

We found that isolated embryos were best cultured on a 2% agar plate bathed in 60 mM sucrose dissolved in M4 media, providing optimal conditions for microinjections. Then, we injected double-stranded (ds)RNA specific to the Distal-less gene (Dll), which is a homeobox transcription factor essential for limb development in invertebrates and vertebrates. Injected embryos presented with defects in the second antenna and appendage development, and dsRNA induced the degradation of Dll mRNAs, indicating that this technique successfully inhibited transcription of the target gene.

Conclusions

We developed a microinjection system for RNAi studies in D. pulex. These techniques add to the growing genomic toolbox and enhance the genetic tractability of this important model for environmental, evolutionary, and developmental genomics.

The role of hybridization in the origin and spread of asexuality in Daphnia

The molecular mechanisms leading to asexuality remain little understood despite their substantial bearing on why sexual reproduction is dominant in nature. Here, we examine the role of hybridization in the origin and spread of obligate asexuality in Daphnia pulex, arguably the best-documented case of contagious asexuality. Obligately parthenogenetic (OP) clones of D. pulex have traditionally been separated into 'hybrid' (Ldh SF) and 'nonhybrid' (Ldh SS) forms because the lactase dehydrogenase (Ldh) locus distinguishes the cyclically parthenogenetic (CP) lake dwelling Daphnia pulicaria (Ldh FF) from its ephemeral pond dwelling sister species D. pulex (Ldh SS). The results of our population genetic analyses based on microsatellite loci suggest that both Ldh SS and SF OP individuals can originate from the crossing of CP female F1 (D. pulex × D. pulicaria) and backcross with males from OP lineages carrying genes that suppress meiosis specifically in female offspring. In previous studies, a suite of diagnostic markers was found to be associated with OP in Ldh SS D. pulex lineages. Our association mapping supports a similar genetic mechanism for the spread of obligate parthenogenesis in Ldh SF OP individuals. Interestingly, our study shows that CP D. pulicaria carry many of the diagnostic microsatellite alleles associated with obligate parthenogenesis. We argue that the assemblage of mutations that suppress meiosis and underlie obligate parthenogenesis in D. pulex originated due to a unique historical hybridization and introgression event between D. pulex and D. pulicaria.

Comparative transcriptome analysis of obligately asexual and cyclically sexual rotifers reveals genes with putative functions in sexual reproduction, dormancy, and asexual egg production

BACKGROUND

Sexual reproduction is a widely studied biological process because it is critically important to the genetics, evolution, and ecology of eukaryotes. Despite decades of study on this topic, no comprehensive explanation has been accepted that explains the evolutionary forces underlying its prevalence and persistence in nature. Monogonont rotifers offer a useful system for experimental studies relating to the evolution of sexual reproduction due to their rapid reproductive rate and close relationship to the putatively ancient asexual bdelloid rotifers. However, little is known about the molecular underpinnings of sex in any rotifer species.

RESULTS

We generated mRNA-seq libraries for obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of the monogonont rotifer, Brachionus calyciflorus, to identify genes specific to both modes of reproduction. Our differential expression analysis identified receptors with putative roles in signaling pathways responsible for the transition from asexual to sexual reproduction. Differential expression of a specific copy of the duplicated cell cycle regulatory gene CDC20 and specific copies of histone H2A suggest that such duplications may underlie the phenotypic plasticity required for reproductive mode switch in monogononts. We further identified differential expression of genes involved in the formation of resting eggs, a process linked exclusively to sex in this species. Finally, we identified transcripts from the bdelloid rotifer Adineta ricciae that have significant sequence similarity to genes with higher expression in CP strains of B. calyciflorus.

CONCLUSIONS

Our analysis of global gene expression differences between facultatively sexual and exclusively asexual populations of B. calyciflorus provides insights into the molecular nature of sexual reproduction in rotifers. Furthermore, our results offer insight into the evolution of obligate asexuality in bdelloid rotifers and provide indicators important for the use of monogononts as a model system for investigating the evolution of sexual reproduction.

Inventory and phylogenetic analysis of meiotic genes in monogonont rotifers

A long-standing question in evolutionary biology is how sexual reproduction has persisted in eukaryotic lineages. As cyclical parthenogens, monogonont rotifers are a powerful model for examining this question, yet the molecular nature of sexual reproduction in this lineage is currently understudied. To examine genes involved in meiosis, we generated partial genome assemblies for 2 distantly related monogonont species, Brachionus calyciflorus and B. manjavacas. Here we present an inventory of 89 meiotic genes, of which 80 homologs were identified and annotated from these assemblies. Using phylogenetic analysis, we show that several meiotic genes have undergone relatively recent duplication events that appear to be specific to the monogonont lineage. Further, we compare the expression of "meiosis-specific" genes involved in recombination and all annotated copies of the cell cycle regulatory gene CDC20 between obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of B. calyciflorus. We show that "meiosis-specific" genes are expressed in both CP and OP strains, whereas the expression of one of the CDC20 genes is specific to cyclical parthenogenesis. The data presented here provide insights into mechanisms of cyclical parthenogenesis and establish expectations for studies of obligate asexual relatives of monogononts, the bdelloid rotifer lineage.

Molecular cloning of doublesex genes of four cladocera (water flea) species

Background

The gene doublesex (dsx) is known as a key factor regulating genetic sex determination in many organisms. We previously identified two dsx genes (DapmaDsx1 and DapmaDsx2) from a freshwater branchiopod crustacean, Daphnia magna, which are expressed in males but not in females. D. magna produces males by parthenogenesis in response to environmental cues (environmental sex determination) and we showed that DapmaDsx1 expression during embryonic stages is responsible for the male trait development. The D. magna dsx genes are thought to have arisen by a cladoceran-specific duplication; therefore, to investigate evolutionary conservation of sex specific expression of dsx genes and to further assess their functions in the environmental sex determination, we searched for dsx homologs in four closely related cladoceran species.

Results

We identified homologs of both dsx genes from, D. pulex, D. galeata, and Ceriodaphnia dubia, yet only a single dsx gene was found from Moina macrocopa. The deduced amino acid sequences of all 9 dsx homologs contained the DM and oligomerization domains, which are characteristic for all arthropod DSX family members. Molecular phylogenetic analysis suggested that the dsx gene duplication likely occurred prior to the divergence of these cladoceran species, because that of the giant tiger prawn Penaeus monodon is rooted ancestrally to both DSX1 and DSX2 of cladocerans. Therefore, this result also suggested that M. macrocopa lost dsx2 gene secondarily. Furthermore, all dsx genes identified in this study showed male-biased expression levels, yet only half of the putative 5’ upstream regulatory elements are preserved in D. magna and D. pulex.

Conclusions

The all dsx genes of five cladoceran species examined had similar amino acid structure containing highly conserved DM and oligomerization domains, and exhibited sexually dimorphic expression patterns, suggesting that these genes may have similar functions for environmental sex determination in cladocerans.

Spindle assembly and spatial distribution of γ-tubulin during abortive meiosis and cleavage division in the parthenogenetic water flea Daphnia pulex

In most animal species, centrosomes, the main microtubule-organizing centers, usually disintegrate in oocytes during meiosis and are reconstructed from sperm-provided centrioles before the first cleavage division. In parthenogenetic oocytes, however, no sperm-derived centrosome-dependent microtubule nucleation is expected, as fertilization does not occur. The water flea Daphnia pulex undergoes parthenogenesis and sexual reproduction differentially in response to environmental cues. We used immunofluorescence microscopy with anti-α-tubulin and anti-γ-tubulin antibodies to examine spindle formation and the occurrence of centrosomes during parthenogenetic oogenesis and the subsequent cleavage division in D. pulex. The spindle formed in abortive meiosis in parthenogenesis is barrel-shaped and lacks centrosomes, whereas the spindle in the subsequent cleavage division is typically spindle-shaped, with centrosomes. During abortive meiosis, γ-tubulin is localized along the spindle, while in the first cleavage division it is localized only at the spindle poles. Thus, D. pulex should provide a useful comparative model system for elucidating mechanisms of spindle formation and improving our understanding of how evolutionary modification of these mechanisms is involved in the switch from sexual to parthenogenetic reproduction.