$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
איור 1 מספק סקירה של ERRBS, הדגשת שלבים עיקריים, אשר הסבירו בכל הפרוטוקול המתואר. ספריות ERRBS הוכנו באמצעות DNA קלט 50 ng.
להעריך את האיכות של הספריות מוכנות. ייצור הספרייה מניב באופן שיגרתי בגדלים שבריר של 150-250 נקודות בסיס ו250-400 נ"ב (איור 3 א-ג). הבדלים קלים בהפצות גודל ספרייה בין הדגימות צפויים. שים לב כי בשני שברים הספרייה נמוכים וגבוהים יש גדלי DNA מאוד אינטנסיביים, מעידים על ההעשרה של רצף מסוים. תוצאות עיכול MspI בהעשרה של משפחה של DNA חוזר על עצמו רצפי נוכחים בגנום האנושי ב 190 נ"ב, 250 נ"ב ו 310 נ"ב בספריות ERRBS. שלוש חזרות אלה מייצגות חתימה אופיינית לספריית ERRBS 20 (ראה איורים 3 א-C ו3G). ספריות נציג היו רצף ברצף של הדור הבאr חד-סוף השימוש קורא. בעת הטעינה בריכוז הספרייה המומלצת על רצף Illumina HiSeq 2500, צפיפות מקבץ של 500,000-700,000 למ"מ 2 צפויות. בצפיפות זה קיבוץ באשכולות, 81.6% ± 3.14% (n = 81) של האשכולות לעבור סינון (איור 4 א). בשל מורכבות של הסוף הנמוך מוסיף הספרייה (אתר ההכרה MspI: C ^ CGG), ערכי עוצמה וציוני איכות שנרשמו ברצף הם משתנים מאוד בשלושת הבסיסים הראשונים (איור 4 ב-ג), לעומת זאת, אם נתיב שליטה עצמאי כלול (ראה דיון), 85% מבסיסים יהיו ציוני איכות של 30 או יותר (ערכי Q30; איור 4D).
יישור נתונים ונחישות מתילציה ציטוזין כפי שמתואר בנתונים ברזולוציה תשואות פרוטוקול בסיס-זוג (לוח 7). לגנום האנושי, ריצת 51-מחזור קריאה יחידה רצף של ספריית ERRBS בנתיב אחד של HiSeq 2500 במצב תפוקה גבוהה רגילly מייצר 153194882 ± 12918302 כולל קורא כי לאחר סינון איכותי ומתאם זמירה תשואות קוראים 152231183 ± 13189678 לקלט לתוך צינור הניתוח. יעילות המיפוי הממוצעת לספריית ERRBS היא בדרך כלל 62.95% ± 5.92% עם ייצוג של 3183594 ± 713547 CpGs עם כיסוי מינימאלי לCPG של 10x וכיסוי לCPG ממוצע של 84.94 ± 16.29 (n = 100).
פרוטוקול ERRBS ניתן לריבוב (ראה קובץ נוסף 1: הסתגלות פרוטוקול עבור סידור ריבוב). נתונים מרצף נציג הרץ מסוכמים באיור 5 נתונים מריצות רצף ריבוב (לקריאה יחידה ריצת רצף 51-מחזור;. N = 128 לשתי ספריות לכל נתיב; n = 11 לשלוש ספריות לכל נתיב; n = 11 לארבע ספריות לכל נתיב) היו בהשוואה לריצוף נתיב מלא של ספריית ERRBS (ריצות רצף קריאה יחידה 51-מחזור; n = 100), כמו גם downsampling נתיב יחיד לsimulatדואר 50%, 33% ו -25% מקורא לכל נתיב (2, 3, ו -4 ריבוב מדגם לכל נתיב בהתאמה; n = 3). ככל שמספר קורא לדגימה יורד עם גורם הריבוב, מספר CpGs המכוסה בכיסוי מינימאלי של 10x והכיסוי לCPG יורד גם כן (איור 5 ולוח 8). שיעורי המרה ממוצעות של אתרים שאינם CPG הצפויים הם 99.85% ± 0.04% (n = 400). שיעורי המרה נמוכים מ -99% עשויים להצביע על פחות מהמרת bisulfite אופטימלית שיכול לגרום לשיעור גבוה של רמות מתילציה שווא.
נתונים מספריית ERRBS הוכנה מהדנ"א הגנומי אנושי נציג נותחו בR 2.15.2 45 באמצעות חבילת methylKit 26 (ראה קובץ קוד נוסף 1 לפרטים פקודה). ניתן דמיינו נתונים בדפדפנים נפוצים הגנום (איור 6 א). נתונים מתילציה ציטוזין נגזר באופן שווה משני הגדילים (איור 6) ונעים כלספקטרום של רמות פוטנציאל מתילציה ציטוזין (איור 6 ג). ניתוח טכני משכפל מגבוהה התאמה תשואות דגימת DNA האנושיות נציג בין תוצאות נתונים (איור 6 ד) ומכסה CpGs בקשת רחבה של לוקוסים הגנומי (איור 6-ו-ו 'וכפי שתואר לעיל 26). בעוד העתקים טכניים יניבו גבוה R 2 ערכים (יותר מ 97%), העתקים ביולוגיים יניבו R 2 ערכים הנעים 0.92-0.96 26, והשוואה בין סוגים שונים של תאי אדם יניב R 2 ערכים נמוכים יותר מאשר 0.86 (מידע לא מוצג).

