Supporting data for "The Long Noncoding RNA Landscape of Neuroendocrine Prostate Cancer and its Clinical Implications"

Treatment induced neuroendocrine prostate cancer (tNEPC) is an aggressive variant of late-stage metastatic castrate resistant (mCRPC) prostate cancer that commonly arises through neuroendocrine transdifferentiation (NEtD). Treatment options are limited, ineffective, and for most patients, results in death in less than a year. We previously developed a first-in-field patient-derived xenograft (PDX) model of NEtD. Longitudinal deep transcriptome profiling of this model enabled monitoring of dynamic transcriptional changes during NEtD and in the context of androgen deprivation. Long non-coding RNA (lncRNA) are implicated in cancer where they can control gene regulation. Until now the expression of lncRNAs during NEtD and their clinical associations were unexplored. We implemented a next-generation sequence analysis pipeline that can detect transcripts at low expression levels and built a genome-wide catalogue (n=37,749) oflncRNAs. We applied this pipeline to 927 clinical samples and our high fidelity NEtD model LTL331 and identified 821 lncRNAs in NEPC. Among these are 122 lncRNAs that robustly distinguish NEPC from prostate adenocarcinoma (AD) patient tumours. The highest expressed lncRNAs within this signature are H19, LINC00617, and SSTR5-AS1. Another 742 are associated with the NEtD process and fall into four distinct patterns of expression (NEtD lncRNA Class I, II, III, and IV) in our PDX model and clinical samples. Each class has significant (z-scores>2) and unique enrichment for transcription factor binding site (TFBS) motifs in their sequences. Enriched TFBS include (1) TP53 and BRN1 in Class I, (2) ELF5, SPIC, and HOXD1 in Class II, (3) SPDEF in Class III, (4) HSF1 and FOXA1 in Class IV, and (5) TWIST1 when merging Class III with IV. Common TFBS in all NEtD lncRNA were also identified and include, E2F, REST, PAX5, PAX9, and STAF. Interrogation of the top deregulated candidates (n=100) in radical prostatectomy adenocarcinoma samples with long-term follow-up (median 18 years) revealed significant clinicopathological associations. Specifically, we identified 25 that are associated with rapid metastasis following androgen deprivation therapy (ADT). Two of these lncRNAs (SSTR5-AS1 and LINC00514) stratified patients undergoing ADT based on patient outcome. To date, a comprehensive characterization of the dynamic landscape of lncRNAs during the NEtD process has not been performed. A temporal analysis of the PDXbased NEtD model has for the first time provide