# Quantifying the removal of red blood cells in *Macaca mulatta* during a *Plasmodium coatneyi* infection

- Luis L. Fonseca
^{1, 4}Email author, - Harnel S. Alezi
^{1}, - Alberto Moreno
^{2, 4}, - John W. Barnwell
^{3, 4}, - Mary R. Galinski
^{2, 4}and - Eberhard O. Voit
^{1, 4}

**Received: **22 March 2016

**Accepted: **29 July 2016

**Published: **12 August 2016

## Abstract

### Background

Malaria is the most deadly parasitic disease in humans globally, and the long-time coexistence with malaria has left indelible marks in the human genome that are the causes of a variety of genetic disorders. Although anaemia is a common clinical complication of malaria, the root causes and mechanisms involved in the pathogenesis of malarial anaemia are unclear and difficult to study in humans. Non-human primate (NHP) model systems enable the mechanistic study and quantification of underlying causative factors of malarial anaemia, and particularly the onset of severe anaemia.

### Methods

Data were obtained in the course of *Plasmodium coatneyi* infections of malaria-naïve and semi-immune rhesus macaques (*Macaca mulatta*), whose red blood cells (RBCs) were labelled in situ with biotin at the time the infections were initiated. The data were used for a survival analysis that permitted, for the first time, an accurate estimation of the lifespan of erythrocytes in macaques. The data furthermore formed the basis for the development and parameterization of a recursive dynamic model of erythrocyte turnover, which was used for the quantification of RBC production and removal in each macaque.

### Results

The computational analysis demonstrated that the lifespan of erythrocytes in macaques is 98 ± 21 days. The model also unambiguously showed that death due to senescence and parasitaemia is not sufficient to account for the extent of infection-induced anaemia. Specifically, the model permits, for the first time, the quantification of the different causes of RBC death, namely, normal senescence, age-independent random loss, parasitization, and bystander effects in uninfected cells. Such a dissection of the overall RBC removal process is hardly possible with experimental means alone. In the infected malaria-naïve macaques, death of erythrocytes by normal physiological senescence processes accounts for 20 % and parasitization for only 4 %, whereas bystander effects are associated with an astonishing 76 % of total RBC losses. Model-based comparisons of alternative mechanisms involved in the bystander effect revealed that most of the losses are likely due to a process of removing uninfected RBCs of all age classes and only minimally due to an increased rate of senescence of the uninfected RBCs.

### Conclusions

A new malaria blood-stage model was developed for the analysis of data characterizing *P. coatneyi* infections of *M. mulatta*. The model used a discrete and recursive framework with age-structure that allowed the quantification of the most significant pathophysiological processes of RBC removal. The computational results revealed that the malarial anaemia caused by this parasite is mostly due to a loss of uninfected RBCs by an age-independent process. The biological identity and complete mechanism of this process is not fully understood and requires further investigation.

### Keywords

*Macaca mulatta*Malarial anaemia Mathematical model

*Plasmodium coatneyi*Red blood cell removal

## Background

Malaria, caused by infection with parasites of the *genus Plasmodium,* is responsible for over half a million deaths per year worldwide, with children being the main victims [1]. The immense severity of the disease and the long coexistence of humans and *Plasmodium* parasites have even led to the emergence and perpetuation of possible protective genetic disorders like thalassaemias and sickle-cell disease, as well as haemoglobin C, haemoglobin E, and G6P dehydrogenase and pyruvate kinase deficiencies [2–7]. In primates, *Plasmodium* sporozoites infect, transform and multiply within parenchymal hepatocytes to form merozoites; these developmental life cycle processes are asymptomatic for the host. Once merozoite forms of the parasite are released from the infected hepatocytes, the parasite begins its cyclical blood-stage development. Merozoites invade and multiply within red blood cells (RBCs) every 24, 48 or 72 h, depending on the species, and new merozoite progeny are released to invade other RBCs [8]. During this cyclical process of invasion and destruction of RBCs, the symptoms and clinical complications associated with malaria emerge. The blood stage of an infection is characterized by repeated rounds of RBC invasion which, if not kept under control by host immune responses or anti-malarial treatment, can lead to exponential growth of the parasite, with a concomitant destruction of the parasitized RBCs. This destruction, however, is not the predominant mechanism of RBC removal that leads to anaemia, and indeed seems to be vastly surpassed by the destruction of uninfected RBCs (uRBCs) [9–11].

Infections with *Plasmodium coatneyi,* a simian malaria species that is closely related to *Plasmodium knowlesi* [12, 13], mirror the biology and pathogenesis of *falciparum* malaria, with severe forms of pathology including anaemia. To study mechanisms of the onset and recovery of anaemia, Moreno et al. [14] established procedures to measure the turnover of in vivo biotinylated RBCs in rhesus macaques (*Macaca mulatta*). These macaques had been experimentally infected for the first time (*i.e.,* when malaria naïve) with *P. coatneyi* infected RBCs and, then again, while partially immune, 9 months after curative anti-malarial drug treatment. Five infected and then re-challenged (semi-immune) animals were compared to five control rhesus with biotinylated RBCs, but no malaria infection. Microscopy-based counts of infected RBCs and haemoglobin levels were monitored daily, and the numbers of biotinylated RBCs were assessed using flow cytometry. This work demonstrated that malarial infections result in an accelerated turnover of uninfected RBCs, and this was most pronounced in the malaria naïve animals. The precise mechanisms causing malarial anaemia are unknown, but have been suggested to be due to multiple possible factors leading to a reduction in circulating RBCs including physicochemical membrane changes, reduced erythrocyte deformability, accelerated erythrocyte senescence, and immunological reactions that cause their removal [14, 15].

