NanoPlasmiQC: Full plasmid sequencing with ONT long-reads and automatic data analysis
تقدم الورقة البحثية NanoPlasmiQC، وهو سير عمل مؤتمت ومنخفض التكلفة يستخدم تقنية تسلسل القراءة الطويلة من Oxford Nanopore لتسلسل وتحليل البلازميدات الكاملة في غضون يوم واحد، مما يوفر بديلاً متفوقاً لطريقة تسلسل سانجر التقليدية.
المؤلفون الأصليون:de Oliveira, J. A. V. S., Ng, V., Wolff, K., Pucker, B.
تخيل أنك طاهٍ يحاول خبز كعكة معقدة للغاية ومحددة جداً. في الماضي، لكي تتأكد من ضبط الوصفة بشكل صحيح، كان عليك تذوق فتات صغير جداً من الأعلى، ثم قطعة أخرى من المنت المنتصف، وأخرى من الأسفل. كان عليك القيام بذلك عشرات المرات، وكتابة ملاحظات لكل لقمة، آملًا ألا تفوتك أي مكونات غريبة مخبأة في المنتصف. كان هذا يشبه تسلسل سانجر (Sanger sequencing): دقيق، ولكنه بطيء ومكلف، ولم يسمح لك إلا بتذوق قطع صغيرة من الكعكة في كل مرة.
الآن، تخيل أن لديك ماصة (شفاطة) سحرية وطويلة جداً يمكنها شرب الكعكة بأكملها في جُرعة واحدة. ستحصل على الوصفة كاملة، من الكريمة حتى الطبقة السفلية، في خطوة واحدة. هذا هو ما يفعله NanoPlasmiQC للعلماء.
إليك قصة هذه الأداة الجديدة، مشروحة ببساطة:
المشكلة: خلطة "بطاقة الوصفة"
في علم الأحياء، يستخدم العلماء حلقات صغيرة من الحمض النووي تسمى البلازميدات (plasmids) كبطاقات وصفة لبناء سمات جديدة في النباتات أو لإنتاج الأدوية. وأحياناً، عندما يقومون بنسخ هذه البطاقات، يحدث خطأ مطبعي (طفرة).
الطريقة القديمة: للتحقق من الأخطاء المطبعية، كانوا يستخدمون طريقة "الفتات الصغير" (تسلسل سانجر). إذا كان البلازميد كبيراً أو يحتوي على أنماط متكررة صعبة، فقد يغفلون عن الخطأ. كما كانت الطريقة مكلفة؛ حيث يمكن أن يكلف فحص بلازميد واحد ما يعادل سعر عشاء فاخر.
الطريقة الجديدة: لدى العلماء تقنية جديدة تسمى أكسفورد نانوبور (Oxford Nanopore - ONT). إنها تشبه تلك الماصة الطويلة جداً. يمكنها قراءة البلازميد بأكرا في قراءة واحدة مستمرة.
الحل: "بار العصائر" الشامل
قام مؤلفو هذه الورقة البحثية، وهم فريق من جامعة بون، بابتكار سير عمل يسمى NanoPlasmiQC. فكر في الأمر كأنه "بار عصائر" (Smoothie Bar) عالي التقنية.
المكونات (تحضير العينة): بدلاً من صنع عصير واحد في كل مرة، يأخذون قطرة صغيرة من العديد من البلازميدات المختلفة ("المكونات") ويخلطونها جميعاً في دلو واحد كبير.
الآلة (التسلسل): يصبون هذا الدلو في آلة تقرأ الحمض النووي. ولأن الآلة حساسة للغاية، يمكنها قراءة الوصفة الكاملة لكل بلازميد في الخليط، على الرغم من أنها جميعاً تسبح معاً.
الفلتر السحري (البرنامج): هذا هو الجزء العبقري حقاً. بمجرد أن تخرج الآلة البيانات، يعمل برنامج كمبيوتر مخصص (مكتوب بلغة بايثون) كأمين مكتبة ذكي للغاية.
يأخذ كومة ضخمة من قراءات الحمض النووي المختلطة.
يقوم بفرزها، واضعاً القراءات الخاصة بـ "البلازميد أ" في كومة، و"البلازميد ب" في كومة أخرى.
يقارن القراءات الجديدة بـ "الوصفة المتوقعة" (ما ظن العالم أنه صنعه).
يسلط الضوء على أي أخطاء مطبعية، أو مكونات مفقودة، أو طفرات غريبة.
لماذا يعد هذا أمراً مهماً؟
السرعة: في الماضي، قد يستغرق فحص بلازميد واحد أياماً أو أسابيع. مع هذه الطريقة الجديدة، يمكنك خلط وتسلسل وتحليل المجموعة بأكملها في يوم واحد فقط.
التكلفة: نظرًا لإمكانية خلط عشرات البلازميدات في عملية واحدة، تنخفض تكلفة البلازميد الواحد بشكل كبير. يمكن أن يكون أرخص من طلب اختبار "فتات" واحد من مختبر تجاري.
الدقة: لأنهم يقرؤون الشيء بأكمله دفعة واحدة، فإنهم لا يفوتون الأجزاء الصعبة حيث ينطوي الحمض النووي على نفسه (التكرارات المعكوسة)، والتي غالباً ما تخفي الطفرات.
"إثبات المفهوم"
لإظهار نجاح الأمر، اختبروه على بلازميد معروف يسمى pBF3038.
مرروه عبر "بار العصائر" الخاص بهم.
قام الكمبيوتر بفرز البيانات ووجد أن البلازميد صحيح في الغالب، لكنه رصد أيضاً بعض الأخطاء المطبعية الصغيرة (الطفرات) التي لم يلاحظها مبتكروه الأصليون.
قاموا حتى بإنشاء خريطة مرئية (مثل المخطط الهندسي) توضح بالضبط أين تصطف خيوط الحمض النووي، مما أثبت أن الهيكل بأكمله تم تجميعه بشكل صحيح.
الخلاصة
NanoPlasmiQC يشبه الترقية من استخدام عدسة مكبرة إلى ماسح ضوئي عالي السرعة. فهو يسمح للعلماء بالتحقق من "بطاقات وصفة" الحمض النووي الخاصة بهم بشكل أسرع، وأرخص، وأكثر شمولاً من أي وقت مضى. وهذا يعني أنهم يمكنهم قضاء وقت أقل في القلق بشأن الأخطاء المطبعية في بياناتهم ووقت أكثر في ابتكار نباتات وأدوية جديدة بالفعل.
باختختصار: إنه نظام مجاني، تلقائي، يعمل في يوم واحد، يتيح لك قراءة دليل تعليمات الحمض النووي الكامل للبلازميد بالكامل في جُرعة واحدة، حتى لو خلطت مائة منها معاً في دلو واحد.
إليك ملخص تقني مفصل لورقة البحث: "NanoPlasmiQC: تسلسل البلازميد الكامل باستخدام قراءات ONT الطويلة والتحليل التلقائي للبيانات."
1. بيان المشكلة
محدودية تسلسل سانجر (Sanger Sequencing): تعتمد عملية التحقق من صحة البلازميدات تقليديًا على تسلسل سانجر، الذي يمتلك حدًا لطول القراءة يبلغ حوالي 1 كيلو قاعدة (1 kb). وهذا يستليزم إجراء تفاعلات وتصميم بادئات (primers) متعددة لتغطية البلازميدات الكاملة (التي تتراوح عادةً بين 5 إلى 20 كيلو قاعدة)، مما يجعل العملية كثيفة العمالة وعرضة لفقدان الطفرات في المناطق المعقدة مثل التكرارات المتناوبة المقلوبة (والتي قد تحتوي على ما يصل إلى 40% من الطفرات غير المتوقعة).
التكلفة والقابلية للتوسع: غالبًا ما تتجاوز خدمات تسلسل البلازميد الكامل التجارية 10 يورو لكل بلازميد. وفي المشاريع التي تتضمن إنتاج والتحقق من عدد كبير من البلازميدات، تتراكم هذه التكاليف بشكل كبير.
المحدودية التاريخية للقراءات الطويلة: بينما توفر تقنية Oxford Nanopore Technologies (ONT) قراءات طويلة قادرة على تغطية البلازميدات بأكملها، إلا أن الإصدارات السابقة كانت غير محبذة للتحقق من البلازميدات بسبب معدلات الخطأ العالية (الإدراج والحذف - indels) في القراءات الخام.
الحاجة إلى الأتمتة: تتطلب سير العمل الحالية لتحليل البلازميدات غالبًا تدخلًا يدويًا أو تفتقر إلى مسار عمل مؤتمت ومبسط مصمم خصيصًا للتحقق من البلازميدات المصنعة مقابل التسلسلات المتوقعة.
2. المنهجية
طور المؤلفون NanoPlasmiQC، وهو سير عمل فعال من حيث التكلفة وشامل (من البداية إلى النهاية) يجمع بين التحضير المخبري للعينات ومسار تحليل معلوماتي حيوي يعتمد على لغة بايثون (Python) ومؤتمت بالكامل.
أ. بروتوكول المختبر الرطب (تحضير العينة والتسلسل)
استراتيجية التجميع (Pooling Strategy): يتم تجميع البلازميدات عن طريق أخذ 1 ميكرولتر من كل عينة (بغض النظر عن الاختلافات في الحجم) لتبسيط التعامل. يتم تخفيف الخليط إلى 10-15 نانوغرام/ميكرولتر، مع استخدام إجمالي 150 نانوغرام لتحضير المكتبة.
تحضير المكتبة: يستخدم مجموعة أدوات التسلسل السريع (SQK-RAD114) من تقنية ONT.
منصة التسلسل: تم الإجراء على خلايا تدفق PromethION (كيمياء R10). والأهم من ذلك، يستخدم المؤلفون خلايا تدفق معاد استخدامها (تم غسلها عبر خطوة تعتمد على DNase، وهي EXP-WSH004) بعد مشاريع جينوم النبات، مما يقلل التكاليف بشكل كبير.
استدعاء القواعد (Basecalling): تم باستخدام Dorado v1.4.0 في وضع الدقة العالية (HAC).
ب. مسار التحليل المعلوماتي الحيوي (NanoPlasmiQC)
يتم تغليف التحليل في سكربت بايثون متاح على GitHub، والذي يؤتمت الخطوات التالية:
الرسم الخرائطي (Mapping): يتم محاذاة القراءات مع التسلسلات المتوقعة باستخدام Minimap2 (الإصدار 2.26) مع تفعيل خيار secondary=no لضمان المحاذاة المحددة.
تقسيم البيانات: يتم تحويل ملفات SAM إلى BAM، وفرزها وتقسيمها حسب مرجع البلازميد باستخدام Samtools (الإصدار 1.19) للسماح بالتحليل المستقل لكل بلازميد.
تصفية القراءات: يتم استخراج القراءات وفحص التكرارات باستخدام Seqkit (الإصدار 2.3.0).
أخذ العينات الفرعية (Subsampling): نظرًا للتغطية العالية، يتم أخذ عينات فرعية للوصوا إلى عمق مستهدف (بهدف تحقيق >100x) لتحسين كفاءة التجميع.
استدعاء المتغيرات (Variant Calling): يتم باستخدام Bcftools (الإصدار 1.19) مع فلاتر جودة محددة (-Q 10, -q 10, QUAL>20, DP>coverage cutoff).
التجميع من الصفر (De Novo Assembly): يتم تجميع القراءات الفرعية لكل بلازميد باستخدام Miniasm (الإصدار 0.3).
الصقل (Polishing): يتم صقل التسلسلات المجمعة باستخدام Racon (الإصدار 1.5.0).
التصور (Visualization): يتم تصور النتائج عبر برنامج Integrative Genomics Viewer (IGV) للفحص اليدوي.
3. المساهمات الرئيسية
خفض التكلفة: يثبت سير العمل أن تسلسل البلازميد الكامل يمكن أن يكون أرخص من تفاعل سانجر واحد عند تجميع عدة بلازميدات، خاصة باستخدام خلايا التدفق المعاد استخدامها.
الأتمتة: يوفر مسار عمل مؤتمت بالكامل وسهل الاستخدام بلغة بايثون، يقوم بمعالجة البيانات الخام إلى تسلسلات بلازميد تم التحقق منها في غض_يوم واحد، مما يلغي الحاجة إلى تصميم بادئات يدوي أو عملية "المشي بالبادئ" (primer walking) المعقدة.
التحقق من دقة القراءات الطويلة: يثبت أن تقنيات ONT الحديثة (R10) مع استدعاء القواعد عالي الدقة، أصبحت الآن دقيقة بما يكفي للتحقق من البلازميد الكامل، متجاوزة المخاوف التاريخية المتعلقة بأخطاء الإدراج والحذف (indels).
المصدر المفتوح: مسار العمل بالكامل مفتوح المصدر ومتاح عبر GitHub (bpucker/NanoPlasmiQC).
4. النتائج
مخرجات التسلسل: أسفرت تجربتان على خلايا تدفق PromethION معاد استخدامها عن 12.75 جيجابايت و0.58 جيجابايت من البيانات، على التوالي.
خصائص القراءات: بلغت قيم N50 لأطوال القراءات 7.2 كيلو قاعدة و2.6 كيلو قاعدة، على التوالي. ورغم أنها أقل من تشغيلات جينوم النبات النموذجية، إلا أن هذه الأطوال كافية لتغطية معظم البلازميدات (5-20 كيلو قاعدة) في قراءات فردية.
التغطية: حقق المؤلفون تغطية كافية (باستهداف >100x) للتحقق من التسلسلات رغم انحياز التجميع.
إثبات المفهوم: تم اختبار سير العمل على البلازميد pBF303、38.
الرسم الخرائطي: أظهر تغطية كاملة لطول البلازميد في قراءات فردية.
كشف المتغيرات: نجح في تحديد الطفرات النقطية في النسخة المتسلسلة مقارنة بالمرجع.
التجميع: نجح التجميع من الصفر (de novo) في إعادة بناء خريطة البلازميد، والتحقق من جميع العناصر الجينية الرئيسية.
5. الأهمية
الكفاءة: يتيح لعلماء الحياة التحقق من بنى البلازميد الكاملة (بما في ذلك التكرارات المعقدة والمدرجات الكبيرة) في يوم واحد دون الحاجة لتصميم بادئات تسلسل.
الأثر الاقتصادي: يخفض بشكل كبير حاجز التحقق عالي الإنتاجية من البلازميدات في مشاريع التكنولوجيا الحيوية والبيولوجيا التخليقية، مما يجعلها متاحة للمختبرات ذات الميزانيات المحدودة.
الموثوقية: يثبت أن القراءات الطويلة لـ ONT، عند دمجها مع مسار NanoPlasmiQC، توفر بديلاً قويًا لتسلسل سانجر، وقادرة على اكتشاف الطفرات غير المتوقعة التي قد تغفل عنها طرق القراءات القصيرة أو تتطلب تفاعلات متعددة لإيجادها.
الاستدامة: من خلال استخدام خلايا التدفق المعاد استخدامها، يعزز سير العمل الاستخدام المستدام لموارد التسلسل.