The evolution of transposable elements in natural populations of self-fertilizing Arabidopsis thaliana and its outcrossing relative Arabidopsis lyrata

Evolutionary divergence and limits of conserved non-coding sequence detection in plant genomes

The discovery of regulatory motifs embedded in upstream regions of plants is a particularly challenging bioinformatics task. Previous studies have shown that motifs in plants are short compared with those found in vertebrates. Furthermore, plant genomes have undergone several diversification mechanisms such as genome duplication events which impact the evolution of regulatory motifs. In this article, a systematic phylogenomic comparison of upstream regions is conducted to further identify features of the plant regulatory genomes, the component of genomes regulating gene expression, to enable future de novo discoveries. The findings highlight differences in upstream region properties between major plant groups and the effects of divergence times and duplication events. First, clear differences in upstream region evolution can be detected between monocots and dicots, thus suggesting that a separation of these groups should be made when searching for novel regulatory motifs, particularly since universal motifs such as the TATA box are rare. Second, investigating the decay rate of significantly aligned regions suggests that a divergence time of ~100 mya sets a limit for reliable conserved non-coding sequence (CNS) detection. Insights presented here will set a framework to help identify embedded motifs of functional relevance by understanding the limits of bioinformatics detection for CNSs.

A brief history of the status of transposable elements: from junk DNA to major players in evolution

The idea that some genetic factors are able to move around chromosomes emerged more than 60 years ago when Barbara McClintock first suggested that such elements existed and had a major role in controlling gene expression and that they also have had a major influence in reshaping genomes in evolution. It was many years, however, before the accumulation of data and theories showed that this latter revolutionary idea was correct although, understandably, it fell far short of our present view of the significant influence of what are now known as "transposable elements" in evolution. In this article, I summarize the main events that influenced my thinking about transposable elements as a young scientist and the influence and role of these specific genomic elements in evolution over subsequent years. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work.

Osvaldo and Isis retrotransposons as markers of the Drosophila buzzatii colonisation in Australia

Background

Transposable elements (TEs) constitute an important source of genetic variability owing to their jumping and regulatory properties, and are considered to drive species evolution. Several factors that are able to induce TE transposition in genomes have been documented (for example environmental stress and inter- and intra-specific crosses) but in many instances the reasons for TE mobilisation have yet to be elucidated. Colonising populations constitute an ideal model for studying TE behaviour and distribution as they are exposed to different environmental and new demographic conditions. In this study, the distribution of two TEs, Osvaldo and Isis, was examined in two colonising populations of D. buzzatii from Australia. Comparing Osvaldo copy numbers between Australian and Old World (reported in previous studies) colonisations provides a valuable tool for elucidating the colonisation process and the effect of new conditions encountered by colonisers on TEs.

Results

The chromosomal distributions of Osvaldo and Isis retrotransposons in two colonising populations of D. buzzatii from Australia revealed sites of high insertion frequency (>10%) and low frequency sites. Comparisons between Osvaldo insertion profiles in colonising populations from the Old World and Australia demonstrate a tendency towards a higher number of highly occupied sites with higher insertion frequency in the Old World than in Australian populations. Tests concerning selection against deleterious TE insertions indicate that Isis is more controlled by purifying selection than Osvaldo. The distribution of both elements on chromosomal arms follows a Poisson distribution and there are non-significant positive correlations between highly occupied sites and chromosomal inversions.

Conclusions

The occupancy profile of Osvaldo and Isis retrotransposons is characterised by the existence of high and low insertion frequency sites in the populations. These results demonstrate that Australian D. buzzatii populations were subjected to a founder effect during the colonisation process. Moreover, there are more sites with high insertion frequency in the Old World colonisation than in the Australian colonisation, indicating a probable stronger bottleneck effect in Australia. The results suggest that selection does not seem to play a major role, compared to demography, in the distribution of transposable elements in the Australian populations.

Consequences of reproductive mode on genome evolution in fungi

An organism's reproductive mode is believed to be a major factor driving its genome evolution. In theory, sexual inbreeding and asexuality are associated with lower effective recombination levels and smaller effective population sizes than sexual outbreeding, giving rise to reduced selection efficiency and genetic hitchhiking. This, in turn, is predicted to result in the accumulation of deleterious mutations and other genomic changes, for example the accumulation of repetitive elements. Empirical data from plants and animals supporting/refuting these theories are sparse and have yielded few conclusive results. A growing body of data from the fungal kingdom, wherein reproductive behavior varies extensively within and among taxonomic groups, has provided new insights into the role of mating systems (e.g., homothallism, heterothallism, pseudohomothallism) and asexuality, on genome evolution. Herein, we briefly review the theoretical relationships between reproductive mode and genome evolution and give examples of empirical data on the topic derived to date from plants and animals. We subsequently focus on the available data from fungi, which suggest that reproductive mode alters the rates and patterns of genome evolution in these organisms, e.g., protein evolution, mutation rate, codon usage, frequency of genome rearrangements and repetitive elements, and variation in chromosome size.

Transposable elements and small RNAs contribute to gene expression divergence between Arabidopsis thaliana and Arabidopsis lyrata

