Schemat blokowy na Rysunku 1 przedstawia ogólne etapy procedury, w tym punkty, w których mierzone są częstości i tempo mutacji oraz wiek komórek. Dokładne określenie częstości i tempa mutacji (lub innej formy niestabilności genomu) w młodych populacjach komórek jest ważnym pierwszym krokiem, ponieważ jest to niezbędne do wyboru odpowiedniej wielkości populacji do znakowania i sortowania magnetycznego. Wartości tempa mutacji można ustalić za pomocą testów fluktuacyjnych9. Przykładowe wyniki dla szczepów drożdży w tle genetycznym BY474120, wyselekcjonowanych pod kątem mutacji w genie CAN1 nadających oporność na kanawaninę lub mutacji w genie URA3 nadających oporność na 5-FOA, przedstawiono na Rysunku 2. Komórki hodowano w temperaturze 20 °C dla reprezentatywnych eksperymentów dotyczących tempa mutacji oraz w 20 °C i 30 °C dla reprezentatywnych eksperymentów dotyczących częstości mutacji. Temperatury te zastosowano, ponieważ początkowe populacje komórek hodowane w 30 °C były następnie hodowane w 20 °C w kolejnych reprezentatywnych eksperymentach sortowania. Wartości tempa mutacji z niezależnych prób były porównywalne dla początkowej populacji komórek i komórek po znakowaniu biotyną dla szczepu, który posiada gen CAN1 w normalnej lokalizacji na chromosomie V, około 32 kilobaz od telomeru lewego ramienia (www.yeastgenome.com) (Rysunek 2A, kolumny lawendowe). Tempo mutacji CAN1 było znacznie wyższe w populacjach początkowych drugiego szczepu drożdży, który posiada delecję genu CAN1 z jego normalnej lokalizacji na chromosomie V oraz insercję CAN1 na prawym ramieniu chromosomu VIII, około 25 kilobaz od telomeru21 (Rysunek 2B). Częstości mutacji genu CAN1 były również porównywalne dla początkowych populacji komórek i komórek znakowanych biotyną w obu temperaturach wzrostu (Rysunek 2A, kolumny niebieskie). Jednakże stwierdzono, że częstości mutacji uzyskane w 20 °C były wyższe niż częstości mutacji w 30 °C. Aby sprawdzić, czy ten efekt temperatury ograniczał się do genu CAN1, zmierzyliśmy mutacje w URA3, selekcjonując oporność na 5-FOA przy użyciu szczepu drożdży z insercją CAN1 na chromosomie VIII. Szczep ten posiada również delecję URA3 z jego normalnego miejsca na chromosomie V oraz insercję URA3 na prawym ramieniu chromosomu VIII, około 40 kilobaz od telomeru21. Częstości mutacji dla URA3 w tej lokalizacji chromosomowej były również wyższe w 20 °C (Rysunek 2C). Ogólnie dowodzi to, że temperatura wzrostu może wpływać na częstość mutacji, natomiast znakowanie biotyną nie wydaje się wpływać na częstość mutacji. W związku z tym tempo i częstość mutacji dla populacji początkowych muszą zostać określone w tej samej temperaturze wzrostu, która będzie stosowana podczas cykli wzrostu i sortowania. Zaleca się zweryfikować, czy znakowanie biotyną nie wpływa na niestabilność genomu przed zastosowaniem protokołu do badania innych rodzajów mutacji lub rearanżacji genomu.
Zgodnie z oczekiwaniami, wartości tempa mutacji przedstawione na Rysunku 2A są kilkukrotnie niższe niż odpowiadające im pomiary częstości mutacji. Różnica ta wynika z faktu, że tempo mutacji mierzy pojawianie się komórek z nowymi mutacjami w każdym cyklu podziału komórkowego, podczas gdy częstość mutacji mierzy skumulowaną liczbę komórek zmutowanych w populacji na koniec okresu wzrostu. Tempa oraz 95% przedziały ufności dla tych prób wskazują, że powtarzalne wyniki można uzyskać, stosując siedem kultur powtórzonych na próbę, co stanowi minimalną liczbę powtórzeń sugerowaną w niniejszym protokole.
Po ustaleniu wstępnych wartości częstotliwości i tempa należy dobrać wielkość populacji, która zapewni obecność odpowiedniej liczby komórek po wielu rundach sortowania, aby wiarygodnie wyznaczyć częstotliwości mutacji. Oznaczenie populacji 108 komórek dla każdego punktu czasowego analizowanego podczas starzenia jest odpowiednie, gdy częstotliwości i tempa są zbliżone do tych przedstawionych na Rysunku 2A. Decyzja ta zależy również od tego, jak wydajnie odzyskuje się oznakowane komórki macierzyste po każdej rundzie sortowania. Wydajność odzyskiwania oznakowanych komórek macierzystych można szybko i prosto ocenić, używając hemocytometru do wyznaczenia liczby komórek wymytych z kolumn po każdym sortowaniu oraz do zbadania morfologii komórek. Przykład danych dotyczących wydajności odzyskiwania (Rysunek 3A) ilustruje dwa główne punkty: odzyskana liczba komórek może wydawać się wyższa niż oczekiwano po pierwszym sortowaniu, a przy dalszym sortowaniu spodziewana jest niewielka, postępująca utrata komórek. Wyższa niż oczekiwana liczba komórek w niektórych eksperymentach po pierwszym sortowaniu może wynikać z oznakowania pączków w komórkach populacji początkowej. Komórka drożdży z pączkiem byłaby zazwyczaj liczona jako pojedyncza komórka podczas wyznaczania liczby komórek do oznakowania biotyną. Jednak oznakowanie pączków obecnych w populacji początkowej pozwoliłoby na zatrzymanie na kolumnach podczas sortowania również komórek, które rozwiną się z tych pączków. Wpływ tego zjawiska na wyznaczenie wydajności odzyskiwania zależy od frakcji komórek pączkujących w początkowej populacji komórek. Drugim potencjalnym powodem wyższej liczby komórek po pierwszym sortowaniu jest tendencja do występowania w tym punkcie czasowym większej liczby dużych komórek pączkujących, które mogą kończyć podział komórkowy podczas sortowania. W rezultacie duże pączki wciąż przyczepione do komórek macierzystych mogą zostać posortowane, a następnie pojawić się jako odrębne komórki podczas badania wymytych komórek. Chociaż przy każdej rundzie wzrostu i sortowania obserwuje się pewną utratę komórek, protokół może zapewnić wiarygodne zatrzymanie 80-90% komórek po każdej rundzie sortowania. Warto poćwiczyć tę procedurę, aby upewnić się, że izolowanych jest 80% lub więcej oznakowanych komórek, co pozwoli uniknąć konieczności oznakowania niepotrzebnie dużej początkowej populacji komórek. Poddanie komórek łagodnemu stresowi po pierwszym sortowaniu (1 mM nadtlenku wodoru w medium YPD przez 30 min) nie zapobiegło wydajnemu odzyskiwaniu komórek w kolejnych rundach wzrostu i sortowania, choć po stresie wystąpiła pewna zmienność w odzyskiwaniu podczas pierwszego sortowania (Rysunek 3A).
Ocena morfologii i żywotności komórek jest przydatna zarówno do potwierdzenia pomyślnej izolacji komórek macierzystych, jak i do dostosowania rozcieńczeń/objętości stosowanych podczas wysiewania komórek na podłoże nieselektywne w celu wyznaczenia gęstości jednostek tworzących kolonie. Komórki macierzyste stają się coraz większe i przyjmują bardziej nieregularne kształty wraz z wiekiem w porównaniu do komórek potomnych (Ryc. 3B). Próbki komórek macierzystych powinny zatem składać się przede wszystkim z dużych komórek o nieregularnym kształcie w miarę postępu eksperymentu przez kolejne rundy sortowania. Obecność wielu małych, owalnych komórek o regularnym kształcie może wskazywać, że komórki macierzyste nie są odpowiednio oddzielone od komórek potomnych. Żywotność komórek macierzystych powinna również stopniowo spadać z każdą rundą hodowli i sortowania, choć w ciągu pierwszych pięciu do dziesięciu generacji komórkowych zmiany mogą być niewielkie. Liczba komórek zdolnych do tworzenia kolonii może spadać znacznie gwałtowniej niż liczba komórek uznanych za żywe na podstawie bezpośredniego barwienia integralności komórkowej, co wynika z powstawania komórek senescentnych. Gdy żywotność mierzona metodą bezpośredniego barwienia zaczyna po raz pierwszy spadać (do około 80% lub mniej), konieczne może być dostosowanie rozcieńczeń i objętości używanych do pomiaru gęstości jednostek tworzących kolonie, przewidując gwałtowniejszy spadek zdolności żywych komórek do tworzenia kolonii (spadek dwukrotny lub wielokrotny).
Przed analizą danych dotyczących częstości mutacji w komórkach macierzystych w celu wykazania zmian w tempie akumulacji mutacji zależnych od wieku, należy wyznaczyć wiek replikacyjny posortowanych populacji oraz populacji kontrolnych. Można to osiągnąć poprzez manualne liczenie blizn po pączkowaniu w komórkach macierzystych (Rysunek 4) lub za pomocą cytometrii przepływowej w celu ilościowego oznaczenia sygnału odczynnika do wykrywania blizn po pączkowaniu (Rysunek 5). Blizny po pączkowaniu można znakować przy stosunkowo niskim sygnale tła, używając konjugatów WGA z fluorochromem (Rysunek 4C). Rysunek 4A pokazuje, że większość komórek z próbek przepływowych procedury sortowania posiada zero lub jedną bliznę po pączkowaniu. Komórki eluowane z kolumn to głównie starsze komórki macierzyste (Rysunki 4B i 5). Zazwyczaj >90% eluowanych komórek to komórki macierzyste, co można łatwiej zaobserwować na wyniku cytometrii przepływowej na Rysunku 5. Komórki z relatywnie niewielką liczbą blizn po pączkowaniu (<6) uzyskane po więcej niż dwóch rundach sortowania mogą reprezentować komórki zanieczyszczające, które nie zostały znakowane biotyną i przeszły kilka rund podziałów komórkowych podczas odpowiedniej fazy wzrostu, w przeciwieństwie do komórek macierzystych dzielących się bardzo powoli. Zmienność wieku komórek może wzrosnąć w kolejnych rundach sortowania, jeśli nie wszystkie komórki w populacji rosną w sposób jednolity.
Większą populację komórek można wykorzystać do szybszej oceny wieku komórek, jeśli zostanie ustanowiona liniowa zależność między intensywnością sygnału fluorescencji WGA a liczbą blizn pączkowania na komórkę. W tej analizie normalizację wszystkich intensywności sygnału fluorescencji WGA przeprowadza się poprzez podzielenie sygnału z barwionych komórek przez sygnał otrzymany dla odpowiedniej niebarwionej populacji komórek (córki lub matki), aby uwzględnić zwiększoną fluorescencję tła w starszych komórkach. Ryciny 4 i 5 przedstawiają analizę tych samych reprezentatywnych populacji komórek. Ręczne liczenie pozwoliło ustalić średni wiek replikacyjny wynoszący odpowiednio 0,95 i 11,4 dla populacji komórek córki (Rycina 4A) oraz komórek matki (Rycina 4B). Znormalizowane sygnały WGA dla tych dwóch populacji oraz trzech dodatkowych populacji (średni wiek 3,0, 6,9 i 14,4) naniesiono na wykres w funkcji średniej liczby blizn pączkowania dla każdej populacji (Rycina 5B). Tę liniową zależność można następnie wykorzystać z tym samym szczepem i odczynnikami w przyszłych eksperymentach w celu określenia średniego wieku komórek na podstawie znormalizowanego sygnału WGA uzyskanego za pomocą cytometrii przepływowej, co pozwala na szybkie i dokładne określenie wieku komórek. Należy nadal uwzględniać populacje kontrolne, aby zweryfikować podobną wydajność barwienia pomiędzy seriami doświadczalnymi.
Wpływ wieku na akumulację mutacji można przeanalizować po wykonaniu wcześniejszych kroków polegających na wyznaczeniu wartości tempa mutacji, efektywnym sortowaniu komórek oraz określeniu wieku komórek. Liczba punktów czasowych, w których można wyznaczyć częstotliwość, zależy od wielkości populacji komórek znakowanych biotyną oraz liczby komórek, które należy wysiać na podłożu selektywnym w celu uzyskania wiarygodnych wyników. Jeśli tempo mutacji zmienia się wraz z wiekiem, to obserwowane częstotliwości mutacji dla starzejących się komórek macierzystych powinny różnić się od tych oczekiwanych wyłącznie na podstawie dodatkowych rund podziałów komórkowych. Porównanie obserwowanej częstotliwości mutacji z przewidywaną częstotliwością mutacji pozwala zidentyfikować różnice w tempie akumulacji mutacji związane z wiekiem. Jak opisano w sekcji 7 protokołu, przewidywaną częstotliwość można otrzymać z iloczynu bazowego tempa mutacji dla początkowej populacji komórek i wzrostu wieku replikacyjnego komórek macierzystych, dodanego do częstotliwości mutacji dla populacji początkowej. Wzrosty częstotliwości mutacji CAN1 zaobserwowano wraz ze wzrostem średniego wieku komórek w szczepie z CAN1 w normalnej lokalizacji na chromosomie V (Rycina 6A) oraz w szczepie z CAN1 na prawym ramieniu chromosomu VIII (Rycina 6B). Ponieważ obserwowane częstotliwości mutacji dla komórek macierzystych na Rycini 6A są podobne do przewidywanych częstotliwości lub znajdują się tuż poniżej nich, dane te nie dostarczają dowodów na specyficzną dla wieku zmianę tempa mutacji. W przeciwieństwie do tego, częstotliwości mutacji dla komórek macierzystych szczepu z CAN1 na chromosomie VIII były wyższe niż przewidywane częstotliwości (Rycina 6B). Należy zauważyć, że przewidywane częstotliwości na wykresie nie wydają się znacząco rosnąć, ponieważ ze względu na duży wzrost obserwowanej częstotliwości mutacji dla osi y należało zastosować skalę logarytmiczną. Zatem wyniki na Rycini 6B dostarczają dowodów wspierających specyficzny dla wieku wzrost tempa akumulacji mutacji. Różnica w wynikach dla zestawów danych na Rycini 6A i 6B wynika prawdopodobnie z różnych lokalizacji genomicznych CAN1. Tego typu obserwacje stanowią punkt wyjścia do formułowania hipotez w celu zbadania mechanizmów wpływających na tempo mutacji w miarę starzenia się komórek oraz do opracowania modeli wyjaśniających akumulację mutacji w zależności od wieku replikacyjnego.

