MysTR: an endogenous retrovirus family in mammals that is undergoing recent amplifications to unprecedented copy numbers

Karyotypic Reshuffling in the Genus Rhipidomys (Rodentia: Cricetidae: Sigmodontinae) Revealed by Zoo-FISH

INTRODUCTION

Rhipidomys is the second most specious and the most widespread genus of the tribe Thomasomyini. Chromosomal data have been an important tool in the taxonomy of the group that presents low variability of diploid number (2n) and highly variable fundamental numbers (FNs). Despite such diversity, the genus has been studied mainly by classical and banding cytogenetic techniques.

METHODS

This study performed a comparative study between R. emiliae (2n = 44, FN = 52), R. macrurus (2n = 44, FN = 49), R. nitela (2n = 50, FN = 71), and R. mastacalis (2n = 44, FN = 72) using chromosome painting probes of two Oryzomyini species.

RESULTS

Our analysis revealed pericentric inversion as the main rearrangement involved in the karyotype evolution of the group, although tandem fusions/fissions were also detected. In addition, we detected eight syntenic associations exclusive of the genus Rhipidomys, and three syntenic associations shared between species of the tribe Thomasomyini and Oryzomyini.

CONCLUSION

Comparative cytogenetic analysis by ZOO-FISH on genus Rhipidomys supports a pattern of chromosomal rearrangement already suggested by comparative G-banding. However, the results suggest that karyotype variability in the genus could also involve the occurrence of an evolutionary new centromere.

Characterization of repetitive DNA on the genome of the marsh rat Holochilus nanus (Cricetidae: Sigmodontinae)

Repetitive DNA are sequences repeated hundreds or thousands of times and an abundant part of eukaryotic genomes. SatDNA represents the majority of the repetitive sequences, followed by transposable elements. The species Holochilus nanus (HNA) belongs to the rodent tribe Oryzomyini, the most taxonomically diverse of Sigmodontinae subfamily. Cytogenetic studies on Oryzomyini reflect such diversity by revealing an exceptional range of karyotype variability. However, little is known about the repetitive DNA content and its involvement in chromosomal diversification of these species. In the search for a more detailed understanding about the composition of repetitive DNA on the genome of HNA and other species of Oryzomyini, we employed a combination of bioinformatic, cytogenetic and molecular techniques to characterize the repetitive DNA content of these species. RepeatExplorer analysis showed that almost half of repetitive content of HNA genome are composed by Long Terminal Repeats and a less significant portion are composed by Short Interspersed Nuclear Elements and Long Interspersed Nuclear Elements. RepeatMasker showed that more than 30% of HNA genome are composed by repetitive sequences, with two main waves of repetitive element insertion. It was also possible to identify a satellite DNA sequence present in the centromeric region of Oryzomyini species, and a repetitive sequence enriched on the long arm of HNA X chromosome. Also, comparative analysis between HNA genome with and without B chromosome did not evidence any repeat element enriched on the supernumerary, suggesting that B chromosome of HNA is composed by a fraction of repeats from all the genome.

An ancient retroviral RNA element hidden in mammalian genomes and its involvement in co-opted retroviral gene regulation

Background

Retroviruses utilize multiple unique RNA elements to control RNA processing and translation. However, it is unclear what functional RNA elements are present in endogenous retroviruses (ERVs). Gene co-option from ERVs sometimes entails the conservation of viral cis-elements required for gene expression, which might reveal the RNA regulation in ERVs.

Results

Here, we characterized an RNA element found in ERVs consisting of three specific sequence motifs, called SPRE. The SPRE-like elements were found in different ERV families but not in any exogenous viral sequences examined. We observed more than a thousand of copies of the SPRE-like elements in several mammalian genomes; in human and marmoset genomes, they overlapped with lineage-specific ERVs. SPRE was originally found in human syncytin-1 and syncytin-2. Indeed, several mammalian syncytin genes: mac-syncytin-3 of macaque, syncytin-Ten1 of tenrec, and syncytin-Car1 of Carnivora, contained the SPRE-like elements. A reporter assay revealed that the enhancement of gene expression by SPRE depended on the reporter genes. Mutation of SPRE impaired the wild-type syncytin-2 expression while the same mutation did not affect codon-optimized syncytin-2, suggesting that SPRE activity depends on the coding sequence.

Conclusions

These results indicate multiple independent invasions of various mammalian genomes by retroviruses harboring SPRE-like elements. Functional SPRE-like elements are found in several syncytin genes derived from these retroviruses. This element may facilitate the expression of viral genes, which were suppressed due to inefficient codon frequency or repressive elements within the coding sequences. These findings provide new insights into the long-term evolution of RNA elements and molecular mechanisms of gene expression in retroviruses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-021-00580-2.

The genome of Peromyscus leucopus, natural host for Lyme disease and other emerging infections

