LTR Retroelements Research Articles

GyDB Mobilomics

The Gypsy Database concerning Mobile Genetic Elements (release 2.0) is a wiki-style project devoted to the phylogenetic classification of LTR retroelements and their viral and host gene relatives characterized from distinct organisms. Furthermore, GyDB 2.0 is concerned with studying mobile elements within genomes. Therefore, an in-progress repository was created for databases with annotations of mobile genetic elements from particular genomes. This repository is called Mobilomics and the first uploaded database contains 549 LTR retroelements and related transposases which have been annotated from the genome of the Pea aphid Acyrthosiphon pisum. Mobilomics is accessible from the GyDB 2.0 project using the URL: http://gydb.org/index.php/Mobilomics.

LTR retroelements in the genome of Daphnia pulex

BackgroundLong terminal repeat (LTR) retroelements represent a successful group of transposable elements (TEs) that have played an important role in shaping the structure of many eukaryotic genomes. Here, we present a genome-wide analysis of LTR retroelements in Daphnia pulex, a cyclical parthenogen and the first crustacean for which the whole genomic sequence is available. In addition, we analyze transcriptional data and perform transposon display assays of lab-reared lineages and natural isolates to identify potential influences on TE mobility and differences in LTR retroelements loads among individuals reproducing with and without sex.ResultsWe conducted a comprehensive de novo search for LTR retroelements and identified 333 intact LTR retroelements representing 142 families in the D. pulex genome. While nearly half of the identified LTR retroelements belong to the gypsy group, we also found copia (95), BEL/Pao (66) and DIRS (19) retroelements. Phylogenetic analysis of reverse transcriptase sequences showed that LTR retroelements in the D. pulex genome form many lineages distinct from known families, suggesting that the majority are novel. Our investigation of transcriptional activity of LTR retroelements using tiling array data obtained from three different experimental conditions found that 71 LTR retroelements are actively transcribed. Transposon display assays of mutation-accumulation lines showed evidence for putative somatic insertions for two DIRS retroelement families. Losses of presumably heterozygous insertions were observed in lineages in which selfing occurred, but never in asexuals, highlighting the potential impact of reproductive mode on TE abundance and distribution over time. The same two families were also assayed across natural isolates (both cyclical parthenogens and obligate asexuals) and there were more retroelements in populations capable of reproducing sexually for one of the two families assayed.ConclusionsGiven the importance of LTR retroelements activity in the evolution of other genomes, this comprehensive survey provides insight into the potential impact of LTR retroelements on the genome of D. pulex, a cyclically parthenogenetic microcrustacean that has served as an ecological model for over a century.

The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences.

BackgroundIn today's age of genomic discovery, no attempt has been made to comprehensively sequence a gymnosperm genome. The largest genus in the coniferous family Pinaceae is Pinus, whose 110-120 species have extremely large genomes (c. 20-40 Gb, 2N = 24). The size and complexity of these genomes have prompted much speculation as to the feasibility of completing a conifer genome sequence. Conifer genomes are reputed to be highly repetitive, but there is little information available on the nature and identity of repetitive units in gymnosperms. The pines have extensive genetic resources, with approximately 329000 ESTs from eleven species and genetic maps in eight species, including a dense genetic map of the twelve linkage groups in Pinus taeda.ResultsWe present here the Sanger sequence and annotation of ten P. taeda BAC clones and Genome Analyzer II whole genome shotgun (WGS) sequences representing 7.5% of the genome. Computational annotation of ten BACs predicts three putative protein-coding genes and at least fifteen likely pseudogenes in nearly one megabase of sequence. We found three conifer-specific LTR retroelements in the BACs, and tentatively identified at least 15 others based on evidence from the distantly related angiosperms. Alignment of WGS sequences to the BACs indicates that 80% of BAC sequences have similar copies (≥ 75% nucleotide identity) elsewhere in the genome, but only 23% have identical copies (99% identity). The three most common repetitive elements in the genome were identified and, when combined, represent less than 5% of the genome.ConclusionsThis study indicates that the majority of repeats in the P. taeda genome are 'novel' and will therefore require additional BAC or genomic sequencing for accurate characterization. The pine genome contains a very large number of diverged and probably defunct repetitive elements. This study also provides new evidence that sequencing a pine genome using a WGS approach is a feasible goal.

