Single-cell RNA-seq reveals invasive trajectory and determines cancer stem cell-related prognostic genes in pancreatic cancer

Identification and validation of stemness-based and ferroptosis-related molecular clusters in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy with an extremely poor prognosis. Cancer stem cells (CSCs) are considered to be responsible for the poor survival, recurrence and therapy resistance of PDAC. Ferroptosis plays a crucial role in the sustain and survival of CSCs. Here, we employed a rigorous evaluation of multiple datasets to identify a novel stemness-based and ferroptosis-related genes (SFRGs) signature to access the potential prognostic application. This work we retrieved RNA-sequencing and clinical annotation data from the TCGA, ICGC, GTEx and GEO database, and acquired 26 stem cell gene sets and 259 ferroptosis genes from StemChecker database and FerrDb database, respectively. Based on consensus clustering and ssGSEA analysis, we identified two expression patterns of CSCs traits (C1 and C2). Then, WGCNA analysis was implemented to screen out hub module genes correlated with stemness. Furthermore, differential expression analysis, Pearson correlation analysis, and the Least absolute shrinkage and selection operator (LASSO) and Cox regression were performed to identify the SFRGs and to construct model. In addition, the differences in prognosis, tumor microenvironment (TME) components and therapy responses were evaluated between two risk groups. Finally, we verified the most influential marker ARNTL2 experimentally by western blot, qRT-PCR, sphere formation assay, mitoscreen assay, intracellular iron concentration determination and MDA determination assays. In conclusion, we developed a stemness-based and ferroptosis-related prognostic model, which could help predict overall survival for PDAC patients. Targeting ferroptosis may be a promising therapeutic strategy to inhibit PDAC progression by suppressing CSCs.

Comprehensive profiling of cell subsets of gastric cancer and liver metastasis based on single cell RNA-sequencing analysis

Background

Liver metastasis (Li) is one of the most common distant metastatic sites for gastric cancer. A deeper understanding of its mechanism of action from a bioinformatics perspective may provide new insights. Therefore, the aim of this study was to use single cell RNA sequencing (scRNA-seq) to evaluate cell subtypes and understand the molecular mechanism of gastric cancer Li.

Methods

The scRNA-seq data GSE163558 of gastric cancer and Li were downloaded from Gene Expression Omnibus (GEO). Single cell data were analyzed by various R packages such as Seurat, CellChat, gene set variation analysis (GSVA), monocle, gene set enrichment analysis (GSEA), and survival analysis and the results were plotted by ggplot2, R4.1.1. Protein expression was confirmed by immunohistochemistry in an additional patient cohort.

Results

The genes APOD, CXCL5, and JUN are involved in epithelial cell metastasis. The infiltration of cytotoxic CD8 T cells was more frequent in gastric primary tumors (PTs) than in Lis. IL7R high natural killer (NK) cells that had high TXNIP and RIPOR2 expression were located at the site of Li and corresponded to a favorable prognosis. NK cells with high TNFAIP3 expression were located at the PT site and corresponded to a poor prognosis. Furthermore, the gene expression of myeloid cells was different depending on their localization in the PT site or Li. MHC-I signaling pathway was found to be activated in the PT compared to MHC-II at the site of Li, as revealed by cell-cell interaction, and HLA-E-CD94/NKG2A of NK cells was only activated in the PT and not in the Li.

Conclusions

The present study revealed the difference between Li and gastric PT by scRNA-seq, demonstrating the impact of partial gene expression on patient prognosis. Our study provides new insights into gastric cancer and Li.

