Prevaccination genomic diversity of human papillomavirus genotype 11: A study on 63 clinical isolates and 10 full-length genome sequences

Study on the detection rate, genetic polymorphism, viral load, persistent infection capacity, and pathogenicity of human papillomavirus type 81

Human papillomavirus (HPV) type 81 has recently become one of the most common low-risk HPV types; however, literature focusing on it is limited. This study aimed to analyze the reasons for the increased detection rate of HPV81 and investigate its evolving pathogenicity. We analyzed the detection rates and trends of HPV81 in 229 061 exfoliated cervical cell samples collected from 2014 to 2023; collected samples of HPV81 single infections from two different time periods; and analyzed the allele frequencies, positive selection, viral load, persistent infection capacity, and pathogenicity of E6 and E7 genotypes. We found that the detection rate of HPV81 ranked first among the low-risk types in exfoliated cervical cells and exhibited a significantly increasing trend (p < 0.001). The frequency of the E6 prototype allele of HPV81 (n = 317) was significantly increased (p = 0.018) and demonstrated the strongest adaptive capacity. The viral load and persistent infection capacity of the E6 prototype were significantly higher than those of the mutants, thus serving as key drivers for increasing the detection rate of HPV81 and enhancing its pathogenicity. The viral load was positively correlated with persistent infection capacity and pathogenicity. Persistent infection was a crucial factor in the pathogenicity of HPV81. Successful adaptive evolution of HPV81 is accompanied by enhanced pathogenicity.

Prevalence of human papillomavirus 6 and 11 variants in recurrent respiratory papillomatosis

Human papillomavirus (HPV) types 6 and 11 are the etiological agents of recurrent respiratory papillomatosis (RRP). We examined the prevalence and distribution of HPVs 6 and 11 genetic variants in juvenile onset (JORRP) and adult onset (AORRP) laryngeal papillomas. Cases of JORRP and AORRP were collected, retrospectively. HPV detection and genotyping were accessed by polymerase chain reaction-sequencing in 67 RRP samples. Overall, the most prevalent HPV-6 variants were from B1 (55.8%) and B3 (27.9%) sublineages, whereas among HPV-11 positive samples A2 (62.5%) variants were predominant. A higher prevalence of HPV-6 B1 was observed in JORRP (83.3% B1 and 16.7% B3), compared with AORRP cases (58.3% B1 and 41.7% B3). HPV-11 A2 variants were more prevalent both in JORRP (57.2%) and in AORRP cases (70.0%). Nevertheless, with the exception that HPV-6 B1 were significantly less likely to recur, there was a lack of association between any particular HPVs 6 or 11 variant and clinicopathological features. Our data do not support an association between HPVs 6 and 11 variability and RRP.

Comparative analysis of human papillomavirus type 6 complete genomes originated from head and neck and anogenital disorders

It is increasingly recognized that fundamental differences exist between high-risk and low-risk human papillomavirus (HPV) genotypes regarding interactions with the host. This study aims to join the recently emerging efforts to uncover these differences at the complete genome level and to study how they may influence the disease caused. Sixteen samples of thirteen patients with various HPV6-mediated benign mucosal disorders (nine recurrent respiratory papillomatoses with 2-8 recurrences, one condyloma acuminatum and three premalignant lesions of the genital mucosa) were sampled to determine the complete virus genomes. We collected the 197 HPV6 complete genomes deposited in the GenBank for cluster analysis to determine (sub)lineages. Genome polymorphisms were determined against the reference sequences of the (sub)lineages. Genome polymorphisms of the long control region (LCR) were tested for putative transcription factor binding sites; their functional analysis was performed by transient transfection of cloned whole LCRs into HEp-2 cells using a luciferase reporter system. Genomes from the same patients were always identical. Three, nine and one patients carried HPV6 lineage A, sublineage B1 and B2 variants, respectively. The three lineage A sequences were highly similar to each other, but distinct from the reference genome. A unique non-synonymous single nucleotide polymorphism (SNP) was found in the E5a open reading frame (ORF). Sublineage B1 genomes were more diverse, exhibited unique non-synonymous SNPs in the LCR and the E2/E4, L1, L2 ORFs. LCR activity of lineage A and sublineage B1 differed significantly; activity of one sublineage B1 LCR exhibiting two unique SNPs was significantly higher than that of other B1 LCR variants, close to the mean of LCR activities of lineage A variants. Different HPV6 lineages showed marked differences in variability patterns of the different genome regions. This may be involved in the differences in their distribution in different diseases or patient populations.

HPV-11 variability, persistence and progression to genital warts in men: the HIM study

HPV-11 and HPV-6 are the etiological agents of about 90 % of genital warts (GWs). The intra-typic variability of HPV-11 and its association with infection persistence and GW development remains undetermined. Here, HPV infection in men (HIM) participants who had an HPV-11 genital swab and/or GW, preceded or not by a normal skin genital swab were analysed. Genomic variants were characterized by PCR-sequencing and classified within lineages (A, B) and sublineages (A1, A2, A3, A4). HPV-11 A2 variants were the most frequently detected in the genital swab samples from controls and in both genital swabs and GW samples from cases. The same HPV-11 variant was detected in the GW sample and its preceding genital swab. There was a lack of association between any particular HPV-11 variant and the increased risk for GW development.

Global Genomic Diversity of Human Papillomavirus 11 Based on 433 Isolates and 78 Complete Genome Sequences

Human papillomavirus 11 (HPV11) is an etiological agent of anogenital warts and laryngeal papillomas and is included in the 4-valent and 9-valent prophylactic HPV vaccines. We established the largest collection of globally circulating HPV11 isolates to date and examined the genomic diversity of 433 isolates and 78 complete genomes (CGs) from six continents. The genomic variation within the 2,800-bp E5a-E5b-L1-upstream regulatory region was initially studied in 181/207 (87.4%) HPV11 isolates collected for this study. Of these, the CGs of 30 HPV11 variants containing unique single nucleotide polymorphisms (SNPs), indels (insertions or deletions), or amino acid changes were fully sequenced. A maximum likelihood tree based on the global alignment of 78 HPV11 CGs (30 CGs from our study and 48 CGs from GenBank) revealed two HPV11 lineages (lineages A and B) and four sublineages (sublineages A1, A2, A3, and A4). HPV11 (sub)lineage-specific SNPs within the CG were identified, as well as the 208-bp representative region for CG-based phylogenetic clustering within the partial E2 open reading frame and noncoding region 2. Globally, sublineage A2 was the most prevalent, followed by sublineages A1, A3, and A4 and lineage B.

IMPORTANCE

This collaborative international study defined the global heterogeneity of HPV11 and established the largest collection of globally circulating HPV11 genomic variants to date. Thirty novel complete HPV11 genomes were determined and submitted to the available sequence repositories. Global phylogenetic analysis revealed two HPV11 variant lineages and four sublineages. The HPV11 (sub)lineage-specific SNPs and the representative region identified within the partial genomic region E2/noncoding region 2 (NCR2) will enable the simpler identification and comparison of HPV11 variants worldwide. This study provides an important knowledge base for HPV11 for future studies in HPV epidemiology, evolution, pathogenicity, prevention, and molecular assay development.

Genomic diversity of human papillomaviruses (HPV) and clinical implications: an overview in adulthood and childhood

During the last years, several researchers have highlighted the importance of characterizing more than one genomic region in order to detect recombination and classify variants of human papillomaviruses (HPVs) properly. HPVs variants differ in their biological, molecular and chemical properties. Therefore, this genomic diversity can present differences in the natural history and pathogenicity of HPVs. Different 'high-risk' HPVs variants of the genotypes HPV 16 and 18 can confer varied risks of viral persistence in the human cervix and influence HPVs progression to cervical cancer. Moreover, different 'low-risk' HPVs variants of the genotypes HPV 6 and 11 can play a unique role in the development of anogenital and cutaneous warts, recurrent respiratory papillomatosis (RRP) and ophthalmic pterygium. In future, the precise impact of genomic HPVs diversity to the clinical course of HPVs-associated diseases as well as to the efficacy of the current HPVs vaccines remains to be elucidated.

Genomic diversity of low-risk human papillomavirus genotypes HPV 40, HPV 42, HPV 43, and HPV 44

In order to investigate the genomic diversity of low-risk human papillomavirus (HPV) genotypes, a total of 108 isolates of HPV 40, HPV 42, HPV 43, or HPV 44, obtained from anal swabs or tissue specimens of patients with anogenital warts, and cervical swabs of women with cervical intraepithelial neoplasia of different grades, were analyzed. The characterization of genomic variants was established by sequencing one third of the viral genome and analysis of three different genomic regions: L1, LCR, and E6. Maximum variant divergence accounted for 0.4-1.1% of the investigated genomic segments. Several novel, potentially important nucleotide substitutions, deletions, and insertions are described. Altogether, among 14 HPV 40 isolates, a total of nine different genomic variants were identified, composed of eight L1, five LCR, and four E6 genomic variants. Among 49 HPV 42 isolates, a total of 30 genomic variants were identified, composed of 20 L1, 18 LCR, and four E6 genomic variants. Among 10 HPV 43 isolates, distributed into two major genomic variant lineages with clearly defined nucleotide signatures, three genomic variants were identified, composed of three L1, two LCR, and two E6 genomic variants. Among 35 HPV 44 isolates, a total of eight HPV 44 and 11 subtype HPV 44 genomic variants were identified, composed of 13 L1, 14 LCR, and 6 E6 genomic variants. A similar level of genomic diversity of HPV 44 and its subtype was identified in our geographic region as has been reported previously on isolates collected worldwide.

