Prevalence of human papillomavirus variants and genetic diversity in the L1 gene and long control region of HPV16, HPV31, and HPV58 found in North-East Brazil

Genetic variability of human papillomavirus type 18 based on E6, E7 and L1 genes in central China

BACKGROUND

High-risk human papillomavirus (HR-HPV) infection is an important factor for the development of cervical cancer. HPV18 is the second most common HR-HPV after HPV16.

METHODS

In this study, MEGA11 software was used to analyze the variation and phylogenetic tree of HPV18 E6-E7 and L1 genes. The selective pressure to E6, E7 and L1 genes was estimated using pamlX. In addition, the B cell epitopes of L1 amino acid sequences and T cell epitopes of E6-E7 amino acid sequences in HPV18 were predicted by ABCpred server and IEDB website, respectively.

RESULTS

A total of 9 single nucleotide variants were found in E6-E7 sequences, of which 2 were nonsynonymous variants and 7 were synonymous variants. Twenty single nucleotide variants were identified in L1 sequence, including 11 nonsynonymous variants and 9 synonymous variants. Phylogenetic analysis showed that E6-E7 and L1 sequences were all distributed in A lineage. In HPV18 E6, E7 and L1 sequences, no positively selected site was found. The nonconservative substitution R545C in L1 affected hypothetical B cell epitope. Two nonconservative substitutions, S82A in E6, and R53Q in E7, impacted multiple hypothetical T cell epitopes.

CONCLUSION

The sequence variation data of HPV18 may lay a foundation for the virus diagnosis, further study of cervical cancer and vaccine design in central China.

Identification and Functional Implications of the E5 Oncogene Polymorphisms of Human Papillomavirus Type 16

The persistence of the human papillomavirus type 16 (HPV16) infection on the cervical epithelium contributes to the progression of cervical cancer. Studies have demonstrated that HPV16 genetic variants may be associated with different risks of developing cervical cancer. However, the E5 oncoprotein of HPV16, which is related to several cellular mechanisms in the initial phases of the infection and thus contributes to carcinogenesis, is still little studied. Here we investigate the HPV16 E5 oncogene variants to assess the effects of different mutations on the biological function of the E5 protein. We detected and analyzed the HPV16 E5 oncogene polymorphisms and their phylogenetic relationships. After that, we proposed a tertiary structure analysis of the protein variants, preferential codon usage, and functional activity of the HPV16 E5 protein. Intra-type variants were grouped in the lineages A and D using in silico analysis. The mutations in E5 were located in the T-cell epitopes region. We therefore analyzed the interference of the HPV16 E5 protein in the NF-kB pathway. Our results showed that the variants HPV16E5_49PE and HPV16E5_85PE did not increase the potential of the pathway activation capacity. This study provides additional knowledge about the mechanisms of dispersion of the HPV16 E5 variants, providing evidence that these variants may be relevant to the modulation of the NF-κB signaling pathway.

Lineage Replacement and Genetic Changes of Four HR-HPV Types during the Period of Vaccine Coverage: A Six-Year Retrospective Study in Eastern China

OBJECTIVE

This study aimed to provide clinical evidence for lineage replacement and genetic changes of High-Risk Human Papillomavirus (HR-HPV) during the period of vaccine coverage and characterize those changes in eastern China.

METHODS

This study consisted of two stages. A total of 90,583 patients visiting the Obstetrics and Gynecology Hospital of Fudan University from March 2018 to March 2022 were included in the HPV typing analysis. Another 1076 patients who tested positive for HPV31, 33, 52, or 58 from November 2020 to August 2023 were further included for HPV sequencing. Vaccination records, especially vaccine types and the third dose administration time, medical history, and cervical cytology samples were collected. Viral DNA sequencing was then conducted, followed by phylogenetic analysis and sequence alignment.

RESULTS

The overall proportion of HPV31 and 58 infections increased by 1.23% and 0.51%, respectively, while infection by HPV33 and 52 decreased by 0.42% and 1.43%, respectively, within the four-year vaccination coverage period. The proportion of HPV31 C lineage infections showed a 22.17% increase in the vaccinated group, while that of the HPV58 A2 sublineage showed a 12.96% increase. T267A and T274N in the F-G loop of HPV31 L1 protein, L150F in the D-E loop, and T375N in the H-I loop of HPV58 L1 protein were identified as high-frequency escape-related mutations.

