Quantitative analysis of mechanisms that govern red blood cell age structure and dynamics during anaemia

Responses of hematopoietic cells after ionizing-irradiation in anemic adult medaka (Oryzias latipes)

PURPOSE

Hematopoietic tissues of vertebrates are highly radiation sensitive and the effects of ionizing radiation on the hematopoiesis have been studied in mammals and teleosts for decades. In this study, radiation responses in the kidney, the main hematopoietic organ in teleosts, were investigated in Japanese medaka (Oryzias latipes), which has been a model animal and a large body of knowledge has been accumulated in radiation biology.

METHODS

Kidney, the main hematopoietic tissue of adult medaka fish, was locally irradiated using proton and carbon ion beams irradiation system of Takasaki Ion Accelerator for Advanced Radiation Application (TIARA), QST, and the effects on peripheral blood cells and histology of the kidney were investigated.

RESULTS

When only kidneys were locally irradiated with proton or carbon ion beam (15 Gy), the hematopoietic cells in the irradiated kidney and cell density in the peripheral blood decreased 7 days after the irradiation in the same manner as after the whole-body irradiation with γ-rays (15 Gy). These results demonstrate that direct irradiation of the hematopoietic cells in the kidney induced cell death and/or cell cycle arrest and stopped the supply of erythroid cells. Then, the cell density in the peripheral blood recovered to the control level within 4 days and 7 days after the γ-ray and proton beam irradiation (15 Gy), respectively, while the cell density in the peripheral blood did not recover after the carbon ion beam irradiation (15 Gy). The hematopoietic cells in the irradiated kidneys temporarily decreased and recovered to the control level within 21 days after the γ-ray or proton beam irradiation (15 Gy), while it did not recover after the carbon ion beam irradiation (15 Gy). In contrast, the recovery of the cell density in the peripheral blood delayed when anemic medaka were irradiated 1 day after the administration of phenylhydrazine. With and without γ-ray irradiation, a large number of hematopoietic cells was still proliferating in the kidney 7 days after the anemia induction.

CONCLUSIONS

The results obtained strongly suggest that the hematopoietic stem cells in medaka kidney prioritize to proliferate and increase peripheral blood cells to eliminate anemia, even when they are damaged by high-dose irradiation.

Linking functional and molecular mechanisms of host resilience to malaria infection

It remains challenging to understand why some hosts suffer severe illnesses, while others are unscathed by the same infection. We fitted a mathematical model to longitudinal measurements of parasite and red blood cell density in murine hosts from diverse genetic backgrounds to identify aspects of within-host interactions that explain variation in host resilience and survival during acute malaria infection. Among eight mouse strains that collectively span 90% of the common genetic diversity of laboratory mice, we found that high host mortality was associated with either weak parasite clearance, or a strong, yet imprecise response that inadvertently removes uninfected cells in excess. Subsequent cross-sectional cytokine assays revealed that the two distinct functional mechanisms of poor survival were underpinned by low expression of either pro- or anti-inflammatory cytokines, respectively. By combining mathematical modelling and molecular immunology assays, our study uncovered proximate mechanisms of diverse infection outcomes across multiple host strains and biological scales.

MULTIPLE PERIODIC SOLUTIONS OF A WITHIN-HOST MALARIA INFECTION MODEL WITH TIME DELAY

In this paper, we study a within-host malaria infection model recently proposed by Schneider et al. in 2018. The stability and Hopf bifurcation analysis at the interior equilibrium are carried out, finding that the basic reproduction number plays a key role in the dynamics of the model, and incrementing the time delay will induce Hopf bifurcation at this equilibrium. The global extension of the local Hopf branch is further tracked numerically by the MatLab package DDE-BIFTOOL. Neimark-Sacker bifurcation of Poincaré map and period-doubling bifurcation of the bifurcated periodic solution are also detected, resulting in the existence of quasi-periodic and multiple periodic solutions, respectively. These results reveal that Hopf bifurcation will indeed bring about the rich dynamics of the model.

Differential drivers of intraspecific and interspecific competition during malaria-helminth co-infection

Various host and parasite factors interact to determine the outcome of infection. We investigated the effects of two factors on the within-host dynamics of malaria in mice: initial infectious dose and co-infection with a helminth that limits the availability of red blood cells (RBCs). Using a statistical, time-series approach to model the within-host ‘epidemiology’ of malaria, we found that increasing initial dose reduced the time to peak cell-to-cell parasite propagation, but also reduced its magnitude, while helminth co-infection delayed peak cell-to-cell propagation, except at the highest malaria doses. Using a mechanistic model of within-host infection dynamics, we identified dose-dependence in parameters describing host responses to malaria infection and uncovered a plausible explanation of the observed differences in single vs co-infections. Specifically, in co-infections, our model predicted a higher background death rate of RBCs. However, at the highest dose, when intraspecific competition between malaria parasites would be highest, these effects of co-infection were not observed. Such interactions between initial dose and co-infection, although difficult to predict a priori, are key to understanding variation in the severity of disease experienced by hosts and could inform studies of malaria transmission dynamics in nature, where co-infection and low doses are the norm.

Ecology of asynchronous asexual replication: the intraerythrocytic development cycle of Plasmodium berghei is resistant to host rhythms

BACKGROUND

Daily periodicity in the diverse activities of parasites occurs across a broad taxonomic range. The rhythms exhibited by parasites are thought to be adaptations that allow parasites to cope with, or exploit, the consequences of host activities that follow daily rhythms. Malaria parasites (Plasmodium) are well-known for their synchronized cycles of replication within host red blood cells. Whilst most species of Plasmodium appear sensitive to the timing of the daily rhythms of hosts, and even vectors, some species present no detectable rhythms in blood-stage replication. Why the intraerythrocytic development cycle (IDC) of, for example Plasmodium chabaudi, is governed by host rhythms, yet seems completely independent of host rhythms in Plasmodium berghei, another rodent malaria species, is mysterious.

METHODS

This study reports a series of five experiments probing the relationships between the asynchronous IDC schedule of P. berghei and the rhythms of hosts and vectors by manipulating host time-of-day, photoperiod and feeding rhythms.

RESULTS

The results reveal that: (i) a lack coordination between host and parasite rhythms does not impose appreciable fitness costs on P. berghei; (ii) the IDC schedule of P. berghei is impervious to host rhythms, including altered photoperiod and host-feeding-related rhythms; (iii) there is weak evidence for daily rhythms in the density and activities of transmission stages; but (iv), these rhythms have little consequence for successful transmission to mosquitoes.

CONCLUSIONS

Overall, host rhythms do not affect the performance of P. berghei and its asynchronous IDC is resistant to the scheduling forces that underpin synchronous replication in closely related parasites. This suggests that natural variation in the IDC schedule across species represents different parasite strategies that maximize fitness. Thus, subtle differences in the ecological interactions between parasites and their hosts/vectors may select for the evolution of very different IDC schedules.

Epidemiology and control of maedi-visna virus: Curing the flock

Maedi-visna (MV) is a complex lentiviral disease syndrome characterised by long immunological and clinical latencies and chronic progressive inflammatory pathology. Incurable at the individual level, it is widespread in most sheep-keeping countries, and is a cause of lost production and poor animal welfare. Culling seropositive animals is the main means of control, but it might be possible to manage virus transmission effectively if its epidemiology was better quantified. We derive a mathematical epidemiological model of the temporal distributions of seroconversion probabilities and estimate susceptibility, transmission rate and latencies in three serological datasets. We demonstrate the existence of epidemiological latency, which has not explicitly been recognised in the SRLV literaure. This time delay between infection and infectiousness apparently exceeds the delay between infection and seroconversion. Poor body condition was associated with more rapid seroconversion, but not with a higher probability of infection. We estimate transmission rates amongst housed sheep to be at about 1,000 times faster than when sheep were at grass, when transmission was negligible. Maternal transmission has only a small role in transmission, because lambs from infected ewes have a low probability of being infected directly by them, and only a small proportion of lambs need be retained to maintain flock size. Our results show that MV is overwhelmingly a disease of housing, where sheep are kept in close proximity. Prevalence of MV is likely to double each year from an initial low incidence in housed flocks penned in typically-sized groups of sheep (c. 50) for even a few days per year. Ewes kept entirely at grass are unlikely to experience transmission frequently enough for MV to persist, and pre-existing infection should die out as older ewes are replaced, thereby essentially curing the flock.

Unraveling Disease Pathophysiology with Mathematical Modeling

Modeling has enabled fundamental advances in our understanding of the mechanisms of health and disease for centuries, since at least the time of William Harvey almost 500 years ago. Recent technological advances in molecular methods, computation, and imaging generate optimism that mathematical modeling will enable the biomedical research community to accelerate its efforts in unraveling the molecular, cellular, tissue-, and organ-level processes that maintain health, predispose to disease, and determine response to treatment. In this review, we discuss some of the roles of mathematical modeling in the study of human physiology and pathophysiology and some challenges and opportunities in general and in two specific areas: in vivo modeling of pulmonary function and in vitro modeling of blood cell populations.

Adaptive phenotypic plasticity in malaria parasites is not constrained by previous responses to environmental change

Background and objectives

Phenotypic plasticity enables organisms to maximize fitness by matching trait values to different environments. Such adaptive phenotypic plasticity is exhibited by parasites, which experience frequent environmental changes during their life cycle, between individual hosts and also in within-host conditions experienced during infections. Life history theory predicts that the evolution of adaptive phenotypic plasticity is limited by costs and constraints, but tests of these concepts are scarce.

Methodology

Here, we induce phenotypic plasticity in malaria parasites to test whether mounting a plastic response to an environmental perturbation constrains subsequent plastic responses to further environmental change. Specifically, we perturb red blood cell resource availability to induce Plasmodium chabaudi to alter the trait values of several phenotypes underpinning within-host replication and between-host transmission. We then transfer parasites to unperturbed hosts to examine whether constraints govern the parasites’ ability to alter these phenotypes in response to their new in-host environment.

Results

Parasites alter trait values in response to the within-host environment they are exposed to. We do not detect negative consequences, for within-host replication or between-host transmission, of previously mounting a plastic response to a perturbed within-host environment.

Conclusions and implications

We suggest that malaria parasites are highly plastic and adapted to adjusting their phenotypes in response to the frequent changes in the within-host conditions they experience during infections. Our findings support the growing body of evidence that medical interventions, such as anti-parasite drugs, induce plastic responses that are adaptive and can facilitate the survival and potentially, drug resistance of parasites.

Lay Summary

Malaria parasites have evolved flexible strategies to cope with the changing conditions they experience during infections. We show that using such flexible strategies does not impact upon the parasites’ ability to grow (resulting in disease symptoms) or transmit (spreading the disease).

Plasticity and genetic variation in traits underpinning asexual replication of the rodent malaria parasite, Plasmodium chabaudi

BACKGROUND

The ability of malaria (Plasmodium) parasites to adjust investment into sexual transmission stages versus asexually replicating stages is well known, but plasticity in other traits underpinning the replication rate of asexual stages in the blood has received less attention. Such traits include burst size (the number of merozoites produced per schizont), the duration of the asexual cycle, and invasion preference for different ages of red blood cell (RBC).

METHODS

Here, plasticity [environment (E) effects] and genetic variation [genotype (G) effects] in traits relating to asexual replication rate are examined for 4 genotypes of the rodent malaria parasite Plasmodium chabaudi. An experiment tested whether asexual dynamics differ between parasites infecting control versus anaemic hosts, and whether variation in replication rate can be explained by differences in burst size, asexual cycle, and invasion rates.

RESULTS

The within-host environment affected each trait to different extents but generally had similar impacts across genotypes. The dynamics of asexual densities exhibited a genotype by environment effect (G×E), in which one of the genotypes increased replication rate more than the others in anaemic hosts. Burst size and cycle duration varied between the genotypes (G), while burst size increased and cycle duration became longer in anaemic hosts (E). Variation in invasion rates of differently aged RBCs was not explained by environmental or genetic effects. Plasticity in burst size and genotype are the only traits making significant contributions to the increase in asexual densities observed in anaemic hosts, together explaining 46.4% of the variation in replication rate.

CONCLUSIONS

That host anaemia induces several species of malaria parasites to alter conversion rate is well documented. Here, previously unknown plasticity in other traits underpinning asexual replication is revealed. These findings contribute to mounting evidence that malaria parasites deploy a suite of sophisticated strategies to maximize fitness by coping with, or exploiting the opportunities provided by, the variable within-host conditions experienced during infections. That genetic variation and genotype by environment interactions also shape these traits highlights their evolutionary potential. Asexual replication rate is a major determinant of virulence and so, understanding the evolution of virulence requires knowledge of the ecological (within-host environment) and genetic drivers of variation among parasites.

