Cette étude impliquant des participants humains, des dossiers cliniques et des échantillons de sérum a été menée conformément à la Déclaration d'Helsinki et aux exigences institutionnelles relatives à la recherche sur des sujets humains13. Le protocole de l'étude a été examiné et approuvé par l'Hôpital provincial de médecine chinoise du Jiangsu sous le numéro d'approbation CJ28J12-45. Pour les données cliniques anonymisées et les échantillons résiduels de sérum collectés de façon rétrospective, l'exigence de consentement éclairé écrit supplémentaire a été levée conformément aux règlements éthiques institutionnels. Pour les échantillons collectés de façon prospective, un consentement éclairé écrit a été obtenu avant l'inclusion. Toutes les informations concernant les participants et les échantillons de sérum ont été anonymisées avant l'analyse. Des identifiants d'étude ont été utilisés pour la liaison des données, et aucune information permettant d'identifier une personne n'a été incluse dans le jeu de données analytique.
Conception de l'étude et répartition de la cohorte
Cette étude observationnelle rétrospective à deux centres a évalué la valeur diagnostique et pronostique d'un panel combiné de biomarqueurs sériques incluant l'antigène carcino-embryonnaire (CEA), l'antigène du carcinome épidermoïde (SCC-Ag) et l'antigène carbonhydraté 125 (CA125) dans le carcinome épidermoïde de l'œsophage (ESCC). Des données rétrospectives ont été extraites des dossiers cliniques éligibles et des échantillons sériques résiduels pré-thérapeutiques, tandis qu'un recrutement prospectif a permis de collecter des échantillons sériques supplémentaires pré-thérapeutiques et d'achever le suivi dans le cadre du même protocole d'étude. Les participants ont été recrutés dans un hôpital tertiaire et un centre de santé communautaire entre janvier 2021 et décembre 2024, avec un suivi poursuivi jusqu'au 31 décembre 2025. Le centre A, un hôpital tertiaire, a constitué la cohorte d'apprentissage destinée au développement du modèle, à l'estimation des coefficients, à la sélection des seuils de coupure et à la validation interne. Le centre B, un centre de santé communautaire, a servi de cohorte de validation externe. L'étude a été rapportée conformément aux recommandations du guide Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD)14.
Un total de 277 participants a été inclus, dont 194 provenant du centre A et 83 du centre B. La cohorte d'entraînement comprenait 114 patients atteints de CÉO, 40 témoins présentant une maladie bénigne de l'œsophage et 40 témoins sains. La cohorte de validation externe comprenait 53 patients atteints de CÉO, 15 témoins présentant une maladie bénigne de l'œsophage et 15 témoins sains. Une répartition par centre a été utilisée à la place d'une division aléatoire afin de réduire les fuites d'information et tester la stabilité du modèle dans un cadre clinique indépendant. L'analyse diagnostique a comparé les taux sériques de CEA, de SCC-Ag et de CA125 chez les patients atteints de CÉO, les patients souffrant d'une maladie bénigne de l'œsophage et les témoins sains. L'analyse pronostique a été limitée aux patients atteints de CÉO et a évalué si la charge combinée de biomarqueurs améliorait la stratification du risque de survie globale (OS) et de survie sans progression (PFS). Le schéma de l'étude, le dépistage des participants, l'attribution des cohortes, le protocole de biomarqueurs sériques, l'élaboration du modèle diagnostique, la validation externe et la comparaison des modèles pronostiques sont résumés dans Figure 1.
Le flux des participants a inclus le dépistage, les exclusions, le regroupement diagnostique et l'affectation par centre. Au total, 326 individus ont été dépistés. Quarante-neuf individus ont été exclus, dont 12 ayant reçu un traitement antitumoral avant le prélèvement sanguin, 8 présentant une autre malignité primaire, 9 souffrant d'une infection grave, d'une maladie auto-immune ou de comorbidités non contrôlées susceptibles d'affecter les marqueurs tumoraux sériques, 11 présentant une interférence due à l'hémolyse, à la lipémie ou à la jaunisse, un volume sérique insuffisant ou un intervalle prolongé de traitement, et 9 ayant des données cliniques, pathologiques ou de suivi essentielles manquantes. La population finale destinée à l'analyse comprenait 277 participants : 167 patients atteints de CÉÉK, 55 témoins présentant une maladie œsophagienne bénigne et 55 témoins sains.
Raisonnement relatif à la taille de l'échantillon
La taille de l'échantillon a été déterminée par la disponibilité consécutive des participants éligibles, des échantillons sériques prétraitement et des informations de suivi pendant la période d'étude prédéfinie. Le modèle diagnostique utilisait trois biomarqueurs sériques continus comme prédicteurs candidats. Dans la cohorte d'apprentissage, 114 cas de CMOE et 80 témoins non CMOE étaient disponibles, soit 38,0 cas de CMOE et 26,7 témoins non CMOE par prédicteur biomarqueur candidat. Ceci permettait une régression logistique en faible dimension et réduisait le risque de surajustement. La cohorte de validation externe comprenait 53 cas de CMOE et 30 témoins non CMOE et a été utilisée uniquement pour la validation, sans réestimation des coefficients. La justification de la taille de l'échantillon était conforme aux principes actuels des modèles de prédiction, qui mettent l'accent sur la fréquence de l'issue, le nombre de paramètres prédicteurs candidats, la performance attendue du modèle et la validation indépendante, plutôt que sur la division aléatoire des données15.
Inscription des participants, critères d'éligibilité et collecte des données
Les participants ont été inclus de manière continue pendant la période d'étude. Dans le Centre A, les patients atteints de CEC ont été recrutés auprès des services de chirurgie thoracique, de gastro-entérologie et d'oncologie. Dans le Centre B, les patients atteints de CEC ont été identifiés par le biais de dépistages communautaires, de dossiers de suivi ou de voies de référence, et ont été inclus uniquement lorsque la confirmation histologique et des échantillons sériques préthérapeutiques étaient disponibles. Les témoins présentant des affections bénignes de l'œsophage ont été recrutés parmi les personnes ayant subi une gastroscopie pour dysphagie, symptômes de reflux ou anomalies œsophagiennes suspectées, et dont les résultats histologiques étaient non cancéreux. Les témoins sains ont été sélectionnés parmi les examens de santé de routine réalisés durant la même période et ne présentaient aucun antécédent de malignité, d'infection active, de dysfonctionnement hépatique ou rénal sévère, ni d'anomalies occupant un espace œsophagien visibles à l'imagerie ou à l'endoscopie.
Les patients atteints de CMOE étaient éligibles s'ils présentaient un CMOE primitif confirmé par biopsie gastroscope ou par examen anatomopathologique postopératoire, n'avaient reçu aucun traitement antitumoral avant le prélèvement sanguin, et disposaient de données cliniques de base complètes, de dossiers histopathologiques et d'échantillons sériques préthérapeutiques. Les participants étaient exclus s'ils présentaient une autre tumeur primitive, un traitement antitumoral antérieur au prélèvement sanguin, une infection aiguë sévère, une maladie auto-immune active, des comorbidités non contrôlées susceptibles d'affecter les taux sériques des marqueurs tumoraux, une hémolyse visible, une lipémie, une interférence due à l'ictère, un volume sérique insuffisant, des variables cliniques essentielles manquantes ou des données de suivi indisponibles. L'analyse pronostique était limitée aux patients atteints de CMOE pour lesquels le stade tumour-nodule-métastase (TNM) était disponible ainsi que des critères de survie ou de progression mesurables. Les patients dont le suivi était inférieur à 6 mois étaient conservés si un décès ou une progression survenait durant cette période.
Les variables de base comprenaient l'âge, le sexe, les antécédents de tabagisme, les antécédents de consommation d'alcool, l'indice de masse corporelle, le centre participant, le groupe diagnostique, le type de lésion et les taux sériques de CEA, de SCC-Ag et de CA125. Chez les patients atteints de CNEE, des variables supplémentaires comprenaient le site tumoral, la longueur de la tumeur, la différenciation histologique, la profondeur d'envahissement, le statut ganglionnaire, le stade TNM, l'envahissement lymphovasculaire, le mode de traitement initial, la récidive, la progression, le statut de survie et la date du dernier suivi. Les antécédents de tabagisme étaient définis comme un tabagisme actuel ou un tabagisme régulier antérieur d'au moins 100 cigarettes au cours de la vie. Les antécédents de consommation d'alcool étaient définis comme une consommation d'alcool d'au moins une fois par semaine pendant plus de 6 mois. Le stade TNM a été déterminé selon la huitième édition du système de stadification du American Joint Committee on Cancer/Union pour le contrôle international du cancer pour le cancer de l'œsophage16. Les coupes histologiques des cas de CNEE et de lésions bénignes ont été examinées indépendamment par deux pathologistes seniors, et toute divergence a été tranchée par un troisième pathologiste.
Deux chercheurs formés ont extrait indépendamment des informations cliniques à partir des dossiers médicaux électroniques, des comptes rendus d'anatomie pathologique, des systèmes d'information de laboratoire, des dossiers de traitement, des comptes rendus de suivi et des registres de déclaration des décès, en utilisant un dictionnaire de données prédéfini. Les divergences ont été résolues par vérification des documents d'origine, un troisième chercheur prenant la décision finale si nécessaire. L'emplacement de la tumeur a été classé comme œsophage thoracique supérieur, moyen ou inférieur. La différenciation histologique a été classée comme bien différenciée, modérément différenciée ou peu différenciée. Les unités des variables continues ont été standardisées entre les deux centres participants avant l'analyse.
Prélèvement sanguin, traitement du sérum et dosage des biomarqueurs
Tous les prélèvements sanguins ont été effectués avant le traitement initial, dans un délai de 7 jours suivant la confirmation histologique, et au moins 3 jours avant le début du traitement. Des échantillons de sang veineux périphérique de 5 à 8 mL ont été recueillis entre 7h00 et 9h00 du matin, après un jeûne nocturne, à l’aide de tubes avec séparation sérique. Le prélèvement était reporté si le participant avait reçu une transfusion sanguine, une perfusion d’albumine en grande quantité ou un geste invasif dans les 72 heures précédant le prélèvement. Les échantillons étaient laissés à coaguler à température ambiante pendant 30 minutes, puis centrifugés à 1 500 × g pendant 10 minutes. Le sérum était séparé et réparti en aliquotes dans 3 à 5 tubes portant une étiquette à code-barres, avec 300 à 500 µL par aliquote. Les aliquotes étaient conservées à −80 °C jusqu’à l’analyse, et aucune aliquote n’était soumise à plus d’un cycle de décongélation-congélation. Les échantillons présentant une hémolyse visible, une lipémie, une interférence due à l’ictère, un volume insuffisant ou un délai supérieur à 2 heures entre le prélèvement et la congélation étaient exclus.
Les taux sériques de CEA, de SCC-Ag et de CA125 ont été mesurés à l’aide de systèmes automatisés de dosage immuno-enzymatique par chimiluminescence, conformément aux instructions des fabricants. Les aliquotes sériques ont été décongelées à 4 °C avant l’analyse, doucement mélangées, puis analysées en une seule série, chaque fois que possible. Si un dosage par lots était nécessaire, l’ordre des échantillons était organisé de manière à éviter un regroupement systématique selon le statut de la maladie, le centre, le stade pathologique ou le résultat du suivi. Les échantillons dont les valeurs dépassaient la plage de mesure analytique ont été dilués et testés à nouveau conformément aux instructions du réactif. Le personnel de laboratoire était aveugle quant au groupe des participants, au stade pathologique, au statut de traitement et au résultat du suivi.
Les seuils cliniques étaient de 5,0 ng/mL pour l'ACE, de 1,5 ng/mL pour l'Ag-CEC et de 35 U/mL pour le CA125, selon les normes locales du laboratoire et les intervalles de référence des fabricants. Ces seuils ont été utilisés pour définir l'élévation des marqueurs et une charge élevée combinée de biomarqueurs dans les analyses d'association clinicopathologique et de stratification de la survie. Pour la classification diagnostique, des seuils déterminés par les données ont été choisis uniquement dans la cohorte d'apprentissage à l'aide de l'indice de Youden, puis appliqués sans modification à la cohorte de validation externe. L'étalonnage de l'appareil et le contrôle interne de qualité ont été effectués avant les analyses. Des contrôles à haut et à bas niveau ont été inclus chaque jour d'analyse, et les tests n'ont été poursuivis que lorsque les résultats du contrôle qualité se situaient dans la plage acceptée.
La reproductibilité inter-essais et inter-centres a été évaluée avant l'analyse finale du modèle à l'aide d'aliquots sériques de raccordement analysés en aveugle et de matériaux de contrôle qualité habituels. Un jeu de raccordement composé de 30 aliquots sériques analysés en aveugle a été testé selon les deux protocoles de laboratoire. Les coefficients de variation inter-centres étaient de 4,8 % pour l'ACE, de 5,6 % pour l'Ag-CCS et de 5,1 % pour le CA125. Les coefficients de corrélation intraclasse correspondants étaient de 0,94 pour l'ACE, de 0,92 pour l'Ag-CCS et de 0,93 pour le CA125. Aucune recalibration spécifique au centre n'a été effectuée, car aucun écart systématique n'a dépassé la plage prédéfinie de contrôle qualité du laboratoire. Les valeurs finales des biomarqueurs ont été harmonisées à l'aide d'une standardisation unitaire et d'une standardisation statistique basée sur le lot d'apprentissage avant la construction du modèle.
Définition combinée des biomarqueurs et construction du modèle
Le panel combiné de biomarqueurs sériques comprenait l'ACE, l'AgCC et le CA125. Pour l'analyse d'association clinicopathologique et la stratification de survie selon Kaplan-Meier, chaque marqueur était classé comme élevé ou non élevé selon le seuil clinique prédéfini. Une forte charge combinée de biomarqueurs était définie par l'élévation d'au moins deux des trois marqueurs. Pour la modélisation diagnostique et pronostique, les trois biomarqueurs ont été conservés comme variables continues, transformés à l'aide du logarithme népérien, puis normalisés à l'aide de la moyenne et de l'écart type du groupe d'apprentissage. Les mêmes paramètres de transformation ont été appliqués au groupe de validation externe.
Les biomarqueurs normalisés et transformés logarithmiquement ont été calculés à l'aide des paramètres suivants de la cohorte d'apprentissage :
zlnCEA = (lnCEA - 1,42)/0,71
zlnSCC - Ag = (lnSCC - Ag - 0,28)/0,57
zlnCA125 = (lnCA125 - 3,05)/0,62
Le panel diagnostique combiné a été construit à l'aide d'une régression logistique multivariée dans la cohorte d'apprentissage. Le prédicteur linéaire diagnostique a été calculé comme suit :
LPdiagnostic = -0,46 + 0,59 × zlnCEA + 1,05 × zlnSCC - Ag+0,42 × zlnCA125
La probabilité prédite individuelle de CECO a été calculée comme suit :
PESCC = 1/1+exp(-LPdiagnostic)
Le seuil diagnostique optimal pour le panel combiné était de 0,57, sélectionné dans la cohorte d'entraînement à l'aide de l'indice de Youden et appliqué sans modification lors de la validation externe.
Le score pronostique basé sur les biomarqueurs a été établi à l'aide d'un modèle de risques proportionnels de Cox. Le modèle reposant uniquement sur le stade TNM a servi de comparateur pour le stadage anatomique. Le modèle basé sur les biomarqueurs utilisait la charge combinée de biomarqueurs. Le modèle clinico-pathologique intégrait l'âge, la longueur de la tumeur, la différenciation histologique, le stade pT, le statut ganglionnaire et le mode de traitement initial. Le modèle intégré ajoutait une forte charge combinée de biomarqueurs au modèle clinico-pathologique. Pour la survie globale (OS), le prédicteur linéaire pronostique intégré a été calculé comme suit :
LPOS = 0,33 × âge ≥65 ans + 0,36 × longueur de la tumeur ≥5 cm + 0,25 × différenciation pauvre + 0,31 × pT3
- 4 étapes + 0,66 × métastase ganglionnaire
-0,58 × chirurgie radicale - traitement basé
+0,81 × highcombinedbiomarkerburden
Pour le SFP, le prédicteur linéaire pronostique intégré a été calculé comme suit :
LPPFS = 0,27 × âge ≥ 65 ans + 0,31 × longueur de la tumeur ≥ 5 cm + 0,24 × différenciation pauvre + 0,25 × pT3
-4 étapes + 0,55 × métastase ganglionnaire
- 0,46 × traitement basé sur la chirurgie radicale
+0,71 × highcombinedbiomarkerburden
Toutes les variables binaires dans les équations pronostiques ont été codées par 1 lorsqu'elles étaient présentes et par 0 lorsqu'elles étaient absentes. Tous les coefficients, paramètres de standardisation et seuils ont été estimés à partir de la cohorte d'apprentissage et appliqués directement à la cohorte de validation externe sans réajustement.
Paramètres diagnostiques et pronostiques
Le diagnostic pathologique a été utilisé comme critère de référence diagnostique. Le critère de jugement diagnostique principal était la discrimination entre le CÉRŒ et les témoins non CÉRŒ, incluant les maladies œsophagiennes bénignes et les témoins sains. Les critères de jugement diagnostiques secondaires comprenaient la discrimination entre le CÉRŒ et les maladies œsophagiennes bénignes, le CÉRŒ et les témoins sains, ainsi que le CÉRŒ précoce et les témoins non malignes. Le CÉRŒ précoce était défini comme une maladie de stade TNM I–II. La performance diagnostique a été évaluée à l'aide de courbes caractéristiques de fonctionnement du récepteur, de la surface sous la courbe, de la sensibilité, de la spécificité, de la valeur prédictive positive, de la valeur prédictive négative, de l'exactitude, du score F1, des courbes précision-rappel, des courbes de calibration, des scores de Brier, des intercepts de calibration, des pentes de calibration et de l'analyse de courbe de décision.
L'analyse pronostique a été réalisée uniquement chez les patients atteints de CMOE. Le critère de jugement pronostique principal était la survie globale (SG), définie comme l'intervalle entre la date du premier traitement antitumoral définitif et le décès pour toute cause. Le critère de jugement pronostique secondaire était la survie sans progression (SSP), définie comme l'intervalle entre le premier traitement et la progression documentée de la maladie, la récidive ou le décès. Chez les patients ayant bénéficié d'une chirurgie radicale, la récidive locale, la métastase à distance ou le décès lié à la tumeur étaient considérés comme un événement de SSP en l'absence de lésion résiduelle après la chirurgie. Chez les patients recevant une chimioradiothérapie définitive ou un traitement palliatif, la progression était déterminée selon les résultats d'imagerie, l'évaluation endoscopique et l'évaluation clinique.
Les données de suivi ont été obtenues par le biais de consultations externes, de dossiers d'hospitalisation, de dossiers de suivi en communauté, de dossiers de référence et de contacts téléphoniques. Les patients ont été suivis tous les 3 mois pendant les 2 premières années suivant le traitement, tous les 6 mois entre la 3e et la 5e année, puis annuellement par la suite. Les patients sans décès, récidive ou progression à la date finale de suivi ont été censurés à la date de leur dernier contact confirmé. Les événements relatifs aux résultats ont été examinés indépendamment par deux investigateurs. Lorsque les dates des événements différaient selon les sources, les rapports d'imagerie, les dossiers de sortie, les rapports endoscopiques, les documents de référence et les registres de déclaration de décès ont été consultés afin de déterminer la date finale de l'événement. La date limite de suivi était fixée au 31 décembre 2025 pour les deux centres participants.
Analyse statistique
L'analyse statistique a été réalisée à l'aide de R 4.3.2 et SPSS 27.0. Les variables continues ont été résumées par la moyenne ± l'écart type ou la médiane (intervalle interquartile), selon la distribution. Les variables catégorielles ont été résumées par des effectifs et des pourcentages. Les comparaisons entre groupes ont été effectuées à l'aide du test t pour échantillons indépendants, du test de Mann-Whitney U, de l'analyse de variance à un facteur, du test de Kruskal-Wallis, du test du chi-deux ou du test exact de Fisher, selon les cas. Les données manquantes ont été évaluées avant la construction des modèles. Les participants présentant des données manquantes concernant le statut diagnostique clé, les valeurs des biomarqueurs avant traitement, le stade TNM ou les données sur la survie/la progression ont été exclus conformément aux critères d'éligibilité prédéfinis. Les covariables non essentielles présentant moins de 5,0 % de données manquantes ont été traitées par imputation multiple par équations chaînées, avec 20 jeux de données imputées et 20 itérations17. Les modèles d'imputation comprenaient le groupe diagnostique, le centre, l'âge, le sexe, les antécédents de tabagisme, la consommation d'alcool, les valeurs des biomarqueurs, le stade TNM, le mode de traitement, le statut de survie globale, le statut de survie sans progression et la durée de suivi. Des analyses sur cas complets ont été menées comme analyses de sensibilité.
Le modèle diagnostique a été développé à partir de la cohorte d'apprentissage. Des courbes caractéristiques de fonctionnement du récepteur ont été utilisées pour calculer la surface sous la courbe, la sensibilité, la spécificité, la valeur prédictive positive, la valeur prédictive négative, l'exactitude et le score F1. Les surfaces sous la courbe ont été comparées à l'aide du test de DeLong18. Des courbes précision-rappel ont été utilisées comme évaluation diagnostique complémentaire. La calibration a été évaluée à l'aide de courbes de calibration, des scores de Brier, des intercepts de calibration et des pentes de calibration. L'utilité clinique a été évaluée par une analyse de courbe de décision19. Une validation interne a été réalisée à l'aide de 1 000 rééchantillonnages bootstrap afin d'estimer la discrimination et la calibration corrigées de l'optimisme. Une validation externe a été effectuée en appliquant directement à partir de la cohorte d'apprentissage l'équation, les paramètres de standardisation et le seuil au Centre B, sans réajustement.
Pour l'analyse pronostique, des courbes de Kaplan-Meier ont été générées afin de comparer la survie globale (OS) et la survie sans progression (PFS) entre les groupes selon la charge de biomarqueurs, et le test du log-rank a été utilisé pour évaluer les différences en matière de survie. Une régression de Cox aux risques proportionnels a été utilisée pour estimer les rapports de risques instantanés (hazard ratios) et les intervalles de confiance à 95 %. Les variables présentant une pertinence clinique ainsi que celles associées aux résultats dans l'analyse univariée ont été prises en compte pour les modèles multivariés. La performance des modèles a été évaluée à l'aide de l'indice C de Harrell, des courbes caractéristiques de fonctionnement du récepteur (ROC) dépendantes du temps, des courbes de calibration, des scores de Brier et de l'analyse de courbe décisionnelle20. La calibration pronostique a été évaluée à 1, 3 et 5 ans. L'hypothèse de proportionnalité des risques a été vérifiée à l'aide des résidus de Schoenfeld et des graphiques de survie log-moins-log21. Les valeurs P globales du test de Schoenfeld étaient de 0,42 pour le modèle de survie globale et de 0,37 pour le modèle de survie sans progression. La multicolinéarité a été évaluée à l'aide des facteurs d'inflation de la variance avant l'interprétation finale du modèle, et toutes les variables candidates retenues dans les modèles diagnostiques et pronostiques finaux présentaient des facteurs d'inflation de la variance inférieurs à 2,1. La valeur pronostique ajoutée par rapport au stade TNM et aux modèles clinicopathologiques conventionnels a été évaluée à l'aide de l'amélioration nette de la reclassification (net reclassification improvement) et de l'amélioration de discrimination intégrée (integrated discrimination improvement). Les paramètres des modèles estimés dans la cohorte d'apprentissage ont été appliqués à la cohorte de validation externe sans réajustement22. Tous les tests étaient bilatéraux, et une valeur de P < 0,05 a été considérée comme statistiquement significative.