A genome-wide screening of BEL-Pao like retrotransposons in Anopheles gambiae by the LTR_STRUC program

The advanced status of assembly of the nematoceran Anopheles gambiae genomic sequence allowed us to perform a wide genome analysis to looking at the presence of Long Terminal Repeats (LTRs) in the range of 10 kb by means of the LTR_STRUC tool. More than three hundred sequences were retrieved and 210 were treated as putative complete retrotransposons that were individually analysed with respect to known retrotransposons of A. gambiae and D. melanogaster. The results show that the vast majority of the retrotransposons analysed belong to the Ty3/gypsy class and only 8% to the Ty1/copia class. In addition, phylogenetic analysis allowed us to characterize in more detail the relationship of a large BEL-Pao lineage in which a single family was shown to harbour an additional env gene.     

家蚕LTR逆转录转座子的鉴定、分类及系统发育分析

转座子是真核生物基因组的重要组成成分。为了研究家蚕Bombyx mori长末端重复序列 (long terminal repeat, LTR)逆转录转座子的分类及进化, 本研究采用de novo预测和同源性搜索相结合的方法, 在家蚕基因组中共鉴定出了38个LTR逆转录转座子家族, 序列长度占整个基因组的0.64%, 远小于先前预测的11.8%, 其中有6个家族为本研究的新发现。38个家族中, 26个家族有表达序列标签 (expression sequence tag, EST)证据, 表明这些家族具有潜在的活性。对有EST证据的6个家族和没有EST证据的5个家族用RT-PCR进行了组织表达谱实验, 结果表明这11个家族在一些组织中有表达, 这进一步证实了这些家族具有转录活性, 基于此我们推测家蚕中大部分的LTR逆转录转座子家族很可能具有潜在活性。对转座子的插入时间进行估计, 结果表明绝大部分元件都是最近1百万年内插入到家蚕基因组中的。我们还比较了黑腹果蝇Drosophila melanogaster、 冈比亚按蚊Anopheles gambiae和家蚕B. mori中Ty3/Gypsy超家族分支的差异, 结果表明不同枝在不同昆虫中有着不同的扩张。家蚕中LTR逆转录转座子的鉴定和系统分析有助于我们理解逆转录转座子在昆虫进化中的作用。     

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.     

Diversity of LTR-retrotransposons and <Emphasis Type="Italic">Enhancer/Suppressor</Emphasis><Emphasis Type="Italic">Mutator</Emphasis>-like transposons in cassava (<Emphasis Type="Italic">Manihot esculenta</Emphasis> Crantz)

Cassava (Manihot esculenta Crantz), though a major world crop with enormous potential, is very under studied. Little is known about its genome structure and organisation. Transposable elements have a key role in the evolution of genome structure, and can be used as important tools in applied genetics. This paper sets out to survey the diversity of members of three major classes of transposable element within the cassava genome and in relation to similar elements in other plants. Members of two classes of LTR-retrotransposons, Ty1/copia-like and Ty3/gypsy-like, and of Enhancer/Suppressor Mutator (En/Spm)-like transposons were isolated and characterised. Analyses revealed 59 families of Ty1/copia, 26 families of Ty3/gypsy retrotransposons, and 40 families of En/Spm in the cassava genome. In the comparative analyses, the predicted amino acid sequences for these transposon classes were compared with those of related elements from other plant species. These revealed that there were multiple lineages of Ty1/copia-like retrotransposons in the genome of cassava and suggested that vertical and horizontal transmission as the source of cassava Mecops may not be mutually exclusive. For the Ty3/gypsy elements network, two groups of cassava Megyps were evident including the Arabidopsis Athila lineage. However, cassava En/Spm-like elements (Meens) constituted a single group within a network of plant En/Spm-like elements. Hybridisation analysis supported the presence of transposons in the genome of cassava in medium (Ty3/gypsy and En/Spm) to high (Ty1/copia) copy numbers. Thus the cassava genome was shown to contain diverse members of three major classes of transposable element; however, the different classes exhibited contrasting evolutionary histories.     

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.     

Ty3/Gypsy retrotransposons: description of new Arabidopsis thaliana elements and evolutionary perspectives derived from comparative genomic data

We performed a comprehensive analysis of the evolution of the Ty3/GYPSY: group of long-terminal-repeat retrotransposons (also known as METAVIRIDAE:). Exhaustive database searches allowed us to detect novel elements of this group. In particular, the Arabidopsis thaliana and Drosophila melanogaster genome sequencing projects have recently disclosed a large number of new Ty3/GYPSY: sequences. So far, elements of three different Ty3/GYPSY: lineages had been described for A. thaliana. Here, we describe six new lineages, which we have called Tit-for-tat1, Tit-for-tat2, Gimli, Gloin, Legolas, and Little Athila. We confirm that plant Ty3/GYPSY: elements form two main monophyletic groups. Moreover, our results suggest that at least four independent ancestral lineages existed before the monocot-dicot split, about 200 MYA. Twelve sequences from D. melanogaster that may correspond to new elements are also described. Some of these sequences are similar to those of OSVALDO: and Ulysses, two elements of the OSVALDO: clade that had never before been described for D. melanogaster. Comparative analyses of multiple organisms, some of them with completely sequenced genomes, show that the number of lineages of Ty3/GYPSY: elements is very variable. Thus, while only 1 lineage is present in Saccharomyces cerevisiae, at least 6 exist in Caenorhabditis elegans, at least 9 are present in the A. thaliana, and perhaps 20 are present in D. melanogaster. Finally, we suggest that the presence of a chromodomain-containing integrase, a feature of some closely related Ty3/GYPSY: elements of fungi, plants, and animals, may be used to define a new METAVIRIDAE: genus.     

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.     

