Partitioning regulatory mechanisms of within-host malaria dynamics using the effective propagation number

Synchrony in malaria infections: how intensifying within-host competition can be adaptive

Malaria parasites exhibit great diversity in the coordination of their asexual life cycle within the host, ranging from asynchronous growth to tightly synchronized cycles of invasion and emergence from red blood cells. Synchronized reproduction should come at a high cost--intensifying competition among offspring--so why would some Plasmodium species engage in such behavior and others not? We use a delayed differential equation model to show that synchronized infections can be favored when (1) there is limited interference among parasites competing for red blood cells, (2) transmission success is an accelerating function of sexual parasite abundance, (3) the target of saturating immunity is short-lived, and (4) coinfections with asynchronous parasites are rare. As a consequence, synchrony may be beneficial or costly, in line with the diverse patterns of synchronization observed in natural and lab infections. By allowing us to characterize diverse temporal dynamics, the model framework provides a basis for making predictions about disease severity and for projecting evolutionary responses to interventions.

Disrupting rhythms in Plasmodium chabaudi: costs accrue quickly and independently of how infections are initiated

BACKGROUND

In the blood, the synchronous malaria parasite, Plasmodium chabaudi, exhibits a cell-cycle rhythm of approximately 24 hours in which transitions between developmental stages occur at particular times of day in the rodent host. Previous experiments reveal that when the timing of the parasite's cell-cycle rhythm is perturbed relative to the circadian rhythm of the host, parasites suffer a (~50%) reduction in asexual stages and gametocytes. Why it matters for parasites to have developmental schedules in synchronization with the host's rhythm is unknown. The experiment presented here investigates this issue by: (a) validating that the performance of P. chabaudi is negatively affected by mismatch to the host circadian rhythm; (b) testing whether the effect of mismatch depends on the route of infection or the developmental stage of inoculated parasites; and, (c) examining whether the costs of mismatch are due to challenges encountered upon initial infection and/or due to ongoing circadian host processes operating during infection.

METHODS

The experiment simultaneously perturbed the time of day infections were initiated, the stage of parasite inoculated, and the route of infection. The performance of parasites during the growth phase of infections was compared across the cross-factored treatment groups (i e, all combinations of treatments were represented).

RESULTS

The data show that mismatch to host rhythms is costly for parasites, reveal that this phenomenon does not depend on the developmental stage of parasites nor the route of infection, and suggest that processes operating at the initial stages of infection are responsible for the costs of mismatch. Furthermore, mismatched parasites are less virulent, in that they cause less anaemia to their hosts.

CONCLUSION

It is beneficial for parasites to be in synchronization with their host's rhythm, regardless of the route of infection or the parasite stage inoculated. Given that arrested cell-cycle development (quiescence) is implicated in tolerance to drugs, understanding how parasite schedules are established and maintained in the blood is important.

Expansion of host cellular niche can drive adaptation of a zoonotic malaria parasite to humans

The macaque malaria parasite Plasmodium knowlesi has recently emerged as an important zoonosis in Southeast Asia. Infections are typically mild but can cause severe disease, achieving parasite densities similar to fatal Plasmodium falciparum infections. Here we show that a primate-adapted P. knowlesi parasite proliferates poorly in human blood due to a strong preference for young red blood cells (RBCs). We establish a continuous in vitro culture system by using human blood enriched for young cells. Mathematical modelling predicts that parasite adaptation for invasion of older RBCs is a likely mechanism leading to high parasite densities in clinical infections. Consistent with this model, we find that P. knowlesi can adapt to invade a wider age range of RBCs, resulting in proliferation in normal human blood. Such cellular niche expansion may increase pathogenesis in humans and will be a key feature to monitor as P. knowlesi emerges in human populations.

Ahead of the curve: next generation estimators of drug resistance in malaria infections

Drug resistance is a major obstacle to controlling infectious diseases. A key challenge is detecting the early signs of drug resistance when little is known about its genetic basis. Focusing on malaria parasites, we propose a way to do this. Newly developing or low level resistance at low frequency in patients can be detected through a phenotypic signature: individual parasite variants clearing more slowly following drug treatment. Harnessing the abundance and resolution of deep sequencing data, our 'selection differential' approach addresses some limitations of extant methods of resistance detection, should allow for the earliest detection of resistance in malaria or other multi-clone infections, and has the power to uncover the true scale of the drug resistance problem.

Disease invasion: impacts on biodiversity and human health

An introduction to the theme issue that includes papers that identify how, where and why infectious diseases in wildlife emerge, while also addressing their possible conservation impacts.

The Cinderella syndrome: why do malaria-infected cells burst at midnight?

An interesting quirk of many malaria infections is that all parasites within a host – millions of them – progress through their cell cycle synchronously. This surprising coordination has long been recognized, yet there is little understanding of what controls it or why it has evolved. Interestingly, the conventional explanation for coordinated development in other parasite species does not seem to apply here. We argue that for malaria parasites, a critical question has yet to be answered: is the coordination due to parasites bursting at the same time or at a particular time? We explicitly delineate these fundamentally different scenarios, possible underlying mechanistic explanations and evolutionary drivers, and discuss the existing corroborating data and key evidence needed to solve this evolutionary mystery.

Mobile phones and malaria: modeling human and parasite travel

Human mobility plays an important role in the dissemination of malaria parasites between regions of variable transmission intensity. Asymptomatic individuals can unknowingly carry parasites to regions where mosquito vectors are available, for example, undermining control programs and contributing to transmission when they travel. Understanding how parasites are imported between regions in this way is therefore an important goal for elimination planning and the control of transmission, and would enable control programs to target the principal sources of malaria. Measuring human mobility has traditionally been difficult to do on a population scale, but the widespread adoption of mobile phones in low-income settings presents a unique opportunity to directly measure human movements that are relevant to the spread of malaria. Here, we discuss the opportunities for measuring human mobility using data from mobile phones, as well as some of the issues associated with combining mobility estimates with malaria infection risk maps to meaningfully estimate routes of parasite importation.

Malaria parasite diversity and transmission intensity affect development of parasitological immunity in a mathematical model

BACKGROUND

The development of parasitological immunity against malaria affects the ability to detect infection, the efficiency of the local human parasite reservoir at infecting mosquitoes, and the response to reintroduction of parasites to previously cleared areas. Observations of similar age-trends in detected prevalence and mean parasitaemia across more than an order-of-magnitude of variation in baseline transmission complicate simple exposure-driven explanations.

METHODS

Mathematical models often employ age-dependent immune factors to match the observed trends, while the present model uses a new detailed mechanistic model of parasite transmission dynamics to explain age-trends through the mechanism of parasite diversity. Illustrative simulations are performed for multiple field sites in Tanzania and Nigeria, and observed age-trends and seasonality in parasite prevalence are recreated in silico, proffering possible mechanistic explanations of the observational data.

RESULTS

Observed temporal dynamics in measured parasitaemia are recreated for each location and age-prevalence outputs are studied. Increasing population-level diversity in malaria surface antigens delays development of broad parasitological immunity. A local parasite population with high diversity can recreate the observed trends in age-prevalence across more than an order of magnitude of variation in transmission intensities.

CONCLUSIONS

Mechanistic models of human immunity and parasite antigen diversity can recreate the observed temporal patterns for the development of parasitological immunity across a wide range of transmission intensities. This has implications for the distribution of disease burden across the population, the human transmission reservoir, design of elimination campaigns, and development and roll-out of potential vaccines.

Quantification of the effects of antibodies on the extra- and intracellular dynamics of Salmonella enterica

Antibodies are known to be essential in controlling Salmonella infection, but their exact role remains elusive. We recently developed an in vitro model to investigate the relative efficiency of four different human immunoglobulin G (IgG) subclasses in modulating the interaction of the bacteria with human phagocytes. Our results indicated that different IgG subclasses affect the efficacy of Salmonella uptake by human phagocytes. In this study, we aim to quantify the effects of IgG on intracellular dynamics of infection by combining distributions of bacterial numbers per phagocyte observed by fluorescence microscopy with a mathematical model that simulates the in vitro dynamics. We then use maximum likelihood to estimate the model parameters and compare them across IgG subclasses. The analysis reveals heterogeneity in the division rates of the bacteria, strongly suggesting that a subpopulation of intracellular Salmonella, while visible under the microscope, is not dividing. Clear differences in the observed distributions among the four IgG subclasses are best explained by variations in phagocytosis and intracellular dynamics. We propose and compare potential factors affecting the replication and death of bacteria within phagocytes, and we discuss these results in the light of recent findings on dormancy of Salmonella.

Ecology of zoonoses: natural and unnatural histories

More than 60% of human infectious diseases are caused by pathogens shared with wild or domestic animals. Zoonotic disease organisms include those that are endemic in human populations or enzootic in animal populations with frequent cross-species transmission to people. Some of these diseases have only emerged recently. Together, these organisms are responsible for a substantial burden of disease, with endemic and enzootic zoonoses causing about a billion cases of illness in people and millions of deaths every year. Emerging zoonoses are a growing threat to global health and have caused hundreds of billions of US dollars of economic damage in the past 20 years. We aimed to review how zoonotic diseases result from natural pathogen ecology, and how other circumstances, such as animal production, extraction of natural resources, and antimicrobial application change the dynamics of disease exposure to human beings. In view of present anthropogenic trends, a more effective approach to zoonotic disease prevention and control will require a broad view of medicine that emphasises evidence-based decision making and integrates ecological and evolutionary principles of animal, human, and environmental factors. This broad view is essential for the successful development of policies and practices that reduce probability of future zoonotic emergence, targeted surveillance and strategic prevention, and engagement of partners outside the medical community to help improve health outcomes and reduce disease threats.

Mechanism-based model of parasite growth and dihydroartemisinin pharmacodynamics in murine malaria

Murine models are used to study erythrocytic stages of malaria infection, because parasite morphology and development are comparable to those in human malaria infections. Mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) models for antimalarials are scarce, despite their potential to optimize antimalarial combination therapy. The aim of this study was to develop a mechanism-based growth model (MBGM) for Plasmodium berghei and then characterize the parasiticidal effect of dihydroartemisinin (DHA) in murine malaria (MBGM-PK-PD). Stage-specific (ring, early trophozoite, late trophozoite, and schizont) parasite density data from Swiss mice inoculated with Plasmodium berghei were used for model development in S-ADAPT. A single dose of intraperitoneal DHA (10 to 100 mg/kg) or vehicle was administered 56 h postinoculation. The MBGM explicitly reflected all four erythrocytic stages of the 24-hour P. berghei life cycle. Merozoite invasion of erythrocytes was described by a first-order process that declined with increasing parasitemia. An efflux pathway with subsequent return was additionally required to describe the schizont data, thus representing parasite sequestration or trapping in the microvasculature, with a return to circulation. A 1-compartment model with zero-order absorption described the PK of DHA, with an estimated clearance and distribution volume of 1.95 liters h(-1) and 0.851 liter, respectively. Parasite killing was described by a turnover model, with DHA inhibiting the production of physiological intermediates (IC(50), 1.46 ng/ml). Overall, the MBGM-PK-PD described the rise in parasitemia, the nadir following DHA dosing, and subsequent parasite resurgence. This novel model is a promising tool for studying malaria infections, identifying the stage specificity of antimalarials, and providing insight into antimalarial treatment strategies.

Gene-gene interaction and functional impact of polymorphisms on innate immune genes in controlling Plasmodium falciparum blood infection level

Genetic variations in toll-like receptors and cytokine genes of the innate immune pathways have been implicated in controlling parasite growth and the pathogenesis of Plasmodium falciparum mediated malaria. We previously published genetic association of TLR4 non-synonymous and TNF-α promoter polymorphisms with P.falciparum blood infection level and here we extend the study considerably by (i) investigating genetic dependence of parasite-load on interleukin-12B polymorphisms, (ii) reconstructing gene-gene interactions among candidate TLRs and cytokine loci, (iii) exploring genetic and functional impact of epistatic models and (iv) providing mechanistic insights into functionality of disease-associated regulatory polymorphisms. Our data revealed that carriage of AA (P = 0.0001) and AC (P = 0.01) genotypes of IL12B 3′UTR polymorphism was associated with a significant increase of mean log-parasitemia relative to rare homozygous genotype CC. Presence of IL12B+1188 polymorphism in five of six multifactor models reinforced its strong genetic impact on malaria phenotype. Elevation of genetic risk in two-component models compared to the corresponding single locus and reduction of IL12B (2.2 fold) and lymphotoxin-α (1.7 fold) expressions in patients'peripheral-blood-mononuclear-cells under TLR4Thr399Ile risk genotype background substantiated the role of Multifactor Dimensionality Reduction derived models. Marked reduction of promoter activity of TNF-α risk haplotype (C-C-G-G) compared to wild-type haplotype (T-C-G-G) with (84%) and without (78%) LPS stimulation and the loss of binding of transcription factors detected in-silico supported a causal role of TNF-1031. Significantly lower expression of IL12B+1188 AA (5 fold) and AC (9 fold) genotypes compared to CC and under-representation (P = 0.0048) of allele A in transcripts of patients' PBMCs suggested an Allele-Expression-Imbalance. Allele (A+1188C) dependent differential stability (2 fold) of IL12B-transcripts upon actinomycin-D treatment and observed structural modulation (P = 0.013) of RNA-ensemble were the plausible explanations for AEI. In conclusion, our data provides functional support to the hypothesis that de-regulated receptor-cytokine axis of innate immune pathway influences blood infection level in P. falciparum malaria.

