(57)【要約】複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することができる判定方法を提供する。複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像と前記スフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定する。
2019176427.tif 000002
【特許請求の範囲】
【請求項１】
複数の細胞の凝集体であるスフェアを撮像したホログラムから前記スフェアの位相差画像を生成し、
前記位相差画像と前記スフェアの形状に応じた形状指標値とに基づいて、前記スフェアの状態を判定する
判定方法。
【請求項２】
前記位相差画像と前記形状指標値とに基づいて、前記スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及び前記スフェアの外形形状の少なくとも１つに関する判定を行う
請求項１に記載の判定方法。
【請求項３】
前記位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、
前記位相差量総和と前記形状指標値との相関性を用いて前記スフェアの状態を判定する
請求項１または請求項２に記載の判定方法。
【請求項４】
前記位相差量総和と前記形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、
前記判定対象のスフェアについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の程度に応じて前記判定対象のスフェアの状態を判定する
請求項３に記載の判定方法。
【請求項５】
前記判定対象のスフェアについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の幅が閾値を超える場合、前記判定対象のスフェアに含まれる細胞の生存率、密度、均質性及び前記スフェアの外形形状の少なくとも１つについて異常ありと判定する
請求項４に記載の判定方法。
【請求項６】
前記位相差量総和と前記形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、
前記判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、前記基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて前記判定対象の培養ロットの良否判定を行う
請求項３に記載の判定方法。
【請求項７】
前記形状指標値は、前記スフェアの体積、断面積、粒径及び周の長さのいずれかである
請求項１から請求項６のいずれか１項に記載の判定方法。
【請求項８】
前記位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、前記スフェアの体積で除算して位相差量密度を導出し、
前記位相差量密度に基づいて、前記スフェアの状態を判定する
請求項１または請求項２に記載の判定方法。
【請求項９】
前記スフェアの状態の判定に用いる位相差画像は、前記ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像である
請求項１から請求項８のいずれか１項に記載の判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
開示の技術は、複数の細胞の凝集体であるスフェアの状態を判定する判定方法に関する。
【背景技術】
【０００２】
細胞の状態を評価または判定する技術として、例えば、以下の技術が知られている。特許文献１には、細胞群に対して照明光を照射し、照明光の照射によって細胞群を透過した透過光を検出して細胞群の画像を検出し、細胞群における散乱光を検出するための検出光を細胞群に対して照射し、検出した散乱光の強度の空間的な広がりまたは時間的な揺らぎに基づいて細胞群を評価する細胞評価方法が記載されている。【０００３】
特許文献２には、神経細胞分化過程の細胞を撮像した第１の撮像画像を入力する画像入力部と、細胞の厚さ方向の位置において焦点の合った状態の全焦点画像を、少なくとも第１の撮像画像に基づく第１の原画像として生成する全焦点画像生成部と、第１の原画像において一定以下の輝度分布を有する領域と、第１の原画像において濃度差が一定以下の領域とで共通する領域を、分化過程に出現するロゼッタとして抽出するロゼッタ抽出部と、抽出されたロゼッタの状態を判定するロゼッタ対応判定部と、を備える細胞評価装置が記載されている。【０００４】
特許文献３には、細胞厚さに基づいて、細胞の分化度を判定する判定ステップを備えることを特徴とする細胞判定方法が記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−００５４３７号公報
【特許文献２】特開２０１７−０００１６３号公報
【特許文献３】特開２０１５−１４６７４７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
細胞の大量生産が可能な培養手法として、細胞の凝集体であるスフェアを培地中に浮遊させた状態で培養する三次元培養法が知られている。三次元培養による細胞の製造工程においては、工程管理の容易性の観点から、細胞の品質をスフェアの状態のまま非破壊且つ簡便に評価する技術が求められる。しかしながら、現時点においては、三次元空間中にランダムに存在する様々なサイズのスフェアを評価する手法が確立されておらず、特にスフェアの内部における細胞の密度及び生存状況について直接観察することは困難である。このため、従来の二次元培養の手法を適用した評価がなされているが、培養する細胞の数の増加に伴って、評価工数が増えるため、多くの人手及び多くの時間が必要となる。従来の二次元培養の手法を適用した評価では、スフェアを単一細胞に分解したり、蛍光色素を添加したりするといった細胞の破壊を伴う処理が必要とされる。【０００７】
開示の技術は、複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することを目的とする。
【課題を解決するための手段】
【０００８】
開示の技術に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定することを含む。【０００９】
開示の技術に係る判定方法によれば、スフェアの状態を非破壊且つ簡便に判定することが可能となる。【００１０】
開示の技術に係る判定方法において、位相差画像と形状指標値とに基づいて、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状の少なくとも１つに関する判定を行ってもよい。【００１１】
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、位相差量総和と形状指標値との相関性を用いてスフェアの状態を判定してもよい。【００１２】
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の程度に応じて判定対象のスフェアの状態を判定してもよい。例えば、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。【００１３】
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて判定対象の培養ロットの良否判定を行ってもよい。【００１４】
開示の技術に係る判定方法において、形状指標値は、スフェアの体積、断面積、粒径及び周の長さのいずれかであってもよい。【００１５】
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、スフェアの体積で除算して位相差量密度を導出し、位相差量密度に基づいて、スフェアの状態を判定してもよい。【００１６】
開示の技術に係る判定方法において、スフェアの状態の判定に用いる位相差画像は、ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像であることが好ましい。
【発明の効果】
【００１７】
開示の技術によれば、複数の細胞の凝集体であるスフェアの状態を非破壊且つ簡便に判定することが可能となる。
【図面の簡単な説明】
【００１８】
【図１】開示の技術の実施形態に係る判定方法の実施に用いる撮像システムの構成の一例を示す図である。
【図２Ａ】開示の技術の実施形態に係る判定方法の実施に用いるホログラムの一例を示す図である。
【図２Ｂ】スフェアのフーリエ変換画像の一例を示す図である。
【図２Ｃ】スフェアのアンラッピング前の位相差画像の一例を示す図である。
【図２Ｄ】スフェアのアンラッピング後の位相差画像の一例を示す図である。
【図３】開示の技術の実施形態に係る位相差画像の概念を示す図である。
【図４】開示の技術の実施形態に係る位相差画像の焦点合わせに関する説明図である。
【図５】開示の技術の実施形態に係るオートフォーカス処理を行うコンピュータのハードウェア構成の一例である。
【図６】開示の技術の実施形態に係るオートフォーカス処理の流れの一例を示すフローチャートである。
【図７】開示の技術の実施形態に係る、スフェアの位相差画像における焦点位置と位相差量のばらつきとの関係の一例を示すグラフである。
【図８】開示の技術の実施形態に係る、スフェアの体積と位相差量総和との相関特性の一例を示すグラフである。
【図９Ａ】開示の技術の実施形態に係る、良ロット及び不良ロットの各々について導出したスフェアの粒径のヒストグラムである。
【図９Ｂ】開示の技術の実施形態に係る、良ロット及び不良ロットの各々について取得した、スフェアの体積と位相差量総和との相関特性の一例を示すグラフである。
【図９Ｃ】開示の技術の実施形態に係る、良ロット及び不良ロットの各々について取得した、位相差量密度のヒストグラムである。
【図１０】開示の技術の実施形態に係る、細胞の生存率と位相差量密度との相関特性の一例を示すグラフである。
【図１１Ａ】Ｈ２Ｏ２を添加した後、０分経過した時点におけるコロニーの蛍光顕微鏡画像である。
【図１１Ｂ】Ｈ２Ｏ２を添加した後、１５分経過した時点におけるコロニーの蛍光顕微鏡画像である。
【図１１Ｃ】Ｈ２Ｏ２を添加した後、３７分経過した時点におけるコロニーの蛍光顕微鏡画像である。
【図１２Ａ】Ｈ２Ｏ２を添加した後、０分経過した時点におけるコロニーの位相差画像である。
【図１２Ｂ】Ｈ２Ｏ２を添加した後、１５分経過した時点におけるコロニーの位相差画像である。
【図１２Ｃ】Ｈ２Ｏ２を添加した後、３７分経過した時点におけるコロニーの位相差画像である。
【図１３】開示の技術の実施形態に係る、ネクローシスに誘導したコロニー、アポトーシスに誘導したコロニー及び健全なコロニーのそれぞれの位相差量総和の相対値を示すグラフである。
【図１４】開示の技術の実施形態に係る、未分化状態を維持したスフェア及び未分化状態を逸脱したスフェアのそれぞれの位相差量総和の平均値を示すグラフである。
【発明を実施するための形態】
【００１９】
以下、本発明の実施形態について図面を参照しつつ説明する。尚、各図面において、実質的に同一又は等価な構成要素又は部分には同一の参照符号を付している。
【００２０】
開示の技術の実施形態に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定する、というものである。
【００２１】
図１は、開示の技術の実施形態に係る判定方法の実施に用いる撮像システム１の構成の一例を示す図である。撮像システム１は、公知のデジタルホログラフィ技術を用いてスフェアのホログラムを取得するためのホログラム光学系１０と、スフェアの蛍光顕微鏡観察を行うための蛍光顕微鏡光学系３０と、を含んで構成されている。
【００２２】
デジタルホログラフィ技術は、物体を透過または反射した物体光と、物体光に対してコヒーレントである参照光との干渉によって生じる像をイメージセンサーを用いて撮像し、撮像によって得られた画像について、光伝搬に基づく数値計算を実施することにより、物体からの光波の波面を復元する技術である。デジタルホログラフィ技術によれば、物体の位相分布を定量化し、また、焦点位置を機械的に移動させることなく、物体の三次元情報を取得することができる。
【００２３】
ホログラム光学系１０は、レーザ光源１１、ビームスプリッタ１２及び１８、コリメートレンズ１３及び２１、対物レンズ１５、結像レンズ１７並びにＣＭＯＳ（Complementary Metal Oxide Semiconductor）カメラ１９を含んで構成されている。サンプルステージにセットされる試料１４としてのスフェアは、コリメートレンズ１３と対物レンズ１５との間に配置される。
【００２４】
レーザ光源１１には、例えば波長６３２．８ｎｍのＨｅＮｅレーザを用いることが可能である。レーザ光源１１から出射されたレーザ光は、ビームスプリッタ１２により、２つのレーザ光に分割される。２つのレーザ光の一方は物体光とされ、他方は参照光とされる。物体光は、コリメートレンズ１３によって平行光とされた後、サンプルステージにセットされた試料１４であるスフェアに照射される。スフェアを透過した物体光による像は、対物レンズ１５によって拡大される。対物レンズ１５を透過した物体光は、結像レンズ１７によって再び平行光とされた後、ビームスプリッタ１８を介してＣＭＯＳカメラ１９の撮像面に結像される。一方、参照光は、光ファイバ２０によってコリメートレンズ２１の手前まで導かれる。光ファイバ２０から出射した参照光は、コリメートレンズ２１によって平行光とされ、ビームスプリッタ１８を介してＣＭＯＳカメラ１９の撮像面に入射する。物体光と参照光との干渉によって生じるホログラムが、ＣＭＯＳカメラ１９によって記録される。なお、ＣＭＯＳカメラ１９の撮像面に入射する物体光と参照光の光軸方向が互いに異なったオフアキシャル光学系が構成されていてもよい。
【００２５】
一方、蛍光顕微鏡光学系３０は、励起光源３１、励起用フィルタ３２、ＯＮ／ＯＦＦ切り替えミラー３３、ダイクロイックミラー３４、対物レンズ１５、分光フィルタ３５及びＣＭＯＳカメラ３６を含んで構成されている。
【００２６】
励起光源３１には、例えば水銀ランプを用いることが可能である。励起用フィルタ３２は、励起光源３１から出射された励起光のうち、例えば４５０ｎｍ〜４９０ｎｍの波長域の光を透過させる。励起用フィルタ３２を透過した励起光は、ＯＮ／ＯＦＦ切り替えミラー３３、ダイクロイックミラー３４及び対物レンズ１５を介してサンプルステージにセットされた試料１４であるスフェアに照射される。スフェアに添加された蛍光色素に励起光が照射されることにより、蛍光色素から発せられる蛍光は、対物レンズ１５、ダイクロイックミラー３４、ＯＮ／ＯＦＦ切り替えミラー３３及び分光フィルタ３５を介してＣＭＯＳカメラ３６の撮像面に入射する。分光フィルタ３５は、蛍光色素から発せられた蛍光のうち、例えば５１５ｎｍ以上の波長域の光を透過させる。
【００２７】
本実施形態に係る撮像システム１によれば、ホログラム光学系１０によるホログラム観察と、蛍光顕微鏡光学系３０による蛍光顕微観察とを並行して行うことが可能である。なお、上記した撮像システム１の構成は、一例にすぎず、上記の構成に限定されない。開示の技術に係る判定方法の実施には、デジタルホログラム技術を用いてホログラムを取得することができる、あらゆる撮像システムを利用することが可能である。
【００２８】
以下に、撮像システム１を用いて取得したスフェアのホログラムから、スフェアの位相差画像を取得する方法の一例について説明する。
【００２９】
はじめに、ＣＭＯＳカメラ１９によって取得した図２Ａに例示するホログラムを、二次元フーリエ変換することにより、物体光のみの複素振幅成分を抜き出す。図２Ｂは、この処理によって得られるスフェアのフーリエ変換画像の一例である。
【００３０】
次に、例えば角スペクトル法を適用して任意の空間位置のスフェアの位相を示す画像を復元する。図２Ｃは、この処理によって得られるスフェアのアンラッピング前の位相差画像の一例である。この時点におけるスフェアの位相は、０〜２πの値に畳みこまれている。そこで、例えば、Unweighted Least Squares（重みなし最小２乗法）またはFlynn's Algorithm（フリンのアルゴリズム）などの位相接続（アンラッピング）手法を適用して２π以上の部分も接合していくことにより、図２Ｄに例示するような最終的なスフェアの位相差画像を得ることができる。なお、アンラッピングの手法は数多く提案されており、位相不整合を生じない適切な手法を適宜選択すればよい。
【００３１】
図３は、位相差画像ＩＰの概念を示す図である。図３の下段には、位相差画像ＩＰの各画素ｋにおける位相差量を３次元表示する。図３の上段には、位相差画像ＩＰの各画素ｋにおける位相差量を平面上にグレースケールで示す。
【００３２】
ここで、位相差画像ＩＰにおける位相差量θは、位相差画像ＩＰの同一焦点面内に存在するバックグラウンド（スフェアの存在しない領域）の位相をθＢとし、スフェアの存在する領域の位相をθＳとしたとき、下記の（１）式によって表わされる。また、本明細書中における「位相」という用語は、光を電磁波とみなした場合の電場振幅の位相であり、より一般的な意味で使用される。
【数１】
2019176427.tif 000003
【００３３】
また、位相差画像ＩＰの各画素ｋにおける位相差量θｋは、下記（２）式によって表わすことができる。但し、ｎｋは位相差画像ＩＰの各画素ｋに対応する部位におけるスフェアの屈折率であり、ｄｋは位相差画像ＩＰの各画素ｋに対応する部位におけるスフェアの厚さであり、λはホログラム光学系１０における物体光の波長である。
【数２】
2019176427.tif 000004
【００３４】
スフェアの位相差画像は、スフェアを透過した物体光の光路長分布を示した画像である。スフェア内における光路長は、スフェアの屈折率とスフェアの厚さの積に相当することから、スフェアの位相差画像は、（２）式にも示されているように、スフェアの屈折率及び厚さ（形状）の情報を含んでいる。
【００３５】
スフェアに対して焦点が合っていない位相差画像からは、回折による広がりの影響によりスフェアの実態に合致した正確な情報が得られない。従って、ＣＭＯＳカメラ１９によって取得したホログラムから位相差画像を取得する際に、スフェアに焦点を合わせることが好ましい。ここで、「スフェアに焦点を合わせる」とは、球形状のスフェアの中央付近においてスライスした位相差画像を得ることを意味する。スフェアに焦点が合った位相差画像を用いてスフェアの状態を判定することにより、より正確な判定結果を得ることができる。なお、ユーザが、取得したスフェアの位相差画像に基づいて、スフェアの状態を判定してもよい。または、判定条件を学習したコンピュータが、取得したスフェアの位相差画像に基づいて、スフェアの状態を判定してもよい。
【００３６】
位相差画像の焦点合わせは、人手によらず、自動化することが好ましい。焦点合わせを自動化することにより、作業者による任意性を排除し、更に、処理時間の短縮を図ることができる。本発明者らは、以下に説明する自動化可能な焦点合わせの手法を見出した。
【００３７】
図４の左のグラフは、スフェアの位相差画像における平面方向の位置と位相差量との関係の一例を示すグラフであり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、位相差画像における特定の位置に急峻なピークが現れる。一方、スフェアに焦点が合っていない場合、焦点が合っている場合と比較してピークが低く且つなだらかになる。
【００３８】
図４の右のグラフは、スフェアの位相差画像における位相差量のヒストグラムの一例であり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、カーブの幅ｗ（位相差量のばらつき）は、相対的に大きくなり、スフェアに焦点が合っていない場合、カーブの幅ｗ（位相差量のばらつき）は、相対的に小さくなる。
【００３９】
従って、互いに異なる焦点位置（スライス位置）毎にスフェアの位相差画像を取得し、取得した位相差画像の各々について、位相差量のヒストグラムにおけるカーブの幅ｗ（位相差量のばらつき）を求め、求めた幅ｗのうち、最大の幅ｗを有する位相差画像を、スフェアに焦点が合った位相差画像として抽出することにより焦点合わせを実現できる。
【００４０】
上記の焦点合わせは、コンピュータを用いて自動化することが可能である。図５は、上記の焦点合わせを自動で行うオートフォーカス処理を行うコンピュータ５００のハードウェア構成の一例である。
【００４１】
コンピュータ５００は、ＣＰＵ（Central Processing Unit）５０１、一時記憶領域としての主記憶装置５０２、不揮発性の補助記憶装置５０３、ＣＭＯＳカメラ１９との間で通信を行うための通信Ｉ／Ｆ（InterFace）５０４、及び液晶ディスプレイ等の表示部５０５を含んで構成されている。ＣＰＵ５０１、主記憶装置５０２、補助記憶装置５０３、通信Ｉ／Ｆ５０４、及び表示部５０５は、それぞれ、バス５０７に接続されている。補助記憶装置５０３には、上記のオートフォーカス処理の手順を記述したオートフォーカスプログラム５０６が格納されている。コンピュータ５００は、ＣＰＵ５０１がオートフォーカスプログラム５０６を実行することにより、オートフォーカス処理を行う。
【００４２】
図６は、コンピュータ５００において実施されるオートフォーカス処理の流れの一例を示すフローチャートである。
【００４３】
ステップＳ１において、ＣＰＵ５０１は、ＣＭＯＳカメラ１９からスフェアのホログラムを取得する。
【００４４】
ステップＳ２において、ＣＰＵ５０１は、取得したホログラムから、焦点位置（スライス位置）が互いに異なる複数の位相差画像を生成する。
【００４５】
ステップＳ３において、ＣＰＵ５０１は、焦点位置（スライス位置）毎の位相差画像の各々について、位相差量のばらつきを導出する。ＣＰＵ５０１は、例えば、位相差画像における位相差量の最大値と最小値との差分を、当該位相差画像における位相差量のばらつきとして導出してもよい。
【００４６】
ステップＳ４において、ＣＰＵ５０１は、焦点位置（スライス位置）が互いに異なる複数の位相差画像のうち、ステップＳ３にて導出した位相差量のばらつきが最大となる位相差画像を、焦点が合った位相差画像として抽出する。
【００４７】
図７は、スフェアの位相差画像における焦点位置（スライス位置）と位相差量のばらつきとの関係の一例を示したグラフである。図７には、焦点位置が、−４００μｍ、−２００μｍ、０μｍ、＋２００μｍ、及び＋４００μｍに対応するスフェアの位相差画像、並びにグラフが例示されている。なお、図７では、位相差量のばらつきが最大となる焦点位置を０μｍとしている。上記のオートフォーカス処理によれば、位相差量のばらつきが最大となる焦点位置０μｍに対応する位相差画像が、焦点が合った位相差画像として抽出される。位相差量のばらつきが最大となる焦点位置０μｍに対応する位相差画像において、スフェアの輪郭が最も鮮明となる。
【００４８】
開示の技術の実施形態に係る判定方法は、スフェアの位相差画像と、スフェアの形状指標値とに基づいてスフェアの状態を判定することを含む。すなわち、スフェアの位相差画像を、スフェアの形状指標値と関連付けて分析することにより、スフェアの状態が判定される。判定されるスフェアの状態には、例えば、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状などが含まれる。
【００４９】
スフェアの形状指標値として、例えば、スフェアの体積、断面積、粒径及び周の長さなどを用いることができる。スフェアの粒径及び周の長さは、スフェアに焦点を合わせた位相差画像（すなわち、スフェアの中央付近においてスライスした位相差画像）におけるスフェアの像から直接的に求めることが可能である。スフェアの断面積は、例えば、スフェアの粒径を直径とする円の断面積として導出することが可能である。スフェアの体積は、例えばスフェアの粒径を直径とする球の体積として導出することが可能である。
【００５０】
開示の技術の実施形態に係る判定方法は、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和θＡを導出し、位相差量総和θＡと形状指標値との相関性を用いてスフェアの状態を判定することを含み得る。
【００５１】
位相差量総和θＡは、下記の（３）式によって表わされる。但し、ｓは位相差画像の各画素ｋの面積であり、ｖｋは位相差画像の各画素ｋに対応する部位におけるスフェアの体積である。（３）式に示されるように、位相差量総和θＡは、スフェアの位相差画像の画素毎の位相差量θｋを、全ての画素ｋについて積算した値に相当する。なお、（３）式におけるｄｋは画素ｋに投影されたスフェア部分の厚みを示し、ｎｋはバックグラウンドである培養液とスフェア内部との屈折率差を表す。（３）式では、ｖｋ＝ｄｋ・sを用いた。ここで、（３）式により、位相差量総和θＡの単位は面積のスケールとなって、例えば[μｍ２]であるが、イメージセンサー間で比較を行わない場合には単に１ピクセルあたりの画素毎の位相差量θｋの和として位相差量総和θＡの単位を[pixel]とし、つまりs ＝ 1 [pixel]として扱ってもよい。
【数３】
2019176427.tif 000005
【００５２】
開示の技術の実施形態に係る判定方法は、位相差量総和θＡとスフェアの形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象のスフェアについての相関性の、基準相関トレンドラインからの乖離の程度に応じて当該判定対象のスフェアの状態を判定することを含み得る。この場合において、判定対象のスフェアについての相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、当該判定対象のスフェアに含まれる細胞の生存率、密度、均質性及び前記スフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
【００５３】
例えば、一方の軸に位相差量総和θＡをとり、他方の軸に形状指標値をとったグラフ上にプロットされた、判定対象のスフェアについてのプロットの、基準相関トレンドラインからの乖離の程度に応じて、当該判定対象のスフェアの状態を判定してもよい。この場合において、上記プロットの、基準相関トレンドラインからの乖離の幅が閾値を超える場合、当該判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
【００５４】
開示の技術の実施形態に係る判定方法は、基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象の培養ロットについての相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて当該判定対象の培養ロットの良否判定を行うことを含み得る。
【００５５】
例えば、一方の軸に位相差量総和θＡをとり、他方の軸に形状指標値をとったグラフ上にプロットされた、判定対象の培養ロットに属する複数のスフェアの各々についてのプロットの、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて当該培養ロットの良否判定を行ってもよい。
【００５６】
例えば、健全なスフェア（標準サンプル）について予め取得した、位相差量総和θＡとスフェアの形状指標値との相関ラインを、基準相関トレンドラインとして用いることが可能である。また、判定対象のスフェアについての各プロットから、例えば最小二乗法を用いて導出される回帰線を、基準相関トレンドラインとして用いてもよい。
【００５７】
開示の技術の実施形態に係る判定方法は、位相差量総和θＡを当該スフェアの体積で除算して位相差量密度ＤＰを導出し、導出した位相差量密度ＤＰに基づいて、当該スフェアの状態を判定することを含み得る。
【００５８】
位相差量密度ＤＰは、下記の（４）式によって表わされる。但し、Ｖはスフェアの体積である。（４）式に示されるように、位相差量密度ＤＰは、位相差量総和θＡを当該スフェアの体積Ｖで除算した値に相当する。健全な細胞は、その恒常性から内部の屈折率は、媒質の屈折率とは異なる一定の値を維持すると考えられる。一方、死細胞は、恒常性を喪失し、内部の屈折率が媒質の屈折率と略同じになると考えられる。従って、位相差量密度ＤＰを、細胞の状態を示す指標として用いることが可能であると考えられる。例えば、判定対象のスフェアについて取得した位相差量密度ＤＰが、閾値以上である場合に当該スフェアの状態が良好であると判定することができ、位相差量密度ＤＰが、閾値未満である場合に当該スフェアの状態が異常であると判定することができる。なお、２π／λは、定数として扱うことができるので、位相差量密度ＤＰの導出に際し、２π／λの乗算を省略してもよい。ここで、スフェアの体積平均屈折率差ＮａｖｅをＮａｖｅ=Σｎｋ ・ （ｖｋ／Ｖ）とすると、（４）式は、ＤＰ = ( ２π／ λ) ×Ｎａｖｅとなることから、位相差密度は体積平均したスフェアの屈折率差を波長の長さで規格化した値である。本明細書において、Ｖはスフェアの位相像の断面画像から球相当径を算出して求めた。より正確に楕円球とすることも可能である。
【数４】
2019176427.tif 000006
【００５９】
開示の技術に係る判定方法の実施例を以下に示す。
【００６０】
［実施例１］
三次元培養法により培養したｉＰＳ細胞（induced pluripotent stem cells）のスフェアを、図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。
【００６１】
得られた各スフェアの位相差画像から、各スフェアの形状指標値として体積を導出した。更に、各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。横軸にスフェアの体積をとり、縦軸に位相差量総和θＡをとったグラフ上に、各スフェアについてプロットを形成することにより、スフェアの体積と位相差量総和θＡの相関特性を取得した。図８は、スフェアの体積と位相差量総和θＡとの相関特性を示すグラフである。図８に示すように、位相差量総和θＡとスフェアの体積とは、比例関係であることが確認された。図８には、スフェアの体積と位相差量総和θＡとの間の相関性の基準を示す基準相関トレンドラインＬＳ及びプロットが示されている。図８に示す各プロットから、導出した回帰線を、基準相関トレンドラインＬＳとして適用した。
【００６２】
図８には、基準相関トレンドラインＬＳ上に存在するプロットａ１及びａ２に対応するスフェアの位相差画像、及び基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ３、ａ４、ａ５に対応するスフェアの位相差画像が示されている。基準相関トレンドラインＬＳ上に存在するプロットａ１及びａ２に対応するスフェアについては、スフェアの全体に亘って輝度が均一な位相差画像が得られた。このことは、スフェアを構成する複数の細胞が均質であること、及びスフェア内における細胞の密度が均一であること等を示している。一方、基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ３及びａ４に対応するスフェアについては、中心部の輝度が他の部位と比較して低下している位相差画像が得られた。このことは、スフェアを構成する複数の細胞が不均質であること、及びスフェア内における細胞の密度が不均一であること等を示している。また、基準相関トレンドラインＬＳから乖離した位置に存在するプロットａ５に対応するスフェアについては、スフェアの輪郭線の凹凸が顕著な位相差画像が得られた。このことは、スフェアを構成する細胞に異常が生じていること等を示している。
【００６３】
以上の結果から、位相差量総和θＡと、スフェアの形状指標値の一例である体積との相関性を用いてスフェアの状態を判定することが可能であるといえる。また、位相差量総和θＡとスフェアの体積との相関性を示す基準相関トレンドラインＬＳと、判定対象のスフェアについての位相差量総和θＡと形状指標値との相関性とを比較し、判定対象のスフェアについての相関性の、基準相関トレンドラインＬＳからの乖離の程度に応じて当該スフェアの状態を判定することが可能であるといえる。具体的には、一方の軸にスフェアの体積をとり、他方の軸に位相差量総和θＡをとったグラフ上にプロットした、判定対象のスフェアについてのプロットの、基準相関トレンドラインＬＳからの乖離の程度に応じて、各スフェアの状態を判定することが可能であるといえる。従って、例えば、位相差量総和θＡの基準相関トレンドラインＬＳからのマイナス幅が閾値以上であるスフェアについては、当該スフェアに含まれる複数の細胞の密度、均質性、及びスフェアの外形形状のうちの少なくとも１つに異常があると判定することができる。
【００６４】
なお、本実施例では、スフェアの形状指標値としてスフェアの体積を用いたが、これに代えて、スフェアの断面積、粒径または周の長さを用いることも可能である。これらいずれの形状指標値を用いる場合でも、位相差量総和θＡとの相関性を用いてスフェアの状態を判定することが可能である。
【００６５】
［実施例２］
本実施例では、複数のスフェアを含む培養ロットを判定対象とした。判定対象の培養ロットは、下記の表１に示す良ロット及び不良ロットである。良ロットに属するスフェア及び不良ロットに属するスフェアは、それぞれ、三次元培養法により培養したｉＰＳ細胞のスフェアである。培養期間の初日から５日目までの細胞の増殖率は、良ロットに属するスフェアおいては２０倍であり、不良ロットに属するスフェアにおいては３．７倍であった。培養５日目における細胞の生存率は、良ロットに属するスフェアおいては９０．２％であり、不良ロットに属するスフェアにおいては６４．１％であった。
【表１】
2019176427.tif 000007
【００６６】
良ロット及び不良ロットの各々について、試料１４であるスフェアを図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。
【００６７】
得られた各スフェアの位相差画像から、各スフェアの粒径を導出し、スフェアの粒径のヒストグラムを取得した。図９Ａは、良ロット及び不良ロットの各々について導出したスフェアの粒径のヒストグラムである。図９Ａに示すように、良ロットと不良ロットとの間でヒストグラムに有意な差はみられなかった。このことは、スフェアの状態は、当該スフェアの粒径に反映されにくく、スフェアの粒径のみに基づいてスフェアの状態を判定することは困難であることを示している。
【００６８】
次に、良ロット及び不良ロットの各々について得られた各スフェアの位相差画像から、各スフェアの形状指標値として体積を導出し、更に、各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。横軸に、スフェアの体積をとり、縦軸に位相差量総和θＡをとったグラフ上に、各培養ロットの各スフェアについてプロットを形成することにより、スフェアの体積と位相差量総和θＡの相関特性を取得した。図９Ｂは、良ロット及び不良ロットの各々について取得した、スフェアの体積と位相差量総和θＡの相関特性を示すグラフである。図９Ｂには、基準相関トレンドラインＬＳ及びプロットが示されている。図９Ｂに示すように、良ロットにおいては、位相差量総和θＡが基準相関トレンドラインＬＳを下回るスフェアが少数であるのに対し、不良ロットにおいては、位相差量総和θＡが基準相関トレンドラインＬＳを下回るスフェアが、良ロットと比べて、多数存在した。
【００６９】
以上の結果から、位相差量総和θＡと、スフェアの形状指標値の一例である体積との相関性を用いて培養ロットの良否判定を行うことが可能であるといえる。また、位相差量総和θＡとスフェアの体積との相関性を示す基準相関トレンドラインＬＳと、判定対象の培養ロットに属するスフェアについての相関性とを比較し、判定対象の培養ロットについての相関性の、基準相関トレンドラインＬＳからの乖離の程度に応じて当該培養ロットの良否判定が可能であるといえる。具体的には、一方の軸にスフェアの体積をとり、他方の軸に位相差量総和θＡをとったグラフ上にプロットした判定対象の培養ロットに属する各スフェアについてのプロットの、基準相関トレンドラインＬＳからの乖離の状態（基準相関トレンドラインＬＳを下回っているか否か）及び乖離の程度の少なくとも一方に基づいて当該培養ロットの状態を判定することが可能であるといえる。従って、例えば、位相差量総和θＡの基準相関トレンドラインＬＳからのマイナス幅が第１の閾値以上であるスフェアの含有率が、第２の閾値以上である場合、当該培養ロットが不良ロットであると判定することができ、上記含有率が、第２の閾値未満である場合、当該培養ロットが良ロットであると判定することができる。
【００７０】
次に、良ロット及び不良ロットの各々について得られた各スフェアの位相差画像から、（４）式によって示される位相差量密度ＤＰを導出した。図９Ｃは、良ロット及び不良ロットの各々について取得した、位相差量密度ＤＰのヒストグラムである。図９Ｃに示すように、良ロットと不良ロットとの間で、ヒストグラムに有意な差が見られた。具体的には、良ロットにおいては、位相差量密度ＤＰが相対的に高い位置に分布のピークが存在し、不良ロットにおいては、位相差量密度ＤＰが相対的に低い位置に分布のピークが存在した。なお、位相差量密度ＤＰのヒストグラムにおける良ロットと不良ロットとの間の有意差の信頼性を示す指標であるｐ値（有意確率）は、０．０５未満であった。
【００７１】
以上の結果から、スフェアの位相差量密度ＤＰに基づいて、培養ロットの良否判定を行うことが可能であるといえる。従って、例えば、判定対象の培養ロットに属する複数のスフェアについて、位相差量密度ＤＰの平均値を求め、上記平均値が閾値以上である場合に、当該培養ロットを良ロットと判定することができ、上記平均値が閾値未満である場合に、当該培養ロットを不良ロットと判定することができる。
【００７２】
なお、本実施例では、位相差量密度ＤＰを培養ロットの良否判定に用いる場合を例示したが、位相差量密度ＤＰをスフェア単体の良否判定に用いることも可能である。例えば、判定対象のスフェアについて取得した位相差量密度ＤＰが閾値以上である場合に、当該スフェアの状態が良好であると判定することができ、位相差量密度ＤＰが閾値未満である場合に、当該スフェアの状態が異常であると判定することができる。
【００７３】
［実施例３］
本実施例では、互いに異なる複数の培養ロットに属する複数のスフェアを判定対象とした。判定対象のスフェアは、いずれも三次元培養法により培養したｉＰＳ細胞のスフェアである。
【００７４】
複数の培養ロットの各々について、試料１４であるスフェアを図１に示す撮像システム１のサンプルステージにセットし、複数のスフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。得られた各スフェアの位相差画像について、（４）式によって示される位相差量密度ＤＰを培養ロット毎に導出した。
【００７５】
次に、複数の培養ロットの各々について、細胞の生存率を取得した。生存率の測定手順は以下のとおりである。遠沈管に収容したスフェアを４００Ｇ（Ｇ：重力加速度）、３分間の遠心分離処理を実施した。遠沈管の底に沈殿したスフェアを回収し、細胞解離酵素であるTrypLE（登録商標） Selectを用いて単一細胞に分解した。分解した細胞に対してTrypan Blueによる死細胞染色処理を施した後、細胞を、市販されている一般的なセルカウンタ（Countess（登録商標）） を用いて計数することにより、染色された死細胞の個数をカウントした。死細胞のカウント数から細胞の生存率を培養ロット毎に導出した。
【００７６】
横軸に細胞の生存率をとり、縦軸に位相差量密度ＤＰをとったグラフ上に、各培養ロットについてプロットを形成することにより、細胞の生存率と位相差量密度ＤＰ（ロット平均値）との相関特性を取得した。図１０は、細胞の生存率と位相差量密度ＤＰ（ロット平均値）との相関特性を示すグラフである。図１０に示すように、細胞の生存率と、位相差量密度ＤＰとは比例関係であることが確認された。
【００７７】
以上の結果より、位相差量密度ＤＰから細胞の生存率を推定できるといえる。位相差量密度ＤＰから細胞の生存率を推定することにより、遠心分離処理及び染色処理等の細胞の破壊を伴う処理を行うことなく（すなわち非破壊で）、細胞の生存率を把握することができる。例えば、細胞の生存率と位相差量密度ＤＰとの相関性を示す相関ラインを予め取得しておき、判定対象の培養ロットに属するスフェアについて取得した位相差量密度ＤＰの平均値と、上記の相関ラインとから当該培養ロットにおける細胞の生存率を推定することが可能である。
【００７８】
なお、本実施例では、位相差量密度ＤＰから培養ロット内における細胞の生存率を推定する場合を例示したが、位相差量密度ＤＰから単一のスフェア内における細胞の生存率を推定することも可能である。例えば、細胞の生存率と位相差量密度ＤＰとの相関性を示す相関ラインを予め取得しておき、判定対象のスフェアについて取得した位相差量密度ＤＰと、上記の相関ラインとから当該スフェア内における細胞の生存率を推定することが可能である。
【００７９】
［実施例４］
細胞の生存率が、位相差画像に与える影響について確認した。二次元培養法（接着培養法）により培養したｉＰＳ細胞をネクローシスに誘導するべく、細胞のコロニーに３ｗｔ％のＨ２Ｏ２（過酸化水素）を添加した後、蛍光色素であるSYTOX（登録商標） Greenによる染色処理を施した。続いて、上記の処理を施したコロニーを図１に示す撮像システム１のサンプルステージにセットし、コロニーのホログラムをＣＭＯＳカメラ１９によって撮影した。取得したホログラムに対してコンピュータによる数値計算を実施することにより、コロニーの位相差画像を取得した。ホログラム撮影と並行して、蛍光顕微鏡画像をＣＭＯＳカメラ３６によって撮影した。
【００８０】
図１１Ａ、図１１Ｂ、及び図１１Ｃは、それぞれ、Ｈ２Ｏ２を添加した後、０分経過した時点、１５分経過した時点、３７経過した時点におけるコロニーの蛍光顕微鏡画像である。図１２Ａ、図１２Ｂ、及び図１２Ｃは、それぞれ、Ｈ２Ｏ２を添加した後、０分経過した時点、１５分経過した時点、３７分経過した時点におけるコロニーの位相差画像である。
【００８１】
図１１Ａ〜図１１Ｃに示すように、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、発光箇所が増加した。このことは、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、死細胞が増加したこと（換言すれば、細胞の生存率が低下したこと）を意味する。また、図１２Ａ〜図１２Ｃに示すように、Ｈ２Ｏ２を添加した時点からの時間経過に伴って、位相差画像における低輝度領域の面積が増加した。このことは、細胞の生存率の低下に伴って、位相差画像における位相差量総和θＡが低下することを意味する。このように、細胞の生存率が、位相差量総和θＡに影響を及ぼすことが確認できた。
【００８２】
［実施例５］
二次元培養法（接着培養法）により培養したｉＰＳ細胞のコロニーにＨ２Ｏ２を添加することにより、ネクローシスに誘導した。また、別のｉＰＳ細胞のコロニーについて、PromoKine社のApoptosis Inducer Setを用いてアポトーシスに誘導した。なお、アポトーシスは、多細胞生物の細胞において増殖制御機構として管理・調節された能動的な細胞死である。一方、ネクローシスは、栄養不足、毒物、外傷などの外的環境要因により起こる受動的な細胞死であり、細胞死に至る過程がアポトーシスと異なる。
【００８３】
上記の処理を施したコロニー及び健全なコロニーを図１に示す撮像システム１のサンプルステージにセットし、各コロニーのホログラムをＣＭＯＳカメラ１９によって撮影した。取得したホログラムに対してコンピュータによる数値計算を実施することにより、各コロニーの位相差画像を取得した。各コロニーの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。
【００８４】
図１３は、ネクローシスに誘導したコロニー（死細胞）、アポトーシスに誘導したコロニー（死細胞）及び健全なコロニー（生細胞）のそれぞれの位相差量総和θＡの相対値を示すグラフである。なお、健全なコロニー（生細胞）の位相差量総和θＡを１とした。図１３に示すように、健全なコロニー（生細胞）と、ネクローシスに誘導したコロニー（死細胞）及びアポトーシスに誘導したコロニー（死細胞）との間で、位相差量総和θＡに有意な差が見られた。有意差の信頼性を示す指標であるｐ値（有意確率）は、０．０５未満であった。
【００８５】
以上の結果から細胞死に至る過程が、ネクローシスであるかアポトーシスであるかにかかわらず、位相差量総和θＡを、細胞の生死の判別に利用できるといえる。例えば、判定対象のスフェアについて取得した位相差量総和θＡが、閾値以上である場合、当該スフェアは、生細胞をより多く含むと判定することができ、判定対象のスフェアについて取得した位相差量総和θＡが、閾値未満である場合、当該スフェアは、死細胞をより多く含むと判定することができる。なお、位相差量総和θＡに代えて位相差量密度ＤＰを、細胞の生死の判別に利用することも可能である。
【００８６】
［実施例６］
未分化状態を維持したｉＰＳ細胞のスフェア及び未分化状態を逸脱したｉＰＳ細胞のスフェアを、図１に示す撮像システム１のサンプルステージにセットし、各スフェアのホログラムをＣＭＯＳカメラ１９によって撮影した。取得した各スフェアのホログラムに対してコンピュータによる数値計算を実施することにより、スフェアの中央付近においてスライスした位相差画像を取得した。各スフェアの位相差画像について、（３）式によって示される位相差量総和θＡを導出した。なお、未分化状態を逸脱したスフェアについては、未分化状態の維持に必要な成分を添加しない基礎培地を用いることにより未分化逸脱状態に誘導した。
【００８７】
図１４は、未分化状態を維持したスフェア及び未分化状態を逸脱したスフェアのそれぞれの位相差量総和θＡの平均値を示すグラフである。図１４に示すように、未分化状態を維持したスフェアと、未分化状態を逸脱したスフェアとの間で、位相差量総和θＡに有意な差が見られた。
【００８８】
以上の結果から、位相差量総和θＡを、幹細胞における未分化状態を維持しているか、未分化状態を逸脱しているかの判定に利用できるといえる。例えば、判定対象のスフェアについて取得した位相差量総和θＡが、閾値以上である場合に当該スフェアが未分化状態を維持していると判定することができ、判定対象のスフェアについて取得した位相差量総和θＡが閾値未満である場合に当該スフェアが未分化状態を逸脱したと判定することができる。なお、位相差量総和θＡに代えて位相差量密度ＤＰを、未分化状態を維持しているか、未分化状態を逸脱しているかの判定に利用することも可能である。
【００８９】
以上の説明から明らかなように、開示の技術の実施形態に係る判定方法によれば、スフェアのホログラムから生成された位相差画像と当該スフェアの形状に応じた形状指標値とに基づいて、スフェアの状態が判定されるので、スフェアの状態を非破壊且つ簡便に判定することが可能である。
【００９０】
すなわち、開示の技術の実施形態に係る判定方法によれば、従来の光学顕微鏡では把握し得ないスフェアの内部における細胞の密度及び生存率などについても判定することができる。また、本判定方法において実施されるホログラムの撮影、画像解析及び統計処理などについては、コンピュータによる自動化が容易であり、従来の手法と比較して、人手及び処理時間を削減することができる。また、作業者による手技のばらつきが影響することもなく、安定した判定結果を得ることができる。また、スフェアを単一細胞に分解したり、蛍光色素を添加したりするといった細胞の破壊を伴う処理を行うことなく、スフェアの状態を非破壊で判定することが可能である。また、ホログラム撮影は、三次元空間中にランダムに存在するスフェアについても容易に行うことが可能であり、三次元培養法による培養中のスフェアの状態を、その場で判定することも可能である。
【００９１】
また、スフェアの状態の判定において、球形状のスフェアの中央付近においてスライスした位相差画像を用いることにより、より正確な判定結果を得ることができる。球形状のスフェアの中央付近においてスライスした１枚の位相差画像を用いることにより、スライス位置が互いに異なる複数の位相差画像を用いる場合と比較して、判定に要する処理を簡略化することができる。
【００９２】
開示の技術に係る判定方法は、複数の細胞の凝集体であるスフェアを撮像したホログラムからスフェアの位相差画像を生成し、位相差画像とスフェアの形状に応じた形状指標値とに基づいて、スフェアの状態を判定することを含む。
【００９３】
開示の技術に係る判定方法によれば、スフェアの状態を非破壊且つ簡便に判定することが可能となる。
【００９４】
開示の技術に係る判定方法において、位相差画像と形状指標値とに基づいて、スフェアに含まれる複数の細胞の生存率、密度、均質性、未分化状態逸脱、及びスフェアの外形形状の少なくとも１つに関する判定を行ってもよい。
【００９５】
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を導出し、位相差量総和と形状指標値との相関性を用いてスフェアの状態を判定してもよい。
【００９６】
スフェアの状態の判定に位相差量総和と形状指標値との相関性を用いることにより、スフェアの状態について適確な判定を行うことが可能となる。
【００９７】
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象のスフェアについての位相差量総和と形状指標値との相関性とを比較し、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の程度に応じて判定対象のスフェアの状態を判定してもよい。例えば、判定対象のスフェアについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の幅が閾値を超える場合、判定対象のスフェアに含まれる細胞の生存率、密度、均質性及びスフェアの外形形状の少なくとも１つについて異常ありと判定してもよい。
【００９８】
この態様によれば、スフェアの状態判定の適確性を更に高めることが可能となる。
【００９９】
開示の技術に係る判定方法において、位相差量総和と形状指標値との相関性の基準を示す基準相関トレンドラインと、判定対象の培養ロットに属する複数のスフェアの各々についての位相差量総和と形状指標値との相関性とを比較し、判定対象の培養ロットについての位相差量総和と形状指標値との相関性の、基準相関トレンドラインからの乖離の状態及び乖離の程度の少なくとも一方に応じて判定対象の培養ロットの良否判定を行ってもよい。
【０１００】
この態様によれば、複数のスフェアを含む培養ロットについての良否判定を適確に行うことが可能となる。
【０１０１】
開示の技術に係る判定方法において、形状指標値は、スフェアの体積、断面積、粒径及び周の長さのいずれかであってもよい。
【０１０２】
開示の技術に係る判定方法において、位相差画像を構成する複数の画素の各々の位相差量を積算した値である位相差量総和を、スフェアの体積で除算して位相差量密度を導出し、位相差量密度に基づいて、スフェアの状態を判定してもよい。
【０１０３】
スフェアの状態を判定に位相差量密度を用いることにより、スフェアの状態について適確な判定を行うことが可能となる。
【０１０４】
開示の技術に係る判定方法において、スフェアの状態の判定に用いる位相差画像は、ホログラムから生成可能な複数の位相差画像のうち、位相差画像を構成する複数の画素間の位相差量のばらつきが最大となる位相差画像であることが好ましい。
【０１０５】
スフェアの状態の判定に用いる位相差画像として、上記の画像を用いることで、当該位相差画像から、スフェアの実態に合致したより正確な情報を得ることができ、スフェアの状態判定の適確性を更に高めることが可能となる。
【符号の説明】
【０１０６】
１  撮像システム
１０  ホログラム光学系
１１  レーザ光源
１２  ビームスプリッタ
１３  コリメートレンズ
１４  試料
１５  対物レンズ
１７  結像レンズ
１８  ビームスプリッタ
１９、３６  ＣＭＯＳカメラ
２０  光ファイバ
２１  コリメートレンズ
３０  蛍光顕微鏡光学系
３１  励起光源
３２  励起用フィルタ
３３  ＯＮ／ＯＦＦ切り替えミラー
３４  ダイクロイックミラー
３５  分光フィルタ
５００  コンピュータ
５０２  主記憶装置
５０３  補助記憶装置
５０４  通信インターフェース
５０５  表示部
５０６  オートフォーカスプログラム
５０７  バス
ＤＰ  位相差量密度
ＩＰ  位相差画像
ＬＳ  基準相関トレンドライン
θ  位相差量
θＢ  バックグランドの位相
θＳ  スフェアの存在する領域の位相
θＡ  位相差量総和
θｋ  １画素の位相差量
Ｖ  スフェアの体積
ｖｋ  位相差画像の各画素ｋに対応する部位におけるスフェアの体積
ａ１、ａ２、ａ３、ａ４、ａ５  プロット
ｋ  画素
ｗ  カーブの幅
【図１】
2019176427.tif 000008
【図２Ａ】
2019176427.tif 000009
【図２Ｂ】
2019176427.tif 000010
【図２Ｃ】
2019176427.tif 000011
【図２Ｄ】
2019176427.tif 000012
【図３】
2019176427.tif 000013
【図４】
2019176427.tif 000014
【図５】
2019176427.tif 000015
【図６】
2019176427.tif 000016
【図７】
2019176427.tif 000017
【図８】
2019176427.tif 000018
【図９Ａ】
2019176427.tif 000019
【図９Ｂ】
2019176427.tif 000020
【図９Ｃ】
2019176427.tif 000021
【図１０】
2019176427.tif 000022
【図１１Ａ】
2019176427.tif 000023
【図１１Ｂ】
2019176427.tif 000024
【図１１Ｃ】
2019176427.tif 000025
【図１２Ａ】
2019176427.tif 000026
【図１２Ｂ】
2019176427.tif 000027
【図１２Ｃ】
2019176427.tif 000028
【図１３】
2019176427.tif 000029
【図１４】
2019176427.tif 000030
【国際調査報告】
2019176427.tif 000031
2019176427.tif 000032
2019176427.tif 000033
2019176427.tif 000034