Analiza dyfrakcyjna oparta na Fourierze żywych Caenorhabditis elegans

Podsumowanie

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.

Streszczenie

Ten manuskrypt opisuje, jak klasyfikować nicienie za pomocą sygnatury dyfrakcyjnej czasowego dalekiego pola. Pojedynczy C. elegans jest zawieszony w słupie wody wewnątrz kuwety optycznej. Laser HeNe o długości fali ciągłej o długości fali 632 nm jest kierowany przez kuwetę za pomocą lusterek na przedniej powierzchni. Znaczna odległość co najmniej 20-30 cm przebyta po przejściu światła przez kuwetę zapewnia użyteczny wzór dyfrakcyjny dalekiego pola (Fraunhofer). Wzór dyfrakcyjny zmienia się w czasie rzeczywistym, gdy nicienie pływają w wiązce laserowej. Fotodioda jest umieszczona poza środkiem we wzorze dyfrakcyjnym. Sygnał napięciowy z fotodiody jest obserwowany w czasie rzeczywistym i rejestrowany za pomocą oscyloskopu cyfrowego. Proces ten powtarza się dla 139 dzikich i 108 "walcowych" C. elegans. Robaki typu dzikiego wykazują szybki wzorzec oscylacji w roztworze. Robaki "rolkowe" mają mutację w kluczowym elemencie naskórka, który zakłóca płynną lokomocję. Interwały czasu, które nie są wolne od nasycenia i bezczynności, są odrzucane. Praktyczne jest podzielenie każdej średniej przez jej maksimum w celu porównania względnych intensywności. Sygnał dla każdego robaka jest przekształcany w Fouriera, tak aby wyłaniał się wzorzec częstotliwości dla każdego robaka. Sygnał dla każdego typu robaka jest uśredniony. Uśrednione widma Fouriera dla typu dzikiego i "rolkowego" C. elegans są wyraźnie różne i ujawniają, że dynamiczne kształty robaków dwóch różnych szczepów robaków można rozróżnić za pomocą analizy Fouriera. Widma Fouriera każdego szczepu robaka pasują do przybliżonego modelu wykorzystującego dwa różne binarne kształty robaków, które odpowiadają momentom lokomotorycznym. Obwiednia uśrednionego rozkładu częstotliwości dla rzeczywistych i modelowanych robaków potwierdza, że model jest zgodny z danymi. Metoda ta może służyć jako punkt odniesienia dla analizy Fouriera dla wielu mikroskopijnych gatunków, ponieważ każdy mikroorganizm będzie miał swoje unikalne widmo Fouriera.

Wprowadzenie

Ta metoda porównuje eksperymentalne i modelowane widma częstotliwości lokomocji C. elegans przy użyciu dwóch szczepów o bardzo różnych wzorcach lokomotorycznych. Wyniki pokazują, że spektrum częstotliwości zależy od zmian czasowych, gdy nicienie pływają w słupie wody, więc do analizy nie są potrzebne wyraźne obrazy mikroskopowe. Metoda ta pozwala na ilościową analizę w czasie rzeczywistym i dostarcza informacji uzupełniających do obrazów/filmów uzyskanych za pomocą tradycyjnych mikroskopów. Dyfrakcja Fraunhofera, zwana również dyfrakcją dalekiego pola, stanowi podstawę do uzyskiwania danych dyfrakcji na żywo^1,2. Natężenie światła w dowolnym pojedynczym punkcie wzoru dyfrakcyjnego jest wynikiem nałożenia światła z każdego punktu w obrysie nicienia³. W rezultacie natężenie światła zebrane w czasie niesie ze sobą informację o lokomocji nicienia. Analiza zależnego od czasu sygnału dyfrakcyjnego może zidentyfikować charakterystyczny ruch odpowiedniego mutanta, ponieważ analiza wszystkich częstotliwości związanych z lokomocją, uzupełnia tradycyjną analizę wideo. W tym przypadku charakterystyczne różnice między lokomocją, jaką uprawia się "walcowo" a dzikim typem C. elegans, potwierdza się poprzez porównanie widm częstotliwości dwóch różnych szczepów nicieni.

