Genetic Drift, Not Life History or RNAi, Determine Long-Term Evolution of Transposable Elements

The curious case of the disappearing piRNAs

Small non-coding RNAs are key regulators of gene expression across eukaryotes. Piwi-interacting small RNAs (piRNAs) are a specific type of small non-coding RNAs, conserved across animals, which are best known as regulators of genome stability through their ability to target transposable elements for silencing. Despite the near ubiquitous presence of piRNAs in animal lineages, there are some examples where the piRNA pathway has been lost completely, most dramatically in nematodes where loss has occurred in at least four independent lineages. In this perspective I will provide an evaluation of the presence of piRNAs across animals, explaining how it is known that piRNAs are missing from certain organisms. I will then consider possible explanations for why the piRNA pathway might have been lost and evaluate the evidence in favor of each possible mechanism. While it is still impossible to provide definitive answers, these theories will prompt further investigations into why such a highly conserved pathway can nevertheless become dispensable in certain lineages. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Quantifying the role of genome size and repeat content in adaptive variation and the architecture of flowering time in Amaranthus tuberculatus

Genome size variation, largely driven by repeat content, is poorly understood within and among populations, limiting our understanding of its significance for adaptation. Here we characterize intraspecific variation in genome size and repeat content across 186 individuals of Amaranthus tuberculatus, a ubiquitous native weed that shows flowering time adaptation to climate across its range and in response to agriculture. Sequence-based genome size estimates vary by up to 20% across individuals, consistent with the considerable variability in the abundance of transposable elements, unknown repeats, and rDNAs across individuals. The additive effect of this variation has important phenotypic consequences-individuals with more repeats, and thus larger genomes, show slower flowering times and growth rates. However, compared to newly-characterized gene copy number and polygenic nucleotide changes underlying variation in flowering time, we show that genome size is a marginal contributor. Differences in flowering time are reflected by genome size variation across sexes and marginally, habitats, while polygenic variation and a gene copy number variant within the ATP synthesis pathway show consistently stronger environmental clines than genome size. Repeat content nonetheless shows non-neutral distributions across the genome, and across latitudinal and environmental gradients, demonstrating the numerous governing processes that in turn influence quantitative genetic variation for phenotypes key to plant adaptation.

Single-worm long-read sequencing reveals genome diversity in free-living nematodes

Obtaining sufficient genetic material from a limited biological source is currently the primary operational bottleneck in studies investigating biodiversity and genome evolution. In this study, we employed multiple displacement amplification (MDA) and Smartseq2 to amplify nanograms of genomic DNA and mRNA, respectively, from individual Caenorhabditis elegans. Although reduced genome coverage was observed in repetitive regions, we produced assemblies covering 98% of the reference genome using long-read sequences generated with Oxford Nanopore Technologies (ONT). Annotation with the sequenced transcriptome coupled with the available assembly revealed that gene predictions were more accurate, complete and contained far fewer false positives than de novo transcriptome assembly approaches. We sampled and sequenced the genomes and transcriptomes of 13 nematodes from early-branching species in Chromadoria, Dorylaimia and Enoplia. The basal Chromadoria and Enoplia species had larger genome sizes, ranging from 136.6 to 738.8 Mb, compared with those in the other clades. Nine mitogenomes were fully assembled, and displayed a complete lack of synteny to other species. Phylogenomic analyses based on the new annotations revealed strong support for Enoplia as sister to the rest of Nematoda. Our result demonstrates the robustness of MDA in combination with ONT, paving the way for the study of genome diversity in the phylum Nematoda and beyond.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

Species-specific chromatin landscape determines how transposable elements shape genome evolution

Transposable elements (TEs) are selfish genetic parasites that increase their copy number at the expense of host fitness. The ‘success’, or genome-wide abundance, of TEs differs widely between species. Deciphering the causes for this large variety in TE abundance has remained a central question in evolutionary genomics. We previously proposed that species-specific TE abundance could be driven by the inadvertent consequences of host-direct epigenetic silencing of TEs—the spreading of repressive epigenetic marks from silenced TEs into adjacent sequences. Here, we compared this TE-mediated local enrichment of repressive marks, or ‘the epigenetic effect of TEs’, in six species in the Drosophila melanogaster subgroup to dissect step-by-step the role of such effect in determining genomic TE abundance. We found that TE-mediated local enrichment of repressive marks is prevalent and substantially varies across and even within species. While this TE-mediated effect alters the epigenetic states of adjacent genes, we surprisingly discovered that the transcription of neighboring genes could reciprocally impact this spreading. Importantly, our multi-species analysis provides the power and appropriate phylogenetic resolution to connect species-specific host chromatin regulation, TE-mediated epigenetic effects, the strength of natural selection against TEs, and genomic TE abundance unique to individual species. Our findings point toward the importance of host chromatin landscapes in shaping genome evolution through the epigenetic effects of a selfish genetic parasite.

All the instructions required for life are encoded in the set of DNA present in a cell. It therefore seems natural to think that every bit of this genetic information should serve the organism. And yet most species carry parasitic ‘transposable’ sequences, or transposons, whose only purpose is to multiply and insert themselves at other positions in the genome.

It is possible for cells to suppress these selfish elements. Chemical marks can be deposited onto the DNA to temporarily ‘silence’ transposons and prevent them from being able to move and replicate. However, this sometimes comes at a cost: the repressive chemical modifications can spread to nearby genes that are essential for the organism and perturb their function.

Strangely, the prevalence of transposons varies widely across the tree of life. These sequences form the majority of the genome of certain species – in fact, they represent about half of the human genetic information. But their abundance is much lower in other organisms, forming a measly 6% of the genome of puffer fish for instance. Even amongst fruit fly species, the prevalence of transposable elements can range between 2% and 25%. What explains such differences?

Huang et al. set out to examine this question through the lens of transposon silencing, systematically comparing how this process impacts nearby regions in six species of fruit flies. This revealed variations in the strength of the side effects associated with transposon silencing, resulting in different levels of perturbation on neighbouring genes. A stronger impact was associated with the species having fewer transposons in its genome, suggesting that an evolutionary pressure is at work to keep the abundance of transposons at a low level in these species. Further analyses showed that the genes which determine how silencing marks are distributed may also be responsible for the variations in the impact of transposon silencing. They could therefore be the ones driving differences in the abundance of transposons between species.

Overall, this work sheds light on the complex mechanisms shaping the evolution of genomes, and it may help to better understand how transposons are linked to processes such as aging and cancer.

piRNA-like small RNAs target transposable elements in a Clade IV parasitic nematode

The small RNA (sRNA) pathways identified in the model organism Caenorhabditis elegans are not widely conserved across nematodes. For example, the PIWI pathway and PIWI-interacting RNAs (piRNAs) are involved in regulating and silencing transposable elements (TE) in most animals but have been lost in nematodes outside of the C. elegans group (Clade V), and little is known about how nematodes regulate TEs in the absence of the PIWI pathway. Here, we investigated the role of sRNAs in the Clade IV parasitic nematode Strongyloides ratti by comparing two genetically identical adult stages (the parasitic female and free-living female). We identified putative small-interfering RNAs, microRNAs and tRNA-derived sRNA fragments that are differentially expressed between the two adult stages. Two classes of sRNAs were predicted to regulate TE activity including (i) a parasite-associated class of 21-22 nt long sRNAs with a 5' uridine (21-22Us) and a 5' monophosphate, and (ii) 27 nt long sRNAs with a 5' guanine/adenine (27GAs) and a 5' modification. The 21-22Us show striking resemblance to the 21U PIWI-interacting RNAs found in C. elegans, including an AT rich upstream sequence, overlapping loci and physical clustering in the genome. Overall, we have shown that an alternative class of sRNAs compensate for the loss of piRNAs and regulate TE activity in nematodes outside of Clade V.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Stochastic Effects in Retrotransposon Dynamics Revealed by Modeling under Competition for Cellular Resources

Transposons are genomic elements that can relocate within a host genome using a ‘cut’- or ‘copy-and-paste’ mechanism. They make up a significant part of many genomes, serve as a driving force for genome evolution, and are linked with Mendelian diseases and cancers. Interactions between two specific retrotransposon types, autonomous (e.g., LINE1/L1) and nonautonomous (e.g., Alu), may lead to fluctuations in the number of these transposons in the genome over multiple cell generations. We developed and examined a simple model of retrotransposon dynamics under conditions where transposon replication machinery competed for cellular resources: namely, free ribosomes and available energy (i.e., ATP molecules). Such competition is likely to occur in stress conditions that a malfunctioning cell may experience as a result of a malignant transformation. The modeling revealed that the number of actively replicating LINE1 and Alu elements in a cell decreases with the increasing competition for resources; however, stochastic effects interfere with this simple trend. We stochastically simulated the transposon dynamics in a cell population and showed that the population splits into pools with drastically different transposon behaviors. The early extinction of active Alu elements resulted in a larger number of LINE1 copies occurring in the first pool, as there was no competition between the two types of transposons in this pool. In the other pool, the competition process remained and the number of L1 copies was kept small. As the level of available resources reached a critical value, both types of dynamics demonstrated an increase in noise levels, and both the period and the amplitude of predator–prey oscillations rose in one of the cell pools. We hypothesized that the presented dynamical effects associated with the impact of the competition for cellular resources inflicted on the dynamics of retrotransposable elements could be used as a characteristic feature to assess a cell state, or to control the transposon activity.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

Comparative analyses of the Hymenoscyphus fraxineus and Hymenoscyphus albidus genomes reveals potentially adaptive differences in secondary metabolite and transposable element repertoires

Background

The dieback epidemic decimating common ash (Fraxinus excelsior) in Europe is caused by the invasive fungus Hymenoscyphus fraxineus. In this study we analyzed the genomes of H. fraxineus and H. albidus, its native but, now essentially displaced, non-pathogenic sister species, and compared them with several other members of Helotiales. The focus of the analyses was to identify signals in the genome that may explain the rapid establishment of H. fraxineus and displacement of H. albidus.

Results

The genomes of H. fraxineus and H. albidus showed a high level of synteny and identity. The assembly of H. fraxineus is 13 Mb longer than that of H. albidus’, most of this difference can be attributed to higher dispersed repeat content (i.e. transposable elements [TEs]) in H. fraxineus. In general, TE families in H. fraxineus showed more signals of repeat-induced point mutations (RIP) than in H. albidus, especially in Long-terminal repeat (LTR)/Copia and LTR/Gypsy elements.

Comparing gene family expansions and 1:1 orthologs, relatively few genes show signs of positive selection between species. However, several of those did appeared to be associated with secondary metabolite genes families, including gene families containing two of the genes in the H. fraxineus-specific, hymenosetin biosynthetic gene cluster (BGC).

Conclusion

The genomes of H. fraxineus and H. albidus show a high degree of synteny, and are rich in both TEs and BGCs, but the genomic signatures also indicated that H. albidus may be less well equipped to adapt and maintain its ecological niche in a rapidly changing environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07837-2.

Movements of transposable elements contribute to the genomic plasticity and species diversification in an asexually reproducing nematode pest

Despite reproducing without sexual recombination, Meloidogyne incognita is an adaptive and versatile phytoparasitic nematode. This species displays a global distribution, can parasitize a large range of plants, and can overcome plant resistance in a few generations. The mechanisms underlying this adaptability remain poorly known. At the whole-genome level, only a few single nucleotide variations have been observed across different geographical isolates with distinct ranges of compatible hosts. Exploring other factors possibly involved in genomic plasticity is thus important. Transposable elements (TEs), by their repetitive nature and mobility, can passively and actively impact the genome dynamics. This is particularly expected in polyploid hybrid genomes such as the one of M. incognita. Here, we have annotated the TE content of M. incognita, analyzed the statistical properties of this TE landscape, and used whole-genome pool-seq data to estimate the mobility of these TEs across twelve geographical isolates, presenting variations in ranges of compatible host plants. DNA transposons are more abundant than retrotransposons, and the high similarity of TE copies to their consensus sequences suggests they have been at least recently active. We have identified loci in the genome where the frequencies of presence of a TE showed substantial variations across the different isolates. Overall, variations in TE frequencies across isolates followed their phylogenetic divergence, suggesting TEs participate in the species diversification. Compared with the M. incognita reference genome, we detected isolate and lineage-specific de novo insertion of some TEs, including within genic regions or in the upstream regulatory regions. We validated by PCR the insertion of some of these TEs inside genic regions, confirming TE movements have possible functional impacts. Overall, we show DNA transposons can drive genomic plasticity in M. incognita and their role in genome evolution of other parthenogenetic animal deserves further investigation.

Network-based visualisation reveals new insights into transposable element diversity

