Single-cell analysis on stromal fibroblasts in the microenvironment of solid tumours

Heterogeneity of cancer-associated fibroblast subpopulations in prostate cancer: Implications for prognosis and immunotherapy

BACKGROUND

Prostate cancer stands as the second most common malignancy among men, notorious for its intricate heterogeneity, especially evident in metastatic disease. This complexity presents substantial challenges in treatment efficacy and patient prognosis.

OBJECTIVE

This study endeavors to elucidate the multifaceted roles of cancer-associated fibroblasts within the tumor microenvironment of prostate cancer, with a focus on their implications for disease prognosis and the potential for novel immunotherapeutic strategies.

METHODS

Leveraging advanced single-cell RNA sequencing technology, we meticulously characterized the diverse CAF subpopulations within prostate cancer samples. Our analysis identified four predominant subsets: C0 IER2+, C1 ABCA8+, C2 ABI3BP+, and C3 MEOX2+. We conducted comprehensive gene expression profiling to construct a robust prognostic model reflecting the clinical relevance of these subpopulations.

RESULTS

C1 ABCA8+ fibroblasts demonstrated heightened proliferative activity, underscoring their pivotal role in fostering tumor growth and metastasis via intricate signaling pathways. In vitro experiments verified that the T transcription factor NFAT5 of C1 ABCA8+ fibroblasts subpopulation was knocked down in LNCaP clone FGC and 22Rv1 cell lines, which was closely related to the proliferation of PC. Moreover, we identified key genes linked to patient outcomes and immune landscape alterations, reinforcing the prognostic significance of CAF characteristics in this context.

CONCLUSION

This investigation illuminates the critical potential of targeting CAFs to augment immunotherapeutic approaches in prostate cancer. Our findings contribute to a deeper understanding of the TME's complexity, advocating for further exploration into CAF-targeted therapies aimed at enhancing treatment responses and ultimately improving patient outcomes.

Desmoglein-2 was a novel cancer-associated fibroblasts-related biomarker for oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is the most common type of oral cancer with alarmingly high morbidity. The cancer-associated fibroblasts (CAFs) play a pivotal role in tumor development, while their specific mechanisms in OSCC remains largely unclear. Our object is to explore a CAFs-related biomarker in OSCC.

METHODS

Single-cell RNA sequencing (ScRNA-seq) analysis was used to pinpoint CAF clusters in OSCC samples. Differentially expressed genes and Cox regression analyses were used to identify candidate genes, and their functions were evaluated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. The prognostic performance of the identified biomarker was evaluated using receiver operating characteristic analysis. The qPCR and western blot were used to assess gene expression. The hub gene related immune characteristics were analyzed in independent cohorts, and gene expression differences between different immunotherapy response groups were investigated using Pearson correlation analysis.

RESULTS

Desmoglein-2 (DSG2) was identified as a CAFs-related biomarker in OSCC exhibiting elevated expression compared to controls and being associated with poor prognosis. Enrichment analyses revealed that DSG2 was involved in signal transduction pathways like focal adhesion. The Area Under Curve values of DSG2 in predicting prognosis exceeded 0.6 in both training-set and validation-set. Furthermore, patients with low DSG2 expression were more likely to benefit from immunotherapy than those DSG2 highly expressed patients.

CONCLUSION

Our study identified DSG2 as a reliable CAFs-related prognostic biomarker in OSCC, providing a new reference for the mechanistic understanding and target therapy of this malignancy.

Revealing a cancer-associated fibroblast-based risk signature for pancreatic adenocarcinoma through single-cell and bulk RNA-seq analysis

PURPOSE

Proliferation of stromal connective tissue is a hallmark of pancreatic adenocarcinoma (PAAD). The engagement of activated cancer-associated fibroblasts (CAFs) contributes to the progression of PAAD through their involvement in tumor fibrogenesis. However, the prognostic significance of CAF-based risk signature in PAAD has not been explored.

METHODS

The single-cell RNA sequencing (scRNA-seq) data sourced from GSE155698 within the Gene Expression Omnibus (GEO) database was supplemented by bulk RNA sequencing data from The Cancer Genome Atlas (TCGA) and microarray data retrieved from the GEO database. The scRNA-seq data underwent processing via the Seurat package to identify distinct CAF clusters utilizing specific CAF markers. Differential gene expression analysis between normal and tumor samples was conducted within the TCGA-PAAD cohort. Univariate Cox regression analysis pinpointed genes associated with CAF clusters, identifying prognostic CAF-related genes. These genes were utilized in LASSO regression to craft a predictive risk signature. Subsequently, integrating clinicopathological traits and the risk signature, a nomogram model was constructed.

RESULTS

Our scRNA-seq analysis unveiled four distinct CAF clusters in PAAD, with two linked to PAAD prognosis. Among 207 identified DEGs, 148 exhibited significant correlation with these CAF clusters, forming the basis of a seven-gene risk signature. This signature emerged as an independent predictor in multivariate analysis for PAAD and demonstrated predictive efficacy in immunotherapeutic outcomes. Additionally, a novel nomogram, integrating age and the CAF-based risk signature, exhibited robust predictability and reliability in prognosticating PAAD. Moreover, the risk signature displayed substantial correlations with stromal and immune scores, as well as specific immune cell types.

CONCLUSIONS

The prognosis of PAAD can be accurately predicted using the CAF-based risk signature, and a thorough analysis of the PAAD CAF signature may aid in deciphering the patient's immunotherapy response and presenting fresh cancer treatment options.

Predicting Survival Signature of Bladder Cancer Related to Cancer-Associated Fibroblast (CAF) Constructed by Intersecting Genes in TCGA and GEO

Bladder cancer was one of the most common carcinomas around the world. However, the mechanism of the disease still remained to be investigated. We expected to establish a prognostic survival model with 9 prognostic genes to predict overall survival (OS) in patients of bladder cancer. The gene expression data of bladder cancer were obtained from TCGA and GEO datasets. TCGA and GEO datasets were used for screening prognostic genes along with developing and validating a 9-gene prognostic survival model by method of weighted gene co-expression network analysis (WGCNA) and LASSO with Cox regression. The relative analysis of evaluate tumor burden mutation (TBM), GO, KEGG, chemotherapy drug and functional pathway were also performed based on CAF-related mRNAs. 151 Overlapping CAF-related genes were distinguished after intersecting differentially expressed genes from 945 genes in TCGA and 491 genes in GEO dataset. 9 Prognostic genes (MSRB2, AGMAT, KLF6, DDAH2, GADD45B, SERPINE2, STMN3, TEAD2, and COMP) were used for construction of prognostic model after LASSO with Cox regression. Receiver operating characteristic (ROC) curves showed a good survival prediction by this model. Functional analysis indicated chemokine, cytokine, ECM interaction, oxidative stress and apoptosis were highly appeared. Potential drugs targeted different CAF-related genes like vemurafenib, irofulven, ghiotepa, and idarubicin were found as well. We constructed a novel 9 CAF-related mRNAs prognostic model and explored the gene expression and potential functional information of related genes, which might be worthy of clinical application. In addition, potential chemotherapy drugs could provide useful insights into the potential clinical treatment of bladder cancer.

Single-cell transcriptome sequencing reveals spatial distribution of IL34+ cancer-associated fibroblasts in hepatocellular carcinoma tumor microenvironment

We utilized scRNA-seq, a well-established technology, to uncover the gene expression characteristics of IL34+ CAFs within HCC. We analyzed the related mechanisms through in vitro and in vivo assays. To begin, we acquired scRNA-seq datasets about HCC, which enabled us to identify distinct cell subpopulations within HCC tissues. We conducted a differential analysis to pinpoint DEGs associated with normal fibroblasts (NFs) and CAFs. Subsequently, we isolated NFs and CAFs, followed by the sorting of IL34+ CAFs. These IL34+ CAFs were then co-cultured with T cells and HCC cells to investigate their potential role in Tregs infiltration, CD8+ T cell toxicity, and the biological processes of HCC cells. We validated our findings in vivo using a well-established mouse model. Our analysis of HCC tissues revealed the presence of seven primary cell subpopulations, with the most significant disparities observed within fibroblast subpopulations. Notably, high IL34 expression was linked to increased expression of receptor proteins and enhanced proliferative activity within CAFs, with specific expression in CAFs. Furthermore, we identified a substantial positive correlation between IL34 expression and the abundance of Tregs. Both in vitro and in vivo experiments demonstrated that IL34+ CAFs promoted Tregs infiltration while suppressing CD8+ T cell toxicity. Consequently, this promoted the growth and metastasis of HCC. In summary, our study affirms that IL34+ CAFs play a pivotal role in augmenting the proliferative activity of CAFs, facilitating Tregs infiltration, and inhibiting CD8+ T cell toxicity, ultimately fostering the growth and metastasis of HCC.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq unravels the heterogeneity of cancer-associated fibroblasts in TNBC

BACKGROUND

The treatment of triple-negative breast cancer (TNBC) is one of the main focuses and key difficulties because of its heterogeneity, and the source of this heterogeneity is unclear.

METHODS

Single-cell RNA (scRNA) and transcriptomics data of TNBC and normal breast samples were retrieved from Gene Expression Omnibus (GEO) database and TCGA-BRCA database. These cells were clustered using the t-SNE and UMAP method, and the marker genes for each cluster were found. We annotated the clusters using the published literature, CellMarker database and "SingleR" R package.

RESULTS

A total of 1535 cells and 21785 genes from 6 TNBC patients and 2068 cells and 15868 genes from 3 normal breast tissues were used for downstream analyses. The scRNA data were divided into 14 clusters labeled into 8 cell types, including epithelial cells, immunocytes, CAFs/fibroblasts and etc. In the TNBC samples, CAFs were divided into three clusters and labelled as prCAFs, myCAFs and emCAFs, and the marker genes were DCN, FAP and RGS5, respectively. The prCAF subgroup is functionally characterized by promoting proliferation and multi drug resistance; myCAF subgroup is involved in constituting the extracellular matrix and collagen production, matrix composition and collagen production, and the emCAF functionally characterized by energy metabolism.

CONCLUSIONS

TNBC has inter- and intra-tumor heterogeneity, and CAF is one of the sources of this heterogeneity. CD74, SASH3, CD2, TAGAP and CCR7 served as significant marker genes with prognostic and therapeutic value.

Effect of activating cancer-associated fibroblasts biomarker TNC on immune cell infiltration and prognosis in breast cancer

BACKGROUND

Cancer-associated fibroblasts (CAFs) are the most important components of the tumor microenvironment (TME). CAFs are heterogeneous and involved in tumor tumorigenesis and drug resistance, contributing to TME remodeling and predicting clinical outcomes as prognostic factors. However, the effect of CAFs the TME and the prognosis of patients with breast cancer (BC) is not fully understood. This study investigated the correlation between CAFs-activating biomarkers immune cell infiltration and survival in patients with breast cancer.

METHODS

RNA sequencing data and survival information for patients with breast cancer were downloaded from The Cancer Genome Atlas (TCGA) using R software. We then analyzed the correlation between CAFs-expressing biomarkers and immune cells using the clusterProfiler package, and evaluated the prognostic role of appealing genes using the Survminer package. Immunohistochemical (IHC) staining was used to determine the expression levels of TNC in 160 breast cancer samples pathologically diagnosed as invasive ductal carcinoma that were not otherwise specified (IDC-NOS).

RESULTS

Data analysis showed that CAFs-expressing genes was higher than in normal tissues (p < 0.05). Pathway enrichment revealed that the overexpression of CAFs-related genes was mainly enriched in the focal adhesion and phosphoinositol-3 kinase-serine/threonine kinase (PI3K-AKT) signaling pathways. Immune infiltration analysis suggested that high expression of CAFs-related genes was significantly positively correlated with the infiltration of naive B cells and resting dendritic cells and inversely correlated with macrophages cell infiltration. In addition, high TNC expression in tumor cells was associated with the most adverse clinicopathological features and reduced metastasis-free survival (MFS) (hazard ratio (HR) 0.574, 95% confidence interval (CI) 0.404-0.815, p = 0.035).

CONCLUSIONS

This study found that CAFs may participate in immunosuppression and regulate tumor cell proliferation and invasion. High TNC expression is associated with several adverse clinicopathological features, and high TNC expression in tumor cells has been identified as an independent prognostic factor for IDC-NOS.

