Gastric cancer risk-scoring system based on analysis of a competing endogenous RNA network

Integrated bioinformatics analysis of noncoding RNAs with tumor immune microenvironment in gastric cancer

In recent years, molecular and genetic research hotspots of gastric cancer have been investigated, including microRNAs, long noncoding RNAs (lncRNAs) and messenger RNA (mRNAs). The study on the role of lncRNAs may help to develop personalized treatment and identify potential prognostic biomarkers in gastric cancer. The RNA-seq and miRNA-seq data of gastric cancer were downloaded from the TCGA database. Differential analysis of RNA expression between gastric cancer samples and normal samples was performed using the edgeR package. The ceRNA regulatory network was visualized using Cytoscape. KEGG pathway analysis of mRNAs in the ceRNA network was performed using the clusterProfiler package. CIBERSORT was used to distinguish 22 immune cell types and the prognosis-related genes and immune cells were determined using Kaplan-Meier and Cox proportional hazard analyses. To estimate these nomograms, we used receiver operating characteristic and calibration curve studies. The ceRNA regulation network of gastric cancer was built in this study, and the genes in the network were analyzed for prognosis. A total of 980 lncRNAs were differentially expressed, of which 774 were upregulated and 206 were downregulated. A survival study identified 15 genes associated with gastric cancer prognosis, including VCAN-AS1, SERPINE1, AL139002.1, LINC00326, AC018781.1, C15orf54, hsa-miR-145. Monocytes and Neutrophils were associated with the survival rate of gastric cancer. Our research uncovers new ceRNA network for the detection, treatment, and monitoring of gastric cancer.

Immune cell death-related lncRNA signature as a predictive factor of clinical outcomes and immune checkpoints in gastric cancer

Background: Immune cell death (ICD) is a type of tumor cell death that has recently been shown to activate and regulate tumor immunity. However, the role of ICD-related long non-coding RNAs (lncRNAs) in gastric cancer remains to be clarified.

Methods: We obtained 375 tumor samples from the Cancer Genome Atlas (TCGA) database and randomly assigned them to training and verification groups. LASSO and Cox regression analysis were utilized to identify ICD-related lncRNAs and establish a risk model. The changes in the immune microenvironment of the two groups were compared by examining the tumor-infiltrating immune cells.

Results: We established a tumor signature based on nine ICD-related lncRNAs. In light of the receiver operating characteristic and Kaplan–Meier curves, the prognostic values of this risk model were verified. Multivariate regression analysis showed that the risk score was an independent risk factor for the prognosis of patients in both the training cohort (HR 2.52; 95% CI: 1.65–3.87) and validation cohort (HR 2.70; 95% CI: 1.54–4.8). A nomogram was developed to predict the 1-, 3-, and 5-year survival of patients with gastric cancer, and the signature was linked to high levels of immunological checkpoint expression (B7-H3, VSIR).

Conclusions: An ICD-related lncRNA signature could predict the immune response and prognosis of patients with gastric cancer. This prognostic signature could be employed to independently monitor the efficacy of immunotherapy for gastric cancer patients.

A novel risk model construction and immune landscape analysis of gastric cancer based on cuproptosis-related long noncoding RNAs

Recent studies have identified cuproptosis, a new mechanism of regulating cell death. Accumulating evidence suggests that copper homeostasis is associated with tumorigenesis and tumor progression, however, the clinical significance of cuproptosis in gastric cancer (GC) is unclear. In this study, we obtained 26 prognostic cuproptosis-related lncRNAs (CRLs) based on 19 cuproptosis-related genes (CRGs) via Pearson correlation analysis, differential expression analysis, and univariate Cox analysis. A risk model based on 10 CRLs was established with the least absolute shrinkage and selection operator (LASSO) Cox regression analysis and multivariate Cox proportional hazards model to predict the prognosis and immune landscape of GC patients from The Cancer Genome Atlas (TCGA). The risk model has excellent accuracy and efficiency in predicting prognosis of GC patients (Area Under Curve (AUC) = 0.742, 0.803, 0.806 at 1,3,5 years, respectively, P < 0.05). In addition, we found that the risk score was negatively correlated with the infiltration of natural killer (NK) cells and helper T cells, while positively correlated with the infiltration of monocytes, macrophages, mast cells, and neutrophils. Moreover, we evaluated the difference in drug sensitivity of patients with different risk patterns. Furthermore, low-risk patients showed higher tumor mutation burden (TMB) and better immunotherapy response than high-risk patients. In the end, we confirmed the oncogenic role of AL121748.1 which exhibited the highest Hazard Ratio (HR) value among 10 CRLs in GC via cellular functional experiments. In conclusion, our risk model shows a significant role in tumor immunity and could be applied to predict the prognosis of GC patients.

Noncoding RNAs-mediated overexpression of KIF14 is associated with tumor immune infiltration and unfavorable prognosis in lung adenocarcinoma

Kinesin family member 14 (KIF14) is potentially oncogenic and acts as a chromokinesin via binding to microtubules and chromatin during the bipolar spindle formation. KIF14 overexpression is a significant prognostic biomarker in various cancers. However, the expression, prognosis, mechanism, and tumor immune regulation of KIF14 in lung adenocarcinoma (LUAD) remain obscure. Our results demonstrated that KIF14 was upregulated in a variety of cancers, including LUAD. High-expression of KIF14 in LUAD was associated with pathological tumor stage, N stage and unfavorable prognosis. Both univariate and multivariate Cox regression results demonstrated that KIF14 was a significant independent risk factor influencing the prognosis of LUAD patients. The most promising upstream ncRNA-associated pathway of KIF14 in LUAD was determined to be GSEC/TYMSOS-hsa-miR-101-3p axis according to the starBase and The Cancer Genome Atlas databases. Furthermore, upregulation of KIF14 in LUAD was positively correlated with tumor mutation burden, microsatellite instability, immune checkpoint-related gene expression, immune cell biomarkers, and tumor immune cell infiltration. This study reveals that ncRNAs-mediated overexpression of KIF14 is associated with tumor immune infiltration and unfavorable prognosis in LUAD.

Identification and validation of a novel cellular senescence-related lncRNA prognostic signature for predicting immunotherapy response in stomach adenocarcinoma

Background: Cellular senescence is a novel hallmark of cancer associated with patient outcomes and tumor immunotherapy. However, the value of cellular senescence-related long non-coding RNAs (lncRNAs) in predicting prognosis and immunotherapy response for stomach adenocarcinoma (STAD) patients needs further investigation.

Methods: The transcriptome and corresponding clinical information of STAD and cellular senescence-related genes were, respectively, downloaded from the Cancer Genome Atlas (TCGA) and CellAge databases. Differential expression analysis and coexpression analysis were performed to obtain cellular senescence-related lncRNAs. Univariate regression analysis and least absolute shrinkage and selection operator (LASSO) Cox analysis were conducted to establish the cellular senescence-related lncRNA prognostic signature (CSLPS). Next, the survival curve, ROC curve, and nomogram were developed to assess the capacity of predictive models. Moreover, principal component analysis (PCA), gene set enrichment analysis (GSEA), tumor microenvironment (TME), tumor mutation burden (TMB), microsatellite instability (MSI), and tumor immune dysfunction and exclusion (TIDE) score analysis were performed between high- and low-risk groups.

Results: A novel CSLPS involving fifteen lncRNAs (REPIN1-AS1, AL355574.1, AC104695.3, AL033527.2, AC083902.1, TYMSOS, LINC00460, AC005165.1, AL136115.1, AC007405.2, AL391152.1, SCAT1, AC129507.1, AL121748.1, and ADAMTS9-AS1) was developed. According to the nomogram, the risk model based on the CSLPS was an independent prognostic factor and could predict 1-, 3-, and 5-year overall survival for STAD patients. GSEA suggested that the high-risk group was mainly associated with Toll-like receptor, JAK/STAT, NOD-like receptor, and chemokine signaling pathways. Further analysis revealed that STAD patients in the low-risk group with better clinical outcomes had a higher TMB, higher proportion of high microsatellite instability (MSI-H), better immune infiltration, and lower TIDE scores.

