Evidence of recombination and positive selection in cetacean papillomaviruses

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.

Viral Diversity and Epidemiology in Critically Endangered Yangtze Finless Porpoises (Neophocaena asiaeorientalis asiaeorientalis)

The Yangtze finless porpoise (YFP) (Neophocaena asiaeorientalis asiaeorientalis) is a critically endangered freshwater cetacean, with about 1,249 individuals thought to be left in the wild. However, viral entities and viral diseases of YFPs remain obscure. In this study, anal swabs for virome analysis were collected during the physical examination of YFPs in the Tian-E-Zhou Oxbow (TEO) ex situ reserve. A total of 19 eukaryotic viral species belonging to 9 families, including Papillomaviridae, Herpesviridae, Picornaviridae, Picobirnaviridae, Caliciviridae, Retroviridae, Parvoviridae, Virgaviridae, and Narnaviridae, and other unclassified viruses were identified based on metasequencing. Among these detected viruses, a novel herpesvirus (NaHV), two different kobuviruses (NaKV1-2), and six different papillomaviruses (NaPV1 to -6) were considered potential risks to YFPs and confirmed by PCR or reverse transcription-PCR (RT-PCR). Most YFPs sampled were found to harbor one or more kinds of detected viral genomes (52/58 [89.7%]). Surveillance results demonstrated that kobuvirus and herpesvirus displayed obvious age distribution and PVs showed significant gender difference in YFPs. According to species demarcation criteria in individual genera in Papillomaviridae, two novel species (referred to as Omikronpapillomavirus 2 and 3) and four novel isolates of PV were identified in YFPs. Further evolutionary analysis suggested that NaPVs would occupy the mucosal niche and that virus-host codivergence mixed with duplications and host-switching events drives the evolution of cetacean PVs. Divergence times of PVs in YFP and other cetacean reflect the incipient speciation of YFPs. In summary, our findings revealed the potential viral entities, their prevalence, and their evolutionary history in YFPs, which raises an important issue regarding effects of viral infection on the fitness of YFPs. IMPORTANCE The Yangtze finless porpoise (YFP) is the only cetacean species in freshwater following the functional extinction of the baiji (Lipotes vexillifer). Health management, disease treatment, and other special measures are important for maintaining the existing YFP populations, especially in in situ and ex situ reserves. The discovery of potential viral entities and their prevalence in YFPs raises an important issue regarding the effects of viral infection on the fitness of YFPs and may contribute to the conservation of YFPs. The evolutionary history of papillomaviruses in YFP and other cetaceans reflects the phylogeny of their hosts and supports the status of incipient species, opening a window to investigate the evolutionary adaptation of cetaceans to freshwater as well as their phylogeny to remedy the deficiency of fossil evidence.

Identification of two novel papillomaviruses in belugas

Introduction

Papillomaviruses (PVs) can cause hyperplasia in the skin and mucous membranes of humans, mammals, and non-mammalian animals, and are a significant risk factor for cervical and genital cancers.

Methods

Using next-generation sequencing (NGS), we identified two novel strains of papillomavirus, PV-HMU-1 and PV-HMU-2, in swabs taken from belugas (Delphinapterus leucas) at Polar Ocean Parks in Qingdao and Dalian.

Results

We amplified the complete genomes of both strains and screened ten belugas and one false killer whale (Pseudorca crassidens) for the late gene (L1) to determine the infection rate. In Qingdao, 50% of the two sampled belugas were infected with PV-HMU-1, while the false killer whale was negative. In Dalian, 71% of the eight sampled belugas were infected with PV-HMU-2. In their L1 genes, PV-HMU-1 and PV-HMU-2 showed 64.99 and 68.12% amino acid identity, respectively, with other members of Papillomaviridae. Phylogenetic analysis of combinatorial amino acid sequences revealed that PV-HMU-1 and PV-HMU-2 clustered with other known dolphin PVs but formed distinct branches. PVs carried by belugas were proposed as novel species under Firstpapillomavirinae.

Conclusion

The discovery of these two novel PVs enhances our understanding of the genetic diversity of papillomaviruses and their impact on the beluga population.

Characterization of a novel papillomavirus identified from a whale (Delphinapterus leucas) pharyngeal metagenomic library

