Evolutionarily conservative and non-conservative regulatory networks during primate interneuron development revealed by single-cell RNA and ATAC sequencing

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Single-cell transcriptomic reveals a cell atlas and diversity of chicken amygdala responded to social hierarchy

Amygdala serves as a highly cellular, heterogeneous brain region containing excitatory and inhibitory neurons and is involved in the dopamine and serotoninergic neuron systems. An increasing number of studies have revealed the underpinned mechanism mediating social hierarchy in mammal and vertebrate, however, there are rare studies conducted on how amygdala on social hierarchy in poultry. In this study, we conducted food competition tests and determined the social hierarchy of the rooster. We performed cross-species analysis with mammalian amygdala, and found that cell types of human and rhesus monkeys were more closely related and that of chickens were more distant. We identified 26 clusters and divided them into 10 main clusters, of which GABAergic and glutamatergic neurons were associated with social behaviors. In conclusion, our results provide to serve the developmental studies of the amygdala neuron system and new insights into the underpinned mechanism of social hierarchy in roosters.

Parvalbumin Interneuron Dysfunction in Neurological Disorders: Focus on Epilepsy and Alzheimer's Disease

Parvalbumin expressing (PV+) GABAergic interneurons are fast spiking neurons that provide powerful but relatively short-lived inhibition to principal excitatory cells in the brain. They play a vital role in feedforward and feedback synaptic inhibition, preventing run away excitation in neural networks. Hence, their dysfunction can lead to hyperexcitability and increased susceptibility to seizures. PV+ interneurons are also key players in generating gamma oscillations, which are synchronized neural oscillations associated with various cognitive functions. PV+ interneuron are particularly vulnerable to aging and their degeneration has been associated with cognitive decline and memory impairment in dementia and Alzheimer's disease (AD). Overall, dysfunction of PV+ interneurons disrupts the normal excitatory/inhibitory balance within specific neurocircuits in the brain and thus has been linked to a wide range of neurodevelopmental and neuropsychiatric disorders. This review focuses on the role of dysfunctional PV+ inhibitory interneurons in the generation of epileptic seizures and cognitive impairment and their potential as targets in the design of future therapeutic strategies to treat these disorders. Recent research using cutting-edge optogenetic and chemogenetic technologies has demonstrated that they can be selectively manipulated to control seizures and restore the balance of neural activity in the brains of animal models. This suggests that PV+ interneurons could be important targets in developing future treatments for patients with epilepsy and comorbid disorders, such as AD, where seizures and cognitive decline are directly linked to specific PV+ interneuron deficits.

PCDHA9 as a candidate gene for amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease. To identify additional genetic factors, we analyzed exome sequences in a large cohort of Chinese ALS patients and found a homozygous variant (p.L700P) in PCDHA9 in three unrelated patients. We generated Pcdhα9 mutant mice harboring either orthologous point mutation or deletion mutation. These mice develop progressive spinal motor loss, muscle atrophy, and structural/functional abnormalities of the neuromuscular junction, leading to paralysis and early lethality. TDP-43 pathology is detected in the spinal motor neurons of aged mutant mice. Mechanistically, we demonstrate that Pcdha9 mutation causes aberrant activation of FAK and PYK2 in aging spinal cord, and dramatically reduced NKA-α1 expression in motor neurons. Our single nucleus multi-omics analysis reveals disturbed signaling involved in cell adhesion, ion transport, synapse organization, and neuronal survival in aged mutant mice. Together, our results present PCDHA9 as a potential ALS gene and provide insights into its pathogenesis.

Investigating the impact of inflammatory response-related genes on renal fibrosis diagnosis: a machine learning-based study with experimental validation