The Challenge of Quantifying Synchrony in Malaria Parasites

Malaria infection is often accompanied by periodic fevers, triggered by synchronous cycles of parasite replication within the host. The degree of synchrony in parasite development influences the efficacy of drugs and immune defenses and is therefore relevant to host health and infectiousness. Synchrony is thought to vary over the course of infection and across different host-parasite genotype or species combinations, but the evolutionary significance - if any - of this diversity remains elusive. Standardized methods are lacking, but the most common metric for quantifying synchrony is the percentage of parasites in a particular developmental stage. We use a heuristic model to show that this metric is often unacceptably biased. Methodological challenges must be addressed to characterize diverse patterns of synchrony and their consequences for disease severity and spread.

Adaptive plasticity in the gametocyte conversion rate of malaria parasites

Sexually reproducing parasites, such as malaria parasites, experience a trade-off between the allocation of resources to asexual replication and the production of sexual forms. Allocation by malaria parasites to sexual forms (the conversion rate) is variable but the evolutionary drivers of this plasticity are poorly understood. We use evolutionary theory for life histories to combine a mathematical model and experiments to reveal that parasites adjust conversion rate according to the dynamics of asexual densities in the blood of the host. Our model predicts the direction of change in conversion rates that returns the greatest fitness after perturbation of asexual densities by different doses of antimalarial drugs. The loss of a high proportion of asexuals is predicted to elicit increased conversion (terminal investment), while smaller losses are managed by reducing conversion (reproductive restraint) to facilitate within-host survival and future transmission. This non-linear pattern of allocation is consistent with adaptive reproductive strategies observed in multicellular organisms. We then empirically estimate conversion rates of the rodent malaria parasite Plasmodium chabaudi in response to the killing of asexual stages by different doses of antimalarial drugs and forecast the short-term fitness consequences of these responses. Our data reveal the predicted non-linear pattern, and this is further supported by analyses of previous experiments that perturb asexual stage densities using drugs or within-host competition, across multiple parasite genotypes. Whilst conversion rates, across all datasets, are most strongly influenced by changes in asexual density, parasites also modulate conversion according to the availability of red blood cell resources. In summary, increasing conversion maximises short-term transmission and reducing conversion facilitates in-host survival and thus, future transmission. Understanding patterns of parasite allocation to reproduction matters because within-host replication is responsible for disease symptoms and between-host transmission determines disease spread.

Malaria parasites in the host replicate asexually and, during each replication cycle, some asexuals transform into sexual stages that enable between-host transmission. It is not understood why the rate of conversion to sexual stages varies during infections despite its importance for the severity and spread of the disease. We combined a mathematical model and experiments to show that parasites adjust conversion rates depending on changes in their in-host population size. When population sizes plummet, between-host transmission is prioritised. However, smaller losses in number elicit reproductive restraint, which facilitates in-host survival and future transmission. We show that increased and decreased conversion in response to a range of in-host environments are actually part of one continuum: a sophisticated reproductive strategy similar to that of multicellular organisms.

Mitochondrial DNA is critical for longevity and metabolism of transmission stage Trypanosoma brucei

The sleeping sickness parasite Trypanosoma brucei has a complex life cycle, alternating between a mammalian host and the tsetse fly vector. A tightly controlled developmental programme ensures parasite transmission between hosts as well as survival within them and involves strict regulation of mitochondrial activities. In the glucose-rich bloodstream, the replicative 'slender' stage is thought to produce ATP exclusively via glycolysis and uses the mitochondrial F1FO-ATP synthase as an ATP hydrolysis-driven proton pump to generate the mitochondrial membrane potential (ΔΨm). The 'procyclic' stage in the glucose-poor tsetse midgut depends on mitochondrial catabolism of amino acids for energy production, which involves oxidative phosphorylation with ATP production via the F1FO-ATP synthase. Both modes of the F1FO enzyme critically depend on FO subunit a, which is encoded in the parasite's mitochondrial DNA (kinetoplast or kDNA). Comparatively little is known about mitochondrial function and the role of kDNA in non-replicative 'stumpy' bloodstream forms, a developmental stage essential for disease transmission. Here we show that the L262P mutation in the nuclear-encoded F1 subunit γ that permits survival of 'slender' bloodstream forms lacking kDNA ('akinetoplastic' forms), via FO-independent generation of ΔΨm, also permits their differentiation into stumpy forms. However, these akinetoplastic stumpy cells lack a ΔΨm and have a reduced lifespan in vitro and in mice, which significantly alters the within-host dynamics of the parasite. We further show that generation of ΔΨm in stumpy parasites and their ability to use α-ketoglutarate to sustain viability depend on F1-ATPase activity. Surprisingly, however, loss of ΔΨm does not reduce stumpy life span. We conclude that the L262P γ subunit mutation does not enable FO-independent generation of ΔΨm in stumpy cells, most likely as a consequence of mitochondrial ATP production in these cells. In addition, kDNA-encoded genes other than FO subunit a are important for stumpy form viability.

Altered life history strategies protect malaria parasites against drugs

Drug resistance has been reported against all antimalarial drugs, and while parasites can evolve classical resistance mechanisms (e.g., efflux pumps), it is also possible that changes in life history traits could help parasites evade the effects of treatment. The life history of malaria parasites is governed by an intrinsic resource allocation problem: specialized stages are required for transmission, but producing these stages comes at the cost of producing fewer of the forms required for within-host survival. Drug treatment, by design, alters the probability of within-host survival, and so should alter the costs and benefits of investing in transmission. Here, we use a within-host model of malaria infection to predict optimal patterns of investment in transmission in the face of different drug treatment regimes and determine the extent to which alternative patterns of investment can buffer the fitness loss due to drugs. We show that over a range of drug doses, parasites are predicted to adopt "reproductive restraint" (investing more in asexual replication and less in transmission) to maximize fitness. By doing so, parasites recoup some of the fitness loss imposed by drugs, though as may be expected, increasing dose reduces the extent to which altered patterns of transmission investment can benefit parasites. We show that adaptation to drug-treated infections could result in more virulent infections in untreated hosts. This work emphasizes that in addition to classical resistance mechanisms, drug treatment generates selection for altered parasite life history. Understanding how any shifts in life history will alter the efficacy of drugs, as well as any limitations on such shifts, is important for evaluating and predicting the consequences of drug treatment.

