High throughput sequencing reveals diversity of Human Papillomaviruses in cutaneous lesions

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.

Comparison of four different human papillomavirus genotyping methods in cervical samples: Addressing method-specific advantages and limitations

Since human papillomavirus (HPV) is recognized as the causative agent of cervical cancer and associated with anogenital non-cervical and oropharyngeal cancers, the characterization of the HPV types circulating in different geographic regions is an important tool in screening and prevention. In this context, this study compared four methodologies for HPV detection and genotyping: real-time PCR (Cobas® HPV test), nested PCR followed by conventional Sanger sequencing, reverse hybridization (High + Low PapillomaStrip® kit) and next-generation sequencing (NGS) at an Illumina HiSeq2500 platform. Cervical samples from patients followed at the Family Health Strategy from Juiz de Fora, Minas Gerais, Brazil, were collected and subjected to the real-time PCR. Of those, 114 were included in this study according to the results obtained with the real-time PCR, considered herein as the gold standard method. For the 110 samples tested by at least one methodology in addition to real-time PCR, NGS showed the lowest concordance rates of HPV and high-risk HPV identification compared to the other three methods (67-75 %). Real-time PCR and Sanger sequencing showed the highest rates of concordance (97-100 %). All methods differed in their sensitivity and specificity. HPV genotyping contributes to individual risk stratification, therapeutic decisions, epidemiological studies and vaccine development, supporting approaches in prevention, healthcare and management of HPV infection.

Prevalence of novel gamma HPV types 223 and 225 in oral cavity and skin of Indian normal and neoplastic participants

Gamma-papillomaviruses, though traditionally classified as cutaneotropic, actual tissue tropism is largely unexplored. This study aimed to evaluate the tissue-specific prevalence of two novel-HPV 223 and 225 in samples of oral mucosa and keratinized epithelium of varied skin parts from 226 female and male subjects, with or without neoplastic/dysplastic lesions in oral cavity or cervix. The gamma-human papillomavirus (gamma-HPV) 223 and 225 DNA presences were determined by polymerase chain reaction (PCR) ursing the HPV type-specific primers and confirmed by Sanger sequencing. Viral load in the HPV 223 and HPV 225 positive samples were determined by absolute real-time quantification method. Alpha-HPV DNA prevalence was also checked in oral mucosa to ascertain coinfection status. Novel HPV 223 was present in 4.4% (10/226) oral mucosal samples of the study population; interestingly all were females with no prevalence in their corresponding skin swab samples. Whereas, the prevalence of HPV 225 was found both in the skin and oral mucosa of 28.2% (N = 37/131) female and 17.9% (N = 17/95) male participants. Alongside, HPV 223 viral load was found to be significantly higher (p = 0.02 < 0.05) in the oral mucosa of diseased participants, whereas, HPV 225 viral load was higher in the oral mucosa of normal participants. Our results suggest that gamma-HPV 223 has its prevalence only in the oral mucosal epithelium, whereas, HPV 225 has its prevalence on both mucosal and keratinized skin epithelium, indicating its dual tropism nature.

The role of HPV in keratinocyte skin cancer development: A systematic review

Keratinocyte skin cancers are the most frequent malignancy, accounting for approximately 30% of all cancers. Although beta genus HPV are the main etiologic agents for squamous cell carcinoma development in patients with epidermodysplasia verruciformis and organ transplant recipients, their role in non-melanoma skin cancer (NMSC) progression in the general population remains controversial. The aim of our review is to summarize current scientific data and to systematically analyse evidence regarding the role of HPV in keratinocyte skin cancers. A total of 2284 patients were included, of which 724 with actinic keratoses, 290 with Bowen's disease, 949 with cutaneous squamous cell carcinomas and 321 with keratoacanthomas. In the case of actinic keratoses, the majority were positive for beta (n = 372, 58.49%) and gamma HPV (n = 256, 40.25%) and only a few (n = 6, 0.94%) were positive for alpha subtypes. Similarly, most of the cutaneous squamous cell carcinomas were positive for beta (n = 248, 55.98%) and gamma HPV (n = 172, 33.82%) and 23 cases (2.42%) were positive for alpha subtypes. Bowen's disease lesions were mostly positive for beta (n = 43, 55.84%) and alpha HPV (n = 30, 38.96%), in contrast to the gamma genus (n = 4, 5.19%). Keratoacanthomas showed a high distribution among beta genus (n = 79, 50.31%) and an equal proportion between alpha (n = 39, 24.84%) and gamma (n = 39, 24.84%) genera. Studies published so far identifying HPV in keratinocyte skin cancers reflect the difference in detection methods rather than a type-specific tendency towards either actinic keratoses, Bowen's disease, squamous cell carcinoma or keratoacanthoma. On the other hand, recent evidence regarding the role of HPV vaccination in patients with non-melanoma skin cancer brings into perspective the idea of a beta-HPV vaccine or a combined alpha and beta-HPV vaccine that could be used as an adjuvant treatment measure in patients with recalcitrant non-melanoma skin cancer.

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.

MinION nanopore sequencing and assembly of a complete human papillomavirus genome

BACKGROUND

The MinION sequencer belongs to the third generation of sequencing technology that allows for the generation of ultra-long reads, representing a potentially more effective approach to characterize entire viral genome sequences than other time-consuming and low-throughput methodologies.

METHODS

We report the use of the MinION nanopore sequencer to sequence the full-length genome of human papillomavirus (HPV)-ICB2 (7441 bp), which was previously characterized in our laboratory. Three independent MinION libraries were prepared and sequenced using either three consecutive 12 -h runs (Protocol A) or a single run of 48 h starting from a pool of three barcoded DNA libraries (Protocol B). A fully automated bioinformatics pipeline was developed for the reconstruction of the viral genome.

RESULTS

Protocols A and B generated 9,354,933 and 3,255,879 reads, respectively. Read length N50 values ranged between 6976 and 7360 nucleotides over the four sequencing runs. Bioinformatics analysis showed that both protocols allowed for the reconstruction of the whole viral genome, with pairwise percentages of identity to HPV-ICB2 of 100 % for protocol A and 99.98 % for protocol B.

CONCLUSION

Our results show that the use of the MinION nanopore sequencer represents an effective strategy for whole-genome sequencing of HPVs with a minimal error rate.

