Single-cell sequencing and bulk RNA sequencing reveal a cell differentiation-related multigene panel to predict the prognosis and immunotherapy response of hepatocellular carcinoma.
Ontology highlight
ABSTRACT:
SUBMITTER: Liu J
PROVIDER: S-EPMC10106201 | biostudies-literature | 2023 Feb
REPOSITORIES: biostudies-literature
Publications
Single-cell sequencing and bulk RNA sequencing reveal a cell differentiation-related multigene panel to predict the prognosis and immunotherapy response of hepatocellular carcinoma.
Liu Jifeng J Yuan Qihang Q Ren Jie J Li Yiming Y Zhang Yunshu Y Shang Dong D
Chinese medical journal 20230220 4
Similar Datasets
Project description:Prostate cancer (PCa) stands as a prominent contributor to morbidity and mortality among males on a global scale. Cancer-associated fibroblasts (CAFs) are considered to be closely connected to tumour growth, invasion, and metastasis. We explored the role and characteristics of CAFs in PCa through bioinformatics analysis and built a CAFs-based risk model to predict prognostic treatment and treatment response in PCa patients. First, we downloaded the scRNA-seq data for PCa from the GEO. We extracted bulk RNA-seq data for PCa from the TCGA and GEO and adopted "ComBat" to remove batch effects. Then, we created a Seurat object for the scRNA-seq data using the package "Seurat" in R and identified CAF clusters based on the CAF-related genes (CAFRGs). Based on CAFRGs, a prognostic model was constructed by univariate Cox, LASSO, and multivariate Cox analyses. And the model was validated internally and externally by Kaplan-Meier analysis, respectively. We further performed GO and KEGG analyses of DEGs between risk groups. Besides, we investigated differences in somatic mutations between different risk groups. We explored differences in the immune microenvironment landscape and ICG expression levels in the different groups. Finally, we predicted the response to immunotherapy and the sensitivity of antitumour drugs between the different groups. We screened 4 CAF clusters and identified 463 CAFRGs in PCa scRNA-seq. We constructed a model containing 10 prognostic CAFRGs by univariate Cox, LASSO, and multivariate Cox analysis. Somatic mutation analysis revealed that TTN and TP53 were significantly more mutated in the high-risk group. Finally, we screened 31 chemotherapeutic drugs and targeted therapeutic drugs for PCa. In conclusion, we identified four clusters based on CAFs and constructed a new CAFs-based prognostic signature that could predict PCa patient prognosis and response to immunotherapy and might suggest meaningful clinical options for the treatment of PCa.
| S-EPMC10511553 | biostudies-literature
Project description:Gastric carcinoma (GC) heterogeneity represents a major barrier to accurate diagnosis and treatment. Here, we established a comprehensive single-cell transcriptional atlas to identify the cellular heterogeneity in malignant epithelial cells of GC using single-cell RNA sequencing (scRNA-seq). A total of 49,994 cells from nine patients with paired primary tumor and normal tissues were analyzed by multiple strategies. This study focused on the malignant epithelial cells, which were divided into three subtypes, including pit mucous cells, chief cells, and gastric and intestinal cells. The trajectory analysis results suggest that the differentiation of the three subtypes could be from the pit mucous cells to the chief cells and then to the gastric and intestinal cells. Lauren's histopathology of GC might originate from various subtypes of malignant epithelial cells. The functional enrichment analysis results show that the three subtypes focused on different biological processes (BP) and pathways related to tumor development. In addition, we generated and validated the prognostic signatures for predicting the OS in GC patients by combining the scRNA-seq and bulk RNA sequencing (bulk RNA-seq) datasets. Overall, our study provides a resource for understanding the heterogeneity of GC that will contribute to accurate diagnosis and prognosis.
| S-EPMC9407012 | biostudies-literature
Project description:Anoikis resistance, a notable factor in osteosarcoma, plays a significant role in tumor invasion and metastasis. This study seeks to identify a distinct gene signature that is specifically associated with the anoikis subcluster in osteosarcoma. Clinical, single-cell, and transcriptional data from TARGET and GEO datasets were used to develop a gene signature for osteosarcoma based on the anoikis subcluster. Univariate Cox and LASSO regression analyses were employed. The signature's predictive value was evaluated using time-dependent ROC and Kaplan-Meier analyses. Functional enrichment analyses and drug sensitivity analyses were conducted. Validation of three modular genes was performed using RT-qPCR and Western blotting. Signature (ZNF583, CGNL1, CXCL13) was developed to predict overall survival in osteosarcoma patients, targeting the anoikis subcluster. The signature demonstrated good performance in external validation. Stratification based on the signature revealed significantly different prognoses. The signature was an independent prognostic factor. The low-risk group showed enhanced immune cell infiltration and improved immune function. Drug sensitivity analysis indicated efficacy of chemotherapy agents. Prognostic nomograms incorporating the signature provided greater predictive accuracy and clinical utility. Signatures related to the anoikis subcluster play a significant role in osteosarcoma progression. Incorporating these findings into clinical decision-making can improve osteosarcoma treatment and patient outcomes.
| S-EPMC10657454 | biostudies-literature
Project description:Myeloid cells are physiologically related to innate immunity and inflammation. Tumor-associated myeloid cells gained increasing interest because of their critical roles in tumor progression and anticancer immune responses in human malignancies. However, the associations between tumor-associated myeloid cell-related genes and hepatocellular carcinoma have yet to be revealed. Here, through the integrating analysis of bulk and single-cell RNA (scRNA) sequencing of public HCC samples, we developed a gene signature to investigate the role of HCC-specific myeloid signature genes in HCC patients. We firstly defined 317 myeloid cell marker genes through analyzing scRNA data of HCC from the GEO dataset. After selecting the differentially expressed genes, eleven genes were also proved prognostic. Then we built a gene signature from the TCGA cohort and verified further with the ICGC dataset by applying the LASSO Cox method. An eight genes signature (FABP5, C15orf48, PABPC1, TUBA1B, AKR1C3, NQO1, AKR1B10, SPP1) was achieved finally. Patients in the high risk group correlated with higher tumor stages and poor survival than those in the low-risk group. The risk score was proved to be an independent risk factor for prognosis. The high risk group had higher infiltrations of dendritic cells, macrophages and Tregs. And the APC co-inhibition, T cell co-inhibition pathways were also activated. Besides, the risk score positively correlated with multidrug resistance proteins. In conclusion, our myeloid cell marker genes related signature can predict patients' survival and may also indicate the levels of immune infiltration and drug resistance.
| S-EPMC10027926 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) frequently recurs even after curative hepatectomy. To develop an integrated multigene expression panel, 144 patients were randomly assigned to either discovery or validation set in a 1:2 ratio. Using surgically resected HCC specimens, expression levels of 12 candidate molecular markers were determined using quantitative reverse-transcriptase PCR. In the discovery set, an expression panel was developed according to the concordance index (C-index) values for overall survival from all 4095 combinations of the 12 candidate molecular markers. Expression scores was determined with weighting according to the coefficient in a Cox regression, and patients were classified into grade 1, 2 and 3. Reproducibility was then tested in the validation set. A panel consisting of four markers, PRMT5, MAGED4, DPYSL3 and AJAP1 was selected as the optimal and most well-balanced set with a C-index value of 0.707. Patient prognosis was clearly stratified by the expression grade using this panel. In the validation set, both overall and disease-free survival rates decreased incrementally with as the grade increased. Higher grades were significantly associated with tumor multiplicity and vessel invasion. The prevalence of extrahepatic recurrences was increased in grade 3 patients. The integrated multigene expression panel clearly stratified HCC patients into low, intermediate and high risk.
| S-EPMC5642618 | biostudies-literature
Project description:Natural killer (NK) cells play essential roles in the tumor development, diagnosis, and prognosis of tumors. In this study, we aimed to establish a reliable signature based on marker genes in NK cells, thus providing a new perspective for assessing immunotherapy and the prognosis of patients with gastric cancer (GC). We analyzed a total of 1560 samples retrieved from the public database. We performed a comprehensive analysis of single-cell RNA-sequencing (scRNA-seq) data of gastric cancer and identified 377 marker genes for NK cells. By performing Cox regression analysis, we established a 12-gene NK cell-associated signature (NKCAS) for the Cancer Genome Atlas (TCGA) cohort, that assigned GC patients into a low-risk group (LRG) or a high-risk group (HRG). In the TCGA cohort, the areas under curve (AUC) value were 0.73, 0.81, and 0.80 at 1, 3, and 5 years. External validation of the predictive ability for the signature was then validated in the Gene Expression Omnibus (GEO) cohorts (GSE84437). The expression levels of signature genes were measured and validated in GC cell lines by real-time PCR. Moreover, NKCAS was identified as an independent prognostic factor by multivariate analysis. We combined this with a variety of clinicopathological characteristics (age, M stage, and tumor grade) to construct a nomogram to predict the survival outcomes of patients. Moreover, the LRG showed higher immune cell infiltration, especially CD8+ T cells and NK cells. The risk score was negatively associated with inflammatory activities. Importantly, analysis of the independent immunotherapy cohort showed that the LRG had a better prognosis and immunotherapy response when compared with the HRG. The identification of NK cell marker genes in this study suggests potential therapeutic targets. Additionally, the developed predictive signatures and nomograms may aid in the clinical management of GC.
| S-EPMC10985121 | biostudies-literature
Project description:The aging microenvironment serves important roles in cancers. However, most studies focus on circumscribed hot spots such as immunity and metabolism. Thus, it is well ignored that the aging microenvironment contributes to the proliferation of tumor. Herein, we established three prognosis-distinctive aging microenvironment subtypes, including AME1, AME2, and AME3, based on aging-related genes and characterized them with "Immune Exclusion," "Immune Infiltration," and "Immune Intermediate" features separately. AME2-subtype tumors were characterized by specific activation of immune cells and were most likely to be sensitive to immunotherapy. AME1-subtype tumors were characterized by inhibition of immune cells with high proportion of Catenin Beta 1 (CTNNB1) mutation, which was more likely to be insensitive to immunotherapy. Furthermore, we found that CTNNB1 may inhibit the expression of C-C Motif Chemokine Ligand 19 (CCL19), thus restraining immune cells and attenuating the sensitivity to immunotherapy. Finally, we also established a robust aging prognostic model to predict the prognosis of patients with hepatocellular carcinoma. Overall, this research promotes a comprehensive understanding about the aging microenvironment and immunity in hepatocellular carcinoma and may provide potential therapeutic targets for immunotherapy.
| S-EPMC9521435 | biostudies-literature
Project description:Background: Globally, hepatocellular carcinoma (HCC) is the sixth most frequent malignancy with a high incidence and a poor prognosis. Immune cell infiltration (ICI) underlies both the carcinogenesis and immunogenicity of tumors. However, a comprehensive classification system based on the immune features for HCC remains unknown. Methods: The HCC dataset from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts was used in this study. The ICI patterns of 571 patients were characterized using two algorithms: the patterns were determined based on the ICI using the ConsensusClusterPlus package, and principal component analysis (PCA) established the ICI scores. Differences in the immune landscape, biological function, and somatic mutations across ICI scores were evaluated and compared, followed by a predictive efficacy evaluation of ICI scores for immunotherapy by the two algorithms and validation using an external immunotherapy cohort. Results: Based on the ICI profile of the HCC patients, three ICI patterns were identified, including three subtypes having different immunological features. Individual ICI scores were determined; the high ICI score subtype was characterized by enhanced activation of immune-related signaling pathways and a significantly high tumor mutation burden (TMB); concomitantly, diminished immunocompetence and enrichment of pathways associated with cell cycle and RNA degradation were found in the low ICI score subtype. Taken together, our results contribute to a better understanding of an active tumor and plausible reasons for its poor prognosis. Conclusion: The present study reveals that ICI scores may serve as valid prognostic biomarkers for immunotherapy in HCC.
| S-EPMC8657761 | biostudies-literature
Project description:BackgroundIt was previously reported that the production of exerkines is positively associated with the beneficial effects of exercise in lung adenocarcinoma (LUAD) patients. This study proposes a novel scoring system based on muscle failure-related genes, to assist in clinical decision making.MethodsA comprehensive analysis of bulk and single cell RNA sequencing (scRNA-seq) of early, advanced and brain metastatic LUAD tissues and normal lung tissues was performed to identify muscle failure-related genes in LUAD and to determine the distribution of muscle failure-related genes in different cell populations. A novel scoring system, named MFI (Muscle failure index), was developed and validated. The differences in biological functions, immune infiltration, genomic alterations, and clinical significance of different subtypes were also investigated.ResultsFirst, we conducted single cell analysis on the dataset GSE131907 and identified eight cell subpopulations. We found that four muscle failure-related genes (BDNF, FNDC5, IL15, MSTN) were significantly increased in tumor cells. In addition, IL15 was widely distributed in the immune cell population. And we have validated it in our own clinical cohort. Then we created the MFI model based on 10 muscle failure-related genes using the LASSO algorithm, and MFI remained an independent prognostic factor of OS in both the training and validation cohorts. Moreover, we generated MFI in the single-cell dataset, in which cells with high MFI received and sent more signals compared to those with low MFI. Biological function analysis of both subtypes revealed stronger anti-tumor immune activity in the low MFI group, while tumor cells with high MFI had stronger metabolic and proliferative activity. Finally, we systematically assessed the immune cell activity and immunotherapy responses in LUAD patients, finding that the low MFI group was more sensitive to immunotherapy.ConclusionOverall, our study can improve the understanding of the role of muscle failure-related genes in tumorigenesis and we constructed a reliable MFI model for predicting prognosis and guiding future clinical decision making.
| S-EPMC9888242 | biostudies-literature
Project description:Colorectal cancer (CRC) is a prevalent malignant tumor characterized by high global incidence and mortality rates. Furthermore, it is imperative to comprehend the molecular mechanisms underlying its development and to identify effective prognostic markers. These efforts are crucial for pinpointing potential therapeutic targets and enhancing patient survival rates. Therefore, we develop a novel prognostic model aimed at providing new theoretical support for clinical prognosis evaluation and treatment. We downloaded data from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Subsequently, we performed single-cell analysis and developed a prognostic model associated with colorectal cancer. We divided the scRNA-seq dataset (GSE221575) into 19 cell clusters and classified these clusters into 11 distinct cell types using marker genes. Using univariate Cox regression and LASSO (Least Absolute Shrinkage and Selection Operator) analyses, we developed a prognostic model consisting of 9 genes. Based on our 9-gene model, we divided patients into high-risk and low-risk groups using the median risk score. The high-risk group demonstrated significant positive correlations with M0 macrophages, CD8+ T cells, and M2 macrophages. The enrichment analyses indicate significant enrichment of immune-related pathways in the high-risk group, including HEDGEHOG_SIGNALING, Wnt signaling pathway, and cell adhesion molecules. Drug sensitivity analysis revealed that the low-risk group was sensitive to 5 chemotherapeutic drugs, while the high-risk group was sensitive to only 1. Additionally, we developed a highly reliable nomogram for clinical application. This suggests that the risk score derived from our modeling analysis is highly effective for stratifying colorectal cancer samples. This study comprehensively applied bioinformatics methods to construct a risk score model. The model showed good predictive performance, offering potential guidance for individualized treatment of colorectal cancer patients. Furthermore, it may provide valuable insights into the disease's pathogenesis and identify potential therapeutic targets for further research.
| S-EPMC11889094 | biostudies-literature