ДолголетиеbioRxivMedina, M. C., Reddy, N. A., Bright, M. G., Sitek, K. R.3 мин чтенияpreprint
Оценка подходов к проверке выводов на основе данных функциональной МРТ одного субъекта с плотной выборкой всего мозга
Evaluating Approaches for Inference Testing of Whole-Brain Densely Sampled Single-Subject Task fMRI Data
Карточка статьи
Рубрика
Долголетие
Источник
bioRxiv
DOI
10.64898/2026.06.29.735344
Дата
30.06.2026
Автор
Medina, M. C., Reddy, N. A., Bright, M. G., Sitek, K. R.
Время чтения
3 мин
Это предварительная публикация, она не прошла научное рецензирование.
Аннотация
Картирование на основе заданий стало многообещающей техникой в функциональной магнитно-резонансной томографии (фМРТ) для надежной характеристики и картирования уникальных паттернов активности отдельных людей. Эти эксперименты включают в себя получение обширных данных визуализации у одного участника, что в конечном итоге улучшает чувствительность и специфичность индивидуальной функциональной локализации. Несмотря на свои преимущества, исследования в основном сосредотачивались на понимании специфической активации коры, что препятствует целостному взгляду на функциональный ответ на уровне систем, и до настоящего времени лучшие подходы к статистическому анализу контролируемых заданий на основе плотно собираемых данных всего мозга еще не были полностью установлены. Поэтому в этом исследовании мы собрали данные многоканальной плотно собранной фМРТ аудиологической системы, которая имеет основные подкорковые компоненты, и оценили чувствительность активации, а также стабильность активации в подмножествах данных, используя стандартные подходы к тестированию выводов на уровне всего мозга и специфичных регионов. Подходы на уровне всего мозга включали стандартные схемы вывода на уровне вокселей и кластеров с различными статистическими порогами, а также непараметрический метод перестановочного вывода. Регионально-специфические подходы включали исследовательские методы топ % t-статистик и непараметрические методы перестановочного вывода. Мы обнаружили, что подход на уровне вокселей с коррекцией по уровню ложных открытий (FDR) (p<0.05) продемонстрировал наивысшую чувствительность среди регионов и субъектов, а также наиболее последовательное обнаружение ожидаемых слуховых регионов, даже при более низкой продолжительности сканирования. Кроме того, мы обнаружили, что подход с использованием т-статистики на уровне регионов может быть полезным исследовательским инструментом функциональной локализации и дополнительным методом к стандартным подходам к тестированию выводов.
Краткое резюме
В исследовании анализируются подходы к проверке выводов на основе данных функциональной МРТ, полученных от одного субъекта. В частности, внимание уделяется чувствительности и стабильности активации при использовании методов, как на уровне всего мозга, так и на уровне регионов.
Практический вывод
Подход на уровне вокселей с коррекцией по уровню ложных открытий показывает наилучшие результаты в выявлении слуховых активностей, что может улучшить точность индивидуальной функциональной локализации.
Ограничения
Это предварительная публикация, она не прошла научное рецензирование. Исследование сосредоточено на одной системе — слуховой — и использует данные только одного участника, что может ограничить обобщаемость полученных результатов на более широкие выборки и разные системы.
Ветераны британской программы ядерных испытаний представляют собой возрастную группу бывших военнослужащих, которые могли подвергнуться воздействию ионизирующего излучения в ходе участия в ядерных испытаниях. Настоящее исследование было направлено на сравнение соматического мутационного ландшафта когорты из 30 ветеранов ядерных испытаний с соответствующей по возрасту когорте из 30 ветеранов-контролей. Варианты были определены на основе общедоступных данных полногеномного секвенирования с использованием биообработочной схемы, разработанной в соответствии с эталонными подходами, определенными Институтом Броуда. Полученный набор сырых SNV и INDEL вариантов для каждого отдельного ветерана подвергся нескольким этапам фильтрации, чтобы уменьшить шум, возникающий от общих мутаций, прежде чем среднее количество и типы были сопоставлены для каждой когорты с использованием теста Груббинса. Геномное распределение этих вариантов также было исследовано путем оценки наличия кластеризации мутаций, характерной для воздействия радиации (SNV и/или INDEL, происходящие в пределах 10 пар оснований) с использованием скользящего окна размером 10 пар оснований, а также отдельным идентификацией мутационных сигнатур путем подгонки SNV к базе данных COSMIC Human Cancer v3.4. При сравнении когорт ветеранов ядерных испытаний и контроля мы не обнаружили статистически значимого увеличения частоты какого-либо типа вариантов или кластеров. Dominant SBS сигнатуры в обеих когортах были характерны для процесс старения. Качественная оценка функционального воздействия наиболее часто наблюдаемых вариантов в каждой когорте показала, что они также ассоциированы с возрастом. Например, в контрольной когорте были обнаружены варианты в LINC02098-ETS1 и RCL1, генах, связанных с классическими возрастными состояниями, такими как потеря волос и остеоартрит. В когорте ветеранов ядерных испытаний мы наблюдали множественные варианты, влияющие на ген CHODL примерно у 40% участников. CHODL кодирует хондролектин, белок, важный для поддержания структурной целостности и функции тканей. В заключение, отсутствие значительных генетических различий между когортами, наряду с преобладанием мутаций, связанных с возрастом, согласуется с тем, что процесс старения является одним из основных факторов наблюдаемого соматического варианта у этих ветеранов, затмевая любые потенциальные воздействия окружающей среды, включая историческое воздействие радиации.
Болезнь Альцгеймера (БА) определяется и стадийно классифицируется на основе стереотипного, прогрессирующего накопления амилоид-бета (A{beta}) бляшек и гиперфосфорилированных тау (pTau) нейрофибриллярных клубков в различных областях мозга. Эти области существенно различаются по своей архитектуре и функциям, но в значительной степени сохраняют клеточный состав, с некоторой региональной спецификацией. Поскольку патология накапливается, определенные нейроны теряются, а не нейронные клетки переходят в состояния, ассоциированные с болезнью, но остается неясным, являются ли типы клеток, пораженные в одной области, такими же в других областях. Здесь мы расширили Атлас клеток мозга при болезни Альцгеймера в Сиэтле (SEA-AD) на десять неокортекальных и аллокортекальных областей, охватывающих кортикальную дугу канонического стадирования БА, профилируя приблизительно семь миллионов ядер от 84 доноров с помощью секвенирования РНК с одиночными ядрами, ATAC-seq и Multiome наряду с количественной нейропатологией и секвенированием всего генома. Ядра были сопоставлены с расширенной референсной таксономией Инициативы BRAIN, содержащей 207 типов клеток, иерархическая псевдо-прогрессия обеспечила непрерывные, донорные уровни оценки патологии БА в каждой области и по всему мозгу, совместно моделируя A{beta} и pTau. Изменения в клетках были как высокоспецифичными, так и поразительно последовательными: лишь около 30% типов клеток изменились в относительной частоте, но те, кто изменился, изменились в согласованном направлении в различных областях. Конкретные подмножества ингибирующих интернейронов Sst, Lamp5, Vip, Sncg и Pvalb, а также миелинизирующих олигодендроцитов были потеряны скорее всего на самых ранних стадиях у донора с минимальной патологией, наряду с первоначальным появлением ассоциированных с БА микроглией; потеря типов возбуждающих нейронов L2/3 и некоторых глубоких слоев, резкий рост микроглии и появление реактивных астроцитов последовали у доноров на более поздних стадиях. Регионально специализированные популяции также оказались уязвимыми, включая ожидаемые аллокортекальные типы, а также, неожиданно, специализированные возбуждающие нейроны слоя 4 (L4 IT) первичной зрительной коры (V1C) и смешанные ингибирующие интернейроны Sst и Pvalb. Ключевые изменения воспроизводились в трех независимых когортах, охватывающих более 700 дополнительных доноров. Мы исследовали два взаимодополняющих уязвимых населения в механистических подробностях: регионально специализированные нейроны V1C L4 IT, которые были потеряны поздно, несмотря на то, что их широко считали устойчивыми, и пан-кортикальные интернейроны Sst, которые потерялись в первые стадии болезни. Применяя многогранный ИИ-рабочий процесс, который конструировал гипотезы, основанные на литературе из дифференциального экспрессии для нейронов L4 IT, были выделены гипервозбудимость, частично опосредованная высоким уровнем экспрессии NMDA-рецепторов, как сходный фенотип уязвимости. Уязвимые интернейроны Sst сходились на гипервозбудимости через частично различные пути и обогащались экспрессией генов, приоритизированных в GWAS для БА, связывая их уязвимость с генетической архитектурой БА. Эти данные, доступные на SEA-AD.org, предоставляют многофункциональный региональный каркас для сообщества, чтобы исследовать молекулярные и клеточные изменения прогрессирования БА.
Обнаружение мутаций митохондриальной ДНК (мтДНК) на уровне отдельных клеток имеет значительный потенциал для определения клональных связей в масштабе, в сочетании с информацией о состоянии клеток у человека. Предыдущие методы фокусировались на мутациях с высокой гетероплазмией, которые, хотя и содержат ценные данные, редки и могут подвергаться функциональному отбору, что ограничивает информацию о родстве и потенциально вводит искажения в отслеживании. Хотя более сложно обнаружить, промежуточные и низкие мутации с гетероплазмией мтДНК имеют ценность из-за своего высокого разнообразия, обилию и меньшей предрасположенности к отбору. Для улучшения обнаружения мутаций мтДНК и содействия точному отслеживанию родословной мы разработали подход, называемый «Регуляторный мультиомикс с глубоким профилированием митохондриальных мутаций» (ReDeeM), представляющий собой интегрированную экспериментальную и вычислительную структуру. Здесь мы конкретно рассматриваем две аналитические проблемы, которые являются центральными для анализа родословной на основе мтДНК: достоверность рабочих процессов вызова вариантов и надежность филогенетического вывода. Мы демонстрируем, что, воспользовавшись ошибками поправки на основе консенсуса, вызовы мутаций мтДНК ReDeeM достигают высокой достоверности, соотносясь с настоящими мутационными сигнатурами даже для мутаций, поддерживаемых одной молекулой на клетку. Мы также разработали улучшенный подход к фильтрации после консенсуса, названный «filter2», который систематически идентифицирует и фильтрует остаточные артефакты, обогащенные краями, хотя они затрагивают лишь меньшинство вызовов мутаций. Чтобы систематически проверить ReDeeM, мы недавно провели эксперимент с лентивирусным баркодированием в человеческих гемопоэтических стволовых клетках, уникально пометив каждую клетку перед ее расширением и дифференциацией, чтобы предоставить основную правду для оценки точности отслеживания родословной. Такая валидация демонстрирует, что оригинальный аналитический подход ReDeeM (filter1) надежно восстанавливает истинную клональную структуру с высокой разрешающей способностью и отзывчивостью. Включение промежуточных и низких вариантов гетероплазмии (<10% на клетку) значительно улучшает выводы о родословной, в то время как исключение мутаций, поддерживаемых одной молекулой на клетку, удаляет истинный клональный сигнал и ухудшает восстановление истинных клонов. Оба фильтра filter1 и filter2 точно восстанавливают истинные клоны и превосходят ранее использовавшиеся подходы на основе мтДНК в отношении точности и отзывчивости, при этом filter2 обеспечивает дополнительные улучшения в производительности. Наконец, хотя ReDeeM продвигает отслеживание родословной на основе мтДНК, восстанавливая истинные клоны с высокой разрешающей способностью, количественная оценка филогенетической неопределенности от митохондриального наследования остается нерешенной. Для решения этой проблемы мы недавно разработали и проверили MitoDrift как платформу следующего поколения, учитывающую дрейф для отслеживания родословной мтДНК, которая дополнительно укрепляет филогенетический вывод и предоставляет интерпретируемые оценки неопределенности.
Интрогрессия от несampled гоминидных линий стала все более важным элементом истории человеческой эволюции. Недавние исследования показывают, что глубоко отличающиеся линии гоминид могли внести предковую составляющую либо непосредственно в современного человека, либо косвенно через интрогрессию от денисовцев, в то время как выводы остаются сложными из-за нехватки эталонных геномов, слабого сигнала и неопределенности в реконструкции глубоких генеалогий. Мы показываем аналитически и с помощью симуляций, что интрогрессия, опосредованная денисовцами, производит предсказуемые изменения в локальной коалиционной глубине, которые можно приблизительно оценить с помощью масштабируемых обобщенных статистик, в частности парной последовательной дивергенции, что указывает на то, что значительная информация о глубоко различающемся предке сохраняется в вариациях последовательностей без явной реконструкции генеалогий. Используя этот инсайт, мы разрабатываем DEEP (Точная оценка глубокого предка с использованием эффективных прокси), нейросетевую структуру без ARG для выявления кандидатных регионов суперархайческого предка. DEEP сохраняет обнаружимую мощность при низких уровнях ложных позитивов в широком диапазоне демографических параметров, остается масштабируемым и восстанавливает сигналы из небольших выборок. Применяя DEEP к океанийцам, тибетцам и ханьцзам, мы определяем примерно 0,4–0,6 % геномных окон с доказательствами суперархайческого предка. Кандидатные регионы показывают как значительное наложение, так и заметную неоднородность между популяциями, с повторным обогащением вблизи локуса HLA во всех популяциях, что указывает на то, что иммунные регионы повторно сохраняют глубоко различающуюся предковую составляющую.
По сравнению с традиционными грубыми объемными показателями, модели, основанные на поверхности, обеспечивают большую пространственную точность для понимания изменений в мозге, связанных с развитием, неврологическими и психиатрическими расстройствах. Глобальные инициативы по исследованию мозга объединяют данные со всего мира, чтобы обнаружить и улучшить маркеры мозга, связанные с болезнями. Здесь мы представляем инструментарий для анализа 3D-геометрии мозга, направленный на решение ключевых задач, стоящих перед крупномасштабными исследованиями нейровизуализации. Наша структура включает масштабируемые методы для интеграции многосайтовых данных, коррекции специфичных для сайта искажений, ускоренного статистического моделирования, интерпретируемого машинного обучения и интерактивной визуализации результатов. Инструментарий был протестирован на данных из 21 независимой выборки, участвующей в Рабочей группе по биполярному расстройству ENIGMA (N = 3,373). По сравнению с традиционными объемными характеристиками мы показываем, как меры формы подкорковых структур могут быть объединены по сайтам исследования для выявления пространственно сложных различий между диагностическими группами и их связью с общими методами лечения. Статистическое моделирование было ускорено с использованием алгоритма быстрого и эффективного смешанного эффекта (FEMA) и достигнуто сокращение времени вычисления в 16 раз по сравнению с традиционными подходами. Модели машинного обучения показали, что характеристики формы могут обеспечить большую предсказательную эффективность по сравнению с традиционными объемными данными как для задач диагностики, так и для прогнозирования лечения, при этом интерпретируемые карты весов предоставляют информацию о локальных признаках, определяющих эффективность модели.
Уточнение взаимосвязи между структурой и функцией важно для понимания работы мозга. В Drosophila melanogaster доступны данные о соединениях, полученные с помощью электронной микроскопии на основе FlyEM и FlyWire, а также данные о кальциевом изображении всего мозга. Мы применили методы предварительной обработки из функциональной магнитно-резонансной томографии (фМРТ) к данным о кальциевом изображении всего мозга и всесторонне исследовали оптимальные параметры. Затем мы обнаружили, что корреляция между функциональной и структурной связностью (FC-SC) линейно снижается с увеличением количества регионов интереса, и эта тенденция была одинаковой как у мушек, так и у людей. Мы также разработали новый, более надежный метод для количественной оценки степени пред- и постсинаптической сегрегации и исследовали это в мозге мух. Это показало, что многие нейроны имеют несегрегированные синапсы. Мы извлекли высоко несегрегированные синапсы и сравнили их с случайно извлеченными нулевыми матрицами SC. Их корреляция FC-SC была значительно выше, что указывает на то, что эти синапсы хорошо способствуют FC. Напротив, высоко сегрегированные синапсы показали значительно более низкую корреляцию FC-SC и меньше способствуют FC. Нейроны с несегрегированными синапсами, подобные неимпульсным нейронам, распределены по всему мозгу и считаются важными для FC.