Badanie to zostało przeprowadzone zgodnie z Deklaracją Helsińską oraz wymogami instytucjonalnymi Szpitala Medycyny Tradycyjnej Chińskiej w Yongkang10. Badanie zostało zatwierdzone przez Komisję Etyczną Szpitala Medycyny Tradycyjnej Chińskiej w Yongkang w 2021 roku. W dokumencie zatwierdzającym nie przypisano oddzielnego numeru zatwierdzenia. Przed screeningiem i rekrutacją od wszystkich uczestników uzyskano pisemną świadomą zgodę. Uczestnicy zostali poinformowani, że udział jest dobrowolny, odmowa nie wpłynie na rutynową opiekę dializacyjną, a wycofanie się z badania jest możliwe w dowolnym momencie. Odczynniki, chemikalia i narzędzia użyte w protokole wymieniono w Tabeli materiałów.
1. Projekt badania i ogólny przebieg prac
Niniejsze prospektywne, jednoośrodkowe, otwarte, randomizowane badanie kontrolowane z grupami równoległymi oceniało, czy połączona interwencja olfaktoryczna z wykorzystaniem preparatu z olejków eterycznych i tradycyjnego chińskiego woreczka ziołowego wiąże się z poprawą subiektywnie ocenianej jakości snu i zmniejszeniem lęku u pacjentów poddawanych hemodializie utrzymującej. Uczestników zrekrutowano z Centrum Oczyszczania Krwi w Szpitalu Medycyny Tradycyjnej w Yongkang w okresie od stycznia 2022 do grudnia 2023 roku. Badanie to nie zostało prospektywnie zarejestrowane w publicznym rejestrze badań klinicznych. Raportowanie z randomizowanego badania oparto na zasadach CONSORT, tam gdzie było to możliwe11.
Potencjalnie kwalifikujących się pacjentów przesiewano podczas rutynowych wizyt dializacyjnych. Po uzyskaniu pisemnej świadomej zgody, przed przydzieleniem do grup, zebrano podstawowe dane demograficzne, informacje związane z dializą, historię przyjmowanych leków, wyniki w Pittsburgh Sleep Quality Index (PSQI) oraz wyniki w podskali lęku Hospital Anxiety and Depression Scale (HADS-A). Kwalifikujący się uczestnicy zostali randomizowani w stosunku 1:1 do grupy kontrolnej lub grupy interwencyjnej. Randomizowano łącznie 88 uczestników, z czego 44 przypisano do każdej z grup. Schemat przepływu uczestników przedstawiono na Rycini 1.

Rycina 1: Schemat przebiegu badania. Przepływ uczestników podczas przesiewu, wykluczania, randomizacji, przydziału, obserwacji i analizy. Oceniono łącznie 104 pacjentów, 16 wykluczono, a 88 zrandomizowano w stosunku 1:1 do grupy kontrolnej lub interwencyjnej. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Wielkość próby obliczono, przyjmując globalny wynik PSQI w 3. miesiącu jako główny punkt końcowy. Za klinicznie istotną uznano różnicę międzygrupową wynoszącą 2,5 punktu. Zakładając wspólne odchylenie standardowe wynoszące 3,5 punktu, dwustronny poziom α wynoszący 0,05, moc testu na poziomie 80% oraz przewidywany odsetek osób rezygnujących z badania wynoszący 10%, wymagano co najmniej 40 uczestników w każdej grupie. Ostateczną docelową wielkość próby ustalono na 44 uczestników w każdej grupie.
2. Rekrutacja i kryteria kwalifikowalności uczestników
Uczestników rekrutowano poprzez ciągły screening podczas rutynowych wizyt na dializach. Personel badawczy omówił cel badania, procedury interwencyjne, harmonogram obserwacji, potencjalny dyskomfort, zasady poufności oraz prawo do wycofania się z badania przed uzyskaniem zgody.
Kryteria włączenia uczestników obejmowały: poddawanie się dializom hemodializacyjnym w trybie podtrzymującym przez co najmniej 12 miesięcy, stabilny harmonogram dializ (2–4 sesje w tygodniu, ok. 4 h na sesję), wiek od 30 do 65 lat, zdolność do współpracy w ramach procedur badania, możliwość samodzielnego wypełnienia kwestionariuszy lub wypełnienia ich przy ustandaryzowanej pomocy, wyjściowy globalny wynik PSQI > 7 oraz brak stosowania leków przeznaczonych specjalnie do poprawy lub modyfikacji snu w ciągu ostatnich 2 tygodni.
Kryteria wykluczenia obejmowały astmę, czynne zapalenie błony śluzowej nosa, zaburzenia węchu, znaną nadwrażliwość na preparaty aromatyczne lub woreczki ziołowe, poważne choroby hematologiczne, ciężkie niestabilne choroby organiczne, niedawne zabiegi chirurgiczne lub urazy, ciężkie infekcje, ostra niestabilność psychiatryczna wymagająca specjalistycznej opieki, znaczne ograniczenie mobilności, poważne upośledzenie słuchu lub wzroku, ciężkie upośledzenie poznawcze lub niekompletne podstawowe dane kliniczne.
Uczestników wycofywano z badania w przypadku zaprzestania hemodializy utrzymującej, zgonu w okresie obserwacji, wycofania zgody, utraty kontaktu w trakcie obserwacji lub wystąpienia poważnego zdarzenia medycznego związanego z dializą, które uniemożliwiało dalszy udział. Odnotowywano datę oraz przyczynę wycofania. Kompletne dane zebrane przed przerwaniem udziału zachowywano w bazie danych badania.
3. Randomizacja, maskowanie alokacji i maskowanie
Po ocenie wyjściowej uczestnicy zostali losowo przydzieleni w stosunku 1:1 przy użyciu tablicy liczb losowych przygotowanej przed rekrutacją. Randomizacja odbywała się w stałych blokach czwórkowych, bez stratyfikacji. Sekwencja przydziału została przygotowana przez członka personelu, który nie brał udziału w przesiewie, realizacji interwencji, ocenie wyników, wprowadzaniu danych ani w analizie statystycznej.
Przydziały do grup umieszczono w ponumerowanych sekwencyjnie, nieprzezroczystych, zapieczętowanych kopertach. Po potwierdzeniu kwalifikowalności i ocenie wyjściowej, otwierano kolejną kopertę zgodnie z kolejnością numeryczną. Formularze wyników wyjściowych wypełniono przed ujawnieniem przydziału.
Ślepa próba w odniesieniu do uczestników i osób prowadzących interwencję nie była możliwa do zastosowania, ponieważ interwencja obejmowała ekspozycję na zapachy oraz zewnętrzne stosowanie substancji aromatycznych. Aby zminimalizować błąd oczekiwań i błąd raportowania, obie grupy otrzymały zestandaryzowane wyjaśnienia, a kwestionariusze były przeprowadzane przy użyciu tej samej treści pytań, kolejności elementów i wytycznych dotyczących odpowiedzi w każdym punkcie czasowym. Formularze wyników były sprawdzane przez personel nieodpowiedzialny za instruktaż interwencji, gdy pozwalały na to zasoby kadrowe. Dane wprowadzano przy użyciu kodów uczestników. Przed analizą etykiety grup zostały zmienione na Grupę A i Grupę B, a statystyk przeanalizował zamaskowany zbiór danych przed przywróceniem tożsamości grup.
4. Zarządzanie warunkami kontrolnymi i współinterwencjami
Grupa kontrolna otrzymała rutynową opiekę w ramach hemodializy utrzymującej oraz standardową edukację zdrowotną. Opieka rutynowa obejmowała dializy zgodnie z istniejącym harmonogramem, obserwacje pielęgniarskie, pielęgnację dostępu naczyniowego, wytyczne dotyczące diety i przyjmowania płynów, przypominanie o lekach oraz edukację w zakresie higieny snu. Edukacja w zakresie higieny snu obejmowała regularne godziny snu, ograniczenie stymulujących aktywności przed pójściem spać, unikanie nadmiernych drzemek w ciągu dnia, ograniczenie spożycia kofeiny lub mocnej herbaty wieczorem oraz utrzymanie cichego otoczenia podczas snu.
Uczestników z obu grup poproszono, aby podczas okresu obserwacji nie rozpoczynali nowej aromaterapii, stosowania saszetek ziołowych, strukturalnego masażu, refleksologii stóp, akupunktury, moksoterapii, programów medytacyjnych ani innych uzupełniających interwencji wspomagających sen. W celach badania nie wdrażano leczenia lekami sedatywnymi, hipnotykami, anksjolitykami, lekami przeciwdepresyjnymi ani przeciwbólowymi. Dopuszczono i odnotowano zmiany w lekach wymaganych klinicznie lub w schematach dializ.
5. Przygotowanie połączonej interwencji węchowej
Interwencja polegała na zastosowaniu preparatu z olejku eterycznego do inhalacji podczas sesji dializ oraz do aplikacji zewnętrznej przed snem, a także tradycyjnego chińskiego woreczka z ziołami w celu nocnej ekspozycji zapachowej.
Preparat z olejków eterycznych zawierał olejki lawendowy, bergamotkowy i kokosowy w stosunku objętościowym 5:3:10. Dla każdej partii 5 mL olejku eterycznego z lawendy, 3 mL olejku eterycznego z bergamotki oraz 10 mL oleju kokosowego wymieszano w szczelnym pojemniku chroniącym przed światłem poprzez delikatne odwracanie przez 2 min. Mieszaninę przechowywano w temperaturze pokojowej, z dala od światła, ciepła i wilgoci, a następnie rozlano do identycznych bursztynowych buteleczek. Każda buteleczka została oznaczona numerem partii, datą przygotowania, stosunkiem składników, numerem uczestnika badania oraz zalecanym okresem stosowania. Dokumentacja wydawania preparatu przypisywała każdego uczestnika do konkretnej partii.
Saszetka ziołowa zawierała po 6 g Albizia julibrissin, róży, miętę, Agastache rugosa, Acorus tatarinowii, goździki, Angelica dahurica oraz 30 g Artemisia argyi. Susze przetworzono na gruby proszek, wymieszano i zapakowano do przepuszczających powietrze bawełnianych saszetek o wymiarach 8 cm × 10 cm. Saszetki oznaczono datą przygotowania i wymiany, a następnie wymieniano co 30 dni.
Przed rozpoczęciem interwencji uczestnicy przeszli instruktaż bezpośredni dotyczący stosowania oleju, rozmieszczenia saszetek, codziennego prowadzenia zapisów oraz zgłaszania dyskomfortu. Pierwsze zastosowanie odbyło się pod nadzorem w oddziale dializ.
6. Przeprowadzenie interwencji podczas dializy i przed snem
W dniach dializ uczestnicy z grupy interwencyjnej stosowali preparat z olejków eterycznych raz podczas sesji dializy. Personel badawczy nanosił dwie krople preparatu na sterylny wacik bawełniany i umieszczał go w odległości około 5–8 cm od nozdrzy, nie dotykając nosa, skóry, miejsca dostępu naczyniowego, drenów dializacyjnych ani błony śluzowej. Każda ekspozycja trwała 45 min. Jeśli zapach stał się wyraźnie niewyczuwalny przed upływem 45 min, wacik bawełniany mógł zostać wymieniony raz; podczas jednej sesji dializy dopuszczano nie więcej niż jedną wymianę. Personel rejestrował czas rozpoczęcia i zakończenia, fakt wymiany wacika, zakończenie sesji oraz wystąpienie dyskomfortu.
Przed snem uczestnicy aplikowali łącznie 4 krople: 1 kroplę za każdym uchem oraz 1 kroplę na punkt Yongquan każdej stopy. Przez 15 min wykonywany był delikatny masaż okrężny. Uczestników poinstruowano, aby nie nakładać olejku na uszkodzoną skórę, obszary zapalne, błony śluzowe ani miejsca wywołujące dyskomfort. W przypadku wystąpienia podrażnienia, pieczenia, dyskomfortu w oddychaniu, zawrotów głowy lub innego wyraźnego dyskomfortu, stosowanie przerwywano i zgłaszano ten fakt.
Przed snem woreczek ziołowy umieszczano obok poduszki, w odległości około 20–30 cm od głowy. Uczestnicy zostali poinstruowani, aby nie zakrywać nosem ani ustami woreczka. Woreczki wymieniano co 30 dni. Schemat interwencji przedstawiono na Rysunku 2A-C.

Rysunek 2: Ścieżka i harmonogram interwencji. Olfaktoryczna interwencja w dwóch środowiskach, prowadzona podczas dializy oraz przed snem w domu. (A) Inhalacja olejków eterycznych podczas sesji dializy z użyciem dwóch kropli na sterylnym waciku bawełnianym umieszczonym 5–8 cm od nozdrzy przez 45 min. (B) Zewnętrzna aplikacja przed snem z użyciem łącznie czterech kropli, a następnie 15 min delikatnego masażu. (C) Umieszczenie ziołowego woreczka przy poduszce na noc, z wymianą raz w miesiącu. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
Uczestnicy prowadzili codzienny dziennik interwencji, w którym dokumentowano inhalacje podczas sesji dializ, aplikację oleju i masaż przed snem, stosowanie saszetek, wymianę saszetek oraz wystąpienie objawów niepożądanych. Logi były sprawdzane podczas wizyt kontrolnych lub telefonicznego monitorowania. Spadek przestrzegania zaleceń skutkował przeprowadzeniem standaryzowanego ponownego instruktażu.
7. Ocena wyników
Głównym punktem końcowym była jakość snu, oceniana za pomocą Pittsburgh Sleep Quality Index12. PSQI zawiera 18 pozycji wypełnianych przez pacjenta, pogrupowanych w siedem wyników składowych, a globalny wynik mieści się w przedziale 0 – 21, gdzie wyższe wyniki wskazują na gorszą jakość snu. Głównym punktem końcowym był globalny wynik PSQI w 3. miesiącu. Wyniki składowe zachowano jako opisowe wtórne wyniki w dziedzinie snu. Do określenia klinicznie istotnych zaburzeń snu w celu kwalifikacji zastosowano wyjściowy globalny wynik PSQI > 7.
Wtórnym punktem końcowym był lęk, oceniany za pomocą podskali lęku w Hospital Anxiety and Depression Scale (HADS-A)13. HADS-A składa się z siedmiu pozycji, a wynik podskali mieści się w zakresie 0–21, przy czym wyższe wyniki wskazują na nasilenie objawów lękowych.
Kwestionariusze PSQI oraz HADS-A były wypełniane w punkcie wyjściowym, w 7. dniu, 15. dniu, w 1. miesiącu oraz w 3. miesiącu przez przeszkolony personel korzystający ze zestandaryzowanego skryptu. Ponieważ PSQI odzwierciedla stan snu z poprzedniego miesiąca, wartości PSQI z 7. i 15. dnia były traktowane jako wczesne oceny monitorowania objawów, natomiast wynik z 3. miesiąca został wykorzystany do głównego porównania PSQI.
W przypadku braku odpowiedzi na pytanie w kwestionariuszu, personel wyjaśniał tę kwestię z uczestnikiem w tym samym dniu, whenever possible. Jeśli z powodu zbyt wielu brakujących odpowiedzi nie można było wiarygodnie obliczyć wyniku całkowitego, wynik dla tego punktu czasowego traktowano jako brakujący i odnotowywano go w rejestrze brakujących danych.
Odpowiedź w zakresie snu została zdefiniowana jako obniżenie całkowitego wyniku w skali PSQI o co najmniej 3 punkty w stosunku do wartości wyjściowej. Odpowiedź w zakresie lęku została zdefiniowana jako obniżenie wyniku w skali HADS-A o co najmniej 2 punkty w stosunku do wartości wyjściowej. Wyniki te w zakresie odpowiedzi pacjentów zostały przeanalizowane jako wtórne punkty końcowe odpowiedzi klinicznej.
8. Monitorowanie przestrzegania protokołu i bezpieczeństwa
Przestrzeganie zaleceń obliczono jako liczbę wykonanych planowanych procedur interwencyjnych podzieloną przez liczbę procedur przewidzianych w okresie obserwacji, pomnożoną przez 100. Za odpowiednie przestrzeganie zaleceń uznano wykonanie co najmniej 80% planowanych procedur. Codzienne zapisy weryfikowano z dokumentacją wydawania leków, rejestrami wymiany saszetek oraz formularzami kontaktu w ramach obserwacji.
Monitorowanie bezpieczeństwa obejmowało zawroty głowy, bóle głowy, nudności, podrażnienia skóry, subiektywny dyskomfort oddechowy oraz inne objawy występujące podczas lub po inhalacji, aplikacji zewnętrznej lub kontakcie z saszetką. Każde zdarzenie zapisywano pod kątem czasu wystąpienia, czasu trwania, nasilenia, związku z interwencją, sposobu postępowania oraz wyniku. Nasilenie oceniano jako łagodne, umiarkowane lub ciężkie. Związek z interwencją był klasyfikowany przez dwóch klinicystów biorących udział w badaniu jako niepowiązany, potencjalnie powiązany, prawdopodobnie powiązany lub jednoznacznie powiązany. Poważne zdarzenia niepożądane zdefiniowano jako zdarzenia skutkujące hospitalizacją, zagrażającym życiu pogorszeniem stanu zdrowia, trwałym przerwaniem interwencji lub innymi poważnymi konsekwencjami medycznymi.
9. Procedury kontrolne i kontrola jakości danych
Oceny kontrolne przeprowadzono w 7. dniu, 15. dniu, 1. miesiącu i 3. miesiącu po rozpoczęciu interwencji. Dopuszczalne okna czasowe wynosiły odpowiednio 7 ± 2 dni, 15 ± 3 dni, 30 ± 5 dni i 90 ± 10 dni. Odnotowano faktyczną datę oceny.
Kwestionariusze były wypełniane na miejscu podczas dializ, whenever było to możliwe. W przypadku braku możliwości przeprowadzenia obserwacji na miejscu, podczas wywiadów telefonicznych stosowano tę samą treści i strukturę odpowiedzi. Odpowiedzi były wprowadzane przez przeszkolony personel i sprawdzane tego samego dnia.
Wszyscy pracownicy badania przeszli ujednolicone szkolenie przed rozpoczęciem rekrutacji. Wypełnione kwestionariusze sprawdzano niezwłocznie po ich ukończeniu. Dane z badania zostały wprowadzone dwukrotnie przez dwóch przeszkolonych pracowników, a rozbieżności uzgodniono na podstawie formularzy źródłowych. Brakujące oceny, dzienniki interwencji, rejestry wydawania leków, rejestry przestrzegania zaleceń, dzienniki bezpieczeństwa i formularze wyników były regularnie weryfikowane krzyżowo. Dane zebrane przed wycofaniem się uczestnika lub utratą kontaktu w ramach obserwacji zachowano w zbiorze danych do analizy.
10. Analiza statystyczna
Analiza statystyczna została przeprowadzona przy użyciu oprogramowania statystycznego. Wszyscy zrandomizowani uczestnicy zostali włączeni do analizy zgodności z przydziałem (intention-to-treat), o ile dostępne były dane dotyczące wyników. Przeprowadzono również analizę zgodności z protokołem (per-protocol) dla uczestników, którzy zrealizowali co najmniej 80% zaplanowanych procedur interwencyjnych.
Zmienne ciągłe podsumowano jako średnia ± odchylenie standardowe w przypadku rozkładu zbliżonego do normalnego oraz jako mediana z rozstępem międzykwartylnym, gdy założenia dotyczące rozkładu nie zostały spełnione. Zmienne kategoryczne podsumowano jako liczby i wartości procentowe. Porównania międzygrupowe parametrów wyjściowych przeprowadzono za pomocą testu t dla prób niezależnych lub testu U Manna-Whitneya dla zmiennych ciągłych oraz testu chi-kwadrat lub dokładnego testu Fishera dla zmiennych kategorycznych, zależnie od potrzeb.
Głównym punktem końcowym była różnica między grupami w całkowitym wyniku PSQI w 3. miesiącu. Podłużne wyniki PSQI oraz HADS-A analizowano przy użyciu liniowych modeli efektów mieszanych. Efekty stałe obejmowały grupę, czas oraz interakcję grupy z czasem. Uczestnika uwzględniono jako losowy punkt przecięcia, aby uwzględnić korelację wewnątrzpacjencką pomiędzy powtórzeniami pomiarów. Czas traktowano jako kategoryczny czynnik powtarzany. W przypadku osiągnięcia zbieżności modelu dla pomiarów powtarzanych stosowano nieustrukturyzowaną macierz kowariancji; w przypadku braku zbieżności, w ramach analizy wrażliwości, stosowano autoregresyjną strukturę kowariancji pierwszego rzędu. Modele szacowano przy użyciu ograniczonej maksymalnej wiarygodności. Modele skorygowane obejmowały wiek, płeć, czas trwania dializ oraz pierwotną chorobę nerek jako współzmienne. Założenia modelu oceniano poprzez analizę rozkładów reszt, wykresów wartości dopasowanych względem reszt oraz obserwacji wpływowych.
Analizy osób odpowiadających na leczenie przeprowadzono w każdym punkcie kontrolnym. Odpowiedź w zakresie snu została zdefiniowana jako obniżenie całkowitego wyniku PSQI o co najmniej 3 punkty względem wartości wyjściowej, natomiast odpowiedź w zakresie lęku zdefiniowano jako obniżenie wyniku HADS-A o co najmniej 2 punkty względem wartości wyjściowej. Wskaźniki odpowiedzi przedstawiono w procentach. Dla różnic w odpowiedzi między grupami obliczono absolutne różnice ryzyka oraz 95% przedziały ufności. Liczbę pacjentów koniecznych do leczenia (NNT) obliczono jako odwrotność absolutnej różnicy ryzyka w przypadkach, gdy grupa interwencyjna wykazała wyższy wskaźnik odpowiedzi niż grupa kontrolna; wartości te zinterpretowano opisowo, ponieważ analizy osób odpowiadających na leczenie miały charakter wtórny.
Przyjęto, że brakujące dane dotyczące wyników są brakujące losowo po uwzględnieniu zaobserwowanych informacji wyjściowych oraz z okresu obserwacji. W analizie głównej zastosowano liniowe modele mieszane, wykorzystując wszystkie dostępne dane z pomiarów powtarzanych, bez wykluczania uczestników z niepełną obserwacją. Przeprowadzono analizę wrażliwości w modelu intention-to-treat z wykorzystaniem wielokrotnej imputacji metodą pełnej specyfikacji warunkowej. Wygenerowano dwadzieścia imputowanych zbiorów danych, wykorzystując jako predyktory grupę, wiek, płeć, czas trwania dializowania, pierwotną chorobę nerek, wyjściowy wynik PSQI, wyjściowy wynik HADS-A oraz dostępne wyniki PSQI i HADS-A z okresu obserwacji. Estymatory ze zbiorów imputowanych połączono zgodnie ze standardowymi procedurami wielokrotnej imputacji.
Analizy w podgrupach przeprowadzono z uwzględnieniem wieku, płci, stażu dializoterapii, pierwotnej choroby nerek, wyjściowego stopnia nasilenia PSQI oraz wyjściowego stopnia nasilenia HADS-A. Wyjściowy stopień nasilenia PSQI podzielono na grupy 8–14 oraz ≥ 15, a wyjściowy stopień nasilenia HADS-A na grupy < 11 oraz ≥ 11. Analizy w podgrupach miały charakter eksploracyjny. Przy opisie efektów w podgrupach badano interakcje typu grupa-czas-podgrupa, zamiast polegać wyłącznie na wartościach P w obrębie poszczególnych podgrup.
Analizy wrażliwości obejmowały analizę per-protocol, analizę intention-to-treat z zastosowaniem wielokrotnej imputacji, dodatkowo skorygowany model uwzględniający bazowe pomiary zmęczenia (tam, gdzie były dostępne) oraz analizę z wyłączeniem uczestników, których bazowe wartości PSQI mieściły się w najwyższych 5% rozkładu próby. Analizy trajektorii i sieci przeprowadzono jako analizy eksploracyjne i nie zostały one wykorzystane do zdefiniowania głównego efektu terapeutycznego.
Wszystkie testy statystyczne były obustronne. Wartość p < 0,05 uznano za statystycznie istotną dla głównego punktu końcowego. Analizy wtórne i eksploracyjne nie były korygowane pod kątem wielokrotnych porównań i zostały zinterpretowane z ostrożnością.