$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
נתוני qRT-PCR ראשוניים הראו כי מוטציה EWS / FLI בשם DAF, עם טירוזין ספציפי למוטציות אלאנין באזור החוזר על עצמו ומופרע של EWS, שמר על היכולת להפעיל גנים היעד EWS / FLI, אבל לא הצליח להדחיק גנים היעד קריטי23. כדי להבין טוב יותר את הקשר בין שאריות אלה בתחום EWS לבין הפונקציה EWS/FLI, נעשה שימוש בפרוטוקול המתואר לעיל ומתואר באיור 1. A673 תאי סרקומה יואינג הועברו באופן ויראלי עם shRNA מיקוד 3'UTR של FLI1, וכתוצאה מכך דלדול של אנדוגני EWS / FLI. לאחר ארבעה ימים של בחירה, פונקציית EWS / FLI חולצה עם טרנסולציה ויראלית של מבנים מוטנטיים שונים מתויגים EWS / FLI, עם וקטור ריק כשליטה ללא הצלה. מוטציה לא פונקציונלית חסרה את תחום EWS, הנקרא Δ22, שימשה כשליטה שלילית וסוג פראי EWS/FLI, הנקרא wtEF, שימשה כשליטה חיובית(איור 2A). DAF שימש כמבנה הבדיקה, אם כי ניתן להשתמש ביותר ממבנה בדיקה אחד במידת הרצון. תאים נבחרו במשך 10 ימים נוספים כדי לאפשר ביטוי מבנה לייצב ולאחר מכן נאספו עבור RNA (עם שלב הסרת gDNA), חלבון ומושבה להרכיב assays. ארבעה העתקים נאספו וייצוג qRT-PCR וכתמים מערביים המציגים נוקאאוט והצלה יעילים מוצגים באיור 2B-D. יש לציין כי תאים שניצלו על ידי DAF לא הצליחו ליצור מושבות כפי שמוצג באיור 2E, דבר המצביע על טרנספורמציה אונקוגנית לקויה.
לאחר השלמת האימות המשכפל ובדיונים פנוטיפיים, RNA הוגש למכון לרפואה גנומית בבית החולים הלאומי לילדים להכנת הספרייה ורצף הדור הבא עם כ -50 מיליון קריאות בסוף זיווג 150 bp שנאספו. הנתונים הוחזרו כקבצים מהירים.gz. קריאות באיכות נמוכה קוצצו מקבצים אלה עם TrimGalore ו- STAR שימש ליישר קריאות לגנום האנושי hg19 ולספור את הקריאות לכל גן. hg19 שימש למטרות תאימות עם ערכות הנתונים האחרות שאוצר עבור EWS/FLI המשמשות בניתוח במורד הזרם. ספירות קריאה אלה שולבו במטריצת ספירה אחת עבור כל הדגימות, ש-6 השורות הראשונות שלהן מוצגות באיור 3.
ספירות בוצעו בתחילה דרך DESeq2 ללא נורמליזציה של אצווה, עם זאת, בדיקה חזותית של המרחק מדגם לדגימה הראתה אפקטי אצווה מבלבלים פוטנציאליים כפי שמוצג עם חצים אדומים באיור 4A. זה ככל הנראה התעורר עקב שונות ביולוגית שהוצגה על ידי המעבר של תאים בתרבית והבדלים בעיבוד של כל אצווה. נורמליזציה עבור אפקטי אצווה בוצעה עם ComBat ומומלץ בדרך כלל. המרחקים לדוגמה לדגימה של הנתונים מנורמלים באצווה מוצגים באיור 4B. לאחר נורמליזציה של אצווה, נעשה שימוש ב- DESeq2 ליצירת פרופילי תמלול עבור שלושת המבנים (wtEF, Δ22 ו- DAF) ביחס לקו הבסיס. שים לב כי בעוד תאי A673 "הורים" (חיקוי מדומה והצלה מדומה, הנקראים "iLuc" כאן) נכללו בניתוח הדיפרנציאלי, ההתייחסות לניסוי זה הם התאים עם תאי IEF מרוקנים EWS/FLI, הנקראים תאי iEF. פרופיל התמלול יכול להיווצר עבור חלבון אנדוגני כאן על ידי השוואת מדגם iLuc ל- iEF, וזה עשוי להיות עניין להבין איך מערכת ההצלה עובדת, אבל זה לא המטרה של ניתוח מסוים זה. פרופילי התמלול שנוצרו עבור המוטנטים כוללים פקדים חיוביים (wtEF) ושליליים (Δ22), ביחס ל- iEF, כך שאלה אמורים לתפקד כמדדים למוטנטים אחרים. זה חשוב, כמו השליטה החיובית בדוגמה זו לא לגמרי recapitulate הפונקציה של אנדוגני EWS / FLI כפי שנדון במקום אחר7,23.
ניתוח הרכיבים העיקרי (PCA) באיור 5 מצביע על כך שפרופיל התמלול של DAF הוא ביניים בין wtEF ל- Δ22, המאשר פונקציה חלקית. יתר על כן, קיבוץ היררכי של 1000 הגנים המשתנים ביותר על פני דגימות הראה כי DAF נכשל בדיכוי גנים היעד EWS / FLI, ורק באופן חלקי שמר על פעילות הפעלת גנים כפי שמוצג איור 6A ו איור S5. ניתוח ToppGene העלה כי מחלקות הגנים ש-DAF מפעילה שונות מבחינה תפקודית מיעדים המופעלים על-ידי EWS/FLI שבהם DAF אינו פונקציונלי (איור 6B). מעניין, הפונקציה של גנים מופעלים שניצלו על ידי wtEF, אבל לא על ידי DAF, נראה קשור בקרת תעתיק ויסות כרומטין. בהתבסס על התוצאות של התקפות היווצרות המושבה, הגנים מחתימת גן הליבה הזו צריכים להיות מנותחים עוד יותר על תפקידם באונקוגנזה בתיווך EWS / FLI. החשיבות של דיכוי גנים בתיווך EWS / FLI תוארה בעבר17.
ידוע כי EWS / FLI בעל זיקה מחייבת ייחודית עבור אלמנטים חוזרים GGAA-microsatellite19,22, וכי כריכה באלמנטים אלה מניעהבמורד הזרם רגולציה גנים 11,15,18,20,22. מיקרו-סטליטים אלה אופיינו כמזוהים עם הפעלה או דיכוי, וקרובים ל- TSS (< 5 kb) או דיסטליים ל- (> 5 kb) TSS25. בנוסף, ישנם גנים מווסתים EWS / FLI עם זיקה גבוהה (HA) מוטיבים ETS פרוקסימלי TSS23. על מנת לנתח עוד יותר את המאפיינים של פונקציית DAF ואת סוגי הגנים המופעלים על ידי EWS / FLI ש- DAF הצליח להציל, נותח ביטוי דיפרנציאלי של גנים הקשורים למעמדות שונים אלה. מעניין לציין ש-DAF הצליחה להציל גנים המופעלים על ידי GGAA-microsatellite, אך לא הצליחה להציל גנים מופעלים ליד אתר HA כפי שניתן לראות באיור 7. כפי שניתן לראות באשכולות היררכיים, DAF נכשל להציל דיכוי בתיווך EWS / FLI על פני מחלקות מוטיב. נתונים אלה מצביעים על כך ש- DAF שומר על תכונות מבניות מספיקות של EWS כדי לאגד ולהפעיל מ- GGAA-microsatellites, הן פרוקסימלי והן דיסטלי ל- TSS. זה נובע ככל הנראה מתחום SYGQ שלם שנחשב חשוב לפעילות EWS / FLI ב GGAA חוזר11. נתונים אלה מצביעים גם על כך שהטירוסינוסים הספציפיים שעברו מוטציה ב- DAF ממלאים תפקידים חשובים, אך לא מובנים היטב, ברגולציה גנטית בתיווך EWS / FLI מאתרי HA, כמו גם בדיכוי גנים, המדגישים תחום חשוב של חקירה נוספת.

איור 1: זרימת עבודה. תיאור של ההליך שלב אחר שלב לביצוע מיפוי פונקציית מבנה על-ידי transcriptomics. התאים היו מוכנים תחילה לבטא את חבילת המבנים הנדרשת למיפוי פונקציית מבנה. בעקבות הביטוי, תאים נקצרו עבור RNA וחלבון ונחשבו פנוטיפים קורלטיביים. ביטוי המבנים אומת, ותהליך זה חזר על עצמו 3-4 פעמים כדי לאסוף שכפולים ביולוגיים עצמאיים. לאחר מכן הוגש RNA לרצף הדור הבא (NGS). כאשר הנתונים התקבלו, הנתונים קוצצו עבור איכות, יישור וספירות לכל תמליל חושבו. אפקטי אצווה נשלטו עבור וחתימות תמלול וביטוי דיפרנציאלי נקבעו באמצעות DESeq2. ניתן לשלב קיבוץ באשכולות הירארכי וניתוח במורד הזרם המשלב ערכות נתונים אחרות של -omics וניתוח מסלול או פונקציונלי שונה. אנא לחץ כאן כדי להציג גירסה גדולה יותר של איור זה.

איור 2: אימות של ביטוי מבנה ובוחות קורלטיביות. (A)סכמטי המתאר את המבנים שנבדקו בדוגמה זו. (B)אימות של נוקאאוט של EWS / FLI אנדוגני וביטוי של מבנים מתויגים 3X-FLAG על ידי אימונובלוט. (C,D) אימות פעילות מבנה בגן יעד המופעל על-ידי EWS/FLI (C), NR0B1ו- (D)מודחק, TGFBR2, על-ידי qRT-PCR. הנתונים מוצגים כסטיית תקן ממוצעת +/. ערכי P חושבו במבחן המשמעות ההגון של טוקי. * p < 0.05, ** p < 0.01, *** p < 0.005 (E) מושבה סופרת מביקורות רכות אגר שבוצעו כדי להעריך את פעילות השינוי של מבנים. ערכי P חושבו במבחן המשמעות ההגון של טוקי. * p < 0.05, ** p < 0.01, *** p < 0.005. נתון זה מותאם מאת Theisen, et al.23אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

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

איור 4: מפת חום של מרחק מדגם לדגימה. (A)התוויית מרחק לדוגמה לדוגמה המציגה את קיבוץ הדוגמאות של נתוני הספירה הגולמית. דוגמאות המקבצות באשכולות הן לפי אצווה והן לפי דוגמה מסומנות בחצים אדומים. (B)התוויית מרחק לדוגמה לדוגמה לאחר נורמליזציה של אצווה עם ComBat. כאן, דגימות מכל המשכפלים מקבצות יחד, ללא תלות באצווה. אנא לחץ כאן כדי להציג גירסה גדולה יותר של איור זה.

איור 5: תוצאות ניתוח ביטוי דיפרנציאלי. (A)נתווה ניתוח רכיבים עקרוני (PCA) של החתימות התמלוליות שנוצרו עבור כל הדגימות מציגות קיבוץ תוך-מדגמי חזק וממחישות כי DAF מתווך בין הפקדים החיוביים (wtEF) ושליליים (Δ22). (B)חלקות הר געש המציגות את ה- -log (p-value) שרטטו כנגד log2FoldChange עבור גנים בכל מבנה. גנים בעלי ערך p מותאם < 0.05 ו- |log2(FoldChange)) | > 1 נחשבים משמעותיים ומוצגים באדום. לוח 5B מותאם מתייסן, ואח'23אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

איור 6: קיבוץ היררכי לזיהוי מחלקות גנים. (A)קיבוץ הירארכי של 1000 הגנים המשתנים ביותר בכל המבנים וקו הבסיס, iEF, מציג ש- DAF מציל חלקית הפעלת גנים בתיווך EWS/FLI. (B)אונטולוגיה גנטית (פונקציה מולקולרית) נובעת מ- ToppGene המציגה את ההעשרה התפקודית של גנים המופעלים על ידי EWS / FLI שניצלו או לא ניצלו על ידי DAF. לוח 6B מותאם מתייסן, ואח'23אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.

איור 7: ניתוח מפורט של רכיבי תגובה שונים של גורם תמלול למבנים שונים: (A) סכמטי המתאר את עיבוד הנתונים המשמש ליצירת לוחות (B) ו- (C) על-ידי שילוב ערכות נתונים זמינות אחרות עם הפרופילים התמלוליים כאן. (B,C) אוסף המציג את ההצלה של סוגים שונים של מטרות ישירות EWS / FLI- (B) מופעל ו (C) מודחק מטרות. גנים כלולים היו רק אותם גנים עם ביטוי דיפרנציאלי ניתן לזיהוי על ידי אנדוגני EWS / FLI. בכל תרשים עוגה, אפור מתאר את חלק הגנים שאינם ניצלים על ידי המבנה. אדום מתאר את חלק הגנים המופעלים באופן דיפרנציאלי, וכחול מתאר את חלק הגנים המודחקים באופן דיפרנציאלי. נתון זה מותאם מאת Theisen, et al.23אנא לחץ כאן כדי להציג גרסה גדולה יותר של איור זה.
איור S1: טעינת הקבצים המהירים.gz לסביבת HPC, לגיזום וליישור. אנא לחץ כאן כדי להוריד נתון זה.
איור S2: איסוף ספירת קריאות בין דוגמאות והפעלת נורמליזציה של אצווה עם ComBat. אנא לחץ כאן כדי להוריד נתון זה.
איור S3: הפעלת DESeq2 וחילוץ תוצאות של ניתוח ביטויים דיפרנציאליים. אנא לחץ כאן כדי להוריד נתון זה.
איור S4: ניתוח פלט. אנא לחץ כאן כדי להוריד נתון זה.
איור S5: קיבוץ הירארכי לזיהוי מחלקות גנים: קיבוץ הירארכי של 1000 הגנים המשתנים ביותר בכל המבנים וקו הבסיס, iEF, ממוין לאשכולות k. במקרה זה k=7, אך פרמטר זה מוגדר על-ידי המשתמש כפי שמוצג באיור S4D. אנא לחץ כאן כדי להוריד נתון זה.
טבלה S1: רשימת גנים (מזהה גן אנסמבל) עם ביאור אשכול. נא לחץ כאן כדי להוריד טבלה זו.