Skip to main content
Download PDF

N6-methyladenosine-related lncRNAs identified as potential biomarkers for predicting the overall survival of Asian gastric cancer patients

BMC Cancer volume 22, Article number: 721 (2022) Cite this article

Abstract

Objective

Gastric cancer (GC) is one of the most prevalent malignant tumors in Asian countries. Studies have proposed that lncRNAs can be used as diagnostic and prognostic indicators of GC due to the high specificity of lncRNAs expression involvement in GC. Recently, N6-methyladenosine (m6A) has also emerged as an important modulator of the expression of lncRNAs in GC. This study aimed at establishing a novel m6A-related lncRNAs prognostic signature that can be used to construct accurate models for predicting the prognosis of GC in the Asian population.

Methods

First, the levels of m6A modification and m6A methyltransferases expression in GC samples were determined using dot blot and western blot analyses. Next, we evaluated the lncRNAs expression profiles and the corresponding clinical data of 88 Asian GC patients retrieved from The Cancer Genome Atlas (TCGA) database. Differential expression of m6A-related lncRNAs between GC and normal tissues was investigated. The relationship between these target lncRNAs and potential immunotherapeutic signatures was also analyzed. Gene set enrichment analysis (GSEA) was performed to identify the malignancy-associated pathways. Univariate Cox regression, LASSO regression, and multivariate Cox regression analyses were performed to establish a novel prognostic m6A-related lncRNAs prognostic signature. Moreover, we constructed a predictive nomogram and determined the expression levels of nine m6A-related lncRNAs in 12 pairs of clinical samples.

Results

We found that m6A methylation levels were significantly increased in GC tumor samples compared to adjacent normal tissues, and the increase was positively correlated with tumor stage. Patients were then divided into two clusters (cluster 1 and cluster 2) based on the differential expression of the m6A-related lncRNAs. Results showed that there was a significant difference in survival probability between the two clusters (p = 0.018). Notably, the low survival rate in cluster 2 may be associated with high expression of immune cells (resting memory CD4+ T cells, p = 0.027; regulatory T cells, p = 0.0018; monocytes, p = 0.00095; and resting dendritic cells, p = 0.015), and low expression of immune cells (resting NK cells, p = 0.033; and macrophages M1, p = 0.045). Enrichment analysis indicated that malignancy-associated biological processes were more common in the cluster 2 subgroup. Finally, the risk model comprising of six m6A-related lncRNAs was identified as an independent predictor of prognoses, which could divide patients into high- or low-risk groups. Time-dependent ROC analysis suggested that the risk score could accurately predict the prognosis of GC patients. Patients in the high-risk group had worse outcomes compared to patients in the low-risk group, and the risk score showed a positive correlation with immune cells (resting memory CD4+ T cells, R = 0.31, P = 0.038; regulatory T cells, R = 0.42, P = 0.0042; monocytes, R = 0.42, P = 0.0043). However, M1 macrophages (R = -0.37, P = 0.012) and resting NK cells (R = -0.31, P = 0.043) had a negative correlation with risk scores. Furthermore, analysis of clinical samples validated the weak positive correlation between the risk score and tumor stage.

Conclusions

The risk model described here, based on the six m6A-related lncRNAs signature, and may predict the clinical prognoses and immunotherapeutic response in Asian GC patients.

Peer Review reports

Introduction

According to worldwide cancer statistics, gastric cancer (GC) is the fifth most common cancer and the third cause of cancer-related death in developing countries, with a very high incidence in Asian countries, especially in China [1]. A recent report indicated that the number of new GC cases reached one million in 2020, of which 769,000 patients succumbed [2]. Specifically, about 150–200,000 people die of GC in China every year according to the 2015 China Cancer Data Report [3]. Notably, there are different factors associated with the occurrence and progression of gastric cancer, including gene polymorphism, lifestyle such as drinking, smoking, high dietary salt intake, and iron deficiency [4], and Helicobacter pylori (H. pylori) infection, which is the strongest known risk factor for GC [5]. Evidence suggests that early detection and accurate non-invasive diagnosis of gastric cancer increases the chances of successful treatment [6], with a high 5-year overall survival (OS) rate (> 60%) in patients with early gastric cancer [7]. For these patients, surgical resection is the best treatment option. In contrast, patients with advanced stages of GC progression require chemotherapy, which is the most important treatment in metastasis [8]. Given the rapid technological advancement, immunotherapy and targeted therapy have been shown to improve the OS of GC patients [9]. Therefore, this calls for identification of the risk factors associated with GC patients’ outcomes, which will aid in classifying GC patients into different risk groups for individualized treatment, with the overarching goal of improving their treatment outcomes.

There are several biomarkers commonly used in GC, including carcinoembryonic antigen (CEA), alpha-fetoprotein (AFP), and carbohydrate antigen (CA) [10]. Considering the low sensitivity and positivity rate of the available biomarkers, many researchers have conducted studies to establish a new, stable, and powerful predictive model for early GC screening. Long non-coding RNAs (lncRNAs) in body fluids have attracted widespread attention as prognostic biomarkers due to their high specificity and sensitivity [11, 12]. Circulating lncRNAs have been associated with many prognosis-related factors (tumor size, stage, depth of invasion, and lymph node metastasis) and have been demonstrated as prognostic biomarkers of gastric cancer [13, 14]. Previous studies have revealed that several circulating exosomal lncRNAs, such as lncRNA PCSK2–2:1 and GNAQ-6:1, are better biomarkers for distinguishing gastric cancer patients from healthy people compared to traditional diagnostic biomarkers, such as CEA and AFP [15, 16]. A recent study explored the landscape of lncRNAs associated with the immunome across 33 cancer types, and found that their expression level was correlated with the level of immune factors and immune cell infiltration [17]. This suggests that the development of an lncRNAs-based prediction model is not only beneficial for the prediction of OS in patients with early gastric cancer, but also for the selection of immunotherapy as the optimal treatment method for patients.

N6-methyladenosine (m6A), one of the most abundant internal modification in RNA, modulates several bioprocesses such as cell differentiation, DNA damage response, and carcinogenesis [18]. Studies have shown that m6A modification is important in the treatment of tumors, especially in targeted therapy [19]. The m6A modification is mediated by an m6A methyltransferase complex, composed of METTL3, METTL14, and WTAP, and is abolished by m6A demethylases, including ALKBH5 and FTO. At the same time, m6A-modification relies on reader proteins such as YTHDC1, YTHDC2, and YTHDF1 [20]. It is well known that, for tumor biology, the level of m6A modification changes significantly, and exerts its role in tumorigenesis, proliferation, metastasis, and immunoregulation, thereby influencing patients’ outcomes [21, 22]. Studies have reported that the level of m6A is significantly higher in hepatocellular carcinoma and gastric cancer patients than in normal individuals, and the level is positively correlated with the clinical outcome [23, 24]. However, in bladder cancer, a lower level of m6A is associated with a worse outcome [21]. This can be attributed to the fact that m6A modification has been identified on mRNA and some non-coding RNA, including mRNA, microRNA, lncRNA, snoRNA, and circRNA [25]. It is worth noting that m6A modification on RNA regulates the progress of oncogenesis and exerts critical roles in innate immunity [26], whereas m6A modification on lncRNA affects RNA-protein interaction [27]. For example, lncRNA MALAT1, associated with lung adenocarcinoma, was highly methylated with m6A [28]. Collectively, these studies suggest that m6A modification plays a critical role in regulating the expression of lncRNAs and may affect cancer development.

