Вестник Военного инновационного технополиса «ЭРА», 2023, T. 4, № 3, стр. 331-336
Особенности оценки данных фМРТ-исследований на групповом уровне
В. А. Орлов 1, *, Д. Г. Малахов 1, С. И. Карташов 1, М. В. Ковальчук 1, Ю. И. Холодный 1
1 Национальный исследовательский центр “Курчатовский институт”
Москва, Россия
* E-mail: Orlov_VA@nrcki.ru
Поступила в редакцию 05.09.2023
После доработки 05.09.2023
Принята к публикации 27.10.2023
Аннотация
фМРТ-исследования различной целевой направленности осуществляются на различных по количеству участников группах лиц, при этом традиционно используется групповая оценка получаемых фМРТ-данных, результаты которой зависят от ряда факторов. В рамках исследований “технологической” направленности, ориентированных на совершенствование системы оценки фМРТ-данных (в том числе на индивидуальном уровне), было признано целесообразным изучить различные варианты групповой оценки таких данных при различных критериях их значимости.
ВВЕДЕНИЕ
В последние годы Национальный исследовательский центр “Курчатовский институт” (НИЦ “КИ”) осуществляет нейрокогнитивные исследования различной целевой направленности, используя метод функциональной магнитно-резонансной томографии (фМРТ), дополненный регистрацией (с помощью МРТ-совместимого полиграфа (МРТсП)) динамики физиологических реакций исследуемого человека. Создание комплексного метода фМРТ-МРТсП-исследования [1, 2] позволило провести цикл экспериментов и заложить основы новой области отечественной нейронауки – нейрокриминалистики. Эти исследования дали возможность получить результаты, которые оказались полезными для фундаментальной науки и для выполнения научно-прикладных работ (например, в области медицины). Выполненные эксперименты наглядно показали необходимость совершенствования методических и аппаратно-программных средств фМРТ-МРТсП-исследований [3, 4], т.е. развития “технологического” направления исследований. При этом был сделан вывод об определенной ограниченности возможностей оценки фМРТ-данных на групповом уровне и указано на необходимость поиска путей и создания методов оценки фМРТ-данных на индивидуальном уровне [5–7], поскольку такая оценка имеет ключевое значение для развития и расширения использования метода фМРТ на практике, например в прикладной нейрокриминалистике и решении ряда задач медицинской практики.
Начатые в рамках исследований “технологического” направления работы по созданию одного из способов оценки на индивидуальном уровне данных, получаемых в ходе фМРТ, привели к созданию метода воксельной зоново-стимульной оценки (ВЗС-оценки) фМРТ-данных [8], который прошел испытание в экспериментальных условиях и показал свою перспективность.
Для дальнейшей разработки системы оценок фМРТ-данных на индивидуальном уровне было признано целесообразным внимательнее изучить итоги групповых оценок фМРТ-данных, полученных от одной и той же выборки исследованных лиц при различных критериях их значимости.
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
При изучении особенностей результатов фМРТ, получаемых на групповом уровне, анализу были подвергнуты фМРТ-данные, полученные в ходе фМРТ-МРТсП-исследования с использованием теста парадигмы скрываемой информации (активно применяемого в криминалистике) – теста со скрываемым именем (ТСИ), который моделирует сокрытие человеком личностно значимой для него информации, сохраняемой в его памяти на протяжении всей жизни. Не останавливаясь на технологии проведения ТСИ (подробное описание порядка его реализации представлено в [9]), отметим лишь, что в условиях фМРТ-МРТсП-исследования этот тест длился около 6–7 мин и содержал 30 стимулов, 25 из которых обеспечивали получение фМРТ-данных, подлежащих дальнейшей обработке.
В ходе экспериментов МРТ-данные регистрировались томографом 3 Tesla SIEMENS Magnetom Verio MR. Для получения структурного T1-взвешенного сагиттального 3D-изображения использовалась последовательность быстрого градиентного эха с параметрами: 176 срезов, TR = 1900 мс, TE = = 2.19 мс, толщина среза = 1 мм, угол поворота = 90°, время инверсии = 900 мс, FOV = 250 мм × × 218 мм2. Данные фМРТ были получены с параметрами: 51 срез, TR = 1110 мс, TE = 24 мс, толщина среза = 2 мм, угол поворота = 90°, FOV = = 192 × 192 мм2.
Эксперименты проводили с использованием сверхбыстрых последовательностей сканирования (TR = 1110 мс), поскольку такой режим способствовал повышению чувствительности методики исследования и обеспечивал приемлемый компромисс между качеством фМРТ-данных и размером области сканирования [5].
Данные функциональной и структурной МРТ были обработаны с использованием программного пакета SPM12 [10]. После конвертирования DICOM-файлов в NIFTI-формат все изображения центрировали вручную в передней комиссуре. Изображения EPI скорректировали на неоднородность магнитного поля с помощью FieldMap-тулбокса для SPM12, а также выполнили временную коррекцию сигналов фМРТ-данных, необходимость применения которой определяли выбором событийно связанной схемы эксперимента.
Анатомические и функциональные данные были нормированы в ICBM-стереотаксическую систему отсчета, а функциональные – сглажены с использованием гауссова фильтра с ядром 6 × 6 × × 6 мм3 FWHM.
Участниками фМРТ-МРТсП-экспериментов являлись студенты (36 мужчин и 36 женщин) в возрасте 21–23 лет, которые сообщили об отсутствии у них каких-либо заболеваний на момент участия в исследовании. Ранее уже было отмечено [9] влияние гендерного фактора на фМРТ-данные, регистрируемые в ходе фМРТ-МРТсП-исследования при реализации тестов в парадигме скрываемой информации. Этот факт определил необходимость формирования для данного исследования двух равных по количеству групп участников экспериментов мужского и женского пола.
Разрешение на проведение фМРТ-МРТсП-экспериментов было предоставлено этическим комитетом НИЦ “КИ”.
РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ
В связи с намеченным планом исследовать особенности результатов, получаемых после обработки фМРТ-данных на групповом уровне, было признано целесообразным рассмотреть результаты, полученные на единой выборке таких данных (72 участника экспериментов), а затем рассмотреть результаты каждой из групп (36 мужчин и 36 женщин).
В ходе этого исследования формировали групповые статистические карты с разными критериями их выделения, которыми были выбраны вероятности выделения: 1) p < 0.001, 5%; 2) p < 0.05 FWE, 1%; 3) p < 0.05 FWE, 2.5%; 4) p < 0.05 FWE, 5%.
Отметим, что выборки мужчин и женщин (каждая по 36 чел.) состояли из двух частей – “предварительной” (13 чел.) и “основной” (23 чел.), которые также были оценены порознь.
Некоторые из результатов оценки указанных выборок фМРТ-данных представлены в табл. 1, они позволили установить следующее.
Таблица 1
Зона | Количество вокселей | Аббревиатуры зон | Группы сравнения | Зо-на | Количество вокселей | Аббревиатуры зон | Группы сравнения | Зо-на | Количество вокселей | Аббревиатуры зон | Группы сравнения | |||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | 2 | 3 | 4 | 5 | 6 | 1 | 2 | 3 | 4 | 5 | 6 | 1 | 2 | 3 | 4 | 5 | 6 | |||||||||
1 | 8125 | FP r | * | # | * | 45 | 2052 | iLOC r | 89 | 96 | SCC l | |||||||||||||||
2 | 7012 | FP l | * | * | 46 | 2056 | iLOC l | 90 | 2541 | OP r | ||||||||||||||||
3 | 1363 | IC r | × | # | * | # | 47 | 752 | ICC r | # | 91 | 2722 | OP l | |||||||||||||
4 | 1319 | IC l | × | # | * | 48 | 655 | ICC l | 2.5 | 92 | 1295 | Thalam r | × | * | # | |||||||||||
5 | 2726 | SFG r | * | # | * | # | 49 | 962 | MedFC | 93 | 1342 | Thalam l | * | * | # | |||||||||||
6 | 2920 | SFG l | × | # | * | # | 50 | 775 | SMA r | * | # | * | # | 94 | 519 | Caudate r | × | * | ||||||||
7 | 2734 | MidFG r | * | # | * | 51 | 721 | SMA L | × | * | # | 95 | 536 | Caudate l | × | * | ||||||||||
8 | 2958 | MidFG l | × | # | * | 52 | 1127 | SubCalC | 96 | 810 | Putam. r | 1 | * | |||||||||||||
9 | 546 | IFG tri r | * | # | * | # | 53 | 1478 | PaCiG r | × | # | * | # | 97 | 869 | Putamen l | 1 | * | * | |||||||
10 | 654 | IFG tri l | * | # | 54 | 1427 | PaCiG l | × | # | * | # | 98 | 268 | Pallid. r | * | * | ||||||||||
11 | 695 | IFGoperr | × | # | * | # | 55 | 2618 | AC | × | # | * | # | 99 | 299 | Pallid. l | 1 | * | * | |||||||
12 | 765 | IFGoper l | × | # | * | # | 56 | 2408 | PC | × | # | * | 100 | 698 | Hippoс r | |||||||||||
13 | 4330 | PreCG r | * | * | 57 | 5571 | Precun. | 2.5 | * | 101 | 758 | Hippoc l | 1 | * | ||||||||||||
14 | 4441 | PreCG l | * | * | 58 | 690 | Cuneal r | 102 | 340 | Amygd r | 1 | * | # | |||||||||||||
15 | 2387 | TP r | * | 59 | 544 | Cuneal l | 103 | 324 | Amygd l | 2.5 | * | |||||||||||||||
16 | 2366 | TP l | * | # | 60 | 1442 | FOrb r | × | # | * | # | 104 | 84 | Accumb r | 2.5 | * | * | |||||||||
17 | 280 | aSTG r | 61 | 1677 | FOrb l | × | # | * | # | 105 | 111 | Accumb l | 1 | * | * | # | ||||||||||
18 | 272 | aSTG l | 62 | 660 | aPaHC r | 106 | 4097 | Brain-St | ||||||||||||||||||
19 | 399 | pSTG r | * | * | 63 | 594 | aPaHC l | 107 | 2292 | Cereb1 l | * | # | * | |||||||||||||
20 | 407 | pSTG l | × | # | * | 64 | 318 | pPaHC r | 108 | 2501 | Cereb1 r | * | # | * | ||||||||||||
21 | 413 | aMTG r | 65 | 382 | pPaHC l | 109 | 1917 | Cereb2 l | * | |||||||||||||||||
22 | 457 | aMTG l | 66 | 1792 | LG r | 110 | 2152 | Cereb2 r | ||||||||||||||||||
23 | 1373 | pMTG r | * | # | * | 67 | 1582 | LG l | 111 | 138 | Cereb3 l | |||||||||||||||
24 | 1363 | pMTG l | × | # | 68 | 294 | aTFusC r | 112 | 190 | Cereb3 r | ||||||||||||||||
25 | 1161 | toMTG r | × | # | * | # | 69 | 319 | aTFusC l | 113 | 815 | Cereb45 l | ||||||||||||||
26 | 865 | toMTG l | * | # | * | # | 70 | 712 | pTFusC r | 114 | 636 | Cereb45 r | ||||||||||||||
27 | 331 | aITG r | 71 | 876 | pTFusC l | 115 | 1248 | Cereb6 l | * | * | ||||||||||||||||
28 | 331 | aITG l | 72 | 817 | TOFusCr | 116 | 1525 | Cereb6 r | * | * | ||||||||||||||||
29 | 960 | pITG r | 73 | 640 | TOFusCl | 117 | 581 | Cereb7 l | * | |||||||||||||||||
30 | 1012 | pITG l | 74 | 903 | OFusG r | 118 | 543 | Cereb7 r | ||||||||||||||||||
31 | 775 | toITG r | 75 | 937 | OFusG l | 119 | 1888 | Cereb8 l | ||||||||||||||||||
32 | 714 | toITG l | 76 | 311 | FO r | × | # | * | # | 120 | 2297 | Cereb8 r | ||||||||||||||
33 | 3237 | PostCG r | 77 | 354 | FO l | × | # | * | 121 | 855 | Cereb9 l | |||||||||||||||
34 | 3692 | PostCG l | 78 | 890 | CO r | 2.5 | * | 122 | 811 | Cereb9 r | ||||||||||||||||
35 | 1467 | SPL r | * | # | * | 79 | 972 | CO l | 2.5 | 123 | 158 | Cereb10 l | ||||||||||||||
36 | 1485 | SPL l | * | * | 80 | 530 | PO r | 2.5 | * | # | 124 | 170 | Cereb10 r | |||||||||||||
37 | 794 | aSMG r | * | # | * | # | 81 | 564 | PO l | * | * | 125 | 47 | Ver12 | 1 | # | ||||||||||
38 | 967 | aSMG l | × | # | * | # | 82 | 376 | PP r | * | 126 | 223 | Ver3 | |||||||||||||
39 | 1239 | pSMGr | × | # | * | # | 83 | 356 | PP l | 2.5 | * | 127 | 635 | Ver45 | 1 | * | * | |||||||||
40 | 1094 | pSMG l | × | # | * | # | 84 | 286 | HG r | 128 | 335 | Ver6 | 1 | * | ||||||||||||
41 | 1451 | AG r | × | # | * | # | 85 | 307 | HG l | 129 | 188 | Ver7 | * | |||||||||||||
42 | 960 | AG l | * | # | * | # | 86 | 438 | PT r | 1 | * | 130 | 248 | Ver8 | ||||||||||||
43 | 4831 | sLOC r | 2.5 | 87 | 566 | PT l | * | * | 131 | 166 | Ver9 | * | ||||||||||||||
44 | 4962 | sLOC l | 1 | 88 | 164 | SCC r | 132 | 101 | Ver10 | |||||||||||||||||
Графа 1 | p < 0.05 FWE, 5% | 72 | м/ж | Графа 3 | p < 0.05 FWE, 5% | × | 36 | м | Графа 5 | p < 0.001, 5% | * | 23 | м | |||||||||||||
Графа 2 | p < 0.001, 5% | 72 | м/ж | Графа 3 | p < 0.001, 5% | * | 36 | м | Графа 6 | p < 0.001, 5% | # | 23 | ж | |||||||||||||
Графа 4 | p < 0.001, 5% | # | 36 | ж |
1. Применение, видимо, наиболее “жесткого” из существующих критериев – p < 0.05 FWE, 5% – к фМРТ-данным, полученным на выборке из 72 чел., позволило (графа 1) выделить 42 из 132 (т.е. 31.8%) зон мозга (от лобных долей до мозжечка), которые были вовлечены в процесс сокрытия значимой для человека информации в ходе ТСИ. Групповые статистические карты активности мозга показаны на рис. 1.
Рис. 1.
Групповые статистические карты (статистика Стьюдента) p < 0.05 с поправкой на множественные сравнения.

