Disentangling the causes of mumps reemergence in the United States.

Over the past two decades, multiple countries with high vaccine coverage have experienced resurgent outbreaks of mumps. Worryingly, in these countries, a high proportion of cases have been among those who have completed the recommended vaccination schedule, raising alarm about the effectiveness of existing vaccines. Two putative mechanisms of vaccine failure have been proposed as driving observed trends: 1) gradual waning of vaccine-derived immunity (necessitating additional booster doses) and 2) the introduction of novel viral genotypes capable of evading vaccinal immunity. Focusing on the United States, we conduct statistical likelihood-based hypothesis testing using a mechanistic transmission model on age-structured epidemiological, demographic, and vaccine uptake time series data. We find that the data are most consistent with the waning hypothesis and estimate that 32.8% (32%, 33.5%) of individuals lose vaccine-derived immunity by age 18 y. Furthermore, we show using our transmission model how waning vaccine immunity reproduces qualitative and quantitatively consistent features of epidemiological data, namely 1) the shift in mumps incidence toward older individuals, 2) the recent recurrence of mumps outbreaks, and 3) the high proportion of mumps cases among previously vaccinated individuals.

The Number and Pattern of Viral Genomic Reassortments are not Necessarily Identifiable from Segment Trees.

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Since large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the "first-infection tree," a potentially counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our observation of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment, but also the four additional situations that can complicate the inference of reassortment.

Markov genealogy processes.

We construct a family of genealogy-valued Markov processes that are induced by a continuous-time Markov population process. We derive exact expressions for the likelihood of a given genealogy conditional on the history of the underlying population process. These lead to a nonlinear filtering equation which can be used to design efficient Monte Carlo inference algorithms. We demonstrate these calculations with several examples. Existing full-information approaches for phylodynamic inference are special cases of the theory.

The contribution of host cell-directed vs. parasite-directed immunity to the disease and dynamics of malaria infections.

Hosts defend themselves against pathogens by mounting an immune response. Fully understanding the immune response as a driver of host disease and pathogen evolution requires a quantitative account of its impact on parasite population dynamics. Here, we use a data-driven modeling approach to quantify the birth and death processes underlying the dynamics of infections of the rodent malaria parasite, Plasmodium chabaudi, and the red blood cells (RBCs) it targets. We decompose the immune response into 3 components, each with a distinct effect on parasite and RBC vital rates, and quantify the relative contribution of each component to host disease and parasite density. Our analysis suggests that these components are deployed in a coordinated fashion to realize distinct resource-directed defense strategies that complement the killing of parasitized cells. Early in the infection, the host deploys a strategy reminiscent of siege and scorched-earth tactics, in which it both destroys RBCs and restricts their supply. Late in the infection, a "juvenilization" strategy, in which turnover of RBCs is accelerated, allows the host to recover from anemia while holding parasite proliferation at bay. By quantifying the impact of immunity on both parasite fitness and host disease, we reveal that phenomena often interpreted as immunopathology may in fact be beneficial to the host. Finally, we show that, across mice, the components of the host response are consistently related to each other, even when infections take qualitatively different trajectories. This suggests the existence of simple rules that govern the immune system's deployment.

The impact of infection-derived immunity on disease dynamics.

When modeling infectious diseases, it is common to assume that infection-derived immunity is either (1) non-existent or (2) perfect and lifelong. However there are many diseases in which infection-derived immunity is known to be present but imperfect. There are various ways in which infection-derived immunity can fail, which can ultimately impact the probability that an individual be reinfected by the same pathogen, as well as the long-run population-level prevalence of the pathogen. Here we discuss seven different models of imperfect infection-derived immunity, including waning, leaky and all-or-nothing immunity. For each model we derive the probability that an infected individual becomes reinfected during their lifetime, given that the system is at endemic equilibrium. This can be thought of as the impact that each of these infection-derived immunity failures have on reinfection. This measure is useful because it provides us with a way to compare different modes of failure of infection-derived immunity.

Differential and enhanced response to climate forcing in diarrheal disease due to rotavirus across a megacity of the developing world.

