Giter VIP home page Giter VIP logo

bcd_nano_cutntag's Introduction

Multimodal chromatin profiling using nanobody-based single-cell CUT&Tag

Marek Bartosovic, Goncalo Castelo-Branco

Code repository related to preprint.

Multimodal chromatin profiling using nanobody-based single-cell CUT&Tag Marek Bartosovic, Gonçalo Castelo-Branco bioRxiv 2022.03.08.483459; doi: https://doi.org/10.1101/2022.03.08.483459

Data availability

Processed files - Seurat objects, fragments file (cellranger), bigwig tracks per cluster and .h5 matrices are available as supplementary files in the GEO repository https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE198467

Reproducing the analysis

Step 1: prepare environment

conda environment is provided in env/environment.yaml

conda env create -f env/environment.yaml

Additional package dependencies that need to be installed:

R
install.packages(c('argparse','ggplot2','funr','Signac','scales','Seurat','rmarkdown','mclust','GGally','BiocManager','patchwork','markdown','UpSetR','pheatmap','viridis','purrr','Rmagic','devtools','raster'))
BiocManager::install(c('ensembldb','EnsDb.Mmusculus.v79','GenomeInfoDb', 'GenomicRanges', 'IRanges', 'Rsamtools','BiocGenerics','rtracklayer','limma','slingshot','BiocGenerics', 'DelayedArray', 'DelayedMatrixStats','limma', 'S4Vectors', 'SingleCellExperiment','SummarizedExperiment', 'batchelor', 'Matrix.utils')))

Have seqtk installed in $PATH

https://github.com/lh3/seqtk

Install papermill for cli for jupyter notebooks

python3 -m pip install papermill

Step 2: Download the data

use fasterq-dump or alternative to download the fastq files

https://github.com/ncbi/sra-tools/wiki/HowTo:-fasterq-dump

GEO repository

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE198467

Step 3: Clone the github repo with analysis code

git clone https://github.com/mardzix/bcd_nano_CUTnTag

Step 4: Modify config

Change config/config.yaml and specify path to

  1. Absolute path to fastq files
  2. Specific path to the tmp folder (Create any folder, e.g. in home)
  3. Specify conda environment to use
  4. Modify path to cellranger reference

Step 5: Run the pipeline

Pipeline is implemented in workflow management software Snakemake

Change the cluster-specific profile to your preference (e.g. slurm, condor etc.), or run without profile

For some example profiles see:

snakemake --snakefile code/workflow/Snakefile_single_modality.smk  --cores 16 --profile htcondor -p                                                                              

Johannes Köster, Sven Rahmann, Snakemake—a scalable bioinformatics workflow engine, Bioinformatics, Volume 28, Issue 19, 1 October 2012, Pages 2520–2522, https://doi.org/10.1093/bioinformatics/bts480

Snakemake

bcd_nano_cutntag's People

Contributors

mardzix avatar

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.