E6 and E7 Gene Polymorphisms in Human Papillomavirus Types-58 and 33 Identified in Southwest China

Distribution of HPV types among women with HPV-related diseases and exploration of lineages and variants of HPV 52 and 58 among HPV-infected patients in China: A systematic literature review

To summarize the distribution of types of human papillomavirus (HPV) associated with HPV-related diseases and investigate the potential causes of high prevalence of HPV 52 and 58 by summarizing the prevalence of lineages, sub-lineages, and mutations among Chinese women. We searched PubMed, EMBASE, CNKI, and WanFang from January, 2012 to June, 2023 to identify all the eligible studies. We excluded patients who had received HPV vaccinations. Data were summarized in tables and cloud/rain maps. A total of 102 studies reporting HPV distribution and 15 studies reporting HPV52/HPV58 variants were extracted. Among Chinese women, the top five prevalent HPV types associated with cervical cancer (CC) were HPV16, 18, 58, 52, and 33. In patients with vaginal cancers and precancerous lesions, the most common HPV types were 16 and 52 followed by 58. For women with condyloma acuminatum (CA), the most common HPV types were 11 and 6. In Chinese women with HPV infection, lineage B was the most prominently identified for HPV52, and lineage A was the most common for HPV58. In addition to HPV types 16, which is prevalent worldwide, our findings revealed the unique high prevalence of HPV 52/58 among Chinese women with HPV-related diseases. HPV 52 variants were predominantly biased toward lineage B and sub-lineage B2, and HPV 58 variants were strongly biased toward lineage A and sub-lineage A1. Further investigations on the association between the high prevalent lineage and sub-lineage in HPV 52/58 and the risk of cancer risk are needed. Our findings underscore the importance of vaccination with the nine-valent HPV vaccine in China.

Phylogenetic analysis and antigenic epitope prediction for E6 and E7 of Alpha-papillomavirus 9 in Taizhou, China

BACKGROUND

Alpha-papillomavirus 9 (α-9) is a member of the human papillomavirus (HPV) α genus, causing 75% invasive cervical cancers worldwide. The purpose of this study was to provide data for effective treatment of HPV-induced cervical lesions in Taizhou by analysing the genetic variation and antigenic epitopes of α-9 HPV E6 and E7.

METHODS

Cervical exfoliated cells were collected for HPV genotyping. Positive samples of the α-9 HPV single type were selected for E6 and E7 gene sequencing. The obtained nucleotide sequences were translated into amino acid sequences (protein primary structure) using MEGA X, and positive selection sites of the amino acid sequences were evaluated using PAML. The secondary and tertiary structures of the E6 and E7 proteins were predicted using PSIPred, SWISS-MODEL, and PyMol. Potential T/B-cell epitopes were predicted by Industrial Engineering Database (IEDB).

RESULTS

From 2012 to 2023, α-9 HPV accounted for 75.0% (7815/10423) of high-risk HPV-positive samples in Taizhou, both alone and in combination with other types. Among these, single-type-positive samples of α-9 HPV were selected, and the entire E6 and E7 genes were sequenced, including 298 HPV16, 149 HPV31, 185 HPV33, 123 HPV35, 325 HPV52, and 199 HPV58 samples. Compared with reference sequences, 34, 12, 10, 2, 17, and 17 nonsynonymous nucleotide mutations were detected in HPV16, 31, 33, 35, 52, and 58, respectively. Among all nonsynonymous nucleotide mutations, 19 positive selection sites were selected, which may have evolutionary significance in rendering α-9 HPV adaptive to its environment. Immunoinformatics predicted 57 potential linear and 59 conformational B-cell epitopes, many of which are also predicted as CTL epitopes.

CONCLUSION

The present study provides almost comprehensive data on the genetic variations, phylogenetics, positive selection sites, and antigenic epitopes of α-9 HPV E6 and E7 in Taizhou, China, which will be helpful for local HPV therapeutic vaccine development.

Cervical carcinogenesis risk association of HPV33 E6 and E7 genetic variations in Taizhou, Southeast China

BACKGROUND

Human papillomavirus (HPV) 33 belongs to the Alphapapillomavirus 9 (α-9 HPV) species group, which also contains types 16, 31, 35, 52, 58 and 67. The purpose of this study was to investigate the genetic variations of HPV33 and to explore its carcinogenicity among women in Taizhou, Southeast China.

