Erik M. Volz

SIR dynamics in random networks with heterogeneous connectivity

Random networks with specified degree distributions have been proposed as realistic models of population structure, yet the problem of dynamically modeling SIR-type epidemics in random networks remains complex. I resolve this dilemma by showing how the SIR dynamics can be modeled with a system of three nonlinear ODE’s. The method makes use of the probability generating function (PGF) formalism for representing the degree distribution of a random network and makes use of network-centric quantities such as the number of edges in a well-defined category rather than node-centric quantities such as the number of infecteds or susceptibles. The PGF provides a simple means of translating between network and node-centric variables and determining the epidemic incidence at any time. The theory also provides a simple means of tracking the evolution of the degree distribution among susceptibles or infecteds. The equations are used to demonstrate the dramatic effects that the degree distribution plays on the final size of an epidemic as well as the speed with which it spreads through the population. Power law degree distributions are observed to generate an almost immediate expansion phase yet have a smaller final size compared to homogeneous degree distributions such as the Poisson. The equations are compared to stochastic simulations, which show good agreement with the theory. Finally, the dynamic equations provide an alternative way of determining the epidemic threshold where large-scale epidemics are expected to occur, and below which epidemic behavior is limited to finite-sized outbreaks.

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.

Epidemic thresholds in dynamic contact networks

The reproductive ratio, R0, is a fundamental quantity in epidemiology, which determines the initial increase in an infectious disease in a susceptible host population. In most epidemic models, there is a specific value of R0, the epidemic threshold, above which epidemics are possible, but below which epidemics cannot occur. As the complexity of an epidemic model increases, so too does the difficulty of calculating epidemic thresholds. Here we derive the reproductive ratio and epidemic thresholds for susceptible–infected–recovered (SIR) epidemics in a simple class of dynamic random networks. As in most epidemiological models, R0 depends on two basic epidemic parameters, the transmission and recovery rates. We find that R0 also depends on social parameters, namely the degree distribution that describes heterogeneity in the numbers of concurrent contacts and the mixing parameter that gives the rate at which contacts are initiated and terminated. We show that social mixing fundamentally changes the epidemiological landscape and, consequently, that static network approximations of dynamic networks can be inadequate.

Edge-based compartmental modelling for infectious disease spread

The primary tool for predicting infectious disease spread and intervention effectiveness is the mass action susceptible–infected–recovered model of Kermack & McKendrick. Its usefulness derives largely from its conceptual and mathematical simplicity; however, it incorrectly assumes that all individuals have the same contact rate and partnerships are fleeting. In this study, we introduce edge-based compartmental modelling, a technique eliminating these assumptions. We derive simple ordinary differential equation models capturing social heterogeneity (heterogeneous contact rates) while explicitly considering the impact of partnership duration. We introduce a graphical interpretation allowing for easy derivation and communication of the model and focus on applying the technique under different assumptions about how contact rates are distributed and how long partnerships last.

Phylodynamics of Infectious Disease Epidemics

We present a formalism for unifying the inference of population size from genetic sequences and mathematical models of infectious disease in populations. Virus phylogenies have been used in many recent studies to infer properties of epidemics. These approaches rely on coalescent models that may not be appropriate for infectious diseases. We account for phylogenetic patterns of viruses in susceptible–infected (SI), susceptible–infected–susceptible (SIS), and susceptible–infected–recovered (SIR) models of infectious disease, and our approach may be a viable alternative to demographic models used to reconstruct epidemic dynamics. The method allows epidemiological parameters, such as the reproductive number, to be estimated directly from viral sequence data. We also describe patterns of phylogenetic clustering that are often construed as arising from a short chain of transmissions. Our model reproduces the moments of the distribution of phylogenetic cluster sizes and may therefore serve as a null hypothesis for cluster sizes under simple epidemiological models. We examine a small cross-sectional sample of human immunodeficiency (HIV)-1 sequences collected in the United States and compare our results to standard estimates of effective population size. Estimated prevalence is consistent with estimates of effective population size and the known history of the HIV epidemic. While our model accurately estimates prevalence during exponential growth, we find that periods of decline are harder to identify.

Effects of Heterogeneous and Clustered Contact Patterns on Infectious Disease Dynamics

The spread of infectious diseases fundamentally depends on the pattern of contacts between individuals. Although studies of contact networks have shown that heterogeneity in the number of contacts and the duration of contacts can have far-reaching epidemiological consequences, models often assume that contacts are chosen at random and thereby ignore the sociological, temporal and/or spatial clustering of contacts. Here we investigate the simultaneous effects of heterogeneous and clustered contact patterns on epidemic dynamics. To model population structure, we generalize the configuration model which has a tunable degree distribution (number of contacts per node) and level of clustering (number of three cliques). To model epidemic dynamics for this class of random graph, we derive a tractable, low-dimensional system of ordinary differential equations that accounts for the effects of network structure on the course of the epidemic. We find that the interaction between clustering and the degree distribution is complex. Clustering always slows an epidemic, but simultaneously increasing clustering and the variance of the degree distribution can increase final epidemic size. We also show that bond percolation-based approximations can be highly biased if one incorrectly assumes that infectious periods are homogeneous, and the magnitude of this bias increases with the amount of clustering in the network. We apply this approach to model the high clustering of contacts within households, using contact parameters estimated from survey data of social interactions, and we identify conditions under which network models that do not account for household structure will be biased.

A discrete-choice approach to modeling social influence on individual decision making

Individual decision making is commonly studied using discrete choice models. Models of this type are applied extensively to the study of travel behavior, residential location, and employment decisions, among other topics of interest. A notable characteristic of the underlying economic theory is the assumption that individuals seek to maximize utility on the basis of their personal attributes and the attributes of the alternatives available to them. This approach ignores the interrelated nature of decision making in social situations—in other words, the role that social structures play in shaping behavior. In this paper we describe a multinomial discrete choice approach to analyzing individual behavior in social situations where position in a social network may encourage or discourage different courses of action. By means of a simulation example, we explore some properties of the model, in particular the effect of network topology.

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.

Viral phylodynamics and the search for an 'effective number of infections'

Information on the dynamics of the effective population size over time can be obtained from the analysis of phylogenies, through the application of time-varying coalescent models. This approach has been used to study the dynamics of many different viruses, and has demonstrated a wide variety of patterns, which have been interpreted in the context of changes over time in the ‘effective number of infections’, a quantity proportional to the number of infected individuals. However, for infectious diseases, the rate of coalescence is driven primarily by new transmissions i.e. the incidence, and only indirectly by the number of infected individuals through sampling effects. Using commonly used epidemiological models, we show that the coalescence rate may indeed reflect the number of infected individuals during the initial phase of exponential growth when time is scaled by infectivity, but in general, a single change in time scale cannot be used to estimate the number of infected individuals. This has important implications when integrating phylogenetic data in the context of other epidemiological data.

Distinguishing epidemic waves from disease spillover in a wildlife population

Serengeti lions frequently experience viral outbreaks. In 1994, one-third of Serengeti lions died from canine distemper virus (CDV). Based on the limited epidemiological data available from this period, it has been unclear whether the 1994 outbreak was propagated by lion-to-lion transmission alone or involved multiple introductions from other sympatric carnivore species. More broadly, we do not know whether contacts between lions allow any pathogen with a relatively short infectious period to percolate through the population (i.e. reach epidemic proportions). We built one of the most realistic contact network models for a wildlife population to date, based on detailed behavioural and movement data from a long-term lion study population. The model allowed us to identify previously unrecognized biases in the sparse data from the 1994 outbreak and develop methods for judiciously inferring disease dynamics from typical wildlife samples. Our analysis of the model in light of the 1994 outbreak data strongly suggest that, although lions are sufficiently well connected to sustain epidemics of CDV-like diseases, the 1994 epidemic was fuelled by multiple spillovers from other carnivore species, such as jackals and hyenas.