Influence of network dynamics on the spread of sexually transmitted diseases

Network epidemiology often assumes that the relationships defining the social network of a population are static. The dynamics of relationships is only taken indirectly into account by assuming that the relevant information to study epidemic spread is encoded in the network obtained, by considering numbers of partners accumulated over periods of time roughly proportional to the infectious period of the disease. On the other hand, models explicitly including social dynamics are often too schematic to provide a reasonable representation of a real population, or so detailed that no general conclusions can be drawn from them. Here, we present a model of social dynamics that is general enough so its parameters can be obtained by fitting data from surveys about sexual behaviour, but that can still be studied analytically, using mean-field techniques. This allows us to obtain some general results about epidemic spreading. We show that using accumulated network data to estimate the static epidemic threshold lead to a significant underestimation of that threshold. We also show that, for a dynamic network, the relative epidemic threshold is an increasing function of the infectious period of the disease, implying that the static value is a lower bound to the real threshold. A practical example is given of how to apply the model to the study of a real population.     

Exogenous re-infection and the dynamics of tuberculosis epidemics: local effects in a network model of transmission

Infection with Mycobacterium tuberculosis leads to tuberculosis (TB) disease by one of the three possible routes: primary progression after a recent infection; re-activation of a latent infection; or exogenous re-infection of a previously infected individual. Recent studies show that optimal TB control strategies may vary depending on the predominant route to disease in a specific population. It is therefore important for public health policy makers to understand the relative frequency of each type of TB within specific epidemiological scenarios. Although molecular epidemiologic tools have been used to estimate the relative contribution of recent transmission and re-activation to the burden of TB disease, it is not possible to use these techniques to distinguish between primary disease and re-infection on a population level. Current estimates of the contribution of re-infection therefore rely on mathematical models which identify the parameters most consistent with epidemiological data; these studies find that exogenous re-infection is important only when TB incidence is high. A basic assumption of these models is that people in a population are all equally likely to come into contact with an infectious case. However, theoretical studies demonstrate that the social and spatial structure can strongly influence the dynamics of infectious disease transmission. Here, we use a network model of TB transmission to evaluate the impact of non-homogeneous mixing on the relative contribution of re-infection over realistic epidemic trajectories. In contrast to the findings of previous models, our results suggest that re-infection may be important in communities where the average disease incidence is moderate or low as the force of infection can be unevenly distributed in the population. These results have important implications for the development of TB control strategies.     

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.     

Integrating stochasticity and network structure into an epidemic model

While the foundations of modern epidemiology are based upon deterministic models with homogeneous mixing, it is being increasingly realized that both spatial structure and stochasticity play major roles in shaping epidemic dynamics. The integration of these two confounding elements is generally ascertained through numerical simulation. Here, for the first time, we develop a more rigorous analytical understanding based on pairwise approximations to incorporate localized spatial structure and diffusion approximations to capture the impact of stochasticity. Our results allow us to quantify, analytically, the impact of network structure on the variability of an epidemic. Using the susceptible–infectious–susceptible framework for the infection dynamics, the pairwise stochastic model is compared with the stochastic homogeneous-mixing (mean-field) model—although to enable a fair comparison the homogeneous-mixing parameters are scaled to give agreement with the pairwise dynamics. At equilibrium, we show that the pairwise model always displays greater variation about the mean, although the differences are generally small unless the prevalence of infection is low. By contrast, during the early epidemic growth phase when the level of infection is increasing exponentially, the pairwise model generally shows less variation.     

Social encounter networks: collective properties and disease transmission

A fundamental challenge of modern infectious disease epidemiology is to quantify the networks of social and physical contacts through which transmission can occur. Understanding the collective properties of these interactions is critical for both accurate prediction of the spread of infection and determining optimal control measures. However, even the basic properties of such networks are poorly quantified, forcing predictions to be made based on strong assumptions concerning network structure. Here, we report on the results of a large-scale survey of social encounters mainly conducted in Great Britain. First, we characterize the distribution of contacts, which possesses a lognormal body and a power-law tail with an exponent of −2.45; we provide a plausible mechanistic model that captures this form. Analysis of the high level of local clustering of contacts reveals additional structure within the network, implying that social contacts are degree assortative. Finally, we describe the epidemiological implications of this local network structure: these contradict the usual predictions from networks with heavy-tailed degree distributions and contain public-health messages about control. Our findings help us to determine the types of realistic network structure that should be assumed in future population level studies of infection transmission, leading to better interpretations of epidemiological data and more appropriate policy decisions.     

Space and contact networks: capturing the locality of disease transmission.

