Document Detail

Foetal haemoglobin and the dynamics of paediatric malaria.
Jump to Full Text
MedLine Citation:
PMID:  23190739     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
BACKGROUND: Although 80% of malaria occurs in children under five years of age, infants under six months of age are known to have low rates of infection and disease. It is not clear why this youngest age group is protected; possible factors include maternal antibodies, unique nutrition (breast milk), and the presence of foetal haemoglobin (HbF). This work aims to gain insight into possible mechanisms of protection, and suggest pathways for focused empirical work, by modelling a range of possible effects of foetal haemoglobin and other red blood cell (RBC) developmental changes on parasite dynamics in infants.
METHODS: A set of ordinary differential equations was created to investigate the leading hypotheses about the possible protective mechanisms of HbF-containing red blood cells, in particular whether HbF suppresses parasite population growth because parasite multiplication in individual RBCs is lower, slower or absent. The model also incorporated the intrinsic changes in blood volume and haematocrit that occur with age, and the possibility of parasite affinities for HbF-containing RBCs or reticulocytes.
RESULTS: The model identified several sets of conditions in which the infant remained protected, or displayed a much slower growth of parasitaemia in the first few months of life, without any intervening immune response. The most protective of the hypothesized mechanisms would be the inhibition of schizont division in foetal RBCs so that fewer merozoites are produced. The model showed that a parasite preference for HbF-containing RBCs increases protective effects for the host, while a preference for reticulocytes has little effect.
CONCLUSIONS: The results from this simple model of haematological changes in infants and their effects on Plasmodium falciparum infection dynamics emphasize the likely importance of HbF and RBC number as an explanatory factor in paediatric malaria, and suggest a framework for organizing related empirical research.
Authors:
Erica M W Billig; Philip G McQueen; F Ellis McKenzie
Related Documents :
3976379 - Effects of electronic fetal monitoring on rates of early neonatal death, low apgar scor...
6991659 - Pediatric follow-up of a randomized controlled trial of intrapartum fetal monitoring te...
1756059 - Birthweight-specific infant mortality risk in cesarean section.
23753649 - Computer programming: quality and safety for neonatal parenteral nutrition orders.
23557509 - Caterpillar cereal as a potential complementary feeding product for infants and young c...
23998449 - Challenges in the practice of human milk nutrition in the neonatal intensive care unit.
7136659 - A planned prospective evaluation of the anti-infective property of varying quantities o...
7567339 - Kernicterus in otherwise healthy, breast-fed term newborns.
8176569 - Increased incidence of bronchopulmonary dysplasia after antenatal administration of ind...
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Intramural     Date:  2012-11-28
Journal Detail:
Title:  Malaria journal     Volume:  11     ISSN:  1475-2875     ISO Abbreviation:  Malar. J.     Publication Date:  2012  
Date Detail:
Created Date:  2013-01-08     Completed Date:  2013-06-06     Revised Date:  2013-07-11    
Medline Journal Info:
Nlm Unique ID:  101139802     Medline TA:  Malar J     Country:  England    
Other Details:
Languages:  eng     Pagination:  396     Citation Subset:  IM    
Affiliation:
National Institutes of Health, Fogarty International Center, Building 16, Room 303, Bethesda, MD 20892, USA. erica.billig@nih.gov
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Age Factors
Blood Volume
Child, Preschool
Erythrocytes / metabolism,  parasitology
Fetal Hemoglobin / metabolism*
Host-Parasite Interactions
Humans
Infant
Infant, Newborn
Malaria, Falciparum / blood*,  parasitology*
Models, Biological
Parasitemia / blood,  parasitology
Plasmodium falciparum / growth & development,  pathogenicity
Reticulocytes / metabolism,  parasitology
Chemical
Reg. No./Substance:
9034-63-3/Fetal Hemoglobin
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Malar J
Journal ID (iso-abbrev): Malar. J
ISSN: 1475-2875
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Billig et al.; licensee BioMed Central Ltd.
open-access:
Received Day: 28 Month: 6 Year: 2012
Accepted Day: 17 Month: 11 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 28 Month: 11 Year: 2012
Volume: 11First Page: 396 Last Page: 396
PubMed Id: 23190739
ID: 3538578
Publisher Id: 1475-2875-11-396
DOI: 10.1186/1475-2875-11-396

Foetal haemoglobin and the dynamics of paediatric malaria
Erica MW Billig1 Email: erica.billig@nih.gov
Philip G McQueen2 Email: mqueenp@mail.nih.gov
F Ellis McKenzie1 Email: ellis.mckenzie@nih.gov
1National Institutes of Health, Fogarty International Center, Building 16, Room 303, Bethesda, MD, 20892, USA
2National Institutes of Health, Center for Information Technology, 12 South Drive, Bethesda, MD, 20892, USA

Background

Malaria disproportionately affects children under five years of age [1]. Repeated infection builds immune responses that protect against severe clinical disease and reduce parasitaemia [2-4], and, presumably, older children in endemic areas have been exposed more often and to a wider range of antigens than younger children. Hence age is often considered a surrogate for cumulative exposure, and age differences in response are attributed accordingly. There may be more to the explanation, however, given that intrinsic developmental differences between children under five and older children, or adults, may affect responses to malaria infection [5,6]. Among the developmental age differences that may affect the course and outcome of a malaria infection are the age-specific abundance and properties of red blood cells (RBCs).

A human malaria infection begins when the parasite enters the bloodstream in the saliva of an infectious, biting mosquito. The parasite invades the liver, multiplies there for approximately 10 days, and then releases thousands of merozoite forms into the bloodstream. Each merozoite invades a RBC, multiplies, and – if Plasmodium falciparum or Plasmodium vivax, the most common species – in about 48 hours bursts the RBC to release eight to 32 new merozoites, which invade new RBCs and continue the cycle [7]. Virtually all of the pathology of human malaria is associated with this cycle in RBCs, which continues until the host mounts an immune response sufficient to control the infection, or the host dies.

