Name: Wiekowa dynamika lokomocji u Caenorhabditis elegans: analiza wykładnika Lyapunova (Wideo) | JoVE Journal
Uploaded: 2025-09-23T14:45:00.000+00:00
Duration: 6 min 44 s
Description: Odkryj analizę wiekowej dynamiki lokomocji u Caenorhabditis elegans za pomocą wykładnika Lyapunova. Metoda ta pozwala badać zmiany w kontroli ruchowej związane ze starzeniem się organizmu.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Dziękujemy Vassar College oraz Lucy Maynard Salmon Research Fund za wsparcie finansowe. Dziękujemy również dr Kathleen Susman, dr Juanowi Merlo oraz dr Susannah Zhang za ich wgląd i wsparcie na wszystkich etapach tych badań.

1. Przygotowanie C. elegans do akwizycji danych
<ol>
	<li>Umieść 0,5 mL E. coli z OD600 na każdej talerzu agarowej Nematode Growth Medium (NGM), aby nicienie mogły je zjeść. Poczekaj, aż E. coli wyschnie na każdym talerzu agaru.</li>
	<li>Zdobyć płytkę kontrolną C. elegans od dowolnego dostawcy materiałów biologicznych, aby stworzyć nowe płytki C. elegans w celu kontroli wieku.</li>
	<li>Wysterylizuj platynowy kostkę za pomocą płomienia z palnika Bunsena lub podobnego. Użyj mikroskopu rozcinającego, aby wybrać 5-10 dorosłych dzikich C. elegans do umieszczenia w każdej płytki.</li>
	<li>Pozwól nicieniom składać jaja przez 4-5 godzin, zanim zbiorzesz dorosłe osobniki z talerza. Wysiaduj pozostałe nicienie i pozwól im je zbierać w wybranym dniu.</li>
	<li>W dniu zbierania danych napełnij optyczną kwarcową kuwetę o pojemności 4,5 mL o wymiarach 10 mm x 10 mm x 45 mm wodą destylowaną w temperaturze pokojowej tuż poniżej górnej części kuwety, aby zapobiec rozlaniu się po dodaniu plastikowej pokrywy. Zbierz 2-3 nicienie i delikatnie umieść je w kuwetcie. Gdy nicienie są już w cuwette, połóż ją na boku, aby łatwo ustawić robaka w wiązce laserowej. C. elegans może unosić się ku dnu, ale będzie kontynuował falowania pływania, dopóki jest całkowicie zanurzonyw wodzie.</li>
</ol>
2. Zbieranie danych
UWAGA: Poniższa procedura rejestrowania szeregu czasu musi być wdrażana każdego dnia o tej samej porze, aby zminimalizować niepewności wieku.
<ol>
	<li>Uruchom eksperyment opisany w poprzednim artykule22. Wprowadź niewielką modyfikację: zbuduj peryskop z dwoma lustrzanymi lustrami i umieścisz kuwettę na boku między nimi, aby ułatwić późniejsze wycentrowanie robaka w wiązce laserowej. Zamiast aparatu, umieść PD w wzorze dyfrakcyjnym (Rysunek 1).<ol>
		<li>Ustawić lustra pionowo, aby zbudować wspomniany peryskop. UWAGA: To jedyne prace związane z laserowym wyrównaniem, które trzeba wykonać.</li>
	</ol>
	</li>
	<li>Włącz laser HeNe (Helium Neon) i pozwól mu się rozgrzać, osiągając równowagę termiczną (~15 min).</li>
	<li>Uruchom cyfrowy oscyloskop, aby rozpocząć zbieranie danych. Ustaw parametry interwału czasowego i bufora pamięci na cyfrowym oscyloskopie. Ustaw interwał czasowy na 100 s i zaimplementuj bufor pamięci o co najmniej 100 kilopróbkach na sekundę (kS).</li>
	<li>Ustaw rozdzielczość na 1 kHz dla szybkości akwizycji danych i 12 bitów dla rozdzielczości amplitud, co pomaga rozróżnić złożoność na małych skalach23.</li>
	<li>Ustaw oscyloskop tak, aby centrował oscylacje intensywności, używając przesunięcia Auto AC , aby wycentrować szeregi czasowe przy zerze woltów.</li>
	<li>Wybierz 2-3 nicienie i delikatnie umieść w kuwetcie wypełnionej wodą destylowaną, aby łatwiej było znaleźć i wycentrować przynajmniej jednego robaka w promieniu laserowym. Nie potrząsaj kuwetą.</li>
	<li>Umieść kuwetę z 2-3 nicieniami w peryskopie i wycentruj 1 nicienie w wiązce laserowej.<ol>
		<li>Gdy C. elegans jest wycentrowany w wiązce laserowej, około 50 cm od peryskopu tworzy się wzór dyfrakcji dalekiego pola.</li>
	</ol>
	</li>
	<li>Umieść fotodiodę w wzorze dyfrakcyjnym dalekiego pola, gdy C. elegans przechodzi przez wiązkę lasera. Upewnij się, że PD jest ustawiony poza środkiem w wzorze dyfrakcyjnym, aby wychwycić światło dyfrakcyjne, a nie centralne maksimum (przepuszczona wiązka laserowa).</li>
	<li>Zbieraj co najmniej 10 sekund danych, aby wiarygodnie obliczyć LLE; tzn. do obliczenia stabilnego LLE potrzebnych jest co najmniej 10 000 punktów danych. UWAGA: Mniej punktów danych powoduje sztuczne wahania LLE.</li>
	<li>Powtarzaj powyższe dane dla każdego dnia 9-15 razy. Zbieraj dane dla osób w wieku 3-12 dni o tej samej porze dnia.</li>
</ol>
3. Analiza danych
UWAGA: Do analizy danych szeregi czasowe są osadzane w przestrzeni fazowej za pomocą wykresulag 16 , a następnie szacują LLE, obliczając dywergencję trajektorii.
<ol>
	<li>Wybierz fragmenty w szeregu czasowym, w którym nicienie swobodnie pływa wewnątrz wiązki laserowej, dokładnie analizując szeregi czasowe. Rysunek 2A pokazuje żywotną serię czasową z ciągłym sygnałem co najmniej 20 s.</li>
	<li>Przefiltruj szeregi czasowe pod kątem niskiego stosunku sygnału do szumu (Rysunek 2B). Ustal, poziom szumów nagrywając szereg czasowy bez robaka w kuwecie. Cechy szeregu czasowego, nie tylko amplituda, muszą mieć ponad dwukrotnie wyższy poziom szumu.</li>
	<li>Wyklucz wszelkie nasycone dane z odpowiedniego szeregu czasowego (rysunek 2C). Jeśli jest zbyt dużo danych, aby wydobyć przydatne informacje, wróć do sekcji Zbieranie danych powyżej, powtórz kroki od 2.7.1 do 2.9 i dostosuj poziom intensywności, odsuwając fotodiodę od centralnego maksimum.</li>
	<li>Określmy średnią częstotliwość daną przez 
	<img alt="Equation 2" src="/files/ftp_upload/68955/68955eq2.jpg" style="margin: auto;" />   (2) gdzie f to częstotliwość, a P to widmo mocy (Rysunek 3). Można to wykonać obliczeniowo, używając szybkiej transformaty Fouriera (FFT), a następnie uśredniając częstotliwości. Wiele programów obliczeniowych posiada wbudowaną funkcję do obliczania średniej częstotliwości zbioru danych.</li>
	<li>Odbuduj topologię przestrzeni fazowej za pomocą metody opóźnienia czasowego. Wybierz szereg czasowy X i opóźnij ten ciąg o czas τ , aby rozstrzygnąć trajektorie (Rysunek 4). Optymalne opóźnienie czasowe określa się poprzez zidentyfikowanie pierwszego lokalnego minimum w informacji wzajemnej (MI)24: 
	<img alt="Equation 3" src="/files/ftp_upload/68955/68955eq3.jpg" style="margin: auto;" />    (3) gdzie N to liczba punktów, Xi to punkt w szeregu czasowym X, a Xi+τ to punkt opóźniony w tym samym szeregu czasowym. p(Xi) to prawdopodobieństwo wystąpienia punktu X i, a p(Xi+τ) to prawdopodobieństwo wystąpienia Xi+τ. p(Xi, Xi+τ) to wspólne prawdopodobieństwo25 zarówno ciągu czasowego Xi, jak i opóźnionego ciągu czasowego Xi+τ pokrywających się.</li>
	<li>Zidentyfikuj pierwsze minimum na wykresie MI między 0,140 a 0,240 s (140 do 240 punktów danych), jak pokazano na rysunku 5. Nie ma dokładnej liczby opóźnienia; Głównym celem jest rozróżnienie trajektorii fazowych na tyle, by określić dywergencję.</li>
	<li>Użyj szeregu czasowego Xi oraz jego opóźnionych wersji przez τ, co ustalano przez minimalizację MI w poprzednim kroku.</li>
	<li>Określ wymiar osadzenia. Najbardziej odpowiednim wymiarem zanurzenia jest najniższy wymiar, dla którego fałszywi najbliżsi sąsiedzi (FNN) stabilizują się przy minimum (Rysunek 6)26. Stosuj metodę opracowaną przez Abarbanela i in.27,28.</li>
	<li>Użyj algorytmu Rosensteina29 do rekonstrukcji układu dynamicznego w przestrzeni fazowej i śledzenia rozbieżności pobliskich trajektorii w czasie. W szczególności wykorzystaj procedurę MATLAB Merve Kizilkaya, opublikowaną na forumMATLAB 30, aby obliczyć LLE. Ta procedura wymaga następujących danych wejściowych: szeregu czasowego, średniej częstotliwości atraktora, opóźnienia czasowego τ, szybkości pozyskiwania danych dla szeregu czasowego oraz przedziału dopasowania do dywergencji.</li>
	<li>Powtórz powyższą procedurę dla zbiorów danych 9-15, aby zmniejszyć niepewność szacunków LLE. Uśrednij LLE dla każdego dnia, a następnie dopasuj do krzywej Moore'a (Rysunek 7).</li>
