JPWO2020071484A1 - 化学物質のスクリーニング方法、プログラム、制御装置、及び培養観察装置 - Google Patents

Abstract

化学物質のスクリーニング方法は、異なる化学物質がそれぞれ異なる濃度で添加された複数の細胞群における細胞の顕微鏡画像に含まれる前記細胞の形態に関する複数の特徴量を取得することと、前記細胞の形状に関する複数の特徴量に基づいて、前記異なる化学物質同士の類似性を判定することと、を有する。

Description

本発明は、化学物質のスクリーニング方法、プログラム、制御装置、及び培養観察装置に関するものである。
本願は、２０１８年１０月５日に、日本に出願された特願２０１８−１８９８１０号に基づき優先権を主張し、その内容をここに援用する。

従来、細胞を用いた化学物質のハイスループットスクリーニングには、単純な系が用いられていた。例えば、抗癌作用を示す低分子化合物のスクリーニング方法の多くには、ＭＴＴ（Ｍｅｔｈｙｌ ｔｈｉａｚｏｌｙｌ ｔｅｔｒａｚｏｒｉｕｍ）アッセイやａｌａｍａｒＢｌｕｅアッセイといった１つのウェル内における細胞の生存数を定量し、生存率を算出する方法が用いられていた。

近年、創薬ターゲットのスクリーニング方法や低分子化合物のスクリーニング方法は、機器やイメージング技術の発展と共に複雑化している。顕微鏡技術がスクリーニングに応用されるようになり、１種類の細胞から複数のパラメータについて解析が行われるようになってきている。このようなスクリーニング方法をＨｉｇｈ ｃｏｎｔｅｎｔ ｓｃｒｅｅｎｉｎｇ（以下、ＨＣＳ）という（非特許文献１参照）。
しかしながら、従来の化学物質のスクリーニング方法では、コストと時間がかかってしまう。

Ｋｒａｕｓｚ、Ｍｏｌ．ＢｉｏＳｙｓｔ．、第３巻、第２３２〜２４０頁、２００７年

上記問題を解決するために、本発明の一態様は、異なる化学物質がそれぞれ異なる濃度で添加された複数の細胞群における細胞の顕微鏡画像に含まれる前記細胞の形態に関する複数の特徴量を取得することと、前記細胞の形状に関する複数の特徴量に基づいて、前記異なる化学物質同士の類似性を判定することと、を有する化学物質のスクリーニング方法である。

上記問題を解決するために、本発明の一態様は、コンピュータに、上記の化学物質のスクリーニング方法に基づくステップを実行させるプログラムである。

上記問題を解決するために、本発明の一態様は、異なる化学物質がそれぞれ異なる濃度で添加された複数の細胞群における細胞の顕微鏡画像に含まれる前記細胞の形態に関する複数の特徴量を算出する特徴量算出部と、前記細胞の形態に関する複数の特徴量に基づいて、前記異なる化学物質同士の類似性を判定する類似性判定部と、を備える制御装置である。

上記問題を解決するために、本発明の一態様は、上記の制御装置を備える培養観察装置である。

本発明に用いられる培養観察装置の概要を示すブロック図である。 本発明に用いられる培養観察装置の正面図である。 本発明に用いられる培養観察装置の平面図である。 本発明に用いられるンキュベータでの観察動作の一例を示すフローチャートである。 本発明の第１実施形態の撮像装置が撮像する顕微鏡画像の一例を示す模式図である。 本発明の第１実施形態の制御装置による化学物質のスクリーニングの動作の一例を示すフローチャートである。 本発明の第１実施形態の主成分分析により処理された特徴量点の一例を示す図である。 本発明の第１実施形態の濃度軌跡の一例を示す図である。 本発明の第１実施形態に係るプロットのコントロール領域の中心からの距離と化学物質濃度との関係の一例を示す図である。 本発明の第１実施形態に係る化学物質ごとの形態変化領域の開始濃度の一例を示す図である。 本発明の第１実施形態に係るプロットの死細胞領域の中心からの距離と化学物質濃度との関係の一例を示す図である。 本発明の第１実施形態に係る化学物質ごとの形態変化領域の一例を示す図である。 本発明の第１実施形態に係る制御装置による化学物質同士の特徴量空間における特徴量点間の比較の動作の一例を示すフローチャートである。 本発明の第１実施形態に係る制御装置による化学物質同士の特徴量空間における特徴量点間の比較の方法の一例を示す図である。 本発明の第２実施形態に用いられる培養観察装置１１ａの概要を示すブロック図である。 本発明の第２実施形態に用いられる制御装置による既知の化学物質の解析及び新規化学物質の解析の動作の一例を示すフローチャートである。 本発明の第２実施形態に係る高濃度側の濃度軌跡及び５０％阻害濃度における時間軌跡の一例を示す図である。 本発明の第２実施形態に係る低濃度側の濃度軌跡及び５０％阻害濃度における時間軌跡の一例を示す図である。 本発明の第２実施形態に係る第１の化学物質の５０％阻害濃度における時間軌跡の一例を示す図である。 本発明の第２実施形態に係る第１の化学物質の低濃度側の線分の和と高濃度側の線分の和との時間変化の一例を示す図である。 本発明の第２実施形態に係る第２の化学物質の５０％阻害濃度における時間軌跡の一例を示す図である。 本発明の第２実施形態に係る第２の化学物質の低濃度側の線分の和と高濃度側の線分の和との時間変化の一例を示す図である。 本発明の第２実施形態に係る制御装置による新規化学物質の判定の動作の一例を示すフローチャートである。 本発明の第２実施形態に係る新規化学物質と第１の化学物質との比較の一例を示す図である。 本発明の第２実施形態に係る新規化学物質と第２の化学物質との比較の一例を示す図である。

（第１実施形態）
≪培養状態評価装置≫
本発明の化学物質のスクリーニング方法には、細胞群に対して、この細胞群に含まれる各細胞に化学物質を添加する第１の工程と、化学物質を添加した細胞の形態変化を観察する第２の工程とが含まれる。ここで細胞群とは、複数の細胞を意味する。第１の工程においては、互いに異なる複数の濃度の化学物質を、細胞に添加する。第２の工程においては、濃度の差による細胞の形態変化の差を観察する。この第２の工程は、培養状態評価装置を用いて行われる。この培養状態評価装置は、培養観察装置に備えられている。培養状態評価装置とは、培養観察装置１１の制御装置４１である。
以下、この培養観察装置について説明する。

図１は、本発明に用いられる培養観察装置の概要を示すブロック図である。また、図２、図３は、本発明に用いられる培養観察装置の正面図および平面図である。

本発明に用いられる培養観察装置１１は、上部ケーシング１２と下部ケーシング１３と、表示部５０とを有している。培養観察装置１１の組立状態において、上部ケーシング１２は下部ケーシング１３の上に載置される。なお、上部ケーシング１２と下部ケーシング１３との内部空間は、ベースプレート１４によって上下に仕切られている。

まず、上部ケーシング１２の構成の概要を説明する。上部ケーシング１２の内部には、細胞の培養を行う恒温室１５が形成されている。この恒温室１５は温度調整装置１５ａおよび湿度調整装置１５ｂを有しており、恒温室１５内は細胞の培養に適した環境（例えば温度３７℃、湿度９０％の雰囲気）に維持されている（なお、図２、図３での温度調整装置１５ａ、湿度調整装置１５ｂの図示は省略する）。

恒温室１５の前面には、大扉１６、中扉１７、小扉１８が配置されている。大扉１６は、上部ケーシング１２および下部ケーシング１３の前面を覆っている。中扉１７は、上部ケーシング１２の前面を覆っており、大扉１６の開放時に恒温室１５と外部との環境を隔離する。小扉１８は、細胞を培養する培養容器１９を搬出入するための扉であって、中扉１７に取り付けられている。この小扉１８から培養容器１９を搬出入することで、恒温室１５の環境変化を抑制することが可能となる。なお、大扉１６、中扉１７、小扉１８は、パッキンＰ１，Ｐ２，Ｐ３によりそれぞれ気密性が維持されている。

また、恒温室１５には、ストッカー２１、観察ユニット２２、容器搬送装置２３、搬送台２４が配置されている。ここで、搬送台２４は、小扉１８の手前に配置されており、培養容器１９を小扉１８から搬出入する。

ストッカー２１は、上部ケーシング１２の前面（図３の下側）からみて恒温室１５の左側に配置される。ストッカー２１は複数の棚を有しており、ストッカー２１の各々の棚には培養容器１９を複数収納することができる。なお、各々の培養容器１９には、培養の対象となる細胞が培地とともに収容されている。

観察ユニット２２は、上部ケーシング１２の前面からみて恒温室１５の右側に配置される。この観察ユニット２２は、培養容器１９内の細胞のタイムラプス観察を実行することができる。

ここで、観察ユニット２２は、上部ケーシング１２のベースプレート１４の開口部に嵌め込まれて配置される。観察ユニット２２は、試料台３１と、試料台３１の上方に張り出したスタンドアーム３２と、位相差観察用の顕微光学系２２３および撮像装置２２４を内蔵した本体部分３３とを有している。そして、試料台３１およびスタンドアーム３２は恒温室１５に配置される一方で、本体部分３３は下部ケーシング１３内に収納される。

試料台３１は透光性の材質で構成されており、その上に培養容器１９を載置することができる。この試料台３１は水平方向に移動可能に構成されており、上面に載置した培養容器１９の位置を調整できる。また、スタンドアーム３２にはＬＥＤ光源２２１が内蔵されている。スタンドアーム３２には照明リング絞りが設けられ、試料台３１に照射されるＬＥＤ光源２２１からの照明光の光強度分布を可変に調整できる。スタンドアーム３２は照明装置２２２として機能する。照明装置２２２と、顕微光学系２２３とは、位相差顕微鏡を構成する。顕微光学系２２３には、コンデンサレンズや対物レンズや位相板などが設けられ、公知の位相差顕微鏡の顕微光学系と同様に構成される。

すなわち、ＬＥＤ光源２２１が照射する照明光は、スタンドアーム３２（照明装置２２２）に設けられた照明リング絞りによって絞られる。照明リング絞りによって絞られた照明光は、コンデンサレンズを通り、培養容器１９の細胞に到達する。培養容器１９の細胞に到達した照明光は、直接光と回折光とに分離する。ここで直接光とは、培養容器１９の細胞の内部を直進する照明光である。回折光とは、位相物体である細胞により回折した照明光である。回折現象は屈折率に違いのある部位において発生するので、回折光は細胞と培養容器１９内において細胞が浸される溶液との境界部分、細胞の内部構造部など、位相物体である細胞の形状の情報を含む。

直接光と回折光とは、対物レンズを通り位相板に到達する。位相板は、４分の１波長板と、ＮＤフィルターとがリング状に形成されており、４分の１波長板及びＮＤフィルターと以外の部分は透明である。ここで４分の１波長板は、光の位相を４分の１波長分だけずらす。ＮＤフィルターは、光を吸収する。
直接光は対物レンズによって集光され、位相板のリング状の４分の１波長板を通過するため位相が４分の１波長分だけずれる。また、４分の１波長板を通過する直接光は明るさが弱められる。一方、回折光は大部分が位相板の透明な部分を通過し、位相及び明るさは変化しない。
直接光と回折光とは顕微光学系２２３の像面に到達して、明暗のコントラストがついた像を結ぶ。

撮像装置２２４は、スタンドアーム３２によって試料台３１の上側から透過照明された培養容器１９の細胞を、顕微光学系２２３を介して撮像することで細胞の顕微鏡画像を取得できる。

容器搬送装置２３は、上部ケーシング１２の前面からみて恒温室１５の中央に配置される。この容器搬送装置２３は、ストッカー２１、観察ユニット２２の試料台３１および搬送台２４との間で培養容器１９の受け渡しを行う。

図３に示すように、容器搬送装置２３は、多関節アームを有する垂直ロボット３４と、回転ステージ３５と、ミニステージ３６と、アーム部３７とを有している。回転ステージ３５は、垂直ロボット３４の先端部に回転軸３５ａを介して水平方向に１８０°回転可能に取り付けられている。そのため、回転ステージ３５は、ストッカー２１、試料台３１および搬送台２４に対して、アーム部３７をそれぞれ対向させることができる。

また、ミニステージ３６は、回転ステージ３５に対して水平方向に摺動可能に取り付けられている。ミニステージ３６には培養容器１９を把持するアーム部３７が取り付けられている。アーム部３７は、把持した培養容器１９をストッカー２１から搬出して観察ユニット２２の試料台３１に載置する。

次に、下部ケーシング１３の構成の概要を説明する。下部ケーシング１３の内部には、観察ユニット２２の本体部分３３や、培養観察装置１１の制御装置４１が収納されている。

