Patterns of Selection Against Transposons Inferred From the Distribution of Tc1, Tc3 and Tc5 Insertions in the mut-7 Line of the Nematode Caenorhabditis elegans

Bioinformatics and Machine Learning Approaches to Understand the Regulation of Mobile Genetic Elements

Simple Summary

Transposable elements (TEs) are DNA sequences that are, or were, able to move (transpose) within the genome of a single cell. They were first discovered by Barbara McClintock while working on maize, and they make up a large fraction of the genome. Transpositions can result in mutations and they can alter the genome size. Cells regulate the activity of TEs using a variety of mechanisms, such as chemical modifications of DNA and small RNAs. Machine learning (ML) is an interdisciplinary subject that studies computer algorithms that can improve through experience and by the use of data. ML has been successfully applied to a variety of problems in bioinformatics and has exhibited favorable precision and speed. Here, we provide a systematic and guided review on the ML and bioinformatic methods and tools that are used for the analysis of the regulation of TEs.

Abstract

Transposable elements (TEs, or mobile genetic elements, MGEs) are ubiquitous genetic elements that make up a substantial proportion of the genome of many species. The recent growing interest in understanding the evolution and function of TEs has revealed that TEs play a dual role in genome evolution, development, disease, and drug resistance. Cells regulate TE expression against uncontrolled activity that can lead to developmental defects and disease, using multiple strategies, such as DNA chemical modification, small RNA (sRNA) silencing, chromatin modification, as well as sequence-specific repressors. Advancements in bioinformatics and machine learning approaches are increasingly contributing to the analysis of the regulation mechanisms. A plethora of tools and machine learning approaches have been developed for prediction, annotation, and expression profiling of sRNAs, for methylation analysis of TEs, as well as for genome-wide methylation analysis through bisulfite sequencing data. In this review, we provide a guided overview of the bioinformatic and machine learning state of the art of fields closely associated with TE regulation and function.

Coevolution between transposable elements and recombination

One of the most striking patterns of genome structure is the tight, typically negative, association between transposable elements (TEs) and meiotic recombination rates. While this is a highly recurring feature of eukaryotic genomes, the mechanisms driving correlations between TEs and recombination remain poorly understood, and distinguishing cause versus effect is challenging. Here, we review the evidence for a relation between TEs and recombination, and discuss the underlying evolutionary forces. Evidence to date suggests that overall TE densities correlate negatively with recombination, but the strength of this correlation varies across element types, and the pattern can be reversed. Results suggest that heterogeneity in the strength of selection against ectopic recombination and gene disruption can drive TE accumulation in regions of low recombination, but there is also strong evidence that the regulation of TEs can influence local recombination rates. We hypothesize that TE insertion polymorphism may be important in driving within-species variation in recombination rates in surrounding genomic regions. Furthermore, the interaction between TEs and recombination may create positive feedback, whereby TE accumulation in non-recombining regions contributes to the spread of recombination suppression. Further investigation of the coevolution between recombination and TEs has important implications for our understanding of the evolution of recombination rates and genome structure.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Natural Variation in the Distribution and Abundance of Transposable Elements Across the Caenorhabditis elegans Species

Transposons are mobile DNA elements that generate both adaptive and deleterious phenotypic variation thereby driving genome evolution. For these reasons, genomes have mechanisms to regulate transposable element (TE) activity. Approximately 12–16% of the Caenorhabditis elegans genome is composed of TEs, of which the majority are likely inactive. However, most studies of TE activity have been conducted in the laboratory strain N2, which limits our knowledge of the effects of these mobile elements across natural populations. We analyzed the distribution and abundance of TEs in 208 wild C. elegans strains to better understand how transposons contribute to variation in natural populations. We identified 3,397 TEs as compared with the reference strain, of which 2,771 are novel insertions and 241 are TEs that have been excised in at least one wild strain. Likely because of their hypothesized deleterious effects, we find that TEs are found at low allele frequencies throughout the population, and we predict functional effects of TE insertions. The abundances of TEs reflect their activities, and these data allowed us to perform both genome-wide association mappings and rare variant correlations to reveal several candidate genes that impact TE regulation, including small regulatory piwi-interacting RNAs and chromatin factors. Because TE variation in natural populations could underlie phenotypic variation for organismal and behavioral traits, the transposons that we identified and their regulatory mechanisms can be used in future studies to explore the genomics of complex traits and evolutionary changes.

