「近赤外線スペクトロスコピー」の版間の差分

編集の要約なし
編集の要約なし
編集の要約なし
1行目: 1行目:
英語名:Near-infrared spectroscopy (NIRS)
英語名:Near-infrared spectroscopy (NIRS)


 近赤外線スペクトロスコピー(near-infrared spectroscopy、 NIRS)は、[[wikipedia:JA:近赤外|近赤外]]領域の光を物質に照射し、透過してきた光の性質(強度など)を解析して、非破壊的に対象物の構成成分を分析する方法で、[[wikipedia:JA:食品科学|食品科学]]や[[wikipedia:JA:農業|農業]]など様々な領域で用いられている。生体への応用は、1977年に[[wikipedia:JA:デュ-ク大学|デュ-ク大学]]の[[wikipedia:Jöbsis|Jöbsis]]が近赤外光を用いて動物の[[wikipedia:JA:心臓|心臓]]や[[脳]]の[[wikipedia:JA:酸素化|酸素化]]状態を非侵襲的に計測したのが始めてで<ref name=ref1><pubmed>929199</pubmed></ref>、以後、生体組織における[[wikipedia:JA:血流|血流]]・酸素代謝モニタ法として研究・開発が進められてきた。さらに、1990年代になってNIRSが神経活動に連動した[[脳血流]]変化に伴う[[wikipedia:JA:ヘモグロビン|ヘモグロビン]](Hb)変化をとらえることができるということが相次いで報告され、本法は新しい脳機能イメージング法(functional NIRS、 fNIRS)としても注目されるようになった。NIRSの応用例については専門誌の特集号<ref name=ref1><pubmed>929199</pubmed></ref><ref name=ref2>Illuminating the future of biomedical optics.<br>Phil Trans R Soc A. 2011, 369</ref>などを参照していただき、ここではfNIRSを中心に基礎的事項を解説する。
 近赤外線スペクトロスコピー(near-infrared spectroscopy、 NIRS)は、[[wikipedia:JA:近赤外|近赤外]]領域の光を物質に照射し、透過してきた光の性質(強度など)を解析して、非破壊的に対象物の構成成分を分析する方法で、[[wikipedia:JA:食品科学|食品科学]]や[[wikipedia:JA:農業|農業]]など様々な領域で用いられている。生体への応用は、1977年に[[wikipedia:JA:デュ-ク大学|デュ-ク大学]]の[[wikipedia:Jöbsis|Jöbsis]]が近赤外光を用いて動物の[[wikipedia:JA:心臓|心臓]]や[[脳]]の[[wikipedia:JA:酸素化|酸素化]]状態を非侵襲的に計測したのが始めてで<ref name=ref1><pubmed>929199</pubmed></ref>、以後、生体組織における[[wikipedia:JA:血流|血流]]・酸素代謝モニタ法として研究・開発が進められてきた。さらに、1990年代になってNIRSが神経活動に連動した[[脳血流]]変化に伴う[[wikipedia:JA:ヘモグロビン|ヘモグロビン]](Hb)変化をとらえることができるということが相次いで報告され、本法は新しい脳機能イメージング法(functional NIRS、 fNIRS)としても注目されるようになった。NIRSの応用例については専門誌の特集号<ref name=ref2>Illuminating the future of biomedical optics.<br>Phil Trans R Soc A. 2011, 369</ref><ref name=ref3>Special Issue on Medical Applications.<br>J Near Infrared Spectrosc. 2012, 20 (1).</ref>などを参照していただき、ここではfNIRSを中心に基礎的事項を解説する。




