Trend of HPV Molecular Epidemiology in the Post-Vaccine Era: A 10-Year Study.

Cervical cancer, a major health concern among women worldwide, is closely linked to human papillomavirus (HPV) infection. This study explores the evolving landscape of HPV molecular epidemiology in Taiwan over a decade (2010-2020), where prophylactic HPV vaccination has been implemented since 2007. Analyzing data from 40,561 vaginal swab samples, with 42.0% testing positive for HPV, we reveal shifting trends in HPV genotype distribution and infection patterns. The 12 high-risk genotypes, in order of decreasing percentage, were HPV 52, 58, 16, 18, 51, 56, 39, 59, 33, 31, 45, and 35. The predominant genotypes were HPV 52, 58, and 16, accounting for over 70% of cases annually. The proportions of high-risk and non-high-risk HPV infections varied across age groups. High-risk infections predominated in sexually active individuals aged 30-50 and were mixed-type infections. The composition of high-risk HPV genotypes was generally stable over time; however, HPV31, 33, 39, and 51 significantly decreased over the decade. Of the strains, HPV31 and 33 are shielded by the nonavalent HPV vaccine. However, no reduction was noted for the other seven genotypes. This study offers valuable insights into the post-vaccine HPV epidemiology. Future investigations should delve into HPV vaccines' effects and their implications for cervical cancer prevention strategies. These findings underscore the need for continued surveillance and research to guide effective public health interventions targeting HPV-associated diseases.

Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer.

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Head-to-head comparison of genotyping of human papillomavirus by real-time multiplex PCR assay using type-specific primers and SPF10-PCR-based line probe assay.

The SPF10-polymerase chain reaction (PCR)-based line probe assay (LiPA-25) with high analytical sensitivity and specificity for human papillomavirus (HPV) genotyping in clinical samples has been widely used in vaccine and epidemiologic studies. A real-time multiplex PCR assay using type-specific primers (Hybribio-23) with low workload and cost has been developed recently. The study aimed to compare the performance of LiPA-25 and Hybribio-23 in selected 1731 cervical swab and 117 tissue samples, with a focus on 20 common HPV types (14 high-risk: 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68/73; 6 low-risk: 6, 11, 42, 43, 44, and 53). The level of agreement of two assays was determined using Cohen's Kappa (κ) statistics. A total of 1296 (74.9%) swab samples were identified as HPV-positive by Hybribio-23 or LiPA-25, of which 814 (62.8%) samples exhibited concordant, 358 (27.6%) showed additional or fewer types (compatible), and 124 (9.6%) were discordant. In addition, the two assays showed a perfect agreement for 20 HPV-combined detection (κ = 0.838) and 17 individual HPV types (all κ > 0.800), a good agreement for HPV31 (κ = 0.792) and 43 (κ = 0.696), and a moderate agreement for HPV42 (κ = 0.504). Hybribio-23 was significantly more sensitive for HPV58, 59, 68/73, 42, 43, and 44, and less sensitive for HPV35 and 66 than LiPA-25 (McNemar's test: all p < 0.05). For 117 HPV-positive tissue specimens, the identification of genotypes was 85.2% identical, 12.2% compatible, and only 2.6% discordant. The agreement for HPV31 (κ = 0.786), 68/73 (κ = 0.742), and HPV53 (κ = 0.742) was good, while for other types (all κ > 0.853) and 20 HPV-combined detection (κ = 0.936) was perfect (all p > 0.05). In conclusion, Hybribio-23 and LiPA-25 are comparable. Hybribio-23 could be used for the detection and genotyping of HPV in cervical samples for epidemiological and vaccine studies worldwide.

SETD6 Regulates E2-Dependent Human Papillomavirus Transcription.

