Multi-species SIR models from a dynamical Bayesian perspective

Estimation of time-varying reproduction numbers underlying epidemiological processes: A new statistical tool for the COVID-19 pandemic

The coronavirus pandemic has rapidly evolved into an unprecedented crisis. The susceptible-infectious-removed (SIR) model and its variants have been used for modeling the pandemic. However, time-independent parameters in the classical models may not capture the dynamic transmission and removal processes, governed by virus containment strategies taken at various phases of the epidemic. Moreover, few models account for possible inaccuracies of the reported cases. We propose a Poisson model with time-dependent transmission and removal rates to account for possible random errors in reporting and estimate a time-dependent disease reproduction number, which may reflect the effectiveness of virus control strategies. We apply our method to study the pandemic in several severely impacted countries, and analyze and forecast the evolving spread of the coronavirus. We have developed an interactive web application to facilitate readers' use of our method.

Individual Specialization and Multihost Epidemics: Disease Spread in Plant-Pollinator Networks

Many parasites infect multiple species and persist through a combination of within- and between-species transmission. Multispecies transmission networks are typically constructed at the species level, linking two species if any individuals of those species interact. However, generalist species often consist of specialized individuals that prefer different subsets of available resources, so individual- and species-level contact networks can differ systematically. To explore the epidemiological impacts of host specialization, we build and study a model for pollinator pathogens on plant-pollinator networks, in which individual pollinators have dynamic preferences for different flower species. We find that modeling and analysis that ignore individual host specialization can predict die-off of a disease that is actually strongly persistent and can badly over- or underpredict steady-state disease prevalence. Effects of individual preferences remain substantial whenever mean preference duration exceeds half of the mean time from infection to recovery or death. Similar results hold in a model where hosts foraging in different habitats have different frequencies of contact with an environmental reservoir for the pathogen. Thus, even if all hosts have the same long-run average behavior, dynamic individual differences can profoundly affect disease persistence and prevalence.

Hand, Foot, and Mouth Disease in China: Modeling Epidemic Dynamics of Enterovirus Serotypes and Implications for Vaccination

BACKGROUND

Hand, foot, and mouth disease (HFMD) is a common childhood illness caused by serotypes of the Enterovirus A species in the genus Enterovirus of the Picornaviridae family. The disease has had a substantial burden throughout East and Southeast Asia over the past 15 y. China reported 9 million cases of HFMD between 2008 and 2013, with the two serotypes Enterovirus A71 (EV-A71) and Coxsackievirus A16 (CV-A16) being responsible for the majority of these cases. Three recent phase 3 clinical trials showed that inactivated monovalent EV-A71 vaccines manufactured in China were highly efficacious against HFMD associated with EV-A71, but offered no protection against HFMD caused by CV-A16. To better inform vaccination policy, we used mathematical models to evaluate the effect of prospective vaccination against EV-A71-associated HFMD and the potential risk of serotype replacement by CV-A16. We also extended the model to address the co-circulation, and implications for vaccination, of additional non-EV-A71, non-CV-A16 serotypes of enterovirus.

METHODS AND FINDINGS

Weekly reports of HFMD incidence from 31 provinces in Mainland China from 1 January 2009 to 31 December 2013 were used to fit multi-serotype time series susceptible-infected-recovered (TSIR) epidemic models. We obtained good model fit for the two-serotype TSIR with cross-protection, capturing the seasonality and geographic heterogeneity of province-level transmission, with strong correlation between the observed and simulated epidemic series. The national estimate of the basic reproduction number, R0, weighted by provincial population size, was 26.63 for EV-A71 (interquartile range [IQR]: 23.14, 30.40) and 27.13 for CV-A16 (IQR: 23.15, 31.34), with considerable variation between provinces (however, predictions about the overall impact of vaccination were robust to this variation). EV-A71 incidence was projected to decrease monotonically with higher coverage rates of EV-A71 vaccination. Across provinces, CV-A16 incidence in the post-EV-A71-vaccination period remained either comparable to or only slightly increased from levels prior to vaccination. The duration and strength of cross-protection following infection with EV-A71 or CV-A16 was estimated to be 9.95 wk (95% confidence interval [CI]: 3.31, 23.40) in 68% of the population (95% CI: 37%, 96%). Our predictions are limited by the necessarily short and under-sampled time series and the possible circulation of unidentified serotypes, but, nonetheless, sensitivity analyses indicate that our results are robust in predicting that the vaccine should drastically reduce incidence of EV-A71 without a substantial competitive release of CV-A16.

CONCLUSIONS

The ability of our models to capture the observed epidemic cycles suggests that herd immunity is driving the epidemic dynamics caused by the multiple serotypes of enterovirus. Our results predict that the EV-A71 and CV-A16 serotypes provide a temporary immunizing effect against each other. Achieving high coverage rates of EV-A71 vaccination would be necessary to eliminate the ongoing transmission of EV-A71, but serotype replacement by CV-A16 following EV-A71 vaccination is likely to be transient and minor compared to the corresponding reduction in the burden of EV-A71-associated HFMD. Therefore, a mass EV-A71 vaccination program of infants and young children should provide significant benefits in terms of a reduction in overall HFMD burden.

Whole Genome Sequencing demonstrates that Geographic Variation of Escherichia coli O157 Genotypes Dominates Host Association

Genetic variation in an infectious disease pathogen can be driven by ecological niche dissimilarities arising from different host species and different geographical locations. Whole genome sequencing was used to compare E. coli O157 isolates from host reservoirs (cattle and sheep) from Scotland and to compare genetic variation of isolates (human, animal, environmental/food) obtained from Scotland, New Zealand, Netherlands, Canada and the USA. Nei’s genetic distance calculated from core genome single nucleotide polymorphisms (SNPs) demonstrated that the animal isolates were from the same population. Investigation of the Shiga toxin bacteriophage and their insertion sites (SBI typing) revealed that cattle and sheep isolates had statistically indistinguishable rarefaction profiles, diversity and genotypes. In contrast, isolates from different countries exhibited significant differences in Nei’s genetic distance and SBI typing. Hence, after successful international transmission, which has occurred on multiple occasions, local genetic variation occurs, resulting in a global patchwork of continental and trans-continental phylogeographic clades. These findings are important for three reasons: first, understanding transmission and evolution of infectious diseases associated with multiple host reservoirs and multi-geographic locations; second, highlighting the relevance of the sheep reservoir when considering farm based interventions; and third, improving our understanding of why human disease incidence varies across the world.