Skip to main content

Advertisement

We’d like to understand how you use our websites in order to improve them. Register your interest.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Abstract

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Background

Non-LTR retrotransposons (nLTR-RTs) are ubiquitous in vertebrate genomes and have profoundly affected the size, structure and function of these genomes [1,2,3,4]. nLTR-RTs constitute a diverse and ancient group of transposable elements whose origin predates the diversification of the main eukaryotic lineages [5]. They can be classified into 28 clades that differ in the number of open-reading frames (ORFs - one or two) and the presence of functional motifs [6]. The mode of mobilization of nLTR-RTs has not been elucidated for all clades but it is likely that, considering their structural similarities, all these elements transpose via a target-primed reverse transcription reaction, as experimentally demonstrated for the R2Bm and L1 elements [7, 8].

Since nLTR-RTs are rarely transmitted horizontally in vertebrates (with the exception of elements of the RTE clade [9,10,11]), the interaction between nLTR-RTs and the genome of their host is among the most intimate and long-lasting co-evolutionary processes found in nature. nLTR-RTs have been a source of evolutionary novelties [3], yet they can also be a source of deleterious alleles [12,13,14]. Thus, vertebrate hosts have evolved a number of repression mechanisms that have in turn affected the evolution of the nLTR-RT [15]. Vertebrate genomes and nLTR-RTs have shaped each other’s evolutionary fate and the signature of these interactions can be seen in the sequence of retrotransposons. For instance, the ORF1 of mammalian L1 carry the signature of adaptive evolution [16,17,18], whereas the recurrent replacement of the promoter region during L1 evolution is indicative of an arms race between L1 and host-encoded repressors of transcription [16, 17, 19, 20].

The base composition of nLTR-RTs may also reflect the nature of the interactions between elements and their hosts. The L1 sequence in mammals is AT-rich with an A bias on the positive strand [16]. From the perspective of L1 transposition, this base composition is not optimal since A-rich sequences are poorly transcribed, due, in part, to the presence of premature poly-adenylation signals causing early termination of transcription [21, 22]. Indeed, a synthetic codon-optimized L1 element replicates much more effectively in a retrotransposition assay [23]. In addition, the A-richness of L1 makes the codon usage of its ORFs poorly adapted for efficient translation [24]. This raises the possibility that the unusual base composition of L1 is a mechanism of self-regulation [25] and may reflect a trade-off between transposition efficiency and the deleterious effect of L1 on its host [21]. For instance, an L1 element with an unbiased base composition will replicate more effectively but this could result in a rate of transposition that is so deleterious for the host that such an element would not persist in the population.

The base composition is one of the most fundamental properties of a DNA sequence because it profoundly affects a number of important functions such as the efficacy of transcription [22, 26], the secondary structure of DNA and RNA molecules, the codon usage [27, 28] as well as the amino acid composition of encoded proteins [29]. All these aspects can potentially affect the reaction of retrotransposition and the overall replicative success of the element [22]. In addition, the base composition of an element can affect the host upon insertion by decreasing the efficacy of transcription of host genes [30], which explains the rarity of AT-rich L1 elements in introns [31], or by modifying the epigenetic status of the region where it inserts [32], for instance by providing novel CpG sites in the case of a GC-rich element. However, since the pioneering work of Lerat et al. [24, 33], the evolution of base composition in transposable elements has not been analyzed in detail, although the number of available genome sequences has drastically increased since these early studies. Here we performed an analysis of the base composition of nLTR-RT in vertebrates. Our goal was to determine if the A-richness of mammalian L1 is a general feature of nLTR-RT or if different clades of nLTR-RT have evolved different nucleotide composition. We also tested if there is a host effect, so that elements belonging to the same clade evolved towards different base composition in different genomes, or if elements from different clades tend to evolve toward similar base composition when in the same genome.

Results

We examined the evolution of base composition in the major clades of nLTR-RT represented in vertebrates (the dataset is available as Supplementary material 1). Our dataset consists of 193 consensus sequences which correspond to the most recently active families (< 5% divergence among sequences within family) in 14 species of mammals (cow, pig, horse, rabbit, human, lemur, armadillo, dog, panda, hyrax, elephant, rat, mouse and opossum), a reptile (the green anole Anolis carolinensis), an amphibian (the African clawed frog Xenopus tropicalis) and five teleost fish (the zebrafish Danio rerio, the three-spined stickleback Gasterosteus aculeatus, the medaka Oryzias latipes, the fugu Takifugu rubripes and the pufferfish Tetraodon nigroviridis). We used the RTclass1 tool [6] to classify the sequences in their respective clade. Our dataset consists of elements belonging to the L1, Tx1, RTE, I, Rex1, CR1, Daphne, Crack, L2 and R4 clades of nLTR-RT. We also included in our study Penelope elements. The number of sequences, and the diversity of clades differed considerably among organisms, from one consensus sequence in each mammalian species, all belonging to the L1 clade, to 56 in the frog, where 6 clades are represented. The genomes analyzed here contain much larger numbers of nLTR-RT families, but we chose to limit our analysis to recently active families to avoid uncertainties when constructing consensus sequences.

Base composition varies among vertebrate lineages and among nLTR-RT clades

We first compared the base composition of the ORF that encodes the reverse transcriptase activity (homologous to the human L1 ORF2) since it is the only region that is shared among all clades of nLTR-RT and because additional ORFs, when present, are not homologous among clades. We first performed a comparison at the level of vertebrate class (mammals, reptiles, amphibians and teleost fish; Fig. 1 and Table 1). There are significant differences among nLTR-RT clades within the same genomes (for instance F2, 33 = 571.02, p < 0.00001 in lizard; F3, 51 = 54.00, p < 0.00001 in fish; using ANOVA) and among organisms for the same clade (for instance F3, 75 = 75.04, p < 0.00001 for L1; using ANOVA).

Fig. 1
figure1

Base composition of ORF2 for the major clades of nLTR-RT in mammals, lizard, frog and teleostean fish. The base composition of the tetrapod CR1 and the fish Daphne elements are presented next to each other since these clades belong to the same Jockey group of elements [6]. For comparison, the base composition of the host’s exomes were obtained for each gene from the RefSeq databases in NCBI

Table 1 Statistics describing the base composition of ORF2 in nLTR-RT

Elements of the L1, CR1, Daphne and Penelope clades show a clear tendency to be enriched in AT when compared with the genes of their host (Fig. 1). In all organisms, the average AT content of the L1 and Penelope clade is significantly higher than in host genes (p < 0.05 for all comparisons using t-test). The difference between element and host genes is significant for Daphne elements in fish (p < 0.0001) but not for CR1 elements in frog (p = 0.06). There are however substantial differences among AT-rich clades within the same genome. For instance in lizard, L1 is on average 67% AT while CR1 is 54% AT. The same clade can also differ substantially in base composition among hosts. For instance, in L1, there is a strong and significant bias in favor of A on the positive strand in lizard and mammals (43 and 42% A versus 24 and 21% T, respectively), a moderate bias in the frog (33% A and 23% T) and an even smaller bias in fish, where A (33%) and T (31%) are represented almost equally. In CR1, there is an A bias in lizard but not in frog. The base composition of Penelope and I elements always showed an A bias on the positive strand (Table 1).

The Rex1, L2 and RTE clades show distinct patterns. The base composition of Rex1 is not significantly different to the base composition of the exome of the source species (p = 0.173 for frog and p = 0.093 for fish; t-test) and does not differ between frog and fish (p = 0.135). The base composition of the RTE clade tends to be GC-rich in fish and lizard, but there are substantial differences among families with proportions of GC ranging from 48 to 61%. Interestingly, RTE families can differ considerably in nucleotide content within the same organism. This is exemplified in medaka, whose genome hosts 3 families of RTE. RTE 2 and 3 have GC content just below 50% (48 and 47%, respectively), while RTE1 contains 60% GC (Supplementary material 1). All fish RTE elements show a G bias on the positive strand (33 to 39% G). The base composition of the L2 clade is equally disparate and shows significant deviation from the composition of host genes. L2 is GC-rich in lizard (55% GC; p < 0.0001 compared with hosts genes using t-test), AT-rich in fish (58% AT; p < 0.0001) and moderately AT-rich in frog (51% AT), C (34%) and then T (30%) being the most represented bases in this organism. It should be noted that the base composition within each vertebrate lineage shows little variation as indicated by small standard deviations on Fig. 1 and Table 1.

We extended our analysis to dinucleotides (Table 2) since genomes exhibit variations in their usage of dinucleotides [34]. We present here the frequencies of dinucleotides in four clades, L1, CR1 (in lizard and frog), Crack (in zebrafish) and L2, since these clades reflect the range of base composition in nLTR-RT. Elements that are AT-rich, such as L1, CR1, Crack and the fish L2, are also enriched in the dinucleotides ApA, ApT, TpA and TpT, while there is a paucity of GC-rich dinucleotides. In the lizard L2, the most abundant dinucleotides are expectedly GC-rich, the four most represented dinucleotides being CpC, CpT, GpG and TpG. The frog L2 is somewhat unusual: the CpT and TpC are abundant, which is consistent with the base composition of the elements, but the next two most common dinucleotides are surprisingly ApG and GpA, although A and G represent respectively only 22 and 14% of nucleotides. Although the frequency of the different dinucleotides generally reflects the base composition of the elements, the observed and expected dinucleotide frequencies differ significantly for all elements and in all organisms (p < 0.0001 for all comparisons using Chi-square test, d.f. = 15). This is explained, in part, by the systematic under-representation of three dinucleotides across species and across clades, CpG, GpT and TpA, suggesting universal selection against these dinucleotides. The most substantial deviation from expectation is found for the frog L2. Here ApC, CpC, CpT, TpC and TpT are substantially under-represented (although C and T are the most abundant nucleotides) but ApG, GpA and GpG are over-represented given the frequency of the constitutive nucleotides.

Table 2 Frequency of nucleotides and ratio between the observed and expected frequencies. For each clade the four most common dinucleotides are highlighted in orange. Ratios that are lower than 0.80 are highlighted in blue

Differences in base composition reflect long-term evolutionary trends

To interpret long-term evolutionary changes in base composition and the possible impact of horizontal gene transfer (HGT), we investigated differences in base composition in a phylogenetic context. To this end, we built phylogenetic trees for the major clades of nLTR-RT (Figs. 2, 3 and 4 and Supplementary material 2) and we estimated the level of similarity among amino acid sequences, within and between clades.

Fig. 2
figure2

Phylogenetic relationships among L1 elements based on the entire ORF2. The base composition of each element is shown above the tree. The tree was built with the maximum likelihood method using the LG + G + I + F model of mutation and its robustness was assessed by 500 bootstrap replicates. The tree was rooted using Tx1 elements. Nodes that are supported by bootstrap values higher than 70% are indicated with an asterisk

Fig. 3
figure3

Phylogenetic relationships among L2 elements based on the reverse transcriptase domains of ORF2. The base composition of each element is shown above the tree. The tree was built with the maximum likelihood method using the LG + G + I + F model of mutation and its robustness was assessed by 500 bootstrap replicates. Nodes that are supported by bootstrap values higher than 70% are indicated with an asterisk. Roman numbers I to VI refer to the L2 sub-clades (see text)

Fig. 4
figure4

Phylogenetic relationships among Tx1 elements based on the entire ORF2. The base composition of each element is shown above the tree. The tree was built with the maximum likelihood method using the LG + G + I + F model of mutation and its robustness was assessed by 500 bootstrap replicates. The tree was rooted using L1 elements as outgroup. Nodes that are supported by bootstrap values higher than 70% are indicated with an asterisk. Elements with comparatively low AT-content are indicated with arrows (see text)

The L1 clade consists of 5 sub-clades: a lizard clade, a mammalian clade, a fish clade (which includes all fish sequences but two) and two frog clades that are not sister to each other (Fig. 2). The topology of the tree suggests the persistence in frog and in fish of ancient and diverse L1 lineages, whose divergence predates the split between teleost fish and tetrapods. The amino acid similarity among the most divergent L1 sequences within the lizard, fish and frog clades are low, with average values of 26.5, 27.5 and 31.2% respectively, suggesting that multiple lineages of L1 can coexist and evolve for extended periods of evolutionary time within the same genome, as previously reported in [16, 35, 36]. Despite this high level of divergence, the base composition remains remarkably constant within vertebrate lineage (Fig. 2), which is consistent with the small standard deviations on Fig. 1. This suggests some long-term selective pressure or functional constraint on L1 to maintain AT-richness (in all vertebrates) and an A bias on the positive strand (in mammals, lizard and frog).

