Cell of origin epigenetic priming determines susceptibility to Tet2 mutation

Hematopoietic stem cell (HSC) mutations can result in clonal hematopoiesis (CH) with heterogeneous clinical outcomes. Here, we investigate how the cell state preceding Tet2 mutation impacts the pre-malignant phenotype. Using an inducible system for clonal analysis of myeloid progenitors, we find that the epigenetic features of clones at similar differentiation status are highly heterogeneous and functionally respond differently to Tet2 mutation. Cell differentiation stage also influences Tet2 mutation response indicating that the cell of origin's epigenome modulates clone-specific behaviors in CH. Molecular features associated with higher risk outcomes include Sox4 that sensitizes cells to Tet2 inactivation, inducing dedifferentiation, altered metabolism and increasing the in vivo clonal output of mutant cells, as confirmed in primary GMP and HSC models. Our findings validate the hypothesis that epigenetic features can predispose specific clones for dominance, explaining why identical genetic mutations can result in different phenotypes.

Epigenetic memory of coronavirus infection in innate immune cells and their progenitors

Inflammation can trigger lasting phenotypes in immune and non-immune cells. Whether and how human infections and associated inflammation can form innate immune memory in hematopoietic stem and progenitor cells (HSPC) has remained unclear. We found that circulating HSPC, enriched from peripheral blood, captured the diversity of bone marrow HSPC, enabling investigation of their epigenomic reprogramming following coronavirus disease 2019 (COVID-19). Alterations in innate immune phenotypes and epigenetic programs of HSPC persisted for months to 1 year following severe COVID-19 and were associated with distinct transcription factor (TF) activities, altered regulation of inflammatory programs, and durable increases in myelopoiesis. HSPC epigenomic alterations were conveyed, through differentiation, to progeny innate immune cells. Early activity of IL-6 contributed to these persistent phenotypes in human COVID-19 and a mouse coronavirus infection model. Epigenetic reprogramming of HSPC may underlie altered immune function following infection and be broadly relevant, especially for millions of COVID-19 survivors.

Systematic benchmarking of single-cell ATAC-sequencing protocols

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Heritable transcriptional defects from aberrations of nuclear architecture

Transcriptional heterogeneity due to plasticity of the epigenetic state of chromatin contributes to tumour evolution, metastasis and drug resistance1-3. However, the mechanisms that cause this epigenetic variation are incompletely understood. Here we identify micronuclei and chromosome bridges, aberrations in the nucleus common in cancer4,5, as sources of heritable transcriptional suppression. Using a combination of approaches, including long-term live-cell imaging and same-cell single-cell RNA sequencing (Look-Seq2), we identified reductions in gene expression in chromosomes from micronuclei. With heterogeneous penetrance, these changes in gene expression can be heritable even after the chromosome from the micronucleus has been re-incorporated into a normal daughter cell nucleus. Concomitantly, micronuclear chromosomes acquire aberrant epigenetic chromatin marks. These defects may persist as variably reduced chromatin accessibility and reduced gene expression after clonal expansion from single cells. Persistent transcriptional repression is strongly associated with, and may be explained by, markedly long-lived DNA damage. Epigenetic alterations in transcription may therefore be inherently coupled to chromosomal instability and aberrations in nuclear architecture.

Positional influence on cellular transcriptional identity revealed through spatially segmented single-cell transcriptomics

Single-cell RNA sequencing (scRNA-seq) is a powerful technique for describing cell states. Identifying the spatial arrangement of these states in tissues remains challenging, with the existing methods requiring niche methodologies and expertise. Here, we describe segmentation by exogenous perfusion (SEEP), a rapid and integrated method to link surface proximity and environment accessibility to transcriptional identity within three-dimensional (3D) disease models. The method utilizes the steady-state diffusion kinetics of a fluorescent dye to establish a gradient along the radial axis of disease models. Classification of sample layers based on dye accessibility enables dissociated and sorted cells to be characterized by transcriptomic and regional identities. Using SEEP, we analyze spheroid, organoid, and in vivo tumor models of high-grade serous ovarian cancer (HGSOC). The results validate long-standing beliefs about the relationship between cell state and position while revealing new concepts regarding how spatially unique microenvironments influence the identity of individual cells within tumors.

