Klasyczna krzywa wzrostu hodowli okresowej ma fazę adaptacji, wzrost wykładniczy, fazę stacjonarną i zamieranie. Ten prosty schemat jest użyteczny, o ile wiadomo, co znajduje się na osi pionowej. Absorbancja, sucha masa, liczba cząstek, liczba komórek zdolnych utworzyć kolonię i ilość ATP mogą dać różne krzywe tej samej hodowli.

Krzywa nie pokazuje też losu „typowej komórki”. Jest wynikiem zbiorczym populacji, w której jedne komórki rosną, inne czekają na pierwszy podział, jeszcze inne umierają lub przechodzą w stan trudny do wyhodowania. Płaski odcinek oznacza brak zmiany mierzonego sygnału netto, a nie bezruch biologiczny.

Najpierw trzeba nazwać mierzoną wielkość

Słowo „wzrost” może oznaczać zwiększenie rozmiaru pojedynczej komórki, przyrost biomasy, wzrost liczby komórek, liczby jednostek tworzących kolonie albo stężenia produktu. Te procesy są powiązane, lecz nie zachodzą w stałych proporcjach w całym doświadczeniu.

Oś pionowa Co rośnie na wykresie Co może zafałszować interpretację
OD lub absorbancja pozorna rozpraszanie światła przez zawiesinę zmiana rozmiaru i kształtu, agregaty, pęcherzyki, osad, nieliniowość aparatu
liczba komórek w mikroskopie lub cytometrze cząstki spełniające kryterium zliczania martwe komórki, dublety, szczątki i trudne bramkowanie
JTK/CFU jednostki zdolne utworzyć kolonię w danych warunkach łańcuchy i agregaty, komórki VBNC, podłoże oraz czas inkubacji
sucha masa masa nielotnego materiału w próbce sole z podłoża, polimery zewnątrzkomórkowe i niedokładne oddzielenie supernatantu
ATP lub aktywność metaboliczna wybrany składnik albo reakcja związana z aktywnością różna zawartość na komórkę i metabolizm bez podziału
sygnał reportera ekspresja lub aktywność konkretnego układu dojrzewanie i trwałość reportera, obciążenie komórki

Pełna nazwa wyniku powinna brzmieć na przykład „OD₆₀₀ w funkcji czasu” albo „JTK/ml w funkcji czasu”, nie po prostu „krzywa liczby bakterii”. Szczegółowe właściwości każdej metody wymagają osobnego porównania; tutaj najważniejsze jest to, że metoda definiuje znaczenie faz.

Ta sama hodowla może rosnąć na jednym wykresie i spadać na drugim

W fazie zamierania liczba JTK może gwałtownie maleć, podczas gdy OD pozostaje wysoka, ponieważ martwe komórki nadal rozpraszają światło. Po przeniesieniu starej hodowli do świeżego podłoża komórki mogą zwiększać rozmiar i masę przed pierwszym podziałem: OD rośnie, ale liczba komórek prawie się nie zmienia.

Tworzenie łańcuchów może zmniejszyć JTK bez śmierci, jeżeli wiele komórek daje jedną kolonię. Rozpad łańcuchów może pozornie zwiększyć JTK bez podziału. Produkcja otoczki lub macierzy zwiększa masę i rozpraszanie, choć liczba genomów nie rośnie w tej samej proporcji.

Nie ma zatem jednej „prawdziwej” krzywej, a pozostałe nie są z definicji błędne. Każda mierzy inną warstwę procesu. Błąd zaczyna się wtedy, gdy wynik jednej metody otrzymuje nazwę wielkości, której nie mierzono.

Hodowla okresowa jest układem zmieniającym się w czasie

W hodowli okresowej, nazywanej też zamkniętą lub batch, na początku wprowadza się inokulum do skończonej porcji podłoża. W typowym układzie nie dopływa świeży substrat ani nie odpływa hodowla. Gaz, ciepło i próbki mogą być wymieniane, więc „zamknięta” nie znaczy hermetyczna.

Komórki zużywają substraty, zmieniają pH i potencjał redoks, wydzielają produkty, pobierają tlen i modyfikują własne otoczenie. Stężenie komórek nie jest jedyną zmienną. Każdy punkt czasowy odpowiada innemu podłożu, nawet jeśli naczynie pozostało to samo.

Klasyczne fazy są jakościowym opisem wyniku tych zmian. Granice nie są osobnymi stanami termodynamicznymi ani uniwersalnymi godzinami. W bogatym podłożu przejścia mogą się nakładać, a w złożonej społeczności kilka populacji może znajdować się w innych fazach jednocześnie.

„Warunki stałe” często oznaczają tylko stałe ustawienia urządzenia

Ustawiona temperatura inkubatora może być stała, lecz temperatura dużej, intensywnie rosnącej hodowli może się różnić. Stała szybkość obrotowa nie gwarantuje stałego natlenienia: wraz z biomasą rośnie zużycie tlenu, zmienia się lepkość i charakter mieszania.

Stałe pH początkowe nie oznacza kontrolowanego pH. Bufor ma skończoną pojemność. W kolbie bez regulacji zmiana kwasowości może zatrzymać wzrost wcześniej niż wyczerpanie głównego źródła węgla.

Dlatego nazwę fazy warto łączyć z pomiarem podłoża: glukozy, tlenu rozpuszczonego, pH, produktu lub właściwego składnika ograniczającego. Sam kształt OD nie wskazuje, który czynnik zakończył wzrost.

