$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Il presente studio è stato approvato dal Comitato Etico del Secondo Ospedale Popolare di Shenzhen (Protocollo n. 20193357024). Il consenso informato è stato ottenuto da tutti i partecipanti o dai loro rappresentanti legali prima dell'iscrizione.
Popolazione e metodologia di studio
Questo studio è stato progettato come uno studio di coorte retrospettivo a centro singolo. Sono stati esaminati i referti medici dei pazienti con diagnosi primaria di BPCO presso il Dipartimento di Medicina Respiratoria del Shenzhen Second People's Hospital, da gennaio 2021 a giugno 2022. Tutti i pazienti adulti (≥18 anni) con una diagnosi primaria di BPCO sono stati inizialmente sottoposti a screening utilizzando la codifica e la revisione delle cartelle cliniche della Classificazione Internazionale delle Malattie (ICD-10) dal sistema elettronico di cartelle cliniche (EMR) dell'ospedale.
Criteri di inclusione
(1) Diagnosi confermata di BPCO in conformità con le linee guida della Global Initiative for chronic obstructive polmonary disease (GOLD); (2) Disponibilità di tomografia computerizzata ad alta risoluzione (HRCT) del torace eseguita entro una settimana dal ricovero; (3) Disponibilità di spirometria completa e dati di laboratorio; e (4) Almeno un anno di dati di follow-up per il monitoraggio delle esacerbazioni acute.
Criteri di esclusione
(1) Infezioni polmonari attive (ad esempio, polmonite o tubercolosi) al momento dell'imaging HRCT; (2) Malignità polmonare convivente; (3) Interventi chirurgici toracici precedenti con potenziale impatto sull'anatomia delle vie aeree; e (4) Dati clinici critici mancanti o immagini non valutabili a causa di artefatti di movimento. Dopo l'applicazione di questi criteri, è stata arruolata una coorte finale di 212 pazienti, composta da 47 pazienti nel gruppo mucoso positivo (MP) e 165 pazienti nel gruppo non-muco plug (NMP). Immagini rappresentative HRCT sono illustrate nella Figura 1. I pazienti del gruppo NMP (n = 165) hanno svolto il ruolo di controlli interni, permettendo il confronto statistico di caratteristiche cliniche, indici di funzione polmonare e biomarcatori di laboratorio tra gruppi. Tutte le analisi sono state condotte su questa coorte controllata internamente per supportare lo sviluppo di modelli guidati da ipotesi.
Raccolta dati
L'estrazione dei dati seguiva un protocollo strutturato e sequenziale. Le variabili demografiche raccolte includevano età, sesso, indice di massa corporea (BMI) e stato di fumo. Le variabili della storia clinica comprendevano la durata della BPCO, la frequenza di esacerbazioni acute e le comorbidità. Parametri di spirometria recuperati tra cui FEV1%, FEV1 a capacità vitale forzata (FVC), capacità vitale (VC), flusso espiratorio forzato (FEF25–75%pred), volume residuo (RV), capacità polmonare totale (TLC) e rapporto RV/TLC. Gli indici di laboratorio includevano immunoglobulina totale E (IgE) sirica, 25-idrossivitamina D(25(OH)D), calcio sierico (Ca2+), fosforo, antigene carboidrato (CA199), ossido nitrico espirato frazionario (FeNO) e ossido nitrico conduttore delle vie aeree (CaNO). Lo screening per la comorbidità includeva sinusite, asma, bronchiectasi, colonizzazione fungina e batterica, e malattie cardiovascolari e metaboliche. Tutti i dati sono stati recuperati dal sistema elettronico di cartelle cliniche (EMR) dell'ospedale. Le immagini HRCT sono state accessibili dall'archivio del sistema di archiviazione e comunicazione delle immagini (PACS) dell'ospedale. I dettagli del software e delle attrezzature utilizzate in questo studio sono forniti nella Tabella dei Materiali. Non venivano utilizzati reagenti fisici né materiali di laboratorio; Tutte le analisi sono state effettuate utilizzando dati clinici e radiologici esistenti. Tutti i dati dei pazienti sono stati esaminati da due investigatori indipendenti. I dati mancanti venivano gestiti utilizzando il metodo di imputazione non parametrica 'missForest' implementato in R, per minimizzare la distorsione nelle analisi multivariate.
Criteri diagnostici HRCT per tappi mucosi
Tutti i pazienti hanno sottoposto una HRCT utilizzando protocolli di imaging istituzionali standardizzati. I tappi di muco sono stati definiti radiologicamente su tacce assiali TAC identificate come strutture tubolari o ramificate di attenuazione dei tessuti molli che occupano un lumen delle vie aeree, visibili su almeno due fette assiali contigue, coerenti con i criteri diagnostici pubblicati. Solo i casi con opacità chiaramente demarcate, segmentarie o subsegmentali delle vie aeree con attenuazione dei tessuti molli simile a quella dei tessuti molli e non attribuibili a artefatti o bronchiectasi sono stati etichettati come positivi al tappi mucosi. L'imaging HRCT è stato eseguito utilizzando uno scanner CT Siemens SOMATOM Definition AS (128-slice) con i seguenti parametri di acquisizione: spessore della fetta 1,0 mm, intervallo di ricostruzione 0,75 mm e uso del kernel ad alta risoluzione B70f. Le immagini sono state esaminate nelle impostazioni standard delle finestre polmonari (larghezza della finestra: 1600 unità di Hounsfield [HU]; Livello finestra: 600 HU. Due radiologi toracici certificati con oltre 8 anni di esperienza hanno esaminato indipendentemente tutte le scansioni. I casi con discrepanze interpretative sono stati risolti tramite discussione consensuale. I criteri diagnostici sono stati applicati in modo uniforme in tutti i casi per garantire la coerenza della classificazione.
Costruzione, valutazione e validazione del nomogramma
È stato sviluppato un nomogramma per prevedere i tappi di muco rilevati dalla TC nei pazienti con BPCO basandosi sui risultati di regressione logistica multivariata. Il modello finale ha incorporato i seguenti predittori indipendenti: bronchiectasi, rinosinusite cronica (CRS), esacerbazioni acute (AE), IMC, FEF25–75% pred, rapporto RV/TLC e livelli sierici di 25(OH)D. A ogni predittore viene assegnato un punteggio su una scala orizzontale; I punteggi individuali vengono sommati per ottenere un punteggio totale, che corrisponde a una probabilità prevista di presenza di tappi mucosi sulla scala di probabilità di output. Il nomogramma è stato sottoposto a validazione interna tramite ricampionamento bootstrap (1000 iterazioni) per valutare l'accuratezza predittiva e la discriminazione utilizzando curve di calibrazione (AUC e ROC).
Analisi statistiche
Tutte le analisi statistiche sono state effettuate utilizzando R versione 4.1.2 e IBM SPSS Statistics versione 25.0. I dati categorici erano espressi come frequenze e percentuali; i confronti tra i gruppi sono stati effettuati utilizzando il test chi-quadrato o il test esatto di Fisher, secondo le opportunità. I dati continui con distribuzione normale sono stati espressi come deviazione media ± standard (SD) e confrontati utilizzando il test t dei campioni indipendenti; i dati continui non normalmente distribuiti sono stati espressi come mediana (intervallo interquartile (IQR)) e confrontati utilizzando il test U di Mann-Whitney. Variabili con P < 0,1 nell'analisi di regressione logistica univariata sono state incluse nel modello, coerente con la pratica standard nello sviluppo di modelli predittivi. I pacchetti R utilizzati erano "rms", "mstate", "data.table", "pROC", "rmada", "rio", "boot" e "missForest". La costruzione del nomogramma veniva implementata utilizzando le funzioni lrm e nomogram del pacchetto rms. Le curve ROC e i valori AUC venivano calcolati utilizzando le funzioni roc e auc dal pacchetto pROC. Le curve di calibrazione venivano generate con la funzione calibrate in RMS. L'analisi della curva di decisione (DCA) veniva eseguita utilizzando la funzione della curva decisionale dal pacchetto RMDA. L'imputazione dei dati mancanti veniva eseguita utilizzando la funzione missForest. La validazione interna di bootstrap (1000 iterazioni) veniva effettuata utilizzando il pacchetto di avvio. All'inizio dell'analisi veniva applicato un seed casuale fisso (set.seed[240708] per garantirne la riproducibilità. Un valore P di < 0,05 è stato considerato statisticamente significativo. La formula del modello di regressione logistica era:
glm(mucus_status ~ bronchiectasi + CRS + BMI + FEF25_75 + RV_TLC + VitD, famiglia = "binomiale")