Conclusion: A fifteen-CSlncRNA prognostic signature could predict survival outcomes, and patients in the low-risk group may be more sensitive to immunotherapy.

Systemic Analysis of the Anticancer Effects of Sijunzi Decoction on Gastric Cancer Based on Network Pharmacology

Objective:

Sijunzi decoction (SJZD) has been used for alleviating peptic ulcer or gastric discomfort, and treating spleen disorders since the Song Dynasty, but its pharmacological effect on human gastric cancer (GC) is still unclear. In this research, a network pharmacology-based strategy was applied to explore active ingredients, potential targets, and molecular mechanisms of SJZD against GC.

Methods:

The active compounds and potential targets of SJZD, as well as GC-associated gene targets, were retrieved from publicly available databases. Bioinformatics approaches were used to assess the network interaction, functional regulation, and signaling pathways between SJZD ingredients and GC targets. The anticancer effects of SJZD against GC were verified in vivo by a mouse subcutaneous model.

Results:

The results of network analysis showed that quercetin was the most active ingredient in SJZD. Several prominent target genes of SJZD were identified, such as AKT1 and STAT3. Gene ontology analysis revealed that the core anti-GC targets of SJZD included transcription factor activity and kinase activity. Pathway enrichment analysis indicated that GC patients could be benefited from SJZD treatment via modulation of signaling pathways related to endocrine system, cancer, and infectious disease. Furthermore, in vivo experiments showed that high-dose SJZD could inhibit GC xenograft tumor growth, reduce GC cell proliferation, induce GC cell apoptosis, and decrease the expression of p-AKT1 and p-STAT3.

Conclusions:

Taken together, our results suggest that SJZD can serve as an effective adjuvant therapeutic agent for GC patients.

Genetic regulation of THBS1 methylation in diabetic retinopathy

BACKGROUND

Diabetic retinopathy (DR) is a common and serious microvascular complication of diabetes mellitus (DM), but its pathological mechanism, especially the formation mechanism of new blood vessels remains unclear. Thrombospondin-1 (THBS1) is a potent endogenous inhibitor of angiogenesis and it was found over expressed in DR in our previous study. Our study aimed to determine whether overexpression of THBS1 is associated with its promoter methylation level, and whether methylation of THBS1 is regulated by genetic variants in DR.

METHODS

Patients diagnosed with DR and DM patients without retinal problems were included in the case-control study. DNA methylation detection of THBS1 by bisulfite sequencing and genotyping of specific SNPs by MassARRAY analysis were performed in the patients recruited from 2019-2020. Real time quantitative PCR was performed to obtain mRNA expression of THBS1 in the patients recruited from August to October 2022. The differentially methylated CpG loci of THBS1 were identified by logistic regression, and associations between 13 SNPs and methylation levels of CpG loci were tested by methylation quantitative trait loci (meQTLs) analysis. Mediation analysis was applied to determine whether CpG loci were intermediate factors between meQTLs and DR.

RESULTS

150 patients diagnosed with DR and 150 DM patients without retinal complications were enrolled in the first recruitment, seven DR patients and seven DM patients were enrolled in the second recruitment. The patients with DR showed promoter hypomethylation of THBS1 (P value = 0.002), and six out of thirty-nine CpG sites within two CpG islands (CGIs) showed hypomethylation(P value < 0.05). THBS1 mRNA expression in peripheral blood was significantly higher in DR patients than in DM patients. Five out of thirteen cis-meQTLs were identified to be associated with CpG sites: rs13329154, rs34973764 and rs5812091 were associated with cis-meQTLs of CpG-4 (P value=0.0145, 0.0095, 0.0158), rs11070177 and rs1847663 were associated with cis-meQTLs of CpG-2 and CpG-3 respectively (P value=0.0201, 0.0275). CpG-4 methylation significantly mediated the effect of the polymorphism rs34973764 on DR (B=0.0535, Boot 95%CI: 0.004~0.1336).

CONCLUSION

THBS1 overexpression is related to THBS1 hypomethylation in patients with DR. DNA methylation may be genetically controlled in DR.