Sortowanie FACS
Strategie bramkowania (gating) przedstawione w niniejszym protokole opierają się na powszechnie akceptowanych metodach12,30,31. W przypadku strategii bramkowania przedstawionej na Rysunku 1, materiałem wyjściowym są progenitorzy z krwi pępowinowej uprzednio oczyszczeni via magnetyczne wzbogacenie CD34+, co wyjaśnia znikomy odsetek komórek Lineage-dodatnich. Niezbędne jest zastosowanie wąskich bram dla czterech kombinacji barwników wewnątrzkomórkowych (np. CTV na rysunku), aby poprawić rozdzielczość pików podczas następującej analizy oraz aby wyznaczyć właściwą populację komórek (Rysunek 1D). W przypadku pokazanym na rysunku, bramki selekcjonują największą i najlepiej zdefiniowaną populację. Obecność wielu, blisko położonych populacji dla każdej kombinacji barwników podziału komórkowego nie jest w naszym doświadczeniu reprezentatywna dla różnic biologicznych. Może ona zamiast tego wskazywać na a) nieoptymalną procedurę barwienia lub b) dużą heterogenność (zwłaszcza w rozmiarze) w wyjściowej puli komórek. Nie jest to zaskakujące w przypadku materiału z krwi pępowinowej lub innych złożonych źródeł biologicznych (np. aspiratów szpiku kostnego, krwi obwodowej). Jeśli brama nie jest ściśle zdefiniowana, postępujące rozcieńczanie różnych kombinacji barwników może prowadzić do zlewania się późniejszych pików, szczególnie w warunkach CV i VC (Rysunek 2D). Inną negatywną konsekwencją nieoptymalnego bramkowania jest brak możliwości efektywnego rozróżnienia poszczególnych pików po hodowli komórkowej, ponieważ heterogenna populacja wyjściowa może prowadzić do powstania niskich i rozmytych pików.

Rysunek 1: Strategia bramkowania dla sortowania komórek. (A) FSC-A względem SSC-A, w celu wykluczenia debris i komórek zanieczyszczających. (B) FSC-A względem FSC-H, w celu wykluczenia dubletów i skupisk komórek. (C) Lin względem FSC-H, w celu wykluczenia komórek Lin+. (D) CTV względem CFSE, aby jednoznacznie zidentyfikować komórki zabarwione kombinacjami barwników CF, CV, VC i VI. Bramki powinny być na tyle rygorystyczne, aby objąć jednorodną populację. (E) CD34 względem CD38, w celu oddzielenia ograniczonych progenitorów CD34+CD38+ (zwanych również HPCs) od komponentu multipotentnego CD34+CD38-. (F) CD45RA względem CD90, z populacji CD34+CD38-, w celu rozróżnienia między najmniej dojrzałymi progenitorami wzbogaconymi w HSC (CD90+CD45RA-), LMPP (CD90midCD45RA+) a bardziej zdeterminowanymi MPP (CD90-CD45RA-). (G) Indeks posortowanych zdarzeń, przedstawiony tutaj pod kątem kombinacji barwników komórkowych oraz (H) ekspresji markerów powierzchniowych CD90 i CD45RA. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Analiza cytometrii przepływowej po hodowli komórkowej
Dane przedstawione na Rycini 2 są reprezentatywne dla ludzkich HSCs z krwi pępowinowej, utrzymywanych w hodowli przez 72 h w obecności wielu cytokin zdolnych do wspierania szerokiej gamy progenitorów i prekursorów mieloidalnych. Panele od 2A do 2D przedstawiają bramkowanie niezbędne do ustalenia pokrewieństwa każdej z poszczególnych komórek, natomiast panele od 2E do 2G umożliwiają fenotypowanie komórkowe. Zmniejszona obecność MEPs na rycinie jest prawdopodobnie wynikiem warunków hodowli zastosowanych w tym reprezentatywnym eksperymencie (Rycina 2F). Zastosowanie innych cytokin i warunków hodowli zmienia względny procent każdego podzbioru, podobnie jak wybór różnych komórek wyjściowych do eksperymentu.