Structural analysis of 83-kb genomic DNA from Thellungiella halophila: Sequence features and microcolinearity between salt cress and Arabidopsis thaliana

Salt cress ( Thellungiella halophila) has become a desirable plant model for molecular analysis of the mechanisms of salt tolerance. Analysis of its physiological action and expressed EST has resulted in better understanding. However, less is known about its genomic features. Here we determined a continuous sequence approximately 83 kb from a salt cress BAC clone, providing the first insight into the genomic feature for this species. The gene density is approximately one gene per 3.6 kb in this sequence. Many types of repetitive sequences are present in this salt cress sequence, including LTR retroelements, DNA transposons and a number of simple sequence repeats. Comparison of sequence similarity indicated that salt cress shares a close relationship with Arabidopsis. Extensive conservation and high-level microcolinearity were uncovered for both genomes. Our study also indicated that genomic DNA alternations (involving chromosome inversion, sequence loss and gene translocation) contributed to the genomic discrepancies between salt cress and Arabidopsis.

Network dynamics of eukaryotic LTR retroelements beyond phylogenetic trees

BackgroundSequencing projects have allowed diverse retroviruses and LTR retrotransposons from different eukaryotic organisms to be characterized. It is known that retroviruses and other retro-transcribing viruses evolve from LTR retrotransposons and that this whole system clusters into five families: Ty3/Gypsy, Retroviridae, Ty1/Copia, Bel/Pao and Caulimoviridae. Phylogenetic analyses usually show that these split into multiple distinct lineages but what is yet to be understood is how deep evolution occurred in this system.ResultsWe combined phylogenetic and graph analyses to investigate the history of LTR retroelements both as a tree and as a network. We used 268 non-redundant LTR retroelements, many of them introduced for the first time in this work, to elucidate all possible LTR retroelement phylogenetic patterns. These were superimposed over the tree of eukaryotes to investigate the dynamics of the system, at distinct evolutionary times. Next, we investigated phenotypic features such as duplication and variability of amino acid motifs, and several differences in genomic ORF organization. Using this information we characterized eight reticulate evolution markers to construct phenotypic network models.ConclusionThe evolutionary history of LTR retroelements can be traced as a time-evolving network that depends on phylogenetic patterns, epigenetic host-factors and phenotypic plasticity. The Ty1/Copia and the Ty3/Gypsy families represent the oldest patterns in this network that we found mimics eukaryotic macroevolution. The emergence of the Bel/Pao, Retroviridae and Caulimoviridae families in this network can be related with distinct inflations of the Ty3/Gypsy family, at distinct evolutionary times. This suggests that Ty3/Gypsy ancestors diversified much more than their Ty1/Copia counterparts, at distinct geological eras. Consistent with the principle of preferential attachment, the connectivities among phenotypic markers, taken as network-represented combinations, are power-law distributed. This evidences an inflationary mode of evolution where the system diversity; 1) expands continuously alternating vertical and gradual processes of phylogenetic divergence with episodes of modular, saltatory and reticulate evolution; 2) is governed by the intrinsic capability of distinct LTR retroelement host-communities to self-organize their phenotypes according to emergent laws characteristic of complex systems.ReviewersThis article was reviewed by Eugene V. Koonin, Eric Bapteste, and Enmanuelle Lerat (nominated by King Jordan)

Bioinformatic flowchart and database to investigate the origins and diversity of Clan AA peptidases