Renal fibrosis plays a crucial role in the progression of renal diseases, yet the lack of effective diagnostic markers poses challenges in scientific and clinical practices. In this study, we employed machine learning techniques to identify potential biomarkers for renal fibrosis. Utilizing two datasets from the GEO database, we applied LASSO, SVM-RFE and RF algorithms to screen for differentially expressed genes related to inflammatory responses between the renal fibrosis group and the control group. As a result, we identified four genes (CCL5, IFITM1, RIPK2, and TNFAIP6) as promising diagnostic indicators for renal fibrosis. These genes were further validated through in vivo experiments and immunohistochemistry, demonstrating their utility as reliable markers for assessing renal fibrosis. Additionally, we conducted a comprehensive analysis to explore the relationship between these candidate biomarkers, immunity, and drug sensitivity. Integrating these findings, we developed a nomogram with a high discriminative ability, achieving a concordance index of 0.933, enabling the prediction of disease risk in patients with renal fibrosis. Overall, our study presents a predictive model for renal fibrosis and highlights the significance of four potential biomarkers, facilitating clinical diagnosis and personalized treatment. This finding presents valuable insights for advancing precision medicine approaches in the management of renal fibrosis.Communicated by Ramaswamy H. Sarma.

Cortical somatostatin long-range projection neurons and interneurons exhibit divergent developmental trajectories

The mammalian cerebral cortex contains an extraordinary diversity of cell types that emerge by implementing different developmental programs. Delineating when and how cellular diversification occurs is particularly challenging for cortical inhibitory neurons because they represent a small proportion of all cortical cells and have a protracted development. Here, we combine single-cell RNA sequencing and spatial transcriptomics to characterize the emergence of neuronal diversity among somatostatin-expressing (SST+) cells in mice. We found that SST+ inhibitory neurons segregate during embryonic stages into long-range projection (LRP) neurons and two types of interneurons, Martinotti cells and non-Martinotti cells, following distinct developmental trajectories. Two main subtypes of LRP neurons and several subtypes of interneurons are readily distinguishable in the embryo, although interneuron diversity is likely refined during early postnatal life. Our results suggest that the timing for cellular diversification is unique for different subtypes of SST+ neurons and particularly divergent for LRP neurons and interneurons.

Spatial and single-nucleus transcriptomics decoding the molecular landscape and cellular organization of avian optic tectum

The avian optic tectum (OT) has been studied for its diverse functions, yet a comprehensive molecular landscape at the cellular level has been lacking. In this study, we applied spatial transcriptome sequencing and single-nucleus RNA sequencing (snRNA-seq) to explore the cellular organization and molecular characteristics of the avian OT from two species: Columba livia and Taeniopygia guttata. We identified precise layer structures and provided comprehensive layer-specific signatures of avian OT. Furthermore, we elucidated diverse functions in different layers, with the stratum griseum periventriculare (SGP) potentially playing a key role in advanced functions of OT, like fear response and associative learning. We characterized detailed neuronal subtypes and identified a population of FOXG1+ excitatory neurons, resembling those found in the mouse neocortex, potentially involved in neocortex-related functions and expansion of avian OT. These findings could contribute to our understanding of the architecture of OT, shedding light on visual perception and multifunctional association.

Loss of Ezh2 in the medial ganglionic eminence alters interneuron fate, cell morphology and gene expression profiles

Introduction

Enhancer of zeste homolog 2 (Ezh2) is responsible for trimethylation of histone 3 at lysine 27 (H3K27me3), resulting in repression of gene expression. Here, we explore the role of Ezh2 in forebrain GABAergic interneuron development.

Methods

We removed Ezh2 in the MGE by generating Nkx2-1Cre;Ezh2 conditional knockout mice. We then characterized changes in MGE-derived interneuron fate and electrophysiological properties in juvenile mice, as well as alterations in gene expression, chromatin accessibility and histone modifications in the MGE.

Results

Loss of Ezh2 increases somatostatin-expressing (SST+) and decreases parvalbumin-expressing (PV+) interneurons in the forebrain. We observe fewer MGE-derived interneurons in the first postnatal week, indicating reduced interneuron production. Intrinsic electrophysiological properties in SST+ and PV+ interneurons are normal, but PV+ interneurons display increased axonal complexity in Ezh2 mutant mice. Single nuclei multiome analysis revealed differential gene expression patterns in the embryonic MGE that are predictive of these cell fate changes. Lastly, CUT&Tag analysis revealed that some genomic loci are particularly resistant or susceptible to shifts in H3K27me3 levels in the absence of Ezh2, indicating differential selectivity to epigenetic perturbation.