Presumably, the basic biological process of parasite invasion and multiplication in RBCs is independent of the literal age of a host (i e, the parasite does not “know” how old its host is), and so, whatever the host’s age, cycles proceed roughly as sketched above, with regular numbers of merozoites released at each cycle, in cycles of regular length. The total blood volume of a host increases very significantly with age and size [8,9]. Therefore, all else equal, overall parasitaemia will increase more rapidly in hosts with smaller total numbers of RBCs – younger children – and the consequent anaemia will become more severe more quickly. That is, while the RBC count of 4–6 × 106 ml-1 of blood may be similar between infants and adults, the blood volume of an adult is 10–15 times that of an infant, so any given population of parasites in an adult host will represent a much lower total parasitaemia than in an infant host.

Another notable difference between infants and older children is that a much higher proportion of RBCs in infants contain foetal haemoglobin (HbF) rather than adult haemoglobin (HbA). HbF is produced only in humans, apes and Old World monkeys: it has a stronger affinity for oxygen than HbA, and so facilitates oxygen transfer across the placenta from mother to foetus during gestation [10]. Beginning in the final trimester of pregnancy, and continuing into the first months of life, haemoglobin production usually largely switches from HbF to HbA, and the fraction of RBCs that contains HbF declines accordingly. The fraction of RBCs that contains HbF at birth varies from about 50% to 100% among neonates [11,12]. Variation may be related to differences in gestation period and/or in the complex polygenic control of haemoglobin production [13-15]. In most adults only about 1% of RBCs contain HbF, but in those with hereditary persistence of foetal haemoglobin (HPFH), a benign condition in which mutations or deletions in the β- or γ-globin genes or regulatory regions alter normal haemoglobin switching, 10-100% of RBCs contain HbF.

Infants under six months of age have a low rate of infection and incidence of severe disease compared to older children [16], despite the relative immaturity of their immune systems. Several factors may contribute to this protection, including passively transferred maternal antibodies, unique nutrition (breast milk), and elevated HbF [17]. Current evidence suggests that HbF retards the expansion of a Plasmodium population because Plasmodium parasites do not survive as well in HbF-containing RBCs as in those containing HbA [18-23], though details and specific mechanisms remain unclear. There could be several reasons for the retardation: the HbF cell may be a “dead end,” meaning that the parasite can invade but cannot replicate, or cannot escape, and so dies within the RBC. Alternatively, the growth process may be suppressed or slower in HbF-containing RBCs, such that fewer merozoites are released per infected RBC, or they take longer to develop to the point of bursting. In addition, P. falciparum may have a higher affinity for HbF-containing RBCs than HbA-containing RBCs. P. falciparum may also have (like P. vivax) a higher affinity for reticulocytes – the youngest age class of RBCs – which may be present in higher proportions in growing infants, and so may affect infection dynamics. In any case, the proportion of HbF-containing RBCs may be a factor in the protection of infants from severe clinical disease and high parasitaemia: by slowing parasite population growth, they may provide extra time for an effective response by the infant’s developing immune system or maternal antibodies.

Although several possible factors underlying variable protection in infants have been suggested, including the effects of maternal antibodies, nutritional differences, and HbF, none has yet been confirmed. Foetal haemoglobin is one of a number of mechanisms, reviewed in an earlier paper [6], that could help to explain the observed neonatal protection from severe disease. Infants have immature immune systems, in which key elements are developing asynchronously [6], so any added protective factors are likely to be of particular importance. This is a difficult set of topics to study, however, given the age group and multiple confounding factors.

No mathematical model to date has investigated the effects of any developmental age difference on a malaria infection. Here previous models of RBC dynamics in malaria infections are extended to focus specifically on the effects of lower RBC counts and higher HbF proportions on the dynamics of P. falciparum infections in infants [24-26].


Methods

A set of ordinary differential equations (ODEs) was used to model the dynamics of circulating RBCs, including the effect of infant growth on total blood volume (via the rate of growth of erythropoietic tissue), the switch in production of HbF- to HbA-containing RBCs, and the change in haematocrit in the first few months of life. The dynamics of two RBC lines were modelled – the uninfected (1) foetal and (2) adult erythrocytes – then the dynamics of malaria infection were incorporated into the dynamics of the circulating RBCs. Three main populations of Plasmodium parasites are involved in the pathology of malaria in an infant: (1) those in infected foetal erythrocytes, (2) those in infected adult erythrocytes, and (3) free merozoites in the blood. (Gametocytes and any cryptic sexual forms were ignored in the model.) Thus, the dynamical model considered has five separate populations, describing the dynamics of uninfected and infected foetal (HbF) RBCs, uninfected and infected adult (HbA) RBCs, and free merozoites within an infant host. For conciseness, vector and matrix notation (Additional file 1) are used to describe the dynamical equations.

Basic structure of the dynamical model

Consider an N × 1 vector P, the components of which may evolve in time. Define T(P) to be the sum of the components of P. Let the following ordinary differential equation (ODE) system in time t determine the evolution of the components:

[Formula ID: bmcM1]
(1) 
dP/dt=Stδ1−ΛDP

where S(t) is greater than or equal to zero for all t, Λ > 0 and dΛ = 0. The components of the vector δ(1) are zero except for the first one. Matrix D is sparse with all diagonal components = 1, and all components just below the diagonal = −1; see Additional file 1 for details. One can show that the contribution to –dT (P)/dt from S at time t - Δt is approximately a Gaussian function of Δt with a mean D = N Λ-1 and standard deviation σ = D N-1/2[24,27]. The ODE system described by equation (1) models the dynamics of individual organisms with total population T(P) and source term S(t). If individuals age with an average lifespan of D with standard deviation σ, then knowing the tangible quantities D and σ sets the abstract quantities N and Λ. (If D = σ, the ODE system reduces to one of simple exponential decay with a source term.) Although the vector–matrix notation is introduced for brevity, P itself can be thought of as containing all the information about the “state” of the population possible in this model. The matrix Λ D operates on the state vector P, incorporating the contribution of the source S(t) into the time evolution of P. This formalism may seem complicated initially, but it incorporates the non-instantaneous propagation of changes in the source S(t) into the population. One could think of other formalisms that incorporate time delays and dispersion in aging, but this particular one allows the use of efficient ODE solvers.

Total blood volume

In uninfected individuals, RBCs containing HbA circulate for ~120 days, while RBCs containing HbF circulate for ~70 days, at which point the RBC is cleared by the spleen [28]. In adults, there are ~5 × 106 RBCs per μl of blood. In neonates, the haematocrit dips within the first few months of life, and then rises again to adult levels [29]. Presumably, this transient dip is related to changes in RBC production in combination with the growing mass, and thus total blood volume, of the infant.

Since RBC and parasite populations are assessed empirically as cell counts per μl of blood, the model dynamics were formulated to consider the population per unit volume of blood rather than the total population in the host. However, the volume growth of the child must be carefully incorporated. If the state of a population of a parasite or RBC line is given by P, then consider ρ = VB-1XP, where VB is the blood volume of the host. (For simplicity, uniform mixing in the blood is assumed.) The time evolution of ρ is given by

[Formula ID: none1]
(none1) 
dρ/dt=dVB−1/dtP+VB−1dP/dt=s(t)δ1–(ΛD+ΛVI)ρ

where

[Formula ID: none2]
(none2) 
st=VB−1St

and

[Formula ID: bmcM2]
(2) 
ΛV=VB−1dVB/dt

The sum of the components of ρ, T(ρ), is VB-1X T(P), the density of the population in the blood. So from the point of view of population density, the growth of blood volume contributes a decay factor ΛV: a growing volume tends to dilute the concentration of cells. Note that ΛV itself is time dependent, as the rate of volume growth slows as the child ages.

Erythrocyte source

Near the time of birth, haemoglobin production usually switches from predominantly HbF to predominantly HbA. The timing of this switch, and the rate and degree of change, vary from infant to infant. The levels of circulating HbF usually reach adult levels of 1 – 2% around one year of age. The production of foetal erythrocytes per kilogram of mass of the host decreases nearly exponentially with age (t) as the host ages, while the production of adult erythrocytes per kilogram of mass of the host is a sigmoid-like function of t that reaches a maximum some months after birth. (For a review, see [13]) The time evolution of the erythrocyte sources in equation (5) is modelled by

[Formula ID: bmcM3]
(3) 
ESat=PamxXtanh((t+tOSa)/tCa)XWtVBt−1ESft=(PfmxXexp(−t/tCf)+Pfmn)XWtVBt−1

Here W(t) is the mass of the host at age t. The time evolution of W(t) and VB(t) is discussed in Additional file 1. The values of the parameters Pamx, tOSa, τCa, Pfmx, τCf and Pfmn are given in Table 1. Three pairs of time constants were used to model different changes in the rates of HbF and HbA production (Table 1). Using these parameters, the HbF:HbA ratio reaches 1:1 at 65, 97, and 163 days of age.

As a measure of the health of the host, the haematocrit hc was tracked, calculated as

[Formula ID: bmcM4]
(4) 
hc=100%XVaTRea+TMa+VfTRef+TMf

Here Vf and Va are the mean volumes of foetal and adult erythrocytes. Their values are given in Table 1. Rea and Ma are the population state vectors for the reticulocytes and mature RBCs for adult RBCs and Ref and Mf are the corresponding state vectors for the foetal RBCs; see next subsection. If the value of hc drops under 24%, it is assumed that the host dies.

Red blood cell populations

Some evidence indicates that P. falciparum preferentially invades young RBCs – reticulocytes – regardless of haemoglobin content [19,30], and also that P. falciparum preferentially invades HbF-containing RBCs [21]. It remains uncertain whether P. falciparum in fact prefers HbF-containing RBCs, or, given the shorter lifespan and changing production patterns in the first few months of life, the key factor is that HbF-containing RBCs have a younger age distribution, skewed toward reticulocytes. The ability to test these possibilities was incorporated into the model by creating a reticulocyte subgroup and a mature subgroup for both the foetal and the adult RBC populations. The model incorporates the effects of the switch from the production of foetal erythrocytes to adult erythrocytes, as well as the increase in blood volume as the host ages. The equations that describe the dynamics of the RBC populations are:

[Formula ID: bmcM5]
(5) 
dRea/dt=ESa(t)δ(1)−(ΛReaD+(ςReaμ+ΛV)I)ReadMa/dt=ΛReaL(Rea)δ(1)−(ΛMaD+(ςMaμ+ΛV)I)MadRef/dt=ESf(t)δ(1)−(ΛRefD+(ςRefμ+ΛV)I)RefdMf/dt=ΛRefL(Ref)δ(1)−(ΛMfD+(ςMfμ+ΛV)I)Mf

Here, I is the identity matrix and ESa(t) and ESf(t) are the source rates of production of adult and foetal reticulocytes, respectively, per unit volume of blood. These quantities change as the host ages, in a manner specified below. ζRea, ζMa, ζRef, and ζMf are the binding affinities of the merozoites to the respective RBCs; the values are defined in the next section. The state vector of the merozoites is μ. The rates ΛRea, ΛMa, ΛRef, ΛMf and the vector lengths were defined as in the subsection “Basic Structure of Dynamical Model” using the values specified in Table 2.

Parasite populations

