Origin and evolution of papillomavirus (onco)genes and genomes

Novel Gammapapillomavirus type in the nasal cavity of a wild red colobus (Piliocolobus tephrosceles)

Papillomaviruses (PVs) are double-stranded, circular, epitheliotropic DNA viruses causing benign warts (papillomas) or inducing dysplasia that can progress to cancer. Although they have been identified in all vertebrate taxa, most classified types are human PVs (HPVs); relatively little is known about PVs in other species. Here we characterize a novel Gammapapillomavirus type, PtepPV1, from a nasal swab of a wild red colobus (Piliocolobus tephrosceles) in Kibale National Park, Uganda. The virus has a genome of 6576 bases, encoding the seven canonical early (E) ORFs (E6, E7, E1, E2, E4, E1^E4 and E8^E2) and two late (L) ORFs (L1 and L2) of the gammapapillomaviruses, and is 81.0% similar to HPV-mSK_118, detected in a cutaneous wart from an immunocompromised human patient, in the L1 gene at the amino acid level. Alphapapillomaviruses (genus Alphapapillomavirus) cause anogenital carcinomas such as cervical cancer and have been described previously in several nonhuman primates. However, the first gammapapillomavirus (genus Gammapapillomavirus), which cause transient cutaneous infections, was not described until 2019 in a healthy rhesus macaque (Macaca mulatta) genital swab. The new virus from red colobus, PtepPV1, has many genomic features encoded by high-risk oncogenic PVs, such as the E7 gene LXSXE and CXXC motifs, suggesting potential for pRb and zinc-finger binding, respectively. To our knowledge, PtepPV1 is also the first reported nonhuman primate PV found in the nasal cavity. PtepPV1 expands the known host range, geographical distribution, tissue tropism and biological characteristics of nonhuman primate PVs.

Determination and Characterization of Novel Papillomavirus and Parvovirus Associated with Mass Mortality of Chinese Tongue Sole (Cynoglossus semilaevis) in China

A massive mortality event concerning farmed Chinese tongue soles occurred in Tianjin, China, and the causative agent remains unknown. Here, a novel Cynoglossus semilaevis papillomavirus (CsPaV) and parvovirus (CsPV) were simultaneously isolated and identified from diseased fish via electron microscopy, virus isolation, genome sequencing, experimental challenges, and fluorescence in situ hybridization (FISH). Electron microscopy showed large numbers of virus particles present in the tissues of diseased fish. Viruses that were isolated and propagated in flounder gill cells (FG) induced typical cytopathic effects (CPE). The cumulative mortality of fish given intraperitoneal injections reached 100% at 7 dpi. The complete genomes of CsPaV and CsPV comprised 5939 bp and 3663 bp, respectively, and the genomes shared no nucleotide sequence similarities with other viruses. Phylogenetic analysis based on the L1 and NS1 protein sequences revealed that CsPaV and CsPV were novel members of the Papillomaviridae and Parvoviridae families. The FISH results showed positive signals in the spleen tissues of infected fish, and both viruses could co-infect single cells. This study represents the first report where novel papillomavirus and parvovirus are identified in farmed marine cultured fish, and it provides a basis for further studies on the prevention and treatment of emerging viral diseases.

Potential of Viruses as Environmental Etiological Factors for Non-Syndromic Orofacial Clefts

In this study, we analyzed the potential of viral infections in the species Homo sapiens as environmental causes of orofacial clefts (OFCs). A scoring system was adapted for qualitatively assessing the potential of viruses to cause cleft lip and/or palate (CL/P). This assessment considered factors such as information from the literature, nucleotide and amino acid similarities, and the presence of Endogenous Viral Elements (EVEs). The analysis involved various algorithm packages within Basic Local Alignment Search Tool 2.13.0 software and databases from the National Center for Biotechnology Information and the International Committee on Taxonomy of Viruses. Twenty significant viral species using different biosynthesis strategies were identified: Human coronavirus NL63, Rio Negro virus, Alphatorquevirus homin9, Brisavirus, Cosavirus B, Torque teno mini virus 4, Bocaparvovirus primate2, Human coronavirus HKU1, Monkeypox virus, Mammarenavirus machupoense, Volepox virus, Souris mammarenavirus, Gammapapillomavirus 7, Betainfluenzavirus influenzae, Lymphocytic choriomeningitis mammarenavirus, Ledantevirus kern, Gammainfluenzavirus influenzae, Betapolyomavirus hominis, Vesiculovirus perinet, and Cytomegalovirus humanbeta5. The evident viral etiological potential in relation to CL/P varies depending on the Baltimore class to which the viral species belongs. Given the multifactorial nature of CL/P, this relationship appears to be dynamic.

