(57)【要約】  本発明の目的は、複数の種類の細胞集団が混在する場合に、細胞集団間の相互作用を効率的に解析するための細胞情報処理方法を提供することである。本発明の細胞情報処理方法は、異なる複数種類の細胞を含む対象細胞群を播種した所定の容器を少なくとも２以上用意し、対象細胞群に対し、所定の細胞刺激を与えて培養し、２以上の時点で、１つの容器内にある全細胞を回収する細胞回収ステップと、各時点で回収した全細胞をシングルセル化するシングルセル化ステップと、各時点のシングルセル化された各細胞からシングルセルデータを取得するシングルセルデータ取得ステップと、各時点のシングルセルデータに基づいて、各時点で回収された全細胞を共通の第１の細胞特徴を有する複数の細胞集団にグループ分けして二次元又は三次元平面上にプロットし、且つ、第２の細胞特徴に基づいて、グループ分けした各細胞集団の細胞種を同定する処理を行い、各時点におけるクラスタリング結果を取得するグルーピングステップと、各時点におけるクラスタリング結果を比較することにより、同じ細胞種の細胞集団の経時的な変化を検出する細胞変化検出ステップと、検出の結果に基づいて、同じ細胞種の細胞集団の経時的な変化の作用機序を解析する作用機序解析ステップと、作用機序解析の結果に基づいて、異なる細胞種の細胞集団の間で起きる相互作用を解析する相互作用解析ステップと、を含む、細胞情報処理方法。
2020071500.tif 000002
【特許請求の範囲】
【請求項１】
異なる複数種類の細胞を含む対象細胞群を播種した所定の容器を少なくとも２以上用意し、前記対象細胞群に対し、所定の細胞刺激を与えて培養し、２以上の時点で、１つの容器内にある全細胞を回収する細胞回収ステップと、
各時点で回収した全細胞をシングルセル化するシングルセル化ステップと、
各時点のシングルセル化された各細胞からシングルセルデータを取得するシングルセルデータ取得ステップと、
各時点の前記シングルセルデータに基づいて、各時点で回収された全細胞を共通の第１の細胞特徴を有する複数の細胞集団にグループ分けして二次元平面又は三次元空間上にプロットし、且つ、第２の細胞特徴に基づいて、グループ分けした各細胞集団の細胞種を同定する処理を行い、各時点におけるクラスタリング結果を取得するグルーピングステップと、
前記各時点におけるクラスタリング結果を比較することにより、同じ細胞種の前記細胞集団の経時的な変化を検出する細胞変化検出ステップと、
前記検出の結果に基づいて、前記同じ細胞種の前記細胞集団の経時的な変化の作用機序を解析する作用機序解析ステップと、
前記作用機序解析の結果に基づいて、異なる細胞種の前記細胞集団の間で起きる相互作用を解析する相互作用解析ステップと
を含む、細胞情報処理方法。
【請求項２】
前記相互作用解析ステップの前記細胞集団間で起きる相互作用には、前記同じ細胞種の前記細胞集団の間で起きる相互作用も含む請求項１に記載の細胞情報処理方法。
【請求項３】
前記細胞変化検出ステップは、前記各時点におけるクラスタリング結果を比較することにより、前記同じ細胞種の前記細胞集団を構成する細胞数の経時的な変化、及び経時的な前記第１の細胞特徴の変化を検出する請求項１または２に記載の細胞情報処理方法。
【請求項４】
前記細胞変化検出ステップは、さらに、前記各時点におけるクラスタリング結果において、前記同じ細胞種の前記細胞集団を構成する細胞数の経時的変化、及び経時的な前記第１の細胞特徴の変化を検出する請求項３に記載の細胞情報処理方法。
【請求項５】
前記細胞回収ステップは、さらに、前記所定の細胞刺激を与えずに培養する前記複数種類の細胞を含む対象細胞群を用意し、１以上の時点で、１つの前記所定の容器内にある全細胞を回収し、
前記細胞変化検出ステップは、前記各時点におけるクラスタリング結果を比較することにより、前記細胞集団の前記所定の細胞刺激の有無による、前記同じ細胞腫の細胞集団の経時的な変化を評価する請求項１〜４のいずれか１項に記載の細胞情報処理方法。
【請求項６】
前記細胞刺激は、化学的刺激および物理的刺激からなる群から選択される少なくとも１種である、請求項１〜５のいずれか１項に記載の細胞情報処理方法。
【請求項７】
前記化学的刺激は、細胞に対して生物学的な反応を誘起する薬剤の添加によるものである、請求項６に記載の細胞情報処理方法。
【請求項８】
前記作用機序は、前記細胞刺激により細胞内外で誘起された生物学的な現象を発揮するための特異的な生化学的反応または相互作用である、請求項１〜７のいずれか１項に記載の細胞情報処理方法。
【請求項９】
前記シングルセルデータは、単一細胞の機能や性質、状態を示す生体物質の情報であり、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾情報、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度からなる群から選択される少なくとも１つである、請求項１〜８のいずれか１項に記載の細胞情報処理方法。
【請求項１０】
前記シングルセルデータは、遺伝子発現量及び遺伝子のＤＮＡ配列である請求項９に記載の細胞情報処理方法。
【請求項１１】
前記グルーピングステップにおいて、前記第１の細胞特徴とは、前記シングルセルデータに含まれるｎ次元の細胞特徴を２次元または３次元に次元削減したものである請求項１〜１０のいずれか１項に記載の細胞情報処理方法。
【請求項１２】
前記次元削減の方法は、主成分分析（PCA）、カーネルあり主成分分析（Kernel-PCA）、多次元尺度構成法（MDS）、t-SNE、及び畳込みニューラルネットワーク（CNN）からなる群から選択される少なくとも１つである、請求項１１に記載の細胞情報処理方法。
【請求項１３】
前記グルーピングステップにおいて、前記第１の細胞特徴とは、前記遺伝子発現量について主成分分析を行い、２次元又は３次元に次元削減して獲得されるものである請求項１０〜１２のいずれか１項に記載の細胞情報処理方法。
【請求項１４】
前記グルーピングステップにおいて、前記第２の細胞特徴とは、細胞の機能または細胞の状態から各細胞種を同定することが可能な少なくとも１つの細胞情報である、請求項１〜１３のいずれか１項に記載の細胞情報処理方法。
【請求項１５】
前記細胞情報とは、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾情報、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度からなる群から選択される少なくとも１つである請求項１４に記載の細胞情報処理方法。
【請求項１６】
前記細胞の機能とは、細胞の増殖、修復、代謝、および細胞間の情報交換から選択される少なくとも１つである、請求項１４または１５に記載の細胞情報処理方法。
【請求項１７】
前記細胞の状態とは、遺伝子の発現状況、タンパク質の発現状況、および酵素活性から選択される少なくとも１つである、請求項１４または１５に記載の細胞情報処理方法。
【請求項１８】
前記相互作用解析ステップにおいて、前記細胞集団の間で起きる相互作用とは、前記異なる細胞種の細胞集団の間で起きる生物学的または物理的な作用、または前記作用による細胞数、前記第１の細胞特徴の変化である、請求項１〜１７のいずれか１項に記載の細胞情報処理方法。
【請求項１９】
前記シングルセルデータから複数の生体物質の量又は状態を表す値を、前記生体物質ごとにそれぞれ複数の時点において取得した時系列データを予め作成し、前記生体物質ごとの時系列データの時間変化と、各生体物質の生物学的機能の類似性に基づいて、前記シングルセルデータを取得した細胞を、前記共通の第１の細胞特徴を有する細胞集団にグループ分けし、
前記相互作用解析ステップにおいて、前記複数の時点の各々について、複数の細胞集団の各々に含まれる１つ以上の第１の細胞特徴から、前記細胞集団の状態を表す値を生成し、生成された、複数時点の、複数の細胞集団の状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データからなるデータセットから、細胞集団間の状態の依存関係を推定する、請求項１〜１８のいずれか１項に記載の細胞情報処理方法。
【請求項２０】
前記生物学的機能の類似性は、共通の遺伝子オントロジーを有すること、共通のカノニカルパスウェイに属すること、共通の上流因子を有すること、共通の表現系に関わること、および、共通の疾患に関わることからなる群から選択される少なくとも１つに基づいて評価されるものである、請求項１９に記載の細胞情報処理方法。
【請求項２１】
前記細胞回収ステップ及び前記シングルセル化ステップにおいて、全細胞を回収し、シングルセル化する方法が、手動、フローサイトメトリー、磁気分離、レーザーキャプチャーマイクロダイセクション、マイクロ流路、マイクロドロップレット、ナノウェル、マイクロピペット微細針吸引、レーザーピンセット、標識アレイ、表面プラズモンレスポンス、およびナノバイオデバイスからなる群から選択される少なくとも１つを用いる方法である、請求項１〜２０のいずれか１項に記載の細胞情報処理方法。
【請求項２２】
前記全細胞を回収し、シングルセル化する方法において、細胞の標識として蛍光標識、ラジオアイソトープ標識、抗体標識、および磁気標識からなる群から選択される少なくとも１つを用いる、請求項２１に記載の細胞情報処理方法。
【請求項２３】
前記複数種類の細胞を含む対象細胞群は、生体組織サンプル、血液サンプル、培養サンプル、および環境サンプルからなる群から選択される少なくとも１つから得られた細胞である、請求項１〜２２のいずれか１項に記載の細胞情報処理方法。
【請求項２４】
前記複数種類の細胞は、動物細胞、植物細胞、真菌細胞および細菌細胞からなる群から選択される少なくとも１つである、請求項２３に記載の細胞情報処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、細胞情報処理方法に関する。特に、異なる複数の細胞種が混在する状況下で培養された細胞集団に対し、薬剤による細胞刺激を与えた場合に生じる異なる細胞種間の相互作用、及び、その作用機序を解析するための方法に関する。
【背景技術】
【０００２】
生体内には多数の細胞種類が存在し、互いに依存しながら共存し、相互刺激により変化する。また、変化する際には、１細胞内で多数の遺伝子が依存関係を持ちながら変化し、遺伝子はこの依存関係を通じて生体内で機能を発揮している。そのため、創薬分野や、再生医療分野において、生命現象および生物現象を理解するために、細胞間または細胞集団間の相互作用および細胞内での分子挙動のメカニズムを解析することは重要である。例えば、複数種類の細胞を混合培養した際に、どのように細胞間で相互に刺激を与えあい、また、その刺激により細胞内での遺伝子発現量または遺伝子の状態が、どのような依存関係をもって時間変化をするかを解析することによって、生命現象の機構を解明するための手がかりを得ることができる。【０００３】
上記メカニズムを解析に利用できる方法として、特許文献１には、「（１）遺伝子の経時的な発現量変化を示すデータから、各データ間の類似度を反映した特徴量を算出するステップと、（２）算出された特徴量から、全ての遺伝子間の組合せについて類似度行列Ｍの固有ベクトルを算出するステップと、（３）類似度行列Ｍを、固有ベクトルの固有値を維持したまま、ブール行列Ｎに変換するステップと、（４）ブール行列Ｎに基づいて各データをクラスタリングするステップと、を少なくとも行う遺伝子クラスタリング方法。」、即ち、経時的な発現量変化の類似性に基づいて、複数の遺伝子をグループ化する方法が提案されている（請求項９）。【０００４】
また、特許文献２には、遺伝子発現プロファイルをクラスター化する方法であって、遺伝子オントロジー（ＧＯ）ツリーから１または２以上のＧＯタームを選択するステップ；遺伝子発現データセットを受け取るステップ；遺伝子発現データセットをＧＯタームに従ってグループに分類するステップ；第１に、遺伝子発現データの類似性に基づいて各群に属する遺伝子発現データをクラスタリングするステップ；および、第２に、第１のクラスタリングの結果をシードとして使用して遺伝子発現データセットをクラスタリングするステップを含む方法」、即ち、発現データの類似性と、それらが関与する生物学的機能の類似性とに基づいて、複数の遺伝子をグループ化する方法が提案されている（クレーム１）。【０００５】
また、従来の薬効及び薬理を解析する方法としては、図８に示すように、各ウェル３６に個々に配置された同種の細胞３０、３２及び３４のみからなる細胞集団に薬剤をそれぞれ与えて培養し、その薬効を各細胞集団の平均化データを用いて解析する方法、又は、臓器から抽出された複数の細胞種を含む細胞集団に薬剤を与えて培養し、薬効を細胞集団の平均化データを用いて解析する方法が用いられてきた。【０００６】
しかし、同種の細胞のみからなる細胞集団でも、個々の細胞の経時的な動的挙動（例えば、細胞の活性化、細胞の阻害、細胞の相互作用、タンパク質発現、タンパク質分泌、細胞増殖、細胞形態の変化、細胞死等）に不均一性があるため、従来の方法では、十分に精度の高い解析を行なえていなかった。そこで、近年では、さらに解析精度を高めるために、同種の細胞のみからなる細胞集団を構成する細胞をシングルセル化し、個々の細胞（単一細胞）レベルでの解析（シングルセル解析）が行われている。このようなシングルセル解析により、細胞集団の平均化データの解析ではなく、個々の細胞を解析することができるようになったため、平均化データによる解析では検出できなかった細胞や、小さな動的挙動の変化を検出することができるようになった。このようなシングルセル解析は、例えば、癌細胞のシングルセルトランスクリプトーム解析により、癌組織中に存在する細胞亜種を分類することに用いられている。また、シングルセル解析を用いて細胞活性を評価する方法として、特許文献３に記載される方法が提案されている。【０００７】
特許文献３には、細胞活性を評価する方法であって、（ａ）領域上に異なる複数種類の細胞を含む細胞集団を設置する工程、即ち、ナノウェルアレイプレート上の各ナノウェルに、異なる複数種類の細胞を含む細胞集団を設置する工程；（ｂ）該細胞集団の動的な挙動を、時間の関数としてアッセイする工程、即ち、単一細胞レベルの動的な挙動を顕微鏡、蛍光顕微鏡等で可視化して（つまり、画像解析により）経時的に測定する工程；（ｃ）該細胞集団から、動的な挙動に基づいて、少なくとも１つの解析対象の細胞を同定する工程；（ｄ）同定された少なくとも１つの解析対象の細胞の分子プロファイルを特徴付ける工程、即ち、同定された少なくとも１つの解析対象細胞について単一細胞レベルで、質量分析、遺伝子分析、タンパク質分析等を行い、転写活性、トランスクリプトームのプロファイル、遺伝子発現活性、ゲノムプロファイル、タンパク質発現活性、プロテオームのプロファイル、タンパク質の相互作用活性、細胞受容体の発現活性、脂質プロファイル、脂質活性、炭水化物プロファイル、微小胞活性、グルコース活性、代謝プロファイル、細胞受容体の発現活性等を取得する工程；および（ｅ）工程（ｂ）および（ｄ）から得られた情報を相関させる工程を含む、上記方法が提案されている（請求項１）。
即ち、特許文献１に記載の方法は、異なる複数種類の細胞を含む細胞集団において、個々の細胞種の細胞活性を評価するものである。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１０−１５７２１４号公報
【特許文献２】米国特許出願公開第２００９／０１１２４８０号明細書
【特許文献３】特表２０１８−５１４１９５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかし、特許文献１及び２に記載の方法はいずれも、細胞内ネットワークを調べる方法であり、細胞間及び細胞種間の相互作用を知ることはできない。【００１０】
また、特許文献３に記載された方法も、各細胞活性を評価するだけの方法であるため、細胞間及び細胞種間の細胞種間の相互作用を知ることはできない。
また、特許文献３に記載された方法は、人間が、画像解析によって少なくとも１００〜２００個の細胞からなる細胞集団の動的な挙動を観察し、複数種類の細胞及び動的な挙動が異なる細胞の中から解析対象の単一細胞の同定及び選択を行い、人間によって選び出された単一細胞の動的な挙動の作用機序および細胞間相互作用を解析するものであるため、つまり、各細胞種の形態や、経時的な各細胞の動的な変化を人間が恣意的に判断して解析を行う細胞を同定及び選択するため、精度の高い解析結果を得ることができないという問題があった。
また、特許文献１の一実施態様として記載されているように、解析対象の細胞集団に免疫細胞が含まれる場合、具体的には、血液サンプル中に免疫細胞が含まれる場合、血液サンプルに薬剤を与えてから約２４時間で細胞が死滅してしまうという事情がある。また、特許文献１に記載の解析方法は、薬剤を与えた細胞集団を経時的に観察し、解析を行うものであるため（つまり、１つのサンプルの観察を解析が終わるまで続けるものであり、また、さらに、解析対象の細胞の選択及び培養を繰り返しながら解析する方法であるため）、特許文献３に記載の方法を採用する場合、薬剤を細胞集団に与えてから少なくとも２４時間以内に全ての解析を終わらせなければならないという非常に厳しい時間的制約がある。【００１１】
また、特許文献３に記載の方法は、細胞同定や対象細胞の選択操作に手間やコストがかかり、また、その解析による評価の算出方法及びその結果も複雑であるため、効率的で、高精度、且つ迅速な評価及び解析を必要とする産業的に用いる方法としては実用的ではないという問題がある。
さらに、異なる複数の種類の細胞を含む細胞集団に薬剤をそれぞれ与えて培養し、シングルセル解析技術を用いて、異なる複数の細胞種を含む細胞集団において実施する各細胞種の同定方法、各細胞種の薬効及び薬理を効率的に解析する方法は現時点で報告されていない。そのため、異なる複数の細胞種が混在する状況下で培養された細胞集団に対し、薬剤による細胞刺激を与えた場合に生じる異なる細胞種間の相互作用、及び、その作用機序をシングルセル解析に基づいて効率的に解析する方法も現時点で報告されていない。【００１２】
また、図８に示すような薬効及び薬理スクリーニング方法を用いた薬剤の探索の初期段階では、薬剤（細胞刺激）による細胞死の無しの判断のみに基づいて、異なる細胞との相互作用解析や、薬効及び薬理のメカニズムの解析が進められるため、膨大な数の候補物質から薬剤となり得る物質を抽出してくる解析初期段階においては、例えば、細胞死が検出されたとしても、それがどのような要因により、細胞死が引き起こされたのか不明なまま、薬剤候補となり得なかったりする可能性があった。【００１３】
本発明は、従来技術の有する上記問題点を鑑みて、異なる複数の種類の細胞集団が混在する場合に、細胞集団間の相互作用を効率的に解析するための細胞情報処理方法を提供することを目的とする。
より具体的に言えば、異なる複数の細胞種が混在する状況下で培養された細胞集団に対し、薬剤または細胞刺激を与えた場合に生じる異なる細胞種間の相互作用、及び、その作用機序をより簡単な操作で、高精度に且つ迅速に評価及び解析することが効率的にでき、産業的な規模でも利用できる方法を提供することを課題とする。
【課題を解決するための手段】
【００１４】
本発明の細胞情報処理方法は、異なる複数種類の細胞を含む対象細胞群を播種した所定の容器を少なくとも２以上用意し、上記対象細胞群に対し、所定の細胞刺激を与えて培養し、２以上の時点で、１つの容器内にある全細胞を回収する細胞回収ステップと、各時点で回収した全細胞をシングルセル化するシングルセル化ステップと、各時点のシングルセル化された各細胞からシングルセルデータを取得するシングルセルデータ取得ステップと、各時点の上記シングルセルデータに基づいて、各時点で回収された全細胞を共通の第１の細胞特徴を有する複数の細胞集団にグループ分けして二次元平面又は三次元空間上にプロットし、且つ、第２の細胞特徴に基づいて、グループ分けした各細胞集団の細胞種を同定する処理を行い、各時点におけるクラスタリング結果を取得するグルーピングステップと、各時点におけるクラスタリング結果を比較することにより、同じ細胞種の上記細胞集団の経時的な変化を検出する細胞変化検出ステップと、検出の結果に基づいて、同じ細胞種の上記細胞集団の経時的な変化の作用機序を解析する作用機序解析ステップと、作用機序解析の結果に基づいて、異なる細胞種の上記細胞集団の間で起きる相互作用を解析する相互作用解析ステップと、を含む。【００１５】
上記相互作用解析ステップの上記細胞集団間で起きる相互作用には、同じ細胞種の細胞集団の間で起きる相互作用も含むことが好ましい。
上記細胞変化検出ステップは、上記各時点におけるクラスタリング結果を比較することにより、同じ細胞種の上記細胞集団を構成する細胞数の経時的な変化、及び経時的な上記第１の細胞特徴の変化を検出することもできる。
上記細胞変化検出ステップは、さらに、上記各時点におけるクラスタリング結果において、同じ細胞種の上記細胞集団を構成する細胞数の経時的変化、及び経時的な上記第１の細胞特徴の変化を検出してもよい。【００１６】
上記細胞回収ステップは、さらに、上記所定の細胞刺激を与えずに培養する上記複数種類の細胞を含む対象細胞群を用意し、１以上の時点で、１つの上記所定の容器内にある全細胞を回収し、上記細胞変化検出ステップは、上記各時点におけるクラスタリング結果を比較することにより、上記細胞集団の上記所定の細胞刺激の有無による、上記同じ細胞腫の細胞集団の経時的な変化を評価することもできる。
上記細胞刺激は、化学的刺激および物理的刺激からなる群から選択される少なくとも１種であるのが好ましい。
上記化学的刺激は、細胞に対して生物学的な反応を誘起する薬剤の添加によるものであることが好ましい。
上記作用機序は、上記細胞刺激により細胞内外で誘起された生物学的な現象を発揮するための特異的な生化学的反応または相互作用であることが好ましい。【００１７】
上記シングルセルデータは、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾情報、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度からなる群から選択される少なくとも１つであることが好ましい。
上記シングルセルデータは、遺伝子発現量及び遺伝子のＤＮＡ配列であることが好ましい。【００１８】
上記グルーピングステップにおいて、上記第１の細胞特徴とは、上記シングルセルデータに含まれるｎ次元の細胞特徴を２次元または３次元に次元削減したものであることが好ましい。
上記次元削減の方法は、主成分分析（PCA）、カーネルあり主成分分析（Kernel-PCA）、多次元尺度構成法（MDS）、t-SNE、及び畳込みニューラルネットワーク（CNN）からなる群から選択される少なくとも１つであることが好ましい。
上記グルーピングステップにおいて、上記第１の細胞特徴とは、上記遺伝子発現量について主成分分析を行い、２次元又は３次元に次元削減して獲得されるものであることが好ましい。
上記グルーピングステップにおいて、上記第２の細胞特徴とは、細胞の機能または細胞の状態から各細胞種を同定することが可能な少なくとも１つの細胞情報であることが好ましい。【００１９】
上記細胞情報とは、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾情報、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度からなる群から選択される少なくとも１つであることが好ましい。
上記細胞の機能とは、細胞の増殖、修復、代謝、および細胞間の情報交換から選択される少なくとも１つであることが好ましい。
上記細胞の状態とは、遺伝子の発現状況、タンパク質の発現状況、および酵素活性から選択される少なくとも１つであることが好ましい。【００２０】
上記相互作用解析ステップにおいて、上記細胞集団の間で起きる相互作用とは、異なる細胞種の細胞集団の間で起きる生物学的または物理的な作用、または上記作用による細胞数、上記第１の細胞特徴の変化であることが好ましい。
上記シングルセルデータから複数の生体物質の量又は状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データを予め作成し、上記生体物質ごとの時系列データの時間変化と、各生体物質の生物学的機能の類似性に基づいて、シングルセルデータを取得した細胞を、共通の第１の細胞特徴を有する細胞集団にグループ分けし、上記相互作用解析ステップにおいて、上記複数の時点の各々について、複数の細胞集団の各々に含まれる１つ以上の第１の細胞特徴から、上記細胞集団の状態を表す値を生成し、生成された、複数時点の、複数の細胞集団の状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データからなるデータセットから、細胞集団間の状態の依存関係を推定することが好ましい。
上記生物学的機能の類似性は、共通の遺伝子オントロジーを有すること、共通のカノニカルパスウェイに属すること、共通の上流因子を有すること、共通の表現系に関わること、および、共通の疾患に関わることからなる群から選択される少なくとも１つに基づいて評価されるものであることが好ましい。【００２１】
上記細胞回収ステップ及び上記シングルセル化ステップにおいて、全細胞を回収し、シングルセル化する方法が、手動、フローサイトメトリー、磁気分離、レーザーキャプチャーマイクロダイセクション、マイクロ流路、マイクロドロップレット、ナノウェル、マイクロピペット微細針吸引、レーザーピンセット、標識アレイ、表面プラズモンレスポンス、およびナノバイオデバイスからなる群から選択される少なくとも１つを用いる方法であることが好ましい。
上記全細胞を回収し、シングルセル化する方法において、細胞の標識として蛍光標識、ラジオアイソトープ標識、抗体標識、および磁気標識からなる群から選択される少なくとも１つを用いることが好ましい。
上記複数種類の細胞を含む対象細胞群は、生体組織サンプル、血液サンプル、培養サンプル、および環境サンプルからなる群から選択される少なくとも１つから得られた細胞であることが好ましい。
上記複数種類の細胞は、動物細胞、植物細胞、真菌細胞および細菌細胞からなる群から選択される少なくとも１つであることが好ましい。
【発明の効果】
【００２２】
本発明によれば、複数の種類の細胞集団が混在する場合に、異なる細胞集団間の相互作用を効率的に解析するための細胞情報処理方法を提供できる。
より具体的に言えば、異なる複数の細胞種が存在する細胞集団において実施する各細胞種の同定、若しくは薬効及び薬理スクリーニングをより簡単な操作で、高精度に且つ迅速に評価及び解析することができ、産業的な規模でも利用できる方法を提供することができる。そのため、同じ細胞種からなる細胞集団間の変化だけでなく、異なる細胞種からなる細胞集団の変化（細胞集団間の相互作用）を効率的、且つ、同時に確認することができる。また、その結果、より広い細胞間及び細胞種間の相互作用（相互ネットワーク）も精度高く、容易に解析することができる。
また、経時的な細胞の動的な挙動等を人間が恣意的に判断することがないため、各細胞種の同定、各種細胞に対する薬剤または細胞刺激の作用機序を、精度高く、解析することができる。
また、画像解析等の手間をかけることなく、対象細胞を同定することができる。
また、免疫細胞に対する薬効スクリーニングを行う場合であっても、特に解析時間に制約を設けることなく解析を行うことができる。
また、同じ細胞種ごとに培養したサンプルそれぞれに薬剤または細胞刺激を与えて解析する必要なく、複数種類の細胞を共培養したサンプルに薬剤または細胞刺激を与えて解析することができるため、解析のために、多量のサンプルを用意したり、解析したりする手間や費用の負担を減らすことができる。
また、解析対象の細胞が、細胞の生育および機能維持のために１種類の細胞だけで培養することができないものであっても、本発明によれば、複数種類の細胞を共培養したサンプルから解析対象の細胞を同定することができる。また、解析対象以外の細胞との共培養による作用機序、例えば、解析対象細胞と免疫細胞との共培養が必要なサンプル薬剤又は細胞刺激の作用機序も適切に解析することができる。
【図面の簡単な説明】
【００２３】
【図１】図１は、本発明の細胞情報処理方法を説明するフローチャートである。
【図２Ａ】図２Ａは、細胞刺激直後（時点０）の細胞を示す図である。
【図２Ｂ】図２Ｂは、細胞刺激後、１時間後（時点１）の細胞を示す図である。
【図２Ｃ】図２Ｃは、細胞刺激後６時間後（時点２）の細胞を示す図である。
【図３Ａ】図３Ａは、図２Ａ（時点０）に係るクラスタリングの結果を示す図である。
【図３Ｂ】図３Ｂは、図２Ｂ（時点１）に係るクラスタリングの結果を示す図である。
【図３Ｃ】図３Ｃは、図２Ｃ（時点２）に係るクラスタリングの結果を示す図である。
【図４】図４は、図３Ｂに基づいて、細胞集団の経時的な変化を示す図である。
【図５】図５は、図３Ｃに基づいて、細胞集団の経時的な変化を示す図である。
【図６】図６は、神経細胞及びグリア細胞の各成分量を経時的に測定した結果を示す図である。
【図７】図７は、図３Ｃ及び図６に基づいて、細胞集団の経時的な変化を示す図である。
【図８】図８は、従来の薬剤スクリーニング方法を説明するための図である。
【発明を実施するための形態】
【００２４】
［細胞情報処理方法］
以下では、本発明の細胞情報処理方法について詳細に説明する。
実施形態１
実施形態１は、神経細胞とグリア細胞間の相互作用の解析、即ち、神経細胞及びグリア細胞からなる対象細胞群に対し薬剤を添加し、培養した際に生じる神経細胞とグリア細胞の相互ネットワークの形成過程（経時的変化）の解析を事例として説明する。
図１は、本発明の実施形態１に係る細胞情報処理方法を示すフローチャートである。
【００２５】
＜細胞回収ステップ（Ｓ１）＞
まず、本ステップにおいては、図２Ａに示すように、異なる複数種類の細胞群（グリア細胞２及び神経細胞４からなる細胞群）を播種した所定の容器１を少なくとも２以上用意する。ここで、所定の容器１に播種された細胞群を対象細胞群という。次に、準備した容器１に播種された対象細胞群に対しては、薬剤添加による細胞刺激を与えて培養し、２以上の時点で（例えば、薬剤添加直後、薬剤添加後、１時間、６時間、１２時間、２４時間、４８時間、７２時間・・・）１つの容器内にある全細胞を回収する。
【００２６】
図２Ａは、細胞刺激直後（時点０）の細胞を示し、図２Ｂは、細胞刺激を与えた後、１時間経過後（時点１）の細胞を示し、図２Ｃは、細胞刺激を与えた後、６時間後（時点２）の細胞を示す。本ステップにおいては、各時点において、容器１内の全細胞を回収する。
なお、時点０は、アストロサイト２及び神経細胞４の形状が変化する前の状態であり、時点１では、薬剤添加による細胞刺激により、アストロサイト２の形状が変化した反応性アストロサイト１０に変化した状態であり、時点２では、変化した反応性アストロサイト１０と神経細胞４が接触したことにより、神経細胞４の形状が変化した神経細胞１２を含む状態を示す。
【００２７】
《対象細胞群》
ここで、対象細胞群は、所定の容器１に播種された細胞群であり、且つ、異なる複数種類の細胞を含む細胞の群を意味する。「異なる複数種類の細胞を含む細胞の群」をより具体的に言えば、例えば、ヒトｉＰＳ細胞及びマウス胎児由来線維芽細胞から構成される細胞群のように、細胞種が異なる複数の細胞から構成される細胞群を意味し、同種由来の細胞で構成される細胞群であっても、異種由来の細胞で構成される細胞群であってもよい。
本実施形態において、対象細胞群を構成する細胞をグリア細胞及び神経細胞として説明しているが、対象細胞群に含まれる細胞の種類は、特に限定されない。
【００２８】
「複数種類の細胞を含む対象細胞群」としては、生体組織サンプル、血液サンプル、培養サンプル、及び環境サンプルが例示される。
上記生体組織サンプルとしては、マウスの脳組織およびヒトの切除腫瘍組織が例示される。
上記血液サンプルとしては、ヒトの採血試料が例示される。
上記培養サンプルとしては、ヒトｉＰＳ細胞とマウス胎児由来線維芽細胞との共培養サンプルが例示される。
上記環境サンプルとしては、土壌サンプルおよび海底熱水噴出孔から採取した水サンプルが例示される。
【００２９】
また、「複数の細胞種を含む対象細胞群」に含まれる「細胞」として、例えば、動物細胞、植物細胞、真菌細胞、および細菌細胞が挙げられる。
【００３０】
上記動物細胞としては、脊椎動物、脊索動物（脊椎動物を除く）または昆虫の細胞が例示される。
上記脊椎動物としては、ヒト、チンパンジー、アカゲザル、イヌ、ブタ、マウス、ラット、チャイニーズハムスター、およびモルモットのような哺乳類、アフリカツメガエル、ゼブラフィッシュ、メダカ、およびトラフグが例示される。
上記哺乳類の細胞（哺乳類細胞）には、腫瘍細胞、肝細胞、繊維芽細胞、幹細胞、及び免疫細胞が含まれるが、これらに限定されない。
上記脊索動物（脊椎動物を除く）としては、ホヤが例示される。
上記昆虫としては、ショウジョウバエ、カイコ、タバコスズメガ、およびミツバチが例示される。
【００３１】
上記植物細胞としては、被子植物の細胞が例示される。
上記被子植物としては、シロイヌナズナ、イネ、コムギ、ミナトカモジグサ、ミヤコグサ、およびタバコが例示される。
【００３２】
上記真菌細胞としては、カビまたは酵母の細胞が例示される。
上記カビとしては、アカパンカビ（Neurospora crassa）、コウジカビ（Aspergillus oryzae）、アスペルギルス・フミガータス（Aspergillus fumigatus）、アスペルギルス・ニジュランス（Aspergillus nidulans）、リゾプス・オリゼ（Rhizopus oryzae）、およびムコール・シルシネロイデス（Mucor circinelloides）が例示される。
上記酵母としては、出芽酵母（Saccharomyces cerevisiae）、分裂酵母（Schizosaccharomyces pombe）、カンジダ・アルビカンス（Candida albicans）、クリプトコッカス・ネオフォルマンス（Cryptococcus neoformans）、およびトリコスポロン・オボイデス（Trichosporon ovoides）が例示される。
【００３３】
上記細菌細胞としては、例えば、エシェリヒア・コリ、サルモネラ・エンテリカ、クロストリジウム・デフィシル、またはバチルス・アンスラシスの細胞が挙げられる。
【００３４】
細胞刺激とは、化学的刺激（化学物質等）および物理的刺激（光、熱、または圧力等）からなる群から選択される少なくとも１種であれば、特に限定されない。
【００３５】
上記化学的刺激の例としては、細胞に対して生物学的な反応を誘起する薬剤の添加によるものが挙げられる。なお、薬剤は生物試料へ添加することにより細胞に対して生物学的な反応を誘起するものであってもよい。
上記生物学的な反応の具体例としては、増殖、細胞死、分化、抗原抗体反応、増殖因子の分泌等が挙げられる。
上述の薬剤とは、身体の構造及び機能に測定可能な効果を有するよう意図される薬剤であれば、特に限定されない。このような薬剤として、抗がん剤等の医薬品、成長因子、サイトカインおよび低分子薬などが挙げられ、成長因子の具体例としては、上皮成長因子（EGF: Epidermal Growth Factor）が挙げられ、サイトカインの具体例としては、腫瘍壊死因子（TNF-α）、インターロイキン１β（IL-1β）、インスリン、グルカゴン様ペプチド−１（GLP-1）、イマチニブ（imatinib）、アセトアミノフェン、アダリムマブ（Adalimumab）、およびニボルマブ（Nivolumab）が挙げられる。
【００３６】
なお、対象細胞群を収容する容器１は、対象細胞群を収容し培養できる細胞培養容器であればよく、特に限定されない。また、使用する培養液も、細胞や解析手法に応じて、適宜好ましいものを用いることができる。
各容器には、複数の細胞種を所定の比率で構成した対象細胞群（例えば、Ａ細胞、Ｂ細胞及びＣ細胞からなり、各細胞の数が、Ａ細胞：Ｂ細胞：Ｃ細胞＝２：１：１で構成される対象細胞群）が播種される。複数の細胞種を含む対象細胞群は、所定の時間培養された後、薬剤が添加され細胞刺激が与えられる。
【００３７】
＜シングルセル化ステップ（Ｓ２）＞
次いで、本ステップにおいて、細胞回収ステップ（Ｓ１）で回収された対象細胞群を、単一細胞化（シングルセル化）する。
【００３８】
対象細胞群をシングルセル化し、単一細胞を回収する方法及び器具は、特に限定されず、公知の方法や器具を使用することができる。例えば、公知の方法としては、手動、フローサイトメトリー、磁気分離、レーザーキャプチャーマイクロダイセクション、マイクロドロップレット法、マイクロピペット微細針吸引法、及び表面プラズモン共鳴法が挙げられ、公知の器具としては、例えば、マイクロ流路、ナノウェル、レーザーピンセット、標識アレイ、およびナノバイオデバイスが挙げられる。
これら公知の方法の中でも、マイクロドロップレット法、マイクロ流路、ナノウェル、フローサイトメトリーを使用することが好ましい。シングルセル化の熟練を必要とせず、大量の細胞を高速に分離・回収できるため、解析精度を高めることができるからである。
【００３９】
単一細胞を回収する際には、蛍光標識、ラジオアイソトープ（ＲＩ）標識、抗体標識、および磁気標識を用いて、各細胞を標識することが好ましい。後述するグルーピングステップ（Ｓ４）において、各細胞集団の細胞種の同定に利用することができるからである。図２Ａ〜２Ｃにおいて、アストロサイト２には、蛍光標識６が付され、神経細胞４には蛍光標識８が付されている。なお、図中では、蛍光標識は、アストロサイト２及び神経細胞４それぞれ、１つの細胞のみに付されているが、全細胞にそれぞれ付されているものとする。
特に、細胞の表面に発現しているタンパク質に結合する抗体と、蛍光、ＲＩ標識、または磁気標識との組み合わせは、抗体の特異性が増すため好ましい。
【００４０】
なお、Ｃ１ＴＭ  Ｓｉｎｇｌｅ−Ｃｅｌｌ  Ａｕｔｏ  Ｐｒｅｐシステム（フリューダイム社製）等のシングルセル・ゲノム研究用自動化ソリューションを用いることもできる。当該ソリューションは、シングルセルの単離、細胞の標識、細胞の溶解、および、後述するシングルセルデータ取得ステップ（Ｓ３）において行うゲノムＤＮＡまたはトータルＲＮＡの抽出までを自動的に行うことができるため、例えば、ゲノムＤＮＡまたはトータルＲＮＡを用いてシングルセルデータを取得しようとする場合に利用すれば、作業効率をより高めることができる。
【００４１】
＜シングルセルデータ取得ステップ（Ｓ３）＞
次に、シングルセルデータ取得ステップ（Ｓ３）において、各時点で回収され、且つシングルセル化された各細胞から、遺伝子のＤＮＡ配列情報、及び遺伝子発現量をシングルセルデータとして取得する。シングルセルデータは、シングルセル化されて回収された全ての単一細胞それぞれについて取得する。
本発明における「遺伝子発現量」とは、遺伝子の転写産物であるｍRNA量であり、遺伝子発現解析により遺伝子の発現状態を調べることにより測定することができる。又は、遺伝子の発現産物であるタンパク質の量の解析を行うようにしてもよい。
【００４２】
《シングルセルデータ》
なお、シングルセルデータとは、単一細胞の機能や性質、状態を示す情報を意味し、上述した遺伝子発現量及びＤＮＡ配列情報に限定されない。例えば、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度などをシングルセルデータとして取得してもよい。
また、単一細胞を回収する際に、各細胞を蛍光物質や、抗体等で標識した場合は、それらの情報もシングルデータとして取得することもできる。
【００４３】
＜グルーピングステップ（Ｓ４）＞
次に、グルーピングステップ（Ｓ４）において、上記シングルセルデータ取得ステップにおいて取得したシングルセルデータ（遺伝子発現量データ、ＤＮＡ配列及び蛍光標識）に基づいて、シングルセルデータを取得した細胞を、共通の第１の細胞特徴を有する細胞集団にグループ分けし、さらに、シングルセルデータに含まれる各細胞種を判別することが可能な第２の細胞特徴に基づいて、第１の細胞特徴でグループ分けした各細胞集団の細胞種を同定する。
具体的には、主成分分析を用いて、遺伝子発現量データに含まれるｎ個（ｎ次元）の細胞特徴を可視化可能な２次元または３次元にまで圧縮処理し、回収された全細胞を複数の細胞集団にグループ分けして２次元平面または３次元空間上にプロットする（つまり、各主成分が「第１の細胞特徴」に該当する）。また、さらに、各細胞集団を構成する細胞の塩基配列及び蛍光標識（第２の細胞特徴）に基づき、グループ分けした各細胞集団の細胞種の同定処理（クラスタリング解析）を行う。
ここで、主成分分析の際の主成分軸は、データ（確率変数）の分散が最大になるような軸を探索することが好ましい。
【００４４】
本ステップにより、各細胞集団の細胞特徴と、細胞数が対応付けられて可視化されるため、細胞数と細胞特徴との関係性を同時に確認又は検出することが容易に、精度高くできる。
【００４５】
なお、回収された全細胞を複数の細胞集団にグループ分けするために使用する細胞特徴の数は、特に限定されない。ｎ個の細胞特徴のうち、１個以上の細胞特徴を、細胞をグループ化するために用いてもよい。
また、本実施形態においては、主成分分析を用いて、２次元または３次元への次元削減により可視化し、シングルセルデータを取得した細胞を共通の第１の細胞特徴を有する細胞集団にグループ分けしているが、これに限定されず、次元削減の方法は、例えば、主成分分析（PCA）、カーネルあり主成分分析（Kernel-PCA）、多次元尺度構成法（MDS）、t-SNE、または、畳込みニューラルネットワーク（CNN）を用いることもできる。
また、各細胞のシングルセルデータ（複数の生体物質に係る情報）から、複数の生体物質の量または状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データを予め用意し、生体物質ごとの時系列データの時間変化と、各生体物質の生物学的機能の類似性に基づいて、シングルセルデータを取得した細胞を、共通の第１の細胞特徴を有する細胞集団にグループ分けしてもよい。
ここで、生物学的機能の類似性は、共通の遺伝子オントロジーを有すること、共通のカノニカルパスウェイに属すること、共通の上流因子を有すること、共通の表現系に関わること、および、共通の疾患に関わることからなる群から選択される少なくとも１つに基づいて評価されるものであることが好ましい。
【００４６】
図３Ａは、図２Ａの細胞刺激直後（時点０）の対象細胞群から取得したシングルセルデータに基づく（主成分及び）クラスタリング結果を示し、図３Ｂは、図２Ｂの細胞刺激後１時間後（時点１）の対象細胞群から取得したシングルセルデータに基づく（主成分及び）クラスタリング結果を示し、図３Ｃは、図２Ｃの細胞刺激後６間後（時点２）の対象細胞群から取得したシングルセルデータに基づく（主成分及び）クラスタリング結果を示す。
図２Ａの細胞刺激直後（時点０）の対象細胞群から取得したシングルセルデータ（遺伝子発現量データ）に主成分分析を行い、２次元に圧縮することで、全細胞が、複数の細胞集団１４及び１６にグループ分けし（図３Ａ）、且つ、シングルセルデータ（ＤＮＡ配列及び蛍光標識６及び８）に基づき、細胞集団１４を構成する細胞種がアストロサイト２であり、細胞集団１６を構成する細胞種が神経細胞４であることを同定する。
同様に、図２Ｂの細胞刺激後、１時間経過後（時点１）の対象細胞群から取得したシングルセルデータ（遺伝子発現量データ）に主成分分析を行い、２次元に圧縮することで、全細胞が、複数の細胞集団１８、２０及び２２にグループ分けし（図３Ｂ）、且つ、シングルセルデータ（ＤＮＡ配列及び蛍光標識６及び８）に基づき、細胞集団１８及び２０を構成する細胞種がアストロサイト２であり、細胞集団２２を構成する細胞種が神経細胞４であることを同定する。
また、図２Ｃの細胞刺激後、６時間経過後（時点２）の対象細胞群から取得したシングルセルデータ（遺伝子発現量データ）に主成分分析を行い、２次元に圧縮することで、全細胞が、複数の細胞集団２４及び２６にグループ分けし（図３Ｃ）、且つ、シングルセルデータ（ＤＮＡ配列及び蛍光標識６及び８）に基づき、細胞集団２４を構成する細胞種がアストロサイト２であり、細胞集団２６を構成する細胞種が神経細胞４であることを同定する。
【００４７】
なお、本実施形態においては、第２の細胞特徴として、ＤＮＡ配列及び蛍光標識を用いて、各細胞集団の細胞種を同定したが、特にこれらに限定されない。細胞の機能（例えば、細胞の増殖、修復、代謝、および細胞間の情報交換）、又は細胞の状態（例えば、遺伝子の発現状況、タンパク質の発現状況、および酵素活性）から各細胞種を判別することが可能な（シングルデータに含まれる）細胞情報、例えば、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度などを利用することができる。
上記各細胞種を判別することが可能な細胞情報には、細胞が本来持っている遺伝子、タンパク質、及び代謝産物等の情報だけでなく、細胞外から導入された遺伝子（例えば、不死化遺伝子）、タンパク質及び代謝物や、有機物等の情報も含む。不死化遺伝子の例としては、ｈＴＥＲＴ遺伝子（ヒトテロメラーゼ逆転写酵素遺伝子）、及びＳＶ４０Ｔ抗原（サルウィルス４０Ｔ抗原遺伝子）が挙げられる。また、シングルセルを回収する際に、各細胞を蛍光物質や、抗体等で標識した場合は、それらの情報も含む。
【００４８】
＜細胞変化検出ステップ（Ｓ５）＞
次に、細胞変化検出ステップ（Ｓ５）は、グルーピングステップ（Ｓ４）で取得されたクラスタリング結果に基づいて、具体的には、細胞刺激直後の細胞に係るクラスタリング結果と、細胞刺激後、所定時間経過した細胞に係るクラスタリング結果との比較することにより、同じ細胞種の細胞集団の経時的な変化（実時間に対する変化、もしくは変化から推定された疑似時間に対する変化）を検出する。なお、細胞集団の経時的変化は、経時的に変化した細胞特徴を抽出し、且つ、その経時的な変化量を数値で検出しておくことが好ましい。
ここで、「細胞集団の経時的変化」とは、細胞数、細胞特徴（第１の細胞特徴）、または、その他細胞特徴（即ち、遺伝子のDNA配列情報（ゲノム）、遺伝子の発現を制御するエピジェネティックな情報（DNAメチル化、ヒストンメチル化、アセチル化、リン酸化）、遺伝子１次転写物（ｍRNA、非翻訳RNA、マイクロRNAなど）情報（トランスクリプトーム）、タンパク質の翻訳量やリン酸化、酸化、糖化等の修飾情報、アミノ酸配列情報（プロテオーム）、代謝産物情報（メタボローム）、細胞内水素イオン濃度指数（pH）、細胞内ATP濃度、イオン濃度（カルシウム、マグネシウム、カリウム、ナトリウムなど）、細胞内温度等）、及び細胞の生体物質の経時的な変化を意味する。
【００４９】
グルーピングステップ（Ｓ４）で取得されたクラスタリング結果である図３Ａ（時点０）、図３Ｂ（時点１）を比較し、同じ細胞種の細胞で構成される細胞集団の細胞特徴（第１の細胞特徴である主成分１及び２）の変化を検出する。具体的には、図３Ａ（時点０）のアストロサイト２と同定された細胞集団１４と、図３Ｂ（時点１）のアストロサイト２と同定された細胞集団１８及び２０との比較、及び、図３Ａ（時点０）の神経細胞４と同定された細胞集団１６の細胞特徴と、図３Ｂ（時点１）の神経細胞４と同定された細胞集団２２との比較により、各細胞集団の細胞特徴の経時的変化の有無を検出する。
具体的には、図３Ａ（時点０）のアストロサイト２と同定された細胞集団１４の細胞特徴（第１の細胞特徴である主成分１及び２）と、図３Ｂ（時点１）のアストロサイト２と同定された細胞集団１８の細胞特徴との比較では、経時的変化がないことを検出し、図３Ａ（時点１）のアストロサイトと同定された細胞集団２０の細胞特徴との比較では、経時的変化があったことを検出する。また、図３Ａ（時点０）の神経細胞４と同定された細胞集団１６の細胞特徴（第１の細胞特徴である主成分１及び２）と、図３Ｂ（時点１）の神経細胞４と同定された細胞集団２２の細胞特徴との比較では、経時的変化がないことを検出する。
【００５０】
また、図３Ｂ（時点１）の細胞集団１８を構成すると同定されたアストロサイト２の細胞数は、図３Ａ（時点０）の細胞集団１４を構成すると同定されたアストロサイト２の細胞数と比べて減少しているが、図３Ｂ（時点１）の細胞集団２０を構成すると同定されたアストロサイト２の細胞数は増加していることを検出する。さらに、図３Ｂ（時点１）の細胞集団１８及び２０を構成するグリア細胞数の合計数は、図３Ａ（時点０）の細胞集団１４を構成すると同定されたアストロサイト２の細胞数に等しいことも検出する。さらに、図３Ｂ（時点１）の細胞集団２２を構成すると同定された神経細胞４の細胞数は、図３Ａ（時点０）の細胞集団１６を構成すると同定された神経細胞４の細胞数は、変化がないことを検出する。
【００５１】
これらの検出結果から、薬剤添加による細胞刺激後、１時時間経過すると、アストロサイト２の細胞特徴及び細胞数に変化がみられるが、神経細胞４には変化が見られないことを容易に検出することができる。さらに、アストロサイト２の細胞特徴及び細胞数の変化から、細胞集団２０が、アストロサイト２の形状が変化した反応性アストロサイト１０から構成されることを検出（推定または確認）することができる。その結果、時点０から時点１にかけて、アストロサイト２の形状が変化し、反応性アストロサイト１０になったことを容易に検出（推定又は確認）することができる（図４実線矢印）。
【００５２】
同様に、グルーピングステップ（Ｓ４）で取得されたクラスタリング結果である図３Ｂ（時点１）及び図３Ｃ（時点２）を比較し、同じ細胞種の細胞で構成される細胞集団の細胞特徴（第１の細胞特徴である主成分１及び２）が時点１から時点２における経時的変化の有無を検出する。
即ち、図３Ｃ（時点２）では、図３Ｂ（時点１）では、アストロサイト２で構成される細胞集団１８が消失していること、及び、図３Ｂ（時点１）で存在する反応性アストロサイト１０で構成されると検出された細胞集団２０だけが残っていることを検出する。また、図３Ｂ（時点１）では存在する神経細胞４で構成される細胞集団２２が消失し、図３Ｃ（時点２）では、神経細胞４で構成されると同定された細胞集団２６が新たに生じていることを検出する。
また、図３Ｂ（時点１）の反応性アストロサイトで構成されると検出された細胞集団２０の細胞特徴（第１の細胞特徴である主成分１及び２）と、時点２のアストロサイト２で構成されると同定された細胞集団２４の細胞特徴とを比較し、経時的変化がないことを検出する。同様に、図３Ｂ（時点１）の神経細胞４で構成されると検出された細胞集団２２の細胞特徴（第１の細胞特徴である主成分１及び２）と、図３Ｃ（時点２）の神経細胞４で構成されると同定された細胞集団２６の細胞特徴とを比較し、経時的変化があることを検出する。
【００５３】
また、さらに、図３Ｂ（時点２）の細胞集団２４を構成すると検出された反応性アストロサイト１０の細胞数は、図３Ａ（時点１）の細胞集団２０及び１８を構成する細胞数の合計数に等しいこと、図３Ｃ（時点２）の細胞集団２６を構成する神経細胞数は、図３Ｂ（時点１）の細胞集団２２の細胞数と同じであることを検出する。
これらの検出結果から、薬剤添加による細胞刺激後、６時時間経過、即ち、時点１から５時間経過すると、反応性アストロサイト１０の細胞特徴及び細胞数には変化が見られないが、神経細胞４には、変化がみられることを容易に検出することができる。さらに、神経細胞４で構成されると検出された細胞集団２２の細胞特徴及び細胞数の変化から、細胞集団２６が神経細胞４の形状が変化した神経細胞１２から構成されることを検出（推定または確認）することができる。その結果、時点１から時点２にかけて神経細胞４の形状が変化し、神経細胞１２になったことを容易に検出（推定または確認）することができる。（図５実線矢印）。
【００５４】
このように、本発明においては、各細胞集団の細胞特徴と、細胞数が対応付けられて表示されるため、細胞死の有無だけでなく、経時的な細胞数や細胞特徴の変化、及び細胞数と細胞特徴との関係性を容易に、精度高く確認又は検出することができる。
【００５５】
＜作用機序解析ステップ（Ｓ６）＞
次に、作用機序解析ステップ（Ｓ６）において、細胞変化検出ステップ（Ｓ５）において取得された細胞変化検出結果に基づいて、同じ細胞種の細胞集団内または細胞集団間の変化の作用機序を解析する。ここで、作用機序とは、薬剤による細胞刺激がその薬理学的効果を発揮するための特異的な作用であり、且つ、同じ細胞種の細胞集団内または細胞集団間でみられる特異的な生化学的反応または相互作用を意味する。
【００５６】
本実施形態において、細胞集団内または細胞集団間の変化の作用機序とは、細胞刺激により細胞内外で誘起された生物学的な現象（例えば、増殖、細胞死、抗原抗体反応、増殖因子の分泌等）であり、より詳細には、細胞刺激により細胞内外で誘起された生物学的な現象を発揮するための特異的な生化学的反応または相互作用（例えば、生体物質の代謝、遺伝子発現、エネルギー代謝、シグナル伝達等）である。
本実施形態においては、細胞変化検出ステップ（Ｓ５）において取得された細胞変化検出結果から、上述したような作用機序を解析することにより、細胞変化検出ステップ（Ｓ５）において取得された各細胞集団の細胞数の変化、及び細胞特徴の変化に関与する因子（例えば、生体物質、遺伝子等）を抽出する。
【００５７】
図６は、シングルセルデータ取得ステップ（Ｓ３）において、先述したＤＮＡ配列及び蛍光標識等とともに、アストロサイト２（及び反応性アストロサイト１０）と神経細胞４及び１２との相互作用に関係すると考えられる成分Ａ〜Ｚの成分量もシングルセルデータとして取得しておき、主成分分析による次元削減後の発現プロファイルの類似度に基づいて推定された疑似的な時間軸（pseudotime）に沿って全細胞をプロットした図である。このように細胞を配置することにより、疑似的な時間的な細胞変化及び遺伝子発現変化がわかる。
図６（１）〜（４）は、アストロサイト２（及び反応性アストロサイト１０）に係る成分Ａ〜Ｃ及びＺの変化を示し、図６（５）〜（８）は、神経細胞４及び１２の成分Ａ〜Ｃ及びＺの変化を示す。
【００５８】
なお、本ステップにおいて使用するアストロサイト２（及び反応性アストロサイト１０）と神経細胞４及び１２との相互作用に関係すると考えられる成分は、特に限定されない。シングルセルデータ取得ステップ（Ｓ３）において、シングルセルデータとして取得可能なものを用いることができる。
【００５９】
本実施形態においては、図６（１）〜（４）から、先述した細胞変化検出ステップ（Ｓ５）で検出された時点０から時点１に係る細胞集団１８（アストロサイト２）から細胞集団２０（反応性アストロサイト１０）への細胞特徴及び細胞数の経時的変化（図４の実線矢印）は、時点０から時点１にかけて成分量が変化しているＢ成分が関係しているという検出結果を獲得する。
また、図６（５）〜（８）から、先述した細胞変化検出ステップ（Ｓ５）で検出された時点１から時点２に係る細胞集団２８、即ち、細胞集団２２（神経細胞４）から細胞集団２６（神経細胞１２）への細胞特徴及び細胞数の経時的変化（図５の実線矢印）は、時点１から時点２にかけて成分量が変化しているＺ成分及びＢ成分が関係しているという検出結果を獲得する。
【００６０】
なお、本実施形態においては、グリア細胞と神経細胞の相互ネットワーク形成の過程（経時的変化）を解析するために使用しているが、これに限定されず、Candida.albicansのような真菌と口腔粘膜上皮細胞との感染過程を解析することにも使用することができる。
【００６１】
上記作用機序解析として、例えば、遺伝子発現プロファイルの解析や、分子ターゲットの解析を用いることもできる。
遺伝子発現プロファイルの解析には、例えば、テンソル分解を用いることができる。
【００６２】
分子ターゲットの解析の例を説明する。
まず、薬剤（細胞刺激）によって発現が変化し、かつ、その変化が疾患によるものと一致していると期待できる遺伝子を特定する。ここで、実際に化合物（薬剤）が結合しているのはタンパクであるが、計測されているのはｍＲＮＡの発現量である。タンパクに化合物（薬剤）が結合することでそのタンパクをコードしている遺伝子のｍＲＮＡの量が変化しているとは思えないので、発現量が変化している遺伝子の中に分子ターゲット（標的タンパク）は無いと推定する。
化合物が実際に結合しているタンパクの他の遺伝子発現プロファイルに対する影響は、そのタンパクがコードされている遺伝子をノックアウトした場合に近いと期待される。そこで遺伝子を網羅的にノックアウトした場合の遺伝子発現プロファイルを参照することで標的遺伝子を推定する。
【００６３】
複数の疾患または症例タイプ候補から細胞刺激の最適な適応例を見出すには、複数の種類の細胞に対して一斉に効果判定および作用機序解析を行うことで、効果の序列を得ることができる。その効果の比較もしくは序列に基づいて、最適な適応症または用途を見出すことができる。ここで、細胞刺激の効果判定は、例えば、細胞刺激（薬剤）処理無のデータと細胞刺激（薬剤）処理有の時間依存データを比較して、同一細胞集団に由来する細胞が生物学的データまたは細胞特徴、細胞数を比較し、変化があった場合は、効果があったと判定できる。
【００６４】
＜相互作用解析ステップ（Ｓ７）＞
次に、相互作用解析ステップ（Ｓ７）において、細胞変化検出ステップ（Ｓ６）において取得された作用機序解析の結果に基づいて、異なる細胞種の細胞集団の間で起きる相互作用を解析する。その結果、同じ細胞種の細胞集団間だけでなく、異なる細胞種の細胞集団の間で起きる相互作用を解析することができるため、より広い細胞間及び細胞種間の相互ネットワークを解析することができる。
【００６５】
上記細胞集団の間で起きる相互作用とは、種類の異なる細胞集団の間で起きる生物学的または物理的な作用、またはそれらの作用により生じる変化である。
上記生物学的または物理的な作用とは、例えば、液性因子の授受、細胞接着、拍動や、または微弱電流であり、上記作用により生じる変化とは、例えば、形態変化、分化、増殖、細胞死、遊走、表面抗原の増加、または液性因子の放出である。
上記液性因子は、例えば、ホルモン、成長因子、サイトカイン、エクソソーム、代謝物質、またはイオン等の無機物である。
【００６６】
上記グルーピングステップ（Ｓ４）において、各細胞のシングルセルデータ（複数の生体物質に係る情報）から、複数の生体物質の量または状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データを予め用意し、生体物質ごとの時系列データの時間変化と、各生体物質の生物学的機能の類似性に基づいて、シングルセルデータを取得した細胞を、共通の第１の細胞特徴を有する細胞集団にグループ分けした場合は、国際公開２０１８/１５０８７８Ａ１に記載の方法を用いて相互作用解析を行うことができる。
本発明においては、相互作用解析ステップ（Ｓ７）において、複数の時点の各々について、複数の細胞集団の各々に含まれる１つ以上の第１の細胞特徴から、細胞集団の状態を表す値を生成し、生成された、複数時点の、複数の細胞集団の状態を表す値を、生体物質ごとにそれぞれ複数の時点において取得した時系列データからなるデータセットから、細胞集団間の状態の依存関係を推定し、そのグループ間（細胞集団間）の依存関係の中において、細胞集団が異なる細胞集団間の依存関係を抜き出すことで異なる細胞種の細胞集団の間で起きる生物学的または物理的な作用、またはそれらの作用（異なる細胞種の細胞集団の間で起きる相互作用）を推定することができる。
【００６７】
細胞集団間の状態の依存関係の推定は、例えば、以下のようにして行うことができる。
細胞を、配列データから同定し、それぞれの遺伝子発現量の時間的な変化パターンが似ている遺伝子を、例えば、７つある隣接２時点における状態値の遷移が、ある閾値に照らして増加した、不変だった、減少した、の３つのうちのどれであるかを判定し、それぞれ３７＝２１８７のグループ（細胞集団）に分類する。
生物学的機能が似ている遺伝子を公共のＷｅｂツールＤＡＶＩＤ（https://david.ncifcrf.gov/）のFunctional Annotation Clusteringを用い、類似した遺伝子オントロジーを有する遺伝子をグループ化する。時間的変化の類似性と、生物学的機能の類似性とに基づいてグループ化し、各グループ（細胞集団）の状態値間の時間的な依存関係を、例えば、ベイジアンネットワークモデルに照らして、または、時系列もしくは生物学的な関係性によってグループ間（細胞集団間）の紐付けをして、推定することができる。
【００６８】
本実施形態においては、先述した細胞変化検出ステップ（Ｓ５）で取得された作用機序解析の結果、即ち、時点０から時点１に係るグリア細胞の細胞特徴及び細胞数の経時的変化（図４の実線矢印）は、時点０から時点１にかけて成分量が変化しているＢ成分が関係していること、及び、時点１から時点２に係る神経細胞の細胞特徴及び細胞数の経時的変化（図５の実線矢印）は、時点１から時点２にかけて成分量が変化しているＺ成分及びＢ成分が関係しているという結果、さらに、図６（１）〜（４）に示されるグリア細胞に係る成分Ａ〜Ｃ及びＺの変化、及び図６（５）〜（８）に示される神経細胞の成分Ａ〜Ｃ及びＺの変化に基づいて、時点１のグリア細胞で構成される細胞集団２０と神経細胞４で構成される細胞集団２２との間で起きる生物学的または物理的な作用（刺激）、またはそれらの作用（刺激）、即ち、異なる細胞種の細胞集団の間で起きる相互作用があること、及び、時点０から時点２にかけて、まず、はじめに、薬剤添加による細胞刺激により、グリア細胞のＢ成分が増加し、グリア細胞の細胞特徴が変化し（図４及び図７実線矢印）、即ち、アストロサイト２が反応性アストロサイト１０への変化し、変化した反応性アストロサイト１０から生物学的または物理的な刺激が神経細胞４なされる（図７太矢印）。続いて、反応性アストロサイト１０から刺激を受けた神経細胞４の成分Ｃが減少することに応じて、神経細胞４のＺ成分が増加することにより、神経細胞の細胞特徴が変化し（図５実線矢印）し、即ち、神経細胞４が神経細胞１２へ変化するという、グリア細胞と神経細胞間の細胞間の相互ネットワークを推定することができる。
【００６９】
実施形態１の方法によれば、異なる複数の細胞種が存在する対象細胞群を構成する全ての細胞について、シングルセル解析を行い、各細胞の細胞種を同定するものであるため、異なる複数の細胞種が存在する対象細胞群から人為的に各細胞の細胞種を同定する従来の方法よりも、各細胞種の同定、各種細胞に対する薬剤または細胞刺激の作用機序を、精度高く、解析することができる。
また、各細胞種の同定、及び各細胞種に対する薬剤又は細胞刺激の作用機序を精度高く解析することができるため、同じ細胞種間の相互作用だけでなく、異なる細胞種間の相互作用も精度高く、容易に解析することができる。また、その結果、より広い細胞間及び細胞種間の相互作用（相互ネットワーク）も精度高く、容易に解析することができる。
また、細胞種の同定及び選択に画像解析等の手間をかけることがなく、また、免疫細胞に対する薬効スクリーニングを行う場合であっても、１つのサンプルの観察を解析が終わるまで続けるものでもないため、特に解析時間に制約を設けることなく解析を行うことができる。
また、シングルセル解析を利用するため、人為的に細胞の変化を確認し、対象細胞を選択する場合よりも、細胞を観察する（回収する）時点数を減らすことができる。
【００７０】
また、実施形態１の方法は、グルーピングステップ（Ｓ４）により、各時点の細胞の状態を可視化できるものであるため、細胞変化検出ステップ（Ｓ５）において、各細胞（細胞集団）の変化を容易に確認することができる。また、その結果、作用機序解析ステップ（Ｓ６）において、さらに解析する必要のある細胞種や、遺伝子等を容易に抽出することができる。
【００７１】
実施形態１の変形例
上述の実施形態１では、上記細胞回収ステップ（Ｓ１）において、細胞刺激直後に回収された細胞と、細胞刺激後、所定時間培養した細胞とを回収し、それら細胞のシングルセルデータの比較に基づく解析を行っているが、これに限定されず、上記細胞回収ステップ（Ｓ１）において、細胞刺激が与えられない対象細胞群と、細胞刺激が与えられる対象細胞群をそれぞれ用意し、それぞれ所定時間培養した後に、細胞を回収し、それら細胞のシングルセルデータの比較に基づくスクリーニング解析を行ってもよい。
【００７２】
具体的には、上記細胞回収ステップ（Ｓ１）において、準備した複数の容器１のうち、一部の容器１に播種された対象細胞群に対しては、細胞刺激を与えずに培養し、残りの容器１に播種された対象細胞群に遺体しては、細胞刺激を与え、それぞれ、２以上の時点で、１つの容器内にある全細胞を回収する作業を行うこと以外は、実施形態１と同様である。
【００７３】
このように、さらに、対照サンプルを用意し解析を行うことにより、薬剤（細胞刺激）による効果であるのか、それとも培養条件等のその他要因による効果であるのかを確認することができる。
より詳細には、作用機序解析ステップ（Ｓ６）において、薬剤（細胞刺激）による個々の細胞集団内または細胞集団間への作用機序を解析することができ、さらに、細胞刺激の分子ターゲットを解析したり、治療を想定した適応例を特定したりすることもできるようになる。
ここで、治療を想定した適応例としては、例えば、細胞刺激により改善が見込める疾患もしくは症状、または分子ターゲットに関連する疾患もしくは症状を挙げることができる。
【実施例】
【００７４】
以下では実施例により本発明をより具体的に説明するが、本発明はこれらの実施例に限定されるものではない。【００７５】
［実施例１］
ヒトｉＰＳ細胞およびマウス胎児由来線維芽細胞を１：１（細胞数比）で混合し、ＲＯＣＫ（Ｒｈｏ−ａｓｓｏｃｉａｔｅｄ  ｃｏｉｌｅｄ−ｃｏｉｌ  ｆｏｒｍｉｎｇ  ｋｉｎａｓｅ／Ｒｈｏ結合キナーゼ）阻害剤Ｙ−２７６３２（１０μＬ；富士フイルム和光純薬社製）を加えて、１０００細胞／ウェル／２００μＬ、３０００細胞／ウェル／２００μＬ、および９０００細胞／ウェル／２００μＬで、６ウェルプレートに播種した。
播種時（０時間）、播種から６時間後、１２時間後、２４時間後、４８時間後、７２時間後、９６時間後、および１２０時間後の細胞を回収した。
細胞分散液を１０００細胞／μＬに希釈し、Ｃ１システム（フリューダイム社製）を用いてシングルセルを捕捉した。次いで、ＳＭＡＲＴｅｒ（Ｒ）  Ｕｌｔｒａ（Ｒ）  Ｌｏｗ  ＲＮＡキットを用いて、細胞の溶解、ｍＲＮＡの逆転写およびｃＤＮＡプレ増幅を行った。
得られたｃＤＮＡを回収し、０．０５ｎｇ／μＬより高い濃度を、ライブラリー調製のために選択した。ライブラリー調製は、Ｎｅｘｔｅｒａ（Ｒ）  ＸＴ  ＤＮＡサンプル調製キット（イルミナ社製）を用いて行った。【００７６】
調製したライブラリーについて、次世代シーケンサー（ＨｉＳｅｑ２５００システム，イルミナ社製）により、２×１００ｂｐの末端読取りを用いて配列決定した。得られたデータから細胞毎の遺伝子発現量を算出し、薬剤毎に時点サンプルをまとめ、主成分分析を用いてクラスタリング解析を行った。
細胞は、配列データからヒトｉＰＳ細胞であるかマウス胎児由来線維芽細胞であるかを同定し、それぞれの遺伝子発現量の時間的な変化パターンが似ている遺伝子を、具体的には、７つある隣接２時点における状態値の遷移が、ある閾値に照らして増加した、不変だった、減少した、の３つのうちのどれであるかを判定し、それぞれ２１８７（＝３の７乗）のグループに分類した。生物学的機能が似ている遺伝子を公共のＷｅｂツールＤＡＶＩＤ(https://david.ncifcrf.gov/)のFunctional Annotation Clusteringを用い、類似した遺伝子オントロジーを有する遺伝子をグループ化した。時間的変化の類似性と、生物学的機能の類似性とに基づいてグループ化した結果、７３０５個のグループが得られた。７３０５個のグループの状態値間の時間的な依存関係を、ベイジアンネットワークモデルに照らし推定したところ、マウス胎児由来線維芽細胞の増殖因子の分泌に伴うヒトｉＰＳ細胞の未分化性維持機構の活性化が時系列的な機序を捉える事ができた。【００７７】
［実施例２］
ヒト初代神経細胞、ＳＶ４０Ｔ抗原を導入した不死化ヒトミクログリア、およびｈＴＥＲＴ遺伝子を導入した不死化アストロサイトを２：１：１（細胞数比）で混合して６ウェルプレートに播種し、２４時間培養した。
生着を確認後に、生理食塩水、ＴＮＦ−α（腫瘍壊死因子α）、ＩＬ−１β（インターロイキン１β）、ＩＬ−１ｒａ（インターロイキン１受容体アンタゴニスト）、ＩＬ−１２（インターロイキン１２）、ＩＦＮγ（インターフェロンγ）、ＩＲＦ１（インターフェロン制御因子１）、ＩＲＦ２（インターフェロン制御因子２）、ＩＲＦ３（インターフェロン制御因子３）、ＩＲＦ４（インターフェロン制御因子４）、ＩＲＦ５（インターフェロン制御因子５）、ＩＲＦ６（インターフェロン制御因子６）、ＩＲＦ７（インターフェロン制御因子７）、ＩＲＦ８（インターフェロン制御因子８）、ＩＲＦ９（インターフェロン制御因子９）、およびＬＰＳ（リポ多糖）を各ウェルに添加した。
添加時（０時間）、添加から６時間後、１２時間後、２４時間後、および４８時間後の時点で細胞をトリプシン処理して回収した。
細胞分散液を１０００細胞／μＬに希釈し、Ｃ１システム（フリューダイム社製）を用いてシングルセルを捕捉した。【００７８】
ｃＤＮＡ調製キット（ＳＭＡＲＴｅｒ（Ｒ）  Ｕｌｔｒａ（Ｒ）  Ｌｏｗ  ＲＮＡキット，クロンテック社製）にＴＰ５３遺伝子の５００−６００および７５０−８７０のターゲットプライマーを加えて、細胞の溶解、ｍＲＮＡの逆転写およびｃＤＮＡプレ増幅を行った。
得られたｃＤＮＡを回収し、０．０５ｎｇ／μＬより高い濃度を、ライブラリー調製のために選択した。ライブラリー調製は、Ｎｅｘｔｅｒａ（Ｒ）  ＸＴ  ＤＮＡサンプル調製キット（イルミナ社製）を用いて行った。
調製したライブラリーについて、次世代シーケンサー（ＨｉＳｅｑ２５００システム，イルミナ社製）により、２×１００ｂｐの末端読取りを用いて配列決定した。
得られたデータから細胞毎の遺伝子発現量を算出し、薬剤毎に時点サンプルをまとめ、主成分分析を用いてクラスタリング解析を行った。【００７９】
細胞は、ＳＶ４０Ｔ抗原遺伝子が検出された細胞を不死化ミクログリアクラスターに、ｈＴＥＲＴ遺伝子が検出された細胞を不死化アストロサイトクラスターに、ＳＶ４０Ｔ抗原遺伝子およびｈＴＥＲＴ遺伝子のいずれもが検出されない細胞をヒト初代神経細胞クラスターに、それぞれグルーピングした。
細胞の種類ごとに、時間経過による細胞数変動および遺伝子発現量変動を見出し、遺伝子の時間変動パターンを取得した。各細胞各遺伝子時間変動パターンをまとめてグループ化およびパターン化し、さらに時間的な依存関係を算出することで、細胞間の依存関係を取得した。【００８０】
以上、本発明の細胞情報処理方法についての種々の実施形態及び実施例を挙げて詳細に説明したが、本発明は、これらの実施形態及び実施例に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良又は変更をしてもよいのはもちろんである。
【図１】
2020071500.tif 000003
【図２Ａ】
2020071500.tif 000004
【図２Ｂ】
2020071500.tif 000005
【図２Ｃ】
2020071500.tif 000006
【図３Ａ】
2020071500.tif 000007
【図３Ｂ】
2020071500.tif 000008
【図３Ｃ】
2020071500.tif 000009
【図４】
2020071500.tif 000010
【図５】
2020071500.tif 000011
【図６】
2020071500.tif 000012
【図７】
2020071500.tif 000013
【図８】
2020071500.tif 000014
【国際調査報告】
2020071500.tif 000015
2020071500.tif 000016
2020071500.tif 000017
2020071500.tif 000018