JP2021121803A - 微小粒子測定システム及び微小粒子測定方法 - Google Patents

Abstract

【課題】測定誤差を簡便な処理によって補正可能な微小粒子測定システム及び微小粒子測定方法を提供する。【解決手段】流路内を通流する微小粒子に対して、レーザー光を照射する照射部と、前記微小粒子からの蛍光を複数の受光素子で検出する検出部と、前記検出部で検出された蛍光強度値に対する補正処理により生成した補正強度値よりスペクトルデータを生成し、前記スペクトルデータに対して基準スペクトルデータを用いたコンペンセーション処理を行う処理部と、を有する微小粒子測定システムであって、前記複数の受光素子は、少なくとも２以上の異なる波長域幅を有する、微小粒子測定システムなどを提供する。【選択図】図１

Description

本技術は、細胞等の微小粒子の特性を光学的に測定する微小粒子測定システム及び微小粒子測定方法に関する。

フローサイトメータは、フローセルを通流する細胞やビーズなどの微小粒子に光を照射し、微小粒子から発せられる蛍光及び散乱光などを検出することによって、各微小粒子の特性を光学的に測定する装置である。

例えば細胞の蛍光を検出する場合、蛍光色素により標識した細胞にレーザー光などの適当な波長かつ強度を有する励起光を照射する。そして、蛍光色素から発せられる蛍光をレンズなどで集光し、フィルタ又はダイクロイックミラーなどの波長選択素子を用いて適当な波長域の光を選択し、選択された光をＰＭＴ（photo multiplier tube）などの受光素子を用いて検出する。このとき、波長選択素子と受光素子とを複数組み合わせることによって、細胞に標識された複数の蛍光色素からの蛍光を同時に検出し、解析することも可能である。更に、複数波長の励起光を組み合わせることで解析可能な蛍光色素の数を増やすこともできる。

フローサイトメータにおける蛍光検出には、フィルタなどの波長選択素子を用いて不連続な波長域の光を複数選択し、各波長域の光の強度を計測する方法の他に、連続した波長域における光の強度を蛍光スペクトルとして計測する方法もある。蛍光スペクトルの計測が可能なスペクトル型フローサイトメータでは、微小粒子から発せられる蛍光を、プリズム又はグレーティングなどの分光素子を用いて分光する。そして、分光された蛍光を、検出波長域幅が異なる複数の受光素子が配列された受光素子アレイを用いて検出する。受光素子アレイには、ＰＭＴ又はフォトダイオードなどの受光素子を一次元に配列したＰＭＴアレイ又はフォトダイオードアレイ、あるいはＣＣＤ又はＣＭＯＳなどの２次元受光素子などの独立した検出チャネルが複数並べられたものが用いられている。

特開２００３−８３８９４号公報

フローサイトメータにおける計測値には様々な要因による誤差が含まれている。測定誤差の補正方法としては、例えば、蛍光特性が既知の標準サンプルを用いた方法が汎用されている。この方法は、標準サンプルを複数計測して得た基準値をもとに各受光素子について電流などの出力値と蛍光強度との関係（キャリブレーション情報）を取得し、この関係に基づいて計測値をキャリブレーションするものである。

上記の補正方法では、フローサイトメータによる測定の都度に標準サンプルを計測し、レーザー光の出力や受光素子の設定値（例えばＰＭＴであれば電圧など）に応じたキャリブレーション情報を取得する必要があり、非常に煩雑であった。

そこで、本技術は、測定誤差を簡便な処理によって補正可能な微小粒子測定システム及び微小粒子測定方法を提供することを主な目的とする。

上記課題解決のため、本技術は、流路内を通流する微小粒子に対して、レーザー光を照射する照射部と、
前記微小粒子からの蛍光を複数の受光素子で検出する検出部と、
前記検出部で検出された蛍光強度値に対する補正処理により生成した補正強度値よりスペクトルデータを生成し、前記スペクトルデータに対して基準スペクトルデータを用いたコンペンセーション処理を行う処理部と、
を有する微小粒子測定システムであって、
前記複数の受光素子は、少なくとも２以上の異なる波長域幅を有する、微小粒子測定システムを提供する。
また、本技術では、流路内を通流する微小粒子に対して、レーザー光を照射する照射工程と、
前記微小粒子からの蛍光を複数の受光素子で検出する検出工程と、
前記検出工程で検出された蛍光強度値に対する補正処理により生成した補正強度値よりスペクトルデータを生成し、前記スペクトルデータに対して基準スペクトルデータを用いたコンペンセーション処理を行う処理工程と、
を有する微小粒子測定方法であって、
前記複数の受光素子は、少なくとも２以上の異なる波長域幅を有する、微小粒子測定方法も提供する。

本技術において、「微小粒子」には、細胞や微生物、リポソームなどの生体関連微小粒子、あるいはラテックス粒子やゲル粒子、工業用粒子などの合成粒子などが広く含まれるものとする。

