Sexual species are separated by larger genetic gaps than asexual species in rotifers

Bdelloid rotifers deploy horizontally acquired biosynthetic genes against a fungal pathogen

Coevolutionary antagonism generates relentless selection that can favour genetic exchange, including transfer of antibiotic synthesis and resistance genes among bacteria, and sexual recombination of disease resistance alleles in eukaryotes. We report an unusual link between biological conflict and DNA transfer in bdelloid rotifers, microscopic animals whose genomes show elevated levels of horizontal gene transfer from non-metazoan taxa. When rotifers were challenged with a fungal pathogen, horizontally acquired genes were over twice as likely to be upregulated as other genes - a stronger enrichment than observed for abiotic stressors. Among hundreds of upregulated genes, the most markedly overrepresented were clusters resembling bacterial polyketide and nonribosomal peptide synthetases that produce antibiotics. Upregulation of these clusters in a pathogen-resistant rotifer species was nearly ten times stronger than in a susceptible species. By acquiring, domesticating, and expressing non-metazoan biosynthetic pathways, bdelloids may have evolved to resist natural enemies using antimicrobial mechanisms absent from other animals.

Ice-inhabiting species of Bdelloidea Rotifera reveal a pre-Quaternary ancestry in the Arctic cryosphere

Historical climate data indicate that the Earth has passed through multiple geological periods with much warmer-than-present climates, including epochs of the Miocene (23-5.3 mya BP) with temperatures 3-4°C above present, and more recent interglacial stages of the Quaternary, for example, Marine Isotope Stage 11c (approx. 425-395 ka BP) and Middle Holocene thermal maximum (7.5-4.2 ka BP), during which continental glaciers may have melted entirely. Such warm periods would have severe consequences for ice-obligate fauna in terms of their distribution, biodiversity and population structure. To determine the impacts of these climatic events in the Nordic cryosphere, we surveyed ice habitats throughout mainland Norway and Svalbard ranging from maritime glaciers to continental ice patches (i.e. non-flowing, inland ice subjected to deep freezing overwinter), finding particularly widespread populations of ice-inhabiting bdelloid rotifers. Combined mitochondrial and nuclear DNA sequencing identified approx. 16 undescribed, species-level rotifer lineages that revealed an ancestry predating the Quaternary (> 2.58 mya). These rotifers also displayed robust freeze/thaw tolerance in laboratory experiments. Collectively, these data suggest that extensive ice refugia, comparable with stable ice patches across the contemporary Norwegian landscape, persisted in the cryosphere over geological time, and may have facilitated the long-term survival of ice-obligate Metazoa before and throughout the Quaternary.

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.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

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.

Diversification rates in Tardigrada indicate a decreasing tempo of lineage splitting regardless of reproductive mode

Understanding the dynamics of speciation and extinction events is one of the most interesting subjects in evolutionary biology that relates to all life forms, even the smallest ones. Tardigrades are microscopic invertebrates that attracted public and scientific attention mostly due to their ability to enter into the diapause stage called cryptobiosis and in such stage resist extremely harsh environmental conditions. However, although recent research solved a considerable number of phylogenetic uncertainties and further uncovered physiological mechanisms of cryptobiosis, not much attention is given to the evolutionary forces shaping tardigrade diversity. Here, we investigated the effect of reproductive mode on diversification rates in tardigrades using three groups: macrobiotids, echiniscids and milnesids, which represent low, moderate and high levels of parthenogenesis, respectively. Our results showed a decreasing tempo of diversification events for each of the studied groups without any differences that could be ascribed to reproductive mode. We discussed the observed lack of effect in tardigrades acknowledging deficiencies in available data sets and encouraging further studies to understand whether our results can be considered reliable.

Thainema gen. nov. (Leptolyngbyaceae, Synechococcales): A new genus of simple trichal cyanobacteria isolated from a solar saltern environment in Thailand

Simple trichal types constitute a group of cyanobacteria with an abundance of novel, often cryptic taxa. Here, we investigated material collected from wet surface-soil in a saline environment in Petchaburi Province, central Thailand. A morphological comparison of the isolated strain with similar known species, as well as its phylogenetic and species delimitation analyses based on the combined datasets of other related organisms, especially simple trichal cyanobacteria, revealed that the material of this study represented an independent taxon. Using a multifaceted method, we propose that this material represents a new genus, Thainema gen. nov., belonging to the family Leptolyngbyaceae, with the type species Thainema salinarum sp. nov. This novel taxon shares similar ecological habitats with strains previously placed in the same lineage.

Parthenogenesis is self-destructive for scaled reptiles

Parthenogenesis is rare in nature. With 39 described true parthenogens, scaled reptiles (Squamata) are the only vertebrates that evolved this reproductive strategy. Parthenogenesis is ecologically advantageous in the short term, but the young age and rarity of parthenogenetic species indicate it is less advantageous in the long term. This suggests parthenogenesis is self-destructive: it arises often but is lost due to increased extinction rates, high rates of reversal or both. However, this role of parthenogenesis as a self-destructive trait remains unknown. We used a phylogeny of Squamata (5388 species), tree metrics, null simulations and macroevolutionary scenarios of trait diversification to address the factors that best explain the rarity of parthenogenetic species. We show that parthenogenesis can be considered as self-destructive, with high extinction rates mainly responsible for its rarity in nature. Since these parthenogenetic species occur, this trait should be ecologically relevant in the short term.

Genomic signatures of recombination in a natural population of the bdelloid rotifer Adineta vaga

Sexual reproduction is almost ubiquitous among extant eukaryotes. As most asexual lineages are short-lived, abandoning sex is commonly regarded as an evolutionary dead end. Still, putative anciently asexual lineages challenge this view. One of the most striking examples are bdelloid rotifers, microscopic freshwater invertebrates believed to have completely abandoned sexual reproduction tens of Myr ago. Here, we compare whole genomes of 11 wild-caught individuals of the bdelloid rotifer Adineta vaga and present evidence that some patterns in its genetic variation are incompatible with strict clonality and lack of genetic exchange. These patterns include genotype proportions close to Hardy-Weinberg expectations within loci, lack of linkage disequilibrium between distant loci, incongruent haplotype phylogenies across the genome, and evidence for hybridization between divergent lineages. Analysis of triallelic sites independently corroborates these findings. Our results provide evidence for interindividual genetic exchange and recombination in A. vaga, a species previously thought to be anciently asexual.

DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots

DNA barcoding identification needs a good characterization of intraspecific genetic divergence to establish the limits between species. Yet, the number of barcodes per species is many times low and geographically restricted. A poor coverage of the species distribution range may hamper identification, especially when undersampled areas host genetically distinct lineages. If so, the genetic distance between some query sequences and reference barcodes may exceed the maximum intraspecific threshold for unequivocal species assignation. Taking a group of Quercus herbivores (moths) in Europe as model system, we found that the number of DNA barcodes from southern Europe is proportionally very low in the Barcoding of Life Data Systems. This geographical bias complicates the identification of southern query sequences, due to their high intraspecific genetic distance with respect to barcodes from higher latitudes. Pairwise intraspecific genetic divergence increased along with spatial distance, but was higher when at least one of the sampling sites was in southern Europe. Accordingly, GMYC (General Mixed Yule Coalescent) single-threshold model retrieved clusters constituted exclusively by Iberian haplotypes, some of which could correspond to cryptic species. The number of putative species retrieved was more reliable than that of multiple-threshold GMYC but very similar to results from ABGD and jMOTU. Our results support GMYC as a key resource for species delimitation within poorly inventoried biogeographic regions in Europe, where historical factors (e.g., glaciations) have promoted genetic diversity and singularity. Future European DNA barcoding initiatives should be preferentially performed along latitudinal gradients, with special focus on southern peninsulas.

Outstanding intraindividual genetic diversity in fissiparous planarians (Dugesia, Platyhelminthes) with facultative sex

BACKGROUND

