Questo studio retrospettivo trasversale di associazione ha arruolato pazienti consecutivi ricoverati presso il Dipartimento di Neurologia dell'Ospedale Popolare N. 2 di Changzhou per capogiri, cefalea, screening di ictus o valutazione cognitiva tra gennaio 2022 e dicembre 2024. Il consenso informato è stato omesso a causa della natura retrospettiva e non interventistica dello studio e dell'utilizzo di dati completamente anonimizzati. Lo studio è stato approvato dal Comitato Etico per le Tecnologie Mediche Cliniche dell'Ospedale Popolare N. 2 di Changzhou (numero di approvazione IRB: [2023] YLJSA069) ed è stato condotto in conformità con le linee guida istituzionali.
Progettazione dello studio e popolazione di pazienti
Inizialmente, sono stati valutati 850 pazienti. I criteri di esclusione sono stati rigorosamente definiti come segue: (1) mancanza di sequenze complete di risonanza magnetica cranica (in particolare l'assenza di imaging ponderato in suscettibilità [SWI], che impediva una valutazione accurata dei microemorragi); (2) assenza di registrazioni dettagliate dei sondaggi parodontali o dati sul numero completo dei denti in bocca; (3) dati incompleti dell'emocromo, che impedivano il calcolo dell'indice sistemico di infiammazione (SII); (4) storia di ictus massivo, tumore cerebrale, trauma cranico o infezione del sistema nervoso centrale; e (5) presenza di infezione acuta attiva, grave malattia autoimmune o malignità negli ultimi trenta giorni, al fine di minimizzare gli effetti confondenti sui marcatori infiammatori sistemici. Per garantire l'integrità dei dati, è stato adottato un approccio di analisi dei casi completi, secondo il quale sono stati esclusi i soggetti privi di uno qualsiasi dei parametri clinico-radiologici fondamentali. Dopo il controllo, 234 partecipanti idonei sono stati inclusi nell'analisi. Questo campione, comprendente 129 eventi ad alto carico, soddisfaceva la regola empirica degli eventi-per-variazione (EPV) >10 necessaria per una modellizzazione multivariata stabile.
Valutazione dello Stato Parodontale e dell'Infiammazione
Due parodontologi calibrati hanno eseguito tutte le valutazioni cliniche orali utilizzando una sonda parodontale standardizzata.
Prima dello studio, entrambi gli esaminatori hanno seguito una formazione di standardizzazione e l'affidabilità tra esaminatori è risultata eccellente (κ di Cohen = 0,82). La profondità della tasca sondabile (PPD) e il livello di attacco clinico (CAL) sono stati misurati in sei siti per dente. In base al quadro classificatorio del 2018 descritto da Tonetti et al.15, la malattia parodontale è stata classificata come lieve/assente, moderata o grave. In particolare, i partecipanti sono stati suddivisi come segue: (1) Parodontite lieve/assente (inclusa la salute parodontale e la parodontite di stadio I), definita come livello di attacco clinico interdentale (CAL) nel sito con la maggiore perdita di attacco ≤2 mm e profondità della tasca sondabile (PPD) ≤4 mm, senza perdita di denti correlata alla parodontite; (2) Parodontite moderata (parodontite di stadio II), definita come CAL interdentale di 3–4 mm, PPD massima ≤5 mm e non più di quattro denti persi a causa della parodontite; e (3) Parodontite grave (parodontite di stadio III/IV), definita come CAL interdentale ≥5 mm, PPD ≥6 mm e/o perdita di quattro o più denti a causa della distruzione parodontale15.
Il numero di denti mantenuti è stato inoltre documentato per ogni partecipante. Poiché il mantenimento di almeno 20 denti funzionali è un parametro clinico ampiamente accettato per preservare la funzione orale di base e un invecchiamento orale di successo16,17,18, il conteggio grezzo dei denti è stato convertito in una variabile binaria: perdita dentale grave (<20 denti) rispetto a perdita dentale non grave (≥20 denti).
Per quanto riguarda l'infiammazione sistemica, l'SII è stato calcolato a partire da campioni di sangue raccolti a digiuno al momento del ricovero. I campioni di sangue venoso sono stati raccolti in provette con EDTA e processati secondo il protocollo standard del laboratorio clinico dell'ospedale, prima dell'analisi mediante un analizzatore ematologico automatico per ottenere i conteggi di neutrofili, linfociti e piastrine. L'indice è stato calcolato come:

Per tenere conto della distribuzione fortemente asimmetrica dei valori di SII, è stata applicata una trasformazione logaritmica naturale (Log_SII) prima dell'inclusione nelle analisi statistiche.
Acquisizione di Risonanza Magnetica e Carico Totale di CSVD
Tutti i soggetti hanno effettuato esami standardizzati di risonanza magnetica cranica mediante uno scanner RM da 3,0 T dotato di una bobina per la testa standard. Il protocollo di imaging includeva sequenze pesate in T1, pesate in T2, fluid-attenuated inversion recovery (FLAIR) e susceptibility-weighted imaging (SWI), acquisite con uno spessore di fetta di 5 mm. I parametri standard di imaging sono stati configurati come segue: imaging pesato in T1 (tempo di ripetizione [TR] = 2.000 ms, tempo d'eco [TE] = 9 ms, campo visivo [FOV] = 230 × 230 mm2, dimensione della matrice = 256 × 256); imaging pesato in T2 (TR = 4.500 ms, TE = 85 ms, FOV = 230 × 230 mm2, dimensione della matrice = 256 × 256); imaging fluid-attenuated inversion recovery (FLAIR) (TR = 8.500 ms, TE = 120 ms, tempo d'inversione [TI] = 2.400 ms, FOV = 230 × 230 mm2, dimensione della matrice = 256 × 256); e susceptibility-weighted imaging (SWI) (TR = 28 ms, TE = 20 ms, angolo di flip = 15°, FOV = 230 × 230 mm2, dimensione della matrice = 256 × 256). Tutte le sequenze sono state acquisite con uno spessore di fetta di 5 mm e un intervallo interslice di 1,0 mm. Lo scanner specifico è indicato nella Tabella dei Materiali. I marcatori di imaging della CSVD sono stati valutati indipendentemente da due medici radiologi specializzati in neuroimaging, in cecità rispetto ai dati clinici dei partecipanti. L'affidabilità inter-osservatore del punteggio totale di carico di CSVD è risultata elevata (κ di Cohen = 0,85). Eventuali discrepanze sono state risolte per consenso con un terzo medico esperto.
In conformità con il consenso internazionale stabilito1,3, il carico cumulativo di CSVD è stato quantificato su una scala da 0 a 4 valutando quattro caratteristiche di neuroimmagine. È stato assegnato un punto per ciascuno dei seguenti reperti alla risonanza magnetica: (1) almeno un infarto lacunare; (2) uno o più microemorragie cerebrali (CMBs); (3) spazi perivascolari dilatati (EPVS) da moderati a gravi nei gangli della base (grado ≥2); e (4) iperintensità severe della sostanza bianca (WMH), definite come un punteggio Fazekas ≥2 nella sostanza bianca profonda o 3 nella regione periventricolare. Le definizioni diagnostiche hanno rigorosamente seguito il consenso internazionale Standards for Reporting Vascular Changes on Neuroimaging (STRIVE-1)1. Le iperintensità della sostanza bianca (WMH) sono state classificate utilizzando la scala di Fazekas19, e gli spazi perivascolari dilatati (EPVS) nei gangli della base sono stati classificati mediante la scala visiva validata a 4 punti descritta da Potter et al.20. Ai fini dello sviluppo del modello di valutazione del rischio, i pazienti con un punteggio totale di carico di CSVD ≥2 sono stati classificati come aventi un elevato carico di CSVD, poiché questa soglia è stata costantemente associata a un rapido declino cognitivo e alla mortalità21,22,23.
Raccolta dei Dati sulle Covariate
Le caratteristiche demografiche e la storia clinica dei pazienti sono state estratte dai record elettronici sanitari istituzionali. Le variabili raccolte includevano età, sesso, indice di massa corporea (BMI), stato di fumatore e consumo di alcol. Inoltre, sono state documentate per ogni partecipante le comorbidità cardiometaboliche e vascolari, tra cui dislipidemia, diabete mellito, ipertensione, malattia coronarica (CAD), storia di infarto miocardico (MI) e precedente ictus ischemico.
Tutte le comorbidità sono state confermate mediante revisione delle diagnosi effettuate dai medici documentate nelle cartelle cliniche dei pazienti ricoverati, della storia farmacologica e delle valutazioni di laboratorio e di imaging eseguite in fase di ricovero, in conformità con linee guida cliniche stabilite. L'ipertensione è stata definita come pressione arteriosa sistolica ≥140 mmHg, pressione arteriosa diastolica ≥90 mmHg o trattamento attuale con farmaci antipertensivi. Il diabete mellito è stato definito come concentrazione di glucosio plasmatico a digiuno ≥7,0 mmol/L, emoglobina glicata (HbA1c) ≥6,5% o uso di farmaci ipoglicemizzanti. La dislipidemia è stata definita sulla base di alterazioni dei lipidi sierici a digiuno o di un trattamento attuale con farmaci ipolipemizzanti. La malattia coronarica (CAD), la storia di infarto miocardico (MI) e l'ictus ischemico pregresso sono stati confermati mediante la storia clinica documentata e i referti di imaging neurovascolare o cardiaco disponibili nel sistema cartella clinica elettronica.
Analisi Statistica e Sviluppo del Modello
Prima dell'analisi statistica, la normalità delle variabili continue è stata valutata mediante il test di Shapiro–Wilk. Le statistiche descrittive per le variabili categoriche sono state riassunte come frequenze (percentuali) e confrontate utilizzando il test chi-quadrato di Pearson o il test esatto di Fisher, secondo opportunità. Le variabili continue con distribuzione normale sono state espresse come medie ± deviazioni standard e confrontate mediante il test t di Student, mentre le variabili con distribuzione non normale sono state espresse come mediane con intervalli interquartili e confrontate mediante il test U di Mann-Whitney. Per indagare potenziali modelli dose-risposta non lineari tra il numero di denti mantenuti e il rischio di elevato carico di CSVD, è stato costruito un modello RCS aggiustato per età, ipertensione e Log_SII, utilizzando quattro nodi posizionati ai percentili 5°, 35°, 65° e 95° della distribuzione del numero di denti. Un test formale per la non linearità è stato effettuato mediante analisi della varianza (ANOVA) per determinare se un effetto soglia o una tendenza lineare continua spiegasse meglio l'associazione osservata.
Per eseguire la selezione delle variabili in presenza di predittori correlati, tutte le variabili alla baseline sono state inserite in un modello di regressione LASSO. Il parametro di regolazione ottimale (λ) è stato determinato mediante validazione incrociata a 10 ripiegamenti e scelto in base al criterio di devianza binomiale minima (λmin). Sono state mantenute per le analisi successive le variabili con coefficienti diversi da zero. Le variabili selezionate sono state successivamente inserite in un modello di regressione logistica multivariabile per stimare i rapporti di probabilità (OR) e gli intervalli di confidenza al 95% (IC) per le loro associazioni con un elevato carico di CSVD. Il modello multivariabile finale è stato quindi utilizzato per costruire un nomogramma clinico esplorativo per la stima del rischio individualizzata.
Per valutare il valore aggiuntivo dei parametri parodontali, un modello di base (inclusi età, ipertensione e Log_SII) è stato confrontato con un modello esteso che includeva il grado di parodontite e il numero di denti. Le prestazioni del modello sono state valutate utilizzando diverse misure complementari. La capacità discriminante è stata valutata mediante l'area sotto la curva ROC (AUC), e le differenze tra i modelli sono state confrontate utilizzando il test di DeLong24. Al valore soglia ottimale determinato dall'indice di Youden massimo, sono state calcolate sensibilità, specificità, valore predittivo positivo (PPV) e valore predittivo negativo (NPV). Sono stati calcolati il miglioramento continuo della riclassificazione netta (NRI) e il miglioramento integrato della discriminazione (IDI) per quantificare il miglioramento nella riclassificazione del rischio fornito dal modello esteso25. La calibrazione è stata valutata mediante grafici di calibrazione generati da 1.000 campionamenti bootstrap effettuati con il pacchetto rms, e l'errore assoluto medio tra i rischi previsti e osservati è stato calcolato. Per tenere conto di un possibile overfitting e fornire una valutazione imparziale delle prestazioni del modello, la validazione bootstrap ha incluso la correzione per l'ottimismo, riportando sia le prestazioni corrette per il bias sia quelle apparenti del modello. È stata eseguita un'analisi della curva decisionale (DCA) per valutare l'utilità clinica potenziale stimando il beneficio netto in un intervallo di probabilità soglia. Tutte le analisi statistiche e le rappresentazioni grafiche sono state effettuate utilizzando software statistici e i pacchetti rms, glmnet, pROC, PredictABEL e dcurves. La significatività statistica è stata definita come un valore P bilaterale <0,05.