MCO Molecular and Clinical Oncology 2049-9450 2049-9469 D.A. Spandidos MCO-17-4-02582 10.3892/mco.2022.2582 Articles Searching for the methylation sites involved in human papillomavirus type 16 and 18-positive women with cervical cancer MaYanyun 1 WangChunxia 1 ShiMengqi 2 LiMingshan 2 LiLin 3 CheTuanjie 3 QuJing 2 Clinical Laboratory, The First People's Hospital of Lanzhou, Lanzhou, Gansu 730050, P.R. China Department of Cell Biology, School of Biology and Basic of Medical Science, Medical College of Soochow University, Suzhou, Jiangsu 215153, P.R. China Laboratory of Precision Medicine and Translational Medicine, Suzhou Science and Technology Town Hospital, Suzhou, Jiangsu 215153, P.R. China Correspondence to: Professor Jing Qu, Department of Cell Biology, School of Biology and Basic of Medical Science, Medical College of Soochow University, 199 Renai Road, Suzhou, Jiangsu 215153, P.R. China qujing@suda.edu.cn

Professor Tuanjie Che, Laboratory of Precision Medicine and Translational Medicine, Suzhou Science and Technology Town Hospital, 1 Lijiang Road, Suzhou, Jiangsu 215153, P.R. China chetj@126.com

 10 2022 02 09 2022 17 4 149 17 03 2022 27 07 2022 Copyright: © Ma et al. 2020 This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

It has been reported that >90% of women with cervical cancer are human papillomavirus (HPV)-positive, with HPV16 and 18 being the most ‘highest-risk’ HPV genotypes. However, in numerous women, HPV infection will not progress to cervical cancer. Accordingly, more appropriate screening markers need to be explored. In the present study, genome-wide DNA methylomic differences between cervical cancer tissues with HPV-16 or HPV-18 infection and normal cervical tissues were detected by using an Illumina Human Methylation 850 K BeadChip. The Gene Ontology functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted in order to define the nearest neighbouring genes of differentiated methylation sites. Moreover, differentiated methylation sites were verified using pyrosequencing. KEGG analyses suggested that the focal adhesion pathway and pathways in cancer were highly enriched. Bioinformatics and statistical analysis indicated that the nine CpG loci had the most significant differences amongst the genes involved in these pathways. Among these, six CpG sites in the CHRM2, LAMA4, COL11A1, FGF10, IGF1 and TEK genes were highly associated with HPV-16-positive cervical cancer, as validated using pyrophosphate sequencing. Additionally, 10 significantly different CpG sites of the HPV-18-positive group were selected and verified in The Cancer Genome Atlas, indicating their possible diagnostic roles in cervical cancer development and determination. In addition, eight hypermethylated CpG island sites that were associated with HPV-16-positive cervical cancer tissues and 10 hypermethylated CpG island sites that were associated with HPV-18-positive cervical cancer tissues were identified, highlighting their potential roles in screening and evaluating targeted therapy efficacy and prognosis. The main focus of the present study was to identify the genetic variability in HPV-16- and HPV-18-positive samples and to elucidate possible methylation biomarkers in HPV-positive women with a risk of developing cervical cancer.

 cervical cancer human papillomavirus 16 human papillomavirus 18 methylation biomarker Funding: The present study was supported by the innovation and entrepreneurship talent project of Lanzhou (2015-RC-68), distinguished professorship program of national research program on prevention and control of major birth defects in reproductive health (2017YFC1000900), special-funded program on national key scientific instruments and equipment development (2016YFF0103800), health and family planning commission research project of Jiangsu province (H201619). Introduction

Cervical cancer is considered one of the leading causes of cancer-associated mortality among women globally (1,2). In particular, >85% of new cases and 90% of cervical cancer-related deaths occur in developing countries (3-5). Despite efforts made to improve cervical cancer therapy, the 5-year survival rate remains <50% (6,7). In China, the incidence rate of cervical cancer is estimated to be ~15.4/100,000, with a relatively high mortality rate (8-10). Due to the huge population and high rate of human papillomavirus (HPV) infection, precise diagnosis and treatment are required for cervical cancer in China.

HPV, a highly prevalent sexually transmitted virus, is a circular dsDNA virus containing six early genes (E1, E2, E4, E5, E6 and E7) and two late genes named L1 and L2 (11-16). Currently, ~200 HPV genotypes have been identified based on the nucleotide diversity of the L1 gene, and 15 of these are regarded as ‘high risk’, contributing to the development of cervical cancer, including HPV-16, -18, -31, -33, -35, -39, -45, -51, -52, -53, -56, -58, -59, -66 and -68 (17-21). It has been reported that HPV16 is the ‘highest-risk’ HPV genotype with HPV-18 being second, according to their oncogenic potential (22-27).

More than 1 year of persistent HPV infection may be an important risk factor for the progression of cervical cancer and its precursors (28,29). HPVs can establish their progeny, spread their viral genes, infect basal cells and further promote epithelial-mesenchymal transition (30,31). However, HPV infection itself is not sufficient for the initiation and establishment of malignant cell transformation (32,33). Accordingly, HPV tests have a poor positive predictive value, since HPV infection will progress to cervical cancer in only a few women (34,35). Clinical data have indicated that HPV infection is self-limiting and regresses in several cases, suggesting that other biomolecular mechanisms are involved in the progression of cervical cancer. A number of studies have demonstrated that DNA methylation is involved in the carcinogenic process of cervical cancer (36-40). The hypermethylation of CpG islands in the promoter regions of specific genes, including tumour suppressor genes, leads to the silencing of the gene and inhibits the downstream pathways. By contrast, disruption of epigenetic processes can lead to the activation of oncogenes, and the accumulation of epigenetic changes is an essential step in the development of cervical cancer (41,42).

An increasing number of studies have indicated that DNA methylation is an early event in tumorigenesis and plays a major role in tumour initiation and the progression of cervical cancer (43). Therefore, it is crucial to identify reliable prognostic and predictive DNA methylation-related biomarkers that may help in the early diagnosis and treatment strategies for cervical cancer. In the present study, to elucidate the effect of the combination of HPV genotypes and DNA methylation, these methylation biomarkers in HPV-16- and HPV-18-positive women with cervical carcinoma were analysed using a human Illumina Human Methylation 850 K BeadChip. The purpose of the present study was to identify the hypermethylation of CpG islands of genetic variability from samples tested positive HPV-16 and HPV-18 and to elucidate possible methylation biomarkers of HPV-positive women with a risk of developing cervical cancer.

 Materials and methods <sec> <title>Human tissue specimens