Novel equid papillomavirus from domestic donkey

BACKGROUND

Papillomaviruses can be of great medical importance as they infect humans and animals such as Equus species, other livestock and pets. They are responsible for several papillomas and benign tumours in their host.

OBJECTIVES

To describe a novel equid papillomavirus detected in oral swab samples collected from donkeys (Equus asinus) found on the Northwest plateau of China.

STUDY DESIGN

Cross-sectional.

METHODS

Swab samples collected from the oral mucosa of 32 donkeys in the Gansu Province of China, were subjected to viral metagenomic analysis to detect the presence of Papillomavirus. After de novo assembly, a novel papillomavirus genome designated as Equus asinus papillomavirus 3 (EaPV3) was identified in the studied samples. Additional bioinformatic analysis of the assembled genome was done using the Geneious prime software (version 2022.0.2).

RESULTS

The complete circular genome of EaPV3 is 7430 bp in length with a GC content of 50.8%. The genome was predicted to contain five ORFs coding for three early proteins (E7, E1, and E2) and two late proteins (L1 and L2). Phylogenic analysis of the nucleotide sequences of the concatenated amino acid sequences of the E1E2L1L2 genes revealed that EaPV3 is most closely related to Equus asinus papillomavirus 1 (EaPV1). The genome analysis of EaPV3 revealed similar genome organisation with other equine papillomavirus and the presence of E7 papillomavirus oncoprotein.

MAIN LIMITATIONS

Since there were no warts in the oral cavity of the donkeys in this study, and no biopsy samples taken, we are unable to conclusively link the novel virus to any clinical condition in the donkeys.

CONCLUSIONS

The Comparative characterisation of EaPV3 and its closest relatives, as well as phylogenetic analysis demonstrated that it is a novel virus specie that clusters within the Dyochipapilloma PV genus.

The importance of combined HPV and CINtec® PLUS genotyping testing for p16 in women with cervical squamous cell carcinoma

INTRODUCTION

Immunohistochemistry (IHC) for p16INK4A (p16) is a reliable surrogate test for the presence of a high-risk, potentially transformative human papillomavirus (HPV) infection in precursor and malignant lesions of the cervix. The purpose of this study was to evaluate changes in cervical cells caused by persistent HPV infection, by IHC (p16 protein) by comparison with HPV genotyping.

PATIENTS, MATERIALS AND METHODS

The study included female patients aged between 26 and 57 years who presented to a public hospital, with complaints related to the genital area, namely vaginal bleeding and dyspareunia. After selecting the patients, samples were subjected to cytological testing and IHC for p16 and for the determination of HPV types.

RESULTS

The relationship between HPV status and p16 status was statistically significant (p=0.0001), of the 41 patients, 53.7% were HPV positive, respectively 56.1% were p16 positive, the agreement relationship between the two indicators was very high (weighted kappa: 0.951). The clinical performance of CINtec® PLUS triage for p16 shows a high positive predictive value (PPV) and a high negative predictive value (NPV) of 95.7% and 100%, respectively, as regards HPV.

CONCLUSIONS

The p16 marker (CINtec® PLUS) can be used as a prognostic biomarker and provides clinical usefulness through increased sensitivity (Se) and specificity (Sp) in the triage of women at risk of developing precancerous lesions, compared to cytology that is based on morphology, but has a rather low Se and high Sp, while HPV testing is very sensitive but slightly more specific.

Diagnosis and phylogenetic analysis of bovine papillomaviruses in cattle papillomatosis cases by different methods

