Epidemic Spread on Weighted Networks

The contact structure between hosts shapes disease spread. Most network-based models used in epidemiology tend to ignore heterogeneity in the weighting of contacts between two individuals. However, this assumption is known to be at odds with the data for many networks (e.g. sexual contact networks) and to have a critical influence on epidemics'' behavior. One of the reasons why models usually ignore heterogeneity in transmission is that we currently lack tools to analyze weighted networks, such that most studies rely on numerical simulations. Here, we present a novel framework to estimate key epidemiological variables, such as the rate of early epidemic expansion () and the basic reproductive ratio (), from joint probability distributions of number of partners (contacts) and number of interaction events through which contacts are weighted. These distributions are much easier to infer than the exact shape of the network, which makes the approach widely applicable. The framework also allows for a derivation of the full time course of epidemic prevalence and contact behaviour, which we validate with numerical simulations on networks. Overall, incorporating more realistic contact networks into epidemiological models can improve our understanding of the emergence and spread of infectious diseases.     

Building epidemiological models from R0: an implicit treatment of transmission in networks

Simple deterministic models are still at the core of theoretical epidemiology despite the increasing evidence for the importance of contact networks underlying transmission at the individual level. These mean-field or 'compartmental' models based on homogeneous mixing have made, and continue to make, important contributions to the epidemiology and the ecology of infectious diseases but fail to reproduce many of the features observed for disease spread in contact networks. In this work, we show that it is possible to incorporate the important effects of network structure on disease spread with a mean-field model derived from individual level considerations. We propose that the fundamental number known as the basic reproductive number of the disease, R0, which is typically derived as a threshold quantity, be used instead as a central parameter to construct the model from. We show that reliable estimates of individual level parameters can replace a detailed knowledge of network structure, which in general may be difficult to obtain. We illustrate the proposed model with small world networks and the classical example of susceptible-infected-recovered (SIR) epidemics.     

Network-based analysis of stochastic SIR epidemic models with random and proportionate mixing

In this paper, we outline the theory of epidemic percolation networks and their use in the analysis of stochastic susceptible-infectious-removed (SIR) epidemic models on undirected contact networks. We then show how the same theory can be used to analyze stochastic SIR models with random and proportionate mixing. The epidemic percolation networks for these models are purely directed because undirected edges disappear in the limit of a large population. In a series of simulations, we show that epidemic percolation networks accurately predict the mean outbreak size and probability and final size of an epidemic for a variety of epidemic models in homogeneous and heterogeneous populations. Finally, we show that epidemic percolation networks can be used to re-derive classical results from several different areas of infectious disease epidemiology. In an Appendix, we show that an epidemic percolation network can be defined for any time-homogeneous stochastic SIR model in a closed population and prove that the distribution of outbreak sizes given the infection of any given node in the SIR model is identical to the distribution of its out-component sizes in the corresponding probability space of epidemic percolation networks. We conclude that the theory of percolation on semi-directed networks provides a very general framework for the analysis of stochastic SIR models in closed populations.     

Pair-level approximations to the spatio-temporal dynamics of epidemics on asymmetric contact networks

The process of infection during an epidemic can be envisaged as being transmitted via a network of routes represented by a contact network. Most differential equation models of epidemics are mean-field models. These contain none of the underlying spatial structure of the contact network. By extending the mean-field models to pair-level, some of the spatial structure can be contained in the model. Some networks of transmission such as river or transportation networks are clearly asymmetric, whereas others such as airborne infection can be regarded as symmetric. Pair-level models have been developed to describe symmetric contact networks. Here we report on work to develop a pair-level model that is also applicable to asymmetric contact networks. The procedure for closing the model at the level of pairs is discussed in detail. The model is compared against stochastic simulations of epidemics on asymmetric contact networks and against the predictions of the symmetric model on the same networks. DEFRA funded project FC1153     

Congruent epidemic models for unstructured and structured populations: analytical reconstruction of a 2003 SARS outbreak

Both the threat of bioterrorism and the natural emergence of contagious diseases underscore the importance of quantitatively understanding disease transmission in structured human populations. Over the last few years, researchers have advanced the mathematical theory of scale-free networks and used such theoretical advancements in pilot epidemic models. Scale-free contact networks are particularly interesting in the realm of mathematical epidemiology, primarily because these networks may allow meaningfully structured populations to be incorporated in epidemic models at moderate or intermediate levels of complexity. Moreover, a scale-free contact network with node degree correlation is in accord with the well-known preferred mixing concept. The present author describes a semi-empirical and deterministic epidemic modeling approach that (a) focuses on time-varying rates of disease transmission in both unstructured and structured populations and (b) employs probability density functions to characterize disease progression and outbreak controls. Given an epidemic curve for a historical outbreak, this modeling approach calls for Monte Carlo calculations (that define the average new infection rate) and solutions to integro-differential equations (that describe outbreak dynamics in an aggregate population or across all network connectivity classes). Numerical results are obtained for the 2003 SARS outbreak in Taiwan and the dynamical implications of time-varying transmission rates and scale-free contact networks are discussed in some detail.     