This paper presents a mathematical model that was developed to study RBC dynamics in circulating blood during malaria infections. The model was parameterized using experimental results from Moreno et al. [14] and implemented as a discrete recursive structure that was previously identified as best suited for this class of problems [16]. The choice of this framework was based on the need to account for the aging of the RBCs rather accurately, which is problematic in delayed differential and integro-differential equations that do not track age directly, but instead approximate the time passed since the cell was generated. Ordinary differential equations with age classes do not model the aging accurately since inspection of age-classes reveals a distribution of ages within each class. Partial differential equation models like the Lotka–McKendrick age-structured population model [17, 18] do accurately model the aging of a population, but are difficult to implement, especially if the model includes variables with and without an age-structure and if the aging process is perturbed by events like a malaria infection. The most effective alternative is a discrete recursive framework with age-classes. In a sense, this structure corresponds to the discretization of a PDE model and hence shares all of its properties, but is more easily implemented and faster to solve, with an accuracy that is readily tuned by stipulating a desired time-step.

The model permits the quantification of the production of newly generated RBCs and of the different processes leading to the removal of RBCs in the absence or presence of a *P. coatneyi* infection in malaria naïve or semi-immune rhesus macaques. Additionally, two alternative mechanisms of uninfected RBC removal, namely accelerated erythrocyte senescence and immunologic removal, are modelled and compared on the basis of their respective predictions. The results demonstrate that the destruction of uninfected RBCs was the dominant process underlying malarial anaemia in the *P. coatneyi* infections reported by Moreno et al. [14], and that the direct destruction of infected RBCs by the parasite accounted for only about 4 % of the total RBC loss. Beyond this specific result, the model can be employed as a tool for predicting and exploring disease severity and evaluating host-directed interventions. This capability includes the study of other species of *Plasmodium* that cause malaria in primates, each with their unique blood-stage biology and pathogenic characteristics [8].

### Red blood cell removal processes

In healthy humans and non-human primates (NHPs), RBCs are produced by the erythropoietic system in the bone marrow, which is under the control of several cytokines including, in particular, erythropoietin (for a review see [19]). Removal of RBCs is a task for the phagocytic arm of the immune system in response to injury, senescence, or other processes. Injured RBCs are usually removed in the spleen, as they fail to deform and can no longer pass through the microcirculation of the red pulp [20].

Senescence-driven removal of RBCs typically occurs due to oxidative stress to haemoglobin, which results from the continuous cycling between normoxia and hypoxia. This cycling between oxidative states eventually triggers the formation of methemoglobin which ultimately denatures into haemichromes [21]. These haemichromes are able to bind to the cytosolic side of AE1 protein (Band 3) and thus to displace ankyrin, which weakens the AE1 connection to the cytoskeleton and eventually results in clustering of AE1 [22, 23]. On the external side, AE1 protein, which is usually found in dimers, has a low affinity for binding naturally occurring IgG antibodies (NAbs). However, once bound to haemichromes on the inside, AE1 dimers are able to aggregate and bind to the NAbs with enhanced affinity [24, 25]. NAbs alone are not efficient at promoting RBC clearance, but they are able to activate the classical pathway of the complement system and thereby induce erythrophagocytosis [26, 27]. It has also been proposed that AE1 recognition by NAbs may depend on the proteolytic degradation of AE1 [28].

## Methods

### Determination of the hazard function of RBCs in *Macaca mulatta*

*c*, was estimated to be 4.36 × 10

^{−5}.

*x*is the age of the RBC in hours. In this case, the parameter values were determined by least-squares optimization (Matlab:

*fminsearch*) as

*a*= 8.488 × 10

^{−45}and

*b*= 12.25.

*Sc*(

*t*)) can then be estimated at time

*t*as the number of living cells older than

*t*:

*Prod*is the hourly production of RBCs. Its value was estimated by least-squares optimization (Matlab:

*fminsearch*) of

*Sc(t)*against experimental data that were obtained by tracking the surviving biotinylated RBCs in a cohort of five healthy uninfected macaques for over 120 days [14].

### Host-pathogen interaction model in *Macaca mulatta*

A host-pathogen interaction model was developed using a discrete recursive framework with age-classes, as this structure allows effective modelling of the exact amount of time a cell survives [16]. Uninfected RBCs exist in two states. Both populations, unlabelled (*RBC*) and biotinylated (*bRBC*) RBCs, are represented with 3840 age-classes (allowing a maximum lifespan of 160 days or 3840 h, which is intentionally chosen longer than a priori expected to ensure total removal). In the absence of any removal processes, this set-up would generate a rectangular age-distribution with the same number of cells in each class.

*H*) is proportional to the total amount of RBCs (Eq. 4) and dependent on the ratio of haemoglobin to RBC (R

_{H/RBC}), which is specific to each macaque:

The malarial merozoites (*M*) are able to infect RBCs, independently of their biotinylation status, and thus produce infected RBCs (*iRBC*). These, in turn, live for another 48 h, thus occupying 48 age classes, at the end of which they burst and release 20 new merozoites [8, 36]. Exposure of the immune system to infected RBCs leads to the up-regulation of the immune response against infected RBCs (*I*), which in turn triggers the removal of infected RBCs by the spleen.

At the beginning of each simulation, all RBCs are biotinylated, while subsequent erythropoiesis produces exclusively non-biotinylated RBCs. Independent of their biotinylation status, all RBCs are equally subject to the different removal processes, as discussed before. According to the literature [32–34], age-independent (random) death accounts for the normal removal of 10 % of all produced RBCs. This process was included in the hazard function as the constant term *c*. Senescent death was modelled with the power-law term of the hazard function shown in Eq. (1), and removal of uninfected RBCs was modelled as an age-independent process and quantified by least-squares optimization using the experimental data for each of the infected macaques.

### Model formulation

The model has five dependent variables (*RBC*, *bRBC*, *iRBC*, *M*, *I*), of which *RBC* and *bRBC* have 3840 age classes each and *iRBC* has 48 age classes. Additionally, an auxiliary dependent variable, *rRBC*, was used to model the pool of cells lost due to the age-independent (random) process. The variables *M*, *I* and *rRBC* do not have an age-structure.

*RBC*,

*bRBC*and

*rRBC*are modelled recursively as follows:

In this formulation, *r*
_{
S
} is the rate of senescence, which has a default value of 1 for the normal rate of senescence where phenotypical age equals chronological age, and a larger value when senescence occurs at a faster rate than the chronological aging of the RBCs. *Sf*
_{
i,t
} is the surviving fraction of RBCs in age-class *i* at time-point *t* that will reach the next time point (*t* + 1). The hazard function *hf*
_{
i
} was estimated in Eq. (1); it determines the relative number of cells lost from age-class *i*, due to age. *RD*
_{
t
} is the total number of cells lost due to random death at time *t* and approximated in Eq. (1) as the age-independent term *c*. *P*
_{
Ery,t
} is the rate of erythropoietic production of RBCs at time *t*; it is calculated in Eq. (19) as the product of the specific RBC production rate (*Prod*; which is a property of the healthy state of each macaque) and the daily relative change in erythropoietic production (*DP*
_{
f
}). *uR*
_{
t
} is the number of RBCs lost during an infection from the peripheral blood at time *t* by processes other than parasitization and is calculated by interpolation of the daily number of uninfected RBCs (*DuR*
_{
f
}) to be removed by these processes. *TRBC*
_{
t
} is the total number of healthy RBCs in circulation at time *t*, while *TiRBC*
_{
t
} is the total number of susceptible RBCs (younger than 50 days) at time *t*. *IA*
_{
t
} is a binary vector of “infectability” for a given age class of RBCs. Equations (5) and (12) allow RBCs in both pools (*bRBC* and *RBC*) to age at a rate faster than their chronological age would normally dictate. This speed-up is achieved by moving cells by *r*
_{
S
} age-classes if *r*
_{
S
} is integer; if it is not an integer, then the cells in one age class are proportionally moved to two age classes *r*
_{
SI
} and *r*
_{
SI
} + 1, where the ratio of cells moved to each of these age classes is calculated based on the decimal part of *r*
_{
S
}, *r*
_{
SD
}. The closer *r*
_{
SD
} is to 1, the larger proportion of cells that are moved by *r*
_{
SI
} + 1 age classes. Since RBCs of any age may die at any time point, *Sf*
_{
i,t
} (Eq. 9) allows the removal of cells by four different processes: random loss [age-independent, Eqs. (16, 17)]; age-dependent, according to the hazard function, Eq. (11); removal of uninfected RBCs [infection dependent, Eq. (18)] and by parasitization, (Eqs. 22–26). Random loss is modelled in a way similar to the proposal by Löffler’s group [33]. Namely, it is assumed that 10 % of all cells produced will eventually be lost by this random process. Unlike Löffler’s model, it is assumed here that the proportion of cells lost randomly is fixed at 10 % rather than being a function of the number of cells being produced. This change was made because a variable ratio caused the model to become unstable under some conditions. The 10 % of cells produced accumulate in an auxiliary variable (*rRBC*) from where they are removed (Eq. 16) with a first-order rate of (800 h)^{−1}. This rate was adjusted from the rate of (1020.4 h)^{−1} for humans [33], based on the fact that the lifespan of rhesus macaque RBCs is approximately one fifth shorter than in humans (100 and 120 days, respectively). The number of cells removed from this auxiliary variable (rRBC) is recorded with Eq. (17), which is then used to estimate the number of cells that are uniformly removed from all age classes of RBCs (Eq. 9).

*iRBC*,

*M*and

*I*pools are recursively modelled as:

*k*is the rate of iRBC clearance by the immune system;

*s*is the maximal rate of the action of the immune response; and \(\varphi\) is the number of iRBCs that produce half-maximal action of the immune response.

### Quantification of the increased removal or senescence rate of uninfected RBCs

The quantification of the production and removal of uninfected RBCs was performed through non-linear regression analysis using the experimental data from Moreno et al. [14] (Fig. 1).

*DuR*

_{ f }) is zero and the daily relative erythropoietic production (

*DP*

_{ f }) is one:

For each control macaque, the specific RBC production rate [*Prod*; Eq. (19)] and haemoglobin-to-RBC ratio (*R*
_{
H/RBC
}) were obtained by least-squares optimization of the trajectories for the total number of biotinylated RBCs \(\left( {TbRBC_{t} = \sum\nolimits_{i = 1}^{3840} {\left( {bRBC_{i,t} } \right)} } \right)\) and the haemoglobin levels \(\left( {H_{t} = R_{H/RBC} \cdot TRBC_{t} } \right)\) against the experimentally obtained time-series data for these variables.

- 1.The infected macaques were simulated first as control macaques. Then the time of infection (
*toi*) and the corresponding initial number of infected cells (*noi*) were determined such thatresults in a simulated trajectory for the iRBCs that approximates the initial growth of the infected RBCs. The control and elimination of the infection were obtained by optimization of the parameters$$M_{t = toi} = noi$$(30)*k*,*s*and \(\varphi\) (Eqs. 24, 27) for the semi-immune, infected macaques. For the naïve macaques, a sub-curative treatment was administered at peak parasitaemic periods post-infection as these macaques were otherwise not able to control their acute parasitaemias and would probably have died. The solutions for the parameters*toi*,*noi*,*k*,*s*and \(\varphi\) were not unique and the actual values were therefore not interpreted biologically. They were used exclusively as a means to simulate the macaque-specific parasitaemia levels. - 2.
Next, the macaque-specific RBC production rate (

*Prod*) and haemoglobin-to-RBC ratio (*R*_{ H/RBC }) were obtained by least-squares optimization against the first five time-points available for the biotinylated RBCs and haemoglobin levels of each macaque. - 3.The mechanism-dependent parameters were then obtained in a mutually exclusive way.
- 3.1.
Age-independent mechanism: The experimentally obtained time series data for the biotinylated RBCs were fitted by optimization of daily levels of uninfected RBC removal (\(DuR_{f} ,\quad f \in \left\{ {1,2, \ldots ,30} \right\}\)) while keeping

*p*_{1}= 0 (Eq. 20). - 3.2.
Increased senescence mechanism: The rate of senescence [

*r*_{ S }:*p*_{1}*,…,p*_{5}in Eq. (20)] was fitted against experimentally obtained time-series data for the biotinylated RBCs, while keeping \(DuR_{f} = 0\,\,( f \in \left\{ {1,2, \ldots ,35} \right\}\)).

- 3.1.
- 4.
The levels of erythropoietic up-regulation of RBC production (\(DP_{f} , \,\,f \in \left\{ {1,2, \ldots ,30} \right\}\)) were optimized by fitting the experimental time-series data for the haemoglobin levels.

- 5.
A re-fitting of

*DuR*_{ f }and*DP*_{ f }(age-independent model) or of*r*_{ S }and*DP*_{ f }(senescence model), using the previous results (2.1 and 3 or 2.2 and 3) as starting points, against the experimental biotinylated RBC numbers and haemoglobin levels was performed to obtain the best overall fit.

## Results

### Modelling the life span of RBCs

*RBC*), while the second pool contains biotinylated RBCs (

*bRBC*). It can be reasonably assumed that the spleen and immune system treat both types in the same manner.

In the beginning of the experiments that are modelled here, all RBCs are labelled (Fig. 1). Therefore, initially only the biotinylated pool contains cells, while production of new RBCs only occurs at the first age-class of the unlabelled pool (Erythropoiesis, Fig. 2). Both pools are assumed to be subject to the same hazard function, which is represented as Eq. (11) in the “Methods” section and with parameters estimated in the “Survival analysis of RBCs in *M. mulatta*
” section. This hazard function (Eq. 1) accounts for the two normal processes of RBC removal present in a healthy host, namely age-dependent death due to senescence, and age-independent death. The latter type of destruction will be referred to as *random*, which follows the terminology of the field [32, 33, 37], even though it is unclear whether the process is truly stochastic or whether it is governed by deterministic processes that are not well understood. This age-independent random death is modelled under the assumption that 10 % of all RBCs produced will eventually die by this process (Eqs. 16, 17). The initial age distribution of bRBCs is calculated by allowing the model to run to steady state, and the age distribution of the variable *RBC* (unlabelled RBCs) is recorded in this state.

To model the blood stage of a malarial infection, a pool of merozoites (*M*) is included with a single 1-h age class since, once released from infected RBCs (*iRBCs*), merozoites target and infect available uninfected RBCs within a matter of minutes. Also taken into account is a pool of infected RBC (*iRBC*) with 48 age classes, corresponding to the 2 days of the intra-erythrocytic lifespan of *P. coatneyi*. To mimic the body’s response to the infection, an immune response (*I*) is introduced as a black box. It acts such that exposure to parasite antigen represented by the *iRBC* leads to an increase in the immune response (*I*), which in turn kills iRBCs in a process that is linearly dependent on the number of iRBCs and on the immune response (Immune-removal, Fig. 2; Eq. 24). Even with this small number of processes representing the immune response, the model does not have a unique solution for each of the macaques in terms of immune-related parameters. For this reason, no other processes, such as immune exhaustion, were included and the numerical values of the parameters will not be discussed, as this immune module is thought to be a utilitarian module without a specific foundation for biological interpretability.

### Survival analysis of RBCs in *Macaca mulatta*

### Quantification of the enhanced removal of uninfected RBCs

The model distinguishes and allows the quantification of RBC loss due to four distinct processes (Eq. 9): age-dependent elimination due to senescence, age-independent random removal, parasitization, and removal of uninfected RBCs by other means. Age-independent removal of RBCs is assumed to be constant, because it is known that 10 % of all newly formed RBCs will be removed at some unpredictable time point before the natural life span is reached (see [16, 32, 33]). Senescent removal is modelled in accordance with the hazard function, which was determined from healthy rhesus macaques (Fig. 3b). Loss of RBCs due to parasitization refers to RBCs that are infected by merozoites and subsequently destroyed upon maturation of the new intracellular parasite progeny. The designation of *removal of uninfected RBCs by other means* excludes all of the other removal processes and is required in order to match the experimental data. For the non-infected control macaques, the model takes into account only the first two processes: age-independent removal and senescence of RBCs.

Control non-infected macaques lost RBCs only due to age-independent (random) and age-dependent (senescence) processes (Fig. 4); in these macaques, the production of RBCs matches exactly the natural loss of cells. The infected malaria-naïve macaques (i.e. not infected previously), by contrast, lost more cells than they produced. Most of their RBC loss is accounted for by the elimination of uninfected RBCs (76 % of total), whereas only a small portion was actually parasitized (4 % of total, or 5 % of the infection-induced loss). Relative to control macaques, infected naïve and semi-immune macaques lose fewer RBCs due to senescence (Fig. 4), simply because so many cells are being removed by other processes that considerably fewer RBCs reach old age. Comparatively however, more cells in both naïve and semi-immune infected macaques are lost through age-independent processes (Fig. 4). The reason is that more cells are being produced in these macaques in response to the infections and 10 % of these cells are removed by age-independent process.

The infected semi-immune macaques had milder parasitaemias (Fig. 1) and, therefore, lost fewer cells to parasitization (1 % of total, *p* < 0.001) than the malaria-naïve infected macaques (Fig. 4). Similarly, they also experienced less cell loss due to the removal of uninfected RBCs, although the difference is not significant (*p* < 0.11). The infected macaques, regardless of whether they were infected when naïve or semi-immune, also exhibited a higher erythropoietic output than the control macaques (*p* < 0.007 and *p* < 0.001, respectively) but the difference between the two infected groups is not significant (*p* < 0.42). The term “Non-infected Production” in Fig. 4 denotes the number of RBCs produced by the infected macaques if they had not been infected; thus, production above this level is required to compensate for the infection-induced destruction. The “non-infected production” of RBCs is not significantly different among the three groups evaluated here (malaria naïve, semi-immune, or control). Strikingly, the infection caused an increase in total RBC production to 290 and 235 % in the naïve and semi-immune infected macaques, respectively, relative to their corresponding production in the absence of infection (Fig. 4).

### Model validation

### The process of removal of uninfected RBCs

## Discussion

This work introduces a mathematical model that allows for the disentanglement of concurrent processes relating to RBC removal and malarial anaemia, based on published experimental data from a longitudinal malaria infection study with NHPs [14]. The main results present a clearer picture of the life span of RBCs in rhesus macaques and the quantification and age-characterization of the enormous proportion of uninfected RBCs that are eliminated coincident with the rise in parasitaemia and progression of the disease.

Developing appropriate models for such purposes requires a flexible computational framework capable of accounting for large numbers of cells in their proper age classes. In particular, a requirement here was that cells could normally be kept in various age classes for more or less fixed amounts of time, but that they could also skip age classes or remain longer within the same classes in certain situations. One approach of addressing such a situation could be based on delay differential equations (DDEs). However, it was shown elsewhere that DDEs are not flexible enough for modelling malarial anaemia at the cellular level [16]. Delays could also be generated with ordinary differential equations (ODEs) with age-classes, but these do not represent the transitions between age classes well. In contrast to these standard approaches, discrete recursive equations with age classes are not only able to generate hard delays but can also represent the age structure correctly, as can be seen by the fact that all cells in a given age class have the same age [16].

Building upon these advantages, a discrete system was constructed with age classes and a hazard function was added to account for death. This structure allowed an effective representation of the age distribution of RBCs. In the absence of cell death according to the hazard function, the cells would have a rectangular age distribution with constant death rate. Instead, use of a data-driven hazard function allowed cell removal to be modelled in a manner that is very close to reality. The hazard function that fitted the experimental results best and was used here is a power-law function, which captures the fact that cells of all ages can die, but that the number of cells dying at a young age is very small. By the same token, the older a cell becomes the more likely it will die, and the power-law hazard function captures this behaviour well (Fig. 3b).

Two mechanisms were tested to explain the removal of uninfected RBCs during a malarial infection by processes beyond the normal physiological age-dependent and age-independent processes and the direct parasitization of RBCs. These mechanisms underlying the loss of uninfected RBCs, sometimes referred to as the bystander effect, corresponded to two hypotheses: (1) that all cells are equally likely to be removed, independent of age; or (2) that the removal occurs by normal senescence, but that during high parasitaemia levels uninfected RBCs age faster. The age-independent mechanism would seem to be the better choice if bystander RBC removal occurs due to age-independent loss of uninfected RBCs. This removal could potentially result from targeting of the RBCs by the immune system. Specifically, upon rupture of the infected RBCs and release of the new brood of merozoites, the intracellular contents of the RBC and the parasites are released into the circulation and can adhere to uninfected RBCs, thereby potentially making them targets for erythrophagocytosis.

By contrast, an increased senescence mechanism would be more appropriate if the processes leading to the removal of uninfected RBCs were in fact age-related, for example, with an increased rate of senescence of the uninfected RBCs in the face of an infection. For instance, the immune response to the infection could lead to an increased level of oxygen radicals, which in turn could trigger an increased rate of RBC senescence.

The age-independent and the increased senescence mechanisms both fit the experimental data equally well (Fig. 7). However, they very clearly predict a different age distribution of the RBCs, in particular during the peak of parasitaemia (Fig. 7d). While the age-independent mechanism predicts almost no change in the age distribution of RBCs during peak parasitaemia, the increased senescence mechanism predicts a population of RBCs that appears to be much older. It has been reported that uninfected human RBCs co-cultured in vitro with *P. falciparum* infected RBCs have a higher proportion of older cells than same-donor control cultures [38]. Although a similar trend has been observed here (Fig. 7) based on in vivo data, the results [38] suggest a different magnitude of increase in older cells, which in turn implies that increased senescence is neither likely the only—nor a major—process leading to the increased removal of uninfected RBCs.

Thus, the age-independent mechanism appears to be more likely an explanation for quantifying the removal processes. Considering this mechanism, the model allows the differentiation between four distinct processes of RBC removal. Two of these processes are normal physiological processes of RBC loss occurring in healthy macaques, namely, age-dependent and age-independent ‘random’ removal of RBCs. The age-dependent process collectively captures the normal processes of RBC senescence, whereas the age-independent processes encompass all RBC losses that occur under normal physiological conditions but are not related to age. This process was modelled as proposed by Löffler’s group [16, 32, 33], where 10 % of all produced RBCs are destined to be lost by this ‘random’ process.

The remaining two distinct processes of RBC removal are directly and indirectly due to the malarial infection. These are the direct parasitization with the consequential destruction of the infected RBCs and the loss of uninfected RBCs as a bystander effect. RBCs removed by direct parasitization are those RBCs that are invaded by a merozoite for the production of the next generation of merozoites, and destroyed concomitantly with their release into the bloodstream. The estimation of the number of infected cells in the current model is based on the experimental parasitaemia measurements by Moreno et al. [14]. Uninfected RBC removal occurs during the malaria infection for reasons that are not well understood and clearly warrant further investigation, because this process dominates the manifestation of malarial anaemia, as exemplified by the case studied here of *P. coatneyi* infecting *M. mulatta*. In both malaria-naïve and semi-immune infected macaques, the loss of uninfected RBCs emerged as the leading process in RBC removal (76 and 67 % of total removal, respectively), whereas the actual destruction of RBCs by parasitization represented only 4 and 1 % (respectively) of the total removal. The remaining losses were 15 % to senescence, and 4 % to age-independent processes, in infected naïve macaques.

RBC production during the 30-day infection of the malaria-naïve macaques was increased by a factor of almost two, which however did not match the total number of cells lost during this period (Fig. 4). The cause of this discrepancy is reflected in the experimental data, as some of the macaques did not fully recover from their anaemic state within the 30-day experimental interval and after 30 days still exhibited decreased haemoglobin levels (Fig. 1). A much smaller difference between the removal and replenishment of RBCs was also inferred for the semi-immune macaques (Fig. 4), as these macaques were mostly able to recover from their anaemias (Fig. 1). Comparison of the inferred RBC production profiles with the experimentally determined percentage of circulating reticulocytes (Fig. 5) showed good agreement, which attests to the model’s ability to predict correctly when erythropoiesis or release of reticulocytes from the bone marrow was upregulated. Yet, closer analysis reveals that this upregulation does not kick in during the time period when RBCs are being actively removed, but only once parasitaemia decreases due to host immune responses or anti-malarial treatment. This discrepancy suggests that the parasite does not reduce erythropoietic output, but that it rather prevents effective erythropoiesis with the normal release of reticulocytes in response to anaemia; only after the parasite is cleared from circulation is the host able to properly respond to the anaemia. Clinically, the infected macaques (whether naïve or semi-immune) showed evidence of erythroid hyperplasia (expansion of erythroid progenitors) in their bone marrow, coincident with the observed peaks of parasitaemia. However, the erythroid hyperplasia was neither reflected by a high number of reticulocytes in the peripheral blood (experimentally measured in [14]) nor by an increased production of RBCs (inferred in this work), corroborating the hypothesis of parasite-induced ineffective erythropoiesis. Consistent with ineffective erythropoiesis, bone marrow biopsies in naïve macaques showed histological signs of dyserythropoiesis in the form of defective development of erythroid cells [14]. Relevantly, macaques that had been malaria naïve when infected exhibited erythropoietin levels that were elevated six-fold during periods of high parasitaemia, indicating that ineffective erythropoiesis is not mediated by the dysregulation (i.e., lack) of compensatory erythropoietin production [14].

It appears that this article is the first to quantify the percentage of removal of uninfected RBCs accurately in an in vivo primate malaria system. It is moreover likely that the results have implications for understanding anaemia in humans that is caused by *P. falciparum* (and perhaps other species of *Plasmodium*). Two earlier, somewhat indirect studies have attempted to quantify the removal of uninfected and parasitized RBCs in humans infected with *P. falciparum* [10, 11]. Jakeman et al. [10] used a discrete mathematical model to calculate the ratio of uninfected erythrocytes destroyed per infected erythrocyte in twelve neurosyphilis patients who were infected with different strains (McLendon and El Limon) of *P. falciparum* as a treatment (malariotherapy). The results showed that on average 8.52 (1.28–19.28) uninfected RBCs are removed per infected RBC [10]; thus, the loss of uninfected RBCs represents about 90 % of the total removal. Price et al. [11] compared the total loss of RBCs inferred from the decrease in haematocrit, with an estimation of the number of RBCs infected and destroyed by parasites inferred from the parasitaemia level. The authors concluded that parasitization was responsible for 7.9 % (6.2–9.6 %) of the RBCs lost in patients with uncomplicated falciparum malaria [11]. These results are in the same range as our average estimates of 95 % (90–97 %) and 99 % (97–99 %) loss of uninfected RBCs in malaria-naïve compared to semi-immune macaques infected with *P. coatneyi*, relative to the total loss of RBCs induced by the infection.

The physiological basis for the high loss of uninfected RBCs is not known. In this paper, two mechanisms are compared: one where any uninfected bystander RBC is prone to removal and one where the infection causes an increased rate of aging and thus the senescence-based removal of uninfected RBCs. If only the second mechanism is in effect, then the model predicts a rather old population of RBCs during high parasitaemia (Fig. 7c, d). This prediction is at odds with experimental evidence obtained in RBC in vitro cultures infected with *P. falciparum* (D10 strain), which suggested a slight (~10 %) decrease in young uninfected RBCs and a small (~5 %) increase in old uninfected RBCs in comparison with a non-infected control culture of RBCs from the same donor [38]. By direct inference, increased senescence would not seem to be the only—or even the predominant—mechanism responsible for the loss of uninfected RBCs during a malarial infection.

Among the age-independent processes that could be responsible for this removal of uninfected RBCs, it has been suggested that reduced deformability, increased oxidation of membranes, inflammatory insults, or surface deposition of parasite proteins like PfRSP-2 could be involved [39]. Also of interest, data obtained in a model of chronic anaemia in semi-immune BALB/c mice infected with the *P. berghei* ANKA strain showed that this increased rate of removal of uninfected RBCs can be delayed by the depletion of macrophages, thus suggesting an immunopathological process where CD4^{+} T cells may be involved [40].

The model put forth here does not account for the microvascular sequestration of parasitized RBCs during the trophozoite and schizont stages of development, a phenomenon that was not experimentally quantifiable with the data, which targeted specifically the dynamic assessment of peripheral blood and bone marrow samples [14]. The daily monitored, circulating ring-stage infected RBCs can be representative of the parasite burden, but this does not reflect the daily parasite load in the host. Including RBC sequestration in the models would slightly change the rate of death by direct parasitization of RBCs and decrease the rate of removal of uninfected RBCs through the bystander effect. However, the overall results would not change much. As an example, consider the case of malaria-naïve macaques where, during a *P. coatneyi* infection, 5 % of RBCs died by parasitization and 95 % were removed as uninfected RBCs. Even if an equal number of infected RBCs were sequestered in the tissue microvasculature and in the peripheral blood circulation, the levels of death by parasitization would merely rise to 10 % and reduce the removal of uninfected RBCs to 90 %. This scenario does not change the main conclusions of this study, and the removal of uninfected RBCs would still greatly outnumber the cell losses due to parasite invasion of the RBCs.

Another assumption of the model, which could be construed as unlikely to occur in humans and in macaques, is the occurrence of random destruction of RBCs as a normal physiological process of RBC loss. This process was here assumed to have a magnitude of 10 %. If the model had been developed without this process, the parameter *c* in the hazard function (Eq. 1) would be zero, and the parameters *a* and *b* would easily compensate for the absence of *c*. As a consequence, the survival curve shown in Fig. 3b (blue) would not exhibit the slight slope present between the RBC ages of 0 and 60 days, but would be essentially flat in this range. Nevertheless, given that the same number of cells would still have to be removed, as it is evident from the biotinylated RBC turnover data in Fig. 3a, the hazard function would not change in magnitude, but just in shape. Therefore, as far as the model calculations are concerned, the current removal by a random process would be subsumed into the removal by senescence, while the inferred extent of removal due to the parasite (directly and through the bystander effect) would remain the same. Hence, the blue and purple bars in Fig. 4 would be merged into one single senescence process which, however, would not change the main results and/or conclusions of this work which are predicated on the large difference between the extent of RBC removal by direct parasitization and by the bystander effect.

## Conclusions

A discrete recursive model of RBC dynamics in *M. mulatta* was developed that allowed, for the first time, the accurate quantification of the processes leading to the production and removal of RBCs in circulation under healthy conditions and during a *P. coatneyi* infection. The results support experimental evidence that removal of uninfected RBCs is the dominant process leading to malarial anaemia. The model suggests that the loss of uninfected RBCs is mostly due to an age-independent process, with only a small number of cells removed due to increased senescence. The mechanisms behind this high loss of uninfected RBCs during malaria are currently unknown and warrant further investigation.

## Declarations

### Authors’ contributions

Conceived and designed the experiments: LLF, EOV. Performed the experiments: LLF. Analysed the data: LLF, HAS, EOV. Wrote the paper: LLF, EOV. Provided expert knowledge, viewpoints, feedback, and manuscript contributions: AM, JWB, MRG. All authors read and approved the final manuscript.

### Acknowledgements

The authors would like to thank Drs. Po-Wei Chen and Sepideh Dolatshahi, Chet Joyner, Mojdeh Faraji, Daniel Olivença, and James Wade for constructive comments and discussions.

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and material

The datasets analysed during the current study are available from the Ref. [14]. Datasets are, however, also available from the authors of this and of Ref. [14] studies upon reasonable request. The data generated during the current study are available from the corresponding author on reasonable request.

### Funding

This project was funded in part by Federal funds from the US National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services under contract # HHSN272201200031C (PI: Mary R. Galinski), which supports the Malaria Host–Pathogen Interaction Center (MaHPIC), as well as the Office of Research Infrastructure Programs/OD P51OD011132 (formerly National Center for Research Resources P51RR000165).

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- WHO. World Malaria Report 2014. Geneva: World Health Organization; 2014. http://www.who.int/malaria/publications/world_malaria_report_2014/en/.
- Mangano VD, Modiano D. An evolutionary perspective of how infection drives human genome diversity: the case of malaria. Curr Opin Immunol. 2014;30C:39–47.View ArticleGoogle Scholar
- Taylor SM, Cerami C, Fairhurst RM. Hemoglobinopathies: slicing the Gordian knot of
*Plasmodium falciparum*malaria pathogenesis. PLoS Pathog. 2013;9:e1003327.View ArticlePubMedPubMed CentralGoogle Scholar - Nagel RL, Roth EF Jr. Malaria and red cell genetic defects. Blood. 1989;74:1213–21.PubMedGoogle Scholar
- Hutagalung R, Wilairatana P, Looareesuwan S, Brittenham GM, Aikawa M, Gordeuk VR. Influence of hemoglobin E trait on the severity of
*falciparum*malaria. J Infect Dis. 1999;179:283–6.View ArticlePubMedGoogle Scholar - May J, Evans JA, Timmann C, Ehmen C, Busch W, Thye T, et al. Hemoglobin variants and disease manifestations in severe
*falciparum*malaria. JAMA. 2007;297:2220–6.View ArticlePubMedGoogle Scholar - Ayi K, Min-Oo G, Serghides L, Crockett M, Kirby-Allen M, Quirt I, et al. Pyruvate kinase deficiency and malaria. N Engl J Med. 2008;358:1805–10.View ArticlePubMedGoogle Scholar
- Coatney GR, Collins WE, Warren M, Contacos PG. The primate malarias. Bethesda; National Institute of Allergy and Infectious Diseases (US); 2003.Google Scholar
- Collins WE, Jeffery GM, Roberts JM. A retrospective examination of anemia during infection of humans with
*Plasmodium vivax*. Am J Trop Med Hyg. 2003;68:410–2.PubMedGoogle Scholar - Jakeman GN, Saul A, Hogarth WL, Collins WE. Anaemia of acute malaria infections in non-immune patients primarily results from destruction of uninfected erythrocytes. Parasitology. 1999;119:127–33.View ArticlePubMedGoogle Scholar
- Price RN, Simpson JA, Nosten F, Luxemburger C, Hkirjaroen L, ter Kuile F, et al. Factors contributing to anemia after uncomplicated
*falciparum*malaria. Am J Trop Med Hyg. 2001;65:614–22.PubMedPubMed CentralGoogle Scholar - Vargas-Serrato E, Corredor V, Galinski MR. Phylogenetic analysis of CSP and MSP-9 gene sequences demonstrates the close relationship of
*Plasmodium coatneyi*to*Plasmodium knowlesi*. Infect Genet Evol. 2003;3:67–73.View ArticlePubMedGoogle Scholar - Singh B, Daneshvar C. Human infections and detection of
*Plasmodium knowlesi*. Clin Microbiol Rev. 2013;26:165–84.View ArticlePubMedPubMed CentralGoogle Scholar - Moreno A, Cabrera-Mora M, Garcia A, Orkin J, Strobert E, Barnwell JW, et al.
*Plasmodium coatneyi*in rhesus macaques replicates the multisystemic dysfunction of severe malaria in humans. Infect Immun. 2013;81:1889–904.View ArticlePubMedPubMed CentralGoogle Scholar - Perkins DJ, Were T, Davenport GC, Kempaiah P, Hittner JB, Ong’echa JM. Severe malarial anemia: innate immunity and pathogenesis. Int J Biol Sci. 2011;7:1427–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Fonseca LL, Voit EO. Comparison of mathematical frameworks for modeling erythropoiesis in the context of malaria infection. Math Biosci. 2015;270:224–36.View ArticlePubMedGoogle Scholar
- McKendrick AG. Applications of the mathematics to medical problems. Proc Edinburgh Math Soc. 1926;44:98.View ArticleGoogle Scholar
- Lotka AJ. Elements of physical biology. Baltimore: Williams and Wilkins; 1925.Google Scholar
- de Back DZ, Kostova EB, van Kraaij M, van den Berg TK, van Bruggen R. Of macrophages and red blood cells; a complex love story. Front Physiol. 2014;5:9.PubMedPubMed CentralGoogle Scholar
- Kipps T. The organization and structure of lymphoid tissues. In: Kaushansky K, Lichtman M, Kipps T, Seligsohn U, Prchal J, editors. Williams Hematology. New York: McGraw-Hill Medical; 2011. p. 75–84.Google Scholar
- Rifkind JM, Nagababu E. Hemoglobin redox reactions and red blood cell aging. Antioxid Redox Signal. 2013;18:2274–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Mohandas N, Gallagher PG. Red cell membrane: past, present, and future. Blood. 2008;112:3939–48.View ArticlePubMedPubMed CentralGoogle Scholar
- Schluter K, Drenckhahn D. Co-clustering of denatured hemoglobin with band 3: its role in binding of autoantibodies against band 3 to abnormal and aged erythrocytes. Proc Natl Acad Sci USA. 1986;83:6137–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Pantaleo A, Giribaldi G, Mannu F, Arese P, Turrini F. Naturally occurring anti-band 3 antibodies and red blood cell removal under physiological and pathological conditions. Autoimmun Rev. 2008;7:457–62.View ArticlePubMedGoogle Scholar
- Arashiki N, Kimata N, Manno S, Mohandas N, Takakuwa Y. Membrane peroxidation and methemoglobin formation are both necessary for band 3 clustering: mechanistic insights into human erythrocyte senescence. Biochemistry. 2013;52:5760–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lutz HU, Nater M, Stammler P. Naturally occurring anti-band 3 antibodies have a unique affinity for C3. Immunology. 1993;80:191–6.PubMedPubMed CentralGoogle Scholar
- Lutz HU, Stammler P, Fasler S. How naturally occurring anti-band 3 antibodies stimulate C3b deposition to senescent and oxidatively stressed red blood cells. Biomed Biochim Acta. 1990;49:S224–9.PubMedGoogle Scholar
- Kay MM. Band 3 and its alterations in health and disease. Cell Mol Biol. 2004;50:117–38.PubMedGoogle Scholar
- Gompertz B. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. Philos Trans R Soc Lond. 1825;115:513–83.View ArticleGoogle Scholar
- Weibull W. A statistical distribution function of wide applicability. J Appl Mech. 1951;18:293–7.Google Scholar
- Makeham WM. On the law of mortality and the construction of annuity tables. Assur Mag J Inst Actuar. 1860;8:301–10.Google Scholar
- Loeffler M, Pantel K, Wulff H, Wichmann HE. A mathematical-model of erythropoiesis in mice and rats. 1. Structure of the model. Cell Tissue Kinet. 1989;22:13–30.PubMedGoogle Scholar
- Schirm S, Engel C, Loeffler M, Scholz M. A biomathematical model of human erythropoiesis under erythropoietin and chemotherapy administration. PLoS ONE. 2013;8:e65630.View ArticlePubMedPubMed CentralGoogle Scholar
- Schirm S, Engel C, Loeffler M, Scholz M. A combined model of human erythropoiesis and granulopoiesis under growth factor and chemotherapy treatment. Theor Biol Med Model. 2014;11:24.View ArticlePubMedPubMed CentralGoogle Scholar
- Gross AJ, Clark V. Survival distributions: reliability applications in the biomedical sciences. New York: Wiley; 1975.Google Scholar
- Eyles DE, Fong YL, Warren M, Guinn E, Sandosham AA, Wharton RH.
*Plasmodium coatneyi*, a new species of primate malaria from Malaya. Am J Trop Med Hyg. 1962;11:597–604.Google Scholar - Beutler E. Destruction of erythrocytes. In: Kaushansky K, Lichtman M, Kipps T, Seligsohn U, Prchal J, editors. Williams hematology. New York: McGraw-Hill Medical; 2011. p. 449–54.Google Scholar
- Omodeo-Sale F, Motti A, Basilico N, Parapini S, Olliaro P, Taramelli D. Accelerated senescence of human erythrocytes cultured with
*Plasmodium falciparum*. Blood. 2003;102:705–11.View ArticlePubMedGoogle Scholar - Lamikanra AA, Brown D, Potocnik A, Casals-Pascual C, Langhorne J, Roberts DJ. Malarial anemia: of mice and men. Blood. 2007;110:18–28.View ArticlePubMedGoogle Scholar
- Evans KJ, Hansen DS, van Rooijen N, Buckingham LA, Schofield L. Severe malarial anemia of low parasite burden in rodent models results from accelerated clearance of uninfected erythrocytes. Blood. 2006;107:1192–9.View ArticlePubMedPubMed CentralGoogle Scholar

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.