Badanie przeprowadzono zgodnie z Deklaracją Helsińską, a protokół został zatwierdzony przez Komisję Etyki Drugiego Szpitala Powiązanego z Uniwersytetem Medycznym w Nankinie (numer zatwierdzenia: 2021-KY-100-01). Od wszystkich uczestników uzyskano pisemną świadomą zgodę.
To retrospektywne obserwacyjne badanie jednomiestresowe zostało przeprowadzone w Katedrze Neurologii w Drugim Szpitalu Afiliacyjnym Uniwersytetu Medycznego w Nanjing. W okresie od stycznia 2022 do grudnia 2023 roku do badania włączono łącznie 32 pacjentów z pierwotną chorobą Parkinsona (PD), którzy spełniali kryteria diagnostyczne choroby Parkinsona (wydanie z 2020 r.). Ostatecznie do dalszej analizy zakwalifikowano sześciu pacjentów, którzy spełnili rygorystyczne kryteria wyboru i przeszli wszystkie procedury badawcze. Szczegóły dotyczące przesiewu, wykluczeń oraz obserwacji przedstawiono na Rysunku uzupełniającym 1.
Pobieranie próbek kału i analiza społeczności mikrobiologicznej
Mikrobiota kałowa została pozyskana od zdrowych ochotników. Proces przesiewowy dawców przeprowadzono ściśle zgodnie z chińskim konsensusem eksperckim dotyczącym klinicznego zastosowania transplantacji mikrobioty kałowej. Kryteria włączenia dawców: wiek 18–40 lat, brak historii chorób przewlekłych; brak historii chorób przewodu pokarmowego, zaburzeń metabolicznych, chorób autoimmunologicznych lub zaburzeń psychiatrycznych; brak stosowania antybiotyków, probiotyków, immunosupresantów lub inhibitorów pompy protonowej w ciągu ostatnich 3 miesięcy; brak podróży do obszarów epidemicznych, przetoczeń krwi lub historii operacji chirurgicznych w ciągu ostatnich 6 miesięcy. Badania przesiewowe krwi i kału w celu wykluczenia chorób zakaźnych i patogenów obejmowały: wirus niedoboru odporności u ludzi (HIV), wirus zapalenia wątroby B (HBV), wirus zapalenia wątroby C (HCV), Treponema pallidum, Clostridium difficile, Salmonella, Shigella, Campylobacter, patogenne Escherichia coli, wirusy jelitowe oraz pasożyty. Objętość pojedynczej prawidłowej próbki kału powinna wynosić co najmniej 50 g. Przy użyciu zautomatyzowanego systemu oczyszczania mikrobioty kałowej (patrz Tabela materiałów), do każdych 100 g świeżego kału dodano 500 mL jałowej soli fizjologicznej, przeprowadzono trzy serie standardowego wirowania i przemywania (700 × g, 3 min), a następnie niezwłocznie wykonano inokulację po przygotowaniu zawiesiny z jałową solą fizjologiczną w stosunku objętości 1:2. Obiekty otrzymały pojedynczą iniekcję 150 mL zawiesiny zawierającej około 5 × 10^13 żywych bakterii13.
Wszyscy biorcy otrzymali FMT drogą kolonoskopową, z wprowadzeniem endoskopu do kątnicy i okrężnicy wstępującej; zawiesinę wprowadzano powoli i równomiernie przez kanał biopsyjny, aby wspomóc preferencyjną kolonizację prawej części okrężnicy. Po infuzji pacjenci pozostawali w bocznej pozycji prawej przez co najmniej 2 h. Przed zabiegiem przez 3 dni stosowano dietę niskoresztkową, a następnie przez 24 h przed procedurą dietę płynną, z zachowaniem postu od pokarmów stałych przez 6 h i od płynów przez 2 h przed infuzją. Wieczorem w dniu poprzedzającym transplantację przeprowadzono standardowe oczyszczanie całego jelita za pomocą elektrolitowego roztworu glikolu polietylenowego do momentu uzyskania przejrzystego, wolnego od treści kałowej wydzieliny. Ograniczenia w stosowaniu leków obejmowały odstawienie antybiotyków systemowych na co najmniej 72 h przed FMT oraz unikanie ich stosowania przez 1 tydzień po zabiegu. Inhibitory pompy protonowej (PPI) odstawiano na 24–48 h przed procedurą, a środki przeczyszczające i probiotyki zawieszano w dniu transplantacji.
W trakcie procedury stale monitorowano tętno, ciśnienie tętnicze oraz saturację tlenem, a także obserwowano wszelkie objawy, takie jak wzdęcia brzucha, ból lub nudności; w przypadku wystąpienia silnego bólu brzucha lub hipotonii natychmiast przerywano infuzję. Przez pierwsze 7 dni po transplantacji przeprowadzano codzienną obserwację w celu zarejestrowania częstotliwości wypróżnień, konsystencji stolca oraz wszystkich zdarzeń niepożądanych, w tym bólu brzucha, biegunki, gorączki i hematochezji. Długoterminowa obserwacja trwała do 5 miesięcy i obejmowała planowe pobieranie próbek stolca do analizy mikrobioty jelitowej oraz oceny kliniczne z wykorzystaniem skali UPDRS (Unified Parkinson’s Disease Rating Scale) i skali Wexner Constipation Score, wraz z dokumentacją zmian objawów i stosowania leków towarzyszących, aby kompleksowo ocenić zarówno skuteczność, jak i bezpieczeństwo.
Próbki kału pobierano w sześciu punktach czasowych: przed FMT oraz w 1., 2., 3., 4. i 5. miesiącu po FMT. Wszyscy uczestnicy otrzymali ustandaryzowane instrukcje dotyczące pobierania próbek kału. Przy użyciu sterylnej pęsety pobrano około 10–20 g środkowej części stolca i umieszczono ją w sterylnych pojemnikach do próbek. Próbki transportowano do laboratorium w ciągu 2–3 h od pobrania, a następnie przeszkoleni badacze w warunkach aseptycznych podzielili je na aliquoty do sterylnych probówek do mikrocentryfugi o pojemności 1,5 mL (patrz Tabela materiałów). Wszystkie aliquoty natychmiast przechowywano w temperaturze −80 °C do czasu dalszej analizy.
Po pobraniu próbek kału wszystkie materiały zostały przetransportowane w suchem lodzie do komercyjnego dostawcy usług sekwencjonowania (patrz Tabela materiałów) w celu przeprowadzenia sekwencjonowania genu 16S rRNA. Przygotowanie biblioteki, sekwencjonowanie oraz wstępna kontrola jakości zostały wykonane zgodnie ze standardowym protokołem dostawcy. Z każdej próbki kału wyekstrahowano genomowe DNA, a następnie za pomocą reakcji łańcuchowej polimerazy (PCR) z użyciem starterów 341F i 805R przeprowadzono amplifikację hiperzmiennych regionów V3–V4 bakteryjnego genu 16S rRNA. Sekwencjonowanie wykonano na laboratoryjnej platformie do sekwencjonowania krótkich odczytów (patrz Tabela materiałów), stosując chemię dla odczytów parzystych 2 × 300 bp, aby scharakteryzować profile mikrobioty jelitowej pacjentów przed FMT oraz w różnych punktach czasowych kontroli po FMT. Po sekwencjonowaniu surowe odczyty zostały przetworzone w celu uzyskania wysokiej jakości sekwencji do dalszej analizy mikrobiologicznej. Klasyfikację taksonomiczną przeprowadzono przy użyciu komercyjnego pakietu do przypisywania taksonomii 16S oraz powiązanej z nim bazy danych identyfikacji mikrobiologicznej (patrz Tabela materiałów). Uzyskane profile taksonomiczne wykorzystano w późniejszych analizach różnorodności mikrobiologicznej, składu społeczności, różnicowej obfitości oraz przewidywanej funkcji mikrobiologicznej.
Oceny kliniczne i etykiety analityczne
Pacjentów z PD podzielono na dwie grupy w zależności od stopnia zaawansowania choroby, określonego przez czas trwania schorzenia: grupę o wysokim stopniu zaawansowania (>6 lat, oznaczona jako „_H”) oraz grupę o niskim stopniu zaawansowania (≤6 lat, oznaczona jako „_L”). Próbki kału pobrane przed FMT zdefiniowano jako próbki wyjściowe i oznaczono jako FB_H lub FB_L zgodnie z odpowiednią grupą zaawansowania. Próbki pobrane w 1., 2., 3., 4. i 5. miesiącu po FMT oznaczono jako FP1M_H, FP2M_H, FP3M_H, FP4M_H i FP5M_H dla grupy o wysokim stopniu zaawansowania oraz FP1M_L, FP2M_L, FP3M_L, FP4M_L i FP5M_L dla grupy o niskim stopniu zaawansowania. Oceniane parametry kliniczne obejmowały skalę zaparć Wexnera, skalę objawów pozaruchowych (NMSS), kwestionariusz objawów pozaruchowych (NMSQ), inwentarz neuropsychiatryczny (NPI), stopień w skali Hoehn–Yahr, skalę czynności życia codziennego (ADL), indeks jakości snu w Pittsburghu (PSQI), oryginalną skalę (UPDRS I–VI) oraz całkowity wynik UPDRS. Do analizy zmian wskaźników klinicznych pacjentów z PD w grupie o wysokim stopniu zaawansowania (grupa _H) i grupie o niskim stopniu zaawansowania (grupa _L) w sześciu punktach czasowych zastosowano kombinację metod statystyki opisowej i indukcyjnej. Aby wstępnie zbadać różnice między grupami we wszystkich punktach czasowych, rozkłady wskaźników między grupami _H i _L porównywano w każdym punkcie czasowym za pomocą testu suma rang Wilcoxona, a następnie zastosowano bardziej rygorystyczny liniowy model mieszany (LMM) z losowymi punktami przecięcia dla poszczególnych pacjentów. Ze względu na ograniczoną liczebność próby, wyniki modelu efektów mieszanych przedstawiono jedynie jako referencje eksploracyjne. Głównym punktem końcowym tego badania eksploracyjnego była zmiana różnorodności i składu mikrobioty jelitowej od stanu wyjściowego do punktów czasowych po przeszczepie mikrobioty kałowej, natomiast zmiany kliniczne oceniano jako wtórne punkty końcowe.
Klastrowanie OTU, analiza różnorodności oraz obfitość taksonomiczna
Do podsumowania składu społeczności mikrobiologicznej wykorzystano jednostki taksonomiczne operacyjne (OTU), a sekwencje o podobieństwie ≥97% pogrupowano w ramach tego samego OTU. Głębokość sekwencjonowania oceniono za pomocą krzywych rarefakcji opartych na indeksie Shannona, przy czym plateau krzywych wskazywało na ogólnie wystarczający stopień pokrycia dla dalszych analiz różnorodności.
Różnorodność alfa oceniano za pomocą indeksów Chao1 oraz ACE (abundance-based coverage estimator) w celu określenia bogactwa mikrobiologicznego, a także indeksów Shannona i Simpsona w celu oceny ogólnej różnorodności i wyrównania. Różnorodność beta ewaluowano przy użyciu ważonych i nieważonych odległości UniFrac, a następnie przeprowadzono analizę głównych współrzędnych (PCoA), aby zwizualizować różnice w składzie społeczności mikrobiologicznej pomiędzy grupami. Różnice między grupami w zakresie różnorodności alfa i beta oceniano odpowiednio za pomocą testów nieparametrycznych oraz permutacyjnej wielowymiarowej analizy wariancji (PERMANOVA), zgodnie z opisem w sekcji analizy statystycznej.
Skład taksonomiczny podsumowano na poziomach typu, klasy, rzędu, rodziny i rodzaju w oparciu o dane z sekwencjonowania genu 16S rRNA. Używając pakietu do analizy ekologii społeczności w programie R (patrz Tabela materiałów), zidentyfikowano i zwizualizowano 10 najliczniejszych taksonów na każdym poziomie taksonomicznym. Profile taksonomiczne przedstawiono za pomocą wykresów słupkowych, a wybrane taksony porównano dodatkowo między grupami przy użyciu wykresów pudełkowych.
Analiza różnic w obfitości taksonów mikrobiologicznych
Taksony o różnej liczności w obrębie zdefiniowanych grup zidentyfikowano za pomocą analizy wielkości efektu dyskryminanta liniowej (LEfSe). LEfSe integruje test Kruskala–Wallisa, test Wilcoxona oraz analizę dyskryminantową liniową (LDA), aby zidentyfikować taksony wykazujące zarówno istotność statystyczną, jak i wartość dyskryminacyjną. Za taksony o różnej liczności uznano te, dla których P < 0,05 oraz wynik LDA (log10) > 3. Względną liczność taksonów dyskryminacyjnych zwizualizowano za pomocą wykresów słupkowych taksonomicznych.
Przewidywanie funkcjonalne i analiza korelacji klinicznej
Przewidywane profile funkcjonalne mikrobioty jelitowej wywnioskowano z danych sekwencjonowania genu 16S rRNA przy użyciu Tax4Fun2 i zanoatowano zgodnie z bazą danych Kyoto Encyclopedia of Genes and Genomes (KEGG). Różnicowo wzbogacone przewidywane szlaki zidentyfikowano za pomocą LEfSe, przyjmując progi P < 0,05 oraz wynik LDA (log10) > 3. Ponieważ profile funkcjonalne te zostały wywnioskowane, a nie zmierzone bezpośrednio metodami metagenomicznymi lub metabolomicznymi, interpretowano je jako przewidywane różnice funkcjonalne. Przeprowadzono analizę korelacji rang Spearmana w celu zbadania powiązań między cechami mikrobiologicznymi a wskaźnikami klinicznymi.
Analiza statystyczna
Analizy statystyczne przeprowadzono przy użyciu oprogramowania R (patrz Tabela Materiałów). Zmienne kliniczne przedstawiono jako mediany i zakresy (minimum–maksimum). Wskaźniki alfa-różnorodności porównano pomiędzy grupami za pomocą testu sum rang Wilcoxona, natomiast beta-różnorodność oceniono za pomocą analizy PERMANOVA z 999 permutacjami. Taksony o różnej obfitości oraz przewidywane ścieżki funkcjonalne zidentyfikowano za pomocą analizy LEfSe, przyjmując progi istotności P < 0,05 oraz wynik LDA (log10) > 3. Do podłużnej analizy różnic w obfitości zastosowano model mieszany z losowymi punktami przecięcia specyficznymi dla pacjenta, a następnie przeprowadzono post-hoc porównania parzyste z korektą stopy odkryć fałszywych Benjamini-Hochberg. Analizę korelacji rang Spearmana wykorzystano do oceny powiązań między cechami mikrobioty jelitowej a wskaźnikami klinicznymi, stosując korektę stopy odkryć fałszywych Benjamini–Hochberg dla wartości P korelacji. Za statystycznie istotne uznano dwustronne P < 0,05. Korektę porównań wielokrotnych zastosowano do analiz korelacji oraz obfitości taksonomicznej. Ze względu na niewielką wielkość próby i projekt powtarzalnego pobierania próbek wzdłużnych, analizy statystyczne uznano za eksploracyjne, a wyniki interpretowano z ostrożnością.