Quantitative Analysis of Mechanisms That Govern Red Blood Cell Age Structure and Dynamics during Anaemia

Mathematical modelling has proven an important tool in elucidating and quantifying mechanisms that govern the age structure and population dynamics of red blood cells (RBCs). Here we synthesise ideas from previous experimental data and the mathematical modelling literature with new data in order to test hypotheses and generate new predictions about these mechanisms. The result is a set of competing hypotheses about three intrinsic mechanisms: the feedback from circulating RBC concentration to production rate of immature RBCs (reticulocytes) in bone marrow, the release of reticulocytes from bone marrow into the circulation, and their subsequent ageing and clearance. In addition we examine two mechanisms specific to our experimental system: the effect of phenylhydrazine (PHZ) and blood sampling on RBC dynamics. We performed a set of experiments to quantify the dynamics of reticulocyte proportion, RBC concentration, and erythropoietin concentration in PHZ-induced anaemic mice. By quantifying experimental error we are able to fit and assess each hypothesis against our data and recover parameter estimates using Markov chain Monte Carlo based Bayesian inference. We find that, under normal conditions, about 3% of reticulocytes are released early from bone marrow and upon maturation all cells are released immediately. In the circulation, RBCs undergo random clearance but have a maximum lifespan of about 50 days. Under anaemic conditions reticulocyte production rate is linearly correlated with the difference between normal and anaemic RBC concentrations, and their release rate is exponentially correlated with the same. PHZ appears to age rather than kill RBCs, and younger RBCs are affected more than older RBCs. Blood sampling caused short aperiodic spikes in the proportion of reticulocytes which appear to have a different developmental pathway than normal reticulocytes. We also provide evidence of large diurnal oscillations in serum erythropoietin levels during anaemia.     

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.     

NMR studies of intracellular free calcium, free magnesium and sodium in the guinea pig reticulocyte and mature red cell

During the maturation process reticulocytes lose their intracellular organelles and undergo changes in membrane lipid composition and ion transport properties. While several reports indicate differences in the levels of magnesium, sodium and calcium in reticulocytes and erythrocytes, controversy remains concerning the actual magnitude and direction of ionic alterations during reticulocyte maturation. One problem with all of these studies is that the techniques used are invasive and are limited to measuring only the total cell ion content. We have used 31P, 23Na and 19F nuclear magnetic resonance (NMR) spectroscopy to compare the intracellular free ion and phosphometabolite levels in guinea pig reticulocytes and mature red blood cells. In contrast to a sharply decreased concentration of ATP in erythrocytes in comparison to reticulocytes, the intracellular free magnesium, measured using 31P-NMR, was increased by about 65% upon maturation (150 mumol/l cell water in reticulocytes in comparison to 250 mumol/l cell water in erythrocytes). Sizeable but opposite changes in intracellular sodium (5.5 mumol/ml cells in reticulocytes vs. 8.5 mumol/ml cells in erythrocytes) and intracellular free calcium (99 nM vs. 31 nM in reticulocytes and mature red cells, respectively) were also observed, suggesting that alterations in the kinetics of membrane ion transport systems, accompanying changes in phospholipid and cholesterol content, occur during the process of red cell maturation. However, in contrast to dog red blood cells, there was no evidence for the presence of a Na+/Ca2+ exchanger in guinea pig reticulocytes or erythrocytes.     

A highly pathogenic simian/human immunodeficiency virus effectively produces infectious virions compared with a less pathogenic virus in cell culture

Background

The host range of human immunodeficiency virus (HIV) is quite narrow. Therefore, analyzing HIV-1 pathogenesis in vivo has been limited owing to lack of appropriate animal model systems. To overcome this, chimeric simian and human immunodeficiency viruses (SHIVs) that encode HIV-1 Env and are infectious to macaques have been developed and used to investigate the pathogenicity of HIV-1 in vivo. So far, we have many SHIV strains that show different pathogenesis in macaque experiments. However, dynamic aspects of SHIV infection have not been well understood. To fully understand the dynamic properties of SHIVs, we focused on two representative strains—the highly pathogenic SHIV, SHIV-KS661, and the less pathogenic SHIV, SHIV-#64—and measured the time-course of experimental data in cell culture.

Methods

