Benjamin L. Ebert

Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles

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.

MicroRNA expression profiles classify human cancers.

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.

The landscape of somatic copy-number alteration across human cancers

A powerful way to discover key genes with causal roles in oncogenesis is to identify genomic regions that undergo frequent alteration in human cancers. Here we present high-resolution analyses of somatic copy-number alterations (SCNAs) from 3,131 cancer specimens, belonging largely to 26 histological types. We identify 158 regions of focal SCNA that are altered at significant frequency across several cancer types, of which 122 cannot be explained by the presence of a known cancer target gene located within these regions. Several gene families are enriched among these regions of focal SCNA, including the BCL2 family of apoptosis regulators and the NF-κΒ pathway. We show that cancer cells containing amplifications surrounding the MCL1 and BCL2L1 anti-apoptotic genes depend on the expression of these genes for survival. Finally, we demonstrate that a large majority of SCNAs identified in individual cancer types are present in several cancer types.

Genome-scale CRISPR-Cas9 knockout screening in human cells.

Improving Whole-Genome Screens Improved methods are needed for the knockout of individual genes in genome-scale functional screens. Wang et al. (p. 80, published online 12 December) and Shalem et al. (p. 84, published online 12 December) used the bacterial CRISPR/Cas9 system to power-screen protocols that avoid several of the pitfalls associated with small interfering RNA (siRNA) screens. Genome editing by these methods completely disrupts target genes, thus avoiding weak signals that can occur when transcript abundance is partially decreased by siRNA. Furthermore, gene targeting by the CRISPR system is more precise and appears to produce substantially fewer off-target effects than existing methods. Genome-editing technology allows improved positive or negative selection screens. The simplicity of programming the CRISPR (clustered regularly interspaced short palindromic repeats)–associated nuclease Cas9 to modify specific genomic loci suggests a new way to interrogate gene function on a genome-wide scale. We show that lentiviral delivery of a genome-scale CRISPR-Cas9 knockout (GeCKO) library targeting 18,080 genes with 64,751 unique guide sequences enables both negative and positive selection screening in human cells. First, we used the GeCKO library to identify genes essential for cell viability in cancer and pluripotent stem cells. Next, in a melanoma model, we screened for genes whose loss is involved in resistance to vemurafenib, a therapeutic RAF inhibitor. Our highest-ranking candidates include previously validated genes NF1 and MED12, as well as novel hits NF2, CUL3, TADA2B, and TADA1. We observe a high level of consistency between independent guide RNAs targeting the same gene and a high rate of hit confirmation, demonstrating the promise of genome-scale screening with Cas9.

MPLW515L Is a Novel Somatic Activating Mutation in Myelofibrosis with Myeloid Metaplasia

Background The JAK2V617F allele has recently been identified in patients with polycythemia vera (PV), essential thrombocytosis (ET), and myelofibrosis with myeloid metaplasia (MF). Subsequent analysis has shown that constitutive activation of the JAK-STAT signal transduction pathway is an important pathogenetic event in these patients, and that enzymatic inhibition of JAK2V617F may be of therapeutic benefit in this context. However, a significant proportion of patients with ET or MF are JAK2V617F-negative. We hypothesized that activation of the JAK-STAT pathway might also occur as a consequence of activating mutations in certain hematopoietic-specific cytokine receptors, including the erythropoietin receptor (EPOR), the thrombopoietin receptor (MPL), or the granulocyte-colony stimulating factor receptor (GCSFR). Methods and Findings DNA sequence analysis of the exons encoding the transmembrane and juxtamembrane domains of EPOR, MPL, and GCSFR, and comparison with germline DNA derived from buccal swabs, identified a somatic activating mutation in the transmembrane domain of MPL (W515L) in 9% (4/45) of JAKV617F-negative MF. Expression of MPLW515L in 32D, UT7, or Ba/F3 cells conferred cytokine-independent growth and thrombopoietin hypersensitivity, and resulted in constitutive phosphorylation of JAK2, STAT3, STAT5, AKT, and ERK. Furthermore, a small molecule JAK kinase inhibitor inhibited MPLW515L-mediated proliferation and JAK-STAT signaling in vitro. In a murine bone marrow transplant assay, expression of MPLW515L, but not wild-type MPL, resulted in a fully penetrant myeloproliferative disorder characterized by marked thrombocytosis (Plt count 1.9–4.0 × 10 12/L), marked splenomegaly due to extramedullary hematopoiesis, and increased reticulin fibrosis. Conclusions Activation of JAK-STAT signaling via MPLW515L is an important pathogenetic event in patients with JAK2V617F-negative MF. The bone marrow transplant model of MPLW515L-mediated myeloproliferative disorders (MPD) exhibits certain features of human MF, including extramedullary hematopoiesis, splenomegaly, and megakaryocytic proliferation. Further analysis of positive and negative regulators of the JAK-STAT pathway is warranted in JAK2V617F-negative MPD.

Clinical effect of point mutations in myelodysplastic syndromes

