Investigating antigenic variation and other parasitehost interactions in Plasmodium falciparum infections in nave hosts

CD36-A Host Receptor Necessary for Malaria Parasites to Establish and Maintain Infection

Plasmodium falciparum-infected erythrocytes (PfIEs) present P. falciparum erythrocyte membrane protein 1 proteins (PfEMP1s) on the cell surface, via which they cytoadhere to various endothelial cell receptors (ECRs) on the walls of human blood vessels. This prevents the parasite from passing through the spleen, which would lead to its elimination. Each P. falciparum isolate has about 60 different PfEMP1s acting as ligands, and at least 24 ECRs have been identified as interaction partners. Interestingly, in every parasite genome sequenced to date, at least 75% of the encoded PfEMP1s have a binding domain for the scavenger receptor CD36 widely distributed on host endothelial cells and many other cell types. Here, we discuss why the interaction between PfIEs and CD36 is optimal to maintain a finely regulated equilibrium that allows the parasite to multiply and spread while causing minimal harm to the host in most infections.

Patient variability in the blood-stage dynamics of Plasmodium falciparum captured by clustering historical data

BACKGROUND

Mathematical models provide an understanding of the dynamics of a Plasmodium falciparum blood-stage infection (within-host models), and can predict the impact of control strategies that affect the blood-stage of malaria. However, the dynamics of P. falciparum blood-stage infections are highly variable between individuals. Within-host models use different techniques to capture this inter-individual variation. This struggle may be unnecessary because patients can be clustered according to similar key within-host dynamics. This study aimed to identify clusters of patients with similar parasitaemia profiles so that future mathematical models can include an improved understanding of within-host variation.

METHODS

Patients' parasitaemia data were analyzed to identify (i) clusters of patients (from 35 patients) that have a similar overall parasitaemia profile and (ii) clusters of patients (from 100 patients) that have a similar first wave of parasitaemia. For each cluster analysis, patients were clustered based on key features which previous models used to summarize parasitaemia dynamics. The clustering analyses were performed using a finite mixture model. The centroid values of the clusters were used to parameterize two established within-host models to generate parasitaemia profiles. These profiles (that used the novel centroid parameterization) were compared with profiles that used individual-specific parameterization (as in the original models), as well as profiles that ignored individual variation (using overall means for parameterization).

RESULTS

To capture the variation of within-host dynamics, when studying the overall parasitaemia profile, two clusters efficiently grouped patients based on their infection length and the height of the first parasitaemia peak. When studying the first wave of parasitaemia, five clusters efficiently grouped patients based on the height of the peak and the speed of the clearance following the peak of parasitaemia. The clusters were based on features that summarize the strength of patient innate and adaptive immune responses. Parameterizing previous within host-models based on cluster centroid values accurately predict individual patient parasitaemia profiles.

CONCLUSION

This study confirms that patients have personalized immune responses, which explains the variation of parasitaemia dynamics. Clustering can guide the optimal inclusion of within-host variation in future studies, and inform the design and parameterization of population-based models.

Mechanistic within-host models of the asexual Plasmodium falciparum infection: a review and analytical assessment

Background

Malaria blood-stage infection length and intensity are important drivers of disease and transmission; however, the underlying mechanisms of parasite growth and the host’s immune response during infection remain largely unknown. Over the last 30 years, several mechanistic mathematical models of malaria parasite within-host dynamics have been published and used in malaria transmission models.

Methods

Mechanistic within-host models of parasite dynamics were identified through a review of published literature. For a subset of these, model code was reproduced and descriptive statistics compared between the models using fitted data. Through simulation and model analysis, key features of the models were compared, including assumptions on growth, immune response components, variant switching mechanisms, and inter-individual variability.

Results

The assessed within-host malaria models generally replicate infection dynamics in malaria-naïve individuals. However, there are substantial differences between the model dynamics after disease onset, and models do not always reproduce late infection parasitaemia data used for calibration of the within host infections. Models have attempted to capture the considerable variability in parasite dynamics between individuals by including stochastic parasite multiplication rates; variant switching dynamics leading to immune escape; variable effects of the host immune responses; or via probabilistic events. For models that capture realistic length of infections, model representations of innate immunity explain early peaks in infection density that cause clinical symptoms, and model representations of antibody immune responses control the length of infection. Models differed in their assumptions concerning variant switching dynamics, reflecting uncertainty in the underlying mechanisms of variant switching revealed by recent clinical data during early infection. Overall, given the scarce availability of the biological evidence there is limited support for complex models.

Conclusions

This study suggests that much of the inter-individual variability observed in clinical malaria infections has traditionally been attributed in models to random variability, rather than mechanistic disease dynamics. Thus, it is proposed that newly developed models should assume simple immune dynamics that minimally capture mechanistic understandings and avoid over-parameterization and large stochasticity which inaccurately represent unknown disease mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03813-z.

Interactions between host immune response and antigenic variation that control Borrelia burgdorferi population dynamics

The population dynamics of pathogens within hosts result from interactions between host immune responses and mechanisms of the pathogen to evade or resist immune responses. Vertebrate hosts have evolved adaptive immune responses to eliminate the infection, while many pathogens evade immune clearance through altering surface antigens. Such interactions can result in a characteristic pattern of pathogen population dynamics within hosts consisting of population growth after infection, rapid population decline following specific immune responses, followed by persistence at low densities during a chronic infection stage. Despite the medical importance of chronic infections, little is known about the conditions of the interactions between variable antigens and the adaptive immune system that cause the characteristic pathogen population dynamics. Using the vls antigenic variation system of the Lyme disease pathogen, Borrelia burgdorferi, as a model system, we investigated conditions of the interaction between the antigenic variation system and the adaptive immune response that can explain the within-host population dynamics of B. burgdorferi using mathematical modelling. This characteristic population dynamic pattern can be explained by models that assume a variable immune removal rate of antibody-bound B. burgdorferi. However, models with a constant immune removal rate could reproduce the rapid population decline of B. burgdorferi populations but not their long-term persistence within hosts using parameter values determined by fitting empirical data. The model predictions, along with the assumptions about the interactions between B. burgdorferi and the immune response, can be tested experimentally to estimate the likelihood that each mechanism affects B. burgdorferi population dynamics in real infections.

Patterns of inflammatory responses and parasite tolerance vary with malaria transmission intensity

Background

In individuals living in malaria-endemic regions, parasitaemia thresholds for the onset of clinical symptoms vary with transmission intensity. The mechanisms that mediate this relationship are however, unclear. Since inflammatory responses to parasite infection contribute to the clinical manifestation of malaria, this study investigated inflammatory cytokine responses in children with malaria from areas of different transmission intensities (ranging from low to high).

Methods

Blood samples were obtained from children confirmed with malaria at community hospitals in three areas with differing transmission intensities. Cytokine levels were assessed using the Luminex^®-based magnetic bead array system, and levels were compared across sites using appropriate statistical tests. The relative contributions of age, gender, parasitaemia and transmission intensity on cytokine levels were investigated using multivariate regression analysis.

Results

Parasite density increased with increasing transmission intensity in children presenting to hospital with symptomatic malaria, indicating that the parasitaemia threshold for clinical malaria increases with increasing transmission intensity. Furthermore, levels of pro-inflammatory cytokines, including tumour necrosis factor alpha (TNF-α), interferon-gamma (IFN-γ), interleukin (IL)-1β, IL-2, IL-6, IL-8, and IL-12, decreased with increasing transmission intensity, and correlated significantly with parasitaemia levels in the low transmission area but not in high transmission areas. Similarly, levels of anti-inflammatory cytokines, including IL-4, IL-7, IL-10 and IL-13, decreased with increasing transmission intensity, with IL-10 showing strong correlation with parasitaemia levels in the low transmission area. Multiple linear regression analyses revealed that transmission intensity was a stronger predictor of cytokine levels than age, gender and parasitaemia.

Conclusion

Taken together, the data demonstrate a strong relationship between the prevailing transmission intensity, parasitaemia levels and the magnitude of inflammatory responses induced during clinical malaria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1796-x) contains supplementary material, which is available to authorized users.

Malaria Parasites in the Asymptomatic: Looking for the Hay in the Haystack

With malaria elimination back on the international agenda, programs face the challenge of targeting all Plasmodium infections, not only symptomatic cases. As asymptomatic individuals are unlikely to seek treatment, they are missed by passive surveillance while remaining infectious to mosquitoes, thus acting as silent reservoirs of transmission. To estimate the risk of asymptomatic infections in various phases of malaria elimination, we need a deeper understanding of the underlying mechanisms favoring carriage over disease, which may involve both pathogen and host factors. Here we review our current knowledge on the determinants leading to Plasmodium falciparum symptomless infections. Understanding the host-pathogen interactions that are most likely to affect transitions between malaria disease states could guide the development of tools to tackle asymptomatic carriers in elimination settings.

A simulation model of the within-host dynamics of Plasmodium vivax infection

Background

The benign reputation of Plasmodium vivax is at odds with the burden and severity of the disease. This reputation, combined with restricted in vitro techniques, has slowed efforts to gain an understanding of the parasite biology and interaction with its human host.

Methods

A simulation model of the within-host dynamics of P. vivax infection is described, incorporating distinctive characteristics of the parasite such as the preferential invasion of reticulocytes and hypnozoite production. The developed model is fitted using digitized time-series’ from historic neurosyphilis studies, and subsequently validated against summary statistics from a larger study of the same population. The Chesson relapse pattern was used to demonstrate the impact of released hypnozoites.

Results

The typical pattern for dynamics of the parasite population is a rapid exponential increase in the first 10 days, followed by a gradual decline. Gametocyte counts follow a similar trend, but are approximately two orders of magnitude lower. The model predicts that, on average, an infected naïve host in the absence of treatment becomes infectious 7.9 days post patency and is infectious for a mean of 34.4 days. In the absence of treatment, the effect of hypnozoite release was not apparent as newly released parasites were obscured by the existing infection.

Conclusions

The results from the model provides useful insights into the dynamics of P. vivax infection in human hosts, in particular the timing of host infectiousness and the role of the hypnozoite in perpetuating infection.

Measurement of the plasma levels of antibodies against the polymorphic vaccine candidate apical membrane antigen 1 in a malaria-exposed population

Background

Establishing antibody correlates of protection against malaria in human field studies and clinical trials requires, amongst others, an accurate estimation of antibody levels. For polymorphic antigens such as apical membrane antigen 1 (AMA1), this may be confounded by the occurrence of a large number of allelic variants in nature.

Methods

To test this hypothesis, plasma antibody levels in an age-stratified cohort of naturally exposed children from a malaria-endemic area in Southern Ghana were determined by indirect ELISA. Titres against four single PfAMA1 alleles were compared with those against three different allele mixtures presumed to have a wider repertoire of epitope specificities. Associations of antibody levels with the incidence of clinical malaria as well as with previous exposure to parasites were also examined.

Results

Antibody titres against PfAMA1 alleles generally increased with age/exposure while antibody specificity for PfAMA1 variants decreased, implying that younger children (≤ 5 years) elicit a more strain-specific antibody response compared to older children. Antibody titre measurements against the FVO and 3D7 AMA1 alleles gave the best titre estimates as these varied least in pair-wise comparisons with titres against all PfAMA1 allele mixtures. There was no association between antibody levels against any capture antigen and either clinical malaria incidence or parasite density.

Conclusions

The current data shows that levels of naturally acquired antigen-specific antibodies, especially in infants and young children, are dependent on the antigenic allele used for measurement. This may be relevant to the interpretation of antibody titre data from measurements against single PfAMA1 alleles, especially in studies involving infants and young children who have experienced fewer infections.

Artemisinin-induced parasite dormancy: a plausible mechanism for treatment failure

BACKGROUND

Artemisinin-combination therapy is a highly effective treatment for uncomplicated falciparum malaria but parasite recrudescence has been commonly reported following artemisinin (ART) monotherapy. The dormancy recovery hypothesis has been proposed to explain this phenomenon, which is different from the slower parasite clearance times reported as the first evidence of the development of ART resistance.

METHODS

In this study, an existing P. falciparum infection model is modified to incorporate the hypothesis of dormancy. Published in vitro data describing the characteristics of dormant parasites is used to explore whether dormancy alone could be responsible for the high recrudescence rates observed in field studies using monotherapy. Several treatment regimens and dormancy rates were simulated to investigate the rate of clinical and parasitological failure following treatment.