We infected HSC-F with SHIV-KS661 and -#64 and measured the concentration of Nef-negative (target) and Nef-positive (infected) HSC-F cells, the total viral load, and the infectious viral load daily for 9 days. The experiments were repeated at two different multiplicities of infection, and a previously developed mathematical model incorporating the infectious and non-infectious viruses was fitted to the full dataset of each strain simultaneously to characterize the infection dynamics of these two strains.

Results and conclusions

We quantified virological indices including virus burst sizes and basic reproduction number of both SHIV-KS661 and -#64. Comparing the burst size of total and infectious viruses (viral RNA copies and TCID₅₀, respectively), we found that there was a statistically significant difference between the infectious virus burst size of SHIV-KS661 and -#64, while there was no significant difference between the total virus burst size. Furthermore, our analyses showed that the fraction of infectious virus among the produced SHIV-KS661 viruses, which is defined as the infectious viral load (TCID₅₀/ml) divided by the total viral load (RNA copies/ml), is more than 10-fold higher than that of SHIV-#64 during overall infection (i.e., for 9 days). Taken together, we conclude that the highly pathogenic SHIV produces infectious virions more effectively than the less pathogenic SHIV in cell culture.

     

Fixation probabilities when generation times are variable: the burst death model

Estimating the fixation probability of a beneficial mutation has a rich history in theoretical population genetics. Typically, to attain mathematical tractability, we assume that generation times are fixed, while the number of offspring per individual is stochastic. However, fixation probabilities are extremely sensitive to these assumptions regarding life history. In this article, we compute the fixation probability for a "burst-death" life-history model. The model assumes that generation times are exponentially distributed, but the number of offspring per individual is constant. We estimate the fixation probability for populations of constant size and for populations that grow exponentially between periodic population bottlenecks. We find that the fixation probability is, in general, substantially lower in the burst-death model than in classical models. We also note striking qualitative differences between the fates of beneficial mutations that increase burst size and mutations that increase the burst rate. In particular, once the burst size is sufficiently large relative to the wild type, the burst-death model predicts that fixation probability depends only on burst rate.     

Preferential invasion of malarial merozoites into young red blood cells

The preferential invasion of malarial merozoites into subpopulations of red blood cells (RBCs) in vivo and in vitro has been the subject of repeated discussions. In this paper, an attempt is made to summarize these discussions and to pinpoint the mechanism by which this preference could arise. The available data suggest that a relatively simple mechanism, related to the capability of the merozoite to rearrange the proteins of the cytoskeleton of the RBC may determine the invasion rate into mature versus very young RBCs (reticulocytes). There is no evidence for significant differences between mature RBCs and reticulocytes in the presence of membrane proteins which might play a role in receptor-ligand binding of merozoites to their host cell. Consequently, the concept of "reticulocyte preference" is left and the ability of penetrating both mature and immature RBCs, versus immature RBCs only, is given as an explanation for the presence of ringforms exclusively in reticulocytes as observed for several species of vivax-type malaria parasites. The possible consequences of preferential invasion for the infection (in vivo) and the culture (in vitro) of different plasmodial species are discussed.     

Plasmodium berghei: relative immunogenicity of infected reticulocytes and infected oxyphilic red blood cells

Hyperbleeding of mice 1 day before and 1 day after infection with Plasmodium berghei resulted in a more aggravated infection. Parasitemia rose significantly faster, but the mean survival time of these mice was not significantly different from control mice. At Day 5 of infection, parasites were almost exclusively in reticulocytes in contrast to control infections in which parasites were found in oxyphilic erythrocytes at Day 5 after infection. Purified parasitized reticulocytes taken from hyperbled mice at Day 5 after infection contained more young developmental parasite stages than purified parasitized oxyphilic erythrocytes taken from normal mice at Day 5 to 7 after infection. Parasitized reticulocytes were more readily opsonized by antibodies from immune serum when compared to parasitized oxyphilic red blood cells and when used to stimulate immune spleen cells the former were better stimulator cells than the latter. Results suggest either that parasitized reticulocytes are more immunogenic then parasitized oxyphilic red blood cells or that suspensions of parasitized reticulocytes contain more immunogenic parasite stages than suspensions of parasitized oxyphilic red blood cells.     

Low red cell production may protect against severe anemia during a malaria infection—Insights from modeling

