The mutational landscape of lethal castration-resistant prostate cancer.

Characterization of the prostate cancer transcriptome and genome has identified chromosomal rearrangements and copy number gains and losses, including ETS gene family fusions, PTEN loss and androgen receptor (AR) amplification, which drive prostate cancer development and progression to lethal, metastatic castration-resistant prostate cancer (CRPC). However, less is known about the role of mutations. Here we sequenced the exomes of 50 lethal, heavily pre-treated metastatic CRPCs obtained at rapid autopsy (including three different foci from the same patient) and 11 treatment-naive, high-grade localized prostate cancers. We identified low overall mutation rates even in heavily treated CRPCs (2.00 per megabase) and confirmed the monoclonal origin of lethal CRPC. Integrating exome copy number analysis identified disruptions of CHD1 that define a subtype of ETS gene family fusion-negative prostate cancer. Similarly, we demonstrate that ETS2, which is deleted in approximately one-third of CRPCs (commonly through TMPRSS2:ERG fusions), is also deregulated through mutation. Furthermore, we identified recurrent mutations in multiple chromatin- and histone-modifying genes, including MLL2 (mutated in 8.6% of prostate cancers), and demonstrate interaction of the MLL complex with the AR, which is required for AR-mediated signalling. We also identified novel recurrent mutations in the AR collaborating factor FOXA1, which is mutated in 5 of 147 (3.4%) prostate cancers (both untreated localized prostate cancer and CRPC), and showed that mutated FOXA1 represses androgen signalling and increases tumour growth. Proteins that physically interact with the AR, such as the ERG gene fusion product, FOXA1, MLL2, UTX (also known as KDM6A) and ASXL1 were found to be mutated in CRPC. In summary, we describe the mutational landscape of a heavily treated metastatic cancer, identify novel mechanisms of AR signalling deregulated in prostate cancer, and prioritize candidates for future study.

Growth factor stimulation induces a distinct ER(alpha) cistrome underlying breast cancer endocrine resistance.

Estrogen receptor α (ERα) expression in breast cancer is predictive of response to endocrine therapy; however, resistance is common in ERα-positive tumors that overexpress the growth factor receptor ERBB2. Even in the absence of estrogen, ERα can be activated by growth factors, including the epidermal growth factor (EGF). EGF induces a transcriptional program distinct from estrogen; however, the mechanism of the stimulus-specific response is unknown. Here we show that the EGF-induced ERα genomic targets, its cistromes, are distinct from those induced by estrogen in a process dependent on the transcription factor AP-1. The EGF-induced ERα cistrome specifically regulates genes found overexpressed in ERBB2-positive human breast cancers. This provides a potential molecular explanation for the endocrine therapy resistance seen in ERα-positive breast cancers that overexpress ERBB2. These results suggest a central role for ERα in hormone-refractory breast tumors dependent on growth factor pathway activation and favors the development of therapeutic strategies completely antagonizing ERα, as opposed to blocking its estrogen responsiveness alone.

Integrative molecular concept modeling of prostate cancer progression.

Despite efforts to profile prostate cancer, the genetic alterations and biological processes that correlate with the observed histological progression are unclear. Using laser-capture microdissection to isolate 101 cell populations, we have profiled prostate cancer progression from benign epithelium to metastatic disease. By analyzing expression signatures in the context of over 14,000 'molecular concepts', or sets of biologically connected genes, we generated an integrative model of progression. Molecular concepts that demarcate critical transitions in progression include protein biosynthesis, E26 transformation-specific (ETS) family transcriptional targets, androgen signaling and cell proliferation. Of note, relative to low-grade prostate cancer (Gleason pattern 3), high-grade cancer (Gleason pattern 4) shows an attenuated androgen signaling signature, similar to metastatic prostate cancer, which may reflect dedifferentiation and explain the clinical association of grade with prognosis. Taken together, these data show that analyzing gene expression signatures in the context of a compendium of molecular concepts is useful in understanding cancer biology.

Integrative analysis of the cancer transcriptome.

DNA microarrays have been widely applied to the study of human cancer, delineating myriad molecular subtypes of cancer, many of which are associated with distinct biological underpinnings, disease progression and treatment response. These primary analyses have begun to decipher the molecular heterogeneity of cancer, but integrative analyses that evaluate cancer transcriptome data in the context of other data sources are often capable of extracting deeper biological insight from the data. Here we discuss several such integrative computational and analytical approaches, including meta-analysis, functional enrichment analysis, interactome analysis, transcriptional network analysis and integrative model system analysis.

Mining for regulatory programs in the cancer transcriptome.

