Abstract
We investigated to what extent a prediction of the ‘ideal’ hematocrit based on individual hemorheological profile with an equation of viscosity is relevant in trained athletes, and how the agreement between theoretical and actual values is modified by changes in training volume and performance. Elite soccer players (national level: 18–32 yr, weight 61–83 kg, body mass index 20.9–25.8 kg/m2) were seen twice at one year interval. Hemorheologic parameters were measured with the MT90 viscometer and the Myrenne aggregometer the theoretical bell-shaped curve of hematocrit/viscosity ratio as a function of hematocrit was reconstructed with Quemada’s equation using actual plasma viscosity and red cell rigidity to predict hematocrit/viscosity at various hematocrit levels. RBC aggregation is correlated at baseline with fat mass (M1 = 0.552 p < 0.02) and changes in aggregation are related to changes in fat mass (M = 0.652, p < 0.05; M1 = 0.647, p < 0.05). Predicted and actual hematocrit are correlated (r = 0.644, p < 0.05) but exhibit discrepancies (mean difference –1% range [3.24 to 1.24]) and those discrepancies are inversely correlated to the level of predicted hematocrit (r = –0.912, p < 0.01), to systolic blood pressure (r = –0.626, p < 0.05), and to the overtraining score (r = –0.693, p < 0.05). After one year changes in hematocrit are a close reflect of the change in training volume (r = –0.877, p < 0.01) but are not correlated to fitness changes. Therefore in these athletes i) systemic hematocrit is close to its predicted ‘ideal value”, suggesting the accuracy of the prediction; ii) red cell aggregation is correlated to fat mass even in nonobese subjects; iii) hematocrit is lower than predicted by the model when markers of sympathetic tone (systolic blood pressure, overtraining score) are increased; iv) weekly training volume appears the main determinant of the reduction of hematocrit.
List of symbols
red blood cell
standard error on the mean
Introduction
Soccer is the sport were was initially described the “paradox of hematocrit”, i.e. the fact that hematocrit is negatively correlated with markers of aerobic capacity and decreases after training, while doping by erythropoietin or autotransfusion is based on the opposite paradigm, i.e. the higher hematocrit, the higher performance [3–5].
The classical concept of hematocrit/viscosity (h/η) ratio is likely in this context to be helpful for understanding this paradox. This parameter is theoretically proportional to the ability of blood to deliver oxygen, if one assumes that hematocrit reflects oxygen carrying by red cells while viscosity impairs blood flow. The ability of vessels to dilate in response to moderate increases in viscosity makes, obviously, this question more complex, but there is a large body of literature supporting the concept that h/η is a statistical determinant of O2 supply to tissues [21].
Our purpose was thus to propose a theoretical prediction of ‘optimal’ hematocrit and h/η in athletes on the basis of a widely admitted equation of blood viscosity, and to compare it to the actual value, and to see how the agreement between theoretical and actual values is modified by changes in training volume and performance.
Subjects and methods
Study subjects
A sample of elite soccer players (national level: 18–32 yr, weight 61–83 kg, body mass index 20.9–25.8 kg/m2). Characteristics of patients are shown on Table 1. Their mean systolic blood pressure 11.4±0.19 cmHg and their mean diastolic blood pressure was 7.25±0.14. Their eating habits collected with an autoquestionnaire and further calculated with a program as previously reported [12] were 2380±149 calories, with 40±4% carbohydrates, 51±4% lipids, 10.4±0.5 proteins, i.e. a daily protein intake evaluated at 0.77±0.05 g/kg/day. Subjects were seen at the same time of the season one year later. Maximal aerobic capacity (VO2max ACSM) was assessed from the submaximal steps according to the American College of Sports Medicine [8]. In addition to a complete clinical examination, the French questionnaire of early symptoms of overtraining developed by the French Society for Sports Medicine (SFMS) downloadable on the website of this society at < http://www.sfmes.org/images/sfmes/pdf/consensus_3.pdf>was employed.
Bioelectrical impedance measurements
Prior to the exercise-test, subjects’ body composition was assessed with bioimpedanceanalysis with a six terminal impedance plethismograph BIACORPUS RX 4000, (Healthnesslink, 8 avenue Jean-Jaurès 92130 Issy-les-Moulineaux, France) with data analysis with the software BodyComp 8.4. This device measures total resistance of the body to an alternative electric current of 50 kHz.
Hemorheological ex vivo measurements
Blood samples for hemorheological measurements (7 ml) were drawn with potassium EDTA as the anticoagulant in a vacuum tube (Vacutainer). Viscometric measurements were done at very high shear rate (1000 s–1) with a falling ball viscometer (MT 90 Medicatest, F-86280 Saint Benoit). The coefficient of variation of this method ranges between 0.6 and 0.8%. We measured with this device apparent viscosity of whole blood at native hematocrit, plasma viscosity, and blood viscosity at corrected hematocrit (45%) according to the equation of Quemada [15].
where φ is hematocrit, ηp is plasma viscosity, and k(γ) is a shear-dependent parameter quantifying the contribution of erythrocyte rheological properties to whole blood viscosity. At the high shear rate used here k(γ) is representative of red cell rigidity (i.e., the lower k(γ), the higher is erythrocyte deformability).
With this equation it is possible to standardize η for hematocrit 45% after calculating k:
This value of k is reintroduced in Equation (1) with φ set at 0.45. Dintenfass’s ‘Tk’ index [6] was also calculated as an index of red cell rigidity.
RBC aggregation was assessed with the Myrenne aggregometer [16] which gives two indices of RBC aggregation: “M” (aggregation during stasis after shearing at 600 s–1) and “M1” (facilitated aggregation at low shear rate after shearing at 600 s–1). Hematocrit was measured with microcentrifuge.
The curve of theoretical optimal h/η plotted vs hematocrit was reconstructed with Quemada’s equation presented above. The equation of h/η as a function of h was thus:
The hematocrit corresponding to the top of this curve was considered as the “theoretical optimal hematocrit”. The highest value of h/η (the top of the curve) was considered as the optimal h/η.
Values are presented as mean±standard error of the mean (SEM). Normality of samples was checked with the Kolmogorov-Smirnov test. After verification of normality, we used ANOVA followed by post-hoc Student’s t test for paired samples. Correlations were assessed with Pearson’s procedure (least square fitting). A value of p < 0.05 was considered as significant.
Results
Evolution of anthropometric and ergometric parameters
Since there was not a training study but a follow-up in stable conditions, there was no overall change in fitness or blood viscosity. After one year non significant tendency to a decrease in the mean training volume (hrs/wk) as assessed by questionnaire. At the beginning they exercised 11±1 hr/wk and one year later they exercised 9±1 hr/wk. This weekly volume of training decreased in 4 subjects and increased in 5.
Hemorheologic parameters
Figure 1 shows the curve of theoretical optimal h/η plotted vs hematocrit reconstructed with Quemada’s equation. This picture shows that on the average the actual values of hematocrit and h/η exhibit a nonsignificant tendency to be lower than predicted with the model.
As shown on Table 2, after one year non significant tendency to a decrease in hematocrit, that was 42.78±0.52% at the beginning of the study and 41.67±1.08% one year later.
Correlations
As shown on Fig. 2, after one year, changes in hematocrit are a close reflect of the change in training volume (r = –0.877, p < 0.01). By contrast (data not shown), the correlation between hematocrit and aerobic capacity is not found in this sample (r = 0.075 before the year of follow-up and r = 0.093 after the year of follow-up).
As shown on Fig. 3, RBC aggregation is correlated at baseline with fat mass (‘M1’ = 0.552, p < 0.02) and changes in aggregation are related to changes in fat mass (‘M’ = 0.652, p < 0.05; ‘M1’ = 0.647, p < 0.05).
Agreement between theoretical and actual values of hematocrit and h/η
As shown on Fig. 4, predicted and actual hematocrit are correlated but exhibit discrepancies (mean difference – 1% range [3.24 to 1.24]). Figure 5 shows that these discrepancies are correlated with two parameters: the value of predicted hematocrit (r = –0.912, p < 0.01), and the clinical score of overtraining (r = –0.693, p < 0.05). In addition discrepancies are also correlated to diastolic blood pressure (r = 0.626, p < 0.05) but neither systolic blood pressure (r = –0.572, NS), pulse pressure (r = –0.409, NS), nor mean blood pressure (r = 0.314, NS).
Discussion
This study shows in a sample of elite athletes that Quemada’s equation predicts an optimal value hematocrit which is very close to the actual one, suggesting that this calculation predicts rather accurately the ‘ideal’ hematocrit, and that this actual hematocrit is set very close to this ‘ideal’ value. However hematocrit is lower on the average than this model-predicted “ideal value” when markers of sympathetic tone (systolic blood pressure, overtraining score) are moderately increased. Weekly training volume rather than fitness appears the main determinant of training-induced reduction of hematocrit. Red cell aggregation is correlated to fat mass within a physiological range.
This study has the advantage of representing a follow-up in real life conditions, with some subjects undergoing an increase in their aerobic performance and some others undergoing a decrease.
Since the prediction is based on the h/η paradigm, the agreement between predicted and actual values supports to some extent this classical concept. Our data are consistent with the hypothesis that there is a homeostatic mechanism aiming at maintaining h/η close from its optimal value defined by the top of the bell-shaped curve of h/η plotted as a function of hematocrit. Two mechanisms of ‘viscoregulation’ have been until now evidenced: the modulation of hematocrit by a loop involving the hormone erythropoietin and triggered by plasma viscosity [17] and the decrease in red cell rigidity as a response to excessive erythrocytosis [20]. Both mechanisms prevent excess blood viscosity. In other terms, not surprisingly, hematocrit which is the most important determinant of blood viscosity is the pivotal parameter in viscoregulation, being both a stimulus and a regulated parameter.
Our modelling approach requires some comments. The equation of Quemada [15] is a prediction of blood viscosity from hematocrit, plasma viscosity and a red cell rheology factor ‘k’ that represents red cell rigidity when the measurements are performed, as done in this study, at high shear rate.
It is well known that hematocrit is not the same in various territories of circulation, and can undergo huge variations according to physiological conditions. For example Vicaut reported in studies on rat myocardium a sub-epicardial hematocrit of 33.1±0.6% while systemic hematocrit was 39.9±2.2% [19]. Okazawa [14] found in human brain 38.3±3.5% while systemic hematocrit was 43.5±4.1%. Lipowsky [11] in 10–15 μm arterioles and venules of cat mesentery reported a microvascular hematocrit of 8.2±4.3% while systemic hematocrit was 35.5±5.6%. Klitzman [10] in 15 μm arterioles of hamster cremaster muscle found a local hematocrit of 13.9±1.2% while systemic hematocrit was 49.4±1.1%. Despite some controversies and the differences among methodologies it is clear that in local circulations hematocrit is lower, and sometimes much lower, than the systemic one, and can undergo very important changes in physiological conditions, such as exercise [9].
Our data deal with systemic hematocrit which is thus higher than these local hematocrits and show that the predicted ‘optimal’ value seems to represent an upper limit of systemic hematocrit in physiological conditions.
We should also point out that our measurements are relevant to venous blood in which the shear rate is lower than that employed in our viscometric measurements. Actually it seems that the concept of h/η is mostly relevant for measurements at high shear rate. With lower shear rates the curve of h/η plotted as a function of hematocrit loses its bell-shaped aspect and it is no longer possible to detect its top [13], probably due to the more complex effect of red cell aggregation. Therefore what we report here is probably not relevant to venous physiology but rather to the arterial bed in which hematocrit is close to the venous one. Most of the time arterial hematocrit is said to be lower than the venous one, as found for example by Yang et al. [22] who reports an average difference of 3.1%. However this difference is not always observed and other studies [2, 18] actually conclude that there is no difference on the average between arterial and venous haematocrit, whose distribution display a large overlap. In a recent study performed in India venous and arterial hematocrit were remarkably correlated with r = 0.99 and a mean difference of 0.3% [1]. Therefore, despite erratic differences of 0.3 to 3% venous hematocrit seems to reflect quite closely on the average arterial hematocrit, which is more likely to be the regulated parameter.
Despite the impressive agreement found between the theoretical prediction of ideal hematocrit and actual values of venous hematocrit, it is interesting to scrutinize the discrepancies between these two values. In this study we evidence three factors correlated with this discrepancy: the predicted hematocrit itself, the clinical score of overtraining, and blood pressure. A correlation between this discrepancy and the predicted optimal hematocrit itself may indicate limits of the model, i.e. when it predicts too high ideal values of hematocrit the actual value is no longer set close of it but at a lower level. We hypothesized that the inclusion of plasma viscosity in our model was not accurate and that a model with fixed plasma viscosity would overcome this effect, but this was not actually the case and a new calculation with fixed plasma viscosity (data not shown) gave a very poor prediction of optimal hematocrit, indicating that the full model including all viscosity factors is better.
Correlations between the score of overtraining and diastolic blood pressure may indicate that changes in the sympathetic tone modifying vascular capacity may interfere with this regulation and increase the gap between ‘optimal’ and ‘actual’ hematocrit. This remains to be more precisely studied.
Finally, this study also suggests that hematocrit in soccer players is under the influence of factors other than the classical viscoregulatory homeostatic loops. The correlation we find with the training volume (while hematocrit in this study is not correlated to aerobic capacity) is in line with the notion that repeated bouts of exercise result in a ‘hemodilution’ prolonged over a couple of days. Therefore the negative correlations repeatedly reported between fitness and hematocrit may be largely explained by the covariable ‘exercise volume’ which is closely related to both of them. This specific issue will require additional studies.
Conclusions
In conclusion, this study shows that in a team of professional soccer players venous hematocrit is set close to its theoretical ‘optimal’ value as predicted with a classical equation of viscosity. The fair agreement between the theoretical prediction and the actual value supports the theory of ‘optimal hematocrit’ and its determination with the bell shaped curve of h/η as a function of hematocrit. However the discrepancies between the theoretical and the actual hematocrit remain to be further analyzed.
Besides, we confirm that RBC aggregation is related to fat mass also in lean subjects and not only in obese, indicating that fat mass is a major regulator of RBC aggregation.
Finally, the strong relationship between changes in training volume and changes in hematocrit is also appealing, indicating that low hematocrit in athletes may be rather a marker of training volume than fitness.
Footnotes
Acknowledgments
The authors affirm that they have no disclosure of interest concerning the issue investigated in this paper. This study was presented at the 15th International Congress of Biorheology and 8th International Conference on Clinical Hemorheology. May 24-28, 2015, Seoul, Korea.