制御装置４１は、温度調整装置１５ａ、湿度調整装置１５ｂ、観察ユニット２２および容器搬送装置２３とそれぞれ接続されている。この制御装置４１は、所定のプログラムに従って培養観察装置１１の各部を統括的に制御する。

一例として、制御装置４１は、温度調整装置１５ａおよび湿度調整装置１５ｂをそれぞれ制御して恒温室１５内を所定の環境条件に維持する。また、制御装置４１は、所定の観察スケジュールに基づいて、観察ユニット２２および容器搬送装置２３を制御して、培養容器１９の観察シーケンスを自動的に実行する。さらに、制御装置４１（培養状態評価装置）は、観察シーケンスで取得した画像に基づいて、細胞の培養状態の評価を行う培養状態評価処理を実行する。

ここで、制御装置４１（培養状態評価装置）は、ＣＰＵ４２および記憶部４３を有している。ＣＰＵ４２は、制御装置４１（培養状態評価装置）の各種の演算処理を実行するプロセッサである。ＣＰＵ４２は、一例としてＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｉｏｎ）などの集積回路で構成される。なお、ＣＰＵ４２の一部、または全部を、ソフトウェアとして構成してもよい。つまり、ＣＰＵ４２は、ハードウェアとソフトウェアとを組み合わせて構成されてもよい。また、ＣＰＵ４２は、プログラムの実行によって、特徴量算出部４４、座標算出部４５、濃度軌跡算出部４６、変化領域判定部４７、及び類似性判定部４８としてそれぞれ機能する。なお、特徴量算出部４４、座標算出部４５、濃度軌跡算出部４６、変化領域判定部４７、及び類似性判定部４８の動作については後述する。

記憶部４３は、ハードディスクや、フラッシュメモリ等の不揮発性の記憶媒体などで構成される。この記憶部４３には、ストッカー２１に収納されている各培養容器１９に関する管理データと、撮像装置で撮像された顕微鏡画像のデータと、既知の化学物質に関するデータとが記憶されている。さらに、記憶部４３には、ＣＰＵ４２によって実行されるプログラムが記憶されている。

なお、上記の管理データには、（ａ）個々の培養容器１９を示すインデックスデータ、（ｂ）ストッカー２１での培養容器１９の収納位置、（ｃ）培養容器１９の種類および形状（ウェルプレート、ディッシュ、フラスコなど）、（ｄ）培養容器１９で培養されている細胞の種類、（ｅ）培養容器１９の観察スケジュール、（ｆ）タイムラプス観察時の撮像条件（対物レンズの倍率、容器内の観察地点等）、などが含まれている。また、ウェルプレートのように複数の小容器で同時に細胞を培養できる培養容器１９については、各々の小容器ごとにそれぞれ管理データが生成される。

表示部５０は、液晶ディスプレイなどを備えており、制御装置４１（培養状態評価装置）が出力する画像を表示する。表示部５０が表示する画像には、制御装置４１（培養状態評価装置）の濃度軌跡算出部４６が算出した軌跡の画像が含まれる。また、表示部５０が表示する画像には、制御装置４１（培養状態評価装置）の類似性判定部４８が判定した結果が含まれる。

≪培養状態評価装置における観察動作の例≫
次に、図４のフローチャートを参照しつつ、本発明に用いられる培養観察装置１１での観察動作の一例を説明する。この図４は、恒温室１５内に搬入された培養容器１９を、登録された観察スケジュールに従ってタイムラプス観察する動作例を示している。

ステップＳ１０１：ＣＰＵ４２は、記憶部４３の管理データの観察スケジュールと現在日時とを比較して、培養容器１９の観察開始時間が到来したか否かを判定する。観察開始時間となった場合（ステップＳ１０１；ＹＥＳ）、ＣＰＵ４２はステップＳ１０２に処理を移行させる。一方、培養容器１９の観察時間ではない場合（ステップＳ１０１；ＮＯ）には、ＣＰＵ４２は次の観察スケジュールの時刻まで待機する。

ステップＳ１０２：ＣＰＵ４２は、観察スケジュールに対応する培養容器１９の搬送を容器搬送装置２３に指示する。そして、容器搬送装置２３は、指示された培養容器１９をストッカー２１から搬出して観察ユニット２２の試料台３１に載置する。なお、培養容器１９が試料台３１に載置された段階で、スタンドアーム３２に内蔵されたバードビューカメラ（不図示）によって培養容器１９の全体観察画像が撮像される。

ステップＳ１０３：ＣＰＵ４２は、観察ユニット２２に対して細胞の顕微鏡画像の撮像を指示する。観察ユニット２２は、ＬＥＤ光源２２１を点灯させて培養容器１９を照明するとともに、撮像装置２２４を駆動させて培養容器１９内の細胞の顕微鏡画像を撮像する。撮像装置２２４が撮像する顕微鏡画像の一例を図５に示す。

図５は、撮像装置２２４が撮像する顕微鏡画像の一例を示す模式図である。この図５の一例においては、顕微鏡画像ＰＩＣには、細胞Ｃｅｌｌ１、細胞Ｃｅｌｌ２、細胞Ｃｅｌｌ３の画像が含まれている。この顕微鏡画像とは、位相差顕微鏡によって撮像される位相差画像である。この顕微鏡画像は、培養容器１９中に培養されている細胞の像を、この細胞を透過する光の細胞の部位ごとの位相差を光の明暗差に変換することによって示す。

図４に戻り、撮像装置２２４は、記憶部４３に記憶されている管理データに基づいて、ユーザーの指定した撮像条件（対物レンズの倍率、容器内の観察地点）に基づいて顕微鏡画像を撮像する。例えば、培養容器１９内の複数の観察地点を観察する場合、観察ユニット２２は、試料台３１の駆動によって培養容器１９の位置を逐次調整し、各々の観察地点でそれぞれ顕微鏡画像を撮像する。なお、ステップＳ１０３で取得された顕微鏡画像のデータは、制御装置４１（培養状態評価装置）に読み込まれるとともに、ＣＰＵ４２の制御によって記憶部４３に記録される。

ステップＳ１０４：ＣＰＵ４２は、観察スケジュールの終了後に培養容器１９の搬送を容器搬送装置２３に指示する。そして、容器搬送装置２３は、指示された培養容器１９を観察ユニット２２の試料台３１からストッカー２１の所定の収納位置に搬送する。その後、ＣＰＵ４２は、観察シーケンスを終了してステップＳ１０１に処理を戻す。以上で、図４のフローチャートの説明を終了する。

≪培養状態評価装置における培養状態評価処理≫
次に、培養状態評価装置における培養状態評価処理の一例を説明する。この一例では、培養容器１９をタイムラプス観察して取得した複数の顕微鏡画像を用いて、培養容器１９の培養細胞における細胞の特徴量を算出する例を説明する。

この一例において、培養状態評価装置には、１種類の細胞について、３種類の化学物質が、それぞれ８段階の濃度にされて添加された、２４個の培養容器１９が用意されている。つまり、この一例において、培養状態評価装置は、３種類の化学物質をスクリーニングする。３種類の化学物質のうち、１種類は作用機序などの性質が未知の化学物質（以下、未知の化学物質と称する）であり、２種類は作用機序などの性質が既知の化学物質（以下、既知の化学物質と称する）である。
本実施形態において、化学物質のスクリーニングとは、複数の化学物質のなかから目的に適合する化学物質を選別することである。ここで目的に適合するとは、例えば、特定の薬効を示すことや、特定の毒性を示すことである。本実施形態においては、化学物質のスクリーニングの一例として、化学物質同士の類似性を判定することにより化学物質を選別する。ここで類似性の判定に用いられる化学物質には、一例として、未知の化学物質と、既知の化学物質とが含まれる。つまり、本実施形態における化学物質のスクリーニングでは、一例として、未知の化学物質と、既知の化学物質との類似性を判定することにより未知の化学物質を選別する。

この培養状態評価処理において、制御装置４１（培養状態評価装置）は、多変量解析を用いて、細胞の形状に関する複数の特徴量の種類毎に座標軸を定めた空間において、顕微鏡画像から算出される特徴量を示す点を、細胞に添加された化学物質の濃度ごとに求める。ここで、特徴量とは、細胞の形状に関する特徴を示す量である。そして、制御装置４１（培養状態評価装置）は、細胞の形状に関する複数の特徴量の種類毎に座標軸を定めた空間において、特徴量を示す点からなるグラフを生成する。本実施形態の化学物質のスクリーニング方法は、異なる化学物質がそれぞれ異なる濃度で添加された細胞の顕微鏡画像に含まれる細胞の形状に関する複数の特徴量に基づいて、異なる化学物質同士の類似性を判定することを有する。

≪計算モデルの生成処理の例≫
以下、図６のフローチャートを参照しつつ、化学物質のスクリーニングの動作の一例を説明する。この図６の例では、細胞群を培養した培養容器１９を培養観察装置１１によって同一視野を同一の撮像条件でタイムラプス観察し、細胞群がタイムラプス観察によって撮像された複数の顕微鏡画像を予め用意する。以下においては、一例として、１種類の細胞について観察する場合について説明する。培養容器１９は、１種類の細胞について２４個用意される。培養容器１９内の細胞には、培養容器１９ごとに種類及び濃度がそれぞれ異なる化学物質が、添加されている。この図６の例では、１回のタイムラプス観察において、化学物質の種類及び濃度がそれぞれ異なる２４個の培養容器１９について、それぞれ１枚の顕微鏡画像が撮像される。なお、培養容器１９毎に撮像する顕微鏡画像の枚数は１枚に限られず、複数枚であってもよい。図６の例では、１回のタイムラプス観察において、２４枚の顕微鏡画像が撮像される。なお、図６の例において、撮像する顕微鏡画像の枚数は、培養容器１９の区画の数や、当該区画毎に設定される培養観察装置１１による観察ポイントの数（ＲＯＩの数）等に基づいて設定してもよい。
つまり、顕微鏡画像とは、異なる化学物質がそれぞれ異なる濃度で添加された細胞の位相差画像である。また、細胞に添加される物質には複数の種類があり、顕微鏡画像とは、物質の種類ごとに撮像された画像である。顕微鏡画像は、物質が細胞に添加されてから経過した時間に基づく複数の時刻において撮像される。
また、添加される化学物質の種類及び濃度がそれぞれ異なる培養容器１９を用いなくてもよく、化学物質の類似性の判定において再現性が得られる（確認できる）ように、添加される化学物質の種類及び／又は濃度が同様の培養容器１９を用いてもよい。また、なお、細胞に添加される物質は、異なる複数の種類や異なる複数の濃度の他、添加する条件を変化させて添加してもよい。

この一例において、細胞とは、肺がん細胞、腎臓がん細胞、皮膚がん細胞、前立腺がん細胞、大腸がん、脳腫瘍のうち、いずれか１種類の細胞である。なお、本実施形態においては、細胞が１種類である場合の一例について説明するが、複数の種類の細胞が用いられて化学物質のスクリーニングが行われてもよい。複数の種類の細胞には、がん細胞だけでなく正常な細胞が含まれてもよい。

また、この一例においては、化学物質とは、Ｐａｃｌｉｔａｘｅｌ、Ｅｔｏｐｏｓｉｄｅ、及び未知の化学物質Ｘの３種類の薬剤である。このＰａｃｌｉｔａｘｅｌは、細胞の微小管に結合して、脱重合を阻害する。Ｅｔｏｐｏｓｉｄｅは、細胞のＤＮＡの複製を阻害する。化学物質Ｘは、化学物質のスクリーニングの対象である。つまり、この一例においては、未知の化学物質Ｘが、既知の化学物質であるＰａｃｌｉｔａｘｅｌとＥｔｏｐｏｓｉｄｅとのいずれと最も類似しているかが判定される。
また、この一例において、化学物質の８段階の濃度は、各化学物質の５０％阻害濃度（ＩＣ５０）が濃度の中央値になるようにして、かつ隣接する濃度間の濃度の比が約３倍になるようにして定められる。ここで５０％阻害濃度（ＩＣ５０）とは、化学物質が添加された細胞群の半数の働きが阻害される化学物質の濃度である。

また、図６の例では、培養容器１９のタイムラプス観察は、培養開始から８時間経過後を初回として、８時間おきに７２時間目まで行う。そのため、図６の例において、１つの培養容器１９の顕微鏡画像が、９枚分（８ｈ、１６ｈ、２４ｈ、３２ｈ、４０ｈ、４８ｈ、５６ｈ、６４ｈ、及び７２ｈ）を１セットとして撮像される。つまり、顕微鏡画像は、化学物質が細胞に添加された後の複数の撮像時刻において撮像される。これら撮像された顕微鏡画像は、記憶部４３に予め記憶されている。以下、撮像時刻とは、タイムラプス観察の各回において顕微鏡画像が撮像された時刻を意味する。

