Warning: fopen(/home/virtual/pediatrics/journal/upload/ip_log/ip_log_2026-11.txt) [function.fopen]: failed to open stream: Permission denied in /home/virtual/pediatrics/journal/ip_info/view_data.php on line 93

Warning: fwrite(): supplied argument is not a valid stream resource in /home/virtual/pediatrics/journal/ip_info/view_data.php on line 94
Hematological alterations in adolescent girls with restrictive anorexia nervosa and secondary amenorrhea at the time of diagnosis

Hematological alterations in adolescent girls with restrictive anorexia nervosa and secondary amenorrhea at the time of diagnosis

Article information

Clin Exp Pediatr. 2026;69(11):891-901
Publication date (electronic) : 2026 October 19
doi : https://doi.org/10.3345/cep.2026.00808
1Institute of Molecular Biomedicine, Faculty of Medicine, Comenius University, Bratislava, Slovakia
2Department of Pediatrics, Faculty of Medicine, National Institute of Children's Diseases, Comenius University, Bratislava, Slovakia
Corresponding author: Katarína Šebeková, MD, DSc. Institute of Molecular Biomedicine, Faculty of Medicine, Comenius University, Sasinkova 4, 811 08 Bratislava, Slovakia Email: sebekova3@uniba.sk, kata.sebekova@gmail.com
Received 2026 April 4; Revised 2026 August 11; Accepted 2026 August 13.

Abstract

Background

Anorexia nervosa (AN) is associated with chronic energy deficiency and endocrine adaptations affecting multiple organ systems involving bone metabolism and hematological homeostasis. However, early alterations in peripheral blood counts of adolescents with a clinically homogeneous AN phenotype remain insufficiently characterized.

Purpose

To examine the differences in blood counts and associated bone and endocrine parameters between adolescent girls with newly diagnosed restrictive AN and secondary amenorrhea and healthy age-matched controls.

Methods

The study included 59 patients and 20 controls. Hematological indices, anthropometric characteristics, body composition, bone mineral status, bone turnover, and endocrine status markers were compared, and a multivariate analysis was used to identify the predictors of blood counts.

Results

Compared to controls, patients had significantly lower leukocyte, erythrocyte, and platelet counts (all P<0.01), and 29% had cytopenia in at least one parameter. Bone formation markers (osteocalcin, amino-terminal propeptide of human procollagen type I, and alkaline phosphatase) and key endocrine parameters (insulin-like growth factor-1, estradiol, and free thyroxine) were suppressed. These alterations were also observed in patients with a normal bone mineral density and adequate vitamin B12 and folate levels. A multivariate analysis identified endocrine- and bone-related predictors of blood counts in controls but no consistent pattern of predictors in patients.

Conclusion

Adolescent girls with restrictive AN and secondary amenorrhea demonstrated hematological alterations at the time of diagnosis despite the absence of a measurable skeletal impairment or micronutrient deficiency These findings suggest that a chronic energy deficiency and endocrine adaptation are associated with early disturbances in hematological homeostasis. Peripheral blood count may represent an accessible clinical marker of systemic disturbances.

Key message

Question: Do adolescent girls with restrictive anorexia nervosa and secondary amenorrhea have low blood counts at diagnosis?

Finding: Patients, even those with bone mineral density z scores similar to those of controls and adequate micronutrient levels, exhibited lower blood counts and suppressed bone formation and endocrine parameters compared with healthy age-matched controls.

Meaning: Hematological alterations may occur early in the disease course and complement routine clinical assessments.

Graphic abstract. AN, anorexia nervosa; IGF-1, insulin-like growth factor-1; BMD, bone mineral density; PINP, aminoterminal propeptide of human procollagen type I; fT4, free thyroxine; OPLS, orthogonal partial least squares.

Introduction

Anorexia nervosa (AN) is a life-threatening eating disorder characterized by distorted body image, intense fear of weight gain, and self-imposed starvation, predominantly affecting adolescent girls and young women [1]. Chronic severe energy deficiency results in profound metabolic and endocrine adaptations that influence multiple organ systems. Because AN frequently develops during adolescence—a critical period for bone mass acquisition—early deterioration of bone mineral status is frequently observed [2]. Malnutrition in AN is accompanied by complex neuroendocrine dysregulation involving the hypothalamic- pituitary-gonadal, hypothalamic-pituitary-thyroid, and hypothalamic-pituitary-somatotropic axes [3,4].

These endocrine adaptations reflect a coordinated response to chronic energy deprivation and contribute to impaired body composition and bone formation. They also affect the bone marrow microenvironment, which plays a central role in both hematopoiesis and skeletal homeostasis.

Bone marrow is a metabolically active organ that integrates hematopoietic, adipogenic, and osteogenic lineages and is particularly sensitive to chronic energy deficiency. Under conditions of prolonged energy deprivation, mesenchymal stem cell differentiation may shift toward adipogenesis at the expense of osteoblastogenesis. Consistent with this concept, despite marked depletion of peripheral fat stores, patients with AN demonstrate paradoxical increases in bone marrow adiposity, meaning that fat accumulates in the bone marrow even as it is lost elsewhere in the body [5,6]. Accumulation of fat in bone marrow may disrupt the hematopoietic niche—the specialized microenvironment supporting blood cell development—leading to abnormalities in peripheral blood counts, including leukopenia, erythropenia, and thrombocytopenia [7].

Previous studies in adolescents with AN have focused primarily on associations with anthropometric severity or disease duration. Investigations often included heterogeneous populations with different menstrual statuses (premenarchal girls, those with primary and secondary amenorrhea), mixed-sex cohorts, patients with different AN subtypes, or did not report the type of AN [8-11]. These factors may introduce potential confounders to the interpretation of hematologic findings.

Thus, it remains unclear whether early hematologic alterations in adolescents with AN occur independently of skeletal impairment or micronutrient deficiency. Moreover, how peripheral blood counts relate to bone metabolism and endocrine adaptations to chronic energy deficiency has not been sufficiently characterized in adolescents in the early stages of the disease.

To address this gap, we evaluated hematologic parameters in adolescent girls at the time of diagnosis of restrictive AN and secondary amenorrhea—a relatively homogeneous phenotype characterized by marked energy deficit and endocrine suppression. We examined the associations of hematologic parameters with anthropometric measures, bone mineral status, bone turnover indices, and endocrine markers. Because these independent variables are biologically interconnected and highly collinear, univariate regression approaches may not capture coordinated physiological patterns in the dataset's relatively high dimensionality. Thus, we applied orthogonal partial least squares (OPLS) regression to explore whether physiological determinants of blood counts differ between healthy adolescents and patients with AN. We hypothesized that patients would exhibit quantitative alterations in peripheral blood counts compared with healthy controls, and that physiological predictors of blood counts observed in healthy adolescents would be less clearly expressed in AN.

Methods

1. Study design

This is a retrospective analysis of a cross-sectional case-control study that included girls with restrictive AN and healthy controls, with the primary outcome being a comparison of bone mineral status and endocrine parameters. This secondary analysis aimed to compare blood counts across groups to investigate the impact of energy deficiency on bone marrow-derived lineages.

2. Study population

We included 59 girls at first diagnosis of restrictive AN and secondary amenorrhea. Restrictive AN was defined according to the Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM-5) criteria [1]. Thus, the patients met the general DSM-5 criteria for AN, which include: (1) significant restriction of energy intake leading to significantly low body weight; (2) intense fear of gaining weight or becoming fat; and (3) disturbance in self-perceived weight or shape. A key feature of a restrictive subtype of AN is the absence of regular bingeing or purging over the past 3 months. Thus, weight reduction is achieved by restricting calorie intake via dieting, fasting, or excessive exercise. Secondary amenorrhea was defined as an absence of regular menstruation for at least 3 months following previously normal cycles [12]. Girls were consecutively recruited at the Department of Pediatrics at the National Institute of Children's Diseases (NICD) in Bratislava between October 2016 and February 2020 among patients hospitalized for the differential diagnosis of severe weight loss. Based on medical history, the onset of the disorder was defined as the time point at which the patient, at her highest attained body weight, initiated intentional weight loss. The duration of the disease was calculated from the onset till hospitalization. AN was diagnosed based on parental information, patient observation, clinical examination, interviews, and consultative evaluations by a psychiatrist or psychologist to exclude other psychiatric disorders.

Twenty healthy age-matched girls recruited via primary care pediatricians practicing within the NICD premises served as controls. All controls and 31 patients with secondary amenorrhea had been included in our former study on bone mineral status in amenorrheic girls with AN [13]. An additional 28 patients with secondary amenorrhea were initially enrolled in a study of oxidative status [14]. The present analysis addresses a distinct research question focused on hematologic parameters and their associations with endocrine and bone-related markers. General exclusion criteria included neurological, gastrointestinal, endocrine, and skeletal disorders; acute or chronic inflammatory diseases; pregnancy; lactation; and current hormonal contraception. The study was conducted in accordance with the Declaration of Helsinki, following approval of the protocol by the Ethics Committee of the NICD, Bratislava, Slovakia, on August 10, 2016. In accordance with Slovak legislation, written informed consent was obtained from the girls' legal guardians, in addition to the participants' verbal assent.

3. Measurements

Clinical and laboratory data were prospectively collected during the initial diagnostic workup, prior to initiation of the re-feeding regimen, in accordance with a predefined institutional protocol. The blood counts were obtained from the same blood draw used for standard blood chemistry, endocrine, and bone metabolism markers.

4. Anthropometry

Height was measured to the nearest 0.1 cm using a wall-mounted stadiometer, and body weight to the nearest 0.1 kg using electronic scales (both: Seca, Germany). Body mass index (BMI) was calculated. Age- and sex-specific reference data from Slovak children and adolescents [15] were used to calculate standard deviation scores.

5. Bone mineral density and body composition

Bone mineral density (BMD, g/cm2), bone mineral content (BMC, g) of the lumbar spine (L1–4), left hip, left femoral neck, lumbar spine trabecular bone score (LSTBS), and body fat percentage were measured using a dual-energy x-ray absorptiometry scanner (Horizon QDR, Hologic Inc., USA), calibrated daily using a manufacturer-supplied standard block. The software (QDR APEXTM v. 5.6.0.5), integrated into the Hologic QDR APEX operating system, was used to compare BMD results against a manufacturer-provided pediatric population reference database and to report them as z scores.

6. Blood sampling and analyses

Blood was collected from the cubital vein in the morning, after at least 8 hours of fasting, into K3EDTA, Liheparin, and gel & clot activator tubes (Vacutainer Plastic Tube, Becton Dickinson, Czech Republic). Blood counts were determined using the Sysmex XN-1000 Hematology Analyzer (Sysmex, Japan). Total erythrocyte, leukocyte, or platelet counts below the lower limit of the appropriate reference ranges of the NICD’s laboratory (Supplementary Table 1) were termed erythropenia, leukopenia, and thrombocytopenia, respectively.

Routine blood chemistry was performed using standard methods on the Cobas c501 analyzer (Roche Diagnostics, Germany). The Vitros 5600 analyzer (Ortho Clinical Diagnostics, USA) was used to determine plasma concentrations of 25-OH vitamin D, vitamin B12, and folate. Vitamin D levels were obtained only from 18 controls and 58 patients. Folate and vitamin B12 concentrations were analyzed only in patients. Adequacy of vitamin levels was assessed according to the NICD’s reference ranges (Supplementary Table 1).

Parathyroid hormone, thyroid-stimulating hormone (TSH), free thyroxine (fT4), growth hormone (GH), cortisol, estradiol, and insulin-like growth factor-1 (IGF-1) were measured using chemiluminescence immunoassays on the Architect i20000SR analyzer (Abbott Laboratories, USA). Osteocalcin, amino-terminal propeptide of human procollagen type I (PINP), and carboxy-terminal telopeptide of type I collagen (CTX) were determined employing the electrochemiluminescence method (Cobas e411, Roche Diagnostics).

