ניתוח ביטוי במאגרי נתונים ציבוריים וברקמות קליניות מזווגות
דפוס הביטוי של miR-192 נבדק מחדש באמצעות נתוני miRNA-seq של TCGA-COAD ו-TCGA-READ שהופקו מ-GDC (איור 1). בקבצים אלו, הרשומה ברמת מסד הנתונים סומנה כ-hsa-mir-192 ולא ישירות כמין ה-hsa-miR-192-5p הבוגר; לפיכך, תוצאה זו ממסד הנתונים הציבורי פורשה כאות miR-192 של TCGA. בקוהורט המשולב של COADREAD, נותחו 616 גידולי ראשוניים ו-11 דגימות רקמה נורמליות מוצקות. דגימות הגידולים הראו ביטוי גבוה יותר של miR-192 מאשר הדגימות הנורמליות (חציונה 16.05 לעומת 13.46 log2[RPM + 1], מבחן Mann-Whitney p = 9.96 x 10-8). אותה מגמה של ביטוי גבוה יותר בגידולים נצפתה ב-COAD לבדו (455 גידולים ו-8 דגימות נורמליות, p = 9.89 x 10-8) וב-READ לבדו (161 גידולים ו-3 דגימות נורמליות, p = 1.94 x 10-5). ניתוח TCGA מזווג הראה גם הוא ביטוי גבוה יותר בגידולים ב-11 מקרים תואמים (הפרש חציוני בין גידול לנורמה = 3.13 log2[RPM + 1], מבחן Wilcoxon p = 0.0029). תוצאות אלו נותרו סותרות את תוצאות ה-qRT-PCR המזווגות המקומיות, שתמכו בביטוי מופחת של miR-192-5p בוגר ברקמות CRC (איור 2).
ביצועים אבחנתיים וקשרים קליניקופתולוגיים
ניתוח ROC מבוסס TCGA חזר על עצמו תוך שימוש בהגדרת ההבחנה בין דגימות גידול לדגימות נורמליות שתוארה לעיל (איור 3). ערך ה-AUC היה 0.968 עבור קוהורט COADREAD המשולב, 0.947 עבור COAD לבדו, ו-0.994 עבור READ לבדו. דגימה חוזרת מאוזנת נוספת הניבה ערכי חציוני AUC של 0.983 עבור COADREAD, 0.984 עבור COAD, ו-1.000 עבור READ; עם זאת, רווחי הסמך של 95% היו רחבים עבור COADREAD (0.868–1.000) ועבור COAD (0.781–1.000), וקבוצת READ כללה רק שלוש דגימות נורמליות. ניתוח "השאר-אחד-בחוץ" (Leave-one-out) של הדגימות הנורמליות הראה גם הוא כי הערכת ה-ROC הייתה רגישה לקבוצה הנורמלית הקטנה. לפיכך, תוצאת ה-ROC נשמרה רק כראיה גששית להפרדה בין גידול לנורמה, ולא כראיה מספקת עבור סמן ביולוגי לאבחון קליני. ניתוח קשר קלינו-פתולוגי המבוסס על ביטוי חציוני של miR-192 ב-TCGA הראה קשרים מובהקים עם שלב N (p = 0.019) וגיל (p = 0.004), בעוד ששלב T (p = 0.110) ושלב M (p = 0.932) לא היו מובהקים (טבלה 1).
ניתוח הישרדות
ניתוח הישרדות כוללת לא תמך בקשר פרוגנוסטי משמעותי עבור ביטוי של TCGA miR-192 ב-CRC (איור 3). בקוהורט המשולב של COADREAD, השוואת קפלן-מייר לפי חציונה (median-cut) לא הייתה מובהקת (log-rank p = 0.857), וגם רגרסיית Cox חד-משתנית לא הראתה קשר מובהק (HR = 0.934 לכל יחידת ביטוי log2, p = 0.374). אותה מסקנה נצפתה ב-COAD (log-rank p = 0.889; Cox p = 0.493) וב-READ (log-rank p = 0.807; Cox p = 0.716). לפיכך, הנתונים אינם תומכים בטענה ל-miR-192 כסמן פרוגנוסטי חזק ב-CRC.
גנים מטרים מועמדים והעשרה תפקודית
ניתוח של גנים מטרת-מועמדים בוצע על ידי שילוב של תחזיות miRWalk עם גנים בעלי ביטוי דיפרנציאלי ממאגרי הנתונים של GEO המסומנים כ-GSE89076 ו-GSE156355. סינון זה זיהה 93 גנים מועמדים לניתוחי העשרה ורשתות בהמשך (איור 4). ניתוח אינטראקציות בין חלבונים השאיר את KIF20A, TPX2, CDCA5, CCNB1, CDK1, PLP1, NRXN1, GRIK3, ו-KIF5C כגנים מרכזיים (hub genes) למחקר (איור 5). ניתוח רגישות של RNA-seq מ-TCGA הראה כי KIF20A, TPX2, CDCA5, CCNB1, ו-CDK1 היו בעלי ביטוי גבוה יותר ברקמת גידול מאשר ברקמה נורמלית, בעוד ש-PLP1, NRXN1, GRIK3, ו-KIF5C היו בעלי ביטוי נמוך יותר ברקמת הגידול. ניתוח מתאם עם ביטוי miR-192 מ-TCGA היה מעורב: CDCA5, CCNB1, ו-CDK1 הראו מתאמים חיוביים, בעוד ש-PLP1, NRXN1, ו-KIF5C הראו מתאמים שליליים. ניתוח העשרה הדגיש מונחים מייצגים הקשורים למטבוליזם של תרופות על ידי cytochrome P450, מטבוליזם של קסנוביוטיקה על ידי cytochrome P450, מטבוליזם של רטינול, מונחים הקשורים לסינפסות גלוטמטרגית וארגון פרה-סינפטי (איור 6). ממצאים אלו תומכים ברלוונטיות הביולוגית של רשת הגנים המרכזיים אך אינם מבססים Targeting ישיר של miR-192-5p.
בדיקת נדידת תאים
תאי HT29 הועברו טרנספקציה עם mimic של miR-192 כדי להעריך את ההשפעות הפונקציונליות של ביטוי יתר (upregulation) של miR-192. מיקרוסקופיה פלואורסצנטית אישרה טרנספקציה מוצלחת בקבוצות הטרנספקציה של miR-192-5p mimic ושל ביקרת השלילה (איור 7). בבדיקת ריפוי פצע (wound-healing assay), קצב הנדידה של תאי HT29 שעברו טרנספקציה של miR-192 היה נמוך יותר מזה של תאי HT29 שלא עברו טרנספקציה ומזה של תאים שעברו טרנספקציה של ביקרת השלילה, הן ב-24 h והן ב-48 h (איור 8). ב-24 h, קצב הנדידה היה 9.75 +/- 2.43% בקבוצה שעברה טרנספקציה של miR-192, בהשוואה ל-15.69 +/- 3.47% בקבוצה שלא עברה טרנספקציה ו-15.45 +/- 3.92% בקבוצת ביקרת השלילה (P = 0.010). ב-48 h, קצב הנדידה היה 20.04 +/- 2.54% בקבוצה שעברה טרנספקציה של miR-192, בהשוואה ל-26.84 +/- 7.65% ו-27.98 +/- 6.78% בשתי קבוצות הביקורת, בהתאמה (P = 0.021). תוצאות אלו מצביעות על כך שביטוי יתר של miR-192 מדכא את נדידת תאי HT29 in vitro.
ניתוח מתאם של חתימה חיסונית
הקשר בין ביטוי miR-192 ב-TCGA לבין חתימות הקשורות למערכת החיסון נבחן מחדש באמצעות נתוני RNA-seq של גידולים ראשוניים (איור 9). סך של 425 דגימות גידול עם נתוני miRNA ו-RNA-seq תואמים נכללו בניתוח המתאם של ציון החתימה. ביטוי miR-192 הראה מתאמים שליליים חלשים עם מספר חתימות של תאי חיסון, כולל תאי Tfh (rho = -0.205, p = 2.15e-05), תאי דנדריטיים מופעלים (rho = -0.186, p = 0.00011), תאי Treg (rho = -0.185, p = 0.00013), תאי דנדריטיים פלסמציטואידיים (rho = -0.169, p = 0.00046), מקרופאגים (rho = -0.149, p = 0.0020), תאי Tcm (rho = -0.136, p = 0.0050), תאי NK CD56bright (rho = -0.132, p = 0.0065), ותאי Tem (rho = -0.129, p = 0.0078). מתאמים חיוביים חלשים נצפו עבור חתימות של אוזינופילים (rho = 0.139, p = 0.0040) ותאי Th17 (rho = 0.121, p = 0.0129). לפיכך, נראה כי הקשר החיסוני של miR-192 ב-CRC הוא הטרוגני ולא следует לתארו אותו כחיובי באופן אחיד בקרב תסנני מערכת החיסון.
זמינות נתונים
מאגרי הנתונים הציבוריים שנבדקו במחקר זה זמינים בפורטל GDC/TCGA ובמאגר הנתונים Gene Expression Omnibus (GEO) תחת מספרי גישה GSE89076 ו- GSE156355ניתוחי המעבדה היבשה (dry-lab) והניתוחים החוזרים שבוצעו במהלך מחקר זה מוצגים ב- טבלה נלווית 1הנתונים הניסויים הבסיסיים התומכים בניתוחי ה-qRT-PCR ובריפוי הפצעים, כולל ערכי Ct גולמיים של qRT-PCR, מדידות של רוחב הפצע, תמונות מיקרוסקופיות מקוריות ונתונים ניסויים קשורים, מופיעים ב- תיקייה משלימה 1.

איור 1: סקירת ביטוי רב-סרטני (Pan-cancer) של miR-192 ב-TCGA. ציר ה-x מציג קיצורי סוגי גידולים, וציר ה-y מציג את הביטוי היחסי של miR-192. כחול מציין רקמה נורמלית, ואדום מציין רקמת גידול. מובהקות סטטיסטית מסומנת באופן הבא: * = p < 0.05, ** = p < 0.01, *** = p < 0.001; NS = לא מובהק. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.

איור 2: ביטוי של miR-192/miR-192-5p במאגרי נתונים של CRC וברקמות זוגיות. (A) ניתוח miRNA-seq של TCGA/GDC המראה אות hsa-mir-192 גבוה יותר ברמת מאגר הנתונים בדגימות גידול של CRC מאשר בדגימות נורמליות. (B) ניתוח qRT-PCR של דגימות רקמה מקומיות זוגיות המראה ביטוי נמוך יותר של miR-192-5p בוגר ברקמות CRC מאשר ברקמות סמוכות שאינן סרטניות. אנא לחצו כאן כדי לצפות בגרסה מורחבת של איור זה.

איור 3: ניתוחי אבחנה והישרדות גששים המבוססים על ביטוי hsa-mir-192 ב-TCGA. (A–C) עקומות הישרדות Kaplan-Meier. ניתוח ההישרדות לא תמך בקשר פרוגנוסטי מובהק. (D) עקומת ROC להבחנה בין גידול לרקמה נורמלית בקוהורט ה-TCGA. יש לפרש את ה-AUC בזהירות בשל המספר הקטן של דגימות נורמליות. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של איור זה.

איור 4: סריקת גנים מועמדים כמטרה וניתוח אינטראקציית חלבונים. (A) דיאגרמת ון המראה את החפיפה בין חיזויי מטרות של miRNA לבין גנים עם ביטוי דיפרנציאלי מ-GSE89076 ו-GSE156355. (B) רשת אינטראקציית חלבון-חלבון של 93 הגנים המועמדים. הצמתים מייצגים גנים, והקשתות מייצגות קשרי גנים חזויים או מאושרים. אנא לחצו כאן כדי לצפות בגרסה מורחבת של איור זה.

איור 5: סריקת גנים מרכזיים (Hub-genes) מרשת מטרות המועמדות. (A,B) מודולי רשת שזוהו באמצעות סריקת מודולים. (C) דירוג גנים מרכזיים על בסיס קישוריות הרשת. (D) מטרות מרכזיות מועמדות שנשמרו לאחר שילוב הסריקה מבוססת המודולים והדרגות. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.

איור 6: מונחי העשרה נבחרים של GO ו-KEGG עבור גנים מועמדים ומודולים משמעותיים ברשת. תרשים הבועות מסכם תפקודי ביולוגיה ומסלולים מועשרים מייצגים. תוצאות ההעשרה הן גששיות ודורשות תיקוף ניסויי. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.

איור 7: תמונות טרנספקציה של תאי HT29. (A) תאי HT29 מייצגים. (B) תמונת פלואורסצנציה של תאי HT29 שעברו טרנספקציה עם miR-192-5p mimic. (C) תמונת פלואורסצנציה של תאי HT29 שעברו טרנספקציה עם ביקורת שלילית con238. Scale bar = 100 µm. אנא לחצו כאן כדי לצפות בגרסה מורחבת של איור זה.

איור 8: בדיקת ריפוי פצע בתאי HT29. (A) תמונות שריטה נבחרות בנקודות זמן של 0 h, 24 h, ו-48 h בתאי HT29 שלא עברו טרנספקציה, תאי HT29 שעברו טרנספקציה ב-miR-192-5p mimic, ותאי HT29 שעברו טרנספקציה ב-con238. (B) כימות של קצבי ריפוי הפצע ב-24 h ו-48 h, מוצגים כאחוזים. סרגל קנה מידה = 100 µm. אנא לחצו כאן כדי לצפות בגרסה מורחבת של איור זה.

איור 9: מתאם בין ביטוי miR-192 ב-TCGA לבין ציוני חתימות הקשורות למערכת החיסון בדגימות גידולי CRC. הניתוח הראה קשרים חיסוניים הטרוגניים, עם מספר מתאמים שליליים חלשים, ובעיקר עבור חתימות של eosinophil ותאי Th17, מתאמים חיוביים. אנא לחץ כאן כדי להציג גרסה מוגדלת של איור זה.
| מאפיינים | ביטוי נמוך של TCGA miR-192 | ביטוי גבוה של TCGA miR-192 | ערך p | שיטה |
| n | 308 | 308 | | |
| שלב T, n (%) | 0.11 | Chisq.test |
| T1 | 9 (1.5%) | 11 (1.8%) | | |
| T2 | 44 (7.2%) | 59 (9.6%) | | |
| T3 | 209 (34.1%) | 211 (34.4%) | | |
| T4 | 43 (7%) | 27 (4.4%) | | |
| שלב N, n (%) | 0.019 | Chisq.test |
| N0 | 156 (25.5%) | 190 (31%) | | |
| N1 | 88 (14.4%) | 63 (10.3%) | | |
| N2 | 61 (10%) | 54 (8.8%) | | |
| שלב M, n (%) | 0.932 | Chisq.test |
| M0 | 214 (39.9%) | 234 (43.7%) | | |
| M1 | 41 (7.6%) | 47 (8.8%) | | |
| גיל, חציון (IQR) | 66 (56, 75) | 69 (60, 77) | 0.004 | Wilcoxon |
טבלה 1: הקשר בין קבוצות הביטוי של TCGA miR-192 לבין מאפיינים קליניקופתולוגיים ב-CRC. קבוצות ביטוי גבוה ונמוך הוגדרו על פי רמת הביטוי החציונית של TCGA miR-192. בטבלה המסכמת, שלב N והגיל היו משמעותיים, בעוד ששלב T ושלב M לא היו משמעותיים.
טבלה נלווית 1: תוצאות ניתוח חוזר של נתוני מעבדה יבשה שנוצרו עבור מחקר זה. חוברת העבודה כוללת השוואות ביטוי של TCGA/GDC, ניתוחי רגישות ROC מאוזנים, ניתוחי leave-one-normal-out, ניתוחי קשר קליני-פתולוגי, ניתוחי הישרדות, סיכומי ביטוי ומתאם של גני מפתח (hub-genes), ותוצאות מתאם של חתימות חיסוניות.אנא לחץ כאן להורדת הקובץ.
תיקיית השלמה 1: נתונים ניסויים גולמיים התומכים בניסויי ה-qRT-PCR, הטרנספקציה של התאים וריפוי הפצעים. קבצים אלו כוללים ערכי Ct מקוריים של qRT-PCR, נתוני ציטומטריה של זרימה/טרנספקציה, מדידות רוחב פצע ותמונות מיקרוסקופיה מייצגות ששימשו לאנליזות המוצגות במאמר.אנא לחצו כאן להורדת קובץ זה.