ステップＳ２０１：ＣＰＵ４２は、予め用意された顕微鏡画像のデータを記憶部４３から読み込む。この一例では、ＣＰＵ４２は、予め用意されている顕微鏡画像のすべてを処理対象として順次読み込む。

ステップＳ２０２：ＣＰＵ４２は、ステップＳ２０１において読み込んだ顕微鏡画像のうちから、処理対象となる画像を指定する。この一例では、ＣＰＵ４２は、ステップＳ２０１において読み込んだ顕微鏡画像のすべてを処理対象として、１枚ずつ順次指定する。

ステップＳ２０３：ＣＰＵ４２の特徴量算出部４４は、ステップＳ２０２において指定した処理対象の顕微鏡画像について、画像内に含まれる細胞の画像を抽出する。例えば、位相差顕微鏡で細胞を撮像すると、細胞膜のように位相差の変化の大きな部位の周辺にはハロが現れる。そのため、特徴量算出部４４は、細胞膜に対応するハロを公知のエッジ抽出手法で抽出するとともに、輪郭追跡処理によってエッジで囲まれた閉空間を細胞と推定する。これにより特徴量算出部４４は、顕微鏡画像から個々の細胞の画像を抽出することができる。なお、特徴量算出部４４が顕微鏡画像から個々の細胞の画像を抽出する方法は、上記のエッジ抽出手法を利用した方法以外の方法であってもよい。例えば、特徴量算出部４４は、公知のコーナー検出法を利用して顕微鏡画像から個々の細胞の画像を抽出してもよい。

ステップＳ２０４：特徴量算出部４４は、ステップＳ２０３において抽出した各細胞の画像について、複数の特徴量をそれぞれ算出する。ここで複数の特徴量とは、異なる化学物質がそれぞれ異なる濃度で添加された細胞の顕微鏡画像に含まれる細胞の形状に関する複数の種類の特徴量である。細胞の形状に関する特徴量には、各細胞の形態についての特徴量、細胞群の集団としての形態についての特徴量、及び細胞群の集団としての形態についての特徴量の時間変化を示す特徴量が含まれる。
各細胞の形態についての特徴量の種類には、細胞の面積、長さ、幅、長さと幅との比率、周囲長などの約５５種類の指標が含まれる。細胞群の集団としての形態についての特徴量には、細胞の分布の広がり、細胞の密集度などの約６０種類の指標が含まれる。ここで細胞の分布の広がりとは、一例として、細胞の位置の分散である。各細胞の位置の分散は、一例として、顕微鏡画像の所定の方向の各細胞の位置に基づいて算出される。細胞群の集団としての形態についての特徴量の時間変化を示す特徴量には、細胞の分布の広がりの時間についての変化率、細胞の密集度の時間についての変化率などの約６００種類の指標が含まれる。細胞群の集団としての形態についての特徴量の時間変化を示す特徴量は、細胞群の集団としての形態についての特徴量の時間変化を示す。ここで特徴量の時間変化とは、初回のタイムラプス観察における撮像時刻から、初回のタイムラプス観察より後の撮像時刻までの時間についての特徴量の時間変化である。

ステップＳ２０５：ＣＰＵ４２は、ステップＳ２０４において特徴量算出部４４が算出した各細胞の特徴量を、記憶部４３に書き込む。このときＣＰＵ４２は、特徴量と、この特徴量の種類、化学物質の種類、化学物質の濃度、及び撮像時刻とを関連付けて、記憶部４３に特徴量を書き込む。

ステップＳ２０６：ＣＰＵ４２は、ステップＳ２０１において記憶部４３から読み出したすべての顕微鏡画像について、特徴量の算出が終了したか否かを判定する。ＣＰＵ４２は、すべての顕微鏡画像について、特徴量の算出が終了したと判定した場合（ステップＳ２０６；ＹＥＳ）には、処理をステップＳ２０７に進める。ＣＰＵ４２は、特徴量の算出が終了していないと判定した場合（ステップＳ２０６；ＮＯ）には、処理をステップＳ２０２に戻す。

ステップＳ２０７：ＣＰＵ４２の座標算出部４５は、ステップＳ２０５において記憶部４３に書き込まれた特徴量に基づいて、細胞の形状に関する複数の特徴量それぞれを示す値の組により指定される点の特徴量空間における座標を算出する。ここで、特徴量空間とは、細胞の形状に関する複数の特徴量の種類毎に座標軸を定めた空間である。特徴量空間は、細胞の形状に関する複数の特徴量の種類毎に座標軸を定めた空間であるため、一般には高次元の空間となる。特徴量空間の次元の数は、細胞の形状に関する複数の特徴量の種類の数に等しい。

細胞の形状に関する複数の特徴量それぞれを示す値の組により指定される点とは、各細胞の形態についての特徴量それぞれの代表値と、細胞群の集団としての形態についての特徴量を示すそれぞれの値と、細胞群の集団としての形態についての特徴量の時間変化を示す特徴量示すそれぞれの値との組により指定される点である。

ここで座標算出部４５は、座標を算出する際に、各細胞の形態についての特徴量については代表値を算出し、この代表値の座標を算出する。特徴量の代表値とは、一例として、各細胞の特徴量の平均値や中央値である。また、細胞群の集団としての形態についての特徴量の一つである細胞の分布の広がりには、各細胞の位置の分散が含まれる。
以下では、特徴量空間における細胞の形状に関する複数の特徴量それぞれを示す値の組により指定される点を、特徴量点と呼ぶ。また、特徴量点からなるグラフを単にグラフと呼ぶ。

座標算出部４５は、化学物質の種類、化学物質の濃度、及び撮像時刻が同じ顕微鏡画像ごとに、特徴量点の特徴量空間における座標を算出する。
ここで、化学物質の種類には、上述したようにＰａｃｌｉｔａｘｅｌ、Ｅｔｏｐｏｓｉｄｅ、及び化学物質Ｘの３種類がある。また、化学物質の濃度には、ＩＣ５０が濃度の中央値になるようにして定められた８種類がある。また、算出対象の顕微鏡画像の撮像時刻には、８時間経過毎の時刻（８ｈ、１６ｈ、２４ｈ、３２ｈ、４０ｈ、４８ｈ、５６ｈ、６４ｈ、及び７２ｈ）の９種類がある。

座標算出部４５は、特徴量点の座標を、記憶部４３に記憶させる。ここで特徴量点の特徴量空間における座標を算出するステップＳ２０７の処理は、多変量解析の一例である。つまり、座標算出部４５は、多変量解析により解析された結果を記憶部４３に記憶させる。

ステップＳ２０８：ＣＰＵ４２の座標算出部４５は、ステップＳ２０７において特徴量空間における座標を算出した特徴量点について、主成分分析を行う。
座標算出部４５は、公知の主成分分析により、各細胞から得られた高次元の特徴量を次元圧縮し、主成分の座標とする。
特徴量点は、視覚化のために２次元や３次元の空間などの低次元の空間に射影されてよい。ここで視覚化とは、例えば、特徴量点をグラフにして表示することである。本実施形態では、一例として、座標算出部４５が特徴量点を、主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面に射影する場合について説明する。以下、主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面を、特徴量平面と呼ぶ。また、視覚化のために特徴量平面に射影された特徴量点を射影特徴量点と呼ぶ。

図７は、主成分分析により処理された射影特徴量点の一例を示す図である。つまり、図７は、特徴量平面における射影特徴量点の一例を示す図である。この図７に示すグラフにおいて、軸ＰＣ１とは、特徴量空間において特徴量点のばらつきが最も大きくなる方向と平行な軸である。また、この図７に示すグラフにおいて、軸ＰＣ２とは、軸ＰＣ１に直交する方向のうち特徴量空間において特徴量点のばらつきが最も大きくなる方向と平行な軸である。図７には、ある細胞に３種類の化学物質、すなわち化学物質Ｃｈ１、化学物質Ｃｈ２、化学物質Ｃｈ３を添加した場合のグラフを示す。この一例では、化学物質Ｃｈ１とは、Ｐａｃｌｉｔａｘｅｌであり、化学物質Ｃｈ２とは、Ｅｔｏｐｏｓｉｄｅであり、化学物質Ｃｈ３とは、化学物質Ｘである。また、この図７に示すグラフにおいて、プロットＰ１−ｍは、化学物質Ｃｈ１が添加された細胞の射影特徴量点を示す。また、このプロットＰ１−ｍの添字ｍは、細胞に添加した化学物質の濃度を、低い方から高い方に順に示す。すなわち、プロットＰ１−１は、最も低い濃度の化学物質Ｃｈ１が添加された細胞の射影特徴量点を示す。また、プロットＰ１−８は、最も高い濃度の化学物質Ｃｈ１が添加された細胞の射影特徴量点を示す。

本実施形態における特徴量には細胞群の集団としての形態についての特徴量の時間変化を示す特徴量が含まれているため、座標算出部４５は、特徴量の時間変化を含めて特徴量の中から主成分の軸（軸ＰＣ１及び軸ＰＣ２）を決定している。これにより、座標算出部４５は、時間変化を含めた特徴量の情報を損なうことなく、多くの種類の特徴量を２軸の指標に要約することができる。すなわち、座標算出部４５によれば、化学物質を添加した後の細胞の形態変化を、その形態変化の特徴を際立たせつつ、視覚化することができる。

図６に戻り、濃度軌跡算出部４６による濃度軌跡の算出について説明する。
ステップＳ２０９：濃度軌跡算出部４６は、座標算出部４５が特徴量空間における座標を算出した特徴量点から、特徴量点の濃度軌跡を算出する。ここで濃度軌跡とは、特徴量点を、細胞に添加された化学物質の濃度の順に並べることにより得られる軌跡である。つまり、濃度軌跡算出部４６は、特徴量算出部４４が算出した特徴量点を、細胞に添加された物質の濃度の順に並べてできる軌跡を算出する。

ここでは、濃度軌跡算出部４６が、特徴量点が主成分分析により２次元平面に射影された特徴量点から濃度軌跡を算出する場合について説明する。濃度軌跡の具体例について図８を参照し説明する。図８は、本実施形態の濃度軌跡の一例を示す図である。
図８に示すように、濃度軌跡算出部４６は、最も低い濃度の化学物質Ｃｈ１が添加された細胞のプロットＰ１−１から、最も高い濃度の化学物質Ｃｈ１が添加された細胞のプロットＰ１−８までを、順に並べることにより、特徴量点の軌跡を算出する。より具体的には、濃度軌跡算出部４６は、プロットＰ１−１と、プロットＰ１−２とを線分Ｌ１によって接続する。また、濃度軌跡算出部４６は、プロットＰ１−２と、プロットＰ１−３とを線分Ｌ２によって接続する。さらに、濃度軌跡算出部４６は、プロットＰ１−３からプロットＰ１−８までを線分Ｌ３から線分Ｌ７によって順に接続する。

ステップＳ２１０：表示部５０は、ステップＳ２０８及びステップＳ２０９において作成された軌跡をグラフとして表示する。この一例においては、表示部５０は、図８に示すグラフを表示する。

なお、特徴量点を２次元や３次元の空間の低次元の空間に射影する処理や、軌跡をグラフとして表示する処理は視覚化のために行われる処理であり省略されてよい。つまり、ステップＳ２０８、ステップＳ２０９、及びステップＳ２１０の処理は省略されてもよい。 以降の処理においては、特徴量空間における特徴量点について処理が行われる。ただし、以降の処理は、特徴量平面などの低次元の空間に射影された射影特徴量点について行われてもよい。例えば、以降の処理は、図７において示した特徴量平面における射影特徴量点について行われてもよい。低次元の空間に射影された射影特徴量点について行われる場合であっても、低次元の空間に射影された射影特徴量点についての処理は、特徴量空間における特徴量点についての処理と同様である。

≪変化領域の特定の例≫
ステップＳ２１１：変化領域判定部４７は、形態変化領域を、細胞に添加された化学物質の濃度の変化に対する特徴量空間における特徴量点の変化の大きさに基づいて判定する。ここで形態変化領域とは、特徴量空間において、細胞に添加された化学物質の濃度の変化に対する特徴量点の変化の大きさが所定の値以上となる化学物質の濃度の範囲である。つまり、変化領域判定部４７は、形態変化領域を、細胞に添加された化学物質の濃度の変化に対する特徴量の変化に基づいて判定する。

ここで、特徴量点とは、特徴量空間における細胞の形状に関する複数の特徴量それぞれの代表値の組により指定される点であるから、特徴量点の変化は細胞の形態変化を示す。 また、細胞に添加された化学物質の濃度の変化に対する特徴量点の変化は、この化学物質の種類に応じて異なる。例えば、強い毒性をもつ化学物質が細胞に添加された場合、添加される化学物質の濃度を僅かに上げただけで細胞は死滅してしまい、その死滅への過程において化学物質の濃度の変化に対し特徴量は急激に変化する。特徴量が急激に変化した後（つまり、細胞が死滅した後）は、添加される化学物質の濃度を上げても細胞は死滅したままであり特徴量に変化はみられない。一方、細胞に作用する効果の小さい（例えば、毒性が弱い）化学物質が細胞に添加された場合、添加される化学物質の濃度をある程度まで上げないと特徴量に変化はみられない。

