Matthew J. Brauer

A comprehensive transcriptional portrait of human cancer cell lines

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

Slow Growth Induces Heat-Shock Resistance in Normal and Respiratory-deficient Yeast

Yeast cells respond to a variety of environmental stresses, including heat shock and growth limitation. There is considerable overlap in these responses both from the point of view of gene expression patterns and cross-protection for survival. We performed experiments in which cells growing at different steady-state growth rates in chemostats were subjected to a short heat pulse. Gene expression patterns allowed us to partition genes whose expression responds to heat shock into subsets of genes that also respond to slow growth rate and those that do not. We found also that the degree of induction and repression of genes that respond to stress is generally weaker in respiratory deficient mutants, suggesting a role for increased respiratory activity in the apparent stress response to slow growth. Consistent with our gene expression results in wild-type cells, we found that cells growing more slowly are cross-protected for heat shock, i.e., better able to survive a lethal heat challenge. Surprisingly, however, we found no difference in cross-protection between respiratory-deficient and wild-type cells, suggesting induction of heat resistance at low growth rates is independent of respiratory activity, even though many of the changes in gene expression are not.

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality.

The programs GMAP and GSNAP, for aligning RNA-Seq and DNA-Seq datasets to genomes, have evolved along with advances in biological methodology to handle longer reads, larger volumes of data, and new types of biological assays. The genomic representation has been improved to include linear genomes that can compare sequences using single-instruction multiple-data (SIMD) instructions, compressed genomic hash tables with fast access using SIMD instructions, handling of large genomes with more than four billion bp, and enhanced suffix arrays (ESAs) with novel data structures for fast access. Improvements to the algorithms have included a greedy match-and-extend algorithm using suffix arrays, segment chaining using genomic hash tables, diagonalization using segmental hash tables, and nucleotide-level dynamic programming procedures that use SIMD instructions and eliminate the need for F-loop calculations. Enhancements to the functionality of the programs include standardization of indel positions, handling of ambiguous splicing, clipping and merging of overlapping paired-end reads, and alignments to circular chromosomes and alternate scaffolds. The programs have been adapted for use in pipelines by integrating their usage into R/Bioconductor packages such as gmapR and HTSeqGenie, and these pipelines have facilitated the discovery of numerous biological phenomena.

The STARD9/Kif16a Kinesin Associates with Mitotic Microtubules and Regulates Spindle Pole Assembly

During cell division, cells form the microtubule-based mitotic spindle, a highly specialized and dynamic structure that mediates proper chromosome transmission to daughter cells. Cancer cells can show perturbed mitotic spindles and an approach in cancer treatment has been to trigger cell killing by targeting microtubule dynamics or spindle assembly. To identify and characterize proteins necessary for spindle assembly, and potential antimitotic targets, we performed a proteomic and genetic analysis of 592 mitotic microtubule copurifying proteins (MMCPs). Screening for regulators that affect both mitosis and apoptosis, we report the identification and characterization of STARD9, a kinesin-3 family member, which localizes to centrosomes and stabilizes the pericentriolar material (PCM). STARD9-depleted cells have fragmented PCM, form multipolar spindles, activate the spindle assembly checkpoint (SAC), arrest in mitosis, and undergo apoptosis. Interestingly, STARD9-depletion synergizes with the chemotherapeutic agent taxol to increase mitotic death, demonstrating that STARD9 is a mitotic kinesin and a potential antimitotic target.

Sustained Brown Fat Stimulation and Insulin Sensitization by a Humanized Bispecific Antibody Agonist for Fibroblast Growth Factor Receptor 1/βKlotho Complex

Dissipating excess calories as heat through therapeutic stimulation of brown adipose tissues (BAT) has been proposed as a potential treatment for obesity-linked disorders. Here, we describe the generation of a humanized effector-less bispecific antibody that activates fibroblast growth factor receptor (FGFR) 1/βKlotho complex, a common receptor for FGF21 and FGF19. Using this molecule, we show that antibody-mediated activation of FGFR1/βKlotho complex in mice induces sustained energy expenditure in BAT, browning of white adipose tissue, weight loss, and improvements in obesity-associated metabolic derangements including insulin resistance, hyperglycemia, dyslipidemia and hepatosteatosis. In mice and cynomolgus monkeys, FGFR1/βKlotho activation increased serum high-molecular-weight adiponectin, which appears to contribute over time by enhancing the amplitude of the metabolic benefits. At the same time, insulin sensitization by FGFR1/βKlotho activation occurs even before the onset of weight loss in a manner that is independent of adiponectin. Together, selective activation of FGFR1/βKlotho complex with a long acting therapeutic antibody represents an attractive approach for the treatment of type 2 diabetes and other obesity-linked disorders through enhanced energy expenditure, insulin sensitization and induction of high-molecular-weight adiponectin.