Niektóre poprzednie cechy zostały potwierdzone za pomocą analizy częstotliwości sygnałów dyfrakcyjnych, takich jak częstotliwości pływania^2,4. Co ważniejsze, metoda ta może być stosowana jako metoda uzupełniająca tradycyjną mikroskopię do obserwacji lokomocji w czasie rzeczywistym na ekranie komputera podczas zbierania danych. Spektrum częstotliwości robaków z wyraźnymi wzorcami lokomotorycznymi można określić ilościowo, biorąc pod uwagę sygnał transformowany przez Fouriera sygnału dyfrakcyjnego.

Multidyscyplinarny charakter dyfrakcji opartej na Fourierze w tej pracy obejmuje dziedziny biologii i fizyki. Dyfrakcja przez niedostateczne próbkowanie jest od dawna stosowana do badania struktur krystalicznych w biologii^{5 i innych dziedzinach}. Jednak w tym eksperymencie nadpróbkowanie^6,7 tworzy wzór dyfrakcji dalekiego pola, tak że organizm jest wyśrodkowany w wiązce laserowej. Nadpróbkowanie jest zwykle używane do obrazowania bez soczewki⁸ w połączeniu z algorytmem odzyskiwania fazy, który rekonstruuje obraz oryginalnego obiektu. Odzyskiwanie faz jest trudne do osiągnięcia, gdy obecne są rozpraszacze, jak ma to miejsce w przypadku nicieni. Sygnatura dyfrakcji czasowej jest wystarczająca do oceny kluczowych częstotliwości ruchu robaka. Ta metoda jest mniej obciążająca obliczeniowo i zapewnia optyczny sposób ilościowego określania lokomocji. Technika ta może być łatwo zaadaptowana do analizy mutacji lub warunków środowiskowych, które zmieniają zachowanie.

Protokół

1. Wzrost i utrzymanie C. elegans

1. Przygotuj potrawy z kultur nicieni.
   1. Napełnij szalki Petriego roztworem agaru i pozostaw je do zestalenia, a następnie wysiewaj je kulturą E. coli szczepu OP50 ^9,10.
2. Przygotuj początkową populację dorosłych nicieni na każdej tablicy, przenosząc kilka dorosłych robaków na świeże szalki Petriego wypełnione agarem z plastrem E. coli. Utrzymywać kultury nicieni w temperaturze 20 °C w inkubatorze.
   UWAGA: Szczepy nicieni można uzyskać w Centrum Genomu Caenorhabditis elegans. W tym badaniu wykorzystano szczep typu dzikiego, N2, i szczep OH7547 (otls199[cat-2::GFP + rgel-1(F25B3.3)::d sRed + rol-6(su1006)]), który wykazuje fenotyp rolkowy.
3. Rozmnażaj robaki dla przyszłych kultur.
   1. Wyjąć szalkę Petriego zawierającą C. elegans i nieużywaną szalkę Petriego z inkubatora o kontrolowanej temperaturze. Umieść je na stoliku mikroskopu preparacyjnego.
   2. Zapal palnik Bunsena i wysterylizuj platynowy nicień, umieszczając metal w płomieniu, aż zaświeci się na czerwono. Pozwól kilofowi ostygnąć do temperatury pokojowej. Nie odkładaj kilofa ani nie pozwól, aby miał on kontakt z zanieczyszczeniami.
   3. Delikatnie dotknij końcówką kilofa krawędzi kręgu bakterii. Substancja ta jest lepka i ułatwi zbieranie pojedynczych dorosłych nicieni.
   4. Przenieść do 4 grawitacyjnych nicieni na pożywkę Petriego wypełnioną agarem z pożywką do wzrostu nicieni (NGM) i inkubować w temperaturze 20 °C. Robaki złożą jaja, które dojrzeją w ciągu czterech dni.
   5. Pozostałe nicienie należy umieścić z powrotem w inkubatorze po przeniesieniu czterech dorosłych nicieni.

2. Układ optyczny (Rysunek 1)

1. Zamocuj laser helowo-neonowy w pobliżu lewego tylnego rogu stołu optycznego i podłącz go do źródła zasilania.
   UWAGA: Wiązka lasera musi spełniać wymagania dotyczące nadpróbkowania. C. elegans mają około 1 mm długości, dlatego wiązka laserowa powinna mieć średnicę większą niż 2 mm, gdy pada na nicienie, ale nie większą niż 5 mm, aby wzór dyfrakcyjny nie był trudny do zlokalizowania.
