Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
Click here for the English version

Genetics

RNA-Seq Analizi ile İnsan Primer Keratinositlerin In Vitro Farklılaşmasının Karakterizasyonu

Published: May 16, 2020 doi: 10.3791/60905
* These authors contributed equally

Summary

Burada sunulan kontak inhibisyonu ile insan primer keratinositlerin in vitro farklılaşması için adım adım bir prosedürdür ve rna-seq analizi ile moleküler düzeyde karakterizasyon uyguluyor.

Abstract

İnsan primer keratinositler genellikle epidermal farklılaşma ve ilgili hastalıklar üzerine çalışmalar için in vitro modeller olarak kullanılır. Çeşitli indüksiyon koşulları kullanılarak iki boyutlu (2D) batık şekillerde kültürlenmiş keratinositlerin in vitro farklılaşması için yöntemler bildirilmiştir. Burada tanımlanan kontak inhibisyonu ve RNA-seq tarafından sonraki moleküler karakterizasyon ile 2D in vitro keratinosit diferansiyasyon yöntemi için bir prosedürdür. Kısacası, keratinositler tanımlanmış keratinosit orta büyüme faktörleri ile tam olarak konca kadar takviye olarak yetiştirilir. Farklılaşma keratinositler arasındaki yakın temaslar tarafından indüklenen ve daha fazla orta büyüme faktörleri hariç uyarılır. RNA-seq analizleri kullanılarak, her iki 1) farklılaştırılmış keratinositlerin farklılaşma sırasında farklı moleküler imzalar gösterdiği ve 2) dinamik gen ekspresyonu deseni epidermal tabakalaşma sırasında hücrelere büyük ölçüde benzediği gösterilmiştir. Normal keratinosit farklılaşmasına karşılık gelince, transkripsiyon faktörü p63 mutasyonlarını taşıyan keratinositler, farklılaşma kusurları ile uyumlu morfoloji ve moleküler imzalar sergilerler. Sonuç olarak, bu protokol 2D in vitro keratinosit farklılaşması ve moleküler karakterizasyonu için adımları ayrıntıları, RNA-seq verilerin biyoinformatik analizi üzerinde durularak. RNA çıkarma ve RNA-seq yordamları iyi belgelenmiş olduğundan, bu protokolün odak noktası değildir. İn vitro keratinosit farklılaşması ve biyoinformatik analiz boru hattının deneysel prosedürü, sağlıklı ve hastalıklı keratinositlerde epidermal farklılaşma sırasında moleküler olayların incelenmesinde kullanılabilir.

Introduction

İnsan derisinden elde edilen insan birincil keratinositler genellikle epidermisin biyolojisi çalışmak için hücresel bir model olarak kullanılır1,2,3,4. Epidermisin tabakalaşması keratinosit farklılaşması ile modellenebilir, ya 2D batık monolayer moda veya 3D hava asansörü organotipik model2,3,5,6,7. Epidermal yapı ve fonksiyonu değerlendirmek için 3D modeller giderek daha önemli hale gelmiş olsa da, 2D farklılaşma modelleri, kolaylıkları ve analizler için çok sayıda hücre üretme olasılığı nedeniyle hala yaygın olarak kullanılmaktadır.

Serum eklenmesi, yüksek kalsiyum konsantrasyonu, düşük sıcaklık ve epidermal büyüme faktörüreseptörlerinininhibisyonu 2 ,3dahil olmak üzere, 2D keratinosit farklılaşmasıinin indüklemesi için çeşitli koşullar uygulanmıştır. Bu yöntemlerin her biri keratinosit farklılaşma belirteç genleri bir dizi tarafından doğrulanmış ve patolojik koşullar altında da dahil olmak üzere keratinosit farklılaşması değerlendirilmesinde etkili olduğu gösterilmiştir. Ancak, bu indüksiyon koşulları da marker genlerin belirli paneller incelendiğinde onların farklılaşma verimliliği ve kinetik farklılıkları göstermektedir2,3.

Bu yöntemlerden biri keratinosit kontak inhibisyonu ve kültür ortamı nda büyüme faktörlerinin tükenmesi içerir8. Bu keratinositler hücreler tam yoğunluğa ulaştığında kendiliğinden ayırt edebilirsiniz gösterilmiştir.  Kültür ortamındaki büyüme faktörlerinin dışında kılması farklılaşmayı daha da artırabilir. Temas inhibisyonu ve büyüme faktörlerinin tükenmesini birleştiren yöntem, çeşitli epidermal belirteçler kullanırken normal tabakalı epidermis benzer gen ekspresyonu desenleri ile farklılaştırılmış keratinositler oluşturmak için gösterilmiştir3, Bu model normal keratinosit farklılaşması için uygun olduğunu düşündürmektedir. Son zamanlarda, bu modeli kullanarak keratinosit farklılaşması iki kapsamlı gen ekspresyonu analizleri bildirilmiştir9,10. Araştırmacılar bu modeli moleküler düzeyde doğruladı lar ve normal ve hastalıklı keratinosit farklılaşmasını incelemek için kullanılabileceğini gösterdiler.

Bu protokol, in vitro diferansiyasyon yöntemi ve RNA-seq kullanarak farklılaşmış hücrelerin moleküler analizi için prosedürü açıklar. Ayrıca, farklılaşma günü 0 (proliferasyon aşaması), 2. gün, 4. Diferansiye keratinositlerin epidermal tabakalaşma sırasında hücrelere büyük ölçüde benzeyen gen ekspresyonu desenlerini gösterdiği gösterilmiştir. Bu yöntemin deri patolojisi üzerinde çalışılıp kullanılamayacağını incelemek için, ektoderaktili, ektodermal displazi ve yarık dudak/damak (EEC) sendromu11,12olan hastalardan elde edilen transkripsiyon faktörü p63 mutasyonlarını taşıyan keratinositleri araştırmak için aynı deneysel ve analiz boru hattını uyguladık. Bu protokol, keratinositlerin in vitro farklılaşmasına ve rna-seq'nin sonraki biyoinformatik analizine odaklanır. RNA çıkarma, RNA-seq numune hazırlama ve kütüphane inşaatı gibi komple yordamdaki diğer adımlar iyi belgelenmiştir ve özellikle yaygın olarak kullanılan birçok ticari kit kullanılırken kolayca takip edilebilir. Bu nedenle, bu adımlar yalnızca kısa bir süre protokolü açıklanır. Veriler, bu boru hattının sağlıklı ve hastalıklı keratinositlerde epidermal farklılaşma sırasında moleküler olayların incelenmesi için uygun olduğunu göstermektedir.

Subscription Required. Please recommend JoVE to your librarian.

Protocol

Deri biyopsileri sağlıklı gönüllülerin veya p63 mutasyonu olan hastaların gövdesinden alınarak primer keratinosit kültürünü nisabı kültürü ne şr; Radboud Üniversitesi Nijmegen Tıp Merkezi ("Komiser Mensgebonden Onderzoek Arnhem-Nijmegen") etik komitesi tarafından insan birincil keratinositlerin kurulmasına ilişkin tüm prosedürler onaylanmıştır. Bilgilendirilmiş onay alındı.

1. Temas inhibisyonu ile insan primer keratinosit farklılaşması

  1. Gerektiğinde Keratinosit Bazal Orta(Malzeme Tablosu)keratinosit büyüme ortamı (KGM) hazırlayın.
  2. Önceden hazırlanmış stoklar kullanılarak 500 mL proliferasyon ortamı olun (KGM-pro; bkz. Tablo 1).
  3. 500 mL farklılaştırma ortamı olun (KGM-dif; bkz. Tablo 2).
    NOT: KgM takviyeleri ile 4 °C'de iki hafta tutulabilir. Daha uzun depolama için -20°C'de saklanabilir.
  4. Hücre proliferatif kapasitesine bağlı olarak 5.0-20 x 103 hücre/cm2yoğunluğunda tohum primer keratinositler. Hücreleri kapsayacak kadar KGM-pro ortam ekleyin.
    NOT: 1) KGM ortamı içinde düzenli hücre kültürünün ayrıntıları Lonza13web sitesinde bulunabilir. 2) Tohumlama yoğunluğu her hücre hattı için test edilmelidir. Örneğin, normal bir primer keratinosit hattı 5.0 x 103 hücre/cm2yoğunlukta tohumlanabilirken, p63 transkripsiyon faktöründe mutasyontaşıyan ve daha az proliferatif olan bir çizgi 20 x 103 hücre/cm2yoğunlukta tohumlanmalıdır. 3) Deneysel kurulumbağlı olarak, birkaç tabak veya hücre kuyuları çoğaltmalar farklı farklılaşma gün örnekleri toplama sağlamak için tohumlu olmalıdır.
  5. KGM-pro ortamı ile hücreleri tohumlamadan sonra iki gün yenileyin (0. günde tohumlama, 3. gün 1saat için ferahlatıcı). Daha sonra, KGM-pro ortamı ile her gün yenileyin. Hücreleri düzenli olarak kontrol edin, en azından her gün.
  6. Hücreler% 90'dan fazla dfluent olduğunda KGM-dif orta değiştirerek hücre farklılaşmasını neden olur. Orta nın KGM-dif olarak değiştirilmesi günü, farklılaşma günü 0 olarak tanımlanır. Daha fazla RNA analizleri için hücreleri ilk olarak DPBS ile 2x yıkayarak, daha sonra lysis arabelleği ekleyerek (RNA ekstraksiyon kitinden) toplayın.
    NOT: Yukarıda bahsedilen tohumlama yoğunluğu ile hücrelerin 7-10 gün içinde biraraya gelmesi gerekir. Hücreler daha erken bir araya ulaşmak için daha yüksek bir başlangıç yoğunluğu ile tohumlanmış olabilir. Hücreler çoğalan değilse, daha uzun süre beklemek tam olarak oluşan hücrelere yol vermeyebilir. Daha yüksek bir tohumlama yoğunluğu gerekebilir.
  7. KGM-dif ortamı ile hücreleri her gün yenileyin ve 2, 4. ve daha fazla RNA çıkarma için 7.

2. RNA çıkarma

  1. Toplam RNA'yı yalıtın. Bu ticari bir RNA izolasyon kiti (Zymo Quick-RNA MicroPrep RNA gibi) veya fenol14kullanılarak yapılabilir. Ancak, DNA kontaminasyonunu ve fenolönlemek için RNA-seq örnekleri için DNA tedavisi içeren bir izolasyon kiti şiddetle tavsiye edilir. Ticari bir RNA izolasyon kiti örneği Malzeme Tablosu'ndaverilmiştir.
  2. RNA konsantrasyonu spektrometre ile ölçün. Hem DNA hem de RNA 260 nm'de bir emilim tepe noktasına sahiptir.
    NOT: 1) RNA'nın saflığını değerlendirmek için 260 nm ve 280 nm'deki absorbanslar arasındaki oran kullanılır. 2.0 civarında bir oran genellikle 'saf' RNA olarak kabul edilir. Oran 2,0'dan önemli ölçüde düşükse, numune proteinler, fenoller veya diğer maddeler gibi 280 nm'de emen kirletici maddelerle kontamine olabilir. 2) 260 nm ve 230 nm'deki absorbanslar arasındaki oran nükleik asit saflığını ölçer. 2.0 ile 2.2 arasında bir oran bekleniyor. Oran önemli ölçüde düşükse, numune EDTA, etanol veya diğer maddeler gibi 230 nm'de emilen kirletici maddelerle kontamine olabilir.

3. RNA kalite kontrolü

  1. RNA Bütünlük Numarası (RIN) nicel doğrulamak için bir biyoanalizör RNA Pico çip veya nitel bir ölçü için bir jel üzerinde RNA çalıştırın. Genel olarak, en az sekiz rin son derece tavsiye edilir. Ancak, dokulardan RNA çıkarılırken RIN daha düşük olabilir.
    NOT: İsteğe bağlı olarak, rna malzeme üzerinde qPCR tarafından kalite kontrol yapılabilir.

4. RNA-seq kütüphane hazırlama

