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

A large percentage of the repetitive elements in mammalian genomes are retroelements, which have been moved primarily by LINE-1 retrotransposons and endogenous retroviruses. Although LINE-1 elements have remained active throughout the mammalian radiation, specific groups of endogenous retroviruses generally remain active for comparatively shorter periods of time. Identification of an unusual extinction of LINE-1 activity in a group of South American rodents has opened a window for examination of the interplay in mammalian genomes between these ubiquitous retroelements. In the course of a search for any type of repetitive sequences whose copy numbers have substantially changed in Oryzomys palustris, a species that has lost LINE-1 activity, versus Sigmodon hispidus, a closely related species retaining LINE-1 activity, we have identified an endogenous retrovirus family differentially amplified in these two species. Analysis of three full-length, recently transposed copies, called mysTR elements, revealed gag, pro, and pol coding regions containing stop codons which may have accumulated either before or after retrotransposition. Isolation of related sequences in S. hispidus and the LINE-1 active outgroup species, Peromyscus maniculatus, by PCR of a pro-pol region has allowed determination of copy numbers in each species. Unusually high copy numbers of approximately 10,000 in O. palustris versus 1,000 in S. hispidus and 4,500 in the more distantly related P. maniculatus leave open the question of whether there is a connection between endogenous retrovirus activity and LINE-1 inactivity. Nevertheless, these independent expansions of mysTR represent recent amplifications of this endogenous retrovirus family to unprecedented levels.     

Larger Mammalian Body Size Leads to Lower Retroviral Activity

Retroviruses have been infecting mammals for at least 100 million years, leaving descendants in host genomes known as endogenous retroviruses (ERVs). The abundance of ERVs is partly determined by their mode of replication, but it has also been suggested that host life history traits could enhance or suppress their activity. We show that larger bodied species have lower levels of ERV activity by reconstructing the rate of ERV integration across 38 mammalian species. Body size explains 37% of the variance in ERV integration rate over the last 10 million years, controlling for the effect of confounding due to other life history traits. Furthermore, 68% of the variance in the mean age of ERVs per genome can also be explained by body size. These results indicate that body size limits the number of recently replicating ERVs due to their detrimental effects on their host. To comprehend the possible mechanistic links between body size and ERV integration we built a mathematical model, which shows that ERV abundance is favored by lower body size and higher horizontal transmission rates. We argue that because retroviral integration is tumorigenic, the negative correlation between body size and ERV numbers results from the necessity to reduce the risk of cancer, under the assumption that this risk scales positively with body size. Our model also fits the empirical observation that the lifetime risk of cancer is relatively invariant among mammals regardless of their body size, known as Peto''s paradox, and indicates that larger bodied mammals may have evolved mechanisms to limit ERV activity.     

The Effect of Life History on Retroviral Genome Invasions

Endogenous retroviruses (ERV), or the remnants of past retroviral infections that are no longer active, are found in the genomes of most vertebrates, typically constituting approximately 10% of the genome. In some vertebrates, particularly in shorter-lived species like rodents, it is not unusual to find active endogenous retroviruses. In longer-lived species, including humans where substantial effort has been invested in searching for active ERVs, it is unusual to find them; to date none have been found in humans. Presumably the chance of detecting an active ERV infection is a function of the length of an ERV epidemic. Intuitively, given that ERVs or signatures of past ERV infections are passed from parents to offspring, we might expect to detect more active ERVs in species with longer generation times, as it should take more years for an infection to run its course in longer than in shorter lived species. This means the observation of more active ERV infections in shorter compared to longer-lived species is paradoxical. We explore this paradox using a modeling approach to investigate factors that influence ERV epidemic length. Our simple epidemiological model may explain why we find evidence of active ERV infections in shorter rather than longer-lived species.     

Cross-Species Transmission and Differential Fate of an Endogenous Retrovirus in Three Mammal Lineages