The role of climate forcing in the population dynamics of infectious diseases has typically been revealed via retrospective analyses of incidence records aggregated across space and, in particular, over whole cities. Here, we focus on the transmission dynamics of rotavirus, the main diarrheal disease in infants and young children, within the megacity of Dhaka, Bangladesh. We identify two zones, the densely urbanized core and the more rural periphery, that respond differentially to flooding. Moreover, disease seasonality differs substantially between these regions, spanning variation comparable to the variation from tropical to temperate regions. By combining process-based models with an extensive disease surveillance record, we show that the response to climate forcing is mainly seasonal in the core, where a more endemic transmission resulting from an asymptomatic reservoir facilitates the response to the monsoons. The force of infection in this monsoon peak can be an order of magnitude larger than the force of infection in the more epidemic periphery, which exhibits little or no postmonsoon outbreak in a pattern typical of nearby rural areas. A typically smaller peak during the monsoon season nevertheless shows sensitivity to interannual variability in flooding. High human density in the core is one explanation for enhanced transmission during troughs and an associated seasonal monsoon response in this diarrheal disease, which unlike cholera, has not been widely viewed as climate-sensitive. Spatial demographic, socioeconomic, and environmental heterogeneity can create reservoirs of infection and enhance the sensitivity of disease systems to climate forcing, especially in the populated cities of the developing world.

Unraveling the Transmission Ecology of Polio.

Sustained and coordinated vaccination efforts have brought polio eradication within reach. Anticipating the eradication of wild poliovirus (WPV) and the subsequent challenges in preventing its re-emergence, we look to the past to identify why polio rose to epidemic levels in the mid-20th century, and how WPV persisted over large geographic scales. We analyzed an extensive epidemiological dataset, spanning the 1930s to the 1950s and spatially replicated across each state in the United States, to glean insight into the drivers of polio's historical expansion and the ecological mode of its persistence prior to vaccine introduction. We document a latitudinal gradient in polio's seasonality. Additionally, we fitted and validated mechanistic transmission models to data from each US state independently. The fitted models revealed that: (1) polio persistence was the product of a dynamic mosaic of source and sink populations; (2) geographic heterogeneity of seasonal transmission conditions account for the latitudinal structure of polio epidemics; (3) contrary to the prevailing "disease of development" hypothesis, our analyses demonstrate that polio's historical expansion was straightforwardly explained by demographic trends rather than improvements in sanitation and hygiene; and (4) the absence of clinical disease is not a reliable indicator of polio transmission, because widespread polio transmission was likely in the multiyear absence of clinical disease. As the world edges closer to global polio eradication and continues the strategic withdrawal of the Oral Polio Vaccine (OPV), the regular identification of, and rapid response to, these silent chains of transmission is of the utmost importance.

Inference for dynamic and latent variable models via iterated, perturbed Bayes maps.

Iterated filtering algorithms are stochastic optimization procedures for latent variable models that recursively combine parameter perturbations with latent variable reconstruction. Previously, theoretical support for these algorithms has been based on the use of conditional moments of perturbed parameters to approximate derivatives of the log likelihood function. Here, a theoretical approach is introduced based on the convergence of an iterated Bayes map. An algorithm supported by this theory displays substantial numerical improvement on the computational challenge of inferring parameters of a partially observed Markov process.

Time-varying, serotype-specific force of infection of dengue virus.

Infectious disease models play a key role in public health planning. These models rely on accurate estimates of key transmission parameters such as the force of infection (FoI), which is the per-capita risk of a susceptible person being infected. The FoI captures the fundamental dynamics of transmission and is crucial for gauging control efforts, such as identifying vaccination targets. Dengue virus (DENV) is a mosquito-borne, multiserotype pathogen that currently infects ∼390 million people a year. Existing estimates of the DENV FoI are inaccurate because they rely on the unrealistic assumption that risk is constant over time. Dengue models are thus unreliable for designing vaccine deployment strategies. Here, we present to our knowledge the first time-varying (daily), serotype-specific estimates of DENV FoIs using a spline-based fitting procedure designed to examine a 12-y, longitudinal DENV serological dataset from Iquitos, Peru (11,703 individuals, 38,416 samples, and 22,301 serotype-specific DENV infections from 1999 to 2010). The yearly DENV FoI varied markedly across time and serotypes (0-0.33), as did daily basic reproductive numbers (0.49-4.72). During specific time periods, the FoI fluctuations correlated across serotypes, indicating that different DENV serotypes shared common transmission drivers. The marked variation in transmission intensity that we detected indicates that intervention targets based on one-time estimates of the FoI could underestimate the level of effort needed to prevent disease. Our description of dengue virus transmission dynamics is unprecedented in detail, providing a basis for understanding the persistence of this rapidly emerging pathogen and improving disease prevention programs.