A fibroblast-associated signature predicts prognosis and immunotherapy in esophageal squamous cell cancer

Background

Current paradigms of anti-tumor therapies are not qualified to evacuate the malignancy ascribing to cancer stroma's functions in accelerating tumor relapse and therapeutic resistance. Cancer-associated fibroblasts (CAFs) has been identified significantly correlated with tumor progression and therapy resistance. Thus, we aimed to probe into the CAFs characteristics in esophageal squamous cancer (ESCC) and construct a risk signature based on CAFs to predict the prognosis of ESCC patients.

Methods

The GEO database provided the single-cell RNA sequencing (scRNA-seq) data. The GEO and TCGA databases were used to obtain bulk RNA-seq data and microarray data of ESCC, respectively. CAF clusters were identified from the scRNA-seq data using the Seurat R package. CAF-related prognostic genes were subsequently identified using univariate Cox regression analysis. A risk signature based on CAF-related prognostic genes was constructed using Lasso regression. Then, a nomogram model based on clinicopathological characteristics and the risk signature was developed. Consensus clustering was conducted to explore the heterogeneity of ESCC. Finally, PCR was utilized to validate the functions that hub genes play on ESCC.

Results

Six CAF clusters were identified in ESCC based on scRNA-seq data, three of which had prognostic associations. A total of 642 genes were found to be significantly correlated with CAF clusters from a pool of 17080 DEGs, and 9 genes were selected to generate a risk signature, which were mainly involved in 10 pathways such as NRF1, MYC, and TGF-Beta. The risk signature was significantly correlated with stromal and immune scores, as well as some immune cells. Multivariate analysis demonstrated that the risk signature was an independent prognostic factor for ESCC, and its potential in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the CAF-based risk signature and clinical stage was developed, which exhibited favorable predictability and reliability for ESCC prognosis prediction. The consensus clustering analysis further confirmed the heterogeneity of ESCC.

Conclusion

The prognosis of ESCC can be effectively predicted by CAF-based risk signatures, and a comprehensive characterization of the CAF signature of ESCC may aid in interpreting the response of ESCC to immunotherapy and offer new strategies for cancer treatment.

Single-cell transcriptome analysis for cancer and biology of the pancreas: A review on recent progress

Single-cell sequencing has become one of the most used techniques across the wide field of biology. It has enabled researchers to investigate the whole transcriptome at the cellular level across tissues, which unlocks numerous potentials for basic and applied studies in future diagnosis and therapy. Here, we review the impact of single-cell RNA sequencing, as the prominent single-cell technique, in pancreatic biology and cancer. We discuss the most recent findings about pancreatic physiology and pathophysiology owing to this technological advancement in the past few years. Using single-cell RNA sequencing, researchers have been able to discover cellular heterogeneity across healthy cell types, as well as cancer tissues of the pancreas. We will discuss the new immunological targets and new molecular mechanisms of progression in the microenvironment of pancreatic cancer studied using single-cell RNA sequencing. The scope is not limited to cancer tissues, and we cover novel developmental, evolutionary, physiological, and heterogenic insights that have also been achieved recently for pancreatic tissues. We cover all biological insights derived from the single-cell RNA sequencing data, discuss the corresponding pros and cons, and finally, conclude how future research can move better by utilizing single-cell analysis for pancreatic biology.

Characterization of cancer-related fibroblasts (CAF) in hepatocellular carcinoma and construction of CAF-based risk signature based on single-cell RNA-seq and bulk RNA-seq data

Background

Cancer-associated fibroblasts (CAFs) are involved in tumor growth, angiogenesis, metastasis, and resistance to therapy. We sought to explore the CAFs characteristics in hepatocellular carcinoma (HCC) and establish a CAF-based risk signature for predicting the prognosis of HCC patients.

Methods

The signal-cell RNA sequencing (scRNA-seq) data was obtained from the GEO database. Bulk RNA-seq data and microarray data of HCC were obtained from the TCGA and GEO databases respectively. Seurat R package was applied to process scRNA-seq data and identify CAF clusters according to the CAF markers. Differential expression analysis was performed to screen differentially expressed genes (DEGs) between normal and tumor samples in TCGA dataset. Then Pearson correlation analysis was used to determine the DEGs associated with CAF clusters, followed by the univariate Cox regression analysis to identify CAF-related prognostic genes. Lasso regression was implemented to construct a risk signature based on CAF-related prognostic genes. Finally, a nomogram model based on the risk signature and clinicopathological characteristics was developed.