Papillomaviruses, known as epitheliotropic, cause proliferation in the skin, mucosa, and different visceral organs. In this study, it was aimed to diagnose bovine papillomavirus (BPV) by using different methods in the lesion taken from twenty cattle with papillomas in different areas of the body and to reveal its molecular characterization. In our study, molecular, immunohistochemistry, and transmission electron microscopy (TEM) methods were used for virus identification. Additionally, sequencing analysis was used to ascertain the phylogenetic relationship between the obtained field strains and other isolates submitted to GenBank. Histopathological analyses of the collected samples were done in addition to diagnostic procedures. Intranuclear virus particles were detected when the papillomas were investigated with TEM. In PCR analyses using degenerate and type-specific primer sets, the presence of BPV nucleic acid was determined in 70% (14/20) and 90% (18/20) of the samples, respectively. No virus could be detected in PCR applications using MY 09/11 degenerate primer sets. Twenty animals of different ages, races, and genders included in the study by random sampling method from different herds were divided into 4 groups according to the body regions where the lesions were located. Sequence analysis was performed on a sample from each group that showed strong positivity in the PCR technique using FAP 59/64 degenerate primer set and type-specific primer set. Sequence analyses were performed using FAP 59/64 degenerate primers of amplicons for phylogenetic research. In these analyses, three of the isolated strains were identified as BPV-1, which is in the Deltapapillomavirus 4 genus, and one as BPV-2. As a result of the study, it was concluded that molecular and phylogenetic studies using type-specific primers are more beneficial in order to fully reveal the etiology of papillomatosis in cattle and it would be correct to determine BPV types before prophylactic (vaccine, etc.) applications.

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

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

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.

Human Papillomavirus: Challenges and Opportunities for the Control of Cervical Cancer

Viruses are the most abundant and genetically diverse entities on the planet, infect all life forms and have evolved with their hosts. To date, 263 viral species have been identified that infect humans, of which only seven are considered type I oncogenic. Human papillomavirus (HPV) is the main virus associated with cancer and is responsible for practically all cases of cervical carcinoma. Screening tests for early detection have been available since the 1960s. Undoubtedly, the entailment between knowledge of HPV biology and the natural history of cervical cancer has contributed to the significant advances that have been made for its prevention since the 21^st century, with the development of prophylactic vaccines and improved screening strategies. Therefore, it is possible to eradicate invasive cervical cancer as a worldwide public health problem, as proposed by the WHO with the 90-70-90 initiative based on vaccination coverage, screening, and treatment, respectively. In addition, the emerging knowledge of viral biology generates opportunities that will contribute to strengthening prevention and treatment strategies in HPV-associated neoplasms.

Revisiting Papillomavirus Taxonomy: A Proposal for Updating the Current Classification in Line with Evolutionary Evidence

Papillomaviruses infect a wide array of animal hosts and are responsible for roughly 5% of all human cancers. Comparative genomics between different virus types belonging to specific taxonomic groupings (e.g., species, and genera) has the potential to illuminate physiological differences between viruses with different biological outcomes. Likewise, extrapolation of features between related viruses can be very powerful but requires a solid foundation supporting the evolutionary relationships between viruses. The current papillomavirus classification system is based on pairwise sequence identity. However, with the advent of metagenomics as facilitated by high-throughput sequencing and molecular tools of enriching circular DNA molecules using rolling circle amplification, there has been a dramatic increase in the described diversity of this viral family. Not surprisingly, this resulted in a dramatic increase in absolute number of viral types (i.e., sequences sharing <90% L1 gene pairwise identity). Many of these novel viruses are the sole member of a novel species within a novel genus (i.e., singletons), highlighting that we have only scratched the surface of papillomavirus diversity. I will discuss how this increase in observed sequence diversity complicates papillomavirus classification. I will propose a potential solution to these issues by explicitly basing the species and genera classification on the evolutionary history of these viruses based on the core viral proteins (E1, E2, and L1) of papillomaviruses. This strategy means that it is possible that a virus identified as the closest neighbor based on the E1, E2, L1 phylogenetic tree, is not the closest neighbor based on L1 nucleotide identity. In this case, I propose that a virus would be considered a novel type if it shares less than 90% identity with its closest neighbors in the E1, E2, L1 phylogenetic tree.

Revealing the uncharacterised diversity of amphibian and reptile viruses

