Poniższe reprezentatywne wyniki podkreślają konieczność wykonywania powtórzeń oraz obrazują różne techniki wizualizacji i analizy danych dla każdego z trzech zestawów próbek. Przedstawiają one wyniki potencjalnych schematów oceny danych, które są odpowiednie do udzielenia odpowiedzi na standardowe pytania badawcze dotyczące poszczególnych zestawów. Przedstawione techniki nie są jednak ograniczone wyłącznie do zestawów próbek, z którymi zostały zaprezentowane. Surowe dane z cytometrii przepływowej zostały udostępnione publicznie pod identyfikatorem FlowRepository FR-FCM-ZY46. Wcześniej przesłane dane ASC są dostępne pod identyfikatorem FR-FCM-ZYD7.
Replikaty biologiczne i techniczne pobrano, przygotowano i przeanalizowano zgodnie z opublikowanymi protokołami11. Replikaty biologiczne z PC i ASC pobrano z trzech równoległych kolb. Replikaty biologiczne BC pobrano jako trzy kolejne próbki z jednego miejsca w instalacji pilotażowej. Trzy alikwoty jednej próbki utrwalono, zabarwiono i przeanalizowano, aby zapewnić powtarzalność techniczną. Powtarzalność metod oceniono zgodnie z wcześniej opublikowanymi procedurami11. Do obliczenia odchylenia standardowego (%) dla każdego bramkowania wykorzystano szablon bramkowania dla wszystkich próbek z jednego zestawu próbek (Supplementary File 1-S6).
W przypadku PC maksymalne odchylenie standardowe względnej liczności wyniosło 3,44%, natomiast średnie odchylenie standardowe wyniosło 2,13% (Rysunek 1). Dla ASC i BC maksymalne odchylenia standardowe względnych liczności wyniosły odpowiednio 1,27% i 0,88%, podczas gdy średnie wartości odchyleń standardowych wyniosły 2,13% i 0,21% (Plik uzupełniający 1-S8).
Standardową procedurą podczas badania czystej kultury szczepu jest sporządzenie krzywej wzrostu w celu analizy fazy lag, czasu generacji i gęstości komórek w określonych warunkach. Połączyliśmy tę procedurę z przebiegiem analizy cytometrycznej, aby uzyskać głębsze zrozumienie systemu (Rycina 2). Wykresy FSC względem fluorescencji DAPI ujawniają stany cyklu komórkowego kultury w różnych punktach hodowli okresowej. Ustanowiono główny szablon bramkowania (plik uzupełniający 1-S6) w celu ilościowego określenia proporcji komórek z jednym (c1n), dwoma (c2n) i wieloma chromosomami (cXn, Rycina 2B). Przeważająca część zaszczepionych komórek posiadała tylko jeden chromosom (93,7% w 0 h). Zmieniło się to drastycznie podczas wzrostu wykładniczego, kiedy niemal wszystkie komórki zawierały więcej niż jeden chromosom (99,6%, 4 h), a ponad połowa populacji (53,1%) zawierała nawet więcej niż dwa chromosomy. Ilustruje to zdolność P. putida do replikacji chromosomów szybciej niż wynika to z czasu generacji.
Analiza cytometryczna z przepływem okazała się bardzo użyteczna w monitorowaniu ewolucji społeczności drobnoustrojów. Pozwala ona śledzić dynamikę społeczności znacznie dokładniej niż bardziej zasobochłonne metody molekularne5,10. Przedstawiliśmy tę dynamikę w filmie, uzyskując przegląd zmian w społeczności osadu czynnego w czasie. Każda jednosekundowa klatka pokazuje wykres fluorescencji FSC względem DAPI dla danego punktu poboru próbek (Supplementary File 2). Bardzo wyraźna zmiana między dniem 0 a dniem 4 była następowana przez ustalenie się rdzennej społeczności po 7. dniu. Dodatkowe podspołeczności pojawiły się w 21. dniu. Stworzyliśmy wzorzec bramy nadrzędnej (Supplementary File 1-S6), który umożliwił ocenę dynamiki drobnoustrojów na poziomie podspołeczności. Dominujące podspołeczności w różnych etapach eksperymentu zostały wyraźnie zidentyfikowane przy użyciu narzędzia CyBar (Rycina 3). Połączenie tego narzędzia z rozkładem częstości względnej obfitości podspołeczności pomogło w wyborze bram, które są interesujące (znacząca zmiana obfitości wywołana w kluczowym punkcie czasowym) i wykonalne (powyżej 5% względnej obfitości) do sortowania i dalszej analizy22.
Zastosowania bioreaktorów w skali przemysłowej mogą wiązać się z potencjalnymi niejednorodnościami przestrzennymi wynikającymi z ograniczeń w mieszaniu. Są one również często narażone na zmieniające się parametry operacyjne, takie jak zmiany w jakości substratów niesyntetycznych. BC reprezentuje taki system i został pobrany z różnych lokalizacji w dynamicznie sterowanym reaktorze przepływowym tłokowego. Wykresy fluorescencji FSC vs. DAPI (Rysunek 4) dla przykładowych punktów poboru próbek wykazują niewielką niejednorodność przestrzenną, lecz wyraźną niejednorodność czasową. BC został następnie zbadany przy użyciu zarówno szybkiej, zautomatyzowanej metody CHIC, jak i bardziej szczegółowego podejścia CyBar opartego na szablonie bramkowania głównego (master gate template).
Automatyczny, szybki i obiektywny charakter sprawia, że metoda CHIC jest szczególnie interesująca do kontroli bioprocesów w warunkach przemysłowych. Pozwala ona na porównanie surowych danych .fcs ze wszystkich dostępnych próbek. Dwa z tych porównań przedstawiono na Rysunku 5. Otrzymane wartości dyssimilarności potwierdzają wcześniejsze obserwacje dotyczące heterogeniczności przestrzennej i czasowej. Wykresy NMDS wygenerowane z macierzy dyssimilarności narzędzi CHIC (Rysunek 6) dodatkowo potwierdzają ten wynik i odpowiednio go wizualizują.
Dodatkowo obliczono względną liczebność 23 subpopulacji BC, wykorzystując główny szablon bramkowania (Supplementary File 1-S6). Przeprowadzono analizę korelacji 48 próbek z okresu jednego roku, aby zrozumieć zależności funkcjonalne w społeczności mikrobiologicznej. Może to pomóc w zidentyfikowaniu bramek interesujących w dalszych badaniach (4 Sortowanie komórek). Silne korelacje dodatnie lub ujemne z parametrami abiotycznymi, takimi jak miana produktów, mogą pomóc w zrozumieniu i optymalizacji ekosystemów oraz procesów biotechnologicznych. Obciążenie organiczne uwzględnione w tej korelacji (Rycina 7) stanowi przykład takiego parametru abiotycznego. G5, G4 i G10 wykazują silne korelacje dodatnie i mogą zawierać gatunki fermentacyjne, natomiast korelacja ujemna w G8 może wskazywać na archeony metanogenne.

