Scholarlink

Background <p>Identity-by-descent (IBD), which describes recent genetic co-ancestry between pairs of genomes, is a fundamental concept in population genomics. It has been used to estimate genetic relatedness, detect selection signals, and understand population demography. The IBD detection method <i>hmmIBD</i> demonstrates high accuracy in inferring IBD segments between haploid genomes, including <i>Plasmodium falciparum</i>, and is widely used in malaria genomic surveillance. However, the current single-threaded implementation of <i>hmmIBD</i> does not utilize the full capacity of multi-processor computers, making it difficult to apply to large data sets, and does not accommodate non-uniform recombination rates across the genome.</p> Methods <p>We developed an enhanced implementation of <i>hmmIBD</i>, named <i>hmmibd-rs</i>, which leverages multi-threaded computing to parallelize IBD inference over genome pairs and which supports optional, user-defined recombination rate maps for more accurate IBD detection and filtration from genomes with non-uniform recombination. We further streamlined large-scale IBD detection by incorporating auxiliary built-in functionalities to preprocess input directly from the standard binary variant call format (BCF) and filter IBD output to reduce disk usage.</p> Results <p>Our new implementation significantly reduces IBD detection computation time nearly linearly with the increased number of CPU threads used; using 128 threads shortens IBD detection time from 5.2&#xa0;days to 1.3&#xa0;h for 220 million pairs of simulated <i>Plasmodium falciparum</i>-like chromosomes, increasing computational speed by approximately 100 × over the single-threaded <i>hmmIBD</i> algorithm. Incorporating non-uniform recombination rates in <i>hmmibd-rs</i> enhances the accuracy of IBD inference by mitigating the overestimation of IBD breakpoints in recombination cold spots and their underestimation in hot spots. Non-uniform rates also improve length filtration of IBD segments, dramatically reducing the rate of false positive in recombination cold spots. When applied to empirical data sets, <i>hmmibd-rs</i> completes the detection of IBD from MalariaGEN Pf7 (n ≈ 10,000 monoclonal samples) within hours, enabling a single-day IBD analysis pipeline for large genomic data sets.</p> Conclusion <p><i>hmmibd-rs</i> builds upon, accelerates, and enhances <i>hmmIBD</i> for efficient and accurate IBD detection, serving as a crucial tool for advancing large-scale malaria genomic surveillance.</p>

错误：搜索内容不能为空，请输入英文关键词
错误：关键词超出字数限制，请精简
高级检索

hmmibd-rs: an enhanced hmmIBD implementation for parallelizable identity-by-descent detection from large-scale Plasmodium genomic data

  • 期刊论文

  • 出版模式：

    开放获取

  • 发表日期：

    2026年2月9日

查看全文
  • Bing Guo,
  • Stephen F. Schaffner,
  • Aimee R. Taylor,
  • Timothy D. O’Connor,
  • Shannon Takala-Harrison

摘要

Background

Identity-by-descent (IBD), which describes recent genetic co-ancestry between pairs of genomes, is a fundamental concept in population genomics. It has been used to estimate genetic relatedness, detect selection signals, and understand population demography. The IBD detection method hmmIBD demonstrates high accuracy in inferring IBD segments between haploid genomes, including Plasmodium falciparum, and is widely used in malaria genomic surveillance. However, the current single-threaded implementation of hmmIBD does not utilize the full capacity of multi-processor computers, making it difficult to apply to large data sets, and does not accommodate non-uniform recombination rates across the genome.

Methods

We developed an enhanced implementation of hmmIBD, named hmmibd-rs, which leverages multi-threaded computing to parallelize IBD inference over genome pairs and which supports optional, user-defined recombination rate maps for more accurate IBD detection and filtration from genomes with non-uniform recombination. We further streamlined large-scale IBD detection by incorporating auxiliary built-in functionalities to preprocess input directly from the standard binary variant call format (BCF) and filter IBD output to reduce disk usage.

Results

Our new implementation significantly reduces IBD detection computation time nearly linearly with the increased number of CPU threads used; using 128 threads shortens IBD detection time from 5.2 days to 1.3 h for 220 million pairs of simulated Plasmodium falciparum-like chromosomes, increasing computational speed by approximately 100 × over the single-threaded hmmIBD algorithm. Incorporating non-uniform recombination rates in hmmibd-rs enhances the accuracy of IBD inference by mitigating the overestimation of IBD breakpoints in recombination cold spots and their underestimation in hot spots. Non-uniform rates also improve length filtration of IBD segments, dramatically reducing the rate of false positive in recombination cold spots. When applied to empirical data sets, hmmibd-rs completes the detection of IBD from MalariaGEN Pf7 (n ≈ 10,000 monoclonal samples) within hours, enabling a single-day IBD analysis pipeline for large genomic data sets.

Conclusion

hmmibd-rs builds upon, accelerates, and enhances hmmIBD for efficient and accurate IBD detection, serving as a crucial tool for advancing large-scale malaria genomic surveillance.