Identification of ONECUT3 as a stemness-related transcription factor regulating NK cell-mediated immune evasion in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) has a dismal response to the current T cell-based immunotherapies, which is attributed to intratumoral heterogeneity caused by PDAC stem cells and lack of major histocompatibility complex class I required for neoantigen presentation. Although this scenario makes natural killer (NK) cells attractive candidates for immunotherapeutic agents targeting MHC-I-deficient cancer stem cells in heterogeneous PDACs, little is known about PDAC stem cell immunology. In our study, PDAC-specific datasets from public databases were collected for in-depth bioinformatic analysis. We found that the abundance of PDAC stemness negatively influenced the infiltration of NK cells and identified the transcription factor ONECUT3 enriched in PDACs with high stemness index scores and Pan-cancer Stemness Signature levels. A series of NK cell-targeted inhibitory immune checkpoints were highly expressed in ONECUT3high PDACs. The patient group with high levels of ONECUT3 expression had a high risk of poor overall survival, even if accompanied by high infiltration of NK cells. Furthermore, the prostanoid metabolic process was enriched in ONECUT3high PDACs with high levels of NK cell-targeted inhibitory immune checkpoints. ONECUT3 enriched in high-stemness PDACs possessed the potential to transcriptionally regulate the prostanoid metabolism-related genes. Our study reveals ONECUT3 as a candidate stemness-related transcription factor regulating NK cell-targeted inhibitory immune checkpoints in PDAC. ONECUT3-mediated prostanoid metabolism may regulate cancer stemness and immune evasion in PDAC. Synergistic inhibition of prostanoid metabolism may improve the efficacy of NK cell-based immunotherapies targeting intratumoral heterogeneity caused by PDAC stem cells.

Exosomes: Another intercellular lipometabolic communication mediators in digestive system neoplasms?

Neoplasms are one of the most concerned public health problems worldwide. Digestive system neoplasms, with a high morbidity and mortality, is one of the most common malignant tumors in human being. It is found that exosomes act as an intercellular communication media to carry the metabolic and genetic information of parental cells to target cells. Likely, exosomes participate in lipid metabolism and regulates multiple processes in digestive system neoplasms, including the information transmission among cancer cells, the formation of neoplastic microenvironment, and the neoplastic biological behaviors like metastasis, invasion, and the chemotherapy resistance. In this review, we firstly introduce the main mechanisms whereas exosomes act as intercellular lipometabolic communication mediator in digestive system neoplasms. Thereafter we introduce the relationship between exosomes lipid metabolism and various type of digestive system neoplasms, including gastric cancer, hepatocellular carcinoma, pancreatic cancer, and colorectal cancer. Eventually, we summarized and prospected the development and implication of exosomes in digestive system neoplasms. The further research of exosomes as intercellular lipid metabolism mediator will contribute to accurate and efficient diagnosis and treatment of digestive system neoplasms.

Role of RNA methylation in the regulation of pancreatic cancer stem cells (Review)

Pancreatic cancer stem cells (CSCs) play a key role in the initiation and progression of pancreatic adenocarcinoma (PDAC). CSCs are responsible for resistance to chemotherapy and radiation, and for cancer metastasis. Recent studies have indicated that RNA methylation, a type of RNA modification, predominantly occurring as m6A methylation, plays an important role in controlling the stemness of cancer cells, therapeutic resistance against chemotherapy and radiation therapy, and their overall relevance to a patient's prognosis. CSCs regulate various behaviors of cancer through cell-cell communication by secreting factors, through their receptors, and through signal transduction. Recent studies have shown that RNA methylation is involved in the biology of the heterogeneity of PDAC. The present review provides an update on the current understanding of RNA modification-based therapeutic targets against deleterious PDAC. Several key pathways and agents that can specifically target CSCs have been identified, thus providing novel insights into the early diagnosis and efficient treatment of PDAC.

Single-cell omics: a new perspective for early detection of pancreatic cancer?

Pancreatic cancer is one of the most lethal cancers, mostly due to late diagnosis and limited treatment options. Early detection of pancreatic cancer in high-risk populations bears the potential to greatly improve outcomes, but current screening approaches remain of limited value despite recent technological advances. This review explores the possible advantages of liquid biopsies for this application, particularly focusing on circulating tumour cells (CTCs) and their subsequent single-cell omics analysis. Originating from both primary and metastatic tumour sites, CTCs provide important information for diagnosis, prognosis and tailoring of treatment strategies. Notably, CTCs have even been detected in the blood of subjects with pancreatic precursor lesions, suggesting their suitability as a non-invasive tool for the early detection of malignant transformation in the pancreas. As intact cells, CTCs offer comprehensive genomic, transcriptomic, epigenetic and proteomic information that can be explored using rapidly developing techniques for analysing individual cells at the molecular level. Studying CTCs during serial sampling and at single-cell resolution will help to dissect tumour heterogeneity for individual patients and among different patients, providing new insights into cancer evolution during disease progression and in response to treatment. Using CTCs for non-invasive tracking of cancer features, including stemness, metastatic potential and expression of immune targets, provides important and readily accessible molecular insights. Finally, the emerging technology of ex vivo culturing of CTCs could create new opportunities to study the functionality of individual cancers at any stage and develop personalised and more effective treatment approaches for this lethal disease.

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.

