Systematic benchmarking of small variant calling pipelines for long-read RNA sequencing data
تقوم هذه الدراسة بقياس أداء خطوط معالجة تحديد المتغيرات الصغيرة وتحديد الأنماط الفردية (haplotype phasing) بشكل منهجي لتسلسل الحمض النووي الريبوزي (RNA) طويل القراءة عبر مجموعات بيانات وتقنيات متنوعة، كاشفةً أن جودة التسلسل هي المحدد الرئيسي للأداء مع تحديد Clair3-RNA وDeepVariant وlongcallR كأفضل أدوات تحديد المتغيرات، وWhatsHap أو HapCUT2 كأدوات مثالية لتحديد الأنماط الفردية اعتماداً على السياق المحدد.
يعد التباين الجيني هو مصدر التنوع البشري، وهي الاختلافات الدقيقة في حمضنا النووي التي تشكل كل شيء، من لون العين إلى القابلية للإصابة بالأمراض. لعقود من الزمن، اعتمد العلماء على تسلسل القراءات القصيرة للعثور على هذه الاختلافات، وهي طريقة تقسم الجينوم إلى قطع صغيرة، وتقرأها، ثم تحاول إعادة تجميع الصورة كأنها لغز. ومع ذلك، فإن هذا النهج غالبًا ما يفقد سياق كيفية ترابط هذه القطع مع بعضها البعض، خاصة عند النظر إلى الحمض النووي الريبوزي (RNA)، وهو الجزيء الذي يحمل التعليمات من الحمض النووي (DNA) لبناء البروتينات. وتعمل تقنية جديدة، وهي تسلسل القراءات الطويلة، على حل هذه المشكلة من خلال قراءة خيوط كاملة من الحمض النووي الريبوزي دفعة واحدة، مما يحافظ على القصة الكاملة لكيفية تعبير الجينات. ولكن بينما توفر هذه التقنية رؤية أوضح، فإنها تجلب أيضًا تحديًا جديدًا: الأدوات المستخدمة للعثور على الاختلافات الجينية بُنيت للبيانات القديمة المجزأة. كان على العلماء معرفة ما إذا كانت هذه الأدوات يمكن أن تتعامل مع الخيوط الطويلة الجديدة دون ارتكاب أخطاء، وأي طرق محددة ستعمل بشكل أفضل لهذه المهمة.
وللإجابة على ذلك، أجرى باحثون في جامعة زيورخ اختبارًا واسع النطاق لمعرفة مدى قدرة البرامج الحاسوبية المختلفة على إيجاد التغيرات الجينية الصغيرة في بيانات الحمض النووي الريبوزي ذات القراءات الطويلة. لقد عاملوا المشكلة كفحص دقيق لمراقبة الجودة، حيث قاموا بتشغيل خمسة برامج متخصصة مختلفة إلى جانب أداة قياسية تُستخدم للبيانات الأقدم. وقد غذوا هذه البرامج ببيانات من ثلاثة خطوط خلوية بشرية معروفة تم تسلسلها باستخدام طرق متنوعة من شركتين رائدتين، وهي "أكسفورد نانوبور" (Oxford Nanopore) و"باسيفيك بايوساينسز" (Pacific Biosciences). لم يكن الهدف مجرد معرفة أي برنامج يجد أكبر عدد من التغييرات، بل فهم أي مجموعات من طرق التسلسل، والبرامج الحاسوبية، وخطوات معالجة البيانات تنتج النتائج الأكثر دقة. كما اختبر الباحثون مدى قدرة هذه البرامج على تجميع التغيرات الجينية معًا في "أنماط فردانية" (haplotypes)، وهي مجموعات من الاختلافات التي تُورث معًا على كروموسوم واحد، وهي خطوة حاسمة لفهم كيفية تأثير تركيبات جينية محددة على الصحة.
وكشفت الدراسة أن العامل الأكثر أهمًا في الحصول على نتائج دقيقة لم يكن البرنامج الحاسوبي نفسه، بل جودة ونوع البيانات القادمة من جهاز التسلسل. ووجد الباحثون أن الكيمياء المحددة المستخدمة في تحضير عينات الحمض النووي الريبوزي كان لها تأثير هائل على النتيجة. فقد أظهرت البيانات الناتلة باستخدام طرق "باسيفيك بايوساينسز" تفوقًا مستمرًا على بيانات "أكسفورد نانوبور"، ويرجع ذلك أساسًا إلى أن بيانات "باسيفيك بايوساينسز" احتوت على أخطاء أقل وسمحت برؤية أكثر اكتمالاً للشفرة الجينية. ومن بين البرامج الحاسوبية المختبرة، برز برنامج يسمى "Clair3-RNA" كالأكثر تعددًا في الاستخدامات، حيث أدى عمله بكفاءة عبر جميع أنواع البيانات المختلفة ووجد كلاً من التغيرات أحادية الحرف والإدراجات أو الحذوفات الصغيرة بدقة عالية. كما قدم برنامج آخر، وهو "DeepVariant"، أداءً جيدًا للغاية، لا سيما عند تحليل بيانات "باسيفيك بايوساينسز" عالية الجودة. ومع ذلك، أظهرت الدراسة أنه لا توجد أداة واحدة مثالية لكل موقف؛ إذ يعتمد الخيار الأفضل بشكل كبير على طريقة التسلسل المحددة المستخدمة.
كما بحث الباحثون فيما إذا كانت خطوة شائعة في معالجة البيانات، والتي تتضمن إعادة تشكيل قراءات الحمض النووي الريبوزي الطويلة لتبدو مثل قراءات الحمض النووي (DNA) القصيرة التي تتوقعها الأدوات القديمة، لا تزال ضرورية. ووجدوا أن هذه الخطوة في إعادة التشكيل كانت غير ضرورية غالبًا للبرامج الجديدة المصممة خصيصًا لقراءات الحمض النووي الريبوزي الطويلة، بل ويمكن أن تضر بالأداء عن طريق تفكيك المعلومات الهامة. وفي الواقع، سمح الإبقاء على البيانات في تنسيقها الأصلي ذي القراءة الطويلة للبرامج المتخصصة بالعمل بشكل أفضل. كما نظرت الدراسة في مدى قدرة هذه الأدوات على التمييز بين التغيرات الجينية الحقيقية وبين "تعديل الحمض النووي الريبوزي" (RNA editing)، وهي عملية طبيعية يقوم فيها الخلية بتغيير الحمض النووي الريبوزي كيميائيًا بعد صنعه. واكتشفوا أن بعض البرامج كانت أفضل في تجاهل هذه التعديلات الطبيعية، بينما أخطأت برامج أخرى واعتبرتها أخطاءً جينية، ولكن يمكن إدارة هذه المشكلة باستخدام قواعد البيانات الموجودة لتصفية مواقع التعديل المعروفة.
أخيرًا، اختبر الفريق مدى قدرة البرامج على ربط التغيرات الجينية معًا في أنماط فردانية. ووجدوا أن الأدوات المختلفة تقدم نقاط قوة مختلفة: فبعضها كان أفضل في ربط عدد كبير من التغييرات معًا، بينما كان بعضها الآخر أكثر دقة وحذرًا ولكنه يربط تغييرات أقل. وقدم برنامج "longcallR" ميزة فريدة، حيث وجد التغيرات الجينية وربطها معًا في خطوة واحدة، مما وفر بديلًا قويًا للباحثين الذين يحتاجون إلى كلتا النتيجتين في آن واحد. وقد أكد الباحثون هذه النتائج باستخدام بيانات من خطوط خلايا سرطانية، مما أظهر أن الدروس المستفادة من الخطوط الخلوية القياسية تظل قائمة في السياقات البيولوجية الأكثر تعقيدًا. وفي نهاية المطاف، توفر هذه الدراسة خارطة طريق واضحة للعلماء، حيث توضح أنه بينما يهم اختيار البرمجيات، فإن جودة بيانات التسلسل هي أساس النجاح، وأن النهج الأفضل يعتمد على مطابقة الأدوات المناسبة مع النوع المحدد من البيانات التي يتم تحليلها.
ملخص تقني: التقييم المرجعي المنهجي لخطوات معالجة استدعاء المتغيرات الصغيرة لبيانات تسلسل الحمض النووي الريبوزي (RNA) طويل القراءة
بيان المشكلة يتيح تسلسل الحمض النووي الريبوزي طويل القراءة (lrRNA-seq) الكشف عن المتغيرات على مستوى النسخ (transcripts)، مما يربط التغيرات الجينية بنماذج نسخ (isoforms) محددة ويسهل تحليل الأحداث المرتبطة بالربط البديل (splicing) والتمايز الأليلي. ومع ذلك، لا تزال التقييمات المنهجية والمحايدة لخطوات معالجة استدعاء المتغيرات الصغيرة (SNVs وindels) لبيانات lrRNA-seq محدودة. كما أن أداء الأدوات الحالية عبر تقنيات التسلسل المختلفة (تقنيات Oxford Nanopore Technologies [ONT] وPacific Biosciences [PacBio])، واستراتيجيات المحاذاة، وخيارات استدعاء المتغيرات، والسياقات الجينومية، وتحديد الأنماط الفردانية (haplotype phasing) اللاحقة، ليس مفهوماً بشكل كامل. علاوة على ذلك، ليس من الواضح ما إذا كان تحويل ملف المحاذاة (تحويل محاذاة RNA المقطعة إلى تمثيلات متصلة تشبه DNA) لا يزال ضرورياً لمستدعي المتغيرات المصمم خصيصاً لـ lrRNA-seq.
المنهجية أجرى المؤلفون تقييماً مرجعياً منهجياً باستخدام مجموعات بيانات Genome in a Bottle (GIAB) التي تضم ثلاثة خطوط خلوية (HG002، HG004، HG005) تم تسلسلها باستخدام أربعة بروتوكولات لتحضير مكتبات ONT (وهي: cDNA R9.4.1، وcDNA R10.4.1، وdRNA002، وdRNA004) وبروتوكولين من PacBio (وهما: Iso-Seq وMas-Seq/Kinnex). كما وسعت الدراسة التقييم ليشمل مجموعات بيانات LongBench الحديثة (H211 وH526 لخطوط خلايا سرطان الرئة).
تضمنت سير عمل التقييم المرجعي ما يلي:
المحاذاة (Alignment): تمت محاذاة القراءات باستخدام minimap2 وpbmm2 (وهو غلاف محسن لـ PacBio). درست الدراسة تأثير تحويل المحاذاة (باستخدام SplitNCAReads لمحاكاة محاذات DNA المتصلة) مقابل استخدام محاذات RNA المقطعة الخام.
استدعاء المتغيرات (Variant Calling): تم تقييم خمسة من مستدعي المتغيرات المخصصة لـ lrRNA-seq وهي: Clair3-RNA، وDeepVariant (باستخدام نموذج MASSEQ لبيانات PacBio ونموذج ONT للـ DNA الخاص بـ ONT)، وisoLASER، وlongcallR، وlongcallR-nn. كما تم إدراج GATK HaplotypeCaller كمرجع للمقارنة مع تقنيات القراءة القصيرة.
تحديد الأنماط الفردانية (Haplotype Phasing): تم تقييم أداء تحديد الأنماط باستخدام WhatsHap، وLongPhase، وHapCUT2، وHiPhase، بالإضافة إلى مخرجات تحديد الأنماط الأصلية من isoLASER وlongcallR.
مقاييس التقييم: تم قياس الدقة كمياً باستخدام مقاييس الدقة (precision)، والاستدعاء (recall)، ودرجة F1 مقابل استدعاءات الحقيقة الأرضية (ground-truth) ضمن المناطق عالية الثقة، مع تصنيفها حسب عمق التسلسل، ونوع المتغير (SNV/indel)، والسياق الجينومي. كما تم تقييم تحديد الأنماط باستخدام معدل خطأ التبديل (switch error rate)، وتغطية تحديد الأنماط، وN50 لكتل تحديد الأنماط (phase-block N50).
المساهمات الرئيسية والنتائج
خصائص التسلسل كمحدد رئيسي: تحدد الدراسة أن خصائص التسلسل (الكيمياء، ومعدل الخطأ، والإنتاجية، وقابلية التعيين) هي المحدد الأقوى لأداء استدعاء المتغيرات، حيث تفسر قدراً أكبر من التباين في درجات F1 وعدد الإيجابيات الحقيقية مقارنة باختيار مستدعي المتغيرات. تفوقت مجموعات بيانات PacBio (خاصة Mas-Seq) باستمرار على مجموعات بيانات ONT بسبب دقة التسلسل الأعلى والعدد الأكبر من القواعد القابلة للاستدعاء. ومن بين كيمياء ONT، حققت dRNA004 وcDNA R10 أفضل أداء، بينما سجلت dRNA002 أسوأ أداء.
أداء مستدعي المتغيرات:
أظهر Clair3-RNA أوسع نطاق وأكثر قدرة على التطبيق عبر كل من كيمياء ONT وPacBio لكل من SNVs وindels.
قدم DeepVariant أداءً قوياً على بيانات PacBio، خاصة في استدعاء الـ indels، لكنه كان ضعيف الأداء على بيانات ONT بسبب غياب نموذج مخصص لـ ONT RNA-seq.
كان longcallR منافساً في استدعاء SNV على بيانات ONT دون الاعتماد على نماذج مدربة مسبقاً، ووفر مخرجات محددة الأنماط بشكل مباشر.
أظهر isoLASER دقة عالية ولكن استدعاءً (recall) منخفضاً، خاصة في بيانات ONT، وكان متحفظاً في استدعاء المتغيرات.
أظهر longcallR-nn أداءً تنافسياً في مجموعات بيانات محددة، لكنه واجه صعوبات مع dRNA002.
استراتيجيات المحاذاة:
اختيار المحاذي (Aligner): حقق minimap2 بشكل عام درجات F1 أعلى من pbmm2؛ حيث مال pbmm2 إلى إحداث المزيد من الإيجابيات الكاذبة، خاصة في مناطق الربط (splice regions)، ويرجع ذلك على الأرجح إلى الاختلافات في الإعدادات الافتية وغياب مدخلات ترميز نقاط الربط (splice-junction annotation).
تحويل ملف BAM: بالنسبة لمستدعي المتغيرات المخصصة لـ lrRNA-seq، لم يؤدِ تحويل المحاذاة (تقسيم القراءات عند وصلات الإكسون-الإنترون) إلى تحسين الأداء، بل أدى في بعض الحالات إلى تقليل الدقة عبر تعطيل ميزات المحاذاة الخاصة بـ RNA (مثل أنماط CIGAR المتعلقة بالربط) ومعلومات الأنماط الفردانية. ظل التحويل مفيداً فقط لـ DeepVariant في بيانات PacBio للكشف عن الـ indels.
السياق الجينومي وتحرير RNA:
تدهور الأداء في السياقات الجينومية الصعبة، لا سيما في مناطق الهوموبوليمر (خاصة A/T) ومناطق GC المتطرفة. كان استدعاء الـ indels أكثر حساسية لطول الهوموبوليمر من استدعاء الـ SNV.
ساهم تحرير RNA (A-to-I) في حدوث إيجابيات كاذبة. وكانت المستدعيات التي دمجت قواعد بيانات تحرير RNA (على سبيل المثال، تمييز مواقع Clair3-RNA عبر REDIportal) أو التي استخدمت تحديد الأنماط التكراري (longcallR) أفضل في التمييز بين المتغيرات الحقيقية وأحداث التحرير.
تحديد الأنماط الفردانية (Haplotype Phasing): لم تكن هناك طريقة واحدة مثالية لجميع المقاييس. فضل WhatsHap تغطية تحديد الأنماط والاتصال (كتل تحديد أنماط أطول) على حساب ارتفاع معدلات خطأ التبديل، بينما أعطى HapCUT2 الأولوية للدقة مع تغطية أقل قليلاً. وفر longcallR بديلاً قوياً للاستدعاء وتحديد الأنماط بشكل متكامل.
القابلية للتعميم: كانت اتجاهات الأداء الملاحظة في مجموعات بيانات GIAB متسقة مع نتائج مجموعات بيانات LongBench لخلايا سرطان الرئة، مما يعزز هيمنة خصائص التسلسل والتراتب النسبي لأداء المستدعيات.
الأهمية يضع هذا البحث حجر الأساس للتقييم المرجعي لعمليات استدعاء المتغيرات في lrRNA-seq، مؤكداً أنه بينما يهم اختيار الأداة، فإن تقنية التسلسل الأساسية وتحضير المكتبات هما العاملان الأكثر حسمًا للنجاح. ويوفر الدراسة توجيهات محددة لبناء خطوط المعالجة:
يُوصى باستخدام Clair3-RNA كأداة الأكثر تنوعاً لاستدعاء متغيرات lrRNA-seq بشكل عام.
يعد DeepVariant خياراً قوياً لبيانات PacBio، خاصة للـ indels.
يجب تجنب تحويل المحاذاة بشكل عام لمستدعي المتغيرات المخصصة لـ lrRNA-seq للحفاظ على المعلومات طويلة المدى، باستثناء حالة استخدام DeepVariant على بيانات PacBio.
يقدم كل من WhatsHap وHapCUT2 مقايضات متكاملة لعمليات تحديد الأنماط اللاحقة، بينما يعد longcallR مناسباً للاستدعاء وتحديد الأنماط المتكاملين.
تخلص الدراسة إلى أن التحسينات المستقبلية في الحساسية قد تتطلب استراتيجيات إثراء مستهدفة للتغلب على القيود المرتبطة بالتعبير الجيني، وأن هناك حاجة لمزيد من العمل لتقييم تعيين المتغيرات على مستوى النسخ (transcript-level) وتطبيقات الخلية الواحدة (single-cell).