$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Mikobakteri bovis Üstel büyüme geçiren BCG (bacilli de Calmette et Guérin) suşu 1173P2, bir zaman serisine (0, 4, 10 ve 20 gün) besin açlığına tabi tutuldu, ardından daha önce Hu ve ark.7'de sunulduğu gibi besin açısından zengin ortamda 6 günlük bir resüsitasyon izledi. Küçük RNA'lar, belirlenen beş zaman noktasının her birinde üç biyolojik kopya ile bakteri kültüründen izole edildi. Illumina kütüphaneleri, yukarıda açıklanan AQRNA-seq kütüphane hazırlama iş akışı (Şekil 1) kullanılarak oluşturuldu ve ardından Massachusetts Institute of Technology'nin BioMicro Center'daki bir sıralayıcı üzerinde dizileme yapıldı. Dizileme verileri daha sonra tRNA bolluğunun nicelleştirilmesi için özelleştirilmiş AQRNA-seq veri analitiği boru hattı (Şekil 2) kullanılarak işlendi.
cDNA kütüphanesinin dizileme primerleri ile PCR amplifikasyonundan sonra, tüm numunelerde 175 baz çifti (bp) boyutunda PCR ürünlerinin varlığı gözlendi (Şekil 3A), bu da primer dimerlerinin oluşumunu düşündürdü. Primer dimerlerin taşınmasını azaltmak için, boyutu 195 bp'yi aşan PCR ürünleri jelden çıkarıldı ve saflaştırıldı (Şekil 3B).
Kalite filtreli ve kırpılmış dizi okumaları, 45 tRNA izoalıcıları, dahili standart ve kontrol dizileri (yani, 23S rRNA, 16S rRNA, 5S rRNA, rnpB ve ssr) dahil olmak üzere özel bir referans dizi kitaplığına eşlendi. tRNA izoseptörleri, belirli bir numunenin toplam eşlenen okumalarının .5 ila @.2'sini oluşturdu ve kontrol dizilerinden çok daha yüksek bir bolluk gösterdi (Şekil 4). Daha da önemlisi, tRNA izoseptörlerinin nispeten düşük okuma oranları, iç standartların daha yüksek nispi bolluğuna bağlanabilir. Bu nedenle, tRNA izoeseptörlerinin iç standartlara göre okuma oranları (Şekil 4, pembe ve yeşil renk blokları), reaksiyona giren iç standart miktarında ince ayar yapılarak operatör tarafından kontrol edilebilir.
Ham tRNA bolluk verileri, R İstatistiksel Programlama Ortamında (bundan böyle R olarak anılacaktır) 1.36.016 paket versiyonu (bundan böyle v olarak anılacaktır) v 4.2.117 ile uygulanan medyan oran yöntemi kullanılarak normalize edildi. Normalleştirmeden sonra, besin açlığı ve resüsitasyon süresi boyunca Mycobacterium bovis BCG'deki tRNA izoseptörlerinin kantitatif bir manzarası elde edilir (Şekil 5).
tRNA izoseptör bolluğundaki modellere dayalı olarak farklı fenotiplere sahip farklı numune kümelerini ortaya çıkarmak için, R'de v 4.2.117 istatistik paketi kullanılarak normalleştirilmiş tRNA bolluğu verileri üzerinde Temel Bileşen Analizi (PCA) gerçekleştirildi (Şekil 6). Analiz, açlık günü 0 ve resüsitasyon günü 6 örneklerini açlık günü 4, 10 ve 20 örneklerinden ayırdı, bu da besin açısından yoksun ortamda ve besin açısından zengin ortamda yetiştirilen Mycobacterium bovis BCG'nin tRNA manzarasında önemli bir fark olduğunu düşündürdü.
Belirlenen beş zaman noktası boyunca her bir tRNA izoeperatörünün bolluğunun dinamiklerini profillemek için, R'de DESeq2 paketi v 1.36.0 kullanılarak normalleştirilmiş tRNA bolluk verileri üzerinde diferansiyel ekspresyon analizi yapıldı (Şekil 7). Analiz, 20 izoedüktör ailesinin 17'sinin, zaman noktalarından en az birinde diferansiyel olarak eksprese edilen (yani önemli ölçüde yukarı veya aşağı regüle edilmiş) izoseptörler içerdiğini ortaya koydu, bu da tüberküloz sırasında Mycobacterium bovis BCG'nin kalıcı durumunda tRNA havuzunun düzenlenmesinin potansiyel bir rolünü düşündürmektedir.