</ol>

<ol><li>Zhen, M., Samuel, A. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((C.+elegans+locomotion%3A+small+circuits%2C+complex+functions%5BTitle%5D)+AND+%22Curr+Opin+Neurobiol%22%5BJournal%5D)">C. elegans locomotion: small circuits, complex functions.</a> Curr Opin Neurobiol. 33, 117-126 (2015).</li>
<li>Corsi, A. K., Wightman, B., Chalfie, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+transparent+window+into+biology%3A+a+primer+on+Caenorhabditis+elegans%5BTitle%5D)+AND+%22WormBook%22%5BJournal%5D)">A transparent window into biology: a primer on Caenorhabditis elegans.</a> WormBook. , http://www.wormbook.org (2015).</li>
<li>Gjorgjieva, J., Biron, D., Haspel, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Neurobiology+of+Caenorhabditis+elegans+locomotion%3A+where+do+we+stand%5BTitle%5D)+AND+%22Bioscience%22%5BJournal%5D)">Neurobiology of Caenorhabditis elegans locomotion: where do we stand.</a> Bioscience. 64 (6), 476-486 (2014).</li>
<li>Boyle, J. H. C. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBoyle%2C+JHC+%22elegans+locomotion%3A+an+integrated+approach%22">elegans locomotion: an integrated approach</a>. , https://core.ac.uk/download/pdf/43059.pdf (2009).</li>
<li>Pierce-Shimomura, J. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Genetic+analysis+of+crawling+and+swimming+locomotory+patterns+in+C.+elegans%5BTitle%5D)+AND+%22Proc+Natl+Acad+Sci+U+S+A%22%5BJournal%5D)">Genetic analysis of crawling and swimming locomotory patterns in C. elegans.</a> Proc Natl Acad Sci U S A. 105, 20982-20987 (2008).</li>
<li>Korta, J., Clark, D. A., Gabel, C. V., Mahadevan, L., Samuel, A. D. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Mechanosensation+and+mechanical+load+modulate+the+locomotory+gait+of+swimming+C.+elegans%5BTitle%5D)+AND+%22J+Exp+Biol%22%5BJournal%5D)">Mechanosensation and mechanical load modulate the locomotory gait of swimming C. elegans.</a> J Exp Biol. 210 (13), 2383-2389 (2007).</li>
<li>Edwards, S. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+novel+molecular+solution+for+ultraviolet+light+detection+in+Caenorhabditis+elegans%5BTitle%5D)+AND+%22PLoS+Biol%22%5BJournal%5D)">A novel molecular solution for ultraviolet light detection in Caenorhabditis elegans.</a> PLoS Biol. 6 (8), e198(2008).</li>
<li>Barbulescu, R., Mestre, G., Oliveira, A. L., Silveira, L. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Learning+the+dynamics+of+realistic+models+of+C.+elegans+nervous+system+with+recurrent+neural+networks%5BTitle%5D)+AND+%22Sci+Rep%22%5BJournal%5D)">Learning the dynamics of realistic models of C. elegans nervous system with recurrent neural networks.</a> Sci Rep. 13 (1), 467(2023).</li>
<li>Sarma, G. P., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((OpenWorm%3A+overview+and+recent+advances+in+integrative+biological+simulation+of+Caenorhabditis+elegans%5BTitle%5D)+AND+%22Philos+Trans+R+Soc+Lond+B+Biol+Sci%22%5BJournal%5D)">OpenWorm: overview and recent advances in integrative biological simulation of Caenorhabditis elegans.</a> Philos Trans R Soc Lond B Biol Sci. 373 (1758), 20170382(2018).</li>
<li>Zanetti, R. F., Canavan, K. L., Zhang, S. G., Magnes, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Multichannel+measurements+of+C.+elegans+largest+Lyapunov+exponents+using+optical+diffraction%5BTitle%5D)+AND+%22Appl+Opt%22%5BJournal%5D)">Multichannel measurements of C. elegans largest Lyapunov exponents using optical diffraction.</a> Appl Opt. 62 (29), 7812-7818 (2023).</li>
<li>Magnes, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Live+C.+elegans+diffraction+at+a+single+point%5BTitle%5D)+AND+%22Open+J+Biophys%22%5BJournal%5D)">Live C. elegans diffraction at a single point.</a> Open J Biophys. 8, 155-162 (2018).</li>
<li>Lorenz, E. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+predictability+of+hydrodynamic+flow%5BTitle%5D)+AND+%22Trans+N+Y+Acad+Sci+Ser+II%22%5BJournal%5D)">The predictability of hydrodynamic flow.</a> Trans N Y Acad Sci Ser II. 25 (4), 409-432 (1963).</li>
<li>Magnes, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chaotic+markers+in+dynamic+diffraction%5BTitle%5D)+AND+%22Appl+Opt%22%5BJournal%5D)">Chaotic markers in dynamic diffraction.</a> Appl Opt. 59 (22), 6642-6650 (2020).</li>
<li>Strogatz, S. H. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aStrogatz%2C+SH+%22Nonlinear+dynamics+and+chaos%3A+with+applications+to+physics%2C+biology%2C+chemistry%2C+and+engineering%22">Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering</a>. , Westview Press. Boulder, CO. (2001).</li>
<li>Nolte, D. D. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aNolte%2C+DD+%22Introduction+to+modern+dynamics%3A+chaos%2C+networks%2C+space+and+time%22">Introduction to modern dynamics: chaos, networks, space and time</a>. , Oxford University Press. Oxford. (2015).</li>
<li>Takens, F. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aTakens%2C+F+%22Detecting+strange+attractors+in+turbulence%22">Detecting strange attractors in turbulence</a>. , Springer. Berlin. (1981).</li>
<li>Murakami, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Manipulation+of+serotonin+signal+suppresses+early+phase+of+behavioral+aging+in+Caenorhabditis+elegans%5BTitle%5D)+AND+%22Neurobiol+Aging%22%5BJournal%5D)">Manipulation of serotonin signal suppresses early phase of behavioral aging in Caenorhabditis elegans.</a> Neurobiol Aging. 29 (7), 1093-1100 (2008).</li>
<li>Marck, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Age-related+changes+in+locomotor+performance+reveal+a+similar+pattern+for+Caenorhabditis+elegans%2C+Mus+domesticus%2C+Canis+familiaris%2C+Equus+caballus%2C+and+Homo+sapiens%5BTitle%5D)+AND+%22J+Gerontol+A+Biol+Sci+Med+Sci%22%5BJournal%5D)">Age-related changes in locomotor performance reveal a similar pattern for Caenorhabditis elegans, Mus domesticus, Canis familiaris, Equus caballus, and Homo sapiens.</a> J Gerontol A Biol Sci Med Sci. 72 (4), 455-463 (2017).</li>
<li>Cohen, N., Sanders, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Nematode+locomotion%3A+dissecting+the+neuronal-environmental+loop%5BTitle%5D)+AND+%22Curr+Opin+Neurobiol%22%5BJournal%5D)">Nematode locomotion: dissecting the neuronal-environmental loop.</a> Curr Opin Neurobiol. 25, 99-106 (2014).</li>
<li>Moore, D. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+study+of+age+group+track+and+field+records+to+relate+age+and+running+speed%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">A study of age group track and field records to relate age and running speed.</a> Nature. 253 (5489), 264-265 (2020).</li>
<li>Varier, S., Kaiser, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Evolution+and+development+of+brain+networks%3A+from+Caenorhabditis+elegans+to+Homo+sapiens%5BTitle%5D)+AND+%22Network%22%5BJournal%5D)">Evolution and development of brain networks: from Caenorhabditis elegans to Homo sapiens.</a> Network. 22 (1-4), 143-147 (2011).</li>
<li>Magnes, J., Susman, K., Eells, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Quantitative+locomotion+study+of+freely+swimming+micro-organisms+using+laser+diffraction%5BTitle%5D)+AND+%22J+Vis+Exp%22%5BJournal%5D)">Quantitative locomotion study of freely swimming micro-organisms using laser diffraction.</a> J Vis Exp. (68), e4412(2012).</li>
<li>Tzepos, D., Trader, O., Magnes, J. 16th Chaotic Modeling and Simulation International Conference, , https://link.springer.com/book/10.1007/978-3-031-60907-7 (2024).</li>
<li>Fraser, A., Swinney, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Independent+coordinates+for+strange+attractors+from+mutual+information%5BTitle%5D)+AND+%22Phys+Rev+A%22%5BJournal%5D)">Independent coordinates for strange attractors from mutual information.</a> Phys Rev A. 33, 1134-1140 (1986).</li>
<li>Starmer, J. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aStarmer%2C+J+%22The+StatQuest+illustrated+guide+to+machine+learning%21%21%21%3A+triple+bam%22">The StatQuest illustrated guide to machine learning!!!: triple bam</a>. , StatQuest Publications. (2022).</li>
<li>Kizilkaya, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aKizilkaya%2C+M+%22MATLAB+routine%22">MATLAB routine</a>. , https://www.mathworks.com/matlabcentral (2020).</li>
<li>Abarbanel, H. D. I., Kennel, M. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Local+false+nearest+neighbors+and+dynamical+dimensions+from+observed+chaotic+data%5BTitle%5D)+AND+%22Phys+Rev+E%22%5BJournal%5D)">Local false nearest neighbors and dynamical dimensions from observed chaotic data.</a> Phys Rev E. 47, 3057-3068 (1993).</li>
<li>Kizilkaya, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aKizilkaya%2C+M+%22False+nearest+neighbor+algorithm%22">False nearest neighbor algorithm</a>. , https://www.mathworks.com/matlabcentral/fileexchange/37239-minimum-embedding-dimension (2020).</li>
<li>Rosenstein, M. T., Collins, J. J., De Luca, C. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+practical+method+for+calculating+largest+Lyapunov+exponents+from+small+data+sets%5BTitle%5D)+AND+%22Physica+D%22%5BJournal%5D)">A practical method for calculating largest Lyapunov exponents from small data sets.</a> Physica D. 65 (1-2), 117-134 (1993).</li>
<li>Kizilkaya, M. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aKizilkaya%2C+M+%22Largest+Lyapunov+exponent+with+Rosenstein%27s+algorithm%22">Largest Lyapunov exponent with Rosenstein's algorithm</a>. , https://www.mathworks.com/matlabcentral/fileexchange/38424-largest-lyapunov-exponent-with-rosenstein-s-algorithm (2020).</li>
<li>Liu, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Functional+aging+in+the+nervous+system+contributes+to+age-dependent+motor+activity+decline+in+C.+elegans%5BTitle%5D)+AND+%22Cell+Metab%22%5BJournal%5D)">Functional aging in the nervous system contributes to age-dependent motor activity decline in C. elegans.</a> Cell Metab. 18 (3), 392-402 (2013).</li>
<li>Zhang, S. G., Singhvi, A., Susman, K. M., Hastings, H. M., Magnes, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+markers+for+chaotic+motion+in+C.+elegans%5BTitle%5D)+AND+%22Nonlinear+Dyn+Psychol+Life+Sci%22%5BJournal%5D)">Dynamic markers for chaotic motion in C. elegans.</a> Nonlinear Dyn Psychol Life Sci. 26 (1), 21-43 (2022).</li>
<li>Giordano, N. J., Nakanishi, H. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aGiordano%2C+NJ+%22Computational+physics%22">Computational physics</a>. , Upper Saddle River, NJ. (2006).</li>
<li>Olsen, A., Vantipalli, M. C., Lithgow, G. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Using+Caenorhabditis+elegans+as+a+model+for+aging+and+age-related+diseases%5BTitle%5D)+AND+%22Ann+N+Y+Acad+Sci%22%5BJournal%5D)">Using Caenorhabditis elegans as a model for aging and age-related diseases.</a> Ann N Y Acad Sci. 1067 (1), 120-128 (2006).</li>
<li>Ahamed, T., Costa, A. C., Stephens, G. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Capturing+the+continuous+complexity+of+behaviour+in+Caenorhabditis+elegans%5BTitle%5D)+AND+%22Nat+Phys%22%5BJournal%5D)">Capturing the continuous complexity of behaviour in Caenorhabditis elegans.</a> Nat Phys. 17 (2), 275-283 (2021).</li>
<li>Zhang, S. G., Singhvi, A., Susman, K. M., Hastings, H. M., Magnes, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Dynamic+markers+for+chaotic+motion+in+C.+elegans%5BTitle%5D)+AND+%22Nonlinear+Dyn+Psychol+Life+Sci%22%5BJournal%5D)">Dynamic markers for chaotic motion in C. elegans.</a> Nonlinear Dyn Psychol Life Sci. 26 (1), 21-43 (2022).</li>
<li>Kennel, M., Brown, R., Abarbanel, H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Determining+embedding+dimension+for+phase-space+reconstruction+using+a+geometrical+construction%5BTitle%5D)+AND+%22Phys+Rev+A%22%5BJournal%5D)">Determining embedding dimension for phase-space reconstruction using a geometrical construction.</a> Phys Rev A. 45, 3403-3411 (1992).</li>
<li>Varier, S., Kaiser, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Evolution+and+development+of+brain+networks%3A+from+Caenorhabditis+elegans+to+Homo+sapiens%5BTitle%5D)+AND+%22Network%22%5BJournal%5D)">Evolution and development of brain networks: from Caenorhabditis elegans to Homo sapiens.</a> Network. 22 (1-4), 143-147 (2011).</li>
<li>Herndon, L. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Stochastic+and+genetic+factors+influence+tissue-specific+decline+in+ageing+C.+elegans%5BTitle%5D)+AND+%22Nature%22%5BJournal%5D)">Stochastic and genetic factors influence tissue-specific decline in ageing C. elegans.</a> Nature. 419, 808-814 (2002).</li>
<li>Rosso, O. A., Larrondo, H. A., Martin, M. T., Plastino, A., Fuentes, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Distinguishing+noise+from+chaos%5BTitle%5D)+AND+%22Phys+Rev+Lett%22%5BJournal%5D)">Distinguishing noise from chaos.</a> Phys Rev Lett. 99 (15), 154102(2007).</li>
</ol>

