Skip to content

danielcgingerich/LOAD_snRNA_snATAC-seq

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

1 Commit
 
 
 
 
 
 
 
 

Repository files navigation

Integrative analysis of parallel single nucleus ATAC and RNA seq data

This repository contains the analysis pipeline used in Gamache, Gingerich & Shwab et al. “Integrative single-nucleus multi-omics analysis prioritizes candidate cis and trans regulatory networks and their target genes in Alzheimer’s disease brains”. The pipeline was created for examining the genomic and epigenomic landscape of Alzheimer’s disease in the temporal cortex, but can be used for a myriad of tissue regions and pathologies. In addition to cross modality integration at the level of cell type, this pipeline uses Seurat's implementation of diagonalized CCA to match snATAC subclusters to snRNA subclusters. The method of subcluster matching was validated using 10X genomics PBMC multiome data and found to be ~70% accurate (see /atac_pipeline/7.cluster_link_validation/). Be sure to keep this in mind when interpreting downstream results.

config.R file

This file is sourced at the beginning of each script. If working on multiple computers or HPCs, the file detects what machine R is running on based on the path to the working directory (grepl(‘/hpc_path_identifier’, getwd()). Then, directory paths and library paths are set accordingly. The config file also silently loads commonly used packages to avoid repetitive use of library() at the beginning of each script.

Directory layout

R scripts in each subdirectory of the pipeline save all output files to a respective user-specified directory. These directories can be specified in the config.R file.
R scripts in this repository are located as follows:

./config.R (this sets all user-specified paths from which to read/write files)

/rna_pipeline/
 1.preprocessing/
 2.integration/
 3.cell_type_annotation/
 4.deg_analysis/

/atac_pipeline/
 1.preprocessing/
 2.cell_type_annotation/
 3.integration/
 4.dap_analysis/
 5.cicero/
 6.motif_enrichment/
 7.cluster_link_validation/

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • R 100.0%