Complexity and anisotropy in host morphology make populations less susceptible to epidemic outbreaks

One of the challenges in epidemiology is to account for the complex morphological structure of hosts such as plant roots, crop fields, farms, cells, animal habitats and social networks, when the transmission of infection occurs between contiguous hosts. Morphological complexity brings an inherent heterogeneity in populations and affects the dynamics of pathogen spread in such systems. We have analysed the influence of realistically complex host morphology on the threshold for invasion and epidemic outbreak in an SIR (susceptible–infected–recovered) epidemiological model. We show that disorder expressed in the host morphology and anisotropy reduces the probability of epidemic outbreak and thus makes the system more resistant to epidemic outbreaks. We obtain general analytical estimates for minimally safe bounds for an invasion threshold and then illustrate their validity by considering an example of host data for branching hosts (salamander retinal ganglion cells). Several spatial arrangements of hosts with different degrees of heterogeneity have been considered in order to separately analyse the role of shape complexity and anisotropy in the host population. The estimates for invasion threshold are linked to morphological characteristics of the hosts that can be used for determining the threshold for invasion in practical applications.     

Estimating epidemic coupling between populations from the time to invasion

Identifying the mechanisms by which diseases spread among populations is important for understanding and forecasting patterns of epidemics and pandemics. Estimating transmission coupling among populations is challenging because transmission events are difficult to observe in practice, and connectivity among populations is often obscured by local disease dynamics. We consider the common situation in which an epidemic is seeded in one population and later spreads to a second population. We present a method for estimating transmission coupling between the two populations, assuming they can be modelled as susceptible–infected–removed (SIR) systems. We show that the strength of coupling between the two populations can be estimated from the time taken for the disease to invade the second population. Confidence in the estimate is low if only a single invasion event has been observed, but is substantially improved if numerous independent invasion events are observed. Our analysis of this simplest, idealized scenario represents a first step toward developing and verifying methods for estimating epidemic coupling among populations in an ever-more-connected global human population.     

Will an outbreak exceed available resources for control? Estimating the risk from invading pathogens using practical definitions of a severe epidemic

Forecasting whether or not initial reports of disease will be followed by a severe epidemic is an important component of disease management. Standard epidemic risk estimates involve assuming that infections occur according to a branching process and correspond to the probability that the outbreak persists beyond the initial stochastic phase. However, an alternative assessment is to predict whether or not initial cases will lead to a severe epidemic in which available control resources are exceeded. We show how this risk can be estimated by considering three practically relevant potential definitions of a severe epidemic; namely, an outbreak in which: (i) a large number of hosts are infected simultaneously; (ii) a large total number of infections occur; and (iii) the pathogen remains in the population for a long period. We show that the probability of a severe epidemic under these definitions often coincides with the standard branching process estimate for the major epidemic probability. However, these practically relevant risk assessments can also be different from the major epidemic probability, as well as from each other. This holds in different epidemiological systems, highlighting that careful consideration of how to classify a severe epidemic is vital for accurate epidemic risk quantification.     

Heterogeneity in susceptible–infected–removed (SIR) epidemics on lattices

The percolation paradigm is widely used in spatially explicit epidemic models where disease spreads between neighbouring hosts. It has been successful in identifying epidemic thresholds for invasion, separating non-invasive regimes, where the disease never invades the system, from invasive regimes where the probability of invasion is positive. However, its power is mainly limited to homogeneous systems. When heterogeneity (environmental stochasticity) is introduced, the value of the epidemic threshold is, in general, not predictable without numerical simulations. Here, we analyse the role of heterogeneity in a stochastic susceptible–infected–removed epidemic model on a two-dimensional lattice. In the homogeneous case, equivalent to bond percolation, the probability of invasion is controlled by a single parameter, the transmissibility of the pathogen between neighbouring hosts. In the heterogeneous model, the transmissibility becomes a random variable drawn from a probability distribution. We investigate how heterogeneity in transmissibility influences the value of the invasion threshold, and find that the resilience of the system to invasion can be suitably described by two control parameters, the mean and variance of the transmissibility. We analyse a two-dimensional phase diagram, where the threshold is represented by a phase boundary separating an invasive regime in the high-mean, low-variance region from a non-invasive regime in the low-mean, high-variance region of the parameter space. We thus show that the percolation paradigm can be extended to the heterogeneous case. Our results have practical implications for the analysis of disease control strategies in realistic heterogeneous epidemic systems.     

