Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 Feb 13;22(1):35.
doi: 10.1186/s12985-025-02650-9.

Genetic variations and carcinogenicity analysis of E6/E7 oncogenes in HPV31 and HPV35 in Taizhou, China

Affiliations

Genetic variations and carcinogenicity analysis of E6/E7 oncogenes in HPV31 and HPV35 in Taizhou, China

Hao-Bo Yuan et al. Virol J. .

Abstract

Background: The purpose of this study was to investigate the genetic variations in the E6 and E7 oncogenes of HPV31 and HPV35, and to explore their potential role in cervical cancer risk among women in Taizhou, China.

Methods: Cervical exfoliated cells were collected for HPV genotyping, and only patients with a single infection of either HPV31 or HPV35 were selected for this study. The ABI 3730xl sequencer was utilized to sequence the E6 and E7 genes, followed by subsequent sequence alignment, analysis of genetic heterogeneity, and construction of maximum likelihood phylogenetic trees for the sequences of HPV31 and HPV35 using BioEdit and MEGA softwares.

Results: From 2013 to 2023, 148 HPV31 E6/E7 gene sequences and 121 HPV35 E6/E7 gene sequences were successfully obtained. We identified 16 distinct HPV31 E6/E7 variants and 5 distinct HPV35 E6/E7 variants, which have been deposited in GenBank under accession numbers OR540563-OR540578 and OR540579-OR540583, respectively. Phylogenetic analysis revealed that most of the HPV31 variants belonged to sublineage A2 (57.4%), followed by sublineages C2 (26.4%), C3 (14.2%) and B1 (2.0%). The proportion of CIN2 + patients in sublineage A2 was greater than that in other HPV31 sublineages, but the difference was not statistically significant (69.2% vs. 30.8%, P > 0.05). The most common variant in A2 was 31CNTZ07, which has a greater risk of CIN2 + than other A2 variants (OR = 3.5, 95% CI = 1.31 to 9.36; P < 0.05). In addition, all the HPV35 variants belonged to lineage A, of which 99.2% belonged to sublineage A1. 35CNTZ01 and 35CNTZ03 were the two most common HPV35 variants in our population, but no significant difference in their carcinogenic ability was observed (P < 0.05).

Conclusion: These data provides a deeper insight into the distribution of geographic/ethnical HPV31 and HPV35 variants, which contribute to the development of multivalent HPV vaccines and diagnostic assays that are suitable for Chinese women.

Keywords: Cervical cancer; E6/E7 oncogenes; Genetic variability; Human papillomavirus 31; Human papillomavirus 35; Phylogenetic analysis.

PubMed Disclaimer

Conflict of interest statement

Declarations. Ethics approval and consent to participate: This study was approved by the Institutional Medical Ethics Review Board of Taizhou Hospital of Zhejiang Province (approval #K20220226), and was carried out in line with the Helsinki Declaration. All participants provided written informed consent for study participation before specimen collection, and the patients’ privacy is strictly protected. Consent for publication: Written informed consent was obtained from all patients for the publication of their medical data. Competing interests: The authors declare no competing interests.

Figures

Fig. 1
Fig. 1
Genetic variability of HPV31 E6 and E7 nucleotide sequences in Taizhou area, Southeast China. Numbering refers to the first nucleotide of the HPV31 prototype reference sequence (GenBank: J04353). Each row indicates the variant identification and the nucleotide sequence alignment compared to the reference
Fig. 2
Fig. 2
Genetic variability of HPV35 E6 and E7 nucleotide sequences in Taizhou area, Southeast China. Numbering refers to the first nucleotide of the HPV35 prototype reference sequence (GenBank: HQ537708). Each row indicates the variant identification and the nucleotide sequence alignment compared to the reference
Fig. 3
Fig. 3
Phylogenetic tree of the HPV31 E6/E7 variants. Maximum-likelihood analysis (with MEGA-X) of E6 and E7 nucleotide sequences was inferred from 16 obtained HPV31 variants and 7 reference sequences. Numbers below branches indicate bootstrap values
Fig. 4
Fig. 4
Phylogenetic tree of the HPV35 E6/E7 variants. Maximum-likelihood analysis (with MEGA-X) of E6 and E7 nucleotide sequences was inferred from 5 obtained HPV35 variants and 3 reference sequences. Numbers below branches indicate bootstrap values

Similar articles

References

    1. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Cancer J Clin. 2022;74(3):229–63. 10.3322/caac.21834 - PubMed
    1. Qiu H, Cao S, Xu R. Cancer incidence, mortality, and burden in China: a time-trend analysis and comparison with the United States and United Kingdom based on the global epidemiological data released in 2020. Cancer Commun. 2021;41(10):1037–48. 10.1002/cac2.12197 - PMC - PubMed
    1. Rahangdale L, Mungo C, O’Connor S, Chibwesha CJ, Brewer NT. Human papillomavirus vaccination and cervical cancer risk. BMJ. 2022;379:e070115. 10.1136/bmj-2022-070115 - PubMed
    1. Chen Z, de Freitas LB, Burk RD. Evolution and classification of oncogenic human papillomavirus types and variants associated with cervical cancer. Methods Mol Biol. 2015;1249:3–26. 10.1007/978-1-4939-2013-6_1 - PMC - PubMed
    1. Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013;445(1–2):232–43. 10.1016/j.virol.2013.07.018 - PMC - PubMed

Publication types

MeSH terms

Substances

Supplementary concepts

LinkOut - more resources