Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Losing identity: structural diversity of transposable elements belonging to different classes in the genome of Anopheles gambiae

Rita D Fernández-Medina1*, José M C Ribeiro2, Claudia M A Carareto3, Luciane Velasque4 and Cláudio J Struchiner15

Author Affiliations

1 Escola Nacional de Saúde Pública Sergio Arouca, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil

2 Section of Vector Biology, Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, MD, 20852, USA

3 Laboratório de Evolução Molecular, Departamento de Biologia, UNESP, Universidade Estadual Paulista, 15054-000, São José do Rio Preto, SP, Brazil

4 Departamento de Matemática e Estatística, Universidade Federal do Estado do Rio de Janeiro, Rio de Janeiro, Brazil

5 Instituto de Medicina Social/Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil

For all author emails, please log on.

BMC Genomics 2012, 13:272  doi:10.1186/1471-2164-13-272


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/272


Received:6 February 2012
Accepted:8 June 2012
Published:22 June 2012

© 2012 Fernandez-Medina et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Transposable elements (TEs), both DNA transposons and retrotransposons, are genetic elements with the main characteristic of being able to mobilize and amplify their own representation within genomes, utilizing different mechanisms of transposition. An almost universal feature of TEs in eukaryotic genomes is their inability to transpose by themselves, mainly as the result of sequence degeneration (by either mutations or deletions). Most of the elements are thus either inactive or non-autonomous. Considering that the bulk of some eukaryotic genomes derive from TEs, they have been conceived as “TE graveyards.” It has been shown that once an element has been inactivated, it progressively accumulates mutations and deletions at neutral rates until completely losing its identity or being lost from the host genome; however, it has also been shown that these “neutral sequences” might serve as raw material for domestication by host genomes.

Results

We have analyzed the sequence structural variations, nucleotide divergence, and pattern of insertions and deletions of several superfamilies of TEs belonging to both class I (long terminal repeats [LTRs] and non-LTRs [NLTRs]) and II in the genome of Anopheles gambiae, aiming at describing the landscape of deterioration of these elements in this particular genome. Our results describe a great diversity in patterns of deterioration, indicating lineage-specific differences including the presence of Solo-LTRs in the LTR lineage, 5′-deleted NLTRs, and several non-autonomous and MITEs in the class II families. Interestingly, we found fragments of NLTRs corresponding to the RT domain, which preserves high identity among them, suggesting a possible remaining genomic role for these domains.

Conclusions

We show here that the TEs in the An. gambiae genome deteriorate in different ways according to the class to which they belong. This diversity certainly has implications not only at the host genomic level but also at the amplification dynamic and evolution of the TE families themselves.

Keywords:
Transposable elements; LTR; Non-LTR; Class II; Deterioration; Anopheles gambiae

Background

Transposable elements (TEs) are genetic elements that share the main characteristic of amplifying their own representation within genomes. Due to their ability to spread in the absence of selection at the host level, they persist in genomes even at the expense of a net negative fitness to the hosts [1]. Under this view, these elements are conceived as genomic parasites. Today, the view of TEs within genomes has changed considerably, and these elements have been shown to be major contributors to both genome evolution and function [2-5]. The relationships between TEs and host genomes where they reside are now considered a continuum from parasitic to beneficial to genomes [6]. TEs show very different genetic structures and transposition strategies and accordingly have been classified into two classes (I and II) and subsequently into orders, superfamilies, and families [7]. Class I elements— also called retrotransposons, which depend on a RNA intermediary for their replication—are further classified into two orders: LTRs (which harbor flanking long terminal repeats [LTRs]) and non-LTRs (NLTRs), which share the transposition mechanism known as “copy and paste.” Both types of elements are first transcribed into mRNA, but LTRs are then retrotranscribed into a DNA copy that is later inserted into the genome, while NLTRs are retrotranscribed at the time they are inserted into the genome. Both processes lead to generation of novel copies. Class II—also called DNA transposons—present a simpler genetic structure that includes a gene encoding for the enzyme transposase (responsible for the element transposition) and flanking terminal inverted repeats (TIRs), which are recognition sites for transposase. Their transposition mechanism is known as “cut and paste” because they are usually excised from one strand at its original place and inserted into a new genomic position. Increasing copy number depends on strand repair mechanisms from homologous chromatids that are dependent on host enzymes. These intrinsic differences imply differences in the way elements belonging to the different classes and orders amplify and degenerate within genomes.

A TE family “life cycle” has been conceived as the successive steps from its “birth” in a given genome until its “death”, and includes the invasion, amplification, inactivation, and elimination from a genome that can occur in geological times. It has been shown that once an element has been inactivated, it progressively accumulates mutations and indels at neutral rates, until completely losing its identity or being lost from the host genome. It has also been shown that these deteriorated elements behave in the genome as “neutral sequences” and might serve as raw material for domestication by host genomes. The concept of “molecular domestication” was used to describe a process whereby a TE sequence is co-opted to perform a different role from the original function that it was selected for and that benefits the host genome [8]. Indeed, TE truncated copies can modulate host gene expression by providing new regulatory sequences, alternative splice sites, polyadenylation signals [9], and new transcription factor binding sites [10], as well as in post transcriptional regulation through RNA editing and translation regulation (reviewed in [11]). In addition, it has been shown that several microRNA genes derive from TEs [12].

An almost universal feature of TEs now in eukaryotic genomes is their inability to transpose by themselves, mainly as the result of sequence degradation although even TEs that have lost their functional transposition machinery can continue to be mobilized by other intact element products. For example, a class II element that has mutated its transposase but not its inverted repeats can be mobilized by the intact transposase of another element of the same family. Accordingly, elements can also be classified as autonomous and non-autonomous according to the nature of their mobilization. TEs thus are also genomically interesting because they can persist and live in the genome as dead elements [13] (after inactivation by mutations, indels, or recombination), and these “fossil sequences” will continue to evolve in the genome, leaving traces of their past history. The relationship between TEs belonging to different families and superfamilies has been recently analyzed as the relationship between different populations and species in community ecology [14,15]. Within this scope, autonomous and non-autonomous elements of the same family are considered as competitors. TEs are simultaneously part of the genome and independent entities living their own life within the genome [16]. TE dynamics can be analyzed both at the intra-genomic level (the different TE families considered as different populations in ecology) and at the intra-population level (considering the complex networks of interactions between autonomous and non-autonomous elements, besides the relationships between different families of elements residing in the same genome). This view is of particular interest when asking about the relative abundances and degrees of activity or level of degradation of TEs in a given genome. Then, in genomes where several TE families coexist, the amplification of one family could have an impact on the dynamics of other families [16]. Frequencies, intragenomic dynamics, distributions, and abundances of retrotransposons and transposons differ considerably between different species.

The description of the “deterioration landscape” of the TE families in a given genome might shed light on the distribution and abundance of the TE families coexisting in a genome. Although the various sequencing projects have allowed the study of the TE content in different genomes, the aspects related to deteriorated elements have been barely studied [17]. Polymorphisms seen at the nucleotide level of sequences belonging to the same family can give information regarding the way that TEs are transposed and regulated as well as how they degenerate (by insertions, deletions, substitutions, or rearrangements).

Differences in these aspects could indicate the existence of regulatory mechanisms acting in a class-specific way. In fact, self-regulation [18], which has been shown to depend on the TE family, regulation by mutant copies, and genomic regulation [19,20] are mechanisms that have all been described in several TE families. On the other hand, host genomes have evolved mechanisms to control or restrict TE replication. Chromatin remodeling [21], methylation of TE DNA sequences [22], and non-coding small RNAs [23,24] or cytidine deaminases [25] are examples of these mechanisms (wee [26] for a review). We have previously analyzed the TE families present in AnoTExcel, an online Anopheles gambiae TE-specific database [27]. In that database, we presented the general features of the TE landscape in the malaria mosquito genome. AnoTExcel presents all the individual sequences belonging to different families in the mosquito genome, which allows analysis of the dynamics and demographic trajectory of certain TE families [28]. Now, to better characterize the An. gambiae mobilome, we have analyzed the sequence structural variations, focusing specifically on the TE deteriorated landscape of elements belonging to classes I and II. We describe a diversity of patterns of deterioration, indicating lineage-specific differences including the presence of Solo-LTRs, 5′-deleted NLTRs, and several non-autonomous and MITEs belonging to class II families.