Five morphologically distinct populations of asexual parasite cells were considered: (1) ring stage, (2) early trophozoites, (3) late trophozoites, (4) schizonts, and (5) merozoites. The possibility is allowed that the parasite can develop in either foetal or adult erythrocytes, but possibly at different rates of development and with different efficiencies of reproduction in each type. The list of populations as well as the average durations of individual residence in those populations, along with the standard deviations in the durations, is given in Table 2.

The primary release of merozoites from the liver apparently involves the release of thousands of merozoites in clusters [31]. For simplicity it is assumed that primary release occurs over a duration τPR = 24 hr with total of 105 merozoites released into the blood volume VB at a constant rate per unit volume κPR = 105 (24 h X VB)-1. The average number of merozoites subsequently released by bursting schizonts of P. falciparum in human blood culture is known from direct microscopy to be >16 [32], so in this report the average number of merozoites released from a bursting schizont formed in an adult erythrocyte is taken to be pa = 16. The possibility is allowed that the average number of merozoites released from a bursting schizont formed in a foetal erythrocyte, pf, may be different from pa. The binding affinity of merozoites to target RBCs is inferred from observations of parasite growth in neurosyphilis patients treated for malaria therapy [33,34], and inoculated volunteers [35], to be 10-6 – 10-5 μl hr-1. Here, ζMa, the binding affinity of a merozoite to the mature adult RBC, is taken to be 3.33333 × 10-6 μl hr-1. The possibility is allowed that the merozoite binding affinity to adult reticulocytes, ζRea, mature foetal RBCs, ζMf, and foetal reticulocytes, ζRef, may differ from ζMa. The value is only changed when directly stated.

Using the information above, letting t be the time since the birth of the infant, and taking AgPR as the age of the infant when primary release starts, the ODE system for the parasite populations is:

[Formula ID: bmcM6]
(6) 
dμ/dt=κPRΘt−AgPRΘ(AgPR−τPR−t)+paΛSaSaNSa+pfΛSfSfNSf−(ςReaT(Rea)+ςMaT(Ma)+ςRefT(Ref)+ςMfT(Mf)+Λμ+ΛV)μdRa/dt=(ςReaT(Rea)+ςMaT(Ma))μδ(1)−(ΛRaD+ΛVI)RadEa/dt=ΛRaL(Ra)δ(1)−(ΛEaD+ΛVI)EadLa/dt=ΛEaL(Ea)δ(1)−(ΛLaD+ΛVI)LadSa/dt=ΛLaL(La)δ(1)−(ΛSaD+ΛVI)SadRf/dt=(ςRefT(Ref)+ςMfT(Mf))μδ(1)−(LRfD+ΛVI)RfdEf/dt=ΛRfL(Rf)δ(1)−(ΛEfD+ΛVI)EfdLf/dt=ΛEfL(Ef)δ(1)−(ΛLfD+ΛVI)LfdSf/dt=ΛLfL(Lf)δ(1)−(ΛSfD+ΛVI)Sf

(Here, Θ (x) = 1 if x > 0, zero otherwise.)

Note that since Dμ = σμ, the state vector for the merozoites, μ, has only one component, so it was treated it as a scalar in equation (6). The rates Λμ, ΛRa, ΛEa, ΛLa, ΛSa, ΛRf, ΛEf, ΛLf, ΛSf and the vector lengths were defined as in the subsection “Basic Structure of Dynamical Model,” using the values specified in Table 2.

For further details on the Methods, see Additional file 1.


Results

Because little is known about the protective effects of HbF at the single-RBC level, several hypotheses about the failure of P. falciparum to develop normally within an HbF-containing RBC were modelled. First, the possibility was investigated that HbF-containing RBCs release 0, 2, or 16 merozoites (pf = 0, 2, or 16), while HbA-containing RBCs always release 16 (pa = 16). Second, the rate at which infants produce HbF and HbA during the first few months of life was varied (τCf =50, 100, or 200; τCa =75, 150 or 300). Third, the possibility was examined that P. falciparum has no preferential affinity, an increased affinity for HbF-containing RBCs, or an increased affinity for reticulocytes (ζRef > ζMf and ζRea > ζMa, or ζRef > ζRea and ζMf > ζMa). Last, the possibility that in HbF-containing RBCs the schizont stage is delayed 60 h (DSf =72), (so that the life cycle is 108 h instead of 48 h), was modelled. Figure 1 shows the changing rates of parasite population growth under these sets of conditions if an infant is infected at 40 days of age, for example. In these plots, it is clear that if HbF-containing RBCs release fewer merozoites per infected cell, the rate increases more slowly than if they release the same number as an HbA-containing cell. Adding other effects either amplifies (when pf = 0, the infection increases more slowly – graphically, a lower slope) or diminishes this effect (i e, the slope of pf = 0 and pf = 2 approaches that of pf = 16). Adding a delay is mildly protective for the host, as it slightly lowers the parasite population growth rate, which increases the time until host death. This change is seen in Figure 1B as a slightly lower slope in the number of infected RBCs. An increased affinity for HbF is very protective for the host: if an infant is infected at 40 days of age and pf = 0, for instance, as in Figure 1C, the infant actually clears the infection, rapidly. An increased affinity for reticulocytes is harmful for the host: parasite population growth increases more rapidly than if all RBCs are infected at equal rates (Figure 1D compared to Figure 1A).

With the parameters set to the values in Tables 1 and 2, the model shows that if HbF-containing RBCs release no new merozoites (pf = 0), the infant remains protected from a fatal outcome of infection through approximately the first 15 days of life, despite the complete absence of an immune response (Figure 2; Additional file 2). If the HbF-containing RBCs release only two merozoites (pf =2), the infant dies but the time until death increases, most significantly in the first two months of life. If the HbF-containing RBCs release 16 merozoites (pf = 16), the same as HbA-containing RBCs, there is no effect; i e, this serves as a control. Now, other parameter modifications will be introduced and compared to this set of parameters.