While an arbitrary level of complexity may be included in simulations of spatial epidemics, computational intensity and analytical intractability mean that such models often lack transparency into the determinants of epidemiological dynamics. Although numerous approaches attempt to resolve this complexity-tractability trade-off, moment closure methods arguably offer the most promising and robust frameworks for capturing the role of the locality of contact processes on global disease dynamics. While a close analogy may be made between full stochastic spatial transmission models and dynamic network models, we consider here the special case where the dynamics of the network topology change on time-scales much longer than the epidemiological processes imposed on them; in such cases, the use of static network models are justified. We show that in such cases, static network models may provide excellent approximations to the underlying spatial contact process through an appropriate choice of the effective neighbourhood size. We also demonstrate the robustness of this mapping by examining the equivalence of deterministic approximations to the full spatial and network models derived under third-order moment closure assumptions. For systems where deviation from homogeneous mixing is limited, we show that pair equations developed for network models are at least as good an approximation to the underlying stochastic spatial model as more complex spatial moment equations, with both classes of approximation becoming less accurate only for highly localized kernels.     

The effect of network mixing patterns on epidemic dynamics and the efficacy of disease contact tracing

In networks, nodes may preferentially contact other nodes with similar (assortatively mixed) or dissimilar (disassortatively mixed) numbers of contacts. Different patterns of contact support different epidemic dynamics, potentially affecting the efficacy of control measures such as contact tracing, which aims to identify and isolate nodes with infectious contacts. We used stochastic simulations to investigate the effects of mixing patterns on epidemic dynamics and contact-tracing efficacy. For uncontrolled epidemics, outbreaks occur at lower infection rates for more assortatively mixed networks, with faster initial epidemic growth rate and shorter epidemic duration than for disassortatively mixed networks. Contact tracing performs better for assortative mixing where epidemic size is large and tracing rate low, but it performs better for disassortative mixing at higher contact rates. For assortatively mixed networks, disease spreads first to highly connected nodes, but this is balanced by contact tracing quickly identifying these same nodes. The converse is true for disassortative mixing, where both disease and tracing are less likely to target highly connected nodes. For small epidemics, contact tracing is more effective on disassortative networks due to the greater resilience of assortative networks to link removal. Multi-step contact tracing is more effective than single-step tracing for assortative mixing, but this effect is smaller for disassortatively mixed networks.     

Stochastic fluctuations in epidemics on networks.

The effects of demographic stochasticity on the long-term behaviour of endemic infectious diseases have been considered for long as a necessary addition to an underlying deterministic theory. The latter would explain the regular behaviour of recurrent epidemics and the former the superimposed noise of observed incidence patterns. Recently, a stochastic theory based on a mechanism of resonance with internal noise has shifted the role of stochasticity closer to the centre stage, by showing that the major dynamic patterns found in the incidence data can be explained as resonant fluctuations, whose behaviour is largely independent of the amplitude of seasonal forcing, and by contrast very sensitive to the basic epidemiological parameters. Here we elaborate on that approach, by adding an ingredient which is missing in standard epidemic models, the 'mixing network' through which infection may propagate. We find that spatial correlations have a major effect on the enhancement of the amplitude and the coherence of the resonant stochastic fluctuations, providing the ordered patterns of recurrent epidemics, whose period may differ significantly from that of the small oscillations around the deterministic equilibrium. We also show that the inclusion of a more realistic, time-correlated recovery profile instead of exponentially distributed infectious periods may, even in the random-mixing limit, contribute to the same effect.     

Rates of coalescence for common epidemiological models at equilibrium

Coalescent theory provides a mathematical framework for quantitatively interpreting gene genealogies. With the increased availability of molecular sequence data, disease ecologists now regularly apply this body of theory to viral phylogenies, most commonly in attempts to reconstruct demographic histories of infected individuals and to estimate parameters such as the basic reproduction number. However, with few exceptions, the mathematical expressions at the core of coalescent theory have not been explicitly linked to the structure of epidemiological models, which are commonly used to mathematically describe the transmission dynamics of a pathogen. Here, we aim to make progress towards establishing this link by presenting a general approach for deriving a model''s rate of coalescence under the assumption that the disease dynamics are at their endemic equilibrium. We apply this approach to four common families of epidemiological models: standard susceptible-infected-susceptible/susceptible-infected-recovered/susceptible-infected-recovered-susceptible models, models with individual heterogeneity in infectivity, models with an exposed but not yet infectious class and models with variable distributions of the infectious period. These results improve our understanding of how epidemiological processes shape viral genealogies, as well as how these processes affect levels of viral diversity and rates of genetic drift. Finally, we discuss how a subset of these coalescent rate expressions can be used for phylodynamic inference in non-equilibrium settings. For the ones that are limited to equilibrium conditions, we also discuss why this is the case. These results, therefore, point towards necessary future work while providing intuition on how epidemiological characteristics of the infection process impact gene genealogies.     