Bilans populacji oddziela przyrost od wyniku netto

Jeżeli N oznacza liczbę komórek, jej zmiana może być zapisana schematycznie:

dN/dt = podziały − śmierć − odpływ + dopływ.

W prostej hodowli okresowej dopływ i odpływ komórek są pomijalne, ale podziały oraz śmierć mogą zachodzić równocześnie. Stałe N nie mówi, czy oba strumienie są bliskie zeru, czy duże i prawie równe.

Dla biomasy X często zapisuje się dX/dt = μX, gdzie μ jest chwilową właściwą szybkością przyrostu biomasy. Równanie wykładnicze działa w przedziale, w którym μ można uznać za stałe i pomiar jest proporcjonalny do X. Poza tym przedziałem dopasowanie jednej prostej ma słabą interpretację.

Gdy osobno uwzględnia się utratę, można zapisać wynikową szybkość jako różnicę tempa wzrostu i tempa zaniku. Parametr wyznaczony z samego N(t) jest wtedy szybkością netto. Nie pozwala odzyskać obu strumieni bez dodatkowego pomiaru.

Wykładniczy wzrost jest lokalnym modelem, nie prawem całej hodowli

Przy stałym μ:

X(t) = X₀e^(μt).

Po zalogarytmowaniu:

ln X(t) = ln X₀ + μt.

Nachylenie prostej na wykresie ln X względem czasu jest więc μ. Dla dwóch punktów:

μ = [ln X₂ − ln X₁] / (t₂ − t₁).

Wybór punktów musi obejmować ten sam odcinek wykładniczy. Użycie pierwszego punktu po zaszczepieniu i ostatniego przed fazą stacjonarną daje średnią przez adaptację i spowolnienie, a nie maksymalną szybkość wzrostu.

Model nie wymaga, by każda komórka dzieliła się dokładnie w tym samym momencie. Wystarczy, że rozkład stanów populacji i średnie tempo są w badanym oknie stabilne. Asynchroniczna populacja może mieć bardzo gładki wzrost wykładniczy.

μ, czas podwojenia i czas generacji

Czas podwojenia biomasy w warunkach wykładniczych wynosi:

t_d = ln(2) / μ.

Jeżeli liczba komórek podwaja się przy każdym podziale i nie ma śmierci, czas podwojenia populacji odpowiada średniemu czasowi generacji. W rzeczywistej hodowli pojęcia mogą się rozchodzić. Komórki mogą rosnąć w masę bez podziału, tworzyć różnej wielkości potomstwo albo równocześnie umierać.

Czas generacji pojedynczej komórki jest odstępem między jej powstaniem a podziałem. Czas podwojenia populacji wynika z tempa wzrostu całego drzewa rodowego. Szybciej dzielące się linie są nadreprezentowane w próbce, dlatego prosta średnia czasów pojedynczych komórek nie zawsze równa się wynikowi populacyjnemu.

W zapisie z logarytmem dwójkowym liczba pokoleń między dwoma punktami wynosi:

n = log₂(N₂/N₁),

a średni czas na pokolenie w tym przedziale:

g = (t₂ − t₁)/n.

Wynik ma sens tylko wtedy, gdy N rzeczywiście reprezentuje mnożące się jednostki i przedział jest wykładniczy.

Ujemne μ nie jest „ujemnym wzrostem komórki”

Spadek liczby żywych jednostek daje ujemne nachylenie logarytmu. Może wynikać ze śmierci, utraty hodowalności, lizy lub agregacji. Nie oznacza, że każda komórka kurczy się wykładniczo.

Podobnie μ = 0 dla OD oznacza brak zmiany sygnału netto, nie brak syntezy. Interpretację mechanistyczną trzeba podeprzeć co najmniej jednym niezależnym pomiarem.

Faza lag jest przygotowaniem populacji, nie pustym czasem

Po przeniesieniu inokulum do nowego środowiska krzywa często ma odcinek o małym przyroście liczby. Komórki naprawiają uszkodzenia, odbudowują rybosomy, zmieniają transportery, dostosowują metabolizm, zwiększają rozmiar i dopiero potem dzielą się. Część może rozpocząć wzrost natychmiast, a inne po wielu godzinach.

Długość lag zależy nie tylko od nowego podłoża, lecz od historii inokulum. Komórki z fazy wykładniczej przeniesione do podobnego podłoża mogą prawie nie mieć laga. Stara kultura, nagła zmiana temperatury, inne źródło węgla, uszkodzenie kwasem lub wysuszenie zwiększają czas potrzebny do pierwszych podziałów.

Wielkość inokulum zmienia widoczny lag populacji. Im więcej komórek, tym większa szansa, że część szybko rozpocznie wzrost i przekroczy próg detekcji. W bardzo małym inokulum czas do sygnału obejmuje zarówno indywidualne opóźnienie, jak i liczbę podwojeń potrzebną do wykrycia.

Krzywizna populacji ukrywa rozkład pierwszych podziałów

Przejście od laga do wzrostu wykładniczego może powstać dlatego, że kolejne komórki włączają się do rosnącej frakcji. Dwie hodowle o podobnym średnim lagu mogą mieć wąski rozkład opóźnień albo małą szybko rosnącą frakcję i długi ogon.

Obserwacje pojedynczych komórek wykazały, że z samej gładkiej krzywej populacyjnej nie da się stabilnie odtworzyć pełnego rozkładu czasów laga.singlelag Ma to znaczenie przy ocenie rzadkich komórek po stresie: średni sygnał może wyglądać prawidłowo, choć część kolonii pojawia się znacznie później.

