Le protocole de l'étude a été examiné et approuvé par le Comité d'éthique de l'Hôpital de médecine traditionnelle chinoise du district de Santai, dans la province du Sichuan (n° d'approbation : 2025003). L'étude a été menée conformément aux principes de la Déclaration d'Helsinki. Un consentement éclairé écrit a été obtenu auprès de chaque patient ou d'un proche parent avant leur inclusion dans l'étude. Cette étude de cohorte rétrospective a passé en revue les dossiers cliniques de patients âgés hospitalisés pour une exacerbation aiguë de la maladie pulmonaire obstructive chronique (AECOPD) dans le service de pneumologie et de soins critiques de l'Hôpital de médecine traditionnelle chinoise du district de Santai entre mai 2023 et mai 2025. L'objectif de l'étude était d'identifier des prédicteurs cliniques couramment disponibles d'hypertension pulmonaire (PH) et de développer un nomogramme permettant une estimation personnalisée du risque de PH. Les instruments de laboratoire, le système d'échographie, les fournitures pour prélèvement sanguin, les outils d'extraction des données et les logiciels statistiques utilisés dans cette procédure sont indiqués dans le Tableau des matériaux.
1. Dépistage des patients et évaluation des critères d'éligibilité
Les patients hospitalisés ayant reçu un diagnostic de rechute aiguë de la BPCO (AECOPD) lors de la période d'étude ont été sélectionnés via le système informatisé des dossiers médicaux des patients hospitalisés. Chaque dossier a été examiné afin de confirmer que le patient avait au moins 60 ans et que le diagnostic de broncho-pneumopathie chronique obstructive était conforme aux critères diagnostiques de l'Initiative mondiale pour la maladie pulmonaire obstructive chronique (GOLD) de 202115. Une rechute aiguë de la BPCO était définie comme une aggravation aiguë des symptômes respiratoires nécessitant une hospitalisation et un traitement complémentaire. Pour les patients ayant plus d'une admission pour rechute aiguë de la BPCO durant la période d'étude, seule la première hospitalisation éligible a été incluse.
Les patients ont été inclus lorsqu'ils répondaient à tous les critères suivants : âge d'au moins 60 ans ; hospitalisation pour une exacerbation aiguë de la BPCO entre mai 2023 et mai 2025 ; échocardiographie transthoracique disponible durant l'hospitalisation initiale ; prélèvement sanguin veineux au premier jeûne disponible après l'admission ; et informations cliniques, biologiques et échocardiographiques complètes nécessaires à la construction du modèle. Les patients ont été exclus s'ils présentaient une insuffisance hépatique sévère, une maladie rénale au stade terminal, une malignité active, une hypertension pulmonaire attribuable à d'autres causes connues, une maladie cardiaque gauche primaire, une maladie du tissu conjonctif telle que la sclérodermie systémique, une tuberculose pulmonaire active, ou des variables essentielles manquantes requises pour la classification du résultat ou la construction du modèle.
Avant que les dossiers ne soient supprimés en raison d'informations incomplètes, l'âge, le sexe, les antécédents de tabagisme, les antécédents de consommation d'alcool, les antécédents d'allergie médicamenteuse, l'hypertension, le diabète sucré, la maladie coronarienne, l'asthme, la bronchiectasie, l'infection pulmonaire, l'emphysème, l'insuffisance respiratoire, l'insuffisance cardiaque, le statut échocardiographique de PH, les paramètres de la numération formule sanguine, l'albumine, la créatinine, la protéine C-réactive de haute sensibilité, le D-dimère, le fibrinogène, le peptide natriurétique de type B et le ratio neutrophiles-lymphocytes ont été vérifiés. Le nombre de dossiers dépistés et exclus, les raisons des exclusions ainsi que les cas finalement inclus ont été enregistrés afin d'établir le diagramme de flux de sélection des patients (Figure 1).

Figure 1 : Sélection des patients et constitution de la cohorte. Le diagramme de flux illustre l'identification des dossiers, l'évaluation des critères d'éligibilité, les exclusions et le regroupement final des patients âgés hospitalisés pour une BPCO-GEA. Les patients ont été classés dans le groupe sans HP ou dans le groupe HP selon les résultats de l'échocardiographie transthoracique. Abréviations : BPCO-GEA, exacerbation aiguë de la maladie pulmonaire obstructive chronique ; HP, hypertension pulmonaire. Veuillez cliquer ici pour afficher une version agrandie de cette figure.
2. Regroupement clinique et définition diagnostique de l'hypertension pulmonaire
Les patients inclus ont été classés dans un groupe sans HP ou dans un groupe HP selon les résultats de l'échocardiographie transthoracique réalisée pendant l'hospitalisation. Le cadre diagnostique était fondé sur les recommandations de 2015 de la Société européenne de cardiologie et de la Société européenne de pneumologie concernant l'hypertension pulmonaire16. L'échocardiographie transthoracique au repos a été utilisée comme base diagnostique dans ce jeu de données rétrospectif. Les examens ont été réalisés par des médecins échocardiographes formés, en utilisant des vues parasternales, apicales et sous-costales standard. L'interprétabilité échocardiographique a été revue avant l'affectation aux groupes, car des fenêtres acoustiques médiocres sont fréquentes chez les patients atteints de bronchopneumopathie chronique obstructive. Un dossier était considéré comme approprié pour la classification de l'HP lorsque le signal Doppler du reflux tricuspidien était mesurable ou lorsque la taille du ventricule droit, l'épaisseur de la paroi et la fonction systolique pouvaient être évaluées à partir de vues cardiaques interprétables.
La pression artérielle pulmonaire systolique a été estimée à l'aide de la vélocité maximale du reflux tricuspidien et de la pression auriculaire droite estimée. La pression auriculaire droite a été évaluée à partir du diamètre de la veine cave inférieure et de son collapsus inspiratoire lorsque les deux mesures étaient disponibles. Une hypertension pulmonaire (HP) a été définie par une pression artérielle pulmonaire systolique estimée d'au moins 40 mmHg. Lorsque la pression artérielle pulmonaire systolique ne pouvait pas être estimée de façon fiable, le diagnostic d'HP nécessitait la présence d'indices échocardiographiques indirects d'atteinte du ventricule droit, notamment un agrandissement du ventricule droit, une hypertrophie du ventricule droit, une dysfonction systolique du ventricule droit, un aplatissage du septum interventriculaire ou un élargissement de l'artère pulmonaire. Les dossiers échocardiographiques comportant des images floues, des signaux Doppler de reflux tricuspidien non traçables ou des anomalies incohérentes du cœur droit ont été revus par un deuxième médecin échocardiographe avant l'attribution définitive au groupe. Les dossiers restant ininterprétables après cette révision ont été exclus de l'analyse finale. Dans la cohorte finale, 230 patients âgés répondaient aux critères d'éligibilité ; 120 patients ont été affectés au groupe sans HP et 110 patients au groupe HP.
3. Extraction des variables démographiques et cliniques
Toutes les données cliniques ont été extraites du système informatisé des dossiers médicaux des patients hospitalisés à l'aide d'un formulaire standardisé de collecte des données. Ont été enregistrés : le sexe, l'âge, les antécédents de tabagisme, les antécédents de consommation d'alcool, les antécédents d'allergie médicamenteuse, l'hypertension artérielle, le diabète sucré, la maladie coronarienne, l'asthme, la bronchiectasie, l'infection pulmonaire, l'emphysème, l'insuffisance respiratoire et l'insuffisance cardiaque. Un antécédent de tabagisme était défini comme une consommation cumulée d'au moins 100 cigarettes au cours de la vie, et un antécédent de consommation d'alcool, comme une consommation régulière d'alcool au moins deux fois par semaine.
L'infection pulmonaire a été diagnostiquée en intégrant les symptômes respiratoires, les signes physiques, les résultats d'imagerie thoracique et, lorsque disponibles, les éléments microbiologiques. L'insuffisance cardiaque a été diagnostiquée selon les manifestations cliniques, les signes physiques, les taux de peptide natriurétique de type B et les résultats échocardiographiques. Les informations sur les comorbidités ont été vérifiées croisées avec les diagnostics de sortie, les notes d'évolution, les comptes rendus d'imagerie, les rapports de laboratoire et les dossiers de consultation, lorsque disponibles.
4. Collecte et traitement des variables de laboratoire
L'échantillon sanguin veineux à jeun prélevé le matin suivant l'admission à l'hôpital a été utilisé comme source de données de laboratoire. Ce moment correspondait au premier prélèvement sanguin matinal après l'admission et avant tout ajustement habituel du traitement diurne, lorsque ces informations étaient disponibles dans le dossier médical. Les paramètres de la numération formule sanguine, les paramètres biochimiques, les marqueurs de la coagulation, les marqueurs inflammatoires et les marqueurs de stress cardiaque ont été extraits du système d'information laboratoriale de l'hôpital.
Le nombre de neutrophiles, le nombre de lymphocytes, le nombre de plaquettes, le volume moyen des plaquettes et le coefficient de variation de la largeur de distribution des hématies ont été enregistrés à partir du compte rendu d'hémogramme. Le ratio neutrophiles/lymphocytes a été calculé comme suit :

L'albumine, la créatinine, le D-dimère, le fibrinogène, la protéine C-réactive à haute sensibilité et le peptide natriurétique de type B ont été relevés à partir des rapports de laboratoire correspondants. Des unités de mesure cohérentes ont été utilisées pour l'ensemble du jeu de données : l'albumine en g/L, la créatinine en µmol/L, la protéine C-réactive à haute sensibilité en mg/L, le D-dimère en µg/mL, le fibrinogène en g/L, le nombre de neutrophiles, le nombre de lymphocytes et le nombre de plaquettes en ×109/L, le volume plaquettaire moyen en fL, et le peptide natriurétique de type B en pg/mL.
5. Nettoyage et préparation des données
Le jeu de données a été examiné avant l'analyse statistique. Chaque enregistrement inclus a été vérifié quant à l'éligibilité, l'affectation au groupe, l'exhaustivité du laboratoire, l'interprétabilité échocardiographique et la disponibilité des prédicteurs candidats. Les identifiants personnels directs ont été supprimés et chaque patient s'est vu attribuer un numéro d'identification d'étude. Les enregistrements en double, les hospitalisations répétées, les étiquettes de groupe incohérentes, les dates impossibles et les valeurs de laboratoire invraisemblables ont été vérifiées par rapport au dossier médical électronique d'origine et au rapport de laboratoire. Aucune imputation n'a été effectuée pour les variables clés d'issue ou de prédiction, car le modèle final était basé sur une analyse des cas complets.
Les variables continues ont été examinées pour leurs caractéristiques de distribution. Les variables normalement distribuées ont été résumées par la moyenne ± l'écart type ; les variables non normalement distribuées par la médiane avec l'intervalle interquartile ; et les variables catégorielles par les effectifs et les pourcentages. Avant la modélisation multivariée, les variables cliniquement apparentées ont été vérifiées pour détecter d'éventuelles redondances. Le nombre de neutrophiles et celui de lymphocytes ont été examinés conjointement avec le ratio neutrophiles/lymphocytes ; l'infection pulmonaire a été examinée conjointement avec le nombre de neutrophiles, la protéine C-réactive de haute sensibilité et le ratio neutrophiles/lymphocytes ; l'insuffisance cardiaque a été examinée conjointement avec le peptide natriurétique de type B et les résultats échocardiographiques cardiaques ; le D-dimère et le fibrinogène ont été examinés comme variables liées à la coagulation ; et le nombre de plaquettes, le volume moyen des plaquettes et le coefficient de variation de la largeur de distribution des globules rouges ont été examinés comme variables hématologiques.
6. Comparaison univariée entre les groupes
Les caractéristiques démographiques, les comorbidités, les indicateurs de laboratoire et les variables de regroupement échocardiographiques ont été comparés entre les groupes sans HP et avec HP. Un test t pour échantillons indépendants a été utilisé pour les variables continues normalement distribuées. Le test U de Mann-Whitney a été utilisé pour les variables continues non normalement distribuées. Le test χ2 ou le test exact de Fisher a été utilisé pour les variables catégorielles, selon les effectifs théoriques par cellule. Une valeur P bilatérale inférieure à 0,05 a été considérée comme statistiquement significative. Les variables présentant des différences significatives entre les groupes ont été prises en compte pour une évaluation multivariée, conjointement avec les prédicteurs cliniquement pertinents.
7. Analyse de régression logistique multivariée
Une régression logistique binaire a été utilisée pour identifier les prédicteurs indépendants associés à l'hypertension pulmonaire (PH) chez les patients âgés atteints d'une exacerbation aiguë de la BPCO (AECOPD). Le statut d'hypertension pulmonaire a été utilisé comme variable dépendante. Les prédicteurs continus ont été saisis dans leurs unités de mesure d'origine, et les prédicteurs catégoriels ont été codés sous forme de variables binaires. Les variables candidates ont été sélectionnées en fonction des résultats univariés et de leur pertinence clinique. Les variables présentant une valeur P < 0,05 dans la comparaison univariée ont d'abord été identifiées, puis examinées pour vérifier leur redondance, leur interprétabilité et leurs chevauchements cliniques.
Les variables candidates initiales comprenaient la protéine C-réactive de haute sensibilité, le D-dimère, le coefficient de variation de la largeur de distribution des globules rouges, le nombre de plaquettes, le volume plaquettaire moyen, le ratio neutrophiles-lymphocytes, le peptide natriurétique de type B, la maladie coronarienne et les antécédents d'asthme. Le nombre de neutrophiles et le nombre de lymphocytes n'ont pas été saisis conjointement avec le ratio neutrophiles-lymphocytes, car ils constituaient les composantes de ce ratio. L'infection pulmonaire et l'insuffisance cardiaque ont été évaluées lors de l'évaluation de la redondance. Afin d'examiner si leur exclusion influençait l'interprétation du modèle, une analyse de sensibilité a été réalisée en ajoutant l'infection pulmonaire et l'insuffisance cardiaque au jeu de prédicteurs principaux. Les estimations de régression, la discrimination, l'étalonnage et le bénéfice net clinique ont été comparés entre les modèles principal et de sensibilité.
La multicolinéarité entre les prédicteurs candidats a été évaluée à l'aide des facteurs d'inflation de la variance. Les variables catégorielles ont été codées avant le calcul, et les variables continues ont été saisies dans les mêmes unités que celles utilisées pour la modélisation par régression. Un facteur d'inflation de la variance inférieur à 5 a été considéré comme indiquant l'absence de multicolinéarité sévère. Les coefficients de régression, les erreurs types, les valeurs de χ2 de Wald, les rapports de cotes, les intervalles de confiance à 95 %, les valeurs P et les facteurs d'inflation de la variance ont été rapportés pour les variables incluses dans le modèle multivarié.
8. Construction du nomogramme
Le nomogramme a été construit à partir des prédicteurs indépendants retenus dans le modèle de régression logistique multivariée. Les prédicteurs finaux étaient le D-dimère, le ratio neutrophiles-lymphocytes, le peptide natriurétique de type B et les antécédents d'asthme. Le modèle de régression logistique a été converti en un nomogramme basé sur un système de points à l'aide du logiciel R et du package rms. Un objet de distribution des données a été généré pour l'ensemble de données de modélisation, le modèle de régression logistique a été ajusté, et les coefficients de régression ont été utilisés pour attribuer des points à chaque prédicteur. Pour chaque patient, les points correspondant au D-dimère, au ratio neutrophiles-lymphocytes, au peptide natriurétique de type B et aux antécédents d'asthme ont été additionnés afin d'obtenir un score total, qui a ensuite été associé à la probabilité estimée d'hypertension pulmonaire (PH).
Le nomogramme a été interprété à l'aide de procédures standard fondées sur des points, décrites dans des études antérieures sur les modèles de prédiction clinique17,18. La valeur de D-dimère du patient, le ratio neutrophiles/lymphocytes, la valeur du peptide natriurétique de type B et le statut asthmatique ont été repérés sur leurs axes respectifs. Une ligne verticale a été tracée depuis chaque valeur prédictive jusqu'à l'échelle de points. Les points ont été additionnés pour obtenir le score total, puis une ligne verticale a été tracée depuis l'échelle du score total jusqu'à l'axe de la probabilité prédite, afin d'obtenir la probabilité estimée individualisée de PH. Le patient à titre d'illustration utilisé dans le manuscrit a été choisi parmi des valeurs situées dans ou proches de la distribution clinique observée dans la cohorte.
9. Évaluation des performances du modèle
La discrimination du modèle a été évaluée à l'aide de l'analyse de la courbe ROC (caractéristique de fonctionnement du récepteur) et de l'aire sous la courbe. L'aire sous la courbe a été indiquée avec un intervalle de confiance à 95 %. Le nomogramme a été comparé aux prédicteurs individuels ainsi qu'à un modèle combinant des marqueurs de laboratoire quantitatifs. La calibration a été évaluée en comparant les probabilités prédites aux résultats observés de PH. Une courbe de calibration a été générée à l'aide d'un rééchantillonnage par bootstrap. L'intercept de calibration, la pente de calibration, l'erreur absolue moyenne, le score de Brier et le test de Hosmer-Lemeshow ont été rapportés. L'intercept de calibration corrigé par bootstrap et la pente de calibration corrigée par bootstrap ont également été calculés. Le test de Hosmer-Lemeshow a été interprété conjointement avec la courbe de calibration et les indices quantitatifs de calibration.
L'utilité clinique a été évaluée à l'aide d'une analyse de courbe décisionnelle. Le bénéfice net a été tracé en fonction des probabilités seuils et comparé aux stratégies de référence « traiter tous » et « ne traiter aucun ». La plage de probabilité seuil dans laquelle le modèle offrait un bénéfice net supérieur à celui des stratégies de référence a été identifiée.
10. Validation interne et reproductibilité
Une validation interne a été réalisée à l'aide d'un rééchantillonnage bootstrap avec 1 000 répétitions. Pour chaque répétition bootstrap, le modèle a été réajusté et testé afin d'estimer le biais optimiste de la performance du modèle. Les indices de discrimination et de calibration corrigés pour le biais optimiste ont été calculés et rapportés conjointement avec la performance apparente du modèle.
Les analyses statistiques ont été réalisées à l'aide de logiciels statistiques standards et de la version 4.5.2 de R. Des packages R ont été utilisés pour la modélisation par régression logistique, la construction de nomogrammes, l'analyse des courbes ROC (caractéristique de fonctionnement du récepteur), l'évaluation de la calibration, la validation interne et l'analyse des courbes de décision. Le nomogramme a été construit à l'aide du package rms dans R.