Pomiary krzywej wzrostu
Krzywe wzrostu dla S. acidocaldarius DSM639 przedstawiono na Rysunku 3A. Stwierdzono, że wzrost był podobny przy porównaniu inkubacji z użyciem termomikserów do inkubacji w konwencjonalnych inkubatorach. Średnie parametry szybkości wzrostu oszacowano poprzez dopasowanie krzywej logistycznej do każdej powtórzonej krzywej wzrostu i obliczenie średniej oraz błędu standardowego. Czas do osiągnięcia środkowej fazy wykładniczej w termomikserze i inkubatorze wynosił odpowiednio 27,2 h ± 1,1 h oraz 31,1 h ± 1,9 h. Szacowany początkowy czas podwojenia dla termomiksera i inkubatora w temperaturze 75 °C wynosił odpowiednio 4,29 h ± 0,28 h i 4,19 h ± 0,44 h, co jest zgodne z wcześniej opublikowanymi wartościami24. Zależność między log10(OD600nm) a log10(CFU) została dobrze opisana modelem liniowym (skorygowane R2 = 0,82, F(1,22) = 104, p < 0,00001, nachylenie = 1,73 ± 0,17, punkt przecięcia = 9,73 ± 0,14). Zależność między OD600nm a CFU jest zatem określona wzorem CFU = 10(1,73 × log10(OD600nm) + 9,73). Wartość OD600nm równa 0,3 odpowiada zatem około 6,7 × 108 CFU/mL (Rysunek 3B).
Eksperyment ewolucji temperaturowej
Trzy warunki temperaturowe: stałe 75 °C, stałe 65 °C oraz spadek temperatury (75-65 °C, obniżanie o 1 °C co dwa transfery) zostały zainicjowane z wykorzystaniem siedmiu niezależnych linii pochodzących od S. acidocaldarius DSM639. Pomiary OD600nm były wykonywane po każdym transferze w ciągu 45 dni eksperymentu (około 150 generacji przy 75 °C) i są przedstawione na Rysunku 4. Pomiary OD600nm wykonywane w odstępach dobowych są z natury szumne, ponieważ mogą podlegać subtelnym różnicom np. w okresie wzrostu, temperaturze itp. Niemniej jednak pomiary dobowe mogą być nadal przydatne do oceny żywotności populacji, a także dawać wskazówkę, czy przystosowanie (fitness) poprawia się w czasie. Linie z warunku stałego 75 °C wykazały wzrost OD600nm z początkowego zakresu 0,125-0,3 do zakresu 0,248-0,471 pod koniec eksperymentu. Sugeruje to, że w tym wariancie nastąpiła poprawa przystosowania. W przeciwieństwie do tego, linie z warunku stałego 65 °C wykazały spadek OD600nm z początkowego zakresu 0,018-0,087 do 0,008-0,04 w końcowym punkcie czasowym. Sugeruje to, że populacje nie były w stanie zaadaptować się do stałej temperatury 65 °C, choć fakt, że udało się odzyskać żywe organizmy, wskazuje, że populacje nie zostały całkowicie usunięte w wyniku sukcesywnych rozcieńczeń, co sugeruje pewien stopień adaptacji. Wreszcie, populacje w wariancie spadku temperatury wzrosły z początkowego zakresu OD600nm 0,099-0,279 do 0,3-0,39 w Tx6 (odpowiadającym 288 h i 73 °C w tym wariancie), po czym nastąpił stały spadek do zakresu 0,003-0,024 w końcowym punkcie czasowym.
Analizy wzrostu/przystosowania
Analizy przystosowania przeprowadzono dla każdej populacji potomnej po zakończeniu eksperymentów ewolucyjnych. OD600nm zmierzono po 48 h wzrostu dla wszystkich siedmiu niezależnych linii, a następnie dopasowano modele liniowe w programie R dla każdej temperatury analizy, przyjmując „środowisko selekcyjne” jako efekt główny oraz „powtórzenie/termomikser” jako efekt blokowy. Wzrost szczepu przodka posłużył jako poziom referencyjny dla kontrastów traktowań. Dane przedstawiono na Rysunku 5.
Podczas analizy w 75 °C odnotowano średnio istotne wzrosty przystosowania w stosunku do szczepu przodka dla linii z grup poddanych stałej temperaturze 75 °C (t-test: t210=3.64, p=0.0003) oraz stałej temperaturze 65 °C (t-test: t210=2.8, p=0.005), natomiast nie zaobserwowano ich dla linii z grupy poddanej obniżeniu temperatury (t-test: t210=−0.87, p=0.38). Podczas analizy w 65 °C średnio wszystkie linie z każdej z grup wykazały wzrost przystosowania (linie ze stałej temperatury 75 °C; t-test: t210=4.68, p<0.0001, linie ze stałej temperatury 65 °C; t-test: t210,=4.24, p<0.0001, linie z obniżenia temperatury; t-test: t210=3.15, p=0.002). Jednakże dla obu temperatur analizy wystąpiły znaczne różnice w przystosowaniu poszczególnych linii (Rycina 5). Niektóre linie nie różniły się znacząco od szczepu przodka lub wykazały spadek przystosowania; było to szczególnie widoczne w przypadku linii z grupy poddanej obniżeniu temperatury.
Warto zauważyć, że zależność między OD600nm a CFU/mL mogła ulec zmianie w trakcie eksperymentu ewolucyjnego. Można to ocenić poprzez wyznaczenie parametrów wzrostu dla linii ewoluujących (zgodnie z krokami 3.1-3.10).
Wyniki sekwencjonowania całego genomu
Analizę całego genomu przeprowadzono przy użyciu programu breseq (wersja 0.38.1)25 dla siedmiu linii pochodzących z warunków stałej temperatury 75 °C, wykorzystując genom referencyjny S. acidocaldarius DSM639 (numer dostępu RefSeq NC_007181.1). W genomach wszystkich linii potomnych ujawniono różnorodne mutacje, w tym insercje, delecje oraz polimorfizmy pojedynczego nukleotydu (SNP) (Tabela 2). Liczne insercje i mutacje niesynonimiczne wystąpiły w genach kodujących białka, a także w regionach międzygenowych, co może wpływać na ekspresję genów ze względu na przesunięcia ramki odczytu w regionach promotorowych genów. Niektóre mutacje były spójne w wielu liniach potomnych w genach zaangażowanych w różne funkcje, takie jak biosynteza ściany komórkowej, transkrypcja, metabolizm, transport komórkowy oraz aktywność katalityczna (Tabela 2). Wśród tych mutacji znajdowała się duża delecja o długości 54 667 par zasad w pięciu z siedmiu linii; potwierdzono to za pomocą wykresu brakującego pokrycia (coverage evidence plot) dla każdej populacji (częstość od 93,2% do 100%). Usunięty region odpowiada utracie 53 genów; rola tych genów w adaptacji zostanie zbadana w przyszłych studiach. Zauważono pewne różnice między użytym izolatem DSM639 a opublikowaną sekwencją referencyjną (przedstawione w Tabeli uzupełniającej 1).
Zużycie energii w inkubatorze wytrząsowym w porównaniu z termomikserem
Zużycie energii przez inkubator z wytrząsaniem porównano z zużyciem energii przez termomikser, wykorzystując dostępną w handlu inteligentną wtyczkę do monitorowania energii w zakresie powszechnie stosowanych temperatur inkubacji oraz 75 °C temperatura zastosowana w tym przypadku. W 75 °Ctermomikser zużył około 1/40th energia tradycyjnego inkubatora z wytrząsarką (Rycina 6), co sugeruje, że termomiksery mogą stanowić potencjalny sposób na zmniejszenie śladu węglowego związanego z ewolucją eksperymentalną.

