Este estudo foi aprovado pelo comitê de ética institucional do Hospital Popular da Região Autônoma Zhuang de Guangxi (número de aprovação IRB: KY-ZC-2023-035). A exigência de consentimento informado foi dispensada porque esta análise retrospectiva utilizou dados clínicos desidentificados coletados durante o atendimento de rotina.
Desenho do estudo e local
Este foi um estudo observacional retrospectivo de coorte único, conduzido na Unidade de Terapia Intensiva Pediátrica (UTIP) do Hospital Popular da Região Autônoma Zhuang de Guangxi. Foram avaliadas todas as admissões consecutivas na UTIP entre 1º de janeiro de 2022 e 31 de dezembro de 2025. Os dados clínicos e laboratoriais foram extraídos do sistema de prontuário eletrônico, do sistema de informações laboratoriais (SIL) e do sistema de informações de terapia intensiva.
O horário de admissão na UTI pediátrica foi definido como tempo zero para todas as janelas de exposição, biomarcadores, suporte orgânico e desfechos. Antes da extração dos dados, os critérios de elegibilidade, janelas de coleta de biomarcadores, definições de desfechos, métricas de gravidade e plano de análise estatística foram previamente especificados e implementados usando um fluxo de trabalho fixo de extração (especificação de extração do LIS versão 1.0; dicionário clínico de dados versão 1.0). Após a extração, todos os registros foram desidentificados. Dois pesquisadores verificaram independentemente os carimbos de tempo e a consistência das fontes para os valores de biomarcadores, variáveis relacionadas à avaliação pediátrica de falência orgânica sequencial (pSOFA), registros de terapia vasoativa e eventos de suporte orgânico. As discrepâncias foram resolvidas mediante revisão conjunta dos registros originais. A montagem da coorte seguiu etapas previamente especificadas, incluindo triagem, confirmação de elegibilidade, confirmação da janela de biomarcadores e revisão da completude dos desfechos. O fluxo final da coorte é apresentado na Figura 1.
Participantes e critérios de elegibilidade
Crianças com idade entre 1 mês e 14 anos admitidas na UTI pediátrica com sepse eram elegíveis. Os casos candidatos foram identificados a partir dos registros de admissão na UTI pediátrica e de diagnósticos relacionados à sepse, e posteriormente confirmados mediante revisão manual dos prontuários. A sepse foi definida como infecção suspeita ou confirmada com disfunção orgânica documentada durante o atendimento na UTI pediátrica e confirmada por avaliação baseada no pSOFA.
Para garantir exposição comparável aos biomarcadores precoces, exigia-se que os pacientes incluídos tivessem a proteína com afinidade por heparina (HBP) e a proteína C-reativa (PCR) medidas nas primeiras 24 horas após a admissão na UTI pediátrica (UTIP) e apresentassem dados suficientes para o cálculo da pSOFA. Os critérios de inclusão foram idade entre 1 mês e 14 anos, admissão na UTIP com sepse, disponibilidade de HBP e PCR entre 0 e 24 horas, e dados completos para a avaliação do desfecho primário, incluindo pSOFA basal, avaliação da pSOFA às 72 horas e acompanhamento de sobrevida. Os critérios de exclusão foram ausência de HBP ou PCR basal, permanência na UTIP inferior a 24 horas, transferência de outra UTIP com mais de 24 horas de internação prévia, admissão pós-transplante ou pós-cirurgia cardíaca, e diagnóstico de sepse secundária a infecção do trato urinário. <24 h, com variáveis-chave ausentes necessárias para a adjudicação do desfecho ou com imunossupressão grave susceptível de alterar substancialmente os perfis dos biomarcadores inflamatórios. Imunossupressão grave foi definida como neoplasia hematológica ativa, uso de quimioterapia ou terapia imunossupressora relacionada a transplante, terapia prolongada com corticosteroides sistêmicos em dose imunossupressora, imunodeficiência primária ou outro tratamento imunossupressor intensivo documentado antes ou durante a admissão na UTI pediátrica. Para admissões repetidas na UTI pediátrica, somente a primeira admissão elegível foi incluída.
Definições dos desfechos e métricas de gravidade
O desfecho primário foi a progressão precoce da disfunção orgânica nas primeiras 72 horas após a admissão na UTI pediátrica, definida como um aumento de 2 ou mais pontos no pSOFA do valor basal para o maior pSOFA registrado nas primeiras 72 horas (ΔpSOFA ≥ 2). O pSOFA basal foi o primeiro escore de pSOFA registrado entre 0 e 24 horas após a admissão na UTI pediátrica; caso múltiplos valores estivessem disponíveis, foi utilizado o escore mais próximo do horário da admissão. O valor de pSOFA às 72 horas foi definido como o pSOFA máximo registrado entre 0 e 72 horas.
A alta gravidade da doença foi definida a priori como um pSOFA máximo em 72 h ≥ 8, e esse limiar foi aplicado uniformemente nas análises categóricas. Esse limiar foi pré-especificado com base na prática de estratificação de gravidade da unidade e na sua interpretabilidade clínica para a identificação precoce de sepse pediátrica de alto risco. O pSOFA também foi analisado como uma medida contínua da carga de disfunção orgânica.
Os desfechos secundários foram choque séptico, ventilação mecânica invasiva, uso de drogas vasoativas, terapia de substituição renal contínua (CRRT), mortalidade na UTI pediátrica e mortalidade por todas as causas aos 28 dias. O choque séptico foi definido como falha circulatória persistente que exigia suporte vasoativo após a ressuscitação volêmica, com hipoperfusão tecidual documentada, incluindo lactato elevado e/ou achados anormais de perfusão periférica. Os desfechos relacionados ao suporte de órgãos foram codificados como variáveis binárias dentro de janelas precoces pré-especificadas. Para fornecer uma medida mais detalhada do suporte hemodinâmico, a pontuação vasoativa-inotrópica (VIS) também foi extraída ou calculada para pacientes que receberam terapia vasoativa nas primeiras 24 horas após a admissão na UTI pediátrica. A VIS foi utilizada como variável quantitativa de intensidade do suporte de órgãos nas análises de correlação e de associação com gravidade, enquanto o uso binário de terapia vasoativa foi mantido como um desfecho clinicamente interpretável de necessidade de tratamento.
Medição de biomarcadores e coleta de dados
HBP e PCR foram medidos no laboratório central utilizando fluxos de trabalho imunoensaios automatizados padronizados, conforme os procedimentos rotineiros de controle de qualidade clínica. Os testes de biomarcadores foram realizados em amostras de sangue venoso coletadas como parte do atendimento de rotina na UTI pediátrica e processadas de acordo com o fluxo pré-analítico padrão do laboratório. Os resultados de HBP foram relatados em ng/mL e os resultados de PCR em mg/L. O HBP foi medido utilizando um ensaio imunoturbidimétrico com reforço por partículas em uma plataforma automatizada de química clínica/imunoensaio. A PCR foi medida utilizando um ensaio imunoturbidimétrico de alta sensibilidade. O fluxo analítico e os procedimentos internos de controle de qualidade seguiram um protocolo de laboratório com controle de versão (versão 2022.1 do POP do laboratório) e permaneceram inalterados na lógica analítica principal durante o período do estudo.
Os registros laboratoriais foram exportados do LIS utilizando um fluxo de trabalho com consulta bloqueada e mapeamento de campos (especificação de extração do LIS versão 1.0), incluindo tanto o horário da coleta sanguínea quanto o horário do laudo. A HBP e a PCR basais foram definidas como os primeiros valores relatados dentro de 0–24 h após a admissão na UTI pediátrica. A HBP e a PCR do dia 2 foram definidas como as medições mais próximas de 48 h dentro de uma janela fixa de 36–60 h. Essa janela do dia 2 foi utilizada de forma consistente em todas as análises dinâmicas.
A variação absoluta de HBP foi calculada como HBP do dia 2 menos o HBP basal. A variação relativa de HBP (%) foi calculada como [(HBP do dia 2 − HBP basal) / HBP basal] × 100. Se o HBP basal estivesse abaixo do limite de quantificação do ensaio ou registrado como zero após a conversão do sistema, a variação relativa de HBP não foi calculada; esses pacientes foram mantidos nas análises utilizando HBP basal e variação absoluta de HBP, mas excluídos das análises de variação percentual. A razão HBP/albumina foi calculada dividindo-se o HBP pela albumina sérica medida dentro da janela de tempo correspondente e foi avaliada como um índice composto exploratório que reflete ativação endotelial e perda proteica relacionada ao vazamento capilar.
Além da HBP e da PCR, foram extraídos marcadores inflamatórios, relacionados à perfusão e à função orgânica, rotineiramente disponíveis, para comparação e modelos multivariáveis. Estes incluíram contagem de glóbulos brancos, procalcitonina, lactato, albumina, creatinina, contagem de plaquetas, variáveis de coagulação e variáveis de função hepática, quando disponíveis. Os valores de interleucina-6 (IL-6) e interleucina-8 (IL-8) foram extraídos do LIS quando tinham sido medidos como parte da avaliação clínica rotineira dentro das mesmas janelas precoces de avaliação. A dosagem de IL-6 e IL-8 não era obrigatória para inclusão na coorte e foi utilizada apenas nas análises em que estavam disponíveis. Essa regra foi aplicada para evitar a exclusão de pacientes que seriam elegíveis, exceto pelo fato de que a dosagem de citocinas não era rotineiramente realizada em todos os casos.
Dados demográficos, comorbidades, fonte da infecção, achados microbiológicos quando disponíveis, sinais vitais na admissão à UTI pediátrica, tratamento antimicrobiano inicial, medidas de ressuscitação e intervenções de suporte de órgãos foram extraídos do prontuário eletrônico e do sistema de informações de terapia intensiva utilizando o dicionário clínico de dados predefinido (versão 1.0). As variáveis de suporte de órgãos incluíram ventilação mecânica invasiva, terapia vasoativa, terapia de substituição renal contínua (CRRT), volume de bolo de fluidos e escore VIS nas primeiras 24 horas, quando aplicável. As variáveis necessárias para o pSOFA foram coletadas dentro de janelas previamente especificadas. Caso um componente do pSOFA gerado automaticamente fosse inconsistente com os valores originais, dois pesquisadores revisaram o prontuário e recalculariam o escore utilizando as regras de pontuação previamente especificadas. O controle de qualidade dos dados incluiu validação de carimbos de tempo, triagem de duplicatas, verificações de intervalo e distribuição, verificação da fonte de valores extremos e registro de auditoria das correções manuais.
Para preservar a reprodutibilidade evitando, ao mesmo tempo, terminologia excessivamente vinculada a marcas no texto principal do protocolo, a documentação dos ensaios foi registrada por classe do método, tipo de plataforma, ambiente de software, regras de faixa relatóvel, uso de calibradores e controles e procedimentos internos de controle de qualidade. Informações detalhadas sobre reagentes, calibradores, controles e software são fornecidas na Tabela de Materiais. O controle interno de qualidade foi realizado diariamente, e somente as séries laboratoriais que atendiam aos critérios de aceitação do laboratório quanto ao coeficiente de variação foram incluídas. Valores abaixo do limite inferior de quantificação foram ajustados para o limite inferior relatóvel. Valores acima do intervalo relatóvel superior foram ajustados para o limite superior relatóvel, exceto quando um resultado repetido e validado da mesma amostra estava disponível. As unidades laboratoriais foram verificadas quanto à consistência antes da análise; nenhuma conversão de unidades foi realizada, a menos que uma divergência de unidade no LIS fosse confirmada em comparação com o relatório original.
Dados ausentes e regras de tratamento de dados
As regras de tratamento de dados foram predefinidas. Pacientes com HBP basal, PCR basal, variáveis pSOFA basais ou variáveis pSOFA às 72 horas ausentes foram excluídos das análises do desfecho primário. Pacientes com HBP do Dia 2 ausente foram mantidos nas análises de biomarcadores basais e nos resumos gerais dos desfechos clínicos, mas excluídos das análises que exigiam índices dinâmicos de HBP. Pacientes sem medições de IL-6 ou IL-8 foram mantidos em todas as análises primárias e excluídos apenas das comparações de modelos que incluíam citocinas.
A análise de casos completos foi utilizada para cada modelo multivariável, e o tamanho da amostra específico para cada análise foi informado. A imputação múltipla não foi utilizada nas análises primárias porque o estudo focou-se nas trajetórias temporais dos biomarcadores observados e na disfunção orgânica. O potencial de viés de seleção causado pela exclusão de pacientes sem HBP basal, PCR basal ou dados completos do desfecho pSOFA foi considerado ao interpretar os resultados.
Valores laboratoriais extremos não foram excluídos unicamente por critérios estatísticos. Os valores atípicos foram verificados em relação aos registros originais e mantidos se confirmados como valores reais. Apenas erros evidentes de documentação ou de entrada de unidades foram corrigidos, e todas as correções foram registradas no log de auditoria com identificadores rastreáveis.
Análise estatística
Todas as análises seguiram um plano de análise estatística pré-especificado (versão 1.0 da especificação de análise). As variáveis contínuas foram avaliadas por inspeção visual de histogramas e pelo teste de Shapiro-Wilk. As variáveis com distribuição normal foram resumidas como média ± desvio padrão e comparadas utilizando testes t para amostras independentes. As variáveis com distribuição não normal foram resumidas como mediana (intervalo interquartílico) e comparadas utilizando o teste de Mann-Whitney U. As variáveis categóricas foram resumidas como contagens (porcentagens) e comparadas utilizando o teste qui-quadrado ou o teste exato de Fisher.
Para análises descritivas, os pacientes foram estratificados por quartis da HBP basal e separadamente pela presença ou ausência de progressão precoce da disfunção orgânica (ΔpSOFA ≥ 2). Os quartis da HBP basal foram utilizados para descrever padrões de dose-resposta ao longo da distribuição observada do biomarcador, sem assumir um ponto de corte clínico universal, pois limiares pediátricos validados de HBP para gravidade da sepse ainda são incertos e podem variar conforme a plataforma do ensaio, distribuição etária e janela de coleta. A idade foi analisada como uma covariável contínua nos modelos principais. Para abordar a heterogeneidade clínica relacionada à idade, análises subgrupadas pré-especificadas também foram realizadas utilizando estratos etários de <1 ano e ≥1 ano, com verificações exploratórias adicionais por categoria etária realizadas quando o tamanho da amostra permitiu.
As associações entre as métricas da PAF, incluindo PAF basal, PAF do dia 2, variação absoluta da PAF, variação relativa da PAF e razão PAF/albumina, e os indicadores de gravidade clínica foram avaliadas por meio de coeficientes de correlação de postos de Spearman. As variáveis que refletem suporte orgânico precoce, incluindo ventilação mecânica invasiva, terapia vasoativa, TRRCC, volume de bolo de fluido, dias de ventilação, tempo de internação na UTIP e VIS, foram consideradas como correlatos da gravidade clínica e do manejo nas análises descritivas e de associação.
O desempenho discriminativo para a progressão precoce da disfunção orgânica e gravidade elevada da doença foi avaliado utilizando curvas de característica de operação do receptor (ROC). A área sob a curva ROC (AUC) e os intervalos de confiança de 95% foram calculados para HBP, PCR, procalcitonina, lactato, albumina, pSOFA, razão HBP/albumina e modelos combinados. As AUCs foram comparadas utilizando o teste de DeLong. Os valores de corte ótimos foram identificados pelo índice de Youden, e sensibilidade, especificidade, valor preditivo positivo, valor preditivo negativo e razões de verossimilhança foram relatados.
A regressão logística multivariável foi utilizada para estimar associações independentes entre as métricas de HBP e os desfechos binários de gravidade, incluindo a progressão precoce da disfunção orgânica e a alta gravidade da doença. As covariáveis foram selecionadas a priori com base na relevância clínica e na disponibilidade de dados, e incluíram idade, categoria da fonte da infecção, pSOFA basal, lactato, albumina, procalcitonina, PCR, choque séptico na admissão e variáveis de suporte orgânico precoce, quando apropriado. Os modelos preditivos principais para os desfechos de gravidade precoce foram restritos às variáveis disponíveis na janela prescrita de avaliação inicial, a fim de reduzir a sobreposição temporal com a evolução dos desfechos. O valor prognóstico incremental foi avaliado mediante a adição sequencial das variáveis de HBP ao modelo clínico basal. Para a mortalidade por todas as causas aos 28 dias, foi utilizada a regressão de riscos proporcionais de Cox, e a suposição de riscos proporcionais foi verificada por meio dos resíduos de Schoenfeld. Antes do ajuste do modelo, foi avaliada a colinearidade entre os preditores candidatos, e o número de covariáveis incluídas em cada modelo foi limitado para preservar uma razão adequada de eventos por variável.
A melhoria do modelo após a adição da HBP foi avaliada utilizando métricas de discriminação, calibração e reclassificação. A calibração foi avaliada por meio de gráficos de calibração, intercepto e inclinação da calibração e escore de Brier, com o teste de Hosmer-Lemeshow relatado como uma verificação suplementar da calibração. A reclassificação foi quantificada usando o índice contínuo de melhoria na reclassificação líquida e a melhoria integrada na discriminação. As análises de sensibilidade pré-especificadas incluíram a redução da janela de coleta do Dia 2 de 36–60 h para 42–54 h, a repetição das análises com uma definição mais rigorosa de progressão da disfunção orgânica (ΔpSOFA ≥ 3), a exclusão de valores extremos de biomarcadores clinicamente atípicos após verificação da fonte e a repetição dos modelos principais após a exclusão dos pacientes sem medidas dinâmicas completas de HBP. As análises de sensibilidade foram relatadas na seção Resultados ou em tabelas suplementares, conforme apropriado.
As análises estatísticas foram realizadas no R versão 4.3.2 utilizando um ambiente de análise fixo (lançamento do ambiente de análise v1.0). Os principais pacotes foram pROC versão 1.18.5 para análises ROC, survival versão 3.5-7 para modelos de tempo até o evento e rms versão 6.8-0 para análises de calibração. A melhoria contínua da reclassificação líquida e a melhoria integrada da discriminação foram calculadas utilizando rotinas validadas implementadas no conjunto predefinido de scripts de análise (conjunto de pacotes de scripts v1.0). Todos os testes foram bicaudais e um valor de p < 0,05 foi considerado estatisticamente significativo.