Photoselective sequencing: microscopically guided genomic measurements with subcellular resolution

In biological systems, spatial organization and function are interconnected. Here we present photoselective sequencing, a new method for genomic and epigenomic profiling within morphologically distinct regions. Starting with an intact biological specimen, photoselective sequencing uses targeted illumination to selectively unblock a photocaged fragment library, restricting the sequencing-based readout to microscopically identified spatial regions. We validate photoselective sequencing by measuring the chromatin accessibility profiles of fluorescently labeled cell types within the mouse brain and comparing with published data. Furthermore, by combining photoselective sequencing with a computational strategy for decomposing bulk accessibility profiles, we find that the oligodendrocyte-lineage-cell population is relatively enriched for oligodendrocyte-progenitor cells in the cortex versus the corpus callosum. Finally, we leverage photoselective sequencing at the subcellular scale to identify features of chromatin that are correlated with positioning at the nuclear periphery. These results collectively demonstrate that photoselective sequencing is a flexible and generalizable platform for exploring the interplay of spatial structures with genomic and epigenomic properties.

Structural and functional properties of mSWI/SNF chromatin remodeling complexes revealed through single-cell perturbation screens

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.

Single-cell multi-scale footprinting reveals the modular organization of DNA regulatory elements

Cis-regulatory elements control gene expression and are dynamic in their structure, reflecting changes to the composition of diverse effector proteins over time1-3. Here we sought to connect the structural changes at cis-regulatory elements to alterations in cellular fate and function. To do this we developed PRINT, a computational method that uses deep learning to correct sequence bias in chromatin accessibility data and identifies multi-scale footprints of DNA-protein interactions. We find that multi-scale footprints enable more accurate inference of TF and nucleosome binding. Using PRINT with single-cell multi-omics, we discover wide-spread changes to the structure and function of candidate cis-regulatory elements (cCREs) across hematopoiesis, wherein nucleosomes slide, expose DNA for TF binding, and promote gene expression. Activity segmentation using the co-variance across cell states identifies "sub-cCREs" as modular cCRE subunits of regulatory DNA. We apply this single-cell and PRINT approach to characterize the age-associated alterations to cCREs within hematopoietic stem cells (HSCs). Remarkably, we find a spectrum of aging alterations among HSCs corresponding to a global gain of sub-cCRE activity while preserving cCRE accessibility. Collectively, we reveal the functional importance of cCRE structure across cell states, highlighting changes to gene regulation at single-cell and single-base-pair resolution.

A transcription factor atlas of directed differentiation

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

Functional inference of gene regulation using single-cell multi-omics

Cells require coordinated control over gene expression when responding to environmental stimuli. Here we apply scATAC-seq and single-cell RNA sequencing (scRNA-seq) in resting and stimulated human blood cells. Collectively, we generate ~91,000 single-cell profiles, allowing us to probe the cis-regulatory landscape of the immunological response across cell types, stimuli, and time. Advancing tools to integrate multi-omics data, we develop functional inference of gene regulation (FigR), a framework to computationally pair scA-TAC-seq with scRNA-seq cells, connect distal cis-regulatory elements to genes, and infer gene-regulatory networks (GRNs) to identify candidate transcription factor (TF) regulators. Utilizing these paired multi-omics data, we define domains of regulatory chromatin (DORCs) of immune stimulation and find that cells alter chromatin accessibility and gene expression at timescales of minutes. Construction of the stimulation GRN elucidates TF activity at disease-associated DORCs. Overall, FigR enables elucidation of regulatory interactions across single-cell data, providing new opportunities to understand the function of cells within tissues.