Autor nie ma nic do ujawnienia.

Kategoryzowaliśmy lokomocję C. elegans, analizując LLE w różnych wiekach, wspierając parametry takie jak częstotliwość pływania, wzajemna informacja (MI) oraz fałszywi najbliżsi sąsiedzi (FNN) – każdy z nich dostarcza dodatkowego wglądu w właściwości dynamiczne systemu. Deterministyczne i nieliniowe właściwości tego systemu badane są za pomocą metod danych zastępczych oraz wykresów rekurencji36.
Wykorzystanie DOD do ilościowego określenia wpływu starzenia się na lokomocję jest techniką uzupełniającą tradycyjną mikroskopię. Nie zastępuje to wizualnej inspekcji mikroskopijnych gatunków; Zapewnia spójny i efektywny sposób ilościowego określania lokomocji wykraczającej poza tradycyjne techniki. Ponieważ DOD omija analizę wideo, część obciążenia obliczeniowego jest z natury zarządzana przez technikę optyczną, ponieważ dyfrakcja bezpośrednio odwzorowuje rozkład intensywności na przestrzeń Fouriera. Superpozycja pola optycznego podczas dyfrakcji śledzi cały mikroskopijny gatunek gatunkowy. Idealnie rozdzielczość jest regulowana przez długość fali λ użytego światła, w tym przypadku 632 nm. Nawet przesunięcie o ułamek długości fali spowoduje zmianę intensywności. Z tego powodu rozdzielczość zależy od zasięgu i rozdzielczości fotodetektora oraz mocy lasera. Na przykład, jeśli detektor wykrywa sygnał szeregu czasowego, który waha się między dwoma napięciami, Vmin i Vmax, a rozdzielczość wynosi n bitów (patrz krok protokołu 2.4), to maksymalna rozdzielczość wynosi λ/n, czyli w naszym przypadku 632 nm/12. Teoretycznie metoda ta może być stosowana do każdego ciała, które zmienia swój kształt; jednak DOD jest szczególnie dobrze przystosowany do gatunków mikroskopijnych, ponieważ tanie lasery o odpowiednich długościach fal są łatwo dostępne i opłacalne.
Aby zapewnić wyraźny i spójny wzór dyfrakcyjny, robaka należy dokładnie śledzić i utrzymywać w centrum wiązki lasera podczas pływania. Jeśli robak wyjdzie z ustawienia, sygnał może się pogorszyć lub całkowicie utracić. Aby złagodzić problemy nasycenia fotodiodami i zmniejszyć szum stochastyczny, stosuje się filtr neutralnej gęstości, który osłabia intensywność lasera. To filtrowanie pomaga utrzymać zakres dynamiczny fotodiody i zapobiega przesyceniu w zarejestrowanym szeregu czasowym. Niemniej jednak, wszelkie segmenty danych wykazujące nadmierne nasycenie lub niskie stosunki sygnału do szumu są wyłączane z końcowej analizy, aby utrzymać jakość danych.
Dryf częstotliwości pływania tworzy ciągłe spektrum częstotliwości, które jest ustalonym wczesnym markerem chaosu13. Odwrotna zależność między średnią częstotliwością a MI sugeruje, że około 1/8 cyklu rozstrzyga trajektorię z wystarczającą dokładnością, aby oszacować LLE, podobnie jak rozdzielczość wymagana do analizy przebiegu sinusoidalnego.
Niski odsetek FNN wokół wymiarów osadzenia 3 i 4 na Rysunku 6 sugeruje niski poziom szumów zarówno w ruchu C. elegans, jak i w systemie optycznym. Jednak u 12-dniowych robaków odsetek FNN nieznacznie wzrasta, ale pozostaje znacznie poniżej 5%, niekoniecznie z powodu wzrostu prawdziwego wymiaru osadzenia, lecz prawdopodobnie ze względu na wzrost szumu w sygnale lokomotorycznym, przypisywany degradacji nerwowej u starzejących się organizmów34,37.
Zmiany w LLE w trakcie rozwoju pokazują, że stopień chaosu w lokomocji zmienia się wraz z wiekiem, co bardzo odpowiada przewidywaniom z równania Moore'a. Szczyt LLE we wczesnym okresie życia jest zgodny z cechami gatunków wybranych R, takich jak C. elegans, które muszą szybko rozwijać funkcje neuromięśniowe bez opieki rodzicielskiej. Co istotne, 70% długodystansowych połączeń neuronalnych powstaje, gdy robak osiąga zaledwie 20% dorosłego rozmiaru37,38. Natomiast gatunki wyselekcjonowane przez K, takie jak ludzie, wykazują wolniejszy rozwój neurologiczny, wspierany przez zaangażowanie rodzicielskie, co skutkuje opóźnionymi, ale bardziej złożonymi wzorcami ruchowymi. Różnice te odzwierciedlają się w czasie szczytowego LLE między gatunkami.
Krótki cykl życiowy C. elegans, około 14 dni, z szybkim rozwojem w ciągu pierwszych dwóch dni w czterech stadium larwalnych2, czyni go szczególnie odpowiednim do badania szczytowej złożoności w lokomocji. Podczas gdy wcześniejsze badania koncentrowały się na spadku struktury neuronów, aby wyjaśnić zmiany w lokomocji związane z wiekiem34,39, nasze badania łączą te fizyczne zmiany z chaotyczną dynamiką systemu. Integracja DOD z dynamiką nieliniową umożliwia precyzyjną ilościową charakterystykę na wielu skalach długości, oferując nowe wglądy w neuronalne czynniki napędzające zachowania. Stale pozytywny LLE na różnych etapach rozwoju silnie wspiera obecność chaotycznych zachowań w układzie lokomotorycznym.
Badanie to oferuje nowatorskie ramy metodologiczne do badania zmian zachowań związanych z wiekiem u C. elegans, poszerzając tym samym obecne rozumienie biologicznych i dynamicznych korelatów spadku ruchu. Dodatkowo potwierdza zastosowanie LLE jako wrażliwego biomarkera dla zmian neurologicznych związanych z wiekiem i wykazuje, że LLE jest wiarygodnym ilościowym wskaźnikiem dynamiki lokomotorycznej wpływającej na starzenie się.
DOD to potężna, nieinwazyjna metoda ilościowego określania lokomocji i chaosu u organizmów mikroskopijnych, najlepiej postrzegana jako uzupełnienie analizy wideo, a nie jej zamiennik. W obecnej wersji eksperymentu wymaga, aby nicienie miały co najmniej trzy dni, ponieważ młodsze robaki są zbyt małe, by dostarczać wiarygodne sygnały dyfrakcyjne. Jednowymiarowy szereg czasowy wyodrębniony z pojedynczego punktu w wzorze dyfrakcyjnym nieuchronnie kompresuje informacje przestrzenne, co może zaciemniać lokalną dynamikę na całym ciele organizmu. Podejścia wielokanałowe pomagają złagodzić to ograniczenie, weryfikując spójność parametrów w całym polu dyfrakcyjnym, choć subtelne heterogeniczności przestrzenne mogą być nadal pomijane. Dokładne oszacowanie LLE zależy dodatkowo od uzyskania wystarczająco długich, wolnych od szumów szeregów czasowych; praktyczne wyzwania, takie jak ruch robaka poza wiązką czy wahania środowiskowe, mogą obniżyć jakość danych i zmniejszyć zaufanie do wartości LLE. Jak we wszystkich metodach eksperymentalnych, szum pomiarowy i subiektywne decyzje podczas dopasowania LLE wprowadzają dodatkową niepewność, która dla pojedynczych pomiarów zazwyczaj utrzymuje się w okolicach 15%.
Przyszłe prace będą koncentrować się na pomiarze LLE w różnych warunkach eksperymentalnych, umożliwiając porównania między modelowanymi a mierzonymi LLE. Ułatwi to rozwój predykcyjnych modeli neurologicznych, które mogą pogłębić nasze zrozumienie, jak złożona kontrola ruchowa powstaje i pogarsza się z czasem. Będziemy również kontynuować badanie spójności naszych ustaleń, badając inne metody obliczeniowe, takie jak obliczanie płaszczyzny entropii-złożoności, aby ponownie zweryfikować deterministyczny charakter tego systemu biologicznego40.