Simultaneous reconstruction of evolutionary history and epidemiological dynamics from viral sequences with the birth–death SIR model

The evolution of RNA viruses, such as human immunodeficiency virus (HIV), hepatitis C virus and influenza virus, occurs so rapidly that the viruses'' genomes contain information on past ecological dynamics. Hence, we develop a phylodynamic method that enables the joint estimation of epidemiological parameters and phylogenetic history. Based on a compartmental susceptible–infected–removed (SIR) model, this method provides separate information on incidence and prevalence of infections. Detailed information on the interaction of host population dynamics and evolutionary history can inform decisions on how to contain or entirely avoid disease outbreaks. We apply our birth–death SIR method to two viral datasets. First, five HIV type 1 clusters sampled in the UK between 1999 and 2003 are analysed. The estimated basic reproduction ratios range from 1.9 to 3.2 among the clusters. All clusters show a decline in the growth rate of the local epidemic in the middle or end of the 1990s. The analysis of a hepatitis C virus genotype 2c dataset shows that the local epidemic in the Córdoban city Cruz del Eje originated around 1906 (median), coinciding with an immigration wave from Europe to central Argentina that dates from 1880 to 1920. The estimated time of epidemic peak is around 1970.     

Prediction of invasion from the early stage of an epidemic

Predictability of undesired events is a question of great interest in many scientific disciplines including seismology, economy and epidemiology. Here, we focus on the predictability of invasion of a broad class of epidemics caused by diseases that lead to permanent immunity of infected hosts after recovery or death. We approach the problem from the perspective of the science of complexity by proposing and testing several strategies for the estimation of important characteristics of epidemics, such as the probability of invasion. Our results suggest that parsimonious approximate methodologies may lead to the most reliable and robust predictions. The proposed methodologies are first applied to analysis of experimentally observed epidemics: invasion of the fungal plant pathogen Rhizoctonia solani in replicated host microcosms. We then consider numerical experiments of the susceptible–infected–removed model to investigate the performance of the proposed methods in further detail. The suggested framework can be used as a valuable tool for quick assessment of epidemic threat at the stage when epidemics only start developing. Moreover, our work amplifies the significance of the small-scale and finite-time microcosm realizations of epidemics revealing their predictive power.     

On spatially explicit models of cholera epidemics

We generalize a recently proposed model for cholera epidemics that accounts for local communities of susceptibles and infectives in a spatially explicit arrangement of nodes linked by networks having different topologies. The vehicle of infection (Vibrio cholerae) is transported through the network links that are thought of as hydrological connections among susceptible communities. The mathematical tools used are borrowed from general schemes of reactive transport on river networks acting as the environmental matrix for the circulation and mixing of waterborne pathogens. Using the diffusion approximation, we analytically derive the speed of propagation for travelling fronts of epidemics on regular lattices (either one-dimensional or two-dimensional) endowed with uniform population density. Power laws are found that relate the propagation speed to the diffusion coefficient and the basic reproduction number. We numerically obtain the related, slower speed of epidemic spreading for more complex, yet realistic river structures such as Peano networks and optimal channel networks. The analysis of the limit case of uniformly distributed population sizes proves instrumental in establishing the overall conditions for the relevance of spatially explicit models. To that extent, the ratio between spreading and disease outbreak time scales proves the crucial parameter. The relevance of our results lies in the major differences potentially arising between the predictions of spatially explicit models and traditional compartmental models of the susceptible–infected–recovered (SIR)-like type. Our results suggest that in many cases of real-life epidemiological interest, time scales of disease dynamics may trigger outbreaks that significantly depart from the predictions of compartmental models.     