Here, using viral metagenomic method, a novel whale papillomavirus (temporarily named wPV, GenBank accession number OP856597) was discovered in a whale (Delphinapterus leucas) pharyngeal metagenomic library. The complete genome size of wPV is 7179 bp, with GC content of 54.4% and a nucleotide composition of 23.4% A, 22.3% T, 28.4% G, and 25.9% C. The viral genome has a typical papillomavirus organization pattern, and five ORFs were predicted, including two late genes encoding L1 and L2, and three early genes encoding E1, E2, and E6. Pairwise sequence comparison and phylogenetic analysis based on the L1 gene sequence indicated that wPV may be a novel species within genus Dyodeltapapillomavirus. In addition, the E2 region of wPV was predicted to have a potential recombination event. The discovery of this novel papillomavirus increases our understanding of the viral ecology of marine mammals, providing insights into possible future infectious diseases.

Cytologic features of oral squamous cell carcinoma in an Indo-Pacific bottlenose dolphin (Tursiops aduncus): Papanicolaou stain and immunocytochemistry using liquid-based cytology

Although oral cytology using Papanicolaou stain is useful for the early detection of oral premalignant lesions and squamous cell carcinoma (SCC) in people, little work has been conducted on this topic in veterinary medicine. This paper describes the features of oral cytology using Papanicolaou stain and immunocytochemistry on liquid-based cytology slides in a case of oral SCC in an Indo-Pacific bottlenose dolphin (Tursiops aduncus). In this case, dysplastic cells with koilocyte-like changes and SCC cells were identified using the Papanicolaou stain. These cells were positive for p53 using an immunocytochemistry analysis. A cytologic diagnosis of SCC was made. We believe that the early detection of premalignant oral lesions and SCC in dolphins can be significantly improved with cytology using liquid-based cytology, Papanicolaou staining, and immunocytochemistry.

First Report of Phodopus sungorus Papillomavirus Type 1 Infection in Roborovski Hamsters (Phodopus roborovskii)

Papillomaviruses (PVs) are considered highly species-specific with cospeciation as the main driving force in their evolution. However, a recent increase in the available PV genome sequences has revealed inconsistencies in virus–host phylogenies, which could be explained by adaptive radiation, recombination, host-switching events and a broad PV host range. Unfortunately, with a relatively low number of animal PVs characterized, understanding these incongruities remains elusive. To improve knowledge of biology and the spread of animal PV, we collected 60 swabs of the anogenital and head and neck regions from a healthy colony of 30 Roborovski hamsters (Phodopus roborovskii) and detected PVs in 44/60 (73.3%) hamster samples. This is the first report of PV infection in Roborovski hamsters. Moreover, Phodopus sungorus papillomavirus type 1 (PsuPV1), previously characterized in Siberian hamsters (Phodopus sungorus), was the only PV detected in Roborovski hamsters. In addition, after a detailed literature search, review and summary of published evidence and construction of a tanglegram linking the cladograms of PVs and their hosts, our findings were discussed in the context of available knowledge on PVs described in at least two different host species.

Anogenital-Associated Papillomaviruses in Animals: Focusing on Bos taurus Papillomaviruses

In contrast to the diverse studies on human papillomaviruses (HPVs), information on animal PVs associated with anogenital lesions is limited. In the animal kingdom, papillomas occur more commonly in cattle than in any other animals, and diverse types of Bos taurus papillomaviruses (BPVs) exist, including the very recently discovered BPV type 29 (BPV29). From this perspective, we will review previous studies describing PV types associated with anogenitals in animals, with a focus on BPVs. To date, two classical BPV types, classified into Deltapapillomavirus (BPV1 and BPV2) and Dyokappapapillomavirus (BPV22), and two novel Xipapillomaviruses (BPV28 and BPV29) have been identified from anogenital lesions and tissues of the domestic cow. Due to the limited reports describing anogenital-associated PVs in animals, the relationships between their phylogenetic and pathogenetic properties are still undiscovered. Animal studies are valuable not only for the veterinary field but also for human medicine, as animal diseases have been shown to mimic human diseases. Studies of anogenital-associated PVs in animals have a positive impact on various research fields.

Genomic and phylogenetic characterization of ChPV2, a novel goat PV closely related to the Xi-PV1 species infecting bovines

BACKGROUND

Papillomaviruses (PVs) infecting artiodactyls are very diverse, and only second in number to PVs infecting primates. PVs associated to lesions in economically important ruminant species have been isolated from cattle and sheep.

METHODS