2. Применение менее “жесткого” и наиболее часто применяемого при оценке фМРТ-данных критерия – p < 0.001, 5% – к тем же данным той же выборки (72 чел.) позволило (графа 2) выделить 72 зоны мозга (т.е. 54.5% зон или на 71.4% больше, чем с помощью первого критерия). Групповые статистические карты активности мозга (при p < 0.001, 5%) представлены на рис. 2.
Рис. 2.
Групповые статистические карты (статистика Стьюдента) p < 0.001 без дополнительных поправок.

3. Помимо оценок, показанных в пп. 1 и 2, в исследовательских целях были применены “промежуточные” критерии (пп. 2 и 3), которые обычно не используют при оценке фМРТ-данных. Если взять за основу оценку фМРТ-данных при p < 0.001 5% (графа 2), то при оценке с критерием p < 0.05 FEW 1% не выделились 11 (отмечены знаком “1”) из 72 зон. При оценке с критерием p < 0.05 FWE 2.5% в дополнение к 11 не выделились еще 9 зон. Таким образом, динамика выделения зон мозга при ужесточении критериев оценки составила 72 (п. 2), 61, 52 и 42 зоны (п. 1).
4. Далее из общего массива фМРТ-данных (72 чел.) были выделены данные, полученные от мужчин (36 чел.) и оценены при двух критериях – p < 0.001 5% и p < 0.05 FWE 5%: итоги оценки представлены в графе 3. Как следует из табл. 1, оценка фМРТ-данных мужчин с критерием p < < 0.001 5% (обозначено звездочкой) в 90% случаев совпадает с оценкой фМРТ-данных с тем же критерием, полученных на выборке в 72 чел.: выделилось 65 зон.
Оценка фМРТ-данных мужчин с критерием p < 0.05 FWE 5% (обозначено крестиком) совпадает с оценкой фМРТ-данных с тем же критерием, полученных на выборке в 72 чел., лишь в 43% случаев: выделилось 18 зон.
5. Оценка фМРТ-данных женщин (36 чел.) с критерием p < 0.05 FWE, 5% не выявила выделившихся зон мозга, а с критерием p < 0.001, 5% (графа 4, обозначено решеткой) удалось выделить 38 зон, т.е. лишь 58% от числа активных зон у мужчин.
6. Оценка фМРТ-данных “основной” выборки мужчин (23 чел.) с критерием p < 0.001, 5% (графа 5, отмечено звездочкой) выделяет 58 зон мозга, что составляет 89% от объема выделившихся зон в выборке из 36 чел.
7. Оценка фМРТ-данных “основной” выборки женщин (23 чел.) с критерием p < 0.001, 5% (графа 6, отмечено решеткой) выделяет 27 зон мозга, что составляет 71% от объема выделившихся зон в женской выборке из 36 чел. и лишь 42% от объема выделившихся зон в мужской выборке из 36 чел.
ЗАКЛЮЧЕНИЕ
Представленные в табл. 1 результаты обработки фМРТ-данных при различных критериях их оценки позволяют сделать следующие выводы.
Таблица 1 наглядно без каких-либо дополнительных вычислений показывает, что активность зон мозга у женщин при выполнении ТСИ явно ниже, чем у мужчин; об этом же свидетельствует сказанное в пп. 4–6.
Ряд зон мозга человека (например, 5–6, 11–12, 25–26, 38–42 и др.) устойчиво проявляют активность в ходе выполнения ТСИ вне зависимости от величины и гендерного состава выборки, в то время как активность иных зон (например, 2, 10, 19, 24, 102 и др.) носит нестабильный, гендернозависимый характер.
Представленная обработка фМРТ-данных подтверждает установленное ранее при проведении нейрокогнитивных исследований требование [11] о строгом, по возможности, соблюдении социально-возрастной и гендерной однородности изучаемой выборки лиц.
Выполненное в итоге данного исследования определение устойчиво реагирующих (с учетом гендерного фактора) зон мозга является весьма полезным для дальнейшего совершенствования создаваемого в НИЦ “КИ” метода ВЗС-оценки фМРТ-данных.
Список литературы
Ковальчук М.В., Холодный Ю.И., Карташов С.И. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2020. Т. 1. № 1. С. 112.
Kovalchuk M.V., Kholodny Y.I. // Adv. Intel. Syst. Comp. 2020. P. 260.https://doi.org/10.1007/978-3-030-25719-4
Малахов Д.Г., Лисицин И.С., Скитева Л.И., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2021. Т. 2. № 2. С. 125.
Орлов В.А., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 1. С. 91.
Ковальчук М.В., Карташов С.И., Орлов В.А., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2021. Т. 2. № 4. С. 103.
Ковальчук М.В., Холодный Ю.И., Орлов В.А. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 2. С. 213.
Орлов В.А., Малахов Д.Г., Карташов С.И. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 2. С. 218.
Орлов В.А., Карташов С.И., Ковальчук М.В., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2023. (в печати).
Холодный Ю.И., Малахов Д.Г., Орлов В.А. и др. // Экспериментальная психология. 2021. Т. 14. № 3. С. 17.
Ashburner J. et al. SPM12 Manual. London: UCL Queen Square Institute of Neurology, 2020. URL: https://www.fil.ion.ucl.ac.uk/spm/doc/manual.pdf (Accessed 9.07.2021).
Bruni T. Cross-cultural variation and fMRI lie-detection // Technologies on the stand: Legal and ethical questions in neuroscience and robotics. 2012. P. 129.
Дополнительные материалы отсутствуют.
Инструменты
Вестник Военного инновационного технополиса «ЭРА»