تم زراعة البكتيريا من البول الذي تم جمعه من النساء بالتراضي كجزء من الدراسات المعتمدة من مجلس المراجعة المؤسسية 19MR0011 (UTD) وS STU 032016-006 (UTSW).
1. تعديل تعزيز ثقافة البول
ملاحظة: يجب تنفيذ كافة الخطوات الثقافية في ظل ظروف معقمة. تعقيم جميع الأدوات والحلول ووسائل الإعلام. تنظيف منطقة العمل مع الإيثانول 70٪، ثم إعداد الموقد بونسن والعمل بعناية على مقربة من اللهب للحد من فرص التلوث. وبالتناوب، يمكن استخدام خزانة السلامة البيولوجية من الفئة الثانية للحفاظ على بيئة معقمة. ارتداء معدات الحماية الشخصية المناسبة (PPE) لتجنب التعرض للميكروبات المسببة للأمراض المحتملة.
- طلاء البول الجلسرين مخزنة وعزل مستعمرة
- إذابة البول المجهز بالجلسرين في درجة حرارة الغرفة (RT). مرة واحدة إذابة، دوامة العينة لمدة 5 ق لخلط. في أنابيب الطرد المركزي الدقيق المعقمة، قم بإعداد 1:3 و 1:30 تمييع البول في محلول ملحي معقم 1x الفوسفات العازل (PBS) إلى حجم نهائي قدره 100 ميكرولتر.
ملاحظة: يتم إعداد البول الجلسرين المخزنة عن طريق خلط 500 ميكرولتر من البول غير المخفف و 500 ميكرولتر من 50٪ الجلسرين المعقم في cryovials وتخزينها في -80 درجة مئوية.
- لوحات أجار قبل الحارة في 37 درجة مئوية لمدة 15 دقيقة قبل الاستخدام. يرجى الاطلاع على الشكل 1 لأنواع الوسائط وظروف الثقافة المناسبة للأجناس البكتيرية البولية الشائعة. خلط البول المخفف جيدا عن طريق pipetting قبل الطلاء، لوحة 100 μL من البول المخفف على لوحة أجار المطلوب ونشر العينة باستخدام الخرز الزجاج العقيم. لوحة 100 ميكرولتر من 1x برنامج تلفزيوني diluent على لوحة منفصلة كما لا تحكم النمو.
ملاحظة: إذا كان محاولة زراعة الأنواع الشائعة المسالك البولية (على سبيل المثال، إشريكيا القولونية، كليبسيلا spp.، Enterococcus faecalis، الخ)، فمن المستحسن استخدام أجار كروموجينيك (جدول المواد) كما أنه يسمح بسهولة تحديد الأنواع البكتيرية المسالك البولية(الشكل 1). حمض Colistin Nalidixic (CNA) أو MRS agar مفيدة لعزل الأنواع الإيجابية للجرام (على سبيل المثال ، Lactobacillus spp.) من البول المعروف أنه يحتوي على الاوروباتوغنز السلبي للجرام ، والذي قد يتفوق على الأنواع السريعة في أجار غير انتقائية.
- احتضان لوحة مقلوب في حالة الغلاف الجوي المطلوب في 35 درجة مئوية لمدة 24 ساعة لuropathogens و 3-5 أيام للبكتيريا fastidious (الشكل 1).
- بعد فترة الحضانة، قم بإزالة الصفائح من الحاضنة. من كل لوحة، اختر المستعمرات التي تظهر لونا فريدا أو مورفولوجيا أو أنماط انحللية.
- إعادة خط مستعمرة بكتيرية باستخدام حلقة معقمة على أجار المقابلة واحتضان لوحة مقلوب لمدة 2-5 أيام في الغلاف الجوي المطلوب للحصول على مستعمرات معزولة جيدا.
ملاحظة: إذا كان استخدام BAP للثقافة الأولية، قد الترقيع المستعمرات على أجار الكروموجينيك توفر معلومات مفيدة حول عدم التجانس من السكان البكتيرية في العينة.
- زراعة في مرق السائل والجلسرين تخزين البكتيرية يعزل
- بمجرد الحصول على المستعمرات المعزولة التي تتطابق مع مورفولوجيا المستعمرة الأم ، اختر مستعمرة واحدة وتلقيحها في 3 مل من مرق السائل باستخدام حلقة تلقيح معقمة. الرجوع إلى الشكل 1 للمرق قادرة على دعم نمو أجناس الميكروبات البولية المشتركة. ختم لوحات أجار مع البارافيلم وتخزينها في 4 درجة مئوية لمدة 2-4 أيام. احتضان الثقافات السائلة في الظروف الجوية المرغوبة لمدة 1-5 أيام حتى تكون الثقافة عكرة بشكل واضح.
- بعد أن لوحظ النمو, دوامة الثقافة, ومن ثم إضافة 1 مل من الثقافة بين عشية وضحاها إلى 500 ميكرولتر من العقيمة 50٪ الجلسيرول في 2 مل cryovial; ختم ومزيج بلطف عن طريق عكس. إعداد اثنين من مخزون الجلسرين لكل مستعمرة (واحد بمثابة نسخة احتياطية) وتخزينها في -80 درجة مئوية.
2. تحديد الأنواع البكتيرية من قبل 16S rRNA جين سانجر التسلسل
ملاحظة: يمكن تأكيد الهوية الميكروبية بدلا من ذلك باستخدام الليزر بمساعدة مصفوفة Desorption وقت تأين قياس الطيف كتلة الطيران (MALDI-TOF)20.
- سلسلة تفاعل Colony-Polymerase المتسلسلة (PCR)
- إعداد 25 ميكرولتر من تفاعل PCR في أنابيب PCR بإضافة 12.5 ميكرولتر من 2x Taq Polymerase ماستر ميكس، 0.5 ميكرولتر من 10 ميكرومتر 8F التمهيدي، 0.5 ميكرولتر من 10 ميكرومتر 1492R التمهيدي(جدول المواد)،و 11.5 ميكرولتر من المياه الخالية من النيوكليز21.
ملاحظة: إذا كان أداء PCR لعينات متعددة، وجعل مزيج رد فعل رئيسي من مزيج Taq بوليمراز، التمهيديات، والمياه الخالية من النيوكليز العقيمة. ثم aliquot 25 ميكرولتر في كل أنبوب PCR.
- لأداء مستعمرة PCR، انتقد مستعمرة معزولة جيدا من إعادة خط باستخدام مسواك معقمة أو تلميح ماصة. Resuspend المستعمرة في مزيج رد فعل PCR أعدت في الخطوة 2.1.1. مزيج بلطف. جمع السائل في الجزء السفلي من الأنبوب عن طريق تدور سريعة في 2000 x ز.
ملاحظة: تأكد من خلو العينة من فقاعات الهواء. تضمين نموذج عنصر تحكم بلا قالب (NTC) يحتوي على مزيج تفاعل PCR وحده.
- ضع أنابيب العينة في المدور الحراري واركض البرنامج التالي: 95 درجة مئوية لمدة 3 دقائق؛ 40 دورة من: 95 درجة مئوية لمدة 30 s، 51 °C لمدة 30 s، و 72 °C لمدة 1 دقيقة 30 s؛ 72 درجة مئوية لمدة 10 دقائق؛ عقد عند 10 °C.
- استخراج الجل وتحديد الأنواع
- عند الانتهاء من تشغيل PCR، تحقق من المنتج PCR على هلام agarose 1٪ أعدت في 0.5x تريس بورات-EDTA (TBE) العازلة. قبل صب الجل، إضافة بروميد الإيثيديوم (EtBr). ثم، يلقي هلام باستخدام أمشاط لل آبار التي تحمل ما لا يقل عن 20 ميكرولتر حجم العينة.
تنبيه: EtBr هو عامل تداخل يشتبه في أنه مسرطن. ارتدي دائما قفازات وPPE عند التعامل معها والتخلص من المواد التي تحتوي على EtBr وفقا لإرشادات المؤسسة.
- عندما يتم تعيين هلام، وضع هلام في خزان الكهربائي مليئة 0.5x TBE العازلة وإزالة المشط. تحميل سلم 1 كيلوبايت في البئر الأول و 10-20 ميكرولتر من تفاعل PCR في الآبار اللاحقة. تشغيل في 100-140 V حتى حلها. تصور هلام تحت ضوء الأشعة فوق البنفسجية وتأكيد وجود نطاق محدد بوضوح في ~ 1.5 كيلو بايت التي هي غائبة في بئر المجلس الانتقالي.
تنبيه: الأشعة فوق البنفسجية ضارة بالجلد والعينين، استخدم حارسا مناسبا عند تصور الجل وارتداء معدات الوقاية الشخصية المناسبة.
ملاحظة: قد تكون مستعمرة PCR غير ناجحة لبعض البكتيريا; المتابعة مع PCR من gDNA معزولة هو خيار بديل22.
- استئصال العصابات ~ 1.5 كيلو بايت باستخدام شفرة الحلاقة ونقل القطع هلام في أنابيب الطرد الدقيق نظيفة. المضي قدما في بروتوكول استخراج هلام وفقا لتعليمات الشركة المصنعة(جدول المواد). قياس تركيز الحمض النووي المنقى بواسطة مطياف ميكروفولوم.
ملاحظة: من المستحسن التركيز >10 نانوغرام/ميكرولتر، و A260/280 بين 1.7-2.0 مقبول.
- إعداد اثنين من ردود الفعل تسلسل سانجر لكل عينة، واحدة باستخدام 8F والآخر باستخدام التمهيدي 1492R في المياه الخالية من النيوكليز وفقا للمبادئ التوجيهية لأي خدمة تسلسل سانجر المختار.
- بمجرد تلقي بيانات التسلسل ، قم بتحميل تسلسل الحمض النووي إلى موقع NCBI Basic Alignment Search Tool (BLAST) (blast.ncbi.nlm.nih.gov/Blast.cgi) ، واختر Nucleotide BLAST (blastn) ، وحدد قاعدة بيانات rRNA / ITS 16S تسلسل الحمض النووي الريبي الريبوسومي (البكتيريا والأرشايا) ، وتشغيل برنامج Megablast. يمكن التعرف على العزل بواسطة أعلى جودة ضرب إلى مرجع من قاعدة البيانات.
ملاحظة: تظهر بعض الأنواع البكتيرية هوية عالية في تسلسلات rRNA 16S الخاصة بها وقد لا يمكن تمييزها عن طريق هذه الطريقة وحدها. سوف يتطلب التحليل تحليل الحمض النووي وتحليلات البيوكيميائية للتمييز بثقة أعضاء من نفس الجنس23.
3. استخراج الحمض النووي الجينومي (gDNA)
ملاحظة: يستخدم هذا القسم الكواشف والأعمدة الدورانية المتوفرة في مجموعة الاستخراج gDNA المشار إليها في جدول المواد لاستخراج الغلة العالية للحمض النووي الجينومي عالي الجودة من الأنواع البكتيرية المتنوعة. فيما يلي التعديلات والتعليمات الموصى بها.
- إعداد الكواشف عدة لكل إرشادات الشركة المصنعة.
- إعداد 3-10 مل الثقافات في مرق العقيمة المناسبة (الشكل 1) عن طريق تلقيح البكتيريا من مستعمرات معزولة جيدا في وسائل الإعلام واحتضان في درجة الحرارة والضغط الجوي المشار إليها في الشكل 1 حتى لوحظ نمو كاف.
- بعد الحضانة، قياس الكثافة البصرية في 600 نانومتر (OD600)من الثقافة باستخدام مطياف24.
- إعداد العينة للقياس الكمي عن طريق تمييع الثقافات بين عشية وضحاها في نسبة 1:10. قم بتضمين فراغ من الوسائط الثقافية العقيمة للقياس أيضا. حساب الكثافة البصرية عن طريق طرح القراءة الفارغة من قراءة العينة والضرب بعامل التخفيف من عشرة.
- باستخدام قياس OD600 ونسبة OD600 إلى CFU/mL المحددة مسبقا للأنواع ، احسب عدد ملليلترات الاستزراع اللازمة للحصول على 2 × 109 خلايا.
- الطرد المركزي حجم الثقافة المطلوبة لمدة 5 دقائق في 5000 × ز إلى بيليه. يستنشق supernatant وإعادة إنفاق بيليه في 200 ميكرولتر العازلة TE الباردة (قبل البرد على الجليد في بداية الإجراء).
- الطرد المركزي العينة لمدة 2 دقيقة في 5000 x g. إزالة supernatant، ومن ثم resuspend بيليه في 180 ميكرولتر من العازلة تحلل الأنزيمية (ELB) وإضافة 20 ميكرولتر من RNase A المسلوقة مسبقا (10 ملغم / مل). للتحلل الفعال للبكتيريا إيجابية الغرام، إضافة 18 ميكرولتر من المتانولين (25 كيلوU/mL). دوامة جيدا، ومن ثم احتضان العينات في 37 درجة مئوية على الدوار لمدة 2 ساعة.
ملاحظة: من المستحسن الاستفادة من ELB الموصوفة في بروتوكول الشركة المصنعة لكل من البكتيريا إيجابية الغرام وغرام سلبية.
- المضي قدما وفقا لتعليمات الشركة المصنعة.
ملاحظة: كرر الخطوات elution مرة أو مرتين للحصول على العائد gDNA إضافية إذا رغبت في ذلك.
- تقييم جودة gDNA المستخرجة وفقا للتعليمات الواردة في القسم 4 وتخزين gDNA في 4 درجة مئوية إذا كان سيتم استخدامه في غضون أسبوع واحد. بدلا من ذلك، احتفظ ب gDNA عند -20 درجة مئوية للتخزين على المدى الطويل.
4. تقييم جودة gDNA المستخرجة
- لتقييم الجودة بواسطة هلام الكهربائي، وإعداد هلام agarose 1٪ كما هو موضح في القسم الفرعي 2.2. إعداد العينة في أنبوب نظيف: مزيج 1-2 ميكرولتر من gDNA المستخرجة و 3 ميكرولتر من 2x تحميل صبغة على parafilm. تشغيل هلام تحميل مرة واحدة، ومن ثم تصور تحت ضوء الأشعة فوق البنفسجية.
ملاحظة: استخراج gDNA ناجحة سوف يكون واضحا من قبل الفرقة منفصلة في الجزء العلوي من هلام والحد الأدنى من تشويه (الشكل 2A). التلطيخ يدل على القص. إذا لم يكن هناك نطاق gDNA واضح و / أو تلطيخ كبيرة، كرر استخراج gDNA. فكر في تقليل أوقات الحضانة في RNase A وبروتيناز K. إذا لوحظت نطاقات اثنين حول 1.5-3 كيلوبايت، وهذا يشير إلى تلوث الحمض النووي الريبي(الشكل 2B). إعداد RNase A الطازجة وتكرار الاستخراج.
- لتقييم الجودة بواسطة مطياف الميكروفولومي، قم بقياس تركيز gDNA ونسبة الامتصاص A260/280 بواسطة مطياف ميكروفولومي. التركيزات >50 نانوغرام/ميكرولتر و A260/280 بين 1.7-2.0 مقبولة.
ملاحظة: قد يكون انخفاض غلة gDNA بسبب انخفاض المدخلات، وارتفاع المدخلات، وتلوث النيوكلياس، وتحلل غير كاف. نسب الامتصاص فوق النطاق تشير إلى تلوث الحمض النووي الريبي. كرر استخراج إذا كانت نوعية gDNA سيئة.
- لتقييم الجودة حسب مقياس الفلور، اتبع تعليمات الشركة المصنعة لتحديد تركيز gDNA باستخدام عدة المقايسة عالية الحساسية وأداة قياس الفلورومتر(جدول المواد). التركيز >50 نانوغرام / ميكرولتر أمر مرغوب فيه.
5. الاقتران نهاية الجيل القادم قصيرة القراءة التسلسل وإعداد المكتبة
ملاحظة: يمكن تنفيذ تسلسل القراءة القصيرة على أدوات مختلفة بأطوال واتجاهات قراءة متميزة. يوصى بتسلسل ثنائي الطرف 150 نقطة أساس (300 دورة) ل WGS البكتيرية. ويمكن الاستعانة بمصادر خارجية لإعداد المكتبات وتسلسلها إلى المرافق الأساسية أو المختبرات التجارية.
- إعداد مكتبة التسلسل وفقا لتعليمات الشركة المصنعة(جدول المواد). اتبع تركيز مكتبة التحميل النهائي الموصى به من قبل الشركة المصنعة؛ ومع ذلك، تعديل الموصى به لتحميل المكتبة المجمعة في 1.8 pM لإنشاء القراءة الأمثل على أدوات NextSeq.
- على الرغم من اختياري، استخدم Bioanalyzer(جدول المواد)لتقييم توزيع جزء المكتبة المجمعة وضمان أن حجم الجزء هو 600 نقطة أساس في المتوسط.
6. نانوبور MinION تسلسل إعداد المكتبة
- إعداد مكتبة التسلسل وفقا لبروتوكول الشركة المصنعة(جدول المواد). باستخدام مجموعتين توسيع الباركود يسمح لتعدد ما يصل إلى 24 عينات على خلية تدفق واحدة. من المستحسن إجراء إعداد المكتبة في جزأين، 12 عينة في وقت تعدد 24 عينات. ويمكن تجميع جميع العينات ال 24 على النحو المبين أدناه.
ملاحظة: قد يتم تخزين عينات عند 4 °C بين عشية وضحاها عند الانتهاء من ربط الباركود الأصلي - وهذا يوفر نقطة توقف في البروتوكول، إذا لزم الأمر. في نهاية قسم ربط الباركود الأصلي من بروتوكول إعداد المكتبة ، يوصى بتجميع كميات متساوية من كل عينة تصل إلى الحد الأقصى لكتلة الحمض النووي (ng) الممكنة.
- للقيام بذلك، قم بتحديد جميع العينات التالية ربط الباركود باستخدام مقياس الفلورومتر(جدول المواد)لكل إرشادات الشركة المصنعة. تقدير حجم العينة مع تركيز dsDNA أدنى، ومن ثم حساب dsDNA الإجمالي وجدت في هذه العينة. استخدم هذا الرقم لتحديد كميات متساوية من كافة العينات الأخرى التي سيتم تجميعها معا.
ملاحظة: لأن حساب equimolar سوف تكبير مقدار dsDNA مجمعة وبالتالي ينتج تجمع عالي الحجم (>65 ميكرولتر) ، تنظيف ضرورية لتركيز التجمع.
- تنظيف تجمع dsDNA وتركيزه
- إضافة حجم 2.5x من الخرز paramagnetic(جدول المواد)إلى تجمع الحمض النووي، ومن ثم نفض الغبار بلطف أنبوب لخلط محتويات. ضع الأنبوب في الدوار لمدة 5 دقائق في RT. سبين أسفل العينة في 2000 × ز والكريه على المغناطيس.
- إضافة 250 ميكرولتر أعدت حديثا 70٪ الإيثانول (في ماء خال من النيوكليز)، مع الحرص على عدم إزعاج بيليه. يستنشق الإيثانول وكرر غسل الإيثانول مرة واحدة.
- بعد الطموح الثاني، تدور أسفل العينة في 2000 x ز ووضعها مرة أخرى على المغناطيس. ماصة قبالة أي الإيثانول المتبقية والسماح للعينة لتجف لمدة 30 ق تقريبا.
- إزالة أنبوب من المغناطيس وإعادة إنفاق بيليه في 60-70 ميكرولتر من الماء الخالي من النيوكليز. احتضان في RT لمدة 2 دقيقة. بيليه العينة على المغناطيس حتى elute واضحة، ومن ثم إزالة elute ونقلها إلى أنبوب الطرد المركزي الدقيق 1.5 مل نظيفة.
- تحديد حجم تجمع مركزة باستخدام مقياس الفلور، ومن ثم إعداد aliquot للمضي قدما في خطوة ربط محول: إعداد 700 نانوغرام من العينة في حجم نهائي 65 ميكرولتر. احتفظ بباقي التجمع عند 4 درجات مئوية للركض ثانية ليتم إكمالها بمجرد انتهاء التشغيل الأول.
- تابع ربط المحول كما هو من قبل الشركة المصنعة ثم قم بتحميل العينة على الخلية التدفق. بدء تشغيل التسلسل.
ملاحظة: أسبيرات الهواء و ~ 200 ميكرولتر من المخزن المؤقت التخزين من منفذ فتيلة الخلية تدفق قبل تحميل العينة. وهذا أمر بالغ الأهمية لنجاح تدفق فتيلة الخلية وتحميل العينة. استخدام ماصة P1000 ونصائح عند رسم وإيداع الحلول من خلال منفذ فتيلة من الخلية تدفق.
- تسلسل المكتبة وفقا لتعليمات الشركة المصنعة.
- افتح برنامج التشغيل للتسلسل وانقر على ابدأ. إدخال اسم للتجربة، تسمية الموصى بها يتضمن تاريخ التشغيل واسم المستخدم. انقر على متابعة لاختيار كيت، حدد مجموعة مناسبة من إعداد المكتبة وحزمة (حزم) توسيع الباركود المستخدمة ، ثم انقر على متابعة تشغيل الخيارات.
- ضبط طول التشغيل إلى 48 ساعة إذا التخطيط لإعداد مكتبة كافية لتشغيل ثانية (وإلا ترك في الافتراضي 72 ساعة). انقر على متابعة إلى Basecalling.
- تحقق من الخيار basecalling التكوين: Basecalling سريع وتأكد من تعيين باركودينج إلى تمكين بحيث سيتم اقتطاع إخراج FASTQ الملفات من تسلسل الباركود و demultiplexed إلى دلائل منفصلة استنادا إلى الباركود. انقر على متابعة الإخراج.
- اختر مكان حفظ بيانات تسلسل الإخراج. توقع ما يقرب من 30-50 غيغابايت من البيانات إذا حفظ إخراج FASTQ فقط و >500 غيغابايت من البيانات إذا أيضا حفظ إخراج FAST5. إلغاء تحديد الخيار تصفية Qscore: 7 | Readlength: غير المصفاة إذا كان التخطيط للمضي قدما في التصفية الموصوفة في القسم 7.2، وإلا ترك فحص وتعديل Readlength إلى 200.
- انقر على متابعة تشغيل الإعداد ومراجعة كافة الإعدادات. إذا كانت الإعدادات صحيحة، انقر على ابدأ،وإلا انقر على العودة وإجراء أي تعديلات ضرورية.
- إذا رغبت في ذلك، قد يتم غسل خلية التدفق وفقا لتعليمات الشركة المصنعة وإعادة تحميلها مع التجمع المتبقي. كرر الخطوات في 6.2 للتجمع المتبقية بمجرد اكتمال التشغيل الأول ثم تم غسل الخلية تدفق.
ملاحظة: عند إعداد التشغيل الثاني، ضبط الجهد التحيز إلى -250 mV لكل توصيات الشركة المصنعة لتدفق الخلايا المستخدمة سابقا في أشواط أكثر من 48 ساعة.
7. تقييم وإعداد القراءات
ملاحظة: يتم تصوير بنية دليل مستحسن في الشكل 4. إنشاء الدلائل الموجودة في سطح المكتب، وهي Long_Reads Short_Reads Trimmed_Reads ، قبل المتابعة مع خطوات الحساب أدناه.
- قراءات قصيرة (الشكل 3)
ملاحظة: يتم إنشاء قراءات قصيرة بتنسيق FASTQ. تحتوي الملفات على 4000 القراءة القصوى لكل FASTQ. وغالبا ما تكون مضغوط هذه (أرشيف .gz) وتنظيمها في ملفات متعددة. اعتمادا على النظام الأساسي، يتم عادة اقتطاع الباركود. بعض البرامج تقبل الملفات في شكل مضغوط ، والبعض الآخر قد تتطلب استخراجها قبل الاستيراد. يجب أن تمر القراءات خطوات مراقبة الجودة (QC) لضمان دقة البيانات أثناء تجميع الجينوم. إذا لم تتوفر منضدة عمل الجينوم CLC، فقد يتم استخدام برامج بديلة للتقليم وقرأات QC القصيرة مثل Trimmomatic25 أو تريم Galore (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) لتقليم وFastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) لتقييم جودة القراءة. ويوصى بأن يكون متوسط تغطية القراءة القصيرة، المقدرة بضرب عدد القراءات بمتوسط طول القراءة وقسمتها على حجم الجينوم، 100x >.- فتح الجينوم برنامج منضدة العمل(جدول المواد)واستيراد جميع الملفات FASTQ المقترنة نهاية قصيرة القراءة. سيتم إنشاء الملفات المقترنة تلقائيا.
- إنشاء مجلد جديد تحت CLC_Data بالنقر فوق جديد في شريط الأدوات العلوي وتحديد المجلد... لتخزين الملفات. اسم المجلد كما هو مطلوب، اصطلاح مستحسن يستخدم نموذج معرف. حفظ كافة الإخراج من الخطوات التالية إلى هذا المجلد.
- في شريط الأدوات العلوي، انقر على زر الاستيراد وحدد Illumina... انتقل إلى كافة الملفات قصيرة القراءة التي تتوافق مع العينة وحددها. تأكد من تحديد خيار القراءات المقترنة ثم قم بإلغاء تحديد الخيار إزالة القراءات الفاشلة. انقر على التالي، حدد حفظ، وانقر على التالي مرة أخرى. اختر حفظ الملفات المستوردة في المجلد الجديد الذي تم إنشاؤه في الخطوة السابقة وانقر فوق إنهاء.
- إنشاء قائمة تسلسل كافة الملفات المقترنة لعزل; سيؤدي ذلك إلى تجميع بيانات القراءة في ملف واحد لبساطة التحليل.
- في شريط الأدوات العلوي، انقر فوق الزر جديد وحدد قائمة التسلسل... في قائمة الدليل على اليسار، حدد الملفات المراد تسلسلها واستخدم الأسهم لنقلها إلى قائمة الملفات المحددة على اليمين. انقر على التالي، حدد حفظ، وانقر على التالي مرة أخرى. اختر حفظ قائمة التسلسل وانقر على إنهاء.
- بمجرد إنشاء قائمة التسلسل، قم بإعادة تسميتها فورا باستخدام نموذج معرف.
- تشغيل أداة QC للتسلسل القراءات على قائمة التسلسل: سيقوم هذا الإجراء بتقييم معلمات الجودة الكلية للقراءات التي تم إنشاؤها بواسطة NGS قصيرة القراءة.
- البحث عن أداة QC ل "قراءة التسلسل" في قائمة مربع الأدوات (الإطار السفلي الأيسر). انقر نقرا مزدوجا فوق الأداة، ثم اختر قائمة التسلسل لتحليلها وانقر على التالي.
- تأكد من فحص جميع خيارات الإخراج واختر حفظ ضمن معالجة النتائج. انقر فوق التالي وحدد لحفظ ملفات الإخراج، ثم انقر فوق إنهاء.
- تشغيل أداة القراءة تريم على قائمة التسلسل: سيتم التشذيب استنادا إلى الجودة والطول والغموض. تفترض هذه العملية أن الباركود المستخدم في التسلسل قد تم اقتطاعه قبل هذه الخطوة.
- البحث عن أداة "القراءات تريم" في مربع الأدوات (الإطار السفلي الأيمن). انقر نقرا مزدوجا على تريم يقرأ، ثم اختر قائمة التسلسل ليتم تحليلها وانقر على التالي.
- التشذيب الجودة: تعيين حد درجة الجودة إلى 0.01 وترك النيوكليوتيدات غامضة في 2. انقر على التالي.
ملاحظة: يمكن تعديل المعلمات حسب تقدير المستخدم؛ هذه هي الإعدادات الموصى بها.
- إلغاء تحديد التلقائي القراءة من خلال محول التشذيب (فقط القيام بذلك إذا تم اقتطاع المحولات من القراءات قبل استيرادها إلى CLC). انقر على التالي وتحقق تجاهل القراءات أدناه طول، استخدم الافتراضي 15.
- انقر على التالي، حدد إنشاء تقرير، ثم اختر حفظ. انقر على التالي وحدد مكان حفظ ملفات الإخراج. انقر على إنهاء.
- تصدير قائمة التسلسل المشذب: سيتم إكمال التجميع والتحليل المختلط اللاحق خارج CLC ويتطلب تصدير الملفات القصيرة القراءة.
- من التنقل الدليل في أعلى اليسار، اختر الملف المشذب الذي تم إنشاؤه في الخطوة 7.1.4، ثم انقر فوق تصدير في شريط الأدوات العلوي. حدد Fastq لنوع ملف التصدير وانقر على التالي. تحقق تصدير قائمة تسلسل المقترنة إلى ملفين. ثم انقر فوق التالي واختر الدليل Trimmed_Reads لتصدير الملفات إليه. انقر على إنهاء. تأكد من أن الملفات القصيرة القراءة التي تم اقتطاعها تم تصديرها بنجاح كملفين (R1 و R2) مع الملحق .fastq.
ملاحظة: يجب تصدير قائمة التسلسل المشذبة إلى ملفين، عادة ما يتم تعيينها بواسطة CLC ك R1 و R2. وهذا أمر بالغ الأهمية لأن التجميع المختلط المتلقين للمعلومات يتطلب إدخال بيانات للقراءة القصيرة لإعدادها على هذا النحو.
- إعادة تسمية الملفات المصدرة، يرجى الامتناع عن استخدام المسافات والأحرف الخاصة في أسماء الملفات. للبساطة تنسيق الموصى بها هو trimmed_short_file. R1.fastq.
- طويل (MinION) يقرأ (الشكل 3)
ملاحظة: خط الأنابيب التالي لإعداد طويل (MinION) تسلسل يقرأ التجميع المختلط يستخدم NanoFilt و Nanostat البرامج26 تنفيذها بواسطة سطر الأوامر. تثبيت الأدوات قبل المتابعة وتكون على دراية أساسيات UNIX من أجل تنفيذ هذه الأوامر. يوصى بالمحطات الطرفية الافتراضية و Bash Shell. تم العثور على دليل الدرس للأوامر الطرفية المشتركة والاستخدام في النجارة البرمجيات27. تفترض الإرشادات أدناه أن الملفات التي تم إنشاؤها سيتم تسميتها بتسمية الباركود (NB01 و NB02 وما إلى ذلك) ويتم حفظها في دليل Long_Reads. بدلا من ذلك، يمكن إنجاز تصفية القراءة باستخدام MinKNOW عند إعداد تشغيل التسلسل. يوصى بأن يكون متوسط تغطية القراءة الطويلة 100x >. متوسط طول القراءة الموصى به هو >2000 نقطة أساس؛ لذلك، عدد القراءات الطويلة المطلوبة أقل من عدد القراءات القصيرة.
- إنشاء أدلة جديدة لكل الباركود المستخدمة في المدى (الباركود01، الباركود02، الخ) داخل الدليل Long_Reads(الشكل 4). نسخ كافة الملفات .fastq التي تتوافق مع كل الباركود في المجلد المناسب. الجمع بين جميع الملفات .fastq لكل الباركود من كل شوط.
- فتح المحطة الطرفية والانتقال إلى دلائل الباركود داخل الدليل Long_Reads باستخدام الأمر cd: cd Desktop/Long_Reads/barcode01
- تجميع كافة ملفات .fastq لكل الباركود في ملف .fastq واحد بتنفيذ الأمر التالي: cat *.fastq > NB01.fastq
ملاحظة: يجمع هذا الأمر كافة القراءات من كل من الملفات FASTQ في واحد FASTQ كبيرة واحدة باسم NB01.fastq.
- استخدم NanoStat لتقييم جودة قراءة العينة بتنفيذ الأمر التالي: NanoStat --fastq NB01.fastq
- تسجيل النتائج عن طريق نسخ الإخراج إلى نص أو ملف Word للرجوع إليها في المستقبل.
- استخدام NanoFilt لتصفية MinION يقرأ تجاهل يقرأ مع Q < 7 وطول < 200 عن طريق تنفيذ الأمر: NanoFilt -q 7 -l 200 bp NB01.fastq | gzip > NB01 _trimmed.fastq.gz
- تشغيل NanoStat على الملف قلصت التي تم إنشاؤها في الخطوة 7.2.6 بتنفيذ الأمر: NanoStat --fastq NB01 _trimmed.fastq.gz
- تسجيل النتائج عن طريق نسخ الإخراج إلى نص أو ملف Word ومقارنة النتائج من الخطوة 7.2.4 للتأكد من أن التصفية كانت ناجحة(الجدول 1).
- كرر الخطوات 7.2.2 إلى 7.2.8 لكل الباركود المستخدمة في تشغيل التسلسل.
ملاحظة: سيتم استخدام الملف NB01_trimmed.fastq.gz الذي تم إنشاؤه في الخطوة 7.2.6 للتجميع المختلط.
8. توليد تجميع الجينوم الهجين
ملاحظة: يستخدم خط أنابيب التجميع التالي Unicycler19و28و29و30 للجمع بين القراءات القصيرة والطويلة المعدة في القسمين 7.1 و 7.2 ( الشكل3). تثبيت Unicycler وتبعياتها وتنفيذ الأوامر أدناه. يفترض أن يتم تسمية الملفات قصيرة القراءة التي تم تصديرها في الخطوة 7.1.5 trimmed_short_file. R1.fastq trimmed_short_file. R2.fastq للبساطة.
- تنظيم الملفات قصيرة القراءة والملفات طويلة القراءة في دليل واحد اسمه Trimmed_Reads. يجب أن يحتوي الدليل على ما يلي:
- ملف .gz.fastq لقراءات طويلة قلصت (التي تم إنشاؤها في الخطوة 7.2.6).
- ملفين .fastq (R1 و R2) لقراءات قصيرة قلصت (التي تم إنشاؤها في الخطوة 7.1.5).
- انتقل إلى Trimmed_Reads الدليل الذي يخزن الملفات المقروءة باستخدام الأمر cd في Terminal: cd Desktop/Trimmed_Reads
- مرة واحدة في الدليل الصحيح، الرمز البريدي اثنين من ملفات القراءة القصيرة بحيث تكون أيضا في تنسيق .gz.fastq عن طريق تنفيذ الأمر التالي: gzip trimmed_short_file. R1.fastq
- كرر الخطوة 8.2 لكل من R1 و R2. تحقق من أن كافة الملفات المقروءة الآن بتنسيق .fastq.gz وتحقق من تطابق كافة الملفات مع نفس العزل.
- بدء التجميع المختلط باستخدام Unicycler عن طريق تشغيل الأمر التالي:
unicycler -1 trimmed_short_file. R1.fastq.gz -2 trimmed_short_file. R2.fastq.gz -l NB01 _trimmed.fastq.gz -o unicycler_output_directory
ملاحظة: -o يحدد الدليل الذي سيتم حفظ الإخراج Unicycler، سيتم إنشاء Unicycler هذا الدليل بمجرد تنفيذ الأمر; لا تقم بإنشاء الدليل مسبقا. وقت التشغيل يختلف حسب القوة الحسابية للكمبيوتر المستخدم وكذلك حجم الجينوم وعدد القراءات. قد يستغرق هذا في أي مكان من 4 ساعة إلى 1 أو 2 أيام. تم تنفيذ هذا البروتوكول على جهاز CentOS Linux 7 مع ذاكرة وصول عشوائي 250 غيغابايت، وحدة المعالجة المركزية Intel Xeon (R) مع 2.5 جيجاهرتز 12 نواة عملية و48 نواة افتراضية. بدلا من ذلك، يمكن لأجهزة الكمبيوتر الشخصية المزودة بذاكرة وصول عشوائي 16 غيغابايت ومعالجات 6-core 2.6 جيجاهرتز حساب هذه التجميعات في وقت معالجة أطول.
- عند اكتمال التشغيل، راجع ملف unicycler.log للتأكد من عدم وجود أخطاء - سجل رقم المتجاورات التي تم إنشاؤها وحجمها وحالتها (كاملة وغير مكتملة).
- إذا تم تعريف المتجاورات غير كاملة (المشار إليها على أنها غير كاملة في سجل Unicycler) ، إعادة تشغيل Unicycler في وضع غامق بإضافة العلامة التالية إلى الأمر في الخطوة 8.4: --وضع غامق.
ملاحظة: وضع غامق سيخفض عتبة الجودة المقبولة الجسور القراءة الطويلة أثناء التجميع; قد ينتج تجميع كامل ولكن قد يتم تقليل جودة التجميع. فمن المستحسن للاستفادة من وضع جريئة فقط عند الضرورة وكأدلة أولية للانضمام إلى الملتق ليتم تأكيده في وقت لاحق من قبل PCR.
9. تقييم جودة التجميع
ملاحظة: يستخدم البروتوكول التالي Bandage31 و QUAST32، برنامجين يجب إعدادهما قبل الاستخدام(الشكل 2 والشكل 4). ضمادة لا تتطلب التثبيت مرة واحدة تحميلها وQUAST يتطلب الألفة مع استخدام سطر الأوامر الأساسية. ومن المستحسن أيضا لتقييم اكتمال الجينوم باستخدام القياس العالمي واحد نسخة تقويم العظام (BUSCO)33.
- ضمادة: انقر على ملف. ثم اختر تحميل Graph وحدد ملف assembly.gfa الذي تم حفظه إلى unicycler_output_directory التي تم إنشاؤها بواسطة Unicycler في الخطوة 8.4. بمجرد تحميلها، انقر على رسم الرسم البياني زر على شريط الأدوات الأيسر وننظر في كيفية اتصال المتجاورة (تسمى العقد) وتنظيمها لتقييم ما إذا كان التجميع هو كاملة (الشكل 5).
ملاحظة: يتم تمثيل التجميعات الكاملة بواسطة متجاورات دائرية مفردة مرتبطة في كلا الطرفين (الشكل 5A, B). التجميعات غير المكتملة لها عدة متجاورات مرتبطة ببعضها البعض أو خطية(الشكل 5C). قد لا تكون المتجاورات الخطية الصغيرة غير مكتملة لأنها قد تشير إلى عناصر خطية خارج نطاق الذاكرة. التغطية، وتسمى أيضا العمق، وسوف يلاحظ في ضمادة ويمثل وفرة النسبية من المتجاورة إلى الكروموسوم، تطبيع في Unicycler إلى 1x.
- كواست
- داخل المحطة الطرفية، انتقل إلى المجلد الذي يخزن إخراج Unicycler باستخدام الأمر cd: cd Desktop/Trimmed_Reads/unicycler_output_directory
ملاحظة: لا يسمح المسافات في المسار إلى حيث يوجد التجميع، أي، لا يمكن أن يكون أية دلائل تؤدي إلى إخراج Unicycler مسافات باسمها. بدلا من ذلك، نسخ ملف assembly.fasta إلى سطح المكتب لسهولة الوصول.
- تشغيل QUAST بتنفيذ الأمر التالي: quast assembly.fasta -o quast_output_directory
- مراجعة التقارير التي تم إنشاؤها بواسطة QUAST في دليل الإخراج quast_output_directory.
10. شرح الجينوم
ملاحظة: خط أنابيب التعليقات التوضيحية أدناه يستخدم Prokka34، أداة سطر الأوامر التي يجب تثبيتها قبل الاستخدام. بدلا من ذلك، استخدم بروكا من خلال الآلي واجهة المستخدم الرسومية K-قاعدة (جدول المواد) أو الجينوم المشروح عبر خادم الويب RAST35. إذا إيداع الجينوم في NCBI، سيتم المشروح تلقائيا باستخدام خط أنابيب الشرح الجينوم Prokaryotic (PGAP)36.
- انتقل داخل المحطة الطرفية إلى المجلد الذي يخزن الإخراج Unicycler باستخدام الأمر cd (راجع الخطوة 9.2.1). ثم قم بتشغيل Prokka بتنفيذ الأمر التالي: prokka --البادئة sample_ID --outdir prokka_output_directory assembly.fasta
ملاحظة: --بادئة سيتم تسمية كافة ملفات الإخراج استنادا إلى sample_ID المحدد.--سيقوم outdir بإنشاء دليل إخراج بالاسم المحدد حيث سيتم حفظ كافة ملفات إخراج Prokka; لا تقم بإنشاء دليل إخراج ل Prokka مسبقا.
- راجع التعليقات التوضيحية عن طريق فتح الجدول .tsv و/أو عن طريق تحميل الملف .gff الذي تم إنشاؤه في برنامج تحليل تسلسلي لتصور وتحليل التعليقات التوضيحية(الشكل 6).
- يمكن إنشاء أنواع محددة من الشروح اعتمادا على العوامل الوراثية ذات الاهتمام. من المستحسن أن تبدأ مع أدوات سهلة الاستخدام على مركز علم الأوبئة الجينومية (www.genomicepidemiology.org/) خادم الويب للتحليل الأولي37،38،39،40،41. تتوفر أدوات إضافية للكشف عن أنظمة CRISPR-cas والسند (الشكل 3)42،43.
11. الممارسات المقترحة لإضفاء الطابع الديمقراطي على البيانات
- عند الإمكان، قم بإيداع جميع بيانات القراءة الأولية وكذلك الجينوم المجمع في مستودع عام مثل أرشيف قراءة تسلسل NCBI (SRA) وGenbank. يتم تعليق الجينوم تلقائيا عبر خط أنابيب PGAP أثناء عملية ترسب NCBI.