Предложенная структура XAI была реализована с использованием набора данных Pima Indians Diabetes Dataset в качестве репрезентативного набора данных из сферы здравоохранения. Pima Indians Diabetes Dataset использовался как единственный набор данных для экспериментальной валидации. Все представленные результаты прогнозирования, объяснимости, статистические данные и результаты воспроизводимости были получены на основе этого набора данных. TCGA-BRCA, TCGA-UCEC, MIMIC-III, ISIC 2019, HAM10, OhioT1DM и BraTS 2021 не подвергались экспериментальной оценке и включены только в качестве примеров, иллюстрирующих применимость общего протокола к различным модальностям медицинских данных. Полный набор данных использовался после удаления недействительных и дублирующих записей, а перед разработкой модели были выполнены указанные операции по предварительной обработке. Набор данных был разделен на независимые обучающую, валидационную и тестовую подвыборки с использованием предопределенной стратегии стратифицированного разделения. Независимость образцов между тремя подвыборками поддерживалась для минимизации возможности утечки данных.
Прогностическая модель была настроена с использованием предопределенной архитектуры и гиперпараметров. Модель обучали с помощью оптимизатора Adam с заданными скоростью обучения и размером пакета в течение максимум 100 эпох. Была применена процедура ранней остановки на основе потерь на валидационной выборке, и была сохранена модель с наилучшими показателями валидации. Для разделения набора данных, инициализации модели и проведения повторных экспериментов были заданы случайные значения (random seeds) для обеспечения воспроизводимости результатов.
Обученная модель была оценена на независимой тестовой выборке с использованием таких показателей, как точность (accuracy), прецизионность (precision), полнота (recall), специфичность, F1-мера, коэффициент корреляции Мэтьюза (MCC), площадь под ROC-кривой (AUC-ROC) и средняя точность (AP). Предложенная модель сравнивали с заранее определенными базовыми моделями, используя идентичные процедуры предобработки, разбиение данных и протоколы оценки. На основе прогнозов для независимой тестовой выборки были построены ROC-кривые (кривые рабочих характеристик приемника), кривые прецизионности-полноты и матрица ошибок.
Для оценки стабильности работы была выполнена пятикратная перекрестная проверка на обучающих данных. Были рассчитаны 95% доверительные интервалы для основных показателей эффективности, а также проведено заранее определенное статистическое сравнение предлагаемой модели с базовыми моделями.
Пятикратная перекрестная проверка показала точность 93.01 ± 0.48%, AUC-ROC 0.954 ± 0.07, прецизионность 92.42 ± 0.87%, полноту 93.02 ± 0.45% и F1-меру 92.72 ± 0.62%. 95% бутстреп-доверительные интервалы, рассчитанные на основе 1 0 ресэмплов, составили 0.897–0.973 для точности, 0.890–0.970 для прецизионности, 0.80–0.963 для полноты, 0.87–0.965 для F1-меры и 0.931–0.984 для AUC-ROC. Статистическое сравнение с базовыми моделями CNN, Random Forest и SVM проводилось с использованием теста Макнемара для оценки точности, теста Делонга для AUC-ROC и парного бутстреп-тестирования с 1 0 ресэмплами для F1-меры.
Полная процедура обучения и оценки была повторена в 10 независимых запусках с использованием различных случайных значений seed. Повторные запуски обеспечили AUC-ROC 0,957 ± 0,008, точность 93,24 ± 0,74% и F1-меру 92,35 ± 0,92%, что указывает на относительно низкую вариабельность между повторными исполнениями. Показатели эффективности, полученные в ходе повторных запусков, были обобщены с использованием среднего значения и стандартного отклонения. Вариабельность между запусками была изучена, чтобы определить, остается ли прогностическая эффективность стабильной при многократном исполнении.
В данном практическом примере внешняя валидация не проводилась, так как независимый внешний набор данных не использовался. Соответственно, все представленные результаты прогнозирования, статистические данные и результаты воспроизводимости были получены на основе набора данных Pima Indians Diabetes с использованием предопределенных разделений на обучающую, валидационную и независимую тестовую выборки. Внешняя валидация оставлена в протоколе в качестве дополнительной процедуры для случаев, когда доступен независимый набор данных из другого учреждения, популяции, географического региона или за другой временной период.
Пояснения к моделям были созданы с использованием методов интерпретируемости, применимых к табличному набору данных Pima Indians Diabetes Dataset, включая SHAP и LIME. Была проведена оценка глобальной значимости признаков, а также сформированы локальные пояснения для конкретных пациентов с использованием отложенных тестовых выборок. Результаты интерпретации были зафиксированы вместе с соответствующими предсказаниями моделей и значениями признаков для обеспечения воспроизводимости и последующего анализа.
Для ясности данный практический пример включал девять основных этапов рабочего процесса, указанных в Таблице 1. Внешняя валидация и оценка клиническими экспертами были оставлены в качестве дополнительных модулей валидации общего протокола. Внешняя валидация не проводилась, так как независимый внешний набор данных не использовался, а оценка клиническими экспертами не проводилась, поскольку в данный практический пример не была включена официальная экспертная комиссия. Соответственно, эти дополнительные модули не считаются частью девятиэтапного экспериментального рабочего процесса и не представлены в качестве экспериментальных результатов.
Процедуры предварительной обработки и разделения данных позволили получить стандартизированный набор данных, пригодный для разработки модели. Дубликаты и противоречивые записи были удалены, пропущенные значения обработаны в соответствии с заранее определенной стратегией, числовые признаки стандартизированы, а набор данных разделен на независимые обучающую, валидационную и тестовую выборки. Характеристики результирующего набора данных и процедуры предварительной обработки обобщены в Таблице 2. Общий рабочий процесс подготовки и предварительной обработки набора данных здравоохранения проиллюстрирован на Рисунке 2, тогда как рабочий процесс экспериментальной подготовки, предварительной обработки, разделения и оценки качества, примененный конкретно к набору данных Pima Indians Diabetes Dataset, показан на Рисунке 3. Итоговая вычислительная среда, архитектура модели и конфигурация гиперпараметров, использованные для получения представленных результатов, обобщены в Таблице 3.
Выполнение разделов 3 и 4 позволило получить стандартизированный набор данных здравоохранения, пригодный для разработки модели. В ходе процедур предобработки были удалены дубликаты и противоречивые записи, заполнены пропущенные значения, стандартизированы числовые признаки, закодированы категориальные переменные (где это было применимо), а также произведено разделение набора данных на обучающую, валидационную и тестовую выборки. Полученные данные обеспечили стандартизированный ввод, необходимый для последующей разработки модели.
После завершения Разделов 5 и 6 сконфигурированная модель продемонстрировала стабильную сходимость в процессе обучения. Кривые обучения показали одновременное снижение потерь при обучении и валидации, а также рост точности при обучении и валидации. Оптимизация гиперпараметров улучшила обобщающую способность модели, при этом признаков существенного переобучения выявлено не было. Для обеспечения воспроизводимости оптимизированная модель была архивирована вместе с окончательной архитектурой, настройками гиперпараметров, версиями программного обеспечения, значениями случайных чисел (seed) и весами обученной модели. Спецификации обучения модели и результаты оптимизации обобщены в Таблице 4, а соответствующие кривые обучения при обучении и валидации представлены на Рисунке 4.
В разделе 7 оценивалась прогностическая эффективность модели с использованием стандартизированных показателей эффективности. Оцениваемые метрики классификации включали точность (accuracy), прецизионность (precision), полноту (recall), специфичность (specificity), F1-меру (F1-score), коэффициент корреляции Мэтьюза (MCC), AUC-ROC и среднюю точность (AP). Предложенная модель сравнивалась с заранее определенными базовыми моделями с использованием тех же разделений наборов данных, процедур предобработки и протокола оценки. Сравнение прогностической эффективности представлено в Таблице 5, в то время как ROC-кривые, кривые прецизионности-полноты, матрица ошибок и сравнительные показатели эффективности показаны на Рисунке 5.
В продолжение Раздела 8 были созданы как глобальные, так и локальные объяснения предсказаний модели для оценки ее интерпретируемости. Объяснения модели были созданы с помощью SHAP и LIME для табличного набора данных Pima Indians Diabetes Dataset, а для иллюстрации изменения предсказания было создано специфичное для конкретного пациента контрфактуальное объяснение. SHAP использовался для генерации глобальной значимости признаков, а также визуализаций в виде диаграмм waterfall и зависимостей признаков для конкретных пациентов, в то время как LIME применялся для создания локальных объяснений на основе отложенных тестовых выборок. Соответствующие результаты интерпретируемости представлены на Рисунке 6.
На заключительном этапе протокола оценивались статистическая достоверность и воспроизводимость рабочего процесса. Полный рабочий процесс был повторен в 10 независимых запусках с использованием различных случайных чисел (random seeds) при сохранении той же стратегии разделения набора данных, параметров предобработки, архитектуры модели, гиперпараметров, программной среды и процедур оценки. Повторные запуски дали значение AUC-ROC 0.957 ± 0.08, точность (accuracy) 93.24 ± 0.74% и F1-меру 92.35 ± 0.92%, что указывает на относительно низкую вариабельность между повторными итерациями.
В настоящей валидационной работе количественная оценка стабильности объяснений не проводилась. Несмотря на то, что для набора данных Pima Indians Diabetes Dataset были созданы объяснения с помощью SHAP и LIME, отдельный анализ ранговой корреляции или перекрытия признаков между запусками не выполнялся. Следовательно, численные показатели стабильности объяснений или согласованности атрибуции не приводятся. Количественная оценка стабильности объяснений сохранена в общем протоколе в качестве необязательной процедуры проверки воспроизводимости для тех случаев, когда доступны результаты многократного получения объяснений.
Статистические сравнения оценивали с использованием заранее определенного порога значимости p < 0.05. В соответствующих случаях применялись двусторонние статистические тесты; для количественной оценки статистической значимости и неопределенности наблюдаемых различий в эффективности приводились соответствующие статистики теста, значения p и 95% доверительные интервалы. Результаты экспериментальной статистической валидации и воспроизводимости, включая показатели перекрестной проверки, доверительные интервалы, статистические сравнения и вариабельность при повторных запусках, обобщены в Таблице 6. Соответствующий статистический анализ и анализ воспроизводимости представлены на Рисунке 7.
Эти результаты предоставили экспериментальные доказательства прогностической стабильности и воспроизводимости при протестированных вычислительных условиях и наборе данных. Вычислительная среда, характеристики набора данных, процедуры предобработки, конфигурация модели, статистический анализ и настройки интерпретируемости, необходимые для независимого воспроизведения, были задокументированы в соответствии с контрольным списком воспроизводимости, представленным в Таблице 7. Экспертная клиническая оценка сгенерированных результатов XAI в приведенном в данной работе примере валидации не проводилась.
В целом, применение данного протокола позволило создать стандартизированный набор данных здравоохранения, подходящий для разработки моделей ИИ, и оптимизированную модель с использованием предопределенных настроек гиперпараметров. Описанный рабочий процесс обеспечил систематическую оценку прогностической способности, генерацию объяснений модели с помощью соответствующих методов XAI, а также статистическую оценку устойчивости и воспроизводимости модели. В совокупности результаты продемонстрировали возможность реализации и воспроизводимости предложенного рабочего процесса XAI в экспериментальных условиях, рассмотренных в примере валидации.