The rate of change of the ratio of HbF- to HbA-containing RBCs is known to vary dramatically from infant to infant. While some infants reach adult levels of HbA at three months of age, others do not reach adult levels until six months. Three different rates of change were modeled by modifying the values of τCf and τCa, such that the host reaches a HbF:HbA ratio of 1:1 at 65, 97, and 163 days of age. Figure 3 shows the rate of change of HbF- to HbA-containing RBCs when τCf = 100, and τCa = 150. Figure 4 shows the effect of these rate variations on the infection dynamics across the first 100 days of life. If HbF-containing RBCs release fewer merozoites than HbA-containing RBCs, then the infants with more HbF are more protected. Figure 4 shows that when τCf = 200 and τCa = 300 the infant is most protected from a fatal outcome.

In normal conditions, the P. falciparum life cycle within a RBC is 48 hours. It is possible that the presence of HbF delays this life cycle. A delay of 60 hours in the schizont phase of HbF-containing RBCs, leading to a total 108-hour life cycle (DSf =72) was modelled. This delay has little effect if the HbF-containing RBCs act as a “dead end” (0 merozoites released; pf = 0) for the parasite (Figure 5; also Figure 1B) compared to the normal life cycle. If pf =2, and the infant is infected at one day of age, the delay increases the time to host death by about 300 hours. This gain in time to host death slowly decreases as the age of the host at the time of infection increases. Similarly, if pf = 16 (the control) in an infant infected at one day of age, the delay increases the time to host death by 400 hours, a gain which slowly decreases as the age at infection increases.

Foetal haemoglobin-containing RBCs may be preferentially invaded due only to their younger age distribution, or there may be some inherent HbF-related RBC property for which P. falciparum has a greater affinity. This possibility was modelled by increasing the parasite’s affinity for HbF RBCs 10-fold (ζRef > ζRea and ζMf > ζMa). Compared to ζRef = ζRea and ζMf = ζMa, the stronger HbF affinity changes infection dynamics in favour of the host, most significantly when pf = 0 (Figure 6). Adding the 60-hour delay to this preferential binding affinity shifts infection dynamics even further in favour of the host, and makes the most significant change in time to host death from the two affinities without delay when pf = 2.

It has been reported that P. falciparum has a higher affinity for reticulocytes than for adult RBCs [30]. This possibility was incorporated into the model by increasing the parasite’s affinity for reticulocytes by 10-fold (ζRea > ζMa and ζRef > ζMf). An increased preference for reticulocytes slightly changes infection dynamics, although the maximum time to host death does not decrease dramatically (Figure 7). Incorporating a 60-hour delay in merozoite release in this scenario does not appreciably change the effects from those without the delay if pf = 0 or 2. If pf = 16 (the control), in an infant infected at one day of age, the time to host death increases by 400 hours, a gain which slowly decreases as the age at infection increases.


Discussion

The model presented here incorporates haematological changes that occur in infants under six months of age, and the results may offer some help toward explaining why neonates are relatively protected from severe malaria. The results are in accord with empirical studies which show that the youngest among the age group in question are protected, but the point at which the shift occurs varies based on transmission: in areas of high transmission, infants are protected only through the first three months of life, while in areas of low transmission, infants seem to be protected through six months [36]. It is possible that in areas of high transmission, HbF-containing RBCs are destroyed more quickly as a result of infection and, therefore, the switch to HbA-containing RBCs occurs more quickly. Thus, in areas of high transmission, infants may be protected for a shorter period of time. In a study of three regions of Togo, parasitaemia prevalence increased from 18.2% in infants zero to two months, to 43.0% in infants three to five months. The prevalence of anaemia was also reported to have increased from zero to two months to three to five months of age, with “anaemia more likely in children with parasitaemia.” Overall, there was no significant difference in parasite burden in the different regions, and, in all regions, most clinical disease was observed in infants under two years of age [37]. Another study examined infants younger than three months of age in Gabon, and found a parasite prevalence of 0.1%, concluding that parasitaemia in infants younger than three months of age is very uncommon [38]. A study in Malawi that measured the parasitaemia of all paediatric admissions (<15 yrs) found that 4.8% of positive smears were in children <6 months and 80.7% were in children between six months and five years of age. They concluded that although the risk of infection is lower in children in the younger age group than in older children, the risk does increase in early infancy and that this age group is at a greater risk for infection than previously thought [16].

In the model, when the parameters are set as in Tables 1 and 2, if pf = 0, the host is protected for the first two months of life, despite the complete absence of an immune response. For the next few months of life, the rate at which parasitaemia increases in the course of an infection is slower than in an adult host, and this difference may be sufficient to allow for effective immune intervention: any host factor that slows down parasite population growth buys time for protective responses by maternal antibody or the developing, immature infant immune system. Whatever the mechanism by which HbF retards the expansion of a Plasmodium population – whether parasite multiplication in an individual RBC is lower, slower or absent – the presence of HbF would appear to be an unmitigated benefit to the malaria-exposed, of any age. Yet the sparse evidence suggests that HPFH is rare, in areas of high malaria endemicity as elsewhere. Sickle cell anaemia and β-thalassaemia are generally accompanied by elevated HbF levels, but the degree of elevation is highly variable [39,40]. Because elevated HbF levels ameliorate the clinical symptoms of sickle cell disease and β-thalassaemia, and the switch from HbF to HbA production is neither complete nor irreversible, therapeutic reactivation is an active area of research [41]: hydroxyurea is the standard agent in current clinical practice [42]. The potential of these approaches for interventions in malaria-endemic regions apparently remains unexplored [43].

It seems possible that there are opportunities for laboratory researchers to gain important insights from in vitro work. For instance, if sufficient quantities of HbF-containing RBCs could be obtained, and parasite cultures established in them, critical observations could be compared to those from the usual cultures of HbA-containing RBCs, on the number of merozoites released per infected HbF-containing RBC, and the time from when the parasite enters the HbF-containing RBC until it bursts.

