Search

tsiR is an open source software package implemented in the R programming language designed to analyze infectious disease time-series data. The software extends a well-studied and widely-applied algorithm, the time-series Susceptible-Infected-Recovered (TSIR) model, to infer parameters from incidence data, such as contact seasonality, and to forward simulate the underlying mechanistic model. The tsiR package aggregates a number of different fitting features previously described in the literature in a user-friendly way, providing support for their broader adoption in infectious disease research. Also included in tsiR are a number of diagnostic tools to assess the fit of the TSIR model. This package should be useful for researchers analyzing incidence data for fully-immunizing infectious diseases.

We analyze the impact of birth seasonality (seasonal oscillations in the birth rate) on the dynamics of acute, immunizing childhood infectious diseases. Previous research has explored the effect of human birth seasonality on infectious disease dynamics using parameters appropriate for the developed world. We build on this work by including in our analysis an extended range of baseline birth rates and amplitudes, which correspond to developing world settings. Additionally, our analysis accounts for seasonal forcing both in births and contact rates. We focus in particular on the dynamics of measles. In the absence of seasonal transmission rates or stochastic forcing, for typical measles epidemiological parameters, birth seasonality induces either annual or biennial epidemics. Changes in the magnitude of the birth fluctuations (birth amplitude) can induce significant changes in the size of the epidemic peaks, but have little impact on timing of disease epidemics within the year. In contrast, changes to the birth seasonality phase (location of the peak in birth amplitude within the year) significantly influence the timing of the epidemics. In the presence of seasonality in contact rates, at relatively low birth rates (20 per 1000), birth amplitude has little impact on the dynamics but does have an impact on the magnitude and timing of the epidemics. However, as the mean birth rate increases, both birth amplitude and phase play an important role in driving the dynamics of the epidemic. There are stronger effects at higher birth rates.

INTRODUCTION: Rotavirus vaccine efficacy ranges from >90% in high socio-economic settings (SES) to 50% in low SES. With the imminent introduction of rotavirus vaccine in low SES countries, understanding reasons for reduced efficacy in these settings could identify strategies to improve vaccine performance. METHODS: We developed a mathematical model to predict rotavirus vaccine efficacy in high, middle and low SES based on data specific for each setting on incidence, protection conferred by natural infection and immune response to vaccination. We then examined factors affecting efficacy. RESULTS: Vaccination was predicted to prevent 93%, 86% and 51% of severe rotavirus gastroenteritis in high, middle and low SES, respectively. Also predicted was that vaccines are most effective against severe disease and efficacy declines with age in low but not high SES. Reduced immunogenicity of vaccination and reduced protection conferred by natural infection are the main factors that compromise efficacy in low SES. DISCUSSION: The continued risk of severe disease in non-primary natural infections in low SES is a key factor underpinning reduced efficacy of rotavirus vaccines. Predicted efficacy was remarkably consistent with observed clinical trial results from different SES, validating the model. The phenomenon of reduced vaccine efficacy can be predicted by intrinsic immunological and epidemiological factors of low SES populations. Modifying aspects of the vaccine (e.g. improving immunogenicity in low SES) and vaccination program (e.g. additional doses) may bring improvements.

Influenza virus tissue tropism defines the host cells and tissues that support viral replication and contributes to determining which regions of the respiratory tract are infected in humans. The location of influenza virus infection along the respiratory tract is a key determinant of virus pathogenicity and transmissibility, which are at the basis of influenza burdens in the human population. As the pathogenicity and transmissibility of influenza virus ultimately determine its reproductive fitness at the population level, strong selective pressures will shape influenza virus tissue tropisms that maximize fitness. At present, the relationships between influenza virus tissue tropism within hosts and reproductive fitness at the population level are poorly understood. The selective pressures and constraints that shape tissue tropism and thereby influence the location of influenza virus infection along the respiratory tract are not well characterized. We use mathematical models that link within-host infection dynamics in a spatially-structured human respiratory tract to between-host transmission dynamics, with the aim of characterizing the possible selective pressures on influenza virus tissue tropism. The results indicate that spatial heterogeneities in virus clearance, virus pathogenicity or both, resulting from the unique structure of the respiratory tract, may drive optimal receptor binding affinity-that maximizes influenza virus reproductive fitness at the population level-towards sialic acids with alpha 2,6 linkage to galactose. The expanding cell pool deeper down the respiratory tract, in association with lower clearance rates, may result in optimal infectivity rates-that likewise maximize influenza virus reproductive fitness at the population level-to exhibit a decreasing trend towards deeper regions of the respiratory tract. Lastly, pre-existing immunity may drive influenza virus tissue tropism towards upper regions of the respiratory tract. The proposed framework provides a new template for the cross-scale study of influenza virus evolutionary and epidemiological dynamics in humans.

