Single-cell RNA-seq reveals dynamic change in tumor microenvironment during pancreatic ductal adenocarcinoma malignant progression.

Pancreatic ductal adenocarcinoma (PDAC) is most aggressive among all gastrointestinal tumors. The complex intra-tumor heterogeneity and special tumor microenvironment in PDAC bring great challenges for developing effective treatment strategies. We aimed to delineate dynamic changes of tumor microenvironment components during PDAC malignant progression utilizing single-cell RNA sequencing. A total of 11 samples (4 PDAC I, 4 PDAC II, 3 PDAC III) were used to construct expression matrix. After identifying distinct cell clusters, subcluster analysis for each cluster was performed. New cancer associated fibroblasts (CAFs) subset was validated by weighted gene co-expression network analysis, RNA in situ hybridization and immunofluorescence. We found that ductal cells were not dominant component while tumor infiltrating immune cells and pancreatic stellate cells gradually accumulated during tumor development. We defined several new Treg and exhausted T cell signature genes, including DUSP4, FANK1 and LAIR2. The analysis of TCGA datasets showed that patients with high expression of DUSP4 had significantly worse prognosis. In addition, we identified a new CAFs subset (complement-secreting CAFs, csCAFs), which specifically expresses complement system components, and constructed csCAFs-related module by weighted gene co-expression network analysis. The csCAFs were located in the tissue stroma adjacent to malignant ductal cells only in early PDAC. We systematically explored PDAC heterogeneity and identified csCAFs as a new CAFs subset special to PDAC, which may be valuable for understanding the crosstalk inside tumor. This study was supported by The Natural Science Foundation of China (NO.81572339, 81672353, 81871954) and the Youth Clinical Research Project of Peking University First Hospital (2018CR28).

Chen K ，Wang Q ，Li M ，Guo H ，Liu W ，Wang F ，Tian X ，Yang Y ... -  《EBioMedicine》

Deciphering the Prognostic Implications of the Components and Signatures in the Immune Microenvironment of Pancreatic Ductal Adenocarcinoma.

Background: The treatment modalities for pancreatic ductal adenocarcinoma (PDAC) are limited and unsatisfactory. Although many novel drugs targeting the tumor microenvironment, such as immune checkpoint inhibitors, have shown promising efficacy for some tumors, few of them significantly prolong the survival of patients with PDAC due to insufficient knowledge on the tumor microenvironment. Methods: A single-cell RNA sequencing (scRNA-seq) dataset and seven PDAC cohorts with complete clinical and bulk sequencing data were collected for bioinformatics analysis. The relative proportions of each cell type were estimated using the gene set variation analysis (GSVA) algorithm based on the signatures identified by scRNA-seq or previous literature. Results: A meta-analysis of 883 PDAC patients showed that neutrophils are associated with worse overall survival (OS) for PDAC, while CD8+ T cells, CD4+ T cells, and B cells are related to prolonged OS for PDAC, with marginal statistical significance. Seventeen cell categories were identified by clustering analysis based on single-cell sequencing. Among them, CD8+ T cells and NKT cells were universally exhausted by expressing exhaustion-associated molecular markers. Interestingly, signatures of CD8+ T cells and NKT cells predicted prolonged OS for PDAC only in the presence of "targets" for pyroptosis and ferroptosis induction. Moreover, a specific state of T cells with overexpression of ribosome-related proteins was associated with a good prognosis. In addition, the hematopoietic stem cell (HSC)-like signature predicted prolonged OS in PDAC. Weighted gene co-expression network analysis identified 5 hub genes whose downregulation may mediate the observed survival benefits of the HSC-like signature. Moreover, trajectory analysis revealed that myeloid cells evolutionarily consisted of 7 states, and antigen-presenting molecules and complement-associated genes were lost along the pseudotime flow. Consensus clustering based on the differentially expressed genes between two states harboring the longest pseudotime span identified two PDAC groups with prognostic differences, and more infiltrated immune cells and activated immune signatures may account for the survival benefits. Conclusion: This study systematically investigated the prognostic implications of the components of the PDAC tumor microenvironment by integrating single-cell sequencing and bulk sequencing, and future studies are expected to develop novel targeted agents for PDAC treatment.

Tang R ，Liu X ，Liang C ，Hua J ，Xu J ，Wang W ，Meng Q ，Liu J ，Zhang B ，Yu X ，Shi S ... -  《Frontiers in Immunology》

Single-cell transcriptome analysis of tumor and stromal compartments of pancreatic ductal adenocarcinoma primary tumors and metastatic lesions.