CONCLUSIONS

Differences in epidemic lineage changes and dominant mutation accumulation may result in a proportional difference in trends of HPV infection. New epidemic lineages and high-frequency escape-related mutations should be noted during the vaccine coverage period, and regional epidemic variants should be considered during the development of next-generation vaccines.

Vaccine escape challenges virus prevention: The example of two vaccine-preventable oncogenic viruses

Over the years, the pace of developing vaccines for HBV and HPV has never stopped. After more than 30 years of application, the HBV vaccine has reduced 80% of hepatocellular carcinoma (HCC). However, vaccine escape variants occur under selective pressure induced by widespread vaccination and antiviral therapy, which results in fulminant infection and horizontal transmission. Several mechanisms have been studied to explain HBV vaccine escape, including vaccine escape mutations (VEMs) in the major hydrophilic region, which leads to a decrease in the binding ability to neutralize antibodies and is the primary escape mechanism, protein conformational and N-linked glycosylation sites changes caused by VEMs, differences in genotype distribution, gene recombination, and some temporarily unknown reasons. However, effective solutions are still being explored. The HPV vaccine has also been proven to prevent 70%-90% of cervical cancer worldwide. Cases of HPV infection after being vaccinated have been observed in clinical practice. However, few researchers have paid attention to the mechanism of HPV vaccine escape. Thus, we reviewed the literature on vaccine escape of both HBV and HPV to discuss the mechanism of the virus escaping from vaccine protection and possible solutions to this problem. We analyzed the gap between studies of HPV and HBV and made prospects for further research in HPV vaccine escape.

E7 Oncogene HPV58 Variants Detected in Northeast Brazil: Genetic and Functional Analysis

Cervical cancer is associated with persistent infections by high-risk Human Papillomavirus (HPV) types that may have nucleotide polymorphisms and, consequently, different oncogenic potentials. Therefore, this study aimed to evaluate the genetic variability and structural effects of the E7 oncogene of HPV58 in cervical scraping samples from Brazilian women. The study was developed with patients from hospitals in the metropolitan area of Recife, PE, Brazil. The most frequent HPV types were, in descending order of abundance, HPV16, 31, and 58. Phylogenetic analysis demonstrated that the isolates were classified into sublineages A2, C1, and D2. Two positively selected mutations were found in E7: 63G and 64T. The mutations G41R, G63D, and T64A in the E7 protein reduced the stability of the protein structure. Utilizing an NF-kB reporter assay, we observed a decrease in the NK-kB pathway activity with the HPV58-E7 variant 54S compared to the WT E7. The other detected E7 HPV58 variants presented similar NF-kB pathway activity compared to the WT E7. In this study, it was possible to identify mutations that may interfere with the molecular interaction between the viral oncoproteins and host proteins.

Genetic diversity and bioinformatic analysis in the L1 gene of HPV genotypes 31, 33, and 58 circulating in women with normal cervical cytology

BACKGROUND

HPV-31, -33, and -58, along with HPV-45 and -52, account for almost 11% of HPV-associated cancers. Our previous studies showed that after HPV-16 and -51, HPV-58 was common and HPV-31 was as frequent as HPV-18 among Iranian women with normal cytology. Hence, in this study, we aimed to investigate the intra-type variations in L1 genes of HPV-58, -31, and -33 to find the predominant lineages circulating in women with normal cytology.

METHODS

Complete coding sequencing of the L1 gene was amplified and nucleotide and amino acid sequences were compared to those of the references. The selective pressure on L1 protein and whether the variations of the L1 genes embed in L1 loops, or N-glycosylated sites were also investigated.

RESULTS

B1, A, and A1 (sub)lineages were common in the HPV-58, -33, and -31 samples, respectively. Ninety nucleotide mutations were observed. Twenty nine nucleotide changes corresponded to nonsynonymous substitutions in which seventeen mutations were located in L1 loops. Only one codon position in HPV-58 sequences was found as the positive selection. No difference was observed in N-glycosylation sites between reference and understudied amino acid sequences.

CONCLUSION