生体関連微小粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラ（細胞小器官）などが含まれる。細胞には、動物細胞（血球系細胞など）及び植物細胞が含まれる。微生物には、大腸菌などの細菌類、タバコモザイクウイルスなどのウイルス類、イースト菌などの菌類などが含まれる。更に、生体関連微小粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。また、工業用粒子は、例えば有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレン・ジビニルベンゼン、ポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカ、磁性体材料などが含まれる。金属には、金コロイド、アルミなどが含まれる。これら微小粒子の形状は、一般には球形であるのが普通であるが、非球形であってもよく、また大きさや質量なども特に限定されない。

本技術により、サンプル解析の都度に標準サンプルを計測しなくとも、測定誤差を補正して正確な解析結果を得ることが可能な微小粒子測定システム及び微小粒子測定方法が提供される。

本技術の第一実施形態に係る微小粒子測定装置Ａの機能的構成を示すブロック図である。 微小粒子測定装置Ａの検出部１０の構成を説明する模式図である。 本技術の第二実施形態に係る微小粒子測定装置Ｂの検出部１０の構成を説明する模式図である。 微小粒子測定装置Ｂに関し、補正処理前の強度値に基づく出力データ（Ａ）と補正処理後の強度値に基づく出力データ（Ｂ）の例を示す図面代用グラフである。 実施例で試作したフローサイトメータにおいて、ＰＭＴアレイの各ＰＭＴの検出波長域幅を決定した結果を示す図面代用グラフである。 実施例で試作したフローサイトメータにおいて、ＰＭＴアレイの各ＰＭＴの相対感度を算出した結果を示す図面代用グラフである。 実施例において蛍光分光光度計を用いた計測により得られた蛍光ビーズのスペクトルチャートを示す図面代用グラフである。 実施例において試作したフローサイトメータを用いた計測により得られた蛍光ビーズFPK505のスペクトルチャートを示す図面代用グラフである。（Ａ）は補正処理前のチャート、（Ｂ）は第一の補正強度値によるチャート、（Ｃ）は第二の補正強度値によるチャートを示す。 実施例において試作したフローサイトメータを用いた計測により得られた蛍光ビーズFPK505のスペクトルチャートを示す図面代用グラフである。（Ａ）は補正処理前のチャート、（Ｂ）は第一の補正強度値によるチャート、（Ｃ）は第二の補正強度値によるチャートを示す。 実施例において試作したフローサイトメータを用いた計測により得られた蛍光ビーズFPK528のスペクトルチャートを示す図面代用グラフである。（Ａ）は補正処理前のチャート、（Ｂ）は第一の補正強度値によるチャート、（Ｃ）は第二の補正強度値によるチャートを示す。 実施例において試作したフローサイトメータを用いた計測により得られた蛍光ビーズFPK549のスペクトルチャートを示す図面代用グラフである。（Ａ）は補正処理前のチャート、（Ｂ）は第一の補正強度値によるチャート、（Ｃ）は第二の補正強度値によるチャートを示す。 実施例において試作したフローサイトメータを用いた計測により得られた蛍光ビーズFPK667のスペクトルチャートを示す図面代用グラフである。（Ａ）は補正処理前のチャート、（Ｂ）は第一の補正強度値によるチャート、（Ｃ）は第二の補正強度値によるチャートを示す。

以下、本技術を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本技術の代表的な実施形態の一例を示したものであり、これにより本技術の範囲が狭く解釈されることはない。説明は以下の順序で行う。
１．第一実施形態に係る微小粒子測定装置
（１）装置の構成
（２）蛍光強度の補正処理
［第一の補正強度値の算出］
［第二の補正強度値の算出］
（３）データ表示
（４）データ解析
２．第二実施形態に係る微小粒子測定装置
（１）装置の構成
（２）蛍光強度の補正処理
［第一の補正強度値の算出］
［第二の補正強度値の算出］
（３）データ表示

１．第一実施形態に係る微小粒子測定装置
（１）装置の構成
図１は、本技術の第一実施形態に係る微小粒子測定装置Ａの機能的構成を説明するブロック図である。また、図２は、微小粒子測定装置Ａの検出部１０の構成を説明する模式図である。

微小粒子測定装置Ａは、微小粒子に対してレーザー光を照射して、微小粒子から発せられる蛍光を検出し、検出された蛍光の強度を電気信号に変換して出力する検出部１０と、ＣＰＵ２０、メモリ３０及びハードディスク（記憶部）４０と、を含む。微小粒子測定装置Ａにおいて、ＣＰＵ２０、メモリ３０及びハードディスク（記憶部）４０は処理部を構成する。また、微小粒子測定装置Ａは、ユーザインターフェースとして、マウス５１及びキーボード５２と、ディスプレイ６１及びプリンタ６２を含んで構成される表示部６０と、を含む。

検出部１０は、従来の微小粒子測定装置と同様の構成とできる。具体的には、光源１０１からのレーザー光を微小粒子Ｐに対して集光・照射する照射系と、微小粒子Ｐから発生する蛍光を分光する分光素子１０２と分光された光を検出する受光素子アレイ１０３とを含む検出系と、から構成される。微小粒子測定装置Ａにおいて、微小粒子Ｐは、フローセル内又はマイクロチップに形成された流路内を一列に配列されて通流される。

照射系は、光源１０１の他、微小粒子Ｐに対してレーザー光を集光・照射するための集光レンズ、ダイクロイックミラー及びバンドパスフィルター等（不図示）からなる。なお光源１０１は、互いに異なる波長を発光する光源を２つ以上組み合わせたものでもよく、その場合２つ以上のレーザー光が微小粒子Ｐを照射する場所は、同一であっても異なっていても構わない。また、検出系は、微小粒子Ｐから発生する蛍光を集光し、分光素子１０２に導光するための集光レンズ等（不図示）を含んでいてもよい。ここでは、受光素子アレイ１０３として、検出波長域幅が異なる３２チャネルのＰＭＴ（photo multiplier tube）を一次元に配列したＰＭＴアレイを用いた構成を例示した。なお、受光素子アレイ１０３には、フォトダイオードアレイ、あるいはＣＣＤ又はＣＭＯＳなどの２次元受光素子などの検出波長域幅が異なる独立した検出チャネルが複数並べられたものも用いることができる。

微小粒子測定装置Ａにおいて、検出部１０は、蛍光の他、レーザー光の照射によって微小粒子Ｐから発生する光であって、例えば前方散乱光、側方散乱光、レイリー散乱及びミー散乱等の散乱光なども検出するように構成してもよい。

（２）蛍光強度の補正処理
ＣＰＵ２０及びメモリ３０は、ハードディスク４０に格納された蛍光強度補正プログラム４１とＯＳ４３と共働して、検出部１０から出力される電気信号に基づき、蛍光の強度値の補正処理を行う。この補正処理は、蛍光の強度値を各受光素子（ここではチャネル１〜３２のＰＭＴ）の検出波長域幅で補正して第一の補正強度値を算出するステップと、第一の補正強度値を各ＰＭＴの感度データを用いて補正して第二の補正強度値を算出するステップとを含む。

［第一の補正強度値の算出］
第一の補正強度値の算出は、各ＰＭＴで取得された蛍光の強度値を、それぞれのＰＭＴの検出波長域幅で除すことにより行われる。

具体的には、チャネル１〜３２のＰＭＴのうち、チャネルｋのＰＭＴで得られたｎ番目の強度値をＩ［ｋ、ｎ］とし、チャネルｋのＰＭＴの検出下限波長をＬ［ｋ］かつ検出上限波長をＨ［ｋ］とする。この場合、第一の補正強度値Ｊ_１［ｋ、ｎ］は、次の式により算出される。ここで、ｋは１〜３２の整数を表す。
Ｊ_１［ｋ、ｎ］＝Ｉ［ｋ、ｎ］／（Ｈ［ｋ］−Ｌ［ｋ］）

チャネル１〜３２のＰＭＴで検出される光の波長域幅は、分光素子１０２を含む検出部１０の光学系が非線形性を有している場合、各ＰＭＴ間で異なった幅となる（後掲の図５参照）。このため、各ＰＭＴで取得される蛍光の強度値は、検出波長域幅が広いチャネルでは相対的に高く、検出波長域幅が狭いチャネルでは相対的に低くなり、このことが測定誤差の要因となる。

特に、分光素子１０２により分光された蛍光を受光素子アレイ１０３を用いて検出し、微小粒子Ｐの蛍光スペクトルを計測する場合、各ＰＭＴで取得される蛍光の強度値をそのまま用いた蛍光スペクトルでは、上記測定誤差によってスペクトル形状に歪みが生じる。すなわち、横軸をチャネル番号、縦軸を強度値とした２次元グラフ（以下、「スペクトルチャート」と称する）を、横軸を検出波長、縦軸を強度値としたスペクトルチャートと比較すると、前者では検出波長域幅が広いチャネルほど強度値が相対的に大きくなる。このため、両者の蛍光スペクトル形状は一致せず、ずれが生じることとなる。

各ＰＭＴで取得された蛍光の強度値をそれぞれのＰＭＴの検出波長域幅で除して得られる第一の補正強度値では、このような光学系の非線形性に起因した測定誤差を補償することが可能となる。

各ＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）は、検出部１０を構成する分光素子１０２、集光レンズ、ダイクロイックミラー及びバンドパスフィルターなどの光学素子の種類や配置によって一意に決定される（後掲の図５参照）。このため、光学素子の選定及び配置を含む装置設計が完了した段階で、各ＰＭＴの検出波長域幅を取得しておくことで、各ＰＭＴで取得された蛍光の強度値から第一の補正強度値を算出することが可能となる。

［第二の補正強度値の算出］
第二の補正強度値の算出は、各ＰＭＴにおける第一の補正強度値を、それぞれのＰＭＴの相対感度で除すことにより行われる。

具体的には、チャネル１〜３２のＰＭＴのうち、チャネルｋのＰＭＴの相対感度をＳ［ｋ］とする。この場合、第二の補正強度値Ｊ_２［ｋ、ｎ］は、次の式により算出される。
Ｊ_２［ｋ、ｎ］＝Ｊ_１［ｋ、ｎ］／Ｓ［ｋ］