Lin W ，Noel P ，Borazanci EH ，Lee J ，Amini A ，Han IW ，Heo JS ，Jameson GS ，Fraser C ，Steinbach M ，Woo Y ，Fong Y ，Cridebring D ，Von Hoff DD ，Park JO ，Han H ... -  《Genome Medicine》

BNIP3+ fibroblasts associated with hypoxia and inflammation predict prognosis and immunotherapy response in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most malignant tumors that lacks effective treatment options. Cancer-associated fibroblasts (CAFs), an important component of the tumor microenvironment, associated with tumor progression, prognosis, and treatment response. This work aimed to explore the novel CAFs-associated target to improve treatment strategies in PDAC. The PDAC single-cell sequencing data (CRA001160, n = 35) were downloaded and integrated based on GSA databases to classify fibroblasts into fine subtypes. Functional enrichment analysis and coexpression regulatory network analysis were used to identify the functional phenotypes and biological properties of the different fibroblast subtypes. Fibroblast differentiation trajectories were constructed using pseudochronological analysis to identify initial and terminally differentiated subtypes of fibroblasts. The changes in the proportions of different fibroblast subtypes before and after PDAC immunotherapy were compared in responsive and nonresponding patients, and the relationships between fibroblast subtypes and PDAC immunotherapy responsiveness were determined based on GSA and GEO database. Using molecular biology methods to confirm the effects of BNIP3 on hypoxia and inflammation in CAFs. CAFs were co cultured with pancreatic cancer cells to detect their effects on migration and invasion of pancreatic cancer. Single-cell data analysis divided fibroblasts into six subtypes. The differentiation trajectory suggested that BNIP3+ Fibro subtype exhibited terminal differentiation, and the expression of genes related to hypoxia and the inflammatory response increased gradually with differentiation time. The specific overexpressed genes in the BNIP3+ Fibro subtype were significantly associated with overall and disease progression-free survival in the patients with PDAC. Interestingly, the greater the proportion of the BNIP3+ Fibro subtype was, the worse the response of PDAC patients to immunotherapy, and the CRTL treatment regimen effectively reduced the proportion of the BNIP3+ Fibro subtype. After knocking out BNIP3, the hypoxia markers and inflammatory factors of CAFs were inhibited. Co-culture of CAFs with pancreatic cancer cells can increase the migration and invasion of pancreatic cancer, but this could be reversed by knocking out BNIP3. This study revealed the BNIP3+ Fibro subtype associated with hypoxia and inflammatory responses, which was closely related to the poor prognosis of patients with PDAC, and identified signature genes that predict the immunotherapy response in PDAC.

Gao B ，Hu G ，Sun B ，Li W ，Yang H ... -  《Journal of Translational Medicine》

Inter- and intra-tumoural heterogeneity in cancer-associated fibroblasts of human pancreatic ductal adenocarcinoma.

Cancer-associated fibroblasts (CAF) are orchestrators of the pancreatic ductal adenocarcinoma (PDAC) microenvironment. Stromal heterogeneity may explain differential pathophysiological roles of the stroma (pro- versus anti-tumoural) in PDAC. We hypothesised that multiple CAF functional subtypes exist in PDAC, that contribute to stromal heterogeneity through interactions with cancer cells. Using molecular and functional analysis of patient-derived CAF primary cultures, we demonstrated that human PDAC-derived CAFs display a high level of inter- and intra-tumour heterogeneity. We identified at least four subtypes of CAFs based on transcriptomic analysis, and propose a classification for human PDAC-derived CAFs (pCAFassigner). Multiple CAF subtypes co-existed in individual patient samples. The presence of these CAF subtypes in bulk tumours was confirmed using publicly available gene expression profiles, and immunostainings of CAF subtype markers. Each subtype displayed specific phenotypic features (matrix- and immune-related signatures, vimentin and α-smooth muscle actin expression, proliferation rate), and was associated with an assessable prognostic impact. A prolonged exposure of non-tumoural pancreatic stellate cells to conditioned media from cancer cell lines (cancer education experiment) induced a CAF-like phenotype, including loss of capacity to revert to quiescence and an increase in the expression of genes related to CAF subtypes B and C. This classification demonstrates molecular and functional inter- and intra-tumoural heterogeneity of CAFs in human PDAC. Our subtypes overlap with those identified from single-cell analyses in other cancers, and pave the way for the development of therapies targeting specific CAF subpopulations in PDAC. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Neuzillet C ，Tijeras-Raballand A ，Ragulan C ，Cros J ，Patil Y ，Martinet M ，Erkan M ，Kleeff J ，Wilson J ，Apte M ，Tosolini M ，Wilson AS ，Delvecchio FR ，Bousquet C ，Paradis V ，Hammel P ，Sadanandam A ，Kocher HM ... -  《-》