Effects of the infectious period distribution on predicted transitions in childhood disease dynamics

The population dynamics of infectious diseases occasionally undergo rapid qualitative changes, such as transitions from annual to biennial cycles or to irregular dynamics. Previous work, based on the standard seasonally forced ‘susceptible–exposed–infectious–removed’ (SEIR) model has found that transitions in the dynamics of many childhood diseases result from bifurcations induced by slow changes in birth and vaccination rates. However, the standard SEIR formulation assumes that the stage durations (latent and infectious periods) are exponentially distributed, whereas real distributions are narrower and centred around the mean. Much recent work has indicated that realistically distributed stage durations strongly affect the dynamical structure of seasonally forced epidemic models. We investigate whether inferences drawn from previous analyses of transitions in patterns of measles dynamics are robust to the shapes of the stage duration distributions. As an illustrative example, we analyse measles dynamics in New York City from 1928 to 1972. We find that with a fixed mean infectious period in the susceptible–infectious–removed (SIR) model, the dynamical structure and predicted transitions vary substantially as a function of the shape of the infectious period distribution. By contrast, with fixed mean latent and infectious periods in the SEIR model, the shapes of the stage duration distributions have a less dramatic effect on model dynamical structure and predicted transitions. All these results can be understood more easily by considering the distribution of the disease generation time as opposed to the distributions of individual disease stages. Numerical bifurcation analysis reveals that for a given mean generation time the dynamics of the SIR and SEIR models for measles are nearly equivalent and are insensitive to the shapes of the disease stage distributions.     

Decreasing stochasticity through enhanced seasonality in measles epidemics

Seasonal changes in the environment are known to be important drivers of population dynamics, giving rise to sustained population cycles. However, it is often difficult to measure the strength and shape of seasonal forces affecting populations. In recent years, statistical time-series methods have been applied to the incidence records of childhood infectious diseases in an attempt to estimate seasonal variation in transmission rates, as driven by the pattern of school terms. In turn, school-term forcing was used to show how susceptible influx rates affect the interepidemic period. In this paper, we document the response of measles dynamics to distinct shifts in the parameter regime using previously unexplored records of measles mortality from the early decades of the twentieth century. We describe temporal patterns of measles epidemics using spectral analysis techniques, and point out a marked decrease in birth rates over time. Changes in host demography alone do not, however, suffice to explain epidemiological transitions. By fitting the time-series susceptible–infected–recovered model to measles mortality data, we obtain estimates of seasonal transmission in different eras, and find that seasonality increased over time. This analysis supports theoretical work linking complex population dynamics and the balance between stochastic and deterministic forces as determined by the strength of seasonality.     

Rubella metapopulation dynamics and importance of spatial coupling to the risk of congenital rubella syndrome in Peru

Rubella is generally a mild childhood disease, but infection during early pregnancy may cause spontaneous abortion or congenital rubella syndrome (CRS), which may entail a variety of birth defects. Consequently, understanding the age-structured dynamics of this infection has considerable public health value. Vaccination short of the threshold for local elimination of transmission will increase the average age of infection. Accordingly, the classic concern for this infection is the potential for vaccination to increase incidence in individuals of childbearing age. A neglected aspect of rubella dynamics is how age incidence patterns may be moulded by the spatial dynamics inherent to epidemic metapopulations. Here, we use a uniquely detailed dataset from Peru to explore the implications of this for the burden of CRS. Our results show that the risk of CRS may be particularly severe in small remote regions, a prediction at odds with expectations in the endemic situation, and with implications for the outcome of vaccination. This outcome results directly from the metapopulation context: specifically, extinction–re-colonization dynamics are crucial because they allow for significant leakage of susceptible individuals into the older age classes during inter-epidemic periods with the potential to increase CRS risk by as much as fivefold.     

Multifractal signatures of infectious diseases

Incidence of infection time-series data for the childhood diseases measles, chicken pox, rubella and whooping cough are described in the language of multifractals. We explore the potential of using the wavelet transform maximum modulus (WTMM) method to characterize the multiscale structure of the observed time series and of simulated data generated by the stochastic susceptible-exposed-infectious-recovered (SEIR) epidemic model. The singularity spectra of the observed time series suggest that each disease is characterized by a unique multifractal signature, which distinguishes that particular disease from the others. The wavelet scaling functions confirm that the time series of measles, rubella and whooping cough are clearly multifractal, while chicken pox has a more monofractal structure in time. The stochastic SEIR epidemic model is unable to reproduce the qualitative singularity structure of the reported incidence data: it is too smooth and does not appear to have a multifractal singularity structure. The precise reasons for the failure of the SEIR epidemic model to reproduce the correct multiscale structure of the reported incidence data remain unclear.     

Modelling the long-term dynamics of pre-vaccination pertussis

The dynamics of strongly immunizing childhood infections is still not well understood. Although reports of successful modelling of several data records can be found in the previous literature, the key determinants of the observed temporal patterns have not yet been clearly identified. In particular, different models of immunity waning and degree of protection applied to disease- and vaccine-induced immunity have been debated in the previous literature on pertussis. Here, we study the effect of disease-acquired immunity on the long-term patterns of pertussis prevalence. We compare five minimal models, all of which are stochastic, seasonally forced, well-mixed models of infection, based on susceptible–infective–recovered dynamics in a closed population. These models reflect different assumptions about the immune response of naive hosts, namely total permanent immunity, immunity waning, immunity waning together with immunity boosting, reinfection of recovered and repeat infection after partial immunity waning. The power spectra of the output prevalence time series characterize the long-term dynamics of the models. For epidemiological parameters consistent with published data, the power spectra show quantitative and even qualitative differences, which can be used to test their assumptions by comparison with ensembles of several-decades-long pre-vaccination data records. We illustrate this strategy on two publicly available historical datasets.     

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.     

Introduction of rubella-containing-vaccine to Madagascar: implications for roll-out and local elimination

Few countries in Africa currently include rubella-containing vaccination (RCV) in their immunization schedule. The Global Alliance for Vaccines Initiative (GAVI) recently opened a funding window that has motivated more widespread roll-out of RCV. As countries plan RCV introductions, an understanding of the existing burden, spatial patterns of vaccine coverage, and the impact of patterns of local extinction and reintroduction for rubella will be critical to developing effective programmes. As one of the first countries proposing RCV introduction in part with GAVI funding, Madagascar provides a powerful and timely case study. We analyse serological data from measles surveillance systems to characterize the epidemiology of rubella in Madagascar. Combining these results with data on measles vaccination delivery, we develop an age-structured model to simulate rubella vaccination scenarios and evaluate the dynamics of rubella and the burden of congenital rubella syndrome (CRS) across Madagascar. We additionally evaluate the drivers of spatial heterogeneity in age of infection to identify focal locations where vaccine surveillance should be strengthened and where challenges to successful vaccination introduction are expected. Our analyses indicate that characteristics of rubella in Madagascar are in line with global observations, with an average age of infection near 7 years, and an impact of frequent local extinction with reintroductions causing localized epidemics. Modelling results indicate that introduction of RCV into the routine programme alone may initially decrease rubella incidence but then result in cumulative increases in the burden of CRS in some regions (and transient increases in this burden in many regions). Deployment of RCV with regular supplementary campaigns will mitigate these outcomes. Results suggest that introduction of RCV offers a potential for elimination of rubella in Madagascar, but also emphasize both that targeted vaccination is likely to be a lynchpin of this success, and the public health vigilance that this introduction will require.     

Stochasticity in staged models of epidemics: quantifying the dynamics of whooping cough

Although many stochastic models can accurately capture the qualitative epidemic patterns of many childhood diseases, there is still considerable discussion concerning the basic mechanisms generating these patterns; much of this stems from the use of deterministic models to try to understand stochastic simulations. We argue that a systematic method of analysing models of the spread of childhood diseases is required in order to consistently separate out the effects of demographic stochasticity, external forcing and modelling choices. Such a technique is provided by formulating the models as master equations and using the van Kampen system-size expansion to provide analytical expressions for quantities of interest. We apply this method to the susceptible–exposed–infected–recovered (SEIR) model with distributed exposed and infectious periods and calculate the form that stochastic oscillations take on in terms of the model parameters. With the use of a suitable approximation, we apply the formalism to analyse a model of whooping cough which includes seasonal forcing. This allows us to more accurately interpret the results of simulations and to make a more quantitative assessment of the predictions of the model. We show that the observed dynamics are a result of a macroscopic limit cycle induced by the external forcing and resonant stochastic oscillations about this cycle.     

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.     

A century of transitions in New York City's measles dynamics

Infectious diseases spreading in a human population occasionally exhibit sudden transitions in their qualitative dynamics. Previous work has successfully predicted such transitions in New York City''s historical measles incidence using the seasonally forced susceptible–infectious–recovered (SIR) model. This work relied on a dataset spanning 45 years (1928–1973), which we have extended to 93 years (1891–1984). We identify additional dynamical transitions in the longer dataset and successfully explain them by analysing attractors and transients of the same mechanistic epidemiological model.     

Implications of spatially heterogeneous vaccination coverage for the risk of congenital rubella syndrome in South Africa

