Identification of TGF-β signaling-related molecular patterns, construction of a prognostic model, and prediction of immunotherapy response in gastric cancer

Identification of a fatty acid metabolism-related gene signature to predict prognosis in stomach adenocarcinoma

BACKGROUND

Fatty acid metabolism (FAM) contributes to tumorigenesis and tumor development, but the role of FAM in the progression of stomach adenocarcinoma (STAD) has not been comprehensively clarified.

METHODS

The expression data and clinical follow-up information were obtained from The Cancer Genome Atlas (TCGA). FAM pathway was analyzed by gene set enrichment analysis (GSEA) and single-sample GSEA (ssGSEA) methods. Univariate Cox regression analysis was conducted to select prognosis genes. Molecular subtypes were classified by consensus clustering analysis. Furthermore, least absolute shrinkage and selection operator (Lasso) analysis was employed to develop a risk model. ESTIMATE and tumour immune dysfunction and exclusion (TIDE) algorithm were used to assess immunity. pRRophetic package was conducted to predict drug sensitivity.

RESULTS

Based on 14 FAM related prognosis genes (FAMRG), 2 clusters were determined. Patients in C2 showed a worse overall survival (OS). Furthermore, a 7-FAMRG risk model was established as an independent predictor for STAD, with a higher riskscore indicating an unfavorable OS. High riskscore patients had higher TIDE score and these patients were more sensitive to anticancer drugs such as Bortezomib, Dasatinib and Pazopanib. A nomogram based on riskscore was an effective prediction tool applicable to clinical settings. The results from pan-cancer analysis supported a prominent application value of riskscore model in other cancer types.

CONCLUSION

The FAMRGs model established in this study could help predict STAD prognosis and offer new directions for future studies on dysfunctional FAM-induced damage and anti-tumor drugs in STAD disease.

Integration of single-cell RNA-seq and bulk RNA-seq to construct liver hepatocellular carcinoma stem cell signatures to explore their impact on patient prognosis and treatment

BACKGROUND

Liver hepatocellular carcinoma (LIHC) is a prevalent form of primary liver cancer. Research has demonstrated the contribution of tumor stem cells in facilitating tumor recurrence, metastasis, and treatment resistance. Despite this, there remains a lack of established cancer stem cells (CSCs)-associated genes signatures for effectively predicting the prognosis and guiding the treatment strategies for patients diagnosed with LIHC.

METHODS

The single-cell RNA sequencing (scRNA-seq) and bulk RNA transcriptome data were obtained based on public datasets and computerized firstly using CytoTRACE package and One Class Linear Regression (OCLR) algorithm to evaluate stemness level, respectively. Then, we explored the association of stemness indicators (CytoTRACE score and stemness index, mRNAsi) with survival outcomes and clinical characteristics by combining clinical information and survival analyses. Subsequently, weighted co-expression network analysis (WGCNA) and Cox were applied to assess mRNAsi-related genes in bulk LIHC data and construct a prognostic model for LIHC patients. Single-sample gene-set enrichment analysis (ssGSEA), Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) and Tumor Immune Estimation Resource (TIMER) analysis were employed for immune infiltration assessment. Finally, the potential immunotherapeutic response was predicted by the Tumor Immune Dysfunction and Exclusion (TIDE), and the tumor mutation burden (TMB). Additionally, pRRophetic package was applied to evaluate the sensitivity of high and low-risk groups to common chemotherapeutic drugs.

RESULTS

A total of four genes (including STIP1, H2AFZ, BRIX1, and TUBB) associated with stemness score (CytoTRACE score and mRNAsi) were identified and constructed a risk model that could predict prognosis in LIHC patients. It was observed that high stemness cells occurred predominantly in the late stages of LIHC and that poor overall survival in LIHC patients was also associated with high mRNAsi scores. In addition, pathway analysis confirmed the biological uniqueness of the two risk groups. Personalized treatment predictions suggest that patients with a low risk benefited more from immunotherapy, while those with a high risk group may be conducive to chemotherapeutic drugs.

CONCLUSION

The current study developed a novel prognostic risk signature with genes related to CSCs, which provides novel ideas for the diagnosis, prognosis and treatment of LIHC.

