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

Genomic diversity of human papillomavirus type 6 from patients with condyloma acuminatum in Eastern China

OBJECTIVE

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts both men and women. However, there is limited data on its genomic characterization in mainland China. The aim of this study was to understand the complete genomic diversity of HPV6 from patients with condyloma acuminatum (CA) and to explore the prevalence of different variant lineages/sublineages in eastern China.

METHODS

CA samples were collected in 3 hospitals in Shandong Province, China from January 2020 to March 2021. DNA extraction, PCR amplification, Sanger sequencing and sequence assembly were performed on HPV6-positive samples. The complete genomes obtained in this study were analyzed phylogenetically with global HPV6 sequences in GenBank database using MEGA 11.

RESULTS

A total of 55 complete genomic sequences of HPV6 were obtained in this study. They were classified as HPV6 variant lineage A (n = 20), sublineage B1 (n = 34) and sublineage B3 (n = 1) by phylogenetic analysis. Sequence alignment showed E1, E5A, E5B, L1, L2, LCR were relatively highly variable regions for sublineage B1 whereas E1, E5A, L2 for lineage A. Both phylogenetic trees of lineage A and sublineage B1 composed of two main branches. Chinese sequences of lineage A segregated into the major branch while those in sublineage B1 belonged to both branches. Genomic divergence between sequences from China and other countries was 0.00% - 0.33% in lineage A and 0.00% - 0.40% in sublineage B1.

CONCLUSIONS

This is the first study on HPV variant lineages circulating in mainland China. The results revealed that lineage A and sublineage B1 were prevalent and they had different highly variable regions. Further surveillance is needed to understand the dynamic change of different variants in the population.

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.

Analysis of Nucleotide Alterations in the E6 Genomic Region of Human Papillomavirus Types 6 and 11 in Condyloma Acuminatum Samples from Brazil

Condyloma acuminata (CA), or genital warts, are benign proliferative epidermal or mucous lesions that are caused by infection with human papillomavirus (HPV), mainly the low-risk types 6 and 11. HPV variants are defined as viral sequences that share identity in the nucleotide sequence of the L1 gene greater than 98%. Based on this criterion, HPV6 and 11 variant lineages have been studied, and there are ongoing attempts to correlate these genetic variants with different clinical findings of infection. Therefore, the aims of this study were to detect variants and nucleotide alterations present in the E6 regions of HPV types 6 and 11 found in CA samples, to correlate the HPV presence with the clinical-pathological data of the patients, and to determine phylogenetic relationships with variants from other places in the world. The E6 regions of 25 HPV6 samples and 7 HPV11 samples from CA were amplified using PCR with specific primers. The products were ligated to a cloning vector and five colonies of each sample were sequenced to observe the nucleotide alterations. Twelve samples were identified as the HPV6B3 variant, presenting the mutation (guanine) G474A (adenine), and one of them also showed the mutation (thymine) T369G. The other 13 patients were positive for HPV6B1 without nucleotide alterations. In the analysis of the HPV11 samples, all patients showed the mutations T137C and (cytosine) C380T. One patient also presented the nucleotide alteration T410C. None of the mutations found in the 32 analyzed samples resulted in amino acid changes. Patient age, local occurrence, and HIV infection did not show significant association with HPV infection. Besides, the data found in this study did not show a relationship with the geographical region of isolation when compared to other data from different regions of the world. In this way, despite the nucleotide alterations found, it was not possible to observe amino acid changes and variants grouping according to geographical region.

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.

Methylation of CpG sites in the upstream regulatory region, physical status and mRNA expression of HPV-6 in adult-onset laryngeal papilloma

The methylation status of HPV-6 upstream regulatory region (URR) in adult-onset laryngeal papillomatosis (AO-LP) remains unclear. The purpose of this study was to investigate the methylation status of URR and the physical status of HPV-6, as well as the dynamic variations of viral load and mRNA expression in AO-LP. We examined 18 specimens from 11 patients with AO-LP by real-time polymerase chain reaction (PCR), bisulfite-sequencing PCR, and amplification of papilloma oncogene transcripts. HPV-6 was identified in 9 of 11 patients (81.8%), and all the 15 specimens derived from 9 HPV-6-positive cases contained only episomal HPV-6 transcripts with intact E2. Three HPV-6-positive patients developed recurrent lesions, and HPV-6 copy numbers and mRNA expression decreased after surgical treatment. Among the 96 CpG sites (16/case), 67 (69.8%) were unmethylated, while 23 (30.2%) were heterogeneous (≥ 1 methylated CpG clone). High viral loads and episomal status of HPV-6 were frequently observed in AO-LP; thus, persistent E6/E7 mRNA expression of LR-HPV-6 may be associated with AO-LP recurrences. Hypomethylation and scattered patterns of methylated CpGs at the URR of HPV-6 were identified in AO-LP.

HPV-6 Molecular Variants Association With the Development of Genital Warts in Men: The HIM Study

Background

Human papillomavirus type 6 (HPV-6) and HPV-11 are the etiological agents of approximately 90% of genital warts (GWs). The impact of HPV-6 genetic heterogeneity on persistence and progression to GWs remains undetermined.

Methods

HPV Infection in Men (HIM) Study participants who had HPV-6 genital swabs and/or GWs preceded by a viable normal genital swab were analyzed. Variants characterization was performed by polymerase chain reaction sequencing and samples classified within lineages (A, B) and sublineages (B1, B2, B3, B4, B5). Country- and age-specific analyses were conducted for individual variants; odds ratios and 95% confidence intervals for the risk of GWs according to HPV-6 variants were calculated.

Results

B3 variants were most prevalent. HPV-6 variants distribution differed between countries and case status. HPV-6 B1 variants prevalence was increased in GWs and genital swabs of cases compared to controls. There was difference in B1 and B3 variants detection in GW and the preceding genital swab. We observed significant association of HPV-6 B1 variants detection with GW development.

Conclusions

HPV-6 B1 variants are more prevalent in genital swabs that precede GW development, and confer an increased risk for GW. Further research is warranted to understand the possible involvement of B1 variants in the progression to clinically relevant lesions.

Determination of the complete genome and functional analysis of HPV6 isolate VBD19/10 from a patient with aggressive recurrent respiratory papillomatosis

Human papillomavirus (HPV) types 6 and 11 are the aetiological agent of recurrent respiratory papillomatosis (RRP). The complete genome of an HPV6 isolate with a 170 base pair (bp) duplication identified within the long control region (LCR) from a patient with aggressive recurrent respiratory papillomatosis was determined. The promoter sequence from the HPV LCR including the 170 bp duplication was placed upstream of a heterologous reporter gene and the activity of the reporter gene product determined using transfected cells. In total, mutations were observed at 157 nucleotide positions of the complete genome and included nucleotide substitutions, deletions and insertions, resulting in amino acid changes at 43 residue positions. Reporter gene activity using an HPV-derived LCR region with a 170 bp duplication was significantly higher than that using an HPV-derived LCR region with no duplication within this region. The results suggest that novel HPV variants warrant further investigation for potential biomarkers of aggressive disease.

Bovine papillomavirus type 7 in Italy: complete genomes and sequence variants

Two novel bovine papillomavirus type 7 (BPV-7) variants have been identified in teat cutaneous papillomas affecting dairy cows in northern Italy. The entire genome sequences of two BPV-7 Italian variants showed major sequence differences in the long control region (LCR) and in the L2 gene compared to the Japanese reference strain. In order to define the stability of these genetic variants, the L2 and LCR sequences of seven further BPV-7 positive isolates were characterized. An insertion of six amino acids in the L2 structural protein has been detected in all samples while different genetic variants have been identified for the LCR. These findings provide new insights on intra-type variability of BPVs and represent a starting point for future studies aimed at establishing the biological role of the different BPV genomic regions and investigating the pathogenic potential of papillomavirus variants.

Differences in Transcriptional Activity of Human Papillomavirus Type 6 Molecular Variants in Recurrent Respiratory Papillomatosis

A significant proportion of recurrent respiratory papillomatosis (RRP) is caused by human papillomavirus type 6 (HPV-6). The long control region (LCR) contains cis-elements for regulation of transcription. Our aim was to characterize LCR HPV-6 variants in RRP cases, compare promoter activity of these isolates and search for cellular transcription factors (TFs) that could explain the differences observed. The complete LCR from 13 RRP was analyzed. Transcriptional activity of 5 variants was compared using luciferase assays. Differences in putative TFs binding sites among variants were revealed using the TRANSFAC database. Chromatin immunoprecipation (CHIP) and luciferase assays were used to evaluate TF binding and impact upon transcription, respectively. Juvenile-onset RRP cases harbored exclusively HPV-6vc related variants, whereas among adult-onset cases HPV-6a variants were more prevalent. The HPV-6vc reference was more transcriptionally active than the HPV-6a reference. Active FOXA1, ELF1 and GATA1 binding sites overlap variable nucleotide positions among isolates and influenced LCR activity. Furthermore, our results support a crucial role for ELF1 on transcriptional downregulation. We identified TFs implicated in the regulation of HPV-6 early gene expression. Many of these factors are mutated in cancer or are putative cancer biomarkers, and must be further studied.

Global genomic diversity of human papillomavirus 6 based on 724 isolates and 190 complete genome sequences

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution.

IMPORTANCE

This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

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.

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.