ここで、相対感度とは、同一強度及び同一波長の光をＰＭＴに照射して各チャネルで得られた強度値を、最も高い強度値が得られたチャネルの強度値に対する相対値により示したものである。相対感度は、同一強度及び同一波長の光をＰＭＴに照射した場合に各チャネルから出力される電気信号量を記録した感度データから予め算出しておくことができる。この感度データは、各ＰＭＴに内在する感度差と、ユーザによって各ＰＭＴに設定された感度差（ゲイン）の両方が反映されたものである。なお、ゲインは、ユーザが印加電圧などの設定値を変化させることによって適宜調整が可能である。

チャネル１〜３２のＰＭＴの感度は、ＰＭＴの個体差及びゲインの設定差によって、各ＰＭＴ間で異なる（後掲の図６参照）。このため、各ＰＭＴで取得される蛍光の強度値は、感度が高いチャネルでは相対的に高く、感度が低いチャネルでは相対的に低くなり、このことも測定誤差の要因となっている。

特に、分光素子１０２により分光された蛍光を受光素子アレイ１０３を用いて検出し、微小粒子Ｐの蛍光スペクトルを計測する場合、各ＰＭＴで取得される蛍光の強度値をそのまま用いた蛍光スペクトルでは、上記測定誤差によってスペクトル形状に歪みが生じる。すなわち、感度が高いチャネルほど強度値が相対的に大きくなっているため、蛍光スペクトル形状がいびつとなり、不正確なものとなる。

各ＰＭＴにおける第一の補正強度値を、それぞれのＰＭＴの相対感度で除して得られる第二の補正強度値では、このような受光素子間の感度差に起因した測定誤差を補償することが可能となる。

（３）データ表示
ＣＰＵ２０、メモリ３０及びハードディスク４０を含んで構成される処理部は、算出された第一の補正強度値又は第二の補正強度値を座標軸としたスペクトルチャートを生成して、表示部６０に出力する。

スペクトルチャートは、各ＰＭＴの検出波長を横軸にとり、第一の補正強度値を縦軸にとったものとできる（後掲の図８〜１２の（Ｂ）参照）。また、スペクトルチャートは、好ましくは、検出波長を横軸にとり、第二の補正強度値を縦軸にとったものとされる（後掲の図８〜１２の（Ｃ）参照）。

スペクトルチャートには、所定の検出波長において所定の蛍光強度値で検出された微小粒子の数（イベント数あるいは密度）に基づいて、強度値を平均値や標準誤差、中央値、四分位点などの統計的な数値で表示することができる（後掲の図９参照）。更に、スペクトルチャートを、イベント数をとった座標軸を追加した３次元グラフとして表示することもでき、この３次元グラフを疑似３Ｄ表示することもできる。また、微小粒子の数（イベント数あるいは密度）に関する情報（頻度情報）を反映させた色相、彩度及び／又は明度によりスペクトルチャートを多色表示することも可能である（後掲の図８参照）。

（４）データ解析
ＣＰＵ２０、メモリ３０及びハードディスク４０を含んで構成される処理部は、算出された第一の補正強度値又は第二の補正強度値を第一のパラメータ、各ＰＭＴの検出波長を第二のパラメータとするスペクトルデータを生成し、これを用いて各種の解析を実行することができる。ここで第一のパラメータは、計測を行った複数の微小粒子Ｐの全てあるいは一部について計算した、第一の補正強度値又は第二の補正強度値の平均値、標準誤差、中央値及び四分位点などの統計的な数値とできる。

また、処理部は、生成されたスペクトルデータと、ハードディスク（記憶部）４０に保持された基準スペクトルデータとを比較して、両者の一致度を評価することができる。更に、その結果を表示部６０に出力することもできる。ここで、基準スペクトルデータとは、既知の蛍光物質を含んでなる微小粒子を事前に微小粒子測定装置Ａで測定し、上述の補正処理を施すことによって得たスペクトルデータか、あるいは当該蛍光物質の蛍光スペクトルを通常の蛍光分光光度計により計測して得たスペクトルデータの、いずれであってもよい。これは、微小粒子測定装置Ａより得られた計測値が、補正によって通常の蛍光分光光度計で得られた計測値と直接比較可能となるためである。基準スペクトルデータは、ハードディスク４０内にリファレンスデータ４２として蓄積されている。基準スペクトルデータとの一致度の評価方法には、例えば各検出波長における差分の和、差分の絶対値の和、差分の二乗和などを利用できる。

微小粒子Ｐを微小粒子測定装置Ａで測定し、これを補正して得たスペクトルデータについて、１ないし２以上の基準スペクトルデータとの比較を行うことにより、測定を行った微小粒子Ｐが発する蛍光が、基準スペクトルデータに記録された蛍光のうちのいずれかと類似するか否かを判定することができる。これにより、例えば測定を行った微小粒子Ｐに含まれる蛍光物質の種類が不明な場合において、記録された基準スペクトルデータの中から類似度が高いものを探索して、微小粒子Ｐに含まれる蛍光物質の種類を予測することが可能となる。

更に、微小粒子Ｐを微小粒子測定装置Ａで測定し、これを補正して得たスペクトルデータについて、同一の微小粒子Ｐについてあらかじめ計測して得ておいた基準スペクトルデータとの比較を行うことにより、微小粒子測定装置Ａの状態診断を行うこともできる。すなわち、微小粒子測定装置Ａは、受光素子アレイの動作異常、フローセルないしマイクロチップ内における微小粒子の流れの乱れ、レンズや分光素子など各部品の温度変化や振動などによるずれなどの影響によってその状態が悪化し、計測の正確性が低下する場合がある。このような装置状態を診断する目的で、同一サンプルを微小粒子測定装置Ａと蛍光分光光度計とで測定し、その結果の比較を行うことは、有力な方法であるといえる。本技術によれば、微小粒子測定装置Ａで得られた計測結果を蛍光分光光度計と直接比較可能となるため、状態診断を簡便かつ正確に行うことが可能となる。

加えて、微小粒子Ｐを微小粒子測定装置Ａで測定し、これを補正して得たスペクトルデータに対して、複数の基準スペクトルデータを用いたコンペンセーション処理を行うことにより、微小粒子Ｐに含まれる複数の蛍光色素の定量を行うことができる。例えば微小粒子Ｐが複数の色素｛D₁, D₂, …, D_n｝で染色されている場合、コンペンセーション処理に用いる基準スペクトルデータには、｛D₁, D₂, …, D_n｝の各蛍光色素で単染色されたサンプルより得られたものが含まれることが望ましい。コンペンセーション処理の方法としては、例えば最小二乗法を採用することができる。

２．第二実施形態に係る微小粒子測定装置
（１）装置の構成
図３は、本技術の第二実施形態に係る微小粒子測定装置Ｂの検出部１０の構成を説明する模式図である。微小粒子測定装置Ｂの機能的構成は、図１に示した第一実施形態に係る微小粒子測定装置Ａと同様であるので説明を省略する。

微小粒子測定装置Ｂは、フィルタなどの波長選択素子を用いて不連続な波長域の光を複数選択し、各波長域の光の強度を計測する点で、連続した波長域における光の強度を蛍光スペクトルとして計測する微小粒子測定装置Ａと異なっている。

微小粒子測定装置Ｂの検出部１０は、光源１０１からのレーザー光を微小粒子Ｐに対して集光・照射する照射系と、微小粒子Ｐから発生する蛍光から所定の波長域の光を選択する波長選択素子１０４〜１０６と選択された光を検出する受光素子１０７〜１１０とを含む検出系と、から構成される。微小粒子測定装置Ｂにおいて、微小粒子Ｐは、フローセル内又はマイクロチップに形成された流路内を一列に配列されて通流される。

照射系は、光源１０１の他、微小粒子Ｐに対してレーザー光を集光・照射するための集光レンズ、ダイクロイックミラー及びバンドパスフィルター等（不図示）からなる。なお光源１０１は、互いに異なる波長を発光する光源を２つ以上組み合わせたものでもよく、その場合２つ以上のレーザー光が微小粒子Ｐを照射する場所は、同一であっても異なっていても構わない。また、検出系は、微小粒子Ｐから発生する蛍光を集光し、波長選択素子１０４に導光するための集光レンズ等（不図示）を含んでいてもよい。ここでは、受光素子１０７〜１１０として、検出波長域幅が異なるＰＭＴ（photo multiplier tube）を用いた構成を例示した。以下、受光素子１０７〜１１０のチャネル番号を１〜４とする。なお、受光素子１０７〜１１０には、フォトダイオードを用いることもできる。

微小粒子測定装置Ｂにおいて、検出部１０を、蛍光に加えて散乱光なども検出するように構成できる点は、微小粒子測定装置Ａと同様である。

（２）蛍光強度の補正処理
微小粒子測定装置Ｂは、検出部１０から出力される電気信号に基づき、蛍光の強度値の補正処理を行う。この補正処理は、蛍光の強度値を各受光素子（ここではチャネル１〜４のＰＭＴ）の検出波長域幅で補正して第一の補正強度値を算出するステップと、第一の補正強度値を各ＰＭＴの感度データを用いて補正して第二の補正強度値を算出するステップとを含む。

［第一の補正強度値の算出］
第一の補正強度値を算出は、各ＰＭＴで取得された蛍光の強度値を、それぞれのＰＭＴの検出波長域幅で除すことにより行われる。

具体的には、チャネル１〜４のＰＭＴのうち、チャネルｋのＰＭＴで得られたｎ番目の強度値をＩ［ｋ、ｎ］とし、チャネルｋのＰＭＴの検出下限波長をＬ［ｋ］かつ検出上限波長をＨ［ｋ］とする。この場合、第一の補正強度値Ｊ_１［ｋ、ｎ］は、次の式により算出される。ここで、ｋは１〜４の整数を表す。
Ｊ_１［ｋ、ｎ］＝Ｉ［ｋ、ｎ］／（Ｈ［ｋ］−Ｌ［ｋ］）

チャネル１〜４のＰＭＴで検出される光の波長域幅は、波長選択素子１０４〜１０６を含む検出部１０の光学系が非線形性を有している場合、各ＰＭＴ間で異なった幅となる。このため、各ＰＭＴで取得される蛍光の強度値は、検出波長域幅が広いチャネルでは相対的に高く、検出波長域幅が狭いチャネルでは相対的に低くなり、このことが測定誤差の要因となる。

各ＰＭＴで取得された蛍光の強度値をそれぞれのＰＭＴの検出波長域幅で除して得られる第一の補正強度値では、このような光学系の非線形性に起因した測定誤差を補償することが可能となる。

各ＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）は、検出部１０を構成する波長選択素子１０４〜１０６、集光レンズ、ダイクロイックミラー及びバンドパスフィルターなどの光学素子の種類や配置によって一意に決定される。このため、光学素子の選定及び配置を含む装置設計が完了した段階で、各ＰＭＴの検出波長域幅を取得しておくことで、各ＰＭＴで取得された蛍光の強度値から第一の補正強度値を算出することが可能となる。

［第二の補正強度値の算出］
第二の補正強度値の算出は、各ＰＭＴにおける第一の補正強度値を、それぞれのＰＭＴの相対感度で除すことにより行われる。

具体的には、チャネル１〜４のＰＭＴのうち、チャネルｋのＰＭＴの相対感度をＳ［ｋ］とする。この場合、第二の補正強度値Ｊ_２［ｋ、ｎ］は、次の式により算出される。
Ｊ_２［ｋ、ｎ］＝Ｊ_１［ｋ、ｎ］／Ｓ［ｋ］

チャネル１〜４のＰＭＴの感度は、ＰＭＴの個体差及びゲインの設定差によって、各ＰＭＴ間で異なる。このため、各ＰＭＴで取得される蛍光の強度値は、感度が高いチャネルでは相対的に高く、感度が低いチャネルでは相対的に低くなり、このことも測定誤差の要因となっている。

各ＰＭＴにおける第一の補正強度値を、それぞれのＰＭＴの相対感度で除して得られる第二の補正強度値では、このような受光素子間の感度差に起因した測定誤差を補償することが可能となる。

（３）データ表示
微小粒子測定装置Ｂは、算出された第一の補正強度値又は第二の補正強度値を座標軸とした２次元グラフを生成して、表示部６０に出力する。

２次元グラフは、各ＰＭＴの検出波長を横軸にとり、第一の補正強度値又は第二の補正強度値を縦軸にとったものとできる。図４（Ａ）に補正処理前のグラフ、（Ｂ）に第二の補正強度値によるグラフを示す。図４（Ａ）は、横軸にチャネル番号ｋ（ｋは１〜４の整数）、縦軸に各チャネルで取得された蛍光の強度値（Ｉ［ｋ］）の対数をとったグラフである。また、図４（Ｂ）は、横軸に検出波長、縦軸に各チャネルで取得された蛍光強度値の補正値（Ｊ_１［ｋ］又はＪ_２［ｋ］）の対数をとったグラフである。

微小粒子測定装置Ｂでは、上記の補正処理によって、波長選択素子を用いて選択した不連続な波長域の光の強度の計測結果を、装置の光学系の非線形性及び受光素子間の感度差に起因した測定誤差を補償した状態で表示させることが可能である。

本技術は、以下のような構成をとることもできる。
（１）微小粒子からの光を検出波長域が異なる複数の受光素子により検出する検出部から取得された前記光の強度値を、各受光素子の検出波長域幅で補正して、第一の補正強度値を算出する処理部、
を少なくとも有する、情報処理装置。
（２）前記検出部は、微小粒子からの光を分光する分光素子と、検出波長域が異なる複数の受光素子が配列された受光素子アレイと、を含む、（１）に記載の情報処理装置。
（３）前記処理部は、前記複数の受光素子により取得された光の強度値を、各受光素子が有する感度データに基づき補正処理した補正強度値を用いてスペクトルデータを生成する、（１）又は（２）に記載の情報処理装置。
（４）前記処理部は、前記第一の補正強度値を、各受光素子の感度データを用いて補正して、第二の補正強度値を算出する、（２）に記載の情報処理装置。
（５）前記処理部は、一の軸を前記検出波長とし、他の一の軸を前記第一の補正強度値又は前記第二の補正強度値とするスペクトルチャートを生成する、（２）又は（４）に記載の情報処理装置。
（６）前記処理部から出力された前記スペクトルチャートを表示する表示部を更に有する、（５）に記載の情報処理装置。
（７）前記処理部は、一のパラメータを前記検出波長とし、他の一のパラメータを前記第一の補正強度値又は前記第二の補正強度値とするスペクトルデータを生成し、基準スペクトルデータと比較し、両者が一致しているかあるいは不一致であるかを前記表示部に出力する、（６）に記載の情報処理装置。
（８）微小粒子からの光を検出波長域が異なる複数の受光素子により検出する検出部から取得された前記光の強度値を、各受光素子の検出波長域幅で補正して、第一の補正強度値を算出する処理部、
を少なくとも有する、情報分析装置。
（９）微小粒子からの光を検出波長域が異なる複数の受光素子により検出する検出部と、
該検出部から取得された前記光の強度値を、各受光素子の検出波長域幅で補正して、第一の補正強度値を算出する処理部と、
を少なくとも有する、フローサイトメータ。
（１０）微小粒子からの光を検出波長域が異なる複数の受光素子で検出し、前記複数の受光素子により取得された前記光の強度値を、各受光素子の検出波長域幅で補正して、第一の補正強度値を算出する、データ解析方法。
（１１）前記複数の受光素子により取得された光の強度値を、各受光素子が有する感度データに基づき補正処理した補正強度値を用いてスペクトルデータを生成する、（１０）に記載のデータ解析方法。
（１２）前記第一の補正強度値を、各受光素子の感度データを用いて補正して、第二の補正強度値を算出する、（１０）に記載のデータ解析方法。
（１３）一の軸を前記検出波長とし、他の一の軸を前記第一の補正強度値又は前記第二の補正強度値とするスペクトルチャートを生成する、（１２）に記載のデータ解析方法。
（１４）一のパラメータを前記検出波長とし、他の一のパラメータを前記第一の補正強度値又は前記第二の補正強度値とするスペクトルデータを生成し、基準スペクトルデータと比較し、両者が一致しているかあるいは不一致であるかを判定する、（１３）に記載のデータ解析方法。

