期刊
MOLECULAR BIOLOGY AND EVOLUTION
卷 36, 期 3, 页码 587-603出版社
OXFORD UNIV PRESS
DOI: 10.1093/molbev/msy242
关键词
SNP; transmission clusters; whole-genome sequencing; public health
资金
- Engineering and Physical Sciences Research Council (EPSRC) [EP/K026003/1]
- EPSRC [EP/N014529/1]
- National Institute of General Medical Sciences [U54GM088558]
- EPSRC [EP/N014529/1, EP/K026003/1] Funding Source: UKRI
Whole-genome sequencing (WGS) is increasingly used to aid the understanding of pathogen transmission. A first step in analyzing WGS data is usually to define transmission clusters, sets of cases that are potentially linked by direct transmission. This is often done by including two cases in the same cluster if they are separated by fewer single-nucleotide polymorphisms (SNPs) than a specified threshold. However, there is little agreement as to what an appropriate threshold should be. We propose a probabilistic alternative, suggesting that the key inferential target for transmission clusters is the number of transmissions separating cases. We characterize this by combining the number of SNP differences and the length of time over which those differences have accumulated, using information about case timing, molecular clock, and transmission processes. Our framework has the advantage of allowing for variable mutation rates across the genome and can incorporate other epidemiological data. We use two tuberculosis studies to illustrate the impact of our approach: with British Columbia data by using spatial divisions; with Republic of Moldova data by incorporating antibiotic resistance. Simulation results indicate that our transmission-based method is better in identifying direct transmissions than a SNP threshold, with dissimilarity between clusterings of on average 0.27bits compared with 0.37bits for the SNP-threshold method and 0.84bits for randomly permuted data. These results show that it is likely to outperform the SNP-threshold method where clock rates are variable and sample collection times are spread out. We implement the method in the R package transcluster.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据