Multiple evolutionary origins of bat papillomaviruses

Viral Metagenomic Data Analyses of Five New World Bat Species from Argentina: Identification of 35 Novel DNA Viruses

Bats are natural reservoirs of a variety of zoonotic viruses, many of which cause severe human diseases. Characterizing viruses of bats inhabiting different geographical regions is important for understanding their viral diversity and for detecting viral spillovers between animal species. Herein, the diversity of DNA viruses of five arthropodophagous bat species from Argentina was investigated using metagenomics. Fecal samples of 29 individuals from five species (Tadarida brasiliensis, Molossus molossus, Eumops bonariensis, Eumops patagonicus, and Eptesicus diminutus) living at two different geographical locations, were investigated. Enriched viral DNA was sequenced using Illumina MiSeq, and the reads were trimmed and filtered using several bioinformatic approaches. The resulting nucleotide sequences were subjected to viral taxonomic classification. In total, 4,520,370 read pairs were sequestered by sequencing, and 21.1% of them mapped to viral taxa. Circoviridae and Genomoviridae were the most prevalent among vertebrate viral families in all bat species included in this study. Samples from the T. brasiliensis colony exhibited lower viral diversity than samples from other species of New World bats. We characterized 35 complete genome sequences of novel viruses. These findings provide new insights into the global diversity of bat viruses in poorly studied species, contributing to prevention of emerging zoonotic diseases and to conservation policies for endangered species.

A Preliminary Study of the Virome of the South American Free-Tailed Bats (Tadarida brasiliensis) and Identification of Two Novel Mammalian Viruses

Bats provide important ecosystem services as pollinators, seed dispersers, and/or insect controllers, but they have also been found harboring different viruses with zoonotic potential. Virome studies in bats distributed in Asia, Africa, Europe, and North America have increased dramatically over the past decade, whereas information on viruses infecting South American species is scarce. We explored the virome of Tadarida brasiliensis, an insectivorous New World bat species inhabiting a maternity colony in Rosario (Argentina), by a metagenomic approach. The analysis of five pooled oral/anal swab samples indicated the presence of 43 different taxonomic viral families infecting a wide range of hosts. By conventional nucleic acid detection techniques and/or bioinformatics approaches, the genomes of two novel viruses were completely covered clustering into the Papillomaviridae (Tadarida brasiliensis papillomavirus type 1, TbraPV1) and Genomoviridae (Tadarida brasiliensis gemykibivirus 1, TbGkyV1) families. TbraPV1 is the first papillomavirus type identified in this host and the prototype of a novel genus. TbGkyV1 is the first genomovirus reported in New World bats and constitutes a new species within the genus Gemykibivirus. Our findings extend the knowledge about oral/anal viromes of a South American bat species and contribute to understand the evolution and genetic diversity of the novel characterized viruses.

Origin and evolution of papillomavirus (onco)genes and genomes

Papillomaviruses (PVs) are ancient viruses infecting vertebrates, from fishes to mammals. Although the genomes of PVs are small and show conserved synteny, PVs display large genotypic diversity and ample variation in the phenotypic presentation of the infection. Most PV genomes contain two small early genes E6 and E7. In a bunch of closely related human papillomaviruses (HPVs), the E6 and E7 proteins provide the viruses with oncogenic potential. The recent discoveries of PVs without E6 and E7 in different fish species place a new root on the PV tree, and suggest that ancestral PVs consisted of the minimal PV backbone E1-E2-L2-L1. Bayesian phylogenetic analyses date the most recent common ancestor of the PV backbone to 424 million years ago (Ma). Common ancestry tests on extant E6 and E7 genes indicate that they share a common ancestor dating back to at least 184 Ma. In AlphaPVs infecting Old World monkeys and apes, the appearance of the E5 oncogene 53-58 Ma concurred with (i) a significant increase in substitution rate, (ii) a basal radiation and (iii) key gain of functions in E6 and E7. This series of events was instrumental to construct the extant phenotype of oncogenic HPVs. Our results assemble the current knowledge on PV diversity and present an ancient evolutionary timeline punctuated by evolutionary innovations in the history of this successful viral family. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Diversity and Evolution of Viral Pathogen Community in Cave Nectar Bats (Eonycteris spelaea)

Bats are unique mammals, exhibit distinctive life history traits and have unique immunological approaches to suppression of viral diseases upon infection. High-throughput next-generation sequencing has been used in characterizing the virome of different bat species. The cave nectar bat, Eonycteris spelaea, has a broad geographical range across Southeast Asia, India and southern China, however, little is known about their involvement in virus transmission. Here we investigate the diversity and abundance of viral communities from a colony of Eonycteris spelaea residing in Singapore. Our results detected 47 and 22 different virus families from bat fecal and urine samples, respectively. Among these, we identify a large number of virus families including Adenoviridae, Flaviviridae, Reoviridae, Papillomaviridae, Paramyxoviridae, Parvoviridae, Picornaviridae, and Polyomaviridae. In most cases, viral sequences from Eonycteris spelaea are genetically related to a group of bat viruses from other bat genera (e.g., Eidolon, Miniopterus, Rhinolophus and Rousettus). The results of this study improve our knowledge of the host range, spread and evolution of several important viral pathogens. More significantly, our findings provide a baseline to study the temporal patterns of virus shedding and how they correlate with bat phenological trends.

Genomic sequencing of a virus representing a novel type within the species Dyopipapillomavirus 1 in an Indian River Lagoon bottlenose dolphin

Fecal samples collected from free-ranging Atlantic bottlenose dolphins (BDs) in the Indian River Lagoon of Florida were processed for viral discovery using a next-generation sequencing (NGS) approach. A 693-bp contig identified in the NGS data was nearly identical to the partial L1 gene sequence of a papillomavirus (PV) previously found in a penile papilloma in a killer whale (Orcinus orca). Based on this partial bottlenose dolphin papillomavirus (BDPV) sequence, a nested inverse PCR and primer-walking strategy was employed to generate the complete genome sequence. The full BDPV genome consisted of 7299 bp and displayed a typical PV genome organization. The BDPV E6 protein contained a PDZ-binding motif, which has been shown to be involved in carcinogenic transformation involving high-risk genital human PVs. Screening of 12 individual fecal samples using a specific endpoint PCR assay revealed that the feces from a single female BD displaying a genital papilloma was positive for the BDPV. Genetic analysis indicated that this BDPV (Tursiops truncatus papillomavirus 8; TtPV8) is a new type of Dyopipapillomavirus 1, previously sequenced from an isolate obtained from a penile papilloma in a harbor porpoise (Phocoena phocoena). Although only a partial L1 sequence has been determined for a PV detected in a killer whale genital papilloma, our finding of a nearly identical sequence in an Atlantic BD may indicate that members of this viral species are capable of host jumping. Future work is needed to determine if this virus is a high-risk PV that is capable of inducing carcinogenic transformation and whether it poses a significant health risk to wild delphinid populations.

Ancient DNA provides evidence of 27,000-year-old papillomavirus infection and long-term codivergence with rodents

The long-term evolutionary history of many viral lineages is poorly understood. Novel sources of ancient DNA combined with phylogenetic analyses can provide insight into the time scale of virus evolution. Here we report viral sequences from ancient North American packrat middens. We screened samples up to 27,000-years old and found evidence of papillomavirus (PV) infection in Neotoma cinerea (Bushy-tailed packrat). Phylogenetic analysis placed the PV sequences in a clade with other previously published PV sequences isolated from rodents. Concordance between the host and virus tree topologies along with a correlation in branch lengths suggests a shared evolutionary history between rodents and PVs. Based on host divergence times, PVs have likely been circulating in rodents for at least 17 million years. These results have implications for our understanding of PV evolution and for further research with ancient DNA from Neotoma middens.

A single bat species in Cameroon harbors multiple highly divergent papillomaviruses in stool identified by metagenomics analysis

A number of PVs have been described in bats but to the best of our knowledge not from feces. Using a previously described NetoVIR protocol, Eidolon helvum pooled fecal samples (Eh) were treated and sequenced by Illumina next generation sequencing technology. Two complete genomes of novel PVs (EhPV2 and EhPV3) and 3 partial sequences (BATPV61, BATPV890a and BATPV890b) were obtained and analysis showed that the EhPV2 and EhPV3 major capsid proteins cluster with and share 60-64% nucleotide identity with that of Rousettus aegyptiacus PV1, thus representing new species of PVs within the genus Psipapillomavirus. The other PVs clustered in different branches of our phylogenetic tree and may potentially represent novel species and/or genera. This points to the vast diversity of PVs in bats and in Eidolon helvum bats in particular, therefore adding support to the current concept that PV evolution is more complex than merely strict PV-host co-evolution.

Novel papillomaviruses in free-ranging Iberian bats: no virus-host co-evolution, no strict host specificity, and hints for recombination

Papillomaviruses (PVs) are widespread pathogens. However, the extent of PV infections in bats remains largely unknown. This work represents the first comprehensive study of PVs in Iberian bats. We identified four novel PVs in the mucosa of free-ranging Eptesicus serotinus (EserPV1, EserPV2, and EserPV3) and Rhinolophus ferrumequinum (RferPV1) individuals and analyzed their phylogenetic relationships within the viral family. We further assessed their prevalence in different populations of E. serotinus and its close relative E. isabellinus. Although it is frequent to read that PVs co-evolve with their host, that PVs are highly species-specific, and that PVs do not usually recombine, our results suggest otherwise. First, strict virus–host co-evolution is rejected by the existence of five, distantly related bat PV lineages and by the lack of congruence between bats and bat PVs phylogenies. Second, the ability of EserPV2 and EserPV3 to infect two different bat species (E. serotinus and E. isabellinus) argues against strict host specificity. Finally, the description of a second noncoding region in the RferPV1 genome reinforces the view of an increased susceptibility to recombination in the E2-L2 genomic region. These findings prompt the question of whether the prevailing paradigms regarding PVs evolution should be reconsidered.

A novel papillomavirus isolated from a nasal neoplasia in an Italian free-ranging chamois (Rupicapra r. rupicapra)

Most amniotes are the hosts of many, distantly related papillomaviruses (PVs). Infection by PVs can be asymptomatic, or lead instead to benign or malignant lesions. However, PVs infecting animals and associated with malignancies are still largely understudied. In the present study, we communicate the complete genome of a novel PV found in a nasal neoplasia of a free-ranging alpine chamois (Rupicapra r. rupicapra) in an Italian national park. Long-PCR and cloning approaches followed for Sanger sequencing were used to identify the first PV found in chamois. The genome of the novel virus - RrupPV1 - of 7256 bp in length, presents the classical PV structure, and lacks the interE2-L2 region that hosts the E5 gene in AlphaPVs and in DeltaPVs. The nucleotide identity percentage of the L1 ORF, places RrupPV1 together with OaPV3 in the same genus. The latter is a PV isolated from a squamous cell carcinoma in sheep in Sardinia. Full-genome phylogenetic reconstructions suggest that these two viruses are sister taxa, and that both of them are very distantly related to any other known PV. Many cetartiodactyl species are infected by non-monophyletic PVs. Our results exemplify further the multiple links between the infection by certain, distantly related PVs and the development of diverse cancers in animals and highlight the need of a systematic search of oncogenic and non-oncogenic animal PVs.