Foamy-like endogenous retroviruses are extensive and abundant in teleosts

Discovery of two novel foamy viruses in sea lions and dolphins provides insight into their evolutionary history

The prevalence and evolution of foamy viruses (FVs) have become the focus of research because of the risk of new zoonotic diseases. FVs have been isolated from various mammals and exhibit long-term co-speciation with their hosts. They also appear to be mild and nonpathogenic to their hosts. However, they may increase the risk of infection by other pathogens or exacerbate the symptoms of other diseases. Based on the data obtained using next-generation sequencing (NGS), we amplified and obtained the complete genomes of the two new FVs discovered in the bottlenose dolphin (Tursiops truncatus) and the South American sea lion (Otaria byronia) at the Qingdao Polar Haichang Ocean Park. Analysis and prediction of the novel FV's genomic structure revealed that it was consistent with that of the known mammalian FVs. The polmerase (pol) genes of the novel OFVoby_1 and DFVttr_1 showed less than 61.87 % and 61.83 % amino acid identity, respectively, with other known FVs belonging to the Retroviridae family. The host was likely to carry the FV for a considerable amount of time, as evidenced by the different times DFVttr_1 was discovered. The phylogenetic analysis revealed that the pol of OFVoby_1 and DFVttr_1 closely clustered with the FVs of Simiispumavirus and Felispumavirus, respectively. However, they both displayed distinct branches. According to the international committee on taxonomy of viruses (ICTV) FV classification criteria, FVs carried by dolphins and sea lions belong to two new genera within the Spumaretrovirinae subfamily. Using Bayesian analysis to simultaneously determine divergence dates and phylogenetic relationships revealed unique FVs with a divergence date of approximately 60 million years. This study helps us understand the FVs evolution and provides a scientific basis for future investigations into animal-borne infectious diseases.

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.

Discovery of an endogenous foamy virus in primitive ruminant chevrotains

Foamy viruses (FVs) are a class of complex retroviruses that could lead to persistent infections in specific species. FVs have two evolutionary characteristics: the first is co-evolution with the host, and the second is difficulty going through the endogenization process. To date, in mammals, only a few species (aye-aye, sloths, and cape golden mole) have been found to have endogenous forms of FVs. In this study, we identified a novel endogenous FV named TraEFV in the genomes of two species of even-toed ungulate chevrotain (genus Tragulus, family Tragulidae) known as mouse-deer. Phylogenetic analysis clustered TraEFVs with an exogenous FV isolated from feline but not with those isolated from cow and horse, and such inconsistent virus-host relationships reflect their complex evolutionary history. Interestingly, TraEFVs could be divided into two lineages, suggesting that TraEFV invaded these hosts at least twice during ancient times. Finally, the molecular clock estimates that TraEFV is approximately 20 million years old, suggesting its ancient nature. Our findings enrich the host taxonomy of spumaretroviruses. IMPORTANCE Foamy viruses (FV) are complex retroviruses that generally codiverge with their hosts. We identified a novel endogenous FV in the genomes of two mouse-deer species, the first endogenous FV found in Artiodactyla. The phylogenetic inconsistency of viruses and hosts suggested that the viruses might have emerged from cross-species transmission in the past. These findings indicate that ancient FVs may have had a wider range of hosts that remain to be expanded.

Identification of Cartilaginous Fish Endogenous Foamy Virus Rooting to Vertebrate Counterparts

Foamy viruses (FVs) are ideal models for studying the long-term evolutionary history between viruses and their hosts. Currently, FVs have been documented in nearly all major taxa of vertebrates, but evidence is lacking for true FV infiltration in cartilaginous fish, the most basal living vertebrates with jaws. Here, we screened 11 available genomes and 10 transcriptome sequence assemblies of cartilaginous fish and revealed a novel endogenous foamy virus, termed cartilaginous fish endogenous foamy virus (CFEFV), in the genomes of sharks and rays. Genomic analysis of CFEFVs revealed feature motifs that were retained among canonical FVs. Phylogenetic analysis using polymerase sequences revealed the rooting nature of CFEFVs to vertebrate FVs, indicating their deep origin. Interestingly, three viral lineages were found in a shark (Scyliorhinus torazame), one of which was clustered with ray-finned fish foamy-like viruses, indicating that multiple episodes of viral infiltrations had occurred in this species. These findings fill a major gap in the Spumaretrovirinae taxon and reveal the aquatic origin of FVs found in terrestrial vertebrates. IMPORTANCE Although foamy viruses (FVs) have been found in major branches of vertebrates, the presence of these viruses in cartilaginous fish, the most basal living vertebrates with jaws, remains to be explored. This study revealed a collection of cartilaginous endogenous FVs in sharks and rays through in silico genomic mining. These viruses were rooted in the polymerase (POL) phylogeny, indicating the ancient aquatic origin of FVs. However, their envelope (ENV) protein grouped with those of amphibian FVs, suggesting different evolutionary histories of different FV genes. Overall, we provide the last missing gap for the taxonomic investigation of Spumaretrovirinae and provide concrete support for the aquatic origin of FVs.

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.

