כל ההליכים הכוללים רקמות אנושיות עמדו בהנחיות המוסדיות ובהצהרת הלסינקי ואושרו על ידי ועדת הביקורת המוסדית של אוניברסיטת פוג'יאן לרפואה (אישור מס' 2021KYB089). התקבלה הסכמה מדעת בכתב מכל המשתתפים לפני רכישת הרקמות.
ביטוי גנים וניתוח הישרדות
נתוני ריצוף ה-RNA והמידע הקליני המתאים התקבלו ממספר מאגרי מידע ציבוריים. 1) קוהורט TCGA: נתוני RNA-seq (FPKM) עבור 175 דגימות גליובלסטומה מולטיפורמה (GBM) ו-534 דגימות גליומה בדרגה נמוכה (LGG) הורדו מאטלס הגנום של הסרטן (https://portal.gdc.cancer.gov/); 2) ביקורת תקינה: פרופילי ביטוי של 211 רקמות מוח תקינות ו-662 רקמות גליומה הורדו ממאגר UCSC Xena (https://xenabrowser.net/datapages/); 3) אימות חיצוני: נתונים מקבוצות CGGA693 ו-CGGA325 התקבלו מאטלס הגנום הסיני לגליומה (http://www.cgga.org.cn); 4) מערך נתונים GEO: מערך הנתונים GSE43378, הכולל נתוני ביטוי וקליניים עבור 50 דגימות גליומה, הורד מ-Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/). כל נתוני הספירה הגולמיים הומרו לתמלילים למיליון (TPM) והוסבו ללוג2. עבור מערכי נתונים שכבר נורמלו, נבדקו מטריצות ביטוי כדי להבטיח התפלגויות דומות. גנים עם ערכי TPM < 1 ביותר מ-80% מהדגימות הוצאו. מידע קליני חסר (גיל, מצב IDH, קודלציה של 1p/19q, מתילציה של MGMT) הוסר באמצעות סינון מקרים מלאים. אפקטי האצווה בין מערכי הנתונים הותאמו באמצעות אלגוריתם ComBat שמיושם בחבילת R sva. ערכי הביטויים תוקנו על ידי טרנספורמציית ציון z בכל מערך נתונים. ניתוחי הישרדות בוצעו באמצעות חבילות R של Survival ו-survminer. המטופלים חולקו לקבוצות ביטוי גבוה ונמוך בהתאם לרמת הביטוי החציונית של IRAIN. נוצרו עקומות הישרדות של קפלן-מאייר, והמשמעות הסטטיסטית הוערכה באמצעות מבחן ה-log-rank. יחסי סיכון (HRs) ורווחי ביטחון של 95% (CIs) הוערכו באמצעות מודלים של רגרסיית סיכונים פרופורציונליים של Cox.
הגדרת קבוצות גנים חיסוניות ומטבוליות
גנים הקשורים למערכת החיסון (IRGs, n = 2,483) התקבלו ממאגר ImmPort (https://www.immport.org/shared/), וגנים הקשורים למטבוליה (MRGs, n = 948) התקבלו ממאגר החתימות המולקולריות (MSigDB, https://www.gsea-msigdb.org/). הסט המשולב של גנים אלו הוגדר כגנים הקשורים לאימונומטבוליה (IMRGs). רשימות גנים אלו שימשו כמקורות לביטוי שונה וניתוחי רשתות מאוחרים יותר.
ביטוי שונה וניתוח רשת קו-ביטוי של גנים משוקללים
גנים מבוטאים באופן שונה (DEGs) בין רקמות המוח והגליומה התקינות זוהו באמצעות חבילת R limma. נתוני הביטוי הותקנו במודל ליניארי ואחריו מתון אמפירי של בייס. גנים עם שינוי קפל |log₂| > 1.5 ושיעור גילוי שגוי (FDR) < 0.05 נחשבו כביטויים באופן מובהק באופן שונה. ניתוח רשת ביטוי משותף של גנים משוקלל (WGCNA) בוצע באמצעות חבילת R WGCNA. מדגמים חריגים הוצאו באמצעות קיבוץ היררכי. עוצמת סף הרכות הוגדר ל-β = 8 כדי להשיג אינדקס התאמה טופולוגי ללא קנה מידה (R2 ≥ 0.85) תוך שמירה על קישוריות ממוצעת מספקת. נבנו מטריצות חפיפה טופולוגית (TOM), והגנים קובצו למודולים בגודל מינימלי של 50 באמצעות אלגוריתם חיתוך עץ דינמי. הגנים העצמיים של המודול היו קשורים לתכונות קליניות, והמודול המזוהה ביותר עם גליומה (Pearson's r > 0.7, P < 1×10-10) נבחר לזיהוי גן המרכזי.
בניית מודלים פרוגנוסטיים מבוססי למידת מכונה
מסגרת LOOCV מקיפה הכוללת עשרה אלגוריתמים של למידת מכונה יושמה לבניית והערכת מודלים פרוגנוסטיים. בסך הכל, 101 תהליכי עבודה קומבינטוריים יושמו באמצעות קבוצת TCGA כמערך הנתונים להדרכה. גנים אימונוממטבוליים הקשורים לפרוגנוזה (IMRGs) זוהו לראשונה באמצעות רגרסיה חד-משתנית של Cox (P < 0.05). המודל האופטימלי נקבע על ידי מקסום מדד הקונקורדנס הממוצע של הארל (מדד C) בשלושה מערכי נתונים לאימות (CGGA693, CGGA325 ו-GSE43378). מודל RSF-Enet שנוצר (α = 0.3) הראה את הביצועים החיזויים הגבוהים ביותר ושמר על הכללה איתנה בין קבוצותעצמאיות 22.
TME וחדירת מערכת חיסון
כדי לאפיין באופן מקיף את הנוף האימונוגנומי, השתמשנו בגישה אנליטית רב-שכבתית. ראשית, רמות החדירה החיסונית והסטרומלית נמדדו באמצעות אלגוריתם ESTIMATE23. ביטוי שונה של מולקולות מפתח בנקודות ביקורת חיסוניות, כולל PDCD1, CTLA4 ו-LAG3, הוערך לאחר מכן באמצעות ניתוח מבוסס לימה, והקורלציות בין גני נקודות הביקורת הוצגו באמצעות מטריצות קורלציה. פרופילי מוטציות סומטיות מ-903 דגימות גליומה בקבוצת TCGA שימשו לחישוב ציוני עומס מוטציה של גידול (TMB), אי-יציבות מיקרו-סווייט (MSI) ודיספונקציה והדרה חיסונית של גידול (TIDE) כדי לחזות תגובות פוטנציאליות לאימונותרפיה. לאחר מכן חולקו המטופלים לארבע קבוצות פרוגנוסטיות לפי ציון משולב של מצב TMB (גבוה/נמוך) וסיכון (גבוה/נמוך), והושוו תוצאות הישרדות באמצעות ניתוח קפלן-מאייר.
ניתוח העשרה פונקציונלית
אונטולוגיית גנים (GO) וניתוחי העשרת מסלולים של אנציקלופדיית קיוטו לגנים וגנומים (KEGG) נערכו באמצעות clusterProfiler של חבילת R. תוצאות העשרה עם ערכי P מותאמים < 0.05 נחשבו מובהקות סטטיסטית. תהליכים ביולוגיים, רכיבים תאיים ותפקודים מולקולריים מיוצגים ביתר על המידה הוצגו באמצעות גרפים של נקודות וגרפים בר. רשתות אינטראקציה בין חלבון לחלבון (PPI) נבנו באמצעות מסד הנתונים STRING (≥ 0.4) והוצגו ב-Cytoscape. מודולים פונקציונליים בתוך רשת ה-PPI זוהו באמצעות אלגוריתם MCODE. אינטראקציה בין גן לגן ורשתות ביטוי משותף נותחו בהמשך באמצעות GeneMANIA (https://string-db.org; ציון ביטחון ≥ 0.4) והוצגו ב-Cytoscape. מודולים פונקציונליים בתוך רשת ה-PPI זוהו באמצעות אלגוריתם MCODE. אינטראקציה בין גנים לרשתות ביטוי משותף וביטוי משותף נותחו בהמשך באמצעות GeneMANIA (https://genemania.org), שמשלב מידע על אינטראקציות פיזיות וגנטיות, מסלולים משותפים ודפוסי ביטוי משותף כדי להסיק אסוציאציות פונקציונליות פוטנציאליות.
דגימות קליניות
רקמות גליומה טריות (n = 6) ורקמות מוח לא גידוליות סמוכות (n = 6; נמצאות לפחות 3 ס"מ משפת הגידול ואושרו היסטולוגית כנטולות גידול) נאספו ממטופלים שעברו כריתה ראשונית של גליומה בבית החולים המסונפ ג'אנגג'ואו של אוניברסיטת פוג'יאן לרפואה. אף אחד מהמטופלים לא קיבל כימותרפיה או רדיותרפיה לפני הניתוח. כל האבחנות הפתולוגיות אושרו באופן עצמאי על ידי שני נוירופתולוגים לפי סיווג ארגון הבריאות העולמי (WHO) לשנת 2021 של גידולים במערכת העצבים המרכזית. מיד לאחר כריתה כירורגית, דגימות רקמה נשטפו במלח קר עם פוספט קר (PBS) להסרת שאריית דם, הוקפאו במהירות בחנקן נוזלי (-196°C), ונשמרו בטמפרטורה של -80°C עד להוצאת ה-RNA.
קווי תאים ותרבות תאים
קווי תאי גליובלסטומה אנושיים SHG44, U251, A172 ו-T98G, וכן תאי גליה אנושיים רגילים (HEB), התקבלו ממאגרים מאומתים ואושרו כי הם נקיים מזיהום מיקופלזמה לפני השימוש. התאים נשמרו במדיום Modified Eagle's Medium (DMEM, גלוקוז גבוה) של Dulbecco, עם תוספת של 10% סרום בקר עוברי (FBS), 2 מ"מ L-גלוטמין ו-1% פניצילין-סטרפטומיצין, בטמפרטורה של 37 °C באינקובטור לח עם 5% CO₂. התאים עברו כל 4-5 ימים עם הגעתם ל-80-90% מפגש מתאים. כדי לקבוע קווי תאים שמביעים ובקרת יתר של IRAIN, תאים הועברו עם וקטורים לנטיויראליים הנושאים את תמלול IRAIN באורך מלא או וקטור ריק כביקורת. שיבוטים יציבים נבחרו באמצעות פורומיצין (2 מיקרוגרם/מ"ל) למשך 14 ימים. יעילות הבעה מופרזת אושרה על ידי PCR שעתוק הפוך כמותי (qRT-PCR) לפני הבדיקות במורד הזרם.
3- (4,5-דימתילתיאזול-2-yl)-2,5-דיפנילטרזוליום ברומיד (MTT)
התאים זרעו בלוחות של 96 בארות בצפיפות של 1 ×-104 תאים לבאר ב-100 מיקרו ליטר של מדיום תרבית מלא. לאחר 24, 48 ו-72 שעות לאחר הזרעה, נוספו לכל באר 20 מיקרוליטר של תמיסת MTT (5 מ"ג/מ"ל במי מלח מפוזר פוספט) ודגרו במשך 4 שעות בטמפרטורה של 37 מעלות צלזיוס. הסופרנטנט הוסר, והוספו 150 מיקרוליטר דימתיל סולפוקסיד (DMSO) כדי להמיס את גבישי הפורמזן. הצלחת נערכה בעדינות במשך 10 דקות כדי להבטיח מסיסות מלאה. הספיגה נמדדה ב-490 ננומטר באמצעות ספקטרופוטומטר מיקרופלט. קריאות רקע מבארות ריקות הוחסרו. חיאום התאים חושב ביחס לקבוצת הביקורת של 24 שעות (נקבע כ-1.0). כל הניסויים בוצעו עם שישה שכפולים טכניים ושלושה שכפולים ביולוגיים עצמאיים. הנתונים מבוטאים כסטיית תקן ממוצעת ± (SD), והמובהקות הסטטיסטית נקבעה באמצעות מבחן t דו-זנב.
ציטומטריית זרימה לאפופטוזיס (Annexin V - FITC/PI TINING)
התאים זרעו בנקודת מפגש של 60-70% וטופלו במשך 24 שעות בתנאים המצוינים. תאים צפים ותאים דביקים נאספו באמצעות טריפסין ללא EDTA, שולבו ונשטפו פעמיים עם PBS קר כקרח. כדורי התאים הושעו מחדש בבופר קשירת Annexin V (10 mM HEPES pH 7.4, 140 mM NaCl, 2.5 mM CaCl2) ב-1 ×10 6 תאים למ"ל. לכל דגימה, 100 מיקרוליטר של תלייה הודגרו יחד עם 5 מיקרוליטר של Annexin V-FITC ו-5 מיקרוליטר יוד פרופידיום (PI; 50 מיקרוגרם/מ"ל) בחושך למשך 15 דקות בטמפרטורת החדר. לאחר הוספת 400 מיקרוליטר של בופר קישור, הדגימות נשמרו על קרח ונותחו תוך שעה על ציטומטר זרימה (488 ננומטר; 530/30 ננומטר ל-FITC ו->585 ננומטר ל-PI). נכללו ביקורות צבע יחיד וביקורת פלואורסצנציה מינוס אחת מתאימות לפיצוי. לפחות 10,000 אירועים נרשמו לכל דגימה. הנתונים נותחו באמצעות גייטינג רבעי: אוכלוסיות חיות (Annexin V⁻/PI⁻), אפופטוטיות מוקדמת (Annexin V⁺/PI⁻), אפופטוטיות מאוחרות (Annexin V⁺/PI⁺), ונמקיות (Annexin V⁻/PI⁺). דווחו אחוזים של תאים אפופטוטיים מוקדמים + מאוחרים (ממוצע ± SD, n = 3).
PCR כמותי בזמן אמת (qRT-PCR)
סך ה-RNA בודד באמצעות תגובת חומצה פנול-גואנידיניום בהתאם לפרוטוקול היצרן. טוהר ה-RNA אושר באמצעות ספקטרופוטומטריה (A₂₆₀/A₂₈₀ = 1.8-2.1), והשלמות אושרה על ידי אלקטרופורזה ג'ל (מספר שלמות RNA ≥ 7). מיקרוגרם אחד של RNA כולל טופל ב-DNase I והועתק לאחור בתגובה של 20 מיקרוליטר באמצעות הקסמרים אקראיים ופריימרים אוליגו(dT). התגובה בוצעה בטמפרטורה של 25 מעלות צלזיוס במשך 10 דקות, 50 מעלות צלזיוס למשך 30 דקות, ו-85 מעלות צלזיוס למשך 5 דקות. PCR כמותי בוצע במערכת של 10 מיקרוליטר שכללה 5 מיקרוליטר של 2× SYBR Green Master Mix, 0.3 מיקרומולר לכל פריימר, ו-1 מיקרוליטר cDNA (≈ 20 ng RNA שקול). תנאי המחזור התרמי היו 95 מעלות צלזיוס למשך 5 דקות, ואחריהם 40 מחזורים של 95 מעלות צלזיוס למשך 15 שניות ו-60 מעלות צלזיוס למשך 30 שניות, ולאחר מכן ניתוח עקומת התכה מ-65 מעלות צלזיוס עד 95 מעלות צלזיוס במרווחים של 0.3 מעלות צלזיוס. כל התגובות בוצעו בשלוש כפילות, יחד עם בקרות ללא תבנית ו-RT מינוס. ערכי Ct > 35 או שכפול טכני של SD > 0.5 הוצאו. הביטוי היחסי חושב באמצעות שיטת 2⁻ΔΔCt, כאשר GAPDH משמש כבקרה פנימית. דווחו ערכי ממוצע ± SD משלושה שכפולים ביולוגיים עצמאיים, והבדלים בקבוצות נותחו באמצעות מבחן t דו-זנב.
ניתוח כתמים מערביים
התאים עברו פירוק על קרח בבופר RIPA (50 מ"מ טריס-HCl, pH 7.4, 150 מ"מ NaCl, 1% NP-40, 0.5% נתרן דאוקסיכולאט, 0.1% SDS) בתוספת מעכבי פרוטאז ופוספטאז. הליזטים דוגרו במשך 30 דקות על קרח עם מערבולת לסירוגין ונוקרו בצנטריפוגה ב-12,000 × גרם למשך 15 דקות בטמפרטורה של 4 מעלות צלזיוס. ריכוזי החלבון נמדדו באמצעות מבחן BCA, הותאמו ל-1-2 מיקרוגרם/מיקרוליטר, וערבבו 1:3 עם 4× בופר Laemmli (בופר סופי 1× המכיל 100 mM DTT). הדגימות עברו שינוי טמפרטורה ב-95 מעלות צלזיוס למשך 5 דקות. כמויות שוות של חלבון (50 מיקרוגרם) נפתרו על ידי 12% SDS-PAGE ב-100 וולט במשך 90 דקות והועברו אלקטרו לממברנות PVDF ב-250 mA למשך 90 דקות. הממברנות נחסמו עם 5% חלב דל שומן ב-TBST (0.1% Tween-20) למשך שעה בטמפרטורת החדר (או 5% BSA לפוספופרוטאנים) ודגרו במהלך הלילה בטמפרטורה של 4°C עם נוגדנים ראשוניים נגד IGF1, IGF1R, JAK2, p-JAK2 (Y1007/1008), STAT3, p-STAT3 (Y705), BIRC5 ו-β-אקטין (דילול טיפוסי 1:1000, β-אקטין 1:5000). לאחר שלוש שטיפות של 10 דקות ב-TBST, הממברנות דגרו בנוגדן משני מצמיד ל-HRP (1:5000) למשך שעה בטמפרטורת החדר, נשטפו שוב ופותחו באמצעות מצע כימילומינצנטי. עוצמות הרצועות נמדדו באמצעות ImageJ, נורמל ל-β-אקטין או חלבון כולל, והובאו לביטוי כממוצע ± SD משלושה ניסויים עצמאיים.
אימונוציטוכימיה
תאים שגדלו על כיסויי זכוכית סטריליים נשטפו פעמיים עם PBS וקיבעו ב-4% פאראפורמלדהיד למשך 15 דקות בטמפרטורת החדר. לאחר שלוש שטיפות PBS, התאים חודרו עם 0.2% Triton X-100 למשך 10 דקות, נחסמו עם אלבומין סרום בקר (BSA) של 5% למשך שעה אחת, ודגרו במהלך הלילה בטמפרטורה של 4 מעלות צלזיוס עם נוגדן אנטי-CD31 ראשוני (דילול של 1:200 ב-1% BSA). לאחר שלוש שטיפות PBS, התאים הוגברו בנוגדן משני מצומד אלקסה פלואור (דילול 1:500) למשך שעה בחושך, צובעו ב-DAPI (1 מיקרוגרם/מ"ל, 5 דקות), והותקנו במדיום אנטיפייד. התמונות צולמו באמצעות מיקרוסקופ פלואורסצנציה באותן הגדרות חשיפה והגבר. אחוז השטח החיובי ל-CD31 נמדד בחמישה שדות לא חופפים שנבחרו באקראי לכל דגימה באמצעות תוכנת ImageJ. בדיקה זו בוצעה במודלים תאיים ולא בחתכי רקמה.
ניתוח סטטיסטי
בוצעו ניתוחים סטטיסטיים באמצעות גרסה 4.3.0 של R יחד עם החבילות הנלוות אליה. להשוואת משתנים קטגוריים, נעשה שימוש במבחן חי-בריבוע, בעוד שמשתנים רציפים הוערכו באמצעות מבחן סכום הדרגה וילקוקסון או מבחן T. הערכת משתנים רציפים בוצעה באמצעות מקדם המתאם של פירסון. ניתוחי הישרדות בוצעו באמצעות חבילת ההישרדות, שכללה מודלים של סיכונים פרופורציונליים של Cox ויצירת עקומות קפלן-מאייר, כאשר ספי הסטרטיפיקציה האופטימליים נקבעו על ידי חבילת הסורוויינר והנוסחה Riskscore =
. חבילת CompareC שימשה להערכת מדדי C של משתנים שונים. עקומת מאפייני ההפעלה של המקלט (ROC), שנועדה לחזות משתנים קטגוריים בינאריים, נוצרה באמצעות חבילת pROC. בנוסף, השטח התלוי בזמן תחת עקומת ROC (AUC) עבור מדדי הישרדות נותח באמצעות חבילת timeROC. כל המבחנים הסטטיסטיים בוצעו בגישה דו-צדדית. רמת מובהקות של P < 0.05 נחשבה מובהקת סטטיסטית.