Why HPV16? Why, now, HPV42? How the discovery of HPV42 in rare cancers provides an opportunity to challenge our understanding about the transition between health and disease for common members of the healthy microbiota

In 2022, a bioinformatic, agnostic approach identified HPV42 as causative agent of a rare cancer, later confirmed experimentally. This unexpected association offers an opportunity to reconsider our understanding about papillomavirus infections and cancers. We have expanded our knowledge about the diversity of papillomaviruses and the diseases they cause. Yet, we still lack answers to fundamental questions, such as what makes HPV16 different from the closely related HPV31 or HPV33; or why the very divergent HPV13 and HPV32 cause focal epithelial hyperplasia, while HPV6 or HPV42 do not, despite their evolutionary relatedness. Certain members of the healthy skin microbiota are associated to rare clinical conditions. We propose that a focus on cellular phenotypes, most often transient and influenced by intrinsic and extrinsic factors, may help understand the continuum between health and disease. A conceptual switch is required towards an interpretation of biology as a diversity of states connected by transition probabilities, rather than quasi-deterministic programs. Under this perspective, papillomaviruses may only trigger malignant transformation when specific viral genotypes interact with precise cellular states. Drawing on Canguilhem's concepts of normal and pathological, we suggest that understanding the transition between fluid cellular states can illuminate how commensal-like infections transition from benign to malignant.

Complete genomic characterization of bovine papillomavirus type 1 and 2 strains infers ongoing cross-species transmission between cattle and horses

Infection with bovine papillomavirus (BPV) types 1 and 2 results in the most common skin tumor of horses, termed equine sarcoid. The persistent and recurrent nature of this tumor stands in contrast to the regressive nature of BPV-1/- 2 induced cutaneous papillomas in cattle. The circulation of horse-specific BPV-1/- 2 variants within equine populations has been suggested as a possible explanation for the difference in clinical presentation of BPV-1/- 2 infection between horses and cattle. In order to investigate this hypothesis, we identified 98 complete BPV-1/- 2 genomes using a Nanopore sequencing approach. Separate BPV-1/- 2 alignments were used to infer Bayesian phylogenetic trees. Phylogeny-trait association concerning host species was investigated using Bayesian Tip-association Significance software (BaTS) Overall, 179 unique BPV-1 and 128 BPV-2 substitutions were found. The E2 coding region in the viral genome exhibited an exceptionally high rate of non-synonymous mutations (81 %, n = 13/16). Interestingly, extensive deletions in the L1/L2 region (up to 1.5 kb) were found exclusively in horse-derived samples. Nevertheless, the most frequently detected single nucleotide polymorphisms were shared between equine and bovine hosts, which is in agreement with BaTS results indicating no phylogeny-host correlation. We found indications that horse-specific mutations might exist in subpopulations of equine derived BPV-1/- 2, but these did not result in horse-adapted genetic variants. Based on these observations, cross-species transmission from cattle to horses seems to be an ongoing process, rather than an ancient occurrence that has been followed by the circulation of horse-adapted BPV variants in the horse population..

Co-infection of distinct papillomavirus types in a captive North American porcupine

BACKGROUND

Only two cases of papillomavirus infections in North American porcupines (Erethizon dorsatum) have been described thus far, and molecular investigation linked these cases to two distinct papillomavirus species.

METHODS

In this report, we present the clinical, histological and molecular investigation of a third case of a porcupine papillomavirus infection. Papillomatous lesions occurred on the upper and lower lip of an otherwise healthy three-year old female that was kept in captivity. Within one month, the lesions progressed into exophytic black nodules, followed by a temporary stabilization and ultimately spontaneous regression within seven months of their initial observation. PCR-based screening using specific primers for Erethizon dorsatum papillomavirus 1 and 2 revealed the presence of both these virus types, after which nanopore sequencing was used to determine the complete sequences of the two virus genomes.

RESULTS

One of the genomes shares 99.9% similarity with the only known sequence for Erethizon dorsatum papillomavirus 1, while the second represents a distinct lineage of Erethizon dorsatum papillomavirus 2, sharing only 93.3% similarity with the previously discovered strain.

CONCLUSIONS

This report marks the first observation of a papillomavirus co-infection in a North American porcupine, although the individual contribution of the two virus types to the clinical presentation was not assessed.

Synonymous nucleotide changes drive papillomavirus evolution

Papillomaviruses have been evolving alongside their hosts for at least 450 million years. This review will discuss some of the insights gained into the evolution of this diverse family of viruses. Papillomavirus evolution is constrained by pervasive purifying selection to maximize viral fitness. Yet these viruses need to adapt to changes in their environment, e.g., the host immune system. It has long been known that these viruses evolved a codon usage that doesn't match the infected host. Here we discuss how papillomavirus genomes evolve by acquiring synonymous changes that allow the virus to avoid detection by the host innate immune system without changing the encoded proteins and associated fitness loss. We discuss the implications of studying viral evolution, lifecycle, and cancer progression.

Farklı işletme koşullarında ineklerin meme başlarında bovine papillomavirus kaynaklı gelişen siğillerin tedavisinde uygulanan yöntemlerin (otolog aşı ve newcastle disease virus La sota suşu) değerlendirilmesi

Burdur-Merkez ve ilçelerinde halk elinde, farklı yapıdaki işletmelerde bulunan 106 adet dişi ineğin bovine papillomavirus (BPV) kaynaklı gelişen meme siğillerinden doku örnekleri toplandı. BPV tip spesifik primerleri açısından PCR testi ile 106 örneğin hepsinde BPV pozitiflik belirlendi. Çalışmada; işletme tipleri, bireysel hayvan bakım ve besleme koşulları, sağım koşulları ve ahır temizliği parametreler yönünden incelenerek değerlendirildi. Ayrıca, meme siğillerinin makroskobik görünümlerine göre (düz ve yuvarlak, pirinç tanesi, saplı-sivri) sınıflandırılmış, BPV tip spesifik primerleri açısından pozitiflik tespit edilmiş 30 örnekte iki farklı tedavi uygulama yöntemleri çalışıldı. Birinci grupta tasnif edilmiş 15 hayvana otojen aşı, ikinci grupta tasnif edilmiş 15 hayvana Newcastle disease virus (NDV) La Sota suş aşı uygulamaları yapıldı. Otojen aşı uygulamalarının yapıldığı, farklı makroskobik görünümlere sahip meme siğillerinde %100 gerileme ve iyileşme sağlandı. Newcastle disease virus La Sota aşı uygulamalarının yapıldığı farklı makroskobik görünümlere sahip meme siğillerinde %53,3 gerileme ve iyileşme sağlandı. İneklerde meme siğillerinin tedavisinde otojen aşı uygulamalarının daha etkin sonuçlar verdiği tespit edildi. Ayrıca, çalışmada meme siğillerinin oluşumunda; işletme, bakım, besleme, sağım, ahır temizliği ve diğer faktörlerin de etkili olabileceği yorumlandı.

