Das Studienprotokoll wurde vom Ethikkomitee des Krankenhauses für Traditionelle Chinesische Medizin der Kreisstadt Santai in der Provinz Mianyang, Sichuan (Genehmigungsnummer 2025003), geprüft und genehmigt. Die Studie wurde gemäß den Grundsätzen der Deklaration von Helsinki durchgeführt. Vor der Aufnahme in die Studie wurde bei jedem Patienten oder einem Verwandten ersten Grades schriftlich informierte Einwilligung eingeholt. Diese retrospektive Kohortenstudie analysierte die klinischen Unterlagen älterer Patienten, die zwischen Mai 2023 und Mai 2025 aufgrund einer akuten Verschlechterung einer chronisch obstruktiven Lungenerkrankung (AECOPD) in der Abteilung für Pneumologie und Intensivmedizin des Krankenhauses für Traditionelle Chinesische Medizin der Kreisstadt Santai hospitalisiert wurden. Ziel der Studie war es, routinemäßig verfügbare klinische Prädiktoren für pulmonale Hypertonie (PH) zu identifizieren und ein Nomogramm zur individuellen Einschätzung des PH-Risikos zu entwickeln. Die in diesem Arbeitsablauf verwendeten Laborgeräte, das Echokardiographiesystem, die Blutentnahmematerialien, die Datengewinnungstools sowie die statistische Software sind in der Tabelle der Materialien aufgeführt.
1. Patienten-Screening und Eignungsbeurteilung
Hospitalisierte Patienten mit einer Entlassungsdiagnose einer akuten Exazerbation der COPD (AECOPD) während des Studienzeitraums wurden über das elektronische Krankenaktensystem der Stationen gescreent. Jeder Datensatz wurde überprüft, um sicherzustellen, dass der Patient mindestens 60 Jahre alt war und die Diagnose einer chronisch obstruktiven Lungenerkrankung den diagnostischen Kriterien der Global Initiative for Chronic Obstructive Lung Disease (GOLD) aus dem Jahr 2021 entsprach15. Eine AECOPD wurde definiert als eine akute Verschlechterung der Atemwegssymptome, die eine Hospitalisierung und zusätzliche Therapie erforderte. Bei Patienten mit mehr als einer AECOPD-Aufnahme während des Studienzeitraums wurde nur die erste berechtigte Hospitalisierung einbezogen.
Patienten wurden eingeschlossen, wenn alle folgenden Kriterien erfüllt waren: Alter von mindestens 60 Jahren; Krankenhauseinweisung wegen einer akuten Exazerbation der COPD (AECOPD) zwischen Mai 2023 und Mai 2025; verfügbare transthorakale Echokardiographie während der Indexhospitalisierung; verfügbare erste nüchterne venöse Blutprobe nach der Aufnahme; sowie vollständige klinische, laborchemische und echokardiographische Informationen, die für die Modellkonstruktion erforderlich sind. Patienten wurden ausgeschlossen, wenn sie eine schwere hepatische Insuffizienz, eine Nierenerkrankung im Endstadium, eine aktive maligne Erkrankung, eine pulmonale Hypertonie aufgrund anderer bekannter Ursachen, eine primäre Linksherzerkrankung, eine Bindegewebserkrankung wie systemische Sklerose, eine aktive pulmonale Tuberkulose oder fehlende Schlüsselvariablen, die für die Outcome-Klassifizierung oder die Modellkonstruktion erforderlich sind, aufwiesen.
Bevor Datensätze aufgrund unvollständiger Angaben, Alter, Geschlecht, Rauchgeschichte, Alkoholkonsumgeschichte, Arzneimittelallergiegeschichte, Hypertonie, Diabetes mellitus, koronare Herzkrankheit, Asthma, Bronchiektasen, pulmonale Infektion, Emphysem, respiratorische Insuffizienz, Herzinsuffizienz, echokardiographischem PH-Status, Parametern des vollständigen Blutbilds, Albumin, Kreatinin, hochsensitivem C-reaktivem Protein, D-Dimer, Fibrinogen, natriuretischem Peptid vom Typ B und Neutrophile-zu-Lymphozyten-Verhältnis entfernt wurden, wurden diese überprüft. Die Anzahl der geprüften und ausgeschlossenen Datensätze, die Ausschlussgründe sowie die endgültig eingeschlossenen Fälle wurden dokumentiert, um das Flussdiagramm zur Patientenauswahl zu erstellen (Abbildung 1).