METHODS

Exfoliated cervical cells were collected for HPV genotyping. Only single HPV33 infection cases were selected, and their E6 and E7 genes were sequenced using the ABI 3730xl sequencer and then analysed using MEGA X.

RESULTS

From 2014 to 2020, a total of 185 single HPV33-positive specimens were successfully amplified. We obtained 15 distinct HPV33 E6/E7 variants, which were published in GenBank under accession numbers OQ672665-OQ672679. Phylogenetic analysis revealed that all HPV33 E6/E7 variants belonged to lineage A, of which 75.7% belonged to lineage A1. Compared with CIN1, the proportion of sublineage A1 in CIN2/3 was higher, but there was no significant difference (76.5% vs. 80.6%, P > 0.05). Altogether, 20 single nucleotide substitutions were identified, of which 6 were novel substitutions, including T196G (C30G), A447T, G458T (R117L), G531A, A704A, and C740T. In addition, no significant trends were observed between the nucleotide substitutions of HPV33 E6/E7 variants and the risk of cervical lesions.

CONCLUSION

This study provides the most comprehensive data on genetic variations, phylogenetics and carcinogenicity of HPV33 E6/E7 variants in Southeast China to date. The data confirmed that cervical lesions among women in Taizhou are attributable to HPV33, which may be due to the high infection rate of sublineage A1 in the population.

Transcriptional activity of the long control region in human papillomavirus type 33 intratype variants

BACKGROUND

High-risk human papillomaviruses (HPVs) are responsible for the development of cervical and other anogenital cancers. Intratype sequence variants of certain high-risk HPV types (e.g. 16, 18 and 31) are thought to have different oncogenic potential, partly due to nucleotide sequence variation in the viral long control region (LCR). The LCR has an important role in the regulation of viral replication and transcription. The purpose of this study was to explore sequence variation in the LCR of HPV 33 intratype variants in Hungary and to see whether there are differences in the transcriptional activities of the variants.

METHODS

The complete HPV 33 LCR was amplified from HPV 33 positive cervical samples. After sequencing the LCR variants, multiple sequence alignment and phylogenetic analyses were carried out. Representative HPV 33 LCR sequence variants were selected for cloning and functional analysis. After transient transfection of HeLa cells, luciferase reporter assays were used to analyse the transcriptional activities of different LCR variants.

RESULTS

Altogether 10 different variants were identified by sequence analysis of the HPV 33 LCR. The results of phylogenetic analysis showed that 3 variants belonged to sublineage A1, while the other 7 variants clustered with sublineage A2. Variants belonging to sublineage A2 had significantly lower transcriptional activities than variants belonging to sublineage A1. Within sublineage A2, the two variants analysed had significantly different transcriptional activities, which was shown to be caused by the A7879G variation.

CONCLUSIONS

Nucleotide variation in the HPV 33 LCR can result in altered transcriptional activity of the intratype variants. Our results can help to understand the correlation between LCR polymorphism and the oncogenic potential of HPV 33 variants.

Characteristics of human papillomavirus infection among females and the genetic variations of HPV18 and HPV58 in Henan province, China

The present study aims to investigate the genotype distribution of Human papillomavirus (HPV) and variations of HPV18 and HPV58 infection among 6538 females in Luoyang city during 2019-2021. The overall positive rate of females with HPV infection was 12.34%, with 9.74% were infected with single HPV and 2.60% with multiple HPVs. The prevalent rate of high-risk HPV (HR-HPV) was 9.85% and the top five HR-HPV genotypes were HPV52 (1.94%), HPV16 (1.93%), HPV58 (1.48%), HPV51 (1.02%) and HPVV39 (0.99%). Two peaks of HPV infections rates were observed in females aged ≤ 20 and 61-65 years old. To characterize mutations, 39 HPV18 and 56 HPV58 L1, E6 and E7 genes were sequenced and submitted to GenBank. In the HPV18 E6-E7-L1 sequences, 38 nucleotides changes were observed with 10/38 were non-synonymous mutations (5 in E6 gene, 1 in E7 gene and 4 in L1 gene). In the HPV58 E6-E7-L1 sequences, 53 nucleotides changes were observed with 23/53 were non-synonymous mutations (3 in E6 gene, 5 in E7 gene and 15 in L1 gene). Phylogenetic analysis based on L1 gene showed that 92.3% (36/39) of HPV18 isolates fell into sublineage A1 and 7.7% (3/39) belonged to A5. For HPV58, 75.0% (42/56) isolates belonged to sublineage A1 and 25.0% (14/56) were sublineage A2. There was no association between amino mutation and cervical lesions. The present study provides basic information about the distribution, genotypes and variations of HPV among females population in Luoyang city, which would assist in the formulation of HPV screening and vaccination programs and preventive strategies for HPV-attributable cancer in this region.