Evidence for horizontal transfer of the LTR retrotransposon mdg3, which lacks an env gene

Horizontal (interspecific) transfer is regarded as a possible strategy for the propagation of transposable elements through evolutionary time. To date, however, conclusive evidence that transposable elements are capable of horizontal transfer from one species to another has been limited to class II or DNA-type elements. We tested the possibility of such transfer for several Drosophila melanogaster LTR retrotransposons of the gypsy group in an experiment in which D. melanogaster and D. virilis somatic cell lines were used as donor and recipient cells, respectively. This approach was chosen in light of the high levels of LTR retrotransposon amplification and expression observed in cultured D. melanogaster cells. In the course of the experiment, parallel analysis for mdg1, mdg3, 17.6, 297, 412 and B104/roo retrotransposons was performed to detect their presence in the genome of recipient cells. Only the mdg3 retrotransposon, which lacks an env gene, was found to be transmitted into recipient cells. This model, based on the use of cultured cells, is a promising system for further investigating the mechanisms of LTR retrotransposon transfer.     

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.     

Increased Length of Long Terminal Repeats Inhibits Ty1 Transposition and Leads to the Formation of Tandem Multimers

The Ty1 retrotransposon of Saccharomyces cerevisiae is bounded by long-terminal repeats (LTRs). We have constructed a variety of Ty1 elements in which the LTR length has been increased from the normal length of 334 bp to >2 kb. Although small insertions in the LTR have minimal effects on transposition frequency, larger insertions dramatically reduce it. Nevertheless, elements with long LTRs are incorporated into the genome at a low frequency. Most of these rare insertion events represent Ty1 tandem (head to tail) multimers.     

Mosquitoes LTR retrotransposons: a deeper view into the genomic sequence of Culex quinquefasciatus

A set of 67 novel LTR-retrotransposon has been identified by in silico analyses of the Culex quinquefasciatus genome using the LTR_STRUC program. The phylogenetic analysis shows that 29 novel and putatively functional LTR-retrotransposons detected belong to the Ty3/gypsy group. Our results demonstrate that, by considering only families containing potentially autonomous LTR-retrotransposons, they account for about 1% of the genome of C. quinquefasciatus. In previous studies it has been estimated that 29% of the genome of C. quinquefasciatus is occupied by mobile genetic elements.The potential role of retrotransposon insertions strictly associated with host genes is described and discussed along with the possible origin of a retrotransposon with peculiar Primer Binding Site region. Finally, we report the presence of a group of 38 retrotransposons, carrying tandem repeated sequences but lacking coding potential, and apparently lacking "master copy" elements from which they could have originated. The features of the repetitive sequences found in these non-autonomous LTR retrotransposons are described, and their possible role discussed.These results integrate the existing data on the genomics of an important virus-borne disease vector.     

A widespread occurrence of extra open reading frames in plant Ty3/gypsy retrotransposons

Long terminal repeat (LTR) retrotransposons make up substantial parts of most higher plant genomes where they accumulate due to their replicative mode of transposition. Although the transposition is facilitated by proteins encoded within the gag-pol region which is common to all autonomous elements, some LTR retrotransposons were found to potentially carry an additional protein coding capacity represented by extra open reading frames located upstream or downstream of gag-pol. In this study, we performed a comprehensive in silico survey and comparative analysis of these extra open reading frames (ORFs) in the group of Ty3/gypsy LTR retrotransposons as the first step towards our understanding of their origin and function. We found that extra ORFs occur in all three major lineages of plant Ty3/gypsy elements, being the most frequent in the Tat lineage where most (77?%) of identified elements contained extra ORFs. This lineage was also characterized by the highest diversity of extra ORF arrangement (position and orientation) within the elements. On the other hand, all of these ORFs could be classified into only two broad groups based on their mutual similarities or the presence of short conserved motifs in their inferred protein sequences. In the Athila lineage, the extra ORFs were confined to the element 3' regions but they displayed much higher sequence diversity compared to those found in Tat. In the lineage of Chromoviruses the extra ORFs were relatively rare, occurring only in 5' regions of a group of elements present in a single plant family (Poaceae). In all three lineages, most extra ORFs lacked sequence similarities to characterized gene sequences or functional protein domains, except for two Athila-like elements with similarities to LOGL4 gene and part of the Chromoviruses extra ORFs that displayed partial similarity to histone H3 gene. Thus, in these cases the extra ORFs most likely originated by transduction or recombination of cellular gene sequences. In addition, the protein domain which is otherwise associated with DNA transposons have been detected in part of the Tat-like extra ORFs, pointing to their origin from an insertion event of a mobile element.     

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.     

Evolutionary Genomics of Transposable Elements in Saccharomyces cerevisiae

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.