Lag obliczony przez przecięcie ekstrapolowanej prostej fazy wykładniczej z poziomem początkowym jest parametrem modelu. Nie musi równać się czasowi pierwszego podziału ani czasowi, w którym wszystkie komórki się zaadaptowały.

Faza przyspieszenia bywa ukryta w czterofazowym schemacie

Pomiędzy lagiem a stałym tempem wykładniczym μ stopniowo rośnie. Komórki odbudowują zdolność biosyntezy, a coraz większa część populacji wchodzi w cykl. Przy rzadkim próbkowaniu odcinek ten wygląda jak ostry narożnik.

Analogicznie przed fazą stacjonarną występuje spowolnienie. Substrat ograniczający może być jeszcze obecny, lecz tlen, pH, produkt hamujący albo zdolność transportu już obniża tempo. Dopasowanie prostej przez ten odcinek zaniża μmax.

Sześciofazowy opis — lag, przyspieszenie, wykładnicza, spowolnienie, stacjonarna i zamieranie — lepiej przypomina ciągłą dynamikę. Nadal jest opisem, nie obowiązkowym zestawem etapów każdej hodowli.

Faza wykładnicza nie zawsze oznacza wzrost z maksymalną możliwą szybkością

Odcinek wykładniczy oznacza w przybliżeniu stałą właściwą szybkość w danych warunkach. Może to być szybkość mniejsza od fizjologicznego maksimum, jeśli substrat, tlen lub inny czynnik już ogranicza komórki w sposób stały.

W bogatej hodowli okresowej maksymalne nachylenie może występować krótko. W minimalnym podłożu stabilny odcinek jest dłuższy, ale wolniejszy. Porównanie szczepów wymaga tego samego podłoża, temperatury, napowietrzania, fazy inokulum i sposobu wyznaczenia zakresu.

„Wzrost zrównoważony” oznacza, że średni skład komórek i rozkład ich stanów są w kolejnych pokoleniach w przybliżeniu stałe, a składniki komórki przyrastają proporcjonalnie. Nie jest synonimem dowolnego fragmentu o wysokiej OD. Po nagłym przesunięciu środowiska biomasa może przez chwilę rosnąć wykładniczo, podczas gdy skład RNA, białek i rozmiar komórek nadal się zmieniają.

Rozmiar komórek zmienia relację między OD a liczbą

W szybkim wzroście bakterie modelowe bywają większe i mają więcej rybosomów niż w głodzie. Po rozcieńczeniu nocnej kultury E. coli do świeżego podłoża szerokość i długość zmieniają się w innym tempie; maksimum rozmiaru może pokrywać się z maksimum szybkości zbiorczej.shape

Stały przelicznik „jedno OD = określona liczba komórek” może więc być poprawny tylko dla konkretnego szczepu, urządzenia i fazy. Kalibracja wykonana na fazie stacjonarnej nie musi pasować do komórek wykładniczych.

Faza stacjonarna jest bilansem, a nie zatrzymanym kadrem

Krzywa wypłaszcza się, gdy wynikowy przyrost mierzonej wielkości zbliża się do zera. Przyczyną może być wyczerpanie źródła węgla, azotu, fosforu lub mikroelementu, brak akceptora elektronów, spadek pH, toksyczny produkt, ograniczenie przestrzeni albo kilka czynników naraz.

Nie każda komórka przestaje się dzielić. Liza części populacji uwalnia składniki wykorzystywane przez inne komórki, co prowadzi do wzrostu ukrytego, zwanego czasem kryptycznym. Mutant lepiej wykorzystujący produkty lizy może zwiększyć udział bez zmiany całkowitej OD.

Populacja przebudowuje proteom, zmniejsza syntezę rybosomów, aktywuje naprawę i mechanizmy stresowe, zmienia powierzchnię oraz może wytwarzać metabolity wtórne, biofilm, formy przetrwalne lub persistery. Płaska linia zawiera zatem intensywną selekcję i różnicowanie.

Pojemność środowiska nie jest stałą gatunku

Maksymalny sygnał zależy od stężenia i jakości substratu, plonu biomasy, tlenu, buforu, objętości, geometrii naczynia oraz produktów. Ten sam szczep osiągnie inną „pojemność” w płytce, kolbie i bioreaktorze.

Parametr K w modelu logistycznym jest dopasowanym poziomem nasycenia w danych warunkach, a nie uniwersalną liczbą komórek, którą środowisko „może pomieścić”. Zmiana mechanizmu ograniczenia w trakcie doświadczenia osłabia jego sens biologiczny.

Faza zamierania zależy od definicji śmierci

Wykres JTK często spada, gdy komórki tracą zdolność tworzenia kolonii. Wykres OD może pozostać płaski do chwili lizy. Barwnik integralności błony może klasyfikować inaczej niż posiew, a komórka uszkodzona może odrosnąć dopiero na podłożu regeneracyjnym.

Spadek JTK bywa w przybliżeniu wykładniczy tylko w ograniczonym przedziale. Później pozostaje długi ogon z persisterami, sporami, mutantami albo komórkami w chronionych agregatach. Jedna stała szybkości śmierci nie musi opisywać całej populacji.

W bardzo starych hodowlach powstaje długotrwała faza stacjonarna. Populacja nie jest niezmienną resztką: warianty o przewadze w wykorzystaniu produktów innych komórek pojawiają się i zastępują poprzednie. Ewolucja zachodzi na tle ciągłego recyklingu biomasy.

