Prognostic value and biological signicance of GUCY1A2 in gastric cancer: a bioinformatics analysis

Background: Guanylate cyclase 1 soluble subunit alpha 2 (sGCα2), also known as GUCY1A2, was reported to be upregulated and promoted tumorigenesis in some cancers. But whether GUCY1A2 was abnormally expressed and the prognostic value in gastric cancer was unknown. The current study aimed to nd out the prognostic value of GUCY1A2 in gastric cancer by analyzing data from The Cancer Genome Atlas (TCGA). Methods: Wilcoxon signed-rank test, cox regression analysis and multivariant analysis were used to analyze the relationship between clinical characteristic and GUCY1A2 expression level. Kaplan-Meier method was used to analyze the association of GUCY1A2 and overall survival. Gene set enrichment analysis (GSEA) was used to identify GUCY1A2-related signaling pathway. Results: Compared to normal tissue, expression of GUCY1A2 was signicantly increased in gastric cancer (p=0.000). Increased GUCY1A2 was associated with advance T stage (p=0.012) and poor survival (p=0.022). Univariate analysis showed that high GUCY1A2 expression was associated with a poor overall survival (HR:1.44, 95% condence interval [CI]: 1.03-2.02, p=0.03). Multivariate analysis indicated that GUCY1A3 remained an independent prognostic predictor of overall survival (HR:1.75, 95% condence interval [CI]: 1.20-2.56, p=0.00). GSEA revealed that calcium signaling pathway, MAPK signaling pathway, TGF-β signaling pathway and Wnt signaling pathway were enriched in GUCY1A2 high expression phenotype. Conclusions: GUCY1A2 maybe a potential prognostic predictor of poor survival in gastric cancer. Calcium signaling pathway, MAPK signaling pathway, TGF-β signaling pathway and Wnt signaling pathway may be the key signaling pathway regulated by GUCY1A2. in expression phenotype. may serve as a predictive biomarker of poor prognosis in gastric


