Publications

2005

Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Although genomewide RNA expression analysis has become a routine tool in biomedical research, extracting biological insight from such information remains a major challenge. Here, we describe a powerful analytical method called Gene Set Enrichment Analysis (GSEA) for interpreting gene expression data. The method derives its power by focusing on gene sets, that is, groups of genes that share common biological function, chromosomal location, or regulation. We demonstrate how GSEA yields insights into several cancer-related data sets, including leukemia and lung cancer. Notably, where single-gene analysis finds little similarity between two independent studies of patient survival in lung cancer, GSEA reveals many biological pathways in common. The GSEA method is embodied in a freely available software package, together with an initial database of 1,325 biologically defined gene sets.
Sweet-Cordero A, Mukherjee S, Subramanian A, You H, Roix JJ, Ladd-Acosta C, Mesirov J, Golub TR, Jacks T. An oncogenic KRAS2 expression signature identified by cross-species gene-expression analysis. Nat Genet. 2005;37:48–55.
Using advanced gene targeting methods, generating mouse models of cancer that accurately reproduce the genetic alterations present in human tumors is now relatively straightforward. The challenge is to determine to what extent such models faithfully mimic human disease with respect to the underlying molecular mechanisms that accompany tumor progression. Here we describe a method for comparing mouse models of cancer with human tumors using gene-expression profiling. We applied this method to the analysis of a model of Kras2-mediated lung cancer and found a good relationship to human lung adenocarcinoma, thereby validating the model. Furthermore, we found that whereas a gene-expression signature of KRAS2 activation was not identifiable when analyzing human tumors with known KRAS2 mutation status alone, integrating mouse and human data uncovered a gene-expression signature of KRAS2 mutation in human lung cancer. We confirmed the importance of this signature by gene-expression analysis of short hairpin RNA-mediated inhibition of oncogenic Kras2. These experiments identified both a pattern of gene expression indicative of KRAS2 mutation and potential effectors of oncogenic KRAS2 activity in human cancer. This approach provides a strategy for using genomic analysis of animal models to probe human disease.
Wang J, Iwasaki H, Krivtsov A, Febbo PG, Thorner AR, Ernst P, Anastasiadou E, Kutok JL, Kogan SC, Zinkel SS, et al. Conditional MLL-CBP targets GMP and models therapy-related myeloproliferative disease. EMBO J. 2005;24:368–81.
Chromosomal translocations that fuse the mixed lineage leukemia (MLL) gene with multiple partners typify acute leukemias of infancy as well as therapy-related leukemias. We utilized a conditional knockin strategy to bypass the embryonic lethality caused by MLL-CBP expression and to assess the immediate effects of induced MLL-CBP expression on hematopoiesis. Within days of activating MLL-CBP, the fusion protein selectively expanded granulocyte/macrophage progenitors (GMP) and enhanced their self-renewal/proliferation. MLL-CBP altered the gene expression program of GMP, upregulating a subset of genes including Hox a9. Inhibition of Hox a9 expression by RNA interference demonstrated that MLL-CBP required Hox a9 for its enhanced cell expansion. Following exposure to sublethal gamma-irradiation or N-ethyl-N-nitrosourea (ENU), MLL-CBP mice developed myelomonocytic hyperplasia and progressed to fatal myeloproliferative disorders. These represented the spectrum of therapy-induced acute myelomonocytic leukemia/chronic myelomonocytic leukemia/myelodysplastic/myeloproliferative disorder similar to that seen in humans possessing the t(11;16). This model of MLL-CBP therapy-related myeloproliferative disease demonstrates the selectivity of this MLL fusion for GMP cells and its ability to initiate leukemogenesis in conjunction with cooperating mutations.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M. Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature. 2005;434:338–45.
Comprehensive identification of all functional elements encoded in the human genome is a fundamental need in biomedical research. Here, we present a comparative analysis of the human, mouse, rat and dog genomes to create a systematic catalogue of common regulatory motifs in promoters and 3' untranslated regions (3' UTRs). The promoter analysis yields 174 candidate motifs, including most previously known transcription-factor binding sites and 105 new motifs. The 3'-UTR analysis yields 106 motifs likely to be involved in post-transcriptional regulation. Nearly one-half are associated with microRNAs (miRNAs), leading to the discovery of many new miRNA genes and their likely target genes. Our results suggest that previous estimates of the number of human miRNA genes were low, and that miRNAs regulate at least 20% of human genes. The overall results provide a systematic view of gene regulation in the human, which will be refined as additional mammalian genomes become available.
Diamond-Blackfan anemia (DBA), a congenital erythroblastopenia, is a model disease for the study of erythroid differentiation but is poorly understood. RPS19 is the only gene yet to have been associated with DBA, but its relevance to erythroid differentiation is unclear. The molecular basis for the stimulation of erythropoiesis by glucocorticoids in patients with DBA has not been identified. We demonstrate that targeted degradation of the RPS19 transcript, through retroviral expression of short hairpin RNAs (shRNAs), blocks the proliferation and differentiation of erythroid progenitor cells in cultured human CD34(+) cells. Treatment of RPS19-deficient cells with dexamethasone restores erythroid differentiation to normal levels. We investigated the molecular basis of pharmacologic therapies for DBA using oligonucleotide microarrays to survey gene expression in CD34(+) cells treated with combinations of dexamethasone, erythropoietin, stem cell factor, and interleukin-3. Dexamethasone did not alter expression of RPS19 but activated a genetic program that includes a set of key hematopoietic regulatory genes. Genes specific to erythroid progenitor cells were up-regulated by dexamethasone, while genes specific to nonerythroid lineages were down-regulated. Deficiency of RPS19 therefore blocks proliferation of immature erythroid progenitor cells, and dexamethasone activates proliferation of the same cell population through mechanisms independent of RPS19.
Febbo PG, Lowenberg M, Thorner AR, Brown M, Loda M, Golub TR. Androgen mediated regulation and functional implications of fkbp51 expression in prostate cancer. J Urol. 2005;173:1772–7.
PURPOSE: Androgen ablation continues to be the most effective therapy for metastatic prostate cancer, although the biologically active androgen receptor (AR) target genes remain largely unknown. Because AR signaling continues in hormone refractory disease, effector AR target genes may have therapeutic import. MATERIALS AND METHODS: We used oligonucleotide microarrays to identify genes with expression induced by androgen and associated with androgen independent growth. The androgen induced expression of FKBP51, a steroid receptor chaperone, was further investigated in LNCaP cells by Northern and Western analysis, and in primary prostate specimens using immunohistochemistry. We used stable clones over expressing FKBP51 to test the functional effects of FKBP51. RESULTS: Many genes had expression that correlates with androgen stimulation in LNCaP cells but relatively few had reproducible, androgen mediated changes in expression across multiple prostate cancer cell lines. FKBP51 had androgen induced RNA and protein expression in LNCaP cells and decreased expression in normal prostate epithelial cells following castration. Further study demonstrated that FKBP51 induction was not a generalized response to cell proliferation, FKBP51 protein physically interacts with AR and LNCaP cells constitutively over expressing FKBP51 have increased ligand mediated AR activation of an exogenous AR reporter construct and endogenous prostate specific antigen. CONCLUSIONS: Taken together these results confirm FKBP51 as an androgen induced gene, demonstrate a physical interaction between FKBP51 and AR and suggest that FKBP51 over expression increases AR transcriptional activity in prostate cancer.
Feuerhake F, Kutok JL, Monti S, Chen W, LaCasce AS, Cattoretti G, Kurtin P, Pinkus GS, Leval L, Harris NL, et al. NFkappaB activity, function, and target-gene signatures in primary mediastinal large B-cell lymphoma and diffuse large B-cell lymphoma subtypes. Blood. 2005;106:1392–9.
Primary mediastinal large B-cell lymphoma (MLBCL) shares important clinical and molecular features with classic Hodgkin lymphoma, including nuclear localization of the nuclear factor kappaB (NFkappaB) subunit c-REL (reticuloendotheliosis viral oncogene homolog) in a pilot series. Herein, we analyzed c-REL subcellular localization in additional primary MLBCLs and characterized NFkappaB activity and function in a MLBCL cell line. The new primary MLBCLs had prominent c-REL nuclear staining, and the MLBCL cell line exhibited high levels of NFkappaB binding activity. MLBCL cells expressing a superrepressor form of inhibitor of kappa B alpha signaling (IkappaB alpha) had a markedly higher rate of apoptosis, implicating constitutive NFkappaB activity in MLBCL cell survival. The transcriptional profiles of newly diagnosed primary MLBCLs and diffuse large B-cell lymphomas (DLBCLs) were then used to characterize the NFkappaB target gene signatures of MLBCL and specific DLBCL subtypes. MLBCLs expressed increased levels of NFkappaB targets that promote cell survival and favor antiapoptotic tumor necrosis factor alpha (TNFalpha) signaling. In contrast, activated B cell (ABC)-like DLBCLs had a more restricted, potentially developmentally regulated, NFkappaB target gene signature. Of interest, the newly characterized host response DLBCL subtype had a robust NFkappaB target gene signature that partially overlapped that of primary MLBCL. In this large series of primary MLBCLs and DLBCLs, NFkappaB activation was not associated with amplification of the cREL locus, suggesting alternative pathogenetic mechanisms.
Garraway LA, Widlund HR, Rubin MA, Getz G, Berger AJ, Ramaswamy S, Beroukhim R, Milner DA, Granter SR, Du J, et al. Integrative genomic analyses identify MITF as a lineage survival oncogene amplified in malignant melanoma. Nature. 2005;436:117–22.
Systematic analyses of cancer genomes promise to unveil patterns of genetic alterations linked to the genesis and spread of human cancers. High-density single-nucleotide polymorphism (SNP) arrays enable detailed and genome-wide identification of both loss-of-heterozygosity events and copy-number alterations in cancer. Here, by integrating SNP array-based genetic maps with gene expression signatures derived from NCI60 cell lines, we identified the melanocyte master regulator MITF (microphthalmia-associated transcription factor) as the target of a novel melanoma amplification. We found that MITF amplification was more prevalent in metastatic disease and correlated with decreased overall patient survival. BRAF mutation and p16 inactivation accompanied MITF amplification in melanoma cell lines. Ectopic MITF expression in conjunction with the BRAF(V600E) mutant transformed primary human melanocytes, and thus MITF can function as a melanoma oncogene. Reduction of MITF activity sensitizes melanoma cells to chemotherapeutic agents. Targeting MITF in combination with BRAF or cyclin-dependent kinase inhibitors may offer a rational therapeutic avenue into melanoma, a highly chemotherapy-resistant neoplasm. Together, these data suggest that MITF represents a distinct class of 'lineage survival' or 'lineage addiction' oncogenes required for both tissue-specific cancer development and tumour progression.
Levine RL, Wadleigh M, Cools J, Ebert BL, Wernig G, Huntly BJ, Boggon TJ, Wlodarska I, Clark JJ, Moore S, et al. Activating mutation in the tyrosine kinase JAK2 in polycythemia vera, essential thrombocythemia, and myeloid metaplasia with myelofibrosis. Cancer Cell. 2005;7:387–97.
Polycythemia vera (PV), essential thrombocythemia (ET), and myeloid metaplasia with myelofibrosis (MMM) are clonal disorders arising from hematopoietic progenitors. An internet-based protocol was used to collect clinical information and biological specimens from patients with these diseases. High-throughput DNA resequencing identified a recurrent somatic missense mutation JAK2V617F in granulocyte DNA samples of 121 of 164 PV patients, of which 41 had homozygous and 80 had heterozygous mutations. Molecular and cytogenetic analyses demonstrated that homozygous mutations were due to duplication of the mutant allele. JAK2V617F was also identified in granulocyte DNA samples from 37 of 115 ET and 16 of 46 MMM patients, but was not observed in 269 normal individuals. In vitro analysis demonstrated that JAK2V617F is a constitutively active tyrosine kinase.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al. MicroRNA expression profiles classify human cancers. Nature. 2005;435:834–8.
Recent work has revealed the existence of a class of small non-coding RNA species, known as microRNAs (miRNAs), which have critical functions across various biological processes. Here we use a new, bead-based flow cytometric miRNA expression profiling method to present a systematic expression analysis of 217 mammalian miRNAs from 334 samples, including multiple human cancers. The miRNA profiles are surprisingly informative, reflecting the developmental lineage and differentiation state of the tumours. We observe a general downregulation of miRNAs in tumours compared with normal tissues. Furthermore, we were able to successfully classify poorly differentiated tumours using miRNA expression profiles, whereas messenger RNA profiles were highly inaccurate when applied to the same samples. These findings highlight the potential of miRNA profiling in cancer diagnosis.