2. Umieść filtr o neutralnej gęstości między laserem helowo-neonowym a próbką w taki sposób, aby wiązka lasera przechodziła przez filtr przed dotarciem do próbki.
3. Korzystając z dwóch aluminiowych lusterek kierowniczych na przedniej powierzchni, zbuduj peryskop, zabezpieczając pierwsze lusterko za filtrem o neutralnej gęstości. Zamocuj drugie lustro około 10 cm poniżej pierwszego lustra, aby dać miejsce na sterowanie wiązką laserową i włożenie kuwety między lustra. Ustaw wiązkę lasera i kuwetę tak, aby wiązka lasera przechodziła pionowo przez kuwetę.
   UWAGA: Odległość od organizmu dyfrakcyjnego do fotodiody musi być znacznie większa niż sam organizm, aby uzyskać dyfrakcję w polu dalekim. W tym eksperymencie odległość od kuwety do fotodiody wynosi 20 cm.
4. Zamocuj fotodiodę bezpośrednio naprzeciwko drugiego lustra, tak aby jej czujnik był skierowany w stronę lustra.
   UWAGA: Kuweta zawierająca nicienie zostanie umieszczona między dwoma lustrami za pomocą zacisków chemicznych. Patrz punkty 4 i 5.
5. Umieść kuwetę wypełnioną wodą na stojaku. Dostosuj wysokość podstawy. Dostosuj wysokości i kąty lustra 1 i lustra 2 tak, aby wiązka laserowa przechodziła przez kuwetę skierowaną w pobliżu, ale nie bezpośrednio na fotodiodę.
6. Użyj poziomicy, aby upewnić się, że stojaki tworzą równą powierzchnię dla kuwety. W razie potrzeby wyreguluj lusterka dalej.
7. Podłącz fotodiodę do oscyloskopu cyfrowego za pomocą USB dostarczonego z oscyloskopem cyfrowym. Podłącz oscyloskop cyfrowy do komputera, który będzie używany do nagrywania i zapisywania danych.

3. Ustawienia oscyloskopu

1. Korzystając z oprogramowania oscyloskopu w komputerze, ustaw częstotliwość próbkowania na co najmniej 8 Hz, aby wystarczająco rozwiązać cykl młócenia robaka.
   UWAGA: Częstotliwość próbkowania powinna być ponad dwukrotnie większa niż oczekiwana częstotliwość młócenia gatunku, aby spełnione było twierdzenie Nyquista¹¹.

4. Przygotowanie robaka i kuwety do zbierania danych

1. Przenieś cztery dorosłe nicienie na świeżą płytkę Petriego wypełnioną agarem NGM¹⁰ za pomocą cienkiego, spłaszczonego szpikulca z drutu platynowego (patrz sekcja 1.3).
2. Wyjmij jedną jednorazową plastikową kuwetę z opakowania, uważając, aby dotykać kuwety tylko jej prążkowanymi bokami.
3. Za pomocą mikropipety odpipetuj wodę destylowaną do kuwety, aż kuweta zostanie wypełniona w około 80% wodą destylowaną.
   UWAGA: Ważne jest, aby podczas pracy z nicieniami używać wyłącznie wody destylowanej lub zjonizowanych, takich jak M9¹⁰ lub sól fizjologiczna buforowana fosforanami (PBS), ponieważ woda z kranu zawiera związki zabijające mikroorganizmy.
4. Umieścić szalkę Petriego zawierającą C. elegans, która ma być użyta, pod lunetą preparacyjną.
5. Używając platynowego kilofa, wyjmij jeden dojrzały C. elegans z szalki Petriego i zanurz go w kuwecie, przesuwając kilkor, jeśli to konieczne, aby usunąć nicienie.
6. Aby zapobiec tworzeniu się pęcherzyków w kuwecie, napełnij kuwetę wodą, aż lekko wybrzuszy się nad górną częścią kuwety. Napełnij całkowicie nasadkę kuwety wodą, a następnie szybko załóż nasadkę na kuwetę.
7. Użyj szmatki do czyszczenia optycznego, aby usunąć kropelki wody, które mogły się rozlać, i papieru do czyszczenia optycznego, aby wyczyścić wszelkie małe pozostałe kropelki.

5. Akwizycja danych w czasie rzeczywistym o zmianach intensywności wzoru dyfrakcyjnego

1. Włącz laser helowo-neonowy i dostosuj ustawienie częstotliwości/koloru tak, aby wytwarzał czerwoną wiązkę. Włącz czujnik.
   ostrożność: Użyj wiązki o niskiej energii przy długości fali 632 nm, ponieważ C. elegans unika światła o wysokiej częstotliwości (niebieskiego)¹².
