$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
1. צור מניות וירוס
הערה: תרשים זרימה של היבט הספסל הרטוב של פרוטוקול זה מתואר באיור 1 הפרטים של ייצור מניית ויראלי והזיהום הבא של תאים בתרבית הרקמה יחולו בדרך כלל סוגים שונים של רטרווירוסים.. בניסויים מסוימים, תא המטרה לא יכול לבטא את קולטן ויראלי אנדוגני (ים), ובמקרים כאלה בניית חלקיקים retroviral Pseudotyped מחסה גליקופרוטאין מעטפת נגיפיים Heterologous, למשל גליקופרוטאין G מנגיף stomatitis שלפוחי (VSV-G), יהיה נדרש עבור זיהום 44,45.
הערה: זהירות יש לנקוט בעת עבודה עם HIV-1. למרות הנחיות ספציפיות תשתנינה ממוסד למוסד, כל העבודה מבוססת הווירוס צריך להתנהל בתוך ארון בטיחות ביולוגית מוגבל ייעודי, מפעיל (בדרך כלל המכונה במנדף בתרבית רקמה). ציוד מגן אישי נכוןהכולל הגנת פן, כיסויי נעליים, שכבת כפפה כפולה, וחליפת סרבל גוף מלא צריכה להיות משוחקת בכל העת. כל פסולת הנוזל המתקבלת מניסויים הקשורות וירוס צריכה להיות מובטלת עם אקונומיקה (10% ריכוז סופי), וכל הפסולת כולל מוצקים צריכה להיות autoclaved לפני הסילוק.
- יום אחד לפני transfection, צלחת 3.3 x 10 6 תאים HEK293T ב 10 מ"ל של מדיום הנשר שונה של Dulbecco (DMEM) בתוספת 10% (v / v) בסרום שור עוברית 1% (v / v) פניצילין / סטרפטומיצין (10,000 U / ml המלאי) בכל אחת מחמש 100 מנות מ"מ.
הערה: השלימו-DMEM המכונה DMEM-FPS מנקודה זו ואילך. - ביום שלאחר מכן, transfect את התאים עם 10 מיקרוגרם של פלסמיד שנשא שיבוטים מולקולריים retroviral באורך מלא או 9 מיקרוגרם של וקטורים יחידים מסביב נמחק מעטפה עם 1 מיקרוגרם של מבנה ביטוי VSV-G באמצעות ריאגנטים transfection או סידן פוספט זמינים מסחרי.
- דגירת גאמות על 37 מעלות צלזיוס חממה תרבית תאים humidified עם 5% CO 2 (להלן מצב זה המכונה "החממה בתרבית רקמה"). לאחר כ 48 שעות, קציר התקשורת התא המכיל וירוס בעזרת פיפטה נפח ולהעביר אותו דרך פילטר 0.45 מיקרומטר ידי זרימת כוח הכבידה.
- לרכז את הנגיף על ידי ultracentrifugation ב 200,000 XG במשך שעה 1 ב 4 ° C.. Resuspend גלולה הנגיף 500 μl DMEM-FPS המכיל 20 U DNase, דגירה במשך שעה 1 ב 37 מעלות צלזיוס.
הערה: שלב DNase עוזר להפחית את ההתאוששות של רצפי פלסמיד לא רצויים על ידי ביטול בניטל של DNA פלסמיד כי נמשך מהליך transfection.
- קביעת ריכוז p24 46 באמצעות ערכת ללכוד אנטיגן p24 HIV-1 לפי הוראות היצרן.
הערה: ריכוז וירוס יכול להיקבע גם על ידי assay פעילות transcriptase הפוכה 47,48. יכולה לחלופין, רמת נגיף פונקציונלילהיקבע על ידי מדידת משרד הפנים. הדבר נעשה ביותר בקלות באמצעות תא מופעל קרינת מיון עם וירוסי המבטאים גני כתב ניאון כגון חלבון פלואורסצנטי ירוק משופר. קביעת משרד הפנים עשויה להיות שימושית במיוחד כאשר עובדים עם תאים ראשוניים שאינו תומכים באותה הרמה של זיהום כמו שורות תאים מותאמות.
2. תאים להדביק עם וירוס
- פלייט 3.0 x 10 5 תאים HEK293T לכל היטב צלחת 6-היטב 2.5 מ"ל DMEM-FPS דגירה לילה חממה בתרבית רקמה.
הערה: מספר אתרי שילוב הייחודיים התאוששו עם פרוטוקול זה עומד ביחס ישר למספר התאים וכמות וירוס פעילים המשמשים בתעשיית הזיהום. - להדביק תאים עם ריכוז p24 ויראלי סופי של 500 ng / ml בנפח סופי של 500 μl טרי DMEM-FPS עבור שעה 2 באינקובטור בתרבית רקמה, ולאחר מכן להוסיף 2 מ"ל DMEM-FPS מחומם מראש ל -37 מעלות צלזיוס לכל היטב להמשיך דגירה.
- בְּ48 שעות לאחר הפגיעה, הסר את התקשורת לשטוף את התאים עם מ"ל 2 פוספט שנאגרו מלוחים (PBS). הוסף 0.5 מיליליטר טריפסין-EDTA מחומם מראש ל 37 מעלות צלזיוס, ולאחר כמה שניות בדוקות חזותי את הבארות עבור שליפת תא.
- הוסף 2 מ"ל מחומם מראש DMEM-FPS ו resuspend את התאים על ידי עדינה למעלה / למטה pipetting עם טפטפת נפח ~ 10 פעמים. מעבירים את הפתרון בקבוק תרבות 75 ס"מ 2 רקמה המכילה 18 מ"ל מראש חימם DMEM-FPS, דגירה התאים בתרבית רקמה חממה.
- לאחר מינימלית חמישה ימים מתחילת הזיהום, לאסוף את התאים על ידי הסרת התקשורת, לשטוף עם 5 מ"ל PBS, להוסיף 2 מ"ל מראש חימם טריפסין-EDTA, ו resuspend עם 5 מ"ל מראש חימם DMEM-FPS ידי pipetting. צנטריפוגה הפתרון עבור 5 דקות בטמפרטורת החדר ב 2500 XG, וזורקים supernatant.
הערה: למרות אינטגרציה תחת מישורי התנאים הללו בסביבות 48 שעות לאחר הדבקה 49,50, בימים 3 הנוספים של התרבות נדרשות sufficiently לדלל את הריכוז של מולקולות DNA משולבות שנובעות רקומבינציה DNA מבוסס תאים או autointegration ויראלי בתיווך. - חלץ הדנ"א הגנומי מן התא גלול באמצעות ערכה זמינה מסחרי (למשל, לראות 51). Elute את ה- DNA מהעמודה חילוף יונים מסופק עם 200 μl של 10 mM Tris-HCl, pH 8.5.
הערה: Aliquot של תאים להתחלק ב 48 שעות לאחר ההדבקה (שלב 2.3) עבור assay infectivity כדי להבטיח זיהום ראוי וירוס לפני NGS.
3. הדנ"א הגנומי שבר ידי Sonication או על ידי הגבלת אנזים Digest
הערה: שברי Sonication הדנ"א הגנומי באופן רצף עצמאי כמעט, ועל כן במצב העדיף פיצול כאשר רצף דגימות עם קצב החלמה צפוי נמוך (למשל, תאי חולה נגועים או זיהומים יזמו בבית יחסית נמוך משרד פנים). יתר על כן, sonication מאפשרת לאדם להבחין כפילויות PCR של Partiרצף אינטגרציה באתר cular מ ואינטגרציות ייחודיות באותו האתר, שהוא קריטי להבחין הרחבת המשובטים של תאים המכיל provirus בחולים נגועים (ראה שלב 11 להלן) 39,52-54.
הערה: ה- DNA צריך להיות בקע מייד במורד זרם מן הזרם LTR לצמצם הגברה של רצפי ויראלי פנימיים במהלך LM-PCR. אנזים הגבלה BglII שנמצאת 43 נ"ב במורד הזרם מן הרצף U5 במעלה או שאינה עולה בקנה אחד עבור קשירת לאחר מכן עם MseI שנוצר DNA מסתיים עובד היטב עם זנים רבים HIV-1 (איור 1 ב). בעת הכנת ה- DNA על ידי sonication, אנזים ההגבלה הפנימי-הביקוע צריך להיות מיושם לאחר קשירת מקשר (ראה תרשים 1C - E ו- שלב 4.3 להלן).
- עבור sonication, לערבב 10 מיקרוגרם של ה- DNA הגנומי במים nuclease ללא לנפח סופי של 120 μl. Sonicate באמצעות פרמטרים של גודל הפסקה ממוצע של 500 נ"ב (שני סבבים של הסעיף הבאמטר: מחזור עבודה: 5%; עוצם: 3; מחזורי פרץ: 200; זמן: 80 שניות).
- לטהר DNA sonicated באמצעות ערכת טיהור PCR. תיקון ה- DNA מסתיים באמצעות ערכת קצה תיקון DNA ולטהר את ה- DNA באמצעות ערכת טיהור PCR. A-זנב את ה- DNA באמצעות exo Klenow - אנזים ולטהר את ה- DNA A-זנב באמצעות ערכת טיהור PCR. עיין 51,52 לפרטים נוספים שימוש בערכה.
- לעיכול endonuclease הגבלה, לקצץ 10 מיקרוגרם של ה- DNA הגנומי לילה בשעה 37 ° C בהיקף של 100 μl עם חיץ שסיפק יצרן וקוקטייל של אנזימים (100 U אחד) שיוצרים המסוכך 5'-ת"א, וכן אנזים בקנה כגון BglII המסלקת במורד זרם מן LTR ויראלי נגד הזרם. לטהר את הדנ"א למחרת באמצעות ערכת טיהור PCR.
הערה: אף אחד אנזימי הגבלה צריך לחתוך בתוך הטרמינל ~ 30 נ"ב של סוף DNA הנגיפי המועצמת על ידי פרוטוקול LM-PCR. פרוטוקול זה במיוחד מגביר את U5סוף DNA HIV-1.
4. לחשל לינקר Oligonucleotides ולקשור ל- DNA הגנומי קטועה
הערה: כן מקשר סימטרי המכיל סככה התואמת את שברי DNA לעיל (ראה טבלת מס '1 עבור הרצפים של oligonucleotides מנוצלים בפרוטוקול זה). והמקשר לשמש עם DNA sonicated חייב לכלול סככה תואמת T-3 ", בעוד המקשר עבור MseI מתעכל DNA חייב לכלול סככת 5'-ת"א תואמת (איור 1). גדיל המקשר הקצר בנוסף חייב לכלול שינוי כימי בלתי להארכה, כגון 3'-אמין, כדי להגביל את תגובות ההגברה הבאות כלפי DNA של עניין.
הערה: בעת הכנת ספריות באתר אינטגרציה שונות מרובות במקביל ו / או כאשר דגימות ייחודיות ריבוב באותו בטווח הרצף, מומלץ להשתמש linkers ייחודי עבור כל דגימה כדי להגביל את פוטנציאל מדגם צולבות contamination במהלך PCR. זה גם מרמז על השימוש פריימרים מקשר ייחודי עבור כל דגימה במהלך חצי מקוננים PCR (כמתואר להלן). גדילי מקשר ייחודיים פריימרים מקשר יכולים להיות מעוצבים על ידי ערבול רצפי oligonucleotide מקשר המפורטות בטבלת 1, תוך שמירה על תוכן GC הכולל% דומה ועמדות סככה ישימות.
- לחשל גדילי מקשר קצר וארוך 35 μl של 10 mM Tris-HCl, pH 8.0-0.1 EDTA מ"מ (ריכוז סופי של 10 מיקרומטר של כל oligonucleotide) על ידי חימום עד 90 מעלות צלזיוס, קירור לאט לטמפרטורת החדר בצעדים של 1 ° C לדקה.
- הכן לפחות ארבע תגובות קשירת במקביל לכל דגימת DNA גנומי, אשר מכילים מקשר 1.5 מיקרומטר ligated, 1 מיקרוגרם DNA מקוטעת, ו -800 האנזים U T4 DNA ב 50 μl. ולקשור לילה בשעה 12 ° C. לטהר למחרת עם ערכת טיהור PCR.
- לקבלת דוגמיות שהכין sonication, לעכל את התגובה קשירת מטוהרים עם 100 U של restricאנזים tion המסלקת במורד זרם מן הזרם LTR (למשל, BglII ל- HIV-1) תחת היצרן מומלץ תנאי לילה. לטהר את הדנ"א בעזרת ערכת טיהור PCR.
5. להגביר ויראלי LTR-Host הדנ"א הגנומי צומת ידי PCR חצי מקוננות
הערה: כדי להבטיח לגיוון ספרייה אופטימלי, לפחות 4-8 PCRs במקביל, תלוי בריכוז ה- DNA של תגובת הקשירה התאוששה, צריכה להיות מוכנה עבור כל דגימה עבור שני סיבובי PCR. ריכוז תבנית ה- DNA יש לכמת ידי לספקטרופוטומטריה. בפרוטוקול זה הסבב הראשון והשני של ה- PCR להעסיק פריימרים LTR ספציפי מקוננות, אבל באותו פריימר ספציפי מקשר משמש הן סיבובים (טבלה 1). פריימר סיבוב השני LTR הספציפי רצפי המתאם לקודד פריימר ספציפי המקשר עבור אשכולות DNA, כמו גם אתרים מחייב פריימר רצף. פריימר LTR הספציפי המקונן גם מקודד רצף מדד 6 NT, WHich יכול להיות מגוון בקרב פריימרים שונים עבור ספריות ריבוב בתוך אותה בטווח הרצף.
- כן PCRs בסיבוב הראשון המכיל את המרכיבים לכל צינור כמפורט בטבלה 2.
הערה: פריימר מקשר ספציפי מטפח 22 NT ההשלמה כדי והמקשר, טמפרטורת ההתכה של 53 מעלות צלזיוס, תוכן GC של 45%, ועל 3 "סוף ממוקם 15-16 נ"ב במעלה הזרם מן 3" טרמיני של הסיבים הארוכים המקשר השונה (טבלה 1). פריימר בסיבוב הראשון 27 NT LTR יש טמפרטורת התכה של 59 מעלות צלזיוס, תוכן GC של 48%, ו -3 שלה "הסוף ממוקם 34 נ"ב זרם ממסוף U5 HIV-1. האזור של פריימר סיבוב שני 26 NT LTR כי הוא משלים LTR HIV-1 יש טמפרטורת ההתכה של 60 מעלות צלזיוס, תוכן GC של 50%, ו -3 שלה "סוף ממוקם 18 נ"ב במעלה הזרם מן U5 ויראלי מָסוּף. מומלץ טמפרטורת ההתכה oligonucleotide ו-תוכן GC צריך לחקות פרמטרים אלה אם משתמשיםפריימרים PCR עיצוב עם רצפי שינו (כולל לשימוש עם רטרווירוסים אחרים) 21. - לרוץ בכל PCR הראשון תחת פרמטרי thermocycler הבאים: מחזור א ': 94 מעלות צלזיוס למשך 2 דקות; 30 מחזורים: 94 מעלות צלזיוס למשך 15 שניות, 55 ° C למשך 30 שניות, 68 מעלות צלזיוס למשך 45 שניות; אחד מחזור: 68 מעלות צלזיוס במשך 10 דקות.
- פינת תגובות ולטהר באמצעות ערכת טיהור PCR. הכן PCRs סיבוב שני המכיל את המרכיבים לכל צינור לפי טבלה 3. הפעל את הסיבוב השני של ה- PCR באמצעות הפרמטרים thermocycler המתואר בשלב 5.2. פינת התגובות ולטהר את ה- DNA באמצעות ערכת טיהור מסחרית PCR בעקבות הוראות היצרן.
הערה: מגוון רצפי מדד מומלץ תואמים NGS אשכולות DNA זמין 71.
6. בצעו QC ו NGS (הושלם בדרך כלל על ידי מתקן רצף)
- (קוו assay # 1) אישור שלב 5.3 ריכוז ה- DNA הספרייה באמצעות fluoroמטר 55. בקיצור, להכין תקני דגימות ניסוי בנפח סופי של 200 מים nuclease ללא μl. צינורות וורטקס למשך 2-3 שניות, דגירה בטמפרטורת החדר למשך 2 דקות, ולאחר מכן לקרוא את דגימות fluorometer.
הערה: דוגמאות צריכות מכילות ריכוז מינימאלי של ה- DNA ספריית 2 ננומטר בהיקף מינימאלי של 15 μl. - (# Assay QC 2) אשר התפלגות גודל DNA שבר באמצעות assay סרט המבוסס על 56.
הערה: התפלגות אידאלית היא שיא DNA רחב יחסית שבמרכזו 500 נ"ב באורך. אם כמות החומר משמעותית היא גדולה מ 1 kb, אז מומלץ לשלב הליך בגודל מבחר לחסל מיני DNA יותר, אשר יעכבו הגברת גשר במהלך אשכולות. לעומת זאת, אם לשיא משמעותי ניכר סביב 100 עד 200 נ"ב, דימר פריימר נוצר במהלך PCR. במקרה זה ההליך צריך להיות מותאם על מנת למזער את ההיווצרות של הדימרים פריימר. - (קוו assay # 3) Confirm התאגדות תקינה של מתאמים לתוך ספריית DNA על ידי PCR כמוני 57.
- בצע NGS בעקבות ספרות היישום של היצרן. לנצל א-ב ספייק של 10% (w / w) ΦX174 DNA, אשר יהיה לייעל את מדדי איכות בזמן אמת על ידי מתן רכב בסיס מאוזן בטווח הרצף.
הערה: ניסויים רצף באתר אינטגרציה חשופים בדרך כלל נ"ב 150 בסוף יחיד (SE150) או-סוף לזווג 150 נ"ב רצף (PE150). PE150 שימושי במיוחד כדי ללכוד את נקודת החיבור מקשר בכל מולקולת הדנ"א (למשל, כאשר בוחנים אתרי אינטגרציה עבור ראיות של שיבוטי התא המארח).
7. השתמש סקריפט פייתון מותאם או PERL כדי לנתח נתונים רצף עבור רצפי LTR המכיל, יבול משם LTR ו לינקר רצפים, ומפת פניית הגנום עם בלאט
- קבצי FASTA סרוקים LTR המכיל רצף קורא, LTR יבול רצפים מקשרים הרחק רצף DNA מארח גנומית, ולייצא רצפים אלה לקובץ FASTA חדש. מפת קצוץ קורא לשני גנום הפניה (למשל hg19 גרסאות הגנום האנושי או GRCh38) ואת הגנום הנגיפי באמצעות בלאט 58, עם אתר אינטגרציה פלט קואורדינטות לייצא לקובץ .txt נפרד, באמצעות ההגדרות הבאות:
stepSize = 6, minIdentity = 97, ו maxIntron = 0 - לנתח את הפלט בלאט .txt קובץ, להסיר autointegrations (עדות כלומר שסוף LTR שלב לתוך באיזור פנימי של הגנום DNA הנגיפי) ומיפוי רצפים אחרים כדי הגנום HIV-1, וליצור פלט נפרד .txt קובץ שבו כל אתרי האינטגרציה הכפולים כבר מרוכזים אחד, להיטים לתאם ייחודיים.
8. יצירת קבצי .bed המכילים מרווח 15-NT סביב ואינטגרציות, להמיר אותם קבצי FASTA, ולבנות רצף לוגו להצגת העדפות מאגר סביב אתרי אינטגרציה
- יצירת קבצי .bed שמפרטים בהפרש של בסיסיםכל אתר אינטגרציה. לפחות 15 בסיסים (5 במעלה או במורד זרם 10) מוצעים עבור דור לוגו רצף. צור קובץ FASTA מקבצי .bed אלה באמצעות פונקצית fastaFromBed מ BEDTools 59 ו בפקודה זו:
-fi / ספרייה fastaFromBed / אל / הפניות / הגנום / -name -s -bed 15_base_pair_file.bed -fo output_file.fasta
הערה: dinucleotide 5'-CA-3 ה"וירלית משתנה מצטרף לארח DNA במהלך אינטגרציה, ואימות לצומת של התחנה סופית LTR ל- DNA של תאים הוא מסנן ראשוני חשוב לזהות אתרי אינטגרציה בתום לב. אנחנו גם לקמפל לוגו רצף מן האוכלוסייה רצף DNA מארח זה כדי לאמת את תוצאות הניסוי. כמו רטרווירוסים להציג העדפות בסיס חתימה סביב האתרים ושילובם 14,15, הלוגוס רצף לשמש כדי לאמת כי באתרי גנומי הממופים התעוררו תוך השתלבות בתיווך לעומת מנגנוני רקומבינציה אחרים כגון DNA שאינו הומולוגיסוף להצטרף 60,61. - השתמש WebLogo 3 (http://weblogo.threeplusone.com/create.cgi) כדי ליצור לוגו רצף מקובצי FASTA. לחץ על 'בחר קובץ' כדי להעלות קובץ FASTA, ולהשתמש את ההגדרות הבאות: פורמט פלט, PDF (וקטור); גודל לוגו, גדול; מספר המיקום הראשון, -5; מגוון לוגו, -5 עד 5; ציר Y מידה, 0.1, ציר Y המרווח טיק, 0.5, ערכת צבע קלאסי (NA).
9. צור מרכזי מאגר זוג .bed קבצים, בדקו בין זיהום לדוגמא, ולמפות את התפלגות אתרי שילוב הייחודי יחסית תכונות הגנום הרלוונטי
- מאז אינטגרצית retroviral מתרחשת בצורה מדורגת על פני גדילי tDNA, להתאים את הקואורדינטות המדויקות של אתרי אינטגרציה על מנת לשקף את נ"ב המרכזי של שכפול אתר היעד למיפוי נכון של יחסי הפצה גנומי לתכונות הגנומי.
- לכן, עבור 5 נ"ב שכפול וירוסים כמו HIV-1, ליצור קובץ .bed עם נ"ב המרכזי לקזז מן integration האתר על ידי שני בסיסים במורד הזרם למיפוי ואינטגרציות לגדיל פלוס, ושני בסיסים upstream למיפוי ואינטגרציות לגדיל מינוס.
- כדי לבדוק אם קיים בין זיהום מדגם, לחשב את מספר אתרי אינטגרציה משותפים בין הספריות השונות על ידי שימוש BEDTools מצטלב פונקציה מצטלב נ"ב מרכזי .bed קבצים עבור שני מדגמים שונים ועל ידי ביצוע פקודה זו:
bedtools מצטלבים -a central_basepair_1.bed -b central_basepair_2.bed -f 1.00 -r -s> overlap1v2.txt - לספור את מספר השורות בתוך קובץ הפלט overlap1v2.txt על מנת לכמת את המספר המדויק של אתרים נפוצים בקרב שתי הספריות באמצעות הפקודה הבאה:
WC -l overlap1v2.txt - הורד את הקובץ RefSeq ביאור .bed עבור גירסת הגנום התייחסות ששימש למיפוי האתר אינטגרציה ממסד הנתונים ביאור הגנום UCSC (למשל http://hgdownload.cse.ucsc.edu/goldenPath/hg38/database) 62.
- לחשב את מספר אתרי אינטגרציה נופלים בתוך גני RefSeq באמצעות BEDTools מצטלב פונקציה כדי לחתוך את קובץ .bed זוג בסיס המרכזי כי נוצר למדגם עם RefSeq .bed קובץ הבא בפקודה זו:
bedtools מצטלבים -a central_basepair_1.bed -b RefSeq_hg38.bed -u> RefSeq_sample1.bed
- לספור את מספר השורות בתוך קובץ הפלט RefSeq_sample1.bed על מנת לכמת את המספר המדויק של אתרים נופלים בגני RefSeq באמצעות הפקודה הבאה:
WC -l RefSeq_sample1.bed - חזור על שלבים 9.3 ו -9.4 עבור אתרי אינטגרציה מיפוי לכל ביאור אחרים בעלי עניין עבורו מרווח .bed הקובץ זמין. הורד את הקובץ .bed ביאור האי CPG העדכני ביותר עבור הגנום התייחסות עניין ממסד הנתונים ביאור הגנום UCSC כמצוין שלב 9.4.
- לחשב את מספר אתרי אינטגרציה נופל בתוך di מסויםעמדה (מאויר בדוגמא זו היא חלון 5 kb) איים CPG באמצעות פונקצית חלון BEDTools ובעקבות בפקודה זו:
bedtools חלון -w 2500 central_basepair_1.bed -b CpG_hg38.bed -u> CpG_sample1.bed
- לספור את מספר השורות בתוך קובץ הפלט CpG_sample1.bed על מנת לכמת את המספר המדויק של אתרים נופלים בתוך 2.5 kb במעלה או במורד הזרם של איי CPG באמצעות הפקודה הבאה:
WC -l CpG_sample1.bed - חזור על שלבי 9.6 ו -9.7 עבור אתרי אינטגרצית מיפוי סמוך TSSs. צור גרסה חלופית של קובץ RefSeq.bed, שבו גנומי קואורדינטות מיפוי ליותר גן אחד הותאם לשקף מתנת גן יחידה רק במיקום זה. הדבר מונע הערכת יתר של צפיפות גן שמסביב אתרי אינטגרציה. חשב את צפיפות הגן באזור 1 Mb סביב כל אתר אינטגרציה באמצעות הפונקציה חלון BEDTools ובעקבות בפקודה זו:
Bedtools חלון -w 500000 central_basepair_1.bed -b RefSeq_hg38_NonRedundant.bed -u> GeneDensity_sample1.bed - חשב את צפיפות הגן הממוצע עבור כל ואינטגרציות של בסיס הנתונים על ידי ביצוע הפקודה:
awk '(sum + = 7 $) END (הדפסה "= ממוצע", סכום / NR)' GeneDensity_sample1.bed
10. סטטיסטית השוואת התפלגויות אינטגרצית אתר בין דגימות באמצעות שני זנב המבחן המדויק של פישר זנב השני ווילקוקסון דרגת סכום בדיקת R
הערה: הבדיקה המדויקת של שימוש פישר להשוואת שיעור אתרי אינטגרציה בתוך גני RefSeq או בתוך חלון של איי CPG או TSSs, אבל להשתמש במבחן סכום דרגה ווילקוקסון עבור השוואת התפלגות צפיפות גן שמסביב באתרי האינטגרציה. תכנית R זמינה בכתובת http://www.r-project.org/.
שני זנב מבחן מדויק של פישר:
- באמצעות המספרים מחושבים בהתאם להוראות צעדים 9.4 ו -9.7, CReate מטריצות עבור כל השוואה במו מופעים ציינו (אינטגרציות בתוך ביאור או בתוך חלון הסובב ביאור) לעומת נותרי אתרים על ידי הפקודה הזאת:
(Annotation_of_interest <- מטריקס (ג (SampleA # ב, SampleA נותרו #, SampleB # ב # SampleB הנותרים), nrow = 2, dimnames = הרשימה (ג ( 'מרכז', 'השאריות'), ג ( 'SampleA', 'SampleB')))) - חשב את ערך P להשוואה ידי בדיקה מדוקדקת של שני זנב פישר עם הפקודה הבאה:
fisher.test (annotation_of_interest, = החלופית" two.sided ') p.value $
שני זנב דרגה ווילקוקסון סכום בדיקה: - יצירת קובץ .txt טאבים שבו כל עמודה מכילה את שם המדגם בתא העליון, ואחריו מתחת מערכיה צפיפות הגן עבור כל אתרי האינטגרציה בספרייה כי (מתוך קובץ .bed שנוצר בשלב 9.9). מופרד באמצעות טאבים זה לייבא קובץ .txt לתוך R באמצעות הפקודה הבאה ו נהvigating אל ספריית הקבצים הנכונה:
- חשב את ערך P להשוואה ידי מבחן סכום בדרגת שני זנב ווילקוקסון עם הפקודה הבאה:
wilcox.test ($ filename SampleA, FILENAME $ SampleB, 'two.sided' = אלטרנטיבית, לזווג = F, מדויק = T) p.value $
הערה: ערכי P ניתן לחשב רק עד גבול מסוים (נמוך מאוד) ב R, שלאחריו אפס יוחזר על ידי התכנית. עבור מדגמים שונים מסיבי שמניבים P = 0 ב R, לאמוד את שווי P כמו <2.2 x 10 -308.
11. בודקים את נתוני רצף גלם עדות השיבוטים של תאים המכילים DNA ויראלי שולב
הערה: פוטנציאל קטן קיים למעלה אינטגרציה אחד בבית NT בדיוק בגנום ההתייחסות. לחלופין, סינגלהאירוע tegration עלול להפוך הנוכחי מיותר בנתונים רצף עקב השימוש של PCR במהלך הכנת הספרייה ו / או על ידי שכפול התא לפני הכנת ה- DNA. ניתוחים אחרונים של הדנ"א הגנומי של נשאי HIV יש להבחין האפשרויות הללו על ידי זיהוי נקודות מצורף מקשר נקודות / גזירה sonication ייחודי (אשר יכול רק להתעורר לפני PCR) בתוך רצפי DNA המכיל 52-54 אתרים אינטגרציה זהה. אין כרגע ויכוח באשר לשאלה האם proviruses טפח בתוך תאים מורחבים clonally לתרום מאגר ויראלי הסמוי, ועל כן הוא בעלת עניין מיוחד לאפיין רמת הרחבתן כאשר לומדים אתרי אינטגרציה בחולים אנושיים.
- בדומה ההליך שתואר בשלב 8.1, ליצור קבצי .bed רישום בהפרש של בסיסי הארכה, במקרה זה, 25 NT במורד זרם מכל אתר שילוב ייחודי (בסיסי upstream מיותרים כאן). צור קובץ FASTA מקבצי .bed אלה (בהתאם להוראותשלב 8.1) באמצעות הפונקציה fastaFromBed מ BEDTools ובעקבות בפקודה זו:
-fi / ספרייה fastaFromBed / אל / הפניות / הגנום / -name -s -bed 25_base_pair_file.bed -fo output_file.fasta
הערה: כדי לשפר את הספציפיות של כל חיפוש מומלץ לחלץ לפחות 25 NT במורד הזרם מכל אתר אינטגרציה עבור שיבוטי מנתח. - רצוי בעזרת סקריפט מותאם אישית, לחפש בקובץ הנתונים הגולמיים רצף FASTA לכל מחרוזות המכילות התאמה מדויקת אל 25 NT במורד הזרם מכל אתר שילוב ייחודי, וכן להפקיד רצפים אלה לתוך קובץ חדש. חתוך LTR רצפים מקשר בין המיתרים גלם. מיזוג רצף PE קורא ידי המרת מקריא המשלים הפוכה, זמירה רצפי LTR ו מקשר, ולאחר מכן הקצאת מחרוזות read2 לזוג read1 שלהם אם מחרוזות לשתף לפחות 20 חופפים NT.
- סרוק את הנקודות מצורפות מקשר של כל בלוק באתר אינטגרציה. לסווג כל אינטגרציה כמו "מורחבת clonally &# 34; אם נקודות מצורפות מקשר הן ≥3 BP בנפרד.
הערה: פרוטוקול לניתוח שיבוטים ללא מיזוג רצף קורא תואר 52.
הערה: פיצול של הגנום במיקום הזהה על ידי sonication מוביל הערכה נמוכה מדי של היקף ההרחבה משובטת, ושיטות לתיקון הטיה הניסיונית המתקבלת תואר 63,64.