اتباع ديناميات المتغيرات الهيكلية في المجموعات السكانية المتطورة تجريبيا

Summary

لقد طورنا طريقة فعالة من حيث التكلفة لمتابعة ديناميكيات أليل تعدد أشكال النوكليوتيدات غير المفردة والتي يمكن تكييفها بسهولة مع أرشيفات التطور التجريبية المجمدة. اقترنت تقنية تفاعل البوليميراز المتسلسل الثلاثي مع الرحلان الكهربائي الشعري المتوازي الآلي لتحديد التردد النسبي لأليل الإدخال على مدار التطور التجريبي.

Abstract

من المعروف الآن أن المتغيرات الهيكلية (SVs) (أي الحذف والإدراج والازدواجية والانقلابات) تلعب دورا مهما في تباين النمط الظاهري ، وبالتالي في عمليات مثل تحديد المرض أو التكيف مع بيئة جديدة. ومع ذلك ، فإن متغيرات النوكليوتيدات المفردة تحظى باهتمام أكبر بكثير من SVs ، ربما لأنها أسهل في اكتشافها ، ومن السهل التنبؤ بآثارها المظهرية. وقد أدى تطوير تكنولوجيات التسلسل العميق القصير والطويل القراءة إلى تحسين الكشف عن التسلسل العميق، ولكن التحديد الكمي لتواترها من بيانات التسلسل المجمع (poolseq) لا يزال معقدا ومكلفا من الناحية التقنية.

هنا ، نقدم طريقة بسيطة وغير مكلفة إلى حد ما ، والتي تسمح للباحثين بمتابعة ديناميات تردد أليل SV. كمثال على التطبيق ، نتبع تكرار إدخال تسلسل الإدراج (IS) في مجموعات التطور التجريبي للبكتيريا. تعتمد هذه الطريقة على تصميم ثلاثة توائم من الاشعال حول حدود المتغير الهيكلي ، بحيث تختلف الأمبليكونات الناتجة عن تضخيم النوع البري (WT) والأليلات المشتقة في الحجم بنسبة 5٪ على الأقل ، وأن كفاءة تضخيمها متشابهة. ثم يتم تحديد كمية كل amplicon عن طريق الكهربائي الشعري المتوازي وتطبيعها إلى منحنى المعايرة. يمكن توسيع هذه الطريقة بسهولة لتشمل التحديد الكمي لتواتر المتغيرات الهيكلية الأخرى (الحذف والازدواجية والانقلاب) ونهج التجميع التسلسلي للسكان الطبيعيين ، بما في ذلك مجموعات مسببات الأمراض داخل المريض.

Introduction

المتغيرات الهيكلية (SVs) هي تغييرات في التسلسل الجيني ، تؤثر بشكل عام على 50 نقطة أساس أو أكثر. الفئات الأربع من SVs الموصوفة هي عمليات إدراج كبيرة وعمليات حذف كبيرة وانعكاسات وازدواجية. حتى وقت قريب ، تم تكريس المزيد من الاهتمام لمتغيرات النوكليوتيدات المفردة (SNVs) مقارنة بالمتغيرات الهيكلية ، من حيث تأثيراتها المظهرية ودورها كمحددات وراثية للمرض ، أو مساهمتها في التكيف. ربما يرجع ذلك إلى أنه من الأسهل اكتشاف SNVs والتنبؤ بآثارها المظهرية. ومع ذلك ، فإن تقنيات التسلسل العميق قصيرة وطويلة القراءة قد حسنت بقوة اكتشاف SVs ، على الأقل في الجينومات الفردية أو النسيلية¹. في موازاة ذلك ، تم توصيف آثارها المظهرية بشكل أفضل ، وتم توثيق العديد من الأمثلة على آثارها كمحددات وراثية للأمراض البشرية ^2,3 أو التكيف مع بيئة جديدة⁴.

