$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Die vorgestellte funktionale Einzelzellplattform ermöglichte die Messung mehrerer Parameter. Zunächst wird ähnlich wie bei Standardtechniken am Ende der Messung die Häufigkeit der sezernierenden Zellen dargestellt (Abbildung 4A). Nach der Stimulation mit 1 μg/ml Lipopolysaccharid (LPS) für 6 Stunden von mononukleären Zellen des peripheren Blutes (PBMC) sezernierten 5,81% der Zellen IL-6 (n= 1270), 4,55% TNFα (n= 995) und 6,06% IL-1β (n= 1326).
Um die Zytokinsekretion zu quantifizieren, wurden Kalibrierkurven mit bekannten Konzentrationen rekombinanter Zytokine erstellt (Abbildung 4B). Diese Kalibrierkurven ermöglichen die Quantifizierung der Zytokinkonzentrationen in Tröpfchen über die Zeit. Beispielhaft erreichte die durchschnittliche IL-6-Konzentration in Tröpfchen nach 90 min für LPS-stimulierte PBMC ein Plateau, während die durchschnittliche IL-1β-Konzentration in Tröpfchen ab 90 min schneller anstieg, was die dynamische Auflösung der Plattform und die Möglichkeit zeigt, Zellsubpopulationen zu extrahieren, die spezifische Zytokine sezernieren (Abbildung 4C). Da sich die Konzentration zwischen den Messpunkten ändert, ist es möglich, dynamische Sekretionsraten pro Zytokin zu berechnen. Unter Berücksichtigung der durchschnittlichen Sekretionsrate für jedes Zytokin (Abbildung 4D) zeigten IL-6-sezernierende Zellen eine konstante Abnahme der durchschnittlichen Sekretionsrate, während TNFα- und IL-1β-sezernierende Zellen beide einen Anstieg der Sekretionsrate nach 90 min Messzeit und einen zweiten Rückgang nach 150 min zeigten.
Darüber hinaus ist es möglich, Zellen in Abhängigkeit von den sezernierten und co-sezernierten Zytokinen in Subpopulationen zu gruppieren (Abbildung 4E). Hier werden IL-6 und TNFα von 30,2 % bzw. 26,4 % der Zellen, die IL-6 bzw. TNFα sezernieren, einfach sezerniert, während einfach sezernierende IL-1β-Zellen 68,8 % aller IL-1β-sezernierenden Zellen ausmachten. Darüber hinaus können die Auswirkungen der Co-Sekretion auf die sezernierten Konzentrationen und Sekretionsraten geklärt werden (Abbildung 4F). Bei der Betrachtung von IL-6-sezernierenden Zellen wurden unterschiedliche Mengen an IL-6 sezerniert, wenn die Zellen zusätzlich TNFα oder IL-1β produzierten. In ähnlicher Weise unterschied sich die Verteilung der gemittelten Sekretionsraten über die Messung statistisch zwischen den Zellen, die neben TNFα (höhere Sekretionsraten) und IL-1β (niedrigere IL-6-Sekretionsraten) nur IL-6 oder IL-6 sezernierten.

