פרטים על כל בעלי החיים, התזונה, המגבים, הערכות, הכלים, התוכנה וספקי השירות המשמשים בפרוטוקול זה מסופקים בטבלת החומרים. כל המחקרים על בעלי חיים בוצעו בהתאם להנחיות המכון הלאומי לבריאות (NIH) לטיפול ושימוש בבעלי חיים במעבדה, וקיבלו אישור מוועדת בעלי החיים הניסיוניים של בית החולים הפרובינציאלי שאנדונג (אישור אתי מס' 2022–007, אושר ב-24 בינואר 2022).
זהירות: דם בעלי חיים, דגימות צואה ורקמות יש לטפל בהם כחומרים שעלולים להיות מסוכנים ביולוגית. כל ההליכים הכוללים בעלי חיים חיים, דגימות ביולוגיות, חדים, ממסים אורגניים וראגנטים לבדיקה כימית צריכים להתבצע בהתאם לבטיחות ביולוגית, טיפול בבעלי חיים ותקנות בטיחות כימית מוסדיות. גופות בעלי חיים, רקמות, חומרים מזוהמים בדם ודגימות צואה יש לאסוף במיכלי פסולת ביולוגית ייעודיים ולהיפנה דרך מוסד לבעלי חיים או משרד בטיחות ביולוגית. מחטים, מחטי גוואג', צינורות נימיים וחומרים חדים אחרים יש להשליך מיד לאחר השימוש במיכלים מאושרים. פסולת ממס אורגנית, כולל מתנול, אצטוניטריל ותמיסות המכילות חומצה פורמית, צריכה להיאסף במיכלי פסולת כימית מסומנים ולהיפנה דרך תוכנית הפסולת הכימית המסוכנת המוסדית. יש לאסוף ולזרוק פסולת מבדיקות עומס ביוכימיות וחמצוניות בהתאם לגיליונות הבטיחות של היצרן ולכללי הבטיחות הכימיים המוסדיים.
בעלי חיים ניסיוניים
עכברים זכרים בני שישה שבועות מסוג C57BL/6J (n = 4 לכל קבוצה) נשמרו בתנאים מסוימים ללא פתוגנים (SPF) עם מחזור אור/חושך של 12 שעות בטמפרטורה של 22 ± 2 מעלות צלזיוס ולחות יחסית של 50–60%. לבעלי חיים הייתה גישה חופשית למזון ולמים, כאשר לכל הכלוב היו מקסימום חמישה עכברים. עכברי הביקורת קיבלו תזונה סטנדרטית שכללה 70% פחמימות, 20% חלבון ו-10% שומן (קלוריות). כדי להקים את מודל העכבר MASLD, העכברים קיבלו תזונה עתירת שומן במשך 8 שבועות שכללה 20% פחמימות, 20% חלבון ו-60% שומן (ק"ל). כל הדיאטות אוחסנו בטמפרטורה של 4 מעלות צלזיוס וחודשו פעמיים בשבוע. העכברים חולקו באקראי לקבוצות ניסוי באמצעות מחולל מספרים אקראיים, וגודל המדגם נקבע (n = 4 לכל קבוצה) בהתבסס על מחקרים קודמים שהשתמשו במודלים דומים של עכברים מבוססי תזונה עתירת שומן וניסויים ראשוניים. כל הטיפולים ניתנו באותו זמן בכל יום כדי למזער את השונות הצירקדית.
שמן המיקרואצות הופק וטוהר על ידי האקדמיה למדעי החקלאות של שאנדונג, עם טוהר של כ-99%. שמן האצות המיקרו אוחסן בטמפרטורה של −30 מעלות צלזיוס, מוגן מאור. לפני ההזנה, השמן הוכן עם אמולסיפייר וערבב היטב כדי להבטיח הומוגניות. במיוחד, שמן מיקרואצות עשיר ב-DHA עבר אמולסיפיקציה בנתרן קרבוקסי מתילצלולוז 0.5% סטרילי המכיל 0.5% טווין-80 במי מלח רגילים. שמן האצות המיקרו ורכב האמולסיפייר התערבבו ביחס של 1:9 (v/v) והועברו לוורטקס למשך 2 דקות מיד לפני ההערכה כדי לקבל מתלה אחיד. עכברים בקבוצת המיקרואצות שטופלו בשמן קיבלו את התכשיר פעם ביום באמצעות מחט האכלה סטרילית במינון של 10 מיקרוליטר/גרם משקל גוף. נפח הגוואג' הותאם מדי יום בהתאם למשקל הגוף כדי להבטיח עקביות במינון. עכברים בקבוצת הביקורת קיבלו נפח דומה של מי מלח סטריליים רגילים באותם תנאים. המלח ששימש היה הזרקת כלוריד נתרן של 0.9%.
טיפול פרוביוטי כלל נוסחה מורכבת שכללה Lactobacillus plantarum (DY-1), Lactobacillus acidophilus (KDB-03), Lactobacillus casei (KDB-LC) ו-Bacillus coagulans (GIM 1.645). ארבעת הזנים התערבבו ביחס תאים חיוניים שווה של 1:1:1:1. כל זן תרם 2.5 × 107 CFU/mL לנוסחה הסופית, מה שהניב ריכוז חיידקים כולל של 1 × 108 CFU/mL. תערובת הפרוביוטיקה סונתזה על ידי האקדמיה למדעי החקלאות של שאנדונג ואינה זמינה מסחרית. טיפול פרוביוטי כלל מרק תסיסה פרוביוטי מורכב המכיל Lactobacillus plantarum DY-1, Lactobacillus acidophilus KDB-03, Lactobacillus casei KDB-LC, ו-Bacillus coagulans GIM 1.645. ריכוז החיידקים של הנוסחה היה 1 ×10 8 CFU/mL. מרק התסיסה הפרוביוטי נשמר בטמפרטורה של 4 מעלות צלזיוס לפני השימוש. עכברים בקבוצת הטיפול בפרוביוטיקה קיבלו את נוסחת הפרוביוטיקה פעם ביום באמצעות גבדה אורלית במינון של 0.01 מ"ל לגרם משקל גוף לאורך תקופת הניסוי.
הזנים הבודדים גודלו בתנאים מתאימים ונקטפו בשלב הגידול הבינוני (OD₆₀₀₀ ≈ 0.6), ואז שולבו ביחס המצוין. תרביות חיידקיות צונטרו כדי להסיר את הסופרנאנט, והכדורים הושעו מחדש בתמלח סטרילי עם בופר פוספט (PBS). המתלה החיידקית הטרייה נשמרה בקפאה לפני המתן כדי לשמור על החיוניות. הריכוז הסופי של החיידקים הותאם ל-1 ×10 8 CFU/mL לפני ההעברה. עכברים בקבוצת הפרוביוטיקה קיבלו את הנוסחה באמצעות גאווג' דרך הפה במינון של 0.01 מ"ל/גרם משקל גוף פעם ביום לאורך תקופת הניסוי. לא הוחל טיפול קדם-אנטיביוטי לפני מתן הפרוביוטיקה. כדי להפחית הטיה ניסויית, בוצעה הערכת תוצאות על ידי חוקרים עיוורים במידת הצורך. MASLD אושר באמצעות ניתוח ביוכימי בדם.
ניתוח דם
בסיום תקופת הניסוי, העכברים היו בצום לילה (12 שעות) עם גישה חופשית למים לפני מדידת גלוקוז ואיסוף דם סופי. העכברים הורדמו עמוקות עם איזופלורן באמצעות מאדה הרדמה עד שאובחן אובדן רפלקס הפדלים. דגימות דם סופיות נאספו באמצעות דימום רטרו-אורביטלי. לאחר איסוף הדם, עכברים הושמתו על ידי פריקה צווארית בזמן הרדמה עמוקה, והמוות אושר על ידי הפסקת הנשימה ופעימות הלב. הדם הורשה להיקרש בטמפרטורת החדר ואז צנטריפוגה ב-3,000 × גרם למשך 10 דקות לקבלת סרום. דגימות הסרום נלקחו ואוחסנו בטמפרטורה של −80 מעלות צלזיוס עד לניתוח. גלוקוז בסרום, כולסטרול כולל (TC), טריגליצרידים (TG), כולסטרול ליפופרוטאין בצפיפות נמוכה (LDL-c), כולסטרול ליפופרוטאין בצפיפות גבוהה (HDL-c), אספרטט אמינוטרנספראז (AST) ואלנין אמינוטרנספראז (ALT) נמדדו באמצעות אנלייזר ביוכימי אוטומטי, בהתאם להוראות היצרן. פרמטרים אלו נמדדו באמצעות בדיקות קולורימטריות אנזימטיות עם ערכות תגובה סטנדרטיות המשולבות במערכת האנלייזרים. כל הפרמטרים הביוכימיים ננורמו לנפח הסרום ובוטאו ביחידות שהוגדרו על ידי היצרן.
בדיקות לחץ חמצוני
סמנים ביולוגיים של לחץ חמצוני נמדדו באמצעות ערכות, כולל אלו של מלונדיאלדהיד, סופראוקסיד דיסמוטאז וקטלז. הבדיקות בוצעו לפי הוראות היצרן. בקצרה, רקמות הכבד עברו הומוגנות בבופר קר וקרח וצנטריפוגה לקבלת סופרנטנטים לניתוח. רמות המלונדיאלדהיד (MDA) נקבעו באמצעות מבחן חומרים תגובתיים בחומצה תיובארבטורית (TBARS); פעילות סופר-חמצן דיסמוטאז (SOD) נמדדה על ידי יכולתו לעכב תגובות מתווכות סופראוקסיד; ופעילות הקטלז (CAT) נמדדה על ידי ניטור קצב הפירוק של מי חמצן. המדידות בוצעו באמצעות קורא מיקרופלט, והערכים חושבו באמצעות עקומות סטנדרטיות או הנוסחאות שסופקו בערכות. סמני לחץ חמצוני ננורמו למשקל רקמת הכבד והובאו לביטוי בהתאם למפרט היצרן.
ריצוף 16S rRNA וניתוח מיקרוביום
דגימות צואה נאספו בסוף תקופת הניסוי ישירות מעכברים בודדים בתנאים סטריליים, הוקפאו מיד בחנקן נוזלי, ואוחסו בטמפרטורה של −80°C עד להוצאת DNA. DNA מיקרוביאלי הופק באמצעות ערכה מסחרית בהתאם לפרוטוקול היצרן. אזור V3–V4 של גן חיידקי 16S rRNA הוגבר באמצעות פריימרים 341F/806R. הגברת PCR בוצעה במערכת תגובה של 25 מיקרוליטר שכללה תבניות DNA, 2x PCR master mix, ו-0.2 מיקרומום מכל פריימר. תוכנית ההגברה כללה דנאטורציה ראשונית ב-95 מעלות צלזיוס למשך 3 דקות, ואחריה 25 מחזורים של 95 מעלות צלזיוס למשך 30 שניות, 55 מעלות צלזיוס ל-30 שניות, ו-72 מעלות צלזיוס למשך 45 שניות, עם הארכה סופית ב-72 מעלות צלזיוס למשך 10 דקות. ספריות Amplicon נבנו ורצפו על פלטפורמת Illumina MiSeq באמצעות אסטרטגיית 2 × 250 bp משולבות על ידי ספק שירותי רצף מסחרי. איכות הספריות הוערכה לפני הרצף, וספריות מאוחדות נטענו בצפיפות האשכול המומלצת על ידי היצרן. קריאות באיכות נמוכה עם בסיסים מעורפלים, אורך לא מספק או ציון איכות ממוצע מתחת ל-Q20 הוסרו לפני ניתוח במעבר.
קריאות רצף גולמיות עובדו באמצעות QIIME 2 גרסה 2019.420. בקצרה, קריאות בקצה זוגי עברו דה-מולטי-פלקס, סינון איכות, דה-נויז, מיזוג וסינון כימרה באמצעות DADA2 עם פרמטרים ברירת מחדל אלא אם צוין אחרת. וריאנטים של רצף אמפליקון לא-סינגלטון (ASVs) נשמרו לניתוח במורד הזרם. הקצאה טקסונומית בוצעה באמצעות גרסת מסד הנתונים Greengenes גרסה 13_8 שמיושמה ב-QIIME221. בהתחשב במגבלות הרזולוציה הטבועות בריצוף גנים של 16S rRNA, במיוחד בשימוש במאגר Greengenes, ההערות הטקסונומיות פורשו אך ורק ברמת הסוג. סיווג ברמת המין לא נחשב אמין ולכן לא שימש בניתוחים במורד הזרם, בהשוואות סטטיסטיות או בפרשנויות ביולוגיות. כל התכונות המיקרוביאליות התכווצו לרמת הסוג לפני ההדמיה והניתוח.
נעשה שימוש במדדי אלפא-גיוון (Chao1, מינים נצפים, שאנון, סימפסון, PD של פיית', ואחידות פיילו). השוואות סטטיסטיות של גיוון אלפא בין קבוצות בוצעו באמצעות מבחנים לא פרמטריים מתאימים (למשל, מבחן קרוסקל–ווליס). מרחק ג'אקארד והשוני בין בריי-קרטיס שימשו כמדדי גיוון בטא. גרפים של ניתוח קואורדינטות עיקרי (PCoA), קנה מידה רב-ממדי לא-מטרי (NMDS) וגודל אפקט ניתוח דיסקרימיננטי ליניארי (LEfSe) שימשו להמחשת שינויים בקהילות מיקרוביום, עם סף ציון LDA של 2.0. ניתוח ורב-משתני פרמוטציוני של שונות (PERMANOVA) שימש להערכת מובהקות סטטיסטית בגיוון בטא. ויזואליזציה של נתונים וניתוחים סטטיסטיים בוצעו באמצעות פלטפורמת GenesCloud (https://www.genescloud.cn).
פרופילים תפקודיים מיקרוביאליים הוסקו מנתוני ריצוף 16S rRNA ולא נמדדו ישירות באמצעות ריצוף מטגנומי או טרנסקריפטומי ברובה ציד. בהתבסס על מסלול מנורמל או טבלת שפע קבוצות פונקציונליות, יחידות פונקציונליות חזויות מופו למסדי נתונים נפוצים, כולל KEGG, MetaCyc ו-COG. שפע מסלולי KEGG סוכם לפי קטגוריות הפונקציונליות ההיררכיות שלו, כאשר סיווגי מסלולים ברמה 2 שימשו להשוואה במורד הזרם. השפע הממוצע של כל קטגוריית מסלולים חושב באמצעות תוכנת R.
ניתוח מטבולומיקה לא ממוקדת של דגימות כבד
רקמות כבד נכרתו במהירות מיד לאחר אישור המוות, הוקפאו מיד בחנקן נוזלי, ונשמרו בטמפרטורה של −80°C לפני הניתוח. כ-50 מ"ג רקמת כבד קפואה עברו הומוגניזציה במתנול:מים קרים (4:1, v/v) בתוספת תקנים פנימיים. יחס הממס בין רקמה להוצאה היה 1:10 (w/v). הדגימות הועברו לוורטקס למשך דקה אחת, סוניקציה על קרח למשך 10 דקות, ואז דגרו בטמפרטורה של −20 מעלות צלזיוס למשך 30 דקות לשיפור משקעי החלבון. ההומוגנטים עורבבו היטב, עברו סוניקציה על קרח, וצנטריפוגה ב-12,000 × גרם למשך 15 דקות בטמפרטורה של 4 מעלות צלזיוס. הסופרנטנטים נאספו בקפידה, התאדו בזרם חנקן קל, והורכבו מחדש ב-50% מתנול לניתוח מאוחר יותר.
דגימות בקרת איכות (QC) נוצרו על ידי איסוף נפחים שווים מכל דגימה ונותחו לסירוגין לאורך כל הריצה להערכת יציבות אנליטית ושחזוריות. ניתוח מטבולומיקה לא ממוקדת בוצע באמצעות מערכת כרומטוגרפיה נוזלית בעלת ביצועים גבוהים במיוחד המחוברת לספקטרומטריית מסה ברזולוציה גבוהה, המופעלת על ידי ספק שירות מסחרי. המטבוליטים הופרדו בעמודת C18 בפאזה הפוכה באמצעות מערכת ממס בינארית המורכבת ממים עם חומצה פורמית של 0.1% ואצטוניטריל עם 0.1% חומצה פורמית. נפח ההזרקה היה 2 מיקרוליטר, טמפרטורת העמודה נשמרה על 40 מעלות צלזיוס, וקצב הזרימה נקבע על 0.30 מ"ל לדקה. תוכנית שחרור גרדיאנט שימשה להשגת הפרדה רחבה בין מטבוליטים. זיהוי ספקטרומטרי מסה בוצע הן במצבי יינון אלקטרוספריי חיובי והן שלילי בטווח m/z של 70–1,000 כדי להשיג כיסוי רחב של מטבוליטים. מתח ריסוס היונים הוגדר ל-3.5 kV במצב חיובי ו-−2.5 kV במצב שלילי. טמפרטורת הנימים נשמרה על 320 מעלות צלזיוס, והנתונים נאספו במצב MS מלא עם רכישת MS/MS תלויה בנתונים לצורך הערות מטבוליטים. הנתונים הגולמיים עובדו באמצעות תוכנת מטבולומיקה ייעודית (MassLynx v4.1) לזיהוי תכונות, יישור ונרמול. זיהוי מטבוליטים בוצע בהתבסס על מסה מדויקת, זמן שימור ודפוסי פיצול בין MS/MS בהשוואה למאגרי מידע ציבוריים, כולל מסד הנתונים של Human Metabolome (HMDB) ואנציקלופדיית קיוטו לגנים וגנומים (KEGG) (v2023), כאשר הדבר רלוונטי.
ניתוח רב-משתני וניתוח מסלולים
נתוני מטבולומיקה מעובדים עברו ניתוחים סטטיסטיים רב-משתניים, כולל ניתוח רכיבים עיקריים (PCA) וניתוח הבחנה חלקית של ריבועים (PLS-DA), כדי להמחיש הבדלים מטבוליים בין קבוצות ניסוי. מטבוליטים שונים זוהו בהתבסס על שילוב של חשיבות משתנה בציוני הקרנה (VIP) (>1.0) ומובהקות סטטיסטית (p < 0.05). מטבוליטים שזוהו מופו עוד למסלולים מטבוליים באמצעות ניתוח העשרת מסלולי KEGG להבהרת תהליכים ביולוגיים הקשורים לתוספי שמן מיקרואצות עשירים ב-DHA.
ניתוח סטטיסטי
הנתונים הוצגו כממוצע ± SEM ונותחו באמצעות ניתוח חד-כיווני של שונות (ANOVA), ואחריו מבחן ההשוואות המרובות של טוקי. המובהקות הסטטיסטית הוגדרה כ-p < 0.05. ניתוח העשרה של מטבוליטים בכבד בוצע באמצעות R (גרסה 4.4.1). אסוציאציות זוגיות בין סוגים מיקרוביאליים למטבוליטים בכבד נותחו באמצעות קורלציות דירוג ספירמן. רק מתאמים משמעותיים (שיעור גילוי שגוי, FDR < 0.05) עם |ρ| ≥ 0.6 נכללו בניתוחי הרשת האינטגרטיבית. עיבוד נתוני מיקרוביום ומטבולומיקה בוצע כפי שתואר לעיל, כולל שימוש בפלטפורמת GenesCloud ונהלי בקרת איכות.