A highly conserved, small LTR retrotransposon that preferentially targets genes in grass genomes

LTR retrotransposons are often the most abundant components of plant genomes and can impact gene and genome evolution. Most reported LTR retrotransposons are large elements (>4 kb) and are most often found in heterochromatic (gene poor) regions. We report the smallest LTR retrotransposon found to date, only 292 bp. The element is found in rice, maize, sorghum and other grass genomes, which indicates that it was present in the ancestor of grass species, at least 50-80 MYA. Estimated insertion times, comparisons between sequenced rice lines, and mRNA data indicate that this element may still be active in some genomes. Unlike other LTR retrotransposons, the small LTR retrotransposons (SMARTs) are distributed throughout the genomes and are often located within or near genes with insertion patterns similar to MITEs (miniature inverted repeat transposable elements). Our data suggests that insertions of SMARTs into or near genes can, in a few instances, alter both gene structures and gene expression. Further evidence for a role in regulating gene expression, SMART-specific small RNAs (sRNAs) were identified that may be involved in gene regulation. Thus, SMARTs may have played an important role in genome evolution and genic innovation and may provide a valuable tool for gene tagging systems in grass.     

Retrotransposon insertion polymorphisms in six rice genes and their evolutionary history

Retrotransposons are abundant in higher plant genomes. Although retrotransposons associated with plant genes have been identified, little is known about their evolutionary conservation at the level of species and subspecies. In the present study, we investigated the phylogenetic distribution of long terminal repeat (LTR) retrotransposon, long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) insertions in six genes in 95 cultivated and wild rice genotypes. These six genes are likely to be functional based on nonsynonymous (Ka) to synonymous (Ks) substitution ratios which were found to be significantly <1. Different conservation patterns of these retrotransposons in genes were observed in cultivated and wild rice species. Four out of seven retrotransposon insertions appear to predate the ancestral Oryza AA genome. Two of these insertions in genes 4 and 5 occurred early in the evolutionary history of Oryza. Two retrotransposon insertions in gene 1 arose after the divergence of Asian cultivated rice from its wild ancestor. Furthermore, the retrotransposon insertion in gene 3 appears to have occurred in the ancestral lineage leading to temperate japonicas. Conservation of retrotransposon insertions in genes in specific groups, species, and lineages might be related to their specific function.     

Characterization of genome-wide long terminal repeat retrotransposons provide insights into trait evolution of four grapevine species

Grapevine is one of the most economically important crops in the world. Although long terminal repeat (LTR) retrotransposons are thought to have played an important role in plants, its distribution in grapevine is not clear. Here, we identified genome-wide intact LTR retrotransposons in a total of six high-quality grapevine genomes from Vitis vinifera L., Vitis sylvestris C.C. Gmel., Vitis riparia Michx. and Vitis amurensis Rupr. with an average of 2938 per genome. Among them, the Copia superfamily (particularly for Ale) is a major component of the LTR retrotransposon in grapevine. Insertion time and copy number analysis revealed that the expansion of 70% LTR retrotransposons concentrating on approximately 2.5 Ma was able to drive genome size variation. Phylogenetic tree and syntenic analyses showed that most LTR retrotransposons in these genomes formed and evolved after species divergence. Furthermore, the function and expression of genes inserted by LTR retrotransposons in V. vinifera (Pinot noir) and V. riparia were explored. The length and expression of genes related to starch metabolism and quinone synthesis pathway in Pinot noir and environmental adaptation pathway in V. riparia were significantly affected by LTR retrotransposon insertion. The results improve the understanding of LTR retrotransposons in grapevine genomes and provide insights for its potential contribution to grapevine trait evolution.     

Structure and evolution of full-length LTR retrotransposons in rice genome

The long terminal repeat (LTR) retrotransposons are the most abundant class of transposable elements in plant genomes and play important roles in genome divergence and evolution. Their accumulation is the main factor influencing genome size increase in plants. Rice (Oryza sativa L.) is a model monocot and is the focus of much biological research due to its economic importance. We conducted a comprehensive survey of full-length LTR retrotransposons based on the completed genome of japonica rice variety Nipponbare (TIGR Release 5), with the newly published tool LTR-FINDER. The elements could be categorized into 29 structural domain categories (SDCs), and their total copy number identified was estimated at >6,000. Most of them were relatively young: more than 90% were less than 10 My. There existed a high level of activity among them as a whole at 0–1 Mya, but different categories possessed distinct amplification patterns. Most recently inserted elements were specific to the rice genome, while a few were conserved across species. This study provides new insights into the structure and evolutionary history of the full-length retroelements in the rice genome.     

High rate of chimeric gene origination by retroposition in plant genomes

Retroposition is widely found to play essential roles in origination of new mammalian and other animal genes. However, the scarcity of retrogenes in plants has led to the assumption that plant genomes rarely evolve new gene duplicates by retroposition, despite abundant retrotransposons in plants and a reported long terminal repeat (LTR) retrotransposon-mediated mechanism of retroposing cellular genes in maize (Zea mays). We show extensive retropositions in the rice (Oryza sativa) genome, with 1235 identified primary retrogenes. We identified 27 of these primary retrogenes within LTR retrotransposons, confirming a previously observed role of retroelements in generating plant retrogenes. Substitution analyses revealed that the vast majority are subject to negative selection, suggesting, along with expression data and evidence of age, that they are likely functional retrogenes. In addition, 42% of these retrosequences have recruited new exons from flanking regions, generating a large number of chimerical genes. We also identified young chimerical genes, suggesting that gene origination through retroposition is ongoing, with a rate an order of magnitude higher than the rate in primates. Finally, we observed that retropositions have followed an unexpected spatial pattern in which functional retrogenes avoid centromeric regions, while retropseudogenes are randomly distributed. These observations suggest that retroposition is an important mechanism that governs gene evolution in rice and other grass species.     

Young but not relatively old retrotransposons are preferentially located in gene‐rich euchromatic regions in tomato (Solanum lycopersicum) plants

Long terminal repeat (LTR) retrotransposons are the major DNA components of flowering plants. They are generally enriched in pericentromeric heterochromatin regions of their host genomes, which could result from the preferential insertion of LTR retrotransposons and the low effectiveness of purifying selection in these regions. To estimate the relative importance of the actions of these two factors on their distribution pattern, the LTR retrotransposons in Solanum lycopersicum (tomato) plants were characterized at the genome level, and then the distribution of young elements was compared with that of relatively old elements. The current data show that old elements are mainly located in recombination‐suppressed heterochromatin regions, and that young elements are preferentially located in the gene‐rich euchromatic regions. Further analysis showed a negative correlation between the insertion time of LTR retrotransposons and the recombination rate. The data also showed there to be more solo LTRs in genic regions than in intergenic regions or in regions close to genes. These observations indicate that, unlike in many other plant genomes, the current LTR retrotransposons in tomatoes have a tendency to be preferentially located into euchromatic regions, probably caused by their severe suppression of activities in heterochromatic regions. These elements are apt to be maintained in heterochromatin regions, probably as a consequence of the pericentromeric effect in tomatoes. These results also indicate that local recombination rates and intensities of purifying selection in different genomic regions are largely responsible for structural variation and non‐random distribution of LTR retrotransposons in tomato plants.     

Unearthing LTR Retrotransposon gag Genes Co-opted in the Deep Evolution of Eukaryotes

LTR retrotransposons comprise a major component of the genomes of eukaryotes. On occasion, retrotransposon genes can be recruited by their hosts for diverse functions, a process formally referred to as co-option. However, a comprehensive picture of LTR retrotransposon gag gene co-option in eukaryotes is still lacking, with several documented cases exclusively involving Ty3/Gypsy retrotransposons in animals. Here, we use a phylogenomic approach to systemically unearth co-option of retrotransposon gag genes above the family level of taxonomy in 2,011 eukaryotes, namely co-option occurring during the deep evolution of eukaryotes. We identify a total of 14 independent gag gene co-option events across more than 740 eukaryote families, eight of which have not been reported previously. Among these retrotransposon gag gene co-option events, nine, four, and one involve gag genes of Ty3/Gypsy, Ty1/Copia, and Bel-Pao retrotransposons, respectively. Seven, four, and three co-option events occurred in animals, plants, and fungi, respectively. Interestingly, two co-option events took place in the early evolution of angiosperms. Both selective pressure and gene expression analyses further support that these co-opted gag genes might perform diverse cellular functions in their hosts, and several co-opted gag genes might be subject to positive selection. Taken together, our results provide a comprehensive picture of LTR retrotransposon gag gene co-option events that occurred during the deep evolution of eukaryotes and suggest paucity of LTR retrotransposon gag gene co-option during the deep evolution of eukaryotes.     

