Canonical AREs are tumor suppressive regulatory elements in the prostate.

The androgen receptor (AR) is the central determinant of prostate tissue identity and differentiation, controlling normal, growth-suppressive prostate-specific gene expression 1 . It is also a key driver of prostate tumorigenesis, becoming "hijacked" to drive oncogenic transcription 2-5 . However, the regulatory elements determining the execution of the growth suppressive AR transcriptional program, and whether this can be reactivated in prostate cancer (PCa) cells remains unclear. Canonical androgen response element (ARE) motifs are the classic DNA binding element for AR 6 . Here, we used a genome-wide strategy to modulate regulatory elements containing AREs to define distinct AR transcriptional programs. We find that activation of these AREs is specifically associated with differentiation and growth suppressive transcription, and this can be reactivated to cause death in AR + PCa cells. In contrast, repression of AREs is well tolerated by PCa cells, but deleterious to normal prostate cells. Finally, gene expression signatures driven by ARE activity are associated with improved prognosis and luminal phenotypes in human PCa patients. This study demonstrates that canonical AREs are responsible for a normal, growth-suppressive, lineage-specific transcriptional program, that this can be reengaged in PCa cells for potential therapeutic benefit, and genes controlled by this mechanism are clinically relevant in human PCa patients.

A novel prostate cancer subtyping classifier based on luminal and basal phenotypes.

BACKGROUND: Prostate cancer (PCa) is a clinically heterogeneous disease. The creation of an expression-based subtyping model based on prostate-specific biological processes was sought. METHODS: Unsupervised machine learning of gene expression profiles from prospectively collected primary prostate tumors (training, n = 32,000; evaluation, n = 68,547) was used to create a prostate subtyping classifier (PSC) based on basal versus luminal cell expression patterns and other gene signatures relevant to PCa biology. Subtype molecular pathways and clinical characteristics were explored in five other clinical cohorts. RESULTS: Clustering derived four subtypes: luminal differentiated (LD), luminal proliferating (LP), basal immune (BI), and basal neuroendocrine (BN). LP and LD tumors both had higher androgen receptor activity. LP tumors also had a higher expression of cell proliferation genes, MYC activity, and characteristics of homologous recombination deficiency. BI tumors possessed significant interferon γactivity and immune infiltration on immunohistochemistry. BN tumors were characterized by lower androgen receptor activity expression, lower immune infiltration, and enrichment with neuroendocrine expression patterns. Patients with LD tumors had less aggressive tumor characteristics and the longest time to metastasis after surgery. Only patients with BI tumors derived benefit from radiotherapy after surgery in terms of time to metastasis (hazard ratio [HR], 0.09; 95% CI, 0.01-0.71; n = 855). In a phase 3 trial that randomized patients with metastatic PCa to androgen deprivation with or without docetaxel (n = 108), only patients with LP tumors derived survival benefit from docetaxel (HR, 0.21; 95% CI, 0.09-0.51). CONCLUSIONS: With the use of expression profiles from over 100,000 tumors, a PSC was developed that identified four subtypes with distinct biological and clinical features. PLAIN LANGUAGE SUMMARY: Prostate cancer can behave in an indolent or aggressive manner and vary in how it responds to certain treatments. To differentiate prostate cancer on the basis of biological features, we developed a novel RNA signature by using data from over 100,000 prostate tumors-the largest data set of its kind. This signature can inform patients and physicians on tumor aggressiveness and susceptibilities to treatments to help personalize cancer management.

Genomic Testing in Localized Prostate Cancer Can Identify Subsets of African Americans With Aggressive Disease.