Divergent Viruses Discovered in Swine Alter the Understanding of Evolutionary History and Genetic Diversity of the Respirovirus Genus and Related Porcine Parainfluenza Viruses

Paramyxoviridae is a rapidly growing family of viruses, whose potential for cross-species transmission makes it difficult to predict the harm of newly emerging viruses to humans and animals. To better understand their diversity, evolutionary history, and co-evolution with their hosts, we analyzed a collection of porcine parainfluenza virus (PPIV) genomes to reconstruct the species classification basis and evolutionary history of the Respirovirus genus. We sequenced 17 complete genomes of porcine respirovirus 1 (also known as porcine parainfluenza virus 1; PPIV-1), thereby nearly tripling the number of currently available PPIV-1 genomes. We found that PPIV-1 was widely prevalent in China with two divergent lineages, PPIV-1a and PPIV-1b. We further provided evidence that a new species, porcine parainfluenza virus 2 (PPIV-2), had recently emerged in China. Our results pointed to a need for revising the current species demarcation criteria of the Respirovirus genus. In addition, we used PPIV-1 as an example to explore recombination and diversity of the Respirovirus genus. Interestingly, we only detected heterosubtypic recombination events between PPIV-1a and PPIV-1b with no intrasubtypic recombination events. The recombination hotspots highlighted a diverse geography-dependent genome structure of paramyxovirus infecting swine in China. Furthermore, we found no evidence of co-evolution between respirovirus and its host, indicating frequent cross-species transmission. In summary, our analyses showed that swine can be infected with a broad range of respiroviruses and recombination may serve as an important evolutionary mechanism for the Respirovirus genus’ greater diversity in genome structure than previously anticipated.

IMPORTANCE Livestock have emerged as critically underrecognized sources of paramyxovirus diversity, including pigs serving as the source of Nipah virus (NiV) and swine parainfluenza virus type 3, and goats and bovines harboring highly divergent viral lineages. Here, we identified a new species of Respirovirus genus named PPIV-2 in swine and proposed to revise the species demarcation criteria of the Respirovirus genus. We found heterosubtypic recombination events and high genetic diversity in PPIV-1. Further, we showed that genetic recombination may have occurred in the Respirovirus genus which may be associated with host range expansion. The continued expansion of Respirovirus genus diversity in livestock with relatively high human contact rates requires enhanced surveillance and ongoing evaluation of emerging cross-species transmission threats.

Identification and evolutionary analysis of papillomavirus sequences in New World monkeys (genera Sapajus and Alouatta) from Argentina

OBJECTIVE

In this study, we investigated the occurrence of papillomavirus (PV) infection in non-human primates (NHPs) in northeastern Argentina. We also explored their evolutionary history and evaluated the co-speciation hypothesis in the context of primate evolution.

METHODS

We obtained DNA samples from 57 individuals belonging to wild and captive populations of Alouatta caraya, Sapajus nigritus, and Sapajus cay. We assessed PV infection by PCR amplification with the CUT primer system and sequencing of 337 bp (112 amino acids) of the L1 gene. The viral sequences were analyzed by phylogenetic and Bayesian coalescence methods to estimate the time to the most common recent ancestor (t_MRCA) using BEAST, v1.4.8 software. We evaluated viral/host tree congruence with TreeMap v3.0.

RESULTS

We identified two novel putative PV sequences of the genus Gammapapillomavirus in Sapajus spp. and Alouatta caraya (SPV1 and AcPV1, respectively). The t_MRCA of SPV1 was estimated to be 11,941,682 years before present (ybp), and that of AcPV1 was 46,638,071 ybp, both before the coalescence times of their hosts (6.4 million years ago [MYA] and 6.8 MYA, respectively). Based on the comparison of primate and viral phylogenies, we found that the PV tree was no more congruent with the host tree than a random tree would be (P > 0.05), thus allowing us to reject the model of virus-host coevolution.

CONCLUSION

This study presents the first evidence of PV infection in platyrrhine species from Argentina, expands the range of described hosts for these viruses, and suggests new scenarios for their origin and dispersal.

Host diversification is concurrent with linear motif evolution in a Mastadenovirus hub protein

Over one hundred Mastadenovirus types infect seven orders of mammals. Virus-host coevolution may involve cospeciation, duplication, host switch and partial extinction events. We reconstruct Mastadenovirus diversification, finding that while cospeciation is dominant, the other three events are also common in Mastadenovirus evolution. Linear motifs are fast-evolving protein functional elements and key mediators of virus-host interactions, thus likely to partake in adaptive viral evolution. We study the evolution of eleven linear motifs in the Mastadenovirus E1A protein, a hub of virus-host protein-protein interactions, in the context of host diversification. The reconstruction of linear motif gain and loss events shows fast linear motif turnover, corresponding a virus-host protein-protein interaction turnover orders of magnitude faster than in model host proteomes. Evolution of E1A linear motifs is coupled, indicating functional coordination at the protein scale, yet presents motif-specific patterns suggestive of convergent evolution. We report a pervasive association between Mastadenovirus host diversification events and the evolution of E1A linear motifs. Eight of 17 host switches associate with the gain of one linear motif and the loss of four different linear motifs, while five of nine partial extinctions associate with the loss of one linear motif. The specific changes in E1A linear motifs during a host switch or a partial extinction suggest that changes in the host molecular environment lead to modulation of the interactions with the retinoblastoma protein and host transcriptional regulators. Altogether, changes in the linear motif repertoire of a viral hub protein are associated with adaptive evolution events during Mastadenovirus evolution.

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.

Ancient Evolutionary History of Human Papillomavirus Type 16, 18 and 58 Variants Prevalent Exclusively in Japan

Human papillomavirus (HPV) is a sexually transmitted virus with an approximately 8-kilo base DNA genome, which establishes long-term persistent infection in anogenital tissues. High levels of genetic variations, including viral genotypes and intra-type variants, have been described for HPV genomes, together with geographical differences in the distribution of genotypes and variants. Here, by employing a maximum likelihood method, we performed phylogenetic analyses of the complete genome sequences of HPV16, HPV18 and HPV58 available from GenBank (n = 627, 146 and 157, respectively). We found several characteristic clusters that exclusively contain HPV genomes from Japan: two for HPV16 (sublineages A4 and A5), one for HPV18 (sublineage A1) and two for HPV58 (sublineages A1 and A2). Bayesian phylogenetic analyses of concatenated viral gene sequences showed that divergence of the most recent common ancestor of these Japan-specific clades was estimated to have occurred ~98,000 years before present (YBP) for HPV16 A4, ~39,000 YBP for HPV16 A5, ~38,000 YBP for HPV18 A1, ~26,000 for HPV58 A1 and ~25,000 YBP for HPV58 A2. This estimated timeframe for the divergence of the Japan-specific clades suggests that the introduction of these HPV variants into the Japanese archipelago dates back to at least ~25,000 YBP and provides a scenario of virus co-migration with ancestral Japanese populations from continental Asia during the Upper Paleolithic period.

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.

The Reservoir of Persistent Human Papillomavirus Infection; Strategies for Elimination Using Anti-Viral Therapies

Human Papillomaviruses have co-evolved with their human host, with each of the over 200 known HPV types infecting distinct epithelial niches to cause diverse disease pathologies. Despite the success of prophylactic vaccines in preventing high-risk HPV infection, the development of HPV anti-viral therapies has been hampered by the lack of enzymatic viral functions, and by difficulties in translating the results of in vitro experiments into clinically useful treatment regimes. In this review, we discuss recent advances in anti-HPV drug development, and highlight the importance of understanding persistent HPV infections for future anti-viral design. In the infected epithelial basal layer, HPV genomes are maintained at a very low copy number, with only limited viral gene expression; factors which allow them to hide from the host immune system. However, HPV gene expression confers an elevated proliferative potential, a delayed commitment to differentiation, and preferential persistence of the infected cell in the epithelial basal layer, when compared to their uninfected neighbours. To a large extent, this is driven by the viral E6 protein, which functions in the HPV life cycle as a modulator of epithelial homeostasis. By targeting HPV gene products involved in the maintenance of the viral reservoir, there appears to be new opportunities for the control or elimination of chronic HPV infections.

K-Mer Analyses Reveal Different Evolutionary Histories of Alpha, Beta, and Gamma Papillomaviruses

Papillomaviruses (PVs) are a heterogeneous group of DNA viruses that can infect fish, birds, reptiles, and mammals. PVs infecting humans (HPVs) phylogenetically cluster into five genera (Alpha-, Beta-, Gamma-, Mu- and Nu-PV), with differences in tissue tropism and carcinogenicity. The evolutionary features associated with the divergence of Papillomaviridae are not well understood. Using a combination of k-mer distributions, genetic metrics, and phylogenetic algorithms, we sought to evaluate the characteristics and differences of Alpha-, Beta- and Gamma-PVs constituting the majority of HPV genomes. A total of 640 PVs including 442 HPV types, 27 non-human primate PV types, and 171 non-primate animal PV types were evaluated. Our analyses revealed the highest genetic diversity amongst Gamma-PVs compared to the Alpha and Beta PVs, suggesting reduced selective pressures on Gamma-PVs. Using a sequence alignment-free trimer (k = 3) phylogeny algorithm, we reconstructed a phylogeny that grouped most HPV types into a monophyletic clade that was further split into three branches similar to alignment-based classifications. Interestingly, a subset of low-risk Alpha HPVs (the species Alpha-2, 3, 4, and 14) split from other HPVs and were clustered with non-human primate PVs. Surprisingly, the trimer-constructed phylogeny grouped the Gamma-6 species types originally isolated from the cervicovaginal region with the main Alpha-HPV clade. These data indicate that characterization of papillomavirus heterogeneity via orthogonal approaches reveals novel insights into the biological understanding of HPV genomes.

Computational Modeling of Immune Response Triggering Immunogenic Peptide Vaccine Against the Human Papillomaviruses to Induce Immunity Against Cervical Cancer

Papillomaviruses are placed within the family Papillomaviride, and the members of this family have a double-stranded circular DNA genome. Every year, ∼30% of cancers are reported to be human papillomavirus (HPV) related, which represents 63,000 cancers of all infectious agent-induced cancers. HPV16 and HPV18 are reported to be associated with 70% of cervical cancers. The quest for an effective drug or vaccine candidate still continues. In this study, we aim to design B cell and T cell epitope-based vaccine using the two structural major capsid protein L1 and L2 as well as other three important proteins (E1, E2, and E6) against HPV strain 16 (HPV16). We used a computational pipeline to design a multiepitope subunit vaccine and tested its efficacy using in silico computational modeling approaches. Our analysis revealed that the multiepitope subunit vaccine possesses antigenic properties, and using in silico cloning method revealed proper expression and downstream processing of the vaccine construct. Besides this, we also performed in silico immune simulation to check the immune response upon the injection. Our results strongly suggest that this vaccine candidate should be tested immediately for the immune response against the cervical cancer-causing agent. The safety, efficacy, expression, and immune response profiling makes it the first choice for experimental and in vivo setup.

Identification of a Novel Papillomavirus Type (MfoiPV1) Associated with Acrochordon in a Stone Marten (Martes foina)

Papillomaviruses (PVs) are an extremely large group of viruses that cause skin and mucosal infections in humans and various domestic and wild animals. Nevertheless, there is limited knowledge about PVs in wildlife hosts, including mustelid species. This study describes a case in stone marten (Martes foina) with a clinical manifestation of skin tumor, which is rather atypical for infections with PVs. The result of the papillomavirus PCR performed on the skin tumor sample was positive, and the complete PV genome was determined in the studied sample using next-generation sequencing technology. The analysis of the PV genome revealed infection of the stone marten with a putative new PV type belonging to the Dyonupapillomavirus genus. The proposed new stone marten PV type was named MfoiPV1.

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.

Primate phageomes are structured by superhost phylogeny and environment