11行目: 11行目:
[[image:図2頭部における光伝播と拡張ベア・ランバート則.png|thumb|300px|'''図2.頭部における光伝播と拡張ベア・ランバート則'''<br>総光路長(L)は光路全体の長さで、散乱により照射・受光間距離ρより長い(βは未知数)。局所脳活動により血流が変化した領域を赤で示した。その領域でのHb濃度変化(ΔC)に由来する吸光度変化(ΔA)は、ΔCと[[wikipedia:JA:モル吸光係数|モル吸光係数]](ε)と部分光路長(partial PL)の積で表わされる。]]
[[image:図2頭部における光伝播と拡張ベア・ランバート則.png|thumb|300px|'''図2.頭部における光伝播と拡張ベア・ランバート則'''<br>総光路長(L)は光路全体の長さで、散乱により照射・受光間距離ρより長い(βは未知数)。局所脳活動により血流が変化した領域を赤で示した。その領域でのHb濃度変化(ΔC)に由来する吸光度変化(ΔA)は、ΔCと[[wikipedia:JA:モル吸光係数|モル吸光係数]](ε)と部分光路長(partial PL)の積で表わされる。]]


 近赤外光(通常用いられているのは波長700~900 nmの光)は、生体に対して高い透過性を示す。これは、近赤外光は可視光に比べて散乱されにくく、生体物質でこの領域の光を吸収するのは、血液中のヘモグロビン(Hb)や筋肉中の[[wikipedia:JA:ミオグロビン|ミオグロビン]](Mb)、そして[[wikipedia:JA:ミトコンドリア|ミトコンドリア]]内にある[[wikipedia:JA:チトクロームCオキシダーゼ|チトクロームCオキシダーゼ]](cyt. ox.)などに限られているためである。NIRSは、近赤外光の吸収される度合がHb、 Mbの酸素化-脱酸素化状態、cyt. ox.の酸化-還元状態によって異なることを利用して、生体に近赤外光を照射して照射点から数センチメートル離れたところで体外に現れた光を検出し、その性状を解析してHb濃度変化などを求めるが、単一の方法ではなく複数の異なる計測方法がある(図1)。その中で最も一般的なのは、連続光(continuous wave light)を用いて拡張ベア・ランバート則<ref name=ref4><pubmed>3070170</pubmed></ref>に基づいて、Hbなど光を吸収する物質の濃度変化を求める方法でありCW計測と呼ばれている。拡張ベア・ランバート則は、A(吸光度)= -logI/I<sub>0</sub> = CL + S で表される(図2)。ここでI0とIは照射ならびに検出光の光量、はモル吸光係数、Cは光吸収物質の濃度、Lは照射された光が検出されるまでに通った生体内における経路の長さ(光路長;個々の光子は違う経路を通るのでそれらの平均光路長)、Sは主として散乱による光の減衰を示す項で通常定数と見なされている。生体計測では、複数の波長を用いて酸素化Hb (oxy-Hb)、脱酸素化Hb (deoxy-Hb)、両者の和である総Hb (t-Hb)の濃度変化を求めるが、光路長を計測することができないため、得られる信号は濃度変化と光路長の積である。しかし、装置によって用いられている波長や演算式が異なり、さらに演算式によってはモル吸光係数も未知数として取り扱っているものもあり、NIRS信号は物理単位をもたない。従って、NIRS信号は[mM・mm](濃度×長さ)の単位で表現されることもあるが、単位をもたないあるいは任意単位(au)として表現されるのが適切である。
 近赤外光(通常用いられているのは波長700~900 nmの光)は、生体に対して高い透過性を示す。これは、近赤外光は可視光に比べて散乱されにくく、生体物質でこの領域の光を吸収するのは、血液中のヘモグロビン(Hb)や筋肉中の[[wikipedia:JA:ミオグロビン|ミオグロビン]](Mb)、そして[[wikipedia:JA:ミトコンドリア|ミトコンドリア]]内にある[[wikipedia:JA:チトクロームCオキシダーゼ|チトクロームCオキシダーゼ]](cyt. ox.)などに限られているためである。NIRSは、近赤外光の吸収される度合がHb、 Mbの酸素化-脱酸素化状態、cyt. ox.の酸化-還元状態によって異なることを利用して、生体に近赤外光を照射して照射点から数センチメートル離れたところで体外に現れた光を検出し、その性状を解析してHb濃度変化などを求めるが、単一の方法ではなく複数の異なる計測方法がある(図1)。その中で最も一般的なのは、連続光(continuous wave light)を用いて拡張ベア・ランバート則<ref name=ref4><pubmed>3070170</pubmed></ref>に基づいて、Hbなど光を吸収する物質の濃度変化を求める方法でありCW計測と呼ばれている。拡張ベア・ランバート則は、A(吸光度)= -logI/I<sub>0</sub> = CL + S で表される(図2)。ここでI<sub>0</sub>とIは照射ならびに検出光の光量、はモル吸光係数、Cは光吸収物質の濃度、Lは照射された光が検出されるまでに通った生体内における経路の長さ(光路長;個々の光子は違う経路を通るのでそれらの平均光路長)、Sは主として散乱による光の減衰を示す項で通常定数と見なされている。生体計測では、複数の波長を用いて酸素化Hb (oxy-Hb)、脱酸素化Hb (deoxy-Hb)、両者の和である総Hb (t-Hb)の濃度変化を求めるが、光路長を計測することができないため、得られる信号は濃度変化と光路長の積である。しかし、装置によって用いられている波長や演算式が異なり、さらに演算式によってはモル吸光係数も未知数として取り扱っているものもあり、NIRS信号は物理単位をもたない。従って、NIRS信号は[mM・mm](濃度×長さ)の単位で表現されることもあるが、単位をもたないあるいは任意単位(au)として表現されるのが適切である。


=== 生体における光伝搬特性 ===  
=== 生体における光伝搬特性 ===  
31行目: 31行目:
=== 神経機能イメージング ===  
=== 神経機能イメージング ===  


 NIRSは脳の表層レベルの情報しか得ることができないが、装置が簡便で測定時に厳しい体動制限を必要とせず、また日常生活と同じような環境で計測することができる。さらに、 NIRSは[[wikipedia:JA:時間分解能|時間分解能]]が高く(1秒以下)、脳血流の変化をリアルタイムかつ経過時的に追跡することができるなどの利点を持つ。このようなNIRS の利点は、 物理的(運動中など)・生理的(乳幼児など)に[[PET]]や[[fMRI]]計測が難しい対象での神経機能イメージング研究を可能にし、認知脳科学にとどまらず、工学、教育学など様々な領域で用いられている。
 NIRSは脳の表層レベルの情報しか得ることができないが、装置が簡便で測定時に厳しい体動制限を必要とせず、また日常生活と同じような環境で計測することができる。さらに、 NIRSは[[wikipedia:JA:時間分解能|時間分解能]]が高く(1秒以下)、脳血流の変化をリアルタイムかつ経時的に追跡することができるなどの利点を持つ。このようなNIRS の利点は、 物理的(運動中など)・生理的(乳幼児など)に[[PET]]や[[fMRI]]計測が難しい対象での神経機能イメージング研究を可能にし、認知脳科学にとどまらず、工学、教育学など様々な領域で用いられている。


