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

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

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.

Geographic origin shapes the adaptive divergences of Rotaria rotatoria (Rotifera, Bdelloidea) to thermal stress: Insights from ecology and transcriptomics

Global warming has raised concerns regarding the potential impact on aquatic biosafety and health. To illuminate the adaptive mechanisms of bdelloid rotifers in response to global warming, the ecological and transcriptomic characteristics of two strains (HX and ZJ) of Rotaria rotatoria were investigated at 25°C and 35°C. Our results showed an obvious genetic divergence between the two geographic populations. Thermal stress significantly reduced the average lifespan of R. rotatoria in both strains, but increased the offspring production in the ZJ strain. Furthermore, the expression levels of genes Hsp70 were significantly upregulated in the HX strain, while GSTo1 and Cu/Zn-SOD were on the contrary. In the ZJ strain, the expression levels of genes Hsp70, CAT2, and GSTo1 were upregulated under thermal stress. Conversely, a significant decrease in the expression level of the Mn-SOD gene was observed in the ZJ strain under thermal stress. Transcriptomic profiling analysis revealed a total of 105 and 5288 differentially expressed genes (DEGs) in the HX and ZJ strains under thermal stress, respectively. The PCA results showed clear differences in gene expression pattern between HX and ZJ strains under thermal stress. Interestingly, compared to the HX strain, numerous downregulated DEGs in the ZJ strain were enriched into pathways related to metabolism under thermal stress, suggesting that rotifers from the ZJ strain prioritize resource allocation to reproduction by suppressing costly metabolic processes. This finding is consistent with the life table results. This study provides new insights into the adaptive evolution of aquatic animals in the context of global climate change.

Ionizing radiation responses appear incidental to desiccation responses in the bdelloid rotifer Adineta vaga

BACKGROUND

The remarkable resistance to ionizing radiation found in anhydrobiotic organisms, such as some bacteria, tardigrades, and bdelloid rotifers has been hypothesized to be incidental to their desiccation resistance. Both stresses produce reactive oxygen species and cause damage to DNA and other macromolecules. However, this hypothesis has only been investigated in a few species.

RESULTS

In this study, we analyzed the transcriptomic response of the bdelloid rotifer Adineta vaga to desiccation and to low- (X-rays) and high- (Fe) LET radiation to highlight the molecular and genetic mechanisms triggered by both stresses. We identified numerous genes encoding antioxidants, but also chaperones, that are constitutively highly expressed, which may contribute to the protection of proteins against oxidative stress during desiccation and ionizing radiation. We also detected a transcriptomic response common to desiccation and ionizing radiation with the over-expression of genes mainly involved in DNA repair and protein modifications but also genes with unknown functions that were bdelloid-specific. A distinct transcriptomic response specific to rehydration was also found, with the over-expression of genes mainly encoding Late Embryogenesis Abundant proteins, specific heat shock proteins, and glucose repressive proteins.

CONCLUSIONS

These results suggest that the extreme resistance of bdelloid rotifers to radiation might indeed be a consequence of their capacity to resist complete desiccation. This study paves the way to functional genetic experiments on A. vaga targeting promising candidate proteins playing central roles in radiation and desiccation resistance.

Horizontal acquisition of a DNA ligase improves DNA damage tolerance in eukaryotes

Bdelloid rotifers are part of the restricted circle of multicellular animals that can withstand a wide range of genotoxic stresses at any stage of their life cycle. In this study, bdelloid rotifer Adineta vaga is used as a model to decipher the molecular basis of their extreme tolerance. Proteomic analysis shows that a specific DNA ligase, different from those usually involved in DNA repair in eukaryotes, is strongly over-represented upon ionizing radiation. A phylogenetic analysis reveals its orthology to prokaryotic DNA ligase E, and its horizontal acquisition by bdelloid rotifers and plausibly other eukaryotes. The fungus Mortierella verticillata, having a single copy of this DNA Ligase E homolog, also exhibits an increased radiation tolerance with an over-expression of this DNA ligase E following X-ray exposure. We also provide evidence that A. vaga ligase E is a major contributor of DNA breaks ligation activity, which is a common step of all important DNA repair pathways. Consistently, its heterologous expression in human cell lines significantly improves their radio-tolerance. Overall, this study highlights the potential of horizontal gene transfers in eukaryotes, and their contribution to the adaptation to extreme conditions.

Parthenogenesis in dipterans: a genetic perspective

Parthenogenesis has been documented in almost every phylum of animals, and yet this phenomenon is largely understudied. It has particular importance in dipterans since some parthenogenetic species are also disease vectors and agricultural pests. Here, we present a catalogue of parthenogenetic dipterans, although it is likely that many more remain to be identified, and we discuss how their developmental biology and interactions with diverse environments may be linked to different types of parthenogenetic reproduction. We discuss how the advances in genetics and genomics have identified chromosomal loci associated with parthenogenesis. In particular, a polygenic cause of facultative parthenogenesis has been uncovered in Drosophila mercatorum, allowing the corresponding genetic variants to be tested for their ability to promote parthenogenesis in another species, Drosophila melanogaster. This study probably identifies just one of many routes that could be followed in the evolution of parthenogenesis. We attempt to account for why the phenomenon has evolved so many times in the dipteran order and why facultative parthenogenesis appears particularly prevalent. We also discuss the significance of coarse genomic changes, including non-disjunction, aneuploidy, and polyploidy and how, together with changes to specific genes, these might relate to both facultative and obligate parthenogenesis in dipterans and other parthenogenetic animals.

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.

Deciphering the Biological Enigma-Genomic Evolution Underlying Anhydrobiosis in the Phylum Tardigrada and the Chironomid Polypedilum vanderplanki

Anhydrobiosis, an ametabolic dehydrated state triggered by water loss, is observed in several invertebrate lineages. Anhydrobiotes revive when rehydrated, and seem not to suffer the ultimately lethal cell damage that results from severe loss of water in other organisms. Here, we review the biochemical and genomic evidence that has revealed the protectant molecules, repair systems, and maintenance pathways associated with anhydrobiosis. We then introduce two lineages in which anhydrobiosis has evolved independently: Tardigrada, where anhydrobiosis characterizes many species within the phylum, and the genus Polypedilum, where anhydrobiosis occurs in only two species. Finally, we discuss the complexity of the evolution of anhydrobiosis within invertebrates based on current knowledge, and propose perspectives to enhance the understanding of anhydrobiosis.

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.

Strategies of tolerance reflected in two North American maple genomes

The first chromosome‐scale assemblies for North American members of the Acer genus, sugar maple (Acer saccharum) and boxelder (Acer negundo), as well as transcriptomic evaluation of the abiotic stress response in A. saccharum are reported. This integrated study describes in‐depth aspects contributing to each species' approach to tolerance and applies current knowledge in many areas of plant genome biology with Acer physiology to help convey the genomic complexities underlying tolerance in broadleaf tree species.

Bacterial N4-methylcytosine as an epigenetic mark in eukaryotic DNA

DNA modifications are used to regulate gene expression and defend against invading genetic elements. In eukaryotes, modifications predominantly involve C5-methylcytosine (5mC) and occasionally N6-methyladenine (6mA), while bacteria frequently use N4-methylcytosine (4mC) in addition to 5mC and 6mA. Here we report that 4mC can serve as an epigenetic mark in eukaryotes. Bdelloid rotifers, tiny freshwater invertebrates with transposon-poor genomes rich in foreign genes, lack canonical eukaryotic C5-methyltransferases for 5mC addition, but encode an amino-methyltransferase, N4CMT, captured from bacteria >60 Mya. N4CMT deposits 4mC at active transposons and certain tandem repeats, and fusion to a chromodomain shapes its “histone-read-DNA-write” architecture recognizing silent chromatin marks. Furthermore, amplification of SETDB1 H3K9me3 histone methyltransferases yields variants preferentially binding 4mC-DNA, suggesting “DNA-read-histone-write” partnership to maintain chromatin-based silencing. Our results show how non-native DNA methyl groups can reshape epigenetic systems to silence transposons and demonstrate the potential of horizontal gene transfer to drive regulatory innovation in eukaryotes.

Eukaryotic DNA can be methylated as 5-methylcytosine and N6-methyladenine, but whether other forms of DNA methylation occur has been controversial. Here the authors show that a bacterial DNA methyltransferase was acquired >60 Mya in bdelloid rotifers that catalyzes N4-methylcytosine addition and is involved in suppression of transposon proliferation.

Convergent consequences of parthenogenesis on stick insect genomes

The shift from sexual reproduction to parthenogenesis has occurred repeatedly in animals, but how the loss of sex affects genome evolution remains poorly understood. We generated reference genomes for five independently evolved parthenogenetic species in the stick insect genus Timema and their closest sexual relatives. Using these references and population genomic data, we show that parthenogenesis results in an extreme reduction of heterozygosity and often leads to genetically uniform populations. We also find evidence for less effective positive selection in parthenogenetic species, suggesting that sex is ubiquitous in natural populations because it facilitates fast rates of adaptation. Parthenogenetic species did not show increased transposable element (TE) accumulation, likely because there is little TE activity in the genus. By using replicated sexual-parthenogenetic comparisons, our study reveals how the absence of sex affects genome evolution in natural populations, providing empirical support for the negative consequences of parthenogenesis as predicted by theory.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Transposable elements promote the evolution of genome streamlining

Eukaryotes and prokaryotes have distinct genome architectures, with marked differences in genome size, the ratio of coding/non-coding DNA, and the abundance of transposable elements (TEs). As TEs replicate independently of their hosts, the proliferation of TEs is thought to have driven genome expansion in eukaryotes. However, prokaryotes also have TEs in intergenic spaces, so why do prokaryotes have small, streamlined genomes? Using an in silico model describing the genomes of single-celled asexual organisms that coevolve with TEs, we show that TEs acquired from the environment by horizontal gene transfer can promote the evolution of genome streamlining. The process depends on local interactions and is underpinned by rock-paper-scissors dynamics in which populations of cells with streamlined genomes beat TEs, which beat non-streamlined genomes, which beat streamlined genomes, in continuous and repeating cycles. Streamlining is maladaptive to individual cells, but improves lineage viability by hindering the proliferation of TEs. Streamlining does not evolve in sexually reproducing populations because recombination partially frees TEs from the deleterious effects they cause. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Genomic Signature of Sexual Reproduction in the Bdelloid Rotifer Macrotrachella quadricornifera

Bdelloid rotifers, common freshwater invertebrates of ancient origin and worldwide distribution have long been thought to be entirely asexual, being the principal exception to the view that in eukaryotes the loss of sex leads to early extinction. That bdelloids are facultatively sexual is shown by a study of allele sharing within a group of closely related bdelloids of the species Macrotrachella quadricornifera, supporting the view that sexual reproduction is essential for long-term success in all eukaryotes.

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.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host

Background

Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensulato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea.

Results

Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria).

Conclusion

Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07857-y.

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.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Evolutionary dynamics of transposable elements in bdelloid rotifers

Transposable elements (TEs) are selfish genomic parasites whose ability to spread autonomously is facilitated by sexual reproduction in their hosts. If hosts become obligately asexual, TE frequencies and dynamics are predicted to change dramatically, but the long-term outcome is unclear. Here, we test current theory using whole-genome sequence data from eight species of bdelloid rotifers, a class of invertebrates in which males are thus far unknown. Contrary to expectations, we find a variety of active TEs in bdelloid genomes, at an overall frequency within the range seen in sexual species. We find no evidence that TEs are spread by cryptic recombination or restrained by unusual DNA repair mechanisms. Instead, we find that that TE content evolves relatively slowly in bdelloids and that gene families involved in RNAi-mediated TE suppression have undergone significant expansion, which might mitigate the deleterious effects of active TEs and compensate for the consequences of long-term asexuality.

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.

Comparative Analysis of ROS Network Genes in Extremophile Eukaryotes

The reactive oxygen species (ROS) gene network, consisting of both ROS-generating and detoxifying enzymes, adjusts ROS levels in response to various stimuli. We performed a cross-kingdom comparison of ROS gene networks to investigate how they have evolved across all Eukaryotes, including protists, fungi, plants and animals. We included the genomes of 16 extremotolerant Eukaryotes to gain insight into ROS gene evolution in organisms that experience extreme stress conditions. Our analysis focused on ROS genes found in all Eukaryotes (such as catalases, superoxide dismutases, glutathione reductases, peroxidases and glutathione peroxidase/peroxiredoxins) as well as those specific to certain groups, such as ascorbate peroxidases, dehydroascorbate/monodehydroascorbate reductases in plants and other photosynthetic organisms. ROS-producing NADPH oxidases (NOX) were found in most multicellular organisms, although several NOX-like genes were identified in unicellular or filamentous species. However, despite the extreme conditions experienced by extremophile species, we found no evidence for expansion of ROS-related gene families in these species compared to other Eukaryotes. Tardigrades and rotifers do show ROS gene expansions that could be related to their extreme lifestyles, although a high rate of lineage-specific horizontal gene transfer events, coupled with recent tetraploidy in rotifers, could explain this observation. This suggests that the basal Eukaryotic ROS scavenging systems are sufficient to maintain ROS homeostasis even under the most extreme conditions.

Improved Reference Genome for Cyclotella cryptica CCMP332, a Model for Cell Wall Morphogenesis, Salinity Adaptation, and Lipid Production in Diatoms (Bacillariophyta)

The diatom, Cyclotella cryptica, is a well-established model species for physiological studies and biotechnology applications of diatoms. To further facilitate its use as a model diatom, we report an improved reference genome assembly and annotation for C. cryptica strain CCMP332. We used a combination of long- and short-read sequencing to assemble a high-quality and contaminant-free genome. The genome is 171 Mb in size and consists of 662 scaffolds with a scaffold N50 of 494 kb. This represents a 176-fold decrease in scaffold number and 41-fold increase in scaffold N50 compared to the previous assembly. The genome contains 21,250 predicted genes, 75% of which were assigned putative functions. Repetitive DNA comprises 59% of the genome, and an improved classification of repetitive elements indicated that a historically steady accumulation of transposable elements has contributed to the relatively large size of the C. cryptica genome. The high-quality C. cryptica genome will serve as a valuable reference for ecological, genetic, and biotechnology studies of diatoms.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

Amino acid sequence associated with bacteriophage recombination site helps to reveal genes potentially acquired through horizontal gene transfer

Background

Horizontal gene transfer, i.e. the acquisition of genetic material from nonparent organism, is considered an important force driving species evolution. Many cases of horizontal gene transfer from prokaryotes to eukaryotes have been registered, but no transfer mechanism has been deciphered so far, although viruses were proposed as possible vectors in several studies. In agreement with this idea, in our previous study we discovered that in two eukaryotic proteins bacteriophage recombination site (AttP) was adjacent to the regions originating via horizontal gene transfer. In one of those cases AttP site was present inside the introns of cysteine-rich repeats. In the present study we aimed to apply computational tools for finding multiple horizontal gene transfer events in large genome databases. For that purpose we used a sequence of cysteine-rich repeats to identify genes potentially acquired through horizontal transfer.

Results

HMMER remote similarity search significantly detected 382 proteins containing cysteine-rich repeats. All of them, except 8 sequences, belong to eukaryotes. In 124 proteins the presence of conserved structural domains was predicted. In spite of the fact that cysteine-rich repeats are found almost exclusively in eukaryotic proteins, many predicted domains are most common for prokaryotes or bacteriophages. Ninety-eight proteins out of 124 contain typical prokaryotic domains. In those cases proteins were considered as potentially originating via horizontal transfer. In addition, HHblits search revealed that two domains of the same fungal protein, Glycoside hydrolase and Peptidase M15, have high similarity with proteins of two different prokaryotic species, hinting at independent horizontal gene transfer events.

Conclusions

Cysteine-rich repeats in eukaryotic proteins are usually accompanied by conserved domains typical for prokaryotes or bacteriophages. These proteins, containing both cysteine-rich repeats, and characteristic prokaryotic domains, might represent multiple independent horizontal gene transfer events from prokaryotes to eukaryotes. We believe that the presence of bacteriophage recombination site inside cysteine-rich repeat coding sequence may facilitate horizontal genes transfer. Thus computational approach, described in the present study, can help finding multiple sequences originated from horizontal transfer in eukaryotic genomes.

Desiccation does not drastically increase the accessibility of exogenous DNA to nuclear genomes: evidence from the frequency of endosymbiotic DNA transfer

Background

Although horizontal gene transfer (HGT) is a widely accepted force in the evolution of prokaryotic genomes, its role in the evolution of eukaryotic genomes remains hotly debated. Some bdelloid rotifers that are resistant to extreme desiccation and radiation undergo a very high level of HGT, whereas in another desiccation-resistant invertebrate, the tardigrade, the pattern does not exist. Overall, the DNA double-strand breaks (DSBs) induced by prolonged desiccation have been postulated to open a gateway to the nuclear genome for exogenous DNA integration and thus to facilitate the HGT process, thereby enhancing the rate of endosymbiotic DNA transfer (EDT).

Results

We first surveyed the abundance of nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs) in five eukaryotes that are highly resistant to desiccation: the bdelloid rotifers Adineta vaga and Adineta ricciae, the tardigrade Ramazzottius varieornatus, and the resurrection plants Dorcoceras hygrometricum and Selaginella tamariscina. Excessive NUMTs or NUPTs were not detected. Furthermore, we compared 24 groups of desiccation-tolerant organisms with their relatively less desiccation-tolerant relatives but did not find a significant difference in NUMT/NUPT contents.

Conclusions

Desiccation may induce DSBs, but it is unlikely to dramatically increase the frequency of exogenous sequence integration in most eukaryotes. The capture of exogenous DNA sequences is possible only when DSBs are repaired through a subtype of non-homologous end joining, named alternative end joining (alt-EJ). Due to the deleterious effects of the resulting insertion mutations, alt-EJ is less frequently initiated than other mechanisms.

The genome of pest Rhynchophorus ferrugineus reveals gene families important at the plant-beetle interface

The red palm weevil, Rhynchophorus ferrugineus, infests palm plantations, leading to large financial losses and soil erosion. Pest-host interactions are poorly understood in R. ferrugineus, but the analysis of genetic diversity and pest origins will help advance efforts to eradicate this pest. We sequenced the genome of R. ferrugineus using a combination of paired-end Illumina sequencing (150 bp), Oxford Nanopore long reads, 10X Genomics and synteny analysis to produce an assembly with a scaffold N50 of ~60 Mb. Structural variations showed duplication of detoxifying and insecticide resistance genes (e.g., glutathione S-transferase, P450, Rdl). Furthermore, the evolution of gene families identified those under positive selection including one glycosyl hydrolase (GH16) gene family, which appears to result from horizontal gene transfer. This genome will be a valuable resource to understand insect evolution and behavior and to allow the genetic modification of key genes that will help control this pest.

HGT in the human and skin commensal Malassezia: A bacterially derived flavohemoglobin is required for NO resistance and host interaction

The skin of humans and animals is colonized by commensal and pathogenic fungi and bacteria that share this ecological niche and have established microbial interactions. Malassezia are the most abundant fungal skin inhabitant of warm-blooded animals and have been implicated in skin diseases and systemic disorders, including Crohn's disease and pancreatic cancer. Flavohemoglobin is a key enzyme involved in microbial nitrosative stress resistance and nitric oxide degradation. Comparative genomics and phylogenetic analyses within the Malassezia genus revealed that flavohemoglobin-encoding genes were acquired through independent horizontal gene transfer events from different donor bacteria that are part of the mammalian microbiome. Through targeted gene deletion and functional complementation in Malassezia sympodialis, we demonstrated that bacterially derived flavohemoglobins are cytoplasmic proteins required for nitric oxide detoxification and nitrosative stress resistance under aerobic conditions. RNA-sequencing analysis revealed that endogenous accumulation of nitric oxide resulted in up-regulation of genes involved in stress response and down-regulation of the MalaS7 allergen-encoding genes. Solution of the high-resolution X-ray crystal structure of Malassezia flavohemoglobin revealed features conserved with both bacterial and fungal flavohemoglobins. In vivo pathogenesis is independent of Malassezia flavohemoglobin. Lastly, we identified an additional 30 genus- and species-specific horizontal gene transfer candidates that might have contributed to the evolution of this genus as the most common inhabitants of animal skin.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

The genome of the marine rotifer Brachionus koreanus sheds light on the antioxidative defense system in response to 2-ethyl-phenanthrene and piperonyl butoxide

BRACHIONUS: spp. (Rotifera: Monogononta) have been introduced as ecotoxicological model-organisms that are widely distributed in aquatic environments. Among the Brachionus spp., the monogonont rotifer Brachionus koreanus has been widely used for ecology, ecotoxicology, and evolution, thus, providing the whole genome data of B. koreanus is important for further understandings of in-depth molecular mechanisms. In this study, the completed assembly and characterization of the B. koreanus genome resulted in a total length of 85.7 Mb with 14,975 annotated genes. The final number of scaffolds was 567 with an N50 value and a GC content of 1.86 Mb and 24.35 %, respectively. Based on the fully constructed genome database, a total of 24 CYPs, 23 GSTs, two SODs, and a single CAT genes were identified and analyzed antioxidant activities (CAT, SOD, and GST), and transcriptional regulation of the entire CYPs, GSTs, SODs, and CAT in response to 2-ethyl-phenanthrene (2-ethyl-PHE) and piperonyl butoxide (PBO), to demonstrate the usefulness of the whole genome library of B. koreanus in response xenobiotic-induced oxidative stress. The assembled B. koreanus genome will provide a better understanding on the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification processes in the rotifer B. koreanus.

GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes

An important assessment prior to genome assembly and related analyses is genome profiling, where the k-mer frequencies within raw sequencing reads are analyzed to estimate major genome characteristics such as size, heterozygosity, and repetitiveness. Here we introduce GenomeScope 2.0 (https://github.com/tbenavi1/genomescope2.0), which applies combinatorial theory to establish a detailed mathematical model of how k-mer frequencies are distributed in heterozygous and polyploid genomes. We describe and evaluate a practical implementation of the polyploid-aware mixture model that quickly and accurately infers genome properties across thousands of simulated and several real datasets spanning a broad range of complexity. We also present a method called Smudgeplot (https://github.com/KamilSJaron/smudgeplot) to visualize and estimate the ploidy and genome structure of a genome by analyzing heterozygous k-mer pairs. We successfully apply the approach to systems of known variable ploidy levels in the Meloidogyne genus and the extreme case of octoploid Fragaria × ananassa.

Prevalence and Implications of Contamination in Public Genomic Resources: A Case Study of 43 Reference Arthropod Assemblies

Thanks to huge advances in sequencing technologies, genomic resources are increasingly being generated and shared by the scientific community. The quality of such public resources are therefore of critical importance. Errors due to contamination are particularly worrying; they are widespread, propagate across databases, and can compromise downstream analyses, especially the detection of horizontally-transferred sequences. However we still lack consistent and comprehensive assessments of contamination prevalence in public genomic data. Here we applied a standardized procedure for foreign sequence annotation to 43 published arthropod genomes from the widely used Ensembl Metazoa database. This method combines information on sequence similarity and synteny to identify contaminant and putative horizontally-transferred sequences in any genome assembly, provided that an adequate reference database is available. We uncovered considerable heterogeneity in quality among arthropod assemblies, some being devoid of contaminant sequences, whereas others included hundreds of contaminant genes. Contaminants far outnumbered horizontally-transferred genes and were a major confounder of their detection, quantification and analysis. We strongly recommend that automated standardized decontamination procedures be systematically embedded into the submission process to genomic databases.

Signatures of the Evolution of Parthenogenesis and Cryptobiosis in the Genomes of Panagrolaimid Nematodes

Most animal species reproduce sexually and fully parthenogenetic lineages are usually short lived in evolution. Still, parthenogenesis may be advantageous as it avoids the cost of sex and permits colonization by single individuals. Panagrolaimid nematodes have colonized environments ranging from arid deserts to Arctic and Antarctic biomes. Many are obligatory meiotic parthenogens, and most have cryptobiotic abilities, being able to survive repeated cycles of complete desiccation and freezing. To identify systems that may contribute to these striking abilities, we sequenced and compared the genomes and transcriptomes of parthenogenetic and outcrossing panagrolaimid species, including cryptobionts and non-cryptobionts. The parthenogens are triploids, most likely originating through hybridization. Adaptation to cryptobiosis shaped the genomes of panagrolaimid nematodes and is associated with the expansion of gene families and signatures of selection on genes involved in cryptobiosis. All panagrolaimids have acquired genes through horizontal gene transfer, some of which are likely to contribute to cryptobiosis.

Genome size and lifestyle in gnesiotrochan rotifers

Gnesiotrochan rotifers display a variety of life styles ranging from taxa with free-swimming larval and sessile adult stages to those with motile adult stages and colonial habits. Several explanations for the C- value enigma posits that genome size is correlated with lifestyle. To investigate this, 13 gnesiotrochan species representing nine genera were measured by flow cytometry. Genome sizes (1C) within Gnesiotrocha ranged from 0.05 pg (Hexarthra mira and Hexarthra fennica) to 0.25 pg (Sinantherina ariprepes). Genome sizes varied within genera and species; e.g., H. fennica (El Huérfano, Mexico) was estimated to be 15% larger than that of H. mira and H. fennica (Keystone Wetland, TX, USA). Gnesiotrochan genome sizes are similar to those reported within Ploima, which range from 0.06 pg (Brachionus rotundiformis, B. dimidiatus) to 0.46 pg (B. asplanchnoidis). Within Gnesiotrocha, genome size was found to be significantly smaller in sessile versus motile species as well as in solitary versus colonial species. To account for phylogenetic background, Linear Mixed Models with hierarchical taxonomic ranks showed that there is a taxonomic component underlying genome size. This study provides the first estimates of genome size within the superorder, providing a baseline for genomic and evolutionary studies within the group.

The genome of the marine monogonont rotifer Brachionus plicatilis: Genome-wide expression profiles of 28 cytochrome P450 genes in response to chlorpyrifos and 2-ethyl-phenanthrene

Brachionus spp. (Rotifera: Monogononta) are globally distributed in aquatic environments and play important roles in the aquatic ecosystem. The marine monogonont rotifer Brachionus plicatilis is considered a suitable model organism for ecology, evolution, and ecotoxicology. In this study, we assembled and characterized the B. plicatilis genome. The total length of the assembled genome was 106.9 Mb and the number of final scaffolds was 716 with an N50 value of 1.15 Mb and a GC content of 26.75%. A total of 20,154 genes were annotated after manual curation. To demonstrate the use of whole genome data, we targeted one of the main detoxifying enzyme of phase I detoxification system and identified in a total of 28 cytochrome P450 s (CYPs). Based on the phylogenetic analysis using the maximum likelihood, 28 B. plicatilis-CYPs were apparently separated into five different clans, namely, 2, 3, 4, mitochondrial (MT), and 46 clans. To better understand the CYPs-mediated xenobiotic detoxification, we measured the mRNA expression levels of 28 B. plicatilis CYPs in response to chlorpyrifos and 2-ethyl-phenanthrene. Most B. plicatilis CYPs were significantly modulated (P < 0.05) in response to chlorpyrifos and 2-ethyl-phenanthrene. In addition, xenobiotic-sensing nuclear receptor (XNR) response element sequences were identified in the 5 kb upstream of promoter regions of 28 CYPs from the genome of B. plicatilis, indicating that these XNR can be associated with detoxification of xenobiotics. Overall, the assembled B. plicatilis genome presented here will be a useful resource for a better understanding the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification in this species.

All Eukaryotes Are Sexual, unless Proven Otherwise: Many So-Called Asexuals Present Meiotic Machinery and Might Be Able to Have Sex

Here a wide distribution of meiotic machinery is shown, indicating the occurrence of sexual processes in all major eukaryotic groups, without exceptions, including the putative "asexuals." Meiotic machinery has evolved from archaeal DNA repair machinery by means of ancestral gene duplications. Sex is very conserved and widespread in eukaryotes, even though its evolutionary importance is still a matter of debate. The main processes in sex are plasmogamy, followed by karyogamy and meiosis. Meiosis is fundamentally a chromosomal process, which implies recombination and ploidy reduction. Several eukaryotic lineages are proposed to be asexual because their sexual processes are never observed, but presumed asexuality correlates with lack of study. The authors stress the complete lack of meiotic proteins in nucleomorphs and their almost complete loss in the fungus Malassezia. Inversely, complete sets of meiotic proteins are present in fungal groups Glomeromycotina, Trichophyton, and Cryptococcus. Endosymbiont Perkinsela and endoparasitic Microsporidia also present meiotic proteins.

Evolutionary dynamics of origin and loss in the deep history of phospholipase D toxin genes

Background

Venom-expressed sphingomyelinase D/phospholipase D (SMase D/PLD) enzymes evolved from the ubiquitous glycerophosphoryl diester phosphodiesterases (GDPD). Expression of GDPD-like SMaseD/PLD toxins in both arachnids and bacteria has inspired consideration of the relative contributions of lateral gene transfer and convergent recruitment in the evolutionary history of this lineage. Previous work recognized two distinct lineages, a SicTox-like (ST-like) clade including the arachnid toxins, and an Actinobacterial-toxin like (AT-like) clade including the bacterial toxins and numerous fungal homologs.

Results

Here we expand taxon sampling by homology detection to discover new GDPD-like SMase D/PLD homologs. The ST-like clade now includes homologs in a wider variety of arthropods along with a sister group in Cnidaria; the AT-like clade now includes additional fungal phyla and proteobacterial homologs; and we report a third clade expressed in diverse aquatic metazoan taxa, a few single-celled eukaryotes, and a few aquatic proteobacteria. GDPD-like SMaseD/PLDs have an ancient presence in chelicerates within the ST-like family and ctenophores within the Aquatic family. A rooted phylogenetic tree shows that the three clades derived from a basal paraphyletic group of proteobacterial GDPD-like SMase D/PLDs, some of which are on mobile genetic elements. GDPD-like SMase D/PLDs share a signature C-terminal motif and a shortened βα1 loop, features that distinguish them from GDPDs. The three major clades also have active site loop signatures that distinguish them from GDPDs and from each other. Analysis of molecular phylogenies with respect to organismal relationships reveals a dynamic evolutionary history including both lateral gene transfer and gene duplication/loss.

Conclusions

The GDPD-like SMaseD/PLD enzymes derive from a single ancient ancestor, likely proteobacterial, and radiated into diverse organismal lineages at least in part through lateral gene transfer.

Electronic supplementary material

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

Cross-Contamination Explains "Inter and Intraspecific Horizontal Genetic Transfers" between Asexual Bdelloid Rotifers

A few metazoan lineages are thought to have persisted for millions of years without sexual reproduction. If so, they would offer important clues to the evolutionary paradox of sex itself [1, 2]. Most "ancient asexuals" are subject to ongoing doubt because extant populations continue to invest in males [3-9]. However, males are famously unknown in bdelloid rotifers, a class of microscopic invertebrates comprising hundreds of species [10-12]. Bdelloid genomes have acquired an unusually high proportion of genes from non-metazoans via horizontal transfer [13-17]. This well-substantiated finding has invited speculation [13] that homologous horizontal transfer between bdelloid individuals also may occur, perhaps even "replacing" sex [14]. In 2016, Current Biology published an article claiming to supply evidence for this idea. Debortoli et al. [18] sampled rotifers from natural populations and sequenced one mitochondrial and four nuclear loci. Species assignments were incongruent among loci for several samples, which was interpreted as evidence of "interspecific horizontal genetic transfers." Here, we use sequencing chromatograms supplied by the authors to demonstrate that samples treated as individuals actually contained two or more highly divergent mitochondrial and ribosomal sequences, revealing cross-contamination with DNA from multiple animals of different species. Other chromatograms indicate contamination with DNA from conspecific animals, explaining genetic and genomic evidence for "intraspecific horizontal exchanges" reported in the same study. Given the clear evidence of contamination, the data and findings of Debortoli et al. [18] provide no reliable support for their conclusions that DNA is transferred horizontally between or within bdelloid species.

Full disclosure: Genome assembly is still hard

Two recent papers highlight the fascinating comparative genomics of anhydrobiosis, the ability to withstand complete desiccation, in bdelloid rotifers and tardigrades. However, both groups had to openly deal with the significant difficulties of generating and interpreting short-read draft assemblies-especially challenging in microscopic species with high sequence polymorphism. These exemplars demonstrate the need to go beyond single draft-quality reference genomes to high-quality multiple species comparative genomics if we are to fully capture the value of genomics.