Введение в количественную МРТ и микроструктуру мозга

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

Основные принципы физики МРТ

Все методы МРТ основаны на поведении ядер водорода (протонов) в сильном внешнем магнитном поле (]B0. Протоны обладают квантово-механическим свойством, называемым спином, которое дает им магнитный момент. В отсутствие внешнего поля эти магнитные моменты случайным образом ориентированы. При размещении в B0 вектор чистой намагниченности (NMV) развивается вдоль направления поля (]z. Протоны также прецируют вокруг B0 на частоте Лармора, которая пропорциональна силе поля: ω0 = γ B0, где γ — гиромагнитное соотношение (приблизительно 42,58 МГц/T для водорода).

Для генерации сигнала на частоту Лармора наносится радиочастотный (РЧ) импульс. Этот импульс вращает НМВ от продольной оси в поперечную плоскость. Угол вращения (угол скольжения) зависит от амплитуды и длительности РЧ импульса. После прекращения РЧ импульса НМВ подвергается двум независимым процессам релаксации:

  • T1 релаксация (спин-решетка релаксации): восстановление продольной намагниченности, как протоны высвобождают энергию в окружающую среду (решетка). T1 является постоянной времени для этого экспоненциального роста.
  • T2 релаксация (спин-спин релаксация): распад поперечной намагниченности вследствие разложения спинов, вызванного взаимодействиями между соседними протонами.Т2 — постоянная времени для этого экспоненциального распада.

Дополнительный распад, T2*, возникает из локальных неоднородностей магнитного поля и короче T2. Излучаемый радиочастотный сигнал (свободный индукционный распад или эхо) обнаруживается приемными катушками и пространственно кодируется с использованием градиентов магнитного поля. Измеренная интенсивность сигнала в вокселе зависит от плотности протона (PD), T1, T2 и T2*. Количественные методы МРТ изолируют эти параметры через преднамеренные вариации импульсной последовательности.

Ключевые количественные методы МРТ

Диффузионная МРТ (dMRI)

Диффузионная МРТ измеряет случайное (броуновское) движение молекул воды в ткани. Кажущийся коэффициент диффузии (ADC) отражает величину диффузии, в то время как диффузионная тензорная визуализация (DTI) моделирует направленность. В белом веществе вода диффундирует преимущественно вдоль аксонов (анизотропная диффузия). Применяя сильные градиенты магнитного поля в нескольких направлениях, dMRI захватывает распределение ориентации волокон. Расширенные модели, такие как диффузионная куртозная визуализация (DKI) и нейритовая ориентация дисперсии и плотности визуализации (NODDI), обеспечивают более конкретные микроструктурные метрики, включая аксональную плотность, ориентационную дисперсию и фракцию свободных и ограниченных водных бассейнов.

Количественное отображение T1 и T2

Количественное отображение времени релаксации T1 и T2 дает карты параметров (так называемые «карты»), которые не зависят от настроек сканера. T1 отображение часто использует восстановление инверсии или переменные последовательности угла поворота. T2 отображение обычно использует последовательности спин-эхо мультиэхо. Эти времена релаксации чувствительны к составу ткани: T1 короче в жирном или богатом миелинами белом веществе, чем в заполненной жидкостью спинномозговой жидкости (CSF); T2 длиннее в отеках и воспалении. Изменения в T1 и T2 коррелируют с демиелинацией, глиозом, осаждением железа и плотностью клеток.

Магнетизация передачи изображений (MTI)

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

Картирование количественной восприимчивости (QSM)

QSM измеряет магнитную восприимчивость ткани, которая зависит от присутствия парамагнитных веществ (например, железа, дезоксигемоглобина) и диамагнитных веществ (например, миелина, кальция). Обрабатывая фазовые изображения из последовательностей градиент-эхо, QSM генерирует карты локальной восприимчивости. Этот метод особенно ценен для количественной оценки железа мозга при нейродегенеративных заболеваниях (например, болезни Паркинсона) и для визуализации глубоких ядер мозга и вен. QSM дополняет методы, основанные на релаксации, поскольку он разделяет эффекты железа и миелина.

Физика, стоящая за диффузионной МРТ

