GJA1 Expression and Its Prognostic Value in Cervical Cancer

Meng, Silu; Fan, Xinran; Zhang, Jianwei; An, Ran; Li, Shuang

doi:https://doi.org/10.1155/2020/8827920

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 8827920 | https://doi.org/10.1155/2020/8827920

GJA1 Expression and Its Prognostic Value in Cervical Cancer

Silu Meng,¹Xinran Fan,¹Jianwei Zhang,¹Ran An,²and Shuang Li¹

Academic Editor: Musthafa M. Essa

Received20 Sept 2020

Accepted30 Oct 2020

Published24 Nov 2020

Abstract

Gap Junction Protein Alpha 1 (GJA1) belongs to the gap junction family and has been widely studied in cancers. We evaluated the role of GJA1 in cervical cancer (CC) using public data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database. The difference of GJA1 expression level between CC and normal tissues was analyzed by the Gene Expression Profiling Interactive Analysis (GEPIA), six GEO datasets, and the Human Protein Atlas (HPA). The relationship between clinicopathological features and GJA1 expression was analyzed by the chi-squared test and the logistic regression. Kaplan–Meier survival analysis and Cox proportional hazard regression analysis were used to assessing the effect of GJA1 expression on survival. Gene set enrichment analysis (GSEA) was used to screen the signaling pathways regulated by GJA1. Immune Cell Abundance Identifier (ImmuCellAI) was chosen to analyze the immune cells affected by GJA1. The expression of GJA1 in CC was significantly lower than that in normal tissues based on the GEPIA, GEO datasets, and HPA. Both the chi-squared test and the logistic regression showed that high-GJA1 expression was significantly correlated with keratinization, hormone use, tumor size, and FIGO stage. The Kaplan–Meier curves suggested that high-GJA1 expression could indicate poor prognosis (). Multivariate analysis showed that high-GJA1 expression was an independent predictor of poor overall survival (HR, 4.084; 95% CI, 1.354-12.320; ). GSEA showed many cancer-related pathways, such as the p53 signaling pathway and the Wnt signaling pathway, were enriched in the high-GJA1-expression group. Immune cell abundance analysis revealed that the abundance of CD8 naive, DC, and neutrophil was significantly increased in the high-GJA1-expression group. In conclusion, GJA1 can be regarded as a potential prognostic marker of poor survival and therapeutic target in CC. Moreover, many cancer-related pathways may be the critical pathways regulated by GJA1. Furthermore, GJA1 can affect the abundance of immune cells.

1. Introduction

Cervical cancer (CC) is the fourth most common cancer in females, with 570,000 new cases and 311,000 deaths estimated for 2018 worldwide [1]. The vast majority of CC is in sub-Saharan Africa and South-Eastern Asia, and there were 98,900 new cases and 30,500 deaths for 2015 in China [1, 2]. Although the incidence and mortality rate of cervical declined due to vaccination, screening, and the control of precancerous lesions, it is even higher in some rural areas in China with the characteristics of young age (peak 45-49) and late stage [3–5]. In clinical, patients with early-stage CC (stages I to IIA) are mainly treated with surgery, whereas those with late-stage CC (stages IIB to IV) are treated with chemoradiotherapy [6]. However, the recurrence rate of CC is approximately 20%–25%, and the 5-year relative survival rate of CC remains poor in China (overall 45.4%), and the data from a retrospective population-based study in Hong Kong also showed that 5-year relative survival rates were 90.9%, 71.0%, 41.7%, and only 7.8% for the stages I, II, III, and IV, respectively [7–9]. Therefore, it is crucial to identify sensitive and specific biomarkers that could predict CC prognosis and serve as a target for CC treatment.

GJA1 maps to 6q22.31. GJA1 is also known as Connexin43 and belongs to the gap junction family [10]. Gap junction family plays an essential role in determining cellular phenotypes [10]. The role of gap junction family in cancers has been widely studied, among which GJA1 is one of the most studied genes. Firstly, GJA1 has been reported to play a suppressive role in tumorigenesis [11–13]. GJA1 could facilitate the transmission of cAMP, leading to increased p27 levels and reduced tumor growth by gap junctional intercellular communication [11]. Moreover, GJA1 could inhibit cell proliferation by inhibiting the Wnt signaling pathway [12, 13]. Secondly, although GJA1 has suppressive effects on cancer, many articles have reported the overexpression of GJA1 in metastatic tumors [14, 15]. Thirdly, GJA1 could also accelerate cancer progression by promoting cell migration via activating p38 [16, 17]. Also, high-GJA1 mRNA expression levels could indicate a more unsatisfactory outcome in cancer [18].

Based on the TCGA database and the GEO database, we firstly explored the difference of GJA1 expression level between CC and normal tissues. Also, we analyzed the relationship between GJA1 expression and the clinicopathological features of patients with CC and its prognostic value in CC. To gain further insight into the biological pathways and immune cell changes involved in CC pathogenesis related to GJA1, GSEA and immune cell abundance analysis were performed.

2. Materials and Methods

2.1. RNA-Sequencing Patient Data in the TCGA Database

The gene expression data (304 cases, workflow type: HTSeq-FPKM-UQ) and corresponding clinical information were downloaded from TCGA database (https://portal.gdc.cancer.gov) by TCGAbiolinks package [19]. The details of clinical information are shown in Table 1. Besides, it should be emphasized that survival information, including vital status and overall survival (OS) time, was directly taken from the articles of Liu et al. published in Cell in 2018 [20].

2.2. Relative GJA1 Expression Level between CC and Normal Tissues

GEPIA (http://gepia.cancer-pku.cn/) was used to access the mRNA expression levels of GJA1 in CC (based on TCGA database) and normal tissues (based on GTEx project) [21]. limma package [22] was used to explore whether the expression of GJA1 is significantly different between CC and normal tissues in six GEO datasets: GSE39001 [23], GSE52903 [24], GSE63514 [25], GSE6791 [26], GSE7803 [27], and GSE9705 [28]. Meeting the criteria of | and adjusted value < 0.05, GJA1 could be considered differentially expressed between CC and normal tissues. Additionally, we explored the expression of GJA1 at the protein level by the HPA database [29].

2.3. Survival Analysis

We matched the GJA1 expression data with survival information derived from Liu et al. [20]. The patients in cohorts were divided into two groups according to the median expression levels of GJA1 (the low-GJA1-expression group and the high-GJA1-expression group). The survival and survminer packages were used for survival analysis and visualization, and then the Kaplan-Meier survival curve was obtained. Also, we conducted a stratified analysis and explored the prognostic value of GJA1 in squamous cell carcinoma (SCC) patients, early-stage patients (FIGO ), and late-stage patients (FIGO ).

2.4. Univariate and Multivariate Cox Regression Analyses

The Cox proportional hazard regression model was used to conduct univariate and multivariate analyses. The hazard ratio (HR) value and 95% confidence interval (CI) were calculated. In univariate analysis, the independent predictive value of the clinicopathological parameters and GJA1 expression on survival was evaluated. Multivariate Cox analysis was conducted to compare the influence of GJA1 expression on survival along with other characteristics as categorical variables. Variables with value ≤ 0.05 or close to 0.05, including age, BMI, tumor size, lymph node, FIGO stage, lymphovascular invasion, and distant metastasis, were entered in the multivariate Cox regression analysis as categorical variables. The cut-off value of GJA1 expression was set based on the median expression value. The data were analyzed by using survival and survminer packages and visualized by forestplot package in R.

2.5. Gene Set Enrichment Analysis (GSEA)

GSEA (version 4.0.3) was used to explore the signaling pathways related to GJA1 in cervical cancer [30]. Gene expression enrichment analysis was carried out between different phenotypes determined by the GJA1 expression level. The annotated gene set was selected (c2.cp.kegg.v7.1.symbols.gmt) as the reference gene set. Gene set permutations were performed 1,000 times for analysis. The normalized enrichment score (NES), nominal value, and false discovery rate (FDR) -value were used to sort the pathways enriched in each group. Pathways with , nominal value < 0.05, and FDR were selected out.

2.6. Immune Cell Abundance Analysis

ImmuCellAI is a tool to estimate the abundance of 24 immune cells from gene expression datasets [31]. We used this tool to survey the infiltration difference of immune cells between the low- and the high-GJA1 expression groups.

2.7. Statistical Analysis

All statistical analyses were performed with IBM SPSS statistical software (version 26.0) and R software (version 3.6.3), and was used to be as the significance level. The differences in GJA1 expression between the two groups were compared by Wilcoxon test. Chi-squared () test and logistic regression analysis were used to evaluate the relationship between GJA1 expression and clinicopathological parameters. Kaplan–Meier analysis and log-rank test were used to explore the relationship between survival rates and the GJA1 expression level. A Cox proportional hazard regression model was used for univariate and multivariate survival analysis.

3. Results

3.1. The Difference of GJA1 Expression between CC and normal tissues

GEPIA was used to explore the mRNA expression levels of GJA1 in CC and normal tissues, showing that GJA1 expression is significantly decreased in CC tissues (Figure S1A). Similar results were observed in five of the six GEO datasets.(Table S1). On the protein level, the expression of GJA1 was also higher in normal tissues than in CC tissues by using the HPA database (Figure S1B). In addition, the expression level of GJA1 was different in groups classified according to histology (, Figure 1(a)) and FIGO stage (, Figure 1(b)). In a stratified analysis of SCC patients, we found the same results that the GJA1 expression level was related to the FIGO stage (, Figure 1(d)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 1

The expression of GJA1 and its association with clinicopathological parameters based on TCGA data: (a) histology type; (b) FIGO stage of all CC patients; (c) impact of GJA1 expression on OS in all CC patients; (d) FIGO stage of SCC patients; (e) impact of GJA1 expression on OS in all SCC patients; (f) impact of GJA1 expression on OS in patients with FIGO ; (g) impact of GJA1 expression on OS in patients with FIGO . TCGA: The Cancer Genome Atlas; FIGO: the International Federation of Gynecology and Obstetrics; CC: cervical cancer; OS: overall survival; SCC: squamous cell carcinoma; ACC: adenomas and adenocarcinomas; GJA1: gap junction protein alpha 1; , , and .

3.2. High Expression of GJA1 in CC Is Related to Poor Overall Survival

We evaluated the relationship between GJA1 expression and prognosis of CC patients by Kaplan–Meier risk estimates. Compared to the low-GJA1 expression, the high-GJA1 expression was significantly associated with poor overall survival (, Figure 1(c)). The median OS of the high-GJA1-expression group was 627.5 days (range: 0-5271 days), while for the low-GJA1-expression group, the median OS was 734.5 days (range: 0-6408 days). The 5-year survival rate of patients in the low-GJA1-expression group (32.9%) was also higher than that in the high-GJA1-expression group (25.0%). In the stratified analysis of SCC patients, early-stage patients (), and late-stage patients (), the high-GJA1 expression group all showed poor prognosis, and the value was 0.004, 0.042, and 0.031, respectively (Figures 1(e)–1(g)).

3.3. The relationship between GJA1 Expression and Clinicopathological Variables

To further explore the relationship between GJA1 expression and clinicopathological parameters, the clinical data of 304 CC samples with GJA1 expression data were analyzed from the TCGA database. Using chi-squared () test, the high expression level of GJA1 was significantly correlated with keratinization (), hormone use (), tumor size (0.042), and FIGO stage () (Table 2). Multiple logistic regression analysis showed that the increased expression of GJA1 in CC was significantly correlated with keratinization (OR, 1.987; 95% CI, 1.011-3.905; ), hormone use (OR, 0.457; 95% CI, 0.236-0.885; ), tumor size (OR, 2.234; 95% CI, 1.013-4.927; ), and FIGO stage (OR, 2.062; 95% CI, 1.274-3.338; ) (Table 3).

3.4. The Effect of GJA1 Expression on Survival Based on Univariate and Multivariate Analyses

The univariate and multivariate Cox proportional hazard regression analyses were performed to investigate whether the high-GJA1 expression is an independent predictor of poor survival in CC patients. After excluding patients with incomplete data, we included 120 patients with CC for the multivariate Cox proportional hazard regression analysis. The univariate analysis showed that BMI (HR, 0.953; 95% CI, 0.911-0.996; ), tumor size (HR, 3.643, 95% CI, 1.922-6.904, ), lymph node (HR, 2.695; 95% CI, 1.358-5.349; ), FIGO stage (HR, 1.861; 95% CI, 1.166-2.970; ), lymphovascular invasion (HR, 10.041; 95% CI, 2.361-42.700; ), distant metastasis (HR, 3.141; 95% CI, 1.866-5.289; ), and high-GJA1 expression (HR, 1.943; 95% CI, 1.202-3.140; ) were important predictors of survival (Table 4). As the value is close to 0.05, age was also included in multivariate analysis. The results showed that tumor size (HR, 6.181; 95% CI, 1.219-31.334; ), lymphovascular invasion (HR, 5.910; 95% CI, 1.124-31.059; ), and the high-GJA1 expression (HR, 4.084; 95% CI, 1.354-12.320; ) were the important independent predictors of poor overall survival of CC (Figure 2 and Table 4).

3.5. Identification of GJA1-Related Signaling Pathways by GSEA

To identify signaling pathways that are differentially activated in CC, we conducted the Gene Set Enrichment Analysis (GSEA) between low- and high-GJA1-expression datasets. We selected out 30 significantly enriched signaling pathways based on the standard that is , nominal value < 0.05, and FDR (Table S2). Figure 3 shows 14 cancer-related pathways enriched in high-GJA1-expression groups, including many well-known pathways, such as the P53 signaling pathway, the Wnt signaling pathway, and the MAPK signaling pathway.

3.6. Immune Cell Abundance Analysis

ImmuCellAI was used to analyze the infiltration difference of immune cells between the low- and the high-GJA1-expression group. We found that the abundance of CD8 naive, DC, and neutrophil is significantly increased in the high-GJA1-expression group, while the abundance of gamma delta and Th1 is significantly decreased in the low-GJA1-expression group (Figure 4).

4. Discussion

The role of GJA1 in cancers has been widely studied in recent years. GJA1 has two opposite effects on cancers as it can act as an oncogene or a tumor suppressor gene. Many studies have unveiled the suppressive roles of GJA1 in tumorigenesis. In the colorectal cancer cell line (HT 29), Sirnes et al. found overexpression of GJA1 could inhibit Wnt signaling by interacting with β-catenin, thus inhibiting the growth of tumor cells [32]. Huang et al. identified that GJA1 could restrain glioblastoma development by reducing the anti-apoptotic Bcl-2 [13]. Similarly, overexpression of GJA1 was also found to inhibit the growth of lung cancer cells (LH7) [33]. Many studies have also reported the role of GJA1 in promoting cancer progression. Tang et al. reported that the overexpression of GJA1 could promote lymph node metastasis and peritoneal metastasis in gastric cancer [34, 35]. Zhao et al. found that high-GJA1 expression could indicate poor prognosis of gastric cancer patients based on five GEO datasets [18]. Ogawa et al. found that the silencing of GJA1 was associated with reduced invasion, migration, and metastasis [36]. From the Affymetrix analysis, Teleki et al. reported that high expression of GJA1 was associated with a reduced relapse-free survival (RFS) and overall survival (OS) in ER-negative breast cancer patients [37]. Moreover, the study by Stoletov et al. confirmed this observation using analysis from the Oncomine database [38]. Also, there are some researches on the impact of GJA1 on CC. Sun et al. found that HPV E6 may regulate GJA1 trafficking in cervical tumor cells, resulting in inhibition of the formation of gap junctions [39]. Using whole-genome microarray data, Cheng et al. showed that GJA1 was a critical gene for CC invasion and metastasis [40].

Using high-throughput RNA-sequencing data from the TCGA database, we aimed to explore the relationship between clinicopathological parameters and GJA1 expression and determine the role of GJA1 expression in CC progression, especially as a prognostic factor for CC. We also tried to screen signaling pathways related to GJA1 in CC to understand the underlying mechanism involved in the regulation of CC development by GJA1. In addition, we tried to explore GJA1 expression’s impact on the abundance of the immune cells. First, we compared the expression of GJA1 in CC and normal tissues by GEPIA. The expression level of GJA1 in CC tissues was significantly lower than that in normal tissues. Six GEO datasets were used to verify this result, and the results from five datasets are the same as those from GEPIA. The HPA database also showed that the expression of GJA1 is higher in normal tissues than that in CC tissues at the protein level. These results suggest that GJA1 may be a tumor suppressor gene and play an important role in the progression of CC. In addition, GJA1 expression levels were different in groups classified by histology and FIGO stage. Further, chi-squared () test showed that the high expression level of GJA1 was significantly correlated with keratinization, hormone use, tumor size, and FIGO stage. Some literatures reported that GJA1 was associated with metastasis in a variety of tumors. Lin et al. found that the expression of GJA1 was higher in metastatic breast tumors than in primary breast tumors [41]. Wang et al. also showed that GJA1 could be used as a marker of metastasis in prostate cancer [42].

The Kaplan–Meier survival analysis showed that the high-GJA1-expression group indicated a worse prognosis than the low-GJA1-expression group. In the stratified analysis of SCC patients, early-stage patients (), and late-stage patients (), we found the same result. Using five GEO datasets, Zhao et al. reported that high expression level of GJA1 was associated with poor prognosis in gastric cancer, which is similar to our result [18].

The univariate analysis showed that BMI, tumor size, lymph node, FIGO stage, lymphovascular invasion, distant metastasis, and high-GJA1expression were important predictors of survival. Importantly, multivariate analysis showed that tumor size, lymphovascular invasion, and the high-GJA1 expression were the important independent predictors of poor overall survival of CC. This result demonstrated the potential of GJA1 to become a prognosis biomarker of CC.

To further investigate the functions of GJA1 in CC, we performed GSEA using the TCGA data and 30 pathways were selected out based on established standards. The result showed that many cancer-related pathways, such as the p53 signaling pathway, the Wnt signaling pathways, and the MAPK signaling pathway, were enriched in the high-GJA1-expression group. Evidence has shown that the proliferation, migration, and invasion of CC cells can be modulated by regulating the p53 signaling pathway [43, 44]. Wnt signaling, which is needed for cell proliferation and differentiation, has been found to play a critical role in CC carcinogenesis [45, 46]. Activation of this pathway was shown to be related to the poor prognosis of CC patients [46]. MAPK signaling pathways regulate a variety of cellular activities, including proliferation, differentiation, survival, and death [47]. Many studies have reported that the interaction between HPV and the MAPK signaling pathway indirectly explains its role in CC [48, 49]. In addition, the TGF-β signaling pathway was enriched in the high-GJA1-expression group, which participates in the regulation of a variety of immune cells, such as NK cell, neutrophils, and macrophage [50]. Thus, we analyzed the differences of immune cells between the low-GJA1-expression group and the high-GJA1-expression group and found that the abundance of CD8 naive cell, DC, and neutrophil were significantly increased in the high-GJA1-expression group.

5. Conclusions

In conclusion, high GJA1 expression may be a potential prognostic molecular marker of poor survival in CC. Moreover, the p53 signaling pathway, Wnt signaling pathway, and MAPK signaling pathways may be the critical pathways regulated by GJA1. Moreover, GJA1 can affect the abundance of immune cells to a certain extent. Further experimental validation should be performed to prove the biologic impact of GJA1 in CC.

Data Availability

Publicly available datasets were analyzed in this study. These can be found in The Cancer Genome Atlas (https://portal.gdc.cancer.gov/); the NCBI Gene Expression Omnibus (GSE39001, GSE52903, GSE63514, GSE6791, GSE7803, GSE9750).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Silu Meng, Xinran Fan, and Jianwei Zhang contributed equally to this work. Silu Meng, Xinran Fan, and Jianwei Zhang designed the overall idea of this study, analyzed the data, prepared the figures and tables, and wrote the drafts of the paper. Ran An checked the statistical results. Shuang Li supervised this study and revised the drafts of the paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant number: 81974410).

Supplementary Materials

Supplementary 1. Supplementary Figure 1: comparison of GJA1 expression between cervical cancer and normal using GEPIA and The Human Protein Atlas.

Supplementary 2. Supplementary Table 1: information of selected six GEO datasets and GJA1 expression in these datasets. Supplementary Table 2: gene sets enriched in the high-GJA1-expression phenotype.

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
View at: Publisher Site | Google Scholar
S. Yu, C. S. Yang, J. Li et al., “Cancer prevention research in China,” Cancer Prevention Research (Philadelphia, Pa.), vol. 8, no. 8, pp. 662–674, 2015.
View at: Publisher Site | Google Scholar
X. Jiang, H. Tang, and T. Chen, “Epidemiology of gynecologic cancers in China,” Journal of Gynecologic Oncology, vol. 29, no. 1, 2018.
View at: Publisher Site | Google Scholar
W. Chen, R. Zheng, S. Zhang et al., “Cancer incidence and mortality in China, 2013,” Cancer Letters, vol. 401, pp. 63–71, 2017.
View at: Publisher Site | Google Scholar
C. Marth, F. Landoni, S. Mahner et al., “Cervical cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up,” Ann Oncol, vol. 28, supplement_4, pp. iv72–iv83, 2017.
View at: Publisher Site | Google Scholar
A. P. M. Heintz, F. Odicino, P. Maisonneuve et al., “Carcinoma of the fallopian tube. FIGO 26th Annual Report on the Results of Treatment in Gynecological Cancer,” International Journal of Gynecology & Obstetrics, vol. 95, pp. S145–S160, 2006.
View at: Publisher Site | Google Scholar
H. Zeng, R. Zheng, Y. Guo et al., “Cancer survival in China, 2003-2005: a population-based study,” International Journal of Cancer, vol. 136, no. 8, pp. 1921–1930, 2015.
View at: Publisher Site | Google Scholar
F. Y. Cheung, O. W. Mang, and S. C. Law, “A population-based analysis of incidence, mortality, and stage-specific survival of cervical cancer patients in Hong Kong: 1997-2006,” Hong Kong Medical Journal, vol. 17, no. 2, pp. 89–95, 2011.
View at: Google Scholar
J. I. Wu and L. H. Wang, “Emerging roles of gap junction proteins connexins in cancer metastasis, chemoresistance and clinical application,” Journal of Biomedical Science, vol. 26, no. 1, p. 8, 2019.
View at: Publisher Site | Google Scholar
Y. W. Zhang, I. Morita, M. Ikeda, K. W. Ma, and S. Murota, “Connexin43 suppresses proliferation of osteosarcoma U2OS cells through post-transcriptional regulation of p27,” Oncogene, vol. 20, no. 31, pp. 4138–4149, 2001.
View at: Publisher Site | Google Scholar
Z. Ai, A. Fischer, D. C. Spray, A. M. C. Brown, and G. I. Fishman, “Wnt-1 regulation of connexin43 in cardiac myocytes,” The Journal of Clinical Investigation, vol. 105, no. 2, pp. 161–171, 2000.
View at: Publisher Site | Google Scholar
R. P. Huang, M. Z. Hossain, R. Huang, J. Gano, Y. Fan, and A. L. Boynton, “Connexin 43 (cx43) enhances chemotherapy-induced apoptosis in human glioblastoma cells,” International Journal of Cancer, vol. 92, no. 1, pp. 130–138, 2001.
View at: Publisher Site | Google Scholar
L. Kanczuga-Koda, S. Sulkowski, A. Lenczewski et al., “Increased expression of connexins 26 and 43 in lymph node metastases of breast cancer,” Journal of Clinical Pathology, vol. 59, no. 4, pp. 429–433, 2006.
View at: Publisher Site | Google Scholar
J. Ming, Y. Zhou, J. du et al., “miR-381 suppresses C/EBPα-dependent Cx43 expression in breast cancer cells,” Bioscience Reports, vol. 35, no. 6, 2015.
View at: Publisher Site | Google Scholar
J. Behrens, P. Kameritsch, S. Wallner, U. Pohl, and K. Pogoda, “The carboxyl tail of Cx43 augments p38 mediated cell migration in a gap junction-independent manner,” European Journal of Cell Biology, vol. 89, no. 11, pp. 828–838, 2010.
View at: Publisher Site | Google Scholar
S. Ghosh, A. Kumar, R. P. Tripathi, and S. Chandna, “Connexin-43 regulates p38-mediated cell migration and invasion induced selectively in tumour cells by low doses of γ-radiation in an ERK-1/2-independent manner,” Carcinogenesis, vol. 35, no. 2, pp. 383–395, 2014.
View at: Publisher Site | Google Scholar
X. Zhao, C. Yu, M. Zheng, and J. Sun, “Prognostic value of the mRNA expression of gap junction α members in patients with gastric cancer,” Oncology Letters, vol. 18, pp. 1669–1678, 2019.
View at: Publisher Site | Google Scholar
A. Colaprico, T. C. Silva, C. Olsen et al., “TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data,” Nucleic Acids Research, vol. 44, no. 8, 2016.
View at: Publisher Site | Google Scholar
J. Liu, T. Lichtenberg, K. A. Hoadley et al., “An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics,” Cell, vol. 173, no. 2, pp. 400–416.e11, 2018.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, p. e47, 2015.
View at: Publisher Site | Google Scholar
A. M. Espinosa, A. Alfaro, E. Roman-Basaure et al., “Mitosis is a source of potential markers for screening and survival and therapeutic targets in cervical cancer,” PLoS One, vol. 8, no. 2, 2013.
View at: Publisher Site | Google Scholar
I. Medina-Martinez, V. Barrón, E. Roman-Bassaure et al., “Impact of gene dosage on gene expression, biological processes and survival in cervical cancer: a genome-wide follow-up study,” PLoS One, vol. 9, no. 5, 2014.
View at: Publisher Site | Google Scholar
J. A. den Boon, D. Pyeon, S. S. Wang et al., “Molecular transitions from papillomavirus infection to cervical precancer and cancer: role of stromal estrogen receptor signaling,” Proceedings of the National Academy of Sciences of the United States of America, vol. 112, no. 25, pp. E3255–E3264, 2015.
View at: Publisher Site | Google Scholar
D. Pyeon, M. A. Newton, P. F. Lambert et al., “Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers,” Cancer Research, vol. 67, no. 10, pp. 4605–4619, 2007.
View at: Publisher Site | Google Scholar
Y. Zhai, R. Kuick, B. Nan et al., “Gene expression analysis of preinvasive and invasive cervical squamous cell carcinomas identifies HOXC10 as a key mediator of invasion,” Cancer Research, vol. 67, no. 21, pp. 10163–10172, 2007.
View at: Publisher Site | Google Scholar
L. Scotto, G. Narayan, S. V. Nandula et al., “Identification of copy number gain and overexpressed genes on chromosome arm 20q by an integrative genomic approach in cervical cancer: potential role in progression,” Genes, Chromosomes & Cancer, vol. 47, no. 9, pp. 755–765, 2008.
View at: Publisher Site | Google Scholar
F. Ponten, K. Jirstrom, and M. Uhlen, “The Human Protein Atlas--a tool for pathology,” The Journal of Pathology, vol. 216, no. 4, pp. 387–393, 2008.
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
Y. R. Miao, Q. Zhang, Q. Lei et al., “ImmuCellAI: a unique method for comprehensive T-cell subsets abundance prediction and its application in cancer immunotherapy,” Advanced Science, vol. 7, no. 7, p. 1902880, 2020.
View at: Publisher Site | Google Scholar
S. Sirnes, J. Bruun, M. Kolberg et al., “Connexin43 acts as a colorectal cancer tumor suppressor and predicts disease outcome,” International Journal of Cancer, vol. 131, no. 3, pp. 570–581, 2012.
View at: Publisher Site | Google Scholar
H. T. Xu, Q. C. Li, Y. X. Zhang et al., “Connexin 43 recruits E-cadherin expression and inhibits the malignant behaviour of lung cancer cells,” Folia Histochemica et Cytobiologica, vol. 46, no. 3, pp. 315–321, 2008.
View at: Publisher Site | Google Scholar
B. Tang, Z. H. Peng, P. W. Yu, G. Yu, and F. Qian, “Expression and significance of Cx43 and E-cadherin in gastric cancer and metastatic lymph nodes,” Medical Oncology, vol. 28, no. 2, pp. 502–508, 2011.
View at: Publisher Site | Google Scholar
B. Tang, Z. H. Peng, P. W. Yu et al., “Aberrant expression of Cx43 is associated with the peritoneal metastasis of gastric cancer and Cx43-mediated gap junction enhances gastric cancer cell diapedesis from peritoneal mesothelium,” PLoS One, vol. 8, no. 9, 2013.
View at: Publisher Site | Google Scholar
K. Ogawa, P. Pitchakarn, S. Suzuki et al., “Silencing of connexin 43 suppresses invasion, migration and lung metastasis of rat hepatocellular carcinoma cells,” Cancer Science, vol. 103, no. 5, pp. 860–867, 2012.
View at: Publisher Site | Google Scholar
I. Teleki, A. M. Szasz, M. E. Maros et al., “Correlations of differentially expressed gap junction connexins Cx26, Cx30, Cx32, Cx43 and Cx46 with breast cancer progression and prognosis,” PLoS One, vol. 9, no. 11, 2014.
View at: Publisher Site | Google Scholar
K. Stoletov, J. Strnadel, E. Zardouzian et al., “Role of connexins in metastatic breast cancer and melanoma brain colonization,” Journal of Cell Science, vol. 126, no. 4, pp. 904–913, 2013.
View at: Publisher Site | Google Scholar
P. Sun, L. Dong, A. MacDonald et al., “HPV16 E6 controls the gap junction protein Cx43 in cervical tumour cells,” Viruses, vol. 7, no. 10, pp. 5243–5256, 2015.
View at: Publisher Site | Google Scholar
Y. Cheng, D. Ma, Y. Zhang, Z. Li, and L. Geng, “Cervical squamous cancer mRNA profiles reveal the key genes of metastasis and invasion,” European Journal of Gynaecological Oncology, vol. 36, no. 3, pp. 309–317, 2015.
View at: Google Scholar
Z. J. Lin, J. Ming, L. Yang, J. Z. du, N. Wang, and H. J. Luo, “Mechanism of regulatory effect of microRNA-206 on connexin 43 in distant metastasis of breast cancer,” Chinese Medical Journal, vol. 129, no. 4, pp. 424–434, 2016.
View at: Publisher Site | Google Scholar
Y. WANG, W. GUO, H. XU et al., “An extensive study of the mechanism of prostate cancer metastasis,” Neoplasma, vol. 65, no. 2, pp. 253–261, 2018.
View at: Publisher Site | Google Scholar
Y. Liu, L. Li, Y. Liu et al., “RECK inhibits cervical cancer cell migration and invasion by promoting p53 signaling pathway,” Journal of Cellular Biochemistry, vol. 119, no. 4, pp. 3058–3066, 2018.
View at: Publisher Site | Google Scholar
V. K. Kashyap, N. Dan, N. Chauhan et al., “VERU-111 suppresses tumor growth and metastatic phenotypes of cervical cancer cells through the activation of p53 signaling pathway,” Cancer Letters, vol. 470, pp. 64–74, 2020.
View at: Publisher Site | Google Scholar
A. Bahrami, F. Amerizadeh, S. ShahidSales et al., “Therapeutic potential of targeting Wnt/β-catenin pathway in treatment of colorectal cancer: rational and progress,” Journal of Cellular Biochemistry, vol. 118, no. 8, pp. 1979–1983, 2017.
View at: Publisher Site | Google Scholar
A. Bahrami, M. Hasanzadeh, S. ShahidSales et al., “Clinical significance and prognosis value of Wnt signaling pathway in cervical cancer,” Journal of Cellular Biochemistry, vol. 118, no. 10, pp. 3028–3033, 2017.
View at: Publisher Site | Google Scholar
E. K. Kim and E. J. Choi, “Pathological roles of MAPK signaling pathways in human diseases,” Biochimica et Biophysica Acta, vol. 1802, no. 4, pp. 396–405, 2010.
View at: Publisher Site | Google Scholar
O. Chakrabarti, K. Veeraraghavalu, V. Tergaonkar et al., “Human papillomavirus type 16 E6 amino acid 83 variants enhance E6-mediated MAPK signaling and differentially regulate tumorigenesis by notch signaling and oncogenic Ras,” Journal of Virology, vol. 78, no. 11, pp. 5934–5945, 2004.
View at: Publisher Site | Google Scholar
A. Contreras-Paredes, E. De la Cruz-Hernandez, I. Martinez-Ramirez, A. Duenas-Gonzalez, and M. Lizano, “E6 variants of human papillomavirus 18 differentially modulate the protein kinase B/phosphatidylinositol 3-kinase (akt/PI3K) signaling pathway,” Virology, vol. 383, no. 1, pp. 78–85, 2009.
View at: Publisher Site | Google Scholar
E. Batlle and J. Massague, “Transforming growth factor-β signaling in immunity and cancer,” Immunity, vol. 50, no. 4, pp. 924–940, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Silu Meng 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

Views

639

Downloads

1107

Citations