A novel prognostic model based on single-cell RNA sequencing data for hepatocellular carcinoma

Identification of a novel FOXO3‑associated prognostic model in hepatocellular carcinoma

Although numerous molecular classifications are available to predict the prognosis of patients with hepatocellular carcinoma (HCC), they are still unsatisfactory. Forkhead box O3 (FOXO3) has been widely reported as a transcription factor involved in human cancers, but its role in HCC remains controversial. The present study aimed to explore the role of FOXO3 in HCC, as well as to identify biomarkers and construct prognostic models based on FOXO3. FOXO3 was highly expressed in HCC and was closely associated with poor prognosis in The Cancer Genome Atlas (the training set) and International Cancer Genome Consortium (the validation set). Subsequently, a co-expression network indicated that the red modules were closely related to FOXO3. Five key FOXO3-related genes [DEAD-box helicase 55 (DDX55), RAB10, member RAS oncogene family (RAB10), RAB7A, TATA-box binding protein associated factor, RNA polymerase I subunit B (TAF1B) and TAF3] were obtained using Cox-least absolute shrinkage and selection operator analyses. The 5-gene signature successfully predicted the prognosis of patients with HCC in both the training and validation sets. Enrichment analysis suggested marked differences in AKT and cell cycle-related (E2F targets and G2/M checkpoints) pathways between HCC subgroups. Furthermore, the tumor microenvironment analysis suggested that the difference in the distribution of M2 macrophages among various subgroups may contribute to the poor prognosis using the CIBERSORTx framework. Furthermore, the mRNA and protein expressions of DDX55, RAB10, RAB7A, TAF1B and TAF3 were found to be higher in HCC tissues compared with paracancerous tissues using RT-qPCR and western blotting. Additionally, knockdown of RAB10, RAB7A and TAF3 inhibited proliferation of Huh7 cells, assessed by a Cell Counting Kit-8 assay. In conclusion, a novel FOXO3-related model was constructed and revealed that RAB10, RAB7A and TAF3 may be potential molecular targets or biomarkers for HCC.

Multi-omics analyses develop and validate the optimal prognostic model on overall survival prediction for resectable hepatocellular carcinoma

Background

Prediction of prognosis in patients with hepatocellular carcinoma (HCC) by single-omics profiling has been widely studied. However, the prognosis related to biomarkers of multiple omics has not been investigated. We aimed to establish and validate a prediction model for prognosis prediction of resectable HCC combining multi-omics and clinicopathological factors.

Methods

The training cohort involved multi-omics data of 330 patients with resectable HCC (stage I-IIIA) at mutational, copy number variation (CNV), transcriptional, and methylation levels from The Cancer Genome Atlas (TCGA) database, along with clinicopathological information. The validation cohort involved samples from 40 HCC patients of Beijing Youan Hospital. Univariate and multivariate analyses were performed in single-omics with clinicopathological variables regarding patient prognosis, and independent risk factors were combined to establish the multi-omics model. The predictive accuracy was assessed by the receiver operating characteristic (ROC) method.

Results

The mutational, copy number, transcriptional, and methylation alterations in HCC were characterized. TP53, CTNNB1, and TTN were among the genes with the top mutational frequency, and FBN1 and MAP1B mutations were independent risk factors for patient overall survival (OS). 1q21.3 and 1q23.3 ranked the highest in copy number amplifications, and 8p12 and 8p23.3 ranked the highest in deletions, and CSMD1, TP53, and RB1 were genes with the most frequent CNVs. AFP, GPC3, and TERT were among genes with the most significant aberrant transcription, and the transcription of CCNJL, FRMD1, and GRPEL2 were independent risk factors for OS. Both hypermethylation and hypomethylation can be observed. The aberrant methylation of CXorf15, DACT2, GP6, KIAA1522, and PDIA3 were independent risk factors. Single-omics models were established with independent risk factors, and were validated by internal and external datasets. A prognostic model for OS with multi-omics independent risk factors and clinicopathlogical information was established. Internal and external validation achieved an optimal maximal area under the curve (AUC) of 0.98 at 1 year and 0.88 at 2 years, respectively.

Conclusions

A multi-omics model combining molecular aberrancies and clinicopathological information was established and proved to be optimal for prognosis prediction of resectable HCC. This model may be helpful for therapeutic strategy selection and survival assessment.

Multiple machine learning algorithms identified SLC6A8 as a diagnostic biomarker of the late stage of Hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a chronic liver disease characterized by persistent tumor growth, contributing significantly to mortality rates worldwide. Consequently, there is an urgent need to develop effective diagnostic and treatment strategies for HCC. This study aims to identify crucial genes for early HCC diagnosis to mitigate disease progression and to investigate differences in immune cell infiltration between early-stage and late-stage HCC. We integrated two published datasets for a comprehensive analysis, identifying 575 DEGs subjected to GSEA to reveal pathways distinguishing early-stage from late-stage HCC. Notably, the gene SLC6A8 emerged as a potential diagnostic biomarker for late-stage HCC through machine learning (LASSO-LR/SVM-RFE/RF-Boruta). ROC curves for SLC6A8 were utilized to evaluate diagnostic accuracy. The ImmuCellAI algorithm assessed immune cell composition differences between early and late-stage HCC, revealing that SLC6A8 expression positively correlates with resting Tfh cells and Th2, while negatively correlating with B cells, indicating its association with immune cell infiltration patterns. To strengthen our results, we further analyzed SLC6A8 expression using single-cell transcriptome data, confirming notably overexpression in late-stage HCC, particularly in key liver cell types such as Hepatocyte cells. Overall, our study nominates SLC6A8 as a dual biomarker for HCC Staging and precision therapy.

Exploring the tumor microenvironment of breast cancer to develop a prognostic model and predict immunotherapy responses

Breast cancer is the most prevalent malignancy in women and exhibits significant heterogeneity. The tumor microenvironment (TME) plays a critical role in tumorigenesis, progression, and response to therapy. However, its impact on the prognosis and immunotherapy responses is incompletely understood. Using public databases, we conducted a comprehensive investigation of transcriptome and single-cell sequencing data. After performing immune infiltration analysis, we conducted consensus clustering, weighted gene co-expression network analysis (WGCNA), Cox regression, and least absolute shrinkage and selection operator (Lasso) regression to identify independent prognostic genes in breast cancer. Subsequently, we developed a prognostic model for patients with breast cancer. Tumor Immune Dysfunction and Exclusion (TIDE) values were used to assess patient's responsiveness to breast cancer. Based on single-cell RNA-sequencing data, we identified various cell types through cluster analysis and investigated the expression of prognostic model genes in each cell type. The drug sensitivity of targeted therapeutic agents for breast cancer treatment was analyzed in different cell types. We identified 12 independent prognostic genes associated with breast cancer and used these genes to construct a prognostic model. The prognostic model accurately discriminated between patients classified as high- and low-risk, providing precise prognostic predictions for individual patients. Additionally, our model exhibited a robust capacity to predict the immunotherapeutic response in breast cancer patients. Our investigation revealed a notable association between the proportion of endothelial cells (ECs) and patient prognosis in breast cancer. A prognostic model for breast cancer was formulated that showed close associations between prognosis and response to immunotherapy. For patients predicted by our model to not respond effectively to immunotherapeutic agents, it may be considered to combine immunotherapeutic agents with targeted therapeutic agents identified through our drug sensitivity analysis, which could potentially enhance treatment efficacy.

