Publications

2014

Moskwa P, Zinn PO, Choi YE, Shukla SA, Fendler W, Chen CC, Lu J, Golub TR, Hjelmeland A, Chowdhury D. A functional screen identifies miRs that induce radioresistance in glioblastomas. Mol Cancer Res. 2014;12:1767–78.
UNLABELLED: The efficacy of radiotherapy in many tumor types is limited by normal tissue toxicity and by intrinsic or acquired radioresistance. Therefore, it is essential to understand the molecular network responsible for regulating radiosensitivity/resistance. Here, an unbiased functional screen identified four microRNAs (miR1, miR125a, miR150, and miR425) that induce radioresistance. Considering the clinical importance of radiotherapy for patients with glioblastoma, the impact of these miRNAs on glioblastoma radioresistance was investigated. Overexpression of miR1, miR125a, miR150, and/or miR425 in glioblastoma promotes radioresistance through upregulation of the cell-cycle checkpoint response. Conversely, antagonizing with antagomiRs sensitizes glioblastoma cells to irradiation, suggesting their potential as targets for inhibiting therapeutic resistance. Analysis of glioblastoma datasets from The Cancer Genome Atlas (TCGA) revealed that these miRNAs are expressed in glioblastoma patient specimens and correlate with TGFbeta signaling. Finally, it is demonstrated that expression of miR1 and miR125a can be induced by TGFbeta and antagonized by a TGFbeta receptor inhibitor. Together, these results identify and characterize a new role for miR425, miR1, miR125, and miR150 in promoting radioresistance in glioblastomas and provide insight into the therapeutic application of TGFbeta inhibitors in radiotherapy. IMPLICATIONS: Systematic identification of miRs that cause radioresistance in gliomas is important for uncovering predictive markers for radiotherapy or targets for overcoming radioresistance.

2013

