Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Dec 11;41(12):2066-2082.e9.
doi: 10.1016/j.ccell.2023.10.009. Epub 2023 Nov 22.

Temporal evolution reveals bifurcated lineages in aggressive neuroendocrine small cell prostate cancer trans-differentiation

Affiliations

Temporal evolution reveals bifurcated lineages in aggressive neuroendocrine small cell prostate cancer trans-differentiation

Chia-Chun Chen et al. Cancer Cell. .

Abstract

Trans-differentiation from an adenocarcinoma to a small cell neuroendocrine state is associated with therapy resistance in multiple cancer types. To gain insight into the underlying molecular events of the trans-differentiation, we perform a multi-omics time course analysis of a pan-small cell neuroendocrine cancer model (termed PARCB), a forward genetic transformation using human prostate basal cells and identify a shared developmental, arc-like, and entropy-high trajectory among all transformation model replicates. Further mapping with single cell resolution reveals two distinct lineages defined by mutually exclusive expression of ASCL1 or ASCL2. Temporal regulation by groups of transcription factors across developmental stages reveals that cellular reprogramming precedes the induction of neuronal programs. TFAP4 and ASCL1/2 feedback are identified as potential regulators of ASCL1 and ASCL2 expression. Our study provides temporal transcriptional patterns and uncovers pan-tissue parallels between prostate and lung cancers, as well as connections to normal neuroendocrine cell states.

Keywords: ASCL1; ASCL2; POU2F3, TFAP4; cancer; lineage plasticity; neuroendocrine; prostate; small cell; stem-like; trans-differentiation.

PubMed Disclaimer

Conflict of interest statement

Declaration of interests E.C. served as a paid consultant to DotQuant and received Institutional sponsored research funding unrelated to this work from AbbVie, Gilead, Sanofi, Zenith Epigenetics, Bayer Pharmaceuticals, Forma Therapeutics, Genentech, GSK, Janssen Research, Kronos Bio, Foghorn Therapeutics, and MacroGenics. P.S.N. has served as a paid consultant for Janssen, Merck, Bristol Myers Squibb and received research funding from Janssen for work unrelated to the present study. O.N.W. currently has consulting, equity, and/or board relationships with Trethera Corporation, Kronos Biosciences, Sofie Biosciences, Breakthrough Properties, Vida Ventures, Nammi Therapeutics, Two River, Iconovir, Appia BioSciences, Neogene Therapeutics, 76Bio, and Allogene Therapeutics. None of these companies contributed to or directed any of the research reported in this article. T.G.G. reports receiving an honorarium from Amgen, having consulting and equity agreements with Auron Therapeutics, Boundless Bio, Coherus BioSciences and Trethera Corporation. The lab of T.G.G. has completed a research agreement with ImmunoActiva. A provisional patent application related to this study was submitted.

Figures

