Diese Studie an menschlichen Probanden, klinischen Aufzeichnungen und Serumperlen wurde gemäß der Deklaration von Helsinki und den institutionellen Anforderungen für die Forschung am Menschen durchgeführt13. Das Studienprotokoll wurde von der Jiangsu Province Hospital of Chinese Medicine unter der Genehmigungsnummer CJ28J12-45 geprüft und genehmigt. Für retrospektiv erhobene anonymisierte klinische Daten und Restserumproben wurde die Anforderung einer zusätzlichen schriftlichen informierten Einwilligung gemäß den ethischen Richtlinien der Institution aufgehoben. Für prospektiv gesammelte Proben wurde vor der Einbindung eine schriftliche informierte Einwilligung eingeholt. Alle Teilnehmerinformationen und Serumproben wurden vor der Analyse anonymisiert. Studienkennungen dienten der Datenverknüpfung, und keine personenbezogenen Identifikationsmerkmale waren im Analyse-Datensatz enthalten.
Studiendesign und Kohortenzuordnung
Diese zweizentrische ambispektive Beobachtungsstudie untersuchte den diagnostischen und prognostischen Wert eines kombinierten Serum-Biomarker-Panels, das das karzinoembryonale Antigen (CEA), das Plattenepithelkarzinom-Antigen (SCC-Ag) und das Kohlenhydrat-Antigen 125 (CA125) für das ösophageale Plattenepithelkarzinom (ESCC) umfasst. Retrospektive Daten wurden aus geeigneten klinischen Befunden und verbliebenen Serumproben vor der Behandlung gewonnen, und eine prospektive Rekrutierung diente der Sammlung zusätzlicher Serumproben vor der Behandlung sowie der Durchführung der Nachbeobachtung im selben Studienrahmen. Die Teilnehmer wurden von einem tertiären Krankenhaus und einem kommunalen Gesundheitszentrum zwischen Januar 2021 und Dezember 2024 rekrutiert, wobei die Nachbeobachtung bis zum 31. Dezember 2025 fortgesetzt wurde. Zentrum A war das tertiäre Krankenhaus und diente als Trainingskohorte für die Modellentwicklung, Schätzung der Koeffizienten, Auswahl der Cutoff-Werte und interne Validierung. Zentrum B war das kommunale Gesundheitszentrum und fungierte als externe Validierungskohorte. Die Studie wurde gemäß den Richtlinien für die transparente Berichterstattung über multivariable Vorhersagemodelle für individuelle Prognose oder Diagnose (TRIPOD) berichtet14.
Insgesamt wurden 277 Teilnehmer eingeschlossen, darunter 194 aus Zentrum A und 83 aus Zentrum B. Der Trainingskohorte gehörten 114 Patienten mit Ösophaguskarzinom (ESCC), 40 Kontrollen mit benignen ösophagealen Erkrankungen und 40 gesunde Kontrollen an. Die externe Validierungskohorte umfasste 53 Patienten mit Ösophaguskarzinom (ESCC), 15 Kontrollen mit benignen ösophagealen Erkrankungen und 15 gesunde Kontrollen. Zur Minimierung von Informationslecks und zur Prüfung der Modellstabilität in einem unabhängigen klinischen Umfeld erfolgte die zentrumsbasierte Zuordnung anstelle einer zufälligen Aufteilung. In der diagnostischen Analyse wurden die Serumspiegel von CEA, SCC-Ag und CA125 bei Patienten mit Ösophaguskarzinom (ESCC), bei Patienten mit benignen ösophagealen Erkrankungen sowie bei gesunden Kontrollen verglichen. Die prognostische Analyse beschränkte sich auf Patienten mit Ösophaguskarzinom (ESCC) und untersuchte, ob die kombinierte Belastung durch Biomarker die Risikostratifizierung für das Gesamtüberleben (OS) und das progressionsfreie Überleben (PFS) verbesserte. Das Studiendesign, die Teilnehmerrekrutierung, die Kohortenzuordnung, der Workflow zur Analyse der serumbasierten Biomarker, die Entwicklung des diagnostischen Modells, die externe Validierung sowie der Vergleich der prognostischen Modelle sind in Abbildung 1 zusammengefasst.
Der Teilnehmerfluss umfasste die Einschreibung, Ausschlüsse, die diagnostische Gruppierung und die zentrumsbasierte Zuordnung. Insgesamt wurden 326 Personen untersucht. 49 Personen wurden ausgeschlossen, darunter 12 mit einer antitumoralen Therapie vor der Blutentnahme, 8 mit einem weiteren primären Malignom, 9 mit schwerer Infektion, Autoimmunerkrankung oder unkontrollierten Komorbiditäten, die wahrscheinlich die serumgebundenen Tumormarker beeinflusst hätten, 11 mit Hämolyse, Lipämie, ikterischer Interferenz, unzureichendem Serumvolumen oder verlängertem Verarbeitungsintervall sowie 9 mit fehlenden wichtigen klinischen, pathologischen oder Nachsorge-Daten. Die endgültige analytische Studienpopulation umfasste 277 Teilnehmer: 167 ESCC-Patienten, 55 Kontrollpersonen mit benignen Speiseröhrenerkrankungen und 55 gesunde Kontrollpersonen.
Begründung zur Stichprobengröße
Die Stichprobengröße wurde anhand der konsekutiven Verfügbarkeit geeigneter Teilnehmer, Vorbehandlungs-Serumproben und Nachverfolgungsinformationen während des vorgegebenen Studienzeitraums bestimmt. Das diagnostische Modell verwendete drei kontinuierliche Serum-Biomarker als Kandidaten-Prädiktoren. In der Trainingskohorte standen 114 ESCC-Fälle und 80 Nicht-ESCC-Kontrollen zur Verfügung, was 38,0 ESCC-Fällen und 26,7 Nicht-ESCC-Kontrollen pro Kandidaten-Biomarker-Prädiktor entspricht. Dies unterstützte eine logistische Regression mit geringer Dimensionalität und verringerte das Überanpassungsrisiko. Die externe Validierungskohorte umfasste 53 ESCC-Fälle und 30 Nicht-ESCC-Kontrollen und wurde ausschließlich zur Validierung verwendet, ohne dass die Koeffizienten neu geschätzt wurden. Die Begründung der Stichprobengröße orientierte sich an aktuellen Grundsätzen für Vorhersagemodelle, die auf der Häufigkeit des Endpunkts, der Anzahl der Kandidaten-Prädiktorparameter, der erwarteten Modellleistung und einer unabhängigen Validierung basieren, anstatt auf zufälliger Datenaufteilung15.
Rekrutierung der Teilnehmer, Eignungskriterien und Datenerhebung
Die Teilnehmer wurden kontinuierlich während des Studienzeitraums rekrutiert. In Zentrum A wurden ESCC-Patienten aus den Abteilungen Thoraxchirurgie, Gastroenterologie und Onkologie rekrutiert. In Zentrum B wurden ESCC-Patienten über screeningsbasierte Programme in der Gemeinschaft, Nachsorgeakten oder Überweisungswegen identifiziert und nur dann eingeschlossen, wenn eine histopathologische Bestätigung und vorbehandlungsbedingte Serumproben vorlagen. Kontrollpersonen mit benignen Speiseröhrenerkrankungen wurden aus Personen rekrutiert, die eine Gastroskopie wegen Dysphagie, Refluxsymptomen oder Verdacht auf Speiseröhrenveränderungen durchlaufen hatten und nicht-maligne pathologische Befunde aufwiesen. Gesunde Kontrollpersonen wurden aus Routine-Gesundheitsuntersuchungen im gleichen Zeitraum ausgewählt und wiesen keine Vorgeschichte von Malignitäten, aktiven Infektionen, schweren Leber- oder Nierenfunktionsstörungen oder bildgebenden/endoskopischen Hinweisen auf raumfordernde Läsionen der Speiseröhre auf.
ESCC-Patienten waren berechtigt, wenn bei ihnen ein primäres Ösophaguskarzinom durch eine gastroskopische Biopsie oder postoperative Pathologie bestätigt wurde, sie vor der Blutentnahme keine antitumorale Therapie erhalten hatten und vollständige klinische Basisdaten, pathologische Befunde sowie vorbehandlungsbezogene Serumproben vorlagen. Teilnehmer wurden ausgeschlossen, wenn sie ein anderes primäres Malignom aufwiesen, zuvor eine antitumorale Therapie vor der Blutentnahme erhalten hatten, eine schwere akute Infektion, eine aktive Autoimmunerkrankung, nicht kontrollierte Komorbiditäten mit wahrscheinlichem Einfluss auf die Serumtumormarkerwerte, sichtbare Hämolyse, Lipämie, ikterische Interferenz, unzureichendes Serumvolumen, fehlende wichtige klinische Variablen oder keine verfügbaren Nachsagedaten hatten. Die prognostische Analyse beschränkte sich auf ESCC-Patienten mit verfügbarer Tumor-Node-Metastase-(TNM-)Stadieneinteilung und messbaren Überlebens- oder Progressionsendpunkten. Patienten mit einer Nachsorge von weniger als 6 Monaten wurden einbezogen, wenn innerhalb dieses Zeitraums ein Tod oder ein Progress auftrat.
Zu den Basisvariablen gehörten Alter, Geschlecht, Rauchgeschichte, Alkoholkonsumgeschichte, Body-Mass-Index, Studienzentrum, Diagnosegruppe, Läsionstyp sowie die Serumspiegel von CEA, SCC-Ag und CA125. Bei Patienten mit Ösophaguskarzinom (ESCC) wurden zusätzliche Variablen erfasst, darunter Tumorlokalisation, Tumorlänge, histologische Differenzierung, Invasions Tiefe, Lymphknotenstatus, TNM-Stadium, lymphovaskuläre Invasion, initiale Behandlungsmodalität, Rezidiv, Progression, Überlebensstatus und Datum der letzten Nachuntersuchung. Die Rauchgeschichte wurde definiert als gegenwärtiges Rauchen oder früheres regelmäßiges Rauchen von mindestens 100 Zigaretten im Laufe des Lebens. Die Alkoholkonsumgeschichte wurde definiert als Alkoholkonsum von mindestens einmal pro Woche über einen Zeitraum von mehr als sechs Monaten. Die TNM-Klassifikation erfolgte gemäß der achten Auflage des Stadiensystems des American Joint Committee on Cancer/Union for International Cancer Control für Ösophaguskarzinom16. Pathologische Präparate von Fällen mit Ösophaguskarzinom (ESCC) und benignen Läsionen wurden unabhängig voneinander von zwei erfahrenen Pathologen überprüft, wobei bei Uneinigkeiten ein dritter Pathologe zur Konsensfindung hinzugezogen wurde.
Zwei geschulte Forscher extrahierten unabhängig voneinander klinische Informationen aus elektronischen Patientenakten, pathologischen Befunden, Labordaten, Behandlungsunterlagen, Nachsorgeakten und Sterberegistern mithilfe eines vorab festgelegten Datenwörterbuchs. Unstimmigkeiten wurden durch Überprüfung der Originaldokumente geklärt, wobei ein dritter Forscher bei Bedarf die endgültige Entscheidung traf. Die Tumoralokalisation wurde in obere, mittlere oder untere Speiseröhre unterteilt. Die histologische Differenzierung wurde als gut, mäßig oder schlecht differenziert klassifiziert. Die Einheiten für kontinuierliche Variablen wurden vor der Analyse in beiden teilnehmenden Zentren standardisiert.
Blutentnahme, Serumverarbeitung und Biomarkerbestimmung
Alle Blutproben wurden vor der initialen Behandlung, innerhalb von 7 Tagen nach der pathologischen Bestätigung und mindestens 3 Tage vor Therapiebeginn entnommen. Zwischen 7:00 Uhr und 9:00 Uhr morgens wurden periphere venöse Blutproben von 5–8 mL nach einer nächtlichen Fastenphase unter Verwendung von Röhrchen mit Serumtrenngel entnommen. Die Probenentnahme wurde verschoben, wenn der Teilnehmer innerhalb von 72 Stunden vor der Entnahme eine Bluttransfusion, eine Albumininfusion mit großem Volumen oder ein invasives Verfahren erhalten hatte. Die Proben durften 30 Minuten lang bei Raumtemperatur gerinnen und wurden anschließend 10 Minuten lang bei 1.500 × g zentrifugiert. Das Serum wurde abgetrennt und in 3–5 mit Barcode beschriftete Röhrchen portioniert, wobei jede Portion 300–500 µL enthielt. Die Aliquote wurden bis zur Analyse bei −80 °C gelagert, wobei kein Aliquot mehr als einem Einfrier- und Auftauzyklus unterzogen wurde. Proben mit sichtbarer Hämolyse, Lipämie, ikterischer Störung, unzureichendem Volumen oder einer Zeitspanne zwischen Probenentnahme und Einfrieren von mehr als 2 Stunden wurden ausgeschlossen.
Serum-CEA, SCC-Ag und CA125 wurden mithilfe automatisierter chemilumineszenter Immunassay-Systeme gemäß den Herstelleranweisungen bestimmt. Serumproben wurden vor der Analyse bei 4 °C aufgetaut, vorsichtig gemischt und nach Möglichkeit in einem einzigen Durchgang analysiert. Falls eine Batch-Analyse erforderlich war, wurde die Reihenfolge der Proben so organisiert, dass keine systematische Gruppierung nach Erkrankungsstatus, Zentrum, pathologischem Stadium oder Follow-up-Ergebnis vorlag. Proben, deren Werte über dem analytischen Messbereich lagen, wurden gemäß den Reagenzianweisungen verdünnt und erneut getestet. Das Laborpersonal war unkenntlich gegenüber der Zugehörigkeit der Teilnehmergruppe, dem pathologischen Stadium, dem Behandlungsstatus und dem Follow-up-Ergebnis.
Die klinischen Schwellenwerte betrugen 5,0 ng/mL für CEA, 1,5 ng/mL für SCC-Ag und 35 U/mL für CA125, basierend auf lokalen Laborstandards und den Referenzintervallen der Hersteller. Diese Schwellenwerte wurden verwendet, um die Markererhöhung und eine hohe kombinierte Biomarkerlast für assoziationsanalytische Untersuchungen bezüglich klinisch-pathologischer Parameter sowie für die Überlebensstratifizierung zu definieren. Für die diagnostische Klassifizierung wurden datenbasierte Cutoffs ausschließlich im Trainingskollektiv anhand des Youden-Index ausgewählt und anschließend unverändert auf die externe Validierungskohorte angewandt. Vor der Durchführung der Tests erfolgten die Kalibrierung des Instruments und die interne Qualitätskontrolle. Hoch- und niedrigkonzentrierte Kontrollen wurden an jedem Testtag mitgeprüft, und die Tests wurden nur durchgeführt, wenn die Ergebnisse der Qualitätskontrolle innerhalb des akzeptierten Bereichs lagen.
Die Zwischenanalysen- und Zwischenzentren-Reproduzierbarkeit wurde vor der endgültigen Modellanalyse anhand verblindeter Überbrückungs-Serumaliquots und routinemäßiger Qualitätskontrollmaterialien beurteilt. Ein Überbrückungssatz aus 30 verblindeten Serumaliquots wurde in beiden Laborabläufen getestet. Die Zwischenzentren-Variationskoeffizienten betrugen 4,8 % für CEA, 5,6 % für SCC-Ag und 5,1 % für CA125. Die entsprechenden Intraklassen-Korrelationskoeffizienten lagen bei 0,94 für CEA, 0,92 für SCC-Ag und 0,93 für CA125. Es erfolgte keine zenterspezifische Nachkalibrierung der Assays, da keine systematischen Abweichungen den vorab definierten Bereich der laborinternen Qualitätskontrolle überschritten. Die endgültigen Biomarkerwerte wurden vor der Modellerstellung durch Einheitenstandardisierung und statistische Standardisierung basierend auf der Trainingskohorte harmonisiert.
Kombinierte Biomarker-Definition und Modellkonstruktion
Das kombinierte Panel aus Serum-Biomarkern bestand aus CEA, SCC-Ag und CA125. Für die Analyse der klinisch-pathologischen Assoziationen und die Kaplan-Meier-Überlebensstratifizierung wurde jeder Marker entsprechend der vordefinierten klinischen Schwelle als erhöht oder nicht erhöht klassifiziert. Ein hoher kombinierter Biomarker-Belastungswert wurde definiert als Erhöhung von mindestens zwei der drei Marker. Für die diagnostische und prognostische Modellierung wurden die drei Biomarker als kontinuierliche Variablen beibehalten, mittels des natürlichen Logarithmus transformiert und mit Hilfe des Mittelwerts und der Standardabweichung aus der Trainingskohorte standardisiert. Die gleichen Transformationsparameter wurden auf die externe Validierungskohorte angewandt.
Die standardisierten log-transformierten Biomarker wurden anhand der folgenden Parameter der Trainingskohorte berechnet:
zlnCEA = (lnCEA - 1,42)/0,71
zlnSCC - Ag = (lnSCC - Ag - 0,28)/0,57
zlnCA125 = (lnCA125 - 3,05)/0,62
Die kombinierte diagnostische Panel wurde mittels multivariater logistischer Regression in der Trainingskohorte erstellt. Der diagnostische lineare Prädiktor wurde wie folgt berechnet:
LPdiagnostic = -0,46 + 0,59 × zlnCEA + 1,05 × zlnSCC - Ag+0,42 × zlnCA125
Die individuelle vorhergesagte Wahrscheinlichkeit von ESCC wurde wie folgt berechnet:
PESCC = 1/1+exp(-LPdiagnostic)
Die optimale diagnostische Grenze für das kombinierte Panel betrug 0,57, ausgewählt in der Trainingskohorte anhand des Youden-Index und unverändert in der externen Validierung angewandt.
Die prognostische Biomarker-Score wurde mithilfe einer Cox-Proportionalitäts-Modellierung konstruiert. Als Vergleichsmodell für die anatomische Stadieneinteilung diente das alleinige TNM-Stadium-Modell. Das Biomarker-Modell verwendete die kombinierte Biomarker-Belastung. Das klinisch-pathologische Modell umfasste Alter, Tumorlänge, histologische Differenzierung, pT-Stadium, Lymphknotenstatus und initiale Behandlungsmodalität. Dem klinisch-pathologischen Modell wurde im integrierten Modell eine hohe kombinierte Biomarker-Belastung hinzugefügt. Für das OS wurde der integrierte prognostische lineare Prädiktor wie folgt berechnet:
LPOS = 0.33 × Alter ≥65 Jahre + 0,36 × Tumorlänge ≥5 cm + 0,25 × schlechte Differenzierung + 0,31 × pT3
- 4-Stadium + 0,66 × Lymphknotenmetastase
-0,58 × Radikaloperation - basierende Behandlung
+0,81 × highcombinedbiomarkerburden
Für PFS wurde der integrierte prognostische lineare Prädiktor wie folgt berechnet:
LPPFS = 0,27 × Alter ≥ 65 Jahre + 0,31 × Tumorlänge ≥ 5 cm + 0,24 × schlechte Differenzierung + 0,25 × pT3
-4. Stadium + 0,55 × Lymphknotenmetastase
– 0,46 × Radikaloperation – basierende Behandlung
+0,71 × hohe kombinierte Biomarkerlast
Alle binären Variablen in den prognostischen Gleichungen wurden mit 1 kodiert, wenn vorhanden, und mit 0, wenn nicht vorhanden. Alle Koeffizienten, Standardisierungsparameter und Schwellenwerte wurden anhand der Trainingskohorte geschätzt und ohne erneute Optimierung direkt auf die externe Validierungskohorte angewandt.
Diagnostische und prognostische Endpunkte
Die pathologische Diagnose wurde als diagnostischer Referenzstandard verwendet. Der primäre diagnostische Endpunkt war die Unterscheidung zwischen ESCC und Nicht-ESCC-Kontrollen, einschließlich benigner Speiseröhrenerkrankungen und gesunder Kontrollen. Zu den sekundären diagnostischen Endpunkten gehörten die Unterscheidung zwischen ESCC und benigner Speiseröhrenerkrankung, ESCC und gesunden Kontrollen sowie zwischen frühen ESCC-Stadien und nicht malignen Kontrollen. Frühe ESCC-Stadien wurden als TNM-Stadium I–II definiert. Die diagnostische Leistung wurde anhand von Receiver-Operating-Characteristic-Kurven, Fläche unter der Kurve, Sensitivität, Spezifität, positivem Vorhersagewert, negativem Vorhersagewert, Genauigkeit, F1-Score, Precision-Recall-Kurven, Kalibrationskurven, Brier-Scores, Kalibrationsintercepts, Kalibrierungssteigungen und der Entscheidungskurvenanalyse bewertet.
Die prognostische Analyse wurde ausschließlich bei Patienten mit Ösophagus-Plattenepithelkarzinom (ESCC) durchgeführt. Der primäre prognostische Endpunkt war das Gesamtüberleben (OS), definiert als der Zeitraum vom Datum der ersten definitiven antitumoralen Behandlung bis zum Tod aus beliebiger Ursache. Der sekundäre prognostische Endpunkt war das progressionsfreie Überleben (PFS), definiert als der Zeitraum von der Erstbehandlung bis zur dokumentierten Krankheitsprogression, Rezidiv oder dem Tod. Bei Patienten, die einer kurativen Operation unterzogen wurden, galt ein lokales Rezidiv, eine Fernmetastasierung oder ein tumorbedingter Tod als PFS-Ereignis, sofern nach der Operation kein Residualläsion vorhanden war. Bei Patienten, die eine definitive Radiochemotherapie oder eine palliative Behandlung erhielten, wurde die Progression anhand bildgebender Verfahren, endoskopischer Beurteilung und klinischer Evaluation festgestellt.
Nachsagedaten wurden durch ambulante Besuche, Krankenhausaufenthaltsakten, Gemeinschaftsnachsorgeakten, Überweisungsakten und telefonischen Kontakt erhoben. Die Patienten wurden in den ersten 2 Jahren nach der Behandlung alle 3 Monate, von Jahr 3 bis 5 alle 6 Monate und danach jährlich nachbetreut. Patienten ohne Tod, Rezidiv oder Progress zum Zeitpunkt des letzten Nachsagekontakts wurden zum Datum ihres letzten bestätigten Kontakts zensiert. Die Ergebnisereignisse wurden unabhängig voneinander von zwei Untersuchern überprüft. Wenn sich die Ereignisdaten zwischen den Quellen unterschieden, wurden Bildgebungsberichte, Entlassungsakten, endoskopische Berichte, Überweisungsdokumente und Sterberegisterakten geprüft, um das endgültige Ereignisdatum festzulegen. Das Nachsage-Stichtag war für beide teilnehmenden Zentren der 31. Dezember 2025.
Statistische Analyse
Die statistische Analyse wurde mit R 4.3.2 und SPSS 27.0 durchgeführt. Stetige Variablen wurden je nach Verteilung als Mittelwert ± Standardabweichung oder Median (Interquartilbereich) zusammengefasst. Kategoriale Variablen wurden als Häufigkeiten und Prozentanteile dargestellt. Gruppenvergleiche wurden je nach Eignung mittels des t-Tests für unabhängige Stichproben, des Mann-Whitney-U-Tests, der einfaktoriellen Varianzanalyse, des Kruskal-Wallis-Tests, des Chi-Quadrat-Tests oder des exakten Fisher-Tests durchgeführt. Fehlende Daten wurden vor der Modellerstellung überprüft. Teilnehmer mit fehlendem wichtigem Diagnosestatus, fehlenden prätherapeutischen Biomarkerwerten, fehlendem TNM-Stadium oder fehlenden Daten zu Überleben/Progression wurden gemäß den vorab festgelegten Einschlusskriterien ausgeschlossen. Nicht wesentliche Kovariaten mit einem Fehlendanteil unter 5,0 % wurden mittels multipler Imputation durch verkettete Gleichungen mit 20 imputierten Datensätzen und 20 Iterationen behandelt17. Die Imputationsmodelle umfassten Diagnosegruppe, Zentrum, Alter, Geschlecht, Rauchgeschichte, Alkoholkonsum, Biomarkerwerte, TNM-Stadium, Behandlungsmodalität, OS-Status, PFS-Status und Nachbeobachtungszeit. Analysen an vollständigen Fällen wurden als Sensitivitätsanalysen durchgeführt.
Das diagnostische Modell wurde an der Trainingskohorte entwickelt. Zur Berechnung der Fläche unter der Kurve, Sensitivität, Spezifität, positiven Vorhersagewerte, negativen Vorhersagewerte, Genauigkeit und des F1-Scores wurden Receiver-Operating-Characteristic-Kurven verwendet. Die Flächen unter den Kurven wurden mittels des DeLong-Tests18 verglichen. Als ergänzende diagnostische Bewertung kamen Precision-Recall-Kurven zum Einsatz. Die Kalibrierung wurde anhand von Kalibrierungskurven, Brier-Scores, Kalibrierungsinterzepten und Kalibrierungssteigungen beurteilt. Die klinische Nützlichkeit wurde mittels Decision-Curve-Analyse19 evaluiert. Die interne Validierung erfolgte durch 1.000 Bootstrap-Stichproben, um die diskriminative Leistung und Kalibrierung unter Korrektur für Optimismus zu schätzen. Die externe Validierung wurde durchgeführt, indem die aus der Trainingskohorte stammende Gleichung, Standardisierungsparameter und der Schwellenwert direkt auf Zentrum B angewandt wurden, ohne erneutes Anpassen.
Für die prognostische Analyse wurden Kaplan-Meier-Kurven erstellt, um das Gesamtüberleben (OS) und das progressionsfreie Überleben (PFS) zwischen den Gruppen mit unterschiedlicher Biomarker-Belastung zu vergleichen, und der Log-Rang-Test wurde zur Beurteilung von Überlebensunterschieden verwendet. Die Cox-Regression mit proportionalen Hazards diente der Schätzung von Hazard-Ratios und 95%-Konfidenzintervallen. Variablen mit klinischer Relevanz sowie Variablen, die in der univariablen Analyse mit den Ergebnissen assoziiert waren, wurden für multivariable Modelle berücksichtigt. Die Modellleistung wurde anhand des Harrellschen C-Index, zeitabhängiger ROC-Kurven (receiver operating characteristic), Kalibrierungskurven, Brier-Scores und der Entscheidungskurvenanalyse bewertet20. Die prognostische Kalibrierung wurde nach 1, 3 und 5 Jahren überprüft. Die Annahme proportionaler Hazards wurde mithilfe von Schoenfeld-Residuen und Log-minus-Log-Überlebensplots geprüft21. Die globalen Schoenfeld-Test-P-Werte betrugen 0,42 für das OS-Modell und 0,37 für das PFS-Modell. Multikollinearität wurde vor der endgültigen Interpretation des Modells anhand von Varianzinflationsfaktoren überprüft, wobei alle Kandidatenvariablen in den endgültigen diagnostischen und prognostischen Modellen Varianzinflationsfaktoren unter 2,1 aufwiesen. Der zusätzliche prognostische Wert jenseits des TNM-Stadiums und konventioneller klinisch-pathologischer Modelle wurde mittels Net Reclassification Improvement und Integrated Discrimination Improvement bewertet. Die im Trainingskollektiv geschätzten Modellparameter wurden ohne Neuanpassung auf die externe Validierungskohorte angewandt22. Alle Tests waren zweiseitig, und ein P < 0,05 galt als statistisch signifikant.