- 1Department of Oncology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
- 2Department of Medical Oncology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, China
Background: Stomach adenocarcinoma (STAD) is one of the most common tumors. Tumor mutation burden (TMB) has been linked to immunotherapy response. We wanted to see if there was any link between TMB and cancer prognosis.
Methods: The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases were used to obtain mutation data, gene expression profiles, and clinical data. We looked at the differences in gene expression and immune markers between low and high TMB groups, built an immune prognostic model, and created a dynamic nomograph App that may be used in the clinic. Simultaneously, We ran the immunotherapy prediction and model comparison at the same time. Finally, model gene mutation and copy number variation (CNV) were displayed. The cellular functional experiments were used to investigate the potential role of GLP2R in gastric cancer.
Results: Firstly, basic mutation information and differences in immune infiltration in STAD are revealed. Secondly, the prognostic model developed by us has good accuracy, and the corresponding dynamic nomograph Apps online and immunotherapy prediction facilitate clinical transformation. Furthermore, GLP2R knockdown significantly inhibited the proliferation, migration of gastric cancer cells in vitro.
Conclusion: Our findings imply that TMB plays a significant role in the prognosis of STAD patients from a biological perspective. GLP2R may serve as a potential target for gastric cancer.
Background
The third most common cancer-related fatality is gastric cancer (GC), among which STAD is the most common type of pathological tissue (Bray et al., 2018). Although many therapeutic ways have been available, the high recurrence rate causes a heavy economic burden in both family and healthcare systems (Mehmedagic et al., 2016; Lee et al., 2017). The cause of STAD is still unclear; many factors, including H. pylori infection, smoking, environmental factors, poison contact history, etc. are associated with the occurrence of STAD (Naumann, 2005). Although clinicians have enriched therapeutic choices in recent years, many clinical obstacles are still unresolved. Immunotherapy has been a newly developed method to treat tumors by targeting PD-1, PD-L1, CTLA-4, etc (Kim and Seo, 2018; Powles and Morrison, 2018; Zhang et al., 2018).
First-line clinical trials of immunotherapy in combination with conventional therapy have shown improved clinical benefit and survival in patients with gastric cancer, especially among pretreated patients. This has contributed to the accelerated approval of some checkpoint inhibitors. However, the therapeutic effects of single-drug immune checkpoint inhibitors do not seem to meet expectations. Therefore, new immunotherapy algorithms should develop more efficient predictive biomarkers to distinguish between gastric tumor subsets with different clinical responses (Lazăr et al., 2018; Coutzac et al., 2019; Liang et al., 2021).
Previous papers have pointed out the correlation between immunotherapy response and TMB (Rizvi et al., 2015; He et al., 2019). Gene mutations in tumor tissues may produce new antigens through transcription and translation, thus being recognized and targeted by the immune system (Matsushita et al., 2012; Riaz et al., 2016). Not all mutations produce immunogenicity, and immune cells can only recognize a few mutated tumor antigens (Snyder and Chan, 2015). The more variations the tumor has, the more antigens it forms. Higher TMB has a tendency to generate more neoantigens, enhancing the immunogenicity of tumors and improving the response to immunotherapy in clinical applications (Rizvi et al., 2015). Therefore, TMB can be used as a biomarker to evaluate neoantigen load in tumors.
TMB can be utilized as a biomarker to predict the survival rate of patients with advanced gastric cancer after immunotherapy, which aids doctors in making the best decisions possible (Wang et al., 2019; Kim et al., 2020). In terms of mechanism, TGFB2 may have a role in the epithelial-mesenchymal transition (EMT) and TMB in gastric cancer, making it a possible therapeutic target (Yang et al., 2020). CXCR4 may also influence gastric cancer growth and prognosis by influencing immune infiltration, TMB, cytolytic activity, tumor purity, and treatment sensitivity (Li et al., 2020). In general, the TMB trial in stomach cancer warrants additional investigation.
TCGA has currently mapped the mutations in the human cancer genome, providing a wealth of mutation and expression profile data to researchers all around the world. We used STAD samples from the TCGA database to find differentially expressed genes (DEGs) across the high and low TMB groups, as well as investigate the relationship between immune cell infiltration features and TMB, creating and verifying an immune prognostic model. In the end, a series of functional experiments were carried out on GLP2R.
Materials and Methods
Data Acquisition
We downloaded various data of STAD patients from the TCGA data portal (https://tcgadata.nci.nih.gov/tcga/) (Hutter and Zenklusen, 2018), including mutation data, gene expression profiles, and clinical data. The STAD probe matrix file (GSE84433 series matrix) and platform file (GPL6947-13512) were also obtained from the GEO database (Barrett et al., 2007). The patient selection criteria for this study was samples of initial gastric adenocarcinoma after primary surgical resection (patients with missing clinical information were not included). To summarize the mutation data, The mutation data were summarized using the Maftools software package (Mayakonda et al., 2018). The relevant parameters of tumor-specific mutant genes were computed to obtain TMB. STAD samples were separated into high and low TMB groups using the median of TMB as the critical value. Then we looked at the relationship between TMB and survival and other clinical factors. Gene expression data from high and low TMB groups were analyzed using limma package (Ritchie et al., 2015). The fold change >2 and the false discovery rate (FDR) < 0.05 were our screening criterion for DEGs.
DEGs Enrichment Analysis and Immune Infiltration Assessment
We used GSEA 4.0.3 software to perform Gene Set Enrichment Analysis (GSEA) and reported the first five Gene Ontology (GO) keywords and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways which enrich most significantly in the high TMB group (compared to the low TMB group) (Subramanian et al., 2005).
The fraction of invading immune cells was computed by applying the CIBERSORT algorithm, and the corresponding violin plot was created. Using the CIBERSORT algorithm with 1,000 permutations and the LM22 signature, an accurate analysis of the immune cells in STAD samples was carried out (Chen et al., 2018).
Visualization of Mutation Data and ssGSEA Analysis
We performed ssGSEA analysis in TCGA-STAD samples using the GSVA package to determine the immune activity of 23 immune-related genesets, and the correction results ranged between 0 and 1. The heat map (pheatmap R package) and violin plots (ggpubr R package) of the tumor microenvironment were created based on the correlation analysis using the ESTIMATE algorithm between high and low risk group. To clarify the relationship between the risk grouping and immune infiltration, we also used the ssGSEA method to calculate the tumor immune cell infiltration score for each STAD sample in the TCGA database, then used the “Bioconductor Limma” R package to run the differential analysis of immune scoring and immune typing, and finally drew the box plot.
Construction and Multiple Validations of Immune Prognostic Model
We obtained differentially expressed immune genes by intersecting the list of immune-related genes from the Immunology Database and Analysis Portal (ImmPort Database) (Bhattacharya et al., 2014) with the preceding DEGs. We randomly separated TCGA samples into two groups: a training group and an internal validation group after correlating the expression of differentially immune genes with survival time. Univariate COX analysis (p < 0.01) identified survival-related immune genes in the training group, and using the R package glmnet, the least absolute shrinkage and selection operator (LASSO) (Tibshirani, 1997) excluded the high correlation genes to avoid the model from overfitting. Finally, the best prognostic model was created using stepwise multivariate Cox regression analysis.
The developed prognostic model formula was used to determine the patient riskScores for the training, internal validation, total, and external validation (GSE84433) groups. We divided the patients in the training, internal validation, total, and external validation groups into a high-risk group and a low-risk group using the median riskScore of the training group as the threshold, plotting the survival curve and Receiver Operating Characteristic (ROC) curve of the training, internal validation, total, and external validation groups with R x64 3.6.3 software. Ultimately, univariate and multivariate prognosis analyses were done on the entire group (p < 0.05) to see if the model’s riskScore may be an independent prognostic factor.
The differences in model gene expression were evaluated between the high and low risk groups, with the ggpubr package drawing model gene boxplots to compute the expression differences. To conduct multiple gene survival analysis for all model genes, we used the PROGgeneV2 online program (Goswami and Nakshatri, 2014) and picked the GSE62254 dataset (Cristescu et al., 2015). We created a dynamic nomograph App online based on the aforesaid prognostic model to aid in the rapid calculation of patient survival rates in clinical practice. In order to demonstrate the application method, we randomly selected a low-risk group sample and a high-risk group sample from the training group, and calculated survival probability by inputting gene expression in the dynamic nomograph App.
Mutation and CNV of Model Genes
We entered the cBioportal website (Cerami et al., 2012) and selected the study (Stomach Adenocarcinoma TCGA PanCancer data) to download the mutation status of model genes. Finally, we investigated the association between CNV of the model genes and immune cells infiltration level, as well as the correlation between immune cells infiltration level and survival of STAD patients by using the Tumor Immune Estimation Resource (TIMER) database (Li et al., 2017).
Immunotherapy Benefit Evaluation and Model Comparison
To explore whether the above prognostic model can evaluate the efficacy of immunotherapy for patients, we analyzed a series of immunotherapy biomarkers. We examined the differences between the high and low risk groups by uploading the TCGA sample expression profile into the Tumor immune dysfunction and exclusion (TIDE) database (Jiang et al., 2018) and obtaining TIDE, Microsatellite instability (MSI), Dysfunction, and Exclusion scores for each sample.
Furthermore, 5-years ROC curves were used to evaluate the riskScore’s prognostic prediction efficiency with TIDE and tumor inflammation signature (TIS) scores (Danaher et al., 2018), demonstrating the correctness of the prognostic model we developed.
Drug Sensitivity Testing
We downloaded transcriptome data from the CellMiner database (https://discover.nci.nih.gov/cellminer/) and FDA-certified drug sensitivity–related data to clarify the influence of model genes in the prognostic model on drug sensitivity and tolerance. To investigate the relationship between gene expression and drug sensitivity, the Pearson correlation test was used. Then, we used the “pRRophetic” R package to draw a box plot and analyze the differences in drug sensitivity between high and low risk groups.
Cell Culture and siRNA Treatment
Human gastric cancer cell lines (BGC823,MKN45) were cultured in RPMI 1640 medium (HyClone, United States)supplemented with 10% fetal bovine serum (FBS, Gibco, United States). These cells were maintained at 37°C under an atmosphere of 5% CO2. GLP2R siRNA were synthesized by RiboBio (Guangzhou, China) and transfected into cells using Lipofectamine 3,000 (Invitrogen, California, United States) according to the manufacturer’s instructions. The cells were cultured in basic RPMI 1640 medium for 12–16 h before the media was for RPMI 1640 medium containing FBS.
Quantitative Real-Time PCR
RNA was extracted from gastric cancer cell lines (BGC823,MKN45) interfered with si-GLP2R and NC as control. cDNA was synthesized for real-time PCR adopting SYBR Green qPCR mix (Vazyme, China). The primers are listed as following: GLP2R -Forward: TCCTGGAAATGTCTCTGTACCC; GLP2R—Reverse: GGCGTTCTCTATCGTCTGCC; GAPDH-Forward:GACCACAGTCCATGCCATCA; GAPDH-reverse:GTCAAAGGTGGAGGAGTGGG.
Western Blot Analysis
RIPA lysis buffer (Servicebio, China) containing PMSF (Servicebio, China) was used to collect proteins from BGC823 and MKN45 cells. 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) was applied in separating protein samples and a polyvinylidene fluoride membrane (PVDF) membrane (Invitrogen, Carlsbad, United States) was used to transfer the separated protein. The membrane was blocked in 5% skim milk at room temperature for 1 h in a shaker, and then incubated with the primary antibody:GLP2R (Abconlal, CAT# A6602, 1:1,000), and GAPDH (Proteintech, CAT# 60004-1-Ig, 1:20,000) at 4°C overnight and subsequent incubation with the secondary antibody for 1 h. Finally, the indicated proteins were visualized by West Pico plus Chemiluminescent Substrate (Thermo Fisher Scientific, United States).
Cell Proliferation and Colony Formation Assays
After transfection with GLP2R siRNA for 48 h, BGC823 and MKN45 cells were cultured in 96-well plates (3,000 cells/well in 100 µl RPMI 1640 medium). The proliferative capacity of treated cells was detected at 24, 48, 72, 96, and 120 h. 10% Cell Counting Kit-8 (CCK8) reagent (Bio-sharp, Hefei, China) was added to each plate according to the kit instructions, and the OD450 value was analyzed by a microplate reader (BioTek, United States). Regarding colony formation experiment, 1,000 cells were seeded in cell culture plates and allowed to grow until visible colonies formed. Then we used methanol to fix clones 15 min, 1% crystal violet to stain clones 20 min, and counted the number of clones (>50 cells).
Transwell Migration and Wound Healing Assays
BGC823 and MKN45 cells were transfected with GLP2R siRNA for 48 h and cultured in 24-well culture plates with 8 mm pore-containing membrane inserts to measure cell migration capacity. 4 × 104 cells were seeded on the upper transwell chambers in 200 µl serum-free culture medium, and 500 µl medium containing 20% FBS was added to the lower chambers. After 24 h incubation, the cells that migrated through membranes were fixed with methanol, stained with 1% crystal violet and counted under light microscope (200×). Additionally, BGC823 and MKN45 cells were cultured in 24-well-plates and scraped with a 10-μl pipette tip. The cells were cultured in RPMI 1640 medium without FBS. Images of wounds were captured at 0, 24, and 48 h, the area of wounds was quantified by ImageJ software (40×).
Statistical Analysis
All data of this study were statistically analyzed by R software 3.6.1 and Prism 9.0. The data were analyzed by two-tailed Student’s t-test and one-way analysis of variance (ANOVA). The difference was considered statistically significant when the p-value was less than 0.05.
Results
Primary Genetic Alterations and Statistical Analysis in STAD
We obtained the clinical information and whole-exome sequencing results of STAD patients from the TCGA database (Table 1). Maftools were applied to summarize the mutation data. We further categorized mutations based on the variant effect predictor, among which the frequency of missense mutations is highest (Figure 1A). Among all mutation types, the occurrence frequency of SNP is the highest (Figure 1B). Similarly, the most frequent type of SNV in STAD is C > T transversion (Figure 1C). Figure D shows the mutation of each sample, which is most relevant to TMB (Figure 1D). Figure E also depicts the mutations of the samples (Figure 1E). The first 10 mutated genes are TTN, MUC16, TP53, LRP1B, ARID1A, SYNE1, FAT4, CSMD3, FLG, and PCLO (Figure 1F). The waterfall diagram graphically shows the gene mutations of the samples (Figure 2). The correlation graph visualizes the correlation between two gene mutations; for example, mutations of PIK3CA and ARID1A co-occur, while mutations of PIK3CA and TP53 are mutually exclusive (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g001.jpg)
FIGURE 1. Primary genetic changes in STAD. (A,B) Gene variant categories in STAD. (C) The SNV classification of STAD. (D,E) The STAD samples’ mutations. (F) STAD mutation profile and first 10 mutant genes.
TMB was calculated by dividing the non-synonymous protein-coding variance by the genome’s overall sequence length. Based on the median TMB, we classified STAD patients into two groups: high TMB and low TMB (Supplementary File S1). The survival rate of patients in the high TMB group is higher than that of patients in the low TMB group, according to Kaplan-Meier survival analysis (Supplementary Figure S1A). The results show that in STAD, TMB is favorably connected with age and negatively correlated with stage, T, and N stages; female patients have higher TMB than male patients. TMB in female patients is higher than in male individuals (Supplementary Figures S1B–H). Between the high and low TMB groups, 816 DEGs were found (Supplementary File S2).
GSEA Enrichment and Immune Infiltration Analysis
The top 5 GO keywords (Figures 4A–E) and KEGG pathways (Figures 4F–J) with the highest significant enrichment in the high TMB group (p < 0.05) were obtained using GSEA enrichment analysis comparing the high and low TMB groups. These KEGG pathways and GO keywords are mostly related to genetic material metabolism and DNA repair.
According to studies, the higher the TMB in tumors, the more neoantigens are produced, making tumors more immunogenic (Rizvi et al., 2015). As a result, we investigated the link between TMB and immune markers in STAD. Using the CIBERSORT method, we determined the fraction of invading immune cells. The inference of immune cell types generated by CIBERSORT was relatively reliable at the p < 0.05 threshold. These findings demonstrated that high TMB tumors included large proportions of follicular helper T cells, activated memory CD4+ T-cells, M1 macrophages, M0 macrophages, and Neutrophils. Resting memory CD4+ T cells, regulatory T cells, monocytes, resting dendritic cells, and resting mast cells were all higher in the low TMB group (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g005.jpg)
FIGURE 5. TMB and immune markers in STAD were found to be correlated. In high (red violin part) and low (blue violin part) TMB groups, there are 22 types of adaptive and innate immune cells (by wilcox test).
Construction and Validations of Prognostic Model
From the overlap of immune-related genes and prior DEGs, differentially expressed immune genes were derived (Supplementary File S3). One-time random grouping was used to create the training (Supplementary File S4) and internal validation groups (Supplementary File S5). Based on univariate COX analysis (p < 0.01) for the training group (Supplementary File S6), we applied the LASSO algorithm to deal with the 21 survival-associated immune genes. Following the stepwise multivariate Cox regression analysis, four immune genes were able to develop a predictive model. The riskScore for each patient is calculated as follows: 0.001*APOD +0.005*APOH +0.039*INHA +0.499*GLP2R (Table 2). Each model gene is considered to with high risk attribute.
We split the patients in the training, internal validation, total, and external validation groups into high and low risk groups (Supplementary Files S7, S8, S9, S10) for later survival analysis after computing their riskScores. The survival curves between the high and low risk groups in the training, internal validation, total, and external validation groups are notably different, and the low risk group’s survival rate is significantly greater than the high risk group’s. The accuracy of the prognostic model we developed is dependable, according to the area under the curve (AUC) values of ROC curves (Figure 6). Ultimately, univariate and multivariate prognostic studies (p < 0.05) show that the riskScore derived from the model is a prognostic factor that is independent of other factors (Supplementary Figure S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g006.jpg)
FIGURE 6. Survival curves and ROC curves. p < 0.05 reveals significant survival differences. The ROC curves show the model prediction efficiency. (A) Training group. (B) Internal validation group. (C) Total group. (D) External validation group GSE84433.
Boxplots of model genes were created to compute the expression level differences between the high and low risk groups (Figures 7A–D). The high-risk group has significantly higher expression of all model genes, which validates the earlier inference that all model genes belong to high risk factors. Furthermore, incorporating all model genes in a multiple gene survival study successfully verifies the efficiency of the prognostic model (Figure 7E). A dynamic nomograph App online (https://u20131050.shinyapps.io/STAD-TMB-Dynamic_nomogram/) was successfully designed to quickly compute patient survival and support clinical decision making in order to improve the clinical translational implications of our study. We take low-risk group sample TCGA-VQ-AA6J and high-risk group sample TCGA-D7-A6F0 as examples to show the calculation results of dynamic nomograph App, which shows that the prediction results are relatively accurate (Supplementary Figure S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g007.jpg)
FIGURE 7. Multiple validations of model genes. (A–D) Boxplots display the differential expression for model genes. *, p < 0.05. **, p < 0.01. ***, p < 0.001. (E) Multiple gene survival validation, from the PROGgeneV2 online tool.
Alteration of Model Genes and Immunity
The cBioportal website provided an overall view of genetic change (Figure 8A) as well as domain mutation plots (Figures 8B–E) for the model genes. These model genes have extremely low mutation frequencies.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g008.jpg)
FIGURE 8. Mutation profiles of model genes. (A) The mutation summary of genetic alteration. Each color refers to corresponding alteration attribute. (B–E) Gene domain mutation. Each color refers to corresponding domain.
According to the TIMER database, CNV of these model genes has little effect on immune cell infiltration level differences (Figures 9A–D), and only Macrophage cell infiltration has a substantial impact on STAD patients’ survival (Figure 9E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g009.jpg)
FIGURE 9. CNV of the model genes and immune cells in STAD. (A–D) The effect of CNV of these model genes on the penetration level difference of immune cells. (E) The correlation between the recruitment of immune cells and the survival of STAD patients.
Visualization of Mutation Data and ssGSEA Analysis
We created a heatmap (Figure 10A) and violin plots (Figures 10B–E) of the tumor microenvironment by analyzing the correlation between risk categorization and tumor microenvironment. In conclusion, high-risk group had higher StromalScore, ImmuneScore, and ESTIMATEScore than low-risk group, while the low-risk group had a higher score of TumorPurity. To investigate the relationship between the risk score and immune status, we used the ssGSEA method to quantify 23 immune cell subsets, discovered that the infiltration of 13 immune cell subsets was significantly higher in the high-risk group than in the low-risk group (p < 0.05). (Figure 10F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g010.jpg)
FIGURE 10. An examination of the immune status of the tumor microenvironment. (A) heatmap. The sample name is represented by the abscissa, and the immune gene set is represented by the ordinate. The upper section contains tumor microenvironment scores, in which Stromal Score, Immune Score, and ESTIMATE Score decrease as risk low, implying that the content of corresponding cells decreases. Tumor purity is diametrically opposed to them. (B–E) Violin Plot. The statistical relationship between the degree of risk and each tumor microenvironmental parameter is investigated. (F) box plot. Immune cell subset differences in high- and low-risk groups. *, p < 0.05. **, p < 0.01. ***, p < 0.001. ns, p > 0.05.
Efficacy Prediction of Immunotherapy and Model Comparison
The high risk group had a higher TIDE score, lower MSI score, higher Dysfunction score, and higher Exclusion score, indicating that the high risk group of STAD has increased immune escape potential, resulting in poor immunotherapy efficacy (Figure 11A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g011.jpg)
FIGURE 11. Efficacy prediction of immunotherapy and model comparison. (A) Violin plots. The comparison of immunotherapy biomarkers in high and low risk categories. *, p < 0.05. **, p < 0.01. ***, p < 0.001. (B) 5-years ROC curves.
In the end, it can be judged from the 5-years ROC curves that the prognostic model constructed by us has the largest AUC value, so its prognostic prediction efficiency is higher than TIDE score and TIS score (Figure 11B). The above results show that the group with high expression of GLP2R has the highest risk score and the worst prognosis.
Drug Sensitivity Testing
By performing a separate drug sensitivity analysis on model genes in the prognostic model, we were able to identify the top 16 drugs with the most statistically significant differences. APOD expression was found to be positively related to the sensitivity of vemurafenib, pd-98059, dabrafenib, hypothemycin, selumetinib, bafetinib, denileukin diftitox ontak and cobimetinib (isomer 1), it indicates that the higher level of APOD expression, the greater sensitivity to the aforementioned drugs, but APOD has a negative correlation with pyrazoloacridine, batracylin, docetaxel and pralatrexate. APOH expression was revealed to be highly linked to the sensitivity of elesclomol, INHA expression was discovered to be strongly tied to the sensitivity of fulvestrant, but inversely correlated with the sensitivity of amonafide. Furthermore, the higher the expression of GLP2R in STAD patients, the greater the patient’s sensitivity to decitabine (Figure 12A). To further improve the clinical value of tumor mutation burden-related prognostic model for the treatment of stomach cancer. We analyzed the commonly used drugs in the clinical treatment of gastric cancer, which include cisplatin, doxorubicin, gemcitabine, lapatinib, and gemcitabine was noticed to be more sensitive in the high-risk group than in the low-risk group (p < 0.001), whereas lapatinib is the opposite (Figures 12B–E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g012.jpg)
FIGURE 12. Gene–drug sensitivity analysis. (A) based on the CellMiner database, the top 16 strongest correlations are shown. (B–E) Used the “pRRophetic” R package to compare differences in drug susceptibility between high and low risk groups, the drug sensitivity of (B) Cisplatin, (C) Doxorubicin, (D) Gemcitabine, (E) Lapatinib in high and low risk groups, respectively displayed.
Downregulation of GLP2R Inhibits STAD Cell Proliferation and Migration
To evaluate the specific role of GLP2R in gastric cancer, the relative mRNA expression level of GLP2R in the eight gastric cancer cell lines (GES1, BGC-823, MKN45, SNU-216, SGC-7901, MGC-803, AGS and N87) was analyzed. The expression level of GLP2R in the BGC-823, MKN45 cell lines was higher than that in the other cell lines (Figure 13A). We first evaluated the transfection efficiency of the cells by qRT-PCR and Western blot, found that the relative expression level of GLP2R was significantly lower after siRNA 3 transfection (Figures 13B,C). To further confirm the role of GLP2R in proliferation, We performed CCK-8 assays to detect the effect of GLP2R knockdown. After GLP2R silencing, BGC-823, MKN45 cell proliferation significantly decreased compared to control cells (Figure 13D). Colony formation assay also indicated that GLP2R silencing significantly suppressed the growth of BGC-823 cell (Supplementary Figures S4A–B). Transwell and wound healing assays were performed to detect migration, Our results showed that the migration rates of BGC-823, MKN45 cells transfected with siRNA were significantly lower than that of the control-transfected cells (Figures 13E–H). These data suggest that GLP2R knockdown repressed the proliferative and migratory abilities of BGC-823, MKN45 cells.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/790920/fcell-10-790920-HTML/image_m/fcell-10-790920-g013.jpg)
FIGURE 13. Knockdown of GLP2R inhibits proliferation, migration of gastric cancer cells. (A) The relative expression level of GLP2R in the GES1, BGC-823, MKN45, SNU-216, SGC-7901, MGC-803, AGS and N87 cell lines detected by RT-PCR. (B) The transfection efficiency of si-GLP2R in the BGC-823 and MKN45 cell lines detected by RT-PCR. (C) Western blot analysis confirmed that the expression of GLP2R was inhibited by GLP2R siRNA administration. (D) The CCK-8 assay was used to detect the effect of si-GLP2R on the proliferation of BGC-823 and MKN45 cell lines. (E) Representative images of the wound healing assay. (F) Statistical analysis of the wound healing assay results after Knockdown of GLP2R. (G) Representative images of the transwell assay. (H) Statistical analysis of the transwell assay results after Knockdown of GLP2R in the BGC-823 and MKN45 cell lines.*p < 0.05, **p < 0.01, and ***p < 0.001.
Discussion
Cancer is an inherited disease, and the mechanisms underlying involve the accumulation of harmful somatic mutations, which leads to a phenotypic consequence. Mutations are caused by faulty repair after DNA replication (especially meiosis) errors or DNA damage (such as exposure to radiation or carcinogens). Mutational processes contribute to different cancers, and the distribution of the rates of different mutational processes in various cancer types also varies (Martincorena and Campbell, 2015). High gene alterations make the tumor more immunogenic, making it a target for immune cell activation and susceptible to programmed cell death -1 (PD-1) immunocheckpoint inhibitors. TMB has been linked to a positive response to checkpoint inhibitor treatment (Snyder et al., 2014; Van Allen et al., 2015). In patients with advanced gastric cancer, a recent study found that pembrolizumab immunotherapy has good antitumor activity and low toxicity (Amatatsu et al., 2018). Due to the generation of immunogenic neoantigens, TMB, an emerging biomarker for response to immunotherapy, has also been a possible predictor of clinical benefits from immunocheckpoint inhibitors across several tumor types (Chalmers et al., 2017; Gong et al., 2019).
Higher TMB, as calculated from the complete exome, is related with greater immunotherapy responsiveness in patients with melanoma or non-small cell lung cancer, according to some studies (Chalmers et al., 2017). TMB levels beyond a certain threshold are more likely to create neoantigens in the tumor microenvironment, boosting antigenic presentation and enticing T cells to infiltrate tumors (Seto et al., 2019). A series of stromal cells, including vascular cells, fibroblasts, and inflammatory cells, come into being the tumor microenvironment (Greten and Grivennikov, 2019). More and more evidence suggests that the tumor microenvironment is important in the beginning and progression of STAD. Meanwhile, tumor-infiltrating lymphocytes have been shown to be effective in predicting GC patient prognosis, demonstrating that immune-related protective and hazardous variables in the GC microenvironment can become prognostic predictors (Yang et al., 2019).
TMB appears to have a considerable impact on the prognosis of STAD, in a study that included TGCA data from 375 gastric cancer patients, 15.7% of the patients had high TMB, and screening 632 up-regulated genes and 979 down-regulated genes were selected (Bai et al., 2020), further pathway analysis showed that patients in the high TMB group were associated with high activated CD4+ memory T cells, follicular helper T cells, quiescent NK cells, M0 and M1 macrophages, and neutrophils Cell infiltration.
TMB has been reported in many articles in tumor patients, as a novel immunotherapy marker, TMB levels are linked to immune cell infiltration in the tumor microenvironment, the higher the TMB, the better the outcome (Li et al., 2016; Kelly, 2017; Wei et al., 2021).
The clinical data of 63 patients with advanced stomach cancer who received immunotherapy were analyzed, assessed whether TMB is associated with response to immunotherapy, the results showed that high TMB was associated with the effectiveness of ICI treatment and chemotherapy, and no progression in the high TMB group extended survival, so TMB may serve as a predictive biomarker for patients with advanced gastric cancer treated with ICI material that aids clinical decision-making (Kim et al., 2020).
We acquired 816 DEGs by comparing the low and high TMB groups. We can identify a close association between genetic material metabolism, DNA repair, and TMB based on the GSEA enrichment data. TMB deficiency impairs genetic material metabolism and is linked to the DNA repair mechanism, an enzyme system that prevents changes in genetic material, as seen in non-small cell lung cancer (Chae et al., 2019). According to our findings, those with a high TMB had more activated memory CD4+ T-cells, follicular helper T cells, M0 macrophages, M1 macrophages, and Neutrophils. High TMB encourages the invasion of these immune cells, as evidenced by studies in different malignancies (Zhang et al., 2019; Wu et al., 2020; Xu et al., 2020). These findings suggested that TMB can alter the features of immune cell infiltration, and that TMB levels beyond a certain threshold can attract immune effector cells. The relative balance between immunosuppressive and anti-tumor immune cells is one of the methods by which tumor cells retain immune-mediated dormancy (Mittal et al., 2014). Activated memory CD4+ T-cells and follicular helper T cells are archetypal anti-tumor immune cells, as we all know. Follicular helper T cells are a subpopulation of CD4 T cells that support B cells in the germinal center of lymphatic tissue. Classically activated M1 macrophages may be implicated in the early removal stage of immunoediting driven by CD8 cytotoxic T lymphocytes and interferons, where they destroy tumor cells and cause tissue death (Quaranta and Schmid, 2019). TMB and the immune microenvironment have a very close interaction, as can be observed. In the growth of malignancies, the interaction of these immune cells is critical.
According to the results of the ssGSEA analysis, the degree of TMB in STAD patients is positively correlated with the level of immune infiltration. According to GSEA enrichment analysis, the risk was positively correlated with StromalScore, ImmuneScore, and ESTIMATEScore, but negatively coupled to TumorPurity, and the higher risk, the stronger immunosuppressive activity. This finding could explain why STAD patients with higher risk levels have a lower survival rate.
In our study, which also indicates that the risk score of the model composed of APOD, APOH, INHA, and GLP2R. At the same time, we demonstrates that increased GLP2R expression is attached to an increase in cancer cell sensitivity to decitabine. Decitabine is a chemotherapeutic pyrimidine nucleoside analogue used for the treatment of myelodysplastic syndromes by inducing DNA hypomethylation and corresponding alterations in gene expression (https://go.drugbank.com/drugs/DB01262). Therefore, we guess that the model gene can be used as an evaluation index of drug efficacy.
It's almost probable that the model genes we discovered are high-risk genes, and that the prognostic model we built is accurate. Through bioinformatics analysis, we found that GLP2R has the highest risk coefficient. And GLP2R was reported to be associated with gastrointestinal cancer in previous studies (Lu et al., 2019; Qiu et al., 2020). Therefore, in order to investigate whether it can be used as a prognostic indicator, we performed in gastric cancer cell lines (BGC823, MKN45) and found that knockdown of GLP2R can significantly inhibit the proliferation and migration of gastric cancer cell lines.
Although there is little study on the significance of these model genes in STAD, they are essential in other malignancies, therefore further mechanism studies are needed. For example, APOD can be employed as a therapeutic tool to promote tumor cell death as a prognostic marker for numerous cancer types (Søiland et al., 2009; Bajo-Grañeras et al., 2013). Similarly, INHA induces tumor angiogenesis and promotes tumor metastasis as a poor survival predictor of various cancer types, and is a potential target for anti-angiogenesis therapy and genetic engineering therapy (Singh et al., 2018; Yoon et al., 2018).
The results of mutation study show that our model genes have low mutation frequency. Perhaps, on the one hand, the model genes in our study do not rely on structural and functional alterations to influence STAD prognosis; on the other hand, they generate high risk effect in STAD as a result of quantitative changes. The judgment could also be explained by the fact that the CNV of these model genes has no effect on the difference in immune cell infiltration levels. Only Macrophage cells, out of the six types of immune cells, can alter STAD patients’ survival. Combined with previous violin plot of immune cell infiltration, to a large extent, Macrophages M0, M1 rather than M2, may play a vitalrole in the prognosis of STAD. At present, studies on the role of Macrophage cells in STAD are still lacking. Therefore, this is a direction worth exploring.
Conclusion
In conclusion, our results imply that STAD patients with a high TMB have a better prognosis. We identified DEGs between high and low TMB groups using STAD samples from the TCGA database, and investigated the relationship between immune cell infiltration features and TMB. In addition, the immune prognostic model was built using a range of bioinformatics methodologies, with several validations. Our study’s clinical translational importance was further strengthened by the establishment of the dynamic nomograph App online and immunotherapy prediction based on the prognostic model. Finally, GLP2R may be expected to be a potential target for gastric cancer.
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
GH and QL conceived of and designed the study. YZ, MF, XP, XL, and NL performed the literature search, generated the figures and tables, and wrote the manuscript. WZ, FY, ZC, SM, and YH supervised the study and reviewed the manuscript. All authors read and approved the final manuscript.
Funding
This study was funded by the National Natural Sciences Foundation of China (Grant No. 82003312).
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 would like to thank the TCGA and GEO databases for the availability of the data. We thank the technical support by the Huazhong University of Science & Technology Analyttical & Testing center.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.790920/full#supplementary-material
Abbreviations
AUC, Area under the curve; CNV, Copy number variation; DEGs, Differentially expressed genes; DNA, Deoxyribonucleic acid; GEO, Gene Expression Omnibus; GC, Gastric cancer; GSEA, Gene Set Enrichment Analysis; GO, Gene Ontology; ICI, Immune checkpoint inhibitor; EMT, Epithelial-mesenchymal transition; FDR, False discovery rate; KEGG, Kyoto Encyclopedia of Genes and Genomes; LASSO, Least absolute shrinkage and selection operator; MSI, Microsatellite instability; PD-1, Programmed cell death -1; ROC, Receiver Operating Characteristic; STAD, Stomach adenocarcinoma; TCGA, The Cancer Genome Atlas; TIMER, Tumor Immune Estimation Resource; TIDE, Tumor immune dysfunction and exclusion; TIS, Tumor inflammation signature; TMB, Tumor mutation burden.
References
Amatatsu, M., Arigami, T., Uenosono, Y., Yanagita, S., Uchikado, Y., Kijima, Y., et al. (2018). Programmed Death-Ligand 1 Is a Promising Blood Marker for Predicting Tumor Progression and Prognosis in Patients with Gastric Cancer. Cancer Sci. 109 (3), 814–820. doi:10.1111/cas.13508
Bai, J., Yang, B., Shi, R., Shao, X., Yang, Y., Wang, F., et al. (2020). Could Microtubule Inhibitors Be the Best Choice of Therapy in Gastric Cancer with High Immune Activity: Mutant DYNC1H1 as a Biomarker. Aging 12 (24), 25101–25119. doi:10.18632/aging.104084
Bajo-Grañeras, R., Crespo-Sanjuan, J., García-Centeno, R. M., Garrote-Adrados, J. A., Gutierrez, G., García-Tejeiro, M., et al. (2013). Expression and Potential Role of Apolipoprotein D on the Death-Survival Balance of Human Colorectal Cancer Cells under Oxidative Stress Conditions. Int. J. Colorectal Dis. 28 (6), 751–766. doi:10.1007/s00384-012-1616-2
Barrett, T., Troup, D. B., Wilhite, S. E., Ledoux, P., Rudnev, D., Evangelista, C., et al. (2007). NCBI GEO: Mining Tens of Millions of Expression Profiles-Ddatabase and Tools Update. Nucleic Acids Res. 35 (Database issue), D760–D765. doi:10.1093/nar/gkl887
Bhattacharya, S., Andorf, S., Gomes, L., Dunn, P., Schaefer, H., Pontius, J., et al. (2014). ImmPort: Disseminating Data to the Public for the Future of Immunology. Immunol. Res. 58 (2-3), 234–239. doi:10.1007/s12026-014-8516-1
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: A Cancer J. Clinicians 68 (6), 394–424. doi:10.3322/caac.21492
Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data: Figure 1. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.cd-12-0095
Chae, Y. K., Davis, A. A., Raparia, K., Agte, S., Pan, A., Mohindra, N., et al. (2019). Association of Tumor Mutational Burden with DNA Repair Mutations and Response to Anti-PD-1/pd-L1 Therapy in Non-small-cell Lung Cancer. Clin. Lung Cancer 20 (2), 88–96. e86. doi:10.1016/j.cllc.2018.09.008
Chalmers, Z. R., Connelly, C. F., Fabrizio, D., Gay, L., Ali, S. M., Ennis, R., et al. (2017). Analysis of 100,000 Human Cancer Genomes Reveals the Landscape of Tumor Mutational burden. Genome Med. 9 (1), 34. doi:10.1186/s13073-017-0424-2
Chen, B., Khodadoust, M. S., Liu, C. L., Newman, A. M., and Alizadeh, A. A. (2018). Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol. Biol. 1711, 243–259. doi:10.1007/978-1-4939-7493-1_12
Coutzac, C., Pernot, S., Chaput, N., and Zaanan, A. (2019). Immunotherapy in Advanced Gastric Cancer, Is it the Future? Crit. Rev. oncology/hematology 133, 25–32. doi:10.1016/j.critrevonc.2018.10.007
Cristescu, R., Lee, J., Nebozhyn, M., Kim, K.-M., Ting, J. C., Wong, S. S., et al. (2015). Molecular Analysis of Gastric Cancer Identifies Subtypes Associated with Distinct Clinical Outcomes. Nat. Med. 21 (5), 449–456. doi:10.1038/nm.3850
Danaher, P., Warren, S., Lu, R., Samayoa, J., Sullivan, A., Pekker, I., et al. (2018). Pan-cancer Adaptive Immune Resistance as Defined by the Tumor Inflammation Signature (TIS): Results from the Cancer Genome Atlas (TCGA). J. Immunotherapy Cancer 6 (1), 63. doi:10.1186/s40425-018-0367-1
Gong, J., Chehrazi-Raffle, A., Placencio-Hickok, V., Guan, M., Hendifar, A., and Salgia, R. (2019). The Gut Microbiome and Response to Immune Checkpoint Inhibitors: Preclinical and Clinical Strategies. Clin. Transl Med. 8 (1), 9. doi:10.1186/s40169-019-0225-x
Goswami, C. P., and Nakshatri, H. (2014). PROGgeneV2: Enhancements on the Existing Database. BMC Cancer 14 (970), 970–2407. doi:10.1186/1471-2407-14-970
Greten, F. R., and Grivennikov, S. I. (2019). Inflammation and Cancer: Triggers, Mechanisms, and Consequences. Immunity 51 (1), 27–41. doi:10.1016/j.immuni.2019.06.025
He, Y., Liu, Z.-X., Jiang, Z.-H., and Wang, X.-S. (2019). Identification of Genomic Features Associated with Immunotherapy Response in Gastrointestinal Cancers. Wjgo 11 (4), 270–280. doi:10.4251/wjgo.v11.i4.270
Hutter, C., and Zenklusen, J. C. (2018). The Cancer Genome Atlas: Creating Lasting Value beyond its Data. Cell 173 (2), 283–285. doi:10.1016/j.cell.2018.03.042
Jiang, P., Gu, S., Pan, D., Fu, J., Sahu, A., Hu, X., et al. (2018). Signatures of T Cell Dysfunction and Exclusion Predict Cancer Immunotherapy Response. Nat. Med. 24 (10), 1550–1558. doi:10.1038/s41591-018-0136-1
Kelly, R. J. (2017). Immunotherapy for Esophageal and Gastric Cancer. Am. Soc. Clin. Oncol. Educ. Book 37, 292–300. doi:10.1200/edbk_175231
Kim, H. S., and Seo, H. K. (2018). Immune Checkpoint Inhibitors for Urothelial Carcinoma. Investig. Clin. Urol. 59 (5), 285–296. doi:10.4111/icu.2018.59.5.285
Kim, J., Kim, B., Kang, S. Y., Heo, Y. J., Park, S. H., Kim, S. T., et al. (2020). Tumor Mutational Burden Determined by Panel Sequencing Predicts Survival after Immunotherapy in Patients with Advanced Gastric Cancer. Front. Oncol. 10, 314. doi:10.3389/fonc.2020.00314
Lazăr, D. C., Avram, M. F., Romoșan, I., Cornianu, M., Tăban, S., and Goldiș, A. (2018). Prognostic Significance of Tumor Immune Microenvironment and Immunotherapy: Novel Insights and Future Perspectives in Gastric Cancer. World J. Gastroenterol. 24 (32), 3583–3616. doi:10.3748/wjg.v24.i32.3583
Lee, I.-S., Ahn, J.-Y., Yook, J.-H., and Kim, B.-S. (2017). Mediastinal Lymph Node Dissection and Distal Esophagectomy Is Not Essential in Early Esophagogastric junction Adenocarcinoma. World J. Surg. Onc 15 (1), 28. doi:10.1186/s12957-016-1088-x
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.can-17-0307
Li, X., Wu, W. K. K., Xing, R., Wong, S. H., Liu, Y., Fang, X., et al. (2016). Distinct Subtypes of Gastric Cancer Defined by Molecular Characterization Include Novel Mutational Signatures with Prognostic Capability. Cancer Res. 76 (7), 1724–1732. doi:10.1158/0008-5472.can-15-2443
Li, Y., Wang, H. C., Wang, J. S., Sun, B., and Li, L. P. (2020). Chemokine Receptor 4 Expression Is Correlated with the Occurrence and Prognosis of Gastric Cancer. FEBS Open Bio 10 (6), 1149–1161. doi:10.1002/2211-5463.12864
Liang, C., Wu, H.-M., Yu, W.-M., and Chen, W. (2021). Research Status on Immunotherapy Trials of Gastric Cancer. Wjcc 9 (21), 5782–5793. doi:10.12998/wjcc.v9.i21.5782
Lu, Y., Kweon, S.-S., Tanikawa, C., Jia, W.-H., Xiang, Y.-B., Cai, Q., et al. (2019). Large-Scale Genome-wide Association Study of East Asians Identifies Loci Associated with Risk for Colorectal Cancer. Gastroenterology 156 (5), 1455–1466. doi:10.1053/j.gastro.2018.11.066
Martincorena, I., and Campbell, P. J. (2015). Somatic Mutation in Cancer and normal Cells. Science 349, 1483–1489. doi:10.1126/science.aab4082
Matsushita, H., Vesely, M. D., Koboldt, D. C., Rickert, C. G., Uppaluri, R., Magrini, V. J., et al. (2012). Cancer Exome Analysis Reveals a T-cell-dependent Mechanism of Cancer Immunoediting. Nature 482 (7385), 400–404. doi:10.1038/nature10755
Mayakonda, A., Lin, D.-C., Assenov, Y., Plass, C., and Koeffler, H. P. (2018). Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res. 28 (11), 1747–1756. doi:10.1101/gr.239244.118
Mehmedagic, I., Hasukic, S., Agic, M., Kadric, N., and Hasukic, I. (2016). Influence of Prognostic Factors for Recurrence of Adenocarcinoma of the Stomach. Med. Arh 70 (6), 441–444. doi:10.5455/medarh.2016.70.441-444
Mittal, D., Gubin, M. M., Schreiber, R. D., and Smyth, M. J. (2014). New Insights into Cancer Immunoediting and its Three Component Phases-Elimination, Equilibrium and Escape. Curr. Opin. Immunol. 27, 16–25. doi:10.1016/j.coi.2014.01.004
Naumann, M. (2005). Pathogenicity Island-dependent Effects of on Intracellular Signal Transduction in Epithelial Cells. Int. J. Med. Microbiol. 295 (5), 335–341. doi:10.1016/j.ijmm.2005.06.007
Powles, T., and Morrison, L. (2018). Biomarker Challenges for Immune Checkpoint Inhibitors in Urothelial Carcinoma. Nat. Rev. Urol. 15 (10), 585–587. doi:10.1038/s41585-018-0056-3
Qiu, X.-T., Song, Y.-C., Liu, J., Wang, Z.-M., Niu, X., and He, J. (2020). Identification of an Immune-Related Gene-Based Signature to Predict Prognosis of Patients with Gastric Cancer. Wjgo 12 (8), 857–876. doi:10.4251/wjgo.v12.i8.857
Quaranta, V., and Schmid, M. C. (2019). Macrophage-Mediated Subversion of Anti-tumour Immunity. Cells 8 (7). 747, doi:10.3390/cells8070747
Riaz, N., Morris, L., Havel, J. J., Makarov, V., Desrichard, A., and Chan, T. A. (2016). The Role of Neoantigens in Response to Immune Checkpoint Blockade. Intimm 28 (8), 411–419. doi:10.1093/intimm/dxw019
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Rizvi, N. A., Hellmann, M. D., Snyder, A., Kvistborg, P., Makarov, V., Havel, J. J., et al. (2015). Mutational Landscape Determines Sensitivity to PD-1 Blockade in Non-small Cell Lung Cancer. Science 348 (6230), 124–128. doi:10.1126/science.aaa1348
Seto, T., Sam, D., and Pan, M. (2019). Mechanisms of Primary and Secondary Resistance to Immune Checkpoint Inhibitors in Cancer. Med. Sci. (Basel) 7 (2), 14. doi:10.3390/medsci7020014
Singh, P., Jenkins, L. M., Horst, B., Alers, V., Pradhan, S., Kaur, P., et al. (2018). Inhibin Is a Novel Paracrine Factor for Tumor Angiogenesis and Metastasis. Cancer Res. 78 (11), 2978–2989. doi:10.1158/0008-5472.can-17-2316
Snyder, A., and Chan, T. A. (2015). Immunogenic Peptide Discovery in Cancer Genomes. Curr. Opin. Genet. Dev. 30, 7–16. doi:10.1016/j.gde.2014.12.003
Snyder, A., Makarov, V., Merghoub, T., Yuan, J., Zaretsky, J. M., Desrichard, A., et al. (2014). Genetic Basis for Clinical Response to CTLA-4 Blockade in Melanoma. N. Engl. J. Med. 371 (23), 2189–2199. doi:10.1056/nejmoa1406498
Søiland, H., Janssen, E. A., Kørner, H., Varhaug, J. E., Skaland, I., Gudlaugsson, E., et al. (2009). Apolipoprotein D Predicts Adverse Outcome in Women >or=70 Years with Operable Breast Cancer. Breast Cancer Res. Treat. 113 (3), 519–528. doi:10.1007/s10549-008-9955-y
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Tibshirani, R. (1997). The Lasso Method for Variable Selection in the Cox Model. Statist. Med. 16 (4), 385–395. doi:10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3
Van Allen, E. M., Miao, D., Schilling, B., Shukla, S. A., Blank, C., Zimmer, L., et al. (2015). Genomic Correlates of Response to CTLA-4 Blockade in Metastatic Melanoma. Science 350, 207–211. doi:10.1126/science.aad0095
Wang, F., Wei, X. L., Wang, F. H., Xu, N., Shen, L., Dai, G. H., et al. (2019). Safety, Efficacy and Tumor Mutational burden as a Biomarker of Overall Survival Benefit in Chemo-Refractory Gastric Cancer Treated with Toripalimab, a PD-1 Antibody in Phase Ib/II Clinical Trial NCT02915432. Ann. Oncol. 30 (9), 1479–1486. doi:10.1093/annonc/mdz197
Wei, X. L., Xu, J. Y., Wang, D. S., Chen, D. L., Ren, C., Li, J. N., et al. (2021). Baseline Lesion Number as an Efficacy Predictive and Independent Prognostic Factor and its Joint Utility with TMB for PD-1 Inhibitor Treatment in Advanced Gastric Cancer. Ther. Adv. Med. Oncol. 13, 1758835921988996. doi:10.1177/1758835921988996
Wu, Z., Wang, M., Liu, Q., Liu, Y., Zhu, K., Chen, L., et al. (2020). Identification of Gene Expression Profiles and Immune Cell Infiltration Signatures between Low and High Tumor Mutation burden Groups in Bladder Cancer. Int. J. Med. Sci. 17 (1), 89–96. doi:10.7150/ijms.39056
Xu, F., Lin, H., He, P., He, L., Chen, J., Lin, L., et al. (2020). A TP53-Associated Gene Signature for Prediction of Prognosis and Therapeutic Responses in Lung Squamous Cell Carcinoma. Oncoimmunology 9 (1), 1731943. doi:10.1080/2162402x.2020.1731943
Yang, B., Bai, J., Shi, R., Shao, X., Yang, Y., Jin, Y., et al. (2020). TGFB2 Serves as a Link between Epithelial-Mesenchymal Transition and Tumor Mutation burden in Gastric Cancer. Int. Immunopharmacology 84, 106532. doi:10.1016/j.intimp.2020.106532
Yang, W., Lai, Z., Lai, Z., Li, Y., Mu, J., Yang, M., et al. (2019). Immune Signature Profiling Identified Prognostic Factors for Gastric Cancer. Chin. J. Cancer Res. 31 (3), 463–470. doi:10.21147/j.issn.1000-9604.2019.03.08
Yoon, W., Yoo, Y., Chae, Y. S., Kee, S.-H., and Kim, B. M. (2018). Therapeutic Advantage of Genetically Engineered Salmonella typhimurium Carrying Short Hairpin RNA against Inhibin Alpha Subunit in Cancer Treatment. Ann. Oncol. 29 (9), 2010–2017. doi:10.1093/annonc/mdy240
Zhang, C., Li, Z., Qi, F., Hu, X., and Luo, J. (2019). Exploration of the Relationships between Tumor Mutation burden with Immune Infiltrates in clear Cell Renal Cell Carcinoma. Ann. Transl. Med. 7 (22), 648. doi:10.21037/atm.2019.10.84
Keywords: bioinformatics, mutation burden, stomach adenocarcinoma, survival, multi-omics, immunity
Citation: Fu M, Huang Y, Peng X, Li X, Luo N, Zhu W, Yang F, Chen Z, Ma S, Zhang Y, Li Q and Hu G (2022) Development of Tumor Mutation Burden-Related Prognostic Model and Novel Biomarker Identification in Stomach Adenocarcinoma. Front. Cell Dev. Biol. 10:790920. doi: 10.3389/fcell.2022.790920
Received: 07 October 2021; Accepted: 21 February 2022;
Published: 23 March 2022.
Edited by:
Weiguo Dong, Renmin Hospital of Wuhan University, ChinaReviewed by:
Liu Le Ping, Central South University, ChinaJinhui Liu, Nanjing Medical University, China
Copyright © 2022 Fu, Huang, Peng, Li, Luo, Zhu, Yang, Chen, Ma, Zhang, Li and Hu. 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: Yuanyuan Zhang, z1731224497@163.com; Qianxia Li, liqianx110@163.com; Guangyuan Hu, h.g.y.121@163.com
†These authors contributed equally to this work