P. falciparum infection durations and infectiousness are shaped by antigenic variation and innate and adaptive host immunity in a mathematical model

Many questions remain about P. falciparum within-host dynamics, immunity, and transmission–issues that may affect public health campaign planning. These gaps in knowledge concern the distribution of durations of malaria infections, determination of peak parasitemia during acute infection, the relationships among gametocytes and immune responses and infectiousness to mosquitoes, and the effect of antigenic structure on reinfection outcomes. The present model of intra-host dynamics of P. falciparum implements detailed representations of parasite and immune dynamics, with structures based on minimal extrapolations from first-principles biology in its foundations. The model is designed to quickly and readily accommodate gains in mechanistic understanding and to evaluate effects of alternative biological hypothesis through in silico experiments. Simulations follow the parasite from the liver-stage through the detailed asexual cycle to clearance while tracking gametocyte populations. The modeled immune system includes innate inflammatory and specific antibody responses to a repertoire of antigens. The mechanistic focus provides clear explanations for the structure of the distribution of infection durations through the interaction of antigenic variation and innate and adaptive immunity. Infectiousness to mosquitoes appears to be determined not only by the density of gametocytes but also by the level of inflammatory cytokines, which harmonizes an extensive series of study results. Finally, pre-existing immunity can either decrease or increase the duration of infections upon reinfection, depending on the degree of overlap in antigenic repertoires and the strength of the pre-existing immunity.

The evolutionary consequences of blood-stage vaccination on the rodent malaria Plasmodium chabaudi

A candidate malaria vaccine promoted the evolution of more virulent malaria parasites in mice.

Malaria vaccine developers are concerned that antigenic escape will erode vaccine efficacy. Evolutionary theorists have raised the possibility that some types of vaccine could also create conditions favoring the evolution of more virulent pathogens. Such evolution would put unvaccinated people at greater risk of severe disease. Here we test the impact of vaccination with a single highly purified antigen on the malaria parasite Plasmodium chabaudi evolving in laboratory mice. The antigen we used, AMA-1, is a component of several candidate malaria vaccines currently in various stages of trials in humans. We first found that a more virulent clone was less readily controlled by AMA-1-induced immunity than its less virulent progenitor. Replicated parasites were then serially passaged through control or AMA-1 vaccinated mice and evaluated after 10 and 21 rounds of selection. We found no evidence of evolution at the ama-1 locus. Instead, virulence evolved; AMA-1-selected parasites induced greater anemia in naïve mice than both control and ancestral parasites. Our data suggest that recombinant blood stage malaria vaccines can drive the evolution of more virulent malaria parasites.

Vaccination can drive the evolution of pathogens. Most obviously, molecules targeted by vaccine-induced immunity can change. Such evolution makes vaccines less effective. A different possibility is that more virulent pathogens are favored in vaccinated hosts. In that case, vaccination would create pathogens that cause more harm to unvaccinated individuals. To test this idea, we studied a rodent malaria parasite in laboratory mice immunized with a component of malaria vaccines currently in human trials. We found that a more virulent parasite clone was less well controlled by vaccine-induced immunity than was its less virulent ancestor. We then passaged parasites through sham- or vaccinated mice to study how the parasites might evolve after multiple rounds of infection of mouse hosts. The parasite molecule targeted by the vaccine did not change during this process. Instead, the parasites became more virulent if they evolved in vaccinated hosts. Our data suggest that some vaccines can drive the evolution of more virulent parasites.

Revealing mechanisms underlying variation in malaria virulence: effective propagation and host control of uninfected red blood cell supply

Malaria parasite clones with the highest transmission rates to mosquitoes also tend to induce the most severe fitness consequences (or virulence) in mammals. This is in accord with expectations from the virulence–transmission trade-off hypothesis. However, the mechanisms underlying how different clones cause virulence are not well understood. Here, using data from eight murine malaria clones, we apply recently developed statistical methods to infer differences in clone characteristics, including induction of differing host-mediated changes in red blood cell (RBC) supply. Our results indicate that the within-host mechanisms underlying similar levels of virulence are variable and that killing of uninfected RBCs by immune effectors and/or retention of RBCs in the spleen may ultimately reduce virulence. Furthermore, the correlation between clone virulence and the degree of host-induced mortality of uninfected RBCs indicates that hosts increasingly restrict their RBC supply with increasing intrinsic virulence of the clone with which they are infected. Our results demonstrate a role for self-harm in self-defence for hosts and highlight the diversity and modes of virulence of malaria.

Pierre TG, Davis TME. Mathematical modeling of malaria infection with innate and adaptive immunity in individuals and agent-based communities

Background

Agent-based modeling of Plasmodium falciparum infection offers an attractive alternative to the conventional Ross-Macdonald methodology, as it allows simulation of heterogeneous communities subjected to realistic transmission (inoculation patterns).

Methodology/Principal Findings

We developed a new, agent based model that accounts for the essential in-host processes: parasite replication and its regulation by innate and adaptive immunity. The model also incorporates a simplified version of antigenic variation by Plasmodium falciparum. We calibrated the model using data from malaria-therapy (MT) studies, and developed a novel calibration procedure that accounts for a deterministic and a pseudo-random component in the observed parasite density patterns. Using the parasite density patterns of 122 MT patients, we generated a large number of calibrated parameters. The resulting data set served as a basis for constructing and simulating heterogeneous agent-based (AB) communities of MT-like hosts. We conducted several numerical experiments subjecting AB communities to realistic inoculation patterns reported from previous field studies, and compared the model output to the observed malaria prevalence in the field. There was overall consistency, supporting the potential of this agent-based methodology to represent transmission in realistic communities.

Conclusions/Significance

Our approach represents a novel, convenient and versatile method to model Plasmodium falciparum infection.

Plasticity in parasite phenotypes: evolutionary and ecological implications for disease

Preventing disease is a major goal of applied bioscience and explaining variation in the harm caused by parasites, and their infectiousness, are major goals of evolutionary biology. The emerging field of evolutionary medicine integrates these two ambitions to inform the development of control strategies that retard or withstand unfavorable parasite evolution. However, as parasites live in hostile and changeable environments – the bodies of other organisms – the success of integrating evolutionary biology with medicine requires a better understanding of how natural selection has solved the problems parasites face. There is increasing appreciation that natural selection shapes parasite strategies to survive in the host and transmit between hosts through facultative (plastic) shifts in parasite traits expressed during infections and in different hosts. This article describes how integrating parasite plasticity into biomedical thinking is central to explaining disease outcomes and transmission patterns, as well as predicting the success of control measures.

The fitness of drug-resistant malaria parasites in a rodent model: multiplicity of infection

Malaria infections normally consist of more than one clonally replicating lineage. Within-host interactions between sensitive and resistant parasites can have profound effects on the evolution of drug resistance. Here, using the Plasmodium chabaudi mouse malaria model, we ask whether the costs and benefits of resistance are affected by the number of co-infecting strains competing with a resistant clone. We found strong competitive suppression of resistant parasites in untreated infections and marked competitive release following treatment. The magnitude of competitive suppression depended on competitor identity. However, there was no overall effect of the diversity of susceptible parasites on the extent of competitive suppression or release. If these findings generalize, then transmission intensity will impact on resistance evolution because of its effect on the frequency of mixed infections, not because of its effect on the distribution of clones per host. This would greatly simplify the computational problems of adequately capturing within-host ecology in models of drug resistance evolution in malaria.