Identification and characterization of genome-wide long terminal repeat retrotransposons provide an insight into elucidating the trait evolution of five Rhododendron species

• Rhododendron is well-known for its beauty and colourful corolla. Although some high-quality whole-genome sequencing of it has been completed, there are few studies on long terminal repeat (LTR) retrotransposons in Rhododendron, which limits our ability to elucidate the causes of genetic variations in Rhododendron species.
• Properties of the intact Rhododendron LTR retrotransposons were investigated at a genome-wide level. Based on available data, the high-quality genomes from five species, i.e. R. griersonianum, R. simsii, R. henanense subsp. lingbaoense, R. mucronatum var. ripense and R. ovatum, were selected as targets with good assembly continuity.
• A total of 17,936 intact LTR retrotransposons were identified; these belong to superfamilies Copia and Gypsy, with 17 clades. The insertion time of these transposons was later than 120 million years ago (Mya), and the outbreak period was concentrated more recently than 30 Mya. Phylogenetic analysis revealed that many LTR retrotransposons might originate from intraspecific duplication. Current evidence also suggests that most LTR retrotransposons were inserted in the interstitial part of genes in R. griersonianum, R. simsii, R. henanense, and R. ovatum, and the functions of the inserted genes mainly involve starch metabolism, proteolysis, etc. The effect of the LTR retrotransposon on gene expression depends on its insertion site and activation. Highly expressed LTR retrotransposons tend to be younger.
• The results herein improve our knowledge of LTR retrotransposons in Rhododendron genomes and facilitate further study of genetic variation and trait evolution in Rhododendron.

     

Temporal dynamics in the evolution of the sunflower genome as revealed by sequencing and annotation of three large genomic regions

Improved knowledge of genome composition, especially of its repetitive component, generates important informations in both theoretical and applied research. In this study, we provide the first insight into the local organization of the sunflower genome by sequencing and annotating 349,380 bp from 3 BAC clones, each including one single-copy gene. These analyses resulted in the identification of 11 putative gene sequences, 18 full-length LTR retrotransposons, 6 incomplete LTR retrotransposons, 2 non-autonomous LTR-retroelements (LINEs), 2 putative DNA transposons fragments and one putative helitron. Among LTR-retrotransposons, non-autonomous elements (the so-called LARDs), which do not carry any protein-encoding sequence, were discovered for the first time in the sunflower. The insertion time of intact retroelements was measured, based on sister LTRs divergence. All isolated elements were inserted relatively recently, especially those belonging to the Gypsy superfamily. Retrotransposon families related to those identified in the BAC clones are present also in other species of Helianthus, both annual and perennial, and even in other Asteraceae. In one of the three BAC clones, we found five copies of a lipid transfer protein (LTP) encoding gene within less than 100,000 bp, four of which are potentially functional. Two of these are interrupted by LTR retrotransposons, in the intron and in the coding sequence, respectively. The divergence between sister LTRs of the retrotransposons inserted within the genes indicates that LTP gene duplication started earlier than 1.749 MYRS ago. On the whole, the results reported in this study confirm that the sunflower is an excellent system to study transposons dynamics and evolution.     

Analysis of genes associated with retrotransposons in the rice genome

Retrotransposons comprise a significant fraction of the rice genome. Despite their prevalence, the effects of retrotransposon insertions are not well understood, especially with regard to how they affect the expression of genes. In this study, we identified one-sixth of rice genes as being associated with retrotransposons, with insertions either in the gene itself or within its putative promoter region. Among genes with insertions in the promoter region, the likelihood of the gene being expressed was shown to be directly proportional to the distance of the retrotransposon from the translation start site. In addition, retrotransposon insertions in the transcribed region of the gene were found to be positively correlated with the presence of alternative splicing forms. Furthermore, preferential association of retrotransposon insertions with genes in several functional classes was identified. Some of the retrotransposons that are part of full-length cDNA (fl-cDNA) contribute splice sites and give rise to novel exons. Several interesting trends concerning the effects of retrotransposon insertions on gene expression were identified. Taken together, our data suggests that retrotransposon association with genes have a role in gene regulation. The data presented in this study provides a foundation for experimental studies to determine the role of retrotransposons in gene regulation.     

Sequence-based alignment of sorghum chromosome 3 and rice chromosome 1 reveals extensive conservation of gene order and one major chromosomal rearrangement

The completed rice genome sequence will accelerate progress on the identification and functional classification of biologically important genes and serve as an invaluable resource for the comparative analysis of grass genomes. In this study, methods were developed for sequence-based alignment of sorghum and rice chromosomes and for refining the sorghum genetic/physical map based on the rice genome sequence. A framework of 135 BAC contigs spanning approximately 33 Mbp was anchored to sorghum chromosome 3. A limited number of sequences were collected from 118 of the BACs and subjected to BLASTX analysis to identify putative genes and BLASTN analysis to identify sequence matches to the rice genome. Extensive conservation of gene content and order between sorghum chromosome 3 and the homeologous rice chromosome 1 was observed. One large-scale rearrangement was detected involving the inversion of an approximately 59 cM block of the short arm of sorghum chromosome 3. Several small-scale changes in gene collinearity were detected, indicating that single genes and/or small clusters of genes have moved since the divergence of sorghum and rice. Additionally, the alignment of the sorghum physical map to the rice genome sequence allowed sequence-assisted assembly of an approximately 1.6 Mbp sorghum BAC contig. This streamlined approach to high-resolution genome alignment and map building will yield important information about the relationships between rice and sorghum genes and genomic segments and ultimately enhance our understanding of cereal genome structure and evolution.     

Natural selection maintains the transcribed LTR retrotransposons in Nosema bombycis

Eight intact LTR retrotransposons (Nbrl-Nbr8) have been previously characterized from the genome of Nosema bombycis, a eu-karyotic parasite with a compact and reduced genome. Here we describe six novel transcribed Nbr elements (Nbr9-Nbrl4) identified through either cDNA library or RT-PCR. Like previously determined ones, all of them belong to the Ty3/Gypsy superfamily. Retro-transposon diversity and incomplete domains with insertions {Nbrll), deletions (Nbrll) and in-frame stop codons in coding regions (Nbr9) were detected, suggesting that both defective and loss events of LTR retrotransposon have happened in N. bombycis genome.Analysis of selection showed that strong purifying selection acts on all elements except Nbrll. This implies that selective pressure keeps both these Nbrs and their functions in genome. Interestingly, Nbrll is under positive selection and some positively selected codons were identified, indicating that new functionality might have evolved in the Nbrll retrotransposon. Unlike other transposable elements, Nbrll has integrated into a conserved syntenic block and probably resulted in the inversion of both flanking regions. This demonstrates that transposable element is an important factor for the reshuffling and evolution of their host genomes, and may be maintained under natural selection.     

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).     

Swine liver methylomes of Berkshire,Duroc and Landrace breeds by MeDIPS

Using a methyl‐DNA immunoprecipitation technique in combination with next‐generation deep sequencing, we conducted comprehensive DNA methylation profiling of liver genomes from three pig breeds: Berkshire, Duroc and Landrace. The profiles revealed that the distribution patterns of methylation signals along the genome are conserved among the three pig breeds. Specifically, many signals in coding genes were found in introns, and most signals in the repetitive elements were identified in non‐long terminal repeat (LTR) retrotransposons such as long and short interspersed repetitive elements, implying a significant association with alternative splicing and expression of retrotransposable elements respectively. Differentially methylated regions among the three pig breeds were identified in the non‐LTR retrotransposons, suggesting that they may lead to differential retrotransposable element activity. Altogether, this study provides advanced swine methylome data and valuable resources for understanding the function of DNA methylation in the evolutionary divergence of different pig breeds.