For vaccine-preventable infections, immunization generally needs to be supplemented by palliative care of individuals missed by the vaccination. Costs and availability of vaccine doses and palliative care vary by disease and by region. In many situations, resources for delivery of palliative care are independent of resources required for vaccination; however we also need to consider the conservative scenario where there is some trade-off between efforts, which is of potential relevance for resource-poor settings. We formulate an SEIR model that includes those two control strategies--vaccination and palliative care. We consider their relative merit and optimal allocation in the context of a highly efficacious vaccine, and under the assumption that palliative care may reduce transmission. We investigate the utility of a range of mixed or pure strategies that can be implemented after an epidemic has started, and look for rule-of-thumb principles of how best to reduce the burden of disease during an acute outbreak over a spectrum of vaccine-preventable infections. Intuitively, we expect the best strategy to initially focus on vaccination, and enhanced palliative care after the infection has peaked, but a number of plausible realistic constraints for control result in important qualifications on the intervention strategy. The time in the epidemic when one should switch strategy depends sensitively on the relative cost of vaccine to palliative care, the available budget, and R0. Crucially, outbreak response vaccination may be more effective in managing low-R0 diseases, while high R0 scenarios enhance the importance of routine vaccination and case management.

Outbreaks of avian influenza in poultry can be devastating, yet many of the basic epidemiological parameters have not been accurately characterised. In 1999-2000 in Northern Italy, outbreaks of H7N1 low pathogenicity avian influenza virus (LPAI) were followed by the emergence of H7N1 highly pathogenic avian influenza virus (HPAI). This study investigates the transmission dynamics in turkeys of representative HPAI and LPAI H7N1 virus strains from this outbreak in an experimental setting, allowing direct comparison of the two strains. The fitted transmission rates for the two strains are similar: 2.04 (1.5-2.7) per day for HPAI, 2.01 (1.6-2.5) per day for LPAI. However, the mean infectious period is far shorter for HPAI (1.47 (1.3-1.7) days) than for LPAI (7.65 (7.0-8.3) days), due to the rapid death of infected turkeys. Hence the basic reproductive ratio, [Formula: see text] is significantly lower for HPAI (3.01 (2.2-4.0)) than for LPAI (15.3 (11.8-19.7)). The comparison of transmission rates and [Formula: see text] are critically important in relation to understanding how HPAI might emerge from LPAI. Two competing hypotheses for how transmission rates vary with population size are tested by fitting competing models to experiments with differing numbers of turkeys. A model with frequency-dependent transmission gives a significantly better fit to experimental data than density-dependent transmission. This has important implications for extrapolating experimental results from relatively small numbers of birds to the commercial poultry flock size, and for how control, including vaccination, might scale with flock size.

Background Influenza viruses show a significant capacity to evade host immunity; this is manifest both as large occasional jumps in the antigenic phenotype of viral surface molecules and in gradual antigenic changes leading to annual influenza epidemics in humans. Recent mouse studies show that avidity for host cells can play an important role in polyclonal antibody escape, and further that electrostatic charge of the hemagglutinin glycoprotein can contribute to such avidity. Methodology/Principal Findings We test the role of glycoprotein charge on sequence data from the three major subtypes of influenza A in humans, using a simple method of calculating net glycoprotein charge. Of all subtypes, H3N2 in humans shows a striking pattern of increasing positive charge since its introduction in 1968. Notably, this trend applies to both hemagglutinin and neuraminidase glycoproteins. In the late 1980s hemagglutinin charge reached a plateau, while neuraminidase charge started to decline. We identify key groups of amino acid sites involved in this charge trend. Conclusions/Significance To our knowledge these are the first indications that, for human H3N2, net glycoprotein charge covaries strongly with antigenic drift on a global scale. Further work is needed to elucidate how such charge interacts with other immune escape mechanisms, such as glycosylation, and we discuss important questions arising for future study.

