An open-access database of infectious disease transmission trees to explore superspreader epidemiology.

Historically, emerging and reemerging infectious diseases have caused large, deadly, and expensive multinational outbreaks. Often outbreak investigations aim to identify who infected whom by reconstructing the outbreak transmission tree, which visualizes transmission between individuals as a network with nodes representing individuals and branches representing transmission from person to person. We compiled a database, called OutbreakTrees, of 382 published, standardized transmission trees consisting of 16 directly transmitted diseases ranging in size from 2 to 286 cases. For each tree and disease, we calculated several key statistics, such as tree size, average number of secondary infections, the dispersion parameter, and the proportion of cases considered superspreaders, and examined how these statistics varied over the course of each outbreak and under different assumptions about the completeness of outbreak investigations. We demonstrated the potential utility of the database through 2 short analyses addressing questions about superspreader epidemiology for a variety of diseases, including Coronavirus Disease 2019 (COVID-19). First, we found that our transmission trees were consistent with theory predicting that intermediate dispersion parameters give rise to the highest proportion of cases causing superspreading events. Additionally, we investigated patterns in how superspreaders are infected. Across trees with more than 1 superspreader, we found preliminary support for the theory that superspreaders generate other superspreaders. In sum, our findings put the role of superspreading in COVID-19 transmission in perspective with that of other diseases and suggest an approach to further research regarding the generation of superspreaders. These data have been made openly available to encourage reuse and further scientific inquiry.

Challenges of COVID-19 Case Forecasting in the US, 2020-2021.

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub (https://covid19forecasthub.org). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1-4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naïve baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making.

A data-driven semi-parametric model of SARS-CoV-2 transmission in the United States.

To support decision-making and policy for managing epidemics of emerging pathogens, we present a model for inference and scenario analysis of SARS-CoV-2 transmission in the USA. The stochastic SEIR-type model includes compartments for latent, asymptomatic, detected and undetected symptomatic individuals, and hospitalized cases, and features realistic interval distributions for presymptomatic and symptomatic periods, time varying rates of case detection, diagnosis, and mortality. The model accounts for the effects on transmission of human mobility using anonymized mobility data collected from cellular devices, and of difficult to quantify environmental and behavioral factors using a latent process. The baseline transmission rate is the product of a human mobility metric obtained from data and this fitted latent process. We fit the model to incident case and death reports for each state in the USA and Washington D.C., using likelihood Maximization by Iterated particle Filtering (MIF). Observations (daily case and death reports) are modeled as arising from a negative binomial reporting process. We estimate time-varying transmission rate, parameters of a sigmoidal time-varying fraction of hospitalized cases that result in death, extra-demographic process noise, two dispersion parameters of the observation process, and the initial sizes of the latent, asymptomatic, and symptomatic classes. In a retrospective analysis covering March-December 2020, we show how mobility and transmission strength became decoupled across two distinct phases of the pandemic. The decoupling demonstrates the need for flexible, semi-parametric approaches for modeling infectious disease dynamics in real-time.

Disasters collide at the intersection of extreme weather and infectious diseases.

Natural disasters interact to affect the resilience and prosperity of communities and disproportionately affect low income families and communities of colour. However, due to lack of a common theoretical framework, these are rarely quantified. Observing severe weather events (e.g. hurricanes and tornadoes) and epidemics (e.g. COVID-19) unfolding in southeastern US communities led us to conjecture that interactions among catastrophic disturbances might be much more considerable than previously recognized. For instance, hurricane evacuations increase human aggregation, a factor that affects the transmission of acute infections like SARS-CoV-2. Similarly, weather damage to health infrastructure can reduce a community's ability to provide services to people who are ill. As globalization and human population and movement continue to increase and weather events are becoming more intense, such complex interactions are expected to magnify and significantly impact environmental and human health.

Season and prey identity mediate the effect of predators on parasites in rodents: a test of the healthy herds hypothesis.

The healthy herds hypothesis (HHH) suggests that predators decrease parasitism in their prey. Repeated tests of this hypothesis across a range of taxa and ecosystems have revealed significant variation in the effect of predators on parasites in prey. Differences in the response to predators (1) between prey taxa, (2) between seasons, and (3) before and after catastrophic disturbance are common in natural systems, but typically ignored in empirical tests of the HHH. We used a predator exclusion experiment to measure the effect of these heterogeneities on the tri-trophic interaction among predators, parasites and prey. We experimentally excluded mammalian predators from the habitats of hispid cotton rats (Sigmodon hispidus) and cotton mice (Peromyscus gossypinus) and measured the effect of exclusion on gastrointestinal parasites in these rodents. Our experiment spanned multiple seasons and before and after a prescribed burn. We found that the exclusion of the same predators had opposite effects on the parasites of small mammal prey species. Additionally, we found that the effect of mammal exclusion on parasitism differed before versus after fire disturbance. Finally, we saw that the effect of predator exclusion was highly dependent on prey capture season. Significant effects of exclusion emerged primarily in the fall and winter months. The presence of so many different effects in one relatively simple system suggests that predator effects on parasites in prey are highly context dependent.

Do predators keep prey healthy or make them sicker? A meta-analysis.

Ecological theory suggests that predators can either keep prey populations healthy by reducing parasite burdens or alternatively, increase parasitism in prey. To quantify the overall magnitude and direction of the effect of predation on parasitism in prey observed in practice, we conducted a meta-analysis of 47 empirical studies. We also examined how study attributes, including parasite type and life cycle, habitat type, study design, and whether predators were able to directly consume prey contributed to variation in the predator-prey-parasite interaction. We found that the overall effect of predation on parasitism differed between parasites and parasitoids and that whether consumptive effects were present, and whether a predator was a non-host spreader of parasites, were the most important traits predicting the parasite response. Our results suggest that the mechanistic basis of predator-prey interactions strongly influences the effects of predators on parasites and that these effects, although context dependent, are predictable.

Overlapping timescales obscure early warning signals of the second COVID-19 wave.

Early warning indicators based on critical slowing down have been suggested as a model-independent and low-cost tool to anticipate the (re)emergence of infectious diseases. We studied whether such indicators could reliably have anticipated the second COVID-19 wave in European countries. Contrary to theoretical predictions, we found that characteristic early warning indicators generally decreased rather than increased prior to the second wave. A model explains this unexpected finding as a result of transient dynamics and the multiple timescales of relaxation during a non-stationary epidemic. Particularly, if an epidemic that seems initially contained after a first wave does not fully settle to its new quasi-equilibrium prior to changing circumstances or conditions that force a second wave, then indicators will show a decreasing rather than an increasing trend as a result of the persistent transient trajectory of the first wave. Our simulations show that this lack of timescale separation was to be expected during the second European epidemic wave of COVID-19. Overall, our results emphasize that the theory of critical slowing down applies only when the external forcing of the system across a critical point is slow relative to the internal system dynamics.

Five approaches to the suppression of SARS-CoV-2 without intensive social distancing.

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.

Socio-demographic, not environmental, risk factors explain fine-scale spatial patterns of diarrhoeal disease in Ifanadiana, rural Madagascar.

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.

Detecting critical slowing down in high-dimensional epidemiological systems.

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.

Assessing temperature-dependent competition between two invasive mosquito species.

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.

Protective population behavior change in outbreaks of emerging infectious disease.

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.

Association between tuberculosis in men and social network structure in Kampala, Uganda.

BACKGROUND: Globally, tuberculosis disease (TB) is more common among males than females. Recent research proposes that differences in social mixing by sex could alter infection patterns in TB. We examine evidence for two mechanisms by which social-mixing could increase men's contact rates with TB cases. First, men could be positioned in social networks such that they contact more people or social groups. Second, preferential mixing by sex could prime men to have more exposure to TB cases. METHODS: We compared the networks of male and female TB cases and healthy matched controls living in Kampala, Uganda. Specifically, we estimated their positions in social networks (network distance to TB cases, degree, betweenness, and closeness) and assortativity patterns (mixing with adult men, women, and children inside and outside the household). RESULTS: The observed network consisted of 11,840 individuals. There were few differences in estimates of node position by sex. We found distinct mixing patterns by sex and TB disease status including that TB cases have proportionally more adult male contacts and fewer contacts with children. CONCLUSIONS: This analysis used a network approach to study how social mixing patterns are associated with TB disease. Understanding these mechanisms may have implications for designing targeted intervention strategies in high-burden populations.

Do synthesis centers synthesize? A semantic analysis of topical diversity in research.

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.

Integrating data mining and transmission theory in the ecology of infectious diseases.

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.

Spatial and temporal variation in nest temperatures forecasts sex ratio skews in a crocodilian with environmental sex determination.

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.

The statistics of epidemic transitions.

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.

Early warning signals of malaria resurgence in Kericho, Kenya.

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.

Anticipating epidemic transitions with imperfect data.

Epidemic transitions are an important feature of infectious disease systems. As the transmissibility of a pathogen increases, the dynamics of disease spread shifts from limited stuttering chains of transmission to potentially large scale outbreaks. One proposed method to anticipate this transition are early-warning signals (EWS), summary statistics which undergo characteristic changes as the transition is approached. Although theoretically predicted, their mathematical basis does not take into account the nature of epidemiological data, which are typically aggregated into periodic case reports and subject to reporting error. The viability of EWS for epidemic transitions therefore remains uncertain. Here we demonstrate that most EWS can predict emergence even when calculated from imperfect data. We quantify performance using the area under the curve (AUC) statistic, a measure of how well an EWS distinguishes between numerical simulations of an emerging disease and one which is stationary. Values of the AUC statistic are compared across a range of different reporting scenarios. We find that different EWS respond to imperfect data differently. The mean, variance and first differenced variance all perform well unless reporting error is highly overdispersed. The autocorrelation, autocovariance and decay time perform well provided that the aggregation period of the data is larger than the serial interval and reporting error is not highly overdispersed. The coefficient of variation, skewness and kurtosis are found to be unreliable indicators of emergence. Overall, we find that seven of ten EWS considered perform well for most realistic reporting scenarios. We conclude that imperfect epidemiological data is not a barrier to using EWS for many potentially emerging diseases.

The global distribution and burden of dengue.

Dengue is a systemic viral infection transmitted between humans by Aedes mosquitoes. For some patients, dengue is a life-threatening illness. There are currently no licensed vaccines or specific therapeutics, and substantial vector control efforts have not stopped its rapid emergence and global spread. The contemporary worldwide distribution of the risk of dengue virus infection and its public health burden are poorly known. Here we undertake an exhaustive assembly of known records of dengue occurrence worldwide, and use a formal modelling framework to map the global distribution of dengue risk. We then pair the resulting risk map with detailed longitudinal information from dengue cohort studies and population surfaces to infer the public health burden of dengue in 2010. We predict dengue to be ubiquitous throughout the tropics, with local spatial variations in risk influenced strongly by rainfall, temperature and the degree of urbanization. Using cartographic approaches, we estimate there to be 390 million (95% credible interval 284-528) dengue infections per year, of which 96 million (67-136) manifest apparently (any level of disease severity). This infection total is more than three times the dengue burden estimate of the World Health Organization. Stratification of our estimates by country allows comparison with national dengue reporting, after taking into account the probability of an apparent infection being formally reported. The most notable differences are discussed. These new risk maps and infection estimates provide novel insights into the global, regional and national public health burden imposed by dengue. We anticipate that they will provide a starting point for a wider discussion about the global impact of this disease and will help to guide improvements in disease control strategies using vaccine, drug and vector control methods, and in their economic evaluation.