2. Zlokalizuj robaka w kuwecie. Trzymając kuwetę po jej prążkowanych bokach, delikatnie przechylaj kuwetę, aż nicienie znajdą się mniej więcej w środku części kuwety.
   UWAGA: Potrząsanie lub przechylanie kuwety powoduje gwałtowne zderzenie robaka ze ściankami kuwety. Może to uszkodzić nicienie.
   1. Umieść kuwetę na stojaku w układzie optycznym, centrując robaka w wiązce lasera odbitej od lustra 1 do lustra 2. Pamiętaj, aby wyeliminować rozproszone światło.
3. Wyśrodkuj robaka w wiązce laserowej.
   1. Umieść fotodiodę we wzorze dyfrakcyjnym tak, aby położenie fotodiody i centralne maksimum wzoru dyfrakcyjnego nie pokrywały się.
   2. Dostosuj filtr o neutralnej gęstości, aby zapobiec nasyceniu fotodiody. Obracanie koła filtrów o neutralnej gęstości reguluje natężenie światła.
      1. Obróć filtr o neutralnej gęstości tak, aby napięcie wyjściowe z fotodiody wzrosło.
         UWAGA: Napięcie wyjściowe jest obserwowane za pomocą oprogramowania fotodiody cyfrowej. Fotodioda jest nasycona, jeśli napięcie się nie zmienia. W takim przypadku obracaj filtr o neutralnej gęstości, aż napięcie spadnie bez spłaszczania się przy minimalnym odczycie. Upewnij się, że voltage sygnał nie spłaszcza się przy szczytowych odczytach, wskazując na nasycenie fotodiody. Zmniejsz natężenie światła, obracając koło filtrów o neutralnej gęstości, jeśli zaobserwowane jest nasycenie.
   3. Gdy poruszający się wzór dyfrakcyjny jest widoczny, zbierz dane za pomocą fotodiody, klikając przycisk start w oprogramowaniu sterującym oscyloskopem, jednocześnie monitorując ruch robaka. Kontynuuj wykonywanie pomiarów, aż robak wyjdzie z wiązki laserowej, a wzór dyfrakcyjny zniknie, co zwykle trwa około 20 sekund.
      1. Zatrzymaj proces zbierania danych, klikając przycisk zatrzymania w oprogramowaniu oscyloskopu. Zapisz dane każdej wersji próbnej w formacie .csv lub .txt.
4. Powtarzaj kroki 5.2 - 5.3, aż zostanie zebranych co najmniej 50 zestawów danych dla każdego fenotypu. Na próbę należy użyć od ośmiu do dziesięciu zwierząt.
5. Jeśli kuweta jest porysowana, wyrzuć ją i robaka, a następnie powtórz krok 4. Jeśli robak zostanie uszkodzony podczas przenoszenia, wyrzuć go i opłucz kuwetę wodą destylowaną przed powtórzeniem kroku 4 przy użyciu tej samej kuwety.
6. Powtórz krok 5 używając szczepu OH7547 "Roller", uważając, aby oznaczyć dane wskazujące szczep robaka.

6. Spektrum danych Fouriera

1. Zaimportuj uzyskane dane do programu do analizy danych, który jest w stanie wykonać dyskretne transformaty Fouriera³.
2. Wykonaj transformacje Fouriera na każdym zestawie danych, korzystając z opcji szybkiej transformaty Fouriera (FFT) oprogramowania.
3. Uśrednić częstotliwości z wyników FFT dla każdej amplitudy dla robaków typu dzikiego N2.
4. Powtórz krok 6.3 przy użyciu FFT z robaków "rolkowych" OH7547.

7. Modelowanie widma Fouriera

1. Zaprogramuj binarny model lokomocji nicieni (patrz program zawarty w materiałach dodatkowych).
   UWAGA: Ten model jest przybliżonym przybliżeniem, które na pierwszy rzut oka pokazuje wyraźne cechy ruchu robaków. Model można udoskonalić, porównując wyniki z rzeczywistymi robakami. Kształty robaków są przybliżeniami przy użyciu obrazów mikroskopowych¹³.
   1. Twórz sekwencyjne klatki modelu binarnego przechodzącego przez co najmniej dwa cykle robaka (Rysunek 2a). Zobacz filmy w materiałach dodatkowych; Film o robaku C (CWorm.avi) i film o robaku W (WWorm.avi).
2. Twórz sekwencyjne wzory dyfrakcyjne. Zobacz filmy w materiałach dodatkowych. Film z dyfrakcją robaka C (CWormDiff.avi) i film z robakiem W (WWormDiff.avi).
   1. Transformacja Fouriera polega na przekształceniu każdej binarnej klatki obrazu robaka. Im większa dopełnienie ramki wokół robaka, tym lepsza będzie rozdzielczość obrazu dyfrakcyjnego.
      UWAGA: Wartość bezwzględna każdej ramki przekształconej Fouriera jest proporcjonalna do odpowiadającego jej wzoru dyfrakcyjnego (Rysunek 2b).
   2. Dostosuj kontrast wzoru dyfrakcyjnego, odwzorowując intensywność wzoru dyfrakcyjnego na skalę logarytmiczną.
      UWAGA: Kamery i oczy mają tendencję do działania w skali nieliniowej. Skala logarytmiczna może symulować, w jaki sposób wzór dyfrakcyjny jest zwykle postrzegany przez ludzkie oko.
3. Wyodrębnij modelowany sygnał dyfrakcyjny.
   1. Wybierz położenie poza środkiem, które odpowiada położeniu fotodiody we wzorze dyfrakcyjnym.
   2. Dodaj sąsiednie elementy matrycy otaczające położenie fotodiody, aby zasymulować rozmiar fotodiody. Rozmiar fotodiody wynosi zazwyczaj 0,1% modelowanego wzoru dyfrakcyjnego.
   3. Zapisać i wykreślić sekwencję wielkości sygnałów dyfrakcyjnych. Sprawdź, czy sygnał jest fizycznie uzasadniony (tzn. w przypadku okresowego młócenia sygnał z fotodiody również powinien być okresowy).
4. Fourier przekształca sygnał dyfrakcyjny uzyskany w 7.3 i porównuje wyniki z danymi doświadczalnymi.
5. Powtórz dla różnych szczepów robaków i porównaj.

Wyniki

Eksperymentalny układ optyczny pokazany w Rysunek 1 pozwala na badanie mikroorganizmów bez konieczności przywiązywania się do płaszczyzny ogniskowej. Sygnał uderzeniowy z fotodiody można obserwować w czasie rzeczywistym na ekranie komputera w miarę zbierania danych. Nietypowe wzory będą widoczne od razu, bez konieczności szczegółowej analizy wideo.

Przykłady modelowanych sekwencyjnych ruchów robaków i odpowiadających im wzorów dyfrakcyjnych są pokazane w Rysunek 2. Modelowane wzory dyfrakcyjne jakościowo przypominają wzorce eksperymentalne¹ i są wstępną wskazówką, że symulacje z powodzeniem modelują nicienie.

Przykładowa sygnatura dyfrakcji czasowej dwóch badanych tutaj typów C. elegans jest pokazana w Rysunek 3. Można zauważyć jakościowo, że każdy nicień miota się w różnym tempie i amplitudzie. Niektóre różnice można określić ilościowo poprzez dopasowanie krzywej, tak jak to miało miejsce w poprzedniej publikacji¹. Dyskretna transformata Fouriera ujawnia jednak więcej szczegółów dotyczących wbudowanych częstotliwości:

, (1)

gdzie F_k to cyfrowa transformata Fouriera (FT), a f_n to zależny od czasu surowy sygnał dyfrakcyjny z dyskretną zmienną czasową n i dyskretną zmienną częstotliwości k. N to całkowita liczba punktów danych. Średnia cyfrowa transformata Fouriera pozwala na identyfikację nicienia na podstawie amplitudy jego widma częstotliwości dyfrakcyjnej (Rysunek 4). Widmo typu dzikiego jest zdominowane przez niższe częstotliwości niż widmo ruchu rolkowego.