Results and discussion

Although different software tools have been used [29-36] for successful identification and/or characterization of repeats in eukaryotic genomes, so far there is no universal TE detector tool, and most softwares are biased toward specific questions (i.e., masking of repeats or identification of a certain class or family of TE by genome blasting). On the other hand, there is no universal database for finding detailed information on the TE families present in a given genome. Repbase is a fairly complete collection of repetitive elements in several eukaryotic genomes [37]; however, it only presents consensus sequences, and the individual sequences used to generate those consensuses are not available. Consequently, there is a lack of reports dealing with the general features of the whole repertoire of TEs in certain genomes, even if there are excellent exceptions as in Lerat et al. 2003 [17]. It is due to this lack of information that we are missing a complete picture of the characteristics of the TE families in An. gambiae, especially in relation to the sequence diversity within given families.

In an attempt to describe the diversity of elements and their deteriorated patterns, we have first compiled all the TEs so far described in the An. gambiae genome and used several families of elements from AnoTExcel—a previously characterized database [38] presenting detailed information on the TEs in An. gambiae[27]—to characterize their degree and patterns of deterioration.

TE in An. gambiae

The TEs in the An. gambiae genome are represented by members of both classes (I and II) with approximately 25 different superfamilies, each composed of several families (127, 121, and 123 for the LTRs, NLTRs, and Class II, respectively) (Figures 1 and 2). They constitute between 12% and 16% of the mosquito genome [39,40]. According to the preliminary genomic analysis [40], LTR elements were the most abundant TEs, followed by SINEs (non-autonomous NLTRs) and MITEs (non-autonomous class II elements) [40]. We have now scrutinized all the TE elements deposited in Repbase [37], TEfam [41], and the publications reporting TEs in An. gambiae so far [27,42-51]. In accordance with previous publications [39,52], we found that NLTRs are the most abundant type of elements in this organism, in terms of both diversity (number of superfamilies) and abundance (number of individual sequences), followed by elements belonging to class II ( 1 and Figure 1); however, the proportion of full-length and truncated elements varies considerably for each class. The former corresponds to 46, 1, and 6% for LTR, NLTR, and class II, respectively (Figure 2). The higher proportion of full-length elements in the Class I, LTR order in comparison to NLTRs and Class II could have important consequences in the dynamic of this particular genome. The activity of LTR elements could have an impact on expression of mosquito genomes, but in addition, these elements could produce important structural consequences through recombination among their LTRs; however, a great part of the LTR elements here analyzed correspond to Solo-LTRs (43%; Figure 2a). Within the LTR order, the Gypsy superfamily has the majority of the families and the higher diversity of elements (Figure 2b) [43]. On the other hand, there are only ten Copia superfamilies in the mosquito genome, and they contain a similar number of full-length and Solo-LTR sequences. The Copia families are composed of very few sequences presenting high degrees of nucleotide identity, indicating that they were active until recently. The Pao-Bel elements have an intermediate number of families and are mainly represented by full-length sequences. As a whole, the LTR elements present a low proportion of fragmented sequences together with a high proportion of Solo LTRs; this seems to indicate that the main mechanism driving the deterioration of LTR elements is through the generation of Solo elements, probably driven by LTR recombination as it was previously described for the Ty3/Gypsy elements [52].

Additional file 1. Tabl S1. Catalog of transposable elements described inAnopheles gambiaeto date. There are three worksheets in the excel file named “LTRs”, “NLTRs,” and “CLASS II”. Each contains information regarding the main TE classes and orders of the respective elements. The source where the specific family was described for the first time is presented. Each worksheet contains a table describing the number of full-length, truncated sequences, Solo-LTRs, and MITEs (for the LTR and class II, respectively) and the total number of families for each superfamily.

Format: XLS Size: 154KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 1. Distribution of transposable elements in the genome ofAnopheles gambiae. Data were compiled from Repbase (RB), Tefam (TF), and a series of publications describing transposable elements in the mosquito genome.

thumbnailFigure 2. A) Proportion of full-length, truncated, Solo-LTR, and MITE-like sequences for i) LTRs, ii) non-LTRs and iii) class II elements in theAn. gambiaegenome.B) Copy number of full-length, truncated, Solo-LTR, and MITE-like sequences for each superfamily analyzed. Data were collected from Repbase, Tefam, and AnoTExcel. Obs: the y-axis for Figure b ii and iii are logarithmic.

On the other hand, NLTRs form a diverse group with all the families presenting heavily deteriorated sequences constituting up to 99% of this order, with some superfamilies presenting only truncated sequences (Figure 2a,b). The retrotransposition mechanism of NLTRs commonly involves generation of 5′ truncated sequences [53-55] by a mechanism that is not clearly understood. These truncated sequences are integrated in the genome and are believed to evolve neutrally [54] and have been designated “dead-on-arrival” (DOA). The mechanism for generation of 5′-truncated sequences seems not to be universal for all NLTR superfamilies, as they are not present in all the species where NLTRs have been described. A great part of the total amount of NLTR elements in this genome is due to the presence of thousands of SINE elements belonging to the SINEX-1 family [56].

The class II elements are also heavily deteriorated (94% correspond to truncated sequences or MITEs) (Figure 2a,b). A great part of them belongs to MITEs (60%), or non-autonomous elements, with several internal deletions and with no coding capacity but able to be mobilized by the action of active elements in trans. This means that the truncated class II elements are still participating—or at least are able to do so—in the transposition dynamics of the families that they belong to.

Deterioration of TEs in An. gambiae

TEs are expected to diverge from their original sequence both in nucleotide composition and in structure during their genomic evolution; however, elements belonging to different classes and orders incorporate errors at different rates and in different ways according to their own mechanisms of replication and to the enzymes involved in their replication. These errors would be lost or fixed, leaving a different landscape of deterioration of the several TE families in different genomes.

The TE content of a genome contributes significantly to differences in the amount of genomic DNA between phylogenetically close species. For instance, the TE content of Aedes aegyptiCulex quinquefasciatus, and An. gambiae differs significantly with DNA genomic size [39]. On the other hand, the great part of the TEs identified in today’s genomes, including An. gambiae, corresponds to deteriorated sequences or remnants of once active elements [57]. It has been previously described [35] that the TEs in Drosophila melanogaster and An. gambiae are, on average, 12 and 24% of the length of their full-length counterparts, respectively, indicating that these elements suffer an important degree of deterioration through deletions in their life cycle. The deterioration process of TEs drives their sequences to a loss of identity, which can be guided by the accumulation of deleterious mutations, incorporation of indels generating frame shifts, or deletion of longer regions in their coding sequences or recognition sites (segmental deletions).

TEs belonging to different classes are known to deteriorate in different ways, in accordance with their different structures and mechanisms of amplification. Class I LTR elements form Solo-LTRs as a by-product of recombination between flanking LTRs of the same element, which is actually believed to be the main force driving the deterioration process of these elements [52,58,59]. The other elements belonging to this class, the NLTRs, have a completely different mechanism of transposition involving a target site-primed reverse transcription (TPRT), which is believed to be related to generation of DOA elements, i.e., sequences lacking different sizes of their 5′ ends [60,61]. Nevertheless, generation and maintenance of 5′-deleted NLTR sequences are not understood completely, and it might be possible that other mechanisms are involved—for instance, previous deterioration of the sequences that are lost during replication, or the presence of common motifs that might serve as templates for (endo) nucleases. Last, class II elements tend to incorporate nucleotide substitutions and internal deletions, which in turn are believed to be related to the formation of MITEs, which are small, non-autonomous class II elements that maintain their terminal inverted repeats, and amplify by using the transposases of active elements in trans.

Here, in an effort to understand the patterns of deterioration of TEs in the genome of An. gambiae, we have analyzed several TE families belonging to the main superfamilies in the genome. We obtained nucleotide sequences belonging to different TE families from AnoTExcel [27]. The family sequences were further aligned to reference sequences from Repbase or TEfam and analyzed by the EMBOSS program Plotcon, which gives a graphic representation of the point by point similarities of sequences in family alignments. This qualitative analysis was complemented by the estimates of the p-distances among different regions of the full-length elements as well as the deteriorated sequences in each family, to present a picture of the main differences in the deterioration process followed by different elements.

