Optymalizacja zautomatyzowanych eksperymentów ChIP-seq dla ośmiu różnych markerów histonowych
W celu pomyślnego opracowania i walidacji zautomatyzowanych protokołów ChIP wybrano przeciwciała klasy ChIP-seq, które zostały wcześniej zwalidowane w manualnych eksperymentach ChIP-Seq (dane nie zostały przedstawione). Do niniejszego badania wybrano następujące przeciwciała klasy ChIP-seq: anty-H3K79me3, -H3K27me3, -H3K4me3, -H3K4me2, -H3K9ac, -H3K9/14ac, -H3K36me3 oraz -H3K9me3. Specyficzność wszystkich przeciwciał klasy ChIP-seq została wcześniej potwierdzona za pomocą techniki dot blot, macierzy peptydowych oraz eksperymentów Western Blot (dane nie zostały przedstawione). Przeprowadzono pilotażowe eksperymenty ChIP-qPCR z zastosowaniem zwiększających się ilości przeciwciał w celu określenia ich czułości (Rysunek 3). Przeanalizowano qPCR z co najmniej dwoma pozytywnymi i dwoma negatywnymi celami kontrolnymi; profile, w których wzbogacenie celu pozytywnego względem negatywnego było ponad pięciokrotne, zakwalifikowano do eksperymentów sekwencjonowania. Ważne jest przeprowadzanie eksperymentów ChIP i ChIP-seq z użyciem wysokiej jakości fragmentowanej chromatyny. Wszystkie eksperymenty ChIP przedstawione w niniejszej publikacji wykonano przy użyciu świeżej chromatyny. Możliwe jest również zamrożenie utrwalonych komórek w temperaturze -80 °C i przystąpienie do przygotowania oraz fragmentacji chromatyny w innym dniu. Jednakże chromatyna przygotowana z zamrożonych utrwalonych komórek może zachowywać się inaczej niż chromatyna przygotowana na świeżo, w związku z czym warunki sonikacji mogą wymagać optymalizacji dla każdej partii chromatyny. Podczas pracy z różnymi typami komórek można stosować bufory do fragmentacji o różnym składzie detergentów (SDS). Typy komórek, takie jak pierwotne linie komórkowe lub komórki hodowane w zawiesinie, są trudne do fragmentacji i wymagają wysokich stężeń SDS (1%), natomiast linie komórkowe łatwe do fragmentacji, takie jak HeLa, wymagają niskich stężeń SDS (0,1%) w buforach do fragmentacji.

Rysunek 3. Walidacja przeciwciał klasy ChIP przy użyciu systemu automatycznego. ChIP przeprowadzono z zastosowaniem króliczych przeciwciał poliklonalnych przeciwko H3K79me3, -H3K27me3, -H3K4me3, -H3K4me2, -H3K9ac, H3K9/14ac, -H3K36me3 oraz -H3K9me3 na fragmentowanej chromatynie z 1 miliona komórek HeLa-S3, w zależności od modyfikacji histonów. W urządzeniu automatycznym do eksperymentów z mianowania przeciwciał wykorzystano zautomatyzowane protokoły ChIP z objętościami roboczymi 200 µl. W każdym eksperymencie ChIP przetestowano ilości przeciwciał wynoszące 1, 2, 5 i 10 μg, a jako kontrolę negatywną zastosowano 2 μg IgG. Wzbogacenia oceniono za pomocą qPCR. Wyniki przedstawiono jako % inputu (względna ilość immunoprecypitowanego DNA w porównaniu do DNA wejściowego po analizie qPCR). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
Po walidacji i wyznaczeniu optymalnych ilości przeciwciał klasy ChIP do wykorzystania w systemie automatycznym, przeprowadzono zautomatyzowane eksperymenty ChIP-seq w celu wygenerowania profili sekwencjonowania dla każdej modyfikacji histonów (Rysunek 4).

Rysunek 4. Profile ChIP-seq histonów wygenerowane przez zautomatyzowane eksperymenty ChIP-seq. Rysunek przedstawia profile ChIP-seq w różnych regionach genomicznych dla H3K4me3, H3K9ac i H3K9/14acH3K4me2 H3K79me3 oraz H3K36me3. 4A pokazuje rozkład piku wzdłuż całego chromosomu X, a 4B rozkład w regionie 75 kb otaczającym gen GAPDH. 4C przedstawia profile H3K27me3, H3K36me3 i H3K4me3 w regionie 500 kb otaczającym gen MYT1, a 4D pokazuje rozkład H3K9me3 w regionie 200 kb otaczającym ZNF12. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Wygenerowano profile epigenetyczne dla sześciu różnych modyfikacji histonów związanych z ekspresją genów (H3K4me3, H3K9ac, H3K9/14ac, H3K4me2, H3K79me3 oraz H3K36me3). Rycina 4A przedstawia profile ChIP-seq wzdłuż chromosomu X dla różnych markerów histonowych. Wysoka korelacja pików zaobserwowana pomiędzy 6 różnymi profilami histonowymi wskazuje na zdolność systemu zautomatyzowanego do generowania dokładnych i wiarygodnych danych. Rycina 4B, 4C i 4D przedstawiają rozkład pików dla różnych modyfikacji histonowych w określonych regionach genomicznych.
Zautomatyzowane eksperymenty immunoprecypitacji chromatyny dla próbek obejmujących nawet 200 komórek
Minimalna liczba komórek, którą można wykorzystać w eksperymentach ChIP, zależy od jakości chromatyny, specyficzności i czułości przeciwciała oraz poziomu badanego białka lub modyfikacji histonów. Wybór wysokiej jakości przeciwciał klasy ChIP-seq jest kluczowy przy pracy z ograniczoną ilością próbki, a dobór optymalnych odczynników i różnych nośników zwiększa wydajność odzysku DNA i przyczynia się do sukcesu eksperymentu ChIP. Aby określić minimalną liczbę komórek, którą może przetworzyć zautomatyzowany protokół ChIP, w systemie zautomatyzowanym IP-Star przetestowano różne ilości chromatyny, przeciwciał oraz kulek magnetycznych, stosując odczynniki do ChIP specjalnie zoptymalizowane pod kątem pracy z niskimi ilościami chromatyny.
W pierwszej kolejności chromatynę z 10 000 komórek poddano sonikacji zgodnie z opisem w protokole. Wyniki ChIP potwierdzono za pomocą qPCR (Rysunek 5A), który wykazuje znaczące wzbogacenie z użyciem przeciwciał H3K4me3 w regionach kontroli pozytywnej oraz pomijalny sygnał w regionach kontroli negatywnej. Dla porównania i w celu potwierdzenia spójności przedstawiono dodatkowe dane uzyskane z zastosowaniem przeciwciał H3K27ac, H3K9me3 oraz H3K27me3 przy użyciu 10 000 komórek.
Następnie przeprowadzono zautomatyzowane eksperymenty ChIP, aby zademonstrować zdolność zautomatyzowanego systemu do pracy z niewielkimi ilościami komórek przy użyciu tego samego przeciwciała przeciwko H3K4me3. Zautomatyzowany ChIP działał sprawnie, co wykazała seria dziesięciu reakcji IP, które były powtarzalne i wysoce porównywalne z wynikami ręcznego ChIP (Rycina 5B). Przeprowadzono eksperymenty ręczne i zautomatyzowane, a korzyści z zastosowania protokołów zautomatyzowanych przejawiły się w zmniejszeniu zmienności między poszczególnymi eksperymentami (Rycina 5C).

Rycina 5. Optymalizacja eksperymentów ChIP i Auto ChIP na 10 000 komórek Manualne eksperymenty ChIP przeprowadzono na 10 000 komórek, stosując 0,25 μg przeciwciał H3K4me3, 0,1 μg H3K27ac, 0,5 μg H3K9me3 oraz 0,25 μg H3K27me3. Jako kontrolę zastosowano identyczne ilości IgG króliczego. qPCR przeprowadzono z użyciem starterów dla dwóch loci pozytywnych i dwóch loci negatywnych dla każdego analizowanego ChIP. Rycina 5A przedstawia odzysk, wyrażony jako procent DNA wejściowego (relatywna ilość DNA immunoprecypitowanego w porównaniu do DNA wejściowego po analizie qPCR). Rycina 5B przedstawia 10 reakcji ChIP przeprowadzonych na systemie IP-Star Compact z użyciem 0,25 μg poliklonalnego przeciwciała H3K4me3 oraz 0,25 μg IgG króliczego jako negatywnego przeciwciała kontrolnego. Następnie wykonano analizę qPCR z użyciem starterów dla loci pozytywnych (promotor EIF4A2 i TSS GAPDH) oraz loci negatywnych (ekson 2 Myoglobiny i Sat2). Rycina przedstawia odzysk, wyrażony jako procent DNA wejściowego (relatywna ilość DNA immunoprecypitowanego w porównaniu do DNA wejściowego po analizie qPCR). Rycina 5C przedstawia dane ChIP dla H3K4me3 z 10 manualnych eksperymentów ChIP w porównaniu do 10 zautomatyzowanych eksperymentów ChIP. Słupki błędów reprezentują odchylenia standardowe dla każdego z dziesięciu powtórzeń. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
W celu oceny czułości zautomatyzowanych protokołów ChIP przeprowadzono eksperymenty z wykorzystaniem liczby komórek w zakresie od 100 000 do 200 komórek na jedną reakcję IP. Zastosowano przeciwciało anty-H3K27me3, ponieważ jest to bardzo powszechna modyfikacja histonowa. Użycie innych przeciwciał antyhistonowych lub nieantyhistonowych może wymagać większej lub mniejszej liczby komórek, w zależności od obfitości epitopu oraz jakości przeciwciała. Eksperymenty zwalidowano za pomocą ilościowego PCR, w wyniku czego zaobserwowano, że zmniejszenie ilości kulek oraz przeciwciał redukuje tło w badaniach, co pozwala na uzyskanie pomyślnych wyników ChIP-qPCR przy użyciu nawet 200 komórek z przeciwciałem (Rysunek 6).

Rysunek 6. Zautomatyzowane analizy ChIP na 200 komórkach. Komórki HeLa-S3 i przeciwciało skierowane przeciwko H3K27me3. Chromatynę poddano fragmentacji z 1 miliona komórek, a następnie do każdej reakcji ChIP wykorzystano rozcieńczenia szeregowe tej chromatyny (od równowartości 100 000 do 200 komórek). W doświadczeniu na 100 000 komórek zastosowano 1 µg H3K27me3 oraz 10 μl magnetycznych kulek pokrytych białkiem A; w przypadku 10 000 i 1 000 komórek zastosowano 0,5 µg H3K27me3 oraz 10 µl kulek; natomiast dla 500 i 200 komórek użyto 0,25 µg H3K27me3 oraz 5 µl kulek. Jako przeciwciało kontrolne negatywne zastosowano odpowiednio 1 µg i 0,5 µg króliczego IgG podczas przeprowadzania doświadczeń z 100 000 i 1 000 komórek. 6A przedstawia stopień zajętości genów TSH2B i GAPDH w % względem inputu. 6B przedstawia względną zajętość TSH2B w stosunku do negatywnego kontrolnego regionu genomicznego GAPDH.
Analiza końcowa wyników ChIP-seq dla 10 000 komórek
W celu oceny ogólnej jakości zautomatyzowanych eksperymentów ChIP-seq przy niskiej początkowej liczbie komórek, przeprowadzono zautomatyzowane analizy ChIP-seq z zastosowaniem 0,25 µg przeciwciała H3K4me3 na 10 000 komórek HeLa, natomiast eksperymenty ChIP na 100 000 komórek HeLa-S3 posłużyły jako kontrola pozytywna eksperymentu. Zautomatyzowane biblioteki przygotowano przy użyciu odczynników z zestawu do przygotowania bibliotek MicroPlex, dostosowanych do bibliotek o niskiej zawartości DNA. Należy zauważyć, że mimo możliwości przeprowadzenia pomyślnych zautomatyzowanych eksperymentów ChIP z liczbą komórek mniejszą niż 10 000, ilość wyizolowanego DNA nie będzie wystarczająca do przygotowania bibliotek przy użyciu odczynników z zestawu. Generowanie klastrów oraz sekwencjonowanie przeprowadzono zgodnie z instrukcjami producenta. Analizy bioinformatyczne po sekwencjonowaniu wykazały wyjątkowe wyniki dla próbek ChIP o niskiej liczbie komórek. Zbiór danych 30 pg (odpowiadający 10 000 komórek materiału wyjściowego) charakteryzuje się niskim szumem tła i wysoce wiarygodnymi pikami wzbogacenia, które zostały potwierdzone zarówno przez zbiór danych 300 pg (odpowiadający 100 000 komórek materiału wyjściowego), jak i zbiorem danych H3K4me3 wygenerowanym przez Broad Institute dla projektu ENCODE, który posłużył jako referencja zewnętrzna. Istotne jest zwrócenie uwagi na dane dotyczące współczynnika nakładania się Top 40, który odnosi się do standardowej metody stosowanej w projekcie ENCODE11; w tej metodzie ChIP-seq uznaje się za powtarzalny, jeśli przy porównaniu dwóch zbiorów danych występuje co najmniej 80% nakładania się najlepszych 40% pików uszeregowanych według wyniku istotności. Zbiór danych 30 pg spełnia te kryteria w porównaniu zarówno ze zbiorem danych 300 pg (uwzględniając wszystkie jego piki, a nie tylko najlepsze 40%), jak i z danymi Broad Institute (Tabela 1). Zbiór danych 300 pg wykazuje niemal identyczne piki jak dane Broad Institute, przy współczynniku nakładania się Top 40 wynoszącym 98% (Rysunek 7).

Rycina 7. Analizy ChIP i generowanie bibliotek z 10 000 komórek Eksperymenty ChIP-seq przeprowadzono na 10 000 i 100 000 komórkach HeLa-S3 z wykorzystaniem przeciwciała H3K4me3 (0,25 µg/µl). Tagi o długości 35 bp zmapowano do genomu ludzkiego za pomocą alignera ELAND. Podczas następczego wyznaczania pików (peak calling) program SICER mógł wiarygodnie zidentyfikować wzbogacenia zarówno z małej liczby komórek, jak i z milionów komórek. Zestawy danych zostały przeanalizowane i porównane ze sobą oraz z danymi referencyjnymi wygenerowanymi przez Broad Institute. Próbki z niską liczbą komórek są spójne i wykazują bardzo wysokie podobieństwo. Próbka 30 pg spełnia kryteria ENCODE11 (min. 80% z górnych 40% pików powinno nakładać się na siebie). Aby zobaczyć powiększoną wersję tej ryciny, kliknij tutaj.

Tabela 1.