Określenie kluczowych czynników sterujących na podstawie transkrypomicznych danych z lotów kosmicznych pomoże NASA w określaniu ryzyk zdrowotnych i opracowywaniu potencjalnych środków przeciwdziałania negatywnym skutkom dla zdrowia astronautów. W naszej niedawnej publikacji zastosowaliśmy powyższe kroki i wykorzystaliśmy zestawy danych GeneLab, aby skutecznie wykazać nowe odkrycie, zgodnie z którym stężenia CO2 na ISS mogą wpływać na zdrowie36. Wykorzystaliśmy również powyższą technikę w innych badaniach, aby pomyślnie określić kluczowe czynniki sterujące badanym układem45,46,47,48,49,50. Poniżej pokażemy, w jaki sposób wyniki uzyskane przy użyciu tego protokołu mogą zostać skutecznie wykorzystane do określenia kluczowych czynników sterujących.
W niniejszym badaniu skupiliśmy się przede wszystkim na różnicach biologicznych występujących u gryzoni utrzymywanych w kontrolnych warunkach naziemnych w klatkach Rodent Habits oraz w kontrolnych warunkach wiwariowych. Jak opisano powyżej, jest to klucz do lepszego zrozumienia tych dwóch siedlisk, co dostarczy nam informacji o możliwych czynnikach zakłócających, które mogą wpływać na stan zdrowia ze względu na środowisko panujące na ISS. W przypadku wszystkich eksperymentów z gryzoniami podczas lotów kosmicznych, kontrole naziemne są również niezbędne do określenia, które czynniki biologiczne są bezpośrednio związane z lotem w kosmos, a które wynikają z warunków środowiskowych na ISS. Zgodnie z protokołem, warunki środowiskowe w siedlisku wiwariowym nie obejmują ekspozycji na wyższy poziom CO2, który występuje w siedlisku Rodent Habitat. W siedlisku wiwariowym panuje normalny poziom CO2 występujący na Ziemi (obecnie od 300 do 380 ppm). Temperatura i wilgotność w obu siedliskach są zbliżone.
Aby określić kluczowe geny odpowiedzialne za różnice między gryzoniami utrzymywanymi w kontrolnych warunkach naziemnych Rodent Habitat a kontrolnymi warunkami naziemnymi w wiwarium, wykorzystaliśmy następujące zestawy danych z platformy GeneLab: GLDS-21, GLDS-111, GLDS-25 oraz GLDS-63. Analizę w celu wyłonienia istotnych genów przeprowadzono zgodnie z powyższym opisem pomiędzy grupą Rodent Habitat (wcześniej AEM) a kontrolą w wiwarium, niezależnie dla każdego zestawu danych. Wykresy PCA wykazały grupowanie powtórzeń biologicznych (Rycina 4 przedstawia wykresy PCA dla GLDS-21). Na podstawie wstępnie przetworzonych danych określiliśmy geny wiodące (leading-edge genes) z różnych zestawów genów GSEA. Wykorzystując geny o zmianie krotności 1,2-fold-change (log2), byliśmy w stanie przewidzieć geny zaangażowane w regulację nadrzędną, szlaki kanoniczne i biofunkcje. Następnie dla każdego zestawu danych znaleźliśmy wspólne/nakładające się geny zaangażowane we wszystkich analizowanych przypadkach (Rycina 5). Obecnie uważa się, że geny te odpowiadają za reakcję gryzoni w Rodent Habitats (lub AEM) w porównaniu z kontrolą w wiwarium. Reprezentacja sieciowa powiązań między tymi kluczowymi genami wskazuje centralne węzły dla każdego analizowanego zestawu danych (Rycina 6). Na przykład MAPK1 jest centralnym węzłem dla tkanek mięśni szkieletowych myszy z misji STS-108 (Rycina 6A). Interpretuje się to jako gen sterujący kluczowymi genami, który najprawdopodobniej odgrywa centralną rolę w powodowaniu różnic biologicznych u myszy utrzymywanych w Rodent Habitats w porównaniu do klatek w wiwarium. W naszej poprzedniej pracy omawiamy na podstawie istniejącej literatury naukowej, w jaki sposób kluczowe geny te są powiązane z odpowiedzią na CO2 i w jaki sposób mogą one odpowiadać za obserwowane u myszy zmiany biologiczne36.
Stosując podejście biologii systemowej, wyznaczyliśmy następnie „główny regulator” (master regulator), który łączy wszystkie zestawy danych/tkanki i jest potencjalnie odpowiedzialny za uniwersalne efekty biologiczne u gryzoni utrzymywanych w AEM w porównaniu do klatek w wiwarium. Dokonano tego poprzez wskazanie genu spośród wszystkich analizowanych zestawów danych, który wykazuje największą liczbę połączeń podczas konstruowania sieci ze wszystkich kluczowych genów. Wykazaliśmy, że MAPK1 jest genem o największej liczbie połączeń i stanowi centralny węzeł spośród wszystkich kluczowych genów (Rysunek 7). Aby potwierdzić, czy MAPK1 może być odpowiedzialny za zmiany biologiczne u myszy wynikające z wyższych poziomów CO2 w AEM, przeszukaliśmy literaturę naukową w poszukiwaniu dowodów wspierających tę hipotezę. Znaleźliśmy kilka badań wskazujących na korelację MAPK1 z CO259 oraz hipoksją19,60,61.