Regulation of transposable elements: Interplay between TE-encoded regulatory sequences and host-specific trans-acting factors in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that can move around the genome, and their expression is one precondition for this mobility. Because the insertion of TEs in new genomic positions is largely deleterious, the molecular mechanisms for transcriptional suppression have been extensively studied. In contrast, very little is known about their primary transcriptional regulation. Here, we characterize the expression dynamics of TE families in Drosophila melanogaster across a broad temperature range (13-29°C). In 71% of the expressed TE families, the expression is modulated by temperature. We show that this temperature-dependent regulation is specific for TE families and strongly affected by the genetic background. We deduce that TEs carry family-specific regulatory sequences, which are targeted by host-specific trans-acting factors, such as transcription factors. Consistent with the widespread dominant inheritance of gene expression, we also find the prevailing dominance of TE family expression. We conclude that TE family expression across a range of temperatures is regulated by an interaction between TE family-specific regulatory elements and trans-acting factors of the host.

Evolution of a transposon in Daphnia hybrid genomes

Background

Transposable elements play a major role in genome evolution. Their capacity to move and/or multiply in the genome of their host may have profound impacts on phenotypes, and may have dramatic consequences on genome structure. Hybrid and polyploid clones have arisen multiple times in the Daphnia pulex complex and are thought to reproduce by obligate parthenogenesis. Our study examines the evolution of a DNA transposable element named Pokey in the D. pulex complex.

Results

Portions of Pokey elements inserted in the 28S rRNA genes from various Daphnia hybrids (diploids and polyploids) were sequenced and compared to sequences from a previous study to understand the evolutionary history of the elements. Pokey sequences show a complex phylogenetic pattern. We found evidence of recombination events in numerous Pokey alleles from diploid and polyploid hybrids and also from non-hybrid diploids. The recombination rate in Pokey elements is comparable to recombination rates previously estimated for 28S rRNA genes in the congener, Daphnia obtusa. Some recombinant Pokey alleles were encountered in Daphnia isolates from multiple locations and habitats.

Conclusions

Phylogenetic and recombination analyses showed that recombination is a major force that shapes Pokey evolution. Based on Pokey phylogenies, reticulation has played and still plays an important role in shaping the diversity of the D. pulex complex. Horizontal transfer of Pokey seems to be rare and hybrids often possess Pokey elements derived from recombination among alleles encountered in the putative parental species. The insertion of Pokey in hotspots of recombination may have important impacts on the diversity and fitness of this transposable element.

Population genomics of transposable elements in Drosophila melanogaster

Transposable elements (TEs) are the primary contributors to the genome bulk in many organisms and are major players in genome evolution. A clear and thorough understanding of the population dynamics of TEs is therefore essential for full comprehension of the eukaryotic genome evolution and function. Although TEs in Drosophila melanogaster have received much attention, population dynamics of most TE families in this species remains entirely unexplored. It is not clear whether the same population processes can account for the population behaviors of all TEs in Drosophila or whether, as has been suggested previously, different orders behave according to very different rules. In this work, we analyzed population frequencies for a large number of individual TEs (755 TEs) in five North American and one sub-Saharan African D. melanogaster populations (75 strains in total). These TEs have been annotated in the reference D. melanogaster euchromatic genome and have been sampled from all three major orders (non-LTR, LTR, and TIR) and from all families with more than 20 TE copies (55 families in total). We find strong evidence that TEs in Drosophila across all orders and families are subject to purifying selection at the level of ectopic recombination. We showed that strength of this selection varies predictably with recombination rate, length of individual TEs, and copy number and length of other TEs in the same family. Importantly, these rules do not appear to vary across orders. Finally, we built a statistical model that considered only individual TE-level (such as the TE length) and family-level properties (such as the copy number) and were able to explain more than 40% of the variation in TE frequencies in D. melanogaster.

