- Open Access
Tracing the history of LINE and SINE extinction in sigmodontine rodents
Mobile DNAvolume 10, Article number: 22 (2019)
L1 retrotransposons have co-evolved with their mammalian hosts for the entire history of mammals and currently compose ~ 20% of a mammalian genome. B1 retrotransposons are dependent on L1 for retrotransposition and span the evolutionary history of rodents since their radiation. L1s were found to have lost their activity in a group of South American rodents, the Sigmodontinae, and B1 inactivation preceded the extinction of L1 in the same group. Consequently, a basal group of sigmodontines have active L1s but inactive B1s and a derived clade have both inactive L1s and B1s. It has been suggested that B1s became extinct during a long period of L1 quiescence and that L1s subsequently reemerged in the basal group.
Here we investigate the evolutionary histories of L1 and B1 in the sigmodontine rodents and show that L1 activity continued until after the L1-extinct clade and the basal group diverged. After the split, L1 had a small burst of activity in the former group, followed by extinction. In the basal group, activity was initially low but was followed by a dramatic increase in L1 activity. We found the last wave of B1 retrotransposition was large and probably preceded the split between the two rodent clades.
Given that L1s had been steadily retrotransposing during the time corresponding to B1 extinction and that the burst of B1 activity preceding B1 extinction was large, we conclude that B1 extinction was not a result of L1 quiescence. Rather, the burst of B1 activity may have contributed to L1 extinction both by competition with L1 and by putting strong selective pressure on the host to control retrotransposition.
LINEs (Long INterspersed Elements) are autonomous non-LTR retrotransposons. They move through an RNA intermediate, but have non-homologous ends and use target-primed reverse transcription . L1 (LINE-1) is the most successful family of LINEs in eutherian mammals  and comprise ~ 20% of a mammalian genome [3,4,5,6,7]. A functional full-length L1 is typically 6000–7000 bp long and composed of a 5′ untranslated region (5’UTR) harboring an RNA polymerase II promoter, two non-overlapping open reading frames (ORFs) known as ORF1 and ORF2 and a 3’UTR followed by a poly-adenosine sequence . The structure of L1 can be diverse among different mammals, particularly in the 5’UTR and ORF1 . The ORF-encoded proteins are strictly required for L1 retrotransposition and are highly cis-preferential [9, 10]. L1s are adenosine rich (~ 40%) on their coding strand, which results in biased codon usage compared to host genes [11, 12], elongation defects , and premature RNA splicing . This A-richness contributes to the inefficiency of L1 retrotransposition and is proposed to regulate the genes in their vicinity .
SINEs (Short INterspersed Elements) are relatively short, non-autonomous, non-LTR retrotransposons. SINEs do not encode proteins for their own retrotransposition and depend on the reverse transcriptase encoded by other transposable elements such as LINEs [15, 16]. Although L1s are highly cis-preferential [9, 10], SINEs can take advantage of L1-encoded proteins for their own retrotransposition [15,16,17], and L1 ORF2 protein is sufficient to drive B1 retrotransposition. Despite their short length, SINEs account for ~ 10% of a typical mammalian genome due to their high copy numbers [3, 4]. Among the ~ 70 SINE families found in mammals , B1 is the most abundant in mouse  and possibly most rodent species . B1s are derived from the RNA component of signal recognition particle 7SL RNA [20, 21] and share features with its ancestors. A functional B1 is ~ 150 bp long and transcribed by RNA polymerase III with the aid of its two transcription factor binding boxes [22, 23]. B1 sequences are rich in CpG sites, which are methylated and thus prone to mutation in mammalian genomes , and the elevated mutation rate is pronounced compared to the A-rich L1s [25, 26].
Both L1 and B1 have long histories of co-evolution with their host genomes. Unlike some LTR retrotransposons [27, 28], there is no known targeted mechanism for L1 excision and thus L1s persist in the genome unless they are removed by non-specific mechanisms. L1s can be found in all placental mammals and marsupials [2, 5, 29]. Mammalian L1s evolve as master lineages so that a single or a few lineages are responsible for the total retrotransposition in a short time window [30,31,32,33]. New master elements replace the old ones, eventually dominating retrotransposition, and this replacement process happens recurrently. B1s are younger than L1s, having arisen just before the divergence of the common ancestor of rodents, ~ 65 MYA , and they are specific to rodents. Other SINEs, including B2, B4 and ID elements, are also present in rodent genomes . SINE families have been interacting with L1s for more than 100 MYA, and fossil remnants of extinct SINE families are detectable in well-characterized mammalian genomes [18, 35]. Despite being targeted by a slew of host restriction mechanisms, L1 and B1 compose approximately a quarter of a typical rodent genome [4, 7]. For example, in the mouse genome, there are ~ 599,000 total copies of L1, responsible for ~ 19% of the genome , of which ~ 3000 copies are potentially functional , and ~ 564,000 copies of B1, responsible for ~ 3% of the genome .
LINEs and SINEs have considerable impact on the mammalian genome, although they were traditionally viewed as “junk DNA”. As LINEs and SINEs, including L1s and B1s, retrotranspose and recombine, they introduce genome instability  and cause disease . These elements may occasionally be co-opted by the host to serve host certain functions, such as their proposed roles in neuro-plasticity [39, 40], X chromosome inactivation [41, 42], regulatory functions [43, 44], and DNA break-repair . Our current picture of the effects of retrotransposition does not take into account what effect silencing of LINEs and SINEs might have on these proposed functions.
Since L1 retrotransposition is under strict control by multiple host defenses , it might seem reasonable for the host to occasionally win the evolutionary arms race with L1s, resulting in loss of L1 activity (L1 extinction). Yet, L1 extinctions are relatively rare. Factors contributing to the rarity of L1 extinction could be but not limited to co-option of L1s for certain essential function to the host. Mammalian L1s are not known to move horizontally, although ancient L1s might have been able to do so given evidence of other non-LTR retrotransposons . The unlikely replacement of L1s by horizontal movement would be evident as a mismatch between L1 phylogeny and host phylogeny. Therefore, mammalian L1 extinctions would affect all derived host species. Two factors are of note here. First, clades with early L1 extinctions could have given rise to large mammalian lineages without L1 activity and be easily detected because of both the number of species affected and the deterioration of the remnant sequences in the genome. Secondly, recent extinctions will be difficult to differentiate from periods of L1 quiescence. To clarify the terms related to loss of L1 activity in this work, we refer to a period of low L1 activity as “quiescence” and complete loss of L1 activity as “extinction”. Given the large phylogenetic impact of early extinctions, one might expect L1s to eventually become extinct in most mammalian genomes, and yet L1s have persisted throughout the entire evolutionary history of their placental mammal and marsupial hosts. Thus, either most L1 extinctions are either recent or rare, or mammalian lineages subject to ancient L1 extinctions do not persist or they give rise to few new species. Understanding the dynamics of L1 extinction will be as important as understanding the dynamics of L1 activity in sorting out the impact of L1s on mammalian genome evolution.
Several cases of L1 extinction have been proposed in the literature [48,49,50,51,52,53,54,55,56] and two of these are deep extinction events that cover major groups of mammals [50,51,52]. One of the major L1 extinctions occurred in a large group of South American rodents and includes most species in Sigmodontinae. Sigmodontinae is a subfamily of the Cricetidae family, including approximately 377 species classified into 74 genera in nine tribes (Fig. 1)  and thus contains 7–8% of the estimated 5000 mammalian species . Given that B1 retrotransposition is dependent on that of L1, it is expected that B1s should lose their activity simultaneously with L1s. However, the B1 extinction in Sigmodontinae appears to have preceded that of L1s based on samples from 14 genera in five tribes [50,51,52], where the basal genus Sigmodon carries inactive B1 and active L1, and the descendant genera carry both inactive L1 and B1 (Fig. 1). It has also been shown that loss of L1 and B1 activity follows the expansion of a group of endogenous retroviruses [60, 61].
It was previously hypothesized that L1 can experience long-term quiescence as a “stealth driver” , and that B1 extinction could happen during this period of L1 quiescence . Since B1s are more prone to mutations than the average sequence due to enriched CpG content , Rinehart et al.  hypothesized that B1 was unable to retrotranspose at a high enough rate during L1 quiescence to replace their active copies, accumulating debilitating mutations more rapidly than L1s. When a more active family of L1 emerged in the Sigmodontini, it was hypothesized that B1 was too degraded to retrotranspose, resulting in B1 extinction even in the presence of high L1 activity.
In this study, we investigate the evolution histories of L1 and B1 spanning the time of their extinctions and the radiation of the extant species in Sigmodontinae (Fig. 1). Since the group carrying extinct L1s and B1s (Oryzomyalia, Fig. 1) shares a common ancestor, we used the marsh rice rat Oryzomys palustris to represent this group, hereafter referred to as the “L1-extinct clade”. We used the hispid cotton rat Sigmodon hispidus to represent the clade carrying active L1 but inactive B1, hereafter referred to as the “basal group”. We used the deer mouse Peromyscus maniculatus to represent a closely related clade carrying both active L1 and B1, hereafter referred to as the “outgroup”.
Using unassembled genome sequences from the species representing the L1-extinct clade and the basal group, we show that the activity of L1 and B1 families preceeding the divergence of the clades is comparable in the current genomes of the two groups. L1 families had been steadily replaced before the split of the two groups and maintained activity after the split of the basal group and the L1-extinct clade. Shortly after this split L1 activity ceased in the L1-extinct clade but became highly active in the basal group. B1s, on the other hand, had a very large increase in activity prior to the split between the L1-extinct clade and the basal group, and there is no evidence of activity in the two groups following their divergence.
To investigate the history of L1 retrotransposition in O. palustris and S. hispidus, we used COSEG  to identify closely related L1 groups based on shared, co-segregating sites as described in Methods. We follow the convention of COSEG to designate these groups as subfamilies. RepeatMasker  was used to assign genomic L1 copies to COSEG-generated subfamilies, and seven subfamilies with no assigned sequences were removed from further consideration, leaving 47 subfamilies for further analysis.
To examine the activity of L1s in O. palustris and S. hispidus, we searched the trace files of both genomes separately with the consensus sequences of the abovementioned 47 subfamilies and identified 19,254 sequences in O. palustris and 90,526 in S. hispidus. The relative age of each sequence was approximated by its percent divergence from the corresponding subfamily consensus - the higher the percent divergence, the older the sequence. The peak of the distribution was used as an approximation of the divergence of the subfamily (Additional file 3 Table S1). Given the possible changes of evolution rate in the detectable range of L1 evolutionary history, a global conversion from percent divergence to time is challenging. However, because of the shared evolutionary history of O. palustris and S. hispidus, we are able to use percent divergence as a reasonably good marker to compare the relative ages of L1 subfamilies of the two species, as is typical in studies of transposable elements .
Subfamily consensus sequences were also used to infer phylogenetic relationships between subfamilies (Additional file 1: Figure S1). Subsequently, phylogenetic relationships and sequence similarities between subfamilies were used to assign subfamilies to families with the stipulation that the pairwise distance between subfamilies within a family be no greater than 3.5%. This distance was determined operationally based on the divergences among phylogenetically clustered subfamilies. Clusters of subfamilies that were similar at the sequence level but differed in divergence were assigned to different families. This process identified five families specific to S. hispidus (L1-S1 to L1-S5), four families shared by O. palustris and S. hispidus (L1-OS1 to L1-OS4) and two shared by P. maniculatus, O. palustris and S. hispidus (L1-OSP1 and L1-OSP2, Additional file 3: Table S1). A distance-based phylogeny reflecting the relationship between L1 families is presented in Fig. 2a. Individual sequences were assigned to the families to which their subfamilies belong; the divergence within a family is based on the distance of each sequence from its subfamily consensus (Fig. 3).
As expected, sequences from L1 families shared by O. palustris and S. hispidus are present in both genomes, and these shared families are fairly synchronized in time and comparable in copy number (Fig. 3a). L1-OS1 is the only shared L1 family between O. palustris and S. hispidus that shows a difference: it is the youngest shared L1 family, the last active L1 family prior to the L1 extinction, and has ~ 1.5-fold higher copy numbers per Gbp of sequence in O. palustris than in S. hispidus. This difference in L1-OS1 deposition between O. palustris and S. hispidus suggests that L1s remained active in the L1-extinct clade after the separation of that group from the basal group. The Sigmodon-specific L1 families (Fig. 3b, families S1–5) experienced substantial amplification after divergence from the L1-extinct clade, whereas no Oryzomys-specific subfamilies were identified by COSEG. The Sigmodon-specific subfamilies had a few sequences from the O. palustris genome assigned to them, but these assignments appear to be anomalous since the sequences are highly divergent from the subfamily consensus sequences (Additional file 3: Table S1). Thus, L1 experienced an expansion (L1-OS1) in the lineage leading to Oryzomyalia immediately before L1 extinction, while the lineage leading to Sigmodontini experienced a delayed but much larger L1 expansion.
In order to study the B1 dynamics in sigmodontine rodents, we performed the analysis on B1 similarly to that done on L1. Because of the short length and CpG-rich nature of B1, we required twice as many sequences to form a subfamily in the second round of COSEG as described in Methods. The analysis revealed 30 subfamilies and five families of B1 in both species (Additional file 4: Table S2). A distance-based phylogeny reflecting the relationships between B1 families is presented in Fig. 2b. One of the families (B1-OS1) is shared by O. palustris and S. hispidus and the subfamilies within B1-OS1 form a polytomy on the distance-based phylogeny (data not shown). The other four B1 families (B1-OSP1–4) are shared by O. palustris, S. hispidus and P. maniculatus. The representation of these families in both O. palustris and S. hispidus genomes is fairly synchronized in time and comparable in copy number (Fig. 4). Since the outgroup, represented by P. maniculatus, carries both active L1s and B1s, we know that B1 extinction happened after the split of the outgroup, yet the point at which B1 lost activity in the basal group is to be determined. Here we show that the peak of the most recent B1 family resides at ~ 11.1% divergence in O. palustris and ~ 10.7% in S. hispidus (Additional file 4: Table S2). These peaks reside in the same divergence window as L1-OS2 (~ 11.1% in O. palustris and ~ 10.3% in S. hispidus, Additional file 3: Table S1), suggesting that B1-OS1 is coincident in time with L1-OS2. Exclusion of CpG sites when calculating percent divergence reduces the variation between the L1 and B1 clocks, but we acknowledge that other minor mutation rate variations between L1s and B1s might still exist. Since L1-OS2 is the youngest L1 family prior to the separation of the basal group and the L1-extinct clade, the last wave of B1 retrotransposition likely preceded the extinction of L1.
In this paper we explore the tempo of L1 and B1 activity surrounding the extinction of both elements that occurred in most species within the rodent subfamily Sigmodontinae. This work is made possible by sequencing methods that allow us to gather large amounts of sequence data and by the availability of a robust species phylogeny for the group (Fig. 1). A recent phylogenetic analysis of muroid rodents  indicates that the tribe Sigmodontini is basal in the Sigmodontinae and sister to the tribe Ichthyomyini. These two tribes are sister to a large, polytomic group (the Oryzomyalia) which includes the remaining six tribes. The subfamily is the result of a rapid radiation of rodents into South America about 5 MYA . Previous work indicated that L1s are extinct in the Oryzomyalia but active in the Sigmodontini, which is composed of 14 species in one genus, Sigmodon. A summary of total L1s and B1s in S. hispidus and O. palustris (Additional file 2: Figure S2) agrees with this pattern of extinction. L1 extinction in the Oryzomyalia has been documented in 13 genera distributed across four tribes spanning this group (Fig. 1). Evidence for this L1 extinction included sequence divergence between L1s cloned by a method that enriches for recently transposed elements , and faint hybridization along with the absence of species- or genus-specific bands in a Southern blot of the 13 genera when probed with L1 .
We reconstructed the shared evolutionary history of L1s and B1s in Sigmodontinae in the period preceding and following extinction of these elements. Our results suggest that L1 master elements were replaced steadily prior to the extinction of both L1 and B1. This is reflected by the consecutive series of L1 families shared by O. palustris and S. hispidus after their divergence from Peromyscus. B1 elements did not appear to take advantage of every wave of L1 activity, but a wave of L1 retrotransposition (family L1-OS2, red color in Figs. 3 and 5) corresponds to the B1 retrotransposition peak just prior to B1 extinction (B1-OS1, red color in Figs. 4 and 5).
By determining the minimum distance of individual L1s and B1s from their consensuses in species spanning 12 genera of the Sigmodontinae, Rinehart et al.  proposed that L1 extinction occurred after the split between the L1-extinct clade and Sigmodon, the basal genus. We further refine this timeframe in a summary diagram showing the higher level of L1-OS1 activity in O. palustris compared to S. hispidus (Fig. 5, magenta color). This supports the conclusion that there was some level of L1 activity in both species after the split, and that the events leading to L1 extinction also happened after the split, rather than via recovery of L1 activity in S. hispidus as previously suggested . The evolutionary history of B1 in O. palustris and S. hispidus is comparable. New B1 deposition into the genome was unremarkable except for the large burst in both species in the period directly preceding B1 extinction (B1-OS1, red color in Figs. 4 and 5). Given the short length of B1s, it is more difficult to identify subfamily clusters, so our estimation of the timing of B1 extinction is weaker than for L1. However, two lines of evidence suggest that the last burst of B1 activity occurred prior to the split between the L1-extinct and basal groups. Firstly, the peak activity of B1-OS1 corresponds most closely to the peak activity of L1-OS2, which appears to precede the split of these two rodent clades (Fig. 5, red color). Secondly, there is no indication of difference in activity for any of the B1 subfamilies in O. palustris and S. hispidus (Additional file 4: Table S2), as was the case for L1 (Additional file 3: Table S1).
Estimation of retrotransposition rate based on historical L1 copy numbers could be affected by the excision rate of the host genome or detection limit of the algorithm. Although no known mechanism specifically targets L1 and B1 for excision, mammalian genomes have been constantly expelling sequences by various mechanisms and the excision rate varies in different groups . Older insertions are exposed to the non-targeted excision mechanisms for longer time, thus fewer copies of the older families are expected to be represented in the genome. Old L1 and B1 copies also suffer from more limited recognition by the available algorithms. The sequences detectable by RepeatMasker decrease drastically beyond 30% divergence. Since the mutation rate in the rodent lineage is one of the highest in all mammals, 30% divergence in L1 and B1 only traces back to the common ancestor of sigmodontine rodents and P. maniculatus, while similar studies on bats  and primates [67, 68] trace back to the common ancestor of mammals. Fortunately, P. maniculatus carries both active L1s and B1s and is close enough to serve as an outgroup in this study. We were able to identify two L1 families shared by O. palustris, S. hispidus and P. maniculatus, L1-OSP1 and L1-OSP2. However, there is an advantage of studying rodents in this type of evolutionary study. Since the mutation rate in the rodent lineage is higher than that of other mammals, evolution of L1 and B1 subfamilies over a given period of time will show greater divergence compared to more slowly evolving species. This gives the divergence distributions of L1s and B1s higher resolution and allows us to discern subtle differences between subfamily divergence.
The patterns of historical L1 and B1 activity reveal the critical time frame at which retrotransposition rates diverged in the ancestral hosts. It is apparent that L1s were still active at the time when the ancestral lineages of Oryzomyalia and Sigmodontinae split. Thus, rather than quiescence of L1 in the common ancestor with resurgence in the sigmodontine lineage, the extinction of L1 appears to have occurred after the split and only in the Oryzomyine lineage. L1s are extinct in Oryzomyalia but active in Sigmodontini. B1s are extinct in Oryzomyalia and Sigmodontini. However, the status of both L1s and B1s in the intermediate tribe, Ichthyomyini, is unknown. Thus, L1 extinction from this single event likely affects between 345 and 362 species, or about 7% of all mammalian species.
Our study also reveals the largest and final wave of B1 expansion indeed occurred in the common ancestor of the host lineages. This huge burst of B1 activity suggests an explanation other than L1 quiescence, and the subsequent deficit of L1 proteins, for its extinction. Given that transposable elements can be involved in an evolutionary arms race with host restriction factors , it is possible that the radical expansion of B1 triggered stronger host defense and eventually led to the extinction of L1 and B1. Therefore, reconstruction of the evolutionary history of L1 and B1 host restriction factors in relevant rodent species could be the key to revealing the mechanism of L1 and B1 extinction.
O. palustris DNA was obtained from the Natural Sciences Research Laboratory in The Museum of Texas Tech University (tissue ID: TK28621), and S. hispidus DNA was obtained from Texas Cooperative Wildlife Collection of the Texas A&M University (tissue ID: MUR15). Whole genome sequencing of each species was done in separate batches using MiSeq (Illumina, Inc., San Diego, CA) at the IBEST Genomic Resources Core (University of Idaho, Moscow, ID). Paired-end libraries were generated with an insert size of 450–550 bp; ~ 13 and 14 million total reads were generated for O. palustris and S. hispidus, respectively. Sequences were processed with SeqyClean (https://bitbucket.org/izhbannikov/seqyclean) and the paired-ends were joined with FLASH . Genome coverage was equivalent to approximately 1.5X; 5.47 Gbp of sequence were generated for O. palustris and 6.06 Gbp for S. hispidus, but we note that genome size within the sigmodontine rodents varies. Although the genome size of O. palustris is not documented to our knowledge, the genome size of sister species in Oryzomys suggest that Sigmodon genomes are 11–16% larger than those of Oryzomys .
L1 reconstructions for both species were generated based on partial genomic sequences generated by 454 Pyrosequencing (Roche Applied Science, Penzberg, Germany) at the IBEST Genomic Resources Core, 203 Mbp of sequence for O. palustris and 214 Mbp for S. hispidus. P. maniculatus genome trace files were obtained from NCBI FTP site (ftp://ftp.ncbi.nlm.nih.gov/pub/TraceDB/peromyscus_maniculatus/). Reconstruction of the 3′ ends of O. palustris and S. hispidus L1s started with a 575 bp consensus seed in the 3′ half of L1 ORF2 generated following Cantrell et al. . A bioinformatic pipeline for reconstructing a full length L1 is described by Yang et al. . Briefly, sequences were acquired from the genome trace files based on percent identity [72, 73]. The overhangs of the found sequences allowed the creation of new seeds at both ends of the L1 fragment and were used to initiate another round of query. In this case, the reconstruction walk was repeated in the 3′ direction until the 3′ end of ORF2 was reached. Percent identity cutoff was set at 92% for O. palustris to capture the most recently active sequences; a higher percent identity (97–99%) and overhang length of at least 100 bp were used for S. hispidus. This assured a satisfactory consensus and exclusion of older L1 elements for each species. The 3′ 300 bp of the reconstructed L1s were then used as the reference sequences for COSEG analysis described below.
B1 sequences from Rinehart et al.  were used as starting seeds for B1 analysis. The PCR-amplified B1s from O. palustris and S. hispidus were aligned with Lasergene MegAlign (DNASTAR, Madison, WI) and the consensus sequence (146 bp) was used as the reference sequence for COSEG analysis.
L1 and B1 subfamilies in O. palustris and S. hispidus were identified and characterized in similar fashion as described below and are summarized in Additional file 3: Table S1 and Additional file 4: Table S2. The alignment of L1 and B1 subfamily sequences are available in Additional file 5 and Additional file 6.
The reconstructed 300 bp sequences from the 3′ end of O. palustris and S. hispidus L1 ORF2 were each used as the initial L1 query sequences, and the full length B1 consensuses from each species, based on Rinehart et al. , were used as the initial B1 query sequences. O. palustris and S. hispidus MiSeq genomic DNA libraries (processed reads without assembly) were queried using RepeatMasker  with default parameters. This initial step was conducted to identify any sequence that is potentially homologous to L1/B1, which saves computation time of the following steps by avoiding the need to query all sequencing reads directly. Hits from RepeatMasker searches were filtered for > 90% coverage of the query sequence and subsequently used for the first COSEG  (http://www.repeatmasker.org/COSEGDownload.html) run to identify subfamilies based on shared, co-segregating sequence variants. All COSEG runs were conducted under default parameters except as noted. Parameters were set such that at least 250 sequences were required to form a L1 subfamily and 1000 were required to form a B1 subfamily. In order to identify older subfamilies, the consensus sequences of the subfamilies identified by the first COSEG run were used as queries to again search the O. palustris and S. hispidus MiSeq libraries using RepeatMasker. The identified sequences from the second RepeatMasker run were filtered for > 90% coverage and extracted. O. palustris and S. hispidus sequences are combined and a second COSEG run was carried out on the combined sequences. To avoid the possible formation of random subfamilies due to the short length of B1 and the high copy number of the detected sequences, the sequences required to form a subfamily was increased from 1000 (for the former separate run) to 2000, whereas this number for L1 remained unchanged at 250. The consensus sequences of the resulting COSEG subfamilies were trimmed to exclude ends that were not common to all subfamilies and the CpG sites were removed and, thus, treated as gaps by RepeatMasker and not counted for the divergence calculation. These modified subfamily consensus sequences were used for a final query of the individual O. palustris and S. hispidus MiSeq libraries using RepeatMasker. Sequences from this third run were assigned to subfamilies based on percent divergence and this information was stored for further analysis.
P. maniculatus genome trace files were data-mined in a similar fashion through a single round of RepeatMasker and COSEG. The O. palustris L1 and B1 sequences described above were used as the initial query seeds for this run. Selected P. maniculatus subfamilies were used to demarcate the divergence of the subfamilies identified in the O. palustris and S. hispidus genomes (Fig. 3).
Subfamily consensus sequences generated by the second COSEG run of the O. palustris and S. hispidus libraries were combined and aligned with MegAlign using the Clustal W method for L1 or Clustal V method for B1 and a distance matrix was calculated based on the alignment. Due to the use of subfamily consensuses, gaps are rare in the alignment and hence kept. CpG sites were manually removed from the alignments. Based on the alignments, maximum likelihood trees were constructed using PhyML  with the GTR + I + G model and 100 bootstrap replicates (Additional file 1: Figure S1). L1 and B1 sequences were then assigned to families based on the topology of the tree and a no more than 3.5% within-family pairwise distance from their subfamily consensuses for L1 and 4.4% for B1. Given that the L1 and B1 masters are constantly being replaced during evolution, perfect designation of large families is not possible. The 3.5% threshold was chosen so as to cluster closely related subfamilies without inflating the number of families. Families are named according to their species-specificity and divergence: “S” indicates Sigmodon-specific families, “OS” for families shared by Sigmodon and Oryzomys and “OSP” for families shared by Sigmodon, Oryzomys and Peromyscus; numbers in family names indicates the divergence of a family within the family group with “1” being the youngest. Family consensus sequences were generated in MegAlign (DNAStar, Madison, WI) using the consensus sequences of subfamilies belonging to each family. Alignment and phylogeny of families were generated as described above for subfamilies. Histograms of L1 and B1 divergence distributions were generated by R  histogram function using a bin size of 1% (Figs. 3 and 4). Percent divergence corresponding to retrotransposition peaks of individual families and subfamilies were determined by R using the kernel smoothing function with 0.4% bandwidth (Additional file 3: Table S1 and Additional file 4: Table S2).
To avoid any bias introduced by using the consensus-based seeds, we performed a similar analysis on the L1s in S. hispidus using RepeatScout , which is a de novo repeat identification method that does not use any priori of known repeats. RepeatScout was run with default parameters on the processed MiSeq library of S. hispidus to find repetitive sequences. Identified repeats were then annotated using RepeatMasker. All L1-like sequences that overlaps with the 3′ 300 bp of S. hispidus L1 ORF2 were used as seeds to perform a COSEG analysis of L1s in S. hispidus following the approach described above. All RepeatScout-based L1 subfamilies (Additional file 7) in S. hispidus were within 3.5% divergence of a COSEG-defined subfamily as described above, and hence can be assigned to the same L1 family.
Long INterspersed Element
Million Years Ago
- O. palustris :
Open Reading Frame
- P. maniculatus :
- S. hispidus :
Short INterspersed Element
Kapitonov VV, Jurka J. A universal classification of eukaryotic transposable elements implemented in Repbase. Nat Rev Genet. 2008;9(5):411–2.
Smit AF. The origin of interspersed repeats in the human genome. Curr Opin Genet Dev. 1996;6(6):743–8.
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409(6822):860–921.
Mouse Genome Sequencing C, Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, et al. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420(6915):520–62.
Boissinot S, Sookdeo A. The evolution of LINE-1 in vertebrates. Genome Biol Evol. 2016;8(12):3485–507.
de Koning AP, Gu W, Castoe TA, Batzer MA, Pollock DD. Repetitive elements may comprise over two-thirds of the human genome. PLoS Genet. 2011;7(12):e1002384.
Platt RN 2nd, Vandewege MW, Ray DA. Mammalian transposable elements and their impacts on genome evolution. Chromosom Res. 2018;26(1–2):25–43.
Furano AV. The biological properties and evolutionary dynamics of mammalian LINE-1 retrotransposons. Prog Nucleic Acid Res Mol Biol. 2000;64:255–94.
Wei W, Gilbert N, Ooi SL, Lawler JF, Ostertag EM, Kazazian HH, et al. Human L1 retrotransposition: cis preference versus trans complementation. Mol Cell Biol. 2001;21(4):1429–39.
Kulpa DA, Moran JV. Cis-preferential LINE-1 reverse transcriptase activity in ribonucleoprotein particles. Nat Struct Mol Biol. 2006;13(7):655–60.
Han JS, Boeke JD. A highly active synthetic mammalian retrotransposon. Nature. 2004;429(6989):314–8.
An W, Dai L, Niewiadomska AM, Yetil A, O'Donnell KA, Han JS, et al. Characterization of a synthetic human LINE-1 retrotransposon ORFeus-Hs. Mob DNA. 2011;2(1):2.
Han JS, Szak ST, Boeke JD. Transcriptional disruption by the L1 retrotransposon and implications for mammalian transcriptomes. Nature. 2004;429(6989):268–74.
Belancio VP, Hedges DJ, Deininger P. LINE-1 RNA splicing and influences on mammalian gene expression. Nucleic Acids Res. 2006;34(5):1512–21.
Dewannieux M, Esnault C, Heidmann T. LINE-mediated retrotransposition of marked Alu sequences. Nat Genet. 2003;35(1):41–8.
Dewannieux M, Heidmann T. L1-mediated retrotransposition of murine B1 and B2 SINEs recapitulated in cultured cells. J Mol Biol. 2005;349(2):241–7.
Wallace N, Wagstaff BJ, Deininger PL, Roy-Engel AM. LINE-1 ORF1 protein enhances Alu SINE retrotransposition. Gene. 2008;419(1–2):1–6.
Vassetzky NS, Kramerov DA. SINEBase: a database and tool for SINE analysis. Nucleic Acids Res. 2013;41(Database issue):D83–9.
Deininger PL, Tiedge H, Kim J, Brosius J. Evolution, expression, and possible function of a master gene for amplification of an interspersed repeated DNA family in rodents. Prog Nucleic Acid Res Mol Biol. 1996;52:67–88.
Weiner AM. An abundant cytoplasmic 7S RNA is complementary to the dominant interspersed middle repetitive DNA sequence family in the human genome. Cell. 1980;22(1 Pt 1):209–18.
Ullu E, Tschudi C. Alu sequences are processed 7SL RNA genes. Nature. 1984;312(5990):171–2.
Geiduschek EP, Kassavetis GA. The RNA polymerase III transcription apparatus. J Mol Biol. 2001;310(1):1–26.
Schramm L, Hernandez N. Recruitment of RNA polymerase III to its target promoters. Genes Dev. 2002;16(20):2593–620.
Bird AP. DNA methylation and the frequency of CpG in animal DNA. Nucleic Acids Res. 1980;8(7):1499–504.
Hodgkinson A, Eyre-Walker A. Variation in the mutation rate across mammalian genomes. Nat Rev Genet. 2011;12(11):756–66.
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.
Magiorkinis G, Belshaw R, Katzourakis A. There and back again': revisiting the pathophysiological roles of human endogenous retroviruses in the post-genomic era. Philos Trans R Soc Lond Ser B Biol Sci. 2013;368(1626):20120504.
Stoye JP. Endogenous retroviruses: still active after all these years? Curr Biol. 2001;11(22):R914–6.
Luo ZX, Yuan CX, Meng QJ, Ji Q. A Jurassic eutherian mammal and divergence of marsupials and placentals. Nature. 2011;476(7361):442–5.
Clough JE, Foster JA, Barnett M, Wichman HA. Computer simulation of transposable element evolution: random template and strict master models. J Mol Evol. 1996;42(1):52–8.
Casavant NC, Hardies SC. The dynamics of murine LINE-1 subfamily amplification. J Mol Biol. 1994;241(3):390–7.
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.
Pascale E, Liu C, Valle E, Usdin K, Furano AV. The evolution of long interspersed repeated DNA (L1, LINE 1) as revealed by the analysis of an ancient rodent L1 DNA family. J Mol Evol. 1993;36(1):9–20.
Kramerov DA, Vassetzky NS. Short retroposons in eukaryotic genomes. Int Rev Cytol. 2005;247:165–221.
Ogiwara I, Miya M, Ohshima K, Okada N. Retropositional parasitism of SINEs on LINEs: identification of SINEs and LINEs in elasmobranchs. Mol Biol Evol. 1999;16(9):1238–50.
Goodier JL, Ostertag EM, Du K, Kazazian HH Jr. A novel active L1 retrotransposon subfamily in the mouse. Genome Res. 2001;11(10):1677–85.
Hedges DJ, Deininger PL. Inviting instability: transposable elements, double-strand breaks, and the maintenance of genome integrity. Mutat Res. 2007;616(1–2):46–59.
Belancio VP, Hedges DJ, Deininger P. Mammalian non-LTR retrotransposons: for better or worse, in sickness and in health. Genome Res. 2008;18(3):343–58.
Muotri AR, Chu VT, Marchetto MC, Deng W, Moran JV, Gage FH. Somatic mosaicism in neuronal precursor cells mediated by L1 retrotransposition. Nature. 2005;435(7044):903–10.
Coufal NG, Garcia-Perez JL, Peng GE, Yeo GW, Mu Y, Lovci MT, et al. L1 retrotransposition in human neural progenitor cells. Nature. 2009;460(7259):1127–31.
Cantrell MA, Carstens BC, Wichman HA. X chromosome inactivation and Xist evolution in a rodent lacking LINE-1 activity. PLoS One. 2009;4(7):e6252.
Chow JC, Ciaudo C, Fazzari MJ, Mise N, Servant N, Glass JL, et al. LINE-1 activity in facultative heterochromatin formation during X chromosome inactivation. Cell. 2010;141(6):956–69.
Sasaki T, Nishihara H, Hirakawa M, Fujimura K, Tanaka M, Kokubo N, et al. Possible involvement of SINEs in mammalian-specific brain formation. Proc Natl Acad Sci U S A. 2008;105(11):4220–5.
Kunarso G, Chia NY, Jeyakani J, Hwang C, Lu X, Chan YS, et al. Transposable elements have rewired the core regulatory network of human embryonic stem cells. Nat Genet. 2010;42(7):631–4.
Morrish TA, Gilbert N, Myers JS, Vincent BJ, Stamato TD, Taccioli GE, et al. DNA repair mediated by endonuclease-independent LINE-1 retrotransposition. Nat Genet. 2002;31(2):159–65.
Ariumi Y. Guardian of the human Genome: host defense mechanisms against LINE-1 Retrotransposition. Front Chem. 2016;4:28.
Gilbert C, Feschotte C. Horizontal acquisition of transposable elements and viral sequences: patterns and consequences. Curr Opin Genet Dev. 2018;49:15–24.
Cantrell MA, Scott L, Brown CJ, Martinez AR, Wichman HA. Loss of LINE-1 activity in the megabats. Genetics. 2008;178(1):393–404.
Yang L, Brunsfeld J, Scott L, Wichman H. Reviving the dead: history and reactivation of an extinct L1. PLoS Genet. 2014;10(6):e1004395.
Casavant NC, Scott L, Cantrell MA, Wiggins LE, Baker RJ, Wichman HA. The end of the LINE?: lack of recent L1 activity in a group of south American rodents. Genetics. 2000;154(4):1809–17.
Grahn RA, Rinehart TA, Cantrell MA, Wichman HA. Extinction of LINE-1 activity coincident with a major mammalian radiation in rodents. Cytogenet Genome Res. 2005;110(1–4):407–15.
Rinehart TA, Grahn RA, Wichman HA. SINE extinction preceded LINE extinction in sigmodontine rodents: implications for retrotranspositional dynamics and mechanisms. Cytogenet Genome Res. 2005;110(1–4):416–25.
Platt RN 2nd, Ray DA. A non-LTR retroelement extinction in Spermophilus tridecemlineatus. Gene. 2012;500(1):47–53.
Boissinot S, Roos C, Furano AV. Different rates of LINE-1 (L1) retrotransposon amplification and evolution in New World monkeys. J Mol Evol. 2004;58(1):122–30.
Waters PD, Dobigny G, Pardini AT, Robinson TJ. LINE-1 distribution in Afrotheria and Xenarthra: implications for understanding the evolution of LINE-1 in eutherian genomes. Chromosoma. 2004;113(3):137–44.
Gallus S, Hallstrom BM, Kumar V, Dodt WG, Janke A, Schumann GG, et al. Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil. Mol Biol Evol. 2015;32(5):1268–83.
Schenk JJ, Rowe KC, Steppan SJ. Ecological opportunity and incumbency in the diversification of repeated continental colonizations by muroid rodents. Syst Biol. 2013;62(6):837–64.
Smith MF, Patton JL. Phylogenetic relationships and the radiation of sigmodontine rodents in South America: evidence from cytochrome b. J Mamm Evol. 1999;6(2):89–128.
Wilson DE. Mammal species of the world: a taxonomic and geographic reference. Baltimore: Johns Hopkins University Press; 2005.
Cantrell MA, Ederer MM, Erickson IK, Swier VJ, Baker RJ, Wichman HA. MysTR: an endogenous retrovirus family in mammals that is undergoing recent amplifications to unprecedented copy numbers. J Virol. 2005;79(23):14698–707.
Erickson IK, Cantrell MA, Scott L, Wichman HA. Retrofitting the genome: L1 extinction follows endogenous retroviral expansion in a group of muroid rodents. J Virol. 2011;85(23):12315–23.
Cordaux R, Batzer MA. The impact of retrotransposons on human genome evolution. Nat Rev Genet. 2009;10(10):691–703.
Smit A, Hubley R. RepeatMasker Open-3.0; 1996-2010.
Marshall LG, Butler RF, Drake RE, Curtis GH, Tedford RH. Calibration of the great american interchange. Science. 1979;204(4390):272–9.
Cantrell MA, Grahn RA, Scott L, Wichman HA. Isolation of markers from recently transposed LINE-1 retrotransposons. Biotechniques. 2000;29(6):1310–6.
Gregory TR. Insertion-deletion biases and the evolution of genome size. Gene. 2004;324:15–34.
Smit AF, Toth G, Riggs AD, Jurka J. Ancestral, mammalian-wide subfamilies of LINE-1 repetitive sequences. J Mol Biol. 1995;246(3):401–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.
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.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.
Gregory TR. Animal Genome Size Database. 2014. http://www.genomesize.com.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinf. 2009;10:421.
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.
R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2013.
Price AL, Jones NC, Pevzner PA. De novo identification of repeat families in large genomes. Bioinformatics. 2005;21(Suppl 1):i351–8.
We thank Dr. Jerzy Jurka at the Genetic Information Research Institute for offering the bioinformatics training, John Brunsfeld and Dr. Celeste Brown on helpful ideas of the L1 reconstruction pipeline design, and Drs. Wenfeng An, Celeste Brown and James Foster for helpful comments and discussions. Materials for this study were provided by the Natural Science Research Laboratory in The Museum of Texas Tech University and the Texas Cooperative Wildlife Collection of Texas A&M University. IBEST Genomics Resources Core assisted with high-throughput sequencing and the IBEST Computer Resources Core provided resources for bioinformatics analyses.
This work was funded by National Institute of Health R01-GM38737 to HAW and National Science Foundation DDIG-1210694 to HAW and LY. Analytical resources were provided by National Institute of Health P20GM103408 and P30GM103324. Support for completing the analysis and manuscript preparation was provided by P20GM102420. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Data generated or analyzed during this study are included in this published article and its additional files.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Maximum likelihood phylogeny of detected L1 subfamilies. Reconstructed O. palustris and S. hispidus L1s, labeled ‘seed’, and P. maniculatus subfamilies 5 and 6 are included as markers. The tree was constructed using PhyML  with the GTR + I + G model and 100 bootstrap replicates. Bootstrap values > 80% are shown. (PDF 39 kb)
Figure S2. Divergence distribution of all detected L1 and B1 sequences. Percent divergence from the corresponding subfamily consensus sequences are plotted in 1% bins. Species and retrotransposon names are indicated at the top of each panel. (PDF 36 kb)
Table S1. Statistics and designation of L1 subfamilies and families. “Ory” stands for O. palustris and “Sig” stands for S. hispidus. “Peak” indicates the peak of the L1 divergence distribution of the subfamily or family identified by kernel smoothing. Copy numbers are normalized as copies per three Gbp of MiSeq sequence used for the search, which approximates the copy number per haploid genome. Designation of families is only shown after the first subfamily that belongs to it; all subsequent subfamilies belong to this family until the demarcation of the next family. Characters in family names: “S” represents S. hispidus-specific, “OS” for shared by O. palustris and S. hispidus and “OSP” for shared by O. palustris, S. hispidus and P. maniculatus. Numbers in the family names reflect their divergence among the family group with “1” being the youngest. Copy numbers of families are rounded sums of subfamily copy numbers per three Gbp of sequences and, thus, are occasionally off by one. (XLSX 14 kb)
Table S2. Statistics and designation of B1 subfamilies and families. “Ory” stands for O. palustris and “Sig” stands for S. hispidus. “Peak” indicates the peak of the B1 divergence distribution of the subfamily or family identified by kernel smoothing. Copy numbers are normalized as copies per three Gbp of MiSeq sequence used for the search. Designation of families is only shown after the first subfamily that belongs to it; all subsequent subfamilies belong to this family until the demarcation of the next family. Characters in family names: “OS” represents families shared by O. palustris and S. hispidus and “OSP” for families shared by O. palustris, S. hispidus and P. maniculatus. Numbers in the family names reflect their divergence within the family group with “1” being the youngest. Copy numbers of families are rounded sums of subfamily copy numbers per three Gbp of sequences.and S. hispidus genome. CpG sites with evidence of methylation-induced mutations were removed from the alignment. (XLSX 13 kb)
Alignment of L1 subfamily sequences used in the tables and figures. CpG sites with evidence of methylation-induced mutations were removed from the alignment. (FA 15 kb)
Alignment of B1 subfamily sequences used in the tables and figures. CpG sites with evidence of methylation-induced mutations were removed from the alignment. (FA 7 kb)
Alignment of L1 subfamily sequences based on RepeatScout-generated seeds and S. hispidus genome. CpG sites with evidence of methylation-induced mutations were removed from the alignment. (FA 21 kb)