Publications

2009

Olson P, Lu J, Zhang H, Shai A, Chun MG, Wang Y, Libutti SK, Nakakura EK, Golub TR, Hanahan D. MicroRNA dynamics in the stages of tumorigenesis correlate with hallmark capabilities of cancer. Genes Dev. 2009;23:2152–65.
While altered expression of microRNAs (miRs) in tumors has been well documented, it remains unclear how the miR transcriptome intersects neoplastic progression. By profiling the miR transcriptome we identified miR expression signatures associated with steps in tumorigenesis and the acquisition of hallmark capabilities in a prototypical mouse model of cancer. Metastases and a rare subset of primary tumors shared a distinct miR signature, implicating a discrete lineage for metastatic tumors. The miR-200 family is strongly down-regulated in metastases and met-like primary tumors, thereby relieving repression of the mesenchymal transcription factor Zeb1, which in turn suppresses E-cadherin. Treatment with a clinically approved angiogenesis inhibitor normalized angiogenic signature miRs in primary tumors, while altering expression of metastatic signature miRs similarly to liver metastases, suggesting their involvement in adaptive resistance to anti-angiogenic therapy via enhanced metastasis. Many of the miR changes associated with specific stages and hallmark capabilities in the mouse model are similarly altered in human tumors, including cognate pancreatic neuroendocrine tumors, implying a generality.
Ong SE, Schenone M, Margolin AA, Li X, Do K, Doud MK, Mani DR, Kuai L, Wang X, Wood JL, et al. Identifying the proteins to which small-molecule probes and drugs bind in cells. Proc Natl Acad Sci U S A. 2009;106:4617–22.
Most small-molecule probes and drugs alter cell circuitry by interacting with 1 or more proteins. A complete understanding of the interacting proteins and their associated protein complexes, whether the compounds are discovered by cell-based phenotypic or target-based screens, is extremely rare. Such a capability is expected to be highly illuminating--providing strong clues to the mechanisms used by small-molecules to achieve their recognized actions and suggesting potential unrecognized actions. We describe a powerful method combining quantitative proteomics (SILAC) with affinity enrichment to provide unbiased, robust and comprehensive identification of the proteins that bind to small-molecule probes and drugs. The method is scalable and general, requiring little optimization across different compound classes, and has already had a transformative effect on our studies of small-molecule probes. Here, we describe in full detail the application of the method to identify targets of kinase inhibitors and immunophilin binders.

2008

Takeyama K, Monti S, Manis JP, Dal Cin P, Getz G, Beroukhim R, Dutt S, Aster JC, Alt FW, Golub TR, et al. Integrative analysis reveals 53BP1 copy loss and decreased expression in a subset of human diffuse large B-cell lymphomas. Oncogene. 2008;27:318–22.
p53-Binding protein 1 (53BP1) encodes a critical checkpoint protein that localizes to sites of DNA double-strand breaks (DSBs) and participates in DSB repair. Mice that are 53bp1 deficient or hemizygous have an increased incidence of lymphoid malignancies. However, 53BP1 abnormalities in primary human tumors have not been described. By combining high-density single nucleotide polymorphism (HD SNP) array data and gene expression profiles, we found 9 of 63 newly diagnosed human diffuse large B-cell lymphomas (DLBCLs) with single copy loss of the chromosome 15q15 region including the 53BP1 locus; these nine tumors also had significantly lower levels of 53BP1 transcripts. 53BP1 single copy loss found with the HD SNP array platform was subsequently confirmed by fluorescence in situ hybridization. These studies highlight the role of 53BP1 copy loss in primary human DLBCLs and the value of integrative analyses in detecting this genetic lesion in human tumors.
Tam WF, Gu TL, Chen J, Lee BH, Bullinger L, Frohling S, Wang A, Monti S, Golub TR, Gilliland DG. Id1 is a common downstream target of oncogenic tyrosine kinases in leukemic cells. Blood. 2008;112:1981–92.
Oncogenic tyrosine kinases, such as BCR-ABL, TEL-ABL, TEL-PDGFbetaR, and FLT3-ITD, play a major role in the development of hematopoietic malignancy. They activate many of the same signal transduction pathways. To identify the critical target genes required for transformation in hematopoietic cells, we used a comparative gene expression strategy in which selective small molecules were applied to 32Dcl3 cells that had been transformed to factor-independent growth by these respective oncogenic alleles. We identified inhibitor of DNA binding 1 (Id1), a gene involved in development, cell cycle, and tumorigenesis, as a common target of these oncogenic kinases. These findings were prospectively confirmed in cell lines and primary bone marrow cells engineered to express the respective tyrosine kinase alleles and were also confirmed in vivo in murine models of disease. Moreover, human AML cell lines Molm-14 and K562, which express the FLT3-ITD and BCR-ABL tyrosine kinases, respectively, showed high levels of Id1 expression. Antisense and siRNA based knockdown of Id1-inhibited growth of these cells associated with increased p27(Kip1) expression and increased sensitivity to Trail-induced apoptosis. These findings indicate that Id1 is an important target of constitutively activated tyrosine kinases and may be a therapeutic target for leukemias associated with oncogenic tyrosine kinases.
Wagner BK, Kitami T, Gilbert TJ, Peck D, Ramanathan A, Schreiber SL, Golub TR, Mootha VK. Large-scale chemical dissection of mitochondrial function. Nat Biotechnol. 2008;26:343–51.
Mitochondrial oxidative phosphorylation (OXPHOS) is under the control of both mitochondrial (mtDNA) and nuclear genomes and is central to energy homeostasis. To investigate how its function and regulation are integrated within cells, we systematically combined four cell-based assays of OXPHOS physiology with multiplexed measurements of nuclear and mtDNA gene expression across 2,490 small-molecule perturbations in cultured muscle. Mining the resulting compendium revealed, first, that protein synthesis inhibitors can decouple coordination of nuclear and mtDNA transcription; second, that a subset of HMG-CoA reductase inhibitors, combined with propranolol, can cause mitochondrial toxicity, yielding potential clues about the etiology of statin myopathy; and, third, that structurally diverse microtubule inhibitors stimulate OXPHOS transcription while suppressing reactive oxygen species, via a transcriptional mechanism involving PGC-1alpha and ERRalpha, and thus may be useful in treating age-associated degenerative disorders. Our screening compendium can be used as a discovery tool both for understanding mitochondrial biology and toxicity and for identifying novel therapeutics.
Antipova AA, Stockwell BR, Golub TR. Gene expression-based screening for inhibitors of PDGFR signaling. Genome Biol. 2008;9:R47.
Here we describe a proof-of-concept experiment designed to explore the possibility of using gene expression-based high-throughput screening (GE-HTS) to find inhibitors of a signaling cascade, using platelet derived growth factor receptor (PDGFR) signaling as the example. The previously unrecognized ability of aurintricarboxylic acid to inhibit PDGFR signaling, discovered through a screen of 1,739 compounds, demonstrates the feasibility and generalizability of GE-HTS for the discovery of small molecule modulators of any signaling pathway of interest.
Ebert BL, Galili N, Tamayo P, Bosco J, Mak R, Pretz J, Tanguturi S, Ladd-Acosta C, Stone R, Golub TR, et al. An erythroid differentiation signature predicts response to lenalidomide in myelodysplastic syndrome. PLoS Med. 2008;5:e35.
BACKGROUND: Lenalidomide is an effective new agent for the treatment of patients with myelodysplastic syndrome (MDS), an acquired hematopoietic disorder characterized by ineffective blood cell production and a predisposition to the development of leukemia. Patients with an interstitial deletion of Chromosome 5q have a high rate of response to lenalidomide, but most MDS patients lack this deletion. Approximately 25% of patients without 5q deletions also benefit from lenalidomide therapy, but response in these patients cannot be predicted by any currently available diagnostic assays. The aim of this study was to develop a method to predict lenalidomide response in order to avoid unnecessary toxicity in patients unlikely to benefit from treatment. METHODS AND FINDINGS: Using gene expression profiling, we identified a molecular signature that predicts lenalidomide response. The signature was defined in a set of 16 pretreatment bone marrow aspirates from MDS patients without 5q deletions, and validated in an independent set of 26 samples. The response signature consisted of a cohesive set of erythroid-specific genes with decreased expression in responders, suggesting that a defect in erythroid differentiation underlies lenalidomide response. Consistent with this observation, treatment with lenalidomide promoted erythroid differentiation of primary hematopoietic progenitor cells grown in vitro. CONCLUSIONS: These studies indicate that lenalidomide-responsive patients have a defect in erythroid differentiation, and suggest a strategy for a clinical test to predict patients most likely to respond to the drug. The experiments further suggest that the efficacy of lenalidomide, whose mechanism of action in MDS is unknown, may be due to its ability to induce erythroid differentiation.
Ebert BL, Pretz J, Bosco J, Chang CY, Tamayo P, Galili N, Raza A, Root DE, Attar E, Ellis SR, et al. Identification of RPS14 as a 5q- syndrome gene by RNA interference screen. Nature. 2008;451:335–9.
Somatic chromosomal deletions in cancer are thought to indicate the location of tumour suppressor genes, by which a complete loss of gene function occurs through biallelic deletion, point mutation or epigenetic silencing, thus fulfilling Knudson's two-hit hypothesis. In many recurrent deletions, however, such biallelic inactivation has not been found. One prominent example is the 5q- syndrome, a subtype of myelodysplastic syndrome characterized by a defect in erythroid differentiation. Here we describe an RNA-mediated interference (RNAi)-based approach to discovery of the 5q- disease gene. We found that partial loss of function of the ribosomal subunit protein RPS14 phenocopies the disease in normal haematopoietic progenitor cells, and also that forced expression of RPS14 rescues the disease phenotype in patient-derived bone marrow cells. In addition, we identified a block in the processing of pre-ribosomal RNA in RPS14-deficient cells that is functionally equivalent to the defect in Diamond-Blackfan anaemia, linking the molecular pathophysiology of the 5q- syndrome to a congenital syndrome causing bone marrow failure. These results indicate that the 5q- syndrome is caused by a defect in ribosomal protein function and suggest that RNAi screening is an effective strategy for identifying causal haploinsufficiency disease genes.
Hahn CK, Ross KN, Warrington IM, Mazitschek R, Kanegai CM, Wright RD, Kung AL, Golub TR, Stegmaier K. Expression-based screening identifies the combination of histone deacetylase inhibitors and retinoids for neuroblastoma differentiation. Proc Natl Acad Sci U S A. 2008;105:9751–6.
The discovery of new small molecules and their testing in rational combination poses an ongoing problem for rare diseases, in particular, for pediatric cancers such as neuroblastoma. Despite maximal cytotoxic therapy with double autologous stem cell transplantation, outcome remains poor for children with high-stage disease. Because differentiation is aberrant in this malignancy, compounds that modulate transcription, such as histone deacetylase (HDAC) inhibitors, are of particular interest. However, as single agents, HDAC inhibitors have had limited efficacy. In the present study, we use an HDAC inhibitor as an enhancer to screen a small-molecule library for compounds inducing neuroblastoma maturation. To quantify differentiation, we use an enabling gene expression-based screening strategy. The top hit identified in the screen was all-trans-retinoic acid. Secondary assays confirmed greater neuroblastoma differentiation with the combination of an HDAC inhibitor and a retinoid versus either alone. Furthermore, effects of combination therapy were synergistic with respect to inhibition of cellular viability and induction of apoptosis. In a xenograft model of neuroblastoma, animals treated with combination therapy had the longest survival. This work suggests that testing of an HDAC inhibitor and retinoid in combination is warranted for children with neuroblastoma and demonstrates the success of a signature-based screening approach to prioritize compound combinations for testing in rare diseases.
Haining WN, Ebert BL, Subrmanian A, Wherry EJ, Eichbaum Q, Evans JW, Mak R, Rivoli S, Pretz J, Angelosanto J, et al. Identification of an evolutionarily conserved transcriptional signature of CD8 memory differentiation that is shared by T and B cells. J Immunol. 2008;181:1859–68.
After Ag encounter, naive lymphocytes differentiate into populations of memory cells that share a common set of functions including faster response to Ag re-exposure and the ability to self-renew. However, memory lymphocytes in different lymphocyte lineages are functionally and phenotypically diverse. It is not known whether discrete populations of T and B cells use similar transcriptional programs during differentiation into the memory state. We used cross-species genomic analysis to examine the pattern of genes up-regulated during the differentiation of naive lymphocytes into memory cells in multiple populations of human CD4, CD8, and B cell lymphocytes as well as two mouse models of memory development. We identified and validated a signature of genes that was up-regulated in memory cells compared with naive cells in both human and mouse CD8 memory differentiation, suggesting marked evolutionary conservation of this transcriptional program. Surprisingly, this conserved CD8 differentiation signature was also up-regulated during memory differentiation in CD4 and B cell lineages. To validate the biologic significance of this signature, we showed that alterations in this signature of genes could distinguish between functional and exhausted CD8 T cells from a mouse model of chronic viral infection. Finally, we generated genome-wide microarray data from tetramer-sorted human T cells and showed profound differences in this differentiation signature between T cells specific for HIV and those specific for influenza. Thus, our data suggest that in addition to lineage-specific differentiation programs, T and B lymphocytes use a common transcriptional program during memory development that is disrupted in chronic viral infection.