$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
W ostatniej dekadzie, metabolomika stała się nauką omiczną dzięki postępom w technologiach analitycznych, takich jak chromatografia gazowa-spektrometria mas (GC-MS) i chromatografia cieczowa-spektrometria mas (LC-MS). Techniki te umożliwiają jednoczesny pomiar setek do tysięcy metabolitów małocząsteczkowych, tworząc złożone, wielowymiarowe zestawy danych. Eksperymenty metabolomiczne mogą być przeprowadzane w trybie ukierunkowanym lub nieukierunkowanym. Ukierunkowane eksperymenty metabolomiczne służą do pomiaru określonych klas metabolitów. Są one zwykle oparte na hipotezach, podczas gdy podejścia nieukierunkowane próbują zmierzyć jak najwięcej metabolitów i mają charakter generujący hipotezy. Testy celowane zwykle obejmują wzorce wewnętrzne, a tym samym pozwalają na bezwzględną kwantyfikację metabolitów będących przedmiotem zainteresowania. W przeciwieństwie do tego, testy nieukierunkowane umożliwiają względną kwantyfikację i obejmują wiele nieznanych metabolitów1.
Analiza danych metabolomicznych to wieloetapowy proces, który wykorzystuje wiele specjalistycznych narzędzi programowych1. Można go podzielić na następujące trzy główne etapy: (1) przetwarzanie danych i kontrola jakości, (2) analiza statystyczna oraz (3) interpretacja danych biologicznych. Opisane tutaj narzędzia zostały zaprojektowane tak, aby umożliwić ostatni etap analizy.
Intuicyjnym i popularnym sposobem interpretacji danych metabolomicznych jest mapowanie pomiarów eksperymentalnych na szlaki metaboliczne. Aby to osiągnąć, zaprojektowano wiele narzędzi2,3,4,5, w tym Metscape, opracowany przez naszą grupę6. Mapowanie ścieżek jest często łączone z analizą wzbogacania, która pomaga zidentyfikować najważniejsze ścieżki7,8. Techniki te po raz pierwszy zyskały na znaczeniu w analizie danych dotyczących ekspresji genów i zostały z powodzeniem zastosowane do analizy danych proteomicznych i epigenomicznych9,10,11,12,13. Analiza danych metabolomicznych wiąże się jednak z szeregiem wyzwań dla podejść opartych na wiedzy. Po pierwsze, oprócz metabolitów endogennych, testy metabolomiczne mierzą związki egzogenne, w tym te, które pochodzą z odżywiania i innych źródeł środowiskowych. Związki te, podobnie jak metabolity wytwarzane przez bakterie, nie mogą być mapowane na szlaki ludzkie lub metaboliczne innych organizmów eukariotycznych. Co więcej, pokrycie szlaku metabolizmu wtórnego i metabolizmu lipidów nie pozwala obecnie na mapowanie w wysokiej rozdzielczości na poziomie, który z łatwością wspierałby biologiczną interpretację danych14,15.
Techniki analizy sieci opartej na danych mogą pomóc w przezwyciężeniu tych wyzwań. Na przykład sieci oparte na korelacjach mogą pomóc w określeniu relacji między znanymi i nieznanymi metabolitami oraz ułatwić dodawanie adnotacji do niewiadomych16. Chociaż obliczanie współczynników korelacji Pearsona jest najprostszym podejściem do ustalania liniowych zależności między metabolitami, wadą jest to, że obejmuje zarówno bezpośrednie, jak i pośrednie powiązania17,18,19. Alternatywą jest obliczenie współczynników korelacji cząstkowej, które umożliwiają rozróżnienie między skojarzeniami bezpośrednimi i pośrednimi. Graficzne modelowanie gaussowskie (GGM) może być używane do szacowania sieci korelacji cząstkowych. GGM wymaga jednak, aby wielkość próby i liczba cech były porównywalne. Warunek ten jest rzadko spełniony w nieukierunkowanych danych LC-MS, które zawierają pomiary dla tysięcy cech metabolicznych. Aby przezwyciężyć to ograniczenie, można zastosować techniki regularyzacji. Graficzne lasso (Glasso) i regresja węzłowa są popularnymi metodami regularnego szacowania sieci korelacji cząstkowej16,20.
Pierwsze z prezentowanych tutaj narzędzi bioinformatycznych, CorrelationCalculator16, opiera się na algorytmie zdetensyjnej rzadkiej korelacji cząstkowej (DSPC). DSPC opiera się na oszczędnym graficznym modelowaniu lasso. Podstawowym założeniem algorytmu jest to, że liczba połączeń między metabolitami jest znacznie mniejsza niż liczba próbek, tj. sieć częściowej korelacji metabolitów jest rzadka. To założenie pozwala DSPC odkryć łączność między dużą liczbą metabolitów przy użyciu mniejszej liczby próbek, wykorzystując regularne techniki regresji. Ponadto, używając kroku odchylania dla unormowanych oszacowań regresji, uzyskuje rozkłady próbkowania dla parametrów krawędzi, które można wykorzystać do konstruowania przedziałów ufności i testowania hipotez będących przedmiotem zainteresowania (np. obecność/brak pojedynczej krawędzi lub grupy krawędzi). Obecność lub brak krawędzi w sieci korelacji cząstkowej można zatem formalnie przetestować przy użyciu obliczonych wartości p.
Kalkulator korelacji okazał się bardzo przydatny do analizy pojedynczej grupy16; jednak celem wielu eksperymentów metabolomicznych jest analiza różnicowa dwóch lub więcej warunków. Chociaż CorrelationCalculator może być stosowany dla każdej z grup osobno w celu wygenerowania sieci korelacji cząstkowych dla każdego warunku, takie podejście ogranicza liczbę próbek, które można wykorzystać do generowania sieci. Ponieważ wystarczająco duża wielkość próby jest jednym z najważniejszych czynników branych pod uwagę w analizie opartej na danych, metody, które mogą wykorzystać wszystkie dostępne próbki w danych do budowy sieci, są wysoce pożądane. Podejście to zostało zaimplementowane w drugim przedstawionym tutaj narzędziu, o nazwie Filigree21. Filigree opiera się na wcześniej opublikowanym algorytmie Differential Network Enrichment Analysis (DNEA), 22. Tabela 1 przedstawia aplikacje i przepływ pracy obu narzędzi.
| Liczba warunków doświadczalnych (k) | k = 1 | k = 2 |
| Narzędzie programowe | Kalkulator korelacji | filigran |
| Dane wejściowe | • Macierz danych metabolity x próbki | • Macierz danych metabolity x próbki
• Grupy eksperymentalne |
Przepływ pracy
• Obróbka wstępna
• Szacowanie sieci
• Klastrowanie sieci
• Analiza wzbogacenia |
• Transformacja logarytmu; autoskalowanie
• DSPC
• Za pośrednictwem zewnętrznych aplikacji
•Nie |
• Transformacja logarytmu; autoskalowanie
• Wspólne szacowanie sieci
• Grupowanie konsensusu
• Sieć NetGSA |
| Wizualizacja danych | Za pomocą zewnętrznej aplikacji, np. Cytoscape | Za pomocą zewnętrznej aplikacji, np. Cytoscape |
| Testowanie modułów metabolicznych pod kątem związku z wynikiem będącym przedmiotem zainteresowania (opcjonalnie) | Za pomocą zewnętrznych aplikacji | Za pomocą zewnętrznych aplikacji |
Tabela 1: Zakres zastosowania i przepływ pracy CorrelationCalculator i Filigran.