Two persistent LINE-1 lineages in Peromyscus have unequal rates of evolution

Flying Around in the Genome: Characterization of LINE-1 in Chiroptera

L1s are transposable elements that move by a copy-and-paste mechanism that continuously increases their copy number in the genome, such that each genome has a record of the L1 history in that host lineage. They make up about 20% of the genomes of eutherian mammals and have played a major role in shaping genome evolution. Chiroptera has the lowest average genome size among mammalian orders and the only documented case of L1 extinction affecting an entire mammalian family. Herein, L1 activity and extinction are characterized in all families of the order Chiroptera using a method that enriches for the youngest lineages of L1s in the genome. In addition to the previously reported L1 extinction in Pteropodidae, L1 extinction was documented to occur in Mormoops blainvilli, but this event did not affect all species of Mormoopidae. Further, there was no evidence of concordance between the evolution of L1s and their chiropteran host. There were two L1 lineages present before the divergence of all extant bats. Both lineages are extinct in the Pteropodidae. One or the other L1 lineage is extinct in almost all bat families, but Taphozous melanopogon maintains active members of both. Most intriguingly, some families within the Rhinolophoidea retain one active L1 lineage whereas other families retain the other, creating a deep discontinuity between L1 phylogeny and chiropteran phylogeny. These results indicate that there have been numerous losses of active L1 lineages over the history of chiropteran evolution, but that all chiropteran families except Pteropodidae have retained L1 activity.

Contrasted patterns of evolution of the LINE-1 retrotransposon in perissodactyls: the history of a LINE-1 extinction

Background

LINE-1 (L1) is the dominant autonomously replicating non-LTR retrotransposon in mammals. Although our knowledge of L1 evolution across the tree of life has considerably improved in recent years, what we know of L1 evolution in mammals is biased and comes mostly from studies in primates (mostly human) and rodents (mostly mouse). It is unclear if patterns of evolution that are shared between those two groups apply to other mammalian orders. Here we performed a detailed study on the evolution of L1 in perissodactyls by making use of the complete genome of the domestic horse and of the white rhinoceros. This mammalian order offers an excellent model to study the extinction of L1 since the rhinoceros is one of the few mammalian species to have lost active L1.

Results

We found that multiple L1 lineages, carrying different 5’UTRs, have been simultaneously active during the evolution of perissodactyls. We also found that L1 has continuously amplified and diversified in horse. In rhinoceros, L1 was very prolific early on. Two successful families were simultaneously active until ~20my ago but became extinct suddenly at exactly the same time.

Conclusions

The general pattern of L1 evolution in perissodactyls is very similar to what was previously described in mouse and human, suggesting some commonalities in the way mammalian genomes interact with L1. We confirmed the extinction of L1 in rhinoceros and we discuss several possible mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0117-4) contains supplementary material, which is available to authorized users.

LINE-1 retrotransposons: from 'parasite' sequences to functional elements

Long interspersed nuclear elements-1 (LINE-1) are the most abundant and active retrotransposons in the mammalian genomes. Traditionally, the occurrence of LINE-1 sequences in the genome of mammals has been explained by the selfish DNA hypothesis. Nevertheless, recently, it has also been argued that these sequences could play important roles in these genomes, as in the regulation of gene expression, genome modelling and X-chromosome inactivation. The non-random chromosomal distribution is a striking feature of these retroelements that somehow reflects its functionality. In the present study, we have isolated and analysed a fraction of the open reading frame 2 (ORF2) LINE-1 sequence from three rodent species, Cricetus cricetus, Peromyscus eremicus and Praomys tullbergi. Physical mapping of the isolated sequences revealed an interspersed longitudinal AT pattern of distribution along all the chromosomes of the complement in the three genomes. A detailed analysis shows that these sequences are preferentially located in the euchromatic regions, although some signals could be detected in the heterochromatin. In addition, a coincidence between the location of imprinted gene regions (as Xist and Tsix gene regions) and the LINE-1 retroelements was also observed. According to these results, we propose an involvement of LINE-1 sequences in different genomic events as gene imprinting, X-chromosome inactivation and evolution of repetitive sequences located at the heterochromatic regions (e.g. satellite DNA sequences) of the rodents' genomes analysed.

Jumping the fine LINE between species: horizontal transfer of transposable elements in animals catalyses genome evolution

