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.

Hidden Markov Models on Variable Blocks with a Modal Clustering Algorithm and Applications

Enhanced Bayesian Gaussian hidden Markov mixture clustering for improved knowledge discovery

Clustering hidden Markov models with variational HEM

Bayesian Bi-clustering Methods with Applications in Computational Biology

Sparse Bayesian Hierarchical Modeling of High-dimensional Clustering Problems

Bayesian mixtures of common factor analyzers: Model, variational inference, and applications

Co-Clustering Multi-View Data Using the Latent Block Model

Hidden Markov Model Variants and their Application

Medical Inpatient Journey Modeling and Clustering: A Bayesian Hidden Markov Model Based Approach.

Variational Bayesian Variable Selection for High-Dimensional Hidden Markov Models

BELMM: Bayesian model selection and random walk smoothing in time-series clustering

DBHC: Discrete Bayesian HMM Clustering

A Bayesian hierarchical hidden Markov model for clustering and gene selection: Application to kidney cancer gene expression data

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

Bayesian hidden Markov tree models for clustering genes with shared evolutionary history

A multi-view latent variable model reveals cellular heterogeneity in complex tissues for paired multimodal single-cell data

Clustering Sparse Binary Data with Hierarchical Bayesian Bernoulli Mixture Model

Clustering of single-cell multi-omics data with a multimodal deep learning method

Efficient Bayesian model selection for coupled hidden Markov models with application to infectious diseases

bmVAE: a variational autoencoder method for clustering single-cell mutation data

Topological Hidden Markov Models