Abbildung 1: Patientenscreening und Kohortenbildung. Das Flussdiagramm zeigt die Identifizierung von Datensätzen, die Beurteilung der Eignung, Ausschlüsse und die endgültige Gruppierung älterer Patienten, die wegen einer akuten Exazerbation einer COPD (AECOPD) hospitalisiert wurden. Die Patienten wurden anhand der Ergebnisse der transthorakalen Echokardiographie in die Gruppe ohne pulmonale Hypertonie (non-PH) oder die Gruppe mit pulmonaler Hypertonie (PH) eingeteilt. Abkürzungen: AECOPD, akute Exazerbation der chronisch obstruktiven Lungenerkrankung; PH, pulmonale Hypertonie. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
2. Klinische Einteilung und diagnostische Definition der pulmonalen Hypertonie
Die eingeschlossenen Patienten wurden anhand der Ergebnisse der transthorakalen Echokardiographie während des Krankenhausaufenthalts in eine Nicht-PH-Gruppe oder eine PH-Gruppe klassifiziert. Der diagnostische Rahmen basierte auf den Leitlinien der Europäischen Gesellschaft für Kardiologie/Europäische Gesellschaft für Pneumologie aus dem Jahr 2015 für pulmonale Hypertonie16. Die ruhende transthorakale Echokardiographie diente als diagnostische Grundlage in diesem retrospektiven Datensatz. Die Untersuchungen wurden von geschulten Echokardiographie-Ärzten unter Verwendung standardmäßiger parasternaler, apikaler und subkostaler Ansichten durchgeführt. Die Interpretierbarkeit der Echokardiographie wurde vor der Gruppeneinteilung überprüft, da schlechte akustische Fenster bei Patienten mit chronisch obstruktiver Lungenerkrankung häufig vorkommen. Eine Aufzeichnung galt als geeignet für die Klassifizierung der pulmonalen Hypertonie, wenn das Doppler-Signal der Trikuspidalklappeninsuffizienz messbar war oder wenn die Größe des rechten Ventrikels, die Wanddicke und die systolische Funktion anhand interpretierbarer Herzansichten beurteilt werden konnten.
Der systolische pulmonalarterielle Druck wurde anhand der maximalen Trikuspidalinsuffizienzgeschwindigkeit und des geschätzten rechtsatrialen Drucks ermittelt. Der rechtsatriale Druck wurde aus dem Durchmesser der unteren Hohlvene und deren inspiratorischem Kollaps geschätzt, sofern beide Messungen verfügbar waren. Eine pulmonale Hypertonie (PH) wurde definiert als ein geschätzter systolischer pulmonalarterieller Druck von mindestens 40 mmHg. Wenn der systolische pulmonalarterielle Druck nicht zuverlässig geschätzt werden konnte, war für die Klassifizierung der PH indirekter echokardiographischer Befundnachweis einer Rechtsherzbeteiligung erforderlich, einschließlich Rechtsherzvergrößerung, Rechtsherzhypertrophie, systolischer Rechtsherzinsuffizienz, Flattening des interventrikulären Septums oder Erweiterung der Pulmonalarterie. Echokardiographische Befunde mit unscharfen Bildern, nicht nachvollziehbaren Doppler-Signalen der Trikuspidalinsuffizienz oder inkonsistenten Befunden der Rechtsherzstrukturen wurden vor der endgültigen Gruppeneinteilung von einem zweiten echokardiographisch tätigen Arzt überprüft. Befunde, die auch nach dieser Überprüfung nicht interpretierbar blieben, wurden von der endgültigen Analyse ausgeschlossen. In der endgültigen Kohorte erfüllten 230 ältere Patienten die Einschlusskriterien; 120 Patienten wurden der Nicht-PH-Gruppe und 110 Patienten der PH-Gruppe zugeordnet.
3. Extraktion demografischer und klinischer Variablen
Alle klinischen Daten wurden mithilfe eines standardisierten Erfassungsbogens aus dem elektronischen Patienten-Aktensystem der stationären Behandlung abgerufen. Es wurden Geschlecht, Alter, Rauchgeschichte, Alkoholkonsumgeschichte, Arzneimittelallergie-Anamnese, Hypertonie, Diabetes mellitus, koronare Herzkrankheit, Asthma, Bronchiektasen, pulmonale Infektion, Emphysem, Atemversagen und Herzinsuffizienz erfasst. Eine Rauchgeschichte wurde definiert als ein lebenslanger Konsum von mindestens 100 Zigaretten, eine Alkoholkonsumgeschichte als regelmäßiger Alkoholkonsum von mindestens zweimal pro Woche.
Eine pulmonale Infektion wurde durch die Integration von respiratorischen Symptomen, körperlichen Befunden, bildgebenden Thoraxbefunden und, falls verfügbar, mikrobiologischen Nachweisen diagnostiziert. Eine Herzinsuffizienz wurde entsprechend den klinischen Manifestationen, körperlichen Zeichen, den Werten des B-Typ-natriuretischen Peptids und echokardiographischen Befunden diagnostiziert. Informationen zu Komorbiditäten wurden anhand von Entlassungsdiagnosen, Verlaufsnotizen, bildgebenden Befunden, Laborbefunden und Konsiliarberichten, falls verfügbar, gegengeprüft.
4. Sammlung und Aufbereitung von Laborvariablen
Die erste nüchterne venöse Blutprobe, die am Morgen nach der Krankenhausaufnahme entnommen wurde, diente als Quelle für die Laborergebnisse. Dieser Zeitpunkt entsprach der ersten morgendlichen Blutentnahme nach der Aufnahme und vor der routinemäßigen Anpassung der Tagesbehandlung, sofern diese Informationen im medizinischen Befund verfügbar waren. Blutbildparameter, biochemische Parameter, Gerinnungsmarker, Entzündungsmarker und Marker für kardiale Belastung wurden aus dem Laborinformationssystem des Krankenhauses extrahiert.
Die Anzahl der Neutrophilen, die Anzahl der Lymphozyten, die Thrombozytenanzahl, das mittlere Thrombozytenvolumen und der Variationskoeffizient der roten Blutkörperchenverteilungsbreite wurden aus dem vollständigen Blutbildbericht erfasst. Das Neutrophilen-Lymphozyten-Verhältnis wurde wie folgt berechnet:

Albumin, Kreatinin, D-Dimer, Fibrinogen, hochsensitives C-reaktives Protein und BNP wurden aus den entsprechenden Laborberichten entnommen. Durchgängig wurden folgende Maßeinheiten im Datensatz verwendet: Albumin in g/L, Kreatinin in µmol/L, hochsensitives C-reaktives Protein in mg/L, D-Dimer in µg/mL, Fibrinogen in g/L, Neutrophile, Lymphozyten und Thrombozyten in ×109/L, mittleres Thrombozytenvolumen in fL sowie BNP in pg/mL.
5. Datenaufbereitung und Vorbereitung
Der Datensatz wurde vor der statistischen Analyse überprüft. Jeder eingeschlossene Datensatz wurde auf Eignung, Gruppenzugehörigkeit, Vollständigkeit der Laborergebnisse, Interpretierbarkeit der Echokardiographie und Verfügbarkeit der Kandidatenprädiktoren geprüft. Direkte personenbezogene Identifikatoren wurden entfernt, und jedem Patienten wurde eine Studien-ID zugewiesen. Doppelte Einträge, wiederholte Krankenhausaufenthalte, inkonsistente Gruppenbezeichnungen, unmögliche Daten und unplausible Laborwerte wurden mit der ursprünglichen elektronischen Patientenakte und dem Laborbericht abgeglichen. Es erfolgte keine Imputation für wesentliche Zielgrößen oder Prädiktorvariablen, da das endgültige Modell auf einer Analyse vollständiger Fälle basierte.
Stetige Variablen wurden hinsichtlich ihrer Verteilungseigenschaften überprüft. Normalverteilte Variablen wurden als Mittelwert ± Standardabweichung zusammengefasst; nicht normalverteilte Variablen als Median mit Interquartilsabstand; und kategoriale Variablen als Anzahlen und Prozentanteile. Vor der multivariablen Modellierung wurden klinisch verwandte Variablen auf Redundanz geprüft. Die Neutrophilenanzahl und die Lymphozytenanzahl wurden gemeinsam mit dem Neutrophilen-zu-Lymphozyten-Verhältnis betrachtet; die pulmonale Infektion wurde gemeinsam mit der Neutrophilenanzahl, hochsensitivem C-reaktivem Protein und dem Neutrophilen-zu-Lymphozyten-Verhältnis bewertet; die Herzinsuffizienz wurde gemeinsam mit dem B-Typ-natriuretischen Peptid und echokardiographischen Herzbefunden analysiert; D-Dimer und Fibrinogen wurden als gerinnungsrelevante Variablen betrachtet; und Thrombozytenanzahl, mittleres Thrombozytenvolumen sowie Koeffizient der Variation der roten Zellverteilungsbreite wurden als hämatologische Variablen ausgewertet.
6. Univariativer Vergleich zwischen Gruppen
Die demografischen Merkmale, Komorbiditäten, Laborparameter und die Gruppierungsechokardiographie-Variablen wurden zwischen der Nicht-PH- und der PH-Gruppe verglichen. Für normalverteilte stetige Variablen wurde ein t-Test für unabhängige Stichproben verwendet. Für nicht normalverteilte stetige Variablen kam der Mann-Whitney-U-Test zum Einsatz. Für kategoriale Variablen wurde je nach erwarteten Zellenhäufigkeiten der χ2-Test oder der exakte Fisher-Test angewendet. Ein zweiseitiger P-Wert unter 0,05 galt als statistisch signifikant. Variablen mit signifikanten Unterschieden zwischen den Gruppen wurden zusammen mit klinisch relevanten Prädiktoren in die multivariable Analyse einbezogen.
7. Multivariable logistische Regressionsanalyse
Die binäre logistische Regression wurde verwendet, um unabhängige Prädiktoren für pulmonale Hypertonie (PH) bei älteren Patienten mit AECOPD zu identifizieren. Der PH-Status wurde als abhängige Variable verwendet. Stetige Prädiktoren wurden in ihren ursprünglichen Messgrößen eingegeben, und kategoriale Prädiktoren wurden als binäre Variablen kodiert. Die Auswahl der Kandidatenvariablen erfolgte basierend auf univariaten Ergebnissen und klinischer Relevanz. Variablen mit einem P-Wert < 0,05 im univariaten Vergleich wurden zunächst identifiziert und anschließend auf Redundanz, Interpretierbarkeit und klinische Überschneidungen überprüft.
Zu den anfänglichen Kandidatenvariablen gehörten hochsensitives C-reaktives Protein, D-Dimer, Koeffizient der Variation der roten Blutkörperchenverteilungsbreite, Thrombozytenzahl, mittleres Thrombozytenvolumen, Neutrophile-zu-Lymphozyten-Quotient, B-Typ natriuretisches Peptid, koronare Herzkrankheit und Asthmaanamnese. Neutrophile und Lymphozyten wurden nicht zusammen mit dem Neutrophile-zu-Lymphozyten-Quotienten berücksichtigt, da sie Bestandteile dieses Quotienten sind. Pulmonale Infektion und Herzinsuffizienz wurden im Rahmen der Redundanzprüfung evaluiert. Um zu untersuchen, ob deren Ausschluss die Modellinterpretation beeinflusste, wurde eine Sensitivitätsanalyse durchgeführt, bei der pulmonale Infektion und Herzinsuffizienz dem primären Prädiktorensatz hinzugefügt wurden. Regressionsabschätzungen, Diskriminierung, Kalibrierung und klinischer Nutzen wurden zwischen dem primären und dem Sensitivitätsmodell verglichen.
Multikollinearität zwischen den Kandidatenprädiktoren wurde mithilfe von Varianzinflationsfaktoren bewertet. Kategoriale Variablen wurden vor der Berechnung codiert, und kontinuierliche Variablen wurden in den gleichen Einheiten eingegeben, die für die Regressionsmodellierung verwendet wurden. Ein Varianzinflationsfaktor unter 5 galt als Hinweis auf keine starke Multikollinearität. Regressionskoeffizienten, Standardfehler, Wald-χ2-Werte, Odds-Ratios, 95-%-Konfidenzintervalle, P-Werte und Varianzinflationsfaktoren wurden für die in das multivariable Modell aufgenommenen Variablen angegeben.
8. Nomogrammkonstruktion
Das Nomogramm wurde aus den unabhängigen Prädiktoren erstellt, die im multivariablen logistischen Regressionsmodell beibehalten wurden. Die endgültigen Prädiktoren waren D-Dimer, Neutrophile-zu-Lymphozyten-Verhältnis, B-Typ natriuretisches Peptid und Asthmaanamnese. Das logistische Regressionsmodell wurde mithilfe der R-Software und des rms-Pakets in ein punktbasiertes Nomogramm umgewandelt. Ein Datenverteilungsobjekt wurde für den Modellierungsdatensatz erstellt, das logistische Regressionsmodell wurde angepasst, und die Regressionskoeffizienten wurden verwendet, um jedem Prädiktor Punkte zuzuweisen. Bei jedem Patienten wurden die Punkte für D-Dimer, Neutrophile-zu-Lymphozyten-Verhältnis, B-Typ natriuretisches Peptid und Asthmaanamnese summiert, um eine Gesamtpunktzahl zu ermitteln, die anschließend der geschätzten Wahrscheinlichkeit für PH zugeordnet wurde.
Das Nomogramm wurde mithilfe standardmäßiger, punktbasierter Verfahren interpretiert, wie sie in früheren Studien zu klinischen Vorhersagemodellen beschrieben wurden17,18. Der D-Dimer-Wert des Patienten, das Verhältnis von Neutrophilen zu Lymphozyten, der Wert des B-Typ-natriuretischen Peptids sowie der Asthmastatus wurden auf den entsprechenden Achsen lokalisiert. Von jedem Prädiktorwert wurde eine senkrechte Linie zur Punkteskala gezogen. Die ermittelten Punkte wurden summiert, um die Gesamtpunktzahl zu erhalten, und von dieser Gesamtpunktzahl wurde eine senkrechte Linie zur Achse der vorhergesagten Wahrscheinlichkeit gezogen, um die individuell geschätzte Wahrscheinlichkeit für PH zu ermitteln. Der im Manuskript dargestellte Beispielpatient wurde anhand von Werten ausgewählt, die innerhalb oder nahe der beobachteten klinischen Verteilung der Kohorte lagen.
9. Bewertung der Modellleistung
Die Modell-Diskriminierung wurde mittels Analyse der Receiver-Operating-Characteristic-Kurve und der Fläche unter der Kurve bewertet. Die Fläche unter der Kurve wurde mit einem 95 %igen Konfidenzintervall angegeben. Das Nomogramm wurde mit einzelnen Prädiktoren sowie mit einem Modell verglichen, das quantitative laborchemische Marker kombiniert. Die Kalibrierung wurde durch den Vergleich vorhergesagter Wahrscheinlichkeiten mit beobachteten PH-Ergebnissen überprüft. Eine Kalibrierungskurve wurde mittels Bootstrap-Resampling erstellt. Kalibrierungsinterzept, Kalibrierungssteigung, mittlerer absoluter Fehler, Brier-Score und der Hosmer-Lemeshow-Test wurden berichtet. Zudem wurden der bootstrap-korrigierte Kalibrierungsinterzept und die bootstrap-korrigierte Kalibrierungssteigung berechnet. Der Hosmer-Lemeshow-Test wurde gemeinsam mit der Kalibrierungskurve und quantitativen Kalibrierungsindizes interpretiert.
Die klinische Nützlichkeit wurde mittels Entscheidungskurvenanalyse bewertet. Der Nettogewinn wurde über die Schwellenwahrscheinlichkeiten aufgetragen und mit den Referenzstrategien „alle behandeln“ und „keinen behandeln“ verglichen. Der Bereich der Schwellenwahrscheinlichkeit, in dem das Modell einen höheren Nettogewinn als die Referenzstrategien lieferte, wurde identifiziert.
10. Interne Validierung und Reproduzierbarkeit
Die interne Validierung erfolgte mittels Bootstrap-Resampling mit 1.000 Wiederholungen. Bei jedem Bootstrap-Durchlauf wurde das Modell erneut angepasst und getestet, um den Optimismus bezüglich der Modellleistung zu schätzen. Optimismus-korrigierte Diskriminierungs- und Kalibrationsindizes wurden berechnet und zusammen mit der scheinbaren Modellleistung angegeben.
Die statistischen Analysen wurden mit herkömmlicher statistischer Software und R-Version 4.5.2 durchgeführt. R-Pakete wurden für die logistische Regressionsmodellierung, die Erstellung von Nomogrammen, die Analyse der Receiver-Operating-Characteristic-Kurve, die Bewertung der Kalibrierung, die interne Validierung sowie die Entscheidungskurvenanalyse verwendet. Das Nomogramm wurde mit dem rms-Paket in R erstellt.