本実施形態では、化学物質の濃度の範囲は、コントロール領域、形態変化領域、及び死細胞領域に分類される。
コントロール領域とは、細胞に添加される化学物質の濃度を最も低い濃度から上げたときに、化学物質の濃度に対する特徴量の変化が小さい濃度の範囲である。コントロール領域では、化学物質の細胞に作用する効果は小さい。細胞に作用する効果の小さい化学物質では、コントロール領域の範囲が大きくなる。
形態変化領域とは、細胞に添加される化学物質の濃度を上げていったときに、化学物質の濃度に対する特徴量の変化が大きくなる濃度の範囲である。形態変化領域では、化学物質の細胞に作用する効果は大きい。以下、形態変化領域に含まれる濃度を形態変化濃度ともいう。形態変化濃度は、細胞の形態が変化する濃度に対応する。
死細胞領域とは、細胞に添加される化学物質の濃度を上げていったときに、細胞が死滅し化学物質の濃度に対する特徴量の変化が小さくなる濃度の範囲である。死細胞領域では細胞が死滅しているため、特徴量は小さくなり、そのため化学物質の濃度に対する特徴量の変化も小さくなる。また、添加する細胞に対して強い毒性をもつ化学物質では、死細胞領域の範囲が大きくなる。

なお、本実施形態では、説明の便宜上、化学物質の濃度を、所定の濃度を単位として表す。例えば、濃度が１とは、所定の濃度の１倍の濃度であり、濃度が８とは、所定の濃度の８倍の濃度である。本実施形態では、細胞に添加される各化学物質の濃度の中で最も低い濃度は１であり、細胞に添加される各化学物質の濃度の中で最も高い濃度は８である。

変化領域判定部４７は、細胞に添加される化学物質の濃度の範囲の中からコントロール領域及び細胞死領域を判定し、判定したコントロール領域及び細胞死領域を除いた化学物質の濃度の範囲を形態変化領域であると判定する。ここでは一例として、変化領域判定部４７は、形態変化領域を、細胞に添加された物質の濃度の変化に対する特徴量の変化に基づいて判定する。図９から図１２を参照し、変化領域判定部４７の形態変化領域の判定の方法について説明する。

図９は、本実施形態に係る特徴量点のコントロール領域中心からの距離と化学物質濃度との関係の一例を示す図である。ここで特徴量点のコントロール領域中心からの距離とは、特徴量空間における距離である。この距離は、特徴量点の座標と、コントロール領域中心の座標とに基づいて変化領域判定部４７により算出される。この距離の単位は無次元である。グラフＣｄ１は、Ｅｔｏｐｏｓｉｄｅについて特徴量点のコントロール領域の中心からの距離と化学物質濃度との関係を示すグラフである。グラフＣｄ２は、Ｐａｃｌｉｔａｘｅｌについて特徴量点のコントロール領域の中心からの距離と化学物質濃度との関係を示すグラフである。グラフＣｄ３は、化学物質Ｘについて特徴量点のコントロール領域の中心からの距離と化学物質濃度との関係を示すグラフである。

コントロール領域を判定するために、変化領域判定部４７は、まずコントロール領域の中心の座標を判定する。ここでコントロール領域の中心とは、一例として、細胞に添加される各化学物質の濃度の中で最も低い濃度に対する特徴量点間の重心である。例えば、コントロール領域の中心は、濃度が１のＥｔｏｐｏｓｉｄｅに対する特徴量点と、濃度が１のＰａｃｌｉｔａｘｅｌに対する特徴量点と、濃度が１の化学物質Ｘに対する特徴量点との重心である。

変化領域判定部４７は、特徴量空間において、各化学物質に対する特徴量点のコントロール領域の中心からの距離が、コントロール距離より小さい特徴量点を、特徴量点の座標と、コントロール領域の中心の座標に基づいて判定する。ここで、コントロール距離は、各化学物質に共通して予め決められた距離である。例えば、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅに対する特徴量点のコントロール領域の中心からの距離が、コントロール距離より小さい特徴量点を判定する。例えば、変化領域判定部４７は、Ｐａｃｌｉｔａｘｅｌに対する特徴量点のコントロール領域の中心からの距離が、コントロール距離より小さい特徴量点を判定する。例えば、変化領域判定部４７は、化学物質Ｘに対する特徴量点のコントロール領域の中心からの距離が、コントロール距離より小さい特徴量点を判定する。コントロール距離を示す情報は、記憶部４３に記憶される。
なお、コントロール距離は、化学物質ごとに予め決められてもよい。

なお、コントロール距離は、細胞に添加される各化学物質の最も低い濃度（１の濃度）に対する特徴量点の特徴量空間における座標に基づいて決められる距離が用いられてもよい。例えば、コントロール距離は、細胞に添加される各化学物質の最も低い濃度（１の濃度）に対する特徴量点とコントロール中心との距離それぞれが、コントロール距離以下となるように決定される。また、コントロール距離は、細胞の形態変化に則して決められる距離が用いられてもよい。コントロール距離として、細胞の形態変化に則して決められる距離が用いられる場合、コントロール距離を大きくすればするほど、コントロール領域に含まれる濃度に対する細胞の形態変化は大きくなる。
本実施形態において距離とは、一例として、ユークリッド距離やマハラビノス距離である。

変化領域判定部４７は、判定した特徴量点に対応する化学物質の濃度をコントロール領域に含まれる濃度であると判定する。変化領域判定部４７は、判定した濃度に基づいてコントロール領域を判定する。

変化領域判定部４７は、グラフＣｄ１から、Ｅｔｏｐｏｓｉｄｅについて化学物質濃度が１及び２に対する特徴量点がコントロール領域の中心からの距離が、コントロール距離より小さいと特徴量点の座標と、コントロール領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅのコントロール領域が濃度１及び２であると判定する。変化領域判定部４７は、グラフＣｄ２から、Ｐａｃｌｉｔａｘｅｌについて化学物質濃度が１、２、３、４及び５に対する特徴量点がコントロール領域の中心からの距離が、コントロール半径より小さいと特徴量点の座標と、コントロール領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、Ｐａｃｌｉｔａｘｅｌのコントロール領域が濃度１、２、３、４及び５であると判定する。変化領域判定部４７は、グラフＣｄ３から、化学物質Ｘについて化学物質濃度が１、２、３、４、５、６及び７に対する特徴量点のコントロール領域の中心からの距離が、コントロール半径より小さいと特徴量点の座標と、コントロール領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、化学物質Ｘのコントロール領域が濃度１、２、３、４、５、６及び７であると判定する。

図１０は、本実施形態に係る化学物質ごとの形態変化領域の開始濃度の一例を示す図である。形態変化領域の開始濃度とは、細胞の形状に関する複数の種類の特徴量が所定の量だけ変化し始めた化学物質の濃度である。ここで所定の量とは、コントロール距離に対応する化学物質の濃度の範囲において、細胞の形状に関する複数の種類の特徴量が変化する量である。例えば、形態変化領域の開始濃度とは、形態変化領域に含まれる濃度のうち最も低い濃度である。変化領域判定部４７は、コントロール領域の中において最も高い化学物質濃度の次に高い濃度を形態変化領域の開始濃度であると判定する。例えば、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅの形態変化領域開始濃度を３であると判定する。変化領域判定部４７は、Ｃｙｔｏｃｈａｌａｓｉｎ Ｂの形態変化領域の開始濃度を７であると判定する。変化領域判定部４７は、化学物質Ｘの形態変化領域の開始濃度を８であると判定する。

このように、変化領域判定部４７は、コントロール領域に基づいて形態変化領域の開始濃度を求める。ここでコントロール領域は、図９において説明したように、細胞に添加される化学物質の濃度の変化に対する特徴量の変化に基づいて求められる。したがって、変化領域判定部４７は、細胞に添加される化学物質の濃度の変化に対する特徴量の変化に基づいて形態変化濃度を求める。さらに、図９において説明したように、形態変化濃度を求めることは、特徴量の変化量と閾値とを比較して形態変化濃度を求める。

Ｅｔｏｐｏｓｉｄｅ、Ｄｏｘｏｒｕｂｉｃｉｎ、Ｌａｔｒｕｎｃｕｌｉｎ Ａ、Ｖｉｎｂｌａｓｔｉｎｅ、及び１７−ＡＧＧは、他の化学物質に比べて形態変化領域の開始濃度が低いことがわかる。５−Ｆｌｕｏｒｏｕｒａｃｉｌ、Ｃｙｔｏｃｈａｌａｓｉｎ Ｂ、及びＰａｃｌｉｔａｘｅｌは、Ｅｔｏｐｏｓｉｄｅ、Ｄｏｘｏｒｕｂｉｃｉｎ、Ｌａｔｒｕｎｃｕｌｉｎ Ａ、Ｖｉｎｂｌａｓｔｉｎｅ、及び１７−ＡＧＧに比べて形態変化領域の開始濃度が高いことがわかる。化学物質Ｘは、他の化学物質に比べて形態変化領域の開始濃度が高いことがわかる。

図１１は、本実施形態に係る特徴量点の死細胞領域の中心からの距離と化学物質濃度との関係の一例を示す図である。ここで特徴量点の死細胞領域の中心からの距離とは、特徴量空間における距離である。この距離は、特徴量点の座標と、死細胞領域の中心の座標とに基づいて変化領域判定部４７により算出される。グラフＤｄ１は、Ｅｔｏｐｏｓｉｄｅについて特徴量点の死細胞領域の中心からの距離と化学物質濃度との関係を示すグラフである。グラフＤｄ２は、Ｐａｃｌｉｔａｘｅｌについて特徴量点の死細胞領域の中心からの距離と化学物質濃度との関係を示すグラフである。グラフＤｄ３は、化学物質Ｘについて特徴量点の死細胞領域の中心からの距離と化学物質濃度との関係を示すグラフである。

変化領域判定部４７は、特徴量空間において、各化学物質に対する特徴量点の死細胞領域の中心からの距離が、死細胞距離より小さい特徴量点を死細胞領域に含まれる特徴量点であると判定する。ここで死細胞領域の中心とは、例えば、細胞に添加される各化学物質の濃度の中で最も高い濃度に対する特徴量点間の重心である。死細胞領域の中心は、例えば、濃度が８のＥｔｏｐｏｓｉｄｅに対する特徴量点と、濃度が８のＰａｃｌｉｔａｘｅｌに対する特徴量点と、濃度が８の化学物質Ｘに対する特徴量点との重心である。死細胞距離とは、化学物質に共通して予め決められた距離である。例えば、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅに対する特徴量点の死細胞領域の中心からの距離が、死細胞距離より小さい特徴量点を判定する。例えば、変化領域判定部４７は、Ｐａｃｌｉｔａｘｅｌに対する特徴量点の死細胞領域の中心からの距離が、死細胞距離より小さい特徴量点を判定する。例えば、変化領域判定部４７は、化学物質Ｘに対する特徴量点の死細胞領域の中心からの距離が、死細胞距離より小さい特徴量点を判定する。死細胞距離を示す情報は、記憶部４３に記憶される。
なお、死細胞距離は、化学物質ごとに予め決められてもよい。

なお、死細胞距離は、細胞に添加される各化学物質の最も高い濃度（８の濃度）に対する特徴量点の特徴量空間における座標に基づいて決められる距離が用いられてもよい。また、死細胞距離は、細胞の形態変化に則して決められる距離が用いられてもよい。死細胞距離として、細胞の形態変化に則して決められる距離が用いられる場合、死細胞距離を大きくすればするほど、死細胞領域に含まれる濃度に対する細胞の形態変化は大きくなる。

変化領域判定部４７は、グラフＤｄ１から、Ｅｔｏｐｏｓｉｄｅについて化学物質濃度が７及び８に対する特徴量点が死細胞領域の中心からの距離が、死細胞距離より小さいと特徴量プロットの座標と、死細胞領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅの死細胞領域が濃度７及び８であると判定する。変化領域判定部４７は、グラフＤｄ２から、Ｐａｃｌｉｔａｘｅｌについて化学物質濃度が６、７及び８に対する特徴量点が死細胞領域の中心からの距離が、死細胞距離より小さいと特徴量点の座標と、死細胞領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、Ｐａｃｌｉｔａｘｅｌの死細胞領域が濃度６、７及び８であると判定する。変化領域判定部４７は、グラフＤｄ３から、化学物質Ｘについて特徴量点が死細胞領域の中心からの距離が、死細胞距離より小さい化学物質濃度はないと特徴量点の座標と、死細胞領域の中心の座標に基づいて判定する。したがって、変化領域判定部４７は、化学物質Ｘの死細胞領域はないと判定する。