Rysunek 1. Schemat blokowy ogólnej procedury badania akumulacji mutacji podczas replikacyjnego starzenia się drożdży. Czarny tekst i strzałki wskazują główne etapy procedury. Zakrzywiona strzałka odzwierciedla fakt, że dodatkowe cykle ponownego wzrostu i sortowania tych samych komórek są wykorzystywane w celu uzyskania komórek o coraz większym wieku. Fioletowe strzałki i tekst wskazują etapy, na których wykonywane są określone pomiary. Freq – częstotliwość.

Rycina 2. Wyznaczenie częstotliwości i tempa mutacji w młodych populacjach komórek. (A) Kolonie mutantów wyselekcjonowano poprzez wysiew komórek na podłoże SC-arg + kanawanina w celu wyselekcjonowania mutacji utraty funkcji w genie CAN1, a następnie porównano je z całkowitą liczbą żywych komórek wysianych na podłoże selekcyjne, aby obliczyć częstotliwości/tempa. Komórki z populacji początkowych przed znakowaniem biotyną (IP) lub po znakowaniu biotyną (PB) hodowano na podłożu YPD, zwiększając gęstość z początkowej wartości 5,000 cells/ml do poziomu bliskiego nasyceniu. Kolumny lawendowe oznaczają pomiary tempa, a kolumny niebieskie reprezentują pomiary częstotliwości wykonane w podanych temperaturach (20 °C i 30 °C). Dla każdej próby hodowano zestawy siedmiu kultur powtórzonych, a wykonano od dwóch do pięciu niezależnych prób. Przedstawiono indywidualne wartości tempa obliczone dla niezależnych prób, a słupki błędów dla tych prób reprezentują 95% przedziały ufności wyznaczone za pomocą kalkulatora online22. Częstotliwości reprezentują średnie i odchylenia standardowe. (B) Tempa mutacji CAN1 uzyskane i przedstawione zgodnie z opisem w części A, z wykorzystaniem szczepu z genem CAN1 na prawym ramieniu chromosomu VIII. (C) Częstotliwości mutacji uzyskane po selekcji kolonii mutantów na 5-FOA z wykorzystaniem szczepu z genem URA3 zlokalizowanym na prawym ramieniu chromosomu VIII. Pozostałe metody były analogiczne do części A; przedstawiono średnie i odchylenia standardowe dla trzech lub czterech niezależnych prób.