Maternal pertussis immunisation: clinical gains and epidemiological legacy.

The increase in whooping cough (pertussis) incidence in many countries with high routine vaccination coverage is alarming, with incidence in the US reaching almost 50,000 reported cases per year, reflecting incidence levels not seen since the 1950s. While the potential explanations for this resurgence remain debated, we face an urgent need to protect newborns, especially during the time window between birth and the first routine vaccination dose. Maternal immunisation has been proposed as an effective strategy for protecting neonates, who are at higher risk of severe pertussis disease and mortality. However, if maternally derived antibodies adversely affect the immunogenicity of the routine schedule, through blunting effects, we may observe a gradual degradation of herd immunity. 'Wasted' vaccines would result in an accumulation of susceptible children in the population, specifically leading to an overall increase in incidence in older age groups. In this Perspective, we discuss potential long-term epidemiological effects of maternal immunisation, as determined by possible immune interference outcomes.

Climate-driven endemic cholera is modulated by human mobility in a megacity.

Although a differential sensitivity of cholera dynamics to climate variability has been reported in the spatially heterogeneous megacity of Dhaka, Bangladesh, the specific patterns of spread of the resulting risk within the city remain unclear. We build on an established probabilistic spatial model to investigate the importance and role of human mobility in modulating spatial cholera transmission. Mobility fluxes were inferred using a straightforward and generalizable methodology that relies on mapping population density based on a high resolution urban footprint product, and a parameter-free human mobility model. In accordance with previous findings, we highlight the higher sensitivity to the El Niño Southern Oscillation (ENSO) in the highly populated urban center than in the more rural periphery. More significantly, our results show that cholera risk is largely transmitted from the climate-sensitive core to the periphery of the city, with implications for the planning of control efforts. In addition, including human mobility improves the outbreak prediction performance of the model with an 11 month lead. The interplay between climatic and human mobility factors in cholera transmission is discussed from the perspective of the rapid growth of megacities across the developing world.

The pertussis enigma: reconciling epidemiology, immunology and evolution.

Pertussis, a highly contagious respiratory infection, remains a public health priority despite the availability of vaccines for 70 years. Still a leading cause of mortality in developing countries, pertussis has re-emerged in several developed countries with high vaccination coverage. Resurgence of pertussis in these countries has routinely been attributed to increased awareness of the disease, imperfect vaccinal protection or high infection rates in adults. In this review, we first present 1980-2012 incidence data from 63 countries and show that pertussis resurgence is not universal. We further argue that the large geographical variation in trends probably precludes a simple explanation, such as the transition from whole-cell to acellular pertussis vaccines. Reviewing available evidence, we then propose that prevailing views on pertussis epidemiology are inconsistent with both historical and contemporary data. Indeed, we summarize epidemiological evidence showing that natural infection and vaccination both appear to provide long-term protection against transmission and disease, so that previously infected or vaccinated adults contribute little to overall transmission at a population level. Finally, we identify several promising avenues that may lead to a consistent explanation of global pertussis epidemiology and to more effective control strategies.

Detecting Adaptive Evolution in Phylogenetic Comparative Analysis Using the Ornstein-Uhlenbeck Model.

Phylogenetic comparative analysis is an approach to inferring evolutionary process from a combination of phylogenetic and phenotypic data. The last few years have seen increasingly sophisticated models employed in the evaluation of more and more detailed evolutionary hypotheses, including adaptive hypotheses with multiple selective optima and hypotheses with rate variation within and across lineages. The statistical performance of these sophisticated models has received relatively little systematic attention, however. We conducted an extensive simulation study to quantify the statistical properties of a class of models toward the simpler end of the spectrum that model phenotypic evolution using Ornstein-Uhlenbeck processes. We focused on identifying where, how, and why these methods break down so that users can apply them with greater understanding of their strengths and weaknesses. Our analysis identifies three key determinants of performance: a discriminability ratio, a signal-to-noise ratio, and the number of taxa sampled. Interestingly, we find that model-selection power can be high even in regions that were previously thought to be difficult, such as when tree size is small. On the other hand, we find that model parameters are in many circumstances difficult to estimate accurately, indicating a relative paucity of information in the data relative to these parameters. Nevertheless, we note that accurate model selection is often possible when parameters are only weakly identified. Our results have implications for more sophisticated methods inasmuch as the latter are generalizations of the case we study.

