המחקר נערך בהתאם להצהרת הלסינקי, והפרוטוקול אושר על ידי ועדת האתיקה של בית החולים השני המסונף לאוניברסיטת ננג'ינג הרפואית (מספר אישור: 2021-KY-100-01). הסכמה מדעת בכתב התקבלה מכל המשתתפים.
מחקר תצפיתי רטרוספקטיבי זה, שנערך במרכז אחד, בוצע במחלקה לנוירולוגיה בבית החולים המשלי שני של אוניברסיטת נאנג'ינג לרפואה. מינואר 2022 עד דצמבר 2023, נכללו בסך הכל 32 חולים עם PD ראשוני שעמדו בקריטריונים האבחנתיים למחלת הפרקינסון (מהדורת 2020). בסופו של דבר, שישה חולים שעמדו בקריטריוני בחירה קפדניים והשלימו את כל הליכי המחקר נכללו בניתוח המשך. פרטי הסינון, ההחרגה והמעקב מוצגים ב-איור משלים 1.
איסוף דגימות צואה וניתוח הקהילה המיקרוביאלית
מיקרוביוטה עיתית הושגה ממתנדבים בריאים. תהליך סינון התורמים בוצע בקפידה בהתאם להסכמת המומחים הסינית בנושא יישום קליני של השתלת מיקרוביוטה עיתית. קריטריונים להכללת תורמים: גיל 18–40 שנים, ללא היסטוריה של מחלות כרוניות; ללא היסטוריה של מחלות במערכת העיכול, הפרעות מטבוליות, מחלות אוטואימוניות או הפרעות פסיכיאטריות; ללא שימוש באנטיביוטיקה, פרוביוטיקה, תרופות מדכאות חיסון או מעכבי משאבת פרוטונים ב-3 החודשים האחרונים; ללא נסיעות לאזורי מגפה, ללא עירויי דם או היסטוריה כירורגית ב-6 החודשים האחרונים. סינון דם ועיתים לשלילת מחלות זיהומיות ופתוגנים כולל: נגיף הכשל החיסוני האנושי (HIV), נגיף הדלקת הכבד B (HBV), נגיף הדלקת הכבד C (HCV), Treponema pallidum, Clostridium difficile, סלמונלה (Salmonella), Shigella, קמפילובקטר (Campylobacter), פתוגני Escherichia coli, נגיפי מעי וטפילים. נפחה של דגימת צואה תקפה אחת צריך להיות לפחות 50 גרם. באמצעות מערכת אוטומטית לטיהור מיקרוביוטה של צואה (ראה טבלת חומרים), 500 מ"ל של תמיסל מלח סטרילי הוספו לכל 100 גרם של צואה טרייה, שלושה סבבים של סנטריפוגציה ושטיפה סטנדרטיים (700 × ג', 3 דקות), ומיד לאחר מכן בוצע חיסון לאחר הכנת תרחיף עם תמיסת מלח סטרילית ביחס נפחים של 1:2. הנבדקים קיבלו הזרקה בודדת של 150 מ"ל תרחיף המכיל כ-5 × 10^13 חיידקים חיים13.
כל המקבלים קיבלו FMT באמצעות קולונוסקופיה, כאשר המכשיר הוכנס עד לעיוור (cecum) ולמעי הגס העולה; התמיה הוזרקה לאט ובאופן אחיד דרך תעלת הביופסיה כדי לעודד יישוב מועדף של המעי הגס הימני. לאחר ההזרקה, נשארו המטופלים בתנוחה שכיבה על צדם הימני למשך 2 h לפחות. לפני הפרוצדורה, נשמרה דיאטה דלה בשיירי מזון למשך 3 ימים, ולאחר מכן דיאטה נוזלית מלאה למשך 24 h לפני הפרוצדורה, עם צום ממוצקים למשך 6 h ומנוזלים למשך 2 h לפני ההזרקה. בערב שלפני הש transplantation, בוצע ניקוי סטנדרטי של כל המעי באמצעות תמיסת אלקטרוליטים של polyethylene glycol עד לקבלת הפרשה צלולה ללא צואה. הגבלות תרופתיות כללו הפסקה של אנטיביוטיקה סיסטמית למשך 72 h לפחות לפני ה-FMT והימנעות ממנה למשך שבוע לאחר מכן. מעכבי פומפת פרוטון (PPIs) הופסקו 24–48 h לפני הפרוצדורה, וחומרים משלשליים ופרוביוטיקה הופסקו ביום ה-transplantation.
במהלך הפרוצדורה, בוצע ניטור רציף של קצב הלב, לחץ הדם ורוויית החמצן, ונבחנו תסמינים כגון נפיחות בבטן, כאב או בחילה; העירוי הופסק באופן מיידי במקרה של כאב בטן חריף או תסמונת של hypotension. במהלך 7 הימים הראשונים לאחר ההשתלה, נערך מעקב יומי לתיעוד תדירות יציאות, עקביות הצואה וכל תופעות הלוואי, כולל כאבי בטן, שלשולים, חום והמטוכזיה (hematochezia). המעקב ארוך הטווח נמשך עד 5 חודשים, וכלל דגימות צואה מתוכננות לניתוח מיקרוביוטת המעי והערכות קליניות באמצעות מדד UPDRS (Unified Parkinson’s Disease Rating Scale) ומדד Wexner Constipation Score, לצד תיעוד של שינויים בתסמינים ושימוש בתרופות נלוות, כדי להעריך באופן מקיף הן את היעילות והן את הבטיחות.
דגימות צואה נאספו בשש נקודות זמן: לפני FMT ובחודשים 1, 2, 3, 4 ו-5 לאחר FMT. כל המשתתפים קיבלו הנחיות סטנדרטיות לאיסוף דגימות צואה. כ-10–20 g של חלק מרכזי מדגימת הצואה נאספו באמצעות פינצטה סטרילית והוכנסו למכלי איסוף סטריליים. הדגימות הועברו למעבדה תוך 2–3 h לאחר האיסוף וחולקו לאליקוטיים בתוך מבחנות מיקרוצנטריפוגה סטריליות של 1.5 mL (ראו טבלת חומרים) על ידי חוקרים מיומנים תחת תנאים אספטיים. כל האליקוטיים הוקפאו מיד בטמפרטורה של −80 °C עד לניתוח נוסף.
לאחר איסוף דגימות הצואה, הובלו כל הדגימות על קרח יבש לספק רצף מסחרי (ראו טבלת החומרים) לצורך רצף של הגן 16S rRNA. הכנת הספרייה, הרצף ובקרת האיכות הראשונית בוצעו בהתאם לפרוטוקול הסטנדרטי של הספק. DNA גנומי הופק מכל דגימת צואה, והאזורים ההיפרו-וריאביליים V3–V4 של הגן הבקטריאלי 16S rRNA הוגברו באמצעות תגובת שרשרת פולימראז (PCR) תוך שימוש בפריימרים 341F ו-805R. הרצף בוצע בפלטפורמת רצף קצר שולחנית (ראו טבלת החומרים) תוך שימוש בכימיה של קצוות זוגיים 2 × 300 bp כדי לאפיין את פרופילי המיקרוביום במעי של המטופלים לפני FMT ובנקודות זמן שונות במעקב לאחר ה-FMT. לאחר הרצף, עובדו הקריאות הגולמיות כדי לקבל רצפים באיכות גבוהה לניתוח מיקרוביאלי המשכי. סיווג טקסונומי בוצע באמצעות חבילת השמה טקסונומית מסחרית של 16S ומסד הנתונים למזהה מיקרוביאלי הנלווה אליה (ראו טבלת החומרים). הפרופילים הטקסונומיים שהתקבלו שימשו לניתוחים עוקבים של גיוון מיקרוביאלי, הרכב קהילה, שפע דיפרנציאלי ותפקוד מיקרוביאלי חזוי.
הערכות קליניות ותוויות אנליטיות
חולים עם PD סווגו לשתי קבוצות חומרה בהתאם למשך המחלה: קבוצת חומרה גבוהה (>6 שנים, סומנה כ-"_H") וקבוצת חומרה נמוכה (≤6 שנים, סומנה כ-"_L"). דגימות צואה שנאספו לפני FMT הוגדרו כדגימות בסיס וסומנו כ-FB_H או FB_L בהתאם לקבוצת החומרה המתאימה. דגימות שנאספו בחודשים 1, 2, 3, 4 ו-5 לאחר ה-FMT סומנו כ-FP1M_H, FP2M_H, FP3M_H, FP4M_H ו-FP5M_H עבור קבוצת החומרה הגבוהה, וכ-FP1M_L, FP2M_L, FP3M_L, FP4M_L ו-FP5M_L עבור קבוצת החומרה הנמוכה. הפרמטרים הקליניים שהוערכו כללו את מדד העצירות של Wexner, מדד תסמינים לא-מוטוריים (NMSS), שאלון תסמינים לא-מוטוריים (NMSQ), סקר נוירופסיכיאטרי (NPI), שלב Hoehn–Yahr, מדד פעילויות היומיום (ADL), מדד איכות השינה של פיטסבורג (PSQI), ה-(UPDRS I–VI) המקורי, וציון UPDRS כולל. שילוב של שיטות סטטיסטיות תיאוריות והסקיות שימש לניתוח השינויים במדדים הקליניים של חולי PD בקבוצת החומרה הגבוהה (קבוצת _H) ובקבוצת החומרה הנמוכה (קבוצת _L) בשש נקודות זמן. כדי לבחון באופן ראשוני הבדלים בין הקבוצות בכל נקודות הזמן, הושוו התפלגויות המדדים בין קבוצות _H ו- _L בכל נקודת זמן באמצעות מבחן Wilcoxon rank-sum, וכן הוחל מודל ליניארי מעורב (LMM) קפדני יותר עם חסימות אקראיות ספציפיות למטופל. בשל גודל המדגם המוגבל, תוצאות המודל עם ההשפעות המעורבות מוצגות כמרחבי רפרונסיים גששיים בלבד. נקודת הקצה הראשית של מחקר גשש זה הייתה השינוי במגוון ובהרכב המיקרוביום במעי מהנקודת הבסיס לנקודות הזמן שלאחר השתלת מיקרוביוטת צואה, כאשר שינויים קליניים הוערכו כנקודות קצה משניות.
הקבצת OTU, ניתוח גיוון ושפע טקסונומי
יחידות טקסונומיות תפעוליות (OTUs) שימשו לסיכום הרכב הקהילה המיקרובסית, ורצפים עם דמיון של ≥97% קובצו לאותו OTU. עומק הרצף הוערך באמצעות עקומות ררפיקציה (rarefaction) מבוססות Shannon, כאשר מישורי העקומה העידו על כיסוי מספיק באופן כללי לצורך ניתוחי שונות במורד הזרם.
מדדי גיוון אלפא (Alpha diversity) הוערכו באמצעות מדדי Chao1 ו-abundance-based coverage estimator (ACE) עבור עושר מיקרוביאלי, ומדדי Shannon ו-Simpson עבור גיוון כללי ומידה של אחידות (evenness). גיוון בטא (Beta diversity) הוערך באמצעות מרחקי UniFrac משוקללים ולא משוקללים, ולאחריהם ניתוח רכיבים ראשיים (PCoA) להדמיית הבדלים בהרכב הקהילה המיקרוביאלית בין הקבוצות. הבדלים בין-קבוצתיים בגיוון אלפא ובגיוון בטא נבחנו באמצעות מבחנים לא-פרמטריים וניתוח שונות רב-משתני permutational (PERMANOVA), בהתאמה, כפי שמתואר בסעיף הניתוח הסטטיסטי.
ההרכב הטקסונומי סוכם ברמות המערכה, המחלקה, הסדרה, המשפחה והסוג, בהתבסס על נתוני ריצוף של גן 16S rRNA. באמצעות חבילת ניתוח אקולוגיה של קהילות עבור R (ראו טבלת חומרים), זהה והצג באופן ויזואלי את 10 הטקסונים השכיחים ביותר בכל רמה טקסונומית. פרופילים טקסונומיים הוצגו באמצעות תרשימי עמודות, וטקסונים נבחרים הושוו בהמשך בין קבוצות באמצעות תרשימי קופסה.
ניתוח שפע דיפרנציאלי של טקסים מיקרוביאליים
טקסונים בעלי שפע שונה בין קבוצות שהוגדרו מראש זוהו באמצעות ניתוח גודל אפקט של אפלימננט ליניארי (LEfSe). LEfSe משלב את מבחן Kruskal–Wallis, מבחן Wilcoxon וניתוח אפלימננט ליניארי (LDA) כדי לזהות טקסונים בעלי מובהקות סטטיסטית וערך מבחין. טקסונים עם P < 0.05 וציון LDA (log10) > 3 נחשבו כבעלי שפע שונה. השפע היחסי של הטקסונים המבחינים הוצג ויזואלית באמצעות תרשימי עמודות טקסונומיים.
ניבוי תפקודי וניתוח קורלציה קלינית
פרופילים תפקודיים חזויים של מיקרוביוטת המעי הוסקו מנתוני רצף של גן 16S rRNA באמצעות Tax4Fun2 ואנוטציה בוצעה על פי מאגר הנתונים Kyoto Encyclopedia of Genes and Genomes (KEGG). מסלולים חזויים שהועשרו באופן מבדל זוהו באמצעות LEfSe, כאשר נקבעו ספי ערכים של P < 0.05 וציון LDA (log10) > 3. מכיוון שפרופילים תפקודיים אלו הוסקו ולא נמדדו באופן ישיר בגישות מטגנומיות או מטבולימיות, הם פורשו כהבדלים תפקודיים חזויים. ניתוח מתאם דרג של Spearman בוצע כדי לבחון קשרים בין מאפיינים מיקרוביאליים למדדים קליניים.
ניתוח סטטיסטי
ניתוחים סטטיסטיים בוצעו באמצעות תוכנת R (ראו טבלת חומרים). משתנים קליניים סוכמו כחציונים וטווחים (מינימום–מקסימום). מדדי alpha diversity הושוו בין קבוצות באמצעות מבחן Wilcoxon rank-sum, בעוד ש-beta diversity הוערך באמצעות PERMANOVA עם 999 פרמוטציות. טקסלומיות בעלות שפע דיפרנציאלי ומסלולים פונקציונליים חזויים זוהו באמצעות LEfSe, כאשר P < 0.05 וציון LDA (log10) > 3 שימשו כספי השמעות. עבור ניתוח שפע דיפרנציאלי לרוחב זמן, הופעל מודל אפקטים מעורבים (mixed-effects model) עם חיתוכי מקרי (random intercepts) ספציפיים למטופל, ולאחריו בוצעו השוואות זוגיות post-hoc עם תיקון Benjamini-Hochberg לקצב הגילוי השגוי (false discovery rate). ניתוח מתאם דרוג של Spearman שימש להערכת קשרים בין מאפייני מיקרוביום המעי למדדים קליניים, תוך יישום תיקון Benjamini–Hochberg לקצב הגילוי השגוי עבור ערכי ה-P של המתאם. ערך P דו-זנבי < 0.05 נחשב למשמעותי סטטיסטית. תיקון להשוואות מרובות הוחל על ניתוחי המתאמים ושפע הטקסלומיות. בהינתן גודל המדגם הקטן ומבנה הדגימה הlongitudinal החוזר, הניתוחים הסטטיסטיים נחשבו כחקרים גישושים, והתוצאות פורשו בזהירות.