Cette étude a été approuvée par le comité d'éthique institutionnel de l'Hôpital populaire de la région autonome zhuang de Guangxi (numéro d'approbation IRB : KY-ZC-2023-035). L'exigence de consentement éclairé a été levée, car cette analyse rétrospective utilisait des données cliniques anonymisées recueillies dans le cadre des soins courants.
Conception de l'étude et cadre d'étude
Il s'agissait d'une étude de cohorte observationnelle rétrospective monocentrique menée en unité de soins intensifs pédiatriques (USIP) au sein de l'Hôpital populaire de la région autonome zhuang du Guangxi. Les admissions consécutives en USIP entre le 1er janvier 2022 et le 31 décembre 2025 ont été examinées. Les données cliniques et de laboratoire ont été extraites du système de dossier médical électronique, du système d'information de laboratoire (LIS) et du système d'information de soins critiques.
Le moment de l'admission en USI pédiatrique a été défini comme le temps zéro pour toutes les fenêtres d'exposition, de biomarqueurs, de soutien organique et de résultats. Avant l'extraction des données, les critères d'éligibilité, les fenêtres de prélèvement des biomarqueurs, les définitions des critères de jugement, les indicateurs de gravité et le plan d'analyse statistique ont été prédéfinis et mis en œuvre selon un protocole d'extraction fixe (spécification d'extraction LIS version 1.0 ; dictionnaire de données cliniques version 1.0). Après extraction, tous les dossiers ont été rendus anonymes. Deux investigateurs ont vérifié indépendamment les horodatages et la cohérence des sources pour les valeurs des biomarqueurs, les variables liées à l'évaluation pédiatrique de l'insuffisance organique séquentielle (pSOFA), les enregistrements de traitement vasactif et les événements de soutien organique. Les divergences ont été résolues par examen conjoint des documents sources. L'assemblage de la cohorte a suivi des étapes prédéfinies, incluant le dépistage, la confirmation de l'éligibilité, la confirmation de la fenêtre de biomarqueur et l'examen de l'exhaustivité des critères de jugement. Le flux final de la cohorte est présenté dans Figure 1.
Participants et critères d'éligibilité
Les enfants âgés de 1 mois à 14 ans admis en USI pédiatrique (PICU) pour un état de sepsis étaient éligibles. Les cas candidats ont été identifiés à partir des dossiers d'admission en PICU et des diagnostics liés à la sepsis, puis confirmés par une analyse manuelle des dossiers. La sepsis était définie comme une infection suspectée ou confirmée accompagnée d'un dysfonctionnement organique documenté pendant les soins en PICU et confirmé par une évaluation basée sur le pSOFA.
Afin de garantir une exposition comparable aux biomarqueurs précoces, les patients inclus devaient avoir fait l'objet d'une mesure de la protéine liant l'héparine (HBP) et de la protéine C-réactive (CRP) dans les 24 heures suivant l'admission en USPIC, ainsi que disposer de données suffisantes pour le calcul du pSOFA. Les critères d'inclusion étaient les suivants : âge compris entre 1 mois et 14 ans, admission en USPIC pour sepsis, dosage de l'HBP et de la CRP disponible entre 0 et 24 heures, et données complètes permettant l'évaluation de l'issue principale, incluant le pSOFA de base, l'évaluation du pSOFA à 72 heures et le suivi de survie. Les critères d'exclusion étaient l'absence de HBP ou de CRP au moment de l'admission, un séjour en USPIC inférieur à 24 heures, ou des données manquantes pour l'un des éléments requis pour l'évaluation de l'issue principale. <24 h, absence de variables clés nécessaires à l'adjudication des critères de jugement, ou immunosuppression sévère susceptible d'altérer de manière importante les profils des biomarqueurs inflammatoires. Une immunosuppression sévère était définie par une hémopathie maligne active, la réception de chimiothérapie ou d'un traitement immunosuppresseur lié à une greffe, un traitement prolongé par corticoïdes systémiques à dose immunosuppressive, une immunodéficience primaire, ou tout autre traitement immunosuppresseur intensif documenté avant ou pendant l'admission en USIP. En cas d'admissions répétées en USIP, seule la première admission éligible était incluse.
Définitions des critères d'évaluation et des métriques de gravité
Le critère d'évaluation principal était la progression précoce d'un dysfonctionnement organique dans les 72 heures suivant l'admission en USPIC, définie par une augmentation d'au moins 2 points du score pSOFA par rapport à la valeur initiale, atteinte au cours des 72 premières heures (ΔpSOFA ≥ 2). Le score pSOFA de base correspondait au premier score pSOFA enregistré entre 0 et 24 heures après l'admission en USPIC ; lorsque plusieurs valeurs étaient disponibles, le score le plus proche de l'heure d'admission était retenu. La valeur pSOFA à 72 heures était définie comme le score pSOFA maximal enregistré entre 0 et 72 heures.
Une gravité élevée de la maladie a été définie a priori comme un pSOFA maximal sur 72 h ≥ 8, et ce seuil a été appliqué uniformément dans les analyses catégorielles. Ce seuil avait été prédéfini en se basant sur la pratique de stratification de la gravité de l'unité et sur son interprétabilité clinique pour l'identification précoce du sepsis pédiatrique à haut risque. Le pSOFA a également été analysé comme une mesure continue de la charge de dysfonctionnement organique.
Les critères d'évaluation secondaires étaient le choc septique, la ventilation mécanique invasive, l'utilisation de médicaments vasoactifs, la thérapie de remplacement rénal continu (CRRT), la mortalité en USI pédiatrique et la mortalité à 28 jours pour toutes causes. Le choc septique était défini comme un échec circulatoire persistant nécessitant un soutien vasoactif après une réanimation liquidienne, avec une hypoperfusion tissulaire documentée, incluant une élévation de la lactatémie et/ou des anomalies des signes de perfusion périphérique. Les critères d'évaluation relatifs au soutien organique ont été codés comme des variables binaires au sein de fenêtres précises prédéfinies. Afin d'obtenir une mesure plus fine du soutien hémodynamique, le score vasoactif-inotrope (VIS) a également été extrait ou calculé pour les patients recevant un traitement vasoactif durant les 24 premières heures suivant l'admission en USI pédiatrique. Le VIS a été utilisé comme variable quantitative d'intensité du soutien organique dans les analyses de corrélation et d'association avec la gravité, tandis que l'utilisation binaire de la thérapie vasoactive a été conservée comme critère d'évaluation cliniquement interprétable reflétant la nécessité d'un traitement.
Mesure des biomarqueurs et collecte des données
L'HBP et la CRP ont été mesurés au laboratoire central à l'aide de protocoles automatisés normalisés d'immunodosage, conformément aux procédures habituelles de contrôle qualité clinique. Les analyses des biomarqueurs ont été réalisées sur des échantillons de sang veineux prélevés dans le cadre des soins de routine en USPIC et traités selon le protocole préanalytique standard du laboratoire. Les résultats d'HBP ont été exprimés en ng/mL et ceux de CRP en mg/L. L'HBP a été dosé à l'aide d'un test immunoturbidimétrique amplifié par particules, sur une plateforme entièrement automatisée d'analyses cliniques et d'immunodosages. La CRP a été mesurée à l'aide d'un test immunoturbidimétrique à haute sensibilité. Le protocole analytique et les procédures internes de contrôle qualité ont suivi un protocole de laboratoire soumis à versionnement (version 2022.1 du protocole opératoire standard du laboratoire) et sont restés inchangés dans leur logique analytique fondamentale pendant toute la durée de l'étude.
Les dossiers de laboratoire ont été extraits du système d'information de laboratoire (SIL) à l'aide d'une requête verrouillée et d'un processus de mappage des champs (spécification d'extraction SIL version 1.0), incluant à la fois l'heure de prélèvement sanguin et l'heure de rapport. Les valeurs initiales de la PNH et de la CRP ont été définies comme les premières valeurs rapportées entre 0 et 24 heures suivant l'admission en USPIC. Les valeurs de PNH et de CRP du jour 2 ont été définies comme les mesures les plus proches de 48 heures, prises dans une fenêtre fixe de 36 à 60 heures. Cette fenêtre du jour 2 a été utilisée de manière constante dans toutes les analyses dynamiques.
La variation absolue de HBP a été calculée comme la valeur de HBP au jour 2 moins la valeur de base de HBP. La variation relative de HBP (%) a été calculée selon la formule suivante : [(valeur de HBP au jour 2 − valeur de base de HBP) / valeur de base de HBP] × 100. Si la valeur de base de HBP était inférieure à la limite de quantification du test ou enregistrée comme nulle après la conversion du système, la variation relative de HBP n’a pas été calculée ; ces patients ont été conservés dans les analyses utilisant la valeur de base de HBP et la variation absolue de HBP, mais exclus des analyses basées sur le pourcentage de variation. Le rapport HBP sur albumine a été calculé en divisant la concentration de HBP par celle de l’albumine sérique mesurée dans l’intervalle de temps correspondant, et a été évalué comme un indice composite exploratoire reflétant l’activation endothéliale et la perte protéique liée aux fuites capillaires.
Outre la protéine de liaison à l'heparine (HBP) et la protéine C-réactive (CRP), des marqueurs inflammatoires, liés à la perfusion et à la fonction organique habituellement disponibles ont été extraits à des fins de comparaison et de modélisation multivariée. Ces marqueurs comprenaient le nombre de globules blancs, la procalcitonine, le lactate, l'albumine, la créatinine, le nombre de plaquettes, les paramètres de coagulation et les paramètres de fonction hépatique, lorsque disponibles. Les concentrations d'interleukine-6 (IL-6) et d'interleukine-8 (IL-8) ont été extraites du système d'information de laboratoire (LIS) lorsqu'elles avaient été mesurées dans le cadre de l'évaluation clinique de routine au cours des mêmes fenêtres temporelles précoces. La mesure de l'IL-6 et de l'IL-8 n'était pas requise pour l'inclusion dans la cohorte et ces paramètres n'ont été utilisés que dans les analyses où les données étaient disponibles. Cette règle a été appliquée afin d'éviter d'exclure des patients autrement éligibles uniquement parce que le dosage des cytokines n'était pas systématiquement réalisé dans tous les cas.
Les données démographiques, les comorbidités, la source de l'infection, les résultats microbiologiques lorsqu'ils étaient disponibles, les signes vitaux au moment de l'admission en USPIC, le traitement antimicrobien précoce, les mesures de réanimation et les interventions de soutien organique ont été extraites du dossier médical électronique et du système d'information de soins critiques à l'aide du dictionnaire clinique de données prédéfini (version 1.0). Les variables relatives au soutien organique comprenaient la ventilation mécanique invasive, le traitement vasactif, l'HTRC, le volume des bolus de fluides et l'indice VIS durant les premières 24 heures, le cas échéant. Les variables nécessaires au pSOFA ont été recueillies dans des fenêtres prédéfinies. Si une composante du pSOFA générée automatiquement était incohérente avec les valeurs sources, deux investigateurs examinaient le dossier et recalculaient le score en utilisant les règles de notation prédéfinies. Le contrôle de qualité des données comprenait la validation des horodatages, le dépistage des doublons, les vérifications de plages et de distribution, la confirmation des valeurs extrêmes à partir des sources et la journalisation des corrections manuelles.
Afin de préserver la reproductibilité tout en évitant un libellé spécifique à une marque dans le texte principal du protocole, la documentation des tests a été enregistrée selon la classe de méthode, le type de plateforme, l'environnement logiciel, les règles de plage de résultats rapportables, l'utilisation des étalons et des contrôles, ainsi que les procédures internes de contrôle qualité. Des informations détaillées sur les réactifs, les étalons, les contrôles et les logiciels sont fournies dans le tableau des matériaux. Le contrôle qualité interne a été effectué quotidiennement, et seules les séries d'analyses satisfaisant aux critères d'acceptation du laboratoire concernant le coefficient de variation ont été retenues. Les valeurs inférieures à la limite de quantification inférieure ont été fixées à la limite inférieure rapportable. Les valeurs supérieures à la plage supérieure rapportable ont été fixées à la limite supérieure rapportable, sauf si un résultat répété validé provenant du même échantillon était disponible. Les unités de laboratoire ont été vérifiées pour assurer leur cohérence avant l'analyse ; aucune conversion d'unité n'a été effectuée, sauf si une inadéquation d'unité du système d'information de laboratoire (LIS) a été confirmée par comparaison avec le rapport d'origine.
Données manquantes et règles de traitement des données
Les règles de traitement des données ont été prédéfinies. Les patients présentant des données manquantes pour la PTH de base, la CRP de base, les variables de pSOFA de base ou les variables de pSOFA à 72 heures ont été exclus des analyses des critères d'évaluation principaux. Les patients dont la PTH du jour 2 était manquante ont été conservés dans les analyses des biomarqueurs de base et dans les résumés globaux des résultats cliniques, mais exclus des analyses nécessitant des indices dynamiques de PTH. Les patients sans mesures d'IL-6 ou d'IL-8 ont été conservés dans toutes les analyses principales et exclus uniquement des comparaisons des modèles incluant des cytokines.
Une analyse des cas complets a été utilisée pour chaque modèle multivarié, et la taille de l'échantillon propre à chaque analyse a été indiquée. L'imputation multiple n'a pas été utilisée dans les analyses principales car l'étude se concentrait sur les trajectoires temporelles des biomarqueurs observés et du dysfonctionnement organique. Le risque de biais de sélection dû à l'exclusion des patients ne disposant pas de données de HBP de base, de CRP de base ou d'évaluation complète des critères pSOFA a été pris en compte lors de l'interprétation des résultats.
Les valeurs extrêmes obtenues en laboratoire n'ont pas été exclues uniquement pour des raisons statistiques. Les valeurs aberrantes ont été vérifiées par rapport aux dossiers d'origine et conservées si elles étaient confirmées comme étant des valeurs réelles. Seules les erreurs manifestes de documentation ou de saisie d'unités ont été corrigées, et toutes les corrections ont été enregistrées dans le journal d'audit avec des identifiants traçables.
Analyse statistique
Toutes les analyses ont suivi un plan d'analyse statistique prédéfini (version 1.0 de la spécification d'analyse). Les variables continues ont été évaluées par inspection des histogrammes et par le test de Shapiro-Wilk. Les variables normalement distribuées ont été résumées par la moyenne ± l'écart-type et comparées à l'aide de tests t pour échantillons indépendants. Les variables non normalement distribuées ont été résumées par la médiane (intervalle interquartile) et comparées à l'aide du test de Mann-Whitney U. Les variables catégorielles ont été résumées par des effectifs (pourcentages) et comparées à l'aide du test du chi-deux ou du test exact de Fisher.
Pour les analyses descriptives, les patients ont été stratifiés selon les quartiles de la pression artérielle diastolique (HBP) de base et séparément selon la présence ou l'absence d'une progression précoce du dysfonctionnement organique (ΔpSOFA ≥ 2). Les quartiles de l'HBP de base ont été utilisés pour décrire les profils dose-réponse à travers la distribution observée du biomarqueur, sans supposer de seuil clinique universel, car les seuils pédiatriques validés de l'HBP pour la gravité du sepsis restent incertains et peuvent varier selon la plateforme d'analyse, la répartition par âge et la fenêtre d'échantillonnage. L'âge a été analysé comme une covariable continue dans les modèles principaux. Afin de tenir compte de l'hétérogénéité clinique liée à l'âge, des analyses préspécifiées par sous-groupe ont également été réalisées en utilisant des strates d'âge de <1 an et ≥1 an, des vérifications exploratoires supplémentaires par catégorie d'âge ayant été effectuées lorsque la taille de l'échantillon le permettait.
Des coefficients de corrélation de Spearman ont été utilisés pour évaluer les associations entre les paramètres de la PTH, notamment la PTH initiale, la PTH du jour 2, la variation absolue de la PTH, la variation relative de la PTH et le rapport PTH/albumine, et les indicateurs de gravité clinique. Les variables reflétant une prise en charge précoce des organes, notamment la ventilation mécanique invasive, le traitement vasoactif, l'épuration extrarénale continue (CRRT), le volume de bolus de liquide, le nombre de jours sous ventilateur, la durée de séjour en USPIC et l'indice VIS, ont été considérées comme des corrélats descriptifs et analytiques de la gravité clinique et de la prise en charge.
La performance discriminante pour la progression précoce du dysfonctionnement organique et la gravité élevée de la maladie a été évaluée à l'aide de courbes caractéristiques de fonctionnement du récepteur (ROC). La surface sous la courbe ROC (AUC) et les intervalles de confiance à 95 % ont été calculés pour le HBP, la CRP, la procalcitonine, le lactate, l'albumine, le pSOFA, le rapport HBP sur albumine, ainsi que pour les modèles combinés. Les AUC ont été comparées à l'aide du test de DeLong. Les valeurs seuils optimales ont été déterminées par l'indice de Youden, et la sensibilité, la spécificité, la valeur prédictive positive, la valeur prédictive négative et les rapports de vraisemblance ont été rapportés.
Une régression logistique multivariée a été utilisée pour estimer les associations indépendantes entre les paramètres de la pression artérielle élevée (HBP) et les résultats binaires de gravité, notamment la progression précoce du dysfonctionnement organique et la gravité élevée de la maladie. Les covariables ont été sélectionnées a priori en fonction de leur pertinence clinique et de la disponibilité des données, et comprenaient l'âge, la catégorie de source d'infection, le score pSOFA de base, le taux de lactate, l'albumine, la procalcitonine, la protéine C-réactive (CRP), la présence d'un choc septique à l'admission et, le cas échéant, les variables de prise en charge organique précoce. Les modèles de prédiction principaux pour les résultats de gravité précoce ont été limités aux variables disponibles dans la fenêtre d'évaluation précoce prédéfinie afin de réduire le chevauchement temporel avec l'évolution du résultat. La valeur pronostique ajoutée a été évaluée en ajoutant progressivement les variables HBP au modèle clinique de base. Pour la mortalité toutes causes confondues à 28 jours, une régression de Cox aux risques proportionnels a été utilisée, et l'hypothèse de proportionnalité des risques a été vérifiée à l'aide des résidus de Schoenfeld. Avant l'ajustement du modèle, la colinéarité entre les prédicteurs candidats a été évaluée, et le nombre de covariables inclus dans chaque modèle a été limité afin de préserver un ratio adéquat événements-par-variable.
L'amélioration du modèle après l'ajout de l'HBP a été évaluée à l'aide de métriques de discrimination, de calibration et de reclassement. La calibration a été évaluée à l'aide de graphiques de calibration, de l'ordonnée à l'origine et de la pente de calibration, ainsi que du score de Brier, le test de Hosmer-Lemeshow étant indiqué comme vérification complémentaire de la calibration. Le reclassement a été quantifié à l'aide de l'amélioration continue du reclassement net et de l'amélioration intégrée de la discrimination. Des analyses de sensibilité préspécifiées comprenaient le resserrement de la fenêtre d'échantillonnage du jour 2, passant de 36 à 60 heures à 42 à 54 heures, la répétition des analyses avec une définition plus stricte de la progression du dysfonctionnement organique (ΔpSOFA ≥ 3), l'exclusion des valeurs extrêmes atypiques des biomarqueurs après vérification initiale, et la répétition des modèles principaux après exclusion des patients ne disposant pas de mesures dynamiques complètes de l'HBP. Les analyses de sensibilité ont été rapportées dans la section Résultats ou dans des tableaux supplémentaires, selon le cas.
Les analyses statistiques ont été réalisées dans R version 4.3.2 à l'aide d'un environnement d'analyse verrouillé (version v1.0 de l'environnement d'analyse). Les principaux packages utilisés étaient pROC version 1.18.5 pour les analyses ROC, survival version 3.5-7 pour les modèles de temps jusqu'à l'événement, et rms version 6.8-0 pour les analyses de calibration. L'amélioration continue de la reclassification nette et l'amélioration intégrée de la discrimination ont été calculées à l'aide de routines validées implémentées dans le jeu de scripts d'analyse prédéfini (jeu de scripts v1.0). Tous les tests étaient bilatéraux, et une valeur de p < 0,05 a été considérée comme statistiquement significative.