$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Określenie kluczowych czynników na podstawie danych transkryptomicznych lotów kosmicznych pomoże NASA w określeniu zagrożeń dla zdrowia i opracowaniu potencjalnych środków zaradczych w celu zwalczania negatywnego wpływu na zdrowie astronautów. W naszej ostatniej publikacji wykonaliśmy powyższe kroki i wykorzystaliśmy zestawy danych GeneLab, aby z powodzeniem pokazać nowe odkrycie, że stężenia CO2 na ISS mogą wpływać na zdrowie36. Wykorzystaliśmy również powyższą technikę w innych badaniach, aby z powodzeniem określić kluczowe czynniki napędzające badany system45,46,47,48,49,50. Tutaj pokażemy, w jaki sposób wyniki z użycia tego protokołu można z powodzeniem wykorzystać do określenia kluczowych czynników.
W tym badaniu skupiliśmy się przede wszystkim na różnicach biologicznych, które występują u gryzoni trzymanych w kontrolach naziemnych dla gryzoni i w kontrolach wiwariów. Jak opisano powyżej, jest to klucz do lepszego zrozumienia tych dwóch siedlisk, co dostarczy nam informacji na temat możliwych czynników zakłócających, które mogą wpływać na zdrowie ze względu na środowisko na ISS. W przypadku wszystkich eksperymentów z lotami kosmicznymi gryzoni, te naziemne kontrole są również niezbędne do określenia, które czynniki biologiczne są bezpośrednio związane z lotami kosmicznymi lub ze względu na warunki środowiskowe na ISS. Jak stwierdzono w protokole, warunki środowiskowe siedliska wiwarium nie są narażone na wyższy poziom CO2 , który występuje w siedlisku gryzoni. Siedlisko wiwarium ma normalny poziom CO2 obecny na Ziemi (obecnie wynosi od 300 do 380 ppm). Temperatura i wilgotność w obu siedliskach są podobne.
Użyliśmy następujących zestawów danych z platformy GeneLab, aby określić kluczowe geny między gryzoniami umieszczonymi w naziemnych kontrolach siedlisk gryzoni a kontrolami naziemnymi wiwariów, które są odpowiedzialne za napędzanie różnic między dwoma siedliskami: GLDS-21, GLDS-111, GLDS-25 i GLDS-63. Analizę w celu określenia istotnych genów przeprowadzono zgodnie z powyższym opisem między siedliskiem gryzoni (dawniej AEM) a grupą kontrolną wiwarium niezależnie dla każdego zestawu danych. Wykresy PCA pokazały grupowanie replik biologicznych (Rysunek 4 pokazuje wykresy PCA dla GLDS-21). Na podstawie wstępnie przetworzonych danych określiliśmy wiodące geny z różnych zestawów genów GSEA. Korzystając z genów z 1,2-krotną zmianą (log2), byliśmy w stanie przewidzieć geny zaangażowane w przewidywania dotyczące regulatorów wyższego szczebla, szlaków kanonicznych i biofunkcji. Dla każdego zestawu danych znaleźliśmy następnie wspólne/nakładające się geny zaangażowane we wszystkie geny (Rysunek 5). Obecnie uważa się, że geny te napędzają reakcję między gryzoniami w siedliskach gryzoni (lub AEM) a kontrolami wiwariów. Sieciowa reprezentacja sposobu, w jaki te kluczowe geny się łączą, pokazuje centralne węzły dla każdego analizowanego zestawu danych (Rysunek 6). Na przykład MAPK1 jest centralnym węzłem dla tkanek mięśni szkieletowych STS-108 od myszy (Figura 6A). Byłoby to interpretowane jako gen, który napędza kluczowe geny i najprawdopodobniej jest głównym graczem powodującym różnice biologiczne u myszy trzymanych w siedliskach gryzoni w porównaniu z klatkami wiwarium. W naszej poprzedniej pracy omówiliśmy, w jaki sposób te kluczowe geny są powiązane z odpowiedzią na CO2 na podstawie istniejącej literatury naukowej i jak te geny mogą być odpowiedzialne za zmiany biologiczne obserwowane u myszy36.
Przyjmując podejście do biologii systemowej, następnie określiliśmy "główny regulator", który łączy wszystkie zestawy danych/tkanki i jest potencjalnie odpowiedzialny za uniwersalne efekty biologiczne u gryzoni trzymanych w AEM w porównaniu do klatek wiwariowych. Dokonano tego poprzez określenie genu ze wszystkich zestawów danych, który jest najbardziej połączony podczas konstruowania sieci ze wszystkich kluczowych genów. Udało nam się wykazać, że MAPK1 jest najbardziej połączonym genem i centralnym węzłem ze 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, przejrzeliśmy literaturę naukową w poszukiwaniu dowodów potwierdzających. Znaleźliśmy kilka badań wskazujących na korelację MAPK1 z CO259 i hypoxia19,60,61.