Competitive numerical analysis for stochastic HIV/AIDS epidemic model in a two‐sex population

This study is an attempt to explain a reliable numerical analysis of a stochastic HIV/AIDS model in a two‐sex population considering counselling and antiretroviral therapy (ART). The authors are comparing the solutions of the stochastic and deterministic HIV/AIDS epidemic model. Here, an endeavour has been made to explain the stochastic HIV/AIDS epidemic model is comparatively more pragmatic in contrast with the deterministic HIV/AIDS epidemic model. The effect of threshold number H ^* holds on the stochastic HIV/AIDS epidemic model. If H ^*  < 1 then condition helps us to control disease in a two‐sex human population while H ^*  > 1 explains the persistence of disease in the two‐sex human population. Lamentably, numerical methods such as Euler–Maruyama, stochastic Euler, and stochastic Runge–Kutta do not work for large time step sizes. The recommended structure preserving framework of the stochastic non‐standard finite difference (SNSFD) scheme conserve all vital characteristics such as positivity, boundedness, and dynamical consistency defined by Mickens. The effectiveness of counselling and ART may control HIV/AIDS in a two‐sex population.Inspec keywords: diseases, stochastic processes, epidemics, patient treatment, finite difference methodsOther keywords: two‐sex human population, antiretroviral therapy, competitive numerical analysis, stochastic HIV‐AIDS epidemic model, structure preserving framework, stochastic nonstandard finite difference scheme, SNSFD scheme, deterministic HIV‐AIDS epidemic model     

The evolutionary dynamics of within-generation immune priming in invertebrate hosts

While invertebrates lack the machinery necessary for ‘acquired immunity’, there is increasing empirical evidence that exposure to low levels of disease may ‘prime’ an invertebrate''s immune response, increasing its defence to subsequent exposure. Despite this increasing empirical data, there has been little theoretical attention paid to immune priming. Here, we investigate the evolution of immune priming, focusing on the role of the unique feedbacks generated by a newly developed susceptible–primed–infected epidemiological model. Contrasting our results with previous models on the evolution of acquired immunity, we highlight that there are important implications to the evolution of immunity through priming owing to these different epidemiological feedbacks. In particular, we find that in contrast to acquired immunity, priming is strongly selected for at high as well as intermediate pathogen virulence. We also find that priming may be greatest at either intermediate or high host lifespans depending on the severity of disease. Furthermore, hosts faced with more severe pathogens are more likely to evolve diversity in priming. Finally, we show when the evolution of priming leads to the exclusion of the pathogens or hosts experiencing population cycles. Overall the model acts as a baseline for understanding the evolution of priming in host–pathogen systems.     

Unrelated toxin–antitoxin systems cooperate to induce persistence

Persisters are drug-tolerant bacteria that account for the majority of bacterial infections. They are not mutants, rather, they are slow-growing cells in an otherwise normally growing population. It is known that the frequency of persisters in a population is correlated with the number of toxin–antitoxin systems in the organism. Our previous work provided a mechanistic link between the two by showing how multiple toxin–antitoxin systems, which are present in nearly all bacteria, can cooperate to induce bistable toxin concentrations that result in a heterogeneous population of slow- and fast-growing cells. As such, the slow-growing persisters are a bet-hedging subpopulation maintained under normal conditions. For technical reasons, the model assumed that the kinetic parameters of the various toxin–antitoxin systems in the cell are identical, but experimental data indicate that they differ, sometimes dramatically. Thus, a critical question remains: whether toxin–antitoxin systems from the diverse families, often found together in a cell, with significantly different kinetics, can cooperate in a similar manner. Here, we characterize the interaction of toxin–antitoxin systems from many families that are unrelated and kinetically diverse, and identify the essential determinant for their cooperation. The generic architecture of toxin–antitoxin systems provides the potential for bistability, and our results show that even when they do not exhibit bistability alone, unrelated systems can be coupled by the growth rate to create a strongly bistable, hysteretic switch between normal (fast-growing) and persistent (slow-growing) states. Different combinations of kinetic parameters can produce similar toxic switching thresholds, and the proximity of the thresholds is the primary determinant of bistability. Stochastic fluctuations can spontaneously switch all of the toxin–antitoxin systems in a cell at once. The spontaneous switch creates a heterogeneous population of growing and non-growing cells, typical of persisters, that exist under normal conditions, rather than only as an induced response. The frequency of persisters in the population can be tuned for a particular environmental niche by mixing and matching unrelated systems via mutation, horizontal gene transfer and selection.     

Age-dependent partnering and the HIV transmission chain: a microsimulation analysis

Efficient planning and evaluation of human immunodeficiency virus (HIV) prevention programmes requires an understanding of what sustains the epidemic, including the mechanism by which HIV transmission keeps pace with the ageing of the infected population. Recently, more detailed population models have been developed which represent the epidemic with sufficient detail to characterize the dynamics of ongoing transmission. Here, we describe the structure and parameters of such a model, called EMOD-HIV v. 0.7. We analyse the chains of transmission that allow the HIV epidemic to propagate across age groups in this model. In order to prevent the epidemic from dying out, the virus must find younger victims faster than its extant victims age and die. The individuals who enable such transmission events in EMOD-HIV v. 0.7 are higher concurrency, co-infected males aged 26–29 and females aged 23–24. Prevention programmes that target these populations could efficiently interrupt the mechanisms that allow HIV to transmit at a pace that is faster than the progress of time.     

Sampling through time and phylodynamic inference with coalescent and birth–death models

Many population genetic models have been developed for the purpose of inferring population size and growth rates from random samples of genetic data. We examine two popular approaches to this problem, the coalescent and the birth–death-sampling model (BDM), in the context of estimating population size and birth rates in a population growing exponentially according to the birth–death branching process. For sequences sampled at a single time, we found the coalescent and the BDM gave virtually indistinguishable results in terms of the growth rates and fraction of the population sampled, even when sampling from a small population. For sequences sampled at multiple time points, we find that the birth–death model estimators are subject to large bias if the sampling process is misspecified. Since BDMs incorporate a model of the sampling process, we show how much of the statistical power of BDMs arises from the sequence of sample times and not from the genealogical tree. This motivates the development of a new coalescent estimator, which is augmented with a model of the known sampling process and is potentially more precise than the coalescent that does not use sample time information.     

Bayesian data assimilation provides rapid decision support for vector-borne diseases

Predicting the spread of vector-borne diseases in response to incursions requires knowledge of both host and vector demographics in advance of an outbreak. Although host population data are typically available, for novel disease introductions there is a high chance of the pathogen using a vector for which data are unavailable. This presents a barrier to estimating the parameters of dynamical models representing host–vector–pathogen interaction, and hence limits their ability to provide quantitative risk forecasts. The Theileria orientalis (Ikeda) outbreak in New Zealand cattle demonstrates this problem: even though the vector has received extensive laboratory study, a high degree of uncertainty persists over its national demographic distribution. Addressing this, we develop a Bayesian data assimilation approach whereby indirect observations of vector activity inform a seasonal spatio-temporal risk surface within a stochastic epidemic model. We provide quantitative predictions for the future spread of the epidemic, quantifying uncertainty in the model parameters, case infection times and the disease status of undetected infections. Importantly, we demonstrate how our model learns sequentially as the epidemic unfolds and provide evidence for changing epidemic dynamics through time. Our approach therefore provides a significant advance in rapid decision support for novel vector-borne disease outbreaks.     

A century of transitions in New York City's measles dynamics

