JP5996625B2 - 単一粒子検出を用いた光分析装置、光分析方法及び光分析用コンピュータプログラム - Google Patents

Description

本発明は、共焦点顕微鏡又は多光子顕微鏡の光学系などの溶液中の微小領域からの光が検出可能な光学系を用いて、溶液中に分散又は溶解した原子、分子又はこれらの凝集体（以下、これらを「粒子」と称する。）、例えば、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞などの粒子状の対象物、或いは、非生物学的な粒子を検出して、それらの状態（相互作用、結合・解離状態など）の分析又は解析に於いて有用な情報を取得することが可能な光分析技術に係り、より詳細には、上記の如き光学系を用いて単一の粒子の存在による光の変化（単一粒子が発する光又は単一粒子が生ずる影）を個別に検出して種々の光分析を可能にする光分析装置、光分析方法並びに光分析用コンピュータプログラムに係る。なお、本明細書に於いて、検出される光は、蛍光、りん光、化学発光、生物発光、散乱光等であってよい。光を発する粒子（以下、「発光粒子」と称する。）は、それ自身が光を発する粒子、又は、任意の発光標識若しくは発光プローブが付加された粒子のいずれであってもよい。

近年の光計測技術の発展により、共焦点顕微鏡の光学系とフォトンカウンティング（１光子検出）も可能な超高感度の光検出技術とを用いて、一光子又は蛍光一分子レベルの微弱光の検出・測定が可能となっている。そこで、そのような微弱光の計測技術を用いて、生体分子等の特性、分子間相互作用又は結合・解離反応の検出を行う光分析技術が種々提案されている。そのような光分析技術としては、例えば、蛍光相関分光分析（Fluorescence Correlation Spectroscopy：ＦＣＳ。例えば、特許文献１−３、非特許文献１−３参照）、蛍光強度分布分析(Fluorescence-Intensity Distribution Analysis：ＦＩＤＡ。例えば、特許文献４、非特許文献４)やフォトンカウンティングヒストグラム（Photon Counting Histogram：ＰＣＨ。例えば、特許文献５）などが知られている。また、特許文献６〜８には、共焦点顕微鏡の光学系を用いて計測される試料溶液の蛍光信号の時間経過に基づいて蛍光性物質を検出する方法が提案されている。

更に、本願出願人は、特許文献９〜１１に於いて、共焦点顕微鏡又は多光子顕微鏡の光学系などの溶液中の微小領域からの光が検出可能な光学系を用いた光分析技術であって、ＦＣＳ、ＦＩＤＡ等の光分析技術とは異なる原理による新規な光分析技術を提案した。上記のＦＣＳ、ＦＩＤＡ等の光分析技術の場合、試料溶液内に於ける光の検出領域である微小領域（以下、「光検出領域」と称する。）に浮遊する蛍光分子からの光を連続的に計測することにより得られた光強度データが統計的なゆらぎを算出する演算処理により解析され、蛍光分子の濃度及び／又はその他の特性が検出される。これに対し、特許文献９〜１１に於いて提案された新規な光分析技術では、試料溶液内に於いて光検出領域の位置を移動させながら、即ち、光検出領域により試料溶液内を走査しながら、光検出領域が試料溶液中に分散してランダムに運動する発光粒子を包含したときに、その発光粒子から発せられる光を個別に検出し、これにより、試料溶液中の発光粒子を一つずつ検出して、発光粒子のカウンティングや試料溶液中の発光粒子の濃度又は数密度に関する情報の取得が可能となる。この光分析技術（以下、「走査分子計数法」と称する。）によれば、測定に必要な試料は、ＦＣＳ、ＦＩＤＡ等の光分析技術と同様に微量（例えば、数十μＬ程度）であってもよく、また、測定時間が短く（一回の測定で秒オーダーの時間の計測が数回繰り返される。）、更に、観測対象となる粒子の濃度が、ＦＣＳ、ＦＩＤＡ等の光分析技術の良好に測定可能なレベル（１ｎＭ程度）よりも低い試料溶液に於いて、発光粒子の存在を検出し、その濃度、数密度又はその他の特性を定量的に検出することが可能となる。

かくして、「走査分子計数法」は、医学・生物学の研究開発の分野でしばしば使用される希少な或いは高価な試料についての分析を行う場合や、病気の臨床診断や生理活性物質のスクリーニングなど、検体数が多い場合に、従前の生化学的方法に比して、低廉に、或いは、迅速に実験又は検査が実行でき、しかも、ＦＣＳ、ＦＩＤＡ等を良好に実行できない程度のより低い濃度の粒子の濃度及び／又は特性の検出が可能な強力なツールとして期待される。

特開２００５−０９８８７６ 特開２００８−２９２３７１ 特開２００９−２８１８３１ 特許第４０２３５２３号 国際公開２００８−０８０４１７ 特開２００７−２０５６５ 特開２００８−１１６４４０ 特開平４−３３７４４６号公報 国際公開第２０１１／１０８３６９ 国際公開第２０１１／１０８３７０ 国際公開第２０１１／１０８３７１

金城政孝、蛋白質 核酸 酵素 Vol.４４、No.９、１４３１−１４３８頁 １９９９年 エフ・ジェイ・メイヤー・アルムス（F.J.Meyer-Alms）、フルオレセンス・コリレーション・スペクトロスコピー（Fluorescence Correlation Spectroscopy）、アール・リグラー編（R.Rigler）、スプリンガー（Springer）、ベルリン、２０００年、２０４−２２４頁 加藤則子外４名、遺伝子医学、Vol.６、No.２、２７１−２７７頁 カスク他３名、米国科学アカデミー紀要 １９９９年、９６巻、13756‐13761頁（P. Kask, K. Palo, D. Ullmann, K. Gall PNAS 96, 13756-13761 (1999)） Anal. Chem. Vol.70 No.3 ６３２頁 １９９８年

上記の特許文献９〜１１に記載の走査分子計数法では、典型的には、光検出領域からの光強度の時系列データ（時系列光強度データ）に於いて、釣鐘状の光強度の変化が観測対象となる単一粒子の存在に対応するとの仮定の下、かかる釣鐘状の光強度の変化の検出が行われる。例えば、観測対象粒子が発光粒子である場合には（後述の如く、観測対象粒子が検出波長帯域の光を発しない粒子（非発光粒子）である場合もあり得る。）、時系列光強度データに観測された釣鐘状又はパルス状の光強度の増大のうち、発光粒子が光検出領域内を通過した際に想定される光強度の増大のプロファイル（ピーク強度、半値全幅等）の条件に適合する光強度の変化が粒子の存在を表す信号として検出される。そして、発光粒子の光強度の増大のプロファイルの条件に適合しない光強度の増大は、ノイズとして判定される。

しかしながら、時系列光強度データに於いて上記の如き釣鐘状の信号を検出する場合、粒子の信号の大きさが弱くなると、粒子の信号とノイズ信号との識別が困難となる。

後述の実施形態の欄に於いてより詳細に説明される如く、光強度の測定は、観測対象粒子による微弱な光強度又はその変化を捉えるために、典型的には、フォトンカウンティングにより実行され、時系列光強度データは、離散的なフォトンカウンティングデータとなる（図２（Ａ）参照）。その場合、釣鐘状の信号の検出を容易にするために、好適には、時系列光強度データが時間について平滑化され、平滑化された時系列光強度データに於いて粒子の光強度の変化のプロファイル条件に適合する釣鐘状の信号が検出される。しかしながら、かかる態様に於いては、粒子の信号の強度変化が小さいとき、その信号は、フォトンカウンティングデータの段階でノイズに埋もれてしまい、平滑化処理後に於いても、粒子の光強度の変化の釣鐘状のプロファイルの特徴を呈さず、粒子の信号として検出されないこととなる。特に、粒子の存在に対応する信号の大きさは、光検出領域内での粒子の通過経路によって変化し、光検出領域の周辺部を通過する粒子の信号の強度変化は小さい。従って、微弱な多くの光検出領域の周辺部を通過する粒子の信号が見落とされ、或いは、ノイズ信号が微弱な粒子の信号として誤判定されることとなり、このことが、走査分子計数法の感度又は精度の向上の妨げの一つとなっていた。

かくして、本発明の主な課題は、上記の走査分子計数法に於いて、その感度及び／又は精度の更なる向上を可能にする新規な粒子信号の検出手段又は手法を提供することである。

また、本発明のもう一つの課題は、上記の走査分子計数法に於いて、微弱な粒子の信号とノイズ信号とのより精度の良い識別を可能にする新規な粒子信号の検出手段又は手法を提供することである。

上記の課題に関して、本発明の発明者は、走査分子計数法の時系列光強度データに於いて、粒子の信号に対応する部分と、ノイズ信号の部分とでは、光強度の時間変化（光子数列）の発生パターンに違いがあることを見出した。即ち、ノイズ信号の場合には、常にランダムに光子が検出されるのに対して、粒子の信号の場合には、光子が時間的に集中して検出される。本発明に於いては、かかる光強度の時間変化の発生パターンの違いを検出可能な新規なアルゴリズムにて走査分子計数法に於ける粒子の信号の検出精度の向上が図られる。また、光強度の時間変化の発生パターンは、粒子の特性（偏光特性や発光波長特性）によっても変化するので、光強度の時間変化の発生パターンの違いに基づいて、異なる特性の粒子の識別も試みられる。

かくして、本発明によれば、上記の課題は、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析装置であって、試料溶液内に於ける顕微鏡の光学系の光検出領域の位置を移動する光検出領域移動部と、光検出領域からの光を検出する光検出部と、試料溶液内に於いて光検出領域の位置を移動させながら光検出部にて検出された光検出領域からの光の時系列の光強度データを生成し、時系列の光強度データに於いて単一粒子の各々の存在を表す信号を個別に検出する信号処理部とを含み、信号処理部が、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出し、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出することを特徴とする装置によって達成される。

上記本発明の構成に於いて、「試料溶液中にて分散しランダムに運動する単一粒子」とは、試料溶液中に分散又は溶解した原子、分子又はそれらの凝集体などの粒子であって、基板などに固定されず、溶液中を自由にブラウン運動している粒子であれば任意の粒子であってよい。観測対象となる単一粒子は、光を発する粒子（発光粒子）であっても、（検出波長帯域に於いて）光を発さない粒子（非発光粒子）であってもよい。観測対象粒子が発光粒子である場合には、「単一粒子の各々の存在を表す信号」は、発光粒子がその光検出領域内の通過時に発する光に対応した時系列光強度データ上に於ける「光強度値の一時的な増大」となる。発光粒子は、典型的には、蛍光性粒子であるが、りん光、化学発光、生物発光、光散乱等により光を発する粒子であってもよい。また、観測対象粒子は非発光粒子であってもよく、その場合には、光検出領域からの光が有意な背景光を含み、「単一粒子の各々の存在を表す信号」が、背景光からの光強度の一時的な低減となる（反転型走査分子計数法）。なお、共焦点顕微鏡又は多光子顕微鏡の光学系の「光検出領域」とは、それらの顕微鏡に於いて光が検出される微小領域であり、対物レンズから照明光が与えられる場合には、その照明光が集光された領域に相当する（共焦点顕微鏡に於いては、特に対物レンズとピンホールとの位置関係により確定される。発光粒子が照明光なしで発光する場合、例えば、化学発光又は生物発光により発光する粒子の場合には、顕微鏡に於いて照明光は要しない。）。また、典型的には、光検出部は、所定の計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光を検出し、その場合、時系列の光強度データが時系列のフォトンカウントデータとなる。

上記から理解される如く、本発明の装置に於いては、基本的には、特許文献９〜１１に記載の「走査分子計数法」と同様に、試料溶液内に於いて光検出領域の位置を移動しながら、即ち、試料溶液内を光検出領域により走査しながら、逐次的に、光の検出及び時系列の光強度値を表す時系列光強度データの生成が為され、時系列光強度データ上に於いて、単一粒子の存在を表す信号が検出される。かかる単一粒子の存在を表す信号の検出に際して、本発明では、単に光強度値の増減量が所定量を越えたか否かに基づいて単一粒子の信号の検出が実行されるのではなく、時系列光強度データ上の光強度値の時間変化のパターンが、光検出領域内に単一粒子がない場合と光検出領域内に単一粒子が有る場合とのいずれに於いて発生し易いパターンであるかが時系列に評価される。即ち、既に触れた如く、時系列光強度データに於いて、粒子の信号に対応する部分とノイズ信号の部分とでは、光強度の時間変化の発生パターンに違いがあるので、時系列光強度データ上の或る部分に於ける光強度値の時間変化のパターンが、光検出領域内に単一粒子がない場合に発生し易いパターンであるときには、その部分は、ノイズ信号のみの部分であると判定でき、光検出領域内に単一粒子がある場合に発生し易いパターンであるときには、その部分は、粒子の信号に対応する部分であると判定できることとなる。そして、或る光強度値の時間変化のパターンが、光検出領域内に単一粒子がある場合又はない場合のいずれに発生し易いパターンであるかは、その光強度値の時間変化のパターンが、光検出領域内に単一粒子がある場合及びない場合のそれぞれに於いて発生する確率によって判定することができる。

