Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci.

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

The long noncoding RNA H19 regulates tumor plasticity in neuroendocrine prostate cancer.

Neuroendocrine (NE) prostate cancer (NEPC) is a lethal subtype of castration-resistant prostate cancer (PCa) arising either de novo or from transdifferentiated prostate adenocarcinoma following androgen deprivation therapy (ADT). Extensive computational analysis has identified a high degree of association between the long noncoding RNA (lncRNA) H19 and NEPC, with the longest isoform highly expressed in NEPC. H19 regulates PCa lineage plasticity by driving a bidirectional cell identity of NE phenotype (H19 overexpression) or luminal phenotype (H19 knockdown). It contributes to treatment resistance, with the knockdown of H19 re-sensitizing PCa to ADT. It is also essential for the proliferation and invasion of NEPC. H19 levels are negatively regulated by androgen signaling via androgen receptor (AR). When androgen is absent SOX2 levels increase, driving H19 transcription and facilitating transdifferentiation. H19 facilitates the PRC2 complex in regulating methylation changes at H3K27me3/H3K4me3 histone sites of AR-driven and NEPC-related genes. Additionally, this lncRNA induces alterations in genome-wide DNA methylation on CpG sites, further regulating genes associated with the NEPC phenotype. Our clinical data identify H19 as a candidate diagnostic marker and predictive marker of NEPC with elevated H19 levels associated with an increased probability of biochemical recurrence and metastatic disease in patients receiving ADT. Here we report H19 as an early upstream regulator of cell fate, plasticity, and treatment resistance in NEPC that can reverse/transform cells to a treatable form of PCa once therapeutically deactivated.

Transient Sox9 Expression Facilitates Resistance to Androgen-Targeted Therapy in Prostate Cancer.

PURPOSE: Patients with metastatic prostate cancer are increasingly presenting with treatment-resistant, androgen receptor-negative/low (AR-/Low) tumors, with or without neuroendocrine characteristics, in processes attributed to tumor cell plasticity. This plasticity has been modeled by Rb1/p53 knockdown/knockout and is accompanied by overexpression of the pluripotency factor, Sox2. Here, we explore the role of the developmental transcription factor Sox9 in the process of prostate cancer therapy response and tumor progression. EXPERIMENTAL DESIGN: Unique prostate cancer cell models that capture AR-/Low stem cell-like intermediates were analyzed for features of plasticity and the functional role of Sox9. Human prostate cancer xenografts and tissue microarrays were evaluated for temporal alterations in Sox9 expression. The role of NF-κB pathway activity in Sox9 overexpression was explored. RESULTS: Prostate cancer stem cell-like intermediates have reduced Rb1 and p53 protein expression and overexpress Sox2 as well as Sox9. Sox9 was required for spheroid growth, and overexpression increased invasiveness and neural features of prostate cancer cells. Sox9 was transiently upregulated in castration-induced progression of prostate cancer xenografts and was specifically overexpressed in neoadjuvant hormone therapy (NHT)-treated patient tumors. High Sox9 expression in NHT-treated patients predicts biochemical recurrence. Finally, we link Sox9 induction to NF-κB dimer activation in prostate cancer cells. CONCLUSIONS: Developmentally reprogrammed prostate cancer cell models recapitulate features of clinically advanced prostate tumors, including downregulated Rb1/p53 and overexpression of Sox2 with Sox9. Sox9 is a marker of a transitional state that identifies prostate cancer cells under the stress of therapeutic assault and facilitates progression to therapy resistance. Its expression may index the relative activity of the NF-κB pathway.

Characterization of transcriptomic signature of primary prostate cancer analogous to prostatic small cell neuroendocrine carcinoma.

Prostatic small cell neuroendocrine carcinoma (SC/NE) is well studied in metastatic castration-resistant prostate cancer; however, it is not well characterized in the primary setting. Herein, we used gene expression profiling of SC/NE prostate cancer (PCa) to develop a 212 gene signature to identify treatment-naïve primary prostatic tumors that are molecularly analogous to SC/NE (SC/NE-like PCa). The 212 gene signature was tested in several cohorts confirming similar molecular profile between prostatic SC/NE and small cell lung carcinoma. The signature was then translated into a genomic score (SCGScore) using modularized logistic regression modeling and validated in four independent cohorts achieving an average AUC >0.95. The signature was evaluated in more than 25,000 primary adenocarcinomas to characterize the biology, prognosis and potential therapeutic response of predicted SC/NE-like tumors. Assessing SCGScore in a prospective cohort of 17,967 RP and 6,697 biopsy treatment-naïve primary tumors from the Decipher Genomic Resource Information Database registry, approximately 1% of the patients were found to have a SC/NE-like transcriptional profile, whereas 0.5 and 3% of GG1 and GG5 patients respectively showed to be SC/NE-like. More than 80% of these patients are genomically high-risk based on Decipher score. Interrogating in vitro drug sensitivity analyses, SC/NE-like prostatic tumors showed higher response to PARP and HDAC inhibitors.

Alternative RNA splicing of the GIT1 gene is associated with neuroendocrine prostate cancer.

Potent androgen receptor pathway inhibition (ARPI) therapies have given rise to a lethal, aggressive subtype of castration-resistant prostate cancer (CRPC) called treatment-induced neuroendocrine prostate cancer (t-NEPC). Now, t-NEPC poses a major clinical problem as approximately 20% of CRPC cases bear this subtype-a rate of occurrence that is predicted to rise with the widespread use of ARPI therapies. Unfortunately, there are no targeted therapies currently available to treat t-NEPC as the origin and molecular underpinnings of t-NEPC development remain unclear. In the present study, we identify that RNA splicing of the G protein-coupled receptor kinase-interacting protein 1 (GIT1) gene is a unique event in t-NEPC patients. Specifically, upregulation of the GIT1-A splice variant and downregulation of the GIT1-C variant expressions are associated with t-NEPC patient tumors, patient-derived xenografts, and cell models. RNA-binding assays show that RNA splicing of GIT1 is directly driven by SRRM4 and is associated with the neuroendocrine phenotype in CRPC cohorts. We show that GIT1-A and GIT1-C regulate differential transcriptomes in prostate cancer cells, where GIT1-A regulates genes associated with morphogenesis, neural function, environmental sensing via cell-adhesion processes, and epigenetic regulation. Consistent with our transcriptomic analyses, we report opposing functions of GIT1-A and GIT1-C in the stability of focal adhesions, whereby GIT1-A promotes focal adhesion stability. In summary, our study is the first to report that alternative RNA splicing of the GIT1 gene is associated with t-NEPC and reprograms its function involving FA-mediated signaling and cell processes, which may contribute to t-NEPC development.

MechRNA: prediction of lncRNA mechanisms from RNA-RNA and RNA-protein interactions.

Motivation: Long non-coding RNAs (lncRNAs) are defined as transcripts longer than 200 nt that do not get translated into proteins. Often these transcripts are processed (spliced, capped and polyadenylated) and some are known to have important biological functions. However, most lncRNAs have unknown or poorly understood functions. Nevertheless, because of their potential role in cancer, lncRNAs are receiving a lot of attention, and the need for computational tools to predict their possible mechanisms of action is more than ever. Fundamentally, most of the known lncRNA mechanisms involve RNA-RNA and/or RNA-protein interactions. Through accurate predictions of each kind of interaction and integration of these predictions, it is possible to elucidate potential mechanisms for a given lncRNA. Results: Here, we introduce MechRNA, a pipeline for corroborating RNA-RNA interaction prediction and protein binding prediction for identifying possible lncRNA mechanisms involving specific targets or on a transcriptome-wide scale. The first stage uses a version of IntaRNA2 with added functionality for efficient prediction of RNA-RNA interactions with very long input sequences, allowing for large-scale analysis of lncRNA interactions with little or no loss of optimality. The second stage integrates protein binding information pre-computed by GraphProt, for both the lncRNA and the target. The final stage involves inferring the most likely mechanism for each lncRNA/target pair. This is achieved by generating candidate mechanisms from the predicted interactions, the relative locations of these interactions and correlation data, followed by selection of the most likely mechanistic explanation using a combined P-value. We applied MechRNA on a number of recently identified cancer-related lncRNAs (PCAT1, PCAT29 and ARLnc1) and also on two well-studied lncRNAs (PCA3 and 7SL). This led to the identification of hundreds of high confidence potential targets for each lncRNA and corresponding mechanisms. These predictions include the known competitive mechanism of 7SL with HuR for binding on the tumor suppressor TP53, as well as mechanisms expanding what is known about PCAT1 and ARLn1 and their targets BRCA2 and AR, respectively. For PCAT1-BRCA2, the mechanism involves competitive binding with HuR, which we confirmed using HuR immunoprecipitation assays. Availability and implementation: MechRNA is available for download at https://bitbucket.org/compbio/mechrna. Supplementary information: Supplementary data are available at Bioinformatics online.

Differentially expressed microRNAs in lung adenocarcinoma invert effects of copy number aberrations of prognostic genes.

In many cancers, significantly down- or upregulated genes are found within chromosomal regions with DNA copy number alteration opposite to the expression changes. Generally, this paradox has been overlooked as noise, but can potentially be a consequence of interference of epigenetic regulatory mechanisms, including microRNA-mediated control of mRNA levels. To explore potential associations between microRNAs and paradoxes in non-small-cell lung cancer (NSCLC) we curated and analyzed lung adenocarcinoma (LUAD) data, comprising gene expressions, copy number aberrations (CNAs) and microRNA expressions. We integrated data from 1,062 tumor samples and 241 normal lung samples, including newly-generated array comparative genomic hybridization (aCGH) data from 63 LUAD samples. We identified 85 "paradoxical" genes whose differential expression consistently contrasted with aberrations of their copy numbers. Paradoxical status of 70 out of 85 genes was validated on sample-wise basis using The Cancer Genome Atlas (TCGA) LUAD data. Of these, 41 genes are prognostic and form a clinically relevant signature, which we validated on three independent datasets. By meta-analysis of results from 9 LUAD microRNA expression studies we identified 24 consistently-deregulated microRNAs. Using TCGA-LUAD data we showed that deregulation of 19 of these microRNAs explains differential expression of the paradoxical genes. Our results show that deregulation of paradoxical genes is crucial in LUAD and their expression pattern is maintained epigenetically, defying gene copy number status.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND: Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS: We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS: Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION: This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING: Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

NBN gain is predictive for adverse outcome following image-guided radiotherapy for localized prostate cancer.

Despite the use of clinical prognostic factors (PSA, T-category and Gleason score), 20-60% of localized prostate cancers (PCa) fail primary local treatment. Herein, we determined the prognostic importance of main sensors of the DNA damage response (DDR): MRE11A, RAD50, NBN, ATM, ATR and PRKDC. We studied copy number alterations in DDR genes in localized PCa treated with image-guided radiotherapy (IGRT; n=139) versus radical prostatectomy (RadP; n=154). In both cohorts, NBN gains were the most frequent genomic alteration (14.4 and 11% of cases, respectively), and were associated with overall tumour genomic instability (p<0.0001). NBN gains were the only significant predictor of 5yrs biochemical relapse-free rate (bRFR) following IGRT (46% versus 77%; p=0.00067). On multivariate analysis, NBN gain remained a significant independent predictor of bRFR after adjusting for known clinical prognostic variables (HR=3.28, 95% CI 1.56-6.89, Wald p-value=0.0017). No DDR-sensing gene was prognostic in the RadP cohort. In vitro studies correlated NBN gene overexpression with PCa cells radioresistance. In conclusion, NBN gain predicts for decreased bRFR in IGRT, but not in RadP patients. If validated independently, Nibrin gains may be the first PCa predictive biomarker to facilitate local treatment decisions using precision medicine approaches with surgery or radiotherapy.

Genetic alterations defining NSCLC subtypes and their therapeutic implications.

Lung cancer is the leading cause of cancer death worldwide, accounting for more deaths than breast, prostate and colon cancer combined. While treatment decisions are determined primarily by stage, therapeutically non small cell lung cancer (NSCLC) has traditionally been treated as a single disease. However, recent findings have led to the recognition of histology and molecular subtypes as important determinants in treatment selection. Identifying the genetic differences that define these molecular and histological subtypes has the potential to impact treatment and as such is currently the focus of much research. Microarray and genomic sequencing efforts have provided unparalleled insight into the genomes of lung cancer subtypes, specifically adenocarcinoma (AC) and squamous cell carcinoma (SqCC), revealing subtype specific genomic alterations and molecular subtypes as well as differences in cell signaling pathways. In this review, we discuss the recurrent genomic alterations characteristic of AC and SqCC (including molecular subtypes), their therapeutic implications and emerging clinical practices aimed at tailoring treatments based on a tumor's molecular alterations with the hope of improving patient response and survival.