Oświadczenie etyczne: Badanie zostało przeprowadzone zgodnie z wytycznymi instytucjonalnymi oraz zasadami etycznymi Deklaracji Helsińskiej. Zgoda etyczna została uzyskana od Komitetu Etycznego Uniwersytetu w Bisha (UBCOM-RELOC), numer zgody H-06-BH-087/(0802.23). Udział w badaniu był dobrowolny, a od wszystkich uczestników przed włączeniem do grupy uzyskano pisemną świadomą zgodę. Materiały, sprzęt i procedury analityczne wykorzystane w badaniu podsumowano w Tabeli materiałów.
Projekt i miejsce badania
W okresie od grudnia 2022 do maja 2023 roku w Szpitalu dla Matek i Dzieci w Bisha w Arabii Saudyjskiej przeprowadzono retrospektywne badanie przekrojowe w celu oceny stanu zapasów żelaza oraz zbadania czynników związanych z niedoborem żelaza u kobiet w ciąży w północnym regionie Asir.
Populacja badawcza i metodologia próbkowania
Do badania zakwalifikowano łącznie 200 uczestniczek, w tym 135 kobiet w ciąży korzystających z rutynowej opieki przedporodowej oraz 65 kobiet niebędących w ciąży, dobranych pod względem wieku, które stanowiły grupę kontrolną. Uczestniczki zostały zrekrutowane z poradni ambulatoryjnych przy użyciu doboru wygodnego w okresie trwania badania. Do kobiet spełniających kryteria kwalifikacji, zgłaszających się w godzinach pracy kliniki, podchodzono sukcesywnie, a rekrutację kontynuowano do momentu osiągnięcia docelowej wielkości próby. Aby zminimalizować potencjalne błędy selekcji w ramach tego modelu, grupę kontrolną dobrano pod względem wieku, stosując proste losowanie spośród kwalifikujących się pacjentek niebędących w ciąży.
Kobiety w ciąży kwalifikowały się do badania, jeśli miały co najmniej 18 lat i wiek ciążowy między 12 a 36 tygodniem. Kwalifikujące się uczestniczki włączano do badania niezależnie od poziomu żelaza, a następnie klasyfikowano jako kobiety w ciąży z IDA lub bez IDA na podstawie wyników laboratoryjnych. Zgodnie z kryteriami WHO, anemię w czasie ciąży definiowano jako stężenie hemoglobiny (Hb) < 11,0 g/dL w pierwszym i trzecim trymestrze oraz < 10,5 g/dL w drugim trymestrze. IDA diagnozowano u kobiet, które spełniały próg Hb specyficzny dla ciąży i wykazywały niedobór żelaza, w tym poziom ferrytyny w surowicy < 15 ng/mL, żelaza w surowicy < 10 µmol/L oraz całkowitą zdolność wiązania żelaza (TIBC) ≥ 68 µmol/L17,18. Kobiety w ciąży, które nie spełniały tych kryteriów diagnostycznych, zaklasyfikowano do grupy kobiet w ciąży bez IDA. Grupę kontrolną stanowiły kobiety w wieku rozrodczym, niebędące w ciąży i bez znanych zaburzeń hematologicznych. Kryteria wykluczenia dla obu grup obejmowały udokumentowaną historię przewlekłych lub systemowych schorzeń związanych z zaburzeniami metabolizmu żelaza lub stanem zapalnym, takich jak przewlekła choroba nerek, choroby autoimmunologiczne, nowotwory, zespół zapalniczy jelit, przewlekła choroba wątroby lub dziedziczne zaburzenia gospodarki żelazem (np. hemochromatoza). Wykluczone zostały również kobiety z ostrymi infekcjami w momencie pobierania próbek, które w ciągu trzech poprzedzających miesięcy otrzymały dożylną terapię żelazem lub transfuzje krwi, a także osoby przyjmujące długoterminową terapię immunosupresyjną. Dla wszystkich uczestniczek odnotowano wskaźnik masy ciała (BMI) i przynależność etniczną; zmienne te nie były jednak stosowane jako kryteria wykluczenia, lecz traktowano je jako potencjalne kowarianty w analizie.
Wyznaczenie wielkości próby
Wielkość próby obliczono przy użyciu oprogramowania do epidemiologicznych obliczeń wielkości próby19, w oparciu o przewidywaną 40-procentową zapadalność na IDA wśród kobiet w ciąży w Arabii Saudyjskiej6, przy 95% poziomie ufności i 5% marginesie błędu. Minimalną wymaganą wielkość próby oszacowano na 182 uczestniczki; w związku z tym zrekrutowano łącznie 200 uczestniczek, aby zwiększyć precyzję statystyczną.
Zbieranie danych
Kwestionariusz
Dane zbierano za pomocą ustrukturyzowanego i zwalidowanego kwestionariusza. Kwestionariusz opracowano po szczegółowym przeglądzie literatury, a następnie został on oceniony przez panel ekspertów z danej dziedziny, w tym położników, dietetyków i epidemiologów, aby zapewnić trafność treści. Przeprowadzono badanie pilotażowe z udziałem 20 uczestników, którzy nie zostali włączeni do końcowej analizy, w celu oceny jasności, niezawodności i wykonalności. Spójność wewnętrzna kwestionariusza była na akceptowalnym poziomie, z alfą Cronbacha wynoszącą 0,81.
Końcowy kwestionariusz składał się z trzech sekcji obejmujących łącznie 35 pozycji (Supplementary File 1):
Sekcja 1 (Dane socjodemograficzne, 10 pozycji): wiek, poziom wykształcenia, zawód, miejsce zamieszkania, narodowość oraz wywiad rodzinny w kierunku anemii.
Sekcja 2 (Nawyki żywieniowe, 12 pozycji): częstotliwość spożycia produktów bogatych w żelazo (czerwone mięso, białe mięso, ryby, warzywa i wzbogacane zboża), spożycie herbaty i kawy oraz reżim dietetyczny. Zmienne dietetyczne wybrano a priori na podstawie ich potwierdzonego związku z przyjęciem żelaza, biodostępnością żelaza oraz czynnikami dietetycznymi, które mogą wpływać na wchłanianie żelaza. Kwestionariusz obejmował główne źródła żelaza w diecie, w tym czerwone mięso, drób, ryby, zielone warzywa liściaste i wzbogacane zboża. Oceniano spożycie produktów bogatych w witaminę C, szczególnie cytrusów, ponieważ witamina C może zwiększać wchłanianie żelaza niehemowego. Z kolei spożycie herbaty i kawy oceniano ze względu na to, że polifenole dietetyczne mogą hamować wchłanianie żelaza niehemowego, jeśli są spożywane w trakcie lub w bliskim odstępie czasu od posiłków. Oceniano również reżim dietetyczny i spożycie fast foodów, aby określić ogólne wzorce żywieniowe i adekwatność żywieniową. Wybrane produkty spożywcze zostały dobrane na podstawie literatury analizowanej podczas opracowywania kwestionariusza, a następnie ocenione pod kątem trafności treści przez ekspertów z zakresu położnictwa, żywienia i epidemiologii.
Sekcja 3 (Wywiad kliniczny i objawy, 13 elementów): wiek ciążowy, liczba porodów, historia transfuzji krwi, zmęczenie, zawroty głowy, gorączka oraz inne objawy powszechnie związane z anemią.
Ocena biochemiczna
Od każdego uczestnika pobrano 5 mL krwi żylnej po nocnym poście. Próbki wirowano z prędkością 2000 × g przez 15 min w celu oddzielenia surowicy, którą następnie poddano natychmiastowej analizie. Stężenie hemoglobiny zmierzono za pomocą automatycznego analizatora hematologicznego i podano w gramach na decylitr (g/dL). Poziom ferrytyny w surowicy oceniono za pomocą testu immunochemicznego z emisją elektrochemiluminescencyjną (ECLIA) na automatycznym analizatorze immunoassay i wyrażono w nanogramach na mililitr (ng/mL). Stężenie żelaza w surowicy, całkowita zdolność wiązania żelaza (TIBC) oraz nienasycona zdolność wiązania żelaza (UIBC) zostały określone metodą kolorymetryczną na automatycznym analizatorze biochemicznym, a wyniki podano w mikromolach na litr (µmol/L). Wysycenie transferyny (TR-SAT) obliczono jako stosunek żelaza w surowicy do TIBC pomnożony przez 100 i wyrażono w procentach.
Statystyka i analiza danych
Dane zakodowano i przeanalizowano przy użyciu standardowego oprogramowania do analizy statystycznej. Do podsumowania zmiennych socjodemograficznych, dietetycznych, klinicznych i biochemicznych wykorzystano statystykę opisową, w tym średnie, odchylenia standardowe i wartości procentowe. Porównania grup przeprowadzono za pomocą nieparzystego testu t Student’a lub jednoczynnikowej analizy wariancji (ANOVA), w zależności od potrzeb. Istotność statystyczną określono jako p < 0,05. Wyniki graficzne wygenerowano przy użyciu oprogramowania do tworzenia wykresów i analizy statystycznej.
W celu zbadania czynników niezależnie związanych z IDA zastosowano wieloczynnikowy model regresji logistycznej. Status IDA określono jako binarną zmienną zależną (IDA = 1, brak IDA/grupa kontrolna = 0). Zmienne niezależne wybrano a priori na podstawie prawdopodobieństwa biologicznego oraz dowodów z poprzednich badań; obejmowały one charakterystykę socjodemograficzną (miejsce zamieszkania, wywiad rodzinny w kierunku anemii), stan fizjologiczny (ciąża) oraz zmienne dietetyczne (spożycie czerwonego mięsa, drobiu i owoców cytrusowych).
Zmienne kategoryczne zakodowano za pomocą zmiennych pozornych, a zmienne dotyczące spożycia pokarmów zaklasyfikowano jako niskie lub odpowiednie na podstawie zgłoszonej częstotliwości spożycia (≤2 razy w tygodniu vs. >2 razy w tygodniu). Wszystkie wybrane zmienne wprowadzono jednocześnie do modelu regresji, aby oszacować ich niezależne powiązania przy jednoczesnym uwzględnieniu potencjalnych czynników zakłócających. Raportowano skorygowane ilorazy szans (OR) wraz z 95% przedziałami ufności (CI), gdzie OR większe od 1 odzwierciedlały wyższe prawdopodobieństwo wystąpienia IDA. Graficzne podsumowanie wartości OR oraz odpowiadających im 95% CI przedstawiono za pomocą wykresu leśnego (forest plot) wygenerowanego w programie arkusza kalkulacyjnego.