تصف الدراسة الحالية سير عمل لتسلسل البلازميدات النانوية عندما يكون تسلسل الكروموسومات معروفا بالفعل أو غير مهم. تضمن خطوة أولية لإثراء الحمض النووي بالبلازميدات عمقا مثاليا لتسلسل البلازميد، وتعظيم عدد العينات التي يمكن تتبعها لكل خلية تدفق.
مقالة منهجية
* These authors contributed equally
تصف الدراسة الحالية سير عمل لتسلسل البلازميدات النانوية عندما يكون تسلسل الكروموسومات معروفا بالفعل أو غير مهم. تضمن خطوة أولية لإثراء الحمض النووي بالبلازميدات عمقا مثاليا لتسلسل البلازميد، وتعظيم عدد العينات التي يمكن تتبعها لكل خلية تدفق.
تصف هذه الدراسة سير عمل لتسلسل النانوفور مصمما لتعظيم عائد الحمض النووي للبلازميدات ودقة التسلسل. يغطي الباحثون كل خطوة من خطوات العملية، والتي تشمل استخراج الحمض النووي بالبلازميد، وتحضير المكتبة باستخدام الترميز الشريطي السريع وربط المحولات، وتحميل المكتبة على خلية تدفق للتسلسل، واستدعاء القاعدة (في الوقت الحقيقي أو بعد التسلسل)، وتجميع التسلسل باستخدام Autocycler. كما يحدد المؤلفون متطلبات الحاسوب، والتي تعتمد بشكل كبير على ما إذا كان الاستدعاء الأساسي يتم في الوقت الحقيقي. قام المؤلفون بتسلسل البلازميدات المنقاة من عشرة سلالات سريرية باستخدام هذا السير العملي. لإزالة الفروق بين سلالات الخلايا المضيفة كمتغير، قام المؤلفون أيضا بربط هذه السلالات السريرية بسلالة متلقية مشتركة وسلسلوا الترانزاجينات. أنتجت هذه التجارب تجميعات بلازميدات تطابق أحجام البلازميدات المستخلصة من حركة الجل، والتي في ثلثي الحالات تتطابق أيضا مع تجميعات مصقولة قصيرة القراءة أو ضمن 0.05٪ من هوية التسلسلات الزوجية. كما أن إثراء تسلسلات الحمض النووي للبلازميدات يعزز كفاءة تسلسل البلازميدات، مما يسمح بالتسلسل المتوازي لما يصل إلى 24 عينة مع تغطية مثالية. هذا البروتوكول كاف للبلازميدات التي تتراوح سعتها بين 4 إلى 174 كيلوبايت ويتحمل تلوث كروموسوم بنسبة تصل إلى 72٪ دون التأثير على الدقة. هذا البروتوكول مثالي للحالات التي يكون فيها تسلسل الكروموسومات معروفا بالفعل أو لا يثير الاهتمام، مثل التحقق من تسلسل البلازميدات، وتحسين دقة تسلسلات البلازميدات التي تم الحصول عليها بالفعل بواسطة WGS، أو تسلسل البلازميدات التي يتم جمعها بالترفق. الكائن الحي المستخدم لتوضيح هذا العمل هو الإشريكية القولونية، والتي تمثل المضيفات المستخدمة في التعبير الجيني المؤتلف، والتقاط البلازميدات من خلال الاقتران.
البلازميدات هي قطع DNA ذاتية التكاثر توجد في كل مكان في البكتيريا. تتميز بأحجام متغيرة للغاية، حيث تتراوح من <800 وحدة توقف لكل بوصة إلى 2.5 ميغابت لكل بوصة، وعادة ما تتبع توزيعا ثنائي النمط، حيث توجد ذروة أصغر تقابل البلازميدات متعددة النسخ غير المرتبطة، وذروة أكبر تقابل البلازميدات المنخفضة النقوش، وغالبا ما تكون مرافق 1,2. تلعب البلازميدات دورا رئيسيا في علم البيئة الميكروبي، حيث تحمل جينات تساعد مضيفيها على التكيف مع التحديات البيئية أو الوصول إلى بيئات بيئية جديدة. تعمل البلازميدات أيضا كمنصات للتطور المتسارع، حيث توفر زيادة المرونة الجينية وتتوسط في نقل الجينات الأفقي عبر السلالات والأنواع والأجناس وحتى العائلات3،4. أكدت التحليلات الشبكية والمقارنة الجينومية دور المقاوم البيئي كمخزن دائم لجينات مقاومة المضادات الحيوية ذات الصلة سريريا (ARGs)، مما يسهل صيانتها وانتشارها عبر البيئات البيئية5. تتبع البلازميدات وحمولتها من جينات مقاومة المضادات الحيوية والضراوة يساعد في تحديد طرق الانتقال6. لذلك، فإن تحديد وتوصيف البلازميدات يحظى باهتمام كبير لمجموعة متنوعة من التطبيقات الصحية العامة، وعلم الأحياء الدقيقة السريري، والتقنيات الحيوية.
فتح إدخال تقنيات التسلسل من الجيل القادم (NGS)، التي تمكن من تسلسل الحمض النووي المجهول عالي الإنتاجية، الباب أمام التوصيف الجينومي للعزلات الميكروبية. يتيح تسلسل الجينوم الكامل (WGS) تتبع ظهور وانتشار وانتقال مسببات الأمراض الجديدة، والسلالات المقاومة للأدوية، والمتحورات التي تتجنب اللقاحات. كما يستخدم بشكل متزايد للمساعدة في تحديد طرق انتقال العدوى أثناء التفشيات السريرية، مما قد يوجه تدخلات للوقاية المستقبلية من التفشي7. أخيرا، يسمح WGS بتنميط مجموعات البدلاء البدائل، والتي يمكن أن تؤثر على القرارات العلاجية8، رغم أن ربط المجموعات البدينية بأنماط مقاومة المضادات الحيوية النوعية والكمية لا يزال قيد التطوير9.
ومع ذلك، فإن بيانات تسلسل القراءة القصيرة (SR) التي تنتجها WGS عادة على نطاق واسع تكون عادة في نطاق 400–500 بيت. هذا الحجم أقصر من حجم معظم العناصر الجينية المتنقلة (MGEs)، وبالتالي لا يمكن حلها10. يتعقد تحديد تسلسلات البلازميدات في بيانات WGS بشكل سيليكي بسبب التغطية غير الكاملة، وبسبب التسلسلات المشتركة بين البلازميدات وبين الكروموسومات البكتيرية والبلازميدات الأخرى10. وبالتالي، غالبا ما تكون التجميعات مجزأة وغير مكتملة، مما يجعل موقعها تحت الجينوم (الكروموسومي أو البلازميدي) مستحيلا نسبه10.
تنتج تقنية تسلسل القراءة الطويلة (LR) تجميعات متجاورة للغاية تمتد عبر عدد أكبر من التكرارات، مما يسهل بشكل كبير إنتاج تجميعات البلازميدات الكاملة. هناك نوعان أساسيان من تقنيات القراءة الطويلة والجزيء الواحد متاحان11,12. أحد هذه الأساليب، وهو تسلسل النانو-مسام، أصبح شائعا في علم الجينوم الميكروبي بسبب تكلفته المنخفضة، وسرعة الاستجابة، وملاءمته للبيئات ذات الموارد المحدودة13،14.
تستخدم أجهزة تسلسل المسام النانوية خلايا تدفق تحتوي على مصفوفات تضم آلاف المسامات النانوية. يتم تحميل الحمض النووي في خلية التدفق، يفتحها بروتين حركي ذو نشاط هليكاز، ويمر عبر النانوسام مدفوعا بفارق تيار. عندما تمر النيوكليوتيدات في سلسلة واحدة عبر مسام نانوية فردية، يتعطل التيار لإنتاج إشارة مميزة تقاس بواسطة حساس إلكتروني متصل بكل مسام. يتم فك تشفير هذه الإشارة باستخدام خوارزميات الاستدعاء الأساسي، ويتم تجميع البيانات الناتجة وتحليلها لاحقا باستخدام أدوات المعلوماتية الحيوية المختلفة.
مقارنة بأساليب تسلسل SR، يميل تسلسل النانوسام إلى إنتاج جودة قراءة أقل (درجات Q أقل) وتوليد إدخالات/حذف في تسلسلات البوليمرات المتجانسة أثناء تحليل بيانات الإشارات الكهربائية (الاستدعاء الأساسي)15. كما لوحظ أن المثيلة القريبة في الحمض النووي الأصلي يمكن أن تسبب أخطاء منهجية في استدعاء القاعدة16. كان من الحلول الشائعة المعتمدة تضمين بيانات قصيرة القراءة لتصحيح الأخطاء بعد التجميع. ومع ذلك، فإن هذه العملية (المعروفة باسم "التلميع") تزيد من التكلفة والتعقيد التشغيلي17. زيادة تغطية التسلسل يمكن أن تحسن الدقة أيضا، لكنها ليست دائما ممكنة تقنيا.
في عام 2024، أظهر أن خلايا تدفق تسلسل النانوبور من أكسفورد (ONT) R10.4.1 وكيمياء V14 تنتج تجميعات عالية الجودة ومتصلة لبكتيريا G، بينما لا تزال التحليلات الدقيقة على مستويات النيوكليوتيدات الأحادية تستفيد من بيانات تسلسل SR18. ومع ذلك، فإن الفجوة بين دقة تسلسل LR وSR تزداد تدريجيا. تشمل التحسينات الحديثة عملية أخذ عينات بيانات أكثر دقة بسرعة تسلسل أعلى، وإدخال مستدعي قاعدة جديد مصمم للحمض النووي البكتيري الأصلي والمثيلة، مما أدى إلى زيادة دقة النداء الأساسي،وأدوات تحليل محدثة مثل خوارزميات جديدة للمرحلة، واستدعاء المتغيرات، واستدعاء الأليل النجمي. لاحظ أيضا أن أحدث خوارزمية دورادو V.0.9.0 ONT (2025) يمكن استخدامها لاستدعاء إعادة تأسيس بيانات الخطوط المتعرجة الخام. في الواقع، أفادت بعض دراسات WGS بالفعل عن توليد جينومات ميكروبية عن طريق تسلسل النانو المسام بجودة لا يتوقع أن يحسن التلميع القصير القراءة تسلسل الإجماع بشكل كبير20,21.
غالبا ما تحتوي البلازميدات على الجزء الأكبر من جينات المقاومة المكتسبة، خاصة في عزلات المعويةالبكتيرية 22. لذلك، يسمح حل تسلسلات البلازميدات بوضع جينات AR في سياقها، وهو أمر أساسي لتتبع انتقال AMR 23,24. تم بالفعل استخدام تسلسل النانوبور لفحص التطور الدقيق ل BLAKPC الذي يحتوي على البلازميدات25. كما تم استخدام أجهزة MinION في التسلسل اللحظي للبلازميدات واكتشاف جينات المقاومة في العازلات السريريةفي المستشفيات. ومع ذلك، فإن استجواب البلازميدات ليس بعد جزءا من الممارسات الروتينية لمكافحة العدوى.
على الرغم من التحسن الكبير في تسلسل LR مقارنة بتسلسل SR في القدرة على إنقاذ البلازميدات من التسلسل الجيني27، إلا أن بعض المشاكل لا تزال قائمة. تنتج المجمعات طويلة القراءة مثل فلاي28 وكانو29 ورافن30 تجميعات كروموسوماتية عالية الجودة لكنها تواجه صعوبة في تجميع تسلسلات البلازميدات. ويرجع ذلك جزئيا إلى نقص عمق التغطية لأن الحمض النووي للبلازميدات ممثل تمثيلا ناقصا في بيانات تسلسل WGS31. لاحظ أنه بسبب تركيبها المعياري، ووفرة البنى الثانوية، ووجود تكرارات متكررة، فإن التجميع عالي الجودة لتسلسل الحمض النووي البلازميدي يتطلب تغطية أعلى من الحمض النووي الكروموسومي (عمق قراءة 50x32).
إثراء عينات الحمض النووي بتسلسلات البلازميدات هو نهج فعال لزيادة عمق التغطية، لكن عزل أو إثراء حمض البلازميد قبل تسلسل الحمض النووي مكلف جدا أو مرهق جدا للتطبيقات التشخيصية السريرية وقد يسبب تحيزا من خلال إثراء تسلسلات الحمض النووي المختارة. في عام 2017، تم استخدام الحمض النووي الغني بالبلازميدات لمحاولة التجميع الكامل للحمض النووي البلازميدي في العينات السريرية للكشف عنARGs 24. تم تسلسل ثماني عزلات من ستة أنواع مختلفة من البكتيريا المعوية لتعظيم تنوع السكان والبلازميدات الهيكلية. بعد تجربة عدة مجمعات، تبين أن كانو/كانو+بيلون تفوق بوضوح على الآخرين؛ كما حاولت هذه المجموعة البحثية نهج التجميع الفوقي الذي حقق دقة هياكل بلازميدات إضافية، مما أدى إلى تجميع كامل ل 78٪ من جميع البلازميدات المعروفة في العينات المخطئة. ومع ذلك، فإن وجود هياكل التكرار الكبير ومشاكل (بدرجة أقل) في استخراج الحمض النووي لا تزال تشكل مشاكل كبيرة في اكتشاف جميع تجمعات البلازميدات ذات الكونتيج الواحد24.
مؤخرا، تم فحص العينة التكيفية كنهج تكميلي لإثراء تسلسلات البلازميدات منخفضة الوفرة من خلال رفض التسلسلات الكروموسومالية في عينات عزل البكتيريا، رغم أن هذه الطريقة تتطلب تسلسلات مرجعية معروفة33.
هنا، يصف المؤلفون بروتوكولات لتسلسل ONT مخصصة للعينات المغنية بحمض نووي البلازميدي. من خلال الجمع بين زيادة تمثيل الحمض النووي للبلازميدات مع بروتوكولات تحضير المكتبة المحدثة، وتقنية خلايا التدفق، وأدوات تجميع التسلسل، أنتج سير العمل تجميعات تطابق الأحجام المستخلصة من حركة الجل، وفي ثلثي الحالات، إما تطابق تلميع تسلسل SR أو كانت ضمن 0.05٪ من هوية التسلسل الزوجي. كما أن إثراء تسلسلات الحمض النووي للبلازميدات يعزز كفاءة تسلسل البلازميدات، مما يسمح بالتسلسل المتوازي لما يصل إلى 24 عينة مع تغطية مثالية. البروتوكولات ليست مصممة لاستبدال أساليب WGS القياسية للمسامات النانوية، بل لتعظيم دقة وكفاءة تسلسل البلازميدات في الحالات التي يكون فيها تسلسل الكروموسومات معروفا بالفعل أو غير مثير للاهتمام.
يشمل سير العمل استخراج الحمض النووي بالبلازميد، وتحضير المكتبة باستخدام الترميز الشريطي، وربط المحول ببروتين محرك، وفحص خلية التدفق، وتحميل المكتبة على خلية تدفق للتسلسل، وتحليل بيانات الإشارات الكهربائية (في الوقت الحقيقي أو بعد التسلسل)، وتجميع التسلسل باستخدام Autocycler34، وتحليلات لاحقة متنوعة (الشكل 1).

