VoxHunt is a package for assessing brain organoid patterning, developmental state, and cell composition through systematic comparisons of single cell transcriptomes to three-dimensional in situ hybridization data from the Allen Brain Atlas and a number of other useful reference datasets.
VoxHunt currently implements the following features:
Presto, one of VoxHunt’s dependencies is not on CRAN and has to be installed from GitHub:
# install.packages('devtools')
devtools::install_github('immunogenomics/presto')
If you want to perform deconvolution, you also have to install EPIC manually:
devtools::install_github('GfellerLab/EPIC')
Once all external dependencies are installed, you can install VoxHunt with
devtools::install_github('quadbiolab/voxhunt')
In addition to the R package itself, you’ll also need to download the the ABA expression data from here.
If you have a seurat_object
with single cell transcriptomic data of your organoid ready, you can start right away with projecting them to the brain:
# Load VoxHunt
library(voxhunt)
# Point VoxHunt to ABA expression data
load_aba_data('~/path/to/data')
# Find 300 most variable genes from the E13.5 mouse brain
genes_use <- variable_genes('E13', 300)$gene
# Calculate the similarity map of a seurat object to the E13.5 mouse brain
vox_map <- voxel_map(seurat_object, genes_use=genes_use)
# Plot the result
plot_map(vox_map)
After loading VoxHunt, we first point it to the directory with the unpacked ABA expression data. Then, we select the 300 most variable features from the E13.5 mouse brain ISH data and use them to calculate similarity maps for your organoid cells. Finally, we plot these maps in the sagittal view.
If you find VoxHunt useful for your research, please consider citing our paper:
@Article{,
title = {Resolving organoid brain region identities by mapping single-cell genomic data to reference atlases},
author = {Jonas Simon Fleck and Fátima Sanchís-Calleja and Zhisong He and Margozata Santel and Michael James Boyle and J. Gray Camp and Barbara Treutlein},
journal = {Cell Stem Cell},
year = {2021},
url = {https://doi.org/10.1016/j.stem.2021.02.015},
doi = {10.1016/j.stem.2021.02.015}
}