Bromodomain-containing protein 4 (Brd4) is a member of the bromodomain and extraterminal domain (BET) family of proteins. Brd4 regulates human papillomavirus (HPV) transcription, genome replication, and segregation by binding to the E2 protein. The SETD6 methyltransferase binds to and methylates Brd4 at lysine 99. We investigated the interactions of SETD6 and Brd4 with E2 and their role in HPV transcription. SETD6 coimmunoprecipitated with the E2 transactivation domain, and its depletion in CIN612 episomal cells reduced human papillomavirus type 31 (HPV-31) transcription, whereas depletion of SETD6 in integrated HPV cell lines had no effect on viral gene expression. The mutant Brd4 K99R (bearing a change of K to R at position 99), which cannot be methylated by SETD6, displayed decreased binding to HPV-31 E2, suggesting that SETD6 methylation of Brd4 also influences E2 association with the Brd4 protein. Using chromatin immunoprecipitation, SETD6 was detected at the enhancer region of the HPV long control region. We propose that methylation of Brd4 at K99 by SETD6 is an important mechanism for E2-Brd4 association and HPV transcriptional activation. IMPORTANCE Human papillomaviruses (HPV) cause cervical, anogenital, and oral cancers. Brd4 plays an important role in the HPV life cycle. SETD6 was recently shown to methylate Brd4. The current study demonstrates that methylation of Brd4 by SETD6 in HPV-episomal cells is required for the activation of viral transcription. This study illustrates a novel regulatory mechanism involving E2, Brd4, and SETD6 in the HPV life cycle and provides insight into the multiple roles of Brd4 in viral pathogenesis.

Analysis of the genomic diversity of human papillomavirus type 31 in cervical samples reveals the presence of novel sublineages in clade C.

Human papillomavirus 31 (HPV31) is the fourth most frequent high-risk HPV (HR-HPV) genotype identified in cervical cancer (CC) worldwide and in Mexico. It has been recently classified into three lineages (A, B, and C) and eight sublineages (A1, A2, B1, B2, and C1 - C4). Here, we report the complete genomic sequences of 14 HPV31 isolates from cervical samples, and these were compared with viral genome sequences from the GenBank database for phylogenetic and genetic distance analysis. The formation of two novel clades within the C lineage (proposed as C5 and C6) was observed, with a well-defined variant-specific mutational pattern. The smallest average pairwise distance was 0.71% for lineages A and B, 0.94% for lineages A and C, and 1.01% for lineages B and C, and between sublineages, these values were 0.21% for clade A, 0.29% for clade B, and 0.24% for clade C. The isolates were grouped into the sublineages A1, B2, C1-C3, and C6. This is the first report on the whole-genome diversity of HPV31 in Mexico.

Phylogenomic Analysis of Human Papillomavirus Type 31 and Cervical Carcinogenesis: A Study of 2093 Viral Genomes.

Human papillomavirus (HPV) type 31 (HPV31) is closely related to the most carcinogenic type, HPV16, but only accounts for 4% of cervical cancer cases worldwide. Viral genetic and epigenetic variations have been associated with carcinogenesis for other high-risk HPV types, but little is known about HPV31. We sequenced 2093 HPV31 viral whole genomes from two large studies, one from the U.S. and one international. In addition, we investigated CpG methylation in a subset of 175 samples. We evaluated the association of HPV31 lineages/sublineages, single nucleotide polymorphisms (SNPs) and viral methylation with cervical carcinogenesis. HPV31 A/B clade was >1.8-fold more associated with cervical intraepithelial neoplasia grade 3 and cancer (CIN3+) compared to the most common C lineage. Lineage/sublineage distribution varied by race/ethnicity and geographic region. A viral genome-wide association analysis identified SNPs within the A/B clade associated with CIN3+, including H23Y (C626T) (odds ratio = 1.60, confidence intervals = 1.17-2.19) located in the pRb CR2 binding-site within the E7 oncogene. Viral CpG methylation was higher in lineage B, compared to the other lineages, and was most elevated in CIN3+. In conclusion, these data support the increased oncogenicity of the A/B lineages and suggest variation of E7 as a contributing risk factor.

Complete genome sequence of a HPV31 isolate from laryngeal squamous cell carcinoma and biological consequences for p97 promoter activity.

Human papillomavirus type 31, although detected less frequently than HPV types 16 and 18, is associated with head and neck squamous cell carcinomas. Previous studies suggest that polymorphisms in the long control region (LCR) may alter the oncogenic potential of the virus. This study reports the first complete genome of a South African HPV31 isolate from a laryngeal squamous cell carcinoma. Sequence variations relative to the HPV31 prototype sequence were identified. The pBlue-Topo® vector, a reporter gene system was used to investigate the possible influence of these variations on the LCR promoter activity in vitro. Using mutagenesis to create two different fragments, ß-galactosidase assays were used to monitor the effect of nucleotide variations on the p97 promoter. Increased ß-galactosidase expression was observed in mutants when compared to the South African HPV31 LCR isolate. Enhanced transcriptional activity was observed with the mutant that possessed a single nucleotide change within the YY1 transcription factor binding site. In conclusion, sequence variation within the LCR of HPV31 isolates may have a functional effect on viral p97 promoter activity.