Polymorphism of E6 and E7 Genes in Human Papillomavirus Types 31 and 33 in Northeast China

Persistent infection with human papillomavirus (HPV) types 31 and 33 is an important causative factor for cervical cancer. The E6/E7 genes are key oncogenes involved in the immortalization and transformation of human epithelial cells. Genetic polymorphism may lead to differences in the virus' carcinogenic potential, the immune reaction of the host, and the potencies of vaccines. Few studies on HPV31/33 E6/E7 genetic polymorphism have been carried out. To study the genetic polymorphism of HPV31 and HPV33 E6/E7 genes in northeast China, these genes (HPV31 E6/E7, n = 151; HPV33 E6/E7, n = 136) were sequenced and compared to reference sequences (J04353.1, M12732.1) using BioEdit. Phylogenetic trees were constructed by the neighbor-joining method using MegaX. The diversity of the secondary structure was estimated using the PSIPred server. The positively selected sites were analyzed using PAML4.9. The major histocompatibility complex (MHC) class I and MHCII epitopes were predicted using the ProPred-I server and ProPredserver. B-cell epitopes were predicted using the ABCpred server. In the 151 HPV31E6 sequences, 25 (25/450) single-nucleotide mutations were found, 14 of which were synonymous mutations and 11 were nonsynonymous. In the 151 HPV31E7 sequences, 8 (8/297) nucleotide mutations were found, 3 of which were synonymous mutations and 5 were nonsynonymous. In the 136 HPV33E6 sequences, 17 (17/450) nucleotide mutations were observed, 7 of which were synonymous mutations and 10 were nonsynonymous. C14T/G (T5I/S) was a triallelic mutation. Finally, in the 136 HPV33E7 sequences, 9 (9/294) nucleotide mutations were observed, 3 of which were synonymous mutations and 6 were nonsynonymous. C134T/A (A45V/E) and C278G/A (T93S/N) were triallelic mutations. Lineage A was the most common lineage in both HPV31 and HPV33. In all of the sequences, we only identified one positively selected site, HPV33 E6 (K93N). Most nonsynonymous mutations were localized at sites belonging to MHC and/or B-cell predicted epitopes. Data obtained in this study should contribute to the development and application of detection probes, targeted drugs, and vaccines.

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.

Mortality in COVID-19 disease patients: Correlating the association of major histocompatibility complex (MHC) with severe acute respiratory syndrome 2 (SARS-CoV-2) variants

Genetic factors such as the HLA type of patients may play a role in regard to disease severity and clinical outcome of patients with COVID-19. Taking the data deposited in the GISAID database, we made predictions using the IEDB analysis resource (TepiTool) to gauge how variants in the SARS-CoV-2 genome may change peptide binding to the most frequent MHC-class I and -II alleles in Africa, Asia and Europe. We caracterized how a single mutation in the wildtype sequence of of SARS-CoV-2 could influence the peptide binding of SARS-CoV-2 variants to MHC class II, but not to MHC class I alleles. Assuming the ORF8 (L84S) mutation is biologically significant, selective pressure from MHC class II alleles may select for viral varients and subsequently shape the quality and quantity of cellular immune responses aginast SARS-CoV-2. MHC 4-digit typing along with viral sequence analysis should be considered in studies examining clinical outcomes in patients with COVID-19.

Bioinformatics analysis of HPV-68 E6 and E7 oncoproteins for designing a therapeutic epitope vaccine against HPV infection

