$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Açıklanan Capture Hi-C protokolü, dört tabanlı bir kesici (DpnII) kullanılarak genom çapında 3C şablonunun hazırlanmasına dayanmaktadır. İlgilenilen genomik bölge boyunca ligasyon fragmanlarının daha sonraki zenginleştirilmesi, bir dizi döşeme RNA probunun hibridizasyonu ve bu çalışmada kullanılan hedef zenginleştirme sistemine göre streptavidin bazlı yakalanması ile elde edilir (Şekil 1). Biyotinile RNA probları, DNA problarına kıyasla hedeflerine daha sıkı bağlanma afinitesi gösterdikleri için seçildi52,60. Yakalanan kitaplıklar daha sonra dizine alınır ve çoğullanmış yüksek verimli sıralama için havuza alınır. Hi-C verilerini yakalamak, yüksek çözünürlüklü Hi-C etkileşim haritaları olarak görselleştirilebilir, ancak aynı zamanda yakalanan bölgenin tamamındaki destekleyiciler veya arttırıcılar gibi daha küçük dizilerin etkileşimlerini özel olarak görselleştirmek için 4C benzeri tek görüş noktası temas haritaları olarak da görselleştirilebilir. Protokolün iş akışı Şekil 4'te gösterilmiştir. Ön sıralama kalite kontrolleri Şekil 2'de gösterilmiştir ve 3C şablonunun uygun sindiriminin ve yeniden ligasyonunun değerlendirilmesini ve protokolün farklı adımlarında etkili kesilmesini ve saflaştırılmasını içerir. Kesilen 3C şablon DNA'nın 150 ila 700 bp arasında çalışması bekleniyor ve parçaların >2 kb zenginleştirilmesi tespit edilmemelidir. Aşağıdaki adımlarda, önce kesmeden sonra, daha sonra yakalama öncesi ve yakalama sonrası PCR'lerden sonra birkaç boncuk bazlı DNA temizleme ve boyut seçimi adımı gerçekleştirilir. Temizlenen kütüphaneler, yüksek hassasiyetli bir DNA biyoanalizöründe görselleştirildiği gibi farklı bir parça zenginleştirme profili gösterir (Şekil 2). Ortalama parça boyutu, bağdaştırıcıların bağlanması, dizileme ve indeksleme primerleri nedeniyle kütüphane hazırlığı boyunca artar. Sıralama sonrası kalite kontrolleri Hi-C Pro ile elde edilir ve Şekil 3'te gösterilmiştir. 3C benzeri veri işleme ve analizi için birçok farklı biyoinformatik yazılım uygulaması önerilmiştir. Bunlar arasında, HiC-Pro boru hattı, ham sıralama verilerinin çeşitli çözünürlüklerde son temas haritalarına işlenmesine izin veren en popüler çözümlerden biridir55. HiC-Pro, referans genom üzerindeki dizileme okumalarını hizalamak için iki adımlı bir haritalama stratejisi kullanır. 3C ürünleri daha sonra bilgilendirici olmayan temas çiftlerini kaldırmak ve temas haritalarını oluşturmak için yeniden yapılandırılır ve filtrelenir. Ek olarak, alele özgü analiz yapmak ve iki ebeveyn alelinden gelen temasları farklı temas haritalarında ayırmak için bilinen polimorfizmlerin bir listesini kullanabilir. Daha yakın zamanlarda, HiC-Pro, nf-core çerçevesine (nf-core-hic) dahil edildi ve genişletildi, bu da yüksek oranda ölçeklenebilir ve tekrarlanabilir topluluk odaklı bir boru hattı61,62 sağladı.
Xic faresini yakalamak için, X kromozomunun 3 Mb'sini döşeyen 28.913 RNA probu dizisi tasarlandı. Bu bölge, XCI'daki kilit oyuncuyu, uzun kodlamayan gen Xist'i ve bilinen ~ 800 kb düzenleyici manzarasını içerir (Şekil 5). Bu ~ 800 kb bölgesi iki TAD'a bölünmüştür: biri Xist promotörü ve bilinen pozitif düzenleyicileri (yani, kodlamayan transkriptler Ftx, Jpx ve Xert ve protein kodlayan gen Rnf12) ve Xist'in negatif cis-düzenleyicilerini kapsayan komşu TAD (yani, antisens transkript Tsix, geliştirici element Xite ve kodlamayan transkript Linx) (inceleme için44, 45).
Tanımlanan Capture Hi-C protokolünü Xic'e uygulayarak, bu lokusun topolojik organizasyonu benzeri görülmemiş bir çözünürlükte elde edildi (Şekil 6 ve Şekil 7). Bu, özellikle Capture Hi-C profili ile daha önce yayınlanmış 5C47 (Şekil 6 ve Şekil 7; Ek Tablo 1) ve Hi-C61 (Şekil 6 ve Şekil 7; Ek Tablo 1) Profil. Örneğin, alt-TAD yapıları daha belirgindir – Xist promotörünü (Xist-TAD) içeren TAD açıkça iki küçük alana bölünmüştür (Şekil 6A, mavi ok ucu). Önceden, bu, yalıtım puanı algoritması kullanılarak bu bölgede bir sınırın tespit edilmesine rağmen, yalnızca 5C profilinden (Şekil 6B) görsel olarak "tahmin edilebilir". Benzer şekilde, Capture Hi-C profilinin çözünürlüğü, Tsix lokusunun (Tsix-TAD) destekleyicisini içeren komşu TAD'de (Şekil 6A, B) iki küçük etki alanının tanımlanmasına izin verir; bu daha önce 5C ile elde edilememişti (Şekil 6B). Not olarak, Capture Hi-C ve 5C verilerinden yalıtım skoru ile belirlenen topolojik sınırlar genellikle biraz farklı konumlarda ve farklı göreceli güçlerde tespit edilir.
Ayrıca, temas döngüleri gibi diğer alt-TAD yapıları, daha önce Capture-C63 ile tanımlanan Xist ve Ftx arasındaki döngü (Şekil 7A) ve yakın zamanda Capture Hi-C48 için benzer bir protokol kullanılarak tanımlanan Xist ve Xert arasındaki döngü (Şekil 7B) gibi Capture Hi-C verilerinden açıkça görülebilir. Linx, Chic1 ve Xite lokusları arasındaki Tsix-TAD içinde bilinen temas noktalarını oluşturanlar gibi Capture Hi-C profillerinin artan çözünürlüğü nedeniyle diğer kontaklar da daha hassas bir şekilde eşlenebilir (Şekil 7A).
Şekil 7'de gösterilen Hi-C verileriyle karşılaştırıldığında, Capture Hi-C çözünürlükte dört kat artışa izin verdi, ancak sıralama derinliğinin yalnızca dörtte birini gerektiriyordu (yani, 126 M okumalara karşı 571 M) (Ek Tablo 1). Çözünürlükteki bu artış, Şekil 6 ve Şekil 7'de gösterilen sıralama derinliğinde Hi-C tarafından tespit edilemeyen alt TAD'lerin ve döngü etkileşimlerinin algılanmasını sağlar. Capture Hi-C için açıklanan protokol, önceki yaklaşımlara kıyasla, geniş bir genomik bölgenin çok daha ayrıntılı, yüksek çözünürlüklü bir karakterizasyonuna izin verir.

Resim 1: Prob tasarımı. Prob tasarımı için kullanılan stratejinin şematik gösterimi. 3 Mb hedef bölge boyunca her bir DpnII kısıtlama bölgesinin 300 bp yukarı ve aşağı akış bölgeleri seçildi ve üst üste binen biyotinile RNA probları ile döşendi. Bu seçilen bölgelerden biri gösterilir, chrX: 102,474,805-102,475,500. Her probda 40'tan fazla baz tekrarlayan diziye izin verilmez. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 2: Hi-C ön sıralama kalite kontrollerini yakalayın . (A) 3C şablon kalite kontrollerinin temsili örneği. 200 ng DNA,% 1'lik bir agaroz jeline yüklendi. Şerit 1: 1 kb merdiven. Şerit 2: Sindirilmemiş, çapraz bağlı ve bozulmamış kromatin, >10 kb'de keskin bir bant olarak çalışır. Şerit 3: DpnII ile sindirilmiş çapraz bağlı kromatin, 1 kb ila 3 kb arasında bir yayma olarak çalışır. Şerit 4: Son 3C kütüphanesi veya şablonu; Sindirilmiş çapraz bağlı DNA fragmanlarının serbest uçları yeniden bağlanır. Düşük moleküler boyuttaki DNA yayması neredeyse tespit edilemez ve ligasyon ürünü >10 kb'lik bir bant olarak tespit edilir. (B) Yüksek hassasiyetli biyoanalizör DNA profillerinin temsili örnekleri. Sol üst: 150 bp ile 700 bp arasında parça boyutunun dağılımını gösteren başarıyla kesilmiş 3C kütüphanesi. Sağ üst: yetersiz kesilmiş 3C kütüphanesi. Kesilmemiş DNA, parçaların geniş zenginleşmesi olarak tespit edilir >2 kb. (C) Sol altta: SPRI boncukları kullanılarak 1:1 sol taraf boyutunda bir seçimi takiben kesilmiş DNA örneği. ~300 bp'lik parçalar zenginleştirilmiştir. Alt orta: Üreticinin protokolüne göre çift uçlu adaptörlerin bağlanmasından sonra PCR profilini önceden yakalayın. Sağ altta: çoğullanmış sıralama için adaptörler, sıralama ve dizin oluşturma primerleri dahil olmak üzere son Capture Hi-C kitaplığı. Kısaltmalar: bp = baz çiftleri, FU = keyfi floresan birimi. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 3: HiC-Pro ile Hi-C dizileme sonrası kalite kontrollerini yakalayın . (A) Dizileme çiftlerinin ilk eşi için referans genom üzerindeki haritalama hızı örneği. Açık mavi fraksiyon, HiC-Pro tarafından hizalanan ve bir ligasyon kavşağını kapsayan okumaları temsil eder. Bu metrik böylece deneysel ligasyon adımını doğrulamak için kullanılabilir. (B) Dizileme eşleri genom üzerinde hizalandıktan sonra, analiz için yalnızca benzersiz şekilde hizalanmış okuma çiftleri tutulur. (C) Sarkan uç, kendi kendine daire veya yeniden ligasyon gibi geçerli olmayan çiftler (kırmızı renkte) analizden atılır. Geçerli çiftlerin oranı, ligasyonun ve aşağı çekme verimliliğinin iyi bir göstergesidir. (D) Geçerli çiftler ayrıca intra/interkromozomal ve kısa/uzun menzilli temaslara ayrılabilir. PCR yapıtlarını temsil etmesi muhtemel yinelenen okuma çiftleri analizden atılır. (E) Alele özgü analiz için HiC-Pro, her ebeveyn genomu için bir veya iki eş tarafından desteklenen allelik okuma sayısını bildirir (yani, C57BL / 6J x CASTEi / J). Anne ve baba aleline atanan okumaların aynı kısmı beklenir. (F) Son olarak, temas haritalarını oluşturmak için yalnızca yakalama bölgesiyle çakışan geçerli çiftler seçilir. Yakalama-yakalama çiftleri hedeflenen bölge içindeki kişileri temsil ederken, yakalama-raporlayıcı çiftleri hedeflenen bölge ile hedef dışı bölge arasındaki etkileşimi içerir. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 4: Capture Hi-C protokolünün iş akışı. Farklı protokol adımlarının şematik gösterimi. Genom çapında 3C şablonunu oluşturmak için, kromatin önce formaldehit ile çapraz bağlanır ve daha sonra DpnII kısıtlama enzimi ile sindirilir. Serbest DNA uçları daha sonra yeniden bağlanır, çapraz bağlanma tersine çevrilir ve DNA saflaştırılır. Hedef bölgeyi kapsayan parçaları zenginleştirmek için, bir dizi biyotinillenmiş RNA probu 3C şablonuna hibridize edilir ve streptavidin aracılı aşağı çekme ile yakalanır. Yakalama kitaplıkları çoğullanmış sıralama için işlenir ve yüksek çözünürlüklü etkileşim haritaları olarak görselleştirilen hedef boyunca kromatin temaslarının sıklığını çıkarmak için geçerli ligasyon parçaları ölçülür. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 5: Fare X kromozomu üzerindeki Xic'i kapsayan bölgeye genel bakış. Fare X kromozomunun şematik gösterimi ve yakalanan 3 Mb bölgesinin yakınlaştırılması (ChrX: 102.475.000-105.475.000). Hedeflenen bölge, XCI'nın ana düzenleyici lokusu olan Xic'e karşılık gelen ~ 800 kb DNA içerir. Xic , uzun kodlamayan genleri, XCI'nın önemli bir oyuncusu olan Xist'i ve düzenleyici manzarasını içerir. Xist'in pozitif düzenleyicileri yeşil, negatif regülatörleri mor renkte gösterilmiştir. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 6: Yakalanan 3 Mb bölge boyunca Hi-C, 5C ve Hi-C etkileşim haritalarını yakalayın. (A) Fare Xic'ini kapsayan 3 Mb hedefin Hi-C etkileşim haritasını 10 kb çözünürlükte yakalayın (bu çalışma). (B) 6 kb çözünürlükte A'dakiyle aynı hedef bölgenin 5C etkileşim haritası (47'den yeniden işlenen veriler). Analizlerde yer almayan tekrarlayan bölgeler beyaz renkte maskelenir. 5C verileri kendi biyoinformatik işlemlerini gerektirir (bkz.47). Temizleme ve hizalamadan sonra, astar çözünürlüğündeki 5C haritaları, 6 kb'lik bir nihai çözünürlüğe ulaşmak için çalışan bir medyan (pencere = 30 kb, adım = 5) kullanılarak bağlanır. (C) 40 kb çözünürlükte A ve B'deki ile aynı genomik bölgenin Hi-C etkileşim haritası (64'ten yeniden işlenen veriler). Tüm etkileşim haritaları fare ESC'lerinden oluşturulmuştur. İzolasyon skoru cooltools kullanılarak hesaplanmıştır ve TAD sınırlarında yalıtım minimleri ile histogramlar olarak temsil edilir. TAD sınırları haritanın altında dikey çizgiler olarak gösterilir. Her çizginin yüksekliği sınır kuvvetini gösterir. Genler, transkripsiyon yönünü gösteren oklar olarak gösterilir. Yalnızca veya daha kesin olarak Capture Hi-C haritalarında algılanan alt TAD sınırları, sırasıyla Tsix ve Xist TAD'lardaki alt-TAD'lar için macenta ve mavi ok uçlarıyla gösterilir. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.

Şekil 7: Yakalanan bölge içinde 1 Mb boyunca Hi-C, 5C ve Hi-C etkileşim haritalarını yakalayın. (A) Fare Xic'ini kapsayan 1 Mb genomik bölgenin Hi-C etkileşim haritasını 5 kb çözünürlükte yakalayın (bu çalışma). (B) A'dakiyle aynı genomik bölgenin 5C etkileşim haritası. 6 KB çözünürlükte (47'den yeniden işlenen veriler). Analizlerde yer almayan tekrarlayan bölgeler beyaz renkte maskelenir. Not olarak, 5C verileri kendi biyoinformatik işlemlerini gerektirir (bkz.47). Temizleme ve hizalamadan sonra, astar çözünürlüğündeki 5C haritaları, 6 kb'lik bir nihai çözünürlüğe ulaşmak için çalışan bir medyan (pencere = 30 kb, adım = 5) kullanılarak bağlanır. (C) Hi-C'nin A ve B'sindeki ile aynı genomik bölgenin Hi-C etkileşim haritası 20 kb çözünürlükte (veriler 64'ten yeniden işlenir). Tüm etkileşim haritaları mESC'lerden oluşturulmuştur. İzolasyon skoru cooltools kullanılarak hesaplanmıştır ve TAD sınırlarında yalıtım minimleri ile histogramlar olarak temsil edilir. TAD sınırları haritanın altında dikey çizgiler olarak gösterilir. Her çizginin yüksekliği sınır kuvvetini gösterir. Genler, transkripsiyon yönüne işaret eden oklar olarak gösterilir. Capture Hi-C'de özel olarak veya daha hassas bir şekilde algılanan temas döngüleri, sırasıyla Tsix ve Xist TAD'lardaki döngüler için macenta ve mavi yıldız işaretleriyle gösterilir. Bu şeklin daha büyük bir versiyonunu görmek için lütfen buraya tıklayın.
Ek Tablo 1: Bu makalede kullanılan veri kümeleri için sıralama sonrası istatistikler: Capture Hi-C (bu çalışma), Hi-C64 ve 5C47. Bu Dosyayı indirmek için lütfen buraya tıklayın.