Metapopulation rescue effects are thought to be key to the persistence of many acute immunizing infections. Yet the enhancement of persistence through spatial coupling has not been previously quantified. Here we estimate the metapopulation rescue effects for four childhood infections using global WHO reported incidence data by comparing persistence on island countries vs all other countries, while controlling for key variables such as vaccine cover, birth rates and economic development. The relative risk of extinction on islands is significantly higher, and approximately double the risk of extinction in mainland countries. Furthermore, as may be expected, infections with longer infectious periods tend to have the strongest metapopulation rescue effects. Our results quantitate the notion that demography and local community size controls disease persistence.

Background: Norovirus (NoV) transmission may be impacted by changes in symptom intensity. Sudden onset of vomiting, which may cause an initial period of hyper-infectiousness, often marks the beginning of symptoms. This is often followed by: a 1–3 day period of milder symptoms, environmental contamination following vomiting, and post-symptomatic shedding that may result in transmission at progressively lower rates. Existing models have not included time-varying infectiousness, though representing these features could add utility to models of NoV transmission. Methods: We address this by comparing the fit of three models (Models 1–3) of NoV infection to household transmission data from a 2009 point-source outbreak of GII.12 norovirus in North Carolina. Model 1 is an SEIR compartmental model, modified to allow Gamma-distributed sojourn times in the latent and infectious classes, where symptomatic cases are uniformly infectious over time. Model 2 assumes infectiousness decays exponentially as a function of time since onset, while Model 3 is discontinuous, with a spike concentrating 50% of transmissibility at onset. We use Bayesian data augmentation techniques to estimate transmission parameters for each model, and compare their goodness of fit using qualitative and quantitative model comparison. We also assess the robustness of our findings to asymptomatic infections. Results: We find that Model 3 (initial spike in shedding) best explains the household transmission data, using both quantitative and qualitative model comparisons. We also show that these results are robust to the presence of asymptomatic infections. Conclusions: Explicitly representing explosive NoV infectiousness at onset should be considered when developing models and interventions to interrupt and prevent outbreaks of norovirus in the community. The methods presented here are generally applicable to the transmission of pathogens that exhibit large variation in transmissibility over an infection.

Early observations from countries that have introduced rotavirus vaccination suggest that there may be indirect protection for unvaccinated individuals, but it is unclear whether these benefits will extend to the long term. Transmission dynamic models have attempted to quantify the indirect protection that might be expected from rotavirus vaccination in developed countries, but results have varied. To better understand the magnitude and sources of variability in model projections, we undertook a comparative analysis of transmission dynamic models for rotavirus. We fit five models to reported rotavirus gastroenteritis (RVGE) data from England and Wales, and evaluated outcomes for short- and long-term vaccination effects. All of our models reproduced the important features of rotavirus epidemics in England and Wales. Models predicted that during the initial year after vaccine introduction, incidence of severe RVGE would be reduced 1.8-2.9 times more than expected from the direct effects of the vaccine alone (28-50% at 90% coverage), but over a 5-year period following vaccine introduction severe RVGE would be reduced only by 1.1-1.7 times more than expected from the direct effects (54-90% at 90% coverage). Projections for the long-term reduction of severe RVGE ranged from a 55% reduction at full coverage to elimination with at least 80% coverage. Our models predicted short-term reductions in the incidence of RVGE that exceeded estimates of the direct effects, consistent with observations from the United States and other countries. Some of the models predicted that the short-term indirect benefits may be offset by a partial shifting of the burden of RVGE to older unvaccinated individuals. Nonetheless, even when such a shift occurs, the overall reduction in severe RVGE is considerable. Discrepancies among model predictions reflect uncertainties about age variation in the risk and reporting of RVGE, and the duration of natural and vaccine-induced immunity, highlighting important questions for future research.