Downregulation of lncRNA EPB41L4A-AS1 promotes gastric cancer cell proliferation, migration and invasion

BACKGROUND

The incidence of gastric cancer ranks the first among digestive tract tumors in China. However, there are no specific symptoms in the early stage of the tumor and the diagnosis process is complex, so more effective detection methods are very needed. In this study, a novel long noncoding RNA (lncRNA) was introduced as a diagnostic biomarker for gastric cancer, which brought new thinking to the exploration of its pathological mechanism and clinical prediction.

METHODS

The level of lncRNA EPB41L4A-AS1 (EPB41L4A-AS1) in gastric cancer serum and cells was verified via real-time quantitative polymerase chain reaction (RT-qPCR). Receiver operating characteristic (ROC) curve was performed based on the EPB41L4A-AS1 level, and the diagnostic possibility of EPB41L4A-AS was analyzed. The chi-square test evaluated the correlation between EPB41L4A-AS expression and clinical information. The cells were cultured and transfected in vitro, and the mediations of abnormal EPB41L4A-AS level on the viability and motility of gastric cancer cells were verified through cell counting kit-8 (CCK-8) and Transwell assay. Furthermore, luciferase activity assay was performed to confirm the sponge molecule microRNA-17-5p (miR-17-5p) of EPB41L4A-AS1.

RESULTS

EPB41L4A-AS1 was decreased in gastric cancer, and low EPB41L4A-AS1 level indicated resultful diagnostic value. Overexpression of EPB41L4A-AS1 inhibited the activity of gastric cancer cells, while knockdown of EPB41L4A-AS1 promoted tumor deterioration. EPB41L4A-AS1 directly targeted and regulated the expression ofmiR-17-5p.

CONCLUSION

This study elaborated that EPB41L4A-AS1 is lowly expressed in gastric cancer. Silencing EPB41L4A-AS1 was beneficial to cell proliferation, migration, and invasion. EPB41L4A-AS1 provides a new possibility for the diagnosis of gastric cancer patients by targeting miR-17-5p.

Single-cell data revealed exhaustion of characteristic NK cell subpopulations and T cell subpopulations in hepatocellular carcinoma

BACKGROUND

The treatment and prognosis of patients with advanced hepatocellular carcinoma (HCC) have been a major medical challenge. Unraveling the landscape of tumor immune infiltrating cells (TIICs) in the immune microenvironment of HCC is of great significance to probe the molecular mechanisms.

METHODS

Based on single-cell data of HCC, the cell landscape was revealed from the perspective of TIICs. Special cell subpopulations were determined by the expression levels of marker genes. Differential expression analysis was conducted. The activity of each subpopulation was determined based on the highly expressed genes. CTLA4+ T-cell subpopulations affecting the prognosis of HCC were determined based on survival analysis. A single-cell regulatory network inference and clustering analysis was also performed to determine the transcription factor regulatory networks in the CTLA4+ T cell subpopulations.

RESULTS

10 cell types were identified and NK cells and T cells showed high abundance in tumor tissues. Two NK cells subpopulations were present, FGFBP2+ NK cells, B3GNT7+ NK cells. Four T cells subpopulations were present, LAG3+ T cells, CTLA4+ T cells, RCAN3+ T cells, and HPGDS+ Th2 cells. FGFBP2+ NK cells, and CTLA4+ T cells were the exhaustive subpopulation. High CTLA4+ T cells contributed to poor prognostic outcomes and promoted tumor progression. Finally, a network of transcription factors regulated by NR3C1, STAT1, and STAT3, which were activated, was present in CTLA4+ T cells.

CONCLUSION

CTLA4+ T cell subsets in HCC exhibited functional exhaustion characteristics that probably inhibited T cell function through a transcription factor network dominated by NR3C1, STAT1, and STAT3.

Integrated single-cell sequencing, spatial transcriptome sequencing and bulk RNA sequencing highlights the molecular characteristics of parthanatos in gastric cancer

BACKGROUND

Parthanatos is a novel programmatic form of cell death based on DNA damage and PARP-1 dependency. Nevertheless, its specific role in the context of gastric cancer (GC) remains uncertain.

METHODS

