Dzięki analizie lokomocji w cieczy (pływania) można wyjaśnić fenotypy, które nie są łatwo zauważalne na podłożach stałych (pełzania). W celu ilościowego określenia lokomocji pływnej opracowaliśmy specjalistyczne oprogramowanie, które mierzy dziesięć nowych parametrów zachowania podczas pływania8. Osiem z tych parametrów, najbardziej użytecznych, zostało szczegółowo opisanych w Tabeli 1. Parametry te to: częstotliwość inicjacji fali (Wave initiation rate), liczba fal ciała (Body wave number), asymetria (Asymmetry), rozciągnięcie (Stretch), zwijanie (Curling), prędkość przemieszczania się (Travel speed), ruch pędzla (Brush stroke) oraz indeks aktywności (Activity index). Badania wykazujące skuteczność oprogramowania pozwoliły określić spadek funkcji u setek starzejących się osobników z genotypem WT, mutantami behawioralnymi lub mutantami długowieczności8, a także przeanalizować dobrze zbadane mutanty długowieczności age-1(hx546) oraz daƒ-16(mgDƒ50), które posiadają mutacje zaburzające prawidłową ścieżkę sygnalizacji insulinowej. Gen age-1 koduje podjednostkę katalityczną kinazy 3-fosfatydyloinozytolu (PIK3) i w przypadku wystąpienia mutacji hx546 powoduje wydłużenie czasu życia oraz odporność na stres27-29. Gen daƒ-16 koduje czynnik transkrypcyjny z rodziny forkhead box O (FOXO), którego delecja skraca czas życia i upośledza odpowiedź stresową30-33.
Niektóre parametry pływania, takie jak częstość inicjacji fal (Wave initiation rate), prędkość przemieszczania się (Travel speed), ruchy szczeczek (Brush stroke) oraz indeks aktywności (Activity index), stopniowo spadały wraz z wiekiem, nawet w sprzyjających liniach genetycznych (Rysunek 1). Zgodnie z aktualną wiedzą, długowieczne mutanty age-1(hx546) wykazywały wyższą sprawność fizyczną niż szczepy WT w wieku zaawansowanym i bardzo podeszłym. Zgodnie z przewidywaniami, krótkożyjące mutanty daƒ-16(mgDƒ50) wykazywały upośledzoną sprawność, szczególnie w bardzo podeszłym wieku. Co niezwykłe, dopiero dzięki zastosowaniu pakietu wizji komputerowej i algorytmów matematycznych CeleST można było wykryć lepszą sprawność pływania mutantów age-1(hx546) już na początku dorosłości. Fakt, że age-1(hx546) prowadzi do zwiększenia sprawności fizycznej w młodej dorosłości, sugeruje, że mutacja ta wpływa na prawidłowy rozwój i/lub fenotyp młodych osobników dorosłych w sposób, który nie był wcześniej dostrzeżony (Rysunek 1).
Parametry liczby fal ciała, asymetrii, rozciągania i zwijania wykazywały tendencję wzrostową wraz z wiekiem u osobników WT oraz dorosłych mutantów starzejących się (Rysunek 2). Co ciekawe, poziom rozdzielczości oprogramowania ujawnił subtelniejsze cechy behawioralne, takie jak utrzymująca się symetria mutantów age-1(hx546) przez cały okres ich życia oraz niezdolność bardzo starych mutantów daƒ-16(mgDƒ50) do rozciągania się i zwijania w takim stopniu, w jakim robią to osobniki WT i dorosłe mutanty age-1(hx546) w tym samym wieku.
Oprócz ogólnego i nieuniknionego spadku sprawności fizycznej związanego z wiekiem, każdy dorosły osobnik wykazuje unikalny wzorzec progresji w procesie starzenia, nawet w przypadku niemal jednorodnego tła genetycznego i środowiskowego7. (Poprzez kontrolę genetyki i środowiska zminimalizowano możliwe zakłócające działanie tych czynników, co ujawniło istotny wkład stochastyczności w degenerację związaną z wiekiem). Zsynchronizowana populacja C. elegans o podobnym tle genetycznym, utrzymywana w kontrolowanym środowisku, nadal zawiera mieszankę różnych klas osobników w zależności od ich specyfiki starzenia się. Chociaż wszystkie zaczynają jako zdrowe osobniki dorosłe, niektóre szybko tracą sprawność fizyczną (tzw. „bad agers”, klasa C), podczas gdy inne zachowują witalność przez dłuższy czas (tzw. „graceful agers”, klasa A). Osobniki z grupy „bad agers” wykazują zatem znacznie krótszy okres zdrowia (healthspan) niż osobniki z grupy „graceful agers”.
Jak szczegółowo opisano w naszym badaniu8, osobniki starzejące się w sposób łagodny (graceful agers) zachowały młodzieńczą sprawność fizyczną, co zaobserwowano poprzez porównanie z profilem pływania znacznie młodszych dorosłych (Rysunki 3, 4 i 5). Ta utrzymująca się sprawność jest porównywalna z wydolnością fizyczną długowiecznych mutantów age-1(hx546) w wieku pozaprodukcyjnym (D 11) (Rysunki 1 i 2). Przeciwnie, osobniki starzejące się w sposób gwałtowny (bad agers) drastycznie utraciły znaczną część swoich możliwości fizycznych wkrótce po reprodukcji, osiągając poziomy zbliżone do ekstremalnie starych oraz progerycznych dorosłych daƒ-16(mgDƒ50) (Rysunki 1 - 4). Podobieństwa te można wywieść z ogólnego porównania, jednak przy dokładniejszej analizie łatwo dostrzec konkretne sygnatury. Na przykład, choć istnieje pewna korelacja między zakresem rozciągania i zwijania u ekstremalnie starych typów dzikich oraz mutantów starzejących się (Rysunek 2), zależność ta nie jest obserwowana u osobników starzejących się w sposób gwałtowny (Rysunki 4 i 5), które w badanych próbkach wykazują większą tendencję do zwijania się, ale nie do rozciągania. Prezentowane przez nas oprogramowanie dodaje zatem nowy wymiar do analizy sprawności fizycznej lub zdolności lokomotorycznych, dostarczając narzędzi do bardziej zaawansowanych badań, które wcześniej nie były wykonalne. Podsumowując, CeleST zapewnia kompleksowe odczyty w formie ośmiu nowych miar omówionych w niniejszej pracy, które definiują behawioralny odcisk palca specyficznych uwarunkowań genetycznych, epigenetycznych i środowiskowych, umożliwiając identyfikację unikalnych i wspólnych wzorców parametrów, które mogą być sygnaturami konkretnych warunków (środowiskowych, farmakologicznych, żywieniowych), procesów biologicznych lub stanów organizmu, takich jak okres zdrowia (healthspan).

