$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Protokół ręcznego QFA (dla szczepów S. cerevisiae)
1. Hodowla szczepów drożdży
- Do 96 niezależnych szczepów drożdży hoduje się w 200 μl bogatego medium płynnego w 96-dołkowej płytce hodowlanej. Szczepy namnaża się do stanu nasycenia w inkubatorze z kontrolowaną temperaturą.
- 96-igłowy manualny replikator o średnicy igieł 1/8" (SIGMA R-2508) sterylizuje się poprzez zanurzenie w 100% etanolu, opalenie nad płomieniem, a następnie pozostawienie do ostygnięcia.
- W 96-dołkowej płytce hodowlanej rozmieszcza się po 200 μl sterylnej wody. Nasycone kultury poddaje się mieszaniu w wortexie (Eppendorf MixMate, 30 sec, 750 rpm) i rozcieńcza, zanurzając wysterylizowane narzędzie igłowe trzy razy w nasyconych szczepach, a następnie raz w płytce zawierającej wodę.
- Narzędzie igłowe sterylizuje się ponownie zgodnie z punktem 1.2.
- Wysterylizowane narzędzie igłowe służy do nanoszenia rozcieńczonej kultury na płytki z agarem stałym poprzez zanurzenie w płytce z rozcieńczoną kulturą trzy razy, a następnie przeniesienie narzędzia na prostokątną płytkę z agarem stałym.
- Prostokątne płytki z agarem stałym inkubuje się w odpowiedniej temperaturze, a następnie obrazuje manualnie za pomocą urządzenia S&P Robotics spImager. Można zastosować alternatywne metody obrazowania, takie jak FujiFilm LAS4000 lub tańszy aparat cyfrowy lustrzanka, pod warunkiem zachowania równomiernego oświetlenia i spójnego ustawienia płytki względem aparatu.
- Metodę można dostosować do mniejszej liczby kultur równoległych, np. 48, które można umieścić na okrągłych szalkach Petriego, dbając o zachowanie tego samego kąta względem aparatu podczas każdego zdjęcia.
W pełni zautomatyzowany protokół QFA (dla szczepów S. cerevisiae)
2. Hodowla szczepów drożdży
- Rozpocznij od prostokątnej macierzy niezależnych szczepów drożdży rosnących na stałym agarze. W tym przykładzie szczepy wyjściowe są wynikiem skrzyżowania temperaturoczułej mutacji zapytania z Kolekcją Delecji Drożdży (YDC) metodą Synthetic Genetic Array (SGA). Końcowe płytki SGA są w formacie 1536. Jest to 16 wierszy i 24 kolumny kolonii/kultur reprezentujących cztery powtórzenia 384 niezależnych szczepów drożdży. Układy płytek przedstawiono na Rysunku 1.
- 384 szczepy spośród 1536 są przenoszone robotycznie przy użyciu robota S&P Robotics Inc BM3-SC z 96-pinowym, sterylnym narzędziem pinowym o średnicy pinu 1 mm do czterech płytek hodowlanych z 96 dołkami zawierających 200 μl pożywki selektywnej (Rysunek 1). Czystość i sterylność narzędzia pinowego uzyskuje się poprzez sekwencyjne mycie w wodzie sterylnej dwukrotnie (raz za pomocą obrotowej szczotki w celu usunięcia zanieczyszczeń), w 70% etanolu (z sonikacją) i na koniec w 100% etanolu.
- Kultury są hodowane do nasycenia (w tym przykładzie przez trzy dni w 20 °C, patrz Rysunek 1).
3. Wysiewanie kultur drożdży
- Przy użyciu urządzenia Beckman Biomek FX, nasycone kultury resuspendsuje się poprzez mieszanie w vortexie Variomag Teleshake (przeciwnie do ruchu wskazówek zegara przy 1000 rpm przez 20 s), a następnie rozcieńcza się je w stosunku około 1:70, przenosząc do 200 μl sterylnej wody za pomocą sterylnego robotycznego narzędzia typu pin tool do płytek 96-dołkowych (magnetyczny pin tool V&P Scientific, średnica pinu 2 mm). Czystość i sterylność narzędzia zapewnia się poprzez płukanie w sterylnej wodzie, płukanie w 70% etanolu (z użyciem szczoteczki w celu usunięcia zanieczyszczeń), a na koniec zanurzenie w 70% etanolu i osuszenie.
- Rozcieńczone kultury nanosi się na stałe pożywki agarowe w formacie 384 za pomocą tego samego robotycznego narzędzia pin tool (patrz Rysunek 1) po jego oczyszczeniu i sterylizacji.
- Po naniesieniu kultur, płytki przenosi się do inkubatora Cytomat z kontrolą temperatury i wilgotności, wyposażonego w automatyczny karuzelowy system transportowy i okno dostępu.
4. Rejestracja obrazu
- Zautomatyzowany system obrazowania S&P Robotics wielokrotnie wyjmuje płytki z inkubatora Cytomat, zdejmuje pokrywki, precyzyjnie umieszcza je w zamkniętej, równomiernie oświetlonej przestrzeni pod aparatem Canon EOS Rebel Ti 35mm DSLR, wykonuje zdjęcie w rozdzielczości 5184 x 3456 px, zakłada pokrywki i zwraca płytki do karuzeli inkubatora. Jedno zdjęcie jest wykonywane bezpośrednio po początkowym przeniesieniu do inkubatora, aby umożliwić pomiar wzrostu w punkcie zerowym (tło). Obsługa robotyczna i fotografia zajmują 2 min na płytkę.
- Przy pełnej karuzeli maksymalna osiągalna częstotliwość przechwytywania obrazów płytek wynosi jedno zdjęcie co 6 hr. Płytki są inkubowane do sześciu dni (do momentu nasycenia wzrostu kolonii). Rysunek 1 przedstawia trzy typowe obrazy z serii czasowej w 20 °C.
- W celach śledzenia płytki są nazywane zgodnie z nazwą inkubatora, temperaturą, unikalnym sekwencyjnym numerem partii oraz pozycją w inkubatorze. Nazwy obrazów stanowią połączenie nazwy płytki i znacznika czasu i są automatycznie tworzone podczas przechwytywania obrazu (np. K000011_030_010_2011-09-22_09-47-54.jpg, typ pliku FT1, Tabela 1).
5. Rejestracja metadanych eksperymentalnych
Pliki metadanych opisane w tej sekcji to pliki tekstowe rozdzielane tabulatorami, które są generowane ręcznie (np. przy użyciu arkuszy kalkulacyjnych). Plik Opisu Eksperymentu (FT2, Tabela 1) jest unikalny dla każdego eksperymentu, natomiast Opis Biblioteki (FT3, Tabela 1) może być wielokrotnie wykorzystywany po jego utworzeniu.
- Eksperyment został opisany w pliku z opisem eksperymentalnym (FT2, Tabela 1), który zawiera kolumny dotyczące kodu kreskowego (lub automatycznie wygenerowanej nazwy płytki), znacznika czasu rozpoczęcia eksperymentu, traktowania płytki, składu stałego podłoża agarowego, nazwy przesiewania, biblioteki, numeru płytki (dla bibliotek obejmujących wiele płytek) oraz numeru kwadrantu powtórzeń (patrz Rysunek 1) w celu przeskalowania z formatu 1536 na format 384.
- Biblioteka szczepów drożdży jest opisana w pliku z opisem biblioteki (FT3, Tabela 1), określającym genotyp hodowany w każdym miejscu kultury na każdej płytce w bibliotece. Zawiera on kolumny: nazwa biblioteki, ORF, numer płytki, wiersz płytki, kolumna płytki oraz opcjonalną kolumnę uwag.
- Można dostarczyć opcjonalny plik ze standardowymi nazwami genów (FT4, Tabela 1) opisujący standardową nazwę genu (np. RAD9) powiązaną z każdym systematycznym numerem y ORF (np. YDR217C) identyfikującym przesiewane szczepy. Plik ten zawiera dwie kolumny: ORF i nazwa genu.
6. Analiza danych
Obliczeniowy schemat pracy QFA wymaga dostępu do odpowiednio wydajnej wielordzeniowej stacji roboczej (np. Dell Precision T3500 z czterordzeniowym procesorem Xeon 2,67 GHz i 12Gb RAM), na której zainstalowane są narzędzie do analizy obrazów Colonyzer 3,4 oraz pakiet R dla QFA7; oba te narzędzia są udokumentowane online, ogólnodostępne i działają na wielu systemach operacyjnych.
- Uruchom program Colonyzer, wykorzystując każde z wykonanych zdjęć płytek jako dane wejściowe, generując jeden plik wynikowy Colonyzer (FT5 Tabela 1) dla każdego wykonanego zdjęcia. Pliki wynikowe Colonyzer określają szacunki gęstości kultury, powierzchnię kultury, kształt i kolor dla każdego z 384 miejsc na obrazowanej płytce. Nazwa pliku wynikowego jest automatycznie kopiowana z nazwy pliku obrazu źródłowego (np. fotografia płytki K000011_030_010_2011-09-22_09-47-54.jpg (FT1, Tabela 1) odpowiada plikowi wynikowemu Colonyzer K000011_030_010_2011-09-22_09-47-54.dat (FT5, Tabela 1)).
- Za pomocą pakietu R QFA załaduj metadane eksperymentalne (FT2, Tabela 1) oraz pliki wynikowe Colonyzer (FT5 Tabela 1). Dane te są łączone w ramki danych R (które mogą zostać wyeksportowane jako tekstowe pliki surowych danych QFA (FT6, Tabela 1)) do dalszej analizy.
- Pakiet R QFA zawiera funkcje służące do gromadzenia przebiegów czasowych gęstości komórek dla każdej kultury, dopasowywania uogólnionych logistycznych modeli populacyjnych do obserwacji oraz ich wykreślania (przykłady znajdują się na Rysunkach 2 i 3). Dopasowane wartości parametrów są zapisywane w plikach parametrów logistycznych QFA (FT7 Tabela 1). Więcej szczegółów znajduje się w dokumentacji pakietu R QFA.
- Pakiet R QFA zawiera również funkcje kontroli jakości: kolonie brzegowe są odrzucane ze względu na większą dostępność składników odżywczych na brzegach płytki i trudności w analizie obrazu w pobliżu ścianek płytki; z analizy usuwane są kultury, które nie przeszły testu SGA, oraz genotypy wykazujące sprzężenie z markerowymi genami specyficznymi dla przesiewu.
- Na podstawie parametrów logistycznego modelu populacyjnego dla każdej kultury wyznacza się kilka ilościowych miar sprawności (fitness). Obejmują one maksymalne tempo podwojenia populacji (MDR) oraz liczbę podziałów od inokulacji do nasycenia (MDP). Dla każdego zestawu kultur powtórzonych (na przykład dla unikalnego genotypu) i dla każdej definicji sprawności oblicza się kilka statystyk podsumowujących szacunki sprawności: średnią i medianę sprawności, odchylenie standardowe sprawności oraz liczbę zaobserwowanych powtórzeń. Statystyki podsumowujące są zapisywane w plikach podsumowania sprawności QFA (FT8, Tabela 1) do dalszej analizy, np. obliczenia wyników interakcji genetycznych (FT9 Tabela 1).
7. Reprezentatywne wyniki
Rysunek 2 przedstawia jedną typową krzywą wzrostu spośród 384 kultur QFA rosnących w temperaturze 20 °C na stałym agarze, zmierzoną za pomocą Colonyzer, w której wzrost gęstości komórek staje się wykrywalny po około 2 dniach. Opóźnienie to jest prawdopodobnie łącznym efektem fazy adaptacyjnej kultury po inokulacji na podłoże stałe oraz dolnej granicy wykrywalności gęstości komórek. Rysunek 3 pokazuje 308 podobnych krzywych wzrostu zarejestrowanych z jednej płytki. Rysunek 4 demonstruje porównanie dwóch genomowych przesiewów QFA w celu określenia siły interakcji genetycznych (adaptowano z 1). Pakiet R dla QFA 7 zawiera również funkcje służące do generowania Rysunku 3 oraz do wyjściowego zestawienia rankingowych list badanych genotypów i szacunków siły interakcji genetycznej (GIS), wraz z wartością q (wartość p skorygowana o stopień fałszywych odkryć (FDR)) dla istotności zaobserwowanego GIS.

Rysunek 1. Schemat robotycznej inokulacji szczepów drożdży w formacie 384 punktów. Procedura ta rozpoczyna się od 1536 niezależnych kultur na płytce (po lewej). W tym typowym przykładzie kolonie w pozycjach 1,1; 1,2; 2,1 i 2,2 (zaznaczone na czerwono) stanowią cztery powtórzenia tego samego genotypu. Kultury his3::KANMX zaznaczone na żółto, rosnące na krawędzi płytki, mają przewagę wzrostu ze względu na brak konkurencji i dlatego nie są badane metodą QFA. Jedno z tych powtórzeń (np. 1,1) jest inokulowane do płynnych mediów wzrostowych w płytkach 96-dołkowych przy użyciu narzędzia 96-igłowca, który za każdym razem inokuluje 96 z 1536 kolonii. Aby zaszczepić jedno powtórzenie dla każdego z 384 delecji genów, cztery różne „kwadranty” (oznaczone jako czerwony, niebieski, zielony i fioletowy) są inokulowane do czterech różnych płytek 96-dołkowych zawierających media wzrostowe. Po wzroście do nasycenia (np. 3 dni w 20 °C), kultury są rozcieńczane w wodzie, a następnie cztery kwadranty z jednego powtórzenia są naniesione w formacie 384 na stałą płytkę agarową (po prawej) w tym samym układzie co na oryginalnej płytce SGA (zgodnie z oznaczeniem kolorystycznym). Proces ten można powtórzyć w celu przetestowania pozostałych powtórzeń: 1,2; 2,1 i 2,2. Przykładowe obrazy time-lapse po prawej stronie zostały wykonane 0,5, 2 i 3,5 dni po inokulacji. Adaptacja z Rysunku uzupełniającego 1 2. Kliknij tutaj, aby wyświetlić większy rysunek.

Rycina 2. Pojedyncza krzywa wzrostu QFA zarejestrowana przez robota. A) Krzywa wzrostu QFA szczepu drożdży rosnących na agarze z galaktozą w temperaturze 20 °C. Obrazy były rejestrowane robotycznie mniej więcej co 2 godziny. Faza wykładnicza w tej temperaturze była obserwowalna przez około 1,5 dnia, a wzrost gęstości hodowli stał się wykrywalny około 2 dni po inokulacji. Do zaobserwowanych danych dopasowano uogólniony model logistyczny (krzywa szara). Przedstawiono parametry modelu automatycznie dopasowane przez pakiet R dla QFA: K (pojemność środowiska (AU)), r (szybkość wzrostu (d-1)), g (gęstość inokulum (AU)), v (symetria wzrostu). B) Podobnie jak w panelu A, wykreślono z gęstością komórek w skali logarytmicznej. Kliknij tutaj, aby wyświetlić powiększoną rycinę.

Rysunek 3. Automatycznie generowane krzywe wzrostu zarejestrowane równolegle z jednej płytki w formacie 384. Każdy panel przedstawia szacunki gęstości komórek (czerwone krzyżyki) i dopasowanie modelu (czarne krzywe) dla niezależnych kultur oznaczonych nazwą genu lub numerem y-ORF hodowanych w procedurze QFA. Ta przykładowa płytka została zobrazowana przez robota i hodowana w automatycznym inkubatorze w temperaturze 20 °C. Metadane eksperymentalne są automatycznie zawarte w tytule rysunku: nazwa płytki, traktowanie płytki, podłoże agarowe i numer płytki bibliotecznej. Dopasowane wartości parametrów modelu są również automatycznie wydrukowane na każdym panelu (patrz Rysunek 2). Kultury krawędziowe zostały usunięte z analizy QFA, co pozostawiło 308 kultur na tym wykresie (patrz Protokół QFA 5.4). Ponieważ kultury krawędziowe nie są analizowane w QFA, miejsca te są zazwyczaj wypełniane neutralnymi szczepami kontrolnymi (w tym przypadku pGAL-HIS3). Oryginalna wersja tego rysunku, wygenerowana przez pakiet R QFA, znajduje się w formacie .pdf i jest zatem nieskończenie powiększalna oraz przeszukiwalna za pomocą zapytań tekstowych. Kliknij tutaj, aby wyświetlić powiększony rysunek.