A total of six paraffin-embedded specimens, including three cases of HPV 18 (HPV-18 group)- and 3 cases of HPV-16 (HPV-16 group)-positive cervical carcinoma tissues, were collected and diagnosed on a pathological basis according to the FIGO (2009) clinical staging criteria (44). Normal cervical tissues were obtained from three women with hysteromyoma who underwent total hysterectomy (normal group) from July, 2014 through December, 2017 at the Department of Obstetrics and Gynecology of the First People's Hospital of Lanzhou and Gansu Provincial Hospital. All experiments performed in the present study were approved by the Ethics Committee of The First People's Hospital of Lanzhou and Gansu Provincial Hospital (approval no. (2016-02). Written informed consent was obtained from all the patients, and the time of signing the agreement was the time of sample collection.

 HPV genotyping and grouping

DNA extraction was performed using a QIAGEN QIAamp DNA Mini kit (cat. no. 51304; Qiagen GmbH) from formalin-fixed and paraffin-embedded tissue sections according to the manufacturer's instructions. The DNA samples were identified and quantified using a NanoDrop™ 8000 spectrophotometer (Thermo Fisher Scientific, Inc.) and agarose gel electrophoresis. Genotyping with positive samples was performed by using the HPV Genotyping Detection kit, the Assay Kit for Genotyping Human Papillomavirus (PCR-reverse dot blot) (cat. no. CP.008.022; Guangzhou, LBP Medicine Science & Technology Co., Ltd.) for subtypes HPV-16 and HPV-18. Samples were tested following the manufacturer's instructions (45). The HPV-positive control and -negative controls were set in each experiment. In total, 12 samples were genotyped with the HPV Genotyping Detection kit, and three HPV-16-positive specimens and three HPV-18-positive specimens were randomly selected (Table SI).

 DNA methylation chip

An Illumina Human Methylation 850K BeadChip (Illumina, Inc.) was used to detect the whole genome methylation status of HPV-16- and HPV-18-positive tissues. Genomic DNA of the normal group (three cases normal cervical tissues), HPV-16 group (HPV-16-positive specimens with cervical cancer), and HPV-18 group (three cases of HPV-18-positive specimens with cervical cancer) was extracted using the QIAamp DNA Mini kit (cat. no. 51304; Qiagen GmbH) and bisulfite-converted using the EZ DNA Methylation kit (cat. no. D5001; Zymo Research Corp.). The converted DNA was hybridized to an Infinium Human Methylation 850 K BeadChip. The subsequent bioinformatics analysis was performed by Genergy Co. The Illumina 850 K methylation chip analysis data have been uploaded in the GEO public database repository (accession no. GSE169622, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE169622).

 Pyrosequencing

Through pyrosequencing, nine candidates of the HPV-16 group screened by the 850K methylation chip were verified. An EZ 96-DNA methylation kit (Zymo Research Corp.) was used for bisulfite conversion in accordance with the manufacturer's standard procedures, with fully methylated and unmethylated samples as test controls. PyroMark Assay Design 2.0 was used for the synthesis of bisulfite-PCR primers, which were synthesized by The Beijing Genomics Institute (BGI). A list of bisulfite PCR primers is presented in Table SII. The bisulfite PCR amplification conditions were as follows: Pre-denaturation at 95˚C for 3 min; 40 cycles at 94˚C for 30 sec, 52˚C for 30 sec, and 72˚C for 1 min; and a final elongation at 72˚C for 7 min. The HotStarTaq DNA polymerase (Qiagen Ltd.) was used for disulfide PCR amplification. Compared to the PyroMark Q96 ID (Qiagen Ltd.), Pyro Q CpG software (version 2.0.6, Qiagen GmbH) automatically analysed the methylation status of each site.

 Bioinformatics analysis

The 850K chip data analysis w implemented in R language. The whole analysis model is based on the highly integrated R analysis package ChAMP (Version: 2.8.9), which inherits the methods of Minfi, Limma, Sva, and IMA analysis packages. The graph was acquired using self-written R script, the basic function in GGplot2 implementation. Microarray data were normalized using BMIQ (Beta MIXture Quantile dilation). SVD (Singular Value Decomposition) was applied to evaluate the major components of variables in the data set, and then a Bayesian model-based Combat method was used to eliminate the batch effect. Quality control is achieved through a set of functions provided by ChAMP, such as CpG.GUI, champ.QC, and QC.GUI, and then MVP (Methylation variable position), DMR (Differently methylated regions), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis are also implemented with ChAMP. Among them, MVP uses the method of Limma R package. Limma firstly establishes multiple linear regression model for all data as a whole, and then applies the regression model to each probe line and calculates the modified T-statistic and P-value in combination with the Bayesian model. The sites with P≤0.05 after the FDR correction were considered as differential methylation sites.

 Cervical cancer data sets

The DNA methylation data from The Cancer Genome Atlas (TCGA)-Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (CESC) were downloaded from TCGA website (https://portal.gdc.cancer.gov/). β-values were extracted to evaluate the DNA methylation level of each probe. The candidate significantly differentially methylated CpG sites of the HPV-18 group were screened using the 850 k methylation chip were verified in (TCGA-CESC).

 Statistical analysis

Statistical analysis was performed using SPSS software (Release 13.0, SPSS Inc.) Statistical analysis data comparisons between two groups were analysed using the Student's t-test. P<0.05 was considered to indicate a statistically significant difference.

 Results <sec> <title>Methylation analysis and general characteristics

The 850K methylation sites in the HPV 16- and HPV 18-positive cervical carcinoma and normal cervical tissues were analysed using the Illumina 850K methylation chip (GSE169622). The data from cervical carcinoma and normal cervical tissue samples were normalized (negative and positive controls were provided by Genergy Co.) and processed. To analyse DNA methylation differences between the cervical cancer group and the normal group, Δ β and P-values were used to construct a volcanic map of CpG sites, in order to reflect the magnitude and statistical significance of differences (P<0.01). As depicted in Fig. 1, the gene methylation profiles of HPV-16 (HPV-16 group) and HPV-18 (HPV-18 group) cervical carcinoma samples differed significantly from those of the control samples (normal group). In total, it was observed that 106,378 sites demonstrated differential expression in HPV 16-positive cervical carcinoma tissues compared to normal cervical tissues, of which 101,152 were hypermethylated and 5,226 were hypomethylated (Fig. 1A-a). In addition, 70,744 sites showed differential expression in HPV 18-positive cervical carcinoma tissues compared to normal cervical tissues, of which 53,168 were hypermethylated and 17,576 were hypomethylated (Fig. 1B-a). In addition, a cluster analysis map was generated to show the methylation status of different groups. As demonstrated in Fig. 1, there were significant differences in methylation patterns and states between cervical cancer and normal samples in the HPV-16 group and HPV-18 group. Thus, there may be biomarkers suitable for cervical cancer screening among these differential methylation sites.

 GO functional analysis

The GO functional annotation analysis of the HPV-16 differentially expressed methylation sites revealed that these enriched genes were mainly involved in Biological Process, Cellular Component and Molecular Function (Fig. 2A). The results of cellular component analyses revealed that molecules distributed in the cell periphery, plasma membrane, cell junction and cell membrane components were significantly enriched. Important functions, such as cell migration, transport and synthesis of substances all occur place at these sites. At the molecular level, functional annotation analysis revealed that the highly enriched genes were related to calcium binding, protein binding, cytoskeletal protein binding, metal ion transmembrane transport activity and phosphotransferase activity (Fig. 2A). The GO functional annotation analysis of the HPV-18-differentially expressed methylation sites demonstrated that these genes were mainly enriched in Biological Processes, Cellular Component and Molecular Function (Fig. 2B). The biological processes of the two groups of methylation differential genes were mainly concentrated in the biological development process and anatomical structure development.

 KEGG signalling pathway analysis

KEGG pathway functional analysis annotates and classifies the functions of pathways in the KEGG database according to whole genes and differential genes. The signalling pathways were further investigated using the KEGG database. According to the criteria of P<0.01 and FDR <0.05, the top 20 related signalling pathways of different methylation sites were selected (Fig. 3). The most prominent major signalling pathways for the HPV-16-positive samples were focal adhesion, pathways in cancer, glutamatergic synapse and the regulation of actin cytoskeleton, suggesting that these pathways are major regulatory factors of cancer behaviours (Fig. 3A). Additionally, focal adhesion, pathways in cancer, glutamatergic synapse and circadian entrainment signalling pathways were significantly upregulated in the HPV-18-positive samples (Fig. 3B). It was observed that focal adhesion and pathways in cancer are among the top pathways in the comparison of the two groups. This suggests that it may be possible to monitor the degree of cervical lesions by detecting gene differential methylation sites in these pathways.

 Identification of different methylation sites in the HPV-16 and HPV-18 groups

Subsequently, the 3,000 methylation variable positions in genes with the most significant difference were further analysed by KEGG according to the selected significantly different signalling pathways, including focal adhesion and pathways in cancer from HPV-16- and HPV-18-positive samples. The results indicated that a total of 10 genes, including CHRM2, GNG4, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1, TEK and COL11A2, were screened from the focal adhesion and PI3K-AKT signalling pathways from the HPV-16 group. Moreover, the Gene network analysis revealed that the most significantly different CpG sites, cg24575234, cg20818778, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548, cg08264401 and cg25459558, can be found in the CHRM2, GNG4, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1, TEK and COL11A2 genes, indicating their possible diagnostic roles in cervical cancer development (Table I). Some of these genes have already exhibited critical roles in several types of cancer, such as thymic, gastric, ovarian, breast, pancreatic and colorectal cancer. For examples, it has been shown that DNA methylation of GNG4 is a common epigenetic alteration in thymic carcinoma (46). LAMA4 and COL11A1 are associated with tumour invasion and metastasis (47,48).

For the HPV-18 group, the 10 most significantly different CpG sites, cg03520644, cg25792518, cg06958829, cg00172849, cg19707040, cg02501779, cg19679123, cg14427009, cg25993718 and cg27423357 was selected. The gene network analysis indicated that they were located in COL11A1, CHAD, CHAD, COL11A1, CTNNA2, CBLN4, SMAD3, PCDH17, CBLN4 and FLT1 (Table II). In addition, the 10 selected methylation significantly different sites were verified in TCGA-CESC (Table III). All 10 sites exhibited statistically significant differences (P<0.05), and seven sites demonstrated highly significant differences (P<0.01, Table III).

 Pyrosequencing verification

In total, 10 candidate significantly differentially methylated CpG sites of the HPV-16 group, including cg24575234, cg20818778, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548, cg08264401 and cg25459558, screened using the 850k methylation chip were verified by pyrosequencing. Among these, cg25459558 was withdrawn from verification due to the failure of primer design. The statistical comparison of the mean value revealed that the other nine sites exhibited significant differences (P<0.05), of which six sites exhibited highly significant differences (P<0.01, Table IV). The statistical comparison of means revealed that the most significantly different CpG sites between the normal and HPV-16 samples cg24575234, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548 and cg08264401, detected in the CHRM2, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1 and TEK. Unexpectedly, the increasing mean value of the CpG site, cg20818778, for the gene GNG4 was not statistically significant (detailed statistical analysis data are contained in Table IV).

The mechanism of cervical cancer development remains unclear, while individual differences are significant. Further investigations of the methylation characteristics of a single gene and a single locus as a biomarker for cancer screening will have a high missed diagnosis rate. Multipoint joint detection is suggested for the future improvement of cancer detection rate.

 Discussion

Cervical cancer formation is affected by several risk factors, including HPV infection. It has been reported that HPV16 and HPV18 contribute to >70% of all cervical cancer cases worldwide and are thus entitled as a ‘high-risk’ HPV genotype (49-51). In some studies, HPV-16/18 genotyping is used as a molecular marker reflecting the underlying carcinogenic process. However, HPV infection is self-limiting and regresses in some clinical cases (30,52,53). It is unclear whether HPV infection acts as a key determinant of the progression to cervical cancer. Consequently, HPV positivity is not a specific diagnostic indicator for cervical cancer or diseases.

DNA methylation is one of the mechanisms that has been closely related to the occurrence and development of cervical cancer. The aberrant DNA methylation of human host cell genes or HPV genomic DNA has been closely associated with the dysfunction of various tumour suppressor genes during persistent high-risk HPV (HR-HPV) infection and cervical carcinogenesis (54-56). Studies have indicated that the tumour suppressor genes, p53 and p73, demonstrate a higher degree of methylation in cervical cancer samples than in normal samples (57,58). The aberrant DNA methylation of CpG islands is comparatively rare in normal cells, suggesting that the differentially methylated CpG sites between cervical cancer and normal samples have the potential to become reliable biomarkers of cervical cancer. In addition, DNA hypermethylation has been associated with long-term HR-HPV infection and is therefore considered a marker of cervical intraepithelial neoplasia lesion severity and invasive cervical cancer risk (59). However, the high heterogeneity of those previously published data renders it difficult to determine the appropriate methylation markers for cervical cancer screening (60). Additionally, the expression levels of E6 oncoprotein have a different effect on the carcinogenic potential of HPV. For example, the enhanced expression of the HPV-16 E6 oncogenes may trigger a neoplastic transformation of squamous epithelial cells at the uterine cervix (61). Thus, the HPV promoter methylation profile could be an easy and measurable biomarker for the examination of the high-risk HPV potential carcinogenicity.

In the present study, the genome-wide methylation level was evaluated by comparing HPV-16-positive or HPV-18-positive cervical cancer cases with normal cervical tissues. The results of the present study indicated that 106,378 and 70,744 sites demonstrated differential expression in HPV-16 and HPV-18 cervical cancer tissues as compared with normal cervical tissues, respectively, indicating that the distribution of methylation sites in cervical tissues varies greatly. It has been reported that hypermethylation at CpG islands (CGIs) of genes acting as tumour suppressors is a common mechanism involved in cancer occurrence (62-64). Other studies have also detected an apparently positive association between the hypomethylation of proto-oncogenes and the progression of cervical cancer. In the present study, 101,152 with higher methylation levels and 5,226 with lower methylation levels CGIs in HPV-16-positive cancer tissues than in normal cervical tissues were identified. By contrast, 53,168 CGIs with increased methylation levels and 17,576 CGIs with decreased methylation levels were identified in HPV-18-positive cancer tissues compared with normal cervical tissues. Genome-wide methylation level evaluation can retrieve additional differential methylation sites that have not been previously discovered.

Moreover, the differentially expressed methylation genes were analysed through GO functional annotation. It has been revealed that a number of methylated genes are closely associated with HPV-positive cervical cancer cases, including SOX1, PAX1, JAM3, EPB41L3, CADM1 and MAL (65,66). For example, expression of SOX1 was shown to be associated with early embryogenesis, central nervous system development, and neural stem cell maintenance. Hypermethylated PAX1 has been detected in cervical carcinoma (67). The aforementioned methylated human gene biomarkers used in combination may be clinically useful for the triage of women with HR-HPV infections. The functional annotation data have previously demonstrated that the highly enriched genes were mainly involved in calcium binding, protein binding, cytoskeletal protein binding, metal ion transmembrane transport activity and phosphotransferase activity (68-71). The results of cellular component analyses revealed that molecules distributed in the cell periphery, plasma membrane, cell junction, and cell membrane components were significantly enriched. Further cluster analysis demonstrated that the differentially methylated genes covered a variety of different functional communities, indicating that there are many types of genes involved in the regulation of the occurrence and progression of cervical cancer (72-74).

Previous data indicate that a variety of cellular pathways can be affected by the methylation status of specific genes. KEGG pathway analysis in the present study revealed that differentially methylated genes were mainly involved in focal adhesion, regulation of actin cytoskeleton, and pathways in cancer. Among these pathways, the most significant pathways were focal adhesion and PI3K-AKT signalling pathways, which are a collection of receptors and ligands on the plasma membrane associated with intracellular and extracellular signalling pathways that regulate cell growth and cell migration. Based on the KEGG pathway analysis results, a total of nine genes, including CHRM2, GNG4, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1 and TEK, associated with nine significantly different CpG sites, cg24575234, cg20818778, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548 and cg08264401, were screened from the focal adhesion and PI3K-AKT signalling pathways of the HPV-16-positive group. However, the pyrosequencing data of the present study indicated that the increasing mean value of the CpG site, cg20818778, for the gene GNG4 was not statistically significant. Thus, the most significantly different CpG sites are cg24575234, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548 and cg08264401, detected in the CHRM2, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1 and TEK, indicating their possible diagnostic roles in cervical carcinoma development.

Additionally, the 10 most significantly different CpG sites of the HPV-18-positive group, cg03520644, cg25792518, cg06958829, cg00172849, cg19707040, cg02501779, cg19679123, cg14427009, cg25993718 and cg27423357 in COL11A1, CHAD, CHAD, COL11A1, CTNNA2, CBLN4, SMAD3, PCDH17, CBLN4 and FLT1 were selected, indicating their future applications as candidate molecular markers of cervical cancer. Furthermore, the 10 selected methylation significantly different sites were verified in TCGA-CESC. All 10 sites exhibited significant differences (P<0.05), and seven sites demonstrated highly significant differences (P<0.01, Table III).

Some of the genes screened in the present study have already demonstrated potential functions in other diseases. Several studies have revealed that CHRM2 is associated with the pathophysiology of schizophrenia (75). GNG4 has been associated with immune infiltration in the tumour microenvironment, which promotes tumour cell migration and proliferation (76-78). LAMA4 can regulate the proliferation and migration of gastric cancer cells (47) and may be a potential gastric cancer prognostic biomarker and therapeutic target (79). ITGB8 silencing inhibits the invasion and migration of lung cancer cells (80). Increased expression of COL11A1 has been detected in several in various types of cancer, such as ovarian, breast, pancreatic and colorectal cancer, and increased levels of COL11A1 are often associated with poor survival, chemoresistance and recurrence (48). Fgf10 induces migration and invasion of pancreatic cancer cells (81). CHAD, however, has been reported to be involved in confronting hepatocellular carcinoma migration and proliferation and predicting good survival (82). In addition, PCDH17 methylation has been reported as a potential prognostic biomarker in some cancer patient markers, including postoperative renal cell carcinoma (83). Validating whether these markers can be used as novel tools for cervical cancer screening and investigating their role in normal cervical and cervical cancer development will be the focus of our future research.

In summary, the host cell gene methylation test may be a promising method for cervical cancer screening. the Illumina Human Methylation 850K BeadChip methylation chip was used for the methylation site detection in HPV 16- and HPV 18-positive cervical carcinoma and normal cervical tissues. The current findings suggested that the methylation modification sites in cervical carcinoma cells may be abnormal. Hypermethylation and hypomethylation sites occur more frequently and are mainly enriched in functional categories, including focal adhesion, regulation of actin cytoskeleton, tumour growth, and pathways in cancer. The eight most significantly different CpG sites, cg24575234, cg14289461, cg06818777, cg08361126, cg06381931, cg08976810, cg20881548 and cg08264401, were screened and verified from HPV-16-positive samples and were associated with the CHRM2, LAMA4, CHAD, ITGA8, COL11A1, FGF10, IGF1 and TEK genes. The 10 most significantly different CpG sites of the HPV-18-positive group, cg03520644, cg25792518, cg06958829, cg00172849, cg19707040, cg02501779, cg19679123, cg14427009, cg25993718 and cg27423357, which are located in COL11A1, CHAD, CHAD, COL11A1, CTNNA2, CBLN4, SMAD3, PCDH17, CBLN4 and FLT1, were selected and verified in TCGA-CESC. It is important to explore and develop DNA methylation assays of improved sensitivity and specificity in order to ameliorate the early detection of cervical cancer (84-86). The findings of the present study may provide fundamental data for the use of methylation biomarkers for cervical cancer diagnosis; however, further research is required.

 Supplementary Material HPV Genotyping and grouping. List of bisulfite PCR and sequencing primers. Acknowledgements

Not applicable.

 Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The Illumina 850K methylation chip analysis data have been uploaded in the GEO public database repository (accession number: GSE169622, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE169622).

 Authors' contributions

JQ and TC conceived and designed this study. YM and CW were responsible for performing the experiments and writing the first draft of the manuscript. JQ was responsible for revising the first draft of the manuscript. MS and ML collected and analysed the data. LL and YM interpreted the data of DNA methylation. JQ and TC confirm the authenticity of all the raw data. All the authors critically reviewed the original manuscript, edited and approved the final version. All the authors have read and approved the final manuscript.

 Ethics approval and consent to participate

The present study was conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all patients. The study was approved by the Medical Ethics Examination of Lanzhou First People's Hospital, China (2016-02).

 Patient consent for publication

Not applicable.

 Competing interests

The authors declare that they have no competing interests.

 References XingBGuoJShengYWuGZhaoYHuman papillomavirus-negative cervical cancer: A comprehensive reviewFront Oncol1060633520213368092810.3389/fonc.2020.606335 BuskwofieADavid-WestGClareCAA review of cervical cancer: Incidence and disparitiesJ Natl Med Assoc11222923220203227847810.1016/j.jnma.2020.03.002 BrayFFerlayJSoerjomataramISiegelRLTorreLAJemalAGlobal cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countriesCA Cancer J Clin6839442420183020759310.3322/caac.21492 SossoSMTchouaketMCTFokamJSimoRKTorimiroJTigaALobeEEAmbadaGNangeASemengueENJHuman immunodeficiency virus is a driven factor of human papilloma virus among women: Evidence from a cross-sectional analysis in Yaoundé, CameroonVirol J176920203243003410.1186/s12985-020-01340-y MapangaWSinghEFeresuSAGirdler-BrownBTreatment of pre- and confirmed cervical cancer in HIV-seropositive women from developing countries: A systematic reviewSyst Rev97920203227667210.1186/s13643-020-01345-2 CeglaPBurchardtERoszakACzepczynskiRKubiakACholewinskiWInfluence of biological parameters assessed in [18F]FDG PET/CT on overall survival in cervical cancer patientsClin Nucl Med4486086320193134808110.1097/RLU.0000000000002733 DingFNGaoBHWuXGongCWWangWQZhangSMmiR-122-5p modulates the radiosensitivity of cervical cancer cells by regulating cell division cycle 25A (CDC25A)FEBS Open Bio91869187920193150510510.1002/2211-5463.12730 SongBDingCChenWSunHZhangMChenWIncidence and mortality of cervical cancer in China, 2013Chin J Cancer Res2947147620172935396910.21147/j.issn.1000-9604.2017.06.01 GuXYZhengRSSunKXZhangSWZengHMZouXNChenWQHeJIncidence and mort ality of cervical cancer in China, 2014Zhonghua Zhong Liu Za Zhi4024124620182973090810.3760/cma.j.issn.0253-3766.2018.04.001(In Chinese) HeRZhuBLiuJZhangNZhangWHMaoYWomen's cancers in China: A spatio-temporal epidemiology analysisBMC Womens Health2111620213374364810.1186/s12905-021-01260-1 JeeBYadavRPankajSShahiSKImmunology of HPV-mediated cervical cancer: Current understandingInt Rev Immunol4035937820203285304910.1080/08830185.2020.1811859 BriantiPDe FlammineisEMercuriSRReview of HPV-related diseases and cancersNew Microbiol408085201728368072 TulayPSerakinciNThe route to HPV-associated neoplastic transformation: A review of the literatureCrit Rev Eukaryot Gene Expr26273920162727888310.1615/CritRevEukaryotGeneExpr.v26.i1.40 ChenJSignaling pathways in HPV-associated cancers and therapeutic implicationsRev Med Virol25 (Suppl 1)S24S5320152575281510.1002/rmv.1823 LototskajaESahharovOPiirsooMKalaMUstavMPiirsooACyclic AMP-dependent protein kinase exhibits antagonistic effects on the replication efficiency of different HPV typesJ Virol10e002512120213385396310.1128/JVI.00251-21 BhattKHNellerMASrihariSCrooksPLekieffreLAftabBTLiuHSmithCKennyLPorcedduSKhannaRProfiling HPV-16-specific T cell responses reveals broad antigen reactivities in oropharyngeal cancer patientsJ Exp Med217e2020038920203271651810.1084/jem.20200389 BouvardVBaanRStraifKGrosseYSecretanBGhissassiFEBenbrahim-TallaaLGuhaNFreemanCGalichetLA review of human carcinogens-Part B: Biological agentsLancet Oncol1032132220091935069810.1016/s1470-2045(09)70096-8 TaoGYalingGZhanGPuLMiaoHHuman papillomavirus genotype distribution among HPV-positive women in Sichuan province, Southwest ChinaArch Virol163657220182898374410.1007/s00705-017-3556-1 CastellsaguéXNatural history and epidemiology of HPV infection and cervical cancerGynecol Oncol110 (3 Suppl 2)S4S720081876071110.1016/j.ygyno.2008.07.045 SchellenbacherCRodenRBSKirnbauerRDevelopments in L2-based human papillomavirus (HPV) vaccinesVirus Res23116617520172788961610.1016/j.virusres.2016.11.020 LippertJBonlokkeSUtkeAKnudsenBRSorensenBSSteinicheTStougaardMTargeted next generation sequencing panel for HPV genotyping in cervical cancerExp Mol Pathol11810456820213317115510.1016/j.yexmp.2020.104568 BordigoniAMotteATissot-DupontHColsonPDesnuesCDevelopment and validation of a multiplex qPCR assay for detection and relative quantification of HPV16 and HPV18 E6 and E7 oncogenesSci Rep11403920213359759210.1038/s41598-021-83489-2 FanZFengXZhangWLiNZhangXLinJMVisual detection of high-risk HPV16 and HPV18 based on loop-mediated isothermal amplificationTalanta21712101520203249891110.1016/j.talanta.2020.121015 LaymanHRickertKWWilsonSAksyukAADuntyJMNatrakulDSwaminathanNDelNagroCJDevelopment and validation of a multiplex immunoassay for the simultaneous quantification of type-specific IgG antibodies to E6/E7 oncoproteins of HPV16 and HPV18PLoS One15e022967220203221436210.1371/journal.pone.0229672 PengSFerrallLGaillardSWangCChiWYHuangCHRodenRBSWuTCChangYNHungCFDevelopment of DNA vaccine targeting E6 and E7 proteins of human papillomavirus 16 (HPV16) and HPV18 for immunotherapy in combination with recombinant vaccinia boost and PD-1 antibodymBio12e03224e0322020213346869810.1128/mBio.03224-20 BertiFCBMathiasCGarciaLEGradiaDFde Araujo SouzaPSCipollaGAde OliveiraJCMalheirosDComprehensive analysis of ceRNA networks in HPV16- and HPV18-mediated cervical cancers reveals XIST as a pivotal competing endogenous RNABiochim Biophys Acta Mol Basis Dis186716617220213404892410.1016/j.bbadis.2021.166172 HammerARositchAQeadanFGravittPEBlaakaerJAge-specific prevalence of HPV16/18 genotypes in cervical cancer: A systematic review and meta-analysisInt J Cancer1382795280320162666188910.1002/ijc.29959 LuveroDLopezSBoganiGRaspagliesiFAngioliRFrom the infection to the immunotherapy in cervical cancer: Can we stop the natural course of the disease?Vaccines (Basel)859720203305048410.3390/vaccines8040597 LinWNiuZZhangHKongYWangZYangXYuanFImbalance of Th1/Th2 and Th17/Treg during the development of uterine cervical cancerInt J Clin Exp Pathol1236043612201931934210 HoGYBiermanRBeardsleyLChangCJBurkRDNatural history of cervicovaginal papillomavirus infection in young womenN Engl J Med3384234281998945964510.1056/NEJM199802123380703 BodilyJLaiminsLAPersistence of human papillomavirus infection: Keys to malignant progressionTrends Microbiol19333920112105076510.1016/j.tim.2010.10.002 GrahamSVThe human papillomavirus replication cycle, and its links to cancer progression: A comprehensive reviewClin Sci (Lond)1312201222120172879807310.1042/CS20160786 HebnerCMLaiminsLAHuman papillomaviruses: Basic mechanisms of pathogenesis and oncogenicityRev Med Virol16839720061628720410.1002/rmv.488 YooSHOckCYKeamBParkSJKimTMKimJHJeonYKChungEJKwonSKHahJHPoor prognostic factors in human papillomavirus-positive head and neck cancer: Who might not be candidates for de-escalation treatment?Korean J Intern Med341313132320193042864610.3904/kjim.2017.397 LindeDSAndersenMSMwaiselageJDManongiRKjaerSKRaschVText messages to increase attendance to follow-up cervical cancer screening appointments among HPV-positive Tanzanian women (Connected2Care): Study protocol for a randomised controlled trialTrials1855520172916214810.1186/s13063-017-2215-x LiXWuXLiYCuiYTianRSinghNDingMYangYGaoYPromoter hypermethylation of SOX11 promotes the progression of cervical cancer in vitro and in vivoOncol Rep412351236020193072013310.3892/or.2019.6993 ChenLQiuXZhangNWangYWangMLiDWangLDuYAPOBEC-mediated genomic alterations link immunity and viral infection during human papillomavirus-driven cervical carcinogenesisBiosci Trends1138338820172871706110.5582/bst.2017.01103 LorinczATCancer diagnostic classifiers based on quantitative DNA methylationExpert Rev Mol Diagn1429330520142464981810.1586/14737159.2014.897610 YangSWuYWangSXuPDengYWangMLiuKTianTZhuYLiNHPV-related methylation-based reclassification and risk stratification of cervical cancerMol Oncol142124214120203240839610.1002/1878-0261.12709 FerteyJHagmannJRuscheweyhHJMunkCKjaerSHusonDHaedicke-JarbouiJStubenrauchFIftnerTMethylation of CpG 5962 in L1 of the human papillomavirus 16 genome as a potential predictive marker for viral persistence: A prospective large cohort study using cervical swab samplesCancer Med91058106820203185641110.1002/cam4.2771 ZhangLHuDWangSZhangYPangLTaoLJiaWAssociation between dense PAX1 promoter methylation and HPV16 infection in cervical squamous epithelial neoplasms of Xin Jiang Uyghur and Han womenGene72314414220203158995710.1016/j.gene.2019.144142 FranzenAVogtTJMullerTDietrichJSchrockAGolletzCBrossartPBootzFLandsbergJKristiansenGDietrichDPD-L1 (CD274) and PD-L2 (PDCD1LG2) promoter methylation is associated with HPV infection and transcriptional repression in head and neck squamous cell carcinomasOncotarget964165020182941664110.18632/oncotarget.23080 El AlianiAEl-AbidHEl MallaliYAttalebMEnnajiMMEl MzibriMAssociation between gene promoter methylation and cervical cancer development: Global distribution and a meta-analysisCancer Epidemiol Biomarkers Prev3045045920213344130810.1158/1055-9965.EPI-20-0833 PecorelliSRevised FIGO staging for carcinoma of the vulva, cervix, and endometriumInt J Gynaecol Obstet10510310420091936768910.1016/j.ijgo.2009.02.012 MuñozNBoschFXde SanjoséSHerreroRCastellsaguéXShahKVSnijdersPJMeijerCJEpidemiologic classification of human papillomavirus types associated with cervical cancerN Engl J Med34851852720031257125910.1056/NEJMoa021641 KishibuchiRKondoKSoejimaSTsuboiMKajiuraKKawakamiYKawakitaNSawadaTTobaHYoshidaMDNA methylation of GHSR, GNG4, HOXD9 and SALL3 is a common epigenetic alteration in thymic carcinomaInt J Oncol5631532620203174637010.3892/ijo.2019.4915 WangMLiCLiuYWangZEffect of LAMA4 on prognosis and its correlation with immune infiltration in gastric cancerBiomed Res Int2021642887320213441423810.1155/2021/6428873 NallanthighalSHeisermanJPCheonDJCollagen type XI alpha 1 (COL11A1): A novel biomarker and a key player in cancerCancers (Basel)1393520213366809710.3390/cancers13050935 XuHHWangKFengXJDongSSLinAZhengLZYanWHPrevalence of human papillomavirus genotypes and relative risk of cervical cancer in China: A systematic review and meta-analysisOncotarget9153861539720182963265210.18632/oncotarget.24169 LinCFranceschiSCliffordGMHuman papillomavirus types from infection to cancer in the anus, according to sex and HIV status: A systematic review and meta-analysisLancet Infect Dis1819820620182915810210.1016/S1473-3099(17)30653-9 MalaryMMoosazadehMHamzehgardeshiZAfshariMMoghaddasifarIAfsharimoghaddamAThe prevalence of cervical human papillomavirus infection and the most at-risk genotypes among iranian healthy women: A systematic review and meta-analysisInt J Prev Med77020162721793610.4103/2008-7802.181756 ChanPKPicconiMACheungTHGiovannelliLParkJSLaboratory and clinical aspects of human papillomavirus testingCrit Rev Clin Lab Sci4911713620122291340510.3109/10408363.2012.707174 KaliternaVBarisicZGenital human papillomavirus infectionsFront Biosci (Landmark Ed)231587161120182929345210.2741/4662 ClarkeMAWentzensenNMirabelloLGhoshAWacholderSHarariALorinczASchiffmanMBurkRDHuman papillomavirus DNA methylation as a potential biomarker for cervical cancerCancer Epidemiol Biomarkers Prev212125213720122303517810.1158/1055-9965.EPI-12-0905 JohannsenELambertPFEpigenetics of human papillomavirusesVirology44520521220132395323010.1016/j.virol.2013.07.016 SzalmásAKónyaJEpigenetic alterations in cervical carcinogenesisSemin Cancer Biol1914415220091942947710.1016/j.semcancer.2009.02.011 JhaAKSharmaVNikbakhtMJainVSehgalACapalashNKaurJA comparative analysis of methylation status of tumor suppressor genes in paired biopsy and serum samples from cervical cancer patients among north Indian populationGenetika5225525920162721504110.7868/s0016675816010070 de la Cruz-HernandezEPerez-CardenasEContreras-ParedesACantuDMoharALizanoMDuenas-GonzalezAThe effects of DNA methylation and histone deacetylase inhibitors on human papillomavirus early gene expression in cervical cancer, an in vitro and clinical studyVirol J41820071732426210.1186/1743-422X-4-18 LendvaiÁJohannesFGrimmCEijsinkJJWardenaarRVoldersHHKlipHGHollemaHJansenRCSchuuringEGenome-wide methylation profiling identifies hypermethylated biomarkers in high-grade cervical intraepithelial neoplasiaEpigenetics71268127820122301886710.4161/epi.22301 KellyHBenaventeYPavonMADe SanjoseSMayaudPLorinczATPerformance of DNA methylation assays for detection of high-grade cervical intraepithelial neoplasia (CIN2+): A systematic review and meta-analysisBr J Cancer12195496520193161603710.1038/s41416-019-0593-4 Hoppe-SeylerKBosslerFBraunJAHerrmannALHoppe-SeylerFThe HPV E6/E7 Oncogenes: Key factors for viral carcinogenesis and therapeutic targetsTrends Microbiol2615816820182882356910.1016/j.tim.2017.07.007 BachmanKEParkBHRheeIRajagopalanHHermanJGBaylinSBKinzlerKWVogelsteinBHistone modifications and silencing prior to DNA methylation of a tumor suppressor geneCancer Cell3899520031255917810.1016/s1535-6108(02)00234-9 ClarkSJMelkiJDNA methylation and gene silencing in cancer: Which is the guilty party?Oncogene215380538720021215440010.1038/sj.onc.1205598 TaharaTShibataTYamashitaHNakamuraMYoshiokaDOkuboMHirataIArisawaTChronic nonsteroidal anti-inflammatory drug (NSAID) use suppresses multiple CpG islands hyper methylation (CIHM) of tumor suppressor genes in the human gastric mucosaCancer Sci1001192119720091943288610.1111/j.1349-7006.2009.01175.x ClarkeMAGradissimoASchiffmanMLamJSollecitoCCFettermanBLoreyTPoitrasNRaine-BennettTRCastlePEHuman papillomavirus DNA methylation as a biomarker for cervical precancer: Consistency across 12 genotypes and potential impact on management of HPV-positive womenClin Cancer Res242194220220182942022210.1158/1078-0432.CCR-17-3251 CuschieriKRoncoGLorinczASmithLOgilvieGMirabelloLCarozziFCubieHWentzensenNSnijdersPEurogin roadmap 2017: Triage strategies for the management of HPV-positive women in cervical screening programsInt J Cancer14373574520182934111010.1002/ijc.31261 ZhaoZZhangXZhaoXCaiJWuNYWangJSOX1 and PAX1 are hypermethylated in cervical adenocarcinoma and associated with better prognosisBiomed Res Int2020398152920203337672210.1155/2020/3981529 MengMSangLWangXS100 calcium binding protein A11 (S100A11) promotes the proliferation, migration and invasion of cervical cancer cells, and activates Wnt/β-catenin signalingOnco Targets Ther128675868520193169542610.2147/OTT.S225248 TianTLiXHuaZMaJWuXLiuZChenHCuiZS100A7 promotes the migration, invasion and metastasis of human cervical cancer cells through epithelial-mesenchymal transitionOncotarget8249642497720172821256410.18632/oncotarget.15329 ChokchaichamnankitDWatcharatanyatipKSubhasitanontPWeeraphanCKeeratichamroenSSritanaNKantathavornNDiskul-Na-AyudthayaPSaharatKChantaraampornJUrinary biomarkers for the diagnosis of cervical cancer by quantitative label-free mass spectrometry analysisOncol Lett175453546820193118676510.3892/ol.2019.10227 TomiyamaNIkedaRNishizawaYMasudaSTajitsuYTakedaYS100A16 up-regulates Oct4 and Nanog expression in cancer stem-like cells of Yumoto human cervical carcinoma cellsOncol Lett159929993320182992836610.3892/ol.2018.8568 LeeHSYunJHJungJYangYKimBJLeeSJYoonJHMoonYKimJMKwonYIIdentification of differesntially-expressed genes by DNA methylation in cervical cancerOncol Lett91691169820152578902510.3892/ol.2015.2917 LiuMYZhangHHuYJChenYWZhaoXNIdentification of key genes associated with cervical cancer by comprehensive analysis of transcriptome microarray and methylation microarrayOncol Lett1247347820162734716710.3892/ol.2016.4658 MaXLiuJWangHJiangYWanYXiaYChengWIdentification of crucial aberrantly methylated and differentially expressed genes related to cervical cancer using an integrated bioinformatics analysisBioscience Rep40BSR2019436520203236878410.1042/BSR20194365 DeanBScarrEPossible involvement of muscarinic receptors in psychiatric disorders: A focus on schizophrenia and mood disordersCurr Mol Med1525326420152581785810.2174/1566524015666150330144821 LiangLHuangJYaoMLiLJinXJCaiXYGNG4 promotes tumor progression in colorectal cancerJ Oncol2021993198420213469117910.1155/2021/9931984 ZhaoHShengDQianZYeSChenJTangZIdentifying GNG4 might play an important role in colorectal cancer TMBCancer Biomark3243545020213427589210.3233/CBM-203009 WenSPengWChenYDuXXiaJShenBZhouGFour differentially expressed genes can predict prognosis and microenvironment immune infiltration in lung cancer: A study based on data from the GEOBMC Cancer2219320223518474810.1186/s12885-022-09296-8 WangXHouQZhouXLAMA4 expression is activated by zinc finger E-box-binding homeobox 1 and independently predicts poor overall survival in gastric cancerOncol Rep401725173320183001586110.3892/or.2018.6564 WangJFWangYZhangSWChenYYQiuYDuanSYLiBPChenJQExpression and prognostic analysis of integrins in gastric cancerJ Oncol2020886222820203333555010.1155/2020/8862228 ItohNOhtaHFgf10: A paracrine-signaling molecule in development, disease, and regenerative medicineCurr Mol Med1450450920142473052510.2174/1566524014666140414204829 DengXWeiWHuangNShiYHuangMYanYLiDYiJWangXTumor repressor gene chondroadherin oppose migration and proliferation in hepatocellular carcinoma and predicts a good survivalOncotarget8602706027920172894796910.18632/oncotarget.19811 LinYLWangYPLiHZZhangXAberrant promoter methylation of PCDH17 (Protocadherin 17) in serum and its clinical significance in renal cell carcinomaMed Sci Monit233318332320172868823210.12659/msm.902077 XuWXuMWangLZhouWXiangRShiYZhangYPiaoYIntegrative analysis of DNA methylation and gene expression identified cervical cancer-specific diagnostic biomarkersSignal Transduct Target Ther45520193187177410.1038/s41392-019-0081-6 BhatSKabekkoduSPVargheseVKChakrabartySMallyaSPRottiHPandeyDKushtagiPSatyamoorthyKAberrant gene-specific DNA methylation signature analysis in cervical cancerTumour Biol39101042831769457320172835129810.1177/1010428317694573 VargheseVKShuklaVKabekkoduSPPandeyDSatyamoorthyKDNA methylation regulated microRNAs in human cervical cancerMol Carcinog5737038220182907723410.1002/mc.22761

The magnitude and statistical significance of the methylation difference between the cervical cancer group and the normal group (A-a and B-a) were reflected by the volcanic map (P<0.01, The Δ β value >0 corresponds to hypermethylation, and the Δ β value <0 to hypomethylation). (A-b and B-b) The hierarchical cluster analysis is depicted by the heatmap of differentially methylated CpG company sites between the cervical cancer group and the normal group (DNA methylation levels are color-coded. Blue represents CpG with the lowest methylation level and red represents CpG with the highest methylation level). HPV, human papillomavirus.

GO functional annotation analysis of the differentially expressed methylation sites from (A) HPV-16- and (B) HPV-18-positive samples. GO, gene ontology; HPV, human papillomavirus.

Methylation difference sites were associated with the top 20 signalling pathways in tested (A) HPV-16- and (B) HPV-18-positive samples (P<0.01 and FDR <0.05). HPV, human papillomavirus.

The top 10 methylation difference sites of HPV-16 group screened using the 850 k methylation chip.

Sequence Probe ID Normal group average methylation level HPV-16 group average methylation level Difference between two groups P-value FDR Gene
1 cg24575234 0.0700±0.0234 0.6615±0.1325 0.592 0.000318 0.00888b CHRM2
2 cg20818778 0.1310±0.0074 0.6913±0.1583 0.560 0.000888 0.01466a GNG4
3 cg14289461 0.1687±0.0852 0.7237±0.1104 0.555 0.000511 0.011113a LAMA4
4 cg06818777 0.0654±0.0250 0.6105±0.1142 0.545 0.000241 0.007846b CHAD
5 cg08361126 0.0901±0.0376 0.6207±0.1280 0.531 0.000514 0.011149a ITGA8
6 cg06381931 0.1690±0.0502 0.6946±0.0608 0.526 4.34E-05 0.004157b COL11A1
7 cg08976810 0.0369±0.0098 0.5484±0.1995 0.511 0.003815 0.033664a FGF10
8 cg20881548 0.2406±0.0908 0.7451±0.0314 0.505 0.000138 0.006232b IGF1
9 cg08264401 0.2143±0.0175 0.7006±0.0315 0.486 1.74E-06 0.002313b TEK
10 cg25459558 0.0932±0.0189 0.5783±0.1569 0.485 0.00171 0.02099a COL11A2

P-values were obtained using a t-test for comparisons between normal control and HPV-16-positive groups. Variance is uneven; following logarithmic conversion, the data meets the parameter test conditions.

aP≤0.05 and

bP≤0.01, compared with the normal group. HPV, human papillomavirus.

The top 10 methylation difference sites of HPV-18 group screened using the 850 k methylation chip.

Sequence Probe ID G5 normal group average methylation level G3 cervical cancer HPV-18-positive group average methylation level Difference between two groups P-value FDR Gene
1 cg03520644 0.1225±0.0342 0.3966±0.0241 0.2741 5.98E-05 0.009034a COL11A1
2 cg25792518 0.3841±0.0564 0.7650±0.0315 0.3809 8.44E-05 0.009827a CHAD
3 cg06958829 0.1182±0.0121 0.6611±0.0797 0.5430 4.17E-05 0.008166a CHAD
4 cg00172849 0.1219±0.0440 0.3997±0.0141 0.2778 8.70E-05 0.009885a COL11A1
5 cg19707040 0.1011±0.0136 0.4713±0.0504 0.3702 3.63E-05 0.007918a CTNNA2
6 cg02501779 0.1815±0.0145 0.7214±0.0734 0.5399 3.00E-05 0.007761a CBLN4
7 cg19679123 0.4993±0.0466 0.8142±0.0246 0.3149 8.42E-05 0.009826a SMAD3
8 cg14427009 0.2460±0.0267 0.7325±0.0381 0.4865 5.43E-06 0.006396a PCDH17
9 cg25993718 0.0971±0.0343 0.5181±0.0515 0.4211 4.20E-05 0.00819a CBLN4
10 cg27423357 0.3920±0.0364 0.7782±0.0226 0.3862 1.20E-05 0.006597a FLT1

P-values are obtained by using t test for comparisons between normal control and HPV-18-positive groups. Variance is uneven; following logarithmic conversion, the data meets the parameter test conditions.

aP≤0.01, compared with the normal group. HPV, human papillomavirus.

The top 10 methylation difference sites of HPV-18 group screened using the 850 k methylation chip were verified in TCGA.

Sequence Probe ID Normal group Cervical cancer group HPV-18-positive t-test P-value Gene
1 cg03520644 0.616±0.0312 0.7381±0.0467 -20.847 <0.001c COL11A1
2 cg25792518 0.2926±0.1046 0.8885±0.0626 -8.464 0.001b CHAD
3 cg06958829 0.1529±0.0501 0.7040±0.0587 -12.366 <0.001c CHAD
4 cg00172849 0.1886±0.0464 0.6807±0.0938 -8.146 0.001b COL11A1
5 cg19707040 0.0350±0.084 0.5199±0.0944 -8.523 0.001b CTNNA2
6 cg02501779 0.1541±0.0400 0.6263±0.1349 -5.814 0.004b CBLN4
7 cg19679123 0.4073±0.1601 0.8318±0.0875 -4.030 0.016a SMAD3
8 cg14427009 0.2136±0.0373 0.6362±0.0527 -11.327 <0.001c PCDH17
9 cg25993718 0.1236±0.0384 0.5085±0.1573 -4.118 0.015a CBLN4
10 cg27423357 0.4072±0.0873 0.7177±0.1584 -2.975 0.041a FLT1

P-values are obtained by using t test for comparisons between normal control and HPV-18-positive groups. Variance is uneven; following logarithmic conversion, the data meets the parameter test conditions.

aP≤0.05,

bP≤0.01 and

cP<0.001, compared with the normal group. HPV, human papillomavirus.

Verification of different top 9 methylation CpG sites in the HPV-16 group using pyrosequencing.

Sequence Probe ID Normal group Cervical cancer group HPV-16 t-test P-value Gene
1 cg24575234 8.68±1.46 33.69±5.71 9.495 <0.001c CHRM2
2 cg20818778 14.92±0.603 28.17±1.33c 1.916 0.092a GNG4
3 cg14289461 10.62±3.91 48.76±16.40 5.059 0.001b LAMA4
4 cg06818777 7.86±2.17 44.13±2.54c 2.775 0.024a CHAD
5 cg08361126 8.46±1.17 24.53±11.36c 2.833 0.025a ITGA8
6 cg06381931 26.16±4.54 59.35±10.52 6.474 <0.001c COL11A1
7 cg08976810 9.39±1.68 52.08±12.63 7.492 <0.001c FGF10
8 cg20881548 36.35±4.50 60.61±10.68 4.684 0.002b IGF1
9 cg08264401 23.79±2.77 44.38±12.60 3.568 0.007b TEK

P-values are obtained by using t test for comparisons between normal control and HPV-16-positive groups. Variance is uneven; following logarithmic conversion, the data meets the parameter test conditions.

aP≤0.05,

bP≤0.01 and

cP<0.001, compared with the normal group. HPV, human papillomavirus.