Threshold parameters for a model of epidemic spread among households and workplaces

The basic reproduction number R ₀ is one of the most important concepts in modern infectious disease epidemiology. However, for more realistic and more complex models than those assuming homogeneous mixing in the population, other threshold quantities can be defined that are sometimes more useful and easily derived in terms of model parameters. In this paper, we present a model for the spread of a permanently immunizing infection in a population socially structured into households and workplaces/schools, and we propose and discuss a new household-to-household reproduction number R _H for it. We show how R _H overcomes some of the limitations of a previously proposed threshold parameter, and we highlight its relationship with the effort required to control an epidemic when interventions are targeted at randomly selected households.     

Stochastic amplification in epidemics

The role of stochasticity and its interplay with nonlinearity are central current issues in studies of the complex population patterns observed in nature, including the pronounced oscillations of wildlife and infectious diseases. The dynamics of childhood diseases have provided influential case studies to develop and test mathematical models with practical application to epidemiology, but are also of general relevance to the central question of whether simple nonlinear systems can explain and predict the complex temporal and spatial patterns observed in nature outside laboratory conditions. Here, we present a stochastic theory for the major dynamical transitions in epidemics from regular to irregular cycles, which relies on the discrete nature of disease transmission and low spatial coupling. The full spectrum of stochastic fluctuations is derived analytically to show how the amplification of noise varies across these transitions. The changes in noise amplification and coherence appear robust to seasonal forcing, questioning the role of seasonality and its interplay with deterministic components of epidemiological models. Childhood diseases are shown to fall into regions of parameter space of high noise amplification. This type of "endogenous" stochastic resonance may be relevant to population oscillations in nonlinear ecological 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.     

Recent network evolution increases the potential for large epidemics in the British cattle population.

Following the foot and mouth disease epidemic in Great Britain (GB) in 2001, livestock movement bans were replaced with mandatory periods of standstill for livestock moving between premises. It was anticipated that these movement restrictions would limit each individual's contact networks, the extent of livestock movements and thus the spread of future disease outbreaks. However, the effect of behaviour changes on the global network in adapting to these restrictions is currently unknown. Here, we take a novel approach using GB cattle movement data to construct week-by-week contact networks between animal holdings (AH) to explore the evolution of the network since this policy was introduced, the first time network theory has been used for this purpose. We show that the number of AH moving cattle as part of the giant strong component (GSC), representing the region of maximal connectivity, has been increasing linearly over time. This is of epidemiological significance as the size of the GSC indicates the number of holdings potentially exposed to disease, thus giving a lower bound of maximum epidemic size. Therefore, despite restriction of cattle movements, emergent behaviour in this self-organizing system has potentially increased the size of infectious disease epidemics within the cattle industry.     

The network of sheep movements within Great Britain: Network properties and their implications for infectious disease spread.

During the 2001 foot and mouth disease epidemic in the UK, initial dissemination of the disease to widespread geographical regions was attributed to livestock movement, especially of sheep. In response, recording schemes to provide accurate data describing the movement of large livestock in Great Britain (GB) were introduced. Using these data, we reconstruct directed contact networks within the sheep industry and identify key epidemiological properties of these networks. There is clear seasonality in sheep movements, with a peak of intense activity in August and September and an associated high risk of a large epidemic. The high correlation between the in and out degree of nodes favours disease transmission. However, the contact networks were largely dissasortative: highly connected nodes mostly connect to nodes with few contacts, effectively slowing the spread of disease. This is a result of bipartite-like network properties, with most links occurring between highly active markets and less active farms. When comparing sheep movement networks (SMNs) to randomly generated networks with the same number of nodes and node degrees, despite structural differences (such as disassortativity and higher frequency of even path lengths in the SMNs), the characteristic path lengths within the SMNs are close to values computed from the corresponding random networks, showing that SMNs have 'small-world'-like properties. Using the network properties, we show that targeted biosecurity or surveillance at highly connected nodes would be highly effective in preventing a large and widespread epidemic.     

On methods for studying stochastic disease dynamics.

