Estimating R and overdispersion in secondary cases from the size of identical sequence clusters of SARS-CoV-2
Emma B Hodcroft,Martin S Wohlfender,Richard A Neher,Julien Riou,Christian L Althaus
DOI: https://doi.org/10.1101/2024.05.26.24307940
2024-09-28
Abstract:The wealth of genomic data that was generated during the COVID-19 pandemic provides an exceptional opportunity to obtain information on the transmission of SARS-CoV-2. Specifically, there is great interest to better understand how the effective reproduction number R and the overdispersion of secondary cases, which can be quantified by the negative binomial dispersion parameter , changed over time and across regions and viral variants. The aim of our study was to develop a Bayesian framework to infer R and from viral sequence data. First, we developed a mathematical model for the distribution of the size of identical sequence clusters, in which we integrated viral transmission, the mutation rate of the virus, and incomplete case-detection. Second, we implemented this model within a Bayesian inference framework, allowing the estimation of R and from genomic data only. We validated this model in a simulation study. Third, we identified clusters of identical sequences in all SARS-CoV-2 sequences in 2021 from Switzerland, Denmark, and Germany that were available on GISAID. We obtained monthly estimates of the posterior distribution of R and , with the resulting R estimates slightly lower than resulting obtained by other methods, and comparable with previous results. We found comparatively higher estimates of in Denmark which suggests less opportunities for superspreading and more controlled transmission compared to the other countries in 2021. Our model included an estimation of the case detection and sampling probability, but the estimates obtained had large uncertainty, reflecting the difficulty of estimating these parameters simultaneously. Our study presents a novel method to infer information on the transmission of infectious diseases and its heterogeneity using genomic data. With increasing availability of sequences of pathogens in the future, we expect that our method has the potential to provide new insights into the transmission and the overdispersion in secondary cases of other pathogens.