Endogenous retroviruses of non-avian/mammalian vertebrates illuminate diversity and deep history of retroviruses

Signature of viral fossils: a comparative genomics approach to understand the diversity of endogenous retroviruses in bats

Endogenous retroviruses (ERVs) are traces of past viral infections commonly found in vertebrate genomes. Many ERVs are tightly regulated by the host genomes and co-opted for various functions within the hosts. Bats are the only true volant mammals, with the smallest mammalian genomes and a high fraction of ERVs within the genomes. They are important hosts for various zoonotic viral pathogens and can effectively modulate their immune response to tolerate viral infections. Integrations of retroviruses have been implicated as one of the mechanisms by which bats have co-evolved strategies to combat viral infections. In this study, we investigated the diversity of ERVs in over 40 publicly available bat genomes to understand the distribution and the evolution of ERVs within bats. We observed all classes of ERVs within bat genomes including even the complex lenti retroviruses. Alpha and spuma retroviruses which are generally considered rare in mammals, were common within bats. We observed a positive correlation between bat genome size and length of ERV elements. Interestingly, nearly 30 % of the ERVs within bats are intact suggesting a recent origin or co-option by the host genome. Future studies focusing on comparative genomic and experimental data will be critical to understand the role of these ERVs in host genome evolution.

A novel approach to exploring the dark genome and its application to mapping of the vertebrate virus fossil record

BACKGROUND

Genomic regions that remain poorly understood, often referred to as the dark genome, contain a variety of functionally relevant and biologically informative features. These include endogenous viral elements (EVEs)-virus-derived sequences that can dramatically impact host biology and serve as a virus fossil record. In this study, we introduce a database-integrated genome screening (DIGS) approach to investigate the dark genome in silico, focusing on EVEs found within vertebrate genomes.

RESULTS

Using DIGS on 874 vertebrate genomes, we uncover approximately 1.1 million EVE sequences, with over 99% originating from endogenous retroviruses or transposable elements that contain EVE DNA. We show that the remaining 6038 sequences represent over a thousand distinct horizontal gene transfer events across 10 virus families, including some that have not previously been reported as EVEs. We explore the genomic and phylogenetic characteristics of non-retroviral EVEs and determine their rates of acquisition during vertebrate evolution. Our study uncovers novel virus diversity, broadens knowledge of virus distribution among vertebrate hosts, and provides new insights into the ecology and evolution of vertebrate viruses.

CONCLUSIONS

We comprehensively catalog and analyze EVEs within 874 vertebrate genomes, shedding light on the distribution, diversity, and long-term evolution of viruses and reveal their extensive impact on vertebrate genome evolution. Our results demonstrate the power of linking a relational database management system to a similarity search-based screening pipeline for in silico exploration of the dark genome.

Discovery and genetic characterization of novel paramyxoviruses from small mammals in Hubei Province, Central China

Paramyxoviruses are a group of single-stranded, negative-sense RNA viruses, some of which are responsible for acute human disease, including parainfluenza virus, measles virus, Nipah virus and Hendra virus. In recent years, a large number of novel paramyxoviruses, particularly members of the genus Jeilongvirus, have been discovered in wild mammals, suggesting that the diversity of paramyxoviruses may be underestimated. Here we used hemi-nested reverse transcription PCR to obtain 190 paramyxovirus sequences from 969 small mammals in Hubei Province, Central China. These newly identified paramyxoviruses were classified into four clades: genera Jeilongvirus, Morbillivirus, Henipavirus and Narmovirus, with most of them belonging to the genus Jeilongvirus. Using Illumina sequencing and Sanger sequencing, we successfully recovered six near-full-length genomes with different genomic organizations, revealing the more complex genome content of paramyxoviruses. Co-divergence analysis of jeilongviruses and their known hosts indicates that host-switching occurred more frequently in the evolutionary histories of the genus Jeilongvirus. Together, our findings demonstrate the high prevalence of paramyxoviruses in small mammals, especially jeilongviruses, and highlight the diversity of paramyxoviruses and their genome content, as well as the evolution of jeilongviruses.

Epigenetic Restriction Factors (eRFs) in Virus Infection

The ongoing arms race between viruses and their hosts is constantly evolving. One of the ways in which cells defend themselves against invading viruses is by using restriction factors (RFs), which are cell-intrinsic antiviral mechanisms that block viral replication and transcription. Recent research has identified a specific group of RFs that belong to the cellular epigenetic machinery and are able to restrict the gene expression of certain viruses. These RFs can be referred to as epigenetic restriction factors or eRFs. In this review, eRFs have been classified into two categories. The first category includes eRFs that target viral chromatin. So far, the identified eRFs in this category include the PML-NBs, the KRAB/KAP1 complex, IFI16, and the HUSH complex. The second category includes eRFs that target viral RNA or, more specifically, the viral epitranscriptome. These epitranscriptomic eRFs have been further classified into two types: those that edit RNA bases-adenosine deaminase acting on RNA (ADAR) and pseudouridine synthases (PUS), and those that covalently modify viral RNA-the N6-methyladenosine (m6A) writers, readers, and erasers. We delve into the molecular machinery of eRFs, their role in limiting various viruses, and the mechanisms by which viruses have evolved to counteract them. We also examine the crosstalk between different eRFs, including the common effectors that connect them. Finally, we explore the potential for new discoveries in the realm of epigenetic networks that restrict viral gene expression, as well as the future research directions in this area.

Comparison of Endogenous Alpharetroviruses (ALV-like) across Galliform Species: New Distant Proviruses

The Genus Alpharetrovirus contains viruses pathogenic mainly for chickens, forming the Avian Sarcoma and Leukosis Virus group (ASLV). Cells of most Galliform species, besides chickens, contain genetic elements (endogenous retroviruses, ERVs) that could recombine with other alpharetroviruses or express proteins, complementing defective ASLV, which may successfully replicate and cause disease. However, they are quite unknown, and only ALV-F, from ring-necked pheasants, has been partially published. Upon scrutiny of 53 genomes of different avian species, we found Alpharetrovirus-like sequences only in 12 different Galliformes, including six full-length (7.4-7.6 Kbp) and 27 partial sequences. Phylogenetic studies of the regions studied (LTR, gag, pol, and env) consistently resulted in five almost identical clades containing the same ERVs: Clade I (presently known ASLVs); Clade II (Callipepla spp. ERVs); Clade IIIa (Phasianus colchicus ERVs); Clade IIIb (Alectoris spp. ERVs); and Clade IV (Centrocercus spp. ERVs). The low pol identity scores suggested that each of these Clades may be considered a different species. ORF analysis revealed that putatively encoded proteins would be very similar in length and domains to those of other alpharetroviruses and thus potentially functional. This will undoubtedly contribute to better understanding the biology of defective viruses, especially in wild Galliformes, their evolution, and the danger they may represent for other wild species and the poultry industry.

Genealogical Diversity of Endogenous Retrovirus in the Jawless Fish Genome

Retroviral integration into ancient vertebrate genomes left traces that can shed light on the early history of viruses. In this study, we explored the early evolution of retroviruses by isolating nine Spuma endogenous retroviruses (ERVs) and one Epsilon ERV from the genomes of Agnatha and Chondrichthyes. Phylogenetic analysis of protein sequences revealed a striking pattern of co-evolution between jawless fish ERV and their host, while shark ERV underwent ancient cross-class viral transmission with jawless fish, ray-finned fish, and amphibians. Nucleotide sequence analysis showed that jawless fish ERV emerged in the Palaeozoic period, relatively later than ray-finned fish ERV. Moreover, codon analysis suggested that the jawless fish ERV employed an infection strategy that mimics the host codon. The genealogical diversity of ERVs in the jawless fish genome highlights the importance of studying different viral species. Overall, our findings provide valuable insights into the evolution of retroviruses and their interactions with their hosts.

Plants acquired mitochondrial linear plasmids horizontally from fungi likely during the conquest of land

Mitochondrial linear plasmids have been sporadically reported in fungi and plants. Yet, much remains obscure about the diversity, distribution, and evolution of mitochondrial linear plasmids. Here, through phylogenomic analyses across 7,163 cellular organisms (including 991 plants), we find that mitochondrial linear plasmids are widely present in land plants and fungi. Phylogenetic analyses indicate that plants are likely to have acquired mitochondrial linear plasmids horizontally from fungi before or during the conquest of terrestrial environments by plants. Gene content analyses show that mitochondrial linear plasmids harbor a highly dynamic and promiscuous repertoire of genes. Our study refines the understanding of the origin and evolution of mitochondrial linear plasmids.

Jawless vertebrates do not escape retrovirus infection

On occasion, retroviruses infect the genome of germline cell, forming endogenous retroviruses (ERVs), which provide molecular fossils for studying the deep evolution of retroviruses. While ERVs have been extensively characterized in the genomes of jawed vertebrates, much remains contentious and unexplored about the diversity and evolution of ERVs within jawless vertebrates. Here, we report the discovery of a novel ERV lineage, designated EbuERVs, in the genome of a hagfish Eptatretus burgeri. Phylogenetic analyses show that EbuERVs pertain to epsilon-retroviruses and might have derived from cross-species transmission from jawed vertebrates. EbuERVs are estimated to have invaded in the hagfish genome at least tens of millions of years ago. Evolutionary dynamics analyses indicate that EbuERVs might have experienced one proliferation peak and have been not active in transposition anymore. However, some EbuERVs can transcribe in embryo and might serve as lncRNA. Overall, these findings expand the distribution of retroviruses from jawed vertebrates to jawless vertebrates.

A Long-Running Arms Race between APOBEC1 Genes and Retroviruses in Tetrapods

Activation-induced cytidine deaminase/apolipoprotein B mRNA editing catalytic polypeptide-like (AID/APOBEC) proteins are cytosine deaminases implicated in diverse biological functions. APOBEC1 (A1) proteins have long been thought to regulate lipid metabolism, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. Endogenous retroviruses (ERVs) document past retroviral infections and are ubiquitous within the vertebrate genomes. Here, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in ERVs, and interrogate the evolutionary arms race between A1 genes and ERVs across vertebrate species. We find that A1 genes are widely present in tetrapods, recurrently amplified and lost in certain lineages, suggesting that A1 genes might have originated during the early evolution of tetrapods. A1-mediated mutation footprints can be detected in ERVs across tetrapods. Moreover, A1 genes appear to have experienced episodic positive selection in many tetrapod lineages. Taken together, we propose that a long-running arms race between A1 genes and retroviruses might have persisted throughout the evolutionary course of tetrapods. IMPORTANCE APOBEC3 (A3) genes have been thought to function in defense against retroviruses, whereas the evolutionary significance of A1 proteins in antiviral defense remains largely obscure. In this study, we identify the A1 gene repertoire, characterize the A1-mediated mutation footprints in endogenous retroviruses (ERVs), and explore the evolutionary arms race between A1 genes and ERVs across vertebrate species. We found A1 proteins originated during the early evolution of tetrapods, and detected the footprints of A1-induced hypermutations in retroviral fossils. A1 genes appear to have experienced pervasive positive selection in tetrapods. Our study indicates a long-running arms race between A1 genes and retroviruses taking place throughout the evolutionary course of tetrapods.

Evolutionary analysis of endogenous intronic retroviruses in primates reveals an enrichment in transcription binding sites associated with key regulatory processes

Background

Endogenous retroviruses (ERVs) are the result of the integration of retroviruses into host DNA following germline infection. Endogenous retroviruses are made up of three main genes: gag, pol, and env, each of which encodes viral proteins that can be conserved or not. ERVs have been observed in a wide range of vertebrate genomes and their functions are associated with viral silencing and gene regulation.

Results

In this work, we studied the evolutionary history of endogenous retroviruses associated with five human genes (INPP5B, DET1, PSMA1, USH2A, and MACROD2), which are located within intron sections. To verify the retroviral origin of the candidates, several approaches were used to detect and locate ERV elements. Both orthologous and paralogous genes were identified by Ensembl and then analyzed for ERV presence using RetroTector. A phylogenetic tree was reconstructed to identify the minimum time point of ERV acquisition. From that search, we detected ERVs throughout the primate lineage and in some other groups. Also, we identified the minimum origin of the ERVs from the parvorder Catarrhini to the Homininae subfamily.

Conclusions

With the data collected, and by observing the transcription factors annotated inside ERVs, we propose that these elements play a relevant role in gene expression regulation and they probably possess important features for tumorigenesis control.

Advances on genetic and genomic studies of ALV resistance

Avian leukosis (AL) is a general term for a variety of neoplastic diseases in avian caused by avian leukosis virus (ALV). No vaccine or drug is currently available for the disease. Therefore, the disease can result in severe economic losses in poultry flocks. Increasing the resistance of poultry to ALV may be one effective strategy. In this review, we provide an overview of the roles of genes associated with ALV infection in the poultry genome, including endogenous retroviruses, virus receptors, interferon-stimulated genes, and other immune-related genes. Furthermore, some methods and techniques that can improve ALV resistance in poultry are discussed. The objectives are willing to provide some valuable references for disease resistance breeding in poultry.

Cross-species transmission of an ancient endogenous retrovirus and convergent co-option of its envelope gene in two mammalian orders

Endogenous retroviruses (ERVs) found in vertebrate genomes are remnants of retroviral invasions of their ancestral species. ERVs thus represent molecular fossil records of ancient retroviruses and provide a unique opportunity to study viral-host interactions, including cross-species transmissions, in deep time. While most ERVs contain the mutated remains of the original retrovirus, on rare occasions evolutionary selection pressures lead to the co-option/exaptation of ERV genes for a host function. Here, we report the identification of two ancient related non-orthologous ERV env genes, ARTenvV and CARenvV, that are preserved with large open reading frames (ORFs) in the mammalian orders Artiodactyla and Carnivora, respectively, but are not found in other mammals. These Env proteins lack a transmembrane motif, but phylogenetic analyses show strong sequence preservation and positive selection of the env surface ORF in their respective orders, and transcriptomic analyses show a broad tissue expression pattern for both ARTenvV and CARenvV, suggesting that these genes may be exapted for a host function. Multiple lines of evidence indicate that ARTenvV and CARenvV were derived from an ancient ancestral exogenous gamma-like retrovirus that was independently endogenized in two mammalian orders more than 60 million years ago, which roughly coincides with the K-Pg mass extinction event and subsequent mammalian diversification. Thus, these findings identify the oldest known retroviral cross-ordinal transmission of a gamma-like retrovirus with no known extant infectious counterpart in mammals, and the first discovery of the convergent co-option of an ERV gene derived from the same ancestral retrovirus in two different mammalian orders.

Origin and Deep Evolution of Human Endogenous Retroviruses in Pan-Primates

Human endogenous retroviruses (HERVs) are viral “fossils” in the human genome that originated from the ancient integration of exogenous retroviruses. Although HERVs have sporadically been reported in nonhuman primate genomes, their deep origination in pan-primates remains to be explored. Hence, based on the in silico genomic mining of full-length HERVs in 49 primates, we performed the largest systematic survey to date of the distribution, phylogeny, and functional predictions of HERVs. Most importantly, we obtained conclusive evidence of nonhuman origin for most contemporary HERVs. We found that various supergroups, including HERVW9, HUERSP, HSERVIII, HERVIPADP, HERVK, and HERVHF, were widely distributed in Strepsirrhini, Platyrrhini (New World monkeys) and Catarrhini (Old World monkeys and apes). We found that numerous HERVHFs are spread by vertical transmission within Catarrhini and one HERVHF was traced in 17 species, indicating its ancient nature. We also discovered that 164 HERVs were likely involved in genomic rearrangement and 107 HERVs were potentially coopted in the form of noncoding RNAs (ncRNAs) in humans. In summary, we provided comprehensive data on the deep origination of modern HERVs in pan-primates.

Evolution and Genetic Diversity of the Retroviral Envelope in Anamniotes

Retroviruses are widely distributed in all vertebrates, as are their endogenous forms, endogenous retroviruses (ERV), which serve as "fossil" evidence to trace the ancient origins and history of virus-host interactions over millions of years. The retroviral envelope (Env) plays a significant role in host range determination, but major information on their genetic diversification and evolution in anamniotes is lacking. Here, by incorporating multiple-round in silico similarity search and phylogenomic analysis, more than 30,000 copies of ERV lineages with gamma-type Env (GTE), covalently associated Env, were discovered by searching against all fish and amphibian genomes and transcriptomic assemblies, but no beta-type Env (BTE), noncovalently associated Env, was found. Furthermore, a nine-type classification system of anamniote GTE was proposed by combining phylogenetic and domain/motif analyses. The elastic genomic organization and overall phylogenetic incongruence between anamniotic Env and its neighboring polymerase (Pol) implied that early retroviral diversification in anamniotic vertebrates was facilitated by frequent recombination. At last, host cellular opioid growth factor receptor (OGFr) gene capturing by anamniotic ERVs with GTE was reported for the first time. Overall, our findings overturn traditional Pol genotyping and reveal a complex evolutionary history of anamniotic retroviruses inferred by Env evolution. IMPORTANCE Although the retroviral envelope (Env) protein in amniotes has been well studied, its evolutionary history in anamniotic vertebrates is ambiguous. By analyzing more than 30,000 copies of ERV lineages with gamma-type Env (GTE) in anamniotes, several important evolutionary features were identified. First, GTE was found to be widely distributed among different amphibians and fish. Second, nine types of GTE were discovered and defined, revealing their great genetic diversity. Third, the incongruence between the Env and Pol phylogenies suggested that frequent recombination shaped the early evolution of anamniote retroviruses. Fourth, an ancient horizontal gene transfer event was discovered from anamniotes to ERVs with GTE. These findings reveal a complex evolution pattern for retroviral Env in anamniotes.

A Missing Link between Retrotransposons and Retroviruses

The origin and deep evolution of retroviruses remain largely unclear. It has been proposed that retroviruses might have originated from a Ty3/Gypsy retrotransposon, but all known Ty3/Gypsy retrotransposons are only distantly related to retroviruses. Retroviruses and some plant Athila/Tat elements (within Ty3/Gypsy retrotransposons) independently evolved a dual RNase H domain and an env/env-like gene. Here, we reported the discovery of a novel lineage of retrotransposons, designated Odin retrotransposons, in the genomes of eight sea anemones (order Actinaria) within the Cnidaria phylum. Odin retrotransposons exhibited unique genome features, encoding a dual RNase H domain (like retroviruses) but no env gene (like most Ty3/Gypsy retrotransposons). Phylogenetic analyses based on reverse transcriptase showed that Odin retrotransposons formed a sister group to lokiretroviruses, and lokiretroviruses and Odin retrotransposons together were sister to canonical retroviruses. Moreover, phylogenetic analyses based on RNase H and integrase also supported the hypothesis that Odin retrotransposons were sisters to lokiretroviruses. Lokiretroviruses and canonical retroviruses did not form a monophyletic group, indicating that lokiretroviruses and canonical retroviruses might represent two distinct virus families. Taken together, the discovery of Odin retrotransposons narrowed down the evolutionary gaps between retrotransposons and canonical retroviruses and lokiretroviruses.

The diversity and evolution of retroviruses: perspectives from viral “fossils”

Retroviruses exclusively infect vertebrates, causing a variety of diseases. The replication of retroviruses requires reverse transcription and integration into host genomes. When infecting germline cells, retroviruses become inherited vertically, forming endogenous retroviruses (ERVs). ERVs document past viral infections, providing molecular fossils for studying the evolutionary history of retroviruses. In this review, we summarize the recent advances in understanding the diversity and evolution of retroviruses from the perspectives of viral fossils, and discuss the effects of ERVs on the evolution of host biology.

•

Recent advances in understanding the diversity and evolution of retroviruses.

•

Methods to analyze ERVs.

•

The effects of ERVs on the evolution of host biology.

Genome-Wide Characterization of Zebrafish Endogenous Retroviruses Reveals Unexpected Diversity in Genetic Organizations and Functional Potentials

Endogenous retroviruses (ERVs) occupy a substantial fraction of mammalian genomes. However, whether ERVs extensively exist in ancient vertebrates remains unexplored. Here, we performed a genome-wide characterization of ERVs in a zebrafish (Danio rerio) model. Approximately 3,315 ERV-like elements (DrERVs) were identified as Gypsy, Copia, Bel, and class I−III groups. DrERVs accounted for approximately 2.3% of zebrafish genome and were distributed in all 25 chromosomes, with a remarkable bias on chromosome 4. Gypsy and class I are the two most abundant groups with earlier insertion times. The vast majority of the DrERVs have varied structural defects. A total of 509 gag and 71 env genes with coding potentials were detected. The env-coding elements were well-characterized and classified into four subgroups. A ERV-E4.8.43-DanRer element shows high similarity with HERV9NC-int in humans and analogous sequences were detected in species spanning from fish to mammals. RNA-seq data showed that hundreds of DrERVs were expressed in embryos and tissues under physiological conditions, and most of them exhibited stage and tissue specificity. Additionally, 421 DrERVs showed strong responsiveness to virus infection. A unique group of DrERVs with immune-relevant genes, such as fga, ddx41, ftr35, igl1c3, and tbk1, instead of intrinsic viral genes were identified. These DrERVs are regulated by transcriptional factors binding at the long terminal repeats. This study provided a survey of the composition, phylogeny, and potential functions of ERVs in a fish model, which benefits the understanding of the evolutionary history of ERVs from fish to mammals.

IMPORTANCE Endogenous retroviruses (ERVs) are relics of past infection that constitute up to 8% of the human genome. Understanding the genetic evolution of the ERV family and the interplay of ERVs and encoded RNAs and proteins with host function has become a new frontier in biology. Fish, as the most primitive vertebrate host for retroviruses, is an indispensable integral part for such investigations. In the present study, we report the genome-wide characterization of ERVs in zebrafish, an attractive model organism of ancient vertebrates from multiple perspectives, including composition, genomic organization, chromosome distribution, classification, phylogeny, insertion time, characterization of gag and env genes, and expression profiles in embryos and tissues. The result helps uncover the evolutionarily conserved and fish-specific ERVs, as well as the immune-relevant ERVs in response to virus infection. This study demonstrates the previously unrecognized abundance, diversification, and extensive activity of ERVs at the early stage of ERV evolution.

Retroviruses of Bats: a Threat Waiting in the Wings?

Bats are infamous reservoirs of deadly human viruses. While retroviruses, such as the human immunodeficiency virus (HIV), are among the most significant of virus families that have jumped from animals into humans, whether bat retroviruses have the potential to infect and cause disease in humans remains unknown. Recent reports of retroviruses circulating in bat populations builds on two decades of research describing the fossil records of retroviral sequences in bat genomes and of viral metagenomes extracted from bat samples. The impact of the global COVID-19 pandemic demands that we pay closer attention to viruses hosted by bats and their potential as a zoonotic threat. Here we review current knowledge of bat retroviruses and explore the question of whether they represent a threat to humans.

Predicted Cellular Interactors of the Endogenous Retrovirus-K Integrase Enzyme

Integrase (IN) enzymes are found in all retroviruses and are crucial in the retroviral integration process. Many studies have revealed how exogenous IN enzymes, such as the human immunodeficiency virus (HIV) IN, contribute to altered cellular function. However, the same consideration has not been given to viral IN originating from symbionts within our own DNA. Endogenous retrovirus-K (ERVK) is pathologically associated with neurological and inflammatory diseases along with several cancers. The ERVK IN interactome is unknown, and the question of how conserved the ERVK IN protein-protein interaction motifs are as compared to other retroviral integrases is addressed in this paper. The ERVK IN protein sequence was analyzed using the Eukaryotic Linear Motif (ELM) database, and the results are compared to ELMs of other betaretroviral INs and similar eukaryotic INs. A list of putative ERVK IN cellular protein interactors was curated from the ELM list and submitted for STRING analysis to generate an ERVK IN interactome. KEGG analysis was used to identify key pathways potentially influenced by ERVK IN. It was determined that the ERVK IN potentially interacts with cellular proteins involved in the DNA damage response (DDR), cell cycle, immunity, inflammation, cell signaling, selective autophagy, and intracellular trafficking. The most prominent pathway identified was viral carcinogenesis, in addition to select cancers, neurological diseases, and diabetic complications. This potentiates the role of ERVK IN in these pathologies via protein-protein interactions facilitating alterations in key disease pathways.

Molecular fossils illuminate the evolution of retroviruses following a macroevolutionary transition from land to water

The ancestor of cetaceans underwent a macroevolutionary transition from land to water early in the Eocene Period >50 million years ago. However, little is known about how diverse retroviruses evolved during this shift from terrestrial to aquatic environments. Did retroviruses transition into water accompanying their hosts? Did retroviruses infect cetaceans through cross-species transmission after cetaceans invaded the aquatic environments? Endogenous retroviruses (ERVs) provide important molecular fossils for tracing the evolution of retroviruses during this macroevolutionary transition. Here, we use a phylogenomic approach to study the origin and evolution of ERVs in cetaceans. We identify a total of 8,724 ERVs within the genomes of 25 cetaceans, and phylogenetic analyses suggest these ERVs cluster into 315 independent lineages, each of which represents one or more independent endogenization events. We find that cetacean ERVs originated through two possible routes. 298 ERV lineages may derive from retrovirus endogenization that occurred before or during the transition from land to water of cetaceans, and most of these cetacean ERVs were reaching evolutionary dead-ends. 17 ERV lineages are likely to arise from independent retrovirus endogenization events that occurred after the split of mysticetes and odontocetes, indicating that diverse retroviruses infected cetaceans through cross-species transmission from non-cetacean mammals after the transition to aquatic life of cetaceans. Both integration time and synteny analyses support the recent or ongoing activity of multiple retroviral lineages in cetaceans, some of which proliferated into hundreds of copies within the host genomes. Although ERVs only recorded a proportion of past retroviral infections, our findings illuminate the complex evolution of retroviruses during one of the most marked macroevolutionary transitions in vertebrate history.

Multiple Infiltration and Cross-Species Transmission of Foamy Viruses across the Paleozoic to the Cenozoic Era

Foamy viruses (FVs) are complex retroviruses that can infect humans and other animals. In this study, by integrating transcriptomic and genomic data, we discovered 412 FVs from 6 lineages in amphibians, which significantly increased the known set of FVs in amphibians. Among these lineages, salamander FVs maintained a coevolutionary pattern with their hosts that could be dated back to the Paleozoic era, while in contrast, frog FVs were much more likely acquired from cross-species (class-level) transmission in the Cenozoic era. In addition, we found that three distinct FV lineages had integrated into the genome of a salamander. Unexpectedly, we identified a lineage of endogenous FVs in caecilians that expressed all complete major genes, demonstrating the potential existence of an exogenous form of FV outside of mammals. Our discovery of rare phenomena in amphibian FVs has significantly increased our understanding of the macroevolution of the complex retrovirus. IMPORTANCE Foamy viruses (FVs) represent, more so than other viruses, the best model of coevolution between a virus and a host. This study represents the largest investigation so far of amphibian FVs and reveals 412 FVs of 6 distinct lineages from three major orders of amphibians. Besides a coevolutionary pattern, cross-species and repeated infections were also observed during the evolution of amphibian FVs. Remarkably, expressed FVs including a potential exogenous form were discovered, suggesting that active FVs might be underestimated in nature. These findings revealed that the multiple origins and complex evolution of amphibian FVs started from the Paleozoic era.

A Sister Lineage of Sampled Retroviruses Corroborates the Complex Evolution of Retroviruses

The origin and deep history of retroviruses remain mysterious and contentious, largely because the diversity of retroviruses is incompletely understood. Here, we report the discovery of lokiretroviruses, a novel major lineage of retroviruses, within the genomes of a wide range of vertebrates (at least 137 species), including lampreys, ray-finned fishes, lobe-finned fishes, amphibians, and reptiles. Lokiretroviruses share a similar genome architecture with known retroviruses, but display some unique features. Interestingly, lokiretrovirus Env proteins share detectable similarity with fusion glycoproteins of viruses within the Mononegavirales order, blurring the boundary between retroviruses and negative sense single-stranded RNA viruses. Phylogenetic analyses based on reverse transcriptase demonstrate that lokiretroviruses are sister to all the retroviruses sampled to date, providing a crucial nexus for studying the deep history of retroviruses. Comparing congruence between host and virus phylogenies suggests lokiretroviruses mainly underwent cross-species transmission. Moreover, we find that retroviruses replaced their ribonuclease H and integrase domains multiple times during their evolutionary course, revealing the importance of domain shuffling in the evolution of retroviruses. Overall, our findings greatly expand our views of the diversity of retroviruses, and provide novel insights into the origin and complex evolutionary history of retroviruses.

Frequent Retroviral Gene Co-option during the Evolution of Vertebrates

Endogenous retroviruses are ubiquitous in the vertebrate genomes. On occasion, hosts recruited retroviral genes to mediate their own biological functions, a process formally known as co-option or exaptation. Much remains unknown about the extent of retroviral gene co-option in vertebrates, although more than ten retroviral gene co-option events have been documented. Here, we use a phylogenomic approach to analyze more than 700 vertebrate genomes to uncover retroviral gene co-option taking place during the evolution of vertebrates. We identify a total of 177 independent retroviral gene co-option events in vertebrates, a majority of which have not been reported previously. Among these retroviral gene co-option events, 93 and 84 involve gag and env genes, respectively. More than 78.0% (138 out of 177) of retroviral gene co-option occurred within mammals. The gag and env co-option events share a generally similar temporal pattern with less frequent retroviral gene co-option identified in the deep branches, suggesting that retroviral gene co-option might have not been maintained for very long time periods. Moreover, we find co-opted retroviral genes are subject to different selection pressure, implying potentially diverse cellular functionality. Our study provides a comprehensive picture of co-opted retroviral genes during the evolution of vertebrates and has implications in understanding the ancient evolution of vertebrate-retrovirus interaction.

Unexpected Discovery and Expression of Amphibian Class II Endogenous Retroviruses

Endogenous retroviruses (ERVs) are the remnants of past retroviral infections. Fossil records of class II retroviruses have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Through genomic mining of all available amphibian genomes, we identified, for the first time, class II ERVs in amphibians. The class II ERVs were found only in Gymnophiona (caecilians) and not in the genomes of the other amphibian orders, Anura (frogs and toads) and Caudata (salamanders and newts), which are phylogenetically closely related. Therefore, the ERV endogenization occurred after the split of Gymnophiona, Anura, and Caudata (323 million years ago). Investigation of phylogenetic relationship and genomic structure revealed that the ERVs may originate from alpha- or betaretroviruses. We offer evidence that class II ERVs infiltrated amphibian genomes recently and may still have infectious members. Remarkably, certain amphibian class II ERVs can be expressed in diverse tissues. This discovery closes the major gap in the retroviral fossil record of class II ERVs and provides important insights into the evolution of class II ERVs in vertebrates.IMPORTANCE Class II retroviruses, largely distributed among mammals and birds, are of particular importance for medicine and economics. Class II ERVs have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Here, for the first time, we discovered class II ERVs in amphibians. We also revealed that the ERVs may originate from alpha- or betaretroviruses. We revealed that class II ERVs were integrated into amphibian genomes recently and certain amphibian class II ERVs can be expressed in diverse tissues. Our discovery closes the major gap in the retroviral fossil record of class II ERVs, and also indicates that amphibians may be still infected by class II retroviruses.

The Unique, the Known, and the Unknown of Spumaretrovirus Assembly

Within the family of Retroviridae, foamy viruses (FVs) are unique and unconventional with respect to many aspects in their molecular biology, including assembly and release of enveloped viral particles. Both components of the minimal assembly and release machinery, Gag and Env, display significant differences in their molecular structures and functions compared to the other retroviruses. This led to the placement of FVs into a separate subfamily, the Spumaretrovirinae. Here, we describe the molecular differences in FV Gag and Env, as well as Pol, which is translated as a separate protein and not in an orthoretroviral manner as a Gag-Pol fusion protein. This feature further complicates FV assembly since a specialized Pol encapsidation strategy via a tripartite Gag-genome–Pol complex is used. We try to relate the different features and specific interaction patterns of the FV Gag, Pol, and Env proteins in order to develop a comprehensive and dynamic picture of particle assembly and release, but also other features that are indirectly affected. Since FVs are at the root of the retrovirus tree, we aim at dissecting the unique/specialized features from those shared among the Spuma- and Orthoretrovirinae. Such analyses may shed light on the evolution and characteristics of virus envelopment since related viruses within the Ortervirales, for instance LTR retrotransposons, are characterized by different levels of envelopment, thus affecting the capacity for intercellular transmission.

Endogenous Retroviruses Drive Resistance and Promotion of Exogenous Retroviral Homologs

Endogenous retroviruses (ERVs) serve as markers of ancient viral infections and provide invaluable insight into host and viral evolution. ERVs have been exapted to assist in performing basic biological functions, including placentation, immune modulation, and oncogenesis. A subset of ERVs share high nucleotide similarity to circulating horizontally transmitted exogenous retrovirus (XRV) progenitors. In these cases, ERV-XRV interactions have been documented and include (a) recombination to result in ERV-XRV chimeras, (b) ERV induction of immune self-tolerance to XRV antigens, (c) ERV antigen interference with XRV receptor binding, and (d) interactions resulting in both enhancement and restriction of XRV infections. Whereas the mechanisms governing recombination and immune self-tolerance have been partially determined, enhancement and restriction of XRV infection are virus specific and only partially understood. This review summarizes interactions between six unique ERV-XRV pairs, highlighting important ERV biological functions and potential evolutionary histories in vertebrate hosts.

Presence of Endogenous Viral Elements Negatively Correlates with Feline Leukemia Virus Susceptibility in Puma and Domestic Cat Cells

While feline leukemia virus (FeLV) has been shown to infect felid species other than the endemic domestic cat host, differences in FeLV susceptibility among species has not been evaluated. Previous reports have noted a negative correlation between endogenous FeLV (enFeLV) copy number and exogenous FeLV (exFeLV) infection outcomes in domestic cats. Since felids outside the genus Felis do not harbor enFeLV genomes, we hypothesized absence of enFeLV results in more severe disease consequences in felid species lacking these genomic elements. We infected primary fibroblasts isolated from domestic cats (Felis catus) and pumas (Puma concolor) with FeLV and quantitated proviral and viral antigen loads. Domestic cat enFeLV env and long terminal repeat (LTR) copy numbers were determined for each individual and compared to FeLV viral outcomes. FeLV proviral and antigen levels were also measured in 6 naturally infected domestic cats and 11 naturally infected Florida panthers (P. concolor coryi). We demonstrated that puma fibroblasts are more permissive to FeLV than domestic cat cells, and domestic cat FeLV restriction was highly related to enFeLV-LTR copy number. Terminal tissues from FeLV-infected Florida panthers and domestic cats had similar exFeLV proviral copy numbers, but Florida panther tissues have higher FeLV antigen loads. Our work indicates that enFeLV-LTR elements negatively correlate with exogenous FeLV replication. Further, Puma concolor samples lacking enFeLV are more permissive to FeLV infection than domestic cat samples, suggesting that endogenization can play a beneficial role in mitigating exogenous retroviral infections. Conversely, presence of endogenous retroelements may relate to new host susceptibility during viral spillover events.IMPORTANCE Feline leukemia virus (FeLV) can infect a variety of felid species. Only the primary domestic cat host and related small cat species harbor a related endogenous virus in their genomes. Previous studies noted a negative association between the endogenous virus copy number and exogenous virus infection in domestic cats. This report shows that puma cells, which lack endogenous FeLV, produce more virus more rapidly than domestic cat fibroblasts following cell culture challenge. We document a strong association between domestic cat cell susceptibility and FeLV long terminal repeat (LTR) copy number, similar to observations in natural FeLV infections. Viral replication does not, however, correlate with FeLV env copy number, suggesting that this effect is specific to FeLV-LTR elements. This discovery indicates a protective capacity of the endogenous virus against the exogenous form, either via direct interference or indirectly via gene regulation, and may suggest evolutionary outcomes of retroviral endogenization.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Six reference-quality genomes reveal evolution of bat adaptations

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

The characterization of two novel neotropical primate papillomaviruses supports the ancient within-species diversity model

Papillomaviruses (PVs) are non-enveloped icosahedral viruses with a circular double-stranded DNA genome of ∼8,000 base pairs (bp). More than 200 different PV types have been identified to date in humans, which are distributed in five genera, with several strains associated with cancer development. Although widely distributed in vertebrates, Neotropical Primates (NP) PV infection was described for the first time only in 2016. Currently, four complete genomes of NP PVs have been characterized, three from Saimiri sciureus (SscPV1 to SscPV3) and one from Alouatta guariba (AgPV1). In this work, we describe two novel PV strains infecting Callithrix penicillata (provisionally named CpenPV1 and CpenPV2), using anal swab samples from animals residing at the Brasilia Primatology Center and next generation sequencing. The genomes of CpenPV1 (7,288 bp; 41.5% guanine-cytosine content - GC) and CpenPV2 (7,250 bp; 40.7% GC) contain the characteristic open reading frames (ORFs) for the early (E6, E7, E1, E2, and E4) and late (L2 and L1) PV genes. The L1 ORFs, commonly used for phylogenetic identification, share 76 per cent similarity with each other and differ 32 per cent from any other known PV, indicating that these new strains meet the criteria for defining novel species. PV genes phylogenetic variance was analyzed and different degrees of saturation revealed similar levels of topological heterogeneity, ruling out saturation as primary etiological factor for this phenomenon. Interestingly, the two CpenPV strains form a monophyletic clade within the Gammapapillomavirus genus (provisionally named gammapapillomavirus 32). Unlike for other NP PV strains, which grouped into a new sister genus of Alphapapillomavirus, this is the first report of NP PV strains grouping into a genus previously considered to exclusively comprise Old World Primates (OWP) PVs, including human PVs. These findings confirm the existence of a common ancestor for Gammapapillomavirus already infecting primates before the split of OWP and NP at ∼40 million years ago. Finally, our findings are consistent with an ancient within-species diversity model and emphasize the importance of increasing sampling to help understanding the PV-primate codivergence dynamics and pathogenic potential.

Cross-Species Transmission of Deltacoronavirus and the Origin of Porcine Deltacoronavirus

Deltacoronavirus is the last identified Coronaviridae subfamily genus. Differing from other coronavirus (CoV) genera, which mainly infect birds or mammals, deltacoronaviruses (δ‐CoVs) reportedly infect both animal types. Recent studies show that a novel δ‐CoV, porcine deltacoronavirus (PDCoV), can also infect calves and chickens with the potential to infect humans, raising the possibility of cross‐species transmission of δ‐CoVs. Here, we explored the deep phylogenetic history and cross‐species transmission of δ‐CoVs. Virus–host cophylogenetic analyses showed that δ‐CoVs have undergone frequent host switches in birds, and sparrows may serve as the unappreciated hubs for avian to mammal transmission. Our molecular clock analyses show that PDCoV possibly originated in Southeast Asia in the 1990s and that the PDCoV cluster shares a common ancestor with Sparrow‐CoV of around 1,810. Our findings contribute valuable insights into the diversification, evolution, and interspecies transmission of δ‐CoVs and the origin of PDCoV, providing a model for exploring the relationships of δ‐CoVs in birds and mammals.

From Ancient to Emerging Infections: The Odyssey of Viruses in the Male Genital Tract

The male genital tract (MGT) is the target of a number of viral infections that can have deleterious consequences at the individual, offspring, and population levels. These consequences include infertility, cancers of male organs, transmission to the embryo/fetal development abnormalities, and sexual dissemination of major viral pathogens such as human immunodeficiency virus (HIV) and hepatitis B virus. Lately, two emerging viruses, Zika and Ebola, have additionally revealed that the human MGT can constitute a reservoir for viruses cleared from peripheral circulation by the immune system, leading to their sexual transmission by cured men. This represents a concern for future epidemics and further underlines the need for a better understanding of the interplay between viruses and the MGT. We review here how viruses, from ancient viruses that integrated the germline during evolution through old viruses (e.g., papillomaviruses originating from Neanderthals) and more modern sexually transmitted infections (e.g., simian zoonotic HIV) to emerging viruses (e.g., Ebola and Zika) take advantage of genital tract colonization for horizontal dissemination, viral persistence, vertical transmission, and endogenization. The MGT immune responses to viruses and the impact of these infections are discussed. We summarize the latest data regarding the sources of viruses in semen and the complex role of this body fluid in sexual transmission. Finally, we introduce key animal findings that are relevant for our understanding of viral infection and persistence in the human MGT and suggest future research directions.

Frequent cross-species transmissions of foamy virus between domestic and wild felids

Emerging viral outbreaks resulting from host switching is an area of continued scientific interest. Such events can result in disease epidemics or in some cases, clinically silent outcomes. These occurrences are likely relatively common and can serve as tools to better understand disease dynamics, and may result in changes in behavior, fecundity, and, ultimately survival of the host. Feline foamy virus (FFV) is a common retrovirus infecting domestic cats globally, which has also been documented in the North American puma (Puma concolor). The prevalent nature of FFV in domestic cats and its ability to infect wild felids, including puma, provides an ideal system to study cross-species transmission across trophic levels (positions in the food chain), and evolution of pathogens transmitted between individuals following direct contact. Here we present findings from an extensive molecular analysis of FFV in pumas, focused on two locations in Colorado, and in relation to FFV recovered from domestic cats in this and previous studies. Prevalence of FFV in puma was high across the two regions, ∼77 per cent (urban interface site) and ∼48 per cent (rural site). Comparison of FFV from pumas living across three states; Colorado, Florida, and California, indicates FFV is widely distributed across North America. FFV isolated from domestic cats and pumas was not distinguishable at the host level, with FFV sequences sharing >93 per cent nucleotide similarity. Phylogenetic, Bayesian, and recombination analyses of FFV across the two species supports frequent cross-species spillover from domestic cat to puma during the last century, as well as frequent puma-to-puma intraspecific transmission in Colorado, USA. Two FFV variants, distinguished by significant difference in the surface unit of the envelope protein, were commonly found in both hosts. This trait is also shared by simian foamy virus and may represent variation in cell tropism or a unique immune evasion mechanism. This study elucidates evolutionary and cross-species transmission dynamics of a highly prevalent multi-host adapted virus, a system which can further be applied to model spillover and transmission of pathogenic viruses resulting in widespread infection in the new host.

Identification and evolution of avian endogenous foamy viruses

A history of long-term co-divergence means that foamy viruses (family Retroviridae) provide an ideal framework to understanding virus-host evolution over extended time periods. Endogenous foamy viruses (EndFVs) are rare, and to date have only been described in a limited number of mammals, amphibians, reptiles and fish genomes. By screening 414 avian genomes we identified EndFVs in two bird species: the Maguari Stork (Ciconia maguari) and the Oriental Stork (Ciconia boyciana). Analyses of phylogenetic relationships, genome structures and flanking sequences revealed a single origin of EndFVs in Ciconia species. In addition, the marked incongruence between the virus and host phylogenies suggested that this integration event occurred independently in birds. In sum, by providing evidence that birds can be infected with foamy viruses, we fill the last major gap in the taxonomic distribution of foamy viruses and their animal hosts.

Convergent Co-option of the Retroviral gag Gene during the Early Evolution of Mammals

Endogenous retroviruses, records of past retroviral infections, are ubiquitous in vertebrate genomes. On occasion, vertebrate hosts have co-opted retroviral genes for their own biological functions. Here, we perform a phylogenomic survey of retroviral gag gene homologs within vertebrate genomes and identify two ancient co-opted retroviral gag genes, designated wucaishi1 (wcs1) and wucaishi2 (wcs2), in mammals. Conserved synteny and evolutionary analyses suggest that the wcs1 and wcs2 co-options occurred before the origin of modern placental mammals (∼100 million years ago) and before the origin of modern marsupials (∼80 million years ago), respectively. We found that the wcs genes were lost or pseudogenized multiple times during the evolutionary course of mammals. While the wcs1 gene is mainly subject to negative selection in placental mammals (except in Perissodactyla), the wcs2 gene underwent positive selection in marsupials. Moreover, analyses of transcriptome-sequencing (RNA-seq) data suggest that the wcs1 and the wcs2 genes are expressed in a wide range of tissues. The convergent wcs co-option in mammals implies the retroviral gag gene might have been repurposed more frequently than previously thought.IMPORTANCE Retroviruses occasionally can infect host germ lines, forming endogenous retroviruses. Vertebrates, in turn, recruited retroviral genes for their own biological functions, a process formally known as co-option or exaptation. To date, co-opted retroviral gag genes have rarely been reported. In this study, we identified two co-opted retroviral gag genes, designated wucaishi1 (wcs1) and wucaishi2 (wcs2), in mammals. The co-option of wcs1 and wcs2 occurred before the origin of modern placentals and before the origin of modern marsupials, respectively. Our study indicates that retroviral gag gene co-option might have occurred more frequently than previously thought during the evolutionary course of vertebrates.

A reptilian endogenous foamy virus sheds light on the early evolution of retroviruses

Endogenous retroviruses (ERVs) represent host genomic 'fossils' of ancient viruses. Foamy viruses, including those that form endogenous copies, provide strong evidence for virus-host co-divergence across the vertebrate phylogeny. Endogenous foamy viruses (EFVs) have previously been discovered in mammals, amphibians, and fish. Here we report a novel endogenous foamy virus, termed ERV-Spuma-Spu, in genome of the tuatara (Sphenodon punctatus), an endangered reptile species endemic to New Zealand. Phylogenetic analyses revealed that foamy viruses have likely co-diverged with their hosts over many millions of years. The discovery of ERV-Spuma-Spu fills a major gap in the fossil record of foamy viruses and provides important insights into the early evolution of retroviruses.

Nomenclature for endogenous retrovirus (ERV) loci

Retroviral integration into germline DNA can result in the formation of a vertically inherited proviral sequence called an endogenous retrovirus (ERV). Over the course of their evolution, vertebrate genomes have accumulated many thousands of ERV loci. These sequences provide useful retrospective information about ancient retroviruses, and have also played an important role in shaping the evolution of vertebrate genomes. There is an immediate need for a unified system of nomenclature for ERV loci, not only to assist genome annotation, but also to facilitate research on ERVs and their impact on genome biology and evolution. In this review, we examine how ERV nomenclatures have developed, and consider the possibilities for the implementation of a systematic approach for naming ERV loci. We propose that such a nomenclature should not only provide unique identifiers for individual loci, but also denote orthologous relationships between ERVs in different species. In addition, we propose that-where possible-mnemonic links to previous, well-established names for ERV loci and groups should be retained. We show how this approach can be applied and integrated into existing taxonomic and nomenclature schemes for retroviruses, ERVs and transposable elements.