A triptych of the evolution of plant transposable elements

Transposable elements (TEs) constitute the majority of angiosperm DNA, but the processes that govern their accumulation remain mysterious. Here we discuss the three major forces that govern the accumulation of TEs, corresponding to the three panels of a triptych. The first force, transposition, creates new copies of TEs, but is regulated by both host- and TE-specific mechanisms. The second force, deletion of TE DNA, is capable of removing vast swaths of genomic regions via recombinational processes, but we still have very little insight into how deletion varies across species and even among TE types. Finally, we focus on the often-ignored third panel of our triptych - the population processes that determine the ultimate evolutionary fate of TE insertions.

Genome-wide analysis of Tol2 transposon reintegration in zebrafish

BACKGROUND

Tol2, a member of the hAT family of transposons, has become a useful tool for genetic manipulation of model animals, but information about its interactions with vertebrate genomes is still limited. Furthermore, published reports on Tol2 have mainly been based on random integration of the transposon system after co-injection of a plasmid DNA harboring the transposon and a transposase mRNA. It is important to understand how Tol2 would behave upon activation after integration into the genome.

RESULTS

We performed a large-scale enhancer trap (ET) screen and generated 338 insertions of the Tol2 transposon-based ET cassette into the zebrafish genome. These insertions were generated by remobilizing the transposon from two different donor sites in two transgenic lines. We found that 39% of Tol2 insertions occurred in transcription units, mostly into introns. Analysis of the transposon target sites revealed no strict specificity at the DNA sequence level. However, Tol2 was prone to target AT-rich regions with weak palindromic consensus sequences centered at the insertion site.

CONCLUSION

Our systematic analysis of sequential remobilizations of the Tol2 transposon from two independent sites within a vertebrate genome has revealed properties such as a tendency to integrate into transcription units and into AT-rich palindrome-like sequences. This information will influence the development of various applications involving DNA transposons and Tol2 in particular.

Epigenetic silencing of transposable elements: a trade-off between reduced transposition and deleterious effects on neighboring gene expression

Transposable elements (TEs) are ubiquitous genomic parasites. The deleterious consequences of the presence and activity of TEs have fueled debate about the evolutionary forces countering their expansion. Purifying selection is thought to purge TE insertions from the genome, and TE sequences are targeted by hosts for epigenetic silencing. However, the interplay between epigenetic and evolutionary forces countering TE expansion remains unexplored. Here we analyze genomic, epigenetic, and population genetic data from Arabidopsis thaliana to yield three observations. First, gene expression is negatively correlated with the density of methylated TEs. Second, the signature of purifying selection is detectable for methylated TEs near genes but not for unmethylated TEs or for TEs far from genes. Third, TE insertions are distributed by age and methylation status, such that older, methylated TEs are farther from genes. Based on these observations, we present a model in which host silencing of TEs near genes has deleterious effects on neighboring gene expression, resulting in the preferential loss of methylated TEs from gene-rich chromosomal regions. This mechanism implies an evolutionary tradeoff in which the benefit of TE silencing imposes a fitness cost via deleterious effects on the expression of nearby genes.

The evolutionary dynamics of the Helena retrotransposon revealed by sequenced Drosophila genomes

BACKGROUND

Several studies have shown that genomes contain a mixture of transposable elements, some of which are still active and others ancient relics that have degenerated. This is true for the non-LTR retrotransposon Helena, of which only degenerate sequences have been shown to be present in some species (Drosophila melanogaster), whereas putatively active sequences are present in others (D. simulans). Combining experimental and population analyses with the sequence analysis of the 12 Drosophila genomes, we have investigated the evolution of Helena, and propose a possible scenario for the evolution of this element.

