Widespread Horizontal Gene Transfer Among Animal Viruses.

The initial objective of this study was to shed light on the evolution of small DNA tumor viruses by analyzing de novo assemblies of publicly available deep sequencing datasets. The survey generated a searchable database of contig snapshots representing more than 100,000 Sequence Read Archive records. Using modern structure-aware search tools, we iteratively broadened the search to include an increasingly wide range of other virus families. The analysis revealed a surprisingly diverse range of chimeras involving different virus groups. In some instances, genes resembling known DNA-replication modules or known virion protein operons were paired with unrecognizable sequences that structural predictions suggest may represent previously unknown replicases and novel virion architectures. Discrete clades of an emerging group called adintoviruses were discovered in datasets representing humans and other primates. As a proof of concept, we show that the contig database is also useful for discovering RNA viruses and candidate archaeal phages. The ancillary searches revealed additional examples of chimerization between different virus groups. The observations support a gene-centric taxonomic framework that should be useful for future virus-hunting efforts.

Two Lineages of Papillomaviruses Identified from Caracals (Caracal caracal) in South Africa.

Papillomaviruses (PV) infect epithelial cells and can cause hyperplastic or neoplastic lesions. In felids, most described PVs are from domestic cats (Felis catus; n = 7 types), with one type identified in each of the five wild felid species studied to date (Panthera uncia, Puma concolor, Leopardus wiedii, Panthera leo persica and Lynx rufus). PVs from domestic cats are highly diverse and are currently classified into three genera (Lambdapapillomavirus, Dyothetapapillomavirus, and Taupapillomavirus), whereas those from wild felids, although diverse, are all classified into the Lambdapapillomavirus genus. In this study, we used a metagenomic approach to identify ten novel PV genomes from rectal swabs of five deceased caracals (Caracal caracal) living in the greater Cape Town area, South Africa. These are the first PVs to be described from caracals, and represent six new PV types, i.e., Caracal caracal papillomavirus (CcarPV) 1-6. These CcarPV fall into two phylogenetically distinct genera: Lambdapapillomavirus, and Treisetapapillomavirus. Two or more PV types were identified in a single individual for three of the five caracals, and four caracals shared at least one of the same PV types with another caracal. This study broadens our understanding of wild felid PVs and provides evidence that there may be several wild felid PV lineages.

Two Novel Geminiviruses Identified in Bees (Apis mellifera and Nomia sp.).

Members of the Geminviridae family are circular single-stranded DNA plant-infecting viruses, some of which impact global food production. Geminiviruses are vectored by sap-feeding insects such as leafhoppers, treehoppers, aphids, and whiteflies. Additionally, geminivirus sequences have also been identified in other insects such as dragonflies, mosquitoes, and stingless bees. As part of a viral metagenomics study on honeybees and solitary bees (Nomia sp.), two geminivirus genomes were identified. These represent a novel citlodavirus (from honeybees collected from Westmoreland, Jamaica) and a mastrevirus-like genome (from a solitary bee collected from Tempe, Arizona, USA). The novel honeybee-derived citlodavirus genome shares ~61 to 69% genome-wide nucleotide pairwise identity with other citlodavirus genome sequences and is most closely related to the passion fruit chlorotic mottle virus identified in Brazil. Whereas the novel solitary bee-derived mastrevirus-like genome shares ~55 to 61% genome-wide nucleotide identity with other mastreviruses and is most closely related to tobacco yellow dwarf virus identified in Australia, based on pairwise identity scores of the full genome, replication-associated protein, and capsid protein sequences. Previously, two geminiviruses in the Begomovirus genus were identified in samples of stingless bee (Trigona spp.) samples. Here, we identify viruses that represent two new species of geminiviruses from a honeybee and a solitary bee, which continues to demonstrate that plant pollinators can be utilized for the identification of plant-infecting DNA viruses in ecosystems.

Diverse papillomaviruses identified from Antarctic fur seals, leopard seals and Weddell seals from the Antarctic.

