26
|
Sarikhani M, Garbern JC, Ma S, Sereda R, Conde J, Krähenbühl G, Escalante GO, Ahmed A, Buenrostro JD, Lee RT. Sustained Activation of AMPK Enhances Differentiation of Human iPSC-Derived Cardiomyocytes via Sirtuin Activation. Stem Cell Reports 2020; 15:498-514. [PMID: 32649901 PMCID: PMC7419706 DOI: 10.1016/j.stemcr.2020.06.012] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Revised: 06/11/2020] [Accepted: 06/12/2020] [Indexed: 12/18/2022] Open
Abstract
Recent studies suggest that metabolic regulation may improve differentiation of cardiomyocytes derived from induced pluripotent stem cells (iPSCs). AMP-activated protein kinase (AMPK) is a master regulator of metabolic activities. We investigated whether AMPK participates in iPSC-derived cardiomyocyte differentiation. We observed that AMPK phosphorylation at Thr172 increased at day 9 but then decreased after day 11 of differentiation to cardiomyocytes. Inhibition of AMPK with compound C significantly reduced mRNA and protein expression of cardiac troponins TNNT2 and TNNI3. Moreover, sustained AMPK activation using AICAR from days 9 to 14 of differentiation increased mRNA and protein expression of both TNNT2 and TNNI3. AICAR decreased acetylation of histone 3 at Lys9 and 56 and histone 4 at Lys16 (known target sites for nuclear-localized sirtuins [SIRT1, SIRT6]), suggesting that AMPK activation enhances sirtuin activity. Sustained AMPK activation during days 9–14 of differentiation induces sirtuin-mediated histone deacetylation and may enhance cardiomyocyte differentiation from iPSCs. iPSC-derived cardiomyocytes transiently increased AMPK phosphorylation at Thr172 Chemical inhibition of AMPK with compound C decreased TNNI3 and TNNT2 expression Sustained activation of AMPK using AICAR increased expression of TNNT2 and TNNI3 AICAR decreased acetylation of histones H3 (at Lys9 and Lys56) and H4 (at Lys16).
Collapse
|
27
|
LaFave LM, Kartha V, Ma S, Meli K, Priore ID, Lareau C, Sanker V, Naranjo S, Westcott P, Chiang Z, Brack A, Law T, Regev A, Buenrostro JD, Jacks T. Abstract PR08: Leveraging single-cell epigenomics to uncover regulatory programs in lung adenocarcinoma. Cancer Res 2020. [DOI: 10.1158/1538-7445.camodels2020-pr08] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Abstract
A hallmark of cancer development is the loss of cellular identity due to the dysregulation of regulatory networks that maintain functional, differentiated states. Chromatin state has been linked to this control of cellular identity and developmental processes; however, the mechanisms by which these regulatory landscapes are disrupted during cancer progression are not well understood. To this end, we leveraged an optimized single-cell ATAC-sequencing (scATAC-seq) strategy using a combinatorial indexing-based technology to assess chromatin state changes that occur in a mouse model of lung adenocarcinoma (LUAD). We utilized a well-established LUAD mouse model in which mice conditionally express the KrasG12D mutation and lose p53 expression in alveolar type II cells following Cre-mediated recombination, termed the KP model. In this model, KP tumors progress from early-stage hyperplasias to metastatic, advanced adenocarcinomas without acquiring additional somatic mutations, suggesting that these phenotypic transitions may be driven in part by epigenetic mechanisms. Using scATAC-seq, we profiled 16,044 normal lung, primary KP tumor, and metastatic cells to assess chromatin state changes across LUAD progression. We found that LUAD tumor evolution was characterized by a heterogeneous epigenomic continuum that was associated with loss of alveolar identity and progression to metastasis. Interestingly, we identified a reproducible metastatic state that was homogeneous, suggesting that cancer cells funnel towards a stable epigenetic state. In addition, we utilized novel computational tools to assess combinatorial transcription factor (TF)-driven programs that delineate stages of cancer progression in the KP model. From these analyses, we identified a clear transition point that corresponded to progressive gain of RUNX2 transcription factor activity. Using CRISPR activation and knockout strategies, we found that RUNX2 regulated the chromatin accessibility of several regulatory programs involved in extracellular matrix remodeling, suggesting that RUNX2 signaling is a critical event in LUAD progression and metastasis. In addition, gene signatures derived from this RUNX-activated state were highly predictive of survival in human LUAD. Together, these results demonstrated the power of single-cell epigenomics to identify TF-driven regulatory programs as key biomarkers of tumor progression.
This abstract is also being presented as Poster A40.
Citation Format: Lindsay M. LaFave, Vinay Kartha, Sai Ma, Kevin Meli, Isabella Del Priore, Caleb Lareau, Venkat Sanker, Santiago Naranjo, Peter Westcott, Zachary Chiang, Alison Brack, Travis Law, Aviv Regev, Jason D Buenrostro, Tyler Jacks. Leveraging single-cell epigenomics to uncover regulatory programs in lung adenocarcinoma [abstract]. In: Proceedings of the AACR Special Conference on the Evolving Landscape of Cancer Modeling; 2020 Mar 2-5; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2020;80(11 Suppl):Abstract nr PR08.
Collapse
|
28
|
Lareau CA, Ma S, Duarte FM, Buenrostro JD. Inference and effects of barcode multiplets in droplet-based single-cell assays. Nat Commun 2020; 11:866. [PMID: 32054859 PMCID: PMC7018801 DOI: 10.1038/s41467-020-14667-5] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2019] [Accepted: 01/23/2020] [Indexed: 12/12/2022] Open
Abstract
A widespread assumption for single-cell analyses specifies that one cell's nucleic acids are predominantly captured by one oligonucleotide barcode. Here, we show that ~13-21% of cell barcodes from the 10x Chromium scATAC-seq assay may have been derived from a droplet with more than one oligonucleotide sequence, which we call "barcode multiplets". We demonstrate that barcode multiplets can be derived from at least two different sources. First, we confirm that approximately 4% of droplets from the 10x platform may contain multiple beads. Additionally, we find that approximately 5% of beads may contain detectable levels of multiple oligonucleotide barcodes. We show that this artifact can confound single-cell analyses, including the interpretation of clonal diversity and proliferation of intra-tumor lymphocytes. Overall, our work provides a conceptual and computational framework to identify and assess the impacts of barcode multiplets in single-cell data.
Collapse
|
29
|
Zhang B, Ma S, Rachmin I, He M, Baral P, Choi S, Gonçalves WA, Shwartz Y, Fast EM, Su Y, Zon LI, Regev A, Buenrostro JD, Cunha TM, Chiu IM, Fisher DE, Hsu YC. Hyperactivation of sympathetic nerves drives depletion of melanocyte stem cells. Nature 2020; 577:676-681. [PMID: 31969699 PMCID: PMC7184936 DOI: 10.1038/s41586-020-1935-3] [Citation(s) in RCA: 149] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2019] [Accepted: 12/13/2019] [Indexed: 01/24/2023]
Abstract
Empirical and anecdotal evidence have associated stress with accelerated hair greying (formation of unpigmented hairs)1,2, but the scientific evidence linking the two is scant. Here, we report that acute stress leads to hair greying through fast depletion of melanocyte stem cells (MeSCs). Combining adrenalectomy, denervation, chemogenetics3,4, cell ablation, and MeSC-specific adrenergic receptor knockout, we found that stress-induced MeSC loss is independent of immune attack or adrenal stress hormones. Rather, hair greying results from activation of the sympathetic nerves that innervate the MeSC niche. Upon stress, sympathetic nerve activation leads to burst release of the neurotransmitter norepinephrine, which drives quiescent MeSCs into rapid proliferation, followed by differentiation, migration, and permanent depletion from the niche. Transient suppression of MeSC proliferation prevents stress-induced hair greying. Our studies demonstrate that acute stress-induced neuronal activity can drive rapid and permanent loss of somatic stem cells, and illustrate an example in which somatic stem cell maintenance is directly influenced by the overall physiological state of the organism.
Collapse
|
30
|
Chen H, Lareau C, Andreani T, Vinyard ME, Garcia SP, Clement K, Andrade-Navarro MA, Buenrostro JD, Pinello L. Assessment of computational methods for the analysis of single-cell ATAC-seq data. Genome Biol 2019; 20:241. [PMID: 31739806 PMCID: PMC6859644 DOI: 10.1186/s13059-019-1854-5] [Citation(s) in RCA: 157] [Impact Index Per Article: 31.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 10/03/2019] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Recent innovations in single-cell Assay for Transposase Accessible Chromatin using sequencing (scATAC-seq) enable profiling of the epigenetic landscape of thousands of individual cells. scATAC-seq data analysis presents unique methodological challenges. scATAC-seq experiments sample DNA, which, due to low copy numbers (diploid in humans), lead to inherent data sparsity (1-10% of peaks detected per cell) compared to transcriptomic (scRNA-seq) data (10-45% of expressed genes detected per cell). Such challenges in data generation emphasize the need for informative features to assess cell heterogeneity at the chromatin level. RESULTS We present a benchmarking framework that is applied to 10 computational methods for scATAC-seq on 13 synthetic and real datasets from different assays, profiling cell types from diverse tissues and organisms. Methods for processing and featurizing scATAC-seq data were compared by their ability to discriminate cell types when combined with common unsupervised clustering approaches. We rank evaluated methods and discuss computational challenges associated with scATAC-seq analysis including inherently sparse data, determination of features, peak calling, the effects of sequencing coverage and noise, and clustering performance. Running times and memory requirements are also discussed. CONCLUSIONS This reference summary of scATAC-seq methods offers recommendations for best practices with consideration for both the non-expert user and the methods developer. Despite variation across methods and datasets, SnapATAC, Cusanovich2018, and cisTopic outperform other methods in separating cell populations of different coverages and noise levels in both synthetic and real datasets. Notably, SnapATAC is the only method able to analyze a large dataset (> 80,000 cells).
Collapse
|
31
|
Chen H, Albergante L, Hsu JY, Lareau CA, Lo Bosco G, Guan J, Zhou S, Gorban AN, Bauer DE, Aryee MJ, Langenau DM, Zinovyev A, Buenrostro JD, Yuan GC, Pinello L. Single-cell trajectories reconstruction, exploration and mapping of omics data with STREAM. Nat Commun 2019; 10:1903. [PMID: 31015418 PMCID: PMC6478907 DOI: 10.1038/s41467-019-09670-4] [Citation(s) in RCA: 126] [Impact Index Per Article: 25.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2018] [Accepted: 03/22/2019] [Indexed: 12/14/2022] Open
Abstract
Single-cell transcriptomic assays have enabled the de novo reconstruction of lineage differentiation trajectories, along with the characterization of cellular heterogeneity and state transitions. Several methods have been developed for reconstructing developmental trajectories from single-cell transcriptomic data, but efforts on analyzing single-cell epigenomic data and on trajectory visualization remain limited. Here we present STREAM, an interactive pipeline capable of disentangling and visualizing complex branching trajectories from both single-cell transcriptomic and epigenomic data. We have tested STREAM on several synthetic and real datasets generated with different single-cell technologies. We further demonstrate its utility for understanding myoblast differentiation and disentangling known heterogeneity in hematopoiesis for different organisms. STREAM is an open-source software package.
Collapse
|
32
|
Ulirsch JC, Lareau CA, Bao EL, Ludwig LS, Guo MH, Benner C, Satpathy AT, Kartha VK, Salem RM, Hirschhorn JN, Finucane HK, Aryee MJ, Buenrostro JD, Sankaran VG. Interrogation of human hematopoiesis at single-cell and single-variant resolution. Nat Genet 2019; 51:683-693. [PMID: 30858613 PMCID: PMC6441389 DOI: 10.1038/s41588-019-0362-6] [Citation(s) in RCA: 109] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2018] [Accepted: 01/28/2019] [Indexed: 11/16/2022]
Abstract
Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here, we perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding for proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations. Fine mapping of blood cell traits in UK Biobank identifies putative causal variants and enrichment of fine-mapped variants in accessible chromatin of hematopoietic progenitor cells. The study provides an analytical framework for single-variant and single-cell analyses of genetic associations.
Collapse
|
33
|
Yoshida H, Lareau CA, Ramirez RN, Rose SA, Maier B, Wroblewska A, Desland F, Chudnovskiy A, Mortha A, Dominguez C, Tellier J, Kim E, Dwyer D, Shinton S, Nabekura T, Qi Y, Yu B, Robinette M, Kim KW, Wagers A, Rhoads A, Nutt SL, Brown BD, Mostafavi S, Buenrostro JD, Benoist C. The cis-Regulatory Atlas of the Mouse Immune System. Cell 2019; 176:897-912.e20. [PMID: 30686579 PMCID: PMC6785993 DOI: 10.1016/j.cell.2018.12.036] [Citation(s) in RCA: 236] [Impact Index Per Article: 47.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Revised: 10/26/2018] [Accepted: 12/20/2018] [Indexed: 01/08/2023]
Abstract
A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.
Collapse
|
34
|
Mor N, Rais Y, Sheban D, Peles S, Aguilera-Castrejon A, Zviran A, Elinger D, Viukov S, Geula S, Krupalnik V, Zerbib M, Chomsky E, Lasman L, Shani T, Bayerl J, Gafni O, Hanna S, Buenrostro JD, Hagai T, Masika H, Vainorius G, Bergman Y, Greenleaf WJ, Esteban MA, Elling U, Levin Y, Massarwa R, Merbl Y, Novershtern N, Hanna JH. Neutralizing Gatad2a-Chd4-Mbd3/NuRD Complex Facilitates Deterministic Induction of Naive Pluripotency. Cell Stem Cell 2018; 23:412-425.e10. [PMID: 30122475 DOI: 10.1016/j.stem.2018.07.004] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Revised: 06/19/2018] [Accepted: 07/10/2018] [Indexed: 01/08/2023]
Abstract
Mbd3, a member of nucleosome remodeling and deacetylase (NuRD) co-repressor complex, was previously identified as an inhibitor for deterministic induced pluripotent stem cell (iPSC) reprogramming, where up to 100% of donor cells successfully complete the process. NuRD can assume multiple mutually exclusive conformations, and it remains unclear whether this deterministic phenotype can be attributed to a specific Mbd3/NuRD subcomplex. Moreover, since complete ablation of Mbd3 blocks somatic cell proliferation, we aimed to explore functionally relevant alternative ways to neutralize Mbd3-dependent NuRD activity. We identify Gatad2a, a NuRD-specific subunit, whose complete deletion specifically disrupts Mbd3/NuRD repressive activity on the pluripotency circuitry during iPSC differentiation and reprogramming without ablating somatic cell proliferation. Inhibition of Gatad2a facilitates deterministic murine iPSC reprogramming within 8 days. We validate a distinct molecular axis, Gatad2a-Chd4-Mbd3, within Mbd3/NuRD as being critical for blocking reestablishment of naive pluripotency and further highlight signaling-dependent and post-translational modifications of Mbd3/NuRD that influence its interactions and assembly.
Collapse
|
35
|
Buenrostro JD, Corces MR, Lareau CA, Wu B, Schep AN, Aryee MJ, Majeti R, Chang HY, Greenleaf WJ. Integrated Single-Cell Analysis Maps the Continuous Regulatory Landscape of Human Hematopoietic Differentiation. Cell 2018; 173:1535-1548.e16. [PMID: 29706549 PMCID: PMC5989727 DOI: 10.1016/j.cell.2018.03.074] [Citation(s) in RCA: 378] [Impact Index Per Article: 63.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2017] [Revised: 01/03/2018] [Accepted: 03/27/2018] [Indexed: 12/26/2022]
Abstract
Human hematopoiesis involves cellular differentiation of multipotent cells into progressively more lineage-restricted states. While the chromatin accessibility landscape of this process has been explored in defined populations, single-cell regulatory variation has been hidden by ensemble averaging. We collected single-cell chromatin accessibility profiles across 10 populations of immunophenotypically defined human hematopoietic cell types and constructed a chromatin accessibility landscape of human hematopoiesis to characterize differentiation trajectories. We find variation consistent with lineage bias toward different developmental branches in multipotent cell types. We observe heterogeneity within common myeloid progenitors (CMPs) and granulocyte-macrophage progenitors (GMPs) and develop a strategy to partition GMPs along their differentiation trajectory. Furthermore, we integrated single-cell RNA sequencing (scRNA-seq) data to associate transcription factors to chromatin accessibility changes and regulatory elements to target genes through correlations of expression and regulatory element accessibility. Overall, this work provides a framework for integrative exploration of complex regulatory dynamics in a primary human tissue at single-cell resolution.
Collapse
|
36
|
Satpathy AT, Saligrama N, Buenrostro JD, Wei Y, Wu B, Rubin AJ, Granja JM, Lareau CA, Li R, Qi Y, Parker KR, Mumbach MR, Serratelli WS, Gennert DG, Schep AN, Corces MR, Khodadoust MS, Kim YH, Khavari PA, Greenleaf WJ, Davis MM, Chang HY. Transcript-indexed ATAC-seq for precision immune profiling. Nat Med 2018; 24:580-590. [PMID: 29686426 PMCID: PMC5948148 DOI: 10.1038/s41591-018-0008-8] [Citation(s) in RCA: 103] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2017] [Accepted: 02/07/2018] [Indexed: 02/04/2023]
Abstract
T cells create vast amounts of diversity in the genes that encode their T cell receptors (TCRs), which enables individual clones to recognize specific peptide-major histocompatibility complex (MHC) ligands. Here we combined sequencing of the TCR-encoding genes with assay for transposase-accessible chromatin with sequencing (ATAC-seq) analysis at the single-cell level to provide information on the TCR specificity and epigenomic state of individual T cells. By using this approach, termed transcript-indexed ATAC-seq (T-ATAC-seq), we identified epigenomic signatures in immortalized leukemic T cells, primary human T cells from healthy volunteers and primary leukemic T cells from patient samples. In peripheral blood CD4+ T cells from healthy individuals, we identified cis and trans regulators of naive and memory T cell states and found substantial heterogeneity in surface-marker-defined T cell populations. In patients with a leukemic form of cutaneous T cell lymphoma, T-ATAC-seq enabled identification of leukemic and nonleukemic regulatory pathways in T cells from the same individual by allowing separation of the signals that arose from the malignant clone from the background T cell noise. Thus, T-ATAC-seq is a new tool that enables analysis of epigenomic landscapes in clonal T cells and should be valuable for studies of T cell malignancy, immunity and immunotherapy.
Collapse
|
37
|
Finucane HK, Reshef YA, Anttila V, Slowikowski K, Gusev A, Byrnes A, Gazal S, Loh PR, Lareau C, Shoresh N, Genovese G, Saunders A, Macosko E, Pollack S, Perry JRB, Buenrostro JD, Bernstein BE, Raychaudhuri S, McCarroll S, Neale BM, Price AL. Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types. Nat Genet 2018; 50:621-629. [PMID: 29632380 PMCID: PMC5896795 DOI: 10.1038/s41588-018-0081-4] [Citation(s) in RCA: 561] [Impact Index Per Article: 93.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 01/29/2018] [Indexed: 02/07/2023]
Abstract
We introduce an approach to identify disease-relevant tissues and cell types by analyzing gene expression data together with genome-wide association study (GWAS) summary statistics. Our approach uses stratified linkage disequilibrium (LD) score regression to test whether disease heritability is enriched in regions surrounding genes with the highest specific expression in a given tissue. We applied our approach to gene expression data from several sources together with GWAS summary statistics for 48 diseases and traits (average N = 169,331) and found significant tissue-specific enrichments (false discovery rate (FDR) < 5%) for 34 traits. In our analysis of multiple tissues, we detected a broad range of enrichments that recapitulated known biology. In our brain-specific analysis, significant enrichments included an enrichment of inhibitory over excitatory neurons for bipolar disorder, and excitatory over inhibitory neurons for schizophrenia and body mass index. Our results demonstrate that our polygenic approach is a powerful way to leverage gene expression data for interpreting GWAS signals.
Collapse
|
38
|
Koh AS, Miller EL, Buenrostro JD, Moskowitz DM, Wang J, Greenleaf WJ, Chang HY, Crabtree GR. Rapid chromatin repression by Aire provides precise control of immune tolerance. Nat Immunol 2018; 19:162-172. [PMID: 29335648 PMCID: PMC6049828 DOI: 10.1038/s41590-017-0032-8] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Accepted: 12/07/2017] [Indexed: 01/23/2023]
Abstract
Aire mediates the expression of tissue-specific antigens in thymic epithelial cells to promote tolerance against self-reactive T lymphocytes. However, the mechanism that allows expression of tissue-specific genes at levels that prevent harm is unknown. Here we show that Brg1 generates accessibility at tissue-specific loci to impose central tolerance. We found that Aire has an intrinsic repressive function that restricts chromatin accessibility and opposes Brg1 across the genome. Aire exerted this repressive influence within minutes after recruitment to chromatin and restrained the amplitude of active transcription. Disease-causing mutations that impair Aire-induced activation also impair the protein's repressive function, which indicates dual roles for Aire. Together, Brg1 and Aire fine-tune the expression of tissue-specific genes at levels that prevent toxicity yet promote immune tolerance.
Collapse
|
39
|
Daugherty AC, Yeo RW, Buenrostro JD, Greenleaf WJ, Kundaje A, Brunet A. Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans. Genome Res 2017. [PMID: 29141961 DOI: 10.1101/088732] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.
Collapse
|
40
|
Schep AN, Wu B, Buenrostro JD, Greenleaf WJ. chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data. Nat Methods 2017; 14:975-978. [PMID: 28825706 PMCID: PMC5623146 DOI: 10.1038/nmeth.4401] [Citation(s) in RCA: 705] [Impact Index Per Article: 100.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2017] [Accepted: 07/21/2017] [Indexed: 12/29/2022]
Abstract
Single cell ATAC-seq (scATAC) yields sparse data that makes application of conventional analysis approaches challenging. We developed chromVAR, an R package for analyzing sparse chromatin accessibility data by estimating gain or loss of accessibility within peaks sharing the same motif or annotation while controlling for technical biases. chromVAR enables accurate clustering of scATAC-seq profiles and enables characterization of known and de novo sequence motifs associated with variation in chromatin accessibility.
Collapse
|
41
|
Soto-Feliciano YM, Bartlebaugh JME, Liu Y, Sánchez-Rivera FJ, Bhutkar A, Weintraub AS, Buenrostro JD, Cheng CS, Regev A, Jacks TE, Young RA, Hemann MT. PHF6 regulates phenotypic plasticity through chromatin organization within lineage-specific genes. Genes Dev 2017; 31:973-989. [PMID: 28607179 PMCID: PMC5495126 DOI: 10.1101/gad.295857.117] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Accepted: 05/15/2017] [Indexed: 12/17/2022]
Abstract
In this study, Soto-Feliciano et al. describe the function of the plant homeodomain finger 6 (PHF6) protein in leukemia and define its role in regulating chromatin accessibility to lineage-specific transcription factors. Their findings suggest that active maintenance of a precise chromatin landscape is essential for sustaining proper leukemia cell identity and that loss of a single factor (PHF6) can cause focal changes in chromatin accessibility and nucleosome positioning that render cells susceptible to lineage transition. Developmental and lineage plasticity have been observed in numerous malignancies and have been correlated with tumor progression and drug resistance. However, little is known about the molecular mechanisms that enable such plasticity to occur. Here, we describe the function of the plant homeodomain finger protein 6 (PHF6) in leukemia and define its role in regulating chromatin accessibility to lineage-specific transcription factors. We show that loss of Phf6 in B-cell leukemia results in systematic changes in gene expression via alteration of the chromatin landscape at the transcriptional start sites of B-cell- and T-cell-specific factors. Additionally, Phf6KO cells show significant down-regulation of genes involved in the development and function of normal B cells, show up-regulation of genes involved in T-cell signaling, and give rise to mixed-lineage lymphoma in vivo. Engagement of divergent transcriptional programs results in phenotypic plasticity that leads to altered disease presentation in vivo, tolerance of aberrant oncogenic signaling, and differential sensitivity to frontline and targeted therapies. These findings suggest that active maintenance of a precise chromatin landscape is essential for sustaining proper leukemia cell identity and that loss of a single factor (PHF6) can cause focal changes in chromatin accessibility and nucleosome positioning that render cells susceptible to lineage transition.
Collapse
|
42
|
Miller EL, Hargreaves DC, Kadoch C, Chang CY, Calarco JP, Hodges C, Buenrostro JD, Cui K, Greenleaf WJ, Zhao K, Crabtree GR. TOP2 synergizes with BAF chromatin remodeling for both resolution and formation of facultative heterochromatin. Nat Struct Mol Biol 2017; 24:344-352. [PMID: 28250416 PMCID: PMC5395302 DOI: 10.1038/nsmb.3384] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2016] [Accepted: 02/01/2017] [Indexed: 12/19/2022]
Abstract
Resolution and formation of facultative heterochromatin is essential to development, reprogramming, and oncogenesis. The mechanisms underlying these changes are poorly understood due to the inability to study heterochromatin dynamics and structure in vivo. We devised an in vivo approach to investigate these mechanisms and found that topoisomerase II (TOP2), but not TOP1, synergizes with BAF (mSWI/SNF) ATP-dependent chromatin remodeling complexes genome-wide to resolve facultative heterochromatin to accessible chromatin independent of transcription, indicating that changes in DNA topology through (de-)catenation rather than release of torsional stress through swiveling is necessary for heterochromatin resolution. In turn, TOP2 and BAF cooperate to recruit pluripotency factors, explaining some of the instructive roles of BAF complexes. Unexpectedly, we found that TOP2, also plays a role in the reformation of facultative heterochromatin, suggesting that facultative heterochromatin and accessible chromatin exist at different states of catenation or other topologies, which may be critical to their structures.
Collapse
|
43
|
Moskowitz DM, Zhang DW, Hu B, Le Saux S, Yanes RE, Ye Z, Buenrostro JD, Weyand CM, Greenleaf WJ, Goronzy JJ. Epigenomics of human CD8 T cell differentiation and aging. Sci Immunol 2017; 2:eaag0192. [PMID: 28439570 PMCID: PMC5399889 DOI: 10.1126/sciimmunol.aag0192] [Citation(s) in RCA: 147] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The efficacy of the adaptive immune response declines dramatically with age, but the cell-intrinsic mechanisms driving immune aging in humans remain poorly understood. Immune aging is characterized by a loss of self-renewing naïve cells and the accumulation of differentiated but dysfunctional cells within the CD8 T cell compartment. Using ATAC-seq, we inferred the transcription factor binding activities correlated with naive and central and effector memory CD8 T cell states in young adults. Integrating our results with RNA-seq, we identified transcription networks associated with CD8 T cell differentiation, with prominent roles implicated for BATF, ETS1, Eomes, and Sp1. Extending our analysis to aged humans, we found that the differences between the memory and naive subsets were largely preserved across age, but that naive and central memory cells from older individuals exhibited a shift toward more differentiated patterns of chromatin openness. Additionally, aged naive cells displayed a loss in chromatin accessibility at gene promoters, largely associated with a decrease in NRF1 binding. This shift was implicated in a marked drop-off in the ability of the aged naive cells to transcribe respiratory chain genes, which may explain the reduced capacity of oxidative phosphorylation in older naïve cells. Our findings identify BATF- and NRF1-driven gene regulation as potential targets for delaying CD8 T cell aging and restoring function.
Collapse
|
44
|
Litzenburger UM, Buenrostro JD, Wu B, Shen Y, Sheffield NC, Kathiria A, Greenleaf WJ, Chang HY. Single-cell epigenomic variability reveals functional cancer heterogeneity. Genome Biol 2017; 18:15. [PMID: 28118844 PMCID: PMC5259890 DOI: 10.1186/s13059-016-1133-7] [Citation(s) in RCA: 72] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2016] [Accepted: 12/12/2016] [Indexed: 12/27/2022] Open
Abstract
Background Cell-to-cell heterogeneity is a major driver of cancer evolution, progression, and emergence of drug resistance. Epigenomic variation at the single-cell level can rapidly create cancer heterogeneity but is difficult to detect and assess functionally. Results We develop a strategy to bridge the gap between measurement and function in single-cell epigenomics. Using single-cell chromatin accessibility and RNA-seq data in K562 leukemic cells, we identify the cell surface marker CD24 as co-varying with chromatin accessibility changes linked to GATA transcription factors in single cells. Fluorescence-activated cell sorting of CD24 high versus low cells prospectively isolated GATA1 and GATA2 high versus low cells. GATA high versus low cells express differential gene regulatory networks, differential sensitivity to the drug imatinib mesylate, and differential self-renewal capacity. Lineage tracing experiments show that GATA/CD24hi cells have the capability to rapidly reconstitute the heterogeneity within the entire starting population, suggesting that GATA expression levels drive a phenotypically relevant source of epigenomic plasticity. Conclusion Single-cell chromatin accessibility can guide prospective characterization of cancer heterogeneity. Epigenomic subpopulations in cancer impact drug sensitivity and the clonal dynamics of cancer evolution. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-1133-7) contains supplementary material, which is available to authorized users.
Collapse
|
45
|
Chen X, Shen Y, Draper W, Buenrostro JD, Litzenburger U, Cho SW, Satpathy AT, Carter AC, Ghosh RP, East-Seletsky A, Doudna JA, Greenleaf WJ, Liphardt JT, Chang HY. ATAC-see reveals the accessible genome by transposase-mediated imaging and sequencing. Nat Methods 2016; 13:1013-1020. [PMID: 27749837 DOI: 10.1038/nmeth.4031] [Citation(s) in RCA: 146] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 09/19/2016] [Indexed: 01/10/2023]
Abstract
Spatial organization of the genome plays a central role in gene expression, DNA replication, and repair. But current epigenomic approaches largely map DNA regulatory elements outside of the native context of the nucleus. Here we report assay of transposase-accessible chromatin with visualization (ATAC-see), a transposase-mediated imaging technology that employs direct imaging of the accessible genome in situ, cell sorting, and deep sequencing to reveal the identity of the imaged elements. ATAC-see revealed the cell-type-specific spatial organization of the accessible genome and the coordinated process of neutrophil chromatin extrusion, termed NETosis. Integration of ATAC-see with flow cytometry enables automated quantitation and prospective cell isolation as a function of chromatin accessibility, and it reveals a cell-cycle dependence of chromatin accessibility that is especially dynamic in G1 phase. The integration of imaging and epigenomics provides a general and scalable approach for deciphering the spatiotemporal architecture of gene control.
Collapse
|
46
|
Corces MR, Buenrostro JD, Wu B, Greenside PG, Chan SM, Koenig JL, Snyder MP, Pritchard JK, Kundaje A, Greenleaf WJ, Majeti R, Chang HY. Lineage-specific and single-cell chromatin accessibility charts human hematopoiesis and leukemia evolution. Nat Genet 2016; 48:1193-203. [PMID: 27526324 PMCID: PMC5042844 DOI: 10.1038/ng.3646] [Citation(s) in RCA: 714] [Impact Index Per Article: 89.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2015] [Accepted: 07/18/2016] [Indexed: 02/07/2023]
Abstract
We define the chromatin accessibility and transcriptional landscapes in 13 human primary blood cell types that span the hematopoietic hierarchy. Exploiting the finding that the enhancer landscape better reflects cell identity than mRNA levels, we enable 'enhancer cytometry' for enumeration of pure cell types from complex populations. We identify regulators governing hematopoietic differentiation and further show the lineage ontogeny of genetic elements linked to diverse human diseases. In acute myeloid leukemia (AML), chromatin accessibility uncovers unique regulatory evolution in cancer cells with a progressively increasing mutation burden. Single AML cells exhibit distinctive mixed regulome profiles corresponding to disparate developmental stages. A method to account for this regulatory heterogeneity identified cancer-specific deviations and implicated HOX factors as key regulators of preleukemic hematopoietic stem cell characteristics. Thus, regulome dynamics can provide diverse insights into hematopoietic development and disease.
Collapse
|
47
|
Mazumdar C, Shen Y, Xavy S, Zhao F, Reinisch A, Li R, Corces MR, Flynn RA, Buenrostro JD, Chan SM, Thomas D, Koenig JL, Hong WJ, Chang HY, Majeti R. Leukemia-Associated Cohesin Mutants Dominantly Enforce Stem Cell Programs and Impair Human Hematopoietic Progenitor Differentiation. Cell Stem Cell 2015; 17:675-688. [PMID: 26607380 PMCID: PMC4671831 DOI: 10.1016/j.stem.2015.09.017] [Citation(s) in RCA: 148] [Impact Index Per Article: 16.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Revised: 08/24/2015] [Accepted: 09/23/2015] [Indexed: 01/14/2023]
Abstract
Recurrent mutations in cohesin complex proteins have been identified in pre-leukemic hematopoietic stem cells and during the early development of acute myeloid leukemia and other myeloid malignancies. Although cohesins are involved in chromosome separation and DNA damage repair, cohesin complex functions during hematopoiesis and leukemic development are unclear. Here, we show that mutant cohesin proteins block differentiation of human hematopoietic stem and progenitor cells (HSPCs) in vitro and in vivo and enforce stem cell programs. These effects are restricted to immature HSPC populations, where cohesin mutants show increased chromatin accessibility and likelihood of transcription factor binding site occupancy by HSPC regulators including ERG, GATA2, and RUNX1, as measured by ATAC-seq and ChIP-seq. Epistasis experiments show that silencing these transcription factors rescues the differentiation block caused by cohesin mutants. Together, these results show that mutant cohesins impair HSPC differentiation by controlling chromatin accessibility and transcription factor activity, possibly contributing to leukemic disease.
Collapse
|
48
|
Maza I, Caspi I, Zviran A, Chomsky E, Rais Y, Viukov S, Geula S, Buenrostro JD, Weinberger L, Krupalnik V, Hanna S, Zerbib M, Dutton JR, Greenleaf WJ, Massarwa R, Novershtern N, Hanna JH. Transient acquisition of pluripotency during somatic cell transdifferentiation with iPSC reprogramming factors. Nat Biotechnol 2015; 33:769-74. [PMID: 26098448 PMCID: PMC4500825 DOI: 10.1038/nbt.3270] [Citation(s) in RCA: 113] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2014] [Accepted: 06/01/2015] [Indexed: 01/20/2023]
Abstract
Somatic cells can be transdifferentiated to other cell types without passing through a pluripotent state by ectopic expression of appropriate transcription factors1,2. Recent reports have proposed an alternative transdifferentiation method in which fibroblasts are directly converted to various mature somatic cell types by brief expression of the induced pluripotent stem cell (iPSC) reprogramming factors Oct4, Sox2, Klf4 and c-Myc (OSKM) followed by cell expansion in media that promote lineage differentiation3–6. Here we test this method using genetic lineage tracing for expression of endogenous Nanog and Oct4 and for X chromosome reactivation, as these events mark acquisition of pluripotency. We show that the vast majority of reprogrammed cardiomyocytes or neural stem cells obtained from mouse fibroblasts by OSKM-induced transdifferentiation pass through a transient pluripotent state, and that their derivation is molecularly coupled to iPSC formation mechanisms. Our findings underscore the importance of defining trajectories during cell reprogramming by different methods.
Collapse
|
49
|
Buenrostro JD, Wu B, Chang HY, Greenleaf WJ. ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide. ACTA ACUST UNITED AC 2015; 109:21.29.1-21.29.9. [PMID: 25559105 DOI: 10.1002/0471142727.mb2129s109] [Citation(s) in RCA: 1826] [Impact Index Per Article: 202.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
This unit describes Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq), a method for mapping chromatin accessibility genome-wide. This method probes DNA accessibility with hyperactive Tn5 transposase, which inserts sequencing adapters into accessible regions of chromatin. Sequencing reads can then be used to infer regions of increased accessibility, as well as to map regions of transcription-factor binding and nucleosome position. The method is a fast and sensitive alternative to DNase-seq for assaying chromatin accessibility genome-wide, or to MNase-seq for assaying nucleosome positions in accessible regions of the genome.
Collapse
|
50
|
Buenrostro JD, Chircus LM, Araya CL, Layton CJ, Chang HY, Snyder MP, Greenleaf WJ. Quantitative analysis of RNA-protein interactions on a massively parallel array reveals biophysical and evolutionary landscapes. Nat Biotechnol 2014; 32:562-8. [PMID: 24727714 PMCID: PMC4414031 DOI: 10.1038/nbt.2880] [Citation(s) in RCA: 155] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2014] [Accepted: 03/19/2014] [Indexed: 01/13/2023]
Abstract
RNA-protein interactions drive fundamental biological processes and are targets for molecular engineering, yet quantitative and comprehensive understanding of the sequence determinants of affinity remains limited. Here we repurpose a high-throughput sequencing instrument to quantitatively measure binding and dissociation of a fluorescently labeled protein to >10(7) RNA targets generated on a flow cell surface by in situ transcription and intermolecular tethering of RNA to DNA. Studying the MS2 coat protein, we decompose the binding energy contributions from primary and secondary RNA structure, and observe that differences in affinity are often driven by sequence-specific changes in both association and dissociation rates. By analyzing the biophysical constraints and modeling mutational paths describing the molecular evolution of MS2 from low- to high-affinity hairpins, we quantify widespread molecular epistasis and a long-hypothesized, structure-dependent preference for G:U base pairs over C:A intermediates in evolutionary trajectories. Our results suggest that quantitative analysis of RNA on a massively parallel array (RNA-MaP) provides generalizable insight into the biophysical basis and evolutionary consequences of sequence-function relationships.
Collapse
|