Rycina 1: Schemat blokowy ilustrujący protokół eksperymentu ewolucyjnego w trzech wariantach temperatury. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 2: Schemat blokowy ilustrujący kolejne etapy protokołu badania wzrostu/sprawności (fitness assay). Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 3: Wyznaczenie kluczowych parametrów wzrostu i porównanie urządzeń inkubacyjnych dla S. acidocaldarius DSM639. Trzy powtórne kultury hodowano w temperaturze 75 °C w 7 oddzielnych probówkach. Do pomiarów zastosowano pobieranie próbek niszczących w celu oznaczenia (A) OD600nm (średnie ± błąd standardowy dla n=3 powtórzeń technicznych; niektóre słupki błędów są mniejsze niż symbole punktów; krzywe reprezentują dopasowane logistyczne modele wzrostu) oraz (B) jednostek tworzących kolonie (CFUs; linie reprezentują dopasowane log-log liniowe modele regresji). Aby wyświetlić powiększoną wersję tego rysunku, kliknij tutaj.

Rycina 4: Reprezentatywne wyniki gęstości optycznych (OD600nm) uzyskane podczas eksperymentów ewolucyjnych. Gęstości optyczne niezależnych linii mierzone podczas eksperymentów ewolucyjnych w trzech wariantach temperatury (stałe 65 °C, stałe 75 °C, spadek 75 °C-65 °C) przeprowadzonych przez około 150 pokoleń. Krzywe przedstawiają wygładzanie Loessa w czasie dla każdej niezależnej linii. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 5: Reprezentatywne wyniki testów wzrostu.Testy wzrostu dla niezależnych linii wywodzących się z S. acidocaldarius DSM639 po eksperymentach z ewolucją temperaturową (stałe 65 °C, stałe 75 °C, spadek 75 °C-65 °C) w porównaniu ze szczepem przodkiem. Dla wszystkich linii wzrost badano w 65 °C i 75 °C. Kolorowe punkty przedstawiają średnią ± błąd standardowy powtórzeń technicznych (zaznaczonych na szaro, n = 12 dla przodka i n = 3 dla każdej wyewoluowanej linii). Szary pasek oznacza średnią ± błąd standardowy przystosowania przodka. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6: Zużycie energii w tradycyjnym inkubatorze w porównaniu do urządzeń typu thermomixer. Zużycie energii rejestrowano przez 2 h przy użyciu dostępnej w handlu inteligentnej wtyczki do monitorowania energii. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.
Tabela 1: Składy pożywek i roztwory zapasowe wymagane do hodowli *S. acidocaldarius*. Wszystkie pożywki i roztwory zapasowe powinny być przygotowane w wodzie dwukrotnie destylowanej H2O (ddH2O)2O) a następnie sterylizowano w autoklawie lub poprzez filtrację przez 0.22 µm przefiltrować zgodnie z instrukcją. Szczegółowy opis przygotowania wszystkich pożywek i roztworów stokowych znajduje się w sekcji 1 protokołu. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Reprezentatywne wyniki sekwencjonowania całego genomu linii potomnych. Mutacje znalezione w liniach potomnych wywodzących się z S. acidocaldarius DSM639 poddanych stałej temperaturze 75 °C. Mutacje wskazują zmiany w stosunku do bezpośredniego przodka linii, który posiada kilka zmian względem sekwencji referencyjnej dla S. acidocaldarius DSM639 (RefSeq NC_007181; mutacje u przodka przedstawiono w Tabeli uzupełniającej 1). n oznacza liczbę linii, w których znaleziono daną mutację. → gen w „ramce odczytu bezpośredniego”; ← gen w „ramce odczytu odwrotnej”. †Zmiany te są określone względem przesunięcia ramki odczytu (A)10→11 występującego w izolacie S. acidocaldarius DSM639 (Tabela uzupełniająca 1). Prosimy kliknąć tutaj, aby pobrać tę tabelę.
Tabela uzupełniająca 1: Mutacje obecne w izolacie S. acidocaldarius DSM639 w stosunku do sekwencji referencyjnej (numer dostępu RefSeq NC_007181.1). → gen w „ramce odczytu bezpośredniego”; ← gen w „ramce odczytu odwrotnego”. ‡SACI_RS04020 jest adnotowany jako pseudogen w NC_007181.1, lecz zaobserwowana tutaj mutacja przesunięcia ramki odczytu Δ1 bp prawdopodobnie przywraca jego funkcję, ponieważ wraz z mutacją koduje on białko o 100% identyczności z genem odwrotnej gyrazy rgy (numer dostępu RefSeq WP_176586667.1). Kliknij tutaj, aby pobrać ten plik.