EISCA and EISTA: Full-Spectrum Pipelines for Single-Cell and Spatial Transcriptomics Analysis
تقدم الورقة البحثية EISCA وEISTA، وهما مسارات عمل (pipelines) مبنية على إطار عمل nf-core، وتتميز بأنها معيارية، ونمطية، وقابلة للتوسع، وتوفر حلولاً شاملة ومتكاملة لتحليل بيانات النسخ الوراثي أحادي الخلية والنسخ الوراثي المكاني عالي الدقة، على التوالي، لتسريع الاكتشاف البيولوجي القابل للتكرار.
المؤلفون الأصليون:Wu, H., Lister, A., Macaulay, I. C., Long, K., Uauy, C., Lan, Y., Wickham, G. J., Swarbreck, D., Videm, P., Stubbs, A., Soranzo, N., de Waard-van Baardwijk, M., Nilchi, A. N., Papatheodorou, I.
المؤلفون الأصليون: Wu, H., Lister, A., Macaulay, I. C., Long, K., Uauy, C., Lan, Y., Wickham, G. J., Swarbreck, D., Videm, P., Stubbs, A., Soranzo, N., de Waard-van Baardwijk, M., Nilchi, A. N., Papatheodorou, I.
تتكون الحياة من خلايا، وهي الوحدات الحية الصغيرة التي تشكل كل نبات وحيوان. لعقود من الزمن، درس العلماء هذه الخلايا عن طريق طحن الأنسجة وتحويلها إلى حساء، وفصل الخلايا الفردية، وقراءة التعليمات الجينية بداخلها. هذا النهج، المعروف باسم تسلسل الخلية الواحدة، كشف أن الأنسجة ليست كتلًا متجانسة من أجزاء متطابقة، بل هي فسيفساء معقدة من أنواع عديدة ومختلفة من الخلايا التي تعمل معًا. ومع ذلك، فإن هذه الطريقة تشوبها ثغرة كبيرة: فمن خلال تفتيت النسيج، يتم تدمير الخريطة. يستطيع العلماء رؤية ماهية الخلايا، لكنهم يفقدون المعلومات الحاسمة حول مكان جلوس تلك الخلايا وكيف كانت مرتبة في الببنية الأصلية. ولإصلاح ذلك، تم تطوير تقنية أحدث تسمى "الترانسكريبتوميات المكانية". تسمح هذه الطريقة للباحثين بقراءة التعليمات الجينية بينما لا تزال الخلايا في مواقعها الأصلية، مما يحافظ على المخطط المعماري للنسيج.
التحدي الآن هو أن هذه التقنيات تولد كميات هائلة من البيانات المعقدة التي يصعب تحليلها. يمكن لتجربة واحدة أن تنتج ملايين النقاط البيانية، وفهمها يتطلب سلسلة طويلة من الخطوات الحاسوبية، بدءًا من تنظيف الإشارات الخام وصولاً إلى إيجاد الأنماط وتجميع الخلايا المتشابهة معًا. وحتى الآن، لم تكن هناك طريقة معيارية واحدة للقيام بذلك لجميع أنواع البيانات، مما ترك الباحثين يبنون أدواتهم الخاصة لكل مشروع جديد. وهذا النقص في وجود نظام موحد يجعل من الصعب مقارنة النتائج بين المختبرات المختلفة أو الوثوق في إمكانية تكرار النتائج.
ولحل هذه المشكلة، قام فريق من البحث في معهد إيرلهام (Earlham Institute) في المملكة المتحدة بإنشاء برنامجين حاسوبيين جديدين، أُطلق عليهما اسم EISCA وEISTA. فكر في هذين البرنامجين كخطوط تجميع مؤتمتة ومعيارية للبيانات البيولوجية. صُمم برنامج EISCA للأسلوب القديم المتمثل في تحليل الخلايا التي تم فصلها عن أنسجتها، بينما بُني برنامج EISTA خصيصًا للأسلوب الأحدث الذي يحافظ على الخلايا في أماكنها الأصلية. كلا البرنامجين مبنيان على إطار عمل مرن يسمح للعلماء بتشغيل التحليل بأكمله من البداية إلى النهاية بأمر واحد، أو التوقف عند أي خطوة لفحص النتائج وتعديل الإعدادات. وهما يتوليان كل شيء، بدءًا من البيانات الخام الأولية وصولًا إلى القائمة النهائية لأنواع الخلايا وتفاعلاتها، وينتجان تقارير واضحة توضح بدقة شكل البيانات في كل مرحلة.
اختبر الباحثون هذه الأدوات على مشكلات بيولوجية حقيقية لإثبات فعاليتها. أولًا، استخدموا EISTA لدراسة أوراق نبات "أرابيدوبسيس ثاليانا" (Arabidopsis thaliana) الشائع بعد إصابته ببكتيريا ضارة. ومن خلال تمرير البيانات عبر المسار البرمجي، نجح البرنامج في تحديد مجموعات محددة من الخلايا التي كانت تقاوم العدوى. وقد أظهر أن هذه الخلايا النشطة مناعيًا لم تكن منتشرة عشوائيًا، بل كانت تتجمع في نقاط محددة داخل نسيج الورقة. كما حدد البرنامج بدقة الجينات التي كانت هذه الخلايا تفعلها للدفاع عن النبات، مما كشف عن استجابة منظمة وواضحة للتهديد، وهو ما كان العلماء يشتبهون فيه ولكن كان من الصعب رسم خريطته بهذه الدقة من قبل.
بعد ذلك، طبق الفريق برنامج EISCA لدراسة عينات دم بشرية من مرضى يعانون من تعفن الدم (Sepsis)، وهو رد فعل مهدد للحياة تجاه العدوى. حلل البرنامج خلايا الدم من المرضى الذين نجوا من المرض وأولئك الذين لم ينجوا. وأكد أن المرض يسبب تحولًا دراماتيكيًا في أنواع الخلايا التي تدور في الدم. وتحديدًا، وجد البرنامج أن المرضى الذين لم ينجوا عانوا من انخفاض حاد في عدد الخلايا البائية (B cells)، وهي خلاسة ضرورية للمناعة طويلة الأمد، وتراكم خطير للصفائح الدموية وخلايا مرتبطة بالإجهاد. وقد طابقت هذه النتائج نتائج الدراسة الأصلية التي ولدت البيانات، مما أثبت أن المسار البرمجي الجديد يمكنه إعادة إنتاج الرؤى البيولوجية المعقدة بشكل موثوق دون الحاجة إلى فريق من الخبراء لتوجيهه يدويًا في كل خطوة.
تكمن القوة الحقيقية لهذه الأدوات في قدرتها على جعل العلوم المتقدمة متاحة للجميع. نظرًا لأن البرامج مؤتمتة ومعيارية، يمكن للباحث الذي يمتلك خبرة أقل في البرمجة الحاسوبية الآن إجراء تحليل متطور والحصول على صورة واضحة لبياناته فورًا. ولا تجبر هذه الأدوات العلماء على الاختيار بين السرعة والمرونة؛ إذ يمكنهم تشغيل التحليل الكامل للحصول على نظرة عامة سريعة أو الغوص في خطوات محددة لاستكشاف أسئلة جديدة. ومن خلال توفير طريقة مشتركة وموثوقة لمعالجة هذه المجموعات الضخمة من البيانات، يزيل EISCA وEISTA عقبة رئيسية في علم الأحياء الحديث. فهما يسمحان للعلماء بالتركيز بشكل أقل على صراعهم مع البيانات، وبشكل أكبر على فهم الأنظمة الحية التي تصفها، مما يضمن أن الاكتشافات التي تتم في مختبر واحد يمكن الوثوق بها والبناء عليها من قبل الآخرين حول العالم.
بيان المشكلة أحدثت تقنيات تسلسل الحمض النووي الريبي أحادي الخلية (scRNA-seq) والتحليل النسخي المكاني ثورة في فهم التباين الخلوي والتنظيم النسيجي. ومع ذلك، لا تزال التعقيدات التحليلية لهذه التقنيات تشكل عائقًا كبيرًا. ورغم وجود العديد من الخوارزميات وأدوات التحليل (مثل Seurat وScanpy وSquidpy)، إلا أن هناك حاجة إلى مسارات عمل معيارية وشاملة (end-to-end) يمكنها التعامل مع كامل نطاق التحليل — بدءًا من معالجة البيانات الأولية وصولاً إلى الرؤى المتقدمة في المراحل النهائية — مع الحفاظ على المرونة اللازمة للبحث الاستكشافي. وتحديدًا، تفتقر منصات التحليل النسخي المكاني القائمة على التصوير مثل Vizgen MERFISH و10x Xenium إلى مسارات عمل شاملة ومتبعة على نطاق واسع تدمج معالجة البيانات الأولية مع التحليل الإحصائي المكاني.
المنهجية يقدم المؤلفون مسارين معياريين، EISCA (تحليل الخلية الواحدة من معهد إيرلهام) وEISTA (التحليل النسخي المكاني من معهد إيرلهام)، واللذين بُنيا على إطار عمل nf-core Nextflow ويعتمدان على منظومة Scanpy/AnnData. تم تصميم كلا المسارين ليكونا نمطيين، وقابلين للتوسع، وقابلين للتكرار، مع دعم التنفيذ على الأجهزة المحلية، وعناقيد الحوسبة عالية الأداء (HPC)، والمنصات السحابية.
بنية مسار العمل:
EISCA: يدعم تحليل scRNA-seq القائم على القطرات (مثل 10x Chromium) والقائم على الأطباق (مثل Smart-seq2).
هيكل سير العمل: ينقسم كلا المسارين إلى ثلاث مراحل:
التحليل الأولي: يحول البيانات الخام (ملفات FASTQ أو ملفات التصوير) إلى مصفوفات عدد الجينات والخلايا. يدعم EISCA أدوات Kallisto/Bustools وSalmon Alevin وSTARsolo للمحاذاة والتقدير الكمي. أما EISTA فيستخدم أداة Vizgen Postprocessing Tool (VPT) لتقسيم الخلايا، وتوزيع النسخ، وحساب البيانات الوصفية، مع إعدادات قابلة للتخصيص لـ Cellpose.
التحليل الثانوي: يركز على ضبط الجودة (QC)، والترشيح، والتطبيع، والتجميع العنقودي (clustering)، والدمج. يستخدم هذا الطور أدوات Scanpy وSquidpy (للبيانات المكانية). وتشمل الميزات الرئيسية الترشيح الآلي للخلايا (الحدود الصارمة، وإزالة القيم المتطرفة بناءً على المدى الربيعي IQR)، وكشف القطرات المزدوجة (Scrublet)، وإزالة القطرات الفارغة (CellBender)، وتصحيح الدفعات (BBKNN أو Harmony أو Scanorama أو scvi-tools). يتميز EISTA بدمج التحليلات الإحصائية المكانية، بما في ذلك مقاييس المركزية، وإثراء الجوار، والارتباط الذاتي المكاني (Moran's I).
التحليل الثالثي: يغطي المهام اللاحقة مثل توصيف أنواع الخلايا (عبر CellTypist أو scvi-tools/scANVI)، وتحليل التعبير التفاضلي (Scanpy أو PyDESeq2 أو scvi-tools)، واستنتاج التواصل بين الخلايا (CellChat أو LIANA+ في نسخة Galaxy).
مبادئ التصميم:
المعيارية مقابل المرونة: يمكن تنفيذ مسارات العمل بشكل كامل للحصول على نتائج أولية سريعة، أو تشغيلها بشكل نمطي للسماوة بالتحليل الاستكشافي المتكرر والمعدل بالمعلمات.
تنفيذ Galaxy: تتوفر نسخة من مسار عمل EISTA عبر نظام Galaxy، باستخدام كائنات SpatialData ودمج LIANA+ لتحليل ارتباط الليغاند والمستقبل.
النتائج الرئيسية تحقق المؤلفون من صحة المسارات باستخدام دراستي حالة متميزتين:
دراسة حالة التحليل النسخي المكاني (EISTA):
مجموعة البيانات: بيانات Vizgen MERFISH لأوراق نبات Arabidopsis thaliana المصابة ببكتيريا Pseudomonas syringae.
النتائج: نجح المسار في معالجة بيانات التصوير الخام، وترشيح الخلايا منخفضة الجودة، وتحديد 11 عنقودًا خلويًا رئيسيًا. كشف التحليل الإحصائي المكاني أن عناقيد محددة (مثل Cluster 0) أظهرت مركزية درجة عالية وأصبحت مركزة مكانيًا بمرور الوقت خلال العدوى (من 4 إلى 24 ساعة).
الرؤية البيولوجية: حدد تحليل التعبير التفاضلي جينات مرتبطة بالمناعة (مثل AT1G19250/FMO1، وAT2G13810/ALD1، وAT1G62300/WRKY6) كجينات متغيرة مكانيًا ومرتفعة التعبير في الخلايا النشطة مناعيًا. نجح نقل الملصقات (Label transfer) باستخدام scANVI في رسم الهويات النسيجية (البشرة، الميزوفيل، الأوعية) والحالات المناعية، مما أعاد إنتاج التوسع الزمني لخلايا الميزوفيل النشطة مناعيًا كما لوحظ في الدراسة الأصلية.
دراسة حالة الخلية الواحدة (EISCA):
مجموعة البيانات: بيانات scRNA-seq من مرضى الإنتان (Sepsis) البشر (الخلايا أحادية النواة في الدم المحيطي - PBMCs)، مع مقارنة الأصحاء، والناجين، وغير الناجين (في المراحل المبكرة والمتأخرة).
النتائج: أعاد EISCA إنتاج النتائج البيولوجية الرئيسية من الدراسة الأصلية، بما في ذلك النقص التدريجي في الخلايا البائية (B cells) وإثراء الصفائح الدموية لدى غير الناجين.
الرؤية البيولوجية: سلط تحليل التعبير التفاضلي الضوء على ارتفاع مستويات الجينات المرتبطة بالالتهاب والجينات المحفزة بالإنترفيرون في الخلايا الوحيدة (monocytes) لدى غير الناجين. كما رصد المسار ظهور خلايا سلفية حمراء دموية (يُستدل عليها بجينات الهيموجلوبين) في حالات الإنتان الشديد. كشف تحليل التواصل بين الخلايا عن تعزيز الإشارات المسببة للالتهاب لدى غير الناجين.
الأهمية والادعاءات يزعم البحث أن EISCA وEISTA يوفران حلولًا فعالة ومرنة وقابلة للتوسع للتحليلات الشاملة للخلية الواحدة والتحليل النسخي المكاني. تكمن أهميتهما الأساسية في:
سهولة الوصول: خفض حواجز الدخول للمستخدمين غير المتخصصين من خلال إزالة الحاجة إلى التثبيت اليدوي للأدوات وتوفير تقارير فورية وقابلة للتفسير.
الشمولية: تقديم إطار عمل موحد يمتد من التحليلات الأولية والثانوية إلى التحليلات الثالثية، وهو أمر ملحوظ بشكل خاص في ESTA كونه أحد مسارات العمل القليلة التي توفر تحليلاً شاملاً للتحليل النسخي المكاني القائم على التصوير.
القابلية للتكرار: الاستفادة من إطار عمل nf-core Nextflow لضمان تكرار مسارات العمل عبر بيئات حوسبية مختلفة.
الفائدة البيولوجية: إثبات القدرة على استعادة رؤى بيولوجية ذات مغزى في سياقات بيولوجية معقدة (مناعة النبات والإنتان البشري) دون الحاجة إلى كتابة أكواد برمجية مخصصة لتدفقات العمل القياسية.
يؤكد المؤلفون أن هذه المسارات توازن بين المعيارية والمرونة المطلوبة للطبيعة الاستكشافية لأبحاث الخلية الواحدة، مما يسمح بالتقييم الأولي السريع والتحقيق العميق القائم على الفرضيات.