Comparing the transmission potential from sequence and surveillance data of 2009 North American influenza pandemic waves

Venkata R Duvvuri; Joseph T Hicks; Lambodhar Damodaran; Martin Grunnill; Thomas Braukmann; Jianhong Wu; Jonathan B Gubbay; Samir N Patel; Justin Bahl

doi:10.1016/j.idm.2023.02.003

Comparing the transmission potential from sequence and surveillance data of 2009 North American influenza pandemic waves

Infect Dis Model. 2023 Feb 16;8(1):240-252. doi: 10.1016/j.idm.2023.02.003. eCollection 2023 Mar.

Authors

Venkata R Duvvuri^{1

2

3

4

5}, Joseph T Hicks^{4

5}, Lambodhar Damodaran^{4

5}, Martin Grunnill^{1

3}, Thomas Braukmann¹, Jianhong Wu³, Jonathan B Gubbay^{1

2}, Samir N Patel^{1

2}, Justin Bahl^{4

5

6}

Affiliations

¹ Public Health Ontario, Toronto, Ontario, Canada.
² Department of Laboratory Medicine and Pathobiology, Faculty of Medicine, University of Toronto, Toronto, Ontario, Canada.
³ Laboratory for Industrial and Applied Mathematics, Department of Mathematics and Statistics, York University, Toronto, Ontario, Canada.
⁴ Center for the Ecology of Infectious Disease, Department of Infectious Diseases, Institute of Bioinformatics, University of Georgia, Athens, Georgia.
⁵ Department of Epidemiology and Biostatistics, Institute of Bioinformatics, University of Georgia, Athens, Georgia.
⁶ Duke-NUS Graduate Medical School, Singapore.

Abstract

Technological advancements in phylodynamic modeling coupled with the accessibility of real-time pathogen genetic data are increasingly important for understanding the infectious disease transmission dynamics. In this study, we compare the transmission potentials of North American influenza A(H1N1)pdm09 derived from sequence data to that derived from surveillance data. The impact of the choice of tree-priors, informative epidemiological priors, and evolutionary parameters on the transmission potential estimation is evaluated. North American Influenza A(H1N1)pdm09 hemagglutinin (HA) gene sequences are analyzed using the coalescent and birth-death tree prior models to estimate the basic reproduction number (R ₀ ). Epidemiological priors gathered from published literature are used to simulate the birth-death skyline models. Path-sampling marginal likelihood estimation is conducted to assess model fit. A bibliographic search to gather surveillance-based R ₀ values were consistently lower (mean ≤ 1.2) when estimated by coalescent models than by the birth-death models with informative priors on the duration of infectiousness (mean ≥ 1.3 to ≤2.88 days). The user-defined informative priors for use in the birth-death model shift the directionality of epidemiological and evolutionary parameters compared to non-informative estimates. While there was no certain impact of clock rate and tree height on the R ₀ estimation, an opposite relationship was observed between coalescent and birth-death tree priors. There was no significant difference (p = 0.46) between the birth-death model and surveillance R ₀ estimates. This study concludes that tree-prior methodological differences may have a substantial impact on the transmission potential estimation as well as the evolutionary parameters. The study also reports a consensus between the sequence-based R ₀ estimation and surveillance-based R ₀ estimates. Altogether, these outcomes shed light on the potential role of phylodynamic modeling to augment existing surveillance and epidemiological activities to better assess and respond to emerging infectious diseases.

Keywords: Birth-death models; Coalescent growth models; Pandemic 2009 H1N1; Pathogen sequence data; Phylodynamics; Public health; Reproduction number.