Regulation of Human Papillomavirus 18 Genome Replication, Establishment, and Persistence by Sequences in the Viral Upstream Regulatory Region.

During persistent human papillomavirus infection, the viral genome replicates as an extrachromosomal plasmid that is efficiently partitioned to daughter cells during cell division. We have previously shown that an element which overlaps the human papillomavirus 18 (HPV18) transcriptional enhancer promotes stable DNA replication of replicons containing the viral replication origin. Here, we perform comprehensive analyses to elucidate the function of this maintenance element. We conclude that no unique element or binding site in this region is absolutely required for persistent replication and partitioning and instead propose that the overall chromatin architecture of this region is important to promote efficient use of the replication origin. These results have important implications for the genome partitioning mechanism of papillomaviruses. IMPORTANCE Persistent infection with oncogenic human papillomaviruses (HPVs) is responsible for ∼5% of human cancers. The viral DNA replicates as an extrachromosomal plasmid and is partitioned to daughter cells in dividing keratinocytes. Using a complementation assay that allows us to separate viral transcription and replication, we provide insight into viral sequences that are required for long-term replication and persistence in keratinocytes. Understanding how viral genomes replicate persistently for such long periods of time will guide the development of antiviral therapies.

Distribution of human papillomavirus genotypes in western China and their association with cervical cancer and precancerous lesions.

The aim of this study was to describe the distribution of human papillomavirus (HPV) genotypes among cervical cancers and pre-cancers in Shaanxi province of western China. A total of 17,341 women who were screened for cervical cancer from January 2014 to December 2016, using HPV genotyping and ThinPrep cytologic test were included. The prevalence and attribution of HPV genotypes were stratified by cervical lesion and age group. Of the subjects, 26.3% were infected with HPV, 28.0% of whom had multiple infections. The crude HPV prevalence increased from atypical squamous cells of undetermined significance/low-grade squamous intraepithelial lesions (ASCUS/LSIL, 64.3%) to high-grade squamous intraepithelial lesions (HSIL, 79.8%) and to invasive cervical cancer (ICC, 89.7%, P < 0.001). The three most prevalent genotypes were HPV 16 (8.0%), 58 (4.2%), and 52 (4.0%), and HPV 16, 31 and 33 were positively correlated with increased severity of cervical lesions. Additionally, the divalent vaccine genotypes HPV 16 and 18 accounted for 68.2% of ICC cases. Although 78.5% of ICC and 60.3% of HSIL cases were attributed to 9-valent vaccine genotypes, the other genotypes not covered by any vaccine still resulted in increases in coverage, with 1.5% for ICC, 5.3% for HSIL, and 13.5% for ASCUS/LSIL. HPV prevalence in western China was consistent with other regions of China. Early vaccination with 9-valent HPV vaccine is recommended in this locality for females younger than 26 years with no prior infection, while divalent the vaccine is more appropriate for women between 26 and 45 years, considering the efficacy, safety and cost-effectiveness of vaccines.

Lineage analysis of human papillomavirus types 31 and 45 in cervical samples of Iranian women.

Knowing the regional lineages/sublineages of human papillomavirus 31 (HPV 31) and 45 would be of great importance for further evolutionary, epidemiological, and biological analysis. In this regard, to characterize more common lineages and sublineages of HPV 31 and 45, the sequence variations of E6 gene were investigated in normal, premalignant, and malignant samples collected from the cervix in Iran. In total, 54 HPV 31- and 24 HPV 45-positive samples were analyzed by hemi-nested polymerase chain reaction (PCR) and nested-PCR, respectively. All PCR products were subjected to direct sequencing analysis. The results indicated that all three lineages A, B, and C were detected in HPV 31-positive samples; among which HPV 31 lineage A was dominant as it was found in 66.7% of all samples. HPV 31 lineages B and C were identified in 5.5% and 27.8% of samples, respectively. In HPV 45-infected samples, lineage B comprised of 62.5% of all samples and the remaining 37.5%  belonged to lineage A. In conclusion, our findings showed that lineage A of HPV 31 was predominant in Iran. Lineage B of HPV 45 was also dominant among Iranian women. However, further studies with larger sample size should be addressed to estimate the pathogenicity risk of HPV 31 or HPV 45 lineages/sublineages in the development of cervical cancer among Iranian women.

Structural and functional impacts of E5 genetic variants of human papillomavirus type 31.

Persistent infections caused by high-risk human papillomavirus (HR-HPV) are important, for the development of cervical lesions, but environmental and genetic factors are also related in the process of carcinogenesis. Among the genetic factors, the genetic variants of HR-HPV appear to be related to the risk of persistent infections. Therefore, the present study investigates variants of HPV31 E5 oncogene in cervical scraping samples from Brazilian women to assess their functional and structural effects, in order to identify possible repercussions of these variants on the infectious and carcinogenic process. Our results detected nucleotide changes previously described in the HPV31 E5 oncogene, which may play a critical role in the development of cancer due to its ability to promote cell proliferation and signal transmission. In our study, the interaction percentage of the 31E5 sequence generated by the Immune Epitope Server database and the Analysis Resource (IEDB) allowed us to include possible immunogenic epitopes with the MHC-I and MHC-II molecules, which may represent a possible relationship between protein suppression of the immune system. In the structural analysis of the HPV31 E5 oncoprotein, the N5D, I48 V, P56A, F80I and V64I polymorphisms can be found inserted within transmembrane regions. The P56A mutation has been predicted to be highly stabilizing and, therefore, can cause a change in protein function. Regarding the interaction of the E5 protein from HPV31 with the signaling of NF-kB pathway, we observed that in all variants of the E5 gene from HPV-31, the activity of the NF-kB pathway was increased compared to the prototype. Our study contributes to a more refined design of studies with the E5 gene from HPV31 and provides important data for a better understanding of how variants can be distinguished under their clinical consequences.

Comprehensive Assessment of the Antigenic Impact of Human Papillomavirus Lineage Variation on Recognition by Neutralizing Monoclonal Antibodies Raised against Lineage A Major Capsid Proteins of Vaccine-Related Genotypes.

Human papillomavirus (HPV) is the causative agent of cervical and other epithelial cancers. Naturally occurring variants of HPV have been classified into lineages and sublineages based on their whole-genome sequences, but little is known about the impact of this diversity on the structure and function of viral gene products. The HPV capsid is an icosahedral lattice comprising 72 pentamers of the major capsid protein (L1) and the associated minor capsid protein (L2). We investigated the potential impact of this genome variation on the capsid antigenicity of lineage and sublineage variants of seven vaccine-relevant, oncogenic HPV genotypes by using a large panel of monoclonal antibodies (MAbs) raised against the L1 proteins of lineage A antigens. Each genotype had at least one variant that displayed a ≥4-fold reduced neutralizing antibody sensitivity against at least one MAb, demonstrating that naturally occurring variation can affect one or more functional antigenic determinants on the HPV capsid. For HPV16, HPV18, HPV31, and HPV45, the overall impact was of a low magnitude. For HPV33 (sublineages A2 and A3 and lineages B and C), HPV52 (lineage D), and HPV58 (lineage C), however, variant residues in the indicated lineages and sublineages reduced their sensitivity to neutralization by all MAbs by up to 1,000-fold, suggesting the presence of key antigenic determinants on the surface of these capsids. These determinants were resolved further by site-directed mutagenesis. These data improve our understanding of the impact of naturally occurring variation on the antigenicity of the HPV capsid of vaccine-relevant oncogenic HPV genotypes.IMPORTANCE Human papillomavirus (HPV) is the causative agent of cervical and some other epithelial cancers. HPV vaccines generate functional (neutralizing) antibodies that target the virus particles (or capsids) of the most common HPV cancer-causing genotypes. Each genotype comprises variant forms that have arisen over millennia and which include changes within the capsid proteins. In this study, we explored the potential for these naturally occurring variant capsids to impact recognition by neutralizing monoclonal antibodies. All genotypes included at least one variant form that exhibited reduced recognition by at least one antibody, with some genotypes affected more than others. These data highlight the impact of naturally occurring variation on the structure of the HPV capsid proteins of vaccine-relevant oncogenic HPV genotypes.