Results

Based on scRNA-seq data, we identified 4 CAF clusters in HCC, 3 of which were associated with prognosis in HCC. A total of 423 genes were identified from 2811 DEGs to be significantly correlated with CAF clusters, and were narrowed down to generate a risk signature with 6 genes. These six genes were primarily connected with 39 pathways, such as angiogenesis, apoptosis, and hypoxia. Meanwhile, the risk signature was significantly associated with stromal and immune scores, as well as some immune cells. Multivariate analysis revealed that risk signature was an independent prognostic factor for HCC, and its value in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the stage and CAF-based risk signature was constructed, which exhibited favorable predictability and reliability in the prognosis prediction of HCC.

Conclusion

CAF-based risk signatures can effectively predict the prognosis of HCC, and comprehensive characterization of the CAF signature of HCC may help to interpret the response of HCC to immunotherapy and provide new strategies for cancer treatment.

Engineering stromal heterogeneity in cancer

Based on our exponentially increasing knowledge of stromal heterogeneity from advances in single-cell technologies, the notion that stromal cell types exist as a spectrum of unique subpopulations that have specific functions and spatial distributions in the tumor microenvironment has significant impact on tumor modeling for drug development and personalized drug testing. In this Review, we discuss the importance of incorporating stromal heterogeneity and tumor architecture, and propose an overall approach to guide the reconstruction of stromal heterogeneity in vitro for tumor modeling. These next-generation tumor models may support the development of more precise drugs targeting specific stromal cell subpopulations, as well as enable improved recapitulation of patient tumors in vitro for personalized drug testing.

Targeting the Stroma in the Management of Pancreatic Cancer

Pancreatic cancer (PC) presents extremely aggressive tumours and is associated with poor survival. This is attributed to the unique features of the tumour microenvironment (TME), which is known to create a dense stromal formation and poorly immunogenic condition. In particular, the TME of PC, including the stromal cells and extracellular matrix, plays an essential role in the progression and chemoresistance of PC. Consequently, several promising agents that target key components of the stroma have already been developed and are currently in multiple stages of clinical trials. Therefore, the authors review the latest available evidence on novel stroma-targeting approaches, highlighting the potential impact of the stroma as a key component of the TME in PC.

Single-Cell Transcriptomic Analysis of Tumor-Derived Fibroblasts and Normal Tissue-Resident Fibroblasts Reveals Fibroblast Heterogeneity in Breast Cancer

Cancer-associated fibroblasts (CAFs) are a prominent stromal cell type in solid tumors and molecules secreted by CAFs play an important role in tumor progression and metastasis. CAFs coexist as heterogeneous populations with potentially different biological functions. Although CAFs are a major component of the breast cancer stroma, molecular and phenotypic heterogeneity of CAFs in breast cancer is poorly understood. In this study, we investigated CAF heterogeneity in triple-negative breast cancer (TNBC) using a syngeneic mouse model, BALB/c-derived 4T1 mammary tumors. Using single-cell RNA sequencing (scRNA-seq), we identified six CAF subpopulations in 4T1 tumors including: 1) myofibroblastic CAFs, enriched for α-smooth muscle actin and several other contractile proteins; 2) ‘inflammatory’ CAFs with elevated expression of inflammatory cytokines; and 3) a CAF subpopulation expressing major histocompatibility complex (MHC) class II proteins that are generally expressed in antigen-presenting cells. Comparison of 4T1-derived CAFs to CAFs from pancreatic cancer revealed that these three CAF subpopulations exist in both tumor types. Interestingly, cells with inflammatory and MHC class II-expressing CAF profiles were also detected in normal breast/pancreas tissue, suggesting that these phenotypes are not tumor microenvironment-induced. This work enhances our understanding of CAF heterogeneity, and specifically targeting these CAF subpopulations could be an effective therapeutic approach for treating highly aggressive TNBCs.