Search

To support decision-making and policy for managing epidemics of emerging pathogens, we present a model for inference and scenario analysis of SARS-CoV-2 transmission in the USA. The stochastic SEIR-type model includes compartments for latent, asymptomatic, detected and undetected symptomatic individuals, and hospitalized cases, and features realistic interval distributions for presymptomatic and symptomatic periods, time varying rates of case detection, diagnosis, and mortality. The model accounts for the effects on transmission of human mobility using anonymized mobility data collected from cellular devices, and of difficult to quantify environmental and behavioral factors using a latent process. The baseline transmission rate is the product of a human mobility metric obtained from data and this fitted latent process. We fit the model to incident case and death reports for each state in the USA and Washington D.C., using likelihood Maximization by Iterated particle Filtering (MIF). Observations (daily case and death reports) are modeled as arising from a negative binomial reporting process. We estimate time-varying transmission rate, parameters of a sigmoidal time-varying fraction of hospitalized cases that result in death, extra-demographic process noise, two dispersion parameters of the observation process, and the initial sizes of the latent, asymptomatic, and symptomatic classes. In a retrospective analysis covering March-December 2020, we show how mobility and transmission strength became decoupled across two distinct phases of the pandemic. The decoupling demonstrates the need for flexible, semi-parametric approaches for modeling infectious disease dynamics in real-time.

Directly observed therapy (DOT) is almost universally used for the treatment of TB. Several meta-analyses using different methods have assessed the effectiveness of DOT compared to self-administered therapy (SAT). The results of these meta-analyses often conflict with some concluding DOT is superior and others that there is little or no difference. Meta-analyses can guide policymaking, but such analyses must be reliable. To assess the validity of a previous meta-analysis, we tried to reproduce it. We encountered problems with the previous analysis that did not allow for a meaningful reproduction. We describe the issues we encountered here. We then performed a new meta-analysis comparing the treatment outcomes of adults given treatment with SAT versus DOT. Outcomes in the new analysis are loss to follow-up, treatment failure, cure, treatment completed, and all-cause mortality. All data, documentation, and code used to generate our results is provided. Our new analysis included four randomized and three observational studies with 1603 and 1626 individuals respectively. The pooled relative risks (RR) are as follows: Lost to follow-up (RR = 1.2, 95% CI 0.9, 1.7), Treatment Failure (RR = 1.1, 95% CI 0.6, 2), Cure (RR = 0.9, 95% CI 0.8, 1.1), Treatment Completion (RR = 1, 95% CI 0.9, 1.1), Mortality (RR = 0.9, 95% CI 0.6, 1.3). Based on data from our new meta-analysis, the magnitude of the difference between DOT and SAT for all reported outcomes is small, and none of the differences are statistically significant.

[This corrects the article DOI: 10.1371/journal.pone.0150568.].

Vaccination is an effective method to protect against infectious diseases. An important consideration in any vaccine formulation is the inoculum dose, i.e., amount of antigen or live attenuated pathogen that is used. Higher levels generally lead to better stimulation of the immune response but might cause more severe side effects and allow for less population coverage in the presence of vaccine shortages. Determining the optimal amount of inoculum dose is an important component of rational vaccine design. A combination of mathematical models with experimental data can help determine the impact of the inoculum dose. We illustrate the concept of using data and models to inform inoculum dose determination for vaccines, wby fitting a mathematical model to data from influenza A virus (IAV) infection of mice and human parainfluenza virus (HPIV) infection of cotton rats at different inoculum doses. We use the model to map inoculum dose to the level of immune protection and morbidity and to explore how such a framework might be used to determine an optimal inoculum dose. We show how a framework that combines mathematical models with experimental data can be used to study the impact of inoculum dose on important outcomes such as immune protection and morbidity. Our findings illustrate that the impact of inoculum dose on immune protection and morbidity can depend on the specific pathogen and that both protection and morbidity do not necessarily increase monotonically with increasing inoculum dose. Once vaccine design goals are specified with required levels of protection and acceptable levels of morbidity, our proposed framework can help in the rational design of vaccines and determination of the optimal amount of inoculum.

Determining the duration of protective immunity requires quantifying the magnitude and rate of loss of antibodies to different virus and vaccine antigens. A key complication is heterogeneity in both the magnitude and decay rate of responses of different individuals to a given vaccine, as well as of a given individual to different vaccines. We analyzed longitudinal data on antibody titers in 45 individuals to characterize the extent of this heterogeneity and used models to determine how it affected the longevity of protective immunity to measles, rubella, vaccinia, tetanus, and diphtheria. Our analysis showed that the magnitude of responses in different individuals varied between 12- and 200-fold (95% coverage) depending on the antigen. Heterogeneity in the magnitude and decay rate contribute comparably to variation in the longevity of protective immunity between different individuals. We found that some individuals have, on average, slightly longer-lasting memory than others-on average, they have higher antibody levels with slower decay rates. We identified different patterns for the loss of protective levels of antibodies to different vaccine and virus antigens. Specifically, we found that for the first 25 to 50 years, virtually all individuals have protective antibody titers against diphtheria and tetanus, respectively, but about 10% of the population subsequently lose protective immunity per decade. In contrast, at the outset, not all individuals had protective titers against measles, rubella, and vaccinia. However, these antibody titers wane much more slowly, with a loss of protective immunity in only 1% to 3% of the population per decade. Our results highlight the importance of long-term longitudinal studies for estimating the duration of protective immunity and suggest both how vaccines might be improved and how boosting schedules might be reevaluated.