Humans harbor diverse communities of microorganisms, the majority of which are bacteria in the gastrointestinal tract. These gut bacterial communities in turn host diverse bacteriophage (hereafter phage) communities that have a major impact on their structure, function, and, ultimately, human health. However, the evolutionary and ecological origins of these human-associated phage communities are poorly understood. To address this question, we examined fecal phageomes of 23 wild nonhuman primate taxa, including multiple representatives of all the major primate radiations. We find relatives of the majority of human-associated phages in wild primates. Primate taxa have distinct phageome compositions that exhibit a clear phylosymbiotic signal, and phage-superhost codivergence is often detected for individual phages. Within species, neighboring social groups harbor compositionally and evolutionarily distinct phageomes, which are structured by superhost social behavior. Captive nonhuman primate phageome composition is intermediate between that of their wild counterparts and humans. Phage phylogenies reveal replacement of wild great ape-associated phages with human-associated ones in captivity and, surprisingly, show no signal for the persistence of wild-associated phages in captivity. Together, our results suggest that potentially labile primate-phage associations have persisted across millions of years of evolution. Across primates, these phylosymbiotic and sometimes codiverging phage communities are shaped by transmission between groupmates through grooming and are dramatically modified when primates are moved into captivity.

A Comprehensive Study of Cutaneous Fibropapillomatosis in Free-Ranging Roe Deer (Capreolus capreolus) and Red Deer (Cervus elaphus): from Clinical Manifestations to Whole-Genome Sequencing of Papillomaviruses

Papillomaviruses (PVs) are an extremely large group of viruses that cause skin and mucosa infections in humans and various animals. In roe deer and red deer, most PVs belong to the Deltapapillomavirus genus and cause neoplastic changes that are generally described as fibropapillomas. Despite the wide distribution of roe and red deer throughout Europe and beyond, the data in the scientific literature regarding the widespread distribution of PVs and the genetic variability of PV genomes in these species are rather scarce. This study describes cutaneous fibropapillomatosis cases in roe and red deer with clinical manifestations that are typical of infections with PVs. In all cases, the presence of PV DNA was confirmed using PCR, followed by Sanger sequencing of the partial L1 gene. The complete PV genomes were determined in all the investigated samples using next-generation sequencing technology, revealing infections of roe deer with the CcaPV1-type and red deer with the CePV1v-type variant. A comparison of the complete CcaPV1-type and CePV1v-type variant genome sequences reported here with already available complete genome sequences in GenBank revealed their great genetic stability across time and space.

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

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

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.

PVAmpliconFinder: a workflow for the identification of human papillomaviruses from high-throughput amplicon sequencing

BACKGROUND

The detection of known human papillomaviruses (PVs) from targeted wet-lab approaches has traditionally used PCR-based methods coupled with Sanger sequencing. With the introduction of next-generation sequencing (NGS), these approaches can be revisited to integrate the sequencing power of NGS. Although computational tools have been developed for metagenomic approaches to search for known or novel viruses in NGS data, no appropriate tool is available for the classification and identification of novel viral sequences from data produced by amplicon-based methods.

RESULTS

We have developed PVAmpliconFinder, a data analysis workflow designed to rapidly identify and classify known and potentially new Papillomaviridae sequences from NGS amplicon sequencing with degenerate PV primers. Here, we describe the features of PVAmpliconFinder and its implementation using biological data obtained from amplicon sequencing of human skin swab specimens and oral rinses from healthy individuals.

CONCLUSIONS

PVAmpliconFinder identified putative new HPV sequences, including one that was validated by wet-lab experiments. PVAmpliconFinder can be easily modified and applied to other viral families. PVAmpliconFinder addresses a gap by providing a solution for the analysis of NGS amplicon sequencing, increasingly used in clinical research. The PVAmpliconFinder workflow, along with its source code, is freely available on the GitHub platform: https://github.com/IARCbioinfo/PVAmpliconFinder.

Risk of cervical lesions in high-risk HPV positive women with normal cytology: a retrospective single-center study in China

Background

To investigate high-risk HPV (hr-HPV) genotype distributions and the association between hr-HPV infection with severity of the cervical lesions in women with normal cytology.

Methods

In this cross-sectional study, the result of the hr-HPV test and biopsy of colposcopy of women with normal cytology from January 2012 to January 2019 were analyzed. The detection rate of high-grade squamous intraepithelial lesion (HSIL) and cervical cancer were calculated among different hr-HPV genotypes, viral load group, and age groups.

Results

Five thousand eight hundred eighty women were enrolled in this study. Overall, 59.97% had normal histological results, 19.32% had HSIL, and 1.07% had cervical cancer. The detection rate of HSIL or worse (HSIL+) in women with single HPV16(34.00%), HPV31(27.50%), HPV33(25.58%), and HPV52(20.88%) infection were higher significantly than single HPV18 (15.59%) infection, respectively. The HSIL+ detection rate between HPV16 single infection and multiple infections (excluding HPV18) was no significant difference (34% vs 35.47%, P = 0.638), contrary to HPV18(12.59% vs 21.67%, P = 0.022). In women without HPV16/18 infections, HSIL+ detection rates for single, double, and triple or more hr-HPV infections were 12.28, 20.31, and 37.50%, the risk of detection of HSIL+ significantly increasing. With the hr-HPV DNA load increases, the risk of detection of HSIL+ (χ2 = 91.01, P < 0.0001) and invasive cervical cancer (χ2 = 5.757, P = 0.016) increase. In age < 30, 31–40, 41–50, 51–60, > 60 group, HSIL+ detection rate were 24.80%、22.10%、19.59%、14.29, and 12.61%, respectively.

Conclusion

Women who have normal cytology with HPV 16/18/31/33/52/58 infections, multiple HPV infections and high viral load, have a higher detection rate of HSIL+.

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.

Leopardus wiedii Papillomavirus type 1, a novel papillomavirus species in the tree ocelot, suggests Felidae Lambdapapillomavirus polyphyletic origin and host-independent evolution

The limited knowledge on Papillomavirus diversity (particularly in wild animal species) influences the accuracy of PVs phylogeny and their evolutionary history, and hinders the comprehension of PVs pathogenicity, especially the mechanism of virus - related cancer progression. This study reports the identification of Leopardus wiedii Papillomavirus type 1 (LwiePV1), the first PV type within Lambdapapillomavirus in a Leopardus host. LwiePV1 full genome sequencing allowed the investigation of its taxonomic position and phylogeny. Based on results, LwiePV1 should be assigned to a novel PV species providing evidence for a polyphyletic origin of feline lambda PVs, and representing an exception to codivergence between feline lambda PVs and their hosts. Results improve our knowledge on PV diversity and pave the way to future studies investigating biological and evolutionary features of animal PVs.