Abbildung 4: Repräsentative Ergebnisse der IL-6-, TNFα- und IL-1β-sezernierenden PBMC nach 6-stündiger Stimulation mit 1 μg/ml LPS. (A) Prozentualer Anteil der PBMC, die IL-6, TNFα und IL-1β am Ende der 4-stündigen Messung sezerniert. (B) Multiplex-Zytokin-Kalibrierungskurven werden mit bekannten Konzentrationen rekombinanter Zytokine erstellt. Dies ermöglicht die Quantifizierung von Zellexperimenten, indem aus dem Verlagerungswert die Zytokinkonzentration im Tröpfchen berechnet wird. Die Punkte wurden mit einer nichtlinearen einphasigen Assoziationskurvenanpassung angepasst, r2=0,9926 (IL-6), 0,9901 (TNFα), 0,9990 (IL-1β). (C) Durchschnittliche sekretierte Konzentrationen von IL-6, TNFα und IL-1β, die durch Seernierung von PBMC über die Messzeit von 4 h freigesetzt wurden. (D) Durchschnittliche Sekretionsraten von IL-6, TNFα und IL-1β über die Messzeit von 4 Stunden. (E) Relativer Prozentsatz der mitsezernierenden Zellen, die IL-6, TNFα oder IL-1β und Kombinationen davon sezernieren. Normalisiert auf alle sekretierenden Zellen, die für jedes Zytokin nachgewiesen wurden. (F) Gemittelte IL-6-Konzentrationen über die Messzeit und Verteilungen der durchschnittlichen Sekretionsrate (log) für IL-6-sezernierende Zellen mit Co-Sekretionsauflösung (n=383 für nur IL-6, n=531 für IL-6 + TNFα, n= 213 für IL-6 + IL-1β und n=143 für IL-6+TNFα+IL-1β). Statistische Unterschiede in den Verteilungen der Sekretionsrate wurden mit zweiseitigen, ungepaarten, nichtparametrischen Kolmogorov-Smirnov-Tests mit einer Konfidenz von 95% bewertet, wobei der p-Wert dargestellt ist. ** (p <0,002) und **** (p <0,0001). Die volle Linie stellt den Median und die gestrichelte Linie die Quartile dar. nZellen insgesamt = 21 866. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.
Um zusätzliche Informationen auf Einzelzellebene zu extrahieren, kann eine Sigmoidfunktion an die Konzentrations-Zeitpunkte jeder Zelle und jedes Zytokins angepasst werden (Abbildung 5). Ein beispielhafter Datensatz der Konzentration über die Zeit für eine Zelle und die entsprechende sigmoidale Anpassung ist in Abbildung 5A dargestellt. Hier liefert das Anpassungsverfahren der kleinsten Quadrate die folgenden Parameter: C, entsprechend dem oberen Plateauwert der Kurve, t50 , der die zeitliche Verschiebung der Kurve von Null quantifiziert, und die Hügelsteigung m, die die Steilheit des ansteigenden Teils der Sigmakurve mit 10 % und 90 % Konzentrationswerten beschreibt, die während der gesamten Messung erreicht wurden. Aus diesen Anpassungsparametern können einige Kurvendeskriptoren extrahiert werden, wie in Schritt 7.12 beschrieben. Dies ergibt Cmax, den höchsten Konzentrationswert der Daten, tstart, die Startzeit der Sekretion, definiert als das Erreichen von 10 % des Konzentrationswerts des oberen Plateaus, und SRlin, die Sekretionsrate während des ansteigenden Teils der Kurve.
Um Zellsubpopulationen zu klassifizieren, wurden die Kurvendeskriptoren, die aus allen Einzelzellanpassungen gewonnen wurden, in jeweils drei Kategorien eingeteilt: DieCmax-Werte wurden in niedrig, mittel und hoch gruppiert,wobei die t-Werte in frühe, mittlere und späteSR-Lin in langsame, mittlere und schnelle Sekretoren eingeteilt wurden. Um diese Klassifizierung zu veranschaulichen, werden vier exemplarische Einzelzellsekretionskurven und ihre entsprechenden Kurvendeskriptoren gezeigt (Abbildung 5A-D), wobei Kurve A die Eigenschaften eines frühen niedrigen Sekretors mittlerer Rate aufweist, Kurve B ein früher, langsamer und hoher Sekretor, Kurve C ein früher schneller Sekretor mit hohem Sekretor und Kurve D eine späte niedrige Sekretion. Es ist wichtig zu beachten, dass die Cutoffs für diese Kriterien zell-, zytokin- und assayparameterspezifisch sind und für jede Forschungsfrage angepasst werden müssen. Darüber hinaus wurde hier nur die IL-6-Sekretion von PBMC nach 1 μg/mL LPS-Stimulation für 6 h berücksichtigt, was bedeutet, dass die meisten Zellen mit 80% bzw. 79% frühe und hohe Sekretoren waren (Abbildung 5E-F). In Bezug auf die Sekretionsrate wurde eine bipolare Reaktion beobachtet, wobei 55% der IL-6-sezernierenden Zellen langsame Sekretoren und 39% schnelle Sekretoren sind (Abbildung 5G).
Um das Sekretionsverhalten weiter zu charakterisieren, wurden die Kurvendeskriptoren für jede Zelle gegeneinander aufgetragen und verschiedene Cluster extrahiert (Abbildung 5H-J). Es gibt keine klare Korrelation zwischent-Start undC-max (Abbildung 5H): Die beiden größten Populationen waren frühe niedrige Sekretoren und hohe Sekretoren unabhängig vom Sekretionsstart. Unter Berücksichtigung der Beziehung zwischent-Start und SRlin (Abbildung 5I) waren die meisten Zellen frühe langsame Sekretoren mit einer klaren Population von frühen hohen Sekretoren und wenigen langsamen/mittleren bis späten Sekretoren. Hinsichtlich der Korrelationen von SRlin und Cmax (Abbildung 5J) waren fast keine schnellen niedrigen bis mittleren Sekretoren vorhanden, mit nur einer größeren Population von schnellen niedrigen Sekretoren. Darüber hinaus gab es eine große Population von schnellen Sekretoren, die nicht von der maximal gemessenen Konzentration abhingen, und zwei Populationen von hohen Sekretoren, die entweder langsam oder schnell sezernierten. Zusammenfassend lässt sich sagen, dass die Untersuchung des Zusammenhangs zwischen den Kurvendeskriptoren für einzelne Zellen zu einer wesentlich detaillierteren Analyse führt und möglicherweise neue biologische Erkenntnisse aus Einzelzellsekretionsmessungen extrahieren kann.
Mit der oben vorgestellten Analyse haben wir die Sekretionsdynamik von co-sekretierenden Zellen extrahiert (Abbildung 6). Zwei Beispielkurven zeigen eine unterschiedliche Dynamik der Co-Sekretion von IL-6 und TNFα aus zwei Einzelzellen mit einem gleichzeitigen Start beider Zytokine (Abbildung 6A) oder einem sequentiellen Sekretionsstart, wobei IL-6 zuerst sezerniert wird (Abbildung 6B). Zur Klassifizierung aller co-sekretierenden Zellen wurde eine Sekretionsverzögerung von 60 min definiert, wobei alle Zellen, die innerhalb dieses Bereichs mit der Sekretion beginnen, als simultane Sekretoren und alle Zellen mit längerer Verzögerung als sequentielle Sekretoren betrachtet werden. Diese Analyse ermöglichte auch die Beobachtung, welches Zytokin zuerst sezerniert wurde. Für IL-6 und TNFα wurde in 76 % der Zellen hauptsächlich eine gleichzeitige Co-Sekretion beobachtet (Abbildung 6C), während für IL-6 und IL-1β in 86 % der Zellen eine sequentielle Co-Sekretion beobachtet wurde, wobei IL-6 in den meisten Fällen das erste Zytokin war, das sezerniert wurde (Abbildung 6D).
Betrachtet man den Startzeitpunkt der Sekretion für die verschiedenen Zytokine für alle einzelnen co-sekretierenden Zellen, so konnte in den durchgeführten Experimenten keine eindeutige Korrelation zwischen den Startzeiten der Sekretion beobachtet werden. Für die Co-Sekretion von IL-6 und TNFα (Abbildung 6E) war ein größerer vertikaler Cluster um 0 min vorhanden, was den co-sezernierenden Zellen entspricht, die häufiger mit IL-6 beginnen. Für die Co-Sekretion von IL-6 und IL-1β (Abbildung 6F) begannen die meisten Zellen zu Beginn der Messung mit der Sekretion von IL-6, während IL-1β hauptsächlich später sezerniert wurde. Zusammenfassend lässt sich sagen, dass die hier vorgestellte Analyse die Identifizierung verschiedener Sekretor-Subpopulationen und komplexer Zytokin-Co-Sekretionsdynamiken ermöglichte.