Explainable deep neural networks for novel viral genome prediction

Viral infection causes a wide variety of human diseases including cancer and COVID-19. Viruses invade host cells and associate with host molecules, potentially disrupting the normal function of hosts that leads to fatal diseases. Novel viral genome prediction is crucial for understanding the complex viral diseases like AIDS and Ebola. While most existing computational techniques classify viral genomes, the efficiency of the classification depends solely on the structural features extracted. The state-of-the-art DNN models achieved excellent performance by automatic extraction of classification features, but the degree of model explainability is relatively poor. During model training for viral prediction, proposed CNN, CNN-LSTM based methods (EdeepVPP, EdeepVPP-hybrid) automatically extracts features. EdeepVPP also performs model interpretability in order to extract the most important patterns that cause viral genomes through learned filters. It is an interpretable CNN model that extracts vital biologically relevant patterns (features) from feature maps of viral sequences. The EdeepVPP-hybrid predictor outperforms all the existing methods by achieving 0.992 mean AUC-ROC and 0.990 AUC-PR on 19 human metagenomic contig experiment datasets using 10-fold cross-validation. We evaluate the ability of CNN filters to detect patterns across high average activation values. To further asses the robustness of EdeepVPP model, we perform leave-one-experiment-out cross-validation. It can work as a recommendation system to further analyze the raw sequences labeled as ‘unknown’ by alignment-based methods. We show that our interpretable model can extract patterns that are considered to be the most important features for predicting virus sequences through learned filters.

Transcription of human papillomaviruses in nonmelanoma skin cancers of the immunosuppressed

Nonmelanoma skin cancer (NMSC) has a greatly increased incidence among the immunosuppressed and the DNA of human papillomavirus (HPV) is commonly found in these tumors. To investigate if there are any actively transcribed HPV infections in these tumors, we identified all skin cancers diagnosed after solid organ transplantation in Sweden during 1964-2011 (n = 7614 NMSCs) and requested the diagnostic tumor blocks from the corresponding pathology archives. For the present study, we selected diagnostic specimens from 345 NMSC and performed whole genome transcriptome analysis using NovaSeq (Illumina), in comparison with three cervical cancers. Although we obtained an abundance of high-quality paired reads per sample (median of 35 million reads), only 15 NMSC specimens contained HPV transcription. Three specimens had transcription of oncogenic anogenital HPVs (HPV16 and 56), six tumors had transcription of HPVs from the beta-2 species (three HPV38, two with HPV23 and one with HPV107) and then there was one observation each of transcription of HPVs 3, 26, 57, 147, 158, 168 and of two nonestablished HPV types belonging to the gamma genus. In conclusion, transcription of specific HPV types can be found in NMSC among the immunosuppressed, but this is not common.

Case control study comparing the HPV genome in patients with oral cavity squamous cell carcinoma to normal patients using metagenomic shotgun sequencing

The aim of this study was to carry out a case control study comparing the HPV genome in patients with oral cavity squamous cell carcinoma (OC-SCC) to normal patients using metagenomic shotgun sequencing. We recruited 50 OC-SCC cases which were then matched with a control patient by age, gender, race, smoking status and alcohol status. DNA was extracted from oral wash samples from all patients and whole genome shotgun sequencing performed. The raw sequence data was cleaned, reads aligned with the human genome (GRCH38), nonhuman reads identified and then HPV genotypes identified using HPViewer. In the 50 patients with OC-SCC, the most common subsite was tongue in 26 (52%). All patients were treated with primary resection and neck dissection. All but 2 tumors were negative on p16 immunohistochemistry. There were no statistically significant differences between the cases and controls in terms of gender, age, race/ethnicity, alcohol drinking, and cigarette smoking. There was no statistically significant difference between the cancer samples and control samples in the nonhuman DNA reads (medians 4,228,072 vs. 5,719,715, P value = 0.324). HPV was detected in 5 cases (10%) of OC-SCC (genotypes 10, 16, 98) but only 1 tumor sample (genotype 16) yielded a high number of reads to suggest a role in the etiology of OC-SCC. HPV was detected in 4 control patients (genotypes 16, 22, 76, 200) but all had only 1-2 HPV reads per human genome. Genotypes of HPV are rarely found in patients with oral cancer.

Immunodiagnosis and Immunotherapeutics Based on Human Papillomavirus for HPV-Induced Cancers

Infection with human papillomavirus (HPV) is one of the main causes of malignant neoplasms, especially cervical, anogenital, and oropharyngeal cancers. Although we have developed preventive vaccines that can protect from HPV infection, there are still many new cases of HPV-related cancers worldwide. Early diagnosis and therapy are therefore important for the treatment of these diseases. As HPVs are the major contributors to these cancers, it is reasonable to develop reagents, kits, or devices to detect and eliminate HPVs for early diagnosis and therapeutics. Immunological methods are precise strategies that are promising for the accurate detection and blockade of HPVs. During the last decades, the mechanism of how HPVs induce neoplasms has been extensively elucidated, and several oncogenic HPV early proteins, including E5, E6, and E7, have been shown to be positively related to the oncogenesis and malignancy of HPV-induced cancers. These oncoproteins are promising biomarkers for diagnosis and as targets for the therapeutics of HPV-related cancers. Importantly, many specific monoclonal antibodies (mAbs), or newly designed antibody mimics, as well as new immunological kits, devices, and reagents have been developed for both the immunodiagnosis and immunotherapeutics of HPV-induced cancers. In the current review, we summarize the research progress in the immunodiagnosis and immunotherapeutics based on HPV for HPV-induced cancers. In particular, we depict the most promising serological methods for the detection of HPV infection and several therapeutical immunotherapeutics based on HPV, using immunological tools, including native mAbs, radio-labelled mAbs, affitoxins (affibody-linked toxins), intracellular single-chain antibodies (scFvs), nanobodies, therapeutical vaccines, and T-cell-based therapies. Our review aims to provide new clues for researchers to develop novel strategies and methods for the diagnosis and treatment of HPV-induced tumors.

Detection of human papillomaviruses in paired healthy skin and actinic keratosis by next generation sequencing

