-
Notifications
You must be signed in to change notification settings - Fork 4
Installation
No installation of the pipeline is required. Everything is downloaded, unzipped and installed for you.
You need a working environment with:
- python (version 3.8.5)
- snakemake (version 5.32.2)
- singularity (version 3.6.3)
It is possible to create this environment with conda by:
source /mnt/beegfs/software/conda/etc/profile.d/conda.sh
conda env create -f /mnt/beegfs/pipelines/single-cell/<version>/envs/conda/single-cell_user.yaml --prefix="<path/MyNameEnvir>"
NB: singularity is not in this conda environment because of a namespace issue for slurm, but singularity is available by module load
.
To activate the complete environnement:
source /mnt/beegfs/software/conda/etc/profile.d/conda.sh
conda activate <path/MyNameEnvir>
module load singularity
To download this workflow you can use 2 methods:
- install git (see: https://git-scm.com/book/en/v2/Getting-Started-Installing-Git)
git clone https://github.com/gustaveroussy/single-cell.git
- go in single-cell/resources/WHITELISTS/ and unzip the 3M-february-2018.txt.gz file:
cd single-cell/resources/WHITELISTS/
gunzip 3M-february-2018.txt.gz
- donwload singularity environments:
cd ../../envs/singularity/ #from WHITELISTS folder
wget "https://zenodo.org/record/5769147/files/single_cell.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_DE.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_integration.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_oldcerebro.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_TCR_BCR.simg?download=1"
Or:
- download this workflow .
- unzip the dowloaded file
- go in single-cell/resources/WHITELISTS/ and unzip the 3M-february-2018.txt.gz file:
cd single-cell/resources/WHITELISTS/
gunzip 3M-february-2018.txt.gz
- donwload singularity environments:
cd ../../envs/singularity/ #from WHITELISTS folder
wget "https://zenodo.org/record/5769147/files/single_cell.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_DE.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_integration.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_oldcerebro.simg?download=1"
wget "https://zenodo.org/record/5769147/files/single_cell_TCR_BCR.simg?download=1"
Note: you can also recompile singularity image with singularity tool (singularity can be installed thanks to the "Installation of working environment" next section):
cd ../../envs/singularity/ #from WHITELISTS folder
singularity build single_cell.simg single_cell.def
singularity build single_cell_TCR_BCR.simg single_cell_TCR_BCR.def
singularity build single_cell_oldcerebro.simg single_cell_oldcerebro.def
You need a working environment with:
- python (version 3.8.5)
- snakemake (version 5.32.2)
- singularity (version 3.6.3)
It is possible to create this environment with conda:
- install conda (see: https://conda.io/projects/conda/en/latest/user-guide/install/index.html)
- install the environment by:
conda env create -f /mnt/beegfs/pipelines/single-cell/envs/conda/single-cell_user_local.yaml -n <MyNameEnvir>
To activate the complete environnement:
conda activate <MyNameEnvir>
Resources of the Theory of single cell RNA-seq
v1.3
Pipeline details
Configuration
-
Parameter file
- Steps
- Alignment_countTable_GE
- Droplets_QC_GE
- Filtering_GE
- Norm_DimRed_Eval_GE
- Clust_Markers_Annot_GE
- Cerebro
- Alignment_countTable_ADT
- Adding_ADT
- Alignment_annotations_TCR_BCR
- Adding_TCR
- Adding_BCR
- Int_Norm_DimRed_Eval_GE
- Int_Clust_Markers_Annot_GE
- Int_Adding_ADT
- Int_Adding_TCR
- Int_Adding_BCR
- Grp_Norm_DimRed_Eval_GE
- Grp_Clust_Markers_Annot_GE
- Grp_Adding_ADT
- Grp_Adding_TCR
- Grp_Adding_BCR
- Additional files
Results help
- Arborescence of all results
-
Observations and weird results
- Not a threshold by emptyDrops
- Large and small cells into the same sample
- emptyDrops does't work well
- More than 15% mitochondrial RNA while I filtered them out at 15%
- Impact of empty droplets on umap
- Choose the right number of dimensions
- Be careful with the colors, they are sometimes misleading
- Impact of bias correction on umap
Complete Examples of school cases
Individual analysis :
1 sample (scRNA-seq + ADT + TCR + BCR)
Grouped/Integrated analysis :
2 samples (scRNA-seq + ADT + TCR + BCR)