Рисунок 1: Девятиэтапный основной рабочий процесс разработки воспроизводимых и объяснимых моделей ИИ в здравоохранении.Рабочий процесс включает следующие этапы: (1) определение задачи прогнозирования в здравоохранении, (2) конфигурация вычислительной среды, (3) сбор и валидация набора данных, (4) предобработка данных, (5) разделение набора данных, (6) обучение прогностической модели, (7) оценка прогностической эффективности, (8) анализ объяснимости и (9) статистическая валидация и оценка воспроизводимости. Внешняя валидация и оценка клиническими экспертами рассматриваются как дополнительные расширения протокола и не включены в основной девятиэтапный рабочий процесс. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 2: Рабочий процесс подготовки и предобработки набора данных для здравоохранения. (A) Сбор данных здравоохранения из клинических записей, медицинских изображений, физиологических сигналов и мультимодальных источников данных. (B) Интеграция и предобработка данных, включая обработку пропущенных значений, нормализацию, удаление шума и аугментацию. (C) Разделение набора данных на обучающую, валидационную и тестовую выборки. (D) Оценка качества, демонстрирующая распределение классов, нормализацию признаков и результаты предобработки перед разработкой модели. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 3: Экспериментальный рабочий процесс подготовки, предварительной обработки, разделения и оценки качества набора данных Pima Indians Diabetes Dataset. Рабочий процесс включает получение набора данных, оценку качества данных, обработку недопустимых и пропущенных значений, стандартизацию численных признаков, разделение набора данных на обучающую, валидационную и независимую тестовую подвыборки, а также окончательную проверку качества перед разработкой модели. Нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 4: Кривые обучения при тренировке и валидации предлагаемой модели с использованием набора данных Pima Indians Diabetes Dataset. (A) Потери при обучении и валидации на протяжении всего периода обучения. (B) Точность при обучении и валидации на протяжении всего периода обучения. (C) Увеличенный вид потерь в течение эпох 50–10. (D) Увеличенный вид точности в течение эпох 50–10. Лучшие показатели валидации были достигнуты на 78-й эпохе с потерями при валидации 0.142 и точностью валидации 94.6%. Ранняя остановка была инициирована на 8-й эпохе при значении patience, равном 10 эпохам. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 5: Экспериментальная оценка прогностической эффективности предлагаемой структуры XAI с использованием независимой тестовой выборки (n = 154). (A) ROC-кривая (кривая рабочих характеристик приемника), демонстрирующая дискриминационную способность предлагаемой модели XAI и базовых моделей. (B) Кривые точности-полноты (PR-кривые), демонстрирующие показатели точности и полноты предлагаемой модели XAI и базовых моделей. (C) Матрица ошибок предлагаемой модели XAI на независимой тестовой выборке с соответствующими значениями точности, чувствительности (полноты) и специфичности. (D) Сравнение точности, прецизионности, полноты, специфичности, F1-меры, коэффициента корреляции Мэтьюза (MCC), площади под ROC-кривой (AUC) и средней точности (AP) для всех оцениваемых моделей. Все количественные результаты, представленные на этом рисунке, относятся к валидационному эксперименту наборе данных по диабету индейцев племени Пима (Pima Indians Diabetes Dataset). Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 6: Результаты интерпретируемости, полученные на основе прогнозов предложенной модели на независимом тестовом наборе данных Pima Indians Diabetes Dataset. (A) Глобальный сводный график SHAP, показывающий распределение вкладов признаков во всем тестовом наборе. (B) График важности признаков SHAP, основанный на средних абсолютных значениях SHAP. (C) Индивидуальный для пациента водопадный график SHAP, показывающий вклад отдельных признаков в прогнозируемую вероятность диабета. (D) Локальное объяснение LIME, показывающее вклад признаков для тестового пациента №125. (E) График зависимости SHAP, показывающий взаимосвязь между значениями глюкозы и их вкладами SHAP. (F) Индивидуальное для пациента контрфактуальное объяснение, показывающее минимальные изменения признаков, связанные с изменением прогноза модели. Сокращения: SHAP = Shapley Additive exPlanations; LIME = Local Interpretable Model-agnostic Explanations. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.

Рисунок 7: Экспериментальная статистическая валидация и анализ воспроизводимости предлагаемой модели с использованием набора данных Pima Indians Diabetes Dataset. (A) Результаты пятикратной перекрестной проверки на обучающей выборке. (B) 95% бутстреп-доверительные интервалы для показателей эффективности на независимой тестовой выборке. (C) Статистическое сравнение с базовыми моделями, включая используемый статистический критерий, значение статистики теста, p-value и уровень значимости. (D) Стабильность работы модели в 10 независимых запусках с различными случайными значениями seed. Для сравнения парной точности классификации использовался критерий Макнемара, для коррелированного ROC-AUC — тест Делона, а для сравнения F1-меры — парный бутстреп-тест с 1 0 повторных выборок. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
| Этап | Действие по протоколу | Ожидаемый результат | Статус реализации |
| 1 | Выбор и валидация набора данных здравоохранения | Проверенный набор данных, цель прогнозирования, распределение классов и информация о качестве данных | Выполнено |
| 2 | Настройка вычислительной среды | Проверенное оборудование, программное обеспечение, библиотеки, версии и конфигурация вычислений | Выполнено |
| 3 | Предварительная обработка и контроль качества данных здравоохранения | Очищенный, преобразованный и стандартизированный набор данных | Выполнено |
| 4 | Разбиение набора данных | Независимые обучающие, валидационные и тестовые наборы данных | Выполнено |
| 5 | Разработка и оптимизация прогностической модели / модели XAI | Окончательная архитектура модели и гиперпараметры | Выполнено |
| 6 | Обучение и сохранение модели | Обученная модель, контрольная точка, конфигурация обучения и логи | Выполнено |
| 7 | Оценка прогностической и вычислительной эффективности | Метрики эффективности на тестовом наборе и сравнение производительности | Выполнено |
| 8 | Генерация и оценка результатов интерпретируемости | Результаты SHAP/LIME и соответствующие пояснения | Выполнено |
| 9 | Статистическая валидация и оценка воспроизводимости | Доверительные интервалы, статистические тесты, результаты повторных запусков и показатели воспроизводимости | Выполнено |
Таблица 1: Резюме основного девятиэтапного рабочего процесса протокола, примененного в ходе валидации с использованием набора данных по диабету индейцев племени Пима (Pima Indians Diabetes Dataset). В таблице девять этапов, реализованных в практическом примере с использованием набора данных Pima Indians Diabetes Dataset, отделены от внешней валидации и оценки клиническими экспертами, которые остаются дополнительными компонентами общего протокола.
| Набор данных | Источник данных | Клиническое применение | Модальность данных | Испытуемые/Записи | Образцы | Классы | Распределение классов | Обучающая/валидационная/тестовая выборки | Роль в настоящем исследовании | Статус валидации | URL-адрес источника |
| Набор данных по диабету у индейцев племени пима | Репозиторий машинного обучения UCI | Прогнозирование диабета | Клиническая таблица | 768 записей | 768 | 2 | 50 лиц без сахарного диабета; 268 лиц с сахарным диабетом | 60% / 20% / 20% | Основной набор данных экспериментальной валидации | Выполнено — основной рабочий процесс из девяти этапов завершен | База данных диабета индейцев племени Пима (Pima Indians Diabetes Database) |
| TCGA-BRCA | Центр общих геномных данных (GDC) Национального института рака (NCI), проект TCGA-BRCA | Классификация рака молочной железы | Клиническая и гистопатологическая характеристика | 1 098 случаев | Зависимость от набора данных в зависимости от типа данных | Зависящий от конечной точки | Отсутствие единого распределения классов по всему набору данных | Неприменимо — экспериментальное разделение не проводилось | Пример только для иллюстрации применимости протокола | Не использовалось для экспериментальной валидации | https://portal.gdc.cancer.gov/projects/TCGA-BRCA |
| TCGA-UCEC | NCI Genomic Data Commons (GDC), проект TCGA-UCEC | Классификация рака эндометрия | Клиническая и гистопатологическая характеристика | Когорта проекта; размер зависит от выбранного типа данных и фильтров | Зависимость набора данных от типа данных | Зависящий от конечной точки | Отсутствие единого распределения классов по всему набору данных | Неприменимо — разделение на экспериментальные группы не проводилось | Пример применимости протокола | Не используется для экспериментальной валидации | https://portal.gdc.cancer.gov/projects/TCGA-UCEC |
| MIMIC-III | PhysioNet, клиническая база данных MIMIC-III v1.4 | Прогнозирование клинических исходов | Клинические временные ряды | 46 520 пациентов | 58 976 госпитализаций в отделения интенсивной терапии | Зависящий от задачи | Отсутствие единого распределения классов по всей базе данных | Неприменимо — разделение на экспериментальные группы не проводилось | Пример применимости протокола | Не используется для экспериментальной валидации | https://physionet.org/content/mimiciii/1.4/ |
| ISIC 2019 | Конкурс International Skin Imaging Collaboration (ISIC) | Классификация кожных поражений | Дерматоскопические изображения | Неприменимо – набор изображений | 25 331 обучающее изображение | 8 категорий обучения | АКИЕЦ 867; БК 323; БКЛ 2624; ДФ 239; МЕЛ 452; НВ 12875; ВАСК 253; ПлК 628 | Неприменимо — разделение на экспериментальные группы не проводилось | Пример только для области применения протокола | Не использовалось для экспериментальной валидации | https://challenge.isic-archive.com/landing/2019/ |
| HAM100 | Harvard Dataverse / Архив ISIC | Классификация поражений кожи | Дермоскопические изображения | 7 470 уникальных поражений | 10 015 изображений | 7 | AKIEC 327; БК 514; БКЛ 1 09; ДФ 115; МЕЛ 1 13; НВ 6 705; ВАСК 142 | Неприменимо — разделение на экспериментальные группы не проводилось | Пример только для применимости протокола | Не используется для экспериментальной валидации | **Оценка влияния ионизирующего излучения на жизнеспособность клеток с использованием анализа MTT**
**Абстракт**
В данной работе представлен подробный протокол проведения анализа MTT (3-(4,5-диметилтиазол-2-ил)-2,5-дифенилтетразолиум бромид) для оценки влияния ионизирующего излучения на жизнеспособность клеток. Анализ MTT является широко используемым колориметрическим методом определения метаболической активности клеток, который служит индикатором их жизнеспособности и пролиферации. В этом эксперименте клетки подвергаются воздействию различных доз ионизирующего излучения, после чего их жизнеспособность оценивается по способности митохондриальных дегидрогеназ восстанавливать тетразолиевую соль MTT до образующихся кристаллов формазана. В протоколе подробно описываются этапы подготовки клеток, процесс облучения, процедура окрашивания MTT и последующее измерение оптической плотности с использованием планшетного ридера. Результаты выражаются в процентах жизнеспособности клеток по сравнению с необлученной контрольной группой. Данный метод предоставляет количественную основу для изучения радиочувствительности различных клеточных линий и оценки эффективности радиопротекторов или радиосенсибилизаторов.
**Ключевые слова**
ионизирующее излучение, жизнеспособность клеток, анализ MTT, радиобиология, метаболическая активность, цитоксичность, пролиферация клеток, митохондриальная активность |
| OhioT1DM | набор данных OhioT1DM, распространяемый в открытом доступе для исследовательских целей | Мониторинг и прогнозирование диабета | Клинические временные ряды | 12 участников | 24 XML-файла, распределенных по обучающей, проверочной и тестовой выборкам | Непрерывное прогнозирование уровня глюкозы; не является задачей фиксированной классификации | Неприменимо | Файлы для обучения/разработки и тестирования, определенные набором данных; экспериментальное разделение здесь не выполнялось | Пример применимости протокола | Не используется для экспериментальной валидации | **Методы анализа экспрессии генов и белков при изучении нейровоспаления**
**Аннотация**
Нейровоспаление характеризуется активацией микроглии и астроцитов, что приводит к высвобождению провоспалительных цитокинов и хемокинов. Понимание молекулярных механизмов, лежащих в основе этих процессов, имеет решающее значение для разработки терапевтических стратегий при нейродегенеративных заболеваниях и травмах центральной нервной системы. В данной статье представлены и подробно описываются общепринятые методы анализа экспрессии генов и белков, используемые для изучения нейровоспаления. Мы рассматриваем количественную полимеразную цепную реакцию с обратной транскрипцией (RT-qPCR) для анализа мРНК, вестерн-блоттинг для определения уровней белков и иммуногистохимическое окрашивание для визуализации пространственной локализации маркеров воспаления в тканях мозга. В работе приводятся детальные протоколы, рекомендации по выбору контрольных групп и методы статистической обработки данных, что позволяет исследователям точно количественно определять и локализовать воспалительные ответы в нейрональных моделях.
**Ключевые слова**
нейровоспаление, микроглия, астроциты, RT-qPCR, вестерн-блоттинг, иммуногистохимия, экспрессия генов, маркеры воспаления
**Введение**
Воспаление в центральной нервной системе (ЦНС) является сложным динамическим процессом, который может играть как защитную, так и патологическую роль. Основными эффекторными клетками нейровоспаления являются микроглия (резидентные макрофаги мозга) и астроциты. При воздействии патогенов, повреждений или белковых агрегатов эти клетки переходят в активированное состояние, изменяя свой профиль экспрессии генов и секретируя широкий спектр медиаторов, включая интерлейкины (например, IL-1β, IL-6), фактор некроза опухоли (TNF-α) и хемокины (например, MCP-1).
Для изучения этих изменений исследователи полагаются на комплексный подход, сочетающий анализ транскрипции и трансляции. Анализ мРНК с помощью RT-qPCR позволяет обнаружить ранние изменения в экспрессии генов и высокую чувствительность к низким уровням транскриптов. Вестерн-блоттинг обеспечивает подтверждение наличия и количества функциональных белков. Иммуногистохимический анализ (ИГХ) критически важен для определения того, какие именно типы клеток экспрессируют данные маркеры и как они распределены в различных анатомических областях мозга.
**Материалы и методы**
**1. Анализ экспрессии мРНК с помощью RT-qPCR**
**Приготовление РНК и синтез кДНК**
1. Выделение общей РНК из тканей мозга или культур клеток проводили с использованием набора RNeasy Mini Kit (Qiagen) в соответствии с инструкциями производителя.
2. Концентрацию и чистоту РНК определяли с помощью спектрофотометра NanoDrop.
3. Синтез первой цепи кДНК осуществляли из 1 мкг общей РНК с использованием High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems).
**Количественная ПЦР в реальном времени**
1. Реакции проводили с использованием SYBR Green PCR Master Mix.
2. Специфические праймеры были разработаны для целевых генов нейровоспаления (например, *Tnf*, *Il1b*, *Il6*, *Gfap*, *Iba1*).
3. В качестве внутреннего контроля (референсного гена) использовали *Gapdh* или *Actb*.
4. Амплификацию проводили в приборе StepOnePlus Real-Time PCR System.
5. Относительный уровень экспрессии генов рассчитывали методом $2^{-\Delta\Delta\text{Ct}}$.
**2. Анализ уровней белков методом вестерн-блоттинга**
**Экстракция белков и электрофорез**
1. Ткани мозга гомогенизировали в лизирующем буфере RIPA с добавлением ингибиторов протеаз и фосфатаз.
2. Общую концентрацию белка определяли методом Брэдфорд.
3. Равные количества белка (20–40 мкг) разделяли с помощью SDS-PAGE (10% или 12% полиакриламидный гель).
**Перенос и иммунодетекция**
1. Белки переносили на PVDF-мембрану с помощью электропереноса.
2. Мембраны блокировали 5% раствором БСА или обезжиренного сухого молока в ТБС-Т (TBS-T) в течение 1 часа при комнатной температуре.
3. Инкубацию с первичными антителами (например, против Iba1, GFAP, NF-κB) проводили в течение ночи при 4°C.
4. После промывки мембраны инкубировали с соответствующими антителами вторичного уровня, конъюгированными с пероксидазой хрена (HRP).
5. Хемилюминесцентный сигнал визуализировали с помощью системы ECL и анализировали с помощью программного обеспечения ImageJ для денситометрии.
**3. Иммуногистохимический анализ (ИГХ)**
**Фиксация и подготовка срезов**
1. Мозг фиксировали перфузией 4% раствором параформальдегида (PFA), после чего ткани погружали в 4% PFA на 24 часа.
2. Срезы толщиной 30–40 мкм готовили с помощью криостата.
**Процедура окрашивания**
1. Срезы инкубировали в растворе Триттона Х-100 (0,3% в PBS) для обеспечения проницаемости мембран.
2. Блокировали нормальной сывороткой соответствующего вида (например, козьей) для минимизации неспецифического связывания.
3. Инкубировали с первичными антителами (например, anti-Iba1 для микроглии, anti-GFAP для астроцитов) в течение 48 часов при 4°C.
4. Специфическое связывание визуализировали с помощью флуоресцентно-меченых вторичных антител или метода иммунопероксидазы с использованием DAB (диаминобензидина).
5. Срезы монтировали и анализировали с помощью флуоресцентного или светового микроскопа.
**Результаты и обсуждение**
Применение данных методов позволяет детально картировать нейровоспалительный ответ. RT-qPCR выявляет быстрый подъем уровней мРНК провоспалительных цитокинов, что свидетельствует об активации транскрипционных факторов, таких как NF-κB. Вестерн-блоттинг подтверждает, что эти транскрипты транслируются в белки, которые затем секретируются или выполняют структурные функции. ИГХ-анализ позволяет отличить общую активацию глии от локализованного ответа в конкретных областях (например, в гиппокампе или коре) и оценить морфологические изменения, такие как утолщение отростков микроглии и гипертрофия астроцитов.
Интеграция этих трех подходов обеспечивает высокую достоверность данных: если повышение уровня мРНК сопровождается увеличением количества белка и соответствующим изменением клеточной морфологии в ткани, это дает веские доказательства наличия нейровоспалительного процесса.
**Заключение**
Описанные методы анализа экспрессии генов и белков представляют собой базовый инструментарий для любого исследования нейровоспаления. Точное следование протоколам, использование соответствующих контролей и комбинирование молекулярных и гистологических методов позволяют глубоко понять механизмы нейроиммунного взаимодействия и оценить эффективность потенциальных противовоспалительных препаратов. |
| BraTS 2021 | RSNA-ASNR-MICCAI BraTS 2021; CBICA/Пенсильванский университет | Сегментация/классификация опухолей головного мозга | МРТ | 2 040 случаев в контрольной когорте | 1251 аннотированного обучающего случая; по четыре модальности МРТ на случай | Зависимый от задачи | Отсутствие единого распределения классов по всему набору данных | Когорты, определенные по результатам испытания; разделение на экспериментальные группы не проводилось | Пример применимости протокола | Не используется для экспериментальной валидации | https://www.med.upenn.edu/cbica/brats2021/ |
Таблица 2: Характеристики наборов данных здравоохранения и применимость предлагаемого протокола. В таблице обобщены источники данных, варианты клинического применения, модальности, характеристики образцов, распределение классов, стратегии разделения наборов данных, статус валидации и информация об источниках для наборов данных здравоохранения, рассмотренных в протоколе. Набор данных Pima Indians Diabetes Dataset служит основным рабочим примером для валидации протокола.
| Элемент конфигурации | Условия фактической рабочей валидации | Статус / Интерпретация |
| Набор данных | Набор данных по диабету индейцев племени пима | Набор данных для фактической экспериментальной валидации |
| Алгоритм / Модель | Табличная прогностическая модель для бинарной классификации диабета | Используйте точное название архитектуры из протокола эксперимента, если оно задокументировано отдельно |
| Скорость обучения | 0.001 | Экспериментальная установка |
| Оптимизатор | Адам | Экспериментальная установка |
| Размер партии | 32 | Экспериментальная установка |
| Максимальное количество эпох | 100 | Экспериментальная установка |
| Функция потерь | Перекрестная энтропия | Экспериментальная установка |
| Функция активации | ReLU (линейный выпрямитель) | Экспериментальная установка |
| Дропаут (исключение нейронов) | 0.5 | Экспериментальная установка |
| Регуляризация весов (Weight Decay) | 1e-4 | Экспериментальная установка |
| Пакетная нормализация | Включено | Экспериментальная установка |
| Аугментация данных / Балансировка классов | Активировано | Указывайте SMOTE только в том случае, если данный метод действительно применялся в представленном прогоне |
| Стратегия валидации | 5-кратная перекрестная проверка | Экспериментальная установка |
| Ранняя остановка | Период ожидания (patience) = 10 эпох | Экспериментальная установка |
| Случайное число (Random Seed) | 42 | Используйте ту же конфигурацию начального значения (seed), которая была зафиксирована для данного эксперимента |
| Повторные запуски | 10 независимых запусков с использованием различных случайных значений (seeds) | Анализ воспроизводимости экспериментов |
| Размер входного изображения | Неприменимо | Набор данных Pima представляет собой табличные данные |
| Размерность признаков | 512 | Используйте только в том случае, если это окончательное реализованное представление признаков |
| Объяснимость | SHAP + LIME; контрфактуальное объяснение показано на рисунке 6 | Фактический отчет об анализе XAI (объяснимого искусственного интеллекта) |
| Метрики оценки | Точность, прецизионность, полнота, специфичность, F1-мера, коэффициент корреляции Мэтьюса (MCC), AUC-ROC, средняя точность (AP) | Приведенные метрики экспериментальной оценки |
| Оборудование | GPU NVIDIA | Замените на точную модель ГП, если она была зафиксирована |
| Программное обеспечение / Фреймворк | Python 3.1; Scikit-learn; компоненты фреймворка, использованные в реализации | Используйте точные версии пакетов, если они зафиксированы |
Таблица 3: Вычислительная среда, архитектура модели и конфигурация гиперпараметров, использованные для реализации протокола. В таблице указаны вычислительная платформа, программная среда, архитектура модели, конфигурация входных данных, параметры оптимизации, настройки регуляризации и параметры воспроизводимости, использованные для реализации предложенного рабочего процесса XAI.
| Параметр | Репрезентативная спецификация / результат |
| Оптимизатор | Адам |
| Скорость обучения | 0.001 |
| Размер партии | 32 |
| Максимальное количество эпох | 100 |
| Период ожидания при ранней остановке | 10 эпох |
| Лучшая эпоха | 78 |
| Эпоха ранней остановки | 88 |
| Наилучшее значение функции потерь на валидационной выборке | 0.142 |
| Наилучшая точность валидации | 94.6% |
| Результаты обучения | Потери при обучении и валидации снижались и достигли сходимости |
| Результаты валидации | Точность обучения и валидации увеличилась и стабилизировалась |
| Результат оптимизации | Наилучшая производительность модели достигнута на 78-й эпохе |
| Контроль переобучения | Ранняя остановка предотвратила дальнейшую оптимизацию после выбранной наилучшей эпохи |
Таблица 4: Спецификации обучения модели и результаты оптимизации. В таблице приведены сведения об оптимизаторе, скорости обучения, размере пакета, максимальном количестве эпох обучения, показателях валидации, сходимости обучения, результате оптимизации и стратегии контроля переобучения, использованных при разработке модели.
| Модель | Точность (%) | Прецизионность (%) | Полнота (%) | Специфичность (%) | F1-мера (%) | MCC | AUC (ROC) | AP (PR) |
| Предлагаемая модель XAI | 93.5 | 94.5 | 92.4 | 94.6 | 93.4 | 0.87 | 0.96 | 0.94 |
| Глубокое обучение (CNN) | 89.1 | 89.8 | 8.1 | 90 | 88.9 | 0.78 | 0.92 | 0.9 |
| Случайный лес (Random Forest) | 84.3 | 85 | 83.2 | 85.7 | 84.1 | 0.67 | 0.86 | 0.81 |
| Метод опорных векторов (SVM) | 78.6 | 79.3 | 7.1 | 80 | 78.2 | 0.54 | 0.79 | 0.72 |
| Базовая модель (мажоритарный класс) | 62.1 | 62.1 | 10 | 0 | 76.6 | 0 | 0.5 | 0.5 |
Таблица 5: Сравнение прогностической эффективности оцениваемых моделей. В таблице приведено сравнение предлагаемой модели XAI с оцениваемыми базовыми моделями по таким показателям, как точность (accuracy), прецизионность (precision), полнота (recall), специфичность (specificity), F1-мера (F1-score), коэффициент корреляции Мэтьюза, площадь под ROC-кривой и средняя точность (average precision).
| Валидационный анализ | Сравнение / Метрика | Статистический критерий | Статистический критерий | pp-значение | Экспериментальный результат / 95% ДИ |
| Пятикратная кросс-валидация | Точность | Описательная (среднее значение ± СО) | — | — | 93.01 ± 0.48% |
| Пятикратная перекрестная проверка | AUC-ROC | Описательная (среднее значение ± СО) | — | — | 0.954 ± 0.007 |
| Пятикратная перекрестная проверка | Точность | Описательная (средн ± СО) | — | — | 92.42 ± 0.87% |
| Пятикратная перекрестная проверка | Полнота | Описательная (среднее значение ± СО) | — | — | 93.02 ± 0.45% |
| Пятикратная кросс-валидация | F-мера | Описательная (среднее значение ± СО) | — | — | 92.72 ± 0.62% |
| 95% бутстрэп-доверительный интервал | Точность | Бутстрэп (10 повторных выборок) | — | — | 0.935 [0.897, 0.973] |
| 95%-й бутстрэп-доверительный интервал | Точность | Бутстрэп (1 00 повторных выборок) | — | — | 0.930 [0.890, 0.970] |
| 95% бутстрэп-доверительный интервал | Полнота (Recall) | Бутстреп (10 повторных выборок) | — | — | 0.922 [0.880, 0.963] |
| 95% бутстреп-доверительный интервал | F-мера | Бутстреп (1 0 повторных выборок) | — | — | 0.926 [0.887, 0.965] |
| 95%-ный бутстреп-доверительный интервал | AUC-ROC (площадь под ROC-кривой) | Бутстреп (10 повторных выборок) | — | — | 0.958 [0.931, 0.984] |
| Предлагаемая модель в сравнении с CNN | Точность (%) | Критерий Макнемара | 9.32 | 0.0023 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении с CNN | AUC-ROC | Критерий Делона | 3.87 | 0.0001 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении с CNN | F1-мера | Парный бутстрэп (1 0 ресемплов) | 2.81 | 0.0044 | Значительный (α = 0.05) |
| Предлагаемая модель в сравнении со случайным лесом | Точность (%) | Критерий Макнемара | 14.27 | 0.0002 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении со случайным лесом | AUC-ROC | Тест Делонга | 5.86 | <0.0001 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении со случайным лесом | F-мера | Парный бутстрап (1 0 повторных выборок) | 4.03 | 0.0003 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении с SVM | Точность (%) | Критерий Макнемара | 16.08 | <0.0001 | Значимый (α = 0.05) |
| Предлагаемая модель в сравнении с SVM | AUC-ROC | Тест Делонга | 6.95 | <0.0001 | Значительный (α = 0.05) |
| Предлагаемая модель в сравнении с SVM | F-мера | Парный бутстреп (1 0 повторных выборок) | 4.62 | <0.0001 | Значимый (α = 0.05) |
| Воспроизводимость при повторных запусках (10 независимых запусков) | AUC-ROC | Описательный (средн ± СО) | — | — | 0.957 ± 0.008 |
| Воспроизводимость при повторных запусках (10 независимых запусков) | Точность (%) | Описательная (среднее значение ± СКО) | — | — | 93.24 ± 0.74% |
| Воспроизводимость при повторных запусках (10 независимых запусков) | F1-мера (%) | Описательная (среднее значение ± СО) | — | — | 92.35 ± 0.92% |
Таблица 6: Результаты статистической валидации и воспроизводимости, полученные в ходе экспериментальной реализации с использованием набора данных Pima Indians Diabetes Dataset. В таблице обобщены показатели эффективности перекрестной проверки по пяти складываниям, 95% доверительные интервалы на основе бутстрэпа, статистические сравнения предлагаемой модели с базовыми моделями, а также воспроизводимость при повторных запусках с использованием 10 независимых случайных значений seed. В данной работе внешняя валидация и оценка клиническими экспертами не проводились.
| Нет. | Пункт контрольного списка | Необходимая документация | Рекомендуемая запись / проверка |
| 1 | Программная среда | Операционная система, язык программирования и программная среда | Зафиксируйте точные версии операционной системы и языка программирования. |
| 2 | Версии программного обеспечения и библиотек | Версии основных программных библиотек и пакетов | Записывайте точные названия и версии пакетов/библиотек. |
| 3 | Технические характеристики оборудования | Характеристики ЦП, ГП, ОЗУ, накопителей и ускорителей | Зафиксируйте используемые вычислительные мощности. |
| 4 | Идентификация набора данных | Название набора данных, источник, версия и информация о доступе | Записывайте точный источник/версию набора данных и дату доступа, где это применимо. |
| 5 | Характеристики набора данных | Размер выборки и распределение классов | Зафиксируйте общее количество образцов и распределение по классам для каждого разделения выборки. |
| 6 | Разбиение набора данных | Стратегия обучения, валидации и независимой тестовой выборки | Задокументируйте пропорции/количество разделения данных и стратификацию. |
| 7 | Предотвращение утечек данных | Процедура, используемая для предотвращения утечки информации | Разделение документов на обучающую и тестовую выборки и оценка параметров предобработки исключительно на обучающих данных. |
| 8 | Параметры предварительной обработки | Настройки очистки, обработки пропущенных значений, нормализации, кодирования и преобразования | Зафиксируйте все операции по предварительной обработке и значения параметров. |
| 9 | Аугментация данных | Методы аугментации и настройки параметров, если применимо | Документируйте методы, вероятности и диапазоны параметров. |
| 10 | Архитектура модели | Итоговая архитектура и конфигурация модели | Задокументируйте тип модели, слои/компоненты, размерности и функции активации. |
| 11 | Гиперпараметры | Обучение и оптимизация гиперпараметров | Зафиксируйте скорость обучения, размер пакета, оптимизатор, количество эпох, параметры ранней остановки и настроенные параметры. |
| 12 | Случайные значения (seeds) | Случайные числовые последовательности (seed), использованные для обеспечения воспроизводимости выполнения | Зафиксируйте значения seed (начальных значений генератора случайных чисел) для разделения данных, инициализации и повторных запусков. |
| 13 | Конфигурация обучения | Процедура обучения и стратегия выбора модели | Валидация документов, создание контрольных точек и критерии выбора модели. |
| 14 | Метрики оценки | Предопределенные метрики прогностической и статистической оценки | Зафиксируйте все зарегистрированные показатели эффективности. |
| 15 | Доверительные интервалы | Интервалы неопределенности для показателей эффективности | Задокументируйте процедуру оценки доверительного интервала (ДИ) и бутстреп-перевыборки, где это применимо. |
| 16 | Статистические критерии | Статистические процедуры для сравнения моделей | Документируйте тест, статистический показатель, значение p, порог значимости и допущения. |
| 17 | Анализ повторяемости результатов | Независимые запуски с использованием различных случайных чисел (сидов) | Зафиксируйте количество прогонов и среднее значение ± СО ключевых показателей. |
| 18 | Конфигурация интерпретируемости | Методы интерпретируемости и настройки параметров | Документируйте SHAP, LIME, Grad-CAM, внимание, контрфактуальные объяснения или иные применимые настройки. |
| 19 | Оценка объяснимости | Объективная оценка надежности и полезности объяснений | Достоверность документа, стабильность, недостоверность, локализация и экспертная оценка, где применимо. |
| 20 | Артефакты модели | Веса обученной модели и конфигурационные файлы | Архивируйте финальную обученную модель, архитектуру/конфигурацию и информацию по ее выбору. |
| 21 | Доступность кода | Код, необходимый для предобработки, обучения, оценки и генерации объяснений | Репозиторий записей или информация о контролируемом доступе к коду, если применимо. |
| 22 | Документация по исполнению | Команды, скрипты, конфигурационные файлы и инструкции | Зафиксируйте команды и конфигурацию, необходимые для воспроизведения рабочего процесса. |
| 23 | Выходная документация | Прогнозы, результаты оценки, рисунки и выходные данные объяснений | Архивируйте полученные результаты и свяжите их с соответствующим экспериментом/запуском. |
| 24 | Проверка воспроизводимости | Проверка согласованности результатов независимых повторных экспериментов | Сравните результаты повторных запусков и представьте предопределенные показатели вариабельности. |
Таблица 7: Контрольный список воспроизводимости для независимого повторения предлагаемого рабочего процесса XAI. В контрольном списке указаны требования к вычислительным ресурсам, наборам данных, предварительной обработке, разработке моделей, оценке, статистической проверке, интерпретируемости и документации, необходимые для воспроизведения экспериментального рабочего процесса.