Infectious diseases spreading in a human population occasionally exhibit sudden transitions in their qualitative dynamics. Previous work has successfully predicted such transitions in New York City''s historical measles incidence using the seasonally forced susceptible–infectious–recovered (SIR) model. This work relied on a dataset spanning 45 years (1928–1973), which we have extended to 93 years (1891–1984). We identify additional dynamical transitions in the longer dataset and successfully explain them by analysing attractors and transients of the same mechanistic epidemiological model.     

Disease control across urban–rural gradients

Controlling the regional re-emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) after its initial spread in ever-changing personal contact networks and disease landscapes is a challenging task. In a landscape context, contact opportunities within and between populations are changing rapidly as lockdown measures are relaxed and a number of social activities re-activated. Using an individual-based metapopulation model, we explored the efficacy of different control strategies across an urban–rural gradient in Wales, UK. Our model shows that isolation of symptomatic cases or regional lockdowns in response to local outbreaks have limited efficacy unless the overall transmission rate is kept persistently low. Additional isolation of non-symptomatic infected individuals, who may be detected by effective test-and-trace strategies, is pivotal to reducing the overall epidemic size over a wider range of transmission scenarios. We define an ‘urban–rural gradient in epidemic size'' as a correlation between regional epidemic size and connectivity within the region, with more highly connected urban populations experiencing relatively larger outbreaks. For interventions focused on regional lockdowns, the strength of such gradients in epidemic size increased with higher travel frequencies, indicating a reduced efficacy of the control measure in the urban regions under these conditions. When both non-symptomatic and symptomatic individuals are isolated or regional lockdown strategies are enforced, we further found the strongest urban–rural epidemic gradients at high transmission rates. This effect was reversed for strategies targeted at symptomatic individuals only. Our results emphasize the importance of test-and-trace strategies and maintaining low transmission rates for efficiently controlling SARS-CoV-2 spread, both at landscape scale and in urban areas.     

Synchronization-induced persistence versus selection for habitats in spatially coupled ecosystems

Critical population phase transitions, in which a persistent population becomes extinction-prone owing to environmental changes, are fundamentally important in ecology, and their determination is a key factor in successful ecosystem management. To persist, a species requires a suitable environment in a sufficiently large spatial region. However, even if this condition is met, the species does not necessarily persist, owing to stochastic fluctuations. Here, we develop a model that allows simultaneous investigation of extinction due to either stochastic or deterministic reasons. We find that even classic birth–death processes in spatially extended ecosystems exhibit phase transitions between extinction-prone and persistent populations. Sometimes these are first-order transitions, which means that environmental changes may result in irreversible population collapse. Moreover, we find that higher migration rates not only lead to higher robustness to stochastic fluctuations, but also result in lower sustainability in heterogeneous environments by preventing efficient selection for suitable habitats. This demonstrates that intermediate migration rates are optimal for survival. At low migration rates, the dynamics are reduced to metapopulation dynamics, whereas at high migration rates, the dynamics are reduced to a multi-type branching process. We focus on species persistence, but our results suggest a unique method for finding phase transitions in spatially extended stochastic systems in general.     

Predictability in a highly stochastic system: final size of measles epidemics in small populations

A standard assumption in the modelling of epidemic dynamics is that the population of interest is well mixed, and that no clusters of metapopulations exist. The well-known and oft-used SIR model, arguably the most important compartmental model in theoretical epidemiology, assumes that the disease being modelled is strongly immunizing, directly transmitted and has a well-defined period of infection, in addition to these population mixing assumptions. Childhood infections, such as measles, are prime examples of diseases that fit the SIR-like mechanism. These infections have been well studied for many systems with large, well-mixed populations with endemic infection. Here, we consider a setting where populations are small and isolated. The dynamics of infection are driven by stochastic extinction–recolonization events, producing large, sudden and short-lived epidemics before rapidly dying out from a lack of susceptible hosts. Using a TSIR model, we fit prevaccination measles incidence and demographic data in Bornholm, the Faroe Islands and four districts of Iceland, between 1901 and 1965. The datasets for each of these countries suffer from different levels of data heterogeneity and sparsity. We explore the potential for prediction of this model: given historical incidence data and up-to-date demographic information, and knowing that a new epidemic has just begun, can we predict how large it will be? We show that, despite a lack of significant seasonality in the incidence of measles cases, and potentially severe heterogeneity at the population level, we are able to estimate the size of upcoming epidemics, conditioned on the first time step, to within reasonable confidence. Our results have potential implications for possible control measures for the early stages of new epidemics in small populations.     

Disease transmission in territorial populations: the small-world network of Serengeti lions

Territoriality in animal populations creates spatial structure that is thought to naturally buffer disease invasion. Often, however, territorial populations also include highly mobile, non-residential individuals that potentially serve as disease superspreaders. Using long-term data from the Serengeti Lion Project, we characterize the contact network structure of a territorial wildlife population and address the epidemiological impact of nomadic individuals. As expected, pride contacts are dominated by interactions with neighbouring prides and interspersed by encounters with nomads as they wander throughout the ecosystem. Yet the pride–pride network also includes occasional long-range contacts between prides, making it surprisingly small world and vulnerable to epidemics, even without nomads. While nomads increase both the local and global connectivity of the network, their epidemiological impact is marginal, particularly for diseases with short infectious periods like canine distemper virus. Thus, territoriality in Serengeti lions may be less protective and non-residents less important for disease transmission than previously considered.     

Multi‐scale approach for simulating time‐delay biochemical reaction systems

This study presents a multi‐scale approach for simulating time‐delay biochemical reaction systems when there are wide ranges of molecular numbers. The authors construct a new efficient approach based on partitioning into slow and fast subsets in conjunction with predictor–corrector methods. This multi‐scale approach is shown to be much more efficient than existing methods such as the delay stochastic simulation algorithm and the modified next reaction method. Numerical testing on several important problems in systems biology confirms the accuracy and computational efficiency of this approach.Inspec keywords: biochemistry, delays, biological techniques, predictor‐corrector methodsOther keywords: multiscale approach, time‐delay biochemical reaction systems, predictor–corrector methods, delay stochastic simulation algorithm, modified next reaction method, numerical testing, systems biology, method accuracy, computational efficiency     

Modular bond‐graph modelling and analysis of biomolecular systems

Bond graphs can be used to build thermodynamically‐compliant hierarchical models of biomolecular systems. As bond graphs have been widely used to model, analyse and synthesise engineering systems, this study suggests that they can play the same rôle in the modelling, analysis and synthesis of biomolecular systems. The particular structure of bond graphs arising from biomolecular systems is established and used to elucidate the relation between thermodynamically closed and open systems. Block diagram representations of the dynamics implied by these bond graphs are used to reveal implicit feedback structures and are linearised to allow the application of control‐theoretical methods. Two concepts of modularity are examined: computational modularity where physical correctness is retained and behavioural modularity where module behaviour (such as ultrasensitivity) is retained. As well as providing computational modularity, bond graphs provide a natural formulation of behavioural modularity and reveal the sources of retroactivity. A bond graph approach to reducing retroactivity, and thus inter‐module interaction, is shown to require a power supply such as that provided by the ATP ⇌ ADP + Pi reaction. The mitogen‐activated protein kinase cascade (Raf–MEK–ERK pathway) is used as an illustrative example.Inspec keywords: molecular biophysics, bond graphs, hierarchical systems, thermodynamics, enzymes, physiological models, biology computingOther keywords: signalling networks, behavioural modularity, Michaelis‐Menten kinetics, Raf‐MEK‐ERK pathway, mitogen‐activated protein kinase cascade, ATP⇌ADP + Pi reaction, intermodule interaction, retroactivity, computational modularity, block diagram representations, thermodynamically‐compliant hierarchical models, biomolecular systems, modular bond‐graph modelling