Wei Zhang1,2,3†
Yutong Shi4†Shumeng Niu5Lintai Li1
Liewen Lin1
Xucan Gao1
Wanxia Cai1
Yumei Chen1Yafang Zhong1
Donge Tang1*
Min Tang4*
Yong Dai1*- 1Clinical Medical Research Center, Shenzhen People’s Hospital, The Second Clinical Medical College of Jinan University, Shenzhen, China
- 2South China Hospital, Health Science Center, Shenzhen University, Shenzhen, China
- 3Medical School, Nanchang Institute of Technology, Nanchang, China
- 4Key Laboratory of Diagnostic Medicine designated by the Chinese Ministry of Education, Chongqing Medical University, Chongqing, China
- 5Laboratory Department, Shanghai Hongkou Jiangwan Hospital, Shanghai, China
According to a recent report by GLOBOCAN, colorectal cancer is the third most common and second most deadly cancer in 2020. In our previous proteomic study, we found that the expression of GSTM2 in colon tissues was significantly lower than that in para-cancer tissues, and its lower expression was associated with reduced overall survival rate of patients, suggesting that this gene might play a role in the occurrence of colon cancer. As a member of the detoxifying enzyme family, GSTM2 is likely to play an important role in the initiation of tumors. Whereas, the functions of GSTM2 in colon cancer are barely known. In this study, using the RNA-Seq datasets of colon cancer patients from public database (ntumor = 457, nnormal = 41), we confirmed the reduced expression of GSTM2 and its prognostic value in colon cancer. Furthermore, we used our own Chinese cohort (ntumor = 100, nnormal = 72) verified the lower GSTM2 expression in colon cancer, and also its effects on patient prognosis. Subsequently, we uncovered two potential reasons for the lower expression of GSTM2 in colon cancer tissues, including the deep deletion of GSTM2 on genome, and the up-regulation of RAD21 or SP1. Moreover, we disclosed that GSTM2 might be involved in several immune-related pathways in colon cancer, such as chemokine signaling and leukocyte transendothelial migration. Finally, we revealed that the GSTM2 expression was closely related to the immune-related scores of colon cancer and the infiltration ratios of various immune cells, suggesting that GSTM2 might regulate the development of colon cancer by modulating immune microenvironment. In conclusion, we uncovered the prognostic value of GSTM2 based on the public data and our own data, revealed its potential regulatory role in tumor immune microenvironment, and disclosed the probable reasons for its lower expression in colon cancer. The findings of our study provide a potential prognostic biomarker and drug target for clinical diagnosis and treatment of colon cancer.
Introduction
According to the latest cancer statistics, the incidence and mortality rate of colorectal cancer is the third and second highest in the world, respectively. Globally, about 1.93 million people have been diagnosed as colorectal cancer, with even more than 930,000 deaths (1). In addition to surgical treatments, radiation, and chemotherapy, the treatments with immune checkpoint inhibitors (ICIs) have achieved success in the therapy of colorectal cancer in recent ten years (2). Many studies have revealed that the use of PD-1 and CTLA-4 inhibitors is effective in advanced colon cancer patients with mismatch repair defects or highly unstable microsatellites. However, most colon cancer patients are in a non-mismatched repair state, or a microsatellite stable state, thereby may not benefit from the ICI treatment (3). Therefore, it is necessary to continue exploring novel immune drug targets for colon cancer treatment.
Recently, numerous studies have revealed that tumorigenesis is strongly associated with the composition of immune microenvironment (TIM), such as immune cells and relevant cytokines (4). There is growing evidence that the TIM is valuable in predicting patient prognosis in multiple cancers (5). Different types of immune cell infiltration may provide a favorable or unfavorable environment for cancer development, and also being one of the key factors for immunotherapy responses.
The glutathione-S-transferase family (GST) is an important antioxidant enzyme family in living organisms, and GSTM belongs to the Mu subfamily of the GSTs (6). The GSTMs encode a series of detoxification enzymes, participating in detoxifying electrophilic compounds through coupling glutathione (carcinogens, environmental toxins, drugs, and oxidative stress products, etc.) (7). Some studies have demonstrated that GSTM1, GSTA4, GSTO1, in addition to their protection against oxidative stress, participate in the processes including immune defense and antiviral response (8–10). However, there are few studies reporting on the GSTM’s functions towards the TIM, and so as their actions on the tumor prognosis.
In our previous proteomic study (ntumor = 8, nnormal = 8), we observed an extreme decrease in the expression of GSTM2 in colon tumor tissues versus normal adjacent tissues (11). Meanwhile, the data set from Clinical Proteomic Tumor Analysis Consortium (CPTAC) further confirmed our finding, and revealed that the reduced GSTM2 protein expression was correlated with lower patient survival rate (Wei Zhang, et al). Therefore, we hypothesized that GSTM2 might play a role in the tumorigenesis of colon cancer. A previous study has proclaimed that the butyrate can induce the GSTM2 expression in colon cancer, implying that GSTM2 may increase the detoxification ability of colon mucosa and play a protective role (12). Meanwhile, the mRNA expression of GSTM2 has also been found to be associated with prognosis of colon cancer patients based on a bioinformatics analysis (13), but its protein expression has not been verified in a self-built cohort, either its potential functions and mechanisms underlying tumorigenesis not explored in depth.
In this study, we studied the mRNA and protein expression of GSTM2 in colon cancer tissues and normal colon tissues using the public data sets. Besides, we analyzed the effects of GSTM2 on patient survival based on three independent data sets. Further, we revealed the possible pathways by which GSTM2 affected the colon cancer tumorigenesis, and investigated the correlation between the GSTM2 expression and the infiltration ratios of immune cells. Moreover, we explored the probable reasons for the decreasing expression of GSTM2 in colon cancer. Finally, we validated the GSTM2 expression and its action on patient survival based on a self-built Chinese cohort (ntumor = 100, nnormal = 72).
Materials and methods
Expression spectrum
Tumor Immune Estimation Resource (TIMER) is a comprehensive database containing RNA-Sequencing data sets from The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx), used for the analysis of the gene expression and immune cell infiltration (14). In our study, the “Diff Exp” module was used to investigate the GSTM2 expression in multiple types of cancers and normal tissues.
The relationship between GSTM2 expression and clinical features
UALCAN is a comprehensive interactive database providing in-depth analysis of data sets from TCGA and Clinical Proteomic Tumor Analysis Consortium (CPTAC) (16). In this study, the mRNA and protein expression of GSTM2 was analyzed using the UALCAN database. Besides, the relationship between the GSTM2 expression and the clinicopathological parameters (sex, age, ethnicity, tumor stage, body weight, histological typing, lymph node metastasis, and TP53 mutation) in colon cancer patients were also explored using the UALCAN database.
Survival analysis
Prognoscan is an integrated database including the publicly available chip data that is used to analyze the actions of gene expression on patient prognosis (17). In this study, the effects of GSTM2 expression on patient survival was investigated using the microarray data sets of colorectal cancer in the Prognoscan database.
Immunological analysis
XCELL, TIMER, CIBERSORT, MCPCOUTER, QUANTISEQ, and EPIC are the algorithms that estimate the proportion of immune cells in tumors based on the gene expression matrices (18). This study used these six algorithms to analyze the correlation between the GSTM2 expression and the ratios of different infiltrating immune cells. Using the “Immune” module on the TIMER 2.0 website, the link between GSTM2 expression and immune infiltration ratios was investigated. SangerBox was used to study the relationship between GSTM2 expression and microsatellite instability (MSI).
Genetic mutations
The cBioportal database is an open platform for exploring gene mutations in multiple cancers (19). In this study, five independent studies were selected to search the GSTM2 mutations in colon cancer patients. The “cancer types summary” module was utilized to explore the frequency of GSTM2 mutations in colorectal cancer. Through the “mRNA vs. Dx” in “Plots” module, the relationship between the GSTM2 expression and different types of GSTM2 mutations in colon cancer was analyzed.
Transcription factors
The hTFtarget database is a large-scale database collecting most human transcription factors and their experimentally confirmed target genes in multiple organs (20). In our study, the transcription factors of GSTM2 in colon tissue was searched using the hTFtarget database.
Correlation analysis
Gene Expression Profiling Interactive Analysis (GEPIA) provides online interaction and customization analyses, including expression profiling of tumor and normal tissues (15). The “Correlation Analysis” function in the GEPIA database was used to investigate the expression correlation between two genes. A p-value < 0.05 was regarded as statistically significant. In this study, the data sets of colon adenocarcinoma (COAD) and normal samples from TCGA and the data sets of sigmoid colon normal tissue samples from GTEx were used.
Gene set enrichment analysis (GSEA)
The co-expressed mRNAs of GSTM2 in colorectal cancer were obtained from the Linkedomics database, and the “LinkInterpreter” function in this database was used to investigate the potential relevant pathways of GSTM2 in colon cancer based on GSEA.
Single sample gene set enrichment analysis (ssGSEA)
For colon cancer, RNAseq data (level 3) and clinical data were retrieved from TCGA database. The R program GSVA package was used for analysis, and the method “ssGSEA” was chosen. Finally, spearman correlation analysis was applied to examine the correlation between GSTM2 expression and the pathway scores. Statistics were judged significant at p value 0.05.
Patients
The paraffin blocks of colon cancer tissues were provided by the National Human Genetic Resources Sharing Service Platform. The inclusion criteria were as follows: age ≤ 90 years old, histologically-confirmed colon cancer, no serious major organ dysfunction, and no chemotherapy or radiotherapy. The exclusion criteria were: age ≥ 90 years, severe major organ dysfunction, prior chemo- or radiation treatment, or a family history of colon cancer. All patients were informed and had signed an informed consent form. This work was carried out under the supervision of the Ethics Committee of Taizhou Hospital and met the requirements of the World Medical Association for the execution of human experiments.
Immunohistochemistry assay
The tissue sections were stayed at room temperature for 60 minutes and then soaked in xylene for 10 minutes. Then, the soaking step was repeated once. Next, the chips were dewaxed in 100%, 95%, and 75% ethanol, respectively, and heated at 95°C for 15 minutes in the 0.01 M lycium buffer (pH = 6.0). This heating step was repeated 1-2 times. The chip was then cleaned three times with the PBS solution for 5 minutes each. Subsequently, 3% H2O2 was dripped onto the chip, and the chip was stayed at room temperature for 10 minutes, and then washed three times with the PBS solution for 5 minutes each. Subsequently, the chip was incubated with the blocking solution for 20 minutes at room temperature. The primary antibody of GSTM2 (Invitrogen, MA5-29311, 1:1500) was then added to the chip and the chip was incubated at 4°C for 12 hours. Then, the chip was cleaned three times with the PBS solution for 5 minutes each. The incubation solution with secondary antibody was then added to the chip, and the chip was incubated at 37°C for one hour. Then, the chip was cleaned three times with the PBS solution for 5 minutes each. Next, the diaminobenzidine (DABs) was used to color the chips, and the staining intensity was observed under a microscope. Subsequently, the chip was rinsed with the tap water for 10 minutes, and incubated with hematoxylin for 2 minutes, and differentiated with hydrochloride alcohol. After washed with the tap water for 10 minutes, the chip was then dehydrated, transparent, and sealed. Next, the stained chip was scanned using the Aperio Scanner (LEICA, Aperio XT), and the staining intensity and the proportion of the stained cells were assessed by two pathologists. 0, 1, 2, and 3 were used to represent the staining strength (0 represented negative, 1 represented weak, 2 represented moderate, and 3 represented strong). The proportion of the stained cells was evaluated based on the percentage of the positively stained area relative to the entire cancer region. Finally, the numbers representing the staining intensity and the proportion were multiplied to obtain the staining score.
Statistical analysis
GraphPad Prism 6.0 was used for the statistical analysis. The GSTM2 expression in cancerous and normal tissues were assessed using the unpaired T-test. The survival analysis was performed using the “Survival” model of GraphPad Prism 6.0. The “Kaplan Meier Survival Curve” module in Sangerbox was used to determine the best cutoff value of GSTM2 expression in clinical samples. The Cox regression analysis was performed using SPSS Statistics 27. A p value < 0.05 was considered statistically significant, and the confidence interval was set as 95% confidence.
Results
GSTM2 is decreasingly expressed in colon cancer tissues versus normal tissues
To explore the potential role of GSTM2 in tumorigenesis, we first used the TIMER database to assess the mRNA expression of GSTM2 in various human tumors. The results showed that GSTM2 was differentially expressed in tumor tissues versus normal tissues in 15 tumors, showing a down-regulated expression in 13 cancers, including COAD, breast invasive carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), kidney chromophobe (KICH), kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), stomach adenocarcinoma (STAD), thyroid carcinoma (THCA), and uterine corpus endometrial carcinoma (UCEC) (Figure 1B). In addition, after analyzing the data sets from TCGA and GTEx, we found that the GSTM2 expression was down-regulated in 23 tumors and only up-regulated in two tumors, and most of the results were consistent with the results of Figure 1B (Figure 1C). Furthermore, we investigated the GSTM2 protein level using the data set from CPTAC and discovered that the GSTM2 protein expression in colon cancer tissues was also reduced compared to normal colon tissues (n tumor = 97, n normal = 100) (Figure 1D).
Figure 1 GSTM2 is Decreasingly Expressed in Colon Cancer Tissues versus Normal Tissues. (A) Flow chart of this research. GSTM2 mRNA expression in different types of cancer was studied using (B) the TIMER database, and (C) the TCGA database and GTEx database. The protein expression of GSTM2 in colon cancer tissues versus normal tissues was investigated using (D) the UALCAN database. *p < 0.05, **p < 0.01, and ***p < 0.001.
The Correlation between the GSTM2 expression and various clinical features
Subsequently, we used the UCLCAN database to investigated the GSTM2 mRNA expression in different patient groups (Supplementary Table S1). In Figure 2A, the GSTM2 expression was significantly decreased in colon cancer tissues of patients at stages 1, 2, 3, and 4 versus normal controls, while there was no significant difference between different stages. In terms of race, the GSTM2 expression was reduced in tumor tissues of Caucasian, African-american, and Asian patients versus normal controls, but there were no significant differences between distinct races (Figure 2B). Regarding to gender, there was no significant difference between male and female patients (Figure 2C). As for age, the GSTM2 expression was decreased in each age group versus normal group, and the GSTM2 expression was significantly higher in patients aged 21-40 years than 41-60 years, 61-80 years, 81-100 years (Figure 2D). In terms of weight, the GSTM2 expression was down-regulated in tumorous tissues from patients of different weights compared to normal tissues, while it was higher in the obese group than in the normal weight group (Figure 2E). As for histological subtypes, the GSTM2 expression was decreased in tumorous tissues of patients with different subtypes of colon cancer versus normal tissues, but there was no statistical difference between each two histological subtypes (Figure 2F). As for lymph node metastases, the GSTM2 expression in cancerous tissues from patients at N0, N1, and N2 stages was lower than in normal tissues. Still, there was no statistical difference in expression between the three groups (Figure 2G). Regarding to TP53 mutations, the GSTM2 was decreasingly expressed in tumor tissues from patients with wild type TP53 and mutated TP53 versus normal tissues, but the TP53 mutations did not affect the GSTM2 expression (Figure 2H).
Figure 2 The Correlation between the GSTM2 Expression and Various Clinical Features. Analysis revealed the correlation between GSTM2 expression and (A) tumor stages, (B) race, (C) gender, (D) age, (E) weight, (F) histological subtypes, (G) metastasis, and (H) TP53 mutation. N0 indicated no regional lymph node metastasis, N1 indicated 1 to 3 axillary lymph node metastasis, N2 indicated 4 to 9 cases of axillary lymph node metastasis. ns: not statistically significant, *p < 0.05, and ***p < 0.001.
At the same time, we also compared the GSTM2 protein expression in different patient subgroups based on data sets from the UALCAN database. Shown in Supplementary Figure S1A, the GSTM2 protein expression in tumor tissues from patients at different tumor stages was lower than normal tissue. Still, there was no significant difference in GSTM2 protein expression between each stage. In terms of races, the expression of GSTM2 protein in cancerous tissues of Caucasian, African-american, and Asian patients was lower than in normal group. The statistical results showed that GSTM2 protein expression in Asian was higher than in Caucasian (Supplementary Figure S1B). As for gender, the GSTM2 protein expression in tumor tissues from both male and female colon cancer patients was lower than in normal tissues, while there was no statistical difference in GSTM2 protein expression between male and female patients (Supplementary Figure S1C). In terms of age, the GSTM2 protein expression was decreasingly expressed in cancerous tissues from patients of different ages than in control group, but there was no statistical different between each subgroup (Supplementary Figure S1D). Regarding to weight, there were no differences in protein expression between distinct groups (Supplementary Figure S1E). As for sub-classification, the GSTM2 protein was decreasingly expressed in mucinous and non-mucinous colon cancer tissues versus normal tissues, but there was no statistical difference between each two sub-classifications (Supplementary Figure S1F).
The lower GSTM2 expression predicts worse prognosis of colon cancer patients
To clarify the functional role and prognostic value of GSTM2 in colon cancer, we further analyzed the relationship between the GSTM2 expression and prognosis of colon cancer patients using four data sets from Gene Expression Omnibus (GEO) collected by the Prognoscan database. Based on the results of three data sets named GSE17536, GSE17537, GSE12945, we demonstrated that colon cancer patients with lower expression of GSTM2 had worse overall survival rates (OS), disease-free survival (DFS) and disease specific survival rates (DSS) (Figures 3A–H). Also, the survival analysis revealed colon cancer patients with low GSTM2 protein expression had a shorter overall survival (Figure 3I). Therefore, the GSTM2 mRNA expression was potentially associated with the survival of colon cancer patients, and GSTM2 might be served as a prognostic marker for colon cancer. To explore the prognosis of GSTM2 mutant group and non-mutant group, we used the data sets from cBioportal database for analysis (nmutant = 31, nnon-mutant = 3376). However, the results demonstrated that there was no statistical difference in survival rate between GSTM2 mutant and non-mutant group (p = 0.127) (Supplementary Figure S2A). We hypothesized that the negative result was due to a lack of samples in the mutant group, or that the GSTM2 mutants were not functional, thus having no effect on patient survival rates.
Figure 3 The Lower GSTM2 Expression Predicts Worse Prognosis of Colon Cancer Patients. (A–C) The correlation between the GSTM2 expression and OS of colon cancer patients using the GSE17536, GSE17537, GSE12945 datasets. (D–F) The correlation between the GSTM2 mRNA expression and DSS of colon cancer patients using the GSE17537, GSE17536, GSE12945 datasets. (G, H) The correlation between the GSTM2 mRNA expression and DFS of colon cancer patients using the GSE17536, GSE17537 datasets. (I) The correlation between the GSTM2 protein expression and OS of colon cancer patients using the dataset from CPTAC.
Two potential reasons for the down-regulation of GSTM2 in human colon cancer
To explore the possible reasons for the down-regulation of GSTM2 in colon cancer, we first studied the DNA mutations of GSTM2. By analyzing the data sets from the cBioportal database (n = 4488), we found GSTM2 had a moderate mutation frequency in colon cancer patients (1%-4%), and about 50% of these mutations were deep deletion (Figure 4A). And, the GSTM2 mRNA expression in cancerous tissues from patients with deep deletion was lower than the median value of all types of mutations (Figure 4B). This result indicated that deep deletion on DNA might be one of the reasons for the down-regulation of GSTM2 expression in colon cancer tissues.
Figure 4 Two Potential Reasons for the Down-regulation of GSTM2 in Human Colon Cancer. (A) The mutation type and frequency of GSTM2 in colorectal cancer. (B) The GSTM2 expression in tumor tissues from colorectal cancer patients with different GSTM2 mutations. (C–E) The protein expression correlation between GSTM2 and CTCF, RAD21, SP1 in colon cancer. (F–H) The protein expression of CTCF, RAD21, SP1 in colon cancer tissues versus normal tissues. ***p < 0.001.
Additionally, we investigated the expression changes of the potential transcription factors of GSTM2 in colon cancer tissues versus normal tissues. We first downloaded 190 experimentally confirmed transcription factors of GSTM2 from the hTFtarget database (Supplementary Table S2), of which CTCF, RAD21, SP1 were identified transcription factors in colon tissues. As we knew that the transcription factors always had a correlated expression with their downstream genes, thus we used the GEPIA database to explore the expression correlation between GSTM2 and CTCF, RAD21, SP1 in colon cancer. The results showed that there were significantly negative expression correlations between the GSTM2 expression with the three transcription factors (Figures 4C–E), which suggested CTCF, RAD21, SP1 as potential transcription factors of GSTM2 in colon cancer. Next, we used the UALCAN database to analyze the expression of these three transcription factors in colon cancer. Shown in Figure 4F, the CTCF protein was expressed lower in tumor tissues than normal adjacent tissues, while RAD21, SP1 protein were expressed more in cancerous tissues than in normal tissues (Figures 4G, H). Therefore, of the three genes, the expression of RAD21 and SP1 in colon cancer matched the results of correlation analysis, so that RAD21 and SP1 might be the potential upstream transcription factors of GSTM2 in colon cancer. Consequently, the up-regulation of RAD21 and SP1 might be the second probable reason for the down-regulated expression of GSTM2 in colon cancer.
The functions of GSTM2 in colon cancer might be associated with immune modulations
To understand the underlying mechanism by which GSTM2 affected the development of colon cancer, we first collected a total of 5538 co-expressing mRNAs of GSTM2 from the linkedomics database (Supplementary Table S3). The first 50 genes positively and negatively correlated with GSTM2 were shown in Figures 5A, B. Then, we performed GSEA on GSTM2 to search the pathways GSTM2’s probably participating in, and the enrichment results showed that GSTM2 were most associated with immune related biological processes and pathways, such as chemokine signaling, leukocyte transendothelial migration, IgA production, Th1 and Th2 cell differentiation, primary immunodeficiency, adaptive immune response, B cell activation, defense response, immunoglobulin binding, and antigen binding (Figures 5C–F). These results suggested that GSTM2 might affect the tumorigenesis of colon cancer via regulating the immune microenvironment. Furthermore, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of GSTM2’s co-expressed genes which were downloaded from Linkedomics database (|correlation coefficient| > 0.3, false discovery rate (FDR) < 0.05). As shown in Supplementary Figures S2B–E, GSTM2 co-expressed genes were not only involved in tumor-related pathways, but also in immune-related pathways, such as chemokine signaling pathway and cytokine binding.
Figure 5 The Functions of GSTM2 in Colon Cancer might be Associated with Immune Modulations. (A) Heat maps of the top 50 genes positively correlated with GSTM2 expression. (B) Heat maps of the top 50 genes negatively correlated with GSTM2 expression. (C–F) GSEA analysis showing the positively and negatively correlated pathways of GSTM2. ns: not statistically significant, *p < 0.05, and ***p < 0.001.
To learn more about GSTM2’s potential function in carcinogenesis, we analyzed the correlation between GSTM2 expression and the activities of key pathways (or gene signatures) involved in tumorigenesis. The results showed that GSTM2 expression was associated with the activities of signatures including cellular response to hypoxia, tumor proliferation signature, apoptosis, DNA repair, G2M checkpoint, MYC targets, TGFB, DNA replication, and degradation of ECM (Supplementary Figure S3).
The GSTM2 expression was potentially related to immune infiltration in colon cancer
To confirm the correlation between the GSTM2 expression and immune microenvironment, we used six distinct algorithms to assess the correlation between the GSTM2 expression and immune cell infiltration in human colon cancer. As shown in Figure 6A, there was a significant correlation between the GSTM2 mRNA expression and stromal score, microenvironment score, and immunity score in colon cancer (p < 0.001). Meanwhile, the XCELL results revealed that the GSTM2 expression was positively correlated with the ratios of 21 types of immune cells, such as regulatory T cells (Tregs), T cell CD8+, monocyte, macrophage, endothelial cell, B cell, and was negatively correlated with the ratios of four types of immune cells, including T cell natural killer cell (NK), T cell CD4+ T helper 2 cell (Th2), T cell CD4+ T helper 1 cell (Th1), and common lymphoid progenitor. Besides, the CIBERSORT results demonstrated that the GSTM2 expression was positively correlated with the levels of T cell CD4+ naïve, mast cell activated, macrophage M2, B cell plasma, B cell memory infiltration, and negatively correlated with mast cell resting (Figure 6B). The results of the TIMER algorithm showed that the GSTM2 expression was highly positively correlated with the infiltration degrees of T cell CD8+, T cell CD4+, neutrophil, myeloid dendritic cell, macrophage, B cell, and other cells (Figure 6C). The MCPCOUTER results uncovered that the GSTM2 expression was significantly positively correlated with the infiltrating ratios of T cell CD8+, T cell, neutrophil, myeloid dendritic cell, monocyte, macrophage, endothelial cell, and B cell (Figure 6D). The QUANTISEQ results showed that the GSTM2 expression was positively correlated with the infiltration levels of Tregs, T cell CD8+, neutrophil, NK cell, myeloid dendritic cell, macrophage M2, and B cell (Figure 6E). Similarly, the EPIC algorithm results showed that the GSTM2 expression was positively correlated with the infiltration degrees of T cell CD8+, T cell CD4+, macrophage, endothelial cell, and B cell (Figure 6F). Taken together, these results disclosed that the GSTM2 mRNA expression was correlated to the immune scores, and universally associated with the infiltration levels of five types of immune cells, including B cell, T cell CD8+, T cell CD4+, macrophage, dendritic cell, which further indicated that GSTM2 was a potentially regulatory gene for immune microenvironment of colon cancer. In colon cancer, MSI was highly associated with immune cell infiltration and patients’ clinical response to immunotherapy (21, 22), thus we evaluated the correlation between GSTM2 expression and MSI in various cancers. As a result, GSTM2 expression was negatively associated with MSI in colon cancer (Figure 6G). This result supports the conclusion of our study that GSTM2 may be involved in the regulation of immune microenvironment in colon cancer, and also indicates that GSTM2 expression may have an impact on immunotherapy. After that, we used TIMER database to analyze the relationship between GSTM2 expression and the infiltration levels of six kinds of classical immune cells in TCGA colon cancer patients. As shown in Figure 6H, GSTM2 expression was significantly positively correlated with the infiltration degrees of B cells, CD4+T cells, macrophages, neutrophils and dendritic cells, which was consistent with Figures 6A–F results calculated by the six algorithms.
Figure 6 The GSTM2 Expression in Colon Cancer was Potentially Related to Immune Infiltration. (A–F) The correlation between the GSTM2 expression and infiltration ratios of immune cells in colon cancer as analyzed separately using the XCELL, CIBERSORT, TIMER, MCPCOUNTER, QUANTISEQ, and EPIC algorithms. (G) The correlation between GSTM2 expression and MSI in various tumors. (H) The correlations between GSTM2 expression and the infiltration levels of six classic types of immune cells in colon cancer. *p < 0.05, **p < 0.01, ***p < 0.001.
The immunohistochemical results confirmed that GSTM2 was decreasingly expressed in tumor tissues versus normal adjacent tissues based on our own self-built Chinese cohort
To confirm the GSTM2 expression in colon cancer, we first analyzed the GSTM2 protein expression in human colon cancer tissues and normal adjacent tissues using our own proteomic data (n = 8, samples of the same stage were pooled into one sample) (23). As a result, the expression of GSTM2 protein was significantly lower in cancerous tissues compared to their paired non-tumor tissues (Figure 7A). Next, we further detected the GSTM2 protein expression in our self-built Chinese cohort (ntumor = 100, nnormal = 72) using the immunohistochemical assay (Supplementary Data), and the patient clinical information were detailed in Table 1. When solely calculating the colonic cells, we observed the staining score of GSTM2 (stain intensity timed the percentage of stained cells) in cancerous cells was significantly lower than normal adjacent cells (p = 0.0015) (Figure 7B). Since GSTM2 is potentially related to the regulation of immune microenvironment in colon cancer, and the previous studies have also shown that the expression changes of functional genes in immune cells may affect the tumorigenesis, we then analyzed the expression changes of GSTM2 in the infiltrating lymphocytes within tumor tissues. As a result, the staining intensity of GSTM2 in the lymphocytes within the cancerous tissues was significantly lower than the lymphocytes within the normal tissues (p = 0.0001) (Figure 7C). The Figures 7D, E were the representative pictures of immunohistochemical detection of GSTM2 in cancer cells or lymphocytes within the tumor tissues and normal adjacent tissues. Furthermore, we also detected GSTM2 level in different grades, T stages, N stages and M stages of colon cancer, and found that the GSTM2 level in tumor cells and lymphocytes decreased with the increase of M stages (Supplementary Figure S4), suggesting that GSTM2 might be associated with metastasis.
Figure 7 The Immunohistochemical Results Disclosed that GSTM2 was Decreasingly Expressed in Tumor Tissues versus Normal tissues Based on our Self-built Chinese Cohort. (A) The GSTM2 expression was evaluated in 4 paired colon cancer tissues and adjacent normal tissues using our own proteomic data. (B) The GSTM2 staining score in colon cancer versus normal tissue. (C) The GSTM2 staining score in lymphocyte within colon cancer tissue versus normal tissue. (D, E) Representative images of immunohistochemical staining of tumor cell and lymphocyte within tumor tissues and normal tissues. The arrows identify representative lymphocytes or tumor cells. (F, G) The correlation between the GSTM2 protein expression with the OS of colon cancer patients based on our self-built Chinese cohort. **p < 0.01 and **** p < 0.0001.
Table 1 The characteristics of 104 colorectal cancer patients in our study.
Also, we analyzed GSTM2 mRNA expression in individual cells using the data set of Li’s work published in Nature Genetics (24). The results revealed that compared with certain cells from normal tissues, GSTM2 expression was statistically reduced in stem cells (p = 0.008), and had a decreasing trend in fibroblasts (p = 0.311) and B cells (p = 0.13) from colorectal cancer tissues (Supplementary Table S4). As is well-known, stem cells are a kind of primitive cells with the self-renewal ability and multidirectional differentiation potential, which is the origin of cell proliferation. Tumor stem cells have the properties of initiating tumors, promoting metastasis, enhancing drug resistance (25). The lower expression of GSTM2 in stem cells further illustrates this gene’s potential importance and functions in tumorigenesis of colon cancer.
At the same time, we also evaluated the effects of GSTM2 protein expression on patient survival after grouping the patients into two groups based on the best cutoff value of GSTM2 expression (the best cutoff value of GSTM2 expression in tumor cell group was 1, and in lymphocyte group was 0.075). When we analyzed the GSTM2 expression in tumor cells, we found that the patients with lower GSTM2 immunohistochemical scores frequently had poorer OS (p = 0.12, HR = 1.487, 95%CI (0.9877 - 2.4707) (Figure 7F). Similarly, when we analyzed the GSTM2 expression in infiltrating lymphocytes, the patients with lower GSTM2 immunohistochemical scores also had worse OS (p < 0.001, HR = 2.607, 95%CI (1.847 - 6.184) (Figure 7G). We also performed univariate and multivariate Cox regression analyses of GSTM2 protein expression and clinical characteristics, and the results showed that the p value of GSTM2 protein level in infiltrating lymphocytes within tumor tissues was less than 0.05, which indicated that GSTM2 was potentially an independent prognosis factor for colon cancer (Table 2).
Table 2 Univariate and multivariate analyses identified GSTM2 level in lymphocytes as an independent prognostic factor for colon cancer.
Additionally, we hoped to find out which medicines GSTM2 expression was linked to their resistance. We downloaded drug-resistance data sets from Genomics of Drug Sensitivity in Cancer (GDSC) database, exploring drug resistance relevant to GSTM2 expression, but no statistical result was found (correlation coefficient >0.2 or <-0.2, p<0.05) (Supplementary Table S5).
To sum up, the above results confirmed the findings we discovered based on the public data, demonstrating that GSTM2 were lower expressed in colon cancer tissues versus normal tissues and its expression were potentially relevant to the patient prognosis (OS).
The mechanism underlying GSTM2 in colon cancer might be associated with the functions of PCDH17 mutations
In addition, we analyzed the differences in gene mutations between GSTM2-high-expression and GSTM2-low-expression groups using the data sets from TCGA database. To further clarify the mechanism of GSTM2 in tumorigenesis, we analyzed genome-wide mutation differences in GSTM2-high-expression and low-expression colon cancer patients. We found that the mutation frequency of ERBB3 were drastically increased in GSTM2 high-expressed patients, while the mutations of MYO15A, LOXL4, DNAH3, CPS1, PCDH17, L3MBTL4, PLCG1, SAFB, OTOG, CNGA3, and DDX60 were significantly elevated in GSTM2-low-expression patients (Figures 8A–C). Among these mutations, PCDH17 has been identified as a driver gene that induces colon cancer tumorigenesis (26). Therefore, we speculated that the functions of GSTM2 in colon cancer carcinogenesis might be related to PCDH17 mutation.
Figure 8 The Mechanism Underlying GSTM2 in Colon Cancer might be Associated with the Functions of PCDH17 mutations. (A) The forest map showing differential mutant genes between GSTM2 high and low expression groups. (B, C) The waterfall maps showing genes’ mutation frequencies in each group. **p < 0.01 and ***p < 0.001.
Discussion
The GSTM family encodes a series of detoxifying enzymes that promote the metabolism of toxic substances, thereby may play an important role in tumor suppression. Some reports have revealed that the GSTM enzymes also participate in the metabolism of tumor chemotherapeutic drugs, and protection of organelle or cells from stress injures (27). GSTM2 encodes a protein with GST structure domain that plays important roles in various biological activities (7). In addition to detoxification, GSTM2 is also a transporter that can carry lipophilic molecules such as cholic acid, bilirubin, hormones, thyroid hormones, steroids, and a series of exogenous compounds and engage in the production of hormones, such as prostaglandins (28). GSTM2 has been discovered to participate in ASK1, JNK, and P38-MAPK pathways in organisms (29). Several studies have proclaimed that the GSTM1 expression is closely related to the poorer prognosis of colorectal cancer, breast cancer, and lung cancer (30–32). However, the available researches on GSTM2 are relatively insufficient.
Guo E, et al. (13)have demonstrated that the low expression of GSTM2 in colon cancer is related to the “better” prognosis of patients using a bioinformatics analysis. Wonderingly, the lower expression of functional genes is usually associated with “poorer” prognosis in general, thus we speculated that some factors, such as sample size, race, tumor subtype, and age, might have influenced the previous result. Using three independent data sets and a data set of self-constructed Chinese cohort, our study revealed that the lower expression of GSTM2 was associated with the poorer patient outcomes. Compared with the previous study, our study has the following innovative, expanded and in-depth discoveries: First, we not only revealed the expression of GSTM2 and its relationship with prognosis through public data, but also verified the results through our self-built Chinese cohort (n tumor = 104, n normal = 74). We confirmed the GSTM2 expression in Chinese patients and its effect on the occurrence of colon cancer. Second, we uncovered the possibility that GSTM2 might regulate tumorigenesis via modulating the immune microenvironment. To our best knowledge, this is the first time that the potential role of GSTM2 in tumor immunity has been disclosed. Third, we identified two potential causes of down-regulation of GSTM2 expression in colon cancer.
The tumor microenvironment of colon cancer plays an extremely important role in tumorigenesis. Our results revealed that the GSTM2 expression was probably positively associated with the infiltration ratios of most immune cells in colon cancer, such as CD8+ T cells. It has been reported that CD8+ T cells can kill tumor cells, affecting the prognosis of colon cancer and immunotherapy responses (33, 34). Therefore, the patients with low GSTM2 expression might have low infiltrating ratios of T cell CD8+, thereby the tumor initiation may be promoted.
We predicted the potential upstream transcription factors of GSTM2 in colon cancer tissues. First, we discovered three GSTM2 transcription factors in colon tissues using the TFtarget database, and we suspected that in colon cancer tissues, these transcription factors might probably launch GSTM2’s transcription. Because transcription factors are frequently connected with the expression of target genes, we used transcriptome data of colon cancer tissues to determine which transcription factors were more likely to co-express with GSTM2. Eventually, we identified RAD21 and SP1 as potential transcription factors of GSTM2 in colon cancer. RAD21 encodes a subunit of cohesin, and mainly involved in modifying the adhesion and separation of sister chromatids, thus regulating gene transcription and DNA damage repairing. Numerous studies have revealed that RAD21 is over expressed in a variety of tumors, such as colorectal cancer, non small cell lung cancer, and breast cancer (35–37). Besides, It has also been proven that the higher RAD21 expression in tumor tissues is associated with shorter DSS of patients with colorectal cancer. Also, the higher expression of RAD21 in patients promotes the resistance of chemoradiotherapy and chemotherapy (36). All in all, our results are consistent with the previous studies, and the function of RAD21 in colon cancer further indicates that GSTM2 may be involved in tumor formation as a suppressor.
Conclusion
To sum up, our findings revealed that GSTM2 was decreasingly expressed in colon cancer tissues versus normal adjacent tissues, and its lower expression was associated with poorer prognosis of patients. Furthermore, we uncovered the potential immunological functions of GSTM2 in colon cancer. Our findings here suggest that GSTM2 may serve as a prognostic biomarker for colon cancer, and as a potential drug target for immunotherapy.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Taizhou Hospital. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Study concept and design: WZ, DT. Acquisition of data: WZ, LinL, YZ. Analysis and interpretation of data: YS, LieL. Statistical analysis: WZ, YC, SN, WC. Drafting of the manuscript: WZ, YS. Critical revision and final approval of the manuscript: all authors. Obtained funding: MT, YD, WZ. Study supervision: XG, MT. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (No. 82003172), the Postdoctoral Science Foundation of China (No. 2020M673065), the Guangdong Basic and Applied Basic Research Foundation (No. 2019A1515111138, No.2021A1515111071).
Acknowledgments
The authors would like to thank Taizhou Hospital for providing patients’ paraffin blocks, and Shenzhen People’s Hospital for providing fund support for the determination of experimental data.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.881906/full#supplementary-material
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2. Ganesh K, Stadler ZK, Cercek A, Mendelsohn RB, Shia J, Segal N H, et al. Immunotherapy in colorectal cancer: Rationale, challenges and potential. Nat Rev Gastroenterol Hepatol (2019) 16(6):361–75. doi: 10.1038/s41575-019-0126-x
3. Chalabi M, Fanchi LF, Dijkstra KK, Van den Berg JG, Aalbers AG, Sikorska K, et al. Neoadjuvant immunotherapy leads to pathological responses in MMR-proficient and MMR-deficient early-stage colon cancers. Nat Med (2020) 26(4):566–76. doi: 10.1038/s41591-020-0805-8
4. Binnewies M, Roberts EW, Kersten K, Chan V, Fearon DF, Merad M, et al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat Med (2018) 24(5):541–50. doi: 10.1038/s41591-018-0014-x
5. Gajewski TF, Schreiber H, Fu YX. Innate and adaptive immune cells in the tumor microenvironment. Nat Immunol (2013) 14(10):1014–22. doi: 10.1038/ni.2703
6. Pearson WR, Vorachek WR, Xu SJ, Berger R, Hart I, Vannais D, et al. Identification of class-mu glutathione transferase genes GSTM1-GSTM5 on human chromosome 1p13. Am J Hum Genet (1993) 53(1):220–33. doi: 10.1016/0378-1119(93)90713-D
7. Nebert DW, Vasiliou V. Analysis of the glutathione s-transferase (GST) gene family. Hum Genomics (2004) 1(6):460–4. doi: 10.1186/1479-7364-1-6-460
8. Wang Y, Wang P, Zhang Y, Xu J, Li Z, Li Z, et al. Decreased expression of the host long-noncoding RNA-GM facilitates viral escape by inhibiting the kinase activity TBK1 via s-glutathionylation. Immunity (2020) 53(6):1168–81. doi: 10.1016/j.immuni.2020.11.010
9. Nadarajapillai K, Liyanage DS, Sellaththurai S, Jeong T, Lee S, Lee J, et al. Glutathione-s-transferase alpha-4 in hippocampus abdominalis (big-belly seahorse): Molecular characterization, antioxidant properties, and its potent immune response. Comp Biochem Physiol C Toxicol Pharmacol (2021) 240:108917. doi: 10.1016/j.cbpc.2020.108917
10. Udayantha H, Liyanage DS, Nadarajapillai K, Omeka WKM, Yang H, Jeong T, et al. Molecular characterization, immune and xenobiotic responses of glutathione s-transferase omega 1 from the big-belly seahorse: Novel insights into antiviral defense. Fish Shellfish Immunol (2021) 109:62–70. doi: 10.1016/j.fsi.2020.12.010
11. Zhang W, Tang D, Lin L, Fan T, Xia L, Cai W, et al. Integrative multiplatform-based molecular profiling of human colorectal cancer reveals proteogenomic alterations underlying mitochondrial inactivation. Am J Cancer Res (2021) 11(6):2893–910.
12. Ebert MN, Klinder A, Peters WH, Schäferhenrich A, Sendt W, Scheele J, et al. Expression of glutathione s-transferases (GSTs) in human colon cells and inducibility of GSTM2 by butyrate. Carcinogenesis (2003) 24(10):1637–44. doi: 10.1093/carcin/bgg122
13. Guo E, Wei H, Liao X, Wu L, Zeng X. Clinical significance and biological mechanisms of glutathione s-transferase mu gene family in colon adenocarcinoma. BMC Med Genet (2020) 21(1):130. doi: 10.1186/s12881-020-01066-2
14. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res (2017) 77(21):e108–10. doi: 10.1158/0008-5472.CAN-17-0307
15. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res (2017) 45(W1):W98–W102. doi: 10.1093/nar/gkx247
16. Chandrashekar DS, Bashel B, Balasubramanya S, Creighton C J, Ponce-Rodriguez I, Chakravarthi BVSK, et al. UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (2017) 19(8):649–58. doi: 10.1016/j.neo.2017.05.002
17. Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: A new database for meta-analysis of the prognostic value of genes. BMC Med Genomics (2009) 2:18. doi: 10.1186/1755-8794-2-18
18. Sturm G, Finotello F, Petitprez F, Zhang JD, Baumbach J, Fridman WH, et al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics (2019) 35(14):i436–45. doi: 10.1093/bioinformatics/btz363
19. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088
20. Zhang Q, Liu W, Zhang HM, Xie GY, Miao YR, Xia M, et al. hTFtarget: A comprehensive database for regulations of human transcription factors and their targets. Genomics Proteomics Bioinf (2020) 18(2):120–8. doi: 10.1016/j.gpb.2019.09.006
21. Lin A, Zhang J, Luo P. Crosstalk between the MSI status and tumor microenvironment in colorectal cancer. Front Immunol (2020) 11:2039. doi: 10.3389/fimmu.2020.02039
22. Laghi L, Negri F, Gaiani F, Cavalleri T, Grizzi F, De' Angelis GL, et al. Prognostic and predictive cross-roads of microsatellite instability and immune response to colon cancer. Int J Mol Sci (2020) 21(24):9680. doi: 10.3390/ijms21249680
23. Zhang W, Lin L, Xia L, Cai W, Dai W, Zou C, et al. Multi-omics analyses of human colorectal cancer revealed three mitochondrial genes potentially associated with poor outcomes of patients. J Transl Med (2021) 19(1):273. doi: 10.1186/s12967-021-02939-7
24. Li H, Courtois ET, Sengupta D, Tan Y, Chen KH, Goh JJL, et al. Reference component analysis of single-cell transcriptomes elucidates cellular heterogeneity in human colorectal tumors. Nat Genet (2017) 49(5):708–18. doi: 10.1038/ng.3818
25. Adorno-Cruz V, Kibria G, Liu X, Doherty M, Junk DJ, Guan D, et al. Cancer stem cells: targeting the roots of cancer, seeds of metastasis, and sources of therapy resistance. Cancer Res (2015) 75(6):924–9. doi: 10.1158/0008-5472.CAN-14-3225
26. Wolff RK, Hoffman MD, Wolff EC, Herrick JS, Sakoda LC, Samowitz WS, et al. Mutation analysis of adenomas and carcinomas of the colon: Early and late drivers. Genes Chromosomes Cancer (2018) 57(7):366–76. doi: 10.1002/gcc.22539
27. Pljesa-Ercegovac M, Savic-Radojevic A, Matic M, Coric V, Djukic T, Radic T, et al. Glutathione transferases: Potential targets to overcome chemoresistance in solid tumors. Int J Mol Sci (2018) 19(12):3785. doi: 10.3390/ijms19123785
28. Beuckmann CT, Fujimori K, Urade Y, Hayaishi O. Identification of mu-class glutathione transferases M2-2 and M3-3 as cytosolic prostaglandin e synthases in the human brain. Neurochem Res (2000) 25(5):733–8. doi: 10.1023/a:1007579507804
29. Jin Y, Tan Y, Zhao P, Guo Y, Chen S, Wu J, et al. Glutathione s-transferase mu 2 inhibits hepatic steatosis via ASK1 suppression. Commun Biol (2022) 5(1):326. doi: 10.1038/s42003-022-03251-w
30. Yan S, Wang Z, Wang Z, Duan Q, Wang X, Li J, et al. Glutathione s-transferase M1 null genotype related to poor prognosis of colorectal cancer. Tumour Biol (2016) 37(8):10229–34. doi: 10.1007/s13277-015-4676-4
31. Almeida M, Soares M, Ramalhinho AC, Moutinho JF, Breitenfeld L. Prognosis of hormone-dependent breast cancer seems to be influenced by KEAP1, NRF2 and GSTM1 genetic polymorphisms. Mol Biol Rep (2019) 46(3):3213–24. doi: 10.1007/s11033-019-04778-8
32. Gonlugur U, Pinarbasi H, Gonlugur TE, Silig Y. The association between polymorphisms in glutathione s-transferase (GSTM1 and GSTT1) and lung cancer outcome. Cancer Invest (2006) 24(5):497–501. doi: 10.1080/07357900600814813
33. Governa V, Trella E, Mele V, Tornillo L, Amicarella F, Cremonesi E, et al. The interplay between neutrophils and CD8(+) T cells improves survival in human colorectal cancer. Clin Cancer Res (2017) 23(14):3847–58. doi: 10.1158/1078-0432.CCR-16-2047
34. Kong JC, Guerra GR, Pham T, Tornillo L, Amicarella F, Cremonesi E, et al. Prognostic impact of tumor-infiltrating lymphocytes in primary and metastatic colorectal cancer: A systematic review and meta-analysis. Dis Colon Rectum (2019) 62(4):498–508. doi: 10.1097/DCR.0000000000001332
35. Deb S, Xu H, Tuynman J, George J, Yan Y, Li J, et al. RAD21 cohesin overexpression is a prognostic and predictive marker exacerbating poor prognosis in KRAS mutant colorectal carcinomas. Br J Cancer (2014) 110(6):1606–13. doi: 10.1038/bjc.2014.31
36. Zhu T, Gao Z, Yuan K, Wang Y. High expression of RAD21 predicts poor survival in patients with operated non-small-cell lung cancer. Tumori (2020) 106(3):223–8. doi: 10.1177/0300891620910805
Keywords: GSTM2, immune microenvironment, prognostic biomarker, colon cancer, RAD21
Citation: Zhang W, Shi Y, Niu S, Li L, Lin L, Gao X, Cai W, Chen Y, Zhong Y, Tang D, Tang M and Dai Y (2022) Integrated computer analysis and a self-built Chinese cohort study identified GSTM2 as one survival-relevant gene in human colon cancer potentially regulating immune microenvironment. Front. Oncol. 12:881906. doi: 10.3389/fonc.2022.881906
Received: 15 April 2022; Accepted: 13 September 2022;
Published: 03 October 2022.
Edited by:
Yu-Chan Chang, National Yang Ming Chiao Tung University, TaiwanCopyright © 2022 Zhang, Shi, Niu, Li, Lin, Gao, Cai, Chen, Zhong, Tang, Tang and Dai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Dai, daiyong22@aliyun.com; dai.yong@szhospital.com; Donge Tang, donge66@126.com; Min Tang, catom@126.com
†These authors have contributed equally to this work