DNA microarrays have been widely applied to cancer transcriptome analysis. The Oncomine database contains a large collection of such data, as well as hundreds of derived gene-expression signatures. We studied the regulatory mechanisms responsible for gene deregulation in these cancer signatures by searching for the coordinate regulation of genes with common transcription factor binding sites. We found that genes with binding sites for the archetypal cancer transcription factor, E2F, were disproportionately overexpressed in a wide variety of cancers, whereas genes with binding sites for other transcription factors, such as Myc-Max, c-Rel and ATF, were disproportionately overexpressed in specific cancer types. These results suggest that alterations in pathways activating these transcription factors may be responsible for the observed gene deregulation and cancer pathogenesis.

Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression.

Molecular profiling of cancer at the transcript level has become routine. Large-scale analysis of proteomic alterations during cancer progression has been a more daunting task. Here, we employed high-throughput immunoblotting in order to interrogate tissue extracts derived from prostate cancer. We identified 64 proteins that were altered in prostate cancer relative to benign prostate and 156 additional proteins that were altered in metastatic disease. An integrative analysis of this compendium of proteomic alterations and transcriptomic data was performed, revealing only 48%-64% concordance between protein and transcript levels. Importantly, differential proteomic alterations between metastatic and clinically localized prostate cancer that mapped concordantly to gene transcripts served as predictors of clinical outcome in prostate cancer as well as other solid tumors.

Validation of Immunotherapy Response Score as Predictive of Pan-solid Tumor Anti-PD-1/PD-L1 Benefit.

Immunotherapy response score (IRS) integrates tumor mutation burden (TMB) and quantitative expression biomarkers to predict anti-PD-1/PD-L1 [PD-(L)1] monotherapy benefit. Here, we evaluated IRS in additional cohorts. Patients from an observational trial (NCT03061305) treated with anti-PD-(L)1 monotherapy were included and assigned to IRS-High (-H) versus -Low (-L) groups. Associations with real-world progression-free survival (rwPFS) and overall survival (OS) were determined by Cox proportional hazards (CPH) modeling. Those with available PD-L1 IHC treated with anti-PD-(L)1 with or without chemotherapy were separately assessed. Patients treated with PD-(L)1 and/or chemotherapy (five relevant tumor types) were assigned to three IRS groups [IRS-L divided into IRS-Ultra-Low (-UL) and Intermediate-Low (-IL), and similarly assessed]. In the 352 patient anti-PD-(L)1 monotherapy validation cohort (31 tumor types), IRS-H versus IRS-L patients had significantly longer rwPFS and OS. IRS significantly improved CPH associations with rwPFS and OS beyond microsatellite instability (MSI)/TMB alone. In a 189 patient (10 tumor types) PD-L1 IHC comparison cohort, IRS, but not PD-L1 IHC nor TMB, was significantly associated with anti-PD-L1 rwPFS. In a 1,103-patient cohort (from five relevant tumor types), rwPFS did not significantly differ in IRS-UL patients treated with chemotherapy versus chemotherapy plus anti-PD-(L)1, nor in IRS-H patients treated with anti-PD-(L)1 versus anti-PD-(L)1 + chemotherapy. IRS associations were consistent across subgroups, including both Europeans and non-Europeans. These results confirm the utility of IRS utility for predicting pan-solid tumor PD-(L)1 monotherapy benefit beyond available biomarkers and demonstrate utility for informing on anti-PD-(L)1 and/or chemotherapy treatment. Significance: This study confirms the utility of the integrative IRS biomarker for predicting anti-PD-L1/PD-1 benefit. IRS significantly improved upon currently available biomarkers, including PD-L1 IHC, TMB, and MSI status. Additional utility for informing on chemotherapy, anti-PD-L1/PD-1, and anti-PD-L1/PD-1 plus chemotherapy treatments decisions is shown.

Development and validation of an integrative pan-solid tumor predictor of PD-1/PD-L1 blockade benefit.

