Cervical cancer is the most common gynecological malignancy and is especially prevalent in young women [1]. While cervical cancer is often curable if detected at an early stage, the overall survival in patients with advanced cervical cancer remains poor. Therefore, better treatments and molecular markers related to carcinogenesis and progression that can help clinical practice are needed urgently. Screening is key for the early detection of precancerous lesions and cancer. ThinPrep cytologic tests (TCT) and human papillomavirus (HPV) detection are 2 important and commonly used methods to screen for cervical cancer [2]. TCTs are inexpensive, but their performance depends on several subjective human factors, which can lead to misdiagnosis and missed diagnosis, and poor accuracy. The sensitivity for detecting HPV is high, but most cases of HPV infection are transient, leading to poor specificity. Cytological screening strategies combined with HPV testing have contributed to a substantial reduction in the incidence and mortality of cervical cancer [3]. However, overall, the positive predictive effects are moderate, as only a small subset of people with HPV will eventually develop cervical cancer. Gene methylation analysis is a nonmorphological molecular detection method, which can provide an objective and potential shunt method for high-risk HPV-positive women. Triage tests using DNA methylation may be helpful to differentiate women who are at high risk of rapidly developing cervical cancer from women with low risk.
DNA methylation of tumor-related genes is closely related to early tumor progression, in a cancer type specific manner, and may provide information that can aid early screening of cervical cancer [4]. Aberrant promoter methylation, and epigenetic changes in tumor suppressor genes, either inactivating or silencing, are associated with tumorigenesis and progression [5,6,7,8]. In humans, DNA methylation occurs almost exclusively at the carbon 5 position on cytosine residues in CpG dinucleotides, which are concentrated in distinct GC-rich regions called “CpG islands” (CGIs) [8]. CpG island hypermethylation of tumor suppressor genes has been linked to the development of numerous human cancers, including cervical cancer [9, 10]. The Wnt/β-catenin signaling pathway is known to be regulated by several secreted antagonists, such as members of the Dickkopf (DKK) protein family (DKK1, DKK2, DKK3, and DKK4), and has a function in embryonic development and tumorigenesis [11].
Although DKK2 plays a significant role in many tumors, its correlation with clinicopathological characteristics such as lymph node metastasis and HPV infection is still unknown in cervical cancer. Here, we evaluated the mRNA expression profiles of
To compare gene expression differences between healthy donors and patients with cervical cancer, all clinicopathological data related to DKK2 expression profiles (GSE6791, GSE7803, GSE9750, and GSE7410) were carefully selected from the Gene Expression Omnibus (GEO) database, The Cancer Genome Atlas (TCGA) database (
Ect1/E6E7 (catalog No. AC39960) normal cervical epithelial cell and HPV-positive HeLa (catalog No. TCHu187), HPV-positive CaSki (catalog No. BNCC338223), HPV-negative HT-3 (catalog No. HTB-32), and HPV-negative C33A (catalog No. BNCC337882) cervical cancer cell lines were screened. Cells were purchased from American Type Culture Collection (ATCC) or BeNa Culture Collection (Beijing, China) and cultured in RPMI-1640 medium (Gibco).
We analyzed a series of cervical specimens with histological diagnoses collected between May 2016 and November 2017 at Liaocheng People's Hospital (Shandong, China). At the time of surgical resection or cervical biopsy, 79 invasive cervical cancer tissues and 63 cervical intraepithelial neoplasia (CIN) tissues were obtained. Of the CIN tissues, 25 were classified as low-grade CIN (low-grade squamous intraepithelial lesions, LSIL) and 38 as high-grade CIN (high-grade squamous intraepithelial lesions, HSIL), and the latter included both CIN2 and CIN3. Normal cervical tissue (29 cases) was retrieved from patients with uterine leiomyomas who underwent a hysterectomy. All samples were stored at −80 °C until analysis. Genotyping was used to detect HPV16 and HPV18 viral DNA [20]. All cases had a confirmed diagnosis, as well as confirmed HPV typing and histopathological type. Exclusion criteria included samples from patients undergoing chemotherapy or radiotherapy, immunocompromised patients, pregnancy, diagnosis of other cancers, uterine cervix operations, and chronic or acute viral infection. Documented written informed consent was provided by all patients for tissue sample collection. The present study was approved by the medical ethics committee of Liaocheng People's Hospital, China (No. 2016022), and was consistent with relevant national regulations and laws, including the People's Republic of China Regulations for the Management of Medical Institutions (promulgated by the Order No. 149 of the State Council on February 26, 1994; and revised in accordance with the Decision of the State Council on Amending Some Administrative Regulations on February 6, 2016), and international medical ethics documents, including the principles of the World Medical Association (WMA) Declaration of Helsinki and its contemporary amendments (2013).
DKK2 mRNA expression was examined by reverse transcription-polymerase chain reaction (RT-PCR) analysis as previously described [18]. Total RNA was extracted from cell lines, CIN, normal cervical tissues, and tumor samples using TRIzol reagent (Invitrogen).
Primer pairs used for mRNA expression and the size of PCR products
5′-GTACCAAGGACTGGCATTCG-3′ (F) | 169 bp | |
5′-CCAGCAAGAGCACAAGAGGAA-3′ (F) | 114 bp | |
5′-AGAGTTAAATCGTCGAGATTTC-3′ (F) | 146 bp | |
5′-GGAGAGTTAAATTGTTGAGATTTT-3′ (F) | 149 bp |
DNA samples were extracted from tissues and cell lines. Genomic DNA (1 μg) was bisulfite-modified and purified following the guidelines described by the manufacturer. After desulfonation with NaOH and precipitation with ethanol, the final products were dissolved in 20 μL of Tris-EDTA buffer and stored at −80 °C.
HeLa, CaSki, HT-3, and C33A cells (2 × 105) were treated with 10 μmol 5-Aza-2′-deoxycytidine (5-Aza-dC) (Sigma) for 96 h. After adding the demethylation treatment, the medium was changed every 24 h. RNA and DNA were isolated according to routine procedures.
SPSS Statistics for Windows (version 17.0) was used for all statistical analyses. Student
We examined DKK2 mRNA expression in 79 normal cervical tissues (50 by data mining, 29 by RT-PCR), 9 cervical cancer lines (HT-3, C4-I, CaSki, MS751, C33A, SiHa, SW756, ME-180, and HeLa cell lines by data mining, HeLa, HT-3, C33A, and CaSki cell lines by RT-PCR), and 498 cervical cancer samples (419 by data mining, 79 by RT-PCR). Based on data from the GEO database (
Figure 1
DKK2 mRNA expression in cervical cells and tissues. (

We then examined DKK2 mRNA expression both in cervical cancer cell lines and in tissue samples. Our RT-PCR analysis indicated that DKK2 mRNA expression was weak in HPV-negative C33A cells, and silenced in HPV-negative HT-3 cells and HPV-positive CaSki and HeLa cells, as compared with the normal cervical epithelial cell line Ect1/E6E7 (
Normal | 29 | 29/29 (100) | 0/29 (0) | ||
LSIL | 25 | 25/25 (100) | 0.00003 | 1/25 (4.0) | 0.005 |
HSIL | 38 | 21/38 (55.2) | 0.001 | 9/38 (23.7) | 0.037 |
Cervical cancer | 79 | 25/79 (31.6) | 0.014 | 41/79 (51.9) | 0.004 |
DKK2, Dickkopf 2; HSIL, high-grade squamous intraepithelial lesions; LSIL, low-grade squamous intraepithelial lesions.
We then examined the status of
Figure 2

We then examined
The correlation between
DKK2 mRNA expression,
Clinical stage | |||||||
Stage IA2 | 28 | 13 (46.4%) | 15 (53.6%) | 0.099 | 12 (42.9%) | 16 (57.1%) | 0.32 |
Stage IB1+IB2 | 41 | 9 (22.0%) | 32 (78.0%) | 22 (53.7%) | 19 (46.3%) | ||
Stage IIA | 10 | 3 (30.0%) | 7 (70.0%) | 7 (70.0%) | 3 (30.0%) | ||
Histology | |||||||
Squamous cell carcinoma | 70 | 21 (30.0%) | 49 (70.0%) | 0.38 | 36 (51.4%) | 34 (48.6%) | 0.55 |
Adenocarcinoma | 9 | 4 (44.4%) | 5 (55.6%) | 5 (55.6%) | 4 (44.4%) | ||
Tumor grade | |||||||
Grade 1 | 15 | 5 (33.3%) | 10 (66.7%) | 0.81 | 6 (40.0%) | 9 (60.0%) | 0.17 |
Grade 2 | 42 | 12 (28.6%) | 30 (71.4%) | 20 (47.6%) | 22 (52.4%) | ||
Grade 3 | 22 | 8 (36.4%) | 14 (63.6%) | 15 (68.2%) | 7 (31.8%) | ||
Lymph node metastasis | |||||||
Negative | 62 | 23 (37.1%) | 39 (62.9%) | 0.047 | 28 (45.2%) | 34 (54.8%) | 0.022 |
Positive | 17 | 2 (11.8%) | 15 (88.2%) | 13 (76.5%) | 4 (23.5%) |
DDK2, Dickkopf 2.
Relationship between HR-HPV infection type and DKK2 methylation rate
Cervical cancer | 0.015 | |||
HPV16/18(+) | 48 | 32 | 66.7 | |
HPV16/18(−) | 25 | 9 | 36 | |
Precancerous lesions | 0.48 | |||
HPV16/18(+) | 29 | 7 | 24.1 | |
HPV16/18(−) | 23 | 3 | 15 |
HPV, human papillomavirus; HR-HPV, high-risk HPV.
In the present study, we investigated
The activated Wnt/β-catenin signaling pathway plays an important role in cervical cancer. CpG island promoter hypermethylation has been shown to inactivate extracellular Wnt antagonists in cervical cancer [22, 23]. CpG island methylation can be widely found in the human genome, but only a subset of loci plays important roles in tumorigenesis. These implicated genes are known to be involved in multiple cellular signaling pathways such as apoptosis, cell cycle regulation, development, differentiation, invasion, and metastasis. As CpG island methylation occurs early and methylated alleles can be detected in a sensitive manner in carcinogenesis, detection of methylation may be a promising tool for early detection of cancer [24]. DKK2, a Wnt antagonist, contributes to tumorigenesis in multiple cancers [12, 13, 25,26,27,28]. In the present study, RT-PCR showed reduced expression of DKK2 mRNA in cervical cancer compared with either normal cervical samples or HSIL samples. Further, through MSP analysis
Methylation of the
HPV infection accounts for over 90% of cervical cancer cases, and the high-risk types of HPV are associated with 87%–88% of squamous cell carcinomas [30]. A triage test using DNA methylation may be helpful to differentiate women who are at high risk of developing cervical cancer rapidly from women with low risk. In the case of HSIL, the combination of HPV genotyping and methylation marker analysis can not only overcome the limitations of cytological examination but also increase diagnosis accuracy. The present study indicated that among HR-HPV-positive cervical cancer patients, the methylation rate of
We acknowledge that the present study also has some limitations. In the precancerous lesion group, the difference was not statistically significant. The correlation between
Figure 1

Figure 2

DKK2 mRNA expression and methylation in cervical neoplasms
Normal | 29 | 29/29 (100) | 0/29 (0) | ||
LSIL | 25 | 25/25 (100) | 0.00003 | 1/25 (4.0) | 0.005 |
HSIL | 38 | 21/38 (55.2) | 0.001 | 9/38 (23.7) | 0.037 |
Cervical cancer | 79 | 25/79 (31.6) | 0.014 | 41/79 (51.9) | 0.004 |
Relationship between HR-HPV infection type and DKK2 methylation rate
Cervical cancer | 0.015 | |||
HPV16/18(+) | 48 | 32 | 66.7 | |
HPV16/18(−) | 25 | 9 | 36 | |
Precancerous lesions | 0.48 | |||
HPV16/18(+) | 29 | 7 | 24.1 | |
HPV16/18(−) | 23 | 3 | 15 |
DKK2 mRNA expression, DKK2 promoter methylation and clinicopathological characteristics of cervical cancer
Clinical stage | |||||||
Stage IA2 | 28 | 13 (46.4%) | 15 (53.6%) | 0.099 | 12 (42.9%) | 16 (57.1%) | 0.32 |
Stage IB1+IB2 | 41 | 9 (22.0%) | 32 (78.0%) | 22 (53.7%) | 19 (46.3%) | ||
Stage IIA | 10 | 3 (30.0%) | 7 (70.0%) | 7 (70.0%) | 3 (30.0%) | ||
Histology | |||||||
Squamous cell carcinoma | 70 | 21 (30.0%) | 49 (70.0%) | 0.38 | 36 (51.4%) | 34 (48.6%) | 0.55 |
Adenocarcinoma | 9 | 4 (44.4%) | 5 (55.6%) | 5 (55.6%) | 4 (44.4%) | ||
Tumor grade | |||||||
Grade 1 | 15 | 5 (33.3%) | 10 (66.7%) | 0.81 | 6 (40.0%) | 9 (60.0%) | 0.17 |
Grade 2 | 42 | 12 (28.6%) | 30 (71.4%) | 20 (47.6%) | 22 (52.4%) | ||
Grade 3 | 22 | 8 (36.4%) | 14 (63.6%) | 15 (68.2%) | 7 (31.8%) | ||
Lymph node metastasis | |||||||
Negative | 62 | 23 (37.1%) | 39 (62.9%) | 0.047 | 28 (45.2%) | 34 (54.8%) | 0.022 |
Positive | 17 | 2 (11.8%) | 15 (88.2%) | 13 (76.5%) | 4 (23.5%) |
Primer pairs used for mRNA expression and the size of PCR products
5′-GTACCAAGGACTGGCATTCG-3′ (F) |
169 bp | |
5′-CCAGCAAGAGCACAAGAGGAA-3′ (F) |
114 bp | |
5′-AGAGTTAAATCGTCGAGATTTC-3′ (F) |
146 bp | |
5′-GGAGAGTTAAATTGTTGAGATTTT-3′ (F) |
149 bp |
A One Health approach to antimicrobial resistance Sialic acid: an attractive biomarker with promising biomedical applications Antibiotic resistance, biofilm forming ability, and clonal profiling of clinical isolates of Staphylococcus aureus from southern and northeastern IndiaAn 85-amino-acid polypeptide from Myrmeleon bore larvae (antlions) homologous to heat shock factor binding protein 1 with antiproliferative activity against MG-63 osteosarcoma cells in vitroLong noncoding and micro-RNA expression in a model of articular chondrocyte degeneration induced by stromal cell-derived factor-1 Promoter methylation analysis of DKK2 may be a potential biomarker for early detection of cervical cancer