Model, który przybliża typ dziki w porównaniu z walcem C. elegans zauważa, że typ dziki ma tendencję do miotania się ruchem przypominającym falę (kształt W lub S) (Rysunek 2a), podczas gdy walec ma tendencję do faworyzowania jednej strony, która z grubsza przypomina oscylujący kształt litery C (Rysunek 5). Daje to pewne wyjaśnienie dla różnych widm. Wałek w większości utworzy C po jednej stronie, podczas gdy oscylację W można traktować jako dwa przeciwstawne ruchy C. Z tego powodu ruch W jest bardziej złożony, ujawniając więcej wtórnych niskich częstotliwości niż ruch C. Wynik ten jest potwierdzony w modelu obliczeniowym. Kształt W ma znacznie większą gęstość częstotliwości niż kształt C (Rysunek 6). Jest to potwierdzone w FFT w Rysunek 4, gdzie częstotliwości rolek są bardziej zgrupowane, ale nie są całkowicie dyskretne. Statystyki rolki są zniekształcone, ponieważ rolka może tymczasowo powrócić do lokomocji typu dzikiego.

Wygładzone widma mocy walcowego typu C. elegans wykazują szeroki szczyt przy ~1,5 Hz, podczas gdy pływający dziki typ C. elegans wykazuje spektrum multimodalne (w tym szczyty przy ~1,0 Hz i 1,75 Hz). Fotodioda (PD) ma skończony rozmiar rozłożony na kilka elementów matrycy. Poszczególne elementy matrycy lub punkty na wzorze dyfrakcyjnym różnią się intensywnością, ponieważ interferencja konstruktywna i destrukcyjna jest różna; niemniej jednak częstotliwości, przy których zmieniają się intensywności, są takie same dla wszystkich elementów macierzy, co można zobaczyć na Rysunek 7. Biorąc pod uwagę pochodną czasową równania 1, można zauważyć, że fluktuacje częstotliwości nie zależą od macierzy fazowej, a jedynie od fluktuacji pierwotnego obiektu:

, (2)

Gdy PD rozprzestrzenia się na kilka elementów matrycy, lokalizacje szczytów uśredniają się do spójnego profilu częstotliwości. Można się spodziewać pewnych zmian, które mogą dać wskazówki co do orientacji robaka. Rozkład częstotliwości zmieni się wraz ze zmianą chodu robaka. Obecny model jest prostym modelem, który pozwala jedynie na ocenę lokalizacji pików, a nie względnych wysokości pików. Różne wzorce lokomotoryczne będą średnio docierać do różnych lokalizacji szczytów.

Rysunek 1. Konfiguracja eksperymentalna. Wiązka lasera o małej mocy przechodzi przez filtr o neutralnej gęstości, jest odbijana przez lustro M1 w dół przez kuwetę zawierającą robaka na zwierciadło M2 i przemieszcza się w kierunku fotodiody. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2. sekwencyjne kształty robaków i odpowiadające im wzory dyfrakcyjne. (a) Niektóre wybierają sekwencyjne obrazy binarne modelowanych nicieni w kształcie litery W i (b) odpowiadające im sekwencyjne wzory dyfrakcyjne. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3. Sygnatury dyfrakcyjne próbek eksperymentalnych. Sygnatury dyfrakcyjne zebrane dla (a) OH7547 "roller" i (b) N2 wild type C. elegans przy użyciu pojedynczej fotodiody we wzorze dyfrakcyjnym. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4. Eksperymentalne uśrednione widma mocy serii dyfrakcyjnej Fraunhofera typu walcowego i dzikiego. Widma pokazują częstotliwości występujące w uśrednionej transformacie Fouriera szeregu czasowego zarejestrowanego za pomocą fotodiody. Do wygładzania używany jest filtr Gaussa o odchyleniu standardowym 0,075 Hz, obcięty na 3 odchylenia standardowe. Zwróć uwagę na szeroki pik widmowy przy ~1,5 Hz w wygładzonym widmie rolkowym, w porównaniu z multimodalnym wygładzonym widmem typu dzikiego (w tym piki przy ~1,0 i 1,75 Hz). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5. Ilustracja formacji dyfrakcyjnej. Wzory dyfrakcyjne można modelować, myśląc o każdym segmencie linii jako o nieskończenie małej linii prostej (po lewej). Nałożenie tych linii (po prawej) demonstruje konstrukcję wzorów dyfrakcyjnych dalekiego pola generowanych przez nicienie w kształcie litery C. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6. Symulowane widma mocy serii dyfrakcyjnej Fraunhofera typu walcowego i dzikiego. (a) robaki w kształcie litery C oraz (b) w kształcie litery W z fotodiodą wyśrodkowaną na elementach matrycy 200 (pionowo) i 175 (poziomo). Kształt W pokazuje większą gęstość częstotliwości ze względu na bardziej złożoną lokomocję. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 7. Symulowane widma mocy serii dyfrakcyjnej Fraunhofera typu walcowego i dzikiego w różnych lokalizacjach fotodiody. (a) ślimak w kształcie litery W i robak w kształcie litery C dla pojedynczych elementów matrycy w różnych miejscach symulujących różne lokalizacje fotodiody. Wysokości szczytów różnią się w zależności od lokalizacji; Jednak lokalizacje pików pozostają takie same dla określonych kształtów. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Dyskusja