BackgroundClan AA of aspartic peptidases relates the family of pepsin monomers evolutionarily with all dimeric peptidases encoded by eukaryotic LTR retroelements. Recent findings describing various pools of single-domain nonviral host peptidases, in prokaryotes and eukaryotes, indicate that the diversity of clan AA is larger than previously thought. The ensuing approach to investigate this enzyme group is by studying its phylogeny. However, clan AA is a difficult case to study due to the low similarity and different rates of evolution. This work is an ongoing attempt to investigate the different clan AA families to understand the cause of their diversity.ResultsIn this paper, we describe in-progress database and bioinformatic flowchart designed to characterize the clan AA protein domain based on all possible protein families through ancestral reconstructions, sequence logos, and hidden markov models (HMMs). The flowchart includes the characterization of a major consensus sequence based on 6 amino acid patterns with correspondence with Andreeva's model, the structural template describing the clan AA peptidase fold. The set of tools is work in progress we have organized in a database within the GyDB project, referred to as Clan AA Reference Database .ConclusionThe pre-existing classification combined with the evolutionary history of LTR retroelements permits a consistent taxonomical collection of sequence logos and HMMs. This set is useful for gene annotation but also a reference to evaluate the diversity of, and the relationships among, the different families. Comparisons among HMMs suggest a common ancestor for all dimeric clan AA peptidases that is halfway between single-domain nonviral peptidases and those coded by Ty3/Gypsy LTR retroelements. Sequence logos reveal how all clan AA families follow similar protein domain architecture related to the peptidase fold. In particular, each family nucleates a particular consensus motif in the sequence position related to the flap. The different motifs constitute a network where an alanine-asparagine-like variable motif predominates, instead of the canonical flap of the HIV-1 peptidase and closer relatives.ReviewersThis article was reviewed by Daniel H. Haft, Vladimir Kapitonov (nominated by Jerry Jurka), and Ben M. Dunn (nominated by Claus Wilke).

The LTR retrotransposon micropia in the cardini group of Drosophila (Diptera: Drosophilidae): a possible case of horizontal transfer

The presence of the micropia retroelement from the Ty1-copia family of LTR retroelements was investigated in three species of the Drosophila cardini group. Southern blot analysis suggested the existence of at least four micropia copies in the genomes of D. cardinoides, D. neocardini and D. polymorpha populations. The high sequence similarity between dhMiF2 and Dm11 clones (micropia retroelements isolated from D. hydei and D. melanogaster, respectively) with micropia sequences amplified from D. cardini group genome supports the hypothesis that this retroelement plays an active role in horizontal transfer events between D. hydei and the D. cardini group.

Relationships of Gag-pol Diversity between Ty3/Gypsy and Retroviridae LTR retroelements and the three kings hypothesis

BackgroundThe origin of vertebrate retroviruses (Retroviridae) is yet to be thoroughly investigated, but due to their similarity and identical gag-pol (and env) genome structure, it is accepted that they evolve from Ty3/Gypsy LTR retroelements the retrotransposons and retroviruses of plants, fungi and animals. These 2 groups of LTR retroelements code for 3 proteins rarely studied due to the high variability – gag polyprotein, protease and GPY/F module. In relation to 3 previously proposed Retroviridae classes I, II and II, investigation of the above proteins conclusively uncovers important insights regarding the ancient history of Ty3/Gypsy and Retroviridae LTR retroelements.ResultsWe performed a comprehensive study of 120 non-redundant Ty3/Gypsy and Retroviridae LTR retroelements. Phylogenetic reconstruction inferred based on the concatenated analysis of the gag and pol polyproteins shows a robust phylogenetic signal regarding the clustering of OTUs. Evaluation of gag and pol polyproteins separately yields discordant information. While pol signal supports the traditional perspective (2 monophyletic groups), gag polyprotein describes an alternative scenario where each Retroviridae class can be distantly related with one or more Ty3/Gypsy lineages. We investigated more in depth this evidence through comparative analyses performed based on the gag polyprotein, the protease and the GPY/F module. Our results indicate that contrary to the traditional monophyletic view of the origin of vertebrate retroviruses, the Retroviridae class I is a molecular fossil, preserving features that were probably predominant among Ty3/Gypsy ancestors predating the split of plants, fungi and animals. In contrast, classes II and III maintain other phenotypes that emerged more recently during Ty3/Gypsy evolution.ConclusionThe 3 Retroviridae classes I, II and III exhibit phenotypic differences that delineate a network never before reported between Ty3/Gypsy and Retroviridae LTR retroelements. This new scenario reveals how the diversity of vertebrate retroviruses is polyphyletically recurrent into the Ty3/Gypsy evolution, i.e. older than previously thought. The simplest hypothesis to explain this finding is that classes I, II and III trace back to at least 3 Ty3/Gypsy ancestors that emerged at different evolutionary times prior to protostomes-deuterostomes divergence. We have called this "the three kings hypothesis" concerning the origin of vertebrate retroviruses.