Identification of CEACAM5 as a stemness-related inhibitory immune checkpoint in pancreatic cancer

BACKGROUND

Immunotherapy has emerged as a new cancer treatment modality. However, tumour heterogeneity can diminish checkpoint blockade response and shorten patient survival. As a source of tumour heterogeneity, cancer stem cells act as an indispensable reservoir for local recurrence and distant metastasis. Thus, precision immunotherapy targeting tumour heterogeneity requires a comprehensive understanding of cancer stem cell immunology. Our study aimed to identify stemness-related inhibitory immune checkpoints and relevant regulatory pathways in pancreatic cancer.

METHODS

Pancreatic cancer-specific datasets in The Cancer Genome Atlas and the Cancer Therapeutics Response Portal were collected for in-depth bioinformatic analysis. Differentially expressed genes between pancreatic cancers with high and low stemness index (mRNAsi) scores were compared to screen out inhibitory immune checkpoints. Survival analysis was used to predict the prognostic value of immune checkpoint plus immune infiltrate in patients with pancreatic cancer. The expression of stemness-related immune checkpoint across immune subtypes of pancreatic cancer was detected and gene set enrichment analysis was performed to figure out the relevant regulatory signallings.

RESULTS

The abundance of cancer stemness predicted a low immunotherapy response to pancreatic cancer. The inhibitory immune checkpoint CEACAM5 that was enriched in pancreatic cancers with high mRNAsi scores also exhibited a strong correlation with invasive cell-enriched signature and Msi⁺ tumour-initiating cell-enriched signature. Levels of CEACAM5 expression were higher in the interferon-γ dominant immune subtype of pancreatic cancers that are characterized by high M1 macrophage infiltration. The patient group with high levels of CEACAM5 expression had a high risk of poor overall survival, even if accompanied by high infiltration of M1 macrophages. Furthermore, prostanoid and long-chain unsaturated fatty acid metabolic processes showed a significant association with cancer stemness and CEACAM5 expression.

CONCLUSIONS

Our findings suggest that CEACAM5 is a candidate stemness-related innate immune checkpoint in pancreatic cancer, and is potentially regulated by prostanoid and long-chain unsaturated fatty acid metabolic processes. Immune checkpoint blockade of CEACAM5, which synergizes with inhibition of those regulatory pathways, may improve the efficacy of precision immunotherapy targeting tumour heterogeneity caused by cancer stem cells.

Estimating Metastatic Risk of Pancreatic Ductal Adenocarcinoma at Single-Cell Resolution

Pancreatic ductal adenocarcinoma (PDAC) is characterized by intra-tumoral heterogeneity, and patients are always diagnosed after metastasis. Thus, finding out how to effectively estimate metastatic risk underlying PDAC is necessary. In this study, we proposed scMetR to evaluate the metastatic risk of tumor cells based on single-cell RNA sequencing (scRNA-seq) data. First, we identified diverse cell types, including tumor cells and other cell types. Next, we grouped tumor cells into three sub-populations according to scMetR score, including metastasis-featuring tumor cells (MFTC), transitional metastatic tumor cells (TransMTC), and conventional tumor cells (ConvTC). We identified metastatic signature genes (MSGs) through comparing MFTC and ConvTC. Functional enrichment analysis showed that up-regulated MSGs were enriched in multiple metastasis-associated pathways. We also found that patients with high expression of up-regulated MSGs had worse prognosis. Spatial mapping of MFTC showed that they are preferentially located in the cancer and duct epithelium region, which was enriched with the ductal cells' associated inflammation. Further, we inferred cell-cell interactions, and observed that interactions of the ADGRE5 signaling pathway, which is associated with metastasis, were increased in MFTC compared to other tumor sub-populations. Finally, we predicted 12 candidate drugs that had the potential to reverse expression of MSGs. Taken together, we have proposed scMetR to estimate metastatic risk in PDAC patients at single-cell resolution which might facilitate the dissection of tumor heterogeneity.

Investigate the stemness of adult adipose-derived stromal cells based on single-cell RNA-sequencing