Transposable elements (TEs) are widespread across eukaryotic genomes, yet their content varies widely between different species. Factors shaping the diversity of TEs are poorly understood. Understanding the evolution of TEs is difficult because their sequences diversify rapidly and TEs are often transferred through non‐conventional means such as horizontal gene transfer. We developed a method to track TE evolution using network analysis to visualise TE sequence and TE content across different genomes. We illustrate our method by first using a monopartite network to study the sequence evolution of Tc1/mariner elements across focal species. We identify a connection between two subfamilies associated with convergent acquisition of a domain from a protein‐coding gene. Second, we use a bipartite network to study how TE content across species is shaped by epigenetic silencing mechanisms. We show that the presence of Piwi‐interacting RNAs is associated with differences in network topology after controlling for phylogenetic effects. Together, our method demonstrates how a network‐based approach can identify hitherto unknown properties of TE evolution across species.

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 Survey of Transposon Landscapes in the Putative Ancient Asexual Ostracod Darwinula stevensoni

How asexual reproduction shapes transposable element (TE) content and diversity in eukaryotic genomes remains debated. We performed an initial survey of TE load and diversity in the putative ancient asexual ostracod Darwinula stevensoni. We examined long contiguous stretches of DNA in clones from a genomic fosmid library, totaling about 2.5 Mb, and supplemented these data with results on TE abundance and diversity from an Illumina draft genome. In contrast to other TE studies in putatively ancient asexuals, which revealed relatively low TE content, we found that at least 19% of the fosmid dataset and 26% of the genome assembly corresponded to known transposons. We observed a high diversity of transposon families, including LINE, gypsy, PLE, mariner/Tc, hAT, CMC, Sola2, Ginger, Merlin, Harbinger, MITEs and helitrons, with the prevalence of DNA transposons. The predominantly low levels of sequence diversity indicate that many TEs are or have recently been active. In the fosmid data, no correlation was found between telomeric repeats and non-LTR retrotransposons, which are present near telomeres in other taxa. Most TEs in the fosmid data were located outside of introns and almost none were found in exons. We also report an N-terminal Myb/SANT-like DNA-binding domain in site-specific R4/Dong non-LTR retrotransposons. Although initial results on transposable loads need to be verified with high quality draft genomes, this study provides important first insights into TE dynamics in putative ancient asexual ostracods.

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.

A genomic survey of Tc1-mariner transposons in nematodes suggests extensive horizontal transposon transfer events

The number of reports concerning horizontal transposon transfers (HTT) in metazoan species is considerably increased, alongside with the exponential growth of genomic sequence data However, our understanding of the mechanisms of such phenomenon is still at an early stage. Nematodes constitute an animal phylum successfully adapted to almost every ecosystem and for this reason could potentially contribute to spreading the genetic information through horizontal transfer. To date, few studies describe HTT of nematode retrotransposons. This is due to the lack of annotation of transposable elements in the sequenced nematode genomes, especially DNA transposons, which are acknowledged as the best horizontal travelers among mobile sequences. We have therefore started a survey of DNA transposons and their possible involvement in HTT in sequenced nematode genomes. Here, we describe 83 new Tc1/mariner elements distributed in 17 nematode species. Among them, nine families were possibly horizontally transferred between nematodes and the most diverse animal species, including ants as preferred partner of HTT. The results obtained suggest that HTT events involving nematodes Tc1/mariner elements are not uncommon, and that nematodes could have a possible role as transposon reservoir that, in turn, can be redistributed among animal genomes. Overall, this could be relevant to understand how the inter-species genetic flows shape the landscape of genetic variation of organisms inhabiting specific environmental communities.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

Long-term experimental evolution reveals purifying selection on piRNA-mediated control of transposable element expression

Background

Transposable elements (TEs) are an almost universal constituent of eukaryotic genomes. In animals, Piwi-interacting small RNAs (piRNAs) and repressive chromatin often play crucial roles in preventing TE transcription and thus restricting TE activity. Nevertheless, TE content varies widely across eukaryotes and the dynamics of TE activity and TE silencing across evolutionary time is poorly understood.

Results

Here, we used experimentally evolved populations of C. elegans to study the dynamics of TE expression over 409 generations. The experimental populations were evolved at population sizes of 1, 10 and 100 individuals to manipulate the efficiency of natural selection versus genetic drift. We demonstrate increased TE expression relative to the ancestral population, with the largest increases occurring in the smallest populations. We show that the transcriptional activation of TEs within active regions of the genome is associated with failure of piRNA-mediated silencing, whilst desilenced TEs in repressed chromatin domains retain small RNAs. Additionally, we find that the sequence context of the surrounding region influences the propensity of TEs to lose silencing through failure of small RNA-mediated silencing.

Conclusions

Our results show that natural selection in C. elegans is responsible for maintaining low levels of TE expression, and provide new insights into the epigenomic features responsible.

Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode

Root-knot nematodes (genus Meloidogyne) are plant parasites causing huge economic loss in the agricultural industry and affecting severely numerous developing countries. Control methods against these plant pests are sparse, the preferred one being the deployment of plant cultivars bearing resistance genes against Meloidogyne species. However, M. enterolobii is not controlled by the resistance genes deployed in the crop plants cultivated in Europe. The recent identification of this species in Europe is thus a major concern. Here, we sequenced the genome of M. enterolobii using short and long-read technologies. The genome assembly spans 240 Mbp with contig N50 size of 143 kbp, enabling high-quality annotations of 59,773 coding genes, 4,068 non-coding genes, and 10,944 transposable elements (spanning 8.7% of the genome). We validated the genome size by flow cytometry and the structure, quality and completeness by bioinformatics metrics. This ensemble of resources will fuel future projects aiming at pinpointing the genome singularities, the origin, diversity, and adaptive potential of this emerging plant pest.

Measurement(s)	genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing assay • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Meloidogyne enterolobii

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12410363

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.

Complex Evolutionary History of Mboumar, a Mariner Element Widely Represented in Ant Genomes

Mboumar-9 is an active mariner-transposable element previously isolated in the ant Messor bouvieri. In this work, a mariner-like element, Mboumar, isolated from 22 species of ants, is analyzed. These species belong to nine different subfamilies, including Leptanillinae, the most primitive ant subfamily, and Myrmicinae and Formicidae, the most derived ones. Consequently, Mboumar-like elements seem to be well-represented in ant genomes. The phylogenetic tree drawn for mariner elements is highly inconsistent with the phylogeny of host ants, with almost identical elements found in clearly distant species and, on the contrary, more variable elements in closely related species. The inconsistency between the two phylogenetic trees indicates that these transposable elements have evolved independently from the speciation events of the ants that host them. Besides, we found closer genetic relationships among elements than among their host ants. We also found potential coding copies with an uninterrupted open reading frame of 345 aa in 11 species. The putative transposase codified by them showed a high sequence identity with the active Mboumar-9 transposase. The results of selection tests suggest the intervention of purifying selection in the evolution of these elements. Overall, our study suggests a complex evolutionary history of the Mboumar-like mariner in ants, with important participation of horizontal transfer events. We also suggest that the evolutionary dynamics of Mboumar-like elements can be influenced by the genetic system of their host ants, which are eusocial insects with a haplodiploid genetic system.

Small RNAs in parasitic nematodes - forms and functions

Small RNAs are important regulators of gene expression. They were first identified in Caenorhabditis elegans, but it is now apparent that the main small RNA silencing pathways are functionally conserved across diverse organisms. Availability of genome data for an increasing number of parasitic nematodes has enabled bioinformatic identification of small RNA sequences. Expression of these in different lifecycle stages is revealed by small RNA sequencing and microarray analysis. In this review we describe what is known of the three main small RNA classes in parasitic nematodes – microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs) and small interfering RNAs (siRNAs) – and their proposed functions. miRNAs regulate development in C. elegans and the temporal expression of parasitic nematode miRNAs suggest modulation of target gene levels as parasites develop within the host. miRNAs are also present in extracellular vesicles released by nematodes in vitro, and in plasma from infected hosts, suggesting potential regulation of host gene expression. Roles of piRNAs and siRNAs in suppressing target genes, including transposable elements, are also reviewed. Recent successes in RNAi-mediated gene silencing, and application of small RNA inhibitors and mimics will continue to advance understanding of small RNA functions within the parasite and at the host–parasite interface.

Asexual reproduction reduces transposable element load in experimental yeast populations