Predicted genetic consequences of asexuality include high intraindividual genetic diversity (i.e., the Meselson effect) and accumulation of deleterious mutations (i.e., Muller's Ratchet), among others. These consequences have been largely studied in parthenogenetic organisms, but studies on fissiparous species are scarce. Differing from parthenogens, fissiparous organisms inherit part of the soma of the progenitor, including somatic mutations. Thus, in the long term, fissiparous reproduction may also result in genetic mosaicism, besides the presence of the Meselson effect and Muller's Ratchet. Dugesiidae planarians show outstanding regeneration capabilities, allowing them to naturally reproduce by fission, either strictly or combined with sex (facultative). Therefore, they are an ideal model to analyze the genetic footprint of fissiparous reproduction, both when it is alternated with sex and when it is the only mode of reproduction.

RESULTS

In the present study, we generate and analyze intraindividual cloned data of a nuclear and a mitochondrial gene of sexual, fissiparous and facultative wild populations of the species Dugesia subtentaculata. We find that most individuals, independently of their reproductive strategy, are mosaics. However, the intraindividual haplotype and nucleotide diversity of fissiparous and facultative individuals is significantly higher than in sexual individuals, with no signs of Muller's Ratchet. Finally, we also find that this high intraindividual genetic diversity of fissiparous and facultative individuals is composed by different combinations of ancestral and derived haplotypes of the species.

CONCLUSIONS

The intraindividual analyses of genetic diversity point out that fissiparous reproduction leaves a very special genetic footprint in individuals, characterized by mosaicism combined with the Meselson effect (named in the present study as the mosaic Meselson effect). Interestingly, the different intraindividual combinations of ancestral and derivate genetic diversity indicate that haplotypes generated during periods of fissiparous reproduction can be also transmitted to the progeny through sexual events, resulting in offspring showing a wide range of genetic diversity and putatively allowing purifying selection to act at both intraindividual and individual level. Further investigations, using Dugesia planarians as model organisms, would be of great value to delve into this new model of genetic evolution by the combination of fission and sex.

Importance of Resolving Fungal Nomenclature: the Case of Multiple Pathogenic Species in the Cryptococcus Genus

Cryptococcosis is a major fungal disease caused by members of the Cryptococcus gattii and Cryptococcus neoformans species complexes. After more than 15 years of molecular genetic and phenotypic studies and much debate, a proposal for a taxonomic revision was made. The two varieties within C. neoformans were raised to species level, and the same was done for five genotypes within C. gattii. In a recent perspective (K. J. Kwon-Chung et al., mSphere 2:e00357-16, 2017, https://doi.org/10.1128/mSphere.00357-16), it was argued that this taxonomic proposal was premature and without consensus in the community. Although the authors of the perspective recognized the existence of genetic diversity, they preferred the use of the informal nomenclature "C. neoformans species complex" and "C. gattii species complex." Here we highlight the advantage of recognizing these seven species, as ignoring these species will impede deciphering further biologically and clinically relevant differences between them, which may in turn delay future clinical advances.

Genetic Exchange among Bdelloid Rotifers Is More Likely Due to Horizontal Gene Transfer Than to Meiotic Sex

Although strict asexuality is supposed to be an evolutionary dead end, morphological, cytogenetic, and genomic data suggest that bdelloid rotifers, a clade of microscopic animals, have persisted and diversified for more than 60 Myr in an ameiotic fashion. Moreover, the genome of bdelloids of the genus Adineta comprises 8%-10% of genes of putative non-metazoan origin, indicating that horizontal gene transfers are frequent within this group and suggesting that this mechanism may also promote genetic exchanges among bdelloids as well. To test this hypothesis, we used five independent sequence markers to study the genetic diversity of 576 Adineta vaga individuals from a park in Belgium. Haplowebs and GMYC analyses revealed the existence of six species among our sampled A. vaga individuals, with strong evidence of both intra- and interspecific recombination. Comparison of genomic regions of three allele-sharing individuals further revealed signatures of genetic exchanges scattered among regions evolving asexually. Our findings suggest that bdelloids evolve asexually but exchange DNA horizontally both within and between species.

Horizontal gene transfer in bdelloid rotifers is ancient, ongoing and more frequent in species from desiccating habitats

Background

Although prevalent in prokaryotes, horizontal gene transfer (HGT) is rarer in multicellular eukaryotes. Bdelloid rotifers are microscopic animals that contain a higher proportion of horizontally transferred, non-metazoan genes in their genomes than typical of animals. It has been hypothesized that bdelloids incorporate foreign DNA when they repair their chromosomes following double-strand breaks caused by desiccation. HGT might thereby contribute to species divergence and adaptation, as in prokaryotes. If so, we expect that species should differ in their complement of foreign genes, rather than sharing the same set of foreign genes inherited from a common ancestor. Furthermore, there should be more foreign genes in species that desiccate more frequently. We tested these hypotheses by surveying HGT in four congeneric species of bdelloids from different habitats: two from permanent aquatic habitats and two from temporary aquatic habitats that desiccate regularly.

Results

Transcriptomes of all four species contain many genes with a closer match to non-metazoan genes than to metazoan genes. Whole genome sequencing of one species confirmed the presence of these foreign genes in the genome. Nearly half of foreign genes are shared between all four species and an outgroup from another family, but many hundreds are unique to particular species, which indicates that HGT is ongoing. Using a dated phylogeny, we estimate an average of 12.8 gains versus 2.0 losses of foreign genes per million years. Consistent with the desiccation hypothesis, the level of HGT is higher in the species that experience regular desiccation events than those that do not. However, HGT still contributed hundreds of foreign genes to the species from permanently aquatic habitats. Foreign genes were mainly enzymes with various annotated functions that include catabolism of complex polysaccharides and stress responses. We found evidence of differential loss of ancestral foreign genes previously associated with desiccation protection in the two non-desiccating species.

Conclusions

Nearly half of foreign genes were acquired before the divergence of bdelloid families over 60 Mya. Nonetheless, HGT is ongoing in bdelloids and has contributed to putative functional differences among species. Variation among our study species is consistent with the hypothesis that desiccating habitats promote HGT.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0202-9) contains supplementary material, which is available to authorized users.

The marbled crayfish as a paradigm for saltational speciation by autopolyploidy and parthenogenesis in animals

The parthenogenetic all-female marbled crayfish is a novel research model and potent invader of freshwater ecosystems. It is a triploid descendant of the sexually reproducing slough crayfish, Procambarus fallax, but its taxonomic status has remained unsettled. By cross-breeding experiments and parentage analysis we show here that marbled crayfish and P. fallax are reproductively separated. Both crayfish copulate readily, suggesting that the reproductive barrier is set at the cytogenetic rather than the behavioural level. Analysis of complete mitochondrial genomes of marbled crayfish from laboratory lineages and wild populations demonstrates genetic identity and indicates a single origin. Flow cytometric comparison of DNA contents of haemocytes and analysis of nuclear microsatellite loci confirm triploidy and suggest autopolyploidisation as its cause. Global DNA methylation is significantly reduced in marbled crayfish implying the involvement of molecular epigenetic mechanisms in its origination. Morphologically, both crayfish are very similar but growth and fecundity are considerably larger in marbled crayfish, making it a different animal with superior fitness. These data and the high probability of a divergent future evolution of the marbled crayfish and P. fallax clusters suggest that marbled crayfish should be considered as an independent asexual species. Our findings also establish the P. fallax-marbled crayfish pair as a novel paradigm for rare chromosomal speciation by autopolyploidy and parthenogenesis in animals and for saltational evolution in general.

Do Species Exist in Asexuals? Theory and Evidence from Bdelloid Rotifers

The possibility for independently evolving entities to form and persist in the absence of sexual recombination in eukaryotes has been questioned; nevertheless, there are organisms that are known to be asexual and that have apparently diversified into multiple species as recognized by taxonomists. These organisms have therefore been identified as an evolutionary paradox. We explore three alternative hypotheses attempting to solve the apparent paradox, focusing on bdelloid rotifers, the most studied group of organisms in which all species are considered asexual: (1) they may have some hidden form of sex; (2) species do not represent biological entities but simply convenient names; and (3) sex may not be a necessary requirement for speciation. We provide ample evidence against the first two hypotheses, reporting several studies supporting (1) bdelloids asexuality from different approaches, and (2) the existence of species from genetics, jaw morphology, ecology, and physiology. Thus, we (3) explore the role of sex in speciation comparing bdelloid and monogonont rotifers, and conclude with some caveats that could still change our understanding of bdelloid species.

Effects of phylogenetic reconstruction method on the robustness of species delimitation using single-locus data

Coalescent-based species delimitation methods combine population genetic and phylogenetic theory to provide an objective means for delineating evolutionarily significant units of diversity. The generalised mixed Yule coalescent (GMYC) and the Poisson tree process (PTP) are methods that use ultrametric (GMYC or PTP) or non-ultrametric (PTP) gene trees as input, intended for use mostly with single-locus data such as DNA barcodes.

Here, we assess how robust the GMYC and PTP are to different phylogenetic reconstruction and branch smoothing methods. We reconstruct over 400 ultrametric trees using up to 30 different combinations of phylogenetic and smoothing methods and perform over 2000 separate species delimitation analyses across 16 empirical data sets. We then assess how variable diversity estimates are, in terms of richness and identity, with respect to species delimitation, phylogenetic and smoothing methods.

The PTP method generally generates diversity estimates that are more robust to different phylogenetic methods. The GMYC is more sensitive, but provides consistent estimates for BEAST trees. The lower consistency of GMYC estimates is likely a result of differences among gene trees introduced by the smoothing step. Unresolved nodes (real anomalies or methodological artefacts) affect both GMYC and PTP estimates, but have a greater effect on GMYC estimates. Branch smoothing is a difficult step and perhaps an underappreciated source of bias that may be widespread among studies of diversity and diversification.

Nevertheless, careful choice of phylogenetic method does produce equivalent PTP and GMYC diversity estimates. We recommend simultaneous use of the PTP model with any model-based gene tree (e.g. RAxML) and GMYC approaches with BEAST trees for obtaining species hypotheses.