Ethical statement: The study was performed in accordance with the institutional guidelines and the ethical principles of the Declaration of Helsinki. Ethical approval was obtained from the Ethical Committee of Bisha University (UBCOM-RELOC), with approval number H-06-BH-087/(0802.23). Participation was voluntary, and written informed consent was obtained from all participants prior to enrollment. The materials, equipment, and analytical procedures used in the study are summarized in the Table of Materials.
Study design and setting
A retrospective cross-sectional study was conducted at the Maternity and Children’s Hospital in Bisha, Saudi Arabia, between December 2022 and May 2023 to assess iron status and examine factors associated with iron deficiency among pregnant women in the Northern Asir Region.
Study population and sampling methodology
A total of 200 participants were enrolled, including 135 pregnant women attending routine antenatal care clinics and 65 age-matched non-pregnant women who served as a control group. Participants were recruited from outpatient clinics via convenience sampling during the study period. Eligible women presenting during clinic hours were consecutively approached, and recruitment continued until the target sample size was reached. To reduce potential selection bias within this framework, controls were age-matched using simple random sampling from eligible non-pregnant clinic attendees.
Pregnant women were eligible for inclusion if they were aged 18 years or older and had a pregnancy with a gestational age between 12 and 36 weeks. Eligible participants were enrolled regardless of their iron status and subsequently classified as pregnant women with IDA or without IDA based on laboratory findings. According to the WHO criteria, anemia during pregnancy was defined as hemoglobin (Hb) concentrations < 11.0 g/dL during the first and third trimesters and < 10.5 g/dL during the second trimester. IDA was diagnosed in women who met the pregnancy-specific Hb threshold and had evidence of iron deficiency, including serum ferritin < 15 ng/mL, serum iron < 10 µmol/L, and total iron-binding capacity (TIBC) ≥ 68 µmol/L17,18. Pregnant women who did not meet these diagnostic criteria were classified as the pregnant without IDA group. Non-pregnant women of reproductive age with no current pregnancy and without known hematological disorders served as the control group. Exclusion criteria for both groups included a documented history of chronic or systemic conditions associated with altered iron metabolism or inflammatory status, such as chronic kidney disease, autoimmune disorders, malignancy, inflammatory bowel syndrome, chronic liver disease, or hereditary iron disorders (e.g., hemochromatosis). Women with acute infections at the time of sampling, those who had received intravenous iron therapy or blood transfusions within the preceding three months, and those receiving long-term immunosuppressive therapy were also excluded. Body mass index (BMI) and ethnicity were recorded for all participants; however, these variables were not used as exclusion criteria and were instead considered as potential covariates in the analysis.
Sample size determination
The sample size was calculated using an epidemiological sample-size calculation software19, based on an anticipated prevalence of IDA among pregnant women in Saudi Arabia of 40%6, with a 95% confidence level and a 5% margin of error. The minimum required sample size was estimated at 182 participants; therefore, a total of 200 participants were recruited to enhance statistical precision.
Data collection
Questionnaire
Data was collected using a structured and validated questionnaire. The questionnaire was developed following an extensive review of the literature and was evaluated by a panel of subject experts, including obstetricians, nutritionists, and epidemiologists, to ensure content validity. A pilot study involving 20 participants, who were not included in the final analysis, was conducted to assess clarity, reliability, and feasibility. The internal consistency of the questionnaire was acceptable, with Cronbach’s alpha of 0.81.
The final questionnaire consisted of three sections comprising a total of 35 items (Supplementary File 1):
Section 1 (Sociodemographic data, 10 items): age, education level, occupation, place of residence, nationality, and family history of anemia.
Section 2 (Dietary habits, 12 items): frequency of consumption of iron-rich foods (red meat, white meat, fish, vegetables, and fortified cereals), tea and coffee intake, and dietary regimen. Dietary variables were selected a priori based on their established relevance to iron intake, iron bioavailability, and dietary factors that may influence iron absorption. The questionnaire included major dietary sources of iron, including red meat, poultry, fish, green leafy vegetables, and fortified cereals. Vitamin C-rich foods, particularly citrus fruits, were assessed because vitamin C can enhance the absorption of non-heme iron. In contrast, tea and coffee consumption were assessed because dietary polyphenols may inhibit non-heme iron absorption when consumed with or near meals. Dietary regimen and fast-food consumption were also assessed to indicate overall dietary patterns and nutritional adequacy. The selected dietary items were informed by literature reviewed during questionnaire development and were subsequently evaluated by obstetric, nutrition, and epidemiology experts for content validity.
Section 3 (Clinical history and symptoms, 13 items): gestational age, parity, history of blood transfusion, fatigue, dizziness, fever, and other symptoms commonly associated with anemia.
Biochemical assessment
For each participant, a 5 mL venous blood sample was obtained after an overnight fast. Samples were centrifuged at 2000 × g for 15 min to separate the serum, which was subsequently analyzed immediately. Hemoglobin concentration was measured using an automated hematology analyzer and reported in grams per deciliter (g/dL). Serum ferritin levels were assessed using an electrochemiluminescence immunoassay (ECLIA) on an automated immunoassay analyzer and expressed in nanograms per milliliter (ng/mL). Serum iron, total iron-binding capacity (TIBC), and unsaturated iron-binding capacity (UIBC) were determined using a colorimetric method on an automated biochemical analyzer, with results reported in micromoles per liter (µmol/L). Transferrin saturation (TR-SAT) was calculated as the ratio of serum iron to TIBC multiplied by 100 and expressed as a percentage.
Statistics and data analyses
Data was coded and analyzed using standardized statistical analysis software. Descriptive statistics, including means, standard deviations, and percentages, were used to summarize sociodemographic, dietary, clinical, and biochemical variables. Group comparisons were conducted using an unpaired Student’s t-test or one-way analysis of variance (ANOVA), as appropriate. Statistical significance was defined as p < 0.05. Graphical outputs were produced using graphing and statistical software.
To examine factors independently associated with IDA, a multivariate logistic regression model was applied. IDA status was specified as the binary dependent variable (IDA = 1, non-IDA/control = 0). Independent variables were selected a priori based on biological plausibility and evidence from previous studies and included sociodemographic characteristics (residence, family history of anemia), physiological status (pregnancy), and dietary variables (intake of red meat, poultry, and citrus fruits).
Categorical variables were coded using dummy variables, and dietary intake variables were classified as low or adequate based on reported consumption frequency (≤2 times per week vs. >2 times per week). All selected variables were entered simultaneously into the regression model to estimate their independent associations while accounting for potential confounding factors. Adjusted odds ratios (ORs) with 95% confidence intervals (CIs) were reported, where ORs greater than 1 reflected a higher odds of IDA. A visual summary of the ORs and corresponding 95% CIs was presented using a forest plot generated in spreadsheet software.