The phylogenetic analysis of L2 elements based on the reverse transcriptase domain shows that this clade consists of six sub-clades (numbered I through VI on Fig. 3): a lizard sub-clade (I), a fish/frog sub-clade (II) and 4 fish sub-clades (III to VI). Elements belonging to the fish/frog sub-clade and to two of the fish sub-clades (sub-clades II to IV) have a single ORF (ORF2) while the other two fish sub-clades (V and VI) have acquired an additional ORF containing an esterase domain. Since these two di-cistronic sub-clades are not sister to each other, it is possible that they have acquired a second ORF independently or that the second ORF was lost in sub-clade IV. The two di-cistronic sub-clades V and VI harbor the strong AT bias observed for other fish elements (~ 61% AT) while elements in the mono-cistronic sub-clades (II and III) contain a large proportion of the nucleotides C and T and elements in clade IV have a T bias, the three other bases being equally represented. Each of the fish sub-clades contains sequences from multiple species and the level of similarity between species in each sub-clade is similar (~ 50% similarity). This suggests that the base composition has been maintained in those genomes since before the species diverged. The lizard L2 sub-clades consist of elements that are all GC-rich. The lizard L2 clade experienced an intense diversification and a number of closely related families are concurrently active in this genome [36]. However, the most divergent subgroups in this sub-clade are on average 36.6% identical at the amino acid level indicating that the GC-richness of L2 in lizard is ancient and has persisted through extended periods of evolutionary time. Most, but not all, lizard L2s have a second ORF, for which we failed to find any conserved domain. It appears that the lizard L2 acquired a second ORF independently twice since the ORF of the most basal elements does not share any homology with the ORF found in the families that experienced a recent diversification.

The Tx1 clade consists of three highly divergent sub-clades (named I through III on Fig. 4). Each of these three sub-clades remains monophyletic in a wider phylogenetic context (data not shown). In all three sub-clades AT-rich elements dominate but four nested sequences, three in sub-clade I and one in sub-clade III (indicated with arrows on Fig. 4), are GC rich (~ 53%). These sequences do not group together and are derived from three different hosts (fugu, stickleback and frog), while all other sequences but one come from zebrafish. This pattern suggests that host-specific forces can cause changes in base composition, or that these four sequences were transferred horizontally from an organism that harbors Tx1 elements with different base composition.

In contrast to the clades described above, the AT-rich base composition of Rex1, I and Penelope has been conserved over extended periods of evolutionary time. For instance, the diversification of Rex1 predates the split between fish and tetrapods and it has persisted in both fish and frog (Supplementary material 2). Yet, despite this ancient history, the base composition of Rex1 has remained constant over a timespan of 525My [37]. It should be noted however that one of the frog sequences is much more similar to its closest fish sequence (78% similarity) than expected, which is suggestive of an ancient event of HGT between frogs and fish, or in both lineages from a common source. Thus, we can’t exclude that the HGT of Rex1 among lineages has contributed to the apparent pattern of homogeneity in base composition among organisms, although the divergences between elements in the rest of the tree are consistent with the host divergence and suggest that vertical transfer is the main mode of transmission of Rex1.

Variation in base composition across the ORFs

In order to determine if the difference in base composition extends outside the ORF containing the reverse transcriptase domain (ORF2), we compared the base composition of ORF2 with the upstream ORFs (ORF1) in those clades that have two ORFs, such as L1, CR1 in tetrapods, Daphne in fish and some L2 elements (Fig. 5). For L1, CR1 and Daphne, ORF1 shows the same nucleotide bias as ORF2 but the bias always tends to be stronger in ORF2 than in ORF1, i.e. ORF2 is always richer in AT than ORF1, and this is true in all vertebrates. In contrast, the two ORFs of L2 differ markedly in nucleotide content. In lizard, ORF2 is GC-rich (~ 55% GC) but ORF1 is AT-rich (54%). In fish, the phylogenetic analysis suggests that ORF1 was independently acquired twice since the two sub-clades containing elements with an ORF1 are not sister to each other. This is supported by the fact that the base composition of their ORFs differs, 55% AT in one of the sub-clades and 45% in the other one, although their ORF2 is equally AT-rich (61% AT). It should be noted that since ORF1 is not homologous among L2 sub-clades, variations in base composition may not result from long term processes acting on base composition but instead may reflect the original nucleotide content of the sequence that was recruited by the element to form a novel ORF.

Fig. 5
figure5

Comparison of the nucleotide content of ORF1 and ORF2 in L1, CR1 (in lizard and frog), Daphne (in fish) and L2

Variation in base composition at the codon level

We then investigated base composition at the codon-level in ORF2 (Fig. 6). The first two positions of codons determine the encoded amino acid but the base composition at the third position, which is mostly neutral, should reflect the mutational process alone. In L1, the base composition shows a strong A bias at the three position of codons, in mammals and lizard, a moderate one in frog and equal contribution of A and T in fish. In tetrapod CR1 and fish Daphne elements, the first and second positions of the codons show similar composition. The two most abundant bases are A and then G at the first position, and A and T at the second position, demonstrating long-term selective constrains with regard to the encoded amino acid. Composition at the mostly neutral third position varies, with similar frequency of AT and GC in lizard CR1, a moderate AT bias in frog CR1 (with T the most abundant base) and a strong AT bias in fish Daphne. In L2, the difference between species is notable. In frog, the three positions are enriched in C and T. In lizard the first and third positions are GC-rich (59% at both positions) while T is the most common base at the second position, followed by C. In fish, L2 tends to be AT-rich at the three positions but there is an A bias at the first position (31% A), a T bias at the second (33%) and third (36%) position. It is interesting to note that the base composition at the third, mostly neutral, position differs among clades between species but also among clades within the same host (ANOVA; F2, 33 = 234.25, p < 0.00001 for lizard; F2, 42 = 14.13, p < 0.00002 for fish). This suggests that a strictly neutral process does not drive the base composition at the third position and that selection is acting on base composition, independently of the protein-coding capability of the sequence.

Fig. 6
figure6

Base composition at the three codon positions in L1, CR1 (in lizard and frog), Daphne (in fish) and L2 ORF2

We then investigated how the base composition at the three positions affects codon usage and the amino acid sequence of the encoded protein. Supplementary figure 3 shows examples of codon usage for L1 and L2 ORF2 in lizard and zebrafish, compared with the codon usage of the hosts’ exomes. For lizard L1, it is always the A-rich codon that is preferred while the opposite trend is found for L2. For instance, codon GAA is used more than 70% of the time in L1 to encode glutamic acid, whereas it is GAG that is used more than 80% of the time in L2, and both codons are used almost equally in the host exome. In contrast, there is very little difference in codon usage between L1 and L2 in fish since both elements are similarly AT-rich and both L1 and L2 show a preference for AT-rich codons compared with the exome. These trends are summarized by the Relative Synonymous Codon Usage (RSCU; Table 3), which shows a preference for codons with an A at the third position in AT-rich elements with an A bias on the positive strand and a preference for an A or a T for the fish AT-rich elements.

Table 3 Relative Synonymous Codon Usage (RSCU) for ORF2 in lizard, frog and fish L1 and L2 elements, lizard and frog CR1 and zebrafish Crack. The highest RSCU for each amino acid is highlighted in grey. Codons with an A at the third position are in bold

The codon usage bias was further investigated using three statistics - Nc, CAI and RCDI. Nc, also called the effective codon usage, can range from 20, when a single codon is used to encode for each amino acid, to 61 when all codons are used equally. We found the lowest values of Nc for elements that have the strongest A bias (mammalian or lizard L1) or AT bias (fish L1, Daphne and Crack) whereas elements that are GC rich (lizard L2) or have a base composition similar to the exome (Rex1) exhibit a higher value (Table 1). This general trend is reflected by the fact that when the GC content at the third position of codon increases, so does the value of Nc (Supplementary material 4).

The Codon Adaptation Index (CAI) is a measure of how closely the synonymous codon usage of a sequence matches that of a reference set, in our case the genome of the host. Table 1 shows the observed and expected average values of CAI given the base composition of the sequence. The values of CAI for the different types of nLTR-RT are remarkably similar to each other and the observed and expected values are almost identical. This indicates that there is no synonymous bias and that the frequency of the different codons fits what is expected given the nucleotide content of the sequence. This is consistent with the weak correlation between the synonymous GC content and the observed values of CAI (r2 = 0.05, p = 0.0016; Supplementary material 4).

The Relative Codon Deoptimization Index (RCDI) is a measure of how different the codon usage in a sequence is relative to a reference set. An RCDI value of 1 indicates that the codon usage of a sequence is identical to the reference and the larger the value of RCDI the larger the difference in codon usage is. The lowest values of RCDI were found for elements with high GC content such as Rex1, the lizard L2, and RTE, and the highest values were found for the high AT rich elements such as L1 and the fish Daphne and Crack (Table 1). In this case, we found a strong negative correlation between the GC content and the RCDI (r2 = 0.57, p < 0.00001; Supplementary material 4), which indicates that the higher the GC content the smaller the difference in the codon usage of the element and the codon usage of the host. Interestingly however, the observed RCDIs are lower than the expected RCDI (z = − 11.44, p < 0.00001), given the base composition of the sequence (Table 1). This suggests that the codon usage of the element is closer to the codon usage of the host than expected given the base composition of its sequence, which indicates a certain level of codon usage adaptation.

Another consequence of nucleotide bias is that it can affect the amino acid composition of the ORFs, which in turn can affect the physico-chemical properties of the proteins as well as their stability. Supplementary figure 5 compares the amino acid composition of L1 and L2 ORF2 in lizard, frog and fish. The A-rich lizard L1 is considerably enriched in amino acids encoded by A-rich codons such as lysine (AAA and AAG) and Isoleucine (ATA, ATT and ATC), which respectively account for 13.6 and 9.8% of ORF2p. In contrast, amino acids encoded by GC rich codons, such as alanine (GCN; 8.6%) and arginine (CGN, AGT, AGC; 8.2%) are more abundant in lizard L2. Similarly, the CT-rich frog L2 encodes a protein enriched in CT-rich encoded amino acids, such as serine (TCN, AGT, AGC; 14.5%) and leucine (CTN, TTA, TTG; 15.9%). As expected, the amino acid composition of L1 and L2 in fish is very similar since these elements have similar nucleotide content.

Impact of base composition differences on transcription

We then examined how biases in base composition can affect the transcription of retrotransposons. Because an AT-rich sequence is more likely to contain premature polyadenylation (polyA) signals, which would result in inefficient transcription [21], we assessed the number of canonical and non-canonical (AATAAA, ATTAAA) polyA signals in ORF2 (Table 1). The number of polyA signals is correlated with the abundance in AT (r2 = 0.55, p < 0.00001; Supplementary material 6). AT rich elements, such as L1 have more polyA signals (up to 23 for ORF2 in lizard) than GC rich elements, which can have zero or one (such as RTE elements). This suggests that the ORFs of some clades may be transcribed much more efficiently than others. It should be noted that the number of polyA signals in L1 seems to exceed the expected number relative to other elements. For instance, the base composition of the lizard L1 and the fish Daphne are almost identical (66.8%AT), yet there are more than two times more polyA signals in lizard L1 (23.0 ± 6.5) than in fish Daphne (10.0 ± 2.9). Similarly, the lizard Penelope has the same base composition as the mammalian L1 (62.7 and 62.3% AT, respectively), yet the mammalian L1 has on average 2.5 times more polyA (13.5 ± 3.5) than the Penelope (5.4 ± 2.6). This suggests either a stronger selection against polyA signals in non-L1 clades or selection in favor of polyA in L1, possibly to tune the level of transcription to a level tolerable by the host.

Mutation pattern

We then investigated if the pattern of mutations in the different elements can account for the difference in base composition. To this end we estimated the relative proportion of mutations in genomic copies relative to the consensus sequence. To make sure that the L1 and L2 datasets were truly comparable we selected subsets of L1 and L2 elements of the same levels of divergence and thus of the same age. We chose to eliminate mutations shared by multiple elements since these mutations are likely inherited from a common progenitor. Thus, we focused our attention on singletons, i.e. mutations that are unique to a sequence. Figure 7 shows the normalized proportion of de novo mutations for L1 and L2 in lizard, frog and fish. Some general trends are apparent, but there are also some differences among elements. In all cases, mutations from C to T and mutations from G to A are the most frequent (except for frog L1 for which we found a large proportion of T to C mutation), and this is true even when mutations in CpG are excluded. Interestingly, the strength of this mutational bias is clade specific. For L1, 42 to 46% of novel mutations are C to T or G to A while this proportion is 33 to 34% for L2. Since the lizard, frog and fish L2 have drastically different base composition, this bias is not related to the nucleotide content of the elements or its CpG content (which is very different among species). It should be noted that mutations at CpG are unlikely to have a substantial effect on the overall base composition of elements since the number of CpG in each sequence represents a tiny fraction of the length (for instance, from 0.9 to 1.8% of all possible dinucleotides in L1 are CpG) and even if all CpG mutated to TpG or CpA this would have a marginal effect on the base composition. In addition, the elements with the highest number of CpG, the lizard L2s (4.1% of all dinucleotides are CpG) have retained a high GC composition over a long period of evolutionary time, although their CpGs are not immune to a high mutation rate toward TpG and CpA (Fig. 7).

Fig. 7
figure7

Normalized percentage of mutations in genomic copies of L1 and L2 in lizard, frog and fish. Only singletons are shown (see text)

Although mutations from C to T are more common than T to C and mutations from G to A are more common than A to G (from 1.2 to 4.8-fold), there are differences among clades and among species. For instance, mutations from T to C and from A to G in lizard L2 and frog L1 account for ~ 30% of all mutations while they account for 14 to 19% for the other elements. Whatever the cause of these differences, it remains that all elements experience a mutation pressure toward an AT-rich nucleotide content, and when all mutation types are combined, we calculated an overall excess of mutations from GC to AT ranging from 1.4 to 2.9.