7. Statistics

The sample size was based on detecting a 15% difference in leukocyte counts at α=0.05 with 80% power and a 1:3 enrollment ratio, yielding groups of 19 controls and 57 patients.

Data distribution was assessed using the Kolmogorov- Smirnov test. Normally distributed data were compared using the 2-sided independent Student t test after checking for equality of variances with Levene test, and are presented as the mean±standard deviation (SD). When equal variances could not be assumed, we used Welch t test. Skewed data were compared using the Mann- Whitney U test and reported as medians (interquartile ranges). Categorical data were expressed as frequency and percentage. Significance was set at P<0.05. In addition to P values, which address only statistical significance, Cohen d effect sizes were calculated to quantify the magnitude of group differences between patients and controls, thereby facilitating interpretation of the potential clinical relevance of statistically significant findings, especially given unequal sample sizes. For normally distributed data, Cohen d is given with a 95% confidence interval. Effect sizes were interpreted according to conventional thresholds as ≥ |0.2| - small, ≥ |0.5| - medium, and ≥ |0.8| - large. Analyses were performed using IBM SPSS Statistics ver. 22.0 (IBM Corp., USA) and Psychometrica [16], and figures were generated using GraphPad Prism v. 11.0.2 (GraphPad Software, Inc., USA).

Multivariate analyses were performed using SIMCA v. 18.2 software (Sartorius Stedim Data Analytics AB, Sweden). Variables with high skewness and low min/max ratio were log-transformed, and all data were mean-centered. Principal component analysis was used to explore the sample distribution without classification. Hotelling T2 test was applied to identify potential outliers. We applied an OPLS model for multivariate regression because it is well-suited to high-dimensional, intercorrelated data and can enhance interpretability. The OPLS analysis was intended as a complementary multivariate approach to support our central hypothesis that chronic energy deficiency disrupts the normal physiological relationships linking nutritional, endocrine, skeletal, and hematopoietic regulation. Independent variables entered into the model included femoral neck z score, body fat percentage, IGF-1, CTX, PINP, alkaline phosphatase (ALKP), fT4, estradiol, and, in patients, duration of amenorrhea or disease duration. Variables with Variable Importance in the Projection (VIP) scores ≥1 were considered relevant contributors. Measures of fit (R²) and predictivity (Q²) ≥0.5 were considered acceptable. Significant models were validated using a permutation test (n=100). The criterion for validity was that the Q2 regression line intersected the vertical axis at or below zero. The detailed descriptions of the model evaluation and outcomes are provided in the Supplementary material, as they primarily support the multivariate analysis, whereas the principal biological findings are summarized in the main manuscript.

Results

Cohort characteristics are given in Table 1. In 58 out of 59 patients, disease duration ranged from 3 to 30 months; in the remaining girl, it was 50 months. Controls had blood counts within the reference ranges. Patients had lower total leukocyte counts (-16%; Cohen d=-0.8; 95% confidence interval [CI], -1.3 to -0.3) (Fig. 1A) than controls, and 10 (17%) had leukopenia. Across all leukocyte lineages, significance was observed only in the monocyte percentage, which was lower in patients (Table 1). They also showed lower erythrocyte counts (-6%; Cohen d=-0.8; 95% CI, -1.4 to -0.3) (Fig. 1B). Nine (15%) had erythropenia (4 of whom concurrently had leukopenia). Patients had lower hematocrit (-5%), while hemoglobin concentration and mean cell volume (MCV) did not differ significantly between the groups. Girls with AN exhibited higher mean cell hemoglobin and red cell distribution width, but the differences were within the reference range (Table 1). Platelet counts (-14%; Cohen d=-0.8; 95% CI, -1.3 to -0.3) (Fig. 1C) and plateletcrit (-19%, Table 1) were also lower in patients than in controls. Two patients (3%) had thrombocytopenia. Mean platelet volume was similar in both groups (Table 1). Girls with AN had adequate folate and vitamin B12 levels (Table 1), although one exhibited a vitamin B12 level below the reference range.

Anthropometrics and blood counts of healthy controls versus patients with anorexia nervosa

Fig. 1.

Blood counts, markers of bone status, and estradiol and IGF-1 concentrations of patients with AN (n=59) versus controls (n=20). Normally distributed data (panels A, B, C, D, F, and H) were compared using the 2-sided independent Student t test after checking for equality of variances with Levene test and are shown as mean and standard deviation. Skewed data (panels E and G) were compared using the Mann-Whitney U test, and medians and interquartile ranges are indicated. Statistical significance was set at P<0.05. IGF-1, insulin-like growth factor-1; AN, anorexia nervosa; BMD, bone mineral density; PINP, amino-terminal propeptide of human procollagen type I.

Hip BMD z scores were lower in patients than in controls (Cohen d=-0.6; 95% CI, -1.2 to -0.1) (Fig. 1D), whereas those of lumbar spine, femoral neck, and LS TBS were lower without reaching significance (Table 1). Six patients (10%) displayed z scores of ≤-2 SD (considered "below the expected range for age") [17] at least at one measured site.

Plasma PINP (Cohen d=-1.8; Fig. 1E) and osteocalcin levels, and ALKP activity were lower in patients. CTX concentrations did not differ significantly between the groups (Table 1). Patients had higher 25-OH vitamin D levels than controls (Cohen d=0.8; 95% CI, 0.3–1.3) (Fig. 1F); 56% of controls and 36% of patients had values below the reference range. Data on anthropometric measurements, BMD, BMC, and plasma mineral concentrations are provided in Supplementary Table 2.

Patients had lower estradiol (Cohen d=-1.2; Fig. 1G) and IGF-1 (Cohen d=-1.4; 95% CI, -1.9 to -0.8) (Fig. 1H) levels than controls. Concentrations of cortisol and TSH were similar between groups, whereas fT4 was lower and GH was higher in patients than in controls (Table 1).

To determine if reduced blood counts were independent of bone density, we analyzed a subset of patients with BMD z scores ≥-1.4 SD at all sites, matching the lower range of our control group. Patients without measurable deterioration of bone mineral status, as shown in Table 2 and Fig. 2D, nevertheless exhibited significantly lower leukocyte (-15%; Cohen d=-0.7; 95% CI, -1.2 to -0.2; Fig. 2A), erythrocyte (-6%; Cohen d=-0.9; 95% CI, -1.4 to -0.3) (Fig. 2B), and platelet (-14%; Cohen d=-0.8; 95% CI, -1.3 to -0.2) (Fig. 2C) counts than controls. Seven of 45 patients (16%) had erythropenia; 5 (11%) had leukopenia (2 of whom had it combined with erythropenia); and 1 (2%) had thrombocytopenia. Compared with controls, these patients exhibited alterations in monocyte percentage, hematocrit, MCV, mean cell hemoglobin (MCH), red cell distribution width (RDW), plateletcrit, and mean platelet volume (MPV) similar to those observed in the whole group, all associated with medium-to-large effect sizes. Concurrently, suppression of bone formation markers and dysregulation of key hormonal axes were also evident (Table 2; Fig. 2E–G). Data on anthropometrics, BMD, BMC, and plasma mineral concentrations are provided in Supplementary Table 3.

Anthropometrics and blood counts of healthy controls versus patients with anorexia nervosa presenting with BMD z scores ≥ -1.4 SD

Fig. 2.

Blood counts, markers of bone status, and estradiol and IGF-1 concentrations in patients with anorexia nervosa presenting with BMD z scores ≥ -1.4 SD (n=45) and healthy controls (n=20). Normally distributed data (panels A, B, C, D, F, and H) were compared using a 2-sided independent Student t test after a check for equality of variances with Levene test and are given as mean and SD. Skewed data (panels E and G) were compared using the Mann-Whitney U test and are shown as median and interquartile range. Statistical significance was set at P<0.05. IGF-1, insulin-like growth factor-1; AN, anorexia nervosa; BMD, bone mineral density; PINP, amino-terminal propeptide of human procollagen type I; SD, standard deviation.

Descriptive statistics for the OPLS models are provided in Supp. Table 4. In healthy controls, the leukocyte model explained 64% of the outcome variance with marginal predictivity (Q²=0.50). FT4, femoral neck z score, and estradiol had VIP values ≥1. The model passed the permutation test. Loading, VIP, and permutation plots are presented in Supplementary Fig. 1A–C. The erythrocyte and platelet models explained 52% and 43% of outcome variance, respectively, but showed limited (Q²=0.24), or poor (Q²=0.03) cross-validated predictivity. In patients with anorexia nervosa, all three models showed poor explanatory performance (25%, 13%, and 11%, respectively) and low or even negative Q² values (0.04, −0.01, and −0.02, respectively). Except for leukocytes in controls, the remaining models did not meet the stated criteria for acceptable predictive performance.

Discussion

The most relevant finding of this study was that adolescent girls with restrictive AN and secondary amenorrhea, whose BMD z scores were similar to those of controls and whose vitamin status was adequate, had, at the time of diagnosis, significantly lower leukocyte, erythrocyte, and platelet counts than healthy controls. Our findings suggested that quantitative hematologic alterations in AN coincide with the disruption of physiological relationships linking nutritional and endocrine status, bone metabolism, and hematopoiesis. Peripheral blood counts may serve as clinical indicators of systemic metabolic adaptation to chronic energy deficiency in adolescent AN.

Leukocyte counts were lower in the AN group, with a medium-to-large effect size. The relative proportions of leukocyte lineages suggest a generalized reduction in production (probably due to chronic undernutrition), although a mild, selective suppression of the monocyte lineage may reflect an impaired innate immune response. Lower erythrocyte counts and hematocrit in girls with AN, with hemoglobin concentrations similar to those of the controls, suggest mild or compensated anemia-related changes, which may be partially explained by the relatively short duration and moderate severity of the disease. Higher MCH and RDW, with a trend toward higher MCV, may indicate early macrocytic adaptation or increased heterogeneity in erythrocyte size, potentially related to nutritional deficiencies. Since reticulocyte counts were not assessed in the standard differential blood count, we can only speculate that a higher RDW may partially reflect impaired erythrocyte maturation. Mean vitamin B12 and folate levels in the patients were not suggestive of severe micronutrient deficiency, implying that the observed alterations may primarily reflect caloric restriction and marrow suppression rather than overt micronutrient depletion. Lower total circulating platelet mass in patients, together with similar MPV across groups, may indicate that platelet production remains qualitatively intact, while overall megakaryocyte activity is suppressed. Thus, our findings likely reflect adaptive hematopoietic and immunologic alterations associated with chronic energy deficiency.

Reduced blood counts are well recognized in AN, and not all 3 hematopoietic lineages are affected simultaneously [18]. In other studies, girls with AN displayed lower leukocyte and platelet counts [9] or lower leukocyte and erythrocyte counts [8] than their healthy peers. In adolescents with AN, a decline in peripheral blood cell counts is generally mild, and the reduction is only partly attributable to cytopenias [8,19-21]. Overt cytopenias are reported particularly in long-lasting (defined as disease duration ranging from 3 to 10 years, with a mode of 7 years), enduring (defined as multiple unsuccessful treatment attempts), and severe disease (BMI <16 kg/m2 in those aged ≥15 years; there is no consensus on percentile for younger patients) [1,22,23], predominantly in adults, with prevalence increasing as malnutrition progresses. In these cases, pancytopenia and morphological changes in the bone marrow, including hypocellularity and gelatinous marrow transformation, may occur [4,24-27]. This suggests that impaired hematopoiesis in AN may be independent of overt skeletal impairment, and not primarily driven by micronutrient deficiency, but rather reflects systemic metabolic adaptations to chronic energy deprivation. Our findings extend these observations by demonstrating that in a relatively homogeneous adolescent cohort, hematologic alterations may be detectable at diagnosis, or early in the disease course, even in patients with preserved BMD and adequate folate and vitamin B12 levels.