عمليات الحذف والإدراج ، غالبا بسبب إدخالات العناصر الجينية المتنقلة (MGE) ، أكثر تعطيلا بكثير من تعدد أشكال النوكليوتيدات المفردة (SNPs) وتؤدي إلى طفرات إزاحة الإطار وتعديلات بنية البروتين. يؤدي الحذف وإدخال MGE داخل الجينات دائما تقريبا إلى تعطيل الجينات ، ويمكن أن يؤدي الإدراج في المناطق غير المشفرة إلى قمع أو تعبير تأسيسي للجينات المجاورة عندما تحتوي تسلسلات الإدراج (ISs) على تسلسلات المروج أو الإنهاء⁵. في حين أن خروج الجينات الأساسية يؤدي إلى آثار ضارة واضحة على اللياقة البكتيرية ، فإن فقدان الجينات غير الأساسية مفيد في بعض الحالات. وعلى الرغم من التكاليف المتأصلة في الازدواجية، فإنها يمكن أن تكون مفيدة أيضا، وتشارك في التكيف لأنها تؤدي إلى تغيير في جرعة الجينات؛ يمكن أن تكون الزيادة في نشاط بروتين معين مفيدة اعتمادا على الظروف⁶.

عادة ما تبدأ مجموعات التطور التجريبي الميكروبي بالحيوانات المستنسخة. هذا الغياب الأولي للتنوع الجيني ، جنبا إلى جنب مع خاصية "البيئة المغلقة" لأنابيب الاختبار ، يؤدي إلى إمكانات محدودة للغاية للتطور عن طريق اكتساب الجينات من خلال نقل الجينات الأفقية وإعادة التركيب. في هذه الظروف المحددة ، تكون المساهمة في تكييف عمليات الحذف والازدواجية وإدخال MGE داخل الجينوم مهمة بشكل خاص ؛ غالبا ما تتكيف البكتيريا من خلال طفرات فقدان الوظيفة (ويرجع ذلك أساسا إلى عمليات الحذف أو إدخال MGE) ، مما يؤثر على الجينات غير المفيدة في البيئات الاصطناعية المستقرة ، الغنية بالمغذيات^{في كثير من الأحيان ، 7}. في أطول تجربة لتطور الإشريكية القولونية ، كانت عمليات إدخال IS150 متكررة بشكل خاص بين المجموعات السكانية التي تطورت بعد 50000 جيل ، حيث تمثل عناصر IS 35٪ من الطفرات التي تصل إلى تردد عال في السكان الذين يحتفظون بمعدل طفرة نقطة أسلافهم⁸.

دراسات التطور وإعادة التسلسل تجمع بين التطور التجريبي وتقنيات التسلسل من الجيل التالي (NGS) للتحقيق في كيفية تكيف البكتيريا ، على المستويين الظاهري والجينومي ، مع الظروف والضغوط البيئية المختلفة ، مثل مصادر الكربون والطاقة المختلفة والمضادات الحيوية والإجهاد التناضحي⁹،^10،11 . تحصل هذه الدراسات عادة على معلومات جينومية عن المجموعات السكانية المتطورة أو الحيوانات المستنسخة فقط عند نقطة النهاية التجريبية ، وفي بعض الحالات ، في عدد من النقاط الزمنية الوسيطة¹²،^13،14. توفر هذه البيانات نظرة ثاقبة للجينات والمسارات المشاركة في التكيف مع بيئة معينة ، ولكنها نادرا ما تسمح للباحثين بمتابعة ديناميكيات الأليلات الناشئة والكاسحة بمرور الوقت.

تتمثل إحدى الطرق التي يجب اتباعها في هذه الديناميكيات في اختيار عدد محدود من الأليلات المنفصلة ذات الأهمية (بسبب وظيفة الجينات التي تؤثر عليها ، لأنها تكتسح بالتوازي في مجموعات سكانية مستقلة ، وما إلى ذلك) واستخدام تسلسل amplicon لتحديد نسبة الأليل ، وتجميع العديد من النقاط الزمنية في نفس تشغيل التسلسل¹⁵. تم استخدام هذه الطريقة بنجاح لمتابعة ديناميكيات المتغيرات صغيرة الحجم (SNPs أو 1 bp indels) في¹⁶ و¹⁷ مجموعة تجريبية من الميكروبات. ومع ذلك ، في حالة إدخالات indels أو MGE الأكبر حجما ، فإن اختلاف حجم الأمبليكونات يؤدي إلى اختلافات كفاءة تفاعل البوليميراز المتسلسل ، مما يشوه العلاقة بين نسب القراءة والأليل . في بعض الحالات ، يكون فرق الحجم بين الأليلين أعلى من الطول الكلاسيكي للأخمبليون. هنا ، قمنا بإقران تقنية تفاعل البوليميراز المتسلسل الثلاثي مع الرحلان الكهربائي الشعري المتوازي الآلي لتحديد التردد النسبي لأليل الإدخال بناء على تمييز الحجم. يسمح هذا النهج باستغلال النقاط الزمنية التجريبية غير المستغلة لتحديد ديناميات أليل طافر ناشئ ومتابعة تواتره للتثبيت أو الفقد ، بطريقة فعالة من حيث التكلفة. طبقنا هذه الطريقة لتتبع الأليلات الطافرة الناشئة ، التي تحورت من خلال إدخال IS10 ، مما يوفر للنمط الجيني المتحور نمطا ظاهريا مفرط التحور.

تتطلب هذه الطريقة أليلين مستهدفين باختلاف ≥5٪ في الحجم. أولا ، تم تصميم ثلاثة توائم أولية لإنتاج شظايا متشابهة الحجم ، والتي تشترك في التمهيدي المشترك. ثانيا ، يتم تحسين ظروف تفاعل البوليميراز المتسلسل ، ويتم إنتاج منحنى معايرة باستخدام مزيج من النوع البري (WT) والحمض النووي الطافر. أخيرا ، يتم تضخيم العينات بواسطة تفاعل البوليميراز المتسلسل ، ويتم قياس التردد النسبي لكل أليل عن طريق الرحلان الكهربائي الشعري الكمي المتوازي.

Protocol

يتطلب إعداد هذا البروتوكول معرفة دقيقة بنقطة الإدراج أو الحذف أو الانعكاس أو الازدواجية داخل تسلسل الأجداد. عادة ما يتم الحصول على هذه المعلومات عن طريق تسلسل الجينوم الكامل (WGS) لعينات النقطة النهائية أو الوسيطة. في البروتوكول التالي ، يتم إعطاء المبدأ العام لحالة طفرة الإدراج لكل خطوة ، جنبا إلى جنب مع حالة تمثيلية حيث يتم اتباع تكرار إدخال IS10 في جين mutS في مجموعة تطور تجريبية من E. coli. في هذه الفئة من السكان، حددت WGS لمجتمع نقطة النهاية إدخال 1,329 bp IS10 بين المواضع 2,463 و 2,471، مما أدى إلى تكرار موقع الإدراج هذا. تنطبق هذه الطريقة على أنواع SV الثلاثة الأخرى ، ويتم إعطاء خصوصيات كل حالة في المناقشة.

1. تصميم الاشعال الثلاثي

1. استخدم ممارسات تصميم التمهيدي الكلاسيكية (18-24 نقطة أساس ، 40٪ -60٪ محتوى GC ، ابدأ / تنتهي بأزواج G / C حيثما أمكن ، فرق T_m < 5 درجات مئوية) لإنتاج البادئات FW1 و RV1. صمم مواد أولية لتضخيم أمبليكون قصير على أليل WT حول موقع إدخال الأليل الطافر (الشكل 1).
   ملاحظة: يمكن أن يتراوح حجم amplicon من 100 bp إلى ما يصل إلى 3000 bp ، بما يتماشى مع سلم حجم الحمض النووي المستخدم في الرحلان الكهربائي الشعري. في هذا المثال ، تم تضخيم amplicon 155 bp. يمنع حجم الشظية الصغير المختار هنا التضخيم خارج الهدف لتسلسل إدخال IS10 بأكمله (انظر القسم 2).