The incidence and mortality of cervical cancer, which mainly results from the infection of human papillomavirus (HPV) is significantly increasing in Xinjiang. According to the previous research, the incidence of HPV-68 in cervical cancer patients in Xinjiang is significantly higher than in other parts of China. HPV E6 and E7 oncoproteins play a crucial role in cervical cancer, and can be used as ideal targets for therapeutic vaccines. Therefore, we analyzed and identified the possible T-cell and B-cell dominant epitopes and various aspects of HPV-68 E6 and E7 oncoproteins, including the physicochemical properties, secondary and tertiary structures using a bioinformatic approach, which provided a basis for designing an effective HPV infection therapeutic vaccine. The results showed that E6 oncoproteins was an unstable and hydrophilic protein, while E7 oncoproteins was unstable and hydrophilic protein. The secondary structure of the E6 oncoproteins consisted of 45.57% alpha helixes, 14.56% extended strands, 4.43% beta turns and 35.44% random coils. The secondary structure of E7 oncoproteins consisted of 35.45% alpha helixes, 17.27% extended strands, 0.91% beta turns and 46.36% random coils. Moreover, our results identified 5 dominant T-cell epitopes and 6 dominant B-cell epitopes in the E6 oncoproteins structure and 5 dominant T-cell epitopes and 3 dominant B-cell epitopes in E7 oncoproteins. In conclusion, this study provides comprehensive biological information about the HPV-68 E6 and E7 oncoproteins, which will lay a theoretical foundation for multi-epitope vaccines against HPV infection.

Phylogeny and polymorphism in the E6 and E7 of human papillomavirus: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44) in women from Shanghai

Background

Persistent infection with human papillomaviruses (HPVs) has been associated with cervical intraepithelial neoplasia (CIN) and cervical cancer. However, why only a fraction of HPV cases progress to cancer is still unclear.

Methods

We focused on the heterogeneity, classification, evolution and dispersal of variants for 14 common HPV types in 262 HPV-positive patients with cervical lesions. The E6 and E7 genes of HPV were sequenced and compared with the HPV reference for sequence analysis. Phylogenetic trees were constructed using the neighbour-joining tree method with MEGA 7.0.

Results

In this study, 233 E6 and 212 E7 sequences were successfully amplified by PCR, and these sequences were divided into 5 species groups: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44). The incidence of high-grade squamous intraepithelial lesion (HSIL) in patients infected with alpha-9 HPV was significantly increased compared with other groups (P < 0.0001), especially HPV16 (P < 0.0001). Strikingly, E7 had significantly fewer nonsynonymous variants in the HSIL compared to <HSIL groups (P = 3.17× 10^- 4). The A388C (K93 N) variation in HPV58 E6 can significantly reduce the risk of HSIL (P = 0.015). However, T7220G (D32E) variation in HPV16 E6 and A7689G (N29S) in HPV16 E7 increased the incidence of HSIL compared to the <HSIL group (P = 0.036 and 0.022).

Conclusions

Strict conservation of E7 is important for HPV carcinogenicity, especially N29 of HPV16. The findings in this work provide preventative/therapeutic interventions for HPV infections and CIN.

Characteristics of human papillomaviruses distribution in Guizhou Province, China

BACKGROUND

Human papillomavirus (HPV) is one of the most common sexually transmitted viruses. Data about HPV infection in Guizhou is limited.

METHODS

56,768 cervical samples were collected and genotyped for 15 main high risk and 6 main low risk HPV types.

RESULTS

16.95% (9623/56768) of samples were HPV positive; 90.70% (8728/9623) of HPV positive women were infected by high risk HPV. High risk and high risk mix infection (1458; 70.85%) was the most common mix HPV infection type. The highest HPV detection rate was found in age group 41-45 years old (detection rate = 17.89%) (χ2 = 204.77; P < 0.001); the highest within-group HPV infection rates were found in the ≤20 (25.62%) and ≥ 61 (24.67%) years old age groups, the lowest within-group HPV infection rate was found in the 31-35 years old age group (15.02%). The highest mix infection proportions were found in the ≥61 (36.06%) and ≤ 20 (33.63%) years old age groups (χ2 = 111.21; P < 0.001), the lowest mix infection proportion was found in the 41-45 (17.42%) years old age group. The highest high risk infection proportions were found in the 26-30 (92.98%), ≥61 (92.68%), and 36-40 (92.16%) years old age groups (χ2 = 31.72; P < 0.001), the lowest high risk infection proportion was found in the ≤20 (84.96%) years old age group. HPV infection rates varied with seasons in Guizhou.

CONCLUSIONS

Characteristics of HPV distribution in Guizhou were identified. There were significant differences in HPV distribution among age groups, prevention strategies should be adjusted according to the characteristics.