Human papillomavirus type 6 and 11 genetic variants found in 71 oral and anogenital epithelial samples from Australia

Genetic variation of 49 human papillomavirus (HPV) 6 and 22 HPV11 isolates from recurrent respiratory papillomatosis (RRP) (n = 17), genital warts (n = 43), anal cancer (n = 6) and cervical neoplasia cells (n = 5), was determined by sequencing the long control region (LCR) and the E6 and E7 genes. Comparative analysis of genetic variability was examined to determine whether different disease states resulting from HPV6 or HPV11 infection cluster into distinct variant groups. Sequence variation analysis of HPV6 revealed that isolates cluster into variants within previously described HPV6 lineages, with the majority (65%) clustering to HPV6 sublineage B1 across the three genomic regions examined. Overall 72 HPV6 and 25 HPV11 single nucleotide variations, insertions and deletions were observed within samples examined. In addition, missense alterations were observed in the E6/E7 genes for 6 HPV6 and 5 HPV11 variants. No nucleotide variations were identified in any isolates at the four E2 binding sites for HPV6 or HPV11, nor were any isolates found to be identical to the HPV6 lineage A or HPV11 sublineage A1 reference genomes. Overall, a high degree of sequence conservation was observed between isolates across each of the regions investigated for both HPV6 and HPV11. Genetic variants identified a slight association with HPV6 and anogenital lesions (p = 0.04). This study provides important information on the genetic diversity of circulating HPV 6 and HPV11 variants within the Australian population and supports the observation that the majority of HPV6 isolates cluster to the HPV6 sublineage B1 with anogenital lesions demonstrating an association with this sublineage (p = 0.02). Comparative analysis of Australian isolates for both HPV6 and HPV11 to those from other geographical regions based on the LCR revealed a high degree of sequence similarity throughout the world, confirming previous observations that there are no geographically specific variants for these HPV types.

HPV associated with recurrent respiratory papillomatosis

Papillomaviruses are members of the Papillomaviridae family. Over 150 HPV types have been identified. Recurrent respiratory papillomatosis (RRP) is a chronic condition caused by HPV characterized by recurrent papillomas of the respiratory tract, mainly the larynx. During the early stages, the condition presents with hoarseness, while more advanced disease presents with stridor and respiratory distress. There is no specific cure and treatment consists of repeated surgical procedures to remove the papillomas. Most patients eventually go into remission, but some suffer for many years with this condition, which may be fatal. HPV-6 and HPV-11 are the HPV types most commonly associated with RRP. Although most studies have found RRP due to HPV-11 to be more aggressive than disease due to HPV-6, the variability in disease aggressiveness is probably multifactorial. Information regarding the current epidemiology, molecular diversity and host immune responses is important for strategizing ways to reduce disease. Data on HPV genotypes associated with RRP would provide valuable information for vaccination programs to reduce the incidence of these genotypes in mothers and, in the long term, reduce the incidence of RRP in children. This review focuses on HPV-6 and HPV-11 as the HPV types that cause RRP, and discusses the viral genome and replication, clinical presentation of RRP, current techniques of diagnosis and genotyping, and the molecular diversity of HPV-6 and HPV-11.

FRET-based detection and genotyping of HPV-6 and HPV-11 causing recurrent respiratory papillomatosis

Recurrent respiratory papillomatosis (RRP) is a potentially life-threatening disease caused by human papillomavirus (HPV), usually HPV types 6 and 11. The conventional method used for detection and typing the RRP isolates in our laboratory is the polymerase chain reaction (PCR) and DNA sequencing method. A real-time PCR assay based on fluorescence resonance energy transfer (FRET) probe technology was developed for the detection and rapid genotyping of HPV-6 and-11 isolates from biopsy material. The primers and probes were designed using multiple alignments of HPV-6 and HPV-11 partial E6 and E7 sequences that included prototypic and non-prototypic variants. Real-time PCR followed by probe-specific melting-curve analysis allowed differentiation of HPV-6 and HPV-11. HPV-6 and HPV-11 amplicons were used to determine detection limits and inter- and intra-assay variability. The detection limit of the assay was 12.8 DNA copies for HPV-6 and 22.5 DNA copies for HPV-11. A total of 60 isolates were genotyped using the FRET real-time PCR assay and a 100% concordance was obtained when results were compared with genotyping based on conventional DNA sequencing. The real-time PCR assay based on FRET technology was able to detect and rapidly genotype HPV from tissue biopsy obtained from patients with RRP. The assay reduces the time required for genotyping from three working days to less than a day.

Nucleotide and phylogenetic analysis of human papillomavirus types 6 and 11 isolated from recurrent respiratory papillomatosis in Brazil

There are few studies about the distribution of natural molecular variants of low-risk HPVs. Our aim was to evaluate the E6 early gene variability among HPV-6 and HPV-11 isolates detected in recurrent respiratory papillomatosis (RRP) samples obtained in a cohort of Brazilian patients. We also performed a phylogenetic analysis in order to compare nucleotide sequences identified in our study with previously reported isolates from different anatomic sites (laryngeal papillomas, genital warts, cervical cancer and anal swabs) obtained from other parts of the world to determine the phylogenetic relationships of variants detected in Brazil. The complete coding region of the E6 gene of 25 samples was cloned and sequenced: 18 isolates of HPV-6 (72%) and 7 isolates of HPV-11 (28%). A total of four different HPV-6 genomic variants and two HPV-11 genomic variants was identified. It was not possible to correlate specific variants with disease severity. Phylogenetic trees for both HPV types were constructed enclosing both E6 sequences detected in our study and formerly published sequences. In both phylogenetic trees, the sequences from Brazil did not group together. We could not establish a geographical association between HPV-6 or HPV-11 variants, unlike HPV-16 and HPV-18.

Identical human papillomavirus (HPV) genomic variants persist in recurrent respiratory papillomatosis for up to 22 years

Seventy initial and 125 follow-up tissue specimens of laryngeal papillomas, obtained from 70 patients who had had recurrent respiratory papillomatosis for from 1-22 years, were investigated for the presence of human papillomavirus (HPV) DNA and HPV E5a, LCR and/or full-length genomic variants. HPV-6 was found in 130/195, HPV-11 in 63/195, and HPV-6/HPV-11 in 2/195 samples. Within 67/70 (95.7%) patients, all follow-up HPV isolates genetically matched completely initial HPV isolate over the highly variable parts of the genome or over the entire genome. Frequent recurrence of laryngeal papillomas is a consequence of long-term persistence of the identical initial HPV genomic variant.

Novel HPV-6 variants of human papillomavirus causing recurrent respiratory papillomatosis in southern Africa

There is currently no information regarding the genetic diversity of HPV-6 variants circulating in South Africa. The aim of this study was to determine the HPV-6 variants affecting patients with recurrent respiratory papillomatosis, to determine whether mutations correlate with disease severity and identify molecular determinants of virulence with prognostic relevance. HPV-6 variants were identified based on genome changes within the 712-991 bp region encompassing the non-coding region (URR) of the genome, with variations in length resulting from insertions and duplications, and the 453-bp gene encoding the E6 protein. Based on manual comparison of sequence data from the URR, the isolates were identified as HPV-6a and HPV-6vc variants. Three novel HPV-6 variants were identified: one based on a mutation in the E6 region; two based on changes in the URR including a unique substitution detected in three isolates and an insertion and 170-bp duplication in the URR genome in one patient, who had clinical features of severe disease.

Evolution and taxonomic classification of human papillomavirus 16 (HPV16)-related variant genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67

Background

Human papillomavirus 16 (HPV16) species group (alpha-9) of the Alphapapillomavirus genus contains HPV16, HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. These HPVs account for 75% of invasive cervical cancers worldwide. Viral variants of these HPVs differ in evolutionary history and pathogenicity. Moreover, a comprehensive nomenclature system for HPV variants is lacking, limiting comparisons between studies.

Methods

DNA from cervical samples previously characterized for HPV type were obtained from multiple geographic regions to screen for novel variants. The complete 8 kb genomes of 120 variants representing the major and minor lineages of the HPV16-related alpha-9 HPV types were sequenced to capture maximum viral heterogeneity. Viral evolution was characterized by constructing phylogenic trees based on complete genomes using multiple algorithms. Maximal and viral region specific divergence was calculated by global and pairwise alignments. Variant lineages were classified and named using an alphanumeric system; the prototype genome was assigned to the A lineage for all types.

Results

The range of genome-genome sequence heterogeneity varied from 0.6% for HPV35 to 2.2% for HPV52 and included 1.4% for HPV31, 1.1% for HPV33, 1.7% for HPV58 and 1.1% for HPV67. Nucleotide differences of approximately 1.0% - 10.0% and 0.5%–1.0% of the complete genomes were used to define variant lineages and sublineages, respectively. Each gene/region differs in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) /noncoding region 2 (NCR2) >upstream regulatory region (URR)> E6/E7 > E2/L2 > E1/L1.

Conclusions

These data define maximum viral genomic heterogeneity of HPV16-related alpha-9 HPV variants. The proposed nomenclature system facilitates the comparison of variants across epidemiological studies. Sequence diversity and phylogenies of this clinically important group of HPVs provides the basis for further studies of discrete viral evolution, epidemiology, pathogenesis and preventative/therapeutic interventions.