Despite severe peripheral fat loss, patients with AN paradoxically exhibit increased bone marrow adiposity [5,28]. Metabolically active bone marrow adipose tissue (BMAT) influences hematopoietic activity and bone remodeling by regulating interactions among osteoblasts, adipocytes, and hematopoietic stem cells [6,29,30]. During periods of caloric restriction, BMAT expands and secretes paracrine and endocrine factors that redirect mesenchymal stem cell differentiation in the bone marrow toward adipocytes rather than osteoblasts. We did not evaluate the bone marrow status. Young females with AN and a history of fractures had higher BMAT but similar BMD or bone microarchitecture compared with those without [31]. Moreover, about 39% of adolescents and young adults with AN exhibit bone marrow hypoplasia or aplasia, and 50% display partial, focal, or complete gelatinous degeneration of the bone marrow [32]. Importantly, peripheral blood counts in AN may underestimate the extent of bone marrow involvement [25,32].

The coexistence of marked energy deficiency with higher vitamin D concentrations and preserved vitamin B12, folate, and mineral status in the patients at the time of diagnosis probably reflects sustained self-supplementation rather than adequate nutritional intake from food sources alone. Regular intake of vitamin supplements and minerals is 10 times more frequent in adolescent girls with AN than in controls [8], probably explaining the lower prevalence of vitamin D deficiency (2% in AN vs. 24% in controls) [33].

Patients with AN exhibited marked endocrine alteration, including lower concentrations of IGF-1, estradiol, fT4, and higher GH, together with suppression of bone formation markers. These hormonal abnormalities are well-recognized adaptive responses to chronic energy deficiency and were accompanied by hematologic alterations that were present even in patients with preserved bone mineral status. Endocrine and bone-derived pathways have been implicated in the regulation of the hematopoietic microenvironment [8,34]. In contrast to healthy controls, multivariate analyses in patients with AN did not identify stable associations among hematologic, endocrine, and bone-related variables. Thus, our results do not support a consistent relationship between endocrine parameters and peripheral blood counts in adolescents with AN. Given the limited predictive performance of the OPLS model and the lack of its validation in patients, the exploratory associations identified in this study should be interpreted with caution. The observed hematologic abnormalities should be considered primarily as concurrent manifestations of malnutrition rather than as direct consequences of specific endocrine alterations.

Although malnutrition-related bone marrow and hormonal abnormalities are largely reversible with nutritional rehabilitation [34], skeletal deficits acquired during adolescence may persist, increasing the risk of premature osteoporosis and fractures [35].

This study has limitations. Its cross-sectional design does not allow causal inference regarding relationships among nutritional status, endocrine alterations, bone metabolism, and hematologic findings. The relatively small control group may have limited power to detect subtle associations; however, controls were carefully agematched and evaluated using the same protocol, enabling reliable comparisons, and effect sizes indicated medium-to-large effects. Our study also lacks a direct evaluation of bone marrow composition. In addition, we did not assess dietary intake, dietary supplement use, physical activity, iron status, or several other factors that may affect hematopoiesis and endocrine-bone interactions. Despite these limitations, the study provides novel insights into hematologic alterations in adolescent girls with newly diagnosed AN of a relatively homogeneous phenotype (restrictive AN with secondary amenorrhea). The well-defined cohort enabled evaluation of early disease-related changes while minimizing potential confounders (such as age, sex, AN subtype, reproductive status), and demonstrated that hematologic abnormalities may already be present even in patients without overt bone mineral impairment.

In conclusion, girls with restrictive AN and secondary amenorrhea showed lower peripheral blood cell counts alongside endocrine suppression and altered bone turnover, despite preserved BMD and the absence of overt micronutrient deficiency. While hematologic abnormalities in AN have been previously described, our findings extend existing evidence by showing that these alterations are also present in a clinically homogeneous subgroup of adolescents with relatively preserved bone status. These findings support the concept that hematologic alterations may reflect early systemic metabolic adaptation in AN, even before measurable bone loss becomes apparent. Future longitudinal studies, ideally including direct assessment of bone marrow characteristics, are needed to clarify underlying mechanisms and potential prognostic implications.

Supplementary materials

Supplementary Tables 1-4 and Supplementary Fig. 1 are available at https://doi.org/10.3345/cep.2026.00808.

Supplementary Table 1.

Age-dependent reference values in females - hematology and blood chemistry

cep-2026-00808-Supplementary-Table-1.pdf
Supplementary Table 2.

Anthropometric characteristics, bone mineral density, bone mineral content, and plasma minerals in healthy controls and girls with anorexia nervosa

cep-2026-00808-Supplementary-Table-2.pdf
Supplementary Table 3.

Anthropometric characteristics, BMD, BMC, and plasma minerals in healthy controls and girls with anorexia nervosa presenting with BMD z scores ≥ -1.4 in all measured locations

cep-2026-00808-Supplementary-Table-3.pdf
Supplementary Table 4.

The multivariate regression of selected variables on leukocytes, erythrocytes, and platelet counts in healthy controls and patients with anorexia using the orthogonal partial least squares model

cep-2026-00808-Supplementary-Table-4.pdf
Supplementary Fig. 1.

The multivariate regression model by OPLS of independent variables on leukocyte counts in healthy controls. (A) Green dots, independent variables; blue dot, dependent variable; IGF-1, insulin-like growth factor-1; FatPerc, body fat percentage (by DEXA); CTX, carboxy-terminal telopeptide of type I collagen; PINP, amino-terminal propeptide of human procollagen type I; ALP, alkaline phosphatase; E2, estradiol; fT4, free thyroxine; FemNeckscore, femoral neck BMD z score; Leu, leukocytes. (B) The VIP plot, variables with Variable Importance in the Projection (VIP) values ≥ 1.00 were considered as important predictors of leukocyte counts. OPLS, orthogonal projections to latent structures; CTRL, controls. (C) Permutation test plot validating the OPLS model. The permutation test results (Q2=[0.0; -0.592]) confirm the model’s accuracy.