Single-cell methods for measuring chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) are rapidly evolving, but tools to integrate data and infer gene-regulatory relationships remain limited. Here we generate multi-omics data of resting and stimulated human blood cells and present a new computational framework for constructing gene-regulatory networks (GRNs). Specifically, we describe functional inference of gene regulation (FigR), a workflow to (1) pair scATAC-seq with scRNA-seq, (2) connect cis-regulatory elements to target genes, and (3) identify TF-gene relationships.

Single-Cell Epigenomics Reveals Mechanisms of Cancer Progression

Cancer initiation is driven by the cooperation between genetic and epigenetic aberrations that disrupt gene regulatory programs critical to maintaining specialized cellular functions. After initiation, cells acquire additional genetic and epigenetic alterations influenced by tumor-intrinsic and -extrinsic mechanisms, which increase intratumoral heterogeneity, reshape the cell's underlying gene regulatory networks and promote cancer evolution. Furthermore, environmental or therapeutic insults drive the selection of heterogeneous cell states, with implications for cancer initiation, maintenance, and drug resistance. The advancement of single-cell genomics has begun to uncover the full repertoire of chromatin and gene expression states (cell states) that exist within individual tumors. These single-cell analyses suggest that cells diversify in their regulatory states upon transformation by co-opting damage-induced and nonlineage regulatory programs that can lead to epigenomic plasticity. Here, we review these recent studies related to regulatory state changes in cancer progression and highlight the growing single-cell epigenomics toolkit poised to address unresolved questions in the field.

Spatial genomics enables multi-modal study of clonal heterogeneity in tissues

The state and behaviour of a cell can be influenced by both genetic and environmental factors. In particular, tumour progression is determined by underlying genetic aberrations^1-4 as well as the makeup of the tumour microenvironment^5,6. Quantifying the contributions of these factors requires new technologies that can accurately measure the spatial location of genomic sequence together with phenotypic readouts. Here we developed slide-DNA-seq, a method for capturing spatially resolved DNA sequences from intact tissue sections. We demonstrate that this method accurately preserves local tumour architecture and enables the de novo discovery of distinct tumour clones and their copy number alterations. We then apply slide-DNA-seq to a mouse model of metastasis and a primary human cancer, revealing that clonal populations are confined to distinct spatial regions. Moreover, through integration with spatial transcriptomics, we uncover distinct sets of genes that are associated with clone-specific genetic aberrations, the local tumour microenvironment, or both. Together, this multi-modal spatial genomics approach provides a versatile platform for quantifying how cell-intrinsic and cell-extrinsic factors contribute to gene expression, protein abundance and other cellular phenotypes.

SMARCA4 inactivation promotes lineage-specific transformation and early metastatic features in the lung

SMARCA4/BRG1 encodes for one of two mutually exclusive ATPases present in mammalian SWI/SNF chromatin remodeling complexes and is frequently mutated in human lung adenocarcinoma. However, the functional consequences of SMARCA4 mutation on tumor initiation, progression, and chromatin regulation in lung cancer remain poorly understood. Here, we demonstrate that loss of Smarca4 sensitizes CCSP+ cells within the lung in a cell-type dependent fashion to malignant transformation and tumor progression, resulting in highly advanced dedifferentiated tumors and increased metastatic incidence. Consistent with these phenotypes, Smarca4-deficient primary tumors lack lung lineage transcription factor activities and resemble a metastatic cell state. Mechanistically, we show that Smarca4 loss impairs the function of all three classes of SWI/SNF complexes, resulting in decreased chromatin accessibility at lung lineage motifs and ultimately accelerating tumor progression. Thus, we propose that the SWI/SNF complex - via Smarca4 - acts as a gatekeeper for lineage-specific cellular transformation and metastasis during lung cancer evolution.

A microRNA expression and regulatory element activity atlas of the mouse immune system