Endogenous retroviruses (ERVs) arise from retroviruses chromosomally integrated in the host germline. ERVs are common in vertebrate genomes and provide a valuable fossil record of past retroviral infections to investigate the biology and evolution of retroviruses over a deep time scale, including cross-species transmission events. Here we took advantage of a catalog of ERVs we recently produced for the bat Myotis lucifugus to seek evidence for infiltration of these retroviruses in other mammalian species (>100) currently represented in the genome sequence database. We provide multiple lines of evidence for the cross-ordinal transmission of a gammaretrovirus endogenized independently in the lineages of vespertilionid bats, felid cats and pangolin ~13–25 million years ago. Following its initial introduction, the ERV amplified extensively in parallel in both bat and cat lineages, generating hundreds of species-specific insertions throughout evolution. However, despite being derived from the same viral species, phylogenetic and selection analyses suggest that the ERV experienced different amplification dynamics in the two mammalian lineages. In the cat lineage, the ERV appears to have expanded primarily by retrotransposition of a single proviral progenitor that lost infectious capacity shortly after endogenization. In the bat lineage, the ERV followed a more complex path of germline invasion characterized by both retrotransposition and multiple infection events. The results also suggest that some of the bat ERVs have maintained infectious capacity for extended period of time and may be still infectious today. This study provides one of the most rigorously documented cases of cross-ordinal transmission of a mammalian retrovirus. It also illustrates how the same retrovirus species has transitioned multiple times from an infectious pathogen to a genomic parasite (i.e. retrotransposon), yet experiencing different invasion dynamics in different mammalian hosts.     

A New Approach to Establish a Cell Line with Reduced Risk of Endogenous Retroviruses

Endogenous retroviruses (ERVs) are integrated as DNA proviruses in the genomes of all mammalian species. Several ERVs are replication-competent and produced as fully infectious viruses from host cell. Thus, live-attenuated vaccines and biological substances have been prepared using the cell lines which may produce ERV. Indeed, we recently reported that several commercial live-attenuated vaccines for pets were contaminated with the infectious feline endogenous retrovirus, RD-114. In this study, to establish a cell line for vaccine manufacture with reduced risk of ERVs, we generated a cell line stably expressing human tetherin (Teth-CRFK cells). The release of infectious ERV from Teth-CRFK cells was suppressed to undetectable levels, while the production of parvovirus in Teth-CRFK cells was similar to that in parental CRFK cells. These observations suggest that Teth-CRFK cells will be useful as a cell line for the manufacture of live-attenuated vaccines or biological substances with reduced risk of ERV.     

Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains

The interspersed repeat content of mammalian genomes has been best characterized in human, mouse and cow. In this study, we carried out de novo identification of repeated elements in the equine genome and identified previously unknown elements present at low copy number. The equine genome contains typical eutherian mammal repeats, but also has a significant number of hybrid repeats in addition to clade-specific Long Interspersed Nuclear Elements (LINE). Equus caballus clade specific LINE 1 (L1) repeats can be classified into approximately five subfamilies, three of which have undergone significant expansion. There are 1115 full-length copies of these equine L1, but of the 103 presumptive active copies, 93 fall within a single subfamily, indicating a rapid recent expansion of this subfamily. We also analysed both interspersed and simple sequence repeats (SSR) genome-wide, finding that some repeat classes are spatially correlated with each other as well as with G+C content and gene density. Based on these spatial correlations, we have confirmed that recently-described ancestral vs. clade-specific genome territories can be defined by their repeat content. The clade-specific Short Interspersed Nuclear Element correlations were scattered over the genome and appear to have been extensively remodelled. In contrast, territories enriched for ancestral repeats tended to be contiguous domains. To determine if the latter territories were evolutionarily conserved, we compared these results with a similar analysis of the human genome, and observed similar ancestral repeat enriched domains. These results indicate that ancestral, evolutionarily conserved mammalian genome territories can be identified on the basis of repeat content alone. Interspersed repeats of different ages appear to be analogous to geologic strata, allowing identification of ancient vs. newly remodelled regions of mammalian genomes.     

Endogenous retroviruses and human evolution

Humans share about 99% of their genomic DNA with chimpanzees and bonobos; thus, the differences between these species are unlikely to be in gene content but could be caused by inherited changes in regulatory systems. Endogenous retroviruses (ERVs) comprise approximately 5% of the human genome. The LTRs of ERVs contain many regulatory sequences, such as promoters, enhancers, polyadenylation signals and factor-binding sites. Thus, they can influence the expression of nearby human genes. All known human-specific LTRs belong to the HERV-K (human ERV) family, the most active family in the human genome. It is likely that some of these ERVs could have integrated into regulatory regions of the human genome, and therefore could have had an impact on the expression of adjacent genes, which have consequently contributed to human evolution. This review discusses possible functional consequences of ERV integration in active coding regions.     

Recombination Marks the Evolutionary Dynamics of a Recently Endogenized Retrovirus

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.     

The first sequenced carnivore genome shows complex host-endogenous retrovirus relationships

Host-retrovirus interactions influence the genomic landscape and have contributed substantially to mammalian genome evolution. To gain further insights, we analyzed a female boxer (Canis familiaris) genome for complexity and integration pattern of canine endogenous retroviruses (CfERV). Intriguingly, the first such in-depth analysis of a carnivore species identified 407 CfERV proviruses that represent only 0.15% of the dog genome. In comparison, the same detection criteria identified about six times more HERV proviruses in the human genome that has been estimated to contain a total of 8% retroviral DNA including solitary LTRs. These observed differences in man and dog are likely due to different mechanisms to purge, restrict and protect their genomes against retroviruses. A novel group of gammaretrovirus-like CfERV with high similarity to HERV-Fc1 was found to have potential for active retrotransposition and possibly lateral transmissions between dog and human as a result of close interactions during at least 10.000 years. The CfERV integration landscape showed a non-uniform intra- and inter-chromosomal distribution. Like in other species, different densities of ERVs were observed. Some chromosomal regions were essentially devoid of CfERVs whereas other regions had large numbers of integrations in agreement with distinct selective pressures at different loci. Most CfERVs were integrated in antisense orientation within 100 kb from annotated protein-coding genes. This integration pattern provides evidence for selection against CfERVs in sense orientation relative to chromosomal genes. In conclusion, this ERV analysis of the first carnivorous species supports the notion that different mammals interact distinctively with endogenous retroviruses and suggests that retroviral lateral transmissions between dog and human may have occurred.     

Classification and nomenclature of endogenous retroviral sequences (ERVs): Problems and recommendations

The genomes of many species are crowded with repetitive mobile sequences. In the case of endogenous retroviruses (ERVs) there is, for various reasons, considerable confusion regarding names assigned to families/groups of ERVs as well as individual ERV loci. Human ERVs have been studied in greater detail, and naming of HERVs in the scientific literature is somewhat confusing not just to the outsider. Without guidelines, confusion for ERVs in other species will also probably increase if those ERVs are studied in greater detail. Based on previous experience, this review highlights some of the problems when naming and classifying ERVs, and provides some guidance for detecting and characterizing ERV sequences. Because of the close relationship between ERVs and exogenous retroviruses (XRVs) it is reasonable to reconcile their classification with that of XRVs. We here argue that classification should be based on a combination of similarity, structural features, (inferred) function, and previous nomenclature. Because the RepBase system is widely employed in genome annotation, RepBase designations should be considered in further taxonomic efforts. To lay a foundation for a phylogenetically based taxonomy, further analyses of ERVs in many hosts are needed. A dedicated, permanent, international consortium would best be suited to integrate and communicate our current and future knowledge on repetitive, mobile elements in general to the scientific community.     

Human endogenous retroviruses in the primate lineage and their influence on host genomes

Primates emerged about 60 million years ago. Since that time various primate-targeting retroviruses have integrated in the germ line of primate species, and some drifted to fixation. After germ line fixation, continued activity of proviruses resulted in intragenomic spread of so-called endogenous retroviruses (ERVs). Variant ERVs emerged, amplified in the genome and profoundly altered genome structures and potentially functionality. Importantly, ERVs are genome modifiers of exogenous origin. The human genome contains about 8% of sequences of retroviral origin. The human ERVs (HERVs) comprise many distinct families that amplified to copy numbers of up to several thousand. We review here the evolution of several well-characterized HERV families in the human lineage since initial germ line fixation. It is apparent that endogenous retroviruses profoundly affected the genomes of species in the evolutionary lineage leading to Homo sapiens.     

Phylogenetic Analysis Reveals That ERVs "Die Young" but HERV-H Is Unusually Conserved

About 8% of the human genome is made up of endogenous retroviruses (ERVs). Though most human endogenous retroviruses (HERVs) are thought to be irrelevant to our biology notable exceptions include members of the HERV-H family that are necessary for the correct functioning of stem cells. ERVs are commonly found in two forms, the full-length proviral form, and the more numerous solo-LTR form, thought to result from homologous recombination events. Here we introduce a phylogenetic framework to study ERV insertion and solo-LTR formation. We then apply the framework to site patterns sampled from a set of long alignments covering six primate genomes. Studying six categories of ERVs we quantitatively recapitulate patterns of insertional activity that are usually described in qualitative terms in the literature. A slowdown in most ERV groups is observed but we suggest that HERV-K activity may have increased in humans since they diverged from chimpanzees. We find that the rate of solo-LTR formation decreases rapidly as a function of ERV age and that an age dependent model of solo-LTR formation describes the history of ERVs more accurately than the commonly used exponential decay model. We also demonstrate that HERV-H loci are markedly less likely to form solo-LTRs than ERVs from other families. We conclude that the slower dynamics of HERV-H suggest a host role for the internal regions of these exapted elements and posit that in future it will be possible to use the relationship between full-length proviruses and solo-LTRs to help identify large scale co-options in distant vertebrate genomes.     

Identification of an ancient endogenous retrovirus,predating the divergence of the placental mammals

The evolutionary arms race between mammals and retroviruses has long been recognized as one of the oldest host–parasite interactions. Rapid evolution rates in exogenous retroviruses have often made accurate viral age estimations highly problematic. Endogenous retroviruses (ERVs), however, integrate into the germline of their hosts, and are subjected to their evolutionary rates. This study describes, for the first time, a retroviral orthologue predating the divergence of placental mammals, giving it a minimum age of 104–110 Myr. Simultaneously, other orthologous selfish genetic elements (SGEs), inserted into the ERV sequence, provide evidence for the oldest individual mammalian-wide interspersed repeat and medium-reiteration frequency interspersed repeat mammalian repeats, with the same minimum age. The combined use of shared SGEs and reconstruction of viral orthologies defines new limits and increases maximum ‘lookback’ times, with subsequent implications for the field of paleovirology.     

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.     

Low diversity,activity, and density of transposable elements in five avian genomes

In this study, we conducted the activity, diversity, and density analysis of transposable elements (TEs) across five avian genomes (budgerigar, chicken, turkey, medium ground finch, and zebra finch) to explore the potential reason of small genome sizes of birds. We found that these avian genomes exhibited low density of TEs by about 10% of genome coverages and low diversity of TEs with the TE landscapes dominated by CR1 and ERV elements, and contrasting proliferation dynamics both between TE types and between species were observed across the five avian genomes. Phylogenetic analysis revealed that CR1 clade was more diverse in the family structure compared with R2 clade in birds; avian ERVs were classified into four clades (alpha, beta, gamma, and ERV-L) and belonged to three classes of ERV with an uneven distributed in these lineages. The activities of DNA and SINE TEs were very low in the evolution history of avian genomes; most LINEs and LTRs were ancient copies with a substantial decrease of activity in recent, with only LTRs and LINEs in chicken and zebra finch exhibiting weak activity in very recent, and very few TEs were intact; however, the recent activity may be underestimated due to the sequencing/assembly technologies in some species. Overall, this study demonstrates low diversity, activity, and density of TEs in the five avian species; highlights the differences of TEs in these lineages; and suggests that the current and recent activity of TEs in avian genomes is very limited, which may be one of the reasons of small genome sizes in birds.     

Phylogenomic analysis of the L1 retrotransposons in Deuterostomia

L1 retrotransposons constitute the largest single component of mammalian genomes. In contrast to the single remaining lineage of L1 retrotransposons in mammalian genomes, some teleost fishes contain a highly diverse L1 retrotransposon repertoire. Major evolutionary changes in L1 retrotransposon repertoires have therefore taken place in the land vertebrates (Tetrapoda). The lack of sequence data for L1 retrotransposons in the basal living Tetrapoda lineages prompted an investigation of their distribution and evolution in the genomes of the key tetrapod lineages, amphibians and reptiles, and in lungfishes. In this study, we combined genome database searches with PCR analysis to demonstrate that L1 retrotransposons are present in the genomes of lungfishes, amphibians, and lepidosaurs. Phylogenomic analysis shows that the genomes of Deuterostomia possess three highly divergent groups of L1 retrotransposons, with distinct distribution patterns. The analysis of L1 diversity shows the presence of a very large number of diverse L1 families, each with very low copy numbers, at the time of the origin of tetrapods. During the evolution of synapsids, all but one L1 lineage have been lost. This study establishes that the loss of L1 diversity and explosion in copy numbers occurred in the synapsid ancestors of mammals, and was most probably caused by severe population bottlenecks.     

Bats with hATs: evidence for recent DNA transposon activity in genus Myotis

Transposable elements make up a significant fraction of many eukaryotic genomes. Although both classes of transposable elements, the DNA transposons and the retrotransposons, show substantial expansion in plants and invertebrates, the DNA transposons are thought to have become inactive in mammalian genomes long ago. Here, we report the first evidence for recent activity of DNA transposons in a mammalian lineage, the bat genus Myotis. Six recently active families of nonautonomous hobo/Activator/TAM transposons were identified in the Myotis lucifugus genome using computational tools. Low sequence divergence among the individual sequences and between individual sequences and their respective consensus sequences suggest their recent expansion in the M. lucifugus genome. Furthermore, amplification and sequencing of polymorphic insertion loci in a related taxon, M. austroriparius, confirms their recent activity. Myotis is one of the largest mammalian genera with 103 species. The discovery of DNA transposon activity in this genus may therefore influence our understanding of genome evolution and diversification in bats and in mammals in general. Furthermore, the identification of a likely autonomous element may lead to new approaches for mammalian genetic manipulation.