Dataset Information

Comprehensive analysis of single cell ATAC-seq data with SnapATAC.

ABSTRACT: Identification of the cis-regulatory elements controlling cell-type specific gene expression patterns is essential for understanding the origin of cellular diversity. Conventional assays to map regulatory elements via open chromatin analysis of primary tissues is hindered by sample heterogeneity. Single cell analysis of accessible chromatin (scATAC-seq) can overcome this limitation. However, the high-level noise of each single cell profile and the large volume of data pose unique computational challenges. Here, we introduce SnapATAC, a software package for analyzing scATAC-seq datasets. SnapATAC dissects cellular heterogeneity in an unbiased manner and map the trajectories of cellular states. Using the Nyström method, SnapATAC can process data from up to a million cells. Furthermore, SnapATAC incorporates existing tools into a comprehensive package for analyzing single cell ATAC-seq dataset. As demonstration of its utility, SnapATAC is applied to 55,592 single-nucleus ATAC-seq profiles from the mouse secondary motor cortex. The analysis reveals ~370,000 candidate regulatory elements in 31 distinct cell populations in this brain region and inferred candidate cell-type specific transcriptional regulators.

SUBMITTER: Fang R

PROVIDER: S-EPMC7910485 | biostudies-literature | 2021 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Comprehensive analysis of single cell ATAC-seq data with SnapATAC.

Fang Rongxin R Preissl Sebastian S Li Yang Y Hou Xiaomeng X Lucero Jacinta J Wang Xinxin X Motamedi Amir A Shiau Andrew K AK Zhou Xinzhu X Xie Fangming F Mukamel Eran A EA Zhang Kai K Zhang Yanxiao Y Behrens M Margarita MM Ecker Joseph R JR Ren Bing B

Nature communications 20210226 1

Identification of the cis-regulatory elements controlling cell-type specific gene expression patterns is essential for understanding the origin of cellular diversity. Conventional assays to map regulatory elements via open chromatin analysis of primary tissues is hindered by sample heterogeneity. Single cell analysis of accessible chromatin (scATAC-seq) can overcome this limitation. However, the high-level noise of each single cell profile and the large volume of data pose unique computational c ...[more]

PMID: 33637727

Similar Datasets

Project description:Single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data provided new insights into the understanding of epigenetic heterogeneity and transcriptional regulation. With the increasing abundance of dataset resources, there is an urgent need to extract more useful information through high-quality data analysis methods specifically designed for scATAC-seq. However, analyzing scATAC-seq data poses challenges due to its near binarization, high sparsity and ultra-high dimensionality properties. Here, we proposed a novel network diffusion-based computational method to comprehensively analyze scATAC-seq data, named Single-Cell ATAC-seq Analysis via Network Refinement with Peaks Location Information (SCARP). SCARP formulates the Network Refinement diffusion method under the graph theory framework to aggregate information from different network orders, effectively compensating for missing signals in the scATAC-seq data. By incorporating distance information between adjacent peaks on the genome, SCARP also contributes to depicting the co-accessibility of peaks. These two innovations empower SCARP to obtain lower-dimensional representations for both cells and peaks more effectively. We have demonstrated through sufficient experiments that SCARP facilitated superior analyses of scATAC-seq data. Specifically, SCARP exhibited outstanding cell clustering performance, enabling better elucidation of cell heterogeneity and the discovery of new biologically significant cell subpopulations. Additionally, SCARP was also instrumental in portraying co-accessibility relationships of accessible regions and providing new insight into transcriptional regulation. Consequently, SCARP identified genes that were involved in key Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to diseases and predicted reliable cis-regulatory interactions. To sum up, our studies suggested that SCARP is a promising tool to comprehensively analyze the scATAC-seq data.

Project description:PurposeCorneal epithelial homeostasis is maintained by coordinated gene expression across distinct cell populations, but the gene regulatory programs underlying this cellular diversity remain to be characterized. Here we applied single-cell multi-omics analysis to delineate the gene regulatory profile of mouse corneal epithelial cells under normal homeostasis.MethodsSingle cells isolated from the cornea epithelium (with marginal conjunctiva) of adult mice were subjected to scRNA-seq and scATAC-seq using the 10×Genomics platform. Cell types were clustered by the graph-based visualization method uniform manifold approximation and projection and unbiased computational informatics analysis. The scRNA-seq and scATAC-seq datasets were integrated following the integration pipeline described in ArchR and Seurat.ResultsWe characterized diverse corneal epithelial cell types based on gene expression signatures and chromatin accessibility. We found that cell type-specific accessibility regions were mainly located at distal regions, suggesting essential roles of distal regulatory elements in determining corneal epithelial cell diversity. Trajectory analyses revealed a continuum of cell state transition and higher coordination between transcription factor (TF) motif accessibility and gene expression during corneal epithelial cell differentiation. By integrating transcriptomic and chromatin accessibility analysis, we identified cell type-specific and shared gene regulation programs. We also uncovered critical TFs driving corneal epithelial cell differentiation, such as nuclear factor I (NFI) family members, Rarg, Elf3. We found that nuclear factor-κB (NF-κB) family members were positive TFs in limbal cells and some superficial cells, but they were involved in regulating distinct biological processes.ConclusionsOur study presents a comprehensive gene regulatory landscape of mouse cornea epithelial cells, and provides valuable foundations for future investigation of corneal epithelial homeostasis in the context of cornea pathologies and regenerative medicine.

Dataset Information

Comprehensive analysis of single cell ATAC-seq data with SnapATAC.

Publications

Comprehensive analysis of single cell ATAC-seq data with SnapATAC.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets