Identification of novel LTR retrotransposons in the genome of Aedes aegypti

We have detected seventy-six novel LTR retrotransposons in the genome of the mosquito Aedes aegypti by a genome wide analysis using the LTR_STRUC program. We have performed a phylogenetic classification of these novel elements and a distribution analysis in the genome of A. aegypti. These mobile elements belong either to the Ty3/gypsy or to the Bel family of retrotransposons and were not annotated in the mosquito LTR retrotransposon database (TEfam). We have found that  1.8% of the genome is occupied by these newly detected retrotransposons that are distributed predominantly in intergenic genomic sequences and introns. The potential role of retrotransposon insertions linked to host genes is described and discussed. We show that a retrotransposon family belonging to the Osvaldo lineage has peculiar structural features, and its presence is likely to be restricted to the A. aegypti and to the Culex pipiens quinquefasciatus genomes. Furthermore we show that the ninja-like group of elements lacks the Primer Binding Site (PBS) sequence necessary for the replication of retrotransposons. These results integrate the knowledge on the complicate genomic structure of an important disease vector.     

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.     

Diversity of LTR retrotransposons and their role in genome reorganization

Current views of retrotransposons possessing long terminal repeats (LTRs) are described. The existing classification and element types isolated by genome organization are considered. Experimental data are summarized to demonstrate that the replicative cycle of a retrotransposon is not restricted to a single cell and that LTR retrotransposons are transferred between somatic cells with a rate comparable with the element transposition rate within the genome of one cell. The major mechanisms mediating the role of LTR retrotransposons in reorganization of the genome are considered with regard to the strategies of their horizontal and vertical transfer.     

Diversity of LTR retrotransposons and their role in genome reorganization

Current views of retrotransposons possessing long terminal repeats (LTRs) are described. The existing classification and element types isolated by genome organization are considered. Experimental data are summarized to demonstrate that the replicative cycle of a retrotransposon is not restricted to a single cell and that LTR retrotransposons are transferred between somatic cells with a rate comparable with the element transposition rate within the genome of one cell. The major mechanisms mediating the role of LTR retrotransposons in reorganization of the genome are considered with regard to the strategies of their horizontal and vertical transfer.__________Translated from Genetika, Vol. 41, No. 4, 2005, pp. 542–548.Original Russian Text Copyright © 2005 by Syomin, Ilyin.     

Evolutionary dynamics of the Ty3/gypsy LTR retrotransposons in the genome of Anopheles gambiae

Ty3/gypsy elements represent one of the most abundant and diverse LTR-retrotransposon (LTRr) groups in the Anopheles gambiae genome, but their evolutionary dynamics have not been explored in detail. Here, we conduct an in silico analysis of the distribution and abundance of the full complement of 1045 copies in the updated AgamP3 assembly. Chromosomal distribution of Ty3/gypsy elements is inversely related to arm length, with densities being greatest on the X, and greater on the short versus long arms of both autosomes. Taking into account the different heterochromatic and euchromatic compartments of the genome, our data suggest that the relative abundance of Ty3/gypsy LTRrs along each chromosome arm is determined mainly by the different proportions of heterochromatin, particularly pericentric heterochromatin, relative to total arm length. Additionally, the breakpoint regions of chromosomal inversion 2La appears to be a haven for LTRrs. These elements are underrepresented more than 7-fold in euchromatin, where 33% of the Ty3/gypsy copies are associated with genes. The euchromatin on chromosome 3R shows a faster turnover rate of Ty3/gypsy elements, characterized by a deficit of proviral sequences and the lowest average sequence divergence of any autosomal region analyzed in this study. This probably reflects a principal role of purifying selection against insertion for the preservation of longer conserved syntenyc blocks with adaptive importance located in 3R. Although some Ty3/gypsy LTRrs show evidence of recent activity, an important fraction are inactive remnants of relatively ancient insertions apparently subject to genetic drift. Consistent with these computational predictions, an analysis of the occupancy rate of putatively older insertions in natural populations suggested that the degenerate copies have been fixed across the species range in this mosquito, and also are shared with the sibling species Anopheles arabiensis.     

A computational study of the dynamics of LTR retrotransposons in the Populus trichocarpa genome