Theory predicts that sexual reproduction can either facilitate or restrain transposable element (TE) accumulation by providing TEs with a means of spreading to all individuals in a population, versus facilitating TE load reduction via purifying selection. By quantifying genomic TE loads over time in experimental sexual and asexual Saccharomyces cerevisiae populations, we provide direct evidence that TE loads decrease rapidly under asexual reproduction. We show, using simulations, that this reduction may occur via evolution of TE activity, most likely via increased excision rates. Thus, sex is a major driver of genomic TE loads and at the root of the success of TEs.

The genetic information of most living organisms contains parasitic invaders known as transposable elements. These genetic sequences multiply by copying and pasting themselves through the genome, but this process can disrupt the activity of important genes and put the organism at risk.

How transposable elements proliferate in a population depends on the way organisms reproduce. If they simply clone themselves asexually, the selfish elements cannot spread between the different clones. If the organisms mate together their respective transposable elements get mixed, which helps the sequences to spread more easily and to potentially become more virulent. However, sexual reproduction also comes with mechanisms that keep transposable elements in check.

Bast, Jaron et al. took advantage of the fact that yeasts can reproduce with or without mating to explore whether sexual or asexual organisms are better at controlling the spread of transposable elements. The number of copies of transposable elements in the genomes of yeast grown sexually or asexually was assessed. The results showed that sexual populations kept constant numbers of selfish elements, while asexual organisms lost these genomic parasites over time. Simulations then revealed that this difference emerged because a defense gene that helps to delete transposable elements was spreading more quickly in the asexual group.

The work by Bast, Jaron et al. therefore suggests that sex is responsible for the evolutionary success of transposable elements, while asexual populations can discard these sequences over time. Sex therefore helps genetic parasites, somewhat similar to sexually transmitted diseases, to spread between individuals and remain virulent.

Diversification of Transposable Elements in Arthropods and Its Impact on Genome Evolution

Transposable elements (TEs) are ubiquitous in arthropods. However, analyses of large-scale and long-term coevolution between TEs and host genomes remain scarce in arthropods. Here, we choose 14 representative Arthropoda species from eight orders spanning more than 500 million years of evolution. By developing an unbiased TE annotation pipeline, we obtained 87 to 2266 TE reference sequences in a species, which is a considerable improvement compared to the reference TEs previously annotated in Repbase. We find that TE loads are diversified among species and were previously underestimated. The highly species- and time-specific expansions and contractions, and intraspecific sequence diversification are the leading driver of long terminal repeat (LTR) dynamics in Lepidoptera. Terminal inverted repeats (TIRs) proliferated substantially in five species with large genomes. A phylogenetic comparison reveals that the loads of multiple TE subfamilies are positively correlated with genome sizes. We also identified a few horizontally transferred TE candidates across nine species. In addition, we set up the Arthropod Transposable Elements database (ArTEdb) to provide TE references and annotations. Collectively, our results provide high-quality TE references and uncover that TE loads and expansion histories vary greatly among arthropods, which implies that TEs are an important driving force shaping the evolution of genomes through gain and loss.

Neutral Theory, Transposable Elements, and Eukaryotic Genome Evolution

Among the multitude of papers published yearly in scientific journals, precious few publications may be worth looking back in half a century to appreciate the significance of the discoveries that would later become common knowledge and get a chance to shape a field or several adjacent fields. Here, Kimura's fundamental concept of neutral mutation-random drift, which was published 50 years ago, is re-examined in light of its pervasive influence on comparative genomics and, more specifically, on the contribution of transposable elements to eukaryotic genome evolution.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

The Emergence of the Tardigrade Hypsibius exemplaris as a Model System

The success of scientists in revealing biological mechanisms has depended in large part on choosing tractable model systems. In 1997, molecular phylogenetics revealed that two of biology's most tractable models-Caenorhabditis elegans and Drosophila-are much more closely related to each other than had been thought previously. I began to explore whether any of the little-studied members of this branch of the tree of life might serve as a new model for comparative biology that could make use of the rich and ongoing sources of information flowing from C. elegans and Drosophila research. Tardigrades, also known as water bears, make up a phylum of microscopic animals. The tardigrade Hypsibius exemplaris (recently disambiguated from a closely related species, Hypsibius dujardini) can be maintained in laboratories and has a generation time of <2 wk at room temperature. Stocks of animals can be stored frozen and revived. The animals and their embryos are optically clear, and embryos are laid in groups, with each synchronous clutch of embryos laid in a clear molt. We have developed techniques for laboratory study of this system, including methods for microinjection of animals, immunolocalization, in situ hybridization, RNA interference, transcriptomics, and methods for identifying proteins that mediate tolerance to extreme environments. Here, I review the development of this animal as an emerging model system, as well as recent molecular studies aimed at understanding the evolution of developmental mechanisms that underpin the evolution of animal form and at understanding how biological materials can survive extreme environments.

Modeling Interactions between Transposable Elements and the Plant Epigenetic Response: A Surprising Reliance on Element Retention

Transposable elements (TEs) compose the majority of angiosperm DNA. Plants counteract TE activity by silencing them epigenetically. One form of epigenetic silencing requires 21-22 nt small interfering RNAs that act to degrade TE mRNA and may also trigger DNA methylation. DNA methylation is reinforced by a second mechanism, the RNA-dependent DNA methylation (RdDM) pathway. RdDM relies on 24 nt small interfering RNAs and ultimately establishes TEs in a quiescent state. These host factors interact at a systems level, but there have been no system level analyses of their interactions. Here, we define a deterministic model that represents the propagation of active TEs, aspects of the host response and the accumulation of silenced TEs. We describe general properties of the model and also fit it to biological data in order to explore two questions. The first is why two overlapping pathways are maintained, given that both are likely energetically expensive. Under our model, RdDM silenced TEs effectively even when the initiation of silencing was weak. This relationship implies that only a small amount of RNAi is needed to initiate TE silencing, but reinforcement by RdDM is necessary to efficiently counter TE propagation. Second, we investigated the reliance of the host response on rates of TE deletion. The model predicted that low levels of deletion lead to few active TEs, suggesting that silencing is most efficient when methylated TEs are retained in the genome, thereby providing one explanation for the large size of plant genomes.

Comparative Genomics of Apomictic Root-Knot Nematodes: Hybridization, Ploidy, and Dynamic Genome Change

The root-knot nematodes (genus Meloidogyne) are important plant parasites causing substantial agricultural losses. The Meloidogyne incognita group (MIG) of species, most of which are obligatory apomicts (mitotic parthenogens), are extremely polyphagous and important problems for global agriculture. While understanding the genomic basis for their variable success on different crops could benefit future agriculture, analyses of their genomes are challenging due to complex evolutionary histories that may incorporate hybridization, ploidy changes, and chromosomal fragmentation. Here, we sequence 19 genomes, representing five species of key root-knot nematodes collected from different geographic origins. We show that a hybrid origin that predated speciation within the MIG has resulted in each species possessing two divergent genomic copies. Additionally, the apomictic MIG species are hypotriploids, with a proportion of one genome present in a second copy. The hypotriploid proportion varies among species. The evolutionary history of the MIG genomes is revealed to be very dynamic, with noncrossover recombination both homogenizing the genomic copies, and acting as a mechanism for generating divergence between species. Interestingly, the automictic MIG species M. floridensis differs from the apomict species in that it has become homozygous throughout much of its genome.

Evolutionary divergence of 3' UTRs in cichlid fishes

Background

Post-transcriptional regulation is crucial for the control of eukaryotic gene expression and might contribute to adaptive divergence. The three prime untranslated regions (3’ UTRs), that are located downstream of protein-coding sequences, play important roles in post-transcriptional regulation. These regions contain functional elements that influence the fate of mRNAs and could be exceptionally important in groups such as rapidly evolving cichlid fishes.

Results

To examine cichlid 3’ UTR evolution, we 1) identified gene features in nine teleost genomes and 2) performed comparative analyses to assess evolutionary variation in length, functional motifs, and evolutionary rates of 3’ UTRs. In all nine teleost genomes, we found a smaller proportion of repetitive elements in 3’ UTRs than in the whole genome. We found that the 3’ UTRs in cichlids tend to be longer than those in non-cichlids, and this was associated, on average, with one more miRNA target per gene in cichlids. Moreover, we provided evidence that 3’ UTRs on average have evolved faster in cichlids than in non-cichlids. Finally, analyses of gene function suggested that both the top 5% longest and 5% most rapidly evolving 3’ UTRs in cichlids tended to be involved in ribosome-associated pathways and translation.

Conclusions

Our results reveal novel patterns of evolution in the 3’ UTRs of teleosts in general and cichlids in particular. The data suggest that 3’ UTRs might serve as important meta-regulators, regulators of other mechanisms governing post-transcriptional regulation, especially in groups like cichlids that have undergone extremely fast rates of phenotypic diversification and speciation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4821-8) contains supplementary material, which is available to authorized users.

