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

HPV Infections-Classification, Pathogenesis, and Potential New Therapies

To date, more than 400 types of human papillomavirus (HPV) have been identified. Despite the creation of effective prophylactic vaccines against the most common genital HPVs, the viruses remain among the most prevalent pathogens found in humans. According to WHO data, they are the cause of 5% of all cancers. Even more frequent are persistent and recurrent benign lesions such as genital and common warts. HPVs are resistant to many disinfectants and relatively unsusceptible to external conditions. There is still no drug available to inhibit viral replication, and treatment is based on removing lesions or stimulating the host immune system. This paper presents the systematics of HPV and the differences in HPV structure between different genetic types, lineages, and sublineages, based on the literature and GenBank data. We also present the pathogenesis of diseases caused by HPV, with a special focus on the role played by E6, E7, and other viral proteins in the development of benign and cancerous lesions. We discuss further prospects for the treatment of HPV infections, including, among others, substances that block the entry of HPV into cells, inhibitors of viral early proteins, and some substances of plant origin that inhibit viral replication, as well as new possibilities for therapeutic vaccines.

Identification and characterization of novel alphacoronaviruses in Tadarida brasiliensis (Chiroptera, Molossidae) from Argentina: insights into recombination as a mechanism favoring bat coronavirus cross-species transmission

Bats are reservoirs of various coronaviruses that can jump between bat species or other mammalian hosts, including humans. This article explores coronavirus infection in three bat species (Tadarida brasiliensis, Eumops bonariensis, and Molossus molossus) of the family Molossidae from Argentina using whole viral metagenome analysis. Fecal samples of 47 bats from three semiurban or highly urbanized areas of the province of Santa Fe were investigated. After viral particle enrichment, total RNA was sequenced using the Illumina NextSeq 550 instrument; the reads were assembled into contigs and taxonomically and phylogenetically analyzed. Three novel complete Alphacoronavirus (AlphaCoV) genomes (Tb1-3) and two partial sequences were identified in T. brasiliensis (Tb4-5), and an additional four partial sequences were identified in M. molossus (Mm1-4). Phylogenomic analysis showed that the novel AlphaCoV clustered in two different lineages distinct from the 15 officially recognized AlphaCoV subgenera. Tb2 and Tb3 isolates appeared to be variants of the same virus, probably involved in a persistent infectious cycle within the T. brasiliensis colony. Using recombination analysis, we detected a statistically significant event in Spike gene, which was reinforced by phylogenetic tree incongruence analysis, involving novel Tb1 and AlphaCoVs identified in Eptesicus fuscus (family Vespertilionidae) from the U.S. The putative recombinant region is in the S1 subdomain of the Spike gene, encompassing the potential receptor-binding domain of AlphaCoVs. This study reports the first AlphaCoV genomes in molossids from the Americas and provides new insights into recombination as an important mode of evolution of coronaviruses involved in cross-species transmission. IMPORTANCE This study generated three novel complete AlphaCoV genomes (Tb1, Tb2, and Tb3 isolates) identified in individuals of Tadarida brasiliensis from Argentina, which showed two different evolutionary patterns and are the first to be reported in the family Molossidae in the Americas. The novel Tb1 isolate was found to be involved in a putative recombination event with alphacoronaviruses identified in bats of the genus Eptesicus from the U.S., whereas isolates Tb2 and Tb3 were found in different collection seasons and might be involved in persistent viral infections in the bat colony. These findings contribute to our knowledge of the global diversity of bat coronaviruses in poorly studied species and highlight the different evolutionary aspects of AlphaCoVs circulating in bat populations in Argentina.

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.

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 biology of papillomavirus PDZ associations: what do they offer papillomaviruses?

The high-risk α-type papillomaviruses have a C-terminal PDZ-binding motif (PBM) on one of the two major oncoproteins E6 or E7; the vast majority on E6. The PBM is essential for the high-risk HPV life cycle, for episomal maintenance of the virus genome, and for maintaining the mitotic stability of the infected cell. The question is why only these viruses have PBMs - are there specific constraints imposed by the mucosal epithelium in which these viruses replicate? However the low-risk α-HPVs, such as HPV-6 and HPV-11 replicate extremely efficiently without a PBM, while viruses of the alpha8 group, such as HPV-40, replicate well with a very primitive PBM. So what does PDZ-binding capacity contribute to the fitness of the virus?

The genetic diversity of human papillomavirus types from the species Gammapapillomavirus 15: HPV135, HPV146, and HPV179

