População do estudo e ética
Este foi um estudo observacional prospectivo de centro único conduzido em um hospital universitário terciário entre janeiro de 2026 e abril de 2026. O protocolo foi aprovado sob o número de aprovação 2026HL-038, e o consentimento informado por escrito foi obtido de todos os participantes antes de qualquer procedimento do estudo.
Adultos com idade entre 18 e 75 anos submetidos à broncoscopia flexível diagnóstica com lavado broncoalveolar (LBA) como parte do atendimento clínico de rotina foram recrutados consecutivamente e alocados em um grupo de infecção ou em um grupo controle sem infecção, com base nos diagnósticos clínicos e microbiológicos finais. Uma visão geral do fluxo de trabalho de visualização integrada, desde o recrutamento dos pacientes até as saídas visuais integradas, é apresentada na Figura 1.

Figura 1: Visão geral do fluxo de trabalho integrado de visualização para o perfilamento baseado em LBA da recrutamento de neutrófilos e respostas inflamatórias de citocinas. Os cinco estágios resumem, da esquerda para a direita, o recrutamento prospectivo dos grupos de infecção e controle com broncoscopia e LBA, o pré-processamento do LBA em um pellet celular e uma fração de sobrenadante, medições celulares paralelas (citologia com Wright–Giemsa e citometria de fluxo multiparamétrica) e solúveis (ensaio imunológico com microesferas 13-plex), o pipeline unificado de visualização (UMAP com metaagrupamento FlowSOM, mapa térmico hierárquico, análise de componentes principais e rede de correlação de Spearman) e a saída visual integrada. Abreviações: LBA = lavado broncoalveolar; LBAF = fluido de LBA; UMAP = aproximação uniforme de variedades e projeção. Clique aqui para visualizar uma versão maior desta figura.
Os critérios de inclusão para o grupo de infecção foram os seguintes: (1) pacientes diagnosticados com pneumonia adquirida na comunidade (PAC), pneumonia associada a atendimento em saúde (PAS) ou pneumonia associada a ventilação mecânica (PAVM), conforme diretrizes da Sociedade de Doenças Infecciosas da América (IDSA) e da Sociedade Torácica Americana (ATS); (2) detecção de um ou mais patógenos por meio de cultura de LBA com ≥104 unidades formadoras de colônia (UFC)/mL, reação em cadeia da polimerase (PCR) respiratória multiplex, sequenciamento metagenômico de nova geração do LBA ou testes específicos de antígeno/anticorpo a partir de espécimes respiratórios/séricos pareados.
Os critérios para o grupo controle foram os seguintes: (1) razões não infecciosas para broncoscopia, consistindo na avaliação de um nódulo pulmonar indeterminado, resultados de radiologia anormais que não sugerem infecção ou exame para tosse crônica; (2) ausência de infecção pulmonar aguda, demonstrada por resultados negativos nos testes microbiológicos e ausência de quaisquer sinais de inflamação ou infecção sistêmica no momento da broncoscopia; e (3) ausência de uso de antibióticos nos 30 dias anteriores à participação. Os controles foram recrutados entre pacientes submetidos à broncoscopia por outras condições pulmonares clinicamente indicadas, e não entre voluntários saudáveis, refletindo a população da qual a LBA clinicamente indicada é rotineiramente obtida no centro do estudo.
Os critérios de exclusão incluíram: (1) neoplasia maligna e/ou quimioterapia citotóxica nos últimos seis meses; (2) soropositividade para HIV; (3) transplante de órgão sólido ou de células-tronco hematopoiéticas; (4) terapia imunossupressora prolongada envolvendo uma dose diária equivalente de prednisona > 0,3 mg/kg por mais de 3 semanas; (5) mulheres grávidas ou em período de amamentação; (6) histórico de diatese hemorrágica ou quaisquer contraindicações à broncoscopia flexível; (7) doença respiratória crônica grave preexistente (doença pulmonar obstrutiva crônica [DPOC] estágio III/IV da Iniciativa Global para Doença Pulmonar Obstrutiva Crônica [GOLD], fibrose pulmonar idiopática, asma persistente grave ou bronquiectasias clinicamente significativas); (8) tuberculose ativa; (9) participação em outro estudo intervencionista nos últimos 30 dias. A lavagem broncoalveolar (LBA) foi realizada durante a primeira broncoscopia, nas primeiras 72 horas após a admissão por infecções e nas primeiras 48 horas após a intubação em pacientes sob ventilação mecânica com infecções.
Os cálculos do tamanho da amostra foram focados no desfecho celular primário: o número de neutrófilos no LBA entre os grupos. Considerando um tamanho do efeito esperado de 0,6 DP entre vias aéreas infectadas e não infectadas e uma margem de 15% para perdas, foi necessário um tamanho mínimo de amostra de 45 indivíduos por grupo, com α bicaudal = 0,05 e poder estatístico de 80%. Dos 156 candidatos elegíveis, 27 não atenderam aos critérios de elegibilidade (18 não cumpriram os critérios de inclusão, principalmente devido à incapacidade de identificar um agente patogênico responsável pela doença dentro do período estipulado ou à falta de concordância com o próprio procedimento; 9 descumpriram os critérios de exclusão, principalmente por malignidade ativa ou histórico de administração de quimioterapia, uso crônico de corticosteroides, doença pulmonar crônica pré-existente grave ou diatese hemorrágica). Os 129 participantes restantes foram submetidos à broncoscopia. Posteriormente, 17 indivíduos foram excluídos do estudo (7 devido à coleta inadequada de LBA < 40%; 4 por retirada do consentimento durante a hospitalização; e 6 por reclassificação dos indivíduos após a avaliação microbiológica final – 4 candidatos do grupo controle diagnosticados com infecção oculta com base nos resultados de cultura; e 2 candidatos do grupo infecção que não apresentaram evidência de infecção após a análise de cultura). O grupo final de participantes do estudo foi composto por 112 indivíduos: 64 com infecção e 48 sem infecção.
Coleta e pré-processamento de LBA
A broncoscopia foi realizada por operadores qualificados seguindo o protocolo institucional padrão, conforme as recomendações técnicas publicadas para lavagem broncoalveolar12. A lidocaína foi usada topicamente, e os pacientes receberam sedação intravenosa moderada enquanto oxigênio suplementar era fornecido para manter a saturação de oxigênio (SpO₂) >92%. O broncoscópio flexível foi inserido e posicionado no segmento brônquico-pulmonar com anormalidades radiológicas. Nos casos em que as anormalidades radiológicas eram difusas ou não localizáveis, os segmentos do lobo médio direito ou da língula esquerda foram escolhidos. Após cada instilação, foi realizada sucção manual suave com pressão negativa < 100 mmHg. Foram administradas um total de três a cinco instilações de solução salina estéril a 37 °C, com 20–40 mL cada, totalizando um volume de 100–150 mL. A proporção de recuperação foi estimada como o volume recuperado em relação ao volume instilado, e somente amostras com recuperação ≥40% foram consideradas para análise.
As amostras coletadas a partir do segundo lavado e dos subsequentes foram agrupadas para reduzir a contaminação brônquica e obter uma amostra enriquecida em alvéolos. O pool de líquido de lavagem broncoalveolar (BALF) foi mantido sob gelo e processado dentro de 2 h após a coleta, levando em consideração descobertas recentes sobre alterações na composição celular e nas contagens diferenciais após armazenamento prolongado e aquecimento de amostras humanas de BALF13. Após a filtração utilizando uma malha de náilon estéril com porosidade de 70 micrômetros para remover o muco, realizada apenas por passagem mecânica, sem adição de agente mucolítico ou anticlumpante como o ditiotreitol, já que testes pilotos demonstraram que isso era suficiente para remover os coágulos de muco sem perda mensurável da integridade dos epítopos superficiais em etapas subsequentes, o fluido foi submetido à centrifugação a 400 × g durante 10 min a 4 °C. O sobrenadante foi decantado e fracionado em criotubos de baixa ligação proteica para análise de citocinas, e as células foram suspensas em solução salina tamponada com fosfato (PBS) contendo 0,5% de albúmina sérica bovina e 2 mM de ácido etilenodiaminotetraacético (EDTA), sendo parte delas mantida para citologia e citometria de fluxo. As amostras de citocinas foram congeladas em nitrogênio líquido, armazenadas a –80 °C e submetidas a não mais do que um ciclo de congelamento-descongelamento antes do uso.
Citologia de neutrófilos e quantificação por citometria de fluxo
Em termos de citopatologia, 100 µL da suspensão celular de LBA foi aplicado em lâmina de vidro utilizando um citocentrífugo (ver Tabela de Materiais) a 300 ×g durante 5 min, secado ao ar e corado com Wright-Giemsa. A contagem diferencial foi realizada por um tecnologista em hematologia sem conhecimento do histórico clínico do paciente, que contou um mínimo de 400 células por lâmina e calculou as porcentagens de macrófagos, linfócitos, neutrófilos e eosinófilos. As porcentagens de neutrófilos no LBA foram consideradas a medida primária de desfecho para citopatologia, dada a evidência atual que sugere sua utilidade na distinção entre doenças pulmonares infecciosas e não infecciosas em pacientes com infiltrados pulmonares14.
O protocolo de quantificação por citometria de fluxo baseou-se em uma validação recente do protocolo de diferenciação de leucócitos em fluido de lavado broncoalveolar (BALF) com quatro cores15 e incluiu marcadores adicionais para identificar neutrófilos. Especificamente, cerca de 5 × 105 células isoladas do BALF foram bloqueadas com um reagente bloqueador de receptor Fc humano (ver Tabela de Materiais) por 10 min a 4 °C, e então coradas por 20 min a 4 °C com anticorpos contra CD45 (clone HI30), CD66b (clone G10F5), CD15 (clone HI98), CD16 (clone 3G8), CD11b (clone ICRF44), CD14 (clone M5E2) e HLA-DR (clone L243), além de um corante fixável de viabilidade (ver Tabela de Materiais). Para a lise de eritrócitos, foi utilizada uma solução comercial de lise de eritrócitos (ver Tabela de Materiais) por 10 min à temperatura ambiente. Após a lise de eritrócitos, as células foram lavadas duas vezes com PBS suplementado com 2% de soro fetal bovino e ressuspendidas em tampão de citometria, sendo então analisadas em um citômetro de fluxo multicolorido (ver Tabela de Materiais). A compensação foi realizada utilizando microesferas de compensação monocromáticas (ver Tabela de Materiais), e foram incluídos tubos de fluorescência menos um e tubos não corados para cada painel.
Foram adquiridos pelo menos 1 × 105 eventos CD45⁺ para cada amostra. Primeiro, foram selecionados os singletos no gráfico FSC-A/FSC-H, seguidos pelos leucócitos viáveis (negativos para corante de viabilidade) CD45⁺, após o que foi realizada a exclusão de monócitos-macrófagos CD14⁺HLA-DR⁺ para selecionar os granulócitos CD66b⁺ com alta dispersão lateral. Os neutrófilos foram então definidos dentro da população de granulócitos como eventos CD15⁺CD16⁺CD11b⁺. Este primeiro gate de granulócitos CD66b⁺ foi desenhado de forma ampla ao redor de todos os eventos com alta dispersão lateral, em vez de ser restrito ao limite CD15⁺CD16⁺CD11b⁺, de modo que neutrófilos ativados CD16 fraco e CD11b alto, bem como outras formas de granulócitos menos maduras, permanecessem incluídos na contagem automatizada total de neutrófilos; a definição mais rigorosa CD15⁺CD16⁺CD11b⁺ foi aplicada apenas para caracterizar o subconjunto de neutrófilos maduros, relatado separadamente nas análises de agrupamento subsequentes. A contagem de neutrófilos no LBA (células/µL) foi estimada com base na porcentagem de neutrófilos e na contagem total de células nucleadas determinada por meio de um analisador hematológico automatizado (ver Tabela de Materiais). As matrizes de compensação, a estratégia de seleção de gates e os relatórios individuais de controle de qualidade das amostras foram arquivados, e os analistas permaneceram cegos quanto à atribuição dos grupos dos participantes durante toda a aquisição e análise.
Perfilagem multiplex de citocinas inflamatórias
Os sobrenadantes de LBA foram descongelados uma vez em gelo úmido, seguidos por vortexação breve e clarificação por centrifugação a 10.000 × g por 5 min a 4 °C imediatamente antes do ensaio. A medição de 13 analitos em três categorias funcionais foi realizada utilizando uma plataforma de ensaio imunoenzimático multiplex baseada em microesferas, com um painel personalizado de microesferas magnéticas (ver Tabela de Materiais). O grupo pró-inflamatório compreendeu IL-1β, IL-6, IL-17A e TNF-α; o grupo de quimiocinas compreendeu CXCL8/IL-8, CXCL1/GRO-α, CXCL10/IP-10, CCL2/MCP-1 e CCL3/MIP-1α; e o grupo regulador compreendeu IL-10, IFN-γ, IL-1RA e G-CSF. A escolha dos analitos baseou-se nos que são mais frequentemente relatados como informativos em análises multiplex de citocinas em LBA ou soro em estudos de pneumonia/infeção respiratória grave.7,11. O fator de diluição (1:2 em tampão de ensaio) foi determinado a partir de um experimento piloto realizado em um subconjunto representativo de amostras infectadas e controles para garantir que >90% das medições estiveram dentro da curva padrão.
As curvas padrão foram determinadas utilizando diluições seriadas de sete pontos dos padrões recombinantes do fabricante, preparadas em duplicata em cada placa, incluindo um poço em branco para subtrair o fundo. O ajuste da curva foi realizado por regressão logística de cinco parâmetros (5PL) por meio do software proprietário de aquisição e análise do instrumento (consulte a Tabela de Materiais). O limite inferior de quantificação (LIQ) foi calculado como a menor quantidade de padrão que gerou uma recuperação por retroajuste entre 80% e 120% e um coeficiente de variação intrarrodada (CV) ≤ 20%. Os valores que ficaram abaixo do LIQ foram substituídos por LIQ/√2. Duas amostras pool de fluido de lavagem broncoalveolar (BALF) para controle de qualidade, cada uma proveniente de um paciente infectado representativo e de um indivíduo não afetado, foram aliquotadas e adicionadas a cada placa para acompanhar o desempenho a longo prazo. Os padrões de garantia de qualidade basearam-se em recomendações de boas práticas derivadas de um programa de proficiência multicêntrico de ensaio multiplex baseado em microesferas com 12 anos de duração16, incluindo requisitos de CV intraensaio ≤ 10% em poços em duplicata, CV interplaca ≤ 20% para controle de qualidade e contagem mínima de microesferas de ≥ 50 por analito por poço. As placas que não atenderam aos critérios acima foram repetidas. A análise foi realizada em um único lote dentro de quatro semanas, utilizando um único lote do kit, com aleatorização intercalada das amostras de infecção e controle.
Fluxo de análise de visualização
Os arquivos no Padrão de Citometria de Fluxo (FCS) foram inicialmente analisados usando software comercial de citometria de fluxo (consulte a Tabela de Materiais) para compensação e seleção de singletos viáveis CD45⁺, e posteriormente exportados para análise imparcial no R v4.4. Para cada amostra, 10.000 eventos CD45⁺ foram subamostrados e concatenados, e as intensidades de sinal foram transformadas utilizando a função arcsinh com um fator de cofator igual a 150. A redução de dimensionalidade e as projeções UMAP foram realizadas usando o pacote uwot com 15 vizinhos mais próximos e uma distância mínima de 0,1, juntamente com incorporações t-SNE paralelas usando o pacote Rtsne com perplexidade 3010. Os agrupamentos fenotípicos foram identificados utilizando o FlowSOM com 100 nós de mapas autoorganizáveis e posteriormente divididos em 12–15 agrupamentos por meio de metaagrupamento, utilizando uma abordagem baseada em otimização que foi recentemente validada para dados de citometria de alta dimensionalidade17. Cada agrupamento foi anotado com base na expressão mediana de marcadores de linhagem e ativação e revisado em comparação com os gates manuais.
Para a análise de dados multiplex de citocinas, os valores brutos foram primeiramente transformados por log₁₀ e, em seguida, padronizados entre analitos por meio de escore z antes das análises adicionais. Tanto os analitos quanto os pacientes passaram por agrupamento hierárquico utilizando distância euclidiana e ligação Ward.D2; os resultados foram exibidos em um mapa de calor com agrupamento em linhas e colunas. Essa abordagem foi utilizada para detectar módulos de citocinas coexpressas e classificar subconjuntos de pacientes em coortes inflamatórias multiplex9. Matrizes de correlação foram geradas para todos os biomarcadores celulares e solúveis utilizando correlações de postos de Spearman pareadas. Mapas de calor das matrizes de correlação foram plotados ao lado de uma rede na qual os nós correspondem aos biomarcadores e as arestas indicam correlações que resistiram à correção para testes múltiplos pelo método de Benjamini–Hochberg (BH). A largura das arestas foi ajustada de acordo com a magnitude do coeficiente de correlação |ρ|, sendo incluídas apenas as arestas com |ρ| ≥ 0,30 e p < 0,05 após a correção. A análise de componentes principais (PCA) foi realizada na matriz padronizada contendo células e citocinas utilizando a função prcomp() com escala de variância unitária. Os participantes foram projetados nos dois primeiros componentes principais, com os grupos de infecção e controle sobrepostos. O subamostragem consistente e a padronização de todas as análises garantiram a comparabilidade entre as leituras de biomarcadores celulares e solúveis em todas as visualizações.
Análise estatística
As variáveis contínuas foram descritas utilizando a mediana (intervalo interquartil, IIQ) quando não normais de acordo com os testes de Shapiro–Wilk, e a média (desvio padrão, DP) caso contrário. As variáveis categóricas foram descritas por meio de frequências e proporções. As comparações entre grupos foram realizadas utilizando o teste de Mann–Whitney U teste quando as variáveis contínuas tinham distribuições não normais, o Welch t-test quando as variáveis contínuas apresentavam distribuição normal e o teste qui-quadrado/teste exato de Fisher quando aplicável para variáveis categóricas. As associações entre variáveis contínuas foram avaliadas por meio do coeficiente de correlação de postos de Spearman. Para controlar a taxa de falsa descoberta (FDR) nas análises multiplex de citocinas e de correlação, os valores de p foram ajustados utilizando o método de Benjamini–Hochberg com um limiar de q < 0,05. Bicaudal p < 0,05 foi considerado estatisticamente significativo em todas as análises primárias. Sempre que possível, os casos com dados ausentes em variáveis individuais foram mantidos e analisados por análise de casos completos para a variável afetada, enquanto imputações múltiplas usando equações encadeadas (m = 20) foram realizadas para conduzir análises de sensibilidade. Curvas de característica operacional do receptor e a área correspondente sob a curva (AUC) foram geradas para citocinas individuais utilizando o pacote pROC na coorte completa do estudo; uma coorte separada de validação externa não foi utilizada. Análises de sensibilidade pré-especificadas repetiram as comparações entre os grupos primários após a exclusão dos participantes que tinham recebido antibióticos nas 72 horas anteriores à broncoscopia, e após estratificar o grupo de infecção conforme o status de ventilação mecânica e o subtipo de pneumonia (adquirida na comunidade, adquirida no hospital ou associada à ventilação). As análises foram realizadas utilizando os pacotes rstatix, FlowSOM, uwot, ComplexHeatmap, igraph, pROC e mice no R v4.4. As figuras foram geradas utilizando software comercial de gráficos científicos (ver Tabela de Materiais).