Toward data-driven genome breeding (digital breeding), we are developing data analysis infrastructure technology essential for genome editing, focusing on gene function analysis using bioinformatics called BioDX.

This is a project specific guide for the Bioiversity Genomics Europe (BGE project use of WorkflowHub.

Gene_fetch

This tool fetches gene sequences from NCBI databases based on taxonomy IDs (taxids) or taxonomic information. It can retrieve both protein and nucleotide sequences for various genes, including protein-coding genes (e.g., cox1, cytb, rbcl, matk) and rRNA genes (e.g., 16S, 18S).

Feature highlight

• Fetch protein and/or nucleotide sequences from NCBI GenBank database.
• Handles both direct nucleotide sequences and protein-linked nucleotide searches (CDS extraction includes fallback ...

gSpreadComp: Streamlining Microbial Community Analysis for Resistance, Virulence, and Plasmid-Mediated Spread

Overview

gSpreadComp is a UNIX-based, modular bioinformatics toolkit designed to streamline comparative genomics for analyzing microbial communities. It integrates genome annotation, gene spread calculation, plasmid-mediated horizontal gene transfer (HGT) detection and resistance-virulence ranking within the analysed microbial community to help researchers identify potential ...

sanger-tol/curationpretext

[![Cite with ...

[![Cite ...

Phasing and imputation pipeline

skim2mito

skim2mito is a snakemake pipeline for the batch assembly, annotation, and phylogenetic analysis of mitochondrial genomes from low coverage genome skims. The pipeline was designed to work with sequence data from museum collections. However, it should also work with genome skims from recently collected samples.

Contents

• Setup
• Example data
• Input
• Output
• Filtering contaminants
• [Assembly and ...

covid-sequence-analysis-workflow

This is the official repository of the SARS-CoV-2 variant surveillance pipeline developed by Danish Technical University (DTU), Eotvos Lorand University (ELTE), EMBL-EBI, Erasmus Medical Center (EMC) under the Versatile Emerging infectious disease Observatory (VEO) project. The project consists of 20 European partners. It is funded by the European Commission.

The ...

[![Cite ...

Assembly Evaluation for ERGA-BGE Reports

One Assembly, HiFi WGS reads + HiC reads

The workflow requires the following:

• Species Taxonomy ID number
• NCBI Genome assembly accession code
• BUSCO Lineage
• WGS accurate reads accession code
• NCBI HiC reads accession code

The workflow will get the data and process it to generate genome profiling (genomescope, smudgeplot -optional-), assembly stats (gfastats), merqury stats (QV, completeness), BUSCO, snailplot, contamination blobplot, and HiC ...

Assembly Evaluation for ERGA-BGE Reports

One Assembly, Illumina WGS reads + HiC reads

The workflow requires the following:

• Species Taxonomy ID number
• NCBI Genome assembly accession code
• BUSCO Lineage
• WGS accurate reads accession code
• NCBI HiC reads accession code

The workflow will get the data and process it to generate genome profiling (genomescope, smudgeplot -optional-), assembly stats (gfastats), merqury stats (QV, completeness), BUSCO, snailplot, contamination blobplot, and ...

The workflow requires the user to provide:

• ENSEMBL link address of the annotation GFF3 file
• ENSEMBL link address of the assembly FASTA file
• NCBI taxonomy ID
• BUSCO lineage
• OMArk database

Thw workflow will produce statistics of the annotation based on AGAT, BUSCO and OMArk.

cfDNA UniFlow is a unified, standardized, and ready-to-use workflow for processing whole genome sequencing (WGS) cfDNA samples from liquid biopsies. It includes essential steps for pre-processing raw cfDNA samples, quality control and reporting. Additionally, several optional utility functions like GC bias correction and estimation of copy number state are included. Finally, we provide specialized methods for extracting coverage derived signals and visualizations comparing cases and controls. ...

deepconsensus 1.2 snakemake pipeline

This snakemake-based workflow takes in a subreads.bam and results in a deepconsensus.fastq

• no methylation calls !

The metadata id of the subreads file needs to be: "m[numeric][numeric][numeric].subreads.bam"

Chunking (how many subjobs) and ccs min quality filter can be adjusted in the config.yaml

the checkpoint model for deepconsensus1.2 should be accessible like this: gsutil cp -r gs://brain-genomics-public/research/deepconsensus/models/v1.2/model_checkpoint/* ...

Annotation of an assembled bacterial genomes to detect genes, potential plasmids, integrons and Insertion sequence (IS) elements.

Short paired-end read analysis to provide quality analysis, read cleaning and taxonomy assignation

Antimicrobial resistance gene detection from assembled bacterial genomes

Assembly of bacterial paired-end short read data with generation of quality metrics and reports

Variant Interpretation Pipeline (VIP) that annotates, filters and reports prioritized causal variants in humans, see https://github.com/molgenis/vip for more information.

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

Pipeline for the identification of circular DNAs

A mapping-based pipeline for creating a phylogeny from bacterial whole genome sequences

GSC (Genotype Sparse Compression)

Genotype Sparse Compression (GSC) is an advanced tool for lossless compression of VCF files, designed to efficiently store and manage VCF files in a compressed format. It accepts VCF/BCF files as input and utilizes advanced compression techniques to significantly reduce storage requirements while ensuring fast query capabilities. In our study, we successfully compressed the VCF files from the 1000 Genomes Project (1000Gpip3), consisting of 2504 samples and 80 ...

GSC (Genotype Sparse Compression)

Genotype Sparse Compression (GSC) is an advanced tool for lossless compression of VCF files, designed to efficiently store and manage VCF files in a compressed format. It accepts VCF/BCF files as input and utilizes advanced compression techniques to significantly reduce storage requirements while ensuring fast query capabilities. In our study, we successfully compressed the VCF files from the 1000 Genomes Project (1000Gpip3), consisting of 2504 samples and 80 ...

Parabricks-Genomics-nf is a GPU-enabled pipeline for alignment and germline short variant calling for short read sequencing data. The pipeline utilises NVIDIA's Clara Parabricks toolkit to dramatically speed up the execution of best practice bioinformatics tools. Currently, this pipeline is configured specifically for NCI's Gadi HPC.

NVIDIA's Clara Parabricks can deliver a significant ...

...

HiC contact map generation

Snakemake pipeline for the generation of .pretext and .mcool files for visualisation of HiC contact maps with the softwares PretextView and HiGlass, respectively.

Prerequisites

This pipeine has been tested using Snakemake v7.32.4 and requires conda for installation of required tools. To run the pipline use the command:

snakemake --use-conda

There are provided a set of configuration and running scripts for exectution on a slurm queueing system. After configuring ...

This is a Nextflow implementaion of the GATK Somatic Short Variant Calling workflow. This workflow can be used to discover somatic short variants (SNVs and indels) from tumour and matched normal BAM files following GATK's Best Practices Workflow. The workflowis currently optimised to run efficiently and at scale on the National Compute Infrastructure, Gadi.

Workflow for Creating a large disease network from various datasets and databases for IBM, and applying the active subnetwork identification method MOGAMUN.

ANNOTATO - Annotation workflow To Annotate Them Oll

• ANNOTATO - Annotation workflow To Annotate Them Oll
• Overview of the workflow
• Input data
• Pipeline steps
• Output data
• Prerequisites
• Installation
• Running ANNOTATO
• Before running the pipeline (IMPORTANT) ...

ERGA Protein-coding gene annotation workflow.

Adapted from the work of Sagane Joye:

https://github.com/sdind/genome_annotation_workflow

Prerequisites

The following programs are required to run the workflow and the listed version were tested. It should be noted that older versions of snakemake are not compatible with newer versions of singularity as is noted here: https://github.com/nextflow-io/nextflow/issues/1659.

conda v 23.7.3 ...

CLAWS (CNAG's Long-read Assembly Workflow in Snakemake)

Snakemake Pipeline used for de novo genome assembly @CNAG. It has been developed for Snakemake v6.0.5.

It accepts Oxford Nanopore Technologies (ONT) reads, PacBio HFi reads, illumina paired-end data, illumina 10X data and Hi-C reads. It does the preprocessing of the reads, assembly, polishing, purge_dups, scaffodling and different evaluation steps. By default it will preprocess the reads, run Flye + Hypo + purge_dups + yahs and evaluate ...

ARA (Automated Record Analysis) : An automatic pipeline for exploration of SRA datasets with sequences as a query

Requirements

• Docker

• Please checkout the Docker installation guide.

or

• Mamba package manager

• Please checkout the mamba or micromamba official installation guide.

• We prefer mamba over conda since it is faster and uses ...

prepareChIPs

This is a simple snakemake workflow template for preparing single-end ChIP-Seq data. The steps implemented are:

1. Download raw fastq files from SRA
2. Trim and Filter raw fastq files using AdapterRemoval
3. Align to the supplied genome using bowtie2
4. Deduplicate Alignments using Picard MarkDuplicates
5. Call Macs2 Peaks using macs2

A pdf of the rulegraph is available here

Full details for each step are given below. Any additional ...

A CWL-based pipeline for calling small germline variants, namely SNPs and small INDELs, by processing data from Whole-genome Sequencing (WGS) or Targeted Sequencing (e.g., Whole-exome sequencing; WES) experiments.

On the respective GitHub folder are available:

• The CWL wrappers and subworkflows for the workflow
• A pre-configured YAML template, based on validation analysis of publicly available HTS data

Briefly, the workflow performs the following steps:

1. Quality control of Illumina reads ...

A CWL-based pipeline for calling small germline variants, namely SNPs and small INDELs, by processing data from Whole-genome Sequencing (WGS) or Targeted Sequencing (e.g., Whole-exome sequencing; WES) experiments.

On the respective GitHub folder are available:

• The CWL wrappers and subworkflows for the workflow
• A pre-configured YAML template, based on validation analysis of publicly available HTS data

Briefly, the workflow performs the following steps:

1. Quality control of Illumina reads ...

Purge dups

This snakemake pipeline is designed to be run using as input a contig-level genome and pacbio reads. This pipeline has been tested with snakemake v7.32.4. Raw long-read sequencing files and the input contig genome assembly must be given in the config.yaml file. To execute the workflow run:

snakemake --use-conda --cores N

Or configure the cluster.json and run using the ./run_cluster command

MoMofy

Module for integrative Mobilome prediction

Bacteria can acquire genetic material through horizontal gene transfer, allowing them to rapidly adapt to changing environmental conditions. These mobile genetic elements can be classified into three main categories: plasmids, phages, and integrons. Autonomous elements are those capable of excising themselves from the chromosome, reintegrating elsewhere, and potentially modifying the host's physiology. Small integrative elements like insertion ...

IGVreport-nf

• Description
• Diagram
• User guide
• Workflow summaries
• Metadata
• Component tools
• Required (minimum) inputs/parameters
• Additional notes
• Help/FAQ/Troubleshooting
• Acknowledgements/citations/credits

Description

Quickly generate [IGV .html ...

GermlineStructuralV-nf is a pipeline for identifying structural variant events in human Illumina short read whole genome sequence data. GermlineStructuralV-nf identifies structural variant and copy number events from BAM files using Manta, Smoove, and TIDDIT. Variants are then merged using SURVIVOR, ...

IndexReferenceFasta-nf

===========

• Description
• Diagram
• User guide
• Benchmarking
• Workflow summaries
• Metadata
• Component tools
• Required (minimum) inputs/parameters
• Additional notes
• Help/FAQ/Troubleshooting
• Acknowledgements/citations/credits ...

Workflow for Metagenomics from bins to metabolic models (GEMs)

Summary

• Prodigal gene prediction
• CarveMe genome scale metabolic model reconstruction
• MEMOTE for metabolic model testing
• SMETANA Species METabolic interaction ANAlysis

Other UNLOCK workflows on WorkflowHub: https://workflowhub.eu/projects/16/workflows?view=default

All tool CWL files and other workflows can be found here: Tools: https://gitlab.com/m-unlock/cwl Workflows: https://gitlab.com/m-unlock/cwl/workflows

**How ...

Workflow for LongRead Quality Control and Filtering

• NanoPlot (read quality control) before and after filtering
• Filtlong (read trimming)
• Kraken2 taxonomic read classification before and after filtering
• Minimap2 read filtering based on given references

Other UNLOCK workflows on WorkflowHub: https://workflowhub.eu/projects/16/workflows?view=default

All tool CWL files and other workflows can be found here: https://gitlab.com/m-unlock/cwl/workflows

**How to setup and use an UNLOCK ...

- Deprecated -

See our updated hybrid assembly workflow: https://workflowhub.eu/workflows/367

And other workflows: https://workflowhub.eu/projects/16#workflows

Workflow for sequencing with ONT Nanopore data, from basecalled reads to (meta)assembly and binning

• Workflow Nanopore Quality
• Kraken2 taxonomic classification of FASTQ reads
• Flye (de-novo assembly)
• Medaka (assembly polishing)
• metaQUAST (assembly quality reports)

When Illumina reads are provided:

• Workflow ...

Workflow for Illumina Quality Control and Filtering

Multiple paired datasets will be merged into single paired dataset.

Summary:

• FastQC on raw data files
• fastp for read quality trimming
• BBduk for phiX and (optional) rRNA filtering
• Kraken2 for taxonomic classification of reads (optional)
• BBmap for (contamination) filtering using given references (optional)
• FastQC on filtered (merged) data

Other UNLOCK workflows on WorkflowHub: https://workflowhub.eu/projects/16/workflows?view=default ...

...

Bootstrapping-for-BQSR @ NCI-Gadi is a pipeline for bootstrapping a variant resource to enable GATK base quality score recalibration (BQSR) for non-model organisms that lack a publicly available variant resource. This implementation is optimised for the National Compute Infrastucture's Gadi HPC. Multiple rounds of bootstrapping can be performed. Users can use Fastq-to-bam @ NCI-Gadi and Germline-ShortV @ NCI-Gadi to ...

Local Cromwell implementation of GATK4 germline variant calling pipeline

See the GATK website for more information on this toolset

Assumptions

• Using hg38 human reference genome build
• Running 'locally' i.e. not using HPC/SLURM scheduling, or containers. This repo was specifically tested on Pawsey Nimbus 16 CPU, 64GB RAM virtual machine, primarily running in the /data volume storage partition.
• Starting from short-read Illumina paired-end fastq ...

Fastq-to-BAM @ NCI-Gadi is a genome alignment workflow that takes raw FASTQ files, aligns them to a reference genome and outputs analysis ready BAM files. This workflow is designed for the National Computational Infrastructure's (NCI) Gadi supercompter, leveraging multiple nodes on NCI Gadi to run all stages of the workflow in parallel, either massively parallel using the scatter-gather approach or parallel by sample. It consists of a number of stages and follows the BROAD Institute's best practice ...

SLURM HPC Cromwell implementation of GATK4 germline variant calling pipeline

See the GATK website for more information on this toolset

Assumptions

• Using hg38 human reference genome build
• Running using HPC/SLURM scheduling. This repo was specifically tested on Pawsey Zeus machine, primarily running in the /scratch partition.
• Starting from short-read Illumina paired-end fastq files as input

Dependencies

The following versions have been ...

Germline-ShortV @ NCI-Gadi is an implementation of the BROAD Institute's best practice workflow for germline short variant discovery. This implementation is optimised for the National Compute Infrastucture's Gadi HPC, utilising scatter-gather parallelism to enable use of multiple nodes with high CPU or memory efficiency. This workflow requires sample BAM files, which can be generated using the Fastq-to-bam @ NCI-Gadi pipeline. Germline-ShortV can be applied ...

ORSON combine state-of-the-art tools for annotation processes within a Nextflow pipeline: sequence similarity search (PLAST, BLAST or Diamond), functional annotation retrieval (BeeDeeM) and functional prediction (InterProScan). When required, BUSCO completness evaluation and eggNOG Orthogroup annotation can be activated. While ORSON results can be analyzed through the command-line, it also offers the possibility to be compatible with BlastViewer or Blast2GO graphical tools.

Pipelines used by the genomes assembly teams part of the Biodiversity Genomics Europe project

https://biodiversitygenomics.eu/

Collection of Galaxy workflows for generating results used for creating ERGA-BGE Reports

For a given genome, two workflows should be run: the assembly evaluation (ASM analyses), and the annotation evaluation (ANNOT analyses)

Depending on the kind of data used for the genome assembly, you should choose HiFi or ONT (Illumina) workflows for ASM analyses

Collection of de-novo genome assembly workflows written for implementation in Galaxy

Input data should be Oxford Nanopore raw reads plus Illumina WGS reads and Illumina 3-dimensional Chromatin Confirmation Capture (HiC) reads

Executing all workflows will output one scaffolded collapsed assembly and the complete QC analyses

Please run the workflows in order: WF0 (there are two, one for ONT, and another one for Illumina that can be used independently for the WGS and HiC reads), WF1, WF2, WF3, WF4

Collection of de-novo genome assembly workflows written for implementation in Galaxy

Input data should be Oxford Nanopore raw reads plus Illumina WGS reads and Illumina 3-dimensional Chromatin Confirmation Capture (HiC) reads

Executing all workflows will output one scaffolded collapsed assembly and the complete QC analyses

Please run the workflows in order: WF0 (there are two, one for ONT, and another one for Illumina that can be used independently for the WGS and HiC reads), WF1, WF2, WF3, WF4

This is a general collection of workflows used by or developed by members of the BGE project.

Collection of de-novo genome assembly workflows written for implementation in Galaxy

Input data should be PacBio HiFi reads and Illumina 3-dimensional Chromatin Confirmation Capture (HiC) reads

Executing all workflows will output two scaffolded haplotype assemblies and the complete QC analyses

Please run the workflows in order: WF0 (there are two, one for HiFi and one for Illumina HiC), WF1, WF2, WF3, WF4

The Vertebrate Genomes Pipelines in Galaxy are intended to allow a user to generate high-quality near error-free assemblies of species from a user's own data or from the GenomeArk database.