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.

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.