Workflow: pipeline.cwl

Fetched 2026-10-08 21:25:02 GMT

The pipeline implemented in this workflow allows for the analysis of paired normal tumor NGS WES BAM files. It runs a comprehensive characterization of the data. In order to ensure the validity of the analysis some QC metrics are included in the pipeline as the normal/tumor genotype distance (to ensure that samples are properly paired) and quality metrics (Picard HsMetrics). Ethnicity is automatically inferred by using EthSEQ. The pipeline then calls SNVs, indels and copy number profiles and next infers purity and ploidy of the tumor sample by using CLONETv2. Upon ploidy and admixture correction, CLONETv2 generates corrected copy number and allele specific copy numbers (the number of copies of each parental allele) calls together with segment clonality (the fraction of the tumor that have a certain aberration). A second purity inference method based on SNVs (TPES) is also run on the data. As a final result the pipeline produces two folders: one with the analysis results (data) and the second that includes the logs of each tool (logs).

children parents
Workflow as SVG
  • Selected
  • Default Values
  • Nested Workflows
  • Tools
  • Inputs/Outputs

Inputs

ID Type Title Doc
threads Integer (Optional)

The number of parallel threads that will be used for computation.

sample_sex https://w3id.org/cwl/view/git/06e4a6f868651c8e26a88db518af1974159ecba9/cwl/types/sex_enum.yaml#sex_enum

The sex of the patient from which the sample was collected. Either \"m\" or \"f\".

log_to_file Boolean

If true, the output generated by each tool will be redirected to a file. Otherwise the output will be printed on the output.

bam_file_tumor File

The BAM file of the tumor sample.

create_reports Boolean

If true enables the creation of graphical reports. Only the tools that support this type of output will include such outputs files.

bam_file_normal File

The BAM file of the normal sample.

kit_bait_bed_file File

The BED file containing the regions captured by the capture kit.

spia_snps_vcf_file File

The file with the SNPs that are used by SPIA to compute the genotype distance.

vep_data_directory Directory

The folder where VEP can find the annotation database.

kit_target_bed_file File

The BED file containing the regions targeted by the capture kit.

ethseq_snps_gds_file File

The GDS file containing the model of the SNPs used for ethnicity inference.

ethseq_snps_vcf_file File

The file with the SNPs that are included in the EthSEQ model.

snps_in_kit_vcf_file File

The list of SNPs that are contained in the regions targeted by the sequencing kit.

kit_bait_interval_file File

The interval list file containing the regions that are captured by the kit.

kit_target_interval_file File

The interval list file containing the regions that are targeted by the kit.

reference_genome_fasta_file File

The file containing the reference genome in FASTA format.

vep_reference_genome_version String

The VEP version of the reference genome to use for annotation.

cnvkit_accessible_regions_bed File (Optional)

An optional file containing the regions of the genome that are accessible (meaning that can be sequenced).

Steps

ID Runs Label Doc
vep
../tools/vep.cwl (CommandLineTool)
Annotates the variants with affected gene and other annotations.

VEP determines the effect of variants (SNPs, insertions, deletions, CNVs or structural variants) on genes, transcripts, and protein sequence, as well as regulatory regions.

spia

Runs SPIA to calculate the genotype distance between normal and tumor samples starting from the pilueps of the SNPs.

tpes
clonet
../tools/clonet.cwl (CommandLineTool)
Infers tumor purity and ploidy and computes corrected copy numbers, allele specific copy numbers and clonality of segments and SNVs.

Analyzes genomic data from next-generation sequencing experiments. CLONETv2 offers a set of functions to compute allele specific copy number and clonality from segmented data leveraging heterozygous SNPs position pileups. The package also calculates the clonality of single nucleotide variants (SNVs) given read counts at mutated positions. The package has been developed at the laboratory of Computational and Functional Oncology, Department of CIBIO, University of Trento (Italy), under the supervision of prof Francesca Demichelis. References: Prandi et al. (2014) <doi:10.1186/s13059-014-0439-6>; Carreira et al. (2014) <doi:10.1126/scitranslmed.3009448>. Prandi et al. (2019) <https://doi.org/10.1002/cpbi.81>

cnvkit
cnvkit_to_seg.cwl (Workflow)

Runs segmentation using CVNKit and exports the copy number segments in a segfile.

ethseq

Creates a genotype file based on the pileup of the SNPs and then runs EthSEQ.

mutect2

Runs mutect2 variant calling and and variant filter.

merge_output
18b35a04f35d0ba81b162e175ce532ea (ExpressionTool)
pileup_tumor
../tools/pacbam_pileup.cwl (CommandLineTool)
Computes the pileup of BAM files at specified positions

PaCBAM is a C command line tool for the complete characterization of genomic regions and single nucleotide positions from next-generation sequencing data. PaCBAM implements a fast and scalable multi-core computational engine, generates exhaustive output files for downstream analysis, introduces an innovative on-the-fly read duplicates filtering strategy and provides comprehensive visual reports.

pileup_normal
../tools/pacbam_pileup.cwl (CommandLineTool)
Computes the pileup of BAM files at specified positions

PaCBAM is a C command line tool for the complete characterization of genomic regions and single nucleotide positions from next-generation sequencing data. PaCBAM implements a fast and scalable multi-core computational engine, generates exhaustive output files for downstream analysis, introduces an innovative on-the-fly read duplicates filtering strategy and provides comprehensive visual reports.

snvs_coverage
snvs_coverage.cwl (Workflow)

Creates a genotype file starting from the pileup of the SNPs and then runs EthSEQ.

hsmetrics_tumor
../tools/picard_hsmetrics.cwl (CommandLineTool)
Computes the capture statistics for hybrid capture sequencing data

Picard is a set of command line tools for manipulating high-throughput sequencing (HTS) data and formats such as SAM/BAM/CRAM and VCF. CollectHsMetrics is used to capture several metrics useful to verify the quality of target-capture sequencing experiments.

merge_log_files
6e86c1f5a85ded5fb439702a57b51de1 (ExpressionTool)
hsmetrics_normal
../tools/picard_hsmetrics.cwl (CommandLineTool)
Computes the capture statistics for hybrid capture sequencing data

Picard is a set of command line tools for manipulating high-throughput sequencing (HTS) data and formats such as SAM/BAM/CRAM and VCF. CollectHsMetrics is used to capture several metrics useful to verify the quality of target-capture sequencing experiments.

merge_vep_snvs_coverage
../tools/merge_vep_snvs_coverage.cwl (CommandLineTool)
Merges SNV coverage data with VEP annotation

Script that allows to merge a file (in tabular format) of variants annotated by VEP with the coverage information of the tumor and normal samples. In particular the script will add the coverage of reference and alternative bases and the allelic fractions at SNV positions in normal and tumor samples.

Outputs

ID Type Label Doc
output Directory (Optional)

The output directory produced by the pipeline.

log_files Directory (Optional)

The log file, if enabled, that captures the output produced by each tool.

output_map https://w3id.org/cwl/view/git/06e4a6f868651c8e26a88db518af1974159ecba9/cwl/types/pipeline_output_map.yaml#pipeline_output_map (Optional)

A data structure that will allow to easily access the outputs produced by the pipeline.

Permalink: https://w3id.org/cwl/view/git/06e4a6f868651c8e26a88db518af1974159ecba9/cwl/workflows/pipeline.cwl