scholarly journals A rapid and robust method for single cell chromatin accessibility profiling

2018 ◽  
Author(s):  
Xi Chen ◽  
Ricardo J Miragaia ◽  
Kedar Nath Natarajan ◽  
Sarah A Teichmann

AbstractThe assay for transposase-accessible chromatin using sequencing (ATAC-seq) is widely used to identify regulatory regions throughout the genome. However, very few studies have been performed at the single cell level (scATAC-seq) due to technical challenges. Here we developed a simple and robust plate-based scATAC-seq method, combining upfront bulk Tn5 tagging with single-nuclei sorting. We demonstrated that our method worked robustly across various systems, including fresh and cryopreserved cells from primary tissues. By profiling over 3,000 splenocytes, we identify distinct immune cell types and reveal cell type-specific regulatory regions and related transcription factors.

2021 ◽  
Vol 12 ◽  
Author(s):  
Zhe Cui ◽  
Ya Cui ◽  
Yan Gao ◽  
Tao Jiang ◽  
Tianyi Zang ◽  
...  

Single-cell Assay Transposase Accessible Chromatin sequencing (scATAC-seq) has been widely used in profiling genome-wide chromatin accessibility in thousands of individual cells. However, compared with single-cell RNA-seq, the peaks of scATAC-seq are much sparser due to the lower copy numbers (diploid in humans) and the inherent missing signals, which makes it more challenging to classify cell type based on specific expressed gene or other canonical markers. Here, we present svmATAC, a support vector machine (SVM)-based method for accurately identifying cell types in scATAC-seq datasets by enhancing peak signal strength and imputing signals through patterns of co-accessibility. We applied svmATAC to several scATAC-seq data from human immune cells, human hematopoietic system cells, and peripheral blood mononuclear cells. The benchmark results showed that svmATAC is free of literature-based markers and robust across datasets in different libraries and platforms. The source code of svmATAC is available at https://github.com/mrcuizhe/svmATAC under the MIT license.


2021 ◽  
Author(s):  
Risa Karakida Kawaguchi ◽  
Ziqi Tang ◽  
Stephan Fischer ◽  
Rohit Tripathy ◽  
Peter K. Koo ◽  
...  

Background: Single-cell Assay for Transposase Accessible Chromatin using sequencing (scATAC-seq) measures genome-wide chromatin accessibility for the discovery of cell-type specific regulatory networks. ScATAC-seq combined with single-cell RNA sequencing (scRNA-seq) offers important avenues for ongoing research, such as novel cell-type specific activation of enhancer and transcription factor binding sites as well as chromatin changes specific to cell states. On the other hand, scATAC-seq data is known to be challenging to interpret due to its high number of zeros as well as the heterogeneity derived from different protocols. Because of the stochastic lack of marker gene activities, cell type identification by scATAC-seq remains difficult even at a cluster level. Results: In this study, we exploit reference knowledge obtained from external scATAC-seq or scRNA-seq datasets to define existing cell types and uncover the genomic regions which drive cell-type specific gene regulation. To investigate the robustness of existing cell-typing methods, we collected 7 scATAC-seq datasets targeting mouse brain for a meta-analytic comparison of neuronal cell-type annotation, including a reference atlas generated by the BRAIN Initiative Cell Census Network (BICCN). By comparing the area under the receiver operating characteristics curves (AUROCs) for the three major cell types (inhibitory, excitatory, and non-neuronal cells), cell-typing performance by single markers is found to be highly variable even for known marker genes due to study-specific biases. However, the signal aggregation of a large and redundant marker gene set, optimized via multiple scRNA-seq data, achieves the highest cell-typing performances among 5 existing marker gene sets, from the individual cell to cluster level. That gene set also shows a high consistency with the cluster-specific genes from inhibitory subtypes in two well-annotated datasets, suggesting applicability to rare cell types. Next, we demonstrate a comprehensive assessment of scATAC-seq cell typing using exhaustive combinations of the marker gene sets with supervised learning methods including machine learning classifiers and joint clustering methods. Our results show that the combinations using robust marker gene sets systematically ranked at the top, not only with model based prediction using a large reference data but also with a simple summation of expression strengths across markers. To demonstrate the utility of this robust cell typing approach, we trained a deep neural network to predict chromatin accessibility in each subtype using only DNA sequence. Through model interpretation methods, we identify key motifs enriched about robust gene sets for each neuronal subtype. Conclusions: Through the meta-analytic evaluation of scATAC-seq cell-typing methods, we develop a novel method set to exploit the BICCN reference atlas. Our study strongly supports the value of robust marker gene selection as a feature selection tool and cross-dataset comparison between scATAC-seq datasets to improve alignment of scATAC-seq to known biology. With this novel, high quality epigenetic data, genomic analysis of regulatory regions can reveal sequence motifs that drive cell type-specific regulatory programs.