かくして、本発明の装置では、時系列光強度データ上で光検出領域内に粒子が存在した部分、即ち、単一粒子の存在を表す信号を検出するために、上記の如く、「時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率」とが算出される。ここで、「解析窓」とは、時系列光強度データ上に於ける任意の時間幅の領域であり、「解析窓」は、時系列光強度データ上にて逐次的に又は時系列に設定される。また、「第一の発生確率」は、光検出領域内に単一粒子が存在していない場合に、「解析窓」に於いて計測された光強度の時間変化が発生する確率であり、「第二の発生確率」は、光検出領域内に単一粒子が存在している場合に、「解析窓」に於いて計測された光強度の時間変化が発生する確率である。そして、光検出領域内に単一粒子が存在していれば、「第二の発生確率」が「第一の発生確率」に対して相対的に増大することとなるので、「第一の発生確率」と「第二の発生確率」とを参照することによって、時系列の光強度データ上にて光検出領域内に単一粒子の存在した時間領域が決定され、これにより、単一粒子の各々の存在を表す信号の検出が可能となる。

なお、上記の構成に於いて、解析窓の幅は、例えば、単一粒子が光検出領域を通過するのに要する時間以上の任意の幅に設定されてよい。また、解析窓は、時系列光強度データ上に於いて、単位時間毎に或いは所定の時間間隔毎に設定されてもよく（この場合、時系列に設定される解析窓は互いに重複する。）、或いは、時系列光強度データを所定の時間幅にて分割することにより設定されてもよい。「単位時間」とは、光計測に於いて一つの光強度値を与える幅の時間であり、フォトンカウンティングの場合であれば、１つ又はそれ以上の数のビンタイムの幅の時間であってよい。

上記の構成に於いて、第一及び第二の発生確率は、それぞれ、解析窓内に於ける第一及び第二の状態を仮定した場合の光強度値の時間変化の平均のパターンからの実際に計測された光強度値の時間変化のパターンのずれに基づいて決定することができる。平均のパターンからの実測値のパターンのずれが小さいほど、その実測値のパターンの発生する確率が高くなる。従って、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の検出された光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出することが可能である。また、時系列光強度データ上に出現する光強度値は、光検出領域内から放出される光子数又はこれに比例する量であるので、単位時間毎の光強度値は、ポアソン分布に従うと考えられる。そこで、上記の構成に於いて、好適には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。

かくして、上記の装置では、簡単には、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されてよい。また、第一及び第二の発生確率の比又はオッズ比を時系列に算出し、かかる確率の比又はオッズ比の値に基づいて光検出領域内に単一粒子が存在した時間が決定されてもよい。

ところで、上記の走査分子計数法に於いて、光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出して成分の各々の時系列の光強度データを生成することが可能である。その場合、検出される成分のデータに於いて、単一粒子の任意の特性が反映されるように、検出される成分を選択することが可能である。例えば、検出される複数の成分として、偏光方向が互いに異なる成分を選択すれば、複数の成分の時系列光強度データに於いて、粒子の偏光特性が反映される。また、検出される複数の成分として、互いに異なる波長帯域の光成分を選択すれば、複数の成分の時系列光強度データに於いて、粒子の発光波長特性が反映される。そして、複数の成分の時系列光強度データに於ける光強度値の時間変化のパターンに於いても、観測対象粒子の前記の如き特性が反映されることとなる。

そこで、上記の本発明の装置は、更に、光検出部が光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出し、信号処理部が、前記の成分の各々の時系列の光強度データを生成し、更に、信号処理部が、成分の各々の第一及び第二の発生確率を算出できるよう構成されていてよい。かかる構成に於いては、検出成分を適宜選択することにより、観測対象となる単一粒子の有する所定の特性値を検出された光の成分の各々の第二の発生確率に反映させること、即ち、第二の発生確率を単一粒子の有する所定の特性値の関数とすることが可能となる。そして、成分毎の第一及び第二の発生確率の発生確率を参照して、所定の特性値の関数である第二の発生確率が相対的に高い領域に於いて、所定の特性値を有する粒子が存在したと判定できることとなる（所定の特性値を有していない粒子が存在していた場合、所定の特性値の関数である第二の発生確率は低くなる。）。即ち、成分毎の第一及び第二の発生確率の発生確率を用いることにより、単に光検出領域内の粒子の存在の有無だけでなく、時系列の光強度データ上にて光検出領域内に所定の特性値を有する粒子が存在している時間が決定され、所定の特性値を有する粒子の存在を検出することが可能となる。

更に、上記の如き、少なくとも二つの互いに異なる成分毎の第一及び第二の発生確率の発生確率を用いる態様によれば、単一粒子として互いに異なる所定の特性値を有する複数の種類の単一粒子が含まれる場合、それぞれの種類の単一粒子について、時系列の光強度データ上にて粒子の存在した時間が決定できることとなる。即ち、単一粒子の種類毎に、それらの互いに異なる所定の特性値の関数である成分の各々の第二の発生確率を算出し（第一の状態では粒子が存在していないので、第一の発生確率には、粒子の特性値は反映されない。）、複数の成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とを参照したとき、相対的に高い値の第二の発生確率を与える種類の粒子が光検出領域内に存在していると推定することができる。かくして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合、上記の本発明の装置は、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率を算出し、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号を検出するよう構成されていてよい。要すれば、かかる構成により、試料溶液中に複数の種類の単一粒子が含まれている場合に、種類を識別しながら、粒子の検出が可能となる。

なお、上記の第二の発生確率に反映させるべき粒子の特性値としては、既に触れた如く、単一粒子の偏光特性を表す指標値、例えば、蛍光異方性や、単一粒子の発光波長特性を表す指標値、例えば、互いに異なる発光波長帯域の発光強度の比などであってよい。

上記の本発明の装置に於ける試料溶液内での光検出領域の位置の移動速度は、観測対象の単一粒子の特性又は試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよい。光検出領域の移動速度が速くなると、単一粒子が発光粒子である場合には、一つの発光粒子から得られる光量は低減することとなり、また、単一粒子が非発光粒子である場合には、一つの非発光粒子の存在による光強度値の低減量が小さくなる。従って、光検出領域の移動速度は、単一粒子による光強度値の変化が精度よく又は感度よく計測できるように、適宜変更可能となっていることが好ましい。また、試料溶液内での光検出領域の位置の移動速度は、好適には、検出対象となる単一粒子の拡散移動速度（ブラウン運動による粒子の平均の移動速度）よりも高く設定される。上記に説明されている如く、本発明の装置では、検出された光強度値の時間変化のパターンを、光検出領域内に粒子が無い場合と有る場合（第一及び第二の状態）を仮定して算出された時間変化パターンの（第一及び第二の）発生確率を用いて評価するところ、粒子の光検出領域内の通過時に粒子がブラウン運動により移動することも考慮すると、第二の発生確率の算出が複雑となる。そこで、本発明では、粒子の光検出領域内の通過中に於ける粒子のブラウン運動による効果が無視できるように、光検出領域の移動速度が検出対象となる単一粒子の拡散移動速度よりも速く設定されていることが好ましい。なお、拡散移動速度は、単一粒子によって変わるので、上記の如く、その特性（特に、拡散係数）に応じて、本発明の装置は、光検出領域の移動速度が適宜変更可能であるよう構成されていることが好ましい。

試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよい。例えば、レーザー走査型光学顕微鏡に於いて採用されているガルバノミラーを用いるなどして、顕微鏡の光学系の光路を変更して光検出領域の位置が変更されるようになっていてよく、或いは、試料溶液の位置を（例えば、顕微鏡のステージを移動するなどして）動かして、光検出領域の試料溶液内に於ける位置を移動するようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。特に、顕微鏡の光学系の光路を変更して光検出領域の位置が変更される場合、光検出領域の移動は、速やかであり、且つ、試料溶液に於いて機械的振動や流体力学的な作用が実質的に発生しないので、検出対象となる単一粒子が力学的な作用の影響を受けることなく安定した状態にて、光の計測が可能である点で有利である。

上記の本発明の実施の態様の一つに於いては、信号の数を計数して光検出領域に包含された単一粒子の数を計数するようになっていてよい（粒子のカウンティング）。その場合、検出された単一粒子の数と光検出領域の位置の移動量と組み合わせることにより、試料溶液中の同定された単一粒子の数密度又は濃度に関する情報が得られることとなる。具体的には、例えば、複数の試料溶液の数密度若しくは濃度の比、或いは、濃度若しくは数密度の基準となる標準試料溶液に対する相対的な数密度若しくは濃度の比が算出されるか、又は、濃度若しくは数密度の基準となる標準試料溶液に対する相対的な数密度若しくは濃度の比を用いて、絶対的な数密度値又は濃度値が決定されてよい。或いは、任意の手法により、例えば、所定の速度にて光検出領域の位置を移動するなどして、光検出領域の位置の移動軌跡の全体積を特定すれば、単一粒子の数密度又は濃度が具体的に算定できることとなる。

上記の本発明の装置に於いて試料溶液内に於ける光検出領域の位置を移動させながら光検出を行い、単一粒子からの信号を個別に検出する光分析技術であって、検出された光強度値の時間変化のパターンが発生する確率（第一及び第二の発生確率）を参照して、単一粒子の有無を決定する光分析技術の処理は、汎用のコンピュータによっても実現可能である。

従って、本発明のもう一つの態様によれば、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出するための光分析用コンピュータプログラムであって、試料溶液内に於いて顕微鏡の光学系の光検出領域の位置を移動する手順と、試料溶液内に於いて光検出領域の位置を移動させながら光検出領域からの光の強度を測定して光強度データを生成する手順と、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度の時間変化値の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する手順と、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出する手順とをコンピュータに実行させることを特徴とするコンピュータプログラムが提供される。この場合に於いても、本発明の装置と同様に、典型的には、光検出領域からの光を検出して時系列の光強度データを生成する手順に於いて、計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光が検出され、時系列の光強度データが時系列のフォトンカウントデータである。「解析窓」は、本発明の装置の場合と同様に設定されてよい。観測対象となる粒子が、発光粒子である場合には、時系列光強度データ上に於ける「光強度値の一時的な増大」が「単一粒子の各々の存在を表す信号」となる。また、観測対象となる粒子が、（検出波長帯域に於いて）非発光粒子である場合には、光検出領域からの光が有意な背景光を含み、背景光からの光強度の一時的な低減が「単一粒子の各々の存在を表す信号」となる。なお、コンピュータプログラムは、コンピュータ読み取り可能な記憶媒体に記憶されて提供される。コンピュータは、記憶媒体に記憶されているプログラムを読み出して、情報の加工・演算処理を実行することにより、上記の手順を実現する。ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等であってよい。更に、上述のプログラムは、通信回線によってコンピュータへ配信され、この配信を受けたコンピュータがプログラムを実行するようにしても良い。

上記の構成に於いても、典型的には、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出され、より具体的には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。そして、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されるか、第二の発生確率と第一の発生確率とのオッズ比に基づいて、光検出領域内に単一粒子が存在した時間領域が決定されてよい。

また、上記のコンピュータプログラムも、光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、かかる成分の各々の時系列の光強度データが生成され、更に、成分の各々について、第一の発生確率と、観測対象である単一粒子の所定の特性値の関数である第二の発生確率が算出され、成分毎の第一及び第二の発生確率の発生確率に基づいて時系列の光強度データ上にて所定の特性値を有する単一粒子の存在を表す信号が検出されるよう構成されていてよい。そして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合には、上記のコンピュータプログラムに於いて、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率が算出され、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号が検出されてよい。

更に、試料溶液内での光検出領域の位置の移動速度は、単一粒子の特性、試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよく、好適には、試料溶液内での光検出領域の位置の移動速度は、検出対象となる単一粒子の拡散移動速度よりも高く設定される。試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよく、好適には、顕微鏡の光学系の光路を変更して或いは試料溶液の位置を動かして光検出領域の位置が変更されるようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。

更に、上記のコンピュータプログラムに於いても、個別に検出された単一粒子からの信号の数を計数して光検出領域の位置の移動中に検出された単一粒子の数を計数する手順及び／又は検出された単一粒子の数に基づいて、試料溶液中の単一粒子の数密度又は濃度を決定する手順が含まれていてよい。

上記の本発明の装置又はコンピュータプログラムによれば、試料溶液内に於ける光検出領域の位置を移動させながら、個々の発光粒子の光の検出を行う光分析方法であって、検出された光強度値の時間変化のパターンが発生する確率（第一及び第二の発生確率）を参照して、単一粒子の有無を決定する新規な方法が実現される。