A novel class III endogenous retrovirus with a class I envelope gene in African frogs with an intact genome and developmentally regulated transcripts in Xenopus tropicalis

BACKGROUND

Retroviruses exist as exogenous infectious agents and as endogenous retroviruses (ERVs) integrated into host chromosomes. Such endogenous retroviruses (ERVs) are grouped into three classes roughly corresponding to the seven genera of infectious retroviruses: class I (gamma-, epsilonretroviruses), class II (alpha-, beta-, delta-, lentiretroviruses) and class III (spumaretroviruses). Some ERVs have counterparts among the known infectious retroviruses, while others represent paleovirological relics of extinct or undiscovered retroviruses.

RESULTS

Here we identify an intact ERV in the Anuran amphibian, Xenopus tropicalis. XtERV-S has open reading frames (ORFs) for gag, pol (polymerase) and env (envelope) genes, with a small additional ORF in pol and a serine tRNA primer binding site. It has unusual features and domain relationships to known retroviruses. Analyses based on phylogeny and functional motifs establish that XtERV-S gag and pol genes are related to the ancient env-less class III ERV-L family but the surface subunit of env is unrelated to known retroviruses while its transmembrane subunit is class I-like. LTR constructs show transcriptional activity, and XtERV-S transcripts are detected in embryos after the maternal to zygotic mid-blastula transition and before the late tailbud stage. Tagged Gag protein shows typical subcellular localization. The presence of ORFs in all three protein-coding regions along with identical 5' and 3' LTRs (long terminal repeats) indicate this is a very recent germline acquisition. There are older, full-length, nonorthologous, defective copies in Xenopus laevis and the distantly related African bullfrog, Pyxicephalus adspersus. Additional older, internally deleted copies in X. tropicalis carry a 300 bp LTR substitution.

CONCLUSIONS

XtERV-S represents a genera-spanning member of the largely env-less class III ERV that has ancient and modern copies in Anurans. This provirus has an env ORF with a surface subunit unrelated to known retroviruses and a transmembrane subunit related to class I gammaretroviruses in sequence and organization, and is expressed in early embryogenesis. Additional XtERV-S-related but defective copies are present in X. tropicalis and other African frog taxa. XtERV-S is an unusual class III ERV variant, and it may represent an important transitional retroviral form that has been spreading in African frogs for tens of millions of years.

Transmission of infectious viruses in the natural setting at human-animal interface

Most viral pathogens causing epidemics and pandemics are zoonotic, emerging from wildlife reservoirs like SARS CoV2 causing the global Covid-19 pandemic, although animal origin of this virus remains a mystery. Cross-species transmission of pathogens from animals to humans is known as zoonosis. However, pathogens are also transmitted from humans to animals in regions where there is a close interaction between animals and humans by 'reverse transmission' (anthroponosis). Molecular evidence for the transmission of two zoonotic RNA viruses at the human-monkey interface in Rajasthan forests is presented here: a) the apathogenic Simian Foamy Viruses (SFV), and b): Influenza A viruses (IAV)-like virus, etiologic agent for human flu infecting wild Indian rhesus monkeys inhabiting Rajasthan forests. The data provide critical information on ecology and evolution of viruses of Public Health relevance. During replication, viral genomes mutate along the transmission route to adapt to the new hosts, generating new variants that are likely to have properties different from the founder viruses. Wild Indian monkeys are under-sampled for monitoring infectious diseases mainly because of the difficulties with sample collection. Monkeys are perceived as religious icons by the Hindus in India. It is extremely difficult to obtain permission from the Forest and Wildlife Department government authorities to collect wild simian blood samples for surveillance of infectious diseases caused by viral pathogens. Reducing animal-human contact and affordable vaccination are two relevant anti-viral strategies to counteract the spread of infectious zoonotic pathogens. Genbank Accession numbers: Indian SFVmac: ADN94420, IAV like virus: MZ298601.

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.

Genomic characterization and distribution of bovine foamy virus in Japan

Bovine foamy virus (BFV) is distributed through worldwide cattle herds. Although the biological features of BFV are not well understood, appearance of clinical manifestation by superinfection with other microorganisms is inferred. In Japan, reports of genomic characterizations and epidemiology of this virus are limited. In this study, we performed whole genomic sequencing of BFV strains Ibaraki and No.43, which were isolated in this country. Additionally, we investigated BFV in geographically distant four daily farms in Japan, to estimate the distribution of BFV and its correlation to bovine leukemia virus (BLV). BFV was distributed throughout Japan; the average positive rate was 12.7%. The nucleotide sequence identities of the isolates were 99.6% when compared with BFV strain isolated in the USA. The phylogenetic tree using env gene sequence showed strains Ibaraki, No.43 and Kagoshima were sorted in the same cluster including the USA and Chinese strains, while Hokkaido strain was in the other cluster including European strains. Although no clear correlation between BFV and BLV could be found, BFV and BLV infections were likely to increase with ages. Our data on epidemiology and characteristics of BFV will provide important information to reveal biological features of BFV.

The Genetic Basis of Morphological Diversity in Domesticated Goldfish

Although domesticated goldfish strains exhibit highly diversified phenotypes in morphology, the genetic basis underlying these phenotypes is poorly understood. Here, based on analysis of transposable elements in the allotetraploid goldfish genome, we found that its two subgenomes have evolved asymmetrically since a whole-genome duplication event in the ancestor of goldfish and common carp. We conducted whole-genome sequencing of 27 domesticated goldfish strains and wild goldfish. We identified more than 60 million genetic variations and established a population genetic structure of major goldfish strains. Genome-wide association studies and analysis of strain-specific variants revealed genetic loci associated with several goldfish phenotypes, including dorsal fin loss, long-tail, telescope-eye, albinism, and heart-shaped tail. Our results suggest that accumulated mutations in the asymmetrically evolved subgenomes led to generation of diverse phenotypes in the goldfish domestication history. This study is a key resource for understanding the genetic basis of phenotypic diversity among goldfish strains.

Avian and serpentine endogenous foamy viruses, and new insights into the macroevolutionary history of foamy viruses

This study reports and characterises two novel distinct lineages of foamy viruses (FVs) in the forms of endogenous retroviruses (ERVs). Several closely related elements were found in the genome of oriental stork (Ciconia boyciana) and other was found in the genome of spine-bellied sea snake (Hydrophis hardwickii), designated ERV-Spuma.N-Cbo (where 'N' runs from one to thirteen) and ERV-Spuma.1-Hha, respectively. This discovery of avian and serpentine endogenous FVs adds snakes, and perhaps more crucially, birds to the list of currently known hosts of FVs, in addition to mammals, reptiles, amphibians, and fish. This indicates that FVs are, or at least were, capable of infecting all major lineages of vertebrates. Moreover, together with other FVs, phylogenetic analyses showed that both of them are most closely related to mammalian FVs. Further examination revealed that reptilian FVs form a deep paraphyletic group that is basal to mammalian and avian FVs, suggesting that there were multiple ancient FV cross-class transmissions among their hosts. Evolutionary timescales of various FV lineages were estimated in this study, in particular, the timescales of reptilian FVs and that of the clade of mammalian, avian, and serpentine FVs. This was accomplished by using the recently established time-dependent rate phenomenon models, inferred using mainly the knowledge of the co-speciation history between FVs and mammals. It was found that the estimated timescales matched very well with those of reptiles. Combined with the observed phylogenetic patterns, these results suggested that FVs likely co-speciated with ancient reptilian animals, but later jumped to a protomammal and/or a bird, which ultimately gave rise to mammalian and avian FVs. These results contribute to our understanding of FV emergence, specifically the emergence of mammalian and avian FVs, and provide new insights into how FVs co-evolved with their non-mammalian vertebrate hosts in the distant past.

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.

Simian Foamy Viruses in Central and South America: A New World of Discovery

Foamy viruses (FVs) are the only exogenous retrovirus to date known to infect neotropical primates (NPs). In the last decade, an increasing number of strains have been completely or partially sequenced, and molecular evolution analyses have identified an ancient co-speciation with their hosts. In this review, the improvement of diagnostic techniques that allowed the determination of a more accurate prevalence of simian FVs (SFVs) in captive and free-living NPs is discussed. Determination of DNA viral load in American primates indicates that oral tissues are the viral replicative site and that buccal swab collection can be an alternative to diagnose SFV infection in NPs. Finally, the transmission potential of NP SFVs to primate workers in zoos and primate centers of the Americas is examined.

Modular nature of simian foamy virus genomes and their evolutionary history

Among all known retroviruses, foamy viruses (FVs) have the most stable virus–host co-speciation history, co-diverging in concert with their vertebrate hosts for hundreds of millions of years. However, detailed molecular analyses indicate that different parts of their genome might have different evolutionary histories. While their polymerase gene displays a robust and straightforward virus–host co-speciation pattern, the evolutionary history of their envelope (env) gene, is much more complicated. Here, we report eleven new FV env sequences in two mandrill populations in Central Africa, geographically separated by the Ogooué River into the North and the South populations. Phylogenetic reconstruction of the polymerase gene shows that the two virus populations are distinct, and each contains two variants of env genes co-existing with one another. The distinction between the two env variants can be mapped to the surface domain, flanked by two recombination hotspots, as previously reported for chimpanzee and gorilla FVs. Our analyses suggest that the two env variants originated during the diversification of Old World monkeys and apes, ∼30 million years ago. We also show that this env gene region forms two phylogenetically distinct clades, each displaying a host co-divergence and geographical separation pattern, while the rest of the genome of the two strains is phylogenetically indistinguishable in each of the host-specific groups. We propose possible evolutionary mechanisms to explain the modular nature of the FV genome.

The First Co-Opted Endogenous Foamy Viruses and the Evolutionary History of Reptilian Foamy Viruses

A recent study reported the discovery of an endogenous reptilian foamy virus (FV), termed ERV-Spuma-Spu, found in the genome of tuatara. Here, we report two novel reptilian foamy viruses also identified as endogenous FVs (EFVs) in the genomes of panther gecko (ERV-Spuma-Ppi) and Schlegel’s Japanese gecko (ERV-Spuma-Gja). Their presence indicates that FVs are capable of infecting reptiles in addition to mammals, amphibians, and fish. Numerous copies of full length ERV-Spuma-Spu elements were found in the tuatara genome littered with in-frame stop codons and transposable elements, suggesting that they are indeed endogenous and are not functional. ERV-Spuma-Ppi and ERV-Spuma-Gja, on the other hand, consist solely of a foamy virus-like env gene. Examination of host flanking sequences revealed that they are orthologous, and despite being more than 96 million years old, their env reading frames are fully coding competent with evidence for strong purifying selection to maintain expression and for them likely being transcriptionally active. These make them the oldest EFVs discovered thus far and the first documented EFVs that may have been co-opted for potential cellular functions. Phylogenetic analyses revealed a complex virus–host co-evolutionary history and cross-species transmission routes of ancient FVs.

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

The deep history and early diversification of retroviruses remains elusive, largely because few retroviruses have been characterized in vertebrates other than mammals and birds. Endogenous retroviruses (ERVs) documented past retroviral infections and thus provide ‘molecular fossils’ for studying the deep history of retroviruses. Here we perform a comprehensive phylogenomic analysis of ERVs within the genomes of 92 non-avian/mammalian vertebrates, including 72 fishes, 4 amphibians, and 16 reptiles. We find that ERVs are present in all the genomes of jawed vertebrates, revealing the ubiquitous presence of ERVs in jawed vertebrates. We identify a total of >8,000 ERVs and reconstruct ~450 complete or partial ERV genomes, which dramatically expands the phylogenetic diversity of retroviruses and suggests that the diversity of exogenous retroviruses might be much underestimated in non-avian/mammalian vertebrates. Phylogenetic analyses show that retroviruses cluster into five major groups with different host distributions, providing important insights into the classification and diversification of retroviruses. Moreover, we find retroviruses mainly underwent frequent host switches in non-avian/mammalian vertebrates, with exception of spumavirus-related viruses that codiverged with their ray-finned fish hosts. Interestingly, ray-finned fishes and turtles appear to serve as unappreciated hubs for the transmission of retroviruses. Finally, we find retroviruses underwent many independent water-land transmissions, indicating the water-land interface is not a strict barrier for retrovirus transmission. Our analyses provide unprecedented insights into and valuable resources for studying the diversification, key evolutionary transitions, and macroevolution of retroviruses.

Retroviruses infect a wide range of vertebrates and cause many diseases, such as AIDS and cancers. To date, retroviruses have been rarely characterized in vertebrates other than mammals and birds, impeding our understanding of the diversity and early evolution of retroviruses. Retroviruses can occasionally integrate into host genomes and become endogenous retroviruses (ERVs), which provide molecular fossils for studying the long-term evolution of retroviruses. Here we performed comparative genomic and evolutionary analyses of ERVs within 92 non-avian/mammalian vertebrates (fishes, amphibians, and reptiles) and uncovered extraordinary diversity of retroviruses in non-avian/mammalian vertebrates. Our analyses reveal an ancient aquatic origin of retroviruses and retroviruses underwent frequent host-switching. Our findings have important implications in understanding the deep history and evolutionary mode of retroviruses.

Spumaretroviruses: Updated taxonomy and nomenclature

Spumaretroviruses, commonly referred to as foamy viruses, are complex retroviruses belonging to the subfamily Spumaretrovirinae, family Retroviridae, which naturally infect a variety of animals including nonhuman primates (NHPs). Additionally, cross-species transmissions of simian foamy viruses (SFVs) to humans have occurred following exposure to tissues of infected NHPs. Recent research has led to the identification of previously unknown exogenous foamy viruses, and to the discovery of endogenous spumaretrovirus sequences in a variety of host genomes. Here, we describe an updated spumaretrovirus taxonomy that has been recently accepted by the International Committee on Taxonomy of Viruses (ICTV) Executive Committee, and describe a virus nomenclature that is generally consistent with that used for other retroviruses, such as lentiviruses and deltaretroviruses. This taxonomy can be applied to distinguish different, but closely related, primate (e.g., human, ape, simian) foamy viruses as well as those from other hosts. This proposal accounts for host-virus co-speciation and cross-species transmission.

Unique genome organization of non-mammalian papillomaviruses provides insights into the evolution of viral early proteins

The family Papillomaviridae contains more than 320 papillomavirus types, with most having been identified as infecting skin and mucosal epithelium in mammalian hosts. To date, only nine non-mammalian papillomaviruses have been described from birds (n = 5), a fish (n = 1), a snake (n = 1), and turtles (n = 2). The identification of papillomaviruses in sauropsids and a sparid fish suggests that early ancestors of papillomaviruses were already infecting the earliest Euteleostomi. The Euteleostomi clade includes more than 90 per cent of the living vertebrate species, and progeny virus could have been passed on to all members of this clade, inhabiting virtually every habitat on the planet. As part of this study, we isolated a novel papillomavirus from a 16-year-old female Adélie penguin (Pygoscelis adeliae) from Cape Crozier, Ross Island (Antarctica). The new papillomavirus shares ∼64 per cent genome-wide identity to a previously described Adélie penguin papillomavirus. Phylogenetic analyses show that the non-mammalian viruses (expect the python, Morelia spilota, associated papillomavirus) cluster near the base of the papillomavirus evolutionary tree. A papillomavirus isolated from an avian host (Northern fulmar; Fulmarus glacialis), like the two turtle papillomaviruses, lacks a putative E9 protein that is found in all other avian papillomaviruses. Furthermore, the Northern fulmar papillomavirus has an E7 more similar to the mammalian viruses than the other avian papillomaviruses. Typical E6 proteins of mammalian papillomaviruses have two Zinc finger motifs, whereas the sauropsid papillomaviruses only have one such motif. Furthermore, this motif is absent in the fish papillomavirus. Thus, it is highly likely that the most recent common ancestor of the mammalian and sauropsid papillomaviruses had a single motif E6. It appears that a motif duplication resulted in mammalian papillomaviruses having a double Zinc finger motif in E6. We estimated the divergence time between Northern fulmar-associated papillomavirus and the other Sauropsid papillomaviruses be to around 250 million years ago, during the Paleozoic-Mesozoic transition and our analysis dates the root of the papillomavirus tree between 400 and 600 million years ago. Our analysis shows evidence for niche adaptation and that these non-mammalian viruses have highly divergent E6 and E7 proteins, providing insights into the evolution of the early viral (onco-)proteins.