cep-2026-00808-Supplementary-Fig-1.pdf

Notes

Conflicts of interest

No potential conflict of interest relevant to this article was reported.

Funding

The Ministry of Health of the Slovak Republic, No. 2018/36-LFUK-10; The Grant Agency of the Ministry of Education, Science, Research and Sport of the Slovak Republic (VEGA), No. 1/0613/17.

Author contribution

Conceptualization: KS, LP; Data curation: KS, AGK, LT; Formal analysis: KS, AGK, LT, LP; Funding acquisition: LP; Methodology: KS, AGK, LT; Project administration: LP; Writing - original draft: KS; Writing - review & editing: KS, AGK, LT, LP

References

1. American Psychiatric Association. Diagnostic and statistical manual of mental disorders (DSM-5) Washington, DC: American Psychiatric Publishing; 2013.
2. Legroux I, Cortet B. Factors influencing bone loss in anorexia nervosa: assessment and therapeutic options. RMD Open 2019;5e001009.
3. Schorr M, Miller KK. The endocrine manifestations of anorexia nervosa: mechanisms and management. Nat Rev Endocrinol 2017;13:174–86.
4. De Filippo E, Marra M, Alfinito F, Di Guglielmo ML, Majorano P, Cerciello G, et al. Hematological complications in anorexia nervosa. Eur J Clin Nutr 2016;70:1305–8.
5. Fazeli PK, Klibanski A. The paradox of marrow adipose tissue in anorexia nervosa. Bone 2019;118:47–52.
6. Liu H, Liu L, Rosen CJ. Bone marrow adipocytes as novel regulators of metabolic homeostasis: clinical consequences of bone marrow adiposity. Curr Obes Rep 2025;14:9.
7. Naveiras O, Nardi V, Wenzel PL, Hauschka PV, Fahey F, Daley GQ. Bone-marrow adipocytes as negative regulators of the haematopoietic microenvironment. Nature 2009;460:259–63.
8. Misra M, Aggarwal A, Miller KK, Almazan C, Worley M, Soyka LA, et al. Effects of anorexia nervosa on clinical, hematologic, biochemical, and bone density parameters in community-dwelling adolescent girls. Pediatrics 2004;114:1574–83.
9. Lambert M, Hubert C, Depresseux G, Vande Berg B, Thissen JP, Nagant de Deuxchaisnes C, et al. Hematological changes in anorexia nervosa are correlated with total body fat mass depletion. Int J Eat Disord 1997;21:329–34.
10. Ünver H, Gökçe Ceylan B, Erdoğdu Yıldırım AB, Perdahlı Fiş N. Serum peripheral markers for inflammation in adolescents with anorexia nervosa. Int J Psychiatry Clin Pract 2024;28:68–72.
11. Oświęcimska J, Malczyk Ż, Szymlak A, Mikołajczak A, Ziora K, Zamlynski J, et al. Changes in platelet count and size indices in adolescent patients with anorexia nervosa. Clin Appl Thromb Hemost 2017;23:562–6.
12. Practice Committee of American Society for Reproductive Medicine. Current evaluation of amenorrhea. Fertil Steril 2008;90(5 Suppl):S219–25.
13. Čagalová A, Tichá Ľ, Gaál Kovalčíková A, Šebeková K, Podracká Ľ. Bone mineral density and oxidative stress in adolescent girls with anorexia nervosa. Eur J Pediatr 2022;181:311–21.
14. Kovalčíková AG, Tichá Ľ, Šebeková K, Celec P, Čagalová A, Sogutlu F, et al. Oxidative status in plasma, urine and saliva of girls with anorexia nervosa and healthy controls: a cross-sectional study. J Eat Disord 2021;9:54.
15. Regecová V, Hamade J, Janechová H, Ševčíková Ľ. Comparison of Slovak reference values for anthropometric parameters in children and adolescents with international growth standards: implications for the assessment of overweight and obesity. Croat Med J 2018;59:313–26.
16. Lenhard W, Lenhard A. Computation of effect sizes [Internet]. Psychometrica; 2022 [cited 2025 Nov. 11]. Available from: https://www.psychometrica.de/effect_size.html.
17. Gordon CM, Leonard MB, Zemel BS, ; International Society for Clinical Densitometry. 2013 Pediatric Position Development Conference: executive summary and reflections. J Clin Densitom 2014;17:219–24.
18. Salari N, Razavizadeh S, Abdolmaleki A, Heidarian P, Rahimi A, Shohaimi S, et al. The global prevalence of anemia in patients with anorexia nervosa: a systematic review and meta-analysis. BMC Psychol 2025;13:343.
19. Alvin P, Zogheib J, Rey C, Losay J. [Severe complications and mortality in mental eating disorders in adolescence. On 99 hospitalized patients]. Arch Fr Pediatr 1993;50:755–62. French.
20. Devuyst O, Lambert M, Rodhain J, Lefebvre C, Coche E. Haematological changes and infectious complications in anorexia nervosa: a case-control study. Q J Med 1993;86:791–9.
21. Palla B, Litt IF. Medical complications of eating disorders in adolescents. Pediatrics 1988;81:613–23.
22. Dalle Grave R. Severe and enduring anorexia nervosa: no easy solutions. Int J Eat Disord 2020;53:1320–1.
23. Hay P, Touyz S. Classification challenges in the field of eating disorders: can severe and enduring anorexia nervosa be better defined? J Eat Disord 2018;6:41.
24. Walsh K, Blalock DV, Mehler PS. Hematologic findings in a large sample of patients with anorexia nervosa and bulimia nervosa. Am J Hematol 2020;95:E98–101.
25. Hütter G, Ganepola S, Hofmann WK. The hematology of anorexia nervosa. Int J Eat Disord 2009;42:293–300.
26. Guinhut M, Melchior JC, Godart N, Hanachi M. Extremely severe anorexia nervosa: hospital course of 354 adult patients in a clinical nutrition-eating disorders-unit. Clin Nutr 2021;40:1954–65.
27. Sabel AL, Gaudiani JL, Statland B, Mehler PS. Hematological abnormalities in severe anorexia nervosa. Ann Hematol 2013;92:605–13.
28. Ecklund K, Vajapeyam S, Feldman HA, Buzney CD, Mulkern RV, Kleinman PK, et al. Bone marrow changes in adolescent girls with anorexia nervosa. J Bone Miner Res 2010;25:298–304.
29. Cunha MC, Lima Fda S, Vinolo MA, Hastreiter A, Curi R, Borelli P, et al. Protein malnutrition induces bone marrow mesenchymal stem cells commitment to adipogenic differentiation leading to hematopoietic failure. PLoS One 2013;8e58872.
30. Li Z, Rosen CJ. The multifaceted roles of bone marrow adipocytes in bone and hematopoietic homeostasis. J Clin Endocrinol Metab 2023;108:e1465–72.
31. Dang T, Faje AT, Meenaghan E, Bredella MA, Bouxsein ML, Klibanski A, et al. Bone marrow adipose tissue is associated with fracture history in anorexia nervosa. Osteoporos Int 2022;33:2619–27.
32. Abella E, Feliu E, Granada I, Millá F, Oriol A, Ribera JM, et al. Bone marrow changes in anorexia nervosa are correlated with the amount of weight loss and not with other clinical findings. Am J Clin Pathol 2002;118:582–8.
33. Haagensen AL, Feldman HA, Ringelheim J, Gordon CM. Low prevalence of vitamin D deficiency among adolescents with anorexia nervosa. Osteoporos Int 2008;19:289–94.
34. Misra M, Klibanski A. Endocrine consequences of anorexia nervosa. Lancet Diabetes Endocrinol 2014;2:581–92.
35. Mumford J, Kohn M, Briody J, Miskovic-Wheatley J, Madden S, Clarke S, et al. Long-term outcomes of adolescent anorexia nervosa on bone. J Adolesc Health 2019;64:305–10.