Origin and spread of the SRY gene on the X and Y chromosomes of the rodent Microtus cabrerae: Role of L1 elements

In the rodent species Microtus cabrerae, males as well as females present several copies of the SRY gene, a single-copy gene located on the Y chromosome in most mammals. Using different PCR approaches, we have characterized the sequence, structure, and organization of the SRY copies and their flanking regions distributed on the X and Y chromosomes of this species. All copies of SRY analyzed, including those from the Y chromosome, proved to be nonfunctional pseudogenes, as they have internal stop codons. In addition, we demonstrated the association of SRY pseudogenes with different fragments of L1 and LTR retroelements in both sex chromosomes of M. cabrerae. Examining the possible origin of SRY pseudogene and retroposons association, we propose that retroposons could have been involved in the mechanism of SRY gene amplification on the Y chromosome and in the transference of the Y-linked SRY copies to the X-chromosome heterochromatin.

Distinct chromosomal distributions of highly repetitive sequences in maize

The majority of genomic DNA in most plant species is made up of repetitive elements including satellites and retrotransposons. The maize genome is intermediate in size and abundance of repetitive elements between small genomes such as Arabidopsis and rice and larger genomes such as wheat. Although repetitive elements are present throughout the maize genome, individual families are non-randomly distributed along chromosomes. In this work we use fluorescence in-situ hybridization (FISH) to examine the distribution of abundant LTR retroelement families and satellites contained in heterochromatic blocks called knobs. Different retroelement families have distinct patterns of hybridization. Prem1 and Tekay, two very closely related elements, both hybridize along the length of all chromosomes but do so with greater intensity near the centromeres, although subtle differences are detectable between the hybridization patterns. Opie, Prem2/Ji, and Huck are enriched away from the centromeres and Grande is distributed uniformly along the chromosomes. Double labeling with proximally and distally enriched elements on pachytene chromosomes produces alternating blocks of element enrichment. The maize elements hybridized in the same general patterns to chromosomes of maize relatives including Zea diploperennis and Tripsacum dactyloides. Additionally, abundant Tripsacum LTR retroelements are enriched in similar chromosomal regions among the different species. The 180 bp knob satellite is present in large blocks at interstitial locations on chromosome arms. With long exposures, smaller sites of hybridization are detected at the ends of chromosomes, adjacent to the telomere tract. This distal position for knob satellites is conserved among Zea and Tripsacum species.

The 5'-end transitional CpGs between the CpG islands and retroelements are hypomethylated in association with loss of heterozygosity in gastric cancers

BackgroundA loss of heterozygosity (LOH) represents a unilateral chromosomal loss that reduces the dose of highly repetitive Alu, L1, and LTR retroelements. The aim of this study was to determine if the LOH events can affect the spread of retroelement methylation in the 5'-end transitional area between the CpG islands and their nearest retroelements.MethodsThe 5'-transitional area of all human genes (22,297) was measured according to the nearest retroelements to the transcription start sites. For 50 gastric cancer specimens, the level of LOH events on eight cancer-associated chromosomes was estimated using the microsatellite markers, and the 5'-transitional CpGs of 20 selected genes were examined by methylation analysis using the bisulfite-modified DNA.ResultsThe extent of the transitional area was significantly shorter with the nearest Alu elements than with the nearest L1 and LTR elements, as well as in the extragenic regions containing a higher density of retroelements than in the intragenic regions. The CpG islands neighbouring a high density of Alu elements were consistently hypomethylated in both normal and tumor tissues. The 5'-transitional methylated CpG sites bordered by a low density of Alu elements or the L1 and LTR elements were hypomethylated more frequently in the high-level LOH cases than in the low-level LOH cases.ConclusionThe 5'-transitional methylated CpG sites not completely protected by the Alu elements were hypomethylated in association with LOH events in gastric cancers. This suggests that an irreversible unbalanced decrease in the genomic dose reduces the spread of L1 methylation in the 5'-end regions of genes.