Functional evaluation of human papillomavirus type 31 long control region variants.

Persistent infections by high-risk human papillomavirus (HR-HPV) are a necessary condition, but not sufficient for cervical cancer development. Genetic variants of HR-HPV appear to be related to the risk of persistent infections. The study performed a functional evaluation of variants of the HPV-31 promoter region (LCR). For this, cloning and subcloning of variants HPV-31/UFPE-21 HPV-31/UFPE-89, HPV-31/UFPE-66, E2 gene and prototype HPV-31 were performed. Transfection with different concentrations of E2 was done and the concentration of 25 ng was determined to be ideal for LCR activation. HPV-31/UFPE-21 and HPV-31/UFPE-89 have a greater ability to alter Nluc reporter gene expression levels and HPV-31/UFPE-66 showed decreased levels of gene expression of Nluc reporter gene compared to control. Statistical analysis showed a significant difference between the polymorphic LCR regions and the control (p < 0.0001). A more refined profile of variants of HPV-31 and its importance for the prognosis of cervical lesions begins to be drawn.

Human Papillomavirus 31 Tyrosine 102 Regulates Interaction with E2 Binding Partners and Episomal Maintenance.

Several serine and threonine residues of the papillomavirus early E2 protein have been found to be phosphorylated. In contrast, only one E2 tyrosine phosphorylation site in BPV-1 (tyrosine 102) and one in HPV-16/31 (tyrosine 138) site have been characterized. Between BPV-1 and HPV-31 E2, 8 of the 11 tyrosines are conserved in the N-terminal domain, suggesting that phosphorylation of tyrosines has an essential role in E2 biology. In this study, we examine the effect of Y102 phosphorylation on HPV-31 E2 biology. Y102 proteins mutated either to the potential phospho-mimetic glutamic acid (Y102E) or to the nonphosphorylated homologue phenylalanine (Y102F) remain nuclear; however, Y102E is more associated with the nuclear matrix fraction. This is consistent with the inability of Y102E to bind TopBP1. Both BPV-1 and HPV-31 Y102E are similar in that neither binds the C terminus of Brd4, but in all other aspects the mutant behaves differently between the two families of papillomaviruses. BPV-1 Y102E was unable to bind E1 and did not replicate in a transient in vitro assay, while HPV-31 Y102E binds E1 and was able to replicate, albeit at lower levels than wild type. To examine the effect of E2 mutations under more native-like infection conditions, a neomycin-selectable marker was inserted into L1/L2 of the HPV-31 genome, creating HPV-31neo. This genome was maintained in every cell line tested for at least 50 days posttransfection/infection. Y102E in both transfection and infection conditions was unable to maintain high episome copy numbers in epithelial cell lines.IMPORTANCE Posttranslational modifications by phosphorylation can change protein activities, binding partners, or localization. Tyrosine 102 is conserved between delta papillomavirus BPV-1 and alpha papillomavirus HPV-31 E2. We characterized mutations of HPV-31 E2 for interactions with relevant cellular binding partners and replication in the context of the viral genome.

Synergistically enhanced colorimetric molecular detection using smart cup: a case for instrument-free HPV-associated cancer screening.

Rationale: Early and accurate detection of disease is crucial for its prevention, identification, and treatment. However, most of disease diagnostics is still limited in clinical laboratories due to the need of complicated instruments and professional personnel. Herein, we reported a smartphone-based synergistically enhanced colorimetric method for molecular diagnostics in our point of care (POC) smart cup platform. Methods: A disposable microfluidic chip was developed for colorimetric loop-mediated isothermal amplification (LAMP) detection of multiple HPV DNA in our POC smart cup platform. The colorimetric detection takes advantage of synergistic effect of PPi4- and H+ ions, two byproducts of LAMP reaction. Color signal of LAMP assay was recorded and analyzed by our custom Android app (dubbed "Hue Analyzer"). Results: Our method not only significantly improves colorimetric readout, but also provides a 10-fold increase in detection sensitivity. It has been successfully applied for HPV-associated cancer screening with spiked saliva and clinical swab samples. Conclusion: The proposed POC diagnostic platform is completely compatible with other nucleic acid biomarkers and has great potential for personalized health monitoring and disease prevention.

