$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Bu iş akışı, tek çekirdek çözünürlüğünde gen ekspresyon profilleri elde etmek için donmuş insan IMAT örneklerinin işlenmesine rehberlik etmek ve hücre tipi tanımlamasını mümkün kılmak üzere tasarlanmıştır. Burada, SOMMA çalışmasına katılan bir katılımcıdan alınan temsili bir IMAT örneği sunulmaktadır.
snRNA-seq verilerinin herhangi bir analizinin ilk adımı, potansiyel olarak veri setinden çıkarılması gereken düşük kaliteli çekirdekleri belirlemek için verilerin kalitesini değerlendirmektir. Daha da önemlisi, yaygın olarak değerlendirilen metrikler dokular ve hücre tipleri arasında farklılık gösterebileceğinden, elinizdeki belirli numune ve veri seti türü için filtreleme adımları ve eşikleri belirlenmelidir22,23. Şekil 4A, oluşturulan snRNA-seq verilerinin kalitesini değerlendirmek için kullanılan bazı temel metriklerin görsellerini sağlar. Çekirdek başına tespit edilen gen sayısı, dizileme derinliğine ve hücre tipine bağlıdır, ancak kaliteli çekirdekler için 200'ün üzerinde olması beklenir18,23. Bu protokol kullanılarak üretilen verilerin, toplam 4662 çekirdekten çekirdek başına ortalama 1134 gen ile beklenen aralıkta olduğu bulundu.
Mitokondriyal okumaların yüzdesi değerlendirilir, çünkü hasarlı çekirdeklerden veya çekirdeğe bağlanan ortam RNA'sından yüksek derecede mitokondriyal kontaminasyon ortaya çıkabilir, bu da düşük kaliteli çekirdekleri gösterir. Burada sunulan veri setinde, literatürde yaygın olarak kullanılan %5- eşiğinin oldukça altında olan 2.65'lik bir medyan mitokondriyal okuma yüzdesi bulunmuştur 24,25,26. Ribozomal okumaların yüzdesi hücre tipleri ve dokular arasında farklılık gösterir. Bununla birlikte, ribozomal genlerin büyük oranları verilerin kümelenmesini etkileyebileceğinden, kümelemeden önce ribozomal okuma yüzdesinin kontrol edilmesi ve potansiyel olarak yüksek düzeyde ribozomal gen içeren ribozomal genlerin veya çekirdeklerin veri kümesinden çıkarılması önerilir. Bu protokolle oluşturulan veriler, medyan %2,46 ve maksimum ,5 ile düşük düzeyde ribozomal okuma gösterdi ve bu nedenle bu metriğe göre filtreleme yapmadık. Son olarak, tespit edilen log (10) gen sayısının, tespit edilen log (10) okuma sayısına bölünmesine dayalı olarak bir hücre karmaşıklık skoru hesaplandı. İyi kaliteli çekirdeklerin 0.8'in üzerinde olması beklenir ve bu çalışmada kullanılan örnekte 0.92'lik bir medyan elde edilmiştir. Bu QC metriklerine dayanarak, veri kümesinden hangi çekirdeğin filtreleneceğine karar verilebilir. Analiz için, çekirdek başına 200'den az veya 10.000'den fazla gen, 'dan fazla mitokondriyal okuma ve 0.8'in altında bir karmaşıklık puanı olan çekirdekleri filtrelemeyi seçtik.
İlk kalite değerlendirmesi ve filtreleme adımını takiben, çekirdeklerin kümelenmesini görselleştirmek için bir UMAP oluşturulabilir. Kümeleme, SCT dönüşümü kullanılarak en değişken 2000 gene dayalı olarak gerçekleştirildi. İlk kümeleme adımları, örneğin yüksek mitokondriyal okumalara sahip çekirdekler gibi QC özelliklerinden herhangi birinin birlikte kümelenip kümelenmediğini kontrol etmek için kullanılabilir. Ayrıca, bu protokolde kullanılan DoubletFinder20 de dahil olmak üzere bazı ikili algılama yöntemleri için kümeleme bilgileri gereklidir. DoubletFinder, damlacık tabanlı platformun sağlayıcıları tarafından önerildiği gibi, beklenen çarpan oranı %4,8 olarak ayarlanmış olarak kullanıldı. Çiftin çıkarılmasından sonra, RNA'nın hücre lizizi üzerine sitoplazmadan salındığı ve emülsiyon içinde jel boncuklara (GEM'ler) dağıtıldığı ve aşağıdaki kütüphane hazırlama adımlarında amplifiye edildiği için özellikle tek çekirdek preparatlarında yaygın olan ortam RNA kontaminasyonu seviyesi tahmin edildi. Bu nedenle, ortamdaki RNA kontaminasyonunun doğasında var olan sorunu düzeltmek için çeşitli araçlar geliştirilmiştir (bakınız Tablo 3). Gen ekspresyon matrisini ayarlamak için ham arka plan matrisinin (yalnızca boş damlacıklar dahil) kullanıldığı ve gerçek gen ekspresyonu imzasını geliştirdiği R paketi decontX21'i kullandık.
Kümelenme ve az miktarda bulunan hücre tiplerini tespit etme yeteneği, çekirdek sayısına bağlıdır. Bu çalışma, QC filtreleme, çift çıkarma ve ortam RNA ayarlamasından sonra toplam 3817 çekirdekten IMAT'ta (Şekil 4B) beklenen tüm ana hücre tiplerini tespit etti. Bunlar arasında kök hücreler, fibro-adipojenik progenitörler (FAP'ler) ve olgun adipositlerin yanı sıra perisitler, düz kas hücreleri, bağışıklık hücreleri, kas progenitör hücreleri ve iskelet kası hücresi kontaminasyonundan kaynaklanan myonukleler vardı.
Genel olarak, bu protokolün, IMAT'ın biyolojisini ve hücresel kökenlerini çözmek için önemli olan hücre tipi açıklamalarının algılanmasını sağlayan yüksek çözünürlüklü tek çekirdek verileri ürettiğini gösterdik.

Şekil 4: Dizileme verilerinin kalite değerlendirmesi, kümeleme ve hücre tipi açıklaması. (A) Çekirdek başına tespit edilen gen sayısı, mitokondriyal okuma yüzdesi, ribozomal okuma yüzdesi ve log (10) gen sayısının tespit edilen log (10) okuma sayısına bölünmesiyle ölçülen hücre karmaşıklığı dahil olmak üzere, numunenin ve dizileme performansının değerlendirilmesi için temel metriklerin keman grafikleri. Her metrik için medyan değerler kapalı kutularda verilmiştir. Toplam çekirdek sayısı: 4662. (B) Tek tek çekirdeklerin kümelenmesini gösteren UMAP ve filtrelemeden sonra her küme için hücre tipi işaretleyici genlerin göreceli gen ekspresyonunu gösteren karşılık gelen DotPlot. Çekirdek sayısı: 3817. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.
Ek Dosya 1: Kalite kontrol ve kümeleme analizi için kod. Bu dosyayı indirmek için lütfen buraya tıklayın.