E6 and E7 gene polymorphisms in human papillomavirus Type-6 identified in Southwest China

BACKGROUND

Human papillomavirus type-6 (HPV6) is the major etiological agent of anogenital warts both men and women. The present study aimed to characterize the genetic diversity among HPV6 in Southwest China, and to investigate the origin of, selective pressure experienced by, and impact of the resultantly identified genetic variants on the HPV6 secondary structure.

METHODS

Phylogenetic trees were constructed by Maximum-likelihood and the Kimura 2-parameters methods by Molecular Evolutionary Genetics Analysis version 6.0. The diversity of secondary structure was analyzed by PSIPred software. The selection pressures acting on the E6/E7 genes were estimated by Phylogenetic Analyses by Maximum Likelihood version 4.8 software.

RESULTS

HPV6 was the most prevalent low risk HPV type in southwest China. In total, 143 E6 and E7 gene sequences of HPV6 isolated from patients were sequenced and compared to GenBank HPV6 reference sequence X00203. The results of these analyses revealed that both the HPV6 E6 and E7 were highly conserved within the analyzed patient samples, and comprised only 3 types of variant sequence, respectively. Furthermore, the analysis of HPV6 E6 and E7 sequences revealed seven/five single-nucleotide mutations, two/four and five/one of which were non-synonymous and synonymous, respectively. The phylogenetic analyses of the E6 and E7 sequences indicated that they belonged to sub-lineage A1 and sub-lineage B1, whereas the selective pressure analyses showed that only the E7 mutation sites 4R, 34E, and 52F were positive selection.

CONCLUSIONS

HPV6 (detection rate = 13.10%) was very prevalent in southwest China, both the HPV6 E6 and E7 sequences were highly conserved within the analyzed patient samples in southwest China, indicating that the low risk HPV6 can adapt to the environment well without much evolution.

Identification of variants and therapeutic epitopes in HPV-33/HPV-58 E6 and E7 in Southwest China

BACKGROUND

Human papillomavirus (HPV) E6 and E7 oncoproteins play a crucial role in HPV-related diseases, such as cervical cancer, and can be used as ideal targets for therapeutic vaccines. Human leukocyte antigen (HLA) participates in the immune response to block HPV infection and invasion by its target/recognition function. HPV-33 and HPV-58 are highly prevalent among Chinese women. Therefore, it is of great significance to study the E6 and E7 region-specific gene polymorphisms of HPV-33 and HPV-58 in Southwest China and to identify ideal epitopes for vaccine design. Both HPV-33 and HPV-58 belong to α-9 genus HPV and are highly homologous, so their correlations are included in our research.

METHODS

To study the E6 and E7 variations and polymorphisms of HPV-33 and HPV-58 in Southwest China, we collected samples, extracted and sequenced DNA, and identified variants. Nucleotide sequences were translated into amino acids by Mega 6.0 software. The physical/chemical properties, amino acid-conserved sequences and secondary structure of protein sequences were analysed by the Protparam server, ConSurf server and PSIPRED software. The T and B cell epitopes of the E6/E7 reference and variant sequences in HPV-33 and HPV-58 were predicted by the Immune Epitope Database (IEDB) analysis server and the ABCpred server, respectively.

RESULTS

Five and seven optimal HLA-I restricted T cell epitopes were selected from HPV-33 and HPV-58 E6, respectively, and these optimal epitopes are mainly located in _41-58EVYDFAFADLTVVYREGN of HPV-33 E6 and _40-60SEVYDFVFADLRIVYRDGNPF of HPV-58 E6. Six optimal HLA-I-restricted T cell epitopes were selected from HPV-33 and HPV-58 E7, and these epitopes are mainly located in _77-90RTIQQLLMGTVNIV of HPV-33 E7 and _78-91RTLQQLLMGTCTIV of HPV-58 E7.

CONCLUSIONS

HPV-33/HPV-58 E6/E7 gene polymorphisms and T/B cell epitopes of their reference and variant sequences were studied, and candidate epitopes were selected by bioinformatics techniques for therapeutic vaccine design for people in Southwest China. This study was the first to investigate the correlation of epitopes between HPV-33 and HPV-58. After experimental validation, these selected epitopes will be employed to induce a wide range of immune responses in heterogeneous HLA populations.