== 課題 ==
== 課題 ==
39行目: 39行目:
 上で述べたように、一般に市販されている装置(CW計測装置)から得られる信号は、 Hb濃度変化に光路長が掛け合わさった値で、脳賦活時のように血流変化が脳内に限局している場合には、血流変化の生じている部分の光路長(部分光路長)との積になる(図2)。CW計測による[[光トポグラフィ]]では、多チャンネル装置を用いて複数の領域を計測し、各計測領域の光路長は照射-受光間距離が同じなら一定であると仮定して、得られる画像はHbの濃度変化のマッピングを示していることになっている。しかし、光路長は照射-受光間距離が同じでも計測部位によって異なり、さらに、部分光路長と総光路長は負の関係にあるため(照射-受光間距離が一定の場合、総光路長が長いほど部分光路長は短い<ref name=ref6><pubmed>16409097</pubmed></ref>)、信号の振幅は必ずしもHb濃度変化の大きさを示しているわけではなく、計測領域全体でHb濃度が一様に変化した場合は、光路長のマッピングを示している可能性もある。しかし、図1に示した[[時間分解計測法]](time-resolved spectroscopy, TRS)や[[位相分解計測法]](phase-resolved spectroscopy, PRS)を用いて総光路長を計測することはできても、現時点では部分光路長を実測することは不可能であるため、拡張ベア・ランバート則に基づくCW計測で、脳活動に連動するHb濃度変化を定量的に算出することはできない。さらに、NIRS計測では照射と受光ファイバペアの位置によって、信号の振幅は異なるため、NIRS信号の振幅の部位間・個体間比較は、血流反応の大小の比較にはならない。
 上で述べたように、一般に市販されている装置(CW計測装置)から得られる信号は、 Hb濃度変化に光路長が掛け合わさった値で、脳賦活時のように血流変化が脳内に限局している場合には、血流変化の生じている部分の光路長(部分光路長)との積になる(図2)。CW計測による[[光トポグラフィ]]では、多チャンネル装置を用いて複数の領域を計測し、各計測領域の光路長は照射-受光間距離が同じなら一定であると仮定して、得られる画像はHbの濃度変化のマッピングを示していることになっている。しかし、光路長は照射-受光間距離が同じでも計測部位によって異なり、さらに、部分光路長と総光路長は負の関係にあるため(照射-受光間距離が一定の場合、総光路長が長いほど部分光路長は短い<ref name=ref6><pubmed>16409097</pubmed></ref>)、信号の振幅は必ずしもHb濃度変化の大きさを示しているわけではなく、計測領域全体でHb濃度が一様に変化した場合は、光路長のマッピングを示している可能性もある。しかし、図1に示した[[時間分解計測法]](time-resolved spectroscopy, TRS)や[[位相分解計測法]](phase-resolved spectroscopy, PRS)を用いて総光路長を計測することはできても、現時点では部分光路長を実測することは不可能であるため、拡張ベア・ランバート則に基づくCW計測で、脳活動に連動するHb濃度変化を定量的に算出することはできない。さらに、NIRS計測では照射と受光ファイバペアの位置によって、信号の振幅は異なるため、NIRS信号の振幅の部位間・個体間比較は、血流反応の大小の比較にはならない。


=== 選択的Hb濃度計測 ===  
=== 選択的脳内Hb濃度計測 ===  


 頭部での計測において、皮膚血流に変化が生じた場合、それがNIRS信号に影響することは必至であるが、顔面・頭皮の皮膚血流制御機構が複雑であるため、血流変化がどのような場合に生じ、その影響がどの程度であるのかについて予測することは難しい。これまでに皮膚血流の影響を取り除く方法は幾つか提案されており、光の照射点を起点にして直線上に複数の異なる照射-受光間距離で計測するマルチディスタンス法は有効で、後述の拡散理論に基づいて光学特性値(吸収係数、換算散乱係数)を求める空間分解計測法を利用して、脳組織のHbの酸素飽和度を算出する装置は市販されている。また、得られたデータから主成分分析や独立成分分析法を用いて体循環などに起因する変化を取り除く方法も考案されているが、装置に組み込まれている方法はまだ少ない。
 頭部での計測において、皮膚血流に変化が生じた場合、それがNIRS信号に影響することは必至であるが、顔面・頭皮の皮膚血流制御機構が複雑であるため、血流変化がどのような場合に生じ、その影響がどの程度であるのかについて予測することは難しい。これまでに皮膚血流の影響を取り除く方法は幾つか提案されており、光の照射点を起点にして直線上に複数の異なる照射-受光間距離で計測するマルチディスタンス法は有効で、後述の拡散理論に基づいて光学特性値(吸収係数、換算散乱係数)を求める空間分解計測法を利用して、脳組織のHbの酸素飽和度を算出する装置は市販されている。また、得られたデータから主成分分析や独立成分分析法を用いて体循環などに起因する変化を取り除く方法も考案されているが、装置に組み込まれている方法はまだ少ない。
48行目: 48行目:
NIRS信号強度の群間比較(たとえば患者群と健常群)に意味がないという事実は、本法が臨床の場に普及する妨げになってきたが、カテゴリカルデザインを用いた群間比較は可能である。たとえば、患者群と健常群において、タスク間比較(コントロールタスクとテストタスクなど)の結果に違いを認めるかどうかという様な比較は行うことができる。また、事象関連デザインによる解析は信号の絶対値を必要とせず、時間分解能がすぐれているNIRSには適した解析法と思われる。しかし、fMRIで用いられている血流動態反応関数(hemodynamic response function、 HRF)を、そのままNIRSの解析に用いることの妥当性などは、今後さらに検討が必要である。
NIRS信号強度の群間比較(たとえば患者群と健常群)に意味がないという事実は、本法が臨床の場に普及する妨げになってきたが、カテゴリカルデザインを用いた群間比較は可能である。たとえば、患者群と健常群において、タスク間比較(コントロールタスクとテストタスクなど)の結果に違いを認めるかどうかという様な比較は行うことができる。また、事象関連デザインによる解析は信号の絶対値を必要とせず、時間分解能がすぐれているNIRSには適した解析法と思われる。しかし、fMRIで用いられている血流動態反応関数(hemodynamic response function、 HRF)を、そのままNIRSの解析に用いることの妥当性などは、今後さらに検討が必要である。


 fMRIやPET計測においては、複数の被験者データを同一空間(標準脳)で解析することがルーチンに行われている。ある個人に特異的な所見ではなく、多くの人に共通に認められる現象を見出すためであるが、マルチチャンネルNIRS装置で得られたデータ解析において、最も困難な解析法の一つである。NIRS計測の前あるいは後で頭部MRI画像を取得して、デジタイザーで光ファイバの位置を脳表に投射してから、解剖学的正規化と空間的平滑化を行って標準脳に変形させることは可能である。脳波電極の配置(たとえば国際10-20法)は、その直下の脳領域を推定するのにしばしば利用されるが<ref name=ref11><pubmed>14741647</pubmed></ref>、たとえ脳波電極設置部位をランドマークにして光ファイバ固定用ホルダーを設置しても、各照射-受光ペアが計測する脳領域は被験者ごとに異なる。これは、ホルダーの大きさは同じでも、頭の大きさは個人で異なるためで、さらにホルダーをとりつける時に多少ずれてしまうこともある。
 fMRIやPET計測においては、複数の被験者データを同一空間(標準脳)で解析することがルーチンに行われている。ある個人に特異的な所見ではなく、多くの人に共通に認められる現象を見出すためであるが、マルチチャンネルNIRS装置で得られたデータ解析において、最も困難な解析法の一つである。NIRS計測の前あるいは後で頭部MRI画像を取得して、デジタイザで光ファイバの位置を脳表に投射してから、解剖学的正規化と空間的平滑化を行って標準脳に変形させることは可能である。脳波電極の配置(たとえば国際10-20法)は、その直下の脳領域を推定するのにしばしば利用されるが<ref name=ref11><pubmed>14741647</pubmed></ref>、たとえ脳波電極設置部位をランドマークにして光ファイバ固定用ホルダーを設置しても、各照射-受光ペアが計測する脳領域は被験者ごとに異なる。これは、ホルダーの大きさは同じでも、頭の大きさは個人で異なるためで、さらにホルダーをとりつける時に多少ずれてしまうこともある。


 そこで、標準脳に変形させてから、計測領域がほぼ同じとみなせる被験者のデータのみを用いて解析することは可能であるが、賢明な解決法とは言いがたい。脳磁場計測(MEG)や脳波計測(EEG)もNIRSと同様に標準的な集団解析の方法はなく、個人解析がまず行われ、その結果が集団解析の対象になっている。NIRS信号の解析においても、直接集団データの中から共通項を見いだすのではなく、MEGやEEGのように個人データ解析から行うのが確実であると考える。
 そこで、標準脳に変形させてから、計測領域がほぼ同じとみなせる被験者のデータのみを用いて解析することは可能であるが、賢明な解決法とは言いがたい。脳磁場計測(MEG)や脳波計測(EEG)もNIRSと同様に標準的な集団解析の方法はなく、個人解析がまず行われ、その結果が集団解析の対象になっている。NIRS信号の解析においても、直接集団データの中から共通項を見いだすのではなく、MEGやEEGのように個人データ解析から行うのが確実であると考える。
54行目: 54行目:
== 選択的・定量的脳内Hb計測 ==
== 選択的・定量的脳内Hb計測 ==


[[image:図3TRSで計測された検出光の時間プロファイル.png|thumb|300px|'''図3.TRSで計測された検出光の時間プロファイル'''<br>重心に対応する時間]]
[[image:図3TRSで計測された検出光の時間プロファイル.png|thumb|300px|'''図3.TRSで計測された検出光の時間プロファイル'''<br>t<sub>m</sub>重心に対応する時間]]


=== 時間分解計測法(time-resolved spectroscopy、 TRS)===  
=== 時間分解計測法(time-resolved spectroscopy、 TRS)===  


 TRSでは、ピコ秒オーダーの極短パルス光を生体に照射し、照射部位から数センチメータ離れた体表に現れた光を時間分解計測することによって、時間に対する検出光の強度分布(検出された[[wikipedia:JA:光子|光子]]数の時間プロファイル)をもとめる(図3)。比較的直進に近い状態で生体内を透過した光子は早い時間成分になり、散乱を強く受けた光子は遅い時間成分になる。つまり、図3は光子の飛行時間分布で、ここから平均飛行時間(t_m:時間プロファイルの重心に対応する時間)が決まり、生体内における光速をt_mで乗ずることによって平均総光路長を求めることができる。また、この時間プロファイルに生体における光伝搬モデルとしてしばしば用いられる光拡散方程式の解析解から得られた検出光強度分布をカーブフィッティングして吸収係数を決定し、吸収係数からHb濃度を算出することができる。このようにして求めたHb濃度は、CW計測に比べてより選択的に脳内Hb濃度変化を計測することができる。しかし、この方法を用いても脳外組織の影響を完全に取り除くことはできない。
 TRSでは、ピコ秒オーダーの極短パルス光を生体に照射し、照射部位から数センチメートル離れた体表に現れた光を時間分解計測することによって、時間に対する検出光の強度分布(検出された[[wikipedia:JA:光子|光子]]数の時間プロファイル)をもとめる(図3)。比較的直進に近い状態で生体内を透過した光子は早い時間成分になり、散乱を強く受けた光子は遅い時間成分になる。つまり、図3は光子の飛行時間分布で、ここから平均飛行時間(t<sub>m</sub>:時間プロファイルの重心に対応する時間)が決まり、生体内における光速をt<sub>m</sub>で乗ずることによって平均総光路長を求めることができる。また、この時間プロファイルに生体における光伝搬モデルとしてしばしば用いられる光拡散方程式の解析解から得られた検出光強度分布をカーブフィッティングして吸収係数を決定し、吸収係数からHb濃度を算出することができる。このようにして求めたHb濃度は、CW計測に比べてより選択的に脳内Hb濃度変化を計測することができる。しかし、この方法を用いても脳外組織の影響を完全に取り除くことはできない。


=== 拡散光トモグラフィ(Diffuse optical tomography、 DOT)<ref name=ref12><pubmed>18163810</pubmed></ref> ===  
=== 拡散光トモグラフィ(Diffuse optical tomography、 DOT)<ref name=ref12><pubmed>18163810</pubmed></ref> ===