Class I, LTR elements

We analyzed the full-length and Solo-LTR elements of several families belonging to the three LTR superfamilies: Pao-Bel, Copia, and Gypsy. The alignments corresponding to full-length sequences were divided into three regions—the 5′ LTR, the internal region, and the 3′ LTR—and p-distances were calculated for each of the three regions. To evaluate the relative age of the LTR elements, we also calculated the p-distances between the 5′ and 3′ LTR for each individual sequence in the alignment. The p-distances of Solo-LTR families were also calculated (i.e., group of Solo-LTRs that share more than 90% of total identity) (Figure 3a–c).

thumbnailFigure 3. Number of base substitutions per site (p-distance) from between sequences for the three LTR superfamilies inAnopheles gambiae. Analyses were conducted using the Kimura 2-parameter model. Rate variation among sites was modeled with a gamma distribution (shape parameter = 1). All ambiguous positions were removed for each sequence pair. a) 5′-LTR (population comparison); b) INTERNAL REGION (population comparison); c) 3′-LTR (population comparison); d) 5′ to 3′-LTR comparison within each full-length sequence within the population; and e) SOLO-LTRs (only for Pao-BELs and Gypsy superfamilies). A) Pao-Bel superfamily (12 full-length and 3 Solo-LTR families). Kruskal-Wallis test: p value < 0.001. Ad hoc, Tukey non-parametric test: all comparisons were significant at p < 0.05 except a–c. B) Copia superfamily (four families). Kruskal-Wallis test: p-value = ns. Not significant differences among the regions and (C) Gypsy superfamily (six full-length and eight Solo-LTR families). Kruskal-Wallis test: p-value < 0.001. Ad hoc, Tukey non-parametric test: all comparisons were significant at p < 0.05, except a–b, a–c, and b–c with p-value > 0.05.

Notably, at the population level, comparison of the 5′ and 3′ LTRs shows higher p-distances than the internal regions for the three LTR superfamilies (Figure 3a–c; boxplots a and c versus b). Although these differences are significant only for the elements in the Pao-Bel superfamily, dispersion of the data for LTR regions (a and c) in the Copia and Gypsy superfamilies is also larger than the distances in the internal regions. It has been previously shown that the noncoding LTR-ULR region, despite its regulatory functional importance, is the most rapidly evolving region of LTR retrotransposons [62]. A similar observation indicating different levels of nucleotide substitutions in different genes was reported for the retroviruses HTLV type I and II [63], suggesting a higher heterogeneity of the LTR regions.

Interestingly, comparison of the 5′ and 3′ LTRs within each full-length sequence shows that these regions are almost identical in all the cases, with a slightly higher dispersion for the Copia elements (Figure 3a–c; boxplots d). This is also evident in the deterioration profiles of these sequences shown in 2, i.e., even if the sequences present heterogeneity in the LTRs at the population level, the deterioration level is quite symmetrical within each individual sequence.

Additional file 2. Figure S1. Graphic representation of the deterioration profiles for the LTR superfamilies and families analyzed in this study. The graphics performed with Plotcon (http://emboss.bioinformatics.nl/cgi-bin/emboss/plotcon) represent the similarity of the sequences along multiple alignments performed for several families belonging to the LTR superfamilies: Pao-Bel: (al) correspond to families of full-length sequences, (mo) correspond to families of Solo-LTR sequences; Copia superfamily (ps) correspond to families of full-length sequences; and the Gypsy superfamily: (ty) correspond to families of full-length sequences, (zac) correspond to families of Solo-LTRs. Each family was aligned to the respective reference sequence as described in Repbase. The family name and the cluster number according to the AnoTExcel numbering (in parenthesis), are indicated in each graph as. Red circles are included to highlight regions where the differences between sequences are mainly due to indels; blue circles indicate regions where the differences are due to nucleotide substitutions; and green circles indicate regions where segmental deletions are present. Red asterisks are included to indicate single deletions and blue asterisks single insertions. The horizontal green bars represent the ORFs in the alignments, with numbers indicating their relative position to the first nucleotide of the alignment. Green arrows at the top of each graph indicated the relative position of the LTRs; the numbers within the arrows indicate the length of the LTRs. The X axis for all plots refers to the relative residue position in each alignment and the Y axis to their similarity indicated as the pairwise scores that are taken from the specified similarity matrix (see Methods section for detailed information). (PDF 331 kb)

Format: PDF Size: 332KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

All of the families containing full-length sequences (with the exception of Gypsy3_AG) ( 2) also contain open reading frames. These previous observations together indicate that these families might still be active or have been recently transposed.

On the other hand, as expected, the families of Solo-LTR sequences belonging to the Pao-Bel and Gypsy superfamilies present highly deteriorated sequences, presenting both nucleotide substitutions and indels or segmental deletions along their alignments ( 2), indicating their neutral way of evolution.

It has been generally assumed that after transposition, the new copies of LTR elements are full length in size and identical in nucleotide composition to the source sequences and that during evolution, it is expected that these sequences will diverge in sequences and structure due to Solo-LTR formation and indel accumulation [52]. In fact, the level of divergence estimated between copies and consensus or canonical sequences has been used to estimate the ages of insertions in several genomes [52][64]. Recently, Arensburger et al. [39] performed a comparative analysis of the TEs in three mosquitoes genomes, Aedes aegyptiAn. gambiae, and Culex quinquefasciatus, and used the divergences among the members of different families to compare the relative ages of TEs in the three genomes. They found that both LTR and NLTR retroelements dominated the most recent relative age classes, consistent with the presence of recently active retrotransposons and a gradual degradation of the sequences. Our data showed that alignments of Pao-Bel families present greater heterogeneity than the Copia or the Gypsy families. This heterogeneity is mainly due to the presence of nucleotide substitutions along their alignments with most of the point or segmental deletions being present in the flanking regions. The Copia families appeared to be the most conservative sequences, which together with the absence of Solo_LTRs sequences in this superfamily could indicate a more recent introduction of this superfamily in the Anopheles genome. Finally, most of the Gypsy families here analyzed present quite homogeneous sequences, with most of the differences between the sequences being due to the presence of nucleotide substitutions. These data are in agreement with a complete analysis of the Ty3-Gypsy [52], which suggested that the main mechanism driving the evolution of Gypsy elements is the formation of Solo LTRs, which in turn must be subjected to lower selective pressure than the full proviral sequences and therefore persisting longer in the genome, allowing for accumulation of mutations and deletions.

The absence of a spectrum of divergence of class I elements and a relatively high homogeneity of LTR sequences found in the genome of Saccharomyces cerevisiae has been attributed to a rapid turnover of copies that become inactivated in that genome by LTR-LTR recombination leading to the formation of Solo-LTRs [17]. On the contrary, Solo-LTRs are mostly absent from the genome of D. melanogaster[17], so the LTR degradation depends on different mechanisms.

In the An. gambiae genome, the LTR elements are represented by members of the three main superfamilies, which present slightly different characteristics. The Copia elements are represented by families of elements with few copies and high homogeneity among their sequences, besides presenting few Solo-LTRs. These features indicated that these elements are young and active in this genome. On the other hand, both the Pao-Bel and Gypsy families here analyzed have quite heterogeneous sequences, even for putative active elements, presenting several nucleotide substitutions and deletions that are far more representative than insertions, as has been previously described for Gypsy elements [52]. Most of the deletions in these sequences are outside the ORFs or inside a few sequences in the alignment, indicating that for the group of sequences analyzed here, the Solo-LTR is a more important source of LTR-family degradation than the deletion of their sequences.

Class I, NLTR elements

Five different superfamilies of NLTR elements were analyzed (CR1, I, Jockey, Outcast, and RTE). The profile of deterioration of the CR1 and RTE superfamilies presents the previously described 5′ truncated sequences [54,65-67], giving a stair-like pattern to the sequence alignments ( 3), while Jockey, I, and Outcast are represented by full-length sequences with no significant differences in the frequencies of indels or nucleotide substitutions along the full-length alignments ( 3). Most of the families analyzed belong to the CR1 superfamily, which is the most abundant superfamily in the mosquito genome.

Additional file 3. Figure S2. Graphic representation of the deterioration profiles for the non-LTR families analyzed in this study. See legend for Additional file 1: Figure S1 for detailed information. (a-e) CR1 families composed of full-length and 5′ truncated sequenced; (fm) CR1 families composed of 5′ truncated sequences; (n,o) CR1 families containing 3′ truncated sequences; (p,q) families of RTE elements composed of full-length and 5′ truncated sequences, (r,s) Jockey families composed of full-length sequences, (t) family of I full-length sequences, (u,v) families of Outcast elements composed of full-length sequences. Colored arrows indicate the relative positions of ORF1 (blue), ORF2 (red), RT domain (yellow), and 3′ conserved region (green). The relative positions according to multiple alignments of the sequences to canonical element are indicated in each region. The number of sequences represented in each region is indicated above the arrows. (PDF 350 kb)

Format: PDF Size: 351KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Analysis of the similarities at the nucleotide level of the sequences in the families of NLTR elements is shown in Figure 4. The p-distances of the families containing full-length sequences corresponding to the Jockey and RTE superfamilies were significantly smaller than their truncated counterparts (Figure 4d,e and h,i), and the Outcast families have very low distances both in their full-length and truncated sequences (Figure 4f, g); however, there were no significant differences between families containing full-length and truncated sequences and those containing only truncated sequences from the CR1 superfamily (Figure 4a, b). CR1 families present a very broad range of p-distances both among full-length and truncated sequences. We have analyzed fifteen CR1 families: five presenting full-length and 5′-truncated sequences ( 3), eight with only 5′-truncated sequences of different lengths ( 3), and two composed only of 3′deleted sequences ( 3). Some of the CR1 families containing both full-length and truncated sequences have large p-distances (more than 6%) (Figure 4a), which might indicate the process of deterioration of the truncated sequences that are also present in the families and included in this analysis. In fact, the deterioration profile of those sequences ( 3) shows that some of the 5′-truncated sequences from each of the CR1 families analyzed are quite deteriorated, containing nucleotide substitutions, insertions, and deletions. Still, all the full-length sequences within these families contain ORFs and no indels at all; however, closer inspection of the individual sequences in these families showed that several of the 5′ truncated sequences keep ORFs covering the whole ORF2, the RT domain, or even a conserved region in the 3′ of ORF2 with unknown function ( 3). In some cases, an elevated number of truncated sequences maintain their ORFs (see 3). On the other hand, some families containing only truncated sequences (Figure 4b) present low p-distances, which might indicate a constraint in the evolution of these truncated sequences. A detailed analysis of truncated sequences in these families also showed that several of the truncated sequences preserved ORFs comprising the regions mentioned before. It is interesting to note that most of the insertions or deletions in the deterioration profiles shown in 3 are present in the 3′ region after the end of ORF2. This is a very unusual pattern, because most commonly the truncated sequences of CR1 elements lack ORFs at all [68].

thumbnailFigure 4. Number of base substitutions per site (p-distance) from between sequences belonging to five superfamilies of NLTR elements. See Figure 3 for detailed information. a) and b) CR1 superfamily: full length (four families) and fragmented sequences (16 families), respectively; c) I elements (one full length family); d) and e) Jockey superfamily, full length (three families) and fragmented sequences (two families), respectively; f) and g) Outcast superfamily, full length (two families) and fragmented sequences (one family), respectively; h) and i) RTE superfamily, full length (two families) and fragmented sequences (one family), respectively. We tested the significance of the difference between families containing full-length sequences and those containing only truncated sequences. Kruskal-Wallis test: p-value < 0,001. Ad hoc, Tukey non-parametric test. p-values: a–b = ns; d–e**; h–i*. P-value: ns = not significant; * < 0.05; ** < 0.001.

