$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Het gepresenteerde functionele single-cell platform maakte het mogelijk om verschillende parameters te meten. Ten eerste, en net als bij standaardtechnieken, wordt de frequentie van de afscheidende cellen aan het einde van de meting weergegeven (Figuur 4A). Na de stimulatie met 1 μg/ml lipopolysaccharide (LPS) gedurende 6 uur van perifere mononucleaire bloedcellen (PBMC), scheidde 5,81% van de cellen IL-6 (n= 1270), 4,55% TNFα (n= 995) en 6,06% IL-1β (n= 1326) uit.
Om de cytokinesecretie te kwantificeren, werden kalibratiecurven gegenereerd met bekende concentraties van recombinante cytokines (Figuur 4B). Deze kalibratiecurven maken het mogelijk om de cytokineconcentraties in de druppel in de loop van de tijd te kwantificeren. Voorbeeldig is dat de gemiddelde IL-6-concentratie in de druppel na 90 minuten een plateau bereikte voor LPS-gestimuleerde PBMC, terwijl de gemiddelde in-druppel IL-1β sneller toenam vanaf 90 minuten, wat de dynamische resolutie van het platform weergeeft en de mogelijkheid om celsubpopulaties te extraheren die specifieke cytokines afscheiden (Figuur 4C). Omdat de concentratie tussen de meetpunten verandert, is het mogelijk om dynamische secretiesnelheden per cytokine te berekenen. Rekening houdend met de gemiddelde secretiesnelheid voor elk cytokine (figuur 4D), vertoonden IL-6-secreterende cellen een constante afname van de gemiddelde secretiesnelheid, terwijl TNFα- en IL-1β-secreterende cellen beide een toename van de secretiesnelheid vertoonden na een meettijd van 90 minuten en een tweede afname na 150 minuten.
Bovendien is het mogelijk om cellen te clusteren in subpopulaties, afhankelijk van de uitgescheiden en gelijktijdig uitgescheiden cytokines (Figuur 4E). Hier worden IL-6 en TNFα enkelvoudig uitgescheiden door respectievelijk 30,2% en 26,4% van de cellen die IL-6 of TNFα afscheiden, terwijl enkelvoudige afscheidende IL-1β-cellen 68,8% van alle IL-1β-uitscheidende cellen uitmaakten. Bovendien kunnen de effecten van co-secretie op de uitgescheiden concentraties en secretiesnelheden worden opgelost (Figuur 4F). Door naar IL-6-afscheidende cellen te kijken, werden verschillende hoeveelheden IL-6 uitgescheiden als de cellen bovendien TNFα of IL-1β produceerden. Evenzo verschilde de verdeling van de gemiddelde secretiesnelheden over de meting statistisch tussen de cellen die alleen IL-6 of IL-6 afscheidden naast TNFα (hogere secretiesnelheden) en IL-1β (lagere IL-6-secretiesnelheden).