Rycina 1: Powtarzalność analizy cytometrycznej czystej kultury. Wykresy FSC względem fluorescencji DAPI z kuleczkami kontrolnymi (A, B) oraz wykresy słupkowe względnej liczebności [%] trzech podspołeczności z słupkami błędów reprezentującymi ± jedno odchylenie standardowe (C, D). Względne liczebności określono przy użyciu szablonu bramkowania przedstawionego w Pliku Uzupełniającym 1-S6. Pobrano trzy powtórzenia biologiczne A, B i C z krzywej wzrostu w 4 h, a z próbki A przygotowano trzy powtórzenia techniczne P1, P2, P3, które zmierzono odpowiednio trzykrotnie: M1, M2, M3; przedstawiono je wraz z odpowiednimi odchyleniami standardowymi. W bramce komórkowej zarejestrowano 50 000 zdarzeń. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 2: Analiza cyklu komórkowego czystej kultury pobranej podczas eksperymentu z krzywą wzrostu w czasie 12 h. (A). Wykresy FSC względem fluorescencji DAPI dla próbek pobieranych co dwie godziny, wykazujące zmianę zawartości DNA podczas fazy wzrostu wykładniczego. Ta seria pomiarów nie obejmuje kulek kalibracyjnych (patrz Supplementary File 1-S3). (B). Krzywa wzrostu oparta na gęstości optycznej z paskami błędów reprezentującymi ± jedno odchylenie standardowe oraz względną liczebność w podspołecznościach w czasie. Prosimy kliknąć tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 3: Ewolucja struktury społeczności osadu czynnego w ciągu 24 dni. (A) Wykres flowCyBar z nałożonymi rozkładami częstości przedstawionymi w formie wykresów pudełkowych dla poszczególnych bram (gates) uporządkowanych według podobieństwa trendów, odpowiadający klucz kolorów (B) oraz analiza podobieństwa na wykresie NMDS z R <0,001 (C). Wykresy bazowe pokazano w sekwencji filmowej w pliku uzupełniającym 2. Względne obfitości określono przy użyciu szablonu bram w pliku uzupełniającym 1- S6. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.

Rysunek 4: Heterogeniczność przestrzenna i czasowa struktury społeczności mikrobowej w przemysłowym reaktorze biogazowym z przepływem tłokowym. (A) Porównanie struktury społeczności mikrobowej w czterech punktach poboru próbek I-IV w dniu 223. (B) Porównanie zależnych od czasu zmian struktury społeczności od dnia 0 do dnia 384 w punkcie II. Szum został usunięty w celu poprawy wizualizacji. W bramce komórkowej zmierzono 250 000 zdarzeń (plik uzupełniający 1-S6). Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 5: Porównanie obrazów histogramów cytometrycznych — CHIC. Zasada działania algorytmu CHIC została zilustrowana poprzez porównanie dwóch próbek reprezentujących odpowiednio heterogenność przestrzenną i czasową. (i) Obrazy CHIC są tworzone z plików .fcs po wybraniu dwóch parametrów do wyświetlenia (FSC vs. fluorescencja DAPI). Skala szarości (0–255) koduje liczbę zdarzeń w pikselu. (ii) Podzbiory CHIC są generowane po określeniu obszaru zawierającego komórki (tutaj: 200 < x < 4000, 200 < y < 2200). (iii) Obraz kombinacji XOR jest tworzony poprzez porównanie pikseli między podzbiorami. Brak różnicy odpowiada wartości 0 i jest wyświetlany jako kolor czarny. Maksymalna różnica odpowiada wartości 200 i jest wyświetlana jako kolor biały. Odpowiadająca wartość XOR jest obliczana poprzez zsumowanie wartości skali szarości wszystkich pikseli w obrazie XOR. (iv) Obraz kombinacji nałożenia (overlay) jest tworzony poprzez dodanie wartości skali szarości pikseli podzbiorów. Odpowiadająca wartość nałożenia jest obliczana poprzez zliczenie pikseli obrazu nałożenia, które nie są równe zero, a zatem zawierają informacje. (v) Wartości odmienności są obliczane poprzez podzielenie wartości XOR przez wartości nałożenia. Stanowią one podstawę wykresu NMDS w Rycini 6. Zarówno wartości odmienności, jak i wykres NMDS wykazują pomijalną heterogenność przestrzenną oraz wyraźną heterogenność czasową wspólnoty. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 6: Analiza podobieństwa próbek społeczności biogazu oparta na metodzie CHIC. Próbka z dnia 0 została wykluczona z tego wykresu ze względu na skrajnie odmienną strukturę społeczności, która zniekształcałaby analizę (Plik uzupełniający 1-S11). Wykres NMDS potwierdza założenie o niskiej heterogeniczności przestrzennej i wyższej heterogeniczności czasowej, sformułowane przy użyciu narzędzia CHIC i wyjaśnione na Rycini 5. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 7: Analiza korelacji w obrębie społeczności biogazu. Macierz została obliczona przy użyciu współczynnika korelacji rang Spearmana i opiera się na względnej liczebności 23 podspołeczności, które zostały określone za pomocą głównego szablonu bramkowania w Pliku Uzupełniającym 1-S6. Zostały one skorelowane z szybkością obciążenia organicznego (OLR) dla 48 próbek z okresu jednego roku. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 8: Analiza podobieństwa społeczności planktonicznych (P) i opartych na osadzie (S) w ściekach. Próbki pobrano z osadnika wstępnego (PCL), komory osadu czynnego (AS) i komory fermentacyjnej (DT) pełnoskalowej oczyszczalni ścieków (wykresy kropkowe w Pliku uzupełniającym 1-S10). Obfitości względne określono przy użyciu szablonu bramkowania pokazanego w Pliku uzupełniającym 1-S6. Obfitości tych podspołeczności przeanalizowano również za pomocą narzędzia flowCyBar w celu utworzenia wykresu CyBar (Plik uzupełniający 1-S10) oraz wykresu NMDS (R <0,01). Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 9: Stabilność fiksacji czystej kultury. Próbki z eksperymentu z krzywą wzrostu pobrane w 0 h przechowywano przez 28 dni w temperaturze -80 °C po fiksacji w 15% glicerolu. Przedstawiono wykresy FSC względem fluorescencji DAPI wraz z koralikami kontrolnymi. Serie pomiarowe nie zawierają koralików (Plik uzupełniający 1-S3). W bramce komórkowej zarejestrowano 50 000 zdarzeń (Plik uzupełniający 1-S6). Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
Plik uzupełniający 1: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Kliknij tutaj, aby pobrać ten plik.