Abstract
Background: G Protein Subunit Gamma 7 (GNG7), an important regulator of cell proliferation and cell apoptosis, has been reported to be downregulated in a variety of tumors including lung adenocarcinoma (LUAD). However, the correlation between low expression of GNG7 and prognosis of LUAD as well as the immune infiltrates of LUAD remains unclear.
Methods: The samples were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). R software was performed for statistical analysis. GNG7 expression and its prognostic value in LUAD were assessed through statistically analyzing the data from different databases. A nomogram was constructed to predict the impact of GNG7 on prognosis. Gene set enrichment analysis (GSEA) and single-sample gene set enrichment analyses GSEA (ssGSEA) were employed to determine the potential signal pathways and evaluated the immune cell infiltration regulated by GNG7. The prognostic significance of GNG7 expression associated with immune cell infiltration was investigated using the Tumor Immune Estimation Resource 2.0 (TIMER2.0) and the Kaplan-Meier plotter database. The UALCAN, cBio Cancer Genomics Portal (cBioPortal) and MethSurv database were used to analyze the correlation between the methylation of GNG7 and its mRNA expression as well as prognostic significance.
Results: GNG7 was demonstrated to be down-regulated in LUAD and its low expression was associated with poor prognosis. A clinical reliable prognostic-predictive model was constructed. Pathway enrichment showed that GNG7 was highly related to the B cell receptor signaling pathway. Further analysis showed that GNG7 was positively associated with B cell infiltration and low levels of B cell infiltration tended to associate with worse prognosis in patients with low GNG7 expression. Moreover, methylation analysis suggested hypermethylation may contribute to the low expression of GNG7 in LUAD.
Conclusion: Decreased expression of GNG7 at least partly caused by hypermethylation of the GNG7 promoter is closely associated with poor prognosis and tumor immune cell infiltration (especially B cells) in LUAD. These results suggest that GNG7 may be a promising prognostic biomarker and a potential immunotherapeutic target for LUAD, which provides new insights into immunotherapy for LUAD.
Introduction
Lung cancer is the most common reason for global cancer-related mortality, of which lung adenocarcinoma (LUAD) is the most common histological subtype (Travis, 2011; Sung et al., 2021). In recent years, although molecular targeted therapies and immunotherapy have significantly improved the prognosis of a small proportion of LUAD patients, the above treatments are ineffective in many patients due to the high heterogeneity and complexity of LUAD (Molina et al., 2008; Saito et al., 2018; Wu and Shih, 2018). The prognosis for many patients, especially those with advanced LUAD, remains poor, with a 5-year survival rate of less than 18% (Singh et al., 2020). Therefore, an in-depth pathogenetic exploration and search for other effective diagnostic and therapeutic approaches as well prognostic markers are essential to improve the prognosis of patients with LUAD.
Accumulating evidence suggests that the immune cells within the tumor microenvironment (TME) play essential roles in tumorigenesis (Hinshaw and Shevde, 2019). Such tumor associated immune cells may exert pro-tumor or anti-tumor function in the initiation and development of tumors (Taube et al., 2018). Studies have shown that immune cell infiltration is an important factor influencing the efficacy of immunotherapy (Martinez and Moon, 2019; Petitprez et al., 2020; Bagchi et al., 2021). In addition, TME is also closely related to the prognosis of patients (Qi et al., 2020; Wu et al., 2020). Therefore, it is crucial to investigate the regulators of immune cell infiltration in the TME to improve the effectiveness of immunotherapy and improve patient prognosis.
G Protein Subunit Gamma 7 (GNG7), a subunit of heterotrimeric G protein, is strongly enriched in the striatum and plays a vital role in the A2A adenosine or D1 dopamine receptor-induced neuro-protective response (Schwindinger et al., 2012). Multiple studies have shown that GNG7 expression is decreased in many cancers, including pancreatic cancer, gastrointestinal tract cancer and renal carcinoma (Shibata et al., 1998; Shibata et al., 1999; Ohta et al., 2008). GNG7 overexpression was shown to inhibit cell growth and tumorigenicity of esophageal carcinoma cells (Hartmann et al., 2012). Also, GNG7 was confirmed as an essential autophagy-inducing agent and participated in inhibiting tumor progression through mTOR pathway (Xu et al., 2019). Recently, GNG7 was reported to be lowly expressed in LUAD and promoted the progression of LUAD through Hedgehog signaling (Liu et al., 2016). These findings indicated that GNG7 may be a potential tumor suppressor implicated in the carcinogenesis and tumor progression. However, the detailed roles and mechanisms of GNG7 especially the effects of GNG7 on immune infiltration in LUAD are largely unknown.
In this study, we aimed to evaluate the clinical significance of GNG7 in LUAD and the possible mechanisms underlying its function through comprehensive bioinformatics analysis. Our results showed that low expression of GNG7 positively correlated with the progression of LUAD, and GNG7 may be an important potential prognostic biomarker for LUAD. We also constructed a reliable clinical prediction model. In addition, we revealed for the first time the underlying mechanisms of GNG7 dysregulation and the correlation between GNG7 expression and immune infiltration in LUAD, which implies that GNG7 could be a potential target for clinical antitumor immunotherapy.
Materials and methods
Data sources and pretreatment
The RNA-seq data of 513 LUAD samples and 59 normal samples were downloaded from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/). The downloaded data format was level 3 HTSeq-fragments per kilobase per million (FPKM) and then was converted to transcripts per million (TPM) format for subsequent analysis. TCGA supplemented prognostic data were obtained from a Cell article (Liu et al., 2018). In addition, three sets of microarray data of LUAD tissues (accession numbers: GSE32665, GSE32863, and GSE43458) were downloaded from the GEO database. The data used in this study were obtained from both the TCGA database and the GEO database, which ensured that all written informed consent was obtained prior to data collection.
Key gene screening
The R package “DESeq2” was used to identify the differentially expressed genes (DEGs) between LUAD tissues and normal tissues (Love et al., 2014). Adjusted p-value <0.05 and | Log2fold change| >1 were set as cut-off criteria. The Survminer R package and the survivor R package were used to screen for genes with better prognostic value. The cut-off threshold was hazard ratio (HR) < 1, and the Cox p-value <0.005. The prognostic indicators included overall survival (OS), disease-specific survival (DSS) and progression-free survival (PFS). The Venn diagram was used to represent the intersection set of DEGs obtained from the four sets.
G protein subunit gamma 7 differential expression analysis in lung adenocarcinoma tissues
Pre-processed TCGA-LUAD data were used for differential expression analysis of GNG7 in LUAD tumor tissues and normal tissues. Three GEO expression profile datasets, GSE32665, GSE32863, and GSE43458, were used to compare the expression of GNG7 between LUAD tissues and normal tissues. Differential protein levels of GNG7 between LUAD tissues and normal tissues were analyzed by UALCAN (http://ualcan.path.uab.edu/) (Chandrashekar et al., 2022). A receiver operating characteristic (ROC) curve was used to evaluate the diagnostic significance of GNG7 using the plotROC R package (Version 1.17.0.1) (Robin et al., 2011).
Clinical statistical analysis on prognosis, model construction and evaluation
The correlation between GNG7 expression and survival in LUAD patients was analyzed by the PrognoScan database (http://www.prognoscan.org/) (Mizuno et al., 2009). Univariate Cox regression analysis, multivariate Cox regression analysis, logistic analysis and Kaplan-Meier (K-M) analysis were employed for prognosis analysis. The independent prognostic factors obtained from multivariate Cox regression analysis were employed to establish nomograms to predict survival probability for 1-, 2-, and 3-year. The calibration curves and nomograms were analyzed and plotted via the rms (version 6.2-0) and survival (version 3.2-10) package of R software. The calibration curves were graphically assessed by mapping the nomogram-predicted probabilities against the observed rates, and the 45-degree line represented the best predictive values. A concordance index (C-index) was used to determine the discrimination of the nomogram. According to the median risk score, patients were divided into a high-risk score group and a low-risk score group. The survival difference between the two groups was assessed by K-M survival curves. The model was compared with the two-by-two model consisting of independent prognostic factors screened from multivariate Cox regression, and ROC curves made by the timeROC R package (version 0.4) were used to assess the accuracy of the model predictions. The risk curve was used to demonstrate the application of the model in predicting clinical outcomes.
Gene set enrichment analysis
DESeq2 package (Version 1.26.0) was employed to identify the DEGs between GNG7-high and GNG7-low expression patients from TCGA datasets. The cut-off threshold was |log fold change (FC)|>1 and adjusted p-value <0.05. All the DEGs were presented in the volcano plots, and the correlation of some representative DEGs with GNG7 was presented in heatmaps. Gene Set Enrichment Analysis (GSEA) is a computational method for determining whether a defined set of genes shows statistically significant differences between two states. In the study, GSEA was performed by using the clusterProfiler R package (version 3.14.3) with c2 (c2. all.v7.0. entrez.gmt) from the Molecular Signatures Database (MSigDB) (Subramanian et al., 2005; Yu et al., 2012). Each analysis procedure was repeated 1000 times. The function or pathway termed with adjusted p-value <0.05 and false discovery rate (FDR) < 0.25 was considered statistically significant enrichment.
Tumor immune infiltration analysis, protein-protein interaction network analyses and the screening of hub genes
ESTIMATE algorithm was used to calculate the immune scores using the “estimate” R package (Yoshihara et al., 2013). Single-sample gene set enrichment analysis (ssGSEA) algorithm was used to assess the relative enrichment of the tumor tissue-infiltrating immune cells in LUAD (Hänzelmann et al., 2013). Based on an immune dataset for the 24 types of immunocytes, the relative enrichment score of every immunocyte was quantified from the gene expression profiles of each tumor sample (Bindea et al., 2013). In addition, we analyzed the differences in the enrichment of these 24 immune cells between the high and low GNG7 expression groups using ssGSEA. The Stat R package (Version 3.6.3) was employed to search for the genes related to GNG7 in LUAD. The correlation results were analyzed by the Spearman coefficient and the cut-off thresholds were |R| >0.4 and p-value <0.05. The list of immune-related genes was obtained from the ImmPort database. Genes associated with GNG7 were intersected with IRGs to obtain IRGs associated with GNG7. To further understand the interactions between IRGs associated with GNG7, we constructed a Protein-Protein Interaction (PPI) using the Search Tool for the Retrieval of Interacting Genes (STRING) (https://string-db.org/). An interaction with a combined score >0.4 was considered statistically significant. Cytoscape (Version 3.7.2) was used to visualize the network, while the cytoHubba plugin was used to rank genes within this network based upon their degree centrality values. Hub genes were considered to be those with the top 10 highest degree values. ClusterProfiler R package was employed to perform Gene Ontology (GO) function enrichment analyses. Furthermore, We used the GEPIA2 database (http://gepia2.cancer-pku.cn/#index) to analyze the correlation between GNG7 and B cell infiltration in LUAD tissues and normal tissues (Tang et al., 2019). Tumor Immune Estimation Resource 2.0 (TIMER2.0) database (http://timer.cistrome.org/) and the Kaplan-Meier plotter database (http://kmplot.com/analysis/) were used for the prognostic analysis of LUAD patients with different GNG7 expression and B cell infiltration (Li T. et al., 2020; Lánczky and Győrffy, 2021). We also used the TIMER database to analyze the correlation between GNG7 and immune cell markers in LUAD (Li et al., 2017).
DNA methylation analysis
To explore the possible mechanism of decreased expression of GNG7 in LUAD, we performed a differential methylation analysis of GNG7 between the normal and LUAD tissues using the UALCAN database. The cBio Cancer Genomics Portal (cBioPortal) (https://www.cbioportal.org), developed based on the TCGA database, was used to perform a correlation analysis between GNG7 mRNA expression and its methylation levels (Cerami et al., 2012; Gao et al., 2013). DNA methylation of GNG7 at CpG sites and the prognostic value of these CpG sites in LUAD were analyzed by MethSurv (https://biit.cs.ut.ee/methsurv/) (Modhukur et al., 2018).
Statistical analysis
All statistical analyses were conducted using R (Version 3.6.3). A part of the figures was plotted using the ggplot2 R package (Version 3.3.3). Dunn’s test, Kruskal–Wallis test, and logistic regression were used to analyze the clinicopathological features of GNG7 in LUAD. Kaplan-Meier survival analysis, univariate and multivariate Cox regression analysis were performed for prognostic analysis. In all analyses, the p-value<0.05 was considered statistically significant. The specific datasets, R packages, software and databases used in each part of this study are detailed in Supplementary Table S4.
Results
G protein subunit gamma 7 is found to be one of the key regulators of lung adenocarcinoma tightly related to the prognosis through large-scale screening
To find key regulators of LUAD, we conducted screening work based on differential expression analysis and prognostic analysis, and performed a series of analytical work based on the target gene. The workflow was shown in Figure 1 we applied the DESeq2 R package to screen for DEGs in LUAD based on TCGA-LUAD datasets. The results showed that there were 7741 up-regulated and 3783 down-regulated genes among the screened 11,524 DEGs (Figure 2A). Combined with further prognostic analysis, we found that among the DEGs, nine genes (HSD17B6, PXMP4, HLF, ADGRD1, CYP17A1, ESYT3, FCAMR, C11orf16, GNG7) were significantly and positively associated with LUAD prognostic indicators including OS, DSS and PFI (Figure 2B). Of note, although GNG7 has been reported to be differentially expressed in a variety of tumors, its roles in the initiation and progression of LUAD remain unclear. In the present study, we focused on GNG7 to explore the underlying mechanism and clinical significance in LUAD.
FIGURE 1
FIGURE 2
G protein subunit gamma 7 expression is downregulated in lung adenocarcinoma
To elucidate the expression pattern of GNG7 in cancers, we first evaluated the expression of GNG7 in 33 types of cancers by a systematic analysis based on the TCGA databases. The results showed that GNG7 expression was significantly down-regulated in 17 different tumors, including Bladder Urothelial Carcinoma (BLCA), Breast invasive carcinoma (BRCA), Colon adenocarcinoma (COAD), Lung adenocarcinoma (LUAD) while it was significantly up-regulated in Cholangiocarcinoma (CHOL), Liver hepatocellular carcinoma (LIHC) and Pheochromocytoma and Paraganglioma (PCPG) (Figure 2C). Then, the low expression of GNG7 in LUAD was further validated by using three GEO datasets (GSE32665, GSE32863, GSE43458) (Figures 2D–F). The paired analysis got the similar results (Figure 2G). In addition, the decreased protein level of GNG7 was also observed in LUAD using the UALCAN database (Figure 2H). Moreover, ROC curve analysis was employed to analyze the distinguishing efficacy of GNG7 between LUAD tissue and normal tissue. The area under the curve (AUC) of GNG7 is 0.871, suggesting that GNG7 may be an ideal biomarker to distinguish LUAD from normal tissue (Figure 2I). Together, these results indicated that GNG7 is lowly expressed in LUAD which may be a potential diagnostic marker for LUAD.
Association between clinicopathological characteristics and GNG7 expression in lung adenocarcinoma
To clarify the correlation between the expression of GNG7 and clinicopathological variables, we collected data from the TCGA database on 535 patients with LUAD. After data preprocessing, the relationship between gene expression profiles and clinicopathological characteristics of 513 LUAD patients was shown in the baseline data table (Supplementary Table S1). The results showed that low expression of GNG7 was positively associated with high T stage, Gender (male sex), poor primary therapy outcome and high pathologic stage of LUAD, while there were no significant associations between GNG7 expression and the other clinical factors such as N stage and M stage (Figures 3A–F). In line with these findings, the logistics regression analysis also revealed that GNG7 expression was significantly associated with T stage (OR = 0.403, 95% CI: 0.274-0.589, p < 0.001), Pathologic stage (OR = 0.538, 95% CI: 0.377-0.766, p < 0.001), Primary therapy outcome (OR = 0.311, 95% CI: 0.133-0.707, p = 0.006) and Gender (OR = 0.502, 95% CI: 0.352-0.713, p < 0.001) (Table 1).
FIGURE 3
TABLE 1
| Characteristics | Total (N) | Odds Ratio (OR) | p value |
|---|---|---|---|
| T stage (T2&T3&T4 vs. T1) | 510 | 0.403 (0.274−0.589) | ***<0.001 |
| N stage (N1&N2&N3 vs. N0) | 501 | 0.700 (0.482−1.014) | 0.060 |
| M stage (M1 vs. M0) | 369 | 0.493 (0.197−1.140) | 0.110 |
| Pathologic stage (Stage II&Stage III&Stage IV vs. Stage I) | 505 | 0.538 (0.377−0.766) | ***<0.001 |
| Primary therapy outcome (PD vs. SD) | 105 | 0.311 (0.133−0.707) | **0.006 |
| Gender (Male vs. Female) | 513 | 0.502 (0.352−0.713) | ***<0.001 |
| Race (Black or African American&White vs. Asian) | 446 | 2.701 (0.576−19.003) | 0.238 |
| Age (>65 vs. <=65) | 494 | 1.050 (0.738−1.496) | 0.785 |
| Residual tumor (R1&R2 vs. R0) | 361 | 0.565 (0.191−1.519) | 0.271 |
| Anatomic neoplasm subdivision (Right vs. Left) | 498 | 0.957 (0.669−1.371) | 0.812 |
| Anatomic neoplasm subdivision2 (Peripheral Lung vs. Central Lung) | 189 | 1.022 (0.555−1.894) | 0.943 |
| number_pack_years_smoked (>=40 vs. <40) | 351 | 0.863 (0.567−1.312) | 0.491 |
| Smoker (Yes vs. No) | 499 | 0.629 (0.378−1.035) | 0.070 |
GNG7 expression association with clinical pathological characteristics (logistic regression).
*p < 0.05; **p < 0.01; ***p < 0.001.
Significance of G protein subunit gamma 7 in clinical prognosis of lung adenocarcinoma and clinical subgroup analysis
We utilized data from the TCGA database to investigate the prognostic significance of GNG7 in LUAD. Kaplan-Meier survival analysis based on the TCGA-LUAD dataset revealed that low expression of GNG7 was associated with poor OS (HR = 0.51, 95% CI: 0.38-0.69, p < 0.001), DSS (HR = 0.56, 95% CI: 0.38-0.82, p = 0.003) and PFI (HR = 0.64, 95% CI: 0.49-0.85, p = 0.002) (Figures 3G–I). To further validate the prognostic value of GNG7 in LUAD, we utilized the PrognoScan database for further study. We included two of the GSE datasets (GSE31210 and GSE13213) in our analysis, where low GNG7 expression was significantly associated with the poorer prognosis (OS, HR = 0.21, 95% CI: 0,08-0.52, Cox p = 0.000748; RFS, HR = 0.25, 95% CI: 0,13-0.49, Cox p = 0.000069 in the GSE31210 dataset; OS, HR = 0.48, 95% CI: 0,34-0.67, Cox p = 0.000023 in the GSE13213 dataset) (Figures 3J–L).
Moreover, the univariate Cox regression analysis model showed that GNG7 expression level was significantly associated with OS (HR: 0.702; 95% CI: 0.599-0.822; p < 0.001) similar to T stage, N stage, M stage and Pathologic stage as well Primary therapy outcome and Residual tumor. Meanwhile, the multivariate Cox regression analysis also revealed that low expression of GNG7, similar to Primary therapy outcome and Residual tumor, was an independent risk factor for the prognosis of LUAD patients (Table 2). Collectively, these results suggest that low expression of GNG7 independently predicts poor prognosis for patients with LUAD.
TABLE 2
| Characteristics | Total (N) | Univariate analysis | Multivariate analysis | ||
|---|---|---|---|---|---|
| Hazard ratio (95% CI) | p value | Hazard ratio (95% CI) | p value | ||
| T stage (T2&T3&T4 vs. T1) | 501 | 1.668 (1.184−2.349) | **0.003 | 1.116 (0.629−1.978) | 0.708 |
| N stage (N1&N2&N3 vs. N0) | 492 | 2.606 (1.939−3.503) | ***<0.001 | 1.689 (0.758−3.764) | 0.200 |
| M stage (M1 vs. M0) | 360 | 2.111 (1.232−3.616) | **0.007 | 1.698 (0.674−4.280) | 0.262 |
| Pathologic stage (Stage II&Stage III&Stage IV vs. Stage I) | 496 | 2.975 (2.188−4.045) | ***<0.001 | 1.109 (0.471−2.610) | 0.812 |
| Primary therapy outcome (PD&SD&PR vs. CR) | 419 | 2.818 (2.004−3.963) | ***<0.001 | 3.662 (2.217−6.049) | ***<0.001 |
| Gender (Male vs. Female) | 504 | 1.060 (0.792−1.418) | 0.694 | ||
| Race (White vs. Asian&Black or African American) | 446 | 1.422 (0.869−2.327) | 0.162 | ||
| Age (>65 vs. <=65) | 494 | 1.228 (0.915−1.649) | 0.171 | ||
| Residual tumor (R1&R2 vs. R0) | 352 | 3.973 (2.217−7.120) | ***<0.001 | 3.670 (1.503−8.964) | **0.004 |
| Anatomic neoplasm subdivision (Right vs. Left) | 490 | 1.024 (0.758−1.383) | 0.878 | ||
| Anatomic neoplasm subdivision2 (Peripheral Lung vs. Central Lung) | 182 | 0.913 (0.570−1.463) | 0.706 | ||
| Number pack years smoked (>=40 vs. <40) | 345 | 1.038 (0.723−1.490) | 0.840 | ||
| Smoker (Yes vs. No) | 490 | 0.887 (0.587−1.339) | 0.568 | ||
| GNG7 | 504 | 0.702 (0.599−0.822) | ***<0.001 | 0.727 (0.561-0.943) | *0.016 |
Univariate analysis and multivariate analysis of the correlation between clinicopathological characteristics and OS in LUAD.
*p < 0.05; **p < 0.01; ***p < 0.001.
Given that multivariate Cox regression analysis identified low expression of GNG7 as an independent risk factor, we investigated the potential prognostic value of GNG7 in LUAD patients with different clinical subgroups. As shown in Figures 4A–C, low expression of GNG7 was associated with poor prognosis in stage N0, including OS (HR = 0.42, 95% CI: 0.27-0.67, p < 0.001), DSS (HR = 0.38, 95% CI: 0.21-0.69, p = 0.001) and PFI (HR = 0.57, 95% CI: 0.40-0.83, p = 0.003). However, there was no statistically significant correlation between GNG7 expression and prognosis in the N1&N2&N3 stage (p > 0.05) (Supplementary Figures S1A–C). In addition, low GNG7 expression was significantly associated with poor prognosis in LUAD patients in M0 stage, including OS (HR = 0.48, 95% CI: 0.33-0.68, p < 0.001), DSS (HR = 0.59, 95% CI: 0.37-0.94, p = 0.026), PFI (HR = 0.66, 95% CI: 0.47-0.92, p = 0.015) (Figures 4D–F). Nevertheless, no significant association was shown between GNG7 expression and prognosis in LUAD patients in the M1 stage (Supplementary Figures S1D–F). These results suggest that low expression of GNG7 is positively associated with the poor prognosis of LUAD patients without lymph node invasion and distal metastasis.
FIGURE 4
Construction and validation of a nomogram based on the independent clinical risk factors
To provide a quantitative approach to predicting the prognosis of LUAD patients, we constructed a prognostic nomogram to predict individual survival probability based on the expression levels of GNG7 and other independent clinical risk factors (Figure 4G). The calibration curve of the nomogram showed that the established lines of 1-, 2-, and 3-y survival highly matched the ideal line (the 45-degree line) (Figure 4H). In addition, the C-index of the prediction model reached 0.690 (0.659–0.720), indicating that the model had a reliable potential to predict the OS of LUAD patients. In addition, ROC curve analysis based on the three time points of 1-, 2-, and 3-Year showed that the Area under the curve (AUC) of this prediction model was higher than the AUC of the two-by-two model consisting of independent prognostic factors screened from multivariate Cox regression, indicating the superiority of the model (Supplementary Figures S2A–D). On the basis of the median risk score, patients were divided into a high-risk score group and a low-risk score group. Survival curve analysis revealed that the high-risk group had a significantly poorer prognosis compared to the low-risk group (HR = 2.59, 95% CI: 2.59 (1.72–3.89), p < 0.001) (Supplementary Figure S2E). Additionally, the risk curve indicated that the high-risk score group had higher mortality and worse prognosis than the low-risk score group (Supplementary Figure S2F).
Functional enrichment and pathway analysis of G protein subunit gamma 7-associated differentially expressed genes in lung adenocarcinoma
To investigate the biological functions and signaling pathways associated with GNG7, we examined the DEGs between GNG7-high and GNG7-low patients which were stratified based on the median GNG7 expression. Resultantly, 1403 mRNAs (492 upregulated and 911 downregulated), 962 lncRNAs (256 upregulated and 706 downregulated), and 21miRNAs (18 upregulated and 3 downregulated) were differently expressed in GNG7-high patients compared to GNG7-low ones (Figure 5A, Supplementary Figures S3A,C). Relative expression values of some representative DEGs between the two cohorts were shown in the form of heatmaps (Figure 5B, Supplementary Figures S3B,D). Strikingly, pathway enrichment analysis showed that the DEGs were most strongly enriched in the B cell receptor signaling pathway, T cell receptor signaling pathway and HIV infection allograft rejection which are highly related to the cellular immune response (Figures 5C–G). These data suggested that GNG7 may play an important role in regulating the tumor immune microenvironment of LUAD.
FIGURE 5
Correlation analysis between the expression of G protein subunit gamma 7 and immune cell infiltration in lung adenocarcinoma
As reported, tumor-associated immune cell infiltration has a close relationship with tumor development and the prognosis of patients. Then, we utilized the ESTIMATE algorithm to assess the correlation between GNG7 and the abundance of immune cell infiltration in LUAD. The results revealed that GNG7 expression was positively correlated with the abundance of immune cell infiltration in LUAD (Figure 6A). Specifically, further Spearman correlation analysis showed that among 24 immune cell subpopulations, GNG7 expression was positively correlated with most immune cell subsets, including Mast cell, DC, B cells, and CD8+ T cells, but negatively correlated with Th2 and Tgd cells (Figure 6B). Consistently, the ssGSEA analysis demonstrated that the infiltration levels of most of the immune cell subsets such as Mast cells, pDCs, B cells, NK cells and CD8+ T cells were remarkably increased in LUAD patients with GNG7 high expression compared to those with GNG7 low expression (Figure 6C). In keeping with this finding, GNG7 was significantly correlated with most immune markers of different immune cells, including CD8+T cell, B cell, Neutrophils, and Dendritic cells (Supplementary Table S2). Moreover, by correlation analysis, we identified immune-related genes (IRG) co-expressed with GNG7 and constructed a PPI network (Supplementary Figure S4A). We screened the top10 of hub genes and showed the correlation between these 10 genes and GNG7 in the form of scatter plots (Figure 6D, Supplementary Figures S4B–K). Additionally, we performed GO enrichment analysis to investigate the possible involvement of GNG7 in the immune response. The terms identified in the BP category showed that aberrantly expressed GNG7 was associated with antigen processing and presentation via MHC class II(MHCII), while in the CC category, the hub genes were significantly enriched in the MHCII protein complex and endoplasmic reticulum-related terms. Furthermore, the MF category revealed significant enrichment in GO terms related to the MHCII protein complex binding, MHCII receptor activity and peptide antigen binding, etc (Figures 6E,F). Together, these results suggest that GNG7 may contribute to the remodeling of the immune microenvironment in LUAD through promoting the infiltration of a variety of tumor-associated immune cells and influencing antigen presentation.
FIGURE 6
G protein subunit gamma 7 high expression with high B cell infiltration predicts a better prognosis of lung adenocarcinoma patients
Of the infiltrated immune cells increased in LUAD with GNG7 high expression, B cell infiltration attracts our attention as the relatively less knowledge of this cell type in tumor immunotherapy currently. Our results showed that GNG7 was closely associated with the level of B cell infiltration in LUAD (Supplementary Figure S5A). Specifically, the level of B cell infiltration was significantly elevated in the GNG7 high expression group compared with that in the GNG7 low expression group (Supplementary Figure S5B). In addition, we respectively investigated the correlation of GNG7 with B cells in the tumor and normal tissues. Strikingly, GNG7 showed a strong positive correlation with the B cell marker genes CD19 and CD79A in LUAD tissues (Figure 7A). In contrast, the correlation of GNG7 with the B cell markers CD19 and CD79A did not reach statistical significance in normal tissues (Figure 7B). These results suggest that GNG7 expression may promote B cell infiltration in the context of LUAD. Meanwhile, through KM plot database analysis, we found that patients with high GNG7 expression tended to predict a better prognosis in the B cell enriched group but not in the B cell decreased group (Figures 7C,D). Such finding was further corroborated by the analysis using the TIMER2.0 database, implying that high GNG7 expression corresponded to a better prognosis for LUAD patients in the context of enriched B cell infiltration. Furthermore, we found that high infiltration levels of B cells in the presence of consistent levels of GNG7 expression corresponded to a good prognosis in patients with LUAD (Figure 7E). Taken together, it is reasonable to suggest that GNG7 may have improved patient prognosis by promoting B cell infiltration.
FIGURE 7
G protein subunit gamma 7 dysregulation is associated with aberrant DNA methylation
Considering the importance of DNA methylation in regulating gene expression, we tested whether aberrant DNA methylation occurs in GNG7 gene in LUAD. By analyzing the data from the UALCAN database, we found that the methylation level of GNG7 was significantly higher in the tumor group compared to that in the normal group (Figure 8A). Next, the correlation analysis based on the cBioPortal database showed that GNG7 expression was significantly negatively correlated with methylation (Figure 8B). To further investigate the methylation of GNG7 in LUAD, we analyzed the methylation levels of different CpG sites of GNG7 in LUAD patients using the MethSurv database and presented them in the form of heat maps (Figure 8C). The results revealed that several CpG sites of GNG7 exhibited high methylation in LUAD patient samples, including cg19477361, cg21462934, and cg27181295. Prognostic analysis showed that the above CpG sites with highly methylated levels were associated with poor prognosis in LUAD (Figures 8D–F). These results suggest that the low expression of GNG7 in LUAD may be partly due to the methylation modification of the abovementioned CpG sites and plays a key role in tumor progression.
FIGURE 8
Discussion
As one of the most common malignancies worldwide, the prognosis of LUAD patients remains very gloomy due to the lack of effective biomarkers for early diagnosis and effective treatment for advanced patients. Thus, intense research has been focused on deciphering the pathogenesis and searching for effective diagnostic and therapeutic approaches as well as prognostic markers to improve the prognosis of patients with LUAD in the last several years. Indeed, a growing number of potential biomarkers for LUAD have been identified, such as PPP1R14D, lncRNA-Ac068228, IFITM1, and so on (Koh et al., 2019; Jiang et al., 2022; Tian Y. et al., 2022). However, most of these biomarkers are associated with the increase in cell numbers resulting from cell division (cell proliferation), programmed cell death (apoptosis), and tumor angiogenesis, while few with tumor immune microenvironment. Accumulating evidence has shown that not only the characteristics of tumor cells but also the tumor microenvironment, especially Tumor infiltrating immune cells (TIICs), plays critical roles in the tumorigenesis and progression of LUAD (Hinshaw and Shevde, 2019). In recent years, with the further understanding of the mechanism of tumor immune infiltration, tumor immunotherapies, such as the immune checkpoint inhibitors (ICIs), have had a revolutionary impact on the treatment of LUAD (Bagchi et al., 2021). However, only a small percentage of patients achieved a durable immune response after treatment. The mechanisms of LUAD development are far from being elucidated. It is of great necessity to further clarify the molecular basis of LUAD and explore contributing factors as well as sensitive diagnosis and prognosis biomarkers of immunotherapy response to improve patient outcomes (Wu and Shih, 2018; Singh et al., 2020).
In the current study, we performed comprehensive bioinformatics analyses to explore the potential key molecules involved in the development of LUAD. Through screening and identification, GNG7 was demonstrated to be lowly expressed in LUAD and had a good diagnostic performance. In addition, low expression of GNG7 was positively associated with the poor clinicopathological characteristics such as poor primary therapy outcome and high pathologic stage of LUAD, implying the tumor suppressive roles of GNG7 in LUAD.
As a subunit of heterotrimeric G protein, GNG7 has been reported to be tightly related to carcinogenesis. GNG7 is frequently downregulated in various cancers including pancreatic cancer, esophageal cancer and clear cell renal cell carcinoma (Shibata et al., 1998; Ohta et al., 2008; Xu et al., 2019). It is worth noting that GNG7 has been identified as one of the hub genes in an eight-gene prognostic signature model and a four-gene panel predicting overall survival for LUAD (Li C. et al., 2020; Ma et al., 2021). However, the role of GNG7 as an independent prognostic factor in LUAD has not been fully elucidated. Here, our study based on GNG7 is similar to, but more distinctive from, the newly identified biomarkers for LUAD in the latest literatures. In our study, we demonstrated GNG7 as an independent prognostic risk factor in LUAD (Supplementary Table S4) (Wan et al., 2021; Tian W. et al., 2022; Zhang et al., 2022; Zhou et al., 2022). Moreover, we constructed a nomogram combined with other clinical independent prognostic risk factors to predict the prognosis of LUAD patients reliably.
Recent studies have reported that GNG7 inhibited the progression of LUAD by inhibiting E2F1 and Hedgehog signaling, but the exact mechanism by which it regulates the development of LUAD is largely unknown (Zhao et al., 2021; Zheng et al., 2021). Our GSEA and GO enrichment analysis found that GNG7 may be involved in regulating the TME of LUAD and antigen processing and presentation via MHCII. Especially, high expression of GNG7 corresponded to increased infiltration level of several immune cells including B cells. Over the past decades, intense research has focused on the roles of T cells in immune regulation in the TME (Joyce and Fearon, 2015). More recently, there is increasing evidence supporting a critical role for B cells in tumor immunology (Bruno, 2020; Cabrita et al., 2020; Helmink et al., 2020). However, there is limited understanding of the biological contributors to the B cell infiltration in the TME. As the key immune cell in humoral immunity, B cells express a large number of MHCII molecules and are important antigen-presenting cells. In the present study, we found that GNG7 may be involved in regulating MHCII-mediated antigen processing and presentation. Furthermore, we found that GNG7 expression may promote B cell infiltration as evidenced by that low GNG7 expression was negatively correlated with B cell infiltration. Strikingly, the results of the prognostic analysis indicated patients with high GNG7 expression tended to predict a better prognosis in the context of enriched B cell infiltration and high infiltration levels of B cells regardless of GNG7-low or GNG7-high expression corresponded to a good prognosis in patients with LUAD. Our results indicated that GNG7 may exert its tumor suppressive roles in LUAD by promoting B cell infiltration and GNG7 expression together with B cell infiltration may be a powerful predictive signature for prognosis and immunotherapy response in LUAD, although the detailed function and mechanism need further in-depth investigation both in vitro and in vivo.
Finally, in this study, we further explored the mechanism of GNG7 low expression in LUAD. To our knowledge, other than miR-19b-3p which was reported to target GNG7 directly and significantly decrease the mRNA level of GNG7, little is known about the mechanism of GNG7 dysregulation in LUAD (Zhao et al., 2021). Given that DNA methylation of CpG islands is known to be a repressive mark of gene expression, we assessed the methylation level of GNG7 in LUAD (Kulis and Esteller, 2010). As expected, elevated methylation level of GNG7 was observed in tumor tissues which may be responsible for the low expression of GNG7 in LUAD. Interestingly, hypermethylation of GNG7 was associated with poor prognosis in patients with LUAD, which is consistent with the prognostic value of GNG7 mRNA expression. These results suggest the importance of DNA methylation in regulating GNG7 expression in LUAD.
In summary, this study revealed for the first time that GNG7 may be involved in regulating the immune microenvironment in LUAD and influence tumor development and patient prognosis at least partly by regulating the B cell infiltration. GNG7 may be not only a potential diagnostic biomarker for LUAD but also a promising predictive signature for prognosis and immunotherapy response for patients with LUAD. Nevertheless, as data from this study were mainly obtained from open databases, more LUAD patient samples are needed to confirm the clinical prognostic value of GNG7. Moreover, the effects of GNG7 on immune cell recruitment and infiltration as well as immunotherapy response are needed to be investigated deeply at the cellular and molecular levels and in future clinical trials.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: TCGA database: https://portal.gdc.cancer.gov/ and GEO database: https://www.ncbi.nlm.nih.gov/geo/.
Author contributions
HY and QW designed the study. QW performed the acquisition and analysis of data. QW wrote the manuscript. PZ, TM, BJ, and HY edited the manuscript. All authors read and approved the final manuscript.
Funding
This work was funded by National Natural Science Foundation of China (Nos 31870781) and Shandong Province Science and Technology Development Plan Foundation for Traditional Chinese Medicine (Nos 2017-047).
Acknowledgments
We would like to thank the National Natural Science Foundation of China (Nos 31870781) and the Shandong Provincial Science and Technology Development Program Fund for Traditional Chinese Medicine (Nos 2017-047) for their financial support, as well as TCGA and GEO databases for providing platforms and contributors to upload their meaningful datasets. It should be noted that this article has been shared as a preprint on the Research Square platform under the DOI number: https://doi.org/10.21203/rs.3.rs-1702007/v1.
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/fgene.2022.984575/full#supplementary-material
References
1
BagchiS.YuanR.EnglemanE. G. (2021). Immune checkpoint inhibitors for the treatment of cancer: clinical impact and mechanisms of response and resistance. Annu. Rev. Pathol.16, 223–249. 10.1146/annurev-pathol-042020-042741
2
BindeaG.MlecnikB.TosoliniM.KirilovskyA.WaldnerM.ObenaufA. C.et al (2013). Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity39 (4), 782–795. 10.1016/j.immuni.2013.10.003
3
BrunoT. C. (2020). New predictors for immunotherapy responses sharpen our view of the tumour microenvironment. Nature577 (7791), 474–476. 10.1038/d41586-019-03943-0
4
CabritaR.LaussM.SannaA.DoniaM.Skaarup LarsenM.MitraS.et al (2020). Tertiary lymphoid structures improve immunotherapy and survival in melanoma. Nature577 (7791), 561–565. 10.1038/s41586-019-1914-8
5
CeramiE.GaoJ.DogrusozU.GrossB. E.SumerS. O.AksoyB. A.et al (2012). The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov.2 (5), 401–404. 10.1158/2159-8290.Cd-12-0095
6
ChandrashekarD. S.KarthikeyanS. K.KorlaP. K.PatelH.ShovonA. R.AtharM.et al (2022). UALCAN: an update to the integrated cancer data analysis platform. Neoplasia25, 18–27. 10.1016/j.neo.2022.01.001
7
GaoJ.AksoyB. A.DogrusozU.DresdnerG.GrossB.SumerS. O.et al (2013). Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal.6 (269), pl1. 10.1126/scisignal.2004088
8
HänzelmannS.CasteloR.GuinneyJ. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma.14, 7. 10.1186/1471-2105-14-7
9
HartmannS.SzaumkesselM.SalaverriaI.SimonR.SauterG.KiwerskaK.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. 10.1007/s13353-011-0079-4
10
HelminkB. A.ReddyS. M.GaoJ.ZhangS.BasarR.ThakurR.et al (2020). B cells and tertiary lymphoid structures promote immunotherapy response. Nature577 (7791), 549–555. 10.1038/s41586-019-1922-8
11
HinshawD. C.ShevdeL. A. (2019). The tumor microenvironment innately modulates cancer progression. Cancer Res.79 (18), 4557–4566. 10.1158/0008-5472.Can-18-3962
12
JiangX.ChenM.DuJ.BiH.GuoX.YangC.et al (2022). LncRNA-AC068228.1 is a novel prognostic biomarker that promotes malignant phenotypes in lung adenocarcinoma. Front. Oncol.12, 856655. 10.3389/fonc.2022.856655
13
JoyceJ. A.FearonD. T. (2015). T cell exclusion, immune privilege, and the tumor microenvironment. Science348 (6230), 74–80. 10.1126/science.aaa6204
14
KohY. W.HanJ. H.JeongD.KimC. J. (2019). Prognostic significance of IFITM1 expression and correlation with microvessel density and epithelial-mesenchymal transition signature in lung adenocarcinoma. Pathol. Res. Pract.215 (7), 152444. 10.1016/j.prp.2019.152444
15
KulisM.EstellerM. (2010). DNA methylation and cancer. Adv. Genet.70, 27–56. 10.1016/b978-0-12-380866-0.60002-2
16
LánczkyA.GyőrffyB. (2021). Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J. Med. Internet Res.23 (7), e27633. 10.2196/27633
17
LiC.LongQ.ZhangD.LiJ.ZhangX. (2020a). Identification of a four-gene panel predicting overall survival for lung adenocarcinoma. BMC Cancer20 (1), 1198. 10.1186/s12885-020-07657-9
18
LiT.FanJ.WangB.TraughN.ChenQ.LiuJ. S.et al (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res.77 (21), e108–e110. 10.1158/0008-5472.Can-17-0307
19
LiT.FuJ.ZengZ.CohenD.LiJ.ChenQ.et al (2020b). TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res.48 (W1), W509–w514. 10.1093/nar/gkaa407
20
LiuJ.JiX.LiZ.YangX.WangW.ZhangX. (2016). G protein γ subunit 7 induces autophagy and inhibits cell division. Oncotarget7 (17), 24832–24847. 10.18632/oncotarget.8559
21
LiuJ.LichtenbergT.HoadleyK. A.PoissonL. M.LazarA. J.CherniackA. D.et al (2018). An integrated TCGA pan-cancer clinical data Resource to drive high-quality survival outcome analytics. Cell173 (2), 400–416. e411. 10.1016/j.cell.2018.02.052
22
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15 (12), 550. 10.1186/s13059-014-0550-8
23
MaW.LiangJ.LiuJ.TianD.ChenZ. (2021). Establishment and validation of an eight-gene metabolic-related prognostic signature model for lung adenocarcinoma. Aging (Albany NY)13 (6), 8688–8705. 10.18632/aging.202681
24
MartinezM.MoonE. K. (2019). CAR T cells for solid tumors: new strategies for finding, infiltrating, and surviving in the tumor microenvironment. Front. Immunol.10, 128. 10.3389/fimmu.2019.00128
25
MizunoH.KitadaK.NakaiK.SaraiA. (2009). PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med. Genomics2, 18. 10.1186/1755-8794-2-18
26
ModhukurV.IljasenkoT.MetsaluT.LokkK.Laisk-PodarT.ViloJ. (2018). MethSurv: a web tool to perform multivariable survival analysis using DNA methylation data. Epigenomics10 (3), 277–288. 10.2217/epi-2017-0118
27
MolinaJ. R.YangP.CassiviS. D.SchildS. E.AdjeiA. A. (2008). Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin. Proc.83 (5), 584–594. 10.4065/83.5.584
28
OhtaM.MimoriK.FukuyoshiY.KitaY.MotoyamaK.YamashitaK.et al (2008). Clinical significance of the reduced expression of G protein gamma 7 (GNG7) in oesophageal cancer. Br. J. Cancer98 (2), 410–417. 10.1038/sj.bjc.6604124
29
PetitprezF.MeylanM.de ReynièsA.Sautès-FridmanC.FridmanW. H. (2020). The tumor microenvironment in the response to immune checkpoint blockade therapies. Front. Immunol.11, 784. 10.3389/fimmu.2020.00784
30
QiX.QiC.QinB.KangX.HuY.HanW. (2020). Immune-Stromal score signature: novel prognostic tool of the tumor microenvironment in lung adenocarcinoma. Front. Oncol.10, 541330. 10.3389/fonc.2020.541330
31
RobinX.TurckN.HainardA.TibertiN.LisacekF.SanchezJ. C.et al (2011). pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinforma.12, 77. 10.1186/1471-2105-12-77
32
SaitoM.SuzukiH.KonoK.TakenoshitaS.KohnoT. (2018). Treatment of lung adenocarcinoma by molecular-targeted therapy and immunotherapy. Surg. Today48 (1), 1–8. 10.1007/s00595-017-1497-7
33
SchwindingerW. F.MirshahiU. L.BaylorK. A.SheridanK. M.StaufferA. M.UsefofS.et al (2012). Synergistic roles for G-protein γ3 and γ7 subtypes in seizure susceptibility as revealed in double knock-out mice. J. Biol. Chem.287 (10), 7121–7133. 10.1074/jbc.M111.308395
34
ShibataK.MoriM.TanakaS.KitanoS.AkiyoshiT. (1998). Identification and cloning of human G-protein gamma 7, down-regulated in pancreatic cancer. Biochem. Biophys. Res. Commun.246 (1), 205–209. 10.1006/bbrc.1998.8581
35
ShibataK.TanakaS.ShiraishiT.KitanoS.MoriM. (1999). G-protein gamma 7 is down-regulated in cancers and associated with p 27kip1-induced growth arrest. Cancer Res.59 (5), 1096–1101.
36
SinghS. S.DahalA.ShresthaL.JoisS. D. (2020). Genotype driven therapy for non-small cell lung cancer: resistance, pan inhibitors and immunotherapy. Curr. Med. Chem.27 (32), 5274–5316. 10.2174/0929867326666190222183219
37
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A.102 (43), 15545–15550. 10.1073/pnas.0506580102
38
SungH.FerlayJ.SiegelR. L.LaversanneM.SoerjomataramI.JemalA.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. 10.3322/caac.21660
39
TangZ.KangB.LiC.ChenT.ZhangZ. (2019). GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res.47 (W1), W556–w560. 10.1093/nar/gkz430
40
TaubeJ. M.GalonJ.ShollL. M.RodigS. J.CottrellT. R.GiraldoN. A.et al (2018). Implications of the tumor immune microenvironment for staging and therapeutics. Mod. Pathol.31 (2), 214–234. 10.1038/modpathol.2017.156
41
TianW.ZhouJ.ChenM.QiuL.LiY.ZhangW.et al (2022a). Bioinformatics analysis of the role of aldolase A in tumor prognosis and immunity. Sci. Rep.12 (1), 11632. 10.1038/s41598-022-15866-4
42
TianY.GuanL.QianY.WuY.GuZ. (2022b). Effect of PPP1R14D gene high expression in lung adenocarcinoma knocked out on proliferation and apoptosis of DMS53 cell. Clin. Transl. Oncol. Online ahead of print. 10.1007/s12094-022-02842-7
43
TravisW. D. (2011). Pathology of lung cancer. Clin. Chest Med.32 (4), 669–692. 10.1016/j.ccm.2011.08.005
44
WanQ.QuJ.LiL.GaoF. (2021). Guanylate-binding protein 1 correlates with advanced tumor features, and serves as a prognostic biomarker for worse survival in lung adenocarcinoma patients. J. Clin. Lab. Anal.35 (2), e23610. 10.1002/jcla.23610
45
WuS. G.ShihJ. Y. (2018). Management of acquired resistance to EGFR TKI-targeted therapy in advanced non-small cell lung cancer. Mol. Cancer17 (1), 38. 10.1186/s12943-018-0777-1
46
WuZ.OuyangC.PengL. (2020). An immune scores-based nomogram for predicting overall survival in patients with clear cell renal cell carcinoma. Med. Baltim.99 (34), e21693. 10.1097/md.0000000000021693
47
XuS.ZhangH.LiuT.ChenY.HeD.LiL. (2019). G Protein γ subunit 7 loss contributes to progression of clear cell renal cell carcinoma. J. Cell. Physiol.234 (11), 20002–20012. 10.1002/jcp.28597
48
YoshiharaK.ShahmoradgoliM.MartínezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun.4, 2612. 10.1038/ncomms3612
49
YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS16 (5), 284–287. 10.1089/omi.2011.0118
50
ZhangY.WangQ.ZhuT.ChenH. (2022). Identification of cigarette smoking-related novel biomarkers in lung adenocarcinoma. Biomed. Res. Int.2022, 9170722. 10.1155/2022/9170722
51
ZhaoX.ZhangX. C.ZangK.YuZ. H. (2021). MicroRNA miR-19b-3p mediated G protein γ subunit 7 (GNG7) loss contributes lung adenocarcinoma progression through activating Hedgehog signaling. Bioengineered12 (1), 7849–7858. 10.1080/21655979.2021.1976896
52
ZhengH.TianH.YuX.RenP.YangQ. (2021). G protein gamma 7 suppresses progression of lung adenocarcinoma by inhibiting E2F transcription factor 1. Int. J. Biol. Macromol.182, 858–865. 10.1016/j.ijbiomac.2021.04.082
53
ZhouH.YaoY.LiY.GuoN.ZhangH.WangZ.et al (2022). Identification of small nucleolar RNA SNORD60 as a potential biomarker and its clinical significance in lung adenocarcinoma. Biomed. Res. Int.2022, 5501171. 10.1155/2022/5501171
Summary
Keywords
immune microenvironment 1, GNG7 2, prognosis 3, lung adenocarcinoma 4, bioinformatics analysis 5
Citation
Wei Q, Miao T, Zhang P, Jiang B and Yan H (2022) Comprehensive analysis to identify GNG7 as a prognostic biomarker in lung adenocarcinoma correlating with immune infiltrates. Front. Genet. 13:984575. doi: 10.3389/fgene.2022.984575
Received
02 July 2022
Accepted
08 August 2022
Published
09 September 2022
Volume
13 - 2022
Edited by
Saurav Mallik, Harvard University, United States
Reviewed by
Loveleen Gaur, Amity University, India
Soumita Seth, Aliah University, India
Updates
Copyright
© 2022 Wei, Miao, Zhang, Jiang and Yan.
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: Hua Yan, 71001643@sdutcm.edu.cn
This article was submitted to Computational Genomics, a section of the journal Frontiers in Genetics
Disclaimer
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.