$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
يعتمد بروتوكول Capture Hi-C الموصوف على إعداد قالب 3C على مستوى الجينوم باستخدام قاطع رباعي القواعد (DpnII). يتم الحصول على الإثراء اللاحق لشظايا الربط عبر المنطقة الجينومية ذات الأهمية عن طريق تهجين مجموعة من مجسات الحمض النووي الريبي التبليطية والتقاطها القائم على الستربتافيدين وفقا لنظام التخصيب المستهدف المستخدم في هذه الدراسة (الشكل 1). تم اختيار مجسات الحمض النووي الريبي البيوتينيلية لأنها تظهر تقاربا أكثر إحكاما لأهدافها مقارنة بمجسات الحمض النووي52,60. ثم يتم فهرسة المكتبات التي تم التقاطها وتجميعها من أجل تسلسل متعدد الإنتاجية عالي الإنتاجية. يمكن تصور بيانات التقاط Hi-C كخرائط تفاعل Hi-C عالية الدقة ولكن أيضا كخرائط اتصال أحادية نقطة عرض تشبه 4C لتصور تفاعلات التسلسلات الأصغر على وجه التحديد مثل المروجين أو المعززات داخل المنطقة الملتقطة بأكملها. يظهر سير عمل البروتوكول في الشكل 4. يتم عرض ضوابط جودة التسلسل المسبق في الشكل 2 وتشمل تقييم الهضم السليم وإعادة ربط قالب 3C والقص والتنقية الفعالة عبر الخطوات المختلفة للبروتوكول. من المتوقع أن يعمل الحمض النووي لقالب 3C المقطوع بين 150 إلى 700 نقطة أساس ، ولا ينبغي الكشف عن أي إثراء للشظايا >2 كيلو بايت. خلال الخطوات التالية ، يتم تنفيذ العديد من خطوات تنظيف الحمض النووي القائمة على الخرزة واختيار الحجم ، أولا بعد القص ، ثم بعد PCRs قبل الالتقاط وما بعد الالتقاط. تظهر المكتبات التي تم تنظيفها ملف تعريف إثراء شظية مميز كما هو موضح في محلل حيوي عالي الحساسية للحمض النووي (الشكل 2). يزداد متوسط حجم الجزء على مدار إعداد المكتبة بسبب ربط المحولات والتسلسل وفهرسة البادئات. يتم الحصول على ضوابط الجودة بعد التسلسل عبر Hi-C Pro وتظهر في الشكل 3. تم اقتراح العديد من تطبيقات برامج المعلوماتية الحيوية المختلفة لمعالجة وتحليل البيانات الشبيهة ب 3C. من بينها ، يعد خط أنابيب HiC-Pro أحد أكثر الحلول شيوعا ، مما يسمح بمعالجة بيانات التسلسل الخام إلى خرائط الاتصال النهائية بدقةمختلفة 55. يستخدم HiC-Pro استراتيجية رسم خرائط من خطوتين لمحاذاة قراءات التسلسل على الجينوم المرجعي. ثم يتم إعادة بناء منتجات 3C وتصفيتها لإزالة أزواج الاتصال غير المفيدة وإنشاء خرائط الاتصال. بالإضافة إلى ذلك ، فإنه قادر على استخدام قائمة من الأشكال المتعددة المعروفة لإجراء تحليل خاص بالأليل وفصل جهات الاتصال القادمة من الأليلين الأبويين في خرائط اتصال متميزة. في الآونة الأخيرة ، تم تضمين HiC-Pro وتوسيعه في إطار عمل nf-core (nf-core-hic) ، مما يوفر خط أنابيب مدفوعا بالمجتمع قابل للتطوير وقابل للتكراربدرجة كبيرة 61,62.
لالتقاط الفأر Xic ، تم تصميم مجموعة من 28,913 مسبار RNA تبليط 3 ميجا بايت من كروموسوم X. تتضمن هذه المنطقة اللاعب الرئيسي في XCI ، الجين الطويل غير المشفر Xist ، والمشهد التنظيمي المعروف ~ 800 كيلو بايت (الشكل 5). تنقسم هذه المنطقة ~ 800 كيلو بايت إلى اثنين من TADs: أحدهما يتضمن مروج Xist ومنظماته الإيجابية المعروفة (أي النصوص غير المشفرة Ftx و Jpx و Xert والجين المشفر للبروتين Rnf12) ، و TAD المجاور الذي يشمل منظمات رابطة الدول المستقلة السلبية ل Xist (أي نسخته المضادة للمعنى Tsix ، وعنصر المحسن Xite ، والنص غير المشفر Linx) (للمراجعة 44 ، 45).
من خلال تطبيق بروتوكول Capture Hi-C الموصوف على Xic ، تم الحصول على التنظيم الطوبولوجي لهذا الموضع بدقة غير مسبوقة (الشكل 6 والشكل 7). هذا واضح بشكل خاص عند مقارنة ملف تعريف Capture Hi-C ب 5C47 المنشور مسبقا (الشكل 6 والشكل 7 ؛ الجدول التكميلي 1) و Hi-C61 (الشكل 6 والشكل 7 ؛ الجدول التكميلي 1) ملفات تعريف. على سبيل المثال ، تكون هياكل TAD الفرعية أكثر وضوحا - من الواضح أن TAD الذي يحتوي على مروج Xist (Xist-TAD) ينقسم إلى مجالين أصغر (الشكل 6A ، رأس السهم الأزرق). في السابق ، كان من الممكن "تخمين" هذا بصريا فقط من ملف تعريف 5C (الشكل 6B) ، وإن كان اكتشاف حدود في هذه المنطقة باستخدام خوارزمية درجة العزل. وبالمثل ، فإن دقة ملف تعريف Capture Hi-C تسمح بتحديد مجالين أصغر في TAD المجاور (الشكل 6A ، B) ، والذي يحتوي على مروج موضع Tsix ( Tsix-TAD ) ؛ لم يتحقق ذلك من قبل مع 5C (الشكل 6B). وتجدر الإشارة إلى أن الحدود الطوبولوجية التي تحددها درجة العزل من بيانات Capture Hi-C و 5C يتم اكتشافها بشكل عام في مواقع مختلفة قليلا وبقوة نسبية مختلفة.
علاوة على ذلك ، فإن هياكل TAD الفرعية الأخرى مثل حلقات الاتصال مرئية بوضوح من بيانات Capture Hi-C ، مثل الحلقة بين Xist و Ftx (الشكل 7A) ، التي تم تحديدها مسبقا مع Capture-C63 ، والحلقة بين Xist و Xert (الشكل 7B) ، تم تحديدها مؤخرا باستخدام بروتوكول مماثل ل Capture Hi-C48. يمكن أيضا تعيين جهات الاتصال الأخرى بشكل أكثر دقة بسبب زيادة دقة ملفات تعريف Capture Hi-C ، مثل تلك التي تشكل نقاط الاتصال المعروفة داخل Tsix-TAD بين مواقع Linx و Chic1 و Xite (الشكل 7A).
بالمقارنة مع بيانات Hi-C الموضحة في الشكل 7 ، سمح Capture Hi-C بزيادة أربعة أضعاف في الدقة ، ومع ذلك فقد تطلب فقط ربع عمق التسلسل (أي 126 M يقرأ مقابل 571 M) (الجدول التكميلي 1). تسمح هذه الزيادة في الدقة باكتشاف subTADs وتفاعلات الحلقات التي لا يمكن اكتشافها بواسطة Hi-C عند عمق التسلسل الموضح في الشكل 6 والشكل 7. وبالتالي ، فإن البروتوكول الموصوف ل Capture Hi-C يسمح بتوصيف أكثر تفصيلا وعالي الدقة لمنطقة جينومية كبيرة ذات أهمية ، مقارنة بالأساليب السابقة.