Papillomaviruses (family Papillomaviridae) are non-enveloped, circular, double-stranded DNA viruses known to infect squamous and mucosal epithelial cells. In the family Papillomaviridae there are 53 genera and 133 viral species whose members infect a variety of mammalian, avian, reptilian, and fish species. Within the Antarctic context, papillomaviruses (PVs) have been identified in Adélie penguins (Pygoscelis adeliae, 2 PVs), Weddell seals (Leptonychotes weddellii, 7 PVs), and emerald notothen (Trematomus bernacchii, 1 PV) in McMurdo Sound and Ross Island in eastern Antarctica. Here we identified 13 diverse PVs from buccal swabs of Antarctic fur seals (Arctocephalus gazella, 2 PVs) and leopard seal (Hydrurga leptonyx, 3 PVs) in western Antarctica (Antarctic Peninsula), and vaginal and nasal swabs of Weddell seals (8 PVs) in McMurdo Sound. These PV genomes group into four genera representing 11 new papillomavirus types, of which five are from two Antarctic fur seals and a leopard seal and six from Weddell seals.

New duck papillomavirus type identified in a mallard in Missouri, USA.

Papillomaviruses are small circular DNA viruses that infect epithelial and mucosal cells and have co-evolved with their hosts. Some papillomaviruses in mammals are well studied (especially those associated with disease). However, there is limited information on papillomaviruses associated with avian hosts. From a cloacal swab sample of a mallard (Anas platyrhynchos) sampled in Missouri, USA (6 Jan 2023), we identified a papillomavirus (7839 nt) that shares ~68% genome-wide nucleotide sequence identity with Anas platyrhynchos papillomavirus 1 (AplaPV1) from a mallard sampled in Newfoundland (Canada) and ~40% with AplaPV2 from a mallard sampled in Minnesota (USA) with mesenchymal dermal tumors. The papillomavirus we identified shares 73.6% nucleotide sequence identity in the L1 gene with that of AplaPV1 and thus represents a new AplaPV type (AplaPV3). The genome sequence of AplaPV3 shares >97% identity with three partial PV genome sequences (1316, 1997, and 4241 nt) identified in a mallard in India, indicating that that virus was also AplaPV3.

Identification of a novel polyomavirus in wild Sonoran Desert rodents of the family Heteromyidae.

Rodents are the largest and most diverse group of mammals. Covering a wide range of structural and functional adaptations, rodents successfully occupy virtually every terrestrial habitat, and they are often found in close association with humans, domestic animals, and wildlife. Although a significant amount of research has focused on rodents' prominence as known reservoirs of zoonotic viruses, there has been less emphasis on the viral ecology of rodents in general. Here, we utilized a viral metagenomics approach to investigate polyomaviruses in wild rodents from the Baja California peninsula, Mexico, using fecal samples. We identified a novel polyomavirus in fecal samples from two rodent species, a spiny pocket mouse (Chaetodipus spinatus) and a Dulzura kangaroo rat (Dipodomys simulans). These two polyomaviruses represent a new species in the genus Betapolyomavirus. Sequences of this polyomavirus cluster phylogenetically with those of other rodent polyomaviruses and two other non-rodent polyomaviruses (WU and KI) that have been identified in the human respiratory tract. Through our continued work on seven species of rodents, we endeavor to explore the viral diversity associated with wild rodents on the Baja California peninsula and expand on current knowledge of rodent viral ecology and evolution.

Identification and characterization of a polyomavirus in the thornback skate (Raja clavata).

Members of the family Polyomaviridae have a circular double-stranded DNA genome that have been identified in various hosts ranging from mammals to arachnids. Here we report the identification and analysis of a complete genome sequence of a novel polyomavirus, Raja clavata polyomavirus (RcPyV1), from a cartilaginous fish, the thornback skate (Raja clavata). The genome sequence was determined using a metagenomics approach with an aim to provide baseline viral data in cartilaginous fish in different ecosystems. The RcPyV1 genome (4,195 nucleotides) had typical organization of polyomavirus, including early antigens (small T; Large T) encoded on one strand and late viral proteins (VP1; VP2) on the complementary strand. Maximum-likelihood phylogenetic analysis of the large T-antigen revealed that RcPyV1 clusters with a polyomavirus obtained from another cartilaginous fish, the guitarfish polyomavirus 1 (GfPyV1). These two share ~ 56% pairwise identity in LT and VP1 protein sequences. These analyses support the hypothesis that cartilaginous fishes have a specific lineage of polyomaviruses.

Viral anti-inflammatory serpin reduces immuno-coagulopathic pathology in SARS-CoV-2 mouse models of infection.

