Stefan Graw

Robust gene expression and mutation analyses of RNA-sequencing of formalin-fixed diagnostic tumor samples

Current genomic studies are limited by the availability of fresh tissue samples. Here, we show that Illumina RNA sequencing of formalin-fixed diagnostic tumor samples produces gene expression that is strongly correlated with matched frozen tumor samples (r > 0.89). In addition, sequence variations identified from FFPE RNA show 99.67% concordance with that from exome sequencing of matched frozen tumor samples. Because FFPE is a routine diagnostic sample preparation, the feasibility results reported here will facilitate the setup of large-scale research and clinical studies in medical genomics that are currently limited by the availability of fresh frozen samples.

Sustained O-GlcNAcylation reprograms mitochondrial function to regulate energy metabolism

Dysfunctional mitochondria and generation of reactive oxygen species (ROS) promote chronic diseases, which have spurred interest in the molecular mechanisms underlying these conditions. Previously, we have demonstrated that disruption of post-translational modification of proteins with β-linked N-acetylglucosamine (O-GlcNAcylation) via overexpression of the O-GlcNAc-regulating enzymes O-GlcNAc transferase (OGT) or O-GlcNAcase (OGA) impairs mitochondrial function. Here, we report that sustained alterations in O-GlcNAcylation either by pharmacological or genetic manipulation also alter metabolic function. Sustained O-GlcNAc elevation in SH-SY5Y neuroblastoma cells increased OGA expression and reduced cellular respiration and ROS generation. Cells with elevated O-GlcNAc levels had elongated mitochondria and increased mitochondrial membrane potential, and RNA-sequencing analysis indicated transcriptome reprogramming and down-regulation of the NRF2-mediated antioxidant response. Sustained O-GlcNAcylation in mouse brain and liver validated the metabolic phenotypes observed in the cells, and OGT knockdown in the liver elevated ROS levels, impaired respiration, and increased the NRF2 antioxidant response. Moreover, elevated O-GlcNAc levels promoted weight loss and lowered respiration in mice and skewed the mice toward carbohydrate-dependent metabolism as determined by indirect calorimetry. In summary, sustained elevation in O-GlcNAcylation coupled with increased OGA expression reprograms energy metabolism, a finding that has potential implications for the etiology, development, and management of metabolic diseases.

Mitochondrial genomic backgrounds affect nuclear DNA methylation and gene expression

Mitochondrial DNA (mtDNA) mutations and polymorphisms contribute to many complex diseases, including cancer. Using a unique mouse model that contains nDNA from one mouse strain and homoplasmic mitochondrial haplotypes from different mouse strain(s)-designated Mitochondrial Nuclear Exchange (MNX)-we showed that mtDNA could alter mammary tumor metastasis. Because retrograde and anterograde communication exists between the nuclear and mitochondrial genomes, we hypothesized that there are differential mtDNA-driven changes in nuclear (n)DNA expression and DNA methylation. Genome-wide nDNA methylation and gene expression were measured in harvested brain tissue from paired wild-type and MNX mice. Selective differential DNA methylation and gene expression were observed between strains having identical nDNA, but different mtDNA. These observations provide insights into how mtDNA could be altering epigenetic regulation and thereby contribute to the pathogenesis of metastasis. Cancer Res; 77(22); 6202-14. ©2017 AACR.

O-GlcNAc Homeostasis Controls Cell Fate Decisions During Hematopoiesis

The addition of O-GlcNAc (a single β-D-N-acetylglucosamine sugar at serine and threonine residues) by O-GlcNAc transferase (OGT) and removal by O-GlcNAcase (OGA) maintains homeostatic levels of O-GlcNAc. We investigated the role of O-GlcNAc homeostasis in hematopoiesis utilizing G1E-ER4 cells carrying a GATA-1 transcription factor fused to the estrogen receptor (GATA-1ER) that undergo erythropoiesis following the addition of β-estradiol (E2) and myeloid leukemia cells that differentiate into neutrophils in the presence of all-trans retinoic acid. During G1E-ER4 differentiation, a decrease in overall O-GlcNAc levels and an increase in GATA-1 interactions with OGT and OGA were observed. Transcriptome analysis on G1E-ER4 cells differentiated in the presence of Thiamet-G (TMG), an OGA inhibitor, identified expression changes in 433 GATA-1 target genes. Chromatin immunoprecipitation demonstrated that the occupancy of GATA-1, OGT, and OGA at Laptm5 gene GATA site was decreased with TMG. Myeloid leukemia cells showed a decline in O-GlcNAc levels during differentiation and TMG reduced the expression of genes involved in differentiation. Sustained treatment with TMG in G1E-ER4 cells prior to differentiation caused a reduction of hemoglobin positive cells during differentiation. Our results show that alterations in O-GlcNAc homeostasis disrupt transcriptional programs causing differentiation errors suggesting a vital role of O-GlcNAcylation in control of cell fate.