NOT: RNA-seq kitaplık hazırlama genellikle ticari bir kit veya ticari ayarları altında yapılır. Açıklanan protokol ticari bir kit uyarlanmıştır, KAPA RNA HyperPrep Kit RiboErase ile (Illumina), gerekli tüm adımları kısa bir açıklama ile: insan ribozomal RNA,RNA parçalanma, birinci iplikçik sentezi, ikinci iplikçik sentezi ve A-kuyruk, ve her adımdan sonra temizlik oligo hibridizasyon ile rRNA tükenmesi15. Diğer kütüphane hazırlama kitleri de bu amaç için kullanılabilir. Oluşturulan cDNA kitaplığın kalitesi genellikle daha tutarlı olduğundan, bu adımı ticari olarak kullanılabilen bir kit kullanarak gerçekleştirmeniz önerilir. 2) 1x kütüphane hazırlığı için aşağıdaki adımlar açıklanmıştır. Birkaç numune hazırlarken, %10 ekstra hacimli ana karışımlar yapın.

  1. Oligo hibridizasyon ve rRNA tükenmesi
    1. Oligo hibridizasyon ana karışımını (toplam = 11 μL, 4,4 μL hibridizasyon tamponu, 4,4 μL hibridizasyon oligosu ve 2,2 μL RNase içermeyen su) ve tükenme ana karışımını (toplam = 5,5 μL, 3,3 μL tükenme tamponu ve 2,2 μL'lik Tükenme seh H) hazırlayın.
    2. PCR reaksiyon programını bir termocycler üzerine ayarlayın: 2 dk için 95 °C; rampa -0.1 °C/s'de 45 °C'ye kadar; 45 °C duraklama; 30 dk için 45 °C; 4 °C sonsuza kadar.
      NOT: Amplifikasyon için gerekenden iki kat daha fazla RNA miktarı ile başlayan düşünün. İlk denemede, aşağıdaki adımları ile devam etmek için miktarın yarısını kullanın. Bu, amplifikasyon döngülerinin yetersiz veya aşırı amplifikatöre dönüşmesi durumunda, bu adımları farklı sayıda döngüyle yinelemek için hala malzeme olduğundan emin olmaktır (bkz.
    3. 10 μL RNAse içermeyen suda 25 ng ile 1 μg arasında toplam RNA kullanın, 10 μL oligo hibridizasyon ana karışımı ekleyin. Numuneleri önceden programlanmış termocycler'a yerleştirin ve programı başlatın.
    4. Program 45 °C'de duraklama adımına ulaştığında, termocycler'dan çıkarmadan 20 μL hibridizasyon reaksiyonuna 5 l'lik tükenme ana karışımı ekleyin. Birkaç kez yukarı ve aşağı boru ile iyice karıştırın.
    5. Termocycler programına devam edin ve tükenme adımı (30 dk için 45 °C) devam edin.
      NOT: 1) Protokolün bir sonraki bölümü için tükenme adımı sırasında rRNA tükenmesi (örneğin, KAPA saf boncuklar) için boncukları oda sıcaklığına (RT) koyun. 2) Reaktifler doğrudan rRNA tükenmesi temizleme sonra gerekli olduğundan, DNase sindirim ana karışımı (bölüm 4.2 için) hazırlamayı düşünün.
    6. 2.2x boncuk bazlı KAPA Pure Boncuk temizliği yapın: RNA karışımı (25 μL) ve KAPA Pure Boncukları (55°L) birleştirerek, RNA karışımındaki boncukları birkaç kez yukarı ve aşağı boruyla ıslatarak iyice askıya alın.
    7. RNA'nın boncuklara bağlanmasını sağlamak için 5 dakika boyunca RT'de kuluçkaya yatırın ve sıvı temizlenene kadar boncukları yakalamak için tüpü mıknatıslı rafa yerleştirin ve 60 μL'lik süpernatantı dikkatlice çıkarın ve atın.
    8. Tüpü mıknatısın üzerinde tutarak, boncukları etanolle 30 ≥ süre kuluçkaya yatırarak 200 μL %80 etanol ile yıkayın ve etanolatın. İkinci yıkamadan sonra boncukrahatsız etmeden tüm artık etanol kaldırmaya çalışın.
    9. 3-5 dakika veya tüm etanol buharlaşana kadar RT'deki boncukları kurutun.
      NOT: Boncukların aşırı kurutulması verimin azalmasına neden olabilir.
  2. DNase sindirim
    1. DNase sindirim ana karışımını hazırlayın (toplam = 22 μL, 2,2 μL DNase tamponu, 2 μL DNase ve 17,8 μL RNase içermeyen su ile) .
    2. Birkaç kez yukarı ve aşağı borulama tarafından DNaAse sindirim ana karışımı (20 μL) boncuk resuspend. RNA'yı boncuklardan çıkarmak için 3 dakika boyunca RT'deki tüpü (ler) kuluçkaya yatırın.
    3. Sıvı temiz olana kadar boncukları yakalamak için tüpü mıknatısrafına yerleştirin ve 20 μL'lik süpernatantı temiz bir tüpe dikkatlice aktarın.
    4. 37 °C'de 30 dakika boyunca tüp(ler) ile tüp(ler) inkübül.
      NOT: Reaktifler Doğrudan DNase sindirim temizliği nden sonra gerekli olduğundan, RNA elüsasyon, parçalanma ve astar ana karışımları (bölüm 4.3 için) hazırlamayı düşünün.
    5. 4.1.6-4.1.9'u takip ederek 2,2 x boncuk bazlı temizlik yapın.
  3. RNA elüsasyonu, parçalanma ve astarlama
    1. Fragment, Prime ve Elute Tampon (1x) ana karışımını 11 μL Fragment, Prime ve Elute Tampon (2x) ve 11 μL RNase içermeyen su ile hazırlayın.
    2. Arıları saflaştırılmış, DNase ile işlenmiş RNA ile 22 μL'lik Parçalanma, Prime ve Elute Tampon (1x) ile birkaç kez yukarı ve aşağı borular oluşturarak iyice yeniden askıya alın.
    3. Oda sıcaklığında 3 dakika boyunca cin s kapalı RNAelute, sonra sıvı temizlenir kadar boncuk yakalamak için bir mıknatıs üzerine tüp (ler) yerleştirin.
    4. 20 μL'lik supernatant'ı dikkatlice bir tüpe aktarın. Tüp(ler) boncukile atın.
      NOT: Numuneler -20 °C'de 24 saat ≤ için depolanabildiği için güvenli bir durma noktasıdır.
    5. Tüpü bir termocycler yerleştirin ve 94 °C'de 6 dk boyunca parçalanma ve astarlama gerçekleştirin. Yaklaşık 200-300 bp uzunluğunda parçalar elde edilir.
      NOT: 100-200 bp parçaları için 94 °C'de 8 dk kuluçka. Kısmen bozulmuş RNA kullanırken, bozulmanın şiddetine bağlı olarak 85 °C'de 1-6 dakika arasında parçabölün.
    6. Tüp(ler) buz üzerine yerleştirin ve hemen ilk iplikçik sentezine devam edin.
  4. İlk iplikçik sentezi
    1. İlk iplik sentezi ana karışımını hazırlayın (toplam = 12°L, 11 μL ilk ipliksentezi tamponu ve 1 μL KAPA komut dosyası).
    2. PCR reaksiyon programını bir termocycler üzerine ayarlayın: 10 dk için 25 °C; 15 dk için 42 °C; 15 dk için 70 °C; 4 °C sonsuza kadar.
    3. Buz üzerinde, 20 μL parçalı astarlanmış RNA ile 10 μL ilk iplik sentezi ana karışımını birleştirin. Tüp(ler) buz üzerinde tutun, yavaşça yukarı ve aşağı birkaç kez reaksiyon borulama tarafından iyice karıştırın.
    4. Tüpü önceden programlanmış termocycler'a yerleştirin ve programı başlatın.
  5. İkinci iplikçik sentezi ve A-kuyruklu
    1. Buz üzerinde ikinci iplikçik sentezi ve A-kuyruklu ana karışımı (toplam = 33°L, 31 μL ikinci iplik sentezi tamponu ve 2 μL ikinci iplikçik sentezi ve A-kuyruklu enzim karışımı) hazırlayın.
    2. BIR termocycler pcr reaksiyon programı ayarlayın: 16 °C 30 dakika için; 10 dk için 62 °C; 4 °C sonsuza kadar.
    3. 30 μL'lik ilk iplik sentezi ürününün 30 μL'lik ikinci ipliksentezi ve A-kuyruklu ana karışımını birleştirin ve buz üzerinde birkaç kez yukarı ve aşağı borular oluşturarak iyice karıştırın.
    4. Tüpü önceden programlanmış termocycler'a yerleştirin ve programı başlatın.
  6. Adaptör ligasyonu ve ligasyon sonrası temizleme
    1. Seyreltik stok adaptörleri (NEXTflex DNA barkodları, 25 μM) RNase içermeyen suda 7 μM'ye 3,57 x.
    2. Adaptör ligasyon ana karışımını hazırlayın (toplam = 50°L, 40 μL ligasyon tamponu ve 10 μL DNA ligaz).
    3. 60 μL ikinci iplik sentezi ürününün, 45 μL adaptör ligasyon ana karışımını ve 5 μL seyreltilmiş adaptörleri birleştirin. Buz üzerinde birkaç kez yukarı ve aşağı boru lar oluşturarak iyice karıştırın.
    4. Tüpleri 20 °C'de 15 dakika kuluçkaya yatırın ve ilk ligasyon sonrası temizleme ile hemen devam edin.
    5. Adaptör ligatlı DNA (110 μL) ve KAPA saf boncukları (70°L) birleştirerek 0,63x boncuk bazlı bir temizlik gerçekleştirin, ardından DNA karışımındaki boncukları birkaç kez yukarı ve aşağı borular oluşturarak iyice askıya alın.
    6. Adımları 4.1.7-4.1.9'u yineleyin.
    7. Mıknatıstan tüpleri çıkarın. 10 mM Tris HCL'nin (pH = 8.0-8.5) 50 μL'lik boncukları iyice askıya alın ve boncukları 2 dakika boyunca RT'de kuluçkaya yatırın.
    8. Plakayı mıknatısın üzerine yerleştirin ve sıvı temize çıkana kadar bekleyin. Açık süpernatantın 50 μL'sini yeni bir tüpe aktarın.
      NOT: 4 °C'de 24 saatten az bir süre için güvenli durma noktası.
    9. 50 μL boncuk ile saflaştırılmış adaptör ligated DNA'yı 35 μL PEG/NaCL çözeltisi ile birleştirerek 0,7 x boncuk bazlı temizlik yapın. Girdap ile iyice karıştırın.
    10. Temizleme adımlarını 4.1.7-4.1.9 gerçekleştirin. 10 mM Tris HCL'nin (pH = 8.0-8.5) 20 μL'lik kısmını iyice askıya alın ve boncukları 2 dakika boyunca RT'de kuluçkaya yatırın.
    11. Plakayı mıknatısın üzerine yerleştirin; sıvı temizlene ne kadar bekleyin. Açık supernatantın 20 μL'sini yeni bir tüpe aktarın ve bölüm 4.7'ye ilerleyin.
      NOT: 4 °C'de 1 haftadan az veya -20 °C'de 1 aydan az bir süre için güvenli durma noktası.
  7. Kütüphane amplifikasyon ve temizleme
    1. Kütüphane amplifikasyon ana karışımını hazırlayın (toplam = 33 μL, 27,5 μL 2x KAPA HiFi HotStart ReadyMix ve 5,5 μL 10x kütüphane amplifikasyon astar karışımı).
    2. Aşağıdaki parametreleri kullanarak BIR termocycler PCR reaksiyon programı ayarlayın. 45 s için 98 °C; N döngüleri (aşağıdaki nota bakınız): 98 °C 15 s; 30 s için 60 °C; 30 s için 72 °C; 72 °C 1 dk için; ve 4 °C sonsuza kadar.
      NOT: Kitaplık amplifikasyonu için döngüler giriş malzemesine bağlıdır. Kabaca şöyle tahmin edilebilir: RNA = 25-100 ng, N = 11-15 döngüleri ile başlayarak; 100-250 ng, N = 9-12 döngüleri; 250-500 ng, N = 7-10 döngüleri.
    3. Güçlendirilmiş kütüphane DNA'sı (50°L) ve KAPA saf boncukları (40°L) birleştirerek 0,8x boncuk bazlı bir temizlik gerçekleştirin, ardından DNA karışımındaki boncukları birkaç kez yukarı ve aşağı borular oluşturarak iyice askıya alın.
    4. Temizleme adımlarını 4.1.7-4.1.9 gerçekleştirin. 10 mM Tris HCL'nin (pH = 8.0-8.5) 50 μL'lik boncukları iyice askıya alın ve boncukları 2 dakika boyunca RT'de kuluçkaya yatırın.
    5. Açık supernatant'ın 50 μL'sini yeni bir tüpe aktarın ve ikinci kütüphane amplifikasyon temizliğine devam edin.
    6. Güçlendirilmiş kütüphane DNA'sı (50°L) ve KAPA saf boncukları (50°L) birleştirerek 1x boncuk bazlı bir temizlik gerçekleştirin, ardından DNA karışımındaki boncukları birkaç kez yukarı ve aşağı borular oluşturarak iyice askıya alın.
    7. Temizleme adımlarını 4.1.7-4.1.9 gerçekleştirin. 10 mM Tris HCL'nin (pH = 8.0-8.5) 22 μL'lik boncukları iyice askıya alın ve boncukları 2 dakika boyunca RT'de kuluçkaya yatırın.
    8. Açık süpernatantın 20 μL'sini yeni bir tüpe aktarın Qbit ve biyoanalizör ölçümlerine ilerleyin.
  8. Konsantrasyon ve parçalanma boyutu
    1. Örneklerin DNA konsantrasyonlarını ölçün. Son derece hassas floresan boya bazlı kiti kullanarak (örneğin, Denovix Qbit, bkz. Malzeme Tablosu),veya konsantrasyon 0,5 ng/μL'nin altında görünüyorsa, daha hassas bir qPCR yaklaşımı kullanarak yeniden ölçün (örneğin, Kappa Quantification, Bkz. Malzeme Tablosu).
    2. Yüksek duyarlı elektroforez (örneğin, bir biyoanalizör) kullanarak kütüphanenin parça boyutunu belirleyin.
      NOT: İsteğe bağlı olarak, qPCR ile kalite kontrolü cDNA yerine seyreltilmiş hazırlanmış bir kitaplık kullanılarak sıralama (Ek) önce yapılabilir.
  9. Sıralama
    1. Sıralama için kütüphaneyi gönderin. RNA-seq diferansiyel gen ekspresyonu analizi için, numune başına 10-25 milyon arasında bir okuma nın sıralama derinliği genellikle yeterlidir.

5. Veri ön işleme

  1. Kalite kontrol Fastq dosyaları
    1. Düşük kaliteli baz çiftleri ve Fastq dosyaları içinde boş okuma kaldırmak için Trimgalore16'yı indirin ve kurun.
      NOT: 1) Burada bahsedilen çoğu yazılım sadece Linux'ta çalışır. Bir Windows bilgisayarıyla sınırlıysa, MobaXterm veya Putty yazılımını kullanarak bir Linux sunucusunda analizler çalıştırmak mümkündür. Genom indeksleme için çok sayıda rastgele erişimbelleği (RAM) gerektiğini unutmayın (yaklaşık 64 GB). 2) Bir conda ortamına gerekli tüm yazılım yükleme şiddetle tavsiye edilir. Bu kolay yazılım yükleme ve paket yönetimi sağlayacaktır. Conda hakkında daha fazla bilgi için adresini ziyaret edin.
    2. 'RNA_seq_KC_diff' klasörü gibi veriler için bir dizin (klasör) oluşturun. "cd /home/user/RNA_seq_KC_diff/" yazarak bunu çalışma dizini olarak ayarlayın.
      NOT: Klasör yapısına genel bir bakış için, ek kodlama dosyalarındaki 'folder_structure.txt' bölümüne bakın.
    3. Çalışma dizininde belirli adlarla birkaç klasör oluşturun: 'Fastq', 'CRCh38_fasta', 'CRCh38', 'komut dosyaları', ve 'Eşleme'.
    4. Sıralı verilerin fasta dosyalarını fastq klasörüne taşıyın. Alternatif olarak, disk alanı kazanmak için Fastq dosyalarına "ln home/user/old_location_fastq/FastqFile home/user/RNAseq_KC_diff/Fastq/Fastqfile" komutunu kullanarak yumuşak bağlantılar oluşturun.
    5. Fastq dosyalarında Trim galore çalıştırın, bash içine yapıştırkopyalamak için kod için TrimGalore.txt dosyasına bakın. Gerektiğinde komut içindeki klasör adlarını ve ayarlarını değiştirin.
  2. Eşleme
    1. Okumaların haritasını çıkarmak için bir genom indirin, örneğin, hg38 ensembl sürüm 97 (maskesiz sürüm): ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/Homo_sapiens.GRCH38.primary_assembly.fa.gz ; ve ilgili gen ek açıklama dosyası: hg 38 gen ek açıklama ensembl sürüm 97, ftp://ftp.ensembl.org/pub/release-97/gtf/homo_sapiens/Homo_sapiens.GRCh38.97.gtf.gz. Her iki dosyayı da CRCh38_fasta klasörüne aktarın.
      NOT: Alternatif olarak genomun farklı bir sürümünü kullanın, örneğin, hg38'in UCSC sürümü, gen iLikleri yerine transkript d'lerine haritalama tercih ediliyorsa.
    2. STAR 2.7.117yükleyin.
    3. Oluşturmak genom komut star 2.7.1 kullanarak bir in-house referans genom oluşturun (Ek Kodlama Dosyalarıbakınız). İş parçacığı miktarını işleyicinin sahip olduğu işlem miktarıyla değiştirin (--runThreadN X). Gerektiğinde bu komut dosyasındaki klasör adlarını ve ayarlarını değiştirin. Kopyalayarak/yapıştırarak bash'a çalıştırın.
    4. Bam dosyaları dizine ekin ve wiggle dosyaları sıkıştırmak için gzip için samtools 1.918 yükleyin.
    5. Trim Galore doğrulanmış sıralama okumalarını insan genomu montajı hg38 (ensembl release 97) star 2.7.1 kullanarak ve ardından samtools ve gzip kullanarak indeksleme ile hizalayın. Bu komut dosyası map_fastq.txt kullanılarak yapılmalıdır (Bkz. Ek Kodlama Dosyaları). Gerektiğinde bu komut dosyasındaki klasör adlarını ve ayarlarını değiştirin. Kopyalayarak/yapıştırarak bash'a çalıştırın.
      NOT: STAR eşleme, bir bam dosyası, bir kıpırdanma dosyası ve deseq2 için giriş olarak kullanmak üzere R tarafından concatenated doğrudan kullanılabilir *_ReadsPerGene.out.tab adlı sayımları tablo okumak dahil olmak üzere çeşitli çıkış dosyaları oluşturur.
  3. Genom tarayıcı görselleştirme
    1. Big2bw script19ile bigwig dosyaları oluşturmak için UCSC genom tarayıcı araçları wigToBigWig yükleyin.
    2. 'convertBigWigChroms.py' ve 'CRCH38_ensembl2UCSC.txt' dosyalarını tamamlayıcı kodlama dosyalarından çalışma dizinindeki bir klasöre 'komut dosyaları' olarak aktarın.
    3. ConvertBigWigChroms.py çalıştırılabilir hale getirin (linux makinesinde 'chmod +x komutları/convertBigWigChroms.py' komutunu kullanarak).
    4. Komut dosyası wig2bwkullanarak, Wiggle dosyaları bigWig dosyaları oluşturun.txt (Ek Kodlama Dosyalarıbakınız). Kopyalayarak/yapıştırarak bash'a çalıştırın.
    5. Görselleştirme için UCSC genom tarayıcısı veya Tümleştirici Genomik Görüntüleyici (IGV) üzerine BigWig dosyalarını girin.

6. RNA-seq veri analizi

  1. Örnek verileri Deseq2'ye girme
    1. Rstudio (sürüm 1.1.456) ve R (sürüm 3.6) indirin ve yükleyin.
    2. Gerekli tüm R paketlerini yükleyin (bkz. Ek kodlama dosyaları).
      NOT: Tüm programlama kodunu ve çıktısını gösteren ayrıntılı bir R-markdown dosyası için, ekteki dosyalara bakın: "RNA_seq_kc_differentiation_wt.html" ve "RNA_seq_kc_differentiation_patient.html". Tüm adımların genel bir açıklaması aşağıda verilmiştir.
    3. ReadsPerGene.out.tab dosyalarından bir sayım tablosu oluşturun.
    4. Tüm dosya adlarını, farklılaşma gününü ve diğer ilgili örnek verileri içeren bir örnek veri dosyası yazın. Örnek bir örnek dosya için, Ek Kodlama Dosyaları'ndaki"sample_data_example.csv" adlı dosyaya bakın.
    5. Hem sayılabilir hem de örnek verileri içeren bir Deseq220 nesnesi oluşturmak için sayım tablosunu ve örnek verileri kullanın.
  2. Gen ekspresyonu normalizasyonu, örnek mesafesi ve PCA
    1. Deseq2 rld veya vst normalleştirme kullanarak Deseq2 nesnesindeki sayım tablosunu normalleştirin. Rld normalleştirme tercih edilir, ancak birçok örnekile vst normalizasyonu çok daha hızlıdır.
    2. R'deki "dist" fonksiyonu kullanılarak normalleştirilmiş okuma sayısı yoğunluklarına dayalı çizim örneklem mesafesi ve daha sonra örnek uzaklığı temel alınarak "hclust" kümeleme gerçekleştirir. Pheatmap paketini kullanarak ısı haritasının kendisini çizin.
    3. Deseq2'nin plotPCA işlevini kullanarak normalleştirilmiş okuma sayısı yoğunluklarının bir PCA çizimini oluşturun.
      NOT: PCA, araştırmacı veri analizinde bir araç olarak hizmet vermektedir ve farklı örnekler arasındaki mesafeyi ve akrabalığı görselleştirmek için kullanılabilir.
  3. Diferansiyel/yüksek değişkenli gen ekspresyonu analizi
    1. Deseq2 sonuç işlevini kullanarak diferansiyel genleri hesaplayın. Ancak, çift yönlü karşılaştırmaların geçerli olmadığı birkaç zaman puanındaki değişiklikleri değerlendirirken, son derece değişken genler kullanın. Bu durumda, rowVars fonksiyonu ile farklı zaman noktalarından örnekler arasındaki varyans sipariş yoluyla ilk 500 son derece değişken genleri ayıklayın.
      NOT: Analizlerimizde, kontrol örnekleri için, daha fazla analiz için en yüksek değişken 500 gen (farklılaşma günlerinde normalleştirilmiş yoğunluklarının en yüksek standart sapması olan genler) kullanılmıştır. Ancak, hastalıklı vs kontrol analizi için hastalık ve zaman içinde sağlıklı kontroller arasında farklı ifade genler deseq2 tarafından 1 x 10-4 bir cutoff olarak düzeltilmiş p-değeri kullanılarak hesaplanmıştır. Başka bir olası filtreleme adımı, diferansiyel genlerin belirli bir kıvrım değişiminden daha az değişen kesimi dışlamaktır.
    2. Diferansiyel veya çok değişken genler üzerinde farklı ifade desenlerine göre kümelemek için kmeans kümeleme gerçekleştirin.
    3. Pheatmap paketini kullanarak bir ısı haritasında diferansiyel veya çok değişken genleri görselleştirin. Isı haritasında çizilen yoğunluk, ortanca değeri çıkarılmış deseq2 normalleştirilmiş yoğunluk.
  4. Gen Ontolojisi (GO) ek açıklama zenginleştirme analizi
    1. Tek bir örnekte 10'dan fazla sayıma sahip tüm genleri alarak ifade edilmiş arka plan genlerinin bir listesini oluşturun.
    2. GOrilla21gibi çevrimiçi araçları kullanarak GO analizi gerçekleştirin. Kümelerdeki diferansiyel/yüksek değişkengenlerin listelerini "gen listesi" olarak ve arka plan genlerini karşılaştırma için arka plan olarak kullanın.
      NOT: Daha gelişmiş R kullanıcıları Go vadeli zenginleştirme analizini otomatikleştirmek için clusterProfiler22 paketini kullanabilir.

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

Normal keratinosit farklılaşması ve RNA-seq analizi
Bu deneyde, beş kişiden elde edilen keratinosit çizgileri farklılaşma ve RNA-seq analizleri için kullanılmıştır. Şekil 1, farklılaşma ve RNA-seq analizi sonuçlarının deneysel prosedürünü özetler. Normal keratinositlerin in vitro farklılaşma prosedürlerine genel bir bakış ve farklılaşma sırasında hücre morfolojisi değişiklikleri Şekil 1A'dagösterilmiştir. İlke komponent analizi (PCA), farklılaşma geçiren keratinositlerin bağlı ancak farklı genel gen ekspresyonu profilleriolduğunu göstermiştir(Şekil 1B). Çok değişken genler farklılaşma sırasında gen ekspresyonu dinamiklerini ve desenlerini görselleştirmek için kmeans tarafından kümelenmiştir(Şekil 1C).

Genlerin her kümesi keratinosit farklılaştırma özelliği genleri (örneğin, proliferasyon için KRT5, erken ve orta farklılaşma için KRT1ve KRT10 ve geç farklılaşma için IVL/LOR/FLG) tarafından temsil edildi. Yüksek değişkengen kümelerinin gen ontolojisi (GO) ek gösterimi analizi(Şekil 1C)bu gen kümelerinin gen fonksiyonlarında (örneğin, farklılaşmanın orta evrelerinde keratinizasyon; ve epidermal hücre farklılaşması, keratinosit farklılaşması ve geç evrelerde peptit çapraz bağlanması) açık bir fark göstermiştir; Şekil 1D). Birkaç farklılaşma belirteçlerinin protein ekspresyonu batı lekeleme ile ölçüldü (Şekil 1E).

P63 mutant keratinosit farklılaşması ve RNA-seq analizi:
İkinci deneyde, hücre morfolojisi ve gen ekspresyonu farklılıkları sağlıklı kontrollerden alınan keratinositler ile p63 mutasyonu taşıyan hastalardan elde edilen üç çizgi (Mutantlar R204W, R279H ve R304W) arasında karşılaştırıldı. Şekil 2A, farklılaşma prosedürleri ve hücre morfolojisi değişikliklerine genel bir bakış gösterir. Mutant keratinositler çanak yüzeyinde düz kaldı ve kalabalık olmadı ya da 7.

PCA analizinde, kontrol hücresi çizgileri Şekil 1Bile karşılaştırıldığında açıkça bir farklılaşma deseni izlemiştir. Ancak, farklılaşma sırasında mutant hücrelerin gen ekspresyonu deseni büyük ölçüde çoğalan / farklılaşmamış hücrelerin benzer kalır. Üç mutant hattı arasında, R279 örnekleri ayırt edici pc1 ve PC2 boyunca bir ölçüde taşındı, onun farklılaşma r204W ve R304W ile karşılaştırıldığında, daha az bozulmuş olduğunu gösteren.

Kümeleme analizinde(Şekil 2C),kontrol hücrelerinde (küme 1) küçültülen genler R204W ve R279W'da kısmen indirgenmiş, ancak gen ekspresyonu R304W'de büyük ölçüde değişmemiştir. Bu genler büyük olasılıkla hücre çoğalmasında rol oynarlar, GO ek açıklamalarında gösterildiği gibi (Şekil 2D). Küme 2'de(Şekil 2C),genler ilk olarak kontrol hücrelerinde indüklendi ve daha sonra indirgendi. Epidermal farklılaşma ve keratinizasyon fonksiyonları bu küme genleri için yüksek oranda zenginleştirilmiş olduğundan, bu genler muhtemelen keratinosit farklılaşmasında rol oynarlar (Şekil 2D). Bu genler R204W ve R304W'da indüklenmezken, R279H hücrelerinde bu genler kontrol hücreleri kadar küçültülmemiştir.

Küme 3'teki genler sadece kontrol hücrelerinde farklılaşma sonunda indüklenmiştir(Şekil 2D). Bununla tutarlı olarak, bu genler epidermisin dış tabakasında rol alabilir, dış uyaranlara ve inflamasyona duyarlı oldukları gösterilmiştir(Şekil 2D). Bu genlerin ekspresyon deseni üç mutant hücre hattında da çok fazla değişmedi. Kontrol ve mutant hücreler arasındaki gen ekspresyonu desenleri görünür farklılıklar mutant hücrelerin bu in vitro farklılaşma modellerinde düzgün ayırt edemez göstermektedir.

Figure 1
Şekil 1: Keratinosit farklılaşması ve analizi. (A) Keratinosit farklılaşma protokolü ve hücre morfolojisi genel bakış. Ölçek çubuğu = 100 μm. (B) Diferansiyel kontrol keratinositlerinin prensip bileşen analizi. (C) Keratinosit farklılaşması sırasında ilk 500 değişken genin ısı haritası. Genler kmean kümeleme kullanılarak üç kümehalinde kümelenir. Yan tarafta her gen kümesi için temsili farklılaşma belirteç genleri gösterilir. (D) GO, yüksek değişkengen kümelerinde genler için aşırı temsil edilen fonksiyonların zenginleştirme analizi, tüm ifade edilen genlerin arka planına kıyasla (>10 sayısı). GOrilla zenginleştirme testi için kullanıldı. (E) Farklılaşma sırasında keratinosit farklılaşma belirteçleri Batı leke. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

Figure 2
Şekil 2: Kontrol ve p63 mutant keratinositlerin karşılaştırılması. (A) Keratinosit farklılaşma protokolü ve kontrol hücre morfolojisi ve p63 mutant keratinositlere genel bakış. Ölçek = 100 μm. (B) Farklılaşma kontrol ve hasta (R204W, R279H & R304W) keratinositler prensibi bileşen analizi. (C) Kontrol ve mutant keratinositler arasındaki diferansiyel genlerin ısı haritası. Genler kmean kümeleme kullanılarak üç kümehalinde kümelenir. Her gen kümesi için temsili farklılaşma belirteç genleri endikedir. (D) GO, yüksek değişkengen kümelerinde genler için aşırı temsil edilen fonksiyonların zenginleştirme analizi, tüm ifade edilen genlerin arka planına kıyasla (>10 sayısı). GOrilla zenginleştirme testi için kullanıldı. Bu rakamın daha büyük bir sürümünü görüntülemek için lütfen buraya tıklayın.

KGM bileşeni Stok Orta Birim
KBM 500 mL
Kalem/Strep 100.000 adet/mL 100 adet/mL 5 mL
BPE ~13 mg/mL 0.4% 2 mL
Etanolamin 0,1 M 0,1 mM 500 μL
o-fosfoetanolamin 0,1 M 0,1 mM 500 μL
Hidrokortizon 0.5 mg/mL 0,5 g/mL 500 μL
Insülin 5 mg/mL 5 μg/mL 500 μL
Egf 10 μg/mL 10 ng/mL 500 μL

Tablo 1. KGM-pro orta takviyeleri.

KGM bileşeni Stok Orta Birim
KBM 500 mL
Kalem/Strep 100.000 adet/mL 100 adet/mL 5 mL
Etanolamin 0,1 M 0,1 mM 500 μL
o-fosfoetanolamin 0,1 M 0,1 mM 500 μL

Tablo 2. KGM-diff orta takviyeleri.

Ek Kodlama Dosyaları: Folder_structure.txt; Generate_genome.txt; Map_fastq.txt; RNA_seq_kc_differentiation_patient.html; RNA_seq_kc_differentiation_wt.html; Sample_data_example.csv; TrimGalore.txt; ve Wig2bw.txt. Bu dosyayı indirmek için lütfen buraya tıklayınız.

Subscription Required. Please recommend JoVE to your librarian.

Discussion

Bu çalışma, RNA-seq analizleri kullanılarak insan keratinosit farklılaşması ve sonraki karakterizasyonu indükleyen bir yöntemi açıklamaktadır. Mevcut literatürde, insan keratinosit farklılaştırma üzerinde birçok çalışma, yüksek kalsiyum konsantrasyonu veya serum ile farklılaşmaiknayöntemleri olarak iki diğer yöntemler kullanın 2,3,23. Bir önceki rapor dikkatle bu üç farklı yöntem3 karşılaştırıldığında ve bu yöntemler keratinosit farklılaşma farklı biyoloji temsil edebilir gösterdi. Aynı raporda yazarlar, serum tarafından indüklenen diferansiye hücrelerin yüksek proliferatif potansiyele ve Sedef liciltte ifade edilen Ancak normal epidermiste ifade edilmeyen KRT16 ve SKALP/ PI3 gibi ekspres genlere sahip olduğunu ve bu nedenle sedef hastalığını incelemek için serum kaynaklı farklılaşma modelinin kullanılabileceğini göstermiştir. Buna karşılık, kontak tinhibisyonu artı büyüme faktörü dışlama yöntemi KRT1 ve KRT10 neden olabilir, normalde epidermis ifade edilir ve normal keratinosit farklılaşması benzer. Yüksek kalsiyum indüksiyonu, serum indüksiyonu ile kontak tinhibisyonu arasında en az spesifik yöntem olarak gen ekspresyonu profiline yol açar. Farklılaşma sırasında RNA-seq analizleri kullanılarak yapılan analizler, kontak inhibisyonu ve büyüme faktörü dışlama yönteminin epidermal farklılaşma için beklenene göre benzer gen ekspresyonu profillerine sahip farklılaştırılmış keratinositlerle sonuçverdiğini doğruladı (örneğin, çoğalan hücrelerde KRT5, erken farklılaşma indüksiyonunda KRT1 ve KRT10, geç farklılaşmada LOR ve FLG; Şekil 1C).

Bu bulgular, bu farklılaşma tekniğinin keratinosit farklılaşmasını incelemek için kullanımı kolay ve güvenilir bir yöntem olduğunu göstermektedir. Ayrıca, p63 mutant keratinositlerden elde edilen keratinositler üzerinde yapılan bu çalışma, hastalığın alabildiği durumlarda etkilenen keratinosit farklılaşmasını incelemek için yöntemin kullanılabileceğini göstermektedir. Bu in vitro farklılaştırma yaklaşımında Lonza'dan satın alınan KGM, tutarlı sonuçlar vermiştir. Termo Fisher Scientific'in keratinosit serumsuz orta (KSFM) gibi benzer bileşimi olan diğer epidermal ortam da büyük olasılıkla uygundur, ancak bunun test edilmesi gerekmektedir.

Bu yöntemin bazı sınırlamaları olduğu unutulmamalıdır. Temas inhibisyonuna dayandığı için hücre yoğunluğunun birleşmesi gerekir. P63 mutant keratinositler ile yapılan deneylerde, hücrelerin biraraya gelmesi için daha yüksek bir başlangıç tohumlama yoğunluğu gereklidir. Buna ek olarak, hücreler aynı hızda büyümez, farklılaşma bazen farklı günlerde indüklenen gerekir. Deneyler ilerlerken bu hususlar göz önünde bulundurulmalıdır. Hücrelerin aynı anda indüklmesi gereken deneysel ortamlarda, serum eklenmesi, yüksek kalsiyum konsantrasyonları ve epidermal büyüme faktörüreseptörlerinininhibisyonu 2,3dahil olmak üzere diğer farklılaşma yöntemleri düşünülmelidir. Yine de, tüm in vitro diferansiyasyon yöntemleri artıları ve eksileri var, muhtemelen in vivo cilt gelişimi sırasında mevcut kısmi farklılaşma indüksiyon sinyalleri temsil olarak2.

Bu farklı yöntemleri karşılaştıran kapsamlı bir moleküler analiz son derece bilgilendirici olacak ve biyolojik süreçler üzerinde farklı çalışmalar için en uygun yöntem seçimi talimat verebilir. Buna ek olarak, transkriptomik düzeyde ölçülen değişikliklerin doğrulanması, örneğin batı lekeleme veya proteomik analizler yoluyla son derece tavsiye edilir. Ayrıca in vitro veriler dikkatli kullanılmalıdır ve bu modellerden elde edilen sonuçlar tercihen insan cilt gelişiminde in vivo olarak doğrulanmalıdır.

RNA çıkarma ve RNA-seq kütüphane hazırlama temel ilkeleri de dahaönce24,25,26,27,28açıklanmıştır . Bu protokolde, RNA çıkarma ve RNA-seq kitaplık hazırlama yordamları ticari olarak kullanılabilen kitlerin iş akışlarına dayanmaktadır. Prensip olarak, RNA kalitesi iyiyse RNA-seq analizleri üzerinde farklı yöntemlerin RNA-seq analizleri üzerinde büyük bir etkisi olmamalıdır. RNA kalitesinin düşük olduğu durumlarda (örneğin, RNA formalin-sabit parafin gömülü dokular çıkarıldığında), RNA-seq hala yapılabilir; ancak, RNA parçalanma adımı ayarlanmalıdır. RNA-seq kitaplık hazırlığı ribozomal RNA (rRNA) çıkarılmasından sıralama için cDNA kitaplığı yapısına kadar kapsar. Ana prosedürler çeşitli kitleri kullanarak, ya da bireysel enzimler ve ev yapımı tamponlar kullanılarak yapılabilir29,30,31,32. RNA-seq analizi farklı temel ilkeler kullanılarak yapılırsa (örneğin, ribozomal RNA tükenmesi veya poliT oligolarına hibridizasyon ile poliA mRNA zenginleştirme), RNA-seq analizlerinin sonucu farklı olabilir. Ayrıca protokol, DNA oligolarının insan, fare ve sıçan rRNA'ya hibridizasyonu ile ribozomal RNA çıkarılmasını kullanır. Farklı türlerle çalışırken alternatif bir oligo seti kullanılmalıdır.

Oluşturulan kitaplığın sıralaması parçaların her iki ucunda (eşleştirilmiş uç) veya parçanın bir ucunda (tek uç) gerçekleştirilebilir. Genel olarak, eşleştirilmiş uç sıralaması büyük ölçüde eşlenebilirliği artırır ve transkript türevleri ile ilgili daha fazla bilgi sağlar. Ancak, burada açıklandığı gibi nispeten basit bir diferansiyel gen analizi için, tek uç sıralama da yeterli bilgi sağlayabilir.

Veri analizinin önemli bir adımı için, fastq dosyalarının kalite kontrolü için nispeten basit bir yöntem tanımlanır ve ardından genoma haritalama okumaları yapılır. Sağlanan bash kodu ideal ölçeklenebilirlik olmamasına rağmen, şeffaflık avantajına sahiptir. Yazılım seçimi, hangi adımları gerçekleştirir ve veri önişleme sırasında kullanılan genom sürümü, tekrarlanabilirlik için gerekli olan iyi belgelenmesi gereken önemsiz olmayan adımlardır.

Daha gelişmiş kullanıcılar için, otomatik bir boru hattı fastQ dosya kalite kontrolü, bağdaştırıcı kırpma ve haritalama, örneğin ARMOR yılan iş akışı33 veya 'RNA-seq' Snakemake iş akışı van Heeringen lab34adımları gerçekleştirmek için kullanılabilir. Ancak, bu tamamen otomatik boru hatları daha az saydam ve değiştirmek daha zordur. Bu araçları kullanırken, bu otomatik boru hatları içindeki işlevleri anlamak için hayati önem taşımaktadır. Son olarak, protokol zaman içinde değişken gen ekspresyonuna vurgu ile RNA-seq veri analizi içerir. Değişken gen ekspresyonu, farklılaşma gibi birden fazla zaman puanıolan süreçlere bakıldığında çift yönlü diferansiyel testlere tercih edilir. Sonuç olarak, bu analiz boru hattı RNAseq verilerinin biyoinformatik analizi için ilgili araçları tanır. Hem normal hem de hastalıklı koşullarda keratinosit farklılaşmasını iyice değerlendirebilen araçlar içerir.

Subscription Required. Please recommend JoVE to your librarian.

Disclosures

Yazarların açıklayacak bir şeyi yok.

Acknowledgments

Bu araştırma Hollanda Bilimsel Araştırma Örgütü (NWO/ALW/MEERVOUD/836.12.010, H.Z.) tarafından desteklenmiştir. (NWO/ALW/Açık Rekabet/ALWOP 376, H.Z., J.G.A.S.); Radboud Üniversitesi bursu (H.Z.); ve Çin Burs Konseyi hibe 201406330059 (J.Q.).

Materials

Name Company Catalog Number Comments
Bioanalyzer 2100 Agilent G2929BA
Bovine pituitary extract (BPE) Lonza Part of the bulletKit
CFX96 Real-Time system Bio-Rad qPCR machine
Dulbecco's Phosphate-Buffered Saline (DPBS) Sigma-Aldrich D8537
Epidermal Growth Factor (EGF) Lonza Part of the bulletKit
Ethanolamine >= 98% Sigma-Aldrich E9508
High Sensitivity DNA chips Agilent 5067-4626
Hydrocortison Lonza Part of the bulletKit
Insulin Lonza Part of the bulletKit
iQ SYBR Green Kit BioRad 170-8886
iScript cDNA synthesis Bio rad 1708890
KAPA Library Quant Kit Roche 07960255001 Low concentration measure kit
KAPA RNA HyperPrep Kit with RiboErase Roche KK8540 RNAseq kit
KGM Gold Keratinocyte Growth Medium BulletKit Lonza 192060
Nanodrop deNovix DS-11 FX (model) Nanodrop and Qbit for DNA and RNA measurements
NEXTflex DNA barcodes -24 Illumnia NOVA-514103 6 bp long primers
Penicillin-Streptomycin Gibco 15140122
RNA Pico Chip Agilent 5067-1513

DOWNLOAD MATERIALS LIST

References

  1. Hardman, J. A. Skin equivalents for studying the secrets of skin: from development to disease. British Journal of Dermatology. 173 (2), 320-321 (2015).
  2. Borowiec, A. S., Delcourt, P., Dewailly, E., Bidaux, G. Optimal differentiation of in vitro keratinocytes requires multifactorial external control. PLoS One. 8 (10), 77507 (2013).
  3. Van Ruissen, F., et al. Induction of normal and psoriatic phenotypes in submerged keratinocyte cultures. Journal of Cellular Physiology. 168 (2), 442-452 (1996).
  4. Ojeh, N., Pastar, I., Tomic-Canic, M., Stojadinovic, O. Stem Cells in Skin Regeneration, Wound Healing, and Their Clinical Applications. International Journal of Molecular Sciences. 16 (10), 25476-25501 (2015).
  5. Ali, N., et al. Skin equivalents: skin from reconstructions as models to study skin development and diseases. British Journal of Dermatology. 173 (2), 391-403 (2015).
  6. Luis, N. M., et al. Regulation of human epidermal stem cell proliferation and senescence requires polycomb- dependent and -independent functions of Cbx4. Cell Stem Cell. 9 (3), 233-246 (2011).
  7. Mulder, K. W., et al. Diverse epigenetic strategies interact to control epidermal differentiation. Nature Cell Biology. 14 (7), 753-763 (2012).
  8. Poumay, Y., Pittelkow, M. R. Cell density and culture factors regulate keratinocyte commitment to differentiation and expression of suprabasal K1/K10 keratins. Journal of Investigative Dermatology. 104 (2), 271-276 (1995).
  9. Kouwenhoven, E. N., et al. Transcription factor p63 bookmarks and regulates dynamic enhancers during epidermal differentiation. EMBO Rep. 16 (7), 863-878 (2015).
  10. Qu, J., et al. Mutant p63 Affects Epidermal Cell Identity through Rewiring the Enhancer Landscape. Cell Reports. 25 (12), 3490-3503 (2018).
  11. Brunner, H. G., Hamel, B. C., Van Bokhoven, H. The p63 gene in EEC and other syndromes. Journal of Medical Genetics. 39 (6), 377-381 (2002).
  12. Browne, G., et al. Differential altered stability and transcriptional activity of DeltaNp63 mutants in distinct ectodermal dysplasias. Journal of Cell Science. 124, Pt 13 2200-2207 (2011).
  13. KGM Gold Keratinocyte Growth Medium BulletKit. , Available from: https://bioscience.lonza.com/Lonza_bs/CH/en/Culture-Media-and-Reagents/p/000000000000192060/KGM-Gold-Keratinocyte-Growth-Medium-BulletKit (2019).
  14. Doyle, K. aM. Protocols and applications guide. , Promega Corporation. Madison. (1996).
  15. Roche. Hyperprep Kit. , Available from: https://sequencing.roche.com/en-us/products-solutions/by-category/Library-preparation/rna-library-preparation/kapa-rna-hyperprep-kit-with-riboerasehmr.html (2019).
  16. Felix Krueger. TrimGalore. , Available from: https://github.com/FelixKrueger/TrimGalore (2019).
  17. Dobin, A., et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 29 (1), 15-21 (2013).
  18. Li, H., et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 25 (16), 2078-2079 (2009).
  19. Ryan, D. Convert chromosome names in a bigWig file. , Available from: https://gist.github.com/dpryan79/39c70b4429dd4559d88fb079b8669721 (2019).
  20. Love, M. I., Huber, W., Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550 (2014).
  21. Eden, E., Navon, R., Steinfeld, I., Lipson, D., Yakhini, Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 10, 48 (2009).
  22. Yu, G., Wang, L. G., Han, Y., He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 16 (5), 284-287 (2012).
  23. Kretz, M., et al. Control of somatic tissue differentiation by the long non-coding RNA TINCR. Nature. 493 (7431), 231-235 (2013).
  24. Mullegama, S. V., et al. Nucleic Acid Extraction from Human Biological Samples. Methods in Molecular Biology. 1897, 359-383 (2019).
  25. Ma, F., et al. A comparison between whole transcript and 3' RNA sequencing methods using Kapa and Lexogen library preparation methods. BMC Genomics. 20 (1), 9 (2019).
  26. Chao, H. P., et al. Systematic evaluation of RNA-Seq preparation protocol performance. BMC Genomics. 20 (1), 571 (2019).
  27. Podnar, J., Deiderick, H., Huerta, G., Hunicke-Smith, S. Next-Generation Sequencing RNA-Seq Library Construction. Current Protocols in Molecular Biology. 106, 21 (2014).
  28. Haque, A., Engel, J., Teichmann, S. A., Lonnberg, T. A practical guide to single-cell RNA-sequencing for biomedical research and clinical applications. Genome Medicine. 9 (1), 75 (2017).
  29. Oyola, S. O., et al. Optimizing Illumina next-generation sequencing library preparation for extremely AT-biased genomes. BMC Genomics. 13, 1 (2012).
  30. Quail, M. A., et al. Optimal enzymes for amplifying sequencing libraries. Nature Methods. 9 (1), 10-11 (2011).
  31. Quail, M. A., et al. A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers. BMC Genomics. 13, 341 (2012).
  32. Ross, M. G., et al. Characterizing and measuring bias in sequence data. Genome Biology. 14 (5), 51 (2013).
  33. Soneson, C. ARMOR. , Available from: https://github.com/csoneson/ARMOR (2019).
  34. Sande, S. F. snakemake-workflows. , Available from: https://github.com/vanheeringen-lab/snakemake-workflows (2019).

Tags

Genetik Sayı 159 insan primer keratinositler 2D batık kültür in vitro farklılaşma RNA-seq biyoinformatik analizi p63
RNA-Seq Analizi ile İnsan Primer Keratinositlerin In Vitro Farklılaşmasının Karakterizasyonu
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Smits, J. G. A., Qu, J., Niehues,More

Smits, J. G. A., Qu, J., Niehues, H., Zhou, H. Characterization of In Vitro Differentiation of Human Primary Keratinocytes by RNA-Seq Analysis. J. Vis. Exp. (159), e60905, doi:10.3791/60905 (2020).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter