Inventory and phylogenetic analysis of meiotic genes in monogonont rotifers

Recombination in bdelloid rotifer genomes: asexuality, transfer and stress

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

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.

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.

Divergent Gene Expression Following Duplication of Meiotic Genes in the Stick Insect Clitarchus hookeri

Some animal groups, such as stick insects (Phasmatodea), have repeatedly evolved alternative reproductive strategies, including parthenogenesis. Genomic studies have found modification of the genes underlying meiosis exists in some of these animals. Here we examine the evolution of copy number, evolutionary rate, and gene expression in candidate meiotic genes of the New Zealand geographic parthenogenetic stick insect Clitarchus hookeri. We characterized 101 genes from a de novo transcriptome assembly from female and male gonads that have homology with meiotic genes from other arthropods. For each gene we determined copy number, the pattern of gene duplication relative to other arthropod orthologs, and the potential for meiosis-specific expression. There are five genes duplicated in C. hookeri, including one also duplicated in the stick insect Timema cristinae, that are not or are uncommonly duplicated in other arthropods. These included two sister chromatid cohesion associated genes (SA2 and SCC2), a recombination gene (HOP1), an RNA-silencing gene (AGO2) and a cell-cycle regulation gene (WEE1). Interestingly, WEE1 and SA2 are also duplicated in the cyclical parthenogenetic aphid Acyrthosiphon pisum and Daphnia duplex, respectively, indicating possible roles in the evolution of reproductive mode. Three of these genes (SA2, SCC2, and WEE1) have one copy displaying gonad-specific expression. All genes, with the exception of WEE1, have significantly different nonsynonymous/synonymous ratios between the gene duplicates, indicative of a shift in evolutionary constraints following duplication. These results suggest that stick insects may have evolved genes with novel functions in gamete production by gene duplication.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

A Reappraisal of the Evolutionary and Developmental Pathway of Apomixis and Its Genetic Control in Angiosperms

Apomixis sensu stricto (agamospermy) is asexual reproduction by seed. In angiosperms it represents an easy byway of life cycle renewal through gamete-like cells that give rise to maternal embryos without ploidy reduction (meiosis) and ploidy restitution (syngamy). The origin of apomixis still represents an unsolved problem, as it may be either evolved from sex or the other way around. This review deals with a reappraisal of the origin of apomixis in order to deepen knowledge on such asexual mode of reproduction which seems mainly lacking in the most basal angiosperm orders (i.e., Amborellales, Nymphaeales and Austrobaileyales, also known as ANA-grade), while it clearly occurs in different forms and variants in many unrelated families of monocots and eudicots. Overall findings strengthen the hypothesis that apomixis as a whole may have evolved multiple times in angiosperm evolution following different developmental pathways deviating to different extents from sexuality. Recent developments on the genetic control of apomixis in model species are also presented and adequately discussed in order to shed additional light on the antagonist theories of gain- and loss-of-function over sexuality.

Embryo bioassays with aquatic animals for toxicity testing and hazard assessment of emerging pollutants: A review

This review article gathers the available information on the use of embryo-tests as high-throughput tools for toxicity screening, hazard assessment and prioritization of new and existing chemical compounds. The approach is contextualized considering the new legal trends for animal experimentation, fostering the 3R policy, with reduction of experimental animals, addressing the potential of embryo-tests as high-throughput toxicity screening and prioritizing tools. Further, the current test guidelines, such as the ones provided by OECD and EPA, focus mainly in a limited number of animal lineages, particularly vertebrates and arthropods. To extrapolate hazard assessment to the ecosystem scale, a larger diversity of taxa should be tested. The use of new experimental animal models in toxicity testing, from a representative set of taxa, was thoroughly revised and discussed in this review. Here, we critically review current tools and the main advantages and drawbacks of different animal models and set researcher priorities.

Did apomixis evolve from sex or was it the other way around?

In angiosperms, there are two pathways of reproduction through seeds: sexual, or amphimictic, and asexual, or apomictic. The essential feature of apomixis is that an embryo in an ovule is formed autonomously. It may form from a cell of the nucellus or integuments in an otherwise sexual ovule, a process referred to as adventitious embryony. Alternatively, the embryo may form by parthenogenesis from an unreduced egg that forms in an unreduced embryo sac. The latter may form from an ameiotic megasporocyte, in which case it is referred to as diplospory, or from a cell of the nucellus or integument, in which case it is referred to as apospory. Progeny of apomictic plants are generally identical to the mother plant. Apomixis has been seen over the years as either a gain- or loss-of-function over sexuality, implying that the latter is the default condition. Here, we consider an additional point of view, that apomixis may be anciently polyphenic with sex and that both reproductive phenisms involve anciently canalized components of complex molecular processes. This polyphenism viewpoint suggests that apomixis fails to occur in obligately sexual eukaryotes because genetic or epigenetic modifications have silenced the primitive sex apomixis switch and/or disrupted molecular capacities for apomixis. In eukaryotes where sex and apomixis are clearly polyphenic, apomixis exponentially drives clonal fecundity during reproductively favorable conditions, while stress induces sex for stress-tolerant spore or egg formation. The latter often guarantees species survival during environmentally harsh seasons.

Repeated Evolution of Asexuality Involves Convergent Gene Expression Changes

Asexual reproduction has evolved repeatedly from sexual ancestors across a wide range of taxa. Whereas the costs and benefits associated with asexuality have received considerable attention, the molecular changes underpinning the evolution of asexual reproduction remain relatively unexplored. In particular, it is completely unknown whether the repeated evolution of asexual phenotypes involves similar molecular changes, as previous studies have focused on changes occurring in single lineages. Here, we investigate the extent of convergent gene expression changes across five independent transitions to asexuality in stick insects. We compared gene expression of asexual females to females of close sexual relatives in whole-bodies, reproductive tracts, and legs. We identified a striking amount of convergent gene expression change (up to 8% of genes), greatly exceeding that expected by chance. Convergent changes were also tissue-specific, and most likely driven by selection for functional changes. Genes showing convergent changes in the reproductive tract were associated with meiotic spindle formation and centrosome organization. These genes are particularly interesting as they can influence the production of unreduced eggs, a key barrier to asexual reproduction. Changes in legs and whole-bodies were likely involved in female sexual trait decay, with enrichment in terms such as sperm-storage and pigmentation. By identifying changes occurring across multiple independent transitions to asexuality, our results provide a rare insight into the molecular basis of asexual phenotypes and suggest that the evolutionary path to asexuality is highly constrained, requiring repeated changes to the same key genes.

Evolutionary diversity and novelty of DNA repair genes in asexual Bdelloid rotifers

Background

Bdelloid rotifers are the oldest, most diverse and successful animal taxon for which males, hermaphrodites, and traditional meiosis are unknown. Their degenerate tetraploid genome, with 2–4 copies of most loci, includes thousands of genes acquired from all domains of life by horizontal transfer. Many bdelloid species thrive in ephemerally aquatic habitats by surviving desiccation at any life stage with no loss of fecundity or lifespan. Their unique genomic diversity and the intense selective pressure of desiccation provide an exceptional opportunity to study the evolution of diversity and novelty in genes involved in DNA repair.

Results

We used genomic data and RNA-Seq of the desiccation process in the bdelloid Adineta vaga to characterize DNA damage reversal, translesion synthesis, and the major DNA repair pathways: base, nucleotide, and alternate excision repair, mismatch repair (MMR), and double strand break repair by homologous recombination (HR) and classical non-homologous end joining (NHEJ). We identify multiple horizontally transferred DNA damage response genes otherwise unknown in animals (AlkD, Fpg, LigK UVDE), and the presence of genes often considered vertebrate specific, particularly in the NHEJ complex and X family polymerases. While 75–100% of genes involved in MMR and HR are present in 0–2 copies, genes involved in NHEJ, which are present in only a single copy in nearly all other animals, are retained in 3–8 copies. We present structural predictions and expression evidence of neo- or sub-functionalization of multiple copy genes involved in NHEJ and other repair processes.

Conclusion

The horizontally-acquired genes and duplicated genes in BER and NHEJ suggest resilience to oxidative damage is conferred in part by increased DNA damage recognition and efficient end repair capabilities. The pattern of gene loss and retention in MMR and HR may facilitate recombination and gene conversion between divergent sequences, thus providing at least some of the benefits of sex. The unique retention and divergence of duplicates genes in NHEJ may be facilitated by the lack of efficient selection in the absence of meiotic recombination and independent assortment, and may contribute to the evolutionary success of bdelloids.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1288-9) contains supplementary material, which is available to authorized users.

Retention of Core Meiotic Genes Across Diverse Hymenoptera

