Данное исследование было рассмотрено и одобрено комитетом по этике при аффилированной первой народной больнице Ханчжоу медицинской школы Университета Уэстлэйк. Номер разрешения — IIT-20230528-0108-01. До забора крови от участников, включенных в исследование проспективно, было получено письменное информированное согласие. Для клинических записей, включенных ретроспективно, и архивных остаточных образцов, собранных в ходе стандартного медицинского обслуживания, комитет по этике освободил исследователей от необходимости получения информированного согласия, так как использовались обезличенные данные и никаких дополнительных вмешательств не проводилось.
Воспроизводимость, контроль систематических ошибок и целостность данных
Критерии соответствия, временные точки забора крови, правила экспертной оценки визуализационных данных, определения подгрупп, потенциальные предикторы, определения пороговых значений, правила обработки пропущенных данных и процедуры валидации были заранее определены до проведения статистического анализа. Образцы крови собирали до подтверждения результатов визуализации с контрастным усилением и до начала любого противоопухолевого лечения. Лабораторные исследования проводились с использованием анонимизированных идентификаторов исследования, а образцы из обучающей и валидационной когорт тестировались в рандомизированном порядке. Лабораторный персонал не был осведомлен о окончательном диагнозе, принадлежности к когорте, стадии по классификации Barcelona Clinic Liver Cancer (BCLC) и статусе AFP-отрицательного анализа. Диагностика по данным визуализации и стадирование проводились независимо от лабораторных исследований. Ввод данных, преобразование единиц измерения, распределение по когортам, кодирование исключений, кодирование пропущенных значений и статистические наборы данных проверялись двумя исследователями перед подбором модели. Обучающая и валидационная когорты оставались разделенными на всех этапах отбора переменных, подбора модели, вывода пороговых значений и валидации.
Популяция и дизайн исследования
В данном диагностическом исследовании с двумя когортами была разработана и валидирована многофакторная модель для выявления ГЦК с использованием показателей аномального протромбина/белка, индуцированного PIVKA-II, АФП и стандартных клинических переменных. Ретроспективно сформированная обучающая когорта включала участников, обследованных в период с 1 января 2021 г. по 31 декабря 2023 г. Эта когорта использовалась для скрининга переменных, подбора модели, определения пороговых значений и внутренней валидации методом бутстрепа. Проспективно набранная валидационная когорта включала участников, обследованных с 1 января 2024 г. по 31 декабря 2025 г. Эта когорта использовалась только для тестирования на независимой выборке. В валидационной когорте не проводилось повторного отбора переменных, переподбора модели, обновления коэффициентов, повторной оптимизации пороговых значений или перекалибровки на основе подгрупп. Поток участников, распределение по когортам, забор крови, обработка плазмы, тестирование биомаркеров, экспертная оценка снимков и статистический рабочий процесс показаны на рисунке 1.
Участники были набраны из группы наблюдения или диагностики заболеваний печени с высоким риском в Первой народной больнице города Ханчжоу, при медицинском факультете Университета Уэстлейк. Статус высокого риска определялся как хроническая инфекция гепатита B, хроническая инфекция гепатита C, цирроз любой этиологии или жировая болезнь печени с подтвержденным фиброзом на основании клинической оценки, лабораторных данных или результатов предшествующей визуализации. Критериями включения были взрослые пациенты в возрасте ≥18 лет, которые прошли забор крови перед визуализацией, а также мультифазную компьютерную томографию (КТ) с контрастным усилением и/или динамическую магнитно-резонансную томографию (МРТ) с контрастным усилением в течение 14 дней после забора крови. Критериями исключения были: предшествующее лечение ГЦК, включая резекцию, абляцию, трансартериальную или системную терапию; введение витамина K в течение 7 дней до забора проб; терапия антагонистами витамина K или другими антикоагулянтными препаратами, которые могли повлиять на интерпретацию коагуляции; острый холангит или обструкция желчных путей на момент забора проб; беременность; отсутствие данных по AFP, PIVKA-II или окончательного диагноза по результатам визуализации; а также неприемлемое качество образцов при отсутствии возможности повторного забора.
Итоговая исследуемая популяция включала 389 участников. Обучающая когорта состояла из 205 участников, включая 112 случаев ГЦК, подтвержденных методами визуализации, и 93 контроля без ГЦК. Валидационная когорта включала 184 участников, включая 58 случаев ГЦК, подтвержденных методами визуализации, и 126 контролей без ГЦК. В обучающей когорте 93 контроля без ГЦК включали 30 здоровых добровольцев и 63 пациента с хроническими заболеваниями печени. В валидационной когорте 126 контролей без ГЦК включали 0 здоровых добровольцев и 126 пациентов с хроническими заболеваниями печени.
Симптоматический статус не определялся только на основании истории заболевания гепатитом или компенсированного цирроза. Участник классифицировался как симптоматический только в том случае, если в записи перед проведением визуализации были задокументированы новые или усиливающиеся боли в правом верхнем квадранте живота, необъяснимая потеря веса, вздутие живота, желтуха, темная моча или клинически подтвержденная декомпенсация, не объясненная другой установленной причиной. Участники без этих признаков классифицировались как бессимптомные или выявленные в ходе скрининга.
Диагностическое заключение и стадирование
Окончательный статус ГЦК определяли с помощью мультифазной КТ с контрастным усилением и/или динамической МРТ с контрастным усилением, выполненных в рамках стандартного медицинского обслуживания. Визуализация должна была включать как минимум артериальную и портальную венозную фазы; данные отсроченной фазы учитывались при их наличии. Очаговое образование классифицировалось как ГЦК, если при визуализации печени в группе риска наблюдалось гиперусиление в артериальную фазу с последующим «вымыванием» контраста в портальную венозную и/или отсроченную фазы, с наличием или отсутствием контрастируемой капсулы. Участники, у которых отсутствовали признаки ГЦК по данным визуализации и не было других злокачественных диагнозов, классифицировались как группа без ГЦК.
Информация о визуализации была извлечена с помощью заранее определенной формы, в которой фиксировались метод визуализации, полученные фазы, гиперусиление в артериальную фазу, вымывание контраста, вид капсулы, количество очагов, максимальный диаметр очага, макрососудистая инвазия и внепеченочное распространение. Два обученных специалиста независимо извлекали характеристики визуализации из окончательных радиологических отчетов. Расхождения разрешались путем консенсуса с использованием тех же предопределенных критериев. Если были доступны и КТ, и МРТ, использовался метод с наиболее полным многофазным описанием. Если оба метода были представлены в полном объеме, использовался тот метод, который был задокументирован как основание для окончательного клинического решения.
Стадия по системе BCLC была определена после окончательного подтверждения статуса ГЦК. ГЦК на ранней стадии определялась как стадия BCLC 0–A. ГЦК с отрицательным уровнем АФП определялась как АФП <20 ng/mL при заборе крови перед визуализацией. Эти метки подгрупп использовались для оценки эффективности в подгруппах, а не для обучения модели или определения соответствия участников критериям включения.
Сбор образцов и обработка плазмы
Периферическую венозную кровь собирали утром после ночного голодания продолжительностью не менее 8 h, до подтверждения диагноза с помощью КТ/МРТ с контрастным усилением и до начала противоопухолевой терапии. Всего 10 mL венозной крови собирали в пробирки с антикоагулянтом K2-EDTA. Сразу после сбора пробирки перемешивали путем переворачивания 8–10 раз и транспортировали в лабораторию при 4 °C. Интервал от венепункции до центрифугирования составлял не более 2 h и фиксировался в журнале образцов.
Плазму отделяли путем центрифугирования при 1,600 × g в течение 10 мин при 4 °C с использованием центрифуги с охлаждением. Супернатант плазмы переносили в микроцентрифужные пробирки объемом 1.5 mL, не содержащие нуклеаз, и разделяли на порции по 0.5–1.0 mL. Аликвоты хранили при -80 °C до проведения анализа. Перед измерением уровня биомаркеров одну аликвоту однократно размораживали при комнатной температуре, 20–25 °C, в течение 20–30 мин, осторожно перемешивали путем переворачивания и тестировали в течение 2 ч после полного размораживания. Образцы с видимым свертыванием, выраженным или умеренным гемолизом или выраженной липемией отбраковывали в соответствии с заранее определенными критериями приемлемости лаборатории. Для когорты проспективной валидации, когда это было клинически возможно, проводили повторный забор крови перед подтверждением данными визуализации.
Измерение уровней AFP и PIVKA-II
Уровни AFP и PIVKA-II измеряли с помощью электрохемилюминесцентного иммуноанализа на автоматическом иммуноанализаторе. AFP регистрировали в ng/mL, а PIVKA-II в mAU/mL. Результаты экспортировали непосредственно из лабораторной информационной системы.
Каждая аналитическая серия включала специфические для данного анализа калибраторы и как минимум два уровня внутреннего контроля качества. Для приемки серии требовалось, чтобы все результаты контроля качества находились в пределах заранее определенного лабораторного контрольного диапазона. Целевые показатели аналитической прецизионности составляли: коэффициент вариации внутри серии ≤10% и коэффициент вариации между сериями ≤15% как для AFP, так и для PIVKA-II. В случае неудовлетворительного контроля качества или обнаружения дрейфа серии, серия отклонялась, анализатор проходил повторную калибровку, и серия повторялась при наличии достаточного объема образца. При получении неправдоподобного или выходящего за пределы диапазона результата для отдельного образца, по возможности выполнялись два повторных измерения из той же размороженной аликвоты, и использовалось среднее значение двух принятых повторов. Если повторное тестирование не дало результата или объем образца был недостаточен, результат помечался как недостоверный и обрабатывался в соответствии с правилом обработки пропущенных данных.
Рутинное извлечение клинических переменных
Стандартные клинические показатели были извлечены из данных того же клинического визита, во время которого проводился забор крови на AFP/PIVKA-II перед визуализацией. В число рассматриваемых переменных вошли возраст, пол, аланинаминотрансфераза, аспартатаминотрансфераза, альбумин, общий билирубин, международное нормализованное отношение (INR), протромбиновое время и количество тромбоцитов. При наличии нескольких значений в окне перед визуализацией выбиралось значение, ближайшее по времени к моменту забора крови на AFP/PIVKA-II. Перед анализом единицы измерения были стандартизированы. Возраст регистрировался в годах; пол кодировался как женский = 0 и мужской = 1; количество тромбоцитов — в ×10⁹/L; альбумин — в g/L; общий билирубин — в μmol/L; AFP — в ng/mL; PIVKA-II — в mAU/mL.
После извлечения данных были проведены проверки диапазонов, единиц измерения и поиск неправдоподобных значений. Отмеченные значения сверялись с исходными записями. Значение исключалось только в том случае, если в исходной записи подтверждалась проблема с качеством образца, ошибка преобразования единиц или отклонение при обработке. Верифицированные экстремальные значения были сохранены. Словарь переменных определял структуру набора данных, названия переменных, единицы измерения, правила кодирования, правила преобразования и коды отсутствующих значений.
Определение предикторов и выбор переменных
Потенциальные предикторы были определены до разработки модели. Основными биомаркерами были AFP и PIVKA-II. В число стандартных клинических предикторов вошли возраст, пол, аланинаминотрансфераза, аспартатаминотрансфераза, альбумин, общий билирубин, МНО или протромбиновое время, а также количество тромбоцитов. Для AFP и PIVKA-II было выполнено натуральное логарифмическое преобразование. Значения ниже предела отчетности анализа перед преобразованием были заменены на половину нижней границы предела отчетности.
Для одномерного сравнения участников с ГЦК и без ГЦК в обучающей когорте использовались t-критерии для независимых выборок для переменных с приблизительно нормальным распределением, U-критерии Манна-Уитни для непрерывных переменных с асимметричным распределением, а также критерии χ2 или точный критерий Фишера для категориальных переменных. Переменные со значением P < 0.10 при одномерном скрининге рассматривались для включения в многофакторное моделирование. Мультиколлинеарность оценивали с помощью фактора инфляции дисперсии (VIF); значение VIF > 5 считалось признаком значимой коллинеарности. В случае сильной коллинеарности двух переменных сохранялась та, которая имела более четкую клиническую интерпретируемость и обеспечивала лучшую приростную эффективность модели.
Окончательная модель многофакторной логистической регрессии была разработана только для обучающей когорты. В модель были принудительно включены ln(AFP) и ln(PIVKA-II). Дополнительные предикторы отбирались методом обратного исключения на основе информационного критерия Акаике. После обратного отбора в модель вошли следующие переменные: ln(AFP), ln(PIVKA-II), возраст, пол, уровень альбумина, МНО и количество тромбоцитов.
Разработка модели и расчет показателя риска
Диагностическая модель была построена с использованием многофакторной логистической регрессии, где бинарным исходом являлся статус ГЦК. Прогнозируемая вероятность ГЦК рассчитывалась на основе линейного предиктора η по следующему уравнению:
p(HCC)=1/[1+exp(-η)]
η = -4.862+0.247×ln(AFP)+0.712×ln(PIVKA-II)+0.018×age+0.331×sex-0.062×albumin+0.554×INR-0.004×plateletcount
Концентрацию AFP измеряли в ng/mL, PIVKA-II в mAU/mL, возраст в годах, пол кодировали как женщина = 0 и мужчина = 1, уровень альбумина измеряли в g/L, INR был безразмерным, а количество тромбоцитов измеряли в ×109/L. Для AFP и PIVKA-II использовались натуральные логарифмы. На основании тех же подогнанных коэффициентов была построена номограмма. Для анализа рабочих характеристик приемника (ROC-анализа), анализа калибровки, анализа кривой принятия решений и валидации использовалась прогнозируемая вероятность, а не балльная шкала номограммы.
Определения пороговых значений и обработка пределов
Пороговые значения PIVKA-II были обозначены в соответствии с аналитической целью. Установленный клинический порог для PIVKA-II составил 40.0 mAU/mL. Порог по индексу Юдена для PIVKA-II в обучающей когорте составил 45.0 mAU/mL. Порог с фиксированной специфичностью был получен путем выбора значения, наиболее близкого к 90% специфичности в обучающей когорте; для PIVKA-II это значение составило 38.0 mAU/mL. Видимый порог только для валидации, равный 37.5 mAU/mL, рассматривался как исследовательский и не использовался в качестве основного порога валидации.
Для комбинированной модели основной бинарный порог был получен на основе обучающей когорты с использованием индекса Юдена и применен без изменений к валидационной когорте. Для определения подгруппы AFP-отрицательных пациентов и сравнения со стандартной клинической практикой использовался установленный клинический порог AFP, составляющий 20.0 ng/mL. Оптимизированные пороги AFP, использованные при ROC-анализе, обозначались отдельно от клинического порога. Таким образом, значения порогов указывались в зависимости от аналитической цели, а не рассматривались как взаимозаменяемые показатели.
Производительность, калибровка и валидация модели
Дискриминационную способность оценивали с помощью кривых рабочих характеристик приемника (ROC-кривых) и площади под кривой (AUC). Значения AUC приводились с 95% доверительными интервалами. Для сравнения AUC использовали тест Делонга. Эффективность классификации оценивали по чувствительности, специфичности, прогностической значимости положительного результата, прогностической значимости отрицательного результата и точности при заранее заданных значениях отсечки.
Калибровку оценивали по точке пересечения калибровочной кривой, наклону калибровочной кривой, показателю Брайера и критерию согласия Хосмера–Лемешоу. На графиках калибровки сравнивали прогнозируемую и наблюдаемую вероятности ГЦК в различных группах риска. Внутренняя валидация проводилась с использованием 1 000 бутстрэп-ресэмплов в обучающей когорте для оценки скорректированного на оптимизм AUC и скорректированного на оптимизм наклона калибровочной кривой. Валидационная когорта использовалась только для внешней оценки эффективности. В валидационной когорте не проводилось обновление коэффициентов, переподгонка модели, повторная оптимизация порога или рекалибровка на основе подгрупп.
Клиническая значимость была оценена с помощью анализа кривых принятия решений при пороговых вероятностях от 0.05 до 0.50. Сравнение чистой выгоды проводилось для AFP отдельно, PIVKA-II отдельно, модели на основе клинических переменных и комбинированной модели.
Обработка пропущенных данных и анализ чувствительности
Для основного анализа использовалось моделирование на основе полных случаев для участников, по которым были сохранены все переменные в окончательной модели. Для каждого потенциального предиктора был проведен обобщающий анализ пропущенных значений. Если доля пропусков для любого сохраненного предиктора превышала 5%, в качестве анализа чувствительности выполнялась многократная импутация методом итерационных уравнений (MICE), при этом в модель импутации включались исходный показатель, все сохраненные предикторы и индикатор когорты. Переменные, доля пропусков в которых превышала 20%, не допускались к многофакторному моделированию. Коэффициенты модели, AUC, калибровочный наклон и показатели классификации, полученные в результате анализа с импутацией, сравнивались с результатами анализа на основе полных случаев.
Программное обеспечение и воспроизводимая отчетность
Все анализы были выполнены с использованием R версии 4.3.2 и RStudio версии 2023.12.1. Логистическая регрессия была проведена с помощью базового пакета stats. ROC-анализ и сравнения по методу DeLong были выполнены с использованием pROC версии 1.18.5. Анализ калибровки проводился с использованием rms версии 6.7-1 и ResourceSelection версии 0.3-6. Для бутстреп-валидации использовалось 1 000 повторных выборок с фиксированным случайным числом (seed) 70558. Анализ кривых принятия решений был выполнен с помощью rmda версии 1.6. Множественная импутация, если она требовалась, выполнялась с помощью mice версии 3.16.0. Скрипт анализа воспроизводил схему потока участников, коэффициенты моделей, ROC-кривые, показатели калибровки, бутстреп-валидацию, анализ кривых принятия решений, анализы на основе пороговых значений и статистические таблицы.