蛍光強度補正方法及び蛍光強度算出装置
【課題】複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示する技術の提供。
【解決手段】蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法を提供する。
【解決手段】蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法を提供する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、蛍光強度補正方法と蛍光強度算出装置に関する。より詳しくは、微小粒子に多重標識された複数の蛍光色素のそれぞれから発生する蛍光の強度を正確に算出するための蛍光強度補正方法等に関する。
【背景技術】
【０００２】
従来、細胞等の微小粒子を蛍光色素を用いて標識し、これにレーザ光を照射して励起された蛍光色素から発生する蛍光の強度やパターンを計測することによって、微小粒子の特性を測定する装置（例えば、フローサイトメータ）が用いられている。近年では、細胞等の特性をより詳細に分析するため、微小粒子を複数の蛍光色素を用いて標識し、各蛍光色素から発生する光を受光波長帯域の異なる複数の光検出器（ＰＤやＰＭＴなど）により計測するマルチカラー測定が行われるようになっている。マルチカラー測定では、用いる蛍光色素の蛍光波長に応じて光検出器側の光学フィルタを選択して蛍光の検出を行っている。
【０００３】
一方、現在利用されている蛍光色素（例えば、ＦＩＴＣ、ＰＥ（フィコエリスリン）、ＡＰＣ（アロフィコシアニン）など）は、蛍光スペクトルに互いに重複する波長帯域が存在する。従って、これらの蛍光色素を組み合わせてマルチカラー測定を行う場合、各蛍光色素から発生する蛍光を光学フィルタにより波長帯域別に分離しても、各光検出器には目的以外の蛍光色素からの蛍光が漏れ込むことがある。蛍光の漏れ込みが生じると、各光検出器で計測される蛍光強度と目的とする蛍光色素からの真の蛍光強度にずれが生じ、測定誤差の原因となる
【０００４】
この測定誤差を補正するため、光検出器で計測された蛍光強度から漏れ込み分の蛍光強度を差し引く蛍光補正（コンペンセーション）が行われている。蛍光補正は、光検出器で計測された蛍光強度が、目的とする蛍光色素からの真の蛍光強度となるように、パルスに電気的あるいは数学的な補正を加えるものである。
【０００５】
数学的に蛍光補正を行う方法として、各光検出器で計測された蛍光強度（検出値）をベクトルとして表し、このベクトルに予め設定した漏れ込み行列の逆行列を作用させることで、目的とする蛍光色素からの真の蛍光強度を算出する方法が用いられている（図１２・１３、特許文献１参照）。この漏れ込み行列は、各蛍光色素を個別に単標識した微小粒子の蛍光波長分布を解析することによって作成されるものであり、各蛍光色素の蛍光波長分布が列ベクトルとして配列されたものである。漏れ込み行列の逆行列は「補正行列」とも称される。図１２・１３には、５種類（ＦＩＴＣ，ＰＥ，ＥＣＤ，ＰＣ５，ＰＣ７）の蛍光色素と５つの光検出器を用いて５カラー測定を行う場合を例に示した。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００３−８３８９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
補正行列を用いた蛍光補正方法では、補正後の蛍光強度が負値になる場合がある。これは、各光検出器の検出値に含まれるノイズが行列演算に影響を与えていることが原因である。しかし、現実には各蛍光色素からの蛍光強度は負値にはなり得ない。さらに、ある蛍光色素からの蛍光強度が負値として算出されるということは、同時に他の蛍光色素の蛍光強度の算出値に正方向への誤差が生じていることを意味する。
【０００８】
解析する微小粒子集団（ポピュレーション）に、ある蛍光色素についての蛍光強度が負値となる小集団（サブポピュレーション）が存在すると、当該蛍光色素の蛍光強度を対数軸（ログスケール）でプロットした二次元相関図（サイトグラム）上に当該サブポピュレーションがプロットされない。このため、二次元相関図上にプロットされたポピュレーションが実際よりも少なくなったような誤解をユーザに与えるおそれがある。
【０００９】
本発明は、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示するための技術を提供することを主な目的とする。
【課題を解決するための手段】
【００１０】
上記課題解決のため、本発明は、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法を提供する。
この蛍光強度補正方法は、前記手順においては、下記式（２）を満足しながら下記式（１）で示される評価関数が最小値となるパラメータｘｊ（ｊ＝１〜Ｍ）を求めることにより、各蛍光色素から発生した蛍光の強度を算出する。式（２）の制約条件を設けて補正演算を行うことにより、蛍光色素からの蛍光強度が負値として算出されることに起因する測定誤差や二次元相関図（サイトグラム）上のポピュレーション減少の問題を解決できる。
【数１】
（式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。）
前記パラメータｘｊ（ｊ＝１〜Ｍ）は、下記式（３）〜（５）の二次計画問題を解くことにより求めることができる。
【数２】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ただし、ＡはＭ×Ｍ次行列、ｂはＭ×１次行列として設定し、Ｕ１〜ＵＭは前記下限値を表すものとする。）
前記式（２）中、前記下限値Ｕｊとしては、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊ、該検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数、あるいは前記検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求めて下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数、のいずれかを用いることができる。
【数３】
【００１１】
また、前記パラメータｘｊ（ｊ＝１〜Ｍ）は、上記式（４）の線形制約条件がない、下記式（１０），（１１）の二次計画問題を解くことにより求めることもできる。
【数４】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。）
さらに、上記式（４）の線形制約条件がない二次計画問題を解く場合には、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めてもよい。
【数５】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。）
あるいは、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めてもよい。
【数６】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
さらに、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【数７】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【００１２】
また、本発明は、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段と、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により前記測定スペクトルを近似する算出手段と、を備え、前記算出手段は、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法により行う蛍光強度算出装置をも提供する。
【００１３】
本発明において、「微小粒子」には、細胞や微生物、リポソームなどの生体関連微小粒子、あるいはラテックス粒子やゲル粒子、工業用粒子などの合成粒子などが広く含まれるものとする。
生体関連微小粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラ(細胞小器官)などが含まれる。細胞には、動物細胞(血球系細胞など)および植物細胞が含まれる。微生物には、大腸菌などの細菌類、タバコモザイクウイルスなどのウイルス類、イースト菌などの菌類などが含まれる。さらに、生体関連微小粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。また、工業用粒子は、例えば有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレン・ジビニルベンゼン、ポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカ、磁性体材料などが含まれる。金属には、金コロイド、アルミなどが含まれる。これら微小粒子の形状は、一般には球形であるのが普通であるが、非球形であってもよく、また大きさや質量なども特に限定されない。
【発明の効果】
【００１４】
本発明により、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示する技術が提供される。
【図面の簡単な説明】
【００１５】
【図１】測定スペクトルを単染色スペクトルの線形和により近似して得られる近似曲線を説明するグラフである。
【図２】実施例で用いた各蛍光色素の単染色スペクトルを示すグラフである。
【図３】蛍光補正を行わずに、仮想出力データから作成した二次元相関図である。
【図４】仮想出力データを、補正行列を用いて蛍光補正して作成した二次元相関図である（比較例）。
【図５】仮想出力データを、下限値を平均値とする第一実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図６】仮想出力データを、線形制約条件なしで、下限値を平均値とする第一実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図７】仮想出力データを、下限値を確率密度関数に従う乱数とする第二実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図８】仮想出力データを、線形制約条件なしで、下限値を確率密度関数に従う乱数とする第二実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図９】仮想出力データを、下限値を正規乱数とする第三実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図１０】仮想出力データを、線形制約条件なしで、下限値を正規乱数とする第三実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図１１】仮想出力データを、下限値を０とする第四実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図１２】従来の補正行列を用いた蛍光補正方法を説明する図である。
【図１３】従来の補正行列の行列要素を説明する図である。
【発明を実施するための形態】
【００１６】
以下、本発明を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本発明の代表的な実施形態の一例を示したものであり、これにより本発明の範囲が狭く解釈されることはない。説明は以下の順序で行う。
１．蛍光強度補正方法
（１）測定手順
（２）算出手順
（２−１）近似曲線
（２−２）制限付き最小二乗法
（２−３）制約条件（下限値条件）
（２−３−１）平均値（第一実施形態とその変形例）
（２−３−２）確率密度関数に従う乱数（第二実施形態とその変形例）
（２−３−３）正規分布に従う乱数（第三実施形態とその変形例）
（２−３−４）０（第四実施形態）
２．蛍光強度算出装置
【００１７】
１．蛍光強度補正方法
本発明に係る蛍光強度補正方法は、以下の２つの手順を含む。
測定手順：蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る手順。
算出手順：制限付き最小二乗法を用いて、測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順。
【００１８】
（１）測定手順
まず、測定対象とする微小粒子を複数の蛍光色素を用いて多重標識する。微小粒子の蛍光色素標識は従来公知の手法によって行うことができる。例えば測定対象を細胞とする場合には、細胞表面分子に対する蛍光標識抗体と細胞とを混合し、細胞表面分子に抗体を結合させる。蛍光標識抗体は、抗体に直接蛍光色素を結合させたものであってよく、ビオチン標識した抗体にアビジンを結合した蛍光色素をアビジン・ビオチン反応によって結合させたものであってもよい。また、抗体は、モノクローナル抗体又はポリクローナル抗体であってよい。
【００１９】
蛍光色素には、従来公知の物質を２以上組み合わせて用いることができる。例えば、フィコエリスリン（ＰＥ）、ＦＩＴＣ、ＰＥ−Ｃｙ５、ＰＥ−Ｃｙ７、ＰＥ−テキサスレッド(PE-Texas red)、アロフィコシアニン（ＡＰＣ）、ＡＰＣ−Ｃｙ７、エチジウムブロマイド（Ethidium bromide）、プロピジウムアイオダイド（Propidium iodide）、ヘキスト（Hoechst）33258/33342、ＤＡＰＩ、アクリジンオレンジ（Acridine orange）、クロモマイシン（Chromomycin）、ミトラマイシン（Mithramycin）、オリボマイシン（Olivomycin）、パイロニン（Pyronin）Ｙ、チアゾールオレンジ（Thiazole orange）、ローダミン（Rhodamine）１０１イソチオシアネート（isothiocyanate）、ＢＣＥＣＦ、ＢＣＥＣＦ−ＡＭ、Ｃ．ＳＮＡＲＦ−１、Ｃ．ＳＮＡＲＦ−１−ＡＭＡ、エクオリン（Aequorin）、Ｉｎｄｏ−１、Ｉｎｄｏ−１−ＡＭ、Ｆｌｕｏ−３、Ｆｌｕｏ−３−ＡＭ、Ｆｕｒａ−２、Ｆｕｒａ−２−ＡＭ、オキソノール（Oxonol）、テキサスレッド(Texas red)、ローダミン（Rhodamine）１２３、１０−Ｎ−ノニ−アクリジンオレンジ（Acridine orange）、フルオレセイン（Fluorecein）、フルオレセインジアセテート（Fluorescein diacetate）、カルボキシフルオレセイン（Carboxyfluorescein）、カルビキシフルオレセインジアセテート（Caboxyfluorescein diacetate）、カルボキシジクロロフルオレセイン（Carboxydichlorofluorescein）、カルボキシジクロロフルオレセインジアセテート（Carboxydichlorofluorescein diacetate）が挙げられる。
【００２０】
次に、複数の蛍光色素により多重標識された微小粒子に対して光を照射し、励起された蛍光色素から発生する蛍光を受光波長帯域の異なる複数の光検出器で受光する。測定手順は、従来公知のマルチカラー測定フローサイトメータを用いた方法と同様にして行うことができる。
【００２１】
（２）算出手順
（２−１）近似曲線
算出手順においては、測定手順で取得された各光検出器の検出値を補正演算して各蛍光色素からの蛍光強度を算出する。この際、本発明に係る蛍光強度補正方法では、制限付き最小二乗法を用いて、測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似することにより、各蛍光色素からの真の蛍光強度を算出する。
【００２２】
ここで、「測定スペクトル」とは、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られるものである。また、「単染色スペクトル」とは、各蛍光色素の蛍光波長分布であり、各蛍光色素を個別に標識した微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を光検出器で受光し、検出値を収集して得られるものである。
【００２３】
図１に基づいて、測定スペクトルを単染色スペクトルの線形和により近似して得られる近似曲線について説明する。
【００２４】
図中、Ｘ軸は観測点を、Ｙ軸は検出値を示す。図では、光検出器１で受光された蛍光の検出値をｐ１、光検出器２で受光された蛍光の検出値をｐ２、光検出器ｎで受光された蛍光の検出値をｐｎで示している。各検出値ｐ１〜ｐｎを結ぶ線が測定スペクトルである。
【００２５】
また、図では、１番目の蛍光色素（蛍光色素１）の単染色スペクトルを表す曲線（基底関数）をＳ１（ｉ），２番目の蛍光色素（蛍光色素２）の単染色スペクトルを表す曲線をＳ２（ｉ），Ｍ番目の蛍光色素（蛍光色素Ｍ）の単染色スペクトルを表す曲線をＳＭ（ｉ）で示している。単染色スペクトルは、測定の都度に各蛍光色素を個別に標識したサンプルを調製して取得してもよく、あるいは予め装置に記憶された標準スペクトルを利用してもよい。
【００２６】
各光検出器では、蛍光色素１から蛍光色素Ｍまでの全ての蛍光色素からの蛍光がそれぞれ所定比率で漏れ込んだ状態で受光される。そのため、各光検出器の検出値ｐｉは、蛍光色素１から蛍光色素Ｍまでの基底関数にそれぞれ所定比率を乗じた値の和として下記式ｙ（ｉ）によって近似することができる。ここで、各光検出器への各蛍光色素からの蛍光の漏れ込み比率ｘｊは、各蛍光色素の発光強度（真の蛍光強度）により規定される。
【００２７】
【数８】
【００２８】
具体的には、例えば、光検出器１の検出値ｐ１は、蛍光色素１の蛍光強度Ｓ１（１）に比率ｘ１を乗じた値から蛍光色素Ｍの基底関数ＳＭ（１）に比率ｘＭを乗じた値までの和ｙ（１）として近似される。そして、光検出器１への蛍光色素１〜Ｍからの蛍光の漏れ込み比率ｘｊ（ｊ＝１〜Ｍ）は、蛍光色素１〜Ｍの発光強度に相応する。
【００２９】
この式で示される近似曲線は、次に説明する制限付き最小二乗法を用いてｘｊを求めることにより得られる。このｘｊは、各蛍光色素の真の蛍光強度に同等であり、物理的に負値にはなり得ない値である。そこで、本発明に係る蛍光強度補正方法では、このｘｊについて、所定値よりも大きいという制約条件（下限値条件）を設けて補正演算を行う。下限値条件を設けて補正演算を行うことにより、蛍光色素からの蛍光強度が負値として算出されることに起因する測定誤差や二次元相関図（サイトグラム）上のポピュレーション減少の問題を解決できる。
【００３０】
（２−２）制限付き最小二乗法
以下に、ｘｊを求めるための手順を説明する。まず、下記式（１）で示される評価関数（カイ二乗）を定義する。そして、この評価関数が、下記式（２）を満足しながら、最小値となるようなパラメータｘｊ（ｊ＝１〜Ｍ）を求める。
【００３１】
【数９】
【００３２】
式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値Ｓｊ（ｉ）を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。
【００３３】
【数１０】
【００３４】
式中、Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。
【００３５】
上記式（２）は、蛍光強度ｘｊが所定値（Ｕｊ）よりも大きいという制約条件（下限値条件）を表す。
【００３６】
ここで、ｓｉｊを要素とするＮ×Ｍ次行列Ｓと、ｘｊを要素とするＭ次行列ｘと、ｐｉを要素とするＮ次ベクトルｐを設定すると、式（２）を満足しながら式（１）で示される評価関数が最小値となるパラメータｘｊを求めることは、次の問題を解くこと同じである。
【００３７】
【数１１】
【００３８】
Ｍ×Ｍ次行列ＡとＭ×１次行列ｂは、蛍光強度ｘが所定値（Ｕ）以上であるように制約する不等式として以下の式（６），（７）のように設定される。
【００３９】
【数１２】
【００４０】
上記式（３）は２乗して展開すると、下記式（２３）のようになる。
【００４１】
【数１３】
【００４２】
式（２３）を最小化するため、最終項のｐＴｐは無視できる。従って、式（２３）を最小化することは、下記式（２４）を最小化することに等しい。
【００４３】
【数１４】
【００４４】
式（２４）及び上記式（４），（５）は、二次計画問題として知られている。「二次計画問題（Quadratic Programming Problem）」は、Ｄをｎ×ｎ次非負定値対称行列、ｃをｎ次ベクトル、Ａをｍ×ｎ次行列、ｂをｍ次ベクトルとするとき、線形制約条「Ａｘ≦ｂ、ｘ≧０」を満たすｎ次ベクトルのうち、２次形式で与えられる目的関数「f（ｘ）=ｘTAｘ／２＋ｃＴｘ」を最小にするものを求める問題である。二次計画問題は、次の式（２５）〜（２７）で表される。二次計画問題は、有限回の反復で厳密な最適解が得られる特殊な非線形計画問題として知られている。
【００４５】
【数１５】
【００４６】
上記式（２４），（４），（５）は、式（２５）中のＤ，ｃを以下のように置き換え、式（２６）のＡ，ｂを式（６），（７）のように置き換えて二次計画問題を解くことに等しい。
【数１６】
【００４７】
（２−３）制約条件（下限値条件）
（２−３−１）平均値（第一実施形態とその変形例）
下記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、蛍光色素を標識していない無染色の微小粒子を測定したとして、ｊ番目の蛍光色素にとって最小の検出値として適当な値を設定する。具体的には、下限値Ｕｊは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊを用いることができる。無染色平均値Ｖｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値が用いられる。あるいは、無染色平均値Ｖｊは、最大の検出値をとる検出器の検出平均値とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の検出平均値のさらに平均をとった値などとすることもできる。
【００４８】
【数１７】
【００４９】
【数１８】
【００５０】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用できる。この場合、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００５１】
【数１９】
【００５２】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。無染色平均値Ｖｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値が用いられる。あるいは、無染色平均値Ｖｊは、最大の検出値をとる検出器の検出平均値とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の検出平均値のさらに平均をとった値などとすることもできる。
【００５３】
（２−３−２）確率密度関数に従う乱数（第二実施形態とその変形例）
また、上記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された光検出器毎の検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いることもできる。無染色確率密度関数Ｆｊ（Ｕｊ）としては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の確率密度関数が用いられる。あるいは、無染色確率密度関数Ｆｊ（Ｕｊ）は、最大の検出値をとる検出器の確率密度関数とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の確率密度関数とから算出した関数などとすることもできる。
【００５４】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用する場合、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００５５】
【数２０】
【００５６】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。無染色確率密度関数Ｆｊ（Ｕｊ）としては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の確率密度関数ｆｉ（ｕｉ）が用いられる。あるいは、無染色確率密度関数Ｆｊ（Ｕｊ）は、最大の検出値をとる検出器の確率密度関数とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の確率密度関数とから算出した関数などとすることもできる。
【００５７】
（２−３−３）正規分布に従う乱数（第三実施形態とその変形例）
さらに、上記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、前記下限値Ｕｊとして、下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いてもよい。この乱数は、正規分布に従った乱数（正規乱数）となる。無染色平均値Ｖｊ及び無染色標準偏差ρｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値及び標準偏差が用いられる。あるいは、無染色平均値Ｖｊ及び無染色標準偏差ρｊは、最大の検出値をとる検出器の平均値及び標準偏差とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の平均値及び標準偏差とから算出した値などとすることもできる。
【００５８】
【数２１】
【００５９】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用する場合、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００６０】
【数２２】
【００６１】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。無染色平均値Ｖｊ及び無染色標準偏差ρｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値ｉ及び標準偏差が用いられる。あるいは、無染色平均値Ｖｊ及び無染色標準偏差ρｊは、最大の検出値をとる検出器の平均値及び標準偏差とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の平均値及び標準偏差とから算出した値などとすることもできる。
【００６２】
（２−３−４）０（第四実施形態）
上記式（４）の線形制約条件がない二次計画法を適用する場合、下限値を０とすれば、下記式（１０），（１１）の二次計画問題を解くことにより、パラメータｘｊ（ｊ＝１〜Ｍ）を求めることができる。
【００６３】
【数２３】
【００６４】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。
【００６５】
２．蛍光強度算出装置
本発明に係る蛍光強度算出装置は、従来のフローサイトメータ等と同様に流体系と光学系、分取系、データ処理系などから構成される。
【００６６】
流体系は、フローセルにおいて測定対象とする微小粒子を含むサンプル液をシース液の層流の中心に流し、フローセル内に微小粒子を一列に配列させる手段である。フローセルに替えて、マイクロチップ上に形成した流路内において、微小粒子を一列に配列させてもよい。
【００６７】
光学系は、蛍光色素により標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段である。光学系では、微小粒子から発生する前方散乱光や側方散乱光、レイリー散乱やミー散乱等の散乱光も検出される。光学系は、具体的には、レーザ光源と、微小粒子に対してレーザ光を集光・照射するための集光レンズやダイクロイックミラー、バンドパスフィルター等からなる照射系と、レーザ光の照射によって微小粒子から発生する蛍光や散乱光を検出する検出系と、によって構成される。検出系は、例えば、PMT（photo multiplier tube）や、CCDやCMOS素子等のエリア撮像素子等によって構成され、受光波長帯域の異なる光検出器が複数配設される。
【００６８】
微小粒子の分取を行う場合には、サンプル液を、微小粒子を含む液滴としてフローセル外の空間に吐出し、液滴の移動方向を制御して所望の特性を備えた微小粒子を分取する。分取系は、サンプル液を液滴化してフローセルから吐出させるピエゾ素子等の振動素子と、吐出される液滴に電荷を付与する荷電手段と、液滴の移動方向に沿って、移動する液滴を挟んで対向して配設された対電極などから構成される。
【００６９】
データ処理系は、光検出器から検出値を電気信号として入力され、電気信号に基づいて微小粒子の光学特性を解析する。データ処理系は、各光検出器から検出値を収集して得られる測定スペクトルを、上述した方法に従って単染色スペクトルの線形和により近似し、各蛍光色素から発生した真の蛍光強度を算出する。このため、データ処理系は、上述した本発明に係る蛍光強度算出方法の各ステップを実行するためのプログラムを記憶するハードディスク等の記録媒体と、プラグラムを実行するCPU及びメモリ等を有する。
【００７０】
データ処理系には、各光検出器からの信号の品質を向上させるため、各光検出器間のデータの平滑化を行うノイズフィルタを設けることが好ましい。平均化処理は、光検出器の数や受光波長帯域、使用する蛍光色素や、測定対象微小粒子の自家蛍光波長、さらに各機器のノイズ周波数成分等を考慮して行われる。データの平滑化を行うことにより、ノイズの影響を抑えることが可能となり、演算精度が向上し、より正確に蛍光強度を算出できる。
【実施例】
【００７１】
３２チャンネルのマルチカラー測定フローサイトメータでの仮想出力データを用いて、従来の補正行列を用いた蛍光補正方法と本発明に係る制限付き最小二乗法を用いた蛍光補正方法との間で演算結果の比較を行った。
【００７２】
蛍光色素には、ＦＩＴＣ（ＣＨ５），Ａｌｅｘａ５００（ＣＨ６），Ａｌｅｘａ５１４（ＣＨ１０），Ａｌｅｘａ５３２（ＣＨ１２），ＰＥ（ＣＨ１５），ＰＥ−ＴＲ（ＣＨ１９），ＰＩ（ＣＨ２１），Ａｌｅｘａ６００（ＣＨ２０），ＰＥ−Ｃｙ５（ＣＨ２４），ＰｅｒＣＰ（ＣＨ２５），ＰｅｒＣＰ−Ｃｙ５．５（ＣＨ２６），ＰＥ−Ｃｙ７（ＣＨ３１）を用いた。括弧内に示すチャネル番号は、補正行列を用いた従来方法における各蛍光色素に対するチャネルの割り当てを示す。
【００７３】
図２に、各蛍光色素の単染色スペクトルを示す。これらの各単染色スペクトルがランダムな強度で排他的に出現し、さらに各検出データには検出器や電気的なノイズが仮想的に付加されたシミュレーション・データをモンテカルロ法にて生成させた。発生させたシミュレーション・データを、図２の単染色スペクトルと、従来方法あるいは本発明に係る方法とを用いて蛍光補正を行って得た典型的な解析結果を以下に示す。
【００７４】
図３は、蛍光補正を行わずに、仮想出力データから作成した二次元相関図を示す。（Ａ）〜（Ｃ）の縦軸はそれぞれＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度を、横軸はＦＩＴＣの蛍光強度を示す。図上段はログスケールプロットした二次元相関図、下段はリニアスケールでプロットした二次元相関図である（以下、図４〜１１について同様）。
【００７５】
図４は、仮想出力データを、比較のために補正行列を用いて蛍光補正して作成した二次元相関図を示す。リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が負値として算出されているサブポピュレーションが存在している。このため、ログスケールの二次元相関図上にプロットされたポピュレーションが、リニアスケールの二次元相関図上にプロットされたポピュレーションに対して顕著に減少してしまっている。
【００７６】
図５は、仮想出力データを、下限値を平均値とする第一実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が正値として算出されている。このため、リニアスケールの二次元相関図上にプロットされたポピュレーションが、全てログスケールの二次元相関図上にもプロットできている。
【００７７】
図６は、仮想出力データを、線形制約条件なしで、下限値を平均値とする第一実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００７８】
図７は、仮想出力データを、下限値を確率密度関数に従う乱数とする第二実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００７９】
図８は、仮想出力データを、線形制約条件なしで、下限値を確率密度関数に従う乱数とする第二実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８０】
図９は、仮想出力データを、下限値を正規乱数とする第三実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８１】
図１０は、仮想出力データを、線形制約条件なしで、下限値を正規乱数とする第三実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８２】
図１１は、仮想出力データを、下限値を０とする第四実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８３】
図６〜１１においても、リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が正値として算出されている。このため、リニアスケールの二次元相関図上にプロットされたポピュレーションが、全てログスケールの二次元相関図上にもプロットできている。
【産業上の利用可能性】
【００８４】
本発明に係る蛍光強度補正方法等によれば、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示できる。従って、本発明に係る蛍光強度補正方法等は、細胞等の微小粒子の特性をより詳細に解析するため寄与し得る。
【技術分野】
【０００１】
本発明は、蛍光強度補正方法と蛍光強度算出装置に関する。より詳しくは、微小粒子に多重標識された複数の蛍光色素のそれぞれから発生する蛍光の強度を正確に算出するための蛍光強度補正方法等に関する。
【背景技術】
【０００２】
従来、細胞等の微小粒子を蛍光色素を用いて標識し、これにレーザ光を照射して励起された蛍光色素から発生する蛍光の強度やパターンを計測することによって、微小粒子の特性を測定する装置（例えば、フローサイトメータ）が用いられている。近年では、細胞等の特性をより詳細に分析するため、微小粒子を複数の蛍光色素を用いて標識し、各蛍光色素から発生する光を受光波長帯域の異なる複数の光検出器（ＰＤやＰＭＴなど）により計測するマルチカラー測定が行われるようになっている。マルチカラー測定では、用いる蛍光色素の蛍光波長に応じて光検出器側の光学フィルタを選択して蛍光の検出を行っている。
【０００３】
一方、現在利用されている蛍光色素（例えば、ＦＩＴＣ、ＰＥ（フィコエリスリン）、ＡＰＣ（アロフィコシアニン）など）は、蛍光スペクトルに互いに重複する波長帯域が存在する。従って、これらの蛍光色素を組み合わせてマルチカラー測定を行う場合、各蛍光色素から発生する蛍光を光学フィルタにより波長帯域別に分離しても、各光検出器には目的以外の蛍光色素からの蛍光が漏れ込むことがある。蛍光の漏れ込みが生じると、各光検出器で計測される蛍光強度と目的とする蛍光色素からの真の蛍光強度にずれが生じ、測定誤差の原因となる
【０００４】
この測定誤差を補正するため、光検出器で計測された蛍光強度から漏れ込み分の蛍光強度を差し引く蛍光補正（コンペンセーション）が行われている。蛍光補正は、光検出器で計測された蛍光強度が、目的とする蛍光色素からの真の蛍光強度となるように、パルスに電気的あるいは数学的な補正を加えるものである。
【０００５】
数学的に蛍光補正を行う方法として、各光検出器で計測された蛍光強度（検出値）をベクトルとして表し、このベクトルに予め設定した漏れ込み行列の逆行列を作用させることで、目的とする蛍光色素からの真の蛍光強度を算出する方法が用いられている（図１２・１３、特許文献１参照）。この漏れ込み行列は、各蛍光色素を個別に単標識した微小粒子の蛍光波長分布を解析することによって作成されるものであり、各蛍光色素の蛍光波長分布が列ベクトルとして配列されたものである。漏れ込み行列の逆行列は「補正行列」とも称される。図１２・１３には、５種類（ＦＩＴＣ，ＰＥ，ＥＣＤ，ＰＣ５，ＰＣ７）の蛍光色素と５つの光検出器を用いて５カラー測定を行う場合を例に示した。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００３−８３８９４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
補正行列を用いた蛍光補正方法では、補正後の蛍光強度が負値になる場合がある。これは、各光検出器の検出値に含まれるノイズが行列演算に影響を与えていることが原因である。しかし、現実には各蛍光色素からの蛍光強度は負値にはなり得ない。さらに、ある蛍光色素からの蛍光強度が負値として算出されるということは、同時に他の蛍光色素の蛍光強度の算出値に正方向への誤差が生じていることを意味する。
【０００８】
解析する微小粒子集団（ポピュレーション）に、ある蛍光色素についての蛍光強度が負値となる小集団（サブポピュレーション）が存在すると、当該蛍光色素の蛍光強度を対数軸（ログスケール）でプロットした二次元相関図（サイトグラム）上に当該サブポピュレーションがプロットされない。このため、二次元相関図上にプロットされたポピュレーションが実際よりも少なくなったような誤解をユーザに与えるおそれがある。
【０００９】
本発明は、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示するための技術を提供することを主な目的とする。
【課題を解決するための手段】
【００１０】
上記課題解決のため、本発明は、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法を提供する。
この蛍光強度補正方法は、前記手順においては、下記式（２）を満足しながら下記式（１）で示される評価関数が最小値となるパラメータｘｊ（ｊ＝１〜Ｍ）を求めることにより、各蛍光色素から発生した蛍光の強度を算出する。式（２）の制約条件を設けて補正演算を行うことにより、蛍光色素からの蛍光強度が負値として算出されることに起因する測定誤差や二次元相関図（サイトグラム）上のポピュレーション減少の問題を解決できる。
【数１】
（式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。）
前記パラメータｘｊ（ｊ＝１〜Ｍ）は、下記式（３）〜（５）の二次計画問題を解くことにより求めることができる。
【数２】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ただし、ＡはＭ×Ｍ次行列、ｂはＭ×１次行列として設定し、Ｕ１〜ＵＭは前記下限値を表すものとする。）
前記式（２）中、前記下限値Ｕｊとしては、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊ、該検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数、あるいは前記検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求めて下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数、のいずれかを用いることができる。
【数３】
【００１１】
また、前記パラメータｘｊ（ｊ＝１〜Ｍ）は、上記式（４）の線形制約条件がない、下記式（１０），（１１）の二次計画問題を解くことにより求めることもできる。
【数４】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。）
さらに、上記式（４）の線形制約条件がない二次計画問題を解く場合には、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めてもよい。
【数５】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。）
あるいは、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めてもよい。
【数６】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
さらに、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【数７】
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【００１２】
また、本発明は、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段と、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により前記測定スペクトルを近似する算出手段と、を備え、前記算出手段は、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法により行う蛍光強度算出装置をも提供する。
【００１３】
本発明において、「微小粒子」には、細胞や微生物、リポソームなどの生体関連微小粒子、あるいはラテックス粒子やゲル粒子、工業用粒子などの合成粒子などが広く含まれるものとする。
生体関連微小粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラ(細胞小器官)などが含まれる。細胞には、動物細胞(血球系細胞など)および植物細胞が含まれる。微生物には、大腸菌などの細菌類、タバコモザイクウイルスなどのウイルス類、イースト菌などの菌類などが含まれる。さらに、生体関連微小粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。また、工業用粒子は、例えば有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレン・ジビニルベンゼン、ポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカ、磁性体材料などが含まれる。金属には、金コロイド、アルミなどが含まれる。これら微小粒子の形状は、一般には球形であるのが普通であるが、非球形であってもよく、また大きさや質量なども特に限定されない。
【発明の効果】
【００１４】
本発明により、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示する技術が提供される。
【図面の簡単な説明】
【００１５】
【図１】測定スペクトルを単染色スペクトルの線形和により近似して得られる近似曲線を説明するグラフである。
【図２】実施例で用いた各蛍光色素の単染色スペクトルを示すグラフである。
【図３】蛍光補正を行わずに、仮想出力データから作成した二次元相関図である。
【図４】仮想出力データを、補正行列を用いて蛍光補正して作成した二次元相関図である（比較例）。
【図５】仮想出力データを、下限値を平均値とする第一実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図６】仮想出力データを、線形制約条件なしで、下限値を平均値とする第一実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図７】仮想出力データを、下限値を確率密度関数に従う乱数とする第二実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図８】仮想出力データを、線形制約条件なしで、下限値を確率密度関数に従う乱数とする第二実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図９】仮想出力データを、下限値を正規乱数とする第三実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図１０】仮想出力データを、線形制約条件なしで、下限値を正規乱数とする第三実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図である。
【図１１】仮想出力データを、下限値を０とする第四実施形態に係る方法により蛍光補正して作成した二次元相関図である。
【図１２】従来の補正行列を用いた蛍光補正方法を説明する図である。
【図１３】従来の補正行列の行列要素を説明する図である。
【発明を実施するための形態】
【００１６】
以下、本発明を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本発明の代表的な実施形態の一例を示したものであり、これにより本発明の範囲が狭く解釈されることはない。説明は以下の順序で行う。
１．蛍光強度補正方法
（１）測定手順
（２）算出手順
（２−１）近似曲線
（２−２）制限付き最小二乗法
（２−３）制約条件（下限値条件）
（２−３−１）平均値（第一実施形態とその変形例）
（２−３−２）確率密度関数に従う乱数（第二実施形態とその変形例）
（２−３−３）正規分布に従う乱数（第三実施形態とその変形例）
（２−３−４）０（第四実施形態）
２．蛍光強度算出装置
【００１７】
１．蛍光強度補正方法
本発明に係る蛍光強度補正方法は、以下の２つの手順を含む。
測定手順：蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る手順。
算出手順：制限付き最小二乗法を用いて、測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順。
【００１８】
（１）測定手順
まず、測定対象とする微小粒子を複数の蛍光色素を用いて多重標識する。微小粒子の蛍光色素標識は従来公知の手法によって行うことができる。例えば測定対象を細胞とする場合には、細胞表面分子に対する蛍光標識抗体と細胞とを混合し、細胞表面分子に抗体を結合させる。蛍光標識抗体は、抗体に直接蛍光色素を結合させたものであってよく、ビオチン標識した抗体にアビジンを結合した蛍光色素をアビジン・ビオチン反応によって結合させたものであってもよい。また、抗体は、モノクローナル抗体又はポリクローナル抗体であってよい。
【００１９】
蛍光色素には、従来公知の物質を２以上組み合わせて用いることができる。例えば、フィコエリスリン（ＰＥ）、ＦＩＴＣ、ＰＥ−Ｃｙ５、ＰＥ−Ｃｙ７、ＰＥ−テキサスレッド(PE-Texas red)、アロフィコシアニン（ＡＰＣ）、ＡＰＣ−Ｃｙ７、エチジウムブロマイド（Ethidium bromide）、プロピジウムアイオダイド（Propidium iodide）、ヘキスト（Hoechst）33258/33342、ＤＡＰＩ、アクリジンオレンジ（Acridine orange）、クロモマイシン（Chromomycin）、ミトラマイシン（Mithramycin）、オリボマイシン（Olivomycin）、パイロニン（Pyronin）Ｙ、チアゾールオレンジ（Thiazole orange）、ローダミン（Rhodamine）１０１イソチオシアネート（isothiocyanate）、ＢＣＥＣＦ、ＢＣＥＣＦ−ＡＭ、Ｃ．ＳＮＡＲＦ−１、Ｃ．ＳＮＡＲＦ−１−ＡＭＡ、エクオリン（Aequorin）、Ｉｎｄｏ−１、Ｉｎｄｏ−１−ＡＭ、Ｆｌｕｏ−３、Ｆｌｕｏ−３−ＡＭ、Ｆｕｒａ−２、Ｆｕｒａ−２−ＡＭ、オキソノール（Oxonol）、テキサスレッド(Texas red)、ローダミン（Rhodamine）１２３、１０−Ｎ−ノニ−アクリジンオレンジ（Acridine orange）、フルオレセイン（Fluorecein）、フルオレセインジアセテート（Fluorescein diacetate）、カルボキシフルオレセイン（Carboxyfluorescein）、カルビキシフルオレセインジアセテート（Caboxyfluorescein diacetate）、カルボキシジクロロフルオレセイン（Carboxydichlorofluorescein）、カルボキシジクロロフルオレセインジアセテート（Carboxydichlorofluorescein diacetate）が挙げられる。
【００２０】
次に、複数の蛍光色素により多重標識された微小粒子に対して光を照射し、励起された蛍光色素から発生する蛍光を受光波長帯域の異なる複数の光検出器で受光する。測定手順は、従来公知のマルチカラー測定フローサイトメータを用いた方法と同様にして行うことができる。
【００２１】
（２）算出手順
（２−１）近似曲線
算出手順においては、測定手順で取得された各光検出器の検出値を補正演算して各蛍光色素からの蛍光強度を算出する。この際、本発明に係る蛍光強度補正方法では、制限付き最小二乗法を用いて、測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似することにより、各蛍光色素からの真の蛍光強度を算出する。
【００２２】
ここで、「測定スペクトル」とは、蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られるものである。また、「単染色スペクトル」とは、各蛍光色素の蛍光波長分布であり、各蛍光色素を個別に標識した微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を光検出器で受光し、検出値を収集して得られるものである。
【００２３】
図１に基づいて、測定スペクトルを単染色スペクトルの線形和により近似して得られる近似曲線について説明する。
【００２４】
図中、Ｘ軸は観測点を、Ｙ軸は検出値を示す。図では、光検出器１で受光された蛍光の検出値をｐ１、光検出器２で受光された蛍光の検出値をｐ２、光検出器ｎで受光された蛍光の検出値をｐｎで示している。各検出値ｐ１〜ｐｎを結ぶ線が測定スペクトルである。
【００２５】
また、図では、１番目の蛍光色素（蛍光色素１）の単染色スペクトルを表す曲線（基底関数）をＳ１（ｉ），２番目の蛍光色素（蛍光色素２）の単染色スペクトルを表す曲線をＳ２（ｉ），Ｍ番目の蛍光色素（蛍光色素Ｍ）の単染色スペクトルを表す曲線をＳＭ（ｉ）で示している。単染色スペクトルは、測定の都度に各蛍光色素を個別に標識したサンプルを調製して取得してもよく、あるいは予め装置に記憶された標準スペクトルを利用してもよい。
【００２６】
各光検出器では、蛍光色素１から蛍光色素Ｍまでの全ての蛍光色素からの蛍光がそれぞれ所定比率で漏れ込んだ状態で受光される。そのため、各光検出器の検出値ｐｉは、蛍光色素１から蛍光色素Ｍまでの基底関数にそれぞれ所定比率を乗じた値の和として下記式ｙ（ｉ）によって近似することができる。ここで、各光検出器への各蛍光色素からの蛍光の漏れ込み比率ｘｊは、各蛍光色素の発光強度（真の蛍光強度）により規定される。
【００２７】
【数８】
【００２８】
具体的には、例えば、光検出器１の検出値ｐ１は、蛍光色素１の蛍光強度Ｓ１（１）に比率ｘ１を乗じた値から蛍光色素Ｍの基底関数ＳＭ（１）に比率ｘＭを乗じた値までの和ｙ（１）として近似される。そして、光検出器１への蛍光色素１〜Ｍからの蛍光の漏れ込み比率ｘｊ（ｊ＝１〜Ｍ）は、蛍光色素１〜Ｍの発光強度に相応する。
【００２９】
この式で示される近似曲線は、次に説明する制限付き最小二乗法を用いてｘｊを求めることにより得られる。このｘｊは、各蛍光色素の真の蛍光強度に同等であり、物理的に負値にはなり得ない値である。そこで、本発明に係る蛍光強度補正方法では、このｘｊについて、所定値よりも大きいという制約条件（下限値条件）を設けて補正演算を行う。下限値条件を設けて補正演算を行うことにより、蛍光色素からの蛍光強度が負値として算出されることに起因する測定誤差や二次元相関図（サイトグラム）上のポピュレーション減少の問題を解決できる。
【００３０】
（２−２）制限付き最小二乗法
以下に、ｘｊを求めるための手順を説明する。まず、下記式（１）で示される評価関数（カイ二乗）を定義する。そして、この評価関数が、下記式（２）を満足しながら、最小値となるようなパラメータｘｊ（ｊ＝１〜Ｍ）を求める。
【００３１】
【数９】
【００３２】
式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値Ｓｊ（ｉ）を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。
【００３３】
【数１０】
【００３４】
式中、Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。
【００３５】
上記式（２）は、蛍光強度ｘｊが所定値（Ｕｊ）よりも大きいという制約条件（下限値条件）を表す。
【００３６】
ここで、ｓｉｊを要素とするＮ×Ｍ次行列Ｓと、ｘｊを要素とするＭ次行列ｘと、ｐｉを要素とするＮ次ベクトルｐを設定すると、式（２）を満足しながら式（１）で示される評価関数が最小値となるパラメータｘｊを求めることは、次の問題を解くこと同じである。
【００３７】
【数１１】
【００３８】
Ｍ×Ｍ次行列ＡとＭ×１次行列ｂは、蛍光強度ｘが所定値（Ｕ）以上であるように制約する不等式として以下の式（６），（７）のように設定される。
【００３９】
【数１２】
【００４０】
上記式（３）は２乗して展開すると、下記式（２３）のようになる。
【００４１】
【数１３】
【００４２】
式（２３）を最小化するため、最終項のｐＴｐは無視できる。従って、式（２３）を最小化することは、下記式（２４）を最小化することに等しい。
【００４３】
【数１４】
【００４４】
式（２４）及び上記式（４），（５）は、二次計画問題として知られている。「二次計画問題（Quadratic Programming Problem）」は、Ｄをｎ×ｎ次非負定値対称行列、ｃをｎ次ベクトル、Ａをｍ×ｎ次行列、ｂをｍ次ベクトルとするとき、線形制約条「Ａｘ≦ｂ、ｘ≧０」を満たすｎ次ベクトルのうち、２次形式で与えられる目的関数「f（ｘ）=ｘTAｘ／２＋ｃＴｘ」を最小にするものを求める問題である。二次計画問題は、次の式（２５）〜（２７）で表される。二次計画問題は、有限回の反復で厳密な最適解が得られる特殊な非線形計画問題として知られている。
【００４５】
【数１５】
【００４６】
上記式（２４），（４），（５）は、式（２５）中のＤ，ｃを以下のように置き換え、式（２６）のＡ，ｂを式（６），（７）のように置き換えて二次計画問題を解くことに等しい。
【数１６】
【００４７】
（２−３）制約条件（下限値条件）
（２−３−１）平均値（第一実施形態とその変形例）
下記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、蛍光色素を標識していない無染色の微小粒子を測定したとして、ｊ番目の蛍光色素にとって最小の検出値として適当な値を設定する。具体的には、下限値Ｕｊは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊを用いることができる。無染色平均値Ｖｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値が用いられる。あるいは、無染色平均値Ｖｊは、最大の検出値をとる検出器の検出平均値とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の検出平均値のさらに平均をとった値などとすることもできる。
【００４８】
【数１７】
【００４９】
【数１８】
【００５０】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用できる。この場合、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００５１】
【数１９】
【００５２】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。無染色平均値Ｖｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値が用いられる。あるいは、無染色平均値Ｖｊは、最大の検出値をとる検出器の検出平均値とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の検出平均値のさらに平均をとった値などとすることもできる。
【００５３】
（２−３−２）確率密度関数に従う乱数（第二実施形態とその変形例）
また、上記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された光検出器毎の検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いることもできる。無染色確率密度関数Ｆｊ（Ｕｊ）としては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の確率密度関数が用いられる。あるいは、無染色確率密度関数Ｆｊ（Ｕｊ）は、最大の検出値をとる検出器の確率密度関数とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の確率密度関数とから算出した関数などとすることもできる。
【００５４】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用する場合、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００５５】
【数２０】
【００５６】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。無染色確率密度関数Ｆｊ（Ｕｊ）としては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の確率密度関数ｆｉ（ｕｉ）が用いられる。あるいは、無染色確率密度関数Ｆｊ（Ｕｊ）は、最大の検出値をとる検出器の確率密度関数とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の確率密度関数とから算出した関数などとすることもできる。
【００５７】
（２−３−３）正規分布に従う乱数（第三実施形態とその変形例）
さらに、上記式（３）〜（５）の二次計画問題を解く際に制約条件とする下限値Ｕｊ（ｊ＝１〜Ｍ）には、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、前記下限値Ｕｊとして、下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いてもよい。この乱数は、正規分布に従った乱数（正規乱数）となる。無染色平均値Ｖｊ及び無染色標準偏差ρｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値及び標準偏差が用いられる。あるいは、無染色平均値Ｖｊ及び無染色標準偏差ρｊは、最大の検出値をとる検出器の平均値及び標準偏差とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の平均値及び標準偏差とから算出した値などとすることもできる。
【００５８】
【数２１】
【００５９】
変形例として、上記式（４）の線形制約条件がない二次計画法を適用する場合、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求めることもできる。
【００６０】
【数２２】
【００６１】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。無染色平均値Ｖｊ及び無染色標準偏差ρｊとしては、例えば、ｊ番目の蛍光色素の単染色スペクトルが最大の検出値をとる検出器における検出値の平均値ｉ及び標準偏差が用いられる。あるいは、無染色平均値Ｖｊ及び無染色標準偏差ρｊは、最大の検出値をとる検出器の平均値及び標準偏差とその前後の検出器（例えば（ｊ−１）番目と（ｊ＋１）番目の検出器）の平均値及び標準偏差とから算出した値などとすることもできる。
【００６２】
（２−３−４）０（第四実施形態）
上記式（４）の線形制約条件がない二次計画法を適用する場合、下限値を０とすれば、下記式（１０），（１１）の二次計画問題を解くことにより、パラメータｘｊ（ｊ＝１〜Ｍ）を求めることができる。
【００６３】
【数２３】
【００６４】
式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。
【００６５】
２．蛍光強度算出装置
本発明に係る蛍光強度算出装置は、従来のフローサイトメータ等と同様に流体系と光学系、分取系、データ処理系などから構成される。
【００６６】
流体系は、フローセルにおいて測定対象とする微小粒子を含むサンプル液をシース液の層流の中心に流し、フローセル内に微小粒子を一列に配列させる手段である。フローセルに替えて、マイクロチップ上に形成した流路内において、微小粒子を一列に配列させてもよい。
【００６７】
光学系は、蛍光色素により標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段である。光学系では、微小粒子から発生する前方散乱光や側方散乱光、レイリー散乱やミー散乱等の散乱光も検出される。光学系は、具体的には、レーザ光源と、微小粒子に対してレーザ光を集光・照射するための集光レンズやダイクロイックミラー、バンドパスフィルター等からなる照射系と、レーザ光の照射によって微小粒子から発生する蛍光や散乱光を検出する検出系と、によって構成される。検出系は、例えば、PMT（photo multiplier tube）や、CCDやCMOS素子等のエリア撮像素子等によって構成され、受光波長帯域の異なる光検出器が複数配設される。
【００６８】
微小粒子の分取を行う場合には、サンプル液を、微小粒子を含む液滴としてフローセル外の空間に吐出し、液滴の移動方向を制御して所望の特性を備えた微小粒子を分取する。分取系は、サンプル液を液滴化してフローセルから吐出させるピエゾ素子等の振動素子と、吐出される液滴に電荷を付与する荷電手段と、液滴の移動方向に沿って、移動する液滴を挟んで対向して配設された対電極などから構成される。
【００６９】
データ処理系は、光検出器から検出値を電気信号として入力され、電気信号に基づいて微小粒子の光学特性を解析する。データ処理系は、各光検出器から検出値を収集して得られる測定スペクトルを、上述した方法に従って単染色スペクトルの線形和により近似し、各蛍光色素から発生した真の蛍光強度を算出する。このため、データ処理系は、上述した本発明に係る蛍光強度算出方法の各ステップを実行するためのプログラムを記憶するハードディスク等の記録媒体と、プラグラムを実行するCPU及びメモリ等を有する。
【００７０】
データ処理系には、各光検出器からの信号の品質を向上させるため、各光検出器間のデータの平滑化を行うノイズフィルタを設けることが好ましい。平均化処理は、光検出器の数や受光波長帯域、使用する蛍光色素や、測定対象微小粒子の自家蛍光波長、さらに各機器のノイズ周波数成分等を考慮して行われる。データの平滑化を行うことにより、ノイズの影響を抑えることが可能となり、演算精度が向上し、より正確に蛍光強度を算出できる。
【実施例】
【００７１】
３２チャンネルのマルチカラー測定フローサイトメータでの仮想出力データを用いて、従来の補正行列を用いた蛍光補正方法と本発明に係る制限付き最小二乗法を用いた蛍光補正方法との間で演算結果の比較を行った。
【００７２】
蛍光色素には、ＦＩＴＣ（ＣＨ５），Ａｌｅｘａ５００（ＣＨ６），Ａｌｅｘａ５１４（ＣＨ１０），Ａｌｅｘａ５３２（ＣＨ１２），ＰＥ（ＣＨ１５），ＰＥ−ＴＲ（ＣＨ１９），ＰＩ（ＣＨ２１），Ａｌｅｘａ６００（ＣＨ２０），ＰＥ−Ｃｙ５（ＣＨ２４），ＰｅｒＣＰ（ＣＨ２５），ＰｅｒＣＰ−Ｃｙ５．５（ＣＨ２６），ＰＥ−Ｃｙ７（ＣＨ３１）を用いた。括弧内に示すチャネル番号は、補正行列を用いた従来方法における各蛍光色素に対するチャネルの割り当てを示す。
【００７３】
図２に、各蛍光色素の単染色スペクトルを示す。これらの各単染色スペクトルがランダムな強度で排他的に出現し、さらに各検出データには検出器や電気的なノイズが仮想的に付加されたシミュレーション・データをモンテカルロ法にて生成させた。発生させたシミュレーション・データを、図２の単染色スペクトルと、従来方法あるいは本発明に係る方法とを用いて蛍光補正を行って得た典型的な解析結果を以下に示す。
【００７４】
図３は、蛍光補正を行わずに、仮想出力データから作成した二次元相関図を示す。（Ａ）〜（Ｃ）の縦軸はそれぞれＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度を、横軸はＦＩＴＣの蛍光強度を示す。図上段はログスケールプロットした二次元相関図、下段はリニアスケールでプロットした二次元相関図である（以下、図４〜１１について同様）。
【００７５】
図４は、仮想出力データを、比較のために補正行列を用いて蛍光補正して作成した二次元相関図を示す。リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が負値として算出されているサブポピュレーションが存在している。このため、ログスケールの二次元相関図上にプロットされたポピュレーションが、リニアスケールの二次元相関図上にプロットされたポピュレーションに対して顕著に減少してしまっている。
【００７６】
図５は、仮想出力データを、下限値を平均値とする第一実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が正値として算出されている。このため、リニアスケールの二次元相関図上にプロットされたポピュレーションが、全てログスケールの二次元相関図上にもプロットできている。
【００７７】
図６は、仮想出力データを、線形制約条件なしで、下限値を平均値とする第一実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００７８】
図７は、仮想出力データを、下限値を確率密度関数に従う乱数とする第二実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００７９】
図８は、仮想出力データを、線形制約条件なしで、下限値を確率密度関数に従う乱数とする第二実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８０】
図９は、仮想出力データを、下限値を正規乱数とする第三実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８１】
図１０は、仮想出力データを、線形制約条件なしで、下限値を正規乱数とする第三実施形態の変形例に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８２】
図１１は、仮想出力データを、下限値を０とする第四実施形態に係る方法により蛍光補正して作成した二次元相関図を示す。
【００８３】
図６〜１１においても、リニアスケールでプロットした二次元相関図において、ＦＩＴＣ及びＡｌｅｘａ５００，Ａｌｅｘａ５１４，Ａｌｅｘａ５３２の蛍光強度が正値として算出されている。このため、リニアスケールの二次元相関図上にプロットされたポピュレーションが、全てログスケールの二次元相関図上にもプロットできている。
【産業上の利用可能性】
【００８４】
本発明に係る蛍光強度補正方法等によれば、複数の蛍光色素により標識された微小粒子を複数の光検出器によってマルチカラー測定する場合に、各蛍光色素からの蛍光強度を正確に算出してユーザに提示できる。従って、本発明に係る蛍光強度補正方法等は、細胞等の微小粒子の特性をより詳細に解析するため寄与し得る。
【特許請求の範囲】
【請求項１】
蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、
前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法。
【請求項２】
前記手順において、下記式（２）を満足しながら下記式（１）で示される評価関数が最小値となるパラメータｘｊ（ｊ＝１〜Ｍ）を求めることにより、各蛍光色素から発生した蛍光の強度を算出する請求項１記載の蛍光強度補正方法。
（式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。）
【請求項３】
前記手順において、下記式（３）〜（５）の二次計画問題を解くことにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ただし、ＡはＭ×Ｍ次行列、ｂはＭ×１次行列として設定し、Ｕ１〜ＵＭは前記下限値を表すものとする。）
【請求項４】
前記式（２）中、前記下限値Ｕｊとして、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊを用いる請求項３記載の蛍光強度補正方法。
【請求項５】
前記式（２）中、前記下限値Ｕｊとして、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いる請求項３記載の蛍光強度補正方法。
【請求項６】
前記式（２）中、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、前記下限値Ｕｊとして、下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いる請求項３記載の蛍光強度補正方法。
【請求項７】
前記手順において、下記式（１０），（１１）の二次計画問題を解くことにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。）
【請求項８】
前記手順において、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。）
【請求項９】
前記手順において、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【請求項１０】
前記手順において、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【請求項１１】
蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段と、
各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により前記測定スペクトルを近似する算出手段と、を備え、
前記算出手段は、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法により行う蛍光強度算出装置。
【請求項１】
蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して得られる測定スペクトルを、各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により近似する手順を含み、
前記手順において、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法を用いて行う蛍光強度補正方法。
【請求項２】
前記手順において、下記式（２）を満足しながら下記式（１）で示される評価関数が最小値となるパラメータｘｊ（ｊ＝１〜Ｍ）を求めることにより、各蛍光色素から発生した蛍光の強度を算出する請求項１記載の蛍光強度補正方法。
（式中、ｓｉｊは、ｊ番目の蛍光色素の単染色スペクトルにおけるｉ番目の光検出器の検出値を表す。ｐｉは、測定スペクトルにおけるｉ番目の光検出器の検出値を表す。σｉは、ｉ番目の光検出器の検出値に対する重みの逆数を表す。Ｕｊは、算出する各蛍光色素の蛍光強度の下限値を表す。）
【請求項３】
前記手順において、下記式（３）〜（５）の二次計画問題を解くことにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ただし、ＡはＭ×Ｍ次行列、ｂはＭ×１次行列として設定し、Ｕ１〜ＵＭは前記下限値を表すものとする。）
【請求項４】
前記式（２）中、前記下限値Ｕｊとして、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値から求めたｊ番目の蛍光色素の無染色平均値Ｖｊを用いる請求項３記載の蛍光強度補正方法。
【請求項５】
前記式（２）中、前記下限値Ｕｊとして、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の確率密度関数から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いる請求項３記載の蛍光強度補正方法。
【請求項６】
前記式（２）中、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得された前記光検出器毎の検出値の平均値及び分散からｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、前記下限値Ｕｊとして、下記式（９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数を用いる請求項３記載の蛍光強度補正方法。
【請求項７】
前記手順において、下記式（１０），（１１）の二次計画問題を解くことにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。）
【請求項８】
前記手順において、下記式（１２），（１３）の二次計画問題を解き、下記式（１４）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉから求めたｊ番目の蛍光色素の無染色平均値Ｖｊを要素とするＭ次ベクトルを表す。）
【請求項９】
前記手順において、下記式（１５），（１６）の二次計画問題を解き、下記式（１７）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の確率密度関数ｆｉ（ｕｉ）に従って発生させた乱数ｕｉを要素とするＮ次ベクトルを表す。Ｕは、確率密度関数ｆｉ（ｕｉ）から求めたｊ番目の蛍光色素の無染色確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【請求項１０】
前記手順において、下記式（２０），（２１）の二次計画問題を解き、下記式（２２）を実行することにより、前記パラメータｘｊ（ｊ＝１〜Ｍ）を求める請求項２記載の蛍光強度補正方法。
（式中、Ｓは、ｓｉｊを要素とするＮ×Ｍ次行列を表す。ｘは、ｘｊを要素とするＭ次行列を表す。ｐは、ｐｉを要素とするＮ次ベクトルを表す。ｕは、ｊ番目の蛍光色素を標識していない微小粒子に光を照射して取得されたｉ番目の光検出器の検出値の平均値ｖｉ及び分散σｉを求め、式（１８）の確率密度関数ｆｊ（ｕｊ）に従って発生させた乱数ｕｊを要素とするＮ次ベクトルを表す。Ｕは、平均値ｖｉ及び分散σｉからｊ番目の蛍光色素の無染色平均値Ｖｊ及び無染色標準偏差ρｊを求め、式（１９）の確率密度関数Ｆｊ（Ｕｊ）に従って発生させた乱数Ｕｊを要素とするＭ次ベクトルを表す。）
【請求項１１】
蛍光波長帯域の重複する複数の蛍光色素により多重標識された微小粒子に光を照射することによって励起された蛍光色素から発生する蛍光を、受光波長帯域の異なる光検出器で受光し、各光検出器から検出値を収集して測定スペクトルを得る測定手段と、
各蛍光色素を個別に標識した微小粒子で得られる単染色スペクトルの線形和により前記測定スペクトルを近似する算出手段と、を備え、
前記算出手段は、前記単染色スペクトルの線形和による前記測定スペクトルの近似を、制限付き最小二乗法により行う蛍光強度算出装置。
【図１２】
【図１３】
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１３】
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１２−１８１０８（Ｐ２０１２−１８１０８Ａ）
【公開日】平成２４年１月２６日（２０１２．１．２６）
【国際特許分類】
【出願番号】特願２０１０−１５６３８２（Ｐ２０１０−１５６３８２）
【出願日】平成２２年７月９日（２０１０．７．９）
【出願人】（０００００２１８５）ソニー株式会社 (34,172)
【Ｆターム（参考）】
【公開日】平成２４年１月２６日（２０１２．１．２６）
【国際特許分類】
【出願日】平成２２年７月９日（２０１０．７．９）
【出願人】（０００００２１８５）ソニー株式会社 (34,172)
【Ｆターム（参考）】
[ Back to top ]