かくして、本発明によれば、更に、共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析方法であって、試料溶液内に於いて顕微鏡の光学系の光検出領域の位置を移動する過程と、試料溶液内に於いて光検出領域の位置を移動させながら光検出領域からの光の強度を測定して光強度データを生成する過程と、時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、光検出領域内に単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、光検出領域内に単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する過程と、第一及び第二の発生確率に基づいて時系列の光強度データ上にて単一粒子の各々の存在を表す信号を検出する過程とを含むことを特徴とする方法が提供される。この場合に於いても、本発明の装置と同様に、典型的には、光検出領域からの光を検出して時系列の光強度データを生成する過程に於いて、計測単位時間（ビンタイム）毎に到来する光子数を計数するフォトンカウンティングにより光検出領域からの光が検出され、時系列の光強度データが時系列のフォトンカウントデータである。「解析窓」は、本発明の装置の場合と同様に設定されてよい。観測対象となる粒子が、発光粒子である場合には、時系列光強度データ上に於ける「光強度値の一時的な増大」が「単一粒子の各々の存在を表す信号」となる。また、観測対象となる粒子が、（検出波長帯域に於いて）非発光粒子である場合には、光検出領域からの光が有意な背景光を含み、背景光からの光強度の一時的な低減が「単一粒子の各々の存在を表す信号」となる。

上記の構成に於いても、典型的には、第一及び第二の発生確率は、それぞれ、解析窓内に於ける単位時間毎の光強度値と第一及び第二の状態を仮定した場合の単位時間毎の期待値とに基づいて算出され、より具体的には、単位時間毎の光強度値が単位時間毎の期待値を有するポアソン分布に従うものとして単位時間毎の光強度値の単位時間発生確率が算出され、第一及び第二の発生確率がそれぞれ対応する単位時間発生確率を用いて算出されてよい。そして、第二の発生確率が第一の発生確率より大きい解析窓の時間に光検出領域内に単一粒子が存在したと判定されるか、第二の発生確率と第一の発生確率とのオッズ比に基づいて、光検出領域内に単一粒子が存在した時間領域が決定されてよい。

また、上記の方法も、光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、かかる成分の各々の時系列の光強度データが生成され、更に、成分の各々について、第一の発生確率と、観測対象である単一粒子の所定の特性値の関数である第二の発生確率が算出され、成分毎の第一及び第二の発生確率の発生確率に基づいて時系列の光強度データ上にて所定の特性値を有する単一粒子の存在を表す信号が検出されるよう構成されていてよい。そして、単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含む場合には、上記の方法に於いて、単一粒子の種類毎に、互いに異なる所定の特性値の関数である成分の各々の第二の発生確率が算出され、成分の各々の第一の発生確率と複数の種類の単一粒子の各々の成分の各々の第二の発生確率とに基づいて時系列の光強度データ上にて単一粒子の種類毎に該単一粒子の存在を表す信号が検出されてよい。

更に、試料溶液内での光検出領域の位置の移動速度は、単一粒子の特性、試料溶液中の数密度又は濃度に基づいて適宜変更可能となっていてよく、好適には、試料溶液内での光検出領域の位置の移動速度は、検出対象となる単一粒子の拡散移動速度よりも高く設定される。試料溶液内での光検出領域の位置の移動は、任意の方式で為されてよく、好適には、顕微鏡の光学系の光路を変更して或いは試料溶液の位置を動かして光検出領域の位置が変更されるようになっていてよい。光検出領域の位置の移動軌跡は、任意に設定されてよく、例えば、円形、楕円形、矩形、直線及び曲線のうちから選択可能であってよい。

更に、上記の方法に於いても、個別に検出された単一粒子からの信号の数を計数して光検出領域の位置の移動中に検出された単一粒子の数を計数する過程及び／又は検出された単一粒子の数に基づいて、試料溶液中の単一粒子の数密度又は濃度を決定する過程が含まれていてよい。

上記の本発明の光分析技術は、典型的には、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞などの粒子状の生物学的な対象物の溶液中の状態の分析又は解析の用途に用いられるが、非生物学的な粒子（例えば、原子、分子、ミセル、金属コロイドなど）の溶液中の状態の分析又は解析に用いられてもよく、そのような場合も本発明の範囲に属することは理解されるべきである。

かくして、本発明によれば、走査分子計数法に於いて、単に、時系列光強度データ上の光強度値の増減を参照するのではなく、光強度値の時間変化のパターンが発生しやすい状態を決定して光検出領域内に単一粒子が存在しているか否かが推定される。かかる構成によれば、単一粒子による光強度値の変化が比較的小さく、光強度値変化の絶対値のみでは、粒子による信号かノイズ信号かが判断が難しい場合でも、より精度良く、粒子の信号とノイズ信号との識別が可能となることが期待される。また、粒子の信号とノイズ信号との識別の精度が向上されることにより、走査分子計数法に於いて検出可能な試料溶液中の単一粒子の濃度範囲の低濃度側への拡大が期待される。更に、観測対象となる単一粒子の所定の特性値を考慮して算出された光強度値の時間変化の発生確率を用いる場合には、観測対象となる単一粒子とそれ以外の粒子との識別も可能となり、不純物等が混在する試料溶液に於いても計測の精度又は感度の向上が期待される。

本発明のその他の目的及び利点は、以下の本発明の好ましい実施形態の説明により明らかになるであろう。

図１（Ａ）は、本発明による走査分子計数法を実行する光分析装置の内部構造の模式図である。図１（Ｂ）は、コンフォーカル・ボリューム（共焦点顕微鏡の観察領域）の模式図である。図１（Ｃ）は、ミラー７の向きを変更して試料溶液内に於いて光検出領域の位置を移動する機構の模式図である。図１（Ｄ）は、マイクロプレートの水平方向位置を移動して試料溶液内に於ける光検出領域の位置を移動する機構の模式図である。 図２（Ａ）、（Ｂ）は、それぞれ、本発明が適用される走査分子計数法に於ける発光粒子の光を検出の原理を説明する模式図及び計測される光強度の時間変化の模式図である。図２（Ｃ）、（Ｄ）は、それぞれ、本発明が適用される反転型走査分子計数法に於ける検出光波長帯域で発光しない単一粒子の存在を検出する原理を説明する模式図及び計測される光強度の時間変化の模式図である。 図３（Ａ）は、時系列光強度データの模式的な例を示している。（左）輝度が大きい発光粒子が存在している場合。（中）輝度が小さい発光粒子が存在している場合。（右）発光粒子が存在していない場合。図３（Ｂ）は、本発明に於いて時系列光強度データ上に設定される解析窓を説明する図である。図３（Ｃ）は、試料溶液内の光検出領域の位置を発光粒子の拡散移動速度よりも速い速度にて移動することにより発光粒子が光検出領域を横切る場合の粒子の運動の態様を表すモデル図である。図３（Ｄ）は、光検出領域に於ける発光粒子から放出され検出される光の光検出領域の半径ｒ方向の強度分布の模式図である。 図４（Ａ）は、説明の目的で模式的に表した時系列光強度データであり、上段は、発光粒子が存在している状態（左：第二の状態）及び発光粒子が存在していない状態（右：第一の状態）をそれぞれ仮定した場合の時系列光強度データ上の光強度値の期待値を表しており、中段は、実際に計測された光強度値を表しており、下段は、発光粒子が存在している状態（左）及び発光粒子が存在していない状態（右）をそれぞれ仮定した場合の計測された光強度値の発生確率を示している。図４（Ｂ）は、発光粒子から放出される光の互いに偏光成分のモデル図（上）と、各成分の光強度値の平均的な時間変化（左下）と、合計の光強度値の平均的な時間変化（右下）とを示している。図４（Ｃ）は、互いに異なる発光波長特性を有する発光粒子Ｄ１、Ｄ２の発光波長スペクトルであり、互いに異なる波長帯域の光の成分を別々に検出した場合に於いて、各波長帯域にて検出される各発光粒子からの成分の光強度値を示している。陰影部分の面積が各成分の光強度値に相当する。 本発明に従って実行される走査分子計数法の処理手順をフローチャートの形式で表した図である。 図６は、本発明による走査分子計数法に従って得られた時系列光強度データ（フォトンカウントデータ）の一部（Ａ）と、（Ａ）の時系列光強度データをスムージングして得られた平滑化された時系列光強度データ（Ｂ）と、（Ａ）の時系列光強度データの、本発明の教示に従って粒子が存在する場合と存在しない場合の各々を仮定して算出された発生確率（第二及び第一の発生確率）のオッズ比（Ｃ）を示している。 図７は、図６と同様の図であって、時系列光強度データの別の部分を示している。 図８は、本発明による走査分子計数法に従って発光粒子の検出及び計数を行った場合の試料溶液中の発光粒子（ＡＴＴＯ６４７Ｎ）濃度と検出された粒子数との関係を示した図である。 図９（Ａ）は、本発明による反転型走査分子計数法に従って得られた時系列光強度データ（フォトンカウントデータ）の一部（上段）と、その時系列光強度データの、本発明の教示に従って粒子が存在する場合と存在しない場合の各々を仮定して算出された発生確率（第二及び第一の発生確率）のオッズ比（下段）を示している。図９（Ｂ）は、（Ａ）の一部を拡大した図である。 図１０は、図９（Ｂ）を更に拡大した図である。 図１１（Ａ）は、本発明による反転型走査分子計数法に従って粒子の検出及び計数を行った場合の試料溶液中の粒子濃度と検出された粒子数との関係を示した図であり、図１１（Ｂ）は、（Ａ）と同一の時系列光強度データに於いて、試料溶液中の粒子濃度と平滑化及び釣鐘型関数のフィッティングにより検出された粒子数との関係を示した図である。図中、エラーバーは、標準偏差である。 図１２は、本発明による走査分子計数法に従って光の計測を互いに異なる偏光方向の成分を別々に検出して得られた時系列光強度データに於いて、発光粒子の蛍光偏光異方性を考慮して算出された発生確率を用いて、発光粒子の検出及び発光粒子の種類の識別を行ったときの結果を示している。図中、左は、プラスミドを含む溶液に於いて検出された発光粒子の種類の判別結果であり、右は、蛍光色素（ＴＡＭＲＡ）を含む溶液に於いて検出された発光粒子の種類の判別結果である。

１…光分析装置（共焦点顕微鏡）
２…光源
３…シングルモードオプティカルファイバー
４…コリメータレンズ
５…ダイクロイックミラー
６、７、１１…反射ミラー
８…対物レンズ
９…マイクロプレート
１０…ウェル（試料溶液容器）
１２…コンデンサーレンズ
１３…ピンホール
１４…バリアフィルター
１４ａ…ダイクロイックミラー又はビームスプリッター
１５…マルチモードオプティカルファイバー
１６…光検出器
１７…ミラー偏向器
１７ａ…ステージ位置変更装置
１８…コンピュータ

以下、本発明の好ましい実施形態について詳細に説明する。

光分析装置の構成
本発明による光分析技術を実現する光分析装置は、基本的な構成に於いて、図１（Ａ）に模式的に例示されている如き、ＦＣＳ、ＦＩＤＡ等が実行可能な共焦点顕微鏡の光学系と光検出器とを組み合わせてなる装置であってよい。同図を参照して、光分析装置１は、光学系２〜１７と、光学系の各部の作動を制御すると共にデータを取得し解析するためのコンピュータ１８とから構成される。光分析装置１の光学系は、通常の共焦点顕微鏡の光学系と同様であってよく、そこに於いて、光源２から放射されシングルモードファイバー３内を伝播したレーザー光（Ｅｘ）が、ファイバーの出射端に於いて固有のＮＡにて決まった角度にて発散する光となって放射され、コリメーター４によって平行光となり、ダイクロイックミラー５、反射ミラー６、７にて反射され、対物レンズ８へ入射される。対物レンズ８の上方には、典型的には、１〜数十μＬの試料溶液が分注される試料容器又はウェル１０が配列されたマイクロプレート９が配置されており、対物レンズ８から出射したレーザー光は、試料容器又はウェル１０内の試料溶液中で焦点を結び、光強度の強い領域（励起領域）が形成される。

観測対象となる単一粒子が発光粒子である場合、試料溶液中には、発光粒子、典型的には、蛍光性粒子又は蛍光色素等の発光標識が付加された粒子が分散又は溶解され、かかる発光粒子が励起領域に進入すると、その間、発光粒子が励起され光が放出される。一方、観測対象となる単一粒子が非発光粒子である場合、試料溶液中には、典型的には、検出波長帯域にて発光しない粒子と、背景光を生ずる任意の発光物質が分散又は溶解され、検出波長帯域にて発光しない粒子が励起領域に進入していないときには、発光物質が励起されて実質的に一定の光が放出されて、背景光となり、検出波長帯域にて発光しない粒子が励起領域に進入すると、背景光が低減することとなる。

