Integrating bulk and single-cell RNA sequencing reveals cellular heterogeneity and immune infiltration in hepatocellular carcinoma

Single-cell analysis extracted CAFs-related genes to established online app to predict clinical outcome and radiotherapy prognosis of prostate cancer

BACKGROUND

Cancer-associated fibroblasts (CAFs) play a significant role in regulating the clinical outcome and radiotherapy prognosis of prostate cancer (PCa). The aim of this study is to identify CAFs-related genes (CAFsRGs) using single-cell analysis and evaluate their potential for predicting the prognosis and radiotherapy prognosis in PCa.

METHODS

We acquire transcriptome and single-cell RNA sequencing (scRNA-seq) results of PCa and normal adjacent tissues from The GEO and TCGA databases. The "MCPcounter" and "EPIC" R packages were used to assess the infiltration level of CAFs and examine their correlation with PCa prognosis. ScRNA-seq and differential gene expression analyses were used to extract CAFsRGs. We also applied COX and LASSO analysis to further construct a risk score (CAFsRS) to assess biochemical recurrence-free survival (BRFS) and radiotherapy prognosis of PCa. The predictive efficacy of CAFsRS was evaluated by ROC curves and subgroup analysis. Finally, we integrated the CAFsRS gene signature with relevant clinical features to develop a nomogram, enhancing the predictive accuracy.

RESULTS

The abundance of CAFs is associated with a poor prognosis of PCa patients. ScRNA-seq and differential gene expression analysis revealed 323 CAFsRGs. After COX and LASSO analysis, we obtained seven CAFsRGs with prognostic significance (PTGS2, FKBP10, ENG, CDH11, COL5A1, COL5A2, and SRD5A2). Additionally, we established a risk score model based on the training set (n = 257). The ROC curve was used to confirm the performance of CAFsRS (The AUC values for 1, 3 and 5-year survival were determined to be 0.732, 0.773, and 0.775, respectively.). The testing set (n = 129), GSE70770 set (n = 199) and GSE116918 set (n = 248) revealed that the model exhibited exceptional predictive performance. This was also confirmed by clinical subgroup analysis. The violin plot demonstrated a statistically significant disparity in the CAFs infiltrations between the high-risk and low-risk groups of CAFsRS. Further analysis confirmed that both CAFsRS and T stage were independent prognostic factors for PCa. The nomogram was then established and its excellent predictive performance was demonstrated through calibration and ROC curves. Finally, we developed an online prognostic prediction app ( https://sysu-symh-cafsnomogram.streamlit.app/ ) to facilitate the practical application of the nomogram.

CONCLUSIONS

The prognostic prediction risk score model we constructed could accurately predict BRFS and radiotherapy prognosis PCa, which can provide new ideas for clinicians to develop personalized PCa treatment and follow-up programs.

Antitumour mechanisms of traditional Chinese medicine elicited by regulating tumour-associated macrophages in solid tumour microenvironments

Tumour-associated macrophages (TAMs), particularly M2-TAMs, constitute the largest proportion of immune cells in the solid tumour microenvironment, playing a crucial role in tumour progression and correlating with poor prognosis. TAMs promote the proliferation, invasion, and metastasis of tumour cells by remodelling the extracellular matrix, inhibiting immunity, promoting immune escape and tumour angiogenesis, and affecting cell metabolism. Traditional Chinese medicine (TCM) has been used clinically in China for millennia. Chinese herbs exhibit potent antitumour effects with minimal to no toxicity, substantially contributing to prolonging the lives of patients with cancer and improving their quality of life. TCM has unique advantages in improving the solid tumour microenvironment, particularly in regulating TAMs to further inhibit tumour angiogenesis, reduce drug resistance, reverse immunosuppression, and enhance antitumour immunity. This review highlights the TAM-associated mechanisms within the solid tumour microenvironment, outlines the recent advancements in TCM targeting TAMs for antitumour effects, emphasises the superiority of combining TCM with standard treatments or new nano-drug delivery systems, and evaluates the safety and efficacy of TCM combined with conventional treatments via clinical trials to provide insights and strategies for future research and clinical treatment.

Signature construction and molecular subtype identification based on liver-specific genes for prediction of prognosis, immune activity, and anti-cancer drug sensitivity in hepatocellular carcinoma

BACKGROUND

Liver specific genes (LSGs) are crucial for hepatocyte differentiation and maintaining normal liver function. A deep understanding of LSGs and their heterogeneity in hepatocellular carcinoma (HCC) is necessary to provide clues for HCC diagnosis, prognosis, and treatment.

METHODS

The bulk and single-cell RNA-seq data of HCC were downloaded from TCGA, ICGC, and GEO databases. Through unsupervised cluster analysis, LSGs-based HCC subtypes were identified in TCGA-HCC samples. The prognostic effects of the subtypes were investigated with survival analyses. With GSVA and Wilcoxon test, the LSGs score, stemness score, aging score, immune score and stromal score of the samples were estimated and compared. The HCC subtype-specific genes were identified. The subtypes and their differences were validated in ICGC-HCC samples. LASSO regression analysis was used for key gene selection and risk model construction for HCC overall survival. The model performance was estimated and validated. The key genes were validated for their heterogeneities in HCC cell lines with quantitative real-time PCR and at single-cell level. Their dysregulations were investigated at protein level. Their correlations with HCC response to anti-cancer drugs were estimated in HCC cell lines.

RESULTS

We identified three LSGs-based HCC subtypes with different prognosis, tumor stemness, and aging level. The C1 subtype with low LSGs score and high immune score presented a poor survival, while the C2 subtype with high LSGs score and immune score indicated an enduring survival. Although no significant survival difference between C2 and C3 HCCs was shown, the C2 HCCs presented higher immune score and stroma score. The HCC subtypes and their differences were confirmed in ICGC-HCC dataset. A five-gene prognostic signature for HCC survival was constructed. Its good performance was shown in both the training and validation datasets. The five genes presented significant heterogeneities in different HCC cell lines and hepatocyte subclusters. Their dysregulations were confirmed at protein level. Furthermore, their significant associations with HCC sensitivities to anti-cancer drugs were shown.

CONCLUSIONS

LSGs-based HCC subtype classification and the five-gene risk model might provide useful clues not only for HCC stratification and risk prediction, but also for the development of more personalized therapies for effective HCC treatment.

Integrated Analysis of Single-Cell and Bulk RNA Sequencing Data Reveals Memory-like NK Cell Subset Associated with Mycobacterium tuberculosis Latency

Natural killer (NK) cells are innate-like lymphocytes that belong to the family of type-1 innate lymphoid cells and rapidly respond to virus-infected and tumor cells. In this study, we have combined scRNA-seq data and bulk RNA-seq data to define the phenotypic and molecular characteristics of peripheral blood NK cells. While the role of NK cells in immune surveillance against virus infections and tumors has been well established, their contribution to protective responses to other intracellular microorganisms, such as Mycobacterium tuberculosis (Mtb), is still poorly understood. In this study, we have combined scRNA-seq data and bulk RNA-seq data to illuminate the molecular characteristics of circulating NK cells in patients with active tuberculosis (TB) disease and subjects with latent Mtb infection (LTBI) and compared these characteristics with those of healthy donors (HDs) and patients with non-TB other pulmonary infectious diseases (ODs). We show here that the NK cell cluster was significantly increased in LTBI subjects, as compared to patients with active TB or other non-TB pulmonary diseases and HD, and this was mostly attributable to the expansion of an NK cell population expressing KLRC2, CD52, CCL5 and HLA-DRB1, which most likely corresponds to memory-like NK2.1 cells. These data were validated by flow cytometry analysis in a small cohort of samples, showing that LTBI subjects have a significant expansion of NK cells characterized by the prevalence of memory-like CD52+ NKG2C+ NK cells. Altogether, our results provide some new information on the role of NK cells in protective immune responses to Mtb.

T-cell infiltration and its regulatory mechanisms in cancers: insights at single-cell resolution

Tumor-infiltrating T cells recognize, attack, and clear tumor cells, playing a central role in antitumor immune response. However, certain immune cells can impair this response and help tumor immune escape. Therefore, exploring the factors that influence T-cell infiltration is crucial to understand tumor immunity and improve therapeutic effect of cancer immunotherapy. The use of single-cell RNA sequencing (scRNA-seq) allows the high-resolution analysis of the precise composition of immune cells with different phenotypes and other microenvironmental factors, including non-immune stromal cells and the related molecules in the tumor microenvironment of various cancer types. In this review, we summarized the research progress on T-cell infiltration and the crosstalk of other stromal cells and cytokines during T-cell infiltration using scRNA-seq to provide insights into the mechanisms regulating T-cell infiltration and contribute new perspectives on tumor immunotherapy.

The counterpart congenital overgrowth syndromes Beckwith-Wiedemann Syndrome in human and large offspring syndrome in bovine involve alterations in DNA methylation, transcription, and chromatin configuration

Beckwith-Wiedemann Syndrome (BWS, OMIM #130650) is a congenital epigenetic disorder in humans which affects approximately 1 in 10,340 children. The incidence is likely an underestimation as the condition is usually recognized based on observable phenotypes at birth. BWS children have up to a 28% risk of developing tumors and currently, only 80% of patients can be corroborated molecularly (epimutations/variants). It is unknown how the subtypes of this condition are molecularly similar/dissimilar globally, therefore there is a need to deeply characterize the syndrome at the molecular level. Here we characterize the methylome, transcriptome and chromatin configuration of 18 BWS individuals together with the animal model of the condition, the bovine large offspring syndrome (LOS). Sex specific comparisons are performed for a subset of the BWS patients and LOS. Given that this epigenetic overgrowth syndrome has been characterized as a loss-of-imprinting condition, parental allele-specific comparisons were performed using the bovine animal model. In general, the differentially methylated regions (DMRs) detected in BWS and LOS showed significant enrichment for CTCF binding sites. Altered chromosome compartments in BWS and LOS were positively correlated with gene expression changes, and the promoters of differentially expressed genes showed significant enrichment for DMRs, differential topologically associating domains, and differential A/B compartments in some comparisons of BWS subtypes and LOS. We show shared regions of dysregulation between BWS and LOS, including several HOX gene clusters, and also demonstrate that altered DNA methylation differs between the clinically epigenetically identified BWS patients and those identified as having DNA variants (i.e. CDKN1C microdeletion). Lastly, we highlight additional genes and genomic regions that have the potential to serve as targets for biomarker development to improve current molecular methodologies. In summary, our results suggest that genome-wide alternation of chromosome architecture, which is partially caused by DNA methylation changes, also contribute to the development of BWS and LOS.

Application of adoptive cell therapy in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) remains a global health challenge. Novel treatment modalities are urgently needed to extend the overall survival of patients. The liver plays an immunomodulatory function due to its unique physiological structural characteristics. Therefore, following surgical resection and radiotherapy, immunotherapy regimens have shown great potential in the treatment of hepatocellular carcinoma. Adoptive cell immunotherapy is rapidly developing in the treatment of hepatocellular carcinoma. In this review, we summarize the latest research on adoptive immunotherapy for hepatocellular carcinoma. The focus is on chimeric antigen receptor (CAR)-T cells and T cell receptor (TCR) engineered T cells. Then tumour-infiltrating lymphocytes (TILs), natural killer (NK) cells, cytokine-induced killer (CIK) cells, and macrophages are briefly discussed. The main overview of the application and challenges of adoptive immunotherapy in hepatocellular carcinoma. It aims to provide the reader with a comprehensive understanding of the current status of HCC adoptive immunotherapy and offers some strategies. We hope to provide new ideas for the clinical treatment of hepatocellular carcinoma.

CDK1-SRC Interaction-Dependent Transcriptional Activation of HSP90AB1 Promotes Antitumor Immunity in Hepatocellular Carcinoma

This study aimed to analyze multiomics data and construct a regulatory network involving kinases, transcription factors, and immune genes in hepatocellular carcinoma (HCC) prognosis. The researchers used transcriptomic, proteomic, and clinical data from TCGA and GEO databases to identify immune genes associated with HCC. Statistical analysis, meta-analysis, and protein-protein interaction analyses were performed to identify key immune genes and their relationships. In vitro and in vivo experiments validated the CDK1-SRC-HSP90AB1 network's effects on HCC progression and antitumor immunity. A prognostic risk model was developed using clinicopathological features and immune infiltration. The immune genes LPA, BIRC5, HSP90AB1, ROBO1, and CCL20 were identified as the key prognostic factors. The CDK1-SRC-HSP90AB1 network promoted HCC cell proliferation and migration, with HSP90AB1 being transcriptionally activated by the CDK1-SRC interaction. Manipulating SRC or HSP90AB1 reversed the effects of CDK1 and SRC on HCC. The CDK1-SRC-HSP90AB1 network also influenced HCC tumor formation and antitumor immunity. Overall, this study highlights the importance of the CDK1-SRC-HSP90AB1 network as a crucial immune-regulatory network in the HCC prognosis.

Immune landscape and immunotherapy of hepatocellular carcinoma: focus on innate and adaptive immune cells

Hepatocellular carcinoma (HCC) is responsible for roughly 90% of all cases of primary liver cancer, and the cases are on the rise. The treatment of advanced HCC is a serious challenge. Immune checkpoint inhibitor (ICI) therapy has marked a watershed moment in the history of HCC systemic treatment. Atezolizumab in combination with bevacizumab has been approved as a first-line treatment for advanced HCC since 2020; however, the combination therapy is only effective in a limited percentage of patients. Considering that the tumor immune microenvironment (TIME) has a great impact on immunotherapies for HCC, an in-depth understanding of the immune landscape in tumors and the current immunotherapeutic approaches is extremely necessary. We elaborate on the features, functions, and cross talk of the innate and adaptive immune cells in HCC and highlight the benefits and drawbacks of various immunotherapies for advanced HCC, as well as future projections. HCC consists of a heterogeneous group of cancers with distinct etiologies and immune microenvironments. Almost all the components of innate and adaptive immune cells in HCC have altered, showing a decreasing trend in the number of tumor suppressor cells and an increasing trend in the pro-cancer cells, and there is also cross talk between various cell types. Various immunotherapies for HCC have also shown promising efficacy and application prospect. There are multilayered interwoven webs among various immune cell types in HCC, and emerging evidence demonstrates the promising prospect of immunotherapeutic approaches for HCC.

A prognostic model related to necrotizing apoptosis of breast cancer based on biorthogonal constrained depth semi-supervised nonnegative matrix decomposition and single-cell sequencing analysis

Breast cancer (BC) is one of the most common malignant tumours in women, and its prognosis is poor. The prognosis of BC patients can be improved by immunotherapy. However, due to the heterogeneity of BC, the identification of new biomarkers is urgently needed to improve the prognosis of BC patients. Necrotic apoptosis has been shown to play an essential role in many cancers. First, this study proposed a novel clustering algorithm called biorthogonal constrained depth semisupervised nonnegative matrix factorization (DO-DSNMF). The DO-DSNMF algorithm added multilayer nonlinear transformation to the coefficient matrix obtained after decomposition, which was used to mine the nonlinear relationship between samples. In addition, we also added orthogonal constraints on the basis matrix and coefficient matrix to reduce the influence of redundant features and samples on the results. We applied the DO-DSNMF algorithm and analysed the differences in survival and immunity between the subtypes. Then, we used prognosis analysis to construct the prognosis model. Finally, we analysed single cells using single-cell sequencing (scRNA-seq) data from the GSE75688 dataset in the GEO database. We identified two BC subtypes based on the BC transcriptome data in the TCGA database. Immune infiltration analysis showed that the necrotizing apoptosis-related genes of BC were related to various immune cells and immune functions. Necrotizing apoptosis was found to play a role in BC progression and immunity. The role of prognosis-related NRGs in BC was also verified by cell experiments. This study proposed a novel clustering algorithm to analyse BC subtypes and constructed an NRG prognostic model for BC. The prognosis and immune landscape of BC patients were evaluated by this model. The cell experiment supported its role in BC, which provides a potential therapeutic target for the treatment of BC.

Iridium(III) complexes inhibit the proliferation and migration of BEL-7402 cells through the PI3K/AKT/mTOR signaling pathway

Iridium(III) complexes are largely studied as anti-cancer complexes due to their excellent anti-cancer activity. In this article, two new iridium(III) complexes [Ir(piq)₂(THPIP)]PF₆ (THPIP = 2,4-di-tert-butyl-6-(1H-imidazo[4,5-f][1,10]phenanthrolin-2-yl)phenol, piq = deprotonated 1-phenylisoquinoline) (Ir1) and [Ir(bzq)₂(THPIP)]PF₆ (bzq = deprotonated benzo[h]quinolone) (Ir2) were synthesized. 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assays showed that complex Ir1 exhibits moderate activity (IC₅₀ = 29.9 ± 4.6 μM) and Ir2 shows high cytotoxicity (IC₅₀ = 9.8 ± 1.8 μM) against BEL-7402 cells. Further studies on the mechanism showed that Ir1 and Ir2 induced apoptosis by changing the mitochondrial membrane potential, Ca²⁺ release, ROS accumulation, and cell cycle arrest at the S phase. The complexes can effectively inhibit cell colony formation and migration. The expression of B-cell lymphoma-2 (Bcl-2) family proteins, PI3K (phosphatidylinositol 3-kinase), AKT (protein kinase B), mTOR (mammalian target of rapamycin), and p-mTOR was studied by immunoblotting. Complexes Ir1 and Ir2 downregulated the expression of anti-apoptotic protein Bcl-2 and increased the expression of autophagy-related proteins of Beclin-1 and LC3-II. Further experiments showed that the complexes inhibited the production of glutathione (GSH) and increased the amounts of malondialdehyde (MDA). Fluorescence of HMGB1 was significantly increased. We also investigated the effect of the complexes on the expression of genes using RNA-sequence analysis, we further calculated the lowest binding energies between the complexes and proteins using molecular docking. Taken together, the above results indicated that complexes Ir1 and Ir2 induce apoptosis in BEL-7402 cells through a ROS-mediated mitochondrial dysfunction and inhibition of the PI3K/AKT/mTOR signaling pathway.

Identification of a CD4+ conventional T cells-related lncRNAs signature associated with hepatocellular carcinoma prognosis, therapy, and tumor microenvironment

Background

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide, and CD4+ T lymphocytes can inhibit hepatocarcinogenesis and mediate tumor regression. However, few studies have focused on the prognostic power of CD4+ Tconv-related lncRNAs in HCC patients.

Method

We obtained data from TCGA and GEO databases and identified CD4+Tconv-related lncRNAs in HCC. The risk score was constructed using lasso regression and the model was validated using two validation cohorts. The RS was also assessed in different clinical subgroups, and a nomogram was established to further predict the patients' outcomes. Furthermore, we estimated the immune cell infiltration and cancer-associated fibroblasts (CAFs) through TIMER databases and assessed the role of RS in immune checkpoint inhibitors response.

Results

We constructed a CD4+ Tconv-related lncRNAs risk score, including six lncRNAs (AC012073.1, AL031985.3, LINC01060, MKLN1-AS, MSC-AS1, and TMCC1-AS1), and the RS had good predictive ability in validation cohorts and most clinical subgroups. The RS and the T stage were included in the nomogram with optimum prediction and the model had comparable OS prediction power compared to the AJCC. Patients in the high-risk group had a poor immune response phenotype, with high infiltrations of macrophages, CAFs, and low infiltrations of NK cells. Immunotherapy and chemotherapy response analysis indicated that low-risk group patients had good reactions to immune checkpoint inhibitors.

Conclusion

We constructed and validated a novel CD4+ Tconv-related lncRNAs RS, with the potential predictive value of HCC patients' survival and immunotherapy response.

A Model of Basement Membrane-Associated Gene Signature Predicts Liver Hepatocellular Carcinoma Response to Immune Checkpoint Inhibitors

Liver hepatocellular carcinoma (LIHC) is a highly lethal malignant tumor originating from the digestive system, which is a serious threat to human health. In recent years, immunotherapy has shown significant therapeutic effects in the treatment of LIHC, but only for a minority of patients. The basement membrane (BM) plays an important role in the occurrence and development of tumors, including LIHC. Therefore, this study is aimed at establishing a risk score model based on basement membrane-related genes (BMRGs) to predict patient prognosis and response to immunotherapy. First, we defined three patterns of BMRG modification (C1, C2, and C3) by consensus clustering of BMRG sets and LIHC transcriptome data obtained from public databases. Survival analysis showed that patients in the C2 group had a better prognosis, and Gene Set Variation Analysis (GSVA) revealed that the statistically significant pathways were mainly enriched in the C2 group. Moreover, we performed Weighted Correlation Network Analysis (WGCNA) on the above three subgroups and obtained 179 intersecting genes. We further applied function enrichment analyses, and the results demonstrated that they were mainly enriched in metabolism-related pathways. Furthermore, we conducted the LASSO regression analysis and obtained 4 BMRGs (MPV17, GNB1, DHX34, and MAFG) that were significantly related to the prognosis of LIHC patients. We further constructed a prognostic risk score model based on the above genes, which was verified to have good predictive performance for LIHC prognosis. In addition, we analyzed the correlation between the risk score and the tumor immune microenvironment (TIM), and the results showed that the high-risk scoring group tended to be in an immunosuppressed status. Finally, we investigated the relationship between the risk score and LIHC immune function. The results demonstrated that the risk score was closely related to the expression levels of multiple immune checkpoints. Patients in the low-risk group had significantly higher IPS scores, and patients in the high-risk group had lower immune escape and TIDE score. In conclusion, we established a novel risk model based on BMRGs, which may serve as a biomarker for prognosis and immunotherapy in LIHC.

Primary Resistance to Immunotherapy-Based Regimens in First Line Hepatocellular Carcinoma: Perspectives on Jumping the Hurdle

Simple Summary

Immune checkpoint inhibitors (ICIs) had been explored extensively in patients affected by unresectable hepatocellular carcinoma. These agents were expected to be the keystones of the disease’s first-line treatment because they were theoretically able to revert the immune suppressive tumor microenvironment of the cancerous liver, and because of their manageable safety profile. However, when used as monotherapies, they showed important activity and efficacy limitations. In this mini-review, we summarize the characteristics of the different ICIs-based regimens which constitute the present gold standard of first-line treatment, then, moving from their shortcomings, we discuss the rationale supporting the strategies currently under investigation: systemic triplets and new paradigms of immune-therapeutic agents such as CAR-T and vaccines.

Abstract

Immune checkpoint inhibitors (ICIs) are a key component of different stages of hepatocellular carcinoma (HCC) treatment, particularly in the first line of treatment. A lesson on the primary resistance which hampers their efficacy and activity was learned from the failure of the trials which tested them as first-line mono-therapies. Despite the combination of anti-PD(L)1 agents with anti-VEGF, anti CTLA4, or TKIs demonstrating relevant improvements in efficacy, the “doublets strategy” still shows room for improvement, due to a limited overall survival benefit and a high rate of progressive disease as best response. In this review, we discuss the results from the currently tested doublet strategies (i.e., atezolizumab+bevacizumab, durvalumab+tremelimumab with a mention to the newly presented ICIs/TKIs combinations), which highlight the need for therapeutic improvement. Furthermore, we examine the rationale and provide an overview of the ongoing trials testing the treatment intensification strategy with triplet drugs: anti-PD1+anti-CTLA4+anti-VEGF/TKIs and anti-PD1+anti-VEGF+alternative immunity targets. Lastly, we report on the alternative strategy to integrate ICIs into the new paradigm of immune therapeutics constituted by CAR-T and anti-cancer vaccines. This review provides up-to-date knowledge of ongoing clinical trials of the aforementioned strategies and critical insight into their mechanistic premises.