Joshua V. Ross

Networks and the epidemiology of infectious disease

The science of networks has revolutionised research into the dynamics of interacting elements. It could be argued that epidemiology in particular has embraced the potential of network theory more than any other discipline. Here we review the growing body of research concerning the spread of infectious diseases on networks, focusing on the interplay between network theory and epidemiology. The review is split into four main sections, which examine: the types of network relevant to epidemiology; the multitude of ways these networks can be characterised; the statistical methods that can be applied to infer the epidemiological parameters on a realised network; and finally simulation and analytical methods to determine epidemic dynamics on a given network. Given the breadth of areas covered and the ever-expanding number of publications, a comprehensive review of all work is impossible. Instead, we provide a personalised overview into the areas of network epidemiology that have seen the greatest progress in recent years or have the greatest potential to provide novel insights. As such, considerable importance is placed on analytical approaches and statistical methods which are both rapidly expanding fields. Throughout this review we restrict our attention to epidemiological issues.

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.

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.

The chlamydia screening studies: rationale and design

Background: Screening has been recommended to reduce the prevalence and morbidity associated with genital chlamydia infection in the United Kingdom. Methods: We describe the rationale and study design of the Chlamydia Screening Studies (ClaSS), a collaborative project designed to evaluate screening outside genitourinary medicine clinics. A non-selective, active screening approach in 16–39 year olds randomly sampled from 27 general practice lists in the Bristol and Birmingham areas formed the basis of interlinked studies: a case-control study was used to investigate factors to improve the targeting of screening; participants with chlamydia were invited to enrol in a randomised controlled trial to evaluate partner notification conducted in primary care; and laboratory based studies were used to assess the best specimens and tests. We also explored psychosocial effects of screening and partner notification and modelled the cost effectiveness of the programme. Conclusion: Results from four pilot practices show that mailing of specimens for chlamydia testing is feasible but that it is difficult to achieve high response rates with postal screening. The high prevalence of asymptomatic infection in men suggests that efforts to screen men for chlamydia should be strengthened.

Measuring social networks in British primary schools through scientific engagement

Primary schools constitute a key risk group for the transmission of infectious diseases, concentrating great numbers of immunologically naive individuals at high densities. Despite this, very little is known about the social patterns of mixing within a school, which are likely to contribute to disease transmission. In this study, we present a novel approach where scientific engagement was used as a tool to access school populations and measure social networks between young (4–11 years) children. By embedding our research project within enrichment activities to older secondary school (13–15) children, we could exploit the existing links between schools to achieve a high response rate for our study population (around 90% in most schools). Social contacts of primary school children were measured through self-reporting based on a questionnaire design, and analysed using the techniques of social network analysis. We find evidence of marked social structure and gender assortativity within and between classrooms in the same school. These patterns have been previously reported in smaller studies, but to our knowledge no study has attempted to exhaustively sample entire school populations. Our innovative approach facilitates access to a vitally important (but difficult to sample) epidemiological sub-group. It provides a model whereby scientific communication can be used to enhance, rather than merely complement, the outcomes of research.

Calculation of disease dynamics in a population of households.

Early mathematical representations of infectious disease dynamics assumed a single, large, homogeneously mixing population. Over the past decade there has been growing interest in models consisting of multiple smaller subpopulations (households, workplaces, schools, communities), with the natural assumption of strong homogeneous mixing within each subpopulation, and weaker transmission between subpopulations. Here we consider a model of SIRS (susceptible-infectious-recovered-susceptible) infection dynamics in a very large (assumed infinite) population of households, with the simplifying assumption that each household is of the same size (although all methods may be extended to a population with a heterogeneous distribution of household sizes). For this households model we present efficient methods for studying several quantities of epidemiological interest: (i) the threshold for invasion; (ii) the early growth rate; (iii) the household offspring distribution; (iv) the endemic prevalence of infection; and (v) the transient dynamics of the process. We utilize these methods to explore a wide region of parameter space appropriate for human infectious diseases. We then extend these results to consider the effects of more realistic gamma-distributed infectious periods. We discuss how all these results differ from standard homogeneous-mixing models and assess the implications for the invasion, transmission and persistence of infection. The computational efficiency of the methodology presented here will hopefully aid in the parameterisation of structured models and in the evaluation of appropriate responses for future disease outbreaks.

Simulation-based Bayesian inference for epidemic models

A powerful and flexible method for fitting dynamic models to missing and censored data is to use the Bayesian paradigm via data-augmented Markov chain Monte Carlo (DA-MCMC). This samples from the joint posterior for the parameters and missing data, but requires high memory overheads for large-scale systems. In addition, designing efficient proposal distributions for the missing data is typically challenging. Pseudo-marginal methods instead integrate across the missing data using a Monte Carlo estimate for the likelihood, generated from multiple independent simulations from the model. These techniques can avoid the high memory requirements of DA-MCMC, and under certain conditions produce the exact marginal posterior distribution for parameters. A novel method is presented for implementing importance sampling for dynamic epidemic models, by conditioning the simulations on sets of validity criteria (based on the model structure) as well as the observed data. The flexibility of these techniques is illustrated using both removal time and final size data from an outbreak of smallpox. It is shown that these approaches can circumvent the need for reversible-jump MCMC, and can allow inference in situations where DA-MCMC is impossible due to computationally infeasible likelihoods.

Epidemiological consequences of household-based antiviral prophylaxis for pandemic influenza.

Antiviral treatment offers a fast acting alternative to vaccination; as such it is viewed as a first-line of defence against pandemic influenza in protecting families and households once infection has been detected. In clinical trials, antiviral treatments have been shown to be efficacious in preventing infection, limiting disease and reducing transmission, yet their impact at containing the 2009 influenza A(H1N1)pdm outbreak was limited. To understand this seeming discrepancy, we develop a general and computationally efficient model for studying household-based interventions. This allows us to account for uncertainty in quantities relevant to the 2009 pandemic in a principled way, accounting for the heterogeneity and variability in each epidemiological process modelled. We find that the population-level effects of delayed antiviral treatment and prophylaxis mean that their limited overall impact is quantitatively consistent (at current levels of precision) with their reported clinical efficacy under ideal conditions. Hence, effective control of pandemic influenza with antivirals is critically dependent on early detection and delivery ideally within 24 h.

Understanding the biological invasion risk posed by the global wildlife trade: propagule pressure drives the introduction and establishment of Nearctic turtles

Biological invasions are a key component of human-induced global change. The continuing increase in global wildlife trade has raised concerns about the parallel increase in the number of new invasive species. However, the factors that link the wildlife trade to the biological invasion process are still poorly understood. Moreover, there are analytical challenges in researching the role of global wildlife trade in biological invasions, particularly issues related to the under-reporting of introduced and established populations in areas with reduced sampling effort. In this work, we use high-quality data on the international trade in Nearctic turtles (1999-2009) coupled with a statistical modelling framework, which explicitly accounts for detection, to investigate the factors that influence the introduction (release, or escape into the wild) of globally traded Nearctic turtles and the establishment success (self-sustaining exotic populations) of slider turtles (Trachemys scripta), the most frequently traded turtle species. We found that the introduction of a species was influenced by the total number of turtles exported to a jurisdiction and the age at maturity of the species, while the establishment success of slider turtles was best associated with the propagule number (number of release events), and the number of native turtles in the jurisdiction of introduction. These results indicate both a direct and indirect association between the wildlife trade and the introduction of turtles and establishment success of slider turtles, respectively. Our results highlight the existence of gaps in the number of globally recorded introduction events and established populations of slider turtles, although the expected bias is low. We emphasize the importance of researching independently the factors that affect the different stages of the invasion pathway. Critically, we observe that the number of traded individuals might not always be an adequate proxy for propagule pressure and establishment success.

How big is an outbreak likely to be? Methods for epidemic final-size calculation

Epidemic models have become a routinely used tool to inform policy on infectious disease. A particular interest at the moment is the use of computationally intensive inference to parametrize these models. In this context, numerical efficiency is critically important. We consider methods for evaluating the probability mass function of the total number of infections over the course of a stochastic epidemic, with a focus on homogeneous finite populations, but also considering heterogeneous and large populations. Relevant methods are reviewed critically, with existing and novel extensions also presented. We provide code in Matlab and a systematic comparison of numerical efficiency.