Actinic keratosis (AK) arises on photo-damaged skin and is considered to be the precursor lesion of cutaneous squamous cell carcinoma (cSCC). Many findings support the involvement of β human papillomaviruses (HPVs) in cSCC, while very little is known on γ HPV types. The objective of this study was to characterize the spectrum of PV types in healthy skin (HS) and AK samples of the same immunocompetent individuals using next generation sequencing (NGS). Viral DNA of 244 AK and 242 HS specimens were amplified by PCR using two different sets of primers (FAP59/64 and FAPM1). Purified amplicons were pooled and sequenced using NGS. The study resulted in the identification of a large number of known β and γ PV types. In addition, 27 putative novel β and 16 γ and 4 unclassified PVs were isolated. HPV types of species γ-1 (e.g. HPV4) appeared to be strongly enriched in AK versus HS. The NGS analysis revealed that a large spectrum of known and novel PVs is present in HS and AK. The evidence that species γ-1 HPV types appears to be enriched in AK in comparison to HS warrants further studies to evaluate their role in development of skin (pre)cancerous lesions.

Perspectives in HPV Secondary Screening and Personalized Therapy Basing on Our Understanding of HPV-Related Carcinogenesis Pathways

As cervical cancer is one of the most common malignancies in women worldwide even with present screening methods, the incidence in most developed countries is not decreasing for the last 15-20 years. A shift has been observed in the age of diagnosis in favour of younger women, and treatment of already developed cervical cancer is a challenge for surgeons. It is imperative to find new diagnostic methods for accurately pointing out patients at high risk of developing malignant disease and developing personalized treatment. Since cervical cancer is almost exclusively associated with HPV infection, understanding changes happening in an infected cell may prove invaluable for search of such methods, but it may also prove helpful in the diagnosis and treatment of other anogenital and nasopharyngeal region cancers. This review follows HPV-related changes in infected cell biology to point what potential markers and targets for therapy are in option when dealing with HPV-related diseases.

Emerging biomarkers and clinical significance of HPV genotyping in prevention and management of cervical cancer

Cervical cancer is a growing and serious problem world-wide in women, but more acute in developing countries especially in Indian subcontinent. The main causative agent for the disease is Human Papilloma Virus (HPV). The history of the cervical cancer goes back to eighteenth century as the HPV infection is reported since 1800s. Presently, the genetic structure of HPV is well defined. Several screening tests including cytology and visual based screening and high risk HPV testing are available. Also available are various clinical and commercial diagnostic tests. However due to the lack of awareness and population-based screening programs, the morbidity and mortality rate is alarmingly high. There are new emerging biomarkers including E6/E7 mRNA, p16^ink4a, markers of aberrant S-phase induction, chromosomal abnormalities and miRNAs along with advanced genotyping methods. These markers have clinical significance and are helpful in disease prevention and management. Further, recent advancement in the field of metagenomics has increased the prospects of identifying newer microbes, viruses hitherto reported thus far in the context of HPV infection. Analysis of HPV cases using modern tools including genotyping using more powerful biomarkers is envisaged to enhance the prospects of early diagnosis, better prognosis, more reliable treatment and eventual management of the disease.

ViraMiner: Deep learning on raw DNA sequences for identifying viral genomes in human samples

Despite its clinical importance, detection of highly divergent or yet unknown viruses is a major challenge. When human samples are sequenced, conventional alignments classify many assembled contigs as "unknown" since many of the sequences are not similar to known genomes. In this work, we developed ViraMiner, a deep learning-based method to identify viruses in various human biospecimens. ViraMiner contains two branches of Convolutional Neural Networks designed to detect both patterns and pattern-frequencies on raw metagenomics contigs. The training dataset included sequences obtained from 19 metagenomic experiments which were analyzed and labeled by BLAST. The model achieves significantly improved accuracy compared to other machine learning methods for viral genome classification. Using 300 bp contigs ViraMiner achieves 0.923 area under the ROC curve. To our knowledge, this is the first machine learning methodology that can detect the presence of viral sequences among raw metagenomic contigs from diverse human samples. We suggest that the proposed model captures different types of information of genome composition, and can be used as a recommendation system to further investigate sequences labeled as "unknown" by conventional alignment methods. Exploring these highly-divergent viruses, in turn, can enhance our knowledge of infectious causes of diseases.

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.

Machine Learning for detection of viral sequences in human metagenomic datasets

BACKGROUND

Detection of highly divergent or yet unknown viruses from metagenomics sequencing datasets is a major bioinformatics challenge. When human samples are sequenced, a large proportion of assembled contigs are classified as "unknown", as conventional methods find no similarity to known sequences. We wished to explore whether machine learning algorithms using Relative Synonymous Codon Usage frequency (RSCU) could improve the detection of viral sequences in metagenomic sequencing data.

RESULTS

We trained Random Forest and Artificial Neural Network using metagenomic sequences taxonomically classified into virus and non-virus classes. The algorithms achieved accuracies well beyond chance level, with area under ROC curve 0.79. Two codons (TCG and CGC) were found to have a particularly strong discriminative capacity.

CONCLUSION

RSCU-based machine learning techniques applied to metagenomic sequencing data can help identify a large number of putative viral sequences and provide an addition to conventional methods for taxonomic classification.

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.

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.

Extension of the viral ecology in humans using viral profile hidden Markov models

When human samples are sequenced, many assembled contigs are "unknown", as conventional alignments find no similarity to known sequences. Hidden Markov models (HMM) exploit the positions of specific nucleotides in protein-encoding codons in various microbes. The algorithm HMMER3 implements HMM using a reference set of sequences encoding viral proteins, "vFam". We used HMMER3 analysis of "unknown" human sample-derived sequences and identified 510 contigs distantly related to viruses (Anelloviridae (n = 1), Baculoviridae (n = 34), Circoviridae (n = 35), Caulimoviridae (n = 3), Closteroviridae (n = 5), Geminiviridae (n = 21), Herpesviridae (n = 10), Iridoviridae (n = 12), Marseillevirus (n = 26), Mimiviridae (n = 80), Phycodnaviridae (n = 165), Poxviridae (n = 23), Retroviridae (n = 6) and 89 contigs related to described viruses not yet assigned to any taxonomic family). In summary, we find that analysis using the HMMER3 algorithm and the "vFam" database greatly extended the detection of viruses in biospecimens from humans.

Detection of novel Betapapillomaviruses and Gammapapillomaviruses in eyebrow hair follicles using a single-tube 'hanging droplet' PCR assay with modified pan-PV CODEHOP primers

A modified pan-PV consensus-degenerate hybrid oligonucleotide primer (CODEHOP) PCR was developed for generic and sensitive detection of a broad-spectrum of human papillomaviruses (HPVs) infecting the cutaneous epithelium. To test the analytical sensitivity of the assay we examined 149 eyebrow hair follicle specimens from immunocompetent male patients. HPV DNA was detected in 60 % (89/149) of analysed eyebrow samples with a total of 48 different HPV sequences, representing 21 previously described HPVs and 27 putative novel HPV types. Evidence for ten novel HPV subtypes and seven viral variants, clustering to three out of five genera containing cutaneous HPVs, was also obtained. Thus, we have shown that the modified pan-PV CODEHOP PCR assay is able to identify multiple HPV types, even from different genera, in the same clinical sample. Overall, these results demonstrate that the pan-PV CODEHOP PCR is an excellent tool for screening and identification of novel cutaneous HPVs, even in samples with low viral loads.

Genotyping of high-risk anal human papillomavirus (HPV): ion torrent-next generation sequencing vs. linear array

Background

Our next generation sequencing (NGS)-based human papillomavirus (HPV) genotyping assay showed a high degree of concordance with the Roche Linear Array (LA) with as little as 1.25 ng formalin-fixed paraffin-embedded-derived genomic DNA in head and neck and cervical cancer samples. This sensitive genotyping assay uses barcoded HPV PCR broad-spectrum general primers 5+/6+ (BSGP)5+/6+ applicable to population studies, but it’s diagnostic performance has not been tested in cases with multiple concurrent HPV infections.

Methods

We conducted a cross-sectional study to compare the positive and negative predictive value (PPV and NPV), sensitivity and specificity of the NGS assay to detect HPV genotype infections as compared to the LA. DNA was previously extracted from ten anal swab samples from men who have sex with men in Nigeria enrolled on the TRUST/RV368 cohort study. Two-sample tests of proportions were used to examine differences in the diagnostic performance of the NGS assay to detect high vs. low-risk HPV type-specific infections.

Results

In total there were 94 type-specific infections detected in 10 samples with a median of 9.5, range (9 to 10) per sample. Using the LA as the gold standard, 84.4% (95% CI: 75.2–91.2) of the same anal type-specific infections detected on the NGS assay had been detected by LA. The PPV and sensitivity differed significantly for high risk (PPV: 90%, 95% CI: 79.5–96.2; sensitivity: 93.1%, 95% CI: 83.3–98.1) as compared to low risk HPV (PPV: 73%, 95% CI: 54.1–87.7; sensitivity: 61.1, 95% CI: 43.5–76.9) (all p < 0.05). The NPV for all types was 92.5% (95% CI: 88.4–95.4). The NPV and specificity were similar for high and low risk HPVs (all p > 0.05). The NGS assay detected 10 HPV genotypes that were not among the 37 genotypes found on LA (30, 32, 43, 44, 74, 86, 87, 90, 91, 114).

Conclusions

The NGS assay accurately detects multiple HPV infections in individual clinical specimens with limited sample volume and has extended coverage compared to LA.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0771-z) contains supplementary material, which is available to authorized users.

Genome Sequence of a Novel Human Gammapapillomavirus Isolated from Skin

A new human gammapapillomavirus (HPV_MTS2) genome was isolated and fully cloned from a skin swab. The L1 open reading frame of HPV_MTS2 was 79% and 80% identical to those of its closest relatives, HPV type 149 (species Gamma-7 of the genus Gammapapillomavirus) and HPV isolate Dysk2 (GenBank accession no. KX781281), respectively, thus qualifying it as a new HPV type.

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.

High diversity of alpha, beta and gamma human papillomaviruses in genital samples from HIV-negative and HIV-positive heterosexual South African men

BACKGROUND

This study investigated the prevalence of HPVs in heterosexual South African men and the impact of HIV co-infection.

METHODS

HPV was detected in penile swabs from 195 HIV-infected and 140 HIV-uninfected men using PCR with FAP59/64 primers and Roche Linear Array HPV genotyping (LA). Genotyping of FAP positive specimens was achieved by high-throughput sequencing of amplicons.

RESULTS

HPV was detected by FAP PCR and LA in 79% (266/335) of the men. Men with HIV co-infection and men with HIV infected sexual partners had a significantly (p<0.0001) higher HPV infection risk (adjusted odds ratio 4.0 (2.1-8.2) and 3.7 (2.1-6.7), respectively). LA genotyping and 454 sequencing of 218 FAP positive specimens detected 45 known α-HPV types, 45 β-HPV types (34 known, 10 putative and 1 novel putative), and 91 γ-HPV types (26 known, 51 putative and 14 novel putative). Alpha, beta and gamma types were detected in 89.8%, 51.4% and 62.4% of the 218 men with HPV-62, HPV-5 and HPV-121 most common in each genus, respectively.

CONCLUSION

A great diversity of known and novel alpha, beta and gamma HPV types were detected with higher prevalence in HIV co-infected men and unknown associations, if any, with genital lesions and cancers.

Molecular tests potentially improving HPV screening and genotyping for cervical cancer prevention

INTRODUCTION

Human papillomavirus (HPV)-related cancers can be averted by type-specific vaccination (primary prevention) and/or through detection and ablation of precancerous cervical lesions (secondary prevention). This review presents current challenges to cervical cancer screening programs, focusing on recent molecular advances in HPV testing and potential improvements on risk stratification. Areas covered: High-risk (HR)-HPV DNA detection has been progressively incorporated into cervix cancer prevention programs based on its increased sensitivity. Advances in next-generation sequencing (NGS) are being rapidly applied to HPV typing. However, current HPV DNA tests lack specificity for identification of cervical precancer (CIN3). HPV typing methods were reviewed based on published literature, with a focus on these applications for screening and risk stratification in the emerging complex clinical scenario post-vaccine introduction. In addition, the potential for NGS technologies to increase specificity is discussed in regards to reflex testing of specimens for emerging biomarkers for cervix precancer/cancer. Expert commentary: Integrative multi-disciplinary molecular tests accurately triaging exfoliated cervical specimens will improve cervical cancer prevention programs while simplifying healthcare procedures in HPV-infected women. Hence, the concept of a 'liquid-biopsy' (i.e., 'molecular' Pap test) highly specific for early identification of cervical precancerous lesions is of critical importance in the years to come.

Seroprevalence of cutaneous human papillomaviruses (HPVs) among men in the multinational HPV Infection in Men study

Data on cutaneous human papillomavirus (HPV) seroprevalence are primarily derived from skin cancer case-control studies. Few studies have reported the seroprevalence of cutaneous HPV among healthy men. This study investigated the seroprevalence of cutaneous HPV types and associated risk factors among men residing in Brazil, Mexico and the USA. Six hundred men were randomly selected from the HPV Infection in Men study. Archived serum specimens were tested for antibodies against 14 cutaneous HPV genotypes, β-HPV types (5/8/12/14/17/22/23/24/38/48), α-HPV 27, γ-HPV 4, µ-HPV1 and ν-HPV 41 using a glutathione S-transferase L1-based multiplex serology assay. Risk factor data were collected by a questionnaire. Binomial proportions were used to estimate seroprevalence, and logistic regression to examine factors associated with seropositivity. Overall, 65.4 % of men were seropositive to ≥1 of the 14 cutaneous HPV types, and 39.0 % were positive for ≥1 β-HPV types. Seroprevalence was 8.9, 30.9, 28.6 and 9.4 % for α-HPV 27, γ-HPV 4, µ-HPV 1 and ν-HPV 41, respectively. In multivariate analyses, seropositivity for any cutaneous HPV type was associated with higher education [adjusted odds ratio (AOR) 1.75; 95 % confidence interval (CI) 1.08-2.83], and seropositivity of any β-HPV type was significantly associated with increasing age (AOR 1.72; 95 % CI 1.12-2.63, for men aged 31-44 years vs men aged 18-30 years). Other factors associated with various type-specific cutaneous HPV seropositivity included country, circumcision and lifetime number of male sexual partners. These data indicate that exposure to cutaneous HPV is common. Future studies are needed to assess the role of cutaneous HPV in diseases.

Identification of a novel human papillomavirus by metagenomic analysis of vaginal swab samples from pregnant women

BACKGROUND

The number of members in the genus Gammapapillomavirus of Family Papillomaviridae has recently been expanding most rapidly. The aim of this study was to characterize a novel human gammapapillomavirus type identified in a vaginal swab from a 25-year-old pregnant woman suffering from vaginitis.

METHODS

Viral metagenomics method was used to detect the viral sequences in 88 vaginal swab samples collected from 88 pregnant women with vaginitis. A novel papillomavirus, named HPV-ZJ01 (GenBank no. KX082661), was detected in one sample and its complete genome sequence was amplified by PCR and sequenced by Sanger walking. Phylogenetic analyses based on the complete genome and the L1 protein of HPV-ZJ01 and other representative human papillomaviruses were done, respectively. Further PCR screening was performed in vaginal swabs (n = 135), cervical smears (n = 40) and cervical cancer tissues (n = 40) using nested-PCR primers designed based on HPV-ZJ01 sequence to investigate the prevalence of HPV-ZJ01.

RESULTS

The genome of HPV-ZJ01 is 7,358 bp in length with a GC content of 37.8 %. HPV-ZJ01 was predicted to contain six open reading frames (E6, E7, E1, E2, L2, and L1) and a non-coding long control region (LCR). The genome shared the highest overall similarity to HPV-166, with 70.6 % nucleotide sequence identity while its L1 gene shared the highest nucleotide similarity to HPV-162, with 71.1 % sequence identity. Phylogenetic analysis suggested that HPV-ZJ01 belongs to a novel HPV type in the Gamma-PV genus, species Gamma-19, already containing HPV161, HPV162 and HPV166. PCR screening results indicated that none of the other samples were positive for HPV-ZJ01 except the original HPV-ZJ01 positive vaginal swab specimen.

CONCLUSION

The genome sequence of a novel type of species Gamma-19 HPV was characterized. The screening PCR results suggested that HPV-ZJ01 is not associated with any of the cervical cancer samples tested. In order to confirm the prevalence and disease association, if any, for HPV-ZJ01, a further study with different sample types and a larger sample size is needed.

Detection of DNA viruses in prostate cancer

We tested prostatic secretions from men with and without prostate cancer (13 cases and 13 matched controls) or prostatitis (18 cases and 18 matched controls) with metagenomic sequencing. A large number (>200) of viral reads was only detected among four prostate cancer cases (1 patient each positive for Merkel cell polyomavirus, JC polyomavirus and Human Papillomavirus types 89 or 40, respectively). Lower numbers of reads from a large variety of viruses were detected in all patient groups. Our knowledge of the biology of the prostate may be furthered by the fact that DNA viruses are commonly shed from the prostate and can be readily detected by metagenomic sequencing of expressed prostate secretions.

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

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

Characterization of Intra-Type Variants of Oncogenic Human Papillomaviruses by Next-Generation Deep Sequencing of the E6/E7 Region

Different human papillomavirus (HPV) types are characterized by differences in tissue tropism and ability to promote cell proliferation and transformation. In addition, clinical and experimental studies have shown that some genetic variants/lineages of high-risk HPV (HR-HPV) types are characterized by increased oncogenic activity and probability to induce cancer. In this study, we designed and validated a new method based on multiplex PCR-deep sequencing of the E6/E7 region of HR-HPV types to characterize HPV intra-type variants in clinical specimens. Validation experiments demonstrated that this method allowed reliable identification of the different lineages of oncogenic HPV types. Advantages of this method over other published methods were represented by its ability to detect variants of all HR-HPV types in a single reaction, to detect variants of HR-HPV types in clinical specimens with multiple infections, and, being based on sequencing of the full E6/E7 region, to detect amino acid changes in these oncogenes potentially associated with increased transforming activity.

HPV Population Profiling in Healthy Men by Next-Generation Deep Sequencing Coupled with HPV-QUEST

Multiple-type human papillomaviruses (HPV) infection presents a greater risk for persistence in asymptomatic individuals and may accelerate cancer development. To extend the scope of HPV types defined by probe-based assays, multiplexing deep sequencing of HPV L1, coupled with an HPV-QUEST genotyping server and a bioinformatic pipeline, was established and applied to survey the diversity of HPV genotypes among a subset of healthy men from the HPV in Men (HIM) Multinational Study. Twenty-one HPV genotypes (12 high-risk and 9 low-risk) were detected in the genital area from 18 asymptomatic individuals. A single HPV type, either HPV16, HPV6b or HPV83, was detected in 7 individuals, while coinfection by 2 to 5 high-risk and/or low-risk genotypes was identified in the other 11 participants. In two individuals studied for over one year, HPV16 persisted, while fluctuations of coinfecting genotypes occurred. HPV L1 regions were generally identical between query and reference sequences, although nonsynonymous and synonymous nucleotide polymorphisms of HPV16, 18, 31, 35h, 59, 70, 73, cand85, 6b, 62, 81, 83, cand89 or JEB2 L1 genotypes, mostly unidentified by linear array, were evident. Deep sequencing coupled with HPV-QUEST provides efficient and unambiguous classification of HPV genotypes in multiple-type HPV infection in host ecosystems.

Identification of a Novel Human Papillomavirus, Type HPV199, Isolated from a Nasopharynx and Anal Canal, and Complete Genomic Characterization of Papillomavirus Species Gamma-12

The novel human papillomavirus type 199 (HPV199) was initially identified in a nasopharyngeal swab sample obtained from a 25 year-old immunocompetent male. The complete genome of HPV199 is 7,184 bp in length with a GC content of 36.5%. Comparative genomic characterization of HPV199 and its closest relatives showed the classical genomic organization of Gammapapillomaviruses (Gamma-PVs). HPV199 has seven major open reading frames (ORFs), encoding five early (E1, E2, E4, E6, and E7) and two late (L1 and L2) proteins, while lacking the E5 ORF. The long control region (LCR) of 513 bp is located between the L1 and E6 ORFs. Phylogenetic analysis additionally confirmed that HPV-199 clusters into the Gamma-PV genus, species Gamma-12, additionally containing HPV127, HV132, HPV148, HPV165, and three putative HPV types: KC5, CG2 and CG3. HPV199 is most closely related to HPV127 (nucleotide identity 77%). The complete viral genome sequence of additional HPV199 isolate was determined from anal canal swab sample. Two HPV199 complete viral sequences exhibit 99.4% nucleotide identity. To the best of our knowledge, this is the first member of Gamma-PV with complete nucleotide sequences determined from two independent clinical samples. To evaluate the tissue tropism of the novel HPV type, 916 clinical samples were tested using HPV199 type-specific real-time PCR: HPV199 was detected in 2/76 tissue samples of histologically confirmed common warts, 2/108 samples of eyebrow hair follicles, 2/137 anal canal swabs obtained from individuals with clinically evident anal pathology, 4/184 nasopharyngeal swabs and 3/411 cervical swabs obtained from women with normal cervical cytology. Although HPV199 was found in 1.4% of cutaneous and mucosal samples only, it exhibits dual tissue tropism. According to the results of our study and literature data, dual tropism of all Gamma-12 members is highly possible.

Does human papillomavirus-negative condylomata exist?

Condylomata acuminata is caused by human papillomavirus (HPV). PCR with consensus primers will typically detect HPV in >96% of condylomata. Metagenomic sequencing has found that some "HPV-negative" condylomata do indeed contain HPV. We wished to perform a renewed evaluation of the "HPV-negative" condylomata using deeper metagenomics sequencing. Sequencing of whole genome amplified DNA from 40 apparently "HPV-negative" condylomata detected HPV in 37/40 specimens. We found 75 different HPV types, out of which 43 represented novel putative HPV types. Three types were cloned and established as HPV types 200, 201 and 202. Molluscum contagiosum virus was detected in 24 of the 40 samples. In summary, deep sequencing enables detection of HPV in almost all condylomata. "HPV-negative" condylomata might largely be explained by clinical misdiagnosis or the presence of viral variants, distantly related HPV types and/or low viral loads.

Viremia during pregnancy and risk of childhood leukemia and lymphomas in the offspring: Nested case-control study

A possible role for infections of the pregnant mother in the development of childhood acute leukemias and lymphomas has been suggested. However, no specific infectious agent has been identified. Offspring of 74,000 mothers who had serum samples taken during pregnancy and stored in a large-scale biobank were followed up to the age of 15 years (750,000 person years) through over-generation linkages between the biobank files, the Swedish national population and cancer registers to identify incident leukemia/lymphoma cases in the offspring. First-trimester sera from mothers of 47 cases and 47 matched controls were retrieved and analyzed using next generation sequencing. Anelloviruses were the most common viruses detected, found in 37/47 cases and in 40/47 controls, respectively (OR: 0.6, 95% CI: 0.2-1.9). None of the detected viruses was associated with leukemia/lymphoma in the offspring. Viremia during pregnancy was common, but no association with leukemia/lymphoma risk in the offspring was found.

Contribution of Beta-HPV Infection and UV Damage to Rapid-Onset Cutaneous Squamous Cell Carcinoma during BRAF-Inhibition Therapy

PURPOSE

BRAF-inhibition (BRAFi) therapy for advanced melanoma carries a high rate of secondary cutaneous squamous cell carcinoma (cSCC) and risk of other cancers. UV radiation and α-genus human papillomavirus (HPV) are highly associated with SCC, but a novel role for β-genus HPV is suspected in BRAFi-cSCC. Cutaneous β-HPV may act in concert with host and environmental factors in BRAFi-cSCC.

EXPERIMENTAL DESIGN

Primary BRAFi-cSCC tissue DNA isolated from patients receiving vemurafenib or dabrafenib from two cancer centers was analyzed for the presence of cutaneous oncogenic viruses and host genetic mutations. Diagnostic specimens underwent consensus dermatopathology review. Clinical parameters for UV exposure and disease course were statistically analyzed in conjunction with histopathology.

RESULTS

Twenty-nine patients contributed 69 BRAFi-cSCC lesions. BRAFi-cSCC had wart-like features (BRAFi-cSCC-WF) in 22% of specimens. During vemurafenib therapy, BRAFi-cSCC-WF arose 11.6 weeks more rapidly than conventional cSCC when controlled for gender and UV exposure (P value = 0.03). Among all BRAFi-cSCC, β-genus HPV-17, HPV-38, HPV-111 were most frequently isolated, and novel β-HPV genotypes were discovered (CTR, CRT-11, CRT-22). Sequencing revealed 63% of evaluated BRAFi-cSCCs harbored RAS mutations with PIK3CA, CKIT, ALK, and EGFR mutations also detected.

CONCLUSIONS

We examined clinical, histopathologic, viral, and genetic parameters in BRAFi-cSCC demonstrating rapid onset; wart-like histomorphology; β-HPV-17, HPV-38, and HPV-111 infection; UV damage; and novel ALK and CKIT mutations. Discovered β-HPV genotypes expand the spectrum of tumor-associated viruses. These findings enhance our understanding of factors cooperating with BRAF inhibition that accelerate keratinocyte oncogenesis as well as broaden the knowledge base of multifactorial mediators of cancer in general.

Molecular methods for identification and characterization of novel papillomaviruses

Papillomaviruses (PV) are a remarkably heterogeneous family of small DNA viruses that infect a wide variety of vertebrate species and are aetiologically linked with the development of various neoplastic changes of the skin and mucosal epithelia. Based on nucleotide similarity, PVs are hierarchically classified into genera, species and types. Novel human PV (HPV) types are given a unique number only after the whole genome has been cloned and deposited with the International HPV Reference Center. As of 9 March 2015, 200 different HPV types, belonging to 49 species, had been recognized by the International HPV Reference Center. In addition, 131 animal PV types identified from 66 different animal species exist. Recent advances in molecular techniques have resulted in an explosive increase in the identification of novel HPV types and novel subgenomic HPV sequences in the last few years. Among PV genera, the γ-PV genus has been growing most rapidly in recent years with 80 completely sequenced HPV types, followed by α-PV and β-PV genera that have 65 and 51 recognized HPV types, respectively. We reviewed in detail the contemporary molecular methods most often used for identification and characterization of novel PV types, including PCR, rolling circle amplification and next-generation sequencing. Furthermore, we present a short overview of 12 and 10 novel HPV types recently identified in Sweden and Slovenia, respectively. Finally, an update on the International Human Papillomavirus Reference Center is provided.

Characterization of two novel gammapapillomaviruses, HPV179 and HPV184, isolated from common warts of a renal-transplant recipient

Gammapapillomavirus (Gamma-PV) is a diverse and rapidly expanding PV-genus, currently consisting of 76 fully characterized human papillomavirus (HPV) types. In this study, DNA genomes of two novel HPV types, HPV179 and HPV184, obtained from two distinct facial verrucae vulgares specimens of a 64 year-old renal-transplant recipient, were fully cloned, sequenced and characterized. HPV179 and HPV184 genomes comprise 7,228-bp and 7,324-bp, respectively, and contain four early (E1, E2, E6 and E7) and two late genes (L1 and L2); the non-coding region is typically positioned between L1 and E6 genes. Phylogenetic analysis of the L1 nucleotide sequence placed both novel types within the Gamma-PV genus: HPV179 was classified as a novel member of species Gamma-15, additionally containing HPV135 and HPV146, while HPV184 was classified as a single member of a novel species Gamma-25. HPV179 and HPV184 type-specific quantitative real-time PCRs were further developed and used in combination with human beta-globin gene quantitative real-time PCR to determine the prevalence and viral load of the novel types in the patient's facial warts and several follow-up skin specimens, and in a representative collection, a total of 569 samples, of HPV-associated benign and malignant neoplasms, hair follicles and anal and oral mucosa specimens obtained from immunocompetent individuals. HPV179 and HPV184 viral loads in patients' facial warts were estimated to be 2,463 and 3,200 genome copies per single cell, respectively, suggesting their active role in the development of common warts in organ-transplant recipients. In addition, in this particular patient, both novel types had established a persistent infection of the skin for more than four years. Among immunocompetent individuals, HPV179 was further detected in low-copy numbers in a few skin specimens, indicating its cutaneous tissue tropism, while HPV184 was further detected in low-copy numbers in one mucosal and a few skin specimens, suggesting its dual tissue tropism.

Papillomaviruses: Viral evolution, cancer and evolutionary medicine

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

International standardization and classification of human papillomavirus types

Established Human Papillomavirus (HPV) types, up to HPV202, belong to 49 species in five genera. International standardization in classification and quality standards for HPV type designation and detection is ensured by the International HPV Reference Center. The center i) receives clones of potentially novel HPV types, re-clones and re-sequences them. If confirmed, an HPV type number is assigned and posted on www.hpvcenter.se. ii) distributes reference clone samples, for academic research, under Material Transfer Agreements agreed with the originator. iii) provides preliminary checking of whether new sequences represent novel types iv) issues international proficiency panels for HPV genotyping. The rate of HPV type discovery is increasing, probably because of metagenomic sequencing. γ-genus today contains 79HPV types and 27 species, surpassing ∝ and β genera with 65 and 51HPV types, respectively. Regular issuing of proficiency panels based on HPV reference clones has resulted in global improvement of HPV genotyping services.

Prevention of human papillomavirus-associated cancers

The oncogenic, anogenital types of human papillomavirus (HPV) are established as causing about 4.8% of all human cancers worldwide, particularly cervical, anal, vulvar, vaginal, penile, and oropharyngeal cancers. Quantitative knowledge of the HPV type-specific risks for these cancers, as well as for the different cervical cancer precursors (cervical intraepithelial neoplasias, CINs), is useful for estimating the effect of elimination of specific HPV types and clinical benefits of screening for specific HPV types. The present review summarizes both the worldwide presence of specific HPV types in cervical cancer precursors and in invasive cervical cancers, and also the long-term follow-up data from a large randomized clinical trial of HPV-based cervical cancer screening. All 12 HPV types classified as class I (established) carcinogens (HPV types 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59) were more common in cervical cancers than among women without cervical lesions. A few rare HPV types also were more common in cervical cancers (eg, HPV26, 67, 68, 69, 73, 82). The follow-up studies found increased long-term risks particularly for HPV types 16, 18, 31, and 33, which had 14-year cumulative incidences for CIN3+above 28%, while HPV35, 45, 52, and 58 had 14-year risks between 14%-18% and HPV39, 51, 56, 59, 66, and 68 had risks<10%. HPV16 contributed to the greatest proportion of CIN2+(first-round population attributable proportion [PAR] 36%), followed by types 31, 52, 45, and 58 (7%-11%). HPV16, 18, 31, 33, 45, 52, and 58 together contributed 73.9% of CIN2+lesions and all high-risk types contributed 86.9%.In summary, the different oncogenic HPV types have substantial differences in their oncogenic potential. These differences are relevant for the design and evaluation of cervical screening tests and programs, as well as for studying the effect of vaccination programs using different HPV vaccines.

High genital prevalence of cutaneous human papillomavirus DNA on male genital skin: the HPV Infection in Men Study

Background

The genital skin of males hosts a diversity of HPV genotypes and uncharacterized HPV genotypes. Previously we demonstrated that a specific viral genotype was not identified in 14% of all genital specimens (i.e., HPV unclassified specimens) using the Roche Linear Array method. Our goal was to identify and assess the prevalence of individual HPV types among genital HPV unclassified specimens collected in the HIM Study population, at enrollment, and examine associations with socio-demographic and behavioral characteristics.

Methods

Genital skin specimens of men that were considered unclassified (HPV PCR positive, no genotype specified) at enrollment were typed by sequencing amplified PGMY09/11 products or cloning of PGMY/GP+ nested amplicons followed by sequencing. PGMY/GP+ negative specimens were further analyzed using FAP primers. HPV type classification was conducted through comparisons with sequences in the GenBank database.

Results

Readable nucleotide sequences were generated for the majority of previously unclassified specimens (66%), including both characterized (77%) and yet uncharacterized (23%) HPV types. Of the characterized HPV types, most (73%) were Beta [β]-HPVs, primarily from β-1 and β-2 species, followed by Alpha [α]-HPVs (20%). Smokers (current and former) were significantly more likely to have an α-HPV infection, compared with any other genus; no other factors were associated with specific HPV genera or specific β-HPV species.

Conclusions

Male genital skin harbor a large number of β-HPV types. Knowledge concerning the prevalence of the diverse HPV types in the men genital is important to better understand the transmission of these viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-014-0677-y) contains supplementary material, which is available to authorized users.

Deep sequencing extends the diversity of human papillomaviruses in human skin

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

Human papillomavirus community in healthy persons, defined by metagenomics analysis of human microbiome project shotgun sequencing data sets

Human papillomavirus (HPV) causes a number of neoplastic diseases in humans. Here, we show a complex normal HPV community in a cohort of 103 healthy human subjects, by metagenomics analysis of the shotgun sequencing data generated from the NIH Human Microbiome Project. The overall HPV prevalence was 68.9% and was highest in the skin (61.3%), followed by the vagina (41.5%), mouth (30%), and gut (17.3%). Of the 109 HPV types as well as additional unclassified types detected, most were undetectable by the widely used commercial kits targeting the vaginal/cervical HPV types. These HPVs likely represent true HPV infections rather than transitory exposure because of strong organ tropism and persistence of the same HPV types in repeat samples. Coexistence of multiple HPV types was found in 48.1% of the HPV-positive samples. Networking between HPV types, cooccurrence or exclusion, was detected in vaginal and skin samples. Large contigs assembled from short HPV reads were obtained from several samples, confirming their genuine HPV origin. This first large-scale survey of HPV using a shotgun sequencing approach yielded a comprehensive map of HPV infections among different body sites of healthy human subjects.

IMPORTANCE

This nonbiased survey indicates that the HPV community in healthy humans is much more complex than previously defined by widely used kits that are target selective for only a few high- and low-risk HPV types for cervical cancer. The importance of nononcogenic viruses in a mixed HPV infection could be for stimulating or inhibiting a coexisting oncogenic virus via viral interference or immune cross-reaction. Knowledge gained from this study will be helpful to guide the designing of epidemiological and clinical studies in the future to determine the impact of nononcogenic HPV types on the outcome of HPV infections.

Genetic variation of human papillomavirus type 16 in individual clinical specimens revealed by deep sequencing

Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled de novo assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis.

Identification and characterization of eleven novel human gamma-papillomavirus isolates from healthy skin, found at low frequency in a normal population

Eleven novel human papillomavirus (HPV) types were isolated and characterized from healthy individuals in China. HPV163 belongs to the γ-1 species, HPV 164 and HPV 168 fit in the γ-8 species, HPV 165 and KC5 belongs to the γ-12 species, HPV 168 is closely allied with the γ-4 species, HPV 169 is closely related to the γ-11 species, and HPV 170 is related to the γ-12 species. In addition, HPV 161, HPV 162, and HPV 166 may form a new HPV species of the γ-PV genus. The prevalence of these HPV types in the normal population is low.

Diversity of human papillomaviruses in skin lesions

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

Cutaneous alpha, beta and gamma human papillomaviruses in relation to squamous cell carcinoma of the skin: a population-based study

Human papillomavirus (HPV) infection is common worldwide and, in immunodeficient populations, may contribute to the pathogenesis of keratinocyte cancers, particularly squamous cell carcinomas (SCC). However, their role in SCC in the general population is less clear. We conducted a comprehensive analysis to investigate the independent effects of seropositivity for cutaneous alpha, beta and gamma HPV types on risk of SCC, and a meta-analysis of the available literature. In a population-based case-control study from New Hampshire, USA (n = 1,408), histologically confirmed SCC cases and controls were tested for L1 antibodies to alpha, beta and gamma cutaneous HPV types 2-5, 7-10, 15, 17, 20, 23, 24, 27b, 36, 38, 48-50, 57, 65, 75-77, 88, 92, 95, 96, 101, 103 and 107 using multiplex serology. An increasing risk of SCC with number of beta HPVs to which an individual tested positive was observed even among those seronegative for gamma types (p for trend = 0.016) with an odds ratio of 1.95 (95% confidence interval (CI) = 1.07-3.56) for four or more beta types positive. In a meta-analysis of six case-control studies, increased SCC risks in relation to beta HPV seropositivity were found across studies (meta odds ratio = 1.45, CI = 1.27-1.66). While the prevalence of gamma HPVs assayed was somewhat higher among SCC cases than controls, the association was only weakly evident among those seronegative for beta HPVs. Overall, the association between cutaneous HPVs and skin cancers appears to be specific to SCC and to genus beta HPVs in a general US population.