Pancreatic neuroendocrine tumors (PNETs) are uncommon malignancies characterized by significant clinical heterogeneity and a pronounced tendency for liver metastasis. Despite this, the cellular mechanisms driving PNET progression remain inadequately elucidated, especially concerning neuroendocrine subpopulations. We analyzed publicly available single-cell RNA sequencing (scRNA-seq) data from 27 samples, including adjacent normal tissues (NT), primary tumors (PT), and hepatic metastases (HM), to explore the heterogeneity of neuroendocrine cells. Our downstream analyses encompassed copy number variation (CNV) inference, pseudotime trajectory modeling using CytoTRACE2 and Monocle3, Single-Cell Regulatory Network Inference and Clustering (SCENIC), cell–cell communication analysis via CellChat, and external validation with bulk RNA-seq datasets. We identified a distinct PCSK1+ neuroendocrine cell subpopulation, predominantly found in HM, which exhibited a high CNV burden, low differentiation potential, and significant transcriptional divergence from the NEUROD1+ neuroendocrine cells prevalent in primary tumors. The trajectory analysis delineated a developmental continuum commencing from the NEUROD1+ subpopulation and culminating in the PCSK1+ subpopulation. SCENIC analysis identified ATF6 as a pivotal transcriptional regulator within the PCSK1+ subpopulation, while KEGG enrichment of its target genes indicated involvement in stress-related signaling pathways. Furthermore, cell–cell communication analysis demonstrated that fibroblasts were the predominant signaling source to the PCSK1+ subpopulation in both primary tumors (PT) and hepatic metastases (HM), with conserved ligand–receptor interactions, including the CD99–CD99 and collagen–integrin axes. Our study identifies a metastasis-enriched, terminally differentiated PCSK1+ subpopulation and elucidates its potential regulatory and microenvironmental characteristics. These findings enhance our understanding of the cellular states linked to the progression of PNETs and lay the groundwork for subsequent mechanistic investigations.
{"title":"Single-cell transcriptomic analysis reveals a metastasis-associated PCSK1+ neuroendocrine subpopulation in pancreatic neuroendocrine tumors","authors":"Xiaolei Yin, Xiaopeng Li, Lili Mi, Jiaojiao Hou, Fei Yin","doi":"10.1111/jne.70084","DOIUrl":"10.1111/jne.70084","url":null,"abstract":"<p>Pancreatic neuroendocrine tumors (PNETs) are uncommon malignancies characterized by significant clinical heterogeneity and a pronounced tendency for liver metastasis. Despite this, the cellular mechanisms driving PNET progression remain inadequately elucidated, especially concerning neuroendocrine subpopulations. We analyzed publicly available single-cell RNA sequencing (scRNA-seq) data from 27 samples, including adjacent normal tissues (NT), primary tumors (PT), and hepatic metastases (HM), to explore the heterogeneity of neuroendocrine cells. Our downstream analyses encompassed copy number variation (CNV) inference, pseudotime trajectory modeling using CytoTRACE2 and Monocle3, Single-Cell Regulatory Network Inference and Clustering (SCENIC), cell–cell communication analysis via CellChat, and external validation with bulk RNA-seq datasets. We identified a distinct PCSK1<sup>+</sup> neuroendocrine cell subpopulation, predominantly found in HM, which exhibited a high CNV burden, low differentiation potential, and significant transcriptional divergence from the NEUROD1<sup>+</sup> neuroendocrine cells prevalent in primary tumors. The trajectory analysis delineated a developmental continuum commencing from the NEUROD1<sup>+</sup> subpopulation and culminating in the PCSK1<sup>+</sup> subpopulation. SCENIC analysis identified ATF6 as a pivotal transcriptional regulator within the PCSK1<sup>+</sup> subpopulation, while KEGG enrichment of its target genes indicated involvement in stress-related signaling pathways. Furthermore, cell–cell communication analysis demonstrated that fibroblasts were the predominant signaling source to the PCSK1<sup>+</sup> subpopulation in both primary tumors (PT) and hepatic metastases (HM), with conserved ligand–receptor interactions, including the CD99–CD99 and collagen–integrin axes. Our study identifies a metastasis-enriched, terminally differentiated PCSK1<sup>+</sup> subpopulation and elucidates its potential regulatory and microenvironmental characteristics. These findings enhance our understanding of the cellular states linked to the progression of PNETs and lay the groundwork for subsequent mechanistic investigations.</p>","PeriodicalId":16535,"journal":{"name":"Journal of Neuroendocrinology","volume":"37 11","pages":""},"PeriodicalIF":4.1,"publicationDate":"2025-08-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144859286","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ulrich Peter Knigge, Magnus Kjellman, Henning Grønbæk, Espen Thiis-Evensen, Camilla Schalin-Jäntti, Staffan Welin, Halfdan Sørbye, Maria del Pilar Schneider, Roger Belusa, The Nordic NET Biomarker Group
There is an unmet need for new methods to predict disease course in patients with neuroendocrine tumors (NET). We investigated 92 putative cancer-related plasma proteins including chromogranin A (CgA) and clinical parameters at the time of diagnosis to identify early factors associated with progressive (PD) or stable disease (SD). Patients with NET grade 1 and 2 of the small intestine (siNET) and pancreas (pNET) were included in this prospective study. Blood samples were obtained at the time of diagnosis before tumor-related therapy was initiated. During 3 years of follow-up, SD or PD was determined according to current clinical practice by each investigator. Association rule mining (ARM) was used to identify combinations of biomarkers and clinical parameters associated with SD or PD. Altogether, 115 patients with siNET and 30 with pNET with complete clinical and biomarker data were included in the analysis representing 3 years of follow-up. Several novel plasma proteins and clinical factors were associated with either PD or SD. In siNET, CgA (>4 upper limits of normal [ULN]) was the most frequent biomarker associated with PD. Females, in contrast to males, with CgA >4 ULN showed a high risk of progression (PPV 100%, NPV 63%). In the siNET cohort, Carboxypeptidase E (CPE) was a discriminating factor between SD and PD. CPE <3.03 was associated with SD, whereas CPE >3.14 was associated with PD (p = 0.003). In the pNET cohort, among clinical variables, only the presence of liver metastasis was associated with PD. CgA was not among the top biomarkers associated with PD. Several parameters, both clinical and biomarker data, as well as combinations of these, were associated with PD or SD 3 years after diagnosis. We identified novel biomarkers improving the association with PD or SD. [Correction added on 28 August 2025, after first online publication: Abstract has been updated.]
{"title":"Association rule mining of clinical and biomarker data in neuroendocrine tumors: A prospective study on disease progression","authors":"Ulrich Peter Knigge, Magnus Kjellman, Henning Grønbæk, Espen Thiis-Evensen, Camilla Schalin-Jäntti, Staffan Welin, Halfdan Sørbye, Maria del Pilar Schneider, Roger Belusa, The Nordic NET Biomarker Group","doi":"10.1111/jne.70069","DOIUrl":"10.1111/jne.70069","url":null,"abstract":"<p>There is an unmet need for new methods to predict disease course in patients with neuroendocrine tumors (NET). We investigated 92 putative cancer-related plasma proteins including chromogranin A (CgA) and clinical parameters at the time of diagnosis to identify early factors associated with progressive (PD) or stable disease (SD). Patients with NET grade 1 and 2 of the small intestine (siNET) and pancreas (pNET) were included in this prospective study. Blood samples were obtained at the time of diagnosis before tumor-related therapy was initiated. During 3 years of follow-up, SD or PD was determined according to current clinical practice by each investigator. Association rule mining (ARM) was used to identify combinations of biomarkers and clinical parameters associated with SD or PD. Altogether, 115 patients with siNET and 30 with pNET with complete clinical and biomarker data were included in the analysis representing 3 years of follow-up. Several novel plasma proteins and clinical factors were associated with either PD or SD. In siNET, CgA (>4 upper limits of normal [ULN]) was the most frequent biomarker associated with PD. Females, in contrast to males, with CgA >4 ULN showed a high risk of progression (PPV 100%, NPV 63%). In the siNET cohort, Carboxypeptidase E (CPE) was a discriminating factor between SD and PD. CPE <3.03 was associated with SD, whereas CPE >3.14 was associated with PD (<i>p</i> = 0.003). In the pNET cohort, among clinical variables, only the presence of liver metastasis was associated with PD. CgA was not among the top biomarkers associated with PD. Several parameters, both clinical and biomarker data, as well as combinations of these, were associated with PD or SD 3 years after diagnosis. We identified novel biomarkers improving the association with PD or SD. [Correction added on 28 August 2025, after first online publication: Abstract has been updated.]</p>","PeriodicalId":16535,"journal":{"name":"Journal of Neuroendocrinology","volume":"37 10","pages":""},"PeriodicalIF":4.1,"publicationDate":"2025-08-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jne.70069","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144855575","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Håkan Ohlsson, Martin Nilsson, Anna Sundlöv, Marlene Malmström, Martin Almquist
Health-related quality of life (HRQoL) has been shown to predict overall survival (OS) in several different malignancies, but not in patients with small intestinal neuroendocrine tumours (siNET). We evaluated the influence of HRQoL on survival in patients with siNET. We included 85 patients with advanced siNET who completed the validated HRQoL instruments, EORTC QLQ-C30 and GI.NET21. We used Cox proportional hazards to calculate the hazard ratio (HR) of survival according to the QLQ-C30 summary score, adjusting for clinical variables selected with causal inference. Flexible parametric modelling using cubic splines was used to illustrate the time-dependent relationship between HRQoL and OS. The QLQ-C30 summary score was correlated with overall survival (OS) with an adjusted HR of 0.62 (95% CI 0.46–0.83, p < .001) for each 10-point increase in summary score. Compared to a model using only clinical variables, the summary score increased predictive accuracy by eight percentage units and improved model fit. The inclusion of GI.NET21 with the summary score yielded similar results. HRQoL predicts overall survival in patients with siNET, providing additional information to clinical variables. Measuring HRQoL might be useful when following patients with siNET.
健康相关生活质量(HRQoL)已被证明可预测几种不同恶性肿瘤的总生存期(OS),但不适用于小肠神经内分泌肿瘤(siNET)患者。我们评估了HRQoL对siNET患者生存的影响。我们纳入了85例晚期siNET患者,他们完成了经验证的HRQoL仪器、EORTC QLQ-C30和GI.NET21。根据QLQ-C30综合评分,采用Cox比例风险法计算生存风险比(HR),并对经因果推理选择的临床变量进行调整。采用三次样条灵活的参数化建模来说明HRQoL和OS之间的时间依赖关系。QLQ-C30总评分与总生存期(OS)相关,调整后风险比为0.62 (95% CI 0.46-0.83, p
{"title":"Quality of life as a predictor for survival in patients with small intestinal neuroendocrine tumours","authors":"Håkan Ohlsson, Martin Nilsson, Anna Sundlöv, Marlene Malmström, Martin Almquist","doi":"10.1111/jne.70081","DOIUrl":"10.1111/jne.70081","url":null,"abstract":"<p>Health-related quality of life (HRQoL) has been shown to predict overall survival (OS) in several different malignancies, but not in patients with small intestinal neuroendocrine tumours (siNET). We evaluated the influence of HRQoL on survival in patients with siNET. We included 85 patients with advanced siNET who completed the validated HRQoL instruments, EORTC QLQ-C30 and GI.NET21. We used Cox proportional hazards to calculate the hazard ratio (HR) of survival according to the QLQ-C30 summary score, adjusting for clinical variables selected with causal inference. Flexible parametric modelling using cubic splines was used to illustrate the time-dependent relationship between HRQoL and OS. The QLQ-C30 summary score was correlated with overall survival (OS) with an adjusted HR of 0.62 (95% CI 0.46–0.83, <i>p</i> < .001) for each 10-point increase in summary score. Compared to a model using only clinical variables, the summary score increased predictive accuracy by eight percentage units and improved model fit. The inclusion of GI.NET21 with the summary score yielded similar results. HRQoL predicts overall survival in patients with siNET, providing additional information to clinical variables. Measuring HRQoL might be useful when following patients with siNET.</p>","PeriodicalId":16535,"journal":{"name":"Journal of Neuroendocrinology","volume":"37 11","pages":""},"PeriodicalIF":4.1,"publicationDate":"2025-08-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1111/jne.70081","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144835382","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}