Prakash Rao

Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification.

Sequencing of DNase I hypersensitive sites (DNase-seq) is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize performance of DNase-seq. Sequencing short fragments of 50–100 base pairs (bp) that accumulate in long internucleosome linker regions was more efficient for identifying transcription factor binding sites compared to sequencing longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy via generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNase I cutting bias, we found a strong effect that varied over more than two orders of magnitude. This indicates that the nucleotide-resolution cleavage patterns at many transcription factor binding sites are derived from intrinsic DNase I cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNase I hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.

The Public Repository of Xenografts Enables Discovery and Randomized Phase II-like Trials in Mice

More than 90% of drugs with preclinical activity fail in human trials, largely due to insufficient efficacy. We hypothesized that adequately powered trials of patient-derived xenografts (PDX) in mice could efficiently define therapeutic activity across heterogeneous tumors. To address this hypothesis, we established a large, publicly available repository of well-characterized leukemia and lymphoma PDXs that undergo orthotopic engraftment, called the Public Repository of Xenografts (PRoXe). PRoXe includes all de-identified information relevant to the primary specimens and the PDXs derived from them. Using this repository, we demonstrate that large studies of acute leukemia PDXs that mimic human randomized clinical trials can characterize drug efficacy and generate transcriptional, functional, and proteomic biomarkers in both treatment-naive and relapsed/refractory disease.

Oncogenic deregulation of EZH2 as an opportunity for targeted therapy in lung cancer

UNLABELLED As a master regulator of chromatin function, the lysine methyltransferase EZH2 orchestrates transcriptional silencing of developmental gene networks. Overexpression of EZH2 is commonly observed in human epithelial cancers, such as non-small cell lung carcinoma (NSCLC), yet definitive demonstration of malignant transformation by deregulated EZH2 remains elusive. Here, we demonstrate the causal role of EZH2 overexpression in NSCLC with new genetically engineered mouse models of lung adenocarcinoma. Deregulated EZH2 silences normal developmental pathways, leading to epigenetic transformation independent of canonical growth factor pathway activation. As such, tumors feature a transcriptional program distinct from KRAS- and EGFR-mutant mouse lung cancers, but shared with human lung adenocarcinomas exhibiting high EZH2 expression. To target EZH2-dependent cancers, we developed a potent open-source EZH2 inhibitor, JQEZ5, that promoted the regression of EZH2-driven tumors in vivo, confirming oncogenic addiction to EZH2 in established tumors and providing the rationale for epigenetic therapy in a subset of lung cancer. SIGNIFICANCE EZH2 overexpression induces murine lung cancers that are similar to human NSCLC with high EZH2 expression and low levels of phosphorylated AKT and ERK, implicating biomarkers for EZH2 inhibitor sensitivity. Our EZH2 inhibitor, JQEZ5, promotes regression of these tumors, revealing a potential role for anti-EZH2 therapy in lung cancer. Cancer Discov; 6(9); 1006-21. ©2016 AACR.See related commentary by Frankel et al., p. 949This article is highlighted in the In This Issue feature, p. 932.

Chromatin immunoprecipitation from fixed clinical tissues reveals tumor-specific enhancer profiles

Extensive cross-linking introduced during routine tissue fixation of clinical pathology specimens severely hampers chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) analysis from archived tissue samples. This limits the ability to study the epigenomes of valuable, clinically annotated tissue resources. Here we describe fixed-tissue chromatin immunoprecipitation sequencing (FiT-seq), a method that enables reliable extraction of soluble chromatin from formalin-fixed paraffin-embedded (FFPE) tissue samples for accurate detection of histone marks. We demonstrate that FiT-seq data from FFPE specimens are concordant with ChIP-seq data from fresh-frozen samples of the same tumors. By using multiple histone marks, we generate chromatin-state maps and identify cis-regulatory elements in clinical samples from various tumor types that can readily allow us to distinguish between cancers by the tissue of origin. Tumor-specific enhancers and superenhancers that are elucidated by FiT-seq analysis correlate with known oncogenic drivers in different tissues and can assist in the understanding of how chromatin states affect gene regulation.

Erratum: The Public Repository of Xenografts Enables Discovery and Randomized Phase II-like Trials in Mice (Cancer Cell (2016) 29 (574–586))

Elizabeth C. Townsend, Mark A. Murakami, Alexandra Christodoulou, Amanda L. Christie, Johannes Köster, Tiffany A. DeSouza, Elizabeth A. Morgan, Scott P. Kallgren, Huiyun Liu, Shuo-Chieh Wu, Olivia Plana, Joan Montero, Kristen E. Stevenson, Prakash Rao, Raga Vadhi, Michael Andreeff, Philippe Armand, Karen K. Ballen, Patrizia Barzaghi-Rinaudo, Sarah Cahill, Rachael A. Clark, Vesselina G. Cooke, Matthew S. Davids, Daniel J. DeAngelo, David M. Dorfman, Hilary Eaton, Benjamin L. Ebert, Julia Etchin, Brant Firestone, David C. Fisher, Arnold S. Freedman, Ilene A. Galinsky, Hui Gao, Jacqueline S. Garcia, Francine Garnache-Ottou, Timothy A. Graubert, Alejandro Gutierrez, Ensar Halilovic, Marian H. Harris, Zachary T. Herbert, Steven M. Horwitz, Giorgio Inghirami, Andrew M. Intlekofer, Moriko Ito, Shai Izraeli, Eric D. Jacobsen, Caron A. Jacobson, Sébastien Jeay, Irmela Jeremias, Michelle A. Kelliher, Raphael Koch, Marina Konopleva, Nadja Kopp, Steven M. Kornblau, Andrew L. Kung, Thomas S. Kupper, Nicole R. LeBoeuf, Ann S. LaCasce, Emma Lees, Loretta S. Li, A. Thomas Look, Masato Murakami, Markus Muschen, Donna Neuberg, Samuel Y. Ng, Oreofe O. Odejide, Stuart H. Orkin, Rachel R. Paquette, Andrew E. Place, Justine E. Roderick, Jeremy A. Ryan, Stephen E. Sallan, Brent Shoji, Lewis B. Silverman, Robert J. Soiffer, David P. Steensma, Kimberly Stegmaier, Richard M. Stone, Jerome Tamburini, Aaron R. Thorner, Paul van Hummelen, Martha Wadleigh, Marion Wiesmann, Andrew P. Weng, Jens U. Wuerthner, David A. Williams, Bruce M. Wollison, Andrew A. Lane, Anthony Letai, Monica M. Bertagnolli, Jerome Ritz, Myles Brown, Henry Long, Jon C. Aster, Margaret A. Shipp, James D. Griffin, and David M. Weinstock* *Correspondence: dweinstock@partners.org http://dx.doi.org/10.1016/j.ccell.2016.06.008

CistromeFinder for ChIP-seq and DNase-seq data reuse

SUMMARY Chromatin immunoprecipitation and DNase I hypersensitivity assays with high-throughput sequencing have greatly accelerated the understanding of transcriptional and epigenetic regulation, although data reuse for the community of experimental biologists has been challenging. We created a data portal CistromeFinder that can help query, evaluate and visualize publicly available Chromatin immunoprecipitation and DNase I hypersensitivity assays with high-throughput sequencing data in human and mouse. The database currently contains 6378 samples over 4391 datasets, 313 factors and 102 cell lines or cell populations. Each dataset has gone through a consistent analysis and quality control pipeline; therefore, users could evaluate the overall quality of each dataset before examining binding sites near their genes of interest. CistromeFinder is integrated with UCSC genome browser for visualization, Primer3Plus for ChIP-qPCR primer design and CistromeMap for submitting newly available datasets. It also allows users to leave comments to facilitate data evaluation and update. AVAILABILITY http://cistrome.org/finder. CONTACT xsliu@jimmy.harvard.edu or henry_long@dfci.harvard.edu.

Embryonic transcription factor SOX9 drives breast cancer endocrine resistance

Significance Resistance to endocrine treatment remains a significant clinical obstacle. ESR1 mutations were found to be the mechanism of endocrine resistance in a substantial number of patients with metastatic ER-positive breast. However, these mutations are primarily linked to aromatase inhibitor resistance and are not strongly associated with tamoxifen resistance. Herein, we show that tamoxifen treatment promotes a RUNX2–ER complex, which mediates an altered ER cistrome that facilitates the up-regulation of SOX9. We show that up-regulation of SOX9, an embryonic transcription factor with key roles in metastases, is a driver of endocrine resistance in the setting of tamoxifen treatment. Our data provide putative targets for the development of new strategies to treat tamoxifen-resistant breast cancer. The estrogen receptor (ER) drives the growth of most luminal breast cancers and is the primary target of endocrine therapy. Although ER blockade with drugs such as tamoxifen is very effective, a major clinical limitation is the development of endocrine resistance especially in the setting of metastatic disease. Preclinical and clinical observations suggest that even following the development of endocrine resistance, ER signaling continues to exert a pivotal role in tumor progression in the majority of cases. Through the analysis of the ER cistrome in tamoxifen-resistant breast cancer cells, we have uncovered a role for an RUNX2–ER complex that stimulates the transcription of a set of genes, including most notably the stem cell factor SOX9, that promote proliferation and a metastatic phenotype. We show that up-regulation of SOX9 is sufficient to cause relative endocrine resistance. The gain of SOX9 as an ER-regulated gene associated with tamoxifen resistance was validated in a unique set of clinical samples supporting the need for the development of improved ER antagonists.

Progesterone receptor isoforms, agonists and antagonists differentially reprogram estrogen signaling

Major roadblocks to developing effective progesterone receptor (PR)-targeted therapies in breast cancer include the lack of highly-specific PR modulators, a poor understanding of the pro- or anti-tumorigenic networks for PR isoforms and ligands, and an incomplete understanding of the cross talk between PR and estrogen receptor (ER) signaling. Through genomic analyses of xenografts treated with various clinically-relevant ER and PR-targeting drugs, we describe how the activation or inhibition of PR differentially reprograms estrogen signaling, resulting in the segregation of transcriptomes into separate PR agonist and antagonist-mediated groups. These findings address an ongoing controversy regarding the clinical utility of PR agonists and antagonists, alone or in combination with tamoxifen, for breast cancer management. Additionally, the two PR isoforms PRA and PRB, bind distinct but overlapping genomic sites and interact with different sets of co-regulators to differentially modulate estrogen signaling to be either pro- or anti-tumorigenic. Of the two isoforms, PRA inhibited gene expression and ER chromatin binding significantly more than PRB. Differential gene expression was observed in PRA and PRB-rich patient tumors and PRA-rich gene signatures had poorer survival outcomes. In support of antiprogestin responsiveness of PRA-rich tumors, gene signatures associated with PR antagonists, but not PR agonists, predicted better survival outcomes. The better patient survival associated with PR antagonists versus PR agonists treatments was further reflected in the higher in vivo anti-tumor activity of therapies that combine tamoxifen with PR antagonists and modulators. This study suggests that distinguishing common effects observed due to concomitant interaction of another receptor with its ligand (agonist or antagonist), from unique isoform and ligand-specific effects will inform the development of biomarkers for patient selection and translation of PR-targeted therapies to the clinic.

Abstract LB-257: Estrogen signaling in mature luminal and luminal progenitor cells of BRCA2 carriers and non-carriers

Background: Carriers of mutations in the breast cancer susceptibility gene BRCA2 have a 50-80% lifetime risk of developing breast cancers, about 70% of which are positive for estrogen receptor (ER), a key oncogenic hormone receptor. Metastatic ER+ tumors also respond to PARP inhibitors associated with their DNA repair defect. There remains much to understand regarding estrogen signaling in BRCA2 carriers and preventing cancers in this high-risk population. Methods: We studied the transcriptomes and in-vitro growth of luminal progenitors and mature luminal cells derived from primary human normal breast tissue obtained from BRCA2 carriers and wild-type patients undergoing prophylactic mastectomy and reduction mammoplasty procedures, respectively. Subsequently, formalin-fixed paraffin-embedded breast cancer tissue samples from twelve case-control matched pairs of BRCA2 carriers and wild-type patients were investigated for ER genomic binding, H3K27Ac enhancer marks, and gene expression changes. Results: Luminal progenitors and mature luminal cells have different gene expression patterns. Unsupervised clustering of transcriptomes from normal tissue of three pairs of BRCA2-mutant and wildtype patients suggested that BRCA2 mutation status can segregate transcriptomes from luminal progenitors but not mature luminal cells. These differential transcriptomes in luminal progenitor cells from BRCA2 carriers were enriched for cell proliferation genes. In agreement, greater proliferative phenotype was observed in BRCA2-mutant luminal progenitor cells in comparison to the wild-type cohort. Interestingly, clustering analyses of twelve pairs of breast tumors indicated no clear differences between ER cistromes, enhancer marks and gene expression between BRCA2 carriers and wildtype patients. Because luminal progenitor cells are more likely to contribute to tumorigenesis, we hypothesize that the identified BRCA2-mutant signature in normal progenitor cells could be retained in tumors obtained from BRCA2 carriers. Conclusions: There are differential gene expression patterns between luminal progenitors and mature luminal cells obtained from normal human breast tissue. The transcriptomes in the luminal progenitor but not mature luminal cells segregate differently based on the BRCA2 mutation status. This identified BRCA2-mutant gene signature in luminal progenitor cells is enriched for pro-proliferation pathways, and BRCA2-mutant luminal progenitor cells exhibit higher cell proliferation, suggesting a potential cell-of-origin in BRCA2-associated breast cancers. Importantly, no significant differences in transcriptomes and ER cistromes are observed in breast tumors from BRCA2 carriers and non-carriers, indicating a need to identify BRCA2-mutant gene signature that may assist in segregating ER+ breast cancers in this high-risk populations. The knowledge from this study may assist in pioneering strategies for cancer prevention in high-risk BRCA2 mutation carriers. Citation Format: Hari Singhal, David Chi, Elgene Lim, Housheng He, Raga Vadhi, Prakash K. Rao, Henry W. Long, Andrea L. Richardson, Judy Garber, Myles Brown. Estrogen signaling in mature luminal and luminal progenitor cells of BRCA2 carriers and non-carriers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-257. doi:10.1158/1538-7445.AM2017-LB-257