$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
הטמעה מוצלחת של פלטפורמת מתיל-Seq עכברוש תלוי מספר קריטריונים. איור 1 מציג את זרימת העבודה הכוללת של המחקר וסימונים צעדים ספציפיים בקרת איכות (QC) הנדרשים לפני שתתקדמו. אחד הגורמים הראשון שיש לקחת בחשבון הוא החוסן של המודל בעלי חיים, משטר מתח, הקובעות את סדר הגודל של שינויים epigenetic המתרחשות על-פני methylome. מאז עבודתנו בעלי חיים מותנה בהשלמת שלנו התצפית הקודם כי חשיפה corticosterone (קורט) יכול להוביל שינויים ב- DNA מתילציה19,20, משטר מתח כרוני משתנה (CVS) שלנו צריכה להיות של התאבנות מספיק כדי לייצר הדגיש חולדות עם רמות מוגברות פלזמה קורט. משטר CVS שבועי טיפוסי שמוצג בטבלה 1 , כללה לחצים יומית בבוקר, אחר הצהריים, בן לילה את זה כל הזמן משתנות כדי למנוע habituation ונחלשו תגובת המתח. לאורך כל הטיפול 3 שבועות, החיות לחוץ הציג רמות גבוהות באופן משמעותי של רשע פלזמה קורט [ימים 4 – 21, לשלוט: 32.7 3.7 ננוגרם למ"ל, מתח: 103.0 11.9 ng/mL (כלומר SEM), P = 2.2 x 10-4, איור 2 א] על אלו של unstressed, לשלוט בבעלי חיים. באופן עקבי, החיות האלה גם הראו התנהגות כמו חרדה גדולה יותר גבוהות פלוס מבוך (EPM), כפי שמציין את הזמן באופן משמעותי יותר בזרועות סגור של EPM ופחות זמן בזרועות פתוחות (איור 2B). תוצאות אלו מדגימים כי החשיפה CVS הוביל אנדוקרינית משמעותיים ושינויי התנהגות, מוביל אותנו לחקור אם שינויים אלה היו קשורים עם חתימות מתילציה DNA ספציפיים.
אנחנו מדגישים מספר מחסומים המהווים מכרעת לבנייה מוצלחת של הספרייה מתיל-Seq. החל בכמות מספקת של ה-DNA הוא הכרחי, sonication, מרובים כביסה/טיהור, היעד העשרה, ביסולפאט המרה צעדים במרוכז להפחית את כמות ה-DNA בספריה סיים. למרות מספר שלבים הגברה PCR להקל על האובדן של תבנית ה-DNA, מספרי מחזור PCR מופרז יכול להציג את קריאות כפולים גבוה יותר. לצורך המחקר הנוכחי עכברוש מתיל-Seq, שימש 2 גר' gDNA דם לכל חולדה. נציין כי ספריות מתיל-Seq יכולה להתבצע באמצעות המתחילים כמות הדנ א 500 ng. מתחיל קטן יותר חומר מאפשר למשתמשים ליצור ספריות מ- DNA מבודדת על ידי FACS (תא לפעיל על-ידי קרינה פלואורסצנטית מיון) או מחט אגרופים, למרות שיש סיכון מוגבר לייצר כמות לא מספקת של ספריות עבור רצף עוקבות. QC מתבצע על ידי אלקטרופורזה של 1 ליטר של המדגם על bioanalyzer, אשר מספק DNA משקל מולקולרי, כמות, ו molarity. שלושה שלבים קריטיים הדורשים את השימוש bioanalyzer הם: 1) בעקבות sonication צעד כדי להבטיח מספיק הטיה של DNA (~ 170 bp, אדום, איור 3); 2) בעקבות צעד מצדו מתאם שצוין על-ידי שינוי בגודל ממוצעת של ה-DNA הוטו (~ 200 bp, כחול, איור 3) כדי להבטיח שלהם הגברה עוקבות מאת PCR; ו 3) בעקבות ספריה הסופי שלב טיהור כדי להבטיח את הכמות והגודל של הספרייה רצף.
R-החבילות של BSSeq BSmooth ב Bioconductor שימשו לניתוח ביסולפאט את רצף נתונים18. החדרים כוללים כלים ושיטות ליישור קוראת את רצף, ביצוע בקרת איכות, זיהוי באופן שונה מפוגל אזורים (DMRs). תוכנה BSmooth הפעלת עניבת הפרפר 2.016,17 כמו רצף פנימי aligner להשיג CpG ברמת מדידה סיכומים על ידי יישור של קריאות קלט גולמי כדי המרה ביסולפאט רצף גנומית. קריאות מיושר ואז מסוננים דרך הליכי בקרת איכות קפדנית המבקשים לזהות רצף שיטתית ושגיאות בסיס-טלפון עשויים להטות ניתוחים במורד הזרם. סדרה של חלקות נוצרות באופן חזותי לסייע בתהליך זה של סינון. מדדי רצף נוצרים גם למידע רלבנטי המסמך, כגון מספר הקריאות מיושר, היעד %, ולפי CpG כיסוי, בין היתר (טבלה 2). לאחר הנתונים מסוננים, באלגוריתם החלקה/נורמליזציה מבוצעת, איפה כל CpG מוקצה ערך מתילציה המשוערת בהתבסס על QC כל קורא מתוך כל דגימה, הערכות מגבולות סחורות ארוזות לצרכן כדי להבטיח קוראים מדויקת יותר של מתילציה מצב אפילו במקרים בהם רצף כיסוי נמוך. ערך זה מספק הערכה מוחלקים של ההסתברות של מתילציה בכל אתר CpG. על ידי השוואת הממוצע של ההערכות מוחלקים מתילציה של כל מדגם בין שתי קבוצות טיפול לבין אזורים גנומית הדירוג מ המעט שונה באופן משמעותי ביותר, נוצרת רשימה של DMRs (טבלה 3).
העליון ש-DMR בין קבוצות לחוץ, unstressed היה ממוקם האמרגן של הגן histocompatibility הגדולות עכברוש Rt1-m4, עם הדגיש חיות מפגין רמות גבוהות יותר של מתילציה מעבר סחורות ארוזות לצרכן כל מבהמות unstressed (איור 4A). כדי לוודא את יישום מוצלח של פלטפורמת מתיל-Seq וניתוח נתונים, תוכננו primers נגד DMR, רמות מתילציה של הדנ א בדם במדגם כולו מתוח, unstressed בעלי חיים (8 וסודרו על ידי מתיל-Seq ו- 8 לא רציף) היו לאומדן ביסולפאט pyrosequencing. תוצאות להדגים עלייה משמעותית מתילציה DNA על פני 10 מתוך 12 סחורות ארוזות לצרכן לבדיקה (שינוי 5.1-10.4% מתילציה, P < 0.037, איור 4B). ניתוח מסלול קג בוצעה על כל אחד DMRs נומינלית משמעותית לזיהוי מסלולים הקשורים עם מתח. באופן עקבי, מסלולים הקשורים DMR מעורבים המחלות הקשורות חשיפה מתח כרוני, כגון סוכרת, מחלות לב וכלי דם, סרטן (טבלה 4). 21 , 22 , 23 להפגין שיוך בין הנתונים epigenetic ומידת החשיפה ללחץ, רמות מתילציה בעשר-CpG הושוו על הרמות קורט כלומר 3 שבועות עבור כל חיה. תוצאות המחקר הראו מתאם צנוע בין הנתונים אנדוקרינית, מתילציה (R2= 0.54, P = 0.001, איור 5).

איור 1: בסך הכל סכמטית זרימת עבודה עבור פלטפורמת עכברוש מתיל-Seq. גרם אחד של ה-DNA גנומי המופק מדם הדגיש, שליטה חולדות מעובד תחילה עבור בניית הספריות של מתיל-Seq עבור רצף, ניתוח, וזיהוי המטרה. עוד 100 ננוגרם של ה-DNA עבור אימות עצמאית של המטרות epigenetic שזוהה משמש ביסולפאט pyrosequencing. אנא לחץ כאן כדי להציג גירסה גדולה יותר של הדמות הזאת.

איור 2: חשיפה מתח כרוני משתנה (CVS) מובילה לשינויים האנדוקרינית והתנהגותיים אצל חולדות. (א) samplings מרובים של corticosterone (קורט) מדגימים את היציבות של השבוע 3 משטר CVS. נאספו דגימות דם בבוקר לפני הטיפול סטרס יומי. (B) חיות Stressed מבלה יותר זמן בזרועות סגור ופחות זמן בזרועות פתוחות גבוהות ועוד מבוך (EPM). Boxplots עם נקודת נתונים עבור כל חיה מוצגים. מבחן T של סטודנט בוצעה מובהקות סטטיסטית. * P < 0.05, * * P < 0.01, ו * * * P < 0.001. אנא לחץ כאן כדי להציג גירסה גדולה יותר של הדמות הזאת.

איור 3: כימות של הוטו, ובין אם-לא מתאם עכברוש דנ א ב bioanalyzer. העקומות אדום וכחול להראות את הכמות והגודל של הדנ א הגנומי (אדום) בעקבות הטיה sonicator איזותרמי, מתאם מצדו, בהתאמה. כל שורה מייצגת דוגמא אחת ואת האדום, כחול עקומות משקפים הן אובדן של DNA במהלך השלבים מספר (end-תיקון, 3'-adenylation וניקוי לדוגמה) להגביר את לחץ הדם גודל בשל מצדו המתאמים. חדות פסגות-25 bp ו- 1500 bp הם סימנים סטנדרטיים נוספו למאגר טעינה. אנא לחץ כאן כדי להציג גירסה גדולה יותר של הדמות הזאת.

איור 4: CVS-induced שינויים epigenetic מזוהים על ידי עכברים מתיל-תת סעיף. (א) ניתוח של העכברוש מתיל-Seq נתונים מעורבים האמרגן של הגן Rt1m4 כאזור מפוגל באופן שונה (DMR) בין לחוץ (אדום) וחולדות שליטה (כחול). פלט גרפי עבור Rt1m4 DMR (אזור מוצל ורוד) מציגה כל CpG (קו אנכי אפור), הדגימות ארבעה לכל קבוצה (קווים אדום או כחול), ואת רמת מתילציה % עבור כל חיה (נקודה אדום או כחול). (B) 12 סחורות ארוזות לצרכן בתוך DMR היו אומת על-ידי ביסולפאט pyrosequencing. מהגרפים מיוצגים אומר SEM, בוצע מבחן T של סטודנט מובהקות סטטיסטית. * P < 0.05. אנא לחץ כאן כדי להציג גירסה גדולה יותר של הדמות הזאת.

איור 5: ניתוח של רגרסיה ליניארית הראו מתאם צנוע בין % DNA מתילציה-CpG-10 מתוך Rt1m4 ובשבוע 3 מתכוון פלזמה רמות קורט של שניהם הדגישו לשלוט בבעלי חיים (N = 16). נתונים מבעלי חיים לחוץ מיוצגים על-ידי עיגולים אדומים. אנא לחץ כאן כדי להציג גירסה גדולה יותר של הדמות הזאת.
| השבוע | יום 1 | יום 2 | יום 3 | יום 4 | יום 5 | יום 6 | יום 7 |
| AM | איפוק | לשחות | חדר קר | לשחות | איפוק | שאכר | לשחות |
| PM | שאכר | כלוב הטיה | איפוק | שאכר | חדר קר | איפוק | חדר קר |
| בן לילה | להגביל את האוכל | רטוב מצעים | בידוד | אור על | הצפיפות | אור על | רטוב מצעים |
טבלה 1: זמנים שבועי טיפוסי של משטר מתח משתנה כרוני (CVS).
| מדדי רצף | לחץ1 | בקרה1 |
| (n = 4) | (n = 4) |
| לזווג סוף קריאות (לפי) | 89,290,397 | 80,165,674 |
| סוף לזווג ממופה באופן ייחודי קורא (UMPER) | 39,200,255 | 35,013,406 |
| יישור קצב/מיפוי יעילות (UMPER /) | 44% | 44% |
| קריאות כפולים (% UMPER) | 73% | 65% |
| UMPER עם מניעת כפילויות | 10,481,031 | 12,306,018 |
| ממוצע לקריאה עומק כיסוי (x) (ARDC) | 6 x | 6 x |
| סחורות ארוזות לצרכן (N) | 12,056,878 | 12,056,878 |
| המרכז לקידום פליטים אפריקאים (x) של סחורות ארוזות לצרכן | 2 x | 2 x |
| סחורות ארוזות לצרכן לפחות 10 קריאות (N) | 481,383 | 595,850 |
| המרכז לקידום פליטים אפריקאים (X) של סחורות ארוזות לצרכן עם קריאות לפחות 10 | 19 | 19 |
| על היעד סחורות ארוזות לצרכן (חפיפה מלאה עם בדיקה אזורי היעד) | 1,923,872 | 2,007,638 |
| המטרה ARDC (x) של סחורות ארוזות לצרכן | 7 x | 8 x |
| על סחורות ארוזות לצרכן היעד לפחות 10 קריאות (N) | 428,249 | 531,419 |
| ב- ARDC היעד (x) של סחורות ארוזות לצרכן לפחות 10 קורא | 18 x | 18 x |
| המטרה (לכל עם חפיפה זוג בסיסים 1 או יותר עם בדיקה אזורי היעד) (UMPER) | 8,277,715 | 9,369,523 |
| % על היעד (של UMPER עם מניעת כפילויות) | 78% | 77% |
| המטרה (סה כ בסיסים ממופה) Mb | 125 mb | 128 מגה-בתים |
| על היעד קריאה ממוצע עומק כיסוי (x) (ARDC) | 9 x | 10 x |
|
1 מדדי רצף המבוסס על ממוצעים על פני נושאים בכל קבוצה | | |
טבלה 2: קביעת רצף מדדים המתקבל פלטפורמת מתיל-Seq חולדה.
| chr | התחל | סוף | ג'ין | מרחק | areaStat | meanDiff | מתח | שליטה | כיוון |
| chr20 | 1,644,246 | 1,644,390 | RT1-M4 | in_gene | 93.03 | 0.22 | 0.33 | 0.11 | רווח |
| chr5 | 160,361,352 | 160,361,564 | LOC690911 | in_gene | -70.75 | -0.19 | 0.72 | 0.91 | אובדן |
| chr3 | 61,138,281 | 61,138,330 | RGD1564319 | 265569 | 61.79 | 0.21 | 0.94 | 0.72 | רווח |
| chr2 | 143,064,811 | 143,065,010 | Ufm1 | 8569 | -59.48 | -0.11 | 0.13 | 0.24 | אובדן |
| chr7 | 30,764,111 | 30,764,284 | Ntn4 | in_gene | 57.04 | 0.21 | 0.94 | 0.73 | רווח |
| chr17 | 12,469,112 | 12,469,218 | Idnk | 41996 | -50.91 | -0.13 | 0.74 | 0.88 | אובדן |
| chr7 | 47,101,725 | 47,101,930 | Pawr | in_gene | -50.54 | -0.12 | 0.64 | 0.76 | אובדן |
| chr5 | 76,111,248 | 76,111,822 | Txndc8 | 151703 | -50.38 | -0.11 | 0.85 | 0.96 | אובדן |
| chr11 | 80,640,132 | 80,640,356 | Dgkg | in_gene | -50.07 | -0.16 | 0.73 | 0.89 | אובדן |
| chr8 | 71,759,248 | 71,759,411 | Mir190 | 210226 | -47.84 | -0.17 | 0.58 | 0.75 | אובדן |
טבלה 3: ראש 10 באופן שונה מפוגל אזורים. עבור כל DMR, טבלת הפלט מראה משמאל לעמודה נכון: מיקום כרומוזומלית (chr), מרכזת (התחלה/סיום), ג'ין שם, מרחק שעתוק להתחיל באתר, אזור דיפרנציאלית סטטיסטיקה הדגישה בין ולשלוט קבוצות (areaStat), כלומר מתילציה דיפרנציאלית (meanDiff), כלומר מתילציה רמות על פני כל DMR עבור הדגיש ולשינוי קבוצות הבקרה (מתח/בקרה), וכיוון של מתילציה פקדים.
| תנאי המעבר קג | ספירת ג'ין | % | ערך-P | בנימיני |
| סוכרת |
| סוכרת סוג II | 12 | 0.1 | 3.6 x 10-4
| 9.8 x 10-3
|
| מחלות לב וכלי דם |
| התכווצות השריר החלק בכלי הדם | 18 | 0.1 | 1.6 x 10-3
| 3.6 x 10-2
|
| Arrhythmogenic חדרית שריר הלב (ARVC) | 13 | 0.1 | 4.0 x 10-3
| 7.1 x 10-2
|
| קרדיומיופתיה מורחבת | 14 | 0.1 | 7.6 x 10-3
| 1.2 x 10-1
|
| נוירון פונקציה |
| הויסעגט | 11 | 0.1 | 1.5 x 10-2
| 1.4 x 10-1
|
| איתות |
| מסלול איתות MAPK | 35 | 0.2 | 2.4 x 10-4
| 9.9 x 10-3
|
| מסלול איתות סידן | 22 | 0.1 | 1.2 x 10-2
| 1.4 x 10-1
|
| מסלול איתות כימוקין | 21 | 0.1 | 1.2 x 10-2
| 1.3 x 10-1
|
| סרטן |
| מסלולים ב סרטן | 42 | 0.3 | 4.1 x 10-5
| 3.4 x 10-3
|
| Glioma | 15 | 0.1 | 4.4 x 10-5
| 2.4 x 10-3
|
| סרטן ריאות תא חדרים-קטן | 10 | 0.1 | 7.9 x 10-3
| 1.1 x 10-1
|
| סרטן המעי הגס | 13 | 0.1 | 8.4 x 10-3
| 1.1 x 10-1
|
| לוקמיה מיאלואידית כרונית | 12 | 0.1 | 1.2 x 10-2
| 1.3 x 10-1
|
טבלה 4: ניתוח קג מסלול של DMRs זיהה מהחולדה מתיל-תת סעיף.