Dataset Information

A model-based clustering method to detect infectious disease transmission outbreaks from sequence variation.

ABSTRACT: Clustering infections by genetic similarity is a popular technique for identifying potential outbreaks of infectious disease, in part because sequences are now routinely collected for clinical management of many infections. A diverse number of nonparametric clustering methods have been developed for this purpose. These methods are generally intuitive, rapid to compute, and readily scale with large data sets. However, we have found that nonparametric clustering methods can be biased towards identifying clusters of diagnosis-where individuals are sampled sooner post-infection-rather than the clusters of rapid transmission that are meant to be potential foci for public health efforts. We develop a fundamentally new approach to genetic clustering based on fitting a Markov-modulated Poisson process (MMPP), which represents the evolution of transmission rates along the tree relating different infections. We evaluated this model-based method alongside five nonparametric clustering methods using both simulated and actual HIV sequence data sets. For simulated clusters of rapid transmission, the MMPP clustering method obtained higher mean sensitivity (85%) and specificity (91%) than the nonparametric methods. When we applied these clustering methods to published sequences from a study of HIV-1 genetic clusters in Seattle, USA, we found that the MMPP method categorized about half (46%) as many individuals to clusters compared to the other methods. Furthermore, the mean internal branch lengths that approximate transmission rates were significantly shorter in clusters extracted using MMPP, but not by other methods. We determined that the computing time for the MMPP method scaled linearly with the size of trees, requiring about 30 seconds for a tree of 1,000 tips and about 20 minutes for 50,000 tips on a single computer. This new approach to genetic clustering has significant implications for the application of pathogen sequence analysis to public health, where it is critical to robustly and accurately identify clusters for the most cost-effective deployment of outbreak management and prevention resources.

SUBMITTER: McCloskey RM

PROVIDER: S-EPMC5703573 | biostudies-literature | 2017 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A model-based clustering method to detect infectious disease transmission outbreaks from sequence variation.

McCloskey Rosemary M RM Poon Art F Y AFY

PLoS computational biology 20171113 11

Clustering infections by genetic similarity is a popular technique for identifying potential outbreaks of infectious disease, in part because sequences are now routinely collected for clinical management of many infections. A diverse number of nonparametric clustering methods have been developed for this purpose. These methods are generally intuitive, rapid to compute, and readily scale with large data sets. However, we have found that nonparametric clustering methods can be biased towards ident ...[more]

PMID: 29131825

Similar Datasets

Project description:Whole-genome sequencing of pathogens from host samples becomes more and more routine during infectious disease outbreaks. These data provide information on possible transmission events which can be used for further epidemiologic analyses, such as identification of risk factors for infectivity and transmission. However, the relationship between transmission events and sequence data is obscured by uncertainty arising from four largely unobserved processes: transmission, case observation, within-host pathogen dynamics and mutation. To properly resolve transmission events, these processes need to be taken into account. Recent years have seen much progress in theory and method development, but existing applications make simplifying assumptions that often break up the dependency between the four processes, or are tailored to specific datasets with matching model assumptions and code. To obtain a method with wider applicability, we have developed a novel approach to reconstruct transmission trees with sequence data. Our approach combines elementary models for transmission, case observation, within-host pathogen dynamics, and mutation, under the assumption that the outbreak is over and all cases have been observed. We use Bayesian inference with MCMC for which we have designed novel proposal steps to efficiently traverse the posterior distribution, taking account of all unobserved processes at once. This allows for efficient sampling of transmission trees from the posterior distribution, and robust estimation of consensus transmission trees. We implemented the proposed method in a new R package phybreak. The method performs well in tests of both new and published simulated data. We apply the model to five datasets on densely sampled infectious disease outbreaks, covering a wide range of epidemiological settings. Using only sampling times and sequences as data, our analyses confirmed the original results or improved on them: the more realistic infection times place more confidence in the inferred transmission trees.

Project description:A better characterization of the early growth dynamics of an epidemic is needed to dissect the important drivers of disease transmission, refine existing transmission models, and improve disease forecasts.We introduce a 2-parameter generalized-growth model to characterize the ascending phase of an outbreak and capture epidemic profiles ranging from sub-exponential to exponential growth. We test the model against empirical outbreak data representing a variety of viral pathogens in historic and contemporary populations, and provide simulations highlighting the importance of sub-exponential growth for forecasting purposes.We applied the generalized-growth model to 20 infectious disease outbreaks representing a range of transmission routes. We uncovered epidemic profiles ranging from very slow growth (p=0.14 for the Ebola outbreak in Bomi, Liberia (2014)) to near exponential (p>0.9 for the smallpox outbreak in Khulna (1972), and the 1918 pandemic influenza in San Francisco). The foot-and-mouth disease outbreak in Uruguay displayed a profile of slower growth while the growth pattern of the HIV/AIDS epidemic in Japan was approximately linear. The West African Ebola epidemic provided a unique opportunity to explore how growth profiles vary by geography; analysis of the largest district-level outbreaks revealed substantial growth variations (mean p=0.59, range: 0.14-0.97). The districts of Margibi in Liberia and Bombali and Bo in Sierra Leone had near-exponential growth, while the districts of Bomi in Liberia and Kenema in Sierra Leone displayed near constant incidences.Our findings reveal significant variation in epidemic growth patterns across different infectious disease outbreaks and highlights that sub-exponential growth is a common phenomenon, especially for pathogens that are not airborne. Sub-exponential growth profiles may result from heterogeneity in contact structures or risk groups, reactive behavior changes, or the early onset of interventions strategies, and consideration of "deceleration parameters" may be useful to refine existing mathematical transmission models and improve disease forecasts.

Project description:BACKGROUND: Infectious disease outbreaks in communities can be controlled by early detection and effective prevention measures. Assessing the relative importance of each individual community member with respect to these two processes requires detailed knowledge about the underlying social contact network on which the disease can spread. However, mapping social contact networks is typically too resource-intensive to be a practical possibility for most communities and institutions. METHODS: Here, we describe a simple, low-cost method - called collocation ranking - to assess individual importance for early detection and targeted intervention strategies that are easily implementable in practice. The method is based on knowledge about individual collocation which is readily available in many community settings such as schools, offices, hospitals, and so on. We computationally validate our method in a school setting by comparing the outcome of the method against computational predictions based on outbreak simulations on an empirical high-resolution contact network. We compare collocation ranking to other methods for assessing the epidemiological importance of the members of a population. To this end, we select subpopulations of the school population by applying these assessment methods to the population and adding individuals to the subpopulation according to their individual rank. Then, we assess how suited these subpopulations are for early detection and targeted intervention strategies. RESULTS: Likelihood and timing of infection during an outbreak are important features for early detection and targeted intervention strategies. Subpopulations selected by the collocation ranking method show a substantially higher average infection probability and an earlier onset of symptoms than randomly selected subpopulations. Furthermore, these subpopulations selected by the collocation ranking method were close to the optimum. CONCLUSIONS: Our results indicate that collocation ranking is a highly effective method to assess individual importance, providing critical low-cost information for the development of sentinel surveillance systems and prevention strategies.

Dataset Information

A model-based clustering method to detect infectious disease transmission outbreaks from sequence variation.

Publications

A model-based clustering method to detect infectious disease transmission outbreaks from sequence variation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets