(57)【特許請求の範囲】
【請求項１】
細胞の局在を評価する画像解析装置であって、
第１の領域内の評価対象細胞群の２次元領域内における細胞群の位置と形状とからなる２次元パターンを抽出するパターン抽出部と、
前記パターン抽出部により抽出された第１の時点における前記評価対象細胞群の位置と形状と前記第１の時点とは異なる第２の時点における前記評価対象細胞群の位置と形状との相関度を求める相関解析を行う相関解析部と、
前記相関解析部により求められた第１の時点における前記評価対象細胞群の位置と形状と第２の時点における評価対象細胞群の位置と形状との相関度に基づいて、前記評価対象細胞群の組織内での自己組織化過程を数値化する位置−時間演算部と、を有することを特徴とする画像解析装置。
【請求項２】
前記相関度に基づいて、前記評価対象細胞群の局在の変化の度合いを判断する判断部を有し、
前記判断部は、前記評価対象細胞群の位置と形状に依存する前記相関度のしきい値に基づいて前記局在の変化の度合いを判断することを特徴とする請求項１に記載の画像解析装置。
【請求項３】
前記相関度に基づいて、前記評価対象細胞群の局在の変化の度合いを判断する判断部を有し、
前記判断部は、第１の試料の前記評価対象細胞群の相関度をしきい値として設定し、前記第１の試料とは異なる第２の試料の評価対象細胞群の相関度を判定する請求項１に記載の画像解析装置。
【請求項４】
前記評価対象細胞群の位置と形状が血管内皮細胞による制御されていないパターンの場合に、前記しきい値を、0.1から0.2に設定する請求項２に記載の画像解析装置。
【請求項５】
前記評価対象細胞群の位置と形状が血管内皮細胞による制御されたパターンの場合に、前記しきい値を、0.6から0.8に設定する請求項２に記載の画像解析装置。
【請求項６】
前記パターン抽出部は、
取得した前記第１の領域内の前記評価対象細胞群の画像を二値化する二値化処理部と、  ノイズを除去するノイズ除去処理部と、を有することを特徴とする請求項１から５までのいずれか１項に記載の画像解析装置。
【請求項７】
さらに、前記第１の時点と前記第２の時点とにおける前記評価対象細胞群の位置と形状と、前記第１及び第２の時点における評価対象細胞群の位置と形状と対応する相関解析結果との対応関係を記憶する記憶部を有することを特徴とする請求項１から６までのいずれか１項に記載の画像解析装置。
【請求項８】
前記評価対象細胞群が、人工的に作製された細胞組織であることを特徴とする請求項１から７までのいずれか１項に記載の画像解析装置。
【請求項９】
細胞の局在を評価する画像解析方法であって、
第１の領域内の評価対象細胞群の２次元領域内における前記評価対象細胞群の位置と形状とからなる２次元パターンを抽出するパターン抽出ステップと、
前記パターン抽出ステップにより抽出された第１の時点における前記評価対象細胞群の位置と形状と前記第１の時点とは異なる第２の時点における前記評価対象細胞群の位置と形状との相関度を求める相関解析を行う相関解析ステップと、
前記相関解析ステップにより求められた第１の時点における前記評価対象細胞群の位置と形状と第２の時点における前記評価対象細胞群の位置と形状との相関度に基づいて、前記評価対象細胞群の組織内での自己組織化過程を数値化する位置−時間演算ステップと、を有することを特徴とする画像解析方法。
【請求項１０】
コンピュータに、請求項９に記載の画像解析方法を実行させるための画像解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像解析技術に関し、特に、評価対象細胞群の位置と形状との局在の変化の指標を求める細胞画像解析技術に関する。
【背景技術】
【０００２】
近年、生体の組織や臓器を模倣した、細胞が高密度に集合している三次元組織モデルを作製し、創薬や再生医療領域へ応用することが期待されている。このような三次元組織内では様々な種類の細胞が適切な位置関係に配置されることで効率的に組織の機能を発現している。そこで、温度応答性表面での細胞パターニングなどの技術を用いて三次元組織内における任意の細胞の局在を制御することで、細胞の配置が三次元的にデザインされた組織を構築する試みが多くなされている。
【０００３】
下記特許文献１による細胞シートを利用した細胞評価システム及びその利用方法では、三次元組織内での三次元的な細胞の遊走および局在の変化を、２次元解析装置により画像を観察しながら追跡している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】国際公開2010/101225号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
一方で、in vitroで作製した三次元組織内において血管内皮細胞などの細胞は非常に活発に組織内で遊走（移動）することが知られている。これにより、組織構築時に任意の細胞の局在を制御した組織を作製したとしても、培養期間を経るごとに組織内で細胞が遊走し局在がランダムに変化してしまう。この三次元化後の細胞の遊走を防ぐために様々な手法の研究が試みられているが、これまで組織内での細胞の局在変化を数値化し定量的に評価することは煩雑で困難な手順を踏む必要があったため、各手法を評価するための研究において障壁となっていた。
【０００６】
特許文献１に記載の技術では、ある平面における細胞群の変化を観察しているが、細胞群の変化に明確な指標がないため、客観的な評価をしにくいという課題がある。
本発明は、細胞群の位置と形状の局在の変化を客観的に評価することを目的とする。
【課題を解決するための手段】
【０００７】
本発明においては、解析可能な領域をあるひとつの平面に絞り、画像相関解析を用いる。これにより、一定時間後の細胞群の形の変化を数値化して捉えることができ、細胞群の変化に明確な指標を得ることができる。
【０００８】
本発明では三次元組織内の血管内皮細胞を始めとする任意の細胞をあらかじめ蛍光標識などによりラベリングした状態でタイムラプス画像を取得し、これらの画像を画像相関解析法により解析することで簡便に数値化し定量的に評価するシステムを構築した。これにより、三次元組織内の任意の細胞の遊走挙動や局在の乱れを定量的に評価可能な汎用性の高い技術を確立することができる。
【０００９】
本発明の一観点によれば、細胞の局在を評価する画像解析装置であって、第１の領域内の評価対象細胞群の２次元パターンを抽出するパターン抽出部と、前記パターン抽出部により抽出された第１の時点における第１の局在パターン形状と前記第１の時点とは異なる第２の時点における第２の局在パターン形状との相関解析を行う相関解析部と、を有することを特徴とする画像解析装置が提供される。
【００１０】
さらに、前記相関解析部により求められた前記第１の局在パターン形状と前記第２の局在パターン形状との相関度を、前記評価対象細胞群の位置と形状との局在の変化の指標として求めることにより前記評価対象細胞群の組織内での自己組織化過程を数値化することができる。
【００１１】
本発明の他の観点によれば、細胞の局在を評価する画像解析方法であって、第１の領域内の評価対象細胞群の２次元パターンを抽出するパターン抽出ステップと、前記パターン抽出ステップにより抽出された第１の時点における第１の局在パターン形状と前記第１の時点とは異なる第２の時点における第２の局在パターン形状との相関解析を行う相関解析ステップと、を有することを特徴とする画像解析方法が提供される。
【００１２】
本発明は、コンピュータに、上記に記載の画像解析方法を実行させるための画像解析プログラムであっても良く、当該プログラムを記録するコンピュータ読み取り可能な記録媒体であっても良い。
【発明の効果】
【００１３】
本発明によれば、細胞群の位置と形状の局在の変化を客観的に評価できる指標を得ることができる。
【図面の簡単な説明】
【００１４】
【図１】図１（ａ）は、評価対象細胞群として、シート状の細胞組織（細胞シート）を作成する方法の一例を示す図である。図１（ｂ）は、細胞シートを含む構成の断面図である。図１（ｃ）は、得られた評価対象試料の平面図である。
【図２】評価対象細胞群の画像の相関解析手順をイメージで示す図である。
【図３Ａ】評価対象細胞群の画像の相関解析処理の流れの一例を示すフローチャート図である。
【図３Ｂ】図３Ｂ（ａ）は、画像解析装置の一構成例を示す機能ブロック図である。図３Ｂ（ｂ）は、パターン抽出部の位置構成例を示す機能ブロック図である。
【図４】ランダムパターン形状を有する異なる試料について相関解析を行った結果の例を示す図である。
【図５】図４の相関解析結果を相関度により比較した図である。
【図６】図４の相関解析結果を相関度の時間依存性として示す図である。
【図７】ラインパターン形状を有する試料の作成方法の一例を示す図である。
【図８】ラインパターン形状を有する同じ試料の画像の経時変化を示す図である。
【図９】図８に示すラインパターン形状を有する同じ試料の相関度の経時変化をグラフで示す図である。
【図１０】ライン間隔の異なるラインパターン形状を有する異なる試料の画像の経時変化を示す図である。
【図１１】ライン間隔の異なるラインパターン形状を有する異なる試料の相関度のライン間隔依存性をグラフで示す図である。
【発明を実施するための形態】
【００１５】
本明細書において、解析対象となる細胞群（細胞集団）は、生体の組織又は器官から直接採取した初代細胞系に由来する細胞でもよく、該初代細胞を何代か継代させた継代細胞系に由来する細胞であってもよい。細胞が由来する動物も特に限定されず、例えば、ヒト、非ヒト霊長類、マウス、ラット、ウサギ、イヌ、ネコ、ウシ、ブタ、ウマ、ヤギ又はヒツジ等の動物に由来する細胞を使用することができる。
【００１６】
また、評価対象細胞群が「パターン形状が制御されていない」とは、例えばランダムパターンのように決められた形状を保持するように制御されていない状態、例えば細胞自身の自己組織化によりランダムな位置に局在している状態を指す。評価対象細胞群が「パターン形状が制御された」とは、例えばラインパターンのように決められた形状を保持するように制御されていないことを指す。
以下に、本発明の実施の形態による細胞の画像解析技術について図面を参照しながら詳細に説明する。
【００１７】
（細胞シート中での血管内皮細胞ネットワークの形成工程）
図１は、評価対象細胞の一例として、シート状の細胞組織（細胞シート）を作製する方法の一例を示す図である。本実施の形態では、緑色蛍光タンパク発現正常ヒト臍帯静脈内皮細胞(GFP-HUVEC) (Angio-proteomie, Boston, USA)２５ａ及び正常ヒト皮膚線維芽細胞(NHDF) (Lonza Inc., Walkersville, USA)２５ｂを用いて細胞シート２５を作成した。細胞シートの作製には温度応答性ポリマーであるPNIPAAm、ポリ-N-イソプロピルアクリルアミドを表面に修飾した温度応答性培養皿(f=35 mm, UpCellTM, CellSeed Inc., Tokyo, Japan)２１を培養液２３で満たし、総細胞数1×106cells/dish、血管内皮細胞：線維芽細胞の比率＝1：7の密度で細胞を播種しコンフルエントに達するまで培養を行なった（図１（ａ）、図１（ｂ））。
【００１８】
次いで、コンフルエントの状態に達した細胞を、20℃で20分以上静置し、GFP-HUVEC２５ａ及びNHDF２５ｂで構成されたシート状の細胞組織を回収した。シート状の細胞組織を細胞シートと称する。
次いで、細胞シートを収縮させて３７℃でインキュベーションを行うことで、GFP-HUVECが組織内でランダムに配置されたシート状の組織２５ｘ（図１（ｃ））を得ることができた。
【００１９】
（タイムラプス蛍光画像の取得方法）
回収したシート状の細胞組織２５ｘを、タイムラプス蛍光顕微鏡装置(BZ-9000, Keyence Corp., Osaka, Japan)を用い、１時間おきに５日間、同一箇所のタイムラプス蛍光画像を取得した。シート状の細胞組織は、endothelial growth medium (EGM-2) (Lonza Inc., Walkersville, USA) 中に50 ng/mL recombinant human vascular endothelial growth factor (rhVEGF) (R and D Systems, Minneapolis, USA)を添加し、抗生物質存在下で、顕微鏡ステージ装着型インキュベータ (INUG2-KI3, Tokai Hit Co. Ltd., Shizuoka, Japan)を用い、37℃、CO2濃度5%環境中で培養した。
【００２０】
図２、図３Ａ、図３Ｂは、上記の方法により作成した血管パターンを例にした評価対象細胞群の画像解析技術を具体的に示す図である。図２は、評価対象細胞群の画像の解析手順のイメージを示す図であり、図３Ａは、評価対象細胞群の画像の解析処理例を示すフローチャート図であり、図３Ｂは、評価対象細胞群の画像解析に用いた画像解析装置の構成例を示す機能ブロック図である。画像解析装置は、ソフトウェアによるかハードウェアによるかを問わない。
【００２１】
発明者らは、２次元領域内における細胞群の位置と形状の局在の変化を客観的に評価できる指標として、画像相関解析を用いることを考案した。相関度に基づいて、前記局在の変化の度合いを判断することができる。
【００２２】
以下に、シート状の細胞組織内の血管内皮細胞による血管ネットワークの形成過程を定量評価する手法を例にして、本実施の形態による細胞の画像相関解析について詳細に説明する。
【００２３】
図３Ｂ（ａ）に示すように、本実施の形態による画像解析装置１は、シート状の細胞組織２５ｘなどの撮像した観察画像を取得する画像取得部１−１と、画像取得部１−１により取得した細胞群のパターン形状を抽出するパターン抽出部１−２と、例えば、同じ領域を異なる時点、すなわち第１の時点と第１の時点とは異なる第２の時点とにおいて抽出した異なる２つのパターン形状の画像相関解析を行う相関解析部１−３と、相関度に基づいて細胞群の局在の時間変化を求める位置−時間演算部１−４と、相関度に基づいて局在の変化の度合いを判断する判断部１−５と、判断結果等を出力する出力部１−６と、を有している。さらに、判断部１−５において局在の変化の度合いを判断するための適切なしきい値を設定するしきい値設定部１−７と、局在パターンと時間との関係等を記憶する記憶部１−８とを備えていても良い。画像解析装置１はこれらすべての構成を一体として備えていなくても良い。画像相関解析には、ImageJ (National Institutes of Health, MD, USA)を用いた。
【００２４】
画像相関解析における相関度の算出は例えば以下のようにして行うことができる。すなわち、ある画素数（ｎ，ｍ）の画像ａと画像ｂの相関度の算出式は、下記の通りである。
ここで、画像ａの画素（ｘ，ｙ）の輝度値をａ（ｘ，ｙ）とする。画像ｂの画素（ｘ，ｙ）の輝度値をｂ（ｘ，ｙ）とする。
【００２５】
【数１】
0006661165.tif 000002
【００２６】
図３Ｂ（ｂ）は、図３Ｂ（ａ）のパターン抽出部１−２の一構成例を示す機能ブロック図である。図３Ｂ（ｂ）に示すように、パターン抽出部１−２は、画像のタイムラプス解析を行うタイムラプス解析部１−２−１と、タイムラプス解析した画像データを取得する画像データ取得部１−２−２と、取得した画像データの二値化を行う二値化部１−２−３と、二値化した画像のノイズ部分を除去するノイズ除去部１−２−４とを有している。このようにして画像解析処理に適した前処理を行って解析に適したパターンを抽出することで、２次元領域内における細胞群の位置と形状の局在の変化を客観的に評価できる指標を精度良く求めることができる。
【００２７】
画像解析処理の流れを図３Ａに示すフローチャートに基づいて説明する。ステップＳ１において処理が開始されると（Ｓｔａｒｔ）、ステップＳ２において、タイムラプス解析部１−２−１はタイムラプス解析を行う、次いで、ステップＳ３において、画像データ取得部１−２−２が２枚の画像を取得する。２枚の画像は、例えば、図２に示すような、評価対象細胞群の第１の領域における蛍光画像（Ａ）Ｐ１と、同じ第１の領域における蛍光画像（Ｂ）Ｐ２であり、これらは、血管パターンの蛍光画像のうち、それぞれ、初期画像と５日後の画像である。画像相関解析を行う２つの画像パターンは、第１の時点と第１の時点とは異なる第２の時点における画像であり、第１の時点は、細胞観察の初期時点であっても良く、それ以降の時点であっても良い。
【００２８】
次いで、ステップＳ４において、二値化部１−２−３が２枚の画像について二値化を行う。次いで、ステップＳ５において、二値化を行った画像のノイズ部分をノイズ除去部１−２−４が除去し、パターン抽出部１−２が全体としてステップＳ６でパターンを抽出する。図２のＰ１１とＰ１２とが、画像相関解析を行うための前処理を行った後の、同じ領域の異なる時間における画像である。
【００２９】
最後に、抽出したパターンを、ステップＳ７において、相関解析部１−３が画像相関解析を行うことで、２つの異なる画像の相関度を求める第１の解析処理が終了する（Ｅｎｄ）。画像相関解析の結果は、評価対象細胞群の位置と形状との局在の変化の指標となる。
【００３０】
尚、画像解析装置１は、第１の時点における第１の画像とは異なる第３の画像と、第２の時点における第２の画像との相関度を求めた第１の画像との相関度を演算する相関度関連性判定部を有していても良い。この相関度関連性判定部の判定結果として大きい（１に近似する）相関度が得られた場合には、第３の画像における異なる評価対象細胞群について、記憶部１−８に記憶された相関度と第１の時点から第２の時点までの経過時間との関係を、異なる評価対象細胞の位置の時間依存として近似するようにしても良い。記憶部を設けることで、類似の試料等における評価対象細胞群の組織内での自己組織化過程の実験を行う前に推定することができる。
【００３１】
［実施例１］
（ランダムパターンの画像解析）
まず、パターン形状が制御されていないランダムパターンについて画像解析を行った結果を示す。
図４は、本実施の形態による相関解析処理を行った結果の一例を示す図であり、図２に示したランダムなパターン形状を有する細胞群の画像解析例を示す図である。図５は、同一試料のタイムラプス画像の比較結果として得られた相関度と、異なる試料のタイムラプス画像の比較結果として得られた相関度とを比較したグラフである。
【００３２】
上記の画像相関解析技術を用いて、血管内皮細胞の初期位置と血管内皮細胞ネットワークの形状の相関性を定量的に評価した。試料Ａについて、血管内皮細胞がランダムな場合の血管内皮細胞の初期位置(Ｐ１：Day0)と、同一視野における血管内皮細胞ネットワークの形状(Ｐ２：Day5)の相関度は0.2246（0.24±0.03）であり、正の相関が認められた。同様に、試料Ｂについても、血管内皮細胞がランダムな場合の血管内皮細胞の初期位置(Ｐ１１：Day0)と、血管内皮細胞ネットワークの形状(Ｐ1２：Day5)の相関度は0.1957）であり、正の相関が認められた。すなわち、試料Ａ及び試料Ｂとのそれぞれにおいて、５日後においても、血管内皮細胞のパターン形状の局在は初期のパターン形状の局在に依存していることが明らかになった。
【００３３】
一方、試料Ａの血管内皮細胞の初期位置(Ｐ１：Day0)と試料Ａとは異なる試料Ｂの血管内皮細胞ネットワークの形状(Ｐ１２：Day5)との相関度と、試料Ｂの血管内皮細胞の初期位置(Ｐ１１：Day0)と試料Ａの血管内皮細胞ネットワークの形状(Ｐ２：Day5)との相関度は0.02±0.01であり、同一試料同士の相関度と有意差がある。
【００３４】
図５に示すように、同一の試料と異なる試料とでは、初期と５日後との相関度に大きな差異がみされることがわかる。すなわち、同一試料内のある視野における一定時間前後の画像の相関度は、異なる試料であるか同じ試料であるかを判断する指標として用いることができる。また、２つの試料を比較する際の相関度のしきい値を、図５等の解析結果に基づいて、例えば、0.1から0.2程度に定めることができる。このしきい値を設定することで、同じ試料か異なる試料かを自動的に判定することも可能である。これにより、培養容器内の試料の取り違えを防ぐシステムを構築することが可能となる。
【００３５】
このようにしきい値設定部１−７が、判断部１−５において局在の変化の度合いを判断するための適切なしきい値を、例えば図５を参照して説明したように、同一の試料であるか否かを切り分けるための相関度の値に基づいて設定する。
【００３６】
また、図６に示すように、判断部１−５において、同じ試料の例えば初期位置と比較した相関度を時間に対してプロットすることで、評価対象細胞群の位置と形状との局在の変化の指標として求めることにより評価対象細胞群の組織内での自己組織化過程を数値化することができる。この処理は位置−時間演算部１−４により行うことができる。
【００３７】
また、図６に示すような、位置−時間演算部１−４により求められた相関度の時間依存性に基づいて、同一又は同様な評価対象細胞群のパターン形状の位置と局在の時間変化を推測することが可能である。
【００３８】
以上の解析結果より、画像相関解析を用いることで、ランダムパターンの血管内皮細胞ネットワーク形状工程は、血管内皮細胞の初期位置に依存することが示唆された。
尚、評価対象細胞群は、人工的に作製された細胞組織であっても良い。
【００３９】
［実施例２］
（ラインパターンの画像解析例）
次に、上記のようなランダムパターンとは異なり、制御されたパターン、例えば同じ方向に延在し互いに隣接する複数のラインパターンの評価結果について説明する。
図７は、このようなラインパターンの作製方法の一例を示す図である。
【００４０】
温度応答性培養皿(35 mmf ,UpCellTM, CellSeed Inc., Tokyo, Japan)上に光リソグラフィ技術を用いて、水溶性カンファーキノンを重合開始剤としてポリアクリルアミド（PAAm）を部分的に修飾した。これにより、ストライプ状の温度応答性領域とその他の細胞非接着性領域からなるパターン化温度応答性培養皿を作製した。作製したパターン化温度応答性培養皿はEOGガスにより滅菌を行い培養に用いた。
【００４１】
図７に示すように、パターンデザインについてはストライプパターンにおける細胞接着領域の幅（ライン）：細胞非接着領域の幅（スペース）が、500 μm: 2000 μm (ストライプパターン中の各血管内皮細胞によるライン間隔Ddish= 2500μm），500 μm: 1000 μm (Ddish= 1500μm）, 500 μm: 500 μm (Ddish= 1000μm）の３通りを作製した。
【００４２】
これを、線維芽細胞シートと伸展した状態を保持したまま積層した後に収縮処理を行い、所定の間隔Ｄｔｉｓｓｕｅを有する血管内皮細胞によるラインパターンを含む三次元組織を作製した。
【００４３】
作製した三次元組織中の血管内皮細胞の動きは、タイムラプス蛍光顕微鏡装置(BZ-9000, Keyence Corp., Osaka, Japan)を用い、5日間同一箇所のタイムラプス蛍光画像を取得することで観察した。シート状の細胞組織は、50 ng/mL recombinant human vascular endothelial growth factor (rhVEGF) (R and D Systems, Minneapolis, USA)を添加したEGM-2中で、抗生物質存在下において顕微鏡ステージ装着型インキュベータ (INUG2-KI3, Tokai Hit Co. Ltd., Shizuoka, Japan) を用いて37 ℃,CO2濃度5%環境下で培養した。
【００４４】
その結果を図８〜図１１までに示す。図８、図９は、同じ試料の画像相関解析結果の経時変化を示すものであり、作製したラインパターンの所定の間隔Ｄｔｉｓｓｕｅが703±45μｍ（もしくは、Ddish2500μｍ）である。図１０、図１１は、ライン間隔の異なる試料の画像相関解析結果の経時変化と、相関度のラインパターンの所定の間隔（Ｄｔｉｓｓｕｅ）依存性存性を示す図である。
【００４５】
図８、図９に示すように、三次元組織内のラインパターンの経時変化をみると、積層後５日間にわたって局在を保つことが可能であることが明らかとなった。
【００４６】
先行研究では伸展した状態の細胞からなる三次元組織内において、血管内皮細胞による２本のラインパターンの所定の間隔（Ｄｔｉｓｓｕｅ）が500 μm程度の場合、積層後１日で血管内皮細胞はランダムに遊走していることが確認されている（参考：M. Muraoka et al. / Biomaterials 34 (2013) 696‐703）。
【００４７】
先行研究ではVEGFやFBSを始めとする血管内皮細胞の遊走や増殖を促す因子をほぼ含まない条件で培養されたものであり、血管内皮細胞の局在を保つという点では有利な条件であったと考えられるが、本結果より組織を収縮させるという簡便な手法で血管内皮細胞を活性化させる因子を加えた条件においてもパターン形状の維持が可能であることが示された。この現象は、組織全体の収縮により血管内皮細胞の足場の状態が変化したことで引き起こされたと考えられる。
【００４８】
また、図１０、図１１に示すように、異なるライン間隔の試料を用いて本実施の形態による相関解析を行うと、試料Ａのようにライン間隔が狭い場合には、相関度の時間依存性が大きくなり、すなわち相関度が小さくなり、試料Ｃのようにライン間隔が大きくなると、相関度の時間依存性が小さくなり、すなわち、相関度が大きくなることがわかる。
【００４９】
図１１に示すように、パターン制御を行わない場合の相関度は０.２程度であり、ラインパターン幅が狭い試料Ａはそれと同程度であることがわかる。また、ラインパターン幅を広くしていくと、相関度が大きくなり、ライン制御が有効になっていることがわかる。
【００５０】
尚、評価対象細胞群のパターン形状が制御された血管内皮細胞によるライン状パターンの場合に、しきい値を、０．６から０．８に設定すると良い。
【００５１】
このように、時間依存を求めることで、評価対象細胞群の組織内での自己組織化過程を数値化することができる。
微細構造をデザインした三次元組織の構築においても、上記のような製造技術が有効であることを、客観的に示すことができる。
【００５２】
尚、本実施の形態においては、評価対象細胞として、血管内皮細胞を例にして説明したが、その他、神経細胞、筋芽細胞、肝細胞、その他各種多能性幹細胞などを対象としても良い。
【００５３】
本発明は、創薬研究のための組織モデル作製および再生医療に利用可能な三次元組織作製において有用な技術となると考えられる。これらの組織においては、血管網を始めとする組織内の微細構造の局在を制御することが生体機能の模倣や効率的な組織構築のために重要であると期待されている。本発明は、三次元組織内での細胞の局在制御や遊走能解析研究において汎用性の高い定量評価用ツールとなりうる。
【００５４】
上記の実施の形態において、添付図面に図示されている構成等については、これらに限定されるものではなく、本発明の効果を発揮する範囲内で適宜変更することが可能である。その他、本発明の目的の範囲を逸脱しない限りにおいて適宜変更して実施することが可能である。
【００５５】
また、本発明の各構成要素は、任意に取捨選択することができ、取捨選択した構成を具備する発明も本発明に含まれるものである。
【００５６】
また、本実施の形態で説明した機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。尚、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
【産業上の利用可能性】
【００５７】
本発明は、細胞群のパターン形状の局在を評価する装置として利用可能である。
【符号の説明】
【００５８】
１…画像解析装置、１−１…画像取得部、１−２…パターン抽出部、１−２−１…タイムラプス解析部、１−２−２…画像データ取得部、１−２−３…二値化部、１−２−４…ノイズ除去部、１−３…相関解析部、１−４…位置−時間演算部、１−５…判断部、１−６…出力部。
【図１】
0006661165.tif 000003
【図２】
0006661165.tif 000004
【図３Ａ】
0006661165.tif 000005
【図３Ｂ】
0006661165.tif 000006
【図４】
0006661165.tif 000007
【図５】
0006661165.tif 000008
【図６】
0006661165.tif 000009
【図７】
0006661165.tif 000010
【図８】
0006661165.tif 000011
【図９】
0006661165.tif 000012
【図１０】
0006661165.tif 000013
【図１１】
0006661165.tif 000014