Association of cervical carcinogenesis risk with HPV16 E6 and E7 variants in the Taizhou area, China

Background Human papillomavirus (HPV) type 16 accounts for a larger share of cervical cancer and has been a major health problem worldwide for decades. The progression of initial infection to cervical cancer has been linked to viral sequence properties; however, the role of HPV16 variants in the risk of cervical carcinogenesis, especially with longitudinal follow-up, is not fully understood in China. Methods We aimed to investigate the genetic variability of HPV16 E6 and E7 oncogenes in isolates from cervical exfoliated cells. Between December 2012 and December 2014, a total of 310 single HPV16-positive samples were selected from women living in the Taizhou area, China. Sequences of all E6 and E7 oncogenes were analysed by PCR-sequencing assay. Detailed sequence comparison, genetic heterogeneity analyses and maximum-likelihood phylogenetic tree construction were performed with BioEdit Sequence Alignment Editor and MEGA X software. Data for cytology tests and histological diagnoses were obtained from our Taizhou Area Study with longitudinal follow-up for at least 5 years. The relationship between HPV16 variants and cervical carcinogenesis risk was analysed by the chi-square test or Fisher’s exact test. Results In this study, we obtained 64 distinct variation patterns with the accession GenBank numbers MT681266-MT681329. Phylogenetic analysis revealed that 98.3% of HPV16 variants belong to lineage A, in which the A4 (Asian) sublineage was dominant (64.8%), followed by A2 (12.1%), A1 (11.4%), and A3 (10.0%). The A4 (Asian) sublineage had a higher risk of CIN2+ than the A1–3 (European) sublineages (OR = 2.69, 95% CI = 1.04–6.97, P < 0.05). Furthermore, nucleotide variation in HPV16 E6 T178G is associated with the development of cervical cancer. Conclusion These data could provide novel insights into the role of HPV16 variants in cervical carcinogenesis risk in China. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-021-08531-y.


Background
Cervical cancer ranks as the fourth most common cancer in women worldwide, and approximately 90% of cervical cancer deaths occur in less developed countries [1]. Persistent human papillomavirus (HPV) infections are the major risk factor for cervical cancer, where nearly 100% of cervical cancer tissues contain HPV DNA [2]. Currently, more than 170 HPV types have been well identified, and approximately 40 of them can be easily transmitted in human mucosa or epithelium [3,4]. Of the known oncogenic (or high-risk [HR]) HPV types, HPV16 is the most frequently (65%) involved in cervical cancer worldwide [5,6]. Our previous epidemiological studies showed that the HPV16 type, with an infection rate of 3.1%, was less common in the Chinese general population but was found most frequently in patients with cervical intraepithelial neoplasia 2 (CIN2), CIN3 and cervical cancer (31.6, 47.9 and 60.6%, respectively) [7,8]. Many lines of evidence have indicated that the genetic variability of the HPV16 genome contributes to HPV-driven cervical carcinogenesis [9,10].

Subject recruitment and HPV genotyping
Cervical exfoliated cell specimens were collected from Chinese patients who underwent cervical cancer screening at our gynaecological clinic from December 2012 to December 2014. HPV genotyping was performed using the GP5+/bioGP6 + −PCR/MPG assay for 27 genotypes using standard procedures (CFDA Certified No. (2017): 3404697) as previously reported [8]. The inclusion criteria were single HPV16-infected patients with no pregnancy and no history of total uterus or cervix resection. Total cellular DNA was extracted using the DNeasy Blood & Tissue Kit (QIAGEN) and stored at − 20°C.
Subsequently, PCR products were purified and sequenced on the ABI 3730xl DNA Genetic Analyzer at BGI (Shanghai, China). To avoid PCR or sequencing errors, all the data were confirmed at least twice by repeating the PCR and sequencing reactions.

Treatment and follow-up
Patients with HPV16 infection were recommended for colposcopy examination immediately, according to the 2012 guidelines of cervical cancer screening [20]. Directed biopsy of suspected lesions was performed using cervical biopsy forceps. The biopsy tissues were placed in a 10% formalin solution for histological examination. Histological diagnoses were performed by pathologists and classified as normal, CIN1, CIN2, CIN3 or invasive cervical cancer, according to the WHO histological criteria.
Loop electrosurgical excision procedure (LEEP) was performed for treatment of CIN2 or worse (CIN2+). After the LEEP operation, re-examinations of HPV genotyping and cytology co-testing were performed every 3 months of follow-up, and colposcopy examination was performed for patients with abnormal followup results. Data for HPV genotyping, cytology tests and histological diagnoses were obtained from our Taizhou Area HPV Study [8] with longitudinal follow-up for at least 5 years.
Statistical analysis SPSS 16.0 statistical software (SPSS Inc., Chicago, IL) was used for this study. The association of cervical carcinogenesis risk with HPV16 variants was analysed using the chi-square test or Fisher's exact test. Odds ratios (ORs) and relative 95% confidence intervals (95% CI) were calculated. All statistical tests were two-sided. P values < 0.05 were accepted as statistically significant.

