Población de estudio y ética
Este fue un estudio observacional prospectivo monocéntrico realizado en un hospital docente terciario entre enero de 2026 y abril de 2026. El protocolo fue aprobado bajo el número de aprobación 2026HL-038, y se obtuvo el consentimiento informado por escrito de todos los participantes antes de cualquier procedimiento del estudio.
Se reclutaron consecutivamente adultos de 18 a 75 años que se sometieron a broncoscopia flexible diagnóstica con lavado broncoalveolar (LBA) como parte de la atención clínica de rutina, y se asignaron a un grupo de infección o a un grupo control sin infección según sus diagnósticos clínicos y microbiológicos finales. Una visión general del flujo de trabajo de visualización integrada, desde el reclutamiento de pacientes hasta las salidas visuales integradas, se muestra en la Figura 1.

Figura 1: Vista general del flujo de trabajo integrado de visualización para el perfilado basado en LBA de la reclutación de neutrófilos y las respuestas inflamatorias de citocinas. Las cinco etapas resumen, de izquierda a derecha, el reclutamiento prospectivo de los grupos de infección y control mediante broncoscopia y LBA, el preprocesamiento del LBA en un sedimento celular y una fracción de sobrenadante, mediciones celulares paralelas (citología con Wright-Giemsa y citometría de flujo multiparamétrica) y solubles (inmunoensayo basado en perlas de 13 plex), la tubería unificada de visualización (UMAP con metaagrupamiento FlowSOM, mapa de calor jerárquico, análisis de componentes principales y red de correlación de Spearman) y la salida visual integrada. Abreviaturas: BAL = lavado broncoalveolar; LBA = líquido de lavado broncoalveolar; UMAP = proyección y aproximación uniforme de variedades. Haga clic aquí para ver una versión más grande de esta figura.
Los criterios de inclusión para el grupo de infección fueron los siguientes: (1) pacientes diagnosticados con neumonía adquirida en la comunidad (NAC), neumonía adquirida en el hospital (NAH) o neumonía asociada al ventilador (NAV), según las guías de la Sociedad Americana de Enfermedades Infecciosas (IDSA) y la Sociedad Torácica Americana (ATS); (2) detección de uno o más patógenos mediante cultivo del lavado broncoalveolar (LBA) con ≥104 unidades formadoras de colonias (UFC)/mL, reacción en cadena de la polimerasa (PCR) respiratoria múltiple, secuenciación genómica de nueva generación del LBA o pruebas específicas de antígeno/anticuerpo a partir de muestras respiratorias y sueros pareadas.
Los criterios para el grupo de control fueron los siguientes: (1) razones no infecciosas para la broncoscopia, que consistieron en la evaluación de un nódulo pulmonar indeterminado, resultados radiológicos anormales que no sugirieran infección, o el examen por tos crónica; (2) ausencia de infección pulmonar aguda, demostrada por resultados negativos en las pruebas microbiológicas y falta de cualquier signo de inflamación o infección sistémica en el momento de la broncoscopia; y (3) no haber utilizado antibióticos en los 30 días previos a la participación. Los controles se reclutaron de pacientes sometidos a broncoscopia por otras afecciones pulmonares clínicamente indicadas, en lugar de voluntarios sanos, reflejando así la población de la cual habitualmente se obtiene BAL clínicamente indicado en el centro del estudio.
Los criterios de exclusión incluyeron: (1) malignidad y/o quimioterapia citotóxica en los últimos seis meses; (2) seropositividad para VIH; (3) trasplante de órganos sólidos o trasplante de células madre hematopoyéticas; (4) terapia inmunosupresora sostenida que implique una dosis diaria equivalente de prednisona > 0,3 mg/kg durante más de 3 semanas; (5) mujeres embarazadas o en lactancia; (6) antecedentes de diatesis hemorrágica o cualquier contraindicación para broncoscopia flexible; (7) enfermedad respiratoria crónica severa preexistente (enfermedad pulmonar obstructiva crónica [EPOC] según la Iniciativa Global para la Enfermedad Pulmonar Obstructiva Crónica [GOLD], estadio III/IV; fibrosis pulmonar idiopática; asma persistente severa o bronquiectasias clínicamente significativas); (8) tuberculosis activa; (9) participación en otro estudio intervencional en los últimos 30 días. Se realizó lavado broncoalveolar durante la primera broncoscopia, dentro de las 72 horas posteriores al ingreso por infecciones, y dentro de las 48 horas posteriores a la intubación en pacientes con ventilación mecánica y enfermedades infecciosas.
Los cálculos del tamaño de la muestra se centraron en el resultado celular principal: el número de neutrófilos en el LBA entre los grupos. Con respecto al tamaño del efecto previsto de 0,6 DE entre vías respiratorias infectadas y no infectadas, y considerando un margen de abandono del 15 %, se requirió un tamaño mínimo de muestra de 45 sujetos por grupo para un α bilateral = 0,05 con un poder estadístico del 80 %. De 156 candidatos elegibles, 27 no cumplieron con los requisitos de elegibilidad (18 no cumplieron con los criterios de inclusión, principalmente por no poder identificar un agente patógeno responsable de su enfermedad dentro del período establecido o por falta de acuerdo con el procedimiento en sí; 9 incumplieron los criterios de exclusión, principalmente debido a malignidad activa o antecedentes de administración de quimioterapia, uso crónico de corticosteroides, enfermedad pulmonar crónica grave preexistente o diatesis hemorrágica). Los 129 participantes restantes se sometieron a broncoscopia. Posteriormente, 17 sujetos fueron retirados del estudio (7 por rendimiento inadecuado de LBA < 40 %; 4 por retiro del consentimiento durante la hospitalización; y 6 por reclasificación de los sujetos tras la evaluación microbiológica final: 4 candidatos del grupo control diagnosticados con infección oculta según los resultados del cultivo; y 2 candidatos del grupo infección que no presentaron evidencia de infección tras el análisis del cultivo). El grupo final de sujetos del estudio comprendió 112 participantes: 64 con infección y 48 sin infección.
Recolección y preprocesamiento del LBA
Se realizó broncoscopia por operadores capacitados siguiendo el protocolo institucional estándar, conforme con las recomendaciones técnicas publicadas para el lavado broncoalveolar12. Se utilizó lidocaína de forma tópica, y los pacientes recibieron sedación intravenosa moderada mientras se administraba oxígeno suplementario para mantener la saturación de oxígeno (SpO₂) >92 %. Se insertó el broncoscopio flexible y se posicionó en el segmento broncopulmonar con anomalías radiológicas. En los casos en que las anomalías radiológicas eran difusas o no localizables, se eligieron los segmentos del lóbulo medio derecho o del segmento lingular izquierdo. Tras cada instilación, se realizó succión manual suave con presión negativa < 100 mmHg. Se administraron en total de tres a cinco instilaciones de solución salina estéril a 37 °C, de 20–40 mL cada una, para un volumen total de 100–150 mL. La proporción de recuperación se estimó como el volumen recuperado en relación con el volumen instilado, y solo se consideraron para el análisis las muestras con una recuperación ≥40 %.
Las muestras recopiladas a partir del segundo lavado y los subsiguientes se agruparon para reducir la contaminación bronquial y obtener una muestra enriquecida en alvéolos. El concentrado de líquido de lavado broncoalveolar (BALF) se mantuvo sobre hielo y se procesó dentro de las 2 h posteriores a su recolección, teniendo en cuenta hallazgos recientes sobre los cambios en la composición celular y los recuentos diferenciales tras el almacenamiento prolongado y el calentamiento de muestras humanas de BALF13. Después de la filtración mediante una malla de nailon estéril con un tamaño de poro de 70 micrones para eliminar el moco, realizada únicamente mediante paso mecánico sin añadir agentes mucolíticos o anticonglomerantes como el ditiotreitol, ya que pruebas piloto mostraron que esto era suficiente para eliminar los coágulos de moco sin pérdida medible de la integridad de los epítopos superficiales en pasos posteriores, el fluido se sometió a centrifugación a 400 × g durante 10 min a 4 °C. Se decantó el sobrenadante y se distribuyó en viales criogénicos de baja unión a proteínas para el análisis de citoquinas, y las células se suspendieron en solución salina tamponada con fosfato (PBS) con 0,5 % de albúmina sérica bovina y 2 mM de ácido etilendiaminotetraacético (EDTA), reservando una parte para citología y citometría de flujo. Las muestras de citoquinas se congelaron en nitrógeno líquido, se mantuvieron a –80 °C y no se sometieron a más de un ciclo de congelación-descongelación antes de su uso.
Citología de neutrófilos y cuantificación mediante citometría de flujo
En cuanto a la citopatología, se aplicaron 100 µL de la suspensión celular de LBA en una portaobjetos utilizando un citocentrífugo (véase la Tabla de Materiales) a 300 ×g durante 5 min, se dejó secar al aire y se tiñó con Wright-Giemsa. Un tecnólogo de hematología realizó los recuentos diferenciales sin conocimiento del historial médico del paciente, contando un mínimo de 400 células por portaobjetos y calculando los porcentajes de macrófagos, linfocitos, neutrófilos y eosinófilos. Se consideró el porcentaje de neutrófilos en el LBA como la medida principal del resultado en citopatología, dada la evidencia actual que sugiere su utilidad para distinguir enfermedades pulmonares infecciosas de las no infecciosas en pacientes con infiltrados pulmonares14.
El protocolo de cuantificación mediante citometría de flujo se basó en una validación reciente del protocolo de diferenciación de leucocitos en líquido de lavado broncoalveolar (BALF) con cuatro colores15 e incluyó marcadores adicionales para identificar neutrófilos. Específicamente, aproximadamente 5 × 105 células aisladas del BALF se bloquearon con un reactivo bloqueador de receptores Fc humanos (ver Tabla de Materiales) durante 10 min a 4 °C, y luego se marcaron durante 20 min a 4 °C con anticuerpos contra CD45 (clon HI30), CD66b (clon G10F5), CD15 (clon HI98), CD16 (clon 3G8), CD11b (clon ICRF44), CD14 (clon M5E2) y HLA-DR (clon L243), además de un colorante vital fijable (ver Tabla de Materiales). Para la lisis de eritrocitos, se utilizó una solución comercial de lisis de eritrocitos (ver Tabla de Materiales) durante 10 min a temperatura ambiente. Tras la lisis de eritrocitos, las células se lavaron dos veces con PBS suplementado con 2% de suero bovino fetal y se resuspendieron en tampón de citometría, luego se analizaron mediante un citómetro de flujo multicolor (ver Tabla de Materiales). La compensación se realizó utilizando perlas de compensación marcadas individualmente (ver Tabla de Materiales), y se incluyeron tubos de control con menos un fluorocromo y tubos sin tinción para cada panel.
Se adquirieron al menos 1 × 105 eventos CD45⁺ por cada muestra. Primero se seleccionaron los eventos únicos en el gráfico FSC-A/FSC-H, seguidos por los leucocitos CD45⁺ viables (negativos para el colorante de viabilidad), y posteriormente se excluyeron los monocitos-macrófagos CD14⁺HLA-DR⁺ para seleccionar los granulocitos CD66b⁺ con alta dispersión lateral. Los neutrófilos se definieron entonces dentro de la población de granulocitos como eventos CD15⁺CD16⁺CD11b⁺. Esta puerta primaria de granulocitos CD66b⁺ se trazó de forma amplia alrededor de todos los eventos con alta dispersión lateral, en lugar de limitarse al contorno CD15⁺CD16⁺CD11b⁺, de modo que los neutrófilos activados CD16-débiles y CD11b-altos, así como otras formas de granulocitos menos maduras, se mantuvieran dentro del recuento total automatizado de neutrófilos; la definición más estricta CD15⁺CD16⁺CD11b⁺ se aplicó únicamente para caracterizar el subconjunto de neutrófilos maduros, que se informó por separado en los análisis de agrupamiento posteriores. La concentración de neutrófilos en el LBA (células/µL) se estimó en función del porcentaje de neutrófilos y del recuento total de células nucleadas determinado mediante un analizador hematológico automatizado (véase Tabla de materiales). Las matrices de compensación, la estrategia de selección de eventos y los informes individuales de control de calidad de las muestras se archivaron, y los analistas permanecieron ciegos respecto a la asignación de los participantes a los grupos durante toda la adquisición y el análisis.
Perfilado multiplex de citocinas inflamatorias
Los sobrenadantes de LBA se descongelaron una vez sobre hielo húmedo, seguido de una breve vortexación y clarificación mediante centrifugación a 10.000 × g durante 5 min a 4 °C inmediatamente antes del ensayo. La medición de 13 analitos en tres categorías funcionales se realizó mediante una plataforma de inmunoensayo multiplex basada en microesferas, utilizando un panel personalizado de microesferas magnéticas (ver Tabla de materiales). El grupo proinflamatorio incluyó IL-1β, IL-6, IL-17A y TNF-α; el grupo de quimiocinas incluyó CXCL8/IL-8, CXCL1/GRO-α, CXCL10/IP-10, CCL2/MCP-1 y CCL3/MIP-1α; y el grupo regulador incluyó IL-10, IFN-γ, IL-1RA y G-CSF. La selección de analitos se basó en aquellos que con mayor frecuencia se han reportado como informativos en análisis multiplex de citocinas en LBA o suero en estudios de neumonía/infección respiratoria grave7,11. El factor de dilución (1:2 en tampón de ensayo) se determinó a partir de un experimento piloto realizado con un subconjunto representativo de muestras infectadas y controles para asegurar que >El 90 % de las mediciones cayó dentro de la curva estándar.
Las curvas estándar se determinaron utilizando diluciones seriadas de siete puntos de los estándares recombinantes del fabricante, preparadas por duplicado en cada placa, incluyendo un pozo en blanco para restar el fondo. El ajuste de la curva se realizó mediante regresión logística de cinco parámetros (5PL) utilizando el software propietario de adquisición y análisis del instrumento (ver Tabla de Materiales). El límite inferior de cuantificación (LLOQ) se calculó como la cantidad más baja del estándar que produjo una recuperación por retroajuste entre el 80 % y el 120 % y un coeficiente de variación intramuestra (CV) ≤ 20 %. Los valores que estuvieron por debajo del LLOQ se sustituyeron por LLOQ/√2. Se prepararon dos muestras agrupadas de control de calidad de líquido de lavado broncoalveolar (BALF), cada una proveniente de un paciente infectado representativo y de un sujeto no afectado, se ali-cuotaron y se añadieron a cada placa para monitorear el rendimiento a largo plazo. Los estándares de aseguramiento de calidad se basaron en recomendaciones de buenas prácticas derivadas de un programa de competencia de ensayos multiplex basados en microesferas en múltiples centros durante 12 años16, incluyendo requisitos de CV intramuestra ≤ 10 % en pozos por duplicado, CV interplaca ≤ 20 % para el control de calidad y un recuento mínimo de microesferas de ≥ 50 por analito por pozo. Las placas que no cumplieron con los criterios anteriores se volvieron a analizar. El análisis se realizó en un solo lote dentro de las cuatro semanas, utilizando un solo lote del kit, con una intercalación aleatoria de las muestras de infección y de control.
Secuencia de análisis de visualización
Los archivos estándar de citometría de flujo (FCS) se analizaron primero utilizando software comercial de citometría de flujo (véase la Tabla de Materiales) para compensar y seleccionar eventos únicos vivos CD45⁺, y luego se exportaron para un análisis imparcial en R v4.4. Para cada muestra, se redujeron aleatoriamente a 10.000 eventos CD45⁺, que se concatenaron, y las intensidades de señal se transformaron utilizando la función arcsinh con un coefactor de 150. La reducción de dimensionalidad y las proyecciones UMAP se realizaron utilizando el paquete uwot con 15 vecinos más cercanos y una distancia mínima de 0,1, junto con embeddings t-SNE paralelos utilizando el paquete Rtsne con una perplejidad de 3010. Los grupos fenotípicos se identificaron mediante FlowSOM con 100 nodos de mapas autoorganizados y se dividieron posteriormente en 12–15 grupos mediante metaagrupamiento, utilizando un enfoque basado en optimización que recientemente ha sido validado para datos de citometría de alta dimensión17. Cada grupo se anotó según su expresión mediana de marcadores de linaje y activación, y se revisó en comparación con las compuertas manuales.
Para el análisis de datos multiplex de citocinas, los valores brutos se transformaron primero mediante log₁₀ y luego se estandarizaron entre analitos mediante puntuaciones z antes de realizar análisis adicionales. Tanto los analitos como los pacientes se sometieron a agrupamiento jerárquico utilizando la distancia euclidiana y el enlace Ward.D2; los resultados se mostraron en un mapa de calor con filas y columnas agrupadas. Este enfoque se utilizó para detectar módulos de citocinas coexpresadas y clasificar subconjuntos de pacientes en cohortes inflamatorias multiplex9. Se generaron matrices de correlación para todos los biomarcadores celulares y solubles utilizando correlaciones de rangos de Spearman por pares. Los mapas de calor de las matrices de correlación se representaron junto con una red en la que los nodos corresponden a biomarcadores y los enlaces indican correlaciones que resistieron la corrección por comparaciones múltiples mediante el método de Benjamini–Hochberg (BH). El ancho de los enlaces se escaló según la magnitud del coeficiente de correlación |ρ|, y solo se incluyeron enlaces si |ρ| ≥ 0,30 y p < 0,05 tras la corrección. Se realizó un análisis de componentes principales (PCA) sobre la matriz estandarizada que contenía células y citocinas utilizando la función prcomp() con escalamiento de varianza unitaria. Los participantes se proyectaron sobre los dos primeros componentes principales, superponiendo los grupos de infección y control. Un submuestreo y estandarización consistentes de todos los análisis garantizaron la comparabilidad entre las lecturas de biomarcadores celulares y solubles en todas las visualizaciones.
Análisis estadístico
Las variables continuas se describieron mediante la mediana (rango intercuartílico, IQR) cuando no presentaban distribución normal según la prueba de Shapiro-Wilk, y mediante la media (desviación estándar, DE) en caso contrario. Las variables categóricas se describieron usando frecuencias y proporciones. Las comparaciones entre grupos se realizaron utilizando la prueba de Mann-Whitney U prueba cuando las variables continuas tenían distribuciones no normales, la de Welch t-prueba cuando las variables continuas presentaban distribución normal, y la prueba de chi-cuadrado o la prueba exacta de Fisher cuando correspondía para variables categóricas. Las asociaciones entre variables continuas se evaluaron mediante el coeficiente de correlación de rangos de Spearman. Para tener en cuenta la tasa de descubrimientos falsos (FDR) en los análisis multiplex de citocinas y de correlación, los valores p se ajustaron utilizando el método de Benjamini-Hochberg con un umbral de q < 0.05. Bilateral p < Se consideró estadísticamente significativo un valor de p < 0,05 en todos los análisis primarios. Siempre que fue posible, se mantuvieron y analizaron los casos con datos faltantes en variables individuales mediante análisis de casos completos para la variable afectada, mientras que se realizaron imputaciones múltiples utilizando ecuaciones encadenadas (m = 20) para llevar a cabo análisis de sensibilidad. Se generaron curvas de característica operativa del receptor y el área correspondiente bajo la curva (AUC) para citocinas individuales utilizando el paquete pROC dentro de la cohorte completa del estudio; no se utilizó una cohorte de validación externa independiente. Los análisis de sensibilidad preespecificados repitieron las comparaciones primarias entre grupos tras excluir a los participantes que habían recibido antibióticos en las 72 horas previas a la broncoscopia, y tras estratificar el grupo de infección según el estado de ventilación mecánica y según el subtipo de neumonía (adquirida en la comunidad, adquirida en el hospital o asociada al ventilador). Los análisis se realizaron utilizando los paquetes rstatix, FlowSOM, uwot, ComplexHeatmap, igraph, pROC y mice en R v4.4. Las figuras se generaron utilizando software comercial para gráficos científicos (véase Tabla de materiales).