図２に示した構成の検出部を備えたスペクトル型フローサイトメータを試作した。光源には、波長４８８ｎｍのレーザーダイオードと波長６３８ｎｍのレーザーダイオードを用いた。また、分光素子には、複数のプリズムを組み合わせたプリズムアレイを用いた。受光素子アレイには、３２チャネルのＰＭＴアレイを使用し、波長５００ｎｍから８００ｎｍの蛍光を分光検出した。

試作した装置において、各ＰＭＴの検出波長域幅を決定したグラフを図５に示す。グラフ中の「×」は各チャネルのＰＭＴの検出下限波長（Ｌ［ｋ］）を、「○」は検出上限波長（Ｈ［ｋ］）を示す。ここで、ｋは１〜３２の整数を表す。各ＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）は、チャネル番号が大きい長波長側のＰＭＴほど広くなっていることが確認される。なお、波長６３８ｎｍ付近の蛍光を検出するチャネル２１番前後のＰＭＴでは、波長６３８ｎｍの光源からのレーザー光の漏れ込みを防止する光学フィルタによって、検出される蛍光も制限されている。

また、各ＰＭＴの相対感度を算出したグラフを図６に示す。相対感度は、同一強度及び同一波長の光を各ＰＭＴに照射して各チャネルで得られた強度値を、最も高い強度値が得られたチャネル３２の強度値を１とした相対値により示した。

初めに、F-4500型蛍光分光光度計（株式会社日立ハイテクノロジーズ）を用いて市販の蛍光ビーズの蛍光スペクトルを計測した。蛍光ビーズには、Sherotech社より入手したFluorescent Particle Kit （FPK）505、FPK528、FPK549、FPK667の４種を用いた。得られたスペクトルチャート（基準スペクトルチャート）を図７に示す。（Ａ）はFPK505、（Ｂ）はFPK528、（Ｃ）はFPK549、（Ｄ）はFPK667の蛍光スペクトルを示す。横軸は蛍光波長（５００〜８００ｎｍ）、縦軸は蛍光強度値（対数表示）である。なお、レーザー光の励起波長は、（Ａ）〜（Ｃ）が波長４８８ｎｍ、（Ｄ）が波長６３８ｎｍである。

次に、試作した装置を用いて蛍光ビーズの蛍光スペクトルを計測した。得られたスペクトルチャートを図８〜１２に示す。図８及び図９はFPK505、図１０はFPK528、図１１はFPK549、図１２はFPK667のチャートを示す。なお、図８では各チャネルでのイベント数をスペクトルの色によって表示した。また、図９では強度値をイベント数に基づく平均値（実線）及び平均値±標準偏差（破線）で表示している。

図８〜１２の（Ａ）は、横軸にチャネル番号、縦軸に各チャネルで取得された蛍光の強度値（Ｉ［ｋ］、ｋは１〜３２の整数を表す）の対数をとったスペクトルチャートである。

図８〜１２の（Ａ）のスペクトルチャートに示されるスペクトル形状は、図７に示した基準スペクトルチャートのスペクトル形状と明らかに異なっている。このことは、ＰＭＴで取得される蛍光の強度値（Ｉ［ｋ］）をそのまま用いた蛍光スペクトルでは、装置の光学系の非線形性及び受光素子間の感度差に起因した測定誤差によってスペクトル形状に歪みが生じていることを示している。

図８〜１２の（Ｂ）には、横軸に検出波長、縦軸に各チャネルで取得された蛍光強度値の第一の補正値（Ｊ_１［ｋ］、ｋは１〜３２の整数を表す）の対数をとったスペクトルチャートである。第一の補正強度値Ｊ_１［ｋ］は、各ＰＭＴで取得された蛍光の強度値（Ｉ［ｋ］）を、図５に示したそれぞれのＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）で除すことにより得た。より具体的には、チャネルｋのＰＭＴで得られたｎ番目の強度値Ｉ［ｋ、ｎ］をＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）で除して第一の補正強度値Ｊ_１［ｋ、ｎ］を得て、Ｊ_１［ｋ、ｎ］の分布を横軸Ｌ［ｋ］〜Ｈ［ｋ］の範囲に描画してスペクトルチャートを作成した。

図８〜１２の（Ｂ）のスペクトルチャートに示されるスペクトル形状は、図７に示した基準スペクトルチャートのスペクトル形状と概ね一致している。このことは、各ＰＭＴで取得された蛍光の強度値（Ｉ［ｋ］）をそれぞれのＰＭＴの検出波長域幅（Ｈ［ｋ］−Ｌ［ｋ］）で除す補正処理によって、装置の光学系の非線形性に起因した測定誤差を補償して、スペクトル形状の歪みを補正できたことを示している。

図８〜１２の（Ｃ）には、横軸に検出波長、縦軸に各チャネルで取得された蛍光強度値の第二の補正値（Ｊ_２［ｋ］、ｋは１〜３２の整数を表す）の対数をとったスペクトルチャートである。第二の補正強度値Ｊ_２［ｋ］は、第一の補正強度値（Ｊ_１［ｋ］）を、図６に示したそれぞれのＰＭＴの相対感度（Ｓ［ｋ］）で除すことにより得た。

図８〜１２の（Ｃ）のスペクトルチャートに示されるスペクトル形状は、図７に示した基準スペクトルチャートのスペクトル形状と良く一致している。特に、図８〜１２の（Ｂ）の第一の補正値（Ｊ_１［ｋ］）に基づくスペクトルチャートでは、波長５００ｎｍ付近の領域でＰＭＴの感度差に起因すると思われるスペクトル形状の歪みがみられたが、図８〜１２の（Ｃ）の第二の補正値（Ｊ_２［ｋ］）に基づくスペクトルチャートでは、この歪みが補正されている。このことは、第一の補正値（Ｊ_１［ｋ］）をそれぞれのＰＭＴの相対感度（Ｓ［ｋ］）で除す補正処理によって、受光素子間の感度差に起因した測定誤差を補償して、スペクトル形状の歪みを補正できることが示された。

以上の結果から、本技術に係る微小粒子測定装置によれば、微小粒子について、通常の蛍光分光光度計により計測して得られる基準スペクトルチャートと良く一致するスペクトル形状を計測可能であることが分かる。

このため、本技術に係る微小粒子測定装置では、標識された蛍光色素の種類や量が未知のサンプルについて得られたスペクトル形状を、既知の蛍光物質のスペクトル形状を記録したデータベース（リフェレンスデータ）から検索することにより、サンプルに標識された蛍光色素の種類や量を推定することも可能となる。

また、本技術に係る微小粒子測定装置では、複数の蛍光色素で標識されたサンプルを解析する際のコンペンセーション計算を、従来必要とされていた各蛍光色素で単染色されたサンプルの計測を行うことなく、予め計測した各蛍光色素の蛍光スペクトルを用いて行うことも可能となる。そして、これにより、サンプル解析のための手間や時間、試薬などのリソースを低減することができる。

本技術に係る微小粒子測定システム及び微小粒子測定方法によれば、サンプル解析の都度に標準サンプルを計測しなくとも、測定誤差を補正して正確な解析結果を得ることができる。従って、本技術に係る微小粒子測定システム及び微小粒子測定方法は、細胞等の微小粒子の光学特性をより詳細に解析するため寄与する。

Ａ：微小粒子測定装置、Ｐ：微小粒子、１０：検出部、１０１：光源、１０２：分光素子、１０３：受光素子アレイ、１０４〜１０６：波長選択素子、１０７〜１１０：受光素子、２０：ＣＰＵ、３０：メモリ、４０：ハードディスク（記憶部）、４１：蛍光強度補正プログラム、４２：リファレンスデータ、４３：ＯＳ、５１：マウス、５２：キーボード、６０：表示部、６１：ディスプレイ、６２：プリンタ

Claims (36)

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. 前記照射工程では、前記２種以上のレーザー光を前記流路内の異なる位置に照射する、請求項２９に記載の微小粒子測定方法。
31. 前記照射工程では、４８８ｎｍを有するレーザー光と第一の蛍光波長を有するレーザー光とを照射し、前記第一の蛍光波長は４８８ｎｍ超である、請求項２９又は３０に記載の微小粒子測定方法。
32. 前記検出工程では、４８８ｎｍ以上の波長域の光を検出し、前記第一の蛍光波長のレーザー光の漏れ込みを防止する、請求項３１に記載の微小粒子測定方法。
33. 前記コンペンセーション処理には、最小二乗法が用いられる、請求項１９〜３２のいずれか一項に記載の微小粒子測定方法。
34. 前記基準スペクトルデータは、前記微小粒子が標識された複数の蛍光色素それぞれで単染色されたサンプルにより取得されたデータである、請求項１９〜３３のいずれか一項に記載の微小粒子測定方法。
35. 前記処理工程では、記憶部から前記基準スペクトルデータを受け取る、請求項１９〜３４のいずれか一項に記載の微小粒子測定方法。
36. 前記微小粒子は、細胞である、請求項１９〜３５のいずれか一項に記載の微小粒子測定方法。
    

Images