Discussion

We performed a comprehensive analysis of the base composition of the major clades of nLTR-RTs active in vertebrates. Our results can be summarized as follows. First, we showed that the nucleotide content differs markedly among clades of nLTR-RTs within the same host and that elements belonging to the same clade can differ in base composition among hosts. Using phylogenetic analysis, we demonstrated that nucleotide content remains constant within the same host over an extended period of evolutionary time. It had been shown that the base composition of TEs differs from the base composition of hosts’ genes [38, 39] and that TEs show a tendency to be enriched in AT at the third position of codons, independently of the genome of origin [24, 40]. Based on this difference, it was proposed that the unusual nucleotide content of TEs can be used to identify them in genomes [39]. Here we showed that this picture needs to be revised and that there are some subtle differences in base composition among elements of the same clade (for instance the A-bias of L1 in mammals, lizard and frog versus the similar abundance of A and T in fish) but also some large differences among elements found in the same genomes (for instance L1 and L2 in lizard).

Base composition variation and horizontal gene transfer

A possible explanation for the heterogeneity in base composition within host is horizontal gene transfer (HGT). There are many cases of HGT across kingdoms of class II elements (e.g. DNA transposons) and some of these cases involve vertebrate hosts [41,42,43]. Such events are however exceedingly rare for nLTR-RTs. A recent investigation of HGT based on 759 eukaryotic genomes [11] showed that, although possible, HGT of nLTR-RTs remains largely limited to elements of the RTE clade [9, 44]. The same study identified six potential instances of HGT of Tx1 elements involving aquatic organisms (mostly invertebrates) but these proposed cases are extremely ancient and most cases are supported by a very small number of sequences. These authors also propose a seventh case, namely the HGT of L1 in mammals but their only argument is the apparent lack of L1 elements in monotremes and they are ignoring other arguments in favor of the vertical transmission of L1 in vertebrates, for instance the fact that the phylogeny of L1 in vertebrates matches perfectly the phylogeny of their hosts [16, 45]. Even if nLTR-RTs, other than RTE, were shown to be transferred horizontally in vertebrates, the data currently at hand support the idea that such cases would be extremely rare, and that vertical transmission is the main mode of transfer of nLTR-RT in vertebrates.

Our dataset however contains some indications that HGT has occurred and could explain some of the differences and similarities in base composition. First, we found within the same genome (Medaka) three RTE families with very different base composition, suggestive of horizontal transmissions from different sources. Second, we found a Rex1 element in frog that showed a higher level of similarity with a fish element than this element had with other fish Rex1. Finally, we found in Tx1, four GC-rich sequences nested among AT-rich fish Tx1, which could be explained by HGT. These different cases will need to be analyzed in more details, however, at the time of our study, BLAST searches of public databases using those sequences did not produce any further support for HGT.

For the other clades analyzed here there is no reason to believe HGT played a role, however we can’t fully exclude it. It is indeed possible that the presence of some nLTR-RT lineages in a genome is the result of an ancient event of HGT from a host with a very different base composition. For instance, it is plausible that the GC-rich L2 element of lizard was acquired a long time ago from an unidentified host. It remains however that this element persisted and diversified in the genome of the lizard for an extended period of evolutionary time and yet retained the same base composition, despite a mutational pressure toward a higher AT content. Thus, even if we hypothesize the ancient transfer of nLTR-RT harboring different base composition, the persistence of the nucleotide content over long periods of evolutionary time remains to be explained.

Base composition is likely maintained by selection

There are two main categories of mechanisms that can account for differences in base composition among clades of nLTR-RTs. First, if elements are exposed to different mutational processes, this can lead to different base composition over time. However, an excess of mutations toward AT is observed for all types of elements, even elements that are not AT-rich, and thus mutation alone does not appear to be the main cause of difference in nucleotide content among nLTR-RTs. This is not to say that mutation does not play a role as we did observe some differences in the pattern of mutations of elements among clades and among hosts. For instance, the lowest bias in favor of AT was found for the GC-rich L2 element of lizard and for the CT-rich L1 element of frog, and could thus contribute to the unusual base composition of these elements. The causes of the differences in mutational bias are unclear and will require further investigations that are beyond the scope of this study. Among possible factors, mutations at the hyper-mutable CpG dinucleotides, whose repair affects the probability of mutations at non-CpG sites [46,47,48], editing by APOBEC proteins [49], that have been shown to differentially affect L1 and L2 elements in lizard [50], or GC-biased gene conversion, which will affect differently elements that reside in regions of different recombination [51, 52] will need to be examined.

Alternatively, selective processes may dictate base composition. There are three lines of evidence suggesting that selection plays a role in the evolution of base composition. First, the base composition of a clade remains stable for extended period of evolutionary time within the same host, even though there are no intrinsic reasons why ORF2 should necessarily be AT-rich, as demonstrated by the diversity in nucleotide content of this ORF. This is exemplified in L1, which has retained an AT-rich composition with an A bias in lizard and an AT-rich composition with no A bias in fish, although L1 experienced intense lineage diversification in these organisms. Second, clades of nLTR-RTs have retained their base composition despite a mutational pressure that should have pushed them all toward an AT-rich composition, yet some clades (L2 in lizard and frog) have retained a GC-rich or CT-rich composition. This means that these mutations towards AT are not recruited in the active lineages, possibly because they are not favorable to the replicative success of the elements. Third, when there are more than one ORF the base composition of the two ORFs differs, while a strict mutational model predicts that the two ORFs should harbor similar nucleotide content. This is exemplified in L1, where the AT content in ORF2 is always higher than in ORF1.

What could be the basis of this selection? A possible explanation is that the base composition reflects selection for transcriptional efficiency or inefficiency, depending on the elements. It had been demonstrated that the AT richness of most nLTR-RTs results in poor transcription, possibly as a means of self-regulation of the elements [21, 23]. The fact that in the most AT-rich elements, the three positions of codons are similarly AT-rich is consistent with a process that is independent of the protein-coding capabilities of the ORF. The base composition of Rex1 and L2 in lizard and frog does not fit a self-regulation model since their ORFs are not enriched in AT and potential premature poly-adenylation signals. This leads to two testable hypotheses. Either these elements are transcribed at a much higher rate than AT-rich elements and are repressed (or self-regulated) by other means, or they are intrinsically less efficiently transcribed because of a weak internal promoter. These two hypotheses will require experimental testing.