The malaria parasite causes lysis of red blood cells, resulting in anemia, a major cause of mortality and morbidity. Intuitively, one would expect the production of red blood cells to increase in order to compensate for this loss. However, it has been observed that this response is weaker than would be expected. Furthermore, iron supplementation for iron deficient children in malaria endemic regions can paradoxically adversely affect the clinical outcome of malaria infection. A possible explanation may lie in the preference that some malaria parasites show for infecting immature red blood cells (reticulocytes). In the presence of a parasite preference for immature red cells, a rise in red cell production can ‘fuel the fire’ of infection by increasing the availability of the parasite's preferred target cell.We present a mathematical model of red blood cell production and infection in order to explore this hypothesis. We assess the effect of varying the reticulocyte replacement rate and preference of the parasite for reticulocytes on four key outcome measures assessing anemia and parasitemia.For a given level of parasite preference for reticulocytes we uncover an optimal erythropoietic response which minimizes disease severity. Increasing red blood cell production much above this optimum confers no benefit to the patient, and in fact can increase the degree of anemia and parasitemia. These conclusions are consistent with epidemiological studies demonstrating that both iron deficiency and anemia are protective against severe malaria, whilst iron supplementation in malaria endemic regions is with an increased number of malaria related adverse effects. Thus, suppression of red blood cell production, rather than being an unfortunate side effect of inflammation, may be a host protective effect against severe malarial anemia.     

Control of red blood cell mass during spaceflight

Data are reviewed from twenty-two astronauts from seven space missions in a study of red blood cell mass. The data show that decreased red cell mass in all astronauts exposed to space for more than nine days, although the actual dynamics of mass changes varies with flight duration. Possible mechanisms for these changes, including alterations in erythropoietin levels, are discussed.     

Cytoskeleton Remodeling Induces Membrane Stiffness and Stability Changes of Maturing Reticulocytes

Reticulocytes, the precursors of erythrocytes, undergo drastic alterations in cell size, shape, and deformability during maturation. Experimental evidence suggests that young reticulocytes are stiffer and less stable than their mature counterparts; however, the underlying mechanism is yet to be fully understood. Here, we develop a coarse-grained molecular-dynamics reticulocyte membrane model to elucidate how the membrane structure of reticulocytes contributes to their particular biomechanical properties and pathogenesis in blood diseases. First, we show that the extended cytoskeleton in the reticulocyte membrane is responsible for its increased shear modulus. Subsequently, we quantify the effect of weakened cytoskeleton on the stiffness and stability of reticulocytes, via which we demonstrate that the extended cytoskeleton along with reduced cytoskeleton connectivity leads to the seeming paradox that reticulocytes are stiffer and less stable than the mature erythrocytes. Our simulation results also suggest that membrane budding and the consequent vesiculation of reticulocytes can occur independently of the endocytosis-exocytosis pathway, and thus, it may serve as an additional means of removing unwanted membrane proteins from reticulocytes. Finally, we find that membrane budding is exacerbated when the cohesion between the lipid bilayer and the cytoskeleton is compromised, which is in accord with the clinical observations that erythrocytes start shedding membrane surface at the reticulocyte stage in hereditary spherocytosis. Taken together, our results quantify the stiffness and stability change of reticulocytes during their maturation and provide, to our knowledge, new insights into the pathogenesis of hereditary spherocytosis and malaria.     

Red blood cells derived from peripheral blood and bone marrow CD34+ human haematopoietic stem cells are permissive to Plasmodium parasites infection

The production of fully functional human red cells in vitro from haematopoietic stem cells (hHSCs) has been successfully achieved. Recently, the use of hHSCs from cord blood represented a major improvement to develop the continuous culture system for Plasmodium vivax. Here, we demonstrated that CD34+hHSCs from peripheral blood and bone marrow can be expanded and differentiated to reticulocytes using a novel stromal cell. Moreover, these reticulocytes and mature red blood cells express surface markers for entrance of malaria parasites contain adult haemoglobin and are also permissive to invasion by P. vivax and Plasmodium falciparum parasites.     

Molecular forms of red blood cell hexokinase

Summary Mammalian red blood cell hexokinase has been shown to exist in two or more distinct molecular forms, which are separable by ion-exchange chromatography. Of these forms just one corresponds to hexokinase type I from other tissues, while the others differ from any previously reported hexokinase isozyme. Analysis of several molecular properties of the three major forms (la, Ib and Ic in the order of their elution from DE-52 columns) of hexokinase prepared from human red cells and of the two forms purified from rabbit reticulocytes, shows significant differences in the isoelectric point. The kinetic and regulatory characteristics, the molecular weight, the temperature and pH-dependence of the various isozymes were similar.The hexokinase isozymic pattern is largely dependent upon red blood cell age. Among all, hexokinase Ib is the predominant form in rabbit reticulocytes and becomes the minor component in the older cells; a similar situation has also been found in the human erythrocyte. At present the molecular basis of hexokinase heterogeneity remains unknown, however preliminary experimental findings indicate a post-translational modification as a possible mechanism.     