Potential PV DNA from teat lesions of a Damascus goat was isolated, cloned and sequenced. The PV genome was analyzed using bioinformatics approaches to detect open reading frames and to predict potential features of encoded proteins as well as putative regulatory elements. Sequence comparison and phylogenetic analyses using the concatenated E1E2L2L1 nucleotide and amino acid alignments was used to reveal the relationship of the new PV to the known PV diversity and its closest relevants.

RESULTS

We isolated and characterized the full-genome of novel Capra hircus papillomavirus. We identified the E6, E7, E1, E2, L2, L1 open reading frames with protein coding potential and putative active elements in the ChPV2 proteins and putative regulatory genome elements. Sequence similarities of L1 and phylogenetic analyses using concatenated E1E2L2L1 nucleotide and amino acid alignments suggest the classification as a new PV type designated ChPV2 with a phylogenetic position within the XiPV genus, basal to the XiPV1 species. ChPV2 is not closely related to ChPV1, the other known goat PV isolated from healthy skin, although both of them belong confidently into a clade composed of PVs infecting cervids and bovids. Interestingly, ChPV2 contains an E6 open reading frame whereas all closely related PVs do not CONCLUSION: ChPV2 is a novel goat PV closely related to the Xi-PV1 species infecting bovines. Phylogenetic relationships and genome architecture of ChPV2 and closely related PV types suggest at least two independent E6 losses within the XiPV clade.

Characterization of Bovine papillomavirus 28 (BPV28) and a novel genotype BPV29 associated with vulval papillomas in cattle

Papillomavirus (PV) infections are associated with the development of cutaneous and mucosal tumors in humans and various animal species. In humans, infection of high-risk human PVs (HPVs) causes anogenital cancers, while in animals, anogenital-associated PVs are not well understood. Among animal PVs, Bos taurus PVs (BPVs) have the most diverse genotypes, up to 28 of them. The present study will report two unique BPVs identified in vulval papilloma lesions from two Holstein Friesian cattle by conventional PCR and sequencing. In the first case, BPV28 harboring two L1 open reading frames (ORFs) due to a five-nucleotide deletion was identified. In the second case, histologically diagnosed as papilloma, an unclassified BPV genotype was detected. However, in both cases, the immunohistochemistry against PV antigen was negative. The full genome of the unclassified BPV was amplified by inverse PCR and analyzed by genome-walking sequencing. The L1 nucleotide sequence was most identical to BPV genotype 6 (BPV6), showing 78 % identity, indicating that this novel BPV should be classified as species Xipapillomavirus 1, genotype BPV29. The mRNA expression of three early genes (E1, E2, E10), but not L1, was confirmed in both BPV28- and BPV29-detected papilloma lesions. The present study suggests the involvement of novel types of BPV in vulval papilloma. The alteration of BPV28 pathogenicity due to the frameshift mutation of L1 needs to be elucidated in the future.

Papillomaviruses infecting cetaceans exhibit signs of genome adaptation following a recombination event

Papillomaviruses (PVs) have evolved through a complex evolutionary scenario where virus-host co-evolution alone is not enough to explain the phenotypic and genotypic PV diversity observed today. Other evolutionary processes, such as host switch and recombination, also appear to play an important role in PV evolution. In this study, we have examined the genomic impact of a recombination event between distantly related PVs infecting Cetartiodactyla (even-toed ungulates and cetaceans). Our phylogenetic analyses suggest that one single recombination was responsible for the generation of extant 'chimeric' PV genomes infecting cetaceans. By correlating the phylogenetic relationships to the genomic content, we observed important differences between the recombinant and non-recombinant cetartiodactyle PV genomes. Notably, recombinant PVs contain a unique set of conserved motifs in the upstream regulatory region (URR). We interpret these regulatory changes as an adaptive response to drastic changes in the PV genome. In terms of codon usage preferences (CUPrefs), we did not detect any particular differences between orthologous open reading frames in recombinant and non-recombinant PVs. Instead, our results are in line with previous observations suggesting that CUPrefs in PVs are rather linked to gene expression patterns as well as to gene function. We show that the non-coding URR of PVs infecting cetaceans, the central regulatory element in these viruses, exhibits signs of adaptation following a recombination event. Our results suggest that also in PVs, the evolution of gene regulation can play an important role in speciation and adaptation to novel environments.

Recombination Between High-Risk Human Papillomaviruses and Non-Human Primate Papillomaviruses: Evidence of Ancient Host Switching Among Alphapapillomaviruses

We use all the currently known 405 Papillomavirus (PV) sequences, 343 curated PV sequences from both humans and animals from the PAVE data base, to analyse the recombination dynamics of these viruses at the whole genome levels. After showing some evidence of human and non-human primate PV recombination, we report a comprehensive recombination analysis of all currently known 82 Alphapapillomaviruses (Alpha-PVs). We carried out an exploratory study and found novel recombination events between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 (PpPV1) Papillomaviruses. This is the first evidence of interactions between PVs from different hosts and hence postulates the likelihood of ancient host switching among Alpha-PVs. Notwithstanding these results should be interpreted with caution because the major and minor parents indicated by RDP4 program are simply the sequences in the alignment that most closely resemble the actual parents. We found statistically significant differences between the phylogenies of the PV sequences with recombination regions and PV sequences without recombination regions using the Shimodaira–Hasegawa phylogenetic incongruence testing. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts before speciation of primate species. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. However, a clearer picture of the PVs evolutionary landscape can only be achieved with the incremental discovery of PV sequences, especially from the animal kingdom.

Genome Plasticity in Papillomaviruses and De Novo Emergence of E5 Oncogenes

The clinical presentations of papillomavirus (PV) infections come in many different flavors. While most PVs are part of a healthy skin microbiota and are not associated to physical lesions, other PVs cause benign lesions, and only a handful of PVs are associated to malignant transformations linked to the specific activities of the E5, E6, and E7 oncogenes. The functions and origin of E5 remain to be elucidated. These E5 open reading frames (ORFs) are present in the genomes of a few polyphyletic PV lineages, located between the early and the late viral gene cassettes. We have computationally assessed whether these E5 ORFs have a common origin and whether they display the properties of a genuine gene. Our results suggest that during the evolution of Papillomaviridae, at least four events lead to the presence of a long noncoding DNA stretch between the E2 and the L2 genes. In three of these events, the novel regions evolved coding capacity, becoming the extant E5 ORFs. We then focused on the evolution of the E5 genes in AlphaPVs infecting primates. The sharp match between the type of E5 protein encoded in AlphaPVs and the infection phenotype (cutaneous warts, genital warts, or anogenital cancers) supports the role of E5 in the differential oncogenic potential of these PVs. In our analyses, the best-supported scenario is that the five types of extant E5 proteins within the AlphaPV genomes may not have a common ancestor. However, the chemical similarities between E5s regarding amino acid composition prevent us from confidently rejecting the model of a common origin. Our evolutionary interpretation is that an originally noncoding region entered the genome of the ancestral AlphaPVs. This genetic novelty allowed to explore novel transcription potential, triggering an adaptive radiation that yielded three main viral lineages encoding for different E5 proteins, displaying distinct infection phenotypes. Overall, our results provide an evolutionary scenario for the de novo emergence of viral genes and illustrate the impact of such genotypic novelty in the phenotypic diversity of the viral infections.

Comparative Innate and Adaptive Immune Responses in Atlantic Bottlenose Dolphins (Tursiops truncatus) With Viral, Bacterial, and Fungal Infections

Free-ranging Atlantic bottlenose dolphins (n = 360) from two southeastern U.S. estuarine sites were given comprehensive health examinations between 2003 and 2015 as part of a multi-disciplinary research project focused on individual and population health. The study sites (and sample sizes) included the Indian River Lagoon (IRL), Florida, USA (n = 246) and Charleston harbor and associated rivers (CHS), South Carolina, USA (n = 114). Results of a suite of clinicoimmunopathologic tests revealed that both populations have a high prevalence of infectious and neoplastic disease and a variety of abnormalities of their innate and adaptive immune systems. Subclinical infections with cetacean morbillivirus and Chlamydiaceae were detected serologically. Clinical evidence of orogenital papillomatosis was supported by the detection of a new strain of dolphin papillomavirus and herpesvirus by molecular pathology. Dolphins with cutaneous lobomycosis/lacaziasis were subsequently shown to be infected with a novel, uncultivated strain of Paracoccidioides brasiliensis, now established as the etiologic agent of this enigmatic disease in dolphins. In this review, innate and adaptive immunologic responses are compared between healthy dolphins and those with clinical and/or immunopathologic evidence of infection with these specific viral, bacterial, and fungal pathogens. A wide range of immunologic host responses was associated with each pathogen, reflecting the dynamic and complex interplay between the innate, humoral, and cell-mediated immune systems in the dolphin. Collectively, these studies document the comparative innate and adaptive immune responses to various types of infectious diseases in free-ranging Atlantic bottlenose dolphins. Evaluation of the type, pattern, and degree of immunologic response to these pathogens provides novel insight on disease immunopathogenesis in this species and as a comparative model. Importantly, the data suggest that in some cases infection may be associated with subclinical immunopathologic perturbations that could impact overall individual and population health.

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'.

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.

Complete Genome Sequencing of a Novel Type of Omikronpapillomavirus 1 in Indian River Lagoon Bottlenose Dolphins (Tursiops truncatus)

The genome sequence of a papillomavirus was determined from fecal samples collected from bottlenose dolphins in the Indian River Lagoon, FL. The genome was 7,772 bp and displayed a typical papillomavirus genome organization. Phylogenetic analysis supported the bottlenose dolphin papillomavirus as being a novel type of Omikronpapillomavirus1.

Diverse papillomaviruses identified in Weddell seals

Papillomaviridae is a diverse family of circular, double-stranded DNA (dsDNA) viruses that infect a broad range of mammalian, avian and fish hosts. While papillomaviruses have been characterized most extensively in humans, the study of non-human papillomaviruses has contributed greatly to our understanding of their pathogenicity and evolution. Using high-throughput sequencing approaches, we identified 7 novel papillomaviruses from vaginal swabs collected from 81 adult female Weddell seals (Leptonychotes weddellii) in the Ross Sea of Antarctica between 2014-2017. These seven papillomavirus genomes were amplified from seven individual seals, and six of the seven genomes represented novel species with distinct evolutionary lineages. This highlights the diversity of papillomaviruses among the relatively small number of Weddell seal samples tested. Viruses associated with large vertebrates are poorly studied in Antarctica, and this study adds information about papillomaviruses associated with Weddell seals and contributes to our understanding of the evolutionary history of papillomaviruses.

Association of papillomavirus E6 proteins with either MAML1 or E6AP clusters E6 proteins by structure, function, and evolutionary relatedness

Papillomavirus E6 proteins bind to LXXLL peptide motifs displayed on targeted cellular proteins. Alpha genus HPV E6 proteins associate with the cellular ubiquitin ligase E6AP (UBE3A), by binding to an LXXLL peptide (ELTLQELLGEE) displayed by E6AP, thereby stimulating E6AP ubiquitin ligase activity. Beta, Gamma, and Delta genera E6 proteins bind a similar LXXLL peptide (WMSDLDDLLGS) on the cellular transcriptional co-activator MAML1 and thereby repress Notch signaling. We expressed 45 different animal and human E6 proteins from diverse papillomavirus genera to ascertain the overall preference of E6 proteins for E6AP or MAML1. E6 proteins from all HPV genera except Alpha preferentially interacted with MAML1 over E6AP. Among animal papillomaviruses, E6 proteins from certain ungulate (SsPV1 from pigs) and cetacean (porpoises and dolphins) hosts functionally resembled Alpha genus HPV by binding and targeting the degradation of E6AP. Beta genus HPV E6 proteins functionally clustered with Delta, Pi, Tau, Gamma, Chi, Mu, Lambda, Iota, Dyokappa, Rho, and Dyolambda E6 proteins to bind and repress MAML1. None of the tested E6 proteins physically and functionally interacted with both MAML1 and E6AP, indicating an evolutionary split. Further, interaction of an E6 protein was insufficient to activate degradation of E6AP, indicating that E6 proteins that target E6AP co-evolved to separately acquire both binding and triggering of ubiquitin ligase activation. E6 proteins with similar biological function clustered together in phylogenetic trees and shared structural features. This suggests that the divergence of E6 proteins from either MAML1 or E6AP binding preference is a major event in papillomavirus evolution.

Molecular characterization of novel mucosotropic papillomaviruses from a Florida manatee (Trichechus manatus latirostris)

We isolated two new manatee papillomavirus (PV) types, TmPV3 and TmPV4, from a Florida manatee (Trichechus manatus latirostris). Two PV types were previously isolated from this species. TmPV1 is widely dispersed amongst manatees and a close-to-root PV; not much is known about TmPV2. The genomes of TmPV3 and TmPV4 were 7622 and 7771 bp in size, respectively. Both PVs had a genomic organization characteristic of all PVs, with one non-coding region and seven ORFs, including the E7 ORF that is absent in other cetacean PVs. Although these PVs were isolated from separate genital lesions of the same manatee, an enlarged E2/E4 ORF was found only in the TmPV4 genome. The full genome and L1 sequence similarities between TmPV3 and TmPV4 were 63.2 and 70.3 %, respectively. These genomes shared only 49.1 and 50.2 % similarity with TmPV1. The pairwise alignment of L1 nucleotide sequences indicated that the two new PVs nested in a monophyletic group of the genus Rhopapillomavirus, together with the cutaneotropic TmPV1 and TmPV2.

Characterization of novel human papillomavirus types 157, 158 and 205 from healthy skin and recombination analysis in genus γ-Papillomavirus

Gammapapillomavirus (γ-PV) is a diverse and rapidly expanding genus, currently consisting of 79 fully characterized human PV (HPV) types. In this study, three novel types, HPV157, HPV158 and HPV205, obtained from healthy sun-exposed skin of two immunocompetent individuals, were amplified by the "Hanging droplet" long PCR technique, cloned, sequenced and characterized. HPV157, HPV158 and HPV205 genomes comprise 7154-bp, 7192-bp and 7298-bp, respectively, and contain four early (E1, E2, E6 and E7) and two late genes (L1 and L2). Phylogenetic analysis of the L1 ORF placed all novel types within the γ-PV genus: HPV157 was classified as a new member of species γ-12 while HPV158 and HPV205 belong to species γ-1. We then explored potential recombination events in genus γ-PV with the RDP4 program in a dataset of 74 viruses (71 HPV types with available full-length genomes and the 3 novel types). Two events, both located in the E1 ORF, met the inclusion criterion (p-values <0.05 with at least four methods) and persisted in different ORF combinations: an inter-species recombination in species γ-8 (major and minor parents: species γ-24 and γ-11, respectively), and an intra-species recombination in species γ-7 (recombinant strain: HPV170; major and minor parents: HPV-109 and HPV-149, respectively). These findings were confirmed by phylogenetic tree incongruence analysis. An additional incongruence was found in members of species γ-9 but it was not detected by the RDP4. This report expands our knowledge of the family Papillomaviridae and provides for the first time in silico evidence of recombination in genus γ-PV.

Oral papillomatosis caused by Enhydra lutris papillomavirus 1 (ElPV-1) in southern sea otters (Enhydra lutris nereis) in California, USA

The southern sea otter (Enhydra lutris nereis) is a threatened marine sentinel. During postmortem investigations of stranded sea otters from 2004 to 2013 in California, US, papillomas were detected in the oral cavity of at least seven otters via necropsy and histopathology. Next-generation sequencing of viral particles purified from a single papilloma revealed a novel papillomavirus, Enhydra lutris papillomavirus 1 (ElPV-1). The genome of ElPV-1 was obtained, representing the first fully sequenced viral genome from southern sea otters. Phylogenetic analysis of the entire L1 gene, as well as a concatenated protein identities plot of all papillomaviral genes revealed that ElPV-1 is a λ-papillomavirus, related to a raccoon papillomavirus (Procyon lotor papillomavirus type 1) and a canine oral papillomavirus. Immunohistochemical staining, using a cross-reactive bovine papillomavirus antibody, suggested that ElPV-1 is present in intranuclear inclusions and intracytoplasmic keratin granules. Virus-infected cells were scattered throughout the stratum granulosum and stratum spinosum of the gingival and buccal papillomas. Using ElPV-1-specific PCR, we confirmed viral DNA in oral papillomas from all seven stranded sea otters, with identical L1 sequences. This virus is associated with the development of oral papillomatosis in southern sea otters.

Papillomaviruses: Viral evolution, cancer and evolutionary medicine

Papillomaviruses (PVs) are a numerous family of small dsDNA viruses infecting virtually all mammals. PVs cause infections without triggering a strong immune response, and natural infection provides only limited protection against reinfection. Most PVs are part and parcel of the skin microbiota. In some cases, infections by certain PVs take diverse clinical presentations from highly productive self-limited warts to invasive cancers. We propose PVs as an excellent model system to study the evolutionary interactions between the immune system and pathogens causing chronic infections: genotypically, PVs are very diverse, with hundreds of different genotypes infecting skin and mucosa; phenotypically, they display extremely broad gradients and trade-offs between key phenotypic traits, namely productivity, immunogenicity, prevalence, oncogenicity and clinical presentation. Public health interventions have been launched to decrease the burden of PV-associated cancers, including massive vaccination against the most oncogenic human PVs, as well as systematic screening for PV chronic anogenital infections. Anti-PVs vaccines elicit protection against infection, induce cross-protection against closely related viruses and result in herd immunity. However, our knowledge on the ecological and intrapatient dynamics of PV infections remains fragmentary. We still need to understand how the novel anthropogenic selection pressures posed by vaccination and screening will affect viral circulation and epidemiology. We present here an overview of PV evolution and the connection between PV genotypes and the phenotypic, clinical manifestations of the diseases they cause. This differential link between viral evolution and the gradient cancer-warts-asymptomatic infections makes PVs a privileged playground for evolutionary medicine research.

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.

Deep sequencing extends the diversity of human papillomaviruses in human skin

Most viruses in human skin are known to be human papillomaviruses (HPVs). Previous sequencing of skin samples has identified 273 different cutaneous HPV types, including 47 previously unknown types. In the present study, we wished to extend prior studies using deeper sequencing. This deeper sequencing without prior PCR of a pool of 142 whole genome amplified skin lesions identified 23 known HPV types, 3 novel putative HPV types and 4 non-HPV viruses. The complete sequence was obtained for one of the known putative types and almost the complete sequence was obtained for one of the novel putative types. In addition, sequencing of amplimers from HPV consensus PCR of 326 skin lesions detected 385 different HPV types, including 226 previously unknown putative types. In conclusion, metagenomic deep sequencing of human skin samples identified no less than 396 different HPV types in human skin, out of which 229 putative HPV types were previously unknown.

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.

Diversity of human papillomaviruses in skin lesions

Pools of frozen biopsies from patients with squamous cell carcinoma (SCC) (n=29) actinic keratosis (AK) (n=31), keratoacanthoma (n=91) and swab samples from 84 SCCs and 91 AKs were analysed with an extended HPV general primer PCR and high-throughput sequencing of amplimers. We found 273 different HPV isolates (87 known HPV types, 139 previously known HPV sequences (putative types) and 47 sequences from novel putative HPV types). Among the new sequences, five clustered in genus Betapapillomavirus and 42 in genus Gammapapillomavirus. Resequencing of the three pools between 21 to 70 times resulted in the detection of 283 different known or putative HPV types, with 156 different sequences found in only one of the pools. Type-specific PCRs for 37 putative types from an additional 296 patients found only two of these putative types. In conclusion, skin lesions contain a large diversity of HPV types, but most appeared to be rare infections.

Animal papillomaviruses

We provide an overview of the host range, taxonomic classification and genomic diversity of animal papillomaviruses. The complete genomes of 112 non-human papillomavirus types, recovered from 54 different host species, are currently available in GenBank. The recent characterizations of reptilian papillomaviruses extend the host range of the Papillomaviridae to include all amniotes. Although the genetically diverse papillomaviruses have a highly conserved genomic lay-out, deviations from this prototypic genome organization are observed in several animal papillomaviruses, and only the core ORFs E1, E2, L2 and L1 are present in all characterized papillomavirus genomes. The discovery of papilloma-polyoma hybrids BPCV1 and BPCV2, containing a papillomaviral late region but an early region encoding typical polyomaviral nonstructural proteins, and the detection of recombination breakpoints between the early and late coding regions of cetacean papillomaviruses, could indicate that early and late gene cassettes of papillomaviruses are relatively independent entities that can be interchanged by recombination.

Identification of human papillomavirus type 156, the prototype of a new human gammapapillomavirus species, by a generic and highly sensitive PCR strategy for long DNA fragments

This study developed a hanging-droplet long PCR, a generic and highly sensitive strategy to facilitate the identification of new human papillomavirus (HPV) genomes. This novel procedure used for the first time the hanging-droplet PCR technique for the amplification of long DNA fragments with generic primers targeting the L1 and E1 regions. It was first applied to the amplification of types belonging to the highly divergent genus Gammapapillovirus (γ-PV). The hanging-droplet long PCR was 100-fold more sensitive than a simple long PCR procedure, detecting as few as ten copies of HPV-4. Nineteen skin samples, potentially containing putative HPV types from the γ-PV genus, were also screened. The method identified four γ-PV genomic halves from new and previously described putative types, and made the full characterization of HPV-156 possible. This novel virus meets the criteria for a new species within the γ-PV genus, with nucleotide identities in the L1 ORF ranging from 58.3 to 67.3 % compared with representative types of the current γ-PV species. HPV-156 showed the highest identity to HPV-60 (67.3 %) from species γ-4, and was consistently closely related to it in both late- and early-gene-derived phylogenies. In conclusion, this report provides a versatile and highly sensitive approach that allowed identification of the prototype of a new species within the γ-PV genus. Its application with primers targeting the different genera in which both human and non-human PVs are distributed may facilitate characterization of the missing members of the family Papillomaviridae.

Recombined sequences between the non-coding control regions of JC and BK viruses found in the urine of a renal transplantation patient

Kidney cells are the common host for JC virus (JCV) and BK virus (BKV). Reactivation of JCV and/or BKV in patients after organ transplantation, such as renal transplantation, may cause hemorrhagic cystitis and polyomavirus-associated nephropathy. Furthermore, JCV and BKV may be shed in the urine after reactivation in the kidney. Rearranged as well as archetypal non-coding control regions (NCCRs) of JCV and BKV have been frequently identified in human samples. In this study, three JC/BK recombined NCCR sequences were identified in the urine of a patient who had undergone renal transplantation. They were designated as JC-BK hybrids 1, 2, and 3. The three JC/BK recombinant NCCRs contain up-stream JCV as well as down-stream BKV sequences. Deletions of both JCV and BKV sequences were found in these recombined NCCRs. Recombination of DNA sequences between JCV and BKV may occur during co-infection due to the relatively high homology of the two viral genomes.

Novel animal papillomavirus sequences and accurate phylogenetic placement

All amniotes are probably infected by specific papillomaviruses (PVs), but knowledge about PV diversity remains sparse. An insufficient taxon sampling, and a focus on humans as hosts, may perturb phylogenetic analyses leading to wrong conclusions about PV evolution. We performed a systematic approach to explore the diversity of PVs combining rolling circle amplification with the use of "universal" primers to search for the presence of novel PV sequences in animal samples. We communicate 12 sequences putatively corresponding to novel PVs gained from 10 host species in eight mammal families: Bovidae, Canidae, Cervidae, Equidae, Hominidae, Phocoenidae, Procyonidae and Pteropodidae. The phylogenetic position of the new sequences was inferred with an evolutionary placement algorithm under a Maximum Likelihood framework using a pre-computed, well-resolved tree constructed with the E1-E2-L1 gene sequences as a backbone. The new sequences were phylogenetically diverse and could be respectively placed with confidence within all four PV crown groups. The prevailing presence of sequences from the crown groups Alpha+Omikron-PVs and Beta+Xi-PVs may correspond to an increased viral diversity in these taxa, or rather reflect a combination of anthropocentric bias and preferential amplification from commonly used "universal" primers. Our results combined with literature data support the view that the number and diversity of animal PVs is overwhelmingly large.

Papillomaviruses and herpesviruses: who is who in genital tumor development of free-ranging Atlantic bottlenose dolphins (Tursiops truncatus)?

The number of studies addressing neoplasia in marine mammals has recently increased, giving rise to concern whether such lesions could be reflective of an emerging infectious disease. Eight species-specific viruses, seven papillomaviruses (PVs) and two herpesviruses (HVs) have separately been shown to be associated with genital tumors in Atlantic bottlenose dolphins (Tursiops truncatus, Tt): TtPV1-6, as well as HVs provisionally assigned the names DeHV4 and -5 (Delphinid HVs). A definite causal role of these viruses in cell transformation remains to be demonstrated. Concurrent PV- and HV-infection has never been reported in marine mammals. DNA extractions from biopsies of genital tumors derived from 15 free-ranging Atlantic bottlenose dolphins were selected for molecular examination. Polymerase chain reaction (PCR) analyses revealed the presence of DeHV4, while a serological screening using an antibody-based TtPV enzyme-linked immunosorbent assay (ELISA) demonstrated previous and/or current infection of the HV-positive dolphins with at least one TtPV type. Therefore, care must be taken when drawing conclusions about viral causalities in tumor development, since the "hit and run" and other mechanisms have been described for types of both viral families. This study presents the first evidence of marine mammals having a history of PV- as well as HV-infection and discusses the disputed effects of viral co-infection.