Retrotransposons are an ubiquitous component of plant genomes, especially abundant in species with large genomes. Populus trichocarpa has a relatively small genome, which was entirely sequenced; however, studies focused on poplar retrotransposons dynamics are rare. With the aim to study the retrotransposon component of the poplar genome, we have scanned the complete genome sequence searching full-length long-terminal repeat (LTR) retrotransposons, i.e., characterised by two long terminal repeats at the 5′ and 3′ ends. A computational approach based on detection of conserved structural features, on building multiple alignments, and on similarity searches was used to identify 1,479 putative full-length LTR retrotransposons. Ty1-copia elements were more numerous than Ty3-gypsy. However, many LTR retroelements were not assigned to any superfamily because lacking of diagnostic features and non-autonomous. LTR retrotransposon remnants were by far more numerous than full-length elements, indicating that during the evolution of poplar, large amplification of these elements was followed by DNA loss. Within superfamilies, Ty3-gypsy families are made of more members than Ty1-copia ones. Retrotransposition occurred with increasing frequency following the separation of Populus sections, with different waves of retrotransposition activity between Ty3-gypsy and Ty1-copia elements. Recently inserted elements appear more frequently expressed than older ones. Finally, different levels of activity of retrotransposons were observed according to their position and their density in the linkage groups. On the whole, the results support the view of retrotransposons as a community of different organisms in the genome, whose activity (both retrotransposition and DNA loss) has heavily impacted and probably continues to impact poplar genome structure and size.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae

The recent availability of the genome of Anopheles gambiae offers an extraordinary opportunity for comparative studies of the diversity of transposable elements (TEs) and their evolutionary dynamics between two related species, taking advantage of the existing information from Drosophila melanogaster. To this goal, we screened the genome of A. gambiae for elements belonging to the Ty3/gypsy group of long-terminal repeat (LTR) retrotransposons. The A. gambiae genome displays a rich diversity of LTR retrotransposons, clearly greater than D. melanogaster. We have characterized in detail 63 families, belonging to five of the nine main lineages of the Ty3/gypsy group. The Mag lineage is the most diverse and abundant, with more than 30 families. In sharp contrast with this finding, a single family belonging to this lineage has been found in D. melanogaster, here reported for the first time in the literature, most probably consisting of old inactive elements. The CsRn1 lineage is also abundant in A. gambiae but almost absent from D. melanogaster. Conversely, the Osvaldo lineage has been detected in Drosophila but not in Anopheles. Comparison of structural characteristics of different families led to the identification of several lineage-specific features such as the primer-binding site (PBS), the gag-pol translational recoding signal (TRS), which is extraordinarily diverse within the Ty3/gypsy retrotransposons of A. gambiae, or the presence/absence of specific amino acid motifs. Interestingly, some of these characteristics, although in general well conserved within lineages, may have evolved independently in particular branches of the phylogenetic tree. We also show evidence of recent activity for around 75% of the families. Nevertheless, almost all families contain a high proportion of degenerate members and solitary LTRs (solo LTRs), indicative of a lower turnover rate of retrotransposons belonging to the Ty3/gypsy group in A. gambiae than in D. melanogaster. Finally, we have detected significant overrepresentations of insertions on the X chromosome versus autosomes and of putatively active insertions on euchromatin versus heterochromatin.     

Correlated evolution of LTR retrotransposons and genome size in the genus <Emphasis Type="Italic">eleocharis</Emphasis>

Background  

Transposable elements (TEs) are considered to be an important source of genome size variation and genetic and phenotypic plasticity in eukaryotes. Most of our knowledge about TEs comes from large genomic projects and studies focused on model organisms. However, TE dynamics among related taxa from natural populations and the role of TEs at the species or supra-species level, where genome size and karyotype evolution are modulated in concert with polyploidy and chromosomal rearrangements, remain poorly understood. We focused on the holokinetic genus Eleocharis (Cyperaceae), which displays large variation in genome size and the occurrence of polyploidy and agmatoploidy/symploidy. We analyzed and quantified the long terminal repeat (LTR) retrotransposons Ty1-copia and Ty3-gypsy in relation to changes in both genome size and karyotype in Eleocharis. We also examined how this relationship is reflected in the phylogeny of Eleocharis.     

B chromosomes and genome size in flowering plants.

B chromosomes are extra chromosomes found in some, but not all, individuals within a species, often maintained by giving themselves an advantage in transmission, i.e. they drive. Here we show that the presence of B chromosomes correlates to and varies strongly and positively with total genome size (excluding the Bs and corrected for ploidy) both at a global level and via a comparison of independent taxonomic contrasts. B chromosomes are largely absent from species with small genomes; however, species with large genomes are studied more frequently than species with small genomes and Bs are more likely to be reported in well-studied species. We controlled for intensity of study using logistic regression. This regression analysis also included effects of degree of outbreeding, which is positively associated with Bs and genome size, and chromosome number, which is negatively associated with Bs and genome size, as well as variable ploidy (more than one ploidy level in a species). Genome size, breeding system and chromosome number all contribute independently to the distribution of B chromosomes, while variable ploidy does not have a significant effect. The genome size correlates are consistent with reduced selection against extra DNA in species with large genomes and with increased generation of B sequences from large A genomes.     

LTR_STRUC: a novel search and identification program for LTR retrotransposons