The role of hypoxia-senescence co-related molecular subtypes and prognostic characteristics in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is known for its high invasiveness, high fatality rate. Both hypoxia and senescence play crucial roles in the initiation and progression of cancer, yet their prognostic implications in HCC are yet to be fully understood. The hypoxia-senescence co-related genes (HSCRGs) were screened from public databases. Transcriptome data and clinical information were obtained from patients with HCC using the Cancer Genome Atlas, GSE76427, and International Cancer Genome Consortium (ICGC). The random forest tree algorithm was used to identify the characteristic genes of the disease, and the genes were verified by related experiments. SVM algorithm was used to classify HCC patients based on HSCRGs. The prediction model based on HSCRGs was established by LASSO, univariate and multivariate COX regression analysis. We used the ICGC for outside validation. The risk score model was analyzed from subgroup analysis, immune infiltration, and functional strength. The expression patterns of key prognostic genes in tumor microenvironment were decoded by single cell analysis. A total of 184 HSCRGs were identified. The expression pattern and functional characteristics of MLH1 gene in HCC were verified. Two HCC subtypes were identified based on HSCRGs. Then, a prediction model based on HSCRGs was established, and risk score was identified as an independent prognostic indicator of HCC. A new nomogram is constructed and shows good prediction ability. We further determined that the level of infiltration of immune cells and the expression of immune checkpoints are significantly affected by the risk score. The immune microenvironment was different between the two risk groups. The high-risk group was dominated by immunosuppressed cells, and the prognosis was poor. Single-cell analysis revealed the expression of seven key prognostic genes in the tumor microenvironment. Finally, qPCR results further verified the expression levels of seven prognostic genes. HSCRGs are of great significance in the prognosis prediction, risk stratification and targeted therapy of patients with HCC.

Integrated analysis of single-cell and bulk RNA-sequencing to predict prognosis and therapeutic response for colorectal cancer

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.

Establishment of a Lactylation-Related Gene Signature for Hepatocellular Carcinoma Applying Bulk and Single-Cell RNA Sequencing Analysis

Background: Lactylation is closely involved in cancer progression, but its role in hepatocellular carcinoma (HCC) is unclear. The present work set out to develop a lactylation-related gene (LRG) signature for HCC. Methods: The lactylation score of tumor and normal groups was calculated using the gene set variation analysis (GSVA) package. The single-cell RNA sequencing (scRNA-seq) analysis of HCC was performed in the "Seurat" package. Prognostic LRGs were selected by performing univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analyses to develop and validate a Riskscore model. Functional enrichment analysis was conducted by gene set enrichment analysis (GSEA) using the "clusterProfiler" package. Genomic characteristics between different risk groups were compared, and tumor mutational burden (TMB) was calculated by the "Maftools" package. Immune cell infiltration was assessed by algorithms of cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT), microenvironment cell populations-counter (MCP-counter), estimating the proportions of immune and cancer cells (EPIC), tumor immune estimation resource (TIMER), and single-sample gene set enrichment analysis (ssGSEA). Immunotherapy response was predicted by the tumor immune dysfunction and exclusion (TIDE) algorithm. Drug sensitivity was analyzed using the "pRRophetic" package. A nomogram was established using the "rms" package. The expressions of the prognostic LRGs in HCC cells were verified by in vitro test, and cell counting kit-8 (CCK-8), wound healing, and transwell assays were carried out to measure the viability, migration, and invasion of HCC cells. Results: The lactylation score, which was higher in the tumor group than in the normal group, has been confirmed as an independent factor for the prognostic evaluation in HCC. Six prognostic LRGs, including two protective genes (FTCD and APCS) and four risk genes (LGALS3, C1orf43, TALDO1, and CCT5), were identified to develop a Riskscore model with a strong prognostic prediction performance in HCC. The scRNA-seq analysis revealed that LGALS3 was largely expressed in myeloid cells, while APCS, FTCD, TALDO1, CCT5, and C1orf43 were mainly expressed in hepatocytes. The high-risk group was primarily enriched in the pathways involved in tumor occurrence and development, with higher T cell infiltration. Moreover, the high-risk group was found to be less responsive to immunotherapy but was more sensitive to chemotherapeutic drugs. By integrating Riskscore and clinical features, a nomogram with a high predictive accuracy was developed. Additionally, C1orf43, CCT5, TALDO1, and LGALS3 were highly expressed in HCC cells. Silencing CCT5 inhibited the viability, migration, and invasion of HCC cells. Conclusion: The present work developed a novel LRG gene signature that could be considered a promising therapeutic target and biomarker for HCC.

Specific BCG-related gene expression levels correlate with immune cell infiltration and prognosis in melanoma

Melanoma, caused by malignant melanocytes, is known for its invasiveness and poor prognosis. Therapies are often ineffective due to their heterogeneity and resistance. Bacillus Calmette-Guérin (BCG), primarily a tuberculosis vaccine, shows potential in treating melanoma by activating immune responses. In this study, data from The Cancer Genome Atlas and the National Center for Biotechnology Information Gene Expression Omnibus database were utilized to determine pivotal DEGs such as DSC2, CXCR1, BOK, and CSTB, which are significantly upregulated in BCG-treated blood samples and are strongly associated with the prognosis of melanoma. We employ tools like edgeR and ggplot2 for functional and pathway analysis and develop a prognostic model using LASSO Cox regression analysis to predict patient survival. A notable finding is the correlation between BCG-related genes and immune cell infiltration in melanoma, highlighting the potential of these genes as both biomarkers and therapeutic targets. Additionally, the study examines genetic alterations in these genes and their impact on the disease. This study highlights the necessity of further exploring BCG-related genes for insights into melanoma pathogenesis and treatment enhancement, suggesting that BCG's role in immune activation could offer novel therapeutic avenues in cancer treatment.

Integrative multi-omics analysis reveals a novel subtype of hepatocellular carcinoma with biological and clinical relevance

Background

Hepatocellular carcinoma (HCC) is a highly heterogeneous tumor, and the development of accurate predictive models for prognosis and drug sensitivity remains challenging.

Methods

We integrated laboratory data and public cohorts to conduct a multi-omics analysis of HCC, which included bulk RNA sequencing, proteomic analysis, single-cell RNA sequencing (scRNA-seq), spatial transcriptomics sequencing (ST-seq), and genome sequencing. We constructed a tumor purity (TP) and tumor microenvironment (TME) prognostic risk model. Proteomic analysis validated the TP-TME-related signatures. Joint analysis of scRNA-seq and ST-seq revealed characteristic clusters associated with TP high-risk subtypes, and immunohistochemistry confirmed the expression of key genes. We conducted functional enrichment analysis, transcription factor activity inference, cell-cell interaction, drug efficacy analysis, and mutation information analysis to identify a novel subtype of HCC.

Results

Our analyses constructed a robust HCC prognostic risk prediction model. The patients with TP-TME high-risk subtypes predominantly exhibit hypoxia and activation of the Wnt/beta-catenin, Notch, and TGF-beta signaling pathways. Furthermore, we identified a novel subtype, XPO1+Epithelial. This subtype expresses signatures of the TP risk subtype and aligns with the biological behavior of high-risk patients. Additional analyses revealed that XPO1+Epithelial is influenced primarily by fibroblasts via ligand-receptor interactions, such as FN1-(ITGAV+ITGB1), and constitute a significant component of the TP-TME subtype. Moreover, XPO1+Epithelial interact with monocytes/macrophages, T/NK cells, and endothelial cells through ligand-receptor pairs, including MIF-(CD74+CXCR4), MIF-(CD74+CD44), and VEGFA-VEGFR1R2, respectively, thereby promoting the recruitment of immune-suppressive cells and angiogenesis. The ST-seq cohort treated with Tyrosine Kinase Inhibitors (TKIs) and Programmed Cell Death Protein 1 (PD-1) presented elevated levels of TP and TME risk subtype signature genes, as well as XPO1+Epithelial, T-cell, and endothelial cell infiltration in the treatment response group. Drug sensitivity analyses indicated that TP-TME high-risk subtypes, including sorafenib and pembrolizumab, were associated with sensitivity to multiple drugs. Further exploratory analyses revealed that CTLA4, PDCD1, and the cancer antigens MSLN, MUC1, EPCAM, and PROM1 presented significantly increase expression levels in the high-risk subtype group.

Conclusions

This study constructed a robust prognostic model for HCC and identified novel subgroups at the single-cell level, potentially assisting in the assessment of prognostic risk for HCC patients and facilitating personalized drug therapy.

Unraveling the Heterogeneity of Tumor Immune Microenvironment in Hepatocellular Carcinoma by SingleCell RNA Sequencing and its Implications for Prognosis and Therapeutic Response

Tumor immune microenvironment (TIME) has become a new hotspot in cancer research over the past few years. Tumor immune microenvironment of hepatocellular carcinoma (HCC) is especially intriguing as HCC is reported to be highly heterogeneous by previous genomic and cytological studies. It is also closely related to patient prognosis and therapeutic outcome. The recently emerged single-cell RNA sequencing (scRNAseq) technique provides a new tool for TIME study, and current studies have made great advances in defining the roles of TIME in HCC pathogenesis and therapy. Current evidence suggests that heterogeneity is a key player influencing therapeutic response, drug resistance, and prognosis. However, our understanding is limited on the roles of TIME heterogeneity in HCC development, prognosis, and therapeutic response, especially in the era of immunotherapy. This review aims to unravel the heterogeneity of TIME in HCC by scRNAseq, with specific focuses on the cellular, transcriptional, and marker perspectives of TIME heterogeneity in HCC, as well as assessing prognostic and therapeutic response by heterogeneity markers. By summarizing current discoveries regarding TIME heterogeneity, we hope to provide clues on the crucial roles of various cellular components in the development and progression of HCC. We also hope to identify potential markers and therapeutic targets for prognosis assessment and personalized treatment to improve patient outcomes. Combined therapies from multiple dimensions regarding heterogeneity may provide new opportunities to treat HCC more effectively.

Identification of Cuproptosis-Associated Prognostic Gene Expression Signatures from 20 Tumor Types

We investigated the mRNA expression of 124 cuproptosis-associated genes in 7489 biopsies from 20 different tumor types of The Cancer Genome Atlas (TCGA). The KM plotter algorithm has been used to calculate Kaplan-Meier statistics and false discovery rate (FDR) corrections. Interaction networks have been generated using Ingenuity Pathway Analysis (IPA). High mRNA expression of 63 out of 124 genes significantly correlated with shorter survival times of cancer patients across all 20 tumor types. IPA analyses revealed that their gene products were interconnected in canonical pathways (e.g., cancer, cell death, cell cycle, cell signaling). Four tumor entities showed a higher accumulation of genes than the other cancer types, i.e., renal clear cell carcinoma (n = 21), renal papillary carcinoma (n = 13), kidney hepatocellular carcinoma (n = 13), and lung adenocarcinoma (n = 9). These gene clusters may serve as prognostic signatures for patient survival. These signatures were also of prognostic value for tumors with high mutational rates and neoantigen loads. Cuproptosis is of prognostic significance for the survival of cancer patients. The identification of specific gene signatures deserves further exploration for their clinical utility in routine diagnostics.

Mechanisms associated with cuproptosis and implications for ovarian cancer

Ovarian cancer, a profoundly fatal gynecologic neoplasm, exerts a substantial economic strain on nations globally. The formidable challenge of its frequent relapse necessitates the exploration of novel cytotoxic agents, efficacious antineoplastic medications with minimal adverse effects, and strategies to surmount resistance to primary chemotherapeutic agents. These endeavors aim to supplement extant pharmacological interventions and elucidate molecular mechanisms underlying induced cytotoxicity, distinct from conventional therapeutic modalities. Recent scientific research has unveiled a novel form of cellular demise, known as copper-death, which is contingent upon the intracellular concentration of copper. Diverging from conventional mechanisms of cellular demise, copper-death exhibits a pronounced reliance on mitochondrial respiration, particularly the tricarboxylic acid (TCA) cycle. Tumor cells manifest distinctive metabolic profiles and elevated copper levels in comparison to their normal counterparts. The advent of copper-death presents alluring possibilities for targeted therapeutic interventions within the realm of cancer treatment. Hence, the primary objective of this review is to present an overview of the proteins and intricate mechanisms associated with copper-induced cell death, while providing a comprehensive summary of the knowledge acquired regarding potential therapeutic approaches for ovarian cancer. These findings will serve as valuable references to facilitate the advancement of customized therapeutic interventions for ovarian cancer.

The cuproptosis-related signature predicts the prognosis and immune microenvironments of primary diffuse gliomas: a comprehensive analysis

BACKGROUND

Evidence has revealed a connection between cuproptosis and the inhibition of tumor angiogenesis. While the efficacy of a model based on cuproptosis-related genes (CRGs) in predicting the prognosis of peripheral organ tumors has been demonstrated, the impact of CRGs on the prognosis and the immunological landscape of gliomas remains unexplored.

METHODS

We screened CRGs to construct a novel scoring tool and developed a prognostic model for gliomas within the various cohorts. Afterward, a comprehensive exploration of the relationship between the CRG risk signature and the immunological landscape of gliomas was undertaken from multiple perspectives.

RESULTS

Five genes (NLRP3, ATP7B, SLC31A1, FDX1, and GCSH) were identified to build a CRG scoring system. The nomogram, based on CRG risk and other signatures, demonstrated a superior predictive performance (AUC of 0.89, 0.92, and 0.93 at 1, 2, and 3 years, respectively) in the training cohort. Furthermore, the CRG score was closely associated with various aspects of the immune landscape in gliomas, including immune cell infiltration, tumor mutations, tumor immune dysfunction and exclusion, immune checkpoints, cytotoxic T lymphocyte and immune exhaustion-related markers, as well as cancer signaling pathway biomarkers and cytokines.

CONCLUSION

The CRG risk signature may serve as a robust biomarker for predicting the prognosis and the potential viability of immunotherapy responses. Moreover, the key candidate CRGs might be promising targets to explore the underlying biological background and novel therapeutic interventions in gliomas.

Exploring the diagnostic value of endothelial cell and angiogenesis-related genes in Hashimoto's thyroiditis based on transcriptomics and single cell RNA sequencing

(1) BACKGROUND: Hashimoto's thyroiditis (HT) can cause angiogenesis in the thyroid gland. However, the molecular mechanism of endothelial cells and angiogenesis related genes (ARGs) has not been extensively studied in HT. (2) METHODS: The HRA001684, GSE29315 and GSE163203 datasets were included in this study. Using single-cell analysis, weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, machine learning algorithms and expression analysis for exploration. And receiver operator characteristic (ROC) curves was draw. Gene set enrichment analysis (GSEA) was utilized to investigate the biological function of the biomarkers. Meanwhile, we investigated into the relationship between biomarkers and different types of immune cells. Additionally, the expression of biomarkers in the TCGA-TC dataset was examined and the mRNA-drug interaction network was constructed. (3) RESULTS: We found 14 cell subtypes were obtained in HT samples after single-cell analysis. A total of 5 biomarkers (CD52, CD74, CD79A, HLA-B and RGS1) were derived, and they had excellent diagnostic performance. Then, 27 drugs targeting biomarkers were predicted. The expression analysis showed that CD74 and HLA-B were significantly up-regulated in HT samples. (4) CONCLUSION: In this study, 5 biomarkers (CD52, CD74, CD79A, HLA-B and RGS1) were screened and their expressions in endothelial cells was compared to offer a new reference for the recognition and management of HT.

Application of artificial intelligence in the prediction of immunotherapy efficacy in hepatocellular carcinoma: Current status and prospects

Artificial Intelligence (AI) has increased as a potent tool in medicine, with promising oncology applications. The emergence of immunotherapy has transformed the treatment terrain for hepatocellular carcinoma (HCC), offering new hope to patients with this challenging malignancy. This article examines the role and future of AI in forecasting the effectiveness of immunotherapy in HCC. We highlight the potential of AI to revolutionize the prediction of therapy response, thus improving patient selection and clinical outcomes. The article further outlines the challenges and future research directions in this emerging field.

Bulk and single-cell transcriptome reveal the immuno-prognostic subtypes and tumour microenvironment heterogeneity in HCC

BACKGROUND & AIMS

Accumulating evidences suggest tumour microenvironment (TME) profoundly influence clinical outcome in hepatocellular carcinoma (HCC). Existing immune subtypes are susceptible to batch effects, and integrative analysis of bulk and single-cell transcriptome is helpful to recognize immune subtypes and TME in HCC.

METHODS

Based on the relative expression ordering (REO) of 1259 immune-related genes, an immuno-prognostic signature was developed and validated in 907 HCC samples from five bulk transcriptomic cohorts, including 72 in-house samples. The machine learning models based on subtype-specific gene pairs with stable REOs were constructed to jointly predict immuno-prognostic subtypes in single-cell RNA-seq data and validated in another single-cell data. Then, cancer characteristics, immune landscape, underlying mechanism and therapeutic benefits between subtypes were analysed.

RESULTS

An immune-related signature with 29 gene pairs stratified HCC samples individually into two risk subgroups (C1 and C2), which was an independent prognostic factor for overall survival. The machine learning models verified the immune subtypes from five bulk cohorts to two single-cell transcriptomic data. Integrative analysis revealed that C1 had poorer outcomes, higher CNV burden and malignant scores, higher sensitivity to sorafenib, and exhibited an immunosuppressive phenotype with more regulators, e.g., myeloid-derived suppressor cells (MDSCs), Mø_SPP1, while C2 was characterized with better outcomes, higher metabolism, more benefit from immunotherapy, and displayed active immune with more effectors, e.g., tumour infiltrating lymphocyte and dendritic cell. Moreover, both two single-cell data revealed the crosstalk of SPP1-related L-R pairs between cancer and immune cells, especially SPP1-CD44, might lead to immunosuppression in C1.

CONCLUSIONS

The REO-based immuno-prognostic subtypes were conducive to individualized prognosis prediction and treatment options for HCC. This study paved the way for understanding TME heterogeneity between immuno-prognostic subtypes of HCC.

Prognostic prediction using a gene signature developed based on exhausted T cells for liver cancer patients

Background

Liver hepatocellular carcinoma (LIHC) is a solid primary malignancy with poor prognosis. This study discovered key prognostic genes based on T cell exhaustion and used them to develop a prognostic prediction model for LIHC.

Methods

SingleR's annotations combined with Seurat was used to automatically annotate the single-cell clustering results of the LIHC dataset GSE166635 downloaded from the Gene Expression Omnibus (GEO) database and to identify clusters related to exhausted T cells. Patients were classified using ConsensusClusterPlus package. Next, weighted gene co-expression network analysis (WGCNA) package was employed to distinguish key gene module, based on which least absolute shrinkage and selection operator (Lasso) and multi/univariate cox analysis were performed to construct a RiskScore system. Kaplan-Meier (KM) analysis and receiver operating characteristic curve (ROC) were employed to evaluate the efficacy of the model. To further optimize the risk model, a nomogram capable of predicting immune infiltration and immunotherapy sensitivity in different risk groups was developed. Expressions of genes were measured by quantitative real-time polymerase chain reaction (qRT-PCR), and immunofluorescence and Cell Counting Kit-8 (CCK-8) were performed for analyzing cell functions.

Results

We obtained 18,413 cells and clustered them into 7 immune and non-immune cell subpopulations. Based on highly variable genes among T cell exhaustion clusters, 3 molecular subtypes (C1, C2 and C3) of LIHC were defined, with C3 subtype showing the highest score of exhausted T cells and a poor prognosis. The Lasso and multivariate cox analysis selected 7 risk genes from the green module, which were closely associated with the C3 subtype. All the patients were divided into low- and high-risk groups based on the medium value of RiskScore, and we found that high-risk patients had higher immune infiltration and immune escape and poorer prognosis. The nomogram exhibited a strong performance for predicting long-term LIHC prognosis. In vitro experiments revealed that the 7 risk genes all had a higher expression in HCC cells, and that both liver HCC cell numbers and cell viability were reduced by knocking down MMP-9.

Conclusion

We developed a RiskScore model for predicting LIHC prognosis based on the scRNA-seq and RNA-seq data. The RiskScore as an independent prognostic factor could improve the clinical treatment for LIHC patients.

SIMarker: Cellular similarity detection and its application to diagnosis and prognosis of liver cancer

BACKGROUND

The emergence of single-cell technology offers a unique opportunity to explore cellular similarity and heterogeneity between precancerous diseases and solid tumors. However, there is lacking a systematic study for identifying and characterizing similarities at single-cell resolution.

METHODS

We developed SIMarker, a computational framework to detect cellular similarities between precancerous diseases and solid tumors based on gene expression at single-cell resolution. Taking hepatocellular carcinoma (HCC) as a case study, we quantified the cellular and molecular connections between HCC and cirrhosis. Core analysis modules of SIMarker is publicly available at https://github.com/xmuhuanglab/SIMarker ("SIM" means "similarity" and "Marker" means "biomarkers).

RESULTS

We found PGA5+ hepatocytes in HCC showed cirrhosis-like characteristics, including similar transcriptional programs and gene regulatory networks. Consequently, the genes constituting the gene expression program of these cirrhosis-like subpopulations were designated as cirrhosis-like signatures (CLS). Strikingly, our utilization of CLS enabled the development of diagnosis and prognosis biomarkers based on within-sample relative expression orderings of gene pairs. These biomarkers achieved high precision and concordance compared with previous studies.

CONCLUSIONS

Our work provides a systematic method to investigate the clinical translational significance of cellular similarities between HCC and cirrhosis, which opens avenues for identifying similar paradigms in other categories of cancers and diseases.

Clathrin light chain A facilitates small extracellular vesicle uptake to promote hepatocellular carcinoma progression

BACKGROUND

Endocytosis is a fundamental process for internalizing small extracellular vesicles (sEVs). The present study aimed to elucidate the role of clathrin light chain A (CLTA) in sEV uptake in hepatocellular carcinoma (HCC).

MATERIALS AND METHODS

CLTA expression was analyzed by bioinformatics, quantitative PCR and immunohistochemistry. The clinical relevance of CLTA was analyzed by Fisher's exact test, Kaplan-Meier analysis, and multivariate cox regression model. The functions of CLTA in sEV uptake and cancerous properties were examined by PKH67-sEV uptake, MTT, colony formation, and transwell assays. Mass spectrometry was used to identify the downstream effectors of CLTA. CLTA inhibitor, Pitstop 2, was tested in a mouse model of patient-derived xenografts (PDXs).

RESULTS

CLTA expression was higher in tumor tissues than in non-tumorous liver tissues and progressively increased from the early to late tumor stage. CLTA overexpression was associated with larger tumor size and poor prognosis in HCC. Cellular CLTA contributed to the sEV uptake, resulting in enhanced cancerous properties. Mechanistically, CLTA increases capping actin protein gelsolin-like (CAPG) expression to facilitate sEV uptake, thereby promoting the proliferation, motility, and invasiveness of HCC cells. What's more, the CLTA inhibitor Pitstop 2 alone or in combination with sorafenib attenuated tumor growth in mice implanted with PDXs.

CONCLUSIONS

The study reveals the role of CLTA in sEV uptake to promote HCC progression. Inhibition of CLTA and its mediated pathway illuminate a new therapeutic strategy for HCC patients.

Multi-omics analysis defines a cuproptosis-related prognostic model for ovarian cancer: Implication of WASF2 in cuproptosis resistance

BACKGROUND

Ovarian cancer (OVC) is one of the deadliest and most aggressive tumors in women, with an increasing incidence in recent years. Cuproptosis, a newly discovered type of programmed cell death, is caused by intracellular copper-mediated lipoylated protein aggregation and proteotoxic stress. However, the role of cuproptosis-related features in OVC remains elusive.

METHODS

The single-cell sequencing data from GSE154600 and bulk transcriptome data of 378 OVC patients from TCGA database. The RNA-seq and clinical data of 379 OVC patients in GSE140082 and 173 OV patients in GSE53963. The PROGENy score was calculated to assess tumor-associated pathways. Based on gene set enrichment analysis (GSEA) of the cuproptosis pathway, the single cells were divided into the cuproptosishigh and cuproptosislow groups. The differentially expressed genes (DEGs) between the two groups were screened, and 47 prognosis-related genes were identified based on univariate cox regression analysis. Randomforest was used to construct a prognostic model. Immuno-infiltration analysis was performed using ssGSEA and xCell algorithms. In vitro and in vivo experiments were used for functional verification.

RESULTS

Six major cell populations was identified, including fibroblast, T cell, myeloid, epithelial cell, endothelial cell, and B cell populations. The PROGENy score which revealed significant activation of the PI3K pathway in T and B cells, and activation of the TGF-β pathway in endothelial cells and fibroblasts. TIMM8B, COX8A, SSR4, HIGD2A, WASF2, PRDX5 and CLDN4 were selected to construct a prognostic model from the identified 47 prognosis-related genes. Furthermore, the cuproptosishigh and cuproptosislow groups showed significant differences in the expression levels of the model genes, immune cell infiltration, and sensitivity to six potential drug candidates. The functional experiments showed that WASF2 is associated with cuproptotic resistance and promotes cancer cell proliferation and resistance to platinum, and its high expression is associated with poor prognosis of OVC patients.

CONCLUSION

A clinically significant cuproptosis-related prognostic model was identified which can accurately predict the prognosis and immune characteristics of OVC patients. WASF2, one of the cuproptosis-related gene in the risk model, promotes the proliferation and platinum resistance of OVC cells, and leads poor prognosis.

Single-cell and bulk RNA-sequencing analysis to predict the role and clinical value of CD36 in lung squamous cell carcinoma

The majority of patients with lung squamous cell carcinoma are diagnosed at an advanced stage, which poses a challenge to the efficacy of chemotherapy. Therefore, the search for an early biomarker needs to be addressed. CD36 is a scavenger receptor expressed in various cell types. It has been reported that it is closely related to the occurrence and development of many kinds of tumours. However, its role in lung squamous cell carcinoma has not been reported. Our research aims to reveal the role of CD36 in lung squamous cell carcinoma by integrating single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing data. We used bioinformatics methods to explore the potential carcinogenicity of CD36 by analysing the data from the cancer genome map (TCGA), gene expression comprehensive map (GEO), human protein map (HPA) comparative toxicology genomics database (CTD) and other resources. Our study dissected the relationship between CD36 and prognosis and gene correlation, functional analysis, mutation of different tumours, infiltration of immune cells and exploring the interaction between CD36 and chemicals. The results showed that the expression of CD36 was heterogeneous. Compared with normal patients, the expression was low in lung squamous cell carcinoma. In addition, CD36 showed early diagnostic value in four kinds of tumours (LUSC, BLCA, BRCA and KIRC) and was positively or negatively correlated with the prognosis of different tumours. The relationship between CD36 and the tumour immune microenvironment was revealed by immunoinfiltration analysis, and many drugs that might target CD36 were identified by the comparative toxicological genomics database (CTD). In summary, through pancancer analysis, we found and verified for the first time that CD36 may play a role in the detection of lung squamous cell carcinoma. In addition, it has high specificity and sensitivity in detecting cancer. Therefore, CD36 can be used as an auxiliary index for early tumour diagnosis and a prognostic marker for lung squamous cell carcinoma.

Integrated single-cell RNA sequencing analysis reveals a mesenchymal stem cell-associated signature for estimating prognosis and drug sensitivity in gastric cancer

BACKGROUND

Mesenchymal stem cells (MSCs) play an important role in regulating all stages of the immune response, angiogenesis, and transformation of matrix components in the tumor microenvironment. The aim of this study was to identify the prognostic value of MSC-related signatures in patients with gastric cancer (GC).

METHODS

MSC marker genes were identified by analyzing single-cell RNA sequencing (scRNA-seq) data for GC from the Gene Expression Omnibus (GEO) database. Using bulk sequencing data from the Cancer Genome Atlas-Stomach adenocarcinoma (TCGA-STAD), as a training cohort, and data from GEO, as a validation cohort, we developed a risk model consisting of MSC prognostic signature genes, and classified GC patients into high- and low-MSC risk subgroups. Multifactorial Cox regression was used to evaluate whether MSC prognostic signature was an independent prognostic factor. An MSC nomogram was constructed combining clinical information and risk grouping. Subsequently, we evaluated the effect of MSC prognostic signature on immune cell infiltration, antitumor drugs and immune checkpoints and verified the expression of MSC prognostic signature by in vitro cellular assays.

RESULTS

In this study, 174 MSC marker genes were identified by analyzing scRNA-seq data. We identified seven genes (POSTN, PLOD2, ITGAV, MMP11, SDC2, MARCKS, ANXA5) to construct MSC prognostic signature. MSC prognostic signature was an independent risk factor in the TCGA and GEO cohorts. GC patients in the high-MSC risk group had worse prognoses. In addition, the MSC nomogram has a high clinical application value. Notably, the MSC signature can induce the development of a poor immune microenvironment. GC patients in the high MSC-risk group were more sensitive to anticancer drugs and tended to have higher levels of immune checkpoint markers. In qRT-PCR assays, the MSC signature was more highly expressed in GC cell lines.

CONCLUSIONS

The MSC marker gene-based risk signature developed in this study can not only be used to predict the prognosis of GC patients, but also has the potential to reflect the efficacy of antitumor therapies.

Clathrin light chain A-enriched small extracellular vesicles remodel microvascular niche to induce hepatocellular carcinoma metastasis

Small extracellular vesicles (sEVs) play a key role in exchanging cargoes between cells in tumour microenvironment. This study aimed to elucidate the functions and mechanisms of hepatocellular carcinoma (HCC) derived sEV-clathrin light chain A (CLTA) in remodelling microvascular niche. CLTA level in the circulating sEVs of HCC patients was analysed by enzyme-linked immunosorbent assay (ELISA). The functions of sEV-CLTA in affecting HCC cancerous properties were examined by multiple functional assays. Mass spectrometry was used to identify downstream effectors of sEV-CLTA in human umbilical vein endothelial cells (HUVECs). Tube formation, sprouting, trans-endothelial invasion and vascular leakiness assays were performed to determine the functions of sEV-CLTA and its effector, basigin (BSG) in HUVECs. BSG inhibitor, SP-8356, was tested in a mouse model of patient-derived xenografts (PDXs). Circulating sEVs of HCC patients had markedly enhanced CLTA levels than control individuals and were reduced in patients after surgery. HCC derived sEV-CLTA enhanced HCC cancerous properties, disrupted endothelial integrity and induced angiogenesis. Mechanistically, CLTA remodels microvascular niche by stabilizing and upregulating BSG. Last, SP-8356 alone or in combination with sorafenib attenuated PDXs growth. The study reveals the role of HCC derived sEV-CLTA in microvascular niche formation. Inhibition of CLTA and its mediated pathway may illuminate a new therapeutic strategy for HCC patients.

The prognostic and biological importance of chromatin regulation-related genes for lung cancer is examined using bioinformatics and experimentally confirmed

BACKGROUND

The pathogenesis and clinical diagnosis of lung adenocarcinoma (LUAD), a malignant illness with substantial morbidity and mortality, are still being investigated. Genes involved in chromatin regulation are crucial in the biological function of LUAD.

METHODS

The prognostic prediction model for LUAD was developed using multivariables and least absolute shrinkage and selection operator (LASSO) regression. It consisted of 10 chromatin regulators. The LUAD has been divided into two groups, high- and low-risk, using a predictive model. The model was shown to be accurate in predicting survival by the nomogram, receiver operating characteristic (ROC) curves, and principal component analysis (PCA). An analysis of differences in immune-cell infiltration, immunologicalfunction, and clinical traits between low- and high-risk populations was conducted. Protein-protein interaction (PPI) networks and Gene Ontology (GO) pathways of differentially expressed genes (DEGs) in the high versus low risk group were also examined to investigate the association between genes and biological pathways. The biological roles of chromatin regulators (CRs) in LUAD were finally estimated using colony formation and cell movement. The important genes' mRNA expression has been measured using real-time polymerase chain reaction (RT-PCR).

RESULTS AND CONCLUSION

Risk score and stage based on the model could be seen as separate prognostic indicators for patients with LUAD. The main signaling pathway difference across various risk groups was in cell cycle. The immunoinfiltration profile of the tumor microenvironment (TME) and individuals with different risk levels were correlated, suggesting that the interaction of immune cells with the tumor led to the creation of a favorable immunosuppressive microenvironment. These discoveries aid in the creation of individualized therapies for LUAD patients.

Identification of fibroblast-related genes based on single-cell and machine learning to predict the prognosis and endocrine metabolism of pancreatic cancer

Background

Single-cell sequencing technology has become an indispensable tool in tumor mechanism and heterogeneity studies. Pancreatic adenocarcinoma (PAAD) lacks early specific symptoms, and comprehensive bioinformatics analysis for PAAD contributes to the developmental mechanisms.

Methods

We performed dimensionality reduction analysis on the single-cell sequencing data GSE165399 of PAAD to obtain the specific cell clusters. We then obtained cell cluster-associated gene modules by weighted co-expression network analysis and identified tumorigenesis-associated cell clusters and gene modules in PAAD by trajectory analysis. Tumor-associated genes of PAAD were intersected with cell cluster marker genes and within the signature module to obtain genes associated with PAAD occurrence to construct a prognostic risk assessment tool by the COX model. The performance of the model was assessed by the Kaplan-Meier (K-M) curve and the receiver operating characteristic (ROC) curve. The score of endocrine pathways was assessed by ssGSEA analysis.

Results

The PAAD single-cell dataset GSE165399 was filtered and downscaled, and finally, 17 cell subgroups were filtered and 17 cell clusters were labeled. WGCNA analysis revealed that the brown module was most associated with tumorigenesis. Among them, the brown module was significantly associated with C11 and C14 cell clusters. C11 and C14 cell clusters belonged to fibroblast and circulating fetal cells, respectively, and trajectory analysis showed low heterogeneity for fibroblast and extremely high heterogeneity for circulating fetal cells. Next, through differential analysis, we found that genes within the C11 cluster were highly associated with tumorigenesis. Finally, we constructed the RiskScore system, and K-M curves and ROC curves revealed that RiskScore possessed objective clinical prognostic potential and demonstrated consistent robustness in multiple datasets. The low-risk group presented a higher endocrine metabolism and lower immune infiltrate state.

Conclusion

We identified prognostic models consisting of APOL1, BHLHE40, CLMP, GNG12, LOX, LY6E, MYL12B, RND3, SOX4, and RiskScore showed promising clinical value. RiskScore possibly carries a credible clinical prognostic potential for PAAD.

Identification of a novel signature based on macrophage-related marker genes to predict prognosis and immunotherapeutic effects in hepatocellular carcinoma

Background

Tumor-related macrophages (TAMs) have emerged as an essential part of the immune regulatory network in hepatocellular carcinoma (HCC). Constructing a TAM-related signature is significant for evaluating prognosis and immunotherapeutic response of HCC patients.

Methods

Informative single-cell RNA sequencing (scRNA-seq) dataset was obtained from the Gene Expression Omnibus (GEO) database, and diverse cell subpopulations were identified by clustering dimension reduction. Moreover, we determined molecular subtypes with the best clustering efficacy by calculating the cumulative distribution function (CDF). The ESTIMATE method, CIBERSORT (cell-type identification by estimating relative subsets of RNA transcripts) algorithm and publicly available tumor immune dysfunction and exclusion (TIDE) tools were used to characterize the immune landscape and tumor immune escape status. A TAM-related gene risk model was constructed through Cox regression and verified in multiple datasets and dimensions. We also performed functional enrichment analysis to detect potential signaling pathways related to TAM marker genes.

Results

In total, 10 subpopulations and 165 TAM-related marker genes were obtained from the scRNA-seq dataset (GSE149614). After clustering 3 molecular subtypes based on TAM-related marker genes, we found significantly different prognostic survival and immune signatures among the three subtypes. Subsequently, a 9-gene predictive signature (TPP1, FTL, CXCL8, CD68, ATP6V1F, CSTB, YBX1, LGALS3, and APLP2) was identified as an independent prognostic factor for HCC patients. Those patients with high RiskScore had a lower survival rate and benefited less from immunotherapy than those with low RiskScore. Moreover, more samples of the Cluster C subtype were enriched in the high-risk group, with higher tumor immune escape incidence.

Conclusions

We constructed a TAM-related signature with excellent efficacy for predicting prognostic survival and immunotherapeutic responses in HCC patients.

An immune-related prognostic model predicts neoplasm-immunity interactions for metastatic nasopharyngeal carcinoma

Background

The prognosis of nasopharyngeal carcinoma (NPC) has been recognized to improve immensely owing to radiotherapy combined with chemotherapy. However, patients with metastatic NPC have a poor prognosis. Immunotherapy has dramatically prolonged the survival of patients with NPC. Hence, further research on immune-related biomarkers is imperative to establish the prognosis of metastatic NPC.

Methods

10 NPC RNA expression profiles were generated from patients with or without distant metastasis after chemoradiotherapy from the Fujian Cancer Hospital. The differential immune-related genes were identified and validated by immunohistochemistry analysis. The method of least absolute shrinkage and selection operator (LASSO)was used to further establish the immune-related prognostic model in an external GEO database (GSE102349, n=88). The immune microenvironment and signal pathways were evaluated in multiple dimensions at the transcriptome and single-cell levels.

Results

1328 differential genes were identified, out of which 520 were upregulated and 808 were downregulated. Notably, most of the immune genes and pathways were down-regulated in the metastasis group. A prognostic immune model involving nine hub genes. Patients in low-risk group were characterized by survival advantage, hot immune phenotype and benefit from immunotherapy. Compared with immune cells, malignant cell exhibited the most active levels of risk score by ssGSEA. Accordingly, intercellular communications including LT, CD70, CD40 and SPP1, and the like, between high-risk and low-risk were explored by the R package “Cellchat”.

Conclusion

We have constructed a model based on immunity of metastatic NPC and determined its prognostic value. The model identified the level of immune cell infiltration, cell-cell communication, along with potential immunotherapy for metastatic NPC.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

Significance of Identifying Key Genes Involved in HBV-Related Hepatocellular Carcinoma for Primary Care Surveillance of Patients with Cirrhosis

Cirrhosis is frequently the final stage of disease preceding the development of hepatocellular carcinoma (HCC) and is one of the risk factors for HCC. Preventive surveillance for early HCC in patients with cirrhosis is advantageous for achieving early HCC prevention and diagnosis, thereby enhancing patient prognosis and reducing mortality. However, there is no highly sensitive diagnostic marker for the clinical surveillance of HCC in patients with cirrhosis, which significantly restricts its use in primary care for HCC. To increase the accuracy of illness diagnosis, the study of the effective and sensitive genetic biomarkers involved in HCC incidence is crucial. In this study, a set of 120 significantly differentially expressed genes (DEGs) was identified in the GSE121248 dataset. A protein-protein interaction (PPI) network was constructed among the DEGs, and Cytoscape was used to extract hub genes from the network. In TCGA database, the expression levels, correlation analysis, and predictive performance of hub genes were validated. In total, 15 hub genes showed increased expression, and their positive correlation ranged from 0.80 to 0.90, suggesting they may be involved in the same signaling pathway governing HBV-related HCC. The GSE10143, GSE25097, GSE54236, and GSE17548 datasets were used to investigate the expression pattern of these hub genes in the progression from cirrhosis to HCC. Using Cox regression analysis, a prediction model was then developed. The ROC curves, DCA, and calibration analysis demonstrated the superior disease prediction accuracy of this model. In addition, using proteomic analysis, we investigated whether these key hub genes interact with the HBV-encoded oncogene X protein (HBx), the oncogenic protein in HCC. We constructed stable HBx-expressing LO2-HBx and Huh-7-HBx cell lines. Co-immunoprecipitation coupled with mass spectrometry (Co-IP/MS) results demonstrated that CDK1, RRM2, ANLN, and HMMR interacted specifically with HBx in both cell models. Importantly, we investigated 15 potential key genes (CCNB1, CDK1, BUB1B, ECT2, RACGAP1, ANLN, PBK, TOP2A, ASPM, RRM2, NEK2, PRC1, SPP1, HMMR, and DTL) participating in the transformation process of HBV infection to HCC, of which 4 hub genes (CDK1, RRM2, ANLN, and HMMR) probably serve as potential oncogenic HBx downstream target molecules. All these findings of our study provided valuable research direction for the diagnostic gene detection of HBV-related HCC in primary care surveillance for HCC in patients with cirrhosis.

Integrating Single-cell RNA-seq to construct a Neutrophil prognostic model for predicting immune responses in non-small cell lung cancer

Non-small cell lung cancer (NSCLC) is the most widely distributed tumor in the world, and its immunotherapy is not practical. Neutrophil is one of a tumor’s most abundant immune cell groups. This research aimed to investigate the complex communication network in the immune microenvironment (TIME) of NSCLC tumors to clarify the interaction between immune cells and tumors and establish a prognostic risk model that can predict immune response and prognosis of patients by analyzing the characteristics of Neutrophil differentiation. Integrated Single-cell RNA sequencing (scRNA-seq) data from NSCLC samples and Bulk RNA-seq were used for analysis. Twenty-eight main cell clusters were identified, and their interactions were clarified. Next, four subsets of Neutrophils with different differentiation states were found, closely related to immune regulation and metabolic pathways. Based on the ratio of four housekeeping genes (ACTB, GAPDH, TFRC, TUBB), six Neutrophil differentiation-related genes (NDRGs) prognostic risk models, including MS4A7, CXCR2, CSRNP1, RETN, CD177, and LUCAT1, were constructed by Elastic Net and Multivariate Cox regression, and patients’ total survival time and immunotherapy response were successfully predicted and validated in three large cohorts. Finally, the causes of the unfavorable prognosis of NSCLC caused by six prognostic genes were explored, and the small molecular compounds targeted at the anti-tumor effect of prognostic genes were screened. This study clarifies the TIME regulation network in NSCLC and emphasizes the critical role of NDRGs in predicting the prognosis of patients with NSCLC and their potential response to immunotherapy, thus providing a promising therapeutic target for NSCLC.

Analysis on heterogeneity of hepatocellular carcinoma immune cells and a molecular risk model by integration of scRNA-seq and bulk RNA-seq

Background

Studies have shown that hepatocellular carcinoma (HCC) heterogeneity is a main cause leading to failure of treatment. Technology of single-cell sequencing (scRNA) could more accurately reveal the essential characteristics of tumor genetics.

Methods

From the Gene Expression Omnibus (GEO) database, HCC scRNA-seq data were extracted. The FindCluster function was applied to analyze cell clusters. Autophagy-related genes were acquired from the MSigDB database. The ConsensusClusterPlus package was used to identify molecular subtypes. A prognostic risk model was built with the Least Absolute Shrinkage and Selection Operator (LASSO)-Cox algorithm. A nomogram including a prognostic risk model and multiple clinicopathological factors was constructed.

Results

Eleven cell clusters labeled as various cell types by immune cell markers were obtained from the combined scRNA-seq GSE149614 dataset. ssGSEA revealed that autophagy-related pathways were more enriched in malignant tumors. Two autophagy-related clusters (C1 and C2) were identified, in which C1 predicted a better survival, enhanced immune infiltration, and a higher immunotherapy response. LASSO-Cox regression established an eight-gene signature. Next, the HCCDB18, GSA14520, and GSE76427 datasets confirmed a strong risk prediction ability of the signature. Moreover, the low-risk group had enhanced immune infiltration and higher immunotherapy response. A nomogram which consisted of RiskScore and clinical features had better prediction ability.

Conclusion

To precisely assess the prognostic risk, an eight-gene prognostic stratification signature was developed based on the heterogeneity of HCC immune cells.

Identification of Diagnostic Genes and Effective Drugs Associated with Osteoporosis Treatment by Single-Cell RNA-Seq Analysis and Network Pharmacology

Background

Osteoporosis is a common bone metabolic disease with increased bone fragility and fracture rate. Effective diagnosis and treatment of osteoporosis still need to be explored due to the increasing incidence of disease.

Methods

Single-cell RNA-seq was acquired from GSE147287 dataset. Osteoporosis-related genes were obtained from chEMBL. Cell subpopulations were identified and characterized by scRNA-seq, t-SNE, clusterProfiler, and other computational methods. “limma” R packages were used to identify all differentially expressed genes. A diagnosis model was build using rms R packages. Key drugs were determined by proteins-proteins interaction and molecular docking.

Results

Firstly, 15,577 cells were obtained, and 12 cell subpopulations were identified by clustering, among which 6 cell subpopulations belong to CD45+ BM-MSCs and the other subpopulations were CD45-BM-MSCs. CD45- BM-MSCs_6 and CD45+ BM-MSCs_5 were consider as key subpopulations. Furthermore, we found 7 genes were correlated with above two subpopulations, and F9 gene had highest AUC. Finally, five compounds were identified, among which DB03742 bound well to F9 protein.

Conclusions

This work discovered that 7 genes were correlated with CD45-BM-MSCs_6 and CD45+ BM-MSCs_5 subpopulations in osteoporosis, among which F9 gene had better research value. Moreover, compound DB03742 was a potential inhibitor of F9 protein.

Single-Cell Transcriptomics of Liver Cancer: Hype or Insights?

Hepatocellular carcinoma (HCC) is characterized by its high degrees of both inter- and intratumoral heterogeneity. Its complex tumor microenvironment is also crucial in promoting tumor progression. Recent advances in single-cell RNA sequencing provide an important highway to characterize the underlying pathogenesis and heterogeneity of HCC in an unprecedented degree of resolution. This review discusses the up-to-date discoveries from the latest studies of HCC with respect to the strength of single-cell RNA sequencing. We discuss its use in the dissection of the landscape of the intricate HCC ecosystem and highlight the major features at cellular levels, including the malignant cells, different immune cell types, and the various cell-cell interactions, which are crucial for developing effective immunotherapies. Finally, its translational applications will be discussed. Altogether, these explorations may give us some hints at the tumor growth and progression and drug resistance and recurrence, particularly in this era of personalized medicine.