Licofelone enhances the efficacy of paclitaxel in ovarian cancer by reversing drug resistance and tumor stem-like properties

Drug development for first-line treatment of epithelial ovarian cancer (EOC) has been stagnant for almost three decades. Traditional cell culture methods for primary drug screening do not always accurately reflect clinical disease. To overcome this barrier, we grew a panel of EOC cell lines in three-dimensional (3D) cell cultures to form multicellular tumor spheroids (MCTS). We characterized these MCTS for molecular and cellular features of EOC and performed a comparative screen with cells grown using two-dimensional (2D) cell culture to identify previously unappreciated anticancer drugs. MCTS exhibited greater resistance to chemotherapeutic agents, showed signs of senescence and hypoxia, and expressed a number of stem cell-associated transcripts including ALDH1A and CD133, also known as PROM1 Using a library of clinically repurposed drugs, we identified candidates with preferential activity in MCTS over 2D cultured cells. One of the lead compounds, the dual COX/LOX inhibitor licofelone, reversed the stem-like properties of ovarian MCTS. Licofelone also synergized with paclitaxel in ovarian MCTS models and in a patient-derived tumor xenograft model. Importantly, the combination of licofelone with paclitaxel prolonged the median survival of mice (>141 days) relative to paclitaxel (115 days), licofelone (37 days), or vehicle (30 days). Increased efficacy was confirmed by Mantel-Haenszel HR compared with vehicle (HR = 0.037) and paclitaxel (HR = 0.017). These results identify for the first time an unappreciated, anti-inflammatory drug that can reverse chemotherapeutic resistance in ovarian cancer, highlighting the need to clinically evaluate licofelone in combination with first-line chemotherapy in primary and chemotherapy-refractory EOC.Significance: This study highlights the use of an in vitro spheroid 3D drug screening model to identify new therapeutic approaches to reverse chemotherapy resistance in ovarian cancer. Cancer Res; 78(15); 4370-85. ©2018 AACR.

Digit-a tool for detection and identification of genomic interchromosomal translocations

Abstract Structural variations (SVs) in genomic DNA can have profound effects on the evolution of living organisms, on phenotypic variations and on disease processes. A critical step in discovering the full extent of structural variations is the development of tools to characterize these variations accurately in next generation sequencing data. Toward this goal, we developed a software pipeline named digit that implements a novel measure of mapping ambiguity to discover interchromosomal SVs from mate-pair and pair-end sequencing data. The workflow robustly handles the high numbers of artifacts present in mate-pair sequencing and reduces the false positive rate while maintaining sensitivity. In the simulated data set, our workflow recovered 96% of simulated SVs. It generates a self-updating library of common translocations and allows for the investigation of patient- or group-specific events, making it suitable for discovering and cataloging chromosomal translocations associated with specific groups, traits, diseases or population structures.

An ensemble-based Cox proportional hazards regression framework for predicting survival in metastatic castration-resistant prostate cancer (mCRPC) patients

From March through August 2015, nearly 60 teams from around the world participated in the Prostate Cancer Dream Challenge (PCDC). Participating teams were faced with the task of developing prediction models for patient survival and treatment discontinuation using baseline clinical variables collected on metastatic castrate-resistant prostate cancer (mCRPC) patients in the comparator arm of four phase III clinical trials. In total, over 2,000 mCRPC patients treated with first-line docetaxel comprised the training and testing data sets used in this challenge. In this paper we describe: (a) the sub-challenges comprising the PCDC, (b) the statistical metrics used to benchmark prediction performance, (c) our analytical approach, and finally (d) our team’s overall performance in this challenge. Specifically, we discuss our curated, ad-hoc, feature selection (CAFS) strategy for identifying clinically important risk-predictors, the ensemble-based Cox proportional hazards regression framework used in our final submission, and the adaptation of our modeling framework based on the results from the intermittent leaderboard rounds. Strong predictors of patient survival were successfully identified utilizing our model building approach. Several of the identified predictors were new features created by our team via strategically merging collections of weak predictors. In each of the three intermittent leaderboard rounds, our prediction models scored among the top four models across all participating teams and our final submission ranked 9 th place overall with an integrated area under the curve (iAUC) of 0.7711 computed in an independent test set. While the prediction performance of teams placing between 2 nd- 10 th (iAUC: 0.7710-0.7789) was better than the current gold-standard prediction model for prostate cancer survival, the top-performing team, FIMM-UTU significantly outperformed all other contestants with an iAUC of 0.7915. In summary, our ensemble-based Cox regression framework with CAFS resulted in strong overall performance for predicting prostate cancer survival and represents a promising approach for future prediction problems.

Abstract POSTER-TECH-1109: Robust gene expression and mutation analyses from RNA-sequencing of formalin-fixed diagnostic tumor samples

Abstracts: 10th Biennial Ovarian Cancer Research Symposium; September 8-9, 2014; Seattle, WA Formalin-Fixed Paraffin Embedding (FFPE) has been the standard sample preparation for pathologists for decades. This created a mass of disease and normal tissue biospecimens with rich clinical annotation and patient follow-up data. If these FFPE archives proved useful for the next-generation sequence analyses, FFPE samples may be sequenced to investigate the complex genetic changes underlying tumor progression, therapy resistance and variability in disease outcome. Previous studies have successfully used FFPE-DNA for copy number analysis and mutation detection even though it is technically difficult to perform sequencing analyses on DNA isolated from FFPE samples. It is even more challenging to perform RNA sequence analysis because routine FFPE processing does not preserve high quality RNA. Aim of this study is to determine the extent to which FFPE RNA expression is correlated with matched frozen tumor samples. Six pairs of matched fresh frozen (FF) and FFPE ovarian tumor samples were processed for RNA sequence analysis. Multiple quality checks were performed before and after aligning the reads to the reference genome (hg19/GRCh37). Differential gene expression analysis was performed utilizing DESeq on aligned RNA-sequencing data, as well as data obtained by NanoString technology (250 selected genes). Both results were compared to each other. A pipeline called SNPiR was applied to detect mutation from the RNA sequencing data, and the mutation calls between the FF and FFPE samples were compared at the RNA and DNA level. For each analysis FF samples from the same cancer patients were used to verify the results from FFPE sample. The quality checks indicate that FFPE samples show uniform coverage of 5’ and 3’ of the transcripts whereas FF samples show 3’ bias in transcript coverage. This difference is likely attributable to differences in library preparation steps for FF and FFPE samples. The sequencing libraries for FF samples were done by poly-A selection while the libraries for FFPE samples were prepared using Ribo- Minus RNA sequencing kit. The differential gene expression analysis showed a good correlation between each matching FF and FFPE samples with one exception. Analyzing the NanoString gene counts show not only a remarkable correlation between the matching FF and FFPE sample but also a correlation between FF and FFPE samples between the two techniques (RNA-sequencing and NanoString). Unfortunately the mutation calling show variable results with 3 paired samples showing high rates of concordant mutation calls between FF and FFPE samples while the remaining 3 paired samples showing low rates of concordant mutation calls. Ongoing studies are focusing on biological and technical variability that may account for the observed differences in concordant calls. In summary, the results indicate that FFPE samples are a suitable replacement for FF samples in differential gene expression analysis. However, mutation analysis from FFPE samples is challenging and results are variable. Citation Format: Stefan Graw, Richard Meier, Kay Minn, Andrew K Godwin, Peter Beyerlein, Jeremy Chien. Robust gene expression and mutation analyses from RNA-sequencing of formalin-fixed diagnostic tumor samples [abstract]. In: Proceedings of the 10th Biennial Ovarian Cancer Research Symposium; Sep 8-9, 2014; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(16 Suppl):Abstract nr POSTER-TECH-1109.