The cellular heterogeneity and genetic features of stemness of adipose-derived stromal cells (ADSCs) remain unclear. Using single-cell RNA sequencing (scRNA-seq), we investigated the genomic features of the stemness gene in ADSCs with genetic variability. We cultured the ADSCs isolated from the fat waste of a healthy adult volunteers undergoing cosmetic plastic surgery to the third generation, used the BD Rhapsody platform to perform scRNA-seq, then used Monocle2 to analyze the growth and development trajectory of ADSCs, Cellular Trajectory Reconstruction Analysis Using Gene Counts and Expression (CytoTRACE) to evaluate the stemness gene characteristics in ADSCs clusters, and Beam to analyze the expression change characteristics of the main stemness related genes of ADSCs. According to the scRNA-seq data of 5325 ADSCs, they could be classified into nine cell clusters. According to CytoTRACE analysis, Cluster 3 of ADSCs had the highest stemness, whereas Cluster 8 had the lowest stemness. Pseudotime analysis revealed that Cluster 3 of ADSCs was primarily dispersed in the middle part of the growth and development trajectory, whereas Cluster 8 was primarily distributed at the end. We summarized the stemness of Cluster 3 in ADSCs with high expression of TPM1 and CCND1 genes in the metaphase of growth and development is the strongest, whereas the stemness of Cluster 8 with high expression of FICD, CREBRF, SDF2L1, HERPUD1, and HYOU1 genes in the telophase of growth and development is the weakest, providing a theoretical basis for screening and improving the therapeutic effect of ADSCs in cell transplantation research.

Single Cell RNA Sequencing: A New Frontier in Pancreatic Ductal Adenocarcinoma

Simple Summary

Pancreatic cancer has a very low survival rate for several reasons. One of those is primarily due to the difficulty in diagnosing it at an early stage. For this reason, it is important to refine our understanding of this disease to guide the development of new diagnostic and therapeutic modalities to combat this fatal illness. Here we attempt to provide a review of current progress in utilizing single-cell RNA sequencing (scRNA-seq) techniques in the molecular profiling of pancreatic ductal adenocarcinoma.

Abstract

Pancreatic ductal adenocarcinoma is a malignancy with a high mortality rate. It exhibits significant heterogeneity in metabolic pathways which are associated with its progression. In this review, we discuss the role of single cell RNA sequencing in unraveling the metabolic and clinical features of these highly malignant tumors.

The clinical significance of circulating tumor cells and T lymphocyte subtypes in pancreatic cancer patients

Circulating tumor cells (CTCs) are sensitive and reliable biomarkers for tracing relapsed and metastatic cancer. Here, we explore the clinical significance of CTCs and T lymphocyte subtypes in patients with pancreatic cancer. A total of 106 patients with the pancreatic cancer were enrolled in this study. The enrichment and identification of CTCs were achieved before treatment by a PatrolCTC detection technique. Flow cytometry (FACS) was used to characterize CD4, CD8, natural killer (NK) cells, and Tregulatory (Treg) lymphocyte subtypes. Interleukin-2 (IL-2), Interleukin-4 (IL-4), Interleukin-17A (IL-17A), Interleukin-10 (IL-10), and Interferon γ (IFN-γ) were measured by meso-scale discovery (MSD) assay. Among these patients, 44 (41.5%) patients with pancreatic ductal adenocarcinoma (PDAC) were female and 62 (58.5%) cases were male. Case numbers with II-IV tumor-node-metastasis (TNM) stages were 32 (30.2%), 50 (47.2%), and 24 (22.6%), respectively. The positive rate of CTCs before surgery was 37.5% (12/32), 88.0% (44/50) and 100% (24/24) in stage II, III, and IV patients, respectively. Total CTCs, mixed CTCs, and mesenchymal CTCs (MCTCs) were strongly relevant to shorter progression-free survival (PFS) of the patients. In addition, total CTCs (≥6) and positive MCTCs were also significantly correlated with recurrence and metastasis. The patients with high CTCs also had low levels of CD4, CD4/CD8 ratio, NK cells, IL-2, and IFNγ. In contrast, Treg cells had significant elevation in PDAC patients. These results indicated that CTCs number in PDAC patients was an independent indicator for worse PFS. High CTCs also had strong correlation with weak cellular immunity functions.