Hartwell KA, Miller PG, Mukherjee S, Kahn AR, Stewart AL, Logan DJ, Negri JM, Duvet M, Jaras M, Puram R, et al. Niche-based screening identifies small-molecule inhibitors of leukemia stem cells. Nat Chem Biol. 2013;9:840–848.
Efforts to develop more effective therapies for acute leukemia may benefit from high-throughput screening systems that reflect the complex physiology of the disease, including leukemia stem cells (LSCs) and supportive interactions with the bone marrow microenvironment. The therapeutic targeting of LSCs is challenging because LSCs are highly similar to normal hematopoietic stem and progenitor cells (HSPCs) and are protected by stromal cells in vivo. We screened 14,718 compounds in a leukemia-stroma co-culture system for inhibition of cobblestone formation, a cellular behavior associated with stem-cell function. Among those compounds that inhibited malignant cells but spared HSPCs was the cholesterol-lowering drug lovastatin. Lovastatin showed anti-LSC activity in vitro and in an in vivo bone marrow transplantation model. Mechanistic studies demonstrated that the effect was on target, via inhibition of HMG-CoA reductase. These results illustrate the power of merging physiologically relevant models with high-throughput screening.
Baca SC, Prandi D, Lawrence MS, Mosquera JM, Romanel A, Drier Y, Park K, Kitabayashi N, MacDonald TY, Ghandi M, et al. Punctuated evolution of prostate cancer genomes. Cell. 2013;153:666–77.
The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.
Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, Carter SL, Stewart C, Mermel CH, Roberts SA, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499:214–218.
Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.
Hoshida Y, Villanueva A, Sangiovanni A, Sole M, Hur C, Andersson KL, Chung RT, Gould J, Kojima K, Gupta S, et al. Prognostic gene expression signature for patients with hepatitis C-related early-stage cirrhosis. Gastroenterology. 2013;144:1024–30.
BACKGROUND & AIMS: Cirrhosis affects 1% to 2% of the world population and is the major risk factor for hepatocellular carcinoma (HCC). Hepatitis C cirrhosis-related HCC is the most rapidly increasing cause of cancer death in the United States. Noninvasive methods have been developed to identify patients with asymptomatic early-stage cirrhosis, increasing the burden of HCC surveillance, but biomarkers are needed to identify patients with cirrhosis who are most in need of surveillance. We investigated whether a liver-derived 186-gene signature previously associated with outcomes of patients with HCC is prognostic for patients with newly diagnosed cirrhosis but without HCC. METHODS: We performed gene expression profile analysis of formalin-fixed needle biopsy specimens from the livers of 216 patients with hepatitis C-related early-stage (Child-Pugh class A) cirrhosis who were prospectively followed up for a median of 10 years at an Italian center. We evaluated whether the 186-gene signature was associated with death, progression of cirrhosis, and development of HCC. RESULTS: Fifty-five (25%), 101 (47%), and 60 (28%) patients were classified as having poor-, intermediate-, and good-prognosis signatures, respectively. In multivariable Cox regression modeling, the poor-prognosis signature was significantly associated with death (P = .004), progression to advanced cirrhosis (P .001), and development of HCC (P = .009). The 10-year rates of survival were 63%, 74%, and 85% and the annual incidence of HCC was 5.8%, 2.2%, and 1.5% for patients with poor-, intermediate-, and good-prognosis signatures, respectively. CONCLUSIONS: A 186-gene signature used to predict outcomes of patients with HCC is also associated with outcomes of patients with hepatitis C-related early-stage cirrhosis. This signature might be used to identify patients with cirrhosis in most need of surveillance and strategies to prevent the development of HCC.
Santagata S, Mendillo ML, Tang YC, Subramanian A, Perley CC, Roche SP, Wong B, Narayan R, Kwon H, Koeva M, et al. Tight coordination of protein translation and HSF1 activation supports the anabolic malignant state. Science. 2013;341:1238303.
The ribosome is centrally situated to sense metabolic states, but whether its activity, in turn, coherently rewires transcriptional responses is unknown. Here, through integrated chemical-genetic analyses, we found that a dominant transcriptional effect of blocking protein translation in cancer cells was inactivation of heat shock factor 1 (HSF1), a multifaceted transcriptional regulator of the heat-shock response and many other cellular processes essential for anabolic metabolism, cellular proliferation, and tumorigenesis. These analyses linked translational flux to the regulation of HSF1 transcriptional activity and to the modulation of energy metabolism. Targeting this link with translation initiation inhibitors such as rocaglates deprived cancer cells of their energy and chaperone armamentarium and selectively impaired the proliferation of both malignant and premalignant cells with early-stage oncogenic lesions.
Dulak AM, Stojanov P, Peng S, Lawrence MS, Fox C, Stewart C, Bandla S, Imamura Y, Schumacher SE, Shefler E, et al. Exome and whole-genome sequencing of esophageal adenocarcinoma identifies recurrent driver events and mutational complexity. Nat Genet. 2013;45:478–86.
The incidence of esophageal adenocarcinoma (EAC) has risen 600% over the last 30 years. With a 5-year survival rate of ~15%, the identification of new therapeutic targets for EAC is greatly important. We analyze the mutation spectra from whole-exome sequencing of 149 EAC tumor-normal pairs, 15 of which have also been subjected to whole-genome sequencing. We identify a mutational signature defined by a high prevalence of A>C transversions at AA dinucleotides. Statistical analysis of exome data identified 26 significantly mutated genes. Of these genes, five (TP53, CDKN2A, SMAD4, ARID1A and PIK3CA) have previously been implicated in EAC. The new significantly mutated genes include chromatin-modifying factors and candidate contributors SPG20, TLR4, ELMO1 and DOCK2. Functional analyses of EAC-derived mutations in ELMO1 identifies increased cellular invasion. Therefore, we suggest the potential activation of the RAC1 pathway as a contributor to EAC tumorigenesis.
Gustafsdottir SM, Ljosa V, Sokolnicki KL, Anthony Wilson J, Walpita D, Kemp MM, Petri Seiler K, Carrel HA, Golub TR, Schreiber SL, et al. Multiplex cytological profiling assay to measure diverse cellular states. PLoS One. 2013;8:e80999.
Computational methods for image-based profiling are under active development, but their success hinges on assays that can capture a wide range of phenotypes. We have developed a multiplex cytological profiling assay that "paints the cell" with as many fluorescent markers as possible without compromising our ability to extract rich, quantitative profiles in high throughput. The assay detects seven major cellular components. In a pilot screen of bioactive compounds, the assay detected a range of cellular phenotypes and it clustered compounds with similar annotated protein targets or chemical structure based on cytological profiles. The results demonstrate that the assay captures subtle patterns in the combination of morphological labels, thereby detecting the effects of chemical compounds even though their targets are not stained directly. This image-based assay provides an unbiased approach to characterize compound- and disease-associated cell states to support future probe discovery.