Horizontal transfer (HT) is the transmission of genetic material between non-mating species, a phenomenon thought to occur rarely in multicellular eukaryotes. However, many transposable elements (TEs) are not only capable of HT, but have frequently jumped between widely divergent species. Here we review and integrate reported cases of HT in retrotransposons of the BovB family, and DNA transposons, over a broad range of animals spanning all continents. Our conclusions challenge the paradigm that HT in vertebrates is restricted to infective long terminal repeat (LTR) retrotransposons or retroviruses. This raises the possibility that other non-LTR retrotransposons, such as L1 or CR1 elements, believed to be only vertically transmitted, can horizontally transfer between species. Growing evidence indicates that the process of HT is much more general across different TEs and species than previously believed, and that it likely shapes eukaryotic genomes and catalyses genome evolution.

X chromosome inactivation and Xist evolution in a rodent lacking LINE-1 activity

Dosage compensation in eutherian mammals occurs by inactivation of one X chromosome in females. Silencing of that X chromosome is initiated by Xist, a large non-coding RNA, whose coating of the chromosome extends in cis from the X inactivation center. LINE-1 (L1) retrotransposons have been implicated as possible players for propagation of the Xist signal, but it has remained unclear whether they are essential components. We previously identified a group of South American rodents in which L1 retrotransposition ceased over 8 million years ago and have now determined that at least one species of these rodents, Oryzomys palustris, still retains X inactivation. We have also isolated and analyzed the majority of the Xist RNA from O. palustris and a sister species retaining L1 activity, Sigmodon hispidus, to determine if evolution in these sequences has left signatures that might suggest a critical role for L1 elements in Xist function. Comparison of rates of Xist evolution in the two species fails to support L1 involvement, although other explanations are possible. Similarly, comparison of known repeats and potential RNA secondary structures reveals no major differences with the exception of a new repeat in O. palustris that has potential to form new secondary structures.

Loss of LINE-1 activity in the megabats

LINE-1 (L1) retrotransposons are the most abundant type of mammalian retroelement. They have profound effects on genome plasticity and have been proposed to fulfill essential host functions, yet it remains unclear where they lie on the spectrum from parasitism to mutualism. Their ubiquity makes it difficult to determine the extent of their effects on genome evolution and gene expression because of the relative dearth of animal models lacking L1 activity. We have isolated L1 sequences from 11 megabat species by a method that enriches for recently inserted L1s and have done a bioinformatic examination of L1 sequences from a 12th species whose genome was recently shotgun sequenced. An L1 extinction event appears to have occurred at least 24 million years ago (MYA) in an ancestor of the megabats. The ancestor was unusual in having maintained two highly divergent long-term L1 lineages with different levels of activity, which appear, on an evolutionary scale, to have simultaneously lost that activity. These megabat species can serve as new animal models to ask what effect loss of L1 activity has on mammalian genome evolution and gene expression.

X accumulation of LINE-1 retrotransposons in Tokudaia osimensis, a spiny rat with the karyotype XO

The observation that LINE-1 transposable elements are enriched on the X in comparison to the autosomes led to the hypothesis that LINE-1s play a role in X chromosome inactivation. If this hypothesis is correct, loss of LINE-1 activity would be expected to result in species extinction or in an alternate pathway of dosage compensation. One such alternative pathway would be to evolve a karyotype that does not require dosage compensation between the sexes. Two of the three extant species of the Ryukyu spiny rat Tokudaia have such a karyotype; both males and females are XO. We asked whether this karyotype arose due to loss of LINE-1 activity and thus the loss of a putative component in the X inactivation pathway. Although XO Tokudaia has no need for dosage compensation, LINE-1s have been recently active in Tokudaia osimensis and show higher density on the lone X than on the autosomes.

Molecular evolution and tempo of amplification of human LINE-1 retrotransposons since the origin of primates

