Studiepopulatie en monsterkenmerken
Er waren 112 proefpersonen in de uiteindelijke analysecohort, met 64 proefpersonen in de infectiegroep en 48 proefpersonen in de niet-infectie controlegroep. De baseline demografie was vergelijkbaar tussen de groepen, terwijl laboratoriumparameters gerelateerd aan infectie hoger waren in de infectiegroep, zoals weergegeven in Tabel 1. Infecties werden primair veroorzaakt door Gram-negatieve bacteriën, waarbij de belangrijkste bacteriële pathogenen Klebsiella pneumoniae, Pseudomonas aeruginosa en Acinetobacter baumannii waren. Andere pathogenen waren Streptococcus pneumoniae, Staphylococcus aureus, atypische agentia, virussen en gemengde infecties.
De herstelsnelheden van de BALF voldeden aan de vooraf vastgestelde drempelwaarde van ≥40% voor alle behouden monsters, met iets lagere percentages in de infectiegroep vergeleken met de controlegroep (mediaan 48,3% vs 53,3%, p = 0,03), wat consistent is met verhoogde ontsteking van de luchtwegen en mucosale congestie in de geïnfecteerde luchtwegen. Omdat het herstelde volume ook verschilde tussen de groepen (58 mL vs 64 mL, p = 0,04), kunnen alle gerapporteerde cytokoneconcentraties onderhevig zijn aan een bescheiden verschil in verdunningsfactor tussen de groepen; een vooraf vastgestelde sensitiviteitsanalyse beperkt tot monsters met een herstelsnelheid van ≥50% (n = 42 infectie, n = 34 controle) reproduceerde de richting en omvang van elke vergelijking tussen de groepen (grootste verandering in effectgrootte < 8%), wat erop wijst dat verdunningsbias waarschijnlijk niet verantwoordelijk is voor de waargenomen verschillen. De tijd tussen collectie en verwerking was 65 min voor de infectiegroep en 62 min voor de controlegroep, wat binnen het venster van 2 u bleef. Kwaliteitscontrolecriteria (QC) op monsterniveau duidden op resultaten van hoge kwaliteit. Cytologische differentiële tellingen werden verkregen op basis van ≥400 getelde cellen in alle 112/112 monsters (100%). De flowcytometrie-acquisities voldeden aan het vooraf vastgestelde criterium van ≥1 × 105 CD45⁺ events in alle 112/112 monsters (100%), met een mediane levensvatbaarheid van levende cellen van 89,4% (IQR 85,6–92,8%). Alle multiplex cytokoneplaten voldeden aan de QC-criteria na één enkele herhaling. Voor alle 13 analyten was de intra-assay variatiecoëfficiënt (CV) 4,7% (2,1–9,3%, alle <10%), de inter-plaat CV van gepoolde QC-monsters was 11,2% (6,8–18,4%, alle <20%), en het mediane aantal beads was 78 events/analyte/well (53–142, alle ≥ 50).
| Variabele | Infectie (n = 64) | Controle (n = 48) | p-waarde |
| Leeftijd, jaren, mediaan [IQR] | 62 [54–71] | 60 [52–68] | 0.41 |
| Man, n (%) | 38 (59.4) | 26 (54.2) | 0.59 |
| BMI, kg/m2, mediaan [IQR] | 23.8 [21.4–26.2] | 23.5 [21.6–25.7] | 0.66 |
| Rookgeschiedenis, n (%) | | | 0.52 |
| Nooit | 28 (43.8) | 25 (52.1) | |
| Voorheen | 22 (34.4) | 16 (33.3) | |
| Huidig | 14 (21.9) | 7 (14.6) | |
| Hypertensie, n (%) | 27 (42.2) | 18 (37.5) | 0.62 |
| Diabetes mellitus, n (%) | 17 (26.6) | 10 (20.8) | 0.49 |
| Coronaire hartziekte, n (%) | 11 (17.2) | 6 (12.5) | 0.51 |
| Chronische nierziekte, n (%) | 5 (7.8) | 3 (6.3) | 0.76 |
| Indicatie voor bronchoscopie, n (%) | | | — |
| Community-acquired pneumonie | 35 (54.7) | — | |
| Hospitaal-acquired pneumonie | 19 (29.7) | — | |
| Ventilator-associated pneumonie | 10 (15.6) | — | |
| Onbepaalde longnodule | — | 23 (47.9) | |
| Aanhoudende radiografische abnormaliteit | — | 15 (31.3) | |
| Evaluatie van chronische hoest | — | 10 (20.8) | |
| Veroorzaker pathogeen, n (%) | | — | |
| Gram-negatieve bacteriën | 28 (43.8) | | |
| Gram-positieve bacteriën | 16 (25.0) | | |
| Atypische pathogenen | 7 (10.9) | | |
| Respiratoire virussen | 8 (12.5) | | |
| Gemengde infectie | 5 (7.8) | | |
| Mechanische ventilatie bij BAL, n (%) | 14 (21.9) | 0 (0) | <0.001 |
| Blootstelling aan antibiotica binnen 72 u vóór BAL, n (%) | 41 (64.1) | 0 (0) | <0.001 |
| Instillatievolume, mL, mediaan [IQR] | 120 [120–140] | 120 [120–140] | 0.84 |
| Herwonnen volume, mL, mediaan [IQR] | 58 [49–68] | 64 [55–73] | 0.04 |
| Herstelpercentage BALF, %, mediaan [IQR] | 48.3 [42.5–55.6] | 53.3 [46.2–60.5] | 0.03 |
| Tijd tot verwerking, min, mediaan [IQR] | 65 [50–85] | 62 [48–80] | 0.55 |
| Witte bloedcellen in bloed, × 10⁹/L, mediaan [IQR] | 11.8 [8.6–14.5] | 6.4 [5.3–7.8] | <0.001 |
| Neutrofielen in bloed, %, mediaan [IQR] | 78.5 [70.2–86.4] | 60.3 [54.6–66.8] | <0.001 |
| C-reactief proteïne, mg/L, mediaan [IQR] | 78.5 [42.6–128.4] | 4.2 [1.8–8.6] | <0.001 |
| Procalcitonine, ng/mL, mediaan [IQR] | 0.86 [0.32–2.45] | 0.05 [0.03–0.08] | <0.001 |
Tabel 1: Demografische, klinische en BALF-verwerkingskenmerken van de onderzoekscohort. De baseline demografische kenmerken, comorbiditeiten, indicaties voor bronchoscopie, veroorzakende pathogenen, procedurele parameters en laboratoriumwaarden van perifeer bloed worden vergeleken tussen de 64 deelnemers in de infectiegroep en de 48 deelnemers in de niet-infectie controlegroep. Continue variabelen worden gerapporteerd als mediaan [interkwartielafstand] en vergeleken met de Mann-Whitney test. U test; categorische variabelen worden gerapporteerd als n (%) en vergeleken met de chi-kwadraattoets of de exacte toets van Fisher, indien van toepassing. Afkortingen: BAL = bronchoalveolaire lavage; BALF = bronchoalveolaire lavagevloeistof; BMI = lichaamsmassaindex; IQR = interkwartielafstand; WBC = aantal witte bloedcellen.
Visualisatie van de rekrutering van neutrofielen in BALF
Een duidelijke verschuiving in de ratio van een door macrofagen gedomineerde populatie naar een door neutrofielen gedomineerde populatie werd waargenomen bij de geïnfecteerde proefpersonen, zoals getoond in Figuur 2. De cytocentrifuge-preparaten werden gekleurd met de Wright-Giemsa-techniek, en de kleuring toonde een duidelijk onderscheid tussen de twee groepen, waarbij de controlegroep meer alveolaire macrofagen, enkele lymfocyten en een paar neutrofielen bevatte, terwijl de geïnfecteerde groep gedomineerd werd door neutrofielen met frequente intracellulaire bacteriën, gedegenereerde kernen en toxische granulaties. Differentiële tellingen op basis van ten minste 400 getelde cellen werden gebruikt om de verschillen te verifiëren, waarbij de neutrofielenratio steeg van een mediaanwaarde van 3,2% (IQR: 1,8–5,6%) in de controlegroep naar 68,4% (IQR: 49,2–82,5%) in de geïnfecteerde groep (p < 0,001), terwijl die van de alveolaire macrofagen afnam van 88,2% (IQR: 82,4–92,5%) naar 22,6% (IQR: 12,3–38,5%) (p < 0,001). De percentages lymfocyten waren vergelijkbaar tussen de groepen (7,4% [4,8–11,3%] vs 6,5% [3,4–11,2%], p = 0,42).
De gating- en dimensionaliteitsreductieanalyses bevestigden en breidden de cellulaire fenotypes uit die via microscopie waren waargenomen. Na het uitsluiten van debris, doubletten en dode cellen, vormden CD45⁺CD66b⁺ granulocyten een mediaan van 65,8% (IQR: 46,2–79,4%) van de CD45⁺ leukocyten in de infectiegroep tegenover 4,1% (IQR: 2,3–6,8%) in de controlegroep (p < 0,001). Binnen de granulocytenpopulatie vormden CD15⁺CD16⁺CD11b⁺ mature neutrofielen 92,4% (IQR: 87,5–95,8%) van alle granulocyten in de infectiegroep tegenover 88,6% (IQR: 82,3–93,5%) in de controlegroep, wat aantoont dat de gerekruteerde pool neutrofielen grotendeels bestond uit mature, volledig gedifferentieerde cellen in plaats van vroege precursoren. De UMAP- en t-SNE-projecties op basis van de arcsinh-getransformeerde oppervlakte-eiwitmatrix vertoonden vrijwel identieke topologieën waarin infectiegevallen een prominent eiland van CD15⁺CD16⁺CD11b⁺ mature neutrofielen bezetten dat in controlemonsters schaars bevolkt was, terwijl het eiland van HLA-DR⁺CD14⁺ alveolaire macrofagen het tegenovergestelde patroon vertoonde: dicht bevolkt in controles en schaars bevolkt bij infectie. Meta-clustering via FlowSOM identificeerde 13 subpopulaties, waaronder een subpopulatie mature neutrofielen (38,4% [22,6–56,2%] vs 2,8% [1,5–4,6%], p < 0,001) en een geactiveerde CD16-dim CD11b-high neutrofielencluster (22,6% [11,4–32,5%] vs 0,4% [0,1–1,2%], p < 0,001) die de grootste expansies tussen de groepen vertoonden.
Wat betreft het absolute aantal neutrofielen was er een ongeveer 50-voudige toename van neutrofielen in de BALF in de infectiegroep vergeleken met de controlegroep. Het totaal aantal gekerneerde cellen per µL BALF steeg van 156 cellen/µL (IQR 88–232) in de controlegroep naar 384 cellen/µL (IQR 180–712) in de infectiegroep (p < 0,001), en het aantal neutrofielen steeg van 5 cellen/µL (IQR 2–13) in de controlegroep naar 262 cellen/µL (IQR 89–578) in de infectiegroep (p < 0,001). Paarsgewijze verschillen in het percentage en het absolute aantal neutrofielen tussen de groepen bleven zeer significant na Benjamini–Hochberg (BH) correctie (q < 0,001).

Figuur 2: Visualisatie van de rekrutering van neutrofielen in BALF. Representatieve Wright-Giemsa-gekleurde cytospin-velden van een controledeelnehmer en een geïnfecteerde deelnemer illustreren de verschuiving van een door macrofagen gedomineerd naar een door neutrofielen gedomineerd cellulair BALF-profiel. Een gating-strategie voor flowcytometrie in vier stappen (FSC singlets → levende CD45⁺ leukocyten → CD66b⁺ granulocyten na CD14⁺HLA-DR⁺ exclusie → CD15⁺CD16⁺CD11b⁺ mature neutrofielen), een UMAP-embedding van arcsinh-getransformeerde flowcytometrie-gegevens met een FlowSOM 13-metacluster overlay gekleurd op basis van biologische annotatie, en boxplots van het BALF-neutrofielenpercentage en de absolute neutrofielenconcentratie (cellen/µL, logaritmische schaal) worden zij aan zij getoond om een directe vergelijking mogelijk te maken van de cellulaire fenotypes die zijn geïdentificeerd in de infectie- en controlegroepen. Afkortingen: BALF = bronchoalveolaire lavagevloeistof; FSC = forward scatter; UMAP = uniform manifold approximation and projection. Klik hier om een grotere versie van deze figuur te bekijken.
Visualisatie van ontstekingscytokineprofielen in BALF
Een patroon van gecoördineerde reacties over meerdere cytokine-assen onderscheidde de infectiegroep van de controlegroep, en de cytokines werden ingedeeld in drie consistente functionele modules, zoals weergegeven in Figuur 3. Clustering van log₁₀-getransformeerde en genormaliseerde cytokineconcentraties met behulp van de Ward.D2-methode resulteerde in een duidelijk tweedelig dendrogram van monsters, waarbij infectiemonsters en controles zich aan tegenovergestelde uiteinden van het dendrogram bevonden en slechts zeven verkeerd geclassificeerde deelnemers werden waargenomen (zes controles geclusterd met infectiemonsters, voornamelijk gevallen van chronische hoest met milde neutrofiele luchtweginflammatie, en één infectiemonster geclusterd met controles, wat een vroege atypische pneumonie vertegenwoordigde). Het dendrogram van de analyten definieerde drie cytokineclusters. Het chemokinecluster, bestaande uit CXCL8, CXCL1, CXCL10, CCL2 en CCL3, vertoonde de grootste fold change tussen de infectie- en controlegroepen en de hoogste intramodulaire correlatie. Het pro-inflammatoire cluster, met IL-1β, IL-6, IL-17A en TNF-α, vertoonde vergelijkbare directionele verhogingen, hoewel deze enigszins minder uitgesproken waren. Het regulatoire cluster, inclusief IL-10, IL-1RA, IFN-γ en G-CSF, was eveneens verhoogd in de infectiemonsters, wat aangeeft dat de respons niet louter pro-inflammatoir was, maar parallelle counter-regulatoire en antivirale mechanismen activeerde.
De medianen (IQR's) van alle geteste analyten waren hoger in de infectiegroep vergeleken met de controlegroep, en de aangepaste p-waarde was kleiner dan 0,001 voor alle 13 analyten na BH-gebaseerde correctie voor de false-discovery-rate (BH-FDR). Van de analyten vertoonde CXCL8/IL-8 de hoogste concentratie en de sterkste fold-change (van 68 [32–142] pg/mL in de controlegroep naar 2.485 [842–6.520] pg/mL in de infectiegroep, een toename van ongeveer 36-voudig). Een parallel patroon werd gezien voor CXCL1/GRO-α, met een toename van ongeveer 21-voudig (78 [36–148] vs 1.672 [524–3.890] pg/mL). IL-1β vertoonde een toename van ongeveer 13-voudig (5,8 [2,6–11,4] vs 78 [29–184] pg/mL). IL-6, IL-1RA en G-CSF vertoonden toenames van 13- tot 16-voudig (respectievelijk 16,4 [7,8–30,5] vs 248 [82–612] pg/mL; 112 [48–258] vs 1.748 [658–3.964] pg/mL; 24 [10–52] vs 384 [125–912] pg/mL). Andere kleinere toenames waren TNF-α (~8-voudig), CXCL10/IP-10 (~9-voudig), CCL2/MCP-1 (~8-voudig), CCL3/MIP-1α (~10-voudig), IL-10 (~10-voudig), IFN-γ (~6-voudig) en IL-17A (~4-voudig).
Receiver operating characteristic-analyse binnen de cohort identificeerde IL-8 (AUC 0,94), IL-1β (AUC 0,92), IL-6 (AUC 0,91), CXCL1 (AUC 0,91) en G-CSF (AUC 0,89) als de analyten met de sterkste discriminatie per enkele marker tussen de infectie- en controlegroepen; deze AUC-waarden vatten de discriminatie binnen de analysecohort samen bij afwezigheid van een aparte externe validatieset. Sensitiviteitsanalyse via multiple imputatie door middel van geketende vergelijkingen (m = 20) leverde schattingen op binnen 5% van de medianen van de complete cases voor elke analyt en veranderde de rangorde van de effectgroottes niet.
Vooraf vastgestelde sensitiviteitsanalyses werden gebruikt om de invloed van eerdere antibioticumbehandeling en het type pneumonie op de waargenomen cytokineverschillen te beoordelen. Wanneer de primaire groepvergelijkingen werden herhaald na uitsluiting van de 41 deelnemers uit de infectiegroep (64,1%) die binnen 72 u sebelum de BAL antibiotica hadden ontvangen, bleven alle 13 analyten significant hoger in de resterende 23 infectiegevallen dan in de controles na BH-FDR-correctie (gecorrigeerde p < 0,001 voor de vijf belangrijkste discriminanten), met puntschattingen die iets hoger waren dan de medianen van de gehele cohort (bijvoorbeeld, CXCL8/IL-8 2.860 [980–7.110] pg/mL en IL-1β 92 [34–210] pg/mL), wat consistent is met een gedeeltelijke demping van luchtwegcytokinen door antibioticumbehandeling voorafgaand aan de BAL. Stratificatie van de infectiegroep naar type pneumonie toonde dezelfde rangorde van chemokinen en pro-inflammatoire cytokinen aan voor community-acquired pneumonia (n = 35), hospital-acquired pneumonia (n = 19) en ventilator-associated pneumonia (n = 10), waarbij de subgroep met ventilator-associated pneumonia de hoogste mediane waarden vertoonde voor CXCL8, CXCL1, IL-6 en G-CSF, en de status van mechanische ventilatie de richting van geen enkele tussen-groepvergelijking veranderde.

Figuur 3: Visualisatie van ontstekingscytokineprofielen in BALF. Hiërarchische clustering heatmap van log₁₀-getransformeerde en z-score-gecorrigeerde concentraties van de 13 gemeten analyten bij alle 112 deelnemers, met gekleurde annotatiebalken boven de heatmap die de groepsindeling (controle, infectie), het subtype pneumonie (CAP, HAP of VAP pneumonie) en de categorie van het veroorzakende pathogeen (Gram-negatief, Gram-positief, atypisch, viraal of mixed infection) aangeven. Boxplots van de belangrijkste differentieel abundante cytokinen (CXCL8/IL-8, IL-1β, IL-6, CXCL1/GRO-α, G-CSF en TNF-α; log₁₀ schaal, pg/mL) vergelijken de infectie- en controlegroepen. *** p < 0.001 na BH false-discovery-rate correctie. Afkortingen: BALF = bronchoalveolair spoelvocht; BH = Benjamini–Hochberg; CAP = community-acquired pneumonia; HAP = hospital-acquired pneumonia; VAP = ventilator-associated pneumonia. Klik hier om een grotere versie van deze figuur te bekijken.
Gekoppelde visualisatie van neutrofielrekrutering en cytokiinprofielen
De koppeling van de cellulaire en oplosbare compartimenten onthulde een sterke biologische verbinding tussen de rekrutering van neutrofielen en het chemokine-systeem, zoals weergegeven in Figuur 4. Een tweedimensionale weergave werd afgeleid uit een principale componentanalyse (PCA) van de gestandaardiseerde cellulaire en cytokinematrix, waarbij PC1 verantwoordelijk was voor 51,1% van de totale variantie en PC2 voor een verdere 12,7%, wat neerkomt op een cumulatief percentage van 63,8%. PC1 scheidde de infectiegroep van de controles, gekenmerkt door hoge positieve loadings voor BALF neutrofielen, CXCL8, CXCL1, IL-6, IL-1β en G-CSF, en hoge negatieve loadings voor macrofagen. PC2 maakte onderscheid tussen een chemokine-rijk traject (hoge positieve loadings voor CXCL8 en CXCL1) en een regulatoire-rijk traject (positieve loadings voor IL-10, IL-1RA en IFN-γ), waardoor deelnemers binnen de infectiegroep geneigd waren zich te scheiden op basis van de ernst, waarbij gevallen van beademingsgebonden pneumonie bij mechanisch beademde patiënten in het chemokine-rijke kwadrant vielen.
Paarsgewijze Spearman-rangcorrelaties tussen 17 cellulaire en oplosbare kenmerken onthulden een sterke positieve correlatie tussen het aantal neutrofielen in de BALF en de algemene chemokine-as (CXCL8 ρ = 0,78; CXCL1 ρ = 0,74; CCL3 ρ = 0,69; CCL2 ρ = 0,61; CXCL10 ρ = 0,54; alle gecorrigeerde p < 0,001). Een tweede niveau van significante correlaties werd gedetecteerd tussen het aantal neutrofielen en de pro-inflammatoire as (G-CSF ρ = 0,65; IL-1β ρ = 0,62; IL-6 ρ = 0,58; TNF-α ρ = 0,51), gevolgd door matige positieve correlaties met de regulatoire as (IL-10 ρ = 0,45; IL-1RA ρ = 0,48). Het percentage macrofagen vertoonde een gespiegelde negatieve correlatie van vergelijkbare sterkte met zowel de chemokine- als de pro-inflammatoire as. Correlaties binnen de as waren aanzienlijk sterker dan correlaties tussen assen (mediaan |ρ| binnen de as = 0,68 versus mediaan |ρ| tussen assen = 0,33), wat het concept ondersteunt van gemeenschappelijke stroomopwaartse regulatoire paden die de chemokine-, pro-inflammatoire en regulatoire assen controleren.
Het correlatienetwerk bracht deze koppeling in kaart op een graaf waarin BALF-neutrofielen de meest centrale hub-node vormden, waarbij alle 13 verbindingen naar cytokines de BH-correctie voor meervoudige toetsing doorstonden. Van de 136 mogelijke paarsgewijze verbindingen tussen de 17 kenmerken voldeden er 99 aan het gezamenlijke criterium van |ρ| ≥ 0,30 en BH-gecorrigeerde p < 0,05 en werden behouden. De graaf van het correlatienetwerk vertoonde een duidelijke modulaire architectuur van drie blokken die de clusters uit de hiërarchische clustering reproduceerde: de dichtste en sterkste verbindingen liepen tussen de chemokine-module en het aantal en percentage BALF-neutrofielen, verbindingen van gemiddeld gewicht koppelden de pro-inflammatoire module aan zowel de neutrofielen als de chemokine-module, en dunne verbindingen vanuit de regulatoire module overbrugden de chemokine- en pro-inflammatoire blokken. De node voor het percentage macrofagen bevond zich tegenover het neutrofiel-chemokineblok en was hiermee verbonden via een waaier van negatief gewogen verbindingen. Collectief onthulden PCA, de correlatiematrix en de netwerkgraaf een samenhangend patroon waarin BALF-neutrofieleninfiltratie en de chemokine-respons sterk geassocieerd zijn op het niveau van de individuele deelnemer, en de pro-inflammatoire en regulatoire paden dienen als goed gedefinieerde maar onderling verbonden modules rond deze neutrofiel-chemokine-as.

Figuur 4: Gekoppelde visualisatie van neutrofielrekrutering en ontstekingscytokineprofielen. Biplot van de hoofdcomponentanalyse van de gestandaardiseerde cellulaire en cytokinematrix met deelnemers ingekleurd per groep, Spearman-correlatieheatmap en netwerkgrafiek van alle 17 cellulaire en oplosbare kenmerken (knopen gegroot naar graad, verbindingen geschaald naar |ρ| en gefilterd op BH-gecorrigeerde p < 0,05 en |ρ| ≥ 0,30), waarbij de neutrofiel-chemokine-hub wordt benadrukt. Afkortingen: BH = Benjamini–Hochberg; PC = hoofdcomponent. Klik hier om een grotere versie van deze figuur te bekijken.
DATABESCHIKBAARHEID:
De geanonimiseerde dataset op deelnemerniveau en de R-analysecode die in deze studie zijn gebruikt, zijn beschikbaar als Aanvullende Tabel S1 en Aanvullend Bestand 1. Andere gegevens ter ondersteuning van de resultaten van deze studie zijn op redelijk verzoek verkrijgbaar via de corresponderende auteur.
Aanvullende Tabel S1: BALF-dataset. Klik hier om dit bestand te downloaden.
Aanvullend bestand 1: Een ZIP-bestand met analysecodes. Klik hier om dit bestand te downloaden.