Estrogeny są ważnymi regulatorami wielu procesów fizjologicznych zarówno u kobiet, jak i mężczyzn, w tym tkanek rozrodczych, tkanek metabolicznych, mózgu i kości1. Oprócz korzystnego działania w tych tkankach, estrogeny napędzają również nowotwory, które powstają w tkankach sutkowych i rozrodczych. Estrogeny działają głównie poprzez ER, aby wywołać efekty specyficzne dla typu komórki. Głębokie sekwencjonowanie transkryptów regulowanych przez ERα za pomocą RNA-Seq i analiza miejsc wiązania ERα DNA w całym genomie za pomocą ChIP-Seq okazało się przydatne do zrozumienia, jak ERα działa w różnych tkankach i nowotworach, które z nich wynikają. My i inni opublikowaliśmy profile ekspresji genów związane z różnymi receptorami (ERα vs. ERβ)2,3, różnymi ligandami3-5 i różnymi koregulatorami2,4,6,7.
RNA-Seq to główna metoda badania transkryptomu, oferująca wyższą precyzję i wydajność w porównaniu do analizy ekspresji genów opartej na mikromacierzyach8. RNA uzyskane z linii komórkowych2-4,7, tkanki lub próbki guza są sekwencjonowane, mapowane do dostępnych zespołów genomowych i identyfikowane są geny regulowane różnicowo. Immunoprecypitacja chromatyny (ChIP) jest stosowana do analizy czynnika transkrypcyjnego i koregulatora wiązania chromatyny ze znanymi miejscami wiązania regulatorowego. ChIP-Seq (ChIP, po którym następuje sekwencjonowanie o wysokiej przepustowości) zapewnia bezstronne wykrywanie globalnych miejsc wiązania. Metabolomika jest kolejnym coraz częściej stosowanym podejściem biologii systemowej, które mierzy ilościowo, dynamiczną wieloparametryczną reakcję żywych systemów na różne bodźce, w tym chemikalia i zaburzenia genetyczne.
Wykonując globalne profilowanie metaboliczne, można uzyskać funkcjonalny odczyt z komórek, tkanek i krwi. Ponadto informacje z eksperymentów transkryptomowych nie zawsze odzwierciedlają rzeczywiste zmiany poziomu enzymów, które przyczyniają się do szlaków biochemicznych. Połączona analiza danych dotyczących transkryptomu i metabolomu pozwala nam zidentyfikować i skorelować zmiany w ekspresji genów z rzeczywistymi zmianami metabolitów. Wykorzystanie informacji ze wszystkich tych wielkoskalowych zbiorów danych dostarcza szczegółów mechanistycznych, aby zrozumieć rolę czynników transkrypcyjnych regulujących złożone systemy biologiczne, zwłaszcza te, które odnoszą się do rozwoju człowieka i chorób takich jak rak i cukrzyca.
Złożona natura genomu ssaków sprawia, że trudno jest zintegrować i w pełni zinterpretować dane uzyskane z eksperymentów transkryptomu, cistromu i metabolomu. Identyfikacja funkcjonalnych zdarzeń wiązania, które prowadziłyby do zmian w ekspresji genów docelowych, jest ważna, ponieważ po zidentyfikowaniu funkcjonalnych miejsc wiązania; Późniejsza analiza, w tym analiza motywu wiązania transkrypcji (TF), może być przeprowadzana z większą dokładnością. Prowadzi to do identyfikacji biologicznie znaczących kaskad i mechanizmów TF. Również bezpośrednie porównanie eksperymentów RNA-Seq i ChIP-Seq nie zawsze jest możliwe, ponieważ dane z każdego eksperymentu mają różne skale i szumy, a w niektórych przypadkach znaczące sygnały są przesłonięte przez szum populacyjny. Nie znamy żadnego badania, w którym zintegrowano informacje z tych trzech niezależnych, ale powiązanych podejść, aby zrozumieć bezpośrednią regulację metaboliczną przez ERα w raku piersi. Dlatego naszym ogólnym celem w tym artykule jest powiązanie produktywnych zdarzeń wiązania z ekspresją genów i zmianami metabolitów. Aby osiągnąć ten cel, połączyliśmy dane z RNA-Seq, ChIP-Seq i eksperymentów metabolomicznych oraz zidentyfikowaliśmy te indukowane estrogenem zdarzenia wiązania ERα, które doprowadziłyby do zmian ekspresji genów w szlakach metabolicznych. Po raz pierwszy udostępniamy kompletny zestaw protokołów (rysunek 1) do generowania ChIP-Seq, RNA-Seq i profilowania metabolomicznego oraz przeprowadzania integracyjnej analizy danych w celu odkrycia nowych obwodów genowych regulujących metabolizm linii komórkowych raka piersi.