To better define the control of immune system regulation, we generated an atlas of microRNA (miRNA) expression from 63 mouse immune cell populations and connected these signatures with assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation followed by sequencing (ChIP-seq) and nascent RNA profiles to establish a map of miRNA promoter and enhancer usage in immune cells. miRNA complexity was relatively low, with >90% of the miRNA compartment of each population comprising <75 miRNAs; however, each cell type had a unique miRNA signature. Integration of miRNA expression with chromatin accessibility revealed putative regulatory elements for differentially expressed miRNAs, including miR-21a, miR-146a and miR-223. The integrated maps suggest that many miRNAs utilize multiple promoters to reach high abundance and identified dominant and divergent miRNA regulatory elements between lineages and during development that may be used by clustered miRNAs, such as miR-99a/let-7c/miR-125b, to achieve distinct expression. These studies, with web-accessible data, help delineate the cis-regulatory elements controlling miRNA signatures of the immune system.

Protocol for single-cell ATAC sequencing using combinatorial indexing in mouse lung adenocarcinoma

Single-cell ATAC sequencing using combinatorial indexing (sciATAC-seq) enables the identification of chromatin accessibility profiles at single-cell resolution with a dual-barcoding approach during transposition and library construction. Unlike commercial droplet-based approaches, sciATAC-seq is a cost-effective, extensible strategy that permits flexibility in the experimental scale via multiplexed barcoding across samples or perturbations. In contrast, droplet-based approaches have higher cell recovery and may be advantageous when cell input is limited. The step-by-step sciATAC-seq protocol described here is amenable to diverse cell types and fixed samples. For complete details on the use and execution of this protocol, please refer to LaFave et al. (2020).

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

Deep learning-based enhancement of epigenomics data with AtacWorks

ATAC-seq is a widely-applied assay used to measure genome-wide chromatin accessibility; however, its ability to detect active regulatory regions can depend on the depth of sequencing coverage and the signal-to-noise ratio. Here we introduce AtacWorks, a deep learning toolkit to denoise sequencing coverage and identify regulatory peaks at base-pair resolution from low cell count, low-coverage, or low-quality ATAC-seq data. Models trained by AtacWorks can detect peaks from cell types not seen in the training data, and are generalizable across diverse sample preparations and experimental platforms. We demonstrate that AtacWorks enhances the sensitivity of single-cell experiments by producing results on par with those of conventional methods using ~10 times as many cells, and further show that this framework can be adapted to enable cross-modality inference of protein-DNA interactions. Finally, we establish that AtacWorks can enable new biological discoveries by identifying active regulatory regions associated with lineage priming in rare subpopulations of hematopoietic stem cells.

ATAC-seq measures chromatin accessibility as a proxy for the activity of DNA regulatory regions across the genome. Here the authors present AtacWorks, a deep learning tool to denoise and identify accessible chromatin regions from low cell count, low-coverage, or low-quality ATAC-seq data.

In situ genome sequencing resolves DNA sequence and structure in intact biological samples

Understanding genome organization requires integration of DNA sequence and three-dimensional spatial context; however, existing genome-wide methods lack either base pair sequence resolution or direct spatial localization. Here, we describe in situ genome sequencing (IGS), a method for simultaneously sequencing and imaging genomes within intact biological samples. We applied IGS to human fibroblasts and early mouse embryos, spatially localizing thousands of genomic loci in individual nuclei. Using these data, we characterized parent-specific changes in genome structure across embryonic stages, revealed single-cell chromatin domains in zygotes, and uncovered epigenetic memory of global chromosome positioning within individual embryos. These results demonstrate how IGS can directly connect sequence and structure across length scales from single base pairs to whole organisms.

Chromatin Potential Identified by Shared Single-Cell Profiling of RNA and Chromatin

Cell differentiation and function are regulated across multiple layers of gene regulation, including modulation of gene expression by changes in chromatin accessibility. However, differentiation is an asynchronous process precluding a temporal understanding of regulatory events leading to cell fate commitment. Here we developed simultaneous high-throughput ATAC and RNA expression with sequencing (SHARE-seq), a highly scalable approach for measurement of chromatin accessibility and gene expression in the same single cell, applicable to different tissues. Using 34,774 joint profiles from mouse skin, we develop a computational strategy to identify cis-regulatory interactions and define domains of regulatory chromatin (DORCs) that significantly overlap with super-enhancers. During lineage commitment, chromatin accessibility at DORCs precedes gene expression, suggesting that changes in chromatin accessibility may prime cells for lineage commitment. We computationally infer chromatin potential as a quantitative measure of chromatin lineage-priming and use it to predict cell fate outcomes. SHARE-seq is an extensible platform to study regulatory circuitry across diverse cells in tissues.

Advances in Chromatin and Chromosome Research: Perspectives from Multiple Fields

Nucleosomes package genomic DNA into chromatin. By regulating DNA access for transcription, replication, DNA repair, and epigenetic modification, chromatin forms the nexus of most nuclear processes. In addition, dynamic organization of chromatin underlies both regulation of gene expression and evolution of chromosomes into individualized sister objects, which can segregate cleanly to different daughter cells at anaphase. This collaborative review shines a spotlight on technologies that will be crucial to interrogate key questions in chromatin and chromosome biology including state-of-the-art microscopy techniques, tools to physically manipulate chromatin, single-cell methods to measure chromatin accessibility, computational imaging with neural networks and analytical tools to interpret chromatin structure and dynamics. In addition, this review provides perspectives on how these tools can be applied to specific research fields such as genome stability and developmental biology and to test concepts such as phase separation of chromatin.

Abstract 5879: Spatially-segmented single-cell transcriptomics by diffusional accessibility to a small-molecule dye

High-throughput single-cell RNA-seq (scRNA-seq) is used to describe complex tissues by characterizing transcriptional states of individual cells. Defining a cell's position, both in regard to tissue margins and its social context, is essential for understanding the intrinsic and extrinsic variables that effect the transcriptional identity of individual cells. Conventional high-throughput scRNA-seq assays, however, decouple cells from their original locations within tissues. In situ hybridization readouts of gene expression and in situ sequencing preserve spatial information in tissues, but currently have a lower total read threshold than NextGen sequencing – imposing a restriction on either cell throughput or transcriptional breadth. Combining the above methods or using regional barcodes to define 2D positions have spatially-reconstructed tissue regions but are seldom employed in unspecialized laboratories. We describe SEgmentation by Exogenous Perfusion (SEEP), a rapid, integrated, method for providing 3D spatial-segmentation to scRNA-seq data. Tissues are divided into layers based on accessibility of a fluorescent dye allowing sorted cells to be characterized by transcriptomic and regional identity. We use SEEP to explore how the transcriptional states of cells in high-grade serous ovarian cancer vary with respect to intratumoral position. We describe an integrated method for correlating 3D radial-spatial cell positions with scRNA-seq data. By employing a basic stain-and-sort method using an off-the-shelf live-dead stain, we are able to correlate transcriptional profiles with radial-spatial positions in both radially symmetric tissues (e.g., spheroids, organoids, spherical tumor masses) and linear tissue samples (e.g., punch biopsies). SEEP uses an imaging-based calibration step to inform the parameters of a FACS and sequencing based measurement step. Across spheroid, organoid, and PDX models of high-grade serous ovarian cancer, we find common positional transcriptomic profiles, ranging from single-gene proclivities to multi-gene-set enrichments across tissue layers.

Citation Format: David B. Morse, Michele Ceribelli, Joachim De Jonghe, Aleksandra Michalowski, Christoph Muus, Maria Vias, Samantha Boyle, David A. Weitz, James Brenton, Jason D. Buenrostro, Craig Thomas, Tuomas P. Knowles. Spatially-segmented single-cell transcriptomics by diffusional accessibility to a small-molecule dye [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5879.

Epigenomic State Transitions Characterize Tumor Progression in Mouse Lung Adenocarcinoma

Regulatory networks that maintain functional, differentiated cell states are often dysregulated in tumor development. Here, we use single-cell epigenomics to profile chromatin state transitions in a mouse model of lung adenocarcinoma (LUAD). We identify an epigenomic continuum representing loss of cellular identity and progression toward a metastatic state. We define co-accessible regulatory programs and infer key activating and repressive chromatin regulators of these cell states. Among these co-accessibility programs, we identify a pre-metastatic transition, characterized by activation of RUNX transcription factors, which mediates extracellular matrix remodeling to promote metastasis and is predictive of survival across human LUAD patients. Together, these results demonstrate the power of single-cell epigenomics to identify regulatory programs to uncover mechanisms and key biomarkers of tumor progression.

Sustained Activation of AMPK Enhances Differentiation of Human iPSC-Derived Cardiomyocytes via Sirtuin Activation

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).