It is however also likely that selection acts at the translational level. Although the synonymous codon usage does not deviate from expectation (as indicated by the CAI analysis), the overall codon usage of ORF2 is always more similar to the codon of the host than expected given the base composition of the element (as suggested by the observed RCDI values which are always lower than the expected one), which is indicative of a certain level of adaptation to the host. Selection for a more optimal codon usage could also explain why the base composition of ORF1 is less biased than ORF2, since ORF1p needs to be produced in larger amount that ORF2p for successful transposition. It is thus likely that the base composition of nLTR-RTs is evolving in response to the joint effect of selection for lower (or higher) transcription and selection for more (or less) efficient translation of the ORFs.

Conclusions

Our analysis on base composition evolution provides some insights on the nature of the interactions between TEs and their host and among TEs within a genome. The persistence over long time scales of base compositions that are not optimal for the replication of elements support a model of co-existence between nLTR-RTs and their hosts. Interactions between TEs and their host can range from an arms race, where hosts evolve repression mechanisms imposing a selective pressure on TEs to evade repression, to domestication, where TEs and the hosts are peacefully co-existing because they are both indispensable to the survival of each other. There are examples of both models in the literature, however it is unclear which model is the most common in nature. In a recent review, Cosby et al. examined in great details the literature on this topic and they proposed that the arms race model may not be the most prevalent one [53], but that instead, strategies that would allow TEs to persist and multiply without jeopardizing the fitness of the host have been overlooked. In this context the stability of sub-optimal base composition may provide an example of self-regulation of nLTR-RTs to maintain a harmonious relationship with their hosts. Another implication of this research is that it supports the idea that the genome is comparable to an ecosystem in which TEs compete for host resources, i.e. the community ecology of the genome [54, 55]. If we push this metaphor a little further, TEs can possibly occupy different “genomic niches” if they don’t use the same resources and thus co-exist in the genome of their hosts. The co-existence of L1 and L2 in lizard illustrates this scenario. Given their base composition, it is likely that these two elements do not use the same pool of tRNA and of amino acids for their translation and therefore occupy different niches in the genome of their host. The idea that nLTR-RTs can coexist because they differ in their use of resources will need to be better studied both theoretically and experimentally, but present an intriguing research direction to understand the mechanisms that account for the diversity of TEs in genomes.

Methods