איור 1: תרשים זרימה של שלבי פרוטוקול ERRBS. תרשים מייצג את שלבים, שיכול להסתיים ביום עבודה מסורתי. * מציין נקודת הפסקה פוטנציאלית (מייד מעקב קשירת ing לנקות ולפני בחירת גודל, צעד פרוטוקול 5) שביכולות להיות קפוא דגימות ב -20 ° C לפני שתמשכנה עם את משך הזמן של פרוטוקול.

איור 2: פרוטוקול בחירת גודל. (א) צילום מסך של הגדרות המשמשות בפרוטוקול ERRBS פיפין Prep (ראה סעיף 5.1.2 פרוטוקול - 5.1.6): (1) סוג הקלטת בחר. (2) הסטנדרטי בחר לשמש. (3) בחר את מצב הגבייה לכל מסלול. (4) הזן את הטווחים נ"ב אוסף. (5) שמור את פרוטוקול שלבים (ב) לחילוץ ג'ל המדריך לשימוש בסעיף 5.2 פרוטוקול:. (1) סולמות ג'ל דמיינו. (2) מסומן גדלים לבחירת גודל באמצעות סכין גילוח. (3) תמונה של דגימות נכרתה (חלק תחתון: 150-250 נקודות בסיס וחלק גבוה יותר: 250-400 נ"ב)."> לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.

איור 3: תוצאות בקרת איכות לספריות ERRBS נציג הוכנו מדגימות DNA אנושיות באמצעות מכונה bioanalyzer. (א) תמונה דמוית ג'ל מראה סולם סטנדרטי (1), חלק תחתון ספרייה (135-240 נ"ב חלק מפיפין Prep); 2) ושבר ספרייה גבוה יותר (240-410 נ"ב חלק מפיפין Prep); . 3) (electropherogram B) Bioanalyzer של החלק תחתון הספרייה הצפויה (C) electropherogram Bioanalyzer של שבריר הספרייה גבוה יותר הצפוי D -.. F) נציגי נתונים מהכנת ספרייה באיכות ירודה. תמונה כמו ג'ל (D) של הסולם הסטנדרטי (1), חלק תחתון ספרייה (2) ושבר הספרייה גבוה יותר (3). הלהקה ב -150 נ"ב מסומנת בarrow מציין כמויות עודפות של מתאם. Electropherogram של נמוך (E) ושברים גבוהים ספרייה (F) עם פסגות המתאם העודפות ב 150 נ"ב (מסומן בחיצים). Electropherogram (G) Bioanalyzer של ספריית ERRBS נקוותה לסידור. עקבות אדומות מייצגת ספרייה ונקוותה באיכות גבוהה עם ייצוג שווה של שברים גבוהים יותר ונמוכים יותר. עקבות כחולות מייצגת ספרייה במאגר לא מתאים לסידור בשל היעדר הייצוג שווה של שברים גבוהים יותר ונמוכים יותר. לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.