The cellular mechanisms of meiosis are critical for proper gamete formation in sexual organisms. Functional studies in model organisms have identified genes essential for meiosis, yet the extent to which this core meiotic machinery is conserved across non-model systems is not fully understood. Moreover, it is unclear whether deviation from canonical modes of sexual reproduction is accompanied by modifications in the genetic components involved in meiosis. We used a robust approach to identify and catalogue meiosis genes in Hymenoptera, an insect order typically characterized by haplodiploid reproduction. Using newly available genome data, we searched for 43 genes involved in meiosis in 18 diverse hymenopterans. Seven of eight genes with roles specific to meiosis were found across a majority of surveyed species, suggesting the preservation of core meiotic machinery in haplodiploid hymenopterans. Phylogenomic analyses of the inventory of meiosis genes and the identification of shared gene duplications and losses provided support for the grouping of species within Proctotrupomorpha, Ichneumonomorpha, and Aculeata clades, along with a paraphyletic Symphyta. The conservation of meiosis genes across Hymenoptera provides a framework for studying transitions between reproductive modes in this insect group.

Comparative genomics of bdelloid rotifers: Insights from desiccating and nondesiccating species

Bdelloid rotifers are a class of microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasised desiccation as the cause of several unusual genomic features in this group. Here, we present high-quality whole-genome sequences of 3 bdelloid species: Rotaria macrura and R. magnacalcarata, which are both desiccation intolerant, and Adineta ricciae, which is desiccation tolerant. In combination with the published assembly of A. vaga, which is also desiccation tolerant, we apply a comparative genomics approach to evaluate the potential effects of desiccation tolerance and asexuality on genome evolution in bdelloids. We find that ancestral tetraploidy is conserved among all 4 bdelloid species, but homologous divergence in obligately aquatic Rotaria genomes is unexpectedly low. This finding is contrary to current models regarding the role of desiccation in shaping bdelloid genomes. In addition, we find that homologous regions in A. ricciae are largely collinear and do not form palindromic repeats as observed in the published A. vaga assembly. Consequently, several features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred nonmetazoan genes in both desiccating and nondesiccating bdelloid species and show that this unusual feature is not shared by other animal phyla, even those with desiccation-tolerant representatives. These comparisons call into question the proposed role of desiccation in mediating horizontal genetic transfer.

Identification of the meiotic toolkit in diatoms and exploration of meiosis-specific SPO11 and RAD51 homologs in the sexual species Pseudo-nitzschia multistriata and Seminavis robusta

BACKGROUND

Sexual reproduction is an obligate phase in the life cycle of most eukaryotes. Meiosis varies among organisms, which is reflected by the variability of the gene set associated to the process. Diatoms are unicellular organisms that belong to the stramenopile clade and have unique life cycles that can include a sexual phase.

RESULTS

The exploration of five diatom genomes and one diatom transcriptome led to the identification of 42 genes potentially involved in meiosis. While these include the majority of known meiosis-related genes, several meiosis-specific genes, including DMC1, could not be identified. Furthermore, phylogenetic analyses supported gene identification and revealed ancestral loss and recent expansion in the RAD51 family in diatoms. The two sexual species Pseudo-nitzschia multistriata and Seminavis robusta were used to explore the expression of meiosis-related genes: RAD21, SPO11-2, RAD51-A, RAD51-B and RAD51-C were upregulated during meiosis, whereas other paralogs in these families showed no differential expression patterns, suggesting that they may play a role during vegetative divisions. An almost identical toolkit is shared among Pseudo-nitzschia multiseries and Fragilariopsis cylindrus, as well as two species for which sex has not been observed, Phaeodactylum tricornutum and Thalassiosira pseudonana, suggesting that these two may retain a facultative sexual phase.

CONCLUSIONS

Our results reveal the conserved meiotic toolkit in six diatom species and indicate that Stramenopiles share major modifications of canonical meiosis processes ancestral to eukaryotes, with important divergences in each Kingdom.

Meiosis gene inventory of four ciliates reveals the prevalence of a synaptonemal complex-independent crossover pathway

To establish which meiosis genes are present in ciliates, and to look for clues as to which recombination pathways may be treaded by them, four genomes were inventoried for 11 meiosis-specific and 40 meiosis-related genes. We found that the set of meiosis genes shared by Tetrahymena thermophila, Paramecium tetraurelia, Ichthyophthirius multifiliis, and Oxytricha trifallax is consistent with the prevalence of a Mus81-dependent class II crossover pathway that is considered secondary in most model eukaryotes. There is little evidence for a canonical class I crossover pathway that requires the formation of a synaptonemal complex (SC). This gene inventory suggests that meiotic processes in ciliates largely depend on mitotic repair proteins for executing meiotic recombination. We propose that class I crossovers and SCs were reduced sometime during the evolution of ciliates. Consistent with this reduction, we provide microscopic evidence for the presence only of degenerate SCs in Stylonychia mytilus. In addition, lower nonsynonymous to synonymous mutation rates of some of the meiosis genes suggest that, in contrast to most other nuclear genes analyzed so far, meiosis genes in ciliates are largely evolving at a slower rate than those genes in fungi and animals.

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.