Studiepopulatie en ethiek
Dit was een prospectieve observationele monocentrische studie die werd uitgevoerd in een tertiair opleidingsziekenhuis tussen januari 2026 en april 2026. Het protocol werd goedgekeurd onder goedkeuringsnummer 2026HL-038, en van elke deelnemer werd schriftelijke geïnformeerde toestemming verkregen voorafgaand aan elke studieprocedure.
Volwassenen in de leeftijd van 18 tot 75 jaar die een diagnostische flexibele bronchoscopie met bronchoalveolaire lavage (BAL) ondergingen als onderdeel van de routinematige klinische zorg, werden opeenvolgend gerekruteerd en toegewezen aan een infectiegroep of een niet-infectie controlegroep op basis van hun definitieve klinische en microbiologische diagnoses. Een overzicht van de geïntegreerde visualisatieworkflow, van patiëntenwerving tot de geïntegreerde visuele output, wordt getoond in Figuur 1.

Figuur 1: Overzicht van de geïntegreerde visualisatieworkflow voor BALF-gebaseerde profilering van neutrofielwerving en inflammatoire cytokinerespons. De vijf fasen vatten van links naar rechts samen: prospectieve werving van de infectie- en controlegroepen met bronchoscopie en BAL, BALF-preprocessen tot een celpellet en een supernatantfractie, parallelle cellulaire (Wright-Giemsa-cytologie en multiparameter flowcytometrie) en oplosbare (13-plex bead-gebaseerde immunoassay) metingen, de uniforme visualisatiepijplijn (UMAP met FlowSOM-meta-clustering, hiërarchische heatmap, principale componentanalyse en Spearman-correlatienetwerk), en de geïntegreerde visuele output. Afkortingen: BAL = bronchoalveolaire lavage; BALF = BAL-vloeistof; UMAP = uniform manifold approximation and projection. Klik hier om een grotere versie van deze figuur te bekijken.
De inclusiecriteria voor de infectiegroep waren als volgt: (1) patiënten gediagnosticeerd met respectievelijk community-acquired pneumonia (CAP), hospital-acquired pneumonia (HAP) of ventilator-associated pneumonia (VAP) volgens de richtlijnen van de Infectious Diseases Society of America (IDSA) en de American Thoracic Society (ATS); (2) detectie van één of meer pathogenen via BALF-kweek ≥104 colony-forming units (CFU)/mL, respiratoire multiplex polymerasekettingreactie (PCR), metagenomische next-generation sequencing van BALF, of specifieke antigeen-/antilichaamtesten van gepaarde respiratoire/serummonsters.
De criteria voor de controlegroep waren als volgt: (1) niet-infectieuze redenen voor bronchoscopie, bestaande uit de beoordeling van een onbepaald longnodulus, abnormale radiologische resultaten die niet wezen op een infectie, of onderzoek naar chronische hoest; (2) geen acute pulmonale infectie, aangetoond door negatieve resultaten van microbiologisch onderzoek en de afwezigheid van tekenen van ontsteking of systemische infectie op het moment van de bronchoscopie; en (3) geen gebruik van antibiotica in de 30 dagen voorafgaand aan deelname. De controles werden geworven uit patiënten die een bronchoscopie ondergingen vanwege andere klinisch geïndiceerde pulmonale aandoeningen in plaats van uit gezonde vrijwilligers, om zo de populatie te weerspiegelen waaruit klinisch geïndiceerde BAL routinematig wordt verkregen in het studiecentrum.
De exclusiecriteria waren: (1) maligniteit en/of cytotoxische chemotherapie in de voorafgaande zes maanden; (2) HIV-seropositiviteit; (3) transplantatie van een solide orgaan of hematopoëtische stamcellen; (4) aanhoudende immunosuppressieve therapie met een dagelijkse equivalente dosis prednisolon > 0,3 mg/kg gedurende meer dan 3 weken; (5) zwangere/lacterende vrouwen; (6) voorgeschiedenis van bloedingsdiathese of enige contra-indicaties voor flexibele bronchoscopie; (7) reeds bestaande ernstige chronische respiratoire ziekte (chronisch obstructieve longziekte [COPD] Global Initiative for Chronic Obstructive Lung Disease [GOLD] stadium III/IV, idiopathische pulmonale fibrose, ernstig persisterend astma of klinisch significante bronchiectasieën); (8) actieve tuberculose; (9) deelname aan een andere interventiestudie binnen de laatste 30 dagen. BAL werd uitgevoerd tijdens de eerste bronchoscopie, binnen 72 u na opname bij infecties en binnen 48 u na intubatie bij mechanisch beademde patiënten met infecties.
Steekproefgrootteberekeningen waren gericht op de primaire cellulaire uitkomst: het aantal neutrofielen in BALF tussen de groepen. Gezien de verwachte effectgrootte van 0,6 SD tussen geïnfecteerde en niet-geïnfecteerde luchtwegen, en rekening houdend met een uitvalpercentage van 15%, was een minimale steekproefgrootte van 45 proefpersonen per groep nodig voor een tweezijdige α = 0,05 met een statistische power van 80%. Van de 156 geschikte kandidaten voldeden er 27 niet aan de geschiktheidseisen (18 voldeden niet aan de inclusiecriteria, voornamelijk door het onvermogen om binnen de gestelde periode een pathogeen agens te identificeren dat verantwoordelijk was voor hun ziekte of door een gebrek aan instemming met de procedure zelf; 9 overtraden de exclusiecriteria, hoofdzakelijk vanwege een actieve maligniteit of een voorgeschiedenis van chemotherapie, chronisch gebruik van corticosteroïden, ernstige pre-existent chronische longziekte of een bloedingsdiathese). De resterende 129 deelnemers ondergingen een bronchoscopie. Vervolgens werden 17 proefpersonen uit de studie verwijderd (7 vanwege een onvoldoende BALF-opbrengst < 40%; 4 vanwege het intrekken van de toestemming tijdens opname in het ziekenhuis; en 6 vanwege herclassificatie van proefpersonen na de definitieve microbiologische evaluatie – 4 controlekandidaten die op basis van cultuurresultaten werden gediagnosticeerd met een occulte infectie; en 2 infectiekandidaten bij wie na cultuuranalyse geen bewijs van infectie was). De uiteindelijke groep studieproefpersonen bestond uit 112 deelnemers: 64 met infectie en 48 zonder.
BALF-verzameling en voorbewerking
Bronchoscopie werd uitgevoerd door ervaren operators volgens het standaard institutionele protocol, in overeenstemming met gepubliceerde technische aanbevelingen voor bronchoalveolaire lavage12. Lidocaïne werd topisch gebruikt en patiënten ontvingen matige intraveneuze sedatie, terwijl supplementaire zuurstof werd toegediend om de zuurstofverzadiging (SpO₂) >92% te houden. De flexibele bronchoscoop werd ingebracht en gepositioneerd in het bronchopulmonale segment met radiologische afwijkingen. In gevallen waarin de radiologische afwijkingen diffuus of niet lokaliseerbaar waren, werden de rechter middenkwab of de linker lingulaire segmenten gekozen. Na elke instillatie werd een zachte handmatige aspiratie uitgevoerd met een negatieve druk < 100 mmHg. In totaal werden drie tot vijf instillaties van steriele zoutoplossing bij 37 °C, telkens 20–40 mL, voor een totaalvolume van 100–150 mL, toegediend. De herstelratio werd geschat als het teruggewonnen volume ten opzichte van het geïnstilleerde volume, en alleen monsters met een herstel van ≥40% werden meegenomen voor analyse.
De verzamelde monsters van de tweede en volgende spoelingen werden samengevoegd om bronchiale contaminatie te verminderen en een alveolair-verrijkt monster te verkrijgen. De pool van bronchoalveolaire lavagevloeistof (BALF) werd op ijs bewaard en binnen 2 u na collectie verwerkt, rekening houdend met recente bevindingen over veranderingen in de cellulaire samenstelling en differentiële tellingen na langdurige bewaring en opwarming van menselijke BALF-monsters13. Na filtratie met een steriel nylon gaas met een poriegrootte van 70 microns om het slijm te verwijderen, uitgevoerd door enkel mechanische passage zonder toevoeging van een mucolyticum of antiklontmiddel zoals dithiothreitol, aangezien pilotonderzoek uitwees dat dit voldoende was om slijmklonters te verwijderen zonder meetbaar verlies van de integriteit van oppervlakte-epitopen, werd de vloeistof gecentrifugeerd bij 400 × g gedurende 10 min bij 4 °C. Het supernatant werd afgeschonken en verdeeld in cryovials met een lage eiwitbinding voor cytokinenanalyse, en de cellen werden gesuspendeerd in fosfaatgebufferde zoutoplossing (PBS) met 0,5% runder-serumalbumine en 2 mM ethyleendiaminetetra-azijnzuur (EDTA), waarbij een deel werd bewaard voor cytologie en flowcytometrie. Cytokine-monsters werden ingevroren in vloeibare stikstof, bewaard bij –80 °C en vóór gebruik onderworpen aan maximaal één vries-dooi-cyclus.
Neutrofielcytologie en flowcytometrische kwantificering
Wat betreft de cytopathologie werd 100 µL van de BALF-cel-suspensie met een cytocentrifuge (zie Materialentabel) op een glazen objectglas aangebracht bij 300 × g gedurende 5 min, aan de lucht gedroogd en gekleurd met Wright-Giemsa. Differentiële tellingen werden uitgevoerd door een hematologie-technoloog zonder kennis van de medische achtergrond van de patiënt, die minimaal 400 cellen per objectglas telde en de percentages macrofagen, lymfocyten, neutrofielen en eosinofielen berekende. BALF-neutrofielenpercentages werden beschouwd als de primaire uitkomstmaat voor de cytopathologie, gezien het huidige bewijs dat suggereert dat deze nuttig zijn bij het onderscheid tussen infectieuze en niet-infectieuze longziekten bij patiënten met pulmonale infiltraten14.
Het protocol voor flowcytometrische kwantificering was gebaseerd op een recente validatie van het protocol voor vierkleurige leukocytendifferentiatie in bronchoalveolair spoelvocht (BALF)15 en bevatte aanvullende markers om neutrofielen te identificeren. Specifiek werden ongeveer 5 × 105 cellen geïsoleerd uit BALF gedurende 10 min bij 4 °C geblokkeerd met een human Fc receptor-blocking reagent (zie Tabel met materialen>) en vervolgens gedurende 20 min bij 4 °C gekleurd met antilichamen tegen CD45 (kloon HI30), CD66b (kloon G10F5), CD15 (kloon HI98), CD16 (kloon 3G8), CD11b (kloon ICRF44), CD14 (kloon M5E2) en HLA-DR (kloon L243), naast een fixeerbare viabiliteitskleurstof (zie Tabel met materialen>). Voor erytrocytlyse werd een commerciële erytrocyt-lyseoplossing (zie Tabel met materialen>) gedurende 10 min bij kamertemperatuur gebruikt. Na erytrocytlyse werden de cellen twee keer gewassen met PBS aangevuld met 2% foetaal kalfserum en geresuspendeerd in flowbuffer, waarna ze werden geanalyseerd met een multicolor flowcytometer (zie Tabel met materialen>). Compensatie werd uitgevoerd met gebruik van enkelvoudig gekleurde compensatiekraaltjes (zie Tabel met materialen>), en fluorescence-minus-one en ongekleurde buisjes werden voor elk paneel meegenomen.
Voor elk monster werden ten minste 1 × 105 CD45⁺ events verzameld. Eerst werden singlets op de FSC-A/FSC-H plot gegated, gevolgd door levende (viability dye⁻) CD45⁺ leukocyten, waarna uitsluiting van CD14⁺HLA-DR⁺ monocyten-macrofagen werd uitgevoerd om high-side scatter CD66b⁺ granulocyten te selecteren. Neutrofielen werden vervolgens binnen de granulocytenpopulatie gedefinieerd als CD15⁺CD16⁺CD11b⁺ events. Deze primaire CD66b⁺ granulocyten-gate werd breed getrokken rond alle high-side-scatter events in plaats van te worden beperkt tot de CD15⁺CD16⁺CD11b⁺ grens, zodat CD16-dim, CD11b-high geactiveerde neutrofielen en andere minder volwassen granulocytenvormen binnen de totale geautomatiseerde neutrofielentelling bleven; de striktere CD15⁺CD16⁺CD11b⁺ definitie werd alleen toegepast om de subset volwassen neutrofielen te karakteriseren die afzonderlijk werd gerapporteerd in de daaropvolgende clusteringanalyses. De neutrofielentelling in BALF (cellen/µL) werd geschat op basis van het percentage neutrofielen en het totale aantal kernhoudende cellen bepaald met een geautomatiseerde hematologie-analyzer (zie Table of Materials). Compensatiematrices, de gating-strategie en individuele kwaliteitscontolerapporten per monster werden gearchiveerd, en analisten bleven gedurende de acquisitie en analyse blind voor de toewijzing van de deelnemersgroepen.
Multiplex profiling van inflammatoire cytokinen
BALF-supernatanten werden één keer ontdooid op nat ijs, gevolgd door kort vortexen en klaring door centrifugatie bij 10.000 × g gedurende 5 min bij 4 °C direct voorafgaand aan de assay. De meting van 13 analyten in drie functionele categorieën werd uitgevoerd met een multiplex bead-gebaseerd immunoassay-platform met een aangepast magnetisch bead-panel (zie Tabel met materialen). De pro-inflammatoire groep bestond uit IL-1β, IL-6, IL-17A en TNF-α; de chemokinegroep bestond uit CXCL8/IL-8, CXCL1/GRO-α, CXCL10/IP-10, CCL2/MCP-1 en CCL3/MIP-1α; en de regulatoire groep bestond uit IL-10, IFN-γ, IL-1RA en G-CSF. De keuze voor de analyten was gebaseerd op diegenen die het vaakst als informatief werden gerapporteerd in multiplex-analyses van BALF- of serumcytokinen in studies naar pneumonie/ernstige respiratoire infecties7,11. De verdunningsfactor (1:2 in assaybuffer) werd bepaald op basis van een pilotexperiment uitgevoerd op een representatieve deelset van zowel geïnfecteerde als controlemonsters om ervoor te zorgen dat >90% van de metingen binnen de standaardcurve viel.
De standaardcurves werden bepaald met behulp van zevenpunts seriële verdunningen van de recombinante standaarden van de fabrikant, die in dubblo replicates op elke plaat werden klaargemaakt, inclusief één blanco put om de achtergrond af te trekken. De curvefitting werd uitgevoerd via een vijfparameter logistieke (5PL) regressie door de eigen acquisitie- en analysesoftware van het instrument (zie Tabel met materialenDe onderste kwantificeringsgrens (LLOQ) werd berekend als de kleinste hoeveelheid standaard die een back-fit recovery tussen 80% en 120% en een intra-run variatiecoëfficiënt (CV) van ≤ 20% opleverde. Waarden die onder de LLOQ vielen, werden vervangen door LLOQ/√2. Twee gepoolde kwaliteitscontrolemonsters van bronchoalveolaire lavagevloeistof (BALF), telkens afkomstig van een representatieve geïnfecteerde patiënt en een niet-aangetast subject, werden in aliquoten verdeeld en aan elke plaat toegevoegd om de prestaties op lange termijn te volgen. De kwaliteitsborgingsnormen waren gebaseerd op best-practice-aanbevelingen afgeleid van een 12-jarig multisite-bekwaamheidsprogramma voor multiplex bead-based assays.16inclusief eisen voor een intra-assay CV ≤ 10% in duplex wells, een inter-plate CV ≤ 20% voor kwaliteitscontrole en een minimaal aantal beads van ≥ 50 per analyt per well. Platen die niet aan de bovenstaande criteria voldeden, werden opnieuw uitgevoerd. De analyse werd in één batch binnen vier weken uitgevoerd, met gebruik van één kit-lot, waarbij de infectie- en controlesamples willekeurig werden afgewisseld.
Pipeline voor visualisatieanalyse
De Flow Cytometry Standard (FCS)-bestanden werden eerst geanalyseerd met commerciële flowcytometrie-software (zie Tabel met Materialen) om levende CD45⁺ singlets te compenseren en te gaten, en werden vervolgens geëxporteerd voor onbevooroordeelde analyse in R v4.4. Voor elk monster werden 10.000 CD45⁺ events gedownsampled en geconcateneerd, en de signaalintensiteiten werden getransformeerd met de arcsinh-functie met een cofactor van 150. Dimensiereductie en UMAP-projecties werden uitgevoerd met het uwot-pakket bij 15 dichtstbijzijnde buren en een minimale afstand van 0,1, samen met parallelle t-SNE-embeddings met het Rtsne-pakket met een perplexiteit van 3010. Fenotypeclusters werden geïdentificeerd met FlowSOM met 100 self-organizing map-nodes en verder onderverdeeld in 12–15 clusters via meta-clustering, waarbij gebruik werd gemaakt van een optimalisatiegebaseerde benadering die onlangs is gevalideerd voor hoogdimensionale cytometriegegevens17. Elk cluster werd geannoteerd op basis van de mediane expressie van lineage- en activatiemarkers en gecontroleerd aan de hand van de handmatige gates.
Voor de analyse van multiplex cytokine-gegevens werden de ruwe waarden eerst log₁₀-getransformeerd en vervolgens gestandaardiseerd tussen analyten via z-scoring voorafgaand aan aanvullende analyses. Zowel analyten als patiënten werden onderworpen aan hiërarchische clustering met behulp van Euclidische afstand en Ward.D2-koppeling; de resultaten werden weergegeven in een heatmap met clustering van rijen en kolommen. Een dergelijke benadering werd gebruikt om co-expressie-modules van cytokines te detecteren en patiëntsubsets in multiplex inflammatoire cohorten te classificeren9. Correlatiematrices werden gegenereerd voor alle cellulaire en oplosbare biomarkers met behulp van paarsgewijze Spearman-rangcorrelaties. Heatmaps van de correlatiematrices werden uitgezet naast een netwerk waarin knooppunten overeenkomen met biomarkers en randen wijzen op correlaties die de correctie voor meervoudig testen via de Benjamini–Hochberg (BH)-methode hebben doorstaan. De breedte van de randen werd geschaald volgens de grootte van de correlatiecoëfficiënt |ρ|, en randen werden alleen opgenomen als |ρ| ≥ 0,30 en p < 0,05 na correctie. Principaalcomponentanalyse (PCA) werd uitgevoerd op de gestandaardiseerde matrix met cellen en cytokines met behulp van de prcomp()-functie met unit variance scaling. De deelnemers werden geprojecteerd op de eerste twee hoofdcomponenten, waarbij infectie- en controlegroepen werden overgelegd. Consistente downsampling en standaardisatie van alle analyses waarborgden de vergelijkbaarheid tussen cellulaire en oplosbare biomarker-uitlezingen in alle visualisaties.
Statistische analyse
Continue variabelen werden beschreven met behulp van de mediaan (interkwartielafstand, IQR) bij een niet-normale verdeling volgens de Shapiro–Wilk-toetsen, en anders met het gemiddelde (standaarddeviatie, SD). Categorische variabelen werden beschreven aan de hand van frequenties en proporties. Vergelijkingen tussen groepen werden uitgevoerd met de Mann–Whitney U-toets wanneer continue variabelen een niet-normale verdeling hadden, de Welch t-toets wanneer continue variabelen normaal verdeeld waren, en de chi-kwadraattoets/Fisher's exact test wanneer dit van toepassing was voor categorische variabelen. De associaties tussen continue variabelen werden beoordeeld via de rangcorrelatiecoëfficiënt van Spearman. Om rekening te houden met de false discovery rate (FDR) bij multiplex cytokine- en correlatieanalyses, werden p-waarden gecorrigeerd met de Benjamini–Hochberg-methode met een afkapwaarde van q < 0,05. Tweezijdige p < 0,05 werd beschouwd als statistisch significant in alle primaire analyses. Waar mogelijk werden gevallen met ontbrekende gegevens over individuele variabelen behouden en geanalyseerd via complete-case-analyse voor de betreffende variabele, terwijl meerdere imputaties met behulp van chained equations (m = 20) werden uitgevoerd om sensitiviteitsanalyses door te voeren. Receiver operating characteristic-curven en de bijbehorende area under the curve (AUC) werden gegenereerd voor individuele cytokines met behulp van het pROC-pakket binnen de volledige onderzoekscohort; er werd geen aparte externe validatiecohort gebruikt. Vooraf gespecificeerde sensitiviteitsanalyses herhaalden de primaire groepvergelijkingen na uitsluiting van deelnemers die binnen 72 h vóór de bronchoscopie antibiotica hadden ontvangen, en na stratificatie van de infectiegroep naar status van mechanische ventilatie en naar pneumoniesubtype (community-acquired, hospital-acquired of ventilator-associated). Analyses werden uitgevoerd met de pakketten rstatix, FlowSOM, uwot, ComplexHeatmap, igraph, pROC en mice in R v4.4. Figuren werden gegenereerd met commerciële wetenschappelijke grafieksoftware (zie Table of Materials).