RESULTS

The model output indicates that following a single treatment with ART parasitological and clinical failures occur in up to 77% and 67% of simulations, respectively. These rates rapidly decline with repeated treatment and are sensitive to the assumed dormancy rate. The simulated parasitological and clinical treatment failure rates after 3 and 7 days of treatment are comparable to those reported from several field trials.

CONCLUSIONS

Although further studies are required to confirm dormancy in vivo, this theoretical study adds support for the hypothesis, highlighting the potential role of this parasite sub-population in treatment failure following monotherapy and reinforcing the importance of using ART in combination with other anti-malarials.

Antigenic variation in Plasmodium falciparum malaria involves a highly structured switching pattern

Many pathogenic bacteria, fungi, and protozoa achieve chronic infection through an immune evasion strategy known as antigenic variation. In the human malaria parasite Plasmodium falciparum, this involves transcriptional switching among members of the var gene family, causing parasites with different antigenic and phenotypic characteristics to appear at different times within a population. Here we use a genome-wide approach to explore this process in vitro within a set of cloned parasite populations. Our analyses reveal a non-random, highly structured switch pathway where an initially dominant transcript switches via a set of switch-intermediates either to a new dominant transcript, or back to the original. We show that this specific pathway can arise through an evolutionary conflict in which the pathogen has to optimise between safeguarding its limited antigenic repertoire and remaining capable of establishing infections in non-naïve individuals. Our results thus demonstrate a crucial role for structured switching during the early phases of infections and provide a unifying theory of antigenic variation in P. falciparum malaria as a balanced process of parasite-intrinsic switching and immune-mediated selection.

The malaria parasite Plasmodium falciparum avoids recognition and clearance by the immune system by sequentially switching between members of the var multi-gene family which encode the immunodominant surface proteins PfEMP1. However, some mechanism must exist to prevent rapid exposure of the pathogen's entire antigenic repertoire as this would quickly terminate the infection. It has previously been shown that the immune system can play an important role in orchestrating the sequential display of variants once an infection is established; however this does not explain how repertoire exhaustion is avoided in the initial phases of infection before an immune response has been established. Here we show that P. falciparum has evolved a highly structured switching pattern to prevent repertoire exhaustion in the early stages of infection without compromising the ability to establish new infections among partially immune individuals.

Interrupting malaria transmission: quantifying the impact of interventions in regions of low to moderate transmission

Malaria has been eliminated from over 40 countries with an additional 39 currently planning for, or committed to, elimination. Information on the likely impact of available interventions, and the required time, is urgently needed to help plan resource allocation. Mathematical modelling has been used to investigate the impact of various interventions; the strength of the conclusions is boosted when several models with differing formulation produce similar data. Here we predict by using an individual-based stochastic simulation model of seasonal Plasmodium falciparum transmission that transmission can be interrupted and parasite reintroductions controlled in villages of 1,000 individuals where the entomological inoculation rate is <7 infectious bites per person per year using chemotherapy and bed net strategies. Above this transmission intensity bed nets and symptomatic treatment alone were not sufficient to interrupt transmission and control the importation of malaria for at least 150 days. Our model results suggest that 1) stochastic events impact the likelihood of successfully interrupting transmission with large variability in the times required, 2) the relative reduction in morbidity caused by the interventions were age-group specific, changing over time, and 3) the post-intervention changes in morbidity were larger than the corresponding impact on transmission. These results generally agree with the conclusions from previously published models. However the model also predicted changes in parasite population structure as a result of improved treatment of symptomatic individuals; the survival probability of introduced parasites reduced leading to an increase in the prevalence of sub-patent infections in semi-immune individuals. This novel finding requires further investigation in the field because, if confirmed, such a change would have a negative impact on attempts to eliminate the disease from areas of moderate transmission.

Generation of humoral immune responses to multi-allele PfAMA1 vaccines; effect of adjuvant and number of component alleles on the breadth of response

There is increasing interest in multi-allele vaccines to overcome strain-specificity against polymorphic vaccine targets such as Apical Membrane Antigen 1 (AMA1). These have been shown to induce broad inhibitory antibodies in vitro and formed the basis for the design of three Diversity-Covering (DiCo) proteins with similar immunological effects. The antibodies produced are to epitopes that are shared between vaccine alleles and theoretically, increasing the number of component AMA1 alleles is expected to broaden the antibody response. A plateau effect could however impose a limit on the number of alleles needed to achieve the broadest specificity. Moreover, production cost and the vaccine formulation process would limit the number of component alleles. In this paper, we compare rabbit antibody responses elicited with multi-allele vaccines incorporating seven (three DiCos and four natural AMA1 alleles) and three (DiCo mix) antigens for gains in broadened specificity. We also investigate the effect of three adjuvant platforms on antigen specificity and antibody functionality. Our data confirms a broadened response after immunisation with DiCo mix in all three adjuvants. Higher antibody titres were elicited with either CoVaccine HT™ or Montanide ISA 51, resulting in similar in vitro inhibition (65–82%) of five out of six culture-adapted P. falciparum strains. The antigen binding specificities of elicited antibodies were also similar and independent of the adjuvant used or the number of vaccine component alleles. Thus neither the four extra antigens nor adjuvant had any observable benefits with respect to specificity broadening, although adjuvant choice influenced the absolute antibody levels and thus the extent of parasite inhibition. Our data confirms the feasibility and potential of multi-allele PfAMA1 formulations, and highlights the need for adjuvants with improved antibody potentiation properties for AMA1-based vaccines.

Problems with continuous-time malaria models in describing gametocytogenesis

Most mathematical models of malaria infection represent parasites as replicating continuously at a constant rate whereas in reality, malaria parasites replicate at a fixed age. The behaviour of continuous-time models when gametocytogenesis is included, in comparison to a more realistic discrete-time model that incorporates a fixed replication age was evaluated. Both the infection dynamics under gametocytogenesis and implications for predicting the amount parasites should invest into gametocytes (level of investment favoured by natural selection) are considered. It is shown that the many malaria models with constant replication rates can be represented by just 3 basic types. For these 3 types, it is then shown that under gametocytogenesis (i) in 2 cases, parasite multiplication and gametocyte production is mostly much too low, (ii) in the third, parasite multiplication and gametocyte production is mostly much too high, (iii) the effect of gametocyte investment on parasite multiplication is mostly too high, (iv) the effect of gametocyte investment on gametocyte production is nearly always too low and (v) with a simple approximation of fitness, the predicted level of gametocyte investment is mostly much too low. However, a continuous model with 48 age-compartments compares well to the discrete model. These findings are a further argument for modelling malaria infections in discrete time.

Pneumocystis murina MSG gene family and the structure of the locus associated with its transcription

Analysis of the Pneumocystis murina MSG gene family and expression-site locus showed that, as in Pneumocystis carinii, P. murina MSG genes are arranged in head-to-tail tandem arrays located on multiple chromosomes, and that a variety of MSG genes can reside at the unique P. murina expression site. Located between the P. murina expression site and attached MSG gene is a block of 132 basepairs that is also present at the beginning of MSG genes that are not at the expression site. The center of this sequence block resembles the 28 basepair CRJE of P. carinii, but the block of conserved sequence in P. murina is nearly five times longer than in P. carinii, and much shorter than in P. wakefieldiae. These data indicate that the P. murina expression-site locus has a distinct structure.

Plasmodium falciparum infection dynamics and transmission potential following treatment with sulfadoxine-pyrimethamine

OBJECTIVES

To investigate the overall efficacy of sulfadoxine-pyrimethamine (SP) treatment and the corresponding transmission potential for patients infected with SP-resistant Plasmodium falciparum.

METHODS

A mathematical model of the in-host dynamics of P. falciparum infections was used to simulate infections with parasites having different numbers of mutations in the dhfr and dhps genes and their responses to SP treatment. The treatment outcome and transmission potential of each simulated infection following SP treatment was assessed by tracking asexual parasite density and the number of days with sufficient mature gametocytes to give a > 95% probability of infecting a mosquito.

RESULTS