BACKGROUND: Anti-PD-1 and PD-L1 (collectively PD-[L]1) therapies are approved for many advanced solid tumors. Biomarkers beyond PD-L1 immunohistochemistry, microsatellite instability, and tumor mutation burden (TMB) may improve benefit prediction. METHODS: Using treatment data and genomic and transcriptomic tumor tissue profiling from an observational trial (NCT03061305), we developed Immunotherapy Response Score (IRS), a pan-tumor predictive model of PD-(L)1 benefit. IRS real-world progression free survival (rwPFS) and overall survival (OS) prediction was validated in an independent cohort of trial patients. RESULTS: Here, by Cox modeling, we develop IRS-which combines TMB with CD274, PDCD1, ADAM12 and TOP2A quantitative expression-to predict pembrolizumab rwPFS (648 patients; 26 tumor types; IRS-High or -Low groups). In the 248 patient validation cohort (248 patients; 24 tumor types; non-pembrolizumab PD-[L]1 monotherapy treatment), median rwPFS and OS are significantly longer in IRS-High vs. IRS-Low patients (rwPFS adjusted hazard ratio [aHR] 0.52, p = 0.003; OS aHR 0.49, p = 0.005); TMB alone does not significantly predict PD-(L)1 rwPFS nor OS. In 146 patients treated with systemic therapy prior to pembrolizumab monotherapy, pembrolizumab rwPFS is only significantly longer than immediately preceding therapy rwPFS in IRS-High patients (interaction test p = 0.001). In propensity matched lung cancer patients treated with first-line pembrolizumab monotherapy or pembrolizumab+chemotherapy, monotherapy rwPFS is significantly shorter in IRS-Low patients, but is not significantly different in IRS-High patients. Across 24,463 molecularly-evaluable trial patients, 7.6% of patients outside of monotherapy PD-(L)1 approved tumor types are IRS-High/TMB-Low. CONCLUSIONS: The validated, predictive, pan-tumor IRS model can expand PD-(L)1 monotherapy benefit outside currently approved indications.

Therapies activating the immune system (checkpoint inhibitors) have revolutionized the treatment of patients with advanced cancer, however new molecular tests may better identify patients who could benefit. Using treatment data and clinical molecular test results, we report the development and validation of Immunotherapy Response Score (IRS) to predict checkpoint inhibitor benefit. Across patients with more than 20 advanced cancer types, IRS better predicted checkpoint inhibitor benefit than currently available tests. Data from >20,000 patients showed that IRS identifies ~8% of patients with advanced cancer who may dramatically benefit from checkpoint inhibitors but would not receive them today based on currently available tests. Our approach may help clinicians to decide which patients should receive checkpoint inhibitors to treat their disease.

AGTR1 overexpression defines a subset of breast cancer and confers sensitivity to losartan, an AGTR1 antagonist.

Breast cancer patients have benefited from the use of targeted therapies directed at specific molecular alterations. To identify additional opportunities for targeted therapy, we searched for genes with marked overexpression in subsets of tumors across a panel of breast cancer profiling studies comprising 3,200 microarray experiments. In addition to prioritizing ERBB2, we found AGTR1, the angiotensin II receptor type I, to be markedly overexpressed in 10-20% of breast cancer cases across multiple independent patient cohorts. Validation experiments confirmed that AGTR1 is highly overexpressed, in several cases more than 100-fold. AGTR1 overexpression was restricted to estrogen receptor-positive tumors and was mutually exclusive with ERBB2 overexpression across all samples. Ectopic overexpression of AGTR1 in primary mammary epithelial cells, combined with angiotensin II stimulation, led to a highly invasive phenotype that was attenuated by the AGTR1 antagonist losartan. Similarly, losartan reduced tumor growth by 30% in AGTR1-positive breast cancer xenografts. Taken together, these observations indicate that marked AGTR1 overexpression defines a subpopulation of ER-positive, ERBB2-negative breast cancer that may benefit from targeted therapy with AGTR1 antagonists, such as losartan.

Development and Validation of StrataNGS, a Multiplex PCR, Semiconductor Sequencing-Based Comprehensive Genomic Profiling Test.

Despite widespread use in targeted tumor testing, multiplex PCR/semiconductor (Ion Torrent) sequencing-based assessment of all comprehensive genomic profiling (CGP) variant classes has been limited. Herein, we describe the development and validation of StrataNGS, a 429-gene, multiplex PCR/semiconductor sequencing-based CGP laboratory-developed test performed on co-isolated DNA and RNA from formalin-fixed, paraffin-embedded tumor specimens with ≥2 mm2 tumor surface area. Validation was performed in accordance with MolDX CGP validation guidelines using 1986 clinical formalin-fixed, paraffin-embedded samples and an in-house developed optimized bioinformatics pipeline. Across CGP variant classes, accuracy ranged from 0.945 for tumor mutational burden (TMB) status to >0.999 for mutations and gene fusions, positive predictive value ranged from 0.915 for TMB status to 1.00 for gene fusions, and reproducibility ranged from 0.998 for copy number alterations to 1.00 for splice variants and insertions/deletions. StrataNGS TMB estimates were highly correlated to those from whole exome- or FoundationOne CDx-determined TMB (Pearson r = 0.998 and 0.960, respectively); TMB reproducibility was 0.996 (concordance correlation coefficient). Limit of detection for all variant classes was <20% tumor content. Together, we demonstrate that multiplex PCR/semiconductor sequencing-based tumor tissue CGP is feasible using optimized bioinformatic approaches described herein.