Genetic variability and oncogenic risk association of human papillomavirus type 58 E6 and E7 genes in Taizhou area, China

It is well recognized that the association of human papillomavirus (HPV) and cervical carcinogenesis is based on the presence of HPV DNA sequence. The E6 and E7 oncoproteins encoded by high-risk HPV types play a key role in carcinogenesis. HPV58 type accounts for a larger share of cervical disease in China, whereas data on HPV58 genetic variability in China is limited. We aimed to evaluate the diversity of HPV58 genetic variants by sequencing the entire E6 and E7 genes. Phylogenetic trees were constructed by Maximum likelihood method by MEGA 5.05 software. In this study, the overall HPV infection rate was 22.6% (2891/12780) in Southeast China and the prevalence of HPV58 infection rate was 2.6% (335/12780). 26 nucleotides substitutions were observed in E6 and E7 genes with 10 novel substitutions and 17 non-synonymous substitutions. We obtained 25 distinct variation patterns which the accession GenBank numbers as MH348918-MH348942. All of HPV58 variants belong to lineage A, while no lineage B, C and D were detected in Taizhou area, Southeast China. The sublineage A1, A2, and A3 variants were found in 136 (68.3%), 39 (19.6%), and 24 (12.1%) of HPV58 isolates, respectively. The sublineage A3 variants with T20I/G63S substitutions at E7 oncoprotein carried a significantly higher risk for high-grade cervical intraepithelial neoplasia (CIN2 or worse, CIN2+) when compared with other HPV58 variants (odds ratio = 4.41, P < 0.05). Nevertheless, there was no association between HPV58 (sub) lineages and cervical lesions. These data provide the critical characteristics of HPV58 variants to assist further investigation of carcinogenic association and the development of next generation vaccines and diagnostic assays in China.

The polymorphisms of LCR, E6, and E7 of HPV-58 isolates in Yunnan, Southwest China

Backgroud

Variations in HPV LCR/E6/E7 have been shown to be associated with the viral persistence and cervical cancer development. So far, there are few reports about the polymorphisms of the HPV-58 LCR/E6/E7 sequences in Southwest China. This study aims to characterize the gene polymorphisms of the HPV-58 LCR/E6/E7 sequences in women of Southwest China, and assess the effects of variations on the immune recognition of viral E6 and E7 antigens.

Methods

Twelve LCR/E6/E7 of the HPV-58 isolates were amplified and sequenced. A neighbor-joining phylogenetic tree was constructed by MEGA 7.0, followed by the secondary structure prediction of the related proteins using PSIPRED v3.3. The selection pressure acting on the HPV-58 E6 and E7 coding regions was estimated by Bayes empirical Bayes analysis of PAML 4.8. Meanwhile, the MHC class-I and II binding peptides were predicted by the ProPred-I server and ProPred server. The transcription factor binding sites in the HPV-58 LCR were analyzed using the JASPAR database.

Results

Twenty nine SNPs (20 in the LCR, 3 in the E6, 6 in the E7) were identified at 27 nucleotide sites across the HPV-58 LCR/E6/E7. From the most variable to the least variable, the nucleotide variations were LCR > E7 > E6. The combinations of all the SNPs resulted in 11 unique sequences, which were clustered into the A lineage (7 belong to A1, 2 belong to A2, and 2 belong to A3). An insertion (TGTCAGTTTCCT) was found between the nucleotide sites 7280 and 7281 in 2 variants, and a deletion (TTTAT) was found between 7429 and 7433 in 1 variant. The most common non-synonymous substitution V77A in the E7 was observed in the sequences encoding the α-helix. 63G in the E7 was determined to be the only one positively selected site in the HPV-58 E6/E7 sequences. Six non-synonymous amino acid substitutions (including S71F and K93 N in the E6, and T20I, G41R, G63S/D, and V77A in the E7) were affecting multiple putative epitopes for both CD4⁺ and CD8⁺ T-cells. In the LCR, C7265G and C7266T were the most variable sites and were the potential binding sites for the transcription factor SOX10.

Conclusion

These results provide an insight into the intrinsic geographical relatedness and biological differences of the HPV-58 variants, and contribute to further research on the HPV-58 epidemiology, carcinogenesis, and therapeutic vaccine development.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-0986-7) contains supplementary material, which is available to authorized users.