The fate of transposable elements in asexual populations

Sexual reproduction and recombination are important for maintaining a stable copy number of transposable elements (TEs). In sexual populations, elements can be contained by purifying selection against host carriers with higher element copy numbers; however, in the absence of sex and recombination, asexual populations could be driven to extinction by an unchecked proliferation of TEs. Here we provide a theoretical framework for analyzing TE dynamics under asexual reproduction. Analytic results show that, in an infinite asexual population, an equilibrium in copy number is achieved if no element excision is possible, but that all TEs are eliminated if there is some excision. In a finite population, computer simulations demonstrate that small populations are driven to extinction by a Muller's ratchet-like process of element accumulation, but that large populations can be cured of vertically transmitted TEs, even with excision rates well below transposition rates. These results may have important consequences for newly arisen asexual lineages and may account for the lack of deleterious retrotransposons in the putatively ancient asexual bdelloid rotifers.     

Analysis of Transposable Elements in the Genome of Asparagus officinalis from High Coverage Sequence Data

Asparagus officinalis is an economically and nutritionally important vegetable crop that is widely cultivated and is used as a model dioecious species to study plant sex determination and sex chromosome evolution. To improve our understanding of its genome composition, especially with respect to transposable elements (TEs), which make up the majority of the genome, we performed Illumina HiSeq2000 sequencing of both male and female asparagus genomes followed by bioinformatics analysis. We generated 17 Gb of sequence (12×coverage) and assembled them into 163,406 scaffolds with a total cumulated length of 400 Mbp, which represent about 30% of asparagus genome. Overall, TEs masked about 53% of the A. officinalis assembly. Majority of the identified TEs belonged to LTR retrotransposons, which constitute about 28% of genomic DNA, with Ty1/copia elements being more diverse and accumulated to higher copy numbers than Ty3/gypsy. Compared with LTR retrotransposons, non-LTR retrotransposons and DNA transposons were relatively rare. In addition, comparison of the abundance of the TE groups between male and female genomes showed that the overall TE composition was highly similar, with only slight differences in the abundance of several TE groups, which is consistent with the relatively recent origin of asparagus sex chromosomes. This study greatly improves our knowledge of the repetitive sequence construction of asparagus, which facilitates the identification of TEs responsible for the early evolution of plant sex chromosomes and is helpful for further studies on this dioecious plant.     

Asexual Experimental Evolution of Yeast Does Not Curtail Transposable Elements

Compared with asexual reproduction, sex facilitates the transmission of transposable elements (TEs) from one genome to another, but boosts the efficacy of selection against deleterious TEs. Thus, theoretically, it is unclear whether sex has a positive net effect on TE’s proliferation. An empirical study concluded that sex is at the root of TE’s evolutionary success because the yeast TE load was found to decrease rapidly in approximately 1,000 generations of asexual but not sexual experimental evolution. However, this finding contradicts the maintenance of TEs in natural yeast populations where sexual reproduction occurs extremely infrequently. Here, we show that the purported TE load reduction during asexual experimental evolution is likely an artifact of low genomic sequencing coverages. We observe stable TE loads in both sexual and asexual experimental evolution from multiple yeast data sets with sufficient coverages. To understand the evolutionary dynamics of yeast TEs, we turn to asexual mutation accumulation lines that have been under virtually no selection. We find that both TE transposition and excision rates per generation, but not their difference, tend to be higher in environments where yeast grows more slowly. However, the transposition rate is not significantly higher than the excision rate and the variance of the TE number among natural strains is close to its neutral expectation, suggesting that selection against TEs is at best weak in yeast. We conclude that the yeast TE load is maintained largely by a transposition–excision balance and that the influence of sex remains unclear.     

Transposable elements in clonal lineages: lethal hangover from sex

Long-term coevolution of transposable elements (TEs) in sexual hosts leads to evolution of extremely active and dangerous mutagens kept in tenuous check by host-derived mechanisms and via natural selection against TE-rich genomes. To the extent that sexual reproduction and recombination are important in maintaining a stable TE copy number and a tolerable mutation load, the switch to clonality from sexual reproduction can be extremely damaging and, generally, should lead to clonal lineage extinction. Surprisingly however, the loss of powerful selective mechanisms constraining TEs can be beneficial in the short-term by immediately eliminating selective load and possibly promoting the early success of clonal lineages. The clonal lineages that do survive in the long-term must find a way to eliminate or domesticate TEs. Indeed bdelloid rotifers, which are ancient asexuals, do appear to have lost most of the otherwise wide-spread TEs and might have domesticated others. The path to this TE-free haven is anything but clear at the moment. We have considered a novel scenario of instantaneous inactivation of TEs by starting off with a genome carrying repressive host alleles for all TEs in the genome. We show that such a scenario appears plausible and provide some limited empirical evidence in its support.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79 , 33–41.     

Evolutionary history of mammalian transposons determined by genome-wide defragmentation

The constant bombardment of mammalian genomes by transposable elements (TEs) has resulted in TEs comprising at least 45% of the human genome. Because of their great age and abundance, TEs are important in comparative phylogenomics. However, estimates of TE age were previously based on divergence from derived consensus sequences or phylogenetic analysis, which can be unreliable, especially for older more diverged elements. Therefore, a novel genome-wide analysis of TE organization and fragmentation was performed to estimate TE age independently of sequence composition and divergence or the assumption of a constant molecular clock. Analysis of TEs in the human genome revealed approximately 600,000 examples where TEs have transposed into and fragmented other TEs, covering >40% of all TEs or approximately 542 Mbp of genomic sequence. The relative age of these TEs over evolutionary time is implicit in their organization, because newer TEs have necessarily transposed into older TEs that were already present. A matrix of the number of times that each TE has transposed into every other TE was constructed, and a novel objective function was developed that derived the chronological order and relative ages of human TEs spanning >100 million years. This method has been used to infer the relative ages across all four major TE classes, including the oldest, most diverged elements. Analysis of DNA transposons over the history of the human genome has revealed the early activity of some MER2 transposons, and the relatively recent activity of MER1 transposons during primate lineages. The TEs from six additional mammalian genomes were defragmented and analyzed. Pairwise comparison of the independent chronological orders of TEs in these mammalian genomes revealed species phylogeny, the fact that transposons shared between genomes are older than species-specific transposons, and a subset of TEs that were potentially active during periods of speciation.     

Wake up of transposable elements following Drosophila simulans worldwide colonization.

Transposable elements (TEs) make up around 10%-15% of the Drosophila melanogaster genome, but its sibling species Drosophila simulans carries only one third as many such repeat sequences. We do not, however, have an overall view of copy numbers of the various classes of TEs (long terminal repeat [LTR] retrotransposons, non-LTR retrotransposons, and transposons) in genomes of natural populations of both species. We analyzed 34 elements in individuals from various natural populations of these species. We show that D. melanogaster has higher average chromosomal insertion site numbers per genome than D. simulans for all TEs except five. The LTR retrotransposons gypsy, ZAM, and 1731 and the transposon bari-1 present similar low copy numbers in both species. The transposon hobo has a large number of insertion sites, with significantly more sites in D. simulans. High variation between populations in number of insertion sites of some elements of D. simulans suggests that these elements can invade the genome of the entire species starting from a local population. We propose that TEs in the D. simulans genome are being awakened and amplified as they had been a long time ago in D. melanogaster.     

Fundamentally different repetitive element composition of sex chromosomes in Rumex acetosa

Background and AimsDioecious species with well-established sex chromosomes are rare in the plant kingdom. Most sex chromosomes increase in size but no comprehensive analysis of the kind of sequences that drive this expansion has been presented. Here we analyse sex chromosome structure in common sorrel (Rumex acetosa), a dioecious plant with XY₁Y₂ sex determination, and we provide the first chromosome-specific repeatome analysis for a plant species possessing sex chromosomes.MethodsWe flow-sorted and separately sequenced sex chromosomes and autosomes in R. acetosa using the two-dimensional fluorescence in situ hybridization in suspension (FISHIS) method and Illumina sequencing. We identified and quantified individual repeats using RepeatExplorer, Tandem Repeat Finder and the Tandem Repeats Analysis Program. We employed fluorescence in situ hybridization (FISH) to analyse the chromosomal localization of satellites and transposons.Key ResultsWe identified a number of novel satellites, which have, in a fashion similar to previously known satellites, significantly expanded on the Y chromosome but not as much on the X or on autosomes. Additionally, the size increase of Y chromosomes is caused by non-long terminal repeat (LTR) and LTR retrotransposons, while only the latter contribute to the enlargement of the X chromosome. However, the X chromosome is populated by different LTR retrotransposon lineages than those on Y chromosomes.ConclusionsThe X and Y chromosomes have significantly diverged in terms of repeat composition. The lack of recombination probably contributed to the expansion of diverse satellites and microsatellites and faster fixation of newly inserted transposable elements (TEs) on the Y chromosomes. In addition, the X and Y chromosomes, despite similar total counts of TEs, differ significantly in the representation of individual TE lineages, which indicates that transposons proliferate preferentially in either the paternal or the maternal lineage.     

The Association Between Breeding System and Transposable Element Dynamics in <Emphasis Type="Italic">Daphnia Pulex</Emphasis>

Transposable elements (TEs) are major sources of genetic variation, and mating systems are believed to play a significant role in their dynamics. For example, insertion number is expected to be higher in sexual than in asexual organisms due to the inability of TEs to colonize new genomes in the absence of sex. The goal of this study was to determine the impact of the loss of sexual reproduction on TE load. Daphnia pulex has two reproductive modes, obligate and cyclical parthenogenesis, which differ with respect to the production of diapausing eggs. Cyclical parthenogens produce them meiotically, while obligate parthenogens produce them clonally. Pokey is a TTAA-specific DNA transposon, and is a stable component of Daphnia genomes. We used a PCR-based approach, TE-Display, to estimate the number of Pokey insertions in 22 cyclic and 22 obligate isolates of D. pulex. As expected, the copy number of Pokey insertions is significantly higher in cyclic than in obligate isolates. However, the distribution of elements among isolates within each breeding system is similar, which is congruent with the recent establishment of obligate lineages from a cyclic ancestor. We also assayed 46 isolates from eight cyclic populations and found that very few Pokey insertions were observed in more than one isolate, suggesting that Pokey has been active recently. Sequencing of PCR products from the TE-Display analysis shows that Pokey inserts into both coding and noncoding regions of the genome. However, there is no obvious similarity among sequences downstream of the TTAA Pokey insertion site.     

No evidence that sex and transposable elements drive genome size variation in evening primroses

Genome size varies dramatically across species, but despite an abundance of attention there is little agreement on the relative contributions of selective and neutral processes in governing this variation. The rate of sex can potentially play an important role in genome size evolution because of its effect on the efficacy of selection and transmission of transposable elements (TEs). Here, we used a phylogenetic comparative approach and whole genome sequencing to investigate the contribution of sex and TE content to genome size variation in the evening primrose (Oenothera) genus. We determined genome size using flow cytometry for 30 species that vary in genetic system and find that variation in sexual/asexual reproduction cannot explain the almost twofold variation in genome size. Moreover, using whole genome sequences of three species of varying genome sizes and reproductive system, we found that genome size was not associated with TE abundance; instead the larger genomes had a higher abundance of simple sequence repeats. Although it has long been clear that sexual reproduction may affect various aspects of genome evolution in general and TE evolution in particular, it does not appear to have played a major role in genome size evolution in the evening primroses.     

Asexual evolution: do intragenomic parasites maintain sex?

Resolving the paradox of sex, with its twofold cost to genic transmission, remains one of the major unresolved questions in evolutionary biology. Counting this genetic cost has now gone genomic. In this issue of Molecular Ecology, Kraaijeveld et al. (2012) describe the first genome‐scale comparative study of related sexual and asexual animal lineages, to test the hypothesis that asexuals bear heavier loads of deleterious transposable elements. A much higher density of such parasites might be expected, due to the inability of asexual lineages to purge transposons via mechanisms exclusive to sexual reproduction. They find that the answer is yes—and no—depending upon the family of transposons considered. Like many such advances in testing theory, more questions are raised by this study than answered, but a door has been opened to molecular evolutionary analyses of how responses to selection from intragenomic parasites might mediate the costs of sex.     

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.     

Transposons but not retrotransposons are located preferentially in regions of high recombination rate in Caenorhabditis elegans

We analyzed the distribution of transposable elements (TEs: transposons, LTR retrotransposons, and non-LTR retrotransposons) in the chromosomes of the nematode Caenorhabditis elegans. The density of transposons (DNA-based elements) along the chromosomes was found to be positively correlated with recombination rate, but this relationship was not observed for LTR or non-LTR retrotransposons (RNA-based elements). Gene (coding region) density is higher in regions of low recombination rate. However, the lower TE density in these regions is not due to the counterselection of TE insertions within exons since the same positive correlation between TE density and recombination rate was found in noncoding regions (both in introns and intergenic DNA). These data are not compatible with a global model of selection acting against TE insertions, for which an accumulation of elements in regions of reduced recombination is expected. We also found no evidence for a stronger selection against TE insertions on the X chromosome compared to the autosomes. The difference in distribution of the DNA and RNA-based elements along the chromosomes in relation to recombination rate can be explained by differences in the transposition processes.     

TEnest: automated chronological annotation and visualization of nested plant transposable elements

Organisms with a high density of transposable elements (TEs) exhibit nesting, with subsequent repeats found inside previously inserted elements. Nesting splits the sequence structure of TEs and makes annotation of repetitive areas challenging. We present TEnest, a repeat identification and display tool made specifically for highly repetitive genomes. TEnest identifies repetitive sequences and reconstructs separated sections to provide full-length repeats and, for long-terminal repeat (LTR) retrotransposons, calculates age since insertion based on LTR divergence. TEnest provides a chronological insertion display to give an accurate visual representation of TE integration history showing timeline, location, and families of each TE identified, thus creating a framework from which evolutionary comparisons can be made among various regions of the genome. A database of repeats has been developed for maize (Zea mays), rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare) to illustrate the potential of TEnest software. All currently finished maize bacterial artificial chromosomes totaling 29.3 Mb were analyzed with TEnest to provide a characterization of the repeat insertions. Sixty-seven percent of the maize genome was found to be made up of TEs; of these, 95% are LTR retrotransposons. The rate of solo LTR formation is shown to be dissimilar across retrotransposon families. Phylogenetic analysis of TE families reveals specific events of extreme TE proliferation, which may explain the high quantities of certain TE families found throughout the maize genome. The TEnest software package is available for use on PlantGDB under the tools section (http://www.plantgdb.org/prj/TE_nest/TE_nest.html); the source code is available from (http://wiselab.org).     

Retrotransposon Sequence Variation in Four Asexual Plant Species

Transposable elements (TEs) can be viewed as genetic parasites that persist in populations due to their capacity for increase in copy number and the inefficacy of selection against them. A corollary of this hypothesis is that TEs are more likely to spread within sexual populations and be eliminated or inactivated within asexual populations. While previous work with animals has shown that asexual taxa may contain less TE diversity than sexual taxa, comparable work with plants has been lacking. Here we report the results of a study of Ty1/copia, Ty3/gypsy, and LINE-like retroelement diversity in four asexual plant species. Retroelement-like sequences, with a high degree of conservation both within and between species, were isolated from all four species. The sequences correspond to several previously annotated retroelement subfamilies. They also exhibit a pattern of nucleotide substitution characterized by an excess of synonymous substitutions, suggestive of a history of purifying selection. These findings were compared with retroelement sequence evolution in sexual plant taxa. One likely explanation for the discovery of conserved TE sequences in the genomes of these asexual taxa is simply that asexuality within these taxa evolved relatively recently, such that the loss and breakdown of TEs is not yet detectable through analysis of sequence diversity. This explanation is examined by conducting stochastic simulation of TE evolution and by using published information to infer rough estimates of the ages of asexual taxa. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Dmitri Petrov     

白菜和甘蓝基因组转座子表达及其对基因调控的潜在影响

转座子或转座元件是大多数真核生物基因组的主要组成成分。甘蓝(Brassica oleracea)基因组比白菜(B. rapa)大主要是转座子的扩增差异造成的。然而, 这两个芸薹属近缘物种转座子表达水平以及对基因的调控和功能的影响目前还不清楚。文章对白菜和甘蓝叶、根、茎3个器官的转录组数据进行了初步分析。结果显示, 转座子的表达量很低, 转录组reads中有1%来自转座子的转录本; 转座子的表达存在器官差异, 且不同类别和家族的转座子表达量相差很大, 相同类别和同一家族的转座子在白菜和甘蓝基因组中的表达活性也不相同。进一步鉴定到转录读出的LTR反转座子, 其与下游基因距离小于2 kb的有41个, 小于100 bp的有9个, 这些LTR的转录读出很可能通过正义或反义的转录本激活或干扰下游基因的表达。同时, 具有转录读出的intact LTR比solo LTR具有更强的读出活性。通过深入分析转座子的插入位点发现, 白菜基因组中转座子插入基因内部的频率比甘蓝基因组中的高; 与反转座子相比, DNA转座子更偏向于插入或保留在基因的内含子当中。这些结果为认识转座子对其他蛋白编码基因的影响提供了基础。     

Sequence composition, organization, and evolution of the core Triticeae genome

We investigated the composition and the basis of genome expansion in the core Triticeae genome using Aegilops tauschii, the D-genome donor of bread wheat. We sequenced an unfiltered genomic shotgun (trs) and a methylation-filtration (tmf) library of A. tauschii, and analyzed wheat expressed sequence tags (ESTs) to estimate the expression of genes and transposable elements (TEs). The sampled D-genome sequences consisted of 91.6% repetitive elements, 2.5% known genes, and 5.9% low-copy sequences of unknown function. TEs constituted 68.2% of the D-genome compared with 50% in maize and 14% in rice. The DNA transposons constituted 13% of the D-genome compared with 2% in maize. TEs were methylated unevenly within and among elements and families, and most were transcribed which contributed to genome expansion in the core Triticeae genome. The copy number of a majority of repeat families increased gradually following polyploidization. Certain TE families occupied discrete chromosome territories. Nested insertions and illegitimate recombination occurred extensively between the TE families, and a majority of the TEs contained internal deletions. The GC content varied significantly among the three sequence sets examined ranging from 42% in tmf to 46% in trs and 52% in the EST. Based on enrichment of genic sequences, methylation-filtration offers one option, although not as efficient as in maize, for isolating gene-rich regions from the large genome of wheat.     

Copy Number Variation in Transcriptionally Active Regions of Sexual and Apomictic Boechera Demonstrates Independently Derived Apomictic Lineages

In asexual (apomictic) plants, the absence of meiosis and sex is expected to lead to mutation accumulation. To compare mutation accumulation in the transcribed genomic regions of sexual and apomictic plants, we performed a double-validated analysis of copy number variation (CNV) on 10 biological replicates each of diploid sexual and diploid apomictic Boechera, using a high-density (>700 K) custom microarray. The Boechera genome demonstrated higher levels of depleted CNV, compared with enriched CNV, irrespective of reproductive mode. Genome-wide patterns of CNV revealed four divergent lineages, three of which contain both sexual and apomictic genotypes. Hence genome-wide CNV reflects at least three independent origins (i.e., expression) of apomixis from different sexual genetic backgrounds. CNV distributions for different families of transposable elements were lineage specific, and the enrichment of LINE/L1 and long term repeat/Copia elements in lineage 3 apomicts is consistent with sex and meiosis being mechanisms for purging genomic parasites. We hypothesize that significant overrepresentation of specific gene ontology classes (e.g., pollen–pistil interaction) in apomicts implies that gene enrichment could be an adaptive mechanism for genome stability in diploid apomicts by providing a polyploid-like system for buffering the effects of deleterious mutations.