Rycina 1: Porównanie siedliska dla gryzoni (wcześniej AEM) z klatkami wiwaryjnymi. (A) Zdjęcie klatki AEM dostarczonej przez NASA (Źródło: NASA/Dominic Hart). (B) Standardowa klatka wiwaryjna obecnie stosowana (zdjęcie wykonane w naszym laboratorium). Rycina została zmodyfikowana na podstawie pracy Beheshti et al.36. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: System sprzętowy Rodent Habitat z trzema różnymi modułami wykorzystywanymi podczas transportu na misje kosmiczne i z powrotem. Moduł po lewej (A) to moduł Rodent Habitat (wcześniej AEM), moduł środkowy (B) to Transporter, a moduł po prawej (C) to Animal Access Unit (AAU). (D)Mouse Transfer Box (MTB). (Źródło: NASA/Dominic Hart). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 3: Przykładowy schemat analizy, który można wykorzystać w interfejsie GeneLab Galaxy do przetwarzania danych RNA-seq. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 4: Analiza głównych składowych (PCA) reprezentatywnego zbioru danych po etapach wstępnego przetwarzania. Przedstawiono zbiór danych GLDS-21 dla AEM vs. klatki w vivarium dla mięśni szkieletowych myszy z misji STS-118. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Diagram Venna przedstawiający kluczowe geny zidentyfikowane przy użyciu różnych narzędzi do przewidywania szlaków. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 6: Kluczowe geny określone dla wszystkich warunków i tkanek mysich w porównaniu klatek AEM vs. klatek w wiwariach. (A-E) Reprezentacja sieciowa kluczowych genów dla każdego zestawu danych/tkanki gryzonia. Log2 fold-change (przy progu 1,2-krotnej zmiany) ekspresji genów posłużył do uzyskania różnych odcieni zieleni dla fold-change w genach downregulowanych, podczas gdy różne odcienie czerwieni obrazują fold-change w genach upregulowanych. Im ciemniejszy odcień zieleni lub czerwieni, tym większa wartość fold-change. Rycina została zmodyfikowana na podstawie pracy Beheshti et al.36. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 7: Wyznaczanie „głównego regulatora” dla gryzoni w klatkach Rodent Habitat w porównaniu do klatek wiwaryjnych. Połączenia między wszystkimi poszczególnymi kluczowymi genami (Rycina 6) zostały określone i przedstawione jako sieć za pomocą programu IPA. Sieć została przedstawiona na wykresie radialnym, z najbardziej połączonym kluczowym genem, MAPK1, w centrum. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Rycina dodatkowa 1: Integracja GeneLab-GenomeSpace z ISACreator w celu usprawnienia operacji przetwarzania danych. Kliknij tutaj, aby pobrać tę rycinę.
Rycina dodatkowa 2: Zrzut ekranu z wyszukiwań w GeneLab z wykorzystaniem federacji/integracji z zewnętrznymi heterogenicznymi bazami bioinformatycznymi (GEO, PRIDE, MG-RAST). Kliknij tutaj, aby pobrać tę rycinę
Rycina uzupełniająca 3: Zrzut ekranu wspólnej przestrzeni roboczej GeneLab przedstawiający zarządzanie kontami użytkowników oraz kontrole dostępu (np., foldery prywatne, współdzielone i publiczne). Kliknij tutaj, aby pobrać tę rycinę