כאן אנו מציגים פרוטוקול המשלב ניתוח מולטי-אומיקס עם אימות in vivo כדי לחקור כיצד Gegen Qinlian Decoction מקל על קוליטיס כיבית על ידי ויסות מטרות הקשורות ל-PANoptosis ומסלולי איתות דלקתיים.
הרשמה ל-JoVE נדרשת לצפייה בתוכן זה. התחברו או התחילו בגרסת ניסיון בחינם.
מאמר מחקר
* These authors contributed equally
כאן אנו מציגים פרוטוקול המשלב ניתוח מולטי-אומיקס עם אימות in vivo כדי לחקור כיצד Gegen Qinlian Decoction מקל על קוליטיס כיבית על ידי ויסות מטרות הקשורות ל-PANoptosis ומסלולי איתות דלקתיים.
אנו מציגים פרוטוקול משולב של ביואינפורמטיקה ופרמקולוגיה רשתית להבהרת המנגנונים הטיפוליים של דקוקציה גגן קינליאן (GQD) בקוליטיס כיבית (UC), עם דגש מיוחד על ויסות הקשור ל-PANoptosis. פרוטוקול זה נועד לזהות באופן שיטתי תרכובות ביואקטיביות פוטנציאליות, מטרות טיפוליות, מסלולי איתות ומנגנונים הקשורים למערכת החיסון באמצעות כרייה רב-מסדי נתונים וניתוח טרנסקריפטומי. פרופילי ביטוי גנים מקבוצות UC, כולל GSE87466 ו-GSE92415 ממאגר GEO, מנותחים באמצעות ניתוח ביטוי שונה וניתוח רשתות ביטוי משולב משוקלל לזיהוי גנים הקשורים למחלות. מטרות מועמדות ל-GQD נאספות ממאגרי מידע ציבוריים של פרמקולוגיה ומצטלבות עם גנים הקשורים ל-UC ו-PANoptosis כדי לקבוע מטרות טיפוליות מרכזיות. לאחר מכן מתבצעים בניית רשת אינטראקציה של חלבון-חלבון, ניתוחי העשרה, ומידול רשת תרכובת-יעד-מסלול כדי לחשוף את הבסיס המולקולרי של פעילות GQD. בנוסף, ניתוח חדירה חיסונית משולב כדי לחקור את הקשר בין גנים מרכזיים למיקרו-סביבה החיסונית במעי. עגינה מולקולרית משמשת גם להערכת פוטנציאל הקישור בין תרכובות פעילות עיקריות לבין מטרות המרכז. פרוטוקול זה מספק מסגרת שיטתית וניתנת לשחזור לחקירת האופן שבו GQD מקל על UC על ידי ויסות PANoptosis ומסלולים דלקתיים נלווים, ויכול לשמש כהפניה מועילה למחקרים מכניים של נוסחאות רפואה סינית מסורתית אחרות במחלות דלקתיות מורכבות.
קוליטיס כיבית (UC), אחת הצורות העיקריות של מחלת מעי דלקתית, היא הפרעה דלקתית כרונית ונסגנת, הכוללת בעיקר את רירית המעי הגס ומתחילה בדרך כלל ברקטום עם הארכה פרוקסימלית רציפה. הפתוגנזה שלו מורכבת וכוללת הפעלה חרירית, תפקוד לקוי של מחסום אפיתלי, הפרעות מיקרוביאליות ואינטראקציות לא מווסתות בין מארח לסביבה, וביחד מטיל נטל בריאותי עולמי משמעותי וגדל 1,2. דלקת רירית מתמשכת לא רק גורמת לתסמינים חוזרים כמו כאבי בטן, שלשולים והמטוצ'זיה, אלא גם מגבירה את הסיכון ארוך הטווח לסרטן המעי הגס הקשור לקוליטיס, במיוחד אצל מטופלים עם מחלה נרחבת וארוכת טווח 3,4.
הניהול הפרמקולוגי הנוכחי של UC מתבסס בעיקר על 5-אמינוסליצילטים, קורטיקוסטרואידים, אימונומודולטורים, ביולוגיים וסוכנים למולקולות קטנות, בהתאם לחומרת המחלה ולתגובה הטיפולית5. למרות שאסטרטגיות אלו שיפרו את התוצאות הקליניות, חלק ניכר מהמטופלים עדיין חווים אי-תגובה ראשונית, אובדן תגובה משני, אי סבילות, תלות בסטרואידים או תופעות לוואי הקשורות לטיפול, מה שמדגיש את הצורך באפשרויות טיפוליות בטוחות וממסוגלות מנגנון6. לכן, זיהוי תוכניות מולקולריות במעלה הדרך המווסתות בו-זמנית את הפגיעה האפיתלית ואת חוסר הוויסות החיסוני של המעי הוא בעל חשיבות תרגומית משמעותית.
מחקרים עדכניים הראו שמוות תאי מתוכנת לא מוסדר מעורב עמוקות בפתוגנזה של UC. מעבר לאפופטוזיס הקלאסי, פירופטוזיס ונקרופטוזיס כולם היו מעורבים בהרס מחסום אפיתל ובהגברה של דלקת רירית7. PANoptosis, מסלול מוות תאי מתוכנת דלקתי שהוגדר לאחרונה ומשלב תכונות מרכזיות של פירופטוזיס, אפופטוזיס ונקרופטוזיס, הפך למסגרת מכנית רלוונטית במיוחד למחלות דלקתיות כרוניות8. תהליך זה מתואם על ידי קומפלקסים מרובי חלבונים של PANoptoזום וכולל מתווכים מרכזיים כמו ZBP1, RIPK1/RIPK3, CASP8, ומולקולות הקשורות לאינפלמזום כולל NLRP3 באופן תלוי הקשר9. חשוב לציין שמצטברות ראיות מצביעות על כך ש-PANoptosis מופעלת בקיבוריס (UC) ועלולה לקדם נזק לאפיתל במעי, קריסת מחסום וחוסר איזון חיסוני, ובכך לתרום להתקדמות המחלה 10,11,12. ממצאים אלו מצביעים על כך ש-PANoptosis אינה רק תוצאה המשך של דלקת, אלא עשויה לייצג צומת פעולה מכנית המקשרת בין פגיעה אפיתלית להגברה חיסונית רירית.
גגן צ'ינליאן דקוקציה (GQD), מרשם קלאסי שנרשם לראשונה בשאנג האן לון, מורכב מפוארריה לובאטה, סקוטלריה בייקלנסיס, קופטיס כיננסיס וגליסיריזה אורלנסיס. הוא משמש זמן רב ברפואה הסינית המסורתית לשלשול והפרעות דלקתיות במעי. מחקרים פרמקולוגיים מודרניים הראו כי GQD יכול להקל על קוליטיס ניסיוני על ידי דיכוי תגובות דלקתיות, שיקום שלמות מחסום המעיים, עיצוב מחדש של מיקרוביוטה במעיים ווויסות הומאוסטזיס חיסוני, כולל איתות ופעילות אינפלמזום הקשורה ל-Th17/Treg13,14. פעולות מרובות מטרות אלו הופכות את GQD למועמד סביר להתערבות בתהליכים פתולוגיים הקשורים ל-PANoptosis ב-UC15. עם זאת, למרות הראיות הגוברות ליעילותו האנטי-קוליטית, עדיין לא ברור אם GQD מפעיל השפעות טיפוליות באמצעות ויסות מתואם של רשתות מולקולריות הקשורות ל-PANoptosis.
נכון לעכשיו, נותרו שני פערים עיקריים. ראשית, מחקרים קיימים על מנגנוני מוות של תאי UC התמקדו בעיקר במסלולים מבודדים כגון אפופטוזיס, פירופטוזיס או פרופטוזיס, בעוד שהתפקיד המשולב של PANoptosis בתיווך ההשפעות הטיפוליות של נוסחאות צמחיות לא הובהר מספיק 7,10,11. שנית, אף שפרמקולוגיה ברשת מתאימה לאפיון האופי של "רב-רכיבים-רב-מטרות-ריבוי-מסלולים" של נוסחאות מסורתיות, תחזיות המבוססות על מסדי נתונים מורכבים-יעדים בלבד עשויות להיות חסרות ספציפיות להקשר מחלה ולדרוש עדיפות טרנסקריפטומית ואימות ניסויי לשיפור הסבירות הביולוגית ואמינות התרגום 16,17,18,19 . לכן, אסטרטגיה משולבת המשלבת טרנסקריפטומיקה של מחלות, פרמקולוגיה רשתית ואימות in vivo עשויה לספק מסגרת קפדנית יותר להבהרת מנגנון GQD ב-UC.
בהשוואה לגישות המבוססות אך ורק על תצפית היסטופתולוגית, בדיקות מולקולריות חד-נתיביות או פרמקולוגיה רשתית מונחית מסדי נתונים בלבד, זרימת העבודה הנוכחית מתאימה יותר לשאלה זו משום שהיא מאפשרת אימות צולב בין אותות טרנסקריפטומיים בהקשר מחלה, חיזוי של נוסחה-יעד ואימות ביולוגי במודל קוליטיס ניסיוני. קריאות היסטולוגיות או מולקולריות לבדן יכולות ללכוד רק ממדים מוגבלים של פתולוגיית UC, בעוד שפרמקולוגיה רשתית לבדה עשויה להיעדר ספציפיות מחלה מספקת ותמיכה ניסויית. על ידי שילוב תעדוף טרנסקריפטומי עם אימות in vivo, תהליך עבודה זה משפר את הסבירות הביולוגית של זיהוי מטרה ומספק מסגרת קוהרנטית יותר לקישור התערבות GQD למנגנונים הקשורים ל-PANoptosis ב-UC.
פרוטוקול זה מתאים במיוחד למחקרים מכניים שמטרתם להעריך כיצד התערבויות רב-רכיביות מווסתות פגיעה אפיתלית, תגובות דלקתיות ונתיבים הקשורים ל-PANoptosis ב-UC ניסיוני. הוא שימושי במיוחד כאשר המטרה היא לחבר בין חיזוי חישובי לאימות ביולוגי ברמת בעלי חיים בזרימת עבודה סטנדרטית. עם זאת, הוא אינו נועד להחליף גישות אבחון קליניות, ואינו מטפל באופן מקיף בכל היבטי הפתוגנזה של ה-UC, כגון הטרוגניות ארוכת טווח בין מטופלים, אינטראקציות סיבתיות כלל-מיקרוביום, או מיפוי חיסון חד-תאי ומרחבי. לכן, היקפו מוגדר בצורה הטובה ביותר כמסגרת פרה-קלינית ניתנת לשחזור להערכה מונחית מנגנונים וטיפולית.
לכן השענו ש-GQD משפר את ה-UC על ידי ויסות רשתות איתות הקשורות ל-PANoptosis, ובכך מפחית את הפגיעה באפיתל המעי ומתכנת מחדש את המיקרו-סביבה הדלקתית. כדי לבדוק השערה זו, שילבנו ניתוח טרנסקריפטומי של UC עם פרמקולוגיה רשתית כדי לזהות מטרות טיפוליות משותפות ונתיבים מרכזיים של GQD, עם דגש מיוחד על רגולטורים הקשורים ל-PANoptosis. אימתנו עוד את המנגנונים החזויים במודל דלקת קוליטיס עכברית מושרה על ידי דסקסטרן סולפט נתרן (DSS). על ידי שילוב עדיפות מטרה חישובית עם אימות ביולוגי, מחקר זה נועד לספק הסבר רלוונטי יותר למחלות ומבוסס מכניסטית יותר להשפעות האנטי-UC של GQD, ולהציע ראיות חדשות התומכות ב-PANoptosis כציר טיפולי פוטנציאלי ב-UC.
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
כל ההליכים בבעלי חיים בוצעו בהתאם להנחיות מוסדיות ואושרו לפני תחילת המחקר על ידי ועדת האתיקה של אוניברסיטת צ'נגדו לרפואה סינית מסורתית (אישור מס' 2022-126). רק מערכי נתונים זמינים לציבור של GEO ו-GWAS שימשו לניתוחי ביואינפורמטיקה; לכן, לא נדרשה הסכמה מדעת נוספת או אישור אתי נוסף עבור הרכיב החישובי של מחקר זה. כל ההליכים הכוללים כימיקלים מסוכנים, דגימות ביולוגיות ופסולת בעלי חיים בוצעו בהתאם לתקנות בטיחות מעבדות מוסדיות ולנהלים מוסדיים להשלכת פסולת מסוכנת.
רכישה ועיבוד מוקדם של מערכי נתונים של GEO
מערכי נתונים טרנסקריפטומיים אנושיים הקשורים לקוליטיס כיבית (UC) נאספו ממאגר הנתונים GEO באמצעות מילת המפתח "קוליטיס כיבית". לאחר סינון למאגרי פרופיל ביטוי עם קבוצות UC וביקורת בריאה מוגדרות היטב וגודל מדגם מספק, נבחרו GSE8746620 ו-GSE92415 21, שניהם שנוצרו בפלטפורמת GPL13158, כקבוצת ההדרכה. מערכי הנתונים הללו כללו 87 דגימות UC ו-21 ביקורות בריאות, ו-53 דגימות UC ו-21 ביקורות בריאות, בהתאמה. לצורך אימות חיצוני, GSE8747322 (106 דגימות UC ו-21 ביקורת בריאה)ו-23 GSE16879 (48 דגימות UC ו-12 ביקורת בריאה; פלטפורמת GPL570) נכללה. קבצים גולמיים או במטריצה סדרתית הורדו ויובאו ל-R לצורך עיבוד מוקדם. מזהי הגשושים הוסבו לסמלי גנים בהתאם לקבצי ההערה של הפלטפורמה המתאימים. כאשר מספר גלאים מופו לאותו סמל גן, ערך הביטוי הממוצע נשמר. דגימות ללא מידע ברור על קבוצות הוצאו לפני הניתוח במורד הזרם. שני מערכי הנתונים של האימון אוחדו על ידי סמלי גנים משותפים. ערכי הביטויים עברו טרנספורמציה ל-log2 בעת הצורך, נורמליזציה באמצעות הליך סטנדרטי בין מערך, ותוקנו לאפקטים אצווה באמצעות שיטת התאמת אצווה לאחר הגדרת מקור מערך הנתונים כמשתנה באץ. ניתוח רכיבים עיקריים ובדיקת קופסאות שימשו לפני ואחרי תיקון הנורמליזציה/אצווה כנקודות ביקורת לבקרת איכות. עיבוד מוקדם מוצלח הוגדר על ידי חפיפה משופרת של התפלגויות דגימות בין מערכי נתונים והפחתת אשכולות מונחי מערכי נתונים. בסך הכול נאספו 902 גנים הקשורים ל-PANoptosis מתוך ספרות שפורסמה ושימשו כמערך גנים ייחוס לניתוח אינטגרטיבי.
זיהוי גנים מבוטאים באופן שונה ובניית רשת WGCNA
מטריצת האימון המאוחדת נותחה ב-R באמצעות זרימת עבודה של ביטויים דיפרנציאליים. גנים עם |log2 fold משתנים| ≥ 0.585 ו-P מותאם < 0.05 הוגדרו כגנים מבוטאים באופן שונה (DEGs). גרפים של הרי געש ומפות חום נוצרו כתוצרים ביניים כדי לאשר כי קריטריוני הסינון הניבו הבדלים ביולוגיים בביטוי בין קבוצות UC לקבוצת הביקורת. בניתוח רשת ביטוי משולבת של גנים משוקללים (WGCNA), הגנים דורגו לפי שונות בין מדגמים, ו-25% הגנים המשתנים ביותר נשמרו כקלט. עץ אשכולות דגימות נבדק תחילה כדי לזהות דגימות חריגות פוטנציאליות; לא נשמרו חריגים ברורים לבניית רשת אלא אם ההסרה מוצדקת לפי קריטריוני בקרת איכות. הספק רך נבחר בהתבסס על מדד התאמה טופולוגי חופשי בקנה מידה, והחזקה הקטנה ביותר שהגיעה לרשת נטולת קנה מידה בערך שימשה לבניית מטריצת השכונות. מטריצת השכנות הוסבה לאחר מכן למטריצת חפיפה טופולוגית, ובוצעה אשכולות היררכיים באמצעות אי-דמיון מבוסס TOM. מודולים זוהו באמצעות כריתת עצים דינמית עם גודל מודול מינימלי של 100 גנים. מודולים קרובים אוחדו כאשר הקורלציה העצמית שלהם עלתה על קריטריון המיזוג המוקדם. חושבו ערכי מורכבות עצמיים של מודול, משמעות מודול וערכי מובהקות גנים כדי לזהות את המודולים הקשורים ביותר לפנוטיפ של UC. שלב מוצלח ב-WGCNA הוגדר על ידי הפרדה יציבה של מודולים, גודל מודולים ביולוגיים סבירים, ומודול אחד או יותר שהראו קורלציה ברורה למצב המחלה.
זיהוי מטרות מועמדות ב-GQD
תרכובות מועמדות ויעדים התואמים לארבעת הרכיבים הצמחיים של Gegen Qinlian Decoction (GQD)—Pueraria lobata, Scutellaria baicalensis, Coptis chinensis ו-Glycyrrhiza uralensis—נאספו מ-TCMSP24 ו-BATMAN-TCM25. ספי הסקר נקבעו כך: זמינות ביולוגית דרך הפה ≥ 30% ודמיון תרופה ≥ 0.18 עבור רשומות שמקורן במאגר נתונים אחד; ציון הביטחון ≥ 0.84 ו-P מותאם <-0.05 עבור רשומות שמקורן במאגר אחר. תרכובות נוספות הושלמו ממאגר תרופות צמחיות כאשר לא נאספו רשומות רלוונטיות בחיפוש הראשוני. שרשראות קנוניות של תרכובות שמורות ב-SMILES הוגשו לפלטפורמת הערכה מקוונת של ADME. התרכובות נשמרו כאשר הראו ספיגה גבוהה במערכת העיכול ועמדו לפחות בשני חוקי דמיון תרופתי בין ליפינסקי, גוזה, ובר, איגן ומואגה. היעדים המשוערים נחזו לאחר מכן באמצעות פלטפורמת חיזוי מטרות עם הסתברות > 0.1 וסטנדרטים לסמלים גנים רשמיים באמצעות מאגר נתונים של אנוטציות חלבון. תרכובות כפולות ומטרות כפולות הוסרו לאחר מיזוג מסדי נתונים. התוצר הצפוי של שלב זה היה מערך נתונים GQD לא מיותר עם יעד מורכב, המתאים לניתוח חפיפה.
זיהוי גנים חופפים וביצוע ניתוח פונקציונלי
מטרות טיפוליות פוטנציאליות הוגדרו כהצטלבות בין מטרות הקשורות ל-GQD, גנים הקשורים ל-UC שזוהו מביטוי שונה ו-WGCNA, וקבוצת הגנים הקשורה ל-PANoptosis. יחסי חפיפה הוצגו באמצעות דיאגרמת ון. הגנים המצטלבים נשלחו למאגר אינטראקציה בין חלבון לחלבון (PPI), וזוגות אינטראקציה עם ציון ביטחון ≥ 0.4 נשמרו. רשת ה-PPI שנוצרה יוצאה והוצגה בפלטפורמת ניתוח רשת. שלב PPI מוצלח הוגדר על ידי נוכחות רשת מחוברת הכוללת את רוב הגנים המצטלבים, ולא רק צמתים מבודדים. ניתוח העשרה פונקציונלית, כולל אנליזה של אונטולוגיית גנים (GO) ואנציקלופדיית קיוטו של גנים וגנומים (KEGG), בוצע ב-R כאשר q < 0.05 כסף המובהקות. תנאים או מסלולים שלא עמדו בסף זה הוצאו. התהליכים הביולוגיים והנתיבי האיתות העשירים ביותר סוכמו כדי להסיק את המנגנון הפוטנציאלי שבאמצעותו GQD עשויה לווסת את UC. גרפים ברים וגרפים בועות נוצרו כתוצרים ביניים כדי לוודא האם תוצאות ההעשרה נשלטות על ידי מסלולים הקשורים לדלקת, חיסון או פציעות אפיתליות, כפי שמצפים מההקשר הביולוגי.
ניתוח דפוסי ביטוי, לוקליזציה כרומוזומלית ומתאמים
רמות הביטוי של הגנים המשותפים הושוו בין קבוצות ה-UC לקבוצת הביקורת באמצעות תהליך הדמיה סטנדרטי של הבעה שונה ב-R. נוצרו תרשימי קופסה ומפות חום, ו-P < 0.05 נחשב מובהק סטטיסטית. רק גנים עם כיוון ביטוי ניתן לפרשנות והתנהגות עקבית לאורך קבוצת ההכשרה נשמרו להערכה מכנית מאוחרת יותר. מיקומי כרומוזומים של גנים ליבה נבחרים מופו לאחר עיצוב קבצי האנווטציה של הגנים למבנה הקלט הנדרש, והמיקומים הגנומיים הוצגו באמצעות תהליך גרוטינג כרומוזומים מעגלי. מטריצות קורלציה בין גנים ליבה חושבו לאחר מכן באמצעות ניתוח מתאם זוגי. התוצאה הצפויה של שלב זה הייתה זיהוי זוגות גנים המציגים דפוסי ביטוי מתואמים שיכולים לתמוך בוויסות משותף או באינטראקציה פונקציונלית.
בנייה והערכה של מודלים של למידת מכונה
פיתוח מודלים אבחנתיים בוצע באמצעות קבוצת ההדרכה המאוחדת (GSE87466 + GSE92415), בעוד GSE16879 ו-GSE87473 שימשו כמאגרי נתונים חיצוניים עצמאיים לאימות חיצוני. שנים-עשר אלגוריתמים ללמידת מכונה - Lasso, Ridge, Stepglm, XGBoost, Random Forest, Enet, plsRglm, GBM, NaiveBayes, LDA, glmBoost ו-SVM - שימשו ליצירת 113 שילובי מודלים, 26. בוצעה אימות צולב פי עשרה בתוך קבוצת ההכשרה. בכל קיפול בוצעו הדרכת מודלים, בחירת תכונות והערכת ביצועים באמצעות אותה לוגיקת חלוקה של אימון/בדיקות כדי למנוע דליפת מידע. ביצועי המודל הוערכו בעיקר לפי השטח שמתחת לעקומת מאפיין ההפעלה של המקלט (AUC). המודל האופטימלי הוגדר כמודל עם ממוצע ה-AUC הגבוה ביותר בין מאגרי הנתונים של האימון והוולידציה, ולא כ-AUC האימון הגבוה ביותר בלבד. מודלים שהראו AUC אימון גבוה מאוד אך ביצועי אימות חיצוניים ירודים לא נשמרו. נוצרו עקומות ROC לכל מודל שנשמר, ונבנתה נומוגרמה מהסמנים הביולוגיים הסופיים שנבחרו. שלב מודל מוצלח הוגדר על ידי ביצועי אימות צולב יציבים ושימור יכולת הבחנה במאגרי הנתונים החיצוניים.
ביצועי GSVA וביטוי גנים מרכזיים בפרופיל
ניתוח שונות של מערך גנים (GSVA) בוצע כדי להשוות העשרת מסלולים בין קבוצות ביטוי גבוה לקבוצות ביטוי נמוך, המוגדרות לפי ערך הביטוי החציוני של כל גן ליבה. P < 0.05 מותאם שימש כחיתוך להבדלים משמעותיים במסלולים. מסלולים שלא עמדו בסף זה הוצאו מפרשנות. נוצרו מפות חום וגרפים של קורלציה כדי להמחיש פעילות מסלולים ודפוסי ביטוי משותף. הצלחה בינונית של שלב זה הועלתה על ידי פרופילי העשרה התואמים לדיסרולציה דלקתית, חיסונית או מטבולית הרלוונטית ל-UC.
ניתוח הנוף החיסוני ובניית רשת ceRNA
הרכב תאי המערכת החיסונית הוערך באמצעות אלגוריתם דה-קונבולוציה עם 1,000 פרמוטציות. דגימות שלא עמדו בקריטריון המשמעות של דה-קונבולוציה פנימית הוצאו מניתוח מתאם חיסוני מאוחר יותר. להשלמת ניתוח זה, נעשה שימוש גם בניתוח העשרת קבוצות גנים חד-דגימית (ssGSEA) להערכת דפוסי חדירה חיסונית. לאחר מכן חושבו והוצגו מתאמים בין חלקי תאי החיסון לביטוי הגנים המרכזי. לצורך בניית רשת ceRNA, זוהו miRNA המיקוד ב-mRNA הליבה ממאגרי אינטראקציה זמינים לציבור, ו-lncRNA מועמדים האינטראקטיביים עם miRNA אלו נבדקו לאחר מכן. נשמרו רק אינטראקציות הנתמכות על ידי רשומות חיזוי או הערות במסד נתונים. רשת הרגולציה mRNA-miRNA-lncRNA הוצגה בפלטפורמת ניתוח רשת. התפוקה הצפויה של שלב זה הייתה רשת ceRNA מובנית עם רגולטורים ביולוגיים סבירים במעלה הזרם המקושרים לגנים המרכזיים.
ביצועי ניתוח תא יחיד
GSE214695 מערך הנתונים של תא יחיד יובא ל-R והומור לאובייקטים חד-תאיים באמצעות תהליך ניתוח תא יחיד סטנדרטי. בקרת איכות בוצעה על ידי סינון תאים לפי מדדים סטנדרטיים, כולל תאים בעלי תכונות נמוכות, טיפות פוטנציאליות ריקות, ותאים עם איכות תמלול ירודה. תאים שעברו בקרת איכות ננורמו, ו-2,000 הגנים המשתנים ביותר נבחרו לניתוח רכיבים עיקריים. המרכיבים העיקריים שתרמו שונות ביולוגית משמעותית נשמרו לאשכולות לא מפוקחים, ו-UMAP שימש להמחשה דו-ממדית. גני סמן אשכולות זוהו באמצעות ניתוח ביטוי שונה בין אשכולות. לאחר מכן בוצעה סיווג מסוג תא באמצעות סיווג מבוסס ייחוס, והאנוטציות הושוו מול ביטוי סמני קנוני כאשר הדבר האפשרי. שלב תא יחיד מוצלח הוגדר על ידי הפרדה ברורה של אוכלוסיות תאים עיקריות בגרף UMAP, גנים של סמנים ניתנים לפרשנות לכל אשכול27, והערות עקביות של סוגי התאים החיסוניים והסטרומליים הרלוונטיים ל-UC.
ביצוע ניתוח אקראיות מנדלי
בוצעה רנדומיזציה מנדלית דו-מדגמית (MR) להערכת קשרים סיבתיים אפשריים בין ביטוי הגנים המרכזיים לרגישות ל-UC. משתנים אינסטרומנטליים הופקו ממאגרי נתוני מוקדי תכונות כמותיות של ביטוי, וסטטיסטיקות סיכום מחקר אסוציאציה כלל-גנומי של UC שימשו כמערך הנתונים לתוצאות. SNPs נבדקו לפי קריטריוני בקרת איכות סטנדרטיים של MR, כולל רלוונטיות לתכונת החשיפה והסרת וריאנטים עמומים או כפולים כאשר נדרש על ידי הרמוניזציה. שיטת השונות ההפוכה שימשה כאומדת סיבתית ראשית. מבחן Q של קוקרן יושם להערכת הטרוגניות. ניתוחי MR-Egger ו-MR-PRESSO שימשו כניתוחי רגישות להערכת וכאשר נדרש, תיקון פליאוטרופיה אופקית. שלב MR מוצלח הוגדר על ידי הרמוניזציה תקפה של מכשירים, היעדר הטרוגניות משמעותית או פליאוטרופיה28, והערכות עקביות כיוונית בין שיטות MR משלימות.
ביצוע עגינה מולקולרית של תרכובות מועמדות
מבנים תלת-ממדיים של חלבוני מטרה התקבלו מבנק נתוני החלבונים, ומבני ליגנדים התקבלו ממאגר מולקולות קטנות ציבורי. לפני העגינה, חלבונים וליגנדים עובדו מראש על ידי הסרת מולקולות מים בעת הצורך, הוספת מימן, הגדרת סוגי אטומים והמרת פורמטים של קבצים הנדרשים לעגינה. העגינה בוצעה לאחר מכן באמצעות תהליך עגינה מולקולרית, ואנרגיות קישור חושבו עבור כל זוג ליגנדה-יעד. זיקות קשירה נמוכות מ-5 קק"ל למול פורשו כקשירה מועדפת, בעוד שערכים נמוכים מ-7 קקלורי/מול פורשו כקשירה יחסית חזקה. מפות חום של עגינה צוירו כדי להשוות דפוסי קישור גלובליים בין מטרות. הקונפורמציות המדורגות ביותר הוצגו בשלושה ממדים כדי לבחון קשרי מימן, אינטראקציות הידרופוביות וקרבה לשאריות. השלמה מוצלחת של שלב זה הוגדרה על ידי קבצי פלט עגינה יציבים, תנוחות ליגנד סבירות באזור הקישור המטרה, ואנרגיות קישור התואמות את תעדף המועמדים.
ביצוע אימות ניסויי in vivo
הכנת GQD
GQD הוכן מארבעה עשבי מרפא סיניים מסורתיים בכמויות גסות של תרופות: 24 גרם Puerariae Lobatae Radix (Gegen), 9 גרם Scutellariae Radix (Huangqin), 9 גרם Coptidis Rhizoma (הואנגליאנית), ו-6 גרם Glycyrrhizae Radix et Rhizoma (גנקאו). פרטים על החומרים הבוטניים ומקורותיהם מופיעים בטבלת החומרים. העשבים נטחנו לאבקה דקה, התערבבו היטב והושרו ב-390 מ"ל מים מזוקקים למשך שעה. התערובת הוטחה פעמיים. לתייבש הראשון, העשבים הושגו על חום גבוה ואז בושלו על אש נמוכה במשך 30 דקות. לאחר סינון דרך גזה, השאריות הורתחו מחדש עם נפח מים פי עשרה. שני הסינון שולבו ורוכזו באמבט מים עד לריכוז תרופה גולמי סופי של 5 גרם למ"ל. התמצית סיננתה, התקררה ונשמרה בטמפרטורה של 4 מעלות צלזיוס עד לשימוש.
הקמת מודל העכבר UC המושר על ידי DSS
במחקר זה שימשו שמונה-עשר עכברים זכרים מסוג BALB/c (בני 8 שבועות), וכל ההליכים אושרו על ידי ועדת האתיקה של אוניברסיטת צ'נגדו לרפואה סינית מסורתית (אישור מס' 2022-126). פרטים על מקור בעלי החיים מופיעים בטבלת החומרים. העכברים אוחסנו בתנאים מסוימים ללא פתוגנים בטמפרטורה של 20–22 מעלות צלזיוס ולחות של 55% במחזור אור/חושך של 12 שעות. לאחר תקופת הסתגלות של שבוע, החיות חולקו באקראי לשלוש קבוצות (n = 6 לכל קבוצה): NC, UC ו-GQD. קבוצת NC קיבלה מים סטריליים לשתייה לאורך כל הניסוי, בעוד שקבוצות UC ו-GQD קיבלו 1.5% (עם משקל דקסטרן סולפט) נתרן (MW 36-50 kDa) במי שתייה במשך 5 ימים רצופים, ואחריהם 2 ימים של מים רגילים. מחזור זה חזר על עצמו שלוש פעמים כדי לבסס קוליטיס כרוני. במהלך תקופת ההתערבות, קבוצות NC ו-UC קיבלו 200 מיקרוליטר של 0.5% נתרן קרבוקסימתיל צלולוז, בעוד שקבוצת GQD קיבלה GQD פעמיים ביום.
המרשם הקליני ל-GQD כלל 24 גרם גגן, 9 גרם חואנגצ'ין, 9 גרם חואנגליאן ו-6 גרם גנקאו, המקביל לכמות תרופה גולמית כוללת של 48 גרם ליום למבוגרים. בהנחה שמשקל גוף מבוגר סטנדרטי של 60 ק"ג, המינון המקביל למבוגרים היה 0.8 גרם/ק"ג/יום. בהתבסס על מקדם ההמרה המקובל בין אדם לעכבר של 1:12, המינון המקביל לעכבר היה 9.6 גרם לק"ג/יום. לפי יחס מינונים נמוכים, בינוניים וגבוהים שנקבע מראש של 1:2:4, המינונים המחושבים היו 9.6, 19.2 ו-38.4 גרם לק"ג/יום, בהתאמה. בהתחשב בחומרת ובהתמדה של דלקת המעי הגס במודל הקוליטיס המושרה על ידי DSS, נבחר משטר המינון הגבוה כדי להבטיח התערבות פרמקולוגית מספקת במהלך אימות in vivo. כדי לתמוך בחשיפה תרופתית מספקת במהלך שלב הדלקת הפעיל, המנה היומית הכוללת של 38.4 גרם/ק"ג/יום חולקה לשתי זריקות שוות. במהלך הניסוי, משקל הגוף, עקביות הצואה ודימום רקטלי נבדקו מדי יום. בסיום הפרוטוקול, עכברים הושמתו ונאספו רקמות מעי גס לניתוחים נוספים.
הערכת מדד משקל גוף ופעילות מחלות
משקל הגוף נרשם מדי יום והוא מוצג כאחוז ממשקל הגוף ההתחלתי ביום 0. חומרת המחלה הוערכה באמצעות מדד פעילות המחלה (DAI), שכלל שלושה פרמטרים: ירידה במשקל, עקביות בצואה ודימום רקטלי. כל פרמטר דורג באופן עצמאי בסולם של 0–4 בהתאם לחומרת התסמינים הקליניים, כאשר ציונים גבוהים יותר הצביעו על פעילות מחלה חמורה יותר. ה-DAI הסופי חושב כממוצע של שלושת הציונים האישיים, והניב ציון כולל שנע בין 0 ל-4. בקצרה, ציון של 0 הצביע על חוסר חריגה, בעוד שציונים של 1–4 שיקפו ירידה הולכת וגדלה במשקל, עקביות צואה רפויה יותר ודימום רקטלי חמור יותר. קריטריוני הניקוד המפורטים מוצגים בטבלה 1.
ביצוע צביעת H&E וחישוב מדד הפעילות ההיסטולוגית
להערכה היסטולוגית, רקמות המעי הגס קיבעו ב-4% פאראפורמלדהיד, יובשו באמצעות אתנול מדורג, הוטמעו בפרפין, וחתכו בעובי של 4 מיקרון. החלקים עברו דה-פרפין, התייבשו מחדש, צובקו בהמטוקסילין, צובעו נגד באאוזין, יובשו, נוקו והותקנו. תמונות היסטולוגיות נצפו באמצעות מיקרוסקופ שדה בהיר ונסרקו באמצעות מערכת סריקת שקופיות. מידע ספציפי על ציוד מופיע בטבלת החומרים. הפגיעה ההיסטופתולוגית הוערכה על פי שלושה פרמטרים, כולל התרוקנות תאי הגביע, נזק מבני הקריפטה וחדירת תאים דלקתיים. כל פרמטר דורג לפי דרגת הפגיעה ברקמה, וסכום שלושת ציוני הפרמטרים הוגדר כמדד פעילות היסטולוגי (HAI), הנע בין 0 ל-10, כפי שמסוכם בטבלה 2. לכימות תאי הגביע, נותחו חמישה שדות לא חופפים שנבחרו באקראי בכל מקטע בהגדלה של 400× מספר תאי הגביע בכל שדה נספר ידנית בצורה עיוורת, והספירה הממוצעת לשדה חושבה עבור כל חלק מעי גס מייצג.
ביצוע צביעה אימונוהיסטוכימית
לאחר דה-פרפיניזציה וההידרציה מחודשת, חתכי הרקמה עברו שאיבת אנטיגן עם 0.1% טריפסין ב-37 מעלות צלזיוס למשך 30 דקות, ולאחר מכן דגירה עם 3% מי חמצן לחסימת פעילות פרוקסידאז אנדוגנית. המקטעים נדגרו במהלך הלילה בטמפרטורה של 4 מעלות צלזיוס עם נוגדנים ראשוניים נגד STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP), ו-HSPA5 (1:1000, טכנולוגיית איתות תאי, #3177), מדוללים ב-1% BSA. פרטים על הנוגדנים ומקורותיהם מופיעים בטבלת החומרים. לאחר שטיפה עם PBS, החתכים הודגרו עם הנוגדן המשני המתאים (1:2000) בטמפרטורת החדר למשך 1.5 שעה. אימונואקטיביות הוצגה באמצעות תגובת זיהוי כרוגני, והגרעינים צוברו נגד בהמטוקסילין. ניתוח תמונה כמותי בוצע באמצעות תוכנת ניתוח תמונה המפורטת בטבלת החומרים.
ביצוע ניתוח qRT-PCR
סך ה-RNA הופק באמצעות תגובת חילוץ RNA מבוסס פנול והועתק הפוך באמצעות ערכת סינתזת cDNA בגדיל ראשון. תמלול הפוך בוצע בטמפרטורה של 25 מעלות צלזיוס למשך 10 דקות, 55 מעלות צלזיוס ל-15 דקות, ו-85 מעלות צלזיוס למשך 5 דקות. לאחר מכן בוצע qRT-PCR באמצעות מערכת PCR בזמן אמת ותערובת SYBR Green qPCR. תנאי ההגברה היו כדלקמן: דנאטורציה ראשונית ב-95 מעלות צלזיוס למשך 60 שניות, ואחריה 40 מחזורים ב-95 מעלות צלזיוס למשך 15 שניות, 60 מעלות צלזיוס ל-15 שניות, ו-72 מעלות צלזיוס למשך 45 שניות. רצפי פריימר עוצבו באמצעות תוכנת עיצוב פריימר ומופיעים בטבלה 3. GAPDH שימש כהתייחסות פנימית, ורמות הביטוי היחסיות של mRNA חושבו באמצעות שיטת 2^-ΔΔCt . כל דגימה נותחה בשלושה עותקים, וכל ניסוי חזר על עצמו לפחות שלוש פעמים. פרטים על המגבים, הכלים והתוכנה מסופקים בטבלת החומרים.
ביצוע ניתוח הכתמים המערביים
רקמות המעי הגס הומוגניזציה בבופר ליזיס המכיל מעכב פרוטאז על קרח ועברו ליזה יסודית במשך שעתיים. ריכוזי החלבון נקבעו באמצעות ערכת בדיקת חלבון. כמויות שוות של חלבון עורבבו עם בופר עומס, עברו דנטור לפני האלקטרופורזה, הופרדו על ידי 10% SDS-PAGE, והועברו לממברנות PVDF ב-100 וולט למשך שעה. לאחר חסימה עם חלב דל שומן של 5% למשך שעה בטמפרטורת החדר, הממברנות נדגרו במהלך הלילה בטמפרטורה של 4°C עם נוגדנים ראשוניים של ארנבים נגד STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP), HSPA5 (1:1000, טכנולוגיית איתות תאים, #3177), ו-GAPDH בדילול העבודה המתאימה. לאחר שטיפה עם TBST (20 mM Tris, 150 mM NaCl, 0.1% Tween-20; pH 7.4), הממברנות נדגרו עם הנוגדן המשני המתאים (1:5000) בטמפרטורת החדר למשך שעה. רצועות החלבון הוצגו באמצעות תגובת גילוי כימילומינסצנטי. עוצמות הרצועות נותחו באמצעות תוכנת ניתוח תמונה, ורמות הביטוי של חלבוני המטרה נורמלו ל-GAPDH. המשקל המולקולרי הצפוי של החלבונים שזוהו היה כ-88 kDa עבור STAT3, 21 kDa עבור TIMP1, 28 kDa עבור SPHK2, 70 kDa עבור HSPA5, ו-36 kDa עבור GAPDH. פרטים על הנוגדנים, הממברנות, המגיבים והתוכנה מסופקים בטבלת החומרים. שמות מוצרים ספציפיים, יצרנים, מספרי קטלוג, דגמי מכשירים ופרטי תוכנה מסופקים בטבלת החומרים, בעוד שתיאורים כלליים משמשים בטקסט הראשי לשיפור הנייטרליות והיישום הכללי.
ביצועי ניתוח סטטיסטי
הנתונים מוצגים כסטיית תקן ממוצעת ±. הבדיקות הסטטיסטיות בוצעו באמצעות תוכנה סטטיסטית. ANOVA חד-כיווני ואחריו מבחן פוסט-הוק של טוקי שימש להשוואות רב-קבוצתיות. ערך P ≤ 0.05 נחשב מובהק סטטיסטית.
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
רכישת מערכי נתונים GEO וזיהוי גנים הקשורים ל-UC
כדי לחקור שינויים בביטוי גנים ב-UC, התקבלו שני מערכי נתונים טרנסקריפטומיים - GSE87466 ו-GSE92415 - ממסד הנתונים GEO תוך שימוש במונח "קוליטיס כיבית" כמונח חיפוש, עם קריטריוני בחירה שכללו סוג מערך נתונים, אנוטציה ברורה של מחלה/ביקורת, וגודל מדגם מספק. מאגרי הנתונים הללו כללו 140 מקרים של UC ו-42 ביקורת בריאה. דגימות ללא תוויות פנוטיפ חד-משמעיות או פרופילי ביטוי תקפים הוצאו לפני האינטגרציה במורד הזר...
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
קוליטיס כיבית (UC) היא מחלה כרונית ואידיופתית דלקתית במעי, המשפיעה בעיקר על רירית המעי הגס, ומאופיינת בדלקת מתמשכת או אפיזודית. המחלה מפחיתה משמעותית את איכות חייו של המטופלים ומגבירה את הסיכון לפתח קרצינומה של המעי הגס. למרות שהשכיחות העולמית שלו עולה בהתמדה, האטיולוגיה המדויקת והמנגנונים הבסיסיים נותרו מובנים רק חלקית. אפשרויות הטיפול הנוכחיות לעיתים מניבות תוצאות לא אופטימליות ולרוב מלוויות בתופעות לוואי, מה שמדגיש את הצורך הקריטי בגישות טיפוליות יעילות יותר. מחקרים מצטברים מדגישים את...
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
המחברים מצהירים כי אין להם אינטרסים פיננסיים מתחרים ידועים או מערכות יחסים אישיות שיכולות להיראות כמשפיעות על העבודה המדוונת במאמר זה.
המחברים גם מודים למאגר הנתונים של GEO על גישה פתוחה לנתונים ולסוקרים על הערותיהם הבנות.
מימון: עבודה זו נתמכה על ידי הקרן הלאומית למדעי הטבע של סין (81973684, 82174358), קרן המדע הטבעי של מחוז סצ'ואן (2023NSFSC1760, 2025ZNSFSC1834), פרויקט מיוחד למחקר ברפואה סינית מסורתית של מנהל הרפואה הסינית המסורתית של מחוז סצ'ואן (2023MS365, 2020JC0094), וקרן החדשנות המשותפת של ועדת הבריאות של צ'נגדו ואוניברסיטת צ'נגדו לרפואה סינית מסורתית (WXLH202402019, WXLH202403012, WXLH202501238, WXLH202501201).
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.
| שם | חברה | מספר קטלוג | הערות |
|---|---|---|---|
| פוראריה לובאטה רדיקס | סצ'ואן נאוטוס לרפואה סינית מסורתית בע"מ. | NTC-001 | חומר צמחי להכנת GQD |
| Scutellariae Radix | סצ'ואן נאוטוס לרפואה סינית מסורתית בע"מ. | NTC-002 | חומר צמחי להכנת GQD |
| קופטידיס ריזומה | סצ'ואן נאוטוס לרפואה סינית מסורתית בע"מ. | NTC-003 | חומר צמחי להכנת GQD |
| Glycyrrhizae Radix et Rhizoma | סצ'ואן נאוטוס לרפואה סינית מסורתית בע"מ. | NTC-004 | חומר צמחי להכנת GQD |
| עכברים BALB/c, זכר, בן 8 שבועות | בייג'ינג HFK Bioscience Co., Ltd. | לא זמין | בעלי חיים ניסיוניים |
| דקסטרן סולפט נתרן (DSS), MW 36-50 kDa | ביומדיקל MP | 160110 | תגובת קוליטיס שמעוררת קוליטיס |
| נתרן קרבוקסימתיל צלולוז | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | C5678 | כלי לניתוח אוראלי |
| פאראפורמלדהיד | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | 158127 | קיבוע רקמות |
| המטוקסילין | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | H9627 | כתם גרעיני |
| אאוסין | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | E4009 | נגד-צבע להיסטולוגיה |
| מיקרוסקופ שדה בהיר | ניקון | אקליפס 80i | תצפית היסטולוגית |
| מערכת סריקת שקופיות | 3DHISTECH בע"מ. | פנורמה 250 | רכישת תמונות שקופית שלמה דיגיטלית |
| טריפסין | גיבקו, קרלסבד, קליפורניה, ארה"ב | 25200-056 | תגובת שחזור אנטיגנים |
| מי חמצן | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | H1009 | חסימת פעילות פרוקסידאז אנדוגנית |
| אלבומין סרום בקר (BSA) | סיגמא-אלדריץ', סנט לואיס, מיזורי, ארה"ב | A7906 | מדלל נוגדנים / מגיב חוסם |
| נוגדן ראשי אנטי-STAT3, ארנב | טכנולוגיית איתות סלולרי | 9139 | נוגדן ראשי ל-IHC/WB |
| נוגדן ראשוני נגד TIMP1, ארנב | אבקאם | AB211925 | נוגדן ראשי ל-IHC/WB |
| נוגדן ראשי נגד SPHK2, ארנב | אבקאם | ab187030 | נוגדן ראשי ל-IHC/WB |
| נוגדן ראשוני נגד HSPA5, ארנב | טכנולוגיית איתות סלולרי | 3177 | נוגדן ראשי ל-IHC/WB |
| נוגדן ראשוני נגד GAPDH, ארנב | אבקאם | ab9485 | בקרת טעינה עבור WB |
| נוגדן משני נגד ארנב עז | Yeasen Biotechnology | 33201 | נוגדן משני ל-IHC/WB |
| מגיב גילוי כרומוגני | Beyotime Biotechnology [Shanghai] Co., Ltd., שנחאי, סין | P0202 | התפתחות אותות לאימונוהיסטוכימיה |
| תוכנת ניתוח תמונות | סייברנטיקה של מדיה | Image-Pro Plus 7.0 | כימות של צביעה ועוצמת רצועה |
| תגובת חילוץ RNA מבוססת פנול | Invitrogen | ריאגנט טריזול, 15596026 | חילוץ RNA כולל |
| ערכת סינתזת cDNA לגדיל ראשון | ביוטכנולוגיה ביוטיים | DRR047A | תעתיק הפוך |
| מערכת PCR בזמן אמת | מערכות ביולוגיות יישומיות | ABI פריזמה 7300 | הגברה כמותית של PCR |
| תערובת SYBR Green qPCR | ביוטכנולוגיה ביוטיים | DRR420A | מגיב גילוי qPCR |
| תוכנת עיצוב ראשוני | Premier Biosoft International, פאלו אלטו, קליפורניה, ארה"ב | פריימר 5.0 | עיצוב פריימר |
| מאגר ניתוק RIPA | ביוטכנולוגיה ביוטיים | P0013B | תגובת חילוץ חלבון |
| מעכב פרוטאז (PMSF) | ביוטכנולוגיה ביוטיים | ST506 | הגנה מפני חלבון במהלך הליזה |
| ערכת בדיקת חלבון BCA | ביוטכנולוגיה ביוטיים | P0010 | כימות חלבונים |
| ממברנת PVDF | מרק מיליפור | IPVH00010 | ממברנת העברת חלבון |
| מחזור / מערכת גילוי כימילומינסצנטי | תרמו פישר סיינטיפיק | 32106 | זיהוי אות עבור כתם מערבי |
| תוכנה סטטיסטית | IBM | SPSS 23.0 | ניתוח סטטיסטי |
הגישה מוגבלת. התחברו או התחילו תקופת ניסיון כדי לצפות בתוכן זה.