The hypothesis that HbF-containing RBCs release fewer merozoites has a greater overall effect on infection dynamics than a delay in the parasite life cycle. Empirical research suggests that the parasite may have a reticulocyte preference, and that this may confound the seeming preference for HbF in vivo through the differing age distributions of HbA- and HbF-containing RBCs [19]. Although a reticulocyte preference seems more beneficial to parasite persistence than host survival, in combination with a reduced merozoite release the overall effect still benefits the host.


Conclusion

It has long been accepted that children under six months of age are protected from severe infection, and this is widely reflected in the literature on paediatric malaria. This topic has recently been analysed more carefully, as noted above, and it appears that the malaria burden in this age group – while still relatively low – may be higher than previously acknowledged [16]. It may be that age-related variation within the overall age group arises from variation in HbF levels and transmission intensity.

The results from this simple model of haematological changes in infants and their effects on P. falciparum infection dynamics seem intuitively sound in a qualitative sense, and they emphasize the likely importance of HbF as an explanatory factor in paediatric malaria. The model provides a framework for examining hypotheses about the protective effects of foetal haemoglobin, organizing empirical observations and making critical quantitative comparisons among them.


Abbreviations

RBC: Red blood cell; P. falciparum: Plasmodium falciparum; HbF: Foetal haemoglobin; HbA: Adult haemoglobin; HPFH: Hereditary persistence of foetal haemoglobin; ODE: Ordinary differential equation.


Competing interests

The authors have declared no competing interests.


Authors' contributions

EB conducted the model runs and compiled the output. All authors contributed to the development of the model and wrote, read and approved the final manuscript.


Supplementary Material Additional file 1

Further Methods; Description: Further details on the methods.


Click here for additional data file (1475-2875-11-396-S1.doc)

Additional file 2

Explanation of Figure 2; Description: Further explanation of the interpretation of Figure2.


Click here for additional data file (1475-2875-11-396-S2.docx)


Acknowledgements

Fogarty International Center, NIH; Dr Gail Weinmann

The research contribution of PGM was supported by the Intramural Research Program of the NIH, Center for Information Technology.


