/mtscatac_protocol

Repository for mtscATAC-seq protocol paper

Primary LanguageR

mtscATAC-seq computational protocol

Last updated: 23 June 2022

Installs

  • Processing of .bcl files to .fastq files and alignment/processing of single-cell data for mtscATAC-seq libraries are directly enabled by CellRanger-ATAC.
    • We recommend version 2.0+ of CellRanger-ATAC for computational efficiency of processing and aligning data.
  • The mgatk python package for efficient genotyping of mitochondrial variants can be installed via PyPi via the following command:
pip3 install mgatk

Execution

We recommend running CellRanger-ATAC and mgatk back to back for each sample. Here is a real example from a shell script that we use!

i="mtscatac_sampleID"
dir="mito_samples/fastq"
ref="software/refdata-cellranger-arc-GRCh38-2020-mtMask"

# Align cellranger-atac samples
cellranger-atac-2.0.0/bin/cellranger-atac count --sample $i --id "${i}_v2-hg38-mtMask" --fastqs $dir --localcores 16 --reference $ref 

# Now run mgatk
mgatk bcall -i "${i}_v2-hg38-mtMask/outs/possorted_bam.bam" -n "${i}_hg19_mask_mgatk" -o "${i}_v2-hg38-mtMask_mgatk" \
	-b "${i}_v2-hg38-mtMask/outs/filtered_peak_bc_matrix/barcodes.tsv" -bt CB -c 16 -jm 12000m -g rCRS -qc

These encompass steps 44 and 45 in the protocol.

Once these commands have been executed, the interpretation of output files in the two output directories are provided in the protocol on steps 46-50.

Resources

  • Pre-computed coordinate files for common reference genomes for blacklisting the NUMTs is available here.
    • A dedicated walk through on masking NUMTs is detailed here.
    • Follow steps 38-41 to generate a new NUMT mask bed file if your reference genome isn't currently supported.
  • Raw mtscATAC-seq in .fastq format is available from our GEO repository here.

Advanced analyses

  • We provide a complete code repository for advanced clonotype calling and integrated analyses from our 2020 paper available here.
  • The Seurat/Signac packages provide compatible interactive workflows for mtDNA variant analysis with mtscATAC-seq
  • Analyses to reproduce bone marrow linage analysis, including clone fate bias, is available here.

Contact