Comprehensive analysis of nine m7G-related lncRNAs as prognosis factors in tumor immune microenvironment of hepatocellular carcinoma and experimental validation

Natural killer (NK) cells-related gene signature reveals the immune environment heterogeneity in hepatocellular carcinoma based on single cell analysis

The early diagnosis of liver cancer is crucial for the treatment and depends on the coordinated use of several test procedures. Early diagnosis is crucial for precision therapy in the treatment of the hepatocellular carcinoma (HCC). Therefore, in this study, the NK cell-related gene prediction model was used to provide the basis for precision therapy at the gene level and a novel basis for the treatment of patients with liver cancer. Natural killer (NK) cells have innate abilities to recognize and destroy tumor cells and thus play a crucial function as the "innate counterpart" of cytotoxic T cells. The natural killer (NK) cells is well recognized as a prospective approach for tumor immunotherapy in treating patients with HCC. In this research, we used publicly available databases to collect bioinformatics data of scRNA-seq and RNA-seq from HCC patients. To determine the NK cell-related genes (NKRGs)-based risk profile for HCC, we isolated T and natural killer (NK) cells and subjected them to analysis. Uniform Manifold Approximation and Projection plots were created to show the degree of expression of each marker gene and the distribution of distinct clusters. The connection between the immunotherapy response and the NKRGs-based signature was further analyzed, and the NKRGs-based signature was established. Eventually, a nomogram was developed using the model and clinical features to precisely predict the likelihood of survival. The prognosis of HCC can be accurately predicted using the NKRGs-based prognostic signature, and thorough characterization of the NKRGs signature of HCC may help to interpret the response of HCC to immunotherapy and propose a novel tumor treatment perspective.

Advancing bladder cancer management: development of a prognostic model and personalized therapy

Background

Bladder cancer (BLCA) was recognized as a significant public health challenge due to its high incidence and mortality rates. The influence of molecular subtypes on treatment outcomes was well-acknowledged, necessitating further exploration of their characterization and application. This study was aimed at enhancing the understanding of BLCA by mapping its molecular heterogeneity and developing a robust prognostic model using single-cell and bulk RNA sequencing data. Additionally, immunological characteristics and personalized treatment strategies were investigated through the risk score.

Methods

Single-cell RNA sequencing (scRNA-seq) data from GSE135337 and bulk RNA-seq data from several sources, including GSE13507, GSE31684, GSE32894, GSE69795, and TCGA-BLCA, were utilized. Molecular subtypes, particularly the basal-squamous (Ba/Sq) subtype associated with poor prognosis, were identified. A prognostic model was constructed using LASSO and Cox regression analyses focused on genes linked with the Ba/Sq subtype. this model was validated across internal and external datasets to ensure predictive accuracy. High- and low-risk groups based on the risk score derived from TCGA-BLCA data were analyzed to examine their immune-related molecular profiles and treatment responses.

Results

Six molecular subtypes were identified, with the Ba/Sq subtype being consistently associated with poor prognosis. The prognostic model, based on basal-squamous subtype-related genes (BSSRGs), was shown to have strong predictive performance across diverse clinical settings with AUC values at 1, 3, and 5 years indicating robust predictability in training, testing, and entire datasets. Analysis of the different risk groups revealed distinct immune infiltration and microenvironments. Generally higher tumor mutation burden (TMB) scores and lower tumor immune dysfunction and exclusion (TIDE) scores were exhibited by the low-risk group, suggesting varied potentials for systemic drug response between the groups. Finally, significant differences in potential systemic drug response rates were also observed between risk groups.

Conclusions

The study introduced and validated a new prognostic model for BLCA based on BSSRGs, which was proven effective in prognosis prediction. The potential for personalized therapy, optimized by patient stratification and immune profiling, was highlighted by our risk score, aiming to improve treatment efficacy. This approach was promised to offer significant advancements in managing BLCA, tailoring treatments based on detailed molecular and immunological insights.

