Metodę tę zastosowano w celu zweryfikowania hipotezy, że curli mogą nadawać sztywność E. coli i Sbiofilmy Typhimurium, ograniczając ruch kulek podczas eksperymentów z wykorzystaniem mikroskopii konfokalnej. Obecny zestaw narzędzi został wykorzystany do porównania właściwości materiałowych Enterococcus faecalis szczep typowy komensalny OG1RF do Salmonella enterica serotyp Typhimurium, E. coli, oraz ich odpowiednie izogeniczne mutanty curli (Rycina 1B i Wideo uzupełniające 1, Wideo uzupełniające 2, Wideo uzupełniające 3, Wideo uzupełniające 4, Wideo uzupełniające 5, Wideo uzupełniające 6Właściwości materiałowe biofilmu mogą potencjalnie różnić się pod względem sztywności (np. curli związane z eDNA) lub oddziaływań elektrostatycznych i hydrofobowych pomiędzy ujemnie naładowanymi koralikami a komórkami biofilmu i materiałami macierzy, a także gęstości komórkowej.
Powtarzalność
Zestaw narzędzi Biofilm został zaprogramowany w językach Groovy30 i Java31 w środowisku VRL-Studio32, co umożliwiło zaprojektowanie modularnego przepływu pracy z automatycznym generowaniem interfejsu użytkownika (UI) dla wszystkich komponentów obliczeniowych. Pozwoliło to na automatyzację procesu, eliminując niezamierzone błędy systematyczne wprowadzone przez eksperymentatora podczas analizy wyników.
Wykorzystanie MSD do potwierdzenia rodzaju ruchu w biofilmach
Do analizy trajektorii za pomocą programu Particle Tracker 2D/3D dostępne są różne ustawienia dynamiki dla analizy różnych typów ruchu kulek. W tych badaniach wybrano ustawienie „Brownian motion” (tj. ruch dyfuzyjny), ponieważ E. faecalis jest bakterią nieporuszającą się, E. coli i Salmonella nie wykazują obecności wici w biofilmach, a eksperymenty przeprowadzono w systemie zamkniętym przy braku przepływu. Ustawienie to można było dodatkowo zweryfikować za pomocą obliczonego średniokwadratowego przemieszczenia (MSD) kulek. Stosując definicję
gdzie m jest liczbą segmentów trajektorii, można obliczyć zmianę MSD w przebiegu każdej trajektorii. Trajektorie liniowe wskazują na dyfuzyjny ruch kulek (Rycina 2A). Przy użyciu kwadratowej metody najmniejszych kwadratów obliczono średni wzorzec ruchu wszystkich kulek w biofilmie, co wykazało dominujący charakter liniowy i potwierdziło dyfuzję pasywną jako siłę napędową (Rycina 2A-2F).
Analiza ramki ograniczającej.
Zestaw narzędzi wykorzystuje ImageJ Mosaic oraz Particle Tracker 2D/3D do generowania trajektorii (Krok 4), a następnie, przy użyciu zautomatyzowanego potoku analizy biofilmu, generuje istotne dane dotyczące trajektorii kulek, które mogą służyć do porównania właściwości materiałowych biofilmu. Objętość ramki ograniczającej w µm3 mierzono poprzez konstrukcję minimalnego prostopadłościanu obejmującego trajektorię i pomiar jego objętości (Rycina 3).
E. faecalis biofilmy wykazują większy ruch kulek, z wartościami bounding box w zakresie 1-6000 µm3 (Rycina 3B, 3C i 3D). Wyniki potwierdzają, że ruch obserwowany przy użyciu szkiełka nakrywkowego zamontowanego na powleczonym preparacie z otworem o grubości powłoki około 25 µm (Rycina 3C) w porównaniu do biofilmów hodowanych na dnie szklanych studzienek optycznych i obrazowanych bezpośrednio (Rycina 3D) daje równoważne wyniki z niewielkimi różnicami. Jedyną różnicą było to, że w pobliżu górnej części zamontowanego szkiełka nakrywkowego w biofilmach E. faecalis można było zarejestrować stabilne trajektorie o czasie trwania przekraczającym 10 minut, ale jednocześnie małe wartości bounding box, podczas gdy na płytce z dnem optycznym można było zarejestrować określoną liczbę kulek o wyższej mobilności. Wspólnie sugeruje to, że montaż szklanego preparatu mógł zmienić napięcie powierzchniowe układu na styku górnej części biofilmu ze szkiełkiem w przypadku biofilmu na zamontowanym szkiełku nakrywkowym, co ostatecznie zmniejszyło mobilność niektórych kulek w obszarach biofilmu o mniejszej lepkości (Rycina 3B, 3C, 3D i 3I). Trajektorie należące do tej kategorii stanowiły bardzo niewielki procent i nawet przy tej małej liczbie uwięzionych kulek średnie MSD dla E. faecalis na zamontowanym szkiełku nakrywkowym było nieco wyższe niż MSD obliczone dla biofilmu na płytce z dnem optycznym (Rycina 3).
S. Typhimurium i E. coli wykazały mniejsze objętości ograniczających pól trajektorii kulek, wynoszące 0-10 µm3 (Rysunek 3A, 3B, 3E i 3F), w porównaniu do izogenicznych mutantów curli, których pola ograniczające wynosiły 1-6000 µm3 dla E. coli oraz 1-5000 µm3 dla S. Typhimurium (Rysunek 3A, 3B, 3F i 3H), co świadczy o większej mobilności kulek. Wyniki te sugerują, że obecność amyloidu korelowała ze zwiększoną sztywnością biofilmów, co było zgodne z brakiem zauważalnego ruchu biofilmu na nagraniach wideo. Objętości pól ograniczających były konsekwentnie małe (0-10 µm3) nawet w obszarach biofilmu o niskiej gęstości. Obserwacja ta jest zgodna z wcześniejszymi doniesieniami, że curli mogą występować w obszarach biofilmu o niskiej gęstości komórek10.
Nie było możliwe porównanie zachowania biofilmu Enterobacteriaceae na płytkach z dnem optycznym, ponieważ bakterie te rosną w formie błon (pellicles) na granicy fazy powietrze-ciecz (Krok 1.2.2). Podczas stosowania szkiełka cover-slip, błona przylegała do szkiełka na granicy faz, a po jego zdjęciu błona spoczywała na szkiełku, tworząc jednolity obszar obrazowania. W przypadku płytki z dnem optycznym hodowanej pod skosem, obrazowanie przeprowadzano przy obecności cieczy w studzience. Oznacza to, że błona wciąż unosiła się nad dnem optycznym, przez co znajdowała się poza głębokością roboczą mikroskopu odwróconego, takiego jak Leica Sp5. Usunięcie wystarczającej ilości medium, aby wprowadzić biofilm w zakres głębokości roboczej mikroskopu, powodowało wysychanie preparatu w trakcie 20-minutowego procesu obrazowania.
Ogólnie rzecz biorąc, wykresy potwierdzają obserwacje wizualne przedstawione w materiałach wideo uzupełniających i są zgodne z zaobserwowanymi różnicami w MSD (Rycina 3I i 3J).
Długość trwania trajektorii
Długość trwania trajektorii mierzono jako liczbę następujących po sobie klatek, w których zarejestrowano kuleczkę (Rycyna 3).
W bardziej lepkich, przypominających ciecz biofilmach E. faecalis wszystkie kulek miały czas trwania trajektorii krótszy niż 10 minut, a większość trajektorii w przypadku biofilmów E. faecalis mieściła się w przedziale od 2 do 5 minut. Jednakże kulki z zarejestrowanym krótkim czasem trwania trajektorii mogły zostać zlokalizowane poprzez inspekcję wizualną w biofilmach E. faecalis w całym oknie czasowym obrazowania (Supplemental Video 1 and 2). Zatem możliwe jest, że kulki poruszają się wzdłuż zarejestrowanej trajektorii, okresowo odseparowując się od biofilmu i kończąc trajektorię, a następnie ponownie łącząc się z biofilmem, w którym momencie inicjowana jest nowa trajektoria. Ostatecznie prowadziłoby to do krótkich czasów trwania trajektorii przy ciągłej obecności kulek w biofilmie. Należy zauważyć, że przy zastosowaniu tej techniki czasy trwania trajektorii, zwłaszcza w lepkim biofilmie, mają tendencję do zaniżania całkowitego czasu, przez który kulka jest związana z biofilmem.
W biofilmach S. Typhimurium, które charakteryzowały się mniejszymi objętościami ograniczających ramek (bounding box), większość kulek (około 80%) wykazywała długie czasy trwania trajektorii wynoszące 16-20 klatek, co odpowiadało około 15-20 minutom czasu rzeczywistego (Rysunek 3A, 3G i 3H). W przeciwieństwie do nich, biofilmy izogenicznych mutantów curli zawierały więcej mobilnych kulek z objętościami ograniczających ramek w zakresie od 1-6000 µm3 (E. coli) i 1-5000 µm3 (S. Typhimurium) (Rysunek 3A, 3B, 3F i 3H). W kontraście do biofilmów E. faecalis, w których >70% trajektorii miało objętości ograniczających ramek większe niż 10 µm3, biofilmy gatunków z rodziny Enterobacteriaceae wykazały jedynie 30% trajektorii kulek o objętościach ograniczających ramek powyżej 10 µm3. Mimo że całkowite czasy trwania trajektorii kulek były krótsze w biofilmach mutantów curli, niektóre trajektorie odzwierciedlały znaczący ruch kulek oraz długie czasy trwania trajektorii (Rysunek 3H). Obserwacja ta może wskazywać, że zmienność ta może odpowiadać różnym właściwościom materiałowym biofilmu, takim jak wiskoelastyczność, i/lub zmianom chemii powierzchni cząstek, np. ładunku.
Analiza długości trajektorii i prędkości kulek
Długość trajektorii jest pomiarem odległości przebytej przez kuleczki w µmPomiar ten jest zgodny z prędkością ruchu kulek w µm/sZgodnie z większymi objętościami ramki ograniczającej, kuleczki w E. faecalis biofilmy wykazywały 10-krotnie dłuższe trajektorie, 5-20 µmkontrastowo z <4 µm w biofilmach zawierających włókna curli. Zgodnie z krótszymi trajektoriami (Rycina 4A). E. faecalis kulki osiągały prędkości nawet 15-krotnie wyższe, przy czym większość z nich poruszała się z prędkościami w zakresie 0,01-0.15 µm/s prędkości względem siebie <0.006 µm/s (Rycina 4B). Niemniej jednak w biofilmach mutantów curli odnotowano ogólnie niższe prędkości i krótsze trajektorie w porównaniu do E. faecalis biofilmy, lecz dłuższe trajektorie i wyższe prędkości niż w przypadku szczepów rodzicielskich zawierających curli (Rycina 4A oraz 4B).
Warto zauważyć, że włóknista, przypominająca sieć struktura curli30 może wpływać na mobilność w sposób anizotropowy, ograniczając ruch w płaszczyźnie xy i umożliwiając zwiększoną mobilność w kierunku z (Rysunek 2G). Duża pula trajektorii (około 800) należąca do w przybliżeniu 50 unikalnych kuleczek w biofilmach zawierających curli byłaby spójna z ograniczeniami metody Mosaic Particle Tracking, która zlicza każdą z tych szybko poruszających się kuleczek jako pojedynczą kuleczkę w osiach x, y i z. Do potwierdzenia tej obserwacji niezbędne będą dalsze badania oraz rozwój oprogramowania.
Analiza zależności ruchu kulek od gęstości komórkowej
Zależność ruchu kulek od gęstości komórkowej określono przy użyciu średnich ważonych prędkości i ich wariancji, a także średnich/średnich ważonych i wariancji objętości ramek ograniczających (bounding box). Do obliczenia lokalnych gęstości komórkowych w obliczeniach prędkości ważonych wykorzystano drugi kanał obrazowania dla bakterii znakowanych Syto9. Gęstość komórkową obliczono poprzez uśrednienie danych wokselowych Syto9 w obrębie ramki ograniczającej każdego krawędziowego odcinka trajektorii (Rysunek 5, prawa strona). W ten sposób prędkość kulek może być ważona lokalną gęstością komórek na każdym odcinku. Istnieje wiele rodzajów barwników, które można wykorzystać do wizualizacji bakterii, w tym barwniki ścian komórkowych, błon oraz zawartości DNA. Do wyznaczenia gęstości komórkowej wybrano Syto9, ponieważ daje on najbardziej spójny sygnał niezależnie od wizualizowanego przekroju optycznego w osi Z. Barwniki otoczkowe (ściany komórkowej i błony) będą dawać różny sygnał w zależności od pozycji przekroju Z. Jeśli przekrój Z obejmuje górną lub dolną część komórki, sygnał będzie silniejszy niż w przypadku przekroju Z przechodzącego przez środek komórki, gdzie zabarwiony jest tylko jej zarys.
Trajektorie kulek z kanału czerwonego śledzono przez 20 klatek, przy czym pojedyncze trajektorie miały minimalny czas życia wynoszący 2 klatki i maksymalny 20 klatek, z 19 segmentami trajektorii łączącymi klatki (Rycina 5). Aby zbadać zależność mobilności kulek od gęstości komórek, wyznaczono intensywność GFP na voksel (każdy z pojedynczych pomiarów w obrazie o rozmiarach 512x512 vokseli). Gęstość komórek wokół każdego segmentu trajektorii kulek obliczono jako średnią lokalną gęstość w ramce otaczającej dany segment.
W przypadku niektórych biofilmów można było odnotować statystycznie istotną zależność od gęstości (Rysunek 6), najsilniej w przypadku biofilmów E. faecalis hodowanych na szklanym zakryciu i odwróconych na szkiełku wieloodworowym (Rysunek 6A). Przeciwnie, biofilmy E. faecalis hodowane na dnie płytki 96-dołkowej (Rysunek 6B) nie wykazały zależności od gęstości. Podsumowując, sugeruje to, że wysoce płynne biofilmy E. faecalis mogły zostać potencjalnie nieznacznie ściśnięte z powodu zamontowania ich na szkiełku wieloodworowym, co jest zgodne z redukcją liczby kulek poruszających się szybciej oraz tych o niewielkich objętościach ramki ograniczającej, które zostały uwięzione u góry biofilmu przy szklanym szkiełku (Rysunek 3C vs Rysunek 3D). Zarówno biofilmy Salmonella i E. coli (Rysunek 6C i 6D), jak i ich mutanty izogeniczne (Rysunek 6E i 6F) wykazały niewielką lub żadną zależność od gęstości komórkowej.

Rycina 1. Schemat obrazowania i analizy (kroki 2-4) (A) Biofilmy obrazowano zgodnie z opisem w sekcji 2.2. Na podstawie obrazów biofilmów (patrz pkt 3) wygenerowano trajektorie kulek zgodnie z opisem w sekcji 4. Wykorzystując trajektorie, obliczono istotne dane za pomocą narzędzi analitycznych (patrz pkt 5) (B) Biofilmy hodowano zgodnie z opisem w kroku 1 na szkiełkach nakrywkowych (E. faecalis, Wideo 1), S. Typhimurium (Wideo 3), E. coli (Wideo 5) oraz izogenicznymi mutantami curli (Wideo 4, Wideo 6) lub w 96-dołkowej płytce z dnem optycznym (E. faecalis dno, Wideo 2). Biały pasek skali odpowiada 20 mm. Rycina powtórzona za zgodą z (31). Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 2. Stwierdzono, że ruch kulek w biofilmach ma charakter dyfuzyjny (ruch Browna, patrz krok 5) (A-F) Dane MSD wykazują zachowanie liniowe (czerwona linia), co potwierdza ruch Browna. Biofilmy hodowano zgodnie z krokiem 1 (G) Przykład eliptycznego ruchu kulek zaobserwowanego w biofilmach E. coli oraz S. Typhimurium, pochodzący z jednej klatki testu 4D biofilmu E. coli (H) Przykład dużych zmian w wzorach rozmieszczenia kulek między klatką 3 a 4, pochodzący z optycznego dołka dna dla E. faecalis. Należy zauważyć, że sam biofilm wywołuje pewien przepływ (Wideo 1 i 2), co sprawia wrażenie, że klatka 3 i 4 są zorientowane inaczej. Rysunek przedrukowano za zgodą z (31). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 3. Analiza różnic w sztywności przy użyciu ramek ograniczających i czasu trwania trajektorii. Biofilmy hodowano zgodnie z opisem w kroku 1 na szkiełkach nakrywkach (E. faecalis, Wideo 1), S. Typhimurium (Wideo 3), E. coli (Wideo 5) i izogenicznych mutantach curli (Wideo 4, Wideo 6) lub w 96-dołkowej płytce z dnem optycznym (dno E. faecalis, Wideo 2). Czas trwania trajektorii przedstawiono jako % całkowitej liczby trajektorii kulek (A) oraz na wykresach rozrzutu (C-H), wraz z objętościami ramek ograniczających (obliczonymi w kroku 5) (I). Porównanie MSD kulek w różnych biofilmach (H) Średnie MSD dla każdego typu biofilmu. Słupki wskazują 95% przedział ufności danych. Rysunek powielony za zgodą z (31). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 4. Analiza różnic w długości trajektorii i prędkości kulek. Biofilmy hodowano zgodnie z opisem w kroku 1. Długości trajektorii podano w µm i przedstawiono jako % całkowitej liczby trajektorii kulek (A). Prędkość podano w µm/s i przedstawiono jako % całkowitej liczby trajektorii kulek (B). Rycina opracowana za zgodą autorów (31). Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 5. Schemat analizy segmentów trajektorii. Ta rycina została opracowana za zgodą autorów z (31). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6Badanie wpływu gęstości biofilmu na prędkość kulek przy użyciu potoku analizy biofilmu (krok 5). Wyniki wskazały, że prędkość kulek nie zależy wyłącznie od gęstości komórek. Biofilmy hodowano zgodnie z opisem w kroku 1. Trajektorie analizowano w skali pojedynczych segmentów trajektorii (Rysunek 5). Dla każdego segmentu prędkość kulek w µm/s wykreślono w zależności od gęstości komórkowej w obrębie ramki ograniczającej (średnia wartość GFP na woxel wewnątrz ramki ograniczającej). Czerwona linia przedstawia regresję liniową, a zielona linia przebieg regresji wykładniczej. Rysunek powielono za zgodą autorów z (31). Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej figury.
Wideo uzupełniające 1. Wideo 4D z 24-godzinnego biofilmu E. faecalis OG1RF hodowanego na szkiełku nakrywkowym z grubością 1,5. Wideo poklatkowe 4D zostało przygotowane przy użyciu mikroskopu o rozdzielczości 512x512. Seria Z składająca się z 40 obrazów została wygenerowana poprzez obrazowanie obszaru biofilmu o grubości około 20 µm w krokach 0,5 µm. Każda seria Z stanowiła jedną klatkę i wymagała 50-60 s na przechwycenie. Przechwycono serię 20 następujących po sobie klatek, aby stworzyć wideo 4D. Prędkość odtwarzania wideo wynosi około 120x. Wideo jest reprezentatywne dla co najmniej 6 niezależnych eksperymentów. Kliknij tutaj, aby pobrać to wideo.
Uzupełniający film 2. Film 4D z 24-godzinnego wzrostu biofilmu E. faecalis OG1RF na płytce 96-dołkowej z optycznym dnem. Film poklatkowy 4D został przygotowany przy użyciu mikroskopu w rozdzielczości 512x512. Seria Z składająca się z 40 obrazów została wygenerowana poprzez obrazowanie obszaru biofilmu o grubości około 20 µm w krokach 0,5 µm. Każda seria Z stanowiła jedną klatkę i wymagała 50-60 s na zarejestrowanie. W celu przygotowania filmu 4D zarejestrowano serię 20 następujących po sobie klatek. Prędkość odtwarzania filmu wynosi około 120x. Film jest reprezentatywny dla 3 niezależnych eksperymentów. Kliknij tutaj, aby pobrać ten film.
Materiał uzupełniający Wideo 3. Wideo 4D biofilmu Salmonella enterica serotyp Typhimurium ATCC 14028 hodowanego przez 6-7 dni na szklanych lamelkach optycznych o grubości 1,5. Wideo time-lapse 4D zostało przygotowane przy użyciu mikroskopu w rozdzielczości 512x512. Serię Z składającą się z 40 obrazów uzyskano poprzez obrazowanie obszaru biofilmu o grubości około 20 µm w krokach 0,5 µm. Każda seria Z stanowiła jedną klatkę i jej rejestracja trwała 50-60 s. W celu przygotowania wideo 4D zarejestrowano serię 20 następujących po sobie klatek. Prędkość odtwarzania wideo wynosi około 120x. Wideo jest reprezentatywne dla 3 niezależnych eksperymentów. Prosimy kliknąć tutaj, aby pobrać to wideo.
Uzupełniające wideo 4. Wideo 4D biofilów Salmonella enterica serotypu Typhimurium ATCC 14028 z mutacją w genach curli (csgBA), hodowanych na szkiełkach przedmiotowych o grubości 1,5 przez 6-7 dni. Wideo poklatkowe 4D zostało wykonane przy użyciu mikroskopu o rozdzielczości 512x512. Seria Z składająca się z 40 obrazów została utworzona poprzez obrazowanie obszaru biofilu o grubości około 20 µm w krokach co 0,5 µm. Każda seria Z stanowiła jedną klatkę i jej zarejestrowanie trwało 50-60 s. W celu uzyskania wideo 4D zarejestrowano serię 20 następujących po sobie klatek. Prędkość odtwarzania wideo wynosi około 120x. Wideo jest reprezentatywne dla 3 niezależnych eksperymentów. Kliknij tutaj, aby pobrać to wideo.
Dodatkowy film 5. Film 4D szczepu E. coli UTI89 hodowanego na szkiełkach dekoporycznych o grubości 1,5 przez 6–7 dni. Film poklatkowy 4D został wykonany przy użyciu mikroskopu w rozdzielczości 512x512. Serię Z składającą się z 40 obrazów uzyskano poprzez obrazowanie obszaru biofilmu o grubości około 20 µm w krokach co 0,5 µm. Każda seria Z stanowiła jedną klatkę i wymagała 50–60 s na zarejestrowanie. W celu przygotowania filmu 4D zarejestrowano serię 20 sąsiednich klatek. Prędkość odtwarzania filmu wynosi około 120x. Film jest reprezentatywny dla 3 niezależnych eksperymentów. Prosimy kliknąć tutaj, aby pobrać ten film.
Uzupełniający film 6. Film 4D mutantów curli (csgBA) szczepu E. coli UTI89 hodowanych na grubych szkiełkach optycznych o grubości 1,5 przez 6-7 dni. Film pokadrowy 4D został utworzony przy użyciu mikroskopu o rozdzielczości 512x512. Seria Z składająca się z 40 obrazów została uzyskana poprzez obrazowanie obszaru biofilmu o grubości około 20 µm w krokach 0,5 µm. Każda seria Z stanowiła jedną klatkę i wymagała 50-60 s na zarejestrowanie. Aby stworzyć film 4D, zarejestrowano serię 20 następujących po sobie klatek. Odtwarzanie filmu odbywa się z prędkością około 120x. Film jest reprezentatywny dla 3 niezależnych eksperymentów. Kliknij tutaj, aby pobrać ten film.