المتفطرة البوفية خضعت سلالة BCG (عصيات كالميت وغيرين) 1173P2 التي تخضع لنمو أسي لسلسلة زمنية (0 و 4 و 10 و 20 يوما) من تجويع المغذيات ، تليها إنعاش لمدة 6 أيام في وسط غني بالمغذيات كما تم تقديمه سابقا في Hu et al.7. تم عزل الحمض النووي الريبي الصغير من الثقافة البكتيرية ، مع ثلاثة تكرارات بيولوجية ، في كل نقطة من النقاط الزمنية الخمس المحددة. تم إنشاء مكتبات Illumina باستخدام سير عمل إعداد مكتبة AQRNA-seq الموصوف أعلاه (الشكل 1) ، متبوعا بالتسلسل على جهاز التسلسل في مركز BioMicro التابع لمعهد ماساتشوستس للتكنولوجيا. ثم تمت معالجة بيانات التسلسل باستخدام خط أنابيب تحليلات بيانات AQRNA-seq (الشكل 2) المخصص للقياس الكمي لوفرة الحمض النووي الريبي.
بعد تضخيم PCR لمكتبة cDNA باستخدام بادئات التسلسل ، لوحظ وجود منتجات PCR بحجم 175 زوجا أساسيا (bp) في جميع العينات (الشكل 3 أ) ، مما يشير إلى تكوين دايمرات التمهيدي. للتخفيف من ترحيل الدايمرات التمهيدية ، تم استئصال منتجات تفاعل البوليميراز المتسلسل التي يتجاوز حجمها 195 نقطة أساس من الجل وتنقيتها (الشكل 3 ب).
تم تعيين قراءات التسلسل المصفاة والمشذبة بجودة الجودة إلى مكتبة تسلسل مرجعية مخصصة بما في ذلك 45 متقبلات أيزو الحمض النووي الريبي tRNA والمعيار الداخلي وتسلسلات التحكم (أي 23S rRNA و 16S rRNA و 5S rRNA و rnpB و ssr). شكلت مستقبلات أيزو الحمض النووي الريبي tRNA 10.5٪ إلى 40.2٪ من إجمالي القراءات المعينة لعينة معينة وأظهرت وفرة أعلى بكثير من تسلسلات التحكم (الشكل 4). الأهم من ذلك ، يمكن أن تعزى نسب القراءة المنخفضة نسبيا لمستقبلات أيزو الحمض النووي الريبي إلى الوفرة النسبية الأعلى للمعايير الداخلية. لذلك ، يمكن التحكم في نسب القراءة لمستقبلات أيزو الحمض النووي الريبي بالنسبة للمعايير الداخلية (الشكل 4 ، كتل اللون الوردي مقابل اللون الأخضر) بواسطة المشغل ، من خلال ضبط كمية المعيار الداخلي المتصاعد في التفاعل.
تمت تسوية بيانات وفرة tRNA الخام باستخدام طريقة متوسط النسب المطبقة مع إصدار حزمة DESeq2 (المشار إليه فيما يلي باسم v) 1.36.016 في بيئة البرمجة الإحصائية R (المشار إليها فيما يلي باسم R) v 4.2.117. بعد التطبيع ، يتم تحقيق مشهد كمي لمستقبلات أيزو الحمض النووي الريبي في المتفطرة البوفية BCG خلال فترة زمنية من تجويع المغذيات وإنعاشها (الشكل 5).
للكشف عن مجموعات متميزة من العينات ذات الأنماط الظاهرية المختلفة بناء على الأنماط في وفرة مستقبلات الحمض النووي الريبي tRNA ، تم إجراء تحليل المكون الرئيسي (PCA) على بيانات وفرة tRNA الطبيعية باستخدام حزمة الإحصائيات v 4.2.117 في R (الشكل 6). ميز التحليل عينات من يوم الجوع 0 ويوم الإنعاش 6 عن عينات أيام الجوع 4 و 10 و 20 ، مما يشير إلى اختلاف كبير في مشهد الحمض النووي الريبي tRNA للمتفطرة البوفية BCG المزروعة في وسط محروم من المغذيات ووسط غني بالمغذيات.
لتحديد ديناميكيات وفرة كل متقبل للحمض النووي الريبي عبر النقاط الزمنية الخمس المحددة ، تم إجراء تحليل التعبير التفاضلي على بيانات وفرة tRNA الطبيعية باستخدام حزمة DESeq2 v 1.36.0 في R (الشكل 7). وكشف التحليل أن 17 من أصل 20 عائلة من عائلات الأيزو مستقبلات تحتوي على متقبلات الأيزو التي تم التعبير عنها بشكل تفاضلي (أي منظمة بشكل كبير صعودا أو هبوطا) في نقطة زمنية واحدة على الأقل، مما يشير إلى دور محتمل لتنظيم تجمع الحمض النووي الريبي في الحالة المستمرة للمتفطرة البوفية BCG أثناء مرض السل.

الشكل 1: رسم تخطيطي لسير عمل إعداد مكتبة AQRNA-seq. يتم سرد الخطوات الرئيسية الموضحة في سير العمل في وسط المخطط ومتصلة بالرسوم التوضيحية الرسومية الخاصة بها بواسطة خطوط منقطة. يمكن العثور على الوصف التفصيلي لكل خطوة في قسم البروتوكول. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 2: رسم تخطيطي لخط أنابيب تحليلات بيانات AQRNA-seq. يتم سرد الخطوات الرئيسية الموضحة في خط الأنابيب في وسط المخطط ومتصلة بالرسوم التوضيحية الرسومية الخاصة بها بواسطة خطوط منقطة. يتوفر الوصف التفصيلي لكل خطوة على GitHub (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git). يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 3: الرحلان الكهربائي لهلام الأغاروز لشظايا cDNA بعد تضخيم تفاعل البوليميراز المتسلسل باستخدام بادئات التسلسل. أ: صورة الجل قبل استخلاص الجل وتنقيته. يحتوي الممران 7 و 14 من الجانب الأيسر على 5 ميكرولتر من سلم الحمض النووي 50 bp ، بينما تحتوي الممرات الأخرى على 20 ميكرولتر من كل عينة من العينات ال 15. يشير توطين حجم منتجات تفاعل البوليميراز المتسلسل إلى أعلى تركيز لها في النطاق من 175 نقطة أساس (دايمرات التمهيدي) إلى 300 نقطة أساس (بادئان + 120 نقطة أساس 5S rRNA). ب: صورة الجل بعد استخراج الجل وتنقيته. لكل عينة ، تم استئصال كتلة الهلام بين 200 نقطة أساس و 400 نقطة أساس لتقليل تلوث دايمرات التمهيدي في مكتبة التسلسل. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 4: عدد قراءات التسلسل التي تم تعيينها بنجاح إلى مكتبة التسلسل المرجعي. يظهر المحور السيني أسماء العينات (على سبيل المثال ، D18-69XX) مجمعة حسب النقطة الزمنية (على سبيل المثال ، يوم الجوع 0). لكل عينة ، يتم تمثيل عدد القراءة المرتبط بفئات الموضوعات المستهدفة المختلفة باستخدام كتل ملونة مكدسة فوق بعضها البعض. تمثل الأرقام الموجودة في وسط كتل الألوان نسب القراءات المقابلة للموضوعات المستهدفة المعنية داخل عينة معينة. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 5: المشهد الكمي لمستقبلات الحمض النووي الريبي الثلاثي ل المتفطرة البوفية BCG في نقاط زمنية مختلفة على طول مسار وقت الجوع والإنعاش. تم تسوية بيانات وفرة الحمض النووي الريبي الخام باستخدام طريقة متوسط النسب. هنا ، يصور كل صف وفرة tRNA الطبيعية (المحور y) كمتوسط ± خطأ قياسي ل 3 مكررات بيولوجية في كل نقطة زمنية. على المحور السيني ، تم تجميع الأيزومستقبلات من نفس العائلة معا ووسمها بالحمض الأميني المقابل. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 6: مخطط جيب التمام التربيعي للعينات المشتقة من تحليل المكون الرئيسي (PCA). تم إجراء PCA على أساس وفرة tRNA الطبيعية. يشير جيب التمام التربيعي إلى أهمية المكونات الرئيسية للعينات، وقد رسمت العينات بالنسبة إلى جيب التمام التربيعي للمكونين الرئيسيين الأولين. تم تصنيف العينات باستخدام معرفات العينة وتم ترميزها بالألوان حسب النقطة الزمنية. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.

الشكل 7: التعبير التفاضلي لمستقبلات الحمض النووي الريبوزي الناقل (tRNA) عبر نقاط زمنية مختلفة. تم تلخيص وفرة tRNA الطبيعية على أنها وسائل (عقدة خطية) ± خطأ معياري (أشرطة خطأ) عبر 3 مكررات بيولوجية. بسبب ضيق المساحة ، تم اختصار الشروط على النحو التالي: S0-S20 = أيام المجاعة 0-20 ؛ R6 = يوم الإنعاش 6. تم إجراء تحليل التعبير التفاضلي لكل مستقبل إيزومتبل للحمض النووي الريبي tRNA ، ومقارنة النقاط الزمنية المختلفة بطريقة زوجية باستخدام اختبار نسبة الاحتمال واختبار والد. تم استخدام الحروف المدمجة لتمثيل الدلالة الإحصائية ، حيث لم تكن وفرة متقبل إيزومت الحمض النووي الريبي tRNA معين في نقاط زمنية تشترك في حرف مشترك واحد على الأقل مختلفة بشكل كبير عن بعضها البعض. على سبيل المثال ، تم تنظيم وفرة tRNA-Lys-CTT-1-1 (في لوحة Lysine) بشكل كبير من S0 إلى S4 ومن S4 إلى S10 ، ولكن ليس من S10 إلى S20. ثم تم تنظيمه بشكل كبير من S20 إلى R6. يرجى النقر هنا لعرض نسخة أكبر من هذا الرقم.
الجدول 1: قليل النيوكليوتيدات المشاركة في سير عمل إعداد مكتبة AQRNA-seq. المعيار الداخلي هو الحمض النووي الريبي ، في حين أن جميع قليل النيوكليوتيدات الأخرى هي الحمض النووي. تعتبر بادئات PCR وبادئات التسلسل المخصصة المدرجة خاصة بمنصات التسلسل. يمكن تصميم بادئات PCR إضافية بتسلسلات فهرس جديدة. الرجاء الضغط هنا لتنزيل هذا الجدول.