איור 4: רצף תרשימים לנציג ERRBS ריצת רצף קריאה יחיד 51-מחזור ברצף HiSeq 2500 בתפוקה גבוההמצב. (א) צפיפות אשכול (K / מ"מ 2 = 1,000 אשכולות לכל מילימטר מרובעים; כחולה). וצפיפויות אשכול עוברות סינון (ירוק) בשני מסלולים עם ספריות ERRBS (B) עוצמות אופייניות ראו ב -30 המחזורים הראשונים במסלול עם ספריית ERRBS. שים לב לחתימת CGG מעיכול MspI בעוצמות של שלושה המחזורים הראשונים. אחוז (C) של בסיסים עם ציון איכות של 30 ומעלה (%> Q30) לכל מחזור בנתיב ERRBS אחד. חלוקת ציון איכות (ד ') ל כל המחזורים בנתיב ERRBS אחד. כחול = פחות מ Q30, ירוק = גדול או שווה לQ30. במסלול זה, 84.7% מבסיסים היו ציוני איכות של 30 ומעלה.

איור 5: רצף תוצאות פלט. חלקות תיבה של נתוני ניסוי ממדגם מרובב ויחיד לכל רצף נתיב ru ns (מוצג כקופסות ירוקות) ושל נתונים שהופקו על ידי downsampling המדומה מריצות רצף של שלוש ERRBS ספריות (מוצג כקופסות כחולות; נדגמו חמש פעמים לכל ריצת רצף) מ -51 מחזור רצף קריאה יחיד פועל גורם הריבוב מתאים. מספר ספריות ERRBS רצף לכל נתיב. 1 = כל נתיב או 100% מקורא ומייצג נתונים מספריית ERRBS אחת לכל נתיב; 2 = 50% מנתיב ומייצגים נתונים משני ERRBS ספריות לכל נתיב; 3 = 33% מנתיב ומייצגים נתונים משלוש ERRBS ספריות לכל נתיב; ו, 4 = 25% מנתיב ומייצגים נתונים מארבע ERRBS ספריות לכל נתיב. (א) סעיפי הקריאה, או מספר הרצפים ניתחו, לגורם ריבוב. (ב) מספר CPG של מכוסה בנתוני הרצף לריבוב גורם. (ג) הכיסוי הממוצע לCPG לגורם ריבוב._blank "> לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.

איור 6: נציגי נתונים מספריית ERRBS הוכנה מהדנ"א הגנומי אנושי () אוניברסיטת קליפורניה, סנטה קרוז (UCSC) דפדפן הגנום 43 תמונה של נתונים נציג מנתיב רצף ERRBS.. סרגל קנה המידה ציר y מייצג מתילציה% 0-100 בכל ציטוזין מכוסה עם מינימום של 10x. המסלול המותאם אישית העליון מייצג את הגדיל קדימה והמסלול המותאם אישית התחתון מייצג את הגדיל ההפוך. מוצג היא chr12:.. 6,489,523-6,802,422 (hg19) כולל של גני refseq ואיי CPG בתוך אזור גנומי זה היסטוגרמות הפצה (ב) לכיסוי CPG לאורך קדימה לאחור גדילים במדגם CD34 אדם + מח עצם נציג היסטוגרמה (C) הפצהרמות מתילציה CPG לאורך שני הגדילים במדגם CD34 אדם + מח עצם נציג. (ד) עלילת מתאם של רמות מתילציה CPG מהעתק טכני נציג של דגימת DNA אנושית. תרשים (E) פאי הממחיש את הפרופורציות של CpGs המכוסה בERRBS ש מבואר לאיי CPG (אור ירוק), CPG חופים (אפור) ואזורים אחרים (לבן) במדגם מייצג שהוכן מהדנ"א הגנומי אנושי. תרשים (F) פאי הממחיש את הפרופורציות של CpGs המכוסה בERRBS שמבואר ליזמי גן (אדום ), אקסונים (ירוק), אינטרונים (כחול) ואזורים חוץ-(סגולים). אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.
| מֵגִיב | נפח | תגובה |
| 10x אנזים T4 DNA | 10 μl | |
| אדנוזין Deoxynucleotide (dNTP) הפתרון Mix | 4 μl | לערבב של 10 מ"מ של כל נוקלאוטיד |
| T4 DNA פולימראז | 5 μl | 3,000 יחידות / מיליליטר |
| קטע DNA פולימראז אני גדול (Klenow) | 1 μl | 5,000 יחידות / מיליליטר |
| T4 polynucleotide קינאז | 5 μl | 10,000 יחידות / מיליליטר |
| ללא מים DNase | 45 μl | |
טבלת 1:. ריאגנטים סוף תגובת תיקון שמות מגיב וכמויות המשמשים בתגובת תיקון הסוף (שלב פרוטוקול 2.1).
| מֵגִיב | נפח | תגובה |
| מאגר תגובת 10x | 5 μl | למשל, NEBuffer 2 |
| 1 מ"מ 2'-deoxyadenosine 5'-טריפוספט (dATP) | 10 μl | |
| Klenow שבר (3 '→ 5' exo-) | 3 μl | 5,000 יחידות / מיליליטר |
טבלת 2:. ריאגנטים-עוקב תגובת שמות מגיב וכמויות המשמשים בתגובה עוקב-A (צעד פרוטוקול 3.1).
| מֵגִיב | נפח | תגובה |
| 15 מיקרומטר מתאמים מרותקים במים ללא DNase | 3 μl | מתאם 1.0 וPE מתאם PE 2.0; ראה טבלה 4 לרצפים והתייחסות |
| מאגר תגובת 10x אנזים T4 DNA | 581; l | |
| אנזים T4 DNA | 1 μl | 2,000,000 יחידות / מיליליטר |
| ללא מים DNase | 31 μl | |
לוח 3: ריאגנטים מתאם קשירת תגובת שמות מגיב וכמויות המשמשות בתגובת קשירת מתאם (צעד פרוטוקול 4.2)..

לוח 4:. Oligos שימוש בפרוטוקול ERRBS רשימה של oligos בשימוש בכל פרוטוקול ERRBS בתגובת קשירה (צעד פרוטוקול 4) וצעדי PCR הגברה (צעד פרוטוקול 7).
| מֵגִיב | נפח | תגובה |
| wi מאגר תגובת 10x FastStart היי פידליטיה -18 מגנזיום כלוריד מ"מ | 20 μl | |
| 10 מ"מ dNTP פתרון Mix | 5 μl | |
| 25 מיקרומטר פריימר PCR PE 1.0 | 4 μl | ראה טבלת 4 |
| 25 מיקרומטר פריימר PCR PE 2.0 | 4 μl | ראה טבלת 4 |
| FastStart נאמנות גבוהה אנזים | 2 μl | 5 יחידות / μl FastStart תקי DNA פולימראז |
| ללא מים DNase | 125 μl | |
לוח 5: ריאגנטים PCR תגובת שמות מגיב וכמויות המשמשות בתגובת PCR ההגברה (צעד פרוטוקול 7.1)..
| צעד פרוטוקול | מגיב / פירוט פרוטוקול | כמות ה- DNA קלט |
| 5-10 ng | 25 ng | 50 ng |
| 1 | אנזים MspI | 1 μl | 2 μl | 2 μl |
| MspI לעכל נפח תגובה | 50 | 100 | 100 |
| 4 | מתאמים בתגובת קשירה | 1 μl | 2 μl | 3 μl |
| קשירת נפח תגובה | 20 μl | 25 μl | 50 μl |
| 5 | פרוטוקול בחירת גודל | ג'ל הידני בלבד | פיפין Prep או ג'ל מדריך ל | פיפין Prep או ג'ל מדריך ל |
| 7 | ריכוז פריימר PCR | 25 מיקרומטר | 25 מיקרומטר | 10 מיקרומטר במשך 14 מחזורים; 25 מיקרומטר במשך 18 מחזורים |
| מספר מחזורי PCR | 18 | 18 | 14-18 |
לוח 6:. שינויי פרוטוקול צעד לכמויות חומר הזנה החל 5-50 ng כמה צעדים לאורך כל הפרוטוקול מחייבים שינוי של כמויות מגיב משמשות ליצירת ספריות באיכות גבוהות מכמויות שונות של חומרי מוצא. שינויים בכמויות מגיב מפתח כלולים כאן. התאם חיץ ונפחי מים בתגובות בהתאם.
| Chr | בסיס | Strand | כיסוי | freqC | freqT |
| chr1 | 10564 | R | 366 | 85.52 | 14.48 |
| chr1 | 10,571 | F | 423 | 91.25 | 8.75 |
| chr1 | 10542 | F | 432 | 91.2 | 8.8 |
| chr1 | 10563 | F | 429 | 94.64 | 5.36 |
| chr1 | 10572 | R | 366 | 96.99 | 3.01 |
| chr1 | 10590 | R | 370 | 88.11 | 11.89 |
| chr1 | 10,526 | R | 350 | 92 | 8 |
| chr1 | 10543 | R | 368 | 92.93 | 7.07 |
| chr1 | 10525 | F | 433 | 91.92 | 8.08 |
| chr1 | 10497 | F | 435 | 88.74 | 11.26 |
לוח 7: נתונים ERRBS נציג. לאחר יישור נתונים ונחישות מתילציה ציטוזין, נתונים זוג הבסיס מתקבל למכוסה כל CPG, פרוטוקול היישור כפי שתואר יקבע לתאם (עמודות: chr = כרומוזום, בסיס וStrand) הגנומי., שיעור הכיסוי של המוקד הספציפי (סיקור ), ושיעור ציטוזין זיהוי לעומת thymidine כאחוז (freqC וfreqT בהתאמה).
| מספר ספריות ERRBS לכל נתיב | אומר מספר מיושר באופן ייחודי קורא | אומר מספר CpGs המכוסה | אומר כיסוי לCPG |
| 1 | 152231184 ± 13189678 | 3183594 ± 713547 | 85 ± 16 |
| 2 | 77680837 ± 7657058 | 2674823± 153494 | 49 ± 9 |
| 3 | 49938156 ± 2436865 | 2552186 ± - 76624 | 39 ± 2 |
| 4 | 34457208 ± 4441686 | 1814461 ± 144339 | 28 ± 4 |
לוח 8:. פרמטרים נציג מרצף ספריות ERRBS אחת ומרובבות מוצג נתונים לכל נתיב מריצות רצף קריאה אחת 51 מחזור: ממוצע וסטיית התקן של מיושר באופן ייחודי קורא, מספר CpGs מכוסה וכיסוי לכל אתר CPG מתקבל מרצף אחד ספריות ERRBS לכל נתיב (n = 100), שתי ERRBS ספריות לכל נתיב (n = 128), שלוש ERRBS ספריות לכל נתיב (n = 11), וארבע ERRBS ספריות לכל נתיב (n = 11).