Rycina 3. Przykład wydajności sortowania określonej poprzez obliczenie liczby eluowanych komórek po każdej rundzie sortowania. (A) Całkowitą liczbę komórek w każdej populacji wyjściowej po znakowaniu biotyną przyjęto za wartość jeden, a całkowitą liczbę komórek obecnych w próbkach eluowanych z każdej rundy magnetycznego sortowania komórek podzielono przez te wartości początkowe, aby uzyskać frakcję komórek znakowanych. Całkowita liczba komórek została określona na podstawie zliczeń wykonanych za pomocą hemocytometru. Pomarańczowe słupki reprezentują średnią i odchylenie standardowe dla trzech niezależnych prób przeprowadzonych zgodnie ze standardowym protokołem. Niebieskie słupki reprezentują średnią i odchylenie standardowe dla dwóch niezależnych prób, w których komórki traktowano 1 mM nadtlenkiem wodoru przez 30 min bezpośrednio po pierwszym sortowaniu. (B) Rozmiar i morfologia komórek po znakowaniu biotyną (post-biotin) oraz komórki macierzyste odzyskane po trzeciej rundzie sortowania magnetycznego (sort 3 mother cells) przedstawione przy użyciu standardowej mikroskopii w polu jasnym i obiektywu 20X. Białe linie w tle to linie ograniczające najmniejsze kwadraty widoczne w standardowym hemocytometrze. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 4. Określanie wieku replikacyjnego za pomocą ręcznego liczenia blizn po pączkowaniu. Mikroskopia konfokalna została wykorzystana do policzenia blizn po pączkowaniu w poszczególnych komórkach znakowanych koniugatem fluorescencyjnym WGA (pobudzenie przy 488 nm i detekcja z kombinacją filtra pasmowo-przepustowego 458/543 nm oraz filtra długoprzepustowego 505 nm). (A) Ręczne liczenie blizn po pączkowaniu w komórkach z frakcji przepływowej (komórki potomne) po trzeciej rundzie sortowania magnetycznego (n = 62). (B) Ręczne liczenie blizn po pączkowaniu w komórkach elowanych (komórki macierzyste) po trzeciej rundzie sortowania magnetycznego (n = 54). (C) Komórki macierzyste zatrzymane po trzeciej rundzie sortowania magnetycznego, sfotografowane przy użyciu obiektywu imersyjnego 63X z 3-krotnym powiększeniem po barwieniu koniugatem fluorescencyjnym WGA. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 5. Określanie wieku replikacyjnego za pomocą cytometrii przepływowej. Próbki po 10 000 komórek zabarwionych koniugatem WGA z fluoresceiną analizowano metodą cytometrii przepływowej przy wzbudzeniu przy 488 nm i detekcji z użyciem zestawu filtrów pasmowoprzepuszczalnego 530/30 oraz długoprzepuszczalnego 505. (A) Histogramy przedstawiające liczbę komórek o określonej intensywności fluorescencji WGA dla populacji komórek potomnych lub eluatu (komórek macierzystych) po trzech rundach sortowania magnetycznego. Populacje odpowiadają tym analizowanym na Rycini 4. Niebieskie linie pionowe wskazują odpowiadającą pozycję dla większości komórek potomnych na histogramie komórek macierzystych. (B) Znormalizowaną średnią geometryczną sygnału fluorescencji dla każdej populacji przedstawionej w A oraz trzech dodatkowych populacji komórek (średni wiek 3,0, 6,9 i 14,4) obliczono jako stosunek średniej geometrycznej komórek zabarwionych do średniej geometrycznej odpowiedniej populacji komórek niezabarwionych. Wartości te naniesiono w porównaniu do średniej liczby blizn pączkowania dla każdej całkowitej populacji (Rycina 4 oraz dane niepokazane). Wartość R2 dla linii trendu tego porównania jest podana na wykresie.