We investigated the evolution of the families of LINE-1 (L1) retrotransposons that have amplified in the human lineage since the origin of primates. We identified two phases in the evolution of L1. From approximately 70 million years ago (Mya) until approximately 40 Mya, three distinct L1 lineages were simultaneously active in the genome of ancestral primates. In contrast, during the last 40 million years (Myr), i.e., during the evolution of anthropoid primates, a single lineage of families has evolved and amplified. We found that novel (i.e., unrelated) regulatory regions (5'UTR) have been frequently recruited during the evolution of L1, whereas the two open-reading frames (ORF1 and ORF2) have remained relatively conserved. We found that L1 families coexisted and formed independently evolving L1 lineages only when they had different 5'UTRs. We propose that L1 families with different 5'UTR can coexist because they don't rely on the same host-encoded factors for their transcription and therefore do not compete with each other. The most prolific L1 families (families L1PA8 to L1PA3) amplified between 40 and 12 Mya. This period of high activity corresponds to an episode of adaptive evolution in a segment of ORF1. The correlation between the high activity of L1 families and adaptive evolution could result from the coevolution of L1 and a host-encoded repressor of L1 activity.

Extinction of LINE-1 activity coincident with a major mammalian radiation in rodents

LINE-1 transposable elements (L1s) are ubiquitous in mammals and are thought to have remained active since before the mammalian radiation. Only one L1 extinction event, in South American rodents in the genus Oryzomys, has been convincingly demonstrated. Here we examine the phylogenetic limits and evolutionary tempo of that extinction event by characterizing L1s in related rodents. Fourteen genera from five tribes within the Sigmodontinae subfamily were examined. Only the Sigmodontini, the most basal tribe in this group, demonstrate recent L1 activity. The Oryzomyini, Akodontini, Phyllotini, and Thomasomyini contain only L1s that appear to have inserted long ago; their L1s lack open reading frames, have mutations at conserved amino acid residues, and show numerous private mutations. They also lack restriction site-defined L1 subfamilies specific to any species, genus or tribe examined, and fail to form monophyletic species, genus or tribal L1 clusters. We determine here that this L1 extinction event occurred roughly 8.8 million years ago, near the divergence of Sigmodon from the remaining Sigmodontinae species. These species appear to be ideal model organisms for studying the impact of L1 inactivity on mammalian genomes.

Teleost fish genomes contain a diverse array of L1 retrotransposon lineages that exhibit a low copy number and high rate of turnover

Retrotransposable elements exhibit a wide range of variation in population dynamics, abundance, and lineage diversity among host genomes across taxa. This range of diversity is illustrated by a single well-defined constituent monophyletic clade of L1 non-LTR retrotransposons that is shared between mammalian and teleost fish genomes. Despite the clear phylogenetic relationships that exist between mammalian and teleost L1 sequences, these elements exhibit markedly different dynamics within their respective taxa. While mammalian genomes typically contain a single, abundant lineage of L1 elements that traces millions of years of evolution, the zebraflsh genome was recently shown to exhibit a high diversity of ancient lineages coexisting at a very low copy number and apparently exhibiting a high rate of turnover. In the present study, a combination of degenerate PCR, lineage-specific PCR, and genomic Southern blot analysis is utilized to demonstrate high L1 lineage diversity, low copy number, and a high proportion of polymorphic inserts in the genomes of the killifish species, Fundulus heteroclitus. Additional species surveyed by degenerate PCR include Cyprinodon variegatus, Rivulus marmoratus, and Menidia beryllina. These results further support the generality of the differences that exist in host-element dynamics between teleost fish and mammalian genomes with regard to L1 retrotransposons.

The Evolution of SINEs and LINEs in the genus Chironomus (Diptera)

Genomic DNA amplification from 51 species of the family Chironomidae shows that most contain relatives of NLRCth1 LINE and CTRT1 SINE retrotransposons first found in Chironomus thummi. More than 300 cloned PCR products were sequenced. The amplified region of the reverse transcriptase gene in the LINEs is intact and highly conserved, suggesting active elements. The SINEs are less conserved, consistent with minimal/no selection after transposition. A mitochondrial gene phylogeny resolves the Chironomus genus into six lineages (Guryev et al. 2001). LINE and SINE phylogenies resolve five of these lineages, indicating their monophyletic origin and vertical inheritance. However, both the LINE and the SINE tree topologies differ from the species phylogeny, resolving the elements into "clusters I-IV" and "cluster V" families. The data suggest a descent of all LINE and SINE subfamilies from two major families. Based on the species phylogeny, a few LINEs and a larger number of SINEs are cladisitically misplaced. Most misbranch with LINEs or SINEs from species with the same families of elements. From sequence comparisons, cladistically misplaced LINEs and several misplaced SINEs arose by convergent base substitutions. More diverged SINEs result from early transposition and some are derived from multiple source SINEs in the same species. SINEs from two species (C. dorsalis, C. pallidivittatus), expected to belong to the clusters I-IV family, branch instead with cluster V family SINEs; apparently both families predate separation of cluster V from clusters I-IV species. Correlation of the distribution of active SINEs and LINEs, as well as similar 3' sequence motifs in CTRT1 and NLRCth1, suggests coevolving retrotransposon pairs in which CTRT1 transposition depends on enzymes active during NLRCth1 LINE mobility.

An ancient retrovirus-like element contains hot spots for SINE insertion

Vertebrate retrotransposons have been used extensively for phylogenetic analyses and studies of molecular evolution. Information can be obtained from specific inserts either by comparing sequence differences that have accumulated over time in orthologous copies of that insert or by determining the presence or absence of that specific element at a particular site. The presence of specific copies has been deemed to be an essentially homoplasy-free phylogenetic character because the probability of multiple independent insertions into any one site has been believed to be nil. Mys elements are a type of LTR-containing retrotransposon present in Sigmodontine rodents. In this study we have shown that one particular insert, mys-9, is an extremely old insert present in multiple species of the genus Peromyscus. We have found that different copies of this insert show a surprising range of sizes, due primarily to a continuing series of SINE (short interspersed element) insertions into this locus. We have identified two hot spots for SINE insertion within mys-9 and at each hot spot have found that two independent SINE insertions have occurred at identical sites. These results have major repercussions for phylogenetic analyses based on SINE insertions, indicating the need for caution when one concludes that the existence of a SINE at a specific locus in multiple individuals is indicative of common ancestry. Although independent insertions at the same locus may be rare, SINE insertions are not homoplasy-free phylogenetic markers.

The biological properties and evolutionary dynamics of mammalian LINE-1 retrotransposons

Mammalian LINE-1 (L1) elements belong to the superfamily of autonomously replicating retrotransposable elements that lack the long terminal repeated (LTR) sequences typical of retroviruses and retroviral-like retrotransposons. The non-LTR superfamily is very ancient and L1-like elements are ubiquitous in nature, having been found in plants, fungi, invertebrates, and various vertebrate classes from fish to mammals. L1 elements have been replicating and evolving in mammals for at least the past 100 million years and now constitute 20% or more of some mammalian genomes. Therefore, L1 elements presumably have had a profound, perhaps defining, effect on the evolution, structure, and function of mammalian genomes. L1 elements contain regulatory signals and encode two proteins: one is an RNA-binding protein and the second one presumably functions as an integrase-replicase, because it has both endonuclease and reverse transcriptase activities. This work reviews the structure and biological properties of L1 elements, including their regulation, replication, evolution, and interaction with their mammalian hosts. Although each of these processes is incompletely understood, what is known indicates that they represent challenging and fascinating biological phenomena, the resolution of which will be essential for fully understanding the biology of mammals.

LINE-1 (L1) lineages in the mouse

Recently, a rapidly amplifying family of mouse LINE-1 (L1) has been identified and named T(F). The evolutionary context surrounding the derivation of the T(F) family was examined through phylogenetic analysis of sequences in the 3' portion of the repeat. The Mus musculus domesticus T(F) family was found to be the terminal subfamily of the previously identified L1Md4 lineage. The L1Md4 lineage joins the other prototypical mouse LINE-1 lineage (the L1MdA2 lineage) approximately 1 MYA at about the time of the common ancestor of M. m. domesticus, Mus spicilegus, and Mus spretus. However, the T(F) family from M. m. domesticus was found to join to the previously reported M. spretus Ms475 and Ms7024 LINE-1 families at just 0.5 MYA, indicating horizontal transfer. The T(F) family from M. m. domesticus was then found to be even more recently related to LINE-1's from another species, M. spicilegus. A separate spretus A2 lineage was found through a directed search of a PCR library. This lineage, in contrast to the spretus T(F) lineage, does join domesticus at about 1 MYA, as would be expected in the absence of horizontal transfer. A third major family was also found that splits off from the L1Md4 lineage shortly after its departure from the L1MdA2 lineage. The new family, named the Z family, was found to contain the de novo LINE-1 inserts causing the beige and med mutations. Whether the split with the Z family was before or after the recombination that introduced the F-type promoters and defined the inception of T(F) as a lineage is unclear. In enumerating copies of the various LINE-1 families, we found that T(F) 3' ends were not much more numerous than the reported number of 5' ends, suggesting that T(F) may not be subjected to the 90% truncation pattern typical of LINE-1 as a whole.

The end of the LINE?: lack of recent L1 activity in a group of South American rodents

L1s (LINE-1: Long Interspersed Nuclear Element 1) are present in all mammals examined to date. They occur in both placental mammals and marsupials and thus are thought to have been present in the genome prior to the mammalian radiation. This unusual conservation of a transposable element family for over 100 million years has led to speculation that these elements provide an advantage to the genomes they inhabit. We have recently identified a group of South American rodents, including rice rats (Oryzomys), in which L1s appear to be quiescent or extinct. Several observations support this conclusion. First, genomic Southern blot analysis fails to reveal genus-specific bands in Oryzomys. Second, we were unable to find recently inserted elements. Procedures to enrich for young elements did not yield any with an intact open reading frame for reverse transcriptase; all elements isolated had numerous insertions, deletions, and stop codons. Phylogenetic analysis failed to yield species-specific clusters among the L1 elements isolated, and all Oryzomys sequences had numerous private mutations. Finally, in situ hybridization of L1 to Oryzomys chromosomes failed to reveal the characteristic L1 distribution in Oryzomys with either a homologous or heterologous probe. Thus, Oryzomys is a viable candidate for L1 extinction from a mammalian host.

Molecular evolution of two lineages of L1 (LINE-1) retrotransposons in the california mouse, Peromyscus californicus

The large number of L1 [long interspersed elements (LINE)-1] sequences found in the genome is due to the insertion of copies of the retrotransposon over evolutionary time. The majority of copies appear to be replicates of a few active, or "master" templates. A continual replacement of master templates over time gives rise to lineages distinguishable by their own unique set of shared-sequence variants. A previous analysis of L1 sequences in deer mice, Peromyscus maniculatus and P. leucopus, revealed two active L1 lineages, marked by different rates of evolution, whose most recent common ancestor predates the expansion of the Peromyscus species. Here we exploit lineage-specific, shared-sequence variants to reveal a paucity of Lineage 2 sequences in at least one species, P. californicus. The dearth of Lineage 2 copies in P. californicus suggests that Lineage 2 may have been unproductive until after the most recent common ancestor of P. californicus and P. maniculatus. We also show that Lineage 1 appears to have a higher rate of evolution in P. maniculatus relative to either P. californicus or P. leucopus. As a phylogenetic tool, L1 lineage-specific variants support a close affinity between P. californicus and P. eremicus relative to the other species examined.

Deer Mice As Laboratory Animals

Although laboratory mice (Mus) and rats (Rattus) are the most widely used research rodents, deer mice (Peromyscus maniculatus) and their congeneric species are favored as nontraditional alternatives for some purposes. Mice of the native genus Peromyscus are the most abundant and widely distributed rodents in North America. They occur in a great diversity of habitats and play a significant role in natural ecosystems. Because of their abundance, peromyscines are commonly hosts for larva of ticks that transmit Lyme disease bacteria, and they are implicated in several other vector-borne diseases. Deer mice also are the principal carriers of the virus that causes hantaviral pulmonary syndrome, or "Four Corners disease." Deer mice are useful as laboratory models for a variety of other types of pure and applied research. They are easily maintained and bred in captivity using the husbandry protocols developed for other small laboratory rodent species. The Peromyscus Genetic Stock Center at the University of South Carolina maintains more than 50 laboratory-bred, well-characterized stocks of deer mice and other peromyscine species for research and educational use.

Nature of B chromosomes in the harvest mouse Reithrodontomys megalotis by fluorescence in situ hybridization (FISH)

Using fluorescence in situ hybridization, we examined the characteristics of two types of B chromosomes in harvest mice of the genus Reithrodontomys. B chromosomes were interrogated with rDNA, telomeric repeat, LINE element and centromeric heterochromatin probes. The two types of B chromosomes share the following features: (a) telomeres present on the ends of both arms; (b) hybridization to LINE probes; (c) absence of hybridization to the ribosomal gene probes; (d) C-band-positive centromeric regions; and (e) euchromatic arms. They differ as follows: (a) the larger B element hybridizes to the centromeric heterochromatin (pMeg-1) probe whereas the smaller B element does not; (b) the amount of C-band-positive material is reduced in the smaller B chromosome relative to that present on the larger B chromosome; and (c) the smaller element is reduced in size by about a third. It is concluded that the larger B chromosome arose as a leftover centromere from centric fusion, whereas the smaller element has a different origin perhaps as an intact fragment or as an amplified region from the A chromosomes. The presence of euchromatic regions on B chromosomes may account for their survival in the karyotype.

Analysis of orthologous retrovirus-like elements in the white-footed mouse, Peromyscus leucopus

Three loci in the genome of the white-footed mouse, Peromyscus leucopus, were examined for the presence or absence of orthologous copies of the retrovirus-like element mys using polymerase chain reaction. We examined these loci in 28 mice collected throughout the P. leucopus species range. Mys insertions were present in only one of the individuals examined at the mys-1 and mys-7 loci. Conversely, the mys-6 element was found in several individuals, but the presence of this element was limited to northern latitudes. Because the long terminal repeats (LTRs) of a given element are expected to be identical at the time of retrotransposition into the genome, and to accumulate changes over evolutionary time, within-element LTR sequence comparisons can be used to estimate the relative age of insertions. Within-element LTR differences are greater in mys-6 than in mys-1 or mys-7. The LTRs from orthologous mys-6 elements of six mice were sequenced. The alignment revealed 13 of the 22 differences between the right and left LTRs that were shared by all orthologous mys-6 sites, suggesting that relative to its time of insertion into the genome, mys-6 has only recently spread across the northern part of the species range.