Introduction
Gastric cancer is one of the most common malignant tumors in the world, with the fth highest incidence rate and the third highest mortality rate (1). The incidence rate of gastric cancer shows geographical differences. It is most frequently occurred in East Asia, South America and Eastern Europe and less common seen in Western Europe, Australia and North America (2). More than 95% of gastric cancers are classi ed as adenocarcinomas (3). Helicobacter pylori infection is the most important cause of gastric cancer (4). Other risk factors include smoking, alcohol consumption, high salt intake, et al (5)(6)(7). The only curative therapy for gastric cancer is surgical resection (8). However, about 50% of patients presents with advanced diseases at diagnosis and about 40-60% of patients relapsed after surgical resection (9). The standard treatment of these patients is systemic chemotherapy and target therapy (10). However, prognosis of metastasis gastric cancer is poor, with the 5-year overall survival (OS) rate of 5-20% (11).
Searching for new prognostic biomarkers and treatment targets may help to improve the prognosis of gastric cancer.
The cyclic GMP (cGMP)/protein Kinase G (PKG) pathway participates in the endogenous apoptotic process in several types of cancers (12)(13)(14). cGMP, which is one of the most important intracellular signaling molecule in mammal cells, is generated from GTP by two classes of guanylate cyclases (GCs): soluble GC (sGC) and particulate GC (pGC) (15). Guanylate cyclase 1 soluble subunit alpha 1 (sGCα1) and alpha 2 (sGCα2) are subunits of sGCs. It is reported that sGCα1 is increased in prostate cancer and overexpression of sGCα1 is related to chemoresistance (16). Another research showed that sGCα1 promoted cell proliferation, migration and survival of cervical cancer cells (17). Inhibition of basal activity of sGC/cGMP pathway induced p53-dependent apoptosis in ovarian cancer cells. Guanylate cyclase soluble subunit alpha-2, also known as GUCY1A2, is found to be upregulated in cervical cancer and correlated with super cial tumor growth (18). However, the expression level and prognostic value of GUCY1A2 in gastric cancer remains unknown.
Thus, the current study aimed at evaluating the expression and prognostic value of GUCY1A2 in gastric cancer by bioinformatics analysis. In addition, we performed the gene set enrichment analysis (GSEA) to study the biological pathway related to the GUCY1 regulation. The data used in this study was obtained from TCGA database (https://portal.gdc.cancer.gov/).

TCGA data download and bioinformatic analysis
The gene expression data and corresponding clinical data of gastric cancer (project ID: TCGA-STAD) were downloaded from TCGA database (https://portal.gdc.cancer.gov/). Expression data of GUCY1A2 was extracted. Then the expression difference of GUCY1A2 between normal and tumor samples were compared. Clinical information was extracted and patients with incomplete clinical information were excluded. Finally, there were 349 patients were included into cox regression and multivariate cox analysis. The media expression value of GUCY1A2 were calculated. Patients with GUCY1A2 expression level higher than media value were distributed to the high expression group and others were distributed to the low expression group. Overall survival in high expression and low expression group were compared. Cox regression and multivariate cox analysis were used to analysis the in uence of clinical characteristics on the survival.

Gene set enrichment analysis
Software GSEA_4.1.0 was used to perform gene set enrichment analysis (GSEA). The expression matrix le and phenotype le were prepared before GSEA. GUCY1A2 expression level was used as a phenotype label. Gene set permutations number was set to 1000. Then the expression matrix le and phenotype le were imported and GSEA was performed.

Statistical analysis
Statistical analyses were conducted by R software (R x64 3.6.2). The comparison of GUCY1A2 expression between normal and tumor sample was analyzed by Wilcoxon signed-rank test. Clinical characteristics associated to survival was analyzed by Cox regression and Kaplan-Meier method. To further compare the effect of GUCY1A2 expression on overall survival along with other clinical factors,

GUCY1A2 expression was increased in gastric carcinoma patients
Data of 343 gastric carcinoma samples and 30 normal samples with GUCY1A2 expression information were downloaded from TCGA database in January 2021. GUCY1A2 expression level were compared in normal samples and tumors samples. In comparison to normal samples, GUCY1A2 expression was signi cantly increased in tumor samples ( Fig. 1A and 1B).

Increased GUCY1A2 was associated with poor overall survival
Patients with complete survival information and GUCY1A2 expression information were included in survival analysis. As shown in Fig. 2A, patients with high GUCY1A2 expression showed a worse prognosis than patients with low GUCY1A2 expression. We also analyzed the association between GUCY1A2 expression and other clinical characteristics. Results showed that expression level of GUCY1A2 was correlated to histological grade, T stage. Patients who with higher histological grade and T stage had higher GUCY1A2 level. But GUCY1A2 expression was independent of age, gender, clinical stage, lymph node invasion or distant metastasis ( Fig. 2B-2H).

Discussion
Soluble guanylate cyclase (sGC) is a catalytic enzyme that catalyzing the conversion of guanosine 5′triphosphate (GTP) into cyclic guanosine 3′,5′-monophosphate (cGMP). sGC is one of the main receptors of nitric oxide (NO). When sGC binds to NO, the catalytic activity was increased by 100-to 200-folds. Thus it plays an important role in the NO-cGMP pathway (19). It has been shown that NO/sGC/cGMP pathway plays an important role in tumor development. Increased cGMP leads to increased protein kinase G (PKG) activity and promotes prostate cancer cell proliferation (20). Inhibition of NO/cGMP/PKG pathway results in suppression of migration and invasion in breast cancer cells (21).
sGC a 150kDa heterodimer consisting of two subunits, subunit α and subunit β (22). Both of the subunits are required for catalytic activity (23). There are two isoforms of both subunits α1, α2 and β1 β2. α1/β1 is the most abundant form. But there is indistinguishable activity between the α1/β1 and α2/β1 heterodimers (24). It has been shown that sGC subunit α1 and α2, also known as GUCY1A1 and GUCY1A2, is increased in some type of cancers and promotes cell proliferation, migration and survival of cancer cells (16)(17)(18). But whether GUCY1A2 was abnormally expressed and the prognostic value in gastric cancer was unknown.
In the current study, we found that GUCY1A2 expression was increased in gastric cancer in comparison to normal tissues. Increased GUCY1A2 expression was associated with advance T stage and poor prognosis in gastric cancer. GSEA indicated that adherens junction, calcium signaling pathway, cell adhesion molecules cams, ECM receptor interaction, focal adhesion, MAPK signaling pathway, TGF-β signaling pathway, Wnt signaling pathway was enriched in GUCY1A2 high expression phenotype. Results suggested that high level of GUCY1A2 may serve as a predictive biomarker of poor prognosis in gastric cancer.
Some research has indicated that upregulation of the calcium signaling pathway, Wnt signaling pathway, MAPK signaling pathway and TGF-β signaling pathway were associated with tumor development and progression (25)(26)(27)(28). A recent study showed that increased cGMP resulted in increased tumor cell stemness and promoted metastasis in breast cancer by activating the cGMP-dependent PKG and MAPK signaling pathway (29). Activation of sGC and increased generation of cGMP also induced activation of TGF-β1 (30). However, little is known about how sGC subunit GUCY1A2 was related to the above signaling pathway and further exploration is needed.

Conclusion
Our current research showed that GUCY1A2 mays serve as a potential prognostic biomarker of poor survival in gastric cancer. Calcium signaling pathway, MAPK signaling pathway, TGF-β signaling pathway and Wnt signaling pathway maybe the critical pathway regulated by GUCY1A2 in gastric cancer. However, GUCY1A2 expression level in this research was mRNA expression rather than protein expression and using mRNA to predict protein expression was not perfect. Further study is needed to validate the results and prove the biological in uence of GUCY1A2. Availability of data and materials

List Of Abbreviations
The datasets generated and analysed during the current study are available in the TCGA repository, http://cancergenome.nih.gov/.

Competing interests
The authors declare that they have no competing interests