Przedstawiona funkcjonalna platforma do analizy pojedynczych komórek umożliwiła pomiar kilku parametrów. Po pierwsze, analogicznie do standardowych technik, na końcu pomiaru przedstawiono częstotliwość występowania komórek wydzielniczych (Rycina 4A). Po stymulacji mononuklearnych komórek krwi obwodowej (PBMC) lipopolisacharydem (LPS) w stężeniu 1 µg/mL przez 6 h, 5,81% komórek wydzielało IL-6 (n= 1270), 4,55% TNFα (n= 995) i 6,06% IL-1β (n= 1326).
Aby określić ilość wydzielanych cytokin, opracowano krzywe kalibracyjne z zastosowaniem znanych stężeń cytokin rekombinowanych (Rysunek 4B). Krzywe kalibracyjne te umożliwiają ilościowe oznaczenie stężeń cytokin wewnątrz kropel w czasie. Przykładowo, średnie stężenie IL-6 wewnątrz kropel osiągnęło plateau po 90 min w przypadku PBMC stymulowanych LPS, podczas gdy średnie stężenie IL-1β wewnątrz kropel rosło szybciej od 90 min, co wykazuje dynamiczną rozdzielczość platformy oraz możliwość wydzielania subpopulacji komórek sekrecyjnych dla konkretnych cytokin (Rysunek 4C). Ponieważ stężenia zmieniają się pomiędzy punktami pomiarowymi, możliwe jest obliczenie dynamicznych wskaźników wydzielania dla poszczególnych cytokin. Biorąc pod uwagę średnią szybkość wydzielania każdej cytokiny (Rysunek 4D), komórki wydzielające IL-6 wykazywały stały spadek średniej szybkości wydzielania, podczas gdy komórki wydzielające TNFα i IL-1β wykazywały wzrost szybkości wydzielania po 90 min pomiaru oraz wtórny spadek po 150 min.
Ponadto możliwe jest pogrupowanie komórek w subpopulacje w zależności od wydzielanych i współwydzielanych cytokin (Rysunek 4E). W tym przypadku IL-6 i TNFα są wydzielane pojedynczo przez odpowiednio 30,2% i 26,4% komórek wydzielających IL-6 lub TNFα, natomiast komórki wydzielające pojedynczo IL-1β stanowiły 68,8% wszystkich komórek wydzielających IL-1β. Dodatkowo można przeanalizować wpływ współwydzielania na stężenia wydzielanych cytokin oraz szybkość ich sekrecji (Rysunek 4F). Analizując komórki wydzielające IL-6, stwierdzono, że wydzielano różne ilości IL-6 w zależności od tego, czy komórki dodatkowo produkowały TNFα lub IL-1β. Podobnie rozkład średnich szybkości sekrecji w czasie pomiaru różnił się statystycznie pomiędzy komórkami wydzielającymi tylko IL-6 a komórkami wydzielającymi IL-6 wraz z TNFα (wyższe szybkości sekrecji) i IL-1β (niższe szybkości sekrecji IL-6).

Rycina 4: Reprezentatywne wyniki sekrecji IL-6, TNFα i IL-1β przez PBMC po 6 h stymulacji 1 µg/mL LPS. (A) Procent PBMC sekrecyjnych dla IL-6, TNFα i IL-1β na koniec 4 h pomiaru. (B) Multiplexowe krzywe kalibracyjne cytokin wygenerowane przy użyciu znanych stężeń cytokin rekombinowanych. Pozwala to na ilościowe określenie wyników eksperymentów komórkowych poprzez obliczenie stężenia cytokiny wewnątrz kropelki na podstawie wartości relokacji. Punkty dopasowano za pomocą nieliniowej jednofazowej krzywej asocjacyjnej, r2=0,9926 (IL-6), 0,9901 (TNFα), 0,9990 (IL-1β). (C) Średnie sekreowane stężenia IL-6, TNFα i IL-1β uwalniane przez sekrecyjne PBMC w ciągu 4 h czasu pomiaru. (D) Średnie szybkości sekrecji IL-6, TNFα i IL-1β w ciągu 4 h czasu pomiaru. (E) Względny procent komórek współsekrecyjnych wydzielających IL-6, TNFα lub IL-1β oraz ich kombinacje. Znormalizowano do wszystkich komórek sekrecyjnych wykrytych dla każdej cytokiny. (F) Średnie stężenia IL-6 w czasie pomiaru oraz rozkłady średniej szybkości sekrecji (log) dla komórek sekrecyjnych IL-6 z rozdzielczością współsekrecji (n=383 dla samych IL-6, n=531 dla IL-6 + TNFα, n=213 dla IL-6 + IL-1β oraz n=143 dla IL-6+TNFα+IL-1β). Różnice statystyczne w rozkładach szybkości sekrecji oceniono za pomocą dwustronnych, niezależnych, nieparametrycznych testów Kołmogorowa-Smirnowa z 95% poziomem ufności; przedstawiono wartości p. ** (p <0,002) i **** (p <0,0001). Linia ciągła reprezentuje medianę, a linia przerywana kwartyle. ncałkowita liczba komórek = 21 866. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
Aby uzyskać dodatkowe informacje na poziomie pojedynczych komórek, do punktów stężenia w czasie dla każdej komórki i cytokiny można dopasować funkcję sigmoidalną (Rycina 5). Przykładowy zestaw danych stężenia w czasie dla jednej komórki oraz odpowiadające mu dopasowanie sigmoidalne przedstawiono na Rycynie 5A. W tym przypadku procedura dopasowania metodą najmniejszych kwadratów pozwala uzyskać następujące parametry: C, odpowiadające wartości górnego plateau krzywej, t50 określający przesunięcie czasowe krzywej od zera oraz nachylenie Hill’a m, opisujące stromość wzrostowej części krzywej sigmoidalnej, przy których osiągnięto wartości stężeń 10% i 90% w trakcie pomiaru. Z tych parametrów dopasowania można wyodrębnić niektóre deskryptory krzywej, zgodnie z wyjaśnieniami w kroku 7.12, otrzymując Cmax (najwyższą wartość stężenia w danych), tstart (czas rozpoczęcia sekrecji, zdefiniowany jako osiągnięcie 10% wartości stężenia górnego plateau) oraz SRlin (szybkość sekrecji podczas wzrostowej części krzywej).
Aby zaklasyfikować subpopulacje komórek, deskryptory krzywych uzyskane ze wszystkich dopasowań dla pojedynczych komórek podzielono na trzy kategorie dla każdego parametru: wartości Cmax pogrupowano jako niskie, średnie i wysokie, tstart na wczesne, średnie i późne, a SRlin na sekretory wolne, średnie i szybkie. Aby zilustrować tę klasyfikację, przedstawiono cztery przykładowe krzywe sekrecji pojedynczych komórek oraz odpowiadające im deskryptory krzywych (Rysunek 5A-D), gdzie krzywa A wykazuje charakterystykę wczesnego sekretora o niskim poziomie i średnim tempie, krzywa B reprezentuje wczesnego, wolnego sekretora o wysokim poziomie, krzywa C wczesnego, szybkiego sekretora o wysokim poziomie, a krzywa D wykazuje późną, niską sekrecję. Należy zauważyć, że wartości progowe dla tych kryteriów są specyficzne dla rodzaju komórek, cytokiny i parametrów testu i muszą być dostosowane do każdego pytania badawczego. Ponadto, w analizie tej uwzględniono jedynie sekrecję IL-6 przez PBMC po stymulacji 1 µg/mL LPS przez 6 h, co oznacza, że większość komórek była wczesnymi i wysokimi sekretorami, odpowiednio w 80% i 79% przypadków (Rysunek 5E-F). W odniesieniu do tempa sekrecji zaobserwowano odpowiedź bimodalną, w której 55% komórek wydzielających IL-6 to sekretory wolne, a 39% to sekretory szybkie (Rysunek 5G).
Aby dalej scharakteryzować zachowanie sekrecyjne, deskryptory krzywych dla każdej komórki naniesiono na wzajemne wykresy i wyodrębniono różne klastry (Rycyna 5H-J). Nie stwierdzono wyraźnej korelacji między tstart a Cmax (Rycina 5H): dwie największe populacje stanowiły wczesne komórki o niskiej sekrecji oraz komórki o wysokiej sekrecji niezależnie od momentu rozpoczęcia procesu. Analizując relację między tstart a SRlin (Rycina 5I), większość komórek okazała się wczesnymi komórkami o powolnej sekrecji, przy wyraźnej populacji wczesnych komórek o wysokiej sekrecji oraz niewielkiej liczbie komórek o powolnej/średniej sekrecji rozpoczynającej się późno. W odniesieniu do korelacji między SRlin a Cmax (Rycina 5J), niemal nie wystąpiły szybkie komórki o niskiej do średniej sekrecji, z wyjątkiem nieco liczniejszej populacji szybkich komórek o niskiej sekrecji. Ponadto odnotowano dużą populację szybkich komórek sekrecyjnych, której liczebność nie zależała od maksymalnego zmierzonego stężenia, oraz dwie populacje komórek o wysokiej sekrecji, które sekretowały albo powolnie, albo szybko. Podsumowując, można wywnioskować, że badanie zależności między deskryptorami krzywych dla poszczególnych komórek pozwala na znacznie bardziej szczegółową analizę i może potencjalnie przynieść nowe odkrycia biologiczne w pomiarach sekrecji na poziomie pojedynczych komórek.
Dzięki przedstawionej powyżej analizie wyekstrahowano dynamikę sekrecji komórek współsekretujących (Rysunek 6). Dwie przykładowe krzywe pokazują różną dynamikę współsekrecji IL-6 i TNFα z dwóch pojedynczych komórek, z jednoczesnym rozpoczęciem sekrecji obu cytokin (Rysunek 6A) lub sekrecją sekwencyjną, w której pierwsza była sekretowana IL-6 (Rysunek 6B). Aby sklasyfikować wszystkie komórki współsekretujące, zdefiniowano opóźnienie sekrecji wynoszące 60 min, gdzie wszystkie komórki rozpoczynające sekrecję w tym przedziale uznano za sekretory jednoczesne, a wszystkie komórki z dłuższym opóźnieniem za sekretory sekwencyjne. Analiza ta pozwoliła również ustalić, która cytokina była sekretowana jako pierwsza. W przypadku IL-6 i TNFα w 76% komórek zaobserwowano głównie jednoczesną współsekrecję (Rysunek 6C), natomiast dla IL-6 i IL-1β współsekrecję sekwencyjną odnotowano w 86% komórek, przy czym w większości przypadków pierwszą sekretowaną cytokiną była IL-6 (Rysunek 6D).
Analizując czas rozpoczęcia sekrecji różnych cytokin dla wszystkich poszczególnych komórek współsekrecyjnych, w przeprowadzonych eksperymentach nie zaobserwowano wyraźnej korelacji między czasami rozpoczęcia sekrecji. W przypadku współsekrecji IL-6 i TNFα (Rysunek 6E) widoczna była większa pionowa grupa wokół 0 min, co odpowiadało częstszemu rozpoczynaniu procesu przez komórki od sekrecji IL-6. W przypadku współsekrecji IL-6 i IL-1β (Rysunek 6F) większość komórek zaczęła wydzielać IL-6 około rozpoczęcia pomiaru, podczas gdy IL-1β była wydzielana głównie później. Podsumowując, zaprezentowana tutaj analiza umożliwiła identyfikację różnych subpopulacji komórek wydzielniczych oraz złożonej dynamiki współsekrecji cytokin.

Rysunek 5: Szczegółowa analiza różnych wzorców dynamiki sekrecji dla krzywych pojedynczych komórek wydzielających IL-6. (A) Reprezentatywne dane dotyczące stężenia cytokin w pojedynczych komórkach w czasie pomiaru wraz z dopasowaną krzywą sigmoidalną i wyznaczonymi parametrami. (B-D) Trzy przykładowe krzywe stężenia cytokin w pojedynczych komórkach dla różnych typów wydzielników cytokin zidentyfikowanych w przypadku sekrecji IL-6 po stymulacji LPS. (E-G) Procenty komórek wydzielających IL-6 zaklasyfikowanych do różnych typów wydzielników na podstawie następujących kryteriów (n=633): E. Cmax: niskie <5 nM, wysokie >19,5 nM, F. tstart: wczesne <30 min, późne >120 min, G. SRlin: powolne <250 cząsteczek/s, szybkie >750 cząsteczek/s. (H-J) Zależność między trzema deskryptorami krzywej sekrecji Cmax, tstart i SRlin dla każdej pojedynczej komórki (n=633). Duża populacja przy Cmax=20 nM wynika z osiągnięcia górnej granicy detekcji testu. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 6: Wyodrębnianie wzorców współwydzielania z krzywych stężeń dla pojedynczych komórek. (A-B) Reprezentatywne krzywe stężeń dla pojedynczych komórek współwydzielających IL-6 i TNFα odpowiednio (A) symultanicznie oraz (B) sekwencyjnie. (C-D) Procent komórek wykazujących symultaniczne i sekwencyjne współwydzielanie odpowiednio IL-6 i TNFα (n=249) oraz IL-6 i IL-1β (n=72). Wydzielanie sekwencyjne zdefiniowano jako opóźnienie rozpoczęcia wydzielania cytokin przekraczające 60 min. Kolory wskazują, która z cytokin zaczęła być wydzielana jako pierwsza. (E-F) Zależność między czasami rozpoczęcia wydzielania różnych cytokin dla każdej wydzielającej komórki (nIL6-TNFα=249, nIL6-IL1β=72). Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.