Research Article  Open Access
Effect of the Polydispersity of RBCs on the Recovery Rate of RBCs during the Removal of CPAs
Abstract
In the process of removing cryoprotectants from cryopreserved blood, the theoretically optimal operating condition, which is based on the assumption that the distribution of red blood cells is uniform, is often used to reduce or even avoid the hypotonic damage to cells. However, due to the polydispersity of cells, the optimal condition is actually not reliable. In this study, based on the discrete concept developed in our previous work, the effect of the polydispersity on the recovery rate of cells in the dilutionfiltration system was statistically investigated by assigning three random parameters, isotonic cell volume, cell surface area, and osmotically inactive cell volume, to cells in small units of blood. The results show that, due to the polydispersity, the real recovery rate deviates from the ideal value that is based on uniform distribution. The deviation significantly increases with the standard errors of cell parameters, and it can be also magnified by high cryoprotectant concentrations. Under the effect of polydispersity, the uniform distributionbased optimized blood or diluent flow rate is not perfect. In practice, one should adopt a more conservative blood or diluent flow rate so that the hypotonic damage to cells can be further reduced.
1. Introduction
The process of removing cryoprotectants (CPAs) is one of five important steps in the cryopreservation of red blood cells (RBCs) [1–8], where the damage to RBCs could happen due to the cell volume excursion induced by osmotic disequilibria. In the past decades, many methods have been proposed to improve the process, such as centrifugation [9–12], dialysis [13–15], and dilutionfiltration [16, 17].
To optimize the methods mentioned above, many theoretical models have also been developed [18–25]. In these models, almost without exception, the three parameters of RBCs, such as the cell volume at isotonic condition, the cell surface area, and the osmotically inactive cell volume, were all assumed to be identical. Thus, the mean values measured from experiments were applied. In the centrifugation method, a single cell was often studied to get optimal conditions for all cells. In the dialysisbased or dilutionfiltration method, the blood was divided into discrete units so that the effect of the randomness of RBCs entering the washing system on the osmotic damage to RBCs and the washing time necessary to remove CPAs could be investigated; however, the cells were still considered the same [20, 24, 25]. In fact, the polydispersity of RBCs, that is, the difference between RBCs, exists widely [26]. For example, the RBCs have different sizes and shapes not only at different growth stages [27, 28], but also at the same stage due to health, gender, and ethnic differences [29, 30]. In abnormal subjects, the RBCs always have various differences compared to the ones in normal subjects (e.g., the volume distribution of RBCs would have a left shift peak in microcytic anemia patients or a right shift peak in pernicious anemia patients [31, 32]). In addition, the osmotically inactive cell volumes are also different in different growth stages for human RBCs [33–35], presenting in the form of random distributions [36, 37]. Therefore, the reliability of these models only using the mean properties of cells to some extent is questionable. In practice, the polydispersity of RBCs could cause the failure in not only the accurate prediction of the recovery rate of RBCs but also the optimization of the washing process of CPAs. Thus, it needs to be noticed in intensive study.
In this work, the effect of the polydispersity of RBCs on the osmotic damage to RBCs or the recovery rate of RBCs in removing CPAs will be focused on. We will perform the study using the dilutionfiltration system (Figure 1) in our previous work as an example to show the effect of polydispersity. The three characteristic parameters of RBCs as random variables will be integrated into the mass transfer equations. Then, the distribution of the cell maximum volumes under the polydispersity will be statistically analyzed, and the change of the distribution caused by the variation of the standard errors will be further studied. Finally, the effects of operating conditions, including cell swelling limits, cryoprotectant concentrations, and blood and diluent flow rates, will be discussed, and the safe operating condition under the effect of polydispersity will be given. To the best of our knowledge, the effect of the polydispersity of RBCs is for the first time considered in the cryopreservation of RBCs. Thus, our study could provide a strategy not only to predict more accurately the recovery rate of cells but also to search safer operating conditions for removing cryoprotectants. Furthermore, the thought in this work could be also extended to other washing systems to analyze the effect of the polydispersity of cells.
2. Modeling
2.1. Discrete Concept Revisited
In our previous work, a discrete concept was developed for the dialysisbased or dilutionfiltration method to trace the volume variations of cells [15, 20, 24]. In the dilutionfiltration method, the blood or red blood cell suspension to be washed is divided into a certain number of units. These units randomly enter the system, extend their volumes due to the dilution in the dilution region, restore their original volumes in the filter, go back to the blood bag, and wait to be selected into the system again (Figure 1). In the circulation of units, the CPA inside cells is transported out due to the decrease of the CPA outside cells induced by the dilution, filtered out of the hemofilter, and removed along with the filtrate, and the cells in units experience swelling and shrinkage.
The discretization makes it possible to trace the detailed information related to all units, including the volume of units, cell volumes, and solute concentrations inside and outside cells in units, to analyze whether the cells suffer osmotic damage and thus to statistically estimate the recovery rate of cells under a given operating condition. Compared to our previous work [17, 24, 25], the novelty of this work is that we consider the effect of the polydispersity of RBCs on the optimized operating condition for safely removing CPAs.
In the dilutionfiltration system, the mass transfer of cryoprotectants from inside to outside cells in any unit can be unified and calculated separately. By means of a similar derivation method in the literature [24], the volume of RBCs and the intracellular concentration of CPAs can be calculated by the classic twoparameter model [40, 41] and the extracellular concentration of CPAs can be obtained according to mass conservation: where and are CPA concentrations inside and outside RBCs, respectively (mol/kg H_{2}O); and are NaCl concentrations inside and outside RBCs, respectively (mol/kg H_{2}O); and are CPA volumes inside and outside RBCs, respectively (μm^{3}). (m/Pa/s) and (m/s) are hydraulic permeability and solute permeability of the RBC membrane, respectively. is the RBC volume (μm^{3}), is the cell membrane area (μm^{2}), and (μm^{3}) refers to the osmotically inactive volume of the cells at the isotonic condition (, , and are random variables). is the universal gas constant (J/mol/K); is the absolute temperature (K); is the CPA partial molar volume (l/mol); and is the cell density. is the volume of the blood units (it equals in the entrance region, in the dilution region, and in the recirculation region and decreases from to in the filter [24]). The superscript 0 denotes the previous time. The subscripts 1 and 2 denote inside and outside cells, respectively. The subscripts and denote blood and diluent, respectively. The subscripts and denote CPA and NaCl, respectively.
In the blood bag, the mixing occurs because the number of blood units is very large. Under this situation, cell volume and intracellular solute concentrations can be described by (1), (2), and (4) whereas extracellular solute concentrations can be calculated by the following equations approximately [24]: In the dilutionfiltration system, at the dilution point (Figure 1), solute concentrations outside cells will change. Based on mass conservation, the relationships for them are as follows [24]: where the superscripts and denote units before and after the dilution point, respectively, and is the NaCl concentration in the dilution solution.
2.2. Polydispersity of RBCs
Once the blood is divided into a certain number of units, we need to assign three parameters of cells for all units (here, we assume that cells are the same in one unit but different in different units). Then, we have to face two questions: the first is what distributions the isotonic volume of cells, the surface area of cells, and the osmotically inactive volume of cells statistically meet, respectively; the second is what the correlation among these parameters is. In practice, the osmotically inactive volume of cells cannot be directly measured and can be only obtained by extrapolation in hypotonic expansion experiments [33]. In this work, to simplify the topic studied, we assume that it is independent of the other two parameters. As for the isotonic volume and surface area of cells, we theoretically confirm that they are subject to bivariate normal distribution, shown in the following paragraphs.
RBCs usually differ from each other in size and shape. This fact results in the disproportion between the isotonic volume of cells () and the surface area of cells (). In experiments, the isotonic volume of cells and the surface area of cells are both subjected to normal distribution; however, they are associated with some certain relationship. If the shape of RBCs is spherical, they should rigorously be subjected to the 2/3 power relationship, but the shape of RBCs varies from sphere to oblate spheroid in different status. In the least squares analysis of volume versus area from the blood sample data of healthy adults, researchers concluded that they meet linear (with a correlation coefficient 0.969 in the literature [42] or 0.943 in the literature [43]) rather than 2/3 power relationship [42, 43].
In our work, based on the pioneer work in the literature [39, 42], we assume that the joint distribution of the isotonic volume of cells and the surface area of cells is subjected to a bivariate normal distribution, as shown in Figure 2(a). The distribution parameters including mean value ( and ), standard error ( and ), and correlation coefficient () are estimated by maximum likelihood estimation, which are presented below (the subscripts and indicate isotonic cell volume and surface area, resp.):where is the number of observations. Then, the joint distribution probability density function (PDF) can be described in vector format as below: where is the parameter vector, represents the mean vector, and represents the covariance matrix.
In this work, the Pearson Chisquare test of goodness of fit is applied for hypothesis test of our bivariate normal distribution [44]. The null hypothesis is that the isotonic cell volume and surface area are subjected to bivariate normal distribution at 5% significance level.
First, we perform coordinate transformation for region division and frequency statistics. The inverse matrix of the covariance matrix () is diagonalized by orthogonal matrix : where is the diagonalized covariance matrix. If we introduce and set , then we can get To calculate the Pearson statistic , we divide the () plane into mutually disjoint regions ( is the frequency statistics in each region, and is the theoretical probability in each region). can be chosen as finite or infinite rectangle to simplify the calculation. Here, we use a serial of equidistant lines parallel to axis to divide the () plane into different regions, which is shown in Figure 2(b) (the multiple points indicate that a same value of isotonic cell volume and surface area is observed multiple times (usually ≥3) in the measurement).
Under the coordinate transformation, is counted in each region (multiple points counted as 3) and is calculated by where represents the cumulative distribution function of the standard normal distribution.
Then, we can finally get the Chisquare critical value (the degree of freedom is calculated by the expression . In this study, is the number of different regions divided from the () plane in Pearson Chisquare test and is the number of estimated parameters): The null hypothesis cannot be rejected at the 5% significance level and, thus, the isotonic cell volume and surface area are subjected to bivariate normal distribution with the estimated parameters. Therefore, the Pearson Chisquare test proves that the bivariate normal distribution is more suitable than the linear regression.
As for the osmotically inactive cell volume (), though it may depend on the cell stage or cell aging, there are not any reliable data or literature to show its direct connections with the other two parameters. Here an independent normaldistributed random parameter is assumed in this work. The distribution of is shown as follows: where the mean value () and standard error () are derived from the several groups of fitted values listed in the literature [33].
2.3. Simulation Strategy
The purpose of this work is to theoretically and statistically study the effect of the polydispersity of RBCs in the process of removing cryoprotectants. The strategy is as follows. Firstly, the blood is divided into a certain number of small units, which are traced and controlled by a program code to randomly enter the washing system. Secondly, the three parameters of RBCs are generated by another program code and then randomly assigned to all units. Then, the above twoparameter equations are used to calculate the volume variations of cells and the concentration variations of CPAs in all units. Finally, the effect of the polydispersity on the recovery rate of RBCs under various situations is studied statistically.
In the simulation, the cryoprotectant was glycerol, the diluent was the isotonic solution only containing NaCl (290 mOsm), and the blood was mimicked by the red blood cell suspension only containing NaCl and glycerol (NaCl: 290 mOsm; glycerol: 6500 mOsm; and hematocrit: 30%). In the dilutionfiltration system [17], the lumen volume of the plasma filter (Plasmflo AP05H/L, ASAHI Co., Japan), the dilution region volume, and the recirculation region volume were 85 mL, 5 mL, and 10 mL, respectively (the diameter of the tubing was 4 mm). Other parameters used are listed in Table 1.
 