Figuur 4: Representatieve resultaten van IL-6, TNFα en IL-1β uitscheidende PBMC na 6 uur stimulatie met 1 μg/ml LPS. (A) Percentage PBMC dat IL-6, TNFα en IL-1β uitscheidt aan het einde van de 4 uurs meting. (B) Gemultiplexte cytokinekalibratiecurven worden gegenereerd met bekende concentraties van recombinante cytokines. Dit maakt het mogelijk om celexperimenten te kwantificeren door uit de verplaatsingswaarde de cytokineconcentratie in de druppel te berekenen. Punten werden aangebracht met behulp van een niet-lineaire eenfasige associatiecurve-fit, r2=0,9926 (IL-6), 0,9901 (TNFα), 0,9990 (IL-1β). (C) Gemiddelde uitgescheiden concentraties van IL-6, TNFα en IL-1β die vrijkomen door PBMC uit te scheiden gedurende de meettijd van 4 uur. (D) Gemiddelde secretiesnelheden van IL-6, TNFα en IL-1β gedurende de meettijd van 4 uur. (E) Relatief percentage gelijktijdig afscheidende cellen die IL-6, TNFα of IL-1β afscheiden en combinaties daarvan. Genormaliseerd naar alle afscheidende cellen die voor elk cytokine zijn gedetecteerd. (F) Gemiddelde IL-6-concentraties gedurende de meettijd en gemiddelde secretiesnelheid (log) verdelingen voor IL-6-secreterende cellen met co-secretieresolutie (n=383 alleen voor IL-6, n=531 voor IL-6 + TNFα, n= 213 voor IL-6 + IL-1β en n=143 voor IL-6+TNFα+IL-1β). Statistische verschillen in secretiesnelheidsverdelingen werden beoordeeld met behulp van tweezijdige, ongepaarde, niet-parametrische Kolmogorov-Smirnov-tests met een betrouwbaarheid van 95%, de p-waarde wordt weergegeven. ** (p <0,002) en **** (p <0,0001). De volledige lijn geeft de mediaan weer en de stippellijn de kwartielen. ntotaal aantal cellen = 21 866. Klik hier om een grotere versie van deze figuur te bekijken.
Om aanvullende informatie op het niveau van een enkele cel te extraheren, kan een sigmoïde functie worden aangepast aan de concentratietijdpunten van elke cel en cytokine (Figuur 5). Een voorbeeldige dataset over concentratie in de tijd voor één cel en de bijbehorende sigmoïdale pasvorm wordt weergegeven in figuur 5A. Hier levert de procedure voor het aanpassen van de kleinste kwadraten de volgende parameters op: C, die overeenkomt met de bovenste plateauwaarde van de curve, t50 die de tijdgewijze verschuiving van de curve van nul kwantificeert, en de Hill-helling m, die de steilheid van het stijgende deel van de sigmoïde curve beschrijft met concentratiewaarden van 10% en 90% die tijdens de meting worden bereikt. Uit deze fitparameters kunnen enkele curve-descriptoren worden geëxtraheerd, zoals uitgelegd in stap 7.12. wat de Cmax oplevert, de hoogste concentratiewaarde van de gegevens, tstart, de starttijd van de secretie, gedefinieerd als het bereiken van 10% van de bovenste plateauconcentratiewaarde, en SRlin, de secretiesnelheid tijdens het stijgende deel van de curve.
Om celsubpopulaties te classificeren, werden de curvedescriptoren verkregen uit alle eencellige fits elk in drie categorieën ingedeeld: Cmax-waarden werden gegroepeerd in laag, gemiddeld en hoog voor tstart in vroeg, gemiddeld en laat en SRlin in langzame, gemiddelde en snelle secretoren. Om deze classificatie te illustreren, worden vier voorbeeldige eencellige secretiecurven en de bijbehorende curvedescriptoren getoond (Figuur 5A-D), waarbij curve A de kenmerken vertoont van een vroege lage secretor met een gemiddelde snelheid, curve B een vroege, langzame en hoge secretor, curve C een vroege snelle hoge secretor en curve D late lage secretie vertoont. Het is belangrijk op te merken dat de cutoffs voor deze criteria cel-, cytokine- en testparameterspecifiek zijn en voor elke onderzoeksvraag moeten worden aangepast. Bovendien werd hier alleen rekening gehouden met IL-6-secretie van PBMC na 1 μg/ml LPS-stimulatie gedurende 6 uur, wat betekent dat de meeste cellen vroege en hoge secretoren waren met respectievelijk 80% en 79% (Figuur 5E-F). Wat de secretiesnelheid betreft, werd een bipolaire respons waargenomen waarbij 55% van de IL-6-secreterende cellen langzame secretoren zijn en 39% snelle secretoren (Figuur 5G).
Om het secretiegedrag verder te karakteriseren, werden de curvedescriptoren voor elke cel tegen elkaar uitgezet en werden verschillende clusters geëxtraheerd (Figuur 5H-J). Er wordt geen duidelijke correlatie gegeven tussen tstart en Cmax (Figuur 5H): de twee grootste populaties waren vroege lage secretoren en hoge secretoren onafhankelijk van de secretiestart. Rekening houdend met de relatie tussen tstart en SRlin (Figuur 5I), waren de meeste cellen vroege langzame secretoren met een duidelijke populatie van vroege hoge secretoren en weinig langzame/gemiddelde tot late secretoren. Wat betreft SRlin en Cmax (Figuur 5J) correlaties, waren er bijna geen snelle lage tot gemiddelde secretoren, met alleen een grotere populatie van snelle lage secretoren. Bovendien was er een grote populatie van snelle secretoren die niet afhankelijk waren van de maximaal gemeten concentratie, en twee populaties van hoge secretoren scheidden langzaam of snel uit. Samenvattend kan worden geconcludeerd dat het onderzoeken van de relatie tussen de curvedescriptoren voor individuele cellen een veel gedetailleerdere analyse oplevert en mogelijk nieuwe biologische bevindingen kan extraheren uit metingen van de secretie van eencellige cellen.
Met de hierboven geïntroduceerde analyse hebben we de secretiedynamiek van co-secreterende cellen geëxtraheerd (Figuur 6). Twee voorbeeldcurves tonen een verschillende dynamiek van co-secretie voor IL-6 en TNFα van twee afzonderlijke cellen met een gelijktijdige start van beide cytokines (Figuur 6A), of een sequentiële secretiestart, waarbij IL-6 eerst wordt uitgescheiden (Figuur 6B). Om alle co-secreterende cellen te classificeren, werd een secretievertraging van 60 minuten gedefinieerd, waarbij alle cellen die binnen dit bereik beginnen met secretie worden beschouwd als gelijktijdige secretoren en alle cellen met langere vertragingen als sequentiële secretoren worden beschouwd. Deze analyse maakte het ook mogelijk om te observeren welk cytokine als eerste werd uitgescheiden. Voor IL-6 en TNFα werd voornamelijk gelijktijdige co-secretie waargenomen in 76% van de cellen (Figuur 6C), terwijl voor IL-6 en IL-1β sequentiële co-secretie werd waargenomen in 86% van de cellen, waarbij IL-6 in de meeste gevallen het eerste cytokine was dat werd uitgescheiden (Figuur 6D).
Kijkend naar de starttijd van secretie voor de verschillende cytokinen voor alle individuele co-secreterende cellen, werd er geen duidelijke correlatie waargenomen tussen de starttijden van de secretie in de uitgevoerde experimenten. Voor de co-secretie van IL-6 en TNFα (figuur 6E) was een groter verticaal cluster aanwezig rond 0 minuten, wat overeenkomt met de co-secreterende cellen die vaker beginnen met IL-6. Voor IL-6 en IL-1β co-secretie (figuur 6F) begonnen de meeste cellen IL-6 uit te scheiden rond het begin van de meting, terwijl IL-1β voornamelijk later werd uitgescheiden. Samenvattend maakte de hier gepresenteerde analyse de identificatie mogelijk van verschillende secretorsubpopulaties en complexe cytokine-co-secretiedynamiek.

Figuur 5: Gedetailleerde analyse van verschillende secretiedynamische patronen voor enkelvoudige IL-6 secreterende celcurves. (A) Representatieve gegevens over de concentratie van eencellige cytokines gedurende de meettijd met de passende sigmoïde curve en de geëxtraheerde parameters. (B-D) Drie voorbeeldige eencellige cytokineconcentratiecurven voor de verschillende soorten cytokinesecretoren die worden gevonden voor IL-6-secretie na LPS-stimulatie. (E-G) Percentages IL-6-secreterende cellen die zijn ingedeeld in de verschillende secretortypen met de volgende criteria (n=633): E. Cmax: laag <5 nM, hoog >19,5 nM, F. tstart: vroeg <30 min, laat >120 min, G. SRlin: langzaam <250 moleculen/s, snel >750 moleculen/s. (H-J) Relatie tussen de drie secretiecurve-descriptoren Cmax, tstart en SRlin voor elke individuele cel (n=633). De grote populatie bij Cmax=20nM is het resultaat van het bereiken van de bovenste detectielimiet van de test. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: Extractie van co-secretiepatronen uit enkelcellige concentratiecurves. (A-B) Representatieve concentratiecurven voor enkelvoudige cellen die respectievelijk IL-6 en TNFα (A) gelijktijdig en (B) achtereenvolgens afscheiden. (C-D) Percentage cellen dat gelijktijdige en sequentiële co-secretie van IL-6 en TNFα (n=249), of respectievelijk IL-6 en IL-1β (n=72) vertoont. Sequentiële secretie wordt gedefinieerd door de vertraging tussen de start van de cytokinesecretie van meer dan 60 minuten. Kleuren geven aan welke van de cytokinen het eerst met de secretie is begonnen. (E-F) Relatie tussen de starttijden van de secretie voor de verschillende cytokinen voor elke secreterende cel (nIL6-TNFα=249, nIL6-IL1β=72). Klik hier om een grotere versie van deze figuur te bekijken.