New insights into Sauropsid Papillomaviridae evolution and epizootiology: discovery of two novel papillomaviruses in native and invasive Island geckos

Papillomaviruses cause persistent infections in skin and mucosal membranes and, in at least one species, are also be able to infect a tissue of mesenchymal origin. Infections may either be subclinical or induce proliferative lesions. Of the known papillomaviruses, the majority that have been characterized are from humans and other mammals. Currently, only fifteen complete bird and reptile papillomavirus genomes have been described, and they have been found in birds (n = 11), turtles (n = 2), and snakes (n = 2). Using next-generation sequencing technologies and virus-specific PCR, we have identified two novel papillomavirus genomes, Hemidactylus frenatus Papillomavirus 1 and 2 (HfrePV1, HfrePV2), in the widely distributed and highly invasive Asian house gecko (H.frenatus) and mute gecko (Gehyra mutilata) on Christmas Island and Cocos (Keeling) Islands. HfrePV1 was also detected in critically endangered Lister’s geckos (Lepidodactylus listeri) in their captive breeding colony on Christmas Island. Tissue-containing virus included epidermis, oral mucosa, and liver (HfrePV1) and epidermis, liver, and colon (HfrePV2). Concurrent infections were found in both H.frenatus and G.mutilata. Invasive mourning geckos (Lepidodactylus lugubris) (n = 4), Sri Lankan house geckos (Hemidactylus parvimaculatus) (n = 3), flat-tailed house geckos (Hemidactylus platyurus) (n = 4) from the Cocos Islands, and blue-tailed skinks (Cryptoblepharus egeriae) (n = 10) from Christmas Island were also screened but were not found to be infected. The novel HfrePV1 and HfrePV2 genomes were 7,378 bp and 7,380 bp in length, respectively, and each contained the early (E1, E2, and E7), and late (L1 and L2) open-reading frames. Phylogenetic analysis of the concatenated E1, E2, and L1 proteins from both papillomaviruses revealed that they clustered with, but were basal to, the Sauropsida clade containing bird and reptile viruses. This study sheds light on the evolution of papillomaviruses and the distribution of pathogens in a highly invasive species impacting endangered populations of geckos.

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.

Novel Polyomaviruses in Mammals from Multiple Orders and Reassessment of Polyomavirus Evolution and Taxonomy

As the phylogenetic organization of mammalian polyomaviruses is complex and currently incompletely resolved, we aimed at a deeper insight into their evolution by identifying polyomaviruses in host orders and families that have either rarely or not been studied. Sixteen unknown and two known polyomaviruses were identified in animals that belong to 5 orders, 16 genera, and 16 species. From 11 novel polyomaviruses, full genomes could be determined. Splice sites were predicted for large and small T antigen (LTAg, STAg) coding sequences (CDS) and examined experimentally in transfected cell culture. In addition, splice sites of seven published polyomaviruses were analyzed. Based on these data, LTAg and STAg annotations were corrected for 10/86 and 74/86 published polyomaviruses, respectively. For 25 polyomaviruses, a spliced middle T CDS was observed or predicted. Splice sites that likely indicate expression of additional, alternative T antigens, were experimentally detected for six polyomaviruses. In contrast to all other mammalian polyomaviruses, three closely related cetartiodactyl polyomaviruses display two introns within their LTAg CDS. In addition, the VP2 of Glis glis (edible dormouse) polyomavirus 1 was observed to be encoded by a spliced transcript, a unique experimental finding within the Polyomaviridae family. Co-phylogenetic analyses based on LTAg CDS revealed a measurable signal of codivergence when considering all mammalian polyomaviruses, most likely driven by relatively recent codivergence events. Lineage duplication was the only other process whose influence on polyomavirus evolution was unambiguous. Finally, our analyses suggest that an update of the taxonomy of the family is required, including the creation of novel genera of mammalian and non-mammalian polyomaviruses.

Evidence for swine and human papillomavirus in pig slurry in Italy

AIMS

The diversity and the geographical distribution of swine papillomaviruses (PVs) are virtually unknown. The occurrence and the diversity of swine PV were therefore investigated in pig slurry collected in Italy, to contribute towards filling this gap in knowledge.

METHODS AND RESULTS

Twenty-two slurry samples underwent analysis by nested PCR and DNA sequencing using published and newly designed specific primer pairs for Sus scrofa papillomavirus (SsPV) type 1 and 2 (SsPV1 and 2), along with degenerate PV-specific primers targeting the major coat protein L1 and the helicase protein E1. Overall, three samples (13·6%) were positive for SsPV1 by specific primers, and nucleotide (nt) sequences showed 99-100% nt identity with SsPV1 variant a (EF395818), while SsPV2 was not found in any sample. Using generic primers, eight samples (36·4%) were tested positive for human papillomavirus (HPV), and were characterized as follows: β1-HPV8, β1-HPV14, β1-HPV206, β2-HPV113, β2-HPV120 and γ1-HPV173. Moreover, one unclassified γ-type was detected.

CONCLUSIONS

Both swine and human PVs were detected in pig slurry in this study. The unexpected presence of HPV in pig waste could be explained as the result of an improper use of the sewage collection pits and/or with improper procedures of the operators.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study reports the first detection of SsPV1 in Italy, along with the first detection of HPVs in pig slurry samples in Italy, and expands our knowledge about PV diversity and geographic distribution.

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

Mucosal and Cutaneous Human Papillomavirus Infections and Cancer Biology

Papillomaviridae is a family of small non-enveloped icosahedral viruses with double-stranded circular DNA. More than 200 different human papillomaviruses (HPVs) have been listed so far. Based on epidemiological data, a subgroup of alphapapillomaviruses (alpha HPVs) was referred to as high-risk (HR) HPV types. HR HPVs are the etiological agents of anogenital cancer and a subset of head and neck cancers. The cutaneous HPV types, mainly from beta and gamma genera, are widely present on the surface of the skin in the general population. However, there is growing evidence of an etiological role of betapapillomaviruses (beta HPVs) in non-melanoma skin cancer (NMSC), together with ultraviolet (UV) radiation. Studies performed on mucosal HR HPV types, such as 16 and 18, showed that both oncoproteins E6 and E7 play a key role in cervical cancer by altering pathways involved in the host immune response to establish a persistent infection and by promoting cellular transformation. Continuous expression of E6 and E7 of mucosal HR HPV types is essential to initiate and to maintain the cellular transformation process, whereas expression of E6 and E7 of cutaneous HPV types is not required for the maintenance of the skin cancer phenotype. Beta HPV types appear to play a role in the initiation of skin carcinogenesis, by exacerbating the accumulation of UV radiation-induced DNA breaks and somatic mutations (the hit-and-run mechanism), and they would therefore act as facilitators rather than direct actors in NMSC. In this review, the natural history of HPV infection and the transforming properties of various HPV genera will be described, with a particular focus on describing the state of knowledge about the role of cutaneous HPV types in NMSC.