Rycina 2: Strategia bramkowania dla analizy cytometrycznej. (A) FSC-A względem SSC-A w celu wykluczenia debris i komórek zanieczyszczających. (B) FSC-A względem FSC-H w celu wykluczenia dubletów i agregatów komórkowych. (C) CTV względem CFSE; bramka Labeled pozwala na wykluczenie zdarzeń autofluorescencyjnych, które mogłyby wpłynąć na rozdzielczość danych. (D) CTV względem CFSE. Niezwykle ważne jest rygorystyczne bramkowanie czterech populacji na podstawie rozcieńczeń barwnika podziałów komórkowych. (E) CD34 względem CD38 w celu rozróżnienia prekursorów zdeterminowanych (CD34-), ograniczonych progenitorów (HPC) (CD34+CD38+) oraz niedojrzałych progenitorów (CD34+CD38-). (F) CD45RA względem CD123 w celu rozróżnienia trzech typów ograniczonych progenitorów: CMP (CD123+CD45RA-), MEP (CD123-CD45RA-) oraz GMP (CD123+CD45RA+). (G) CD45RA względem CD90 z populacji CD34+CD38- w celu identyfikacji HSCs, LMPPs oraz MPPs. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Etapy definiowania i przypisywania pików (Rycina 3 i Rycina 4) są kluczowymi aspektami protokołu i wymagają zdefiniowania rygorystycznych bramek (gates). Dla definicji piku (Rycina 3) niezbędne jest co najmniej 1 000 zdarzeń, aby identyfikacja była wiarygodna. W związku z tym korzystne może być wyizolowanie większej liczby komórek podczas etapu sortowania komórek do studzienek typu „Bulk”. Rycina 4 przedstawia cztery przykłady pojedynczych studzienek zawierających wiele rodzin. Rycina ta wyjaśnia znaczenie bramkowania z Ryciny 2D i Ryciny 3, szczególnie w odniesieniu do identyfikacji każdej rodziny i każdego piku. Rycina 4A ilustruje prosty przykład, w którym wszystkie komórki w bramce CF znajdują się bardzo blisko siebie i mogą zostać łatwo przypisane do jednego piku. Rycina 4C pokazuje inny przykład rodziny rozłożonej jednoznacznie na dwa wyraźnie oddzielone piki, co jest wyraźnie widoczne na histogramie na Rycynie 4D. Rycina 4E,G ujawniają znaczenie rygorystycznego bramkowania opartego na dużej liczbie zdarzeń; obie pokazują kilka zdarzeń, które są bliskie, ale znajdują się poza bramkami dla kombinacji barwników. Zdarzenia te mogłyby zostać błędnie włączone do bramek VI i CF, gdyby opierano się wyłącznie na analizie pojedynczej studzienki. Na koniec Rycina 4F,H przedstawiają dwa różne przykłady rodzin rozproszonych na wielu pikach: jeden przykład z dwoma pikami o podobnej intensywności (Rycina 4F) oraz jeden z dwoma pikami o różnej intensywności (Rycina 4H).