Rysunek 1: Siedlisko gryzoni (poprzednio AEM) w porównaniu do klatek wiwariowych.(A) Zdjęcie klatki AEM dostarczone przez NASA (Credits: NASA/Dominic Hart). (B) Standardowa klatka wiwarium, która jest obecnie używana (zdjęcie wykonane przez nasze laboratorium). Ten rysunek został zmodyfikowany z Beheshti et al.36. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: System sprzętowy dla gryzoni z trzema różnymi modułami zaangażowanymi podczas transportu do i z misji kosmicznych. Lewy moduł (A) to moduł siedliska gryzoni (dawniej AEM), moduł środkowy (B) to transporter, a prawy moduł (C) to jednostka dostępu do zwierząt (AAU). orazSkrzynka transferowa myszy (MTB). (Kredyty: NASA/Dominic Hart). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Przykładowy proces analizy, który może być wykorzystany w interfejsie GeneLab Galaxy do przetwarzania danych sekwencyjnych RNA. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Analiza głównych składowych (PCA) reprezentatywnego zbioru danych po wstępnym przetwarzaniu. Zestaw danych GLDS-21 dla AEM i klatki wiwarium jest pokazany dla mysich mięśni szkieletowych z misji STS-118. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Diagram Venna przedstawiający, jakie kluczowe geny są określane za pomocą różnych narzędzi do przewidywania ścieżek. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Kluczowe geny określone dla wszystkich schorzeń i tkanek myszy między AEM a AEM. klatki wiwariowe. (A-E) Sieciowa reprezentacja kluczowych genów dla każdego zestawu danych/tkanki gryzonia. Logarytmiczne zmiany krotności (z odcięciem 1,2-krotnej zmiany) ekspresji genów wykorzystano do uzyskania różnych odcieni zieleni dla zmiany fałdowania w genach regulowanych w dół, podczas gdy różne odcienie czerwieni obrazują zmianę fałdowania w genach regulowanych w górę. Im ciemniejszy odcień zieleni lub czerwieni, tym większa zmiana fałdowania. Ten rysunek został zmodyfikowany z Beheshti et al.36. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7: Określanie "głównego regulatora" dla gryzoni w pomieszczeniach dla gryzoni w porównaniu do klatek dla gryzoni. Połączenia między wszystkimi poszczególnymi kluczowymi genami (Ryc. 6) zostały określone i wyświetlone jako sieć za pomocą IPA. Sieć jest reprezentowana jako wykres promienisty z najbardziej połączonym kluczowym genem, MAPK1, w środku. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Uzupełniający Rysunek 1: Integracja GeneLab-GenomeSpace z ISACreator w celu usprawnienia operacji przetwarzania danych. Kliknij tutaj, aby pobrać ten rysunek.
Uzupełniające Rysunek 2: Zrzut ekranu wyszukiwania GeneLab przy użyciu federacji/integracji z heterogenicznymi zewnętrznymi bazami danych bioinformatycznych (GEO, PRIDE, MG-RAST). Kliknij tutaj, aby pobrać ten rysunek
Uzupełniające Rysunek 3: Zrzut ekranu przedstawiający wspólną przestrzeń roboczą GeneLab pokazującą zarządzanie kontami użytkowników i kontrolę dostępu (np. foldery prywatne, współdzielone, publiczne). Kliknij tutaj, aby pobrać ten rysunek