Discussion

Thus, loss of Ezh2 in the MGE alters interneuron fate, morphology, and gene expression and regulation. These findings have important implications for both normal development and potentially in disease etiologies.

Patterning ganglionic eminences in developing human brain organoids using a morphogen-gradient-inducing device

In early neurodevelopment, the central nervous system is established through the coordination of various neural organizers directing tissue patterning and cell differentiation. Better recapitulation of morphogen gradient production and signaling will be crucial for establishing improved developmental models of the brain in vitro. Here, we developed a method by assembling polydimethylsiloxane devices capable of generating a sustained chemical gradient to produce patterned brain organoids, which we termed morphogen-gradient-induced brain organoids (MIBOs). At 3.5 weeks, MIBOs replicated dorsal-ventral patterning observed in the ganglionic eminences (GE). Analysis of mature MIBOs through single-cell RNA sequencing revealed distinct dorsal GE-derived CALB2+ interneurons, medium spiny neurons, and medial GE-derived cell types. Finally, we demonstrate long-term culturing capabilities with MIBOs maintaining stable neural activity in cultures grown up to 5.5 months. MIBOs demonstrate a versatile approach for generating spatially patterned brain organoids for embryonic development and disease modeling.

Apically localized PANX1 impacts neuroepithelial expansion in human cerebral organoids

Dysfunctional paracrine signaling through Pannexin 1 (PANX1) channels is linked to several adult neurological pathologies and emerging evidence suggests that PANX1 plays an important role in human brain development. It remains unclear how early PANX1 influences brain development, or how loss of PANX1 alters the developing human brain. Using a cerebral organoid model of early human brain development, we find that PANX1 is expressed at all stages of organoid development from neural induction through to neuroepithelial expansion and maturation. Interestingly, PANX1 cellular distribution and subcellular localization changes dramatically throughout cerebral organoid development. During neural induction, PANX1 becomes concentrated at the apical membrane domain of neural rosettes where it co-localizes with several apical membrane adhesion molecules. During neuroepithelial expansion, PANX1-/- organoids are significantly smaller than control and exhibit significant gene expression changes related to cell adhesion, WNT signaling and non-coding RNAs. As cerebral organoids mature, PANX1 expression is significantly upregulated and is primarily localized to neuronal populations outside of the ventricular-like zones. Ultimately, PANX1 protein can be detected in all layers of a 21-22 post conception week human fetal cerebral cortex. Together, these results show that PANX1 is dynamically expressed by numerous cell types throughout embryonic and early fetal stages of human corticogenesis and loss of PANX1 compromises neuroepithelial expansion due to dysregulation of cell-cell and cell-matrix adhesion, perturbed intracellular signaling, and changes to gene regulation.

Multi-omics approach reveals dysregulated genes during hESCs neuronal differentiation exposure to paracetamol

Prenatal paracetamol exposure has been associated with neurodevelopmental outcomes in childhood. Pharmacoepigenetic studies show differences in cord blood DNA methylation between unexposed and paracetamol-exposed neonates, however, causality and impact of long-term prenatal paracetamol exposure on brain development remain unclear. Using a multi-omics approach, we investigated the effects of paracetamol on an in vitro model of early human neurodevelopment. We exposed human embryonic stem cells undergoing neuronal differentiation with paracetamol concentrations corresponding to maternal therapeutic doses. Single-cell RNA-seq and ATAC-seq integration identified paracetamol-induced chromatin opening changes linked to gene expression. Differentially methylated and/or expressed genes were involved in neurotransmission and cell fate determination trajectories. Some genes involved in neuronal injury and development-specific pathways, such as KCNE3, overlapped with differentially methylated genes previously identified in cord blood associated with prenatal paracetamol exposure. Our data suggest that paracetamol may play a causal role in impaired neurodevelopment.

Innovative breakthroughs facilitated by single-cell multi-omics: manipulating natural killer cell functionality correlates with a novel subcategory of melanoma cells

Background