Examining the Reticulocyte Preference of Two Plasmodium berghei Strains during Blood-Stage Malaria Infection

The blood-stage of the Plasmodium parasite is one of the key phases within its life cycle that influences disease progression during a malaria infection. The efficiency of the parasite in infecting red blood cells (RBC) determines parasite load and parasite-induced hemolysis that is responsible for the development of anemia and potentially drives severe disease progression. However, the molecular factors defining the infectivity of Plasmodium parasites have not been completely identified so far. Using the Plasmodium berghei mouse model for malaria, we characterized and compared the blood-stage infection dynamics of PbANKA WT and a mutant parasite strain lacking a novel Plasmodium antigen, PbmaLS_05, that is well conserved in both human and animal Plasmodium parasite strains. Infection of mice with parasites lacking PbmaLS_05 leads to lower parasitemia levels and less severe disease progression in contrast to mice infected with the wildtype PbANKA strain. To specifically determine the effect of deleting PbmaLS_05 on parasite infectivity we developed a mathematical model describing erythropoiesis and malarial infection of RBC. By applying our model to experimental data studying infection dynamics under normal and drug-induced altered erythropoietic conditions, we found that both PbANKA and PbmaLS_05 (-) parasite strains differed in their infectivity potential during the early intra-erythrocytic stage of infection. Parasites lacking PbmaLS_05 showed a decreased ability to infect RBC, and immature reticulocytes in particular that are usually a preferential target of the parasite. These altered infectivity characteristics limit parasite burden and affect disease progression. Our integrative analysis combining mathematical models and experimental data suggests that deletion of PbmaLS_05 affects productive infection of reticulocytes, which makes this antigen a useful target to analyze the actual processes relating RBC preferences to the development of severe disease outcomes in malaria.

Phenotypic plasticity in reproductive effort: malaria parasites respond to resource availability

The trade-off between survival and reproduction is fundamental in the life history of all sexually reproducing organisms. This includes malaria parasites, which rely on asexually replicating stages for within-host survival and on sexually reproducing stages (gametocytes) for between-host transmission. The proportion of asexual stages that form gametocytes (reproductive effort) varies during infections-i.e. is phenotypically plastic-in response to changes in a number of within-host factors, including anaemia. However, how the density and age structure of red blood cell (RBC) resources shape plasticity in reproductive effort and impacts upon parasite fitness is controversial. Here, we examine how and why the rodent malaria parasite Plasmodium chabaudi alters its reproductive effort in response to experimental perturbations of the density and age structure of RBCs. We show that all four of the genotypes studied increase reproductive effort when the proportion of RBCs that are immature is elevated during host anaemia, and that the responses of the genotypes differ. We propose that anaemia (counterintuitively) generates a resource-rich environment in which parasites can afford to allocate more energy to reproduction (i.e. transmission) and that anaemia also exposes genetic variation to selection. From an applied perspective, adaptive plasticity in parasite reproductive effort could explain the maintenance of genetic variation for virulence and why anaemia is often observed as a risk factor for transmission in human infections.

Modelling primaquine-induced haemolysis in G6PD deficiency

Primaquine is the only drug available to prevent relapse in vivax malaria. The main adverse effect of primaquine is erythrocyte age and dose-dependent acute haemolytic anaemia in individuals with glucose-6-phosphate dehydrogenase deficiency (G6PDd). As testing for G6PDd is often unavailable, this limits the use of primaquine for radical cure. A compartmental model of the dynamics of red blood cell production and destruction was designed to characterise primaquine-induced haemolysis using a holistic Bayesian analysis of all published data and was used to predict a safer alternative to the currently recommended once weekly 0.75 mg/kg regimen for G6PDd. The model suggests that a step-wise increase in daily administered primaquine dose would be relatively safe in G6PDd. If this is confirmed, then were this regimen to be recommended for radical cure patients would not require testing for G6PDd in areas where G6PDd Viangchan or milder variants are prevalent.

DOI:http://dx.doi.org/10.7554/eLife.23061.001

Malaria is the most important parasitic disease that affects humans. Over half of the malaria cases in Asia and South America are caused by a species of malaria parasite called Plasmodium vivax (known as vivax malaria). This form of malaria results in repeated illness because dormant parasites in the liver wake at intervals to infect the blood. The only available drug that can stop these relapses is a drug called primaquine, which was developed seventy years ago.

Unfortunately, primaquine causes dangerous side effects in certain individuals who are deficient in an enzyme called G6PD, which helps defend red blood cells against stresses. Primaquine damages these cells so that they burst, leading to anaemia. This is a major problem because G6PD deficiency is common in regions where malaria is present: in some areas up to 30% of the population may be G6PD deficient. Since G6PD testing is not widely available, doctors often avoid prescribing primaquine to treat malaria, which results in more cases of disease relapse. Failing to prevent vivax relapses causes extensive illness and hinders efforts to eliminate malaria.

Is there a way to give this drug to patients that would be safer for people with G6PD deficiency? Primaquine destroys older rather than younger red blood cells. Watson et al. used mathematical modelling to see whether it is possible to develop a primaquine treatment strategy that would allow a gradual destruction of older red blood cells in individuals with G6PD deficiency, which would be safer. The mathematical model incorporates data from previous studies in malaria patients and healthy volunteers with G6PD deficiency and combines this with knowledge of how red blood cells are produced and destroyed. Watson et al. predicted that giving primaquine over 20 days in a steadily increasing dose was safer than current recommendations.

Mathematical models are simplifications of real world processes. The only way to test these findings properly will be to run a clinical trial that gives healthy volunteers who are G6PD deficient a course of primaquine treatment with a steadily increasing dose.

DOI:http://dx.doi.org/10.7554/eLife.23061.002

Gamma radiation at a human relevant low dose rate is genotoxic in mice

Even today, 70 years after Hiroshima and accidents like in Chernobyl and Fukushima, we still have limited knowledge about the health effects of low dose rate (LDR) radiation. Despite their human relevance after occupational and accidental exposure, only few animal studies on the genotoxic effects of chronic LDR radiation have been performed. Selenium (Se) is involved in oxidative stress defence, protecting DNA and other biomolecules from reactive oxygen species (ROS). It is hypothesised that Se deficiency, as it occurs in several parts of the world, may aggravate harmful effects of ROS-inducing stressors such as ionising radiation. We performed a study in the newly established LDR-facility Figaro on the combined effects of Se deprivation and LDR γ exposure in DNA repair knockout mice (Ogg1^−/−) and control animals (Ogg1^+/−). Genotoxic effects were seen after continuous radiation (1.4 mGy/h) for 45 days. Chromosomal damage (micronucleus), phenotypic mutations (Pig-a gene mutation of RBC^CD24−) and DNA lesions (single strand breaks/alkali labile sites) were significantly increased in blood cells of irradiated animals, covering three types of genotoxic activity. This study demonstrates that chronic LDR γ radiation is genotoxic in an exposure scenario realistic for humans, supporting the hypothesis that even LDR γ radiation may induce cancer.

Facilitation through altered resource availability in a mixed-species rodent malaria infection

A major challenge in disease ecology is to understand how co-infecting parasite species interact. We manipulate in vivo resources and immunity to explain interactions between two rodent malaria parasites, Plasmodium chabaudi and P. yoelii. These species have analogous resource-use strategies to the human parasites Plasmodium falciparum and P. vivax: P. chabaudi and P. falciparum infect red blood cells (RBC) of all ages (RBC generalist); P. yoelii and P. vivax preferentially infect young RBCs (RBC specialist). We find that: (1) recent infection with the RBC generalist facilitates the RBC specialist (P. yoelii density is enhanced ~10 fold). This occurs because the RBC generalist increases availability of the RBC specialist's preferred resource; (2) co-infections with the RBC generalist and RBC specialist are highly virulent; (3) and the presence of an RBC generalist in a host population can increase the prevalence of an RBC specialist. Thus, we show that resources shape how parasite species interact and have epidemiological consequences.

Mathematical model of susceptibility, resistance, and resilience in the within-host dynamics between a Plasmodium parasite and the immune system

We developed a coupled age-structured partial differential equation model to capture the disease dynamics during blood-stage malaria. The addition of age structure for the parasite population, with respect to previous models, allows us to better characterize the interaction between the malaria parasite and red blood cells during infection. Here we prove that the system we propose is well-posed and there exist at least two global states. We further demonstrate that the numerical simulation of the system coincides with clinically observed outcomes of primary and secondary malaria infection. The well-posedness of this system guarantees that the behavior of the model remains smooth, bounded, and continuously dependent on initial conditions; calibration with clinical data will constrain domains of parameters and variables to physiological ranges.

Serum pantetheinase/vanin levels regulate erythrocyte homeostasis and severity of malaria

Tissue pantetheinase, encoded by the VNN1 gene, regulates response to stress, and previous studies have shown that VNN genes contribute to the susceptibility to malaria. Herein, we evaluated the role of pantetheinase on erythrocyte homeostasis and on the development of malaria in patients and in a new mouse model of pantetheinase insufficiency. Patients with cerebral malaria have significantly reduced levels of serum pantetheinase activity (PA). In mouse, we show that a reduction in serum PA predisposes to severe malaria, including cerebral malaria and severe anemia. Therefore, scoring pantetheinase in serum may serve as a severity marker in malaria infection. This disease triggers an acute stress in erythrocytes, which enhances cytoadherence and hemolysis. We speculated that serum pantetheinase might contribute to erythrocyte resistance to stress under homeostatic conditions. We show that mutant mice with a reduced serum PA are anemic and prone to phenylhydrazine-induced anemia. A cytofluorometric and spectroscopic analysis documented an increased frequency of erythrocytes with an autofluorescent aging phenotype. This is associated with an enhanced oxidative stress and shear stress-induced hemolysis. Red blood cell transfer and bone marrow chimera experiments show that the aging phenotype is not cell intrinsic but conferred by the environment, leading to a shortening of red blood cell half-life. Therefore, serum pantetheinase level regulates erythrocyte life span and modulates the risk of developing complicated malaria.

The role of spatial organization of cells in erythropoiesis

Erythropoiesis, the process of red blood cell production, occurs mainly in the bone marrow. The functional unit of mammalian erythropoiesis, the erythroblastic island, consists of a central macrophage surrounded by adherent erythroid progenitor cells (CFU-E/Pro-EBs) and their differentiating progeny, the erythroblasts. Central macrophages display on their surface or secrete various growth or inhibitory factors that influence the fate of the surrounding erythroid cells. CFU-E/Pro-EBs have three possible fates: (a) expansion of their numbers without differentiation, (b) differentiation into reticulocytes that are released into the blood, (c) death by apoptosis. CFU-E/Pro-EB fate is under the control of a complex molecular network, that is highly dependent upon environmental conditions in the erythroblastic island. In order to assess the functional role of space coupled with the complex network behavior in erythroblastic islands, we developed hybrid discrete-continuous models of erythropoiesis. A model was developed in which cells are considered as individual physical objects, intracellular regulatory networks are modeled with ordinary differential equations and extracellular concentrations by partial differential equations. We used the model to investigate the impact of an important difference between humans and mice in which mature late-stage erythroblasts produce the most Fas-ligand in humans, whereas early-stage erythroblasts produce the most Fas-ligand in mice. Although the global behaviors of the erythroblastic islands in both species were similar, differences were found, including a relatively slower response time to acute anemia in humans. Also, our modeling approach was very consistent with in vitro culture data, where the central macrophage in reconstituted erythroblastic islands has a strong impact on the dynamics of red blood cell production. The specific spatial organization of erythroblastic islands is key to the normal, stable functioning of mammalian erythropoiesis, both in vitro and in vivo. Our model of a simplified molecular network controlling cell decision provides a realistic functional unit of mammalian erythropoiesis that integrates multiple microenvironmental influences within the erythroblastic island with those of circulating regulators of erythropoiesis, such as EPO and glucocorticosteroids, that are produced at remote sites.

Immune-mediated competition in rodent malaria is most likely caused by induced changes in innate immune clearance of merozoites

Malarial infections are often genetically diverse, leading to competitive interactions between parasites. A quantitative understanding of the competition between strains is essential to understand a wide range of issues, including the evolution of virulence and drug resistance. In this study, we use dynamical-model based Bayesian inference to investigate the cause of competitive suppression of an avirulent clone of Plasmodium chabaudi (AS) by a virulent clone (AJ) in immuno-deficient and competent mice. We test whether competitive suppression is caused by clone-specific differences in one or more of the following processes: adaptive immune clearance of merozoites and parasitised red blood cells (RBCs), background loss of merozoites and parasitised RBCs, RBC age preference, RBC infection rate, burst size, and within-RBC interference. These processes were parameterised in dynamical mathematical models and fitted to experimental data. We found that just one parameter , the ratio of background loss rate of merozoites to invasion rate of mature RBCs, needed to be clone-specific to predict the data. Interestingly, was found to be the same for both clones in single-clone infections, but different between the clones in mixed infections. The size of this difference was largest in immuno-competent mice and smallest in immuno-deficient mice. This explains why competitive suppression was alleviated in immuno-deficient mice. We found that competitive suppression acts early in infection, even before the day of peak parasitaemia. These results lead us to argue that the innate immune response clearing merozoites is the most likely, but not necessarily the only, mediator of competitive interactions between virulent and avirulent clones. Moreover, in mixed infections we predict there to be an interaction between the clones and the innate immune response which induces changes in the strength of its clearance of merozoites. What this interaction is unknown, but future refinement of the model, challenged with other datasets, may lead to its discovery.

Malaria infections often consist of more than one strain of the same parasitic species. Understanding the within-host competition between these various strains is essential to understanding the evolution and epidemiology of drug resistance in malarial infections. The infection process and the competition between strains involve complicated biological processes that are explained by various hypotheses. Mathematical models tested against experimental data provide quantitative measures to compare these hypotheses and enable us to discern the actual biological processes that contribute to the observed dynamics. We use a group of models against experimental data on rodent malaria to test various hypotheses. Such quantitative measures, in understanding rodent malaria, can be considered as a step towards understanding within-host parasite dynamics. Our work presented here demonstrates how confronting mathematical models with data allows the discovery of subtle and novel interactions between hosts and parasites that would be impractical to do in an experiment and allows the rejection of hypotheses that are incorrect. It is our contention that understanding the forces controlling within-host parasite dynamics in well-defined experimental model is a necessary step towards understanding these features in natural infections.

Quantitative analysis of porcine reproductive and respiratory syndrome (PRRS) viremia profiles from experimental infection: a statistical modelling approach

Porcine reproductive and respiratory syndrome (PRRS) is one of the most economically significant viral diseases facing the global swine industry. Viremia profiles of PRRS virus challenged pigs reflect the severity and progression of infection within the host and provide crucial information for subsequent control measures. In this study we analyse the largest longitudinal PRRS viremia dataset from an in-vivo experiment. The primary objective was to provide a suitable mathematical description of all viremia profiles with biologically meaningful parameters for quantitative analysis of profile characteristics. The Wood's function, a gamma-type function, and a biphasic extended Wood's function were fit to the individual profiles using Bayesian inference with a likelihood framework. Using maximum likelihood inference and numerous fit criteria, we established that the broad spectrum of viremia trends could be adequately represented by either uni- or biphasic Wood's functions. Three viremic categories emerged: cleared (uni-modal and below detection within 42 days post infection(dpi)), persistent (transient experimental persistence over 42 dpi) and rebound (biphasic within 42 dpi). The convenient biological interpretation of the model parameters estimates, allowed us not only to quantify inter-host variation, but also to establish common viremia curve characteristics and their predictability. Statistical analysis of the profile characteristics revealed that persistent profiles were distinguishable already within the first 21 dpi, whereas it is not possible to predict the onset of viremia rebound. Analysis of the neutralizing antibody(nAb) data indicated that there was a ubiquitous strong response to the homologous PRRSV challenge, but high variability in the range of cross-protection of the nAbs. Persistent pigs were found to have a significantly higher nAb cross-protectivity than pigs that either cleared viremia or experienced rebound within 42 dpi. Our study provides novel insights into the nature and degree of variation of hosts' responses to infection as well as new informative traits for subsequent genomic and modelling studies.

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

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

A biomathematical model of human erythropoiesis under erythropoietin and chemotherapy administration

Anaemia is a common haematologic side effect of dose-dense multi-cycle cytotoxic polychemotherapy requiring erythrocyte transfusions or erythropoietin (EPO) administration. To simulate the effectiveness of different EPO application schedules, we performed both modelling of erythropoiesis under chemotherapy and pharmacokinetic and dynamic modelling of EPO applications in the framework of a single comprehensive biomathematical model. For this purpose, a cell kinetic model of bone marrow erythropoiesis was developed that is based on a set of differential compartment equations describing proliferation and maturation of erythropoietic cell stages. The system is regulated by several feedback loops comprising those mediated by EPO. We added a model of EPO absorption after injection at different sites and a pharmacokinetic model of EPO derivatives to account for the effects of external EPO applications. Chemotherapy is modelled by a transient depletion of bone marrow cell stages. Unknown model parameters were determined by fitting the predictions of the model to data sets of circulating erythrocytes, haemoglobin, haematocrit, percentage of reticulocytes or EPO serum concentrations derived from the literature or cooperating clinical study groups. Parameter fittings resulted in a good agreement of model and data. Depending on site of injection and derivative (Alfa, Beta, Delta, Darbepoetin), nine groups of EPO applications were distinguished differing in either absorption kinetics or pharmacokinetics. Finally, eight different chemotherapy protocols were modelled. The model was validated on the basis of scenarios not used for parameter fitting. Simulations were performed to analyze the impact of EPO applications on the risk of anaemia during chemotherapy. We conclude that we established a model of erythropoiesis under chemotherapy that explains a large set of time series data under EPO and chemotherapy applications. It allows predictions regarding yet untested EPO schedules. Prospective clinical studies are needed to validate model predictions and to explore the feasibility and effectiveness of the proposed schedules.

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

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

Novel methods for quantifying individual host response to infectious pathogens for genetic analyses

We propose two novel approaches for describing and quantifying the response of individual hosts to pathogen challenge in terms of infection severity and impact on host performance. The first approach is a direct extension of the methodology for estimating group tolerance (the change in performance with respect to changes in pathogen burden in a host population) to the level of individuals. The second approach aims to capture the dynamic aspects of individual resistance and tolerance over the entire time course of infections. In contrast to the first approach, which provides a means to disentangle host resistance from tolerance, the second approach focuses on the combined effects of both characteristics. Both approaches provide new individual phenotypes for subsequent genetic analyses and come with specific data requirements. In particular, both approaches rely on the availability of repeated performance and pathogen burden measurements of individuals over the time course of one or several episodes of infection. Consideration of individual tolerance also highlights some of the assumptions hidden within the concept of group tolerance, indicating where care needs to be taken in trait definition and measurement.

Plasticity in parasite phenotypes: evolutionary and ecological implications for disease

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

Causes of variation in malaria infection dynamics: insights from theory and data

Parasite strategies for exploiting host resources are key determinants of disease severity (i.e., virulence) and infectiousness (i.e., transmission between hosts). By iterating the development of theory and empirical tests, we investigated whether variation in parasite traits across two genetically distinct clones of the rodent malaria parasite, Plasmodium chabaudi, explains differences in within-host infection dynamics and virulence. First, we experimentally tested key predictions of our earlier modeling work. As predicted, the more virulent genotype produced more progeny parasites per infected cell (burst size), but in contrast to predictions, invasion rates of red blood cells (RBCs) did not differ between the genotypes studied. Second, we further developed theory by confronting our earlier model with these new data, testing a new set of models that incorporate more biological realism, and developing novel theoretical tools for identifying differences between parasite genotypes. Overall, we found robust evidence that differences in burst sizes contribute to variation in dynamics and that differential interactions between parasites and host immune responses also play a role. In contrast to previous work, our model predicts that RBC age structure is not important for explaining dynamics. Integrating theory and empirical tests is a potentially powerful way of progressing understanding of disease biology.

Large-scale transcriptional profiling and functional assays reveal important roles for Rho-GTPase signalling and SCL during haematopoietic differentiation of human embryonic stem cells

Understanding the transcriptional cues that direct differentiation of human embryonic stem cells (hESCs) and human-induced pluripotent stem cells to defined and functional cell types is essential for future clinical applications. In this study, we have compared transcriptional profiles of haematopoietic progenitors derived from hESCs at various developmental stages of a feeder- and serum-free differentiation method and show that the largest transcriptional changes occur during the first 4 days of differentiation. Data mining on the basis of molecular function revealed Rho-GTPase signalling as a key regulator of differentiation. Inhibition of this pathway resulted in a significant reduction in the numbers of emerging haematopoietic progenitors throughout the differentiation window, thereby uncovering a previously unappreciated role for Rho-GTPase signalling during human haematopoietic development. Our analysis indicated that SCL was the 11th most upregulated transcript during the first 4 days of the hESC differentiation process. Overexpression of SCL in hESCs promoted differentiation to meso-endodermal lineages, the emergence of haematopoietic and erythro-megakaryocytic progenitors and accelerated erythroid differentiation. Importantly, intrasplenic transplantation of SCL-overexpressing hESC-derived haematopoietic cells enhanced recovery from induced acute anaemia without significant cell engraftment, suggesting a paracrine-mediated effect.

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

Immune clearance and resource limitation (via red blood cell depletion) shape the peaks and troughs of malaria parasitemia, which in turn affect disease severity and transmission. Quantitatively partitioning the relative roles of these effects through time is challenging. Using data from rodent malaria, we estimated the effective propagation number, which reflects the relative importance of contrasting within-host control mechanisms through time and is sensitive to the inoculating parasite dose. Our analysis showed that the capacity of innate responses to restrict initial parasite growth saturates with parasite dose and that experimentally enhanced innate immunity can affect parasite density indirectly via resource depletion. Such a statistical approach offers a tool to improve targeting of drugs or vaccines for human therapy by revealing the dynamics and interactions of within-host regulatory mechanisms.

Energy production and redox status of rat red blood cells after reticulocytosis induced by various treatments

Stimulated erythropoiesis and reticulocytosis can be induced by daily bleeding, or by phenylhydrazine (PHZ) treatment. We compared the in vivo effects of PHZ and bleeding treatment on haematological, energy and redox status parameters in red blood cells (RBC) of rats. The results showed that all followed haematological parameters were significantly lower in bleeding, compared to PHZ-treated rats. PHZ induced even 2.58-fold higher reticulocytosis as compared to bleeding treatment. Although PHZ induced higher reticulocytosis, respiration intensity and energy production was lower than in bleeding-induced reticulocytes. These alterations were the consequence of increased superoxide anion and peroxynitrite concentrations in PHZ-treated rats. Bleeding treatment resulted in increased activity of an antioxidative enzyme, superoxide dismutase. In conclusion, differences in these two experimental models for reticulocytosis may be used as tools for appropriate pharmacological testing of redox-active substances considering energy and redox processes, as well as apoptosis pathways.

Ribavirin-induced anemia in hepatitis C virus patients undergoing combination therapy

The current standard of care for hepatitis C virus (HCV) infection – combination therapy with pegylated interferon and ribavirin – elicits sustained responses in only ∼50% of the patients treated. No alternatives exist for patients who do not respond to combination therapy. Addition of ribavirin substantially improves response rates to interferon and lowers relapse rates following the cessation of therapy, suggesting that increasing ribavirin exposure may further improve treatment response. A key limitation, however, is the toxic side-effect of ribavirin, hemolytic anemia, which often necessitates a reduction of ribavirin dosage and compromises treatment response. Maximizing treatment response thus requires striking a balance between the antiviral and hemolytic activities of ribavirin. Current models of viral kinetics describe the enhancement of treatment response due to ribavirin. Ribavirin-induced anemia, however, remains poorly understood and precludes rational optimization of combination therapy. Here, we develop a new mathematical model of the population dynamics of erythrocytes that quantitatively describes ribavirin-induced anemia in HCV patients. Based on the assumption that ribavirin accumulation decreases erythrocyte lifespan in a dose-dependent manner, model predictions capture several independent experimental observations of the accumulation of ribavirin in erythrocytes and the resulting decline of hemoglobin in HCV patients undergoing combination therapy, estimate the reduced erythrocyte lifespan during therapy, and describe inter-patient variations in the severity of ribavirin-induced anemia. Further, model predictions estimate the threshold ribavirin exposure beyond which anemia becomes intolerable and suggest guidelines for the usage of growth hormones, such as erythropoietin, that stimulate erythrocyte production and avert the reduction of ribavirin dosage, thereby improving treatment response. Our model thus facilitates, in conjunction with models of viral kinetics, the rational identification of treatment protocols that maximize treatment response while curtailing side effects.

The treatment of HCV infection poses a major global health-care challenge today. The current standard of care, combination therapy with interferon and ribavirin, works in only about half of the patients treated. Because no alternatives are available yet for patients in whom combination therapy fails, identifying ways to improve response to combination therapy is critical. Increasing exposure to ribavirin does improve response but is associated with the severe side-effect, anemia. One way to maximize treatment response therefore is to increase ribavirin exposure to levels just below where anemia becomes intolerable. A second way is to supplement combination therapy with growth hormones, such as erythropoietin, that increase the production of red blood cells (erythrocytes) and compensate for ribavirin-induced anemia. Rational optimization of combination therapy thus relies on a quantitative description of ribavirin-induced anemia, which is currently lacking. Here, we develop a model of the population dynamics of erythrocytes in individuals exposed to ribavirin that quantitatively describes ribavirin-induced anemia. Model predictions capture several independent observations of ribavirin-induced anemia in HCV patients undergoing combination therapy, estimate the threshold ribavirin exposure beyond which anemia becomes intolerable, suggest guidelines for the usage of growth hormones, and facilitate rational optimization of therapy.

Red cell survival: relevance and mechanism involved

Red blood cells (RBCs) often have a short circulating half-life in hemodialysis patients, which increases the difficulty of achieving a stable hemoglobin level. Fluctuations in erythropoietin (EPO) levels contribute to this increased RBC turnover because a decline in the level of EPO triggers the preferential destruction of newly-formed RBC, a process termed neocytolysis. The RBCs that are released during the treatment of renal anemia are often hypochromic, with a low content of iron; these RBCs are vulnerable to rapid turnover because iron-deficiency affects RBCs in several ways, such as, increased exposure of the phagocytic signaling molecule phosphatidylserine, loss of deformability, and increased oxidative stress. Both EPO fluctuation and the release of iron-deficient RBCs are characteristic events occurring during the management of renal anemia, and the shorter RBC lifetime is a component of the large fluctuations in hemoglobin level seen in patients on hemodialysis.

Fitness costs of disrupting circadian rhythms in malaria parasites

Circadian biology assumes that biological rhythms maximize fitness by enabling organisms to coordinate with their environment. Despite circadian clocks being such a widespread phenomenon, demonstrating the fitness benefits of temporal coordination is challenging and such studies are rare. Here, we tested the consequences--for parasites--of being temporally mismatched to host circadian rhythms using the rodent malaria parasite, Plasmodium chabaudi. The cyclical nature of malaria infections is well known, as the cell cycles across parasite species last a multiple of approximately 24 h, but the evolutionary explanations for periodicity are poorly understood. We demonstrate that perturbation of parasite rhythms results in a twofold cost to the production of replicating and transmission stages. Thus, synchronization with host rhythms influences in-host survival and between-host transmission potential, revealing a role for circadian rhythms in the evolution of host-parasite interactions. More generally, our results provide a demonstration of the adaptive value of circadian rhythms and the utility of using an evolutionary framework to understand parasite traits.

Quantitative analysis of immune response and erythropoiesis during rodent malarial infection

Malarial infection is associated with complex immune and erythropoietic responses in the host. A quantitative understanding of these processes is essential to help inform malaria therapy and for the design of effective vaccines. In this study, we use a statistical model-fitting approach to investigate the immune and erythropoietic responses in Plasmodium chabaudi infections of mice. Three mouse phenotypes (wildtype, T-cell-deficient nude mice, and nude mice reconstituted with T-cells taken from wildtype mice) were infected with one of two parasite clones (AS or AJ). Under a Bayesian framework, we use an adaptive population-based Markov chain Monte Carlo method and fit a set of dynamical models to observed data on parasite and red blood cell (RBC) densities. Model fits are compared using Bayes' factors and parameter estimates obtained. We consider three independent immune mechanisms: clearance of parasitised RBCs (pRBC), clearance of unparasitised RBCs (uRBC), and clearance of parasites that burst from RBCs (merozoites). Our results suggest that the immune response of wildtype mice is associated with less destruction of uRBCs, compared to the immune response of nude mice. There is a greater degree of synchronisation between pRBC and uRBC clearance than between either mechanism and merozoite clearance. In all three mouse phenotypes, control of the peak of parasite density is associated with pRBC clearance. In wildtype mice and AS-infected nude mice, control of the peak is also associated with uRBC clearance. Our results suggest that uRBC clearance, rather than RBC infection, is the major determinant of RBC dynamics from approximately day 12 post-innoculation. During the first 2-3 weeks of blood-stage infection, immune-mediated clearance of pRBCs and uRBCs appears to have a much stronger effect than immune-mediated merozoite clearance. Upregulation of erythropoiesis is dependent on mouse phenotype and is greater in wildtype and reconstitited mice. Our study highlights the informative power of statistically rigorous model-fitting techniques in elucidating biological systems.

Deep-coverage rhesus red blood cell proteome: a first comparison with the human and mouse red blood cell

BACKGROUND

Macaques are the closest evolutionary relatives of humans routinely used in basic and applied biomedical research. Their genetic, physiological, immunological and metabolic similarity to humans, second only to that of the great apes, makes them invaluable models of human disease. These similarities also mean that macaques are often the only experimental models available for evaluating increasingly specific drugs in development, and as a proof-of-concept bridge can help reduce the numbers of compounds that fail in clinical pharmaceutical research. In vertebrates, red blood cells (RBCs) diseases are frequently severe as their role as sole gas transporter makes them indispensable to survival; much research has therefore focused on an in-depth understanding of the functioning of the RBC. RBCs also host malaria, babesia and other parasites. Recently, we presented an in-depth proteome for the human RBC and a comparative human/mouse RBC proteome.

MATERIAL AND METHODS

Here, we present directly comparable data for the human, mouse and rhesus RBC proteomes. All proteins were identified, validated and categorized in terms of sub-cellular localization, protein family and function and, in comparison with the human and mouse RBC, were classified as orthologues, family-related or unique. Splice isoforms were identified and polypeptides migrating with anomalous apparent molecular weights were grouped into putatively ubiquitinylated or partially degraded complexes.

RESULTS AND DISCUSSION

Overall there was close concordance between mouse, human and rhesus proteomes, confirming the unexpected RBC complexity. Several novel findings in the human and mouse proteomes have been confirmed here. This comparison sheds light on several open issues in RBC biology and provides a departure point for more comprehensive understanding of RBC function.