Multiple sources and routes of information transmission: Implications for epidemic dynamics

In a recent paper [20], we proposed and analyzed a compartmental ODE-based model describing the dynamics of an infectious disease where the presence of the pathogen also triggers the diffusion of information about the disease. In this paper, we extend this previous work by presenting results based on pairwise and simulation models that are better suited for capturing the population contact structure at a local level. We use the pairwise model to examine the potential of different information generating mechanisms and routes of information transmission to stop disease spread or to minimize the impact of an epidemic. The individual-based simulation is used to better differentiate between the networks of disease and information transmission and to investigate the impact of different basic network topologies and network overlap on epidemic dynamics. The paper concludes with an individual-based semi-analytic calculation of R₀ at the non-trivial disease free equilibrium.     

Disease contact tracing in random and clustered networks

The efficacy of contact tracing, be it between individuals (e.g. sexually transmitted diseases or severe acute respiratory syndrome) or between groups of individuals (e.g. foot-and-mouth disease; FMD), is difficult to evaluate without precise knowledge of the underlying contact structure; i.e. who is connected to whom? Motivated by the 2001 FMD epidemic in the UK, we determine, using stochastic simulations and deterministic 'moment closure' models of disease transmission on networks of premises (nodes), network and disease properties that are important for contact tracing efficiency. For random networks with a high average number of connections per node, little clustering of connections and short latency periods, contact tracing is typically ineffective. In this case, isolation of infected nodes is the dominant factor in determining disease epidemic size and duration. If the latency period is longer and the average number of connections per node small, or if the network is spatially clustered, then the contact tracing performs better and an overall reduction in the proportion of nodes that are removed during an epidemic is observed.     

Social encounter networks: characterizing Great Britain

A major goal of infectious disease epidemiology is to understand and predict the spread of infections within human populations, with the intention of better informing decisions regarding control and intervention. However, the development of fully mechanistic models of transmission requires a quantitative understanding of social interactions and collective properties of social networks. We performed a cross-sectional study of the social contacts on given days for more than 5000 respondents in England, Scotland and Wales, through postal and online survey methods. The survey was designed to elicit detailed and previously unreported measures of the immediate social network of participants relevant to infection spread. Here, we describe individual-level contact patterns, focusing on the range of heterogeneity observed and discuss the correlations between contact patterns and other socio-demographic factors. We find that the distribution of the number of contacts approximates a power-law distribution, but postulate that total contact time (which has a shorter-tailed distribution) is more epidemiologically relevant. We observe that children, public-sector and healthcare workers have the highest number of total contact hours and are therefore most likely to catch and transmit infectious disease. Our study also quantifies the transitive connections made between an individual''s contacts (or clustering); this is a key structural characteristic of social networks with important implications for disease transmission and control efficacy. Respondents'' networks exhibit high levels of clustering, which varies across social settings and increases with duration, frequency of contact and distance from home. Finally, we discuss the implications of these findings for the transmission and control of pathogens spread through close contact.     

The effect of contact heterogeneity and multiple routes of transmission on final epidemic size

Heterogeneity in the number of potentially infectious contacts amongst members of a population increases the basic reproduction ratio (R(0)) and markedly alters disease dynamics compared to traditional mean-field models. Most models describing transmission on contact networks only account for one specific route of transmission. However, for many infectious diseases multiple routes of transmission exist. The model presented here captures transmission through a well defined network of contacts, complemented by mean-field type transmission amongst the nodes of the network that accounts for alternative routes of transmission. The impact of these combined transmission mechanisms on the final epidemic size is investigated analytically. The analytic predictions for the purely mean-field case and the transmission through the network-only case are confirmed by individual-based network simulations. There is a critical transmission potential above which an increased contribution of the mean-field type transmission increases the final epidemic size while an increased contribution of the transmission through the network decreases it. Below the critical transmission potential the opposite effect is observed.     

Susceptible-infected-recovered epidemics in dynamic contact networks