The results show treatment failure only occurring in patients infected with parasites having two mutations in dihydrofolate reductase (DHFR) combined with at least two mutations in dihydropteroate synthetase (DHPS) or with parasites having a triple mutation in DHFR. Highly mutated parasites (three mutations in each gene) caused the highest clinical failure rate, while moderately mutated parasites (three mutations in DHFR plus one mutation in DHPS) produced a high rate of asymptomatic parasitological failure following SP treatment. This high rate of asymptomatic recrudescence caused the transmission potential of infections with moderately resistant parasites to exceed that of highly resistant parasites.

CONCLUSIONS

The model output suggests that infection dynamics following SP treatment and the overall transmission potential are inherently linked. The combination of prolonged asymptomatic parasitaemia and increased transmission potential allows parasites having three mutations in DHFR, but fewer mutations in DHPS, to expand largely unnoticed.

Parasite immune escape: new views into host-parasite interactions

For parasites of humans and animals that rely on vectors or on sexual contact for transmission, it is particularly important that infection does not to terminate before the occurrence of the crucial event that completes its lifecycle (e.g. another mosquito bite). For chronic infection to occur, it is essential that the parasite avoids clearance by the host immune system. Much progress has been made in elucidating the immunological interactions and the molecular mechanisms involved in the process of immune evasion. Mathematical models have also been invaluable in understanding these processes, particularly in the generation of new ideas about a complex form of immune evasion known as antigenic variation whereby a major target of the host immune response is varied during the course of a single infection to avoid recognition.

Modeling the development of acquired clinical immunity to Plasmodium falciparum malaria

Individuals living in regions where malaria is endemic develop an acquired immunity to malaria which enables them to remain asymptomatic while still carrying parasites. Field studies indicate that cumulative exposure to a variety of diverse Plasmodium parasites is required for the transition from symptomatic to asymptomatic malaria. This study used a simulation model of the within-host dynamics of P. falciparum to investigate the development of acquired clinical immunity under different transmission conditions and levels of parasite diversity. Antibodies developed to P. falciparum erythrocyte membrane protein 1 (PfEMP1), a clonally variant molecule, were assumed to be a key human immunological response to P. falciparum infection, along with responses to clonally conserved but polymorphic antigens. The time to the development of clinical immunity was found to be proportional to parasite diversity and inversely proportional to transmission intensity. The effect of early termination of symptomatic infections by chemotherapy was investigated and found not to inhibit the host's ability to develop acquired immunity. However, the time required to achieve this state was approximately double that compared to when no treatment was administered. This study demonstrates that an immune response primarily targeted against PfEMP1 has the ability to reduce clinical symptoms of infections irrespective of whether treatment is administered, supporting its role in the development of acquired clinical immunity. The results also illustrate a novel use for simulation models of P. falciparum infections, investigation of the influence of intervention strategies on the development of naturally acquired clinical immunity.

Evolution of resistance to sulfadoxine-pyrimethamine in Plasmodium falciparum

The development of resistance to sulfadoxine-pyrimethamine by Plasmodium parasites is a major problem for the effective treatment of malaria, especially P. falciparum malaria. Although the molecular basis for parasite resistance is known, the factors promoting the development and transmission of these resistant parasites are less clear. This paper reports the results of a quantitative comparison of factors previously hypothesized as important for the development of drug resistance, drug dosage, time of treatment, and drug elimination half-life, with an in-host dynamics model of P. falciparum malaria in a malaria-naïve host. The results indicate that the development of drug resistance can be categorized into three stages. The first is the selection of existing parasites with genetic mutations in the dihydrofolate reductase or dihydropteroate synthetase gene. This selection is driven by the long half-life of the sulfadoxine-pyrimethamine combination. The second stage involves the selection of parasites with allelic types of higher resistance within the host during an infection. The timing of treatment relative to initiation of a specific anti-P. falciparum EMP1 immune response is an important factor during this stage, as is the treatment dosage. During the third stage, clinical treatment failure becomes prevalent as the parasites develop sufficient resistance mutations to survive therapeutic doses of the drug combination. Therefore, the model output reaffirms the importance of correct treatment of confirmed malaria cases in slowing the development of parasite resistance to sulfadoxine-pyrimethamine.