Figure 1.
Figure 1.. Temporal gene expression programs of the PARCB transformation model reveal SCNPC trans-differentiation pathways.
(A) Schematic summary of PARCB time course study and representative Hematoxylin and eosin (H&E) staining and immunohistochemistry (IHC) staining of neuroendocrine markers (SYP and NCAM1) on sequential tumors from the tissue microarray. Time point (TP1-6) samples were sequenced using bulk RNA sequencing (green circle), bulk ATAC-sequencing (red circle) and/or single cell RNA sequencing (blue circle, tumors only). (B) Projection of the PARCB time course samples onto the PCA framework defined by pan-cancer clinical tumor datasets ,,–. LUAD: Lung adenocarcinoma. LUAD norm: lung adenocarcinoma adjacent normal tissue. SCLC: small cell lung cancer. PRAD: prostate adenocarcinoma. PRAD norm: prostate adenocarcinoma adjacent normal tissue. CRPC: castration resistant prostate cancer. SCNPC: small cell neuroendocrine prostate cancer. (C) Arage gene expression of selected SCNPC-associated proteins and markers. (D) Heatmap of hierarchical clusters (HC) of samples (columns) and corresponding differentially upregulated gene modules (rows). Differential expression defined by one HC vs all other HCs). (E) PCA of the PARCB time course samples and trans-differentiation trajectories including primary arc and secondary bifurcation. A 3-dimensional rotatable version of this figure is available on the PARCB Multi-omics Explorer website. [For review, a 3D rotatable version is included as Data S1.] (F) Selected enriched GO terms across HC. See also Figure S1.
Figure 2.
Figure 2.. Sequential transcription regulators modulate reprogramming and neuroendocrine programs through a highly entropic and accessible chromatin state.
(A) Overall differential chromatin accessibility across HC. (B) PCA of chromatin accessibility of PARCB time course samples with entropy analysis using ATAC sequencing. (C) Overall mean accessible peaks near TSS of each HC in PARCB time course study. (D) Enriched motifs from suites of transcription factors in each HC using ATAC-sequencing. Top 5 motif suites for each comparison are shown, with additional analysis in Figure S2B, and full results in Table S1D. (E) Top ranked transcription factors and known neuroendocrine transcription factors across PARCB time course using bulk RNA sequencing. HOXC TFs avg: Average expression of HOXC4, HOXC5, HOXC6, HOXC8, HOXC9, HOXC10, HOXC11, HOXC12 and HOXC13. (F) Expression of ASCL1, ASCL2, NEURDO1 and POU2F3 in each HC. See also Figure S2.
Figure 3.
Figure 3.. Transcription factor-defined cell populations contribute to lineage divergence and tumor heterogeneity.
(A) Dimension reduction UMAP analysis of four patient series (P2, P5, P6 and P7) over time (TP3-6) using single cell RNA sequencing. (B) Temporal UMAP analysis of all the samples. (C) Expression of selected markers and transcription factors. KRT5 marks basal cells. KRT15 marks luminal cells. The expression is presented in log normalized counts. (D) Top enriched inferred cell types from the Human Cell Type Database using SingleR . (E) Projection of single cell RNA-sequencing samples on PCA framework by bulk RNA-sequencing samples (top panel) and the expression of selected markers and transcription factors (bottom panel). Each data point is a single cell colored by their corresponding HC. (F) Expression of ASCL1 (top) and ASCL2 (middle) and percentage of ASCL1/2 positive cells (cells with expression value >0) (bottom) in human biopsy and GEMM model tumors from five single cell RNA-sequencing datasets ,–. Other: prostatic intraepithelial neoplasia. NMYC_RB_M: Ptenf/f; Rb1f/f;MYCN + (PRN) and RB_M: Ptenf/f; Rb1f/f (PR) mouse model in Brady et al. See also Figure S3.
Figure 4.
Figure 4.. ASCL1 and ASCL2 specify independent transcriptional programs and sub-lineages in SCNPC.
(A) Inferred clonal tracing analysis of the PARCB time course samples using Monocle 2 . (B) Relative expression of KRT5, ASCL1 and ASCL2 in the inferred clonal tracing analysis (pseudo-time). (C) Percentages of ASCL1 or ASCL2 positive, double positive and double negative cell populations over time. (D) Volcano plot of differential gene expression in high ASCL1+ vs high ASCL2+ cell populations. (E) Representative genes from the predicted transcriptional programs of ASCL1 and ASCL2 trained on data from patient and model prostate cancer tumors (,, including TCGA), as determined by the ARACNE algorithm . (F) Western blot of panel of genes in the PARCB tumor derived cell lines from different tissue of origin (prostate, bladder and lung) ,. (G) Representative images of in situ hybridization of ASCL1 and ASCL2 mRNA analysis on transitional tumors (P7-TP5 and P9-TP4). See also Figure S4.
Figure 5.
Figure 5.. ASCL1 and ASCL2 as pan-cancer classifiers.
See also Figure S5. (A) Projection of the PARCB time course samples on the PCA framework defined by the CRPC subtypes using RNA sequencing (left) and ATAC-sequencing (right) . SCL: stem-cell like. NEPC: Neuroendocrine prostate cancer. 3-dimensional rotatable versions of these figures are available on the PARCB Multi-omics Explorer website. [For review, 3D rotatable versions are included as Data S2 and Data S3.] (B) Projection of the PARCB time course samples on the PCA framework defined by the SCLC subtypes ,. (C) mRNA expression of ASCL1 and ASCL2 in the PARCB time course samples and multiple sets of clinical CRPC-PRAD and SCNPC samples including TCGA and different research groups ,–. (D) Representative images of in situ RNA hybridization of ASCL1 and ASCL2 in clinical SCNPC tissues. (E) mRNA expression of ASCL1 and ASCL2 in pan cancer cell lines (CCLE). (F) mRNA expression of ASCL1 and ASCL2 in pan cancer tumors from TCGA.
Figure 6.
Figure 6.. Alternating ASCL1 and ASCL2 expression through reciprocal interaction and TFAP4 epigenetic regulation.
See also Figure S6. (A) Western blot analysis of exogenously expressing either V5 tagged ASCL2 in ASCL1+ cell lines (left) or V5-tagged ASCL1 in ASCL2+ cell lines (right). (B) Schematic of putative cis regulatory elements (CREs) of ASCL1 and ASCL2 (top) and the heatmap of open chromatin accessibility across CREs of ASCL1 and ASCL2 using the PARCB time course ATAC-sequencing (bottom). Red box: CREs containing predicted TFAP4 binding sites by HOMER motif enrichment analysis . (C) Top 8 ranked transcription factor motifs in ASCL1 promoter and ASCL2 enhancer regions, ranked by p-values. (D) Western blot analysis of doxycycline-inducible knockout of TFAP4 and proteins of interest in P7-TP6 (ASCL1+) and P3-TP5 (ASCL2+) cell lines. DOX: doxycycline. (E) Cell proliferation analysis of P7-TP6 (ASCL1+) and P3-TP5 (ASCL2+) cell lines with doxycycline-inducible knockout of TFAP4. Ctrl: no addition of doxycycline. TFAP4: with addition of doxycycline. (G) Schematic summary of the PARCB time course study.

References

    1. Nadal R, Schweizer M, Kryvenko ON, Epstein JI, and Eisenberger MA (2014). Small cell carcinoma of the prostate. Nat Rev Urol 11, 213–219. 10.1038/nrurol.2014.21. - DOI - PMC - PubMed
    1. Pietanza MC, Byers LA, Minna JD, and Rudin CM (2015). Small cell lung cancer: will recent progress lead to improved outcomes? Clin Cancer Res 21, 2244–2255. 10.1158/1078-0432.CCR-14-2958. - DOI - PMC - PubMed
    1. Yamada Y, and Beltran H (2021). Clinical and Biological Features of Neuroendocrine Prostate Cancer. Curr Oncol Rep 23, 15. 10.1007/s11912-020-01003-9. - DOI - PMC - PubMed
    1. Balanis NG, Sheu KM, Esedebe FN, Patel SJ, Smith BA, Park JW, Alhani S, Gomperts BN, Huang J, Witte ON, and Graeber TG (2019). Pan-cancer Convergence to a Small-Cell Neuroendocrine Phenotype that Shares Susceptibilities with Hematological Malignancies. Cancer Cell 36, 17–34 e17. 10.1016/j.ccell.2019.06.005. - DOI - PMC - PubMed
    1. Cejas P, Xie Y, Font-Tello A, Lim K, Syamala S, Qiu X, Tewari AK, Shah N, Nguyen HM, Patel RA, et al. (2021). Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer. Nat Commun 12, 5775. 10.1038/s41467-021-26042-z. - DOI - PMC - PubMed

Publication types

MeSH terms

Substances