ABSTRACT
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.
