「機能的磁気共鳴画像法」の版間の差分

ナビゲーションに移動 検索に移動
編集の要約なし
編集の要約なし
編集の要約なし
10行目: 10行目:
</div>
</div>
英:functional magnetic resonance imaging 英略称:fMRI 独:Funktionelle Magnetresonanztomographie 仏:Imagerie par résonance magnétique fonctionnelle
英:functional magnetic resonance imaging 英略称:fMRI 独:Funktionelle Magnetresonanztomographie 仏:Imagerie par résonance magnétique fonctionnelle
{{box|text= 抄録をお願いいたします。}}
{{box|text= 1段落程度の抄録をお願いいたします。}}


== はじめに ==
== はじめに ==
 機能的磁気共鳴画像(fMRI)とは、MRIを用いて生体の脳や脊髄を一定時間連続的に撮像し、脳活動(神経活動とシナプス活動等の総和)と相関するMRI信号の変動を非侵襲的に計測する技術である。1990年代の初頭に開発されるやいなや、当時ヒト脳機能イメージング研究手法の主流であったポジトロン断層像(PET)による血流・代謝測定を置き換えた。現在では、脳機能イメージング研究の代名詞として、健常脳の機能分離や機能連関の理解、あるいは精神・神経疾患の病態生理の解明のため、欠かすことのできないツールとなっている。ただしfMRIは、PETと同様、脳活動の本態である神経細胞の電気化学的活動そのものを測定しているのではなく、脳活動の代用マーカー(surrogate marker)としての局所酸素代謝・血流動態を画像化していることには留意が必要である。また、脳活動に由来するfMRI信号の変動は、脳活動以外の要因による信号変動と比較して必ずしも大きくないため、興味のある脳活動を抽出するために適切な画像・信号処理を行うことも重要である。本項目では、脳機能を解明するツールとしてのfMRIの原理、解析法とそれらを応用した脳科学研究の潮流を概説する。
 機能的磁気共鳴画像とは、磁気共鳴画像 (magnetic resonance imaging; MRI)を用いて生体の脳や脊髄を一定時間連続的に撮像し、脳活動(神経活動とシナプス活動等の総和)と相関するMRI信号の変動を非侵襲的に計測する技術である。


[[File:Hanakawa_fMRI_Fig1.png|thumb|right|'''図1. 神経血管連関の模式図'''<br>脳血流(CBF)は動脈〜小動脈(arteriole)から流入し、動脈血中では赤血球のヘモグロビンは酸素化(oxy-Hb)されている。刺激がない安静時'''(左)'''であっても、酸素は脳の基礎代謝要求により消費される。酸素は毛細血管で脳組織に供給され、酸素化ヘモグロビン(oxy-Hb)は、常磁性体である還元ヘモグロビン(deoxy-Hb)に変わる。外的刺激などによりシナプス入力と神経活動が増加すると'''(右)'''、局所の酸素・エネルギー代謝要求は安静時と比較して増加する。酸素消費により還元ヘモグロビン(常磁性体)が増加するはずだから、局所磁場が乱れてT2*[脳科学辞典wiki:T2*]が短縮するように思われる(陰性BOLD信号)。しかし、神経血管単位は基礎代謝要求の増加を検知して動脈血の流入を要求量以上に増加させるらしい。この過程にはプロスタグランジン(PG)や一酸化窒素(NO)が関わっているとされる。これらの影響の総和として、脳活動が増加する部分ではdeoxy-Hbが相対的に薄まって局所磁場が安定し、T2*延長が観察される。多くのfMRI法ではこのT2*の延長を陽性BOLD信号として計測している。]]
 1990年代の初頭に開発されるやいなや、当時ヒト脳機能イメージング研究手法の主流であったポジトロン断層像(PET)による血流・代謝測定を置き換えた。現在では、脳機能イメージング研究の代名詞として、健常脳の機能分離や機能連関の理解、あるいは精神・神経疾患の病態生理の解明のため、欠かすことのできないツールとなっている。ただしfMRIは、PETと同様、脳活動の本態である神経細胞の電気化学的活動そのものを測定しているのではなく、脳活動の代用マーカー(surrogate marker)としての局所酸素代謝・血流動態を画像化していることには留意が必要である。また、脳活動に由来するfMRI信号の変動は、脳活動以外の要因による信号変動と比較して必ずしも大きくないため、興味のある脳活動を抽出するために適切な画像・信号処理を行うことも重要である。本項目では、脳機能を解明するツールとしてのfMRIの原理、解析法とそれらを応用した脳科学研究の潮流を概説する。
 
[[File:Hanakawa_fMRI_Fig1.png|thumb|right|'''図1. 神経血管連関の模式図'''<br>脳血流(CBF)は動脈〜小動脈(arteriole)から流入し、動脈血中では赤血球のヘモグロビンは酸素化(oxy-Hb)されている。刺激がない安静時'''(左)'''であっても、酸素は脳の基礎代謝要求により消費される。酸素は毛細血管で脳組織に供給され、酸素化ヘモグロビン(oxy-Hb)は、常磁性体である還元ヘモグロビン(deoxy-Hb)に変わる。外的刺激などによりシナプス入力と神経活動が増加すると'''(右)'''、局所の酸素・エネルギー代謝要求は安静時と比較して増加する。酸素消費により還元ヘモグロビン(常磁性体)が増加するはずだから、局所磁場が乱れて[[磁気共鳴画像法#T2*|T2*]]が短縮するように思われる(陰性BOLD信号)。しかし、神経血管単位は基礎代謝要求の増加を検知して動脈血の流入を要求量以上に増加させるらしい。この過程にはプロスタグランジン(PG)や一酸化窒素(NO)が関わっているとされる。これらの影響の総和として、脳活動が増加する部分ではdeoxy-Hbが相対的に薄まって局所磁場が安定し、T2*延長が観察される。多くのfMRI法ではこのT2*の延長を陽性BOLD信号として計測している。]]
== 原理 ==
== 原理 ==
=== BOLD信号の発見 ===
=== BOLD信号の発見 ===
 1990年に小川誠二らはグラジエントエコー(GE)法MRIでマウス生体脳内の血管を撮像し、血管近傍のMRI信号が酸素飽和度に依存して変化することを発見した<ref name=Ogawa1990><pubmed>2161986</pubmed></ref>。小川はこの観察について、酸素飽和度が減少する部位では、赤血球内の還元ヘモグロビン(deoxy-Hb)が増加し、近傍の脳実質の信号が変化するためであると考えた。なぜなら還元ヘモグロビンは不対電子による常磁性(paramagnetic)を示し、還元ヘモグロビンの増加は局所磁場を乱し横磁化の減衰を増強するからである。小川はこの原理をblood oxygen-level-dependent (BOLD) コントラストと名付けた。さらに小川らは、ヒト健常ボランティアに視覚刺激を与えると一次視覚野のMRI信号が増加することを示した。このMRI信号の増加は、後述するように脳活動の増加に続く局所脳血流の過剰流入が酸素飽和度を上昇させ、還元ヘモグロビン濃度が薄まるために生じると説明されている <ref><pubmed>1631079</pubmed></ref>。なお、還元ヘモグロビンによる横磁化減衰の程度は、静磁場強度の2乗に比例して増加する<ref><pubmed>6275909</pubmed></ref>ため、高磁場MRIに測定上の優位性がある。
 1990年に小川誠二らはグラジエントエコー(GE)法MRIでマウス生体脳内の血管を撮像し、血管近傍のMRI信号が酸素飽和度に依存して変化することを発見した<ref name=Ogawa1990><pubmed>2161986</pubmed></ref>。小川はこの観察について、酸素飽和度が減少する部位では、赤血球内の還元ヘモグロビン(deoxy-Hb)が増加し、近傍の脳実質の信号が変化するためであると考えた。なぜなら還元ヘモグロビンは不対電子による常磁性(paramagnetic)を示し、還元ヘモグロビンの増加は局所磁場を乱し横磁化の減衰を増強するからである。小川はこの原理をblood oxygen-level-dependent (BOLD) コントラストと名付けた。
 
 さらに小川らは、ヒト健常ボランティアに視覚刺激を与えると一次視覚野のMRI信号が増加することを示した。このMRI信号の増加は、後述するように脳活動の増加に続く局所脳血流の過剰流入が酸素飽和度を上昇させ、還元ヘモグロビン濃度が薄まるために生じると説明されている <ref><pubmed>1631079</pubmed></ref>。なお、還元ヘモグロビンによる横磁化減衰の程度は、静磁場強度の2乗に比例して増加する<ref><pubmed>6275909</pubmed></ref>ため、高磁場MRIに測定上の優位性がある。


=== 神経血管連関 ===
=== 神経血管連関 ===
25行目: 29行目:
 2001年のLogothetisらによるサルを対象としたfMRIと神経活動の同時計測は、刺激や課題に伴う脳活動の増加とBOLD信号の関係性の理解に大きく貢献した<ref><pubmed>11449264</pubmed></ref>。ある視覚刺激条件において、マルチユニット神経活動(複数の神経細胞の出力)は一過性にのみ増加を示したが、BOLD信号と局所電場電位(local field potential, LFP)は一過性の増加に引き続く持続性の増加を示した。すなわちBOLD信号変化は、マルチユニット神経活動よりもLFPと良く相関する。LFPはシナプス活動、すなわち神経細胞への情報入力を反映するから、BOLD信号変化は皮質からの出力量よりも皮質への入力量との関係が強いことになる。最近では、超高磁場・高解像度のfMRIにより、皮質層ごとの情報処理を反映するBOLD信号変化の計測技術が精力的に開発されている<ref><pubmed> 26832438</pubmed></ref>。
 2001年のLogothetisらによるサルを対象としたfMRIと神経活動の同時計測は、刺激や課題に伴う脳活動の増加とBOLD信号の関係性の理解に大きく貢献した<ref><pubmed>11449264</pubmed></ref>。ある視覚刺激条件において、マルチユニット神経活動(複数の神経細胞の出力)は一過性にのみ増加を示したが、BOLD信号と局所電場電位(local field potential, LFP)は一過性の増加に引き続く持続性の増加を示した。すなわちBOLD信号変化は、マルチユニット神経活動よりもLFPと良く相関する。LFPはシナプス活動、すなわち神経細胞への情報入力を反映するから、BOLD信号変化は皮質からの出力量よりも皮質への入力量との関係が強いことになる。最近では、超高磁場・高解像度のfMRIにより、皮質層ごとの情報処理を反映するBOLD信号変化の計測技術が精力的に開発されている<ref><pubmed> 26832438</pubmed></ref>。


 脳活動とBOLD信号の相関関係は、外的刺激や課題の無い安静時にも観察される。そもそも課題遂行では脳の酸素代謝は数%しか増えず、脳のエネルギーは課題の無い(task free)安静時(resting state)の活動に大半が消費されている。これは脳が安静時にも組織的かつ活発な自発活動を示すことによる。1990年代後半のPET研究により、内側前頭前野、後部帯状回や両側外側頭頂葉などは、課題遂行時と比べて安静時にむしろ脳血流が増加することが知られていた<ref><pubmed> 25938726</pubmed></ref>。安静時に著明な自発的神経活動を示すこれらの領域は、Raichleによりデフォルトモードネットワーク(default mode network, DMN)と名づけられ、基底状態の脳の統合性に関わる内因性機構として提唱された。一方で、Biswalは1995年に安静状態のfMRIを解析し、両側運動感覚野の信号が主に0.1Hz以下の低い周波数(f)帯域において1/fのパターンで同期していることを見出していた<ref name=Biswal1995><pubmed> 8524021</pubmed></ref>。このようなfMRI信号同期は、安静状態神経ネットワーク(resting-state network, RSN)が有する機能結合(functional connectivity, FC)を反映すると考えられた。さらに興味深いことに、安静時fMRIにより、Raichleの提唱したDMN内の脳領域間には強い機能結合が存在することがわかった<ref><pubmed> 12506194</pubmed></ref>。このように2つの独立した研究の潮流が融合したことで、安静時fMRIを用いて局所の自発脳活動と領域間の機能結合状態を評価できる可能性に大きな注目が集まった。その後、安静時fMRIの信号同期性がサル脳における神経連絡性に対応していることも判明した<ref name=Ogawa1990></ref>。現在、安静時fMRIの同期現象は、神経連絡を持つ遠隔領域間で同期して発生する自発性のシナプス・神経活動(及びこれらに伴うBOLD効果)に基づいていると考えられている。領域間の同期の詳細を知るための解析手法の改善、覚醒時の基底状態としての意識との関連、精神疾患や認知症などの病態との関連、神経連絡性との対応などについて研究が進んでいる。
 脳活動とBOLD信号の相関関係は、外的刺激や課題の無い安静時にも観察される。そもそも課題遂行では脳の酸素代謝は数%しか増えず、脳のエネルギーは課題の無い(task free)安静時(resting state)の活動に大半が消費されている。これは脳が安静時にも組織的かつ活発な自発活動を示すことによる。1990年代後半のPET研究により、内側前頭前野、後部帯状回や両側外側頭頂葉などは、課題遂行時と比べて安静時にむしろ脳血流が増加することが知られていた<ref><pubmed> 25938726</pubmed></ref>。安静時に著明な自発的神経活動を示すこれらの領域は、Raichleによりデフォルトモードネットワーク(default mode network, DMN)と名づけられ、基底状態の脳の統合性に関わる内因性機構として提唱された。一方で、Biswalは1995年に安静状態のfMRIを解析し、両側運動感覚野の信号が主に0.1Hz以下の低い周波数(f)帯域において1/fのパターンで同期していることを見出していた<ref name=Biswal1995><pubmed> 8524021</pubmed></ref>。このようなfMRI信号同期は、安静状態神経ネットワーク(resting-state network, RSN)が有する機能結合(functional connectivity, FC)を反映すると考えられた。さらに興味深いことに、安静時fMRIにより、Raichleの提唱したDMN内の脳領域間には強い機能結合が存在することがわかった<ref><pubmed> 12506194</pubmed></ref>。このように2つの独立した研究の潮流が融合したことで、安静時fMRIを用いて局所の自発脳活動と領域間の機能結合状態を評価できる可能性に大きな注目が集まった。その後、安静時fMRIの信号同期性がサル脳における神経連絡性に対応していることも判明した<ref name=Ogawa1990></ref>
 
 現在、安静時fMRIの同期現象は、神経連絡を持つ遠隔領域間で同期して発生する自発性のシナプス・神経活動(及びこれらに伴うBOLD効果)に基づいていると考えられている。領域間の同期の詳細を知るための解析手法の改善、覚醒時の基底状態としての意識との関連、精神疾患や認知症などの病態との関連、神経連絡性との対応などについて研究が進んでいる。


[[File:Hanakawa_fMRI_Fig2.png|thumb|right|'''図2. 持続時間のごく短い感覚入力事象(t=0)に応答する血流動態応答関数(HRF)'''<br>MRI撮像の繰り返し時間(TR)が一秒の場合を示す。横軸の単位は秒、縦軸の単位は任意で、点線はfMRI信号強度の基準線(baseline)を示す。]]
[[File:Hanakawa_fMRI_Fig2.png|thumb|right|'''図2. 持続時間のごく短い感覚入力事象(t=0)に応答する血流動態応答関数(HRF)'''<br>MRI撮像の繰り返し時間(TR)が一秒の場合を示す。横軸の単位は秒、縦軸の単位は任意で、点線はfMRI信号強度の基準線(baseline)を示す。]]
33行目: 39行目:
== fMRIを用いた脳機能解明法 ==
== fMRIを用いた脳機能解明法 ==
=== fMRI撮像・実験デザイン ===
=== fMRI撮像・実験デザイン ===
 fMRIの撮影には、BOLD信号のコントラスト源であるT2*減衰[脳科学辞典wiki:T2*緩和]を鋭敏にとらえ、かつ時間分解能も高いエコープラナー画像(EPI)法[脳科学辞典wiki:エコープラナー法]を用いることがほとんどである。fMRIの性能を上げるためには、空間と時間両次元の分解能を高めることが理想的である。従って、撮像する際の時・空間解像度としては、MRI装置の性能を最大限に活かした高解像度が望ましいが、一方で解像度を上げすぎると信号雑音(S/N)比が低下することには注意が必要である。空間解像度と時間解像度もトレードオフの関係にあるが、マルチバンドEPI[脳科学辞典wiki:Multi-band/Simultaneous Multi-Slice撮像])などの高速撮像法の開発により長足の進歩があった。マルチバンドEPIは、従来のEPIに比べて空間解像度を1.2-1.5倍程度あるいは時間分解能を6-8倍近く向上させることができ、脳科学研究に汎用され始めている。
 fMRIの撮影には、BOLD信号のコントラスト源である[[T2*緩和|T2*減衰]]を鋭敏にとらえ、かつ時間分解能も高い[[エコープラナー法|エコープラナー画像(EPI)法]]を用いることがほとんどである。fMRIの性能を上げるためには、空間と時間両次元の分解能を高めることが理想的である。従って、撮像する際の時・空間解像度としては、MRI装置の性能を最大限に活かした高解像度が望ましいが、一方で解像度を上げすぎると信号雑音(S/N)比が低下することには注意が必要である。空間解像度と時間解像度もトレードオフの関係にあるが、[[磁気共鳴画像法#Multi-band/Simultaneous Multi-Slice撮像|マルチバンドEPI]]などの高速撮像法の開発により長足の進歩があった。マルチバンドEPIは、従来のEPIに比べて空間解像度を1.2-1.5倍程度あるいは時間分解能を6-8倍近く向上させることができ、脳科学研究に汎用され始めている。


 fMRIは撮像中の研究参加者(あるいは動物)の状態により2種類に大別される。安静状態の脳活動を計測する安静時(resting-stateまたはtask-free) fMRIと、研究参加者(あるいは動物)が指示された課題を遂行中に計測する課題(task) fMRIである。安静時fMRIでは、閉眼条件あるいは眼前の十字を固視した開眼条件で、「何も考えない」ように指示して5-10分程度のfMRI撮像を行うことが多い。課題の練習を行う必要がない楽な検査法であり、健常者のみならず患者や小児の脳研究にも使いやすく、上述のDMNの評価にも応用可能である。一方、課題fMRIでは、数秒から数十秒単位の課題遂行を間歇的に行うよう指示をして実行してもらう。通常、コンピュータを用いて感覚刺激を提示したり応答を記録したりする。そのため、課題fMRIを施行するためには、MRI撮像に影響を与えない感覚刺激・応答記録用のシステムを整備する必要がある。課題fMRIでは、課題遂行の時間情報と前述のHRFを畳み込んだBOLD信号変化モデルを作成し統計解析を行う(後述「fMRIの統計解析」参照)。課題遂行時間がやや長い場合(10-20秒以上)の場合は前述のDMNの周期や呼吸変動ノイズと周波数が近くなり、数秒以内の課題遂行の場合には体動や心拍の周波数帯と近くなる。そのためこうした関心の無い実験要因による解析精度の劣化の可能性を考慮した課題デザインを設計することが望ましいが、課題デザインについての詳細は本項目の範囲を超えるため別項目に譲りたい。
 fMRIは撮像中の研究参加者(あるいは動物)の状態により2種類に大別される。安静状態の脳活動を計測する安静時(resting-stateまたはtask-free)fMRIと、研究参加者(あるいは動物)が指示された課題を遂行中に計測する課題(task) fMRIである。安静時fMRIでは、閉眼条件あるいは眼前の十字を固視した開眼条件で、「何も考えない」ように指示して5-10分程度のfMRI撮像を行うことが多い。課題の練習を行う必要がない楽な検査法であり、健常者のみならず患者や小児の脳研究にも使いやすく、上述のDMNの評価にも応用可能である。一方、課題fMRIでは、数秒から数十秒単位の課題遂行を間歇的に行うよう指示をして実行してもらう。通常、コンピュータを用いて感覚刺激を提示したり応答を記録したりする。そのため、課題fMRIを施行するためには、MRI撮像に影響を与えない感覚刺激・応答記録用のシステムを整備する必要がある。課題fMRIでは、課題遂行の時間情報と前述のHRFを畳み込んだBOLD信号変化モデルを作成し統計解析を行う(後述「fMRIの統計解析」参照)。課題遂行時間がやや長い場合(10-20秒以上)の場合は前述のDMNの周期や呼吸変動ノイズと周波数が近くなり、数秒以内の課題遂行の場合には体動や心拍の周波数帯と近くなる。そのためこうした関心の無い実験要因による解析精度の劣化の可能性を考慮した課題デザインを設計することが望ましいが、課題デザインについての詳細は本項目の範囲を超えるため別項目に譲りたい。


=== fMRIデータの前処理 ===
=== fMRIデータの前処理 ===
55行目: 61行目:
==== 単変量解析 ====
==== 単変量解析 ====
univariate analysis
univariate analysis
 課題fMRIの解析には、課題に対応する脳活動とそれに応答するBOLD信号変化を合わせてモデル化し、この理論的な信号変化モデルと実際に測定したBOLD信号との適合度を体積画素/体素(voxel)ごとに計算し、課題に関連する脳活動を反映する統計画像を得る方法が広く用いられている。信号変化モデルは、複数の説明変数(課題遂行と相関する「脳活動」、頭部の動きなど脳活動以外の要因など)や脳活動に対するBOLD反応の時間的遅れやバラつき(血流動態関数、HRF)を考慮して構築する。多重回帰を含む一般線形モデル(general linear model, GLM)により、信号変化モデルとデータとの適合程度を調べ、各説明変数の重み付け係数値画像や誤差値画像を算出する(第1段階解析)('''図5A''')。臨床検査としての課題fMRI(術前の言語機能評価など)のように、個人の結果に興味がある場合は、係数値画像を誤差値画像で割って得られるt値画像などを用いた個人内要因の統計検定を行う。例えば、課題Aと課題B遂行時の脳活動の差を見るためには、係数の引き算に相当する重みづけ操作を行い(”コントラストを立てる”)、得られたt値画像を用いて統計検定を行う。医学研究や神経科学研究では、測定対象とした複数の個人を母集団からランダムに抽出したサンプルと考え、母集団についての推定を行いたいことが一般的である。その際は、研究者が関心のある実験要因に対する重みづけ係数画像について、個人間でのバラつきの違い「変量効果(random effect)」を考慮した統計検定(一般化線形混合モデル[generalized linear mixed model, GLMM])を行い、課題と相関する脳活動とその群間差を示す脳部位を検出する(2段階解析)('''図5B''')。
 課題fMRIの解析には、課題に対応する脳活動とそれに応答するBOLD信号変化を合わせてモデル化し、この理論的な信号変化モデルと実際に測定したBOLD信号との適合度を体積画素/体素(voxel)ごとに計算し、課題に関連する脳活動を反映する統計画像を得る方法が広く用いられている。信号変化モデルは、複数の説明変数(課題遂行と相関する「脳活動」、頭部の動きなど脳活動以外の要因など)や脳活動に対するBOLD反応の時間的遅れやバラつき(血流動態関数、HRF)を考慮して構築する。多重回帰を含む一般線形モデル(general linear model, GLM)により、信号変化モデルとデータとの適合程度を調べ、各説明変数の重み付け係数値画像や誤差値画像を算出する(第1段階解析)('''図5A''')
 
 臨床検査としての課題fMRI(術前の言語機能評価など)のように、個人の結果に興味がある場合は、係数値画像を誤差値画像で割って得られるt値画像などを用いた個人内要因の統計検定を行う。例えば、課題Aと課題B遂行時の脳活動の差を見るためには、係数の引き算に相当する重みづけ操作を行い(”コントラストを立てる”)、得られたt値画像を用いて統計検定を行う。
 
 医学研究や神経科学研究では、測定対象とした複数の個人を母集団からランダムに抽出したサンプルと考え、母集団についての推定を行いたいことが一般的である。その際は、研究者が関心のある実験要因に対する重みづけ係数画像について、個人間でのバラつきの違い「変量効果(random effect)」を考慮した統計検定(一般化線形混合モデル[generalized linear mixed model, GLMM])を行い、課題と相関する脳活動とその群間差を示す脳部位を検出する(2段階解析)('''図5B''')。
   
   
 ここで問題となるのが統計学的有意性の閾値の決定の方法である。画像統計では数十~百万単位の体素を同時に扱うため、適切な多重比較補正を行わないと、脳科学的に意味のない信号変化を脳活動として検出されてしまう危険(偽陽性)を伴う。そのため、古典的には脳の中に脳活動領域としての塊(cluster)を想定し、そのような塊が偶然形成される確率から閾値を設定する理論が採用されている<ref><pubmed> 9788073</pubmed></ref>。しかし、この統計閾値法の妥当性については議論が続いており<ref><pubmed> 27357684</pubmed></ref>、課題fMRIの再現性が低い主要因になっているとの見方がある<ref><pubmed> 32483374</pubmed></ref>。再現性は画像の平滑化処理や分画化などの前処置法にも依存するため、偽陽性率を適切に制御できる統計法の開発が行われている(例: 順列統計Permutation testなど<ref><pubmed> 24530839</pubmed></ref>。
 ここで問題となるのが統計学的有意性の閾値の決定の方法である。画像統計では数十~百万単位の体素を同時に扱うため、適切な多重比較補正を行わないと、脳科学的に意味のない信号変化を脳活動として検出されてしまう危険(偽陽性)を伴う。そのため、古典的には脳の中に脳活動領域としての塊(cluster)を想定し、そのような塊が偶然形成される確率から閾値を設定する理論が採用されている<ref><pubmed> 9788073</pubmed></ref>。しかし、この統計閾値法の妥当性については議論が続いており<ref><pubmed> 27357684</pubmed></ref>、課題fMRIの再現性が低い主要因になっているとの見方がある<ref><pubmed> 32483374</pubmed></ref>。再現性は画像の平滑化処理や分画化などの前処置法にも依存するため、偽陽性率を適切に制御できる統計法の開発が行われている(例: 順列統計Permutation testなど<ref><pubmed> 24530839</pubmed></ref>。
68行目: 78行目:


 システムとしての脳の仕組みを探るにあたり、機能局在を追求することだけでは限界がある。そこで脳部位間の繋がりが注目され、2005年には脳の網羅的な結線図を意味するコネクトーム(connectome)という概念が提唱された<ref><pubmed> 16201007</pubmed></ref>。脳コネクトームには様々な空間スケールがありうる。全脳を観察できるfMRIでは、脳領域間機能結合の正方行列を計算することで、全脳レベル(macroscale)の機能的コネクトーム解析が可能である。fMRIでは各脳領域間の連絡の強さの指標を機能結合(functional connetivity, FC)と呼ぶ。FCの産出法としては領域間の時系列信号のピアソン相関係数が良く使われる<ref name=Biswal1995></ref>。ある特定の脳領域を関心領域(シード)として設定し、相手方の関心領域(または体素)を、全脳に渡って隙間なく設定した相関解析を行えば、特定のシードに対する脳全体の機能結合マップが得られる。例えば、安静時fMRI解析において後部帯状回に関心領域を設定すれば、PETで見られたDMNがRSNとして抽出できる<ref><pubmed> 12506194</pubmed></ref><ref><pubmed> 17476267</pubmed></ref>。さらにシードを全脳の関心領域すべてに設定することで、多対多の網羅的な機能結合、すなわち機能的コネクトームが算出できる。このような機能的コネクトーム解析を行う際には、通常のピアソンの相関係数ではなく偏相関分析(partial correlation)を用いることで、2領域間に固有性の高いFCを評価できることが示唆されている<ref><pubmed> 22248579</pubmed></ref>。また機能的コネクトームを用いた新たな機能結合性モデルの提唱や、より高次な解析への展開(グラフ理論、独立成分分析、機械学習など)、他の指標や測定(脳波や臨床兆候など)との関連付け、疾患バイオマーカーの探索等の様々な研究も行われている。一方でヒトの脳コネクトーム研究単独では結果の妥当性の検証が難しいため、動物脳での検証も重要であり、例えば神経連絡トレーサーと機能的結合の比較は重要な課題である。
 システムとしての脳の仕組みを探るにあたり、機能局在を追求することだけでは限界がある。そこで脳部位間の繋がりが注目され、2005年には脳の網羅的な結線図を意味するコネクトーム(connectome)という概念が提唱された<ref><pubmed> 16201007</pubmed></ref>。脳コネクトームには様々な空間スケールがありうる。全脳を観察できるfMRIでは、脳領域間機能結合の正方行列を計算することで、全脳レベル(macroscale)の機能的コネクトーム解析が可能である。fMRIでは各脳領域間の連絡の強さの指標を機能結合(functional connetivity, FC)と呼ぶ。FCの産出法としては領域間の時系列信号のピアソン相関係数が良く使われる<ref name=Biswal1995></ref>。ある特定の脳領域を関心領域(シード)として設定し、相手方の関心領域(または体素)を、全脳に渡って隙間なく設定した相関解析を行えば、特定のシードに対する脳全体の機能結合マップが得られる。例えば、安静時fMRI解析において後部帯状回に関心領域を設定すれば、PETで見られたDMNがRSNとして抽出できる<ref><pubmed> 12506194</pubmed></ref><ref><pubmed> 17476267</pubmed></ref>。さらにシードを全脳の関心領域すべてに設定することで、多対多の網羅的な機能結合、すなわち機能的コネクトームが算出できる。このような機能的コネクトーム解析を行う際には、通常のピアソンの相関係数ではなく偏相関分析(partial correlation)を用いることで、2領域間に固有性の高いFCを評価できることが示唆されている<ref><pubmed> 22248579</pubmed></ref>。また機能的コネクトームを用いた新たな機能結合性モデルの提唱や、より高次な解析への展開(グラフ理論、独立成分分析、機械学習など)、他の指標や測定(脳波や臨床兆候など)との関連付け、疾患バイオマーカーの探索等の様々な研究も行われている。一方でヒトの脳コネクトーム研究単独では結果の妥当性の検証が難しいため、動物脳での検証も重要であり、例えば神経連絡トレーサーと機能的結合の比較は重要な課題である。
[[File:Hanakawa_fMRI_Fig6.png|thumb|'''図6. 安静時fMRI画像を前処置後、独立成分分析により得たDMN'''<br>後部帯状回を中心とし頭頂葉、前頭前野前内側部を含むデフォルトモードネットワーク(DMN)。左上は膨らました皮質表面にマッピングしたもの、右上は軸断像に表示、および同ネットワークの信号変化(左下)および周波数分析結果(右下)。周波数0.01-0.1Hzに高いパワーを持つネットワーク活動である。図3と同じfMRIデータで、動き補正・ノイズ処理を行ったあとに得られたもの。]]
[[File:Hanakawa_fMRI_Fig6.png|thumb|'''図6. 安静時fMRI画像を前処置後、独立成分分析により得たDMN'''<br>後部帯状回を中心とし頭頂葉、前頭前野前内側部を含むデフォルトモードネットワーク(DMN)。左上は膨らました皮質表面にマッピングしたもの、右上は軸断像に表示、および同ネットワークの信号変化(左下)および周波数分析結果(右下)。周波数0.01-0.1Hzに高いパワーを持つネットワーク活動である。図3と同じfMRIデータで、動き補正・ノイズ処理を行ったあとに得られたもの。]]


==== 独立成分分析 ====
==== 独立成分分析 ====
independent component analysis
independent component analysis
 ICAはデータ駆動型の解析でありブラインド信号分離の技術として多くの分野で活用されている。主成分分析(principle component analysis, PCA)と比べ、データマイニングの手法として柔軟であると考えられている。主成分分析では無相関、つまりベクトルとして直交する成分のセットが解析的に一意に定まる。これに対しICAでは、互いに他の情報を持たないという特性(独立性)が大きくなるよう成分を分離する。無相関であっても独立ではないことがあるから、独立性は相関性よりも柔軟な条件である。確率変数としてデータの空間的な分布を採るか、経時変化を採るかによって、空間的ICAと時間的ICAが可能である。fMRIデータが空間的に高次元であることから、初期のfMRIでの活用は空間的ICAにより数十個の独立成分を求めるものであった。神経解剖学的に解釈可能な皮質ネットワークを得ることができDMNの同定にも使用される('''図6''')。ICAは従前から脳波など電気生理データのノイズ除去に利用されてきたが、機械学習判別によって個人レベルのfMR元データから構造ノイズ成分を除去する前処理過程にもICAが利用されている(fMRIデータの前処理の項参照)。
 ICAはデータ駆動型の解析でありブラインド信号分離の技術として多くの分野で活用されている。主成分分析(principle component analysis, PCA)と比べ、データマイニングの手法として柔軟であると考えられている。主成分分析では無相関、つまりベクトルとして直交する成分のセットが解析的に一意に定まる。これに対しICAでは、互いに他の情報を持たないという特性(独立性)が大きくなるよう成分を分離する。無相関であっても独立ではないことがあるから、独立性は相関性よりも柔軟な条件である。確率変数としてデータの空間的な分布を採るか、経時変化を採るかによって、空間的ICAと時間的ICAが可能である。fMRIデータが空間的に高次元であることから、初期のfMRIでの活用は空間的ICAにより数十個の独立成分を求めるものであった。神経解剖学的に解釈可能な皮質ネットワークを得ることができDMNの同定にも使用される('''図6''')。ICAは従前から脳波など電気生理データのノイズ除去に利用されてきたが、機械学習判別によって個人レベルのfMR元データから構造ノイズ成分を除去する前処理過程にもICAが利用されている(fMRIデータの前処理の項参照)。
   
   
80行目: 92行目:


 1990年後半からfMRI研究は飛躍的に世界中に広まり多くの成果をあげてきた。その間に基礎となるMRI技術も大きく進展し、質の高いデータが得られるようになった。初期の研究の興味は脳の機能局在・分離の同定にあったが、2010年代以後はネットワークとしての脳の解明に興味がシフトした<ref><pubmed> 22481337</pubmed></ref>。一つの研究が扱うfMRIデータの数も、数10人から数百人、数万人の規模へと拡大し、種としてのヒトの脳機能構築やその複雑性の理解、個体差の理解へと興味が広がった。2010年から2016年まで、米国NIHの支援により大規模研究プロジェクトHuman Connectome Project(HCP)が行われた。このプロジェクトでは、ワシントン大学セントルイス医学校およびオックスフォード大学が中心となり約1200名の健康な若年成人被験者を対象として、安静時fMRI、標準的な課題を用いた課題fMRI、構造MRI(T1・T2強調画像)、拡散強調画像を撮像し、脳内の機能的結合と構造的結合について統合的に解析を進めた。このプロジェクトは(1)質が高い画像データを大量に取得する、(2)空間分解能を犠牲にすることなくデータを処理する手法を開発する、(3)FACT、すなわち皮質機能(Function)・連絡性(Connectivity)、構造(Architecture)、位置(Topography)の情報を統合した脳の領域分割(parcellation)を行う、(4)取得した生データ、解析データ、データ解析に必要なプログラム・コードを無料で公開する、といった多くの野心的な目標を達成し、現在も世界のMRI脳科学研究に大きな影響を与え続けている。特にFACT法は片側半球を180領域に分割することに成功し<ref name= Glasser2016b><pubmed> 27437579</pubmed></ref>、100年以上続くヒト脳の機能分画の歴史の中で初めて非侵襲手法により生きた個人の脳地図を作成することに成功した('''図7''')。こうした技術をさらに発展・拡張させ、脳の個体差の検出や、疾患の診断技術の開発が期待されており、更に規模の大きな研究プロジェクトが各国で推進されている(UKバイオバンク、ABCDなど)。また動物モデルへの技術展開も進めることで動物種を超えた脳のシステムの解明や、適切な動物疾患モデルの開発や評価法検証も期待される。
 1990年後半からfMRI研究は飛躍的に世界中に広まり多くの成果をあげてきた。その間に基礎となるMRI技術も大きく進展し、質の高いデータが得られるようになった。初期の研究の興味は脳の機能局在・分離の同定にあったが、2010年代以後はネットワークとしての脳の解明に興味がシフトした<ref><pubmed> 22481337</pubmed></ref>。一つの研究が扱うfMRIデータの数も、数10人から数百人、数万人の規模へと拡大し、種としてのヒトの脳機能構築やその複雑性の理解、個体差の理解へと興味が広がった。2010年から2016年まで、米国NIHの支援により大規模研究プロジェクトHuman Connectome Project(HCP)が行われた。このプロジェクトでは、ワシントン大学セントルイス医学校およびオックスフォード大学が中心となり約1200名の健康な若年成人被験者を対象として、安静時fMRI、標準的な課題を用いた課題fMRI、構造MRI(T1・T2強調画像)、拡散強調画像を撮像し、脳内の機能的結合と構造的結合について統合的に解析を進めた。このプロジェクトは(1)質が高い画像データを大量に取得する、(2)空間分解能を犠牲にすることなくデータを処理する手法を開発する、(3)FACT、すなわち皮質機能(Function)・連絡性(Connectivity)、構造(Architecture)、位置(Topography)の情報を統合した脳の領域分割(parcellation)を行う、(4)取得した生データ、解析データ、データ解析に必要なプログラム・コードを無料で公開する、といった多くの野心的な目標を達成し、現在も世界のMRI脳科学研究に大きな影響を与え続けている。特にFACT法は片側半球を180領域に分割することに成功し<ref name= Glasser2016b><pubmed> 27437579</pubmed></ref>、100年以上続くヒト脳の機能分画の歴史の中で初めて非侵襲手法により生きた個人の脳地図を作成することに成功した('''図7''')。こうした技術をさらに発展・拡張させ、脳の個体差の検出や、疾患の診断技術の開発が期待されており、更に規模の大きな研究プロジェクトが各国で推進されている(UKバイオバンク、ABCDなど)。また動物モデルへの技術展開も進めることで動物種を超えた脳のシステムの解明や、適切な動物疾患モデルの開発や評価法検証も期待される。
==関連項目==
* [[磁気共鳴画像法]]


==参考文献==
==参考文献==
<references />
<references />

案内メニュー