This study identified lncRNAs associated with m6A regulators, and analyzed their expression and related pathways in Asian gastric cancer patients. The six identified lncRNAs were then used to construct a prognosis model which could accurately predict the prognosis of patients, and may be potential biomarkers for classifying GC patients for targeted immunotherapy.

Methods

Clinical specimens

This study included GC patients admitted at the Tongling People’s Hospital from January 2018 to December 2019. Tissue specimens were collected from patients and immediately stored in liquid nitrogen, including 12 paired GC tumors and adjacent non-tumor tissues. Notably, signed informed consent was obtained from all study participants before enrollment. All participants had not received chemotherapy or radiotherapy before surgery. The study was approved by the Institutional Review Board of Nanjing Medical University and the Ethical committee of the Tongling People’s Hospital (ethical review No.2019–008) and was conducted in accordance with the Declaration of Helsinki.

m6A dot blot analysis

Total RNA was separately extracted from human tissues using TRIzol (Invitrogen) reagent and the RNA concentration was determined using a NanoDrop 2000 (Thermo Fisher Scientific). The RNA concentration was then adjusted to 250, 500, and 1000 ng/μL. Two microliters of Poly(A) + RNAs were first denatured by heating at 65 °C for 5 min and spotted onto a nitrocellulose membrane (Amersham, GE Healthcare, USA). Next, the membrane was cross-linked by 1500 J UV, blocked with 5% non-fat milk twice, and incubated with m6A antibody (1:1000, Thermo Fisher Scientific, USA) overnight at 4 °C. Subsequently, the membrane was incubated with HRP-conjugated goat anti-mouse IgG dilution (1:3000, Proteintech, USA). Finally, the membrane was visualized using the chemiluminescence imaging analysis system (Tanon, Shanghai, China). The other membrane was dyed in 0.02% methylene blue (MB) for 30 min, followed by rinsing twice with ultrapure water and photographing. The membrane was visualized with chemiluminescence image analysis system (Tanon, Shanghai, China) and the grayscale was measured by ImageJ (NIH, Bethesda, MD) software.

Western blot analysis

Western blotting experiments were performed in accordance with our previous study [29]. Briefly, tissues were harvested and lysed in RIPA buffer on ice, and the extracted protein was quantified by bicinchoninic acid (BCA) analysis (Beyotime, China). Next, proteins were resolved using 10% SDS-PAGE and transferred to polyvinylidene difluoride (PVDF) membranes (Millipore, USA). Nonspecific binding sites were then blocked by immersing the membranes in 5% bovine serum albumin in PBS at room temperature. Membranes were then incubated with the following high affinity primary antibodies: anti-WTAP antibody (1:1000, Cell Signaling Technology, USA), anti-METTL3 antibody (1:1000, Cell Signaling Technology, USA), anti-METTL14 antibody (1:1000, Cell Signaling Technology, USA), and anti-GAPDH antibody (1:1000, Cell Signaling Technology, USA). After washing, membranes were incubated with peroxidase (HRP)-conjugated secondary antibody (1:1000, Cell Signaling Technology, USA). Finally, protein signals were detected using a chemiluminescence imaging analysis system (Tanon, Shanghai, China).

Acquisition of information on Asian GC patients

RNA sequence transcriptome data 88 Asian GC patients and their clinicopathological information were downloaded from the Cancer Genome Atlas (TCGA) data portal (https://cancergenome.nih.gov/). In addition, expression matrixes of 23 m6A-related genes were extracted from the TCGA database, including expression data on writers (METTL3, METTL14, METTL16, WTAP, VIRMA, ZC3H13, RBM15, and RBM15B), erasers (FTO and ALKBH5), and readers (YTHDC1, YTHDC2, YTHDF1, YTHDF2, YTHDF3, HNRNPC, FMR1, LRPPRC, HNRNPA2B1, IGF2BP1, IGF2BP2, IGF2BP3, and RBMX). The clinicopathological information for these samples is shown in Table 1 and supplementary Table S1.

Table 1 Clinical pathological parameters of Asian patients with gastric cancer in TCGA database
Full size table

Prognostic signature construction

The long non-coding RNA annotation file of Genome Reference Consortium Human Build 38 (GRCh38) was acquired from the GENCODE website for annotation of the lncRNAs in the TCGA dataset. Pearson’s correlation analysis (with the |Pearson R| > 0.4 and p < 0.001) was used to screen m6A-related lncRNAs, which resulted in identification of 1054 m6A-related lncRNAs (supplementary Table S1). We used univariate Cox regression analysis combined with GC overall survival information in TCGA to explore the prognosis of m6A-related lncRNAs. In addition, limma and corrplot packages in R software were used to analyze the relationship between the differential expressions of lncRNAs, cancer-related factors, and immune factors in all samples. Gene Set Enrichment Analysis (GSEA) was then performed according to our previous study [30].

LncRNA expression profiles were normalized by [log2 (data + 1)]. Notably, the entire set included 72 patients since 16 patients with missing clinical data were excluded. The entire set was then randomly separated into the training set and the test set in a ratio of 5 to 5. Next, LASSO Cox regression and multivariate Cox proportional hazards regression models were performed to identify and analyze the prognostic genes in order to establish a predictive model. Notably, an m6A-related lncRNAs risk model was ultimately established. The following formula was used to calculate the risk score: Risk score = expression of lncRNA1 × β1lncRNA1 + expression of lncRNA2 × β2lncRNA2 + …expression of lncRNAn × βnlncRNAn (β represents the regression coefficient of each lncRNA). Two subgroups (low- and high- risk groups) were then established based on the median risk score. Subsequently, Kaplan-Meier curves and log-rank methods were used to evaluate the prognostic significance of the risk score. The receiver operating characteristic (ROC) curve analysis was used to evaluate the prognosis-prediction accuracy and risk score of each lncRNA, whereas the area under the curve (AUC) with cutoffs was used to measure prognosis prediction accuracy. Moreover, survival curves of clinic-pathologic characteristics and model validation between the two groups were generated using the Kaplan-Meier method. CIBERSORT algorithm was used to calculate the fractions of infiltrating immune cells using the ‘CIBERSORT’ R package (CIBERSORT R script v1.03). In addition, immune scores and stromal scores were obtained by the ESTIMATE algorithm using the ‘ESTIMATE’ R package (R package, v1.1.0). Tumor purity and infiltrating immune cells in Asian gastric cancer patients were analyzed in low- and high-risk groups. Finally, a model of the nomogram based on the final prediction model was constructed.

Quantitative reverse-transcriptase PCR

Total RNA was extracted from tissues using TRIzol (Invitrogen) reagent. The integrity of isolated RNA was assessed by NanoDrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Next, RNA was reversely transcribed into cDNA using (Vazyme, Nanjing, China) according to the manufacturer’s protocol. RT-qPCR was then performed in the LightCycler96 (Roche) using SYBR Green technology. Relative lncRNA expression level was normalized to GAPDH and calculated using the 2-∆∆Ct method. The primers for GAPDH and lncRNAs were purchased from Generay (Nanjing, China) and their sequences are shown in supplementary Table S2.

Statistical analysis

All statistical analyses were performed using the R 3.6.3 and GraphPad Prism 7 software. Survival curves were plotted by Kaplan-Meier method, and differences in the survival rates between high-risk and low-risk groups were evaluated using the log-rank test. ROC curve and AUC were plotted by ‘SurvivalROC’ package in R. P-value < 0.05 was considered statistically significant.

Results

The level of m6A modification and expression of m6A regulators

We first evaluated the level of m6A modification in GC tissues using dot blot analysis. Results showed that expression of m6A level were higher in 11/12 patient carcinoma tissues compared to paracarcinoma tissues (Fig. 1 A). Then we quantified the m6A levels with densitometry from Fig. 1 A for correlation analysis. Correlation analysis also indicated that the tumor grade and clinical stage of patients were higher with increasing m6A level (R = 0.485, P = 0.111), especially in clinical stage (R = 0.581, P = 0.048) (Fig. 1 B). Given that m6A writers are important for cancer progression, usually as a tumor promoter during tumor progression, this study explored the expression of m6A writers, including METTL3, METTL14, and WTAP, in tumor and adjacent normal tissues. In 11/12 patient tumors, the expression of METTL3 and WTAP were higher compared to adjacent normal samples, and in 10/12 patient tumors, the expression of METTL14 were higher (Fig. 1 C). Overall, these results showed that the expression of METTL3, METTL14, and WTAP were higher in tumor samples than in the adjacent normal samples (supplementary Fig. S1).

Fig. 1
figure 1

Expression profiles of m6A writer protein, its methylation levels in gastric cancer samples and correlation with clinical information. A Levels of m6A methylation in from gastric tumor samples alongside normal controls. B Correlation between levels of m6A methylation with grade, and tumor stages in gastric cancer samples. C Expression levels of METTL3, METTL14 and WTAP proteins in twelve paired gastric cancer tissue samples

Full size image

Identification of m6A-related lncRNAs in Asian GC patients

The matrix expression of 23 m6A genes and 1054 m6A-related lncRNAs was retrieved from the TCGA database (Fig. 2 A, supplementary Fig. S2 and supplementary Table S1). Univariate Cox regression analysis was the used to screen and identify the m6A-related lncRNAs associated with overall survival. The candidate lncRNAs were classified into two types: a risk type (AC022034.1, AC022034.4, AL133346.1, AC135012.3, AC005586.2, AC104083.1, AC026691.1, AC008808.2, AP001528.1, LINC01140, AC090825.1, and AL133355.1) with HR > 1, which was associated with poor prognosis, and a protective type (SNHG3, AC099850.4, MIR17HG, AL033527.3, AC091057.1, AC026333.4, AL512506.1 and SCAT2) with HR < 1, which was associated with good prognosis (Fig. 2 B, supplementary Table S1). Among the lncRNAs, eight were overexpressed in tumor tissues (SNHG3, AC099850.4, MIR17HG, AL033527.3, AC091057.1, AC026333.4, AL512506.1, and SCAT2) (Fig. 2 C and D). We clustered the 72 patients who were included in this study into two clusters using the ‘ConsensusClusterPlus’ package in R. When k = 2, two molecular subtypes were conducted (Fig. 2 E). The overall survival results showed that the patients’ prognosis was significantly different between the two subtypes (Fig. 2 E), with patients in cluster 2 exhibiting a significantly worse prognosis than those in cluster 1 (P = 0.018; Fig. 2 E). This suggested that the clinical results of the two subtypes were significantly different. However, the heatmap and clinic pathologic features showed that there were no significant differences in other clinical factors, including age, gender, grade, and TNM stage (Fig. 2 F).

Fig. 2
figure 2

Screening and genotyping survival analysis of m6A-related lncRNAs related to the prognosis of gastric cancer. A A co-expression network of LncRNA-mRNAs. Blue and yellow nodes represent m6A-related mRNAs and lncRNAs, respectively. B Forest plot showing m6A-related lncRNAs significantly associated with prognosis of gastric cancer patients. Cyan and magenta fonts denote favorable and adverse prognostic factors, respectively. C Expression profiles of candidate lncRNAs in TCGA gastric cancer tissues. Red and blue boxes represent tumor and normal tissues, respectively, whereas cyan and magenta fonts denote downregulated and upregulated lncRNAs, respectively. D Heatmap showing differentially expressed lncRNAs between normal and tumor tissues. Magenta and cyan color show highly and lowly expressed lncRNAs, respectively. *P < 0.05, **P < 0.01, ***P < 0.001, relative to corresponding normal tissues. E Identification of two molecular subtypes based on the candidate lncRNAs and Kaplan-Meier survival curve analysis. F Heatmap and clinicopathologic features of the two clusters defined by the prognosis-related m6A-related lncRNAs expression

Full size image

Distribution of tumor-infiltrating immune cells in Asian gastric cancer patients

Next, we compared the expressions of four important immune checkpoints (PD-L1, CTLA4, IL4I1, and IDO1) between normal and tumor tissues, and between the two clusters. Results indicated that the four immune checkpoints exhibited increased expression in GC tumors. However, there was no significant difference between the two molecular subtypes (Fig. 3). Results also revealed a positive correlation between the expression of PD-L1 and AC099850.4 (P < 0.05). There was a positive correlation between the expression of CTLA4 and AC090825.1 and AL133355.1 (P < 0.05), and a negative correlation between the expression of IL4I1 and AL512506.1 (P < 0.05). Furthermore, the expression of IDO1 was significantly positively correlated with AC099850.4 and AC026333.4 (P < 0.05). Collectively, the results suggest that these lncRNAs were all tightly associated with tumor-immune microenvironment.

Fig. 3
figure 3

Expression profiles of four important immune checkpoints in gastric cancer tissues, and their correlation with prognosis-related m6A-related lncRNAs. Differences in the expression of PD-L1, CTLA4, IL4I1 and IDO1 in different tissue types (A, D, G, J) and various subtypes (B, E, H, K). Correlation analysis of the relationship between the expression of the target genes PD-L1, CTLA4, IL4I1 and IDO1 and that of prognosis-related m6A-related lncRNAs (C, F, I, L). Red indicates positive relationships, while blue indicates negative relationships. *P < 0.05, ***P < 0.001

Full size image

In addition, the proportions of 22 types of immune cells were calculated using the CIBERSORT algorithm (Fig. 4 A). Patients in cluster 2 exhibited a worse prognosis, with increased immune cells (resting memory CD4+ T cells, P = 0.027; regulatory T cells, P = 0.0018; monocytes, P = 0.00095; and resting dendritic cells, P = 0.015), and decreased immune cells (resting NK cells, P = 0.033; and macrophages M1, P = 0.045) (Fig. 4 B). After ESTIMATE algorithm was processed, individuals in cluster 2 exhibited higher estimate score (P = 1.1e-6), higher stromal score (P = 6.3e-9), and higher immune score (P = 0.0012) than individuals in cluster 1 (Fig. 4 C).

Fig. 4
figure 4

Proportions of 22 types of immune cells analyzed by the CIBERSORT algorithm, and estimate or stromal scores determined by the ESTIMATE algorithm. A Distribution and expression of 22 types of immune cells across the 2 molecular clusters. B Significantly differentially expressed immune cells in the two clusters. C ESTIMATE, stromal and immune scores in two clusters

Full size image

Functional annotation was further performed through gene set enrichment analysis (GSEA). Results demonstrated that the genes involved in ECM_RECEPTOR_INTERACTION (NES = 2.10; FDR = 0.003), FOCAL_ADHESION (NES = 2.11; FDR = 0.004), CELL_ADHESION_MOLECULES_CAMS (NES = 1.93; FDR = 0.026) and MAPK_SIGNALING_PATHWAY (NES = 1.89; FDR = 0.027) were significantly enriched in cluster 2 individuals (Fig. 5 A). Meanwhile, the genes involved in BASE_EXCISION_REPAIR (NES = -2.15; FDR = 0.000), SPLICEOSOME (NES = -2.20; FDR = 0.000), HOMOLOGOUS_RECOMBINATION (NES = -1.96; FDR = 0.009), PROTEASOME (NES = -2.09; FDR = 0.002), and OXIDATIVE_PHOSPHORYLATION (NES = -1.99; FDR = 0.007) were significantly enriched in cluster 1 individuals (Fig. 5 B).

Fig. 5
figure 5

Gene set enrichment analysis (GSEA) results for Asian GC patients from TCGA datasets. A Tumor hallmarks were significantly enriched in the cluster 2 subgroup. B Significant enrichment of pathways related to cell cycle and cell damage repair in cluster 1 subgroup

Full size image

Construction of the six m6A-related lncRNAs signature as an indicator for prognosis

LASSO COX regression analysis was performed to optimize the prognostic model and prevent overfitting (supplementary Fig. S3). Finally, six m6A-related lncRNAs were screened as independent potential prognostic factors for OS and their risk coefficients are listed in Table 2. A prognostic model was constructed based on the expression of SNHG3, AC026333.4, AL133346.1, AL512506.1, AC026691.1, and AC090825.1 (Table 2, supplementary Fig. S4). The risk score of gastric cancer prognosis was quantified using the following formula: risk score = SNHG3× (− 0.0203) + AC026333.4× (− 0.5629) + AL133346.1× (0.2539) + AL512506.1× (− 0.2812) + AC026691.1× (1.1868) + AC090825.1× (1.0361). Patients in the training set were then stratified into high-risk and low-risk groups based on the median risk score. Similarly, patients in the test set were also divided into low-risk and high-risk groups based on the cut-offs in the training set. Results showed that low-risk patients had longer OS than high-risk patients in both training set and test set (Fig. 6 A and D). ROC curve analysis showed that sensitivity was highest when AUC was 0.828 in the training set and 0.785 in the test set (Fig. 6 B and E). Moreover, the risk score and survival status predicted by the prognostic model were displayed in Fig. 6 C and F. From the results, it was evident that this model presented a good prediction power.

Table 2 Information on six prognostic m6A-related lncRNAs weighted by their multivariable LASSO regression coefficients
Full size table
Fig. 6
figure 6

Construction and prognostic value of a nine m6A-related lncRNAs signature. A The overall survival rates of high- and low-risk groups in train group (N = 36). B Time-dependent ROC curves, based on risk scores, for predicting overall survival in the training dataset. C Distribution of risk scores, survival status and a heatmap of nine prognostic m6A-related lncRNAs in the training group. D Overall survival of patients in the high- and low-risk groups based on test dataset (N = 36). E Time-dependent ROC curves, based on risk scores, for predicting overall survival in the test group. F Distribution of risk scores, survival status and a heatmap of nine prognostic m6A-related lncRNAs in the test group

Full size image

Relationships between the risk scores and clinical factors ​

Herein, we built a prognostic model based on the entire set to evaluate the risk score in different stratification. The model showed a good performance in the stratification of individuals aged less than 65 or more than 65, gender, grade G1–2&G3, T3–4, M0, N0&N1–3, and clinical stage I–II & III–IV (Fig. 7). These results showed a good performance in stratifying age ≤ 65 (P = 0.059) and > 65 (P < 0.001), male (P = 0.003) and female (P = 0.185), grade G1–2 (P = 0.053) and G3 (P = 0.008), T3–4 (P < 0.001), M0 (P < 0.001), N0 (P = 0.021) and N1–3 (P = 0.002), clinical-stage I-II (P = 0.006), and III-IV (P = 0.013) (Fig. 7). The high-risk group in both subgroups was associated with worse OS than the low-risk group in Asian gastric cancer patients, which was consistent with the above-described results. Figure 8 A shows the heatmap of the correlation analysis between risk scores and clinical traits. It was found that the gender, cluster, and immunescore had a significant positive correlation with the risk score. In addition, patients in the high-risk score group had higher immune scores (P = 0.0075), whereas all patients in cluster 2 had high-risk scores (P = 8.6e-07), which indicates worse prognostic results and suggests that patients in cluster 2 subgroup might benefit from immunotherapy. Meanwhile, results revealed that male patients had significantly higher risk scores than female patients (Fig. 8 B). Although, there was only difference in the expression of IL4I1 (P = 0.019), and no difference in the expression of other immune checkpoints (PD-L1, CTLA4 and IDO1) between the high- and low-risk score groups (Fig. 8 C), the results showed that the level of resting memory CD4+ T cells (R = 0.31, P = 0.038), regulatory T cells (Tregs) (R = 0.42, P = 0.0042), monocytes (R = 0.42, P = 0.0043) had a positive correlation with risk scores, whereas M1 macrophages (R = -0.37, P = 0.012) and resting NK cells (R = -0.31, P = 0.043) had a negative correlation with risk scores (Fig. 8 D).

Fig. 7
figure 7

Kaplan-Meier survival curves for Asian GC patients in the TCGA dataset. Stratified analyses of clinicopathological factors in gastric cancer: age, gender, grade G1–2&G3, T3–4, M0, N0 & N1–3, and stage I-II & III-IV

Full size image
Fig. 8
figure 8

A heat map and box plot showing the relationship between risk scores with clinical features and immune cells. A A heatmap showing the correlation between risk scores and clinical characteristics. B Box plot showing the correlation between risk scores and clinical characteristics. C Profiles of PD-L1, CTLA4, IL4I1, and IDO1 expression in the high- and low-risk groups. D The correlation between risk scores with expression of immune cells, including resting memory CD4+ T cells, regulatory T cells, Monocytes, M1 macrophages and resting NK cells

Full size image

Construction of nomogram and its clinical utility

A nomogram was constructed to incorporate the risk score and clinical factors, and its accuracy was validated using the calibration curve of 1-year and 2-year survival (Fig. 9 A). Basic information and clinicopathological parameters of gastric cancer patients included in the study are shown in supplementary Table S3. qRT-PCR results validated the expression of the six predicted m6A-related lncRNAs (SNHG3, AC026333.4, AL133346.1, AL512506.1, AC026691.1, and AC090825.1) in clinical tissue samples (Fig. 9 B). Furthermore, a weak positive correlation was found between the risk score calculated by our model and stage of gastric cancer patients (R = 0.384, P = 0.218) (Fig. 9 C).

Fig. 9
figure 9

A novel nomogram for predicting overall survival rates of Asian GC patients based on independent prognostic factors and the nine m6A-related lncRNAs expression between normal and tumor tissues. A Nomogram for predicting overall survival probability of GC patients. B Expression profiles of six m6A-related lncRNAs in twelve paired gastric cancer tissue samples. C. Correlation between risk scores calculated by model with GC stage in patients

Full size image

Discussion

Currently, gastric cancer is a major gastrointestinal malignancy worldwide, with high incidence rates in East Asia [31]. It is worth noting that distinguishing the risk factors before treatment is important in improving patients’ outcomes. Therefore, this calls for identification of reliable biomarkers and establishing powerful prognostic models, with the overarching goal of improving the OS for GC patients. Although some models for predicting patient survival rates have been identified in recent years, they all have several limitations as described in our previous study [30]. Given the better stability and sensitivity of lncRNAs, they have become key biomarkers or drug treatment targets for GC patients [32, 33]. In many tumors, the lncRNAs signature is often used to predict the prognosis of tumor patients [34]. M6A methylation is the most prevalent post-transcriptional modification of lncRNAs, which affects RNA stability and translational activation [35]. Despite the elucidation of lncRNAs functions in gastric cancer, there is still less findings on the prognostic prediction model of lncRNAs’ m6A methylation modification.

Changes in the m6A modification have previously been implicated in the progression of oncogenesis, especially in gastric cancer [36]. Several studies have reported that progression of GC is promoted by m6A writers, for example, METTL3-mediated m6A modification of some mRNA (HDGF [24], SPHK2 [37], and MYC [38]), which ultimately promotes gastric cancer progression, thus, exhibiting a prognostic significance. Meanwhile, many of the current studies are focusing on the role of m6A readers in lung cancer. One study found that m6A reader YTHDF1 facilitates the tumorigenesis and metastasis of GC via USP14 translation [39]. Knockdown of m6A reader IGF2BP3 inhibited hypoxia-induced cell migration and angiogenesis by regulating hypoxia inducible factor-1alpha in GC [40]. These findings suggest that abnormally high levels of m6A modification in GC patients are often positively correlated with poor prognosis, which is associated with high expression of writers and readers. This fact was also confirmed by our clinical tissue protein results and clinical correlation tests. M6A modification not only affects the functions of mRNAs, but also the lncRNAs. Notably, the functional link between lncRNA GAS5 and m6A modification in colorectal cancer has already been proven, where the lncRNA GAS5, YAP signaling, and YTHDF3 formed a negative feedback loop that promoted cancer progression [41]. The m6A-mediated upregulation of LINC01320 promotes the proliferation, migration, and invasion of gastric cancer via the miR495-5p/RAB19 axis [42]. However, although recent studies have paid more attention on the correlation between m6A modification and lncRNAs, further research is still needed.

Herein, it was found that higher level of m6A modification and selected lncRNAs in human gastric tumor were associated with poor OS. Notably, we screened out six lncRNAs, including SNHG3, AC026333.4, AL133346.1, AL512506.1, AC026691.1, and AC090825.1. LncRNA SNHG3 is a oncogenic and an abnormally expressed lncRNA in various tumors such as gastric cancer [43] and lung cancer [44]. Aberrant expression of SNHG3 is correlated with worse outcomes of patients through contributing to tumor cell proliferation, migration, and invasion [45]. A previous study proposed that SNHG3 could be used as an independent prognostic biomarker of intrahepatic cholangiocarcinoma [46]. AL133346.1 and neighboring protein CCN2 were high-expressed in leukemia samples and could be used to predict the prognosis of patients with B-cell acute lymphoblastic leukemia [47]. AL512506.1 had been found as m6A-related lncRNA and used to establish risk assessment model with other lncRNAs in gastric cancer patients [48]. These findings are in line with our results, which provides further evidence on the potential functions of those lncRNAs and their potential capacities as biomarkers in cancer treatment. However, to date, studies about their molecular mechanism in Asian gastric cancer patients are rare.

After analyzing the relationship between lncRNAs and immune checkpoints, it was found that lncRNAs were closely associated with immune scores. Several immune factors, including PD-L1, CTLA4, IL4I1, and IDO1, were significantly enriched in GC tumors and were correlated with the six identified lncRNAs. The enrichment of CTLA4 and IL4I1 were consistent with their functions in cancer [49]. For example, programmed cell death 1 (PD-L1), cytotoxic T lymphocyte antigen 4 (CTLA4), and the indoleamine 2, 3-dioxygenase 1 (IDO1) have achieved impressive success in the treatment of different cancer types, especially the combined treatment of PD-L1 and CTLA4 [49]. In a previous study, CTLA4, PD-1, and PD-L1, as the checkpoint molecules, could be targeted by antibodies to slow cancer progression [50]. CTLA4 could affect diversity of T cells through increased T cell priming, and acted on both cytotoxic T cells and helper T cells [51]. In gastric cancer, significant benefits were observed after inhibiting PD-L1 and CTLA4 [52]. On the other side, the expression of interleukin-4induced-1 (IL4I1) in tumors has also been observed very frequently, especially in tumor-associated macrophages [53]. IL4I1 can activate aryl hydrocarbon receptor (AHR) and promote progression of cancer by enhancing chronic lymphocytic leukemia progression [54]. Similarly, IDO1 can promote AHR-driven processes and is associated with worse outcomes for tumor patients [55]. However, the function of these genes in cancer is still controversial, and thus more studies should be conducted to clarify their functions. Meanwhile, the poorer prognosis may be associated with the abnormally high expression of three immune cells (resting memory CD4+ T cells, regulatory T cells and monocytes) and down expression of two immune cells (M1 macrophages and resting NK cells), especially in Asian males with gastric cancer. Activated regulatory T cells accumulated in tumor microenvironment and correlated with tumor progression in patients with colorectal cancer [56]. Monocytes and macrophages are the major sources of numerous cytokines, including tumor necrosis factor alpha (TNF), interleukins, colony-stimulating factors, chemokines and cytotoxic mediators, which appear to play an important role in the regulation of tumor growth [57]. Notably, tumor-infiltrating NK cells are mostly “resting” non-activated and a component of innate immunity [58]. Studies report that reduced NK cell activity in gastric cancer patients is associated with poorer prognosis [59]. Collectively, the findings of this study suggest that our prognosis prediction model can identify Asian gastric cancer patients who may benefit from immunotherapy, thereby prolonging their overall survival. Further validation of the six m6A-related lncRNAs signature in a large sample population is necessary and more data is be needed to validate the general applicability of this signature in clinical decisions.

Conclusions

This study has combined clinical characteristics with the expression of six m6A-related lncRNAs, including SNHG3, AC026333.4, AL133346.1, AL512506.1, AC026691.1, and AC090825.1, to construct an accurate prognosis prediction model for predicting the overall survival of Asian GC patients.

Availability of data and materials

The data sets used and/or analyzed during the current study are publicly available data from The Cancer Genome Atlas (TCGA, http://cancergenome.nih.gov/). The figures and materials supporting the conclusions of this article are included within the article.

Abbreviations

GC:

Gastric cancer

TCGA:

The cancer genome atlas

m6A:

N6-methyladenosine

GSEA:

Gene set enrichment analysis

LncRNAs:

Long noncoding RNAs

CEA:

Carcinoembryonic antigen

AFP:

Alpha-fetoprotein

CA:

Carbohydrate antigen

ROC:

The receiver operating characteristic curve

AUC:

The area under the curve

OS:

Overall survival

PD-L1:

Programmed cell death 1

CTLA4:

Cytotoxic T lymphocyte antigen 4

IDO1:

Indoleamine 2, 3-dioxygenase 1

IL4I1:

Interleukin-4 induced-1

References

  1. GlobalSurg C. National Institute for Health Research Global Health research unit on global S: global variation in postoperative mortality and complications after cancer surgery: a multicentre, prospective cohort study in 82 countries. Lancet. 2021;397(10272):387–97.

    Article  Google Scholar 

  2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  3. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.

    Article  PubMed  Google Scholar 

  4. Mills JC, Samuelson LC. Past questions and current understanding about gastric Cancer. Gastroenterology. 2018;155(4):939–44.

    Article  PubMed  Google Scholar 

  5. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  6. Song Z, Zou S, Zhou W, Huang Y, Shao L, Yuan J, et al. Clinically applicable histopathological diagnosis system for gastric cancer detection using deep learning. Nat Commun. 2020;11(1):4294.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Thrift AP, El-Serag HB. Burden of gastric Cancer. Clin Gastroenterol Hepatol. 2020;18(3):534–42.

    Article  PubMed  Google Scholar 

  8. Wagner AD, Syn NL, Moehler M, Grothe W, Yong WP, Tai BC, et al. Chemotherapy for advanced gastric cancer. Cochrane Database Syst Rev. 2017;8:CD004064.

    PubMed  Google Scholar 

  9. Joshi SS, Badgwell BD. Current treatment and recent progress in gastric cancer. CA Cancer J Clin. 2021;71(3):264–79.

    Article  PubMed  Google Scholar 

  10. Shimada H, Noie T, Ohashi M, Oba K, Takahashi Y. Clinical significance of serum tumor markers for gastric cancer: a systematic review of literature by the task force of the Japanese gastric Cancer association. Gastric Cancer. 2014;17(1):26–33.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang E, He X, Zhang C, Su J, Lu X, Si X, et al. A novel long noncoding RNA HOXC-AS3 mediates tumorigenesis of gastric cancer by binding to YBX1. Genome Biol. 2018;19(1):154.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Shuai Y, Ma Z, Liu W, Yu T, Yan C, Jiang H, et al. TEAD4 modulated LncRNA MNX1-AS1 contributes to gastric cancer progression partly through suppressing BTG2 and activating BCL2. Mol Cancer. 2020;19(1):6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Zhang G, Chi N, Lu Q, Zhu D, Zhuang Y. LncRNA PTCSC3 is a biomarker for the treatment and prognosis of gastric Cancer. Cancer Biother Radiopharm. 2020;35(1):77–81.

    Article  CAS  PubMed  Google Scholar 

  14. Peng C, Li X, Yu Y, Chen J. LncRNA GASL1 inhibits tumor growth in gastric carcinoma by inactivating the Wnt/beta-catenin signaling pathway. Exp Ther Med. 2019;17(5):4039–45.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Li S, Zhang M, Zhang H, Hu K, Cai C, Wang J, et al. Exosomal long noncoding RNA lnc-GNAQ-6:1 may serve as a diagnostic marker for gastric cancer. Clin Chim Acta. 2020;501:252–7.

    Article  CAS  PubMed  Google Scholar 

  16. Cai C, Zhang H, Zhu Y, Zheng P, Xu Y, Sun J, et al. Serum Exosomal long noncoding RNA pcsk2-2:1 as a potential novel diagnostic biomarker for gastric Cancer. Onco Targets Ther. 2019;12:10035–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Li Y, Jiang T, Zhou W, Li J, Li X, Wang Q, et al. Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers. Nat Commun. 2020;11(1):1000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Rauch S, Jones KA, Dickinson BC. Small molecule-inducible RNA-targeting Systems for Temporal Control of RNA regulation. ACS Cent Sci. 2020;6(11):1987–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hu BB, Wang XY, Gu XY, Zou C, Gao ZJ, Zhang H, et al. N(6)-methyladenosine (m(6)a) RNA modification in gastrointestinal tract cancers: roles, mechanisms, and applications. Mol Cancer. 2019;18(1):178.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Fang R, Chen X, Zhang S, Shi H, Ye Y, Shi H, et al. EGFR/SRC/ERK-stabilized YTHDF2 promotes cholesterol dysregulation and invasive growth of glioblastoma. Nat Commun. 2021;12(1):177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Gu C, Wang Z, Zhou N, Li G, Kou Y, Luo Y, et al. Mettl14 inhibits bladder TIC self-renewal and bladder tumorigenesis through N(6)-methyladenosine of Notch1. Mol Cancer. 2019;18(1):168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yin H, Zhang X, Yang P, Zhang X, Peng Y, Li D, et al. RNA m6A methylation orchestrates cancer growth and metastasis via macrophage reprogramming. Nat Commun. 2021;12(1):1394.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lin X, Chai G, Wu Y, Li J, Chen F, Liu J, et al. RNA m(6)a methylation regulates the epithelial mesenchymal transition of cancer cells and translation of snail. Nat Commun. 2019;10(1):2065.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Wang Q, Chen C, Ding Q, Zhao Y, Wang Z, Chen J, et al. METTL3-mediated m(6)a modification of HDGF mRNA promotes gastric cancer progression and has prognostic significance. Gut. 2020;69(7):1193–205.

    Article  CAS  PubMed  Google Scholar 

  25. Pan T. N6-methyl-adenosine modification in messenger and long non-coding RNA. Trends Biochem Sci. 2013;38(4):204–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Chen YG, Chen R, Ahmad S, Verma R, Kasturi SP, Amaya L, et al. N6-Methyladenosine modification controls circular RNA immunity. Mol Cell. 2019;76(1):96–109 e109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Huang H, Weng H, Chen J. M(6)a modification in coding and non-coding RNAs: roles and therapeutic implications in Cancer. Cancer Cell. 2020;37(3):270–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Liu N, Parisien M, Dai Q, Zheng G, He C, Pan T. Probing N6-methyladenosine RNA modification status at single nucleotide resolution in mRNA and long noncoding RNA. RNA. 2013;19(12):1848–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Wang C, Zhu J, Zhang Z, Chen H, Ji M, Chen C, et al. Rno-miR-224-5p contributes to 2,2′,4,4′-tetrabromodiphenyl ether-induced low triiodothyronine in rats by targeting deiodinases. Chemosphere. 2020;246:125774.

    Article  CAS  PubMed  Google Scholar 

  30. Xia R, Tang H, Shen J, Xu S, Liang Y, Zhang Y, et al. Prognostic value of a novel glycolysis-related gene expression signature for gastrointestinal cancer in the Asian population. Cancer Cell Int. 2021;21(1):154.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kong R, Zhang EB, Yin DD, You LH, Xu TP, Chen WM, et al. Long noncoding RNA PVT1 indicates a poor prognosis of gastric cancer and promotes cell proliferation through epigenetically regulating p15 and p16. Mol Cancer. 2015;14:82.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Cai C, Yang L, Zhuang X, He Y, Zhou K. A five-lncRNA model predicting overall survival in gastric cancer compared with normal tissues. Aging (Albany NY). 2021;13(21):24349–59.

    Article  CAS  Google Scholar 

  33. Ding Z, Li R, Han J, Sun D, Shen L, Wu G. Identification of an immune-related LncRNA signature in gastric Cancer to predict survival and response to immune checkpoint inhibitors. Front Cell Dev Biol. 2021;9:739583.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Cao W, Liu JN, Liu Z, Wang X, Han ZG, Ji T, et al. A three-lncRNA signature derived from the atlas of ncRNA in cancer (TANRIC) database predicts the survival of patients with head and neck squamous cell carcinoma. Oral Oncol. 2017;65:94–101.

    Article  CAS  PubMed  Google Scholar 

  35. Patil DP, Pickering BF, Jaffrey SR. Reading m(6)a in the transcriptome: m(6)A-binding proteins. Trends Cell Biol. 2018;28(2):113–27.

    Article  CAS  PubMed  Google Scholar 

  36. Jaffrey SR, Kharas MG. Emerging links between m(6)a and misregulated mRNA methylation in cancer. Genome Med. 2017;9(1):2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Huo FC, Zhu ZM, Zhu WT, Du QY, Liang J, Mou J. METTL3-mediated m(6)a methylation of SPHK2 promotes gastric cancer progression by targeting KLF2. Oncogene. 2021;40(16):2968–81.

    Article  CAS  PubMed  Google Scholar 

  38. Yang Z, Jiang X, Li D, Jiang X. HBXIP promotes gastric cancer via METTL3-mediated MYC mRNA m6A modification. Aging (Albany NY). 2020;12(24):24967–82.

    Article  CAS  Google Scholar 

  39. Chen XY, Liang R, Yi YC, Fan HN, Chen M, Zhang J, et al. The m(6)a reader YTHDF1 facilitates the tumorigenesis and metastasis of gastric Cancer via USP14 translation in an m(6)A-dependent manner. Front Cell Dev Biol. 2021;9:647702.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Jiang L, Li Y, He Y, Wei D, Yan L, Wen H. Knockdown of m6A reader IGF2BP3 inhibited hypoxia-induced cell migration and angiogenesis by regulating hypoxia inducible factor-1alpha in stomach Cancer. Front Oncol. 2021;11:711207.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ni W, Yao S, Zhou Y, Liu Y, Huang P, Zhou A, et al. Long noncoding RNA GAS5 inhibits progression of colorectal cancer by interacting with and triggering YAP phosphorylation and degradation and is negatively regulated by the m(6) a reader YTHDF3. Mol Cancer. 2019;18(1):143.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Hu N, Ji H. N6-methyladenosine (m6A)-mediated up-regulation of long noncoding RNA LINC01320 promotes the proliferation, migration, and invasion of gastric cancer via miR495-5p/RAB19 axis. Bioengineered. 2021;12(1):4081–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Xuan Y, Wang Y. Long non-coding RNA SNHG3 promotes progression of gastric cancer by regulating neighboring MED18 gene methylation. Cell Death Dis. 2019;10(10):694.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Liu L, Ni J, He X. Upregulation of the long noncoding RNA SNHG3 promotes lung adenocarcinoma proliferation. Dis Markers. 2018;2018:5736716.

    PubMed  PubMed Central  Google Scholar 

  45. Xu B, Mei J, Ji W, Bian Z, Jiao J, Sun J, et al. LncRNA SNHG3, a potential oncogene in human cancers. Cancer Cell Int. 2020;20(1):536.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Tian D, Wei X, Zhu H, Zhu L, Li T, Li W. LncRNA-SNHG3 is an independent prognostic biomarker of intrahepatic cholangiocarcinoma. Int J Clin Exp Pathol. 2019;12(7):2706–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Cuadros M, Garcia DJ, Andrades A, Arenas AM, Coira IF, Balinas-Gavira C, et al. LncRNA-mRNA co-expression analysis identifies AL133346.1/CCN2 as biomarkers in pediatric B-cell acute lymphoblastic leukemia. Cancers (Basel). 2020;12(12):3803.

    Article  CAS  Google Scholar 

  48. Wang JM, Li X, Yang P, Geng WB, Wang XY. Identification of a novel m6A-related lncRNA pair signature for predicting the prognosis of gastric cancer patients. BMC Gastroenterol. 2022;22(1):76.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Havel JJ, Chowell D, Chan TA. The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy. Nat Rev Cancer. 2019;19(3):133–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Kalbasi A, Ribas A. Tumour-intrinsic resistance to immune checkpoint blockade. Nat Rev Immunol. 2020;20(1):25–39.

    Article  CAS  PubMed  Google Scholar 

  51. Gide TN, Quek C, Menzies AM, Tasker AT, Shang P, Holst J, et al. Distinct Immune Cell Populations Define Response to Anti-PD-1 Monotherapy and Anti-PD-1/Anti-CTLA-4 Combined Therapy. Cancer Cell. 2019;35(2):238–255 e236.

    Article  CAS  PubMed  Google Scholar 

  52. Smyth E, Thuss-Patience PC. Immune checkpoint inhibition in gastro-Oesophageal Cancer. Oncol Res Treat. 2018;41(5):272–80.

    Article  CAS  PubMed  Google Scholar 

  53. Carbonnelle-Puscian A, Copie-Bergman C, Baia M, Martin-Garcia N, Allory Y, Haioun C, et al. The novel immunosuppressive enzyme IL4I1 is expressed by neoplastic cells of several B-cell lymphomas and by tumor-associated macrophages. Leukemia. 2009;23(5):952–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Sadik A, Somarribas Patterson LF, Ozturk S, Mohapatra SR, Panitz V, Secker PF, et al. IL4I1 Is a Metabolic Immune Checkpoint that Activates the AHR and Promotes Tumor Progression. Cell. 2020;182(5):1252–1270 e1234.

    Article  CAS  PubMed  Google Scholar 

  55. Mitchell TC, Hamid O, Smith DC, Bauer TM, Wasser JS, Olszanski AJ, et al. Epacadostat plus Pembrolizumab in patients with advanced solid tumors: phase I results from a multicenter, open-label phase I/II trial (ECHO-202/KEYNOTE-037). J Clin Oncol. 2018;36(32):3223–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Lin YC, Mahalingam J, Chiang JM, Su PJ, Chu YY, Lai HY, et al. Activated but not resting regulatory T cells accumulated in tumor microenvironment and correlated with tumor progression in patients with colorectal cancer. Int J Cancer. 2013;132(6):1341–50.

    Article  CAS  PubMed  Google Scholar 

  57. Siedlar M, Szaflarska A, Szczepanik A, Ruggiero I, Frankenberger M, Szatanek R, et al. Depressed tumor necrosis factor alpha and interleukin-12p40 production by peripheral blood mononuclear cells of gastric cancer patients: association with IL-1R-associated kinase-1 protein expression and disease stage. Int J Cancer. 2005;114(1):144–52.

    Article  CAS  PubMed  Google Scholar 

  58. Fiore PF, Vacca P, Tumino N, Besi F, Pelosi A, Munari E, et al. Wilms' tumor primary cells display potent Immunoregulatory properties on NK cells and macrophages. Cancers (Basel). 2021;13(2):224.

    Article  CAS  Google Scholar 

  59. Saito H, Takaya S, Osaki T, Ikeguchi M. Increased apoptosis and elevated Fas expression in circulating natural killer cells in gastric cancer patients. Gastric Cancer. 2013;16(4):473–9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Most of the data in this study were obtained through the Cancer Genome Atlas Database (TCGA). We are grateful for the source of data used in our research.

Funding

This work was supported by the National Natural Science Foundation of China (81903353), the Technology Development Fund of Nanjing Medical University [NMUB2018001] and the project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions and the Collegiate Natural Science Foundation of Jiangsu Province [19KJB330003], and the Scientific Research Funding of Tongling Municipal Health Commission (201912).

Author information

Author notes

  1. Shuyu Xu, Wenlong Chen and Yiwen Wang contributed equally to this work.

Authors and Affiliations

  1. Key Lab of Modern Toxicology of Ministry of Education, Center for Global Health, School of Public Health, Nanjing Medical University, 101 Longmian Avenue, Nanjing, 211166, People’s Republic of China

    Shuyu Xu, Rong Xia, Jiemiao Shen, Xing Gong, Yinyin Liang, Jiayi Xu & Chao Wang

  2. State Key Lab of Reproductive Medicine, Institute of Toxicology, Nanjing Medical University, 101 Longmian Avenue, Nanjing, 211166, People’s Republic of China

    Shuyu Xu, Rong Xia, Jiemiao Shen, Xing Gong, Yinyin Liang, Jiayi Xu & Chao Wang

  3. Department of Colorectal Surgery, the First Affiliated Hospital of Nanjing Medical University, Nanjing, 210000, People’s Republic of China

    Wenlong Chen, Yiwen Wang, Yi Zhang & Tao Chen

  4. The first clinical medical college of Nanjing Medical University, Nanjing, 211166, People’s Republic of China

    Yuxin Zhang

  5. Department of General Surgery, Yueqing People’s Hospital, 338 Qingyuan road, Yueqing, 325600, Zhejiang province, People’s Republic of China

    Hua Tang

  6. Department of General Surgery, Tongling People’s Hospital, 468 Bijiashan Road, Tongling, Anhui Province, 244000, People’s Republic of China

    Tie Zhao

Authors
  1. Shuyu Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Wenlong Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yiwen Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yuxin Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Rong Xia
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jiemiao Shen
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Xing Gong
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Yinyin Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jiayi Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Hua Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Tie Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Yi Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Tao Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Chao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CW, TC and YZ conceived and designed the model. SX, WC and YW analyzed the data and wrote the manuscript. YZ, RX, JS, XG, YL, JX, HT and TZ contributed significantly to the data analyses and experiment. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yi Zhang, Tao Chen or Chao Wang.

Ethics declarations

Ethics approval and consent to participate

This research was approved by the Institutional Review Board of Nanjing Medical University and the Ethical committee of the Tongling People’s Hospital (ethical review No.2019–008) and was conducted in accordance with the Declaration of Helsinki. Informed consent was obtained from all subjects.

Consent for publication

None.

Competing interests

The authors declare no competing financial interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary Figure S1-S4.

Additional file 2:

 Supplementary Table S1.

Additional file 3:

 Supplementary Table S2.

Additional file 4:

 Supplementary Table S3.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, S., Chen, W., Wang, Y. et al. N6-methyladenosine-related lncRNAs identified as potential biomarkers for predicting the overall survival of Asian gastric cancer patients. BMC Cancer 22, 721 (2022). https://doi.org/10.1186/s12885-022-09801-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09801-z

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us