Differential effect of neuraminidase-treatment on the surface charge-associated properties of rat reticulocytes and erythrocytes. Studies by partitioning in two-polymer aqueous phases

Rat reticulocytes undergo charge-associated surface changes, detectable by cell partitioning in charged dextran-poly(ethylene glycol) aqueous phase systems, as they become mature erythrocytes. Young reticulocytes have a lower partition coefficient, i.e., quantity of cells in the top phase as a percentage of total cells added, than do mature erythrocytes. Sialic acid is the main charge-bearing group on red blood cells and, in the case of the rat, most of the sialic acid can be removed by treatment of the cells with neuraminidase (Vibrio cholerae). By combining isotopic 59Fe-labeling of reticulocytes with countercurrent distribution of the entire red blood cell population in charged dextran-poly(ethylene glycol) aqueous phases we have now studied the relative effect of neuraminidase-treatment on rat reticulocytes and mature erythrocytes. It was found that neuraminidase-treatment (a) does not eliminate surface differences, detectable by partitioning, between rat reticulocytes and erythrocytes and (b) reduces the partition coefficient of mature erythrocytes to a greater extent than the partition coefficient of reticulocytes indicating a differential effect of this enzyme on the two cell populations.     

A reticulocyte-binding protein complex of Plasmodium vivax merozoites.

Plasmodium vivax merozoites primarily invade reticulocytes. The basis of this restricted host cell preference has been debated. Here we introduce two novel P. vivax proteins that comigrate on reducing SDS-polyacrylamide gels, colocalize at the apical pole of merozoites, and adhere specifically to reticulocytes. The genes encoding these proteins, P. vivax reticulocyte-binding proteins 1 and 2 (PvRBP-1 and PvRBP-2), have been cloned and analyzed. Homologous genes are evident in the closely related simian malaria parasite, P. cynomolgi, which also prefers to invade reticulocytes, but are not evident in the genome of another related simian malaria parasite, P. knowlesi, which invades all red blood cell subpopulations. Native PvRBP-1 is likely a transmembrane-anchored disulfide-linked protein, and along with PvRBP-2 may function as an adhesive protein complex. We propose that the RBPs of P. vivax, and homologous proteins of P. cynomolgi, function to target the reticulocyte subpopulation of red blood cells for invasion.     

Differential effect of neuraminidase-treatment on the surface charge-associated properties of rat reticulocytes and erythrocytes: Studies by partitioning in two-polymer aqueous phases

Rat reticulocytes undergo charge-associated surface changes, detectable by cell partitioning in charged dextran-poly(ethylene glycol) aqueous phase systems, as they become mature erythrocytes. Young reticulocytes have a lower partition coefficient, i.e., quantity of cells in the top phase as a percentage of total cells added, than do mature erythrocytes. Sialic acid is the main charge-bearing group on red blood cells and, in the case of the rat, most of the sialic acid can be removed by treatment of the cells with neuraminidase (Vibrio cholerae). By combining isotopic ⁵⁹Fe-labeling of reticulocytes with countercurrent distribution of the entire red blood cell population in charged dextran-poly(ethylene glycol) aqueous phases we have now studied the relative effect of neuraminidase-treatment on rat reticulocytes and mature erythrocytes. It was found that neuraminidase-treatment (a) does not eliminate surface differences, detectable by partitioning, between rat reticulocytes and erythrocytes and (b) reduces the partition coefficient of mature erythrocytes to a greater extent than the partition coefficient of reticulocytes indicating a differential effect of this enzyme on the two cell populations.     

Exosomes from Plasmodium yoelii-infected reticulocytes protect mice from lethal infections

Exosomes are 30-100-nm membrane vesicles of endocytic origin that are released after the fusion of multivesicular bodies (MVBs) with the plasma membrane. While initial studies suggested that the role of exosomes was limited to the removal of proteins during the maturation of reticulocytes to erythrocytes, recent studies indicate that they are produced by different types of cells and are involved in promoting inter-cellular communication and antigen presentation. Here, we describe the isolation and characterization of exosomes from peripheral blood of BALB/c mice infected with the reticulocyte-prone non-lethal Plasmodium yoelii 17X strain. Importantly, proteomic analysis revealed the presence of parasite proteins in these vesicles. Moreover, immunization of mice with purified exosomes elicited IgG antibodies capable of recognizing P. yoelii-infected red blood cells. Furthermore, lethal challenge of immunized mice with the normocyte-prone lethal P. yoelii 17XL strain caused a significant attenuation in the course of parasitaemia, increased survival time, and altered the cell tropism to reticulocytes. These results were obtained also when the exosomes were isolated from a P. yoelii-infected reticulocyte culture indicating that reticulocyte-derived exosomes carry antigens and are involved in immune modulation. Moreover, inclusion of CpG ODN 1826 in exosome immunizations elicited IgG2a and IgG2b antibodies and promoted survival, clearance of parasites and subsequent sterile protection of 83% of the animals challenged with P. yoelli 17XL. To our knowledge, this is the first report of immune responses elicited by exosomes derived from reticulocytes opening new avenues for the modulation of anti-malaria responses.     

Fixation probabilities depend on life history: fecundity, generation time and survival in a burst-death model

The burst-death model has been developed to describe the life history of organisms with variable generation times and a burst of a fixed number of offspring. The model also includes an optional constant clearance rate, such as washout from a chemostat, and the possibility of sustained periods of population growth followed by severe bottlenecks, as in serial passaging. In this model, a beneficial mutation can either increase the burst rate or the burst size, or reduce the clearance rate, thus increasing survival. In this article we examine the effects of these three possible mechanisms on both the Malthusian fitness and the fixation probability of the lineage. We find that equivalent relative increases in the burst rate or burst size confer equivalent increases in the Malthusian fitness of a lineage, whereas increasing survival typically has a more moderate effect on Malthusian fitness. In contrast, for beneficial mutations that confer the same increase in fitness, mutations that increase survival are the most likely to fix, followed by mutations that increase the burst rate. Mutations that increase the burst size are the least likely to fix. These results imply that mutant lineages with the highest Malthusian fitness are not, in many cases, the most likely to escape extinction.     

Macroxyproteinase (M.O.P.): a 670 kDa proteinase complex that degrades oxidatively denatured proteins in red blood cells

Erythrocytes and reticulocytes are shown to undergo rapid rates of protein degradation following exposure to oxidative stress. Experiments with ATP depletion revealed that, unlike the proteolysis of many other abnormal proteins, the degradation of oxidatively modified proteins is an ATP-independent process. Ion exchange chromatography (DEAE Sepharose CL-6B), ammonium sulfate precipitation, gel filtration chromatography (Sephacryl S-300 or Sepharose CL-6B), and a second ion exchange step were used to resolve the activity responsible for degrading oxidatively modified proteins from (dialyzed) cell-free extracts of erythrocytes and reticulocytes. Gel filtration studies revealed that some 70-80% of the activity in erythrocytes, and some 60-70% of the activity in reticulocytes, is expressed by a 670 kDa proteinase complex that is not stimulated by ATP (in fact, ATP is slightly inhibitory). This proteinase complex is inhibited by sulfhydryl reagents, serine reagents, and transition metal chelators, and has a pH optimum of 7.8. We propose the trivial name "macroxyproteinase" or "M.O.P." (abbreviated from Macro-Oxy-Proteinase) for the complex because of its large size, substrate preference (oxidatively modified proteins), and inhibitor profile (which indicates multiple catalytic sites). Electrophoresis studies of the 670 kDa M.O.P. complex revealed the presence of 8 distinct polypeptide subunits with the following apparent molecular sizes: 21.5, 25.3, 26.2, 28.1, 30.0, 31.9, 33.3, and 35.7 kDa. The large molecular size of the M.O.P. complex, its ATP- and ubiquitin-independence, its inhibitor profile, its distinctive subunit banding pattern in denaturing electrophoresis gels, its pH optimum, and its proteolytic profile with fluorogenic peptide substrates all indicate that M.O.P. is identical to 600-700 kDa neutral/alkaline proteinase complexes that have been isolated from a wide variety of eucaryotic cells and tissues, but for which no function has previously been clear. We propose that macroxyproteinase is responsible for catalyzing most of the selective degradation of oxidatively denatured proteins in red blood cells. We further suggest that M.O.P. may perform the same function in other eucaryotic cells and tissues.