In this study, we integrated multi-omics algorithms to investigate the molecular characteristics of parthanatos in GC. A series of bioinformatics algorithms were utilized to explore clinical heterogeneity of GC and further predict the clinical outcomes.

RESULTS

Firstly, we conducted a comprehensive analysis of the omics features of parthanatos in various human tumors, including genomic mutations, transcriptome expression, and prognostic relevance. We successfully identified 7 cell types within the GC microenvironment: myeloid cell, epithelial cell, T cell, stromal cell, proliferative cell, B cell, and NK cell. When compared to adjacent non-tumor tissues, single-cell sequencing results from GC tissues revealed elevated scores for the parthanatos pathway across multiple cell types. Spatial transcriptomics, for the first time, unveiled the spatial distribution characteristics of parthanatos signaling. GC patients with different parthanatos signals often exhibited distinct immune microenvironment and metabolic reprogramming features, leading to different clinical outcomes. The integration of parthanatos signaling and clinical indicators enabled the creation of novel survival curves that accurately assess patients' survival times and statuses.

CONCLUSIONS

In this study, the molecular characteristics of parthanatos' unicellular and spatial transcriptomics in GC were revealed for the first time. Our model based on parthanatos signals can be used to distinguish individual heterogeneity and predict clinical outcomes in patients with GC.

Novel application of the ferroptosis-related genes risk model associated with disulfidptosis in hepatocellular carcinoma prognosis and immune infiltration

Hepatocellular carcinoma (HCC) stands as the prevailing manifestation of primary liver cancer and continues to pose a formidable challenge to human well-being and longevity, owing to its elevated incidence and mortality rates. Nevertheless, the quest for reliable predictive biomarkers for HCC remains ongoing. Recent research has demonstrated a close correlation between ferroptosis and disulfidptosis, two cellular processes, and cancer prognosis, suggesting their potential as predictive factors for HCC. In this study, we employed a combination of bioinformatics algorithms and machine learning techniques, leveraging RNA sequencing data, mutation profiles, and clinical data from HCC samples in The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC) databases, to develop a risk prognosis model based on genes associated with ferroptosis and disulfidptosis. We conducted an unsupervised clustering analysis, calculating a risk score (RS) to predict the prognosis of HCC using these genes. Clustering analysis revealed two distinct HCC clusters, each characterized by significantly different prognostic and immune features. The median RS stratified HCC samples in the TCGA, GEO, and ICGC cohorts into high-and low-risk groups. Importantly, RS emerged as an independent prognostic factor in all three cohorts, with the high-risk group demonstrating poorer prognosis and a more active immunosuppressive microenvironment. Additionally, the high-risk group exhibited higher expression levels of tumor mutation burden (TMB), immune checkpoints (ICs), and human leukocyte antigen (HLA), suggesting a heightened responsiveness to immunotherapy. A cancer stem cell infiltration analysis revealed a higher similarity between tumor cells and stem cells in the high-risk group. Furthermore, drug sensitivity analysis highlighted significant differences in response to antitumor drugs between the two risk groups. In summary, our risk prognostic model, constructed based on ferroptosis-related genes associated with disulfidptosis, effectively predicts HCC prognosis. These findings hold potential implications for patient stratification and clinical decision-making, offering valuable theoretical insights in this field.

The role of molecular subtypes and immune infiltration characteristics based on disulfidptosis-associated genes in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common type of lung cancer which accounts for about 40% of all lung cancers. Early detection, risk stratification and treatment are important for improving outcomes for LUAD. Recent studies have found that abnormal accumulation of cystine and other disulfide occurs in the cell under glucose starvation, which induces disulfide stress and increases the content of disulfide bond in actin cytoskeleton, resulting in cell death, which is defined as disulfidptosis. Because the study of disulfidptosis is in its infancy, its role in disease progression is still unclear. In this study, we detected the expression and mutation of disulfidptosis genes in LUAD using a public database. Clustering analysis based on disulfidptosis gene was performed and differential genes of disulfidptosis subtype were analyzed. 7 differential genes of disulfidptosis subtype were used to construct a prognostic risk model, and the causes of prognostic differences were investigated by immune-infiltration analysis, immune checkpoint analysis, and drug sensitivity analysis. qPCR was used to verify the expression of 7 key genes in lung cancer cell line (A549) and normal bronchial epithelial cell line (BEAS-2B). Since G6PD had the highest risk factor of lung cancer, we further verified the protein expression of G6PD in lung cancer cells by western blot, and confirmed through colony formation experiment that interference with G6PD was able to significantly inhibit the proliferation ability of lung cancer cells. Our results provide evidence for the role of disulfidptosis in LUAD and provide new ideas for individualized precision therapy of LUAD.

