Deze studie met menselijke deelnemers, klinische dossiers en serummonsters werd uitgevoerd in overeenstemming met de Verklaring van Helsinki en de institutionele vereisten voor onderzoek met menselijke proefpersonen13. Het studieprotocol is beoordeeld en goedgekeurd door het Jiangsu Province Hospital of Chinese Medicine onder goedkeuringsnummer CJ28J12-45. Voor retrospectief verzamelde geanonimiseerde klinische gegevens en residuele serummonsters is de vereiste voor aanvullende schriftelijke geïnformeerde toestemming kwijtgescholden in overeenstemming met de institutionele ethische regelgeving. Voor prospectief verzamelde monsters is voorafgaand aan de inclusie schriftelijke geïnformeerde toestemming verkregen. Alle deelnemersinformatie en serummonsters werden geanonimiseerd vóór de analyse. Studie-identificatoren werden gebruikt voor datakoppeling, en er is geen persoonlijk identificeerbare informatie opgenomen in de analytische dataset.
Studieopzet en cohorttoewijzing
Deze ambispectieve observationele studie in twee centra evalueerde de diagnostische en prognostische waarde van een gecombineerd serumbiomarkerpaneel bestaande uit carcinoembryonaal antigeen (CEA), plaveiselcelcarcinoom-antigeen (SCC-Ag) en koolhydraatantigeen 125 (CA125) bij plaveiselcelcarcinoom van de slokdarm (ESCC). Retrospectieve gegevens werden verkregen uit geschikte klinische dossiers en resterende serummonsters van vóór de behandeling, en prospectieve werving werd gebruikt om aanvullende serummonsters van vóór de behandeling te verzamelen en de follow-up binnen hetzelfde studiekader te voltooien. Deelnemers werden geworven uit één tertiair ziekenhuis en één gemeentelijk gezondheidscentrum tussen januari 2021 en december 2024, waarbij de follow-up voortduurde tot 31 december 2025. Centrum A was het tertiaire ziekenhuis en diende als trainingscohort voor modelontwikkeling, coëfficiëntschatting, selectie van afkapwaarden en interne validatie. Centrum B was het gemeentelijk gezondheidscentrum en diende als extern validatiecohort. De studie is gerapporteerd in overeenstemming met de Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD)-richtlijnen14.
In totaal werden 277 deelnemers geïncludeerd, waarvan 194 afkomstig waren van Centrum A en 83 van Centrum B. De trainingscohort bestond uit 114 ESCC-patiënten, 40 controles met goedaardige slokdarmziekten en 40 gezonde controles. De externe validatiecohort bestond uit 53 ESCC-patiënten, 15 controles met goedaardige slokdarmziekten en 15 gezonde controles. Om informatielekkage te verminderen en de stabiliteit van het model in een onafhankelijke klinische setting te testen, werd toewijzing op basis van het centrum gebruikt in plaats van een willekeurige verdeling. De diagnostische analyse vergeleek de serumspiegels van CEA, SCC-Ag en CA125 tussen ESCC-patiënten, patiënten met goedaardige slokdarmziekten en gezonde controles. De prognostische analyse was beperkt tot ESCC-patiënten en beoordeelde of de gecombineerde biomarkerlast de risicostratificatie voor algehele overleving (OS) en progressievrije overleving (PFS) verbeterde. Het studieontwerp, de screening van deelnemers, de cohorttoewijzing, de workflow voor serumbiomarkers, de ontwikkeling van het diagnostische model, de externe validatie en de vergelijking van het prognostische model zijn samengevat in Figuur 1.
De stroom van deelnemers omvatte screening, exclusies, diagnostische groepering en toewijzing per centrum. In totaal werden 326 personen gescreend. Negenenveertig personen werden uitgesloten, waaronder 12 met een antitumorbehandeling voorafgaand aan de bloedafname, 8 met een andere primaire maligniteit, 9 met een ernstige infectie, auto-immuunziekte of ongecontroleerde comorbiditeiten die waarschijnlijk de serumtumormarkers beïnvloeden, 11 met hemolyse, lipemie, geelzuchtinterferentie, onvoldoende serumvolume of een verlengd verwerkingsinterval, en 9 met ontbrekende cruciale klinische, pathologische of follow-upgegevens. De uiteindelijke analytische populatie bestond uit 277 deelnemers: 167 ESCC-patiënten, 55 controles met goedaardige slFockziekten en 55 gezonde controles.
Rationale voor steekproefomvang
De steekproefomvang werd bepaald door de opeenvolgende beschikbaarheid van geschikte deelnemers, serummonsters vóór de behandeling en follow-upinformatie gedurende de vooraf gedefinieerde studieperiode. Het diagnostische model maakte gebruik van drie continue serumbiomarkers als potentiële predictoren. In de trainingscohort waren 114 ESCC-gevallen en 80 non-ESCC-controles beschikbaar, wat overeenkomt met 38,0 ESCC-gevallen en 26,7 non-ESCC-controles per potentiële biomarker-predictor. Dit ondersteunde laag-dimensionale logistieke regressie en verminderde het risico op overfitting. De externe validatiecohort bestond uit 53 ESCC-gevallen en 30 non-ESCC-controles en werd uitsluitend gebruikt voor validatie zonder herberekening van de coëfficiënten. De onderbouwing van de steekproefomvang was in lijn met huidige principes voor predictiemodellen, waarbij de nadruk ligt op de frequentie van de uitkomst, het aantal parameters van de potentiële predictoren, de verwachte modelprestaties en onafhankelijke validatie in plaats van willekeurige splitsing van de gegevens15.
Werving van deelnemers, geschiktheid en gegevensverzameling
Deelnemers werden continu gedurende de studieperiode geïncludeerd. In Centrum A werden ESCC-patiënten geworven vanuit de afdelingen thoraxchirurgie, gastroenterologie en oncologie. In Centrum B werden ESCC-patiënten geïdentificeerd via gemeenschapsgebaseerde screening, follow-upgegevens of verwijzingspaden, en zij werden alleen geïncludeerd wanneer pathologische bevestiging en serummonsters van vóór de behandeling beschikbaar waren. Controles met benigne slokdarmaandoeningen werden geworven uit personen die een gastroscopie ondergingen vanwege dysfagie, refluxsymptomen of vermoedelijke afwijkingen aan de slokdarm en bij wie non-maligne pathologische bevindingen waren vastgesteld. Gezonde controles werden geselecteerd uit routineonderzoeken in dezelfde periode en hadden geen voorgeschiedenis van maligniteit, actieve infectie, ernstige lever- of nierdysfunctie, of beeldvormend/endoscopisch bewijs van ruimte-innemende laesies in de slokdarm.
ESCC-patiënten kwamen in aanmerking als er sprake was van primair ESCC bevestigd door gastroscopisch biopt of postoperatieve pathologie, als er geen antitumorbehandeling was ondergaan vóór de bloedafname, en als er volledige baseline klinische gegevens, pathologische verslagen en serummonsters van vóór de behandeling beschikbaar waren. Deelnemers werden uitgesloten bij de aanwezigheid van een andere primaire maligniteit, eerdere antitumorbehandeling vóór de bloedafname, een ernstige acute infectie, een actieve auto-immuunziekte, ongecontroleerde comorbiditeiten die waarschijnlijk de serumwaarden van tumormarkers beïnvloeden, zichtbare hemolyse, lipemie, geelzuchtinterferentie, onvoldoende serumvolume, ontbrekende belangrijke klinische variabelen of niet beschikbare follow-upgegevens. De prognostische analyse was beperkt tot ESCC-patiënten met een bekend tumor-node-metastasis (TNM) stadium en meetbare overlevings- of progressie-uitkomsten. Patiënten met een follow-up van minder dan 6 maanden werden behouden als er binnen die periode overlijden of progressie optrad.
De basisvariabelen omvatten leeftijd, geslacht, rookgeschiedenis, geschiedenis van alcoholconsumptie, bodymassindex, deelnemend centrum, diagnostische groep, type laesie en serumspiegels van CEA, SCC-Ag en CA125. Voor ESCC-patiënten omvatten aanvullende variabelen de tumorlocatie, tumorlengte, histologische differentiatie, invasiediepte, lymfekliertstatus, TNM-stadium, lymfovasculaire invasie, initiële behandelmethode, recidief, progressie, overlevingsstatus en de datum van de laatste follow-up. Rookgeschiedenis werd gedefinieerd als huidige rookstatus of eerdere regelmatige consumptie van ten minste 100 sigaretten gedurende het leven. Geschiedenis van alcoholconsumptie werd gedefinieerd als alcoholinname ten minste één keer per week gedurende meer dan 6 maanden. De TNM-stadiëring werd bepaald volgens de achtste editie van het stadiëringssysteem voor slokdarmkanker van het American Joint Committee on Cancer/Union for International Cancer Control16. Pathologische coupes van ESCC-gevallen en gevallen met goedaardige laesies werden onafhankelijk beoordeeld door twee senior pathologen, waarbij eventuele meningsverschillen werden opgelost door een derde patholoog.
Twee getrainde onderzoekers extraheerden onafhankelijk klinische informatie uit elektronische medische dossiers, pathologierapporten, laboratoriuminformatiesystemen, behandelingsverslagen, follow-upverslagen en overlijdensregistraties met behulp van een vooraf gedefinieerd datadictionary. Verschillen werden opgelost door de originele dossiers te controleren, waarbij een derde onderzoeker indien nodig de definitieve beslissing nam. De tumorlocatie werd gecategoriseerd als de bovenste, middelste of onderste thoracale slokdarm. Histologische differentiatie werd gecategoriseerd als goede, matige of slechte differentiatie. Eenheden van continue variabelen werden voor de analyse gestandaardiseerd over beide deelnemende centra.
Bloedafname, serumverwerking en biomarkermeting
Alle bloedmonsters werden afgenomen vóór de initiële behandeling, binnen 7 dagen na pathologische bevestiging en ten minste 3 dagen vóór het begin van de behandeling. Perifere veneuze bloedmonsters van 5–8 mL werden verzameld tussen 7:00 AM en 9:00 AM na een nacht vasten met behulp van serumscheidingsbuizen. De afname werd uitgesteld indien de deelnemer binnen 72 h vóór de afname een bloedtransfusie, een albumininfusie van groot volume of een invasieve procedure had ondergaan. De monsters kregen 30 min de tijd om te stollen bij kamertemperatuur en werden vervolgens 10 min gecentrifugeerd bij 1.500 × g . Het serum werd gescheiden en verdeeld over 3–5 met barcodes gelabelde buisjes, met 300–500 µL per aliquot. De aliquoten werden bewaard bij −80°C tot aan de analyse, waarbij geen enkele aliquot aan meer dan één vries-dooi-cyclus werd blootgesteld. Monsters met zichtbare hemolyse, lipemie, geelzuchtinterferentie, onvoldoegend volume of een interval tussen afname en invriezen van langer dan 2 h werden uitgesloten.
Serum CEA, SCC-Ag en CA125 werden gemeten met geautomatiseerde chemiluminescente immunoassay-systemen volgens de instructies van de fabrikant. Serumaliquots werden voor het testen ontdooid bij 4°C, voorzichtig gemengd en indien mogelijk in één run geanalyseerd. Indien batch-testen noodzakelijk was, werd de volgorde van de monsters zo bepaald dat systematische groepering op basis van ziektestatus, centrum, pathologisch stadium of follow-upresultaat werd vermeden. Monsters die het analytische meetbereik overschreden, werden verdund en opnieuw getest volgens de instructies van het reagens. Het laboratoriumpersoneel was geblinderd voor de deelnemersgroep, het pathologische stadium, de behandelstatus en het follow-upresultaat.
De klinische drempelwaarden waren 5,0 ng/mL voor CEA, 1,5 ng/mL voor SCC-Ag en 35 U/mL voor CA125, gebaseerd op lokale laboratoriumnormen en de referentie-intervallen van de fabrikanten. Deze drempelwaarden werden gebruikt om marker-elevatie en een hoge gecombineerde biomarkerlast te definiëren voor clinicopathologische associatie- en overlevingsstratificatie-analyses. Voor de diagnostische classificatie werden datagestuurde afkapwaarden uitsluitend in de trainingscohort geselecteerd met behulp van de Youden-index en vervolgens ongewijzigd toegepast op de externe validatiecohort. Instrumentkalibratie en interne kwaliteitscontrole werden uitgevoerd vóór de tests. Controles op hoog en laag niveau werden tijdens elke testdag meegenomen, en het testen ging alleen door wanneer de kwaliteitscontroleresultaten binnen het geaccepteerde bereik vielen.
De reproduceerbaarheid tussen assays en tussen centra werd vóór de uiteindelijke modelanalyse beoordeeld met behulp van geblindeerde bridging-serumaliquots en routinematige kwaliteitscontrolematerialen. Een bridging-set van 30 geblindeerde serumaliquots werd getest over de twee laboratoriumworkflows. De variatiecoëfficiënten tussen de centra waren 4,8% voor CEA, 5,6% voor SCC-Ag en 5,1% voor CA125. De overeenkomstige intraclass-correlatiecoëfficiënten waren 0,94 voor CEA, 0,92 voor SCC-Ag en 0,93 voor CA125. Er werd geen centrumspecifieke herkalibratie van de assay uitgevoerd, omdat geen enkele systematische afwijking het vooraf gedefinieerde laboratorium-kwaliteitscontrolegebied overschreed. De uiteindelijke biomarkerwaarden werden geharmoniseerd door middel van unit-standaardisatie en statistische standaardisatie op basis van de trainingscohort vóór de constructie van het model.
Gecombineerde biomarkerdefinitie en modelconstructie
Het gecombineerde serum-biomarkerpanel bestond uit CEA, SCC-Ag en CA125. Voor de clinicopathologische associatieanalyse en Kaplan-Meier overlevingsstratificatie werd elke marker geclassificeerd als verhoogd of niet-verhoogd op basis van de vooraf gedefinieerde klinische drempelwaarde. Een hoge gecombineerde biomarkerlast werd gedefinieerd als de verhoging van ten minste twee van de drie markers. Voor diagnostische en prognostische modellering werden de drie biomarkers behouden als continue variabelen, getransformeerd met behulp van de natuurlijke logaritme en gestandaardiseerd met het gemiddelde en de standaarddeviatie van de trainingscohort. Dezelfde transformatieparameters werden toegepast op de externe validatiecohort.
De gestandaardiseerde log-getransformeerde biomarkers werden berekend met behulp van de volgende parameters van de trainingscohort:
zlnCEA = (lnCEA - 1.42)/0.71
zlnSCC - Ag = (lnSCC - Ag - 0.28)/0.57
zlnCA125 = (lnCA125 - 3.05)/0.62
Het gecombineerde diagnostische panel werd opgesteld met behulp van multivariabele logistische regressie in de trainingscohort. De diagnostische lineaire voorspeller werd als volgt berekend:
LPdiagnostic = -0,46 + 0,59 × zlnCEA + 1,05 × zlnSCC - Ag+0,42 × zlnCA125
De individuele voorspelde waarschijnlijkheid van ESCC werd als volgt berekend:
PESCC = 1/1+exp(-LPdiagnostic)
De optimale diagnostische afkapwaarde voor het gecombineerde panel was 0.57, geselecteerd in de trainingscohort met behulp van de Youden-index en ongewijzigd toegepast bij de externe validatie.
De prognostische biomarkerscore werd opgesteld met behulp van Cox proportional hazards-modellering. Het model op basis van enkel het TNM-stadium werd gebruikt als vergelijking voor de anatomische stadiëring. Het biomarkermodel maakte gebruik van de gecombineerde biomarkerbelasting. Het clinicopathologische model maakte gebruik van leeftijd, tumorlengte, histologische differentiatie, pT-stadium, lymfekliersituatie en de initiële behandelmethode. Het geïntegreerde model voegde een hoge gecombineerde biomarkerbelasting toe aan het clinicopathologische model. Voor de OS werd de geïntegreerde prognostische lineaire predictor als volgt berekend:
LPOS = 0,33 × leeftijd ≥65 jaar + 0,36 × tumorlengte ≥5cm + 0,25 × slechte differentiatie + 0,31 × pT3
- stadium 4 + 0,66 × lymfekliermetastase
-0.58 × behandeling op basis van radicale chirurgie
+0,81 × hoge gecombineerde biomarkerlast
Voor de PFS werd de geïntegreerde prognostische lineaire predictor als volgt berekend:
LPPFS = 0,27 × leeftijd ≥ 65 jaar + 0,31 × tumorgrootte ≥ 5 cm + 0,24 × slechte differentiatie + 0,25 × pT3
-4stadium + 0,55 × lymfekliermetastase
- 0,46 × behandeling op basis van radicale chirurgie
+0,71 × hoge gecombineerde biomarkerlast
Alle binaire variabelen in de prognostische vergelijkingen werden gecodeerd als 1 bij aanwezigheid en 0 bij afwezigheid. Alle coëfficiënten, standaardisatieparameters en afkapwaarden werden geschat in de trainingscohort en direct toegepast op de externe validatiecohort zonder heroptimalisatie.
Diagnostische en prognostische eindpunten
Pathologische diagnose werd gebruikt als de diagnostische referentiestandaard. Het primaire diagnostische eindpunt was het onderscheid tussen ESCC en non-ESCC controles, inclusief benigne slziekte en gezonde controles. Secundaire diagnostische eindpunten omvatten het onderscheid tussen ESCC en benigne slziekte, ESCC en gezonde controles, en ESCC in een vroeg stadium en non-maligne controles. ESCC in een vroeg stadium werd gedefinieerd als ziekte in TNM-stadium I–II. De diagnostische prestaties werden geëvalueerd met behulp van receiver operating characteristic-curven, area under the curve, sensitiviteit, specificiteit, positief voorspellende waarde, negatief voorspellende waarde, nauwkeurigheid, F1-score, precision-recall-curven, kalibratiecurven, Brier-scores, kalibratie-intercepten, kalibratieshellingen en decision curve-analyse.
De prognostische analyse werd uitsluitend uitgevoerd bij patiënten met ESCC. Het primaire prognostische eindpunt was OS, gedefinieerd als het interval vanaf de datum van de eerste definitieve antitumorbehandeling tot het overlijden door een willekeurige oorzaak. Het secundaire prognostische eindpunt was PFS, gedefinieerd als het interval vanaf de eerste behandeling tot gedocumenteerde ziekteprogressie, recidief of overlijden. Voor patiënten die een radicale chirurgische ingreep ondergingen, werden lokaal recidief, metastase op afstand of tumorgerelateerd overlijden geteld als PFS-event wanneer er na de operatie geen residuele laesie aanwezig was. Voor patiënten die definitieve chemoradiotherapie of palliatieve behandeling ontvingen, werd progressie bepaald op basis van beeldvormingsresultaten, endoscopische beoordeling en klinische evaluatie.
Follow-upgegevens werden verzameld via poliklinische bezoeken, ziekenhuisgegevens, follow-upgegevens uit de gemeenschap, verwijzingsgegevens en telefonisch contact. Patiënten werden elke 3 maanden gevolgd gedurende de eerste 2 jaar na de behandeling, elke 6 maanden van jaar 3 tot 5, en daarna jaarlijks. Patiënten bij wie geen overlijden, recidief of progressie was vastgesteld op de datum van de laatste follow-up, werden gecensureerd bij hun laatste bevestigde contact. Uitkomstgebeurtenissen werden onafhankelijk beoordeeld door twee onderzoekers. Wanneer gebeurtenisdata tussen bronnen verschilden, werden beeldvormingsverslagen, ontslagdocumenten, endoscopieverslagen, verwijzingsdocumenten en overlijdensaktes gecontroleerd om de definitieve datum van de gebeurtenis te bepalen. De uiterste datum voor de follow-up was 31 december 2025 voor beide deelnemende centra.
Statistische analyse
Statistische analyses werden uitgevoerd met R 4.3.2 en SPSS 27.0. Continue variabelen werden samengevat als gemiddelde ± standaarddeviatie of mediaan (interkwartielafstand), afhankelijk van de distributie. Categorische variabelen werden samengevat als aantallen en percentages. Vergelijkingen tussen groepen werden uitgevoerd met de t-toets voor onafhankelijke steekproeven, de Mann-Whitney U-toets, eenweg-variantieanalyse, de Kruskal-Wallis-toets, de chi-kwadraattoets of de exacte toets van Fisher, naar gelang van de situatie. Ontbrekende gegevens werden beoordeeld vóór de constructie van het model. Deelnemers bij wie essentiële gegevens over de diagnostische status, biomarkerwaarden vóór de behandeling, het TNM-stadium of de overleving/progressie ontbraken, werden uitgesloten op basis van de vooraf gedefinieerde geschiktheidscriteria. Niet-essentiële covariabelen met een percentage ontbrekende gegevens lager dan 5,0% werden behandeld middels multiple imputatie via geketende vergelijkingen met 20 geïmputeerde datasets en 20 iteraties17. De imputatiemodellen omvatten de diagnostische groep, het centrum, leeftijd, geslacht, rookgeschiedenis, alcoholgebruik, biomarkerwaarden, TNM-stadium, behandelmethode, OS-status, PFS-status en follow-up tijd. Complete-case analyses werden uitgevoerd als sensitiviteitsanalyses.
Het diagnostische model werd ontwikkeld in de trainingscohort. Receiver operating characteristic-curves werden gebruikt om de oppervlakte onder de curve (AUC), sensitiviteit, specificiteit, positief voorspellende waarde, negatief voorspellende waarde, accuratesse en de F1-score te berekenen. De oppervlaktes onder de curve werden vergeleken met behulp van de DeLong-test18. Precision-recall-curves werden gebruikt als aanvullende diagnostische evaluatie. De kalibratie werd beoordeeld met behulp van kalibratiecurves, Brier-scores, kalibratie-intercepten en kalibratieshellingen. Het klinisch nut werd geëvalueerd met behulp van decision curve-analyse19. Interne validatie werd uitgevoerd met 1.000 bootstrap-resamples om de optimisme-gecorrigeerde discriminatie en kalibratie te schatten. Externe validatie werd uitgevoerd door de vergelijking van de trainingscohort, de standaardisatieparameters en de afkapwaarde direct toe te passen op Centrum B zonder herfitting.
Voor de prognostische analyse werden Kaplan-Meier-curves gegenereerd om de OS en PFS tussen groepen met verschillende biomarker-belasting te vergelijken, en de log-ranktest werd gebruikt om verschillen in overleving te beoordelen. Cox proportional hazards-regressie werd gebruikt om hazard ratio's en 95% betrouwbaarheidsintervallen te schatten. Variabelen met klinische relevantie en variabelen die in de univariabele analyse geassocieerd waren met uitkomsten, werden meegenomen in multivariabele modellen. De prestaties van het model werden geëvalueerd met behulp van de C-index van Harrell, tijdsafhankelijke receiver operating characteristic-curves, kalibratiecurves, Brier-scores en decision curve-analyse20. Prognostische kalibratie werd beoordeeld na 1, 3 en 5 jaar. De aanname van proportionele hazards werd getoetst met behulp van Schoenfeld-residuen en log-minus-log survival-plots21. De P-waarden van de globale Schoenfeld-test waren 0,42 voor het OS-model en 0,37 voor het PFS-model. Multicollineariteit werd beoordeeld met behulp van variance inflation factors vóór de uiteindelijke modelinterpretatie, en alle kandidaatvariabelen in de definitieve diagnostische en prognostische modellen hadden variance inflation factors onder de 2,1. De incrementele prognostische waarde bovenop het TNM-stadium en conventionele clinicopathologische modellen werd beoordeeld met behulp van net reclassification improvement en integrated discrimination improvement. Modelparameters die in de trainingscohort waren geschat, werden toegepast op de externe validatiecohort zonder herfitting22. Alle tests waren tweezijdig, en P < 0,05 werd beschouwd als statistisch significant.