Rubella is generally a mild childhood disease, but infection during early pregnancy may cause spontaneous abortion or congenital rubella syndrome (CRS), which may entail a variety of birth defects. Since vaccination at levels short of those necessary to achieve eradication may increase the average age of infection, and thus potentially the CRS burden, introduction of the vaccine has been limited to contexts where coverage is high. Recent work suggests that spatial heterogeneity in coverage should also be a focus of concern. Here, we use a detailed dataset from South Africa to explore the implications of heterogeneous vaccination for the burden of CRS, introducing realistic vaccination scenarios based on reported levels of measles vaccine coverage. Our results highlight the potential impact of country-wide reductions of incidence of rubella on the local CRS burdens in districts with small population sizes. However, simulations indicate that if rubella vaccination is introduced with coverage reflecting current estimates for measles coverage in South Africa, the burden of CRS is likely to be reduced overall over a 30 year time horizon by a factor of 3, despite the fact that this coverage is lower than the traditional 80 per cent rule of thumb for vaccine introduction, probably owing to a combination of relatively low birth and transmission rates. We conclude by discussing the likely impact of private-sector vaccination.     

Accuracy of the Michaelis–Menten approximation when analysing effects of molecular noise

Quantitative biology relies on the construction of accurate mathematical models, yet the effectiveness of these models is often predicated on making simplifying approximations that allow for direct comparisons with available experimental data. The Michaelis–Menten (MM) approximation is widely used in both deterministic and discrete stochastic models of intracellular reaction networks, owing to the ubiquity of enzymatic activity in cellular processes and the clear biochemical interpretation of its parameters. However, it is not well understood how the approximation applies to the discrete stochastic case or how it extends to spatially inhomogeneous systems. We study the behaviour of the discrete stochastic MM approximation as a function of system size and show that significant errors can occur for small volumes, in comparison with a corresponding mass-action system. We then explore some consequences of these results for quantitative modelling. One consequence is that fluctuation-induced sensitivity, or stochastic focusing, can become highly exaggerated in models that make use of MM kinetics even if the approximations are excellent in a deterministic model. Another consequence is that spatial stochastic simulations based on the reaction–diffusion master equation can become highly inaccurate if the model contains MM terms.     

A dynamic dose–response model to account for exposure patterns in risk assessment: a case study in inhalation anthrax

The most commonly used dose–response models implicitly assume that accumulation of dose is a time-independent process where each pathogen has a fixed risk of initiating infection. Immune particle neutralization of pathogens, however, may create strong time dependence; i.e. temporally clustered pathogens have a better chance of overwhelming the immune particles than pathogen exposures that occur at lower levels for longer periods of time. In environmental transmission systems, we expect different routes of transmission to elicit different dose–timing patterns and thus potentially different realizations of risk. We present a dose–response model that captures time dependence in a manner that incorporates the dynamics of initial immune response. We then demonstrate the parameter estimation of our model in a dose–response survival analysis using empirical time-series data of inhalational anthrax in monkeys in which we find slight dose–timing effects. Future dose–response experiments should include varying the time pattern of exposure in addition to varying the total doses delivered. Ultimately, the dynamic dose–response paradigm presented here will improve modelling of environmental transmission systems where different systems have different time patterns of exposure.     

The interaction of seasonal forcing and immunity and the resonance dynamics of malaria

Theory has emphasized the importance of both intrinsic factors such as host immunity and extrinsic drivers such as climate in determining disease dynamics. In particular, seasonality may lead to multi-annual cycles in prevalence, but the likelihood of this depends on the role of acquired immunity. Some diseases including malaria have immunity that falls between the classic susceptible–infectious–removed and susceptible–infectious–susceptible models. Here, we investigate the general conditions promoting the subharmonic resonance behaviour that may lead to multi-annual cycles in a general malaria dynamical model. Utilizing two complementary approaches to bifurcation analyses, we show that resonance is promoted by processes shortening the length of the infectious period and that subharmonic cycles are favoured in situations with strong seasonality in transmission but at intermediate levels of endemicity. We discuss the implications of our results for understanding prevalence patterns in long-term malaria datasets from Kenya that show multi-annual cycles and one from Thailand that does not and discuss the possible implications of treatment.     

Parameterizing state–space models for infectious disease dynamics by generalized profiling: measles in Ontario

