{"title":"ChoCallate: a computational pipeline for consensus germline variant calling in genomic sequencing data.","authors":"Aleksey Ermolaev, Mikhail Divashuk","doi":"10.7717/peerj.21758","DOIUrl":null,"url":null,"abstract":"<p><p>Accurate germline variant calling in non-model and polyploid organisms remains challenging due to sequencing errors, genomic complexity, and caller-specific biases. Existing ensemble methods of variant calling are constrained by species-specific training, limited portability, or an inability to process multiple samples concurrently. To address these limitations we developed ChoCallate, a Nextflow-based pipeline for consensus germline variant calling. The pipeline integrates three species-independent callers, bcftools, FreeBayes, and Genome Analysis Toolkit (GATK) HaplotypeCaller, within a flexible voting mechanism. A distinctive feature of ChoCallate is its incorporation of coverage-aware consensus generation <i>via</i> Browser Extensible Data (BED) file, enabling reliable variant identification even in reduced-representation sequencing data such as Restriction site Associated DNA sequencing (RAD-seq) or RNA-seq. In a simulated dataset containing 10,000 single-nucleotide polymorphisms (SNPs) and 1,000 InDels of <i>Arabidopsis thaliana</i> across coverage depths ranging from 3× to 50×, the two-vote consensus strategy provided a favourable precision-recall balance, particularly at lower coverage (3×-10×), where SNP precision was enhanced without substantial loss of sensitivity. False-negative variants were found to concentrate in centromeric and nucleolus organizer regions (NORs) across all strategies, indicating that repetitive genomic structure, rather than consensus logic, constitutes the primary limitation. In a case study of 200 bread wheat (<i>Triticum aestivum</i>) samples, ChoCallate identified 11,178 high-confidence variants, and a genome-wide association study (GWAS) for awnedness successfully recovered the known B1/ALI-1 locus on chromosome 5A, replicating prior findings with reduced genomic inflation (<i>λ</i> (gc) = 1.06) by decreasing the number of false-positive variants through a consensus-based approach. ChoCallate is portable, species-agnostic, and freely available.</p>","PeriodicalId":19799,"journal":{"name":"PeerJ","volume":"14 ","pages":"e21758"},"PeriodicalIF":2.9000,"publicationDate":"2026-09-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13632679/pdf/","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"PeerJ","FirstCategoryId":"99","ListUrlMain":"https://doi.org/10.7717/peerj.21758","RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"2026/1/1 0:00:00","PubModel":"eCollection","JCR":"Q2","JCRName":"MULTIDISCIPLINARY SCIENCES","Score":null,"Total":0}
引用次数: 0
Abstract
Accurate germline variant calling in non-model and polyploid organisms remains challenging due to sequencing errors, genomic complexity, and caller-specific biases. Existing ensemble methods of variant calling are constrained by species-specific training, limited portability, or an inability to process multiple samples concurrently. To address these limitations we developed ChoCallate, a Nextflow-based pipeline for consensus germline variant calling. The pipeline integrates three species-independent callers, bcftools, FreeBayes, and Genome Analysis Toolkit (GATK) HaplotypeCaller, within a flexible voting mechanism. A distinctive feature of ChoCallate is its incorporation of coverage-aware consensus generation via Browser Extensible Data (BED) file, enabling reliable variant identification even in reduced-representation sequencing data such as Restriction site Associated DNA sequencing (RAD-seq) or RNA-seq. In a simulated dataset containing 10,000 single-nucleotide polymorphisms (SNPs) and 1,000 InDels of Arabidopsis thaliana across coverage depths ranging from 3× to 50×, the two-vote consensus strategy provided a favourable precision-recall balance, particularly at lower coverage (3×-10×), where SNP precision was enhanced without substantial loss of sensitivity. False-negative variants were found to concentrate in centromeric and nucleolus organizer regions (NORs) across all strategies, indicating that repetitive genomic structure, rather than consensus logic, constitutes the primary limitation. In a case study of 200 bread wheat (Triticum aestivum) samples, ChoCallate identified 11,178 high-confidence variants, and a genome-wide association study (GWAS) for awnedness successfully recovered the known B1/ALI-1 locus on chromosome 5A, replicating prior findings with reduced genomic inflation (λ (gc) = 1.06) by decreasing the number of false-positive variants through a consensus-based approach. ChoCallate is portable, species-agnostic, and freely available.
期刊介绍:
PeerJ is an open access peer-reviewed scientific journal covering research in the biological and medical sciences. At PeerJ, authors take out a lifetime publication plan (for as little as $99) which allows them to publish articles in the journal for free, forever. PeerJ has 5 Nobel Prize Winners on the Board; they have won several industry and media awards; and they are widely recognized as being one of the most interesting recent developments in academic publishing.