Article information Continued

Fig. 1.

Blood counts, markers of bone status, and estradiol and IGF-1 concentrations of patients with AN (n=59) versus controls (n=20). Normally distributed data (panels A, B, C, D, F, and H) were compared using the 2-sided independent Student t test after checking for equality of variances with Levene test and are shown as mean and standard deviation. Skewed data (panels E and G) were compared using the Mann-Whitney U test, and medians and interquartile ranges are indicated. Statistical significance was set at P<0.05. IGF-1, insulin-like growth factor-1; AN, anorexia nervosa; BMD, bone mineral density; PINP, amino-terminal propeptide of human procollagen type I.

Fig. 2.

Blood counts, markers of bone status, and estradiol and IGF-1 concentrations in patients with anorexia nervosa presenting with BMD z scores ≥ -1.4 SD (n=45) and healthy controls (n=20). Normally distributed data (panels A, B, C, D, F, and H) were compared using a 2-sided independent Student t test after a check for equality of variances with Levene test and are given as mean and SD. Skewed data (panels E and G) were compared using the Mann-Whitney U test and are shown as median and interquartile range. Statistical significance was set at P<0.05. IGF-1, insulin-like growth factor-1; AN, anorexia nervosa; BMD, bone mineral density; PINP, amino-terminal propeptide of human procollagen type I; SD, standard deviation.

Table 1.

Anthropometrics and blood counts of healthy controls versus patients with anorexia nervosa

Variable Controls (n=20) AN (n=59) P value Cohen d (95% CI)
Age (yr) 15.7±1.7 15.5±1.8 0.660 ND
Height (SDS) 0.36±1.04 0.45±1.10 0.750 ND
Body weight (SDS) 0.26±0.75 -1.03±0.97 <0.001 -1.4 (-1.9 to -0.9)
Body mass index (SDS) -0.16±0.67 -1.78±1.02 <0.001 -1.7 (-2.3 to -1.1)
Fat (%) 31.7±4.0 24.0±7.9 <0.001 -1.1 (-1.6 to -0.5)
Neutrophils (%) 50.2±8.8 49.8±13.0 0.884 ND
Lymphocytes (%) 37.8±7.4 40.2±12.1 0.440 ND
Monocytes (%) 9.0±1.6 7.5±1.6 0.001 -0.9 (-1.5 to -0.4)
Eosinophils (%) 2.4±1.3 2.1±1.9 0.456 ND
Basophils (%) 0.6±0.3 0.6±0.2 0.161 ND
Hematocrit (%) 40.4±1.3 38.5±3.5 0.001b) -0.6 (-1.1 to -0.1)
Hemoglobin (g/dL) 13.3±0.5 13.2±1.0 0.531 ND
MCV (fL) 86±4 88±3 0.061 ND
MCH (pg) 29±2 30±2 <0.001 1.0 (0.5–1.5)
RDW (%) 12.4±0.5 13.2±1.1 <0.001 0.8 (0.2–1.3)
Plateletcrit (%) 0.27±0.04 0.22±0.05 <0.001 -1.0 (-1.5 to -0.5)
MPV (fL) 10.4±0.7 10.5±1.0 0.663 ND
PDW (%) 12.0±1.3 12.2±2.2 0.518 ND
LS BMD z score -0.08±0.74 -0.20±1.17 0.611 ND
FN BMD z score 0.11±0.77 -0.27±1.13 0.098 ND
LS TBS 1.44±0.06 1.43±0.07 0.320 ND
Osteocalcin (ng/mL) 57.9±25.8 23.3±18.1 <0.001 -1.7 (-2.3 to -1.1)
CTX (ng/mL) 1.19±0.48 1.34±0.77 0.391 ND
ALKP (µkat/L) 1.7±0.9 1.0±0.4 0.002 -1.2 (-1.8 to -0.7)
TSH (mIU/L) 2.8±1.1 2.7±2.0 0.866 ND
fT4 (pmol/L) 16.4±1.4 13.3±2.9 <0.001 -1.2 (-1.7 to -0.6)
Cortisol (nmol/L) 403±85 444±121 0.176 ND
Growth hormone (µg/L) 0.9 (0.5–7.0) 6.7 (2.4–18.2) 0.001 0.9
Folate (nmol/L) ND 17.9±9.6 NA NA
Vitamin B12 (pmol/L) ND 336±220a) NA NA
Duration of illness (mo) NA 13.3±8.8 NA NA
Duration of amenorrhea (mo) NA 11.3±8.7 NA NA