Parameter estimation for infectious disease models is important for basic understanding (e.g. to identify major transmission pathways), for forecasting emerging epidemics, and for designing control measures. Differential equation models are often used, but statistical inference for differential equations suffers from numerical challenges and poor agreement between observational data and deterministic models. Accounting for these departures via stochastic model terms requires full specification of the probabilistic dynamics, and computationally demanding estimation methods. Here, we demonstrate the utility of an alternative approach, generalized profiling, which provides robustness to violations of a deterministic model without needing to specify a complete probabilistic model. We introduce novel means for estimating the robustness parameters and for statistical inference in this framework. The methods are applied to a model for pre-vaccination measles incidence in Ontario, and we demonstrate the statistical validity of our inference through extensive simulation. The results confirm that school term versus summer drives seasonality of transmission, but we find no effects of short school breaks and the estimated basic reproductive ratio ℛ₀ greatly exceeds previous estimates. The approach applies naturally to any system for which candidate differential equations are available, and avoids many challenges that have limited Monte Carlo inference for state–space models.     

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.     

Spatio-temporal waves and targeted vaccination in recurrent epidemic network models

The success of an infectious disease to invade a population is strongly controlled by the population''s specific connectivity structure. Here, a network model is presented as an aid in understanding the role of social behaviour and heterogeneous connectivity in determining the spatio-temporal patterns of disease dynamics. We explore the controversial origins of long-term recurrent oscillations believed to be characteristic of diseases that have a period of temporary immunity after infection. In particular, we focus on sexually transmitted diseases such as syphilis, where this controversy is currently under review. Although temporary immunity plays a key role, it is found that, in realistic small-world networks, the social and sexual behaviour of individuals also has a great influence in generating long-term cycles. The model generates circular waves of infection with unusual spatial dynamics that depend on focal areas that act as pacemakers in the population. Eradication of the disease can be efficiently achieved by eliminating the pacemakers with a targeted vaccination scheme. A simple difference equation model is derived, which captures the infection dynamics of the network model and gives insights into their origins and their eradication through vaccination. Illustrative videos may be found in the electronic supplementary material.     

Role of cell polarity dynamics and motility in pattern formation due to contact-dependent signalling

A key challenge in biology is to understand how spatio-temporal patterns and structures arise during the development of an organism. An initial aggregate of spatially uniform cells develops and forms the differentiated structures of a fully developed organism. On the one hand, contact-dependent cell–cell signalling is responsible for generating a large number of complex, self-organized, spatial patterns in the distribution of the signalling molecules. On the other hand, the motility of cells coupled with their polarity can independently lead to collective motion patterns that depend on mechanical parameters influencing tissue deformation, such as cellular elasticity, cell–cell adhesion and active forces generated by actin and myosin dynamics. Although modelling efforts have, thus far, treated cell motility and cell–cell signalling separately, experiments in recent years suggest that these processes could be tightly coupled. Hence, in this paper, we study how the dynamics of cell polarity and migration influence the spatiotemporal patterning of signalling molecules. Such signalling interactions can occur only between cells that are in physical contact, either directly at the junctions of adjacent cells or through cellular protrusional contacts. We present a vertex model which accounts for contact-dependent signalling between adjacent cells and between non-adjacent neighbours through long protrusional contacts that occur along the orientation of cell polarization. We observe a rich variety of spatiotemporal patterns of signalling molecules that is influenced by polarity dynamics of the cells, relative strengths of adjacent and non-adjacent signalling interactions, range of polarized interaction, signalling activation threshold, relative time scales of signalling and polarity orientation, and cell motility. Though our results are developed in the context of Delta–Notch signalling, they are sufficiently general and can be extended to other contact dependent morpho-mechanical dynamics.     

Predators indirectly control vector-borne disease: linking predator–prey and host–pathogen models

Pathogens transmitted by arthropod vectors are common in human populations, agricultural systems and natural communities. Transmission of these vector-borne pathogens depends on the population dynamics of the vector species as well as its interactions with other species within the community. In particular, predation may be sufficient to control pathogen prevalence indirectly via the vector. To examine the indirect effect of predators on vectored-pathogen dynamics, we developed a theoretical model that integrates predator–prey and host–pathogen theory. We used this model to determine whether predation can prevent pathogen persistence or alter the stability of host–pathogen dynamics. We found that, in the absence of predation, pathogen prevalence in the host increases with vector fecundity, whereas predation on the vector causes pathogen prevalence to decline, or even become extinct, with increasing vector fecundity. We also found that predation on a vector may drastically slow the initial spread of a pathogen. The predator can increase host abundance indirectly by reducing or eliminating infection in the host population. These results highlight the importance of studying interactions that, within the greater community, may alter our predictions when studying disease dynamics. From an applied perspective, these results also suggest situations where an introduced predator or the natural enemies of a vector may slow the rate of spread of an emerging vector-borne pathogen.