BACKGROUND Myelodysplastic syndromes are clinically heterogeneous disorders characterized by clonal hematopoiesis, impaired differentiation, peripheral-blood cytopenias, and a risk of progression to acute myeloid leukemia. Somatic mutations may influence the clinical phenotype but are not included in current prognostic scoring systems. METHODS We used a combination of genomic approaches, including next-generation sequencing and mass spectrometry-based genotyping, to identify mutations in samples of bone marrow aspirate from 439 patients with myelodysplastic syndromes. We then examined whether the mutation status for each gene was associated with clinical variables, including specific cytopenias, the proportion of blasts, and overall survival. RESULTS We identified somatic mutations in 18 genes, including two, ETV6 and GNAS, that have not been reported to be mutated in patients with myelodysplastic syndromes. A total of 51% of all patients had at least one point mutation, including 52% of the patients with normal cytogenetics. Mutations in RUNX1, TP53, and NRAS were most strongly associated with severe thrombocytopenia (P<0.001 for all comparisons) and an increased proportion of bone marrow blasts (P<0.006 for all comparisons). In a multivariable Cox regression model, the presence of mutations in five genes retained independent prognostic significance: TP53 (hazard ratio for death from any cause, 2.48; 95% confidence interval [CI], 1.60 to 3.84), EZH2 (hazard ratio, 2.13; 95% CI, 1.36 to 3.33), ETV6 (hazard ratio, 2.04; 95% CI, 1.08 to 3.86), RUNX1 (hazard ratio, 1.47; 95% CI, 1.01 to 2.15), and ASXL1 (hazard ratio, 1.38; 95% CI, 1.00 to 1.89). CONCLUSIONS Somatic point mutations are common in myelodysplastic syndromes and are associated with specific clinical features. Mutations in TP53, EZH2, ETV6, RUNX1, and ASXL1 are predictors of poor overall survival in patients with myelodysplastic syndromes, independently of established risk factors. (Funded by the National Institutes of Health and others.).

Age-Related Clonal Hematopoiesis Associated with Adverse Outcomes

BACKGROUND The incidence of hematologic cancers increases with age. These cancers are associated with recurrent somatic mutations in specific genes. We hypothesized that such mutations would be detectable in the blood of some persons who are not known to have hematologic disorders. METHODS We analyzed whole-exome sequencing data from DNA in the peripheral-blood cells of 17,182 persons who were unselected for hematologic phenotypes. We looked for somatic mutations by identifying previously characterized single-nucleotide variants and small insertions or deletions in 160 genes that are recurrently mutated in hematologic cancers. The presence of mutations was analyzed for an association with hematologic phenotypes, survival, and cardiovascular events. RESULTS Detectable somatic mutations were rare in persons younger than 40 years of age but rose appreciably in frequency with age. Among persons 70 to 79 years of age, 80 to 89 years of age, and 90 to 108 years of age, these clonal mutations were observed in 9.5% (219 of 2300 persons), 11.7% (37 of 317), and 18.4% (19 of 103), respectively. The majority of the variants occurred in three genes: DNMT3A, TET2, and ASXL1. The presence of a somatic mutation was associated with an increase in the risk of hematologic cancer (hazard ratio, 11.1; 95% confidence interval [CI], 3.9 to 32.6), an increase in all-cause mortality (hazard ratio, 1.4; 95% CI, 1.1 to 1.8), and increases in the risks of incident coronary heart disease (hazard ratio, 2.0; 95% CI, 1.2 to 3.4) and ischemic stroke (hazard ratio, 2.6; 95% CI, 1.4 to 4.8). CONCLUSIONS Age-related clonal hematopoiesis is a common condition that is associated with increases in the risk of hematologic cancer and in all-cause mortality, with the latter possibly due to an increased risk of cardiovascular disease. (Funded by the National Institutes of Health and others.).

Bone progenitor dysfunction induces myelodysplasia and secondary leukaemia

Mesenchymal cells contribute to the ‘stroma’ of most normal and malignant tissues, with specific mesenchymal cells participating in the regulatory niches of stem cells. By examining how mesenchymal osteolineage cells modulate haematopoiesis, here we show that deletion of Dicer1 specifically in mouse osteoprogenitors, but not in mature osteoblasts, disrupts the integrity of haematopoiesis. Myelodysplasia resulted and acute myelogenous leukaemia emerged that had acquired several genetic abnormalities while having intact Dicer1. Examining gene expression altered in osteoprogenitors as a result of Dicer1 deletion showed reduced expression of Sbds, the gene mutated in Schwachman–Bodian–Diamond syndrome—a human bone marrow failure and leukaemia pre-disposition condition. Deletion of Sbds in mouse osteoprogenitors induced bone marrow dysfunction with myelodysplasia. Therefore, perturbation of specific mesenchymal subsets of stromal cells can disorder differentiation, proliferation and apoptosis of heterologous cells, and disrupt tissue homeostasis. Furthermore, primary stromal dysfunction can result in secondary neoplastic disease, supporting the concept of niche-induced oncogenesis.

Identification of RPS14 as a 5q- syndrome gene by RNA interference screen.

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.

Rational design of highly active sgRNAs for CRISPR-Cas9–mediated gene inactivation

Components of the prokaryotic clustered, regularly interspaced, short palindromic repeats (CRISPR) loci have recently been repurposed for use in mammalian cells. The CRISPR-associated (Cas)9 can be programmed with a single guide RNA (sgRNA) to generate site-specific DNA breaks, but there are few known rules governing on-target efficacy of this system. We created a pool of sgRNAs, tiling across all possible target sites of a panel of six endogenous mouse and three endogenous human genes and quantitatively assessed their ability to produce null alleles of their target gene by antibody staining and flow cytometry. We discovered sequence features that improved activity, including a further optimization of the protospacer-adjacent motif (PAM) of Streptococcus pyogenes Cas9. The results from 1,841 sgRNAs were used to construct a predictive model of sgRNA activity to improve sgRNA design for gene editing and genetic screens. We provide an online tool for the design of highly active sgRNAs for any gene of interest.