2020 ◽  
Author(s):  
Alexandre P. Marand ◽  
Zongliang Chen ◽  
Andrea Gallavotti ◽  
Robert J. Schmitz

ABSTRACTCis-regulatory elements (CREs) encode the genomic blueprints for coordinating spatiotemporal gene expression programs underlying highly specialized cell functions. To identify CREs underlying cell-type specification and developmental transitions, we implemented single-cell sequencing of Assay for Transposase Accessible Chromatin in an atlas of Zea mays organs. We describe 92 distinct states of chromatin accessibility across more than 165,913 putative CREs, 56,575 cells, and 52 known cell-types in maize using a novel implementation of regularized quasibinomial logistic regression. Cell states were largely determined by combinatorial accessibility of transcription factors (TFs) and their binding sites. A neural network revealed that cell identity could be accurately predicted (>0.94) solely based on TF binding site accessibility. Co-accessible chromatin recapitulated higher-order chromatin interactions, with distinct sets of TFs coordinating cell type-specific regulatory dynamics. Pseudotime reconstruction and alignment with Arabidopsis thaliana trajectories identified conserved TFs, associated motifs, and cis-regulatory regions specifying sequential developmental progressions. Cell-type specific accessible chromatin regions were enriched with phenotype-associated genetic variants and signatures of selection, revealing the major cell-types and putative CREs targeted by modern maize breeding. Collectively, our analysis affords a comprehensive framework for understanding cellular heterogeneity, evolution, and cis-regulatory grammar of cell-type specification in a major crop species.


2022 ◽  
Author(s):  
Matthew T Buckley ◽  
Eric Sun ◽  
Benson M. George ◽  
Ling Liu ◽  
Nicholas Schaum ◽  
...  

Aging manifests as progressive dysfunction culminating in death. The diversity of cell types is a challenge to the precise quantification of aging and its reversal. Here we develop a suite of 'aging clocks' based on single cell transcriptomic data to characterize cell type-specific aging and rejuvenation strategies. The subventricular zone (SVZ) neurogenic region contains many cell types and provides an excellent system to study cell-level tissue aging and regeneration. We generated 21,458 single-cell transcriptomes from the neurogenic regions of 28 mice, tiling ages from young to old. With these data, we trained a suite of single cell-based regression models (aging clocks) to predict both chronological age (passage of time) and biological age (fitness, in this case the proliferative capacity of the neurogenic region). Both types of clocks perform well on independent cohorts of mice. Genes underlying the single cell-based aging clocks are mostly cell-type specific, but also include a few shared genes in the interferon and lipid metabolism pathways. We used these single cell-based aging clocks to measure transcriptomic rejuvenation, by generating single cell RNA-seq datasets of SVZ neurogenic regions for two interventions - heterochronic parabiosis (young blood) and exercise. Interestingly, the use of aging clocks reveals that both heterochronic parabiosis and exercise reverse transcriptomic aging in the niche, but in different ways across cell types and genes. This study represents the first development of high-resolution aging clocks from single cell transcriptomic data and demonstrates their application to quantify transcriptomic rejuvenation.


2020 ◽  
Vol 25 (7) ◽  
pp. 755-769 ◽  
Author(s):  
Clinton Willis ◽  
Johanna Nyffeler ◽  
Joshua Harrill

Cell Painting is a high-throughput phenotypic profiling assay that uses fluorescent cytochemistry to visualize a variety of organelles and high-content imaging to derive a large number of morphological features at the single-cell level. Most Cell Painting studies have used the U-2 OS cell line for chemical or functional genomics screening. The Cell Painting assay can be used with many other human-derived cell types, given that the assay is based on the use of fluoroprobes that label organelles that are present in most (if not all) human cells. Questions remain, however, regarding the optimization steps required and overall ease of deployment of the Cell Painting assay to novel cell types. Here, we used the Cell Painting assay to characterize the phenotypic effects of 14 phenotypic reference chemicals in concentration–response screening mode across six biologically diverse human-derived cell lines (U-2 OS, MCF7, HepG2, A549, HTB-9 and ARPE-19). All cell lines were labeled using the same cytochemistry protocol, and the same set of phenotypic features was calculated. We found it necessary to optimize image acquisition settings and cell segmentation parameters for each cell type, but did not adjust the cytochemistry protocol. For some reference chemicals, similar subsets of phenotypic features corresponding to a particular organelle were associated with the highest-effect magnitudes in each affected cell type. Overall, for certain chemicals, the Cell Painting assay yielded qualitatively similar biological activity profiles among a group of diverse, morphologically distinct human-derived cell lines without the requirement for cell type–specific optimization of cytochemistry protocols.


2020 ◽  
Vol 52 (9) ◽  
pp. 1428-1442 ◽  
Author(s):  
Jeongwoo Lee ◽  
Do Young Hyeon ◽  
Daehee Hwang

Abstract Advances in single-cell isolation and barcoding technologies offer unprecedented opportunities to profile DNA, mRNA, and proteins at a single-cell resolution. Recently, bulk multiomics analyses, such as multidimensional genomic and proteogenomic analyses, have proven beneficial for obtaining a comprehensive understanding of cellular events. This benefit has facilitated the development of single-cell multiomics analysis, which enables cell type-specific gene regulation to be examined. The cardinal features of single-cell multiomics analysis include (1) technologies for single-cell isolation, barcoding, and sequencing to measure multiple types of molecules from individual cells and (2) the integrative analysis of molecules to characterize cell types and their functions regarding pathophysiological processes based on molecular signatures. Here, we summarize the technologies for single-cell multiomics analyses (mRNA-genome, mRNA-DNA methylation, mRNA-chromatin accessibility, and mRNA-protein) as well as the methods for the integrative analysis of single-cell multiomics data.


Author(s):  
Ryan S. Ziffra ◽  
Chang N. Kim ◽  
Amy Wilfert ◽  
Tychele N. Turner ◽  
Maximilian Haeussler ◽  
...  

AbstractDynamic changes in chromatin accessibility coincide with important aspects of neuronal differentiation, such as fate specification and arealization and confer cell type-specific associations to neurodevelopmental disorders. However, studies of the epigenomic landscape of the developing human brain have yet to be performed at single-cell resolution. Here, we profiled chromatin accessibility of >75,000 cells from eight distinct areas of developing human forebrain using single cell ATAC-seq (scATACseq). We identified thousands of loci that undergo extensive cell type-specific changes in accessibility during corticogenesis. Chromatin state profiling also reveals novel distinctions between neural progenitor cells from different cortical areas not seen in transcriptomic profiles and suggests a role for retinoic acid signaling in cortical arealization. Comparison of the cell type-specific chromatin landscape of cerebral organoids to primary developing cortex found that organoids establish broad cell type-specific enhancer accessibility patterns similar to the developing cortex, but lack many putative regulatory elements identified in homologous primary cell types. Together, our results reveal the important contribution of chromatin state to the emerging patterns of cell type diversity and cell fate specification and provide a blueprint for evaluating the fidelity and robustness of cerebral organoids as a model for cortical development.


2020 ◽  
Vol 8 (Suppl 2) ◽  
pp. A4.2-A5
Author(s):  
CM Schürch ◽  
DJ Phillips ◽  
M Matusiak ◽  
B Rivero Gutierrez ◽  
SS Bhate ◽  
...  

BackgroundImmunotherapies have induced long-lasting remissions in countless advanced-stage cancer patients, but many more patients have not benefitted. Therefore, novel predictive markers are needed to stratify patients before treatment and select those who will most likely benefit from immunotherapy, while avoiding potentially devastating adverse effects and high treatment costs for those who will not. We reasoned that thoroughly characterizing the architecture of the tumor microenvironment (TME) at the single-cell level by highly multiplexed tissue imaging should reveal novel spatial biomarkers of immunotherapy response.Materials and MethodsWe used CODEX (CO-Detection by indEXing) highly multiplexed fluorescence microscopy to investigate the TME of cutaneous T cell lymphoma (CTCL) in samples from patients treated with pembrolizumab. 55 protein markers were visualized simultaneously using a tissue microarray of matched pre- and post-treatment skin biopsies from 7 pembrolizumab responders and 7 non-responders. After computational image processing and extraction of single-cell information, cell types were identified by unsupervised clustering followed by supervised curation, and cell-cell distances and ‘cellular neighborhoods’ were computed. We also performed RNA sequencing on laser-capture microdissected tissue microarray cores to extract cell-type specific gene expression profiles by CIBERSORTx analysis.ResultsCODEX enabled the identification and characterization of malignant CD4+ tumor cells and reactive immune cells in the CTCL TME at the single-cell level, resulting in 21 different cell type clusters with spatial information. Cluster frequencies were not significantly different between responders and non-responders pre- and post-treatment. However, advanced computational analysis of the tumor architecture revealed cellular neighborhoods (CNs) that dynamically changed during pembrolizumab therapy and were correlated with response. Effector-type CNs enriched in tumor-infiltrating CD4+ T cells and dendritic cells were significantly increased after treatment in responders. In contrast, a regulatory T cell-enriched CN was significantly increased in non-responders before and after therapy. Furthermore, a spatial signature of cell-cell distances between tumor cells and effector/regulatory immune cells predicted therapy outcome. In addition, CIBERSORTx analysis revealed that tumor cells in responders, but not in non-responders, increased their expression of immune-activating genes.ConclusionsHigh-dimensional spatial analysis of CTCL tumors revealed a pre-existing immunosuppressive state in pembrolizumab non-responders. Thorough analysis of the TME therefore enables the discovery of novel spatial biomarkers in a concept that accounts for both cell type information and higher-order tumor architecture. Combining highly multiplexed microscopy with CIBERSORTx allows for the discovery of novel, predictive spatial biomarkers of immunotherapy response and will pave the way for future studies that functionally address these cell types and their interactions.Disclosure InformationC.M. Schürch: F. Consultant/Advisory Board; Modest; Enable Medicine, LLC. D.J. Phillips: None. M. Matusiak: None. B. Rivero Gutierrez: None. S.S. Bhate: None. G.L. Barlow: None. M.S. Khodadoust: B. Research Grant (principal investigator, collaborator or consultant and pending grants as well as grants already received); Significant; Corvus Pharmaceuticals. R. West: None. Y.H. Kim: B. Research Grant (principal investigator, collaborator or consultant and pending grants as well as grants already received); Significant; Merck, Horizon, Soligenix, miRagen, Forty Seven Inc., Neumedicine, Trillium, Galderma, Elorac. D. Speakers Bureau/Honoraria (speakers bureau, symposia, and expert witness); Significant; Innate Pharma, Eisai, Kyowa Hakko Kirin, Takeda, Seattle Genetics, Medivir, Portola Pharmaceuticals, Corvus Pharmaceuticals. G.P. Nolan: E. Ownership Interest (stock, stock options, patent or other intellectual property); Significant; Akoya Biosciences. F. Consultant/Advisory Board; Significant; Akoya Biosciences.


2020 ◽  
Author(s):  
Paola Benaglio ◽  
Jacklyn Newsome ◽  
Jee Yun Han ◽  
Joshua Chiou ◽  
Anthony Aylward ◽  
...  

AbstractGene regulation is highly cell type-specific and understanding the function of non-coding genetic variants associated with complex traits requires molecular phenotyping at cell type resolution. In this study we performed single nucleus ATAC-seq (snATAC-seq) and genotyping in peripheral blood mononuclear cells from 10 individuals. Clustering chromatin accessibility profiles of 66,843 total nuclei identified 14 immune cell types and sub-types. We mapped chromatin accessibility QTLs (caQTLs) in each immune cell type and sub-type which identified 6,248 total caQTLs, including those obscured from assays of bulk tissue such as with divergent effects on different cell types. For 3,379 caQTLs we further annotated putative target genes of variant activity using single cell co-accessibility, and caQTL variants were significantly correlated with the accessibility level of linked gene promoters. We fine-mapped loci associated with 16 complex immune traits and identified immune cell caQTLs at 517 candidate causal variants, including those with cell type-specific effects. At the 6q15 locus associated with type 1 diabetes, in line with previous reports, variant rs72928038 was a naïve CD4+ T cell caQTL linked to BACH2 and we validated the allelic effects of this variant on regulatory activity in Jurkat T cells. These results highlight the utility of snATAC-seq for mapping genetic effects on accessible chromatin in specific cell types and provide a resource for annotating complex immune trait loci.


2021 ◽  
Author(s):  
Tianxing Ma ◽  
Haochen Li ◽  
Xuegong Zhang

eQTL studies are essential for understanding genomic regulation. Effects of genetic variations on gene regulation are cell-type-specific and cellular-context-related, so studying eQTLs at a single-cell level is crucial. The ideal solution is to use both mutation and expression data from the same cells. However, current technology of such paired data in single cells is still immature. We present a new method, eQTLsingle, to discover eQTLs only with single cell RNA-seq (scRNA-seq) data, without genomic data. It detects mutations from scRNA-seq data and models gene expression of different genotypes with the zero-inflated negative binomial (ZINB) model to find associations between genotypes and phenotypes at single-cell level. On a glioblastoma and gliomasphere scRNA-seq dataset, eQTLsingle discovered hundreds of cell-type-specific tumor-related eQTLs, most of which cannot be found in bulk eQTL studies. Detailed analyses on examples of the discovered eQTLs revealed important underlying regulatory mechanisms. eQTLsingle is a unique powerful tool for utilizing the huge scRNA-seq resources for single-cell eQTL studies, and it is available for free academic use at https://github.com/horsedayday/eQTLsingle.


Sign in / Sign up

Export Citation Format

Share Document