Amphibians and non-avian reptiles represent a significant proportion of terrestrial vertebrates, however knowledge of their viruses is not proportional to their abundance. Many amphibians and reptiles have strict habitual environments and localised populations and are vulnerable to viral outbreaks and potential elimination as a result. We sought to identify viruses that were hidden in amphibian and reptile metatranscriptomic data by screening 235 RNA-sequencing datasets from a 122 species covering 25 countries. We identified 26 novel viruses and eight previously characterised viruses from fifteen different viral families. Twenty-five viruses had RNA genomes with identity to Arteriviridae, Tobaniviridae, Hantaviridae, Rhabdoviridae, Astroviridae, Arenaviridae, Hepeviridae, Picornaviridae, Orthomyxoviridae, Reoviridae, Flaviviridae and Caliciviridae. In addition to RNA viruses, we also screened datasets for DNA viral transcripts, which are commonly excluded from transcriptomic analysis. We identified ten DNA viruses with identity to Papillomaviridae, Parvoviridae, Circoviridae and Adomaviridae. With the addition of these viruses, we expand the global amphibian and reptile virome and identify new potentially pathogenic viruses that could challenge populations. We speculate that amphibian viruses often have simpler genomes than those in amniotes, as in the case of the Secondpapillomavirinae and Orthomyxoviridae viruses identified in this study. In addition, we find evidence of inter-family recombination in RNA viruses, and we also identify new members of the recombinant Adomaviridae family. Overall, we provide insights into the uncharacterised diversity of amphibian and reptile viruses with the aim of improving population management, treatment and conservation into the future.

The Drivers, Mechanisms, and Consequences of Genome Instability in HPV-Driven Cancers

Simple Summary

Cells infected with high-risk human papillomaviruses (HPV) can accumulate DNA damage and eventually transform into HPV-driven cancers. Genome instability, or the progressive accumulation of DNA alterations (e.g., mutations), in HPV-infected cells is directly induced by the HPV genes and indirectly promoted by HPV infection through the consequences of chronic infection maintenance, increased cell growth, and accumulation of damaging mutations in genes that themselves affect genome instability. While the HPV genome typically exists as a separate entity within cells, genome instability increases the chances of HPV integrating within the host (human) genome, which is common in HPV-induced cancers. The DNA regions surrounding HPV integrations are unstable and can undergo complex alterations that affect both human and HPV genes. This review discusses HPV-dependent and -independent drivers and mechanisms of genome instability in HPV-driven cancers, both globally and around sites of HPV integration, and describes the changes induced in the tumour genome.

Abstract

Human papillomavirus (HPV) is the causative driver of cervical cancer and a contributing risk factor of head and neck cancer and several anogenital cancers. HPV’s ability to induce genome instability contributes to its oncogenicity. HPV genes can induce genome instability in several ways, including modulating the cell cycle to favour proliferation, interacting with DNA damage repair pathways to bring high-fidelity repair pathways to viral episomes and away from the host genome, inducing DNA-damaging oxidative stress, and altering the length of telomeres. In addition, the presence of a chronic viral infection can lead to immune responses that also cause genome instability of the infected tissue. The HPV genome can become integrated into the host genome during HPV-induced tumorigenesis. Viral integration requires double-stranded breaks on the DNA; therefore, regions around the integration event are prone to structural alterations and themselves are targets of genome instability. In this review, we present the mechanisms by which HPV-dependent and -independent genome instability is initiated and maintained in HPV-driven cancers, both across the genome and at regions of HPV integration.

Differences in integration frequencies and APOBEC3 profiles of five high-risk HPV types adheres to phylogeny

Persistent infection with Human Papillomavirus (HPV) is responsible for almost all cases of cervical cancers, and HPV16 and HPV18 associated with the majority of these. These types differ in the proportion of viral minor nucleotide variants (MNVs) caused by APOBEC3 mutagenesis as well as integration frequencies. Whether these traits extend to other types remains uncertain. This study aimed to investigate and compare genomic variability and chromosomal integration in the two phylogenetically distinct Alpha-7 and Alpha-9 clades of carcinogenic HPV types. The TaME-seq protocol was employed to sequence cervical cell samples positive for HPV31, HPV33 or HPV45 and combine these with data from a previous study on HPV16 and HPV18. APOBEC3 mutation signatures were found in Alpha-9 (HPV16/31/33) but not in Alpha-7 (HPV18/45). HPV45 had significantly more MNVs compared to the other types. Alpha-7 had higher integration frequency compared to Alpha-9. An increase in integration frequency with increased diagnostic severity was found for Alpha-7. The results highlight important differences and broaden our understanding of the molecular mechanisms behind cervical cancer induced by high-risk HPV types from the Alpha-7 and Alpha-9 clades.