Rycina 6. Porównanie przewidywanej i obserwowanej częstości mutacji podczas starzenia replikacyjnego. Komórki z mutacjami w genie CAN1 wyselekcjonowano zgodnie z opisem dla Ryciny 2, wykorzystując szczepy z genem CAN1 w normalnej lokalizacji na chromosomie V (A) lub na prawym ramieniu chromosomu VIII (B). Komórki hodowano na podłożu stałym w temperaturze 30 °C przed znakowaniem biotyną, a następnie w 20 °C. Obserwowane wartości częstości przedstawiono za pomocą pomarańczowych słupków (Obs), a przewidywane częstości dla komórek starszych za pomocą brązowych słupków (Pre). Pierwszy słupek na każdym wykresie przedstawia średnią i odchylenie standardowe początkowej częstości mutacji po znakowaniu biotyną dla czterech niezależnych prób (PB). Liczby pod słupkami wskazują średnią liczbę blizn po pączkowaniu dla każdej populacji (Wiek komórki). Niezależne wartości przewidywane określono zgodnie z opisem w sekcji 7 protokołu, wykorzystując każdą z początkowych wartości szybkości przedstawionych na Rycini 2A (słupki IP) oraz 2B. Następnie wyciągnięto średnią z tych niezależnych wartości. Wartości obserwowane są średnimi z trzech prób. Paski błędów oznaczają odchylenie standardowe.