図１２は、本実施形態に係る化学物質ごとの死細胞領域の一例を示す図である。変化領域判定部４７は、各化学物質について、形態変化領域の開始濃度より高い濃度の範囲から死細胞領域を除いた濃度の範囲を形態変化領域と判定する。つまり、変化領域判定部４７は、形態変化領域の開始濃度より高い濃度の範囲から、細胞の形状に関する複数の種類の特徴量の変化が所定の量以下となる濃度の範囲を除いた範囲を形態変化領域と判定する。例えば、変化領域判定部４７は、Ｅｔｏｐｏｓｉｄｅについて、形態変化領域の開始濃度３より高い濃度の範囲から死細胞領域である濃度７及び８の範囲を除いた領域を形態変化領域と判定する。例えば、変化領域判定部４７は、Ｐａｃｌｉｔａｘｅｌについて、形態変化領域の開始濃度である濃度６より高い濃度の範囲から死細胞領域である濃度６、７及び８の濃度の範囲を除くと濃度の範囲は残らないため、形態変化領域は存在しないと判定する。変化領域判定部４７は、化学物質Ｘについて、死細胞領域は存在しないため、形態変化領域の開始濃度より高い濃度の範囲である８の濃度の範囲を形態変化領域と判定する。

このように、形態変化濃度は、細胞の形態が変化してから細胞が細胞死するまでの濃度を含む。
また、このように、変化領域判定部４７は、死細胞領域に基づいて形態変化領域を求める。ここで死細胞領域は、図１１において説明したように、細胞に添加される化学物質の濃度の変化に対する特徴量の変化に基づいて求められる。したがって、変化領域判定部４７は、細胞に添加される化学物質の濃度の変化に対する特徴量の変化に基づいて形態変化濃度を求める。さらに、図９において説明したように、形態変化濃度を求めることは、特徴量の変化量と閾値とを比較して形態変化濃度を求める。

再び、図６に戻り化学物質のスクリーニングの動作の説明を続ける。
ステップＳ２１２：類似性判定部４８は、ステップＳ２１１において変化領域判定部４７により判定された形態変化領域を、細胞に添加された化学物質の間において比較することにより、化学物質同士の類似性を判定する。つまり、類似性判定部４８は、異なる特徴量の内、細胞の形態が変化する形態変化濃度に対応する特徴量を異なる化学物質同士で比較して、化学物質同士の類似性を判定する。したがって、類似性判定部４８は、ステップＳ２０７及びステップＳ２１１の多変量解析により解析された結果が示す濃度の変化に対する特徴量の変化に基づいて化学物質同士の類似性を判定する。つまり、類似性判定部４８は、複数の特徴量が多変量解析により解析された結果に基づいて化学物質同士の類似性を判定する。

ここで化学物質同士の類似性とは、化学物質が細胞に添加されたときに、化学物質が変化させる細胞の形状の特徴量の種類や変化の大きさが、化学物質同士において類似していることである。本実施形態において化学物質同士の類似性を判定することで、化学物質同士の作用機序などの性質の類似性を判定できる。すなわち、既知の化学物質に対する未知の化学物質の作用機序などの性質の類似性が判定できるため、効率的に化学物質のスクリーニングを行うことができる。

ここで図１３を参照し図１２におけるステップＳ２１２において類似性判定部４８が化学物質同士の類似性を判定する処理の詳細について説明する。
図１３は、本実施形態に係る制御装置による化学物質同士の特徴量空間における特徴量点間の比較の動作の一例を示すフローチャートである。

ステップＳ３０１：類似性判定部４８は、類似性を判定する対象の化学物質の特徴量空間における特徴量点について、一方の形態変化領域に含まれる特徴量点のそれぞれについて、他方の形態変化領域に含まれる特徴量点のうち最も近い距離にある特徴量点を特徴量点の座標に基づいて探索する。

ここでステップＳ３０１の類似性判定部４８の動作について図１４を参照し説明する。 図１４は、本実施形態に係る制御装置による化学物質同士の特徴量空間における特徴量点間の比較の方法の一例を示す図である。類似性判定部４８の処理は高次元の特徴量空間における特徴量点について行われる。ただし、図１４では、視覚化のために主成分分析の結果得られた特徴量平面に、化学物質Ａ、化学物質Ｂ、及び化学物質Ｘについての射影特徴量点が示されている。つまり、図１４に示す特徴量平面は、特徴量空間のある平面が抜き出された平面であり、図１４に示す射影特徴量点は、特徴量空間における特徴量点がこの平面に射影された点である。

図１４に示す例では、類似性判定部４８は、化学物質Ｘについて、化学物質Ａとの類似性、及び化学物質Ｂとの類似性を判定する。
化学物質Ａについての形態変化領域ＣｈＡは、プロットＰＡ−１〜プロットＰＡ−５からなる。化学物質Ｂについての形態変化領域ＣｈＢは、プロットＰＢ−１〜プロットＰＢ−５からなる。化学物質Ｘについての形態変化領域ＣｈＸは、プロットＰＸ−１〜プロットＰＸ−３からなる。

類似性判定部４８は、形態変化領域ＣｈＸに含まれるプロットＰＸ−１からの距離が最も近い特徴量点を、形態変化領域ＣｈＡに含まれるプロットＰＡ−１〜プロットＰＡ−５の中から探索する。類似性判定部４８は、プロットＰＸ−１からの距離が距離ｄ１であるプロットＰＡ−３を、プロットＰＸ−１からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。同様に、類似性判定部４８は、プロットＰＸ−２からの距離が距離ｄ２であるプロットＰＡ−４を、プロットＰＸ−２からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。類似性判定部４８は、プロットＰＸ−３からの距離が距離ｄ３であるプロットＰＡ−５を、プロットＰＸ−３からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。
ただし、形態変化領域に含まれる特徴量点に対応する化学物質の濃度と、この特徴量点からの距離が最も近いと判定される特徴量点に対応する化学物質の濃度とは同じであるとは限らない。例えば、プロットＰＸ−１に対応する化学物質の濃度と、プロットＰＡ−３に対応する化学物質の濃度とは同じであるとは限らない。

類似性判定部４８は、形態変化領域ＣｈＸに含まれるプロットＰＸ−１からの距離が最も近い特徴量点を、形態変化領域ＣｈＢに含まれるプロットＰＢ−１〜プロットＰＢ−５の中から探索する。類似性判定部４８は、プロットＰＸ−１からの距離が距離ｄ４であるプロットＰＢ−２を、プロットＰＸ−１からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。同様に、類似性判定部４８は、プロットＰＸ−２からの距離が距離ｄ５であるプロットＰＢ−３を、プロットＰＸ−２からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。類似性判定部４８は、プロットＰＸ−３からの距離が距離ｄ６であるプロットＰＢ−３を、プロットＰＸ−３からの距離が最も近い特徴量点であると特徴量点の座標に基づいて判定する。

図１３に戻り、制御装置による化学物質同士の特徴量平面における特徴量点間の比較の動作の説明を続ける。
ステップＳ３０２：類似性判定部４８は、化学物質Ｘの形態変化領域の特徴量点と、ステップＳ３０１において判定した化学物質Ａの形態変化領域の最も近い特徴量点とのそれぞれの距離の平均を特徴量点の座標に基づいて算出する。類似性判定部４８は、例えば、形態変化領域ＣｈＸと形態変化領域ＣｈＡとに対して、距離ｄ１、距離ｄ２、及び距離ｄ３の３つの距離の平均を特徴量点の座標に基づいて算出する。
一方、類似性判定部４８は、化学物質Ｘの形態変化領域の特徴量点と、ステップＳ３０１において判定した化学物質Ｂの形態変化領域の最も近い特徴量点とのそれぞれの距離の平均を特徴量点の座標に基づいて算出する。類似性判定部４８は、例えば、形態変化領域ＣｈＸと形態変化領域ＣｈＢとに対して、距離ｄ４、距離ｄ５、及び距離ｄ６の３つの距離の平均を特徴量点の座標に基づいて算出する。

ステップＳ３０３：類似性判定部４８は、類似性の判定対象である化学物質のなかから、未知の化学物質Ｘに最も類似する既知の化学物質を判定する。類似性判定部４８は、ステップＳ３０２において算出した形態変化領域の特徴量点間の距離の平均が最も小さい化学物質を、未知の化学物質Ｘに最も類似する既知の化学物質であると判定する。
図１４に示す例においては、類似性判定部４８は、形態変化領域ＣｈＡをもつ化学物質Ａを、化学物質Ｘに最も類似する既知の化学物質であると判定する。

制御装置４１は、判定結果に基づいて、細胞に生じさせる形態変化について、未知の化学物質Ｘを分類することができる。
なお、本実施形態では、制御装置４１は、未知の化学物質Ｘと既知の化学物質との類似性を判定する場合について説明したが、これに限らない。制御装置４１は、既知の化学物質間において類似性を判定してもよい。

≪まとめ≫
以上説明したように、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法は、異なる化学物質がそれぞれ異なる濃度で添加された細胞の顕微鏡画像に含まれる細胞の形状に関する複数の特徴量に基づいて、異なる化学物質同士の類似性を判定することと、判定した結果を出力することとを有する。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、細胞の形状に関する複数の特徴量に基づいて異なる化学物質同士の類似性を判定できるため、化学物質のスクリーニングにおけるコストと時間を削減することができる。 本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、化学物質同士の類似性を判定することで、化学物質同士の作用機序などの性質の類似性を判定できる。すなわち、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、既知の化学物質に対する未知の化学物質の作用機序などの性質の類似性が判定できるため、効率的に化学物質のスクリーニングを行うことができる。

また、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、複数の特徴量が多変量解析により解析された結果に基づいて化学物質同士の類似性を判定する。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、複数の特徴量が多変量解析により解析された結果に基づいて化学物質同士の類似性を判定することができるため、多変量解析により解析された結果に基づいて化学物質同士の類似性を判定しない場合に比べて化学物質のスクリーニングの精度を高くすることができる。

また、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、異なる化学物質同士の類似性を判定することは、異なる特徴量の内、細胞の形態が変化する形態変化濃度に対応する特徴量を異なる化学物質同士で比較して、異なる化学物質同士の類似性を判定する。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、化学物質が細胞に生じさせる形態変化が大きい化学物質の濃度の範囲に限定して化学物質同士の類似性を判定できるため、化学物質が細胞に生じさせる形態変化が大きい化学物質の濃度の範囲に限定せずに化学物質同士の類似性を判定する場合に比べて、化学物質のスクリーニングにおけるコストと時間を削減することができる。

また、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、細胞に添加される化学物質の濃度の変化に対する特徴量の変化に基づいて形態変化濃度を求めることを有する。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、細胞の形態の変化が開始される濃度より高い濃度の範囲に限定して化学物質同士の類似性を判定できるため、細胞の形態の変化が開始される濃度より高い濃度の範囲に限定しない場合に比べ、化学物質のスクリーニングにおけるコストと時間を削減することができる。

また、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、形態変化濃度は、細胞の形態が変化してから細胞が細胞死するまでの濃度を含む。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、細胞の細胞死を生じさせる濃度より低い濃度の範囲に限定して化学物質同士の類似性を判定できるため、細胞の細胞死を生じさせる濃度より低い濃度の範囲に限定しない場合に比べ、化学物質のスクリーニングにおけるコストと時間を削減することができる。

本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、形態変化濃度を求めることは、特徴量の変化量と閾値とを比較して形態変化濃度を求める。
この構成により、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、形態変化濃度は特徴量の変化量と閾値とを比較して求められるため、形態変化濃度を簡便に求めることができる。

また、本実施形態の制御装置４１が提供する化学物質のスクリーニング方法では、異なる化学物質は、既知の化学物質と、未知の化学物質とが含まれ、既知の化学物質と、未知の化学物質との類似性を判定する。
この構成により、既知の化学物質と、未知の化学物質との類似性を判定できるため、未知の化学物質が既知の化学物質のいずれと類似しているかを判定することができる。

なお、制御装置４１は、ステップＳ２０７及びステップＳ２１１での複数の特徴量を多変量解析により解析する処理を行う代わりに、ステップＳ２０７及びステップＳ２１１での複数の特徴量を多変量解析により解析する処理を外部装置に行わせてもよい。ここで外部装置とは、制御装置４１の外部に設けられた装置である。この外部装置は、制御装置４１とは別に下部ケーシング１３に設けられてもよいし、培養観察装置１１とは独立して設けられてもよい。この外部装置は、多変量解析の演算処理を実行するプロセッサ、及び多変量解析の解析結果を記憶する記憶部を備える。

制御装置４１がステップＳ２０７及びステップＳ２１１での複数の特徴量を多変量解析により解析する処理を外部装置に行わせる場合、ステップＳ２０５においてＣＰＵ４２は、ステップＳ２０４において特徴量算出部４４が算出した各細胞の特徴量を、外部装置に供給する。ステップＳ２１２では、ＣＰＵ４２は、外部装置が行った複数の特徴量を多変量解析により解析する処理の結果を取得し、類似性判定部４８に供給する。類似性判定部４８は、外部装置が行った複数の特徴量を多変量解析により解析する処理の結果に基づいて、異なる化学物質同士の類似性を判定する。

また、制御装置４１は、ステップＳ２０１からステップＳ２０５までの、顕微鏡画像に含まれる細胞の形状に関する複数の特徴量を算出し記憶させる処理を、外部装置に行わせてもよい。この外部装置は、上述の外部装置と同じ装置であってもよいし、上述の外部装置とは別の装置であってもよい。制御装置４１がステップＳ２０１からステップＳ２０５までの処理を外部装置に行わせる場合、制御装置４１は、外部装置が算出した複数の特徴量を取得し、ステップＳ２０６以降の処理を行う。

表示部５０は、類似性判定部４８が判定した異なる化学物質同士の類似性の判定結果を表示する。つまり、本実施形態の化学物質のスクリーニング方法は、判定した結果を出力することとを有する。

なお、本実施形態においては、ステップＳ２１１において変化領域判定部４７が形態変化領域を判定する際に、コントロール領域及び細胞死領域を判定する場合について説明したが、これに限らない。コントロール領域を示す情報、及び細胞死領域を示す情報が、予め記憶部４３に記憶されていてもよい。コントロール領域を示す情報、及び細胞死領域を示す情報が予め記憶部４３に記憶されている場合、変化領域判定部４７は、記憶部４３からコントロール領域を示す情報、及び細胞死領域を示す情報を取得し、取得したコントロール領域を示す情報が示すコントロール領域、及び取得した細胞死領域を示す情報が示す細胞死領域を除いた化学物質の濃度の範囲を形態変化領域であると判定してよい。

また、本実施形態においては、ステップＳ２１１において変化領域判定部４７は、形態変化領域を判定するために、コントロール領域、細胞死領域、及び形態変化領域の開始濃度を判定したが、これに限らない。ステップＳ２１１において変化領域判定部４７は、形態変化領域を判定する代わりに、コントロール領域、細胞死領域、及び形態変化領域の開始濃度のうち少なくも１つを判定してもよい。つまり、変化領域判定部４７は、多変量解析により解析された結果に基づいて、細胞の形態の変化が開始される濃度を判定してもよい。また、変化領域判定部４７は、多変量解析により解析された結果に基づいて、前細胞の細胞死を生じさせる濃度を判定してもよい。

また、変化領域判定部４７は、コントロール領域及び細胞死領域を求めずに、形態変化領域を求めてもよい。例えば、変化領域判定部４７は、距離の上限と下限に含まれる濃度を形態変化領域と判定してもよい。変化領域判定部４７は、多変量解析により解析された結果に基づいて、細胞の形態の変化が開始される濃度、及び細胞の形態の変化が終了する濃度を判定してもよい。変化領域判定部４７は、判定した細胞の形態の変化が開始される濃度から細胞の形態の変化が終了する濃度の範囲を、形態変化領域と判定してよい。

ステップＳ２１１において変化領域判定部４７が、コントロール領域、細胞死領域、及び形態変化領域の開始濃度のうち少なくも１つを判定する場合、ステップＳ２１２において類似性判定部４８は、コントロール領域、細胞死領域、及び形態変化領域の開始濃度のうち少なくも１つに基づいて、化学物質の濃度の範囲を限定し、化学物質同士の類似性を判定してもよい。

例えば、類似性判定部４８は、コントロール領域を除いた化学物質の濃度の範囲を、細胞に添加された化学物質の間において比較することにより、化学物質同士の類似性を判定してよい。別の一例では、類似性判定部４８は、細胞死領域を除いた化学物質の濃度の範囲を、細胞に添加された化学物質の間において比較することにより、化学物質同士の類似性を判定してよい。また、別の一例では、類似性判定部４８は、形態変化領域の開始濃度より高い化学物質の濃度の範囲を、細胞に添加された化学物質の間において比較することにより、化学物質同士の類似性を判定してよい。また、類似性判定部４８は、コントロール領域において化学物質同士の類似性を判定してもよい。類似性判定部４８は、細胞死領域において化学物質同士の類似性を判定してもよい。

なお、上記の実施形態においては、コントロール領域は、細胞に添加される各化学物質の最も低い濃度に対する特徴量点間の重心からの所定の距離に基づいて判定され、死細胞領域は、細胞に添加される各化学物質の最も高い濃度に対する特徴量点間の重心からの所定の距離に基づいて判定される場合について説明したが、コントロール領域や死細胞領域は、予め取得された既知の化学物質のデータを用いて決定されてもよい。

なお、上記の実施形態においては、既知の化学物質がＰａｃｌｉｔａｘｅｌ及びＥｔｏｐｏｓｉｄｅである場合について説明したが、これに限らない。既知の化学物質は、存在するいずれの化学物質であってもよい。例えば、既知の化学物質は、ＰＤ−１８０９７０や、図１０及び図１２において示したＤｏｘｏｒｕｂｉｃｉｎ、５−Ｆｌｕｏｒｏｕｒａｃｉｌ、Ｃｙｔｏｃｈａｌａｓｉｎ Ｂ、Ｌａｔｒｕｎｃｕｌｉｎ Ａ、Ｖｉｎｂｌａｓｔｉｎｅ、及び１７−ＡＡＧであってもよい。

また、上記の実施形態においては、既知の化学物質が２種類である場合について説明したが、これに限らない。上記の実施形態における化学物質のスクリーニングにおいて、１種類、または３種類以上の既知の化学物質が用いられてもよい。

また、上記の実施形態においては、コントロール領域を判定する際に、コントロール領域の中心が、細胞に添加される各化学物質の濃度の中で最も低い濃度に対する特徴量点間の重心である場合について説明したが、これに限らない。コントロール領域の中心は、細胞の形態の変化に基づいて決定されてもよい。つまり、コントロール領域の中心は、特徴量の変化に基づいて決定されてもよい。

また、上記の実施形態においては、死細胞領域を判定する際に、死細胞領域の中心が、一例として、細胞に添加される各化学物質の濃度の中で最も高い濃度に対する特徴量点間の重心である場合について説明したが、これに限らない。死細胞領域の中心は、例えば、生細胞に対する死細胞の割合に基づいて決定されてもよい。ここで死細胞領域の中心が生細胞に対する死細胞の割合に基づいて決定される場合、死細胞領域の中心は、生細胞に対する死細胞の割合が所定の割合以上となる化学物質の濃度に対する特徴量点に基づいて決定される。

また、上記の実施形態においては、一例として、形態変化領域において異なる化学物質同士の類似性を判定する場合について説明したが、これに限らない。例えば、添加した化学物質の全濃度範囲において類似性を判定してもよい。また、添加した化学物質の一部の濃度範囲において類似性を判定してもよい。

（第２実施形態）
以下、図面を参照しながら本発明の第２の実施形態について詳しく説明する。
上記第１の実施形態では、制御装置が、細胞に生じさせる形態変化が大きい化学物質の濃度を判定してから化学物質同士の類似性について判定する場合について説明した。本実施形態では、制御装置は、５０％阻害濃度（ＩＣ５０）の化学物質が添加された場合の細胞の形態の時間変化に基づいて、化学物質同士の類似性について判定する場合について説明する。
本実施形態において用いられる培養観察装置を培養観察装置１１ａという。

図１５は、本実施形態に用いられる培養観察装置１１ａの概要を示すブロック図である。本実施形態に係る培養観察装置１１ａ（図１５）と第１の実施形態に係る培養観察装置１１（図１）とを比較すると、制御装置４１ａが異なる。しかし、他の構成要素（上部ケーシング１２、下部ケーシング１３、観察ユニット２２、及び表示部５０）が持つ機能は第１の実施形態と同じである。第１の実施形態と同じ機能の説明は省略し、第２の実施形態では、第１の実施形態と異なる部分を中心に説明する。

制御装置４１ａは、ＣＰＵ４２ａおよび記憶部４３を有している。ＣＰＵ４２ａは、制御装置４１ａの各種の演算処理を実行するプロセッサである。ＣＰＵ４２ａは、一例としてＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｉｏｎ）などの集積回路で構成される。なお、ＣＰＵ４２ａの一部、または全部を、ソフトウェアとして構成してもよい。つまり、ＣＰＵ４２ａは、ハードウェアとソフトウェアとを組み合わせて構成されてもよい。また、ＣＰＵ４２は、プログラムの実行によって、特徴量算出部４４ａ、座標算出部４５ａ、時間軌跡算出部４６ａ、時間軌跡定量化部４７ａ、及び類似化学物質判定部４８ａとしてそれぞれ機能する。

≪化学物質のスクリーニングの動作の例≫
以下、図１６のフローチャートを参照しつつ、化学物質のスクリーニングの動作の一例を説明する。図１６は、本実施形態に用いられる制御装置４１ａによる既知の化学物質の解析及び未知の化学物質の解析の動作の一例を示すフローチャートである。

ステップＳ４０１：制御装置４１ａは、既知の化学物質の解析を行う。ここでステップＳ４０１の処理は、図６のステップＳ２０１からステップＳ２０９までの処理と同様の処理を経た後に行われる。特徴量算出部４４ａは、抽出した各細胞の画像について、特徴量をそれぞれ算出する。特徴量算出部４４ａは、細胞が撮像された経過時間に対応する顕微鏡画像ごと特徴量を算出する。特徴量算出部４４ａが算出する特徴量の種類には、各細胞の形態についての特徴量、細胞群の集団としての形態についての特徴量が含まれる。
各細胞の形態についての特徴量の種類には、細胞の面積、長さ、幅、長さと幅との比率、周囲長などの約５５種類の指標が含まれる。細胞群の集団としての形態についての特徴量には、細胞の分布の広がり、細胞の密集度などの約６０種類の指標が含まれる。

座標算出部４５ａは、特徴量算出部４４ａが算出した特徴量に関連付けられている、特徴量の種類、化学物質の種類、化学物質の濃度、及び算出対象の顕微鏡画像の撮像時刻ごとに、特徴量空間における特徴量点の座標を算出する。

ここで時間軌跡算出部４６ａは、座標算出部４５ａが座標を算出した特徴量点から、特徴量点の時間軌跡及び濃度軌跡を算出してもよい。ここで時間軌跡とは、特徴量空間内の特徴量点について、細胞に添加された化学物質の濃度ごとの特徴量点同士を、顕微鏡画像の撮像時刻の順に線分または曲線を用いて接続して得られる軌跡である。つまり、時間軌跡算出部４６ａは、特徴量算出部４４ａが算出した特徴量点を、細胞が撮像された経過時間の順に並べてできる軌跡である時間軌跡を算出する。また、濃度軌跡とは、特徴量空間内の特徴量点について、顕微鏡画像の撮像時刻ごとの特徴量点同士を、細胞に添加された化学物質の濃度の順に線分または曲線を用いて接続して得られる軌跡である。

ここで図１７を参照し、濃度軌跡及び時間軌跡の具体例について説明する。濃度軌跡及び時間軌跡は、高次元の特徴量空間において算出される。ただし、図１７では、視覚化のために、高次元の特徴量空間における特徴量点が、主成分分析により特徴量平面における射影特徴量点に射影された一例について示している。後述する図１８、図１９、及び図２１においても同様に、高次元の特徴量空間における特徴量点が、主成分分析により特徴量平面における射影特徴量点に射影された一例について示している。
また、図１７、図１８、図１９、及び図２１においては、時間軌跡算出部４６ａが特徴量点の時間軌跡及び濃度軌跡を算出する場合の一例について示しているが、時間軌跡算出部４６ａは時間軌跡及び濃度軌跡を算出しなくてもよい。つまり、図１７、図１８、図１９、及び図２１において、特徴量点同士は、線分または曲線を用いて接続されなくてもよい。

図１７は、本実施形態に係る高濃度側の濃度軌跡及び５０％阻害濃度における時間軌跡の一例を示す図である。図１７には、ある細胞に５０％阻害濃度（ＩＣ５０）の１種類の化学物質が添加された場合の時間軌跡が示されている。この時間軌跡は、主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面（特徴量平面）に示されている。この時間軌跡は、プロットＩＣ５０−ｔからなる。ただし、プロットＩＣ５０−ｔの添字ｔは、顕微鏡画像の撮像時刻を、時間の順に示す。すなわち、プロットＩＣ５０−１は、最も古い撮像時刻に対する５０％阻害濃度（ＩＣ５０）の濃度に対応する射影特徴量点を示す。また、プロットＩＣ５０−５は、最も新しい撮像時刻に対する５０％阻害濃度（ＩＣ５０）の濃度に対応する射影特徴量点を示す。