HPV Integration Site Mapping: A Rapid Method of Viral Integration Site (VIS) Analysis and Visualization Using Automated Workflows in CLC Microbial Genomics

Human papillomavirus (HPV) integration within the host genome may contribute to carcinogenesis through various disruptive mechanisms. With next-generation sequencing (NGS), identification of viral and host genomic breakpoints and chimeric sequences are now possible. However, a simple, streamlined bioinformatics workflow has been non-existent until recently. Here, we tested two new, automated workflows in CLC Microbial Genomics, i.e., Viral Hybrid Capture (VHC) Data Analysis and Viral Integration Site (VIS) Identification for software performance and efficiency. The workflows embedded with HPV and human reference genomes were used to analyze a publicly available NGS dataset derived from pre- and cancerous HPV+ cervical cytology of 21 Gabonese women. The VHC and VIS workflow median runtimes were 19 and 7 min per sample, respectively. The VIS dynamic graphical outputs included read mappings, virus-host genomic breakpoints, and virus-host integration circular plots. Key findings, including disrupted and nearby genes, were summarized in an auto-generated report. Overall, the VHC and VIS workflows proved to be a rapid and accurate means of localizing viral-host integration site(s) and identifying disrupted and neighboring human genes. Applying HPV VIS-mapping to pre- or invasive tumors will advance our understanding of viral oncogenesis and facilitate the discovery of prognostic biomarkers and therapeutic targets.

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.

Viruses Infecting the European Catfish (Silurus glanis)

In aquaculture, disease management and pathogen control are key for a successful fish farming industry. In past years, European catfish farming has been flourishing. However, devastating fish pathogens including limiting fish viruses are considered a big threat to further expanding of the industry. Even though mainly the ranavirus (Iridoviridea) and circovirus (Circoviridea) infections are considered well- described in European catfish, more other agents including herpes-, rhabdo or papillomaviruses are also observed in the tissues of catfish with or without any symptoms. The etiological role of these viruses has been unclear until now. Hence, there is a requisite for more detailed information about the latter and the development of preventive and therapeutic approaches to complete them. In this review, we summarize recent knowledge about viruses that affect the European catfish and describe their origin, distribution, molecular characterisation, and phylogenetic classification. We also highlight the knowledge gaps, which need more in-depth investigations in the future.

HPV DeepSeq: An Ultra-Fast Method of NGS Data Analysis and Visualization Using Automated Workflows and a Customized Papillomavirus Database in CLC Genomics Workbench

Next-generation sequencing (NGS) has actualized the human papillomavirus (HPV) virome profiling for in-depth investigation of viral evolution and pathogenesis. However, viral computational analysis remains a bottleneck due to semantic discrepancies between computational tools and curated reference genomes. To address this, we developed and tested automated workflows for HPV taxonomic profiling and visualization using a customized papillomavirus database in the CLC Microbial Genomics Module. HPV genomes from Papilloma Virus Episteme were customized and incorporated into CLC “ready-to-use” workflows for stepwise data processing to include: (1) Taxonomic Analysis, (2) Estimate Alpha/Beta Diversities, and (3) Map Reads to Reference. Low-grade (n = 95) and high-grade (n = 60) Pap smears were tested with ensuing collective runtimes: Taxonomic Analysis (36 min); Alpha/Beta Diversities (5 s); Map Reads (45 min). Tabular output conversion to visualizations entailed 1–2 keystrokes. Biodiversity analysis between low- (LSIL) and high-grade squamous intraepithelial lesions (HSIL) revealed loss of species richness and gain of dominance by HPV-16 in HSIL. Integrating clinically relevant, taxonomized HPV reference genomes within automated workflows proved to be an ultra-fast method of virome profiling. The entire process named “HPV DeepSeq” provides a simple, accurate and practical means of NGS data analysis for a broad range of applications in viral research.

Complete genome sequence of a novel fish papillomavirus detected in farmed wels catfish (Silurus glanis)

A novel papillomavirus (PV) was detected in farmed wels catfish (Silurus glanis) in Hungary showing clinical signs resembling those of wels catfish herpesvirus disease. The whole genome of Silurus glanis papillomavirus 1 (SgPV1) was identified using next-generation sequencing. The 5,612-bp complete genome contains four predicted protein coding regions (E1, E2, L1, and L2), which seem to have homologues in every PV genome sequenced to date. Five complete fish PV genome sequences are available in the GenBank database. Their genomes range between 5,748 and 6,086 bp and contain the minimal PV backbone genes E1, E2, L2, and L1, unlike PVs of higher vertebrates, which have larger genomes (6.8-8.6 kbp) and additional (onco)genes. Considering the current species demarcation criteria for the family Papillomaviridae, the establishment of a novel species named "Nunpapillomavirus siluri" is proposed for the SgPV1 in a novel genus, "Nunpapillomavirus", in the subfamily Secondpapillomavirinae.

Recurrent infantile digital fibromatosis with HPV infection: a case report

Infantile digital fibromatosis is a rare, benign fibrous tumor among infants usually limited to fingers and toes. Most cases present themselves with round eosinophilic inclusion bodies of different sizes, similar to erythrocyte in the parakarytoplasm. Although infantile digital fibromatosis had been observed with a tendency of spontaneous regression after a rapid phase of growth in some reports, the recurrence rate of early surgical intervention remains high. And the mechanism of recurrence is still unknown. Human papillomavirus (HPV), as a circular icosahedral double-stranded DNA virus, is famous for its remarkable significant correlation with cervical cancer. However, the reports about the possible relationship of recurrent infantile digital fibromatosis and HPV infection are rare and inconsistent. Here, we report a recurrent case of infantile digital fibromatosis after surgical resection. Pathological biopsy of the recurrent site not only identified the diagnosis of infantile digital fibromatosis again, but found the sign of HPV infection. Family history indicated that the patient's grandmother had a history of verruca plana. After complete resection of recurrence, the tumor recurred again and the case is still being followed up. The unique case may serve as a clue to the pathogenesis of the relationship between recurrent infantile digital fibromatosis and HPV infection.

The oncogenic pathways of papillomaviruses

Papillomaviruses are oncogenic DNA viruses and induce hyperplastic benign lesions of both cutaneous and mucosal tissues in their various hosts, including many domestic and wild animals as well as humans. There are some Papillomavirus genotypes that can infect hosts different from their own, such as BPV 1 and BPV 2 originated from cattle, which can also infect horses and are responsible for fibroblastic tumours in horses. This review article summarizes the origin and evolution of papillomaviruses as an etiological agent in the historical process. The main focus in this review is the evaluation of the interactions between high-risk papillomavirus oncoproteins and programmed cell-death pathways. It further exemplifies the role of these interactions in the malignant cell transformation process. In parallel with this, the use and importance of the bovine model system to enlighten the papillomavirus-associated cancers is discussed with an in-depth examination. Furthermore, it focuses on the epidemiological situation of BPV infections in Turkey in the cattle herds.

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.

APOBEC: A molecular driver in cervical cancer pathogenesis

Cervical cancer is one of the foremost common cancers in women. Human papillomavirus (HPV) infection remains a major risk factor of cervical cancer. In addition, numerous other genetic and epigenetic factors also are involved in the underlying pathogenesis of cervical cancer. Recently, it has been reported that apolipoprotein B mRNA editing enzyme catalytic polypeptide like (APOBEC), DNA-editing protein plays an important role in the molecular pathogenesis of cancer. Particularly, the APOBEC3 family was shown to induce tumor mutations by aberrant DNA editing mechanism. In general, APOBEC3 enzymes play a pivotal role in the deamination of cytidine to uridine in DNA and RNA to control diverse biological processes such as regulation of protein expression, innate immunity, and embryonic development. Innate antiviral activity of the APOBEC3 family members restrict retroviruses, endogenous retro-element, and DNA viruses including the HPV that is the leading risk factor for cervical cancer. This review briefly describes the pathogenesis of cervical cancer and discusses in detail the recent findings on the role of APOBEC in the molecular pathogenesis of cervical cancer.

•

APOBEC enzymes deaminate cytidine to uridine and control diverse biological processes including viral restriction.

•

APOBEC3, DNA/RNA-editing enzyme plays an important role in the molecular pathogenesis of cervical cancer.

•

APOBEC3-mediated DNA editing leads to the accumulation of somatic mutations in tumors and HPV genome.

•

