期刊
BMC BIOINFORMATICS
卷 23, 期 1, 页码 -出版社
BMC
DOI: 10.1186/s12859-022-05100-3
关键词
SARS-CoV-2; Wastewater analysis; Variant proportion estimation; Probabilistic modeling; Weighted mixture model
类别
资金
- Slovak Research and Development Agency
- Slovak grant agency [APVV-18-0239]
- Operational Program Integrated Infrastructure [VEGA 1/0463/20, 1/0538/22]
- European Union's Horizon 2020 Research and Innovation Staff Exchange programme under the Marie Sklodowska-Curie grant [ITMS:313011ATL7]
- [872539]
This study proposes a new method for estimating the proportions of SARS-CoV-2 variants in wastewater samples. The method is based on a probabilistic model that captures sequence diversity and sequencing errors, and it has been evaluated for accuracy using different sequencing platforms.
Background: The genomes of SARS-CoV-2 are classified into variants, some of which are monitored as variants of concern (e.g. the Delta variant B.1.617.2 or Omicron variant B.1.1.529). Proportions of these variants circulating in a human population are typically estimated by large-scale sequencing of individual patient samples. Sequencing a mixture of SARS-CoV-2 RNA molecules from wastewater provides a cost-effective alternative, but requires methods for estimating variant proportions in a mixed sample. Results: We propose a new method based on a probabilistic model of sequencing reads, capturing sequence diversity present within individual variants, as well as sequencing errors. The algorithm is implemented in an open source Python program called VirPool. We evaluate the accuracy of VirPool on several simulated and real sequencing data sets from both Illumina and nanopore sequencing platforms, including wastewater samples from Austria and France monitoring the onset of the Alpha variant. Conclusions: VirPool is a versatile tool for wastewater and other mixed-sample analysis that can handle both short- and long-read sequencing data. Our approach does not require pre-selection of characteristic mutations for variant profiles, it is able to use the entire length of reads instead of just the most informative positions, and can also capture haplotype dependencies within a single read.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据