76
|
Mina MJ, Grenfell BT, Metcalf CJE. Response to Comment on "Long-term measles-induced immunomodulation increases overall childhood infectious disease mortality". SCIENCE (NEW YORK, N.Y.) 2019; 365:365/6449/eaax6498. [PMID: 31296742 DOI: 10.1126/science.aax6498] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 06/11/2019] [Indexed: 11/02/2022]
Abstract
Thakkar and McCarthy suggest that periodicity in measles incidence artifactually drives our estimates of a 2- to 3-year duration of measles "immune-amnesia." We show that periodicity has a negligible effect relative to the immunological signal we detect, and demonstrate that immune-amnesia is largely undetectable in small populations with large fluctuations in mortality of the type they use for illustration.
Collapse
|
77
|
Lau MSY, Grenfell BT, Worby CJ, Gibson GJ. Model diagnostics and refinement for phylodynamic models. PLoS Comput Biol 2019; 15:e1006955. [PMID: 30951528 PMCID: PMC6469796 DOI: 10.1371/journal.pcbi.1006955] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2018] [Revised: 04/17/2019] [Accepted: 03/18/2019] [Indexed: 11/29/2022] Open
Abstract
Phylodynamic modelling, which studies the joint dynamics of epidemiological and evolutionary processes, has made significant progress in recent years due to increasingly available genomic data and advances in statistical modelling. These advances have greatly improved our understanding of transmission dynamics of many important pathogens. Nevertheless, there remains a lack of effective, targetted diagnostic tools for systematically detecting model mis-specification. Development of such tools is essential for model criticism, refinement, and calibration. The idea of utilising latent residuals for model assessment has already been exploited in general spatio-temporal epidemiological settings. Specifically, by proposing appropriately designed non-centered, re-parameterizations of a given epidemiological process, one can construct latent residuals with known sampling distributions which can be used to quantify evidence of model mis-specification. In this paper, we extend this idea to formulate a novel model-diagnostic framework for phylodynamic models. Using simulated examples, we show that our framework may effectively detect a particular form of mis-specification in a phylodynamic model, particularly in the event of superspreading. We also exemplify our approach by applying the framework to a dataset describing a local foot-and-mouth (FMD) outbreak in the UK, eliciting strong evidence against the assumption of no within-host-diversity in the outbreak. We further demonstrate that our framework can facilitate model calibration in real-life scenarios, by proposing a within-host-diversity model which appears to offer a better fit to data than one that assumes no within-host-diversity of FMD virus. Integrated modelling of conventional epidemiological data and modern genomic data (i.e. phylodynamics) has made significant progress in recent years, due to the ever-increasing availability of genomic data and development of statistical methods. However, there is a lack of tools for carrying out effective diagnostics for phylodynamic models. We propose a novel model diagnostic framework that involves a latent residual process which is a priori independent of model assumptions and which can be used to quantify, and reveal the nature of, model inadequacy. Our results suggest that our framework may systematically detect deviation from a particular model assumption and greatly facilitate subsequent model calibration.
Collapse
|
78
|
Abstract
Healthcare-associated infections (HAIs) pose a significant burden to patient safety. Institutions can implement hospital infection control (HIC) measures to reduce the impact of HAIs. Since patients can carry pathogens between institutions, there is an economic incentive for hospitals to free ride on the HIC investments of other facilities. Subsidies for infection control by public health authorities could encourage regional spending on HIC. We develop coupled mathematical models of epidemiology and hospital behavior in a game-theoretic framework to investigate how hospitals may change spending behavior in response to subsidies. We demonstrate that under a limited budget, a dollar-for-dollar matching grant outperforms both a fixed-amount subsidy and a subsidy on uninfected patients in reducing the number of HAIs in a single institution. Additionally, when multiple hospitals serve a community, funding priority should go to the hospital with a lower transmission rate. Overall, subsidies incentivize HIC spending and reduce the overall prevalence of HAIs.
Collapse
|
79
|
Dalziel BD, Kissler S, Gog JR, Viboud C, Bjørnstad ON, Metcalf CJE, Grenfell BT. Urbanization and humidity shape the intensity of influenza epidemics in U.S. cities. Science 2019; 362:75-79. [PMID: 30287659 PMCID: PMC6510303 DOI: 10.1126/science.aat6030] [Citation(s) in RCA: 175] [Impact Index Per Article: 35.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 08/10/2018] [Indexed: 01/14/2023]
Abstract
Influenza epidemics vary in intensity from year to year, driven by climatic conditions and by viral antigenic evolution. However, important spatial variation remains unexplained. Here we show predictable differences in influenza incidence among cities, driven by population size and structure. Weekly incidence data from 603 cities in the United States reveal that epidemics in smaller cities are focused on shorter periods of the influenza season, whereas in larger cities, incidence is more diffuse. Base transmission potential estimated from city-level incidence data is positively correlated with population size and with spatiotemporal organization in population density, indicating a milder response to climate forcing in metropolises. This suggests that urban centers incubate critical chains of transmission outside of peak climatic conditions, altering the spatiotemporal geometry of herd immunity.
Collapse
|
80
|
Morris SE, Yates AJ, de Swart RL, de Vries RD, Mina MJ, Nelson AN, Lin WHW, Kouyos RD, Griffin DE, Grenfell BT. Modeling the measles paradox reveals the importance of cellular immunity in regulating viral clearance. PLoS Pathog 2018; 14:e1007493. [PMID: 30592772 PMCID: PMC6310241 DOI: 10.1371/journal.ppat.1007493] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2018] [Accepted: 11/29/2018] [Indexed: 12/15/2022] Open
Abstract
Measles virus (MV) is a highly contagious member of the Morbillivirus genus that remains a major cause of childhood mortality worldwide. Although infection induces a strong MV-specific immune response that clears viral load and confers lifelong immunity, transient immunosuppression can also occur, leaving the host vulnerable to colonization from secondary pathogens. This apparent contradiction of viral clearance in the face of immunosuppression underlies what is often referred to as the 'measles paradox', and remains poorly understood. To explore the mechanistic basis underlying the measles paradox, and identify key factors driving viral clearance, we return to a previously published dataset of MV infection in rhesus macaques. These data include virological and immunological information that enable us to fit a mathematical model describing how the virus interacts with the host immune system. In particular, our model incorporates target cell depletion through infection of host immune cells-a hallmark of MV pathology that has been neglected from previous models. We find the model captures the data well, and that both target cell depletion and immune activation are required to explain the overall dynamics. Furthermore, by simulating conditions of increased target cell availability and suppressed cellular immunity, we show that the latter causes greater increases in viral load and delays to MV clearance. Overall, this signals a more dominant role for cellular immunity in resolving acute MV infection. Interestingly, we find contrasting dynamics dominated by target cell depletion when viral fitness is increased. This may have wider implications for animal morbilliviruses, such as canine distemper virus (CDV), that cause fatal target cell depletion in their natural hosts. To our knowledge this work represents the first fully calibrated within-host model of MV dynamics and, more broadly, provides a new platform from which to explore the complex mechanisms underlying Morbillivirus infection.
Collapse
|
81
|
Gadroen K, Dodd CN, Masclee GMC, de Ridder MAJ, Weibel D, Mina MJ, Grenfell BT, Sturkenboom MCJM, van de Vijver DAMC, de Swart RL. Impact and longevity of measles-associated immune suppression: a matched cohort study using data from the THIN general practice database in the UK. BMJ Open 2018; 8:e021465. [PMID: 30413497 PMCID: PMC6231568 DOI: 10.1136/bmjopen-2017-021465] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
OBJECTIVE To test the hypothesis that measles infection increases the incidence of non-measles infectious diseases over a prolonged period of time. DESIGN A population-based matched cohort study. DATA SOURCES This study examined children aged 1-15 years in The Health Improvement Network UK general practice medical records database. Participants included 2228 patients diagnosed with measles between 1990 and 2014, which were matched on age, sex, general practitioner practice and calendar year with 19 930 children without measles. All controls had received at least one measles vaccination. Children with a history of immune-compromising conditions or with immune-suppressive treatment were excluded. PRIMARY OUTCOME MEASURES Incidence rate ratio (IRR) of infections, anti-infective prescriptions and all-cause hospitalisations following measles in predetermined periods using multivariate analysis to adjust for confounding variables. RESULTS In children with measles, the incidence rate for non-measles infectious disease was significantly increased in each time period assessed up to 5 years postmeasles: 43% in the first month (IRR: 1.43; 95% CI 1.22 to 1.68), 22% from month one to the first year (IRR: 1.22; 95% CI 1.14 to 1.31), 10% from year 1 to 2.5 years (IRR: 1.10; 95% CI 1.02 to 1.19) and 15% (IRR: 1.15; 95% CI 1.06 to 1.25) in years 2.5 to 5 years of follow-up. Children with measles were more than three times as likely to receive an anti-infective prescription in the first month and 15%-24% more likely between the first month and 5 years. The rate of hospitalisation in children with measles was increased only in the month following diagnosis but not thereafter (IRR: 2.83; 95% CI 1.72 to 4.67). CONCLUSION Following measles, children had increased rates of diagnosed infections, requiring increased prescribing of antimicrobial therapies. This population-based matched cohort study supports the hypothesis that measles has a prolonged impact on host resistance to non-measles infectious diseases.
Collapse
|
82
|
Kissler SM, Gog JR, Viboud C, Charu V, Bjørnstad ON, Simonsen L, Grenfell BT. Geographic transmission hubs of the 2009 influenza pandemic in the United States. Epidemics 2018; 26:86-94. [PMID: 30327253 DOI: 10.1016/j.epidem.2018.10.002] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2017] [Revised: 10/05/2018] [Accepted: 10/08/2018] [Indexed: 10/28/2022] Open
Abstract
A key issue in infectious disease epidemiology is to identify and predict geographic sites of epidemic establishment that contribute to onward spread, especially in the context of invasion waves of emerging pathogens. Conventional wisdom suggests that these sites are likely to be in densely-populated, well-connected areas. For pandemic influenza, however, epidemiological data have not been available at a fine enough geographic resolution to test this assumption. Here, we make use of fine-scale influenza-like illness incidence data derived from electronic medical claims records gathered from 834 3-digit ZIP (postal) codes across the US to identify the key geographic establishment sites, or "hubs", of the autumn wave of the 2009 A/H1N1pdm influenza pandemic in the United States. A mechanistic spatial transmission model is fit to epidemic onset times inferred from the data. Hubs are identified by tracing the most probable transmission routes back to a likely first establishment site. Four hubs are identified: two in the southeastern US, one in the central valley of California, and one in the midwestern US. According to the model, 75% of the 834 observed ZIP-level outbreaks in the US were seeded by these four hubs or their epidemiological descendants. Counter-intuitively, the pandemic hubs do not coincide with large and well-connected cities, indicating that factors beyond population density and travel volume are necessary to explain the establishment sites of the major autumn wave of the pandemic. Geographic regions are identified where infection can be statistically traced back to a hub, providing a testable prediction of the outbreak's phylogeography. Our method therefore provides an important way forward to reconcile spatial diffusion patterns inferred from epidemiological surveillance data and pathogen sequence data.
Collapse
|
83
|
Takahashi S, Metcalf CJE, Arima Y, Fujimoto T, Shimizu H, Rogier van Doorn H, Le Van T, Chan YF, Farrar JJ, Oishi K, Grenfell BT. Epidemic dynamics, interactions and predictability of enteroviruses associated with hand, foot and mouth disease in Japan. J R Soc Interface 2018; 15:rsif.2018.0507. [PMID: 30209044 PMCID: PMC6170776 DOI: 10.1098/rsif.2018.0507] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Accepted: 08/20/2018] [Indexed: 12/28/2022] Open
Abstract
Outbreaks of hand, foot and mouth disease have been documented in Japan since 1963. This disease is primarily caused by the two closely related serotypes of Enterovirus A71 (EV-A71) and Coxsackievirus A16 (CV-A16). Here, we analyse Japanese virologic and syndromic surveillance time-series data from 1982 to 2015. As in some other countries in the Asia Pacific region, EV-A71 in Japan has a 3 year cyclical component, whereas CV-A16 is predominantly annual. We observe empirical signatures of an inhibitory interaction between the serotypes; virologic lines of evidence suggest they may indeed interact immunologically. We fit the time series to mechanistic epidemiological models: as a first-order effect, we find the data consistent with single-serotype susceptible–infected–recovered dynamics. We then extend the modelling to incorporate an inhibitory interaction between serotypes. Our results suggest the existence of a transient cross-protection and possible asymmetry in its strength such that CV-A16 serves as a stronger forcing on EV-A71. Allowing for asymmetry yields accurate out-of-sample predictions and the directionality of this effect is consistent with the virologic literature. Confirmation of these hypothesized interactions would have important implications for understanding enterovirus epidemiology and informing vaccine development. Our results highlight the general implication that even subtle interactions could have qualitative impacts on epidemic dynamics and predictability.
Collapse
|
84
|
Lau MSY, Grenfell BT. Vaccination under uncertainty. Nat Ecol Evol 2018; 2:1350-1351. [DOI: 10.1038/s41559-018-0652-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
85
|
Goldlust SM, Thuan PD, Giang DDH, Thang ND, Thwaites GE, Farrar J, Thanh NV, Nguyen TD, Grenfell BT, Boni MF, Hien TT. The decline of malaria in Vietnam, 1991-2014. Malar J 2018; 17:226. [PMID: 29880051 PMCID: PMC5992833 DOI: 10.1186/s12936-018-2372-8] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Accepted: 05/28/2018] [Indexed: 02/04/2023] Open
Abstract
Background Despite the well-documented clinical efficacy of artemisinin-based combination therapy (ACT) against malaria, the population-level effects of ACT have not been studied thoroughly until recently. An ideal case study for these population-level effects can be found in Vietnam’s gradual adoption of artemisinin in the 1990s. Methods and results Analysis of Vietnam’s national annual malaria reports (1991–2014) revealed that a 10% increase in artemisinin procurement corresponded to a 32.8% (95% CI 27.7–37.5%) decline in estimated malaria cases. There was no consistent national or regional effect of vector control on malaria. The association between urbanization and malaria was generally negative and sometimes statistically significant. Conclusions The decline of malaria in Vietnam can largely be attributed to the adoption of artemisinin-based case management. Recent analyses from Africa showed that insecticide-treated nets had the greatest effect on lowering malaria prevalence, suggesting that the success of interventions is region-specific. Continuing malaria elimination efforts should focus on both vector control and increased access to ACT. Electronic supplementary material The online version of this article (10.1186/s12936-018-2372-8) contains supplementary material, which is available to authorized users.
Collapse
|
86
|
Morris SE, Freiesleben de Blasio B, Viboud C, Wesolowski A, Bjørnstad ON, Grenfell BT. Analysis of multi-level spatial data reveals strong synchrony in seasonal influenza epidemics across Norway, Sweden, and Denmark. PLoS One 2018; 13:e0197519. [PMID: 29771952 PMCID: PMC5957349 DOI: 10.1371/journal.pone.0197519] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Accepted: 05/03/2018] [Indexed: 12/02/2022] Open
Abstract
Population structure, spatial diffusion, and climatic conditions mediate the spatiotemporal spread of seasonal influenza in temperate regions. However, much of our knowledge of these dynamics stems from a few well-studied countries, such as the United States (US), and the extent to which this applies in different demographic and climatic environments is not fully understood. Using novel data from Norway, Sweden, and Denmark, we applied wavelet analysis and non-parametric spatial statistics to explore the spatiotemporal dynamics of influenza transmission at regional and international scales. We found the timing and amplitude of epidemics were highly synchronized both within and between countries, despite the geographical isolation of many areas in our study. Within Norway, this synchrony was most strongly modulated by population size, confirming previous findings that hierarchical spread between larger populations underlies seasonal influenza dynamics at regional levels. However, we found no such association when comparing across countries, suggesting that other factors become important at the international scale. Finally, to frame our results within a wider global context, we compared our findings from Norway to those from the US. After correcting for differences in geographic scale, we unexpectedly found higher levels of synchrony in Norway, despite its smaller population size. We hypothesize that this greater synchrony may be driven by more favorable and spatially uniform climatic conditions, although there are other likely factors we were unable to consider (such as reduced variation in school term times and differences in population movements). Overall, our results highlight the importance of comparing influenza spread at different spatial scales and across diverse geographic regions in order to better understand the complex mechanisms underlying disease dynamics.
Collapse
|
87
|
Pons-Salort M, Oberste MS, Pallansch MA, Abedi GR, Takahashi S, Grenfell BT, Grassly NC. The seasonality of nonpolio enteroviruses in the United States: Patterns and drivers. Proc Natl Acad Sci U S A 2018; 115:3078-3083. [PMID: 29507246 PMCID: PMC5866597 DOI: 10.1073/pnas.1721159115] [Citation(s) in RCA: 59] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Nonpolio enteroviruses are diverse and common viruses that can circulate year-round but tend to peak in summer. Although most infections are asymptomatic, they can result in a wide range of neurological and other diseases. Many serotypes circulate every year, and different serotypes predominate in different years, but the drivers of their geographical and temporal dynamics are not understood. We use national enterovirus surveillance data collected by the US Centers for Disease Control and Prevention during 1983-2013, as well as demographic and climatic data for the same period, to study the patterns and drivers of the seasonality of these infections. We find that the seasonal pattern of enterovirus cases is spatially structured in the United States and similar to that observed for historical prevaccination poliomyelitis (1931-1954). We identify latitudinal gradients for the amplitude and the timing of the peak of cases, meaning that those are more regularly distributed all year-round in the south and have a more pronounced peak that arrives later toward the north. The peak is estimated to occur between July and September across the United States, and 1 month earlier than that for historical poliomyelitis. Using mixed-effects models, we find that climate, but not demography, is likely to drive the seasonal pattern of enterovirus cases and that the dew point temperature alone explains ∼30% of the variation in the intensity of transmission. Our study contributes to a better understanding of the epidemiology of enteroviruses, demonstrates important similarities in their circulation dynamics with polioviruses, and identifies potential drivers of their seasonality.
Collapse
|
88
|
Lee EC, Arab A, Goldlust SM, Viboud C, Grenfell BT, Bansal S. Deploying digital health data to optimize influenza surveillance at national and local scales. PLoS Comput Biol 2018. [PMID: 29513661 PMCID: PMC5858836 DOI: 10.1371/journal.pcbi.1006020] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
The surveillance of influenza activity is critical to early detection of epidemics and pandemics and the design of disease control strategies. Case reporting through a voluntary network of sentinel physicians is a commonly used method of passive surveillance for monitoring rates of influenza-like illness (ILI) worldwide. Despite its ubiquity, little attention has been given to the processes underlying the observation, collection, and spatial aggregation of sentinel surveillance data, and its subsequent effects on epidemiological understanding. We harnessed the high specificity of diagnosis codes in medical claims from a database that represented 2.5 billion visits from upwards of 120,000 United States healthcare providers each year. Among influenza seasons from 2002-2009 and the 2009 pandemic, we simulated limitations of sentinel surveillance systems such as low coverage and coarse spatial resolution, and performed Bayesian inference to probe the robustness of ecological inference and spatial prediction of disease burden. Our models suggest that a number of socio-environmental factors, in addition to local population interactions, state-specific health policies, as well as sampling effort may be responsible for the spatial patterns in U.S. sentinel ILI surveillance. In addition, we find that biases related to spatial aggregation were accentuated among areas with more heterogeneous disease risk, and sentinel systems designed with fixed reporting locations across seasons provided robust inference and prediction. With the growing availability of health-associated big data worldwide, our results suggest mechanisms for optimizing digital data streams to complement traditional surveillance in developed settings and enhance surveillance opportunities in developing countries. Influenza contributes substantially to global morbidity and mortality each year, and epidemiological surveillance for influenza is typically conducted by sentinel physicians and health care providers recruited to report cases of influenza-like illness. While population coverage and representativeness, and geographic distribution are considered during sentinel provider recruitment, systems cannot always achieve these standards due to the administrative burdens of data collection. We present spatial estimates of influenza disease burden across United States counties by leveraging the volume and fine spatial resolution of medical claims data, and existing socio-environmental hypotheses about the determinants of influenza disease disease burden. Using medical claims as a testbed, this study adds to literature on the optimization of surveillance system design by considering conditions of limited reporting and spatial aggregation. We highlight the importance of considering sampling biases and reporting locations when interpreting surveillance data, and suggest that local mobility and regional policies may be critical to understanding the spatial distribution of reported influenza-like illness.
Collapse
|
89
|
Morris DH, Gostic KM, Pompei S, Bedford T, Łuksza M, Neher RA, Grenfell BT, Lässig M, McCauley JW. Predictive Modeling of Influenza Shows the Promise of Applied Evolutionary Biology. Trends Microbiol 2018; 26:102-118. [PMID: 29097090 PMCID: PMC5830126 DOI: 10.1016/j.tim.2017.09.004] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Revised: 09/06/2017] [Accepted: 09/19/2017] [Indexed: 01/16/2023]
Abstract
Seasonal influenza is controlled through vaccination campaigns. Evolution of influenza virus antigens means that vaccines must be updated to match novel strains, and vaccine effectiveness depends on the ability of scientists to predict nearly a year in advance which influenza variants will dominate in upcoming seasons. In this review, we highlight a promising new surveillance tool: predictive models. Based on data-sharing and close collaboration between the World Health Organization and academic scientists, these models use surveillance data to make quantitative predictions regarding influenza evolution. Predictive models demonstrate the potential of applied evolutionary biology to improve public health and disease control. We review the state of influenza predictive modeling and discuss next steps and recommendations to ensure that these models deliver upon their considerable biomedical promise.
Collapse
|
90
|
Dalziel BD, Lau MSY, Tiffany A, McClelland A, Zelner J, Bliss JR, Grenfell BT. Unreported cases in the 2014-2016 Ebola epidemic: Spatiotemporal variation, and implications for estimating transmission. PLoS Negl Trop Dis 2018; 12:e0006161. [PMID: 29357363 PMCID: PMC5806896 DOI: 10.1371/journal.pntd.0006161] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2017] [Revised: 02/09/2018] [Accepted: 12/11/2017] [Indexed: 11/18/2022] Open
Abstract
In the recent 2014–2016 Ebola epidemic in West Africa, non-hospitalized cases were an important component of the chain of transmission. However, non-hospitalized cases are at increased risk of going unreported because of barriers to access to healthcare. Furthermore, underreporting rates may fluctuate over space and time, biasing estimates of disease transmission rates, which are important for understanding spread and planning control measures. We performed a retrospective analysis on community deaths during the recent Ebola epidemic in Sierra Leone to estimate the number of unreported non-hospitalized cases, and to quantify how Ebola reporting rates varied across locations and over time. We then tested if variation in reporting rates affected the estimates of disease transmission rates that were used in surveillance and response. We found significant variation in reporting rates among districts, and district-specific rates of increase in reporting over time. Correcting time series of numbers of cases for variable reporting rates led, in some instances, to different estimates of the time-varying reproduction number of the epidemic, particularly outside the capital. Future analyses that compare Ebola transmission rates over time and across locations may be improved by considering the impacts of differential reporting rates. Epidemics are defined by a surge of cases of a disease, yet often a significant number of cases in an epidemic are never reported, for example because not all infected individuals have access to medical care. This underreporting can introduce bias into analyses of disease spread, by distorting patterns in where and when the most cases are observed. Conversely, quantifying underreporting can improve epidemic forecasts and containment strategies. In this study, we analyze data from the recent Ebola epidemic in West Africa, including the time, location and Ebola status of 6491 individual community burials, conducted over 25 weeks in four districts in Sierra Leone. We quantify how reporting rates varied over space and time, and show that estimates of transmission rates that are corrected for dynamic underreporting diverge significantly from uncorrected estimates, particularly earlier in the epidemic and outside the capital.
Collapse
|
91
|
Lau MSY, Gibson GJ, Adrakey H, McClelland A, Riley S, Zelner J, Streftaris G, Funk S, Metcalf J, Dalziel BD, Grenfell BT. A mechanistic spatio-temporal framework for modelling individual-to-individual transmission-With an application to the 2014-2015 West Africa Ebola outbreak. PLoS Comput Biol 2017; 13:e1005798. [PMID: 29084216 PMCID: PMC5679647 DOI: 10.1371/journal.pcbi.1005798] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2017] [Revised: 11/09/2017] [Accepted: 09/28/2017] [Indexed: 11/18/2022] Open
Abstract
In recent years there has been growing availability of individual-level spatio-temporal disease data, particularly due to the use of modern communicating devices with GPS tracking functionality. These detailed data have been proven useful for inferring disease transmission to a more refined level than previously. However, there remains a lack of statistically sound frameworks to model the underlying transmission dynamic in a mechanistic manner. Such a development is particularly crucial for enabling a general epidemic predictive framework at the individual level. In this paper we propose a new statistical framework for mechanistically modelling individual-to-individual disease transmission in a landscape with heterogeneous population density. Our methodology is first tested using simulated datasets, validating our inferential machinery. The methodology is subsequently applied to data that describes a regional Ebola outbreak in Western Africa (2014-2015). Our results show that the methods are able to obtain estimates of key epidemiological parameters that are broadly consistent with the literature, while revealing a significantly shorter distance of transmission. More importantly, in contrast to existing approaches, we are able to perform a more general model prediction that takes into account the susceptible population. Finally, our results show that, given reasonable scenarios, the framework can be an effective surrogate for susceptible-explicit individual models which are often computationally challenging. Availability of individual-level, spatio-temporal disease data (e.g. GPS locations of infected individuals) has been increasing in recent years, primarily due to the increased use of modern communication devices such as mobile phones. Such data create invaluable opportunities for researchers to study disease transmission on a more refined individual-to-individual level, facilitating the designs of potentially more effective control measures. However, the growing availability of such precise data has not been accompanied by development of statistically sound mechanistic frameworks. Developing such frameworks is an essential step for systematically extracting maximal information from data, in particular, evaluating the efficacy of individually-targeted control strategies and enabling forward epidemic prediction at the individual level. In this paper we develop a novel statistical framework that overcomes a few key limitations of existing approaches, enabling a machinery that can be used to infer the history of partially observed outbreaks and, more importantly, to produce a more comprehensive epidemic prediction. Our framework may also be a good surrogate for more computationally challenging individual-based models.
Collapse
|
92
|
Pitzer VE, Aguas R, Riley S, Loeffen WLA, Wood JLN, Grenfell BT. High turnover drives prolonged persistence of influenza in managed pig herds. J R Soc Interface 2017; 13:rsif.2016.0138. [PMID: 27358277 PMCID: PMC4938081 DOI: 10.1098/rsif.2016.0138] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Accepted: 06/08/2016] [Indexed: 11/16/2022] Open
Abstract
Pigs have long been hypothesized to play a central role in the emergence of novel human influenza A virus (IAV) strains, by serving as mixing vessels for mammalian and avian variants. However, the key issue of viral persistence in swine populations at different scales is ill understood. We address this gap using epidemiological models calibrated against seroprevalence data from Dutch finishing pigs to estimate the ‘critical herd size’ (CHS) for IAV persistence. We then examine the viral phylogenetic evidence for persistence by comparing human and swine IAV. Models suggest a CHS of approximately 3000 pigs above which influenza was likely to persist, i.e. orders of magnitude lower than persistence thresholds for IAV and other acute viruses in humans. At national and regional scales, we found much stronger empirical signatures of prolonged persistence of IAV in swine compared with human populations. These striking levels of persistence in small populations are driven by the high recruitment rate of susceptible piglets, and have significant implications for management of swine and for overall patterns of genetic diversity of IAV.
Collapse
|
93
|
Van Boeckel TP, Glennon EE, Chen D, Gilbert M, Robinson TP, Grenfell BT, Levin SA, Bonhoeffer S, Laxminarayan R. Reducing antimicrobial use in food animals. Science 2017; 357:1350-1352. [PMID: 28963240 PMCID: PMC6510296 DOI: 10.1126/science.aao1495] [Citation(s) in RCA: 335] [Impact Index Per Article: 47.9] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
|
94
|
Saad NJ, Bowles CC, Grenfell BT, Basnyat B, Arjyal A, Dongol S, Karkey A, Baker S, Pitzer VE. The impact of migration and antimicrobial resistance on the transmission dynamics of typhoid fever in Kathmandu, Nepal: A mathematical modelling study. PLoS Negl Trop Dis 2017; 11:e0005547. [PMID: 28475605 PMCID: PMC5435358 DOI: 10.1371/journal.pntd.0005547] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2016] [Revised: 05/17/2017] [Accepted: 04/03/2017] [Indexed: 11/21/2022] Open
Abstract
Background A substantial proportion of the global burden of typhoid fever occurs in South Asia. Kathmandu, Nepal experienced a substantial increase in the number of typhoid fever cases (caused by Salmonella Typhi) between 2000 and 2003, which subsequently declined but to a higher endemic level than in 2000. This epidemic of S. Typhi coincided with an increase in organisms with reduced susceptibility against fluoroquinolones, the emergence of S. Typhi H58, and an increase in the migratory population in Kathmandu. Methods We devised a mathematical model to investigate the potential epidemic drivers of typhoid in Kathmandu and fit this model to weekly data of S. Typhi cases between April 1997 and June 2011 and the age distribution of S. Typhi cases. We used this model to determine if the typhoid epidemic in Kathmandu was driven by heightened migration, the emergence of organisms with reduced susceptibility against fluoroquinolones or a combination of these factors. Results Models allowing for the migration of susceptible individuals into Kathmandu alone or in combination with the emergence of S. Typhi with reduced susceptibility against fluoroquinolones provided a good fit for the data. The emergence of organisms with reduced susceptibility against fluoroquinolones organisms alone, either through an increase in disease duration or increased transmission, did not fully explain the pattern of S. Typhi infections. Conclusions Our analysis is consistent with the hypothesis that the increase in typhoid fever in Kathmandu was associated with the migration of susceptible individuals into the city and aided by the emergence of reduced susceptibility against fluoroquinolones. These data support identifying and targeting migrant populations with typhoid immunization programmes to prevent transmission and disease. Typhoid fever is endemic in Nepal, with Kathmandu coined “the typhoid capital of the world”. We developed a mathematical model to assess the importance of migration and antimicrobial resistance on the transmission of typhoid fever in Kathmandu, Nepal from April 1997 to June 2011. During this period, the burden of typhoid fever increased markedly from January 2000 to December 2003, after which the epidemic declined, but to a higher endemic level than in 2000. Our findings are consistent with the hypothesis that migration of susceptible individuals into Kathmandu played an important role in the epidemic, and may have been further aided by the emergence of typhoid fever with reduced susceptibility against fluoroquinolones. This study showed that identifying and targeting migrant populations with control efforts could be an important avenue to prevent typhoid transmission and disease.
Collapse
|
95
|
Charu V, Zeger S, Gog J, Bjørnstad ON, Kissler S, Simonsen L, Grenfell BT, Viboud C. Human mobility and the spatial transmission of influenza in the United States. PLoS Comput Biol 2017; 13:e1005382. [PMID: 28187123 PMCID: PMC5349690 DOI: 10.1371/journal.pcbi.1005382] [Citation(s) in RCA: 115] [Impact Index Per Article: 16.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2016] [Revised: 03/14/2017] [Accepted: 01/26/2017] [Indexed: 11/18/2022] Open
Abstract
Seasonal influenza epidemics offer unique opportunities to study the invasion and re-invasion waves of a pathogen in a partially immune population. Detailed patterns of spread remain elusive, however, due to lack of granular disease data. Here we model high-volume city-level medical claims data and human mobility proxies to explore the drivers of influenza spread in the US during 2002–2010. Although the speed and pathways of spread varied across seasons, seven of eight epidemics likely originated in the Southern US. Each epidemic was associated with 1–5 early long-range transmission events, half of which sparked onward transmission. Gravity model estimates indicate a sharp decay in influenza transmission with the distance between infectious and susceptible cities, consistent with spread dominated by work commutes rather than air traffic. Two early-onset seasons associated with antigenic novelty had particularly localized modes of spread, suggesting that novel strains may spread in a more localized fashion than previously anticipated. The underlying mechanisms dictating the spatial spread of seasonal influenza remain poorly understood, in part because of the lack of spatially resolved disease data to quantify patterns of spread. In this paper, we address this issue by analyzing fine-grain insurance claims data on influenza-like-illnesses over eight seasons in ~300 locations throughout the United States. Using statistical methods, we found that seven of eight epidemics likely originated in the Southern US, that influenza spatial transmission is dominated by local traffic between cities, and that seasons marked by novel influenza virus circulation had a particularly radial, localized spatial structure. These findings are in stark contrast to prevailing theories of influenza spatial transmission that suggest that transmission is favored in low humidity environments and that spread is a dominated by air traffic between populous hubs.
Collapse
|
96
|
Metcalf CJE, Mina MJ, Winter AK, Grenfell BT. Opportunities and challenges of a World Serum Bank - Authors' reply. Lancet 2017; 389:252. [PMID: 28118913 DOI: 10.1016/s0140-6736(17)30054-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/20/2016] [Accepted: 09/26/2016] [Indexed: 11/17/2022]
|
97
|
Klepac P, Megiddo I, Grenfell BT, Laxminarayan R. Self-enforcing regional vaccination agreements. J R Soc Interface 2016; 13:20150907. [PMID: 26790996 PMCID: PMC4759795 DOI: 10.1098/rsif.2015.0907] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
In a highly interconnected world, immunizing infections are a transboundary problem, and their control and elimination require international cooperation and coordination. In the absence of a global or regional body that can impose a universal vaccination strategy, each individual country sets its own strategy. Mobility of populations across borders can promote free-riding, because a country can benefit from the vaccination efforts of its neighbours, which can result in vaccination coverage lower than the global optimum. Here we explore whether voluntary coalitions that reward countries that join by cooperatively increasing vaccination coverage can solve this problem. We use dynamic epidemiological models embedded in a game-theoretic framework in order to identify conditions in which coalitions are self-enforcing and therefore stable, and thus successful at promoting a cooperative vaccination strategy. We find that countries can achieve significantly greater vaccination coverage at a lower cost by forming coalitions than when acting independently, provided a coalition has the tools to deter free-riding. Furthermore, when economically or epidemiologically asymmetric countries form coalitions, realized coverage is regionally more consistent than in the absence of coalitions.
Collapse
|
98
|
Bharti N, Djibo A, Tatem AJ, Grenfell BT, Ferrari MJ. Measuring populations to improve vaccination coverage. Sci Rep 2016; 5:34541. [PMID: 27703191 PMCID: PMC5050518 DOI: 10.1038/srep34541] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2016] [Accepted: 09/14/2016] [Indexed: 11/09/2022] Open
Abstract
In low-income settings, vaccination campaigns supplement routine immunization but often fail to achieve coverage goals due to uncertainty about target population size and distribution. Accurate, updated estimates of target populations are rare but critical; short-term fluctuations can greatly impact population size and susceptibility. We use satellite imagery to quantify population fluctuations and the coverage achieved by a measles outbreak response vaccination campaign in urban Niger and compare campaign estimates to measurements from a post-campaign survey. Vaccine coverage was overestimated because the campaign underestimated resident numbers and seasonal migration further increased the target population. We combine satellite-derived measurements of fluctuations in population distribution with high-resolution measles case reports to develop a dynamic model that illustrates the potential improvement in vaccination campaign coverage if planners account for predictable population fluctuations. Satellite imagery can improve retrospective estimates of vaccination campaign impact and future campaign planning by synchronizing interventions with predictable population fluxes.
Collapse
|
99
|
Lessler J, Metcalf CJE, Cutts FT, Grenfell BT. Impact on Epidemic Measles of Vaccination Campaigns Triggered by Disease Outbreaks or Serosurveys: A Modeling Study. PLoS Med 2016; 13:e1002144. [PMID: 27727285 PMCID: PMC5058560 DOI: 10.1371/journal.pmed.1002144] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/13/2016] [Accepted: 09/02/2016] [Indexed: 11/23/2022] Open
Abstract
BACKGROUND Routine vaccination supplemented by planned campaigns occurring at 2-5 y intervals is the core of current measles control and elimination efforts. Yet, large, unexpected outbreaks still occur, even when control measures appear effective. Supplementing these activities with mass vaccination campaigns triggered when low levels of measles immunity are observed in a sample of the population (i.e., serosurveys) or incident measles cases occur may provide a way to limit the size of outbreaks. METHODS AND FINDINGS Measles incidence was simulated using stochastic age-structured epidemic models in settings conducive to high or low measles incidence, roughly reflecting demographic contexts and measles vaccination coverage of four heterogeneous countries: Nepal, Niger, Yemen, and Zambia. Uncertainty in underlying vaccination rates was modeled. Scenarios with case- or serosurvey-triggered campaigns reaching 20% of the susceptible population were compared to scenarios without triggered campaigns. The best performing of the tested case-triggered campaigns prevent an average of 28,613 (95% CI 25,722-31,505) cases over 15 y in our highest incidence setting and 599 (95% CI 464-735) cases in the lowest incidence setting. Serosurvey-triggered campaigns can prevent 89,173 (95% CI, 86,768-91,577) and 744 (612-876) cases, respectively, but are triggered yearly in high-incidence settings. Triggered campaigns reduce the highest cumulative incidence seen in simulations by up to 80%. While the scenarios considered in this strategic modeling exercise are reflective of real populations, the exact quantitative interpretation of the results is limited by the simplifications in country structure, vaccination policy, and surveillance system performance. Careful investigation into the cost-effectiveness in different contexts would be essential before moving forward with implementation. CONCLUSIONS Serologically triggered campaigns could help prevent severe epidemics in the face of epidemiological and vaccination uncertainty. Hence, small-scale serology may serve as the basis for effective adaptive public health strategies, although, in high-incidence settings, case-triggered approaches are likely more efficient.
Collapse
|
100
|
Morris SE, Zelner JL, Fauquier DA, Rowles TK, Rosel PE, Gulland F, Grenfell BT. Partially observed epidemics in wildlife hosts: modelling an outbreak of dolphin morbillivirus in the northwestern Atlantic, June 2013-2014. J R Soc Interface 2016; 12:rsif.2015.0676. [PMID: 26577594 PMCID: PMC4685842 DOI: 10.1098/rsif.2015.0676] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Morbilliviruses cause major mortality in marine mammals, but the dynamics of transmission and persistence are ill understood compared to terrestrial counterparts such as measles; this is especially true for epidemics in cetaceans. However, the recent outbreak of dolphin morbillivirus in the northwestern Atlantic Ocean can provide new insights into the epidemiology and spatio-temporal spread of this pathogen. To deal with uncertainties surrounding the ecology of this system (only stranded animals were observed), we develop a statistical framework that can extract key information about the underlying transmission process given only sparse data. Our self-exciting Poisson process model suggests that individuals are infectious for at most 24 days and can transfer infection up to two latitude degrees (220 km) within this time. In addition, the effective reproduction number is generally below one, but reaches 2.6 during a period of heightened stranding numbers near Virginia Beach, Virginia, in summer 2013. Network analysis suggests local movements dominate spatial spread, with seasonal migration facilitating wider dissemination along the coast. Finally, a low virus transmission rate or high levels of pre-existing immunity can explain the lack of viral spread into the Gulf of Mexico. More generally, our approach illustrates novel methodologies for analysing very indirectly observed epidemics.
Collapse
|