$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Ön qRT-PCR verileri, EWS'nin tekrarlayan ve düzensiz bölgesindeki alanin mutasyonlarına özgü tirozin içeren DAF adlı bir EWS / FLI mutantının, EWS / FLI hedef genlerini aktive etme yeteneğini koruduğunu, ancak kritik hedef genleri bastırmayı başaramadığını öne sürdü23. EWS etki alanındaki bu kalıntılar ile EWS/FLI işlevi arasındaki ilişkiyi daha iyi anlamak için yukarıda açıklanan ve Şekil 1'de özetlenen protokol kullanılmıştır. A673 Ewing sarkom hücreleri, FLI1'in3'UTR'sını hedefleyen bir shRNA ile viral olarak transdüklendi ve endojen EWS / FLI'nin tükenmesiyle sonuçlandı. Dört günlük seçimden sonra, EWS/FLI fonksiyonu, farklı 3XFLAG etiketli EWS/FLI mutant yapılarının viral transdüksiyonu ile kurtarıldı ve kurtarma için kontrol olarak boş vektör vardı. Δ22 adı verilen EWS etki alanından yoksun işlevsel olmayan bir mutant negatif kontrol olarak, wtEF adı verilen vahşi tip EWS/FLI ise pozitif kontrol olarak kullanılmıştır (Şekil 2A). DAF test yapısı olarak kullanılmıştır, ancak istenirse birden fazla test yapısı kullanılabilir. Hücreler, yapı ekspresyonlarının stabilize olmasını sağlamak için ek 10 gün boyunca seçildi ve daha sonra RNA (gDNA kaldırma adımı ile), protein ve koloni oluşturan tahliller için toplandı. Dört çoğaltma toplandı ve etkili devirme ve kurtarma gösteren temsili qRT-PCR ve batı lekeleri Şekil 2B-D'degösterilmiştir. DAF tarafından kurtarılan hücrelerin Şekil 2E'degösterildiği gibi koloniler oluşturamadığı belirtilmelidir , bu da onkojenik dönüşümün bozulacağını düşündürmektedir.
Çoğaltma doğrulama ve fenotipik testlerin tamamlanmasının ardından RNA, toplanan ~50 milyon 150 bp eşleştirilmiş uç okuma ile kütüphane hazırlama ve yeni nesil sıralama için Nationwide Çocuk Hastanesi Genomik Tıp Enstitüsü'ne gönderildi. Veriler fastq.gz dosyaları olarak döndürüldü. TrimGalore ile bu dosyalardan düşük kaliteli okumalar kesildi ve STAR, okumaları insan genomu hg19'a hizalamak ve gen başına okumaları saymak için kullanıldı. HG19, aşağı akış analizinde kullanılan EWS/FLI için diğer seçilmiş veri kümeleriyle uyumluluk amacıyla kullanılmıştır. Bu okuma sayıları, ilk 6 satırı Şekil 3'tegösterilen tüm örnekler için tek bir sayım matrisinde birleştirildi.
Sayımlar başlangıçta toplu normalleştirme olmadan DESeq2 üzerinden çalıştırıldı, ancak örnekten örneğe mesafenin görsel incelemesi, Şekil 4A'dakırmızı oklarla vurgulandığı gibi potansiyel kafa karıştırıcı toplu iş efektleri gösterdi. Bu muhtemelen kültürdeki hücrelerin geçişi ve her partinin işlenmesindeki farklılıklar nedeniyle ortaya çıktı. Toplu iş efektleri için normalleştirme ComBat ile gerçekleştirildi ve genellikle önerilir. Toplu olarak normalleştirilen verilerin örneklemden örneğe mesafeleri Şekil 4B'de gösterilmiştir. Toplu normalleştirmeden sonra, DESeq2 taban çizgisine göre üç yapı (wtEF, Δ22 ve DAF) için transkripsiyon profilleri oluşturmak için kullanıldı. "Ebeveyn" A673 hücreleri (burada "iLuc" olarak adlandırılan sahte nakavt ve sahte kurtarma) diferansiyel analize dahil edilirken, bu deney için referansın iEF hücreleri adı verilen EWS / FLI tükenmiş hücreler olduğunu unutmayın. Transkripsiyonal profil, iLuc örneğini iEF ile karşılaştırarak buradaki endojen protein için oluşturulabilir ve bu, kurtarma sisteminin nasıl çalıştığını anlamakla ilgilenebilir, ancak bu özel analizin amacı bu değildir. Mutantlar için oluşturulan transkripsiyon profilleri, iEF ile ilgili olarak pozitif (wtEF) ve negatif (Δ22) kontrolleri içerir, böylece bunlar diğer mutantlar için ölçüt olarak işlev görmelidir. Bu, bu örnekteki pozitif kontrol, başka bir yerde tartışıldığı gibi endojen EWS / FLI işlevini tamamen yeniden yakalamadığı için önemlidir7,23.
Şekil 5'teki ana bileşen analizi (PCA), DAF'ın transkripsiyon profilinin wtEF ve Δ22 arasında orta olduğunu ve kısmi işlevi doğruladığını göstermektedir. Ayrıca, örnekler arasında en değişken 1000 genin hiyerarşik kümelenmesi, DAF'ın EWS / FLI hedef genlerini bastıramadığını ve Şekil 6A ve Şekil S5'te gösterildiği gibi gen aktivasyon aktivitesini sadece kısmen koruduğunu göstermiştir. ToppGene analizi, DAF'ın aktive ettiği gen sınıflarının, DAF'ın işlevsel olmadığı EWS/FLI ile etkinleştirilen hedeflerden işlevsel olarak farklı olduğunu ileri sürmektedir (Şekil 6B). İlginçtir ki, wtEF tarafından kurtarılan aktif genlerin işlevi, ancak DAF tarafından değil, transkripsiyon kontrolü ve kromatin regülasyonu ile ilgili görünüyor. Koloni oluşumunun sonuçlarına dayanarak, bu çekirdek gen imzasından elde edilen genler EWS/FLI aracılı onkogenezdeki rolleri için daha fazla analiz edilmelidir. EWS/FLI aracılı gen baskısının önemi daha önce tanımlanmıştır17.
EWS/FLI'nin GGAA-mikrosatellit tekrar elemanları19,22için benzersiz bir bağlayıcı benzeşime sahip olduğu ve bu elemanlardaki bağlamanın aşağı akış gen regülasyon 11 ,15,18,20,22' yi yönlendirdiği bilinmektedir. Bu mikrosatellitler aktivasyon veya baskı ile ilişkili ve proksimal ila (< 5 kb) TSS veya distal (> 5 kb) TSS 25 olarak karakterizeedilmiştir. Ek olarak, TSS23'eproksimal olarak yüksek benzeşimli (HA) ETS motiflerine sahip EWS / FLI düzenlemeli genler vardır. DAF fonksiyonunun özelliklerini ve DAF'ın ne tür EWS/FLI ile aktive genleri kurtarabildiğini daha fazla analiz etmek için, bu farklı sınıflarla ilişkili genlerin diferansiyel ekspresyöz ifadesi analiz edildi. İlginçtir ki, DAF en çok GGAA-mikrosatellit aktive genlerini kurtarabildi, ancak Şekil 7'degörüldüğü gibi bir HA sitesinin yakınında aktif genleri kurtaramadı. Hiyerarşik kümelemede görüldüğü gibi, DAF motif sınıfları arasında EWS/FLI aracılı baskıyı kurtaramaz. Bu veriler, DAF'ın GGAA-mikrosatellitlerden TSS'ye hem proksimal hem de distal olarak bağlanmak ve etkinleştirmek için EWS'nin yeterli yapısal özelliklerini koruduğunu göstermektedir. Bu büyük olasılıkla GGAA'daki EWS / FLI etkinliği için önemli olduğu düşünülen bozulmamış SYGQ etki alanından kaynaklanmaktadır11. Bu veriler ayrıca, DAF'ta mutasyona uğrayan spesifik tirozinlerin HA sitelerinden EWS/ FLI aracılı gen düzenlemesinde ve gen baskısında önemli rol oynadığını ve daha fazla araştırmanın önemli bir alanını vurguladığını göstermektedir.

Şekil 1: İş Akışı. Transkriptomik olarak yapı işlevi eşlemesi gerçekleştirmek için adım adım yordamın tasviri. Hücreler ilk olarak yapı-işlev eşlemesi için gerekli yapı paketini ifade etmek için hazırlandı. İfadenin ardından, hücreler RNA ve protein için hasat edildi ve korelatif fenotipler için test edildi. Yapıların ekspresys doğrulandı ve bu süreç bağımsız biyolojik kopyaları toplamak için 3-4 kez tekrarlandı. RNA daha sonra yeni nesil sıralama (NGS) için gönderildi. Veriler alındığında, veriler kalite için kırpıldı, hizalandı ve transkript başına sayımlar hesaplandı. Deseq2 kullanılarak parti efektleri kontrol edildi ve transkriptomik imzalar ve diferansiyel ekspresyonu belirlendi. Hiyerarşik kümeleme ve aşağı akış analizi entegre diğer -omics veri kümeleri ve farklı yol veya fonksiyonel analiz dahil edilebilir. 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: Yapı ifadesinin ve korelatif tahlillerin doğrulanması. (A) Bu örnekte sınanan yapıları gösteren şematik. (B) Endojen EWS/FLI'nin devrilmesinin doğrulanması ve 3X-FLAG etiketli yapıların immünoblot tarafından ifade edilerek. (C,D) EWS/FLI (C) aktif hedef gen, NR0B1ve (D) bastırılmış hedef gen olan TGFBR2'deyapı aktivitesinin qRT-PCR ile doğrulanması. Veriler ortalama +/- standart sapma olarak sunulur. P değerleri Tukey'in dürüst önem testi ile hesaplandı. * p < 0.05, ** p < 0.01, *** p < 0.005 (E) Yapıların dönüşüm aktivitesini değerlendirmek için gerçekleştirilen yumuşak agar tahlillerinden koloni sayısı. P değerleri Tukey'in dürüst önem testi ile hesaplandı. * s < 0.05, ** p < 0.01, *** p < 0.005. Bu rakam Theisen,vd.23'tenuyarlanmıştırBu rakamın daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 3: Analiz için son harmanlanmış sayım verileri. Toplu olarak normalleştirilecek ve analiz edilecek tüm örnekler için gen sayısına sahip sayım dosyasının ilk 6 satırının ekran görüntüsü. 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: Örnek-numune mesafesi ısı eşlemleri. (A) Ham sayım verilerinin örnek kümelenini gösteren örnek-örnek mesafe çizimi. Hem toplu iş hem de örnekle kümelenen örnekler kırmızı oklarla gösterilir. (B) ComBat ile toplu normalleşmenin ardından örnek-örnek mesafe çizimi. Burada, tüm çoğaltma kümesinden örnekler toplu işlemden bağımsız olarak birlikte çoğalır. 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: Diferansiyel ekspresyon analizinin sonuçları. (A) Tüm numuneler için oluşturulan transkriptomik imzaların prensip bileşen analizi (PCA) grafiği güçlü numune içi kümeleme gösterir ve DAF'ın pozitif (wtEF) ve negatif (Δ22) kontroller arasında ara olduğunu gösterir. (B) Her yapıdaki genler için log2FoldChange'e çizilen -log(p-değerini) gösteren volkan çizimleri. Ayarlanmış p değeri 0,05 < ve |log2(FoldChange)| > 1 önemli olarak kabul edilir ve kırmızı renkle gösterilir. Panel 5B Theisen'den uyarlanmıştır, ve diğerleri23Bu rakamın daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 6: Gen sınıflarını tanımlamak için hiyerarşik kümeleme. (A) Tüm yapılarda en değişken 1000 genin hiyerarşik kümelenmesi ve temel olan iEF, DAF'ın EWS/FLI aracılı gen aktivasyonunu kısmen kurtar ettiğini gösterir. (B) Gen ontolojisi (moleküler fonksiyon) ToppGene'den elde edilen ve DAF tarafından kurtarılan veya kurtarılmayan EWS/FLI ile aktive genlerin fonksiyonel zenginleştirilmesini gösteren sonuçlar. Panel 6B Theisen'den uyarlanmıştır, ve diğerleri23Bu rakamın daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 7: Farklı transkripsiyon faktörü yanıt öğelerinin farklı yapılara ayrıntılı analizi: (A) Buradaki transkriptomik profillerle mevcut diğer veri kümelerini birleştirerek panel (B) ve (C) oluşturmak için kullanılan veri işlemeyi gösteren şematik. (B,C) Doğrudan EWS/FLI- (B) etkinleştirilmiş ve (C) bastırılmış hedeflerin farklı sınıflarının kurtarılmasını gösteren derleme. Genler sadece endojen EWS/FLI ile sapkınabilir diferansiyel ekspresyonu olan genlerdi. Her pasta grafikte gri, genlerin yapı tarafından kurtarılmayan kısmını tasvir eder. Kırmızı, genlerin farklı olarak aktive olan kısmını, mavi ise genlerin farklı olarak bastırılan kısmını tasvir eder. Bu rakam Theisen,vd.23'tenuyarlanmıştırBu rakamın daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Şekil S1: Fastq.gz dosyalarını HPC ortamına yükleme, kırpma ve hizalama. Bu rakamı indirmek için lütfen tıklayınız.
Şekil S2: Örnekler arasında okuma sayılarını harmanlama ve ComBat ile toplu normalleştirme çalıştırma. Bu rakamı indirmek için lütfen tıklayınız.
Şekil S3: DESeq2 çalıştırıyor ve diferansiyel ifade analizinin sonuçlarını ayıklıyor. Bu rakamı indirmek için lütfen tıklayınız.
Şekil S4: Çıktı analiz. Bu rakamı indirmek için lütfen tıklayınız.
Şekil S5: Gen sınıflarını tanımlamak için hiyerarşik kümeleme: Tüm yapılarda en değişken ilk 1000 genin hiyerarşik kümelenmesi ve taban çizgisi olan iEF, k kümeleri halinde sıralanır. Bu örnekte k=7, ancak bu parametre Şekil S4D'de gösterildiği gibi kullanıcı tarafından ayarlanır. Bu rakamı indirmek için lütfen tıklayınız.
Tablo S1: Küme ek açıklamalı genlerin (Ensembl gen kimliği) listesi. Bu tabloyu indirmek için lütfen tıklayınız.