Het studieprotocol is beoordeeld en goedgekeurd door de ethische commissie van het Santai County Hospital of Traditional Chinese Medicine in Mianyang, provincie Sichuan (goedkeuringsnummer 2025003). De studie is uitgevoerd in overeenstemming met de principes van de Verklaring van Helsinki. Voorafgaand aan de inclusie in de studie is schriftelijke geïnformeerde toestemming verkregen van elke patiënt of van een eerstegraads familielid. Deze retrospectieve cohortstudie bekeek de klinische dossiers van oudere patiënten die tussen mei 2023 en mei 2025 waren opgenomen voor acute exacerbatie van chronisch obstructieve longziekte (AECOPD) op de afdeling Respiratoire en Intensive Care Geneeskunde van het Santai County Traditional Chinese Medicine Hospital. De studie was erop gericht om routinematig beschikbare klinische voorspellers van pulmonale hypertensie (PH) te identificeren en een nomogram te ontwikkelen voor geïndividualiseerde PH-risicoschatting. De laboratoriuminstrumenten, het echocardiografiesysteem, de bloedafnamebenodigdheden, de instrumenten voor data-extractie en de statistische software die in deze workflow zijn gebruikt, staan vermeld in de Tabel met Materialen.
1. Patiëntscreening en beoordeling van geschiktheid
Gehospitaliseerde patiënten met een ontslagdiagnose van AECOPD gedurende de studieperiode werden gescreend via het elektronisch medisch dossier van de kliniek. Elk dossier werd beoordeeld om te bevestigen dat de patiënt ten minste 60 jaar oud was en dat de diagnose chronisch obstructieve longziekte overeenkwam met de diagnostische criteria van de Global Initiative for Chronic Obstructive Lung Disease uit 202115. AECOPD werd gedefinieerd als een acute verslechtering van de respiratoire symptomen die hospitalisatie en aanvullende behandeling vereiste. Voor patiënten met meer dan één AECOPD-opname tijdens de studieperiode werd alleen de eerste geschikte hospitalisatie meegenomen.
Patiënten werden geïncludeerd wanneer aan al de volgende criteria werd voldaan: een leeftijd van ten minste 60 jaar; hospitalisatie vanwege AECOPD tussen mei 2023 en mei 2025; beschikbare transthoracale echocardiografie tijdens de indexhospitalisatie; een beschikbare eerste nuchtere veneuze bloedafname na opname; en volledige klinische, laboratorium- en echocardiografische informatie die vereist was voor de modelconstructie. Patiënten werden geëxcludeerd indien zij leden aan ernstige leverinsufficiëntie, nierfalen in het eindstadium, een actieve maligniteit, pulmonale hypertensie die toeschrijfbaar was aan andere bekende oorzaken, primaire linkerhartziekte, bindweefselziekten zoals systemische sclerose, actieve pulmonale tuberculose, of wanneer cruciale variabelen ontbraken die vereist waren voor de uitkomstclassificatie of modelconstructie.
Voordat dossiers werden verwijderd vanwege onvolledige informatie, werden leeftijd, geslacht, rookgeschiedenis, drinkgeschiedenis, geschiedenis van drugallergie, hypertensie, diabetes mellitus, coronaire hartaandoening, astma, bronchiectasie, longinfectie, emfyseem, respiratoir falen, hartfalen, echocardiografische PH-status, indices van het volledig bloedbeeld, albumine, creatinine, hooggevoelig C-reactief proteïne, D-dimeer, fibrinogeen, B-type natriuretisch peptide en de neutrofiel-lymfocyt ratio gecontroleerd. Het aantal gescreende en uitgesloten dossiers, de redenen voor uitsluiting en de uiteindelijke geïncludeerde gevallen werden genoteerd om het stroomdiagram voor de selectie van patiënten op te stellen (Figuur 1).

