An open-source analysis pipeline to detect germline or somatic variants from whole genome sequencing.

Introduction

Previously known as the Cancer Analysis Workflow (CAW), Sarek is a workflow tool designed to run analyses on WGS data from regular samples or tumour / normal pairs, including relapse samples if required.

It's built using Nextflow, a bioinformatics domain specific language for workflow building. Software dependencies are handled using Docker or Singularity - container technologies that provide excellent reproducibility and ease of use. Singularity has been designed specifically for high-performance computing environments. This means that although Sarek has been primarily designed for use with the Swedish UPPMAX HPC systems, it should be able to run on any system that supports these two tools.

Sarek was developed at the National Genomics Infastructure and National Bioinformatics Infastructure Sweden which are both platforms at SciLifeLab. It is listed on the Elixir - Tools and Data Services Registry.

Workflow steps

Sarek is built with several workflow scripts. A wrapper script contained within the repository makes it easy to run the different workflow scripts as a single job. To test your installation, follow the tests documentation.

Raw FastQ files or aligned BAM files (with or without realignment & recalibration) can be used as inputs. You can choose which variant callers to use, plus the pipeline is capable of accommodating additional variant calling software or CNV callers if required.

The worflow steps and tools used are as follows:

1. Preprocessing - main.nf (based on GATK best practices)
   • Read alignment
     • BWA
   • Read realignment and recalibration of short-read data
     • GATK
2. Germline variant calling - germlineVC.nf
   • SNVs and small indels
     • GATK HaplotyeCaller
     • Strelka
   • Structural variants
     • Manta
3. Somatic variant calling - somaticVC.nf (optional)
   • SNVs and small indels
     • MuTect1
     • MuTect2
     • Freebayes
     • Strelka
   • Structural variants
     • Manta
   • Sample heterogeneity, ploidy and CNVs
     • ASCAT
4. Annotation - annotate.nf (optional)
   • Variant annotation
     • SnpEff
     • VEP (Variant Effect Predictor)
5. Reporting - runMultiQC.nf
   • Reporting
     • MultiQC

Documentation

The Sarek pipeline comes with documentation in the doc/ directory:

1. Installation documentation
2. Installation documentation specific for UPPMAX rackham
3. Installation documentation specific for UPPMAX bianca
4. Tests documentation
5. Reference files documentation
6. Configuration and profiles documentation
7. Intervals documentation
8. Running the pipeline
9. Examples
10. TSV file documentation
11. Processes documentation
12. Documentation about containers
13. Documentation about building
14. More information about ASCAT
15. Folder structure

Contributions & 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 Gitter or contact us: maxime.garcia@scilifelab.se, szilveszter.juhos@scilifelab.se

CHANGELOG

• CHANGELOG

Authors

• Sebastian DiLorenzo
• Jesper Eisfeldt
• Phil Ewels
• Maxime Garcia
• Szilveszter Juhos
• Max Käller
• Malin Larsson
• Marcel Martin
• Björn Nystedt
• Pall Olason

---