Integrating single-cell and bulk sequencing data to identify glycosylation-based genes in non-alcoholic fatty liver disease-associated hepatocellular carcinoma

Background

The incidence of non-alcoholic fatty liver disease (NAFLD) associated hepatocellular carcinoma (HCC) has been increasing. However, the role of glycosylation, an important modification that alters cellular differentiation and immune regulation, in the progression of NAFLD to HCC is rare.

Methods

We used the NAFLD-HCC single-cell dataset to identify variation in the expression of glycosylation patterns between different cells and used the HCC bulk dataset to establish a link between these variations and the prognosis of HCC patients. Then, machine learning algorithms were used to identify those glycosylation-related signatures with prognostic significance and to construct a model for predicting the prognosis of HCC patients. Moreover, it was validated in high-fat diet-induced mice and clinical cohorts.

Results

The NAFLD-HCC Glycogene Risk Model (NHGRM) signature included the following genes: SPP1, SOCS2, SAPCD2, S100A9, RAMP3, and CSAD. The higher NHGRM scores were associated with a poorer prognosis, stronger immune-related features, immune cell infiltration and immunity scores. Animal experiments, external and clinical cohorts confirmed the expression of these genes.

Conclusion

The genetic signature we identified may serve as a potential indicator of survival in patients with NAFLD-HCC and provide new perspectives for elucidating the role of glycosylation-related signatures in this pathologic process.

M7G-related tumor immunity: novel insights of RNA modification and potential therapeutic targets

RNA modifications play a pivotal role in regulating cellular biology by exerting influence over distribution features and molecular functions at the post-transcriptional level. Among these modifications, N7-methylguanosine (m7G) stands out as one of the most prevalent. Over recent years, significant attention has been directed towards understanding the implications of m7G modification. This modification is present in diverse RNA molecules, including transfer RNAs, messenger RNAs, ribosomal RNAs, and other noncoding RNAs. Its regulation occurs through a series of specific methyltransferases and m7G-binding proteins. Notably, m7G modification has been implicated in various diseases, prominently across multiple cancer types. Earlier studies have elucidated the significance of m7G modification in the context of immune biology regulation within the tumor microenvironment. This comprehensive review culminates in a synthesis of findings related to the modulation of immune cells infiltration, encompassing T cells, B cells, and various innate immune cells, all orchestrated by m7G modification. Furthermore, the interplay between m7G modification and its regulatory proteins can profoundly affect the efficacy of diverse adjuvant therapeutics, thereby potentially serving as a pivotal biomarker and therapeutic target for combinatory interventions in diverse cancer types.

Tumor Mutational Burden for Predicting Prognosis and Therapy Outcome of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), the primary hepatic malignancy, represents the second-highest cause of cancer-related death worldwide. Many efforts have been devoted to finding novel biomarkers for predicting both patients' survival and the outcome of pharmacological treatments, with a particular focus on immunotherapy. In this regard, recent studies have focused on unravelling the role of tumor mutational burden (TMB), i.e., the total number of mutations per coding area of a tumor genome, to ascertain whether it can be considered a reliable biomarker to be used either for the stratification of HCC patients in subgroups with different responsiveness to immunotherapy, or for the prediction of disease progression, particularly in relation to the different HCC etiologies. In this review, we summarize the recent advances on the study of TMB and TMB-related biomarkers in the HCC landscape, focusing on their feasibility as guides for therapy decisions and/or predictors of clinical outcome.

Comprehensive multi-omics analysis of the m7G in pan-cancer from the perspective of predictive, preventive, and personalized medicine

BACKGROUND

The N7-methylguanosine modification (m7G) of the 5' cap structure in the mRNA plays a crucial role in gene expression. However, the relation between m7G and tumor immune remains unclear. Hence, we intended to perform a pan-cancer analysis of m7G which can help explore the underlying mechanism and contribute to predictive, preventive, and personalized medicine (PPPM / 3PM).

METHODS