, , and are random variables (mean ± STD). 
In this work, 200 mL of blood was divided into 5000 units so that the error caused by the division was less than 0.5% (the volume of units was 40 μL and cells were considered identical in one unit but different in different units). Five thousand sets of random parameters (i.e., isotonic cell volume, cell surface area, and osmotically inactive volume) were generated by a program code and assigned to cells in 5000 units, subjecting to the distributions mentioned above. In calculation, units were randomly selected to enter the system and the cell volume changes in all units were traced, using another program code. If the cell maximum volume in one unit exceeds the upper tolerance limit (here, the upper tolerance limit for completely avoiding the hypotonic damage was set to 1.53 according to the literature [45]), all cells in this unit will be marked and regarded as dead. Then, the recovery rate of cells was statistically obtained. In this study, the program codes, developed in FORTRAN, were used to simulate the washing process and perform the statistical analysis. For each operating condition, the simulation was repeated 6 times and a mean value was used.
3. Results and Discussion
3.1. Statistical Distribution of the Maximum Volumes of the RBCs
In the dilutionfiltration system, if all cells are assumed to be the same (i.e., the parameter distributions of cells are assumed to be uniform), when the optimized diluent flow rate is used, the maximum volumes of all cells are below the upper tolerance limit and the hypotonic damage to cells can be avoided completely (Figure 3(a)). However, due to the effect of the cell polydispersity, there are still some cells suffering from the hypotonic damage in practice (Figure 3(b); the standard errors of cell parameters are listed in Table 2). Moreover, the damage will increase with the standard error of no matter isotonic volume, surface area, or inactive volume (Figure 4 and Table 3). The effect of the osmotically inactive cell volume is relatively more significant. In addition, our results also show that the distribution of the cell maximum volume is normal under the condition of the polydispersity (Figure 4). Here, the optimized diluent flow rate is an allowable value set by system, which makes the maximum volume of the RBCs very close to the upper tolerance limit and is kept constant in all cycles.


(a)
(b)
(a)
(b)
(c)
Having established that the polydispersity results in the optimized condition deviating from the actual expectation, we wonder how to improve the condition to minimize the effect of the polydispersity. For the dilutionfiltration system, an effective method is to set a small or conservative cellswelling limit in searching the optimized diluent flow rate (Figure 5(a)). By doing so, the distribution of the cell maximum volume will move left, as shown in Figure 5(b), and the recovery rate of cells will increase. For example, when the blood flow rate is 100 mL/min, if the three standard errors are 12.7 μm^{3}, 13.8 μm^{2}, and 2.28 μm^{3}, respectively, and the expected recovery rate is higher than 95%, the cell swelling limit should be smaller than 1.50 , instead of 1.53 (Figure 5(a)). Corresponding to the change, the optimized diluent flow rate should be smaller than 19 mL/min (Figure 5(b)).
(a)
(b)
3.2. Effects of the RBC Polydispersity under Various Initial CPA Concentrations
Since the polydispersity can cause the decrease in the recovery rate of cells, even if the optimized condition is used, it is important to understand in practice when the effect of the polydispersity becomes significant. Our results show that the polydispersity decreases the recovery rates of cells, especially when the initial CPA concentration is high (the high concentration of CPAs causes the large volume change of cells in cycles and then the effect of the polydispersity is accumulated). Moreover, the difference of the recovery rates of cells between uniform and random distributions also increases with the initial CPA concentration (Figure 6). For example, when the three standard errors are 12.7 μm^{3}, 13.8 μm^{2}, and 2.28 μm^{3}, respectively, if the CPA concentration added to blood at the beginning of cryopreservation is increased from 5.5 mol/kg H_{2}O to 6.5 mol/kg H_{2}O, the difference will increase from 0.38% to 13.14%. Therefore, a more conservative cellswelling limit should be set for the dilutionfiltration system in practice when a higher CPA concentration needs to be removed.
3.3. Effects of the RBC Polydispersity under Various Blood and Diluent Flow Rates
In the dilutionfiltration system, diluent and blood flow rates are two controllable factors to reduce the hypotonic damage to cells. Under the assumption of uniform distribution, for a fixed diluent flow rate, the ideal recovery rate of cells first remains 0, then increases approximately linearly, and finally keeps 1 as the blood flow rate increases (Figure 7(a); the upper inflection point is called the uniform distributionbased optimal blood flow rate for a fixed diluent flow rate). However, due to the polydispersity, the change trend of the real recovery rate of cells presents a smooth curve rather than a broken line: when the blood flow rate is larger, the real recovery rate is less than the ideal one; when the blood flow rate is lower, the result is just the opposite (Figure 7). Similar to the effect of polydispersity under various blood flow rates, the deviation between ideal and real recovery rates of cells also exists under various diluent flow rates (Figure 7(b); the lower inflection point is called the uniform distributionbased optimal diluent flow rate for a fixed blood flow rate). Our results indicate that the uniform distributionbased optimal blood or diluent flow rate is inappropriate and can still cause cell loss because of the polydispersity of cells. Therefore, in practice, to reduce the effect of the polydispersity of cells, a more conservative condition (a higher blood flow rate or a lower diluent flow rate) should be used so that the expected recovery rate of cells can be obtained. In this work, a basic picture on the recovery rates of cells under various blood and diluent flow rates is shown in Figure 8, taking into account the polydispersity of cells. The picture clearly shows the nonlinear trend of the recovery rate of cells for a fixed blood or diluent flow rate and it can be used to practically guide the washing process of CPAs using the dilutionfiltration system.
(a)
(b)
4. Conclusions
In this work, we provide a strategy to statistically study the effect of the polydispersity of RBCs on the recovery rate of RBCs in removing CPAs from cryopreserved blood by assigning three random parameters to all cells, isotonic volume, surface area, and inactive volume. To confirm the effect of the polydispersity, the deviation between ideal and real recovery rates of RBCs in the dilutionfiltration system proposed in our previous work is investigated. The results show that due to the polydispersity, the real recovery rate deviates from the ideal one; moreover, the deviation significantly increases with the standard errors of cell parameters (if the three standard errors vary from 6.35 μm^{3}, 6.9 μm^{2}, and 1.14 μm^{3} to 19.05 μm^{3}, 20.7 μm^{2}, and 3.42 μm^{3}, resp., the deviation will increase from 6.36% to 19.72%). The high concentration of CPAs added at the beginning of the cryopreservation process can magnify the effect of polydispersity. If the concentration is increased from 5.5 mol/kg H_{2}O to 6.5 mol/kg H_{2}O, when the three standard errors are 12.7 μm^{3}, 13.8 μm^{2}, and 2.28 μm^{3}, respectively, the deviation will increase from 0.38% to 13.14%. Under the effect of polydispersity, the uniform distributionbased optimized blood or diluent flow rate is not perfect. In practice, one still needs to adopt a more conservative condition, that is, a blood flow rate higher than the optimized blood flow rate or a diluent flow rate lower than the optimized diluent flow rate, so that the expected recovery rate of cells can be obtained.
In this study, we only theoretically study the effect of the polydispersity of RBCs. To validate the simulation with the measured recovery rate of RBCs, one should first experimentally determine the distributions of the three parameters of RBCs, then measure the recovery rate of RBCs after using the dilutionfiltration system under a given operating condition, and finally compare the measured recovery rate of RBCs with the theoretical value. In practice, to establish a customized filtrationdilution system by considering the polydispersity of RBCs as well as the associated optimal blood and diluent flow rates, an extra storage for the table that includes the precalculated optimal operating conditions needs to be embedded into the system. In addition, it should be noted that the strategy developed here could be used for not only the cryoprotectant unloading process but also the cryoprotectant loading process.
Conflict of Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the Fundamental Research Funds for the Central Universities of China (WK2100000001), the Specialized Research Fund for the Doctoral Program of Higher Education of China (WJ2100230004), and the Natural Science Foundation of Anhui Province (BJ2100230008).
References
 B. Calmels, P. Houzé, J.C. Hengesse, T. Ducrot, C. Malenfant, and C. Chabannon, “Preclinical evaluation of an automated closed fluid management device: cytomate, for washing out DMSO from hematopoietic stem cell grafts after thawing,” Bone Marrow Transplantation, vol. 31, no. 9, pp. 823–828, 2003. View at: Publisher Site  Google Scholar
 R. Syme, M. Bewick, D. Stewart, K. Porter, T. Chadderton, and S. Glück, “The role of depletion of dimethyl sulfoxide before autografting: on hematologic recovery, side effects, and toxicity,” Biology of Blood and Marrow Transplantation, vol. 10, no. 2, pp. 135–141, 2004. View at: Publisher Site  Google Scholar
 C. R. Valeri, G. Ragno, P. van Houten et al., “Automation of the glycerolization of red blood cells with the highseparation bowl in the Haemonetics ACP 215 instrument,” Transfusion, vol. 45, no. 10, pp. 1621–1627, 2005. View at: Publisher Site  Google Scholar
 H. T. Meryman, “Cryopreservation of living cells: principles and practice,” Transfusion, vol. 47, no. 5, pp. 935–945, 2007. View at: Publisher Site  Google Scholar
 G. J. RuizDelgado, C. MancíasGuerra, E. L. TamezGómez et al., “Dimethyl sulfoxideinduced toxicity in cord blood stem cell transplantation: report of three cases and review of the literature,” Acta Haematologica, vol. 122, no. 1, pp. 1–5, 2009. View at: Publisher Site  Google Scholar
 G. B. Quan, Y. Han, M. X. Liu et al., “Addition of oligosaccharide decreases the freezing lesions on human red blood cell membrane in the presence of dextran and glucose,” Cryobiology, vol. 62, no. 2, pp. 135–144, 2011. View at: Publisher Site  Google Scholar
 K. A. Almansoori, V. Prasad, J. F. Forbes et al., “Cryoprotective agent toxicity interactions in human articular chondrocytes,” Cryobiology, vol. 64, no. 3, pp. 185–191, 2012. View at: Publisher Site  Google Scholar
 V. Pallotta, G. M. D'Amici, A. D'Alessandro, R. Rossetti, and L. Zolla, “Red blood cell processing for cryopreservation: from fresh blood to deglycerolization,” Blood Cells, Molecules, and Diseases, vol. 48, no. 4, pp. 226–232, 2012. View at: Publisher Site  Google Scholar
 C. R. Valeri, G. L. E. Ragno, R. Pivacek et al., “A multicenter study of in vitro and in vivo values in human RBCs frozen with 40percent (wt/vol) glycerol and stored after deglycerolization for 15 days at 4 degrees C in AS3: assessment of RBC processing in the ACP 215 (vol 41, pg 933, 2001),” Transfusion, vol. 42, no. 12, p. 1618, 2001. View at: Google Scholar
 C. C. M. Lelkens, F. Noorman, J. G. Koning et al., “Stability after thawing of RBCs frozen with the high and lowglycerol method,” Transfusion, vol. 43, no. 2, pp. 157–164, 2003. View at: Publisher Site  Google Scholar
 N. Bandarenko, J. Cancelas, E. L. Snyder et al., “Successful in vivo recovery and extended storage of additive solution (AS)5 red blood cells after deglycerolization and resuspension in AS3 for 15 days with an automated closed system,” Transfusion, vol. 47, no. 4, pp. 680–686, 2007. View at: Publisher Site  Google Scholar
 R. E. Lusianti, J. D. Benson, J. P. Acker, and A. Z. Higgins, “Rapid removal of glycerol from frozenthawed red blood cells,” Biotechnology Progress, vol. 29, no. 3, pp. 609–620, 2013. View at: Publisher Site  Google Scholar
 A. Zelman, D. Gisser, D. Smith, and R. Stephen, “Dialysis processing of cryopreserved red blood cells,” in Synthetic Membranes: Hyper and Ultrafiltration Uses, vol. 2, pp. 109–118, 1981. View at: Google Scholar
 F. Castino and S. R. Wickramasinghe, “Washing frozen red blood cell concentrates using hollow fibres,” Journal of Membrane Science, vol. 110, no. 2, pp. 169–180, 1996. View at: Publisher Site  Google Scholar
 W. Ding, J. Yu, E. Woods, S. Heimfeld, and D. Gao, “Simulation of removing permeable cryoprotective agents from cryopreserved blood with hollow fiber modules,” Journal of Membrane Science, vol. 288, no. 12, pp. 85–93, 2007. View at: Publisher Site  Google Scholar
 F. Arnaud, E. Kapnik, and H. T. Meryman, “Use of hollow fiber membrane filtration for the removal of DMSO from platelet concentrates,” Platelets, vol. 14, no. 3, pp. 131–137, 2003. View at: Publisher Site  Google Scholar
 X. Zhou, Z. Liu, Z. Shu et al., “A dilutionfiltration system for removing cryoprotective agents,” Journal of Biomechanical Engineering, vol. 133, no. 2, Article ID 021007, 2011. View at: Publisher Site  Google Scholar
 D. Y. Gao, J. Liu, C. Liu et al., “Prevention of osmotic injury to human spermatozoa during addition and removal of glycerol,” Human Reproduction, vol. 10, no. 5, pp. 1109–1122, 1995. View at: Google Scholar
 J. A. Gilmore, J. Liu, D. Y. Gao, and J. K. Critser, “Determination of optimal cryoprotectants and procedures for their addition and removal from human spermatozoa,” Human Reproduction, vol. 12, no. 1, pp. 112–118, 1997. View at: Publisher Site  Google Scholar
 W. Ding, S. Sun, C. Liu, S. Heimfeld, and D. Gao, “A random method for theoretical estimation of RBC osmotic damage in removing CPAs from cryopreserved blood with hollow fiber modules: closedloop blood flow mode,” Journal of Membrane Science, vol. 360, no. 12, pp. 17–25, 2010. View at: Publisher Site  Google Scholar
 W. Ding, X. Zhou, S. Heimfeld, J.A. Reems, and D. Gao, “A steadystate mass transfer model of removing CPAs from cryopreserved blood with hollow fiber modules,” Journal of Biomechanical Engineering, vol. 132, no. 1, 2010. View at: Publisher Site  Google Scholar
 J. D. Benson, C. C. Chicone, and J. K. Critser, “A general model for the dynamics of cell volume, global stability, and optimal control,” Journal of Mathematical Biology, vol. 63, no. 2, pp. 339–359, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 J. D. Benson, A. J. Kearsley, and A. Z. Higgins, “Mathematical optimization of procedures for cryoprotectant equilibration using a toxicity cost function,” Cryobiology, vol. 64, no. 3, pp. 144–151, 2012. View at: Publisher Site  Google Scholar
 L. Gong, W. Ding, Y. Ma, S. Sun, G. Zhao, and D. Gao, “Statistical estimation of red blood cell osmotic damage during cryoprotective agent removal from cryopreserved blood,” Biopreservation and Biobanking, vol. 11, no. 5, pp. 299–308, 2013. View at: Publisher Site  Google Scholar
 H. Y. Qiao, W. P. Ding, S. J. Sun, L. Q. Gong, and D. Y. Gao, “Theoretical optimization of the removal of cryoprotective agents using a dilutionfiltration system,” Biomedical Engineering Online, vol. 13, article 120, 2014. View at: Publisher Site  Google Scholar
 Y. M. Alomari, S. N. H. Sheikh Abdullah, R. Zaharatul Azma, and K. Omar, “Automatic detection and quantification of WBCs and RBCs using iterative structured circle detection algorithm,” Computational and Mathematical Methods in Medicine, vol. 2014, Article ID 979302, 17 pages, 2014. View at: Publisher Site  Google Scholar
 T. Shinozuka, Y. Miyata, S. Takei et al., “Changes in En(a) human red blood cell membranes during in vivo ageing,” Mechanisms of Ageing and Development, vol. 86, no. 1, pp. 27–37, 1996. View at: Publisher Site  Google Scholar
 S. C. Gifford, T. Yoshida, S. S. Shevkoplyas, and M. W. Bitensky, “A highresolution, doublelabeling method for the study of in vivo red blood cell aging,” Transfusion, vol. 46, no. 4, pp. 578–588, 2006. View at: Publisher Site  Google Scholar
 I. Sincer, A. Zorlu, M. B. Yilmaz et al., “Relationship between red cell distribution width and right ventricular dysfunction in patients with chronic obstructive pulmonary disease,” Heart and Lung: Journal of Acute and Critical Care, vol. 41, no. 3, pp. 238–243, 2012. View at: Publisher Site  Google Scholar
 S. K. Zalawadiya, V. Veeranna, S. S. Panaich, L. Afonso, and J. K. Ghali, “Gender and ethnic differences in red cell distribution width and its association with mortality among low risk healthy united state adults,” The American Journal of Cardiology, vol. 109, no. 11, pp. 1664–1670, 2012. View at: Publisher Site  Google Scholar
 J. D. Bessman and R. K. Johnson, “Erythrocyte volume distribution in normal and abnormal subjects,” Blood, vol. 46, no. 3, pp. 369–379, 1975. View at: Google Scholar
 J. Guitard, D. Ribes, N. Kamar et al., “Predictive factors for anemia within the first year after orthotopic liver transplantation,” Transplantation Proceedings, vol. 38, no. 7, pp. 2342–2345, 2006. View at: Publisher Site  Google Scholar
 O. Linderkamp and H. J. Meiselman, “Geometric, osmotic, and membrane mechanical properties of densityseparated human red cells,” Blood, vol. 59, no. 6, pp. 1121–1127, 1982. View at: Google Scholar
 A. K. Solomon, M. R. Toon, and J. A. Dix, “Osmotic properties of human red cells,” Journal of Membrane Biology, vol. 91, no. 3, pp. 259–273, 1986. View at: Publisher Site  Google Scholar
 S. C. Gifford, J. Derganc, S. S. Shevkoplyas, T. Yoshida, and M. W. Bitensky, “A detailed study of timedependent changes in human red blood cells: from reticulocyte maturation to erythrocyte senescence,” British Journal of Haematology, vol. 135, no. 3, pp. 395–404, 2006. View at: Publisher Site  Google Scholar
 C. J. Hunt, S. E. Armitage, and D. E. Pegg, “Cryopreservation of umbilical cord blood: 1. Osmotically inactive volume, hydraulic conductivity and permeability of CD34^{+} cells to dimethyl sulphoxide,” Cryobiology, vol. 46, no. 1, pp. 61–75, 2003. View at: Publisher Site  Google Scholar
 I. L. Lisovskaya, I. M. Shcherbachenko, R. I. Volkova, and V. P. Tikhonov, “Modulation of RBC volume distributions by oxidants (phenazine methosulfate and tertbutyl hydroperoxide): role of Gardos channel activation,” Bioelectrochemistry, vol. 73, no. 1, pp. 49–54, 2008. View at: Publisher Site  Google Scholar
 T. H. Papanek, The Water Permeability of the Human Erythrocyte in the Temperature Range +25°C to −10°C, Massachusetts Institute of Technology, Department of Mechanical Engineering, 1978.
 A. W. L. Jay, “Geometry of the human erythrocyte. I. Effect of albumin on cell geometry,” Biophysical Journal, vol. 15, no. 3, pp. 205–222, 1975. View at: Publisher Site  Google Scholar
 F. W. Kleinhans, “Membrane permeability modeling: kedemkatchalsky vs a twoparameter formalism,” Cryobiology, vol. 37, no. 4, pp. 271–289, 1998. View at: Publisher Site  Google Scholar
 I. I. Katkov, “A twoparameter model of cell membrane permeability for multisolute systems,” Cryobiology, vol. 40, no. 1, pp. 64–83, 2000. View at: Publisher Site  Google Scholar
 P. B. Canham and A. C. Burton, “Distribution of size and shape in populations of normal human red cells.,” Circulation Research, vol. 22, no. 3, pp. 405–422, 1968. View at: Publisher Site  Google Scholar
 E. K. Hoffmann and P. B. Dunham, “Membrane mechanisms and intracellular signalling in cell volume regulation,” International Review of Cytology, vol. 161, pp. 173–262, 1995. View at: Publisher Site  Google Scholar
 C. B. Bates and Naval Weapons Laboratory Dahlgren, The ChiSquare Test of Goodness of Fit for a Bivariate Normal Distribution, Defense Technical Information Center, 1966.
 A. K. Parpart, P. B. Lorenz, E. R. Parpart, J. R. Gregg, and A. M. Chase, “The osmotic resistance (fragility) of human red cells,” Journal of Clinical Investigation, vol. 26, no. 4, pp. 636–640, 1947. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Heyuan Qiao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.