2. صمم برايمر أمامي ثان FW2 ضمن تسلسل الإدراج ، لإنتاج أمبليكون ثان أكبر أو أصغر بنسبة 5٪ تقريبا من أمبليكون WT (الشكل 1). هذا الاختلاف في الحجم بنسبة 5٪ هو الحد الأدنى لفرق الحجم الذي يمكن لجهاز الرحلان الكهربائي الشعري المتوازي تمييزه بشكل موثوق. لذلك ، قم بتصميم الاشعال بحيث يكون للمضخمين اختلاف في الحجم ، وهو أعلى ولكن أقرب ما يمكن إلى العتبة النسبية.
   ملاحظة: تأكد من تقليل فرق T_m وتشكيل dimer التمهيدي. في هذا المثال ، تم تصميم التمهيدي الأمامي الثاني لإنتاج 226 bp amplicon ، 71 bp أكبر من WT amplicon. في هذا المثال التمثيلي ، تكون تسلسلات التمهيدي كما يلي:
   FW1:AAAGCATTTCCCGAACGCC
   RV1: GCGATAAATCCACTCCAGCGCC
   FW2: AGTTCGCTTAGGCATGGAAG

الشكل 1: مخطط تصميم التمهيدي الثلاثي على إدخال mutS WT وإدخال mutS IS10 الطافر. يمثل المثلث الأسود موقع إدخال IS10 في جين mutS. جين WT باللون الأزرق ، و IS10 باللون البرتقالي. تشير البادئات FW1 و RV1 إلى موقع الإدراج IS10 وتنتج مضخم 155 نقطة أساس WT. ينتج التمهيدي RV1 والتمهيدي داخل IS10 FW2 أمبليكون ثان 226 bp. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

2. تحسين ظروف PCR

1. تنمو ثقافة بين عشية وضحاها من WT الثابتة واستنساخ أليل متحولة.
2. استخراج الحمض النووي باستخدام أي مجموعة.
3. تحديد كمية الحمض النووي.
4. تحضير عينة من الحمض النووي عن طريق تخفيف استخراج WT والحمض النووي الطافر إلى 5 نانوغرام / ميكرولتر. امزج عينتي الحمض النووي بنسبة 50/50.
5. تضخيم 10 نانوغرام من عينات الحمض النووي الثلاثة (WT ، متحولة ، مزيج 50/50) باستخدام مزيج رئيسي PCR جاهز للاستخدام 2x ، و 0.5 μM FW1 التمهيدي ، و 1 μM RV1 التمهيدي ، و 0.5 μM FW2 التمهيدي في حجم تفاعل 20 ميكرولتر. استخدم جل الأغاروز بنسبة 2٪ لترحيل منتج PCR عن طريق الرحلان الكهربائي الكلاسيكي ، وتحديد ظروف PCR المثلى.
   ملاحظة: يجب تقليل أوقات الاستطالة لمنع تكوين أمبليكون FW1 و RV1 على الأليل الطافر. يجب ضبط درجة حرارة التلدين لتقليل التضخيم المتحيز للأليلات والتضخيم غير المحدد.
   1. لمتابعة البرنامج في هذا المثال ، استخدم الإعدادات التالية: 98 درجة مئوية لمدة 10 ثوان ، متبوعة ب 25 دورة من 98 درجة مئوية لمدة 1 ثانية ، 58 درجة مئوية لمدة 15 ثانية ، 72 درجة مئوية لمدة 8 ثوان ، وخطوة استطالة نهائية عند 72 درجة مئوية لمدة 1 دقيقة.
      ملاحظة: تم تقليل وقت الاستطالة إلى 8 ثوان لمنع تضخيم منتج >1000 bp من التمهيدي الأمامي والتمهيدي RV1 على الأليل الطافر (mutS مع إدخال IS10).

3. منحنى المعايرة

1. امزج عينتي الحمض النووي ، WT و mutant ، بنسب 10/90 و 25/75 و 40/60 و 50/50 و 60/40 و 75/25 و 90/10.
   ملاحظة: يتم تحضير النسخ المتماثلة البيولوجية من مزارع بكتيرية مستقلة بين عشية وضحاها.
2. تضخيم باستخدام ظروف PCR المحسنة (انظر القسم 2).
3. تحديد كمية منتج الأمبليكون.
4. تمييع منتجات PCR إلى 0.1 نانوغرام / ميكرولتر.
5. إعداد أداة الكهربائي الشعري المتوازي.
   1. امزج الجل الطازج والصبغة (مجموعة التحليل الكمي NGS (22 أو 33 أو 55) ؛ جزء HS NGS 1-6000 نقطة أساس لهذا المثال).
      ملاحظة: انظر دليل شظية الرحلان الكهربائي الشعري المتوازي HS NGS للحصول على تعليمات مفصلة (انظر جدول المواد).
6. استبدل محلول التخزين الشعري والمخزن المؤقت للمدخل ، وضع لوحة الشطف العازلة في مواقع الأدراج الصحيحة لأداة الرحلان الكهربائي الشعري المتوازي.
7. أضف 2 ميكرولتر من علامة مخفف HS إلى 22 ميكرولتر من كل عينة مخففة في لوحة 96 بئرا.
8. أضف سلم حجم (سلم حجم الحمض النووي ؛ النطاق 1-6000 نقطة أساس) من مجموعة التحليل الكمي HS NGS إلى بئر واحد من لوحة 96 بئرا.
9. ضع لوحة 96 بئرا في الدرج الصحيح لأداة الرحلان الكهربائي الشعري المتوازي ، وحدد تشغيل على برنامج أداة الرحلان الكهربائي الشعري المتوازي.
10. تحليل النتائج باستخدام برنامج تحليل البيانات ، الذي يكتشف ويحدد كل قمة من سلم الحجم ، ويخصص كل قمة من العينات لحجمها الفعلي المعروف.
    1. عند استخدام مجموعة كمية ، استخدم البرنامج لتحديد كمية الحمض النووي لكل جزء من خلال دمج المنطقة تحت الذروة ، كما هو الحال في تحليل بيانات الكروماتوغرافيا. مرة أخرى ، قارن العينات بكميات معروفة في المعايير لتحديد كل قمة من عينة ، وحساب النسب بين القمم المختلفة المكتشفة في العينات.
    2. بناء منحنى معايرة (الشكل 2) ، يربط النسبة المعروفة للأليل الطافر (مزيج الحمض النووي) بتلك المقاسة باستخدام أداة الرحلان الكهربائي الشعري المتوازي. يسمح منحنى المعايرة هذا بتقييم موثوقية الطريقة وتصحيحها لتحيز التضخيم الطفيف.

4. إعداد عينة

1. تنمو عينات النقطة الزمنية بين عشية وضحاها في الظروف القياسية.
2. استخراج الحمض النووي.
3. تحديد كمية الحمض النووي.
4. تضخيم العينات باستخدام ظروف تفاعل البوليميراز المتسلسل المحسنة (انظر القسم 2).
5. قم بتشغيل العينات في أداة الرحلان الكهربائي الشعري المتوازي (انظر الخطوات 3.5-3.10).

5. قياس كمية الأليل

1. استخراج كميات الأليل الطافرة من بيانات أداة الرحلان الكهربائي الشعري المتوازي مع البرنامج ، وحساب النسب الفعلية عن طريق رسم هذه القيم على منحنى المعايرة.

Results

باستخدام الحمض النووي المستخرج من استنساخ الأسلاف واستنساخ مفرط الحركة المعزول من مجموعة S2.11 عند الجيل 1000 ، أنشأنا منحنى المعايرة الموضح في الشكل 2. تم ربط النسب الطافرة الفعلية من خلطات الحمض النووي المعدة في المختبر والمقاسة بواسطة أداة الرحلان الكهربائي الشعري المتوازي بعلاقة خطية من الميل 1.0706 ، مع R² من 0.9705. بالإضافة إلى ذلك ، كان هناك اتفاق جيد بين النسخ البيولوجية. كان الانحراف المعياري بين 0.61 و 17.74 عبر النقاط التسع للمنحنى المعياري.

Discussion

هنا ، اقترحنا طريقة فعالة من حيث التكلفة تسمح باتباع ديناميكيات أليلات SV التكيفية الناشئة في مجموعات التطور التجريبية. تجمع هذه الطريقة بين تقنيات تفاعل البوليميراز المتسلسل الكلاسيكية والرحلان الكهربائي الشعري المتوازي الآلي ، مما يسمح بتحديد الكميات النسبية لألينتين. بمجرد إعداده ، فإنه يسمح بتحديد نسب الأليل في العديد من العينات بالتوازي ، وهو أقل تكلفة بكثير من WGS. يمكن اعتبار هذه الطريقة مكافئة لتسلسل amplicon للطفرات غير SNP وكحل للقيود التقنية لتسلسل amplicon ل SVs الكبيرة. لقد أظهرنا أن الطريقة تحدد بدقة نسبة ا...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
96 لوحة PCR ذات التفاف الجيد	4titude	4Ti - 0740	PCR
Agarose Molecular biology grade	Eurogentec	EP-0010-05	Agarose gel electrophorese
Agilent DNF-474 HS NGS Fragment Kit دليل سريع لأنظمة محلل الشظايا	دليل تعليمات Agilent		PDF
Buffer TBE	Panreac appliChem	A4228،5000Pc	Agarose gel electrophoreses
حلقات التلقيح والإبر المعايرة التي يمكن التخلص	LABELIANS	8175CSR40H	الثقافة البكتيرية
Dneasy الدم والأنسجة كيت	Qiagen	69506	استخراج الحمض
النووي الرحلان الكهربائي مزود الطاقة	Amilabo	ST606T	الرحلان الكهربائي جل الاغاروز   ؛
محلل الشظايا نظام CE الآلي	الرحلان الكهربائي الشعري المتوازي	الرشيقة
جزيء سلم الحمض النووي	Agilent	DNF-396 ، النطاق 1-6000bp	الرحلان الكهربائي الشعري المتوازي
GENTAMICIN SUL SALTREAGENT	SIGMA-Aldrich	G1264-1G	ثقافة بكتيرية
علامة مخففة عالية	الحساسية Agilent	DNF-373	الشعيرات الدموية المتوازية الرحلان الكهربائي
عالي الحساسية NGS طقم التحليل الكمي	Agilent	DNF-474	الرحلان الكهربائي الشعري المتوازي
سلم تحميل سريع 1 كيلو بايت بالإضافة إلى سلم الحمض النووي	NEB	N0469S	الرحلان الكهربائي لهلام الاغاروز
مرق LB ، VegitoneNutriSelect Plus	Millipore	28713	ثقافة بكتيرية
ماستر ميكس PCR عالية الدقة Phusion فلاش	Thermo Fisher Scientific	F548L	PCR
Primers	Eurogentec		PCR
Prosize تحليل البيانات برنامج v.4	Agilent	V.4	الرحلان الكهربائي الشعري المتوازي
Qubit فحوصات	Invitrogen	MAN0010876	القياس الكمي للحمض النووي
Qubit dsDNA HS ASE KSE	LIFE TECHNOLOGIES SAS	Q32854	القياس الكمي للحمض النووي
Thermocycler	Eppendorf	Ep تدرجات	PCR
UVbox ، eBOX VX5	Vilber Lourmat		Agarose هلام التصور الكهربائي
الماء للتحضير عن طريق الحقن	Aguettant	PROAMP	PCR