The rodent Peromyscus leucopus is the natural reservoir of several tick-borne infections, including Lyme disease. To expand the knowledge base for this key species in life cycles of several pathogens, we assembled and scaffolded the P. leucopus genome. The resulting assembly was 2.45 Gb in total length, with 24 chromosome-length scaffolds harboring 97% of predicted genes. RNA sequencing following infection of P. leucopus with Borreliella burgdorferi, a Lyme disease agent, shows that, unlike blood, the skin is actively responding to the infection after several weeks. P. leucopus has a high level of segregating nucleotide variation, suggesting that natural resistance alleles to Crispr gene targeting constructs are likely segregating in wild populations. The reference genome will allow for experiments aimed at elucidating the mechanisms by which this widely distributed rodent serves as natural reservoir for several infectious diseases of public health importance, potentially enabling intervention strategies.

Tracing the history of LINE and SINE extinction in sigmodontine rodents

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0164-5) contains supplementary material, which is available to authorized users.

Birth, School, Work, Death, and Resurrection: The Life Stages and Dynamics of Transposable Element Proliferation

Transposable elements (TEs) can be maintained in sexually reproducing species even if they are harmful. However, the evolutionary strategies that TEs employ during proliferation can modulate their impact. In this review, I outline the different life stages of a TE lineage, from birth to proliferation to extinction. Through their interactions with the host, TEs can exploit diverse strategies that range from long-term coexistence to recurrent movement across species boundaries by horizontal transfer. TEs can also engage in a poorly understood phenomenon of TE resurrection, where TE lineages can apparently go extinct, only to proliferate again. By determining how this is possible, we may obtain new insights into the evolutionary dynamics of TEs and how they shape the genomes of their hosts.

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.

Comparative Genomic In Situ Hybridization and the Possible Role of Retroelements in the Karyotypic Evolution of Three Akodontini Species

South American Akodontini rodents are characterized by a large number of chromosome rearrangements. Among them, the genus Akodon has been extensively analyzed with classical and molecular cytogenetics, which allowed the identification of a large number of intra- and interspecific chromosomal variation due to Robertsonian rearrangements, pericentric inversions, and heterochromatin additions/deletions. In order to shed some light on the cause of these rearrangements, we comparatively analyzed the karyotypes of three Akodontini species, Akodon cursor (2n = 14, FN = 19), A. montensis (2n = 24, FN = 42), and Necromys lasiurus (2n = 34, FN = 34), after GTG- and CBG-banding. The karyotypes differed by Robertsonian rearrangements, pericentric inversions, centromere repositioning, and heterochromatin variation. Genome comparisons were performed through interspecific fluorescent in situ hybridization (FISH) with total genomic DNAs of each species as probes (GISH). Our results revealed considerable conservation of the euchromatic portions among the three karyotypes suggesting that they mostly differ in their heterochromatic regions. FISH was also performed to assess the distribution of telomeric sequences, long and short interspersed repetitive elements (LINE-1 and B1 SINE) and of the endogenous retrovirus mysTR in the genomes of the three species. The results led us to infer that transposable elements have played an important role in the enormous chromosome variation found in Akodontini.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Positive selection and multiple losses of the LINE-1-derived L1TD1 gene in mammals suggest a dual role in genome defense and pluripotency

Mammalian genomes comprise many active and fossilized retroelements. The obligate requirement for retroelement integration affords host genomes an opportunity to 'domesticate' retroelement genes for their own purpose, leading to important innovations in genome defense and placentation. While many such exaptations involve retroviruses, the L1TD1 gene is the only known domesticated gene whose protein-coding sequence is almost entirely derived from a LINE-1 (L1) retroelement. Human L1TD1 has been shown to play an important role in pluripotency maintenance. To investigate how this role was acquired, we traced the origin and evolution of L1TD1. We find that L1TD1 originated in the common ancestor of eutherian mammals, but was lost or pseudogenized multiple times during mammalian evolution. We also find that L1TD1 has evolved under positive selection during primate and mouse evolution, and that one prosimian L1TD1 has 'replenished' itself with a more recent L1 ORF1 from the prosimian genome. These data suggest that L1TD1 has been recurrently selected for functional novelty, perhaps for a role in genome defense. L1TD1 loss is associated with L1 extinction in several megabat lineages, but not in sigmodontine rodents. We hypothesize that L1TD1 could have originally evolved for genome defense against L1 elements. Later, L1TD1 may have become incorporated into pluripotency maintenance in some lineages. Our study highlights the role of retroelement gene domestication in fundamental aspects of mammalian biology, and that such domesticated genes can adopt different functions in different lineages.

Prototype endogenous avian retroviruses of the genus Gallus

Ancient endogenous retroviruses (ERVs), designated endogenous avian retrovirus (EAVs), are present in all Gallus spp. including the chicken, and resemble the modern avian sarcoma and leukosis viruses (ASLVs). The EAVs comprise several distinct retroviruses, including EAV-0, EAV-E51 and EAV-HP, as well as a putative member previously named the avian retrotransposon of chickens (ART-CH). Thus far, only the EAV-HP elements have been well characterized. Here, we determined sequences of representative EAV-0 and EAV-E51 proviruses by cloning and data mining of the 2011 assembly of the Gallus gallus genome. Although the EAV-0 elements are primarily deleted in the env region, we identified two complete EAV-0 env genes within the G. gallus genome and prototype elements sharing identity with an EAV-E51-related clone previously designated EAV-E33. Prototype EAV-0, EAV-E51 and EAV-E33 gag, pol and env gene sequences used for phylogenetic analysis of deduced proteins showed that the EAVs formed three distinct clades, with EAV-0 sharing the last common ancestor with the ASLVs. The EAV-E51 clade showed the greatest level of divergence compared with other EAVs or ASLVs, suggesting that these ERVs represented exogenous retroviruses that evolved and integrated into the germline over a long period of time. Moreover, the degree of divergence between the chicken and red jungle fowl EAV-E51 sequences suggested that they were more ancient than the other EAVs and may have diverged through mutations that accumulated post-integration. Finally, we showed that the ART-CH elements were chimeric defective ERVs comprising portions of EAV-E51 and EAV-HP rather than authentic retrotransposons.

Transposable elements and viruses as factors in adaptation and evolution: an expansion and strengthening of the TE-Thrust hypothesis

In addition to the strong divergent evolution and significant and episodic evolutionary transitions and speciation we previously attributed to TE-Thrust, we have expanded the hypothesis to more fully account for the contribution of viruses to TE-Thrust and evolution. The concept of symbiosis and holobiontic genomes is acknowledged, with particular emphasis placed on the creativity potential of the union of retroviral genomes with vertebrate genomes. Further expansions of the TE-Thrust hypothesis are proposed regarding a fuller account of horizontal transfer of TEs, the life cycle of TEs, and also, in the case of a mammalian innovation, the contributions of retroviruses to the functions of the placenta. The possibility of drift by TE families within isolated demes or disjunct populations, is acknowledged, and in addition, we suggest the possibility of horizontal transposon transfer into such subpopulations. “Adaptive potential” and “evolutionary potential” are proposed as the extremes of a continuum of “intra-genomic potential” due to TE-Thrust. Specific data is given, indicating “adaptive potential” being realized with regard to insecticide resistance, and other insect adaptations. In this regard, there is agreement between TE-Thrust and the concept of adaptation by a change in allele frequencies. Evidence on the realization of “evolutionary potential” is also presented, which is compatible with the known differential survivals, and radiations of lineages. Collectively, these data further suggest the possibility, or likelihood, of punctuated episodes of speciation events and evolutionary transitions, coinciding with, and heavily underpinned by, intermittent bursts of TE activity.

Retrofitting the genome: L1 extinction follows endogenous retroviral expansion in a group of muroid rodents

Long interspersed nuclear element 1 (LINE-1; L1) retrotransposons are the most common retroelements in mammalian genomes. Unlike individual families of endogenous retroviruses (ERVs), they have remained active throughout the mammalian radiation and are responsible for most of the retroelement movement and much genome rearrangement within mammals. They can be viewed as occupying a substantial niche within mammalian genomes. Our previous demonstration that L1s and B1 short interspersed nuclear elements (SINEs) are inactive in a group of South American rodents led us to ask if other elements have amplified to fill the empty niche. We identified a novel and highly active family of ERVs (mysTR). To determine whether loss of L1 activity was correlated with expansion of mysTR, we examined mysTR activity in four South American rodent species that have lost L1 and B1 activity and four sister species with active L1s. The copy number of recent mysTR insertions was extremely high, with an average of 4,200 copies per genome. High copy numbers exist in both L1-active and L1-extinct species, so the mysTR expansion appears to have preceded the loss of both SINE and L1 activity rather than to have filled an empty niche created by their loss. It may be coincidental that two unusual genomic events--loss of L1 activity and massive expansion of an ERV family--occur in the same group of mammals. Alternatively, it is possible that this large ERV expansion set the stage for L1 extinction.

Recent amplification of the kangaroo endogenous retrovirus, KERV, limited to the centromere

Mammalian retrotransposons, transposable elements that are processed through an RNA intermediate, are categorized as short interspersed elements (SINEs), long interspersed elements (LINEs), and long terminal repeat (LTR) retroelements, which include endogenous retroviruses. The ability of transposable elements to autonomously amplify led to their initial characterization as selfish or junk DNA; however, it is now known that they may acquire specific cellular functions in a genome and are implicated in host defense mechanisms as well as in genome evolution. Interactions between classes of transposable elements may exert a markedly different and potentially more significant effect on a genome than interactions between members of a single class of transposable elements. We examined the genomic structure and evolution of the kangaroo endogenous retrovirus (KERV) in the marsupial genus Macropus. The complete proviral structure of the kangaroo endogenous retrovirus, phylogenetic relationship among relative retroviruses, and expression of this virus in both Macropus rufogriseus and M. eugenii are presented for the first time. In addition, we show the relative copy number and distribution of the kangaroo endogenous retrovirus in the Macropus genus. Our data indicate that amplification of the kangaroo endogenous retrovirus occurred in a lineage-specific fashion, is restricted to the centromeres, and is not correlated with LINE depletion. Finally, analysis of KERV long terminal repeat sequences using massively parallel sequencing indicates that the recent amplification in M. rufogriseus is likely due to duplications and concerted evolution rather than a high number of independent insertion events.

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.