Avoidable errors in the modelling of outbreaks of emerging pathogens, with special reference to Ebola.

As an emergent infectious disease outbreak unfolds, public health response is reliant on information on key epidemiological quantities, such as transmission potential and serial interval. Increasingly, transmission models fit to incidence data are used to estimate these parameters and guide policy. Some widely used modelling practices lead to potentially large errors in parameter estimates and, consequently, errors in model-based forecasts. Even more worryingly, in such situations, confidence in parameter estimates and forecasts can itself be far overestimated, leading to the potential for large errors that mask their own presence. Fortunately, straightforward and computationally inexpensive alternatives exist that avoid these problems. Here, we first use a simulation study to demonstrate potential pitfalls of the standard practice of fitting deterministic models to cumulative incidence data. Next, we demonstrate an alternative based on stochastic models fit to raw data from an early phase of 2014 West Africa Ebola virus disease outbreak. We show not only that bias is thereby reduced, but that uncertainty in estimates and forecasts is better quantified and that, critically, lack of model fit is more readily diagnosed. We conclude with a short list of principles to guide the modelling response to future infectious disease outbreaks.

Highly localized sensitivity to climate forcing drives endemic cholera in a megacity.

The population dynamics of endemic cholera in urban environments--in particular interannual variation in the size and distribution of seasonal outbreaks--remain poorly understood and highly unpredictable. In part, this situation is due to the considerable demographic, socioeconomic, and environmental heterogeneity of large and growing urban centers. Despite this heterogeneity, the influence of climate variability on the population dynamics of infectious diseases is considered a large-scale, regional, phenomenon, and as such has been previously addressed for cholera only with temporal models that do not incorporate spatial structure. Here we show that a probabilistic spatial model can explain cholera dynamics in the megacity of Dhaka, Bangladesh, and afford a basis for cholera forecasts at lead times of 11 mo. Critically, we find that the action of climate variability (El Niño southern oscillation and flooding) is quite localized: There is a climate-sensitive urban core that acts to propagate risk to the rest of the city. The modeling framework presented here should be applicable to cholera in other cities, as well as to other infectious diseases in urban settings and other biological systems with spatiotemporal interactions.

Human birth seasonality: latitudinal gradient and interplay with childhood disease dynamics.

More than a century of ecological studies have demonstrated the importance of demography in shaping spatial and temporal variation in population dynamics. Surprisingly, the impact of seasonal recruitment on infectious disease systems has received much less attention. Here, we present data encompassing 78 years of monthly natality in the USA, and reveal pronounced seasonality in birth rates, with geographical and temporal variation in both the peak birth timing and amplitude. The timing of annual birth pulses followed a latitudinal gradient, with northern states exhibiting spring/summer peaks and southern states exhibiting autumn peaks, a pattern we also observed throughout the Northern Hemisphere. Additionally, the amplitude of United States birth seasonality was more than twofold greater in southern states versus those in the north. Next, we examined the dynamical impact of birth seasonality on childhood disease incidence, using a mechanistic model of measles. Birth seasonality was found to have the potential to alter the magnitude and periodicity of epidemics, with the effect dependent on both birth peak timing and amplitude. In a simulation study, we fitted an susceptible-exposed-infected-recovered model to simulated data, and demonstrated that ignoring birth seasonality can bias the estimation of critical epidemiological parameters. Finally, we carried out statistical inference using historical measles incidence data from New York City. Our analyses did not identify the predicted systematic biases in parameter estimates. This may be owing to the well-known frequency-locking between measles epidemics and seasonal transmission rates, or may arise from substantial uncertainty in multiple model parameters and estimation stochasticity.

Inapparent infections and cholera dynamics.

In many infectious diseases, an unknown fraction of infections produce symptoms mild enough to go unrecorded, a fact that can seriously compromise the interpretation of epidemiological records. This is true for cholera, a pandemic bacterial disease, where estimates of the ratio of asymptomatic to symptomatic infections have ranged from 3 to 100 (refs 1-5). In the absence of direct evidence, understanding of fundamental aspects of cholera transmission, immunology and control has been based on assumptions about this ratio and about the immunological consequences of inapparent infections. Here we show that a model incorporating high asymptomatic ratio and rapidly waning immunity, with infection both from human and environmental sources, explains 50 yr of mortality data from 26 districts of Bengal, the pathogen's endemic home. We find that the asymptomatic ratio in cholera is far higher than had been previously supposed and that the immunity derived from mild infections wanes much more rapidly than earlier analyses have indicated. We find, too, that the environmental reservoir (free-living pathogen) is directly responsible for relatively few infections but that it may be critical to the disease's endemicity. Our results demonstrate that inapparent infections can hold the key to interpreting the patterns of disease outbreaks. New statistical methods, which allow rigorous maximum likelihood inference based on dynamical models incorporating multiple sources and outcomes of infection, seasonality, process noise, hidden variables and measurement error, make it possible to test more precise hypotheses and obtain unexpected results. Our experience suggests that the confrontation of time-series data with mechanistic models is likely to revise our understanding of the ecology of many infectious diseases.

Natural immune boosting in pertussis dynamics and the potential for long-term vaccine failure.

Incidence of whooping cough, unlike many other childhood diseases for which there is an efficacious vaccine, has been increasing over the past twenty years despite high levels of vaccine coverage. Its reemergence has been particularly noticeable among teenagers and adults. Many hypotheses have been put forward to explain these two patterns, but parsimonious reconciliation of clinical data on the limited duration of immunity with both pre- and postvaccine era age-specific incidence remains a challenge. We consider the immunologically relevant, yet epidemiologically largely neglected, possibility that a primed immune system can respond to a lower dose of antigen than a naive one. We hypothesize that during the prevaccine era teenagers' and adults' primed immunity was frequently boosted by reexposure, so maintaining herd immunity in the face of potentially eroding individual immunity. In contrast, low pathogen circulation in the current era, except during epidemic outbreaks, allows immunity to be lost before reexposure occurs. We develop and analyze an age-structured model that encapsulates this hypothesis. We find that immune boosting must be more easily triggered than primary infection to account for age-incidence data. We make age-specific and dynamical predictions through bifurcation analysis and simulation. The boosting model proposed here parsimoniously captures four key features of pertussis data from highly vaccinated countries: (i) the shift in age-specific incidence, (ii) reemergence with high vaccine coverage, (iii) the possibility for cyclic dynamics in the pre- and postvaccine eras, and (iv) the apparent shift from susceptible-infectious-recovered (SIR)-like to susceptible-infectious-recovered-susceptible (SIRS)-like phenomenology of infection and immunity to Bordetella pertussis.

EXACT PHYLODYNAMIC LIKELIHOOD VIA STRUCTURED MARKOV GENEALOGY PROCESSES.

We consider genealogies arising from a Markov population process in which individuals are categorized into a discrete collection of compartments, with the requirement that individuals within the same compartment are statistically exchangeable. When equipped with a sampling process, each such population process induces a time-evolving tree-valued process defined as the genealogy of all sampled individuals. We provide a construction of this genealogy process and derive exact expressions for the likelihood of an observed genealogy in terms of filter equations. These filter equations can be numerically solved using standard Monte Carlo integration methods. Thus, we obtain statistically efficient likelihood-based inference for essentially arbitrary compartment models based on an observed genealogy of individuals sampled from the population.

Red blood cell dynamics during malaria infection violate the assumptions of mathematical models of infection dynamics.

For decades, mathematical models have been used to understand the course and outcome of malaria infections (i.e., infection dynamics) and the evolutionary dynamics of the parasites that cause them. A key conclusion of these models is that red blood cell (RBC) availability is a fundamental driver of infection dynamics and parasite trait evolution. The extent to which this conclusion holds will in part depend on model assumptions about the host-mediated processes that regulate RBC availability i.e., removal of uninfected RBCs and supply of RBCs. Diverse mathematical functions have been used to describe host-mediated RBC supply and clearance, but it remains unclear whether they adequately capture the dynamics of RBC supply and clearance during infection. Here, we use a unique dataset, comprising time-series measurements of erythrocyte (i.e., mature RBC) and reticulocyte (i.e., newly supplied RBC) densities during Plasmodium chabaudi malaria infection, and a quantitative data-transformation scheme to elucidate whether RBC dynamics conform to common model assumptions. We found that RBC clearance and supply are not well described by mathematical functions commonly used to model these processes. Furthermore, the temporal dynamics of both processes vary with parasite growth rate in a manner again not captured by existing models. Together, these finding suggest that new model formulations are required if we are to explain and ultimately predict the within-host population dynamics and evolution of malaria parasites.