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.     

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.     

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.     

Epidemics in networks of spatially correlated three-dimensional root-branching structures

Using digitized images of the three-dimensional, branching structures for root systems of bean seedlings, together with analytical and numerical methods that map a common susceptible–infected–recovered (‘SIR’) epidemiological model onto the bond percolation problem, we show how the spatially correlated branching structures of plant roots affect transmission efficiencies, and hence the invasion criterion, for a soil-borne pathogen as it spreads through ensembles of morphologically complex hosts. We conclude that the inherent heterogeneities in transmissibilities arising from correlations in the degrees of overlap between neighbouring plants render a population of root systems less susceptible to epidemic invasion than a corresponding homogeneous system. Several components of morphological complexity are analysed that contribute to disorder and heterogeneities in the transmissibility of infection. Anisotropy in root shape is shown to increase resilience to epidemic invasion, while increasing the degree of branching enhances the spread of epidemics in the population of roots. Some extension of the methods for other epidemiological systems are discussed.     

Demographic buffering: titrating the effects of birth rate and imperfect immunity on epidemic dynamics

Host demography can alter the dynamics of infectious disease. In the case of perfectly immunizing infections, observations of strong sensitivity to demographic variation have been mechanistically explained through analysis of the susceptible–infected–recovered (SIR) model that assumes lifelong immunity following recovery from infection. When imperfect immunity is incorporated into this framework via the susceptible–infected–recovered–susceptible (SIRS) model, with individuals regaining full susceptibility following recovery, we show that rapid loss of immunity is predicted to buffer populations against the effects of demographic change. However, this buffering is contrary to the dependence on demography recently observed for partially immunizing infections such as rotavirus and respiratory syncytial virus. We show that this discrepancy arises from a key simplification embedded in the SIR(S) framework, namely that the potential for differential immune responses to repeat exposures is ignored. We explore the minimum additional immunological information that must be included to reflect the range of observed dependencies on demography. We show that including partial protection and lower transmission following primary infection is sufficient to capture more realistic reduced levels of buffering, in addition to changes in epidemic timing, across a range of partially and fully immunizing infections. Furthermore, our results identify key variables in this relationship, including R₀.     

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.     

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.     

Transmission network of the 2014–2015 Ebola epidemic in Sierra Leone

Understanding the growth and spatial expansion of (re)emerging infectious disease outbreaks, such as Ebola and avian influenza, is critical for the effective planning of control measures; however, such efforts are often compromised by data insufficiencies and observational errors. Here, we develop a spatial–temporal inference methodology using a modified network model in conjunction with the ensemble adjustment Kalman filter, a Bayesian inference method equipped to handle observational errors. The combined method is capable of revealing the spatial–temporal progression of infectious disease, while requiring only limited, readily compiled data. We use this method to reconstruct the transmission network of the 2014–2015 Ebola epidemic in Sierra Leone and identify source and sink regions. Our inference suggests that, in Sierra Leone, transmission within the network introduced Ebola to neighbouring districts and initiated self-sustaining local epidemics; two of the more populous and connected districts, Kenema and Port Loko, facilitated two independent transmission pathways. Epidemic intensity differed by district, was highly correlated with population size (r = 0.76, p = 0.0015) and a critical window of opportunity for containing local Ebola epidemics at the source (ca one month) existed. This novel methodology can be used to help identify and contain the spatial expansion of future (re)emerging infectious disease outbreaks.     

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.     

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     

Utility of R0 as a predictor of disease invasion in structured populations

Early theoretical work on disease invasion typically assumed large and well-mixed host populations. Many human and wildlife systems, however, have small groups with limited movement among groups. In these situations, the basic reproductive number, R0, is likely to be a poor predictor of a disease pandemic because it typically does not account for group structure and movement of individuals among groups. We extend recent work by combining the movement of hosts, transmission within groups, recovery from infection and the recruitment of new susceptibles into a stochastic model of disease in a host metapopulation. We focus on how recruitment of susceptibles affects disease invasion and how population structure can affect the frequency of superspreading events (SSEs). We show that the frequency of SSEs may decrease with the reduced movement and the group sizes due to the limited number of susceptible individuals available. Classification tree analysis of the model results illustrates the hierarchical nature of disease invasion in host metapopulations. First, the pathogen must effectively transmit within a group (R0>1), and then the pathogen must persist within a group long enough to allow for movement among the groups. Therefore, the factors affecting disease persistence--such as infectious period, group size and recruitment of new susceptibles--are as important as the local transmission rates in predicting the spread of pathogens across a metapopulation.     

Joint modelling of serological and hospitalization data reveals that high levels of pre-existing immunity and school holidays shaped the influenza A pandemic of 2009 in The Netherlands

Obtaining a quantitative understanding of the transmission dynamics of influenza A is important for predicting healthcare demand and assessing the likely impact of intervention measures. The pandemic of 2009 provides an ideal platform for developing integrative analyses as it has been studied intensively, and a wealth of data sources is available. Here, we analyse two complementary datasets in a disease transmission framework: cross-sectional serological surveys providing data on infection attack rates, and hospitalization data that convey information on the timing and duration of the pandemic. We estimate key epidemic determinants such as infection and hospitalization rates, and the impact of a school holiday. In contrast to previous approaches, our novel modelling of serological data with mixture distributions provides a probabilistic classification of individual samples (susceptible, immune and infected), propagating classification uncertainties to the transmission model and enabling serological classifications to be informed by hospitalization data. The analyses show that high levels of immunity among persons 20 years and older provide a consistent explanation of the skewed attack rates observed during the pandemic and yield precise estimates of the probability of hospitalization per infection (1–4 years: 0.00096 (95%CrI: 0.00078–0.0012); 5–19 years: 0.00036 (0.00031–0.0044); 20–64 years: 0.0015 (0.00091–0.0020); 65+ years: 0.0084 (0.0028–0.016)). The analyses suggest that in The Netherlands, the school holiday period reduced the number of infectious contacts between 5- and 9-year-old children substantially (estimated reduction: 54%; 95%CrI: 29–82%), thereby delaying the unfolding of the pandemic in The Netherlands by approximately a week.     

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.     

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.     

Epidemic cycles driven by host behaviour

Host immunity and demographics (the recruitment of susceptibles via birthrate) have been demonstrated to be a key determinant of the periodicity of measles, pertussis and dengue epidemics. However, not all epidemic cycles are from pathogens inducing sterilizing immunity or are driven by demographics. Many sexually transmitted infections are driven by sexual behaviour. We present a mathematical model of disease transmission where individuals can disconnect and reconnect depending on the infectious status of their contacts. We fit the model to historic syphilis (Treponema pallidum) and gonorrhea (Neisseria gonorrhoeae) incidence in the USA and explore potential intervention strategies against syphilis. We find that cycles in syphilis incidence can be driven solely by changing sexual behaviour in structured populations. Our model also explains the lack of similar cycles in gonorrhea incidence even if the two infections share the same propagation pathways. Our model similarly illustrates how sudden epidemic outbreaks can occur on time scales smaller than the characteristic demographic time scale of the population and that weaker infections can lead to more violent outbreaks. Behaviour also appears to be critical for control strategies as we found a bigger sensitivity to behavioural interventions than antibiotic treatment. Thus, behavioural interventions may play a larger role than previously thought, especially in the face of antibiotic resistance and low intervention efficacies.     

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.     

Modelling cholera epidemics: the role of waterways,human mobility and sanitation

We investigate the role of human mobility as a driver for long-range spreading of cholera infections, which primarily propagate through hydrologically controlled ecological corridors. Our aim is to build a spatially explicit model of a disease epidemic, which is relevant to both social and scientific issues. We present a two-layer network model that accounts for the interplay between epidemiological dynamics, hydrological transport and long-distance dissemination of the pathogen Vibrio cholerae owing to host movement, described here by means of a gravity-model approach. We test our model against epidemiological data recorded during the extensive cholera outbreak occurred in the KwaZulu-Natal province of South Africa during 2000–2001. We show that long-range human movement is fundamental in quantifying otherwise unexplained inter-catchment transport of V. cholerae, thus playing a key role in the formation of regional patterns of cholera epidemics. We also show quantitatively how heterogeneously distributed drinking water supplies and sanitation conditions may affect large-scale cholera transmission, and analyse the effects of different sanitation policies.     

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.     

Estimation of measles vaccine efficacy and critical vaccination coverage in a highly vaccinated population

Measles is a highly infectious disease that has been targeted for elimination from four WHO regions. Whether and under which conditions this goal is feasible is, however, uncertain since outbreaks have been documented in populations with high vaccination coverage (more than 90%). Here, we use the example of a large outbreak in a German public school to show how estimates of key epidemiological parameters such as the basic reproduction number (R₀), vaccine efficacy (VE_S) and critical vaccination coverage (p_c) can be obtained from partially observed outbreaks in highly vaccinated populations. Our analyses rely on Bayesian methods of inference based on the final size distribution of outbreak size, and use data which are easily collected. For the German public school the analyses indicate that the basic reproduction number of measles is higher than previously thought (, 95% credible interval: 23.6–40.4), that the vaccine is highly effective in preventing infection (, 95% credible interval: 0.993–0.999), and that a vaccination coverage in excess of 95 per cent may be necessary to achieve herd immunity (, 95% credible interval: 0.961–0.978). We discuss the implications for measles elimination from highly vaccinated populations.