OBJECTIVES

To determine the prevalence, viral load, tissue tropism, and genetic variability of novel human papillomavirus (HPV) type 179, which is etiologically associated with sporadic cases of common warts in immunocompromised patients, and phylogenetically related HPV types 135 and 146.

METHODS

The representative collection of 850 HPV-associated clinical samples (oral/nasopharyngeal/anal, archival specimens of oral/oropharyngeal/conjunctival/cervical/skin cancer, benign lesions of the larynx/conjunctiva/skin, and eyebrows), obtained from immunocompetent individuals, was tested for the presence of HPV179, HPV135, and HPV146 using type-specific real-time PCRs. To assess the genetic diversity of the HPVs investigated in the non-coding long control region (LCR), several highly sensitive nested PCR protocols were developed for each HPV type. The genetic diversity of HPV179 was additionally determined in 12 HPV179 isolates from different anatomical sites of an only immunocompromised patient included in the study.

RESULTS

HPV179, HPV135, and HPV146 were detected in 1.4, 2.0, and 1.5% of the samples tested, respectively, with no preference for cutaneous or mucosal epithelial cells. One (with five single nucleotide polymorphisms; SNPs), four (with one to six SNPs), and four (with one to eight SNPs) genetic variants of HPV179, HPV135, and HPV146, respectively, were identified among eligible samples. HPV179 isolates from the immunocompromised patient exhibited the identical LCR nucleotide sequence, suggesting that HPV179 can cause generalized HPV infections.

CONCLUSIONS

HPV179, HPV135, and HPV146 have a mucocutaneous tissue tropism and are associated with sporadic infections in immunocompromised and immunocompetent individuals. Because the majority of mutations were found outside the major functional domains of the respective LCRs, we assume that HPV179, HPV135, and HPV146 genetic variants pathogenically do not differ from their prototypes. In addition, no association was found between specific HPV179, HPV135, and HPV146 genetic variants and anatomical sites of infection and/or specific neoplasms.

Whole Genomic Analysis and Comparison of Two Canine Papillomavirus Type 9 Strains in Malignant and Benign Skin Lesions

Papillomaviruses (PVs) usually cause benign proliferative lesions in the stratified epithelium of various animal species. However, some high-risk types of PVs have been proven to lead to malignant transformations. In dogs, several canine papillomaviruses (CPVs) have been identified in malignant lesions and are suggested as one of the risk factors for the development of squamous cell carcinomas (SCCs). In the present study, the full genomes of two CPV9 strains from recurrent SCCs of Dog 1 and skin viral papilloma (viral plaque) of Dog 2 were sequenced. Alignment of the two CPV9 sequences with the genome of the reference CPV9 strain (accession no. JF800656.1) derived from a solitary pigmented plaque was performed. Compared with the reference strain, a 27 bp in-frame insertion in the E1 gene was identified in both CPV9 strains in this study. In comparison with the CPV9 strains derived from benign lesions, the CPV9 from the SCCs of Dog 1 exhibited a 328 bp deletion at the 3′ end of the E2 and spacer sequence, which encoded a truncated deduced E2 protein and a chimeric E8^E2 protein. However, there was no difference in the mRNA expression levels of viral oncoproteins of E6 and E7 between the two CPV9 cases, suggesting that the oncogenesis of CPV9 for malignant transformation might be different from that of human papillomaviruses. The roles of E2 and E8^E2 deleted CPV9 in the oncogenesis of benign and malignant lesions should be further investigated.

Genomic organization of a Gamma-6 papillomavirus metagenomic discovered from vaginal swab samples of Chinese pregnant women

A complete genome sequence of human papillomaviruses (HPV) named as HPV-ujs-21015 was determined by viral metagenomic and PCR methods. The complete genome is 7354 bp in length with GC content of 41.7%, of which the genome was predicted to contain six ORFs (Open Reading Frame, ORF) coding for four early proteins (E7, E1, E4, and E2) and two late proteins (L1 and L2). Phylogenetic analysis based on the complete genome and the L1 protein showed that HPV-ujs-21015 belongs to a type 214 member within genus Gamma-6 papillomavirus. It is the first complete genome of Gamma-6 papillomavirus discovered from pregnant women in China.

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.

Discovery, characterisation and genomic variation of six novel Gammapapillomavirus types from penile swabs in South Africa

Six novel human papillomaviruses from penile swabs were characterised. Multiple full genome clones for each novel type were generated, and complete genome sizes were: HPV211 (7253bp), HPV212 (7208bp), HPV213 (7096bp), HPV214 (7357), HPV215 (7186bp) and HPV216 (7233bp). Phylogenetically the novel papillomaviruses all clustered with Gammapapillomaviruses: HPV211 is most closely related to HPV168 (72% identity in the L1 nucleotide sequence) of the Gamma-8 species, HPV212 is most closely related to HPV144 (82.9%) of the Gamma-17 species, HPV213 is most closely related to HPV153 (71.8%) of the Gamma-13 species, HPV214 is most closely related to HPV103 (75.3%) of the Gamma-6 species, HPV215 and HPV216 are most closely related to HPV129 (76.8% and 79.2% respectively) of the Gamma-9 species. The novel HPV types demonstrated the classical genomic organisation of Gammapapillomavirusess, with seven open reading frames (ORFs) encoding five early (E1, E2, E4, E6 and E7) and two late (L1 and L2) proteins. Typical of Gammapapillomavirusess the novel types all lacked the E5 ORF and HPV214 also lacked the E6 ORF. HPV212 had nine unique variants, HPV213 had five and HPV215 had four variants. Conserved domains observed among the novel types are the Zinc finger Binding Domain and PDZ domains. A retinoblastoma binding domain (pRB) binding domain in E7 protein was additionally identified in HPV214. This study expands the knowledge of the rapidly growing Gammapapillomavirus genus.

High prevalence of Gammapapillomaviruses (Gamma-PVs) in pre-malignant cutaneous lesions of immunocompetent individuals using a new broad-spectrum primer system, and identification of HPV210, a novel Gamma-PV type

Genus Gammapapillomavirus (Gamma-PV) is the most diverse and largest clade within the Papillomaviridae family. A novel set of degenerate primers targeting the E1 gene was designed and further used in combination with the well-known CUT PCR assay to assess HPV prevalence and genus distribution in a variety of cutaneous samples from 448 immunocompetent individuals. General HPV, Gamma-PV and mixed infections prevalence were significantly higher in actinic keratosis with respect to benign and malignant neoplasms, respectively (p = 0.0047, p = 0.0172, p = 0.00001). Gamma-PVs were significantly more common in actinic keratosis biopsies than Beta- and Alpha-PVs (p = 0.002). The full-length genome sequence of a novel putative Gamma-PV type was amplified by 'hanging droplet' long-range PCR and cloned. The novel virus, designated HPV210, clustered within species Gamma-12. This study provides an additional tool enabling detection of HPV infections in skin and adds new insights about possible early roles of Gamma-PVs in the development of cutaneous malignant lesions.

Complete Genome Sequences of Four Novel Human Gammapapillomavirus Types, HPV-219, HPV-220, HPV-221, and HPV-222, Isolated from Penile Skin Swabs from South African Men

Four novel human gammapapillomaviruses were characterized from penile specimens using genome amplification, cloning, and sequencing. The HPV-219 L1 gene showed 87% nucleotide identity to that of HPV-213 of species gamma-13, HPV-220 had 72% identity to L1 of HPV-212 (gamma-17), HPV-221 had 80% identity to L1 of HPV-142 (gamma-10), and HPV-222 had 73% nucleotide identity to L1 of HPV-162 (gamma-19).

Four novel human gammapapillomaviruses were characterized from penile specimens using genome amplification, cloning, and sequencing. The HPV-219 L1 gene showed 87% nucleotide identity to that of HPV-213 of species gamma-13, HPV-220 had 72% identity to L1 of HPV-212 (gamma-17), HPV-221 had 80% identity to L1 of HPV-142 (gamma-10), and HPV-222 had 73% nucleotide identity to L1 of HPV-162 (gamma-19).

Identification and characterization of two novel Gammapapillomavirus genomes in skin of an immunosuppressed Epidermodysplasia Verruciformis patient

Two novel human gamma-papillomavirus genomes (HPV_MTS3, and HPV_MTS4) were isolated from the skin of an immunosuppressed, late-onset Epidermodysplasia Verruciformis patient and fully cloned. The L1 open reading frames of HPV_MTS3 and HPV_MTS4 were 77% and 91% identical to their closest HPV full genome isolates w18c39 and EV03c60, which belong to the species gamma-22and gamma-7 of the genus Gammapapillomavirus, respectively.

Mouse papillomavirus infections spread to cutaneous sites with progression to malignancy

