Code to replicate analyses by Brann et al.
- Create a new conda environment.
conda create -n dv_score python=3.9
- Activate that env
conda activate dv_score
. - Clone and enter this repo:
git clone git@github.com:dattalab/Brann_olfactory_dorsoventral.git && cd Brann_olfactory_dorsoventral
- Install the code in this directory via
pip install -e .
- To install the specific versions of packages used when this repo was created do
pip install -r requirements.txt
. The additional requirements for running the notebooks in this repo are topip install numpy seaborn scikit-learn jupyter notebook
. The scripts also rely onpip install pysam scanpy
.
- Processed data is available on the NCBI GEO at Accession number GSE173947 and raw fastq files can be found on the SRA (accession SRP318630).
- Data were preprocessed by running the Nextflow pipeline in the scripts folder.
- Additional instructions for how to work with the raw data and to work with the olfactory gene expression programs (GEPs) can be found in the follow GitHub repo: Tsukahara_Brann_OSN.
Code to generate key results, focusing on the dorsoventral (DV) score.
- Open a new jupyter notebook with
jupyter notebook
. - Run the notebooks.
- Additional stand-alone scripts demonstrate the Nextflow pipelines that were used to uniformly preprocess scRNA-seq data, as well the scVI and scANVI models that were used for data integration.
For more details, please post an issue here or contact the authors.