Infection-induced increases to population size during cycles in a discrete-time epidemic model.

One-dimensional discrete-time population models, such as those that involve Logistic or Ricker growth, can exhibit periodic and chaotic dynamics. Expanding the system by one dimension to incorporate epidemiological interactions causes an interesting complexity of new behaviors. Here, we examine a discrete-time two-dimensional susceptible-infectious (SI) model with Ricker growth and show that the introduction of infection can not only produce a distinctly different bifurcation structure than that of the underlying disease-free system but also lead to counter-intuitive increases in population size. We use numerical bifurcation analysis to determine the influence of infection on the location and types of bifurcations. In addition, we examine the appearance and extent of a phenomenon known as the 'hydra effect,' i.e., increases in total population size when factors, such as mortality, that act negatively on a population, are increased. Previous work, primarily focused on dynamics at fixed points, showed that the introduction of infection that reduces fecundity to the SI model can lead to a so-called 'infection-induced hydra effect.' Our work shows that even in such a simple two-dimensional SI model, the introduction of infection that alters fecundity or mortality can produce dynamics can lead to the appearance of a hydra effect, particularly when the disease-free population is at a cycle.

Prior exposure to pathogens augments host heterogeneity in susceptibility and has key epidemiological consequences.

Pathogen epidemics are key threats to human and wildlife health. Across systems, host protection from pathogens following initial exposure is often incomplete, resulting in recurrent epidemics through partially-immune hosts. Variation in population-level protection has important consequences for epidemic dynamics, but whether acquired protection influences host heterogeneity in susceptibility and its epidemiological consequences remains unexplored. We experimentally investigated whether prior exposure (none, low-dose, or high-dose) to a bacterial pathogen alters host heterogeneity in susceptibility among songbirds. Hosts with no prior pathogen exposure had little variation in protection, but heterogeneity in susceptibility was significantly augmented by prior pathogen exposure, with the highest variability detected in hosts given high-dose prior exposure. An epidemiological model parameterized with experimental data found that heterogeneity in susceptibility from prior exposure more than halved epidemic sizes compared with a homogeneous population with identical mean protection. However, because infection-induced mortality was also greatly reduced in hosts with prior pathogen exposure, reductions in epidemic size were smaller than expected in hosts with prior exposure. These results highlight the importance of variable protection from prior exposure and/or vaccination in driving host heterogeneity and epidemiological dynamics.

Why Similar Policies Resulted In Different COVID-19 Outcomes: How Responsiveness And Culture Influenced Mortality Rates.

In the first two years of the COVID-19 pandemic, per capita mortality varied by more than a hundredfold across countries, despite most implementing similar nonpharmaceutical interventions. Factors such as policy stringency, gross domestic product, and age distribution explain only a small fraction of mortality variation. To address this puzzle, we built on a previously validated pandemic model in which perceived risk altered societal responses affecting SARS-CoV-2 transmission. Using data from more than 100 countries, we found that a key factor explaining heterogeneous death rates was not the policy responses themselves but rather variation in responsiveness. Responsiveness measures how sensitive communities are to evolving mortality risks and how readily they adopt nonpharmaceutical interventions in response, to curb transmission. We further found that responsiveness correlated with two cultural constructs across countries: uncertainty avoidance and power distance. Our findings show that more responsive adoption of similar policies saves many lives, with important implications for the design and implementation of responses to future outbreaks.

Extraordinary parasite multiplication rates in human malaria infections.

For pathogenic organisms, faster rates of multiplication promote transmission success, the potential to harm hosts, and the evolution of drug resistance. Parasite multiplication rates (PMRs) are often quantified in malaria infections, given the relative ease of sampling. Using modern and historical human infection data, we show that established methods return extraordinarily - and implausibly - large PMRs. We illustrate how inflated PMRs arise from two facets of malaria biology that are far from unique: (i) some developmental ages are easier to sample than others; (ii) the distribution of developmental ages changes over the course of infection. The difficulty of accurately quantifying PMRs demonstrates a need for robust methods and a subsequent re-evaluation of what is known even in the well-studied system of malaria.

Eos Promotes TH2 Differentiation by Interacting with and Propagating the Activity of STAT5.

