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.

Review 1: "estimating the Reproduction Number and Transmission Heterogeneity from the Size Distribution of Clusters of Identical Pathogen Sequences"

Estimating the reproduction number and transmission heterogeneity from the size distribution of clusters of identical pathogen sequences

Deterministic epidemic models overestimate the basic reproduction number of observed outbreaks

Monitoring Real-Time Transmission Heterogeneity from Incidence Data

Estimation of the Basic Reproductive Number and Mean Serial Interval of a Novel Pathogen in a Small, Well-Observed Discrete Population

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

Characterizing the reproduction number of epidemics with early sub-exponential growth dynamics

Distributed Reproduction Numbers of Networked Epidemics

Exploring the nuances of R0: Eight estimates and application to 2009 pandemic influenza

Evaluating undercounts in epidemics: response to Maruotti et al. 2022

Estimating the overdispersion in COVID-19 transmission using outbreak sizes outside China

Improved estimation of the effective reproduction number with heterogeneous transmission rates and reporting delays

Unveiling the reproduction number scaling in characterizing social contagion coverage

Estimating the basic reproduction number at the beginning of an outbreak

Estimating R and overdispersion in secondary cases from the size of identical sequence clusters of SARS-CoV-2

Using an Epidemiological Model for Phylogenetic Inference Reveals Density Dependence in HIV Transmission

Impacts and shortcomings of genetic clustering methods for infectious disease outbreaks

Type reproduction number for epidemic models on heterogeneous networks

A Comparative Analysis of Statistical Methods to Estimate the Reproduction Number in Emerging Epidemics, With Implications for the Current Coronavirus Disease 2019 (COVID-19) Pandemic

Refining Estimation of the Instantaneous Reproduction Number During Early Pandemic Stages: Addressing Case-Reporting Variability and Serial Interval Uncertainty