Host immune constraints on malaria transmission: insights from population biology of within-host parasites
© McQueen et al.; licensee BioMed Central Ltd. 2013
Received: 4 February 2013
Accepted: 30 May 2013
Published: 15 June 2013
Plasmodium infections trigger complex immune reactions from their hosts against several life stages of the parasite, including gametocytes. These immune responses are highly variable, depending on age, genetics, and exposure history of the host as well as species and strain of parasite. Although the effects of host antibodies that act against gamete stages in the mosquito (due to uptake in the blood meal) are well documented, the effects of host immunity upon within-host gametocytes are not as well understood. This report consists of a theoretical population biology-based analysis to determine constraints that host immunity impose upon gametocyte population growth. The details of the mathematical models used for the analysis were guided by published reports of clinical and animal studies, incorporated plausible modalities of immune reactions to parasites, and were tailored to the life cycl es of the two most widespread human malaria pathogens, Plasmodium falciparum and Plasmodium vivax.
For the same ability to bind and clear a target, the model simulations suggest that an antibody attacking immature gametocytes would tend to lower the overall density of transmissible mature gametocytes more than an antibody attacking the mature forms directly. Transmission of P. falciparum would be especially vulnerable to complete blocking by antibodies to its immature forms since its gametocytes take much longer to reach maturity than those of P. vivax. On the other hand, antibodies attacking the mature gametocytes directly would reduce the time the mature forms can linger in the host. Simulation results also suggest that varying the standard deviation in the time necessary for individual asexual parasites to develop and produce schizonts can affect the efficiency of production of transmissible gametocytes.
If mature gametocyte density determines the probability of transmission, both Plasmodium species, but especially P. falciparum, could bolster this probability through evasion or suppression of host immune responses against the immature gametocytes. However, if the long term lingering of mature gametocytes at low density in the host is also important to ensure transmission, then evasion or suppression of antibodies against the mature stages would bolster probability of transmission as well.
Although the pathology of malaria is mostly due to asexually reproducing forms of the Plasmodium parasites within their hosts [1–3], these forms cannot survive in the mosquito vectors of the diseases. Instead, by a process not fully understood [4, 5], specialized sexual forms called gametocytes arise within the host, as progeny of asexual forms, and it is these cells which propagate the infection. The mating gametes emerge within the midgut of the mosquito from gametocytes harbored in the blood meal.
The resultant zygote becomes an ookinete and then an oocyst that produces thousands of sporozoites which invade the salivary glands of the vector and thereby carry the infection to a new host when the mosquito feeds [6–8]. Since Plasmodium species infect vertebrates, gametocytes along with the asexual intrahost forms have to contend with vertebrate immune responses. Many studies have shown that antibodies elicited by intrahost stages of the parasite can interfere with the mating of the gametes in the bolus of the blood meal in the mosquito gut, thus blocking transmission of the pathogen. Production of such antibodies (“reaching” from the host into the vector) has been demonstrated during Plasmodium infections in chickens [9, 10], mice , Rhesus monkeys , and human infections with Plasmodium falciparum and Plasmodium vivax[14, 15]. Other studies have shown that monoclonal antibodies against specific surface antigens on the gametes of P. falciparum can reduce the number of oocysts in mosquitoes [16, 17].
The effects of a host’s acquired immunity upon the intrahost gametocytes before uptake by mosquitoes are less clear. A 1977 study  of eleven Gambian children who carried P. falciparum gametocytes found that four subjects had antibodies against gametocytes, but these antibodies did not interact with the surface of erythrocytes parasitized by gametocytes. The seven subjects without antibodies to gametocytes still were able to eliminate them. Since this study, though, newer evidence suggests a reconsideration of transmission blocking due to acquired immunity to eliminate or otherwise cripple gametocytes before uptake by the vector . In particular a 2008 study of Gambian children showed that 34% of them had antibodies to antigens on the surface of erythrocytes parasitized with mature gametocytes (referred to gametocyte surface antigens or GSAs) of the 3D7 strain of P. falciparum. This study also showed that those patients treated with chloroquine and/or sulphadoxine-pyrimethamine for four weeks had statistically significantly lower gametocyte density after treatment if they had antibodies to GSAs than if they did not. (Chloroquine and sulphadoxine-pyrimethamine act against asexual forms, and do not block gametocytemia in patients infected with P. falciparum parasites sensitive to these drugs ).
In addition to activating acquired immune responses, gametocytes, along with other intrahost stages, interact with host innate immune responses as well . Studies of Plasmodium cynomolgi infections in toque monkeys showed that cytokines TNF and INF- γ were needed for the killing of gametocytes . Another experiment demonstrated that white blood cells need nitric oxide to kill gametocytes of P. falciparum and P. vivax. A recent study on mice infected with Plasmodium chabaudi showed that blockage of the TNF receptor (by antibodies to this receptor) increased the transmissibility of this parasite . Finally, epidemiological evidence suggest species-dependent differences in human immune responses. Gametocytemia was found to correlate with high fever in P. vivax infection, but not in P. falciparum infection, in two studies: one of malaria patients in Peru and Thailand , and another of neurosyphilis patients undergoing malariatherapy . Thus, gametocyte-host immune system interactions remain an active area of research for many reasons.
Since much about gametocytogenesis is still uncertain, a “non-cryptic sexual” (non-CS) modality of gametocytogenesis was examined as a test of robustness of the model results. In this modality 5% of merozoites from each bursting schizont go directly into gametocyte development after parasitizing an erythrocyte with no cryptic sexual forms. Results for models that incorporated the non-CS modality were found to be very similar to results form simulations using the CS modality. For that reason, this report focuses on results from models using the CS modality, while results based on non-CS models are given in Additional file 1.
The two major Plasmodium species responsible for the bulk of public health burden of human malarias are P. falciparum and P. vivax[3, 38], which share the same basic lifecycle and development stages. Gametocytogenesis during P. falciparum infection apparently is modulated by factors in the blood as infection progresses [39–42]; for simplicity it is assumed that commitment of some merozoites to the sexual pathway starts when the density of intraerythrocytic asexual forms reaches 0.01 μL−1 for both species. Despite the similarities of P. falciparum and P. vivax, their last common ancestor was ≈108 years ago [43, 44], so they have evolved major physiological differences. The mature gametocytes of the two species have different physical appearance: those of P. falciparum have an elongated or crescent shape, while those of P. vivax tend to be round [45, 46]. Immature gametocytes of P. falciparum sequester during that period in bone marrow , and perhaps other organs  while those of P. vivax remain visible in the peripheral blood, (although possible sequestration of P. vivax in the human host remains an area of investigation) . In this report, the two major differences of interest are the time required for gametocytes to reach maturity, and age specificity of the erythrocytes attacked. The two species differ in the time required for gametocytes to mature into the forms that can be taken up by mosquitoes: gametocytes of P. vivax mature in 2-4 days, while those of P. falciparum require 6-10 days [50, 51]. P. vivax invades mainly the very youngest RBCs, the reticulocytes, while P. falciparum can attack RBCs of any age [52–54]. Thus, the mathematical models were tailored to take into account these differences between P. vivax and P. falciparum. It will be seen that the difference in maturity time can lead to very different gametocyte population dynamics for the two species. The focus of the Results and discussion section is on the characteristic effects host immune responses have upon the population dynamics of transmissible gametocytes as predicted by the model formalism.
Results and discussion
The outcomes of simulated malaria infections were determined for thousands of values of model parameters, allowing us to examine a large spectrum of behavior. All simulations incorporated the following two host immune responses: (1) a quickly acting, quickly decaying innate-like response that is triggered by merozoite production and acts against all the asexual intracellular parasites and their cryptic sexual counterparts (a “pyrogenic” response), and (2) a slowly acting, slowly decaying response to emulate an antibody response triggered by and acting against schizonts. For each immune responses, including any against gametocytes, the parasitemia levels that trigger the response and the maximum rates of target clearance were parameters varied from simulation to simulation; see “Simulating immune response dynamics” in Methods below for details.
Severe and potentially fatal cases of both P. falciparum malaria  and P. vivax malaria  manifest as respiratory distress, cerebral complications, severe anemia, and other symptoms, but for this report it is assumed that the host dies if the uninfected red blood cell counts drops below 60% of normal, either by direct parasitization of erythrocytes or dyserythropoiesis or both. This reports gives results only for simulated infections for which intracellular asexual forms exceed the threshold to trigger sexual forms (0.01 μL−1), and infection duration is less than three years.(Here, infection duration means the time from primary release of merozoites from the liver until either the clearance of all forms of the parasite from the host or death of host.) Over 65% of the simulations generated met these criteria. Of those that did, ≈5% resulted in death of host, with a median survival time of 152 days for P. falciparum and 161 days for P. vivax. Of the simulations considered for this report in which the host clears the parasite, the duration of the asexual forms in the host ranged from 2 days to over 1080 days, for both P. falciparum and P. vivax with a median ≈11 days for both species. The results of over 105 simulations were used for this report. Immune constraints on the asexual parasite population dynamics were discussed in a previous publication , so this report focuses mainly on the behavior of the gametocytes. (Additional file 1: Figures S1–S5 give a summary of the behavior of the asexual populations as a function of model parameters.)
Example time series
Parameters for the two simulated infections shown in Figure 2
Value for cleared infection
Value of fatal infection
Model results if gametocytes are invisible to acquired immunity
If the model innate response can remove gametocytes, almost all simulated P. vivax infections still make mature gametocytes if PD Asx ≥ 1 μ L−1day, although there is some reduction in F MG for P. falciparum infections if PD Asx ≥ 1 μ L−1day. For the subset of simulated infections that make mature gametocytes, the range in PD MG for a given PD Asx is somewhat larger, but PD MG still grows almost linearly with PD Asx , although at a reduced level: for P. falciparum, , and for P. vivax, . The exception is for when the host dies, especially for P. falciparum: then PD MG is suppressed by one to three orders of magnitude if the innate response affects gametocytes (versus when it does not). Investigation of model infections that end in death of host showed that the innate immune response is triggered to be active most of the time due to the failure of innate and acquired immunity combined to control the asexual forms (as for the example illustrated in Figure 2C and 2D). A useful quantity to quantify the activeness of the innate immune response is the time integral of the clearance rate of the innate response integrated over the course of infection, I χ Inn . For the P. falciparum simulations that ended in host death, I χ Inn =108.4, while I χ Inn =15.0 for simulations in which the parasite is cleared with PD Asx in the interval 104.5μ L−1day−105μ L−1day. (Note that I χ Inn is a dimensionless quantity.) Similarly, for the P. vivax simulations that ended in host death, I χ Inn = 125.6, while I χ Inn =4.53 for simulations in which the parasite is cleared with P D Asx in the interval 104μ L−1day−104.5μ L−1day. Since gametocytes for P. falciparum and P. vivax have a longer duration for development than the asexual stages, the gametocytes would be vulnerable if a sustained immune response could affect them. The case of infections ending with host death is analyzed in more detail in the Discussion section below.
Similar results were found for simulations with the non-CS models; see Additional file 1: Figures S6 and S7.
Antibodies against immature gametoctyes tend to be more effective in reducing the density of transmissible gametocytes than antibodies directly against the transmissible forms
Two modalities of acquired response against gametocytes were considered: (1) an antibody response to the mature gametocytes, and (2) an antibody response to the immature gametocytes. It is assumed that the vulnerable stage makes the triggering antigen during the entire duration of the stage. As a measure of the strength of the antibody response to the gametocyte, the quantity IS Gcy =χAb,Mx/Tar Th was defined, where χAb,Mx is the maximum clearance rate of the targeted stage, and Tar Th is the density of the targeted stage which triggers the response. See “Model for immune response dynamics” in the Methods section below.
For the model antibody modalities, the non-CS models showed very similar behavior; see Additional file 1: Figures S8 and S9.
When infections make mature gametocytes, antibodies against mature gametocytes are more effective in reducing the window of transmission than antibodies against the immature forms
Mature gametocyte production is affected by the standard deviation in the development time of the asexual forms
The most important insight from the simulations studied in this report is that differing modalities of host immune response against gametocytes would affect their population dynamics in fundamentally different ways that may have implications for transmission. First, if there is no acquired immunity against gametocytes, the density of the asexual forms determines the density of the mature gametocytes, assuming that that after some triggering event, a set proportion of merozoites develop into sexual forms. The parasite-days PD MG and maximum density Max MG of mature gametocytes grow almost linearly with the asexual parasite days PD Asx and maximum asexual density Max Asx , respectively. This is true even if the model innate response could remove gametocytes at the same clearance rate as the asexual forms, although in infections in which the host clears the parasite, PD MG and Max MG are reduced by the effects of the innate immunity. In simulated infections in which the host dies by anemia, the model innate response suppresses PD MG and Max MG by a factor of 10 to 100 in P. vivax infection, and a factor of 1000 in P. falciparum infections compared to the case of gametocytes being invisible to all immunity. However, for the model examined here, one should realized that death by severe anemia takes nearly three weeks at a minimum from primary release, and the model innate response is activated more on average than during infections in which the host clears the parasites. If innate responses are inhibited or exhausted, or if severe malaria is induced by a mechanism other than anemia with hyper-parasitemia, then the gametocyte levels in severe malaria might differ from the behavior suggested by the models. In addition, some researchers have suggested that the sequestration of immature gametocytes of P. falciparum (as discussed in the Background above) helps the immature forms evade host TNF pyrogenic response, thus moderating suppression of P. falciparum gametocytes by this mechanism of innate immunity .
For the immune modalities in which gametocytes are attacked by an antibody response, antibodies capable of clearing immature gametocytes would suppress the density of the transmissible mature forms due to “choking off” the source of future mature gametocytes. For the same value of the immune strength measure IS Gcy , such antibodies are much more effective at reducing overall numbers of the mature forms than antibodies directly against the mature forms. This was true for both Plasmodium species; in particular, P. falciparum is extremely vulnerable to having its potential transmission completely blocked by this immune modality in simulated infections due to the long time needed for its immature gametocytes to mature. Even when mature gametocytes are made, their density is suppressed much more than in P. vivax infections for the same values of PD Asx and IS Gcy . Yet clinical measurements in malaria patients suggest that the densities of mature gametocytes of the two species are comparable , which leads to a question: how does P. falciparum produce enough mature forms to ensure transmission? One possibility is that the proportion of merozoites of this species that go down the sexual pathway is so large that despite enormous depletion of the immature forms, enough of them reach maturity at a level capable of causing transmission. A survey of blood from the bone marrow from a group of P. falciparum infected Gambian children found only ≈0.05% of erythrocytes in that compartment actually contained immature gametocytes, while ≈8% of the peripheral erythrocytes were infected with asexual forms , so a colossal rate of gametocytogenesis in P. falciparum infection seems unlikely. A second way that an antibody response could be mitigated is by sequestering of mature gametocytes in peripheral capillaries, thus boosting their density in the blood that a mosquito would uptake. Experimental evidence suggests that this happens in the rodent malaria caused by P. chabaudi. Of course, some immature forms have to survive to maturity in order for this mechanism to be effective. The results above suggests a third way: perhaps the immature gametocytes of P. falciparum can evade host antibody responses by not showing surface antigens, or are vulnerable only for a fraction of their duration. A recent study found that immature gametocytes of P. falciparum do not form the “adhesive knob structures” typical seen on the surface of erythrocytes parasitize by the asexual forms of P. falciparum and express much less erythrocyte membrane protein PfEMP1 then the asexual forms . As mentioned in the introduction, a study in Gambia found antibodies to antigens on the surface of erythrocytes parasitized with mature gametocytes of P. falciparum, but not to the immature gametocytes . And there is some evidence that antigenic variation of a surface antigen may aid the survival of immature gametocytes , but more needs to be done to determine this.
A study on neurosyphilis patients treated with P. falciparum suggested that the probability that mosquitoes become infected after feeding on a patient is proportional to the gametocyte density , although the authors of that study noted that some patients with high gametocyte density failed to infect, while others with gametocyte density <10 μL−1 infected. Another study on malaria patients in Ghana did not indicate a clear correlation between gametocyte density and transmission to mosquitoes . Furthermore, a study on neurosyphilis patients with P. vivax infection suggested that transmission success depends on only the density of male gametocytes . These conflicting results suggest that while the probability that a patient can transmit the infection may be affected by PD MG there must be other factors influencing transmission. The window of time when transmission is possible, Dur MG may be just as important. The results reported above indicate that once mature gametocytes are produced, antibodies to the mature forms would be important in reducing Dur MG . Of course, such antibodies would still reduce PD MG , even if not to as large a degree as antibodies to the immature forms. Thus if a threshold density of mature gametocytes is needed to ensure transmission, or if mature gametocytes concentrate in subdermal capillaries, antibodies to the mature forms could still interfere with transmission.
Another insight from this study is that details of development of the asexual forms can combine with a host’s immune system dynamics and affect the production of mature gametocytes. Theoretical studies suggest that the presence of a quickly acting, quickly decaying innate immune response forces the synchronization of schizont rupturing , and synchronization of schizont rupture occurs in models examined here as well . In addition, the circadian melatonin cycle of the host may also regulate the overall cycle of asexual development . The synchronization of merozoite release is not directly set by σ Asx , the standard deviation in the population-averaged duration of the asexual intracellular stages, but σ Asx sets the “sharpness” of this release. However, the value of σ Asx in real infections is not clear. It was previously reported that the larger σ Asx is, the larger the density of the asexual forms, and that this effect seems to be due to the presence of a quickly-acting innate response triggered by merzoite releases . It was shown above that if the gametocytes are invisible to host immunity, PD MG /PD Asx weakly decreases with increasing σ Asx for a given PD Asx . However, if the model innate immunity could remove gametocytes, smaller σ Asx tends to suppress PD MG /PD Asx , especially for P. falciparum infections with PD Asx ≥ 300 μL−1day and σ Asx ≤ 0.7 hr. If the rupturing of schizonts releases immune suppression factors , and if these factors help suppress damage to gametocytes, then the mass release of such factors when schizonts rupture in unison would preserve the efficiency of gametocyte production enough to ensure transmission even if σ Asx ≤ 1 hr. Synchronization of schizont rupture and release of immune suppression factors during rupture may be processes evolutionarily shaped by intense host-parasite interactions.
In order to elucidate the consequences of different host immune modalities against the gametocytes, this report consists of a theoretical population biological study of Plasmodium populations interacting with their human hosts. It is shown that antibodies against the immature gametocytes would suppress the density of mature gametocytes to a greater degree than antibodies that act directly against the mature forms for the same ability of antibody to clear their targets. This would be true for both P. vivax and P. falciparum infections; in particular, P. falciparum is very vulnerable to having its potential transmission completely blocked by this immune modality. Antibodies directly against mature forms would still interfere with transmission by reducing the window of transmissibility during and after a bout of malaria. In addition to the direct effects of immune modalities upon gametocytes, it was also found that the standard deviation in the development time of individual asexual intraerythrocytic parasites can affect the efficiency at which tranmissible gametocytes are produced due to the interplay between the development of the asexual forms and host innate immunity.
Basic algorithm for modeling dynamics of aging populations
In this formalism, the total population is given by the sum of all P n , and s is the rate (source term) at which new individuals enter the population. The tangible quantities D and σ set N, the total number of components. (If D=σ, then the system of Equations 1 reduces to a single equation describing exponential decay.) As mentioned above, this formalism has been used to describe within-host populations of the asexual Plasmodium cells interacting with various systems of the host [29–31]. For this report, this formalism has been adapted to include gametocytes as well. The populations considered in this report are shown schematically in Figure 1, along with the corresponding values of D and σ. (All population sizes are stated as number per μ L of blood.) Because the full system of differential equations for the modal formalism involves many hundreds of components, it is given in detail in Additional file 2. A qualitative description is given below for each of the major parts of the model.
Model for asexual parasite population dynamics
As shown in Figure 1, five populations of asexual parasite cells which are morphologically distinct from each other were considered: (1) ring stage, (2) early trophozoites, (3) late trophozoites, (4) schizonts, and (5) merozoites. The dynamics of each of these populations is described by a set of equations similar to Eq. 1 with the addition of terms due to host immune responses described below. (See equation 6 in Additional file 2). The duration and variance of each of the intracellular stages (1-4) were taken to be the same, with duration D Asx = 12 hr, but with the variance varied from simulation to simulation. (The sampling of parameters varied from simulation to simulation as explained and tabulated in Subsection “Sampling the parameter space” below). Since merozoites are short lived in blood, with duration D μ =0.1 hr, just one compartment was used for them, μ, and take σ μ =D μ . There are two sources for the asexual population: (1) the primary release of merozoites from the liver, thus triggering the blood infection phase of the disease, the (2) merozoite that are released from the bursting schizonts, with p number of merozoites released on average from each schizont. The simulated infection begins with the primary release of merozoites, an event that takes just a few hours to complete . For every simulated infection, it was assumed that the rate of primary release was 0.002 (μL hr)−1, maintained for the one hour. This is equivalent to 104 merozoites being released into an adult human with blood volume 5×106μL. Whether from primary release or from bursting schizonts, the rate of creation for new trophozoites is μ ζ V, where μ is the density of merozoites, ζ is the binding affinity of merozoites to their target erythrocyte, and V is the density of the targeted erythroyctes (reticulocytes only for P. vivax, all red blood cells for P. falciparam). (Remember that μ is also time dependent.) For this report, asexual density Asx means the sum of the ring stage, early and late trophozoite, and schizont blood densities. (Although the later stages of the intracellular asexual forms sequester in P. falciparum infection, Asx was taken as a measure of the full load of the asexual forms on the host.) The parasite days for the asexuals mentioned in the discussion, PD Asx , would be the average of Asx during infection in host times the duration of infection in host. IfAsx ≥ 0.01 μL−1, then it is assumed that a fraction r=0.05 of new infected erythrocytes are committed to the cryptic sexual pathway.
(Here V0 is the initial density of vulnerable red blood cells.) Observations of parasite growth in neurosyphilis patients treated with malaria therapy [64, 73] as well as in inoculated volunteers  indicate that R0≈15 for both P. vivax and P. falciparum. Photographs of bursting P. falciparum schizonts indicate p ≥ 16 . Thus, for this report p = 16, R0=15. Equation 2 fixes ζ to 3×10−5μLhr−1 for P. falciparum and 2.4×10−3μlhr−1 for P. vivax.
It is assumed that the four intracellular stages are attacked by the innate immune response of the host, with time-dependent clearance rate χ Inn . Host acquired responses in malaria are complicated, (see  and  for review), but for simplicity, it is assumed for this report that the main antibody response is against the schizont stage with time-dependent clearance rate χSch,Ab. (Immune dynamics models are explained in more detail in Subsection “Model for immune response dynamics” below.) Also for simplicity, the variation in surface antigens which P. falciparum employs for immune evasion was not considered .
Model for population dynamics of sexual forms
As mentioned above, two very different models of gametocytogenesis were considered, the “cryptic sexual” (CS) model and the “non-cryptic sexual” (non-CS) model. As illustrated in Figure 1, in the CS model there is a set of parasite stages that morphologically resemble the asexual stages, but are committed to eventual development into overt gametocytes. The ring stage, early trophozoites, late trophozoites, schizonts, and merozoites each have cryptic sexual counterparts. (See equation 7 in Additional file 2.) An erythrocyte invaded by a cryptic sexual merozoite becomes a gametocyte. For each simulation, it was assumed that the values corresponding to D Asx , σ Asx , N Asx , D μ , p, and ζ are the same as for the asexual populations. Based on studies of gametocyte duration [50, 51]. the immature gametocyte duration D IG was set to 216hr for P. falciparum, with σ IG = 24 hr. For P. vivax, D IG was set to 72 hr with σ IG = 12 hr. For simplicity, the mature gametocyte population was represented with a single compartment, with exponential decay, D MG = σ MG = 156hr, based on a study of gametocyte dynamics in malaria therapy patients . (The actual point at which a cryptic sexual form emerges is not yet clear, but the CS model is a extreme formulation of the cryptic sexual thesis.) It was assumed that the cryptic sexual forms are subject to the same innate clearance rate χ Inn and antibody clearance rate χSc,Ab. In the Non-CS model, there are no cryptic sexual stages at all; a certain proportion (r=0.05) of erythrocytes invaded by regular merozoites directly become gametocytes once Asx ≥ 0.01 μL−1 (See equation 8 in Additional file 2.)
Four immune modalities were considered for gametocytes in this report for both P. falciparum and P. vivax, and for both the CS and Non-CS models: (1) no host immune responses to gametocytes at all, (2) no antibody response to gametocytes, but an innate response with the same clearance rate χ Inn as for the asexuals, (3) an innate response to all gametocytes with clearance rate χ Inn , and an antibody response to the immature gametocytes but not to mature gametocytes (and different from the one to schizonts) with clearance rate χIG,Ab, and (4) an innate response to all gametocytes with clearance rate χ Inn , and an antibody response on the mature but not immature gametocytes, with clearance rate χMG,A b. The parameters for antibody responses to gametocytes are chosen independently to the response against schizonts. See “Model for immune response dynamics” below.
Model for red blood cell dynamics
Three populations were used to describe the red blood cells: (1) reticulocytes, the youngest of the erythrocytes, (2) mature red blood cells, and (3) senescent red blood cells ready to be removed by phagocytosis in the spleen, liver or bone marrow . (See equation 9 in Additional file 2.) Based on hematological studies , the respective durations of the stages are taken as D Re =36 hr (with σ Re =6 hr),), D Ma =2796 hr (with σ Ma =148 hr), and D Se =48 hr (with σ Se =12 hr). As mentioned above, if V is the density of the erythrocyte stage or stages vulnerable to the merozoite invasion, then the total number of erythrocytes loss to infection by merozoites is ζ μ V. The basal erythrocyte density for a healthy adult was assumed to be 5×106μL−1, and if the total density of uninfected erythrocytes drop to under 60% of this value, it was assumed that the host died.
The rate of production of new reticulocytes from the bone marrow, , is the source term for the erythrocytes and has its own dynamics. (See equation 10 in Additional file 2). In a healthy human, the erythropoietic system makes new cells at a basal rate to maintain the red blood cell density at 5×106μ L−1, (≈1736 (μ L hr)−1). In the event of blood loss, the erythropoietic system can boost the rate of RBC production up to within a couple of days . However, both P. vivax and P. falciparum can modulate erythropoiesis, causing a dysfunction of erythropoiesis  probably through several mechanisms [81–84]. To account for these effects, erythropoietic source dynamics of the model were designed to be controlled by feedback from the dynamics of the RBC populations, so if there were no dyserythropoeisis, the rate of increase itself is proportional to λ E S ×ζ μ V. Here . To account for dyserythropoiesis, an offset δ Dys ζ μ V was subtracted from the growth rate in . For simplicity, it was assumed that the factor δ Dys had no dynamics of its own although it varied form from simulation to simulation; see Subsection “Sampling the parameter space” below. A value δ Dys =0 means there is no dyserythropoiesis.
Model for immune response dynamics
As indicated in Figure 1, host immune responses were modeled in a very phenomenological manner: the presence of some stage of the parasite triggers the actuator stage, which is self-amplifying and eventually triggers the attacker component that removes some stage of the parasite. The entire response becomes self-limiting. The equation of dynamics are modifications of Eq. 1 above which incorporate the self-amplification and self-limiting. Although a highly simplified model of real immune processes, it does capture real aspects of human immune reactions in malaria , and is similar to models used to describe immune responses against viruses . Up to three immune responses were considered for each simulation: (1) a quickly acting, quickly decaying innate response, (2) a slow-to-start but long lasting antibody response to intracellular asexual parasites, (3) a slow-to-start but long lasting antibody response to gametocytes. (As discussed below, in some simulations gametocytes were invisible to the host immune systems.) Consider the innate response first: the actuator is present at some background level when the response is quiet. Various studies indicate that malarias can trigger a rapid and strong cytokine response in the host upon bursting of schizonts, which then decays on a time scale of ≈1 h r until triggered again [22, 86–88]. Thus, in models reported here the presence of merozoites triggers the production of more actuator, which in turn triggers the attacker component that clears its targets at rate χ Inn . A buildup of or χ Inn then limits the response. The limit on growth of the actuator is set by parameter . The limit on growth of the attacker is set by maximum clearance rate χInn,M x, which was varied from simulation to simulation. The level of merozoites that triggers the model innate response was varied from simulation to simulation; see Subsection “Sampling the parameter space” below. The values of D and σ used for the actuator and attacker components, as well as the self-amplification and self-limiting terms were chosen to emulate the dynamics seen in the cytokine responses reported in malaria patients [87, 88]. (See equation 11 in Additional file 2.)
The actuator-attack formalism was adapted to model host antibody responses as well, since antibody responses also have feedback . Acquired immunity responses take some days to activate once the triggering antigen is sensed, so a delay stage was included in the antibody models with its own D and σ. The values of the D and σ used for actuator and delay stages (as shown in Figure 1) as well as the self-amplification and self-limiting terms were chosen to emulate the T-cell dynamics seen acquired immunity responses . (See equation 12 in Additional file 2.) Every simulation incorporated an antibody response to the schizont stage of the asexuals. As discussed above, some simulations incorporated an additional antibody response to the gametocytes, either against the immature forms or the mature forms. For each antibody response, the density of the targeted parasite stage, and the maximum clearance rate of the response were parameters varied from simulation to simulation. In addition, studies indicate that antibodies to malaria antigens can last from a couple of months to years [90, 91], so D for the attacker stage of the antibodies was varied from simulation to simulation as well; see Subsection“Sampling the parameter space” below. For simplicity, the effects of memory B cells were not considered.
Sampling the parameter space
Several model parameters values were varied from simulation to simulation, either because the parameters vary strongly from patient to patient, or else their values are not known. For a full listing of the parameters that were varied from simulation to simulation, the range of values used, and the equations in which they appear, see Table S1 in Additional file 2. For a given class of sexual pathway, species, and immune structure, the Latin hypercube algorithm  was used to sample among the relevant parameters in the following manner: let be the parameters varied for given model class. The plausible range for each P n is divided into ten equal intervals. A 10×M matrix M was defined in which the columns consist of a random ordering of the integers 1 through 10, with no integer repeated in a column. The integers are associated with the parameter intervals as follows: integer k=Mi,n labels interval k for parameter P n . The first simulation uses values of chosen randomly within the intervals labeled by . The second simulation uses values chosen randomly within the intervals labeled by . This procedure is repeated until values in the intervals labeled by are used. The order of the integers in each column is scrambled to repeat the procedure again. Thus, for a given class of sexual pathway, species, and immune structure, 2000 randomized versions of M were used so that there would be 20000 simulations altogether. The Latin hypercube algorithm was used to attempt uniform sampling of the parameter space for a class of models, although with so many variable parameters the sampling could not be perfect.
Solving the model equations
To solve the ODE system, a fifth-order Runge-Kutta-Fehlberg algorithm with adaptive step-size control was used for time integration [93, 94]so that the difference between the fourth- and fifth-order solutions for each component of the ODE systems was less than one part in 106. Since cells are discrete entities the following constraint was enforced: if any density of a cell population indicates less than one cell in a normal blood volume of 5×106μ l, then all components associated with that population are set to zero. See Additional file 2 for more details on how this was enforced. A simulation stopped when (1) the host died, (2) all parasite stages are cleared, (3) simulated time reached 3 years, or (4), the adaptive time stepping algorithm itself could not converge with the desired precision. With the choice of parameters used, only for P. vivax CS models did (4) occur, affecting 39 out of 8×104 simulations.
The study involved no experimental research on humans or animals.
This study utilized the high-performance computational capabilities of the Biowulf Linux cluster at the National Institutes of Health, Bethesda, MD http://biowulf.nih.gov. PGM was supported by the Intramural Research Program of the NIH, Center for Information Technology. KCW was supported by the Intramural Research Program of the NIH, National Institute of Allergy and Infectious Diseases.
- Miller L, Baruch D, Marsh K, Doumbo O: The pathogenic basis of malaria. Nature. 2002, 415: 673-679. 10.1038/415673a.View ArticlePubMed
- English M, Newton C: Malaria: pathogenicity and disease. Chem Immunol. 2002, 80: 50-69.View ArticlePubMed
- White N, Breman J: Chapter 210: Malaria. Harrison’s Principles of Internal Medicine. Edited by: Loscalzo J, Jameson J, Hauser S, Kasper D, Fauci A, Longo D, Longo D, Fauci A, Kasper D, Hauser S, Jameson J, Loscalzo J. 2010, New York: McGraw-Hill, [http://www.accessmedicine.com]
- Talman A, Domarle O, McKenzie F, Ariey F, Robert V: Gametocytogenesis: the puberty of Plasmodium falciparum. Malar J. 2004, 3: 24-10.1186/1475-2875-3-24.PubMed CentralView ArticlePubMed
- Alano P: Plasmodium falciparum gametocytes: still many secrets of a hidden life. Mol Microbiol. 2007, 66: 291-302. 10.1111/j.1365-2958.2007.05904.x.View ArticlePubMed
- Sinden R: Sexual development of malarial parasites. Adv Parasitol. 1983, 22: 153-216.View ArticlePubMed
- Billker O, Lindo V, Panico M, Etienne A, Paxton T, Dell A, Rogers M, Sinden R, Morris H: Identification of xanthurenic acid as the putative inducer of malaria development in the mosquito. Nature. 1998, 392: 289-292. 10.1038/32667.View ArticlePubMed
- Alavi Y, Arai M, Mendoza J, Tufet-Bayona M, Sinha R, Fowler K, Billker O, Franke-Fayard B, Janse C, Waters A, Sinden R: The dynamics of interactions between Plasmodium and the mosquito: a study of the infectivity of Plasmodium berghei and Plasmodium gallinaceum, and their transmission by Anopheles stephensi, Anopheles gambiae and Aedes aegypti. Int J Parasitol. 2003, 33: 933-943. 10.1016/S0020-7519(03)00112-7.View ArticlePubMed
- Carter R, Chen D: Malaria transmission blocked by immunisation with gametes of the malaria parasite. Nature. 1976, 263: 57-60. 10.1038/263057a0.View ArticlePubMed
- Gwadz R: Successful immunization against the sexual stages of Plasmodium gallinaceum. Science. 1976, 193: 1150-1151. 10.1126/science.959832.View ArticlePubMed
- Mendis K, Targett G: Immunisation against gametes and asexual erythrocytic stages of a rodent malaria parasite. Nature. 1979, 277: 389-391. 10.1038/277389a0.View ArticlePubMed
- Gwadz R, Green I: Malaria immunization in Rhesus monkeys. A vaccine effective against both the sexual and asexual stages of Plasmodium knowlesi. J Exp Med. 1978, 148: 1311-1323. 10.1084/jem.148.5.1311.View ArticlePubMed
- Graves P, Carter R, Burkot T, Quakyi I, Kumar N: Antibodies to Plasmodium faciparum gamete surface antigens in Papua New Guinea sera. Parasite Immunol. 1988, 10: 209-218. 10.1111/j.1365-3024.1988.tb00215.x.View ArticlePubMed
- Mendis K, Munesinghe Y, de Silva Y, Keragalla I, Carter R: Malaria transmission-blocking immunity induced by natural infections of Plasmodium vivax in humans. Infect Immun. 1987, 55: 369-372.PubMed CentralPubMed
- Snewin V, Premawansa S, Kapilananda G, Ratnayaka L, Udagama P, Mattei D, Khouri E, Del Giudice G, Peiris J, Mendis K, David P: Transmission blocking immunity in Plasmodium vivax malaria: antibodies raised against a peptide block parasite development in the mosquito vector. J Exp Med. 1995, 181: 357-362. 10.1084/jem.181.1.357.View ArticlePubMed
- Rener J, Graves P, Carter R, Williams J, Burkot T: Target antigens of transmission-blocking immunity on gametes of Plasmodium falciparum. J Exp Med. 1983, 158: 976-981. 10.1084/jem.158.3.976.View ArticlePubMed
- Vermeulen A, Ponnudurai T, Beckers P, Verhave J, Smits M, Meuwissen J: Sequential expression of antigens on sexual stages of Plasmodium falciparum accessible to transmission-blocking antibodies in the mosquito. J Exp Med. 1985, 162: 1460-1476. 10.1084/jem.162.5.1460.View ArticlePubMed
- Smalley M, Sinden R: Plasmodium falciparum gametocytes: their longevity and infectivity. Parasitol. 1977, 74: 1-8. 10.1017/S0031182000047478.View Article
- Sutherland C: Surface antigens of Plasmodium falciparum gametocytes–a new class of transmission-blocking vaccine targets?. Mol Biochem Parasitol. 2009, 166: 93-98. 10.1016/j.molbiopara.2009.03.007.View ArticlePubMed
- Saeed M, Roeffen W, Alexander N, Drakeley C, Targett G, Sutherland C: Plasmodium falciparum antigens on the surface of the gametocyte-infected erythrocyte. PLoS One. 2008, 3: e2280-10.1371/journal.pone.0002280.PubMed CentralView ArticlePubMed
- Hogh B, Thompson R, Hetze C, Fleck S, Kruse N, Jones I, Dgedge M, Barreto J, Sinden R: Specific and nonspecific responses to Plasmodium falciparum blood-stage parasites and observations on the gametocytemia in schoolchildren living in a malaria-endemic area of Mozambique. Am J Trop Med Hyg. 1995, 52: 50-59.PubMed
- Mendis K, Naotunne T, Karunaweera N, Del Giudice G, Grau G, Carter R: Anti-parasite effects of cytokines in malaria. Immunol Lett. 1990, 25: 217-220. 10.1016/0165-2478(90)90118-A.View ArticlePubMed
- Naotunne T, Karunaweera N, Del Giudice G, Kularatne M, Grau G, Carter R, Mendis K: Cytokines kill malaria parasites during infection crisis: extracellular complementary factors are essential. J Exp Med. 1991, 173: 523-529. 10.1084/jem.173.3.523.View ArticlePubMed
- Naotunne T, Karunaweera N, Mendis K, Carter R: Cytokine-mediated inactivation of malarial gametocytes is dependent on the presence of white blood cells and involves reactive nitrogen intermediates. Immunol. 1993, 78: 555-562.
- Long G, Chan B, Allen J, Read A, Graham A: Blockade of TNF receptor 1 reduces disease severity but increases parasite transmission during Plasmodium chabaudi chabaudi infection. Int J Parasitol. 2008, 38: 1073-1081. 10.1016/j.ijpara.2007.12.001.View ArticlePubMed
- McKenzie F, Wongsrichanalai C, Magill A, Forney J, Permpanich B, Lucas C, Erhart L, O’Meara W, Smith D, Sirichaisinthop J, Gasser R: Gametocytemia in Plasmodium vivax and Plasmodium falciparum infections. J Parasitol. 2006, 92: 1281-1285. 10.1645/GE-911R.1.PubMed CentralView ArticlePubMed
- McKenzie F, Jeffery G, Collins W: Gametocytemia and fever in human malaria infections. J Parasitol. 2007, 93: 627-633. 10.1645/GE-1052R.1.PubMed CentralView ArticlePubMed
- Wickramasinghe S, Abdalla S: Blood and bone marrow changes in malaria. Baillieres Best Pract Res Clin Haematol. 2000, 13: 277-299. 10.1053/beha.1999.0072.View ArticlePubMed
- McQueen P, McKenzie F: Age-structured red blood cell susceptibility and the dynamics of malaria infections. Proc Natl Acad Sci USA. 2004, 101: 9161-9166. 10.1073/pnas.0308256101.PubMed CentralView ArticlePubMed
- McQueen P, McKenzie F: Host control of malaria infections: constraints on immune and erythropoeitic response kinetics. PLoS Comp Biol. 2008, 4: e100149-View Article
- McQueen P: Population dynamics of a pathogen: the conundrum of vivax malaria. Biophys Rev. 2010, 2: 111-120. 10.1007/s12551-010-0034-3.PubMed CentralView ArticlePubMed
- Inselburg J: Gametocyte formation by the progeny of single Plasmodium falciparum schizonts. J Parasitol. 1983, 69: 584-591. 10.2307/3281374.View ArticlePubMed
- Bruce M, Alano P, Duthie S, Carter R: Commitment of the malaria parasite Plasmodium falciparum to sexual and asexual development. Parasitol. 1990, 100: 191-200. 10.1017/S0031182000061199.View Article
- Eksi S, Haile Y, Furuya T, Ma L, Su X, Williamson K: Identification of a subtelomeric gene family expressed during the asexual-sexual stage transition in Plasmodium falciparum. Mol Biochem Parasitol. 2005, 143: 90-99. 10.1016/j.molbiopara.2005.05.010.View ArticlePubMed
- Eksi S, Suri A, Williamson K: Sex- and stage-specific reporter gene expression in Plasmodium falciparum. Mol Biochem Parasitol. 2008, 160: 148-151. 10.1016/j.molbiopara.2008.04.005.PubMed CentralView ArticlePubMed
- Smith T, Lourenço P, Carter R, Walliker D, Ranford-Cartwright L: Commitment to sexual differentiation in the human malaria parasite, Plasmodium falciparum. Parasitol. 2000, 121: 127-133. 10.1017/S0031182099006265.View Article
- Silvestrini F, Alano P, Williams J: Commitment to the production of male and female gametocytes in the human malaria parasite Plasmodium falciparum. Parasitol. 2000, 121: 465-471. 10.1017/S0031182099006691.View Article
- Mueller I, Galinski M, Baird J, Carlton J, Kochar D, Alonso P, del Portillo HA: Key gaps in the knowledge of Plasmodium vivax, a neglected human malaria parasite. Lancet Infect Dis. 2009, 9: 555-566. 10.1016/S1473-3099(09)70177-X.View ArticlePubMed
- Smalley M: Plasmodium falciparum gametocytogenesis in vitro. Nature. 1976, 264: 271-272. 10.1038/264271a0.View ArticlePubMed
- Carter R, Miller L: Evidence for environmental modulation of gametocytogenesis in Plasmodium falciparum in continuous culture. Bull World Health Organ. 1979, 57 (Suppl 1): 37-52.PubMed CentralPubMed
- Smalley M, Brown J: Plasmodium falciparum gametocytogenesis stimulated by lymphocytes and serum from infected Gambian children. Trans R Soc Trop Med Hyg. 1981, 75: 316-317. 10.1016/0035-9203(81)90348-5.View ArticlePubMed
- Williams J: Stimulation of Plasmodium falciparum gametocytogenesis by conditioned medium from parasite cultures. Am J Trop Med Hyg. 1999, 60: 7-13.PubMed
- Carter R: Speculations on the origins of Plasmodium vivax malaria. Trends Parasitol. 2003, 19: 214-219. 10.1016/S1471-4922(03)00070-9.View ArticlePubMed
- Duval L, Fourment M, Nerrienet E, Rousset D, Sadeuh S, Goodman S, Andriaholinirina N, Randrianarivelojosia M, Paul R, Robert V, Ayala F, Ariey F: African apes as reservoirs of Plasmodium falciparum and the origin and diversification of the Laverania subgenus. Proc Natl Acad Sci USA. 2010, 107: 10561-10566. 10.1073/pnas.1005435107.PubMed CentralView ArticlePubMed
- Brumpt E: The Human parasites of the genus Plasmodium. Malariology. Edited by: Boyd M. 1949, Philadelphia: WB Saunders Company, 65-121.
- Garnham P: Malaria parasites and Other Haemosporidia. 1966, Oxford: Blackwell
- Smalley M, Abdalla S, Brown J: The distribution of Plasmodium falciparum in the peripheral blood and bone marrow of Gambian children. Trans R Soc Trop Med Hyg. 1980, 75: 103-105.View Article
- Thomson J, Robertson A: The structure and development of Plasmodium falciparum gametocytes in the internal organs and peripheral circulation. Trans R Soc Trop Med Hyg. 1935, 29: 31-40. 10.1016/S0035-9203(35)90015-3.View Article
- Anstey N, Russell B, Yeo T, Price R: The pathophysiology of vivax malaria. Trends Parasitol. 2009, 25: 220-227. 10.1016/j.pt.2009.02.003.View ArticlePubMed
- James S: Epidemiological results of a laboratory study of malaria in England. Tr Roy Soc Trop Med Hyg. 1926, 20: 143-165. 10.1016/S0035-9203(05)80002-1.View Article
- Boyd M, Kitchen S: On the infectiousness of patients infected with Plasmodium vivax and Plasmodium falciparum. Am J Trop Med. 1936, 17: 253-262.
- Craik R: A note on the erythrocytes in malaria. Lancet. 1920, 195: 1110-View Article
- Kitchen S: The infection of reticulocytes by Plasmodium vivax. Am J Trop Med Hyg. 1938, 18: 347-359.
- Shushan M, Blitz O, Adams C: The role of reticulocytes in malaria. J Lab Clin Med. 1937, 22: 364-370.
- Kochar D, Das A, Kochar S, Saxena V, Sirohi P, Garg S, Kochar A, Khatri M, Gupta V: Severe Plasmodium vivax malaria: a report on serial cases from Bikaner in northwestern India. Am J Trop Med Hyg. 2009, 80: 194-198.PubMed
- Eichner M, Diebner H, Molineaux L, Collins W, Jeffery G, Dietz K: Genesis, sequestration and survival of Plasmodium falciparum gametocytes: parameter estimates from fitting a model to malariatherapy data. Trans R Soc Trop Med Hyg. 2001, 95: 497-501. 10.1016/S0035-9203(01)90016-1.View ArticlePubMed
- Drakeley C, Sutherland C, Bousema J, Sauerwein R, Targett G: The epidemiology of Plasmodium falciparum gametocytes: weapons of mass dispersion. Trends Parasitol. 2006, 22: 424-430. 10.1016/j.pt.2006.07.001.View ArticlePubMed
- Taylor L, Read A: Why so few transmission stages? Reproductive restraint by malaria parasites. Parasitol Today. 1997, 13: 135-140. 10.1016/S0169-4758(97)89810-9.View ArticlePubMed
- Gautret P, Miltgen F, Gantier J, Chabaud A, Landau I: Enhanced gametocyte formation by Plasmodium chabaudi in immature erythrocytes: pattern of production, sequestration, and infectivity to mosquitoes. J Parasitol. 1996, 82: 900-906. 10.2307/3284196.View ArticlePubMed
- Tibúrcio M, Silvestrini F, Bertuccini L, Sander A, Turner L, Lavstsen T, Alano P: Early gametocytes of the malaria parasite Plasmodium falciparum specifically remodel the adhesive properties of infected erythrocyte surface. Cell Microbiol. 2012,http://dx.doi.org/10.1111/cmi.12062,
- Piper K, Hayward R, Cox M, Day K: Malaria transmission and naturally acquired immunity to PfEMP-1. Infect Immun. 1999, 67: 6369-6374.PubMed CentralPubMed
- Jeffery G, Eyles D: Infectivity to mosquitoes of Plasmodium falciparum as related to gametocyte density and duration of infection. Am J Trop Med Hyg. 1955, 4: 781-789.PubMed
- Muirhead-Thomson R: Low gametocyte thresholds of infection of Anopheles with Plasmodium falciparum; a significant factor in malaria epidemiology. Br Med J. 1954, 1: 68-70. 10.1136/bmj.1.4853.68.PubMed CentralView ArticlePubMed
- McKenzie F, Jeffery G, Collins W: Plasmodium vivax blood-stage dynamics. J Parasitol. 2002, 88: 521-535.PubMed CentralView ArticlePubMed
- Kwiatkowksi D, Nowak M: Periodic and chaotic host-parasite interactions in human malaria. Proc Natl Acad Sci USA. 1991, 88: 5111-5113. 10.1073/pnas.88.12.5111.View Article
- Hotta C, Gazarini M, Beraldo F, Varotti F, Lopes C, Markus R, Pozzan T, Garcia C: Calcium-dependent modulation by melatonin of the circadian rhythm in malarial parasites. Nat Cell Biol. 2000, 2: 466-468. 10.1038/35017112.View ArticlePubMed
- Kubata B, Eguchi N, Urade Y, Ymashita K, Mitamure T, Kumiko T, Osamu H, Toshihiro H: Plasmodium falciparum produces prostaglandins that are pyrogenic, somnogenic, and immunosuppressive substances in humans. J Exp Med. 1998, 188: 1197-1202. 10.1084/jem.188.6.1197.View Article
- MacDonald N: Biological Delay Systems: Linear Stability Theory. 1989, Cambridge: Cambridge University Press
- Lloyd A: The dependence of viral parameter estimates on the assumed viral life cycle: limitations of studies of viral load data. Proc R Soc B Biol Sci. 2001, 268: 847-854. 10.1098/rspb.2000.1572.View Article
- Lloyd A: Destabilization of epidemic models with the inclusion of realistic distributions of infectious periods. Proc R Soc B Biol Sci. 2001, 268: 985-993. 10.1098/rspb.2001.1599.View Article
- Johnson J, Epstein N, Shiroishi T, Miller L: Factors affecting the ability of isolated Plasmodium knowlesi merozoites to attach to and invade erythrocytes. Parasitol. 1980, 80: 539-550. 10.1017/S0031182000000998.View Article
- Baer K, Klotz C, Kappe S, Schnieder T, Frevert U: Release of hepatic Plasmodium yoelii merozoites into the pulmonary microvasculature. PLoS Pathog. 2007, 3: e171-10.1371/journal.ppat.0030171.PubMed CentralView ArticlePubMed
- Collins W, Jeffery G: A retrospective examination of sporozoite- and trophozoite-induced infections with Plasmodium falciparum: development of parasitologic and clinical immunity during primary infection. Am J Trop Med Hyg. 1999, 61 (Suppl 1): 4-19.View ArticlePubMed
- Fairley N: Sidelights on malaria in man obtained by sub-inoculation experiments. Trans R Soc Trop Med Hyg. 1947, 40: 621-676. 10.1016/0035-9203(47)90025-4.View ArticlePubMed
- Glushakova S, Yin D, Li T, Zimmerberg J: Membrane transformation during malaria parasite release from human red blood cells. Curr Biol. 2005, 15: 1645-1650. 10.1016/j.cub.2005.07.067.View ArticlePubMed
- Engwerda C, Good M: Interactions between malaria parasites and the host immune system. Curr Opin Immunol. 2005, 17: 381-387. 10.1016/j.coi.2005.05.010.View ArticlePubMed
- Plebanski M, Hill A: The immunology of malaria infection. Curr Opin Immunol. 2000, 12: 437-441. 10.1016/S0952-7915(00)00117-5.View ArticlePubMed
- Hommel M, David P, Oligino L: Surface alterations of erythrocytes in Plasmodium falciparum malaria. Antigenic variation, antigenic diversity, and the role of the spleen. J Exp Med. 1983, 157: 1137-1148. 10.1084/jem.157.4.1137.View ArticlePubMed
- Knutson M, Wessling-Resnick M: Iron metabolism in the reticuloendothelial system. Crit Rev Biochem Mol Biol. 2003, 38: 61-88. 10.1080/713609210.View ArticlePubMed
- Rapaport S: Introduction to Hematology. 1987, Philadelphia: J.B. Lippincott
- Chang K, Tam M, Stevenson M: Inappropriately low reticulocytosis in severe malarial anemia correlates with suppression in the development of late erythroid precursors. Blood. 2004, 103: 3727-3735. 10.1182/blood-2003-08-2887.View ArticlePubMed
- Layez C, Nogueira P, Combes V, Costa F, Juhan-Vague I, da Silva L, Gysin J: Plasmodium falciparum rhoptry protein RSP2 triggers destruction of the erythroid lineage. Blood. 2005, 106: 3632-3638. 10.1182/blood-2005-04-1574.View ArticlePubMed
- Casals-Pascual C, Kai O, Cheung J, Williams S, Lowe B, Nyanoti M, Williams T, Maitland K, Molyneux M, Newton C, Peshu N, Watt S, Roberts D: Suppression of erythropoiesis in malarial anemia is associated with hemozoin in vitro and in vivo. Blood. 2006, 108: 2569-2577. 10.1182/blood-2006-05-018697.View ArticlePubMed
- Ru Y, Mao B, Zhang F, Pang T, Zhao S, Liu J, Wickramasinghe S: Invasion of erythroblasts by Pasmodium vivax: A new mechanism contributing to malarial anemia. Ultrastruct Pathol. 2009, 33: 236-242. 10.3109/01913120903251643.View ArticlePubMed
- de Boer R, Oprea M, Antia R, Murali-Krishna K, Ahmed R, Perelson A: Recruitment times, proliferation, and apoptosis rates during the CD8(+) T-cell response to lymphocytic choriomeningitis virus. J Virol. 2001, 75: 10663-10669. 10.1128/JVI.75.22.10663-10669.2001.PubMed CentralView ArticlePubMed
- Kwiatkowski D: Malarial toxins and the regulation of parasite density. Parasitol Today. 1995, 11: 206-212. 10.1016/0169-4758(95)80079-4.View ArticlePubMed
- Karunaweera N, Grau G, Gamage P, Carter R, Mendis K: Dynamics of fever and serum levels of tumor necrosis factor are closely associated during clinical paraoxysms in Plasmodium vivax malaria. Proc Natl Acad Sci USA. 1992, 89: 3200-3203. 10.1073/pnas.89.8.3200.PubMed CentralView ArticlePubMed
- Karunaweera N, Wijesekera S, Wanasekera D, Mendis K, Carter R: The paroxysms of Plasmodium vivax malaria. Trends Parasitol. 2003, 19: 188-193. 10.1016/S1471-4922(03)00036-9.View ArticlePubMed
- Hjelm F, Carlsson F, Getahun A, Heyman B: Antibody-mediated regulation of the immune response. Scand J Immunol. 2006, 64: 177-184. 10.1111/j.1365-3083.2006.01818.x.View ArticlePubMed
- Cavanagh D, Elhassan I, Roper C, Robinson V, Giha H, Holder A, Hviid L, Theander T, Arnot D, McBride J: A longitudinal study of type-specific antibody responses to Plasmodium falciparum merozoite surface protein-1 in an area of unstable malaria in Sudan. J Immunol. 1998, 161: 347-359.PubMed
- Wipasa J, Suphavilai C, Okell L, Cook J, Corran P, Thaikla K, Liewsaree W, Riley E, Hafalla J: Long-lived antibody and B Cell memory responses to the human malaria parasites, Plasmodium falciparum and Plasmodium vivax. PLoS Pathog. 2010, 6: e1000770-10.1371/journal.ppat.1000770.PubMed CentralView ArticlePubMed
- McKay M, Beckman R, Conover W: A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics. 1979, 21: 239-245.
- Cash J, Karp A: A variable order Runge-Kutta method for initial-value problems with rapidly varying right-hand sides. ACM Trans Math Soft. 1990, 16: 201-222. 10.1145/79505.79507.View Article
- Press W, Teukolsky S, Vettering W, Flannery B: Numerical Recipes in C. 1992, Cambridge: Cambridge University Press
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.