Diese retrospektive, querschnittsanalytische Assoziationsstudie schloss aufeinanderfolgende Patienten ein, die zwischen Januar 2022 und Dezember 2024 zur Abklärung von Schwindel, Kopfschmerzen, Schlaganfallvorsorge oder kognitiver Beurteilung in die Abteilung für Neurologie des Zweiten Volkskrankenhauses von Changzhou aufgenommen wurden. Die informierte Einwilligung wurde aufgrund des retrospektiven, nicht-interventionellen Charakters der Studie sowie der Verwendung vollständig anonymisierter Daten nicht eingeholt. Die Studie wurde vom Ethikkomitee für klinische medizinische Technologien des Zweiten Volkskrankenhauses von Changzhou genehmigt (IRB-Zulassungs-ID: [2023] YLJSA069) und gemäß den institutionellen Richtlinien durchgeführt.
Studienaufbau und Patientenpopulation
Zunächst wurden 850 Patienten evaluiert. Die Ausschlusskriterien wurden streng wie folgt definiert: (1) unvollständige Schädel-MRT-Sequenzen (insbesondere fehlende susceptibility-gewichtete Bildgebung [SWI], was eine genaue Beurteilung von Mikroblutungen verhinderte); (2) fehlende detaillierte parodontale Sondierungsbefunde oder Angaben zur Vollgebiss-Zahlanzahl; (3) unvollständige Blutbildbefunde, die eine Berechnung des SII verhinderten; (4) Anamnese eines massiven Schlaganfalls, eines Hirntumors, einer traumatischen Hirnverletzung oder einer Infektion des zentralen Nervensystems; sowie (5) Vorliegen einer aktiven akuten Infektion, einer schweren Autoimmunerkrankung oder einer Malignität innerhalb des letzten Monats, um störende Einflüsse auf systemische inflammatorische Marker zu minimieren. Zur Sicherstellung der Datenintegrität wurde ein Complete-Case-Analyseansatz verwendet, bei dem Personen, bei denen irgendein zentrales klinisch-radiologisches Parameter fehlte, ausgeschlossen wurden. Nach der Selektion wurden 234 berechtigte Teilnehmer in die Analyse eingeschlossen. Diese Stichprobe, bestehend aus 129 Hochlast-Ereignissen, erfüllte das Verhältnis von Ereignissen pro Variable (EPV) >10 Faustregeln für stabile multivariable Modellierung.
Beurteilung des parodontalen Status und der Entzündung
Zwei kalibrierte Parodontologen führten alle klinischen oralen Beurteilungen mit einer standardisierten Parodontaltaschenmesssonde durch.
Vor der Studie absolvierten beide Untersucher eine Standardisierungsschulung, und die Übereinstimmung zwischen den Untersuchern war ausgezeichnet (Cohens κ = 0,82). Die Taschentiefe bei der Befundung (PPD) und das klinische Attachmentniveau (CAL) wurden an sechs Stellen pro Zahn gemessen. Basierend auf dem von Tonetti et al.15 beschriebenen Klassifikationsrahmen aus dem Jahr 2018 wurde die Parodontalerkrankung als leicht/unauffällig, moderat oder schwer eingestuft. Konkret wurden die Teilnehmer wie folgt kategorisiert: (1) Leichte/unauffällige Parodontitis (einschließlich parodontaler Gesundheit und Parodontitis Stadium I), definiert als interdentales klinisches Attachmentniveau (CAL) an der Stelle des größten Attachmentverlusts ≤2 mm und Befundungstaschentiefe (PPD) ≤4 mm, ohne zahnverlustbedingte Folgen der Parodontitis; (2) Moderate Parodontitis (Parodontitis Stadium II), definiert als interdentales CAL von 3–4 mm, maximale PPD ≤5 mm und nicht mehr als vier Zähne verloren aufgrund der Parodontitis; und (3) Schwere Parodontitis (Parodontitis Stadium III/IV), definiert als interdentales CAL ≥5 mm, PPD ≥6 mm und/oder Verlust von vier oder mehr Zähnen aufgrund parodontaler Zerstörung15.
Die Anzahl der erhaltenen Zähne wurde ebenfalls für jeden Teilnehmer dokumentiert. Da das Vorhandensein von mindestens 20 funktionellen Zähnen als weit verbrecheter klinischer Maßstab für die Erhaltung der grundlegenden Mundfunktion und eines erfolgreichen oralen Alterns gilt16,17,18, wurde die absolute Zahl der Zähne in eine binäre Variable umgewandelt: schwerer Zahnverlust (<20 Zähne) versus nicht-schwerer Zahnverlust (≥20 Zähne).
Bezüglich der systemischen Entzündung wurde der SII aus nüchternen Blutproben berechnet, die bei der Aufnahme gesammelt wurden. Venöse Blutproben wurden in EDTA-Röhrchen gesammelt und gemäß dem standardmäßigen klinischen Laborprotokoll des Krankenhauses vor der Analyse mit einem automatisierten Hämatologieanalysator verarbeitet, um die Anzahlen von Neutrophilen, Lymphozyten und Thrombozyten zu ermitteln. Der Index wurde wie folgt berechnet:

Um die stark schief verteilten SII-Werte zu berücksichtigen, wurde vor der Einbeziehung in die statistischen Analysen eine natürliche logarithmische Transformation (Log_SII) angewendet.
MRT-Akquise und gesamte CSVD-Belastung
Alle Probanden unterzogen sich standardisierten zerebralen MRT-Untersuchungen mit einem 3,0-T-MRT-Scanner, der mit einer Standard-Kopfspule ausgestattet war. Das Bildgebungsprotokoll umfasste T1-gewichtete, T2-gewichtete, fluid-attenuierte Inversionsrecovery-(FLAIR)- und suszeptibilitätsgewichtete Bildgebungssequenzen (SWI) mit einer Schichtdicke von 5 mm. Die standardmäßigen Bildgebungsparameter wurden wie folgt festgelegt: T1-gewichtete Bildgebung (Wiederholungszeit [TR] = 2.000 ms, Echozeit [TE] = 9 ms, Bildfeld [FOV] = 230 × 230 mm2, Matrixgröße = 256 × 256); T2-gewichtete Bildgebung (TR = 4.500 ms, TE = 85 ms, FOV = 230 × 230 mm2, Matrixgröße = 256 × 256); FLAIR-Bildgebung (TR = 8.500 ms, TE = 120 ms, Inversionszeit [TI] = 2.400 ms, FOV = 230 × 230 mm2, Matrixgröße = 256 × 256); und SWI-Bildgebung (TR = 28 ms, TE = 20 ms, Kippwinkel = 15°, FOV = 230 × 230 mm2, Matrixgröße = 256 × 256). Alle Sequenzen wurden mit einer Schichtdicke von 5 mm und einem zwischen den Schichten liegenden Abstand von 1,0 mm aufgenommen. Der spezifische Scanner ist in der Materialtabelle aufgeführt. Die bildgebenden Marker der CSVD wurden unabhängig voneinander von zwei auf Neuroradiologie spezialisierten Ärzten bewertet, die keine Kenntnis von den klinischen Daten der Teilnehmer hatten. Die Übereinstimmung zwischen den Beurteilern hinsichtlich der Gesamtbelastungsskala für CSVD war hoch (Cohens κ = 0,85). Unstimmigkeiten wurden im Konsens mit einem dritten, erfahrenen Arzt gelöst.
Entsprechend dem etablierten internationalen Konsensus1,3 wurde die kumulative CSVD-Belastung anhand von vier neurobildgebenden Merkmalen auf einer Skala von 0 bis 4 quantifiziert. Für jeden der folgenden Befunde im MRT erhielt der Patient einen Punkt: (1) mindestens ein lakunärer Infarkt; (2) eine oder mehrere zerebrale Mikroblutungen (CMBs); (3) mäßig bis schwergradig erweiterte perivaskuläre Räume (EPVS) im Bereich der Basalganglien (Grad ≥2); und (4) schwere weiße Substanz-Hyperintensitäten (WMH), definiert als Fazekas-Score von ≥2 in der tiefen weißen Substanz oder 3 im periventrikulären Bereich. Die diagnostischen Definitionen folgten strikt dem internationalen Konsensus „Standards for Reporting Vascular Changes on Neuroimaging“ (STRIVE-1)1. Die weiße Substanz-Hyperintensitäten (WMH) wurden anhand der Fazekas-Skala19 eingestuft, und die erweiterten perivaskulären Räume (EPVS) in den Basalganglien wurden mit der validierten 4-Punkte-Visuell-Bewertungsskala nach Potter et al.20 klassifiziert. Zur Entwicklung des Risikoabschätzungsmodells wurden Patienten mit einem Gesamt-CSVD-Belastungsscore von ≥2 als hoch belastet eingestuft, da diese Schwelle konsistent mit einer beschleunigten kognitiven Verschlechterung und erhöhter Mortalität assoziiert ist21,22,23.
Erhebung von Kovariablen-Daten
Die demografischen Merkmale der Patienten und die klinische Anamnese wurden aus den elektronischen Gesundheitsakten der Institution extrahiert. Zu den erfassten Variablen gehörten Alter, Geschlecht, Body-Mass-Index (BMI), Raucherstatus und Alkoholkonsum. Zusätzlich wurden kardiometabolische und vaskuläre Komorbiditäten, einschließlich Dyslipidämie, Diabetes mellitus, Hypertonie, koronare Herzkrankheit (CAD), Anamnese eines Myokardinfarkts (MI) und vorangegangener ischämischer Schlaganfall, für jeden Teilnehmer dokumentiert.
Alle Komorbiditäten wurden durch Überprüfung der ärztlichen Diagnosen in den stationären Krankenakten, der Medikationsanamnese sowie routinemäßiger Labor- und bildgebender Untersuchungen bei Aufnahme gemäß etablierter klinischer Leitlinien bestätigt. Hypertonie wurde definiert als systolischer Blutdruck ≥140 mmHg, diastolischer Blutdruck ≥90 mmHg oder aktuelle antihypertensive Therapie. Diabetes mellitus wurde definiert als nüchternen Plasma-Glukosekonzentration ≥7,0 mmol/L, glykiertes Hämoglobin (HbA1c) ≥6,5 % oder Einnahme von blutzuckersenkenden Medikamenten. Dyslipidämie wurde anhand nüchterner serumlipider Auffälligkeiten oder aktueller lipid-senkender Therapie festgestellt. KHK, Vorgeschichte eines Myokardinfarkts und vorangegangener ischämischer Schlaganfall wurden anhand dokumentierter klinischer Anamnese und historischer neurovaskulärer oder kardiologischer Bildgebungsbefunde, die im elektronischen Gesundheitsakten-System verfügbar waren, bestätigt.
Statistische Analyse und Modellentwicklung
Vor der statistischen Analyse wurde die Normalverteilung der stetigen Variablen mithilfe des Shapiro-Wilk-Tests überprüft. Deskriptive Statistiken für kategoriale Variablen wurden als Anzahlen (Prozentanteile) zusammengefasst und, soweit angemessen, mittels Chi-Quadrat-Test nach Pearson oder exaktem Test nach Fisher verglichen. Stetige Variablen mit normaler Verteilung wurden als Mittelwerte ± Standardabweichungen angegeben und mittels Students t-Test verglichen, während nicht normalverteilte Variablen als Mediane mit Interquartilsabständen dargestellt und mittels Mann-Whitney-U-Test verglichen wurden. Um mögliche nichtlineare Dosis-Wirkungs-Muster zwischen der Anzahl der erhaltenen Zähne und dem Risiko einer hohen CSVD-Belastung zu untersuchen, wurde ein RCS-Modell, das für Alter, Hypertonie und Log_SII adjustiert war, mit vier Knoten an den 5., 35., 65. und 95. Perzentilen der Zahnanzahlverteilung erstellt. Eine formelle Prüfung auf Nichtlinearität erfolgte mittels Varianzanalyse (ANOVA), um zu bestimmen, ob ein Schwelleneffekt oder ein kontinuierlicher linearer Trend den beobachteten Zusammenhang am besten erklärt.
Um eine Variablenselektion bei Vorliegen korrelierter Prädiktoren durchzuführen, wurden alle Basisvariablen in ein LASSO-Regressionsmodell eingegeben. Der optimale Tuning-Parameter (λ) wurde mittels 10-facher Kreuzvalidierung bestimmt und gemäß dem Kriterium der minimalen binomialen Abweichung ausgewählt (λmin). Variablen mit von null verschiedenen Koeffizienten wurden für nachfolgende Analysen beibehalten. Die ausgewählten Variablen wurden anschließend in ein multivariates logistisches Regressionsmodell eingegeben, um Odds Ratios (ORs) und 95%-Konfidenzintervalle (CIs) bezüglich ihres Zusammenhangs mit einer hohen CSVD-Belastung zu schätzen. Anschließend wurde das endgültige multivariable Modell verwendet, um ein exploratives klinisches Nomogramm zur individuellen Risikoschätzung zu erstellen.
Um den zusätzlichen Nutzen der parodontalen Parameter zu bewerten, wurde ein Basismodell (bestehend aus Alter, Hypertonie und Log_SII) mit einem erweiterten Modell, das den Parodontitisgrad und die Anzahl der Zähne einschloss, verglichen. Die Modellleistung wurde anhand mehrerer komplementärer Maßzahlen beurteilt. Die Diskriminationsfähigkeit wurde mittels der Fläche unter der Receiver-Operating-Characteristic-Kurve (AUC) ermittelt, und die Unterschiede zwischen den Modellen wurden mit dem DeLong-Test24 verglichen. Bei der optimalen Grenze, bestimmt durch den maximalen Youden-Index, wurden zusätzlich Sensitivität, Spezifität, positiver Vorhersagewert (PPV) und negativer Vorhersagewert (NPV) berechnet. Die kontinuierliche netto-Reklassifizierungsverbesserung (NRI) und die integrierte Diskriminierungsverbesserung (IDI) wurden berechnet, um die zusätzliche Risikoreklassifizierung durch das erweiterte Modell zu quantifizieren25. Die Kalibrierung wurde anhand von Kalibrierungsdiagrammen beurteilt, die aus 1.000 Bootstrap-Stichproben mithilfe des rms-Pakets erstellt wurden, und der mittlere absolute Fehler zwischen den vorhergesagten und beobachteten Risiken wurde berechnet. Um eine mögliche Überanpassung zu berücksichtigen und eine unverzerrte Bewertung der Modellleistung zu gewährleisten, wurde bei der Bootstrap-Validierung eine Korrektur für Optimismus vorgenommen, wobei sowohl die bias-korrigierte als auch die scheinbare Modellleistung angegeben wurden. Eine Entscheidungskurvenanalyse (DCA) wurde durchgeführt, um den potenziellen klinischen Nutzen zu bewerten, indem der Nettogewinn über einen Bereich von Schwellenwahrscheinlichkeiten hinweg geschätzt wurde. Alle statistischen Analysen und Visualisierungen wurden mit statistischer Software sowie den Paketen rms, glmnet, pROC, PredictABEL und dcurves durchgeführt. Statistische Signifikanz wurde als zweiseitiger P-Wert von <0,05 definiert.