Comprehensive Genomic and Epigenomic Analyses on Transcriptomic Regulation in Stomach Adenocarcinoma

A novel endoplasmic reticulum stress-related lncRNA signature for prognosis prediction and immune response evaluation in Stomach adenocarcinoma

BACKGROUND

Stomach adenocarcinoma (STAD) is a significant contributor to cancer-related mortality worldwide. Although previous research has identified endoplasmic reticulum stress (ERS) as a regulator of various tumor-promoting properties of cancer cells, the impact of ERS-related long non-coding RNAs (lncRNAs) on STAD prognosis has not yet been investigated. Therefore, our study aims to develop and validate an ERS-related lncRNA signature that can accurately predict the prognosis of STAD patients.

METHODS

We collected RNA expression profiles and clinical data of STAD patients from The Cancer Genome Atlas (TCGA) and identified ERS-related genes from the Molecular Signature Database (MSigDB). Co-expression analysis enabled us to identify ERS-related lncRNAs, and we applied univariate Cox, least absolute shrinkage, and selection operator (LASSO), and multivariate Cox regression analyses to construct a predictive signature comprising of 9 ERS-related lncRNAs. We assessed the prognostic accuracy of our signature using Kaplan-Meier survival analysis, and validated our predictive signature in an independent gene expression omnibus (GEO) cohort. We also performed tumor mutational burden (TMB) and tumor immune microenvironment (TIME) analyses. Enrichment analysis was used to investigate the functions and biological processes of the signature, and we identified two distinct STAD patient subgroups through consensus clustering. Finally, we performed drug sensitivity analysis and immunologic efficacy analysis to explore further insights.

RESULTS

The 9 ERS related-lncRNAs signature demonstrated satisfactory predictive performance as an independent prognostic marker and was significantly associated with STAD clinicopathological characteristics. Furthermore, patients in the high-risk group displayed a worse STAD prognosis than those in the low-risk group. Notably, gene set enrichment analysis (GSEA) revealed significant enrichment of extracellular matrix pathways in the high-risk group, indicating their involvement in STAD progression. Additionally, the high-risk group exhibited significantly lower TMB expression levels than the low-risk group. Consensus clustering revealed two distinct STAD patient subgroups, with Cluster 1 exhibiting higher immune cell infiltration and more active immune functions. Drug sensitivity analysis suggested that the low-risk group was more responsive to oxaliplatin, epirubicinl, and other drugs.

CONCLUSION

Our study highlights the crucial regulatory roles of ERS-related lncRNAs in STAD, with significant clinical implications. The 9-lncRNA signature we have constructed represents a reliable prognostic indicator that has the potential to inform more personalized treatment decisions for STAD patients. These findings shed new light on the pathogenesis of STAD and its underlying molecular mechanisms, offering opportunities for novel therapeutic strategies to be developed for STAD patients.

Multiomics Integration at Single-Cell Resolution Using Bayesian Networks: A Case Study in Hepatocellular Carcinoma

Multiomics data integration is one of the leading frontiers of complex disease research and integrative biology. The advances in single-cell sequencing technologies offer yet another crucial dimension in multiomics research. The single-cell studies enable the study and integration of multiomics data simultaneously in the same cell. We report in this study multiomics data integration in single-cell resolution using Bayesian networks (BNs) in a case study of hepatocellular carcinoma (HCC). A BN encodes the conditional dependencies/independencies of variables using a graphical model with an accompanying joint probability. RNA-seq and Reduced Representation Bisulfite Sequencing data were analyzed separately, and copy number variations were estimated by the hidden Markov model method. Several BN models were constructed to reveal omics' causal and associational relationships. These methods were subjected to a validation study using an independent data set. We show the heterogeneity of the multiple cellular layers of HCC at single-cell omics resolution by identifying best-fitted BN models of 295 genes. We also provide novel insights into the multiomics mechanistic relationships in the human lymphocyte antigen class I genes in HCC. To the best of our knowledge, this is the first study to focus on integrating omics data using a machine learning algorithm, BNs, at the single-cell resolution using a case study of HCC.