Abstract PR08: Leveraging single-cell epigenomics to uncover regulatory programs in lung adenocarcinoma

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.

Inference and effects of barcode multiplets in droplet-based single-cell assays

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.

Hyperactivation of sympathetic nerves drives depletion of melanocyte stem cells

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, chemogenetics^3,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.

Assessment of computational methods for the analysis of single-cell ATAC-seq data

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).

Single-cell trajectories reconstruction, exploration and mapping of omics data with STREAM

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.

Interrogation of human hematopoiesis at single-cell and single-variant resolution

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.

The cis-Regulatory Atlas of the Mouse Immune System

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.

Neutralizing Gatad2a-Chd4-Mbd3/NuRD Complex Facilitates Deterministic Induction of Naive Pluripotency

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.

Integrated Single-Cell Analysis Maps the Continuous Regulatory Landscape of Human Hematopoietic Differentiation

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.

Transcript-indexed ATAC-seq for precision immune profiling

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.

Rapid chromatin repression by Aire provides precise control of immune tolerance

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.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

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.

chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data

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.

PHF6 regulates phenotypic plasticity through chromatin organization within lineage-specific genes

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, Phf6^KO 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.

TOP2 synergizes with BAF chromatin remodeling for both resolution and formation of facultative heterochromatin

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.

Epigenomics of human CD8 T cell differentiation and aging

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.

Single-cell epigenomic variability reveals functional cancer heterogeneity

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.

ATAC-see reveals the accessible genome by transposase-mediated imaging and sequencing

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.

Lineage-specific and single-cell chromatin accessibility charts human hematopoiesis and leukemia evolution

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.

Leukemia-Associated Cohesin Mutants Dominantly Enforce Stem Cell Programs and Impair Human Hematopoietic Progenitor Differentiation

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.

Transient acquisition of pluripotency during somatic cell transdifferentiation with iPSC reprogramming factors

Somatic cells can be transdifferentiated to other cell types without passing through a pluripotent state by ectopic expression of appropriate transcription factors^1,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 differentiation^3–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.

ATAC-seq: A Method for Assaying Chromatin Accessibility Genome-Wide

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.

Exome sequencing identifies a DNAJB6 mutation in a family with dominantly-inherited limb-girdle muscular dystrophy

Limb-girdle muscular dystrophy primarily affects the muscles of the hips and shoulders (the "limb-girdle" muscles), although it is a heterogeneous disorder that can present with varying symptoms. There is currently no cure. We sought to identify the genetic basis of limb-girdle muscular dystrophy type 1 in an American family of Northern European descent using exome sequencing. Exome sequencing was performed on DNA samples from two affected siblings and one unaffected sibling and resulted in the identification of eleven candidate mutations that co-segregated with the disease. Notably, this list included a previously reported mutation in DNAJB6, p.Phe89Ile, which was recently identified as a cause of limb-girdle muscular dystrophy type 1D. Additional family members were Sanger sequenced and the mutation in DNAJB6 was only found in affected individuals. Subsequent haplotype analysis indicated that this DNAJB6 p.Phe89Ile mutation likely arose independently of the previously reported mutation. Since other published mutations are located close by in the G/F domain of DNAJB6, this suggests that the area may represent a mutational hotspot. Exome sequencing provided an unbiased and effective method for identifying the genetic etiology of limb-girdle muscular dystrophy type 1 in a previously genetically uncharacterized family. This work further confirms the causative role of DNAJB6 mutations in limb-girdle muscular dystrophy type 1D.