Rysunek 3: Definicja pików dla analizy cytometrii przepływowej. (A-D) Piki powinny być definiowane przy rejestracji co najmniej 500 zdarzeń, aby zapewnić dobrą reprezentację każdego pojedynczego piku. (A) Histogram dla intensywności CFSE-A. Można zidentyfikować kilka pików, z których każdy odpowiada innej populacji dzielących się komórek. (B,C) Histogramy dla intensywności parametru pochodnego, reprezentującego mieszaninę CFSE-CTV, odpowiednio CV (B) i VC (C). (D) Histogram dla intensywności CTV-A. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 4: Przypisanie piku. (A,B) Dla tej studzienki można wykryć tylko jeden pik, w bramce CF. (C,D) W tej studni można wykryć dwa piki o niemal równej intensywności, w bramce VI. Piki są dobrze rozdzielone. (E,F) W tej studni można wykryć dwa piki o zbliżonej intensywności, w bramce VI. Uwzględniono jedynie zdarzenia w bramce, zgodnie ze strategią opracowaną dla studni zbiorczych. (G-H) W tej studni można wykryć dwa piki o różnej intensywności, w bramce CF. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Reprezentacja danych i testy statystyczne
Rysunek 5 przedstawia różne typy reprezentacji danych z dwóch oddzielnych eksperymentów, oba przeprowadzonych po 72 h hodowli komórkowej. Komórki HSC i MPP były hodowane w dwóch różnych mediach hodowlanych, które miały na celu zmianę właściwości podziału i różnicowania komórek. Media te nazwano „Diff” (Differentiation – różnicowanie)32 oraz „GT”33; pierwsze promuje różnicowanie mieloidalne i erytroidalne, ponieważ zawiera erytropoetynę (EPO) oraz czynnik stymulujący kolonie granulocytowo-monocytarne (GM-CSF), natomiast drugie zostało opracowane w kontekście klinicznych badań nad terapią genową, z celem utrzymania i namnożenia wysokiego odsetka HSPC. Rysunek 5A to reprezentatywna mapa ciepła dla warunku „Diff”, przedstawiająca różnorodne rodziny komórek, zarówno pod kątem ich losów, jak i podziałów. Na tej mapie ciepła każdy wiersz reprezentuje pojedynczą rodzinę, każdy kwadrat pojedynczą komórkę, a kolumny grupują wszystkie komórki znajdujące się w tym samym pokoleniu (np. komórki w 2. pokoleniu podzieliły się co najmniej dwukrotnie). Można wyróżnić rodziny wysoce homogeniczne, składające się z jednego typu komórek i wykazujące tę samą liczbę podziałów (np. rodzina nr 63) oraz rodziny heterogeniczne, obejmujące trzy typy komórek w obrębie dwóch pokoleń (np. rodzina nr 84). Ponieważ wskaźnik odzyskiwania komórek dla tej analizy wynosi około 70%, rzadko obserwuje się rodziny kompletne, definiowane jako takie, w których odzyskano wszystkie komórki, potencjalnie w różnych pokoleniach (np. rodzina składająca się z jednej komórki w 1. pokoleniu i dwóch komórek w 2. pokoleniu) (oznaczone znakiem kratki obok numeru ID na Rysunku 5A). Istnieje wiele wyjaśnień niepełnej detekcji, które mogą być techniczne (problem z barwieniem, utrata komórek wynikająca z protokołu) lub biologiczne (śmierć komórek i/lub apoptoza). Ograniczenia techniczne można pokonać, stosując analizator zaprojektowany w celu zmniejszenia objętości martwej związanej z pojedynczą próbką oraz wykonując barwienie komórek bezpośrednio w płytce hodowlanej, aby zminimalizować pipetowanie objętości. Z kolei ortogonalne metody określania poziomu śmierci komórek (np. poprzez eksperymenty z obrazowaniem żywych komórek) mogą pomóc w odróżnieniu czynników technicznych i biologicznych prowadzących do niepełnej detekcji.
Rysunek 5Bi pokazuje, jak wizualizować wpływ warunków hodowli na skład typów komórek, analogicznie do wykonania analizy bulk. W tym przypadku warunek Diff sprzyja powstaniu większej liczby losów komórkowych oraz wyższemu odsetkowi komórek CD34+ (zdefiniowanych jako wszystkie typy komórek z wyjątkiem CD34-). Przedziały ufności zostały obliczone w skrypcie via podstawowy bootstrap, przy użyciu 250 000 zestawów danych z bootstrapu34. Warto zauważyć, że wszystkie pozostałe histogramy na Rysunku 5 prezentują przedziały ufności obliczone w ten sam sposób. Tabela 5 podsumowuje wszystkie informacje dotyczące liczby rodzin oraz liczby komórek w każdym pokoleniu.
Rysunek 5Bii graficznie przedstawia wyniki testów statystycznych przeprowadzonych w skrypcie „2_bar_plot”. Formalny opis ram statystycznych jest dostępny26. W skrócie, ramy te umożliwiają testowanie hipotez statystycznych przy założeniu, że komórki z tej samej rodziny są od siebie zależne (założenie to jest samo w sobie testowalne), w przeciwieństwie do statystyki klasycznej, która wymagałaby niezależności między wszystkimi obserwowanymi komórkami. W konkretnym przypadku przedstawionym na rysunku, test statystyczny weryfikuje hipotezę, że wybory losu komórkowego MPPs, mierzone jako częstości różnych typów komórek obecnych w hodowli, są niezależne od zastosowanych warunków hodowli komórkowej. Najpierw statystyka testu G jest wykorzystywana do oceny rozbieżności między częstościami typów komórek z różnych mediów hodowlanych (dla przykładu w Bii statystyka ta jest reprezentowana przez czerwony słupek). Następnie przeprowadzana jest randomizacja danych poprzez permutację, polegająca na zamianie całych rodzin komórek między dwoma warunkami hodowli. Ma to na celu zachowanie zależności między komórkami spokrewnionymi rodzinnie, przy jednoczesnym utrzymaniu liczby rodzin w każdym zbiorze zgodnie z danymi oryginalnymi. Statystyka testu G jest obliczana z zrandomizowanego zbioru danych. Wartości niebieskie przedstawione w 5Bii to statystyka testu G dla 250 000 permutacji. Na koniec obliczana jest wartość p, aby ocenić stopień, w jakim oryginalny zbiór danych odbiega od rozkładu zbiorów permutowanych. W tym przykładzie oryginalna statystyka znacząco odbiega od rozkładu, co skutkuje małą wartością p i tym samym odrzuceniem hipotezy, że los komórkowy MPPs jest niezależny od warunków hodowli.
Rycina 5C przedstawia procent rodzin komórek w zależności od maksymalnej generacji, aby zbadać, jak różne warunki wpływają na podziały komórkowe w obrębie pojedynczej rodziny. Wykres ten pokazuje, że po 72 h komórki hodowane w warunkach Diff wykonują większą liczbę podziałów niż komórki w warunkach GT. Przedstawiono liczbę maksymalnych generacji dla każdej rodziny, zatem rodzina, w której występują komórki z generacji 1 i 2, jest uznawana za generację 2. Do statystycznego sprawdzenia niezależności między podziałami komórkowymi a warunkami hodowli można zastosować ten sam model statystyczny, który wykorzystano w Rycini 5B.
Rysunek 5D analizuje rodzaj symetrii/asymetrii pierwszego podziału dla różnych typów komórek przodków (HSCs lub MPPs). Dla pełnych rodzin komórek w pokoleniu 1 – jedynym pokoleniu, w którym można jednoznacznie ustalić, że dwie komórki potomne są komórkami siostrzanymi – można zdefiniować cztery różne typy symetrii/asymetrii: etykieta „Sym Undiff” opisuje rodziny, w których obie komórki potomne zachowują fenotyp komórki macierzystej. „Sym Diff” oznacza, że obie komórki potomne mają ten sam fenotyp, który różni się od fenotypu komórki macierzystej. „Asym Undiff” oznacza, że tylko jedna komórka potomna zachowuje fenotyp komórki macierzystej. Wreszcie, „Asym Diff” opisuje rodziny, w których obie komórki potomne mają różne fenotypy i żadna z nich nie jest taka sama jak komórka macierzysta. Aby uzyskać odpowiednią moc statystyczną w ocenie tych symetrycznych/asymetrycznych losów, wskazane jest przeprowadzenie analizy MultiGen we wczesnych punktach czasowych, aby zaobserwować więcej rodzin, których potomstwo znajduje się w pokoleniu 1.
Na koniec, Rycina 5E przedstawia procentowy udział typów komórek w funkcji liczby podziałów, aby uzyskać wgląd w progresję wzorca różnicowania w kolejnych podziałach. Na przykład dane przedstawione na rycinie sugerują, że komórki przechodzą w stan CD34-, przy czym ponad 50% wykrytych komórek w tej klasie występuje już po trzech podziałach. Ponadto można wywnioskować, że MPPs nie preferują podziałów prowadzących do samoodnowienia, ponieważ tylko niewielki procent komórek zachowuje pierwotny fenotyp. Niektóre z tych wniosków można następnie zweryfikować, korzystając z ram statystycznych przedstawionych na poprzednich rycinach.

Rysunek 5: Przykład reprezentacji danych dla jednego eksperymentu trwającego 72 h z użyciem HSPC z krwi pępowinowej. (A) Mapy ciepła dla wybranego zestawu danych (HSC, w medium „Diff”, po 72 h hodowli). Wykresy przedstawiają wszystkie pojedyncze komórki (kwadraty) zgodnie z ich pokrewieństwem (wiersze), liczbą wykonanych podziałów (kolumny, nazwane generacjami) oraz fenotypem (kolory). (Bi) Histogram porównujący proporcje typów komórek w potomstwie komórek HSC i MPP pomiędzy warunkiem GT a warunkiem Diff. (Bii) Wykres przedstawia testy statystyczne wykonane w skrypcie „2_bar_plot” dla MPP po 72 h hodowli, porównujące koktajle cytokinowe „Diff” i „GT”. Wartość eksperymentalna jest wyświetlana na czerwono, a wartości wygenerowane poprzez 250 000 permutacji na niebiesko. W prawym górnym rogu wskazano wartość p. testu G wraz z liczbą rodzin użytych do testu. (C) Histogram porównujący procent rodzin (łącznie 314 rodzin) w każdej generacji (oznaczonej kolorami) dla HSC i MPP w zależności od warunków hodowli. Przedziały ufności obliczono na podstawie 250 000 zestawów danych metodą bootstrap. (D) Histogram przedstawiający typ symetrii/asymetrii pomiędzy losem komórek potomnych dla rodzin z dwiema komórkami w generacji 1: Sym Undiff (obie komórki potomne zachowują fenotyp komórki macierzystej), Sym Diff (obie komórki potomne mają ten sam fenotyp, który różni się od fenotypu komórki macierzystej), Asym Undiff (tylko jedna komórka potomna zachowuje fenotyp komórki macierzystej) oraz Asym Diff (obie komórki potomne mają różne fenotypy i żadna z nich nie przypomina komórki macierzystej). (E) Histogramy wkładu typów komórek zaklasyfikowanych według generacji dla MPP hodowanych z koktajlem „Diff”; n = 204 komórki i 97 rodzin. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Tabela 5: Opis liczby rodzin i komórek analizowanych dla każdego warunku eksperymentalnego (komórka macierzysta i pożywka hodowlana). Kliknij tutaj, aby pobrać tę tabelę.
Plik uzupełniający 1: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 2: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 3: Kliknij tutaj, aby pobrać ten plik.
Plik uzupełniający 4: Kliknij tutaj, aby pobrać ten plik.