We report secondary cutaneous infections in the mouse papillomavirus (MmuPV1)/mouse model. Our previous study demonstrated that cutaneous MmuPV1 infection could spread to mucosal sites. Recently, we observed that mucosal infections could also spread to various cutaneous sites including the back, tail, muzzle and mammary tissues. The secondary site lesions were positive for viral DNA, viral capsid protein and viral particles as determined by in situ hybridization, immunohistochemistry and transmission electron microscopy analyses, respectively. We also demonstrated differential viral production and tumour growth at different secondarily infected skin sites. For example, fewer viral particles were detected in the least susceptible back tissues when compared with those in the infected muzzle and tail, although similar amounts of viral DNA were detected. Follow-up studies demonstrated that significantly lower amounts of viral DNA were packaged in the back lesions. Lavages harvested from the oral cavity and lower genital tracts were equally infectious at both cutaneous and mucosal sites, supporting the broad tissue tropism of this papillomavirus. Importantly, two secondary skin lesions on the forearms of two mice displayed a malignant phenotype at about 9.5 months post-primary infection. Therefore, MmuPV1 induces not only dysplasia at mucosal sites such as the vagina, anus and oral cavity but also skin carcinoma at cutaneous sites. These findings demonstrate that MmuPV1 mucosal infection can be spread to cutaneous sites and suggest that the model could serve a useful role in the study of the viral life cycle and pathogenesis of papillomavirus.

Natural history of human papillomavirus infection of sun-exposed healthy skin of immunocompetent individuals over three climatic seasons and identification of HPV209, a novel betapapillomavirus

We present the first longitudinal study reporting the natural history of human papillomavirus (HPV) infection in sun-exposed skin of healthy individuals living in a geographical area in which solar UV radiation is influenced by the ozone content of the atmosphere. During three climatic seasons, skin swab samples were obtained from 78 healthy individuals and the prevalence of cutaneous HPVs was assessed with broad-spectrum FAP and CUT primers and determined at 54, 45 and 47 % in spring, summer and winter, respectively. Frequencies of mixed HPV infections were significantly higher in spring with respect to summer and winter (P=0.02). Seventy-one different HPV types/putative types were identified. While 62 volunteers were HPV-infected in at least one season, 23 had persistent infections. β-PVs (β-1) were the most prevalent and persistent. Age was associated with both the infection status (P=0.01) and the type of HPV infection (no infection, indeterminate/transient, persistent P=0.02). The molecular/phylogenetic analysis of the newly identified β-PV, officially designated as HPV209, showed that the virus has a typical genomic organization of cutaneous HPVs with five early (E6, E7, E1, E2 and E4) and two late genes (L2 and L1), which clusters to the species β-2. This provides useful data on cutaneous HPV infections in high UV-exposed regions.

Molecular characterization, tissue tropism, and genetic variability of the novel Mupapillomavirus type HPV204 and phylogenetically related types HPV1 and HPV63

HPV204 is the only newly identified Mupapillomavirus (Mu-PV) type in more than a decade. To comprehensively characterize HPV204, we performed a detailed molecular analysis of the viral genome and evaluated its clinical relevance in comparison to the other Mu-PVs, HPV1 and HPV63. The 7,227-bp long genome of HPV204 exhibits typical genomic organization of Mu-PVs with eight open reading frames (ORFs) (E6, E7, E1, E2, E8, E4, L2, and L1). We developed three type-specific quantitative real-time PCRs and used them to test a representative collection (n = 1,006) of various HPV-associated benign and malignant neoplasms, as well as samples of clinically normal cutaneous, mucosal, and mucocutaneous origins. HPV204, HPV1, and HPV63 were detected in 1.1%, 2.7%, and 1.9% of samples tested, respectively, and were present in skin and mucosa, suggesting dual tissue tropism of all Mu-PVs. To evaluate the etiological role of Mu-PVs in the development of HPV-associated neoplasms, Mu-PV viral loads per single cell were estimated. HPV1 and HPV63 were present in high viral copy numbers in 3/43 and 1/43 cutaneous warts, respectively, and were identified as the most likely causative agents of these warts. HPV204 viral load was extremely low in a single HPV204-positive cutaneous wart (7.4 × 10⁻⁷ viral copies/cell). Hence, etiological association between HPV204 and the development of cutaneous warts could not be established. To the best of our knowledge, this is the first study to evaluate the genetic variability of Mu-PVs by sequencing complete LCR genomic regions of HPV204, HPV1, and HPV63. We detected several nucleotide substitutions and deletions within the LCR genomic regions of Mu-PVs and identified two genetic variants of HPV204 and HPV63 and five genetic variants of HPV1.