{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,6,18]],"date-time":"2026-06-18T07:12:13Z","timestamp":1781766733367,"version":"3.54.5"},"reference-count":69,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2023,11,7]],"date-time":"2023-11-07T00:00:00Z","timestamp":1699315200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Bioinform."],"abstract":"<jats:p><jats:bold>Background:<\/jats:bold> Automating data analysis pipelines is a key requirement to ensure reproducibility of results, especially when dealing with large volumes of data. Here we assembled automated pipelines for the analysis of High-throughput Sequencing (HTS) data originating from RNA-Seq, ChIP-Seq and Germline variant calling experiments. We implemented these workflows in Common workflow language (CWL) and evaluated their performance by: i) reproducing the results of two previously published studies on Chronic Lymphocytic Leukemia (CLL), and ii) analyzing whole genome sequencing data from four Genome in a Bottle Consortium (GIAB) samples, comparing the detected variants against their respective golden standard truth sets.<\/jats:p><jats:p><jats:bold>Findings:<\/jats:bold> We demonstrated that CWL-implemented workflows clearly achieved high accuracy in reproducing previously published results, discovering significant biomarkers and detecting germline SNP and small INDEL variants.<\/jats:p><jats:p><jats:bold>Conclusion:<\/jats:bold> CWL pipelines are characterized by reproducibility and reusability; combined with containerization, they provide the ability to overcome issues of software incompatibility and laborious configuration requirements. In addition, they are flexible and can be used immediately or adapted to the specific needs of an experiment or study. The CWL-based workflows developed in this study, along with version information for all software tools, are publicly available on GitHub (<jats:ext-link>https:\/\/github.com\/BiodataAnalysisGroup\/CWL_HTS_pipelines<\/jats:ext-link>) under the MIT License. They are suitable for the analysis of short-read (such as Illumina-based) data and constitute an open resource that can facilitate automation, reproducibility and cross-platform compatibility for standard bioinformatic analyses.<\/jats:p>","DOI":"10.3389\/fbinf.2023.1275593","type":"journal-article","created":{"date-parts":[[2023,11,7]],"date-time":"2023-11-07T09:01:30Z","timestamp":1699347690000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":1,"title":["Software pipelines for RNA-Seq, ChIP-Seq and germline variant calling analyses in common workflow language (CWL)"],"prefix":"10.3389","volume":"3","author":[{"given":"Konstantinos A.","family":"Kyritsis","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Nikolaos","family":"Pechlivanis","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Fotis","family":"Psomopoulos","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2023,11,7]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"W537","DOI":"10.1093\/nar\/gky379","article-title":"The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update","volume":"46","author":"Afgan","year":"2018","journal-title":"Nucleic Acids Res."},{"key":"B2","doi-asserted-by":"publisher","first-page":"9354","DOI":"10.1038\/s41598-019-45839-z","article-title":"The ENCODE blacklist: identification of problematic regions of the genome","volume":"9","author":"Amemiya","year":"2019","journal-title":"Sci. Rep."},{"key":"B3","volume-title":"FastQC","author":"Andrews","year":"2012"},{"key":"B4","doi-asserted-by":"publisher","first-page":"68","DOI":"10.1038\/nature15393","article-title":"A global reference for human genetic variation","volume":"526","author":"Auton","year":"2015","journal-title":"Nature"},{"key":"B5","doi-asserted-by":"publisher","first-page":"452","DOI":"10.1038\/533452a","article-title":"1,500 scientists lift the lid on reproducibility","volume":"533","author":"Baker","year":"2016","journal-title":"Nature"},{"key":"B6","doi-asserted-by":"publisher","first-page":"622","DOI":"10.1038\/s41597-022-01710-x","article-title":"Introducing the FAIR Principles for research software","volume":"9","author":"Barker","year":"2022","journal-title":"Sci. Data"},{"key":"B7","doi-asserted-by":"publisher","first-page":"2114","DOI":"10.1093\/bioinformatics\/btu170","article-title":"Trimmomatic: a flexible trimmer for Illumina sequence data","volume":"30","author":"Bolger","year":"2014","journal-title":"Bioinformatics"},{"key":"B8","doi-asserted-by":"publisher","first-page":"75","DOI":"10.3389\/fgene.2014.00075","article-title":"Impact of artifact removal on ChIP quality metrics in ChIP-seq and ChIP-exo data","volume":"5","author":"Carroll","year":"2014","journal-title":"Front. Genet."},{"key":"B9","volume-title":"Common Workflow Language, v1.0","author":"Chapman","year":"2016"},{"key":"B10","doi-asserted-by":"publisher","first-page":"62","DOI":"10.1145\/2812803","article-title":"Repeatability in computer systems research","volume":"59","author":"Collberg","year":"2016","journal-title":"Commun. ACM"},{"key":"B11","doi-asserted-by":"publisher","first-page":"54","DOI":"10.1109\/MCC.2016.100","article-title":"To docker or not to docker: a security perspective","volume":"3","author":"Combe","year":"2016","journal-title":"IEEE Cloud Comput."},{"key":"B12","doi-asserted-by":"publisher","first-page":"1840","DOI":"10.1182\/blood.v94.6.1840","article-title":"Ig V gene mutation status and CD38 expression as novel prognostic indicators in chronic lymphocytic leukemia","volume":"94","author":"Damle","year":"1999","journal-title":"Blood"},{"key":"B13","doi-asserted-by":"publisher","first-page":"giab008","DOI":"10.1093\/gigascience\/giab008","article-title":"Twelve years of SAMtools and BCFtools","volume":"10","author":"Danecek","year":"2021","journal-title":"GigaScience"},{"key":"B14","doi-asserted-by":"publisher","first-page":"491","DOI":"10.1038\/ng.806","article-title":"A framework for variation discovery and genotyping using next-generation DNA sequencing data","volume":"43","author":"DePristo","year":"2011","journal-title":"Nat. Genet."},{"key":"B15","doi-asserted-by":"publisher","first-page":"316","DOI":"10.1038\/nbt.3820","article-title":"Nextflow enables reproducible computational workflows","volume":"35","author":"Di Tommaso","year":"2017","journal-title":"Nat. Biotechnol."},{"key":"B16","doi-asserted-by":"publisher","first-page":"3047","DOI":"10.1093\/bioinformatics\/btw354","article-title":"MultiQC: summarize analysis results for multiple tools and samples in a single report","volume":"32","author":"Ewels","year":"2016","journal-title":"Bioinformatics"},{"key":"B17","doi-asserted-by":"publisher","first-page":"1371","DOI":"10.1182\/blood.v116.21.1371.1371","article-title":"Unmutated IGHV1-69\/D3-16\/J3 stereotyped HCDR3 rearrangements (subset 6) are associated with indolent disease course and have outcome independent of mutational status in early stage CLL (rai 0)","volume":"116","author":"Forconi","year":"2010","journal-title":"BLOOD"},{"key":"B18","doi-asserted-by":"publisher","first-page":"243","DOI":"10.1038\/nbt.3172","article-title":"Ballgown bridges the gap between transcriptome assembly and expression analysis","volume":"33","author":"Frazee","year":"2015","journal-title":"Nat. Biotechnol."},{"key":"B19","doi-asserted-by":"publisher","first-page":"108","DOI":"10.1162\/dint_a_00033","article-title":"FAIR computational workflows","volume":"2","author":"Goble","year":"2020","journal-title":"Data Intell."},{"key":"B20","volume-title":"Implementing FAIR digital objects in the EOSC-life workflow collaboratory","author":"Goble","year":"2021"},{"key":"B21","doi-asserted-by":"publisher","first-page":"84","DOI":"10.1016\/j.cpc.2018.05.021","article-title":"Enabling rootless Linux Containers in multi-user environments: the udocker tool","volume":"232","author":"Gomes","year":"2018","journal-title":"Comput. Phys. Commun."},{"key":"B22","doi-asserted-by":"publisher","first-page":"1848","DOI":"10.1182\/blood.v94.6.1848.418k05_1848_1854","article-title":"Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia","volume":"94","author":"Hamblin","year":"1999","journal-title":"Blood"},{"key":"B23","unstructured":"FASTX-Toolkit. FASTQ\/A short-reads pre-processing tools\n            HannonG.\n          2010"},{"key":"B24","doi-asserted-by":"publisher","first-page":"83","DOI":"10.1186\/s13059-017-1215-1","article-title":"Multi-omics approaches to disease","volume":"18","author":"Hasin","year":"2017","journal-title":"Genome Biol."},{"key":"B25","doi-asserted-by":"publisher","first-page":"907","DOI":"10.1038\/s41587-019-0201-4","article-title":"Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype","volume":"37","author":"Kim","year":"2019","journal-title":"Nat. Biotechnol."},{"key":"B26","doi-asserted-by":"publisher","first-page":"giz014","DOI":"10.1093\/gigascience\/giz014","article-title":"Common workflow language (CWL)-based software pipeline for de novo genome assembly from long- and short-read data","volume":"8","author":"Korhonen","year":"2019","journal-title":"GigaScience"},{"key":"B27","doi-asserted-by":"publisher","first-page":"1125","DOI":"10.1080\/15592294.2019.1633867","article-title":"EZH2 upregulates the PI3K\/AKT pathway through IGF1R and MYC in clinically aggressive chronic lymphocytic leukaemia","volume":"14","author":"Kosalai","year":"2019","journal-title":"Epigenetics"},{"key":"B28","doi-asserted-by":"publisher","first-page":"2520","DOI":"10.1093\/bioinformatics\/bts480","article-title":"Snakemake\u2014a scalable bioinformatics workflow engine","volume":"28","author":"K\u00f6ster","year":"2012","journal-title":"Bioinformatics"},{"key":"B29","doi-asserted-by":"publisher","first-page":"giz084","DOI":"10.1093\/gigascience\/giz084","article-title":"CWL-Airflow: a lightweight pipeline manager supporting Common Workflow Language","volume":"8","author":"Kotliar","year":"2019","journal-title":"GigaScience"},{"key":"B30","volume-title":"FelixKrueger\/TrimGalore: v0.6.7 - DOI via Zenodo (0.6.7)","author":"Krueger","year":"2021"},{"key":"B31","doi-asserted-by":"publisher","first-page":"555","DOI":"10.1038\/s41587-019-0054-x","article-title":"Best practices for benchmarking germline small-variant calls in human genomes","volume":"37","author":"Krusche","year":"2019","journal-title":"Nat. Biotechnol."},{"key":"B32","doi-asserted-by":"publisher","first-page":"e0177459","DOI":"10.1371\/journal.pone.0177459","article-title":"Singularity: scientific containers for mobility of compute","volume":"12","author":"Kurtzer","year":"2017","journal-title":"PLOS ONE"},{"key":"B33","volume-title":"Supporting data for \u201csoftware pipelines for RNA-seq, ChIP-seq and germline variant calling analyses in common Workflow Language (CWL)","author":"Kyritsis","year":"2023"},{"key":"B34","doi-asserted-by":"publisher","first-page":"525","DOI":"10.1038\/nature15395","article-title":"Mutations driving CLL and their evolution in progression and relapse","volume":"526","author":"Landau","year":"2015","journal-title":"Nature"},{"key":"B35","doi-asserted-by":"publisher","first-page":"D28","DOI":"10.1093\/nar\/gkq967","article-title":"The European Nucleotide archive","volume":"39","author":"Leinonen","year":"","journal-title":"Nucleic Acids Res."},{"key":"B36","doi-asserted-by":"publisher","first-page":"D19","DOI":"10.1093\/nar\/gkq1019","article-title":"The sequence read archive","volume":"39","author":"Leinonen","year":"","journal-title":"Nucleic Acids Res."},{"key":"B37","volume-title":"Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM","author":"Li","year":"2013"},{"key":"B38","doi-asserted-by":"publisher","first-page":"2078","DOI":"10.1093\/bioinformatics\/btp352","article-title":"The sequence alignment\/map format and SAMtools","volume":"25","author":"Li","year":"2009","journal-title":"Bioinformatics"},{"key":"B39","doi-asserted-by":"publisher","first-page":"923","DOI":"10.1093\/bioinformatics\/btt656","article-title":"featureCounts: an efficient general purpose program for assigning sequence reads to genomic features","volume":"30","author":"Liao","year":"2014","journal-title":"Bioinformatics"},{"key":"B40","doi-asserted-by":"publisher","first-page":"550","DOI":"10.1186\/s13059-014-0550-8","article-title":"Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2","volume":"15","author":"Love","year":"2014","journal-title":"Genome Biol."},{"key":"B41","doi-asserted-by":"publisher","first-page":"320","DOI":"10.1016\/j.cell.2013.03.036","article-title":"Selective inhibition of tumor oncogenes by disruption of super-enhancers","volume":"153","author":"Lov\u00e9n","year":"2013","journal-title":"Cell"},{"key":"B42","doi-asserted-by":"publisher","first-page":"1297","DOI":"10.1101\/gr.107524.110","article-title":"The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data","volume":"20","author":"McKenna","year":"2010","journal-title":"Genome Res."},{"key":"B43","article-title":"Docker: lightweight linux containers for consistent development and deployment","volume":"2","author":"Merkel","year":"2014","journal-title":"Linux J."},{"key":"B44","doi-asserted-by":"publisher","first-page":"53","DOI":"10.1186\/s13059-016-0917-0","article-title":"The real cost of sequencing: scaling computation to keep pace with data generation","volume":"17","author":"Muir","year":"2016","journal-title":"Genome Biol."},{"key":"B45","doi-asserted-by":"publisher","first-page":"2695","DOI":"10.1002\/ijc.31999","article-title":"Integrated epigenomic and transcriptomic analysis reveals TP63 as a novel player in clinically aggressive chronic lymphocytic leukemia","volume":"144","author":"Papakonstantinou","year":"2019","journal-title":"Int. J. Cancer"},{"key":"B46","doi-asserted-by":"publisher","first-page":"1650","DOI":"10.1038\/nprot.2016.095","article-title":"Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown","volume":"11","author":"Pertea","year":"2016","journal-title":"Nat. Protoc."},{"key":"B47","volume-title":"Broad Institute, GitHub repository","author":"Picard","year":"2021"},{"key":"B48","doi-asserted-by":"publisher","first-page":"201178","DOI":"10.1101\/201178","article-title":"Scaling accurate genetic variant discovery to tens of thousands of samples","author":"Poplin","year":"2018","journal-title":"bioRxiv"},{"key":"B49","doi-asserted-by":"publisher","first-page":"841","DOI":"10.1093\/bioinformatics\/btq033","article-title":"BEDTools: a flexible suite of utilities for comparing genomic features","volume":"26","author":"Quinlan","year":"2010","journal-title":"Bioinformatics"},{"key":"B50","doi-asserted-by":"publisher","first-page":"W160","DOI":"10.1093\/nar\/gkw257","article-title":"deepTools2: a next generation web server for deep-sequencing data analysis","volume":"44","author":"Ram\u00edrez","year":"2016","journal-title":"Nucleic Acids Res."},{"key":"B51","doi-asserted-by":"publisher","first-page":"279","DOI":"10.1111\/j.1365-2354.2004.00489.x","article-title":"The clinical and epidemiological burden of chronic lymphocytic leukaemia","volume":"13","author":"Redaelli","year":"2004","journal-title":"Eur. J. Cancer Care (Engl)"},{"key":"B52","doi-asserted-by":"publisher","first-page":"4415","DOI":"10.1158\/1078-0432.CCR-08-3266","article-title":"Stereotyped B-cell receptor is an independent risk factor of chronic lymphocytic leukemia transformation to Richter syndrome","volume":"15","author":"Rossi","year":"2009","journal-title":"Clin. Cancer Res."},{"key":"B53","doi-asserted-by":"publisher","first-page":"389","DOI":"10.1038\/nature10730","article-title":"Differential oestrogen receptor binding is associated with clinical outcome in breast cancer","volume":"481","author":"Ross-Innes","year":"2012","journal-title":"Nature"},{"key":"B54","doi-asserted-by":"publisher","first-page":"e1003285","DOI":"10.1371\/journal.pcbi.1003285","article-title":"Ten simple rules for reproducible computational research","volume":"9","author":"Sandve","year":"2013","journal-title":"PLOS Comput. Biol."},{"key":"B55","doi-asserted-by":"publisher","first-page":"136","DOI":"10.5731\/pdajpst.2016.006734","article-title":"Biocompute objects\u2014a step towards evaluation and validation of biomedical scientific computations","volume":"71","author":"Simonyan","year":"2017","journal-title":"PDA J. Pharm. Sci. Technol."},{"key":"B56","doi-asserted-by":"publisher","first-page":"97","DOI":"10.3233\/DS-210053","article-title":"Packaging research artefacts with RO-Crate","volume":"5","author":"Soiland-Reyes","year":"2022","journal-title":"Data Sci."},{"key":"B57","doi-asserted-by":"publisher","first-page":"2","DOI":"10.1038\/s41746-019-0079-z","article-title":"The reproducibility crisis in the age of digital medicine","volume":"2","author":"Stupple","year":"2019","journal-title":"npj Digit. Med."},{"key":"B58","doi-asserted-by":"publisher","first-page":"11","DOI":"10.1002\/0471250953.bi1110s43","article-title":"From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline","volume":"43","author":"Van der Auwera","year":"2013","journal-title":"Curr. Protoc. Bioinforma."},{"key":"B59","volume-title":"Genomics in the cloud: using docker, GATK, and WDL in terra","author":"Van der Auwera","year":"2020"},{"key":"B60","doi-asserted-by":"publisher","first-page":"21","DOI":"10.1146\/annurev-pharmtox-010818-021511","article-title":"Systems pharmacology: defining the interactions of drug combinations","volume":"59","author":"van Hasselt","year":"2019","journal-title":"Annu. Rev. Pharmacol. Toxicol."},{"key":"B61","doi-asserted-by":"publisher","first-page":"314","DOI":"10.1038\/nbt.3772","article-title":"Toil enables reproducible, open source, big biomedical data analyses","volume":"35","author":"Vivian","year":"2017","journal-title":"Nat. Biotechnol."},{"key":"B62","doi-asserted-by":"publisher","first-page":"1122","DOI":"10.1038\/s41592-021-01205-4","article-title":"DOME: recommendations for supervised machine learning validation in biology","volume":"18","author":"Walsh","year":"2021","journal-title":"Nat. Methods"},{"key":"B63","doi-asserted-by":"publisher","first-page":"e164","DOI":"10.1093\/nar\/gkq603","article-title":"ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data","volume":"38","author":"Wang","year":"2010","journal-title":"Nucleic Acids Res."},{"key":"B64","doi-asserted-by":"publisher","first-page":"1113","DOI":"10.1038\/ng.2764","article-title":"The cancer genome Atlas pan-cancer analysis project","volume":"45","author":"Weinstein","year":"2013","journal-title":"Nat. Genet."},{"key":"B65","doi-asserted-by":"publisher","first-page":"307","DOI":"10.1016\/j.cell.2013.03.035","article-title":"Master transcription factors and mediator establish super-enhancers at key cell identity genes","volume":"153","author":"Whyte","year":"2013","journal-title":"Cell"},{"key":"B66","doi-asserted-by":"publisher","first-page":"160018","DOI":"10.1038\/sdata.2016.18","article-title":"The FAIR Guiding Principles for scientific data management and stewardship","volume":"3","author":"Wilkinson","year":"2016","journal-title":"Sci. Data"},{"key":"B67","doi-asserted-by":"publisher","first-page":"R137","DOI":"10.1186\/gb-2008-9-9-r137","article-title":"Model-based analysis of ChIP-seq (MACS)","volume":"9","author":"Zhang","year":"2008","journal-title":"Genome Biol."},{"key":"B68","doi-asserted-by":"publisher","first-page":"160025","DOI":"10.1038\/sdata.2016.25","article-title":"Extensive sequencing of seven human genomes to characterize benchmark reference materials","volume":"3","author":"Zook","year":"2016","journal-title":"Sci. Data"},{"key":"B69","doi-asserted-by":"publisher","first-page":"246","DOI":"10.1038\/nbt.2835","article-title":"Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls","volume":"32","author":"Zook","year":"2014","journal-title":"Nat. Biotechnol."}],"container-title":["Frontiers in Bioinformatics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbinf.2023.1275593\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,11,7]],"date-time":"2023-11-07T09:01:38Z","timestamp":1699347698000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fbinf.2023.1275593\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,11,7]]},"references-count":69,"alternative-id":["10.3389\/fbinf.2023.1275593"],"URL":"https:\/\/doi.org\/10.3389\/fbinf.2023.1275593","relation":{},"ISSN":["2673-7647"],"issn-type":[{"value":"2673-7647","type":"electronic"}],"subject":[],"published":{"date-parts":[[2023,11,7]]},"article-number":"1275593"}}