図１７には、５０％阻害濃度（ＩＣ５０）の濃度よりも高い濃度に対応する射影特徴量点である高濃度プロットＰＨ１〜高濃度プロットＰＨ３が示されている。高濃度プロットＰＨ１〜高濃度プロットＰＨ３に対応する顕微鏡画像の撮像時刻は、プロットＩＣ５０−３に対応する顕微鏡画像の撮像時刻と同じ撮像時刻である。高濃度プロットＰＨ１〜高濃度プロットＰＨ３は、高濃度プロットＰＨ１、高濃度プロットＰＨ２、高濃度プロットＰＨ３の順に、対応する化学物質の濃度は高くなる。
高濃度プロットＰＨ１〜高濃度プロットＰＨ３に対しても、プロットＩＣ５０−３と同様に、プロットＩＣ５０−１、プロットＩＣ５０−２、プロットＩＣ５０−４、及びプロットＩＣ５０−５に対応する顕微鏡画像の撮像時刻に対応する時間軌跡（高濃度プロットＰＨ１―ｔ〜高濃度プロットＰＨ３―ｔ）が算出される。
プロットＩＣ５０−３、高濃度プロットＰＨ１、高濃度プロットＰＨ２、及び高濃度プロットＰＨ３は、高濃度側線分ＬＨ１〜高濃度側線分ＬＨ３により接続され濃度軌跡を構成する。

図１８は、本実施形態に係る低濃度側の濃度軌跡及び５０％阻害濃度における時間軌跡の一例を示す図である。図１８には、図１７において示した時間軌跡（プロットＩＣ５０−ｔ）に加え、５０％阻害濃度（ＩＣ５０）の濃度よりも低い濃度に対応するプロットである低濃度プロットＰＬ１〜低濃度プロットＰＬ３が示されている。低濃度プロットＰＬ１〜低濃度プロットＰＬ３は、主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面（特徴量平面）に示されている。低濃度プロットＰＬ１〜低濃度プロットＰＬ３に対応する顕微鏡画像の撮像時刻は、プロットＩＣ５０−３に対応する顕微鏡画像の撮像時刻と同じ撮像時刻である。低濃度プロットＰＬ１〜低濃度プロットＰＬ３は、低濃度プロットＰＬ１、低濃度プロットＰＬ２、低濃度プロットＰＬ３の順に、対応する化学物質の濃度は低くなる。
低濃度プロットＰＬ１〜低濃度プロットＰＬ３に対しても、プロットＩＣ５０−３と同様に、プロットＩＣ５０−１、プロットＩＣ５０−２、プロットＩＣ５０−４、及びプロットＩＣ５０−５に対応する顕微鏡画像の撮像時刻に対応する時間軌跡（低濃度プロットＰＬ１―ｔ〜低濃度プロットＰＬ３―ｔ）が算出される。
プロットＩＣ５０−３、低濃度プロットＰＬ１、低濃度プロットＰＬ２、及び低濃度プロットＰＬ３は、低濃度側線分ＬＬ１〜低濃度側線分ＬＬ３により接続され濃度軌跡を構成する。

図１６に戻って化学物質のスクリーニングの動作の一例の説明を続ける。
ステップＳ４０２：時間軌跡定量化部４７ａは、座標算出部４５ａが算出した特徴量空間における特徴量点について、細胞に添加された化学物質の濃度ごとの特徴量点の定量化を行う。また、時間軌跡定量化部４７ａは、座標算出部４５ａが算出した特徴量空間における特徴量点について、顕微鏡画像の撮像時刻ごとの特徴量点の定量化を行う。ここで、細胞に添加された化学物質の濃度ごとの特徴量点の定量化とは、一例として、５０％阻害濃度における特徴量点間の距離の算出である。顕微鏡画像の撮像時刻ごとの特徴量点の定量化とは、一例として、５０％阻害濃度より高濃度側の特徴量点間の距離の和の算出、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和の算出である。５０％阻害濃度より高濃度側の特徴量点間の距離の和、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和の算出は、それぞれ顕微鏡画像の撮像時刻ごとに算出される。特徴量点間の距離とは、一例として、ユークリッド距離やマハラビノス距離である。

図１９は、本実施形態に係る第１の化学物質の５０％阻害濃度における時間軌跡の一例を示す図である。この時間軌跡は、主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面（特徴量平面）に示されている。第１の化学物質とは、一例として、Ｐａｃｌｉｔａｘｅｌである。第１の化学物質の５０％阻害濃度における時間軌跡は、プロットＰＴ１〜プロットＰＴ５の５つの特徴量点からなる。

図２０は、本実施形態に係る第１の化学物質の低濃度側の特徴量点間の距離の和と高濃度側の特徴量点間の距離の和との時間変化の一例を示す図である。図２０では、低濃度側の線分の和と高濃度側の線分の和とが撮像時刻に対して示されている。時間軌跡定量化部４７ａは、低濃度側線分和ＳＬ１及び高濃度側線分和ＳＨ１を、５つの時刻ｔ１〜時刻ｔ５ごとに算出する。

図２１は、本実施形態に係る第２の化学物質の５０％阻害濃度における時間軌跡の一例を示す図である。図２１では、時間軌跡が主成分分析の結果得られた軸ＰＣ１及び軸ＰＣ２を座標軸とする２次元平面（特徴量平面）に示されている。第２の化学物質とは、一例として、Ｅｔｏｐｏｓｉｄｅである。第２の化学物質の５０％阻害濃度における時間軌跡は、プロットＥＴ１〜プロットＥＴ５の５つのプロットからなる。

図２２は、本実施形態に係る第２の化学物質の低濃度側の特徴量点間の距離の和と高濃度側の特徴量点間の距離の和との時間変化の一例を示す図である。図２２では、低濃度側の特徴量点間の距離の和と高濃度側の特徴量点間の距離の和とが撮像時刻に対して示されている。時間軌跡定量化部４７ａは、低濃度側線分和ＳＬ２及び高濃度側線分和ＳＨ２を、５つの時刻ｔ１〜時刻ｔ５ごとに算出する。

図１６に戻って化学物質のスクリーニングの動作の一例の説明を続ける。
ステップＳ４０３：時間軌跡定量化部４７ａは、ステップＳ４０２において算出した定量化されたデータをデータベースに登録する。ステップＳ４０２において算出した定量化されたデータとは、５０％阻害濃度における特徴量点間の距離の和、５０％阻害濃度より高濃度側の特徴量点間の距離の和、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和である。データベースとは、例えば、記憶部４３である。

ステップＳ４０４：制御装置４１ａは、ステップＳ４０１において既知の化学物質に対して行ったのと同様に、未知の化学物質の解析を行う。具体的には、制御装置４１ａは、未知の化学物質について、特徴量の分布に基づいて特徴量空間に特徴量点をプロットする。また制御装置４１は、解析結果から得られる時間軌跡及び濃度軌跡の算出を行ってもよい。
ステップＳ４０５：時間軌跡定量化部４７ａは、ステップＳ４０２において既知の化学物質に対して行ったのと同様に、未知の化学物質について座標算出部４５ａが算出した特徴量空間における特徴量点について、細胞に添加された化学物質の濃度ごとの特徴量点の定量化を行う。

ステップＳ４０６：類似化学物質判定部４８ａは、未知の化学物質が既知の化学物質のいずれと最も類似しているかについての判定を行う。本実施形態では、類似化学物質判定部４８ａは、未知の化学物質が、既知の化学物質である第１の化学物質（Ｐａｃｌｉｔａｘｅｌ）と、既知の化学物質である第２の化学物質（Ｅｔｏｐｏｓｉｄｅ）とのいずれと最も類似しているかについての判定を行う。つまり、本実施形態の化学物質のスクリーニング方法では、異なる化学物質は、既知の化学物質と、未知の化学物質とが含まれ、既知の化学物質と、未知の化学物質との類似性を判定する。類似化学物質判定部４８ａの未知の化学物質の判定の動作の詳細は、図２３のフローチャートにおいて説明する。

図２３は、本実施形態に係る制御装置４１ａによる未知の化学物質の判定の動作の一例を示すフローチャートである。
ステップＳ６１：類似化学物質判定部４８ａは、ステップＳ４０５において算出した未知の化学物質の定量化されたデータと、ステップＳ４０３においてデータベースに登録された既知の化学物質の定量化されたデータとを比較する。類似化学物質判定部４８ａは、３つの指標を算出することにより比較を行う。ここで３つの指標とは、ＩＣ５０距離、低濃度側差分、及び高濃度側差分である。

ＩＣ５０距離とは、特徴量空間において、未知の化学物質に対する５０％阻害濃度における特徴量点と、既知の化学物質に対する５０％阻害濃度における特徴量点との距離を、顕微鏡画像の撮像時刻が同じ特徴量点間において算出したものである。
低濃度側差分ＰＤＬとは、未知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和との差分が顕微鏡画像の各撮像時刻ごとにおいて算出されたものである。
高濃度側差分ＰＤＨとは、未知の化学物質に対する５０％阻害濃度より高濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より濃度軌跡の高濃度側の特徴量点間の距離の和との差分が顕微鏡画像の各撮像時刻ごとにおいて算出されたものである。

図２４は、本実施形態に係る未知の化学物質と第１の化学物質との比較の一例を示す図である。図２４では、ＩＣ５０距離ＰＤＩ、低濃度側差分ＰＤＬ、及び高濃度側差分ＰＤＨが撮像時刻に対して示されている。類似化学物質判定部４８ａは、ＩＣ５０距離ＰＤＩ、低濃度側差分ＰＤＬ、及び高濃度側差分ＰＤＨを算出する。

図２５は、本実施形態に係る未知の化学物質と第２の化学物質との比較の一例を示す図である。図２５では、ＩＣ５０距離ＥＤＩ、低濃度側差分ＥＤＬ、及び高濃度側差分ＥＤＨが撮像時刻に対して示されている。類似化学物質判定部４８ａは、ＩＣ５０距離ＥＤＩ、低濃度側差分ＥＤＬ、及び高濃度側差分ＥＤＨを算出する。

図２３に戻って、制御装置４１ａによる未知の化学物質の判定の動作の説明を続ける。 ステップＳ６２：類似化学物質判定部４８ａは、ステップＳ６１において算出した３つの指標の重みづけの選択を行う。類似化学物質判定部４８ａは３つの指標の各々に対する重みを選択する。例えば、類似化学物質判定部４８ａは、３つの指標のうちいずれか１つの指標のみを未知の化学物質の判定に用いる場合、この１つの指標の重みを１とし、他の２つの指標の重みをゼロとする。また、類似化学物質判定部４８ａは、３つの指標の平均を未知の化学物質の判定に用いる場合、３つの指標の重みを全て１とする。類似化学物質判定部４８ａは、３つの指標の重みを各々任意に選択してよい。

ステップＳ６３：類似化学物質判定部４８ａは、既知の化学物質の中から未知の化学物質と最も類似する化学物質を判定する。ここで類似化学物質判定部４８ａは、多変量解析により解析された結果に基づく、所定の濃度における特徴量の複数の時刻毎の大きさの変化と、所定の濃度よりも高濃度である場合の特徴量の複数の時刻毎の大きさの変化と、所定の濃度よりも低濃度である場合の特徴量の複数の時刻毎の大きさの変化とにそれぞれに重みづけされた結果に基づく異なる化学物質同士の類似性を判定する。

具体的には、類似化学物質判定部４８ａは、ステップＳ６１において算出した３つの指標の各々について、各撮像時刻ごとに値の二乗を算出し、算出した二乗の値の和を算出する。類似化学物質判定部４８ａは、３つの指標の各々について、ステップＳ６２において選択した重みを、算出した二乗の値に乗じ、和を算出する。類似化学物質判定部４８ａは、算出した和が最も小さい化学物質を、未知の化学物質と最も類似する化学物質であると判定する。なお、類似化学物質判定部４８ａは、二乗の値の和の代わりに、二乗以外の冪乗の和を算出してもよいし、絶対値の和を算出してもよい。
図２４及び図２５に示す例では、類似化学物質判定部４８ａは、第１の化学物質（Ｐａｃｌｉｔａｘｅｌ）と第２の化学物質（Ｅｔｏｐｏｓｉｄｅ）のうち第１の化学物質（Ｐａｃｌｉｔａｘｅｌ）を未知の化学物質と最も類似する化学物質であると判定する。

