Genetic variability in E6, E7, and L1 genes of human papillomavirus genotype 52 from Southwest China

Genetic variation of E6 and E7 genes of human papillomavirus type 16 from central China

BACKGROUND

Persistent high-risk human papillomavirus (HR-HPV) infection is an important factor in the development of cervical cancer, and human papillomavirus type 16 (HPV-16) is the most common HR-HPV type worldwide. The oncogenic potential of HPV-16 is closely related to viral sequence variation.

METHODS

In order to clarify the variant characteristics of HPV-16 E6 and E7 genes in central China, E6 and E7 sequences of 205 HPV-16 positive samples were amplified by polymerase chain reaction. PCR products of E6 and E7 genes were further sequenced and subjected to variation analysis, phylogenetic analysis, selective pressure analysis and B-cell epitope prediction.

RESULTS

Twenty-six single nucleotide variants were observed in E6 sequence, including 21 non-synonymous and 5 synonymous variants. Twelve single nucleotide variants were identified in E7 sequence, including 6 non-synonymous and 6 synonymous variants. Four new variants were found. Furthermore, nucleotide variation A647G (N29S) in E7 was significantly related to the higher risk of HSIL and cervical cancer. Phylogenetic analysis showed that the E6 and E7 sequences were all distributed in A lineage. No positively selected site was found in HPV-16 E6 and E7 sequences. Non-conservative substitutions in E6, H31Y, D32N, D32E, I34M, L35V, E36Q, L45P, N65S and K75T, affected multiple B-cell epitopes. However, the variation of E7 gene had little impact on the corresponding B-cell epitopes (score < 0.85).

CONCLUSION

HPV-16 E6 and E7 sequences variation data may contribute to HR-HPV prevention and vaccine development in Jingzhou, central China.

Epidemiological analysis of HPV in Sichuan during 2014-2021

INTRODUCTION

Cervical cancer is a common malignancy among woman, strong molecular epidemiological data show that high risk HPV infection is the main cause of cervical cancer.

MATERIAL AND METHODS

Samples were collected from Sichuan women's and children's hospital based on the relevant guidelines and regulations, HPV DNA was extracted and evaluated by Human Papillomavirus Genotyping Kit for 21 types, according to the manufacturer's guidelines to analyze the epidemic age, mixed infection types, variation trend of HPV types in Sichuan from 2014 to 2021; Results: Out of 51174 samples11165 (21.82 %) HPV positive samples were detected, all belonging to alpha family, 53.32 % HPV positive samples and 61.51 % high-risk (HR) HPV positive samples are alpha-9 genus; The three commonest HR were HPV-52, HPV-16, HPV-58, and the low-risk (LR) HPV were HPV-81, HPV-6, HPV-11; Single infection was absolutely predominant and the age group with the highest HPV detection rate was 26-30 years old. During 2014-2021, HPV-16, HPV-6 and HPV-11 decline, while HPV-58 and HPV-52 increased; Conclusions: The most prevalent age group of HPV in this region was 26-30 years old. The detection rate of HPV-52 increased in the region, overtaking HPV-16 as the commonest type of HPV. α-9 genus HPV with strong pathogenicity is the commonest HR HPV. HPV prevalence systematic comparison in certain areas and continuous time can accurately and intuitively understand its distribution changes, achieve analysis of the epidemic trend, and provide guidance for the prevention, treatment and scientific research of HPV in Sichuan.

Genetic variability in the E6 and E7 oncogenes of HPV52 and its prevalence in the Taizhou area, China

Background

Human papillomavirus (HPV) 52 is one of the prevalent oncogenic HPV genotypes in East Asia. Chinese women have the highest susceptibility to the HPV52 type, but research data on HPV52 genetic variability and its carcinogenicity in China is lacking.

Methods

The present study aimed to investigate the genetic variability of HPV52 currently circulating among Chinese women by PCR sequencing the entire E6 and E7 oncogenes. HPV52 sequence alignment, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed by BioEdit software and MEGA X software.

Results

Between 2016 and 2018, the overall HPV infection rate was 21.3%, of which HPV52 was the most prevalent high-risk type (17.2%) in the Taizhou area, China. A total of 339 single HPV52-positive samples were included in this study. We obtained 27 distinct variation patterns of HPV52 with the accession GenBank numbers ON529577-ON529603. Phylogenetic analysis showed that 96.6% of HPV52 variants belonged to lineage B, which seemed to be uniquely defined by G350T, A379G (K93R) in the E6 gene and C751T, A801G in the E7 gene. Due to the dominance of lineage B in our study population, the results could not be used to assess the association of the HPV52 (sub)lineage with the risk of cervical lesions. In addition, no significant trends were observed between the nucleotide substitutions of HPV52 variants and the risk of cervical carcinogenesis.

Conclusion