Melanoma is typically regarded as the most dangerous form of skin cancer. Although surgical removal of in situ lesions can be used to effectively treat metastatic disease, this condition is still difficult to cure. Melanoma cells are removed in great part due to the action of natural killer (NK) and T cells on the immune system. Still, not much is known about how the activity of NK cell-related pathways changes in melanoma tissue. Thus, we performed a single-cell multi-omics analysis on human melanoma cells in this study to explore the modulation of NK cell activity.

Materials and methods

Cells in which mitochondrial genes comprised > 20% of the total number of expressed genes were removed. Gene ontology (GO), gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and AUCcell analysis of differentially expressed genes (DEGs) in melanoma subtypes were performed. The CellChat package was used to predict cell-cell contact between NK cell and melanoma cell subtypes. Monocle program analyzed the pseudotime trajectories of melanoma cells. In addition, CytoTRACE was used to determine the recommended time order of melanoma cells. InferCNV was utilized to calculate the CNV level of melanoma cell subtypes. Python package pySCENIC was used to assess the enrichment of transcription factors and the activity of regulons in melanoma cell subtypes. Furthermore, the cell function experiment was used to confirm the function of TBX21 in both A375 and WM-115 melanoma cell lines.

Results

Following batch effect correction, 26,161 cells were separated into 28 clusters and designated as melanoma cells, neural cells, fibroblasts, endothelial cells, NK cells, CD4+ T cells, CD8+ T cells, B cells, plasma cells, monocytes and macrophages, and dendritic cells. A total of 10137 melanoma cells were further grouped into seven subtypes, i.e., C0 Melanoma BIRC7, C1 Melanoma CDH19, C2 Melanoma EDNRB, C3 Melanoma BIRC5, C4 Melanoma CORO1A, C5 Melanoma MAGEA4, and C6 Melanoma GJB2. The results of AUCell, GSEA, and GSVA suggested that C4 Melanoma CORO1A may be more sensitive to NK and T cells through positive regulation of NK and T cell-mediated immunity, while other subtypes of melanoma may be more resistant to NK cells. This suggests that the intratumor heterogeneity (ITH) of melanoma-induced activity and the difference in NK cell-mediated cytotoxicity may have caused NK cell defects. Transcription factor enrichment analysis indicated that TBX21 was the most important TF in C4 Melanoma CORO1A and was also associated with M1 modules. In vitro experiments further showed that TBX21 knockdown dramatically decreases melanoma cells' proliferation, invasion, and migration.

Conclusion

The differences in NK and T cell-mediated immunity and cytotoxicity between C4 Melanoma CORO1A and other melanoma cell subtypes may offer a new perspective on the ITH of melanoma-induced metastatic activity. In addition, the protective factors of skin melanoma, STAT1, IRF1, and FLI1, may modulate melanoma cell responses to NK or T cells.

Construction of a telomere-related gene signature to predict prognosis and immune landscape for glioma

Background

Glioma is one of the commonest malignant tumors of the brain. However, glioma present with a poor clinical prognosis. Therefore, specific detection markers and therapeutic targets need to be explored as a way to promote the survival rate of BC patients. Therefore, we need to search for quality immune checkpoints to support the efficacy of immunotherapy for glioma.

Methods

We first recognized differentially expressed telomere-related genes (TRGs) and accordingly developed a risk model by univariate and multivariate Cox analysis. The accuracy of the model is then verified. We evaluated the variations in immune function and looked at the expression levels of immune checkpoint genes. Finally, to assess the anti-tumor medications often used in the clinical treatment of glioma, we computed the half inhibitory concentration of pharmaceuticals.

Results

We finally identified nine TRGs and built a risk model. Through the validation of the model, we found good agreement between the predicted and observed values. Then, we found 633 differentially expressed genes between various risk groups to identify the various molecular pathways between different groups. The enrichment of CD4+ T cells, CD8+ T cells, fibroblasts, endothelial cells, macrophages M0, M1, and M2, mast cells, myeloid dendritic cells, and neutrophils was favorably correlated with the risk score, but the enrichment of B cells and NK cells was negatively correlated with the risk score. The expression of several immune checkpoint-related genes differed significantly across the risk groups. Finally, in order to create individualized treatment plans for diverse individuals, we searched for numerous chemotherapeutic medications for patients in various groups.

Conclusion

The findings of this research provide evidence that TRGs may predict a patient's prognosis for glioma, assist in identifying efficient targets for glioma immunotherapy, and provide a foundation for an efficient, customized approach to treating glioma patients.

A novel signature predicts prognosis and immunotherapy in lung adenocarcinoma based on cancer-associated fibroblasts

Background

Extensive research has established the significant correlations between cancer-associated fibroblasts (CAFs) and various stages of cancer development, including initiation, angiogenesis, progression, and resistance to therapy. In this study, we aimed to investigate the characteristics of CAFs in lung adenocarcinoma (LUAD) and develop a risk signature to predict the prognosis of patients with LUAD.

Methods

We obtained single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the public database. The Seurat R package was used to process the scRNA-seq data and identify CAF clusters based on several biomarkers. CAF-related prognostic genes were further identified using univariate Cox regression analysis. To reduce the number of genes, Lasso regression was performed, and a risk signature was established. A novel nomogram that incorporated the risk signature and clinicopathological features was developed to predict the clinical applicability of the model. Additionally, we conducted immune landscape and immunotherapy responsiveness analyses. Finally, we performed in vitro experiments to verify the functions of EXO1 in LUAD.

Results

We identified 5 CAF clusters in LUAD using scRNA-seq data, of which 3 clusters were significantly associated with prognosis in LUAD. A total of 492 genes were found to be significantly linked to CAF clusters from 1731 DEGs and were used to construct a risk signature. Moreover, our immune landscape exploration revealed that the risk signature was significantly related to immune scores, and its ability to predict responsiveness to immunotherapy was confirmed. Furthermore, a novel nomogram incorporating the risk signature and clinicopathological features showed excellent clinical applicability. Finally, we verified the functions of EXP1 in LUAD through in vitro experiments.

Conclusions

The risk signature has proven to be an excellent predictor of LUAD prognosis, stratifying patients more appropriately and precisely predicting immunotherapy responsiveness. The comprehensive characterization of LUAD based on the CAF signature can predict the response of LUAD to immunotherapy, thus offering fresh perspectives into the management of LUAD patients. Our study ultimately confirms the role of EXP1 in facilitating the invasion and growth of tumor cells in LUAD. Nevertheless, further validation can be achieved by conducting in vivo experiments.

Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development

Inhibitory GABAergic interneurons originate in the embryonic medial ganglionic eminence (MGE) and control network activity in the neocortex. Dysfunction of these cells is believed to lead to runaway excitation underlying seizure-based neurological disorders such as epilepsy, autism, and schizophrenia. Despite their importance in heath and disease, our knowledge about the development of this diverse neuronal population remains incomplete. Here we conducted single-cell RNA sequencing (scRNA-seq) of human foetal MGE from 10 to 15 weeks post conception. These MGE tissues are composed of largely cycling progenitors and immature post-mitotic interneurons with characteristic regional marker expression. Analysis of integrated human and mouse MGE data revealed species-conserved transcriptomic profiles and regulatory programs. Moreover, we identified novel candidate transcription regulators for human interneuron differentiation. These findings provide a framework for in vitro modelling of interneuron development and a strategy for potentially enhancing interneuron production from human pluripotent stem cells.

Deciphering inhibitory neuron development: The paths to diversity

The regulation of fate decisions in progenitor cells lays the foundation for the generation of neuronal diversity and the formation of specialized circuits with remarkable processing capacity. Since the discovery more than 20 years ago that inhibitory (GABAergic) neurons originate from progenitors in the ventral part of the embryonic brain, numerous details about their development and function have been unveiled. GABAergic neurons are an extremely heterogeneous group, comprising many specialized subtypes of local interneurons and long-range projection neurons. Clearly distinguishable types emerge during postmitotic maturation, at a time when precursors migrate, morphologically mature, and establish synaptic connections. Yet, differentiation begins at an earlier stage within their progenitor domains, where a combination of birthdate and place of origin are key drivers. This review explains how new insights from single-cell sequencing inform our current understanding of how GABAergic neuron diversity emerges.

Cortical interneuron specification and diversification in the era of big data

Inhibition in the mammalian cerebral cortex is mediated by a small population of highly diverse GABAergic interneurons. These largely local neurons are interspersed among excitatory projection neurons and exert pivotal regulation on the formation and function of cortical circuits. We are beginning to understand the extent of GABAergic neuron diversity and how this is generated and shaped during brain development in mice and humans. In this review, we summarise recent findings and discuss how new technologies are being used to further advance our knowledge. Understanding how inhibitory neurons are generated in the embryo is an essential pre-requisite of stem cell therapy, an evolving area of research, aimed at correcting human disorders that result in inhibitory dysfunction.

The uses of 3D human brain organoids for neurotoxicity evaluations: A review

Neurotoxicity studies aim at understanding the toxic effects and mechanisms of toxicants to human central nervous systems (CNS). However, human brains are the most complex organs, whereas the most commonly used models, such as 2D cell cultures and animal brains, are probably too simple to predict the responses of human brains. Embryonic stem cells (ESCs) or induced pluripotent stem cells (iPSCs)-based 3D human brain organoids hold unprecedented promise for the understanding of neurodevelopment and brain disease development. This review summarizes recent advances of using 3D human brain organoids for neurotoxicity studies. Comparative studies showed that 3D human brain organoids could support the findings obtained by animal or cohort studies, indicating that 3D human brain organoids are reliable models to evaluate the developmental neurotoxicity. 3D human brain organoids have been used to understand the toxicological mechanisms by using both conventional toxicological methods to investigate the signaling pathway changes as well as single cell RNA-sequencing to understand the neuron diversity. Some studies also used brain organoids carrying gene mutations or with virus infections to understand the toxicological responses of brains under diseased conditions. Although there are still limitations associated, 3D human brain organoids are promising tools for future neurotoxicity studies.

Characterization of neuroendocrine regulation- and metabolism-associated molecular features and prognostic indicators with aid to clinical chemotherapy and immunotherapy of patients with pancreatic cancer

The worldwide prevalence of pancreatic cancer has been rising in recent decades, and its prognosis has not improved much. The imbalance of substance and energy metabolism in tumour cells is among the primary causes of tumour formation and occurrence, which is often controlled by the neuroendocrine system. We applied Cox and LASSO regression analysis to develop a neuroendocrine regulation- and metabolism-related prognostic risk score model with three genes (GSK3B, IL18 and VEGFA) for pancreatic cancer. TCGA dataset served as the training and internal validation sets, and GSE28735, GSE62452 and GSE57495 were designated as external validation sets. Patients classified as the low-risk population (category, group) exhibited considerably improved survival duration in contrast with those classified as the high-risk population, as determined by the Kaplan-Meier curve. Then, we combined all the samples, and divided them into three clusters using unsupervised clustering analysis. Unsupervised clustering, t-distributed stochastic neighbor embedding (t-SNE), and principal component analysis (PCA) were further utilized to demonstrate the reliability of the prognostic model. Moreover, the risk score was shown to independently function as a predictor of pancreatic cancer in both univariate and multivariate Cox regression analyses. The results of gene set enrichment analysis (GSEA) illustrated that the low-risk population was predominantly enriched in immune-associated pathways. "ESTIMATE" algorithm, single-sample GSEA (ssGSEA) and the Tumor Immune Estimation Resource (TIMER) database showed immune infiltration ratings were enhanced in the low-risk category in contrast with the high-risk group. Tumour immune dysfunction and exclusion (TIDE) database predicted that immunotherapy for pancreatic cancer may be more successful in the high-risk than in the low-risk population. Mutation analysis illustrated a positive link between the tumour mutation burden and risk score. Drug sensitivity analysis identified 44 sensitive drugs in the high- and low-risk population. GSK3B expression was negatively correlated with Oxaliplatin, and IL18 expression was negatively correlated with Paclitaxel. Lastly, we analyzed and verified gene expression at RNA and protein levels based on GENPIA platform, HPA database and quantitative real-time PCR. In short, we developed a neuroendocrine regulation- and metabolism-associated prognostic model for pancreatic cancer that takes into account the immunological microenvironment and drug sensitivity.