The majority of the sequences analyzed here had previously been described [16, 36, 56] or were obtained from Repbase (https://www.girinst.org/repbase/). Classification of sequences was determined using the RTclass1 tool [6]. For all consensi we verified that the ORFs were intact. In the few cases they were not intact, we collected genomic copies and refined the consensus sequences. Note that we only analyzed families with a divergence from consensus lower than 5% to reduce the uncertainty in building consensus sequences. Sequences were aligned at the DNA and protein level using Geneious version 8.1.5 (www.geneious.com), which was used to estimate nucleotide content. Geneious 8.1.5 was also used to estimate the level of similarity among amino acid sequences. For comparison, we estimated the base composition of the genes in each host genome. This was accomplished by downloading the corresponding RefSeq databases from NCBI and by calculating for each gene the nucleotide composition. We first verified that the distribution of base composition values for host genes and for each nLTR-RT clade within each genome was consistent with normality using the Kolmogorov-Smirnov test of normality. Then, base composition between host genes and nLTR_RT clades and among nLTR-RT clades within and between genomes were statistically compared using ANOVA or t-test.

Phylogenetic reconstructions were performed on protein sequences using the maximum likelihood method in MEGA 6.06 [57] or in phyML [58]. The best model of sequence evolution was selected using the AIC criteria with SMS [59]. The robustness of the trees was assessed using 500 bootstrap resamplings. Trees were reconstructed either using the RT domain alone, the entire ORF2 or a concatenation of the endonuclease (from conserved motifs I to IX) and RT domain (from conserved motifs 0 to 8) as defined in [60].

Analyses of the base composition at the codon level were performed using the CAIcal platform at http://genomes.urv.es/CAIcal [61]. This site calculates the base composition at the three positions of codons and estimates statistics used to assess codon bias. For each codon we estimated the Relative Synonymous Codon Usage (RSCU), defined as the number of times a codon is used, divided by the number of synonymous codons encoding the same amino acid [62]. We also calculated three estimators of codon usage bias, Nc, CAI and RCDI. Nc is the effective number of codons and quantifies how much the use of a specific codon in a gene deviates from equal use of all synonymous codons [63]. Its value ranges from 20, when each amino acid is encoded by a single synonymous codon, through 61, when all synonymous codons are equally represented. The Codon Adaptation Index, CAI [64], estimates codon bias given the codon usage of an organism and the nucleotide content of the gene. The CAI ranges from 0 to 1, a value of 1 indicating that it is the most common synonymous codon that is used, which is suggestive of a low codon bias. Significance of CAI is determined by comparing the observed values of CAI with the expected CAI (eCAI), which is an estimator of the random codon usage assuming the base composition of the sequence studied [65]. Finally, we calculated the Relative Codon Deoptimization Index (RCDI), which is a measure of how different the codon usage in a gene is relative to a reference set [66]. The higher the similarity between codon usage of the host and the sequence of interest, the closer the value of RCDI is to 1. Statistical significance of RCDI was assessed by calculating the expected RCDI (eRCDI), which is determined by generating random sequences with similar nucleotide content and amino acid composition to the input sequence [67].

We determined the pattern of mutation of each clade in each species by collecting at least 8 genomic copies for each family. The alignments consisted of the entire second open-reading frame (ORF2). The genomic copies were aligned to each other and to the consensus of their respective family. The different types of mutations were tabulated using the “Find variations/SNPs” in Geneious 8.1.5. Differences among elements may result from de novo mutations or from differences they inherited from their progenitor. Since we were interested specifically in the type of mutations that elements experience after, or at the time of insertion, we excluded mutations that were shared between elements since those were likely inherited from a common progenitor, and may have been filtered by selective processes (due to the constraints acting on the ORF of the elements). Thus, only singletons were compared among sequences. Note that the families analyzed here are very young and, thus, determining the ancestral and derived state of the nucleotide is trivial. For mutations from C to T and from G to A, we distinguished mutations in CpG dinucleotides from mutations in non-CpG context. This is because C and G in a CpG context mutate at a rate 10 to 50 times higher than in non-CpG [68,69,70].

Availability of data and materials

The consensus sequences analyzed here were obtained from RepBase or as supplementary material in [16, 36, 56]. The raw dataset is available as an excel spreadsheet in supplementary material.

References

  1. 1.

    Tollis M, Boissinot S. The evolutionary dynamics of transposable elements in eukaryote genomes. Genome Dyn. 2012;7:68–91.

  2. 2.

    Chalopin D, Naville M, Plard F, Galiana D, Volff JN. Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates. Genome Biol Evol. 2015;7(2):567–80.

  3. 3.

    Warren IA, Naville M, Chalopin D, Levin P, Berger CS, Galiana D, et al. Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates. Chromosom Res. 2015;23(3):505–31.

  4. 4.

    Sotero-Caio CG, Platt RN 2nd, Suh A, Ray DA. Evolution and diversity of transposable elements in vertebrate genomes. Genome Biol Evol. 2017;9(1):161–77.

  5. 5.

    Malik HS, Burke WD, Eickbush TH. The age and evolution of non-LTR retrotransposable elements. Mol Biol Evol. 1999;16(6):793–805.

  6. 6.

    Kapitonov VV, Tempel S, Jurka J. Simple and fast classification of non-LTR retrotransposons based on phylogeny of their RT domain protein sequences. Gene. 2009;448(2):207–13.

  7. 7.

    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(4):595–605.

  8. 8.

    Cost GJ, Feng Q, Jacquier A, Boeke JD. Human L1 element target-primed reverse transcription in vitro. EMBO J. 2002;21(21):5899–910.

  9. 9.

    Kordis D, Gubensek F. Unusual horizontal transfer of a long interspersed nuclear element between distant vertebrate classes. Proc Natl Acad Sci U S A. 1998;95(18):10704–9.

  10. 10.

    Gentles AJ, Wakefield MJ, Kohany O, Gu W, Batzer MA, Pollock DD, et al. Evolutionary dynamics of transposable elements in the short-tailed opossum Monodelphis domestica. Genome Res. 2007;17(7):992–1004.

  11. 11.

    Ivancevic AM, Kortschak RD, Bertozzi T, Adelson DL. Horizontal transfer of BovB and L1 retrotransposons in eukaryotes. Genome Biol. 2018;19(1):85.

  12. 12.

    Boissinot S, Davis J, Entezam A, Petrov D, Furano AV. Fitness cost of LINE-1 (L1) activity in humans. Proc Natl Acad Sci U S A. 2006;103:9590–4.

  13. 13.

    Boissinot S, Entezam A, Furano AV. Selection against deleterious LINE-1-containing loci in the human lineage. Mol Biol Evol. 2001;18(6):926–35.

  14. 14.

    Hancks DC, Kazazian HH Jr. Roles for retrotransposon insertions in human disease. Mob DNA. 2016;7:9.

  15. 15.

    Goodier JL. Restricting retrotransposons: a review. Mob DNA. 2016;7:16.

  16. 16.

    Boissinot S, Sookdeo A. The evolution of LINE-1 in vertebrates. Genome Biol Evol. 2016;8(12):3485–507.

  17. 17.

    Khan H, Smit A, Boissinot S. Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates. Genome Res. 2006;16(1):78–87.

  18. 18.

    Boissinot S, Furano AV. Adaptive evolution in LINE-1 retrotransposons. Mol Biol Evol. 2001;18(12):2186–94.

  19. 19.

    Adey NB, Schichman SA, Graham DK, Peterson SN, Edgell MH, Hutchison CA 3rd. Rodent L1 evolution has been driven by a single dominant lineage that has repeatedly acquired new transcriptional regulatory sequences. Mol Biol Evol. 1994;11(5):778–89.

  20. 20.

    Jacobs FM, Greenberg D, Nguyen N, Haeussler M, Ewing AD, Katzman S, et al. An evolutionary arms race between KRAB zinc-finger genes ZNF91/93 and SVA/L1 retrotransposons. Nature. 2014;516(7530):242–5.

  21. 21.

    Perepelitsa-Belancio V, Deininger PL. RNA truncation by premature polyadenylation attenuates human mobile element activity. Nat Genet. 2003;35:363–6.

  22. 22.

    Han JS, Szak ST, Boeke JD. Transcriptional disruption by the L1 retrotransposon and implications for mammalian transcriptomes. Nature. 2004;429:268–74.

  23. 23.

    Han JS, Boeke JD. A highly active synthetic mammalian retrotransposon. Nature. 2004;429(6989):314–8.

  24. 24.

    Lerat E, Capy P, Biemont C. Codon usage by transposable elements and their host genes in five species. J Mol Evol. 2002;54(5):625–37.

  25. 25.

    Charlesworth B, Langley CH. The evolution of self-regulated transposition of transposable elements. Genetics. 1986;112(2):359–83.

  26. 26.

    Zhou Z, Dang Y, Zhou M, Li L, Yu CH, Fu J, et al. Codon usage is an important determinant of gene expression levels largely through its effects on transcription. Proc Natl Acad Sci U S A. 2016;113(41):E6117–25.

  27. 27.

    Plotkin JB, Kudla G. Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2011;12(1):32–42.

  28. 28.

    Hershberg R, Petrov DA. Selection on codon bias. Annu Rev Genet. 2008;42:287–99.

  29. 29.

    Sueoka N. Correlation between base composition of deoxyribonucleic acid and amino acid composition of protein. Proc Natl Acad Sci U S A. 1961;47(8):1141–9.

  30. 30.

    Han JS, Boeke JD. LINE-1 retrotransposons: modulators of quantity and quality of mammalian gene expression? BioEssays. 2005;27(8):775–84.

  31. 31.

    Medstrand P, van de Lagemaat LN, Mager DL. Retroelement distributions in the human genome: variations associated with age and proximity to genes. Genome Res. 2002;12:1483–95.

  32. 32.

    Slotkin RK, Martienssen R. Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007;8(4):272–85.

  33. 33.

    Lerat E, Capy P, Biemont C. The relative abundance of dinucleotides in transposable elements in five species. Mol Biol Evol. 2002;19(6):964–7.

  34. 34.

    Burge C, Campbell AM, Karlin S. Over- and under-representation of short oligonucleotides in DNA sequences. Proc Natl Acad Sci U S A. 1992;89(4):1358–62.

  35. 35.

    Furano AV, Duvernell D, Boissinot S. L1 (LINE-1) retrotransposon diversity differs dramatically between mammals and fish. Trends Genet. 2004;20(1):9–14.

  36. 36.

    Novick PA, Basta H, Floumanhaft M, McClure MA, Boissinot S. The evolutionary dynamics of autonomous non-LTR retrotransposons in the lizard Anolis carolinensis shows more similarity to fish than mammals. Mol Biol Evol. 2009;26(8):1811–22.

  37. 37.

    Blair JE, Hedges SB. Molecular phylogeny and divergence times of deuterostome animals. Mol Biol Evol. 2005;22(11):2275–84.

  38. 38.

    Shields DC, Sharp PM. Evidence that mutation patterns vary among drosophila transposable elements. J Mol Biol. 1989;207(4):843–6.

  39. 39.

    Andrieu O, Fiston AS, Anxolabehere D, Quesneville H. Detection of transposable elements by their compositional bias. BMC Bioinformatics. 2004;5:94.

  40. 40.

    Jia J, Xue Q. Codon usage biases of transposable elements and host nuclear genes in Arabidopsis thaliana and Oryza sativa. Genomics Proteomics Bioinformatics. 2009;7(4):175–84.

  41. 41.

    Novick P, Smith J, Ray D, Boissinot S. Independent and parallel lateral transfer of DNA transposons in tetrapod genomes. Gene. 2010;449(1–2):85–94.

  42. 42.

    Pace JK 2nd, Gilbert C, Clark MS, Feschotte C. Repeated horizontal transfer of a DNA transposon in mammals and other tetrapods. Proc Natl Acad Sci U S A. 2008;105(44):17023–8.

  43. 43.

    Schaack S, Gilbert C, Feschotte C. Promiscuous DNA: horizontal transfer of transposable elements and why it matters for eukaryotic evolution. Trends Ecol Evol. 2010;25(9):537–46.

  44. 44.

    Walsh AM, Kortschak RD, Gardner MG, Bertozzi T, Adelson DL. Widespread horizontal transfer of retrotransposons. Proc Natl Acad Sci U S A. 2013;110(3):1012–6.

  45. 45.

    Kordis D, Lovsin N, Gubensek F. Phylogenomic analysis of the L1 retrotransposons in Deuterostomia. Syst Biol. 2006;55(6):886–901.

  46. 46.

    Chen J, Miller BF, Furano AV. Repair of naturally occurring mismatches can induce mutations in flanking DNA. eLife. 2014;3:e02001.

  47. 47.

    Walser JC, Furano AV. The mutational spectrum of non-CpG DNA varies with CpG content. Genome Res. 2010;20(7):875–82.

  48. 48.

    Walser JC, Ponger L, Furano AV. CpG dinucleotides and the mutation rate of non-CpG DNA. Genome Res. 2008;18(9):1403–14.

  49. 49.

    Carmi S, Church GM, Levanon EY. Large-scale DNA editing of retrotransposons accelerates mammalian genome evolution. Nat Commun. 2011;2:519.

  50. 50.

    Lindic N, Budic M, Petan T, Knisbacher BA, Levanon EY, Lovsin N. Differential inhibition of LINE1 and LINE2 retrotransposition by vertebrate AID/APOBEC proteins. Retrovirology. 2013;10:156.

  51. 51.

    Duret L, Galtier N. Biased gene conversion and the evolution of mammalian genomic landscapes. Annu Rev Genomics Hum Genet. 2009;10:285–311.

  52. 52.

    Mugal CF, Weber CC, Ellegren H. GC-biased gene conversion links the recombination landscape and demography to genomic base composition: GC-biased gene conversion drives genomic base composition across a wide range of species. BioEssays. 2015;37(12):1317–26.

  53. 53.

    Cosby RL, Chang NC, Feschotte C. Host-transposon interactions: conflict, cooperation, and cooption. Genes Dev. 2019;33(17–18):1098–116.

  54. 54.

    Brookfield JF. The ecology of the genome - mobile DNA elements and their hosts. Nat Rev Genet. 2005;6(2):128–36.

  55. 55.

    Venner S, Feschotte C, Biemont C. Dynamics of transposable elements: towards a community ecology of the genome. Trends Genet. 2009;25(7):317–23.

  56. 56.

    Blass E, Bell M, Boissinot S. Accumulation and rapid decay of non-LTR retrotransposons in the genome of the three-spine stickleback. Genome Biol Evol. 2012;4(5):687–702.

  57. 57.

    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–9.

  58. 58.

    Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.

  59. 59.

    Lefort V, Longueville JE, Gascuel O. SMS: smart model selection in PhyML. Mol Biol Evol. 2017;34(9):2422–4.

  60. 60.

    Metcalfe CJ, Filee J, Germon I, Joss J, Casane D. Evolution of the Australian lungfish (Neoceratodus forsteri) genome: a major role for CR1 and L2 LINE elements. Mol Biol Evol. 2012;29(11):3529–39.

  61. 61.

    Puigbo P, Bravo IG, Garcia-Vallve S. CAIcal: a combined set of tools to assess codon usage adaptation. Biol Direct. 2008;3:38.

  62. 62.

    Sharp PM, Tuohy TM, Mosurski KR. Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res. 1986;14(13):5125–43.

  63. 63.

    Wright F. The 'effective number of codons' used in a gene. Gene. 1990;87(1):23–9.

  64. 64.

    Sharp PM, Li WH. The codon adaptation index--a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987;15(3):1281–95.

  65. 65.

    Puigbo P, Bravo IG, Garcia-Vallve S. E-CAI: a novel server to estimate an expected value of codon adaptation index (eCAI). BMC Bioinformatics. 2008;9:65.

  66. 66.

    Mueller S, Papamichail D, Coleman JR, Skiena S, Wimmer E. Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity. J Virol. 2006;80(19):9687–96.

  67. 67.

    Puigbo P, Aragones L, Garcia-Vallve S. RCDI/eRCDI: a web-server to estimate codon usage deoptimization. BMC Res Notes. 2010;3:87.

  68. 68.

    Duncan BK, Miller JH. Mutagenic deamination of cytosine residues in DNA. Nature. 1980;287(5782):560–1.

  69. 69.

    Sved J, Bird A. The expected equilibrium of the CpG dinucleotide in vertebrate genomes under a mutation model. Proc Natl Acad Sci U S A. 1990;87(12):4692–6.

  70. 70.

    Hwang DG, Green P. Bayesian Markov chain Monte Carlo sequence analysis reveals varying neutral substitution patterns in mammalian evolution. Proc Natl Acad Sci U S A. 2004;101(39):13994–4001.

Download references

Acknowledgements

We thank Imtiyaz Hariyani, Yann Bourgeois, Justin Wilcox, Sebastian Kirchhof and two anonymous reviewers for their helpful comments on the manuscript.

Funding

This research was supported by New York University Abu Dhabi research funds (AD180) to S.B.

Author information

Affiliations

Authors

Contributions

RR and SB generated the datasets, analyzed the data and wrote the manuscript. All authors are taking full responsibility for the work described in this article. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Stéphane Boissinot.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ruggiero, R.P., Boissinot, S. Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons. Mobile DNA 11, 14 (2020). https://doi.org/10.1186/s13100-020-00209-9

Download citation

Keywords

  • Base composition
  • Codon usage
  • Retrotransposons
  • LINE
  • Vertebrate