Skip to content

A bioinformatics pipeline for SARS-CoV-2 phylogenetic analysis

License

Notifications You must be signed in to change notification settings

nhhaidee/scovtree

Repository files navigation

nhhaidee/scovtree

SARS-CoV-2 phylogenetic analysis pipeline.

Nextflow install with bioconda

Introduction

nhhaidee/scovtree is a bioinformatics pipeline for SARS-CoV-2 phylogenetic analysis. Given an input FASTA file with SARS-CoV-2 sequences, this workflow will generate a maximum-likelihood phylogenetic tree (using IQ-TREE from a MAFFT or Nextalign multiple sequence alignment) and interactive HTML tree visualization (shiptv).

This pipeline also allows you to visualize your sequences along with the most closely SARS-CoV-2 sequences from GISAID (if both the GISAID sequences and metadata .tar.xz files are provided). Amino acid mutations can also be determined using Nextclade and shown in the tree visualization.

The pipeline is built using Nextflow, a workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It comes with docker containers making installation trivial and results highly reproducible.

Quick Start

  1. Install nextflow

  2. Install any of Docker, Singularity for full pipeline reproducibility (please only use Conda as a last resort; see docs)

  3. Download the pipeline and test it on a minimal dataset with a single command:

    nextflow run nhhaidee/scovtree -profile test,<docker/singularity/conda>
  4. Start running your own analysis!

    • Typical command for phylogenetic analysis:

      nextflow run nhhaidee/scovtree -profile <docker/singularity/conda> \
          --input your-sars-cov-2-sequences.fasta
    • Typical command for phylogenetic analysis with GISAID SARS-CoV-2 data:

      nextflow run nhhaidee/scovtree -profile <docker/singularity/conda> \
          --input your-sars-cov-2-sequences.fasta \
          --gisaid_sequences sequences_fasta_2021_06_14.tar.xz \
          --gisaid_metadata metadata_tsv_2021_06_14.tar.xz

Credits

nhhaidee/scovtree was originally written by Hai Nguyen.

Contributors:

  • Peter Kruczkiewicz for workflow conceptualization and software development.
  • jts/ncov-tools for bin/align2alleles.py and bin/phylogenetic_tree_snps.r to create phylogenetic tree with SNPs highlighted using R ggtree

Contributions and Support

If you would like to contribute to this pipeline, please see the contributing guidelines.

For further information or help, don't hesitate to get in touch on the Slack #scovtree channel (you can join with this invite).

Citations

You can cite the nf-core publication as follows:

The nf-core framework for community-curated bioinformatics pipelines.

Philip Ewels, Alexander Peltzer, Sven Fillinger, Harshil Patel, Johannes Alneberg, Andreas Wilm, Maxime Ulysse Garcia, Paolo Di Tommaso & Sven Nahnsen.

Nat Biotechnol. 2020 Feb 13. doi: 10.1038/s41587-020-0439-x.

In addition, references of tools and data used in this pipeline are as follows:

License

Copyright 2021 Canadian Food Inspection Agency of Canada, Government of Canada.

Distributed under the MIT license.