(57)【要約】【課題】染色等の細胞に影響を与える処理を行うことなく細胞の活性度を測定する。【解決手段】複数の細胞Ａを含む細胞Ａ群の位相情報を測定する位相情報測定部９と、該位相情報測定部９により測定された位相情報から各細胞Ａの有する特徴量を抽出する特徴量抽出部１０と、該特徴量抽出部１０により抽出された特徴量に基づいて細胞Ａ群全体の活性度を算出する活性度算出部１２とを備える活性度測定装置１を提供する。【選択図】図１
2010276585.tif 000002
【特許請求の範囲】
【請求項１】
複数の細胞を含む細胞群の位相情報を測定する位相情報測定部と、
該位相情報測定部により測定された位相情報から各前記細胞の有する特徴量を抽出する特徴量抽出部と、
該特徴量抽出部により抽出された特徴量に基づいて前記細胞群全体の活性度を算出する活性度算出部とを備える活性度測定装置。
【請求項２】
前記特徴量抽出部が、特徴量として前記細胞の最大位相値を抽出する請求項１に記載の活性度測定装置。
【請求項３】
前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、全体頻度に対する所定の閾値以上の最大位相値を有する細胞の出現頻度の割合を活性度として算出する請求項２に記載の活性度測定装置。
【請求項４】
前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、細胞活性度毎に予め計算された少なくとも２つ以上の前記最大位相値の正規分布に基づいて、該正規分布にそれぞれ係数をかけて重畳させたものとヒストグラムとの差分が最小となる各正規分布の前記係数の割合を活性度として算出する請求項２に記載の活性度測定装置。
【請求項５】
前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、該ヒストグラムに現れる極小値以上の最大位相値を有する細胞の出現頻度の割合を活性度として算出する請求項２に記載の活性度測定装置。
【請求項６】
前記極小値が複数存在する場合に、各極小値となる最大位相値における出現頻度とその両側に隣接する最大位相値における出現頻度との差分の小さい方の値が、前記ヒストグラム内において最大となる最大位相値における出現頻度を、前記活性度の算出のための極小値とする請求項５に記載の活性度測定装置。
【請求項７】
複数の細胞を含む細胞群の位相情報を測定する位相情報測定ステップと、
該位相情報測定ステップにより測定された位相情報から各前記細胞の有する特徴量を抽出する特徴量抽出ステップと、
該特徴量抽出ステップにより抽出された特徴量に基づいて前記細胞群全体の活性度を算出する活性度算出ステップとを含む活性度測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、活性度測定装置および活性度測定方法に関するものである。
【背景技術】
【０００２】
従来、細胞のような測定対象の生存率（バイアビリティ）等の活性度を測定する方法が知られている。この方法は、細胞自体または特定の分子を蛍光物質によって染色しておき、励起光を照射して発生する蛍光の強度を測定するものである（例えば、特許文献１〜３参照。）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平７−１９８７１４号公報
【特許文献２】特開２００６−２３０３３３号公報
【特許文献３】特開２００６−３１４２１４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、従来の生存率測定方法は、生存率を測定する測定対象である細胞を染色する必要があるため、一度測定に用いた細胞は、他の用途、例えば、移植等に利用することができず、廃棄しなければならないという不都合がある。
【０００５】
本発明は上述した事情に鑑みてなされたものであって、染色等の細胞に影響を与える処理を行うことなく細胞の活性度を測定することができる活性度測定装置および活性度測定方法を提供することを目的としている。
【課題を解決するための手段】
【０００６】
上記目的を達成するために、本発明は以下の手段を提供する。
本発明は、複数の細胞を含む細胞群の位相情報を測定する位相情報測定部と、該位相情報測定部により測定された位相情報から各前記細胞の有する特徴量を抽出する特徴量抽出部と、該特徴量抽出部により抽出された特徴量に基づいて前記細胞群全体の活性度を算出する活性度算出部とを備える活性度測定装置を提供する。
【０００７】
本発明によれば、位相情報測定部によって測定された細胞群の位相情報から、特徴量抽出部によって各細胞の有する特徴量が抽出され、その特徴量に基づいて活性度算出部により細胞群全体の活性度が算出される。
細胞群の位相情報は、細胞群に対して焦点位置をずらした複数の画像に基づいて測定することができる。
【０００８】
位相情報は細胞内部の液体の屈折率によって変化する。そして、生細胞は外部の培地と内部の液体とを所定の濃度を保ちながら交換しているため、培地より高い屈折率が維持されているが、死細胞は細胞壁が破れることによって内部の液体が外部に放出されて培地が浸入するために、培地の屈折率に近い屈折率に変化する。したがって、位相情報に基づいて抽出された各細胞の特徴量、例えば、位相値φを任意の屈折率ｎで除算することによって算出された見かけの最大厚さ寸法等は生細胞か死細胞かによって変化しており、これによって細胞群の活性度を精度よく算出することができる。すなわち、本発明によれば、細胞を蛍光染色する等の物理的な処理を施すことなく、細胞の活性度を測定することができる。
【０００９】
上記発明においては、前記特徴量抽出部が、特徴量として前記細胞の最大位相値を抽出してもよい。
このようにすることで、特徴量抽出部によって抽出された最大位相値を用いて、細胞の活性度を簡易に判定することができる。
【００１０】
また、上記発明においては、前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、全体頻度に対する所定の閾値以上の最大位相値を有する細胞の出現頻度の割合を活性度として算出してもよい。
細胞の有する特徴量は、細胞の大きさや構造等の細胞の個性によってばらつきがあるので、細胞単独では測定誤差が大きいが、細胞群全体においてその分布を示すヒストグラムを生成することで、生細胞と死細胞との割合を求めることができ、活性度を精度よく測定することができる。
【００１１】
また、上記発明においては、前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、細胞活性度毎に予め計算された少なくとも２つ以上の前記最大位相値の正規分布に基づいて、該正規分布にそれぞれ係数をかけて重畳させたものとヒストグラムとの差分が最小となる各正規分布の前記係数の割合を活性度として算出することとしてもよい。
【００１２】
また、上記発明においては、前記活性度算出部が、前記特徴量抽出部により抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成し、該ヒストグラムに現れる極小値以上の最大位相値を有する細胞の出現頻度の割合を活性度として算出してもよい。
このようにすることで、生細胞と死細胞との特徴量が離れた分布を示している場合に、間に極小値が現れるので、これを用いて生細胞と死細胞との割合を求めることができ、活性度を精度よく測定することができる。
【００１３】
また、上記発明においては、前記極小値が複数存在する場合に、各極小値となる最大位相値における出現頻度とその両側に隣接する最大位相値における出現頻度との差分の小さい方の値が、前記ヒストグラム内において最大となる最大位相値における出現頻度を、前記活性度の算出のための極小値としてもよい。
【００１４】
また、本発明は、複数の細胞を含む細胞群の位相情報を測定する位相情報測定ステップと、該位相情報測定ステップにより測定された位相情報から各前記細胞の有する特徴量を抽出する特徴量抽出ステップと、該特徴量抽出ステップにより抽出された特徴量に基づいて前記細胞群全体の活性度を算出する活性度算出ステップとを含む活性度測定方法を提供する。
【発明の効果】
【００１５】
本発明によれば、染色等の細胞に影響を与える処理を行うことなく細胞の活性度を測定することができるという効果を奏する。
【図面の簡単な説明】
【００１６】
【図１】本発明の一実施形態に係る活性度測定装置を示す模式図である。
【図２】図１の活性度測定装置を用いた本実施形態に係る活性度測定方法を示すフローチャートである。
【図３】図１の活性度測定装置の位相情報測定部により生成された位相分布画像の一例を示す図である。
【図４】図１の活性度測定装置の特徴量抽出部により最大位相値が抽出された位相分布画像の一例を示す図である。
【図５】図１の活性度測定装置のヒストグラム生成部により生成されたヒストグラムであって、所定の閾値を基準として活性度を算出する場合を示す図である。
【図６】図５の変形例であって、発生頻度が最小値となる最大位相差を基準として活性度を算出する場合を示す図である。
【図７】図５の変形例であって、正規分布の合成によって活性度を算出する場合を示す図である。
【図８】（ａ）図５の方法、（ｂ）図７の方法でそれぞれ活性度を算出した結果を示す図である。
【発明を実施するための形態】
【００１７】
本発明の一実施形態に係る活性度測定装置１および活性度測定方法について、図面を参照して以下に説明する。
本実施形態に係る活性度測定装置は、図１に示されるように、細胞Ａ群を培地Ｂとともに収容したシャーレ２を搭載するステージ３と、該ステージ３の鉛直下方に配置され、所定の波長の光、例えば、近赤外光を発生する光源４と、該光源４から発せられた近赤外光を鉛直方向に配される光軸Ｃに沿って細胞群Ａに照射する照明光学系５と、細胞群Ａを上方に透過した透過光を集光する集光光学系６と、該集光光学系６によって集光された透過光を撮影するＣＣＤのような撮像素子７と、該撮像素子７により取得された画像に基づいて細胞Ａの活性度を判定する細胞活性度判定部８とを備えている。
【００１８】
ステージ３は、集光光学系６の光軸Ｃ方向に移動可能に設けられており、ステージ３を移動させた複数の位置において撮影することにより、１枚のフォーカス画像と２枚のデフォーカス画像とを取得することを可能にしている。なお、ステージ３を固定して集光光学系６を光軸Ｃ方向に移動可能に設けてもよい。
【００１９】
細胞活性度判定部８は、撮像素子７により取得された１枚のフォーカス画像と２枚のデフォーカス画像とを用いた公知の方法によって、位相情報の２次元的な分布を示す位相分布画像を取得する位相情報測定部９と、取得された位相分布画像に基づいて、該位相分布画像内に含まれる複数の細胞Ａの最大位相値を抽出する特徴量抽出部１０と、抽出された最大位相値の出現頻度の分布を示すヒストグラムを生成するヒストグラム生成部１１と、生成されたヒストグラムの中で、予め定められた閾値を超える頻度で出現する最大位相値の割合に基づいて細胞Ａ群の活性度を算出する活性度算出部１２とを備えている。
【００２０】
このように構成された本実施形態に係る活性度測定装置１を用いた活性度測定方法について以下に説明する。
まず、図２に示されるように、複数の細胞Ａを含む細胞Ａ群を収容したシャーレ２をステージ３に載置し、光源４から発生した近赤外光を照明光学系５を介してステージ３上の細胞Ａ群に照射する。細胞Ａ群を透過した透過光を集光光学系６によって集光し、撮像素子７によって撮影する。この場合には、ステージ３を光軸Ｃ方向に移動して１枚のフォーカス画像と２枚のデフォーカス画像とを取得する（ステップＳ１）。
【００２１】
次に、取得されたフォーカス画像およびデフォーカス画像に基づいて、位相情報測定部９が、位相情報の２次元的な分布を示す位相分布画像を生成する（ステップＳ２）。そして、特徴量抽出部１０が、位相分布画像に基づいて、各細胞Ａの有する最大位相値を抽出する（ステップＳ３）。特徴量抽出部１０による最大位相値の抽出は、画像処理により各細胞Ａの輪郭形状を抽出するとともに、抽出された各細胞Ａの輪郭形状内において位相値が最大となる箇所を探索することにより行われる。
【００２２】
この後に、ヒストグラム生成部１１が、抽出された各細胞Ａの最大位相値の大きさと出現頻度との関係を示すヒストグラムを生成する（ステップＳ４）。そして、生成されたヒストグラムに基づいて活性度算出部１２が細胞Ａ群全体の活性度を算出する（ステップＳ５）。
【００２３】
具体的には、図３に示されるように、ステップＳ２において位相分布画像が生成され、図４に示されるように、ステップＳ３において該位相分布画像を処理して各細胞Ａの有する最大位相値が抽出される。図中、符号＋は、抽出された各細胞Ａの最大位相値の位置を示している。
そして、図５に示されるように、ステップＳ４において、抽出された最大位相値を横軸とし、その出現頻度を縦軸としたヒストグラムが生成される。
【００２４】
生細胞は細胞膜内に細胞内液を保有し、外部の培地Ｂとは異なる屈折率を保持しているので、細胞を通過する近赤外光と細胞以外の場所を通過する近赤外光とでは、屈折率の相違に基づく位相差が発生し、その位相差の分布が位相分布画像として生成される。
死細胞も外部の培地とは異なる屈折率を有しているものの、細胞膜が破壊されることにより、細胞内液が放出され外部の培地が入り込んで培地に近い屈折率になっている。したがって、生細胞とはその位相差が異なっている。
【００２５】
このため、生細胞と死細胞とが同じ厚さ寸法を有していても、屈折率の相違により位相分布が異なり、最大厚さ寸法に対応する最大位相値の値も相違してくる。各細胞Ａは、その個体差によって最大厚さ寸法が異なるが、その分布は所定の平均値の周囲に正規分布をなして分散している。したがって、多数の細胞Ａについて最大位相値の出現頻度を示すヒストグラムを作成することにより、生細胞と死細胞とが混在した細胞Ａ群の最大位相値の分布を得ることができる。
【００２６】
細胞内液の屈折率は培地Ｂの屈折率よりも高いので、生細胞の屈折率は死細胞の屈折率より高く、最大位相値の大きい側に生細胞が分布している。したがって、所定の閾値より大きな最大位相値の出現頻度の割合が多い程、生細胞が多いと考えることができる。
【００２７】
そこで、生成されたヒストグラムにおいて、図５に示すように、所定の閾値を超える最大位相値の出現頻度の全データ数に対する割合を求めることにより、細胞Ａの活性度を算出することができる。逆に、所定の閾値以下の最大位相値の出現頻度の全データ数に対する割合を求めることにより、細胞Ａの死亡率を算出することができる。ここで、所定の閾値は、細胞Ａの種類毎に予めデータを採取することにより求めておくことにすればよい。
【００２８】
本実施形態に係る活性度測定装置１および活性度測定方法によれば、細胞Ａ群に近赤外光を照射して生成された位相分布画像を用いて細胞Ａの活性度を算出するので、従来、染色することによって測定していた細胞Ａの活性度を無線色のまま測定することができる。したがって、活性度の測定に使用した細胞Ａを廃棄することなくそのまま用いることができ、貴重な細胞Ａを無駄にせずに済むという利点がある。
【００２９】
なお、本実施形態においては、予めデータを採取して求めておいた閾値を用いて、当該閾値を超える最大位相値を有する細胞Ａの割合によって細胞Ａ群の活性度を評価することとしたが、これに代えて、図６に示されるように、ヒストグラムに現れる極小値の最大位相値を基準として活性度を評価することにしてもよい。生細胞と死細胞との間の屈折率の差が大きい場合には、ヒストグラムにおける生細胞と死細胞との分布の間に極小値が現れるので、それを基準として活性度を評価することができる。
【００３０】
ヒストグラムの中に複数の極小値が現れる場合には、注目極小値の前のデータとの差分と後のデータとの差分を計算し、差分が小さい方を記録する。同様の処理をヒストグラム内の全データで実施する。得られた差分データ内で値が最大となる極小値を閾値とする。
【００３１】
また、生細胞と死細胞のそれぞれの正規分布を予め用意しておき、それらの正規分布にそれぞれ係数をかけたものを合成した分布が、取得されたヒストグラムに近似するように係数を調節し、得られた係数の割合によって活性度を評価することにしてもよい。
例えば、図７に示す例では、（ａ）に示されるヒストグラムに対し、（ｂ）に示す生細胞の正規分布ＸをＰ倍し、死細胞の正規分布ＹをＱ倍することで、（ｃ）に示されるように取得されたヒストグラムの分布に近似する分布が得られた場合に、細胞Ａ群の活性度はＰ／（Ｐ＋Ｑ）によって算出することができ、細胞群の死亡率は、Ｑ／（Ｐ＋Ｑ）によって算出することができる。
【００３２】
ここで、本実施形態に係る活性度測定装置１を用いた細胞Ａ群の活性度測定の実施例について説明する。
まず、シャーレ２内に脂肪由来幹細胞Ａを培地Ｂとともに収容したものを６組用意し、３組の脂肪由来幹細胞Ａには紫外線刺激を４時間行い、他の３組の脂肪由来幹細胞Ａは常温で放置した。
【００３３】
各シャーレ２内の脂肪由来幹細胞Ａに対して、本実施形態に係る活性度測定装置１を用いて、脂肪由来幹細胞Ａの位相分布画像を、視野を変えながら１０枚取得した。次いで、各位相分布画像から各細胞Ａの最大位相値を抽出した。
最大位相値を横軸、出現頻度を縦軸としてヒストグラムを生成し、予め定めた閾値に基づく方法による細胞死亡率と、正規分布Ｘ，Ｙの係数Ｐ，Ｑを調整する方法による細胞死亡率とを求めた結果、図８（ａ），（ｂ）に示す通りの結果が得られた。
【００３４】
各組の脂肪由来細胞Ａに対しては、比較のために従来と同様のNucleoCounterによる細胞死亡率の測定も合わせて行った。
その結果、本実施形態に係る活性度測定装置１を用いたいずれの方法によっても、従来の測定方法と同様の結果を得ることができた。
【符号の説明】
【００３５】
Ａ  細胞
１  活性度測定装置
９  位相情報測定部
１０  特徴量抽出部
１２  活性度算出部
Ｓ２  位相分布画像生成ステップ（位相情報測定ステップ）
Ｓ３  最大位相値抽出ステップ（特徴量抽出ステップ）
Ｓ５  活性度算出ステップ
【図１】
2010276585.tif 000003
【図２】
2010276585.tif 000004
【図５】
2010276585.tif 000005
【図６】
2010276585.tif 000006
【図７】
2010276585.tif 000007
【図８】
2010276585.tif 000008
【図３】
2010276585.tif 000009
【図４】
2010276585.tif 000010