الشكل 1: سير عمل تسلسل البلازميدات باستخدام تقنيات أكسفورد نانوبور. الخطوة 1: استخراج الحمض النووي؛ الخطوة 2: تحضير المكتبة؛ الخطوة 3: التسلسل؛ الخطوة 4: التجميع والتعليق. تم إنشاؤه في BioRender. كورتيس، ج. (2026) https://BioRender.com/u3b15st. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
1. استخراج الحمض النووي من البلازميدات
ملاحظة: استخدم هذا البروتوكول لاستخراج الحمض النووي البلازميدي من سلالة الإشريكية القولونية السريرية والمتجانسات المقابلة لها باستخدام مجموعة استخراج الحمض النووي البلازميديةالتجارية رقم 35. يمكن العثور على بروتوكولات مفصلة للتصريف، فيالقسم 36. البروتوكول مناسب لكل من البلازميدات عالية النصوص والمنخفضة النسخ، وتتراوح من ~3 كيلوبايت إلى ~200 كيلوبايت. تتطلب البلازميدات الكبيرة اهتماما خاصا لأنها غالبا ما يكون لديها عدد نسخ منخفض (ما قد لا يقل عن 1–2 نسخة لكل خلية)، مما يحد من إنتاج الحمض النووي. للتعويض، يوصي البروتوكول بزيادة كتلة الخلية وحجم مخازن التحلل والتعادلية لضمان كمية كافية من حمض نووي البلازميدي لربط الأعمدة. لا يزال بإمكان البلازميدات التي تصل إلى 200 كيلوبايت باستخدام هذا البروتوكول؛ ومع ذلك، فإن قابليتها للتطبيق على البلازميدات الضخمة (>200 كيلوبايت) محدودة.



الشكل 2: طبقة التحليل السطحية المعبدة بعد خطوة طرد مركزي واحدة. يتطلب الأمر الطرد المركزي التسلسلي حتى يصبح السائل الفائق نظيفا لتجنب انسداد عمود الترشيح. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
| أنابيب الفحص القياسية | أنابيب اختبار عينات المستخدم | |
| محلول يعمل (صبغة ومخزن مؤقت، 1:200) | 190 ميكرولتر | 180–199 ميكرولتر |
| القياسي (من المجموعة) | 10 ميكرولتر | – |
| عينة مستخدم | – | 1–20 ميكرولتر |
| الحجم الكلي في كل أنبوب اختبار | 200 ميكرولتر | 200 ميكرولتر |
الجدول 1: أحجام وكميات المواد للفلورومتر. تظهر الكميات لعينة واحدة؛ يمكن إعداد تجمع بناء على عدد العينات التي سيتم تحليلها.
2. الطريقة 2: التحضير للمكتبة
ملاحظة: الاعتبارات. تم إجراء تحضير المكتبة لتسلسل البلازميدات المستخرجة وفقا للتعليمات المقدمة مع مجموعة الترميزالشريطي 40. يستخدم هذا البروتوكول مجموعة تعتمد على نهج التحويل البديل لتجزئة الحمض النووي والرمز الشريطي في نفس الوقت، مما يوفر سير عمل سريع وسهل يقلل من وقت الممارسة ويبسط تحضير المكتبة.
| المادة الكيميائية | الذوبان في درجة حرارة الغرفة | انزلق | الخلط عن طريق التوصيل بالماصة |
| الرموز الشريطية السريعة | ليس مجمدا | ✓ | ✓ |
| المحول السريع | ليس مجمدا | ✓ | ✓ |
| خرز AMPure XP | ✓ | ✓ | اخلط عن طريق الماصة أو الدوامة مباشرة قبل الاستخدام |
| مخزن الترشيح | ✓ | ✓ | ✓ |
| مخزن المحول | ✓ | ✓ | المزج عن طريق الدوامة |
الجدول 2: حجم وكميات المواد المستخدمة في تحضير مكتبة النانوبور. يتم وضع علامة على الإجراءات التي يجب تنفيذها مع كل مادة ب (✓).
| المادة الكيميائية | الحجم/العينة |
| حمض نووي قالب 10 ميكرولتر (200 نانوغرام من الخطوة السابقة) | 10 ميكرولتر |
| 1.5 ميكرولتر من الرموز الشريطية السريعة، واحدة لكل عينة | 1.5 ميكرولتر |
| الحجم الكلي | 11.5 ميكرولتر |
الجدول 3: الكواشف والحجم للترميز الشريطي السريع. يتم عرض الأحجام لكل عينة من الحمض النووي؛ يمكن تحضير خليط مجمع، وإضافة الحمض النووي بشكل منفصل.
| 1-24 باركود | 48 باركود | 72 باركود | 96 باركود | |
| حجم مخزن التفريغ | 15 ميكرولتر | 30 ميكرولتر | 45 ميكرولتر | 60 ميكرولتر |
الجدول 4: حجم مخزن الاستشفاع المستخدم لكل 24 باركود. يتم تحديد الأحجام حسب عدد الرموز الشريطية.
| المادة الكيميائية | الحجم |
| المحول السريع (RA) | 1.5 ميكرولتر |
| مخزن المحول (ADB) | 3.5 ميكرولتر |
| الحجم الكلي | 5 ميكرولتر |
الجدول 5: المواد والكواشف للمحول الكهربائي. يمكن تعديل الأحجام حسب عدد العينات.
3. تشغيل خلايا التدفق
ملاحظة: اختر خلية تدفق مخرجة قياسية عند تسلسل عدد قليل من البلازميدات ذات الاحتياجات المتوسطة للتغطية المتوسطة. توصي نانوبور بإجراء فحص خلية تدفق لتحديد عدد المسام النشطة؛ يجب أن تحتوي خلايا التدفق في أجهزة نانوبور على ما لا يقل عن 800 مسام نشطة كما هو محدد في معايير الضمان. توفر خلية تدفق MinION واحدة مخرجا كافيا لتسلسل عينات DNA متعددة غنية بحمض نووي بلازميد (انظر المزيد من التفاصيل حول التعدد في قسم النقاش)، مع توفير تكلفة أقل ومرونة أكبر من خلية تدفق GRIDion أو Promethean.
| المواد الكيميائية | الحجم لكل خلية تدفق |
| تدفق خلية التدفق (FCF) | 1,170 ميكرولتر |
| مصل البقر الألبومين (BSA) بجرعة 50 ملغ/مل | 5 ميكرولتر |
| حبل خلية التدفق (FCT) | 30 ميكرولتر |
| الحجم الكلي النهائي في الأنبوب | 1,205 ميكرولتر |
الجدول 6: الكواشف والحجم لتجهيز وتحميل خلية التدفق. يتم الإشارة إلى أحجام المواد الكيميائية لكل خلية تدفق.
| المواد الكيميائية | الحجم لكل خلية تدفق |
| مخزن التسلسل (SB) | 37.5 ميكرولتر |
| خرزات المكتبة (LIB) أو حل المكتبة (LIS)، إذا كنت تستخدم | 25.5 ميكرولتر |
| مكتبة الحمض النووي | 12 ميكرولتر |
| المجموع | 75 ميكرولتر |
الجدول 7: المواد الكيميائية والمجلدات لتحميل المكتبة. تعرض المجلدات لمكتبة الحمض النووي؛ يجب خلط خرزات المكتبة (المزايدة) مباشرة قبل الاستخدام.
4. الطريقة 4: التجميع باستخدام جهاز إعادة تدوير
ملاحظة: هناك العديد من الخيارات لتجميع القراءات في كونتيغ. لكل مجمع إيجابياته وتحيزاته. يوصي هذا البروتوكول باستخدام Autocycler، الذي يستخدم العديد من التجميعات المختلفة ويجمعها في تجميع واحد. هذا يتطلب جهدا أكبر، ويمكنك استخدام مجمع واحد لتوفير الوقت. يرجى تثبيت Conda في الطرفية لتحميل البرنامج المناسب. يحتوي مستند الأوامر التكميلية على رابط لملف بيئي لإعداد النظام.
5. الطريقة 5: التعليق
ملاحظة: هناك العديد من أدوات التعليق الفعالة للبكتيريا. البرامج التالية موصى بها للحصول على تعليق شامل.
سلالات الإشريكية القولونية السريرية
لتوضيح استخدام ONT لتسلسل البلازميدات من السلالات السريرية، قام المؤلفون بتسلسل 10 سلالات، وكانت هدية سخية من الدكتور ستيفن ساليبانتي في جامعة واشنطن. كانت هذه السلالات جزءا من مجموعة من 312 عزلة مشتقة من الدم أو البول من الإشريكية القولونية الممرضة في الأمعاء (ExPEC)، تم الحصول عليها خلال فترة الرعاية السريرية الروتينية في مركز جامعة واشنطن الطبي45.
تصريف البلازميدات وتوصيفه بواسطة اختبار حركة الجل
كما أجرى الباحثون اختبار اقترافا باستخدام سلالة الإشريكية القولونية LMB100 كمستلم. التصريف هو عملية تؤدي إلى نقل بلازميد موجود في خلية (المتبرع) إلى خلية أخرى (تسمى المستلم)46. يتطلب النقل اتصالا جسديا بين الخليتين وغالبا ما يتم اكتشافه باستخدام علامات تشير إلى وجود البلازميد الاقترابي في الخلايا المستقبلة. سمح ذلك للمؤلفين بإزالة الفروق في سلالة المضيف كمتغير وتوضيح استخدام البروتوكول لتسلسل البلازميدات التي تم الحصول عليها من خلال التقاط البلازميدات (وهي طريقة تستخدم لتحليل البلازميدات من معقل بيئي معين). جميع اختبارات التصريفات أنتجت متوافقات متقاطعة، مؤكدة وجود بلازميد مرافق في المتبرعين. البروتوكول الذي اتبعه المؤلفون موصوف فيالكتاب 36. حددت المجموعة البلازميدات الموجودة في سلالات المتبرع والمتلقي وقدرت حركتها باستخدام الرحلان الكهربائي الهلامي النبضي (PFGE). تظهر قائمة بعشرة سلالات مانحة ممثلة في الجدول 8، وأحجام البلازميدات المقدرة الموجودة فيها مدرجة في العمود الثالث. البلازميدات خارج نطاق دقة هلام PFGE مستثنية من هذا الجدول، لكن الباحثين تمكنوا من تجميع عدة بلازميدات أصغر حول ~4 كيلو بايت (غير موضح).
| سلالة المتبرعين | البلازميد | الحجم المقدر ل PFGE | طول التجميع | العمق المتوسط | عمق الصفر الكبير | نسبة التغطية | المناصب التي تغيرت | نسبة الدقة قبل SR |
| blood_08 _0081 | pblood_08_ 0081_1 | 170000 | 171829 | 100 | 99.941803 | 54 | 0.0314 | 99.9686 |
| blood_08 _0081 | pBlood_08_ 0081_2 | 82000 | 85576 | 46.9 | 0 | 100 | 2 | 99.9977 |
| blood_d-08 _0094 | pBlood_d-08 _0094 | 179000 | 173334 | 382.6 | 3121 | 98.19943 | 2,347 | 98.646 |
| blood_08 _1447 | pBlood_08 _1447 | 97000 | 99624 | 36.2 | 0 | 100 | 0 | 100 |
| blood_10 _0913 | pBlood_10 _0913_1 | 70000 | 71788 | 49 | 4 | 99.994428 | 188 | 99.7381 |
| blood_10 _0913 | pBlood_10 _0913_2 | 112000 | 118727 | 41.2 | 126 | 99.893874 | 276 | 99.7675 |
| blood_1 1_184 | pBlood_11 _184_1 | 150000 | 146583 | 65.1 | 287 | 99.804206 | 134 | 99.9086 |
| blood_11_184 | pBlood_11 _184_2 | 75000 | 74263 | 82.1 | 0 | 100 | 32 | 99.9569 |
| blood_2011 _0238 | pBlood_ 2011_0238 | 170000 | 167286 | 110.9 | 0 | 100 | 0 | 100 |
| upec_108 | pUPEC_108 | 130500 | 131039 | 103.9 | 0 | 100 | 0 | 100 |
| upec_134 | pUPEC_134 | 160500 | 157829 | 169 | 0 | 100 | 0 | 100 |
| upec_271 | pUPEC_271 | 145500 | 144519 | 86 | 33311 | 76.950436 | 1,557 | 98.9226 |
| upec_90 | pUPEC_90 | 70000 | 74624 | 139.8 | 0 | 100 | 6 | 99.992 |
الجدول 8: تقديرات حجم البلازميد ومقاييس دقة التسلسل للجنة مكونة من 10 سلالات متبرعة. يتم سرد تقديرات حجم البلازميدات بناء على تسلسلات PFGE أو تسلسلات بلازميدات المتبرع ب ONT المجمعة. كما تظهر المقارنات الزوجية بين تجميعات البلازميدات المصقولة بتسلسل ONT فقط وSR كمؤشر على الدقة. تشمل المقاييس: طول التجميع، متوسط العمق، مناطق العمق الصفري، نسبة التغطية، المواقع، ونسبة الدقة قبل SR.
مقاييس تسلسل البلازميدات المعتمدة على ONT
تم ضبط مسار التسلسل على 24 ساعة. تم تعدد 24 عينة إجمالا. البيانات الفنية لكل عينة، بما في ذلك درجة الجودة (Q)، والطول، والتغطية، مدرجة في الجدول 9. أنتج تسلسل النانوبور مخرجا قراءة وجودة كافية لجميع العينات التي تم تحليلها. عبر جميع المكتبات، أنتج التسلسل متوسط 134,115 قراءة لكل عينة، وهو ما يعادل متوسط عائد 576.2 ميجابايت. كان متوسط طول القراءات 8,797 بيتا مع متوسط N50 يبلغ 3,693 قطعة قطع، مما وفر استمرارية قراءة كافية لتجميع البلازميدات. كان متوسط درجة جودة القراءة Q16.0 مع متوسط Q16.6، وتراوحت الجولات الفردية من Q13.7 إلى Q18.0.
| الانفعال المتسلسل | قراءات (كيلوبايت) | القواعد (MB) | متوسط الطول | أقصى طول | الحد الأدنى للطول | N50 | الطول الوسيط | متوسط Q | الوسيط Q | |
| Blood_08_0081 المتبرع | 109.8 | 261 | 5713 | 508418 | 78 | 2377 | 1351 | 14.5 | 15 | |
| Blood_08_0094 المتبرع | 89.3 | 229 | 4177 | 114948 | 81 | 2565 | 1065 | 13.7 | 14 | |
| Blood_08_1447 المتبرع | 87.5 | 246 | 5041 | 135147 | 75 | 2816 | 1609 | 15.1 | 15.7 | |
| Bblood_10_0913 المتبرع | 247.7 | 1026 | 9260 | 174231 | 70 | 4143 | 1661 | 15.5 | 16.1 | |
| Blood_11_154 المتبرع | 51.5 | 110 | 4172 | 91148 | 73 | 2127 | 1387 | 16.1 | 16.5 | |
| Blood_2011_0238 المتبرع | 16.8 | 42 | 6707 | 104848 | 62 | 2534 | 845 | 15.3 | 15.2 | |
| UPEC_108 المتبرع | 50.6 | 77 | 1915 | 90234 | 67 | 1525 | 1103 | 15.8 | 16.2 | |
| UPEC_134 المتبرع | 50 | 143 | 7110 | 158393 | 72 | 2865 | 1133 | 16 | 16.2 | |
| UPEC_217 المتبرع | 64 | 236 | 7512 | 144275 | 78 | 3678 | 1895 | 17.4 | 17.9 | |
| UPEC_90 المتبرع | 52.9 | 157 | 5666 | 95944 | 72 | 2969 | 1593 | 16.8 | 17.2 | |
| Blood_08_0081 LMB100 عبر الاقتران | 53.5 | 119 | 5162 | 225856 | 1 | 2232 | 1009 | 16 | 16.8 | |
| Blood_08_0094 LMB100 عبر التزامن | 160 | 442 | 7257 | 252176 | 1 | 2761 | 1083 | 17.3 | 18.2 | |
| Blood_08_1447 LMB100 عبر الاقتراجات | 238 | 1029 | 10249 | 209221 | 1 | 4320 | 1857 | 17.7 | 18.7 | |
| Blood_10_0913 LMB100 عبر الاقتران | 792.3 | 4341 | 11508 | 474325 | 1 | 5479 | 2731 | 17.7 | 18.8 | |
| Blood_11_184 LMB100 عبر الاقتران | 42.5 | 281 | 17931 | 146626 | 87 | 6615 | 2448 | 16 | 16.7 | |
| Blood_2011_0238 LMB100 عبر الاقتراجات | 147.7 | 585 | 9000 | 344203 | 1 | 3962 | 1781 | 18 | 18.9 | |
| UPEC_108 LMB100 ترانزكوجانت | 173.4 | 778 | 10127 | 168756 | 80 | 4485 | 2079 | 15.3 | 16.1 | |
| UPEC_134 LMB100 عبر التزامن | 105.4 | 667 | 17568 | 349230 | 71 | 6323 | 2232 | 15.5 | 16.3 | |
| UPEC_271 LMB100 عبر الاقتران | 74.2 | 336 | 14521 | 342195 | 83 | 4520 | 1306 | 14.6 | 14.8 | |
| UPEC_90 LMB100 عبر الاقتران | 75.2 | 418 | 15348 | 225661 | 83 | 5561 | 2016 | 15.7 | 16.4 | |
الجدول 9: تقرير NanoQ لمقاييس مخرجات التسلسل للجنة مكونة من 10 سلالات متبرعة. مقاييس مخرجات التسلسل مدرجة كما هي. تشمل هذه القراءات، والقواعد، والطول (المتوسط، الوسيط، الأقصى، والحد)، والأطوال وقيم N50، وقيم Q (المتوسط والوسيط) من المتبرعين والمتبادلين المترافقين.
تجميع البلازميدات ودقتها
تم تجميع تسلسلات البلازميدات باستخدام أوتوسايكلر34. أحجام تسلسلات البلازميدات المجمعة ل 10 عينات مدرجة في الجدول 8، العمود الرابع، وكانت في جميع الحالات متوافقة مع التقديرات المبنية على حركة PFGE (المدرجة في العمود الثالث). تم تقييم الدقة من خلال مقارنة مجموعات القراءة الطويلة المصقولة بميداكا مع تجميعات ميداكا (v2.2.0) بالإضافة إلى تجميعات بولي بوليشية (v0.6.1). تم حساب هوية التسلسل باستخدام dnadiff من MUMmer (الإصدار 3.2.4) (الجدول 8). عبر جميع البلازميدات التي تم تحليلها (n = 12)، كان متوسط طول التجميعات 120.4 كيلو بارت، وكانت مدعومة بمتوسط عمق تسلسل يبلغ 109 أضعاف، مع تغطية متوسطة 97.9٪ عبر تسلسلات البلازميدات. قبل التلميع، كانت التجميعات تظهر متوسط هوية تسلسلية بنسبة 99.74٪ مقارنة بالتجميع المصقول بتسلسل SR.
تقدير تلوث الكروموسومات
لتقدير تغطية البلازميدات والتلوث الكروموسومي، تم ربط القراءات بتسلسلات البلازميدات المجمعة وتسلسل الكروموسومات LMB100 في الإشريكية القولونية . تم محاذاة قراءات نانوبور الأساسية باستخدام minimap2 (الإصدار 2.2.6). تمت معالجة المحاذاة الناتجة باستخدام أدوات SAMtools (الإصدار 1.1.7) لإنشاء ملفات BAM مرتبة وحساب إحصائيات الخرائط. تم تقدير تغطية البلازميدات والعمق المتوسط باستخدام تغطية SAMtools. تفاوت عمق تسلسل البلازميدات بين العينات لكنه ظل أعلى بكثير من الحد الأدنى المطلوب للتجميع الدقيق (50–70 x). جزء من القراءات تم تعيينها على الكروموسوم، مما يشير إلى انتقال الحمض النووي للمضيف أثناء استخراج البلازميدات (الجدول 10).
| البلازميدات | إجمالي القراءات | قراءات الكروموسومات | ٪ قراءة الكروموسومات | متوسط تغطية الكروموسومات | متوسط تغطية البلازميدات | نسبة P:C |
| pBlood_08_0081_2 | 53491 | 20116 | 37.61% | 5 | 1068 | 236 |
| pBlood_d-08_0094 | 160024 | 74827 | 46.76% | 37 | 1618 | 44 |
| pBlood_08_1447 | 238031 | 127967 | 53.76% | 62 | 7842 | 127 |
| pBlood_10_0913_1 | 792263 | 221610 | 27.97% | 87 | 7801 | 90 |
| pBlood_10_0913_2 | 792263 | 221610 | 27.97% | 87 | 7772 | 89 |
| pBlood_11_184_1 | 42494 | 15919 | 37.46% | 7 | 1637 | 250 |
| pBlood_11_184_2 | 42494 | 15919 | 37.46% | 7 | 2188 | 335 |
| pBlood_2011_0238 | 147680 | 80526 | 54.53% | 56 | 2965 | 53 |
| pUPEC_108 | 173409 | 107884 | 62.21% | 61 | 4008 | 65 |
| pUPEC_134 | 105418 | 68645 | 65.12% | 29 | 3640 | 127 |
| pUPEC_271 | 74232 | 53554 | 72.14% | 34 | 1202 | 35 |
| pUPEC_90 | 75193 | 30052 | 39.97% | 16 | 4488 | 278 |
الجدول 10: تقديرات تلوث الحمض النووي الكروموسومي. تم قياس التلوث الكروموسومي في البلازميدات المنقاة من السلالات المستقبلة كنسبة القراءات المرتبطة بتسلسل الكروموسومات LMB100 مقارنة بقراءات البلازميدات. لكل بلازميد، يتم إدراج قراءات البلازميد، والقراءات الكلية، وقراءات الكروموسومات، ومتوسط تغطية الكروموسومات، ومتوسط تغطية البلازميد، ونسبة P:C. لاحظ أن pBlood_08_0081_1 غائب لأنه لم يكن بالإمكان استعادته في المتلقي، وقد أجريت هذه التقديرات على المتلقين بسبب خلفيتهم الكروموسومالية الموحدة.
تم قياس التلوث الكروموسومي كنسبة القراءات التي تربط بتسلسل الكروموسومات المجمع LMB100 مقارنة بقراءات البلازميد. تراوحت معدلات التلوث الكروموسومي بين 72.1٪ و28.0٪، بمتوسط 46.9٪ (الجدول 10). من المرجح أن تبالغ هذه الطريقة في تقدير التلوث الكروموسومي لأن النماذج التي تربط كل من البلازميد والكروموسوم تحتسب ككروموسومية فقط، لكنها لا تزال تشير إلى انتقال كبير للكروموسومات. على الرغم من هذا الانتقال الكروموسومي، نجح سير العمل في التجميع في إعادة بناء تسلسلات البلازميدات الكاملة، والتي كانت لا تزال متوافقة مع الحجم المتوقع من حركة PFGE عبر جميع العينات (الجدول 8، العمودان الثالث والرابع).
أزواج البلازميدات المرافقة الممثلة بين المتبرع والمتلقي
تم عزل السلالة upec271 من عينة بولية بين عامي 2011 و2013 وتم تعيينها إلى مجموعة التطور B2 ونوع التسلسل ST131. كشف تحليل PFGE عن وجود شريط يتوافق مع بلازميد يبلغ حوالي 145 كيلو بيت [الشكل 3A، (المانح)]. حدد تحليل التسلسل قصير القراءة وجود ردود IncF.

