$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Prezentowana funkcjonalna platforma jednokomórkowa pozwoliła na pomiar kilku parametrów. Po pierwsze, podobnie jak w przypadku standardowych technik, częstotliwość wydzielania komórek jest przedstawiona na końcu pomiaru (Rysunek 4A). Po stymulacji 1 μg/ml lipopolisacharydu (LPS) przez 6 godzin komórek jednojądrzastych krwi obwodowej (PBMC), 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ściowo wydzielanie cytokin, wygenerowano krzywe kalibracyjne ze znanymi stężeniami rekombinowanych cytokin (Rysunek 4B). Te krzywe kalibracyjne umożliwiają ilościowe określenie stężeń cytokin w kroplach w czasie. Na przykład średnie stężenie IL-6 w kropli osiągnęło plateau po 90 minutach dla PBMC stymulowanego przez LPS, podczas gdy średnie stężenie IL-1β w kropli rosło szybciej od 90 minut, co pokazuje dynamiczną rozdzielczość platformy i możliwość ekstrakcji subpopulacji komórek wydzielających określone cytokiny (Rysunek 4C). Ponieważ stężenie zmienia się między punktami pomiarowymi, możliwe jest obliczenie dynamicznych wskaźników wydzielania na cytokinę. Biorąc pod uwagę średnią szybkość wydzielania dla każdej cytokiny (Ryc. 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β wykazały wzrost szybkości wydzielania po 90 minutach pomiaru i drugi spadek po 150 minutach.
Ponadto, możliwe jest grupowanie 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 30,2% i 26,4% komórek wydzielających odpowiednio IL-6 lub TNFα, podczas gdy komórki IL-1β wydzielające pojedynczo stanowiły 68,8% wszystkich komórek wydzielających IL-1β. Dodatkowo, wpływ współwydzielania na wydzielane stężenia i szybkość wydzielania można rozwiązać (Rysunek 4F). Patrząc na komórki wydzielające IL-6, różne ilości IL-6 były wydzielane, jeśli komórki dodatkowo wytwarzały TNFα lub IL-1β. Podobnie, rozkład uśrednionych szybkości wydzielania w trakcie pomiaru różnił się statystycznie między komórkami wydzielającymi tylko IL-6 lub IL-6 wraz z TNFα (wyższe wskaźniki wydzielania) i IL-1β (niższe wskaźniki wydzielania IL-6).

Rysunek 4: Reprezentatywne wyniki wydzielania PBMC przez IL-6, TNFα i IL-1β po 6 godzinnej stymulacji 1 μg/ml LPS. (A) Procent PBMC wydzielających IL-6, TNFα i IL-1β pod koniec 4-godzinnego pomiaru. (B) Krzywe kalibracji cytokin multipleksowanych są generowane ze znanymi stężeniami cytokin rekombinowanych. Pozwala to na ilościowe określenie eksperymentów komórkowych poprzez obliczenie na podstawie wartości relokacji stężenia cytokin w kropli. Punkty dopasowano za pomocą nieliniowej jednofazowej krzywej asocjacyjnej, r2=0,9926 (IL-6), 0,9901 (TNFα), 0,9990 (IL-1β). (C) Średnie wydzielane stężenia IL-6, TNFα i IL-1β uwalniane przez wydzielanie PBMC w ciągu 4-godzinnego czasu pomiaru. (D) Średnie tempo wydzielania IL-6, TNFα i IL-1β w czasie pomiaru 4 godzin. (E) Względny odsetek komórek współwydzielających IL-6, TNFα lub IL-1β i ich kombinacje. Znormalizowane do wszystkich komórek wydzielających wykrytych dla każdej cytokiny. (F) Uśrednione stężenia IL-6 w czasie pomiaru i rozkłady średniej szybkości wydzielania (log) dla komórek wydzielających IL-6 z rozdzielczością współwydzielania (n = 383 tylko dla IL-6, n = 531 dla IL-6 + TNFα, n = 213 dla IL-6 + IL-1β i n = 143 dla IL-6 + TNFα + IL-1β). Statystyczne różnice w rozkładach szybkości wydzielania oceniano za pomocą dwustronnych, niesparowanych, nieparametrycznych testów Kołmogorowa-Smirnowa z 95% pewnością, reprezentowane są wartości p. ** (p <0,002) i **** (p <0,0001). Pełna linia reprezentuje medianę, a linia przerywana kwartyle. nkomórek ogółem = 21 866. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.
Aby uzyskać dodatkowe informacje na poziomie pojedynczej komórki, funkcja esicy może być dopasowana do punktów stężenia-czasu każdej komórki i cytokiny (Rysunek 5). Przykładowy zestaw danych stężenia w czasie dla jednej komórki i odpowiadającego jej dopasowania sigmoidalnego przedstawiono na rysunku Rysunek 5A. W tym przypadku procedura dopasowania metodą najmniejszych kwadratów daje następujące parametry: C, odpowiadający górnej wartości plateau krzywej, t50 określający ilościowo przesunięcie krzywej w czasie od zera oraz nachylenie wzgórza m, opisujące stromość wznoszącej się części krzywej esicy z wartościami stężeń 10% i 90% osiągniętymi podczas pomiaru. Z tych parametrów dopasowania można wyodrębnić niektóre deskryptory krzywych, jak wyjaśniono w kroku 7.12. uzyskując Cmax, najwyższą wartość stężenia w danych, tstart, czas rozpoczęcia wydzielania, zdefiniowany jako osiągnięcie 10% wartości stężenia na górnym platau, oraz SRlin, szybkość wydzielania w rosnącej części krzywej.
Aby sklasyfikować subpopulacje komórek, deskryptory krzywych uzyskane ze wszystkich dopasowań pojedynczych komórek zostały podzielone na trzy kategorie:Maksymalne wartości C zostały pogrupowane na niskie, średnie i wysokie dlat start na wczesny, średni i późny, a SRlin na wolne, średnie i szybkie wydzielacze. Aby zilustrować tę klasyfikację, pokazano cztery przykładowe krzywe wydzielania pojedynczych komórek i odpowiadające im deskryptory krzywych (Rysunek 5A-D), gdzie krzywa A wykazuje cechy wczesnego niskiego wydzielacza o średniej szybkości, krzywa B jest wczesnym, wolnym i wysokim wydzielaczem, krzywa C jest wczesnym szybkim wysokim wydzielaczem, a krzywa D pokazuje późne niskie wydzielanie. Ważne jest, aby pamiętać, że punkty odcięcia dla tych kryteriów są specyficzne dla komórki, cytokin i parametrów testu i muszą być dostosowane do każdego pytania badawczego. Co więcej, uwzględniono tutaj tylko wydzielanie IL-6 PBMC po stymulacji 1 μg / ml LPS przez 6 godzin, co oznacza, że większość komórek była wczesna i wysoko wydzielająca odpowiednio 80% i 79% (Figura 5E-F). Jeśli chodzi o szybkość wydzielania, zaobserwowano odpowiedź dwubiegunową, przy czym 55% komórek wydzielających IL-6 to powolne wydzielacze, a 39% to szybkie wydzielacze (Figura 5G).
Aby dokładniej scharakteryzować zachowanie wydzielania, deskryptory krzywych dla każdej komórki zostały wykreślone względem siebie i wyodrębniono różne klastry (Rysunek 5H-J). Nie podano wyraźnej korelacji międzypoczątkiem t a Cmax (Ryc. 5H): dwie największe populacje były wczesnymi niskimi wydzielaczami i wysokimi wydzielaczami, niezależnie od początku wydzielania. Biorąc pod uwagę związek międzypoczątkiem t a SRlin (Ryc. 5I), większość komórek była wczesnymi powolnymi wydzielaczami z wyraźną populacją wczesnych wysokich wydzielaczy i kilkoma wolnymi/średnimi do późnych wydzielaczy. Jeśli chodzi o korelacje SRlin i Cmax (Rysunek 5J), prawie nie występowały szybkie niskie i średnie wydzielacze, a tylko większa populacja szybkich niskich wydzielaczy. Co więcej, istniała duża populacja szybkich wydzielaczy, które nie zależały od maksymalnego zmierzonego stężenia, oraz dwie populacje wysokich wydzielaczy wydzielanych albo wolno, albo szybko. Podsumowując, można stwierdzić, że badanie związku między deskryptorami krzywych dla poszczególnych komórek daje znacznie bardziej szczegółową analizę i może potencjalnie wydobyć nowe odkrycia biologiczne z pomiarów wydzielania pojedynczych komórek.
Za pomocą analizy przedstawionej powyżej, wyodrębniliśmy dynamikę wydzielania komórek współwydzielających (Rysunek 6). Dwie przykładowe krzywe pokazują różną dynamikę współwydzielania IL-6 i TNFα z dwóch pojedynczych komórek z jednoczesnym początkiem obu cytokin (Figura 6A) lub sekwencyjnym początkiem wydzielania, przy czym IL-6 jest wydzielana jako pierwsza (Figura 6B). Aby sklasyfikować wszystkie komórki współwydzielające, zdefiniowano opóźnienie wydzielania wynoszące 60 minut, gdzie wszystkie komórki rozpoczynające wydzielanie w tym zakresie są uważane za równoczesne wydzielacze, a wszystkie komórki z dłuższymi opóźnieniami są uważane za sektory sekwencyjne. Analiza ta pozwoliła również zaobserwować, która cytokina została wydzielona jako pierwsza. W przypadku IL-6 i TNFα głównie jednoczesne współwydzielanie zaobserwowano w 76% komórek (Figura 6C), podczas gdy w przypadku IL-6 i IL-1β sekwencyjne współwydzielanie zaobserwowano w 86% komórek, przy czym IL-6 była pierwszą cytokiną, która była wydzielana w większości przypadków (Figura 6D).
Patrząc na czas rozpoczęcia wydzielania dla różnych cytokin dla wszystkich indywidualnych komórek współwydzielających, nie zaobserwowano wyraźnej korelacji między czasami rozpoczęcia wydzielania w przeprowadzonych eksperymentach. W przypadku współwydzielania IL-6 i TNFα (Rysunek 6E) występowało większe pionowe skupisko około 0 minut, odpowiadające komórkom współwydzielającym częściej począwszy od IL-6. W przypadku jednoczesnego wydzielania IL-6 i IL-1β (Figura 6F), większość komórek zaczęła wydzielać IL-6 na początku pomiaru, podczas gdy IL-1β była głównie wydzielana później. Podsumowując, przedstawiona analiza umożliwiła identyfikację różnych subpopulacji wydzielaczy i złożonej dynamiki współwydzielania cytokin.

Rysunek 5: Szczegółowa analiza różnych wzorców dynamicznych wydzielania dla krzywych pojedynczych komórek wydzielających IL-6. (A) Reprezentatywne dane dotyczące stężenia cytokin w pojedynczej komórce w czasie pomiaru z dopasowaną krzywą esicy i wyekstrahowanymi parametrami. (B-D) Trzy przykładowe krzywe stężenia cytokin jednokomórkowych dla różnych typów wydzielaczy cytokin znalezione dla wydzielania IL-6 po stymulacji LPS. (E-G) Procentowy udział komórek wydzielających IL-6, które są klasyfikowane do różnych typów wydzielaczy według następujących kryteriów (n = 633): E. Cmax: niskie <5 nM, wysokie >19,5 nM, F. tpoczątek: wcześnie <30 min, późno >120 min, G. SRlin: wolno <250 cząsteczek/s, szybko >750 molekuł/s. (H-J) Zależność między trzema deskryptorami krzywej wydzielania Cmax, tstart i SRlin dla każdej pojedynczej komórki (n=633). Duża populacja przy Cmax=20nM wynika z osiągnięcia górnej granicy wykrywalności testu. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 6: Ekstrakcja wzorców współwydzielania z krzywych koncentracji pojedynczej komórki. (A-B) Reprezentatywne krzywe stężenia dla pojedynczych komórek współwydzielających IL-6 i TNFα (A) odpowiednio jednocześnie i (B) sekwencyjnie. (C-D) Odsetek komórek wykazujących jednoczesne i sekwencyjne współwydzielanie IL-6 i TNFα (n = 249) lub odpowiednio IL-6 i IL-1β (n = 72). Wydzielanie sekwencyjne definiuje się jako opóźnienie między rozpoczęciem wydzielania cytokin większe niż 60 minut. Kolory wskazują, które z cytokin zaczęły wydzielać jako pierwsze. (E-F) Zależność między czasami rozpoczęcia wydzielania dla różnych cytokin dla każdej komórki wydzielającej (nIL6-TNFα=249, nIL6-IL1β=72). Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.