かくして、励起領域から放出された光（Ｅｍ）は、対物レンズ８、ダイクロイックミラー５を通過し、ミラー１１にて反射してコンデンサーレンズ１２にて集光され、ピンホール１３を通過し、バリアフィルター１４を透過して（ここで、特定の波長帯域の光成分のみが選択される。）、マルチモードファイバー１５に導入されて、光検出器１６に到達し、時系列の電気信号に変換された後、コンピュータ１８へ入力され、後に説明される態様にて光分析のための処理が為される。なお、当業者に於いて知られている如く、上記の構成に於いて、ピンホール１３は、対物レンズ８の焦点位置と共役の位置に配置されており、これにより、図１（Ｂ）に模式的に示されている如きレーザー光の焦点領域、即ち、光検出領域内から発せられた光のみがピンホール１３を通過し、光検出領域以外からの光は遮断される。図１（Ｂ）に例示されたレーザー光の焦点領域は、通常、１〜１０ｆＬ程度の実効体積を有する本光分析装置に於ける光検出領域であり（典型的には、光強度が領域の中心を頂点とするガウス様分布となる。実効体積は、光強度が中心光強度の１／ｅ^２となる面を境界とする略楕円球体の体積である。）、コンフォーカル・ボリュームと称される。また、本発明では、１つの発光粒子からの光、例えば、一つの蛍光色素分子からの微弱光又は非発光粒子の存在による背景光の低減が検出されるので、光検出器１６としては、好適には、フォトンカウンティングに使用可能な超高感度の光検出器が用いられる。光の検出がフォトンカウンティングによる場合、光強度の測定は、所定時間に亘って、逐次的に、計測単位時間（ＢＩＮ ＴＩＭＥ）毎に、光検出器に到来するフォトンの数を計測する態様にて実行される。従って、この場合、時系列の光強度のデータは、時系列のフォトンカウントデータである。また、顕微鏡のステージ（図示せず）には、観察するべきウェル１０を変更するべく、マイクロプレート９の水平方向位置を移動するためのステージ位置変更装置１７ａが設けられていてよい。ステージ位置変更装置１７ａの作動は、コンピュータ１８により制御されてよい。かかる構成により、検体が複数在る場合にも、迅速な計測が達成可能となる。

更に、上記の光分析装置の光学系に於いては、光検出領域により試料溶液内を走査する、即ち、試料溶液内に於いて焦点領域、即ち、光検出領域の位置を移動するための機構が設けられる。かかる光検出領域の位置を移動するための機構としては、例えば、図１（Ｃ）に模式的に例示されている如く、反射ミラー７の向きを変更するミラー偏向器１７が採用されてよい（光検出領域の絶対的な位置を移動する方式）。かかるミラー偏向器１７は、通常のレーザー走査型顕微鏡に装備されているガルバノミラー装置と同様であってよい。或いは、別の態様として、図１（Ｄ）に例示されている如く、試料溶液が注入されている容器１０（マイクロプレート９）の水平方向の位置を移動し、試料溶液内に於ける光検出領域の相対的な位置を移動するべくステージ位置変更装置１７ａが作動されてもよい（試料溶液の絶対的な位置を移動する方式）。いずれの方式による場合も、所望の光検出領域の位置の移動パターンを達成するべく、ミラー偏向器１７又はステージ位置変更装置１７ａは、コンピュータ１８の制御の下、光検出器１６による光検出と協調して駆動される。光検出領域の位置の移動軌跡は、円形、楕円形、矩形、直線、曲線又はこれらの組み合わせから任意に選択されてよい（コンピュータ１８に於けるプログラムに於いて、種々の移動パターンが選択できるようになっていてよい。）。なお、図示していないが、対物レンズ８又はステージを上下に移動することにより、光検出領域の位置が上下方向に移動されるようになっていてもよい。

観測対象物となる発光粒子又は背景光を与える物質が多光子吸収により発光する場合には、上記の光学系は、多光子顕微鏡として使用される。その場合には、励起光の焦点領域（光検出領域）のみで光の放出があるので、ピンホール１３は、除去されてよい。また、発光粒子又は背景光を与える物質が化学発光や生物発光現象により励起光によらず発光する場合には、励起光を生成するための光学系２〜５が省略されてよい。発光粒子又は背景光を与える物質がりん光又は散乱により発光する場合には、上記の共焦点顕微鏡の光学系がそのまま用いられる。更に、光分析装置１に於いては、図示の如く、複数の励起光源２が設けられていてよく、発光粒子又は背景光を与える物質の励起波長によって適宜、励起光の波長が選択できるようになっていてよい。

また更に、光検出器１６も複数個備えられ、各々の光検出器１６が光検出領域からの光のうちの互いに異なる成分を別々に検出するよう構成されていてよい。後により詳細に述べる如く、検出する成分を適宜選択することにより、特定の発光特性を有する単一粒子を選択的に検出することが可能となる。そのような光検出領域からの光のうちの互いに異なる成分を検出する場合には、ピンホール１３の通過後の検出光光路に於いては、光路を任意の態様にて分割するための機構が設けられる。例えば、光検出領域からの光を互いに異なる偏光成分に分けて検出する場合には、検出光光路の符号１４ａが付された部位に偏光ビームスプリッタ１４ａが挿入される。この場合、励起光光路には、ポーラライザ（図示せず）が挿入される。また、検出光光路の部位１４ａに、検出光光路に特定の波長帯域の光を反射し、別の波長帯域を透過するダイクロイックミラー１４ａが挿入されることにより、互いに異なる波長帯域の光成分が別々に検出可能となる。

コンピュータ１８は、ＣＰＵおよびメモリを備え、ＣＰＵが各種演算処理を実行することにより、本発明の手順を実行する。なお、各手順は、ハードウェアにより構成するようにしてもよい。本実施形態で説明される処理の全て或いは一部は、それらの処理を実現するプログラムを記憶したコンピュータ読み取り可能な記憶媒体を用いて、コンピュータ１８により実行されてよい。即ち、コンピュータ１８は、記憶媒体に記憶されているプログラムを読み出して、情報の加工・演算処理を実行することにより、本発明の処理手順を実現するようになっていてよい。ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等であってよく、或いは、上記のプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータがプログラムを実行するようにしても良い。

本発明の光分析技術の原理
「発明の概要」の欄に記載されている如く、本発明の光分析技術に於いては、端的に述べれば、走査分子計数法又は反転型走査分子計数法に於いて、光検出領域の位置を移動しながら計測して得られた時系列光強度データから観測対象となる単一粒子の信号を検出する際に、単に、時系列光強度データ上の光強度値の増減を参照するのではなく、計測された光強度値の時間変化のパターンが、光検出領域内に粒子が存在する場合と存在しない場合とのいずれに於いて、より高い確率で発生するパターンであるかが推定される。そして、その推定結果に基づいて、時系列光強度データ上で観測対象となる単一粒子の信号の存在の有無及び／又はその数が検出される。以下、本発明の走査分子計数法及び観測対象となる単一粒子の信号の検出の原理について説明する。

１.走査分子計数法の原理
走査分子計数法に於いて実行される基本的な処理に於いては、特許文献９〜１１に記載されている如く、端的に述べれば、光検出領域の位置を移動するための機構（ミラー偏向器１７）を駆動して光路を変更し、或いは、試料溶液が注入されている容器１０（マイクロプレート９）の水平方向の位置を移動して、図２（Ａ）にて模式的に描かれているように、試料溶液内に於いて光検出領域ＣＶの位置を移動しながら、即ち、光検出領域ＣＶにより試料溶液内を走査しながら、光検出が実行される。そうすると、観測対象粒子が発光粒子である場合、光検出領域ＣＶが移動する間（図中、時間ｔｏ〜ｔ２）に於いて１つの発光粒子の存在する領域を通過する際（ｔ１）には、発光粒子から光が放出され、図２（Ｂ）に描かれている如き時系列の光強度データ上に有意な光強度（Ｅｍ）のパルス状の信号が出現することとなる。そこで、上記の光検出領域ＣＶの位置の移動と光検出を実行し、その間に出現する図２（Ｂ）に例示されている如きパルス状の信号（有意な光強度）を一つずつ検出することによって、発光粒子が個別に検出され、その数をカウントすることにより、計測された領域内に存在する発光粒子の数、或いは、濃度若しくは数密度に関する情報が取得できることとなる。

また、観測対象となる粒子が非発光粒子である場合、上記の走査分子計数法による光計測に於いて、光検出領域から背景光を放出させ（或いは、光検出領域を照明光にて照明し）、観測対象となる粒子が光検出領域に進入した際の検出される背景光の低下を捉えることにより、非発光粒子の存在を検出することが可能となる（反転型走査分子計数法）。この場合、より具体的には、図２（Ｃ）にて模式的に描かれているように、試料溶液内に於いて光検出領域ＣＶの位置を移動しながら、光検出が実行される。ここで、試料溶液中には発光物質が分散させておくと、光検出領域ＣＶ内には、多数の発光物質が存在し、光検出領域ＣＶの移動中（図中、時間ｔｏ〜ｔ２）、それらの発光物質からの光が略一様に検出される。しかしながら、光検出領域ＣＶが移動する間に単一の非発光粒子の存在する領域を通過する際（ｔ１）には、発光物質の存在領域の体積が低減するので、発光物質の放出する光の総量が低減し、図２（Ｄ）に描かれている如く、時系列の光強度データ上に、釣鐘型のパルス状の有意な光強度（Ｅｍ）の低下が出現することとなる。かくして、上記の光検出領域ＣＶの位置の移動と光検出を実行し、その間に出現する図２（Ｄ）に例示されている如きパルス状の有意な光強度の低下、即ち、非発光粒子の存在を表す信号を一つずつ検出することによって、非発光粒子が個別に検出され、その数をカウントすることにより、計測された領域内に存在する単一粒子の数、或いは、濃度若しくは数密度に関する情報が取得できることとなる。

ところで、上記の如き走査分子計数法に於いて、実際に計測された時系列光強度データ上には、発光粒子の信号の他に、熱ノイズ、迷光、水のラマン散乱等による光強度値の増大（ノイズ信号）が存在する。また、反転型走査分子計数法の場合には、背景光の強度のゆらぎによる光強度値の低減（この場合も、以下、ノイズ信号と称する。）が存在する。そこで、時系列光強度データ上に於いて単一粒子の信号の検出を行う際には、まず、パルス状の信号を検出し、しかる後に、検出されたパルス状の信号の強度、時間幅及び形状が検査され、単一粒子が光検出領域を通過した場合に得られるべき信号の強度、時間幅及び形状の条件に適合している信号のみが、単一粒子の信号として判定され、それ以外の信号はノイズ信号として判定される。

しかしながら、単一粒子の信号の強度値が弱くなると、信号の強度、時間幅及び形状に基づく単一粒子の信号とノイズ信号との識別は、困難となる。特に、観測対象粒子が発光粒子である場合、発光粒子から放出される光子数が微量であり、且つ、光子は確率的に放出されるので、光強度値のプロファイルは、実際には、図２（Ｂ）に描かれている如き滑らかな釣鐘状にはならず、図３（Ａ）左にて例示されている如く、離散的になる。従って、発光粒子から放出される光子数が少なくなると、その場合の光強度データ上の信号の強度、時間幅及び形状（図３（Ａ）中）は、粒子が存在せずノイズ信号のみ存在する光強度データ（図３（Ａ）右）の場合と区別がつき難くなる。また、光検出領域に於ける発光粒子から放出され検出される光の強度値は、発光粒子の位置が光検出領域の略中心から離れるほど低下し（図３（Ｄ）参照）、且つ、光検出領域の略中心から離れるほど発光粒子の数は多くなるので、ノイズ信号との識別が困難な粒子の信号の数が大きくなる。同様の状況は、観測対象粒子が非発光粒子である反転型走査分子計数法に於ける背景光の光子数の低下の場合にも生ずる。

走査分子計数法及び反転型走査分子計数法に於いて、その精度及び／又は感度は、光検出領域の或る移動距離当たりの粒子信号の検出精度と検出数が高いほど、向上される。しかしながら、上記の如く、粒子信号の検出数を多くすべく、信号の強度の判定基準値を低減すると、ノイズ信号を誤って粒子信号と判定する可能性が増大し、粒子信号の検出精度を向上すべく、信号の強度の判定基準値を増大すると、多くの輝度の低い粒子信号が検出されないこととなる。そこで、本発明に於いては、信号の強度、時間幅及び形状によってはノイズ信号との識別が困難な単一粒子の信号の検出がより精度良く達成可能とする新規なアルゴリズムが提案される。

