ПРИМЕЧАНИЕ: Все программное обеспечение, используемое здесь, показано в таблице материалов. Данные, использованные в этом исследовании в демонстрационных целях, были взяты из проекта Human Connectome Project (HCP: http://www. humanconnectome.org)15. Все экспериментальные процедуры были одобрены Институциональным наблюдательным советом (IRB) при Вашингтонском университете. Данные визуализации в наборе данных HCP были получены с помощью модифицированного 3T сканера Siemens Skyra с 32-канальной головной катушкой. Другие параметры получения изображений подробно описаны в более ранней статье16. Для демонстрации были загружены минимальные предварительно обработанные данные, которые завершили следующие этапы предварительной обработки: коррекция градиентных искажений, коррекция движения, предварительная обработка карты поля, коррекция пространственных искажений, пространственная нормализация пространства Монреальского неврологического института (MNI), нормализация интенсивности и удаление поля смещения. Также могут быть использованы данные фМРТ в состоянии покоя, полученные в проектах исследователей.
1. Предварительная обработка данных
- Проверьте качество данных и исключите участников с отсутствующими данными повторного тестирования и чрезмерными движениями головы (3 мм при перемещении и 3° при повороте).
ПРИМЕЧАНИЕ: Пять участников были исключены, а 38 молодых людей (22-35 лет) были включены в основной анализ.
- Откройте набор инструментов теоретико-сетевого анализа графов (GRETNA)17 в MATLAB, чтобы выполнить дальнейшие шаги предварительной обработки. Щелкните на пакете FC Matrix Construction. Выберите путь к функциональному набору данных для загрузки документов NIFTI и выполните следующие шаги, как показано в опции конвейера на рисунке 1:
- Удалите первые 10 изображений, дважды щелкнув Номер временной точки для удаления в разделе Удалить первые изображения и введя 10.
- Пространственно сглаживание (полная ширина при половинном максимуме [FWHM] = [4 4 4] двойным щелчком FWHM (мм) в поле Пространственно сглаженное и вводом [4 4 4]).
- Регрессия ковариат. Выберите Сигналы белого вещества, Сигналы спинномозговой жидкости и Движение головы как TRUE. Выберите подходящую маску в соответствии с фактическим размером воксела, например, здесь маска с 2 мм, и выберите параметры Friston-24 для Head Motion.
- Временной фильтр. Введите значение TR в соответствии со временем повторения МРТ-сканирования (например, 720 мс) и удалите высокочастотный и низкочастотный шум, дважды щелкнув Полоса (Гц) и введя [0.01 0.1].
ПРИМЕЧАНИЕ: Результаты с регрессией и без регрессии сигналов всего мозга представлены ниже. При использовании непредварительно обработанных данных также рекомендуются хорошо зарекомендовавшие себя конвейеры, такие как fMRI-prep18 и Data Processing Assistant for Resting-State fMRI (DPARSF)19 .

Рисунок 1: Предварительный процесс Rs-fMRI и оценка функциональной связности сети. Изменены настройки предварительной обработки (удаление первых 10 изображений, пространственное сглаживание с FWHM 4 мм, линейное временное устранение тренда, регрессия сигналов белого вещества, сигналов спинномозговой жидкости (СМЖ) и движения головы с 24 параметрами, фильтрация полосы 0,01-0,1 Гц) и статической корреляции с Z-фактором Фишера. Сокращения: Rs-фМРТ = функциональная магнитно-резонансная томография в состоянии покоя; FWHM = полная ширина в половине максимума; Спинномозговая жидкость = спинномозговая жидкость. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этого рисунка.
2. Построение и анализ сети
ПРИМЕЧАНИЕ: Общий рабочий процесс построения и анализа навигационной сети сводится к трем основным этапам (рис. 2).

Рисунок 2: Общий рабочий процесс построения и анализа навигационной сети. (A) Выберите навигацию в качестве термина для поиска в базе данных Neurosynth. (B) Может быть сгенерирован список координат активации. (C) Запустите мета-анализ с использованием функций из Neurosynth, чтобы получить несколько карт мозга. (Д,Э) Используя метааналитическую карту и атлас парцелляции всего мозга (AICHA), можно создавать узлы (ROI). (F) Построение навигационной сети с использованием полученных навигационных узлов и их функциональной связности (оценка связности и сетевой анализ). Сокращения: ROI = область интереса; AICHA = атлас внутренней связности гомотопических областей. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этого рисунка.
- Определение узла сети
- Загрузите последнюю версию базы данных Neurosynth (neurosynth.org)20, набрав команду на Python:
>импортировать neurosynth как ns
>ns.dataset.download (path='./', unpack = True)
ПРИМЕЧАНИЕ: Архив набора данных ('current_data.tar.gz') содержит два файла: 'database.txt' и 'features.txt'. Они содержат все координаты активации из статей о нейровизуализации и теги мета-анализа, которые встречаются в этой статье с высокой частотой соответственно.
- Сгенерируйте новый экземпляр набора данных из database.txt и добавьте объекты к этим данным, введя команду:
> из neurosynth.base.dataset import Dataset
> dataset = Dataset('data/database.txt')
> dataset.add_features('данные/features.txt')
- Запустите мета-анализ с интересующим вас термином (т.е. 'navigation'), введя команду:
> ids = dataset.get_ids_by_features ('navigation', threshold=0.01)
> ma = мета. Метаанализ (набор данных, идентификаторы)
> ma.save_results('.', 'навигация')
ПРИМЕЧАНИЕ: Результатом мета-анализа является несколько карт мозга в формате NIFTI. Для контроля частоты ложных срабатываний применялся порог частоты ложных срабатываний (FDR), равный 0,01. На этом этапе необходимы соответствующие знания, чтобы гарантировать, что регионы, о которых часто сообщают, включены в метааналитическую карту. Аналогичные шаги могут быть применены для проведения мета-анализа других когнитивных функций, таких как язык и память.
- Определите интересующие вас кластеры, включив метааналитическую карту и атлас парцелляции всего мозга, введя команду из FSL:
>fslmaths navigation_0.01.nii.gz -bin navi_bin.nii.gz
>fslmaths navi_bin.nii.gz -mul AICHA/AAL.nii.gz navi_label_aicha/aal.nii.gz
>fslmaths navi_label_aicha/aal.nii.gz -thr n -uthr n label _n.nii.gz
>кластер -i метка _n.nii.gz -t 0.2 -o cluster_n.nii.gz
>fslmaths cluster_n.nii.gz -thr m -uthr m cluster_n_m.nii.gz
>fslmaths cluster_n_m.nii.gz -bin -mul x node_x.nii.gz
>fslmaths node_1.nii.gz -добавить ... -добавить node_x.nii.gz navi_AICHA/AAL_mask.nii.gz
ПРИМЕЧАНИЕ: Здесь были использованы два атласа: AAL и AICHA. AAL — это атлас, который использовался в первоначальном исследовании для определения узла12. Этот атлас был создан на основе анатомических профилей. Атлас внутренней связности гомотопических областей (AICHA)21 имеет гораздо более высокое разрешение и более высокую функциональную однородность. Мы определили интересующие нас регионы с помощью каждого из атласов.
- Наберите скрипты на Python для проверки размера каждого региона на карте:
>для i в np.arange(n)+1:
>____region_list.append(i)
>____size1_list.append(np.sum(img_dat==i))
>____size2_list.append(np.sum(aicha_img_dat==i))
>____pct_list.append(np.sum(img_dat==i)/np.sum(aicha_img_dat==i))
ПРИМЕЧАНИЕ: Целое число n в скрипте указывает общее количество регионов в рамках парцелляции AICHA и AAL (384 и 128 соответственно). Чтобы избежать влияния ложных кластеров, предполагается, что кластеры с относительно небольшими размерами (например, 100 вокселей) могут быть удалены. Используемый здесь атлас AICHA создан на основе данных функциональной связности, при этом каждый регион показывает однородность функциональной временной активности внутри себя.
- Оценка сетевого подключения
ПРИМЕЧАНИЕ: Набор инструментов GRETNA используется для оценки связности и анализа сети.- Щелкните на пакете FC Matrix Construction. Загрузите предварительно обработанные данные rs-fMRI, выбрав путь к функциональному набору данных. Выберите опцию статической корреляции . Загрузите узел, полученный на предыдущем шаге, в виде атласа для вычисления статической корреляции сигналов rs-fMRI каждой пары областей и переведите их в z-баллы Фишера для улучшения нормальности.
ПРИМЕЧАНИЕ: Подробная операция показана на рисунке 1. Матрицы навигационной сети N × N (N — количество узлов) для каждого участника будут получены в .txt формате.
- Получите положительную и взвешенную сеть, выполнив следующие шаги, как показано на рисунке 3.
- Щелкните пакет Сетевой анализ. Добавьте сетевые матрицы в окно Brain Connectivity Matrix и выберите один выходной каталог для подготовки.
- Для параметра конвейера в Network Configuration (Конфигурация сети) выберите positive (Положительный ) в матрице Sign of (Знак признака), который установит отрицательные соединения в матрице соединений функции равным 0 и устранит неоднозначные соединения22. Выберите тип сети как взвешенный , чтобы получить сеть с ненаправленным взвешенным коэффициентом.
ПРИМЕЧАНИЕ: Помимо взвешенных сетей, можно также бинаризовать сети для создания бинарных сетей для последующего анализа (с различными подходами), но взвешенная сеть часто считается более надежной23,24.
- Сетевой анализ
- Добавьте малый мир, глобальную эффективность, коэффициент кластеризации, длину кратчайшего пути, степень центральности и локальную эффективность в конвейер сетевого метрического анализа GRETNA, как показано на рисунке 3.
ПРИМЕЧАНИЕ: Малый мир и глобальная эффективность — это две глобальные сетевые метрики. В частности, сеть с малым миром может максимизировать эффективность передачи информации при сравнительно низких затратах на проводку. Глобальная эффективность отражает эффективность передачи параллельной информации в транспортной сети. Для метрик узловой сети степень центральности измеряет количество связей, подключенных к узлу. Кратчайшая длина пути, как и ее название, является основой для измерения интегрирования. Коэффициент кластеризации показывает степень взаимосвязи соседей узлов друг с другом. Локальная эффективность – это эффективность связи с узлом и его соседями (подробная формула и использование приведены в этих статьях) 17,25. Набор инструментов для подключения мозга (BCT)25 и другие наборы инструментов также могут быть использованы для вычисления сетевых метрик.
- Выберите Разреженность сети в методе определения пороговых значений, чтобы исключить искажающие эффекты паразитных соединений, и введите набор пороговых последовательностей (т. е. 0,05, 0,1, 0,15, 0,2, 0,25, 0,3, 0,35, 0,4, 0,45, 0,5) для дальнейшего определения соответствующего порога в соответствии со статистическими результатами.
ПРИМЕЧАНИЕ: Отношение ребер к максимальному числу ребер в сети с выборочным числом узлов известно как порог разреженности. Порог разреженности гарантирует, что разные особи имеют одинаковое количество ребер. Мы решили изучить различные пороговые значения для валидации, которые могли бы предоставить полезные данные для выбора оптимального порога в будущих исследованиях.
- Задайте случайное число сети равным 1 000 для создания случайных сетей с использованием алгоритма связиМаркова 26. Нажмите кнопку Выполнить , чтобы запустить конвейер в GRETNA после настройки всех шагов.
Примечание: Подобно настоящим мозговым сетям, случайные сети поддерживают одинаковое количество узлов, ребер и распределение степеней. Чтобы определить, являются ли они в значительной степени неслучайными топологически построенными, их сравним с сетями мозга. После запуска конвейера будет получена группа оценок для сетевых метрик для каждого из пороговых значений для дальнейшего статистического анализа.
- Определите оптимальное количество модулей в сети в четыре шага.
- Рассчитайте усредняющую навигационную сеть. Нажмите на кнопку Сравнение метрик и выберите Подключение. Загрузите сетевые матрицы, полученные выше, и выберите операцию Усредненное (Функциональное). Выберите направление вывода, чтобы сохранить усредненную сетевую матрицу; Подробнее см. рисунок 4 .
- Разделите среднюю сеть, полученную на предыдущем шаге, на 2, 3, 4 и 5 модулей, используя функцию spectralcluster в MATLAB.
- Вычислите долю узлов, разделенных на один модуль в REST 1 и REST 2 после выравнивания делений модуля с помощью скрипта procrustes_alignment.m. Используйте пропорцию в качестве показателя повторяемости раздела модуля.
- Выберите количество модулей с наибольшей повторяемостью.
- Статистический анализ
ПРИМЕЧАНИЕ: Следующие анализы предназначены в основном для валидации и не являются необходимыми при применении этого протокола к отдельным вариационным исследованиям.- Изучите сходство этих сетевых метрик между двумя сетями с различными типами стратегий определения узлов (т.е. новой, сгенерированной в настоящем исследовании, названной NaviNet_AICHA, и более ранней от Kong et al., названной NaviNet_AAL)12. Вычислите корреляцию Пирсона, используя функцию corrcoef в MATLAB, и повторите анализ для каждого порога разреженности.
ПРИМЕЧАНИЕ: После извлечения сетевых метрик можно проводить любой интересующий вас статистический анализ.
- Проверяем тестово-ретестовую надежность этих сетевых метрик с помощью функции ICC в MATLAB27,28, реализующей расчет внутриклассового коэффициента корреляции.
ПРИМЕЧАНИЕ: Исходные нескорректированные значения p были указаны в разделе репрезентативных результатов. 0,2 < ICC < 0,4 интерпретируется как показатель справедливой надежности при повторных испытаниях, а ICC > 0,4 интерпретируется как надежность от умеренной до хорошей29,30. Отрицательные баллы ICC были установлены на нуле, учитывая тот факт, что наличие отрицательных ICC бессмысленно и трудно интерпретируется31.

Рисунок 3: Анализ сетевых метрик. Этот анализ определяет взвешенные положительные сети с 10 пороговыми значениями. Вычислите две метрики глобальной сети (малые слова и эффективность), четыре узловые метрики сети (коэффициент кластеризации), кратчайшую длину пути, эффективность и степень центральности. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этого рисунка.

Рисунок 4: Расчет средних навигационных сетей. Усредненная (функциональная) операция помогает рассчитать средние сети всех участников. Пожалуйста, нажмите здесь, чтобы увидеть увеличенную версию этого рисунка.