Lokomomocja Caenorhabditis elegans (C. elegans), mikroskopijnego robaka, była badana w celu pogłębienia zrozumienia układu neuronów ruchowych, ponieważ nicienie to jest neurologicznie proste, zawiera tylko 302 neurony1. C. elegans to organizm modelowy, który łatwo się utrzymuje, żyje tylko 14 dni2. Do poruszania się służy tylko około 72 neuronów z grupy 302 u C. elegans i znajdują się w całym dorosłym ciele nicienia1. Ruch niskowymiarowy (czyli na boki, do przodu i do tyłu) C. elegans sprawia, że jest to łatwy okaz do śledzenia odległości w LLE. Układy nerwowe nicieni są rozumiane w niesamowitym stopniu4.
Analiza wideo pomogła określić ilość ruchu C. elegans poprzez pomiar wielkości takich jak promień krzywizny, częstotliwości falowania oraz długość fali5. Pomiary te ustaliły zmienne kontrolne umożliwiające porównanie właściwości lokomotorycznych w różnych środowiskachi warunkach 6,7. Te informacje tworzą ewoluujący model obwodów napędzających lokomocję C. elegans8, a nawet umożliwiają tworzenie wirtualnych, jak i fizycznych symulacji dynamicznych robaków9.
Dynamiczna dyfrakcja optyczna (DOD)10 została również wykorzystana do ilościowego określenia lokomocji C. elegans. Podczas DOD niskonatężene światło laserowe zakrzywia się wokół żywego robaka, tworząc wzór dyfrakcji dalekiego pola znany jako wzór dyfrakcji Fraunhofera. Rozkład intensywności w żywym wzorcu zmienia się wraz z ruchem nicienia. Jednym z punktów dyfrakcji jest superpozycja wszystkich punktów na robaku, tak że zależna od czasu intensywność w wzorcu dyfrakcyjnym tworzy jednowymiarowy szereg czasowy zawierający informacje o dynamice lokomotorycznej11. Skupiamy się na zależnym od wieku Największym Wykładniku Ljapunowa (LLE) lokomocji, który obliczamy na podstawie eksperymentalnego jednowymiarowego szeregu czasowego. LLE na różnych etapach rozwoju jest porównywalny z innymi badaniami dotyczącymi wieku, pokazując, że obwody nerwowe przejawiają się w wzorcach lokomotorycznych, które można mierzyć za pomocą różnych narzędzi, takich jak analiza wideo czy DOD5.
DOD zapewnia szeregi czasowe wrażliwe na wiele skal, ponieważ wzór interferencji może rozłożyć ruch do ułamka używanej długości fali, a jednocześnie uwzględniać duże zmiany plastyczności (kształtu) gatunku. Ta cecha jest szczególnie przydatna, gdy układ jest wrażliwy na warunki początkowe, jak w systemach chaotycznych, gdzie niewielkie zmiany prowadzą do wykładniczych zmian trajektorii, znanych również jako efekt motyla12. Z tego powodu kluczowe jest rejestrowanie szeregów czasowych w tempie umożliwiającym rejestrowanie drobnych zmian, ponieważ te zmiany mogą powodować znaczące przesunięcia w miarę ewolucji systemu. Fotodioda (PD) może rejestrować szeregi czasowe szybciej niż wiele drogich kamer o dużej prędkości. Połączenie czułości przestrzennej w wzorze dyfrakcyjnym i wysokiej szybkości gromadzenia danych może oddać istotę chaotycznego układu13,14.
LLE opisuje wykładniczą dywergencję trajektorii przestrzeni fazowej. Przestrzeń fazowa zawiera wszystkie możliwe stany układu, opisane przez ich zmienne stanowe (lub współrzędne) oraz powiązane pochodne (lub pędy)15. Różne trajektorie opisują ewolucję układu dla różnych warunków początkowych. W układzie chaotycznym dwie pobliskie trajektorie będą się wykładniczo rozbiegać w czasie, a ta dywergencja jest kwantyfikowana przez dodatnią LLE. W układach fizycznych ta dywergencja objawia się propagowaniem niepewności w warunkach początkowych układu, a także niemożnością wiarygodnego przewidzenia stanu systemu po pewnym czasie, jak daje LLE14.
W systemach eksperymentalnych zmienne stanu są często nieznane; jednak topologię przestrzeni fazowej można odtworzyć tylko z jednej zmierzonej zmiennej, korzystając z twierdzenia16 osadzonego Takensa, konstruując jej wersje opóźnione czasowo (równoważne pochodnym czasowym). Każda wersja z opóźnieniem czasowym jest reprezentowana na osi; Razem opóźnienia czasowe tworzą wykres opóźniający, który jest topologicznie identyczny z wykresem fazowym, jeśli szereg czasowy jest nieskończenie długi. Eksperymentalne szeregi czasowe są ograniczone długością i dlatego mogą generować jedynie szacowane LLE ograniczone eksperymentalnymi niepewnościami.
Ruch C. elegans niesie dodatni LLE10, co jest wskaźnikiem chaosu. Młody organizm przechodzi zmiany lokomotoryczne w wieku17 lat. Większy LLE oznacza mniejszą przewidywalność w porównaniu do mniejszego LLE. LLE dostarcza wiarygodnej metody ilościowego określania wzorców ruchowych organizmu 10,13,14. Cechy lokomotoryczne nicieni są powiązane z neuronami w obrębie nicieni18. Według Cohena i in.19, lokomocyja nicieni zależy od neuronów ruchowych, co sugeruje, że badanie złożoności w poruszaniu nicieni jest również powiązane z jego obwodami neuronalnymi.
Kalibrujemy LLE na podstawie wieku, aby zminimalizować niepewność na potrzeby przyszłych badań. Wcześniejsze badania uśredniły LLE u nicieni13 lat w wieku od trzech do sześciu dni. Badanie to starannie monitorowało wiek nicieni, aby ułatwić rygorystyczną analizę spadków związanych ze starzeniem. Dodatkowo badane są zmiany w poruszaniu się w celu identyfikacji innych biologicznych zmian u C. elegans18. Równanie Moore'a śledzi parametry wieku związane z lokomocją20. Tutaj używamy zmodyfikowanej wersji równania Moore'a, która pozwala na elastyczność w określaniu mechanizmu lokomotorycznego do oceny, wprowadzonej jako marker czasowy T:
<img alt="Equation 1" src="/files/ftp_upload/68955/68955eq1.jpg" style="margin: auto;" />, (1)
gdzie P(t) oznacza wydajność, t to czas, a i c to parametry skalowania, natomiast b i d to charakterystyczne czasy odpowiednio wykładniczego wzrostu i spadku. Wykazano, że równanie Moore'a przewiduje trajektorie pozornie niezwiązanych ze sobą wielkości, takich jak szybkość i aktywność dobrowolna, w miarę starzenia się organizmów. Gatunki można scharakteryzować na podstawie kształtu krzywej określonego przez równanie Moore'a; t. j. niektóre gatunki mogą osiągnąć szczyt wcześniej niż inne21.