الشكل 3: توصيف pUPEC271. (أ) PFGE للمتبرع والمترافق. تظهر أشرطة PFGE للمتبرع (D) وعبر الزوجية (T)؛ تشير الأسهم الخضراء إلى نطاق البلازميد المقابل ل pUPEC271. تم وصف البروتوكول المستخدم ل PFGE في دراسة36. (ب) خريطة pUPEC271. تم توجيه الخاتم المزدوج ضد المتبرع، والإجماع (الذي هو متطابق 100٪) يظهر كحلقة خضراء. تم توليد هذا التمثيل باستخدام بروكسي48. الحلقة الخارجية تظهر موقع CDS، وحلقتان داخليتان تعرضان محتوى GC والميل، والحلقة الداخلية تظهر موقع جينات tRNA. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
تم تأكيد نقل هذا البلازميد، الذي أطلق عليه المؤلفون pUPEC271، بشكل إضافي من خلال تحليل PFGE باستخدام البروتوكول الموصوف سابقا في36 [الشكل 3A، T(ترانزاجينت)]. تم تسلسل الحمض النووي المنقى من سلالة upec271 الممانحة (pUPEC271_donor) باستخدام النانومسام. تماشيا مع حجم نطاقه على هلام PFGE، كشف تسلسل pUPEC271_donor عن بلازميد بحجم 144,519 بيطا (الشكل 3B). كان لهذا البلازميد مترجم من نوع IncF ويحمل جينات مقاومة لست فئات مختلفة من المضادات الحيوية (انظر الجدول التكميلي 1). كما قام المؤلفون بتنقية الحمض النووي من أحد المرافقين (pUPEC271_TC) وحصلوا على تسلسل مطابق لتسلسل المتبرع. يتم رسم التجميع النهائي، مقارنة ب pUPEC271_donor، في الشكل 3B. أكد هذا أن البلازميد التوافقي الموجود في المتبرع تم نقله إلى السلالة المتجاوزة. البيانات الفنية موضحة في الجدول 9.
Blood_10_0913 الإجهاد هو مثال على سلالة تحتوي على أكثر من بلازميد. تم جمع هذه العينة من الدم بين عامي 2008 و2013، وأظهرت القراءة القصيرة وجود اثنين من الربيليكون، أحدهما من نوع IncF والآخر من نوع IncY. حدد تحليل PFGE لهذا السلالة بلازميديتين بحوالي 70 و110 كيلو بايت على التوالي (الشكل 4A). وقد أكد ذلك تسلسل النانوبور، الذي أنتج مجموعتين، بحجم 71,182 قطعة بوصة و111,226 بيتن على التوالي. أطلق المؤلفون على هذين البلازميدين اسم pBlood_10_0913_1 و pBlood_10_0913_2 (الشكل 4B). ومع ذلك، فشلت التجميعات في البداية في الدوران بسبب التشابه التسلسل الواسع في قسم البلازميدات الحامل لآلية الاقتران (تم تسليط الضوء على قسم التجانس الواسع في الشكل 4B). مخطط التجميع الفاشل، الذي تم إنشاؤه باستخدام الضمادة47 ، موضح في الشكل 4C. كان التجميع الكامل للبلازميدتين ممكنا فقط بعد التخلص العدواني من القراءة لتسلسلات القراءة القصيرة (< 3000 نقطة حركة).

الشكل 4: تحديد وتجميع بلازميدات Blood_10_0913 السلالة. (أ) PFGE. تشير الأسهم الخضراء إلى وجود أشرطة بلازميد. البروتوكول المستخدم ل PFGE موصوف في36 (B) تجميع التسلسل. تجميع البلازميدتين الصغيرتين الموجودتين في العينة pBlood_10_0913، pBlood_10_0913_2 و pBlood_10_0913_2. (ج) فشل في التجميع. تمثيل فشل التجميع بسبب وجود تسلسل مشترك بين تجميعين مختلفين. تم توليد هذا التمثيل باستخدام Bandage47. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.
الجدول التكميلي 1: جينات المقاومة التي يحملها pUPEC271. تظهر حجم البلازميدات، والريبليكونات، وجينات المقاومة المرتبطة ب pUPEC271. يرجى الضغط هنا لتحميل هذا الملف.
أوامر إضافية. يرجى الضغط هنا لتحميل هذا الملف.
يعد تحديد وتوصيف البلازميدات ذا أهمية كبيرة للصحة العامة، وعلم الأحياء الدقيقة السريري، وعلم البيئة الميكروبي، حيث تسهل صيانة ونشر الجينات التكيفية عبر البيئات البيئية 1,5 وتعمل كمنصات لتطور متسارع4. سمح إدخال تقنيات NGS بتسلسل عالي الإنتاجية لجينومات العزلات الميكروبية والتجمعات، رغم أن تحديد البلازميدات الكاملة أو الكونتيجات المشتقة من البلازميدات بعد التجميع يتطلب عادة تسلسلات قراءة طويلة لأن بيانات التسلسل القصيرة القراءة غير قادرة على تحديد المناطق المتكررة التي تتجاوز حجم القراءات. أصبح ONT نهجا شائعا جدا للتسلسل طويل القراءة، مع تطبيقات تتجاوز تسلسل الجينوم الكامل وعلم الميتاجينوميات، وتشمل التحقق من تسلسل البلازميدات، والكشف وتوصيف التعديلات اللاجينية، وتسلسل الحمض النووي الريبي، وتقييم عدم الاستقرار الجوهري لمصفوفات الجينات المتتالية (تمت مراجعته في12,49).
هنا، يقدم الباحثون سير عمل لتسلسل البلازميدات باستخدام ONT مصمم للحالات التي يكون فيها تسلسل الكروموسومات معروفا بالفعل أو لا يثير الاهتمام. السبب هو أن استخراج البلازميدات يتطلب جهدا كبيرا، مكلفا، ويستغرق وقتا طويلا (لأن مجموعة المضيف تحتاج إلى زراعة ومعالجة)، وغالبا ما يكون تسلسل الجينوم الكامل باستخدام تقنية النانو ساوما كافيا للحصول على تجميعات بلازميدات كاملة مع توفير السياق الجينومي الخاص بها. تشمل الأمثلة التي يكون فيها البروتوكول مفيدا التحقق من تسلسل البلازميدات في سياق الهندسة الوراثية للبلازميدات 50,51، وتحسين دقة تسلسلات البلازميدات التي تم الحصول عليها بالفعل من قبلWGS 24، أو توصيف البلازميدات التي تم التقاطها من عينات سريرية أو بيئية عن طريق الترفق.
قد يحتاج استخراج الحمض النووي من العينات السريرية والبيئية إلى بعض التخصيص. هنا، يذكر المؤلفون أن بعض السلالات السريرية التي استخدموها أنتجت مصفوفة بيوغشية تسد المرشحات، مما يتطلب عدة خطوات للحبيبات وأحيانا يتطلب إضافة مرشحين. أيضا، لتسلسل البلازميدات المرتبطة، هناك حاجة إلى حجم أكبر من الزراعة مقارنة بالبلازميدات متعددة النسخ المتسلسلة، حيث تميلان البلازميدات المتوسطة والكبيرة الطبيعية إلى أن يكون لديها نسخة بلازميدية منخفضةرقم 2.
سمح البروتوكول بتسلسل البلازميدات بحجم يتراوح بين 4 كيلوبايت و173 كيلوبايت. هذا النطاق يقع ضمن مواصفات الشركة المصنعة (بين 2 و200 كيلوبايت). مؤخرا، نظرت دراسة إلى 23,000 تسلسل بلازميد تم ترسيبه في NCBI ووجدت فروقا كبيرة في الحجم حسب المصدر، حيث أن البلازميدات ذات الأصل البشري لها أصغر حجم وسيط (76 كيلو بل)، والبلازميدات في التربة والنباتات هي الأكبر (215 و427 كيلو بارب على التوالي)؛ البلازميدات من جميع المصادر الأخرى كان حجمها الوسيط بين 79 و147 كيلوبايت، لذا يجب أن يغطي هذا البروتوكول معظم البلازميدات، والتي قد تكون استثناء ميغابلازميدات52 وبلازميدات التربة أو النبات3.
أدرج المؤلفون خطوة إثراء الحمض النووي بالبلازميدات لتحسين دقة تجميعات الحمض النووي للبلازميدات. وجود كميات كبيرة من الحمض النووي الكروموسومي يغمر مسام الخلايا على خلية تدفق التسلسل بقراءات غير مستهدفة، مما يقلل بشكل كبير من عمق حمض البلازميد. تم بالفعل تنفيذ خطوة إثراء الحمض النووي بالبلازميدات لمحاولة التجميع الكامل للحمض النووي البلازميدي في العينات السريرية للكشف عنARGs 24. أظهر هذا تحسنا كبيرا في فعالية إكمال تجميع البلازميدات باستخدام تسلسل النانو-مسام فقط (78٪ من التجميعات المكتملة مقارنة بالمعيار المعياري). تبع ذلك أمثلة أخرى50، 53، 54.
لتحسين دقة القراءة، استخدمت اثنتان من هذه الدراسات تقنية طورتها نانوبور تسمى "استدعاء الأساس المزدوج على القراءات الزائفة الزوجية"50,53. يتضمن هذا النهج محاذاة الإشارات الكهربائية الخام من خيوط الحس والمضاد للحاسة لنفس جزيء الحمض النووي، مما يؤدي إلى تحسين في الدقة بحوالي مرتبة مقدار. الاستدعاء المزدوج على القراءات الزائفة متوافق مع التعدد، طالما أن الاستراتيجية تسمح بتحديد الخيوط الأمامية والخلفية لنفس الجزيء. لذا، لا يمكن استخدامه في البلازميدات المنقاة من عينات ذات مكمل بلازميدي غير معروف أو تحتوي على مزيج من البلازميدات. لاحظ أن عدد البلازميدات في الجينوم عادة ما يتراوح بين 0 و7 بلازميدات (خصوصا في العزلات السريرية)2,55.
لاحظ الباحثون مستوى عاليا من التلوث الكروموسومي في العينات (الجدول 10). تراوح مستوى التلوث في 10 عينات تم تحليلها بين 28 و72٪ بناء على رسم الخرائط القراءة. ومع ذلك، عندما قام الباحثون بتطبيع هذا التسلسل إلى حجم الكروموسوم مقارنة بحجم البلازميد أو البلازميدات الموجودة، تبين أن تسلسلات البلازميدات كانت في المتوسط أكثر وفرة بمقدار 144 مرة من تلك الكروموسوماتية. ترجمة هذا الرقم إلى مستوى إثراء الحمض النووي بالبلازميدات يتطلب معرفة عدد نسخ البلازميدات في كل عينة.
عادة ما تنتج خلايا تدفق MinION مخرج تسلسل ~30 جيجابايت، مما يتيح نظريا تسلسل أكثر من 96 بلازميدا في كل جولة. في الواقع، يتم تقييد التعدد أيضا بحد 96 باركود، مما يجعل 96 بلازميدا الحد الأعلى العملي لكل خلية تدفق. ومع ذلك، فإن تلوث الحمض النووي الكروموسومي واختلال التوازن الطبيعي في الباركود يقلل من عمق القراءة الفعال لكل عينة. لضمان عمق قراءة كاف لكل عينة، فإن النهج الأكثر تحفظا هو مضاعفة حوالي 24 بلازميدا لكل خلية تدفق. العمق الموصى به للتجميع الدقيق لتسلسل البلازميدات يتراوح بين 50 و70×32.56. في الأمثلة التوضيحية المعروضة في هذا المقال، أجرى المؤلفون 21 عينة متعددة الإرسال في خلية تدفق واحدة وحصلوا على عمق متوسط (109x)، وهو يتوافق جيدا مع المعيار الأمثل الموصى به لكل خلية تدفق.
لاحظ أنه مع زيادة عمق القراءة، تكون أقسام القراءة الناجحة/الرسوب أقل فائدة لتقييم جودة القراءة بشكل عام. بدلا من ذلك، توفر المقاييس المستمدة من توزيع الطول تمثيلا أوضح لجودة مجموعة البيانات. أحد هذه المقاييس هو N50، الذي يعرف بأنه طول القراءة الذي تشكل فيه القراءات بهذا الطول أو أطول 50٪ من القواعد الكلية. هذه القيمة حساسة للغاية لخطوات مراقبة الجودة. على سبيل المثال، تصفية القراءات منخفضة الجودة أو القصيرة باستخدام أدوات مثل Filtlong تزيل جزءا كبيرا من الذيل الأقصر للتوزيع، مما يؤدي إلى نفخ N50. هذا التأثير واضح في الجدول 9.
استخدم المؤلفون مجموعة الباركود السريع من ONT، والتي تستخدم ترانسبوسوم لتفكيك حمض نووي البلازميد، وبالتالي لا تتطلب معرفة مسبقة بنمط تقييد البلازميد. وجدت دراسة سابقة أنه بالنسبة للبلازميدات الصغيرة (<20 كيلو بايت)، فإن التجزئة الإنزيمية ضرورية لأن البلازميدات الصغيرة تبقى دائرية أثناء معظم استخراج الحمض النووي، وبالتالي لا تحتوي على نهايات حرة لربط المحول54. تلك الدراسة لم تستخدم الرموز الشريطية، التي تولد نهايات مجانية إضافية من خلال تضخيم PCR، لذا في هذه الحالة، من المحتمل أن يكون هذا نقطة غير مهمة. يتجنب OnRamp، أحد الإجراءات المتقدمة للتحقق الروتيني من البلازميدات، استخدام الرموز الشريطية من خلال الاستفادة من قراءات البلازميد الكاملة للتجميع، مما يبسط تحضير العينة54. في هذه الحالة، فضل المؤلفون الباركود لثلاثة أسباب. الأول هو بناء تجميع بغض النظر عن طول البلازميد. الثانية هي تحمل تلوث الحمض النووي الكروموسومي. الثالث هو أن هذا يفتح إمكانية تعدد عدة عينات، كما ذكر أعلاه.
أشارت دراسة سابقة إلى أن تجميع المجمعات المختلفة حسن النتيجة النهائية من خلال إلغاء التحيزات المدمجةفي منصات التجميع الفردية. هذا التحليل التلوي هو ميزة تم دمجها في برنامج التجميع الذي يستخدمه المؤلفون، Autocycler34. مقارنة بنسخته السابقة المسماة Trycycler57، فإن Autocycler أكثر أتمتة قليلا؛ يستخدم عدة مجموعات فرعية من القراءات، وكما ذكر سابقا، يقوم بعمليات التجميع باستخدام عدة مجمعات لإنشاء تجميع توافقي أكثر دقة. ومع ذلك، يواجه معظم المجمعين، حتى الأفضل، مثل فلاي28، كانو29، ورافين30، صعوبة في وجود البلازميدات ذات التجانس التسلسلي الواسع، أو دمج تسلسلات البلازميدات بشكل خاطئ، أو توليد نسخ متعددة من بلازميد واحد، وهم عرضة لفقدان البلازميدات الصغيرة58,59 (انظر أيضا الشكل 4C). عندما يحدث ذلك، يمكن أن يؤدي حذف القراءات ذات الجودة الأقل وإزالة القراءات القصيرة إلى حل المشاكل الكيميرية.
لتحديد الدقة العامة لبروتوكولات التسلسل، قارن المؤلفون التجميعات التي تعتمد فقط على ONT مع التجميعات المصقولة من خلال الجمع بين القراءة الطويلة والقصيرة القراءة، والتي (باستثناء بيانات PacBio) هي المعيار الذهبي الحالي. كما كان لدى المؤلفين بيانات هلام تظهر عدد البلازميدات الموجودة وأحجامها المقدرة (الجدول 8، العمود الثالث). نجح الباحثون في تجميع جميع البلازميدات الثلاثة عشر الموجودة في العينات العشر وحصلوا على أحجام تتوافق مع التسلسلات المجمعة. وهذا نجاح كبير، حيث غالبا ما تفشل مجمعات تسلسل LR في تجميع البلازميدات ذات التسلسلات الهيكلية المعقدة. على مستوى أكثر تفصيلا، أشارت المقارنات الزوجية بين مجموعتي التسلسلات إلى مستوى دقة متغير (الجدول 8). خمسة من تجميعات البلازميدات كان لها مناطق بدون تغطية في تجميعات ONT، ويفترض أن ذلك نتيجة إضافة تتابع اصطناعي أثناء تجميع التسلسل الخاص ب ONT فقط. من بين السبعة الأخرى، أظهرت ثلاثة منها تعدد أشكال النيوكليوتيد الواحد، مع اختلافات زوجية تراوحت بين 0.0023 و0.05٪، وأظهرت أربعة توافقا مثاليا. لا يبدو أن مستوى الدقة يرتبط بمستوى إثراء أو عمق التسلسل في الحمض النووي بالبلازميد، لكن وجود أكثر من بلازميد واحد في عينة معينة يرتبط، ربما لأن وجود أكثر من بلازميد واحد يميل إلى تضخيم مستوى تعقيد التسلسل (الجدول 8). باختصار، تم تصميم سير العمل لتحسين عائد ودقة تسلسل البلازميدات. بينما يمكنه تجميع تسلسلات البلازميدات، فإنه ينتج مستوى متغيرا من الدقة على مستوى النيوكليوتيدات. يمكن استخدام البروتوكول لتسلسل البلازميدات بأحجام متنوعة، ويتحمل تلوث DNA كروموسومي بنسبة لا تقل عن 72٪، ويمكنه استيعاب (حسب حجم البلازميد وحالة خلية التدفق) ما لا يقل عن 24 بلازميدا في خلية تدفق واحدة.
المؤلفون ليس لديهم ما يكشفون عنه.
نحن ممتنون جدا لستيف ساليبانتي (جامعة واشنطن) لمشاركته بسخاء مجموعة من سلالات الإشريكية القولونية السريرية التي نشرت فيالعدد 45. وقد تم دعم هذا العمل جزئيا من خلال صناديق دعم البحث لتعيين مدير الابتكار وتطوير ريادة الأعمال من جامعة كاليفورنيا الجنوبية إلى كلية السل، ومن خلال منحة جامعة كاليفورنيا أليانزا MX بعنوان "التنميط الجيني لعزلات السل بين الفئات عالية العبء والمعرض للخطر في المكسيك: رؤى حول انتقال السل ومقاومة الأدوية" إلى MC.
| الاسم | الشركة | رقم فهرسي | التعليقات |
|---|---|---|---|
| 1.5 mL Eppendorf DNA LoBind tube | Sigma Aldrich | EP022431081 | يمكن استبدالها بعلامة تجارية أخرى |
| 1.5 mL tube racks | Thermo Fisher Scientific | 22-313630 | يمكن استبدالها بعلامة تجارية أخرى |
| 15 mL conical tubes | FALCON | 352096 | يمكن استبدالها بعلامة تجارية أخرى |
| 15 mL tube racks | Thermo Fisher Scientific | 8850 | يمكن استبدالها بعلامة تجارية أخرى |
| 2.0 mL microcentrifuge tubes | Fisherbrand | 05-408-138 | يمكن استبدالها بعلامة تجارية أخرى |
| 250 mL sterile flasks | PYREX | 5320 | يمكن استبدالها بعلامة تجارية أخرى |
| Carbenicillin | GOLD BIOTECHNOLOGY | C-103-50 | يمكن استبدالها بعلامة تجارية أخرى |
| Centrifuge | Eppendorf | 5424R-PR-R | يمكن استبدالها بعلامة تجارية أخرى |
| Centrifuge | Eppendorf | 5810R-P-R | يمكن استبدالها بعلامة تجارية أخرى |
| Centrifuge bottles | Sigma Aldrich | B1283-4EA | يمكن استبدالها بعلامة تجارية أخرى |
| Disposable Inoculating Loops: 1 µL | Fisherbrand | 22-363-595 | يمكن استبدالها بعلامة تجارية أخرى |
| Electronic Digital Scale | Denver Instrument | APX-2001 | يمكن استبدالها بعلامة تجارية أخرى |
| Ethanol | Fisher Scientific | 04-355-720 | يمكن استبدالها بعلامة تجارية أخرى |
| Ethanol-resistant markers | Sharpie | 37001; 37002; 37003 | يمكن استبدالها بعلامة تجارية أخرى |
| Flow cell (R10.4.1) | Oxford Nanopore Technologies | FLO-MIN114 | |
| Flow Cell Wash Kit | Oxford Nanopore Technologies | EXP-WSH004 | |
| Gloves | X-GEN | 44-100M | أي علامة تجارية للقفازات النيتريلية تعمل |
| Ice bucket | Thermo Fisher Scientific | 432128 | يمكن استبدالها بعلامة تجارية أخرى |
| Isopropanol | Fisher Scientific | BP26184 | يمكن استبدالها بعلامة تجارية أخرى |
| MinION Sequencing Device | Oxford Nanopore Technologies | MIN-101B | |
| Native Barcoding Kit 96 V14 | Oxford Nanopore Technologies | SQK-NBD114.96 | |
| NucleoBond Xtra MIDI Plasmid DNA Extraction Kit | Takara | 740410.5 | |
| PCR tube racks | AXYGEN | R96PCRFSP | يمكن استبدالها بعلامة تجارية أخرى |
| PCR tubes: (0.2 mL) | AXYGEN | PCR-02-C | يمكن استبدالها بعلامة تجارية أخرى |
| Program autocycler v0.5.2 | Ryan Wick / University of Melbourne | ||
| Program bcftools v1.22 | Wellcome Sanger Institute (HTSlib / Samtools team) | ||
| Program bwa v0.7.19 | Heng Li, originally at the Broad Institute | ||
| Program canu v2.3 | University of California, Davis & Pacific Biosciences collaborators | ||
| Program flye v2.9.6 | University of California, San Diego (Pavel Pevzner lab) | ||
| Program medaka v2.1.1 | Oxford Nanopore Technologies | ||
| Program metamdbg v1.2 | INRIA | ||
| Program miniasm v0.3 | Heng Li | ||
| Program minimap2 v2.28 | Heng Li | ||
| Program minipolish v0.2.0 | Ryan Wick / University of Melbourne | ||
| Program porechop v0.2.4 | Ryan Wick / University of Melbourne | ||
| Program racon v1.5.0 | Genome Institute of Singapore | ||
| Program samtools v1.22.1 | Wellcome Sanger Institute (HTSlib / Samtools project) | ||
| Qubit™ 4 Fluorometer | Invitrogen | Q33226 | |
| Qubit™ Assay Tubes | Invitrogen | Q32856 | |
| Rapid Barcoding Kit 24 and 96 V14 | Oxford Nanopore Technologies | SQK-RBK114.96 | |
| Set of micropipettes that dispense between 1 – 10 μL (P10), 2 – 20 μL (P20), 20 – 200 μL (P200) and 200 – 1000 μL (P1000) | RAININ | 17008648; 17008650; 17008652; 17008653 | يجب معايرة الماصات الدقيقة |
| Spectrophotometer | Thermo Scientific | 335905P | يمكن استبدالها بعلامة تجارية أخرى |
| Sterilized tips for 10 μL, 200 μL and 1000 μL | Eclipse | 1011-260-000-9; 1018-260-000; 1019-260-000-9 | |
| Thermal cycler C1000 Touch | Bio-Rad | 1851196 | يمكن استبدالها بعلامة تجارية أخرى |
| Water bath | Thermo Scientific | 51221052 | يمكن استبدالها بعلامة تجارية أخرى |