sanger-tol/genomeassembly is a bioinformatics pipeline for a genome assembly for HiFi, Illumina 10x (optional), and HiC data. It performs the following steps: raw assembly, purging from haplotigs, optional polishing, and scaffolding.
Original assembly of HiFi reads is performed using hifiasm assembler in two modes - original and using HiC data (optional). Then assembly is purged from alternative haplotigs using purge_dups. Next optional step is polishing of the purged assembly using Illumina 10X read sequencing. 10X reads are mapped to the full assembly (purged + haplotigs) using Longranger and polishing is implemented using Freebayes. HiC reads are further mapped with bwamem2 to the primary contigs, which are further scaffolded with YaHS using the provided Hi-C data. Polished and scaffolded assemblies are evaluated using GFASTATS, BUSCO and MERQURY.FK
The pipeline is built using Nextflow, a workflow tool to run tasks across multiple compute infrastructures in a very portable manner. It uses Docker/Singularity containers making installation trivial and results highly reproducible. The Nextflow DSL2 implementation of this pipeline uses one container per process which makes it much easier to maintain and update software dependencies. Where possible, these processes have been submitted to and installed from nf-core/modules in order to make them available to all nf-core pipelines, and to everyone within the Nextflow community!
On release, automated continuous integration tests run the pipeline on a full-sized dataset on the LSF infrastructure. This ensures that the pipeline runs on LSF, has sensible resource allocation defaults set to run on real-world datasets, and permits the persistent storage of results to benchmark between pipeline releases and other analysis sources.
While the steps are described in a sequential order, many of them can be executed as parallel jobs.
hifiasm_hic_on
option is set
polishing_on
organelles_on
Note If you are new to Nextflow and nf-core, please refer to this page on how to set-up Nextflow. Make sure to test your setup with
-profile test
before running the workflow on actual data.
Currently, it is advised to run the pipeline with docker or singularity as some modules do not have a conda env associated with them.
To run the pipeline use a command-line:
nextflow run sanger-tol/genomeassembly -profile singularity,YOURPROFILE --outdir <OUTDIR>
For more details on how to run the pipeline and interprete the results see usage and output sections of the documentation.
sanger-tol/genomeassembly was originally written by @ksenia-krasheninnikova.
We thank the following people for their extensive assistance in the development of this pipeline:
@priyanka-surana for the code review, very helpful coding suggestions, and assistance with pushing this pipeline forward through development.
@mcshane and @c-zhou for the design and implementation of the original pipelines for purging (@mcshane), polishing (@mcshane) and scaffolding (@c-zhou).
TreeVal team Damon-Lee Pointon (@DLBPointon), Yumi Sims (@yumisims) and William Eagles (@weaglesBio) for implementation of the hic-mapping pipeline.
@muffato for help with nf-core integration, dealing with infrastructure and troubleshooting, for the code reviews and valuable suggestions at the different stages of the pipeline development.
@gq1 for the code review, valuable suggestions to the code improvement and contributions to the full test setup.
@mahesh-panchal for nextflow implementation of the purging pipeline, code review and valuable suggestions to the nf-core modules implementation.
If you would like to contribute to this pipeline, please see the contributing guidelines.
If you use sanger-tol/genomeassembly for your analysis, please cite it using the following doi: 10.5281/zenodo.10391851.
An extensive list of references for the tools used by the pipeline can be found in the CITATIONS.md
file.
This pipeline uses code and infrastructure developed and maintained by the nf-core community, reused here under the MIT license.
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.