<table><tbody><tr><td>2 aluminiowe lusterka z przednią powierzchnią</td><td>Thorlabs</td><td>PF10-03-F01</td><td></td></tr><tr><td>Laser HeNe 632 nm</td><td>Newport</td><td>LGX1</td><td>Każdy laser z czerwonym światłem</td></tr><tr><td>Pusta miska Petriego</td><td>Carolina</td><td>971632</td><td>Plastikowe szalki Petriego, do których wlewamy nicienie wzrostowe agaru</td></tr><tr><td>Escherichia coli K12, żywy, bakteriofagowy gospodarz</td><td>Carolina</td><td>124500</td><td>Wykorzystywany jako źródło pożywienia dla C. elegans; OD600</td></tr><tr><td>Mikroskop Leica S9i</td><td>Leica Microsystems</td><td>LED2500</td><td>Mikroskop rozcinający</td></tr><tr><td>Zapalniczka</td><td>Bic Lighters</td><td></td><td>Każde narzędzie sterylizujące używane do sterlizacji kilofa przed i po zerwaniu każdego robaka</td></tr><tr><td>MATLAB</td><td>MathWorks</td><td></td><td>Algorytm-Rutyna Rosensteina na forum MATLAB stworzony przez Merve Kizilkaya</td></tr><tr><td>Nicienie Growth Agar</td><td>Carolina</td><td>173520</td><td>Przygotowana butelka mediów, 135 mL</td></tr><tr><td>Fotodiody</td><td>Thorlabs</td><td>DET36A</td><td>Detektor Si o napięciu 350-1100 nm</td></tr><tr><td>Picoscope5204</td><td>Technologia Pico</td><td>PP376</td><td>Oscyloskop PC www.picotech.com</td></tr><tr><td>Platinum Pick</td><td></td><td></td><td>Wykorzystywany do wyławiania C. elegans; To mały, ręcznie robiony kostek z szklanym elementem trzymanym w ręku i platynową łyżką&nbsp;</td></tr><tr><td>Kwarcowa kuweta</td><td>Komórki Starna</td><td>21/G/5</td><td>Wypełnione wodą destylowaną, aby umieścić C. elegans w środku</td></tr></tbody></table>

wiekowa dynamika lokomocji u caenorhabditis elegans: analiza wykładnika lyapunova

Age-dependent Dynamics of Locomotion in Caenorhabditis elegans: A Lyapunov Exponent Analysis

Niniejsze badanie bada wpływ wieku na mobilność Caenorhabditis elegans (C. elegans) za pomocą dynamicznej dyfrakcji optycznej (DOD) do oszacowania największego wykładnika Lyapunova (LLE). LLE, kluczowa miara w układach dynamicznych, ilościowo określa szybkość dywergencji lub zbieżności trajektorii w przestrzeni fazowej, wskazując na przewidywalność i chaos w dynamice układu, w tym przypadku zachowania lokomotywy robaka. Światło laserowe o długości 632 nm odbija się od pływającego robaka w kolumnie wodnej, tworząc wzór dyfrakcyjny. Fotodioda wykrywa światło w jednym punkcie w obrębie wzoru dyfrakcyjnego, rejestrując jednowymiarowy szereg czasowy podczas falowania robaka. Ten szereg czasowy pełni funkcję złożonej reprezentacji całego ruchu robaka, obejmując jego dynamikę lokomocji, ponieważ jeden punkt w wzorze dyfrakcyjnym jest superpozycją wszystkich punktów na robaku. Szeregi czasowe są następnie osadzane w przestrzeni fazowej o wyższych wymiarach, aby obliczyć LLE. C. elegans zazwyczaj żyje około 14 dni, podążając wzorcem wzrastającej i malejącej kontroli motorycznej wraz z wiekiem. Aby wyizolować efekty specyficzne dla wieku, robaki były przenoszone na świeże płytki agarowe zawierające E. coli co dwa dni, zapewniając odpowiednie starzenie (od 3 do 12 dni). Analiza kohorty 13 C. elegans wykazała, że LLE osiągnęła szczyt pięć dni po wykluciu, z tempem 1,34 ± 0,03 1/s. Ten szczyt oznacza krytyczny punkt rozwoju, w którym ruch robaków wykazuje największą złożoność i chaotyczne zachowanie. Zaobserwowane wartości LLE są zgodne z równaniem Moore'a, dobrze ugruntowanym modelem opisującym zmiany aktywności dobrowolnej związane z wiekiem, łączące spadek kontroli ruchowej i poziomu aktywności z rosnącym wiekiem u C. elegans.

Mierząc LLE względem wieku, mierzyliśmy także inne wielkości często związane ze złożonymi układami dynamicznymi, takie jak częstotliwość atraktora, MI i FNN. Zebraliśmy dane z badań wiekowych od 3 do 12 dni. Robaki są zbyt małe i niedojrzałe, zanim skończą 3 dni, by ręcznie się nimi zajmować. Po 12 dniach robaki są starzejące się i prawie się nie ruszają, ponieważ ich obwody neuronalne się pogarszają, a układ nerwowy pogarsza sięo 31 dni.
Algorytm Rosensteina wybiera sąsiedni punkt w odstępie co najmniej jednego średniego okresu dla każdego punktu na osadzonej trajektorii. Zaczynając od wektora separacji między dwoma punktami na dwóch sąsiednich trajektoriach (Rysunek 8), algorytm śledzi ewolucję tej dywergencji. Odległość trajektorii początkowo rośnie wykładniczo z powodu chaotycznej natury układu32, a następnie ustępuje plateau, ponieważ trajektoria jest ograniczona. Wykreślenie logarytmu średniej rozbieżności w czasie śledzi rozdzielenie trajektorii (Rysunek 8 i Rysunek 9). Nachylenie liniowych najmniejszych kwadratów dopasowuje się do pierwszej części krzywej przed jej wyprostowaniem, co daje wiarygodne oszacowanie dla LLE33 układu. Należy zauważyć, że wybór zakresu dopasowania najmniejszych kwadratów oraz samego dopasowania jest niepewny. Czułość na warunki początkowe powoduje, że liniowy obszar dopasowania wykazuje oscylacyjne zachowanie podczas śledzenia dywergencji na atraktorze. Biologiczna zmienność między robakami, które były uśredniane dla każdej grupy wiekowej, przewyższa zmienność w liniowym dopasowaniu LLE dla jednego badania, co pokazano w poprzedniej publikacji10.
Seria czasowa (Rysunek 2A) pokazuje, że szczyty i doliny pokrywają się z aspektami poruszania się robaka, takimi jak zmiany częstotliwości, kształtu i orientacji pływania. Szeregi czasowe są aperiodyczne; Oscylacje w szeregu czasowym nigdy się nie powtarzają, a jednak pozostają ograniczone intensywnością. Częstotliwość jest stała, choć nieznacznie się waha, co wskazuje na złożoność i ograniczenie częstotliwości. Rysunek 10 pokazuje spadek częstotliwości pływania wraz z wiekiem robaka, który zanika się rozkłada 12 dnia. Średnia częstotliwość skacze po 12. dniu u większości robaków. To prawdopodobnie oznaka awarii obwodu nerwowego. Natomiast średnie opóźnienie τ, określone przez pierwsze minimum MI, rośnie wraz z wiekiem populacji, ponieważ cykl się wydłuża (Rysunek 11).
W naszym przypadku MI reprezentuje statystyczne nakładanie się dwóch opóźnionych szeregów czasowych. Jeśli MI jest minimalne, to nakładanie się trajektorii jest minimalne. Teoretycznie chaotyczne trajektorie nigdy się nie nakładają; jednak w praktyce, przy ograniczonych znaczących cyfrach i pewnym szumie eksperymentalnym, niektóre punkty będą się nakładać w granicach tolerancji24. Dążymy do minimalizacji możliwości nakładania się poprzez minimalizację MI. Różnica między zminimalizowanym MI a znacznie większym MI jest zilustrowana na Rysunku 4.
Wymiar osadzenia jest określany przez fałszywych najbliższych sąsiadów (FNN). Liczba FNN spłaszcza się o 5% lub mniej FNN wokół wymiarów osadzenia 3 lub 4. Na Rysunku 6 wiek nicieni nie wpływa na wymiar osadzenia, z wyjątkiem 12. dnia, gdy robak jest już prawie u końca swojego życia.
Trajektoria LLE podąża za krzywą Moore'a (Rysunek 7) z szczytem po 5 dniach, co wskazuje, że najbardziej nieprzewidywalna (i być może złożona) lokomocja występuje, gdy C. elegans są świeżo dojrzałe. Na Rysunku 7 wykres LLE od dni 3 do 12 pokazuje liniowy wzrost z szczytem po 5 dniach, a następnie spadek po 5 dniach. Słupki błędu na Rysunku 10 reprezentują odchylenie standardowe średniej 4,34 i odzwierciedlają kilka czynników, w tym różnorodność gatunków biologicznych oraz szacunki LLE. Różnice między robakami zwykle przeważają nad niepewnością w rutynie dopasowania LLE. Zmienność między dniami, jak dni 5 i 7, nie pokrywa się, więc dane są wyraźnie rozróżnialne. Trend LLE z wiekiem bardzo odpowiada równaniu Moore'a z wcześniejszego badania18, które opisuje wpływ wieku na pamięć i lokomocję u innych organizmów żywych.
LLE w Tabeli 1 są spójne, pokazując powolny wzrost, a następnie malejący trend wraz z wiekiem. Wartości są zgodne z wcześniej opublikowanymi wynikami na podstawie DOD10 oraz analizy wideo34,35.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/68955/68955fig1.jpg" /> 
Rysunek 1: Eksperymentalny układ dyfrakcji dalekiego pola (nie narysowany w skali). Lusterka sterujące tworzą peryskop. Kuweta z C. elegans umieszczana jest między lusterkami sterującymi. Wiązka laserowa jest dyfrakcyjna przez robak i przemieszcza się w kierunku PD przez drugie lustro sterujące. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig1large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/68955/68955fig2v2.jpg" /> 
Rysunek 2: Szeregi czasowe intensywności. Szeregi czasowe DOD wykazują wahania intensywności podczas poruszania nicienia w wiązce laserowej, z (A) odpowiednim szeregiem czasowym do analizy danych. (B) Ten szereg czasowy wskazuje, że odstęp między 11 a 30 s nie wykazuje sygnału; Pokazuje jedynie poziom szumów systemu, ponieważ krótsze amplitudy stałe pochodzą głównie z rozproszonego światła. (C) Ta seria czasowa zawiera kilka przypadków nasyconych danych. Każdy punkt danych stanowi znaczący skok na wykresie, który jest odcięty i spłaszczany na szczytach i dołkach. W tym przykładzie kolce spłaszczają się, gdy intensywność waha się od -100 do 100 AU w zakresie od 45 do 50 s. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig2large.jpg" target="_blank">Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/68955/68955fig3.jpg" /> 
Rysunek 3: Spektrum mocy. Widmo mocy pokazuje maksimum około 0,95 Hz. Częstotliwości są rozproszone, ponieważ częstotliwość trajektorii atraktora się przesuwa i nigdy dokładnie się nie powtarza. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig3large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/68955/68955fig4.jpg" /> 
Rysunek 4: Dane eksperymentalne osadzone w wykresach 3D z opóźnieniami z dwoma różnymi opóźnieniami tej samej serii czasowej 9-dniowego robaka. Ta seria czasowa pokazuje, że lokomocja jest ograniczonym atraktorem zgodnym z teorią chaosu. (A) Trajektoria jest rozwiązywana przy użyciu pierwszego lokalnego minimum (MI ≈ 2,11), co daje opóźnienie około 0,183 s (183 punkty danych). Widoczne przecięcia trajektorii są wynikiem rzutowania na przestrzeń dwuwymiarową. (B) Ten wykres opóźnień pozostaje nierozwiązany przez niewłaściwe opóźnienie 0,002 s (2 punkty danych, (MI > 7)), ponieważ punkty są zbyt blisko siebie, by je rozróżnić, i nie wykazują rozbieżności. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig4large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/68955/68955fig5.jpg" /> 
Rysunek 5: Wykładniczy spadek wzajemnej informacji. Pierwsze minimum w MI określa opóźnienie zbioru danych w 0,161 s (161 punktów danych) w tym konkretnym przypadku, aby rozstrzygnąć trajektorię fazową. To opóźnienie przesuwa wartość szeregu czasowego, aby odtworzyć atraktor w przestrzeni fazowej. Opóźnienie musi być tylko blisko pierwszego minimum, wystarczające, by rozstrzygnąć trajektorie. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig5large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/68955/68955fig6.jpg" /> 
Rysunek 6: Średnie wymiary osadzenia odpowiednio w dniach 3, 5, 9 i 12. Powiększony widok wymiaru osadzenia służy do pokazania konkretnej różnicy między wymiarami osadzenia każdego dnia. Tylko dzień 12 wykazuje zauważalną różnicę. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig6large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/68955/68955fig7.jpg" /> 
Rysunek 7: Średnia szacowanych LLE między dniem 3 a 12. LLE osiąga szczyt po 5 dniach. Przez maksymalnie 5 dni nicienie wykazują trajektorię wzrostu zgodną z przejściem od niedojrzałości do dojrzałości. Po tych 5 dniach nicienie doświadczają spadku rozbieżności. Niepewności są miarą odchylenia standardowego średniej, które jest zdominowane przez zmiany między robakami. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig7large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/68955/68955fig8.jpg" /> 
Rysunek 8: Przedstawienie rozbieżności obecnej w teorii chaosu. Dodatnia LLE spowoduje, że początkowo bliskie trajektorie będą się z czasem rozchodzić. t to przedział czasowy, d to dywergencja, a x(t) to punkt w przestrzeni fazowej. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig8large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/68955/68955fig9.jpg" /> 
Rysunek 9: Dywergencja trajektorii fazowej w skali logarytmicznej. Szacowany LLE to nachylenie liniowego dopasowania (1,08 1/s) w obszarze wznoszenia przed wypłaszczeniem się z powodu ograniczenia atraktora między 0,96 a 1,01 s. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig9large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tego wykresu.</a>
<img alt="Figure 10" class="xfigimg" src="/files/ftp_upload/68955/68955fig10.jpg" /> 
Rysunek 10: Średnia częstotliwość występowania każdego dnia. Średnia z dnia 12 oznacza istotną zmianę sygnałów neurologicznych napędzających lokomocję. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig10large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<img alt="Figure 11" class="xfigimg" src="/files/ftp_upload/68955/68955fig11.jpg" /> 
Rysunek 11: Średnie opóźnienie dla każdego dnia. MI prowadzi średnio do niższych opóźnień między dniem 3 a 6, podczas gdy wyższe średnie opóźnienie pojawia się między dniem 7 a 12. <a href="https://www.jove.com/files/ftp_upload/68955/68955fig11large.jpg" target="_blank">Proszę kliknąć tutaj, aby zobaczyć większą wersję tej figurki.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Wiek (dni)</td>
			<td>Liczba zbiorów danych</td>
			<td>LLE (1/s)</td>
			<td>Odchylenie standardowe średniej (1/s)</td>
		</tr>
		<tr>
			<td>3</td>
			<td>11</td>
			<td>1.14</td>
			<td>0.02</td>
		</tr>
		<tr>
			<td>4</td>
			<td>15</td>
			<td>1.20</td>
			<td>0.03</td>
		</tr>
		<tr>
			<td>5</td>
			<td>13</td>
			<td>1.34</td>
			<td>0.03</td>
		</tr>
		<tr>
			<td>6</td>
			<td>11</td>
			<td>1.23</td>
			<td>0.02</td>
		</tr>
		<tr>
			<td>7</td>
			<td>10</td>
			<td>1.16</td>
			<td>0.03</td>
		</tr>
		<tr>
			<td>8</td>
			<td>12</td>
			<td>1.13</td>
			<td>0.02</td>
		</tr>
		<tr>
			<td>9</td>
			<td>12</td>
			<td>1.07</td>
			<td>0.02</td>
		</tr>
		<tr>
			<td>10</td>
			<td>9</td>
			<td>1.00</td>
			<td>0.03</td>
		</tr>
		<tr>
			<td>11</td>
			<td>9</td>
			<td>0.96</td>
			<td>0.03</td>
		</tr>
		<tr>
			<td>12</td>
			<td>9</td>
			<td>0.92</td>
			<td>0.02</td>
		</tr>
	</tbody>
