Robustness of differential gene expression analysis of RNA-seq.

RNA-sequencing (RNA-seq) is a relatively new technology that lacks standardisation. RNA-seq can be used for Differential Gene Expression (DGE) analysis, however, no consensus exists as to which methodology ensures robust and reproducible results. Indeed, it is broadly acknowledged that DGE methods provide disparate results. Despite obstacles, RNA-seq assays are in advanced development for clinical use but further optimisation will be needed. Herein, five DGE models (DESeq2, voom + limma, edgeR, EBSeq, NOISeq) for gene-level detection were investigated for robustness to sequencing alterations using a controlled analysis of fixed count matrices. Two breast cancer datasets were analysed with full and reduced sample sizes. DGE model robustness was compared between filtering regimes and for different expression levels (high, low) using unbiased metrics. Test sensitivity estimated as relative False Discovery Rate (FDR), concordance between model outputs and comparisons of a 'population' of slopes of relative FDRs across different library sizes, generated using linear regressions, were examined. Patterns of relative DGE model robustness proved dataset-agnostic and reliable for drawing conclusions when sample sizes were sufficiently large. Overall, the non-parametric method NOISeq was the most robust followed by edgeR, voom, EBSeq and DESeq2. Our rigorous appraisal provides information for method selection for molecular diagnostics. Metrics may prove useful towards improving the standardisation of RNA-seq for precision medicine.

Clinically Actionable Insights into Initial and Matched Recurrent Glioblastomas to Inform Novel Treatment Approaches.

Glioblastoma is the most common primary adult brain tumour, and despite optimal treatment, the median survival is 12-15 months. Patients with matched recurrent glioblastomas were investigated to try to find actionable mutations. Tumours were profiled using a validated DNA-based gene panel. Copy number variations (CNVs) and single nucleotide variants (SNVs) were examined, and potentially pathogenic variants and clinically actionable mutations were identified. The results revealed that glioblastomas were IDH-wildtype (IDH WT; n = 38) and IDH-mutant (IDH MUT; n = 3). SNVs in TSC2, MSH6, TP53, CREBBP, and IDH1 were variants of unknown significance (VUS) that were predicted to be pathogenic in both subtypes. IDH WT tumours had SNVs that impacted RTK/Ras/PI(3)K, p53, WNT, SHH, NOTCH, Rb, and G-protein pathways. Many tumours had BRCA1/2 (18%) variants, including confirmed somatic mutations in haemangioblastoma. IDH WT recurrent tumours had fewer pathways impacted (RTK/Ras/PI(3)K, p53, WNT, and G-protein) and CNV gains (BRCA2, GNAS, and EGFR) and losses (TERT and SMARCA4). IDH MUT tumours had SNVs that impacted RTK/Ras/PI(3)K, p53, and WNT pathways. VUS in KLK1 was possibly pathogenic in IDH MUT. Recurrent tumours also had fewer pathways (p53, WNT, and G-protein) impacted by genetic alterations. Public datasets (TCGA and GDC) confirmed the clinical significance of findings in both subtypes. Overall in this cohort, potentially actionable variation was most often identified in EGFR, PTEN, BRCA1/2, and ATM. This study underlines the need for detailed molecular profiling to identify individual GBM patients who may be eligible for novel treatment approaches. This information is also crucial for patient recruitment to clinical trials.

Transcriptional subtyping and CD8 immunohistochemistry identifies poor prognosis stage II/III colorectal cancer patients who benefit from adjuvant chemotherapy.

PURPOSE: Transcriptomic profiling of colorectal cancer (CRC) has led to identification of four consensus molecular subtypes (CMS1-4), which have prognostic value in stage II/III disease. More recently, the Colorectal Cancer Intrinsic Subtypes (CRIS) classification system has helped to define the biology specific to the epithelial component of colorectal tumors. However, the clinical value of these classifications in predicting response to standard-of-care adjuvant chemotherapy remains unknown. PATIENTS AND METHODS: Using samples from 4 European sites, we assembled a novel stage II/III CRC patient cohort and performed transcriptomic profiling on 156 samples, targeted sequencing and generated a tissue microarray to enable integrated "multi-omics" analyses. We also accessed data from 2 published stage II/III CRC patient cohorts: GSE39582 and GSE14333 (479 and 185 samples respectively). RESULTS: The epithelial-rich CMS2 subtype of CRC benefitted significantly from adjuvant chemotherapy treatment in both stage II and III disease (p=0.02 and p<0.0001 respectively), while the CMS3 subtype significantly benefitted in stage III only (p=0.00073). Following CRIS sub-stratification of CMS2, we observed that only the CRIS-C subtype significantly benefitted from adjuvant chemotherapy in stage II and III disease (p=0.0081 and p<0.0001 respectively), while CRIS-D significantly benefitted in stage III only (p=0.0034). We also observed that CRIS-C patients with low levels of CD8+ tumor-infiltrating lymphocytes were most at risk of relapse in both stage II and III disease (p=0.0031). CONCLUSION: Patient stratification using a combination of transcriptional subtyping and CD8 immunohistochemistry analyses is capable of identifying poor prognostic stage II/III patients who benefit from adjuvant standard-of-care chemotherapy. These findings are particularly relevant for stage II disease, where the overall benefit of adjuvant chemotherapy is marginal.

Validation of a Metastatic Assay using biopsies to improve risk stratification in patients with prostate cancer treated with radical radiation therapy.

Background: Radiotherapy is an effective treatment of intermediate/high-risk locally advanced prostate cancer, however, >30% of patients relapse within 5 years. Clinicopathological parameters currently fail to identify patients prone to systemic relapse and those whom treatment intensification may be beneficial. The purpose of this study was to independently validate the performance of a 70-gene Metastatic Assay in a cohort of diagnostic biopsies from patients treated with radical radiotherapy and androgen deprivation therapy. Patients and methods: A bridging cohort of prostate cancer diagnostic biopsy specimens was profiled to enable optimization of the Metastatic Assay threshold before further independent clinical validation in a cohort of diagnostic biopsies from patients treated with radical radiotherapy and androgen deprivation therapy. Multivariable Cox proportional hazard regression analysis was used to assess assay performance in predicting biochemical failure-free survival (BFFS) and metastasis-free survival (MFS). Results: Gene expression analysis was carried out in 248 patients from the independent validation cohort and the Metastatic Assay applied. Ten-year MFS was 72% for Metastatic Assay positive patients and 94% for Metastatic Assay negative patients [HR = 3.21 (1.35-7.67); P = 0.003]. On multivariable analysis the Metastatic Assay remained predictive for development of distant metastases [HR = 2.71 (1.11-6.63); P = 0.030]. The assay retained independent prognostic performance for MFS when assessed with the Cancer of the Prostate Assessment Score (CAPRA) [HR = 3.23 (1.22-8.59); P = 0.019] whilst CAPRA itself was not significant [HR = 1.88, (0.52-6.77); P = 0.332]. A high concordance [100% (61.5-100)] for the assay result was noted between two separate foci taken from 11 tumours, whilst Gleason score had low concordance. Conclusions: The Metastatic Assay demonstrated significant prognostic performance in patients treated with radical radiotherapy both alone and independent of standard clinical and pathological variables. The Metastatic Assay could have clinical utility when deciding upon treatment intensification in high-risk patients. Genomic and clinical data are available as a public resource.

Prognostic and therapeutic relevance of FLIP and procaspase-8 overexpression in non-small cell lung cancer.

Non-small cell lung carcinoma remains by far the leading cause of cancer-related deaths worldwide. Overexpression of FLIP, which blocks the extrinsic apoptotic pathway by inhibiting caspase-8 activation, has been identified in various cancers. We investigated FLIP and procaspase-8 expression in NSCLC and the effect of HDAC inhibitors on FLIP expression, activation of caspase-8 and drug resistance in NSCLC and normal lung cell line models. Immunohistochemical analysis of cytoplasmic and nuclear FLIP and procaspase-8 protein expression was carried out using a novel digital pathology approach. Both FLIP and procaspase-8 were found to be significantly overexpressed in tumours, and importantly, high cytoplasmic expression of FLIP significantly correlated with shorter overall survival. Treatment with HDAC inhibitors targeting HDAC1-3 downregulated FLIP expression predominantly via post-transcriptional mechanisms, and this resulted in death receptor- and caspase-8-dependent apoptosis in NSCLC cells, but not normal lung cells. In addition, HDAC inhibitors synergized with TRAIL and cisplatin in NSCLC cells in a FLIP- and caspase-8-dependent manner. Thus, FLIP and procaspase-8 are overexpressed in NSCLC, and high cytoplasmic FLIP expression is indicative of poor prognosis. Targeting high FLIP expression using HDAC1-3 selective inhibitors such as entinostat to exploit high procaspase-8 expression in NSCLC has promising therapeutic potential, particularly when used in combination with TRAIL receptor-targeted agents.