To better understand the possible significance of the above findings, we further analyzed the CR1 families and classified them into three groups according to the level of sequence deterioration in each family as i) only full-length, ii) truncated sequences belonging to families where full-length sequences are also present, and iii) families composed only of truncated sequences. The multiple sequence alignments analyzed were divided into four regions comprising ORF1, ORF2, the RT domain, and the 3′ conserved region. Analysis of the similarities among the different regions of the full-length and truncated sequences is shown in Figure 5. When considering only the full-length sequences, ORF1 presents significantly less variability than ORF2, including the RT domain, but surprisingly have no significant differences with the 3′ conserved region of ORF2 (Figure 5). This conserved region within ORF2 has no known function, although it appears to be much conserved in NLTRs. There were no significant differences in any of the regions for the two groups of truncated sequences analyzed.

thumbnailFigure 5. Number of base substitutions per site (p-distance) from sequences belonging to CR1 elements. See Figure 3 for detailed information. Lines ad represent the p-distances among different regions (ORF1, ORF2, RT, 3′-conserved region) of full-length sequences; lines eh represent the p-distances from truncated sequences belonging to the same families analyzed in lines ad (truncated I); and lines il represent distances from truncated sequences from families were no full-length sequences were identified (truncated II). The mean p-distances and standard deviation are shown on top of each boxplot. Kruskal-Wallis test: p-value < 0.001. Ad hoc, Tukey non-parametric test: for full length sequences: a–b**, a–c**, d–c*; the rest are not significant. For ORF1 comparisons: a–e***; a–i***; e–I = ns; for ORF2: b–f = *; b–j = *; f–j = ns; for RT: c–g = **; c–k = **; g–k = *; for 3′ conserved region: d–h = ***; d–l = ***; h–l = ***. p-value: ns = not significant; * < 0.05; ** < 0.005, *** < 0.0001.

When comparing the different regions for the three groups of sequences, both ORF1 and ORF2 presented significantly lower distances in the full-length than in the truncated sequences from both groups. On the other hand, there were no significant differences in ORFs 1 and 2 for the truncated sequences. For the RT region and the 3′ conserved region, the full-length sequences presented significant lower distances than their truncated counterparts, but in both cases, the truncated sequences belonging to families with full-length sequences presented significantly lower distances among their sequences than the other group of truncated sequences.

We further analyzed the truncated sequences corresponding to RT and 3′ conserved regions of two families of CR1 elements that contain complete ORFs: CR1-Ele7(clu36) and CR1-Ele13(clu16) contain sequences covering only the complete RT domain and the “3-common” region. Two truncated sequences from family CR1-Ele7, both located in chromosome 2L, are very similar, presenting distances between them both in the RT and in the 3-common region, 0.00539 and 0.00335, respectively (Figure 6a). Five sequences of the CR1-Ele13(clu16) family present a very small distance among them only in the RT region, 0.00743, while the distance of the 3-common region is one order of magnitude higher (Figure 6b). The high similarity among these sequences is surprising, as they are not part of full-length elements. They are all located in different chromosomes. These observations are compatible either with a recent origin of these sequences by an incomplete transcription of full-length elements or with a direct constraint in the evolution of these sequences, such as domestication of the element. We further analyzed expression of these sequences by comparing them by blastn to expression libraries, but they did not show any significant result. It has been previously suggested that genomes might have recurrently recruited TE-derived enzymatic or structural functions for their benefit [69]. In the case of the CR1 families here analyzed, it is surprising that ORF preservation in truncated sequences appears in most of the families. It has been previously noted that DOA elements present low nucleotide polymorphism together with a large number of internal deletions, a fact that has been related to a high rate of DNA loss in Drosophila[70]. On the other hand, other NLTR elements previously analyzed as t F, docs, or Jockey elements do not present the same pattern of DNA loss [17].

thumbnailFigure 6. Schematic representation of truncated sequences spanning the RT domain and 3′ conserved regions for two CR1 families: A) family CR1-ele7 (cluster 36) and B) family CR1-ele13 (cluster 16). Top sequences are canonical sequences for reference. The numbers in white represent the relative positions of each truncated fragment. At the bottom of each panel, the p-distances for each of the regions indicated are shown. The exact chromosomal position of each truncated sequence is shown at the left.