Rycina 4. Porównanie dostosowania z dwóch przesiewów QFA w celu wywnioskowania interakcji genetycznej.Wykres ten przedstawia porównanie dostosowania (fitness) delecji drożdży pączkujących (yfgΔ) w połączeniu z neutralną mutacją ura3Δ lub mutacją wrażliwą na temperaturę cdc13-1, hodowanych w temperaturze półprzepuszczalnej dla cdc13-1 (27 °C). Dostosowanie obliczono jako iloczyn maksymalnej szybkości podwajania i maksymalnego potencjału podwajania 1. Większość szczepów z delecjami rośnie dobrze w połączeniu z neutralną mutacją ura3Δ, co było oczekiwane, i wykazuje wysokie dostosowanie, natomiast delecje w połączeniu z cdc13-1 wykazują szeroki zakres dostosowania. Jasnoniebieskie linie przecinają się w miejscu neutralnej mutacji his3Δ; linia ciągła stanowi model regresji liniowej niezależności genetycznej (oczekiwane dostosowanie mutantów cdc13-1 yfgΔ przy danym dostosowaniu mutantów ura3Δ yfgΔ), a linia przerywana jest linią równości dostosowania. Delecje zaznaczone na zielono znacząco wzmacniają defekt dostosowania szczepów cdc13-1. Delecje zaznaczone na czerwono znacząco supresują ten sam defekt. Odległość pionowa dowolnej mutacji cdc13-1 yfgΔ od linii ciągłej wskazuje siłę interakcji genetycznej między cdc13-1 a yfgΔ w warunkach hodowli na płytkach. Należy zauważyć, że niektórymi z najsilniejszych supresorów cdc13-1 są wcześniej zidentyfikowane geny punktu kontrolnego RAD17, DDC1 i RAD24 oraz nukleaza EXO1 w prawym górnym rogu. Geny, których delecja zmniejsza dostosowanie, obejmują YKU80 (kodujący białko naprawy DNA i cappingu telomerów Yku80) w prawym dolnym rogu. Warto również zauważyć, że niektóre grupy genów funkcjonujących wspólnie mają tendencję do współlokalizacji na tym wykresie. Trzy przykładowe klastry są zaznaczone na niebiesko: SPE1, SPE2 i SPE3: geny syntezy espermidyny; RAD17, DDC1 i RAD24: kodujące składniki przesuwnej obręczy punktu kontrolnego (checkpoint sliding clamp) i ładowarki obręczy punktu kontrolnego (checkpoint clamp loader); MRE11, RAD50 i XRS2: kompleks MRX zaangażowany w DDR oraz EST1 i EST3: geny zaangażowane w regulację długości telomerów. Adaptowano z Supplementary Figure 1B 1; surowe dane, na podstawie których wygenerowano ten wykres, można pobrać stąd: http://research.ncl.ac.uk/colonyzer/AddinallQFA/. Kliknij tutaj, aby wyświetlić powiększoną rycinę.
| Id pliku | Typ pliku | Utworzony ręcznie | Częstotliwość aktualizacji | Wymóg wstępny |
| FT1 | Fotografia płytki | Nie | Na każdą płytkę i punkt czasowy | Tak |
| FT2 | Opis eksperymentu | Tak | Na każdy eksperyment | Tak |
| FT3 | Opis biblioteki | Tak | Na każdą bibliotekę | Tak |
| FT4 | Standardowe nazwy genów | Tak | Na każdą bibliotekę | Nie |
| FT5 | Wynik Colonyzer | Nie | Na każdy punkt czasowy na płytkę | Tak |
| FT6 | Surowe dane QFA | Nie | Na każdy eksperyment | - |
| FT7 | Parametry logistyczne QFA | Nie | Na każdy eksperyment | - |
| FT8 | Podsumowanie sprawności QFA | Nie | Na każdy eksperyment | - |
| FT9 | Lista trafień interakcji genetycznych QFA | Nie | Na każde badanie GIS (2 eksperymenty) | - |
Tabela 1. Elektroniczne pliki QFA. Wszystkie wymienione pliki (z wyjątkiem FT1) są plikami tekstowymi rozdzielanymi tabulatorami i mogą być tworzone lub odczytywane przy użyciu dowolnego edytora tekstu lub programu do arkuszy kalkulacyjnych. Aby uzyskać więcej szczegółów na temat tworzenia plików z metadanymi QFA i dokładnej interpretacji wyników QFA, należy zapoznać się z dokumentacją pakietu R dla QFA 7 oraz dokumentacją oprogramowania Colonyzer 3.

Tabela 2. Metody oceny kondycji hodowli. Podsumowanie cech i wymagań szeregu metod oceny kondycji hodowli drobnoustrojów. Przegląd autorstwa Blomberga 8 zawiera bardziej szczegółowe porównanie niektórych z tych metod. Kliknij tutaj, aby wyświetlić większą grafikę.