For most pathogens, vaccination reduces the spread of the infection and total number of cases; thus, public policy usually advocates maximizing vaccination coverage. We use simple mathematical models to explore how this may be different for pathogens, such as influenza, which exhibit strain variation. Our models predict that the total number of seasonal influenza infections is minimized at an intermediate (rather than maximal) level of vaccination, and, somewhat counter-intuitively, further increasing the level of the vaccination coverage may lead to higher number of influenza infections and be detrimental to the public interest. This arises due to the combined effects of: competition between multiple co-circulating strains; limited breadth of protection afforded by the vaccine; and short-term strain-transcending immunity following natural infection. The study highlights the need for better quantification of the components of vaccine efficacy and longevity of strain-transcending cross-immunity in order to generate nuanced recommendations for influenza vaccine coverage levels.

The World Health Organization identifies influenza as a major public health problem. While the strains commonly circulating in humans usually do not cause severe pathogenicity in healthy adults, some strains that have infected humans, such as H5N1, can cause high morbidity and mortality. Based on the severity of the disease, influenza viruses are sometimes categorized as either being highly pathogenic (HP) or having low pathogenicity (LP). The reasons why some strains are LP and others HP are not fully understood. While there are likely multiple mechanisms of interaction between the virus and the immune response that determine LP versus HP outcomes, we focus here on one component, namely macrophages (MP). There is some evidence that MP may both help fight the infection and become productively infected with HP influenza viruses. We developed mathematical models for influenza infections which explicitly included the dynamics and action of MP. We fit these models to viral load and macrophage count data from experimental infections of mice with LP and HP strains. Our results suggest that MP may not only help fight an influenza infection but may contribute to virus production in infections with HP viruses. We also explored the impact of combination therapies with antivirals and anti-inflammatory drugs on HP infections. Our study suggests a possible mechanism of MP in determining HP versus LP outcomes, and how different interventions might affect infection dynamics.

CD8(+) T cells directed against conserved viral regions elicit broad immunity against distinct influenza viruses, promote rapid virus elimination and enhanced host recovery. The influenza neuraminidase inhibitor, oseltamivir, is prescribed for therapy and prophylaxis, although it remains unclear how the drug impacts disease severity and establishment of effector and memory CD8(+) T cell immunity. We dissected the effects of oseltamivir on viral replication, inflammation, acute CD8(+) T cell responses and the establishment of immunological CD8(+) T cell memory. In mice, ferrets and humans, the effect of osteltamivir on viral titre was relatively modest. However, prophylactic oseltamivir treatment in mice markedly reduced morbidity, innate responses, inflammation and, ultimately, the magnitude of effector CD8(+) T cell responses. Importantly, functional memory CD8(+) T cells established during the drug-reduced effector phase were capable of mounting robust recall responses. Moreover, influenza-specific memory CD4(+) T cells could be also recalled after the secondary challenge, while the antibody levels were unaffected. This provides evidence that long-term memory T cells can be generated during an oseltamivir-interrupted infection. The anti-inflammatory effect of oseltamivir was verified in H1N1-infected patients. Thus, in the case of an unpredicted influenza pandemic, while prophylactic oseltamivir treatment can reduce disease severity, the capacity to generate memory CD8(+) T cells specific for the newly emerged virus is uncompromised. This could prove especially important for any new influenza pandemic which often occurs in separate waves.

Tuberculosis (TB) is an infectious disease with a peculiar feature: Upon infection with the causative agent, Mycobacterium Tuberculosis (MTB), most hosts enter a latent state during which no transmission of MTB to new hosts occurs. Only a fraction of latently infected hosts develop TB disease and can potentially infect new hosts. At first glance, this seems like a waste of transmission potential and therefore an evolutionary suboptimal strategy for MTB. It might be that the human immune response keeps MTB in check in most hosts, thereby preventing it from achieving its evolutionary optimum. Another possible explanation is that long latency and progression to disease in only a fraction of hosts are evolutionary beneficial to MTB by allowing it to persist better in small host populations. Given that MTB has co-evolved with human hosts for millenia or longer, it likely encountered small host populations for a large share of its evolutionary history and had to evolve strategies of persistence. Here, we use a mathematical model to show that indeed, MTB persistence is optimal for an intermediate duration of latency and level of activation. The predicted optimal level of activation is above the observed value, suggesting that human co-evolution has lead to host immunity, which keeps MTB below its evolutionary optimum.

Successful replication within an infected host and successful transmission between hosts are key to the continued spread of most pathogens. Competing selection pressures exerted at these different scales can lead to evolutionary trade-offs between the determinants of fitness within and between hosts. Here, we examine such a trade-off in the context of influenza A viruses and the differential pressures exerted by temperature-dependent virus persistence. For a panel of avian influenza A virus strains, we find evidence for a trade-off between the persistence at high versus low temperatures. Combining a within-host model of influenza infection dynamics with a between-host transmission model, we study how such a trade-off affects virus fitness on the host population level. We show that conclusions regarding overall fitness are affected by the type of link assumed between the within- and between-host levels and the main route of transmission (direct or environmental). The relative importance of virulence and immune response mediated virus clearance are also found to influence the fitness impacts of virus persistence at low versus high temperatures. Based on our results, we predict that if transmission occurs mainly directly and scales linearly with virus load, and virulence or immune responses are negligible, the evolutionary pressure for influenza viruses to evolve toward good persistence at high within-host temperatures dominates. For all other scenarios, influenza viruses with good environmental persistence at low temperatures seem to be favored.

If repeated interventions against multiple outbreaks are not feasible, there is an optimal level of control during the first outbreak. Any control measures above that optimal level will lead to an outcome that may be as sub-optimal as that achieved by an intervention that is too weak. We studied this scenario in more detail.An age-stratified ordinary-differential-equation model was constructed to study infectious disease outbreaks and control in a population made up of two groups, adults and children. The model was parameterized using influenza as an example. This model was used to simulate two consecutive outbreaks of the same infectious disease, with an intervention applied only during the first outbreak, and to study how cumulative attack rates were influenced by population composition, strength of inter-group transmission, and different ways of triggering and implementing the interventions. We assumed that recovered individuals are fully immune and the intervention does not confer immunity.The optimal intervention depended on coupling between the two population sub-groups, the length, strength and timing of the intervention, and the population composition. Population heterogeneity affected intervention strategies only for very low cross-transmission between groups. At more realistic values, coupling between the groups led to synchronization of outbreaks and therefore intervention strategies that were optimal in reducing the attack rates for each subgroup and the population overall coincided. For a sustained intervention of low efficacy, early intervention was found to be best, while at high efficacies, a delayed start was better. For short interventions, a delayed start was always advantageous, independent of the intervention efficacy. For most scenarios, starting the intervention after a certain cumulative proportion of children were infected seemed more robust in achieving close to optimal outcomes compared to a strategy that used a specified duration after an outbreak's beginning as the trigger.

Small populations are thought to be adaptively handicapped, not only because they suffer more from deleterious mutations but also because they have limited access to new beneficial mutations, particularly those conferring large benefits.Here, we test this widely held conjecture using both simulations and experiments with small and large bacterial populations evolving in either a simple or a complex nutrient environment. Consistent with expectations, we find that small populations are adaptively constrained in the simple environment; however, in the complex environment small populations not only follow more heterogeneous adaptive trajectories, but can also attain higher fitness than the large populations. Large populations are constrained to near deterministic fixation of rare large-benefit mutations. While such determinism speeds adaptation on the smooth adaptive landscape represented by the simple environment, it can limit the ability of large populations from effectively exploring the underlying topography of rugged adaptive landscapes characterized by complex environments.Our results show that adaptive constraints often faced by small populations can be circumvented during evolution on rugged adaptive landscapes.

Neuraminidase Inhibitors (NI) are currently the most effective drugs against influenza. Recent cases of NI resistance are a cause for concern. To assess the danger of NI resistance, a number of studies have reported the fraction of treated patients from which resistant strains could be isolated. Unfortunately, those results strongly depend on the details of the experimental protocol. Additionally, knowing the fraction of patients harboring resistance is not too useful by itself. Instead, we want to know how likely it is that an infected patient can generate a resistant infection in a secondary host, and how likely it is that the resistant strain subsequently spreads. While estimates for these parameters can often be obtained from epidemiological data, such data is lacking for NI resistance in influenza. Here, we use an approach that does not rely on epidemiological data. Instead, we combine data from influenza infections of human volunteers with a mathematical framework that allows estimation of the parameters that govern the initial generation and subsequent spread of resistance. We show how these parameters are influenced by changes in drug efficacy, timing of treatment, fitness of the resistant strain, and details of virus and immune system dynamics. Our study provides estimates for parameters that can be directly used in mathematical and computational models to study how NI usage might lead to the emergence and spread of resistance in the population. We find that the initial generation of resistant cases is most likely lower than the fraction of resistant cases reported. However, we also show that the results depend strongly on the details of the within-host dynamics of influenza infections, and most importantly, the role the immune system plays. Better knowledge of the quantitative dynamics of the immune response during influenza infections will be crucial to further improve the results.

Pathogens that evolve resistance to drugs usually have reduced fitness. However, mutations that largely compensate for this reduction in fitness often arise. We investigate how these compensatory mutations affect population-wide resistance emergence as a function of drug treatment. Using a model of gonorrhea transmission dynamics, we obtain generally applicable, qualitative results that show how compensatory mutations lead to more likely and faster resistance emergence. We further show that resistance emergence depends on the level of drug use in a strongly nonlinear fashion. We also discuss what data need to be obtained to allow future quantitative predictions of resistance emergence.