We also found two families of CR1 elements (CR1-2 (clu46) and CR1-3 (clu17) presenting only 3′ deleted sequences, which is a very unusual pattern for NLTR elements ( 3). The CR1-2_AG (clu46) family is composed of 20 sequences that were aligned to the respective canonical sequence in Repbase. Twelve of these sequences consist of three blocks of homology with canonical CR1-2_AG covering positions 1–1762, 2150–2600, and 4511–4674; three of the sequences cover the positions 2150–2600 and 4511–4674, and six sequences cover the first 545 nucleotides of the canonical element. The overall similarity of the sequences is high (p-dist = 0.0141, sd = 0.0012) and they present few indels. Eight of these truncated sequences present full-length ORF1 sequences, with p-distances among them of 0.0114 (sd = 0.0015). The 50 sequences belonging to the CR1-3_AG(clu17) family, on the other hand, correspond to the first 820 nucleotides of the canonical sequence, which is 5515 nucleotides long and contain several insertions, deletions, and point substitutions that are evident in the Plotcon analysis ( 3). A different mechanism than that proposed for the generation of 5′-deleted elements of the NLTR order (DOA)—involving the reverse transcription of these elements—needs to be invoked to explain this finding. Analyses of diverse NLTRs in several species have shown that the 5′ deletion is common among these elements [71-73]. The reverse transcription starts at the 3′ end and is believed to fail to proceed to the 5′ end, generating a large number of copies with varying lengths in the genome [54,66,67]. On the other hand, complete absence of 5′-deleted NLTRs in other species has also been described. In the rotifer Adineta vaga, NLTRs have been found to be inactivated by internal deletions, while no 5′ truncated elements were found [74]. A mechanism related to the target site-primed reverse transcription of NLTR elements is normally used to explain the 5′-deleted elements; however, considering the total absence of these types of truncated elements in other NLTR families, together with the fact that reverse transcription is the mechanism used to replicate these families—which in many species are shown to be extremely successful families (e.g., Alu sequences in humans)—indicate the possible presence of a different, yet unknown, mechanism involved in this process.

Class II

Class II elements clearly represent the most heterogeneous set of sequences of all the TE families in An. gambiae. For all the families analyzed, the profiles of deterioration show heavily deleted sequences ( 4), and the distances among their sequences present broad differences (Figure 7a–h). Most of the families are represented by putative MITE elements, i.e., groups of sequences sharing TIRs with known or even unknown class II elements that are amplified by the transposases of active elements in the genome.

Additional file 4. Figure S3. Graphic representation of the deterioration profiles for the class II families analyzed in this study. See legend for Additional file 1: Figure S1 for detailed information. The blue arrows at the top of each graph represent the position and length of the TIRs in the canonical full-length element. The red arrows indicate the actual TIRs in the MITE-like elements. An “I” after the number indicates incomplete TIRs. Note the different scales on the y-axis. (PDF 201 kb)

Format: PDF Size: 201KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

thumbnailFigure 7. Number of base substitutions per site (p-distance) from between sequences belonging to eight superfamilies of class II elements. See Figure 3 for detailed information. A) Gambol (seven families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: p-value = ** for all comparisons except b–d that was ns. B) Harbinger (three families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: p-value a–b = **; a–c = **; a–d = ns; b–c = ns; b–d = ns; c–d = **. C) hAT (four families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: all comparisons were significant at p < 0.001 except a–c = ns. D) Ikirara (two families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: a–b = **; a–c = ns; a–d = ns; b–c = **; b–d = ns, c–d = ns. E) Mariner (seven families). Kruskal-Wallis test: p-value = **, Ad hoc, Tukey non-parametric test: integer overflow. F) P (five families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: integer overflow. G) Tc1 (six families). Kruskal-Wallis test: p-value = ** Ad hoc, Tukey non-parametric test: integer overflow. H) Transib (Two families). Kruskal-Wallis test: p-value = **. Ad hoc, Tukey non-parametric test: all comparisons were significant at p < 0.001 except b–d = ns. P-value: ns = not significant; * < 0.05; ** < 0.001. a) 5′-TIR (population comparison), b) internal region (population comparison), c) 3′-TIR (population comparison), d) 5′-TIR to 3′-TIR within element (performed for each element within each of the populations analyzed). On top of each boxplot, the mean and standard deviation for the p-distances in each group of sequences are shown.

Some of the MITE-like sequences analyzed share TIRs with recognizable elements present in the mosquito genome, as is the case of the Gambol elements, hAT and Harbinger; however, for some families—although they have a MITE-like structure due to their small size, absence of coding capacity, and presence of TIRs—there is no obvious full-length counterpart described in the An. gambiae genome so far. This is the case for all the marinerN elements harbingerN1 and N2, hatN1, as well as of a series of MITE-like described in the genome that do not show any similarity with known class II elements in their TIRs or their internal regions [27].

All the Gambol elements analyzed correspond to MITE-like sequences, which have been previously described [27]. Interestingly, they present shorter TIRs than the canonical elements, which in one case have different sizes (Gambol_ele3). The TIRs of the canonical Gambol families are larger than typical TIRs of class II elements, while the TIRs of the MITE like sequences are shorter. This might indicate that a smaller portion of the TIR is necessary for proper recognition by the active transposase. On the other hand, TIRs are very similar at the population level with slightly higher distances for the internal regions (Figure 7A:a–c). The 5′ to 3′ TIR comparison reveals differences at the nucleotide composition of the flanking TIRs in certain sequences (Figure 7A:d).

The Harbinger and hAT elements also present MITE-like structures with families presenting imperfect TIRs, which is revealed by the larger p-distances in the 3-5′ nucleotide similarities comparison (Figure 7B:d and C:d). The 5′and 3′ TIRs for the hAT families are identical at the population level, favoring the idea that they might have transposed recently.

The Ikirara, Mariner, and P elements (Figure 7D–F) show larger similarity distances in their 5′ and 3′ TIRs. The Mariner sequences seem to be composed of MITE-like families and have been previously characterized as MarinerN families (Repbase). There are 21 sequences of this type of element described in An. gambiae (Repbase reports). Most of them contain imperfect TIRs; two of the families present sequences lacking either the 3′ or the 5′ region and therefore lack one TIR, which probably makes them dead elements. They have been compared to the MarinerN elements to which they belong ( 4), as they do not share identical TIRs with any full-length Mariner. The P MITE elements [35] have also been compared to the MITE sequences already described as such ( 4). The Tc1 families present heterogeneous 5′ TIRs (Figure 7G:a) but conserved 3′ TIRs (Figure 7G:c) at the population level, and both TIRs are quite different within the same sequences (Figure 7G:d). Finally, the Transib elements (Figure 7H) have almost identical 5′and 3′ TIRs at the population level, although they have imperfect TIRs ( 4).

Some of the families analyzed are quite heterogeneous in relation to the canonical sequences of the families to which they belong, giving very irregular profiles on the Plotcon analysis ( 4), but they are quite homogeneous within their own families. This probably means that they constitute subfamilies of MITEs that are different from the sequences used to create the consensus sequences deposited in Repbase.

Finally, we have based our analysis on the families of deteriorated elements that we found in our previous study and that are present in AnoTExcel [27]. We did not find any member of the Herves family and therefore did not include it in the present analysis.

Conclusions

We described an in silico analysis of elements representing the main superfamilies of TEs identified in the genome of An. gambiae. Our study concentrates on the deterioration pattern of the TEs in this genome and is based on families presenting both full-length and truncated sequences.

It is not an easy task to represent the deterioration profile of sequences that are repetitive and that have undergone nucleotide substitutions, insertions, deletions, and recombinations during their genomic evolution. Although there have been previous works on the well characterized and model organism, D. melanogaster[17,64], there are to our knowledge few reports that have presented a description of the different structural variations and nucleotide polymorphisms of the various TEs in other genomes. This information is of great importance when trying to understand how these elements replicate, amplify, degrade, and are eventually eliminated from a genome.

One such study is that of Lerat et al. [17] that presented an analysis of sequence divergence of 23 transposable element families in the D. melanogaster genome. These authors found a high degree of homogeneity and a lack of divergent elements between sequences of TEs within a given family. On the other hand, they found divergent elements displaying very low percentage of similarities to the full-length sequences. These findings suggested that TEs are highly active within the genome and that the highly divergent copies reflect relics that have been degenerated and rearranged. In the An. gambiae genome, several reports have described the evolution and dynamics of different superfamilies of TEs (e.g., [27,42,43,52]); however, as far as we are concerned, there is no previous report dealing with TEs belonging to all the main superfamilies in this genome.

Here we present the deterioration landscape of most of the TE superfamilies described in An. gambiae. We chose to show the patterns of deterioration as a graphic representation of the different elements present in this particular genome. This approach has not been used before, and we believe that it helps to reveal the overall differences identified in the TE families of the An. gambiae genome.

In any given genome, TEs might present significant diversity regarding the level of sequence deterioration of a given family. The relationship between TEs and the genome where they reside produces particular conditions for evolution of different TE families, and the success or failure in settling and amplifying within a genome would be the result of interactions and evolutionary dynamics between the element and the specific genome. What different strategies can be adopted by different classes of TEs—and what evolutionary forces are involved in this process— are still unanswered questions.

We have not found truncated LTR elements other than the Solo-LTRs; on the other hand, a great part of the NLTRs appears as truncated sequences. In addition, most of the class II are MITE-like elements with highly degraded copies. Overall, we found less divergence for class I (both LTR and NLTR families) than for class II families. A similar finding was described in the Drosophila genome [17]. This result is somewhat surprising, as it is well known that replication based on reverse transcriptases (as is the case of retroelements and retroviruses as well) are quite error prone mechanisms [75,76]. Indels are present in several LTR families and are evident in the deterioration profiles in 2, tending to be more frequent at the 5′ and 3′ ends of the elements. The absence of divergence of the LTR copies could be interpreted as a rapid turnover of the elements once they have been inactivated.

For NLTRs, besides the 5′ truncated sequences previously described, the elements tend to be quite similar within families. They contain indels and point mutations throughout their lengths, but the internal deletions do not seem to be a driving force in their deterioration process as has been previously observed in Drosophila[77-79]. It is also interesting to note the high frequency of truncated sequences preserving ORFs. The presence of truncated sequences belonging to biologically important domains—such as the reverse transcriptase in the NLTR elements with very small nucleotide distances among some of the sequences—is intriguing. One would expect that the truncated sequences would have the same time of evolution than other truncated sequences of the same family, and therefore, according to a molecular clock, these neutral sequences might have accumulated mutations at more or less the same rate. This is not the case for families CR1-ele7(clu36) and ele13(clu16). It might well be that they have been domesticated by the host genome in a certain manner.

On the other hand, the distances among the sequences of the same families for the class II elements tend to be larger than the NLTRs and LTRs. The high sequence variation and the indels between different copies of the same family, mainly for the MITE-like elements, indicate they are ancient sequences in the An. gambiae genome.

Although it is possible to compare the relative ages of TEs belonging to the same family found in different genomes based on their level of nucleotide divergence [17,39,52], comparison of ages of elements belonging to different orders/classes or even families is more complicated, because the mechanisms responsible for their replication involves quite different processes. For instance, it is well known that the replication of retroviruses is a very error-prone process due in particular to the lack of proof reading repair activity of the RNA polymerase and reverse transcriptase enzymes (review in [80]). As retrotransposition of class I elements resembles retroviral replication, it is reasonable to think that this process might also be error prone. In fact, the mutation rate during a single transposition cycle of the yeast Ty1 element has been estimated to be 2.5 × 10–5 substitution/nucleotide [81]—as high as that for retroviruses. We do not have data regarding mutation rates associated with different retrotransposons in An gambiae. Given the lack of information regarding the dynamics and rates of evolution of these elements, we consider it safer not to make comparisons on the relative ages of the different classes of elements in this genome.

In summary, we show here that the transposable elements in the An. gambiae genome deteriorate in different ways according to the class to which they belong. This diversity certainly has implications not only at the host genomic level, but also at the amplification dynamic and evolution of the TE families themselves.

Methods

Pipeline description

The TE families analyzed in this work were extracted from AnoTExcel, a TE-specific database from the An. gambiae genome [27]. Sequence alignments belonging to the main classes of TEs in An. gambiae were further analyzed. We used families containing elements at different deterioration levels. We analyzed 33 families of LTR elements (12 Pao-Bel with full length sequences, and 3 with Solo-LTRs; 4 Copia with full-length sequences; and 14 Gypsy families—6 with full-length elements and 8 with Solo-LTRs), 24 families of NLTRs (16 CR1; 1 I; 3 Jockey; 2 Outcast; and 2 RTE), and 21 class II (4 Gambol; 3 Harbinger; 2 hAT; 1 Tsessebeii; 5 Mariners; 2 P; 1 Tc1; and 2 Transib).

Sequence similarity (Plotcon)

Similarity along the sequences was calculated by moving a window of different lengths along the aligned sequences ([82]). Within the window, the similarity of each position is taken to be the average of all the possible pairwise scores of the bases or residues at that position. The pairwise scores are taken from the specified similarity matrix. The average of the position similarities within the window is plotted against the positions in the alignment. The average similarity is calculated as:

<a onClick="popup('http://www.biomedcentral.com/1471-2164/13/272/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/13/272/mathml/M1">View MathML</a>

(1)

sum, over column*window size. w, sequence weighting; M, matrix comparison table; i,j, with respect to residue i or j; Nseq, number of sequences in the alignment; Wsize, window size.

p-Distances

The number of base substitutions per site between sequences were calculated using MEGA 5.0 [83]. Analyses were conducted using the Kimura 2-parameter model. Rate variation among sites was modeled with a gamma distribution (shape parameter = 1). All ambiguous positions were removed for each sequence pair.

Statistics

Data analysis was carried out with the R System (R version 2.14) [84]. We used non-parametric tests Kruskal-Wallis [85] and non-parametric ad-hoc Tukey test [86]. Differences were considered significant at p values < 0.05.

Abbreviations

Clu: Cluster; DOA: Dead-on-arrival; EST: Expressed sequence tag; LTR: Long terminal repeat; Mb: Megabase; MITE: Miniature inverted TE; NLTR: Non-LTR element; ORF: Open reading frame; RB: Repbase; RT: Reverse transcriptase; TE: Transposable element; TF: TEfam; TIR: Terminal inverted repeat; TRIM: Terminal-repeat retrotransposons in miniature.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RDFM conceived the report, participated in its design, performed data analysis, interpreted results, and drafted the manuscript. CMAC, JMR, and CJS participated in intellectual discussion and manuscript drafting. JMR: helped with bioinformatics. LV helped with statistics. All authors read and approved the final manuscript.

Acknowledgements

This work was supported in part by the Intramural Research Program of the Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health. We thank Brenda Rae Marshall, DPSS, NIAID, for editing. RDFM and CJS were funded by FAPERJ (Fundação de Amparo a Pesquisa do Estado de Rio de Janeiro) and CNPq (Conselho Nacional de Pesquisa). CMCA was funded by CNPq.

Because JMC Ribeiro is a government employee and this is a government work, the work is in the public domain in the United States. Notwithstanding any other agreements, the NIH reserves the right to provide the work to PubMedCentral for display and use by the public, and PubMedCentral may tag or modify the work consistent with its customary practices. You can establish rights outside of the U.S. subject to a government use license.

References

  1. Werren JH: Selfish genetic elements, genetic conflict, and evolutionary innovation.

    Proc Natl Acad Sci U S A 2011, 108(Suppl 2):10863-10870. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Britten RJ: DNA sequence insertion and evolutionary variation in gene regulation.

    Proc Natl Acad Sci U S A 1996, 93:9374-9377. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Britten RJ: Mobile elements inserted in the distant past have taken on important functions.

    Gene 1997, 205:177-182. PubMed Abstract | Publisher Full Text OpenURL

  4. Jordan IK, Rogozin IB, Glazko GV, Koonin EV: Origin of a substantial fraction of human regulatory sequences from transposable elements.

    Trends Genet 2003, 19:68-72. PubMed Abstract | Publisher Full Text OpenURL

  5. van de Lagemaat LN, Landry JR, Mager DL, Medstrand P: Transposable elements in mammals promote regulatory variation and diversification of genes with specialized functions.

    Trends Genet 2003, 19:530-536. PubMed Abstract | Publisher Full Text OpenURL

  6. Kidwell MG, Lisch DR: Perspective: transposable elements, parasitic DNA, and genome evolution.

    Evolution 2001, 55:1-24. PubMed Abstract OpenURL

  7. Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O, et al.: A unified classification system for eukaryotic transposable elements.

    Nat Rev Genet 2007, 8:973-982. PubMed Abstract | Publisher Full Text OpenURL

  8. Miller WJ, McDonald JF, Pinsker W: Molecular domestication of mobile elements.

    Genetica 1997, 100:261-270. PubMed Abstract | Publisher Full Text OpenURL

  9. Marino-Ramirez L, Lewis KC, Landsman D, Jordan IK: Transposable elements donate lineage-specific regulatory sequences to host genomes.

    Cytogenet Genome Res 2005, 110:333-341. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Polavarapu N, Marino-Ramirez L, Landsman D, McDonald JF, Jordan IK: Evolutionary rates and patterns for human transcription factor binding sites derived from repetitive DNA.

    BMC Genomics 2008, 9:226. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  11. Muotri AR, Marchetto MC, Coufal NG, Gage FH: The necessary junk: new functions for transposable elements.

    Hum Mol Genet 2007, 16(Spec No. 2):R159-R167. PubMed Abstract | Publisher Full Text OpenURL

  12. Piriyapongsa J, Marino-Ramirez L, Jordan IK: Origin and evolution of human microRNAs from transposable elements.

    Genetics 2007, 176:1323-1337. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Bureau TE, Ronald PC, Wessler SR: A computer-based systematic survey reveals the predominance of small inverted-repeat elements in wild-type rice genes.

    Proc Natl Acad Sci U S A 1996, 93:8524-8529. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Venner S, Feschotte C, Biemont C: Dynamics of transposable elements: towards a community ecology of the genome.

    Trends Genet 2009, 25:317-323. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Brookfield JF: Host-parasite relationships in the genome.

    BMC Biol 2011, 9:67. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Hua-Van A, Le Rouzic A, Boutin TS, Filee J, Capy P: The struggle for life of the genome’s selfish architects.

    Biol Direct 2011, 6:19. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  17. Lerat E, Rizzon C, Biemont C: Sequence divergence within transposable element families in the Drosophila melanogaster genome.

    Genome Res 2003, 13:1889-1896. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Charlesworth B, Langley CH: The evolution of self-regulated transposition of transposable elements.

    Genetics 1986, 112:359-383. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Lozovskaya ER, Hartl DL, Petrov DA: Genomic regulation of transposable elements in Drosophila.

    Curr Opin Genet Dev 1995, 5:768-773. PubMed Abstract | Publisher Full Text OpenURL

  20. Lippman Z, May B, Yordan C, Singer T, Martienssen R: Distinct mechanisms determine transposon inheritance and methylation via small interfering RNA and histone modification.

    PLoS Biol 2003, 1:E67. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Lisch D: Epigenetic regulation of transposable elements in plants.

    Annu Rev Plant Biol 2009, 60:43-66. PubMed Abstract | Publisher Full Text OpenURL

  22. Bestor TH, Bourc’his D: Transposon silencing and imprint establishment in mammalian germ cells.

    Cold Spring Harb Symp Quant Biol 2004, 69:381-387. PubMed Abstract | Publisher Full Text OpenURL

  23. Malone CD, Hannon GJ: Small RNAs as guardians of the genome.

    Cell 2009, 136:656-668. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Golden DE, Gerbasi VR, Sontheimer EJ: An inside job for siRNAs.

    Mol Cell 2008, 31:309-312. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Chiu YL, Greene WC: The APOBEC3 cytidine deaminases: an innate defensive network opposing exogenous retroviruses and endogenous retroelements.

    Annu Rev Immunol 2008, 26:317-353. PubMed Abstract | Publisher Full Text OpenURL

  26. Levin HL, Moran JV: Dynamic interactions between transposable elements and their hosts.

    Nat Rev Genet 2011, 12:615-627. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Fernandez-Medina RD, Struchiner CJ, Ribeiro JM: Novel transposable elements from Anopheles gambiae.

    BMC Genomics 2011, 12:260. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Struchiner CJ, Massad E, Tu Z, Ribeiro JM: The tempo and mode of evolution of transposable elements as revealed by molecular phylogenies reconstructed from mosquito genomes.

    Evolution 2009, 63:3136-3146. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Kennedy RC, Unger MF, Christley S, Collins FH, Madey GR: An automated homology-based approach for identifying transposable elements.

    BMC Bioinforma 2011, 12:130. BioMed Central Full Text OpenURL

  30. Britten R: Transposable elements have contributed to thousands of human proteins.

    Proc Natl Acad Sci U S A 2006, 103:1798-1803. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Han Y, Burnette JM, Wessler SR: TARGeT: a web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences.

    Nucleic Acids Res 2009, 37:e78. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Huda A, Jordan IK: Analysis of transposable element sequences using CENSOR and RepeatMasker.

    Methods Mol Biol 2009, 537:323-336. PubMed Abstract | Publisher Full Text OpenURL

  33. Meglecz E, Costedoat C, Dubut V, Gilles A, Malausa T, Pech N, Martin JF: QDD: a user-friendly program to select microsatellite markers and design primers from large sequencing projects.

    Bioinformatics 2009, 26:403-404. PubMed Abstract | Publisher Full Text OpenURL

  34. Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, et al.: The TIGR Rice Genome Annotation Resource: improvements and new features.

    Nucleic Acids Res 2007, 35:D883-D887. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Quesneville H, Nouaud D, Anxolabehere D: Detection of new transposable element families in Drosophila melanogaster and Anopheles gambiae genomes.

    J Mol Evol 2003, 57(Suppl 1):S50-S59. PubMed Abstract | Publisher Full Text OpenURL

  36. McCarthy EM, McDonald JF: LTR_STRUC: a novel search and identification program for LTR retrotransposons.

    Bioinformatics 2003, 19:362-367. PubMed Abstract | Publisher Full Text OpenURL

  37. Jurka J: Repbase update: a database and an electronic journal of repetitive elements.

    Trends Genet 2000, 16:418-420. PubMed Abstract | Publisher Full Text OpenURL

  38. AnoTExcel. , ;  .

    [(http://exon.niaid.nih.gov/transcriptome/TE/A_gambiae/AnoTExcel-WEB.zip]

    OpenURL

  39. Arensburger P, Megy K, Waterhouse RM, Abrudan J, Amedeo P, Antelo B, Bartholomay L, Bidwell S, Caler E, Camara F, et al.: Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics.

    Science 2010, 330:86-88. PubMed Abstract | Publisher Full Text OpenURL

  40. Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, et al.: The genome sequence of the malaria mosquito Anopheles gambiae.

    Science 2002, 298:129-149. PubMed Abstract | Publisher Full Text OpenURL

  41. Tefam. , ;  .

    [http://tefam.biochem.vt.edu/tefam/search.php]

    OpenURL

  42. Tubio JM, Costas JC, Naveira HF: Evolution of the mdg1 lineage of the Ty3/gypsy group of LTR retrotransposons in Anopheles gambiae.

    Gene 2004, 330:123-131. PubMed Abstract | Publisher Full Text OpenURL

  43. Tubio JM, Naveira H, Costas J: Structural and evolutionary analyses of the Ty3/gypsy group of LTR retrotransposons in the genome of Anopheles gambiae.

    Mol Biol Evol 2005, 22:29-39. PubMed Abstract | Publisher Full Text OpenURL

  44. Besansky NJ, Bedell JA, Benedict MQ, Mukabayire O, Hilfiker D, Collins FH: Cloning and characterization of the white gene from Anopheles gambiae.

    Insect Mol Biol 1995, 4:217-231. PubMed Abstract | Publisher Full Text OpenURL

  45. Besansky NJ, Mukabayire O, Bedell JA, Lusz H: Pegasus, a small terminal inverted repeat transposable element found in the white gene of Anopheles gambiae.

    Genetica 1996, 98:119-129. PubMed Abstract | Publisher Full Text OpenURL

  46. Biedler J, Tu Z: Non-LTR retrotransposons in the African malaria mosquito, Anopheles gambiae: unprecedented diversity and evidence of recent activity.

    Mol Biol Evol 2003, 20:1811-1825. PubMed Abstract | Publisher Full Text OpenURL

  47. Grossman GL, Cornel AJ, Rafferty CS, Robertson HM, Collins FH: Tsessebe, Topi and Tiang: three distinct Tc1-like transposable elements in the malaria vector, Anopheles gambiae.

    Genetica 1999, 105:69-80. PubMed Abstract | Publisher Full Text OpenURL

  48. Kojima KK, Fujiwara H: Evolution of target specificity in R1 clade non-LTR retrotransposons.

    Mol Biol Evol 2003, 20:351-361. PubMed Abstract | Publisher Full Text OpenURL

  49. Oliveira de Carvalho M, Silva JC, Loreto EL: Analyses of P-like transposable element sequences from the genome of Anopheles gambiae.

    Insect Mol Biol 2004, 13:55-63. PubMed Abstract | Publisher Full Text OpenURL

  50. Quesneville H, Nouaud D, Anxolabehere D: P elements and MITE relatives in the whole genome sequence of Anopheles gambiae.

    BMC Genomics 2006, 7:214. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  51. Rohr CJ, Ranson H, Wang X, Besansky NJ: Structure and evolution of mtanga, a retrotransposon actively expressed on the Y chromosome of the African malaria vector Anopheles gambiae.

    Mol Biol Evol 2002, 19:149-162. PubMed Abstract | Publisher Full Text OpenURL

  52. Tubio JM, Tojo M, Bassaganyas L, Escaramis G, Sharakhov IV, Sharakhova MV, Tornador C, Unger MF, Naveira H, Costas J, Besansky NJ: Evolutionary dynamics of the Ty3/gypsy LTR retrotransposons in the genome of Anopheles gambiae.

    PLoS One 2011, 6:e16328. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Finnegan DJ: Eukaryotic transposable elements and genome evolution.

    Trends Genet 1989, 5:103-107. PubMed Abstract OpenURL

  54. Hutchison CAI, Hardies SC, Loeb DD, Shehee WR, Edgell MH (Eds): LINEs and related retroposons: long interspersed repeated sequences in the eukaryotic genome. American Society for Microbiology, Washington, D.C.; 1989. OpenURL

  55. Luan DD, Korman MH, Jakubczak JL, Eickbush TH: Reverse transcription of R2Bm RNA is primed by a nick at the chromosomal target site: a mechanism for non-LTR retrotransposition.

    Cell 1993, 72:595-605. PubMed Abstract | Publisher Full Text OpenURL

  56. Kapitonov VV, Jurka J: SINEX-1_AG, a family of nonautonomous non-LTR retrotransposons from African malaria mosquito.

    Repbase reports 2004, 44-44. OpenURL

  57. Feschotte C, Pritham EJ: Computational analysis and paleogenomics of interspersed repeats in eukaryotes. In Computational genomics: current methods. Edited by Stojanovic . Horizon Bioscience, Norfolk, U.K.; 2007:31-53. OpenURL

  58. Jordan IK, McDonald JF: Comparative genomics and evolutionary dynamics of Saccharomyces cerevisiae Ty elements.

    Genetica 1999, 107:3-13. PubMed Abstract | Publisher Full Text OpenURL

  59. Jordan IK, McDonald JF: Tempo and mode of Ty element evolution in Saccharomyces cerevisiae.

    Genetics 1999, 151:1341-1351. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Petrov DA, Chao YC, Stephenson EC, Hartl DL: Pseudogene evolution in Drosophila suggests a high rate of DNA loss.

    Mol Biol Evol 1998, 15:1562-1567. PubMed Abstract | Publisher Full Text OpenURL

  61. Petrov DA, Lozovskaya ER, Hartl DL: High intrinsic rate of DNA loss in Drosophila.

    Nature 1996, 384:346-349. PubMed Abstract | Publisher Full Text OpenURL

  62. Arkhipova IR, Lyubomirskaya NV, Ilyin YV: Drosophila retrotransposons. Landes, R.G., Austin, Texas; 1995. OpenURL

  63. Slattery JP, Franchini G, Gessain A: Genomic evolution, patterns of global dissemination, and interspecies transmission of human and simian T-cell leukemia/lymphotropic viruses.

    Genome Res 1999, 9:525-540. PubMed Abstract | Publisher Full Text OpenURL

  64. Kaminker JS, Bergman CM, Kronmiller B, Carlson J, Svirskas R, Patel S, Frise E, Wheeler DA, Lewis SE, Rubin GM, et al.: The transposable elements of the Drosophila melanogaster euchromatin: a genomics perspective.

    Genome Biol 2002, 3:RESEARCH0084. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. Abrusan G, Krambeck HJ, Junier T, Giordano J, Warburton PE: Biased distributions and decay of long interspersed nuclear elements in the chicken genome.

    Genetics 2008, 178:573-581. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Eickbush TH: Origin and evolutionary relationships of retroelements. In Evolutionary biology of viruses. Edited by Morse SS. Raven, New York, NY; 1994:121-157. OpenURL

  67. Kajikawa M, Ohshima K, Okada N: Determination of the entire sequence of turtle CR1: the first open reading frame of the turtle CR1 element encodes a protein with a novel zinc finger motif.

    Mol Biol Evol 1997, 14:1206-1217. PubMed Abstract | Publisher Full Text OpenURL

  68. Wicker T, Robertson JS, Schulze SR, Feltus FA, Magrini V, Morrison JA, Mardis ER, Wilson RK, Peterson DG, Paterson AH, Ivarie R: The repetitive landscape of the chicken genome.

    Genome Res 2005, 15:126-136. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Hua-Van A, Le Rouzic A, Boutin TS, Filee J, Capy P: The struggle for life of the genome’s selfish architects.

    Biol Direct 2011, 6:19. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  70. Petrov DA, Hartl DL: Trash DNA is what gets thrown away: high rate of DNA loss in Drosophila.

    Gene 1997, 205:279-289. PubMed Abstract | Publisher Full Text OpenURL

  71. Besansky NJ: A retrotransposable element from the mosquito Anopheles gambiae.

    Mol Cell Biol 1990, 10:863-871. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Felger I, Hunt JA: Detection and characterization of transposable elements.

    Methods Enzymol 1993, 224:322-334. PubMed Abstract OpenURL

  73. Leeton PR, Smyth DR: An abundant LINE-like element amplified in the genome of Lilium speciosum.

    Mol Gen Genet 1993, 237:97-104. PubMed Abstract OpenURL

  74. Gladyshev EA, Arkhipova IR: A subtelomeric non-LTR retrotransposon Hebe in the bdelloid rotifer Adineta vaga is subject to inactivation by deletions but not 5′ truncations.

    Mob DNA 2010, 1:12. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  75. Domingo E, Holland JJ: RNA virus mutations and fitness for survival.

    Annu Rev Microbiol 1997, 51:151-178. PubMed Abstract | Publisher Full Text OpenURL

  76. Casacuberta JM, Vernhettes S, Audeon C, Grandbastien MA: Quasispecies in retrotransposons: a role for sequence variability in Tnt1 evolution.

    Genetica 1997, 100:109-117. PubMed Abstract | Publisher Full Text OpenURL

  77. Granzotto A, Lopes FR, Lerat E, Vieira C, Carareto CM: The evolutionary dynamics of the Helena retrotransposon revealed by sequenced Drosophila genomes.

    BMC Evol Biol 2009, 9:174. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  78. Rebollo R, Lerat E, Kleine LL, Biemont C, Vieira C: Losing helena: the extinction of a drosophila line-like element.

    BMC Genomics 2008, 9:149. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  79. Petrov DA, Hartl DL: High rate of DNA loss in the Drosophila melanogaster and Drosophila virilis species groups.

    Mol Biol Evol 1998, 15:293-302. PubMed Abstract | Publisher Full Text OpenURL

  80. Domingo E, Holland JJ: Mutation rates and rapid evolution of RNA viruses. Raven, New York, NY; 1994. OpenURL

  81. Gabriel A, Willems M, Mules EH, Boeke JD: Replication infidelity during a single cycle of Ty1 retrotransposition.

    Proc Natl Acad Sci U S A 1996, 93:7767-7771. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  82. Plotcon. , ;  .

    [http://emboss.bioinformatics.nl/cgi-bin/emboss/help/plotcon]

    OpenURL

  83. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.

    Mol Biol Evol 2011, 28:2731-2739. PubMed Abstract | Publisher Full Text OpenURL

  84. R: A Language and Environment for Statisitcal Computing. , ;  .

    [http://www.R-project.org]

    OpenURL

  85. Hollander M, Wolfe DA: Nonparametric statistical methods. 2nd edition. John Wiley & Sons, Inc, New York; 1973. OpenURL

  86. Munzel U, Hothorn LA: A unified approach to simultaneous rank tests procedures in the unbalanced one-way layout.

    Biometric Journal 2001, 43:553-569. Publisher Full Text OpenURL