RESULTS

We show that six species of Drosophila have the Helena transposable element at different stages of its evolution. The copy number is highly variable among these species, but most of them are truncated at the 5' ends and also harbor several internal deletions and insertions suggesting that they are inactive in all species, except in D. mojavensis in which quantitative RT-PCR experiments have identified a putative active copy.

CONCLUSION

Our data suggest that Helena was present in the common ancestor of the Drosophila genus, which has been vertically transmitted to the derived lineages, but that it has been lost in some of them. The wide variation in copy number and sequence degeneration in the different species suggest that the evolutionary dynamics of Helena depends on the genomic environment of the host species.

Identification of expressed transposable element insertions in the sequenced genome of Drosophila melanogaster

Transposable elements (TEs) are major components of most genomes, and their impact on genome evolution is now well documented. However, the way they affect the transcriptome is still not clearly established. Using the sequenced genome of Drosophila melanogaster and EST libraries, we describe here the TE insertions that are unequivocally transcribed, and we have determined their location in the sequenced genome. We show that most TE families are transcribed, and we have specifically identified 69 expressed TE insertions, half of which are located inside genes, mostly within introns and 5'UTRs.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

Population frequencies of transposable elements in selfing and outcrossing Caenorhabditis nematodes

Population genetics theory predicts that differences in breeding systems should be an important factor in the dynamics of selfish genetic elements, because of different intensities of selection on both hosts and elements. We examined population frequencies of transposable elements (TEs) in natural populations of the self-fertilizing nematode Caenorhabditis elegans and its outcrossing relative Caenorhabditis remanei. We identified a Tc1-like class of elements in the C. remanei genome with homology to the terminal inverted repeats of the C. elegans Tc1 transposon, which we name mTcre1. We measured levels of insertion polymorphism for all 32 Tc1 elements present in the genome sequence of the C. elegans N2 strain, and 16 mTcre1 elements from the genome sequence of the C. remanei PB4641 strain. We show that transposons are less polymorphic and segregate at higher frequencies in C. elegans compared with C. remanei. Estimates of the intensity of selection based on the population frequencies of polymorphic elements suggest that transposons are selectively neutral in C. elegans, but subject to purifying selection in C. remanei. These results are consistent with a reduced efficacy of natural selection against TEs in selfing populations, but may in part be explained by non-equilibrium TE dynamics.

Genome-wide screening and characterization of transposable elements and their distribution analysis in the silkworm, Bombyx mori

To elucidate the contribution of transposable elements (TEs) to the silkworm genome structure and evolution, we have conducted genome-wide analysis of TEs using the newly released genome assembly. The TEs made up 35% of the genome and contributed greatly to the genome size. Non-long terminal repeat retrotransposons (non-LTRs) and short interspersed nuclear elements (SINEs) were the predominant TE classes. From characterization of the TE distribution in the genome, it was revealed that non-LTRs, especially R1 clade elements, are frequently inserted into GC-rich regions. The GC content of non-LTRs themselves was over 40%, which indicate their contribution to the GC content of the insertion region. TEs accumulated in regions with low gene density, and there were relatively strong positive correlations between TE density and chromosomal recombination rate. We also characterized the clade distribution of the non-LTRs. The silkworm non-LTRs represented 10 of the 16 previously defined clades, which had the most variety than that reported for other genomes. Two partial CRE clade elements were found, which is one of the most ancient lineages of non-LTRs, and have been only found in Trypanosoma and fungi before. This analysis suggests that Bombyx genome is influenced by numerous amounts and variety of TEs.

DNA transposons and the evolution of eukaryotic genomes

Transposable elements are mobile genetic units that exhibit broad diversity in their structure and transposition mechanisms. Transposable elements occupy a large fraction of many eukaryotic genomes and their movement and accumulation represent a major force shaping the genes and genomes of almost all organisms. This review focuses on DNA-mediated or class 2 transposons and emphasizes how this class of elements is distinguished from other types of mobile elements in terms of their structure, amplification dynamics, and genomic effect. We provide an up-to-date outlook on the diversity and taxonomic distribution of all major types of DNA transposons in eukaryotes, including Helitrons and Mavericks. We discuss some of the evolutionary forces that influence their maintenance and diversification in various genomic environments. Finally, we highlight how the distinctive biological features of DNA transposons have contributed to shape genome architecture and led to the emergence of genetic innovations in different eukaryotic lineages.

The heterochromatic copies of the LTR retrotransposons as a record of the genomic events that have shaped the Drosophila melanogaster genome

Transposable elements, which are major components of most genomes, are known to accumulate in heterochromatic regions in which they have progressively diverged in sequence by mutations and internal deletions and insertions (indels) during the course of evolution. They therefore provide a record of the genomic events that have shaped the genomes, some of which could correspond to speciation events. Using the sequence divergence between the long terminal repeats (LTRs), we estimated the date of the insertion events of the LTR retrotransposon copies embedded within the heterochromatin regions of the Drosophila melanogaster genome. We did not detect traces of any specific waves of mobilization of retrotransposons within heterochromatin, apart from a very recent wave, which corresponds to the numerous LTR retrotransposon copies found in euchromatin.

The distribution of copia-type retrotransposons and the evolutionary history of tomato and related wild species

Retrotransposons are mobile genetic elements that amplify throughout the genome and may be important contributors of genetic diversity. Their distribution is influenced by element behaviour and host-driven controls. We analysed the distribution of three copia-type retrotransposons, ToRTL1, T135 and Tnt1 using sequence-specific amplification polymorphism in self-compatible (SC) and incompatible (SI) species of Solanum subsection Lycopersicon, and genetically mapped polymorphic insertions in S. lycopersicum (tomato). The majority of polymorphic insertions (61%) are located in centromeric regions of the tomato genome. A significant positive relationship was detected between insertion polymorphisms and mating system, independent of selection as most insertions were found to be neutral. As insertion patterns successfully inferred interspecific relationships of Solanum subsection Lycopersicon, our results suggest that the distribution of ToRTL1, T135 and Tnt1 may essentially be determined by selection removing strongly deleterious insertions, with genetic drift and mating system, but not recombination rate, playing important roles.

Recombination: an underappreciated factor in the evolution of plant genomes

Our knowledge of recombination rates and patterns in plants is far from being comprehensive. However, compelling evidence indicates a central role for recombination, through its influences on mutation and selection, in the evolution of plant genomes. Furthermore, recombination seems to be generally higher and more variable in plants than in animals, which could be one of the primary reasons for differences in genome lability between these two kingdoms. Much additional study of recombination in plants is needed to investigate these ideas further.

Transposable element orientation bias in the Drosophila melanogaster genome

Nonrandom distributions of transposable elements can be generated by a variety of genomic features. Using the full D. melanogaster genome as a model, we characterize the orientations of different classes of transposable elements in relation to the directionality of genes. DNA-mediated transposable elements are more likely to be in the same orientation as neighboring genes when they occur in the nontranscribed region's that flank genes. However, RNA-mediated transposable elements located in an intron are more often oriented in the direction opposite to that of the host gene. These orientation biases are strongest for genes with highly biased codon usage, probably reflecting the ability of such loci to respond to weak positive or negative selection. The leading hypothesis for selection against transposable elements in the coding orientation proposes that transcription termination poly(A) signal motifs within retroelements interfere with normal gene transcription. However, after accounting for differences in base composition between the strands, we find no evidence for global selection against spurious transcription termination signals in introns. We therefore conclude that premature termination of host gene transcription due to the presence of poly(A) signal motifs in retroelements might only partially explain strand-specific detrimental effects in the D. melanogaster genome.

Differences in genome size between closely related species: the Drosophila melanogaster species subgroup

Genome size varies considerably among organisms due to differences in the amplification, deletion, and divergence of various kinds of repetitive sequences, including the transposable elements, which constitute a large fraction of the genome. However, while the changes in genome size observed at a wide taxonomic level have been thoroughly investigated, we still know little about the process involved in closely related species. We estimated genome sizes and the reverse transcriptase-related sequence (RTRS) content in the nine species of the Drosophila melanogaster species subgroup. We showed that the species differ with regard to their genome size and that the RTRS content is correlated with genome size for all species except Drosophila orena. The genome of D. orena, which is 1.6-fold as big as that of D. melanogaster, has in fact not undergone any major increase in its RTRS content.

High-resolution genome-wide mapping of transposon integration in mammals

The Sleeping Beauty (SB) transposon is an emerging tool for transgenesis, gene discovery, and therapeutic gene delivery in mammals. Here we studied 1,336 SB insertions in primary and cultured mammalian cells in order to better understand its target site preferences. We report that, although widely distributed, SB integration recurrently targets certain genomic regions and shows a small but significant bias toward genes and their upstream regulatory sequences. Compared to those of most integrating viruses, however, the regional preferences associated with SB-mediated integration were much less pronounced and were not significantly influenced by transcriptional activity. Insertions were also distinctly nonrandom with respect to intergenic sequences, including a strong bias toward microsatellite repeats, which are predominantly enriched in noncoding DNA. Although we detected a consensus sequence consistent with a twofold dyad symmetry at the target site, the most widely used sites did not match this consensus. In conjunction with an observed SB integration preference for bent DNA, these results suggest that physical properties may be the major determining factor in SB target site selection. These findings provide basic insights into the transposition process and reveal important distinctions between transposon- and virus-based integrating vectors.

Deleterious transposable elements and the extinction of asexuals

The genomes of virtually all sexually reproducing species contain transposable elements. Although active elements generally transpose more rapidly than they are inactivated by mutation or excision, their number can be kept in check by purifying selection if its effectiveness becomes disproportionately greater as their copy number increases. In sexually reproducing species, such synergistic selection can result from ectopic crossing-over or from homologous recombination under negative epistasis. In addition, there may be controls on transposon activity that are associated with meiosis. Because a sexual lineage that abandons sex must lack such mechanisms, it may be driven to extinction by the unchecked proliferation of deleterious transposons inherited from its sexual progenitor. An important component of the evolutionary advantage of sex over asex may therefore lie in the ability of sex, despite facilitating the spread of deleterious elements within interbreeding populations, also to restrain their intragenomic proliferation.

Epigenetic control of CACTA transposon mobility in Arabidopsis thaliana

Epigenetic mutation, heritable developmental variation not based on a change in nucleotide sequence, is widely reported in plants. However, the developmental and evolutionary significance of such mutations remains enigmatic. On the basis of our studies of the endogenous Arabidopsis transposon CACTA, we propose that the inheritance of epigenetic gene silencing over generations can function as a transgenerational genome defense mechanism against deleterious movement of transposons. We previously reported that silent CACTA1 is mobilized by the DNA hypomethylation mutation ddm1 (decrease in DNA methylation). In this study, we report that CACTA activated by the ddm1 mutation remains mobile in the presence of the wild-type DDM1 gene, suggesting that de novo silencing is not efficient for the defense of the genome against CACTA movement. The defense depends on maintenance of transposon silencing over generations. In addition, we show that the activated CACTA1 element transposes throughout the genome in DDM1 plants, as reported previously for ddm1 backgrounds. Furthermore, the CACTA1 element integrated into both the ddm1-derived and the DDM1-derived chromosomal regions in the DDM1 wild-type plants, demonstrating that this class of transposons does not exhibit targeted integration into heterochromatin, despite its accumulation in the pericentromeric regions in natural populations. The possible contribution of natural selection as a mechanism for the accumulation of transposons and evolution of heterochromatin is discussed.