Normally distributed data are shown as mean±SD and were compared using the independent Student t test when equal variances could be assumed and Welch t tests otherwise, based on Levene test. Skewed data are shown as median (interquartile range) and were compared using the Mann-Whitney U test. Effect size is shown as Cohen d (with 95% confidence intervals for normally distributed data).

AN, anorexia nervosa; CI, confidence interval; SDS, standard deviation score; MCV, mean cell volume; MCH, mean cell hemoglobin; RDW, red cell distribution width; MPV, mean platelet volume; PDW, platelet distribution width; LS, lumbar spine; FN, femoral neck; BMD, bone mineral density; TBS, trabecular bone score; CTX, carboxy-terminal telopeptide of type I collagen; ALKP, alkaline phosphatase; TSH, thyroid-stimulating hormone; IU, international units; fT4, free thyroxine; ND, not determined; NA, not applicable.

a)

One patient exhibited levels below the reference range.

b)

Welch t test.

Boldface indicates a statistically significant difference with P<0.05.

Table 2.

Anthropometrics and blood counts of healthy controls versus patients with anorexia nervosa presenting with BMD z scores ≥ -1.4 SD

Variable Control (n=20) AN (n=45) P value Cohen d (95% CI)
Age (yr) 15.7±1.7 15.3±1.8 0.402 ND
Height (SDS) 0.36±1.04 0.54±1.14 0.549 ND
Body weight (SDS) 0.26±0.75 -0.85±0.87 <0.001 -1.3 (-1.9 to -0.7)
Body mass index (SDS) -0.16±0.67 -1.57±1.00 <0.001 -1.5 (-2.1 to -1.0)
Fat (%) 31.7±4.0 24.1±8.3 <0.001 -1.0 (-1.7 to -0.5)
Neutrophils (%) 50.2±8.8 48.2±12.3 0.511 ND
Lymphocytes (%) 37.8±7.4 41.6±11.5 0.181 ND
Monocytes (%) 9.0±1.6 7.5±1.7 0.001 -0.9 (-1.5 to -0.4)
Eosinophils (%) 2.4±1.3 2.1±1.9 0.432 ND
Basophils (%) 0.6±0.3 0.7±0.3 0.121 ND
Hematocrit (%) 40.4±1.3 38.8±3.3 0.001 -0.6 (-1.1 to -0.0)
Hemoglobin (g/dL) 13.3±0.5 13.3±1.0 0.630 ND
MCV (fL) 86±4 88±3 0.087 ND
MCH (pg) 29±2 30±2 0.001 0.9 (0.4–1.5)
RDW (%) 12.4±0.5 13.3±1.1 <0.001 0.8 (0.3–1.3)
Plateletcrit (%) 0.27±0.04 0.23±0.05 0.001 -0.9 (-1.5 to -0.4)
MPV (fL) 10.4±0.7 10.5±1.0 0.647 ND
PDW (%) 12.0±1.3 12.2±2.2 0.518 ND
LS BMD z score -0.08±0.74 0.18±0.93 0.276 ND
FN BMD z score 0.11±0.77 0.12±0.96 0.967 ND
LS TBS 1.44±0.06 1.44±0.07 0.933 ND
Osteocalcin (ng/mL) 57.9±25.8 25.1±20.3 <0.001 -1.5 (-2.1 to -0.9)
CTX (ng/mL) 1.19±0.48 1.29±0.75 0.569 ND
ALKP (µkat/L) 1.7±0.9 1.1±0.4 <0.001 -1.0 (-1.6 to -0.5)
TSH (mIU/L) 2.8±1.1 2.9±2.2 0.900 ND
fT4 (pmol/L) 16.4±1.4 13.2±2.7 <0.001 -1.3 (-1.9 to -0.8)
Cortisol (nmol/L) 403±85 442±127 0.227 ND
Growth hormone (µg/L) 0.9 (0.5; 7.0) 4.4 (1.8; 15.7) 0.004 0.8
Folate (nmol/L) ND 18.3±9.9 NA NA
Vit. B12 (pmol/L) ND 350±243a) NA NA
Duration of illness (mo) NA 13.5±9.2 NA NA
Duration of amenorrhea (mo) NA 10.5±7.8 NA NA

Normally distributed data are shown as mean±standard deviation and were compared using the independent Student t test. Skewed data are shown as median (interquartile range) and were compared using the Mann-Whitney U test. Effect size is shown as Cohen d (with 95% confidence intervals for normally distributed data).

SDS, standard deviation score; AN, anorexia nervosa; CI, confidence interval; MCV, mean cell volume; MCH, mean cell hemoglobin; RDW, red cell distribution width; MPV, mean platelet volume; PDW, platelet distribution width; LS, lumbar spine; FN, femoral neck; BMD, bone mineral density; TBS, trabecular bone score; CTX, carboxy-terminal telopeptide of type I collagen; ALKP, alkaline phosphatase; TSH, thyroid-stimulating hormone; IU, international units; fT4, free thyroxine; ND, not determined; NA, not applicable.

a)

One patient exhibited levels below the reference range.

Boldface indicates a statistically significant difference with P<0.05.