Transposable elements and polyploid evolution in animals

Polyploidy in animals is much less common than in plants, where it is thought to be pervasive in all higher plant lineages. Recent studies have highlighted the impact of polyploidization and the associated process of diploidy restoration on the evolution and speciation of selected taxonomic groups in the animal kingdom: from vertebrates represented by salmonid fishes and African clawed frogs to invertebrates represented by parasitic root-knot nematodes and bdelloid rotifers. In this review, we focus on the unique and diverse roles that transposable elements may play in these processes, from marking and diversifying subgenome-specific chromosome sets before hybridization, to influencing genome restructuring during rediploidization, to affecting subgenome-specific regulatory evolution, and occasionally providing opportunities for domestication and gene amplification to restore and improve functionality. There is still much to be learned from the future comparative genomic studies of chromosome-sized and haplotype-aware assemblies, and from postgenomic studies elucidating genetic and epigenetic regulatory phenomena across short and long evolutionary distances in the metazoan tree of life.

The gene regulatory program of Acrobeloides nanus reveals conservation of phylum-specific expression

The evolution of development has been studied through the lens of gene regulation by examining either closely related species or extremely distant animals of different phyla. In nematodes, detailed cell- and stage-specific expression analyses are focused on the model Caenorhabditis elegans, in part leading to the view that the developmental expression of gene cascades in this species is archetypic for the phylum. Here, we compared two species of an intermediate evolutionary distance: the nematodes C. elegans (clade V) and Acrobeloides nanus (clade IV). To examine A. nanus molecularly, we sequenced its genome and identified the expression profiles of all genes throughout embryogenesis. In comparison with C. elegans, A. nanus exhibits a much slower embryonic development and has a capacity for regulative compensation of missing early cells. We detected conserved stages between these species at the transcriptome level, as well as a prominent middevelopmental transition, at which point the two species converge in terms of their gene expression. Interestingly, we found that genes originating at the dawn of the Ecdysozoa supergroup show the least expression divergence between these two species. This led us to detect a correlation between the time of expression of a gene and its phylogenetic age: evolutionarily ancient and young genes are enriched for expression in early and late embryogenesis, respectively, whereas Ecdysozoa-specific genes are enriched for expression during the middevelopmental transition. Our results characterize the developmental constraints operating on each individual embryo in terms of developmental stages and genetic evolutionary history.

Dynamics of genome size evolution in birds and mammals

Genome size in mammals and birds shows remarkably little interspecific variation compared with other taxa. However, genome sequencing has revealed that many mammal and bird lineages have experienced differential rates of transposable element (TE) accumulation, which would be predicted to cause substantial variation in genome size between species. Thus, we hypothesize that there has been covariation between the amount of DNA gained by transposition and lost by deletion during mammal and avian evolution, resulting in genome size equilibrium. To test this model, we develop computational methods to quantify the amount of DNA gained by TE expansion and lost by deletion over the last 100 My in the lineages of 10 species of eutherian mammals and 24 species of birds. The results reveal extensive variation in the amount of DNA gained via lineage-specific transposition, but that DNA loss counteracted this expansion to various extents across lineages. Our analysis of the rate and size spectrum of deletion events implies that DNA removal in both mammals and birds has proceeded mostly through large segmental deletions (>10 kb). These findings support a unified "accordion" model of genome size evolution in eukaryotes whereby DNA loss counteracting TE expansion is a major determinant of genome size. Furthermore, we propose that extensive DNA loss, and not necessarily a dearth of TE activity, has been the primary force maintaining the greater genomic compaction of flying birds and bats relative to their flightless relatives.

Ecological networks to unravel the routes to horizontal transposon transfers

Transposable elements (TEs) represent the single largest component of numerous eukaryotic genomes, and their activity and dispersal constitute an important force fostering evolutionary innovation. The horizontal transfer of TEs (HTT) between eukaryotic species is a common and widespread phenomenon that has had a profound impact on TE dynamics and, consequently, on the evolutionary trajectory of many species' lineages. However, the mechanisms promoting HTT remain largely unknown. In this article, we argue that network theory combined with functional ecology provides a robust conceptual framework and tools to delineate how complex interactions between diverse organisms may act in synergy to promote HTTs.