Liza może chwilowo zwiększyć dostępność substratów

Śmierć jednych komórek tworzy nisze dla innych. DNA, aminokwasy i lipidy uwolnione do podłoża mogą być przyswajane. Jeżeli mierzy się tylko całkowity węgiel organiczny lub OD, nie widać przepływu między martwą i żywą frakcją.

Bilans izotopowy, pomiar metabolitów i śledzenie linii pozwalają rozdzielić pierwotny substrat od nekromasy. Jest to ważne w środowisku, gdzie klasyczna faza śmierci jednej populacji może być fazą wzrostu innej.

Diauksja tworzy dwie fazy wzrostu

Gdy podłoże zawiera dwa użyteczne źródła węgla, komórki mogą preferencyjnie zużywać jedno, a dopiero potem uruchomić szlaki dla drugiego. Krzywa ma pierwszy odcinek wzrostu, przejściowe spowolnienie i drugi odcinek. Klasycznym przykładem jest regulacja wykorzystania glukozy i laktozy przez E. coli.

Przerwa diauksyjna nie oznacza, że wszystkie komórki równocześnie całkowicie przestają rosnąć. Śledzenie pojedynczych komórek podczas zmiany z glukozy na laktozę ujawniło zróżnicowane trajektorie: część mocno zwalniała, część utrzymywała mierzalny wzrost, a średnia kolonii tworzyła klasyczny lag.diauxie

W innych układach podpopulacje mogą specjalizować się w różnych substratach albo produkty jednej frakcji żywić drugą. Dwie fazy na wykresie są wskazówką zmiany metabolizmu, nie dowodem konkretnego operonu bez analizy substratów i ekspresji.

Wyczerpanie jednego składnika może ujawnić ukrytą reakcję

Zmiana pH, tlenu lub produktu może wyprzedzać spadek OD. Próbkowanie metabolitów pokazuje, czy przerwa pokrywa się z wyczerpaniem pierwszego substratu i początkiem zużycia drugiego. Bez tego podobny kształt może wynikać z chwilowej zmiany temperatury, awarii mieszania lub agregacji.

Nagłe przesunięcie warunków nie jest nową krzywą od zera

Po dodaniu substratu, zmianie temperatury albo natlenienia populacja ma pamięć poprzedniego stanu: białka, metabolity, uszkodzenia i rozkład rozmiaru. Odpowiedź na przesunięcie w górę do bogatszego podłoża różni się od odpowiedzi na zaszczepienie świeżego podłoża komórkami z głodu.

Przy nutrient upshift synteza RNA i wzrost objętości mogą przyspieszyć przed zmianą częstości podziałów. Przy downshift istniejące zasoby pozwalają przez pewien czas kontynuować część procesów. Chwilowe μ wyznaczone z OD może nie reprezentować ustalonego stanu.

Porównanie przesunięć wymaga zdefiniowania punktu odniesienia: czasu dodania, czasu rzeczywistego wymieszania i opóźnienia pomiaru. W mikropłytce odparowanie podczas długiego eksperymentu może naśladować wzrost sygnału.

Populacja zbiorcza nie jest przeciętną trajektorią komórki

W rosnącym drzewie rodowym szybkie linie tworzą więcej potomków i silniej wpływają na wynik końcowy. Próbka pobrana w danym czasie jest obciążona w stronę komórek o krótszych generacjach. Średnia z losowo wybranych linii śledzonych od początku odpowiada innemu pytaniu.

Pojedyncza komórka zwiększa objętość i nagle dzieli się, podczas gdy całkowita biomasa populacji może rosnąć gładko. Faza stacjonarna całej hodowli może zawierać nisze o wzroście wykładniczym i nisze bez wzrostu.

Mikroprzepływowe „mother machine”, komory przepływowe i mikrokolonie na podłożu stałym zapewniają inne środowisko niż wymieszana kolba. Ułatwiają śledzenie linii, ale powierzchnia, stały dopływ i usuwanie potomstwa zmieniają ekologię. Wyniku pojedynczych komórek nie należy bezpośrednio nazywać krzywą hodowli okresowej.

Synchronizacja jest eksperymentalnym zaburzeniem

W asynchronicznej populacji fazy cyklu komórkowego są wymieszane. Synchronizacja przez selekcję rozmiaru, głodzenie lub zatrzymanie podziału pozwala obserwować okresowe zdarzenia, ale sama metoda może zmieniać fizjologię.

Oscylacja liczby lub rozmiaru po uwolnieniu synchronizacji stopniowo zanika, bo czasy generacji nie są identyczne. Brak oscylacji w zwykłej krzywej nie znaczy, że komórki nie mają uporządkowanego cyklu.

Natlenienie może określić kształt bardziej niż skład podłoża

Tlen rozpuszcza się słabo, a jego dopływ zależy od powierzchni gaz–ciecz, mieszania, lepkości i geometrii. W kolbie zbyt duża objętość robocza lub słabe wytrząsanie tworzą ograniczenie tlenu. W mikropłytce studzienki brzegowe i środkowe mogą różnić się odparowaniem oraz wymianą gazową.

Przejście do niedotlenienia zmienia oddychanie, wydajność biomasy i produkty. Krzywa może nadal rosnąć, lecz z innym μ i plonem. U organizmu względnie beztlenowego nie jest to po prostu „brak wzrostu”, tylko zmiana strategii energetycznej.