MOTIVATION: Long terminal repeat (LTR) retrotransposons constitute a substantial fraction of most eukaryotic genomes and are believed to have a significant impact on genome structure and function. Conventional methods used to search for LTR retrotransposons in genome databases are labor intensive. We present an efficient, reliable and automated method to identify and analyze members of this important class of transposable elements. RESULTS: We have developed a new data-mining program, LTR_STRUC (LTR retrotransposon structure program) which identifies and automatically analyzes LTR retrotransposons in genome databases by searching for structural features characteristic of such elements. LTR_STRUC has significant advantages over conventional search methods in the case of LTR retrotransposon families having low sequence homology to known queries or families with atypical structure (e.g. non-autonomous elements lacking canonical retroviral ORFs) and is thus a discovery tool that complements established methods. LTR_STRUC finds LTR retrotransposons using an algorithm that encompasses a number of tasks that would otherwise have to be initiated individually by the user. For each LTR retrotransposon found, LTR_STRUC automatically generates an analysis of a variety of structural features of biological interest. AVAILABILITY: The LTR_STRUC program is currently available as a console application free of charge to academic users from the authors.     

Mechanisms of recent genome size variation in flowering plants

BACKGROUND AND AIMS: Plant nuclear genomes vary tremendously in DNA content, mostly due to differences in ancestral ploidy and variation in the degree of transposon amplification. These processes can increase genome size, but little is known about mechanisms of genome shrinkage and the degree to which these can attenuate or reverse genome expansion. This research focuses on characterizing DNA removal from the rice and Arabidopsis genomes, and discusses whether loss of DNA has effectively competed with amplification in these species. METHODS: Retrotransposons were analyzed for sequence variation within several element families in rice and Arabidopsis. Nucleotide sequence changes in the two termini of individual retrotransposons were used to date their time of insertion. KEY RESULTS: An accumulation of small deletions was found in both species, caused by unequal homologous recombination and illegitimate recombination. The relative contribution of unequal homologous recombination compared to illegitimate recombination was higher in rice than in Arabidopsis. However, retrotransposons are rapidly removed in both species, as evidenced by the similar apparent ages of intact elements (most less than 3 million years old) in these two plants and all other investigated plant species. CONCLUSIONS: Differences in the activity of mechanisms for retrotransposon regulation or deletion generation between species could explain current genome size variation without any requirement for natural selection to act on this trait, although the results do not preclude selection as a contributing factor. The simplest model suggests that significant genome size variation is generated by lineage-specific differences in the molecular mechanisms of DNA amplification and removal, creating major variation in nuclear DNA content that can then serve as the substrate for fitness-based selection.     

Efficient algorithms and software for detection of full-length LTR retrotransposons

LTR retrotransposons constitute one of the most abundant classes of repetitive elements in eukaryotic genomes. In this paper, we present a new algorithm for detection of full-length LTR retrotransposons in genomic sequences. The algorithm identifies regions in a genomic sequence that show structural characteristics of LTR retrotransposons. Three key components distinguish our algorithm from that of current software--(i) a novel method that preprocesses the entire genomic sequence in linear time and produces high quality pairs of LTR candidates in run-time that is constant per pair, (ii) a thorough alignment-based evaluation of candidate pairs to ensure high quality prediction, and (iii) a robust parameter set encompassing both structural constraints and quality controls providing users with a high degree of flexibility. We implemented our algorithm into a software program called LTR_par, which can be run on both serial and parallel computers. Validation of our software against the yeast genome indicates superior results in both quality and performance when compared to existing software. Additional validations are presented on rice BACs and chimpanzee genome.     

Fine-grained annotation and classification of de novo predicted LTR retrotransposons

Long terminal repeat (LTR) retrotransposons and endogenous retroviruses (ERVs) are transposable elements in eukaryotic genomes well suited for computational identification. De novo identification tools determine the position of potential LTR retrotransposon or ERV insertions in genomic sequences. For further analysis, it is desirable to obtain an annotation of the internal structure of such candidates. This article presents LTRdigest, a novel software tool for automated annotation of internal features of putative LTR retrotransposons. It uses local alignment and hidden Markov model-based algorithms to detect retrotransposon-associated protein domains as well as primer binding sites and polypurine tracts. As an example, we used LTRdigest results to identify 88 (near) full-length ERVs in the chromosome 4 sequence of Mus musculus, separating them from truncated insertions and other repeats. Furthermore, we propose a work flow for the use of LTRdigest in de novo LTR retrotransposon classification and perform an exemplary de novo analysis on the Drosophila melanogaster genome as a proof of concept. Using a new method solely based on the annotations generated by LTRdigest, 518 potential LTR retrotransposons were automatically assigned to 62 candidate groups. Representative sequences from 41 of these 62 groups were matched to reference sequences with >80% global sequence similarity.