$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Sıçan metil-Seq platformu başarılı bir uygulama çeşitli ölçütlere bağlıdır. Şekil 1 çalışmanın genel iş akışı gösterilmektedir ve ileri taşımadan önce gerekli olan belirli kalite kontrol (QC) adımlar vurgular. Dikkat edilmesi gereken ilk faktör biridir methylome arasında oluşan epigenetik değişiklikleri büyüklüğünü belirlemek hayvan modeli ve stres rejimi sağlamlık. Hayvan çalışmalarımız corticosterone (CORT) maruz kalma DNA metilasyonu19,20değişikliklere neden olabilir bizim önceki gözlem esas beri bizim kronik değişken stres (CVS) rejimi üretmek için yeterli rigor gerekiyordu yüksek plazma CORT düzeyleri ile Rats vurguladı. Tipik bir haftalık CVS rejimi Tablo 1 ' de gösterilen ve günlük stres sabah, öğleden sonra oluşuyordu ve gecede sürekli habituation önlemek için değiştirilir ve stres tepkisi azalmış. 3 haftalık rejim stresli hayvanlar ortalama plazma CORT önemli ölçüde yüksek seviyede sergilenen [gün 4 – 21, kontrol: 32.7 3.7 ng/mL, stres: 103.0 11,9 ng/mL (ortalama SEM), P 2.2 x 10-4, Şekil 2A=] Bu olarak vurgusuz, hayvan kontrol. Sürekli olarak, bu hayvanlar da yükseltilmiş durumda büyük anksiyete gibi davranış artı labirent (EPM), EPM ve kollarını açarak (2B rakam) zamanında daha az kapalı kollarında geçirdiğim önemli ölçüde daha fazla gösterildiği gibi gösterdi. Bu sonuçlar CVS pozlama önemli endokrin ve bize bu değişiklikler belirli DNA metilasyonu imzaları ile ilişkili olup olmadığını araştırmak için önde gelen davranış değişiklikleri için yol göstermektedir.
Biz metil-Seq Kütüphane başarılı yapımı için çok önemli birkaç denetim noktaları vurgulamak. Yeterli miktarda DNA ile başlayan gerekli, sonication, birden fazla yıkama/arıtma, hedef zenginleştirme ve bisülfit dönüşüm adımları sırayla DNA bitmiş kütüphanede bulunan miktarı azaltın. Her ne kadar birçok PCR güçlendirme adımları DNA şablon kaybı hafifletmek, aşırı PCR döngüsü numaraları yüksek yinelenen okuma neden olabilir. Sıçan metil-Seq çalışmada için kan gDNA fare başına 2 g kullanıldı. Biz metil-Seq kitaplıkları DNA miktarı düşük 500 olarak başlayan ile yapılabilir notu ng. Sonraki sıralama için kitaplıklarına yetersiz miktarda üreten riski olsa daha küçük malzeme başlayan (floresan aktif hücre sıralama) FACS tarafından izole DNA kitaplıkları oluşturmak veya yumruk, iğne kullanıcılara izin verir. QC Elektroforez örneği DNA molekül ağırlığı, miktar ve derişim sağlar bir bioanalyzer üzerinde 1 m tarafından gerçekleştirilir. Bioanalyzer kullanımını gerektiren üç kritik adım vardır: 1) sonication adım yeterli DNA'sı (~ 170 bp, kırmızı, Şekil 3); kesme emin olmak için takip 2) bir vardiyada yamultulmuş DNA (~ 200 bp, mavi, Şekil 3) onların sonraki amplifikasyon PCR tarafından; emin olmak için Ortalama büyüklüğü bağdaştırıcısı ligasyonu adım takip simgesiyle ve 3) miktar ve sıralama için Kütüphane boyutu sağlamak için son Kütüphane arıtma adım takip.
R-paketleri BSSeq ve BSmooth Bioconductor içinde veri18sıralama bisülfit analiz etmek için kullanılmıştır. Onlar araçları ve yöntemleri, kalite kontrol, gerçekleştirme sırası okur hizalamak için içerir ve differentially tanımlayan metillenmiş bölgeleri (DMRs). BSmooth yazılım papyon 2.016,17 CpG düzeyi ölçüm özetleri, hizalama genomik sıralarını bisülfit dönüştürülmüş ham giriş okur tarafından elde etmek için bir iç sıra aligner olarak çağırır. Hizalanmış okuma sonra sistematik sıralama ve aşağı akım analizleri çarpık Bankası-arama hataları tanımlamak için aramak sıkı kalite kontrol işlemleri yoluyla filtre uygulanır. Araziler bir dizi görsel olarak filtreleme bu süreç içinde yardım etmek için oluşturulur. Sıralama ölçütleri de belge numarası hizalanmış okur, % hedef ve diğerleri (Tablo 2) CpG kapsama başına gibi ilgili bilgileri için oluşturulur. Bir kez veri süzülür, yumuşatma/normalleştirme algoritması gerçekleştirildiğinde, nerede her CpG tüm QC dayalı bir tahmini metilasyonu değer atanmış her örnekten okur ve CpGs metilasyon daha doğru arama emin olmak için komşu tahminleri durumu bile sıra kapsamı düşük olduğu durumlarda. Bu değer her CpG sitede metilasyonu olasılığının Düzgünleştirilmiş bir tahmin sağlar. Sıralaması için en azından en önemli ölçüde farklı genomik bölgelerine ve iki tedavi grupları arasında her örneğinin düzgünleştirilmiş metilasyonu tahminleri ortalaması karşılaştırarak, DMRs listesi oluşturulur (Tablo 3).
DMR vurguladı ve vurgusuz gruplar arasında fare MHC gen Rt1-m4, organizatörü ile yer oldu üst vurgusuz hayvanların (Şekil 4A) daha tüm CpGs arasında metilasyon katlara sergilenmesi hayvanlar vurguladı. Metil-Seq platformu ve veri analizi başarılı uygulanması onaylamak için astar DMR ve kan DNA metilasyonu düzeyleri ve vurgusuz stresli hayvanlar (metil-Seq tarafından sıralı 8 ve 8 değil sıralı) tüm kohort karşı tasarlanmıştır bisülfit pyrosequencing tarafından değerlendirildi. Sonuçları 10 dışarı-in denetlesinler 12 CpGs arasında DNA metilasyonu önemli artış göstermektedir (% metilasyonu, P 5.1-10,4 değişiklik < 0,037, Şekil 4B). KEGG yolu analiz tüm sözde önemli DMRs stres ile ilgili yollar tanımlamak için gerçekleştirildi. Sürekli olarak, DMR ilişkili yolları hastalıkları kronik stres maruz kalma, diyabet, kardiyovasküler hastalık ve kanser (Tablo 4) gibi ilişkili karıştığı. 21 , 22 , 23 , her hayvan için ortalama 3 haftalık CORT düzeylere KSY-10 düzeylerinde metilasyonu strese maruz kalma derecesi ve epigenetik veri arasındaki ilişkiyi göstermek için karşılaştırıldı. Sonuçlar gösterdi endokrin ve metilasyonu veri arasında mütevazı bir korelasyon (R2= 0.54, P = 0.001, Şekil 5).

Şekil 1: genel olarak şematik iş akışı sıçan metil-Seq platformu için. Bir g kandan ayıklanan genomik DNA'ın vurguladı ve kontrol fareler ilk sıralama, analiz ve hedef tanımlama için metil-Seq kitaplıkları oluşturmak için işlenmiş. Başka bir 100 ng DNA'ın bisülfit pyrosequencing tarafından tanımlanan epigenetik hedefleri bağımsız doğrulamak için kullanılır. Bu rakam daha büyük bir versiyonunu görüntülemek için buraya tıklayınız.

Şekil 2: kronik değişken stres (CVS) maruz açar sıçanlarda endokrin ve davranış değişiklikleri için. (A)corticosterone (CORT) birden fazla örnekleme göstermek sağlamlık 3 hafta CVS rejimi. Kan örnekleri günlük stres alay önce sabahları toplanmıştır. (B) vurguladı hayvanlar kapalı silah ve yükseltilmiş kollarını açarak zamanında artı labirent (EPM) daha az daha fazla zaman geçirdim. Öğretici Boxplots her hayvan için veri noktası ile gösterilir. Öğrenci T-testi için istatistiksel anlamlılık gerçekleştirildi. * P < 0,05, ** P < 0,01, ve *** P < 0,001. Bu rakam daha büyük bir versiyonunu görüntülemek için buraya tıklayınız.

Şekil 3: Nefelometri yamultulmuş ve bağdaştırıcı bakmaksızın sıçan DNA üzerinde bir bioanalyzer. Kırmızı ve mavi eğrileri miktar ve genomik DNA (kırmızı) bir izotermal sonicator ve bağdaştırıcı ligasyonu, sırasıyla kesme takip boyutunu gösterir. Her satır bir örnek ve kırmızı temsil eder ve mavi eğrileri DNA ' (bitiş-tamir, 3'-adenylation ve örnek temizleme) birkaç adım sırasında her iki kaybı yansıtmak ve bağdaştırıcıları ligasyonu nedeniyle bp boyutunda artırmak. Keskin doruklarına 25 bp ve 1500 bp yükleme arabelleğe eklendi standart işaretleri vardır. Bu rakam daha büyük bir versiyonunu görüntülemek için buraya tıklayınız.

Şekil 4: epigenetik değişiklikleri CVS kaynaklı sıçan metil-Seq. tarafından tespit Sıçan metil-Seq veri analizi(a)karıştığı gen Rt1m4 organizatörü differentially metillenmiş bölgesine (DMR) stresli (kırmızı) ve denetim (mavi) fareler arasında. Grafik çıktı Rt1m4 DMR (pembe gölgeli bölgesi) için her CpG (dikey gri çizgi), dört örnekleri her grup (kırmızı veya mavi çizgiler) ve her hayvan (kırmızı veya mavi nokta) için % metilasyonu düzeylerini görüntüler. (B) on iki CpGs içinde DMR bisülfit pyrosequencing tarafından onaylanmış. Çubuk grafikler SEM demek olarak temsil edilir ve bir öğrenci T-testi için istatistiksel anlamlılık gerçekleştirildi. * P < 0,05. Bu rakam daha büyük bir versiyonunu görüntülemek için buraya tıklayınız.

Şekil 5: doğrusal regresyon çözümlemesi gösterdi mütevazı bir korelasyon arasında % DNA metilasyonu KSY-10 / Rt1m4 ve 3 hafta plazma CORT düzeyleri vurguladı her ikisi anlamına gelir ve kontrol hayvanlar (N = 16). Stresli hayvanlar verilerden kırmızı daireler tarafından temsil edilen. Bu rakam daha büyük bir versiyonunu görüntülemek için buraya tıklayınız.
| Hafta | Gün 1 | 2. gün | 3. gün | 4. gün | Gün 5 | Gün 6 | Gün 7 |
| AM | Kısıtlama | Yüzmek | Soğuk Oda | Yüzmek | Kısıtlama | Shaker | Yüzmek |
| PM | Shaker | Kafes Tilt | Kısıtlama | Shaker | Soğuk Oda | Kısıtlama | Soğuk Oda |
| Bir gecede | Gıda kısıtlamak | Odada mevcut yataklar kullanılıyorsa ıslak | Yalıtım | Işık | Kalabalık | Işık | Odada mevcut yataklar kullanılıyorsa ıslak |
Tablo 1: Bir normal haftalık zamanlama, kronik değişken stres rejimi (CVS).
| Sıralama ölçütleri | Stres1 | Kontrol1 |
| (n = 4) | (n = 4) |
| Son Okunma (başı) eşleştirilmiş | 89,290,397 | 80,165,674 |
| Benzersiz olarak eşlenen eşleştirilmiş sonunda (UMPER) okur | 39,200,255 | 35,013,406 |
| Hizalama oranı/eşleme verimlilik (UMPER / başı) | % 44 | % 44 |
| Yinelenen okuma (UMPER %) | % 73 | % 65 |
| Deduplicated UMPER | 10,481,031 | 12,306,018 |
| Ortalama okuma derinlik kapsama (x) (ARDC) | 6 x | 6 x |
| CpGs (N) | 12,056,878 | 12,056,878 |
| ARDC (x) CpGs | 2 x | 2 x |
| En az 10 okuma (N) ile CpGs | 481,383 | 595,850 |
| ARDC (X) CpGs ile en az 10 defa okundu | 19 | 19 |
| Hedef CpGs (sonda hedef bölgeleriyle tam örtüşme) | 1,923,872 | 2,007,638 |
| Hedef ARDC (x) CpGs | 7 x | 8 x |
| Hedef CpGs üzerinde en az 10 okuma (N) ile | 428,249 | 531,419 |
| Hedef ARDC (x) CpGs en az 10 ile okur | 18 x | 18 x |
| Hedef (hedef bölgeleriyle Probe 1 veya daha fazla baz çifti örtüşme ile PER) (UMPER) | 8,277,715 | 9,369,523 |
| Hedef (Deduplicated UMPER) Toplam % | % 78 | % 77 |
| Hedefe (Toplam üsleri eşlenen) Mb | 125 mb | 128 mb |
| Açık hedef ortalama okuma derinlik kapsama (x) (ARDC) | 9 x | 10 x |
|
1 Her grup bireylerde genelinde ortalamalar göre sıralama ölçütleri | | |
Tablo 2: fare metil-Seq platformu elde edilen ölçümler sıralanıyor.
| Chr | Başlat | bitiş | Gen | mesafe | areaStat | meanDiff | stres | Denetim | Yön |
| chr20 | 1,644,246 | 1,644,390 | RT1-M4 | in_gene | 93.03 | 0,22 | 0,33 | 0,11 | kazanç |
| chr5 | 160,361,352 | 160,361,564 | LOC690911 | in_gene | -70.75 | -0.19 | 0,72 | 0,91 | kaybı |
| chr3 | 61,138,281 | 61,138,330 | RGD1564319 | 265569 | 61.79 | 0.21 | 0,94 | 0,72 | kazanç |
| chr2 | 143,064,811 | 143,065,010 | Ufm1 | 8569 | -59.48 | -0.11 | 0,13 | 0,24 | kaybı |
| chr7 | 30,764,111 | 30,764,284 | Ntn4 | in_gene | 57.04 | 0.21 | 0,94 | 0,73 | kazanç |
| chr17 | 12,469,112 | 12,469,218 | Idnk | 41996 | -50.91 | -0.13 | 0,74 | 0,88 | kaybı |
| chr7 | 47,101,725 | 47,101,930 | Pato | in_gene | -50.54 | -0.12 | 0.64 | 0,76 | kaybı |
| chr5 | 76,111,248 | 76,111,822 | Txndc8 | 151703 | -50.38 | -0.11 | 0.85 | 0,96 | kaybı |
| chr11 | 80,640,132 | 80,640,356 | Dgkg | in_gene | -50.07 | -0.16 | 0,73 | 0,89 | kaybı |
| chr8 | 71,759,248 | 71,759,411 | Mir190 | 210226 | -47.84 | -0.17 | 0,58 | 0,75 | kaybı |
Tablo 3: en iyi 10 differentially bölgeleri metillenmiş. Her DMR çıkış tablosu için sağ sütun soldan gösterir: kromozom konumu (chr), koordinatları (başlangıç/bitiş), gen adı, transkripsiyon mesafeden sitesini başlatın, vurguladı ve kontrol grupları (areaStat) arasındaki fark alanı istatistikleri, demek Ortalama metilasyon düzeyleri için her DMR arasında fark metilasyonu (meanDiff), vurguladı ve kontrol grubu (stres/denetimi) ve metilasyonu yönünü denetimlerden değiştirin.
| KEGG yol şartları | Gen sayısı | % | P-değeri | Benjamini |
| Diyabet |
| Tip II diyabet | 12 | 0,1 | 3.6 x 10-4
| 9,8 10-3 x
|
| Kardiyovasküler hastalık |
| Vasküler düz kas kasılma | 18 | 0,1 | 1.6 x 10-3
| 3.6 x 10-2
|
| Aritmojenik Sağ ventrikül kardiyomiyopati (R.V.C) | 13 | 0,1 | 4.0 x 10-3
| 7.1 x 10-2
|
| Dilate kardiyomiyopati | 14 | 0,1 | 7,6 10-3 x
| 1.2 x 10-1
|
| Nöron işlevi |
| Uzun vadeli kullanılmasının muhtemelen | 11 | 0,1 | 1.5 x 10-2
| 1.4 x 10-1
|
| Sinyal verme |
| MAPK sinyal yolu | 35 | 0,2 | 2.4 x 10-4
| 10-3 x 9,9
|
| Kalsiyum sinyal yolu | 22 | 0,1 | 1.2 x 10-2
| 1.4 x 10-1
|
| Kemokin sinyal yolu | 21 | 0,1 | 1.2 x 10-2
| 1.3 x 10-1
|
| Kanser |
| Yolları kanser | 42 | 0,3 | 4.1 x 10-5
| 3.4 x 10-3
|
| Tümörü | 15 | 0,1 | 4.4 x 10-5
| 2.4 x 10-3
|
| Sigara-küçük hücreli akciğer kanseri | 10 | 0,1 | 7,9 10-3 x
| 10-1 x 1.1
|
| Kolorektal kanser | 13 | 0,1 | 8,4 10-3 x
| 10-1 x 1.1
|
| Kronik miyeloid lösemi | 12 | 0,1 | 1.2 x 10-2
| 1.3 x 10-1
|
Tablo 4: KEGG yolu analizini sıçan metil-Seq. tespit DMRs