Este estudio fue aprobado por el comité de ética institucional del Hospital Popular de la Región Autónoma Zhuang de Guangxi (número de aprobación IRB: KY-ZC-2023-035). Se dispensó el requisito de consentimiento informado porque este análisis retrospectivo utilizó datos clínicos desidentificados recopilados durante la atención rutinaria.
Diseño del estudio y entorno
Este fue un estudio observacional retrospectivo de cohorte unicéntrico realizado en la Unidad de Cuidados Intensivos Pediátricos (UCIP) del Hospital Popular de la Región Autónoma Zhuang de Guangxi. Se revisaron todos los ingresos consecutivos a la UCIP entre el 1 de enero de 2022 y el 31 de diciembre de 2025. Los datos clínicos y de laboratorio se extrajeron del sistema de historias clínicas electrónicas, del sistema de información de laboratorio (LIS) y del sistema de información de cuidados críticos.
El momento de ingreso en la UCI pediátrica se definió como el tiempo cero para todas las ventanas de exposición, biomarcadores, soporte orgánico y resultados. Antes de la extracción de datos, se especificaron previamente y se implementaron los criterios de elegibilidad, las ventanas temporales de muestreo de biomarcadores, las definiciones de los puntos finales, las métricas de gravedad y el plan de análisis estadístico, utilizando un protocolo fijo de extracción (especificación de extracción LIS versión 1.0; diccionario clínico de datos versión 1.0). Tras la extracción, todos los registros fueron desidentificados. Dos investigadores verificaron de forma independiente las marcas de tiempo y la coherencia con las fuentes para los valores de biomarcadores, las variables relacionadas con la evaluación pediátrica secuencial de fallo orgánico (pSOFA), los registros de terapia vasoactiva y los eventos de soporte orgánico. Las discrepancias se resolvieron mediante revisión conjunta de los registros originales. La conformación de la cohorte siguió pasos previamente especificados, incluyendo cribado, confirmación de elegibilidad, confirmación de la ventana de biomarcadores y revisión de la completitud de los puntos finales. El flujo final de la cohorte se muestra en Figura 1.
Participantes y criterios de elegibilidad
Eran elegibles los niños de 1 mes a 14 años admitidos en la UCI pediátrica con sepsis. Los casos candidatos se identificaron a partir de los registros de ingreso a la UCI pediátrica y de diagnósticos relacionados con sepsis, y luego se confirmaron mediante revisión manual de historias clínicas. La sepsis se definió como infección sospechada o confirmada con disfunción orgánica documentada durante la atención en la UCI pediátrica y respaldada por una evaluación basada en pSOFA.
Para garantizar una exposición comparable a biomarcadores tempranos, se requirió que los pacientes incluidos tuvieran medidos tanto la proteína con afinidad por heparina (HBP) como la proteína C reactiva (PCR) dentro de las 24 horas posteriores al ingreso en la UCI pediátrica (UCIP) y que contaran con datos suficientes para el cálculo de la puntuación pSOFA. Los criterios de inclusión fueron edad entre 1 mes y 14 años, ingreso en UCIP con sepsis, disponibilidad de HBP y PCR entre 0 y 24 horas, y datos completos para la evaluación del resultado principal, incluyendo la puntuación basal de pSOFA, la evaluación de pSOFA a las 72 horas y el seguimiento de supervivencia. Los criterios de exclusión fueron la ausencia de HBP o PCR basal, estancia en UCIP < 24 h, muerte o alta dentro de las primeras 24 h del ingreso, y falta de consentimiento informado. <24 h, falta de variables clave necesarias para la adjudicación del desenlace, o inmunosupresión grave que probablemente altere sustancialmente los perfiles de biomarcadores inflamatorios. Se consideró inmunosupresión grave la presencia de una neoplasia hematológica activa, administración de quimioterapia o tratamiento inmunosupresor relacionado con trasplante, terapia prolongada con corticosteroides sistémicos a dosis inmunosupresoras, inmunodeficiencia primaria, u otro tratamiento inmunosupresor intensivo documentado antes o durante la admisión en la UCI pediátrica. En caso de múltiples admisiones en la UCI pediátrica, solo se incluyó la primera admisión elegible.
Definiciones de resultados y métricas de gravedad
El punto final primario fue la progresión temprana de disfunción orgánica dentro de las 72 h posteriores al ingreso en la UCI pediátrica, definida como un aumento de 2 o más puntos en el pSOFA desde el valor basal hasta el pSOFA más alto registrado durante las primeras 72 h (ΔpSOFA ≥ 2). El pSOFA basal fue la primera puntuación de pSOFA registrada entre 0 y 24 h tras el ingreso en la UCI pediátrica; si había múltiples valores disponibles, se utilizó la puntuación más cercana al momento del ingreso. El valor de pSOFA a las 72 h se definió como el pSOFA máximo registrado entre las 0 y 72 h.
Se definió a priori una alta gravedad de la enfermedad como un pSOFA máximo a las 72 h ≥ 8, y este umbral se aplicó de forma uniforme en los análisis categóricos. Este umbral se especificó previamente basándose en la práctica de estratificación de gravedad de la unidad y en su interpretabilidad clínica para la identificación temprana de sepsis pediátrica de alto riesgo. El pSOFA también se analizó como una medida continua de la carga de disfunción orgánica.
Los resultados secundarios fueron shock séptico, ventilación mecánica invasiva, uso de fármacos vasoactivos, terapia de reemplazo renal continuo (CRRT), mortalidad en la UCI pediátrica y mortalidad por todas las causas a los 28 días. Se definió shock séptico como una falla circulatoria persistente que requiere soporte vasoactivo tras la reanimación con líquidos, con hipoperfusión tisular documentada, incluyendo niveles elevados de lactato y/o hallazgos anormales de perfusión periférica. Los resultados relacionados con el soporte orgánico se codificaron como variables binarias dentro de ventanas tempranas preespecificadas. Para obtener una medida más detallada del soporte hemodinámico, también se extrajo o calculó la puntuación vasoactiva-inotrópica (VIS) en pacientes que recibieron terapia vasoactiva durante las primeras 24 h tras el ingreso a la UCI pediátrica. La VIS se utilizó como una variable cuantitativa de intensidad del soporte orgánico en los análisis de correlación y asociación con la gravedad, mientras que el uso binario de la terapia vasoactiva se mantuvo como un resultado clínicamente interpretable de necesidad de tratamiento.
Medición de biomarcadores y recopilación de datos
La HBP y la PCR se midieron en el laboratorio central utilizando flujos de trabajo inmunoenzimáticos automatizados y estandarizados, de acuerdo con procedimientos rutinarios de control de calidad clínica. Las pruebas de biomarcadores se realizaron con muestras de sangre venosa recogidas como parte de la atención rutinaria en la UCI pediátrica y procesadas según el flujo de trabajo preanalítico estándar del laboratorio. Los resultados de HBP se informaron en ng/mL y los de PCR en mg/L. La HBP se midió mediante un ensayo inmunoturbidimétrico potenciado con partículas en una plataforma automatizada completa de química clínica/inmunoensayo. La PCR se midió utilizando un ensayo inmunoturbidimétrico de alta sensibilidad. El flujo de trabajo analítico y los procedimientos internos de control de calidad siguieron un protocolo de laboratorio con control de versiones (SOP de laboratorio versión 2022.1) y no cambiaron en su lógica analítica fundamental durante el período del estudio.
Los registros de laboratorio se exportaron desde el LIS mediante una consulta bloqueada y un flujo de trabajo de mapeo de campos (especificación de extracción del LIS versión 1.0), que incluyó tanto la hora de extracción sanguínea como la hora del informe. La HBP y la PCR basales se definieron como los primeros valores informados dentro de las 0–24 h posteriores al ingreso a la UCI pediátrica. La HBP y la PCR del día 2 se definieron como las mediciones más cercanas a las 48 h dentro de una ventana fija de 36–60 h. Esta ventana del día 2 se utilizó de forma consistente en todos los análisis dinámicos.
El cambio absoluto en HBP se calculó como el HBP del día 2 menos el HBP basal. El cambio relativo en HBP (%) se calculó como [(HBP del día 2 − HBP basal) / HBP basal] × 100. Si el HBP basal estaba por debajo del límite de cuantificación del ensayo o se registraba como cero tras la conversión del sistema, no se calculó el cambio relativo en HBP; estos pacientes se mantuvieron en los análisis que utilizaban el HBP basal y el cambio absoluto en HBP, pero se excluyeron de los análisis de cambio porcentual. La relación HBP/albumina se calculó dividiendo el HBP por la albúmina sérica medida dentro del intervalo de tiempo correspondiente y se evaluó como un índice compuesto exploratorio que refleja la activación endotelial y la pérdida de proteínas relacionada con la fuga capilar.
Además de HBP y PCR, se extrajeron marcadores inflamatorios, relacionados con la perfusión y de función orgánica, disponibles de forma rutinaria, para comparación y modelos multivariables. Estos incluyeron recuento de glóbulos blancos, procalcitonina, lactato, albúmina, creatinina, recuento de plaquetas, variables de coagulación y variables de función hepática cuando estaban disponibles. Los valores de interleucina-6 (IL-6) e interleucina-8 (IL-8) se extrajeron del sistema de información de laboratorio (LIS) cuando se habían medido como parte de la evaluación clínica rutinaria dentro de las mismas ventanas tempranas de evaluación. No se requirieron IL-6 ni IL-8 para la inclusión en la cohorte y solo se utilizaron en los análisis donde estuvieron disponibles. Esta regla se aplicó para evitar excluir a pacientes que de otro modo serían elegibles únicamente porque la prueba de citocinas no se realizaba de forma rutinaria en todos los casos.
Se extrajeron del registro médico electrónico y del sistema de información de cuidados críticos, mediante el uso del diccionario clínico de datos predefinido (versión 1.0), datos demográficos, comorbilidades, fuente de infección, hallazgos microbiológicos cuando estuvieron disponibles, signos vitales al ingreso en la UCI pediátrica, tratamiento antimicrobiano temprano, medidas de resucitación e intervenciones de soporte orgánico. Las variables de soporte orgánico incluyeron ventilación mecánica invasiva, terapia vasoactiva, terapia de reemplazo renal continuo (CRRT), volumen de bolo de fluidos y el índice de soporte vasoactivo (VIS) durante las primeras 24 horas cuando fue aplicable. Las variables necesarias para el pSOFA se recopilaron dentro de ventanas preespecificadas. Si un componente del pSOFA generado automáticamente era inconsistente con los valores originales, dos investigadores revisaron la historia clínica y recalculaban la puntuación utilizando las reglas de puntuación preestablecidas. El control de calidad de los datos incluyó validación de marcas de tiempo, detección de duplicados, verificación de rangos y distribuciones, confirmación de la fuente de valores extremos y registro de auditoría de correcciones manuales.
Para preservar la reproducibilidad y evitar al mismo tiempo formulaciones innecesarias específicas de una marca en el texto principal del protocolo, la documentación del ensayo se registró según la clase del método, el tipo de plataforma, el entorno de software, las reglas del rango de cuantificación, el uso de calibradores y controles, y los procedimientos internos de control de calidad. Se proporciona información detallada sobre reactivos, calibradores, controles y software en la Tabla de Materiales. El control interno de calidad se realizó diariamente, y solo se incluyeron las series de trabajo del laboratorio que cumplían con los criterios de aceptación del coeficiente de variación del laboratorio. Los valores por debajo del límite inferior de cuantificación se establecieron en el límite inferior reportable. Los valores por encima del rango reportable superior se establecieron en el límite superior reportable, a menos que estuviera disponible un resultado repetido validado de la misma muestra. Las unidades del laboratorio se verificaron para garantizar su consistencia antes del análisis; no se realizó ninguna conversión de unidades a menos que se confirmara una discrepancia en la unidad del LIS frente al informe original.
Datos faltantes y reglas de manejo de datos
Las reglas de manejo de datos se especificaron previamente. Se excluyeron de los análisis del punto final primario los pacientes con datos faltantes de HBP basal, PCR basal, variables pSOFA basales o variables pSOFA a las 72 h. Se mantuvieron en los análisis de biomarcadores basales y en los resúmenes generales de resultados clínicos a los pacientes con HBP del día 2 faltante, pero se excluyeron de los análisis que requerían índices dinámicos de HBP. Se mantuvieron en todos los análisis primarios a los pacientes sin mediciones de IL-6 o IL-8, y solo se excluyeron de las comparaciones de modelos que incluían citocinas.
Se utilizó un análisis de casos completos para cada modelo multivariable, y se informó el tamaño muestral específico del análisis. No se utilizó imputación múltiple en los análisis principales porque el estudio se centró en las trayectorias temporales de los biomarcadores observados y la disfunción orgánica. Se consideró el riesgo de sesgo de selección provocado por la exclusión de pacientes sin HBP basal, PCR basal o datos completos del punto final pSOFA al interpretar los hallazgos.
No se excluyeron valores de laboratorio extremos únicamente por motivos estadísticos. Los valores atípicos se verificaron frente a los registros originales y se mantuvieron si se confirmaron como valores reales. Únicamente se corrigieron errores evidentes de documentación o de entrada de unidades, y todas las correcciones se registraron en el historial de auditoría con identificadores trazables.
Análisis estadístico
Todos los análisis siguieron un plan de análisis estadístico preespecificado (versión 1.0 de la especificación de análisis). Las variables continuas se evaluaron mediante inspección visual de histogramas y la prueba de Shapiro-Wilk. Las variables con distribución normal se resumieron como media ± desviación estándar y se compararon utilizando pruebas t para muestras independientes. Las variables con distribución no normal se resumieron como mediana (rango intercuartílico) y se compararon utilizando la prueba de Mann-Whitney U. Las variables categóricas se resumieron como frecuencias (porcentajes) y se compararon mediante la prueba de chi-cuadrado o la prueba exacta de Fisher.
Para los análisis descriptivos, los pacientes se estratificaron por cuartiles de HBP basal y por separado según la presencia o ausencia de progresión temprana de disfunción orgánica (ΔpSOFA ≥ 2). Los cuartiles basales de HBP se utilizaron para describir patrones de dosis-respuesta a lo largo de la distribución observada del biomarcador sin asumir un punto de corte clínico universal, ya que los umbrales pediátricos validados de HBP para la gravedad de la sepsis aún son inciertos y pueden variar según la plataforma del ensayo, la distribución por edad y la ventana de muestreo. La edad se analizó como una covariable continua en los modelos principales. Para abordar la heterogeneidad clínica relacionada con la edad, también se realizaron análisis prespecificados de subgrupos utilizando estratos de edad de <1 año y ≥1 año, con comprobaciones exploratorias adicionales por categorías de edad cuando el tamaño de la muestra lo permitió.
Se evaluaron las asociaciones entre las métricas de PTH, incluyendo PTH basal, PTH del día 2, cambio absoluto de PTH, cambio relativo de PTH y relación PTH-albumina, y los indicadores de gravedad clínica mediante coeficientes de correlación de rangos de Spearman. Las variables que reflejan soporte orgánico temprano, incluyendo ventilación mecánica invasiva, terapia vasoactiva, TCRRT, volumen de bolo de fluidos, días de ventilador, duración de la estancia en la UCI pediátrica y el índice VIS, se consideraron como correlatos de gravedad clínica y manejo en los análisis descriptivos y de asociación.
Se evaluó el rendimiento discriminativo para la progresión temprana de disfunción orgánica y la gravedad elevada de la enfermedad mediante curvas de característica operativa del receptor (ROC). Se calcularon el área bajo la curva ROC (AUC) e intervalos de confianza del 95 % para HBP, PCR, procalcitonina, lactato, albúmina, pSOFA, la relación HBP-albúmina y modelos combinados. Las AUC se compararon utilizando la prueba de DeLong. Los valores óptimos de corte se identificaron mediante el índice de Youden, y se informaron sensibilidad, especificidad, valor predictivo positivo, valor predictivo negativo y razones de verosimilitud.
Se utilizó una regresión logística multivariable para estimar las asociaciones independientes entre las métricas de HBP y los resultados binarios de gravedad, incluyendo la progresión temprana de disfunción orgánica y la alta severidad de la enfermedad. Las covariables se seleccionaron a priori según su relevancia clínica y disponibilidad de datos, e incluyeron edad, categoría de fuente de infección, pSOFA basal, lactato, albúmina, procalcitonina, PCR, shock séptico al ingreso y, cuando correspondía, variables de soporte orgánico temprano. Los modelos de predicción principales para los resultados de gravedad temprana se limitaron a las variables disponibles dentro de la ventana preestablecida de evaluación temprana, con el fin de reducir la superposición temporal con la evolución del resultado. El valor pronóstico incremental se evaluó mediante la adición secuencial de variables de HBP al modelo clínico basal. Para la mortalidad por todas las causas a los 28 días, se utilizó una regresión de riesgos proporcionales de Cox, y se verificó el supuesto de riesgos proporcionales mediante los residuos de Schoenfeld. Antes del ajuste del modelo, se evaluó la colinealidad entre los predictores candidatos, y se limitó el número de covariables incluidas en cada modelo para mantener una relación adecuada de eventos por variable.
La mejora del modelo tras la adición de HBP se evaluó utilizando métricas de discriminación, calibración y reclasificación. La calibración se evaluó mediante gráficos de calibración, intercepto y pendiente de calibración y la puntuación de Brier, con la prueba de Hosmer-Lemeshow informada como verificación complementaria de calibración. La reclasificación se cuantificó mediante la mejora continua neta de reclasificación y la mejora integrada de discriminación. Los análisis de sensibilidad preespecificados incluyeron reducir la ventana de muestreo del día 2 de 36–60 h a 42–54 h, repetir los análisis con una definición más estricta de progresión de disfunción orgánica (ΔpSOFA ≥ 3), excluir valores extremos de biomarcadores clínicamente atípicos tras la verificación de la fuente y repetir los modelos principales tras excluir a los pacientes sin mediciones dinámicas completas de HBP. Los análisis de sensibilidad se informaron en la sección de Resultados o en tablas complementarias, según correspondiera.
Los análisis estadísticos se realizaron en R versión 4.3.2 utilizando un entorno de análisis bloqueado (versión del entorno de análisis v1.0). Los paquetes principales fueron pROC versión 1.18.5 para los análisis ROC, survival versión 3.5-7 para los modelos de tiempo hasta el evento y rms versión 6.8-0 para los análisis de calibración. La mejora continua neta de la reclasificación y la mejora integrada de la discriminación se calcularon utilizando rutinas validadas implementadas en el conjunto prescrito de guiones de análisis (conjunto de paquetes de guiones v1.0). Todas las pruebas fueron bilaterales y se consideró estadísticamente significativo un valor de p < 0,05.