The Ikaros zinc-finger transcription factor Eos has largely been associated with sustaining the immunosuppressive functions of regulatory T cells. Paradoxically, Eos has more recently been implicated in promoting proinflammatory responses in the dysregulated setting of autoimmunity. However, the precise role of Eos in regulating the differentiation and function of effector CD4+ T cell subsets remains unclear. In this study, we find that Eos is a positive regulator of the differentiation of murine CD4+ TH2 cells, an effector population that has been implicated in both immunity against helminthic parasites and the induction of allergic asthma. Using murine in vitro TH2 polarization and an in vivo house dust mite asthma model, we find that EosKO T cells exhibit reduced expression of key TH2 transcription factors, effector cytokines, and cytokine receptors. Mechanistically, we find that the IL-2/STAT5 axis and its downstream TH2 gene targets are one of the most significantly downregulated pathways in Eos-deficient cells. Consistent with these observations, we find that Eos forms, to our knowledge, a novel complex with and supports the tyrosine phosphorylation of STAT5. Collectively, these data define a regulatory mechanism whereby Eos propagates STAT5 activity to facilitate TH2 cell differentiation.

Multiple models for outbreak decision support in the face of uncertainty.

Policymakers must make management decisions despite incomplete knowledge and conflicting model projections. Little guidance exists for the rapid, representative, and unbiased collection of policy-relevant scientific input from independent modeling teams. Integrating approaches from decision analysis, expert judgment, and model aggregation, we convened multiple modeling teams to evaluate COVID-19 reopening strategies for a mid-sized United States county early in the pandemic. Projections from seventeen distinct models were inconsistent in magnitude but highly consistent in ranking interventions. The 6-mo-ahead aggregate projections were well in line with observed outbreaks in mid-sized US counties. The aggregate results showed that up to half the population could be infected with full workplace reopening, while workplace restrictions reduced median cumulative infections by 82%. Rankings of interventions were consistent across public health objectives, but there was a strong trade-off between public health outcomes and duration of workplace closures, and no win-win intermediate reopening strategies were identified. Between-model variation was high; the aggregate results thus provide valuable risk quantification for decision making. This approach can be applied to the evaluation of management interventions in any setting where models are used to inform decision making. This case study demonstrated the utility of our approach and was one of several multimodel efforts that laid the groundwork for the COVID-19 Scenario Modeling Hub, which has provided multiple rounds of real-time scenario projections for situational awareness and decision making to the Centers for Disease Control and Prevention since December 2020.

Declines in prevalence alter the optimal level of sexual investment for the malaria parasite Plasmodium falciparum.

Successful infectious disease interventions can result in large reductions in parasite prevalence. Such demographic change has fitness implications for individual parasites and may shift the parasite's optimal life history strategy. Here, we explore whether declining infection rates can alter Plasmodium falciparum's investment in sexual versus asexual growth. Using a multiscale mathematical model, we demonstrate how the proportion of polyclonal infections, which decreases as parasite prevalence declines, affects the optimal sexual development strategy: Within-host competition in multiclone infections favors a greater investment in asexual growth whereas single-clone infections benefit from higher conversion to sexual forms. At the same time, drug treatment also imposes selection pressure on sexual development by shortening infection length and reducing within-host competition. We assess these models using 148 P. falciparum parasite genomes sampled in French Guiana over an 18-y period of intensive intervention (1998 to 2015). During this time frame, multiple public health measures, including the introduction of new drugs and expanded rapid diagnostic testing, were implemented, reducing P. falciparum malaria cases by an order of magnitude. Consistent with this prevalence decline, we see an increase in the relatedness among parasites, but no single clonal background grew to dominate the population. Analyzing individual allele frequency trajectories, we identify genes that likely experienced selective sweeps. Supporting our model predictions, genes showing the strongest signatures of selection include transcription factors involved in the development of P. falciparum's sexual gametocyte form. These results highlight how public health interventions impose wide-ranging selection pressures that affect basic parasite life history traits.

Modeling the effects of Aedes aegypti's larval environment on adult body mass at emergence.