Şekil 1: AQRNA-seq kütüphanesi hazırlama iş akışının şeması. İş akışında özetlenen temel adımlar şemanın ortasında listelenir ve ilgili grafik çizimlerine noktalı çizgilerle bağlanır. Her adımın ayrıntılı açıklaması Protokol bölümünde bulunabilir. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 2: AQRNA-seq veri analizi boru hattının şeması. Boru hattında özetlenen temel adımlar şemanın ortasında listelenir ve noktalı çizgilerle ilgili grafik çizimlerine bağlanır. Her adımın ayrıntılı açıklaması GitHub'da (https://github.com/Chenrx9293/AQRNA-seq-JoVE.git) mevcuttur. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 3: Dizileme primerleri ile PCR amplifikasyonundan sonra cDNA fragmanlarının agaroz jel elektroforezi. (A) Jel ekstraksiyonu ve saflaştırmasından önce jelin görüntüsü. Sol taraftaki 7 ve 14 numaralı şeritler 50 bp DNA merdiveninin 5 μL'sini içerirken, diğer şeritler 15 örneğin her birinden 20 μL içerir. PCR ürünlerinin boyut lokalizasyonu, 175 bp (primer dimerler) ila 300 bp (iki primer + 120 bp 5S rRNA) aralığında en yüksek konsantrasyonlarını gösterir. (B) Jel ekstraksiyonu ve saflaştırılmasından sonra jelin görüntüsü. Her numune için, 200 bp ile 400 bp arasındaki jel blok, sıralama kütüphanesindeki primer dimerlerin kontaminasyonunu en aza indirmek için eksize edildi. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 4: Referans dizi kitaplığına başarıyla eşlenen dizi okumalarının sayısı. X ekseni, örneklerin adlarını (örneğin, D18-69XX) zaman noktasına göre gruplandırılmış olarak gösterir (örneğin, Açlık Günü 0). Her örnek için, çeşitli hedef konu kategorileriyle ilişkili okuma sayısı, üst üste yığılmış renk blokları kullanılarak temsil edilir. Renk bloklarının merkezinde yer alan sayılar, belirli bir örneklem içindeki ilgili hedef konulara karşılık gelen okuma oranlarını temsil eder. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 5: Mycobacterium bovis BCG'nin tRNA izoalıcılarının açlık ve resüsitasyon zaman seyri boyunca çeşitli zaman noktalarında kantitatif manzarası. Ham tRNA bolluk verileri, medyan oran yöntemi kullanılarak normalize edildi. Burada, her satır, normalleştirilmiş tRNA bolluklarını (y ekseni) her zaman noktasında 3 biyolojik kopya için ortalama ± standart hata olarak gösterir. X ekseninde, aynı aileden izoeperatörler birlikte gruplandırıldı ve karşılık gelen amino asit ile etiketlendi. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 6: Temel bileşen analizinden (PCA) türetilen örneklerin kare-kosinüs grafiği. PCA, normalleştirilmiş tRNA bolluğuna dayalı olarak gerçekleştirildi. Kare kosinüs, ana bileşenlerin numuneler için önemini gösterir ve numuneler, ilk iki ana bileşenin kare kosinüsüne göre çizilmiştir. Numuneler, numune kimlikleri kullanılarak etiketlendi ve zaman noktasına göre renk kodlandı. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Şekil 7: Farklı zaman noktalarında tRNA izoepertörlerinin diferansiyel ifadesi. Normalleştirilmiş tRNA bollukları, 3 biyolojik replikat boyunca ortalamalar (çizgi düğümleri) ± standart hata (hata çubukları) olarak özetlendi. Alan sınırlaması nedeniyle, koşullar şu şekilde kısaltıldı: S0-S20 = açlık günleri 0-20; R6 = resüsitasyon günü 6. Her bir tRNA izoeseptörü için, olasılık oranı testi ve Wald testi kullanılarak çeşitli zaman noktalarını ikili bir şekilde karşılaştırarak diferansiyel ekspresyon analizi yapıldı. Kompakt harfler, en az bir ortak harfi paylaşan zaman noktalarında belirli bir tRNA izoalıcısının bolluğunun birbirinden önemli ölçüde farklı olmadığı istatistiksel anlamlılığı temsil etmek için kullanıldı. Örneğin, tRNA-Lys-CTT-1-1'in (Lizin panelinde) bolluğu, S0'dan S4'e ve S4'ten S10'a önemli ölçüde aşağı regüle edildi, ancak S10'dan S20'ye değil. Daha sonra S20'den R6'ya önemli ölçüde yukarı regüle edildi. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Tablo 1: AQRNA-seq kütüphane hazırlama iş akışında yer alan oligonükleotidler. İç standart RNA'dır, diğer tüm oligonükleotidler ise DNA'dır. Listelenen PCR primerleri ve özel dizileme primerleri, dizileme platformlarına özgüdür. Yeni indeks dizileri ile ek PCR primerleri tasarlanabilir. Bu Tabloyu indirmek için lütfen buraya tıklayın.