High DHCR7 Expression Predicts Poor Prognosis for Cervical Cancer

Zou, Jing; Liu, Sha; Long, Jian; Yan, Bingbing

doi:https://doi.org/10.1155/2022/8383885

Computational and Mathematical Methods in Medicine

On this page

Abstract Background Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 8383885 | https://doi.org/10.1155/2022/8383885

High DHCR7 Expression Predicts Poor Prognosis for Cervical Cancer

Jing Zou,¹Sha Liu,¹Jian Long,¹and Bingbing Yan²

Academic Editor: Kathiravan Srinivasan

Received04 Dec 2021

Revised19 Jul 2022

Accepted26 Jul 2022

Published16 Sept 2022

Abstract

DHCR7 is a rate-limiting enzyme in cholesterol synthesis. The expression pattern and prognostic value of DHCR7 in cervical cancer are unknown. We investigated the relationship between DHCR7 expression and clinicopathological features of cervical cancer patients. The dataset was acquired from TCGA database. The Wilcoxon rank sum test was used to explore DHCR7 expression level in cervical cancer. The Kruskal-Wallis test and the logistic regression were performed to estimate the association between the DHCR7 and clinical features. The Kaplan-Meier and Cox regression analyses were used to evaluate factors that affect cervical cancer prognosis. GSEA was used to screen the DHCR7-related pathways. We found that DHCR7 was increased in cervical cancer samples and increased DHCR7 was correlated with advanced T stage, lymph node invasion, and clinical stage (). Patients with elevated DHCR7 levels had poorer overall survival (), progression-free interval (), and disease-specific survival (). Cox analysis revealed that DHCR7 was an independent prognostic factor in cervical cancer (). WNT activated receptor activity, G2/M checkpoint, mTORC1 signaling, KRAS signaling, regulation of cholesterol biosynthetic, FGF signaling, T-cell receptor signaling, JAK/STAT signaling cascade T cell activation, and macrophage migration were enriched in high DHCR7 phenotype. Our data also showed that DHCR7 moderately correlates with T-cell infiltration, including CD8⁺ T-cells. Conclusion. Increased DHCR7 expression is associated with poor survival in cervical cancer.

1. Background

Cervical cancer is one of the most common cancers worldwide and is associated with high mortality rate. It is a major health issue for women, especially in developing countries [1]. Surgery, radiotherapy, and chemotherapy are the main treatment strategies for cervical cancer [2]. Although the survival of cervical cancer has been improved with the combination of the above treatments, the prognosis of high-risk patients remains poor, especially those with metastasis and recurrent diseases [3]. Discovery of biomarkers for diagnosis and therapy may help to improve the prognosis of cervical cancer.

DHCR7 (7-dehydrocholesterol reductase) is a rate-limiting enzyme in cholesterol synthesis, which is mainly involved in the conversion of 7-dehydrocholesterol (7-DHC) to cholesterol [4]. It has been pointed out that cholesterol has an impact on the prognosis of several types of cancers, including cervical cancer [5]. DHCR7 deficiency causes Smith-Lemli-Opitz syndrome (SLOS), which is an autosomal recessive syndrome charactered by multiple congenital abnormalities, including morphogenic, mental retardation and behavioral disorders, due to lack of cholesterol-derived steroid hormones [6]. DHCR7 is also considered to be a vitamin D-related gene which may have an influence on calcium homeostasis, bone health, and various cancers [7–9]. Carvalho et al. [10] showed that DHCR7 polymorphisms might increase the risk of thyroid cancer via its effects on circulating vitamin D levels. However, there is no report about the prognostic value of DHCR7 in cervical cancer.

Here, we used the RNA-Seq data and clinical information from TCGA database to explore the prognostic value of DHCR7 in cervical cancer. Firstly, we compared the DHCR7 expression between cervical cancer tissues and normal tissues. Then, we studied the relationship between DHCR7 expression and clinicopathological characteristics. Cervical cancer patients were clustered into 2 groups based on the DHCR7 expression level. The Kaplan-Meier method and Cox regression analyses were then utilized to compare the survival of patients with different DHCR7 expression levels. Gene ontology (GO) and KEGG enrichment analyses were used to explore the biological role of DHCR7. Furthermore, DHCR7-related biological pathways involved in cervical cancer were screened by gene set enrichment analysis (GSEA). Additionally, we also explored its role in tumor immunity.

2. Materials and Methods

2.1. RNA Sequencing Data Source and Data Preprocessing

RNA sequencing data were obtained from TCGA (uhttps://http://portal.gdc.cancer.gov/) (normal cervical tissues and cervical cancer tissues) and GTEx (normal cervical tissues). RNA-seq data were analyzed to determine DHCR7 expression level via UCSC’s Xena (https://xena.ucsc.edu/). UCSC’s Toil algorithm [11] was used to diminish the batch effects. Patients with incomplete clinical information or gene expression data were excluded from the analysis. Since the present study met the publication guidelines stated by TCGA (https://cancergenome.nih.gov/publications/publicationguidelines) and the data were obtained from TCGA, ethical approval and informed consent were not required.

2.2. Identification of Differentially Expressed Genes (DEGs)

Tumor samples were grouped into DHCR7-high and DHCR7-low expression groups based on the median DHCR7 expression level. The expression profiles of DHCR7-high and DHCR7-low expression groups were compared using DESeq2 (3.8) package of R [12]. Genes with |log2 fold change (FC)| >1 and adjust (adj.) -value <0.05 were considered DEGs. Gene expression differences were visualized using volcano plot and heat map.

2.3. Functional Enrichment Analysis

Gene ontology (GO) [13] functional analysis and KEGG [14] pathway analysis were performed using the clusterProfiler method (3.14.3) [15]. FDR (false discovery rate) <0.05 and adj. -value <0.05 were considered statistically significant.

2.4. Gene Set Enrichment Analysis (GSEA)

Gene set enrichment analysis (GSEA) [16] of DHCR7-high and DHCR7-low groups were conducted using the clusterProfiler package of R software. The analysis involved 1000 permutations for each analysis. Adj. -value <0.05, FDR -value <0.25, and normalized enrichment score (|NES|) ≥1 were employed to identify enriched pathways.

2.5. Inference of Infiltrating Immune Analysis

Immunocyte signatures involving 509 genes were used to predict the infiltration level of 24 immune cells in each tissue sample [17]. Immune infiltration level of immune cells was analyzed via single-sample gene set enrichment analysis (ssGSEA) [18] using GSVA (http://bioconductor.org/packages/release/bioc/html/GSVA.html) package on R. To assess the correlation between DHCR7 expression and immune cells infiltration levels, Spearman’s correlation analysis was performed. Immune cell infiltration in DHCR7-high and DHCR7-low expression groups was compared using the Wilcoxon rank sum test.

2.6. Protein-Protein Interaction (PPI) Network

To further study the functional interaction between DEGs, Metascape [19] online analysis (http://metascape.org) (update date: 2020-3-20) was used to construct a PPI network of the DEGs. Terms with , a minimum count of 3, and the enrichment factor >1.5 were regarded as significant. Hub genes and significant modules of this PPI network were extracted by Molecular Complex Detection (MCODE) [20] algorithm.

2.7. Detection of DHCR7 Expression by Western Blotting

Cervical cancer cells were cultured in RPMI 1640 medium (Gibco, Thermo Fisher Scientific, Inc., MA, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco, Thermo Fisher Scientific, Inc., MA, USA). Cells were incubated at 37°C in a 5% CO2 humidified incubator. Protein lysis buffer (Beyotime, Shanghai, China) was used to extract the cell lysates and bicinchoninic acid (BCA) protein assay kit (Beyotime, Shanghai, China) was used to determine the protein concentration. The protein samples were resolved by SDS-PAGE and then transferred to a polyvinylidene fluoride membrane (EMD Millipore, MA, USA). 5% bovine serum albumin (Beyotime, Shanghai, China) was used to block the membranes. Rabbit anti-DHCR7 antibody (1 : 2000, Invitrogen) and rabbit anti-β-actin were used to react with the protein samples at 4°C overnight. Then, goat anti-rabbit antibody (1 : 1000, Beyotime, Shanghai, China) was used to incubate the membranes at room temperature for 1 h. ECL plus reagent (Beyotime, Shanghai, China) was used to detect the brands.

2.8. Downregulation of DHCR7 Expression

About 2 × 10⁵ cells were seeded into each well of the 6-well plates and incubated at 37°C in a 5% CO2 humidified incubator overnight so that the cells were adherent to the plates. 100 pmol of siRNA of DHCR7 (si-DHCR7) or negative control and 7.5uL lipo3000 transfection reagent (Invitrogen) were added to each well and the cells were incubated for 48 h. Then, DHCR7 expression was detected by western blotting. Sequences of si-DHCR7: sense: 5- CUAUAUGAUGGGAAUUGAG -3, antisense: 5- GAUAUACUACCCUUAACUC -3.

2.9. Cell Proliferation Was Determined by CCK8 Assay

Cells were cultured in RPMI 1640 medium supplemented with 10% FBS by 96-well plates for 24 h, 48 h, and 72 h, respectively. Then, the medium was replaced by RPMI 1640 medium not containing serum and 20μl CCK8 reagent (Beyotime, Shanghai, China) was added to each well and incubated in dark for 2 h. The optical density (OD) at 490 nm was detected.

2.10. Statistical Analysis

Statistical analyses were performed using R software (v3.6.2). The Wilcoxon rank sum test was employed to evaluate DHCR7 expression in 304 cervical cancer and 13 normal samples. DHCR7 protein expression in normal vs. cervical cancer tissues was assessed using immunohistochemistry (IHC) images on Human Protein Atlas (HPA, https://www.proteinatlas.org). The “pROC” package on R was used for ROC (receiver operating characteristic) analysis to determine the diagnostic power of DHCR7 in cervical cancer [21]. Relationship between DHCR7 expression and clinicopathological features was tested using the Kruskal-Wallis test, Wilcoxon signed-rank test, and logistic regression analysis. Survival was analyzed using the Kaplan-Meier method and tested by log-rank test. Univariate and multivariate Cox regression analyses were used to analyze the impact of different characteristics on patient prognosis. Overall survival (OS), progression-free interval (PFI), and disease-specific survival (DSS) were used to evaluate the clinical outcomes. The HR (95% confidence interval; CI) was used to evaluate hazard risk for individual factors. A nomogram was used to predict patient prognosis using rms package on R based on multivariate analysis results. Calibration plots were used to assess the nomogram’s prediction accuracy. indicated statistically significant differences.

3. Results

3.1. Clinical Characteristics

Toil was used to unify TPM format of RNA-seq data from 317 samples. 304 samples were cervical cancer samples from TCGA and 13 samples were normal tissues from TCGA and GTEx. Patient characteristics, including age, metastasis (M stage), lymph node status (N stage), clinical stage, T stage (primary tumors), histological type, histological grade, radiation therapy, menopause status, primary therapy outcome, birth control pill history, and smoking history, were recorded. Median age at diagnosis was 46 years old. 52 (17.1%) adenosquamous cases and 252 (82.9%) squamous cell carcinoma cases were included. There were 162 (53.2%), 69 (22.6%), 45 (14.8%), and 21 (6.9%) cases of patients with stages I, II, III, and IV diseases, respectively. Sixty (19.7%) patients had lymph node invasion (N1) and 10 (3.3%) patients had distant metastasis (M1) (Table 1).

3.2. DHCR7 Expression and Clinicopathological Characteristics

The Wilcoxon rank sum test analysis of DHCR7 expression levels in cervical cancer samples vs. normal samples revealed that DHCR7 level was significantly higher in cervical cancer tissues (; Figure 1(a)). Assessment of DHCR7 protein levels on Human Protein Atlas showed that expression intensity of DHCR7 was weak in non-tumor cervical tissues and strong in cervical cancer tissues (Figure 1(d)). To obtain a more comprehensive evaluation of DHCR7 expression in cancer, we compared DHCR7 mRNA expression in the TCGA pan-cancer cohort which included 33 types of cancers. Data analysis showed that DHCR7 expression was elevated in 27 types of cancer types in comparison with normal samples (Figure 1(c)). ROC curve analysis indicated that DHCR7 expression efficiently discriminated between tumor and non-tumor cervical tissues, with an area under the ROC curve (AUC) of 0.821 (95% CI: 0.673-0.968) (Figure 1(b)). The ROC curve suggested that DHCR7 had strong diagnostic power for cervical cancer. Comparison of DHCR7 expression in cervical cancer patients with different clinicopathological features revealed that DHCR7 expression increased with advanced clinical stage () and T stage (primary tumors) (). DHCR7 was also upregulated in patients with lymph nodes invasion (N1) (). Patients with worsen primary therapy outcome (PD) also showed higher DHCR7 expression level () (Figures 2(a)–2(d)). However, DHCR7 expression showed no difference in patients with different metastatic status (M stage) or different histological types () (Figures 2(e) and 2(f)). Logistic regression analysis revealed that high DHCR7 levels were correlated with advanced T stage (OR =1.01 (1.00-1.02), ), lymph node invasion (OR =1.01 (1.00-1.02), ), metastasis (OR =1.01 (1.00-1.03), ), advance clinical stage (OR =1.01 (1.00-1.01) for stage II and stage III and stage IV vs. stage I, ), and primary therapy outcome (OR =1.01 (1.00-1.01) for PR and SD and PD vs. CR, ) (Table 2).

(a)

(b)

(c)

(d)

Figure 1

DHCR7 expression in tumor and normal tissues. (a) DHCR7 mRNA expression in cervical cancer samples and normal cervical samples. (b) The ROC curve of DHCR7 in cervical cancer. (c) Expression of DHCR7 in 33 types of cancers and non-cancer samples. ns, . , . , . , . (d) Representative immunohistochemistry images of DHCR7 expression in cervical cancer tissues (intensitive, strong. Quantity, >75% positive) and normal tissues (intensitive, weak. Quantity, <25% positive) (Human Protein Atlas). DHCR7 protein expression was weak in non-tumor tissues and strong in cervical cancer tissues.

(a)

(b)

(c)

(d)

(e)

(f)

3.3. Survival Outcomes and Multivariate Analysis

304 cervical cancer cases with clinical outcome information were grouped into DHCR7-low and DHCR7-high groups based on the median DHCR7 expression level. The Kaplan-Meier analyses showed that the OS (P =0.021), PFI (P =0.002), and DSS () (Figures 3(a), 3(f), and 3(k)) of patients in DHCR7-low group were significantly better than patients in DHCR7-high group. We further conducted a subgroup analysis in patients with different clinical characteristics. In squamous cell carcinoma subgroup, patients with high DHCR7 expression showed obviously worse OS (), PFI (), and DSS () (Figures 3(b), 3(g), and 3(l)). In adenosquamous subgroup, patients with high DHCR7 also showed shorter PFI () (Figure S1B). However, there were no OS or DSS differences between patients with different DHCR7 expression levels (Figure S1A, C). And we also found that in lymph nodes invasion subgroup (N1), high DHCR7 patients had worse OS, PFI, and DSS than low DHCR7 patients () (Figures 3(c), 3(h), and 3(m)). In patients who received radiation therapy, OS, PFI, and DSS were significantly shorter for high-DHCR7 group () (Figures 3(d), 3(i), and 3(n)). The hazard ratio of different clinicopathological features to OS, PFI, and DSS was showed in the forest plots (Figures 3(e), 3(j), and 3(o)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

(n)

(o)

Figure 3

The prognostic value of DHCR7 for cervical cancer. (a) Kaplan-Meier analysis of OS in low- and high-DHCR7 groups. (b–d) OS in clinicopathological subgroups of squamous cell carcinoma, lymph node metastasis, and radiation therapy. (e) The hazard ratio of different clinicopathological features to OS in different clinicopathological subgroups. (f) Kaplan-Meier analysis of PFI in low- and high-DHCR7 groups. (g–i) PFI in clinicopathological subgroups of squamous cell carcinoma, lymph node metastasis, and radiation therapy. (j) The hazard ratio of different clinicopathological features to PFI in different clinicopathological subgroups. (K) Kaplan-Meier analysis of DSS in low- and high-DHCR7 groups. (l–n) DSS in clinicopathological subgroups of squamous cell carcinoma, lymph node metastasis, and radiation therapy. (o) The hazard ratio of different clinicopathological features to DSS in different clinicopathological subgroups. OS: overall survival; DSS: disease-specific survival; PFI: progression-free interval; DFI: disease-free interval.

Univariate Cox regression analysis showed that high-DHCR7 patients had worse OS (HR 1.758 (1.088-2.841), ), PFI (HR 2.224 (1.352-3.657), ), and DSS (HR = 2.258 (1.271-4.013), ) than low-DHCR7 patients (Table 3(a)). Multivariate Cox regression analysis was further performed to investigate the independent prognostic factors for cervical cancer. It was indicated that high DHCR7 expression level was independently associated with unfavorable OS, PFI, and DSS () (Table 3(b)). DHCR7 and other independent clinical risk factors (primary therapy outcome and lymph nodes (N)) were included in the nomogram analysis, which was aimed to predict cervical cancer prognosis. The C-index value was 0.760 (95% CI, 0.698-0.822) for OS, 0.703 (95% CI, 0.645-0.760) for PFI, and 0.802 (95% CI, 0.751-0.852) for DSS (Figure S2A, C, E). Calibration curve analysis presented good agreement between fit and actual observation for 1-, 3-, and 5-year OS, PFI, and DSS (Figure S2B, D, F). These results suggested that DHCR7 expression level effectively predicted the OS, FPI, and DSS for cervical squamous cell carcinoma patients, with reliable performance. High DHCR7 levels may be a biomarker of poor prognosis for cervical cancer.

3.4. Identifying Differentially Expressed Genes in Cervical Cancer

Using the cut-off threshold of adj. and |log2FC| >1.0, 554 differentially expressed genes (DEGs) were identified between high-DHCR7 and low-DHCR7 cervical cancer patients. Among which, 261 genes were upregulated and 293 genes were downregulated (Figures 4(a) and 4(b)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 4

554 differentially expressed genes (DEGs) were identified between high-DHCR7 and low-DHCR7 expression cervical cancer patients. (a) Volcano plot of DEGs. (b) Heat map DEGs. (c–f) GO analysis and KEGG enrichment analysis of DEGs. (c) Enriched GO terms in the “Biological process” category. (d) Enriched GO terms in the “Cellular Component” category. (e) Enriched GO terms in the “Molecular Function” category. (f) Enriched KEGG pathways. (g and h) PPI network and MCODE components of downregulated DEGs. (i and j) PPI network and the MCODE components of upregulated DEGs. DEGs: differential expressed genes.

3.5. Functional Enrichment Analysis of DEGs

To estimate the functional significance of DHCR7 in cervical cancer, GO classification and KEGG enrichment analyses were conducted based on the 554 DEGs using clusterProfiler package. In the BP category, 148 enriched GO terms were identified. They were mainly enriched in T cell activation (GO: 0042110), regulation of T cell activation (GO: 0050863), regulation of lymphocyte activation (GO: 0051249), regulation of lymphocyte proliferation (GO: 0050670), and regulation of mononuclear cell proliferation (GO: 0032944) (Figure 4(c)). The results suggested there may be a link between aberrant DHCR7 expression and immunity. Results on the CC category revealed 18 enriched GO terms, including external side of plasma membrane (GO: 0009897), anchored component of membrane (GO: 0031225), intrinsic component of synaptic membrane (GO: 0099240), and integral component of synaptic membrane (GO: 0099699) (Figure 4(d)). MF category results showed 13 enriched GO terms associated with carbohydrate binding (GO: 0030246), receptor ligand activity (GO: 0048018), lipase activity (GO: 0016298), and channel activity (GO: 0015267) (Figure 4(e)). KEGG pathway analysis identified 15 enriched pathways, including Staphylococcus aureus infection (hsa05150), cell adhesion molecules (CAMs) (hsa04514), natural killer cell mediated cytotoxicity(hsa04650), and T cell receptor signaling pathway (hsa04660) (Figure 4(f)).

3.6. DHCR7 Related Signaling Pathways Based on GSEA

GSEA was used to screen the potential biological processes and signaling pathway enriched using the 554 DEGs. This analysis revealed significant differences (FDR <0.05, adjusted -value <0.05) in MSigDB Collection (c2.cp.v6.2.symbols).The most significantly enriched biological processes and signaling pathways were selected based on Normalized Enrichment Score (NES). DHCR7 was found to be correlated with WNT activated receptor activity, PCG protein complex, G2/M checkpoint, E2/F targets, mTORC1 signaling, KRAS signaling, regulation of cholesterol biosynthetic, fibroblast growth factor receptor signaling pathway, T-cell receptor signaling pathway, JAK/STAT signaling pathway, T cell activation, and macrophage migration (Figure 5).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

3.7. The Correlation between DHCR7 Expression and Immune Infiltration

To identify the correlation between immune infiltration levels of tumor-infiltrating lymphocytes (TILs) and DHCR7 expression in the cervical cancer, Spearman’s correlation analysis was applied to 304 cervical cancer samples. Immune cell infiltration levels were quantified using ssGSEA. The result showed that DHCR7 expression was negatively correlated with cytotoxic cells, T cell, and CD8+T cell ( = -0.436, -0.413, and -0.384, respectively (, Figures 6(b)–6(d))). The Wilcoxon rank sum test indicated that infiltrating level of cytotoxic cells and T-cell was significantly higher in DHCR7-low samples in comparison with DHCR7-high ones (, Figures 6(e)–6(h)). These data showed that DHCR7 expression was correlated with reduced infiltration of cytotoxic cells, CD8+T cell, and T cell.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

3.8. Protein-Protein Interaction (PPI) Network Construction and Hub Gene Selection

PPI analysis of cervical cancer DEGs was conducted using Metascape. The PPI network and MCODE components identified using the DEGs in the gene lists are shown on Figures 4(g)–4(j). The top 3 functional clusters of modules in the downregulated DEGs and 4 MCODE components in the upregulated DEGs detected from the PPI network by MCODE were involved in G alpha, T cell receptor, and CD8 TCR signaling pathways.

3.9. DHCR7 Expression Was Increased in Cervical Cancer Cell Lines and Downregulation of DHCR7 Suppressed Cell Proliferation

Western blot was performed to validate DHCR7 expression in cancer cell lines. Compared with non-tumor epithelial cells, DHCR7 expression level in Hela and C33-A cervical cancer cell lines was increased obviously (Figure 7(a)). We used siRNA of DHCR7 to knock down DHCR7 expression in Hela and C33-A (Figure 7(b)) and detected cell proliferation. After downregulation of DHCR7, cell proliferation was significantly suppressed in Hela and C33-A cells (; Figures 7(c) and 7(d)).

(a)

(b)

(c)

(d)

4. Discussion

DHCR7 (7-dehydrocholesterol reductase) is the terminal enzyme of the Kandutsch-Russell (MK–R) and Bloch pathways of cholesterol synthesis [22]. The MK–R pathway generates 7DHC for vitamin D3 synthesis [23]. DHCR7 and EBP serve as regulatory and catalytic subunits of antiestrogen binding site (AEBS), which is a target of tamoxifen [24]. DHCR7 inhibitors and 7-DHC are considered antiviral therapies against emerging or highly pathogenic viruses [25]. Gerrick et al. [26] found that DHCR7 was a key regulator to M2 macrophage polarization and suppresses IL-10 and TARC secretion. These findings highlight the cholesterol pathway as a potential target for macrophage reprogramming therapies.

We found that DHCR7 was differentially expressed between cancer and normal tissues. Elevated DHCR7 expression in cervical carcinoma was correlated with poor clinicopathological features including advanced clinical stage, advanced T stage, and lymph node invasion. ROC curve indicated that DHCR7 had a strong diagnostic power for cervical cancer. Patients with elevated DHCR7 levels had poor OS, DSS, and PFI. Though there were no OS or DSS differences between patients with different DHCR7 expression levels in adenosquamous subgroup, it should be noticed that only 52 patients were included in this subgroup. Multivariate analysis indicated that DHCR7 was an independent prognostic biomarker for clinical outcome. Increased DHCR7 predicted poor prognosis. Taken together, the current study highlighted the potential role of DHCR7 as a prognostic biomarker for cervical cancer.

For a more accurate prognosis prediction, nomograms combining DHCR7 with other clinical features were developed to predict the 1-, 3-, and 5-year of OS, PFI, and DSS in cervical cancer. This analysis revealed C-index values of 0.760, 0.703, and 0.803 for OS, PFI, and DSS, respectively (Figure S2), indicating reliable predictive performance.

It should be noticed that we also used the GEO cohort (GSE30760) to validate whether DHCR7 expression was associated with survival of cervical cancer. However, the result showed that survival of patients in the high DHCR7 and low DHCR7 expression groups was not significantly different (Figure S3). One reason may be that the sample size was very small. Only 48 patients were included in the cohort. Another reason may be that patient characteristic in GSE30760 was different with patient characteristics in TCGA database. There were 6.9% patients in TCGA database while no patients in GSE30760 were diagnosed with stage IV disease.

To investigate the potential biological function of DHCR7 and other co-expression genes, we performed functional enrichment and GSEA analysis. It was shown that DHCR7 was significantly involved in multiple signaling pathways, including G2/M checkpoint, mTORC1 signaling, KRAS signaling, WNT signaling, PCG protein complex, and regulation of cholesterol biosynthetic. The above pathways have also been previously implicated in other cancers. The results suggested that DHCR7 may promote tumor development and progression via these pathways in cervical cancer.

mTOR signaling plays critical roles in the tumor microenvironment (TME) and tumor angiogenesis [27]. Inhibition of tumor angiogenesis is an effective strategy for limiting tumor growth and preventing tumor metastasis [28]. As it is shown by other research, G2/M checkpoint abrogation prevents cancer cells from repairing DNA damage, especially radiation-induced DNA damage, which has emerged as an anticancer target [29]. In cancer cells, DNA damage is triggered by radioactive rays and chemical toxins and causes G2/M cell cycle arrest. Similarly, cancer cells become sensitive to radiotherapy after the abrogation of the radiation-induced G2 arrest [29, 30]. DHCR7 may modulate G2/M checkpoint in response to DNA damage. Interestingly, we found that DHCR7-high patients undergoing radiotherapy had worse prognosis. DHCR7 is a potential anticancer target to improve radiosensitivity and chemosensitivity in cervical cancer. However, this hypothesis needs experimental validation.

Immune infiltration and tumor microenvironment have emerged as modulators of cancer progression and may affect clinical outcomes. In this study, we used ssGSEA and Spearman’s correlation analysis to unveil connections between DHCR7 expression and immune infiltration level. Our data demonstrated that DHCR7 was negatively correlated with the infiltration level of T cells, CD8+ T cells, and cytotoxic cells. Recent studies showed that the presence of tumor-infiltrating lymphocytes (TILs), particularly CD8+T cells, indicates favorable prognosis in several solid tumor types [31]. Patients with large T or B- cells infiltration or with T-cell or B-cell gene activation signature exhibit better survival [32, 33]. Additionally, primary tumor growth and metastasis are associated with decreased intratumoral immune T-cell densities [34]. High DHCR7 expression may suppress the function of immune system in cervical cancer by reducing CTLs and T-cell infiltration.

Biological progresses or pathways like T-cell receptor signaling, JAK/STAT, T cell activation, and macrophage migration were also detected in our study. T-cells are reported to have critical anti-tumor responses and have exhibited anticancer efficacy [35, 36]. Huang W et al. found that TCR-based immunotherapies may reduce solid tumor viability, including immune-checkpoint inhibitor refractory cancers [37]. Additionally, the PPI network analysis revealed that DHCR7 may participate in G alpha signaling events, T cell receptor signaling pathway, and the CD8 TCR pathway in cancer. Thus, we speculate that DHCR7 may affect cervical cancer development and progression by regulating the immune-related pathways.

Some limitations of our study should be noticed. Firstly, the study was based on bioinformatic analysis and without clinical validation. Secondly, it was a single-gene analysis. Thus, the analysis was not comprehensive. What’s more, the size of control samples was small. Experimental studies using clinical specimens, in vitro and in vivo approaches are warranted.

5. Conclusion

In summary, we found that DHCR7 was significantly upregulated in cervical cancers and was correlated with unfavorable outcomes. DHCR7 may be a potential biomarker for poor survival in cervical cancer.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

JZ, BY, SL, and JL had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. JZ, BY, SL, and JL designed the study. JZ, BY, and JL contributed to the acquisition of data. JZ, BY, SL, and JL contributed to the analysis and interpretation of data. BY, SL, and JL contributed to the statistical analysis. JZ, BY, SL, and JL contributed to the drafting and revising of the manuscript.

Supplementary Materials

Supplementary 1. Figure S1 showed the Kaplan-Meier survival curves of patients with high and low DHCR7 expression in the adenosquamous cervical cancer subgroup.

Supplementary 2. Figure S2 showed the nomogram and calibration curves of DHCR7.

Supplementary 3. Figure S3. The GEO cohort (GSE30760) was used to validate whether DHCR7 expression was associated with survival of cervical cancer. However, the result showed that survival of patients in the high DHCR7 and low DHCR7 expression groups was not significantly different.

References

R. L. Siegel, K. D. Miller, H. E. Fuchs, and A. Jemal, “Cancer statistics, 2022,” CA: a Cancer Journal for Clinicians, vol. 72, no. 1, pp. 7–33, 2022.
View at: Publisher Site | Google Scholar
W. J. Koh, N. R. Abu-Rustum, S. Bean et al., “Cervical cancer, version 3.2019, NCCN clinical practice guidelines in oncology,” Journal of the National Comprehensive Cancer Network, vol. 17, no. 1, pp. 64–84, 2019.
View at: Publisher Site | Google Scholar
J. Y. Kim, S. J. Byun, Y. S. Kim, and J. H. Nam, “Disease courses in patients with residual tumor following concurrent chemoradiotherapy for locally advanced cervical cancer,” Gynecologic Oncology, vol. 144, no. 1, pp. 34–39, 2017.
View at: Publisher Site | Google Scholar
X. Li, R. Roberti, and G. Blobel, “Structure of an integral membrane sterol reductase from Methylomicrobium alcaliphilum,” Nature, vol. 517, no. 7532, pp. 104–107, 2015.
View at: Publisher Site | Google Scholar
F. Lin, R. Zheng, C. Yu, Y. Su, X. Yan, and F. Qu, “Predictive role of serum cholesterol and triglycerides in cervical cancer survival,” International Journal of Gynecological Cancer, vol. 31, no. 2, pp. 171–176, 2021.
View at: Publisher Site | Google Scholar
Y. Peng, R. Myers, W. Zhang, and E. Alexov, “Computational investigation of the missense mutations in DHCR7 gene associated with Smith-Lemli-Opitz syndrome,” International Journal of Molecular Sciences, vol. 19, no. 1, p. 141, 2018.
View at: Publisher Site | Google Scholar
E. Theodoratou, T. Palmer, L. Zgaga et al., “Instrumental variable estimation of the causal effect of plasma 25-hydroxy-vitamin D on colorectal cancer risk: a mendelian randomization analysis,” PLoS One, vol. 7, no. 6, article e37662, 2012.
View at: Publisher Site | Google Scholar
C. Carlberg and A. Munoz, “An update on vitamin D signaling and cancer,” Seminars in Cancer Biology, vol. 79, pp. 217–230, 2020.
View at: Publisher Site | Google Scholar
W. B. Grant, “Review of recent advances in understanding the role of vitamin D in reducing cancer risk: breast, colorectal, prostate, and overall cancer,” Anticancer Research, vol. 40, no. 1, pp. 491–499, 2020.
View at: Publisher Site | Google Scholar
I. S. Carvalho, C. I. Gonçalves, J. T. Almeida et al., “Association of vitamin D pathway genetic variation and thyroid cancer,” Genes (Basel), vol. 10, no. 8, p. 572, 2019.
View at: Publisher Site | Google Scholar
J. Vivian, A. A. Rao, F. A. Nothaft et al., “Toil enables reproducible, open source, big biomedical data analyses,” Nature Biotechnology, vol. 35, no. 4, pp. 314–316, 2017.
View at: Publisher Site | Google Scholar
M. I. Love, W. Huber, and S. Anders, “Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2,” Genome Biology, vol. 15, no. 12, p. 550, 2014.
View at: Publisher Site | Google Scholar
C. Gene Ontology, “Gene Ontology Consortium: going forward,” Nucleic Acids Research, vol. 43, no. D1, pp. D1049–D1056, 2015.
View at: Publisher Site | Google Scholar
M. Kanehisa, Y. Sato, M. Furumichi, K. Morishima, and M. Tanabe, “New approach for understanding genome variations in KEGG,” Nucleic Acids Research, vol. 47, no. D1, pp. D590–D595, 2019.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 43, pp. 15545–15550, 2005.
View at: Publisher Site | Google Scholar
G. Bindea, B. Mlecnik, M. Tosolini et al., “Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer,” Immunity, vol. 39, no. 4, pp. 782–795, 2013.
View at: Publisher Site | Google Scholar
D. A. Barbie, P. Tamayo, J. S. Boehm et al., “Systematic RNA interference reveals that oncogenic KRAS -driven cancers require TBK1,” Nature, vol. 462, no. 7269, pp. 108–112, 2009.
View at: Publisher Site | Google Scholar
Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, p. 1523, 2019.
View at: Publisher Site | Google Scholar
G. D. Bader and C. W. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, p. 2, 2003.
View at: Publisher Site | Google Scholar
X. Robin, N. Turck, A. Hainard et al., “pROC: an open-source package for R and S+ to analyze and compare ROC curves,” BMC Bioinformatics, vol. 12, no. 1, p. 77, 2011.
View at: Publisher Site | Google Scholar
A. V. Prabhu, W. Luu, L. J. Sharpe, and A. J. Brown, “Phosphorylation regulates activity of 7-dehydrocholesterol reductase (DHCR7), a terminal enzyme of cholesterol synthesis,” The Journal of Steroid Biochemistry and Molecular Biology, vol. 165, no. Part B, pp. 363–368, 2017.
View at: Publisher Site | Google Scholar
M. A. Mitsche, J. G. McDonald, H. H. Hobbs, and J. C. Cohen, “Flux analysis of cholesterol biosynthesis in vivo reveals multiple tissue and cell-type specific pathways,” eLife, vol. 4, article e07999, 2015.
View at: Publisher Site | Google Scholar
R. L. Sutherland, L. C. Murphy, M. S. Foo, M. D. Green, A. M. Whybourne, and Z. S. Krozowski, “High-affinity anti-oestrogen binding site distinct from the oestrogen receptor,” Nature, vol. 288, no. 5788, pp. 273–275, 1980.
View at: Publisher Site | Google Scholar
J. Xiao, W. Li, X. Zheng et al., “Targeting 7-dehydrocholesterol reductase integrates cholesterol metabolism and IRF3 activation to eliminate infection,” Immunity, vol. 52, no. 1, pp. 109–122.e6, 2020.
View at: Publisher Site | Google Scholar
K. Y. Gerrick, E. R. Gerrick, A. Gupta, S. J. Wheelan, S. Yegnasubramanian, and E. M. Jaffee, “Transcriptional profiling identifies novel regulators of macrophage polarization,” PLoS One, vol. 13, no. 12, article e0208602, 2018.
View at: Publisher Site | Google Scholar
L. C. Kim, R. S. Cook, and J. Chen, “mTORC1 and mTORC2 in cancer and the tumor microenvironment,” Oncogene, vol. 36, no. 16, pp. 2191–2201, 2017.
View at: Publisher Site | Google Scholar
S. Q. Liang, E. D. Bührer, S. Berezowska et al., “mTOR mediates a mechanism of resistance to chemotherapy and defines a rational combination strategy to treat KRAS -mutant lung cancer,” Oncogene, vol. 38, no. 5, pp. 622–636, 2019.
View at: Publisher Site | Google Scholar
M. T. Dillon, J. S. Good, and K. J. Harrington, “Selective targeting of the G2/M cell cycle checkpoint to improve the therapeutic index of radiotherapy,” Clinical Oncology (Royal College of Radiologists), vol. 26, no. 5, pp. 257–265, 2014.
View at: Publisher Site | Google Scholar
N. Bucher and C. D. Britten, “G2 checkpoint abrogation and checkpoint kinase-1 targeting in the treatment of cancer,” British Journal of Cancer, vol. 98, no. 3, pp. 523–528, 2008.
View at: Publisher Site | Google Scholar
J. C. Kong, G. R. Guerra, T. Pham et al., “Prognostic impact of tumor-infiltrating lymphocytes in primary and metastatic colorectal cancer: a systematic review and meta-analysis,” Diseases of the Colon and Rectum, vol. 62, no. 4, pp. 498–508, 2019.
View at: Publisher Site | Google Scholar
M. Schmidt, B. Hellwig, S. Hammad et al., “A comprehensive analysis of human gene expression profiles identifies stromal immunoglobulin κ C as a compatible prognostic marker in human solid tumors,” Clinical Cancer Research, vol. 18, no. 9, pp. 2695–2703, 2012.
View at: Publisher Site | Google Scholar
J. S. Nielsen, R. A. Sahota, K. Milne et al., “CD20+ tumor-infiltrating lymphocytes have an atypical CD27- memory phenotype and together with CD8+ T cells promote favorable prognosis in ovarian cancer,” Clinical Cancer Research, vol. 18, no. 12, pp. 3281–3292, 2012.
View at: Publisher Site | Google Scholar
B. Mlecnik, M. Tosolini, A. Kirilovsky et al., “Histopathologic-based prognostic factors of colorectal cancers are associated with the state of the local immune reaction,” Journal of Clinical Oncology, vol. 29, no. 6, pp. 610–618, 2011.
View at: Publisher Site | Google Scholar
A. Ohno, T. Iwata, Y. Katoh et al., “Tumor-infiltrating lymphocytes predict survival outcomes in patients with cervical cancer treated with concurrent chemoradiotherapy,” Gynecologic Oncology, vol. 159, no. 2, pp. 329–334, 2020.
View at: Publisher Site | Google Scholar
D. E. Speiser, P. C. Ho, and G. Verdeil, “Regulatory circuits of T cell function in cancer,” Nature Reviews. Immunology, vol. 16, no. 10, pp. 599–611, 2016.
View at: Publisher Site | Google Scholar
W. Huang and A. August, “The signaling symphony: T cell receptor tunes cytokine-mediated T cell differentiation,” Journal of Leukocyte Biology, vol. 97, no. 3, pp. 477–485, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Jing Zou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies