1
|
Otto JE, Ursu O, Wu AP, Winter EB, Cuoco MS, Ma S, Qian K, Michel BC, Buenrostro JD, Berger B, Regev A, Kadoch C. Structural and functional properties of mSWI/SNF chromatin remodeling complexes revealed through single-cell perturbation screens. Mol Cell 2023; 83:1350-1367.e7. [PMID: 37028419 DOI: 10.1016/j.molcel.2023.03.013] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2022] [Revised: 02/07/2023] [Accepted: 03/10/2023] [Indexed: 04/09/2023]
Abstract
The mammalian SWI/SNF (mSWI/SNF or BAF) family of chromatin remodeling complexes play critical roles in regulating DNA accessibility and gene expression. The three final-form subcomplexes-cBAF, PBAF, and ncBAF-are distinct in biochemical componentry, chromatin targeting, and roles in disease; however, the contributions of their constituent subunits to gene expression remain incompletely defined. Here, we performed Perturb-seq-based CRISPR-Cas9 knockout screens targeting mSWI/SNF subunits individually and in select combinations, followed by single-cell RNA-seq and SHARE-seq. We uncovered complex-, module-, and subunit-specific contributions to distinct regulatory networks and defined paralog subunit relationships and shifted subcomplex functions upon perturbations. Synergistic, intra-complex genetic interactions between subunits reveal functional redundancy and modularity. Importantly, single-cell subunit perturbation signatures mapped across bulk primary human tumor expression profiles both mirror and predict cBAF loss-of-function status in cancer. Our findings highlight the utility of Perturb-seq to dissect disease-relevant gene regulatory impacts of heterogeneous, multi-component master regulatory complexes.
Collapse
Affiliation(s)
- Jordan E Otto
- Department of Pediatric Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA; Broad Institute of MIT and Harvard, Cambridge, MA, USA; Chemical Biology Program, Harvard University, Cambridge, MA, USA
| | - Oana Ursu
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Department of Biology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Alexander P Wu
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA, USA; Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Evan B Winter
- Department of Pediatric Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA; Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | | | - Sai Ma
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA, USA
| | - Kristin Qian
- Department of Pediatric Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA; Broad Institute of MIT and Harvard, Cambridge, MA, USA; Biological and Biomedical Sciences Program, Harvard Medical School, Boston, MA, USA
| | - Brittany C Michel
- Department of Pediatric Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA; Broad Institute of MIT and Harvard, Cambridge, MA, USA; Biological and Biomedical Sciences Program, Harvard Medical School, Boston, MA, USA
| | - Jason D Buenrostro
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA, USA
| | - Bonnie Berger
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA, USA; Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Aviv Regev
- Broad Institute of MIT and Harvard, Cambridge, MA, USA; Department of Biology, Massachusetts Institute of Technology, Cambridge, MA, USA; Howard Hughes Medical Institute, Chevy Chase, MD, UA.
| | - Cigall Kadoch
- Department of Pediatric Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, MA, USA; Broad Institute of MIT and Harvard, Cambridge, MA, USA; Chemical Biology Program, Harvard University, Cambridge, MA, USA; Howard Hughes Medical Institute, Chevy Chase, MD, UA.
| |
Collapse
|
2
|
Al-Jazrawe M, Molnar C, Rindtorff N, Eser P, Misek S, Alimova M, Ursu O, Colgan W, Attari A, Tsang N, Keskula P, Rios C, Tseng M, Carpenter A, McFarland J, Bass A, Klempner S, Boehm J. Abstract PO-093: Evaluating dependencies by rapid image-based ex vivo cancer biosensors. Clin Cancer Res 2021. [DOI: 10.1158/1557-3265.adi21-po-093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Abstract
Successful mapping of cancer dependencies requires conducting genetic and drug screens on a diversity of models. However, the difficulty in generating long-term models of many cancers limits the share of patient samples that can be studied. Such long-term models have likely also lost the cellular heterogeneity present in the original tumor due to in vitro propagation. To overcome these limitations, we are developing image-based ex vivo cancer biosensors from early patient material. Using freshly received gastroesophageal cancer ascites, we are optimizing perturbation methods and utilizing single-cell transcriptomics and label-free microscopy to infer a subpopulation-specific vulnerability profile. We show that label-free microscopy can infer cell identity and viability in heterogeneous early patient samples. Additionally, early drug perturbation recapitulates observations made in established gastroesophageal cancer organoids. Successful implementation of ex vivo biosensors will expand the cancer dependency space by making perturbational studies accessible to more diverse samples, and by identifying and validating hits in a more immediate setting to the original tumor.
Citation Format: Mushriq Al-Jazrawe, Csaba Molnar, Niklas Rindtorff, Pinar Eser, Sean Misek, Maria Alimova, Oana Ursu, William Colgan, Adel Attari, Natalie Tsang, Paula Keskula, Carmen Rios, Moony Tseng, Anne Carpenter, James McFarland, Adam Bass, Samuel Klempner, Jesse Boehm. Evaluating dependencies by rapid image-based ex vivo cancer biosensors [abstract]. In: Proceedings of the AACR Virtual Special Conference on Artificial Intelligence, Diagnosis, and Imaging; 2021 Jan 13-14. Philadelphia (PA): AACR; Clin Cancer Res 2021;27(5_Suppl):Abstract nr PO-093.
Collapse
Affiliation(s)
| | - Csaba Molnar
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | | | - Pinar Eser
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | - Sean Misek
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | | | - Oana Ursu
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | | | - Adel Attari
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | | | | | - Carmen Rios
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | - Moony Tseng
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| | | | | | - Adam Bass
- 2Dana-Farber Cancer Institute, Boston, MA,
| | | | - Jesse Boehm
- 1Broad Institute of MIT and Harvard, Cambridge, MA,
| |
Collapse
|
3
|
Ursu O, Boley N, Taranova M, Wang YXR, Yardimci GG, Stafford Noble W, Kundaje A. GenomeDISCO: a concordance score for chromosome conformation capture experiments using random walks on contact map graphs. Bioinformatics 2019; 34:2701-2707. [PMID: 29554289 DOI: 10.1093/bioinformatics/bty164] [Citation(s) in RCA: 51] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2017] [Accepted: 03/15/2018] [Indexed: 02/04/2023] Open
Abstract
Motivation The three-dimensional organization of chromatin plays a critical role in gene regulation and disease. High-throughput chromosome conformation capture experiments such as Hi-C are used to obtain genome-wide maps of three-dimensional chromatin contacts. However, robust estimation of data quality and systematic comparison of these contact maps is challenging due to the multi-scale, hierarchical structure of chromatin contacts and the resulting properties of experimental noise in the data. Measuring concordance of contact maps is important for assessing reproducibility of replicate experiments and for modeling variation between different cellular contexts. Results We introduce a concordance measure called DIfferences between Smoothed COntact maps (GenomeDISCO) for assessing the similarity of a pair of contact maps obtained from chromosome conformation capture experiments. The key idea is to smooth contact maps using random walks on the contact map graph, before estimating concordance. We use simulated datasets to benchmark GenomeDISCO's sensitivity to different types of noise that affect chromatin contact maps. When applied to a large collection of Hi-C datasets, GenomeDISCO accurately distinguishes biological replicates from samples obtained from different cell types. GenomeDISCO also generalizes to other chromosome conformation capture assays, such as HiChIP. Availability and implementation Software implementing GenomeDISCO is available at https://github.com/kundajelab/genomedisco. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Oana Ursu
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
| | - Nathan Boley
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
| | - Maryna Taranova
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
| | - Y X Rachel Wang
- Department of Statistics, Stanford University, Stanford, CA, USA
| | | | - William Stafford Noble
- Department of Genome Sciences, University of Washington, WA, USA.,Department of Computer Science and Engineering, University of Washington, WA, USA
| | - Anshul Kundaje
- Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA.,Department of Computer Science, Stanford University, Stanford, CA, USA
| |
Collapse
|
4
|
Tycko J, Wainberg M, Marinov GK, Ursu O, Hess GT, Ego BK, Aradhana, Li A, Truong A, Trevino AE, Spees K, Yao D, Kaplow IM, Greenside PG, Morgens DW, Phanstiel DH, Snyder MP, Bintu L, Greenleaf WJ, Kundaje A, Bassik MC. Mitigation of off-target toxicity in CRISPR-Cas9 screens for essential non-coding elements. Nat Commun 2019; 10:4063. [PMID: 31492858 PMCID: PMC6731277 DOI: 10.1038/s41467-019-11955-7] [Citation(s) in RCA: 75] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2019] [Accepted: 08/07/2019] [Indexed: 12/26/2022] Open
Abstract
Pooled CRISPR-Cas9 screens are a powerful method for functionally characterizing regulatory elements in the non-coding genome, but off-target effects in these experiments have not been systematically evaluated. Here, we investigate Cas9, dCas9, and CRISPRi/a off-target activity in screens for essential regulatory elements. The sgRNAs with the largest effects in genome-scale screens for essential CTCF loop anchors in K562 cells were not single guide RNAs (sgRNAs) that disrupted gene expression near the on-target CTCF anchor. Rather, these sgRNAs had high off-target activity that, while only weakly correlated with absolute off-target site number, could be predicted by the recently developed GuideScan specificity score. Screens conducted in parallel with CRISPRi/a, which do not induce double-stranded DNA breaks, revealed that a distinct set of off-targets also cause strong confounding fitness effects with these epigenome-editing tools. Promisingly, filtering of CRISPRi libraries using GuideScan specificity scores removed these confounded sgRNAs and enabled identification of essential regulatory elements.
Collapse
Affiliation(s)
- Josh Tycko
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Michael Wainberg
- Department of Computer Science, Stanford University, Stanford, CA, 94305, USA
| | - Georgi K Marinov
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Oana Ursu
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Gaelen T Hess
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Braeden K Ego
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Aradhana
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Amy Li
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Alisa Truong
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Alexandro E Trevino
- Center for Personal Dynamic Regulomes, Stanford University, Stanford, CA, 94305, USA
- Department of Bioengineering, Stanford University, Stanford, CA, 94305, USA
| | - Kaitlyn Spees
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - David Yao
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Irene M Kaplow
- Department of Computer Science, Stanford University, Stanford, CA, 94305, USA
- Department of Biology, Stanford University, Stanford, CA, 94305, USA
| | - Peyton G Greenside
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
- Program in Biomedical Informatics, Stanford University School of Medicine, Stanford, CA, 94305, USA
| | - David W Morgens
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Douglas H Phanstiel
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
- Department of Cell Biology and Physiology, University of North Carolina, Chapel Hill, NC, 27599, USA
- Thurston Arthritis Research Center, University of North Carolina, Chapel Hill, NC, 27599, USA
| | - Michael P Snyder
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA
| | - Lacramioara Bintu
- Department of Bioengineering, Stanford University, Stanford, CA, 94305, USA
| | - William J Greenleaf
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA.
- Department of Applied Physics, Stanford University, Stanford, CA, 94305, USA.
- Chan Zuckerberg Biohub, San Francisco, CA, 94158, USA.
| | - Anshul Kundaje
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA.
- Department of Computer Science, Stanford University, Stanford, CA, 94305, USA.
| | - Michael C Bassik
- Department of Genetics, Stanford University, Stanford, CA, 94305, USA.
- Chemistry, Engineering, and Medicine for Human Health (ChEM-H), Stanford University, Stanford, CA, 94305, USA.
| |
Collapse
|
5
|
Abstract
Chromosome conformation capture experiments such as Hi-C are used to map the three-dimensional spatial organization of genomes. One specific feature of the 3D organization is known as topologically associating domains (TADs), which are densely interacting, contiguous chromatin regions playing important roles in regulating gene expression. A few algorithms have been proposed to detect TADs. In particular, the structure of Hi-C data naturally inspires application of community detection methods. However, one of the drawbacks of community detection is that most methods take exchangeability of the nodes in the network for granted; whereas the nodes in this case, that is, the positions on the chromosomes, are not exchangeable. We propose a network model for detecting TADs using Hi-C data that takes into account this nonexchangeability. in addition, our model explicitly makes use of cell-type specific CTCF binding sites as biological covariates and can be used to identify conserved TADs across multiple cell types. The model leads to a likelihood objective that can be efficiently optimized via relaxation. We also prove that when suitably initialized, this model finds the underlying TAD structure with high probability. using simulated data, we show the advantages of our method and the caveats of popular community detection methods, such as spectral clustering, in this application. Applying our method to real Hi-C data, we demonstrate the domains identified have desirable epigenetic features and compare them across different cell types.
Collapse
|
6
|
Ursu O, Gosline SJC, Beeharry N, Fink L, Bhattacharjee V, Huang SSC, Zhou Y, Yen T, Fraenkel E. Network modeling of kinase inhibitor polypharmacology reveals pathways targeted in chemical screens. PLoS One 2017; 12:e0185650. [PMID: 29023490 PMCID: PMC5638242 DOI: 10.1371/journal.pone.0185650] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2017] [Accepted: 09/15/2017] [Indexed: 01/22/2023] Open
Abstract
Small molecule screens are widely used to prioritize pharmaceutical development. However, determining the pathways targeted by these molecules is challenging, since the compounds are often promiscuous. We present a network strategy that takes into account the polypharmacology of small molecules in order to generate hypotheses for their broader mode of action. We report a screen for kinase inhibitors that increase the efficacy of gemcitabine, the first-line chemotherapy for pancreatic cancer. Eight kinase inhibitors emerge that are known to affect 201 kinases, of which only three kinases have been previously identified as modifiers of gemcitabine toxicity. In this work, we use the SAMNet algorithm to identify pathways linking these kinases and genetic modifiers of gemcitabine toxicity with transcriptional and epigenetic changes induced by gemcitabine that we measure using DNaseI-seq and RNA-seq. SAMNet uses a constrained optimization algorithm to connect genes from these complementary datasets through a small set of protein-protein and protein-DNA interactions. The resulting network recapitulates known pathways including DNA repair, cell proliferation and the epithelial-to-mesenchymal transition. We use the network to predict genes with important roles in the gemcitabine response, including six that have already been shown to modify gemcitabine efficacy in pancreatic cancer and ten novel candidates. Our work reveals the important role of polypharmacology in the activity of these chemosensitizing agents.
Collapse
Affiliation(s)
- Oana Ursu
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
| | - Sara J. C. Gosline
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
| | - Neil Beeharry
- Fox Chase Cancer Center, Philadelphia, Pennsylvania, United States of America
| | - Lauren Fink
- Cancer Biology Program, Fox Chase Cancer Center; Philadelphia, Pennsylvania, United States of America
| | | | - Shao-shan Carol Huang
- Plant Biology Laboratory, The Salk Institute for Biological Studies, La Jolla, California, United States of America
- Genomic Analysis Laboratory, The Salk Institute for Biological Studies, La Jolla, California, United States of America
- Howard Hughes Medical Institute, The Salk Institute for Biological Studies, La Jolla, California, United States of America
| | - Yan Zhou
- Fox Chase Cancer Center, Philadelphia, Pennsylvania, United States of America
| | - Tim Yen
- Fox Chase Cancer Center, Philadelphia, Pennsylvania, United States of America
| | - Ernest Fraenkel
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America
- * E-mail:
| |
Collapse
|
7
|
Morgens DW, Wainberg M, Boyle EA, Ursu O, Araya CL, Tsui CK, Haney MS, Hess GT, Han K, Jeng EE, Li A, Snyder MP, Greenleaf WJ, Kundaje A, Bassik MC. Genome-scale measurement of off-target activity using Cas9 toxicity in high-throughput screens. Nat Commun 2017; 8:15178. [PMID: 28474669 PMCID: PMC5424143 DOI: 10.1038/ncomms15178] [Citation(s) in RCA: 206] [Impact Index Per Article: 29.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2016] [Accepted: 03/07/2017] [Indexed: 12/24/2022] Open
Abstract
CRISPR-Cas9 screens are powerful tools for high-throughput interrogation of genome function, but can be confounded by nuclease-induced toxicity at both on- and off-target sites, likely due to DNA damage. Here, to test potential solutions to this issue, we design and analyse a CRISPR-Cas9 library with 10 variable-length guides per gene and thousands of negative controls targeting non-functional, non-genic regions (termed safe-targeting guides), in addition to non-targeting controls. We find this library has excellent performance in identifying genes affecting growth and sensitivity to the ricin toxin. The safe-targeting guides allow for proper control of toxicity from on-target DNA damage. Using this toxicity as a proxy to measure off-target cutting, we demonstrate with tens of thousands of guides both the nucleotide position-dependent sensitivity to single mismatches and the reduction of off-target cutting using truncated guides. Our results demonstrate a simple strategy for high-throughput evaluation of target specificity and nuclease toxicity in Cas9 screens.
Collapse
Affiliation(s)
- David W. Morgens
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Michael Wainberg
- Department of Computer Science, Stanford University, Stanford, California 94305, USA
| | - Evan A. Boyle
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Oana Ursu
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Carlos L. Araya
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - C. Kimberly Tsui
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Michael S. Haney
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Gaelen T. Hess
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Kyuho Han
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Edwin E. Jeng
- Department of Genetics, Stanford University, Stanford, California 94305, USA
- Program in Cancer Biology, Stanford University, Stanford, California 94305, USA
| | - Amy Li
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | - Michael P. Snyder
- Department of Genetics, Stanford University, Stanford, California 94305, USA
| | | | - Anshul Kundaje
- Department of Genetics, Stanford University, Stanford, California 94305, USA
- Department of Computer Science, Stanford University, Stanford, California 94305, USA
| | - Michael C. Bassik
- Department of Genetics, Stanford University, Stanford, California 94305, USA
- Stanford University Chemistry, Engineering, and Medicine for Human Health (ChEM-H), Stanford, California 94305, USA
| |
Collapse
|
8
|
Liu T, Oprea T, Ursu O, Hasselgren C, Altman RB. Estimation of Maximum Recommended Therapeutic Dose Using Predicted Promiscuity and Potency. Clin Transl Sci 2016; 9:311-320. [PMID: 27736015 PMCID: PMC5161261 DOI: 10.1111/cts.12422] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2016] [Accepted: 09/01/2016] [Indexed: 01/08/2023] Open
Abstract
We report a simple model that predicts the maximum recommended therapeutic dose (MRTD) of small molecule drugs based on an assessment of likely protein-drug interactions. Previously, we reported methods for computational estimation of drug promiscuity and potency. We used these concepts to build a linear model derived from 238 small molecular drugs to predict MRTD. We applied this model successfully to predict MRTDs for 16 nonsteroidal antiinflammatory drugs (NSAIDs) and 14 antiretroviral drugs. Of note, based on the estimated promiscuity of low-dose drugs (and active chemicals), we identified 83 proteins as "high-risk off-targets" (HROTs) that are often associated with low doses; the evaluation of interactions with HROTs may be useful during early phases of drug discovery. Our model helps explain the MRTD for drugs with severe adverse reactions caused by interactions with HROTs.
Collapse
Affiliation(s)
- T Liu
- Department of Genetics, Stanford University, Stanford, California, USA
| | - T Oprea
- Department of Internal Medicine, Translational Informatics Division, University of New Mexico, Albuquerque, New Mexico, USA
| | - O Ursu
- Department of Internal Medicine, Translational Informatics Division, University of New Mexico, Albuquerque, New Mexico, USA
| | - C Hasselgren
- PureInfo Discovery Corp, Albuquerque, New Mexico, USA
| | - R B Altman
- Departments of Bioengineering and Genetics, Stanford University, Stanford, California, USA
| |
Collapse
|
9
|
Grubert F, Zaugg JB, Kasowski M, Ursu O, Spacek DV, Martin AR, Greenside P, Srivas R, Phanstiel DH, Pekowska A, Heidari N, Euskirchen G, Huber W, Pritchard JK, Bustamante CD, Steinmetz LM, Kundaje A, Snyder M. Genetic Control of Chromatin States in Humans Involves Local and Distal Chromosomal Interactions. Cell 2015; 162:1051-65. [PMID: 26300125 PMCID: PMC4556133 DOI: 10.1016/j.cell.2015.07.048] [Citation(s) in RCA: 221] [Impact Index Per Article: 24.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2014] [Revised: 03/05/2015] [Accepted: 07/21/2015] [Indexed: 01/12/2023]
Abstract
Deciphering the impact of genetic variants on gene regulation is fundamental to understanding human disease. Although gene regulation often involves long-range interactions, it is unknown to what extent non-coding genetic variants influence distal molecular phenotypes. Here, we integrate chromatin profiling for three histone marks in lymphoblastoid cell lines (LCLs) from 75 sequenced individuals with LCL-specific Hi-C and ChIA-PET-based chromatin contact maps to uncover one of the largest collections of local and distal histone quantitative trait loci (hQTLs). Distal QTLs are enriched within topologically associated domains and exhibit largely concordant variation of chromatin state coordinated by proximal and distal non-coding genetic variants. Histone QTLs are enriched for common variants associated with autoimmune diseases and enable identification of putative target genes of disease-associated variants from genome-wide association studies. These analyses provide insights into how genetic variation can affect human disease phenotypes by coordinated changes in chromatin at interacting regulatory elements.
Collapse
Affiliation(s)
- Fabian Grubert
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Judith B Zaugg
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA; The European Molecular Biology Laboratory Heidelberg, 69117 Heidelberg, Germany
| | - Maya Kasowski
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Oana Ursu
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Damek V Spacek
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Alicia R Martin
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Peyton Greenside
- Biomedical Informatics Graduate Training Program, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Rohith Srivas
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Doug H Phanstiel
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Aleksandra Pekowska
- The European Molecular Biology Laboratory Heidelberg, 69117 Heidelberg, Germany
| | - Nastaran Heidari
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Ghia Euskirchen
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Wolfgang Huber
- The European Molecular Biology Laboratory Heidelberg, 69117 Heidelberg, Germany
| | - Jonathan K Pritchard
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA; Department of Biology, Stanford University, Stanford, CA 94305, USA; Howard Hughes Medical Institute, Stanford University, Stanford, CA 94305, USA
| | - Carlos D Bustamante
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Lars M Steinmetz
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA; The European Molecular Biology Laboratory Heidelberg, 69117 Heidelberg, Germany
| | - Anshul Kundaje
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA; Department of Computer Science, Stanford University, Stanford, CA 94305, USA
| | - Michael Snyder
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA.
| |
Collapse
|
10
|
Klingel K, Sauter M, Ettischer N, Kandolf R, Ursu O. Heme oxygenase-1 mediates ROS production and ongoing injury in CVB3 myocarditis. Eur Heart J 2013. [DOI: 10.1093/eurheartj/eht309.p3867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
|
11
|
Gosline SJC, Spencer SJ, Ursu O, Fraenkel E. SAMNet: a network-based approach to integrate multi-dimensional high throughput datasets. Integr Biol (Camb) 2013; 4:1415-27. [PMID: 23060147 DOI: 10.1039/c2ib20072d] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The rapid development of high throughput biotechnologies has led to an onslaught of data describing genetic perturbations and changes in mRNA and protein levels in the cell. Because each assay provides a one-dimensional snapshot of active signaling pathways, it has become desirable to perform multiple assays (e.g. mRNA expression and phospho-proteomics) to measure a single condition. However, as experiments expand to accommodate various cellular conditions, proper analysis and interpretation of these data have become more challenging. Here we introduce a novel approach called SAMNet, for Simultaneous Analysis of Multiple Networks, that is able to interpret diverse assays over multiple perturbations. The algorithm uses a constrained optimization approach to integrate mRNA expression data with upstream genes, selecting edges in the protein-protein interaction network that best explain the changes across all perturbations. The result is a putative set of protein interactions that succinctly summarizes the results from all experiments, highlighting the network elements unique to each perturbation. We evaluated SAMNet in both yeast and human datasets. The yeast dataset measured the cellular response to seven different transition metals, and the human dataset measured cellular changes in four different lung cancer models of Epithelial-Mesenchymal Transition (EMT), a crucial process in tumor metastasis. SAMNet was able to identify canonical yeast metal-processing genes unique to each commodity in the yeast dataset, as well as human genes such as β-catenin and TCF7L2/TCF4 that are required for EMT signaling but escaped detection in the mRNA and phospho-proteomic data. Moreover, SAMNet also highlighted drugs likely to modulate EMT, identifying a series of less canonical genes known to be affected by the BCR-ABL inhibitor imatinib (Gleevec), suggesting a possible influence of this drug on EMT.
Collapse
Affiliation(s)
- Sara J C Gosline
- Dept. of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | | | | | | |
Collapse
|
12
|
Liu L, Ursu O, Ross J. Microtubule Motility in Crowded Conditions in vitro. Biophys J 2011. [DOI: 10.1016/j.bpj.2010.12.2654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022] Open
|
13
|
Bakkour S, Kolhatkar N, Karki S, Herman N, Ursu O, Raghavan V, Stranford S. Differential gene expression patterns associated with disease resistance in the MAIDS model (39.29). The Journal of Immunology 2010. [DOI: 10.4049/jimmunol.184.supp.39.29] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Abstract
Infection of susceptible strains of mice such as C57BL/6 (B6) with the LP-BM5 isolate of murine leukemia virus (MuLV) leads to an immunodeficiency syndrome similar to HIV-induced AIDS, known as MAIDS. However, other MAIDS-resistant strains such as BALB/c develop protective immunity against the virus after infection. Recently, we used DNA microarray analysis to compare lymphoid organ gene expression in these two mouse strains during the first week post infection. We identified a specific set of proteases and certain chemokines that are highly differentially expressed in BALB/c versus B6 mice after infection. This suggested to us that early changes in lymphoid organ architecture and/or cellular recruitment might help shape the early immune response to the virus. In follow-up experiments using real-time PCR and protein-based assays, we now show that neither protease expression nor lymphatic vascular architecture changes appear to be specifically associated with resistance to MAIDS. This despite the fact that these animals very rapidly control the viral load, showing significantly reduced viremia by five days after infection. In B6 mice, the viral load exponentially rises over the same time-course, with strongly divergent patterns apparent by day 5. We do find that certain chemokines are upregulated in the MAIDS-resistant mice at approximately one week after infection. These findings suggest that differences in cell trafficking might influence disease in this AIDS model system.
Collapse
Affiliation(s)
- Sonia Bakkour
- 1Biological Sciences, Mount Holyoke College, South Hadley, MA
| | | | - Sophiya Karki
- 1Biological Sciences, Mount Holyoke College, South Hadley, MA
| | - Nicole Herman
- 1Biological Sciences, Mount Holyoke College, South Hadley, MA
| | - Oana Ursu
- 1Biological Sciences, Mount Holyoke College, South Hadley, MA
| | - Vidya Raghavan
- 1Biological Sciences, Mount Holyoke College, South Hadley, MA
| | | |
Collapse
|