:脳科学辞典(BOLD信号から転送)これは、このページの承認済み版であり、最新版でもあります。MRI)を用いて1990年代から脳機能を解明するための研究する技術として使用され、ヒト脳の機能解明に大きく貢献した。開発が急速に進んでいる。臨床診断法として確立した技術ではなく、あくまで研究用として使用されている。MRI)を用いて生体の脳や脊髄を一定時間連続的に撮像し、脳活動(神経活動とシナプス活動等の総和)と相関するMRI信号の変動を非侵襲的に計測する技術である。1990年代の初頭に開発されるやいなや、当時ヒト脳機能イメージング研究手法の主流であったポジトロン断層法(PET)による血流・代謝測定を置き換えた。現在では、脳機能イメージング研究の代名詞として、健常脳の機能分離や機能連関の理解、あるいは精神・神経疾患の病態生理の解明のため、欠かすことのできないツールとなっている。画像化していることには留意が必要である。脳活動に由来するfMRI信号の変動は、脳活動以外の要因による信号変動と比較して必ずしも大きくないため、興味のある脳活動を抽出するために適切な画像・信号処理を行うことも重要である。本項目では、脳機能を解明するツールとしてのfMRIの原理、解析法とそれらを応用した脳科学研究の潮流を概説する。脳血流(CBF)は動脈〜小動脈(arteriole)から流入し、動脈血中では赤血球のヘモグロビンは酸素化(oxy-Hb)されている。刺激がない安静時(左)であっても、酸素は脳の基礎代謝要求により消費される。酸素は毛細血管で脳組織に供給され、酸素化ヘモグロビン(oxy-Hb)は、常磁性体である還元ヘモグロビン(deoxy-Hb)に変わる。外的刺激などによりシナプス入力と神経活動が増加すると(右)、局所の酸素・エネルギー代謝要求は安静時と比較して増加する。酸素消費により還元ヘモグロビン(常磁性体)が増加するはずだから、局所磁場が乱れてT2*が短縮するように思われる(陰性BOLD信号)。神経血管単位は基礎代謝要求の増加を検知して動脈血の流入を要求量以上に増加させるらしい。この過程にはプロスタグランジン(PG)や一酸化窒素(NO)が関わっているとされる。これらの影響の総和として、脳活動が増加する部分ではdeoxy-Hbが相対的に薄まって局所磁場が安定し、T2*延長が観察される。多くのfMRI法ではこのT2*の延長を陽性BOLD信号として計測している。原理BOLD信号の発見1990年に小川誠二らはグラジエントエコー(GE)法MRIでマウス生体脳内の血管を撮像し、血管近傍のMRI信号が酸素飽和度に依存して変化することを発見した[1]。小川はこの観察について、酸素飽和度が減少する部位では、赤血球内の還元ヘモグロビン(deoxy-Hb)が増加し、近傍の脳実質の信号が変化するためであると考えた。還元ヘモグロビンは不対電子による常磁性(paramagnetic)を示し、還元ヘモグロビンの増加は局所磁場を乱し横磁化の減衰を増強するからである。コントラストと名付けた。小川らは、ヒト健常ボランティアに視覚刺激を与えると一次視覚野のMRI信号が増加することを示した。hyperemia)をもたらす。coupling)に依存している(図1)。2001年のLogothetisらによるサルを対象としたfMRIと神経活動の同時計測は、刺激や課題に伴う脳活動の増加とBOLD信号の関係性の理解に大きく貢献した[4]。LFP)は一過性の増加に引き続く持続性の増加を示した。BOLD信号変化は、マルチユニット神経活動よりもLFPと良く相関する。LFPはシナプス活動、すなわち神経細胞への情報入力を反映するから、BOLD信号変化は皮質からの出力量よりも皮質への入力量との関係が強いことになる。fMRIにより、皮質層ごとの情報処理を反映するBOLD信号変化の計測技術が精力的に開発されている[5]。脳活動とBOLD信号の相関関係は、外的刺激や課題の無い安静時にも観察される。state)の活動に大半が消費されている。これは脳が安静時にも組織的かつ活発な自発活動を示すことによる。両側外側頭頂葉などは、課題遂行時と比べて安静時にむしろ脳血流が増加することが知られていた[6]。DMN)と名づけられ、基底状態の脳の統合性に関わる内因性機構として提唱された。一方で、Biswalは1995年に安静状態のfMRIを解析し、両側運動感覚野の信号が主に0.1Hz以下の低い周波数(f)帯域において1/fのパターンで同期していることを見出していた[7]。FC)を反映すると考えられた。興味深いことに、安静時fMRIにより、Raichleの提唱したDMN内の脳領域間には強い機能結合が存在することがわかった[8]。このように2つの独立した研究の潮流が融合したことで、安静時fMRIを用いて局所の自発脳活動と領域間の機能結合状態を評価できる可能性に大きな注目が集まった。その後、安静時fMRIの信号同期性がサル脳における神経連絡性に対応していることも判明した[1]。現在、安静時fMRIの同期現象は、神経連絡を持つ遠隔領域間で同期して発生する自発性のシナプス・神経活動(及びこれらに伴うBOLD効果)に基づいていると考えられている。領域間の同期の詳細を知るための解析手法の改善、覚醒時の基底状態としての意識との関連、精神疾患や認知症などの病態との関連、神経連絡性との対応などについて研究が進んでいる。持続時間のごく短い感覚入力事象(t=0)に応答する血流動態応答関数(HRF)MRI撮像の繰り返し時間(TR)が一秒の場合を示す。感覚刺激後に生じるBOLD信号を経時的に計測すると、信号強度は刺激呈示から約2秒後に基準線(baseline)を超え、約5秒後に最大値を示す。HRF)と呼ぶ(図2)。dip)が観測されることもある。神経活動増加の後、信号強度は基準線より低下し、しばらくその状態が続く(undershoot)。dip)に関しては、血管拡張が最も早い皮質最深部の層では見られず、最も遅い最表層の皮質で見られることから、血流動態応答に先行して生じる酸素消費と還元ヘモグロビンの増加を反映する、という仮説が支持されている[9]。dip)のほうが、陽性BOLD反応よりも空間選択性が高いことを示している[10]。血流動態応答は、1)前述の神経血管連関に基づく現象であり、2)後述のようにBOLD信号の変化から課題遂行中の神経活動を推定する逆問題を解く場合に必要な応答関数であることから、その生物学的機構や応答関数の妥当性などについてさらなる解明が望まれる。fMRIを用いた脳機能解明法fMRI撮像・実験デザインfMRIの撮影には、BOLD信号のコントラスト源であるT2*減衰を鋭敏にとらえ、かつ時間分解能も高いエコープラナー画像(EPI)法を用いることがほとんどである。fMRIの性能を上げるためには、空間と時間両次元の分解能を高めることが理想的である。従って、撮像する際の時・空間解像度としては、MRI装置の性能を最大限に活かした高解像度が望ましいが、一方で解像度を上げすぎると信号雑音(S/N)比が低下することには注意が必要である。空間解像度と時間解像度もトレードオフの関係にあるが、マルチバンドEPIなどの高速撮像法の開発により長足の進歩があった。マルチバンドEPIは、従来のEPIに比べて空間解像度を1.2-1.5倍程度あるいは時間分解能を6-8倍近く向上させることができ、脳科学研究に汎用され始めている。fMRIは撮像中の研究参加者(あるいは動物)の状態により2種類に大別される。安静時fMRI閉眼条件あるいは眼前の十字を固視した開眼条件で、「何も考えない」ように指示して5-10分程度のfMRI撮像を行うことが多い。課題の練習を行う必要がない楽な検査法であり、健常者のみならず患者や小児の脳研究にも使いやすく、上述のデフォルトモードネットワークの評価にも応用可能である。課題fMRI数秒から数十秒単位の課題遂行を間歇的に行うよう指示をして実行してもらう。通常、コンピュータを用いて感覚刺激を提示したり応答を記録したりする。システムを整備する必要がある。課題fMRIでは、課題遂行の時間情報と前述の血流動態応答関数を畳み込んだBOLD信号変化モデルを作成し統計解析を行う(後述「fMRIの統計解析」参照)。課題遂行時間がやや長い場合(10-20秒以上)の場合は前述のデフォルトモードネットワークの周期や呼吸変動ノイズと周波数が近くなり、数秒以内の課題遂行の場合には体動や心拍の周波数帯と近くなる。そのためこうした関心の無い実験要因による解析精度の劣化の可能性を考慮した課題デザインを設計することが望ましいが、課題デザインについての詳細は本項目の範囲を超えるため別項目に譲りたい。fMRIデータの前処理3テスラMRI装置でのfMRIにおいて、神経活動による信号変化は全信号変動の約7パーセント程度とされる。生体や装置の熱ノイズが約50%を占め、残りはMRI装置の不完全性による信号ドリフトや、非撮像体の動き等に由来するアーチファクトとされる[11]。correction/realignment)することが通常であるが、この処理では補正しきれないMRI信号変動が残存する。加えて心拍や呼吸などの生理的要因による信号変動も存在するため、これら非神経性の信号を可能な限り分離・除去することが重要である(図3)。filter)をかけることでノイズ除去が行われてきた。この方法では、関心ある神経活動信号も除去してしまうため、独立成分分析によるノイズ同定・軽減への期待が高い[12](後述「fMRIの統計解析」参照)。correction/realignment)のみを行ったfMRI画像。動きはある程度補正されているものの、撮像中の頭部の動きは大きなMRI信号の変動を引き起こすためにノイズ軽減処理はまだ十分ではないことが見て取れる。動き補正に続き、独立成分分析によるノイズ軽減処理と脳外組織のマスク除去を行った後の画像。頭部の動きによる画像アーチファクトは目立たない。いずれも脳標準空間への位置合わせを行い、30倍速で表示している。機能的MRI画像の歪み補正の効果fMRI画像で用いるEPI画像は、主に撮像時の位相方向の設定と静磁場不均一性に依存して歪みが生じる。位相方向を前⇒後(矢印)、および後⇒前(矢印)に設定した場合のEPI画像。それぞれの画像が前後軸に沿ってに逆方向に歪んでいることがわかる。同じ位置での高解像度T1強調画像の断面像。歪み補正により画像内の各脳部位の位置がほぼ一致している。歪み補正により推定した位置ズレ(shift)の大きさを示す画像。歪み補正する前と後では同断面内で最大15mmの位置ズレがみられた。上述のようにfMRIデータにおける熱ノイズ量はBOLD信号変動に対して非常に大きい。smoothing)が行われる。この方法は位置情報の正確さを低下させるため、位置情報を犠牲にしないノイズ軽減法(例:ウィシャートフィルター)の開発も現在進んでいる([13]の図S7を参照)。EPI画像法には、生体との干渉により生じる静磁場(B0)の空間的不均一性により画像の歪み・位置ずれが生じる。そこで静磁場不均一性の空間分布を示す画像を別に撮像し、この情報を用いてfMRI画像の歪みを補正することで本来の位置を復元する(図4)[14]。mm程度発生し、fMRI画像を解剖画像に位置合わせする際に誤差を引き起こす。mm程度しかないことを考えると、正確な機能マッピングを行うためには位置ずれ補正が重要である。個人間の脳の形や大きさの違いによらない機能マッピングを行うため、3次元位置合わせ法を適用することで脳形状の標準化が行われる。z)上に脳活動を位置同定(マッピング)することで脳の機能構築の解明が進んできた。これまでの膨大なfMRI研究の成果は、MNI座標での機能マッピングデータベースとして構築され[15]、web上の視覚ツールも公開されている(neurosynth.org)。MNI座標系のような3次元空間での位置合わせ技術も、高次元の非線形法が進歩したことで精度が向上している。複雑な脳回のパターンを保ったまま3次元位置合わせを行うことは原理的に困難である。そこで、大脳皮質が2次元のシートが折れ曲がった構造をしていることに注目した皮質表面解析法が生まれた[16][17]。構造MRI画像から皮質を分画化し、その外表面(皮質と脳軟膜境界)と内表面(白質・皮質境界)を抽出し「皮質表面」を得る。registration)が提案され[18]、個人間の皮質機能構築を2次元座標上で標準化できるようになった。皮質表面の2次元座標系と、皮質下構造の3次元MNI座標系を組み合わせた新しい座標(灰白質座標)と専用フォーマットCIFTIも設計された。脳構造画像が必要であり、得られた表面境界とfMRI画像との正確な位置合わせを行うことで[19]、fMRI画像を灰白質画像に正確に重ね合わせ、個人間解析を高い位置精度で行うことが可能になってきている。脳回情報だけでなく機能や構造に関わる値(安静状態神経ネットワークやミエリンマップ等)を複数組み合わせた、さらに高い精度の表面位置合わせ法も開発され、応用も進んでいる[20]。課題fMRIの課題デザイン・解析・結果の例ヒト脳コネクトームプロジェクト(HCP)の課題fMRIデータを使用。)を動かす課題fMRIの結果。(HRF)を畳み込んでBOLD信号変化モデル(緑線)を作成する。データへのあてはめ(GLM)を行い、モデルとデータが有意に相関した部位(体素)の信号変化(赤線)を示している。運動課題のBOLD信号変化モデルを30例の被験者のデータにあてはめた際のt検定の結果(t値)を、大脳皮質表面に投影して表示。2次体性感覚野の神経活動が示唆される。analysis課題fMRIの解析には、課題に対応する脳活動とそれに応答するBOLD信号変化を合わせてモデル化し、この理論的な信号変化モデルと実際に測定したBOLD信号との適合度を体積画素/体素(voxel)ごとに計算し、課題に関連する脳活動を反映する統計画像を得る方法が広く用いられている。信号変化モデルは、複数の説明変数(課題遂行と相関する「脳活動」、頭部の動きなど脳活動以外の要因など)や脳活動に対するBOLD反応の時間的遅れやバラつき(血流動態関数、HRF)を考慮して構築する。GLM)により、信号変化モデルとデータとの適合程度を調べ、各説明変数の重み付け係数値画像や誤差値画像を算出する(第1段階解析)(図5A)。臨床検査としての課題fMRI(術前の言語機能評価など)のように、個人の結果に興味がある場合は、係数値画像を誤差値画像で割って得られるt値画像などを用いた個人内要因の統計検定を行う。課題Aと課題B遂行時の脳活動の差を見るためには、係数の引き算に相当する重みづけ操作を行い(”コントラストを立てる”)、得られたt値画像を用いて統計検定を行う。医学研究や神経科学研究では、測定対象とした複数の個人を母集団からランダムに抽出したサンプルと考え、母集団についての推定を行いたいことが一般的である。GLMM])を行い、課題と相関する脳活動とその群間差を示す脳部位を検出する(2段階解析)(図5B)。ここで問題となるのが統計学的有意性の閾値の決定の方法である。画像統計では数十~百万単位の体素を同時に扱うため、適切な多重比較補正を行わないと、脳科学的に意味のない信号変化を脳活動として検出されてしまう危険(偽陽性)を伴う。そのため、古典的には脳の中に脳活動領域としての塊(cluster)を想定し、そのような塊が偶然形成される確率から閾値を設定する理論が採用されている[21]。この統計閾値法の妥当性については議論が続いており[22]、課題fMRIの再現性が低い主要因になっているとの見方がある[23]。testなど[24]。MVPA一般線形モデル(GLM)では体素ごとの統計検定を行うのに対し、MVPAでは複数の体素の信号値の関係性(パターン)が課題の差や、疾患と健常の差を表現すると考えた多変量解析を行う。MVPAは「課題に対する脳活動や疾患に特有の脳活動の時空間的なパターンを判別」するアプローチとも言える。初期にはGLMによる重み付け画像同士の相関性がパターンの類似性指標として使われた[25]。SVM)などの機械学習技術が判別法に採用されている。脳活動の空間的なパターンを情報として、与えられた課題(例えば視覚刺激)の逆問題を推定または解読(decoding)する研究も行われている[26]。安静時fMRIから得られた多変量の機能結合パターンを、機械学習技術を応用して判別し、精神神経疾患の診断補助を目指す先駆的な医療技術の開発も進んでいる[27]。analysisシステムとしての脳の仕組みを探るにあたり、機能局在を追求することだけでは限界がある。そこで脳部位間の繋がりが注目され、2005年には脳の網羅的な結線図を意味するコネクトーム(connectome)という概念が提唱された[28]。脳コネクトームには様々な空間スケールがありうる。正方行列を計算することで、全脳レベル(macroscale)の機能的コネクトーム解析が可能である。FCの産出法としては領域間の時系列信号のピアソン相関係数が良く使われる[7]。ある特定の脳領域を関心領域(シード)として設定し、相手方の関心領域(または体素)を、全脳に渡って隙間なく設定した相関解析を行えば、特定のシードに対する脳全体の機能結合マップが得られる。安静時fMRI解析において後部帯状回に関心領域を設定すれば、PETで見られたデフォルトモードネットワークが安静状態神経ネットワークとして抽出できる[29][30]。シードを全脳の関心領域すべてに設定することで、多対多の網羅的な機能結合、すなわち機能的コネクトームが算出できる。correlation)を用いることで、2領域間に固有性の高いFCを評価できることが示唆されている[31]。)、他の指標や測定(脳波や臨床兆候など)との関連付け、疾患バイオマーカーの探索等の様々な研究も行われている。一方でヒトの脳コネクトーム研究単独では結果の妥当性の検証が難しいため、動物脳での検証も重要であり、例えば神経連絡トレーサーと機能的結合の比較は重要な課題である。含むデフォルトモードネットワーク(DMN)。左上は膨らました皮質表面にマッピングしたもの、右上は軸断像に表示、および同ネットワークの信号変化(左下)および周波数分析結果(右下)。周波数0.01-0.1Hzに高いパワーを持つネットワーク活動である。図3と同じfMRIデータで、動き補正・ノイズ処理を行ったあとに得られたもの。analysis独立成分分析はデータ駆動型の解析でありブラインド信号分離の技術として多くの分野で活用されている。PCA)と比べ、データマイニングの手法として柔軟であると考えられている。主成分分析では無相関、つまりベクトルとして直交する成分のセットが解析的に一意に定まる。これに対し独立成分分析では、互いに他の情報を持たないという特性(独立性)が大きくなるよう成分を分離する。無相関であっても独立ではないことがあるから、独立性は相関性よりも柔軟な条件である。より数十個の独立成分を求めるものであった。神経解剖学的に解釈可能な皮質ネットワークを得ることができデフォルトモードネットワークの同定にも使用される(図6)。独立成分分析は従前から脳波など電気生理データのノイズ除去に利用されてきたが、機械学習判別によって個人レベルのfMR元データから構造ノイズ成分を除去する前処理過程にも独立成分分析が利用されている(fMRIデータの前処理の項参照)。適応することで、個人間に共通のネットワークを同定し(グループ独立成分分析)、グループとしての安静状態神経ネットワークを検出することも可能である。このグループRSNの各個人のデータへのあてはめ(一般線形モデル)を、時間と空間について2回行うことで、個人毎のRSNが得られる。この個人毎の安静状態神経ネットワークを用いた集団解析が可能であり[32]、病態解明への応用研究も進んでいる。最近は時間分解能の高いfMRIデータが取得されるようになってきたことで、時間的ICAの応用も進み、空間的に重なり合う成分や神経連絡性とも合致する成分の抽出[33]や、脳全体の広範なノイズの検出に有用であることが分かってきている。独立成分分析の結果の解釈には数学的な仮定(例えば次元数や線形性等)と脳のシステム・神経-血管連関・画像ノイズの妥当性や検証に留意する必要がある[34]。)を用いて、皮質表面上で機能構築を位置合わせし、皮質表面上で急激に変化する(傾斜が高い)部位を分画の境界線として設定して作成された[35]。本分画は一般公開されており誰でも利用できる[1]1990年後半からfMRI研究は飛躍的に世界中に広まり多くの成果をあげてきた。その間に基礎となるMRI技術も大きく進展し、質の高いデータが得られるようになった。初期の研究の興味は脳の機能局在・分離の同定にあったが、2010年代以後はネットワークとしての脳の解明に興味がシフトした[36]。一つの研究が扱うfMRIデータの数も、数10人から数百人、数万人の規模へと拡大し、種としてのヒトの脳機能構築やその複雑性の理解、個体差の理解へと興味が広がった。Project(HCP)が行われた。このプロジェクトでは、ワシントン大学セントルイス医学校およびオックスフォード大学が中心となり約1200名の健康な若年成人被験者を対象として、安静時fMRI、標準的な課題を用いた課題fMRI、構造MRI(T1・T2強調画像)、拡散強調画像を撮像し、脳内の機能的結合と構造的結合について統合的に解析を進めた。このプロジェクトは(1)質が高い画像データを大量に取得する、(2)空間分解能を犠牲にすることなくデータを処理する手法を開発する、(3)FACT、すなわち皮質機能(Function)・連絡性(Connectivity)、構造(Architecture)、位置(Topography)の情報を統合した脳の領域分割(parcellation)を行う、(4)取得した生データ、解析データ、データ解析に必要なプログラム・コードを無料で公開する、といった多くの野心的な目標を達成し、現在も世界のMRI脳科学研究に大きな影響を与え続けている。FACT法は片側半球を180領域に分割することに成功し[35]、100年以上続くヒト脳の機能分画の歴史の中で初めて非侵襲手法により生きた個人の脳地図を作成することに成功した(図7)。国際のライセンスで利用できます。脳科学辞典について免責事項モバイルビュ。