Диффузионная МРТ опирается на явление самодиффузии воды в присутствии градиентов магнитного поля. Базовой импульсной последовательностью является импульсное градиентное спиновое эхо (PGSE), первоначально разработанное Стейскалем и Таннером. Два идентичных градиентных импульса (длительность δ, разделение Δ) размещены по обе стороны от перефокусирующего импульса 180°. Вращения, которые перемещаются между градиентными импульсами, накапливают чистый фазовый сдвиг, пропорциональный их смещению вдоль градиентного направления. Этот фазовый сдвиг приводит к затуханию сигнала, описанному уравнением Стейскала-Таннера:

S = S0 exp(−b · ADC)

где S0 является сигналом без градиентов, ADC является кажущимся коэффициентом диффузии, а b является коэффициентом диффузионного взвешивания, который зависит от силы градиента (G), δ и Δ: b = γ2 G2 δ2 (Δ− δ/3).

В анизотропных тканях, таких как белое вещество, одного АЦП недостаточно. ДТИ приобретает диффузионно-взвешенные изображения по меньшей мере в шести неколлинеарных направлениях для реконструкции диффузионного тензора 3×3. Собственные значения тензора (λ1, λ2, λ3) и собственные векторы описывают величину и ориентацию диффузии. Производные метрики включают:

  • Фракционная анизотропия (FA): мера направленной когерентности (0 = изотропная, 1 = идеально анизотропная)
  • Средняя диффузионность (MD): среднее значение собственных значений, представляющее общую величину диффузии
  • Осевая диффузивность (AD): λ1, чувствительная к аксональной целостности
  • Радиальная диффузивность (RD): среднее значение λ2 и λ3, чувствительное к целостности миелина

Однако ДТИ предполагает гауссову диффузию, которая является упрощением в сложных тканевых микросредах. Более продвинутые модели учитывают негауссово поведение. Например, диффузионная куртозная визуализация (ДКИ) оценивает избыточный куртоз (отклонение от гауссовской), обеспечивая дополнительную чувствительность к тканевой гетерогенности. NODDI использует многокамерную модель (внутриклеточную, внеклеточную и CSF) для оценки плотности нейрита, дисперсии ориентации и обмена водой. Эти методы требуют приобретения при множественных b-значениях и направлениях, увеличивая время сканирования, но давая более богатую микроструктурную информацию.

Выбор b-значения имеет решающее значение. Для визуализации мозга типичные b-значения варьируются от 1000 с/мм2 (DTI) до 2000-3000 с/мм2 (DKI) и даже выше для NODDI. Высокие b-значения улучшают чувствительность к ограниченной диффузии, но уменьшают отношение сигнал-шум (SNR). Градиентное оборудование - максимальная градиентная амплитуда и скорость убоя - ограничивает достижимые b-значения и способность разрешать очень короткие времена диффузии.

Время релаксации и микроструктура тканей

Релаксация T1 и T2 регулируется локальной молекулярной средой. Для релаксации T1 (спин-решетка) эффективность передачи энергии зависит от скорости падения молекул воды относительно частоты Лармора. Вода, связанная с макромолекулами или ограниченная в узких пространствах (например, между слоями миелина), показывает более быстрое релаксацию T1 (более короткий T1) из-за более медленного падения и усиленных диполярных взаимодействий. И наоборот, свободная вода (CSF) быстро падает, делая T1 неэффективным и, следовательно, длинным (≈4 секунды при 3T). T1 также подвержен влиянию парамагнитных центров, таких как железо, которые ускоряют релаксацию.

Расслабление T2 (спин-спин) является результатом снижения поперечной намагниченности из-за местных изменений поля. Неоднородности мелкомасштабного поля от молекулярного движения вызывают необратимый распад T2. В ткани компартментализация воды - например, внутри-аксоновая, вне-аксоновая и миелиновая вода - производит многоэкспоненциальный распад T2. Миелиновая вода имеет очень короткий T2 (10-30 мс), в то время как внутри- и вне-аксоновая вода имеют промежуточные значения T2 (70-100 мс). Многоэхо-Т2 визуализация может разрешать эти компоненты, обеспечивая карты фракции миелиновой воды (MWF) - прямая мера содержания миелина.

T1 и T2 также зависят от силы поля. При более высоких полях (7T против 3T) T1 увеличивается, а T2 немного уменьшается, влияя на контраст. Коррекция неоднородностей поля необходима при количественном отображении. Такие методы, как DESPOT1 (управляемое равновесным одноимпульсным наблюдением T1) или переменным методом угла поворота, позволяют быстрое отображение T1. Для T2 распространены последовательности многоэхо-спина с стимулированной компенсацией эха или многоэхо-градиент-эхо (для T2*).

Количественное отображение времени релаксации нашло широкое применение при изучении изменений белого вещества. Например, при рассеянном склерозе удлинение T1 и T2 наблюдается при демиелинизирующих поражениях, тогда как укорочение T1 может происходить в богатых железом структурах, таких как базальные ганглии. При болезни Альцгеймера повышенный T2 в гиппокампальных областях коррелирует с тау-патологией. Сочетание показателей T1, T2 и диффузии обеспечивает мультипараметрический отпечаток микроструктуры ткани.

Передача магнитизации и химический обмен

Изображение передачи магнитизации (МТ) является мощным методом исследования макромолекулярного содержимого. Эффект МТ измеряется путем сравнения изображений, полученных с импульсом насыщения вне резонанса и без него. Соотношение передачи намагниченности (МТР) рассчитывается как:

MTR = (S0 − Ssat) / S0

где S0 — сигнал без насыщения, а Ssat — сигнал после насыщения. Более низкий MTR указывает на снижение макромолекулярной целостности, часто из-за демиелинизации. В мозге белое вещество имеет более высокий MTR (≈40—50%), чем серое вещество (≈30—35%), отражая его обильное миелиновое содержание. МТ-изображение не обеспечивает абсолютных значений содержания миелина, поскольку MTR также зависит от содержания воды и обменных курсов. Однако оно остается надежным суррогатным маркером.

Передача насыщения химического обмена (CEST) является расширением MT, которое нацелено на конкретные обмениваемые протоны (например, амидные, аминные, гидроксильные группы) на метаболиты. Например, визуализация переноса амидных протонов (APT) обнаруживает мобильные белки и пептиды в опухолях, помогая в дифференцировке глиом высокого класса от более низкого уровня. Методы CEST требуют тщательного моделирования рН и температурной зависимости, но они предлагают окно в метаболические процессы при разрешении МРТ.

Применение в неврологических заболеваниях

Количественные методы МРТ преобразовали оценку неврологических расстройств. В рассеянном склерозе (MS), DTI и МТ-визуализация обнаруживают микроструктурные повреждения в нормальном появлении белого вещества (NAWM) задолго до появления обычных поражений. Сокращение MTR в NAWM предсказывает прогрессирование инвалидности. Картирование T1 и T2 может дифференцировать острые воспалительные поражения от хронических. QSM обнаруживает накопление железа в глубоком сером веществе, которое коррелирует с нейродегенерацией.

При болезни Альцгеймера (AD) диффузионная МРТ показывает повышенную среднюю диффузивность и снижение фракционной анизотропии в гиппокампе и туловище, отражая ранние дегенеративные изменения. Количественное картирование T1 и T2 было связано с бета-амилоидным и тау-отложением на животных моделях. МТ-визуализация демонстрирует снижение МТР в задней поясной и прекунеусной областях, пораженных в начале AD. Методы qMRI на основе перфузии, такие как маркировка артериального спина (ASL), также показывают снижение мозгового кровотока в AD.

В опухолях головного мозга диффузионно-взвешенная визуализация (DWI) с картированием ADC помогает отличить высокосортную от низкосортной глиомы — более низкий ADC в высокосортных опухолях из-за повышенной клеточности. Перфузионная МРТ (динамический контраст восприимчивости, DSC) дает относительные карты объема мозговой крови (rCBV), которые направляют биопсию и оценивают ответ на лечение. Количественное картирование T1 с контрастными агентами на основе гадолиния может измерять проницаемость гематоэнцефалического барьера более точно, чем полуколичественное усиление.

Stroke is another area where qMRI plays a critical role. Perfusion-diffusion mismatch (using DWI and PWI) identifies salvageable penumbra. Quantitative maps of T₂ and T₂* help distinguish hemorrhagic from ischemic stroke and evaluate edema progression. Advanced diffusion models like DKI have been shown to detect microstructural damage in chronic stroke beyond the infarct core.

Наконец, невродегенеративные заболевания, такие как болезнь Паркинсона, болезнь Хантингтона и боковой амиотрофический склероз (ALS) выигрывают от qMRI. QSM показывает накопление железа в черной субстанции и путамене при болезни Паркинсона, помогая диагностировать. Показатели ДТИ в кортикоспинальных трактах коррелируют с двигательной инвалидностью в БАС. По мере того, как эти методы становятся более стандартизированными, они включаются в клинические протоколы и испытания лекарств.

Будущие направления и вызовы

Поле количественной МРТ головного мозга продолжает быстро развиваться. Одной из основных тенденций является толчок к сверхвысокопольной МРТ (7T и выше). Более высокие силы поля улучшают SNR и пространственное разрешение, позволяя визуализировать корковые слои, подкорковые ядра и небольшие пучки белого вещества. Однако они также создают проблемы: увеличение неоднородностей B0 и B1, более сильные артефакты восприимчивости и более крупное осаждение мощности RF. Разрабатываются расширенные последовательности импульсов (например, параллельная передача, адиабатические импульсы) и методы переключения для смягчения этих проблем.

Машинное обучение (ML) всё чаще применяется к количественной МРТ. Модели глубокого обучения могут ускорить отображение параметров из непроверенных данных, обесценить приобретения с низким SNR и даже предсказать микроструктурные свойства непосредственно из необработанных изображений. Сегментация на основе ML карт qMRI позволяет автоматически определять структуры мозга и количественную оценку нагрузки поражения. Также растет интерес к синтетической МР-визуализации , где одно многоэхо-приобретение используется для создания карт T1, T2 и PD, а затем синтезирует любой желаемый контраст — экономя время сканирования при сохранении количественной точности.

Стандартизация остаётся препятствием.Количественные значения МРТ различаются в разных сканерах и участках из-за различий в аппаратном обеспечении, импульсных последовательностях и реконструкции. Такие инициативы, как Альянс биомаркеров количественной визуализации (QIBA) и Международное общество магнитного резонанса в медицине (ISMRM) работают над протоколами гармонизации, фантомными калибровками и консенсусными последовательностями. Для клинического принятия эти методы должны обеспечивать воспроизводимые и диагностически функциональные метрики.

Другим перспективным направлением является мультипараметрическая визуализация, где несколько параметров qMRI объединяются в один классификатор или биомаркер. Например, данные T1, T2, MT и диффузии могут быть интегрированы для генерации «подписи» нормального старения в белом веществе. При болезни Альцгеймера композит объемных, диффузионных и MT метрик превосходит любой один параметр в прогнозировании когнитивного снижения. По мере увеличения вычислительной мощности становится осуществимым многомерный анализ по вокселю.

Наконец, исследуются новые контрастные механизмы. МРТ натрия и МРТ фосфора обеспечивают прямую метаболическую информацию, но страдают от низкого СНР. Химические агенты обмена и гиперполяризованные газы (например, 129Xe) открывают новые окна в рН ткани, метаболизм и микроструктуру. Хотя они остаются в основном исследовательскими инструментами, они указывают на еще более богатое будущее для количественной визуализации микроструктуры мозга.

Заключение

Количественные методы МРТ построены на строгих физических принципах - расслабление, диффузия, передача намагниченности и магнитная восприимчивость. Преобразуя эти принципы в точные измерения, qMRI позволяет исследователям и клиницистам исследовать микроструктурную среду живого мозга без инвазивных процедур. Сочетание диффузионной визуализации, картирования релаксации, МТИ и QSM предлагает всеобъемлющий взгляд на здоровье тканей при заболеваниях, начиная от рассеянного склероза и заканчивая нейродегенерацией. Текущие достижения в области аппаратного обеспечения, дизайна импульсной последовательности и анализа данных обещают сделать количественную МРТ еще более точной, быстрой и более широко доступной, в конечном итоге улучшая диагностику, мониторинг лечения и наше фундаментальное понимание микроструктуры мозга.