Characteristic of molecular subtypes based on PANoptosis-related genes and experimental verification of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a type of liver cancer that originates from liver cells. It is one of the most common types of liver cancer and a leading cause of cancer-related death worldwide. Early detection and treatment can improve the HCC prognosis. Therefore, it is necessary to further improve HCC markers and risk stratification. PANoptosome is a cytoplasmic polymer protein complex that regulates a proinflammatory programmed cell death pathway called "PANoptosis". The role of PANoptosis in HCC remains unclear. In this study, the molecular changes of PANoptosis related genes (PAN-RGs) in HCC were systematically evaluated. We characterized the heterogeneity of HCC by using consensus clustering to identify two distinct subtypes. The two subtypes showed different survival rate, biological function, chemotherapy drug sensitivity and immune microenvironment. After identification of PAN-RG differential expression genes (DEGs), a prognostic model was established by Cox regression analysis using minimum absolute contraction and selection operator (LASSO), and its prognostic value was verified by Cox regression analysis, Kaplan-Meier curve and receiver operating characteristic (ROC) curve. Our own specimens were also used to further validate the prognostic significance and possible clinical value of the selected targets. Subsequently, we conducted a preliminary discussion on the reasons for the influence of the model on the prognosis through TME analysis, drug resistance analysis, TMB analysis and other studies. This study provides a new idea for individualized and precise treatment of HCC.

Characteristic of molecular subtype based on lysosome-associated genes reveals clinical prognosis and immune infiltration of gastric cancer

Background

Lysosome are involved in nutrient sensing, cell signaling, cell death, immune responses and cell metabolism, which play an important role in the initiation and development of multiple tumors. However, the biological function of lysosome in gastric cancer (GC) has not been revealed. Here, we aim to screen lysosome-associated genes and established a corresponding prognostic risk signature for GC, then explore the role and underlying mechanisms.

Methods

The lysosome-associated genes (LYAGs) were obtained from MSigDB database. Differentially expressed lysosome-associated genes (DE-LYAGs) of GC were acquired based on the TCGA database and GEO database. According to expression profiles of DE-LYAGs, we divided the GC patients into different subgroups and then explored tumor microenvironment (TME) landscape and immunotherapy response in LYAG subtypes using GSVA, ESTIMATE and ssGSEA algorithms. Univariate Cox regression analysis, LASSO algorithm and multivariate Cox regression analysis were adopted to identify the prognostic LYAGs and then establish a risk model for patients with GC. The Kaplan-Meier analysis, Cox regression analysis and ROC analysis were utilized to evaluate the performance of the prognostic risk model. Clinical GC specimens were also used to verify the bioinformatics results by qRT-PCR assay.

Results

Thirteen DE-LYAGs were obtained and utilized to distinguish three subtypes in GC samples. Expression profiles of the 13 DE-LYAGs predicted prognosis, tumor-related immunological abnormalities and pathway dysregulation in these three subtypes. Furthermore, we constructed a prognostic risk model for GC based on DEG in the three subtypes. The Kaplan-Meier analysis suggested that higher risk score related to short OS rate. The Cox regression analysis and ROC analysis indicated that risk model had independent and excellent ability in predicting prognosis of GC patients. Mechanistically, a remarkable difference was observed in immune cell infiltration, immunotherapy response, somatic mutation landscape and drug sensitivity. qRT-PCR results showed that compared with corresponding adjacent normal tissues, most screened genes showed significant abnormal expressions and the expression change trends were consistent with the bioinformatics results.

Conclusions

We established a novel signature based on LYAGs which could be served as a prognostic biomarker for GC. Our study might provide new insights into individualized prognostication and precision treatment for GC.