$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Pomiary krzywej wzrostu
Krzywe wzrostu S. acidocaldarius DSM639 pokazano na Rysunek 3A. Stwierdzono, że wzrost jest podobny, porównując inkubację przy użyciu termomieszalników z inkubatorami konwencjonalnymi. Parametry średniej szybkości wzrostu oszacowano, dopasowując krzywą logistyczną do każdej powtórzonej krzywej wzrostu i obliczając błąd średni i standardowy. Czasy do środkowej fazy wykładniczej na termomikserze i inkubatorze wynosiły odpowiednio 27,2 h ± 1,1 h i 31,1 h ± 1,9 h. Szacowane początkowe czasy podwajania dla termomiksera i inkubatora w temperaturze 75 °C wynosiły odpowiednio 4,29 h ± 0,28 h i 4,19 h ± 0,44 h, co jest zgodne z wcześniej publikowanymi wartościami24. Zależność między logarytmem10 (OD600 nm) a logarytmem10 (CFU) została dobrze scharakteryzowana przez model liniowy (skorygowany 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). OD600 nm 0,3 odpowiada zatem około 6,7 × 108 CFU/ml (Rysunek 3B).
Eksperyment z ewolucją temperatury
Trzy warunki temperaturowe, stała 75 °C, stała 65 °C i spadek temperatury (75-65 °C, zmniejszający się o 1 °C co dwa transfery), zostały zainicjowane przy użyciu siedmiu niezależnych linii wywodzących się od S. acidocaldarius DSM639. Pomiary OD600 nm zostały wykonane po każdym transferze w ciągu 45 dni eksperymentu (około 150 pokoleń w temperaturze 75 °C) i są pokazane na Rysunek 4. Pomiary OD600 nm wykonywane w ciągu dni są z natury hałaśliwe, ponieważ mogą podlegać subtelnym różnicom, na przykład w okresie wzrostu, temperaturze itp. Jednak pomiary wykonane w ciągu dni mogą być nadal przydatne do oceny żywotności populacji, a także dać wskazówkę, czy sprawność poprawia się z czasem. Linie od stałego warunku 75 °C wzrosły w OD600 nm z początkowego zakresu 0,125-0,3 do zakresu 0,248-0,471 pod koniec eksperymentu. Sugeruje to, że sprawność fizyczna poprawiła się dzięki temu zabiegowi. Natomiast linie ze stałego leczenia w temperaturze 65 °C wykazywały spadek OD600 nm, 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 przystosować się do stałej temperatury 65 °C, chociaż fakt, że organizmy zdolne do życia mogły zostać odzyskane, pokazuje, że populacje nie zostały wypłukane przez kolejne rozcieńczenia, co sugeruje pewien stopień adaptacji. Wreszcie, populacje w leczeniu spadkiem temperatury wzrosły z początkowego zakresu OD600 nm wynoszącego 0,099-0,279 do 0,3-0,39 w Tx6 (co odpowiada 288 godzinom i 73 °C w tym zabiegu), po czym nastąpił stały spadek do zakresu 0,003-0,024 w końcowym punkcie czasowym.
Testy wzrostu/sprawności
Testy dopasowania przeprowadzono dla każdej populacji potomków po eksperymentach ewolucyjnych. OD600 nm oznaczono po 48 godzinach wzrostu dla wszystkich siedmiu niezależnych linii, a następnie dopasowano modele liniowe w R dla każdej temperatury testu, z "środowiskiem selekcji" jako głównym efektem i "replikat/termomikser" jako efektem blokowym. Wzrost szczepu przodków został wykorzystany jako poziom odniesienia dla kontrastów leczenia. Dane są pokazane w Rysunek 5.
Gdy testowano w temperaturze 75 °C, zaobserwowano średnio znaczący wzrost dopasowania w stosunku do szczepu przodków dla linii ze stałej 75 °C (test t: t210=3,64, p=0,0003) i stałej 65 °C (test t: t210=2,8, p=0,005), ale nie z leczenia spadkiem temperatury (test t: t210=-0,87, p=0,38). Podczas oznaczania w temperaturze 65 °C, średnio, linie ze wszystkich metod leczenia wykazywały wzrost dopasowania (stałe linie o temperaturze 75 °C; test t: t210=4,68, p<0,0001, stałe linie 65 °C; Test t: t210,=4,24, p<0,0001, linie spadku temperatury; Test t: t210 = 3,15, p = 0,002). Jednak w przypadku obu temperatur testu wystąpiły znaczne różnice między liniami pod względem ich dopasowania (Ryc. 5). Niektóre linie nie różniły się znacząco od rodowodu lub miały słabsze przystosowanie; Było to szczególnie widoczne w przypadku linii pochodzących z leczenia spadkiem temperatury.
Warto zauważyć, że związek między OD600nm a CFU/mL mógł ulec zmianie podczas eksperymentu ewolucyjnego. Można to ocenić, określając parametry wzrostu dla wyewoluowanych linii (następujące kroki 3.1-3.10).
Wyniki sekwencjonowania całego genomu
Analizę całego genomu przeprowadzono przy użyciu breseq (wersja 0.38.1)25 na siedmiu liniach od stałego warunku 75 °C przy użyciu genomu referencyjnego S. acidocaldarius DSM639 (akces RefSeq NC_007181.1). Ujawniono różnorodne mutacje obejmujące insercje, delecje i polimorfizm pojedynczego nukleotydu (SNP) w genomach wszystkich linii potomnych (Tabela 2). Wielokrotne mutacje insercyjne i niesynonimiczne występowały w genach kodujących białka, a także w regionach międzygenowych, które mogą wpływać na ekspresję genów z powodu przesunięć ramki w regionach promotorowych genów. Niektóre z mutacji 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 i aktywność katalityczna (Tabela 2). Wśród tych mutacji była duża delecja 54 667 par zasad w pięciu z siedmiu linii; Zostało to potwierdzone przez wykres brakujących dowodów pokrycia dla każdej populacji (częstość waha się od 93,2% do 100%). Usunięty region jest równoznaczny z utratą 53 genów; Rola tych genów w adaptacji zostanie zbadana w przyszłych badaniach. Stwierdzono pewne różnice między zastosowanym izolatem DSM639 a opublikowaną sekwencją referencyjną (przedstawioną w tabeli uzupełniającej 1).
Zużycie energii przez inkubator z wytrząsaniem w porównaniu z termomikserem
Zużycie energii przez inkubator z wytrząsaniem porównano z termomikserem przy użyciu dostępnej na rynku inteligentnej wtyczki do monitorowania energii w zakresie typowych temperatur inkubacji i temperatury 75 °C stosowanej w tym przypadku. W temperaturze 75 °C termomikser zużywał około 1/40 energii tradycyjnego inkubatora z wytrząsaniem (Rysunek 6), co sugeruje termomiksery jako potencjalny sposób na zmniejszenie śladu węglowego związanego z ewolucją eksperymentalną.

Rysunek 1: Schemat blokowy ilustrujący protokół eksperymentu ewolucyjnego w trzech temperaturach. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 2: Schemat blokowy ilustrujący etapy protokołu testu wzrostu/sprawności. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 3: Określenie kluczowych parametrów wzrostu i porównanie urządzeń do inkubacji dla S. acidocaldarius DSM639. Trzy powtórzone kultury hodowano w temperaturze 75 °C w 7 oddzielnych probówkach. Do pomiaru (A) OD600nm (oznacza ± błąd standardowy n=3 kontrprób technicznych; niektóre słupki błędów są mniejsze niż symbole na wykresie); krzywe reprezentują dopasowane logistyczne modele wzrostu i (B) jednostki tworzące kolonie (CFU); reprezentują dopasowane modele regresji liniowej logarytmicznej). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 4: Reprezentatywne wyniki dla gęstości optycznej (OD600nm) uzyskane podczas eksperymentów ewolucyjnych. Gęstości optyczne niezależnych linii genetycznych zmierzone podczas eksperymentów ewolucyjnych w trzech temperaturach (stała 65 °C, stała 75 °C, spadek 75 °C-65 °C) prowadzonych przez około 150 pokoleń. Krzywe przedstawiają wygładzenia lessu w czasie dla każdej niezależnej linii. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 5: Reprezentatywne wyniki dla testów wzrostu. Testy wzrostu dla niezależnych linii wywodzących się od S. acidocaldarius DSM639 po eksperymentach z ewolucją temperatury (stała 65 °C, stała 75 °C, spadek 75 °C-65 °C) w porównaniu ze szczepem przodka. Dla wszystkich linii wzrost oznaczano w temperaturze 65 °C i 75 °C. Kolorowe punkty pokazują średni błąd standardowy ± powtórzeń technicznych (pokazanych na szaro, n = 12 dla przodka i n = 3 dla każdej wyewoluowanej linii). Szary pasek oznacza średni ± błąd standardowy dopasowania przodków. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Zużycie energii przez tradycyjne inkubatory w porównaniu z urządzeniami termomiksującymi. Zużycie energii było rejestrowane za pomocą dostępnej na rynku inteligentnej wtyczki do monitorowania energii przez 2 godziny. Kliknij tutaj, aby zobaczyć większą wersję tego wykresu.
Tabela 1: Receptury pożywek i roztwory podstawowe wymagane do uprawy S. acidocaldarius. Wszystkie pożywki i roztwory podstawowe powinny być wykonane z podwójnie destylowanegoH2O(ddH2O), a następnie wysterylizowane przez autoklawowanie lub sterylizację filtrem przez filtr 0,22 μm, jak wskazano. Proszę zapoznać się z sekcją 1 protokołu w celu uzyskania szczegółowego opisu sposobu przygotowania wszystkich pożywek i roztworów podstawowych. Kliknij tutaj, aby pobrać tę tabelę.
Tabela 2: Reprezentatywne wyniki dla sekwencjonowania całego genomu linii potomnych. Mutacje znalezione w liniach potomnych wywodzących się od S. acidocaldarius DSM639 ze stałej obróbki w temperaturze 75 °C. Mutacje wskazują na zmiany w stosunku do bezpośredniego przodka linii, który posiada kilka zmian w stosunku do sekwencji referencyjnej dla S. acidocaldarius DSM639 (RefSeq NC_007181; mutacje u przodka przedstawiono w Tabeli Uzupełniającej 1). n wskazuje liczbę linii, w których stwierdzono mutację. → gen na "ramce odczytu do przodu"; ← gen w "ramce odczytu wstecznego". †Zmiany te są względne w stosunku do przesunięcia ramki ramki (A)10→11 obecnego w izolacie S. acidocaldarius DSM639 (Tabela uzupełniająca 1). Kliknij tutaj, aby pobrać tę tabelę.
Tabela uzupełniająca 1: Mutacje obecne w izolacie S. acidocaldarius DSM639 w stosunku do sekwencji referencyjnej (akces RefSeq NC_007181.1). → gen na "ramce odczytu do przodu"; ← gen w "ramce odczytu wstecznego". ‡SACI_RS04020 jest opisany jako pseudogen w NC_007181.1, ale obserwowana tutaj mutacja przesunięcia ramki Δ1 bp przypuszczalnie przywraca jego funkcję, tak jak w przypadku mutacji, koduje białko o 100% identyczności z genem odwrotnej gyrazy rgy (akcesium RefSeq WP_176586667.1). Kliknij tutaj, aby pobrać ten plik.