Uwzględnienie odcinków danych z brakiem aktywności spowoduje zniekształcenie wyników, ponieważ sztuczne niższe częstotliwości zostaną uśrednione w wynikach. Nasycenie fotodiody można rozpoznać po płaskich pikach lub "odciętych" pikach w surowych danych. Podzielenie każdego zestawu surowych danych przez szczytowe intensywności pomoże w uwzględnieniu wahań intensywności lasera.

Częstotliwości szczytowe są wskaźnikiem ogólnej częstotliwości młócenia; Jednak skomplikowany ruch powoduje zakłócenia przy częstotliwościach dudnień we wzorze dyfrakcyjnym i należy go dokładnie zbadać.

Metodę tę można wykorzystać do zbadania lokomocji innych nicieni. Środowisko można zmienić na inne medium. Długości fal również mogą ulec zmianie. Praca w widzialnym zakresie widma elektromagnetycznego jest najłatwiejsza i najbezpieczniejsza.

Bardziej wyrafinowany model będzie w przyszłości symulować widma dyfrakcyjne w bardziej realistyczny sposób. Przyszły model może zawierać robaka, który może zmieniać orientację, co nie wpłynęłoby na lokalizacje częstotliwości, ale na względne wysokości szczytów. Bardziej realistyczny model pozwoliłby na probabilistyczny rozkład częstotliwości młócenia, co poszerzyłoby szczyty, tak jak w danych eksperymentalnych. Rozpiętość częstotliwości odpowiadałaby za zmiany w częstotliwościach młócenia.

Obecny kształt ślimaka jest surowy, szczególnie w okolicy głowy i ogona, które powinny być bardziej zwężające się niż w obecnym modelu. Interesujące może być przeprowadzenie szczegółowej analizy szeregów czasowych sygnału, ponieważ może to dostarczyć wskazówek na temat złożoności lokomocji u różnych mutantów.

Warto zastanowić się nad praktycznością rozszerzenia tej techniki na charakteryzowanie wielu nicieni jednocześnie. Metodę tę należy rozumieć jako metodę uzupełniającą istniejące metody z wykorzystaniem tradycyjnych mikroskopów. Ta metoda ma tę zaletę, że nie wymaga mikroskopu podczas pozyskiwania danych, dzięki czemu robak może wydostać się poza płaszczyznę ogniskową. Uśrednione widma częstotliwości pokazują wyraźne różnice w ruchu robaków i mogą być określone ilościowo za pomocą przeważających pików częstotliwości, co jest nowatorską metodą ilościowego określania lokomocji robaków. Analiza danych dotyczących sygnatur dyfrakcyjnych jest w dalszym rozwoju i miejmy nadzieję, że doprowadzi do zautomatyzowanego procesu identyfikacji wielu mutantów i osobników.

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
Przestrajalny laser helowo-neonowy	Elektrooptyka badawcza	30602	Można wybrać cztery długości fal od 543 nm do 633 nm.
2 Aluminiowe lusterka z przednią powierzchnią	Thorlabs	PF10-03-F01
Fotodioda: SI Ampliified Detector	Thorlabs	PDA 100A
Kwarcowe kuwety	Starna Cells	21 / G / 5	Można również użyć plastikowych ogniw
MatLab (Oprogramowanie)	MathWorks	R2016b (9.1.0.441655)	Użyj polecenia fft, aby symulować dyfrakcję
Excel	Microsoft	14.7.1	Służy do analizy danych Rysunek 4
Caenorhabditis elegans Roller	>University of Minnesota Caenorhabditis elegans Center (CGC)	Szczep: OH7547 Genotyp: otIs199.	https://cbs.umn.edu/cgc/home
Caenorhabditis elegans Typ dziki	University of Minnesota Caenorhabditis elegans Center (CGC)	Szczep: N2 Genotyp: C. elegans izolat dziki	https://cbs.umn.edu/cgc/home