Pomiar tlenu w jednej lokalizacji dużego naczynia nie gwarantuje jednorodności. Biofilm, granula i agregat mają gradienty na skali mikrometrów, choć ciecz zbiorcza jest natleniona.

Krzywa w płytce różni się od krzywej w kolbie

Mała objętość ułatwia równoległe powtórzenia, ale zwiększa wpływ parowania, menisku, kondensacji i osiadania. Odczyt od dołu oraz od góry ma inne artefakty. Krótkie mieszanie przed pomiarem może nie rozbić kłaczków, a intensywne mieszanie może napowietrzać lub uszkadzać delikatne agregaty.

W kolbie pobieranie kolejnych próbek zmniejsza objętość i zmienia stosunek powierzchni do objętości. Osobne naczynia dla punktów czasowych eliminują ten efekt, lecz zwiększają zmienność między naczyniami.

Bioreaktor pozwala rejestrować pH, tlen, gazy wylotowe i dodatki regulatorów. Kontrola usuwa część zmian środowiska, ale dodanie zasady lub kwasu zwiększa objętość i zasolenie. Każdy układ ma własny bilans.

Hodowla ciągła nie przechodzi klasycznej sekwencji faz

W chemostacie świeże podłoże dopływa, a taka sama objętość hodowli odpływa. Szybkość rozcieńczenia D = F/V wynika z przepływu F i objętości V. W idealnym stanie ustalonym właściwa szybkość wzrostu populacji równa się szybkości rozcieńczenia: μ = D.

Nie oznacza to, że komórki rosną zawsze z μmax. Chemostat celowo utrzymuje ograniczający substrat na małym stężeniu, a operator ustala D. Gdy D przekracza zdolność populacji do wzrostu w danych warunkach, następuje wypłukanie.

Turbidostat steruje dopływem na podstawie zmętnienia, utrzymując zadany sygnał biomasy. Często pozwala rosnąć blisko szybkości wynikającej z podłoża. Nazwy nie są zamienne: chemostat kontroluje przepływ i ograniczenie składnikiem, turbidostat używa sprzężenia zwrotnego z mętnością.

Stan ustalony nie oznacza, że każda komórka jest taka sama

Stałe średnie stężenia mogą współistnieć z mutacjami, oscylacjami, przywieraniem do ścian i podpopulacjami. Biofilm w przewodzie lub reaktorze nie podlega temu samemu odpływowi co komórki planktoniczne i może zasilać ciecz potomkami.

Przed uznaniem układu za ustalony trzeba obserwować kilka czasów retencji oraz bilans substratu i produktu. Płaska OD po zmianie D może być przejściem maskowanym przez opóźnienie aparatu.

Modele logistyczny, Gompertza i Baranyiego opisują kształt na różne sposoby

Model logistyczny tworzy symetryczną krzywą esowatą z poziomem nasycenia. Zmodyfikowany model Gompertza elastyczniej opisuje asymetrię i jest często używany do laga, maksymalnego tempa oraz poziomu końcowego. Model Baranyiego wprowadza stan fizjologicznego przygotowania populacji.

Dobre dopasowanie nie potwierdza mechanizmu. Kilka modeli może przechodzić przez te same punkty, a parametry o nazwach „lag” i „maksimum” mogą mieć inne definicje. Ekstrapolacja poza zakres temperatury, pH lub inokulum jest ryzykowna.

Modele pierwotne opisują zmianę w czasie w jednym środowisku. Modele wtórne wiążą parametry z temperaturą, pH lub aktywnością wody. Błąd obu poziomów oraz zmienność biologiczna powinny zostać przeniesione do prognozy.

R² nie wystarcza do wyboru modelu

Punkty czasowe są zależne, wariancja rośnie z sygnałem, a reszty mogą mieć strukturę. Trzeba obejrzeć reszty, niepewność parametrów, korelacje i walidację na niezależnej serii. Model z większą liczbą parametrów może poprawić dopasowanie bez lepszej prognozy.

Logarytmowanie stabilizuje względny błąd, ale nie można logarytmować zera bez arbitralnej poprawki. Punkty poniżej granicy oznaczalności są cenzurowane, nie równe zeru.

Jak zaprojektować wiarygodną krzywą

Najpierw formułuje się pytanie. Do μmax potrzebne są gęste punkty w fazie wykładniczej. Do długości laga — bardzo wczesne punkty i dobrze zdefiniowana historia inokulum. Do śmierci — metoda żywotności, odpowiedni odzysk i długi ogon czasowy.

Następnie ustala się:

  • szczep, numer pasażu i sposób przygotowania inokulum;
  • podłoże z numerem partii, pH i składnik ograniczający;
  • objętość, naczynie, temperatura, mieszanie i natlenienie;
  • wielkość inokulum oraz sposób jej wyrażenia;
  • próbę ślepą, kontrole dodatnie i kontrolę sterylności;
  • częstotliwość pomiaru i wpływ pobierania próbek;
  • zakres liniowy aparatu i regułę rozcieńczania;
  • co najmniej niezależne powtórzenia biologiczne;
  • kryterium wyboru odcinka i modelu przed obejrzeniem wyniku.

Powtórzenia techniczne oceniają odczyt tej samej próbki. Powtórzenia biologiczne zaczynają się od niezależnych hodowli i obejmują zmienność przygotowania. Dziesięć odczytów jednej studzienki nie zastępuje trzech niezależnych inokulów.

Częstotliwość próbkowania powinna odpowiadać dynamice

Jeżeli czas podwojenia wynosi kilkadziesiąt minut, punkt co kilka godzin nie wyznaczy wiarygodnego nachylenia. Zbyt częste otwieranie kolby zmienia temperaturę i napowietrzenie. Automatyczny odczyt rozwiązuje część problemu, ale częste światło może wpływać na organizmy fotosyntetyczne lub reportery.

Wstępna krzywa szerokiego zakresu pozwala zaplanować właściwe okna. W kolejnej serii można gęsto próbkować przejścia, zamiast równomiernie zużywać pomiary w mało informacyjnej części.

Jak wyznaczać fazę wykładniczą bez wybierania „ładnej prostej”

Warto wykreślić logarytm sygnału po korekcie tła i obliczyć lokalne nachylenie w przesuwającym się oknie. Stabilny plateau μ wskazuje zakres wykładniczy. Okno musi zawierać dość punktów, by szum nie dominował, ale nie tak wiele, by mieszać przejścia.

Zakres liniowy OD należy sprawdzić osobno. Spłaszczenie przez nieliniowość detektora może wyglądać jak biologiczne spowolnienie. Rozcieńczenie próbki i przeliczenie do wartości pierwotnej jest wiarygodne tylko wtedy, gdy rozcieńczalnik nie tworzy dodatkowego tła i zawiesina pozostaje jednorodna.

Regułę wyboru punktów dobrze ustalić przed porównaniem szczepów. Ręczne wybieranie innych okien dla każdej krzywej zwiększa ryzyko dopasowania oczekiwanego wyniku.

Niepewność μ pochodzi z więcej niż regresji

Błąd nachylenia jednej prostej opisuje rozrzut punktów tej hodowli. Nie obejmuje różnic między dniami, inokulami i partiami podłoża. Wniosek biologiczny powinien opierać się na rozkładzie μ z niezależnych hodowli.

Jeżeli model jest dopasowany do wszystkich powtórzeń wspólnie, hierarchiczna struktura danych powinna zostać zachowana. Traktowanie każdego punktu czasowego jako niezależnego powtórzenia sztucznie zwiększa liczebność próby.

Minimalny zestaw do interpretacji mechanizmu

Krzywa sygnału biomasy mówi, kiedy zmienia się wynik zbiorczy. Aby odpowiedzieć dlaczego, warto równolegle mierzyć przynajmniej:

  1. wskaźnik liczby lub żywotności, niezależny od rozpraszania;
  2. substrat przypuszczalnie ograniczający;
  3. istotny produkt i pH;
  4. tlen lub inny akceptor, jeśli ma znaczenie;
  5. morfologię oraz agregację w wybranych punktach.

W badaniu diauksji dodaje się oba substraty. W badaniu śmierci — integralność, lizę i warunki odzysku. W badaniu różnic szczepów — kontrolę komplementacji lub niezależne allele, jeśli przyczyna ma być genetyczna.

Wielowymiarowy eksperyment nie musi mierzyć wszystkiego w każdej minucie. Kluczowe są punkty przed i po zmianie nachylenia oraz próbki zachowane do analizy metabolitów.

Szybkość wzrostu i plon biomasy odpowiadają na inne pytania

Szczep może szybko zużyć substrat i osiągnąć mało biomasy albo rosnąć wolniej, lecz wykorzystać większą część węgla do budowy komórek. Maksymalne μ opisuje tempo w sprzyjającym przedziale. Plon Y_X/S opisuje ilość biomasy X powstałą na jednostkę zużytego substratu S:

Y_X/S = ΔX / (−ΔS).

Minus przypomina, że stężenie substratu maleje. W praktyce plon obserwowany obejmuje utrzymanie komórki, produkty uboczne, straty gazowe i zmianę innych zasobów. Nie jest stałą absolutną, jeżeli zmienia się metabolizm lub ograniczenie.

W oddychaniu z dobrym natlenieniem większa część energii może wspierać biomasę niż w fermentacji, lecz szczegół zależy od organizmu i produktu. Dwie krzywe o takim samym nachyleniu OD mogą mieć inny pobór tlenu, wydzielanie kwasów i końcową masę.

Wysokie μ może obniżać wydajność procesu

W biotechnologii celem nie zawsze jest najszybszy przyrost komórek. Wysokie tempo może sprzyjać produktom ubocznym, przepełnieniu metabolizmu, stresowi tlenowemu albo niestabilności białka rekombinowanego. Optimum produktu może leżeć poniżej μmax.

Produktywność objętościowa łączy ilość produktu z czasem i objętością. Wydajność produktu względem substratu mówi o ekonomii węgla, a swoista produktywność — o pracy jednostki biomasy. Wysoka końcowa OD nie wystarcza do wyboru procesu.

Krzywa populacji powinna być więc zestawiona z bilansem. Bez stężenia substratu nie wiadomo, czy niższy poziom końcowy wynika z gorszego plonu, niepełnego zużycia, toksyczności czy wcześniejszego ograniczenia innym składnikiem.

Czas do progu detekcji nie jest bezpośrednio długością laga

Automatyczne systemy często podają czas, w którym sygnał przekroczył ustalony próg. Obejmuje on indywidualną adaptację, liczbę początkowych żywych jednostek, szybkość wykładniczą i czułość urządzenia. Ten sam lag przy dziesięciokrotnie mniejszym inokulum da późniejszy próg o około liczbę podwojeń potrzebną do odrobienia różnicy.

W idealnym modelu wykładniczym czas do progu X_T można zapisać:

t_T = [ln(X_T) − ln(X₀)] / μ.

Jeżeli istnieje lag λ, przybliżenie staje się t_T ≈ λ + [ln(X_T) − ln(X₀)]/μ. Jedna obserwacja t_T nie pozwala oddzielić λ, X₀ i μ. Potrzebna jest znana dawka początkowa albo pełna krzywa.

Brak wykrycia jest informacją przedziałową

Jeżeli sygnał nie przekroczył progu do końca eksperymentu, nie wiadomo, czy nie było zdolnych do wzrostu komórek, lag był dłuższy, czy tempo było zbyt małe. Wynik jest prawostronnie cenzurowany w czasie. Przypisanie mu wartości równej końcowi pomiaru zaniża opóźnienie.

Analogicznie początkowe punkty poniżej progu aparatu nie są zerami biomasy. Zastąpienie ich zerem lub połową limitu może obciążyć dopasowanie. Metoda statystyczna powinna jawnie uwzględniać granicę detekcji.

W testach niskiego inokulum duże znaczenie ma rozkład Poissona. Nie każda studzienka otrzymuje jedną komórkę tylko dlatego, że średnia wynosi jeden obiekt na objętość. Część jest pusta, część ma dwie lub więcej jednostek, a agregat łamie założenie niezależnych cząstek.

Krzywa może zmieniać skład genetyczny populacji

Kilka–kilkanaście pokoleń wystarcza, aby rzadki wariant o przewadze zwiększył udział. W hodowli okresowej przewaga może dotyczyć krótszego laga, większego μ, lepszego przeżycia fazy stacjonarnej albo wykorzystania produktu lizy. Są to różne składowe dostosowania.

Mutant szybciej rosnący wykładniczo może przegrywać w cyklach częstego głodu, jeśli ma długi lag lub małą przeżywalność. Z kolei wariant dominujący w starej hodowli nie musi rosnąć szybciej na świeżym podłożu; może lepiej używać nekromasy albo hamować konkurentów.

Powtarzane pasażowanie wybiera cechy zależnie od momentu transferu i rozcieńczenia. Transfer w połowie fazy wykładniczej selekcjonuje inny zestaw niż przenoszenie po długiej fazie stacjonarnej. „Adaptacja do podłoża” wymaga opisania pełnego cyklu selekcyjnego.

Plazmid i konstrukcja reportera mogą zmieniać krzywą

Utrzymanie plazmidu, ekspresja białka i antybiotyk selekcyjny obciążają komórkę. Populacja może tracić plazmid, mutować konstrukcję albo różnicować poziom ekspresji. Wzrost OD może wtedy pochodzić głównie od komórek nieprodukujących oczekiwanego białka.

Porównanie wymaga szczepu z pustym wektorem, kontroli stabilności i oznaczenia odsetka komórek zachowujących konstrukcję. Sama zgodność początkowego inokulum nie gwarantuje zgodności składu po kilkunastu pokoleniach.

Antybiotyk może degradować się albo być inaktywowany przez populację, osłabiając selekcję w późnej fazie. Zwiększenie dawki nie jest neutralnym rozwiązaniem, bo zmienia fizjologię i lag.

Krzywa mieszanej społeczności sumuje różne populacje

W kokulturze wzrost jednego organizmu może pokrywać się ze spadkiem drugiego. Całkowita OD pozostaje płaska, choć skład zmienia się radykalnie. Sekwencjonowanie względnych udziałów również nie daje bezwzględnej liczebności: udział taksonu może wzrosnąć tylko dlatego, że inny zanikł.

Do bilansu potrzebne są miary taksonowo swoiste połączone z całkowitą biomasą lub standardem ilościowym. qPCR, cytometria z sondami, posiewy wybiórcze i znakowane szczepy mają różne ograniczenia. Współczynnik względny bez mianownika nie odtwarza krzywych poszczególnych partnerów.

Krzyżowe odżywianie tworzy przesunięte fazy

Pierwszy organizm może przekształcać substrat w produkt wykorzystywany przez drugi. Partner rośnie dopiero po nagromadzeniu metabolitu, więc krzywa zbiorcza ma ramię lub drugi odcinek. Nie jest to klasyczna diauksja jednej komórki, choć kształt może być podobny.

Producent i konsument mogą stabilizować pH albo usuwać toksyczny produkt. Izolowany szczep osiąga wtedy inną fazę stacjonarną niż w społeczności. Mechanizm rozstrzyga się przez osobne krzywe, metabolity i doświadczenia z supernatantem lub kontrolowaną kokulturą.

W biofilmie odłączenie fragmentu zwiększa liczbę komórek planktonicznych bez miejscowego podziału w cieczy. Pobranie tylko supernatantu nie opisuje całkowitej populacji naczynia.

Pole pod krzywą nie ma jednej biologicznej definicji

AUC, pole pod krzywą, bywa używane do szybkiego porównania warunków. Łączy lag, tempo i poziom końcowy w jedną liczbę. Dwie hodowle mogą mieć to samo AUC: jedna rośnie szybko i wcześnie zamiera, druga wolno osiąga umiarkowany poziom.

Wartość zależy od czasu zakończenia i odstępów próbkowania. Dodanie długiego plateau zwiększa pole bez nowego przyrostu. AUC z surowej OD zawiera tło podłoża, a AUC po odjęciu tła może być ujemne przy szumie.