Mosquitoes vector harmful pathogens that infect millions of people every year, and developing approaches to effectively control mosquitoes is a topic of great interest. However, the success of many control measures is highly dependent upon ecological, physiological, and life history traits of mosquito species. The behavior of mosquitoes and their potential to vector pathogens can also be impacted by these traits. One trait of interest is mosquito body mass, which depends upon many factors associated with the environment in which juvenile mosquitoes develop. Our experiments examined the impact of larval density on the body mass of Aedes aegypti mosquitoes, which are important vectors of dengue, Zika, yellow fever, and other pathogens. To investigate the interactions between the larval environment and mosquito body mass, we built a discrete time mathematical model that incorporates body mass, larval density, and food availability and fit the model to our experimental data. We considered three categories of model complexity informed by data, and selected the best model within each category using Akaike's Information Criterion. We found that the larval environment is an important determinant of the body mass of mosquitoes upon emergence. Furthermore, we found that larval density has greater impact on body mass of adults at emergence than on development time, and that inclusion of density dependence in the survival of female aquatic stages in models is important. We discuss the implications of our results for the control of Aedes mosquitoes and on their potential to spread disease.

Parasitism of Aedes albopictus by Ascogregarina taiwanensis lowers its competitive ability against Aedes triseriatus.

BACKGROUND: Mosquitoes are vectors for diseases such as dengue, malaria and La Crosse virus that significantly impact the human population. When multiple mosquito species are present, the competition between species may alter population dynamics as well as disease spread. Two mosquito species, Aedes albopictus and Aedes triseriatus, both inhabit areas where La Crosse virus is found. Infection of Aedes albopictus by the parasite Ascogregarina taiwanensis and Aedes triseriatus by the parasite Ascogregarina barretti can decrease a mosquito's fitness, respectively. In particular, the decrease in fitness of Aedes albopictus occurs through the impact of Ascogregarina taiwanensis on female fecundity, larval development rate, and larval mortality and may impact its initial competitive advantage over Aedes triseriatus during invasion. METHODS: We examine the effects of parasitism of gregarine parasites on Aedes albopictus and triseriatus population dynamics and competition with a focus on when Aedes albopictus is new to an area. We build a compartmental model including competition between Aedes albopictus and triseriatus while under parasitism of the gregarine parasites. Using parameters based on the literature, we simulate the dynamics and analyze the equilibrium population proportion of the two species. We consider the presence of both parasites and potential dilution effects. RESULTS: We show that increased levels of parasitism in Aedes albopictus will decrease the initial competitive advantage of the species over Aedes triseriatus and increase the survivorship of Aedes triseriatus. We find Aedes albopictus is better able to invade when there is more extreme parasitism of Aedes triseriatus. Furthermore, although the transient dynamics differ, dilution of the parasite density through uptake by both species does not alter the equilibrium population sizes of either species. CONCLUSIONS: Mosquito population dynamics are affected by many factors, such as abiotic factors (e.g. temperature and humidity) and competition between mosquito species. This is especially true when multiple mosquito species are vying to live in the same area. Knowledge of how population dynamics are affected by gregarine parasites among competing species can inform future mosquito control efforts and help prevent the spread of vector-borne disease.

Interactions Between Opioids and Dextroamphetamine on Locomotor Activity: Influence of an Opioid's Relative Efficacy at the Mu Receptor.

Opioids and stimulants are often used in combination for both recreational and non-recreational purposes. High-efficacy mu opioid agonists generally increase the behavioral effects of stimulants, whereas opioid receptor antagonists generally attenuate the behavioral effects of stimulants; however, less is known regarding the interactions between stimulants and opioids possessing low to intermediate efficacy at the mu receptor. The purpose of this study was to examine the role of an opioid's relative efficacy at the mu receptor in altering the behavioral effects of dextro(d-)amphetamine. To this end, opioids possessing a range of relative efficacy at the mu receptor were examined alone and in combination with cumulative doses of d-amphetamine on a test of open-field, locomotor activity in male rats. Levorphanol, buprenorphine, butorphanol, nalbuphine, (-)-pentazocine, (-)-metazocine, (-)-cyclazocine, (-)-NANM, and nalorphine increased the locomotor effects of d-amphetamine in either an additive or greater-than-additive manner according to an effect-additive model. Only the selective, high-efficacy kappa agonist, spiradoline, and the non-selective opioid receptor antagonist, naloxone, failed to increase the effects of d-amphetamine under the conditions examined. These data indicate that opioids possessing a large range of relative efficacy at the mu receptor, including those possessing very low relative efficacy, significantly increase the locomotor effects of d-amphetamine.

Multiple blood feeding in mosquitoes shortens the Plasmodium falciparum incubation period and increases malaria transmission potential.

Many mosquito species, including the major malaria vector Anopheles gambiae, naturally undergo multiple reproductive cycles of blood feeding, egg development and egg laying in their lifespan. Such complex mosquito behavior is regularly overlooked when mosquitoes are experimentally infected with malaria parasites, limiting our ability to accurately describe potential effects on transmission. Here, we examine how Plasmodium falciparum development and transmission potential is impacted when infected mosquitoes feed an additional time. We measured P. falciparum oocyst size and performed sporozoite time course analyses to determine the parasite's extrinsic incubation period (EIP), i.e. the time required by parasites to reach infectious sporozoite stages, in An. gambiae females blood fed either once or twice. An additional blood feed at 3 days post infection drastically accelerates oocyst growth rates, causing earlier sporozoite accumulation in the salivary glands, thereby shortening the EIP (reduction of 2.3 ± 0.4 days). Moreover, parasite growth is further accelerated in transgenic mosquitoes with reduced reproductive capacity, which mimic genetic modifications currently proposed in population suppression gene drives. We incorporate our shortened EIP values into a measure of transmission potential, the basic reproduction number R0, and find the average R0 is higher (range: 10.1%-12.1% increase) across sub-Saharan Africa than when using traditional EIP measurements. These data suggest that malaria elimination may be substantially more challenging and that younger mosquitoes or those with reduced reproductive ability may provide a larger contribution to infection than currently believed. Our findings have profound implications for current and future mosquito control interventions.

The impact of within-vector parasite development on the extrinsic incubation period.

Mosquito-borne diseases, in particular malaria, have a significant burden worldwide leading to nearly half a million deaths each year. The malaria parasite requires a vertebrate host, such as a human, and a vector host, the Anopheles mosquito, to complete its full life cycle. Here, we focus on the parasite dynamics within the vector to examine the first appearance of sporozoites in the salivary glands, which indicates a first time of infectiousness of mosquitoes. The timing of this period of pathogen development in the mosquito until transmissibility, known as the extrinsic incubation period, remains poorly understood. We develop compartmental models of within-mosquito parasite dynamics fitted with experimental data on oocyst and sporozoite counts. We find that only a fraction of oocysts burst to release sporozoites and bursting must be delayed either via a time-dependent function or a gamma-distributed set of compartments. We use Bayesian inference to estimate distributions of parameters and determine that bursting rate is a key epidemiological parameter. A better understanding of the factors impacting the extrinsic incubation period will aid in the development of interventions to slow or stop the spread of malaria.

COVID-19 reopening strategies at the county level in the face of uncertainty: Multiple Models for Outbreak Decision Support.

Policymakers make decisions about COVID-19 management in the face of considerable uncertainty. We convened multiple modeling teams to evaluate reopening strategies for a mid-sized county in the United States, in a novel process designed to fully express scientific uncertainty while reducing linguistic uncertainty and cognitive biases. For the scenarios considered, the consensus from 17 distinct models was that a second outbreak will occur within 6 months of reopening, unless schools and non-essential workplaces remain closed. Up to half the population could be infected with full workplace reopening; non-essential business closures reduced median cumulative infections by 82%. Intermediate reopening interventions identified no win-win situations; there was a trade-off between public health outcomes and duration of workplace closures. Aggregate results captured twice the uncertainty of individual models, providing a more complete expression of risk for decision-making purposes.

Comparative Impact of Individual Quarantine vs. Active Monitoring of Contacts for the Mitigation of COVID-19: a modelling study.

Background: Voluntary individual quarantine and voluntary active monitoring of contacts are core disease control strategies for emerging infectious diseases, such as COVID-19. Given the impact of quarantine on resources and individual liberty, it is vital to assess under what conditions individual quarantine can more effectively control COVID-19 than active monitoring. As an epidemic grows, it is also important to consider when these interventions are no longer feasible, and broader mitigation measures must be implemented. Methods: To estimate the comparative efficacy of these case-based interventions to control COVID-19, we fit a stochastic branching model to reported parameters for the dynamics of the disease. Specifically, we fit to the incubation period distribution and each of two sets of the serial interval distribution: a shorter one with a mean serial interval of 4.8 days and a longer one with a mean of 7.5 days. To assess variable resource settings, we consider two feasibility settings: a high feasibility setting with 90% of contacts traced, a half-day average delay in tracing and symptom recognition, and 90% effective isolation; and low feasibility setting with 50% of contacts traced, a two-day average delay, and 50% effective isolation. Findings: Our results suggest that individual quarantine in high feasibility settings where at least three-quarters of infected contacts are individually quarantined contains an outbreak of COVID-19 with a short serial interval (4.8 days) 84% of the time. However, in settings where this performance is unrealistically high and the outbreak continues to grow, so too will the burden of the number of contacts traced for active monitoring or quarantine. When resources are prioritized for scalable interventions such as social distancing, we show active monitoring or individual quarantine of high-risk contacts can contribute synergistically to mitigation efforts. Interpretation: Our model highlights the urgent need for more data on the serial interval and the extent of presymptomatic transmission in order to make data-driven policy decisions regarding the cost-benefit comparisons of individual quarantine vs. active monitoring of contacts. To the extent these interventions can be implemented they can help mitigate the spread of COVID-19.

Individual quarantine versus active monitoring of contacts for the mitigation of COVID-19: a modelling study.

BACKGROUND: Voluntary individual quarantine and voluntary active monitoring of contacts are core disease control strategies for emerging infectious diseases such as COVID-19. Given the impact of quarantine on resources and individual liberty, it is vital to assess under what conditions individual quarantine can more effectively control COVID-19 than active monitoring. As an epidemic grows, it is also important to consider when these interventions are no longer feasible and broader mitigation measures must be implemented. METHODS: To estimate the comparative efficacy of individual quarantine and active monitoring of contacts to control severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), we fit a stochastic branching model to reported parameters for the dynamics of the disease. Specifically, we fit a model to the incubation period distribution (mean 5·2 days) and to two estimates of the serial interval distribution: a shorter one with a mean serial interval of 4·8 days and a longer one with a mean of 7·5 days. To assess variable resource settings, we considered two feasibility settings: a high-feasibility setting with 90% of contacts traced, a half-day average delay in tracing and symptom recognition, and 90% effective isolation; and a low-feasibility setting with 50% of contacts traced, a 2-day average delay, and 50% effective isolation. FINDINGS: Model fitting by sequential Monte Carlo resulted in a mean time of infectiousness onset before symptom onset of 0·77 days (95% CI -1·98 to 0·29) for the shorter serial interval, and for the longer serial interval it resulted in a mean time of infectiousness onset after symptom onset of 0·51 days (95% CI -0·77 to 1·50). Individual quarantine in high-feasibility settings, where at least 75% of infected contacts are individually quarantined, contains an outbreak of SARS-CoV-2 with a short serial interval (4·8 days) 84% of the time. However, in settings where the outbreak continues to grow (eg, low-feasibility settings), so too will the burden of the number of contacts traced for active monitoring or quarantine, particularly uninfected contacts (who never develop symptoms). When resources are prioritised for scalable interventions such as physical distancing, we show active monitoring or individual quarantine of high-risk contacts can contribute synergistically to mitigation efforts. Even under the shorter serial interval, if physical distancing reduces the reproductive number to 1·25, active monitoring of 50% of contacts can result in overall outbreak control (ie, effective reproductive number <1). INTERPRETATION: Our model highlights the urgent need for more data on the serial interval and the extent of presymptomatic transmission to make data-driven policy decisions regarding the cost-benefit comparisons of individual quarantine versus active monitoring of contacts. To the extent that these interventions can be implemented, they can help mitigate the spread of SARS-CoV-2. FUNDING: National Institute of General Medical Sciences, National Institutes of Health.

Modeling the waning and boosting of immunity from infection or vaccination.

Immunity following natural infection or immunization may wane, increasing susceptibility to infection with time since infection or vaccination. Symptoms, and concomitantly infectiousness, depend on residual immunity. We quantify these phenomena in a model population composed of individuals whose susceptibility, infectiousness, and symptoms all vary with immune status. We also model age, which affects contact, vaccination and possibly waning rates. The resurgences of pertussis that have been observed wherever effective vaccination programs have reduced typical disease among young children follow from these processes. As one example, we compare simulations with the experience of Sweden following resumption of pertussis vaccination after the hiatus from 1979 to 1996, reproducing the observations leading health authorities to introduce booster doses among school-aged children and adolescents in 2007 and 2014, respectively. Because pertussis comprises a spectrum of symptoms, only the most severe of which are medically attended, accurate models are needed to design optimal vaccination programs where surveillance is less effective.

Mathematical model of broadly reactive plasma cell production.

Strain-specific plasma cells are capable of producing neutralizing antibodies that are essential for clearance of challenging pathogens. These neutralizing antibodies also function as a main defense against disease establishment in a host. However, when a rapidly mutating pathogen infects a host, successful control of the invasion requires shifting the production of plasma cells from strain-specific to broadly reactive. In this study, we develop a mathematical model of germinal center dynamics and use it to predict the events that lead to improved breadth of the plasma cell response. We examine scenarios that lead to germinal centers that are composed of B-cells that come from a single strain-specific clone, a single broadly reactive clone or both clones. We find that the initial B-cell clonal composition, T-follicular helper cell signaling, increased rounds of productive somatic hypermutation, and B-cell selection strength are among the mechanisms differentiating between strain-specific and broadly reactive plasma cell production during infections. Understanding the contribution of these factors to emergence of breadth may assist in boosting broadly reactive plasma cells production.

Evolutionary consequences of feedbacks between within-host competition and disease control.

Lay Summary: Competition often occurs among diverse parasites within a single host, but control efforts could change its strength. We examined how the interplay between competition and control could shape the evolution of parasite traits like drug resistance and disease severity.

Linked within-host and between-host models and data for infectious diseases: a systematic review.

The observed dynamics of infectious diseases are driven by processes across multiple scales. Here we focus on two: within-host, that is, how an infection progresses inside a single individual (for instance viral and immune dynamics), and between-host, that is, how the infection is transmitted between multiple individuals of a host population. The dynamics of each of these may be influenced by the other, particularly across evolutionary time. Thus understanding each of these scales, and the links between them, is necessary for a holistic understanding of the spread of infectious diseases. One approach to combining these scales is through mathematical modeling. We conducted a systematic review of the published literature on multi-scale mathematical models of disease transmission (as defined by combining within-host and between-host scales) to determine the extent to which mathematical models are being used to understand across-scale transmission, and the extent to which these models are being confronted with data. Following the PRISMA guidelines for systematic reviews, we identified 24 of 197 qualifying papers across 30 years that include both linked models at the within and between host scales and that used data to parameterize/calibrate models. We find that the approach that incorporates both modeling with data is under-utilized, if increasing. This highlights the need for better communication and collaboration between modelers and empiricists to build well-calibrated models that both improve understanding and may be used for prediction.

Exposing Anopheles mosquitoes to antimalarials blocks Plasmodium parasite transmission.

Bites of Anopheles mosquitoes transmit Plasmodium falciparum parasites that cause malaria, which kills hundreds of thousands of people every year. Since the turn of this century, efforts to prevent the transmission of these parasites via the mass distribution of insecticide-treated bed nets have been extremely successful, and have led to an unprecedented reduction in deaths from malaria1. However, resistance to insecticides has become widespread in Anopheles populations2-4, which has led to the threat of a global resurgence of malaria and makes the generation of effective tools for controlling this disease an urgent public health priority. Here we show that the development of P. falciparum can be rapidly and completely blocked when female Anopheles gambiae mosquitoes take up low concentrations of specific antimalarials from treated surfaces-conditions that simulate contact with a bed net. Mosquito exposure to atovaquone before, or shortly after, P. falciparum infection causes full parasite arrest in the midgut, and prevents transmission of infection. Similar transmission-blocking effects are achieved using other cytochrome b inhibitors, which demonstrates that parasite mitochondrial function is a suitable target for killing parasites. Incorporating these effects into a model of malaria transmission dynamics predicts that impregnating mosquito nets with Plasmodium inhibitors would substantially mitigate the global health effects of insecticide resistance. This study identifies a powerful strategy for blocking Plasmodium transmission by female Anopheles mosquitoes, which has promising implications for efforts to eradicate malaria.