</table>
Tabela 1: Liczba zbiorów danych i średnia LLE dla każdego dnia badania wiekowego. LLE są zgodne z wcześniejszymi danymi zebranymi przez DOD i pokazują rosnącą, a następnie malejącą trajektorię krzywej Moore'a na Rysunku 7. Zebrano 9-15 zbiorów danych na każdy dzień danych. Zbiory danych spadały wraz z wiekiem robaków ze względu na charakter degradacji mięśni i spadku nerwów, co powodowało zmniejszoną lokomocję.

Badanie to analizuje wpływ wieku na lokomocję u C. elegans poprzez pomiar największego wykładnika Lyapunova (LLE). Wraz z wiekiem C. elegans wykazują wzrost i późniejszy spadek kontroli motorycznej. Wyniki pokazują szczyt LLE po pięciu dniach, a następnie spadek wraz z wiekiem robaków.

Wiekowa dynamika lokomocji u Caenorhabditis elegans: analiza wykładnika Lyapunova

This study investigates the influence of age on the mobility of Caenorhabditis elegans (C. elegans) by employing Dynamic Optical Diffraction (DOD) to ...

W ramach tych badań dążymy do scharakteryzowania lokomocji mikroskopijnego gatunku, aby lepiej zrozumieć układ nerwowy tego ruchu. Dokładniej ustaliłyśmy, że ruch C. elegans jest chaotyczny i doświadczalny jest niski szum eksperymentalny. Eksperymentalny chaos może być trudny do analizy w kontekście eksperymentalnego szumu.Dynamiczna dyfrakcja optyczna, czyli DOD, jest stosunkowo nowa. Nie ma odpowiedników w tej dziedzinie, jest łatwy do wdrożenia, opłacalny i dostarcza dane wieloskalowe. Wieloskalowa informacja dynamiczna sprawdza modelowany obwód neuronowy.Gdy zbudujemy układ nerwowy z dopasowanymi informacjami dynamicznymi, wiemy, że go zrozumieliśmy. W przyszłości nasze badania będą koncentrować się na eksploracji, ale także innych wielkości biologicznych, które możemy scharakteryzować największym wykładnikiem Ljapunowa. Dodatkowo zbudujemy model neuronów lokomotorycznych.Na początek pipetuj 0,5 mililitra Escherichia coli z OD600 na każdą płytkę agarową z ośrodkiem wzrostu nicieni, aby zapewnić pożywienie nicieniom. Poczekaj, aż Escherichia coli całkowicie wyschnie na powierzchni każdej płytki agarowej. Pobierz płytkę kontrolną Caenorhabditis elegans od dostawcy materiałów biologicznych, aby wykorzystać ją do generowania populacji nicieni zsynchronizowanych z wiekiem.Używając palnika Bunsena, wysterylizuj platynowy kostkę w płomieniu, aż zacznie świecić pomarańczowo. Pod mikroskopem rozcinającym użyj sterylnego kilofa do zebrania od pięciu do dziesięciu dorosłych dzikich Caenorhabditis elegans. Przełóż zebrane robaki na przygotowaną płytkę agarową.Po okresie składania jaj trwającym cztery do pięciu godzin usuń dorosłe robaki z talerza za pomocą platynowego kilku. Inkubuj pozostałą płytkę, aby jaja mogły się rozwijać do pożądanego dnia użycia. W dniu zbierania danych napełnij kwarcową kuwetę optycznej o średnicy 4,5 mililitra wodą destylowaną w temperaturze pokojowej tuż poniżej krawędzi, aby zapobiec rozlaniu się.Użyj platynowego kilofa do delikatnego przenieśienia dwóch do trzech nicieni. Gdy nicienie są już w środku, połóż kuwetę na boku, aby ułatwić wyrównanie robaka z wiązką lasera. Obserwuj ruch pływania w pełni zanurzonych nicieni, które mogą unosić się ku dnu, ale pozostają aktywne.Uruchomił system eksperymentalny z modyfikacją polegającą na konstrukcji peryskopu z użyciem dwóch luster. Umieść kuwetę na boku między dwoma lustrami, aby ułatwić ustawienie robaka w wiązce laserowej. Zamiast aparatu, ustaw detektor zdjęć w wzorze dyfrakcyjnym, aby zbierać sygnały optyczne.Następnie włącz neon z helem i pozwól mu się ogrzać przez około 15 minut, aż osiągnie równowagę termiczną. Uruchom cyfrowy oscyloskop, aby rozpocząć zbieranie danych. Aby skonfigurować oscyloskop, ustaw interwał czasowy na 100 sekund i wybierz bufor pamięci o wydajności co najmniej 100 kilo próbek na sekundę.Ustaw rozdzielczość akwizycji danych na jeden kiloherc. Dostosuj rozdzielczość bitów do 12 bitów, aby zapewnić czułość na małe zmiany amplitudy. Użyj funkcji automatycznego przesunięcia AC na oscyloskopie, aby wycentrować oscylacje intensywności na zerze woltów.Następnie użyj sterylnego platynowego kilofa delikatnie zerwij dwa do trzech nicieni i przełóż je do kuwety wypełnionej wodą destylowaną. Włóż kuwetę zawierającą nicienie do zespołu peryskopowego i dostosuj ustawienie tak, aby jeden nicienie był dokładnie wycentrowany na ścieżce wiązki laserowej. Obserwuj powstawanie wzoru dyfrakcyjnego w dalekim polu około 50 centymetrów za peryskopem, gdy Caenorhabditis elegans zostanie wycentrowany w wiązce lasera.Umieść fotodiodę w wzorze dyfrakcyjnym dalekiego pola, gdy nicienie przechodzi przez wiązkę lasera. Zapisuj dane przez co najmniej 10 sekund, aby zebrać co najmniej 10 000 punktów danych, co jest niezbędne do dokładnego obliczenia największego wykładnika Liapunowa. Przeanalizuj szeregi czasowe, aby zidentyfikować segmenty, w których nicienie swobodnie pływa w promieniu lasera.Przeanalizuj szeregi czasowe, aby ocenić stosunek sygnału do szumu. Ustalić poziom szumu bazowego, nagrywając szereg czasowy bez nicieni w kuwecie. Potwierdź, że zarówno amplituda, jak i cechy sygnału eksperymentalnego przekraczają dwukrotność poziomu szumu.Usuń dowolną część szeregu czasowego zawierającą nasycone dane. Zidentyfikuj pierwsze minimum w wykresie informacji wzajemnej w przedziale czasowym od 0,140 do 0,240 sekundy. Wybierz wartość opóźnienia, która skutecznie oddziela trajektorie fazowe, aby umożliwić wyraźne pomiary dywergencji.Określ wymiar osadzenia za pomocą metody fałszywych najbliższych sąsiadów. Wybierz najniższy wymiar, przy którym procent fałszywych najbliższych sąsiadów stabilizuje się na minimum. Największy wykładnik Ljapunowa osiągnął szczyt piątego dnia, a następnie stopniowo spadał od szóstego do dwunastego dnia.Rozbieżność pobliskich trajektorii początkowo rosła wykładniczo, zanim osiągnęła płaskowyż, potwierdzając ograniczone, chaotyczne zachowanie układu. Na początku krzywej logarytmicznej zaobserwowano nachylenie liniowe, co daje wiarygodne oszacowanie największego wykładnika Ljapunowa. Średnia częstotliwość pływania spadła od trzeciego do dwunastego dnia, z zauważalnym załamaniem trendu w dwunastym dniu.