Stratified risk of high-grade cervical disease using onclarity HPV extended genotyping in women, ≥25 years of age, with NILM cytology.

OBJECTIVES: Increasing evidence suggests that extended human papillomavirus (HPV) genotyping (beyond 16/18) is effective for risk stratification in women with normal cytology. This report provides extended genotyping results, using the BD Onclarity HPV Assay, for individual genotypes HPV16, 18, 31, 45, 51, and 52 ̶ and three pooled genotype results for HPV33/58, 35/39/68, and 56/59/66. METHODS: 27,037 women with normal cytology, ≥25 years, were enrolled into the Onclarity HPV trial during routine screening. Women positive for any HPV genotype were referred to colposcopy/biopsy. Hierarchical-ranked prevalence and risk values, associated with cervical intraepithelial neoplasia, grade 2 or worse (≥CIN2) or ≥CIN3, were calculated based on extended genotyping results. RESULTS: HPV 16 and 31 carried the highest risk for ≥CIN2 (11.6% and 12.1%, respectively) and ≥CIN3 (8.1% and 7.5%, respectively); these genotypes were the most prevalent in both ≥CIN2 (29.6% and 19.3%, respectively) and ≥CIN3 (43.7% and 22.5%, respectively). Of the other 12 genotypes, HPV 18, 33/58, and 52 comprised an intermediate risk band (≥CIN2 risk range: 4.9-6.8%; ≥CIN3 risk range: 3.9-5.0%). Genotypes 45, 51, 35/39/68, and 56/59/66 constituted the lowest risk band for both disease grades (≥CIN2 value risk range: 1.7-3.0%; ≥CIN3 value risk range: 1.2-3.6%). CONCLUSIONS: Extended genotyping stratifies risk for ≥CIN2/3 in the ≥25 year-old, normal cytology population. While baseline HPV 16/31 values exceeded the risk threshold for colposcopy referral, the management of women with normal cytology who were positive for the intermediate- or lower-risk genotypes may evolve based on refined risk estimates as well as clinical factors.

Displaying 31RG-1 peptide on the surface of HPV16 L1 by use of a human papillomavirus chimeric virus-like particle induces cross-neutralizing antibody responses in mice.

Current available human papillomavirus (HPV) vaccines are based on the major capsid protein L1 virus-like particles (VLPs), which mainly induce type-specific neutralizing antibodies against vaccine types. Continuing to add more types of VLPs in a vaccine raises the complexity and cost of production which remains the principal impediment to achieve broad implementation of HPV vaccines, particularly in developing regions. In this study, we constructed 16L1-31L2 chimeric VLP (cVLP) by displaying HPV31 L2 aa.17-38 on the h4 coil surface region of HPV16 L1, and assessed its immunogenicity in mouse model. We found that the cVLP adjuvanted with alum plus monophosphoryl lipid A could induce cross-neutralizing antibody responses against 16 out of 17 tested HPV pseudoviruses, and the titer against HPV16 was as high as that was induced by HPV16 L1VLP (titer > 105), more importantly, titers over 103 were observed against two HR-HPVs including HPV31 (titer, 2,200) and -59 (titer, 1,013), among which HPV59 was not covered by Gardasil-9, and medium or low titers of cross-neutralizing antibodies against other 13 tested HPV pseudoviruses were also observed. Our data demonstrate that 16L1-31L2 cVLP is a promising candidate for the formulation of broader spectrum HPV vaccines.

Genetic variability in E5, E6, E7 and L1 genes of human papillomavirus type 31.

Human papillomavirus (HPV) type 31 is an important pathogenic subtype associated with cervical cancer. The aims of the present study were to analyze E5, E6, E7 and L1 gene mutations of HPV­31 among females, and to elucidate the evolutionary associations between them. In total, 87 positive samples were collected. The E5, E6, E7 and L1 genes were amplified by polymerase chain reaction and sequenced. Subsequently, two phylogenetic trees were constructed from the nucleotide sequences of the E5, E6 and E7 and the L1 variants of HPV­31. In total, 31 mutation sites of E5, E6 and E7 genes were identified, of which 16 were non­synonymous. T4053A (F80I), C285T (H60Y), C520T (A138V) and A743G (K62E) were the most common non­synonymous mutations. A total of 30 mutation sites of L1 genes were identified, of which four were non­synonymous. The most common non­synonymous mutations of L1 genes were A6350G (T29A) and C6372A (T36N). By phylogenetic analysis, A and C variants were most frequently detected, while B variants were less frequently detected in this population. The sequence variation data obtained in the present study provides a foundation for future research regarding HPV­induced oncogenesis, and may prove valuable for developing diagnostic probes and in the design of HPV vaccines for targeted populations.

Human Papillomavirus 16 (HPV-16), HPV-18, and HPV-31 E6 Override the Normal Phosphoregulation of E6AP Enzymatic Activity.

The human papillomavirus (HPV) E6 oncoproteins recruit the cellular ubiquitin ligase E6AP/UBE3A to target cellular substrates for proteasome-mediated degradation, and one consequence of this activity is the E6 stimulation of E6AP autoubiquitination and degradation. Recent studies identified an autism-linked mutation within E6AP at T485, which was identified as a protein kinase A phosphoacceptor site and which could directly regulate E6AP ubiquitin ligase activity. In this study, we have analyzed how T485-mediated regulation of E6AP might affect E6 targeting of some of its known substrates. We show that modulation of T485 has no effect on the ability of E6 to direct either p53 or Dlg for degradation. Furthermore, T485 regulation has no effect on HPV-16 or HPV-31 E6-induced autodegradation of E6AP but does affect HPV-18 E6-induced autodegradation of E6AP. In cells derived from cervical cancers, we find low levels of both phosphorylated and nonphosphorylated E6AP in the nucleus. However, ablation of E6 results in a dramatic accumulation of phospho-E6AP in the cytoplasm, whereas nonphosphorylated E6AP accumulates primarily in the nucleus. Interestingly, E6AP phosphorylation at T485 confers association with 14-3-3 proteins, and this interaction seems to be important, in part, for the ability of E6 to recruit phospho-E6AP into the nucleus. These results demonstrate that HPV E6 overrides the normal phosphoregulation of E6AP, both in terms of its enzymatic activity and its subcellular distribution.IMPORTANCE Recent reports demonstrate the importance of phosphoregulation of E6AP for its normal enzymatic activity. Here, we show that HPV E6 is capable of overriding this regulation and can promote degradation of p53 and Dlg regardless of the phosphorylation status of E6AP. Furthermore, E6 interaction with E6AP also significantly alters how E6AP is subject to autodegradation and suggests that this is not a simple stimulation of an already-existing activity but rather a redirection of E6AP activity toward itself. Furthermore, E6-mediated regulation of the subcellular distribution of phospho-E6AP appears to be dependent, in part, upon the 14-3-3 family of proteins.

Stop codon mutagenesis for homogenous expression of human papillomavirus L1 protein in Escherichia coli.

Human papillomavirus (HPV) is widely accepted to be the major causative pathogen of cervical cancer, warts, and other epithelial tumors. Virus infection and subsequent disease development can be prevented by vaccination with HPV vaccines derived from eukaryotic expression systems. Here, we report the soluble expression of the major capsid protein L1 of HPV31, a dominant carcinogenic HPV genotype, in Escherichia coli. HPV31 L1 protein and its elongated form (L1+) were observed in SDS-PAGE and CE-SDS analysis, generated by the native HPV31 L1 gene with a TAA stop codon. Replacing the TAA with TAG but not TGA could completely terminate protein translation. Mass spectrometry sequencing showed that L1+ comprised L1 with a C-terminal extension of 38 amino acids (aa). RNA folding analysis revealed that the unfaithful L1+ expression may result from translational read-through, as TAG is more stable and accessible than the other stop codons. The 38-aa elongated fragment perturbs self-assembly of HPV31 L1+, as shown in size and morphology analyses. By 3D cryo-electron microscopy structure determination, we show self-assembly of purified HPV31 L1 (TAG) VLPs into T = 7 icosahedral symmetry particles, resembling the native HPV virion. Finally, through additional characterization and antigenicity/immunogenicity assays, we verified that the E.coli-derived HPV31 VLPs are an ideal immunogen for HPV vaccine development. Our findings outline a codon optimization stratagem for protein expression and provide a method for the in-depth investigation of prokaryotic translation regulation.