このように、類似化学物質判定部４８ａは、所定の濃度における、化学物質が細胞に添加された後の複数の撮像時刻の変化に対する特徴量の変化に基づいて、異なる化学物質同士の類似性を判定する。ここで所定の濃度は、５０％阻害濃度、５０％阻害濃度よりも高い濃度、５０％阻害濃度よりも低い濃度の少なくとも一方の濃度を含む。

また、類似化学物質判定部４８ａは、５０％阻害濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも高い濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも低い濃度における撮像時刻の変化に対する特徴量の変化とにそれぞれに重み付けされた結果に基づいて異なる化学物質同士の類似性を判定する。

なお、本実施形態においては、ステップＳ６１において類似化学物質判定部４８ａが、低濃度側差分ＰＤＬとして、顕微鏡画像の各撮像時刻ごとにおいて、未知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和との差分を算出する場合の例について説明したが、差分に限らない。例えば、類似化学物質判定部４８ａは、未知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より低濃度側の特徴量点間の距離の和との割合を算出してもよいし、その他の既知の演算を行ってもよい。また、同様に類似化学物質判定部４８ａは、未知の化学物質に対する５０％阻害濃度より高濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より高濃度側の特徴量点間の距離の和と差分に限らず、未知の化学物質に対する５０％阻害濃度より高濃度側の特徴量点間の距離の和と、既知の化学物質に対する５０％阻害濃度より高濃度側の特徴量点間の距離の和との割合を算出してもよいし、その他の既知の演算を行ってもよい。
また、本実施形態においては、ステップＳ６２において類似化学物質判定部４８ａが３つの指標の重みづけの選択を行う場合の例について説明したが、ステップＳ６２の処理は省略されてよい。ステップＳ６２の処理が省略される場合、ステップＳ６２において、類似性判定部４８は、例えば、３つの指標の各々について二乗の値の和を算出してよい。 また、本実施形態においては、ステップＳ６３において類似化学物質判定部４８ａが、算出した和が最も小さい化学物質を、未知の化学物質と最も類似する化学物質であると判定する場合について説明したが、これに限らない。類似化学物質判定部４８ａは、算出した和が所定の閾値よりも大きい場合、判定対象の化学物質は未知の化学物質と類似していないと判定してもよい。

なお、本実施形態においては、ステップＳ４０２及びステップＳ４０５において時間軌跡定量化部４７ａが特徴量空間における特徴量点の定量化として、５０％阻害濃度より高濃度側の特徴量点間の距離の和の算出、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和の算出の各処理を行う場合の例について説明したが、これに限らない。ステップＳ４０２において時間軌跡定量化部４７ａは、特徴量空間における特徴量点の定量化として、５０％阻害濃度より高濃度側の特徴量点間の距離の和の算出、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和の算出の少なくとも１つの処理を行ってもよい。ステップＳ４０２及びステップＳ４０５において時間軌跡定量化部４７ａが、特徴量空間における特徴量点の定量化として、５０％阻害濃度より高濃度側の特徴量点間の距離の和の算出、及び５０％阻害濃度より低濃度側の特徴量点間の距離の和の算出の少なくとも１つの処理を行う場合、時間軌跡定量化部４７ａが行った特徴量空間における特徴量点の定量化の処理に応じて、ステップＳ４０６において、類似化学物質判定部４８ａは、ＩＣ５０距離、低濃度側差分、及び高濃度側差分の３つの指標のうち少なくとも１つの指標に基づいて、未知の化学物質が既知の化学物質のいずれと最も類似しているかについての判定を行う。また、特徴量点間の距離の算出において５０％阻害濃度を基準としなくてもよい。例えば、時間軌跡定量化部４７ａは、所定の濃度を基準として、当該所定の濃度より高濃度側の特徴量点間の距離の和、及び／または低濃度側の特徴量点間の距離の和を算出し、特徴量点の定量化を行ってもよい。また、ステップＳ４０２及びステップＳ４０５において時間軌跡定量化部４７ａは、５０％阻害濃度における特徴量点間の距離（ＩＣ５０距離）を算出する場合の例について説明したが、これに限らない。ステップＳ４０２及びステップＳ４０５において時間軌跡定量化部４７ａは、５０％阻害濃度以外の濃度について顕微鏡画像の撮像時刻ごとの特徴量点間の距離を算出してもよい。

≪まとめ≫
以上説明したように、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法は、顕微鏡画像は、化学物質が細胞に添加された後の複数の撮像時刻において撮像され、所定の濃度における、撮像時刻の変化に対する特徴量の変化に基づいて、異なる化学物質同士の類似性を判定する。
この構成により、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法では、細胞を固定、染色することなく、安価で迅速な解析により化学物質が細胞に生じさせる時間についての形態変化を所定の濃度である場合に定量化できるため、細胞を固定、染色することなく、安価で迅速な解析により細胞に所定の濃度である場合に形態変化を生じさせる化学物質同士の類似性を判断することができる。

また、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法は、所定の濃度は、５０％阻害濃度、前記５０％阻害濃度よりも高い濃度、前記５０％阻害濃度よりも低い濃度の少なくとも一方の濃度を含む。
この構成により、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法では、細胞を固定、染色することなく、安価で迅速な解析により化学物質が細胞に生じさせる時間についての形態変化を５０％阻害濃度、前記５０％阻害濃度よりも高い濃度、前記５０％阻害濃度よりも低い濃度の少なくとも一方の濃度を含む所定の濃度に基づいて定量化できるため、安価で迅速な解析により、細胞に添加される化学物質の濃度が５０％阻害濃度、前記５０％阻害濃度よりも高い濃度、前記５０％阻害濃度よりも低い濃度の少なくとも一方の濃度である場合に形態変化を生じさせる化学物質同士の類似性を判断することができる。

また、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法は、５０％阻害濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも高い濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも低い濃度における撮像時刻の変化に対する特徴量の変化とにそれぞれに重み付けされた結果に基づいて異なる化学物質同士の類似性を判定する。
この構成により、本実施形態の制御装置４１ａが提供する化学物質のスクリーニング方法では、細胞に形態変化を生じさせる濃度に応じて、５０％阻害濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも高い濃度における撮像時刻の変化に対する特徴量の変化と、５０％阻害濃度よりも低い濃度における撮像時刻の変化に対する特徴量の変化とにそれぞれ重みづけをして化学物質同士の類似性を判断することができるため、重みづけをしない場合に比べ化学物質のスクリーニングの精度を上げることができる。

本発明の化学物質のスクリーニング方法は、非破壊の系であるため、タイムラプス解析に好適に用いられ、かつスクリーニングに用いたサンプルそのものを別の評価系に供することができ、スクリーニング結果を別の系から評価することができる。

なお、本発明の実施形態における培養観察装置１１及び培養観察装置１１ａの各処理を実行するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、当該記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、上述した種々の処理を行ってもよい。

なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものであってもよい。また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、フラッシュメモリ等の書き込み可能な不揮発性メモリ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。

さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（例えばＤＲＡＭ（Ｄｙｎａｍｉｃ Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ））のように、一定時間プログラムを保持しているものも含むものとする。また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。

以上、本発明の実施形態について図面を参照して詳述したが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。

１１、１１ａ…培養観察装置、４１、４１ａ…制御装置、４４、４４ａ…特徴量算出部、４５、４５ａ…座標算出部、４６…濃度軌跡算出部、４７…変化領域判定部、４８…類似性判定部、４６ａ…時間軌跡算出部、４７ａ…時間軌跡定量化部、４８ａ…類似化学物質判定部、５０…表示部

Claims (22)

1. 異なる化学物質がそれぞれ異なる濃度で添加された複数の細胞群における細胞の顕微鏡画像に含まれる前記細胞の形態に関する複数の特徴量を取得することと、
   前記細胞の形態に関する複数の特徴量に基づいて、前記異なる化学物質同士の類似性を判定することと、
   を有する化学物質のスクリーニング方法。
2. 前記判定することは、前記異なる特徴量の内、前記細胞の形態が変化する形態変化濃度に対応する特徴量を前記異なる化学物質同士で比較して、前記類似性を判定する請求項１に記載の化学物質のスクリーニング方法。
3. 前記形態変化濃度は、前記細胞の形態が変化してから前記細胞が細胞死するまでの濃度を含む請求項２に記載の化学物質のスクリーニング方法。
4. さらに、前記濃度の変化に対する前記特徴量の変化に基づいて前記形態変化濃度を求めることを有する請求項２または請求項３に記載の化学物質のスクリーニング方法。
5. 前記形態変化濃度を求めることは、前記特徴量の変化量と閾値とを比較して前記形態変化濃度を求める請求項４に記載の化学物質のスクリーニング方法。
6. 前記顕微鏡画像は、前記化学物質が前記細胞に添加された後の複数の撮像時刻において撮像され、
   所定の前記濃度における、前記撮像時刻の変化に対する前記特徴量の変化に基づいて、前記異なる化学物質同士の前記類似性を判定する請求項１から請求項５のいずれか一項に記載の化学物質のスクリーニング方法。
7. 前記所定の濃度は、５０％阻害濃度、前記５０％阻害濃度よりも高い濃度、前記５０％阻害濃度よりも低い濃度の少なくとも一方の濃度を含む請求項６に記載の化学物質のスクリーニング方法。
8. 前記５０％阻害濃度における前記撮像時刻の変化に対する前記特徴量の変化と、前記５０％阻害濃度よりも高い濃度における前記撮像時刻の変化に対する前記特徴量の変化と、前記５０％阻害濃度よりも低い濃度における前記撮像時刻の変化に対する前記特徴量の変化とにそれぞれに重み付けされた結果に基づいて前記類似性を判定する請求項７に記載の化学物質のスクリーニング方法。
9. 前記異なる化学物質は、既知の化学物質と、未知の化学物質とが含まれ、
   前記既知の化学物質と、前記未知の化学物質との前記類似性を判定する
   請求項１から請求項８のいずれか一項に記載の化学物質のスクリーニング方法。
10. 前記複数の特徴量が多変量解析により解析された結果に基づいて前記類似性を判定する請求項１から請求項９のいずれか一項に記載の化学物質のスクリーニング方法。
11. コンピュータに、
    請求項１から請求項１０のいずれか一項に記載の化学物質のスクリーニング方法に基づくステップを実行させるプログラム。
12. 異なる化学物質がそれぞれ異なる濃度で添加された複数の細胞群における細胞の顕微鏡画像に含まれる前記細胞の形態に関する複数の特徴量を算出する特徴量算出部と、
    前記細胞の形態に関する複数の特徴量に基づいて、前記異なる化学物質同士の類似性を判定する類似性判定部と、
    を備える制御装置。
13. 前記類似性判定部は、前記異なる特徴量の内、前記細胞の形態が変化する形態変化濃度に対応する特徴量を前記異なる化学物質同士で比較して、前記類似性を判定する請求項１２に記載の制御装置。
14. 前記形態変化濃度は、前記細胞の形態が変化してから前記細胞が細胞死するまでの濃度を含む請求項１３に記載の制御装置。
15. さらに、前記濃度の変化に対する前記特徴量の変化に基づいて前記形態変化濃度を求める濃度判定部を備える請求項１３または請求項１４に記載の制御装置。
16. 前記濃度判定部は、前記特徴量の変化量と閾値とを比較して前記形態変化濃度を求める請求項１５に記載の制御装置。
17. 前記顕微鏡画像は、前記化学物質が前記細胞に添加された後の複数の撮像時刻において撮像され、
    前記類似性判定部は、所定の前記濃度における、前記撮像時刻の変化に対する前記特徴量の変化に基づいて、前記異なる化学物質同士の前記類似性を判定する請求項１２から請求項１６のいずれか一項に記載の制御装置。
18. 前記所定の濃度は、５０％阻害濃度、前記５０％阻害濃度よりも高い濃度、前記５０％阻害濃度よりも低い濃度の少なくとも一方の濃度を含む請求項１７に記載の制御装置。
19. 前記類似性判定部は、前記５０％阻害濃度における前記撮像時刻の変化に対する前記特徴量の変化と、前記５０％阻害濃度よりも高い濃度における前記撮像時刻の変化に対する前記特徴量の変化と、前記５０％阻害濃度よりも低い濃度における前記撮像時刻の変化に対する前記特徴量の変化とにそれぞれに重み付けされた結果に基づいて前記類似性を判定する請求項１８に記載の制御装置。
20. 前記異なる化学物質は、既知の化学物質と、未知の化学物質とが含まれ、
    前記既知の化学物質と、前記未知の化学物質との前記類似性を判定する
    請求項１２から請求項１９のいずれか一項に記載の制御装置。
21. 前記類似性判定部は、前記複数の特徴量が多変量解析により解析された結果に基づいて前記類似性を判定する請求項１２から請求項２０のいずれか一項に記載の制御装置。
22. 請求項１２から請求項２１のいずれか一項に記載の制御装置を備える培養観察装置。

Images