Contact patterns in populations fundamentally influence the spread of infectious diseases. Current mathematical methods for epidemiological forecasting on networks largely assume that contacts between individuals are fixed, at least for the duration of an outbreak. In reality, contact patterns may be quite fluid, with individuals frequently making and breaking social or sexual relationships. Here, we develop a mathematical approach to predicting disease transmission on dynamic networks in which each individual has a characteristic behaviour (typical contact number), but the identities of their contacts change in time. We show that dynamic contact patterns shape epidemiological dynamics in ways that cannot be adequately captured in static network models or mass-action models. Our new model interpolates smoothly between static network models and mass-action models using a mixing parameter, thereby providing a bridge between disparate classes of epidemiological models. Using epidemiological and sexual contact data from an Atlanta high school, we demonstrate the application of this method for forecasting and controlling sexually transmitted disease outbreaks.     

Large-scale properties of clustered networks: implications for disease dynamics

We consider previously proposed procedures for generating clustered networks and investigate how these procedures lead to differences in network properties other than clustering. We interpret our findings in terms of the effect of the network structure on the disease outbreak threshold and disease dynamics. To generate null-model networks for comparison, we implement an assortativity-conserving rewiring algorithm that alters the level of clustering while causing minimal impact on other properties. We show that many theoretical network models used to generate networks with a particular property often lead to significant changes in network properties other than that of interest. For high levels of clustering, different procedures lead to networks that differ in degree heterogeneity and assortativity, and in broader scale measures such as ?(0) and the distribution of shortest path lengths. Hence, care must be taken when investigating the implications of network properties for disease transmission or other dynamic process that the network supports.     

Inferring social network structure from bacterial sequence data

Using DNA sequence data from pathogens to infer transmission networks has traditionally been done in the context of epidemics and outbreaks. Sequence data could analogously be applied to cases of ubiquitous commensal bacteria; however, instead of inferring chains of transmission to track the spread of a pathogen, sequence data for bacteria circulating in an endemic equilibrium could be used to infer information about host contact networks. Here, we show--using simulated data--that multilocus DNA sequence data, based on multilocus sequence typing schemes (MLST), from isolates of commensal bacteria can be used to infer both local and global properties of the contact networks of the populations being sampled. Specifically, for MLST data simulated from small-world networks, the small world parameter controlling the degree of structure in the contact network can robustly be estimated. Moreover, we show that pairwise distances in the network--degrees of separation--correlate with genetic distances between isolates, so that how far apart two individuals in the network are can be inferred from MLST analysis of their commensal bacteria. This result has important consequences, and we show an example from epidemiology: how this result could be used to test for infectious origins of diseases of unknown etiology.     

Using network properties to predict disease dynamics on human contact networks

Recent studies have increasingly turned to graph theory to model more realistic contact structures that characterize disease spread. Because of the computational demands of these methods, many researchers have sought to use measures of network structure to modify analytically tractable differential equation models. Several of these studies have focused on the degree distribution of the contact network as the basis for their modifications. We show that although degree distribution is sufficient to predict disease behaviour on very sparse or very dense human contact networks, for intermediate density networks we must include information on clustering and path length to accurately predict disease behaviour. Using these three metrics, we were able to explain more than 98 per cent of the variation in endemic disease levels in our stochastic simulations.     

Moment Approximation of Infection Dynamics in a Population of Moving Hosts

The modelling of contact processes between hosts is of key importance in epidemiology. Current studies have mainly focused on networks with stationary structures, although we know these structures to be dynamic with continuous appearance and disappearance of links over time. In the case of moving individuals, the contact network cannot be established. Individual-based models (IBMs) can simulate the individual behaviours involved in the contact process. However, with very large populations, they can be hard to simulate and study due to the computational costs. We use the moment approximation (MA) method to approximate a stochastic IBM with an aggregated deterministic model. We illustrate the method with an application in animal epidemiology: the spread of the highly pathogenic virus H5N1 of avian influenza in a poultry flock. The MA method is explained in a didactic way so that it can be reused and extended. We compare the simulation results of three models: 1. an IBM, 2. a MA, and 3. a mean-field (MF). The results show a close agreement between the MA model and the IBM. They highlight the importance for the models to capture the displacement behaviours and the contact processes in the study of disease spread. We also illustrate an original way of using different models of the same system to learn more about the system itself, and about the representation we build of it.     

Infectious disease transmission and contact networks in wildlife and livestock

The use of social and contact networks to answer basic and applied questions about infectious disease transmission in wildlife and livestock is receiving increased attention. Through social network analysis, we understand that wild animal and livestock populations, including farmed fish and poultry, often have a heterogeneous contact structure owing to social structure or trade networks. Network modelling is a flexible tool used to capture the heterogeneous contacts of a population in order to test hypotheses about the mechanisms of disease transmission, simulate and predict disease spread, and test disease control strategies. This review highlights how to use animal contact data, including social networks, for network modelling, and emphasizes that researchers should have a pathogen of interest in mind before collecting or using contact data. This paper describes the rising popularity of network approaches for understanding transmission dynamics in wild animal and livestock populations; discusses the common mismatch between contact networks as measured in animal behaviour and relevant parasites to match those networks; and highlights knowledge gaps in how to collect and analyse contact data. Opportunities for the future include increased attention to experiments, pathogen genetic markers and novel computational tools.     

Targeted pandemic containment through identifying local contact network bottlenecks

Decision-making about pandemic mitigation often relies upon simulation modelling. Models of disease transmission through networks of contacts–between individuals or between population centres–are increasingly used for these purposes. Real-world contact networks are rich in structural features that influence infection transmission, such as tightly-knit local communities that are weakly connected to one another. In this paper, we propose a new flow-based edge-betweenness centrality method for detecting bottleneck edges that connect nodes in contact networks. In particular, we utilize convex optimization formulations based on the idea of diffusion with p-norm network flow. Using simulation models of COVID-19 transmission through real network data at both individual and county levels, we demonstrate that targeting bottleneck edges identified by the proposed method reduces the number of infected cases by up to 10% more than state-of-the-art edge-betweenness methods. Furthermore, the proposed method is orders of magnitude faster than existing methods.     

A Simulation Study Comparing Epidemic Dynamics on Exponential Random Graph and Edge-Triangle Configuration Type Contact Network Models

We compare two broad types of empirically grounded random network models in terms of their abilities to capture both network features and simulated Susceptible-Infected-Recovered (SIR) epidemic dynamics. The types of network models are exponential random graph models (ERGMs) and extensions of the configuration model. We use three kinds of empirical contact networks, chosen to provide both variety and realistic patterns of human contact: a highly clustered network, a bipartite network and a snowball sampled network of a “hidden population”. In the case of the snowball sampled network we present a novel method for fitting an edge-triangle model. In our results, ERGMs consistently capture clustering as well or better than configuration-type models, but the latter models better capture the node degree distribution. Despite the additional computational requirements to fit ERGMs to empirical networks, the use of ERGMs provides only a slight improvement in the ability of the models to recreate epidemic features of the empirical network in simulated SIR epidemics. Generally, SIR epidemic results from using configuration-type models fall between those from a random network model (i.e., an Erdős-Rényi model) and an ERGM. The addition of subgraphs of size four to edge-triangle type models does improve agreement with the empirical network for smaller densities in clustered networks. Additional subgraphs do not make a noticeable difference in our example, although we would expect the ability to model cliques to be helpful for contact networks exhibiting household structure.     

Revealing mechanisms of infectious disease spread through empirical contact networks

The spread of pathogens fundamentally depends on the underlying contacts between individuals. Modeling the dynamics of infectious disease spread through contact networks, however, can be challenging due to limited knowledge of how an infectious disease spreads and its transmission rate. We developed a novel statistical tool, INoDS (Identifying contact Networks of infectious Disease Spread) that estimates the transmission rate of an infectious disease outbreak, establishes epidemiological relevance of a contact network in explaining the observed pattern of infectious disease spread and enables model comparison between different contact network hypotheses. We show that our tool is robust to incomplete data and can be easily applied to datasets where infection timings of individuals are unknown. We tested the reliability of INoDS using simulation experiments of disease spread on a synthetic contact network and find that it is robust to incomplete data and is reliable under different settings of network dynamics and disease contagiousness compared with previous approaches. We demonstrate the applicability of our method in two host-pathogen systems: Crithidia bombi in bumblebee colonies and Salmonella in wild Australian sleepy lizard populations. INoDS thus provides a novel and reliable statistical tool for identifying transmission pathways of infectious disease spread. In addition, application of INoDS extends to understanding the spread of novel or emerging infectious disease, an alternative approach to laboratory transmission experiments, and overcoming common data-collection constraints.     

Disease evolution on networks: the role of contact structure

Owing to their rapid reproductive rate and the severe penalties for reduced fitness, diseases are under immense evolutionary pressure. Understanding the evolutionary response of diseases in new situations has clear public-health consequences, given the changes in social and movement patterns over recent decades and the increased use of antibiotics. This paper investigates how a disease may adapt in response to the routes of transmission available between infected and susceptible individuals. The potential transmission routes are defined by a computer-generated contact network, which we describe as either local (highly clustered networks where connected individuals are likely to share common contacts) or global (unclustered networks with a high proportion of long-range connections). Evolution towards stable strategies operates through the gradual random mutation of disease traits (transmission rate and infectious period) whenever new infections occur. In contrast to mean-field models, the use of contact networks greatly constrains the evolutionary dynamics. In the local networks, high transmission rates are selected for, as there is intense competition for susceptible hosts between disease progeny. By contrast, global networks select for moderate transmission rates because direct competition between progeny is minimal and a premium is placed upon persistence. All networks show a very slow but steady rise in the infectious period.