Transposable elements (TEs) are often the primary determinant of genome size differences among eukaryotes. In plants, the proliferation of TEs is countered through epigenetic silencing mechanisms that prevent mobility. Recent studies using the model plant Arabidopsis thaliana have revealed that methylated TE insertions are often associated with reduced expression of nearby genes, and these insertions may be subject to purifying selection due to this effect. Less is known about the genome-wide patterns of epigenetic silencing of TEs in other plant species. Here, we compare the 24-nt siRNA complement from A. thaliana and a closely related congener with a two- to threefold higher TE copy number, Arabidopsis lyrata. We show that TEs--particularly siRNA-targeted TEs--are associated with reduced gene expression within both species and also with gene expression differences between orthologs. In addition, A. lyrata TEs are targeted by a lower fraction of uniquely matching siRNAs, which are associated with more effective silencing of TE expression. Our results suggest that the efficacy of RNA-directed DNA methylation silencing is lower in A. lyrata, a finding that may shed light on the causes of differential TE proliferation among species.

Strategies for silencing and escape: the ancient struggle between transposable elements and their hosts

Over the past several years, there has been an explosion in our understanding of the mechanisms by which plant transposable elements (TEs) are epigenetically silenced and maintained in an inactive state over long periods of time. This highly efficient process results in vast numbers of inactive TEs; indeed, the majority of many plant genomes are composed of these quiescent elements. This observation has led to the rather static view that TEs represent an essentially inert portion of plant genomes. However, recent work has demonstrated that TE silencing is a highly dynamic process that often involves transcription of TEs at particular times and places during plant development. Plants appear to use transcripts from silenced TEs as an ongoing source of information concerning the mobile portion of the genome. In contrast to our understanding of silencing pathways, we know relatively little about the ways in which TEs evade silencing. However, vast differences in TE content between even closely related plant species suggest that they are often wildly successful at doing so. Here, we discuss TE activity in plants as the result of a constantly shifting balance between host strategies for TE silencing and TE strategies for escape and amplification.

T-lex: a program for fast and accurate assessment of transposable element presence using next-generation sequencing data

Transposable elements (TEs) are repetitive DNA sequences that are ubiquitous, extremely abundant and dynamic components of practically all genomes. Much effort has gone into annotation of TE copies in reference genomes. The sequencing cost reduction and the newly available next-generation sequencing (NGS) data from multiple strains within a species offer an unprecedented opportunity to study population genomics of TEs in a range of organisms. Here, we present a computational pipeline (T-lex) that uses NGS data to detect the presence/absence of annotated TE copies. T-lex can use data from a large number of strains and returns estimates of population frequencies of individual TE insertions in a reasonable time. We experimentally validated the accuracy of T-lex detecting presence or absence of 768 previously identified TE copies in two resequenced Drosophila melanogaster strains. Approximately 95% of the TE insertions were detected with 100% sensitivity and 97% specificity. We show that even at low levels of coverage T-lex produces accurate results for TE copies that it can identify reliably but that the rate of 'no data' calls increases as the coverage falls below 15×. T-lex is a broadly applicable and flexible tool that can be used in any genome provided the availability of the reference genome, individual TE copy annotation and NGS data.

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.

Mobile element scanning (ME-Scan) by targeted high-throughput sequencing

Background

Mobile elements (MEs) are diverse, common and dynamic inhabitants of nearly all genomes. ME transposition generates a steady stream of polymorphic genetic markers, deleterious and adaptive mutations, and substrates for further genomic rearrangements. Research on the impacts, population dynamics, and evolution of MEs is constrained by the difficulty of ascertaining rare polymorphic ME insertions that occur against a large background of pre-existing fixed elements and then genotyping them in many individuals.

Results

Here we present a novel method for identifying nearly all insertions of a ME subfamily in the whole genomes of multiple individuals and simultaneously genotyping (for presence or absence) those insertions that are variable in the population. We use ME-specific primers to construct DNA libraries that contain the junctions of all ME insertions of the subfamily, with their flanking genomic sequences, from many individuals. Individual-specific "index" sequences are designed into the oligonucleotide adapters used to construct the individual libraries. These libraries are then pooled and sequenced using a ME-specific sequencing primer. Mobile element insertion loci of the target subfamily are uniquely identified by their junction sequence, and all insertion junctions are linked to their individual libraries by the corresponding index sequence. To test this method's feasibility, we apply it to the human AluYb8 and AluYb9 subfamilies. In four individuals, we identified a total of 2,758 AluYb8 and AluYb9 insertions, including nearly all those that are present in the reference genome, as well as 487 that are not. Index counts show the sequenced products from each sample reflect the intended proportions to within 1%. At a sequencing depth of 355,000 paired reads per sample, the sensitivity and specificity of ME-Scan are both approximately 95%.

Conclusions

Mobile Element Scanning (ME-Scan) is an efficient method for quickly genotyping mobile element insertions with very high sensitivity and specificity. In light of recent improvements to high-throughput sequencing technology, it should be possible to employ ME-Scan to genotype insertions of almost any mobile element family in many individuals from any species.