Identification of transposons, retroelements, and a gene family predominantly expressed in floral tissues in chromosome 3DS of the hexaploid wheat progenitor Aegilops tauschii

A multigene family expressed during early floral development was identified on the short arm of wheat chromosome 3D in the region of the Ph2 locus, a locus controlling homoeologous chromosome pairing in allohexaploid wheat. Physical, genetic and molecular characterisation of the Wheat Meiosis 1 (WM1) gene family identified seven members that localised within a region of 173-kb. WM1 gene family members were sequenced and they encode mainly type Ia plasma membrane-anchored leucine rich repeat-like receptor proteins. In situ expression profiling suggests the gene family is predominantly expressed in floral tissue. In addition to the WM1 gene family, a number of other genes, gene fragments and pseudogenes were identified. It has been predicted that there is approximately one gene every 19-kb and that this region of the wheat genome contains 23 repetitive elements including BARE-1 and Wis2-1 like sequences. Nearly 50% of the repetitive elements identified were similar to known transposons from the CACTA superfamily. Ty1-copia, Ty3-gypsy and Athila LTR retroelements were also prevalent within the region. The WM1 gene cluster is present on 3DS and on barley 3HS but missing from the A and B genomes of hexaploid wheat. This suggests either recent generation of the cluster or specific deletion of the cluster during wheat polyploidisation. The evolutionary significance of the cluster, its possible roles in disease response or floral and early meiotic development and its location at or near the Ph2 locus are discussed.

APOBEC3A Is a Potent Inhibitor of Adeno-Associated Virus and Retrotransposons

APOBEC3 proteins constitute a family of cytidine deaminases that provide intracellular resistance to retrovirus replication and transposition of endogenous retroelements. One family member, APOBEC3A (hA3A), is an orphan, without any known antiviral activity. We show that hA3A is catalytically active and that it, but none of the other family members, potently inhibits replication of the parvovirus adeno-associated virus (AAV). hA3A was also a potent inhibitor of the endogenous LTR retroelements, MusD, IAP, and the non-LTR retroelement, LINE-1. Its function was dependent on the conserved amino acids of the hA3A active site, consistent with a role for cytidine deamination, although mutations in retroelement sequences were not found. These findings demonstrate the potent activity of hA3A, an APOBEC3 family member with no previously identified function. They also highlight the functional differences between APOBEC3 proteins. The APOBEC3 family members have distinct functions and may have evolved to resist various classes of genetic elements.

The VIPER elements of trypanosomes constitute a novel group of tyrosine recombinase-enconding retrotransposons

VIPER was initially characterized as a 2326 bp LTR-like retroelement associated to SIRE, a short interspersed repetitive element specific of Trypanosoma cruzi. It carried a single ORF that coded for a putative reverse transcriptase-RNAse H protein, suggesting that it could be a truncated copy of a longer retroelement. Herein we report the identification and characterization of a complete 4480 bp long VIPER in the T. cruzi genome. The complete VIPER harbored three non-overlapped domains encoding for a GAG-like, a tyrosine recombinase and a reverse transcriptase-RNAse H proteins. VIPER elements were also found in the genomes of Trypanosoma brucei and Trypanosoma vivax, but not in Leishmania sp. On the basis of its reverse transcriptase phylogeny, VIPER was classified as an LTR retroelement. However, VIPER was structurally related to the tyrosine recombinase encoding retroelements, DIRS and Ngaro. Phylogenetic analysis showed that VIPER's tyrosine recombinase grouped with the transposases RCI1 of Escherichia coli and Ye24 and Ye72 of Haemophilus influenzae within a major branch of prokaryotic recombinases. Taken together, VIPER's structure, the nature of its tyrosine recombinase, the unique features of its reverse transcriptase catalytic consensus motif and the fact that it was found in Trypanosomes, an early branching eukaryote, suggest that VIPER may be the closest relative of the founder element of the tyrosine recombinase encoding retrotransposons known up to date. Our analysis revealed that tyrosine recombinase-encoding retroelements were originated as early in evolution as non-LTR retroelements and suggests that VIPER, Ngaro and DIRS elements may constitute a third group of retrotransposons, distinct from both LTR and non-LTR retroelements.

Expression of the Idefix retrotransposon in early follicle cells in the germarium of Drosophila melanogaster is determined by its LTR sequences and a specific genomic context.

Retrotransposons are transcriptionally activated in different tissues and cell types by a variety of genomic and environmental factors. Transcription of LTR retrotransposons is controlled by cis-acting regulatory sequences in the 5' LTR. Mobilization of two LTR retroelements, Idefix and ZAM, occurs in the unstable RevI line of Drosophila melanogaster, in which their copy numbers are high, while they are low in all other stocks tested. Here we show that both a full-length and a subgenomic Idefix transcript that are necessary for its mobilization are present in the Rev1 line, but not in the other lines. Studies on transgenic strains demonstrate that the 5' LTR of Idefix contains sequences that direct the tissue-specific expression of the retroelement in testes and ovaries of adult flies. In ovaries, expression occurs in the early follicle and in other somatic cells of the germarium, and is strictly associated with the unstable genetic context conferred by the RevI line. Control of tissue-specific Idefix expression by interactions between cis-acting sequences of its LTR and trans-acting genomic factors provides an opportunity to use this retroelement as a tool for the study of the early follicle cell lineage in the germarium.

Moose, a new family of LTR-retrotransposons in the mosquito Anopheles gambiae.

A novel LTR retrotransposable element called moose has been cloned and characterized from the malaria transmitting mosquito, Anopheles gambiae. This element has all the characteristic features of LTR retroelements and is related to retroelements from other insects and nematodes, belonging to a subgroup of retroelements distinct from the copia/Ty1 and gypsy/Ty3 groups. The moose element appears to be active in A. gambiae, and strong RNA expression is detected in the male and female gonads. The use of this retroelement as a potential vector for germ line transformation is discussed.

A hotspot for the Drosophila gypsy retroelement in the ovo locus.

The Drosophila retroelement gypsy has a number of unusual features including an unusual LTR terminal sequence and an apparent target sequence preference. The ovo locus is a known hotspot for gypsy insertion. We examined the target sequence preference of gypsy within ovo by isolating 26 new insertions and sequencing the gypsy/ovo junctions. Insertions were found at multiple sites within the ovo locus. The insertions clustered within an approximately 150 bp region in the non-translated region of the ovo beta transcript, with most insertions falling within the first intron. There were seven sites of insertion within this region and these mostly conform to the consensus sequence YRYRYR (where Y = pyrimidine and R = purine). However, this target sequence is at best necessary but not sufficient to specify a hotspot, as there were several other sequences conforming to this consensus in the ovo locus that were not hit. The results indicate that gypsy may have a higher degree of target specificity than most infectious LTR retroelements.

LTR retrotransposons and the evolution of eukaryotic enhancers.

Since LTR retrotransposons and retroviruses are especially prone to regional duplications and recombination events, these viral-like systems may be especially conducive to the evolution of closely spaced combinatorial regulatory motifs. Using the Drosophila copia LTR retrotransposon as a model, we show that a regulatory region contained within the element's untranslated leader region (ULR) consists of multiple copies of an 8 bp motif (TTGTGAAA) with similarity to the core sequence of the SV40 enhancer. Naturally occurring variation in the number of these motifs is correlated with the enhancer strength of the ULR. Our results indicate that inter-element selection may favor the evolution of more active enhancers within permissive genetic backgrounds. We propose that LTR retroelements and perhaps other retrotransposons constitute drive mechanisms for the evolution of eukaryotic enhancers which can be subsequently distributed throughout host genomes to play a role in regulatory evolution.