Данное исследование с участием людей, использованием клинических записей и образцов сыворотки крови проводилось в соответствии с Хельсинкской декларацией и требованиями учреждения по проведению исследований с участием людей13. Протокол исследования был рассмотрен и одобрен больницей традиционной китайской медицины провинции Цзянсу (Jiangsu Province Hospital of Chinese Medicine) под номером одобрения CJ28J12-45. Для ретроспективно собранных анонимизированных клинических данных и остаточных образцов сыворотки крови требование о получении дополнительного письменного информированного согласия было отменено в соответствии с правилами институциональной этики. Для проспективно собранных образцов письменное информированное согласие было получено перед включением в исследование. Вся информация об участниках и образцы сыворотки были анонимизированы перед анализом. Для связывания данных использовались идентификаторы исследования; персонально идентифицируемая информация в аналитический набор данных не включалась.
Дизайн исследования и распределение когорт
В данном двухцентровом амбиспективном обсервационном исследовании оценивалась диагностическая и прогностическая значимость комбинированной панели сывороточных биомаркеров, включающей раково-эмбриональный антиген (CEA), антиген плоскоклеточного рака (SCC-Ag) и углеводный антиген 125 (CA125), при плоскоклеточном раке пищевода (ESCC). Ретроспективные данные были получены из соответствующих клинических записей и остаточных образцов сыворотки крови, взятых до начала лечения; проспективный набор участников использовался для сбора дополнительных дотерапевтических образцов сыворотки и завершения последующего наблюдения в рамках того же исследования. Участники набирались в одной многопрофильной больнице и одном общественном центре здравоохранения в период с января 2021 года по декабрь 2024 года, при этом наблюдение продолжалось до 31 декабря 2025 года. Центр А (многопрофильная больница) служил обучающей когортой для разработки модели, оценки коэффициентов, выбора пороговых значений и внутренней валидации. Центр Б (общественный центр здравоохранения) служил когортой для внешней валидации. Отчет об исследовании был составлен в соответствии с рекомендациями Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD)14.
Всего было включено 277 участников, из них 194 из центра А и 83 из центра B. Обучающая когорта включала 114 пациентов с ESCC, 40 человек из группы контроля с доброкачественными заболеваниями пищевода и 40 здоровых добровольцев. Когорта внешней валидации включала 53 пациента с ESCC, 15 человек из группы контроля с доброкачественными заболеваниями пищевода и 15 здоровых добровольцев. Вместо случайного разделения использовалось распределение по центрам, чтобы уменьшить утечку данных и проверить стабильность модели в независимых клинических условиях. В рамках диагностического анализа сравнивались уровни CEA, SCC-Ag и CA125 в сыворотке крови пациентов с ESCC, пациентов с доброкачественными заболеваниями пищевода и здоровых добровольцев. Прогностический анализ проводился только для пациентов с ESCC с целью оценки того, улучшает ли комбинированная нагрузка биомаркеров стратификацию риска общей выживаемости (OS) и выживаемости без прогрессирования (PFS). Дизайн исследования, скрининг участников, распределение по когортам, рабочий процесс анализа биомаркеров сыворотки, разработка диагностической модели, внешняя валидация и сравнение прогностических моделей обобщены на Рисунке 1.
Поток участников включал скрининг, исключение, диагностическую группировку и распределение по центрам. Всего был проведен скрининг 326 человек. Сорок девять человек были исключены, в том числе 12 — из-за проведения противоопухолевой терапии до забора крови, 8 — из-за наличия другого первичного злокачественного новообразования, 9 — из-за тяжелой инфекции, аутоиммунного заболевания или неконтролируемых сопутствующих патологий, которые могли повлиять на уровень сывороточных опухолевых маркеров, 11 — из-за гемолиза, липемии, интерференции при желтухе, недостаточного объема сыворотки или длительного интервала обработки, и 9 — из-за отсутствия ключевых клинических, патологических данных или данных последующего наблюдения. Итоговая аналитическая выборка составила 277 участников: 167 пациентов с плоскоклеточным раком пищевода (ESCC), 55 человек из группы контроля с доброкачественными заболеваниями пищевода и 55 здоровых добровольцев из группы контроля.
Обоснование размера выборки
Размер выборки определялся последовательной доступностью подходящих участников, образцов сыворотки до начала лечения и данных последующего наблюдения в течение заранее определенного периода исследования. В диагностической модели в качестве потенциальных предикторов использовались три непрерывных серологических биомаркера. В обучающей когорте были доступны 114 случаев ESCC и 80 контрольных образцов без ESCC, что соответствует 38,0 случаям ESCC и 26,7 контрольным образцам без ESCC на один биомаркер-предиктор. Это обеспечило возможность применения низкоразмерной логистической регрессии и снизило риск переобучения модели. Когорта внешней валидации включала 53 случая ESCC и 30 контрольных образцов без ESCC и использовалась только для валидации без переоценки коэффициентов. Обоснование размера выборки соответствовало современным принципам построения прогностических моделей, в которых приоритет отдается частоте исходов, количеству параметров предикторов, ожидаемой эффективности модели и независимой валидации, а не случайному разделению данных15.
Набор участников, критерии соответствия и сбор данных
Участники набирались непрерывно в течение периода исследования. В центре А пациенты с плоскоклеточным раком пищевода (ПКРП) привлекались из отделений торакальной хирургии, гастроэнтерологии и онкологии. В центре Б пациенты с ПКРП выявлялись с помощью скрининга на уровне местных сообществ, записей последующего наблюдения или по направлениям и включались в исследование только при наличии патологического подтверждения и образцов сыворотки крови, взятых до начала лечения. В качестве контроля с доброкачественными заболеваниями пищевода привлекались лица, которым проводилась гастроскопия по поводу дисфагии, симптомов рефлюкса или подозрений на патологии пищевода, и у которых по результатам патоморфологического исследования не было выявлено злокачественных новообразований. Здоровые добровольцы отбирались из числа лиц, проходивших плановые медицинские осмотры в тот же период; у них отсутствовал анамнез злокачественных заболеваний, активных инфекций, тяжелой печеночной или почечной недостаточности, а также данных визуализации или эндоскопии, свидетельствующих о наличии объемных образований в пищеводе.
Пациенты с плоскоклеточным раком пищевода (ПКРП) считались подходящими для исследования, если у них был подтвержден первичный ПКРП по результатам гастроскопической биопсии или послеоперационного патоморфологического исследования, отсутствовало противоопухолевое лечение до забора крови, а также имелись полные исходные клинические данные, патологические записи и образцы сыворотки до начала лечения. Участники исключались, если у них имелась другая первичная злокачественная опухоль, проводилось противоопухолевое лечение до забора крови, наблюдалась тяжелая острая инфекция, активное аутоиммунное заболевание, неконтролируемые сопутствующие заболевания, которые могли повлиять на уровни сывороточных опухолевых маркеров, а также при наличии видимого гемолиза, липемии, иктеричности, недостаточного объема сыворотки, отсутствия ключевых клинических переменных или недоступности данных последующего наблюдения. Прогностический анализ был ограничен пациентами с ПКРП, для которых были доступны стадия по системе TNM (опухоль-узел-метастаз) и измеримые показатели выживаемости или прогрессирования. Пациенты с периодом наблюдения менее 6 месяцев были включены в анализ, если смерть или прогрессирование наступили в течение этого периода.
Исходные переменные включали возраст, пол, стаж курения, историю употребления алкоголя, индекс массы тела, центр участия, диагностическую группу, тип поражения, а также уровни CEA, SCC-Ag и CA125 в сыворотке крови. Для пациентов с ESCC дополнительными переменными были локализация опухоли, длина опухоли, гистологическая дифференцировка, глубина инвазии, состояние лимфатических узлов, стадия по системе TNM, лимфоваскулярная инвазия, метод первичного лечения, рецидив, прогрессирование, статус выживаемости и дата последнего последующего наблюдения. Стаж курения определялся как текущее курение или регулярное курение в прошлом (не менее 100 сигарет за всю жизнь). История употребления алкоголя определялась как прием алкоголя не менее одного раза в неделю в течение более 6 месяцев. Стадирование по TNM определялось в соответствии с восьмым изданием системы стадирования рака пищевода Объединенного комитета по раку США/Международного агентства по борьбе с раком16. Патологоанатомические стекла случаев ESCC и доброкачественных поражений были независимо изучены двумя ведущими патологоанатомами, а любые разногласия разрешались третьим патологоанатомом.
Двое квалифицированных исследователей независимо извлекали клиническую информацию из электронных медицинских карт, патологоанатомических отчетов, лабораторных информационных систем, записей о лечении, данных последующего наблюдения и свидетельств о смерти, используя заранее определенный словарь данных. Расхождения устранялись путем сверки с исходными записями, при этом в случае необходимости окончательное решение принимал третий исследователь. Локализацию опухоли разделяли на верхний, средний или нижний отделы грудного отдела пищевода. Гистологическую дифференцировку классифицировали как высокую, умеренную или низкую. Единицы измерения непрерывных переменных были стандартизированы в обоих участвующих центрах перед началом анализа.
Забор крови, обработка сыворотки и измерение уровней биомаркеров
Все образцы крови собирали до начала первичного лечения, в течение 7 дней после патологоанатомического подтверждения и как минимум за 3 дня до начала терапии. Образцы периферической веночной крови объемом 5–8 mL собирали с 7:00 до 9:00 утра после ночного голодания с использованием пробирок для разделения сыворотки. Забор крови откладывали, если участнику проводили переливание крови, инфузию альбумина в большом объеме или инвазивную процедуру в течение 72 h до сбора. Образцам давали свернуться при комнатной температуре в течение 30 min, затем центрифугировали при 1,500 × g в течение 10 min. Сыворотку отделяли и распределяли по 3–5 пробиркам с штрихкодами, по 300–500 µL на одну аликвоту. Аликвоты хранили при −80°C до проведения анализа; ни одна аликвота не подвергалась более чем одному циклу замораживания-оттаивания. Образцы с видимым гемолизом, липемией, желтушной интерференцией, недостаточным объемом или интервалом от сбора до замораживания более 2 h исключали из исследования.
Уровни РЭА, SCC-Ag и CA125 в сыворотке крови измеряли с помощью автоматизированных систем хемилюминесцентного иммуноанализа в соответствии с инструкциями производителей. Аликвоты сыворотки размораживали при 4°C перед тестированием, осторожно перемешивали и по возможности анализировали в одном цикле. Если требовалось серийное тестирование, порядок образцов определяли таким образом, чтобы избежать систематической группировки по статусу заболевания, центру, патологической стадии или исходу наблюдения. Образцы, значения которых превышали диапазон аналитического измерения, разбавляли и тестировали повторно согласно инструкциям к реагентам. Лабораторный персонал не был осведомлен о группе участников, патологической стадии, статусе лечения и исходе наблюдения.
Клинические пороговые значения составили 5,0 ng/mL для CEA, 1,5 ng/mL для SCC-Ag и 35 U/mL для CA125, согласно местным лабораторным стандартам и референсным интервалам производителей. Эти пороги использовались для определения повышения уровня маркеров и высокой комбинированной нагрузки биомаркеров при проведении анализа клинико-патологической корреляции и стратификации выживаемости. Для диагностической классификации значения отсечки, основанные на данных, подбирались только в обучающей когорте с использованием индекса Юдена, а затем применялись без изменений к внешней валидационной когорте. Перед началом тестирования проводились калибровка приборов и внутренний контроль качества. В каждый день тестирования включались контроли высокого и низкого уровней, и анализы проводились только в том случае, если результаты контроля качества находились в допустимых пределах.
Воспроизводимость между анализами и между центрами была оценена перед окончательным анализом модели с использованием заслепленных связующих аликвот сыворотки и стандартных материалов контроля качества. Связующий набор из 30 заслепленных аликвот сыворотки был протестирован в двух рабочих процессах лабораторий. Коэффициенты вариации между центрами составили 4,8% для CEA, 5,6% для SCC-Ag и 5,1% для CA125. Соответствующие коэффициенты внутриклассовой корреляции составили 0,94 для CEA, 0,92 для SCC-Ag и 0,93 для CA125. Перекалибровка анализа для конкретных центров не проводилась, так как ни одно систематическое отклонение не превысило заранее определенный диапазон лабораторного контроля качества. Окончательные значения биомаркеров были гармонизированы путем стандартизации единиц измерения и статистической стандартизации на основе обучающей когорты перед построением модели.
Определение комбинированного биомаркера и построение модели
Комбинированная панель сывороточных биомаркеров включала CEA, SCC-Ag и CA125. Для анализа клинико-патологических взаимосвязей и стратификации выживаемости по методу Каплана-Мейера каждый маркер классифицировали как повышенный или не повышенный в соответствии с заранее определенным клиническим порогом. Высокую комбинированную нагрузку биомаркеров определяли как повышение уровней как минимум двух из трех маркеров. Для диагностического и прогностического моделирования три биомаркера использовали в качестве непрерывных переменных, преобразованных с помощью натурального логарифма и стандартизированных с использованием среднего значения и стандартного отклонения обучающей когорты. Те же параметры преобразования применяли к когорте внешней валидации.
Стандартизированные логарифмированные биомаркеры были рассчитаны с использованием следующих параметров обучающей когорты:
zlnCEA = (lnCEA - 1.42)/0.71
zlnSCC - Ag = (lnSCC - Ag - 0.28)/0.57
zlnCA125 = (lnCA125 - 3.05)/0.62
Комбинированная диагностическая панель была сформирована с использованием многофакторной логистической регрессии в обучающей когорте. Диагностический линейный предиктор рассчитывался следующим образом:
LPdiagnostic = -0,46 + 0,59 × zlnCEA + 1,05 × zlnSCC - Ag+0,42 × zlnCA125
Индивидуальная прогнозируемая вероятность ESCC рассчитывалась следующим образом:
PESCC = 1/1+exp(-LPdiagnostic)
Оптимальным диагностическим пороговым значением для комбинированной панели было 0,57; оно было определено в обучающей когорте с использованием индекса Юдена и применялось без изменений при внешней валидации.
Прогностическая оценка биомаркеров была построена с использованием модели пропорциональных рисков Кокса. Модель, основанная только на стадиях TNM, использовалась в качестве эталона для анатомического стадирования. В модели биомаркеров использовалась комбинированная нагрузка биомаркеров. В клинико-патологической модели учитывались возраст, размер опухоли, гистологическая дифференцировка, стадия pT, статус лимфатических узлов и метод первичного лечения. Интегрированная модель дополняла клинико-патологическую модель показателем высокой комбинированной нагрузки биомаркеров. Для OS интегрированный прогностический линейный предиктор рассчитывался следующим образом:
LPOS = 0,33 × возраст ≥65 лет + 0,36 × длина опухоли ≥5см + 0,25 × низкая степень дифференцировки + 0,31 × pT3
- 4стадия + 0,66 × метастазы в лимфатических узлах
-0,58 × лечение на основе радикальной хирургии
+0,81 × highcombinedbiomarkerburden
Для PFS интегрированный прогностический линейный предиктор рассчитывали следующим образом:
LPPFS = 0,27 × возраст ≥ 65 лет + 0,31 × длина опухоли ≥ 5 см + 0,24 × низкая степень дифференцировки + 0,25 × pT3
-4стадия + 0,55 × метастазы в лимфоузлах
- 0,46 × лечение на основе радикальной хирургии
+0,71 × highcombinedbiomarkerburden
Все бинарные переменные в прогностических уравнениях были закодированы как 1 при наличии признака и 0 при его отсутствии. Все коэффициенты, параметры стандартизации и пороговые значения были рассчитаны на обучающей когорте и применены непосредственно к когорте внешней валидации без повторной оптимизации.
Диагностические и прогностические конечные точки
В качестве диагностического эталона использовался патологоанатомический диагноз. Первичной диагностической конечной точкой была дифференциация между плоскоклеточным раком пищевода (ESCC) и контрольной группой без ESCC, включая доброкачественные заболевания пищевода и здоровых добровольцев. Вторичные диагностические конечные точки включали дифференциацию между ESCC и доброкачественными заболеваниями пищевода, ESCC и здоровыми добровольцами, а также между ESCC на ранней стадии и контрольной группой без злокачественных новообразований. ESCC на ранней стадии определялся как заболевание стадии I–II по системе TNM. Диагностическая эффективность оценивалась с помощью ROC-кривых, площади под кривой, чувствительности, специфичности, прогностической ценности положительного результата, прогностической ценности отрицательного результата, точности, F1-меры, кривых точности-полноты, калибровочных кривых, показателей Брайера, интерцептов калибровки, наклонов калибровки и анализа кривых принятия решений.
Прогностический анализ проводился только среди пациентов с ESCC. Основной прогностической конечной точкой была общая выживаемость (OS), определенная как интервал от даты первого окончательного противоопухолевого лечения до смерти от любой причины. Вторичной прогностической конечной точкой была выживаемость без прогрессирования (PFS), определенная как интервал от первого лечения до задокументированного прогрессирования заболевания, рецидива или смерти. Для пациентов, перенесших радикальную операцию, местный рецидив, отдаленный метастаз или смерть от опухоли расценивались как событие PFS при отсутствии остаточных очагов после операции. Для пациентов, получавших окончательную химиолучевую терапию или паллиативное лечение, прогрессирование определялось на основании данных визуализации, эндоскопической оценки и клинического обследования.
Данные последующего наблюдения были получены в ходе амбулаторных визитов, из записей о госпитализации, отчетов о наблюдении в местных сообществах, документов о направлении к специалистам и посредством телефонной связи. Пациенты проходили обследование каждые 3 месяца в течение первых 2 лет после лечения, каждые 6 месяцев с 3-го по 5-й год и далее ежегодно. Пациенты, у которых к дате последнего визита не было зафиксировано смерти, рецидива или прогрессирования заболевания, были цензурированы по дате последнего подтвержденного контакта. События исхода независимо оценивались двумя исследователями. В случае расхождения дат событий в разных источниках, для определения окончательной даты события проверялись отчеты о визуализации, выписные эпикризы, протоколы эндоскопии, документы о направлении и записи о регистрации смерти. Датой окончания наблюдения для обоих участвующих центров было 31 декабря 2025 года.
Статистический анализ
Статистический анализ проводился с использованием R 4.3.2 и SPSS 27.0. Непрерывные переменные представлялись в виде среднего значения ± стандартного отклонения или медианы (межквартильный размах) в зависимости от распределения. Категориальные переменные представлялись в виде абсолютного количества и процентов. Межгрупповые сравнения проводились с использованием t-критерия для независимых выборок, U-критерия Манна–Уитни, однофакторного дисперсионного анализа, критерия Краскела–Уоллиса, критерия хи-квадрат или точного критерия Фишера, в зависимости от ситуации. Оценка пропущенных данных проводилась перед построением модели. Участники, у которых отсутствовали основные данные о диагностическом статусе, значениях биомаркеров до начала лечения, стадии по системе TNM или данные о выживаемости/прогрессировании, были исключены в соответствии с заранее определенными критериями приемлемости. Для второстепенных ковариат с долей пропусков менее 5.0% применялся метод множественной импутации с помощью цепочек уравнений с 20 импутированными наборами данных и 20 итерациями17. Модели импутации включали диагностическую группу, центр, возраст, пол, историю курения, употребление алкоголя, значения биомаркеров, стадию TNM, метод лечения, статус OS, статус PFS и время последующего наблюдения. Анализы с использованием только полных случаев (complete-case analyses) проводились в качестве анализа чувствительности.
Диагностическая модель была разработана в обучающей когорте. Для расчета площади под кривой, чувствительности, специфичности, прогностической ценности положительного результата, прогностической ценности отрицательного результата, точности и показателя F1 использовались ROC-кривые (рабочие характеристики приемника). Сравнение площадей под кривыми проводилось с помощью теста DeLong18. В качестве дополнительной диагностической оценки использовались кривые точности-полноты. Калибровка оценивалась с помощью калибровочных кривых, показателей Брайера, интерцептов и наклонов калибровки. Клиническая полезность оценивалась с помощью анализа кривых принятия решений19. Внутренняя валидация выполнялась с использованием 1 000 бутстрэп-выборок для оценки дискриминации и калибровки с поправкой на оптимизм. Внешняя валидация проводилась путем прямого применения уравнения обучающей когорты, параметров стандартизации и порогового значения к данным Центра B без повторной подгонки модели.
Для прогностического анализа были построены кривые Каплана — Мейера для сравнения OS и PFS между группами с различной нагрузкой биомаркеров, а для оценки различий в выживаемости использовался логранговый критерий. Регрессия пропорциональных рисков Кокса применялась для оценки отношений рисков и 95% доверительных интервалов. Для многофакторных моделей рассматривались переменные, имеющие клиническую значимость, а также переменные, связанные с исходами в однофакторном анализе. Эффективность модели оценивали с помощью C-индекса Харрелла, зависящих от времени ROC-кривых, калибровочных кривых, показателей Брайера и анализа кривых принятия решений20. Прогностическая калибровка оценивалась через 1, 3 и 5 лет. Допущение о пропорциональности рисков проверялось с помощью остатков Шенфельда и графиков log-minus-log выживаемости21. Значения P глобального теста Шенфельда составили 0,42 для модели OS и 0,37 для модели PFS. Перед окончательной интерпретацией модели оценивали мультиколлинеарность с помощью факторов инфляции дисперсии; для всех переменных-кандидатов в итоговых диагностических и прогностических моделях факторы инфляции дисперсии были ниже 2,1. Дополнительная прогностическая ценность по сравнению со стадией TNM и традиционными клинико-патологическими моделями оценивалась с помощью показателей чистого улучшения реклассификации и интегрированного улучшения дискриминации. Параметры модели, рассчитанные в обучающей когорте, применялись к внешней валидационной когорте без повторной подгонки22. Все тесты были двусторонними, и значение P < 0,05 считалось статистически значимым.