Adaptation of Alpha-Papillomavirus over Millennia

Papillomaviruses (PVs) are a group of small DNA viruses that, with around 350 million years of evolution, acquired the capacity of infecting a broad range of vertebrates, including humans. To date, more than 300 PV types have been isolated. Viruses that have a long common evolutionary history with their host typically cause unapparent infections. However, in some Alpha-PV infections, lesions become apparent and may cause benign proliferative disorders or even malignant proliferative lesions of the cervix, vulva, vagina, anus, penis, and oropharynx. The incongruence observed between the topology of the phylogenetic tree of Alpha-PVs and that of their hosts suggests that virus-host codivergence is not the only evolutionary force that has driven the progression of PVs. The integration of the precursors of E5, E6, and E7 on the genome of the ancestral Alpha-PV was important and made the colonization of new niches and the emergence of carcinogenic types possible.

Niche adaptation and viral transmission of human papillomaviruses from archaic hominins to modern humans

Recent discoveries on the origins of modern humans from multiple archaic hominin populations and the diversity of human papillomaviruses (HPVs) suggest a complex scenario of virus-host evolution. To evaluate the origin of HPV pathogenesis, we estimated the phylogeny, timing, and dispersal of HPV16 variants using a Bayesian Markov Chain Monte Carlo framework. To increase precision, we identified and characterized non-human primate papillomaviruses from New and Old World monkeys to set molecular clock models. We demonstrate specific host niche adaptation of primate papillomaviruses with subsequent coevolution with their primate hosts for at least 40 million years. Analyses of 212 HPV16 complete genomes and 3582 partial sequences estimated ancient divergence of HPV16 variants (between A and BCD lineages) from their most recent common ancestors around half a million years ago, roughly coinciding with the timing of the split between archaic Neanderthals and modern Homo sapiens, and nearly three times longer than divergence times of modern Homo sapiens. HPV16 A lineage variants were significantly underrepresented in present African populations, whereas the A sublineages were highly prevalent in European (A1-3) and Asian (A4) populations, indicative of viral sexual transmission from Neanderthals to modern non-African humans through multiple interbreeding events in the past 80 thousand years. Remarkably, the human leukocyte antigen B*07:02 and C*07:02 alleles associated with increased risk in cervix cancer represent introgressed regions from Neanderthals in present-day Eurasians. The archaic hominin-host-switch model was also supported by other HPV variants. Niche adaptation and virus-host codivergence appear to influence the pathogenesis of papillomaviruses.

Genomic characterization of Erethizon dorsatum papillomavirus 2, a new papillomavirus species marked by its exceptional genome size

We report here the complete sequence and genome organization of a new papillomavirus, Erethizon dorsatum papillomavirus 2 (EdPV2), which was isolated from cutaneous lesions observed on the muzzle of a North American porcupine. The complete genome is 8809 nucleotides long and encodes five early (E6-E7-E1-E2-E4) and two late proteins (L2-L1). In addition to the upstream regulatory region, the EdPV2 genome contains an exceptionally large secondary non-coding region with no apparent functional relevance. EdPV2 is strongly divergent from the previously described porcupine papillomavirus EdPV1 and phylogenetic analysis shows EdPV2 clustering near members of the genus Pipapillomavirus, a group of rodent papillomaviruses. Pairwise sequence comparison based on the L1 open reading frame identifies Rattus norvegicus papillomavirus 1 as the closest related virus (59.97 % similarity). Based on its low sequence similarity to other known papillomaviruses, EdPV2 is thought to represent a new genus in the family Papillomaviridae.

Spatio-temporal dynamics and aetiology of proliferative leg skin lesions in wild British finches

Proliferative leg skin lesions have been described in wild finches in Europe although there have been no large-scale studies of their aetiology or epizootiology to date. Firstly, disease surveillance, utilising public reporting of observations of live wild finches was conducted in Great Britain (GB) and showed proliferative leg skin lesions in chaffinches (Fringilla coelebs) to be widespread. Seasonal variation was observed, with a peak during the winter months. Secondly, pathological investigations were performed on a sample of 39 chaffinches, four bullfinches (Pyrrhula pyrrhula), one greenfinch (Chloris chloris) and one goldfinch (Carduelis carduelis) with proliferative leg skin lesions and detected Cnemidocoptes sp. mites in 91% (41/45) of affected finches and from all species examined. Fringilla coelebs papillomavirus (FcPV1) PCR was positive in 74% (23/31) of birds tested: a 394 base pair sequence was derived from 20 of these birds, from all examined species, with 100% identity to reference genomes. Both mites and FcPV1 DNA were detected in 71% (20/28) of birds tested for both pathogens. Histopathological examination of lesions did not discriminate the relative importance of mite or FcPV1 infection as their cause. Development of techniques to localise FcPV1 within lesions is required to elucidate the pathological significance of FcPV1 DNA detection.

Detection and Characterization of Okapi (Okapia johnstoni)-specific Papillomavirus type 1 (OjPV1)

Papillomavirus-specific DNA was detected in skin lesions collected from an okapi (Okapia johnstoni) in the Zoo Basel. According to the nucleotide sequence analysis, the virus belongs to the genus Deltapapillomavirus. Based on bioinformatics analysis, we propose to designate the newly identified virus as Okapia johnstoni Papillomavirus type 1 (OjPV1). OjPV1 is genetically most closely related to a recently described giraffe (Giraffa camelopardalis) -specific papillomavirus (GcPV1). Of note, the putative oncogenic E5 proteins from OjPV1 and GcPV1 are more conserved than the L1 proteins. This indicates, that the selection pressure on E5 may be more pronounced than that on the otherwise most conserved major capsid protein L1.

Molecular detection and immunopathological examination of Deltapapillomavirus 4 in skin and udder of Egyptian cattle

Aim

Bovine papillomaviruses (BPVs) are the main cause of bovine papillomatosis resulting in cutaneous and/or mucosal benign tumors that could be transformed to malignant ones with marked economic importance, especially in the dairy farms. Molecular, pathological, and immunohistochemical (IHC) diagnosis of bovine papillomatosis cases was conducted to identify and characterize the circulating BPV genotype in some Egyptian governorates.

Materials and Methods

Wart-like lesions in skin, udder, and teats were collected from 123 infected cases in Giza, Beni Suef, and El Menoufia Governorates, Egypt, during 2016-2017. Pathological and IHC characterization, molecular identification, genotyping, and phylogenetic analysis based on the conserved late (L1) gene of the all samples were carried out.

Results

89 of the 123 collected samples (72.3%) were positively detected by polymerase chain reaction (PCR). The sequence analysis of the obtained PCR amplicons was identical revealing identification and genotyping of only one type (Deltapapillomavirus 4 isolate EGY 2017) with accession number (MG547343) which found to be closely related to the recently detected Deltapapillomavirus 4 isolate 04_asi_UK (accession no. MF384288.1) and isolate Deltapapillomavirus 4 isolate 25_equ_CH (accession no. MF384286.1) with 99% nucleotide sequence identity. Histopathological examination revealed severe hyperkeratosis in stratum corneum and acanthosis in most of the cases. These tissue changes were confirmed by the presence of golden brown stained proliferating cell nuclear antigen which was localized intranuclear and perinuclear in other cells using IHC Technique.

Conclusion

It is the first time to detect and genotype the BPVs in these areas with no record of previous genotyping in the whole country. The obtained results will highlight the importance of this disease.

Generation of a novel next-generation sequencing-based method for the isolation of new human papillomavirus types

With the advent of new molecular tools, the discovery of new papillomaviruses (PVs) has accelerated during the past decade, enabling the expansion of knowledge about the viral populations that inhabit the human body. Human PVs (HPVs) are etiologically linked to benign or malignant lesions of the skin and mucosa. The detection of HPV types can vary widely, depending mainly on the methodology and the quality of the biological sample. Next-generation sequencing is one of the most powerful tools, enabling the discovery of novel viruses in a wide range of biological material. Here, we report a novel protocol for the detection of known and unknown HPV types in human skin and oral gargle samples using improved PCR protocols combined with next-generation sequencing. We identified 105 putative new PV types in addition to 296 known types, thus providing important information about the viral distribution in the oral cavity and skin.

The Role of aDNA in Understanding the Coevolutionary Patterns of Human Sexually Transmitted Infections

Analysis of pathogen genome data sequenced from clinical and historical samples has made it possible to perform phylogenetic analyses of sexually transmitted infections on a global scale, and to estimate the diversity, distribution, and coevolutionary host relationships of these pathogens, providing insights into pathogen emergence and disease prevention. Deep-sequenced pathogen genomes from clinical studies and ancient samples yield estimates of within-host and between-host evolutionary rates and provide data on changes in pathogen genomic stability and evolutionary responses. Here we examine three groups of pathogens transmitted mainly through sexual contact between modern humans to provide insight into ancient human behavior and history with their pathogens. Exploring ancient pathogen genomic divergence and the ancient viral-host parallel evolutionary histories will help us to reconstruct the origin of present-day geographical distribution and diversity of clinical pathogen infections, and will hopefully allow us to foresee possible environmentally induced pathogen evolutionary responses. Lastly, we emphasize that ancient pathogen DNA research should be combined with modern clinical pathogen data, and be equitable and provide advantages for all researchers worldwide, e.g., through shared data.

Two Novel, Distantly Related Papillomaviruses Isolated from Healthy Skin of the Timor Deer (Rusa timorensis)

We report the complete genome sequences of Rusa timorensis papillomavirus 1 (RtimPV1) and Rusa timorensis papillomavirus 2 (RtimPV2), isolated from hair follicles of asymptomatic skin from the same Timor deer specimen. RtimPV1 and RtimPV2 are evolutionarily only distantly related. RtimPV1 lacks a canonical E2-binding site, and RtimPV2 does not carry an E6 gene.

Molecular identification of papillomavirus in ducks

Papillomaviruses infect many vertebrates, including birds. Persistent infections by some strains can cause malignant proliferation of cells (i.e. cancer), though more typically infections cause benign tumours, or may be completely subclinical. Sometimes extensive, persistent tumours are recorded-notably in chaffinches and humans. In 2016, a novel papillomavirus genotype was characterized from a duck faecal microbiome, in Bhopal, India; the sixth papillomavirus genotype from birds. Prompted by this finding, we screened 160 cloacal swabs and 968 faecal samples collected from 299 ducks sampled at Ottenby Bird Observatory, Sweden in 2015, using a newly designed real-time PCR. Twenty one samples (1.9%) from six individuals (2%) were positive. Eighteen sequences were identical to the published genotype, duck papillomavirus 1. One additional novel genotype was recovered from three samples. Both genotypes were recovered from a wild strain domestic mallard that was infected for more than 60 days with each genotype. All positive individuals were adult (P = 0.004). Significantly more positive samples were detected from swabs than faecal samples (P < 0.0001). Sample type data suggests transmission may be via direct contact, and only infrequently, via the oral-faecal route. Infection in only adult birds supports the hypothesis that this virus is sexually transmitted, though more work is required to verify this.

Genome Sequences of Two Novel Papillomaviruses Isolated from Healthy Skin of Pudu puda and Cervus elaphus Deer

We report the complete genome sequences of Pudu puda papillomavirus1 (PpudPV1) and Cervus elaphus papillomavirus2 (CelaPV2), isolated from healthy skin hair follicles of a Southern pudu and a red deer, respectively. PpudPV1 is basal to the DyokappaPVs, whereas CelaPV2 is basal to the XiPVs (Beta-XiPV crown group).

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.

Rodent Papillomaviruses

Preclinical infection model systems are extremely valuable tools to aid in our understanding of Human Papillomavirus (HPV) biology, disease progression, prevention, and treatments. In this context, rodent papillomaviruses and their respective infection models are useful tools but remain underutilized resources in the field of papillomavirus biology. Two rodent papillomaviruses, MnPV1, which infects the Mastomys species of multimammate rats, and MmuPV1, which infects laboratory mice, are currently the most studied rodent PVs. Both of these viruses cause malignancy in the skin and can provide attractive infection models to study the lesser understood cutaneous papillomaviruses that have been frequently associated with HPV-related skin cancers. Of these, MmuPV1 is the first reported rodent papillomavirus that can naturally infect the laboratory strain of mice. MmuPV1 is an attractive model virus to study papillomavirus pathogenesis because of the ubiquitous availability of lab mice and the fact that this mouse species is genetically modifiable. In this review, we have summarized the knowledge we have gained about PV biology from the study of rodent papillomaviruses and point out the remaining gaps that can provide new research opportunities.

Genomic comparison of bovine papillomavirus 1 isolates from bovine, equine and asinine lesional tissue samples

Several attempts have been made to categorize equid- and bovid-specific bovine papillomavirus 1 (BPV1) isolates based on sequence tags. This study includes newly determined sequence information from 33 BPV1 isolates of equine, asinine and bovine origin and investigates sequence bias due to host species. Twenty of the viral genomes were sequenced over their entire length and a further thirteen were sequenced, including flanking sequences, at two specific sites, the LCR and the E5 ORF. Alignment and analyses of the sequences did not reveal statistically significant site differences between the sequences of bovine and equid origin. None of the proposed sites of divergence noted by other authors demonstrated significant species-specific characteristics. Our results suggest that BPV1 is shared between equine, asinine and bovine host species, and that viral transfer between bovines and equids is a repeated and ongoing phenomenon.

Ancient Evolution and Dispersion of Human Papillomavirus 58 Variants

Human papillomavirus 58 (HPV58) is found in 10 to 18% of cervical cancers in East Asia but is rather uncommon elsewhere. The distribution and oncogenic potential of HPV58 variants appear to be heterogeneous, since the E7 T20I/G63S variant is more prevalent in East Asia and confers a 7- to 9-fold-higher risk of cervical precancer and cancer. However, the underlying genomic mechanisms that explain the geographic and carcinogenic diversity of HPV58 variants are still poorly understood. In this study, we used a combination of phylogenetic analyses and bioinformatics to investigate the deep evolutionary history of HPV58 complete genome variants. The initial splitting of HPV58 variants was estimated to occur 478,600 years ago (95% highest posterior density [HPD], 391,000 to 569,600 years ago). This divergence time is well within the era of speciation between Homo sapiens and Neanderthals/Denisovans and around three times longer than the modern Homo sapiens divergence times. The expansion of present-day variants in Eurasia could be the consequence of viral transmission from Neanderthals/Denisovans to non-African modern human populations through gene flow. A whole-genome sequence signature analysis identified 3 amino acid changes, 16 synonymous nucleotide changes, and a 12-bp insertion strongly associated with the E7 T20I/G63S variant that represents the A3 sublineage and carries higher carcinogenetic potential. Compared with the capsid proteins, the oncogenes E7 and E6 had increased substitution rates indicative of higher selection pressure. These data provide a comprehensive evolutionary history and genomic basis of HPV58 variants to assist further investigation of carcinogenic association and the development of diagnostic and therapeutic strategies.IMPORTANCE Papillomaviruses (PVs) are an ancient and heterogeneous group of double-stranded DNA viruses that preferentially infect the cutaneous and mucocutaneous epithelia of vertebrates. Persistent infection by specific oncogenic human papillomaviruses (HPVs), including HPV58, has been established as the primary cause of cervical cancer. In this work, we reveal the complex evolutionary history of HPV58 variants that explains the heterogeneity of oncogenic potential and geographic distribution. Our data suggest that HPV58 variants may have coevolved with archaic hominins and dispersed across the planet through host interbreeding and gene flow. Certain genes and codons of HPV58 variants representing higher carcinogenic potential and/or that are under positive selection may have important implications for viral host specificity, pathogenesis, and disease prevention.

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.

HPV8-E6 Interferes with Syntenin-2 Expression through Deregulation of Differentiation, Methylation and Phosphatidylinositide-Kinase Dependent Mechanisms

The E6 oncoproteins of high-risk human papillomaviruses (HPV) of genus alpha contain a short peptide sequence at the carboxy-terminus, the PDZ binding domain, with which they interact with the corresponding PDZ domain of cellular proteins. Interestingly, E6 proteins from papillomaviruses of genus beta (betaPV) do not encode a comparable PDZ binding domain. Irrespective of this fact, we previously showed that the E6 protein of HPV8 (betaPV type) could circumvent this deficit by targeting the PDZ protein Syntenin-2 through transcriptional repression (Lazic et al., 2012). Despite its high binding affinity to phosphatidylinositol-4,5-bisphosphate (PI(4,5)P₂), very little is known about Syntenin-2. This study aimed to extend the knowledge on Syntenin-2 and how its expression is controlled. We now identified that Syntenin-2 is expressed at high levels in differentiating and in lower amounts in keratinocytes cultured in serum-free media containing low calcium concentration. HPV8-E6 led to a further reduction of Syntenin-2 expression only in cells cultured in low calcium. In the skin of patients suffering from Epidermodysplasia verruciformis, who are predisposed to betaPV infection, Syntenin-2 was expressed in differentiating keratinocytes of non-lesional skin, but was absent in virus positive squamous tumors. Using 5-Aza-2′-deoxycytidine, which causes DNA demethylation, Syntenin-2 transcription was profoundly activated and fully restored in the absence and presence of HPV8-E6, implicating that E6 mediated repression of Syntenin-2 transcription is due to promoter hypermethylation. Since Syntenin-2 binds to PI(4,5)P₂, we further tested whether the PI(4,5)P₂ metabolic pathway might govern Syntenin-2 expression. PI(4,5)P₂ is generated by the activity of phosphatidylinositol-4-phosphate-5-kinase type I (PIP5KI) or phosphatidylinositol-5-phosphate-4-kinase type II (PIP4KII) isoforms α, β and γ. Phosphatidylinositide kinases have recently been identified as regulators of gene transcription. Surprisingly, transfection of siRNAs directed against PIP5KI and PIP4KII resulted in higher Syntenin-2 expression with the highest effect mediated by siPIP5KIα. HPV8-E6 was able to counteract siPIP4KIIα, siPIP4KIIβ and siPIP5KIγ mediated Syntenin-2 re-expression but not siPIP5KIα. Finally, we identified Syntenin-2 as a key factor regulating PIP5KIα expression. Collectively, our data demonstrates that Syntenin-2 is regulated through multiple mechanisms and that downregulation of Syntenin-2 expression may contribute to E6 mediated dedifferentiation of infected skin cells.