Skip to main content
Download PDF

CACNA1C is a prognostic predictor for patients with ovarian cancer

Journal of Ovarian Research volume 14, Article number: 88 (2021) Cite this article

Abstract

Background

CACNA1C, as a type of voltage-dependent calcium ion transmembrane channel, played regulatory roles in the development and progress of multiple tumors. This study was aimed to analyze the roles of CACNA1C in ovarian cancer (OC) of overall survival (OS) and to explore its relationships with immunity.

Methods

Single gene mRNA sequencing data and corresponding clinical information were obtained from The Cancer Genome Atlas Database (TCGA) and the International Cancer Genome Consortium (ICGC) datasets. Gene set enrichment analysis (GSEA) was used to identify CACNA1C-related signal pathways. Univariate and multivariate Cox regression analyses were applied to evaluate independent prognostic factors. Besides, associations between CACNA1C and immunity were also explored.

Results

CACNA1C had a lower expression in OC tumor tissues than in normal tissues (P < 0.001), with significant OS (P = 0.013) and a low diagnostic efficiency. We further validated the expression levels of CACNA1C in OC by means of the ICGC dataset (P = 0.01), qRT-PCR results (P < 0.001) and the HPA database. Univariate and multivariate Cox hazard regression analyses indicated that CACNA1C could be an independent risk factor of OS for OC patients (both P < 0.001). Five significant CACNA1C-related signaling pathways were identified by means of GSEA. As for genetic alteration analysis, altered CACNA1C groups were significantly associated with OS (P = 0.0169), progression-free survival (P = 0.0404), disease-free survival (P = 0.0417) and disease-specific survival (P = 9.280e-3), compared with unaltered groups in OC. Besides, CACNA1C was dramatically associated with microsatellite instability (MSI) and immunity.

Conclusions

Our results shed light on that CACNA1C could be a prognostic predictor of OS in OC and it was closely related to immunity.

Introduction

Ovarian cancer (OC) is the third most common gynecologic malignancy worldwide, with approximately 21,750 newly estimated cases and 13,940 newly estimated death in USA, 2020 [1]. Despite other cancers such as endometrial cancer (EC) having higher rates of incidence, OC is the deadliest of all female reproductive cancers [2]. This disease is rare in young women, especially under the age of 30 and its risks will increase with age, sharply after the age of 50. So, the average diagnosis age of OC is between 50 and 70 [3]. Because the clinical manifestations of early OC are hidden or unspecific, OC is often called a silent killer, and approximately 70 % of them are diagnosed with an advanced stage [4]. As reported, the 5-year survival rate for FIGO stage I is 90 %, stage II is 65 %, stage III is 34 %, and stage IV is 15 % [5, 6]. However, current improvements in treatment methods could only slightly improve OC survival rates [7]. Therefore, there was an urgently need to clarify the potential molecular mechanisms of OC and identify novel biomarkers for early diagnosis and prognosis assessment [8, 9].

Calcium channels could be generally classified into two categories containing the voltage-gated calcium channels (VGCCs) and the ligand-gated calcium channels (LGCCs) [10]. As a type of VGCCs, calcium voltage-gated channel subunit alpha1C (CACNA1C) was coded by the α1 subunit Cav1.2 and it was reported to be involved in regulating cell-matrix adhesion, collagen fibril organization, cell adhesion, cellular response to amino acid stimulus, and negative regulation of cell proliferation [11,12,13]. Previous meta-analysis results showed that CACNA1C was up-regulated in brain tumors, leukemia, breast cancer and other tumors, suggesting its regulatory roles in cancer progression [14]. Bioinformatics analysis also revealed that mutation of CACNA1C was significantly associated with longer overall survival (OS) in EC patients [15]. All of these indicated that CACNA1C could be served as a prognostic biomarker for human cancers. Therefore, this study was aimed to explore the prognostic values of CACNA1C in OC, to identify CACNA1C-related signal pathways and to further explore its associations with immunity, thus helping to provide some references for future basic CACNA1C-related researches in OC.

Materials and methods

Data acquisition and processing

CACNA1C single gene expression matrix and corresponding clinical data of 379 OC tumor samples were obtained from The Cancer Genome Atlas (TCGA) database (http://cancergenome.nih.gov/), while 88 ovary normal tissue samples were getting from the Genetype-Tissue Expression (GTEx) data portal (https://gtexportal.org/). In addition, another OC dataset from the International Cancer Genome Consortium (ICGC) data portal (https://icgc.org/) was used as external verification database, involving 111 tumor samples. All data were preprocessed by normalization or log2 transformation and analyzed by R software (https://www.rproject.org/) [16]. The “limma” package was used to calculate differentially expressed genes (DEGs), with |log2 fold change (FC)| >1 and adjusted P-values (FDR) < 0.05 as the cut-off criteria.

Quantitative real-time PCR (qRT‑PCR)

The qRT‑PCR analysis was performed to detect the expression levels of CACNA1C in OC tissues, according to the manufacturer’s instructions (TAKARA) [17]. All reactions were carried out in triplicate and Actin was used as the internal controls for mRNA. The relative gene expression in tissues was calculated using the comparative delta-delta CT (2-ΔΔCt) method. Primers were synthesized by TSINGKE (Beijing, China), including Actin (F:5’ ATGACTTAGTTGCGTTACACC 3’, R:5’ GACTTCCTGTAACAACGCATC 3’) and CACNA1C (F:5’ GAAGCGGCAGCAATATGGGA 3’, R:5’ TTGGTGGCGTTGGAATCATCT 3’). Six paired tumor and adjacent normal tissues were obtained from OC patients from Shengjing Hospital of China Medical University. Ethical approval was obtained from the Institutional Research Ethics Committees of Shengjing Hospital of China Medical University and informed written consent was obtained from all subjects.

The Human Protein Atlas (HPA) database and protein-protein interaction (PPI) network analysis

The Human Protein Atlas (HPA, http://www.proteinatlas.org/) online database was utilized to validate the protein expression of CACNA1C in OC by immunohistochemical staining in this study. With the help of online STRING (https://string-db.org/) database, PPI full network was conducted to explore the potential relationships between CACNA1C and other related genes in OC, with the medium confidence > 0.40.

The receiver operating characteristic (ROC) curves and logistic regression analysis

By utilizing the R “survivalROC” package, ROC curves and the area under the curve (AUC) values were performed to assess the specificity and sensitivity of CACNA1C for OC. Moreover, logistic regression analysis was utilized to explore the associations between CACNA1C expression and 10 clinicopathological parameters (stage, grade, age, location, race, status, diagnosis, venous-invasion, chemotherapy, outcome) in OC patients and P-values below 0.05 were regarded as statistical significance [18].

Univariate and multivariate Cox hazard regression analysis

We performed univariate and multivariate Cox hazard regression analysis to determine whether or not CACNA1C and 10 clinicopathological parameters (stage, grade, age, location, race, status, diagnosis, venous-invasion, chemotherapy, outcome) could be independent factors related to OS on OC patients in the TCGA dataset by means of R package.

Gene set enrichment analysis (GSEA) analysis

GSEA, as a computational method determining whether or not a priori defined set of genes shows statistically significant, concordant differences between two biological states, was carried out to explore the significant survival differences between high- and low-CACNA1C groups [19]. The expression level of CACNA1C was used as a phenotype label and gene set permutations were performed 1000 times for each analysis with the consideration of the nominal p value < 0.05 and normalized enrichment score (NES) > 1.5 as the threshold.

Genetic alteration analysis

We utilized the cBioPortal for Cancer Genomics (http://cbioportal.org) to explore the CACNA1C alteration frequency, mutation type and CNA (copy number alteration) across all TCGA tumors by means of the “Cancer Types Summary” module. Furthermore, the overall survival, progression-free survival, disease-free survival, and disease-specific survival differences for OC with or without CACNA1C genetic alteration were also showed by Kaplan-Meier (K-M) survival plots and log-rank p-values [20, 21].

The evaluation of microsatellite instability (MSI), tumor mutational burden (TMB), tumor neoantigen burden (TNB) and immunity

In order to explore the associations between CACNA1C gene expression and MSI or TMB or TNB, correlation analyses was performed by the Spearman’s method and visualized by radar maps [22, 23]. As previously described, we further explored the relationships between CACNA1C gene expression and immunity. Four aspects were analyzed including tumor immune infiltration, tumor microenvironment, immune checkpoint molecules and immune cells pathway by R packages“estimate”, “limma”, “reshape2” and “RColorBrewer” [24, 25]. All of these above mentioned analyses were performed by using the Sangerbox tool, which is a free online data analysis platform (http://www.sangerbox.com/tool).

Results

The expression level of CACNA1C in OC based on TCGA and GTEx datasets

As showed in Fig. 1 A, we obtained the expression level of CACNA1C in various normal tissues and organs by GTEx database and noticed that there were significant differences between female and male in adipose tissues, breast tissues, and heart tissues (all P < 0.01). Anatomical heatmap of CACNA1C expression in female indicated that it was lowly expressed in ovary (Fig. 1B). We combined TCGA and GTEx datasets to investigate the different expression levels of CACNA1C between pan-cancers and adjacent normal tissues (Fig. 1 C). As displayed in Fig. 1D, the expression level of CACNA1C in OC tissues was obviously lower than that in normal tissues (P < 0.001). Based on the median expression as the cut-off value, OC patients were classified into low- and high-CACNA1C groups and significant expression differences were displayed in these two groups (P < 0.0001). K-M survival analysis presented that high-CACNA1C groups had a shorted OS than those in the low-CACNA1C groups (P = 0.013; Fig. 1E). ROC curves were carried out to demonstrate the performance of the CACNA1C expression level for survival prediction and its AUC values for 1-, 3-, and 5-year survival were 0.622, 0.609, 0.572, merely having a low diagnostic efficiency (Fig. 1 F).

Fig. 1
figure 1

The expression level of CACNA1C in OC based on TCGA and GTEx datasets; (A) CACNA1C expression levels in various normal tissues and organs by GTEx database; (B) Anatomical heatmap of CACNA1C expression in female by GTEx database; (C) CACNA1C mRNA expression levels in pan-cancers from TCGA and GTEx datasets; (D) Boxplot of CACNA1C expression between the OC tumor and normal tissues from TCGA and GTEx datasets (Normal = 88 and Tumor = 379); (E) K-M survival analysis of CACNA1C in TCGA; (F) ROC curves and its AUCs for 1-, 3-, and 5-year survival of CACNA1C; **P < 0.01; ***P < 0.001

Full size image

Validation of CACNA1C expression in OC

We utilized the ICGC data portal (https://icgc.org/) as external verification database. Based on the median expression, OC patients in ICGC dataset were classified into low- and high-risk subgroups and high-CACNA1C groups also had a lower OS than those in the low-CACNA1C groups (P = 0.01, Fig. 2 A). Moreover, qRT-PCR results verified that the expression levels of CACNA1C were significantly down-regulated in the OC tissues compared with normal tissues (Normal = 6; Tumor = 6; P < 0.001; Fig. 2B). Besides, immunohistochemical pictures from the HPA database (https://www.proteinatlas.org/) showed that CACNA1C utilizing HPA039796 antibody had a medium expression in 3/3 normal ovary tissues, while its expression was not detected in 10/12, low in 1/12, medium in 1/12 OC tumor tissues. Obviously, CACNA1C was also lowly expressed in OC tumor tissues compared with normal tissues (Fig. 2 C-D).

Fig. 2
figure 2

Validation of CACNA1C expression in OC; (A) K-M survival analysis of CACNA1C in OC from ICGC dataset; (B) qRT-PCR verification of CACNA1C expression in normal and OC tissues (Normal = 6; Tumor = 6); (C-D) Immunohistochemical staining from the HPA database for CACNA1C in OC; ***P < 0.001

Full size image

Associations between CACNA1C expression and clinicopathologic characteristics in OC

Logistic regression analysis results indicated that P value(s) for the associations between CACNA1C expression and age was 0.43; chemotherapy was 0.64; grade was 0.59; stage was 0.35; status was 0.27; venous-invasion was 0.96; location were all above 0.05; race were all above 0.05; diagnosis were all above 0.05; outcome were all above 0.05 (Fig. 3). Obviously, all P values were above 0.05. In other words, there were no significant associations between the CACNA1C expression and stage, grade, age, location, race, status, diagnosis, venous-invasion, chemotherapy or outcome in TCGA OC samples.

Fig. 3
figure 3

Associations between CACNA1C expression and clinicopathologic characteristics; (A) Age; (B) Chemotherapy; (C) Grade; (D) Stage; (E) Status; (F) Venous-invasion; (G) Location; (H) Race; (I) Diagnosis; (J) Outcome

Full size image

CACNA1C could serve as an independent prognostic factor

To identify independent prognostic factors, univariate and multivariate Cox hazard regression analyses were utilized to exclude CACNA1C and 10 clinicopathological parameters (stage, grade, age, location, race, status, diagnosis, venous-invasion, chemotherapy, outcome) with little OS values. According to our results, the CACNA1C expression was significantly associated with OS in both univariate and multivariate Cox hazard regression analyses (both P < 0.001), indicating that CACNA1C could be an independent risk factor of OS for OC prognosis (Fig. 4; Table 1).

Fig. 4
figure 4

CACNA1C could serve as an independent prognostic factor of OS in ccRCC; (A) Univariate Cox hazard regression analysis of CACNA1C and clinicopathologic variables of ccRCC in TCGA cohort; (B) Multivariate Cox hazard regression analysis of CACNA1C and clinicopathologic variables of ccRCC in TCGA cohort

Full size image
Table 1 Univariate and multivariate cox analysis of CACNA1C and clinicopathologic characteristics of overall survival in TCGA OV cohort
Full size table

GSEA identified CACNA1C-related signaling pathways

In order to explore how CACNA1C participates in the pathogenesis of OC, GSEA analysis was conducted to select related signaling pathways. GSEA analysis was performed between the high- and low-CACNA1C expression datasets, with the consideration of the nominal p value < 0.05 and NES > 1.5 as the threshold. Finally, five pathways that exhibited significant differential enrichment in the high-CACNA1C expression phenotype were identified, containing Base excision repair signaling pathway, Glyoxylate and dicarboxylate metabolism signaling pathway, Mismatch repair signaling pathway, Nucleotide excision repair signaling pathway and Oxidative phosphorylation signaling pathway (Fig. 5; Table 2), thus helping to understand the pathogenesis of OC.

Fig. 5
figure 5

GSEA identified CACNA1C-related signaling pathways; (A) Base excision repair signaling pathway; (B) Glyoxylate and dicarboxylate metabolism signaling pathway; (C) Mismatch repair signaling pathway; (D) Nucleotide excision repair signaling pathway; (E) Oxidative phosphorylation signaling pathway; (F) All of the five significant signaling pathways

Full size image
Table 2 Gene sets enriched in phenotype high
Full size table

Genetic alteration analysis of CACNA1C in OC by cBioPortal database

As displayed in Fig. 6 A, we noticed that the genetic alteration frequency of CACNA1C was about 10 % in OC by means of the cBioPortal tool from TCGA dataset. Moreover, CACNA1C mutation sites in OC were also detailed in Fig. 6B. Further K-M survival analyses indicated that altered CACNA1C groups were significantly associated with OS (P = 0.0169), progression-free survival (P = 0.0404), disease-free survival (P = 0.0417) and disease-specific survival (P = 9.280e-3), compared with unaltered groups in OC (Fig. 6 C-F). Besides, altered CACNA1C-related signaling pathways were also identified in Figure S1, including Cell-cycle-signaling-pathway, RTK-RAS-PI(3)K-pathway and Notch-signaling-pathway. All of these indicated that genetic alteration of CACNA1C might play vital roles in OC.

Fig. 6
figure 6

Genetic alteration analysis of CACNA1C in OC by cBioPortal database; (A) The alteration frequency with mutation type of CACNA1C in different tumor samples from TCGA cohorts; (B) CACNA1C mutation sites in OC from TCGA cohort; (C) K-M survival analysis of OS with CACNA1C altered or unaltered groups; (D) K-M survival analysis of progression-free survival with CACNA1C altered or unaltered groups; (E) K-M survival analysis of disease-free survival with CACNA1C altered or unaltered groups; (F) K-M survival analysis of disease-specific survival with CACNA1C altered or unaltered groups

Full size image

Associations between the CACNA1C expression and PPI, MSI, TMB, TNB in OC

PPI network indicated that ten genes (CALM1、CALM2、CALM3、CACNA2D1、CACNA2D2、CACNA2D4、CACNB1、CACNB2、CACNB3 and CACNB4) were significantly associated with the CACNA1C expression (Fig. 7 A). Based on OC samples from the TCGA cohort, we explored the associations between the CACNA1C expression and MSI, TMB, TNB. Our results shed light on that the CACNA1C expression was significantly related to MSI (P = 0.022) in OC, while it had nothing to do with TMB (P = 0.75) or TNB (P = 0.29) (Fig. 7B and D).

Fig. 7
figure 7

Associations between the CACNA1C expression and PPI, MSI, TMB, TNB in OC; (A) PPI network; (B) Associations between the CACNA1C expression and MSI; (C) Associations between the CACNA1C expression and TNB; (D) Associations between the CACNA1C expression and TMB

Full size image

Associations between CACNA1C and tumor immune infiltrations, tumor microenvironment, immune checkpoint molecules, immune cells pathway in OC

To further explore the relationships between the CACNA1C expression and immunity, four aspects were analyzed respectively, including tumor immune infiltration, tumor microenvironment, immune checkpoint molecules and immune cells pathway. We assessed the correlations between CACNA1C and the infiltration levels of six immune cells through online analysis of TIMER, and found that the CACNA1C expression was significantly associated with B cell infiltration, CD4 + T cell infiltration, neutrophil infiltration, macrophage infiltration, and dendritic cell infiltration (all P < 0.0; Fig. 8 A). As for tumor microenvironment, the CACNA1C expression was markedly related to ESTIMATEScore and StromalScore (both P < 0.001; Fig. 8B). Correlations between the CACNA1C expression and immune checkpoint molecules or immune cells pathway indicated that this gene was remarkably linked to ADORA2A, CD276, CD28, NRP1, TNFRSF8, TNFSF4, Central memory CD8 T cells pathway and Type 17 T helper cells pathway (Fig. 8 C-D).

Fig. 8
figure 8

Associations between CACNA1C and tumor immune infiltrations, tumor microenvironment, immune checkpoint molecules, immune cells pathway in OC; (A) Associations between the CACNA1C expression and tumor immune infiltration; (B) Associations between the CACNA1C expression and immune microenvironment; (C) Correlations between the CACNA1C expression and immune checkpoint molecules; (D) Correlations between the CACNA1C expression and immune cells pathway

Full size image

Discussion

As a common gynecological cancer, OC was a highly malignant disease with a high recurrence rate and a poor 5-year survival rate [26]. Although various treatments such as chemotherapy, radiotherapy, immunotherapy, and targeted therapy had been applied, the 5-year survival rate of advanced OC patients was still less than 30 % [27]. At present, it was difficult to diagnose early and there were also no accurate strategies to predict the prognosis and recurrence of OC [28]. Moreover, the classical TNM staging system, relying mainly on anatomical information without any molecular biological features, had difficulties in accurately stratifying indolent and aggressive cancers to a certain extent [29, 30]. Therefore, it was urgent to identify novel and effective biomarkers for OC prognosis and to further explore its related mechanisms [31].

In this study, we explored the prognostic values of CACNA1C in OC, identified CACNA1C-related signal pathways and further explored its associations with immunity. Our results indicated that CACNA1C had a lower expression in OC tumor tissues than in normal tissues, with significant OS and a low diagnostic efficiency. We further validated the expression levels of CACNA1C in OC by means of the ICGC dataset, qRT-PCR results and the HPA database. Furthermore, univariate and multivariate Cox hazard regression analyses indicated that CACNA1C could be an independent risk factor of OS affecting the prognosis of OC patients. In line with previously published articles, CACNA1C was up-regulated in brain tumors, leukemia, breast cancer and other tumors [14]. Besides, variants of CACNA1C were associated with bladder cancer, endometrial cancer and breast cancer risk [15, 32, 33]. All of these indicated that CACNA1C could also be a prognostic predictor for patients with OC.

GSEA, as a useful tool, had been widely applied in OC to identify related signaling pathways [34, 35]. Shi et al. found that knocking down CerS6 could significantly affect cell cycle in serous OC cells by GSEA and RNA sequencing analyses [36]. Zhang et al. performed a GSEA of SUCNR1 based on TCGA OC dataset and identified immune-related pathways [37]. In this article, GSEA analysis was performed between the high- and low-CACNA1C expression datasets to explore how CACNA1C participates in the pathogenesis of OC. Finally, five pathways that exhibited significant differential enrichment in the high-CACNA1C expression phenotype were detected, containing Base excision repair signaling pathway, Glyoxylate and dicarboxylate metabolism signaling pathway, Mismatch repair signaling pathway, Nucleotide excision repair signaling pathway and Oxidative phosphorylation signaling pathway.

In terms of genetic alteration, it had been found to be associated with the development and progress of various cancers [38,39,40]. Shen et al. revealed that genetic and epigenetic alterations of p33ING1b were related to the pathogenesis of OC [41]. Pongstaporn et al. found that the GSTO2 gene polymorphism might be easily associated with OC risks [42]. Our results shed light on that the genetic alteration frequency of CACNA1C was about 10 % in OC by means of the cBioPortal tool from TCGA dataset. Moreover, K-M survival analyses showed that altered CACNA1C groups were significantly associated with OS, progression-free survival, disease-free survival and disease-specific survival, compared with unaltered groups in OC. All of these indicated that genetic alterations of CACNA1C might play vital roles in OC. Furthermore, we also identified altered CACNA1C-related signaling pathways, including Cell-cycle-signaling-pathway, RTK-RAS-PI(3)K-pathway and Notch-signaling-pathway.

A growing number of studies reported that MSI, TMB and TNB played important roles in tumorigenesis and they could serve as important predictive biomarkers for the application of immunotherapy [43,44,45]. Fan et al. developed and verified a TMB-related prognostic signature, which could be a promising prognostic signature for OC prognosis [46]. Xiao et al. suggested that mismatch repair (MMR) deficient OC patients might be good candidates for anti-PD-1/PD-L1 therapy, characterized by increased tumor-infiltrating lymphocytes and MSI phenotype [47]. Based on OC samples from the TCGA cohort, we also explored the associations between the CACNA1C expression and MSI, TMB, TNB. Our results shed light on that the CACNA1C expression was significantly related to MSI in OC, while it had nothing to do with TMB or TNB.

As for immunity, four aspects were analyzed respectively, including tumor immune infiltration, tumor microenvironment, immune checkpoint molecules and immune cells pathway. Previous articles reported that tumor immune infiltration and tumor microenvironment might be involved in OC tumorigenesis and response to immunotherapy [48, 49]. Our results indicated that the CACNA1C expression was significantly associated with B cell infiltration, CD4 + T cell infiltration, neutrophil infiltration, macrophage infiltration, and dendritic cell infiltration for tumor immune infiltration and markedly related to ESTIMATEScore and StromalScore for tumor microenvironment. Correlations between gene expression and immune checkpoint molecules or immune cells pathway were also analyzed and various genes had been identified in OC [50, 51]. Our results indicated that the CACNA1C expression was remarkably linked to ADORA2A, CD276, CD28, NRP1, TNFRSF8, TNFSF4, Central memory CD8 T cells pathway and Type 17 T helper cells pathway.

As for its potential mechanism, our results indicated that CACNA1C could be a prognostic predictor of OS in OC and it was closely related to immunity, having an effect on oxidative phosphorylation and DNA repair activities including base excision repair, mismatch repair and nucleotide excision repair signaling pathways. As reported by previous article, oxidative phosphorylation could serve as a target for novel therapeutic strategies against OC [52]. DNA repair activities were also revealed to be the most common cause of hereditary OC [53,54,55]. Further experimental researches were required to verify these findings.

The strength of this article was that it was the first time for us to explore the roles of CACNA1C in OC and we further validated this gene in the ICGC dataset, qRT-PCR results and the HPA database. This study had several limitations too. Firstly, due to the absence of normal OC samples in TCGA, we had to analyze CACNA1C single gene data in combination of TCGA and GTEx datasets. This inevitably would lead to some biases. Secondly, our study was a retrospective study based on public databases without further experimental verification from our own samples, except for qRT-PCR. Last but not least, a larger sample size of patients was required to verify the prognostic roles of CACNA1C in OC and to confirm its associations with immunity.

Conclusions

Taken together, our results shed light on that CACNA1C could be a prognostic predictor for patients with OC and five signaling pathways regulated by CACNA1C were also identified. Moreover, our analyses of CACNA1C indicated statistical correlations of CACNA1C expression with clinical prognosis, genetic alteration, MSI, tumor immune infiltration, tumor microenvironment, immune checkpoint molecules and immune cells pathway, helping to understand its role in OC from the perspective of clinical tumor samples. Further experimental validations were still required to mine the molecular mechanisms associated with CACNA1C in OC.

Availability of data and materials

The RNA-sequencing data and corresponding clinical information were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) and the International Cancer Genome Consortium (ICGC) database (https://icgc.org/).

References

  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.

    Article  PubMed  Google Scholar 

  2. Kuroki L, Guntupalli SR. Treatment of epithelial ovarian cancer. Bmj. 2020;371:m3773.

    Article  PubMed  Google Scholar 

  3. Menon U, Griffin M, Gentry-Maharaj A. Ovarian cancer screening–current status, future directions. Gynecol Oncol. 2014;132(2):490–5.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Liu HD, et al. Organoid of ovarian cancer: genomic analysis and drug screening. Clin Transl Oncol. 2020;22(8):1240–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Nash Z, Menon U. Ovarian cancer screening: Current status and future directions. Best Pract Res Clin Obstet Gynaecol. 2020;65:32–45.

    Article  PubMed  Google Scholar 

  6. Gogineni V, et al. Current Ovarian Cancer Maintenance Strategies and Promising New Developments. J Cancer. 2021;12(1):38–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Irusta G. Roads to the strategic targeting of ovarian cancer treatment. Reproduction. 2021;161(1):R1-r11.

    Article  PubMed  Google Scholar 

  8. Foo T, George A, Banerjee S. PARP inhibitors in ovarian cancer: An overview of the practice-changing trials. Genes Chromosomes Cancer; 2020.

  9. Tyagi K, Roy A. Evaluating the current status of protein kinase C (PKC)-protein kinase D (PKD) signalling axis as a novel therapeutic target in ovarian cancer. Biochim Biophys Acta Rev Cancer. 2020;1875(1):188496.

    Article  PubMed  CAS  Google Scholar 

  10. Phan NN, et al. Voltage-gated calcium channels: Novel targets for cancer therapy. Oncol Lett. 2017;14(2):2059–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Nie F, et al. Genetic analysis of SNPs in CACNA1C and ANK3 gene with schizophrenia: A comprehensive meta-analysis. Am J Med Genet B Neuropsychiatr Genet. 2015;168(8):637–48.

    Article  CAS  PubMed  Google Scholar 

  12. Endres D, et al. New Cav1.2 channelopathy with high-functioning autism, affective disorder, severe dental enamel defects, a short QT interval, and a novel CACNA1C loss-of-function mutation. Int J Mol Sci. 2020;21(22):8611.

  13. Lazary J, et al. Genetic analyses of the endocannabinoid pathway in association with affective phenotypic variants. Neurosci Lett. 2021;744:135600.

  14. Wang CY, et al. Meta-Analysis of Public Microarray Datasets Reveals Voltage-Gated Calcium Gene Signatures in Clinical Cancer Patients. PLoS One. 2015;10(7):e0125766.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Qiao Z, et al. Mutations in KIAA1109, CACNA1C, BSN, AKAP13, CELSR2, and HELZ2 Are Associated With the Prognosis in Endometrial Cancer. Front Genet. 2019;10:909.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhang N, et al. A six-gene-based prognostic model predicts complete remission and overall survival in childhood acute myeloid leukemia. Onco Targets Ther. 2019;12:6591–604.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Zhao X, et al. Linc00511 acts as a competing endogenous RNA to regulate VEGFA expression through sponging hsa-miR-29b-3p in pancreatic ductal adenocarcinoma. J Cell Mol Med. 2018;22(1):655–67.

    Article  CAS  PubMed  Google Scholar 

  18. Chen YN, et al. EZH2 is a potential prognostic predictor of glioma. J Cell Mol Med. 2021;25(2):925–36.

    Article  CAS  PubMed  Google Scholar 

  19. Subramanian A, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Cerami E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.

    Article  PubMed  Google Scholar 

  21. Gao J, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Miao Y, et al., Prognostic value and immunological role of PDCD1 gene in pan-cancer. Int Immunopharmacol, 2020. 89(Pt B): p. 107080.

  23. Hundal J, et al. pVAC-Seq: A genome-guided in silico approach to identifying tumor neoantigens. Genome Med. 2016;8(1):11.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612.

    Article  PubMed  CAS  Google Scholar 

  25. Newman AM, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li Y, et al. Prognostic Significance and Related Mechanisms of Hexokinase 1 in Ovarian Cancer. Onco Targets Ther. 2020;13:11583–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. DeSantis CE, et al. Cancer statistics for African Americans, 2019. CA Cancer J Clin. 2019;69(3):211–33.

    Article  PubMed  Google Scholar 

  28. Eisenhauer EA. Real-world evidence in the treatment of ovarian cancer. Ann Oncol. 2017;28(suppl_8):viii61–5.

    Article  CAS  PubMed  Google Scholar 

  29. Choudhary D, et al. Isolated abdominal wall metastasis in early ovarian cancer: a true systemic metastasis or local tumour cells implantation. J Obstet Gynaecol. 2019;39(5):724–6.

    Article  PubMed  Google Scholar 

  30. McIntyre JB, et al. Molecular alterations in indolent, aggressive and recurrent ovarian low-grade serous carcinoma. Histopathology. 2017;70(3):347–58.

    Article  PubMed  Google Scholar 

  31. Zhao X, He M. Comprehensive pathway-related genes signature for prognosis and recurrence of ovarian cancer. PeerJ. 2020;8:e10437.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Chen M, et al. Pathway analysis of bladder cancer genome-wide association study identifies novel pathways involved in bladder cancer development. Genes Cancer. 2016;7(7–8):229–39.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Al Amri WS, et al. Identification of candidate mediators of chemoresponse in breast cancer through therapy-driven selection of somatic variants. Breast Cancer Res Treat. 2020;183(3):607–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Chen Y, et al. Identification of pathological grade and prognosis-associated lncRNA for ovarian cancer. J Cell Biochem. 2019;120(9):14444–54.

    Article  CAS  PubMed  Google Scholar 

  35. Ruan X, et al. Silencing LGR6 Attenuates Stemness and Chemoresistance via Inhibiting Wnt/β-Catenin Signaling in Ovarian Cancer. Mol Ther Oncolytics. 2019;14:94–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Shi Y, et al. Ceramide synthase 6 predicts poor prognosis and activates the AKT/mTOR/4EBP1 pathway in high-grade serous ovarian cancer. Am J Transl Res. 2020;12(9):5924–39.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Zhang J, et al. Association Between Succinate Receptor SUCNR1 Expression and Immune Infiltrates in Ovarian Cancer. Front Mol Biosci. 2020;7:150.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Fois SS, et al. Molecular epidemiology of the main druggable genetic alterations in non-small cell lung cancer. Int J Mol Sci. 2021;22(2):612.

  39. Hsu LW, et al. Genetic alterations in gastric cancer patients according to sex. Aging (Albany NY). 2020;13(1):376–88.

  40. Karunakaran K, Muniyan R. Genetic alterations and clinical dimensions of oral cancer: a review. Mol Biol Rep. 2020;47(11):9135–48.

    Article  PubMed  CAS  Google Scholar 

  41. Shen DH, et al. Epigenetic and genetic alterations of p33ING1b in ovarian cancer. Carcinogenesis. 2005;26(4):855–63.

    Article  CAS  PubMed  Google Scholar 

  42. Pongstaporn W, et al. Genetic alterations in chromosome 10q24.3 and glutathione S-transferase omega 2 gene polymorphism in ovarian cancer. J Exp Clin Cancer Res. 2006;25(1):107–14.

    CAS  PubMed  Google Scholar 

  43. Pećina-Šlaus N, et al. Mismatch Repair Pathway, Genome Stability and Cancer. Front Mol Biosci. 2020;7:122.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Wong RL, Yu EY. Refining Immuno-Oncology Approaches in Metastatic Prostate Cancer: Transcending Current Limitations. Curr Treat Options Oncol. 2021;22(2):13.

    Article  PubMed  Google Scholar 

  45. Boyero L, et al. Primary and acquired resistance to immunotherapy in lung cancer: unveiling the mechanisms underlying of immune checkpoint blockade therapy. Cancers (Basel). 2020;12(12):3729.

  46. Fan S, et al., Association between tumor mutation burden and immune infiltration in ovarian cancer. Int Immunopharmacol, 2020. 89(Pt A): p. 107126.

  47. Xiao X, et al. Mismatch repair deficiency is associated with MSI phenotype, increased tumor-infiltrating lymphocytes and PD-L1 expression in immune cells in ovarian cancer. Gynecol Oncol. 2018;149(1):146–54.

    Article  CAS  PubMed  Google Scholar 

  48. Wang J, et al. Modulation of immune infiltration of ovarian cancer tumor microenvironment by specific subpopulations of fibroblasts. Cancers (Basel). 2020;12(11):3184.

  49. Wall JA, et al. Manipulating the Wnt/β-catenin signaling pathway to promote anti-tumor immune infiltration into the TME to sensitize ovarian cancer to ICB therapy. Gynecol Oncol. 2021;160(1):285–94.

    Article  CAS  PubMed  Google Scholar 

  50. Kim MS, et al. Enhanced immunotherapy with LHRH-R targeted lytic peptide in ovarian cancer. Mol Cancer Ther. 2020;19(11):2396–406.

  51. Zhang QF, et al. CDK4/6 inhibition promotes immune infiltration in ovarian cancer and synergizes with PD-1 blockade in a B cell-dependent manner. Theranostics. 2020;10(23):10619–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Nayak AP, et al. Oxidative phosphorylation: a target for novel therapeutic strategies against ovarian cancer. Cancers (Basel). 2018;10(9):337.

  53. Xiao X, Melton DW, Gourley C. Mismatch repair deficiency in ovarian cancer -- molecular characteristics and clinical implications. Gynecol Oncol. 2014;132(2):506–12.

    Article  CAS  PubMed  Google Scholar 

  54. Saldivar JS, et al. Nucleotide excision repair pathway review I: implications in ovarian cancer and platinum sensitivity. Gynecol Oncol. 2007;107(1 Suppl 1):S56–71.

    Article  CAS  PubMed  Google Scholar 

  55. Zhang M, et al. The Association of Polymorphisms in Base Excision Repair Genes with Ovarian Cancer Susceptibility in Chinese Women: A Two-Center Case-Control Study. J Cancer. 2021;12(1):264–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank the researchers and study participants for their contributions.

Funding

None declared.

Author information

Authors and Affiliations

  1. Department of Obstetrics and Gynecology, Shengjing Hospital of China Medical University, No. 36 Sanhao street, Liaoning Province, 110004, Shenyang, P.R. China

    Xiaohan Chang

  2. Department of Anesthesiology, Shengjing Hospital of China Medical University, No. 36 Sanhao street, Liaoning Province, 110004, Shenyang, P.R. China

    Yunxia Dong

Authors
  1. Xiaohan Chang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yunxia Dong
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YX.D: Protocol/project development, data analysis and manuscript revision; XH.C: Manuscript writing/editing, data collection or management; The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Yunxia Dong.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the Institutional Research Ethics Committees of Shengjing Hospital of China Medical University and informed written consent was obtained from all subjects.

Consent for publication

Not applicable.

Competing interests

None declared.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Figure S1. Altered CACNA1C-related signaling pathways; (A) Cell-cycle-signaling-pathway; (B) RTK-RAS-PI(3)K-pathway; (C) Notch-signaling-pathway.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chang, X., Dong, Y. CACNA1C is a prognostic predictor for patients with ovarian cancer. J Ovarian Res 14, 88 (2021). https://doi.org/10.1186/s13048-021-00830-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13048-021-00830-z

Keywords

Journal of Ovarian Research

ISSN: 1757-2215

Contact us