Identification of Mutant K-Ras-dependent Phenotypes Using a Panel of Isogenic Cell Lines

Background: Many studies analyzing K-Ras function rely on overexpression. Results: Knock-out and knock-in of endogenous mutant K-ras have modest effects on downstream signaling but strong effects on gene expression and transformation. Conclusion: Genomic manipulation allows physiological determination of WT and mutant K-Ras consequences. Significance: Gene expression patterns can be used to monitor inhibition of mutant K-Ras. To assess the consequences of endogenous mutant K-Ras, we analyzed the signaling and biological properties of a small panel of isogenic cell lines. These include the cancer cell lines DLD1, HCT116, and Hec1A, in which either the WT or mutant K-ras allele has been disrupted, and SW48 colorectal cancer cells and human mammary epithelial cells in which a single copy of mutant K-ras was introduced at its endogenous genomic locus. We find that single copy mutant K-Ras causes surprisingly modest activation of downstream signaling to ERK and Akt. In contrast, a negative feedback signaling loop to EGFR and N-Ras occurs in some, but not all, of these cell lines. Mutant K-Ras also had relatively minor effects on cell proliferation and cell migration but more dramatic effects on cell transformation as assessed by growth in soft agar. Surprisingly, knock-out of the wild type K-ras allele consistently increased growth in soft agar, suggesting tumor-suppressive properties of this gene under these conditions. Finally, we examined the effects of single copy mutant K-Ras on global gene expression. Although transcriptional programs triggered by mutant K-Ras were generally quite distinct in the different cell lines, there was a small number of genes that were consistently overexpressed, and these could be used to monitor K-Ras inhibition in a panel of human tumor cell lines. We conclude that there are conserved components of mutant K-Ras signaling and phenotypes but that many depend on cell context and environmental cues.

Identification and Analysis of In Vivo VEGF Downstream Markers Link VEGF Pathway Activity with Efficacy of Anti- VEGF Therapies

Purpose: The aim of this study was to identify conserved pharmacodynamic and potential predictive biomarkers of response to anti-VEGF therapy using gene expression profiling in preclinical tumor models and in patients. Experimental Design: Surrogate markers of VEGF inhibition [VEGF-dependent genes or VEGF-dependent vasculature (VDV)] were identified by profiling gene expression changes induced in response to VEGF blockade in preclinical tumor models and in human biopsies from patients treated with anti-VEGF monoclonal antibodies. The potential value of VDV genes as candidate predictive biomarkers was tested by correlating high or low VDV gene expression levels in pretreatment clinical samples with the subsequent clinical efficacy of bevacizumab (anti-VEGF)-containing therapy. Results: We show that VDV genes, including direct and more distal VEGF downstream endothelial targets, enable detection of VEGF signaling inhibition in mouse tumor models and human tumor biopsies. Retrospective analyses of clinical trial data indicate that patients with higher VDV expression in pretreatment tumor samples exhibited improved clinical outcome when treated with bevacizumab-containing therapies. Conclusions: In this work, we identified surrogate markers (VDV genes) for in vivo VEGF signaling in tumors and showed clinical data supporting a correlation between pretreatment VEGF bioactivity and the subsequent efficacy of anti-VEGF therapy. We propose that VDV genes are candidate biomarkers with the potential to aid the selection of novel indications as well as patients likely to respond to anti-VEGF therapy. The data presented here define a diagnostic biomarker hypothesis based on translational research that warrants further evaluation in additional retrospective and prospective trials. Clin Cancer Res; 19(13); 3681–92. ©2013 AACR.

Constitutive NOTCH3 Signaling Promotes the Growth of Basal Breast Cancers

Notch ligands signal through one of four receptors on neighboring cells to mediate cell-cell communication and control cell fate, proliferation, and survival. Although aberrant Notch activation has been implicated in numerous malignancies, including breast cancer, the importance of individual receptors in distinct breast cancer subtypes and the mechanisms of receptor activation remain unclear. Using a novel antibody to detect active NOTCH3, we report here that NOTCH3 signals constitutively in a panel of basal breast cancer cell lines and in more than one third of basal tumors. Selective inhibition of individual ligands revealed that this signal does not require canonical ligand induction. A NOTCH3 antagonist antibody inhibited growth of basal lines, whereas a NOTCH3 agonist antibody enhanced the transformed phenotype in vitro and in tumor xenografts. Transcriptomic analyses generated a Notch gene signature that included Notch pathway components, the oncogene c-Myc, and the mammary stem cell regulator Id4 This signature drove clustering of breast cancer cell lines and tumors into the common subtypes and correlated with the basal classification. Our results highlight an unexpected ligand-independent induction mechanism and suggest that constitutive NOTCH3 signaling can drive an oncogenic program in a subset of basal breast cancers. Cancer Res; 77(6); 1439-52. ©2017 AACR.

A high-content cellular senescence screen identifies candidate tumor suppressors, including EPHA3.

Activation of a cellular senescence program is a common response to prolonged oncogene activation or tumor suppressor loss, providing a physiological mechanism for tumor suppression in premalignant cells. The link between senescence and tumor suppression supports the hypothesis that a loss-of-function screen measuring bona fide senescence marker activation should identify candidate tumor suppressors. Using a high-content siRNA screening assay for cell morphology and proliferation measures, we identify 12 senescence-regulating kinases and determine their senescence marker signatures, including elevation of senescence-associated β-galactosidase, DNA damage and p53 or p16INK4a expression. Consistent with our hypothesis, SNP array CGH data supports loss of gene copy number of five senescence-suppressing genes across multiple tumor samples. One such candidate is the EPHA3 receptor tyrosine kinase, a gene commonly mutated in human cancer. We demonstrate that selected intracellular EPHA3 tumor-associated point mutations decrease receptor expression level and/or receptor tyrosine kinase (RTK) activity. Our study therefore describes a new strategy to mine for novel candidate tumor suppressors and provides compelling evidence that EPHA3 mutations may promote tumorigenesis only when key senescence-inducing pathways have been inactivated.

Development of a robust flow cytometry-based pharmacodynamic assay to detect phospho- protein signals for phosphatidylinositol 3-kinase inhibitors in multiple myeloma

BackgroundThe phosphatidylinositol 3-kinase (PI3K) pathway plays an important role in multiple myeloma (MM), a blood cancer associated with uncontrolled proliferation of bone marrow plasma cells. This study aimed to develop a robust clinical pharmacodynamic (PD) assay to measure the on-target PD effects of the selective PI3K inhibitor GDC-0941 in MM patients.MethodsWe conducted an in vitro drug wash-out study to evaluate the feasibility of biochemical approaches in measuring the phosphorylation of S6 ribosomal protein (S6), one of the commonly used PD markers for PI3K pathway inhibition. We then developed a 7-color phospho-specific flow cytometry assay, or phospho flow assay, to measure the phosphorylation state of intracellular S6 in bone marrow aspirate (BMA) and peripheral blood (PB). Integrated mean fluorescence intensity (iMFI) was used to calculate fold changes of phosphorylation. Assay sensitivity was evaluated by comparing phospho flow with Meso Scale Discovery (MSD) and immunohistochemistry (IHC) assays. Finally, a sample handling method was developed to maintain the integrity of phospho signal during sample shipping and storage to ensure clinical application.ResultsThe phospho flow assay provided single-cell PD monitoring of S6 phosphorylation in tumor and surrogate cells using fixed BMA and PB, assessing pathway modulation in response to GDC-0941 with sensitivity similar to that of MSD assay. The one-shot sample fixation and handling protocol herein demonstrated exceptional preservation of protein phosphorylation. In contrast, the IHC assay was less sensitive in terms of signal quantification while the biochemical approach (MSD) was less suitable to assess PD activities due to the undesirable impact associated with cell isolation on the protein phosphorylation in tumor cells.ConclusionsWe developed a robust PD biomarker assay for the clinical evaluation of PI3K inhibitors in MM, allowing one to decipher the PD response in a relevant cell population. To our knowledge, this is the first report of an easily implemented clinical PD assay that incorporates an unbiased one-shot sample handling protocol, all (staining)-in-one (tube) phospho flow staining protocol, and an integrated modified data analysis for PD monitoring of kinase inhibitors in relevant cell populations in BMA and PB. The methods described here ensure a real-time, reliable and reproducible PD readout, which can provide information for dose selection as well as help to identify optimal combinations of targeted agents in early clinical trials.