In the current study, we reported, for the first time, the (sub) lineages, amino acid, and genetic diversity in the L1 gene of circulating HPV-58, -33, and -31, in women with normal cytology, in Iran. Such studies can not only have epidemiological values, but also aid to set vaccination programs.

Genetic variation within the human papillomavirus type 16 genome is associated with oropharyngeal cancer prognosis

PURPOSE

A significant barrier to adoption of de-escalated treatment protocols for human papillomavirus-driven oropharyngeal cancer (HPV-OPC) is that few predictors of poor prognosis exist. We conducted the first large whole-genome sequencing (WGS) study to characterize the genetic variation of the HPV type 16 (HPV16) genome and to evaluate its association with HPV-OPC patient survival.

PATIENTS AND METHODS

A total of 460 OPC tumor specimens from two large United States medical centers (1980-2017) underwent HPV16 whole-genome sequencing. Site-specific variable positions [single nucleotide polymorphisms (SNPs)] across the HPV16 genome were identified. Cox proportional hazards model estimated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall survival by HPV16 SNPs. Harrell C-index and time-dependent positive predictive value (PPV) curves and areas under the PPV curves were used to evaluate the predictive accuracy of HPV16 SNPs for overall survival.

RESULTS

A total of 384 OPC tumor specimens (83.48%) passed quality control filters with sufficient depth and coverage of HPV16 genome sequencing to be analyzed. Some 284 HPV16 SNPs with a minor allele frequency ≥1% were identified. Eight HPV16 SNPs were significantly associated with worse survival after false discovery rate correction (individual prevalence: 1.0%-5.5%; combined prevalence: 15.10%); E1 gene position 1053 [HR for overall survival (HRos): 3.75, 95% CI 1.77-7.95; Pfdr = 0.0099]; L2 gene positions 4410 (HRos: 5.32, 95% CI 1.91-14.81; Pfdr = 0.0120), 4539 (HRos: 6.54, 95% CI 2.03-21.08; Pfdr = 0.0117); 5050 (HRos: 6.53, 95% CI 2.34-18.24; Pfdr = 0.0030), and 5254 (HRos: 7.76, 95% CI 2.41-24.98; Pfdr = 0.0030); and L1 gene positions 5962 (HRos: 4.40, 95% CI 1.88-10.31; Pfdr = 0.0110) and 6025 (HRos: 5.71, 95% CI 2.43-13.41; Pfdr = 0.0008) and position 7173 within the upstream regulatory region (HRos: 9.90, 95% CI 3.05-32.12; Pfdr = 0.0007). Median survival time for patients with ≥1 high-risk HPV16 SNPs was 3.96 years compared with 18.67 years for patients without a high-risk SNP; log-rank test P < 0.001. HPV16 SNPs significantly improved the predictive accuracy for overall survival above traditional factors (age, smoking, stage, treatment); increase in C-index was 0.069 (95% CI 0.019-0.119, P < 0.001); increase in area under the PPV curve for predicting 5-year survival was 0.068 (95% CI 0.015-0.111, P = 0.008).

CONCLUSIONS

HPV16 genetic variation is associated with HPV-OPC prognosis and can improve prognostic accuracy.

Worldwide genetic variations in high-risk human papillomaviruses capsid L1 gene and their impact on vaccine efficiency

BACKGROUND

Human papillomavirus is the most common sexually transmitted infection. It is associated with different cancers, mainly cervical cancer, which remains the fourth most frequent cancer among women worldwide; it is also related to anogenital (anus, vulvar, vagina, and penis) and oropharyngeal cancers. Vaccination against HPV infection is the major way of prevention, and it has demonstrated impressive efficacy in reducing cervical cancer incidence. Nowadays, all the licensed HPV recombinant vaccines were designed based on HPV major capsid L1 protein. However, some variations in the HPV L1 gene sequence may induce structural changes within the L1 protein, which may alter the affinity and interaction of monoclonal antibodies (MAbs) with L1 protein epitopes, and influence host immune response and recognition. Hence, the importance of accuracy in delineating epitopes relevant to vaccine design and defining genetic variations within antigenic regions in the L1 gene to predict its impact on prophylactic vaccine efficiency. The present review reports the sequence variations in HR-HPV L1 gene isolates from different countries around the world, which may help to understand the effect of HPV L1 gene variations on vaccine efficiency.

METHODS

Research studies were retrieved from PubMed, Google Scholar, Science direct, and the National Center for Biotechnology Information (NCBI) database. A total of 31 articles describing genetic variations within the major capsid L1 gene and conducted in Africa, Europe, America and Asia were found. Only 26 studies conducted on HPV16, 18, 31, 33, 58, 45 and 52 which are the targets of HPV prophylactic vaccines, and which reported genetic variations within the L1 gene, were selected and evaluated in this review.

FINDINGS

We found a total of 87, 49, 11, 7, 22, 3, and 17 non-synonymous single nucleotide polymorphisms (SNPs) within HPV16, HPV18, HPV31, HPV58, HPV45, and HPV52 L1 gene, respectively. Four mutations were frequently observed in HPV16 L1 sequences: T353P in the HI loop, H228D in the EF loop, T266A in the FG loop, and T292A in the FG loop. Two mutations in HPV58 L1 sequences: T375N in the HI loop and L150F in the DE loop. Three mutations in HPV33 L1 sequences: T56N in the BC loop, G133S in the DE loop, T266K in the FG loop. Other mutations were found in HPV18, HPV45, and HPV52 L1 sequences. Some were found in different countries, and others were specific to a given population. Furthermore, some variations were located on peptide binding epitopes and lead to a modification of epitopes, which may influence MAbs interactions. Others need further investigations due to the lack of studies.

CONCLUSION

This study investigated the major capsid L1 genetic diversity of HPV16, 18, 31, 33, 58, 45, and 52 circulating in different populations around the world. Further investigations should be conducted to confirm their effect on immunogenicity and prophylactic vaccine efficiency.

Epidemiology and Burden of Human Papillomavirus and Related Diseases, Molecular Pathogenesis, and Vaccine Evaluation

Diagnosed in more than 90% of cervical cancers, the fourth deadliest cancer in women, human papillomavirus (HPV) is currently the most common pathogen responsible for female cancers. Moreover, HPV infection is associated with many other diseases, including cutaneous and anogenital warts, and genital and upper aerodigestive tract cancers. The incidence and prevalence of these pathologies vary considerably depending on factors including HPV genotype, regional conditions, the study population, and the anatomical site sampled. Recently, features of the cervicovaginal microbiota are found to be associated with the incidence of HPV-related diseases, presenting a novel approach to identify high-risk women through both blood and cervical samples. Overall, the HPV repartition data show that HPV infection and related diseases are more prevalent in developing countries. Moreover, the available (2-, 4-, and 9-valent) vaccines based on virus-like particles, despite their proven effectiveness and safety, present some limitations in terms of system development cost, transport cold chain, and oncogenic HPV variants. In addition, vaccination programs face some challenges, leading to a considerable burden of HPV infection and related diseases. Therefore, even though the new (9-valent) vaccine seems promising, next-generation vaccines as well as awareness programs associated with HPV vaccination and budget reinforcements for immunization are needed.

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.

Prevalence of oncogenic Human papillomavirus and genetic diversity in the L1 gene of HPV16 HPV 18 HPV31 and HPV33 found in women from Vojvodina Province Serbia

The main purpose of this paper is to estimate the pre-vaccination prevalence of 12 hrHPV types among 564 women from Vojvodina province (Serbia). The corrected contingency coefficient (Ccorr) was used to estimate the importance of association of examined HPV types and cytological diagnosis. The highest association with the abnormal cytology was observed for HPV 16 (Ccorr = 0.493) in all age groups of participants. The effect of HPV 16 was especially clear within the group of women older than 35 years (Ccorr = 0.691), compared with women younger than 35 (Ccorr = 0.333). The molecular characterization at the level of L1 gene of HPV 16, 18, 31 and 33 variants was for the first time assessed in our region. Nearly all HPV 16 isolates cluster with variant lineage A (96.4%) the remaining isolates clustering with variant lineage D. All of HPV 18 and HPV 33 isolates are clustering within the lineage A while isolates of HPV 31 group with lineages A and C. This contributes to understanding of intrinsic geographical and biological differences of examined HPV types and could be useful for development of cervical cancer screening strategies in Vojvodina (Serbia) and diagnosis of HPV related cervical cancer in general.

Identification of human papillomavirus type 16 variants circulating in the Calabria region by sequencing and phylogenetic analysis of HPV16 from cervical smears

Sequence analysis of HPV16 isolates reveals the presence of genome variants with characteristic mutations. The HPV16 variants have different geographical distribution and diverge into four phylogenetic lineages (A, B, C and D) and 16 sub-lineages: A1, A2, A3 (previously known as European variants), A4 (Asian variant), B1, B2, B3, B4, C1, C2, C3, and C4 (African variants), D1 (North-American variant), D2, D3 (Asian-American variants) and D4. Population studies showed that infections with viruses belonging to specific HPV16 sublineages confer different risks of viral persistence and cancer. In this study, 39 HPV16-positive cervical smears from European women living in Calabria (Italy) were analyzed for the presence of HPV16 variants. Cervical DNA extracts were processed by PCR to amplify L1, the Long Control Region (LCR), E6 and E7, which were sequenced. The sequences were concatenated and the 3169 nucleotides long fragments were characterized by BLAST and phylogenetic analysis. A total of 96 Single Nucleotide Polymorphism (SNPs) were detected, 29 of which mapping in the L1, 45 in the LCR, 15 in the E6 and 7 in the E7. The most common SNP was the T350G (29/39 samples, 74.4%), causing the L83 V amino acid change in the E6. Most of the HPV16 isolates (89.7%) had 99% of nucleotide (nt) identity to members of the A1 and A2 sublineages, while 4 isolates had 99% nt identity to members of the B2, B4, C1 and D4 sublineages. In conclusion, viruses belonging to the A1, A2, B2, B4, C1 and D4 HPV16 sublineages were found to circulate in the Calabria region.

Follow-up study of HPV58 variants in women with incident HPV58 infection from a Colombian cohort

Certain variants of human papillomavirus (HPV)type 58 are associated with an increased risk of high grade squamous intraepithelial lesions and cervical cancer. However, little is known about the persistence of HPV58 E6/E7 variants in women with incident HPV58 infections. The aim of the present study was to evaluate the presence and persistence of HPV58 E6/E7 variants in 71 women with incident HPV58 infection throughout their follow-up. These women belonged to a cohort examined in a longitudinal study of 1,610 Colombian women, who were HPV-negative and had normal baseline cytology. E6/E7 DNA regions of HPV58-positive samples were amplified and sequenced using automated direct sequencing. A total of 639 samples were analyzed from the 71 women, and 117 samples (18.3%) were HPV58-positive. HPV58 E6/E7 variants were detected in 85.5% of the samples. The T307/A694/G744/A761 variant was identified in 88% of the samples, the T307/G744 variant was identified in 9% of samples and the T187/T307/A367/G744/G793/T798/A801/T840/C852 was identified in 3% of the samples. Overall, 50% of the HPV58 infections were present after 1 year of follow-up and all infections were cleared after 7 years. Women who had first sexual intercourse at >15 years of age had a lower clearance rate than those who had sexual intercourse for the first time at ≤15 years of age [hazard ratio (HR)=0.29; 95% confidence interval (CI)=0.09-0.92]. Likewise, parous women had a higher clearance rate than nulliparous women (HR=3.43, 95% CI=1.23-9.60). There was no difference in clearance rates between HPV58 E6/E7 variants. In conclusion, HPV58 variants were not associated with persistence of the infection in this group of women.

Human papillomavirus detection using PCR and ATR-FTIR for cervical cancer screening

The human papillomavirus (HPV) genital infection is considered one of the most common sexually transmitted diseases worldwide, and has been associated with cervical cancer. The objective of this study was to investigate the efficacy of the diagnostic methods: polymerase chain reaction (PCR) and Fourier transform infrared (FTIR) equipped with an ATR (Attenuated Total Reflectance) unit (Pike Tech) spectroscopy, to diagnose HPV infection in women undergoing gynecological examination. Seventeen patients (41.46%) of the 41 patients analyzed were diagnosed with exophytic/condyloma acuminate lesions by clinical analysis, 29 patients (70.7%) (G1 group) of the 41 patients, showed positive result for HPV cell injury by oncotic colpocitology and 12 patients (29.3%) (G2 group), presented negative result for cellular lesion and absence of clinical HPV lesion. Four samples were obtained per patient, which were submitted oncotic colpocitology analysis (Papanicolau staining, two samples), PCR (one sample) and ATR-FTIR analysis (one sample). L1 gene was amplified by PCR technique with specific GP5+/GP6+ and MY09/MY11 primers. PCR results were uniformly positive for presence of HPV in all analyzed samples. Multivariate analysis of ATR-FTIR spectra suggests no significant biochemical changes between groups and no clustering formed, concurring with results of PCR. This study suggests that PCR and ATR-FTIR are highly sensitive technique for HPV detection.

Analysis of human papillomavirus 16 E6, E7 genes and Long Control Region in cervical samples from Uruguayan women

OBJECTIVE

This study aims to investigate the HPV16 variant distribution by sequence analyses of E6, E7 oncogenes and the Long Control Region (LCR), from cervical cells collected from Uruguayan women, and to reconstruct the phylogenetic relationships among variants.

METHODS

Forty-seven HPV16 variants, obtained from women with HSIL, LSIL, ASCUS and NILM cytological classes were analyzed for LCR and 12 were further studied for E6 and E7. Detailed sequence comparison, genetic heterogeneity analyses and phylogenetic reconstruction were performed.

RESULTS

A high variability was observed among LCR sequences, which were distributed in 18 different variants. E6 and E7 sequences exhibited novel non-synonymous substitutions. Uruguayan sequences mainly belonged to the European lineage, and only 5 sequences clustered in non-European branches; 3 of them in the Asian-American and North-American linage and 2 in an African branch. Additionally, 6 new variants from European and African clusters were identified.

CONCLUSIONS

HPV16 isolates mainly belonged to the European lineage, though strains from African and Asian-American lineages were also identified. Herein is reported for the first time the distribution and molecular characterization of HPV16 variants from Uruguay, providing novel insights on the molecular epidemiology of this infectious disease in the South America.

SYNOPSIS

A high variability among HPV 16 isolates mainly belonged to European lineage, provides an extensive sequence dataset from a country with high burden of cervical cancer.

Different Association of Human Papillomavirus 16 Variants with Early and Late Presentation of Cervical Cancer

The median age of cervical cancer (CC) presentation coincides with the mean age of menopause presentation (49 years) in Mexico. Here, we investigated the association between different HPV16 variants and early (≤ 49 years) or delayed (≥ 50 years) CC presentation. We conducted a case-case study that included 462 CCs, 386 squamous cell carcinomas (SCC), 63 adenocarcinomas (ACC), and 13 additional cell types. Variants were identified by PCR and DNA sequencing. The risk conferred by each variant for developing CC earlier than 50 years was analyzed using a univariate logistic regression model considering old-aged patients (≥ 50 years) and non-HPV16 cases as the reference variables. Overall, the frequency of HPV16 was 50.9%, and the only identified variants were the European A1/2 (31.2%) and the Asian-American D2 (10.8%), and D3 (8.9%). D2 was mainly associated with ≤ 49-year-old patients (15.9%); A1/2 was uniformly distributed between the two age groups (~31%), whereas D3 increased with age to a frequency of 11.8% in the older group. Only the D2 variant conferred a 3.3-fold increase in the risk of developing CC before 50 years of age (OR = 3.3, 95% CI = 1.7-6.6, p < 0.001) in relation with non-HPV16 cases. Remarkably, this risk was higher for ACC (OR = 6.0, 95% CI = 1.1-33, p < 0.05) than for SCC (OR = 2.8, 95% CI = 1.3-5.9, p < 0.01). Interestingly, when analyzing only the HPV16-positive CC, D2 increases (OR = 2.5, 95% CI = 1.2-5, p < 0.05) and D3 decreases (OR = 0.45, 95% CI 0.2-0.9, p < 0.05) the risk to develop CC before 50 years old in relation with A1/2 variant. These results indicated that D2 variant is associated with early and D3 with delayed CC presentation, whereas A1/2 variant was uniformly distributed between the two age groups.

L1 and L2 gene polymorphisms in HPV-58 and HPV-33: implications for vaccine design and diagnosis

Background

Cervical cancer is associated with infection by certain subtypes of human papillomavirus (HPV). The L1 protein comprising HPV vaccine formulations elicits high-titre neutralizing antibodies and confers protection against specific HPV subtypes. HPV L2 protein is an attractive candidate for cross-protective vaccines. HPV-33 and HPV-58 are very prevalent among Chinese women.

Methods

To study the gene intratypic variations and polymorphisms of HPV-33 and HPV-58 L1/L2 in Sichuan China, HPV-33 and HPV-58 L1 and L2 genes were sequenced and compared with other genes submitted to GenBank. Phylogenetic trees were constructed by maximum-likelihood and the Kimura 2-parameters methods (MEGA 6). The secondary structure was analyzed by PSIPred software, and HPV-33 and HPV-58 L1 homology models were created by SWISS-MODEL software. The selection pressures acting on the L1/L2 genes were estimated by PAML 4.8.

Results

Among 124 HPV-33 L1 sequences 20 single nucleotide mutations were observed included 8/20 non-synonymous and 12/20 synonymous mutations. The 101 HPV-33 L2 sequences included 12 single nucleotide mutations comprising 7/12 non-synonymous and 5/12 synonymous mutations. The 223 HPV-58 L1 sequences included 32 single nucleotide mutations comprising 9/32 non-synonymous and 23/32 synonymous mutations. The 201 HPV-58 L2 sequences comprised 26 single nucleotide mutations including 9/26 non-synonymous and 17/26 synonymous mutations. Selective pressure analysis showed that most of the common non-synonymous mutations showed a positive selection. HPV-33 and HPV-58 L2 were more stable than HPV-33 and HPV-58 L1.

Conclusions

HPV-33 and HPV-58 L2 were better candidates as clinical diagnostic targets compared with HPV-33 and HPV-58 L1. Clinical diagnostic probes and second-generation polyvalent vaccines should be designed on the basis of the unique sequence of HPV-33 and 58 L1/L2 variations in Sichuan, to improve the accuracy of clinical detection and the protective efficiency of vaccines.

Electronic supplementary material

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

Association Study between Cervical Lesions and Single or Multiple Vaccine-Target and Non-Vaccine Target Human Papillomavirus (HPV) Types in Women from Northeastern Brazil

We performed an association between high-grade squamous intraepithelial lesions (HSIL), low-grade squamous intraepithelial lesions (LSIL) and single or multiple vaccine-target as well as non-vaccine target Human papillomavirus (HPV) types. Using bead-based HPV genotyping, 594 gynecological samples were genotyped. An association between squamous intraepithelial lesion (SIL) and presence of HPV16, 18, 31, 58 and 56 types were calculated. The risk was estimated by using odds ratio (OR) and 95% of confidence intervals (CI). A total of 370 (62.3%) women were HPV positive. Among these, 157 (42.7%) presented a single HPV infection, and 212 (57.3%) were infected by more than one HPV type. HPV31 was the most prevalent genotype, regardless single and multiple HPV infections. Single infection with HPV31 was associated with LSIL (OR=2.32; 95%CI: 1.01 to 5.32; p=0.04); HPV31 was also associated with LSIL (OR=3.28; 95%CI: 1.74 to 6.19; p= 0.0002) and HSIL (OR=3.82; 95%CI: 2.10 to 6.97; p<0.001) in multiple HPV infections. Risk to harbor cervical lesions was observed in multiple HPV infections with regard to the HPV56 (OR=5.39; 95%CI: 2.44 to 11.90; p<0.001for LSIL; OR=5.37; 95%CI: 2.71 to 10.69; p<0.001) and HPV58 (OR=3.29; 95%CI: 1.34 to 8.09; p=0.0091 for LSIL; OR=3.55; 95%CI: 1.56 to 8.11; p=0.0026) genotypes. In addition, women coinfected with HPV16/31/56 types had 6 and 5-fold increased risk of HSIL (OR=6.46; 95%CI: 1.89 to 22.09; p=0.002) and LSIL (OR=5.22; 95%CI: 1.10 to 24.70; p=0.03), respectively. Multiple HPV infections without HPV16/18 has 2-fold increased risk of HSIL (OR=2.57; 95%CI: 1.41 to 4.70; p=0.002) and LSIL OR=2.03; 95%CI: 1.08 to 3.79; p=0.02). The results of this study suggest that single and multiple vaccine target as well as non-vaccine target HPV types are associated with LSIL and HSIL. These finding should be taken into consideration in the design of HPV vaccination strategies.