Can Fang
Bing-bing Zou- Department of General Surgery, First Affiliated Hospital of Anhui Medical University, Hefei, China
Background: G Protein Subunit Gamma 7 (GNG7) is an important gene that regulates cell proliferation and induces apoptosis. However, the correlation between GNG7 expression and immune infiltration as well as patient prognosis of colorectal cancer (CRC) remains unclear.
Methods: The GNG7 expression differences between tumor tissues and normal tissues were explored via the Oncomine database, Tumor Immune Estimation Resource (TIMER) site and UALCAN database. Then, the influence of GNG7 on clinical prognosis were evaluated, using the PrognoScan database. In addition, the relationship between GNG7 and tumor-related immune infiltration as well as gene marker sets of immune infiltration was investigated via TIMER, TISIDB and GEPIA.
Results: We found that GNG7 expression was down-regulated in multiple malignant tumors including colorectal cancer (CRC) and the GNG7 expression was associated with tumor stage, histology subtype, lymph node metastasis and poor prognosis in colorectal cancer (CRC). In addition, the expression of GNG7 was significantly associated with infiltration level of multiple immune cells, immunomodulatory factors as well as part of the immune cell markers.
Conclusion: GNG7 displays validated prognostic value in CRC and was associated with its immune cell infiltration and immunoregulation. These results suggest that GNG7 is a potential prognostic marker and is associated with tumor immune infiltration, thus providing a new perspective for the immunotherapy of CRC.
Introduction
In 2020, there were approximately more than 1.9 million new cases of colorectal cancer (CRC) and 935,000 deaths from colorectal cancer globally. Overall, colorectal cancer accounting for the third most common cancer in the world, but ranks second in terms of mortality (Sung et al., 2021). At present, surgery is still the cornerstone of the treatment of CRC, but its shortcomings are also obvious. For example, the choice of surgical methods and the effectiveness of surgical treatment are largely affected by the depth of tumor invasion and anatomical location (Brenner et al., 2014). Although significant progress has been made in the screening, diagnosis, surgery and chemotherapy for patients with colorectal cancer in recent decades, the patients’ clinical benefits acquired from traditional treatment approaches are limited and the patients’ prognosis still needs to be improved. It is reported that about 25% of patients with CRC have metastases at the time of diagnosis, and approximately 50% of them will eventually metastasize in their lifetime (Vatandoust et al., 2015).
Multiple immune cells are vital components of the tumor microenvironment (TME) and deeply influence tumor occurrence and development (Hinshaw and Shevde 2019). Nowadays, immunotherapy has achieved satisfactory effects in dealing with some malignant tumors such as small cell lung cancer, renal cancer and breast cancer (Motzer et al., 2015; Gandhi et al., 2018; Schmid et al., 2018). However, there are still many unknown areas worth exploring in the field of immunotherapy for colorectal cancer. For example, the application of ICIS (immune checkpoint inhibitors) in CRC has epoch-making significance for the immunotherapy of CRC. Nevertheless, less than 15% of colorectal cancer patients exhibit this predictive biomarker, so the vast majority of the remaining CRC patients currently have no opportunity to extend survival with immunotherapy (Lichtenstern et al., 2020). Therefore, to determine immunotherapy targets by studying the influencing factors of tumor-related immune infiltration will enrich the immunotherapy approaches of colorectal cancer and help improve the clinical outcome of cancer patients.
GNG7 belongs to the large guanine nucleotide binding protein γ family and is located on chromosome 19 (Shibata et al., 1998). Related genes of the G protein family are involved in a variety of transmembrane signaling transduction pathways. In addition, studies have shown that GNG7 gene may also be associated with transmembrane signaling pathways and is involved in cell contact-induced growth arrest, thereby inhibiting uncontrolled cell proliferation in multicellular organisms (Hartmann et al., 2012). GNG7 is ubiquitously expressed in multiple tissues but is down-regulated in various cancers (Liu et al., 2016). Several studies have shown that its expression is decreased in pancreatic cancer (Liu et al., 2019), esophagus cancer, gastrointestinal tract cancer (Shibata et al., 1999), and clear cell renal cell carcinoma (ccRCC) (Xu et al., 2019). Combined with the research conclusions listed above, it can therefore be regarded as a potential suppressor of malignant tumor.
The expression profiles and prognostic value of GNG7 in CRC are still unclear. In addition, the association between GNG7 and immune infiltration in CRC still remains largely unexplored. Elucidating these issues will pave the way for the development of new treatments for CRC.
In this study, we analyzed the expression of GNG7 in various cancers through Oncomine, UALCAN, and TIMER databases. Furthermore, the prognostic significance of GNG7 in different tumors was investigated by the PrognoScan web servers. Moreover, we investigated the correlation between GNG7 expression and immune infiltration as well as immunoregulation level in colorectal cancer via TIMER, TISIDB, and GEPIA databases. Our findings revealed the prognostic value of GNG7 in CRC, and provide new insights into the correlation and potential mechanism between GNG7 expression and CRC immune infiltration.
Materials and Methods
Tissue Samples
We collected 20 samples of colorectal cancer and para-cancerous tissues as well as related clinical information from the First Affiliated Hospital of Anhui Medical University (Supplementary Table S2). None of the subjects received radiotherapy or chemotherapy before surgery. All patients’ data were obtained with informed consent, and this study was approved by the Clinical Research Ethics Committees. (Tissue specimen collection time is from 10 June 2021 to 20 July 2021).
Quantitative Real-Time PCR
TRIzol reagent (Invitrogen, Carlsbad, CA, United States) was used for total RNA extraction from tissues. The assay was conducted in triplicate and performed using 2× TB Green Premix Ex Taq (Takara) in a Bio-Rad detection system. All reactions were repeated. Relative gene expression was calculated by the 2–ΔΔCt method. GAPDH served as the internal reference. The primer sequences were as follows (“F” represents “forward”; “R”represents“reverse”). GAPDH: 5′-AGAAGGCTGGGGCTCATTTG-3′(F), 5′AGGGGCCATCCACAGTCTTC-3′(R); GNG7:5′-CCCGAGCGCAGGGAGT-3′(F), 5′-CGCTCAATCCCGGCTTCTAT-3′(R).
Oncomine Database Analysis
The oncomine database (http://www.oncomine.org/) was designed for the mining and analysis of tumor-related data (Rhodes et al., 2007). We employed it to analyze the expression differences of GNG7 in tumor tissues and adjacent normal tissues on it. The threshold was defined as a p-value < 0.01, a fold change>2, and a gene ranking of 10%.
PrognoScan Database Analysis
PrognoScan (http://www.prognoscan.org/) database can be used to analyze the relationship between gene expression and different prognostic outcomes (including OS, DFS, etc.) of patients (Mizuno et al., 2009). We analyzed the correlation between GNG7 expression and survival outcomes in different types of cancers and generated corresponding survival curves.
UALCAN Database Analysis
UALCAN (http://ualcan.path.uab.edu/index.html) is a practical online analysis database with multiple function modules, which is widely used in TCGA data analysis (Chandrashekar et al., 2017). We used UALCAN to analyze the mRNA and protein expression levels of GNG7 in tumor and adjacent normal tissues, as well as its relationship with clinicopathological features.
Tumor Immune Estimation Resource Database Analysis
Tumor Immune Estimation Resource (https://cistrome.shinyapps.io/timer/) (TIMER) platform is a comprehensive resource database which can be used to analyze differential expression of genes and the relationship between specific gene expression and the abundance of immune cell infiltration in various cancers, and the abundance of tumor-infiltrating immune cell (TIIC) is calculated use specific algorithms (Li et al., 2017). We used TIMER to explore the GNG7 expression in different types of cancers and the correlation between GNG7 expression and the infiltration abundance of six kinds of immune cells (TIICs, including B cells, CD4 + T cells, CD8 + T cells, neutrophils, macrophages, and dendritic cells (DCs)). In addition, the relationship between GNG7 expression and immune cell marker genes was further analyzed, and the results were presented in the form of a scatter plot with the Spearman correlation coefficient and p-value. Gene expression levels were shown as log2 TPM.
Gene Expression Profiling Interactive Analysis Database Analysis
The Gene Expression Profiling Interactive Analysis (http://gepia.cancer-pku.cn/index.html) (GEPIA) database integrates wide bioinformatics data, which can be used for a variety of bioinformatics research including genetic difference analysis, prognosis analysis and gene correlation analysis (Tang et al., 2017). The correlation between GNG7 expression and specific immune cell marker genes identified in TIMER database previously was further proved in GEPIA.
TISIDB Database Analysis
TISIDB database (https://cis.hku.hk/TISIDB) is a web portal for conducting tumor and immune system interactions analysis that integrates multiple heterogeneous data types (Ru et al., 2019). We used this database to help investigate the correlation between GNG7 expression and tumor-infiltrating lymphocytes as well as immunomodulators.
Schematic Diagram and Flow Chart
Well-designed diagrams help readers better understand the researcher’s research content (Behzadi and Gajdács 2021). For this reason, we present the schematic diagram (Figure 1) and technical flow chart (Figure 2) of the entire research process in the article for easy reference.
FIGURE 1. Illustration of the main principle of this study.
FIGURE 2. The basic analysis process of this study.
Statistical Analysis
Survival curves were generated from the PrognoScan and the results displayed with HR and P or Cox p-value from a log-rank test. The results generated from Oncomine were displayed with p-value, fold changes and ranks. The correlation of gene expression with the lymphocytes, immunomodulators and gene markers of tumor-infiltrating immune cells was evaluated by Spearman’s correlation coefficient and statistical significance. p-value <0.05 is considered statistically significant.
Results
The mRNA Expression Levels of GNG7 in Different Types of Human Cancers
The differences of GNG7 expression between tumor and normal tissues were analyzed in Oncomine database. It was observed that the GNG7 expression was up-regulated in kidney cancer, leukemia and lung cancer compared with the normal tissues, while it was down-regulated significantly in some datasets of bladder cancer, brain and CNS cancer, breast cancer, colorectal cancer, gastric cancer, head and neck cancer, leukemia, lymphoma, ovarian cancer and pancreatic cancer (Figure 3A, Supplementary Table S1).
FIGURE 3. GNG7 mRNA expression in a variety of human cancers. (A) GNG7 expression was downregulated or upregulated in diverse cancers in the Oncomine database. (B) Human GNG7 expression in multiple types of cancer was determined using the TIMER database (*p < 0.05, **p < 0.01, ***p < 0.001).
To further verify the outcome from Oncomine database, we examined the GNG7 expression levels in multiple malignancies from TCGA database in TIMER. The differential expression between the tumor and adjacent normal tissues of GNG7 was shown in Figure 3B, and as we can see, GNG7 was down-regulated in various cancers, especially in colon cancer and rectal cancer, which was consistent with the results in Oncomine. In addition, we found that GNG7 expression was not only down-regulated in COAD and READ, but it was also associated with tumor histology, stage, lymph node metastasis according to the result in UALCAN (Figure 4).
FIGURE 4. The relationship between GNG7 expression and individual clinical parameters of COAD and READ in UALCAN database. (A) Differential expression of GNG7 in COAD tissues. Association of GNG7 mRNA expression with different COAD patients’ individual cancer histologic subtypes (B), tumor stages (C), nodal metastasis (D). (E) Differential expression of GNG7 in READ tissues. Association of GNG7 mRNA expression with different READ patients’ individual cancer histologic subtypes (F), tumor stages (G), nodal metastasis (H).
We used qRT-PCR to identify the expression of GNG7 in CRC tissues and para-cancerous tissues. The results are displayed in the form of a histogram. It is confirmed that the GNG7 mRNA level was decreased in CRC tissues than para-cancerous tissues (Figure 5, In addition, experimental related pictures are shown in Supplementary File S1). These results suggests that GNG7 was obviously lowly expressed in various cancers, especially in CRC.
FIGURE 5. The results of RT-qPCR experiments are consistent with the analysis results of Oncomine and TIMER databases. The expression level of GNG7 mRNA in colorectal cancer tumor tissues is significantly down-regulated compared with normal para-cancerous tissues. ***p < 0.0001.
Prognostic Potential of GNG7 Across Different Cancers
In order to investigate the prognostic value of GNG7 expression in different types of cancers, we further analyzed the relationship between GNG7 expression and prognosis of cancer patients via PrognoScan database. The analysis uncovered that GNG7 expression was significantly related to the prognosis of eight types of tumors including lung cancer, breast cancer and colorectal cancer (Figure 6). For example, it can be observed that in the GSE6532-GPL570 cohort which included 87 breast cancer samples, the discrepancies of DMFS (Distant Metastasis Free Survival) and RFS (Relapse Free Survival) between two groups were significant. In the GSE31210 cohort which included 204 lung cancer samples, high GNG7 expression was associated with better OS (Overall Survival, p = 0.0005, HR = 0.29,95%CI 0.14–0.58) and RFS (p = 0.0002, HR = 0.25,95%CI 0.13–0.49). Notably, In the GSE17536 and GSE17537 cohorts, which included 177 and 55 colorectal cancer samples respectively, high GNG7 expression was associated with better OS (p = 0.034, HR = 0.08,95%CI 0.01–0.83) and DSS (Disease Specific Survival, p = 0.031, HR = 0.28,95%CI 0.09–0.89). Hence, these findings indicated that GNG7 can be regarded as a potential indicator for predicting the prognosis of certain types of cancer including CRC.
FIGURE 6. Association of GNG7 mRNA expression with prognosis in different human cancers in the PrognoScan. (A–B) Survival curves of DMFS and RFS in a breast cancer cohort (GSE6532-GPL570, n = 87). (C–D) Survival curves of OS and RFS in a lung cancer cohort (GSE31210, n = 204). (E–F) Survival curves of DSS and OS in two colorectal cancer cohorts [GSE17536 (n = 177) and GSE17537 (n = 55)].
GNG7 Expression is Correlated With Immune infiltration Level in Colorectal Cancer
The infiltration of tumor-related immune cells was regarded as an independent factor affecting the prognosis of patients with colorectal cancer (Ohtani 2007). Therefore, we explored the relationship between GNG7 expression and the abundance of immune infiltration in different types of tumors via TIMER database. Tumor purity is an important reference factor in analyzing the degree of tumor immune infiltration (Yoshihara et al., 2013). The results showed that the expression of GNG7 was significantly correlated with tumor purity in 21 types of cancer, including COAD and READ. In addition, the expression of GNG7 had significant correlations with the infiltration level of B cells in 18 tumors, CD8+ T cells in 17 tumors, CD4+ T cells in 27 tumors, macrophages in 22 tumors, neutrophils in 14 tumors, and dendritic cells in 21 tumors, respectively. Notably, we found that the GNG7 expression level was correlated with poor prognosis and high immune infiltration in COAD and READ. GNG7 expression level had significant positive correlations with infiltrating levels of B cells (r = 0.282, p = 8.32e-09), CD4+ T cells (r = 0.504, p = 2.47e-27), macrophages (r = 0.320, p = 4.35e-11), neutrophils (r = 0.215, p = 1.36e-05) and DCs (r = 0.337, p = 4.13e-12) in COAD (Figure 7A). Moreover, GNG7 expression levels were also positively associated with infiltration levels of B cells (r = 0.305, p = 2.56e-04), CD4+ T cells (r = 0.435, p = 8.67e-08), macrophages (r = 0.307, p = 2.37e-04) and DCs (r = 0.380, p = 3.96e-06) in READ (Figure 7B). In contrast, GNG7 expression had no significant correlation with tumor purity and infiltrating levels of B cells, CD8+ T cells, macrophages, neutrophils, or dendritic cells in LIHC (Figure 7C).
FIGURE 7. Correlation of GNG7 expression with immune infiltration level in COAD (colon adenocarcinoma), READ (rectum adenocarcinoma), and LIHC (liver hepatocellular carcinoma). (A) GNG7 expression is significantly negatively related to tumor purity and has significant positive correlations with infiltrating levels of B cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells in COAD, other than CD8+ T cells. (B) GNG7 expression is significantly negatively related to tumor purity and has significant positive correlations with infiltrating levels of B cells, CD4+ T cells, macrophages, and dendritic cells in READ but no significant correlation with infiltrating level of CD8+ T cells and neutrophils. (C) GNG7 expression has no significant correlations with tumor purity and infiltrating levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells in LIHC.
The Function of GNG7 in Regulating Immune Molecules
Then, in order to study the regulation mechanism of GNG7 on immune molecules in CRC, we further conducted an integrated analysis on the correlation between GNG7 expression and 28 tumor immune infiltrating cell subtypes as well as three subtypes of immunomodulators in TISIDB database. The result showed that GNG7 was associated with twenty-four immune cell subtypes in COAD and READ (Table 1). Notably, Activated B cell (rho = 0.63, p < 2.2e-16), Immature B cell (rho = 0.536, p < 2.2e-16) and Eosinophil (rho = 0.445, p < 2.2e-16) occupied the top three TILs in terms of the strength of the correlation with GNG7 expression in COAD (Figures 8A–C), while the top three TILs with the strongest correlation with GNG7 expression in READ were Activated B cell (rho = 0.544, p < 2.2e-16), Immature B cell (rho = 0.482, p = 5.23e-11) and Macrophage (rho = 0.473, p < 1.5e-10) (Figure 8D–F). Immunoinhibitors, immunostimulators and major histocompatibility complex (MHC) molecules comprised all immunomodulators (Gou et al., 2019). The correlation heat maps of GNG7 expression with immunoinhibitors, immunostimulators and major histocompatibility complex (MHC) molecules in different cancers were displayed respectively in (Figures 8B–D). In addition, the correlation scatter plots of the three immunomodulators which had the strongest correlation with GNG7 expression in COAD and READ were displayed near the related heat maps. The specific coefficients and p value were also shown in Supplementary Table S3. These findings provided compelling evidence for demonstrating that GNG7 plays a specific role in immune infiltration in colorectal cancer.
TABLE 1. The correlation between GNG7 expression and tumor lymphocyte infiltration in COAD, READ and LIHC(TISIDB).
FIGURE 8. Spearman’s correlation of GNG7 with lymphocytes and immunomodulators in TISIDB database. (A) Relationship between the abundance of TILs and GNG7 expression in different human cancers. Top 3 TILs displaying the greatest correlation with GNG7 expression in COAD (A–C) and READ (D–F). (B) Relationship between the abundance of immunoinhibitors and GNG7 expression in different human cancers. Top 3 immunoinhibitors displaying the greatest correlation with GNG7 expression in COAD (A–C) and READ (D–F). (C) Relationship between the abundance of immunostimulators and GNG7 expression in different human cancers. Top 3 immunostimulators displaying the greatest correlation with GNG7 expression in COAD (A–C) and READ (D–F). (D) Relationship between the abundance of MHC molecules and GNG7 expression in different human cancers. Top 3 MHC molecules displaying the greatest correlation with GNG7 expression in COAD (A–C) and READ (D–F).
Correlation Analysis Between GNG7 Expression and Immune Markers
To further investigate the relationships between GNG7 and immune infiltrating cells, we used TIMER and GEPIA database to analyze the correlations between GNG7 expression and immune cell markers of diverse immune cells in COAD and READ. These immune cell markers were used to characterize immune cells, including B cells, CD8+ T cells, T cells (general), M1/M2 macrophages, tumor-associated macrophages (TAMs), monocytes, NK cells, neutrophils and DCs. We also analyzed different functional T cells, such as Th1 cells, Th2 cells, Tfh cells, Th17 cells, and Tregs, as well as exhausted T cells. The results revealed that GNG7 was significantly correlated with most immune marker sets of multiple immune cells and different T cells in COAD and READ (Table 2). In contrast, only 13 immune cell gene markers in LIHC were significantly correlated with the GNG7 expression.
TABLE 2. Correlation analysis between GNG7 and relate genes and markers of immune infiltration cells in TIMER.
Interestingly, we found that compared with COAD, the expression of GNG7 in READ is associated with more immune infiltrating cell markers and the correlation coefficient is higher. As shown in Table 2, among these TILs, B cells, TAM, M2 macrophages, dendritic cells and Treg cells displayed stronger correlation with GNG7 expression than others. Therefore, we further explored the correlation between GNG7 expression in COAD, READ and these gene markers by using GEPIA database. The analysis results showed that the degree of correlation is basically the same as the results we obtained in TIMER Table 3.
TABLE 3. Correlation analysis between GNG7 and relate genes and markers of immune infiltration cells in GEPIA.
Consequently, we have reason to believe that GNG7 is closely correlated to immune cell infiltration in CRC, which guides us to speculate that GNG7 plays an important role in the tumor microenvironment and immune response of CRC.
Discussion
Colorectal cancer is a kind of cancer with a high level of morbidity and mortality worldwide and its clinical diagnosis and treatment progress has always been focused on (Arnold et al., 2017). Moreover, due to the lack of effective early detection methods, combined with its low penetration, many patients have already reached advanced stage at the time they received clinical diagnosis (Dekker et al., 2019). Therefore, the discovery of new, relatively convenient, and non-invasive detection markers will help improve the early screening rate of CRC, thereby increasing the chance of patients being treated and improving the prognosis of patients.
In recent years, with the progress of tumor molecular biology research, more and more new and effective markers have been discovered, for example, Ritwik Patra’s research found that COL11A1 gene mutation profoundly affects the occurrence and development of colorectal cancer, and confirmed the qualification of COL11A1 gene as a potential biomarker for colorectal cancer prognosis (Patra et al., 2021). In addition, it is reported that the use of combined detection of SDC2 and TFPI2 methylation as a new biomarker for a non-invasive tool of CRC will help improve the current dilemma of low early screening rates for CRC (Zhang et al., 2021). Meanwhile, as a cell cycle-related protein, SEPT9 has been confirmed to be involved in tumorigenesis and development in a variety of malignant tumors, and the methylation detection of SERT9 gene in peripheral blood has high sensitivity and specificity for the diagnosis of CRC, and can be used in early CRC screening, it may be a useful biomarker for the prognosis of CRC (Mostowy and Cossart 2012; Sun et al., 2019; Yang et al., 2019).
Immune cells are an important part of tumor stroma. Accumulating evidences suggests that innate immune cells (including macrophages, neutrophils, dendritic cells, NK cells, innate lymphocytes and myeloid-derived suppressor cells) as well as adaptive immune cells (T cells and B cells) participate in the formation of the tumor microenvironment (TME). The crosstalk between cancer cells and proximal immune cells will eventually create an environment that fosters tumor growth and metastasis (Gajewski et al., 2013; Hinshaw and Shevde 2019). Moreover, the high level of tumor immune infiltration may be related to poor prognosis in CRC (Guo et al., 2020).
At present, with the in-depth understanding of the mechanism of tumor immune infiltration, the tumor immunotherapy has also been rapidly developed. In particular, tumor immunotherapy has yielded significant effects in a variety of cancer types (Rosenberg et al., 2019; Brewer et al., 2020; Farid and Liu 2020).
However, the progress of immunotherapy for colorectal cancer is still relatively limited. The most important reason is that most CRCs are microsatellite stable, but it is disappointing that currently immunotherapy only produces a good response in a small subset of CRCs that exhibit microsatellite instability (MSI), and the difference in the sensitivity of the two types of CRC to immunotherapy is mainly related to the immune microenvironment at the primary site of the tumor (Mlecnik et al., 2016). Therefore, discovering new immune-related markers and therapeutic targets become the main content of research on CRC immunotherapy.
This study aims to predict the role of GNG7 in CRC by studying the relationship between the expression level of GNG7 and the infiltration level of different immune cells in colorectal cancer and the prognosis of patients with CRC through multi-channel bioinformatics data mining and analysis.
GNG7 is located on chromosome 19 and is a member of the guanine nucleotide-binding proteins (G proteins) family which is involved as a modulator or transducer in various transmembrane signaling systems (Shibata et al., 1998; Hartmann et al., 2012). Previous studies have reported that the loss of GNG7 is related to the tumor volume of esophageal cancer and affects tumor invasiveness, which may attributed to loss of heterozygosity (LOH) (Ohta et al., 2008). Zheng H et al.'s study on LUAD found that GNG7 can exert a tumor suppressor effect in LUAD cells by inhibiting E2F1, moreover, low GNG7 expression is significantly associated with a poorer LUAD prognosis and a higher tumor grade (Zheng et al., 2021). However, the role of GNG7 in colorectal cancer remains to be elucidated. In order to acquire more detailed information about the potential role and biological mechanism of GNG7 expression in CRC, we conducted a stepwise bioinformatics analysis. First, we performed analyses in Oncomine database and TIMER database and found that compared to normal tissues, GNG7 expression in CRC was significantly down-regulated. In addition, the findings from UALCAN indicated that GNG7 expression is associated with tumor histology, stage, lymph node metastasis in colorectal cancer which means that GNG7 may play an important role in the occurrence and progression of colorectal cancer. In addition, we further carried out the qRT-PCR experiment of clinical colorectal tissue samples, and the results are consistent with the results in the database as well. Meanwhile, it is worth noting that analysis of the Oncomine and TIMER databases revealed that the expression levels of GNG7 in various cancers were different relative to normal tissues, however, GNG7 was down-regulated in some cancers and up-regulated in others. Unfortunately, The specific mechanism causing this result has not yet been fully elucidated. We speculate that this may be related to differences in cancer cell tissue types, or that the changes in signaling pathways vary in different cancers, which may also lead to increased or decline expression of downstream proteins. which needs to be determined by further basic experimental research in the future.
Cancer cells are distinguished from normal cells by two major characteristics: immortal proliferation and decreased apoptosis (Hanahan and Weinberg 2011). It is reported that GNG7 has the ability to induce apoptosis and arrest the cell cycle, which may be related to the reduction of actin polymers in cells. In addition, GNG7 induces autophagy by inhibiting the MTOR pathway, thereby limiting the number of cells, reducing the proliferation ability to a certain extent, and playing the role of tumor suppressor genes (Liu et al., 2016; Lai and Ding 2019).
Our survival analysis results in the database showed that GNG7 expression is correlated to the prognosis of different cancers, including CRC. Low GNG7 level means poor colorectal cancer DSS and OS. Therefore, we speculated that GNG7 exhibits significance as a prognostic marker in CRC.
Another important part of this study is to explore the correlation between the expression of GNG7 and the level of infiltration of different immune cells as well as immunoregulators in cancer, especially in colorectal cancer. Our analysis showed that the relationship between GNG7 expression and the level of CD4+ T cell infiltration in CRC is moderate to strong. Simultaneously, there is a positive correlation between GNG7 expression and the infiltration level of B cells, macrophages, neutrophils and dendritic cells. In addition, the GNG7 expression has a significant correlation with numerous immunoregulators as well as some immune cell marker genes in COAD and READ, which implies that GNG7 may be involved in the tumor immune process of COAD and READ, and may play a significant role in tumor-related immune activation or suppression. Hence, we hypothesized that GNG7 is a potential marker for predicting tumor immune microenvironment homeostasis.
Macrophages belong to a differentiated subtype of immune cells, which can be divided into M1 macrophages and M2 macrophages (Jayasingam et al., 2020). They play an important role in the development, homeostasis and immune response of organism (Wynn et al., 2013). Increased macrophage infiltration in colorectal cancer tissues was significantly associated with increased chemoresistance and poor prognosis (Yin et al., 2017).
It has now been fully confirmed that a large number of Treg cell infiltrations in tumor tissues of various cancers are ubiquitous, and their abundant presence often indicates a poor clinical prognosis. Among the various mechanisms of immunological self-tolerance, endogenous Treg cell-mediated immunosuppression is essential (Tanaka and Sakaguchi 2017). Dendritic cell is a quintessential APC of the immune system, responsible for combining innate immunity and adaptive immunity, including activation of anti-tumor T cells (Broz et al., 2014; Gardner et al., 2020). And dendritic cell immunotherapy has received good feedback in animal experiments, so it is regarded as a new cancer therapy which is worth looking forward to (Cranmer et al., 2004). Additionally, DC can assist tumor metastasis by increasing Treg cells and reducing the toxicity of CD8+ T cells (Sawant et al., 2012). These underlying mechanisms may be the basis for GNG7 to recruit and regulate immune cells. This provides a perspective for us to understand how the low expression of GNG7 affects the prognosis of CRC by regulating immune infiltration.
Therefore, we can conclude that GNG7 may be a potential biomarker and is related to the immune infiltration of colorectal cancer, which means that GNG7 can be used as a potential target for regulating anti-tumor immunotherapy. This study provided insights for further research on immunotherapy of CRC. Nevertheless, this study existed limitations, the data of this study is mainly based on data obtained in open access databases, lacking multi-center and large sample data support. In addition, the relevant mechanism of GNG7 in recruiting and regulating immune cells as well as immunoregulators needs to be further explored at the cellular and molecular level and in animal experiments in the future.
Conclusion
Our study confirmed that GNG7 expression is down-regulated in colorectal cancer, and GNG7 expression is related to the prognosis of colorectal cancer. moreover, the correlation between GNG7 gene and multiple immune cell infiltration levels, immunomodulatory molecules, and immune cell related genes in colorectal cancer revealed that GNG7 has a potential role in regulating immune cell function in colorectal cancer, which may provide a new insight to individualized treatment for CRC. In addition, further research should aim to elucidate the molecular mechanisms that may help determine the immune target of CRC.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
Study conception and design: B-bZ and CF. Data collection: CF, RZ and CQ. Statistical analysis: CF, RZ and CQ. Manuscript draft: CF. Manuscript review and editing: B-bZ.
Funding
This study was supported from Education Department of Anhui province fund projects which item number is: NO.KJ2019A0249.
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.
Acknowledgments
The authors thank the participating patients for the source of clinical tissue samples.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.833013/full#supplementary-material.
References
Arnold, M., Sierra, M. S., Laversanne, M., Soerjomataram, I., Jemal, A., and Bray, F. (2017). Global Patterns and Trends in Colorectal Cancer Incidence and Mortality. Gut 66, 683–691. doi:10.1136/gutjnl-2015-310912
Behzadi, P., and Gajdács, M. (2021). Writing a strong Scientific Paper in Medicine and the Biomedical Sciences: a Checklist and Recommendations for Early Career Researchers. Biol. Futura 72 (4), 395–407. doi:10.1007/s42977-021-00095-z
Brenner, H., Kloor, M., and Pox, C. P. (2014). Colorectal Cancer. Lancet 383 (9927), 1490–1502. doi:10.1016/S0140-6736(13)61649-9
Brewer, M., Angioli, R., Scambia, G., Lorusso, D., Terranova, C., Panici, P. B., et al. (2020). Front-line Chemo-Immunotherapy with Carboplatin-Paclitaxel Using Oregovomab Indirect Immunization in Advanced Ovarian Cancer: A Randomized Phase II Study. Gynecol. Oncol. 156 (3), 523–529. doi:10.1016/j.ygyno.2019.12.024
Broz, M. L., Binnewies, M., Boldajipour, B., Nelson, A. E., Pollack, J. L., Erle, D. J., et al. (2014). Dissecting the Tumor Myeloid Compartment Reveals Rare Activating Antigen-Presenting Cells Critical for T Cell Immunity. Cancer Cell 26 (5), 638–652. doi:10.1016/j.ccell.2014.09.007
Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B. V. S. K., et al. (2017). UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 19 (8), 649–658. doi:10.1016/j.neo.2017.05.002
Cranmer, L. D., Trevor, K. T., and Hersh, E. M. (2004). Clinical Applications of Dendritic Cell Vaccination in the Treatment of Cancer. Cancer Immunol. Immunother. 53 (4), 275–306. doi:10.1007/s00262-003-0432-5
Dekker, E., Tanis, P. J., Vleugels, J. L. A., Kasi, P. M., and Wallace, M. B. (2019). Colorectal Cancer. Lancet 394, 1467–1480. doi:10.1016/s0140-6736(19)32319-0
Farid, S., and Liu, S. V. (2020). Chemo-immunotherapy as First-Line Treatment for Small-Cell Lung Cancer. Ther. Adv. Med. Oncol. 12, 1758835920980365. doi:10.1177/1758835920980365
Gajewski, T. F., Schreiber, H., and Fu, Y.-X. (2013). Innate and Adaptive Immune Cells in the Tumor Microenvironment. Nat. Immunol. 14 (10), 1014–1022. doi:10.1038/ni.2703
Gandhi, L., Rodríguez-Abreu, D., Gadgeel, S., Esteban, E., Felip, E., De Angelis, F., et al. (2018). Pembrolizumab Plus Chemotherapy in Metastatic Non-Small-Cell Lung Cancer. N. Engl. J. Med. 378 (22), 2078–2092. doi:10.1056/NEJMoa1801005
Gardner, A., de Mingo Pulido, Á., and Ruffell, B. (2020). Dendritic Cells and Their Role in Immunotherapy. Front. Immunol. 11, 924. doi:10.3389/fimmu.2020.00924
Gou, R., Zhu, L., Zheng, M., Guo, Q., Hu, Y., Li, X., et al. (2019). Annexin A8 Can Serve as Potential Prognostic Biomarker and Therapeutic Target for Ovarian Cancer: Based on the Comprehensive Analysis of Annexins. J. Transl. Med. 17 (1), 275. doi:10.1186/s12967-019-2023-z
Guo, L., Wang, C., Qiu, X., Pu, X., and Chang, P. (2020). Colorectal Cancer Immune Infiltrates: Significance in Patient Prognosis and Immunotherapeutic Efficacy. Front. Immunol. 11, 1052. doi:10.3389/fimmu.2020.01052
Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of Cancer: the Next Generation. Cell 144 (5), 646–674. doi:10.1016/j.cell.2011.02.013
Hartmann, S., Szaumkessel, M., Salaverria, I., Simon, R., Sauter, G., Kiwerska, K., et al. (2012). Loss of Protein Expression and Recurrent DNA Hypermethylation of the GNG7 Gene in Squamous Cell Carcinoma of the Head and Neck. J. Appl. Genet. 53 (2), 167–174. doi:10.1007/s13353-011-0079-4
Hinshaw, D. C., and Shevde, L. A. (2019). The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res. 79 (18), 4557–4566. doi:10.1158/0008-5472.CAN-18-3962
Jayasingam, S. D., Citartan, M., Thang, T. H., Mat Zin, A. A., Ang, K. C., and Ch'ng, E. S. (2020). Evaluating the Polarization of Tumor-Associated Macrophages into M1 and M2 Phenotypes in Human Cancer Tissue: Technicalities and Challenges in Routine Clinical Practice. Front. Oncol. 9, 1512. doi:10.3389/fonc.2019.01512
Lai, W. S., and Ding, Y. L. (2019). GNG7 Silencing Promotes the Proliferation and Differentiation of Placental Cytotrophoblasts in Preeclampsia Rats through Activation of the mTOR Signaling Pathway. Int. J. Mol. Med. 43 (5), 1939–1950. doi:10.3892/ijmm.2019.4129
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472
Lichtenstern, C. R., Ngu, R. K., Shalapour, S., and Karin, M. (2020). Immunotherapy, Inflammation and Colorectal Cancer. Cells 9 (3), 618. doi:10.3390/cells9030618
Liu, J., Ji, X., Li, Z., Yang, X., Wang, W., and Zhang, X. (2016). G Protein γ Subunit 7 Induces Autophagy and Inhibits Cell Division. Oncotarget 7 (17), 24832–24847. doi:10.18632/oncotarget.8559
Liu, L., Wang, S., Cen, C., Peng, S., Chen, Y., Li, X., et al. (2019). Identification of Differentially Expressed Genes in Pancreatic Ductal Adenocarcinoma and normal Pancreatic Tissues Based on Microarray Datasets. Mol. Med. Rep. 20 (2), 1901–1914. doi:10.3892/mmr.2019.10414
Mizuno, H., Kitada, K., Nakai, K., and Sarai, A. (2009). PrognoScan: a New Database for Meta-Analysis of the Prognostic Value of Genes. BMC Med. Genomics 2, 18. doi:10.1186/1755-8794-2-18
Mlecnik, B., Bindea, G., Angell, H. K., Maby, P., Angelova, M., Tougeron, D., et al. (2016). Integrative Analyses of Colorectal Cancer Show Immunoscore Is a Stronger Predictor of Patient Survival Than Microsatellite Instability. Immunity 44, 698–711. doi:10.1016/j.immuni.2016.02.025
Mostowy, S., and Cossart, P. (2012). Septins: the Fourth Component of the Cytoskeleton. Nat. Rev. Mol. Cel. Biol. 13 (3), 183–194. doi:10.1038/nrm3284
Motzer, R. J., Escudier, B., McDermott, D. F., George, S., Hammers, H. J., Srinivas, S., et al. (2015). Nivolumab versus Everolimus in Advanced Renal-Cell Carcinoma. N. Engl. J. Med. 373 (19), 1803–1813. doi:10.1056/NEJMoa1510665
Ohta, M., Mimori, K., Fukuyoshi, Y., Kita, Y., Motoyama, K., Yamashita, K., et al. (2008). Clinical Significance of the Reduced Expression of G Protein Gamma 7 (GNG7) in Oesophageal Cancer. Br. J. Cancer 98 (2), 410–417. doi:10.1038/sj.bjc.6604124
Ohtani, H. (2007). Focus on TILs: Prognostic Significance of Tumor Infiltrating Lymphocytes in Human Colorectal Cancer. Cancer Immun. 7, 4.
Patra, R., Das, N. C., and Mukherjee, S. (2021). Exploring the Differential Expression and Prognostic Significance of the COL11A1 Gene in Human Colorectal Carcinoma: An Integrated Bioinformatics Approach. Front. Genet. 12, 608313. doi:10.3389/fgene.2021.608313
Rhodes, D. R., Kalyana-Sundaram, S., Mahavisno, V., Varambally, R., Yu, J., Briggs, B. B., et al. (2007). Oncomine 3.0: Genes, Pathways, and Networks in a Collection of 18,000 Cancer Gene Expression Profiles. Neoplasia 9 (2), 166–180. doi:10.1593/neo.07112
Rosenberg, A. R., Tabacchi, M., Ngo, K. H., Wallendorf, M., Rosman, I. S., Cornelius, L. A., et al. (2019). Skin Cancer Precursor Immunotherapy for Squamous Cell Carcinoma Prevention. JCI Insight 4 (6), e125476. doi:10.1172/jci.insight.125476
Ru, B., Wong, C. N., Tong, Y., Zhong, J. Y., Zhong, S. S. W., Wu, W. C., et al. (2019). TISIDB: an Integrated Repository portal for Tumor-Immune System Interactions. Bioinformatics 35 (20), 4200–4202. doi:10.1093/bioinformatics/btz210
Sawant, A., Hensel, J. A., Chanda, D., Harris, B. A., Siegal, G. P., Maheshwari, A., et al. (2012). Depletion of Plasmacytoid Dendritic Cells Inhibits Tumor Growth and Prevents Bone Metastasis of Breast Cancer Cells. J. Immunol. 189 (9), 4258–4265. doi:10.4049/jimmunol.1101855
Schmid, P., Adams, S., Rugo, H. S., Schneeweiss, A., Barrios, C. H., Iwata, H., et al. (2018). Atezolizumab and Nab-Paclitaxel in Advanced Triple-Negative Breast Cancer. N. Engl. J. Med. 379 (22), 2108–2121. doi:10.1056/NEJMoa1809615
Shibata, K., Mori, M., Tanaka, S., Kitano, S., and Akiyoshi, T. (1998). Identification and Cloning of Human G-Protein γ 7, Down-Regulated in Pancreatic Cancer. Biochem. Biophys. Res. Commun. 246 (1), 205–209. doi:10.1006/bbrc.1998.8581
Shibata, K., Tanaka, S., Shiraishi, T., Kitano, S., and Mori, M. (1999). G-protein Gamma 7 Is Down-Regulated in Cancers and Associated with P 27kip1-Induced Growth Arrest. Cancer Res. 59 (5), 1096–1101.
Sun, J., Fei, F., Zhang, M., Li, Y., Zhang, X., Zhu, S., et al. (2019). The Role of mSEPT9 in Screening, Diagnosis, and Recurrence Monitoring of Colorectal Cancer. BMC Cancer 19 (1), 450. doi:10.1186/s12885-019-5663-8
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Tanaka, A., and Sakaguchi, S. (2017). Regulatory T Cells in Cancer Immunotherapy. Cell Res. 27 (1), 109–118. doi:10.1038/cr.2016.151
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45 (W1), W98–W102. doi:10.1093/nar/gkx247
Vatandoust, S., Price, T. J., and Karapetis, C. S. (2015). Colorectal Cancer: Metastases to a Single Organ. World J. Gastroenterol. 21 (41), 11767–11776. doi:10.3748/wjg.v21.i41.11767
Wynn, T. A., Chawla, A., and Pollard, J. W. (2013). Macrophage Biology in Development, Homeostasis and Disease. Nature 496 (7446), 445–455. doi:10.1038/nature12034
Xu, S., Zhang, H., Liu, T., Chen, Y., He, D., and Li, L. (2019). G Protein γ Subunit 7 Loss Contributes to Progression of clear Cell Renal Cell Carcinoma. J. Cel. Physiol. 234 (11), 20002–20012. doi:10.1002/jcp.28597
Yang, X., Xu, Z.-J., Chen, X., Zeng, S.-S., Qian, L., Wei, J., et al. (2019). Clinical Value of Preoperative Methylated Septin 9 in Chinese Colorectal Cancer Patients. World J. Gastroenterol. 25 (17), 2099–2109. doi:10.3748/wjg.v25.i17.2099
Yin, Y., Yao, S., Hu, Y., Feng, Y., Li, M., Bian, Z., et al. (2017). The Immune-Microenvironment Confers Chemoresistance of Colorectal Cancer through Macrophage-Derived IL6. Clin. Cancer Res. 23 (23), 7375–7387. doi:10.1158/1078-0432.CCR-17-1283
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun. 4, 2612. doi:10.1038/ncomms3612
Zhang, W., Yang, C., Wang, S., Xiang, Z., Dou, R., Lin, Z., et al. (2021). SDC2 and TFPI2 Methylation in Stool Samples as an Integrated Biomarker for Early Detection of Colorectal Cancer. Cancer Manag. Res. 13, 3601–3617. doi:10.2147/CMAR.S300861
Zheng, H., Tian, H., Yu, X., Ren, P., and Yang, Q. (2021). G Protein Gamma 7 Suppresses Progression of Lung Adenocarcinoma by Inhibiting E2F Transcription Factor 1. Int. J. Biol. Macromol. 182, 858–865. doi:10.1016/j.ijbiomac.2021.04.082
Glossary
ccRCC clear cell renal cell carcinoma
CNS central nervous system
COAD colon adenocarcinoma
CRC colorectal cancer
dcs dendritic cells
DMFS distant metastasis free survival
DSS disease specific survival
GNG7 G protein subunit gamma 7
HR hazard ratio
ICIS immune checkpoint inhibitors
LIHC liver hepatocellular carcinoma
MHC molecule major histocompatibility complex molecule
NK cells natural killer cells
NK T cells natural killer T cells
OS overall survival
READ rectum adenocarcinoma
RFS relapse-free survival
tams tumor-associated macrophages
TCGA the cancer genome atlas
tems Tie2-expressing monocytes
tfh T follicular helper cell
th 1 cells type I helper T cells
th 2 cells type II helper T cells
Th17 cells type 17 T helper cells
TIIC tumor infiltrating immune cell
tils tumor infiltrating lymphocytes
TIMER tumor immune estimation resource
TME tumor microenvironment
treg regulatory T cells
Keywords: GNG7(G Protein Subunit Gamma 7), prognosis, CRC, colorecta cancer, immune infiltration
Citation: Fang C, Zhong R, Qiu C and Zou B-b (2022) The Prognostic Value of GNG7 in Colorectal Cancer and Its Relationship With Immune Infiltration. Front. Genet. 13:833013. doi: 10.3389/fgene.2022.833013
Received: 10 December 2021; Accepted: 31 January 2022;
Published: 23 February 2022.
Edited by:
Suprabhat Mukherjee, Kazi Nazrul University, IndiaReviewed by:
Parth Sarthi Sen Gupta, Indian Institute of Science Education and Research Berhampur (IISER), IndiaPayam Behzadi, Islamic Azad University, Iran
Copyright © 2022 Fang, Zhong, Qiu and Zou. 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: Bing-bing Zou, amy0550@163.com