２．単一粒子の信号の検出
時系列光強度データがフォトンカウントデータである場合に於いて、光強度値は、ビンタイム当たりの検出光子数となる。従って、図３（Ａ）に模式的に例示されているように、光強度値は、時間軸方向に離散的に分布する。その場合、輝度が大きい粒子信号は、ノイズ信号の影響が少なく、略釣鐘状のプロファイル（左図）を有するが、輝度の小さい粒子信号は、強度値がノイズ信号と略同様となり、更に、ノイズ信号が重畳するので（中図）、略釣鐘状のプロファイルを抽出することが困難となり、粒子のないノイズ信号のみが存在する時間領域（右図）との識別が困難となる。しかしながら、輝度の小さい粒子信号の存在する時間領域とノイズ信号のみが存在する時間領域とでは、光子が検出される事象（光子数が１以上となる事象）の発生頻度とパターンに違いを有する。即ち、図からも理解される如く、ノイズ信号の場合は、常にランダムに光子検出の事象が生ずるのに対し、粒子信号の場合は、光子検出の事象が時間的に集中し、特に、その中央付近の強度値が高くなる傾向を有する。同様の現象は、非発光粒子の信号と背景光のゆらぎとの間にも観察される。

かくして、本発明では、かかる光子検出の事象の発生頻度とパターンの違いに着目して、粒子信号を選択的に検出することが試みられる。そのために、本発明による新規な単一粒子信号の検出のためのアルゴリズムに於いては、時系列光強度データ上の所定の時間幅毎の光強度値の時間変化（光子数列）について、光検出領域内に粒子が存在していたと仮定した場合に前記の光強度値の時間変化が発生する確率（発生確率（粒子存在時）、第二の発生確率）と、光検出領域内に粒子が存在していなかったと仮定した場合に前記の光強度値の時間変化が発生する確率（発生確率（粒子不在時）、第一の発生確率）とが算出される。そして、高い発生確率を与える状態が実際の状態であると推定される。

より具体的には、まず、図３（Ｂ）を参照して、時系列光強度データに於いて、任意の時間幅の区間（以下、解析窓と称する。）が設定される。この解析窓は、単位時間（通常は、ビンタイムであってよい。）ｔｉ（ｉ＝１，２，… 以下、同様）毎に、検出された光子数Ｃｉを有する。ところで、各単位時間内に於いて生ずる光子検出の事象の数は、その単位時間に於ける期待値を有するポアソン分布に従うと考えられるので、任意の単位時間ｔｉに於いて、検出光子数Ｃｉとなる確率（単位時間発生確率）は、

により与えられる。ここで、Ｅｉは、単位時間ｔｉに於ける光子数の期待値である。そして、解析窓にｎ＋１個の単位時間が含まれるとき、解析窓内に於ける検出光子数列Ｃｉが発生する確率Ｐ（発生確率）は、

により与えられる。

上記の各単位時間ｔｉに於ける光子検出事象の発生数の期待値Ｅｉは、光計測中の解析窓に対応する時間領域に於ける光検出領域内の粒子の有無によって決定される。かくして、光検出領域内に粒子が存在していない場合、光子検出の事象は、常にランダムに生ずるので、各計測単位時間ｔｉに於ける期待値Ｅｎｉは、
Ｅｎｉ＝Ｂｇ …（３）
と設定されてよい。ここで、Ｂｇは、観測対象粒子が発光粒子である場合は、ノイズ信号の時間平均値であり、観測対象粒子が非発光粒子である場合は、背景光の時間平均値である。従って、式（３）の値を式（１）に代入して、各計測単位時間ｔｉに於ける単位時間発生確率Ｐｎｉが時系列に算定され、更に、式（２）により、光検出領域内に粒子が存在していない状態を仮定した場合に実際の検出光子数列が発生する確率Ｐｎ（第一の発生確率）が算出される。

一方、光検出領域内に粒子が存在していた場合、光検出領域ＣＶの位置は移動しているので、粒子は、図３（Ｃ）に模式的に描かれている如く、光検出領域ＣＶ内を通過することとなる。かかる過程に於いて、光検出領域に於ける粒子から放出され検出される光の強度値又は光検出領域内に粒子が存在することによって低減する背景光の低減量は、図３（Ｄ）に示されている如く、粒子の位置が光検出領域の略中心から離れるほど低下する。従って、光検出領域内に粒子が存在していた場合の各単位時間ｔｉに於ける期待値Ｅｐｉも、時間を変数とする釣鐘状の関数として表される。ここで、その釣鐘状の関数がガウス関数で近似されるとすると、期待値Ｅｐｉは、

により与えられる。ここで、ガウス関数は、解析窓内の任意の時間ｔｃ（例えば、解析窓の中心）にてピーク強度Ｑを有すると仮定されている。また、式（４）のガウス関数の半値全幅は、移動速度ｖの光検出領域が、図３（Ｄ）に例示されている如き、光検出領域に於ける粒子から放出され検出される光の強度値又は光検出領域内に粒子が存在することによって低減する背景光の低減量の半径ｒ方向の分布の半値全幅ｄを通過する時間ｄ／ｖに等しく、この条件から、ｗは、

により与えられる。なお、図３（Ｄ）の半値全幅ｄは、光学系から決定可能である。

式（４）中のピーク強度Ｑは、観測対象粒子が発光粒子の場合、解析窓内の光子数の総計が式（４）の期待値の総計に一致すると仮定すると、

により与えられる。また、観測対象粒子が非発光粒子の場合には、光子数の低減量の絶対値の期待値が式（４）に従うと設定して、ピーク強度Ｑは、

により与えられる。かくして、式（４）の値を式（１）に代入して、各計測単位時間ｔｉに於ける単位時間発生確率Ｐｐｉが時系列に算定され、更に、式（２）により、光検出領域内に粒子が存在していた状態を仮定した場合に実際の検出光子数列が発生する確率Ｐｐ（第二の発生確率）が算出される。そして、粒子が存在していた状態を仮定した場合の検出光子数列の発生確率Ｐｐが粒子が存在していなかった状態を仮定した場合の検出光子数列の発生確率Ｐｎを所定の程度上回るとき、その解析窓に於いて、粒子の信号が存在していると判定される。

図４（Ａ）は、上記の一連の処理の例を示している。同図を参照して、図の例に於いて、図中中段の如く、時間軸（矢印）上の単位時間（長方形）の各々に順に光子数が［０００１１１００］となる光子数列が検出されているとき、背景光の平均が０．１であるとすると、上段右図の如く、粒子が存在していないと仮定した場合の各単位時間の期待値は全域に亘って０．１となる。一方、粒子が存在していると仮定した場合の各単位時間の期待値は、検出光子数を用いた式（６）と、更にその結果を用いた式（４）とから、上段左図の如き釣鐘状のプロファイルを有する値となる。かくして、単位時間毎の検出光子数と期待値を用いた式（１）により、単位時間毎に、検出光子数が発生する確率（ビンタイム毎の発生確率）が、粒子が存在していないと仮定した場合（下段右）と粒子が存在していたと仮定した場合（下段左）とのそれぞれについて算出される。そして、かかる検出光子数が発生する確率をそれぞれ式（２）に従って総乗して、検出された光子数列について、粒子が存在していないと仮定した場合の光子数列の発生確率Ｐｎと粒子が存在していたと仮定した場合光子数列の発生確率Ｐｐとが算出される。図示の例では、発生確率Ｐｐ＞発生確率Ｐｎとなるので、検出光子数列が粒子の信号として判定される。

なお、実施の態様に於いては、解析窓内に粒子の信号が存在するか否かの判定は、発生確率Ｐｐと発生確率Ｐｎとのオッズ比ＯＲ
ＯＲ＝Ｐｐ（１−Ｐｎ）／（１−Ｐｐ）Ｐｎ …（８）
を算出し、その大きさが所定値を越えたときに、解析窓内に於ける粒子の信号の存在が判定されてよい。

上記の粒子の信号の検出処理過程に於いて、解析窓は、好適には、単一粒子が光検出領域内を通過するのに要する時間以上に設定される。半径ｒの光検出領域が速度ｖにて移動しているとすれば、解析窓の時間幅は、
２ｒ／ｖ …（９）
よりも長く設定される。また、解析窓は、好適には、時系列光強度データ上にて単位時間毎に順に設定されてよい。かかる設定によれば、発生確率Ｐｐ、発生確率Ｐｎ及び／又はオッズ比ＯＲが時系列光強度データに沿って算出されることとなる。しかしながら、その場合、演算量が多くなるので、数個の単位時間毎に解析窓が設定されるようになっていてもよい。更に、時系列光強度データ上を解析窓の時間幅にて分割して解析窓を設定してもよい。この場合、解析窓は、互いに重複せずに設定されることとなる。

解析窓が単位時間毎に設定される場合には、一つの粒子信号が存在するとき、粒子の信号の存在の判定が連続する解析窓に於いて継続することとなる。即ち、一つの粒子の信号は、粒子の信号の存在の判定が継続する一つの区間に対応することとなる。従って、粒子の信号の存在の判定が継続する区間の数を計数することにより粒子の信号の計数が可能となる。また、上記の粒子の信号の検出処理過程に於いて、ビンタイムは、単一粒子が光検出領域内を通過するのに要する時間（式９）以下に設定される。これは、単一粒子が通過する際の信号を複数のビンタイムに亘って捉えて、単一粒子が通過する間の光強度値の時間変化のパターンが検出できるようにするためである。（もしビンタイムが式９の時間よりも長いと単一粒子が通過する間の光強度値の時間変化のパターンが捕捉できない。）

３．複数の光成分の計測を用いた単一粒子信号の検出
上記の発生確率Ｐｐ、Ｐｎを用いた単一粒子信号の検出手法は、互いに異なる複数の光成分を別々に計測して成分毎に時系列光強度データを生成する場合に拡張することが可能である。その場合、観測対象となる単一粒子の発光特性が、成分毎の時系列光強度データに反映されるように、検出成分を選択することにより、かかる観測対象となる単一粒子の発光特性の反映された信号を選択的に検出することが可能となる。以下、単一粒子の発光特性の反映された信号を選択的に検出する例のいくつかについて説明する。

（ｉ）特定の偏光特性を有する単一粒子の信号の検出
一定の方向に偏光した光を励起光とし、検出光として、発光粒子の蛍光のｐ波とｓ波を別々に検出する場合、特定の偏光特性を有する粒子の信号を選択的に検出することが可能となる（図４（Ｂ）上参照）。まず、或る発光粒子の蛍光異方性Ｒは、蛍光のｐ波とｓ波の光強度Ｃｐ、Ｃｓを用いて下記の如く与えられる。
Ｒ＝（Ｃｐ−Ｃｓ）／（Ｃｐ＋Ｃｓ） …（１０）
従って、全蛍光強度（Ｃｐ＋Ｃｓ）に対する蛍光のｐ波とｓ波の光強度の割合は、それぞれ
Ｃｐ／（Ｃｐ＋Ｃｓ）＝１／２＋Ｒ／２ …（１１ａ）
Ｃｓ／（Ｃｐ＋Ｃｓ）＝１／２−Ｒ／２ …（１１ｂ）
により与えられる。ところで、図４（Ｂ）下段に模式的に描かれている如く、蛍光のｓ波とｐ波の光強度の和について、粒子が存在している場合の或る解析窓内に於ける光強度の各単位時間の期待値Ｅｐｉは、式（４）と同様であるから、ピーク強度Ｑは、

により与えられる。ここで、Ｃｐｉ、Ｃｓｉ、Ｂｇｐ、Ｂｇｓは、それぞれ、単位時間毎のｐ波とｓ波の光子数、ｓ波とｐ波の背景光の強度である。また、蛍光のｐ波とｓ波のそれぞれの各単位時間の期待値Ｅｐｐｉ、Ｅｐｓｉは、
Ｅｐｐｉ＝（１／２＋Ｒ／２）Ｅｐｉ＋Ｂｇｐ …（１３ａ）
Ｅｐｓｉ＝（１／２−Ｒ／２）Ｅｐｉ＋Ｂｇｓ …（１３ｂ）
により与えられる。ここで、Ｅｐｉは、式（４）に於いて、Ｑを式（１２）により与え、ＢｇをＢｇｐ＋Ｂｇｓにより与えて得られた値である。かくして、粒子が存在していると仮定した場合の蛍光のｐ波とｓ波に於ける単位時間毎の検出光子数の発生確率Ｐｐｐｉ、Ｐｐｓｉは、式（１）と同様に

により与えられる。従って、蛍光異方性Ｒを有する発光粒子が存在していたと仮定した場合の解析窓内の検出光子数列の発生確率Ｐｐは、発生確率Ｐｐｐｉ、Ｐｐｓｉの総乗となるので、

により与えられる。一方、粒子が存在していないと仮定した場合、解析窓内の検出光子数列の発生確率Ｐｎは、

となる。ここで、Ｐｎｐｉ、Ｐｎｓｉは、各成分の単位時間毎の発生確率であり、それぞれ、背景光Ｂｇｐ、Ｂｇｓを期待値として用いて、式（１）と同様に

となる。

かくして、発光粒子の蛍光のｐ波とｓ波を別々に検出する場合、ｐ波及びｓ波の時系列光強度データのそれぞれに於いて、式（１４）及び（１７）による単位時間毎の発生確率を算出した後、式（１５）、（１６）による解析窓に於ける発生確率Ｐｐ、Ｐｎを算出して、これらの発生確率を比較することにより、蛍光異方性Ｒを有する粒子の信号の有無が検出できることとなる。なお、この場合、蛍光異方性Ｒは、既知数である。従って、蛍光異方性Ｒの値は、任意の手法にて実験的又は理論的に決定した値が用いられてよい。理論的に蛍光異方性Ｒの値を算出する場合、例えば、非特許文献５に記載されている如く、粒子の分子量Ｍを用いて、下記の如く推定することができる。
Ｒ＝Ｒｏ／（１＋τθ） …（１８ａ）
ここで、Ｒｏは、粒子が運動していないときの異方性の値（＝０．３）であり、τは、回転相関時間である。θは、下記の式にて与えられる。
θ＝ηＭ（Ｖ＋ｈ）／ＲoＴ …（１８ｂ）
ここで、η、Ｖ＋ｈ、Ｒo、Ｔは、それぞれ、粘度、水和体積、気体定数、温度である。

また、上記の如く、特定の蛍光異方性を有する単一粒子の信号を検出する態様によれば、試料溶液内に複数の互いに蛍光異方性の異なる複数種類の発光粒子が混在している場合に、種類を識別しつつ粒子の信号の検出が可能となる。即ち、粒子の種類の識別を行う場合には、各解析窓の光子数列について、上記の如き粒子が存在していると仮定した場合の発生確率が、蛍光異方性の値毎に算出される。そして、互いに異なる蛍光異方性の値を用いて得られた発生確率のうち、最も高い発生確率を与える蛍光異方性を有する粒子が、解析されている光子数列上に存在していると推定できることとなる。

（ii）特定の発光波長特性を有する単一粒子の信号の検出
光検出領域からの光を検出する際に、互いに異なる波長帯域の光成分を別々に計測する場合、各波長帯域の成分の光強度値の大きさは、発光粒子の発光波長特性に依存して変化することとなる。例えば、図４（Ｃ）に示されている如く、互いに異なる発光波長スペクトルを有する発光粒子Ｄ１、Ｄ２について、図中、点線枠にて示された波長帯域の光成分をそれぞれ別々にＣｈ１、Ｃｈ２として検出する場合、Ｃｈ１及びＣｈ２の各々に於いて検出される各発光粒子の光は、各発光粒子の発光波長スペクトルのうち、Ｃｈ１及びＣｈ２の検出波長帯域と重複する部分（陰影のつけられた部分）となり、それぞれの光量は、陰影のつけられた部分の面積となる。従って、図から理解される如く、Ｃｈ１及びＣｈ２に於ける光量の比は、発光粒子の発光波長スペクトルのプロファイルに依存して変化することとなる。

上記の如き発光粒子の発光波長スペクトルに依存した互いに異なる検出波長帯域にて検出された光強度（光子数）の比は、走査分子計数法による時系列光強度データ上の光子数に反映される。そこで、上記の光子数列の発生確率の算出に於いて、かかる光子数の比を所定の特性値として考慮することにより、特定の発光波長特性を有する発光粒子の信号を選択的に検出することが可能となる。

具体的には、まず、或る発光粒子について、Ｃｈ１及びＣｈ２の検出波長帯域にて検出される光の成分の強度比α、βは、それぞれ、
α＝Ｓα／（Ｓα＋Ｓβ）；β＝Ｓβ／（Ｓα＋Ｓβ）…（１９）
により与えられる。Ｓα、Ｓβは、それぞれ、Ｃｈ１及びＣｈ２の検出波長帯域に於ける発光粒子の発光波長スペクトル強度の積分値（面積）である。

Ｃｈ１及びＣｈ２の検出光子数の総和について、粒子が存在している場合の或る解析窓内に於ける光強度の各単位時間の期待値Ｅｐｉは、式（４）と同様であり、ピーク強度Ｑは、

により与えられる。ここで、Ｃαｉ、Ｃβｉ、Ｂｇαｐ、Ｂｇβは、それぞれ、単位時間毎のＣｈ１とＣｈ２の光子数、Ｃｈ１とＣｈ２の背景光の強度である。従って、Ｃｈ１とＣｈ２の光子数のそれぞれの期待値Ｅｐα、Ｅｐβは、

により与えられ、粒子が存在していると仮定した場合のＣｈ１及びＣｈ２に於ける単位時間毎の検出光子数の発生確率Ｐｐαｉ、Ｐｐβｉは、式（１）と同様に

により与えられる。従って、Ｃｈ１とＣｈ２の光強度比がα：βとなる発光波長特性を有する発光粒子が存在していたと仮定した場合の解析窓内の検出光子数列の発生確率Ｐｐは、発生確率Ｐｐαｉ、Ｐｐβｉの総乗となるので、

により与えられる。一方、粒子が存在していないと仮定した場合、解析窓内の検出光子数列の発生確率Ｐｎは、

となる。ここで、Ｐｎαｉ、Ｐｎβｉは、各成分の単位時間毎の発生確率であり、それぞれ、背景光Ｂｇα、Ｂｇβを期待値として用いて、式（１）と同様に

となる。

かくして、発光粒子の互いに異なる発光波長帯域の成分を別々に検出する場合、各検出波長帯域の時系列光強度データのそれぞれに於いて、式（２２）及び（２５）による単位時間毎の発生確率を算出した後、式（２３）、（２４）による解析窓に於ける発生確率Ｐｐ、Ｐｎを算出して比較することにより、特定の発光波長特性を有する粒子の信号の有無が検出できることとなる。なお、この場合、Ｃｈ１とＣｈ２の光強度比α：βは、既知数である。かかる強度比は、任意の手法にて実験的又は理論的に決定した値が用いられてよい。また、上記の手法によれば、特定の発光波長特性を有する粒子の信号を選択的に検出できることとなるので、試料溶液内に複数の互いに発光波長特性の異なる複数種類の発光粒子が混在している場合に、種類を識別しつつ粒子の信号の検出が可能となる。即ち、粒子の種類の識別を行う場合には、各解析窓の光子数列について、上記の如き粒子が存在していると仮定した場合の発生確率が、Ｃｈ１とＣｈ２の光強度比（α：β）毎に算出される。そして、互いに異なる光強度比を用いて得られた発生確率のうち、最も高い発生確率を与える光強度比を有する粒子が、解析されている光子数列上に存在していると推定できることとなる。

走査分子計数法の処理操作過程
図１（Ａ）に例示の光分析装置１を用いた本発明に従った走査分子計数法の実施形態に於いては、具体的には、（１）単一粒子を含む試料溶液の調製、（２）試料溶液の光強度の測定処理、（３）単一粒子信号の検出処理が実行される。図５は、フローチャートの形式にて表した本実施形態に於ける処理を示している。

（１）試料溶液の調製
本発明の光分析技術の観測対象物となる粒子は、溶解された分子等の、試料溶液中にて分散し溶液中にてランダムに運動する粒子であれば、任意のものであってよく、例えば、タンパク質、ペプチド、核酸、脂質、糖鎖、アミノ酸若しくはこれらの凝集体などの生体分子、ウイルス、細胞、或いは、金属コロイド、その他の非生物学的分子などであってよい。観測対象物となる粒子が発光粒子であり、その粒子が、元来、光を発する粒子でない場合には、発光標識（蛍光分子、りん光分子、化学・生物発光分子）が観測対象物となる粒子に任意の態様にて付加されたものが用いられる。試料溶液は、典型的には水溶液であるが、これに限定されず、有機溶媒その他の任意の液体であってよい。

反転型走査分子計数法を行う場合、観測対象物となる非発光粒子であるが、発光粒子の場合と同様に任意のものであってよい。なお、粒径が光検出領域の径の好適には１５％以上、より好適には、３５％以上であることが見出されている。また、背景光を与える発光物質としては、任意の発光分子、例えば、蛍光分子、りん光分子、化学・生物発光分子であってよく、発光物質は、光検出領域内に常に数分子以上存在する濃度にて試料溶液中に溶解又は分散される。試料溶液は、典型的には水溶液であるが、これに限定されず、有機溶媒その他の任意の液体であってもよい。

（２）試料溶液の光強度の測定（図５−ステップ１００）
本実施形態の走査分子計数法又は反転型走査分子計数法による光分析に於ける光強度の測定は、測定中にミラー偏向器１７又はステージ位置変更装置１７ａを駆動して、試料溶液内での光検出領域の位置の移動（試料溶液内の走査）を行う他は、ＦＣＳ又はＦＩＤＡに於ける光強度の測定過程と同様の態様にて実行されてよい。操作処理に於いて、典型的には、マイクロプレート９のウェル１０に試料溶液を注入して顕微鏡のステージ上に載置した後、使用者がコンピュータ１８に対して、測定の開始の指示を入力すると、コンピュータ１８は、記憶装置（図示せず）に記憶されたプログラム（試料溶液内に於いて光検出領域の位置を移動する手順と、光検出領域の位置の移動中に光検出領域からの光を検出して時系列の光強度データを生成する手順）に従って、試料溶液内の光検出領域に於ける励起光の照射及び光強度の計測が開始される。かかる計測中、コンピュータ１８のプログラムに従った処理動作の制御下、ミラー偏向器１７又はステージ位置変更装置１７ａは、ミラー７（ガルバノミラー）又は顕微鏡のステージ上のマイクロプレート９を駆動して、ウェル１０内に於いて光検出領域の位置の移動を実行し、これと同時に光検出器１６は、逐次的に検出された光を電気信号に変換してコンピュータ１８へ送信し、コンピュータ１８では、任意の態様にて、送信された信号から時系列の光強度データを生成して保存する。なお、典型的には、光検出器１６は、一光子の到来を検出できる超高感度光検出器であるので、光の検出が、フォトンカウンティングによる場合、時系列光強度データは、時系列のフォトンカウントデータであってよい。また、互いに異なる光成分を別々に検出する場合には、複数の光検出器１６の各々が、それらの光成分の光強度値（光子数）をそれぞれ同時に検出し検出された互いに異なる成分毎に時系列光強度データが生成される。

ところで、光検出領域の位置の移動速度は、好適には、単一粒子のランダムな運動、即ち、ブラウン運動による移動速度よりも速い値に設定される。光検出領域の位置の移動速度が粒子のブラウン運動による移動に比して遅い場合には、粒子が光検出領域内をランダムに移動することとなる。そうすると、粒子が存在している状態を仮定して算定する発生確率のための期待値の演算が複雑となり、また、精度も低下し得る。そこで、好適には、図３（Ｃ）に描かれている如く、粒子が光検出領域を略直線に横切り、光強度の変化のプロファイルの期待値が、図３（Ｄ）に描かれている如き粒子の位置が光検出領域の略中心から離れるほど低下するプロファイルと同様の釣鐘状となるように、光検出領域の位置の移動速度は、粒子のブラウン運動による平均の移動速度（拡散移動速度）よりも速く設定される。

具体的には、拡散係数Ｄを有する単一粒子がブラウン運動によって半径ｒの光検出領域（コンフォーカルボリューム）を通過するときに要する時間Δτは、平均二乗変位の関係式
（２ｒ）^２＝６Ｄ・Δτ …（２６）
から、
Δτ＝（２ｒ）^２／６Ｄ …（２７）
となるので、単一粒子がブラウン運動により移動する速度（拡散移動速度）Ｖdifは、概ね、
Ｖdif＝２ｒ／Δτ＝３Ｄ／ｒ …（２８）
となる。そこで、光検出領域の位置の移動速度は、かかるＶdifを参照して、それよりも十分に早い値に設定されてよい。例えば、単一粒子の拡散係数が、Ｄ＝２．０×１０^−１０ｍ^２／ｓ程度であると予想される場合には、ｒが、０．６２μｍ程度だとすると、Ｖdifは、１．０×１０^−３ｍ／ｓとなるので、光検出領域の位置の移動速度は、その１０倍以上の１５ｍｍ／ｓと設定されてよい。なお、単一粒子の拡散係数が未知の場合には、光検出領域の位置の移動速度を種々設定して光強度の変化のプロファイルが、予想されるプロファイル（典型的には、励起光強度分布と略同様）となる条件を見つけるための予備実験を繰り返し実行して、好適な光検出領域の位置の移動速度が決定されてよい。

（３）単一粒子の信号の個別検出（ステップ１１０〜１６０）
時系列光強度データが生成されると、まず、時系列光強度データ上にて背景光の強度値の算出が行われる（ステップ１１０）。背景光の強度値は、時系列光強度データ上に於いて、粒子の信号が存在していない領域の強度値（光子数）の平均値であってよい。かくして、背景光の強度値の算出の一つの手法に於いては、通常の走査分子計数法の場合（観測対象粒子が発光粒子である場合）には、得られた時系列光強度データ上の光強度値のうち、高い方から所定の割合（例えば、２０％）のデータを除いた全強度値の平均値が背景光の強度値として採用されてよい。これは、光強度値の高い方から所定の割合のデータが粒子の信号であると考えられるためである。また、反転型走査分子計数法の場合（観測対象粒子が非発光粒子である場合）には、得られた時系列光強度データ上の光強度値のうち、低い方から所定の割合（例えば、２０％）のデータを除いた全強度値の平均値が背景光の強度値として採用されてよい。これは、光強度値の低い方から所定の割合のデータが粒子の信号であると考えられるためである。時系列光強度データが複数の成分のそれぞれに生成されているときは、成分毎に背景光の強度値の算出が行われる。なお、背景光の強度値は、観測対象の粒子を含まない試料溶液を用いて得られた時系列光強度データ上の光強度値の平均であってもよい。

次に、本実施形態の処理過程に於いては、時系列光強度データ上に於いて解析窓の設定が実行される（ステップ１２０）。既に触れた如く、一つの解析窓の長さは、光検出領域の大きさと移動速度とから決定されてよい（式（９）参照）。また、解析窓は、好適には、ビンタイム毎に時系列に設定されてよい。しかしながら、演算量を低減する目的で、数個のビンタイム毎に設定されてもよく、或いは、解析窓が互いに重複しないように設定されてもよい。

かくして、解析窓の設定が為されると、上記に説明された原理に従って、光検出領域内に粒子が無いと仮定した場合の解析窓内の光強度値列又は光子数列の発生確率Ｐｎ（ステップ１３０）、光検出領域内に粒子が有ると仮定した場合の解析窓内の光強度値列又は光子数列の発生確率Ｐｐ（ステップ１４０）がそれぞれ算出され、発生確率Ｐｎと発生確率Ｐｐとが比較され、解析窓内に粒子が存在しているか否かが判定される（ステップ１５０）。かかる判定に於いては、発生確率Ｐｎと発生確率Ｐｐのオッズ比が算出され（式（８）参照）、オッズ比が所定値以上のとき、粒子が存在していると判定されてよい。互いに異なる複数の成分について時系列光強度データが得られ、粒子の所定の特性値を考慮して発生確率Ｐｐが得られている場合、所定の特性値を有する粒子の信号の有無が選択的に検出されることが理解されるべきである。

なお、光検出領域内に粒子が有ると仮定した場合の発生確率Ｐｐの算出に於いて、強度値の期待値は、解析窓の中心に強度のピークが存在するものと仮定されてよい。実際には、粒子の信号のピークは、解析窓の中心に存在しないことが殆どであるが、実際の粒子の信号のピークの位置が解析窓の中心から離れるほど、発生確率Ｐｐの値が低減するが、実際の粒子の信号が存在していないときの発生確率Ｐｐの値よりは、高い値となる。

上記のステップ１３０〜１５０の処理に於ける解析窓内の発生確率Ｐｐ、Ｐｎの算出及び粒子信号の検出は、光強度データに設定された全ての解析窓に於いて実行されてよい（ステップ１６０）。

（５）粒子濃度その他の分析
かくして、上記の処理により、各解析窓に於ける粒子信号の有無が判定されると、光強度データ上に於ける粒子信号の計数や濃度の算出等の分析が行われてよい（ステップ１７０）。既に触れた如く、各解析窓に於ける粒子信号の有無が判定されると、一つの粒子信号が存在する解析窓が連続することとなる。従って、連続する粒子信号が存在する解析窓の組を計数することにより、粒子の数が得られることとなる。

また、粒子の数の決定が為されている場合、任意の手法にて、光検出領域の通過した領域の総体積が算定されれば、その体積値と単一粒子の数とから試料溶液中の単一粒子の数密度又は濃度が決定される。光検出領域の通過した領域の総体積は、励起光又は検出光の波長、レンズの開口数、光学系の調整状態に基づいて理論的に算定されてよい。また、実験的に、例えば、粒子濃度が既知の溶液（対照溶液）について、検査されるべき試料溶液の測定と同様の条件にて、上記に説明した光強度の測定、単一粒子の検出及び計数を行うことにより検出された粒子の数と、対照溶液の粒子の濃度とから決定されるようになっていてよい。具体的には、例えば、粒子濃度Ｃの対照溶液について、対照溶液の粒子の検出数がＮであったとすると、光検出領域の通過した領域の総体積Ｖｔは、
Ｖｔ＝Ｎ／Ｃ …（２９）
により与えられる。また、対照溶液として、粒子の複数の異なる濃度の溶液が準備され、それぞれについて測定が実行されて、算出されたＶｔの平均値が光検出領域の通過した領域の総体積Ｖｔとして採用されるようになっていてよい。そして、Ｖｔが与えられると、粒子の計数結果がｎの試料溶液の発光粒子の濃度ｃは、
ｃ＝ｎ／Ｖｔ …（３０）
により与えられる。なお、光検出領域の体積、光検出領域の通過した領域の総体積は、上記の方法によらず、任意の方法にて、例えば、ＦＣＳ、ＦＩＤＡを利用するなどして与えられるようになっていてよい。また、本実施形態の光分析装置に於いては、想定される光検出領域の移動パターンについて、種々の標準的な粒子についての濃度Ｃと粒子の数Ｎとの関係（式（２９））の情報をコンピュータ１８の記憶装置に予め記憶しておき、装置の使用者が光分析を実施する際に適宜記憶された関係の情報を利用できるようになっていてよい。なお、互いに異なる複数の成分について時系列光強度データが得られ、粒子の所定の特性値を考慮して発生確率Ｐｐが得られている場合、所定の特性値を有する粒子の信号の計数と濃度の算出が可能であることが理解されるべきである。従って、試料溶液中に互いに異なる特性値を有する複数の種類の粒子が混在している場合、粒子の種類を識別しながら、種類毎に粒子の計数と濃度の算出が可能である

上記に説明した本発明の有効性を検証するために、以下の如き実験を行った。なお、以下の実施例は、本発明の有効性を例示するものであって、本発明の範囲を限定するものではないことは理解されるべきである。

走査分子計数法に於いて、本発明による時系列光強度データ上の光子数列の発生確率Ｐｐ、Ｐｎに基づいて粒子信号の検出を行った。

試料溶液として、トリス緩衝液（１０ｍＭ Tris-HCl(pH8.0)、０．０５％ Tween20を含む）に発光粒子としてＡＴＴＯ６４７Ｎ（シグマ）を種々の濃度にて含む溶液を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。

光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０〜１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を検出し計数した。なお、本実施例に於いては、時系列光強度データは、一つの検出波長帯域の成分のみであるので、発生確率Ｐｎ、発生確率Ｐｐは、「２．単一粒子の信号の検出」の欄に記載された式（１）、（２）、（６）等を用いて算出した。

図６は、時系列光強度データの一部（Ａ）と、（Ａ）を平滑化して得られたデータ（Ｂ）と、（Ａ）から本発明の教示に従って算出された発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）の一つの例を示している。同図のデータの取得に於いては、励起光は、６３３ｎｍのレーザー光を用い、バンドパスフィルターを用いて、６６０−７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、３０００ｒｐｍ（７．５ｍｍ／秒）とし、ＢＩＮ ＴＩＭＥを１０μ秒とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、図示の時間領域の周辺１ｍ秒の全データの光子数のうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、８０μ秒とした。また、解析窓の時間幅は、３００μ秒とした。試料溶液中のＡＴＴＯ６４７Ｎの濃度は、１ｐＭとした。

図６を参照して、平滑化して得られたデータ（Ｂ）に於いて、特許文献９〜１１に記載されている如き、平滑化データに於ける信号強度の大きさが所定の閾値による粒子信号の判定処理では、パルス状の信号の強度が所定の閾値を越えたか否かにより、粒子の信号か否かが判定される。所定の閾値は、典型的には、光子数１に設定されるので、図の略中央の（ａ）にて示した光子数の増大の除き、全てノイズとして判定される。（ａ）にて示したパルス状信号のピーク値は、光子数〜１であるので、粒子信号であると思料されるところ、判定結果は誤差の影響で変わってしまい、ノイズ信号として判定される場合もあり得る（即ち、結果が不安定となる。）。

一方、発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）に於いては、（ａ）にて示した信号に対応する領域に於いて、値の顕著な増大が見られ、従って、安定的に粒子信号であると判定できることとなる。

図７は、図６と同様の時系列光強度データの別の一部を示した図である。同図のデータの取得に於いては、試料溶液中に於ける光検出領域の位置の移動速度は、６０００ｒｐｍ（１５ｍｍ／秒）とした。図７を参照すると、平滑化して得られたデータ（Ｂ）に於いて、（ａ）にて示したパルス状信号の他に有意なパルス状の信号が観察されるが、発生確率Ｐｎ、Ｐｐのオッズ比（Ｃ）に於いて（ａ）にて示したパルス状信号に対応する領域以外には、有意な増大は見られなかった。このことは、本発明の教示に従って算出された発生確率Ｐｎ、Ｐｐに基づく判定によれば、粒子信号とノイズ信号との識別が明確となり、粒子判定の安定性が向上されることを示唆している。

図８は、種々の濃度の発光粒子を含む試料溶液について本発明による発生確率Ｐｎ、Ｐｐのオッズ比を用いて粒子信号の判定を行った場合の発光粒子濃度と検出粒子数との関係を示している。なお、本実験に際しては、励起光は、６４２ｎｍのレーザー光（３ｍＷ）を用い、バンドパスフィルターを用いて、５６０−６２０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、１２０００ｒｐｍ（９０ｍｍ／秒）とし、ＢＩＮ ＴＩＭＥを１０μ秒とし、光の測定時間を６００秒間とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、５０μ秒とした。また、解析窓の時間幅は、３５０μ秒とした。そして、オッズ比＞１０^８を越える領域に粒子信号が存在するものと判定した。図８から理解される如く、１ｆＭ以上の濃度域に於いて、検出粒子数は、発光粒子濃度と共に増大した。このことは、本発明の粒子信号の検出手法により、時系列フォトンカウントデータ上の粒子の信号を検出できること、また、試料溶液中の粒子濃度を１ｆＭ以上の濃度域に於いて決定できることを示している。

反転型走査分子計数法に於いて、本発明による時系列光強度データ上の光子数列の発生確率Ｐｐ、Ｐｎに基づいて粒子信号の検出を行った。

試料溶液として、背景光を生成する発光物質として、１ｍＭ ＡＴＴＯ６４７Ｎを含むトリス緩衝液（１０ｍＭ Tris-HCl(pH8.0)、０．０５％ Tween20を含む）にポリスチレンビーズ（直径４μｍ）を分散した溶液を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。

光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０〜１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を計数した。なお、本実施例に於いては、時系列光強度データは、一つの検出波長帯域の成分のみであるので、発生確率Ｐｎ、発生確率Ｐｐは、「２．単一粒子の信号の検出」の欄に記載された式（１）、（２）、（７）等を用いて算出した。

図９（Ａ）、（Ｂ）、図１０は、得られた時系列光強度データの一部を順に拡大した図（上段）とそのオッズ比を示している。同図のデータの取得に於いては、励起光は、６３３ｎｍのレーザー光（５０μＷ）を用い、バンドパスフィルターを用いて、６６０−７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、６０００ｒｐｍ（１５ｍｍ／秒）とし、ＢＩＮ ＴＩＭＥを５０μ秒とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１００μ秒とした。また、解析窓の時間幅は、３００μ秒とした。ビーズの濃度は、１００ｆＭとした。

これらの図を参照して、オッズ比の増大は、フォトンカウントデータ上の背景光の低減部分と対応していることが理解される。これにより、本発明の粒子信号検出の手法は、反転型走査分子計数法にも適用可能であることが示唆された。

図１１（Ａ）は、種々の濃度のビーズを含む試料溶液について本発明による発生確率Ｐｎ、Ｐｐのオッズ比を用いて粒子信号の判定を行った場合の粒子濃度と検出粒子数との関係を示している。なお、本実験に際しては、励起光は、６３３ｎｍのレーザー光（５０μＷ）を用い、バンドパスフィルターを用いて、６６０−７１０ｎｍの波長帯域の光を測定し、時系列フォトンカウントデータを生成した。試料溶液中に於ける光検出領域の位置の移動速度は、９０００ｒｐｍ（９０ｍｍ／秒）とし、ＢＩＮ ＴＩＭＥを１０μ秒とし、光の測定時間を２００秒間とした。発生確率Ｐｎ、Ｐｐの算出に於いては、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１００μ秒とした。また、解析窓の時間幅は、３００μ秒とした。そして、オッズ比＞１０^１０を越える領域に粒子信号が存在するものと判定した。なお、比較のため、図１１（Ｂ）は、同一のデータについて、特許文献９〜１１に記載の如く時系列フォトンカウントデータを平滑化した後、ガウス関数をフィッティングしてピーク光子数≧５であるパルス状信号を検出して得られた検出粒子数と粒子濃度との関係を示している。

図１１（Ａ）、（Ｂ）を参照して、図示の濃度域（１ａＭ〜１ｆＭ）に於いて、平滑化された時系列フォトンカウントデータに於いてガウス関数をフィッティングして得られた結果（Ｂ）の場合、濃度の増大に対して検出粒子数の増大が観察されるが、ばらつきが比較的大きかった（近似直線の相関係数ｒ^２＝０．９３）。これに対し、本発明による粒子信号の検出手法による結果（Ａ）に於いては、検出粒子数は、粒子濃度に概ね比例した（近似直線の相関係数ｒ^２＝０．９９）。また、同図から、試料溶液中の粒子濃度を５０ａＭ以上の濃度域に於いて決定できることを示している。（（Ｂ）の場合は、決定可能な濃度範囲は、５００ａＭ以上である。）

蛍光異方性の異なる二種類の発光粒子について、本発明による走査分子計数法により種類毎に粒子の検出が可能であることを検証した。

試料溶液として、リン酸緩衝液（０．０５％ Tween20を含む）に蛍光色素TAMRA(Ｍ．Ｗ．４３０．４５ Sigma-Aldrich Cat.No.C2734)を１００ｆＭにて含む溶液（蛍光色素溶液）と、リン酸緩衝液に１ｐＭのプラスミド（pbr322 ２．９ＭＤａ タカラバイオ株式会社 Cat.No.3035）と１０ｎＭのＤＮＡインターカレータ蛍光色素SYTOX Orange（インビトロジェン社 Cat.No.S-11368）を含む溶液（プラスミド溶液、SYTOX Orangeが単一のプラスミドに結合して単一の発光粒子となる。）を調製した。光の測定に於いては、光分析装置として、共焦点蛍光顕微鏡の光学系とフォトンカウンティングシステムを備えた１分子蛍光測定装置ＭＦ２０（オリンパス株式会社）を用い、上記の「（２）試料溶液の光強度の測定」にて説明した要領に従って、上記の試料溶液について、ｓ偏光成分及びｐ偏光成分の時系列光強度データ（フォトンカウントデータ）を同時に且つ別々に取得した。その際、ＴＡＭＲＡ溶液及びプラスミド溶液の両方について、励起光は、５４３ｎｍのレーザー光を用い、バンドパスフィルターを用いて、５６０−６２０ｎｍの波長帯域の光を検出した。また、励起光の偏光方向は、検出光のｐ偏光成分と同じ方向となるよう設定した。試料溶液中に於ける光検出領域の移動速度は、６０００rpm（１５ｍｍ／秒）と設定し、ＢＩＮ ＴＩＭＥを１０μ秒とし、測定時間は、２秒間とした。上記の試料溶液について、時系列光強度データ（フォトンカウントデータ）を取得した。

光計測後のデータ処理に於いては、まず、取得された時系列フォトンカウントデータに於いて、「（３）単一粒子の信号の個別検出」及び図５のステップ１１０〜１６０に記載された要領に従って、背景光強度の算出、解析窓の設定、光検出領域内に粒子が無いと仮定した場合の解析窓内の光子数列の発生確率Ｐｎ及び光検出領域内に粒子が有ると仮定した場合の解析窓内の光子数列の発生確率Ｐｐの算出、発生確率のオッズ比の算出を順に実行し、単一粒子の信号を計数した。なお、本実施例に於いては、ｐ波とｓ波の二成分の時系列光強度データから、蛍光異方性を既知のパラメータとして含む発生確率Ｐｎ、発生確率Ｐｐを、「３．複数の光成分の計測を用いた単一粒子信号の検出（ｉ）特定の偏光特性を有する単一粒子の信号の検出」の欄に記載された式（１５）、（１６）等を用いて算出した。その際、背景光の強度は、各解析窓の周辺１ｍ秒のデータのうちの光子数の高い方と低い方の２０％を除いた全データの光子数の平均値とした。粒子が存在していると仮定した場合の光子数の期待値の算出に於いて、期待値の分布を表すガウス関数（式（４））の半値全幅（ｄ／ｖ）は、１２０μ秒とした。また、解析窓の時間幅は、４００μ秒とした。また、発生確率Ｐｐに関しては、プラスミドの蛍光異方性の０．４を用いた値と、ＴＡＭＲＡの蛍光異方性の０．３２を用いた値とを算出し（即ち、一つの解析窓について、二つの発生確率Ｐｐが算出される。）、それぞれについて、発生確率Ｐｎとのオッズ比が２０を越えたとき、その解析窓に対応する種類の粒子の信号が存在するものと判定した。（プラスミドの発生確率Ｐｐのオッズ比とＴＡＭＲＡの発生確率Ｐｐのオッズ比とが共に２０を越えたときには、オッズ比の高い種類の粒子信号が存在するものとした。）

図１２は、プラスミド溶液と蛍光色素溶液について上記の走査分子計数法により得られた時系列光強度データに於いて、発生確率Ｐｎと互いに異なる蛍光異方性を用いて得られた発生確率Ｐｐとのオッズ比を用いて、粒子信号の存在を判定し、粒子の種類毎に粒子数を計数した結果を示している。同図を参照して、図から理解される如く、それぞれの溶液に於いて、大半の検出された粒子信号がそれぞれの種類に対応して検出された。（プラスミド溶液中では、４２１個の粒子信号が検出され、そのうち、ＴＡＭＲＡとして誤検出された数は、８２個であった。一方、蛍光色素溶液中では、２７９個の粒子信号が検出され、そのうち、プラスミドとして誤検出された数は、５８個であった。）この結果は、本発明の手法により、種類を識別しつつ試料溶液中の単一粒子の検出が可能であることを示唆している。

かくして、上記の実施例の結果から理解される如く、本発明の教示に従って、走査分子計数法に於いて時系列光強度データに於ける光強度値列の発生確率を、粒子が存在したと仮定した場合と粒子が存在していなかったと仮定した場合とのそれぞれについて、算出し、算出された発生確率に基づいて、粒子の信号の有無を判定することにより、粒子の信号の検出精度又はＳ／Ｎ比が向上され、また、検出感度（許容可能な精度にて粒子の信号が検出可能な試料溶液の粒子濃度の範囲）が改善できることが示された。

Claims (30)

1. 共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析装置であって、
   前記試料溶液内に於ける前記光学系の光検出領域の位置を移動する光検出領域移動部と、
   前記光検出領域からの光を検出する光検出部と、
   前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出部にて検出された前記光検出領域からの光の時系列の光強度データを生成し、前記時系列の光強度データに於いて前記単一粒子の各々の存在を表す信号を個別に検出する信号処理部とを含み、
   前記信号処理部が、前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出し、前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出することを特徴とする装置。
2. 請求項１の装置であって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とする装置。
3. 請求項２の装置であって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とする装置。
4. 請求項１乃至３の装置であって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とする装置。
5. 請求項１の装置であって、前記単一粒子が所定の特性値を有し、前記光検出部が前記光検出領域からの光の少なくとも二つの互いに異なる成分を別々に検出し、前記信号処理部が、前記成分の各々の時系列の光強度データを生成し、更に、前記信号処理部が、前記成分の各々の前記第一及び第二の発生確率を算出し、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とする装置。
6. 請求項５の装置であって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とする装置。
7. 請求項１乃至６の装置であって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とする装置。
8. 請求項５又は６を引用する請求項７の装置であって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とする装置。
9. 請求項５又は６を引用する請求項７の装置であって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とする装置。
10. 請求項１乃至４の装置であって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とする装置。
11. 共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出する光分析方法であって、
    前記試料溶液内に於いて前記光学系の光検出領域の位置を移動する過程と、
    前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出領域からの光の強度を測定して光強度データを生成する過程と、
    前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度値の時間変化の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する過程と、
    前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出する過程と
    を含むことを特徴とする方法。
12. 請求項１１の方法であって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とする方法。
13. 請求項１２の方法であって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とする方法。
14. 請求項１１乃至１３の方法であって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とする方法。
15. 請求項１１の方法であって、前記単一粒子が所定の特性値を有し、前記光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、前記成分の各々の時系列の光強度データが生成され、更に、前記成分の各々の前記第一及び第二の発生確率が算出され、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とする方法。
16. 請求項１５の方法であって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とする方法。
17. 請求項１１乃至１６の方法であって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とする方法。
18. 請求項１５又は１６を引用する請求項１７の方法であって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とする方法。
19. 請求項１５又は１６を引用する請求項１７の方法であって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とする方法。
20. 請求項１１乃至１４の方法であって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とする方法。
21. 共焦点顕微鏡又は多光子顕微鏡の光学系を用いて試料溶液中にて分散しランダムに運動する単一粒子を検出するための光分析用コンピュータプログラムであって、
    前記試料溶液内に於いて前記光学系の光検出領域の位置を移動する手順と、
    前記試料溶液内に於いて前記光検出領域の位置を移動させながら前記光検出領域からの光の強度を測定して光強度データを生成する手順と、
    前記時系列の光強度データ上に時系列に設定される解析窓の各々に於ける光強度の時間変化値の、前記光検出領域内に前記単一粒子が存在していない第一の状態を仮定した場合に於ける第一の発生確率と、前記光検出領域内に前記単一粒子が存在している第二の状態を仮定した場合に於ける第二の発生確率とを算出する手順と、
    前記第一及び第二の発生確率に基づいて前記時系列の光強度データ上にて前記単一粒子の各々の存在を表す信号を検出する手順と
    をコンピュータに実行させることを特徴とするコンピュータプログラム。
22. 請求項２１のコンピュータプログラムであって、前記第一及び第二の発生確率が、それぞれ、前記解析窓内に於ける単位時間毎の光強度値と前記第一及び第二の状態を仮定した場合の前記単位時間毎の期待値とに基づいて算出されることを特徴とするコンピュータプログラム。
23. 請求項２２のコンピュータプログラムであって、前記単位時間毎の光強度値が前記単位時間毎の前記期待値を有するポアソン分布に従うものとして前記単位時間毎の光強度値の単位時間発生確率が算出され、前記第一及び第二の発生確率がそれぞれ対応する前記単位時間発生確率を用いて算出されることを特徴とするコンピュータプログラム。
24. 請求項２１乃至２３のコンピュータプログラムであって、前記第二の発生確率が前記第一の発生確率より大きい解析窓の時間に前記光検出領域内に前記単一粒子が存在したと判定されることを特徴とするコンピュータプログラム。
    
25. 請求項２１のコンピュータプログラムであって、前記単一粒子が所定の特性値を有し、前記光検出領域からの光の少なくとも二つの互いに異なる成分が別々に検出されて、前記成分の各々の時系列の光強度データが生成され、更に、前記成分の各々の前記第一及び第二の発生確率が算出され、前記成分の各々の前記第二の発生確率が前記所定の特性値の関数であり、前記成分毎の前記第一及び第二の発生確率の発生確率に基づいて前記時系列の光強度データ上にて前記所定の特性値を有する前記単一粒子の存在を表す信号が検出されることを特徴とするコンピュータプログラム。
26. 請求項２５のコンピュータプログラムであって、前記単一粒子が互いに異なる所定の特性値を有する複数の種類の単一粒子を含み、前記単一粒子の種類毎に、前記互いに異なる所定の特性値の関数である前記成分の各々の前記第二の発生確率が算出され、前記成分の各々の前記第一の発生確率と前記複数の種類の単一粒子の各々の前記成分の各々の前記第二の発生確率とに基づいて前記時系列の光強度データ上にて前記単一粒子の種類毎に該単一粒子の存在を表す信号が検出されることを特徴とするコンピュータプログラム。
27. 請求項２１乃至２６のコンピュータプログラムであって、前記単一粒子が発光粒子であり、前記単一粒子の各々の存在を表す信号が光強度の一時的な増大であることを特徴とするコンピュータプログラム。
28. 請求項２５又は２６を引用する請求項２７のコンピュータプログラムであって、前記所定の特性値が前記単一粒子の偏光異方性であることを特徴とするコンピュータプログラム。
29. 請求項２５又は２６を引用する請求項２７のコンピュータプログラムであって、前記所定の特性値が前記単一粒子の互いに異なる発光波長帯域の発光強度の比であることを特徴とするコンピュータプログラム。
30. 請求項２１乃至２４のコンピュータプログラムであって、前記単一粒子が検出波長帯域にて発光しない粒子であり、前記光検出領域からの光が背景光を含み、前記単一粒子の各々の存在を表す信号が前記背景光からの光強度の一時的な低減であることを特徴とするコンピュータプログラム。

Images