$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
1. توليد مخزونات فيروس
ملاحظة: وصفت مخطط تدفق الجانب مقاعد البدلاء الرطب من هذا البروتوكول في الشكل 1 وتفاصيل إنتاج الأسهم الفيروسية والعدوى لاحق من خلايا نسيج الثقافة ينطبق عموما على أنواع مختلفة من الفيروسات. بالنسبة لبعض التجارب، والخلية المستهدفة قد لا تعبر عن مستقبلات الفيروسية الذاتية (ق)، وفي مثل هذه الحالات بناء جزيئات فيروس pseudotyped إيواء مغاير الفيروسي بروتين سكري مغلف، مثل بروتين سكري G من فيروس التهاب الفم الحويصلي (VSV-G)، سيكون مطلوب للعدوى 44،45.
ملاحظة: يجب أن تؤخذ الاحتياطات عند العمل مع HIV-1. على الرغم من مبادئ توجيهية محددة تختلف من مؤسسة لأخرى، ينبغي أن تتم جميع الأعمال القائمة على الفيروس في مخصص، مشغل تقتصر خزانة السلامة البيولوجية (عادة يشار إلى نسيج الثقافة هود). معدات الوقاية الشخصية المناسبةيتضمن حماية الوجه، ويغطي الحذاء، وطبقة قفاز مزدوجة، ودعوى المعطف كامل الجسم يجب أن ترتديه في جميع الأوقات. وينبغي العمل جميع النفايات السائلة الناتجة من التجارب ذات الصلة الفيروس مع التبييض (تركيز النهائي 10٪)، وجميع النفايات بما في ذلك المواد الصلبة يجب تعقيمها قبل التخلص منها.
- قبل يوم واحد ترنسفكأيشن، لوحة 3.3 × 10 6 خلايا HEK293T في 10 مل من Dulbecco لتعديل النسر متوسطة (DMEM) تستكمل مع 10٪ (ت / ت) مصل بقري جنيني و 1٪ (ت / ت) البنسلين / الستربتومايسين (10000 U / مل الأسهم) في كل من خمسة أطباق 100 ملم.
ملاحظة: تستكمل-DMEM ويشار إلى DMEM-FPS من هذه النقطة. - وفي اليوم اللاحق، transfect الخلايا مع 10 ميكروغرام البلازميد تحمل كامل طول استنساخ الجزيئية فيروسات أو 9 ميكروغرام من ناقلات جولة واحدة في حذف المغلف مع 1 ميكروغرام من التعبير بناء VSV-G باستخدام الكواشف ترنسفكأيشن المتاحة تجاريا أو فوسفات الكالسيوم.
- احتضان جملتعلمي اللغة اإلنكليزية عند 37 درجة مئوية في حاضنة الثقافة خلية ترطيب مع 5٪ CO 2 (هذا الشرط الآخرة المشار إليها باسم "نسيج الثقافة الحاضنة"). بعد ما يقرب من 48 ساعة، حصاد وسائل الإعلام الخلية التي تحتوي على الفيروس باستخدام ماصة حجمية وتمر من خلال 0.45 ميكرون تصفية حسب تدفق الجاذبية.
- تركيز الفيروس عن طريق تنبيذ فائق في 200000 x ج لمدة 1 ساعة على 4 درجات مئوية. resuspend الكرية الفيروس في 500 ميكرولتر DMEM-FPS تحتوي على 20 U الدناز، واحتضان لمدة 1 ساعة على 37 درجة مئوية.
ملاحظة: الخطوة الدناز يساعد على الحد من استعادة تسلسل البلازميد غير المرغوب فيها من قبل القضاء وطأة DNA البلازميد الذي لا يزال قائما من الإجراء ترنسفكأيشن.
- تحديد تركيز P24 46 باستخدام HIV-1 P24 مستضد عدة أسر حسب تعليمات الشركة الصانعة.
ملاحظة: يمكن أيضا تركيز الفيروسات يحددها النشاط الناسخ العكسي فحص 47،48. بدلا من ذلك، مستوى الفيروس وظيفية يمكنيتحدد من خلال قياس وزارة الداخلية. ويتم ذلك بسهولة أكبر باستخدام مضان تنشيط الخلايا الفرز مع الفيروسات التي تعبر عن الجينات مراسل الفلورسنت مثل تعزيز البروتين الفلوري الأخضر. قد يكون تقرير وزارة الداخلية مفيدة بشكل خاص عند العمل مع الخلايا الأولية التي قد لا تدعم نفس المستوى من العدوى كما خطوط الخلايا الأمثل.
2. الخلايا تعدى مع الفيروسات
- لوحة 3.0 × 10 5 خلايا HEK293T لكل بئر في لوحة 6 جيدا في 2.5 مل DMEM-FPS واحتضان بين عشية وضحاها في نسيج الثقافة الحاضنة.
ملاحظة: عدد من المواقع التكامل فريدة تعافى مع هذا البروتوكول يتناسب طرديا مع عدد الخلايا وكمية الفيروس النشط المستخدمة في العدوى. - تصيب الخلايا مع تركيز P24 الفيروسي النهائي من 500 نانوغرام / مل في الحجم النهائي من 500 ميكرولتر الطازجة 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 مل برنامج تلفزيوني، إضافة 2 مل قبل تحسنت التربسين-EDTA، و resuspend مع 5 مل قبل تحسنت DMEM-FPS قبل pipetting. أجهزة الطرد المركزي في حل لمدة 5 دقائق في درجة حرارة الغرفة في 2500 x ج، وتجاهل طاف.
ملاحظة: على الرغم من التكامل في ظل هذه الظروف الهضاب في حوالي 48 ساعة بعد العدوى 49،50، يطلب من 3 أيام إضافية من ثقافة إلى ذلك كافياذ تخفيف تركيز جزيئات الحمض النووي غير متكاملة التي تنتج عن خلية القاعدة إعادة التركيب الحمض النووي أو autointegration الفيروسية التي تتوسط فيها. - استخراج الحمض النووي الجيني من بيليه الخلية باستخدام مجموعة المتاحة تجاريا (على سبيل المثال، انظر 51). أزل الحمض النووي من العمود التبادل الأيوني المتوفرة مع 200 ميكرولتر من 10 ملي تريس، حمض الهيدروكلوريك، ودرجة الحموضة 8.5.
ملاحظة: قسامة من الخلايا يجب أن توزع في 48 ساعة بعد العدوى (الخطوة 2.3) لفحص العدوى لضمان عدوى فيروس الصحيح قبل خ ع.
3. جزء الجينوم الحمض النووي عن طريق صوتنة أو تقييد انزيم دايجست
ملاحظة: شظايا صوتنة الحمض النووي الجيني بطريقة مستقلة عن التسلسل تقريبا، وبالتالي هو الوضع الأفضل للتجزئة عندما تسلسل العينات مع معدل الانتعاش المتوقع انخفاض (على سبيل المثال، الخلايا المصابة المريض أو التهابات بدأت في الانخفاض النسبي وزارة الداخلية). وعلاوة على ذلك، صوتنة يسمح احد لتمييز التكرارات PCR للحزبدائرية تسلسل موقع التكامل من التكامل فريدة من نوعها في نفس الموقع، وهو أمر حاسم للتمييز بين التوسع نسيلي من الخلايا التي تحتوي على طليعة الفيروس لدى المرضى المصابين (راجع الخطوة 11 أدناه) 39،52-54.
ملاحظة: يجب أن يكون المشقوق الحمض النووي فورا المصب من LTR المنبع لتقليل التضخيم من متواليات الفيروسية الداخلية خلال LM-PCR. الانزيم تقييد BglII التي تقع 43 نقطة أساس المصب من تسلسل U5 المنبع ويتعارض ذلك لربط لاحق مع الحمض النووي ولدت MseI-ينتهي يعمل بشكل جيد مع العديد من HIV-1 سلالات (الشكل 1B). عند إعداد الحمض النووي عن طريق صوتنة، ينبغي أن تطبق على انزيم التقييد الشق الداخلي بعد ربط رابط (انظر الشكل 1C - E والخطوة 4.3 أدناه).
- لصوتنة، خلط 10 ميكروغرام من الحمض النووي الجيني في nuclease خالية من المياه إلى الحجم النهائي من 120 ميكرولتر. يصوتن باستخدام معلمات لمتوسط حجم كسر 500 سنة مضت (جولتين من الفقرة التاليةمتر: دورة العمل: 5٪. كثافة: 3؛ دورة في انفجار: 200. الوقت: 80 ثانية).
- تنقية الحمض النووي sonicated باستخدام طقم تنقية PCR. إصلاح الحمض النووي ينتهي باستخدام الحمض النووي عدة النهائي إصلاح وتنقية الحمض النووي باستخدام طقم تنقية PCR. A-ذيل الحمض النووي باستخدام Klenow إكسو - انزيم وتنقية الحمض النووي الذيل وباستخدام طقم تنقية PCR. الرجوع إلى 51،52 لمزيد من التفاصيل من استخدام عدة.
- لهضم تقييد نوكلياز داخلية، وقطع 10 ميكروغرام من الحمض النووي الجيني بين عشية وضحاها في 37 درجة مئوية في حجم 100 ميكرولتر مع العازلة التي قدمتها الشركة المصنعة ومزيج من الأنزيمات (100 U لكل منهما) التي تولد يتدلى 5'-TA، فضلا عن انزيم يتعارض مثل BglII أن يشق مجرى النهر من المنبع LTR الفيروسي. تنقية الحمض النووي في اليوم التالي باستخدام طقم تنقية PCR.
ملاحظة: أي من الانزيمات التقييد يجب أن تقطع داخل محطة ~ 30 نقطة أساس في نهاية الحمض النووي الفيروسي الذي يتم تضخيمه من قبل بروتوكول LM-PCR. هذا البروتوكول تضخيم تحديدا U5نهاية HIV-1 الحمض النووي.
4. يصلب رابط [أليغنوكليوتيد وLigate إلى تجزئتها الجينوم الحمض النووي
ملاحظة: إعداد لرابط غير المتماثلة التي تحتوي على المتراكمة والتي تتوافق مع شظايا الحمض النووي أعلاه (انظر الجدول رقم 1 لتسلسل أليغنوكليوتيد] المستخدمة في هذا البروتوكول). رابط ليتم استخدامها مع الحمض النووي sonicated يجب أن يحتوي على متوافق تي 3 "فيض، في حين رابط لالحمض النووي يهضم MseI-يجب أن يحتوي على متوافق 5'-TA تراكم (الشكل 1). يجب أن تحتوي على حبلا رابط قصير بالإضافة إلى تعديل الكيميائية غير قابلة للتمديد، مثل 3'-أمين، لكبح ردود الفعل التضخيم لاحقة نحو الحمض النووي من الفائدة.
ملاحظة: عند إعداد عدة مكتبات موقع التكامل المختلفة بصورة متوازية و / أو عندما عينات فريدة مضاعفة على نفس المدى التسلسل، فمن المستحسن استخدام linkers فريدة لكل عينة للحد من إمكانية عينة عبر كونتامانأوجه خلال PCR. وهذا يعني بالإضافة إلى استخدام بادئات رابط فريدة لكل عينة خلال شبه متداخلة PCR (هو موضح أدناه). يمكن تصميم مسارات رابط فريدة من نوعها والاشعال رابط من قبل الهرولة متواليات النوكليوتيد رابط المدرجة في الجدول 1 مع الحفاظ مماثل العام المحتوى٪ GC والمواقف المتراكمة المعمول بها.
- يصلب فروع رابط قصيرة وطويلة في 35 ميكرولتر من 10 ملي تريس، حمض الهيدروكلوريك، ودرجة الحموضة 8،0-0،1 مم EDTA (تركيز النهائي من 10 ميكرومتر من كل النوكليوتيد) عن طريق التسخين إلى 90 درجة مئوية والتبريد ببطء إلى درجة حرارة الغرفة في الخطوات من 1 ° C لكل دقيقة.
- إعداد أربعة على الأقل ردود الفعل ربط موازية لكل عينة الحمض النووي الجيني، والتي تحتوي على 1.5 ميكرومتر ligated رابط، 1 ميكروغرام الحمض النووي مجزأة، و 800 يغاز U T4 الحمض النووي في 50 ميكرولتر. Ligate بين عشية وضحاها في 12 درجة مئوية. تنقية اليوم التالي مع عدة تنقية PCR.
- للحصول على عينات من صوتنة أعدت، هضم رد فعل ربط تنقيته مع 100 U من restricانزيم نشوئها أن يشق مجرى النهر من المنبع LTR (على سبيل المثال، BglII لHIV-1) تحت الشركة المصنعة أوصت الظروف بين عشية وضحاها. تنقية الحمض النووي باستخدام طقم تنقية PCR.
5. أسهب الفيروسية LTR-المضيف الجينوم المفارق الحمض النووي عن طريق PCR متداخلة نصف
ملاحظة: لضمان التنوع مكتبة الأمثل، 4-8 على الأقل تقارير إنجاز المشروعات موازية، اعتمادا على تركيز الحمض النووي للتفاعل ربط استردادها، يجب أن تكون مستعدة لكل عينة لالجولتين PCR. يجب أن يكون كميا تركيز قالب الحمض النووي عن طريق القياس الطيفي. في هذا البروتوكول في الجولتين الأولى والثانية من PCR تستخدم بادئات LTR محددة متداخلة، ولكن يتم استخدام نفس التمهيدي رابط معين لالجولتين (الجدول 1). والثاني التمهيدي LTR محددة مستديرة وتسلسل محول التمهيدي ترميز رابط معين لتجميع الحمض النووي وكذلك التسلسل مواقع ملزم التمهيدي. أيضا بترميز التمهيدي LTR محددة متداخلة تسلسل المؤشر 6 الإقليم الشمالي، ذوي الخوذات البيضاءالتراث الثقافي غير المادي يمكن أن تختلف بين الاشعال مختلفة للمكتبات المتنوعة داخل نفس المدى التسلسل.
- إعداد تقارير إنجاز المشروعات الجولة الاولى التي تحتوي على المكونات في أنبوب كما هو موضح في الجدول رقم (2).
ملاحظة: التمهيدي رابط معين المرافئ 22 الإقليم الشمالي التكامل على رابط، ودرجة حرارة انصهار من 53 درجة مئوية، وهو GC-المحتوى من 45٪، ولها 3 'نهاية يقع 15-16 سنة مضت المنبع من 3' ترميني من فروع رابط مختلفة طويلة (الجدول 1). الجولة الأولى 27 الإقليم الشمالي LTR التمهيدي لديه درجة حرارة انصهار من 59 درجة مئوية، وهو GC-المحتوى من 48٪، ويقع "نهايته 3 34 سنة مضت المنبع من U5 محطة HIV-1. منطقة الدور الثاني التمهيدي 26 الإقليم الشمالي LTR التي هي مكملة لHIV-1 LTR لديها درجة حرارة انصهار 60 درجة مئوية، وهو GC-المحتوى من 50٪، ونهايته 3 "يقع 18 سنة مضت المنبع من U5 الفيروسي محطة. فمن المستحسن أن درجة الحرارة النوكليوتيد ذوبان وGC-المحتوى يجب أن تحاكي هذه المعايير إذا المستخدمينالاشعال تصميم PCR مع تسلسل المتغيرة (بما في ذلك للاستخدام مع الفيروسات الأخرى) 21. - تشغيل أول جولة PCR تحت المعلمات thermocycler التالية: دورة واحدة: 94 درجة مئوية لمدة 2 دقيقة. 30 دورات: 94 درجة مئوية لمدة 15 ثانية و 55 درجة مئوية لمدة 30 ثانية، 68 درجة مئوية لمدة 45 ثانية؛ دورة واحدة: 68 درجة مئوية لمدة 10 دقيقة.
- تجميع ردود الفعل وتنقية باستخدام طقم تنقية PCR. إعداد تقارير إنجاز المشروعات الدور الثاني يحتوي على المكونات في أنبوب كما في الجدول 3. تشغيل الجولة الثانية من PCR باستخدام المعلمات thermocycler هو موضح في الخطوة 5.2. تجميع ردود الفعل وتنقية الحمض النووي باستخدام PCR عدة تنقية التجارية باتباع إرشادات الشركة المصنعة.
ملاحظة: مجموعة متنوعة من متواليات مؤشر أوصى متوافقة مع خ ع الحمض النووي تجميع تتوفر 71.
6. تنفيذ مراقبة الجودة وخ ع (مكتمل عادة من قبل مرفق التسلسل)
- (QC فحص رقم 1) تأكيد الخطوة 5.3 تركيز الحمض النووي المكتبة باستخدام الفلوريةمتر 55. لفترة وجيزة، وإعداد معايير ونماذج تجريبية في الحجم النهائي من 200 nuclease خالية من المياه ميكرولتر. أنابيب دوامة لمدة 2-3 ثانية، احتضان في درجة حرارة الغرفة لمدة 2 دقيقة، ثم قراءة العينات في التألق.
ملاحظة: يجب أن العينات تحتوي على الحد الأدنى من تركيز 2 نانومتر الحمض النووي مكتبة في حد أدنى لحجم 15 ميكرولتر. - (QC فحص رقم 2) تأكيد توزيع حجم تفتيت الحمض النووي باستخدام الفحص القائمة على الشريط 56.
ملاحظة: إن التوزيع المثالي هو ذروة DNA واسعة نسبيا تتمحور حول 500 شركة بريتيش بتروليوم في طول. إذا كمية كبيرة من المواد أكبر من 1 كيلو بايت، ثم فمن المستحسن أن تتضمن إجراء الحجم اختيار للقضاء على الأنواع الحمض النووي لفترة أطول، والتي سوف تعيق جسر التضخيم خلال المجموعات. على النقيض من ذلك، إذا كان ذروة كبيرة هو واضح حول 100-200 نقطة أساس، وهو ديمر التمهيدي قد شكلت خلال PCR. في هذه الحالة يجب أن يكون الأمثل الإجراء للحد من تشكيل dimers التمهيدي. - (QC فحص رقم 3) شركةnfirm التأسيس السليم للمحولات في مكتبة الحمض النووي عن طريق الكمي PCR 57.
- أداء خ ع التالية الأدب تطبيق الشركة المصنعة. الاستفادة من ارتفاع في ومن 10٪ (ث / ث) ΦX174 الحمض النووي، والتي سوف تحسين مقاييس الجودة في الوقت الحقيقي من خلال توفير تكوين قاعدة متوازنة لتشغيل التسلسل.
ملاحظة: عادة ما يتعرضون تجارب موقع التكامل التسلسل إلى واحد نهاية 150 سنة مضت (SE150) أو إقران نهاية 150 سنة مضت (PE150) التسلسل. PE150 مفيدة بشكل خاص للقبض على نقطة التعلق رابط على كل جزيء الحمض النووي (على سبيل المثال، عند التدقيق في مواقع التكامل عن أدلة على التوسع نسيلي الخلية المضيفة).
7. استخدام بيثون أو برل مخصص لتحليل البيانات تسلسل لاحتواء LTR-متواليات، والمحاصيل بعيدا LTR و رابط متواليات، وتعيين إلى المرجعي الجينوم مع BLAT
- ملفات FASTA المسح الضوئي لLTR التي تحتوي على تسلسل يقرأ، LTR المحاصيل وتسلسل رابط بعيدا عن المضيف تسلسل الحمض النووي الجيني، وتصدير هذه السلاسل إلى ملف FASTA الجديد. خريطة اقتصاص يقرأ على كل من الجينوم المرجعية (مثل الإصدارات الجينوم البشري hg19 أو GRCh38) والجينوم الفيروسي باستخدام بلاط 58، مع موقع التكامل الناتج إحداثيات تصديرها إلى ملف txt منفصل، وذلك باستخدام الإعدادات التالية:
stepSize = 6، minIdentity = 97، وmaxIntron = 0 - تحليل الإخراج BLAT. TXT الملف، إزالة autointegrations (أي أدلة على أن نهاية LTR وإدماجها في المنطقة الداخلية للجينوم الحمض النووي الفيروسي) ورسم تسلسل الأخرى لجينوم HIV-1، وإنشاء الناتج منفصل. TXT الملف الذي تم تكثيف جميع المواقع التكامل مكررة في واحدة، وتنسيق الزيارات فريدة من نوعها.
8. إنشاء ملفات .bed يحتوي على فترات 15-NT التكاملات المحيطة، تحويل هذه الملفات إلى FASTA، وإنشاء تسلسل شعارات لعرض تفضيلات قاعدة محيطة مواقع التكامل
- إنشاء ملفات .bed تلك القائمة فاصل من قواعد لكل موقع التكامل. واقترح لا يقل عن 15 قواعد (5 المنبع والمصب 10) لتوليد تسلسل الشعار. إنشاء ملف FASTA من هذه الملفات .bed باستخدام الدالة fastaFromBed من BEDTools 59 وهذا الأمر:
fastaFromBed -فاي / الدليل / إلى / إشارة / الجينوم / -اسم -s -bed 15_base_pair_file.bed -fo output_file.fasta
ملاحظة: انضم ثابتة الفيروسي 5'-CA-3 "ثنائي النوكليوتيد لاستضافة الحمض النووي خلال التكامل، والتحقق من تقاطع محطة LTR إلى الحمض النووي الخلوي هو مرشح الأولي المهم تحديد بحسن نية مواقع التكامل. نحن بالإضافة إلى تجميع الشعارات تسلسل من هذا المضيف تسلسل الحمض النووي السكان للتحقق من النتائج التجريبية. كما تعرض الفيروسات القهقرية تفضيلات قاعدة توقيع المحيطة المواقع اندماجهم 14،15، والشعارات تسلسل تعمل على التحقق من أن المواقع الجينومية معين نشأت من خلال التكامل بوساطة IN-بالمقارنة مع آليات إعادة التركيب أخرى مثل الحمض النووي غير المتجانسإنهاء الانضمام 60،61. - استخدام WebLogo 3 (http://weblogo.threeplusone.com/create.cgi) لإنشاء الشعارات تسلسل من ملفات FASTA. انقر على "اختيار ملف" لتحميل ملف FASTA، واستخدام الإعدادات التالية: تنسيق الإخراج وقوات الدفاع الشعبي (ناقلات)؛ حجم الشعار، عموما؛ عدد المركز الأول، -5. مجموعة الشعار، -5 إلى 5؛ على نطاق والعمودي، 0.1، Y-محور تباعد عرة، 0.5، نظام الألوان الكلاسيكية (NA).
9. خلق الوسطى قاعدة زوج .bed الملفات، التحقق من وجود نموذج عبر التلوث، وخريطة توزيع مواقع تكامل فريد نسبة إلى ذات الصلة الجينوم الميزات
- منذ يحدث التكامل فيروسات بطريقة متداخلة في الفروع tDNA، وضبط الإحداثيات الدقيقة للمواقع التكامل لتعكس بي بي المركزي للازدواجية لرسم الخرائط الصحيح لتوزيع الجينومية بالنسبة للميزات الجينومية الموقع المستهدف.
- لذلك، لمدة 5 بي بي بتكرار الفيروسات مثل فيروس نقص المناعة البشرية-1، إنشاء ملف .bed مع شركة بريتيش بتروليوم المركزي تعويض من أناموقع ntegration من قبل اثنين من قواعد المصب لرسم الخرائط التكامل إلى زائد حبلا، وقاعدتين المنبع لرسم الخرائط التكامل لحبلا ناقص.
- للتحقق من عينة انتقال التلوث، وحساب عدد من المواقع التكامل المشتركة بين المكتبات المختلفة باستخدام BEDTools تتقاطع وظيفة لتتقاطع بي بي المركزي .bed الملفات لمدة عينات مختلفة وباتباع هذا الأمر:
bedtools تتقاطع -a central_basepair_1.bed -b central_basepair_2.bed -f 1.00 -r -s> overlap1v2.txt - حساب عدد خطوط داخل ملف الإخراج overlap1v2.txt من أجل تحديد العدد الدقيق للمواقع المشتركة بين المكتبتين باستخدام الأمر التالي:
مرحاض -l overlap1v2.txt - تحميل ملف .bed RefSeq الشرح لإصدار الجينوم المرجعية التي تم استخدامها لرسم خرائط الموقع التكامل من قاعدة بيانات 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 باستخدام الأمر التالي:
مرحاض -l RefSeq_sample1.bed - كرر الخطوات من 9.3 و 9.4 لمواقع التكامل رسم الخرائط إلى أي شرح آخر من الفائدة التي فاصل .bed يتوفر الملف. تحميل الدليل السياسي الشامل ملف .bed الشرح أكثر الجزر الحالي للجينوم إشارة من الاهتمام من قاعدة البيانات UCSC الجينوم الشرح كما هو مشروح في الخطوة 9.4.
- حساب عدد من المواقع التكامل التي تقع ضمن دي معينموقف (موضح في هذا المثال هو نافذة 5 كيلو بايت) من الجزر الدليل السياسي الشامل باستخدام وظيفة BEDTools النافذة وبعد هذا الأمر:
bedtools نافذة -w 2500 central_basepair_1.bed -b CpG_hg38.bed -u> CpG_sample1.bed
- حساب عدد خطوط داخل ملف الإخراج CpG_sample1.bed من أجل تحديد العدد الدقيق للمواقع التي تدخل في نطاق 2.5 كيلوبايت المنبع أو المصب من الجزر الدليل السياسي الشامل باستخدام الأمر التالي:
مرحاض -l CpG_sample1.bed - كرر الخطوات من 9.6 و 9.7 لمواقع التكامل رسم الخرائط TSSs قريب. توليد صيغة بديلة للملف RefSeq.bed، حيث ينسق الجينومية رسم الخرائط لأكثر من جين واحد قد تم تعديلها لتعكس فقط هدية جين واحد في هذا الموقف. هذا يمنع المبالغة في تقدير كثافة الجينات المحيطة مواقع التكامل. حساب كثافة الجينات في المنطقة ميجا بايت 1 المحيطة بكل موقع التكامل باستخدام الدالة BEDTools النافذة وبعد هذا الأمر:
bedtools نافذة -w 500000 central_basepair_1.bed -b RefSeq_hg38_NonRedundant.bed -u> GeneDensity_sample1.bed - حساب متوسط كثافة الجينات لجميع التكامل في مجموعة البيانات باتباع هذا الأمر:
AWK '(المبلغ + = 7 $) END (المطبوعة "متوسط =" مبلغ / NR) "GeneDensity_sample1.bed
10. قارن إحصائيا التكامل الموقع التوزيعات بين العينات عن طريق اثنين من الذيل اختبار دقيق فيشر والذيل اثنين من يلكوكسون الرتبة مجموع اختبار في R
ملاحظة: اختبار استخدام فيشر الدقيق لمقارنة نسبة المواقع التكامل داخل الجينات RefSeq أو ضمن إطار من الجزر الدليل السياسي الشامل أو TSSs، ولكن استخدام يلكوكسون رتبة اختبار مبلغ لمقارنة التوزيع في كثافة الجينات المحيطة مواقع التكامل. يتوفر برنامج R في http://www.r-project.org/.
اثنين من الذيل اختبار فيشر بالضبط:
- باستخدام الأرقام محسوبة وفقا للتعليمات في خطوات 9.4 و 9.7، كرeate المصفوفات لكل المقارنة في مجال البحث من الحوادث لاحظ (التكامل ضمن الشرح أو ضمن إطار المحيطة تعليق توضيحي) مقابل المتبقية المواقع عن طريق تتبع هذا الأمر:
(annotation_of_interest <- مصفوفة (ج (SampleA # في، SampleA # المتبقي، SampleB # في، SampleB # المتبقي)، nrow = 2، dimnames = قائمة (ج ( 'المركز'، 'الباقي')، ج ( 'SampleA، "SampleB ')))) - حساب قيمة P للمقارنة من قبل ثنائي الطرف اختبار فيشر بالضبط مع الأمر التالي:
fisher.test (annotation_of_interest، بديل = 'two.sided') $ p.value
اثنين من الذيل اختبار Wilcoxon رتبة المبلغ: - إنشاء ملف txt المفصول فيها كل عمود يحتوي على اسم عينة في الخلية الأولى، تليها أقل من القيم كثافة الجينات لجميع المواقع التكامل في تلك المكتبة (تم الحصول عليها من ملف .bed ولدت في الخطوة 9.9). استيراد هذا الملف txt المفصول إلى R باستخدام الأمر التالي وغvigating إلى الدليل الملف الصحيح:
FILENAME <- as.data.frame (read.delim (file.choose ()، رأس = T، check.names = FALSE، وملء = TRUE، سبتمبر = ' ر')) - حساب قيمة P للمقارنة من قبل اختبار Wilcoxon رتبة مبلغ ثنائي الطرف مع الأمر التالي:
wilcox.test (FILENAME $ SampleA، FILENAME $ SampleB، = البديلة "two.sided"، وإرفاقها = F، بالضبط = T) $ p.value
ملاحظة: يمكن حساب قيم P فقط لتصل إلى (منخفضة للغاية) حد معين في R، سيتم إرجاع بواسطة برنامج بعدها الصفر. للحصول على عينات مختلفة على نطاق واسع أن تسفر عن P = 0 في البحث، وتقدير قيمة ف ك <2.2 × 10 -308.
11. دراسة أولية البيانات التسلسل لدليل على توسع نسيلي من الخلايا التي تحتوي على الحمض النووي الفيروسي المتكاملة
ملاحظة: هناك إمكانية صغيرة لدمج أكثر من واحد في نفس الإقليم الشمالي المحدد في الجينوم المرجعية. بدلا من ذلك، واحدة فيقد يصبح الحدث tegration الحالي بوفرة في تسلسل البيانات نتيجة لاستخدام PCR خلال إعداد مكتبة و / أو عن طريق خلية الازدواجية قبل إعداد الحمض النووي. وقد تميز التحليلات الأخيرة من الحمض النووي الجيني من المرضى المصابين بفيروس نقص المناعة البشرية هذه الاحتمالات من خلال تحديد نقاط صوتنة القص / رابط مرفق فريدة من نوعها (والتي يمكن أن تنشأ إلا قبل PCR) ضمن تسلسل الحمض النووي التي تحتوي على مواقع التكامل متطابقة 52-54. حاليا هناك نقاش حول ما إذا كان proviruses آوى داخل الخلايا توسيع متطابق بسلالة تساهم في الخزان الفيروسي الكامنة، وبالتالي فإنه من الأهمية الخاصة للتميز مستواها من التوسع عند دراسة مواقع التكامل في المرضى من البشر.
- على غرار الإجراء المذكورة في الخطوة 8.1، إنشاء ملفات .bed إدراج فاصل من قواعد تمتد، في هذه الحالة، 25 الإقليم الشمالي المصب من كل موقع التكامل الفريد (قواعد المنبع ليست ضرورية هنا). إنشاء ملف FASTA من هذه الملفات .bed (كما أشرنا فيالخطوة 8.1) باستخدام وظيفة fastaFromBed من BEDTools وبعد هذا الأمر:
fastaFromBed -فاي / الدليل / إلى / إشارة / الجينوم / -اسم -s -bed 25_base_pair_file.bed -fo output_file.fasta
ملاحظة: لتحسين خصوصية كل بحث فمن المستحسن لاستخراج لا يقل عن 25 الإقليم الشمالي المصب من كل موقع التكامل ليحلل التوسع نسيلي. - ويفضل استخدام برنامج نصي حسب الطلب، بحث ملف FASTA تسلسل البيانات الخام لجميع السلاسل التي تحتوي على تطابق تام إلى 25 الإقليم الشمالي المصب من كل موقع التكامل فريدة من نوعها، وتودع هذه المتتاليات في ملف جديد. تقليم LTR وتسلسل رابط من السلاسل الخام. دمج سلسلة PE يقرأ عن طريق تحويل يقرأ لتكملة العكسي، وتقليم LTR ورابط متواليات، ومن ثم تعيين سلاسل read2 إلى زوج من read1 إذا تشترك في سلاسل لا يقل عن 20 متداخلة الإقليم الشمالي.
- مسح نقاط التعلق رابط كل كتلة موقع التكامل. تصنيف كل التكامل ب "توسيع متطابق بسلالة و# 34؛ إذا نقاط التعلق رابط هي ≥3 بي بي بعيدا.
ملاحظة: يقرأ وقد وصفت بروتوكول للتحليل توسع نسيلي دون دمج تسلسل 52.
ملاحظة: تجزئة من الجينوم في نفس المكان المحدد من قبل صوتنة يؤدي إلى التقليل من مدى توسع نسيلي، وأساليب لتصحيح وقد وصفت 63،64 التحيز التجريبي الناتج.