Single-cell analysis of castration-resistant prostate cancers to identify potential biomarkers for diagnosis and prognosis of neuroendocrine prostate cancer

The high heterogeneity and low percentage of neuroendocrine cells in prostate cancer limit the utility of traditional bulk RNA sequencing and even single-cell RNA sequencing to find better biomarkers for early diagnosis and stratification. Re-clustering of specific cell-type holds great promise for identification of intra-cell-type heterogeneity. However, this has not yet been used in studying neuroendocrine prostate cancer heterogeneity. Neuroendocrine cluster(s) were individually identified in each castration-resistant prostate cancer specimen and combined for trajectory analysis. Three neuroendocrine states were identified. Neuroendocrine state 2 with the highest AR score was considered the initial starting state of neuroendocrine transdifferentiation. State 1 and state 3 with distinct high neuroendocrine scores and marker genes enriched in N-Myc and REST target genes, respectively, were considered as two different types of neuroendocrine differentiated cancer cells. These two states contained distinct groups of prostate cancer biomarkers and a strong distinguishing ability of normal versus cancerous prostate across different pathological grading was found in the N-Myc-associated state. Our data highlight the central role of N-Myc and REST in mediating lineage plasticity and classifying neuroendocrine phenotypes.