Opisany tutaj protokół został wykorzystany do stworzenia map przystosowania dla wariantów genetycznych S. cerevisiae w dwuwymiarowym gradiencie stresu obejmującym temperaturę i zasolenie. W szczególności dla każdego szczepu wariantowego wykorzystano system eVOLVER do przeprowadzenia eksperymentów z konkurencją parzystą przeciwko fluorescencyjnie znakowanemu szczepowi referencyjnemu (szczep S. cerevisiae BY4741 z zintegrowanym konstytutywnym reporterem mNeonGreen) w 16 różnych kombinacjach tych stresów (Rysunek 1). Jako warianty wybrano dwa szczepy z nokautem, z których każdy był przewidywany jako wrażliwy na jedną z osi warunków stresowych: ΔPRX1, o którym zgłaszano defekty przystosowania w wysokich temperaturach20, oraz ΔPBS2 z defektami przystosowania przy wysokich stężeniach soli21. Jako wariant kontrolny wykorzystano nieznakowany szczep dziki BY4741, od którego oczekiwano wyników podobnych do szczepu referencyjnego. Łącznie te trzy 72-godzinne eksperymenty konkurencji przeprowadzono jednocześnie w 48 fiolkach w trzech urządzeniach eVOLVER (po 16 fiolek w każdym), sterowanych z jednego komputera.
Interaktywny panel sterowania eVOLVER służy do nawigacji po dużej ilości danych generowanych podczas każdego eksperymentu (Rysunek 2A). Reprezentatywne krzywe OD dla wszystkich 16 fiolek w jednym urządzeniu eVOLVER przedstawiono na Rysunku 2B. Każda krzywa wykazuje charakterystyczny wzór ząbkowania wynikający z algorytmu sterowania turbidostatem, który wykorzystuje sprzężenie zwrotne OD do inicjowania rozcieńczeń i utrzymywania kultur w wąskim zakresie gęstości (w tym przypadku 0,2-0,3). Dane OD i temperatury są mierzone w sposób ciągły, przesyłane do bazy danych w chmurze i aktualizowane w czasie rzeczywistym w interaktywnym panelu sterowania eVOLVER. Z danych przesyłanych w sposób ciągły można obliczyć wskaźniki wyższego rzędu (np. tempo wzrostu, skumulowaną liczbę pokoleń) i wykreślić je w panelu sterowania (Rysunek 2C), a nawet wykorzystać je jako parametry sprzężenia zwrotnego w różnych schematach selekcji (np. morbidostat).
Aby wygenerować mapy sprawności (fitness maps), mierzymy tempo, w jakim szczep referencyjny wypiera szczep wariantowy, łącząc pomiary cytometrii przepływowej z danymi o wzroście z systemu eVOLVER. W szczególności frakcje populacyjne oblicza się jako logarytm stosunku szczepu wariantowego do szczepu referencyjnego, wykorzystując dane z cytometrii przepływowej dla próbek z każdego punktu czasowego. Dla każdej hodowli i punktu czasowego liczba pokoleń jest interpolowana na podstawie danych o tempie wzrostu z eVOLVER w każdym okresie rozcieńczenia. Po naniesieniu logarytmów stosunków na liczbę pokoleń zaczynają pojawiać się jakościowe różnice między warunkami (Rysunek 3A). Aby obliczyć sprawność, przez liniową część każdego wykresu dopasowuje się linię prostą18,19, co pozwala uzyskać ilościową wartość sprawności (zarówno znak, jak i tempo), opisującą sposób konkurowania szczepu wariantowego ze szczepem referencyjnym (Rysunek 3B). W tym eksperymencie ujemne wartości sprawności wskazują, że szczep wariantowy jest wypierany przez szczep referencyjny. Konstrukcja map ciepła na podstawie wszystkich obliczeń sprawności umożliwia wizualizację dwuwymiarowego krajobrazu sprawności, ujawniając subtelne ilościowe różnice w wydajności między szczepami a warunkami (Rysunek 3C).
Z map ciepła dostosowania wynika, że każdy szczep wariantu wykazuje deficyty dostosowania, które manifestują się w różny sposób. ΔPRX1 jest przede wszystkim wrażliwy na wysokie temperatury, jednak stres solny wydaje się mieć efekt addytywny, zwiększając deficyt dostosowania. Z kolei ΔPBS2 wykazuje deficyty dostosowania głównie w odpowiedzi na wysokie stężenia soli, przy minimalnych różnicach wzdłuż osi temperatury. Wyniki te podkreślają użyteczność eksperymentów selekcyjnych o wysokiej rozdzielczości, szczególnie w celu rozszyfrowania interakcji wielu stresorów środowiskowych, które mogą mieć efekty addytywne, synergistyczne lub epistatyczne.
Wreszcie należy zauważyć, że w niektórych przypadkach, szczególnie w warunkach silnego stresu, obliczenia przystosowania mogą prowadzić do przesadzonych lub zniekształconych wyników. Na przykład ΔPBS2 wydaje się wykazywać stromy dodatni nachylenie w warunkach 39 °C/1 M NaCl w stosunku do szczepu referencyjnego (Rycina 3A). Przypisuje się to słabemu wzrostowi obu szczepów, co odzwierciedlają dane dla typu dzikiego i ogranicza użyteczność tego konkretnego parametru w tych warunkach. Niedostateczna liczba generacji skutkuje 1) słabym rozdzieleniem punktów czasowych, co zakłóca dopasowanie nachylenia, oraz 2) wrażliwością na efekty stochastyczne wynikające z fazy adaptacyjnej. Chociaż nie zdecydowaliśmy się na to w niniejszym badaniu, dane o wzroście w czasie rzeczywistym gromadzone przez eVOLVER podczas eksperymentu mogłyby zostać wykorzystane do podjęcia decyzji o przedłużeniu badania i zebraniu dodatkowych punktów czasowych dla tych warunków przy niewielkim nakładzie pracy. W kolejnych eksperymentach można by również zmodyfikować stosunki początkowe, aby wydłużyć obszar liniowy przed wystąpieniem nasycenia.

Rycina 1: Przegląd platformy eVOLVER i projektu eksperymentalnego. (A) eVOLVER to zautomatyzowana platforma do ciągłych kultur, umożliwiająca wieloparametrową kontrolę warunków hodowli. Platforma składa się z tzw. Smart Sleeves, w których znajdują się czujniki i elementy wykonawcze do kontroli warunków hodowli, zestawu pomp perystaltycznych do doprowadzania pożywki i odprowadzania odpadów z kultur, ekranu dotykowego do ręcznej interakcji z urządzeniem oraz komputera służącego do komunikacji z infrastrukturą chmurową, do której strumieniowane są dane z platformy. (B) eVOLVER można konfigurować na wiele sposobów: fizycznie, poprzez dopływ komórek i pożywek, oraz programowo, poprzez dostosowanie algorytmów sterujących modułami fiolek hodowlanych Smart Sleeve. Można to wykorzystać do programowania wielowymiarowych gradientów selekcyjnych/środowiskowych o wysokiej rozdzielczości, przy czym wynikami są fizycznie pobrane komórki w określonych punktach czasowych oraz strumieniowanie danych o hodowli w czasie rzeczywistym. (C) Konfiguracja eVOLVER w celu przeprowadzenia eksperymentu dopasowania wzrostu (growth fitness) w dwuwymiarowym gradiencie selekcyjnym, składającym się z temperatury i stresu osmotycznego. Fiolki hodowlane eVOLVER zostały zaszczepione równymi proporcjami szczepu drożdży referencyjnego i wariantowego. Zaprogramowano rutynę turbidostatu, aby utrzymać wszystkie kultury w fazie wykładniczej w określonym oknie OD (OD 0.2-0.3). Warunki temperatury i pożywek były zróżnicowane w obrębie zestawu Smart Sleeve, aby utworzyć dwa niezależne gradienty stresu. Pożywka podstawowa składała się z YPD (2% glukozy) + 100 µg/mL karbenicyliny + 25 µg/mL chloramfenikolu jako środka zapobiegającego kontaminacji bakteryjnej. Skład pożywek zmieniano poprzez dodanie suplementacyjnego NaCl do stężeń 0 M, 0.6 M, 0.8 M lub 1.0 M. Temperatury fiolek zaprogramowano na jedną z czterech wartości: 30 °C, 35 °C, 37 °C lub 39 °C. (D) Przykłady surowych wyników dla OD kultury i frakcji populacji z trzech testowanych warunków. Eksperyment konkurencji był prowadzony przez 72 h, z pobieraniem próbek w 24 h, a następnie co 12 h aż do zakończenia eksperymentu. Aby wyświetlić większą wersję tej ryciny, kliknij tutaj.

Rysunek 2: Panel analizy danych w czasie rzeczywistym eVOLVER. (A) Dane z eVOLVER są przetwarzane i przesyłane w czasie rzeczywistym za pomocą oprogramowania w chmurze do interaktywnego panelu sterowania. Za pomocą panelu użytkownicy mogą definiować i uruchamiać własne rutyny hodowli na podłączonym urządzeniu eVOLVER, monitorować aktualnie prowadzone eksperymenty we wszystkich jednostkach eVOLVER oraz ręcznie zmieniać warunki eksperymentalne w pojedynczej fiolce lub zestawie fiolek, jeśli jest to konieczne. (B) Krzywe OD dla każdej fiolki w macierzy badanych warunków przez cały czas trwania eksperymentu. Krzywe można obserwować w czasie rzeczywistym, aby upewnić się, że eksperyment przebiega zgodnie z planem, oraz wykorzystać je po zakończeniu eksperymentu do dalszej analizy. (C) Tempa wzrostu i liczba generacji komórek mogą być obliczane na bieżąco z krzywych OD w celu śledzenia postępu eksperymentu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3: Konstrukcja powierzchni dostosowania. (A) Logarytm naturalny frakcji populacji komórkowej (szczep wariantowy/szczep referencyjny) wykreślony w funkcji generacji komórkowych. Kolor wskazuje temperaturę selekcji, a kreski rozróżniają stężenia soli. (B) Aby ilościowo określić względne dostosowanie wariantu do szczepu referencyjnego, wskaźnik dostosowania jest wyprowadzany z szybkości, z jaką frakcja populacji komórkowej zmienia się w kolejnych generacjach. Wskaźnik ten jest obliczany z liniowego obszaru wykresów frakcji komórkowych przy użyciu minimum trzech punktów. (C) Względne wskaźniki dostosowania są przedstawione na mapie ciepła w celu skonstruowania powierzchni dostosowania w gradientach stresu środowiskowego. Warunek w prawym górnym rogu (39 °C/1,0 M NaCl) dla typu dzikiego i ΔPBS2 nie został uwzględniony w tej analizie, ponieważ kultury przeszły przez niewystarczającą liczbę generacji (patrz Reprezentatywne wyniki). Prosimy kliknąć tutaj, aby zobaczyć powiększoną wersję tej ryciny.