Skip to content

noobCoding/CellEnrich

 
 

Repository files navigation

CellEnrich

Pathway enrichment analysis and visualization for Single Cell Data

Installation

NOTE: on a fresh installation, users may need to install some required interpreter compilers for the system to install other R packages:

  • C++ compiler
  • gfortran compiler (FYI: tips for MAC users or other OS)
  • Seurat >= 5.0.0 is REQUIRED

Install Dependent Packages:

  • Some packages require to be installed from sources that need compilation and a proper version of RTools.
  • It’s recommended to install devtools, and BiocManager packages first, to install the following packages in Github / Bioconductor (not in CRAN).
# install required packages
install.packages('Seurat') # RStudio may need a RESTART for Seurat v5.0.1 to be activated
install.packages('remotes')
install.packages('waiter')
install.packages('farver')
remotes::install_github('vqv/ggbiplot')

# install CellEnrich
remotes::install_github('noobCoding/CellEnrich')

Example with PBMC_3K data

# Download data, if not downloaded
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/pbmcData.RData','pbmcData.RData', mode = 'wb')
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/pbmcClustInfo.RData','pbmcClustInfo.RData', mode = 'wb')
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/Human_Reactome.RData', 'Human_Reactome.RData', mode = 'wb')

# Load library and data
library(CellEnrich)
library(Seurat)

load("pbmcData.RData")
load("pbmcClustInfo.RData")

CountData <- pbmcData
GroupInfo <- pbmcClustInfo

# CellEnrich uses normalized count data as input
CountData <- NormalizeData(CountData)

# This will run CellEnrich
CellEnrich(CountData, GroupInfo)

Example with Alzheimer's data

# Download Alzheimer's data to the working directory
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/Alzheimer_Counts_sampled.RDS','Alzheimer_Counts_sampled.RDS', mode = 'wb')
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/Alzheimer_CellType_sampled.RDS','Alzheimer_CellType_sampled.RDS', mode = 'wb')
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/Human_Reactome.RData', 'Human_Reactome.RData', mode = 'wb')

# Load library and data
library(CellEnrich)
library(Seurat)

GroupInfo <- readRDS("Alzheimer_CellType_sampled.RDS")
CountData <- readRDS("Alzheimer_Counts_sampled.RDS")

# CellEnrich uses normalized count data as input
CountData <- NormalizeData(CountData)

# Run cellenrich
CellEnrich(CountData, GroupInfo)

Example with primary mouse dendritic cells (DCs) stimulated with lipopolysaccharide (LPS)

# download minimal data to the working directory
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/LPS_exp.rds','LPS_exp.rds', mode = 'wb')
download.file('https://github.com/noobcoding/CellEnrich/raw/master/data/Mouse_WikiPathways.RData', 'Mouse_WikiPathways.RData', mode = 'wb')

# Load library and data
library(Seurat)
library(CellEnrich)

LPS_exp <- readRDS("LPS_exp.rds")
CountData <- LPS_exp$counts
GroupInfo <- LPS_exp$meta

# CellEnrich uses normalized count data as input
CountData <- NormalizeData(CountData)

# Run cellenrich
CellEnrich(CountData, GroupInfo)

Dependency

Authors

License

This project is MIT licensed