Identifying Robust Single-Cell 3D Genome Domain Boundaries with CellRoB
يُعد CellRoB إطارًا حوسبيًا يتغلب على ندرة بيانات الجينوم ثلاثي الأبعاد أحادية الخلية لتحديد حدود نطاقات التفاعل الطولي (TAD) الشائعة والنادرة بقوة، مما يكشف عن عدم التجانس الهيكلي ذي الأهمية البيولوجية ويسمح بتمييز المجموعات الفرعية من الخلايا التي لا يمكن تمييزها من خلال البيانات النسخية وحدها.
داخل كل خلية من خلايا جسمك، لا يكون الحمض النووي (DNA) مجرد خيط متشابك وفوضوي، بل هو بنية مطوية بإحكام، منظمة تماماً مثل مكتبة حيث تُرتّب الكتب على رفوف محددة. هذا الترتيب ثلاثي الأبعاد أمر بالغ الأهمية لأنه يحدد أي الجينات يتم تفعيلها أو إيقافها، ليعمل فعلياً كدليل تعليمات الخلية. لقد عرف العلماء منذ زمن طويل أن هذا الجينوم ينقسم إلى أحياء متميزة تسمى "نطاقات" (domains)، تفصل بينها أسوار غير مرئية تُعرف بـ "الحدود" (boundaries). تضمن هذه الأسوار عدم اختلاط تعليمات جزء من الخلية مع تعليمات جزء آخر عن طريق الخطأ. ومع ذلك، عندما ينظر الباحثون إلى خلية واحدة بدلاً من حشد من الخلايا، تصبح الصورة ضبابية للغاية؛ إذ تكون البيانات شحيحة ومليئة بالفجوات لدرجة تجعل من المستحيل تقريباً رؤية أين تقف هذه الأسوار بالفعل، مما يجعل من الصعب فهم كيف تختلف الخلايا الفردية عن بعضها البعض أو كيف تقرر التحول إلى أنواع معينة من الأنسجة.
قام فريق من الباحثين في جامعة شيديان (Xidian University) بتطوير أداة حوسبة جديدة تسمى "CellRoB" لحل مشكلة عدم الوضوح هذه. فبدلاً من محاولة فرض صورة واضحة من البيانات الضوضائية، بنى الباحثون طريقة تحترم التباين الطبيعي بين الخلايا مع استمرارها في إيجاد الأنماط المشتركة. بدأوا بملء القطع المفقودة لخرائط الاتصال الجيني للخلايا الفردية، باستخدام معلومات من خلايا مماثلة لتقديم تخمينات مدروسة حول مكان الروابط المفترض وجودها. وبمجرد حصولهم على رؤية أوضح للهيكل الداخلي لكل خلية، حددوا الحدود المحتملة. وجاء الابتكار الحقيقي في الخطوة التالية: حيث قاموا بمحاذاة هذه الحدود عبر آلاف الخلايا، ليس من خلال إجبارها على التطابق التام، بل عبر إيجاد الموقع الأكثر احتمالاً لوجود حدٍ ما لمجموعة من الخلايا. سمح هذا النهج لهم بالتمييز بين الحدود الموجودة في معظم الخلايا، وتلك التي تظهر في عدد قليل فقط، وتلك الفريدة لخلية واحدة.
اختبر الباحثون طريقتهم على بيانات من خلايا دماغ الفئران وخلايا بشرية، وقارنوا نتائجهم بالتقنيات الحالية. ووجدوا أن "CellRoB" يمكنه تحديد الحدود بدقة واتساق أكبر بكثير، حتى عندما كانت البيانات تفتقر إلى عدد كبير من الاتصالات أو عندما كانت الدقة عالية للغاية، وصولاً إلى 5,000 زوج قاعدي. وعلى عكس الطرق الأخرى التي إما كانت تفقد الحدود النادرة أو تخلق الكثير من الحدود الزائفة كلما أصبحت البيانات أكثر تفصيلاً، ظل "CellRoB" مستقراً. لقد نجح في تحديد الحدود المشتركة بين معظم الخلايا، والتي غالباً ما ترتبط بالبروتينات الهيكلية الأساسية للخلية، لكنه كشف أيضاً عن طبقة خفية من التباين: وهي الحدود التي تظهر في جزء صغير فقط من الخلايا. لم تكن هذه الحدود النادرة مجرد ضجيج؛ بل كانت تحمل معلومات محددة حول هوية الخلية وارتبطت بجينات مسؤولة عن وظائف متخصصة، مثل تطور الجهاز العصبي.
ربما كان الاكتشاف الأكثر إثارة هو عندما طبق الفريق أداة "CellRoB" على مجموعة بيانات لأجنة الفئران تضمنت بيانات الهيكل الجيني والتعبير الجيني معاً. في هذه المجموعة من البيانات، ظهرت مجموعة من الخلايا تُعرف باسم خلايا "الأديم الظاهر" (epiblast) و"الخط البدائي" (primitive streak) كمجموعة واحدة متجانسة عند تحليل نشاطها الجيني وحده. ومع ذلك، عندما رسم الباحثون خرائط لهذه الخلايا بناءً على قوة نطاقات الهيكل، انقسمت المجموعة إلى مجموعتين فرعيتين متميزتين. لم تكن هاتان المجموعتان منفصلتين بسبب وجود أو غياب الحدود، بل بسبب اختلافات دقيقة في مدى قوة تلك الحدود. أظهرت الخلايا في إحدى المجموعتين الفرعيتين سمات جينية تشير إلى أنها تتجه لتصبح نسيجاً عصبياً، بينما بدت المجموعة الأخرى منحازة نحو سلالات غير عصبية. يشير هذا الاكتشاف إلى أن قوة الأسوار الهيكلية للجينوم هي عامل حاسم، وإن كان مغفولاً عنه سابقاً، في تحديد مصير الخلية. ومن خلال الكشف عن هذه الاختلافات الهيكلية منخفضة التكرار وقياس قوة البنية الجينية، يوفر "CellRoB" طريقة جديدة لرؤية التنوع الخفي داخل الخلايا الحية، مما يقدم رؤية أوضح لكيفية تطور الكائنات المعقدة من نقطة بداية واحدة.
ملخص تقني: تحديد حدود نطاقات الجينوم ثلاثي الأبعاد القوية في الخلية الواحدة باستخدام CellRoB
بيان المشكلة توفر تقنيات الجينوم ثلاثي الأبعاد في الخلية الواحدة القدرة على استقصاء تنظيم الكروماتين بما يتجاوز متوسطات المجموعات، مما يكشف عن التباين الهيكلي. ومع ذلك، فإن الندرة الشديدة والضجيج التقني المتأصل في خرائط الاتصال في الخلية الواحدة (مثل scHi-C، وDip-C، وscMicro-C) يحجبان حدود النطاقات المرتبطة بالارتباط الطوبولوجي (TAD). تؤدي هذه الندرة إلى تحديين رئيسيين:
موثوقية الكشف: غالبًا ما تظل إشارات الحدود الضعيفة غير مكتشفة، كما أن عدم اليقين الموضعي يعقد المقارنات بين الخلايا.
القابلية للمقارنة: حتى عند اكتشاف الحدود، فإن التفاوتات الموضعية بين الخلايا تعيق المقارنة المنهجية لظهور الحدود وقوتها. تركز النهج الحسابية الحالية (مثل scHiCluster، وHigashi، وDeTOKI، وDeDoc2) غالبًا على استكمال الاتصال أو كشف النطاقات، لكنها تفشل في إنشاء تمثيل ثابت للحدود عبر الخلايا. علاوة على ذلك، فإن تجميع البيانات في خرائط "شبه مجمعة" (pseudo-bulk) يعزز الإشارات المتكررة ولكنه يحجب الحدود منخفضة التكرار التي تشترك فيها نسبة ضئيلة فقط من الخلايا، مما يحد من توصيف الاختلافات الهيكلية النادرة.
المنهجية: إطار عمل CellRob CellRoB هو إطار عمل حسابي مصمم لتحديد ومقارنة حدود نطاقات الكروماتين في الخلية الواحدة بقوة مع الحفاظ على التباين الخاص بكل خلية. يتكون خط الإنتاج من أربع مراحل رئيسية:
استكمال الاتصال (Contact Imputation):
يستخدم الأسلوب أولاً Higashi لتوليد مصفوفة تشابه بين الخلايا (W) بناءً على خرائط اتصال مقسمة إلى فئات بحجم 500 كيلوبايت.
لكل خلية، يتم الاحتفاظ فقط بأفضل k=15 من أقرب الجيران.
يتم استكمال خريطة الاتصال الخام (40 كيلوبايت) (Hiraw) عبر المتوسط المرجح للخرائط الخام لجيرانها: Hiimp=1+∑j=1NWijHiraw+∑j=1NWijHjraw
الكشف الأولي لحدود الخلية الواحدة:
مستوحى من التجميع القائم على الكثافة، يحدد CellRoB مرشحي الحدود في الخلايا الفردية باستخدام مؤشرين لكل فئة جينومية i:
قوة العزل (ρ(i)): تُحسب كفرق طبيعي بين مجموع الاتصالات داخل النطاق وخارج النطاق ضمن نافذة (w=15 فئة عند 40 كيلوبايت).
المسافة الدنيا (δ(i)): المسافة إلى أقرب فئة ذات قوة عزل أعلى ضمن نطاق M (20 ميجابايت).
يتم تطبيع هذه المؤشرات ودمجها في درجة η(i)=ρ′(i)δ′(i).
يتم تحديد نقطة انعكاس r لموازنة الدرجات العالية مع الرتب المنخفضة. الفئات التي تحقق η(i)>αη(r) (حيث α يتناسب مع حجم الفئة) يتم اختيارها كحدود أولية.
المحاذاة عبر الخلايا (Cross-Cell Alignment):
لمعالاة التفاوتات الموضعية، يتم تجميع إشارات الحدود الخام من جميع الخلايا وتنعيمها باستخدام نواة موحدة.
تقوم الأودية المحلية في الإشارة المنعمة بتقسيم الفضاء الجينومي إلى مجموعات.
لكل مجموعة، يتم تحديد فئة إجماع (consensus bin) عن طريق تقليل خسارة درجة العزل التربيعية عبر جميع الخلايا التي لديها حد في تلك المجموعة. هذا يؤدي إلى محاذاة الحدود إلى نظام إحداثيات مشترك دون فرض قرار ثنائي "موجود/غير موجود" لكل خلية.
تصنيف الحفظ (Conservation Classification):
يتم تصنيف الحدود المحاذية بناءً على تكرار حدوثها (fi) عبر مجتمع الخلايا:
نادرة (Rare):1%≤fi<10%
محفوظة ضعيفًا (Weakly Conserved):10%≤fi<50%
محفوظة بقوة (Strongly Conserved):fi≥50%
النتائج الرئيسية
الأداء على بيانات Dip-C: عند الاختبار على بيانات Dip-C لـ mouse brain (كروموسوم 6، 40 كيلوبايت)، أظهر CellRoB أداءً فائقًا مقارنة بـ HiCS، وHigashi، وDeTOKI، وDeDoc2.
الصلاحية البيولوجية: أظهرت الحدود المتوقعة أعلى إثراء للعوامل المعيارية (CTCF، وH3K4me3، والجينات التدبيرية/housekeeping genes).
الاتساق: حقق CellRoB أعلى مؤشر جاكارد (JI) ومقياس التوافق (MoC) عند مقارنة الحدود بين الخلايا المتشابهة هيكليًا.
المتانة: حافظ على أعلى درجات F1 عبر مستويات الضجيج (من 5% إلى 40% من تساقط الاتصال) وأظهر تكرارية عالية في مجموعات البيانات المكررة.
الأداء عالي الدقة (scMicro-C): تم التقييم على بيانات scMicro-C البشرية لـ GM12878 عبر دقات تتراوح من 40 كيلوبايت وصولًا إلى 5 كيلوبايت.
على عكس الطرق الأخرى (HiCS، وDeDoc2، وDeTOKI) حيث يزداد عدد الحدود المتوقعة بشكل أسي مع تناقص حجم الفئة، ظل CellRoB مستقرًا.
حقق أعلى مقياس توافق (MoC) وأقوى إثراء للميزات البيولوجية حتى عند دقة 5 كيلوبايت.
كشف تحليل إثراء المصطلحات الجينية (GO enrichment) للجينات المرتبطة بفئات الحدود المختلفة في خلايا طور G1 عن ملفات تعريف وظيفية متميزة: ارتبطت الحدود النادرة بتطور الجهاز العصبي، بينما ارتبطت الحدود المحفوظة بقوة بانقسام الخلايا.
توصيف الحدود النادرة:
نجح CellRoB في اكتشاف الحدود الموجودة في مجموعة فرعية صغيرة من الخلايا، والتي عادة ما يتم إغفالها في التحليلات المجمعة.
حملت الحدود النادرة معلومات مميزة لنوع الخلية (مثل الإثراء لجينات التمايز العصبي) وأظهرت بصمات فوق جينية مختلفة مقارنة بالحدود المحفوظة.
حل المجموعات الفرعية التنموية (مجموعة HiRES):
عند تطبيقه على بيانات الجنين للفئرون متعددة الأنماط (Hi-C + RNA-seq)، استطاع CellRoB حل مجموعتين فرعيتين داخل مجموعة "الابللاست" (Epiblast) و"الخط البدائي" (Primitive Streak) اللتين لم يكن بالإمكان التمييز بينهما في التضمينات المستمدة من RNA.
كشف التحليل التفاضلي أن هذه الأنواع الفرعية لم تكن متميزة بوجود/غياب الحدود (ميزات ثنائية)، بل من خلال التباينات الكمية في قوة الحدود (درجات العزل).
أظهر النوع الفرعي 1 إثراءً لجينات إغلاق الأنبوب العصبي، بينما أظهر النوع الفرعي 2 إثراءً لعمليات نمو غير عصبية، مما يشير إلى مسارات تنموية متميزة.
الأهمية والادعاءات يزعم البحث أن CellRoB يوفر إطار عمل قويًا لربط التحديد الموثوق لحدود الخلية الواحدة باكتشاف الاختلاف الهيكلي منخفض التكرار. وتتمثل مساهماته الرئيسية في:
تحسين القابلية للمقارنة: من خلال نمذجة مناطق الحدود المتكررة ومحاذاتها عبر تقليل درجة العزل، فإنه يتيح المقارنة المنهجية لحدوث الحدود وقوتها عبر الخلايا.
استعادة الميزات النادرة: يتيح تعريف وتوصيف الحدود النادرة التي تحجبها عمليات التجميع السكاني، مما يكشف عن التباين الهيكلي ذي المعلومات البيولوجية.
الميزات الهيكلية الكمية: توضح الدراسة أن قوة الحدود الكمية، وليس مجرد وجود الحدود الثنائي، يمكن أن تعمل كميزة هيكلية لحل التباين الخلوي والحالات التنموية التي لا تستطيع البيانات النسخية (transcriptomic) فصلها وحدها.
يشير المؤلفون إلى أنه بينما يحسن CellRoB الاتساق، فإنه يعتمد على تغطية المدخلات وخيارات المعالجة المسبقة (الاستكمال، وحجم الفئة). كما يقرون بأن تصنيف الحدود النادرة يعتمد على عتبات التكرار، وأن هناك حاجة للتحقق التجريبي المستقل لتحديد الأهمية الفيزيائية للحدود المستنتجة، نظراً لعدم وجود معيار ذهبي نهائي لتعليقات حدود الخلية الواحدة.