Models that deal with the individual level of populations have shown the importance of stochasticity in ecology, epidemiology and evolution. An increasingly common approach to studying these models is through stochastic (event-driven) simulation. One striking disadvantage of this approach is the need for a large number of replicates to determine the range of expected behaviour. Here, for a class of stochastic models called Markov processes, we present results that overcome this difficulty and provide valuable insights, but which have been largely ignored by applied researchers. For these models, the so-called Kolmogorov forward equation (also called the ensemble or master equation) allows one to simultaneously consider the probability of each possible state occurring. Irrespective of the complexities and nonlinearities of population dynamics, this equation is linear and has a natural matrix formulation that provides many analytical insights into the behaviour of stochastic populations and allows rapid evaluation of process dynamics. Here, using epidemiological models as a template, these ensemble equations are explored and results are compared with traditional stochastic simulations. In addition, we describe further advantages of the matrix formulation of dynamics, providing simple exact methods for evaluating expected eradication (extinction) times of diseases, for comparing expected total costs of possible control programmes and for estimation of disease parameters.     

The scaling of human interactions with city size

The size of cities is known to play a fundamental role in social and economic life. Yet, its relation to the structure of the underlying network of human interactions has not been investigated empirically in detail. In this paper, we map society-wide communication networks to the urban areas of two European countries. We show that both the total number of contacts and the total communication activity grow superlinearly with city population size, according to well-defined scaling relations and resulting from a multiplicative increase that affects most citizens. Perhaps surprisingly, however, the probability that an individual''s contacts are also connected with each other remains largely unaffected. These empirical results predict a systematic and scale-invariant acceleration of interaction-based spreading phenomena as cities get bigger, which is numerically confirmed by applying epidemiological models to the studied networks. Our findings should provide a microscopic basis towards understanding the superlinear increase of different socioeconomic quantities with city size, that applies to almost all urban systems and includes, for instance, the creation of new inventions or the prevalence of certain contagious diseases.     

Perspectives on the basic reproductive ratio.

The basic reproductive ratio, R0, is defined as the expected number of secondary infections arising from a single individual during his or her entire infectious period, in a population of susceptibles. This concept is fundamental to the study of epidemiology and within-host pathogen dynamics. Most importantly, R0 often serves as a threshold parameter that predicts whether an infection will spread. Related parameters which share this threshold behaviour, however, may or may not give the true value of R0. In this paper we give a brief overview of common methods of formulating R0 and surrogate threshold parameters from deterministic, non-structured models. We also review common means of estimating R0 from epidemiological data. Finally, we survey the recent use of R0 in assessing emerging diseases, such as severe acute respiratory syndrome and avian influenza, a number of recent livestock diseases, and vector-borne diseases malaria, dengue and West Nile virus.     

Likelihood-based estimation of continuous-time epidemic models from time-series data: application to measles transmission in London

We present a new statistical approach to analyse epidemic time-series data. A major difficulty for inference is that (i) the latent transmission process is partially observed and (ii) observed quantities are further aggregated temporally. We develop a data augmentation strategy to tackle these problems and introduce a diffusion process that mimics the susceptible-infectious-removed (SIR) epidemic process, but that is more tractable analytically. While methods based on discrete-time models require epidemic and data collection processes to have similar time scales, our approach, based on a continuous-time model, is free of such constraint. Using simulated data, we found that all parameters of the SIR model, including the generation time, were estimated accurately if the observation interval was less than 2.5 times the generation time of the disease. Previous discrete-time TSIR models have been unable to estimate generation times, given that they assume the generation time is equal to the observation interval. However, we were unable to estimate the generation time of measles accurately from historical data. This indicates that simple models assuming homogenous mixing (even with age structure) of the type which are standard in mathematical epidemiology miss key features of epidemics in large populations.     

Spread of infectious disease through clustered populations

Networks of person-to-person contacts form the substrate along which infectious diseases spread. Most network-based studies of this spread focus on the impact of variations in degree (the number of contacts an individual has). However, other effects such as clustering, variations in infectiousness or susceptibility, or variations in closeness of contacts may play a significant role. We develop analytic techniques to predict how these effects alter the growth rate, probability and size of epidemics, and validate the predictions with a realistic social network. We find that (for a given degree distribution and average transmissibility) clustering is the dominant factor controlling the growth rate, heterogeneity in infectiousness is the dominant factor controlling the probability of an epidemic and heterogeneity in susceptibility is the dominant factor controlling the size of an epidemic. Edge weights (measuring closeness or duration of contacts) have impact only if correlations exist between different edges. Combined, these effects can play a minor role in reinforcing one another, with the impact of clustering the largest when the population is maximally heterogeneous or if the closer contacts are also strongly clustered. Our most significant contribution is a systematic way to address clustering in infectious disease models, and our results have a number of implications for the design of interventions.