Characteristics of the study population
Between December 2012 and December 2014, single HPV16 positivity was detected in 310 women (median age 41.6 years; range 19-77) who were selected for this study. The flow diagram of the present study is shown Fig. 1. A total of 298 (96.1%) sequences of the E6 and E7 genes from HPV16 isolates were obtained. Twelve (3.9%) sequences were excluded due to the small number of HPV copies. Among 298 HPV16-infected women, 63 (21.1%) had normal and adequate cervix by colposcopy examination, so no further biopsy diagnosis was performed. 41 (13.8%) women refused further colposcopy examination. The remaining 194 (65.1%) underwent colposcopy biopsy for diagnosis which were further used for risk association analysis, and 65 were diagnosed with normal cervices after biopsy, 22 with CIN1, 31 with CIN2, 61 with CIN3, and 15 with squamous cell carcinoma (SCC). Characteristics of the study population categorized by HPV16 (sub)lineage are shown in Table 1.

Variations in the E6 and E7 genes
Compared with the prototype HPV16 reference sequence K02718, 97.7% (291/298) of the HPV16 isolates showed nucleotide variation in our study. A summary of nucleotide and amino acid sequence variation throughout the E6 and E7 fragments is shown Fig. 2. We obtained 64 distinct variation patterns denoted as 16CNTZ01-16CNTZ64, which were published with the GenBank accession codes MT681266 to MT681329. In Fig. 1 Flow diagram of the present study this study, 43 (67.2%, 43/64) novel HPV16 variants were detected, which are highlighted in bold in Fig. 2. It is worth noting that the insertion of ATAATC between nt561 and nt562 was first detected in the 16CNTZ64 variant, which accounted for 0.7% (2/298) of HPV16 isolates.
At follow-up, 56 women with histologically confirmed CIN2/3 underwent the LEEP operation, among whom, 10 women were diagnosed as the residual/ recurrent disease during the follow-up visit. Seven women with positive margin: 3 women with HPV16 positive and cytology  D32E  38  11  18  46 10 123