References
WHOWorld Malaria ReportYear: 2010World Health Organization, Geneva
Carneiro I,Roca-Feltrer A,Griffin JT,Smith L,Tanner M,Schellenberg JA,Greenwood B,Schellenberg D,Age-patterns of malaria vary with severity, transmission intensity and seasonality in sub-Saharan Africa: a systematic review and pooled analysisPLoS OneYear: 20105e898810.1371/journal.pone.000898820126547
Reyburn H,Mbatia R,Drakeley C,Bruce J,Carneiro I,Olomi R,Cox J,Nkya WM,Lemnge M,Greenwood BM,Riley EM,Association of transmission intensity and age with clinical manifestations and case fatality of severe Plasmodium falciparum malariaJAMAYear: 20052931461147010.1001/jama.293.12.146115784869
Okiro EA,Al-Taiar A,Reyburn H,Idro R,Berkley JA,Snow RW,Age patterns of severe paediatric malaria and their relationship to Plasmodium falciparum transmission intensityMalar JYear: 20098410.1186/1475-2875-8-419128453
Baird JK,Age-dependent characteristics of protection v. susceptibility to Plasmodium falciparumAnn Trop Med ParasitolYear: 19989236739010.1080/000349898593669683890
Billig E,O'Meara W,Riley E,McKenzie FE,Developmental allometry and paediatric malariaMalar JYear: 2012116410.1186/1475-2875-11-6422394452
Miller LH,Good MF,Milon G,Malaria pathogenesisScienceYear: 18781994264
Raes A,Van Aken S,Craen M,Donckerwolcke R,Vande Walle J,A reference frame for blood volume in children and adolescentsBMC PediatrYear: 20066310.1186/1471-2431-6-316503982
Boer P,Estimated lean body mass as an index for normalization of body fluid volumes in humansAm J PhysiolYear: 1984247F6326366496691
Sankaran VG,Xu J,Orkin SH,Advances in the understanding of haemoglobin switchingBritish J HaematolYear: 201014918119410.1111/j.1365-2141.2010.08105.x
Cook CD,Brodie HR,Allen DW,Measurement of fetal hemoglobin in newborn infantsPediatricsYear: 19572027227813452668
Stamatoyannopoulos G,Molecular and cellular basis of hemoglobin switchingDisorders of hemoglobin: Genetics, pathophysiology, and clinical managementYear: 2001
Karlsson S,Nienhuis AW,Developmental regulation of human globin genesAnn Rev BiochemYear: 1985541071110810.1146/annurev.bi.54.070185.0052312411209
Forget BG,Progress in understanding the hemoglobin switchN Engl J MedYear: 201136585285410.1056/NEJMe110696921879905
Oneal PA,Gantt NM,Schwartz JD,Bhanu NV,Lee YT,Moroney JW,Reed CH,Schechter AN,Luban NL,Miller JL,Fetal hemoglobin silencing in humansBloodYear: 20061082081208610.1182/blood-2006-04-01585916735596
Larru B,Molyneux E,ter Kuile F,Taylor T,Molyneux M,Terlouw D,Malaria in infants below six months of age: retrospective surveillance of hospital admission records in Blantyre, MalawiMalar JYear: 2009831010.1186/1475-2875-8-31020038299
Riley E,Wagner G,Ofori M,Wheeler J,Akanmori B,Tetteh K,McGuinness D,Bennett S,Nkrumah F,Anders R,Lack of association between maternal antibody and protection of African infants from malaria infectionInfect ImmunYear: 2000685856586310.1128/IAI.68.10.5856-5863.200010992495
Nagel RL,Innate resistance to malaria: the intraerythrocytic cycleBlood CellsYear: 199016321339 discussion 340–329. 2257317
Pasvol G,Weatherall DJ,Wilson RJM,Effects of foetal haemoglobin on susceptibility of red cells to Plasmodium falciparumNatureYear: 197727017117310.1038/270171a0337159
Pasvol G,Wilson RJ,The interaction of malaria parasites with red blood cellsBr Med BullYear: 1982381331407052193
Pasvol G,Weatherall DJ,Wilson RJ,Smith DH,Gilles HM,Fetal haemoglobin and malariaLancetYear: 197611269127273695
Wilson RJ,Pasvol G,Weatherall DJ,Invasion and growth of Plasmodium falciparum in different types of human erythrocyteBull World Health OrganYear: 197755179186338178
Amaratunga C,Lopera-Mesa TM,Brittain NJ,Cholera R,Arie T,Fujioka H,Keefer JR,Fairhurst RM,A role for fetal hemoglobin and maternal immune IgG in infant resistance to Plasmodium falciparum malariaPLoS OneYear: 20116e1479810.1371/journal.pone.001479821532754
McQueen PG,McKenzie FE,Age-structured red blood cell susceptibility and the dynamics of malaria infectionsProc Natl Acad Sci U S AYear: 20041019161916610.1073/pnas.030825610115178766
McQueen PG,McKenzie FE,Competition for red blood cells can enhance Plasmodium vivax parasitemia in mixed-species malaria infectionsAm J Trop Med HygYear: 20067511212516837717
McQueen PG,McKenzie FE,Host control of malaria infections: constraints on immune and erythropoeitic response kineticsPLoS Comput BiolYear: 20084e100014910.1371/journal.pcbi.100014918725923
Lloyd AL,Destabilization of epidemic models with the inclusion of realistic distributions of infectious periodsProc Biol SciYear: 200126898599310.1098/rspb.2001.159911370974
Pearson HA,Life-span of the fetal red blood cellJ PediatrYear: 19677016617110.1016/S0022-3476(67)80410-45334979
Yamashita H,Kukita J,Ohga S,Nakayama H,Akazawa K,Ueda K,Serum erythropoietin levels in term and preterm infants during the first year of lifeJ Pediatr Hematol OncolYear: 19941621310.1097/00043426-199408000-00005
Pasvol G,Weatherall DJ,Wilson RJM,The increased susceptibility of young red cells to invasion by the malarial parasite Plasmodium falciparumBrit J HaematolYear: 19804528529510.1111/j.1365-2141.1980.tb07148.x7002199
Baer K,Klotz C,Kappe SHI,Schnieder T,Frevert U,Release of hepatic Plasmodium yoelii merozoites into the pulmonary microvasculaturePLoS PathogYear: 20073e17110.1371/journal.ppat.003017117997605
Glushakova S,Yin D,Li T,Zimmerberg J,Membrane transformation during malaria parasite release from human red blood cellsCurr BiolYear: 2005151645165010.1016/j.cub.2005.07.06716169486
McKenzie FE,Jeffery GM,Collins WE,Plasmodium vivax blood-stage dymanicsJ ParasitolYear: 20028852153512099421
Collins W,Jeffery G,A retrospective examination of the patterns of recrudescence in patients infected with Plasmodium falciparumAm J Trop Med HygYear: 199961444810432044
Hamilton Fairley BN,Sidelights on malaria in man obtained by subinoculation experimentsTrans R Soc Trop Med HygYear: 19474062167610.1016/0035-9203(47)90025-420243883
Snow RW,Nahlen B,Palmer A,Donnelly CA,Gupta S,Marsh K,Risk of severe malaria among African infants: direct evidence of clinical protection during early infancyJ Infect DisYear: 199817781982210.1086/5178189498474
Eliades MJ,Wolkon A,Morgah K,Crawford SB,Dorkenoo A,Sodahlon Y,Hawley WA,Hightower AW,Kuile FOT,Terlouw DJ,Burden of malaria at community level in children less than 5 years of age in TogoAm J Trop Med HygYear: 20067562262917038683
Klein Klouwenberg PMC,Oyakhirome S,Schwarz NG,Gläser B,Issifou S,Kiessling G,Klöpfer A,Kremsner PG,Längin M,Lassmann B,et al. Malaria and asymptomatic parasitaemia in Gabonese infants under the age of 3 monthsActa TropYear: 200595818510.1016/j.actatropica.2005.05.00315950165
Weatherall DJ,Phenotype-genotype relationships in monogenic disease: lessons from the thalassaemiasNat Rev GenetYear: 2001224525510.1038/3506604811283697
Thein SL,Menzel S,Lathrop M,Garner C,Control of fetal hemoglobin: new insights emerging from genomics and clinical implicationsHum Mol GenYear: 200918R216R22310.1093/hmg/ddp40119808799
Wilber A,Nienhuis AW,Persons DA,Transcriptional regulation of fetal to adult hemoglobin switching: new therapeutic opportunitiesBloodYear: 20111173945395310.1182/blood-2010-11-31689321321359
McGann PT,Ware RE,Hydroxyurea for sickle cell anemia: what have we learned and what questions still remain?Curr Opin HematolYear: 20111815816510.1097/MOH.0b013e32834521dd21372708
Makani J,Williams TN,Marsh K,Sickle cell disease in Africa: burden and research prioritiesAnn Trop Med ParasitolYear: 200710131410.1179/136485907X15463817244405
Miller MW,Brayman AA,Sherman TA,Abramowicz JS,Cox C,Comparative sensitivity of human fetal and adult erythrocytes to hemolysis by pulsed 1 MHz ultrasoundUltrasound Med BiolYear: 20012741942510.1016/S0301-5629(00)00350-111369128
White JBN,AS Fauci EB, Kasper DLMalariaHarrison's Principles of Internal MedicineYear: 2008McGraw-Hill, New York
Johnson JG,Epstein N,Shiroishi T,Miller LH,Factors affecting the ability of isolated Plasmodium knowlesi merozoites to attach to and invade erythrocytesParasitologyYear: 19808053955010.1017/S00311820000009986771738
Rapaport S,Introduction to HematologyYear: 1987J.B. Lippincott, Philadelphia

Figures

[Figure ID: F1]
Figure 1 

The number of uninfected and infected HbA- and HbF-containing RBCs (y-axis) in an infant infected at 40 days of age (x-axis: 1000 h). Each graph shows the rate of parasite population growth with pf = 0, 2, or 16. The number of uninfected RBCs is shown for pf = 2: there is no evident difference in these numbers with pf = 0 or pf = 16. A: Parameters set to values in Tables 1 and 2; B: 108 h life cycle (DSf =72); C: Parasite has a stronger affinity for HbF-containing RBCs (ζRef = 10ζRea and ζMf = 10ζMa); D: Parasite has a stronger affinity for reticulocytes (ζRef = 10ζMf and ζRea = 10ζMa).



[Figure ID: F2]
Figure 2 

The time until the host clears infection (y-axis: below h = 0) or the host dies (above h = 0; haematocrit < 24%), when infected, occurs between one and 100 days of age (x-axis). Parameters set to values in Tables 1 and 2. Results are shown for three different values of pf, the number of merozoites released per bursting schizont in foetal RBCs. The number of merozoites released from bursting schizonts in adult RBCs, pa, is 16 for all cases. For more explanation of this figure, see Additional file 2 .



[Figure ID: F3]
Figure 3 

The number of circulating HbF- and HbA-containing RBCs in the model, in an uninfected infant, with parameters set to the values in Tables 1 and 2 .



[Figure ID: F4]
Figure 4 

The effect of different rates of HbF:HbA production. Three distinct sets of parameters were used: τCf =50 and τCa =75; τCf =100 and τCa =150; or τCf =200 and τCa =300. As expected, when τCf and τCa are larger (i.e. there is more HbF for longer), the host is protected for a longer period of time. The solid lines are repeated from Figure 2 for comparison.



[Figure ID: F5]
Figure 5 

The time until the host clears infection (y-axis: below h = 0) or the host dies (above h = 0; haematocrit < 24%) when infected between one and 100 days of age (x-axis). Note that this figure contains the same plot as Figure 2, and the effect of pf = 0, 2, or 16 when there is a delayed, 108 h life cycle (ζRef = 10ζRea and ζMf = 10ζMa).



[Figure ID: F6]
Figure 6 

The time until the host clears infection or the host dies when infected between one and 100 days of age. In this figure, the parasite has a stronger affinity for HbF-containing RBCs, and either has a normal 48 h life cycle, or a delayed 108 h life cycle (ζRef = 10ζRea and ζMf = 10ζMa).



[Figure ID: F7]
Figure 7 

The time until the host clears infection or the host dies when infected between one and 100 days of age. In this figure, the parasite has a stronger affinity for reticulocytes, and either has a normal 48 h life cycle, or a delayed 108 h life cycle (ζRef = 10ζMf and ζRea = 10ζMa).



Tables
[TableWrap ID: T1] Table 1 

Hematological parameters


Parameter Symbol Value
Maximum adult reticulocyte production rate per kilogram of tissue
Pamx,
1.21 × 108 (kg hr) -1
Maximum foetal reticulocyte production rate per kilogram of tissue
Pfmx
1.7 × 108 (kg hr) -1
Minimum foetal reticulocyte production rate per kilogram of tissue
Pfmn
2.48 × 106 (kg hr) -1
Time offset for production of adult erythrocytes
tOSa
50 days
Time constant for production of adult erythrocytes
τCa
75, 150, and 300 days
Time constant for decay of production of foetal reticulocytes
τCf
50, 100, and 200 days
Volume of adult red blood cell
Va
8.0 × 10-8 μm3
Volume of foetal red blood cell Vf 1.25 × 10-7 μm3

Baseline parameters for adult (HbA) and foetal (HbF) RBCs [44].


[TableWrap ID: T2] Table 2 

Populations used in the model


Population State vector Average duration (h) Standard deviation (h)
Merozoites
μ
Dμ = 0.1
σ μ = 0.1 = Dμ
Ring stage in adult RBC
Ra
DRa = 12
σRa = 1.2
Early trophozoite in adult RBC
Ea
DEa = 12
σEa = 1.2
Late trophozoite in adult RBC
La
DLa = 12
σLa = 1.2
Schizont in adult RBC
Sa
DSa = 12
σSa = 1.2
Ring stage in foetal RBC
Rf
DRf = 12
σRf = 1.2
Early trophozoite in foetal RBC
Ef
DEf = 12
σEf = 1.2
Late trophozoite in foetal RBC
Lf
DEf = 12
σLf = 1.2
Schizont in foetal RBC
Sf
DSf = 12 or 72
σSf = 0.1
Adult Reticulocyte
Rea
DRea = 36
σRea = 6
Adult mature RBC
Ma
DMa = 2844
σMa = 168
Foetal Reticulocyte
Ref
DRef = 36
σRef = 6
Foetal mature RBC Mf DMf = 1644 σMf = 120

Durations of intracellular parasite stages in adult (HbA) and foetal (HbF) RBCs [45]. Merozoite duration [46], Adult erythrocyte durations: [47], Foetal erythrocyte durations [28]. The standard deviations are plausible guesses.



Article Categories:
  • Research

Keywords: Malaria, Plasmodium falciparum, Foetal haemoglobin.

Previous Document:  High-Throughput Screen Identifies Cyclic Nucleotide Analogs That Inhibit Prostatic Acid Phosphatase.
Next Document:  Ultrasound detects occult entheseal involvement in early psoriatic arthritis independently of clinic...