Deregulation of APOBEC3 family genes cause genomic instability and result in drug resistance, and immune-evasion in tumors.

Genomic characterization of a novel bovine papillomavirus type 28

Infection of bovine papillomavirus (BPV) has been associated with mucosal and/or cutaneous tumor development in bovids. To date, up to 27 genotypes of BPVs have been identified and classified based on the nucleotide sequence identity of L1 open reading frame. In the present study, the complete sequence of a novel BPV concurrently identified with BPV1 and BPV2 in the facial cutaneous papilloma lesion of a domestic cattle was characterized. The whole genome of the unclassified BPV was 7263 base pairs in full length with GC ratio of 42.9%. In comparison with published BPV sequences, L1 nucleotide sequence of the novel BPV shared 75% identity with BPV15, and was suggested to be classified in the genus, Xipapillomavirus. According to the criteria established by the International Committee on the Taxonomy of Viruses, the novel BPV was designated as BPV type 28.

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.

Towards a multi-level and a multi-disciplinary approach to DNA oncovirus virulence

One out of 10 cancers is estimated to arise from infections by a handful of oncogenic viruses. These infectious cancers constitute an opportunity for primary prevention through immunization against the viral infection, for early screening through molecular detection of the infectious agent, and potentially for specific treatments, by targeting the virus as a marker of cancer cells. Accomplishing these objectives will require a detailed understanding of the natural history of infections, the mechanisms by which the viruses contribute to disease, the mutual adaptation of viruses and hosts, and the possible viral evolution in the absence and in the presence of the public health interventions conceived to target them. This issue showcases the current developments in experimental tissue-like and animal systems, mathematical models and evolutionary approaches to understand DNA oncoviruses. Our global aim is to provide proximate explanations to the present-day interface and interactions between virus and host, as well as ultimate explanations about the adaptive value of these interactions and about the evolutionary pathways that have led to the current malignant phenotype of oncoviral infections. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

Modelling the evolution of viral oncogenesis

Most human oncogenic viruses share several characteristics, such as being DNA viruses, having long (co)evolutionary histories with their hosts and causing either latent or chronic infections. They can reach high prevalences while causing relatively low case mortality, which makes them quite fit according to virulence evolution theory. After analysing the life histories of DNA oncoviruses, we use a mathematical modelling approach to investigate how the virus life cycle may generate selective pressures favouring or acting against oncogenesis at the within-host or at the between-host level. In particular, we focus on two oncoprotein activities, namely extending cell life expectancy and increasing cell proliferation rate. These have immediate benefits (increasing viral population size) but can be associated with fitness costs at the epidemiological level (increasing recovery rate or risk of cancer) thus creating evolutionary trade-offs. We interpret the results of our nested model in light of the biological features and identify future perspectives for modelling oncovirus dynamics and evolution. This article is part of the theme issue 'Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses'.

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.

Evolutionary dynamics of ten novel Gamma-PVs: insights from phylogenetic incongruence, recombination and phylodynamic analyses

Background

Human papillomaviruses (HPVs) are genetically diverse, belonging to five distinct genera: Alpha, Beta, Gamma, Mu and Nu. All papillomaviruses have double stranded DNA genomes that are thought to evolve slowly because they co-opt high-fidelity host cellular DNA polymerases for their replication. Despite extensive efforts to catalogue all the HPV species that infect humans, it is likely that many still remain undiscovered. Here we use the sequences of ten novel Gammapapillomaviruses (Gamma-PVs) characterized in previous studies and related HPVs to analyse the evolutionary dynamics of these viruses at the whole genome and individual gene scales.

Results

We found statistically significant incongruences between the phylogenetic trees of different genes which imply gene-to-gene variation in the evolutionary processes underlying the diversification of Gamma-PVs. We were, however, only able to detect convincing evidence of a single recombination event which, on its own, cannot explain the observed incongruences between gene phylogenies. The divergence times of the last common ancestor (LCA) of the Alpha, Beta, Mu, Nu and Gamma genera was predicted to have existed between 49.7–58.5 million years ago, before splitting into the five main lineages. The LCA of the Gamma-PVs at this time was predicted to have existed between 45.3 and 67.5 million years ago: approximately at the time when the simian and tarsier lineages of the primates diverged.

Conclusion

Consequently, we report here phylogenetic tree incongruence without strong evidence of recombination.