Discussion
China is one of the top contributors to the global burden of cervical cancer, of which HPV16 causes 70% of all cervical cancer worldwide, and some HPV16 variants are more oncogenic than others [13,21]. It has been clear that the distribution of HPV16 (sub)lineages depends on geographical origin and ethnicity, combined with more recent migration patterns (especially migration from Europe and Africa to the Americas) [22,23]. Globally, the data showed that variants in the HPV16 sublineages A1-3 were the most widespread in Europe, lineage D in Central-South America, and sublineage A4 in Asia, and lineages B and C variants were mostly restricted to Africa [24]. Several epidemiological studies have suggested that HPV16 non-European variants are associated with the persistence of HPV infection and its progression to cervical cancer, especially sublineages A4, C, D2, and D3 [11, 14,16,17,24,25].
In this study, we aimed to shed light on the HPV16 variants currently circulating in the Taizhou area, China. We obtained 298 complete sequences of the E6 and E7 genes from HPV16 isolates. Most of the HPV16 variants belong to the A lineage (98.3%), of which the A4 (Asian) sublineage was dominant (64.8%). The results of our study are in agreement with previous data in Central (60.3%) and South China (65.5%) [26][27][28][29]. In China, HPV16 variants in the A1-3 (European) sublineages are common in Xinjiang, which lies in Northwest China [30,31]. Our results indicated that the A4 (Asian) variants had an increased risk for CIN2+ compared to the A1-3 (European) variants (P < 0.05), which is in accordance with the previous findings reported in a pooled worldwide analysis (adjusted by country) [14]. Thus, the higher contribution of HPV16 to cervical cancer in China may be due to the higher oncogenicity of sublineage A4 (Asian) variants [25]. To our knowledge, the number of HPV16 isolates analysed in this study represents the largest reported collection sampled in Southeast China.
As the HPV16 oncoproteins play an important role in the development of cervical cancer, nucleotide variations in the HPV16 E6 and E7 genes could also be associated with the progression of cervical carcinogenesis. Here, we identified 54 single nucleotide substitutions in HPV16 isolates in the Taizhou population, including 37 in the E6 gene and 17 in the E7 gene (Fig. 2). In particular, we found that the three most common variations in HPV16 isolates were the T178G (D32E) in E6 (64.1%) and A647G (N29S) and T846C in E7 (65.4 and 64.4%, respectively), and these three nucleotide substitutions are apparently linked because of their simultaneous occurrence in 96.4% of A4 (Asian) variants in this study. It has been reported that the prevalence of both nonsynonymous substitutions D32E and N29S is much higher in Asia (65.5% in China, 68.0% in Korea, 44.2% in Japan and 73.9% in Thailand) than in Europe (2%), North America (3%), and Africa (0%) [25,[32][33][34][35]. Our results showed that the oncogenicity of E6 D32E and E7 N29S variations was associated with the development of cervical cancer. It has been confirmed that the nonsynonymous substitution A647G(N29S) may block the physiological function of Rb, thereby maintaining longterm infection of HPV and increasing the likelihood of persistent viral infection and cervical cancer progression [31,36,37]. Therefore, these data may help to explain the higher cervical cancer burden observed in China.
T350G (L90V) has been shown to be associated with the progression of cervical lesions [38][39][40]. Many studies have reported that 350 T/G polymorphism is a common variation located in the E6 oncogene, but it has significant heterogeneity by the world region [41]. In Europe/ Central Asia and East Asia, cervical cancer risk was significantly associated with the 350 T. However, the opposite was true in South Asia and South/Central America [40,41]. In our study, we did not find that HPV16 E6 T350G (L90V) was associated with cervical cancer risk, which is in good agreement with previous data in China [28]. Notably, 350 T/G polymorphism has been shown to influence the cervical cancer risk of European lineages, and also occurs in non-European lineages [40]. Additionally, HPV16 E6 H85Y and E120D belong to the A2 sublineage in this study but are not associated with cervical cancer. Therefore, the carcinogenicity of 350 T vs 350G, 85H vs 85 V, 120E vs 120D might be population-dependent. E6 variations might contribute to oncogenesis by disrupting p53 degradation by affecting the interaction between E6 and p53/E6AP.
Our follow-up data showed that 56 women with HPV16-positive CIN2/3 were treated by the LEEP operation, and 10 women (17.9%) were diagnosed as the residual/ recurrent disease during the follow-up visit. Among 5 women with rCIN2/3, 80% (4/5) of HPV16 isolates belong to the A4(Asian) sublineage. The reasons for the residual/ recurrent after the LEEP operation might be incomplete excision of baseline CIN2/3, highrisk HPV DNA persistence, or multifocal disease. Our follow-up data also showed that 65 women with histologically confirmed normal cervix, among whom, 6 women (9.2%) were progressed to CIN grade. 66.7%(4/6) of HPV16 isolates belong to the A4(Asian) sublineage. The persistence HPV16 DNA is the root cause of the recurrent disease for the women treated for the CIN2/3, and the same HPV16 variants were found in both lesions [42]. Therefore, our follow-up data showed that 24.7% (48/194) women with persistent infection of HPV16 genotype, which may involves a persistent infection with the same variant.
The largest limitation of our study was the time of follow-up. Since the time between HPV16 infections and CIN2/3 lesions was calculated to be 5-15 years postinfection, our longitudinal follow-up for HPV genotyping, cytology tests and histological diagnoses may affected the outcome in the present study.

Conclusions
In summary, the present study reported genetic variants in the HPV16 E6 and E7 genes in the Taizhou area, Southeast China, and provided evidence for their involvement in the increased risk of cervical carcinogenesis. Data about different HPV16 variants in the population of specific regions have significance in uncovering the carcinogenic mechanism of HPV16 and in developing preventive and therapeutic vaccines against HPV.
Additional file 1. Clinical data for HPV16 study in Taizhou area, China.

Acknowledgments
We appreciate all patients for their contribution to this study.
Authors' contributions HHX, MZD, and WWS designed the experiments, performed analysis and drafted the manuscript. YQ, XHD, and HHX carried out the sample collection, performed laboratory testing. All authors read and approved the final manuscript.

Funding
This work was supported by grants from National Natural Science Foundation of China (81901625), National Natural Science Foundation of Zhejiang province (LY20H100004), and by Science and Technology Bureau of Taizhou (1901ky05). None of the funders had any influence on the study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.

Availability of data and materials
All data generated during this study are included in this published article. The supplementary materials included the nucleotide variations of the E6 and E7 genes from HPV16 isolates and the follow-up data of patients. In addtion, these sequences have been released to GenBank database with the accession codes of MT681266 to MT681329. The links are https://www.ncbi. nlm.nih.gov/nuccore/MT681266~https://www.ncbi.nlm.nih.gov/nuccore/ MT681329.

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 # MERB-2017-020). All patients provided written informed consent before specimen collection.

Consent for publication
Written informed consent was obtained from all patients for the publication of their medical data.