2012

Straussman R, Morikawa T, Shee K, Barzily-Rokni M, Qian ZR, Du J, Davis A, Mongare MM, Gould J, Frederick DT, et al. Tumour micro-environment elicits innate resistance to RAF inhibitors through HGF secretion. Nature. 2012;487:500–4.
Drug resistance presents a challenge to the treatment of cancer patients. Many studies have focused on cell-autonomous mechanisms of drug resistance. By contrast, we proposed that the tumour micro-environment confers innate resistance to therapy. Here we developed a co-culture system to systematically assay the ability of 23 stromal cell types to influence the innate resistance of 45 cancer cell lines to 35 anticancer drugs. We found that stroma-mediated resistance is common, particularly to targeted agents. We characterized further the stroma-mediated resistance of BRAF-mutant melanoma to RAF inhibitors because most patients with this type of cancer show some degree of innate resistance. Proteomic analysis showed that stromal cell secretion of hepatocyte growth factor (HGF) resulted in activation of the HGF receptor MET, reactivation of the mitogen-activated protein kinase (MAPK) and phosphatidylinositol-3-OH kinase (PI(3)K)-AKT signalling pathways, and immediate resistance to RAF inhibition. Immunohistochemistry experiments confirmed stromal cell expression of HGF in patients with BRAF-mutant melanoma and showed a significant correlation between HGF expression by stromal cells and innate resistance to RAF inhibitor treatment. Dual inhibition of RAF and either HGF or MET resulted in reversal of drug resistance, suggesting RAF plus HGF or MET inhibitory combination therapy as a potential therapeutic strategy for BRAF-mutant melanoma. A similar resistance mechanism was uncovered in a subset of BRAF-mutant colorectal and glioblastoma cell lines. More generally, this study indicates that the systematic dissection of interactions between tumours and their micro-environment can uncover important mechanisms underlying drug resistance.
Wei G, Margolin AA, Haery L, Brown E, Cucolo L, Julian B, Shehata S, Kung AL, Beroukhim R, Golub TR. Chemical genomics identifies small-molecule MCL1 repressors and BCL-xL as a predictor of MCL1 dependency. Cancer Cell. 2012;21:547–62.
MCL1, which encodes the antiapoptotic protein MCL1, is among the most frequently amplified genes in human cancer. A chemical genomic screen identified compounds, including anthracyclines, that decreased MCL1 expression. Genomic profiling indicated that these compounds were global transcriptional repressors that preferentially affect MCL1 due to its short mRNA half-life. Transcriptional repressors and MCL1 shRNAs induced apoptosis in the same cancer cell lines and could be rescued by physiological levels of ectopic MCL1 expression. Repression of MCL1 released the proapoptotic protein BAK from MCL1, and Bak deficiency conferred resistance to transcriptional repressors. A computational model, validated in vivo, indicated that high BCL-xL expression confers resistance to MCL1 repression, thereby identifying a patient-selection strategy for the clinical development of MCL1 inhibitors.