26
|
Evans MV, Drake JM, Jones L, Murdock CC. Assessing temperature-dependent competition between two invasive mosquito species. ECOLOGICAL APPLICATIONS : A PUBLICATION OF THE ECOLOGICAL SOCIETY OF AMERICA 2021; 31:e02334. [PMID: 33772946 DOI: 10.1002/eap.2334] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Revised: 12/10/2020] [Accepted: 01/14/2021] [Indexed: 06/12/2023]
Abstract
Invasive mosquitoes are expanding their ranges into new geographic areas and interacting with resident mosquito species. Understanding how novel interactions can affect mosquito population dynamics is necessary to predict transmission risk at invasion fronts. Mosquito life-history traits are extremely sensitive to temperature, and this can lead to temperature-dependent competition between competing invasive mosquito species. We explored temperature-dependent competition between Aedes aegypti and Anopheles stephensi, two invasive mosquito species whose distributions overlap in India, the Middle East, and North Africa, where An. stephensi is currently expanding into the endemic range of Ae. aegypti. We followed mosquito cohorts raised at different intraspecific and interspecific densities across five temperatures (16-32°C) to measure traits relevant for population growth and to estimate species' per capita growth rates. We then used these growth rates to derive each species' competitive ability at each temperature. We find strong evidence for asymmetric competition at all temperatures, with Ae. aegypti emerging as the dominant competitor. This was primarily because of differences in larval survival and development times across all temperatures that resulted in a higher estimated intrinsic growth rate and competitive tolerance estimate for Ae. aegypti compared to An. stephensi. The spread of An. stephensi into the African continent could lead to urban transmission of malaria, an otherwise rural disease, increasing the human population at risk and complicating malaria elimination efforts. Competition has resulted in habitat segregation of other invasive mosquito species, and our results suggest that it may play a role in determining the distribution of An. stephensi across its invasive range.
Collapse
|
27
|
Lodge EK, Schatz AM, Drake JM. Protective population behavior change in outbreaks of emerging infectious disease. BMC Infect Dis 2021; 21:577. [PMID: 34130652 PMCID: PMC8205197 DOI: 10.1186/s12879-021-06299-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Accepted: 06/09/2021] [Indexed: 11/25/2022] Open
Abstract
BACKGROUND During outbreaks of emerging and re-emerging infections, the lack of effective drugs and vaccines increases reliance on non-pharmacologic public health interventions and behavior change to limit human-to-human transmission. Interventions that increase the speed with which infected individuals remove themselves from the susceptible population are paramount, particularly isolation and hospitalization. Ebola virus disease (EVD), Severe Acute Respiratory Syndrome (SARS), and Middle East Respiratory Syndrome (MERS) are zoonotic viruses that have caused significant recent outbreaks with sustained human-to-human transmission. METHODS This investigation quantified changing mean removal rates (MRR) and days from symptom onset to hospitalization (DSOH) of infected individuals from the population in seven different outbreaks of EVD, SARS, and MERS, to test for statistically significant differences in these metrics between outbreaks. RESULTS We found that epidemic week and viral serial interval were correlated with the speed with which populations developed and maintained health behaviors in each outbreak. CONCLUSIONS These findings highlight intrinsic population-level changes in isolation rates in multiple epidemics of three zoonotic infections with established human-to-human transmission and significant morbidity and mortality. These data are particularly useful for disease modelers seeking to forecast the spread of emerging pathogens.
Collapse
|
28
|
Sánchez CA, Venkatachalam-Vaz J, Drake JM. Spillover of zoonotic pathogens: A review of reviews. Zoonoses Public Health 2021; 68:563-577. [PMID: 34018336 DOI: 10.1111/zph.12846] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Revised: 02/22/2021] [Accepted: 04/03/2021] [Indexed: 11/30/2022]
Abstract
Zoonotic spillover and subsequent disease emergence cause significant, long-lasting impacts on our social, economic, environmental and political systems. Identifying and averting spillover transmission is crucial for preventing outbreaks and mitigating infectious disease burdens. Investigating the processes that lead to spillover fundamentally involves interactions between animals, humans, pathogens and the environments they inhabit. Accordingly, it is recognized that transdisciplinary approaches provide a more holistic understanding of spillover phenomena. To characterize the discourse about spillover within and between disciplines, we conducted a review of review papers about spillover from multiple disciplines. We systematically searched and screened literature from several databases to identify a corpus of review papers from ten academic disciplines. We performed qualitative content analysis on text where authors described either a spillover pathway, or a conceptual gap in spillover theory. Cluster analysis of pathway data identified nine major spillover processes discussed in the review literature. We summarized the main features of each process, how different disciplines contributed to them, and identified specialist and generalist disciplines based on the breadth of processes they studied. Network analyses showed strong similarities between concepts reviewed by 'One Health' disciplines (e.g. Veterinary Science & Animal Health, Public Health & Medicine, Ecology & Evolution, Environmental Science), which had broad conceptual scope and were well-connected to other disciplines. By contrast, awas focused on processes that are relatively overlooked by other disciplines, especially those involving food behaviour and livestock husbandry practices. Virology and Cellular & Molecular Biology were narrower in scope, primarily focusing on concepts related to adaption and evolution of zoonotic viruses. Finally, we identified priority areas for future research into zoonotic spillover by studying the gap data.
Collapse
|
29
|
Drake JM, Dahlin K, Rohani P, Handel A. Five approaches to the suppression of SARS-CoV-2 without intensive social distancing. Proc Biol Sci 2021; 288:20203074. [PMID: 33906405 PMCID: PMC8080008 DOI: 10.1098/rspb.2020.3074] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Accepted: 03/27/2021] [Indexed: 02/06/2023] Open
Abstract
Initial efforts to mitigate transmission of SARS-CoV-2 relied on intensive social distancing measures such as school and workplace closures, shelter-in-place orders and prohibitions on the gathering of people. Other non-pharmaceutical interventions for suppressing transmission include active case finding, contact tracing, quarantine, immunity or health certification, and a wide range of personal protective measures. Here we investigate the potential effectiveness of these alternative approaches to suppression. We introduce a conceptual framework represented by two mathematical models that differ in strategy. We find both strategies may be effective, although both require extensive testing and work within a relatively narrow range of conditions. Generalized protective measures such as wearing face masks, improved hygiene and local reductions in density are found to significantly increase the effectiveness of targeted interventions.
Collapse
|
30
|
Evans MV, Bonds MH, Cordier LF, Drake JM, Ihantamalala F, Haruna J, Miller AC, Murdock CC, Randriamanambtsoa M, Raza-Fanomezanjanahary EM, Razafinjato BR, Garchitorena AC. Socio-demographic, not environmental, risk factors explain fine-scale spatial patterns of diarrhoeal disease in Ifanadiana, rural Madagascar. Proc Biol Sci 2021; 288:20202501. [PMID: 33653145 PMCID: PMC7934917 DOI: 10.1098/rspb.2020.2501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
Precision health mapping is a technique that uses spatial relationships between socio-ecological variables and disease to map the spatial distribution of disease, particularly for diseases with strong environmental signatures, such as diarrhoeal disease (DD). While some studies use GPS-tagged location data, other precision health mapping efforts rely heavily on data collected at coarse-spatial scales and may not produce operationally relevant predictions at fine enough spatio-temporal scales to inform local health programmes. We use two fine-scale health datasets collected in a rural district of Madagascar to identify socio-ecological covariates associated with childhood DD. We constructed generalized linear mixed models including socio-demographic, climatic and landcover variables and estimated variable importance via multi-model inference. We find that socio-demographic variables, and not environmental variables, are strong predictors of the spatial distribution of disease risk at both individual and commune-level (cluster of villages) spatial scales. Climatic variables predicted strong seasonality in DD, with the highest incidence in colder, drier months, but did not explain spatial patterns. Interestingly, the occurrence of a national holiday was highly predictive of increased DD incidence, highlighting the need for including cultural factors in modelling efforts. Our findings suggest that precision health mapping efforts that do not include socio-demographic covariates may have reduced explanatory power at the local scale. More research is needed to better define the set of conditions under which the application of precision health mapping can be operationally useful to local public health professionals.
Collapse
|
31
|
Hackett EJ, Leahey E, Parker JN, Rafols I, Hampton SE, Corte U, Chavarro D, Drake JM, Penders B, Sheble L, Vermeulen N, Vision TJ. Do synthesis centers synthesize? A semantic analysis of topical diversity in research. RESEARCH POLICY 2021; 50:104069. [PMID: 33390628 PMCID: PMC7695893 DOI: 10.1016/j.respol.2020.104069] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Revised: 11/26/2019] [Accepted: 07/01/2020] [Indexed: 11/18/2022]
Abstract
Synthesis centers are a form of scientific organization that catalyzes and supports research that integrates diverse theories, methods and data across spatial or temporal scales to increase the generality, parsimony, applicability, or empirical soundness of scientific explanations. Synthesis working groups are a distinctive form of scientific collaboration that produce consequential, high-impact publications. But no one has asked if synthesis working groups synthesize: are their publications substantially more diverse than others, and if so, in what ways and with what effect? We investigate these questions by using Latent Dirichlet Analysis to compare the topical diversity of papers published by synthesis center collaborations with that of papers in a reference corpus. Topical diversity was operationalized and measured in several ways, both to reflect aggregate diversity and to emphasize particular aspects of diversity (such as variety, evenness, and balance). Synthesis center publications have greater topical variety and evenness, but less disparity, than do papers in the reference corpus. The influence of synthesis center origins on aspects of diversity is only partly mediated by the size and heterogeneity of collaborations: when taking into account the numbers of authors, distinct institutions, and references, synthesis center origins retain a significant direct effect on diversity measures. Controlling for the size and heterogeneity of collaborative groups, synthesis center origins and diversity measures significantly influence the visibility of publications, as indicated by citation measures. We conclude by suggesting social processes within collaborations that might account for the observed effects, by inviting further exploration of what this novel textual analysis approach might reveal about interdisciplinary research, and by offering some practical implications of our results.
Collapse
|
32
|
Evans MV, Garchitorena A, Rakotonanahary RJL, Drake JM, Andriamihaja B, Rajaonarifara E, Ngonghala CN, Roche B, Bonds MH, Rakotonirina J. Reconciling model predictions with low reported cases of COVID-19 in Sub-Saharan Africa: insights from Madagascar. Glob Health Action 2020; 13:1816044. [PMID: 33012269 PMCID: PMC7580764 DOI: 10.1080/16549716.2020.1816044] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 08/24/2020] [Indexed: 12/14/2022] Open
Abstract
COVID-19 has wreaked havoc globally with particular concerns for sub-Saharan Africa (SSA), where models suggest that the majority of the population will become infected. Conventional wisdom suggests that the continent will bear a higher burden of COVID-19 for the same reasons it suffers from other infectious diseases: ecology, socio-economic conditions, lack of water and sanitation infrastructure, and weak health systems. However, so far SSA has reported lower incidence and fatalities compared to the predictions of standard models and the experience of other regions of the world. There are three leading explanations, each with different implications for the final epidemic burden: (1) low case detection, (2) differences in epidemiology (e.g. low R 0 ), and (3) policy interventions. The low number of cases have led some SSA governments to relaxing these policy interventions. Will this result in a resurgence of cases? To understand how to interpret the lower-than-expected COVID-19 case data in Madagascar, we use a simple age-structured model to explore each of these explanations and predict the epidemic impact associated with them. We show that the incidence of COVID-19 cases as of July 2020 can be explained by any combination of the late introduction of first imported cases, early implementation of non-pharmaceutical interventions (NPIs), and low case detection rates. We then re-evaluate these findings in the context of the COVID-19 epidemic in Madagascar through August 2020. This analysis reinforces that Madagascar, along with other countries in SSA, remains at risk of a growing health crisis. If NPIs remain enforced, up to 50,000 lives may be saved. Even with NPIs, without vaccines and new therapies, COVID-19 could infect up to 30% of the population, making it the largest public health threat in Madagascar for the coming year, hence the importance of clinical trials and continually improving access to healthcare.
Collapse
|
33
|
O'Regan SM, O'Dea EB, Rohani P, Drake JM. Transient indicators of tipping points in infectious diseases. J R Soc Interface 2020; 17:20200094. [PMID: 32933375 DOI: 10.1098/rsif.2020.0094] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The majority of known early warning indicators of critical transitions rely on asymptotic resilience and critical slowing down. In continuous systems, critical slowing down is mathematically described by a decrease in magnitude of the dominant eigenvalue of the Jacobian matrix on the approach to a critical transition. Here, we show that measures of transient dynamics, specifically, reactivity and the maximum of the amplification envelope, also change systematically as a bifurcation is approached in an important class of models for epidemics of infectious diseases. Furthermore, we introduce indicators designed to detect trends in these measures and find that they reliably classify time series of case notifications simulated from stochastic models according to levels of vaccine uptake. Greater attention should be focused on the potential for systems to exhibit transient amplification of perturbations as a critical threshold is approached, and should be considered when searching for generic leading indicators of tipping points. Awareness of this phenomenon will enrich understanding of the dynamics of complex systems on the verge of a critical transition.
Collapse
|
34
|
Drake JM, Rohani P, Dahlin K, Handel A. Five approaches to the suppression of SARS-CoV-2 without intensive social distancing. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2020:2020.07.30.20165159. [PMID: 32766603 PMCID: PMC7402061 DOI: 10.1101/2020.07.30.20165159] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Initial efforts to mitigate transmission of SARS-CoV-2 relied on intensive social distancing measures such as school and workplace closures, shelter-in-place orders, and prohibitions on the gathering of people. Other non-pharmaceutical interventions for suppressing transmission include active case finding, contact tracing, quarantine, immunity or health certification, and a wide range of personal protective measures. Here we investigate the potential effectiveness of these alternative approaches to suppression. We introduce a conceptual framework represented by two mathematical models that differ in strategy. We find both strategies may be effective, although both require extensive testing and work within a relatively narrow range of conditions. Generalized protective measures such as wearing face masks, improved hygiene, and local reductions in density are found to significantly increase the effectiveness of targeted interventions.
Collapse
|
35
|
Han BA, O'Regan SM, Paul Schmidt J, Drake JM. Integrating data mining and transmission theory in the ecology of infectious diseases. Ecol Lett 2020; 23:1178-1188. [PMID: 32441459 PMCID: PMC7384120 DOI: 10.1111/ele.13520] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Revised: 01/21/2020] [Accepted: 03/27/2020] [Indexed: 01/07/2023]
Abstract
Our understanding of ecological processes is built on patterns inferred from data. Applying modern analytical tools such as machine learning to increasingly high dimensional data offers the potential to expand our perspectives on these processes, shedding new light on complex ecological phenomena such as pathogen transmission in wild populations. Here, we propose a novel approach that combines data mining with theoretical models of disease dynamics. Using rodents as an example, we incorporate statistical differences in the life history features of zoonotic reservoir hosts into pathogen transmission models, enabling us to bound the range of dynamical phenomena associated with hosts, based on their traits. We then test for associations between equilibrium prevalence, a key epidemiological metric and data on human outbreaks of rodent-borne zoonoses, identifying matches between empirical evidence and theoretical predictions of transmission dynamics. We show how this framework can be generalized to other systems through a rubric of disease models and parameters that can be derived from empirical data. By linking life history components directly to their effects on disease dynamics, our mining-modelling approach integrates machine learning and theoretical models to explore mechanisms in the macroecology of pathogen transmission and their consequences for spillover infection to humans.
Collapse
|
36
|
Walker JW, Kittur N, Binder S, Castleman JD, Drake JM, Campbell CH, King CH, Colley DG. Environmental Predictors of Schistosomiasis Persistent Hotspots following Mass Treatment with Praziquantel. Am J Trop Med Hyg 2020; 102:328-338. [PMID: 31889506 PMCID: PMC7008331 DOI: 10.4269/ajtmh.19-0658] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Schistosomiasis control programs rely heavily on mass drug administration (MDA) campaigns with praziquantel for preventative chemotherapy. Areas where the prevalence and/or intensity of schistosomiasis infection remains high even after several rounds of treatment, termed "persistent hotspots" (PHSs), have been identified in trials of MDA effectiveness conducted by the Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) in Kenya, Mozambique, Tanzania, and Côte d'Ivoire. In this analysis, we apply a previously developed set of criteria to classify the PHS status of 531 study villages from five SCORE trials. We then fit logistic regression models to data from SCORE and publically available georeferenced datasets to evaluate the influence of local environmental and population features, pre-intervention infection burden, and treatment scheduling on PHS status in each trial. The frequency of PHS in individual trials ranged from 35.3% to 71.6% in study villages. Significant relationships between PHS status and MDA frequency, distance to freshwater, rainfall, baseline schistosomiasis burden, elevation, land cover type, and village remoteness were each observed in at least one trial, although the strength and direction of these relationships was not always consistent among study sites. These findings suggest that PHSs are driven in part by environmental conditions that modify the risk and frequency of reinfection.
Collapse
|
37
|
Drake JM, O’Regan SM, Dakos V, Kéfi S, Rohani P. Alternative stable states, tipping points, and early warning signals of ecological transitions. THEOR ECOL-NETH 2020. [DOI: 10.1093/oso/9780198824282.003.0015] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Ecological systems are prone to dramatic shifts between alternative stable states. In reality, these shifts are often caused by slow forces external to the system that eventually push it over a tipping point. Theory predicts that when ecological systems are brought close to a tipping point, the dynamical feedback intrinsic to the system interact with intrinsic noise and extrinsic perturbations in characteristic ways. The resulting phenomena thus serve as “early warning signals” for shifts such as population collapse. In this chapter, we review the basic (qualitative) theory of such systems. We then illustrate the main ideas with a series of models that both represent fundamental ecological ideas (e.g. density-dependence) and are amenable to mathematical analysis. These analyses provide theoretical predictions about the nature of measurable fluctuations in the vicinity of a tipping point. We conclude with a review of empirical evidence from laboratory microcosms, field manipulations, and observational studies.
Collapse
|
38
|
Bock SL, Lowers RH, Rainwater TR, Stolen E, Drake JM, Wilkinson PM, Weiss S, Back B, Guillette L, Parrott BB. Spatial and temporal variation in nest temperatures forecasts sex ratio skews in a crocodilian with environmental sex determination. Proc Biol Sci 2020; 287:20200210. [PMID: 32345164 DOI: 10.1098/rspb.2020.0210] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Species displaying temperature-dependent sex determination (TSD) are especially vulnerable to the effects of a rapidly changing global climate due to their profound sensitivity to thermal cues during development. Predicting the consequences of climate change for these species, including skewed offspring sex ratios, depends on understanding how climatic factors interface with features of maternal nesting behaviour to shape the developmental environment. Here, we measure thermal profiles in 86 nests at two geographically distinct sites in the northern and southern regions of the American alligator's (Alligator mississippiensis) geographical range, and examine the influence of both climatic factors and maternally driven nest characteristics on nest temperature variation. Changes in daily maximum air temperatures drive annual trends in nest temperatures, while variation in individual nest temperatures is also related to local habitat factors and microclimate characteristics. Without any compensatory nesting behaviours, nest temperatures are projected to increase by 1.6-3.7°C by the year 2100, and these changes are predicted to have dramatic consequences for offspring sex ratios. Exact sex ratio outcomes vary widely depending on site and emission scenario as a function of the unique temperature-by-sex reaction norm exhibited by all crocodilians. By revealing the ecological drivers of nest temperature variation in the American alligator, this study provides important insights into the potential consequences of climate change for crocodilian species, many of which are already threatened by extinction.
Collapse
|
39
|
Evans MV, Hintz CW, Jones L, Shiau J, Solano N, Drake JM, Murdock CC. Microclimate and Larval Habitat Density Predict Adult Aedes albopictus Abundance in Urban Areas. Am J Trop Med Hyg 2020; 101:362-370. [PMID: 31190685 PMCID: PMC6685558 DOI: 10.4269/ajtmh.19-0220] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
The Asian tiger mosquito, Aedes albopictus, transmits several arboviruses of public health importance, including chikungunya and dengue. Since its introduction to the United States in 1985, the species has invaded more than 40 states, including temperate areas not previously at risk of Aedes-transmitted arboviruses. Mathematical models incorporate climatic variables in predictions of site-specific Ae. albopictus abundances to identify human populations at risk of disease. However, these models rely on coarse resolutions of environmental data that may not accurately represent the climatic profile experienced by mosquitoes in the field, particularly in climatically heterogeneous urban areas. In this study, we pair field surveys of larval and adult Ae. albopictus mosquitoes with site-specific microclimate data across a range of land use types to investigate the relationships between microclimate, density of larval habitat, and adult mosquito abundance and determine whether these relationships change across an urban gradient. We find no evidence for a difference in larval habitat density or adult abundance between rural, suburban, and urban land classes. Adult abundance increases with increasing larval habitat density, which itself is dependent on microclimate. Adult abundance is strongly explained by microclimate variables, demonstrating that theoretically derived, laboratory-parameterized relationships in ectotherm physiology apply to the field. Our results support the continued use of temperature-dependent models to predict Ae. albopictus abundance in urban areas.
Collapse
|
40
|
Harris MJ, Hay SI, Drake JM. Early warning signals of malaria resurgence in Kericho, Kenya. Biol Lett 2020; 16:20190713. [PMID: 32183637 PMCID: PMC7115183 DOI: 10.1098/rsbl.2019.0713] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Accepted: 02/19/2020] [Indexed: 11/12/2022] Open
Abstract
Campaigns to eliminate infectious diseases could be greatly aided by methods for providing early warning signals of resurgence. Theory predicts that as a disease transmission system undergoes a transition from stability at the disease-free equilibrium to sustained transmission, it will exhibit characteristic behaviours known as critical slowing down, referring to the speed at which fluctuations in the number of cases are dampened, for instance the extinction of a local transmission chain after infection from an imported case. These phenomena include increases in several summary statistics, including lag-1 autocorrelation, variance and the first difference of variance. Here, we report the first empirical test of this prediction during the resurgence of malaria in Kericho, Kenya. For 10 summary statistics, we measured the approach to criticality in a rolling window to quantify the size of effect and directions. Nine of the statistics increased as predicted and variance, the first difference of variance, autocovariance, lag-1 autocorrelation and decay time returned early warning signals of critical slowing down based on permutation tests. These results show that time series of disease incidence collected through ordinary surveillance activities may exhibit characteristic signatures prior to an outbreak, a phenomenon that may be quite general among infectious disease systems.
Collapse
|
41
|
Brett T, Ajelli M, Liu QH, Krauland MG, Grefenstette JJ, van Panhuis WG, Vespignani A, Drake JM, Rohani P. Detecting critical slowing down in high-dimensional epidemiological systems. PLoS Comput Biol 2020; 16:e1007679. [PMID: 32150536 PMCID: PMC7082051 DOI: 10.1371/journal.pcbi.1007679] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Revised: 03/19/2020] [Accepted: 01/23/2020] [Indexed: 01/05/2023] Open
Abstract
Despite medical advances, the emergence and re-emergence of infectious diseases continue to pose a public health threat. Low-dimensional epidemiological models predict that epidemic transitions are preceded by the phenomenon of critical slowing down (CSD). This has raised the possibility of anticipating disease (re-)emergence using CSD-based early-warning signals (EWS), which are statistical moments estimated from time series data. For EWS to be useful at detecting future (re-)emergence, CSD needs to be a generic (model-independent) feature of epidemiological dynamics irrespective of system complexity. Currently, it is unclear whether the predictions of CSD-derived from simple, low-dimensional systems-pertain to real systems, which are high-dimensional. To assess the generality of CSD, we carried out a simulation study of a hierarchy of models, with increasing structural complexity and dimensionality, for a measles-like infectious disease. Our five models included: i) a nonseasonal homogeneous Susceptible-Exposed-Infectious-Recovered (SEIR) model, ii) a homogeneous SEIR model with seasonality in transmission, iii) an age-structured SEIR model, iv) a multiplex network-based model (Mplex) and v) an agent-based simulator (FRED). All models were parameterised to have a herd-immunity immunization threshold of around 90% coverage, and underwent a linear decrease in vaccine uptake, from 92% to 70% over 15 years. We found evidence of CSD prior to disease re-emergence in all models. We also evaluated the performance of seven EWS: the autocorrelation, coefficient of variation, index of dispersion, kurtosis, mean, skewness, variance. Performance was scored using the Area Under the ROC Curve (AUC) statistic. The best performing EWS were the mean and variance, with AUC > 0.75 one year before the estimated transition time. These two, along with the autocorrelation and index of dispersion, are promising candidate EWS for detecting disease emergence.
Collapse
|
42
|
O'Dea EB, Park AW, Drake JM. Estimating the distance to an epidemic threshold. J R Soc Interface 2019; 15:rsif.2018.0034. [PMID: 29950512 PMCID: PMC6030631 DOI: 10.1098/rsif.2018.0034] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2018] [Accepted: 05/31/2018] [Indexed: 12/20/2022] Open
Abstract
The epidemic threshold of the susceptible-infected-recovered model is a boundary separating parameters that permit epidemics from those that do not. This threshold corresponds to parameters where the system's equilibrium becomes unstable. Consequently, we use the average rate at which deviations from the equilibrium shrink to define a distance to this threshold. However, the vital dynamics of the host population may occur slowly even when transmission is far from threshold levels. Here, we show analytically how such slow dynamics can prevent estimation of the distance to the threshold from fluctuations in the susceptible population. Although these results are exact only in the limit of long-term observation of a large system, simulations show that they still provide useful insight into systems with a range of population sizes, environmental noise and observation schemes. Having established some guidelines about when estimates are accurate, we then illustrate how multiple distance estimates can be used to estimate the rate of approach to the threshold. The estimation approach is general and may be applicable to zoonotic pathogens such as Middle East respiratory syndrome-related coronavirus (MERS-CoV) as well as vaccine-preventable diseases like measles.
Collapse
|
43
|
Schmidt JP, Maher S, Drake JM, Huang T, Farrell MJ, Han BA. Ecological indicators of mammal exposure to Ebolavirus. Philos Trans R Soc Lond B Biol Sci 2019; 374:20180337. [PMID: 31401967 PMCID: PMC6711296 DOI: 10.1098/rstb.2018.0337] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Much of the basic ecology of Ebolavirus remains unresolved despite accumulating disease outbreaks, viral strains and evidence of animal hosts. Because human Ebolavirus epidemics have been linked to contact with wild mammals other than bats, traits shared by species that have been infected by Ebolavirus and their phylogenetic distribution could suggest ecological mechanisms contributing to human Ebolavirus spillovers. We compiled data on Ebolavirus exposure in mammals and corresponding data on life-history traits, movement, and diet, and used boosted regression trees (BRT) to identify predictors of exposure and infection for 119 species (hereafter hosts). Mapping the phylogenetic distribution of presumptive Ebolavirus hosts reveals that they are scattered across several distinct mammal clades, but concentrated among Old World fruit bats, primates and artiodactyls. While sampling effort was the most important predictor, explaining nearly as much of the variation among hosts as traits, BRT models distinguished hosts from all other species with greater than 97% accuracy, and revealed probable Ebolavirus hosts as large-bodied, frugivorous, and with slow life histories. Provisionally, results suggest that some insectivorous bat genera, Old World monkeys and forest antelopes should receive priority in Ebolavirus survey efforts. This article is part of the theme issue ‘Dynamic and integrative approaches to understanding pathogen spillover’.
Collapse
|
44
|
Kramer AM, Teitelbaum CS, Griffin A, Drake JM. Multiscale model of regional population decline in little brown bats due to white-nose syndrome. Ecol Evol 2019; 9:8639-8651. [PMID: 31410268 PMCID: PMC6686297 DOI: 10.1002/ece3.5405] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2019] [Accepted: 05/11/2019] [Indexed: 01/26/2023] Open
Abstract
The introduced fungal pathogen Pseudogymnoascus destructans is causing decline of several species of bats in North America, with some even at risk of extinction or extirpation. The severity of the epidemic of white-nose syndrome caused by P. destructans has prompted investigation of the transmission and virulence of infection at multiple scales, but linking these scales is necessary to quantify the mechanisms of transmission and assess population-scale declines.We built a model connecting within-hibernaculum disease dynamics of little brown bats to regional-scale dispersal, reproduction, and disease spread, including multiple plausible mechanisms of transmission.We parameterized the model using the approach of plausible parameter sets, by comparing stochastic simulation results to statistical probes from empirical data on within-hibernaculum prevalence and survival, as well as among-hibernacula spread across a region.Our results are consistent with frequency-dependent transmission between bats, support an important role of environmental transmission, and show very little effect of dispersal among colonies on metapopulation survival.The results help identify the influential parameters and largest sources of uncertainty. The model also offers a generalizable method to assess hypotheses about hibernaculum-to-hibernaculum transmission and to identify gaps in knowledge about key processes, and could be expanded to include additional mechanisms or bat species.
Collapse
|
45
|
Drake JM, Brett TS, Chen S, Epureanu BI, Ferrari MJ, Marty É, Miller PB, O’Dea EB, O’Regan SM, Park AW, Rohani P. The statistics of epidemic transitions. PLoS Comput Biol 2019; 15:e1006917. [PMID: 31067217 PMCID: PMC6505855 DOI: 10.1371/journal.pcbi.1006917] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Emerging and re-emerging pathogens exhibit very complex dynamics, are hard to model and difficult to predict. Their dynamics might appear intractable. However, new statistical approaches-rooted in dynamical systems and the theory of stochastic processes-have yielded insight into the dynamics of emerging and re-emerging pathogens. We argue that these approaches may lead to new methods for predicting epidemics. This perspective views pathogen emergence and re-emergence as a "critical transition," and uses the concept of noisy dynamic bifurcation to understand the relationship between the system observables and the distance to this transition. Because the system dynamics exhibit characteristic fluctuations in response to perturbations for a system in the vicinity of a critical point, we propose this information may be harnessed to develop early warning signals. Specifically, the motion of perturbations slows as the system approaches the transition.
Collapse
|
46
|
Abstract
Second-order statistics such as the variance and autocorrelation can be useful indicators of the stability of randomly perturbed systems, in some cases providing early warning of an impending, dramatic change in the system’s dynamics. One specific application area of interest is the surveillance of infectious diseases. In the context of disease (re-)emergence, a goal could be to have an indicator that is informative of whether the system is approaching the epidemic threshold, a point beyond which a major outbreak becomes possible. Prior work in this area has provided some proof of this principle but has not analytically treated the effect of imperfect observation on the behavior of indicators. This work provides expected values for several moments of the number of reported cases, where reported cases follow a binomial or negative binomial distribution with a mean based on the number of deaths in a birth-death-immigration process over some reporting interval. The normalized second factorial moment and the decay time of the number of reported cases are two indicators that are insensitive to the reporting probability. Simulation is used to show how this insensitivity could be used to distinguish a trend of increased reporting from a trend of increased transmission. The simulation study also illustrates both the high variance of estimates and the possibility of reducing the variance by averaging over an ensemble of estimates from multiple time series.
Collapse
|
47
|
Chen S, O'Dea EB, Drake JM, Epureanu BI. Eigenvalues of the covariance matrix as early warning signals for critical transitions in ecological systems. Sci Rep 2019; 9:2572. [PMID: 30796264 PMCID: PMC6385210 DOI: 10.1038/s41598-019-38961-5] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 12/27/2018] [Indexed: 11/23/2022] Open
Abstract
Many ecological systems are subject critical transitions, which are abrupt changes to contrasting states triggered by small changes in some key component of the system. Temporal early warning signals such as the variance of a time series, and spatial early warning signals such as the spatial correlation in a snapshot of the system's state, have been proposed to forecast critical transitions. However, temporal early warning signals do not take the spatial pattern into account, and past spatial indicators only examine one snapshot at a time. In this study, we propose the use of eigenvalues of the covariance matrix of multiple time series as early warning signals. We first show theoretically why these indicators may increase as the system moves closer to the critical transition. Then, we apply the method to simulated data from several spatial ecological models to demonstrate the method's applicability. This method has the advantage that it takes into account only the fluctuations of the system about its equilibrium, thus eliminating the effects of any change in equilibrium values. The eigenvector associated with the largest eigenvalue of the covariance matrix is helpful for identifying the regions that are most vulnerable to the critical transition.
Collapse
|
48
|
Park AW, Farrell MJ, Schmidt JP, Huang S, Dallas TA, Pappalardo P, Drake JM, Stephens PR, Poulin R, Nunn CL, Davies TJ. Characterizing the phylogenetic specialism-generalism spectrum of mammal parasites. Proc Biol Sci 2019. [PMID: 29514973 DOI: 10.1098/rspb.2017.2613] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
The distribution of parasites across mammalian hosts is complex and represents a differential ability or opportunity to infect different host species. Here, we take a macroecological approach to investigate factors influencing why some parasites show a tendency to infect species widely distributed in the host phylogeny (phylogenetic generalism) while others infect only closely related hosts. Using a database on over 1400 parasite species that have been documented to infect up to 69 terrestrial mammal host species, we characterize the phylogenetic generalism of parasites using standard effect sizes for three metrics: mean pairwise phylogenetic distance (PD), maximum PD and phylogenetic aggregation. We identify a trend towards phylogenetic specialism, though statistically host relatedness is most often equivalent to that expected from a random sample of host species. Bacteria and arthropod parasites are typically the most generalist, viruses and helminths exhibit intermediate generalism, and protozoa are on average the most specialist. While viruses and helminths have similar mean pairwise PD on average, the viruses exhibit higher variation as a group. Close-contact transmission is the transmission mode most associated with specialism. Most parasites exhibiting phylogenetic aggregation (associating with discrete groups of species dispersed across the host phylogeny) are helminths and viruses.
Collapse
|
49
|
Walker JW, Han BA, Ott IM, Drake JM. Transmissibility of emerging viral zoonoses. PLoS One 2018; 13:e0206926. [PMID: 30403733 PMCID: PMC6221319 DOI: 10.1371/journal.pone.0206926] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2018] [Accepted: 10/17/2018] [Indexed: 01/23/2023] Open
Abstract
Effective public health research and preparedness requires an accurate understanding of which virus species possess or are at risk of developing human transmissibility. Unfortunately, our ability to identify these viruses is limited by gaps in disease surveillance and an incomplete understanding of the process of viral adaptation. By fitting boosted regression trees to data on 224 human viruses and their associated traits, we developed a model that predicts the human transmission ability of zoonotic viruses with over 84% accuracy. This model identifies several viruses that may have an undocumented capacity for transmission between humans. Viral traits that predicted human transmissibility included infection of nonhuman primates, the absence of a lipid envelope, and detection in the human nervous system and respiratory tract. This predictive model can be used to prioritize high-risk viruses for future research and surveillance, and could inform an integrated early warning system for emerging infectious diseases.
Collapse
|
50
|
|