Rysunek 1: Raporty oprogramowania CeleST dotyczące szybkości inicjacji fal (A), indeksu aktywności (B), ruchu pędzla (C) i prędkości przemieszczania się (D) dla osobników WT, age-1(hx546), oraz daƒ-16(mgDƒ50) dorosłych w dniu 4 (wczesna dorosłość), 11 (po okresie reprodukcyjnym) i 20 (skrajna starość). Znak '#' na osi y oznacza 'liczbę'. WT zaznaczono kolorem szarym, age-1 zielonym, a daƒ-16 czerwonym. Słupki błędów reprezentują standardowy błąd średniej (SEM). W celu oceny istotności statystycznej porównano osobniki WT i mutanty starzejące się w tym samym wieku, stosując jednoczynnikową analizę wariancji (ANOVA), a następnie test wielokrotnych porównań Dunnetta. **, p = 0,001 - <0,01; ***, p = 0,0001 - <0,001. n = 62 dla każdego punktu danych z czterech niezależnych prób. Należy zauważyć, że tutaj oraz na Rysunku 2 każde pojedyncze 30 s nagranie wideo obejmuje 4 zwierzęta; dla każdej próby oceniamy łącznie 16 zwierząt z 4 filmów z pływania, co wykonuje się dla 4 powtórzeń biologicznych dla każdego przedstawionego punktu danych. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 2: Raporty programistyczne dotyczące liczby fal ciała (A), asymetrii (B), rozciągnięcia (C) i zwinięcia (D) dla osobników dorosłych WT, age-1(hx546) oraz daƒ-16(mgDƒ50) w dniu 4 (wczesna dorosłość), 11 (po okresie rozrodczym) i 20 (skrajna starość). Znak '#' na osi y oznacza 'liczbę'. WT zaznaczono kolorem szarym, age-1 zielonym, a daf-16 czerwonym. Słupki błędów przedstawiają błąd standardowy średniej (SEM). WT oraz mutanty starzejące się w tym samym wieku porównano pod kątem istotności statystycznej za pomocą jednoczynnikowej analizy wariancji (ANOVA), a następnie testu wielokrotnych porównań Dunnetta. *, p = 0,01 - <0,05; **, p = 0,001 - <0,01; ***, p = 0,0001 - <0,001. n = 62 dla każdego punktu danych z czterech niezależnych, 30 s prób pływania. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Raporty programowe dotyczące częstotliwości inicjacji fal (A), indeksu aktywności (B), ruchu szczecin (C) oraz prędkości pływania (D) dla młodych dorosłych WT (D 4) oraz osób w tym samym wieku z grupy „graceful agers” i „bad agers” (D10 i 11). Znak '#' na osi y oznacza 'liczbę'. Młode osobniki WT zaznaczono kolorem szarym, „graceful agers” klasy A kolorem zielonym, a „bad agers” klasy C kolorem czerwonym. Słupki błędów reprezentują błąd standardowy średniej (SEM). „Graceful agers” klasy A oraz „bad agers” klasy C porównywano z młodymi dorosłymi D 4 za pomocą jednoczynnikowej analizy wariancji (ANOVA), a następnie testu wielokrotnych porównań Dunnetta. ****, p <0,0001. n = 27 dla każdego punktu danych z dwóch niezależnych 30-sekundowych prób pływania. Wykres został lekko zmodyfikowany z publikacji Restif et al. (2014)8, która została opublikowana na licencji Creative Commons Attribution (CC BY) http://creativecommons.org/licenses/by/4.0/. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Raporty oprogramowania dotyczące liczby fal ciała (A), asymetrii (B) oraz zwijania się (C) dla młodych osobników WT (D 4) oraz osobników w tym samym wieku z grupy „graceful agers” i „bad agers” (D 10 i 11). Znak '#' na osi y oznacza „liczbę”. Młode osobniki WT zaznaczono kolorem szarym, „graceful agers” klasy A kolorem zielonym, a „bad agers” klasy C kolorem czerwonym. Słupki błędów przedstawiają standardowy błąd średniej (SEM). Grupy „graceful agers” klasy A oraz „bad agers” klasy C porównano do młodych osobników D 4 za pomocą jednoczynnikowej analizy wariancji (ANOVA), a następnie testu wielokrotnych porównań Dunnetta. **, p = 0,001 - <0,01; ****, p <0,0001; n/a, nie dotyczy, ponieważ tylko jedno zwierzę z całej próby zwijało się. n = 27 dla każdego punktu danych z dwóch niezależnych, 30 s prób pływania. Wykres został nieznacznie zmodyfikowany z pracy Restif et al. (2014)8, która została opublikowana na licencji Creative Commons Attribution (CC BY) http:/creativecommons.org/licenses/by/4.0/. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Raport programistyczny dotyczący rozciągliwości dla młodych osób z grupy WT (D 4) oraz osób w tym samym wieku z grupy „graceful agers” (D 10) i „bad agers” (D 11). Młode osoby WT zaznaczono kolorem szarym, „graceful agers” klasy A kolorem zielonym, a „bad agers” klasy C kolorem czerwonym. Słupki błędów przedstawiają standardowy błąd średniej (SEM). Osoby „graceful agers” klasy A oraz „bad agers” klasy C porównano z młodymi dorosłymi z grupy D 4 przy użyciu jednoczynnikowej analizy wariancji (ANOVA), a następnie testu wielokrotnych porównań Dunnetta. n = 27 dla każdego punktu danych z dwóch niezależnych serii doświadczeń. Wykres został nieznacznie zmodyfikowany z pracy Restif et al. (2014)8, która została opublikowana na licencji Creative Commons Attribution (CC BY) http:/creativecommons.org/licenses/by/4.0/. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Wideo 1: Pływanie reprezentatywnej grupy dorosłych osobników C. elegans . Kliknij tutaj, aby obejrzeć to wideo. (Kliknij prawym przyciskiem myszy, aby pobrać).*

Wideo 2: Obliczenia indywidualnych map krzywizny wyników pływania badanych zwierząt za pomocą oprogramowania CeleST. Mapy krzywizny są obliczane w tle; nie są one widoczne w interfejsie użytkownika oprogramowania. Kliknij tutaj, aby obejrzeć ten film. (Kliknij prawym przyciskiem myszy, aby pobrać).

Wideo 3: Obliczanie parametrów pływania przez oprogramowanie na podstawie indywidualnych map krzywizny. Kliknij tutaj, aby obejrzeć to wideo. (Kliknij prawym przyciskiem myszy, aby pobrać).*

Wideo 4: Obliczenia programowe parametrów pływania niewymagające map krzywizny. Kolej kliknij tutaj, aby obejrzeć to wideo. (Kliknij prawym przyciskiem myszy, aby pobrać).*