Single-cell mutation identification via phylogenetic inference (SCIPhI) is a new approach to mutation detection in individual tumor cells by leveraging the evolutionary relationship among cells. SCIPhI, jointly calls mutations in each individual cell and estimates the tumor phylogeny on these cells. Employing a Markov Chain Monte Carlo scheme we robustly account for the various sources of noise in single-cell sequencing data. Our approach enables us to reliably call mutations in each single cell even in experiments with high dropout rates and missing data
SCIPhI is freely available under a GNU General Public License v3.0 at https://github/cbg-ETHZ/SCIPhI
SCIPhI has the following dependencies which need to be installed:
Boost >= 1.6.x
,SeqAn >= 2.3.2
,DLIB >= 19.9
zlib
Clone SCIPhI with submodules from the soucre code:
git clone --recurse-submodules https://github.com/cbg-ethz/SCIPhI.git
In order to install SCIPhI issue the following commands in the github directory:
autoreconf -vif
./configure --with-boost=BOOST_PATH SEQAN_INCLUDEDIR=SEQAN_PATH DLIB_INCLUDEDIR=DLIB_PATH
make
An alternative is to use CMake (tested with version >= 3.11.1). Simply issue the commands:
mkdir build
cd build
cmake ..
cmake --build . --config Release
Afterwards you will find SCIPhI in the build directory.
SCIPhI expects the sequencing information to be passed in form of the well known mpileup format (http://www.htslib.org/doc/samtools.html). In order to generate such a file you need to align your fastq files to a reference and post process the result (e.g., following the instuctions here: https://gatk.broadinstitute.org/hc/en-us/articles/360035535912-Data-pre-processing-for-variant-discovery).
In order to see all available options type
sciphi -h
Executing
sciphi -o result --in cellNames.txt --seed 42 example.mpileup
will run SCIPhI using the cell names provided in cellNames.txt (same order as in the mpileup file). Note that cellNames.txt is a tab delimited file with the cell name in the first column and a cell type identifier in the second column. The cell type can be either CT (tumor cell), CN (control normal cell), or BN (control bulk normal). Note that SCIPhI assumes a pileup against a reference and ignores positions with 'N' as reference.