SARS-CoV-2 acute respiratory distress syndrome (ARDS) induces uncontrolled lung inflammation and coagulopathy with high mortality. Anti-viral drugs and monoclonal antibodies reduce early COVID-19 severity, but treatments for late-stage immuno-thrombotic syndromes and long COVID are limited. Serine protease inhibitors (SERPINS) regulate activated proteases. The myxoma virus-derived Serp-1 protein is a secreted immunomodulatory serpin that targets activated thrombotic, thrombolytic, and complement proteases as a self-defense strategy to combat clearance. Serp-1 is effective in multiple animal models of inflammatory lung disease and vasculitis. Here, we describe systemic treatment with purified PEGylated Serp-1 as a therapy for immuno-coagulopathic complications during ARDS. Treatment with PEGSerp-1 in two mouse-adapted SARS-CoV-2 models in C57Bl/6 and BALB/c mice reduced lung and heart inflammation, with improved outcomes. PEGSerp-1 significantly reduced M1 macrophages in the lung and heart by modifying urokinase-type plasminogen activator receptor (uPAR), thrombotic proteases, and complement membrane attack complex (MAC). Sequential changes in gene expression for uPAR and serpins (complement and plasminogen inhibitors) were observed. PEGSerp-1 is a highly effective immune-modulator with therapeutic potential for severe viral ARDS, immuno-coagulopathic responses, and Long COVID.

Cervicovaginal Human Papillomavirus Genomes, Microbiota Composition and Cytokine Concentrations in South African Adolescents.

The interaction between cervicovaginal virome, bacteriome and genital inflammation has not been extensively investigated. We assessed the vaginal DNA virome from 33 South African adolescents (15-19 years old) using shotgun DNA sequencing of purified virions. We present analyses of eukaryote-infecting DNA viruses, with a focus on human papillomavirus (HPV) genomes and relate these to the vaginal bacterial microbiota (assessed by 16S rRNA gene sequencing) and cytokines (assessed by Luminex). The DNA virome included single-stranded (Anelloviridae, Genomoviridae) and double-stranded DNA viruses (Adenoviridae, Alloherpesviridae, Herpesviridae, Marseilleviridae, Mimiviridae, Polyomaviridae, Poxviridae). We identified 110 unique, complete HPV genomes within two genera (Alphapapillomavirus and Gammapapillomavirus) representing 40 HPV types and 12 species. Of the 40 HPV types identified, 35 showed positive co-infection patterns with at least one other type, mainly HPV-16. HPV-35, a high-risk genotype currently not targeted by available vaccines, was the most prevalent HPV type identified in this cohort. Bacterial taxa commonly associated with bacterial vaginosis also correlated with the presence of HPV. Bacterial vaginosis, rather than HPV, was associated with increased genital inflammation. This study lays the foundation for future work characterizing the vaginal virome and its role in women's health.

Four principles to establish a universal virus taxonomy.

A universal taxonomy of viruses is essential for a comprehensive view of the virus world and for communicating the complicated evolutionary relationships among viruses. However, there are major differences in the conceptualisation and approaches to virus classification and nomenclature among virologists, clinicians, agronomists, and other interested parties. Here, we provide recommendations to guide the construction of a coherent and comprehensive virus taxonomy, based on expert scientific consensus. Firstly, assignments of viruses should be congruent with the best attainable reconstruction of their evolutionary histories, i.e., taxa should be monophyletic. This fundamental principle for classification of viruses is currently included in the International Committee on Taxonomy of Viruses (ICTV) code only for the rank of species. Secondly, phenotypic and ecological properties of viruses may inform, but not override, evolutionary relatedness in the placement of ranks. Thirdly, alternative classifications that consider phenotypic attributes, such as being vector-borne (e.g., "arboviruses"), infecting a certain type of host (e.g., "mycoviruses," "bacteriophages") or displaying specific pathogenicity (e.g., "human immunodeficiency viruses"), may serve important clinical and regulatory purposes but often create polyphyletic categories that do not reflect evolutionary relationships. Nevertheless, such classifications ought to be maintained if they serve the needs of specific communities or play a practical clinical or regulatory role. However, they should not be considered or called taxonomies. Finally, while an evolution-based framework enables viruses discovered by metagenomics to be incorporated into the ICTV taxonomy, there are essential requirements for quality control of the sequence data used for these assignments. Combined, these four principles will enable future development and expansion of virus taxonomy as the true evolutionary diversity of viruses becomes apparent.

Novel polyomaviruses identified in fecal samples from four carnivore species.

Polyomaviruses are oncogenic viruses that are generally thought to have co-evolved with their hosts. While primate and rodent polyomaviruses are increasingly well-studied, less is known about polyomaviruses that infect other mammals. In an effort to gain insight into polyomaviruses associated with carnivores, we surveyed fecal samples collected in the USA from bobcats (Lynx rufus), pumas (Puma concolor), Canada lynxes (Lynx canadensis), and grizzly bears (Ursus arctos). Using a viral metagenomic approach, we identified six novel polyomavirus genomes. Surprisingly, four of the six genomes showed a phylogenetic relationship to polyomaviruses found in prey animals. These included a putative rabbit polyomavirus from a bobcat fecal sample and two possible deer-trophic polyomaviruses from Canada lynx feces. One polyomavirus found in a grizzly bear sample was found to be phylogenetically distant from previously identified polyomaviruses. Further analysis of the grizzly bear fecal sample showed that it contained anelloviruses that are known to infect pigs, suggesting that the bear might have preyed on a wild or domestic pig. Interestingly, a polyomavirus genome identified in a puma fecal sample was found to be closely related both to raccoon polyomavirus 1 and to Lyon-IARC polyomavirus, the latter of which was originally identified in human saliva and skin swab specimens but has since been found in samples from domestic cats (Felis catus).

Genome Sequences of Anelloviruses, a Genomovirus, Microviruses, Polyomaviruses, and an Unclassified Caudovirus Identified in Vaginal Secretions from South African Adolescents.

Other than for papillomaviruses, there is a paucity of whole-genome sequences for bacteriophages and eukaryote-infecting viruses isolated from the female genital tract. Here, we report the genome sequences of 16 microviruses, 3 anelloviruses, 2 polyomaviruses, 1 genomovirus, and 1 caudovirus that were identified in vaginal secretion samples from adolescents in South Africa.

Interspecies Papillomavirus Type Infection and a Novel Papillomavirus Type in Red Ruffed Lemurs (Varecia rubra).

The Papillomaviridae are a family of vertebrate-infecting viruses of oncogenic potential generally thought to be host species- and tissue-specific. Despite their phylogenetic relatedness to humans, there is a scarcity of data on papillomaviruses (PVs) in speciose non-human primate lineages, particularly the lemuriform primates. Varecia variegata (black-and-white ruffed lemurs) and Varecia rubra (red ruffed lemurs), two closely related species comprising the Varecia genus, are critically endangered with large global captive populations. Varecia variegata papillomavirus (VavPV) types -1 and -2, the first PVs in lemurs with a fully identified genome, were previously characterized from captive V. variegata saliva. To build upon this discovery, saliva samples were collected from captive V. rubra with the following aims: (1) to identify PVs shared between V. variegata and V. rubra and (2) to characterize novel PVs in V. rubra to better understand PV diversity in the lemuriform primates. Three complete PV genomes were determined from V. rubra samples. Two of these PV genomes share 98% L1 nucleotide identity with VavPV2, denoting interspecies infection of V. rubra by VavPV2. This work represents the first reported case of interspecies PV infection amongst the strepsirrhine primates. The third PV genome shares <68% L1 nucleotide identity with that of all PVs. Thus, it represents a new PV species and has been named Varecia rubra papillomavirus 1 (VarPV1). VavPV1, VavPV2, and VarPV1 form a new clade within the Papillomaviridae family, likely representing a novel genus. Future work diversifying sample collection (i.e., lemur host species from multiple genera, sample type, geographic location, and wild populations) is likely to uncover a world of diverse lemur PVs.

Identification of diverse papillomaviruses in captive black-and-white ruffed lemurs (Varecia variegata).

Papillomaviruses (PVs) are host-species-specific and tissue-specific viruses that infect a diverse array of vertebrate hosts, including humans and non-human primates, with varying pathogenic outcomes. Although primate PVs have been studied extensively, no complete genome sequences of PVs from lemurs have been determined to date. Saliva samples from three critically endangered, captive black-and-white ruffed lemurs (Varecia variegata variegata) at the Duke Lemur Center (USA) were analyzed, using high-throughput sequencing, for the presence of oral papillomaviruses. We identified three PVs from two individuals, one of which had a coinfection with two different PVs. Two of the three PVs share 99.6% nucleotide sequence identity, and we have named these isolates "Varecia variegata papillomavirus 1" (VavPV1). The third PV shares ~63% nucleotide sequence identity with VavPV1, and thus, we have named it "Varecia variegata papillomavirus 2" (VavPV2). Based on their E1 + E2 + L1 protein sequence phylogeny, the VavPVs form a distinct clade. This clade likely represents a novel genus, with VavPV1 and VavPV2 belonging to two distinct species. Our findings represent the first complete genome sequences of PVs found in lemuriform primates, with their presence suggesting the potential existence of diverse PVs across the over 100 species of lemurs.

Impact of sample clarification by size exclusion on virus detection and diversity in wastewater-based epidemiology.

The use of wastewater-based epidemiology (WBE) for early detection of virus circulation and response during the SARS-CoV-2 pandemic increased interest in and use of virus concentration protocols that are quick, scalable, and efficient. One such protocol involves sample clarification by size fractionation using either low-speed centrifugation to produce a clarified supernatant or membrane filtration to produce an initial filtrate depleted of solids, eukaryotes and bacterial present in wastewater (WW), followed by concentration of virus particles by ultrafiltration of the above. While this approach has been successful in identifying viruses from WW, it assumes that majority of the viruses of interest should be present in the fraction obtained by ultrafiltration of the initial filtrate, with negligible loss of viral particles and viral diversity. We used WW samples collected in a population of ~700,000 in southwest USA between October 2019 and March 2021, targeting three non-enveloped viruses (enteroviruses [EV], canine picornaviruses [CanPV], and human adenovirus 41 [Ad41]), to evaluate whether size fractionation of WW prior to ultrafiltration leads to appreciable differences in the virus presence and diversity determined. We showed that virus presence or absence in WW samples in both portions (filter trapped solids [FTS] and filtrate) are not consistent with each other. We also found that in cases where virus was detected in both fractions, virus diversity (or types) captured either in FTS or filtrate were not consistent with each other. Hence, preferring one fraction of WW over the other can undermine the capacity of WBE to function as an early warning system and negatively impact the accurate representation of virus presence and diversity in a population.

Genome Sequences of Anelloviruses, Genomovirus, and Papillomavirus Isolated from Nasal Pharyngeal Swabs.

The genome sequences of three anelloviruses (genus Alphatorquevirus), a genomovirus (genus Gemykolovirus), and an unclassified papillomavirus were identified in four human nasopharyngeal swabs, and one was positive for influenza A and one for influenza B virus. The influenza B virus-positive sample had a coinfection with an anellovirus and a papillomavirus.

A parasite outbreak in notothenioid fish in an Antarctic fjord.

Climate changes can promote disease outbreaks, but their nature and potential impacts in remote areas have received little attention. In a hot spot of biodiversity on the West Antarctic Peninsula, which faces among the fastest changing climates on Earth, we captured specimens of two notothenioid fish species affected by large skin tumors at an incidence never before observed in the Southern Ocean. Molecular and histopathological analyses revealed that X-cell parasitic alveolates, members of a genus we call Notoxcellia, are the etiological agent of these tumors. Parasite-specific molecular probes showed that xenomas remained within the skin but largely outgrew host cells in the dermis. We further observed that tumors induced neovascularization in underlying tissue and detrimentally affected host growth and condition. Although many knowledge gaps persist about X-cell disease, including its mode of transmission and life cycle, these findings reveal potentially active biotic threats to vulnerable Antarctic ecosystems.

RNA Virus Gene Signatures Detected in Patients With Cardiomyopathy After Chemotherapy; A Pilot Study.

Background: Viral infections are pervasive and leading causes of myocarditis. Immune-suppression after chemotherapy increases opportunistic infections, but the incidence of virus-induced myocarditis is unknown. Objective: An unbiased, blinded screening for RNA viruses was performed after chemotherapy with correlation to cardiac function. Methods: High-throughput sequencing of RNA isolated from blood samples was analyzed following chemotherapy for hematological malignancies (N = 28) and compared with left ventricular ejection fraction (LVEF). Results: On initial rigorous analysis, low levels of influenza orthomyxovirus and avian paramyxovirus sequences were detectable, but without significant correlation to LVEF (r = 0.208). A secondary broad data mining analysis for virus sequences, without filtering human sequences, detected significant correlations for paramyxovirus with LVEF after chemotherapy (r = 0.592, P < 0.0096). Correlations were similar for LVEF pre- and post- chemotherapy for orthomyxovirus (R = 0.483, P < 0.0421). Retrovirus detection also correlated with LVEF post (r = 0.453, p < 0.0591), but not pre-chemotherapy, but is suspect due to potential host contamination. Detectable phage and anellovirus had no correlation. Combined sequence reads (all viruses) demonstrated significant correlation (r = 0.621, P < 0.0078). Reduced LVEF was not associated with chemotherapy (P = NS). Conclusions: This is the first report of RNA virus screening in circulating blood and association with changes in cardiac function among patients post chemotherapy, using unbiased, blinded, high-throughput sequencing. Influenza orthomyxovirus, avian paramyxovirus and retrovirus sequences were detectable in patients with reduced LVEF. Further analysis for RNA virus infections in patients with cardiomyopathy after chemotherapy is warranted.

Identification of a Novel Myxoma Virus C7-Like Host Range Factor That Enabled a Species Leap from Rabbits to Hares.

Myxoma virus (MYXV) is naturally found in rabbit Sylvilagus species and is known to cause lethal myxomatosis in European rabbits (Oryctolagus cuniculus). In 2019, an MYXV strain (MYXV strain Toledo [MYXV-Tol]) causing myxomatosis-like disease in Iberian hares (Lepus granatensis) was identified. MYXV-Tol acquired a recombinant region of ∼2.8 kb harboring several new genes, including a novel host range gene (M159) that we show to be an orthologous member of the vaccinia virus C7 host range family. Here, to test whether M159 alone has enabled MYXV to alter its host range to Iberian hares, several recombinant viruses were generated, including an MYXV-Tol ΔM159 (knockout) strain. While MYXV-Tol underwent fully productive infection in hare HN-R cells, neither the wild-type MYXV-Lau strain (lacking M159) nor vMyxTol-ΔM159 (deleted for M159) was able to infect and replicate, showing that the ability of MYXV-Tol to infect these cells and replicate depends on the presence of M159. Similar to other C7L family members, M159 was shown to be expressed as an early/late gene but was translocated into the nucleus at later time points, indicating that further studies are needed to elucidate its role in the nucleus. Finally, in rabbit cells, the M159 protein did not contribute to increased replication but was able to upregulate the replication levels of MYXV in nonpermissive and semipermissive human cancer cells, suggesting that the M159-targeted pathway is conserved across mammalian species. Altogether, these observations demonstrate that the M159 protein plays a critical role in determining the host specificity of MYXV-Tol in hare and human cells by imparting new host range functions. IMPORTANCE The coevolution of European rabbit populations and MYXV is a textbook example of an arms race between a pathogen and a host. Recently, a recombinant MYXV (MYXV-Tol) crossed the species barrier by jumping from leporid species to another species, causing lethal myxomatosis-like disease. Given the highly pathogenic nature of this new virus in hares and the incidences of other poxvirus cross-species spillovers into other animals, including humans, it is important to understand how and why MYXV-Tol was able to become virulent in a new host species. The results presented clearly demonstrate that M159 is the key factor allowing MYXV-Tol replication in hare cells by imparting new host range functions. These results have the potential to improve current knowledge about the virulence of poxviruses and provide a platform to better understand the new MYXV-Tol, rendering the virus capable of leaping into a new host species.

Coevolutionary Analysis Implicates Toll-Like Receptor 9 in Papillomavirus Restriction.

Upon infection, DNA viruses can be sensed by pattern recognition receptors (PRRs), leading to the activation of type I and III interferons to block infection. Therefore, viruses must inhibit these signaling pathways, avoid being detected, or both. Papillomavirus virions are trafficked from early endosomes to the Golgi apparatus and wait for the onset of mitosis to complete nuclear entry. This unique subcellular trafficking strategy avoids detection by cytoplasmic PRRs, a property that may contribute to the establishment of infection. However, as the capsid uncoats within acidic endosomal compartments, the viral DNA may be exposed to detection by Toll-like receptor 9 (TLR9). In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we showed that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts. Finally, we demonstrated that the cancer-associated human papillomaviruses show a reduction in CpG dinucleotides within a TLR9 recognition complex. IMPORTANCE Viruses must avoid detection by the innate immune system. In this study, we characterized two new papillomaviruses from bats and used molecular archeology to demonstrate that their genomes altered their nucleotide compositions to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we demonstrated that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for coevolution between papillomaviruses and their hosts.