1. Przygotowanie kultur bakteryjnych i płytek do ogólnej hodowli nicieni
- Przygotować 10 mL pożywki LB (lysogeny broth) w szklanej butelce o pojemności 250 mL i wysterylizować w autoklawie.
- Zaszczepić butelkę pojedynczą kolonią szczepu Escherichia coli OP50 i inkubować przez noc w temperaturze 37 °C, a następnie przechowywać kulturę w temperaturze 4 °C do momentu użycia.
- Przygotować 5 L sterylnego agaru NGM (nematode growth medium)15 i rozlać go w objętościach 12 mL do sterylnych, plastikowych szalek Petriego o średnicy 6 cm.
- Pozostawić agar do zastygnięcia przez noc, a następnie przechowywać w temperaturze 4 °C do momentu użycia.
- Zasiać płytki NGM co najmniej 3 dni przed użyciem, nakładając 25 μL schłodzonej kultury OP50. Przechowywać płytki w temperaturze 20 °C do momentu użycia.
2. Przygotowanie kultur bakteryjnych i rozcieńczeń do doświadczeń
- Przygotuj dwie porcje po 50 ml pożywki LB w kolbie Erlenmeyera o pojemności 2 l i wysterylizuj w autoklawie.
- Zaszczep każdą kolbę pojedynczą kolonią szczepu OP50 i inkubuj w 37 °C inkubacja z wytrząsaniem przy 20 obr./min przez około 14 h.
- Uzupełnić hodowle antybiotykiem streptomycyną do stężenia końcowego wynoszącego 50 μg/mLKontynuować wstrząsanie przez 30 min w temperaturze 37 °C następnie umieścić kolby w lodzie na 15 min.
- Przenieś 450 mL każdej kultury do oddzielnych sterylnych butelek do wirowania (najlepiej o pojemności 50 mL, aby uniknąć rozdzielania kultur). Pozostałą część kultury przechowuj na lodzie, ponieważ zostanie ona wykorzystana do oznaczenia stężenia bakterii.
- Odwirować zawartość probówek z prędkością 450 x g g w wirówce chłodzącej ustawionej na 4 °C przez 25 min. Odrzucić supernatant i przechowywać probówki w lodzie.
- Rozcieńczyć 100 μL pozostałej kultury z każdej kolby do 900 μL sterylnego LB. Należy użyć 1 ml sterylnego LB do wyzerowania spektrofotometru, a następnie wyznaczyć wartość OD600 z 10-krotnego rozcieńczenia każdej hodowli. Jeśli na przykład wartość OD600 Jeśli wartość gęstości optycznej (OD) dla 10-krotnego rozcieńczenia hodowli nocnej wynosi 0,28, to rzeczywista wartość OD hodowli wynosiła600 wynoszącego 2,8.
- Użyj roztworu roboczego bakterii o stężeniu 1,12 x 1010 komórek OP50/ml, co odpowiada gęstości optycznej (OD)600 z 56. W związku z tym, wykorzystując przykład OD600 z punktu 2.8 powyżej, resuspender osad z 450-ml hodowli w 20th objętości pierwotnej, która w tym przypadku wynosi 2,50 ml. Ponownie zawiesić bakterie w sterylnym roztworze podstawowym S z dodatkiem streptomycyny do 50 μg/mL (SB + strep).
- Resuspenduj osad w odpowiedniej objętości sterylnego medium SB + strep.
- Przygotować wszystkie kolejne stężenia wykorzystywane w eksperymentach poprzez rozcieńczenia szeregowe zapasu roboczego w pożywce SB + strep, stosując współczynniki rozcieńczenia podane w Tabeli 1.
3. Przygotowanie eksperymentów dotyczących długości życia
UWAGA: Badania długości życia przeprowadza się na specjalnie przygotowanych 6 cm płytkach do hodowli tkankowej wypełnionych 12 mL agaru NGM z dodatkiem streptomycyny i karbenicyliny (NSC), oba w końcowym stężeniu 50 μg/mL. Płytki do hodowli tkankowej są bardzo dobrze dostosowane do badań długości życia przy zerowym lub bardzo niskim stężeniu bakterii, ponieważ nicienie znacznie rzadziej przywierają do ścianek płytki i ulegają wysuszeniu. Zastosowanie dwóch różnych antybiotyków zapobiega wykształceniu oporności na leki u szczepu OP50, co jest kluczowe dla kontroli stężenia bakterii na płytce. Ponadto, poprzez inaktywację bakterii antybiotykami, minimalizujemy wpływ infekcji patogennych na długość życia nicieni16. Pozwala nam to brać pod uwagę w tych eksperymentach wyłącznie komponenty stężenia bakterii związane z pożywieniem. W wielu badaniach nad starzeniem się C. elegans płytki NGM uzupełnia się chemicznie fluoro-2′-deoksyurydyną (FUdR), aby uczynić dorosłe osobniki sterylnymi. Jednak stosowanie FUdR może być problematyczne, ponieważ może powodować szereg czynników zakłócających w badaniach długowieczności17,18,19,20,21. Entchev i wsp. (2015)12 ominęli ten problem, poddając larwy L4 działaniu RNAi genu egg-522,23, co hamuje przejście oocytu w embrion poprzez blokowanie formowania się skorupki jajowych zapłodnionych oocytów C. elegans, prowadząc do ich obumarcia.
- Hodować szczepy C. elegans przeznaczone do analizy na wysianych płytkach NGM w temperaturze 20 °C. Pozostawić płytki do momentu wyczerpania pokarmu, aby zapewnić jednorodny wzrost wszystkich zwierząt oraz potencjalnie wyeliminować transgeneracyjne efekty warunków rozwojowych.
- Przenieść głodujące larwy w stadium L1 na nowe wysiane płytki NGM, wycinając mały kawałek agaru (5 mm x 5 mm) z głodnej płytki za pomocą sterylnego skalpela i kładąc go na nowej płytce; proces ten jest określany w społeczności badawczej C. elegans jako „chunking”15. Hodować nicienie, aż osiągną stadium L4. Następnie pięć larw L4 z każdego szczepu przenosi się na oddzielne wysiane płytki NGM i utrzymuje w temperaturze 20 °C. Zwierzęta te stanowią pokolenie P0.
- Wykorzystać potomstwo L4 pokolenia P0 do utworzenia pokolenia F1. Pojedynczo przenieść larwy F1 L4 szczepu dzikiego typu N2 na 30 wysianych płytek NGM.
UWAGA: W przypadku szczepów mutantów liczba wymaganych płytek zależy od tempa ich wzrostu. Na przykład zwierzęta daf-7(ok3125) przechodzą przejściowe zatrzymanie w stadium dauer w temperaturze 20 °C. W związku z tym szczep ten produkuje mniej larw L4 niż płytka z nicniami dzikiego typu N2.
- Aby zrekompensować tę różnicę, przygotować płytki dla daf-7(ok3125) dzień wcześniej niż płytki dla N2 i w większej liczbie; odpowiedni stosunek to 3:1.
- Przenieść zsynchronizowane potomstwo rodziców F1 w stadium L4 na płytki NGM uzupełnione o 1 mM IPTG i 50 µg/mL karbenicyliny (płytki RNAi), które zostały wysiane 25 µL bakterii HT15 wykazujących ekspresję dsRNA celującego w gen egg-5 i utrzymywane w temperaturze 20 °C przez 24 h. Na każdy eksperyment wymagane jest co najmniej 360 nicieni na szczep, przy gęstości 15 zwierząt na płytkę.
- Po traktowaniu RNAi egg-5 przenieść nicienie na płytki NSC wysiane 25 µL streptomycyną traktowanego szczepu OP50 w stężeniu 2 x 109 komórek/mL (Tabela 1) na dodatkowe 24 godziny w temperaturze 20 °C. Aby uniknąć uszkodzeń mechanicznych nicieni podczas tego i wszystkich kolejnych przeniesień, delikatnie podnieść zwierzęta od spodu, używając bardzo cienkiej, lekko wygiętej igły do nicieni. Zwierzęta zdjąć z igły poprzez zanurzenie jej w 10 µL kropli SB + strep umieszczonej na powierzchni nowej płytki.
- Następnego dnia ponownie rozdzielić nicienie na płytki NSC reprezentujące każdy z poziomów pokarmu wymienionych w Tabeli 1, a także na zestaw płytek NSC niezawierających bakterii. Wszystkie płytki NSC wysiać 225 µL bakterii o odpowiednim stężeniu, a płytki NSC wolne od bakterii wysiać 25 µL SB + strep. Utrzymać gęstość 15 zwierząt na płytkę, co daje co najmniej cztery płytki na każdy warunek pokarmowy dla każdego szczepu. Przenieść płytki do żądanej temperatury eksperymentalnej.
- Przenosić nicienie na świeże płytki NSC wysiane odpowiednim stężeniem pokarmu zgodnie z harmonogramem przedstawionym w Tabeli 2.
- Oceniać ruch zwierząt poprzez delikatne nakłuwanie drucianą igłą. Brak reakcji oceniać jako śmierć. Oceniać przeżywalność zwierząt przy każdym przeniesieniu, a następnie codziennie po ostatnim punkcie przeniesienia.
4. Konfiguracja eksperymentów obrazowania
UWAGA: Kroki opisane w tej sekcji są wystarczające do wygenerowania odpowiedniej liczby nicieni dla każdego szczepu, aby przeprowadzić obrazowanie jednego eksperymentalnego warunku żywieniowego w danej temperaturze. Protokół ten można przeskalować w zależności od liczby warunków, które będą obrazowane w danym dniu. Należy jednak zadbać o taki projekt eksperymentalny, aby ramy czasowe były rozsądne, a różnica wieku zwierząt między różnymi szczepami nie przekraczała 12 godzin w dniu obrazowania. Wysoce zaleca się, aby obrazowane reportery transkrypcyjne były transgenikami pojedynczej kopii, ponieważ będzie to bardziej przypominać naturalną regulację genu. Protokół opisany poniżej i podsumowany w Tabeli 3 dostarcza również uproszczonej metody otrzymywania dużych, zsynchronizowanych i dopasowanych wiekowo populacji szczepów, które mogą być wykorzystane w innych schematach doświadczalnych.
- Hoduj zwierzęta do eksperymentów obrazowych na 10 cm modyfikowanych płytkach do hodowli tkanek, aby umożliwić większą gęstość osobników (∼10 zwierząt na płytkę) niż w testach długości życia.
- Wypełnij płytki 10 cm 30 mL agaru NGM i zaszczep pięcioma alikwotami po 25 µL schłodzonej kultury OP50, rozmieszczając je w kształcie krzyża.
5. Wstępna hodowla szczepów reporterowych
- Hodować badane szczepy C. elegans na wysianych płytkach NGM w temperaturze 20 °C. Pozostawić płytki do momentu wyczerpania pożywienia, aby zapewnić jednolity wzrost wszystkich zwierząt i potencjalnie uwzględnić transgeneracyjne efekty warunków rozwojowych.
- Przenieść zagłodzone larwy w stadium L1 na wysiane płytki NGM i hodować do momentu osiągnięcia stadium L4. Przenieść po trzy larwy L4 każdego szczepu na dwie wysiane płytki NGM w temperaturze 20 °C. Zwierzęta te stanowią generację P0.
- Wykorzystać potomstwo L4 z generacji P0 do utworzenia generacji F1. Umieścić larwy L4 generacji F1 szczepu reporterowego typu dzikiego na 4 wysianych płytkach NGM. W przypadku szczepów mutantów liczba wymaganych płytek zależy od tempa ich wzrostu. Odwołując się do poprzedniego przykładu daf-7(ok3125), szczepy reporterowe na tym tle wymagają trzykrotnej liczby płytek i muszą być przygotowane dwa dni wcześniej niż szczepy reporterowe typu dzikiego, aby uwzględnić różnice w tempie wzrostu.
6. Zbieranie jaj i synchronizacja szczepów reporterowych
UWAGA: Niektóre geny będące przedmiotem badań w społeczności zajmującej się starzeniem C. elegans w przypadku mutacji prowadzą do wad wzrostu i składania jaj, co utrudnia generowanie dużych, synchronicznych populacji zwierząt o różnych genotypach. Na przykład mutacja daf-7(ok3125) powoduje poważne zaburzenia składania jaj w porównaniu ze szczepem dzikim N2. W związku z tym, aby uzyskać wystarczającą liczbę synchronicznych larw L4 różnych szczepów do ilościowych eksperymentów obrazowych, wymagana jest bardziej skuteczna metodologia niż ręczne wyłapywanie nicieni. Z tego powodu szczepy z reporterami transkrypcyjnymi poddano traktowaniu roztworem podchlorynu sodu (NaClO)/wodorotlenku sodu (NaOH) w przypadku dorosłych osobników ciężarnych, aby rozbić ciała zwierząt i uwolnić ich jaja; proces ten jest powszechnie określany w społeczności badawczej C. elegans jako „bleaching” (wybielanie)15.
- Uwzględnij różnice w tempie wzrostu poszczególnych szczepów i oblicz, kiedy należy zebrać i poddać wybielaniu płyt z każdym szczepem. Przykład czasu wybielania szczepów typu dzikiego oraz szczepów reporterowych daf-7(ok3125) znajduje się w Tabeli 3.
- W odpowiednim dniu sukcesywnie zmyj nicienie danego szczepu z płyt, używając 15 mL SB, i zbierz je do sterylnej probówki 15 mL. Pozostaw zwierzęta do naturalnego osadzenia, a następnie dostosuj objętość cieczy do 7 mL.
- Do probówki zawierającej dorosłe osobniki z jajami dodaj 2 mL 5% NaClO oraz 1 mL 5 M NaOH. Delikatnie kołysz mieszaniną w temperaturze pokojowej nie dłużej niż przez 3 min. NaClO zabije bakterie i nicienie, natomiast NaOH spowoduje rozpad organizmów nicieni, uwalniając zawarte w nich zapłodnione jaja do cieczy. Chitynowa skorupka jaj otaczająca embriony chroni je przed działaniem odczynników, pod warunkiem, że czas ekspozycji pozostanie stosunkowo krótki. Po 3-minutowej inkubacji wymieszaj zawartość probówki w wirówce typu vortex przez ~30 s, aby ułatwić dalszy rozpad ciał nicieni.
- Po 3-minutowej inkubacji odwiruj mieszaninę przy 1 0 x g przez 1 min, aby uzyskać pellet jaj. Odessaj większość supernatantu za pomocą sterylnej pipety szklanej podłączonej do kolby próżniowej, pozostawiając ~0,5 mL w każdej probówce, aby nie zakłócić osadu jaj.
- Resuspenduj pellet w 9,5 mL SB, a następnie powtórz etap wirowania i resuspensji jeszcze dwa razy, tak aby jaja zostały łącznie trzykrotnie przemyte SB.
- Po ostatnim płukaniu oddziel jaja przez wirowanie, a następnie odlej cały supernatant z wyjątkiem 0,5 mL. Resuspenduj jaja w pozostałych 0,5 mL SB, a następnie nanieś aliquoty po 100 µL na trzy posiane płytki NGM o średnicy 10 cm. Rozdziel 100 µL równomiernie na wszystkie pięć obszarów wzrostu bakterii na każdej płytce. W przypadku szczepów typu dzikiego nie należy nanosić więcej niż 20 jaj na płytkę.
- Przechowuj płytki w temperaturze 20 °C przez 48 h, aby uzyskać wysoce homogenną populację larw w stadium L4. W przypadku szczepów z fenotypem opóźnionego wzrostu zaleca się naniesienie jak największej liczby dostępnych jaj i wydłużenie czasu inkubacji. Na przykład szczepy reporterowe zawierające daf-7(ok3125) należy inkubować w temperaturze 20 °C przez ~64 h, aby umożliwić wyklucie się jaj i osiągnięcie stadium L4.
7. Traktowanie szczepów reporterowych za pomocą RNAi egg-5
- Kolejno zmyj nicienie z trzech płytek, używając 15 mL SB, i zbierz ciecz do sterylnej probówki 15 ml. Pozostaw larwy L4 do naturalnego osadzenia, a następnie odessij całą ciecz z wyjątkiem około 0,5 mL. W szczepach reporterowych typu dzikiego krok ten usuwa wszystkie larwy młodsze niż L4. W przypadku mutantów krok ten wspomaga usuwanie larw z zahamowanym rozwojem, takich jak dauery w przypadku daf-7(ok3125).
- Resuspenduj nicienie w 9 mL SB. Ponownie monitoruj tempo sedymentacji larw L4, a następnie odessij cały nadsącz z wyjątkiem około 0,5 mL, gdy większość larw L4 utworzy osad. Powtórz ten proces jeszcze raz. Następnie odessij całą ciecz z wyjątkiem około 0,5 mL, gdy większość larw L4 opadnie na dno.
- Dodaj 10 µL sterylnego SB uzupełnionego o 0,1% Pluronic F-127 (SB + Plu) do cieczy zawierającej larwy L4. Substancja ta działa jako surfaktant i zapobiega przywieraniu larw do wewnętrznej powierzchni plastikowych końcówek pipet.
- Delikatnie resuspenduj larwy, używając końcówki do pipety P20 o niskiej retencji, a następnie nanieś alikwoty po 150 µL na trzy płytki RNAi 10 cm, zasiające 5 x 25 µL bakterii RNAi egg-5. Upewnij się, że nicienie są równomiernie rozłożone na wszystkich pięciu trawnikach bakteryjnych.
- Po wchłonięciu cieczy przez agar, usuń z płytek wszystkie larwy inne niż L4, których nie wyeliminowano podczas procedury płukania, zbierając je manualnie. Następnie przechowuj płytki w temperaturze 20 °C przez 24 h.
8. Inicjacja szerokozakresowej DR
UWAGA: Po 24 h traktowania RNAi egg-5, larwy w stadium L4, które początkowo naniesiono na płytki, staną się jednodniowymi dorosłymi osobnikami.
- Na tym etapie usuń wszystkie młode larwy, które przeżyły poprzedni etap ręcznego usuwania, tak aby na płytkach pozostały jedynie 1-dniowe osobniki dorosłe.
- Dla każdego szczepu zmyj 1-dniowe osobniki dorosłe z trzech płytek za pomocą 15 mL sterylnego SB + strep do probówki o pojemności 15 mL. Pozostaw nicie do naturalnego osadzenia, a następnie odessaj cały nadsącz z wyjątkiem 0,5 mL. Resuspenduj nicie w 9,5 mL SB + strep i powtórz etapy sedymentacji oraz płukania.
- Po ostatnim płukaniu pozwól niciom osąść, a następnie odessaj cały nadsącz z wyjątkiem 0,5 mL. Dodaj 10 µL SB + Plu i delikatnie resuspenduj larwy za pomocą końcówki pipety P20, a następnie nanieś 10 µL na płytkę NSC zasiającą 5 x 25 µL bakterii w stężeniu 2 x 109 komórek/mL.
- Rozdziel 10 µL równomiernie na wszystkie pięć kolonii bakteryjnych. Pod mikroskopem oszacuj liczbę zwierząt obecnych na płytce: celem jest uzyskanie od 10 do 150 nicni na płytce.
- Określ objętość cieczy wymaganą do uzyskania gęstości nicni w tym zakresie, a następnie nanieś ją na dwie dodatkowe płytki. Dostosuj liczbę nicni na pierwszej płytce, aby również mieściła się w tym zakresie, a następnie przechowuj płytki w temperaturze 20 °C przez 24 h.
- Następnego dnia zbierz 2-dniowe osobniki dorosłe i przenieś je na nowe płytki NSC zasiające pożądaną eksperymentalną stężenie pożywki (Tabela 1), stosując metody opisane w krokach 2 i 3. Po wchłonięciu cieczy przez agar przenieś płytki do pożądanej temperatury eksperymentalnej na 24 h.
- Następnego dnia zbierz 3-dniowe osobniki dorosłe i przenieś je na świeże płytki NSC zasiające to samo eksperymentalne stężenie pożywki, powtarzając metody opisane w krokach 8.3 i 8.4. Po wchłonięciu cieczy przez agar przywróć płytki do temperatury eksperymentalnej na 48 h.
- Zbierz 5-dniowe osobniki dorosłe i przenieś je na świeże płytki NSC zasiające to samo eksperymentalne stężenie pożywki, powtarzając metody opisane w krokach 8.3 i 8.4. Po wchłonięciu cieczy przez agar przywróć płytki do temperatury eksperymentalnej na 24 h.
9. Mikrofluidycyjne obrazowanie szczepów reporterowych
- W 6. dniu dorosłości przeprowadź obrazowanie zwierząt przy użyciu niestandardowej platformy mikroprzepływowej24,25. Ręcznie pobierz zwierzęta z trzech płytek i przenieś je do probówki kriogenicznej o pojemności 5 mL zawierającej 4,5 mL pożywki SB + strep.
- Po osadzeniu się nicieni odessij cały nadsącz z wyjątkiem ok. 0,5 mL i resuspenduj zwierzęta w 4 mL pożywki SB + strep. Płukanie to pozwala usunąć nadmiar bakterii, które mogłyby zakłócić obrazowanie. Nicienie są wprowadzane do niestandardowego urządzenia mikroprzepływowego za pomocą przepływu wymuszonego ciśnieniem24,25. Wewnątrz urządzenia poszczególne nicienie są kierowane do kanału obrazowania i uwięzione w nim za pomocą zaworów na chipie sterowanych ciśnieniem26 pod kontrolą dedykowanego oprogramowania.
- Po uwięzieniu nicienia głową w kanale obrazowania, wykonaj fluorescencyjny stos obrazów z-stack, składający się z 50 sekcji w odstępach 2 µm, używając standardowego mikroskopu epifluorescencyjnego z obiektywem olejowym 40X (1,3 NA) oraz kamery. Obrazy z fluorescencją czerwoną i zieloną dla każdego z reporterów transkrypcyjnych zbieraj jednocześnie przy użyciu rozdzielacza emisji i przechowuj do analizy. Akwizycja obrazów jest zautomatyzowana za pomocą dedykowanego oprogramowania.
- Przetwórz obraz automatycznie przy użyciu niestandardowych skryptów MATLAB27 (dostępnych pod adresem https://github.com/meizhan/SVMelegans). Stosy Z-stack są ładowane do programu MATLAB i analizowane w celu identyfikacji par neuronów oraz ich lokalizacji w płaszczyźnie obrazowania. Następnie obliczane są projekcje maksymalne, a algorytm progowania jest wykorzystywany do zlokalizowania poszczególnych fluorescencyjnych komórek. Identyfikacja komórek jest następnie obliczana na podstawie względnych odległości i lokalizacji w głowie nicienia.
- Aby ilościowo określić fluorescencję reportera, wyodrębnij trójwymiarową objętość wokół każdej lokalizacji komórkowej ze stosu z-stack. Zintegruj intensywność dla stałej liczby najjaśniejszych pikseli, które w każdym przypadku w pełni obejmują całą komórkę.
- Aby wyeliminować zakłócenia wynikające ze zmian w autofluorescencji jelit każdego zwierzęcia, specyficznych dla warunków doświadczalnych lub szczepu, oblicz intensywność tła dla par komórek położonych najbliżej jelit (ADF i ASI) poprzez oszacowanie mody rozkładu intensywności w objętości wokół neuronu. Odejmij tę wartość intensywności tła od zintegrowanej fluorescencji, aby uzyskać wynik końcowy.
10. Montaż danych
UWAGA: Intensywności fluorescencji wszystkich neuronów analizowanych przez oprogramowanie do przetwarzania obrazów są łączone w pliku przefiltrowanych danych ekspresji (FED), który służy do szacowania profili rozkładu ekspresji genów (szablony skryptów w językach R i C++ są dostępne pod adresem https://github.com/giovannidiana/templates).
- Ręcznie sprawdź każdy przetworzony obraz, aby potwierdzić prawidłową identyfikację wszystkich komórek. Obrazy z błędną identyfikacją komórek należy odnotować w pliku wykluczeń. W pracy Diana et al. (2017)13 plik wykluczeń stanowi tablica binarna, w której „0” oznacza prawidłową, a „1” błędną identyfikację komórki. Liczba wierszy w tablicy odpowiada liczbie obrazowanych nicieni, z osobną kolumną dla każdej obrazowanej komórki (tj. ASI, ADF i NSM).
- Wygeneruj plik FED:
- Uruchom skrypt bash „gen_data+time”, aby wygenerować pliki specyficzne dla komórek, łączące wartości ekspresji uzyskane z każdego obrazowanego nicienia po przefiltrowaniu za pomocą pliku wykluczeń. Dla każdego folderu wygenerowanego przez oprogramowanie do przetwarzania obrazów, skrypt bash odczytuje plik adnotacji _EXP.txt w celu wyodrębnienia warunków eksperymentalnych oraz plik wykluczeń _X.csv, aby wybrać tylko prawidłowo zidentyfikowane neurony. Wartości ekspresji są odczytywane z plików _data_.csv.
- Połącz wszystkie pliki w jeden plik „FED_split.dat” i posortuj go według kodu partii eksperymentalnej, etykiety nicznia oraz tożsamości komórki.
- Uruchom program „sort” (program C++, użycie: ./sort FED_split.dat FED_merged.dat), aby wybrać wpisy z niezerową fluorescencją dla każdej komórki i połączyć je w pojedyncze wiersze w pliku FED_merged.dat.
1. Szacowanie ilości zakodowanych informacji
UWAGA: Poniższa procedura opisuje sposób ilościowego określania informacji o konkretnych warunkach środowiskowych zakodowanych przez zestaw ekspresji genów. W pracy Diana et al. (2017)13 badano informacje zakodowane w odniesieniu do dostępności pożywienia w środowisku, jednak sama metoda jest możliwa do zastosowania dla dowolnej dyskretnej liczby stanów środowiskowych. Kluczowym elementem niezbędnym do ilościowego określenia zmiennych teorii informacji, takich jak entropie informacji czy redundancja, jest wspólny rozkład prawdopodobieństwa odpowiedzi neuronalnych przy rozważanym zestawie bodźców środowiskowych. Aby przeprowadzić taką estymację, kluczowe jest posiadanie wystarczającej próby odpowiedzi w populacjach nicieni. Rozkłady Gaussa można estymować na podstawie stosunkowo małych próbek; jednak ważne jest posiadanie wiedzy o oczekiwanym kształcie rozkładu ekspresji, aby określić odpowiednią wielkość próby dla wiarygodnej estymacji gęstości. Ze względu na nieuniknioną zmienność między różnymi próbami, niezbędne jest sprawdzenie, czy wartości centralne rozkładów uzyskanych z różnych powtórzeń tego samego eksperymentu nie są systematycznie przesunięte oraz czy żadna z cech statystycznych rozkładu ekspresji nie uległa znaczącej zmianie między próbami. W przypadku, gdy zmienność między próbami jest porównywalna ze zmiennością wewnątrz każdej próby, kluczowe jest zbalansowanie liczby prób i liczby nicieni w obrębie prób, aby uśrednić te czynniki środowiskowe/biologiczne, które wpływają na zmienność między próbami. Niedostateczne próbkowanie tych czynników mogłoby znacząco obciążyć analizę teoretyczno-informacyjną.
- Ponieważ skrypt R „code3D.R” stanowi szablon do generowania gęstości trójwymiarowych na podstawie pliku „FED_merged.dat”, należy zmodyfikować ten szablon zgodnie z konkretnym formatem nagłówka FED. Lista „HeaderNames” zawiera nazwy poszczególnych pól w pliku FED, natomiast „RONames” jest listą odczytów. Skrypt wykorzystuje pakiet R „ks”.28,29do szacowania rozkładów wielowymiarowych w obrębie siatki hiperkostki z podziałem na GS przedziały w każdym wymiarze, dzieląc zakres między wartościami minimalną a maksymalną w zbiorze danych dla każdego odczytu. Gdy wewnętrzna zmienna „group” ma wartość 0, do szacowania gęstości wykorzystuje się wszystkie dane. Gdy etykieta „group” przyjmuje wartości od 1 do 5, zbiór danych zostaje podzielony na 5 rozłącznych zbiorów, a gęstości ekspresji są szacowane na podstawie 80% danych, wynikających z wykluczenia jednego z pięciu zbiorów. Funkcja ta jest wykorzystywana w dalszej kolejności do szacowania niepewności.
UWAGA: Uruchomienie pliku code3D.R: Rscript code3D.R <żywność> gdzie GT oznacza genotyp, food jest warunkiem środowiskowym, outfolder to istniejący folder, w którym będą przechowywane rozkłady, label to prefiks nazwy pliku, a frac to ułamek wykorzystanego zbioru danych.
- Dla każdego warunku środowiskowego wygenerowano rozkłady wielowymiarowe o różnych rozmiarach siatki GS (np. (20, 30, 40 przedziały). W celu zmniejszenia obciążenia obliczeniowego preferowane są mniejsze wartości GS, gdy wyższa rozdzielczość nie zmienia znacząco szacowania informacji. Rozkłady ekspresji genów zapisywane są w folderze określono podczas uruchamiania skryptu code3D.R jako jednokolumnowe pliki tekstowe o strukturze nazw plików __<żywność>_GSgrupa.dat.
UWAGA: Poprzednie kroki dostarczają szacunki rozkładów prawdopodobieństwa warunkowego
gdzie g oznacza wektor wszystkich odczytów, a f warunki środowiskowe. Jednak w celu obliczenia informacji wzajemnej między ekspresją genów a środowiskiem30,31.

Potrzebujemy rozkładu sygnału wejściowego
co również wyznacza średnią (wejściową) ekspresję genów
W przypadku braku bezpośredniego dostępu do rozkładu wejściowego, istotną wielkością pozwalającą na charakterystykę cech kodowania jest pojemność kanału, którą można wyznaczyć poprzez maksymalizację informacji wzajemnej dla wszystkich możliwych rozkładów wejściowych.
- Wykorzystując rozkład ekspresji genów, zastosuj algorytm Arimoto-Blahuta32 algorytm służący do szacowania pojemności kanału systemu oraz rozkładu danych wejściowych z otoczenia, który maksymalizuje ilość informacji. Przykład implementacji algorytmu znajduje się w programie w języku C++ „Ccap3D.cpp” dla trójwymiarowego kodu ekspresji genów analizowanego przez Dianę i wsp. (2017)13Przykład: jeśli rozkłady dla genotypu GT123 zostały uzyskane z 80% danych (grupa 3) przy rozmiarze siatki równym 30 i zapisane w folderze „./pdf/” z prefiksem etykiety „PDF”, komenda służąca do obliczenia informacji zakodowanej w tle GT123 to: ./Ccap3D GT123 pdf PDF 30 3
UWAGA: Na szacowanie informacji zakodowanych przez system wpływa wiele źródeł niepewności, w tym wybór algorytmu szacowania gęstości oraz błąd wynikający z wielkości próby. Kroki 1.1–1.2 należy powtórzyć, stosując różne metody szacowania gęstości, aby ocenić wielkość błędu systematycznego wprowadzanego przez każdy z algorytmów.
- Aby oszacować niepewność wynikającą z wielkości próby, należy ponownie obliczyć informacje z kroków 1.1-1.2 dla grup od 1 do 5. Zmienność w obrębie tych pięciu niezależnych prób pobranych z 80% danych odzwierciedli wpływ wielkości próby na szacowanie informacji.
- Oblicz informacje zgodnie z krokami 1.1–1.2 dla zwiększającej się frakcji danych (podobnie jak w kroku 5), aby skorygować błąd wynikający z wielkości próby (metoda jackknife).33.
12. Obliczanie redundancji, szumu i korelacji sygnałów
- Aby obliczyć informację wzajemną, należy wykorzystać optymalne rozkłady wejściowe uzyskane z oszacowania maksymalnej informacji.
między sygnałem wejściowym a odpowiedzią w ekspresji genów każdego neuronu NPlik źródłowy C++ „GetMI1D.c” jest programem przykładowym służącym do wyznaczania brzegowej informacji wzajemnej z łącznych rozkładów prawdopodobieństwa.
- Oblicz redundancję, sumując wyznaczone powyżej wartości informacji wzajemnej dla każdego neuronu, a następnie odejmując przepustowość kanału.
- Oblicz wyraz informacyjny „shuffle”
13,34Szablon pliku źródłowego C++ znajduje się w pliku „GetShuffle.c”.
- Użyj terminu „shuffle” do obliczenia korelacji sygnału i szumu. W przypadku korelacji sygnału mamy
a w przypadku korelacji szumu otrzymujemy
.
- Niepewności dotyczące redundancji; należy wyznaczyć korelację szumu i sygnału analogicznie jak w przypadku całkowitej informacji (krok 1.3 poprzedniej sekcji), wykonując wielokrotne próbkowania 80% danych i obliczając odchylenie standardowe.