Jeżeli pytanie dotyczy mechanizmu, lepiej raportować osobno lag, μ, maksimum, czas przejścia i wskaźnik śmierci. AUC jest przydatnym wynikiem przesiewowym po zdefiniowaniu okna, lecz nie zastępuje parametrów.

Normalizacja do punktu początkowego może ukryć różne inokula

Podzielenie każdej krzywej przez jej wartość t=0 ułatwia porównanie względne, ale usuwa informację o dawce początkowej. Przy nieliniowym OD mały błąd tła w mianowniku silnie zmienia wynik.

Odejmowanie pierwszego punktu zakłada, że całe początkowe rozpraszanie jest niezmiennym tłem. W rzeczywistości pochodzi z komórek, które zmieniają rozmiar lub ulegają lizie. Sposób normalizacji powinien odpowiadać modelowi, nie tylko poprawiać wygląd wykresu.

Najczęstsze błędy interpretacyjne

  1. „OD to liczba żywych bakterii.” OD mierzy rozpraszanie przez cząstki i może pozostawać wysokie po śmierci.
  2. „Lag oznacza brak aktywności.” Komórki mogą rosnąć, naprawiać się i przebudowywać bez zwiększenia liczby.
  3. „Każda komórka wchodzi w fazę log jednocześnie.” Gładka krzywa może ukrywać szeroki rozkład laga.
  4. „Odcinek prosty na wykresie liniowym jest wykładniczy.” Wykładniczość ocenia się po logarytmie właściwej wielkości.
  5. „Najwyższa OD wyznacza pojemność gatunku.” Poziom zależy od podłoża, aparatu i geometrii.
  6. „Faza stacjonarna to brak podziałów.” Podziały i śmierć mogą się równoważyć.
  7. „Spadek CFU jest równy lizie.” Utrata hodowalności może poprzedzać rozpad komórki.
  8. „Czas podwojenia zawsze jest czasem generacji.” Zgodność wymaga określonych założeń o podziale i śmierci.
  9. „μ z dowolnych dwóch punktów opisuje szczep.” Punkty muszą pochodzić ze stabilnego odcinka i konkretnych warunków.
  10. „Diauksyjny lag zatrzymuje wszystkie komórki.” Średnia może powstać z różnych trajektorii.
  11. „Stałe obroty zapewniają stały tlen.” Zapotrzebowanie rośnie wraz z biomasą.
  12. „Dobre R² potwierdza mechanizm.” Dopasowanie opisuje kształt, nie przyczynę.
  13. „Powtórzenia studzienek są niezależnymi hodowlami.” Mogą być tylko powtórzeniami technicznymi.
  14. „Brak sygnału równa się zero komórek.” Istnieje granica detekcji i czas do progu.
  15. „Jedna krzywa wystarczy do porównania szczepów.” Potrzebna jest zmienność biologiczna i jawna reguła analizy.

Praktyczny tok czytania gotowego wykresu

Najpierw odczytuje się dokładną oś, jednostkę, skalę i korektę tła. Następnie sprawdza się metodę inokulacji, naczynie, warunki i liczbę niezależnych hodowli. Bez tych danych nie warto przypisywać fazom mechanizmu.

Kolejny krok to zaznaczenie przedziałów o zmiennym nachyleniu, a nie mechaniczne dopisanie czterech etykiet. Dwa plateau, brak śmierci lub kilka załamań mogą być prawidłowym wynikiem. Potem zestawia się krzywą z substratami, produktami, pH i żywotnością.

Na końcu rozdziela się obserwację od wniosku. „OD przestała rosnąć” jest obserwacją. „Wyczerpała się glukoza” jest hipotezą, dopóki glukozy lub bilansu nie zmierzono.

Granice tego artykułu

Szczegółowe przeliczenia OD, JTK, liczby komórek i suchej masy należą do osobnego artykułu o miarach biomasy. Tam zostaną omówione zakres liniowy, rozcieńczenia, granice detekcji i propagacja błędu.

Wpływ temperatury, pH, tlenu, zasolenia i ciśnienia na μ wymaga osobnego porównania czynników wzrostu. Chemostat, kinetyka Monoda i bilanse bioreaktora zostaną rozwinięte w kontekście hodowli ciągłej i biotechnologii.

Stan stacjonarny, persystencja i śmierć pojawiają się tutaj jako fazy krzywej, ale ich mechanizmy fizjologiczne wymagają osobnego tekstu.

Podsumowanie

Krzywa wzrostu jest wykresem określonego sygnału populacji w określonym układzie. Jej kształt zależy od tego, czy mierzono rozpraszanie, biomasę, cząstki, JTK czy metabolizm. Fazy nie są uniwersalnymi godzinami ani stanami każdej komórki.

Lag obejmuje przygotowanie i rozkład czasów pierwszego podziału. Odcinek wykładniczy pozwala wyznaczyć μ oraz czas podwojenia, jeżeli sygnał jest proporcjonalny do badanej wielkości i warunki są wystarczająco stabilne. Faza stacjonarna jest bilansem przyrostu i strat, a zamieranie zależy od definicji żywotności.

Diauksja, przesunięcia środowiska, ograniczenie tlenu i heterogeniczność pojedynczych komórek pokazują, dlaczego sam kształt nie wskazuje mechanizmu. Najlepsza interpretacja łączy krzywą z substratem, produktem, pH, tlenem, żywotnością i morfologią oraz jasno oddziela obserwację od hipotezy.

Opis hodowli okresowej, szybkości wzrostu i hodowli ciągłej oparto również na zakupionych podręcznikach mikrobiologii.baj