Das Studienprotokoll wurde vom Ethikkomitee des Allgemeinen Krankenhauses des nördlichen Befehlsgebiets (Genehmigungsnummer: Ethik Y (2026) 75) geprüft und genehmigt. Obwohl die klinischen Daten von Patienten stammen, die zwischen Januar 2024 und Dezember 2025 behandelt wurden, wurde die vorliegende Studie als retrospektive Analyse bereits vorhandener klinischer Unterlagen durchgeführt. Die im Jahr 2026 erhaltene ethische Genehmigung umfasste die retrospektive Datenerhebung, Anonymisierung, Analyse und Veröffentlichung dieser zuvor gesammelten klinischen Daten. Es wurde keine prospektive Intervention oder Patientenrekrutierung vor der ethischen Genehmigung durchgeführt. Die Anforderung einer schriftlichen informierten Einwilligung wurde aufgrund des retrospektiven Beobachtungsstudien-Designs aufgehoben. Alle Verfahren wurden gemäß der Deklaration von Helsinki durchgeführt, und alle klinischen Daten wurden vor der Analyse anonymisiert. Die vorliegende Analyse wurde als Teil eines genehmigten retrospektiven Forschungsprojekts zur neurointerventionellen klinischen Datenauswertung durchgeführt. Die für die unten beschriebenen Verfahren benötigten Materialien und Geräte sind in der Tabelle der Materialien zusammengefasst.
Patientenrekrutierung und Bestimmung der CIE
Möglicherweise berechtigte Patienten wurden identifiziert, indem das elektronische Krankenaktensystem nach Patienten durchsucht wurde, die zwischen Januar 2024 und Dezember 2025 in der Abteilung für Neurochirurgie des Allgemeinen Krankenhauses des Nordtheaterkommandos behandelt wurden. Die anfängliche Suche wurde durch Angaben zum Aufnahmedatum, zur Abteilung, zur Diagnose einer zerebrovaskulären Erkrankung und zu neurointerventionellen Verfahren eingeschränkt. Der detaillierte Ablauf zur Patientenrekrutierung und Zuteilung der Kohorte ist in Abbildung 2 dargestellt. Patienten kamen für die Einschlussnahme in Betracht, wenn sie alle folgenden Kriterien erfüllten: (1) Alter zwischen 30 und 80 Jahren; (2) Diagnose einer zerebrovaskulären Erkrankung, die eine neurointerventionelle Beurteilung oder Behandlung erforderte; (3) Vorhandensein vollständiger klinischer Unterlagen und Laborergebnisse; sowie (4) Durchführung einer postoperativen Schädel-CT-Untersuchung innerhalb von 3 Tagen nach der Operation.
Patienten wurden ausgeschlossen, wenn sie unvollständige klinische oder bildgebende Daten aufwiesen, keine endovaskuläre Intervention erhielten, schwere präoperative neurologische Defizite (modifizierter Rankin-Skala-Score ≥4) hatten, eine akute zerebrale Infarktion im präoperativen diffusionsgewichteten MRT zeigten, eine chronische Nierenerkrankung im Stadium G4 oder höher gemäß den Kriterien der Kidney Disease: Improving Global Outcomes (KDIGO) aufwiesen oder eine schwere kardiovaskuläre Erkrankung hatten. Die Einschluss- und Ausschlusskriterien wurden nacheinander angewendet. Nach der initialen elektronischen Suche wurde jeder Kandidatendatensatz manuell überprüft, um die Eignung zu bestätigen und gegebenenfalls den Ausschlussgrund zu dokumentieren. Nach der Selektion wurden die in Frage kommenden Patienten im Verhältnis 8:2 zufällig den Trainings- und Testkohorten zugewiesen, bestehend aus einer Trainingskohorte (n = 128) und einer Testkohorte (n = 33). Es wurde eine geschichtete Zufallsstichprobe verwendet, um eine gleichmäßige Verteilung von CIE- und Nicht-CIE-Fällen zwischen den beiden Kohorten sicherzustellen.
Die Diagnose einer CIE basierte hauptsächlich auf dem zeitlichen Zusammenhang zwischen dem Auftreten der Symptome und zerebrovaskulären interventionsmedizinischen Verfahren, kombiniert mit einer ausschließenden bildgebenden Beurteilung. Zwei unabhängige Spezialisten für zerebrovaskuläre Interventionen überprüften die postoperativen klinischen Befunde, den Zeitpunkt des Symptombeginns, neurologische Manifestationen sowie die postoperativen bildgebenden Ergebnisse bei allen Verdachtsfällen. Unstimmigkeiten wurden durch Diskussion bis zur Einigung behoben. Patienten wurden als kompatibel mit einer CIE angesehen, wenn innerhalb von 24 Stunden nach Abschluss des neurointerventionellen Eingriffs neu aufgetretenes kortikales Erblinden, Bewusstseinsstörungen (einschließlich Somnolenz, Delirium oder Koma), Krampfanfälle oder fokale neurologische Defizite wie Hemiparese oder Aphasie auftraten.
Nach dem Auftreten der Symptome wurden alternative Ursachen, die ähnliche neurologische Manifestationen hervorrufen können, durch die Überprüfung von postoperativen CT- oder MRT-Bildern zusammen mit dem klinischen Verlauf ausgeschlossen. Die postoperative Schädel-CT wurde gemäß dem institutionellen Standardprotokoll für nicht kontrastmittelverstärkte kraniale CT durchgeführt. Zu den wesentlichen Aufnahmeparametern gehörten eine Röhrenspannung von 80 kVp, ein Röhrenstrom von 260 mA, eine Schichtdicke von 0,5 mm und eine standardmäßige axiale Rekonstruktion. Die Bilder wurden von erfahrenen Ärzten in Hirn- und Knochenfenstern analysiert. Bei der Differenzialdiagnose wurden akute intrakranielle Blutung, subarachnoidale Blutung, intrazerebrale Blutung, neu aufgetretene großflächige zerebrale Infarkte, anfallsassoziierte Veränderungen, Infektionen und metabolische Störungen berücksichtigt. Postsymptomatische CT-Untersuchungen zeigten typischerweise fokale oder diffuse zerebrale Ödeme sowie Schatten mit hoher Dichte, die die Großhirnrinde, subkortikale Bereiche oder den subarachnoidalen Raum betrafen und hämorrhagische Läsionen nachahmten. MRT-Befunde zeigen gewöhnlich eine Schwellung der Großhirnrinde, insbesondere im Bereich der Temporo-Parieto-Okzipital-Kortex. Bei einigen Patienten zeigten CT- oder MRT-Untersuchungen jedoch keine offensichtlichen Auffälligkeiten, obwohl die klinischen Symptome entsprechend waren.
Erhebung perioperativer Variablen und Erstellung des CGR
Allgemeine demografische Daten, einschließlich Alter, Geschlecht und Körpergewicht, wurden aus strukturierten Feldern des elektronischen Krankenakten-Systems extrahiert. Variablen zur klinischen Anamnese, darunter Rauchgeschichte, Alkoholkonsumgeschichte, Hypertonie, Diabetes mellitus und koronare Herzkrankheit, wurden aus Eintrittsnotizen, vorbestehenden Krankengeschichten und Entlassungsdiagnosen entnommen und auf Konsistenz überprüft. Laborwerte wurden aus dem Laborinformationssystem extrahiert. Für jeden Patienten wurde die erste venöse Blutprobe verwendet, die nach der Aufnahme und vor dem neurointerventionellen Eingriff entnommen wurde. Zu den extrahierten Laborvariablen gehörten Serumkreatinin, die geschätzte glomeruläre Filtrationsrate (eGFR), Gesamtcholesterin und Triglyceridspiegel. Die eGFR wurde anhand der Gleichung der Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI), die auf Kreatinin basiert, berechnet und dem Laborinformationssystem des Krankenhauses entnommen.
Da die Verfahrenskomplexität, die Exposition gegenüber Kontrastmittel und die Läsionslage die Blut-Hirn-Schranke und die Retention von Kontrastmittel während neurointerventioneller Eingriffe beeinflussen können, wurden auch perioperative Verfahrensvariablen in das Vorhersagemodell integriert. Zu den verfahrensbezogenen Variablen gehörten die Läsionslage, der Eingriffstyp, die Dauer des Eingriffs, der Typ des Kontrastmittels und das gesamte Kontrastmittelvolumen. Die Läsionslage und der Eingriffstyp wurden anhand der Operationsberichte und angiographischen Aufzeichnungen bestimmt. Die Dauer des Eingriffs wurde als Zeit von der arteriellen Punktion bis zum Abschluss des neurointerventionellen Eingriffs definiert. Das gesamte Kontrastmittelvolumen wurde als kumulatives Volumen des verabreichten jodhaltigen Kontrastmittels vom Beginn bis zum Ende des Eingriffs definiert. Die in dieser Studie verwendeten Kontrastmittel waren nichtionische IKM, darunter Iodixanol (320 mg Jod/mL), die gemäß der klinischen Praxis des Instituts verabreicht wurden. Alle in diese Studie eingeschlossenen neurointerventionellen Eingriffe wurden unter digitaler Subtraktionsangiographie-Begleitung durch ein erfahrenes neurointerventionelles Team am Studienzentrum durchgeführt. Alle Eingriffe wurden über den transfemoralen arteriellen Zugang vorgenommen. Während des Eingriffs erhielten die Patienten eine intraoperative Antikoagulation mit Heparin, und die aktivierte Gerinnungszeit wurde im Zielbereich von 250–300 s aufrechterhalten. Die Vitalzeichen wurden während des gesamten Eingriffs kontinuierlich überwacht.
Therapeutische Maßnahmen, einschließlich der Embolisierung von Aneurysmen, Angioplastie und Stentplatzierung, wurden entsprechend den Läsionscharakteristika und dem Ermessen des Operateurs durchgeführt. Entsprechende Mikrokatheter, Coils, Stents oder Ballons wurden je nach den Anforderungen des Eingriffs ausgewählt. Während der angiographischen Bildgebung und Intervention wurden Kontrastmittel entweder mittels Hochdruckinjektion oder manueller Infusion mit konstanter Rate verabreicht. Art des Kontrastmittels und Gesamtkontrastmittelvolumen wurden aus den Prozedurunterlagen entnommen und, sofern verfügbar, mit den Aufzeichnungen der Anästhesie oder Pflegepersonal abgeglichen. Unterlagen mit widersprüchlichen oder unvollständigen Angaben zum Kontrastmittelvolumen wurden vor der Aufnahme in den endgültigen analytischen Datensatz manuell überprüft.
Alle verfahrensbezogenen Variablen, einschließlich Dauer des Eingriffs, Art der therapeutischen Intervention, Nutzung von Geräten und Aufzeichnungen zur Kontrastmittelverabreichung, wurden systematisch überprüft und hinsichtlich Vollständigkeit und Konsistenz vor der Analyse verifiziert. Um die Diskrepanz zwischen Kontrastmittelbelastung und renaler Ausscheidungskapazität quantitativ zu bewerten, wurde das CGR als zentrale prädiktive Variable festgelegt und gemäß der folgenden Gleichung berechnet:
(1)
Das gesamte Kontrastmittelvolumen wurde in mL erfasst und das eGFR in mL/min/1,73 m2. Der für die CGR-Berechnung verwendete eGFR-Wert war der präoperative eGFR-Wert, der aus der ersten venösen Blutprobe nach der Aufnahme und vor dem neurointerventionellen Eingriff gewonnen wurde. Für jeden Patienten wurde das CGR berechnet, nachdem das gesamte Kontrastmittelvolumen und das eGFR überprüft worden waren. Der Zähler entsprach dem gesamten Kontrastmittelvolumen in mL, der Nenner dem eGFR in mL/min/1,73 m2. Dieselbe Berechnungsmethode wurde auf alle Patienten vor der Modellentwicklung angewandt.
Merkmalsauswahl und ML-Workflow
Um die analytische Robustheit sicherzustellen, wurden alle klinischen und laborbezogenen Variablen vor der Analyse systematischen Qualitätskontrollverfahren unterzogen. Der Anteil fehlender Daten wurde für jede Kandidatenvariable vor der Modellerstellung berechnet. Im endgültigen analytischen Datensatz wurden keine fehlenden Werte bei den eingeschlossenen Prädiktoren oder Ergebnisvariablen beobachtet; daher wurde keine Variable aufgrund fehlender Werte ausgeschlossen, und eine multiple Imputation war nicht erforderlich. Die detaillierte Zusammenfassung der fehlenden Daten ist in Supplementary Table 1 enthalten.
Alle potenziellen Prädiktorvariablen wurden zunächst mittels univariater Analysen hinsichtlich ihres Zusammenhangs mit dem Auftreten von CIE bewertet. Die statistischen Methoden wurden entsprechend der Verteilungsmerkmale der Daten ausgewählt. Variablen mit normaler Verteilung wurden als Mittelwert ± Standardabweichung angegeben und mit dem t-Test für unabhängige Stichproben verglichen. Variablen mit nichtnormaler Verteilung wurden als Median (interquartiler Bereich, IQR) dargestellt und mittels des Mann-Whitney-U-Tests analysiert. Kategoriale Variablen wurden als Häufigkeiten und Prozentanteile angegeben und mit dem Chi-Quadrat-Test oder dem exakten Fisher-Test analysiert. Eine statistische Signifikanz wurde bei p < 0,05 festgelegt.
Der gesamte Datensatz wurde vor der Modellentwicklung zufällig im Verhältnis 8:2 in Trainings- und Testkohorten aufgeteilt. Um Überanpassung zu minimieren und Multikollinearität zu reduzieren, wurde eine Regression mit dem Verfahren der kleinsten absoluten Schrumpfung und Auswahl (LASSO) in Kombination mit einer 10-fachen Kreuzvalidierung zur Merkmalsauswahl angewandt. Die Merkmalsauswahl, die Optimierung der Hyperparameter sowie die Modellentwicklung erfolgten ausschließlich anhand der Trainingskohorte, um Datenlecks zu vermeiden. Die interne Testkohorte wurde während der Imputationsparameter-Schätzung, Merkmalsauswahl, Hyperparameter-Optimierung oder Modelltrainierung nicht verwendet und erst einmalig zur abschließenden Leistungsbewertung herangezogen. Die Optimierung der Hyperparameter wurde mittels einer Gitter-Suchstrategie mit 10-facher Kreuzvalidierung innerhalb der Trainingskohorte durchgeführt. Konkret wurde die Trainingskohorte zufällig in 10 sich gegenseitig ausschließende Teilmengen unterteilt. In jeder Iteration dienten neun Teilmengen zum Modelltraining, während die verbleibende Teilmenge zur Validierung genutzt wurde. Dieser Vorgang wurde 10-mal wiederholt, um sicherzustellen, dass jede Teilmenge genau einmal als Validierungskohorte fungierte.
Alle computergestützten Analysen wurden in einer Python-Umgebung durchgeführt. Die Entwicklung und Bewertung der ML-Modelle erfolgte mit scikit-learn. Die Analyse-Skripte für die Datenvorverarbeitung, Merkmalsauswahl, Modelltraining, Leistungsbewertung und DCA wurden als ergänzende Code-Dateien bereitgestellt. Für die Modellimplementierung wurden die aus der LASSO-Regression ausgewählten Merkmale als Eingabevariablen für alle in Frage kommenden Klassifikatoren verwendet. Stetige Variablen wurden mittels z-Score-Normalisierung innerhalb des Trainingskollektivs standardisiert, und dieselben Skalierungsparameter wurden auf das interne Testkollektiv angewandt. Kategoriale Variablen wurden gemäß dem vordefinierten Codierungsschema kodiert; die Variablendefinitionen und das Codierungsschema sind in Supplementary Table 2 angegeben. Fünf Klassifikatoren wurden implementiert: Gaussian Naive Bayes, SVM, KNN, LightGBM und MLP. Bei SVM wurden der Kernel-Typ, der Regularisierungsparameter C und der Kernel-Koeffizient gamma optimiert. Bei KNN wurden die Anzahl der Nachbarn und die Distanzgewichtungsstrategie angepasst. Bei LightGBM wurden die Anzahl der Schätzer, die Lernrate, die maximale Baumtiefe und die Anzahl der Blätter optimiert. Bei MLP wurden die Struktur der versteckten Schichten, die Aktivierungsfunktion, der Regularisierungsparameter und die Lernrate angepasst. Gaussian Naive Bayes wurde mit der vorgegebenen Varianzglättungseinstellung implementiert. Die optimale Kombination von Hyperparametern für jedes Modell wurde basierend auf der Kreuzvalidationsleistung innerhalb des Trainingskollektivs ausgewählt. Die endgültig ausgewählten optimierten Hyperparameterkombinationen für jedes Modell sind in Supplementary Table 3 zusammengefasst. Nach der Hyperparameterauswahl wurde jedes finale Modell erneut auf dem gesamten Trainingskollektiv angepasst und anschließend einmalig im internen Testkollektiv evaluiert.
Vergleiche der Modellleistung
Basierend auf den ausgewählten Prädiktoren wurden fünf ML-Modelle entwickelt und miteinander verglichen. Die Modellleistung wurde sowohl in der Trainings- als auch in der Testkohorte anhand der Fläche unter der Receiver-Operating-Characteristic-Kurve (AUC), des 95 % Konfidenzintervalls (95 % CI), der Sensitivität, Spezifität und Genauigkeit bewertet. Receiver-Operating-Characteristic-Kurven wurden erstellt, um die Unterscheidungsleistung der verschiedenen Modelle zu vergleichen. Der optimale Klassifizierungsschwellenwert für jedes Modell wurde ausschließlich anhand der Trainingskohorte mithilfe des Youden-Index bestimmt. Der ermittelte Schwellenwert wurde anschließend fixiert und unverändert auf die interne Testkohorte angewandt, um schwellenwertabhängige Leistungsmetriken zu berechnen.
Decision-Curve-Analyse
Die DCA wurde durchgeführt, indem der Nettovorteil jedes Modells über einen Bereich von Schwellenwahrscheinlichkeiten berechnet und die modellbasierte Strategie mit zwei Standardstrategien verglichen wurde: Behandlung aller Patienten und Behandlung keiner Patienten. Die DCA wurde gewählt, da sie eine gleichzeitige Beurteilung der Diskriminierungsleistung und der potenziellen klinischen Anwendbarkeit in Entscheidungssituationen im perioperativen Bereich ermöglicht. Zur Verbesserung der Reproduzierbarkeit wurde der gesamte Arbeitsablauf in der folgenden Reihenfolge durchgeführt: Identifizierung der Patienten aus dem elektronischen Krankenakten-System, Anwendung der Einschluss- und Ausschlusskriterien, Adjudikation von CIE durch zwei unabhängige Fachärzte, Extraktion und Gegenprüfung von demografischen, anamnestischen, laborchemischen und verfahrensbezogenen Variablen, Berechnung des CGR, Bewertung fehlender Daten sowie Aufteilung in Trainings- und Testkohorte, Merkmalsauswahl basierend auf LASSO innerhalb der Trainingskohorte, Optimierung der Hyperparameter mittels 10-facher Kreuzvalidierung, erneutes Anpassen des endgültigen Modells in der gesamten Trainingskohorte, interne Evaluierung anhand der Testkohorte, leistungsbezogene Bewertung basierend auf der ROC-Kurve, Bestimmung des Schwellenwerts nach dem Youden-Index und DCA.