The gene expression, genetic variation, clinical information, methylation, and digital pathological section from 33 cancer types were downloaded from the TCGA database. Immunohistochemistry (IHC) was used to validate the expression of the m7G regulator genes (m7RGs) hub-gene. The m7G score was calculated by single-sample gene-set enrichment analysis. The association of m7RGs with copy number variation, clinical features, immune-related genes, TMB, MSI, and tumor immune dysfunction and exclusion (TIDE) was comprehensively assessed. CellProfiler was used to extract pathological section characteristics. XGBoost and random forest were used to construct the m7G score prediction model. Single-cell transcriptome sequencing (scRNA-seq) was used to assess the activation state of the m7G in the tumor microenvironment.

RESULTS

The m7RGs were highly expressed in tumors and most of the m7RGs are risk factors for prognosis. Moreover, the cellular pathway enrichment analysis suggested that m7G score was closely associated with invasion, cell cycle, DNA damage, and repair. In several cancers, m7G score was significantly negatively correlated with MSI and TMB and positively correlated with TIDE, suggesting an ICB marker potential. XGBoost-based pathomics model accurately predicts m7G scores with an area under the ROC curve (AUC) of 0.97. Analysis of scRNA-seq suggests that m7G differs significantly among cells of the tumor microenvironment. IHC confirmed high expression of EIF4E in breast cancer. The m7G prognostic model can accurately assess the prognosis of tumor patients with an AUC of 0.81, which was publicly hosted at https://pan-cancer-m7g.shinyapps.io/Panca-m7g/.

CONCLUSION

The current study explored for the first time the m7G in pan-cancer and identified m7G as an innovative marker in predicting clinical outcomes and immunotherapeutic efficacy, with the potential for deeper integration with PPPM. Combining m7G within the framework of PPPM will provide a unique opportunity for clinical intelligence and new approaches.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13167-022-00305-1.

Constructing and validating of m7G-related genes prognostic signature for hepatocellular carcinoma and immune infiltration: potential biomarkers for predicting the overall survival

Background

To investigate the prognostic significance of N7-methylguanosine (m7G) regulators and immune infiltration in liver hepatocellular carcinoma (LIHC).

Methods

The research measured predictive m7G genes in LIHC samples from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets. Data on the stemness index based on mRNA expression (mRNAsi), gene mutations, and corresponding clinical characteristics were obtained from TCGA and ICGC. Lasso regression was used to construct the prediction model to assess the m7G prognostic signals in LIHC. Based on these genes, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to identify key biological functions and pathways. The correlation between m7G RNA methylation regulators and the prognosis and immune infiltration of LIHC was evaluated.

Results

There were 21 m7G-related differentially expressed genes (DEGs) in LIHC and healthy tissues, and LIHC patients could be divided into two categories by consensus clustering of these DEGs. A five-gene predictive approach was employed using least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Patients in the low-risk group showed a significantly higher survival rate compared with those in the high-risk group (P=0.001). Validations using the ICGC database. Also, univariate and multivariate Cox regression analyses suggested that the risk score produced by the predictive model is an independent predictor for LIHC [hazard ratio (HR): 1.848, 95% confidence interval (CI): 1.286-2.656; HR: 2.597, 95% CI: 1.358-4.965]. The ROC curves of the ICGC cohort revealed that the five-gene prediction model performed well [area under the curve (AUC) =0.642 at 1 year, AUC =0.686 at 2 years, and AUC =0.667 at 3 years]. Immuno-oncology scoring revealed that in the high-risk group, among 16 immune cells, the expressions of neutrophils and natural killer (NK) cells were low and that of regulatory T-cells (Tregs) was high.

Conclusions

LIHC occurrence and progression are linked to m7G-related genes. Corresponding prognostic models help forecast the prognosis of LIHC patients. m7G-related genes and associated immune cell infiltration in the TME may serve as potential therapeutic targets in LIHC, which requires further trials. In addition, the m7G-related gene signature offers a viable alternative to predict LIHC, and these m7G-related genes show a prospective research area for LIHC targeted treatment in the future.