Our data showed that HPV52 variants were strongly biased towards lineage B. These results confirmed that cervical lesions in the Taizhou area are highly attributable to HPV52, which may be due to the high infection rate of lineage B in the population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-022-01929-5.

Prevalence and distribution of human papillomavirus (HPV) in Luoyang city of Henan province during 2015-2021 and the genetic variability of HPV16 and 52

Background

Persistent high-risk Human papillomavirus (HPV) subtypes infection has been implicated as a causative of cervical cancer. Distribution and genotypes of HPV infection among females and their variations would assist in the formulation of preventive strategy for cervical cancer. The purpose of the present study is to investigate the prevalence of HPV among females in central China.

Methods

The distribution and genotypes of HPV among 9943 females attending the gynecological examinations in central of China during 2015–2021 were investigated. HPV genotypes were detected using a commercial kit. Nucleotides sequences of L1, E6 and E7 genes in HPV16 or HPV52 positive samples collected in 2021 were amplified by polymerase chain reaction (PCR). Variations of L1, E6 and E7 in HPV16 and HPV52 were gained by sequencing and compared with the reference sequence. Sublineages of HPV16 and HPV52 were determined by the construction of phylogenetic tree based on L1 gene.

Results

The overall prevalence of HPV infection was 22.81%, with the infection rate of high-risk human papillomavirus (HR-HPV) was 19.02% and low-risk human papillomavirus (LR-HPV) was 6.40%. The most top five genotypes of HPV infection were HPV16 (7.49%), HPV52 (3.04%), HPV58 (2.36%), HPV18 (1.65%) and HPV51 (1.61%). Plots of the age-infection rate showed that the single HPV, multiple HPV, HR-HPV, LR-HPV infection revealed the same tendency with two peaks of HPV infection were observed among females aged ≤ 20 year-old and 60–65 year-old. The predominant sublineage of HPV16 was A1 and B2 for HPV52. For HPV16, The most prevalent mutations were T266A (27/27) and N181T (7/27) for L1, D32E for E6 and S63F for E7 in HPV16. For HPV52, all of the nucleotide changes were synonymous mutation in L1 (except L5S) and E7 genes. The K93R mutation was observed in most HPV52 E6 protein.

Conclusions

The present study provides basic information about the distribution, genotypes and variations of HPV among females population in Henan province, which would assist in the formulation of preventive strategies and improvements of diagnostic probe and vaccine for HPV in this region.

Genetic variation of E6 and E7 genes of human papillomavirus 52 from Central China

Cervical cancer is the fourth most common malignant tumor in women worldwide and is closely related to human papillomavirus (HPV). Women have the highest susceptibility to HPV-52 type in Jingzhou, China. In this study, E6-E7 sequences of 183 HPV-52 positive samples were amplified by a polymerase chain reaction and sequenced. HPV-52 E6-E7 gene variations were analyzed. The phylogenetic tree was constructed using the Kimura 2-parameter method. The secondary structure of the protein was analyzed. The selective pressure to E6-E7 genes was estimated using PAML. In addition, the B cell epitopes of the E6-E7 sequences in HPV-52 were predicted by the ABCpred server. In E6 sequences, 15 single nucleotide variants were observed, including 6 nonsynonymous variants and 9 synonymous variants. In E7 sequences, 19 single nucleotide variants occurred, including 10 nonsynonymous variants and 9 synonymous variants. Six amino acid variants, including 3 nonconservative substitutions, were found in sequences encoding the alpha helix. Eight amino acid variants, including three nonconservative substitutions, occurred in sequences encoding the strand. Through phylogenetic analysis, the E6-E7 sequences were mainly distributed in B lineage. In HPV-52 E6-E7 sequences, no positively selected site was found. The nonconservative substitutions, such as K93R, K93E in E6, T37I, and D38N in E7, affected multiple hypothetical epitopes in the B cell. This study provides information for the investigation of HPV epidemic characters. The discovery of new variants of HPV-52 may lay the basis for the development of the virus diagnosis, further study of cervical cancer, and vaccine design in Central China.

The genetic variability, phylogeny and functional significance of E6, E7 and LCR in human papillomavirus type 52 isolates in Sichuan, China

Background

Variations in human papillomavirus (HPV) E6 and E7 have been shown to be closely related to the persistence of the virus and the occurrence and development of cervical cancer. Long control region (LCR) of HPV has been shown multiple functions on regulating viral transcription. In recent years, there have been reports on E6/E7/LCR of HPV-16 and HPV-58, but there are few studies on HPV-52, especially for LCR. In this study, we focused on gene polymorphism of the HPV-52 E6/E7/LCR sequences, assessed the effects of variations on the immune recognition of viral E6 and E7 antigens, predicted the effect of LCR variations on transcription factor binding sites and provided more basic date for further study of E6/E7/LCR in Chengdu, China.

Methods

LCR/E6/E7 of the HPV-52 were amplified and sequenced to do polymorphic and phylogenetic analysis. Sequences were aligned with the reference sequence by MEGA 7.0 to identify SNP. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED 4.0. The selection pressure of E6 and E7 coding regions were estimated by Bayes empirical Bayes analysis of PAML 4.9. The HLA class-I and II binding peptides were predicted by the Immune Epitope Database server. The B cell epitopes were predicted by ABCpred server. Transcription factor binding sites in LCR were predicted by JASPAR database.

Results

50 SNP sites (6 in E6, 10 in E7, 34 in LCR) were found. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. Two deletions were found between the nucleotide sites 7387–7391 (TTATG) and 7698–7700 (CTT) in all samples. A deletion was found between the nucleotide sites 7287–7288 (TG) in 97.56% (40/41) of the samples. The combinations of all the SNP sites and deletions resulted in 12 unique sequences. As shown in the neighbor-joining phylogenetic tree, except for one belonging to sub-lineage C2, others sequences clustered into sub-lineage B2. No positive selection was observed in E6 and E7. 8 non-synonymous amino acid substitutions (including E3Q and K93R in the E6, and T37I, S52D, Y59D, H61Y, D64N and L99R in the E7) were potential affecting multiple putative epitopes for both CD4+ and CD8+ T-cells and B-cells. A7168G was the most variable site (100%) and the binding sites for transcription factor VAX1 in LCR. In addition, the prediction results showed that LCR had the high probability binding sites for transcription factors SOX9, FOS, RAX, HOXA5, VAX1 and SRY.

Conclusion

This study provides basic data for understanding the relation among E6/E7/LCR mutations, lineages and carcinogenesis. Furthermore, it provides an insight into the intrinsic geographical relatedness and biological differences of the HPV-52 variants, and contributes to further research on the HPV-52 therapeutic vaccine development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01565-5.

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.

Intra-host sequence variability in human papillomavirus

Human papillomaviruses (HPVs) co-evolve slowly with the human host and each HPV genotype displays epithelial tropisms. We assessed the evolution of intra HPV genotype variants within samples, and their association to anogenital site, cervical cytology and HIV status. Variability in the L1 gene of 35 HPV genotypes was characterized phylogenetically using maximum likelihood, and portrayed by phenotype. Up to a thousand unique variants were identified within individual samples. In-depth analyses of the most prevalent genotypes, HPV16, HPV18 and HPV52, revealed that the high diversity was dominated by a few abundant variants. This suggests high intra-host mutation rates. Clades of HPV16, HPV18 and HPV52 were associated to anatomical site and HIV co-infection. Particularly, we observed that one HPV16 clade was specific to vaginal cells and one HPV52 clade was specific to anal cells. One major HPV52 clade, present in several samples, was strongly associated with cervical neoplasia. Overall, our data suggest that tissue tropism and HIV immunosuppression are strong shapers of HPV evolution.

Genetic variability and lineage phylogeny of human papillomavirus type-16 and -53 based on the E6, E7, and L1 genes in Southwest China

Human papillomaviruses (HPVs) are circular double-stranded DNA viruses that are highly prevalent in the general population, and account for the cervical cancer burden in women worldwide. In this study, we analyzed HPV-16, the most prevalent type worldwide, and HPV-53, a possible high-risk type from infected women in Southwest China. To characterize mutations, intratypic variants, and genetic variability in the E6, E7, and L1 genes of HPV-16 (n=97) and HPV-53 (n=15), these genes were sequenced and submitted to GenBank. Phylogenetic trees were constructed using Bayesian trees, followed by secondary structure analysis and B-cell epitope prediction. Moreover, the selection pressures of the E6, E7, and L1 genes were estimated. In total, 27 novel variants of HPV-16 and 11 novel variants of HPV-53 were identified. In the HPV-16 E6-E7-L1 sequences, 73 nucleotide changes were observed with 40/73 being non-synonymous mutations (two in the alpha helix and five in the beta sheet) and 33/73 being synonymous. In the HPV-53 E6-E7-L1 sequences, 64 nucleotide changes were observed with 26/64 being non-synonymous mutations (three in the alpha helix and one in the beta sheet) and 38/64 being synonymous. Selective pressure analysis showed that most of these mutations did not reflect positive selection. The maximal divergence between any two variants within each gene of these two HPV types ranging from 0.94%(HPV-16 L1 gene)to 2.80%(HPV-53 E6 gene). Identifying new variants of HPV-16 and -53 from women in Southwest China may be helpful to design vaccines specifically for women in Southwest China and testing methods specifically for this region. The results of our study may contribute to future researches in diagnostic probes, vaccines improvement, or screening methods for a particular population.