Search for: All records

Abstract The evolutionary fate of multi-strain pathogens is shaped by host-pathogen ecological interactions. In bacterial pathogens of plants, enhanced strain characterization and advances in our understanding of molecular mechanisms underlying defense pathways open the door for revisiting the role of negative frequency-dependent selection (NFDS) in strain structure, including its interplay with genetic exchange. NFDS arising from specific defense is one potential mechanism for generating, maintaining, and structuring pathogen diversity. In plants, specific protection against microbial pathogens involves Resistance proteins (R-proteins) that recognize virulence factors (effectors) secreted by pathogens, typically to subvert the initial line of host defense. Here we formulate a stochastic computational co-evolution model that explicitly incorporates variable length R-gene and effector repertoires, and migration from their regional pools. We use this model to understand potential mechanisms shaping effector repertoire structure and associated strain coexistence in the generalist plant pathogenP. syringae. The demonstration of a modular structure in our numerical simulations motivates the analysis of genome sequences from 76 strains collected in the Midwestern US and 1104 strains from global sources. We find that effector repertories both locally and globally exhibit a modular structure, with higher similarity within than between clusters. The observed modules are consistent with the core genome phylogeny and are unexplained by plant host species, location of isolation, and genetic linkage between effectors. An extension of the model is needed to take into account the evidence for genetic exchange and the phylogenetic congruence of effector modules. We initialize the system with a phylogenetically congruent modular structure and include recombination rates decreasing as a function of phylogenetic distance. We show that NFDS can counter-balance the effects of mixing due to recombination and in so doing, contributes to the maintenance of strain structure. These findings indicate that the observed similarity clusters may constitute, in part, emergent niches arising from eco-evolutionary dynamics that contribute to strain coexistence. 

The coevolutionary dynamics of lytic viruses and microbes with CRISPR-Cas immunity exhibit alternations between sustained host control of viral proliferation and major viral epidemics in previous computational models. Thesealternatingdynamics have yet to be observed in other host–pathogen systems. Here, we address the breakdown of control and transition to large outbreaks with a stochastic eco-evolutionary model. We establish the role of host density-dependent competition in punctuated virus-driven succession and associated diversity trends that concentrate escape pathways during control phases. Using infection and escape networks, we derive the viral emergence probability whose fluctuations of increasing size and frequency characterize the approach to large outbreaks. We explore alternation probabilities as a function of non-dimensional parameters related to the probability of viral escape and host competition. Our results demonstrate how emergent feedbacks between host competition and viral diversification render the host immune structure fragile, potentiating a dynamical transition to large epidemics. 

Mathematical models in ecology and epidemiology must be consistent with observed data in order to generate reliable knowledge and evidence-based policy. Metapopulation systems, which consist of a network of connected sub-populations, pose technical challenges in statistical inference owing to nonlinear, stochastic interactions. Numerical difficulties encountered in conducting inference can obstruct the core scientific questions concerning the link between the mathematical models and the data. Recently, an algorithm has been proposed that enables computationally tractable likelihood-based inference for high-dimensional partially observed stochastic dynamic models of metapopulation systems. We use this algorithm to build a statistically principled data analysis workflow for metapopulation systems. Via a case study of COVID-19, we show how this workflow addresses the limitations of previous approaches. The COVID-19 pandemic provides a situation where mathematical models and their policy implications are widely visible, and we revisit an influential metapopulation model used to inform basic epidemiological understanding early in the pandemic. Our methods support self-critical data analysis, enabling us to identify and address model weaknesses, leading to a new model with substantially improved statistical fit and parameter identifiability. Our results suggest that the lockdown initiated on 23 January 2020 in China was more effective than previously thought. 

The role of climate factors on transmission of mosquito-borne infections within urban landscapes must be considered in the context of the pronounced spatial heterogeneity of such environments. Socio-demographic and environmental variation challenge control efforts for emergent arboviruses transmitted via the urban mosquitoAedes aegypti. We address at high resolution, the spatial heterogeneity of dengue transmission risk in the megacity of Delhi, India, as a function of both temperature and the carrying-capacity of the human environment for the mosquito. Based on previous results predicting maximum mosquitoes per human for different socio-economic typologies, and on remote sensing temperature data, we produce a map of the reproductive number of dengue at a resolution of 250m by 250m. We focus on dengue risk hotspots during inter-epidemic periods, places where chains of transmission can persist for longer. We assess the resulting high-resolution risk map of dengue with reported cases for three consecutive boreal winters. We find that both temperature and vector carrying-capacity per human co-vary in space because of their respective dependence on population density. The synergistic action of these two factors results in larger variation of dengue’s reproductive number than when considered separately, with poor and dense locations experiencing the warmest conditions and becoming the most likely reservoirs off-season. The location of observed winter cases is accurately predicted for different risk threshold criteria. Results underscore the inequity of risk across a complex urban landscape, whereby individuals in dense poor neighborhoods face the compounded effect of higher temperatures and mosquito carrying capacity. Targeting chains of transmission in inter-epidemic periods at these locations should be a priority of control efforts. A better mapping is needed of the interplay between climate factors that are dominant determinants of the seasonality of vector-borne infections and the socio-economic conditions behind unequal exposure. 

Identifying climate drivers is essential to understand and predict epidemics of mosquito-borne infections whose population dynamics typically exhibit seasonality and multiannual cycles. Which climate covariates to consider varies across studies, from local factors such as temperature to remote drivers such as the El Niño–Southern Oscillation. With partial wavelet coherence, we present a systematic investigation of nonstationary associations between mosquito-borne disease incidence and a given climate factor while controlling for another. Analysis of almost 200 time series of dengue and malaria around the globe at different geographical scales shows a systematic effect of global climate drivers on interannual variability and of local ones on seasonality. This clear separation of time scales of action enhances detection of climate drivers and indicates those best suited for building early-warning systems. 

ABSTRACT Microbial host populations evolve traits conferring specific resistance to viral predators via various defence mechanisms, while viruses reciprocally evolve traits to evade these defences. Such coevolutionary dynamics often involve diversification promoted by negative frequency‐dependent selection. However, microbial traits conferring competitive asymmetries can induce directional selection, opposing diversification. Despite extensive research on microbe–virus coevolution, the combined effect of both host trait types and associated selection remains unclear. Using a CRISPR‐mediated coevolutionary system, we examine how the co‐occurrence of both trait types impacts viral evolution and persistence, previously shown to be transient and nonstationary in computational models. A stochastic model incorporating host competitive asymmetries via variation of intrinsic growth rates reveals that competitively advantaged host clades generate the majority of immune diversity. Greater asymmetries extend viral extinction times, accelerate viral adaptation locally in time and augment long‐term local adaptation. These findings align with previous experiments and provide further insights into long‐term coevolutionary dynamics. 

Abstract The spread of dengue and other arboviruses constitutes an expanding global health threat. The extensive heterogeneity in population distribution and potential complexity of movement in megacities of low and middle-income countries challenges predictive modeling, even as its importance to disease spread is clearer than ever. Using surveillance data at fine resolution following the emergence of the DENV4 dengue serotype in Rio de Janeiro, we document a pattern in the size of successive epidemics that is invariant to the scale of spatial aggregation. This pattern emerges from the combined effect of herd immunity and seasonal transmission, and is strongly driven by variation in population density at sub-kilometer scales. It is apparent only when the landscape is stratified by population density and not by spatial proximity as has been common practice. Models that exploit this emergent simplicity should afford improved predictions of the local size of successive epidemic waves. 

The contributions of asymptomatic infections to herd immunity and community transmission are key to the resurgence and control of COVID-19, but are difficult to estimate using current models that ignore changes in testing capacity. Using a model that incorporates daily testing information fit to the case and serology data from New York City, we show that the proportion of symptomatic cases is low, ranging from 13 to 18%, and that the reproductive number may be larger than often assumed. Asymptomatic infections contribute substantially to herd immunity, and to community transmission together with presymptomatic ones. If asymptomatic infections transmit at similar rates as symptomatic ones, the overall reproductive number across all classes is larger than often assumed, with estimates ranging from 3.2 to 4.4. If they transmit poorly, then symptomatic cases have a larger reproductive number ranging from 3.9 to 8.1. Even in this regime, presymptomatic and asymptomatic cases together comprise at least 50% of the force of infection at the outbreak peak. We find no regimes in which all infection subpopulations have reproductive numbers lower than three. These findings elucidate the uncertainty that current case and serology data cannot resolve, despite consideration of different model structures. They also emphasize how temporal data on testing can reduce and better define this uncertainty, as we move forward through longer surveillance and second epidemic waves. Complementary information is required to determine the transmissibility of asymptomatic cases, which we discuss. Regardless, current assumptions about the basic reproductive number of severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) should be reconsidered. 

Predicting arbovirus re-emergence remains challenging in regions with limited off-season transmission and intermittent epidemics. Current mathematical models treat the depletion and replenishment of susceptible (non-immune) hosts as the principal drivers of re-emergence, based on established understanding of highly transmissible childhood diseases with frequent epidemics. We extend an analytical approach to determine the number of ‘skip’ years preceding re-emergence for diseases with continuous seasonal transmission, population growth and under-reporting. Re-emergence times are shown to be highly sensitive to small changes in low R 0 (secondary cases produced from a primary infection in a fully susceptible population). We then fit a stochastic Susceptible–Infected–Recovered (SIR) model to observed case data for the emergence of dengue serotype DENV1 in Rio de Janeiro. This aggregated city-level model substantially over-estimates observed re-emergence times either in terms of skips or outbreak probability under forward simulation. The inability of susceptible depletion and replenishment to explain re-emergence under ‘well-mixed’ conditions at a city-wide scale demonstrates a key limitation of SIR aggregated models, including those applied to other arboviruses. The predictive uncertainty and high skip sensitivity to epidemiological parameters suggest a need to investigate the relevant spatial scales of susceptible depletion and the scaling of microscale transmission dynamics to formulate simpler models that apply at coarse resolutions.