الشكل 1: تصميم المسبار. تمثيل تخطيطي للاستراتيجية المستخدمة لتصميم المسبار. تم اختيار مناطق 300 نقطة أساس في المنبع والمصب لكل موقع تقييد DpnII عبر المنطقة المستهدفة 3 ميجا بايت وتم تبليطها بمجسات الحمض النووي الريبي البيوتينيلية المتداخلة. تظهر إحدى هذه المناطق المحددة، chrX: 102,474,805-102,475,500. لا يسمح بأكثر من 40 قاعدة من التسلسلات المتكررة في كل مسبار. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 2: التقاط ضوابط جودة التسلسل المسبق Hi-C . (أ) مثال تمثيلي لضوابط جودة قالب 3C. تم تحميل 200 نانوغرام من الحمض النووي على هلام أغاروز 1٪. حارة 1: 1 كيلو بايت سلم. الحارة 2: يعمل الكروماتين غير المهضوم والمتشابك والسليم كشريط حاد عند >10 كيلو بايت. الحارة 3: يعمل الكروماتين المتصالب المهضوم DpnII كمسحة يتراوح حجمها بين 1 كيلو بايت إلى 3 كيلو بايت. الحارة 4: مكتبة أو قالب 3C النهائي ؛ يتم إعادة ربط الأطراف الحرة لشظايا الحمض النووي المتشابكة المهضومة. لا يمكن اكتشاف مسحة الحمض النووي ذات الحجم الجزيئي المنخفض تقريبا ، ويتم اكتشاف منتج الربط كنطاق >10 كيلو بايت. (ب) أمثلة تمثيلية لملامح الحمض النووي للمحلل الحيوي عالي الحساسية. أعلى اليسار: مكتبة 3C تم قصها بنجاح تظهر توزيعا لحجم الجزء بين 150 نقطة أساس و 700 نقطة أساس. أعلى اليمين: مكتبة 3C غير مرضية. تم الكشف عن الحمض النووي غير المقطوع على أنه إثراء واسع للشظايا >2 كيلو بايت. (ج) أسفل اليسار: عينة الحمض النووي المنفصمة بعد اختيار حجم الجانب الأيسر 1: 1 باستخدام حبات SPRI. يتم إثراء شظايا ~ 300 bp. الوسط السفلي: ملف تعريف PCR قبل الالتقاط بعد ربط المحولات ذات النهاية المزدوجة وفقا لبروتوكول الشركة المصنعة. أسفل اليمين: مكتبة Capture Hi-C النهائية بما في ذلك المحولات والتسلسل وبادئات الفهرسة للتسلسل متعدد الإرسال. الاختصارات: bp = أزواج القواعد ، FU = وحدة مضان تعسفية. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 3: التقاط ضوابط الجودة بعد التسلسل Hi-C باستخدام HiC-Pro . (أ) مثال على معدل رسم الخرائط على الجينوم المرجعي للرفيق الأول لأزواج التسلسل. يمثل الجزء الأزرق الفاتح القراءات المحاذية بواسطة HiC-Pro وتمتد عبر تقاطع ربط. وبالتالي يمكن استخدام هذا المقياس للتحقق من صحة خطوة الربط التجريبية. (ب) بمجرد محاذاة زملاء التسلسل على الجينوم ، يتم الاحتفاظ بأزواج القراءة المحاذاة بشكل فريد فقط للتحليل. (ج) يتم استبعاد الأزواج غير الصالحة (باللون الأحمر) مثل الطرف المتدلي أو الدائرة الذاتية أو إعادة الربط من التحليل. يعد جزء الأزواج الصالحة مؤشرا جيدا على كفاءة الربط والقائمة المنسدلة. (د) يمكن تقسيم الأزواج الصالحة إلى ملامسات داخل الكروموسومات وبينها وبين ملامسات قصيرة/طويلة المدى. يتم تجاهل أزواج القراءة المكررة التي من المحتمل أن تمثل عناصر تفاعل البوليميراز المتسلسل من التحليل. (ه) بالنسبة للتحليل الخاص بالأليل ، يبلغ HiC-Pro عن عدد القراءات الأليلية التي يدعمها إما واحد أو اثنين من الشركاء لكل جينوم أبوي (أي C57BL / 6J x CASTEi / J). من المتوقع نفس الجزء من القراءات المخصصة لأليل الأم والأب. (F) أخيرا ، يتم تحديد الأزواج الصالحة المتداخلة مع منطقة الالتقاط فقط لإنشاء خرائط الاتصال. تمثل أزواج الالتقاط والالتقاط جهات اتصال داخل المنطقة المستهدفة ، بينما تتضمن أزواج الالتقاط والمراسل التفاعل بين المنطقة المستهدفة والمنطقة خارج الهدف. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 4: سير عمل بروتوكول التقاط Hi-C. تمثيل تخطيطي لخطوات البروتوكول المختلفة. لإنشاء قالب 3C على مستوى الجينوم ، يتم ربط الكروماتين أولا بالفورمالديهايد ثم يتم هضمه باستخدام إنزيم تقييد DpnII. ثم يتم إعادة ربط نهايات الحمض النووي الحرة ، ويتم عكس الارتباط المتقاطع ، ويتم تنقية الحمض النووي. لإثراء الشظايا التي تشمل المنطقة المستهدفة ، يتم تهجين مجموعة من مجسات الحمض النووي الريبي البيوتينيلاتي إلى قالب 3C ويتم التقاطها بواسطة السحب بوساطة الستربتافيدين لأسفل. تتم معالجة مكتبات الالتقاط للتسلسل المتعدد ، ويتم تحديد أجزاء الربط الصالحة لاستنتاج تكرار اتصالات الكروماتين عبر الهدف ، والتي يتم تصورها كخرائط تفاعل عالية الدقة. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 5: نظرة عامة على المنطقة التي تضم Xic على كروموسوم X للفأر. تمثيل تخطيطي لكروموسوم الماوس X وتكبير المنطقة الملتقطة 3 ميجا بايت (ChrX: 102,475,000-105,475,000). تتضمن المنطقة المستهدفة ~ 800 كيلو بايت من الحمض النووي المقابل ل Xic ، الموضع التنظيمي الرئيسي ل XCI. يتضمن Xic الجينات الطويلة غير المشفرة ، Xist ، وهو لاعب رئيسي في XCI ، والمشهد التنظيمي. تظهر المنظمات الإيجابية ل Xist باللون الأخضر ، والمنظمين السالبين باللون الأرجواني. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 6: التقط خرائط تفاعل Hi-C و 5C و Hi-C عبر المنطقة الملتقطة 3 ميجا بايت. (أ) التقاط خريطة تفاعل Hi-C لهدف 3 ميجا بايت يشمل الماوس Xic بدقة 10 كيلو بايت (هذه الدراسة). (B) خريطة تفاعل 5C لنفس المنطقة المستهدفة كما في A بدقة 6 كيلو بايت (البيانات المعاد معالجتها من47). المناطق المتكررة غير المدرجة في التحليلات مقنعة باللون الأبيض. تتطلب بيانات 5C معالجة المعلوماتية الحيوية الخاصة بها (انظر47). بعد التنظيف والمحاذاة ، يتم تثبيت خرائط 5C بدقة التمهيدي باستخدام وسيط قيد التشغيل (نافذة = 30 كيلو بايت ، الخطوة = 5) للوصول إلى دقة نهائية تبلغ 6 كيلوبايت. (C) خريطة تفاعل Hi-C لنفس المنطقة الجينومية كما في A و B بدقة 40 كيلو بايت (البيانات المعاد معالجتها من64). تم إنشاء جميع خرائط التفاعل من ESCs الماوس. تم حساب درجة العزل باستخدام أدوات التبريد ويتم تمثيلها كرسوم بيانية مع الحد الأدنى للعزل عند حدود TAD. تظهر حدود TAD كخطوط عمودية أسفل الخريطة. يشير ارتفاع كل خط إلى قوة الحدود. تظهر الجينات في صورة أسهم تشير إلى اتجاه النسخ. تتم الإشارة إلى حدود TAD الفرعية التي يتم اكتشافها بشكل حصري أو أكثر دقة في خرائط Capture Hi-C بواسطة رؤوس الأسهم الأرجواني والأزرق ل TADs الفرعية في Tsix و Xist TADs ، على التوالي. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 7: التقط خرائط تفاعل Hi-C و 5C و Hi-C عبر 1 ميجا بايت داخل المنطقة الملتقطة. (أ) التقاط خريطة تفاعل Hi-C للمنطقة الجينومية 1 ميجا بايت التي تشمل الفأر Xic بدقة 5 كيلو بايت (هذه الدراسة). (ب) خريطة تفاعل 5C لنفس المنطقة الجينومية كما في A. بدقة 6 كيلوبايت (تمت إعادة معالجة البيانات من47). المناطق المتكررة غير المدرجة في التحليلات مقنعة باللون الأبيض. وتجدر الإشارة إلى أن بيانات 5C تتطلب معالجة المعلوماتية الحيوية الخاصة بها (انظر47). بعد التنظيف والمحاذاة ، يتم تثبيت خرائط 5C بدقة التمهيدي باستخدام وسيط قيد التشغيل (نافذة = 30 كيلو بايت ، الخطوة = 5) للوصول إلى دقة نهائية تبلغ 6 كيلوبايت. (C) خريطة تفاعل Hi-C لنفس المنطقة الجينومية كما في A و B من Hi-C بدقة 20 كيلو بايت (البيانات المعاد معالجتها من64). تم إنشاء جميع خرائط التفاعل من mESCs. تم حساب درجة العزل باستخدام أدوات التبريد ويتم تمثيلها كرسوم بيانية مع الحد الأدنى للعزل عند حدود TAD. تظهر حدود TAD كخطوط عمودية أسفل الخريطة. يشير ارتفاع كل خط إلى قوة الحدود. تظهر الجينات في صورة أسهم تشير إلى اتجاه النسخ. يشار إلى حلقات التلامس التي يتم اكتشافها بشكل حصري أو أكثر دقة في Capture Hi-C بعلامات نجمية أرجوانية وزرقاء للحلقات في Tsix و Xist TADs ، على التوالي. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
الجدول التكميلي 1: إحصائيات ما بعد التسلسل لمجموعات البيانات المستخدمة في هذه المخطوطة: Capture Hi-C (هذه الدراسة) ، Hi-C64 ، و 5C47. الرجاء الضغط هنا لتنزيل هذا الملف.