BACKGROUND: Personalized genomic classifiers have transformed the management of prostate cancer (PCa) by identifying the most aggressive subsets of PCa. Nevertheless, the performance of genomic classifiers to risk classify African American men is thus far lacking in a prospective setting. METHODS: This is a prospective study of the Decipher genomic classifier for National Comprehensive Cancer Network low- and intermediate-risk PCa. Study-eligible non-African American men were matched to African American men. Diagnostic biopsy specimens were processed to estimate Decipher scores. Samples accrued in NCT02723734, a prospective study, were interrogated to determine the genomic risk of reclassification (GrR) between conventional clinical risk classifiers and the Decipher score. RESULTS: The final analysis included a clinically balanced cohort of 226 patients with complete genomic information (113 African American men and 113 non-African American men). A higher proportion of African American men with National Comprehensive Cancer Network-classified low-risk (18.2%) and favorable intermediate-risk (37.8%) PCa had a higher Decipher score than non-African American men. Self-identified African American men were twice more likely than non-African American men to experience GrR (relative risk [RR] = 2.23, 95% confidence interval [CI] = 1.02 to 4.90; P = .04). In an ancestry-determined race model, we consistently validated a higher risk of reclassification in African American men (RR = 5.26, 95% CI = 1.66 to 16.63; P = .004). Race-stratified analysis of GrR vs non-GrR tumors also revealed molecular differences in these tumor subtypes. CONCLUSIONS: Integration of genomic classifiers with clinically based risk classification can help identify the subset of African American men with localized PCa who harbor high genomic risk of early metastatic disease. It is vital to identify and appropriately risk stratify the subset of African American men with aggressive disease who may benefit from more targeted interventions.

Naturally occurring combinations of receptors from single cell transcriptomics in endothelial cells.

VEGF inhibitor drugs are part of standard care in oncology and ophthalmology, but not all patients respond to them. Combinations of drugs are likely to be needed for more effective therapies of angiogenesis-related diseases. In this paper we describe naturally occurring combinations of receptors in endothelial cells that might help to understand how cells communicate and to identify targets for drug combinations. We also develop and share a new software tool called DECNEO to identify them. Single-cell gene expression data are used to identify a set of co-expressed endothelial cell receptors, conserved among species (mice and humans) and enriched, within a network, of connections to up-regulated genes. This set includes several receptors previously shown to play a role in angiogenesis. Multiple statistical tests from large datasets, including an independent validation set, support the reproducibility, evolutionary conservation and role in angiogenesis of these naturally occurring combinations of receptors. We also show tissue-specific combinations and, in the case of choroid endothelial cells, consistency with both well-established and recent experimental findings, presented in a separate paper. The results and methods presented here advance the understanding of signaling to endothelial cells. The methods are generally applicable to the decoding of intercellular combinations of signals.

Digital Cell Sorter (DCS): a cell type identification, anomaly detection, and Hopfield landscapes toolkit for single-cell transcriptomics.

MOTIVATION: Analysis of singe cell RNA sequencing (scRNA-seq) typically consists of different steps including quality control, batch correction, clustering, cell identification and characterization, and visualization. The amount of scRNA-seq data is growing extremely fast, and novel algorithmic approaches improving these steps are key to extract more biological information. Here, we introduce: (i) two methods for automatic cell type identification (i.e., without expert curator) based on a voting algorithm and a Hopfield classifier, (ii) a method for cell anomaly quantification based on isolation forest, and (iii) a tool for the visualization of cell phenotypic landscapes based on Hopfield energy-like functions. These new approaches are integrated in a software platform that includes many other state-of-the-art methodologies and provides a self-contained toolkit for scRNA-seq analysis. RESULTS: We present a suite of software elements for the analysis of scRNA-seq data. This Python-based open source software, Digital Cell Sorter (DCS), consists in an extensive toolkit of methods for scRNA-seq analysis. We illustrate the capability of the software using data from large datasets of peripheral blood mononuclear cells (PBMC), as well as plasma cells of bone marrow samples from healthy donors and multiple myeloma patients. We test the novel algorithms by evaluating their ability to deconvolve cell mixtures and detect small numbers of anomalous cells in PBMC data. AVAILABILITY: The DCS toolkit is available for download and installation through the Python Package Index (PyPI). The software can be deployed using the Python import function following installation. Source code is also available for download on Zenodo: DOI 10.5281/zenodo.2533377. SUPPLEMENTARY INFORMATION: Supplemental Materials are available at PeerJ online.

Modeling disease progression in Multiple Myeloma with Hopfield networks and single-cell RNA-seq.

Associative memories in Hopfield's neural networks are mapped to gene expression pattern to model different paths of disease progression towards Multiple Myeloma (MM). The model is built using single cell RNA-seq data from bone marrow aspirates of MM patients as well as patients diagnosed with Monoclonal Gammopathy of Undetermined Significance (MGUS) and Smoldering Multiple Myeloma (SMM), two medical conditions that often progress to full MM. Results: We identify different clusters of MGUS, SMM, and MM cells, map them to Hopfield associative memory patterns, and model the dynamics of transition between the different patterns. The model is then used to identify genes that are differentialy expressed across different MM stages and whose simultaneous inhibition is associated to a delayed disease progression.