Вестник Военного инновационного технополиса «ЭРА», 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

Полный текст (PDF)

Аннотация

фМРТ-исследования различной целевой направленности осуществляются на различных по количеству участников группах лиц, при этом традиционно используется групповая оценка получаемых фМРТ-данных, результаты которой зависят от ряда факторов. В рамках исследований “технологической” направленности, ориентированных на совершенствование системы оценки фМРТ-данных (в том числе на индивидуальном уровне), было признано целесообразным изучить различные варианты групповой оценки таких данных при различных критериях их значимости.

ВВЕДЕНИЕ

В последние годы Национальный исследовательский центр “Курчатовский институт” (НИЦ “КИ”) осуществляет нейрокогнитивные исследования различной целевой направленности, используя метод функциональной магнитно-резонансной томографии (фМРТ), дополненный регистрацией (с помощью МРТ-совместимого полиграфа (МРТсП)) динамики физиологических реакций исследуемого человека. Создание комплексного метода фМРТ-МРТсП-исследования [1, 2] позволило провести цикл экспериментов и заложить основы новой области отечественной нейронауки – нейрокриминалистики. Эти исследования дали возможность получить результаты, которые оказались полезными для фундаментальной науки и для выполнения научно-прикладных работ (например, в области медицины). Выполненные эксперименты наглядно показали необходимость совершенствования методических и аппаратно-программных средств фМРТ-МРТсП-исследований [3, 4], т.е. развития “технологического” направления исследований. При этом был сделан вывод об определенной ограниченности возможностей оценки фМРТ-данных на групповом уровне и указано на необходимость поиска путей и создания методов оценки фМРТ-данных на индивидуальном уровне [57], поскольку такая оценка имеет ключевое значение для развития и расширения использования метода фМРТ на практике, например в прикладной нейрокриминалистике и решении ряда задач медицинской практики.

Начатые в рамках исследований “технологического” направления работы по созданию одного из способов оценки на индивидуальном уровне данных, получаемых в ходе фМРТ, привели к созданию метода воксельной зоново-стимульной оценки (ВЗС-оценки) фМРТ-данных [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] о строгом, по возможности, соблюдении социально-возрастной и гендерной однородности изучаемой выборки лиц.

Выполненное в итоге данного исследования определение устойчиво реагирующих (с учетом гендерного фактора) зон мозга является весьма полезным для дальнейшего совершенствования создаваемого в НИЦ “КИ” метода ВЗС-оценки фМРТ-данных.

Список литературы

  1. Ковальчук М.В., Холодный Ю.И., Карташов С.И. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2020. Т. 1. № 1. С. 112.

  2. Kovalchuk M.V., Kholodny Y.I. // Adv. Intel. Syst. Comp. 2020. P. 260.https://doi.org/10.1007/978-3-030-25719-4

  3. Малахов Д.Г., Лисицин И.С., Скитева Л.И., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2021. Т. 2. № 2. С. 125.

  4. Орлов В.А., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 1. С. 91.

  5. Ковальчук М.В., Карташов С.И., Орлов В.А., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2021. Т. 2. № 4. С. 103.

  6. Ковальчук М.В., Холодный Ю.И., Орлов В.А. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 2. С. 213.

  7. Орлов В.А., Малахов Д.Г., Карташов С.И. и др. // Вестник Военного инновационного технополиса “ЭРА”. 2022. Т. 3. № 2. С. 218.

  8. Орлов В.А., Карташов С.И., Ковальчук М.В., Холодный Ю.И. // Вестник Военного инновационного технополиса “ЭРА”. 2023. (в печати).

  9. Холодный Ю.И., Малахов Д.Г., Орлов В.А. и др. // Экспериментальная психология. 2021. Т. 14. № 3. С. 17.

  10. 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).

  11. Bruni T. Cross-cultural variation and fMRI lie-detection // Technologies on the stand: Legal and ethical questions in neuroscience and robotics. 2012. P. 129.

Дополнительные материалы отсутствуют.