Figuur 1: Patiëntscreening en cohortconstructie. Het stroomdiagram toont de identificatie van dossiers, de beoordeling van geschiktheid, de uitsluiting en de uiteindelijke groepering van oudere patiënten die waren opgenomen met AECOPD. Patiënten werden ingedeeld in de non-PH- of PH-groep op basis van transthoracale echocardiografische bevindingen. Afkortingen: AECOPD, acute exacerbatie van chronische obstructieve longziekte; PH, pulmonale hypertensie. Klik hier om een grotere versie van deze figuur te bekijken.
2. Klinische groepering en diagnostische definitie van pulmonale hypertensie
De geïncludeerde patiënten werden ingedeeld in een non-PH-groep of een PH-groep op basis van transthoracale echocardiografische bevindingen tijdens de ziekenhuisopname. Het diagnostisch kader was gebaseerd op de 2015-richtlijnen van de European Society of Cardiology/European Respiratory Society voor pulmonale hypertensie16. In deze retrospectieve dataset werd transthoracale echocardiografie in rust gebruikt als diagnostische basis. De onderzoeken werden uitgevoerd door getrainde echocardiografie-artsen met behulp van standaard parasternale, apicale en subcostale views. De interpreteerbaarheid van de echocardiografie werd beoordeeld vóór de groepsindeling, aangezien slechte akoestische vensters veel voorkomen bij patiënten met chronisch obstructieve longziekte. Een dossier werd als geschikt beschouwd voor PH-classificatie wanneer het Doppler-signaal van de tricuspidalisinsufficiëntie meetbaar was of wanneer de grootte van de rechterventrikel, de wanddikte en de systolische functie konden worden geëvalueerd op basis van interpreteerbare hartbeelden.
De systolische pulmonale arteriële druk werd geschat op basis van de piek-tricuspidalisregurgitatiesnelheid en de geschatte druk in de rechterboezem. De druk in de rechterboezem werd geschat aan hand van de diameter van de vena cava inferior en de inspiratoire collaps wanneer beide metingen beschikbaar waren. PH werd gedefinieerd als een geschatte systolische pulmonale arteriële druk van ten minste 40 mmHg. Wanneer de systolische pulmonale arteriële druk niet betrouwbaar kon worden geschat, vereiste de PH-classificatie indirect echocardiografisch bewijs van betrokkenheid van de rechterventrikel, waaronder vergroting van de rechterventrikel, hypertrofie van de rechterventrikel, systolische dysfunctie van de rechterventrikel, afvlakking van het interventriculaire septum of vergroting van de arteria pulmonalis. Echocardiografische verslagen met onduidelijke beelden, niet-traceerbare Doppler-signalen van tricuspidalisregurgitatie of inconsistente bevindingen in de rechterharthelft werden beoordeeld door een tweede echocardiograaf vóór de definitieve groepsindeling. Verslagen die na beoordeling oninterpreteerbaar bleven, werden uitgesloten van de definitieve analyse. In de uiteindelijke cohort voldeden 230 oudere patiënten aan de geschiktheidscriteria; 120 patiënten werden toegewezen aan de non-PH-groep en 110 patiënten werden toegewezen aan de PH-groep.
3. Extractie van demografische en klinische variabelen
Alle klinische gegevens werden opgehaald uit het elektronische medische dossier van de inpatient-afdeling met behulp van een gestandaardiseerd gegevensverzamelingsformulier. Geslacht, leeftijd, rookgeschiedenis, drinkgeschiedenis, geschiedenis van geneesmiddelenallergieën, hypertensie, diabetes mellitus, coronaire hartziekte, astma, bronchiectasie, longinfectie, emfyseem, respiratoir falen en hartfalen werden vastgelegd. Rookgeschiedenis werd gedefinieerd als een levenslange consumptie van ten minste 100 sigaretten, en drinkgeschiedenis werd gedefinieerd als regelmatige alcoholconsumptie van ten minste twee keer per week.
Longinfecties werden gediagnosticeerd door respiratoire symptomen, fysieke tekenen, bevindingen uit thoraxbeeldvorming en, indien beschikbaar, microbiologisch bewijs te integreren. Hartfalen werd gediagnosticeerd op basis van klinische manifestaties, fysieke tekenen, spiegels van B-type natriuretisch peptide en echocardiografische bevindingen. Informatie over comorbiditeit werd gecontroleerd aan de hand van ontslagdiagnoses, voortgangsverslagen, beeldvormingsrapporten, laboratoriumrapporten en consultatierapporten, indien beschikbaar.
4. Verzameling en verwerking van laboratoriumvariabelen
Het eerste nuchtere veneuze bloedmonster, afgenomen op de ochtend na opname in het ziekenhuis, werd gebruikt als bron voor de laboratoriumgegevens. Dit tijdstip kwam overeen met de eerste ochtendbloedafname na opname en vóór de routinematige aanpassing van de dagbehandeling, mits deze informatie beschikbaar was in het medisch dossier. Indices van het volledige bloedbeeld, biochemische indices, stollingsmarkers, ontstekingsmarkers en markers voor cardiale stress werden geëxtraheerd uit het laboratoriuminformatiesysteem van het ziekenhuis.
Het aantal neutrofielen, het aantal lymfocyten, het aantal bloedplaatjes, het gemiddelde plaatjesvolume en de variatiecoëfficiënt van de breedte van de rode bloedcelverdeling werden genoteerd op basis van het rapport van het volledig bloedbeeld. De neutrofiel-lymfocyt-ratio werd als volgt berekend:

Albumin, creatinine, D-dimeer, fibrinogeen, hooggevoelig C-reactief proteïne en B-type natriuretisch peptide werden genoteerd op basis van de betreffende laboratoriumrapporten. Er werden consistente meeteenheden gebruikt in de gehele dataset: albumin in g/L, creatinine in µmol/L, hooggevoelig C-reactief proteïne in mg/L, D-dimeer in µg/mL, fibrinogeen in g/L, neutrofielenaantal, lymfocytenaantal en bloedplaatjesaantal in ×109/L, gemiddeld bloedplaatjesvolume in fL en B-type natriuretisch peptide in pg/mL.
5. Opschoning en voorbereiding van gegevens
De dataset werd gecontroleerd voorafgaand aan de statistische analyse. Elk opgenomen record werd gecontroleerd op geschiktheid, groepsindeling, volledigheid van laboratoriumgegevens, interpreteerbaarheid van de echocardiografie en de beschikbaarheid van potentiële predictoren. Directe persoonlijke identificatiekenmerken werden verwijderd en elke patiënt kreeg een studienummer toegewezen. Dubbele records, herhaalde hospitalisaties, inconsistente groepslabels, onmogelijke data en onwaarschijnlijke laboratoriumwaarden werden gecontroleerd aan de hand van het oorspronkelijke elektronisch medisch dossier en het laboratoriumverslag. Er werd geen imputatie uitgevoerd voor belangrijke uitkomst- of predictormatige variabelen, aangezien het definitieve model gebaseerd was op een complete-case analyse.
Continue variabelen werden geëvalueerd op hun distributiekarakteristieken. Normaal verdeelde variabelen werden samengevat als gemiddelde ± standaarddeviatie; niet-normaal verdeelde variabelen als mediaan met interkwartielafstand; en categorische variabelen als aantallen en percentages. Voorafgaand aan de multivariabele modellering werden klinisch gerelateerde variabelen onderzocht op redundantie. Het aantal neutrofielen en het aantal lymfocyten werden samen met de neutrofiel-lymfocyt ratio geëvalueerd; pulmonale infectie werd geëvalueerd samen met het aantal neutrofielen, hooggevoelig C-reactief proteïne en de neutrofiel-lymfocyt ratio; hartfalen werd geëvalueerd samen met B-type natriuretisch peptide en echocardiografische bevindingen van het hart; D-dimeer en fibrinogeen werden geëvalueerd als stollingsgerelateerde variabelen; en het aantal bloedplaatjes, het gemiddelde bloedplaatjesvolume en de variatiecoëfficiënt van de breedte van de erytrocytendistributie werden geëvalueerd als hematologische variabelen.
6. Univariante vergelijking tussen groepen
Demografische kenmerken, comorbiditeiten, laboratoriumindicatoren en echocardiografische groeperingsvariabelen werden vergeleken tussen de non-PH- en PH-groepen. Voor normaal verdeelde continue variabelen werd een t-toets voor onafhankelijke steekproeven gebruikt. Voor niet-normaal verdeelde continue variabelen werd de Mann-Whitney U-toets gebruikt. Voor categorische variabelen werd de χ2-toets of de exacte toets van Fisher gebruikt, afhankelijk van de verwachte celaantallen. Een tweezijdige P-waarde kleiner dan 0,05 werd beschouwd als statistisch significant. Variabelen met significante verschillen tussen de groepen werden, samen met klinisch relevante predictoren, meegenomen voor multivariabele evaluatie.
7. Multivariabele logistische regressieanalyse
Binaire logistische regressie werd gebruikt om onafhankelijke predictoren te identificeren die geassocieerd zijn met PH bij oudere patiënten met AECOPD. De PH-status werd gebruikt als de afhankelijke variabele. Continue predictoren werden ingevoerd in hun oorspronkelijke meeteenheden, en categorische predictoren werden gecodeerd als binaire variabelen. Kandidaatvariabelen werden geselecteerd op basis van univariaat resultaten en klinische relevantie. Variabelen met P < 0.05 in de univariate vergelijking werden eerst geïdentificeerd en vervolgens beoordeeld op redundantie, interpreteerbaarheid en klinische overlap.
De initiële kandidaatvariabelen omvatten hooggevoelig C-reactief proteïne, D-dimeer, de variatiecoëfficiënt van de breedte van de rode bloedceldistributie, het aantal bloedplaatjes, het gemiddelde bloedplaatjesvolume, de neutrofiel-lymfocytenratio, B-type natriuretisch peptide, coronaire hartziekten en astma in de voorgeschiedenis. Het aantal neutrofielen en het aantal lymfocyten werden niet samen met de neutrofiel-lymfocytenratio ingevoerd, omdat zij componenten van deze ratio waren. Longinfecties en hartfalen werden geëvalueerd tijdens de redundantiebeoordeling. Om te onderzoeken of hun uitsluiting de interpretatie van het model beïnvloedde, werd een sensitiviteitsanalyse uitgevoerd door longinfecties en hartfalen toe te voegen aan de primaire set voorspellers. Regressieschattingen, discriminatie, kalibratie en het klinische netto voordeel werden vergeleken tussen het primaire model en het sensitiviteitsmodel.
Multicollineariteit tussen kandidaat-predictoren werd beoordeeld aan de hand van variantie-inflatiefactoren. Categorische variabelen werden gecodeerd vóór de berekening, en continue variabelen werden ingevoerd in dezelfde eenheden als gebruikt voor de regressiemodellering. Een variantie-inflatiefactor lager dan 5 werd beschouwd als een indicatie dat er geen sprake was van ernstige multicollineariteit. Regressiecoëfficiënten, standaardfouten, Wald χ2-waarden, odds ratio's, 95% betrouwbaarheidsintervallen, P-waarden en variantie-inflatiefactoren werden gerapporteerd voor de variabelen die waren opgenomen in het multivariabele model.
8. Constructie van het nomogram
Het nomogram werd opgesteld op basis van de onafhankelijke predictoren die waren behouden in het multivariabele logistische regressiemodel. De uiteindelijke predictoren waren D-dimeer, de neutrofiel-lymfocytenratio, B-type natriuretisch peptide en astma in de medische anamnese. Het logistische regressiemodel werd omgezet in een op punten gebaseerd nomogram met behulp van R-software en het rms-pakket. Er werd een datadistributie-object gegenereerd voor de modelleringsdataset, het logistische regressiemodel werd gefit en de regressiecoëfficiënten werden gebruikt om punten toe te kennen aan elke predictor. Voor elke patiënt werden de punten voor D-dimeer, de neutrofiel-lymfocytenratio, B-type natriuretisch peptide en astma in de anamnese opgeteld om een totaalscore te genereren, die vervolgens werd gekoppeld aan de geschatte waarschijnlijkheid van PH.
Het nomogram werd geïnterpreteerd aan de hand van standaard puntgebaseerde procedures die zijn beschreven in eerdere studies naar klinische predictiemodellen17,18. De D-dimeerwaarde, de neutrofiel-lymfocytenratio, de B-type natriuretisch peptide-waarde en de astmastatus van de patiënt werden op de overeenkomstige assen gelokaliseerd. Van elke predictorwaarde werd een verticale lijn getrokken naar de puntenschaal. De punten werden bij elkaar opgeteld om de totale score te verkrijgen, waarna een verticale lijn van de schaal voor de totale score naar de as voor de voorspelde waarschijnlijkheid werd getrokken om de geïndividualiseerde geschatte waarschijnlijkheid van PH te bepalen. De illustratieve patiënt in het manuscript is geselecteerd op basis van waarden binnen of nabij de geobserveerde klinische distributie van de cohort.
9. Evaluatie van modelprestaties
Modeldiscriminatie werd geëvalueerd met behulp van een receiver operating characteristic (ROC)-curve-analyse en de oppervlakte onder de curve (area under the curve). De oppervlakte onder de curve werd gerapporteerd met een 95% betrouwbaarheidsinterval. Het nomogram werd vergeleken met individuele predictoren en met een model waarin kwantitatieve laboratoriummarkers waren gecombineerd. De kalibratie werd beoordeeld door voorspelde waarschijnlijkheden te vergelijken met waargenomen PH-uitkomsten. Er werd een kalibratiecurve gegenereerd middels bootstrap-resampling. Het kalibratie-intercept, de kalibratieshelling, de gemiddelde absolute fout, de Brier-score en de Hosmer-Lemeshow-test werden gerapporteerd. Ook het bootstrap-gecorrigeerde kalibratie-intercept en de bootstrap-gecorrigeerde kalibratieshelling werden berekend. De Hosmer-Lemeshow-test werd geïnterpreteerd in samenhang met de kalibratiecurve en de kwantitatieve kalibratie-indices.
Het klinisch nut werd geëvalueerd met behulp van decision curve-analyse. Het netto voordeel werd uitgezet tegenover drempelwaarschijnlijkheden en vergeleken met referentiestrategieën voor 'iedereen behandelen' en 'niemand behandelen'. Het bereik van de drempelwaarschijnlijkheid waarin het model een hoger netto voordeel bood dan de referentiestrategieën werd vastgesteld.
10. Interne validatie en reproduceerbaarheid
Interne validatie werd uitgevoerd met bootstrap-resampling met 1.000 replicaat. In elk bootstrap-replicaat werd het model opnieuw gefit en getest om het optimisme in de modelprestaties te schatten. Optimisme-gecorrigeerde discriminatie- en kalibratie-indexen werden berekend en gerapporteerd samen met de schijnbare modelprestaties.
Statistische analyses werden uitgevoerd met standaard statistische software en R versie 4.5.2. R-pakketten werden gebruikt voor logistische regressiemodellering, de constructie van het nomogram, receiver operating characteristic curve-analyse, kalibratiebeoordeling, interne validatie en decision curve-analyse. Het nomogram werd geconstrueerd met het rms-pakket in R.