age-dependent-dynamics-locomotion-caenorhabditis-elegans-lyapunov

Research

JoVE Journal

Bioengineering

Bioinżynieria

Opracowaliśmy system mikroskopu śledzącego w czasie rzeczywistym, który umożliwia rejestrację i kwantyfikację zachowania C. elegans z wysoką rozdzielczością i przy dużych prędkościach. Opracowaliśmy również metody obliczeniowe pozwalające na redukcję wymiarowości obrazów nicieni do podstawowego zestawu pomiarów, które w pełni opisują kształt organizmu.

Śledzenie i pomiar zachowania C. elegans

Badacze mogą uruchomić ruch na żądanie u C. elegans, stosując pole elektryczne, które robaki wyczuwają i na które reagują. Ten film przedstawia zachowanie elektrotaktyczne i pokazuje przykładowy protokół, który jest wykonywany w urządzeniu mikroprzepływowym.

Test elektrotaksji: metoda obserwacji lokomocji u C. elegans

Ten film opisuje użycie oprogramowania do nagrywania ruchu robaków na arenie. Przykładowy protokół mierzy wpływ ekspozycji na etanol na poruszanie się robaków.

C. elegans (C. elegans) Śledzenie ruchu: metoda oceny lokomocji u robaków

Ten film demonstruje zastosowanie testu lokomocji promieniowej do oceny dysfunkcji motorycznej u C. elegans wykazujących ekspresję ludzkiego białka TDP-43 typu dzikiego lub z mutacją ALS. Robaki o różnym poziomie ekspresji TDP-43 typu dzikiego lub zmutowanego mogą pełzać, a ich ruch jest mierzony odległością promieniową przebytą od centralnego punktu na płytce testowej. Test rozróżnia łagodne, umiarkowane i ciężkie upośledzenie motoryczne na podstawie odległości przebytej przez różne szczepy C.

Ocena upośledzenia motorycznego w modelach C. elegans za pomocą testu lokomocji radialnej

Opisano półautomatyczną metodę mikro-elektro-fluidyczną do wywoływania lokomocji na żądanie u Caenorhabditis elegans. Metoda ta opiera się na neurofizjologicznym zjawisku robaków reagujących na łagodne pola elektryczne ("elektrotaksję") wewnątrz kanałów mikroprzepływowych. Elektrotaksja mikroprzepływowa służy jako szybka, czuła, tania i skalowalna technika badań przesiewowych pod kątem czynników wpływających na zdrowie neuronów.

Elektrotaksja oparta na mikroprzepływach do analizy ilościowej na żądanie lokomocji Caenorhabditis elegans

Przedstawiona tutaj technika mierzy ścieżkę swobodnie pływających mikroskopijnych gatunków za pomocą ekspozycji na pojedynczej długości fali.  C. elegans są używane do demonstracji obrazowania cieni jako niedrogiej alternatywy dla kosztownych mikroskopów.  Technikę tę można dostosować do różnych orientacji, środowisk i gatunków w celu pomiaru kierunku, prędkości, przyspieszenia i sił.

Obrazowanie cieni o pojedynczej długości fali lokomocji Caenorhabditis elegans, w tym oszacowania siły

Opisana jest skuteczna i prosta metodologia komputerowej analizy zachowania nicieni podczas pływania w cieczy. Metoda ta wymaga niewielkich lub żadnych inwestycji dla laboratoriów C. elegans. Używany sprzęt jest standardowy, a oprogramowanie komputerowe do analizy behawioralnej (CeleST) jest oprogramowaniem typu open source.

Zautomatyzowana analiza zachowania C. elegans podczas pływania przy użyciu oprogramowania CeleST

Ten manuskrypt opisuje, jak rozróżniać różne nicienie za pomocą sygnatur dyfrakcyjnych dalekiego pola. Porównujemy lokomocję 139 dzikich i 108 "rolkowych" C. elegans, uśredniając częstotliwości związane z czasową sygnaturą dyfrakcyjną Fraunhofera w jednym miejscu za pomocą lasera o fali ciągłej.

Analiza dyfrakcyjna oparta na Fourierze żywych Caenorhabditis elegans

W tym badaniu demonstrujemy wykorzystanie kinematycznej analizy chodu opartej na obrazowaniu w płaszczyźnie brzusznej do monitorowania subtelnych zmian w koordynacji ruchowej, jak również postępu neurodegeneracji wraz z wiekiem w modelach mysich (np. linie myszy z mutacją endofiliny).

Analiza chodu zależnych od wieku zaburzeń motorycznych u myszy z neurodegeneracją

To badanie przedstawia protokoły dla dwóch półautomatycznych metod analizy aktywności lokomotorycznej u robaków C. elegans complex I disease gas-1(fc21), a mianowicie ZebraLab (test o średniej przepustowości) i WormScan (test o wysokiej przepustowości) i dostarcza analizy porównawczej wśród szerokiego wachlarza metod badawczych do ilościowego określenia zachowania nicieni i zintegrowanej funkcji nerwowo-mięśniowej.

Analiza porównawcza metod eksperymentalnych do ilościowego określania aktywności zwierząt w modelach choroby mitochondrialnej Caenorhabditis elegans

Ten protokół opisuje dwa czułe testy do rozróżniania między łagodnym, umiarkowanym i ciężkim upośledzeniem motorycznym w modelach stwardnienia zanikowego bocznego C. elegans, z ogólną użytecznością dla szczepów C. elegans, ze zmienioną ruchliwością.

Ocena upośledzenia motorycznego w modelach stwardnienia zanikowego bocznego C. elegans

Odkryj analizę wiekowej dynamiki lokomocji u Caenorhabditis elegans za pomocą wykładnika Lyapunova. Metoda ta pozwala badać zmiany w kontroli ruchowej związane ze starzeniem się organizmu.

Wiekowa dynamika lokomocji u Caenorhabditis elegans: analiza wykładnika Lyapunova

Przeglądaj więcej filmów