Abbildung 5: Detaillierte Analyse verschiedener dynamischer Sekretionsmuster für einzelne IL-6-sezernierende Zellkurven. (A) Repräsentative Einzelzell-Zytokin-Konzentrationsdaten über die Messzeit mit der angepassten Sigmoidkurve und den extrahierten Parametern. (B-D) Drei beispielhafte Einzelzell-Zytokin-Konzentrationskurven für die verschiedenen Zytokinsekretortypen, die für die IL-6-Sekretion nach LPS-Stimulation gefunden wurden. (E-G) Prozentuale Anteile von IL-6-sezernierenden Zellen, die mit folgenden Kriterien (n=633) in die verschiedenen Sekretortypen eingeteilt werden: E. Cmax: niedrig <5 nM, hoch >19,5 nM, F. tStart: früh <30 min, spät >120 min, G. SRlin: langsam <250 Moleküle/s, schnell >750 Moleküle/s. (H-J) Beziehung zwischen den drei Sekretionskurvendeskriptoren Cmax, tstart und SRlin für jede einzelne Zelle (n=633). Die große Population bei Cmax=20nM resultiert aus dem Erreichen der oberen Nachweisgrenze des Assays. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.

Abbildung 6: Extraktion von Co-Sekretionsmustern aus Einzelzell-Konzentrationskurven. (A-B) Repräsentative Konzentrationskurven für Einzelzellen, die IL-6 und TNFα (A) gleichzeitig bzw. (B) sequentiell co-sezernieren. (C-D) Prozentsatz der Zellen, die eine gleichzeitige und sequentielle Co-Sekretion von IL-6 und TNFα (n=249) bzw. IL-6 und IL-1β (n=72) aufweisen. Die sequentielle Sekretion wird durch die Verzögerung zwischen den Zytokinsekretionsstarts von mehr als 60 Minuten definiert. Farben zeigen an, welches der Zytokine zuerst mit der Sekretion begonnen hat. (E-F) Beziehung zwischen den Startzeiten der Sekretion für die verschiedenen Zytokine für jede sekretierende Zelle (nIL6-TNFα=249, nIL6-IL1β=72). Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzuzeigen.