JP7548672B2 - 細胞を分析する方法 - Google Patents
細胞を分析する方法 Download PDFInfo
- Publication number
- JP7548672B2 JP7548672B2 JP2021501001A JP2021501001A JP7548672B2 JP 7548672 B2 JP7548672 B2 JP 7548672B2 JP 2021501001 A JP2021501001 A JP 2021501001A JP 2021501001 A JP2021501001 A JP 2021501001A JP 7548672 B2 JP7548672 B2 JP 7548672B2
- Authority
- JP
- Japan
- Prior art keywords
- cellular
- cell
- expression
- cells
- perturbation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N15/1429—Signal processing
- G01N15/1433—Signal processing using image recognition
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6881—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for tissue or cell typing, e.g. human leukocyte antigen [HLA] probes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/28—Databases characterised by their database models, e.g. relational or object models
- G06F16/284—Relational databases
- G06F16/288—Entity relationship models
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional [2D] or three-dimensional [3D] molecular structures, e.g. structural or functional relations or structure alignment
- G16B15/30—Drug targeting using structural data; Docking or binding prediction
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/10—Signal processing, e.g. from mass spectrometry [MS] or from PCR
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/30—Unsupervised data analysis
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N2015/1006—Investigating individual particles for cytology
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N2015/1402—Data analysis by thresholding or gating operations performed on the acquired signals or stored data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/10—Investigating individual particles
- G01N15/14—Optical investigation techniques, e.g. flow cytometry
- G01N2015/1497—Particle shape
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Biotechnology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Analytical Chemistry (AREA)
- Molecular Biology (AREA)
- Databases & Information Systems (AREA)
- Genetics & Genomics (AREA)
- Immunology (AREA)
- Data Mining & Analysis (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Biochemistry (AREA)
- Bioethics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Epidemiology (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Microbiology (AREA)
- Signal Processing (AREA)
- Dispersion Chemistry (AREA)
- Pathology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Pharmacology & Pharmacy (AREA)
Description
本出願は、2018年7月16日に出願された米国仮特許出願第62/698,701号、2019年2月14日に出願された同第62/805,884号、及び2019年2月14日に出願された同第62/805,888号に対する優先権の利益を主張するものであり、これらの仮特許出願の全内容が、参照により本明細書に援用される。
本開示のいくつかの態様の概説及び本開示において使用されるいくつかの定義が提供されてきたが、例示的なシステムの詳細が、図1とともに記載される。
・様々な基本システムサービスを処理するため及びハードウェアに依存するタスクを実行するための手順を含む、オペレーティングシステム116(例えば、ANDROID、iOS、DARWIN、RTXC、LINUX、UNIX、OS X、WINDOWS、又はVxWorkなどの組み込みオペレーティングシステム);
・システム100を、他のデバイス及び/又は通信ネットワーク104と接続するための任意選択のネットワーク通信モジュール(又は命令)118;
・複数のデータセット122を保存するデータセットストア120(各データセットは、1つ以上の識別子(例えば、サンプル識別子124及び/又は細胞/データセット識別子126)、関連する期間128、及び1つ以上の細胞成分132を含む細胞成分ベクトル130を含む);並びに
・1つ以上の単一細胞の移行シグネチャー142及び1つ以上の摂動シグネチャー150を保存するシグネチャーストア140。
細胞の運命の原因を同定するための本明細書に記載される技術を行う際、単一細胞から得られる細胞成分測定に関するデータセットを生成することが有用である。これらのデータセットを生成する(例えば、図3の単一細胞の測定310によって図1のデータセット122-1を生成する)ために、対象とする細胞の集団が、インビトロで培養される。対象とする1つ以上の細胞成分132の単一細胞の測定が、培養中に1つ以上の期間で行われて、データセット122を生成する。(例えば、図3の単一細胞の測定310)。ある実施形態において、対象とする細胞成分としては、DNA、修飾(例えば、メチル化)DNAを含む核酸、コーディング(例えば、mRNA)若しくはノンコーディングRNA(例えば、sncRNA)を含むRNA、転写後修飾タンパク質(例えば、リン酸化、グリコシル化、ミリストイル化タンパク質など)を含むタンパク質、脂質、炭水化物、環状アデノシン一リン酸(cAMP)及び環状グアノシン一リン酸(cGMP)などの環状ヌクレオチドを含むヌクレオチド(例えば、アデノシン三リン酸(ATP)、アデノシン二リン酸(ADP)及びアデノシン一リン酸(AMP))、ニコチンアミドアデニンジヌクレオチド(NADP/NADPH)の酸化及び還元形態などの他の小分子細胞成分、及びそれらの任意の組合せが挙げられる。ある実施形態において、細胞成分測定は、RNAレベルなどの遺伝子発現測定を含む。
0日目:第1の細胞状態のいくつかの細胞を、細胞の増殖に好適な培地中でプレート上に解凍する。
1日目:第1の細胞状態のある数の細胞を、マルチウェルプレート中に播種する。該当する場合、細胞の細胞成分に影響を与えるためのさらなる工程を行う。例えば、対象とする細胞成分をノックアウトするために1つ以上のウイルスに同時に感染させる。
ウェル中の細胞について細胞成分の発現測定の反復t1を行う。
1+l日目:何らかのさらなるプロセスが行われた場合、必要に応じて培地を交換する。
該当する場合、ウェル中の細胞について細胞成分の発現測定の反復tlを行う。
1+m日目:変化した細胞状態の細胞の増殖を補助するのに適切な培地に、培地を交換する。
該当する場合、ウェル中の細胞について細胞成分の発現測定の反復tmを行う。
1+n、o、p日目など:第1の細胞状態から変化した細胞状態へのさらなる細胞状態の移行を補助するために必要に応じて培地を交換する。該当する場合、第1の細胞状態から変化した細胞状態へのさらなる移行に影響を与えるためのさらなる工程を行う。例えば、細胞を変化した細胞状態にするための、対象とする摂動を加える。
該当する場合、ウェル中の細胞について、細胞成分の発現測定の反復tn、to、tpなどを行う。
q日目:ウェル中の変化した細胞状態の細胞について細胞成分の発現測定の反復tqを行う。
プレートを固定し、対象とする細胞成分/タンパク質に適合した抗体で染色して、測定されるべきそれらを溶解/破壊させずに細胞をソート/同定する。それはまた、細胞質の状況におけるほどの分解を伴って見られない表面マーカーを同定するのに使用され得る。各ウェルを走査することによるMolecular Devices HCI IXM4などの細胞イメージングシステムを用いでイメージングする。所望の変化した細胞状態にあるウェル当たりの細胞の数を定量する。
III.A.概説
細胞状態の移行(すなわち、第1の細胞状態から変化した細胞状態への細胞の状態の移行)は、細胞における細胞成分132の発現の変化によって特徴付けられる。例えば、移行は、細胞における細胞成分の発現132の変化、ひいては細胞によって産生される細胞成分(例えば、mRNA、転写因子)の同定及び定量によって特徴付けられ得る。しかしながら、少なくとも現在、細胞状態の移行は、細胞内活性の複雑性のため、完全に決定性ではない。この複雑性を理解しようとして、本明細書は、細胞状態の移行の様々な段階での対象とする1つ以上の測定される細胞成分の可変の存在、非存在又は量に関連する可変の細胞成分の発現が、有意義な知識が抽出され得る高次元データセット(例えば、図1の細胞成分ベクトル130)を提供するという理論の下で、細胞の集団の細胞中の細胞成分132を定量する単一細胞データセット122に統計的技術を適用する。本明細書において、データの高次元は、データセット122に含まれる前駆細胞成分の測定に由来する。各細胞成分132は、次元を表し、集合的に、各細胞成分についての細胞成分測定データセット122が、「前駆」細胞の、異なる細胞型への移行に関する生物学的過程についての潜在情報を符号化する形状を有し得る。実際には、細胞成分132の数は、およそ数万であり得るため、本明細書に記載される計算は、知的に又は手動によって行うことが不可能ではないにしても実用的ではない。
どの種類の方法が使用されるかにかかわらず、異なって発現された細胞成分の決定は、どのような結果が求められるかに応じて変化し得る。例えば、使用される方法が特定の細胞を系統内(on-lineage)又は系統外(off-lineage)であると同定する場合、どの細胞成分が異なって発現されるかの決定は、系統内(on-lineage)であると決定された細胞の細胞成分の発現レベルを、系統外(off-lineage)であると決定された細胞の細胞成分と比較することによって行われ得る。それらの細胞成分の相対的発現は、どの細胞成分が、個別に又は組み合わされて、1つの型又は別の型の細胞中で活性であるかを示す。上記のように、この発現データは、異なって発現されるとフラグを付けられる細胞成分のサブセットを同定するのに使用され得る。次に、因果関係が、同定された細胞成分をインビトロでノックアウトし、実験細胞集団の細胞の運命が、どの細胞成分が活性であるかの変化によって影響されるか否かを評価することによって決定され得る。
図3は、一実施形態に係る、異なって発現された細胞成分132のセットを決定するための細胞成分の発現の差のアッセイの第1の例についてのフローチャートである。図3は、細胞成分の発現の差を用いて記載される一般的な場合の非限定的な例示的な実施形態を提供することに留意されたい。工程310において、単一細胞発現測定が、上記の第II節に記載されるように行われて、細胞の集団についての複数のデータセット122を生成する。上記のように、各細胞についての各データセット122が、細胞成分(例えば、図1の細胞成分132-1-1~132-1-B)のそれぞれの量を含む細胞成分のベクトルri(例えば、図1の細胞成分ベクトル130)として表され得る。単一細胞発現測定310から得られるデータセット122は、一般に、計算装置(例えば、図1のシステム100)の永続的メモリ(例えば、図1の永続的メモリ112)中にデジタル形式で保存されるが、それらは、本明細書に記載される残りの工程を行うために、必要に応じて能動メモリ(例えば、図1の非永続的メモリ111)中にロードされ得る。一般に、図3のプロセスの残りの工程は、1つ以上の計算装置(例えば、図1のシステム100)によって行われる。計算装置の例が、図1に関して説明される。しかしながら、実際には、図3のプロセスは、本明細書に記載される工程の結果に基づいて行われるさらなるインビトロ試験又は臨床決定などの、コンピュータの外で行われ得るさらなる中間の又は後続の工程を含み得る。
上で紹介されるように、細胞成分132のそれぞれが、データの異なる次元を表す際、データセット122は、全体で、高次元を有する。工程320において、次元削減が、計算装置(例えば、システム100)によって行われて、データセット122の細胞成分132の量中に存在する任意の潜在的パターンの構造を保ちながらデータの次元を削減する。
行列M中の削減された次元データ(例えば、次元削減成分ストア146)は、単一細胞発現データセット122から元の高次元データと比べて次元が有意に削減される。しかしながら、得られる行列Mは、非線形多様体(例えば、図1の多様体149)を埋め込む。工程330において、多様体学習技術は、多様体を抽出するために行列Mに適用される。多様体149自体が、擬似時間にわたって細胞の間で細胞成分の発現の差についての有用な情報を提供するだけでなく、それはまた、その情報を視覚化するのにも使用され得る。
工程340において、次元削減320によって提供される低次元空間中の点の位置におけるパターン(例えば、関連する複数の次元削減ベクトル146のサブセットに対応する)を同定するために、クラスタリングが行われて、jクラスターのセットCjを生成する。これらのクラスターを用いて、同様の点(細胞/データセット)を統合して、低次元空間において互いに類似している点の群(例えば、第1のクラスター、第2のクラスターなど)についての統計的に関連する情報を引き出す。以下の表4は、ラスタリング340の出力であり得る点のクラスタリングの例を示す。
次元削減320、任意選択の多様体学習330、及びクラスタリング340の工程は、一般に、集団の細胞、及びそれらの対応する単一細胞発現データセット122を、削減された次元空間内のクラスターへと組織化するように動作し、それによって、基本的な細胞成分ごとの発現測定データが、有意義な情報を抽出するために統合及び分析され得るようになっている。ある実施形態において、この削減された次元空間は、本開示の方法を完成させるのに必要な時間及び/又は処理能力をさらに減少させる。
異なって発現された細胞成分Ekのセットは、それ自体で有用な出力を表す。しかしながら、そのセットのサブセットを同定するために異なって発現された細胞成分のセットをさらに分析する(360)ことが有用であり得る。
特定の細胞の移行を特徴付ける細胞成分の発現の差を、摂動への細胞の曝露によって引き起こされる細胞成分の発現の差にマッチさせることによって、特定の細胞状態の移行に影響を与える摂動が予測され得る。細胞の摂動は、1つ以上の化合物による細胞の任意の処理を含む。1つ以上の化合物は、例えば、小分子、生物製剤、タンパク質、小分子と組み合わされたタンパク質、ADC、siRNA若しくは干渉RNAなどの核酸、cDNA過剰発現野生型及び/又は突然変異型shRNA、cDNA過剰発現野生型及び/又は突然変異型ガイドRNA(例えば、Cas9システム又は他の細胞成分編集システム)、又は上記のいずれかの任意の組合せを含み得る。特定の細胞の移行について異なって発現された細胞成分は、摂動への細胞の曝露によって引き起こされる異なって発現された細胞成分と比較され得る。次に、特定の細胞の移行の細胞成分の発現の差にマッチする細胞成分の発現の差を引き起こす摂動は、特定の細胞の移行に影響を与えることが予測され得る。
ある実施形態において、開示される方法は、摂動についての生物学的有用性を同定するために使用される。これらの方法は、異なる状態又は表現型、例えば、病変及び正常表現型を有する細胞内に異なって存在することが示され得る任意の細胞成分(又は異なる細胞成分の組合せ)の測定を包含する。すなわち、細胞成分の存在、非存在、又は量は、細胞状態又は表現型に関連する。一実施形態において、本方法は、複数の細胞を摂動に曝す工程と;第1の細胞成分の発現の差のアッセイを行う工程であって、アッセイが、摂動への細胞の曝露の前及び後に複数の細胞から得られた第1の複数の単一細胞発現データセットにアクセスする工程を含み、データセットのそれぞれが、細胞成分のベクトルriを含み、ベクトル中の各エントリが、複数の細胞成分の1つに関連し、各エントリの値が、細胞についての細胞成分の量を表す、工程と;統計技術を、第1の複数のデータセットに適用して、摂動への曝露に応答して異なって発現された細胞成分Ekのセットを生成する工程と;摂動への曝露に応答して異なって発現された細胞成分Ekのセットと、病変細胞表現型と正常細胞表現型との間に差に関連する異なって発現された細胞成分Elのセットとの間の類似度を決定する工程とを含み、ここで、EkとElとの間の有意な類似度が、病変細胞表現型と正常細胞表現型との間で細胞を移行させる摂動についての有用性を示す。
実施形態1.複数の単一細胞細胞成分の発現データセットにアクセスする工程であって、各データセットが、同じ「前駆」細胞型から移行された複数の細胞の1つの細胞から得られ、各データセットが、細胞成分のベクトルriを含み、ベクトル中の各エントリが、複数の細胞成分の1つに関連し、各エントリの値が、細胞についての細胞成分の量を表す、工程と;データセットに対して次元削減を行って、行列Mを生成する工程であって、行列Mが、第1の次元における行及び第2の次元における列を含み、各行が、複数の細胞の1つの細胞に対応し、行列Mの値が、第1及び第2の次元空間中のその点に位置する細胞成分の量から生成される値を含む、工程と;クラスタリングを行って、クラスターのセットCjを生成する工程であって、各クラスターが、行列M中の行のサブセットに対応する複数の点及びそれらの対応する細胞を含む、工程と;クラスターのセットCjを用いて、細胞についての、異なって発現された細胞成分Ekのセットを決定する工程とを含む方法。
以下の実施例は、上記の第II及びIII節において紹介される方法を実証する。さらに詳細には、実施例は、細胞状態の移行の軌道に影響を与えることが知られている遺伝子及び/又は摂動を正確に同定する第II及びIII節の方法の能力を実証する。さらに、以下に説明される実施例は、細胞状態の移行の軌道を制御するのに使用され得る新規な生物学的洞察を生成する第II及びIII節の方法の能力を実証する。具体的には、実施例は、以前に知られていない、細胞状態の移行に影響を与える因子(例えば、遺伝子及び摂動)を同定する第II及びIII節の方法の能力を実証する。
この節は、上述されるインビトロ実験のためのプロトコルを説明する。このインビトロ実験からのデータを、公表されているデータとともにプールして、図4B-5A及び7A-9を生成し、独立して使用して、図6を生成した。
0日目:100万個のMEF細胞を、MEF培地中で10cmプレート中に解凍する。
1日目:24ウェルプレート中に20K/ウェルで播種する。
該当する場合、Ascl1ウイルスに同時にスピン感染させる(spin infect)(感染多重度(MOI)8)。MEF培地(250μl/ウェル)及びポリブレン8μg/mlの存在下で、32℃で1時間にわたって2000rpmで回転させる。
単一細胞リボ核酸(RNA)シーケンシング(scRNA-seq)を行って、各細胞についてのd2データセットを得る。
2日目:ウイルス実験のためにポリブレン(MEF培地)を洗い落とすために培地を交換する。
摂動実験のために、小分子(ジメチルスルホキシド(DMSO)又はエタノール中で再懸濁された)を加える。
3日目:培地を神経細胞用培地に交換する。
摂動実験のために、分子(DMSO又はエタノール中で再懸濁された)を加える。
5日目:半分培地交換(該当する場合、小分子を加える)
8日目:半分培地交換(該当する場合、小分子を加える)
9日目:半分培地交換(該当する場合、小分子を加える)
11日目:半分培地交換(該当する場合、小分子を加える)
13日目:半分培地交換(該当する場合、小分子を加える)
15日目:プレートを固定し、Map2及びTuj1抗体で染色する。各ウェルを走査することによって、Molecular Devices HCI IXM4又は他のハイコンテントイメージング顕微鏡においてイメージングする。ウェルごとのMap2/Tuj1陽性ニューロンの数を定量する。
図4Aは、一実施形態に係る、ある期間にわたって誘導された細胞状態の移行の軌道を辿るタイムラインを示す。より具体的には、図4Aは、0日目から22日目までの23日の期間にわたるMEFの誘導された移行の軌道を辿るタイムラインを示す。
上述されるように、Ascl1のみが過剰発現されたMEFから2、5、及び22日目に得られた遺伝子発現測定を、BAM因子の全てが過剰発現されたMEFから22日目に取られた公表されている遺伝子発現測定とともにプールした。第II節において上述される方法を用いて、細胞における遺伝子発現を測定した日のそれぞれについて、細胞のそれぞれについての遺伝子発現測定を用いて、転写物のベクトルriのデータセットを生成した。転写物の各ベクトルriは、転写物のベクトルriに含まれる遺伝子発現測定が得られた特定の日における特定の細胞に関連していた。転写物のベクトルriにおける各転写物は、細胞のゲノム中の特定の遺伝子に関連しており、転写物のベクトルriにおける各エントリの値は、転写物のベクトルriに関連していた特定の日における転写物のシーケンシング深さ(転写物数)を表した。
第III.C.節において上述されるように、次元削減による行列Mの生成の後、クラスタリングを行って、行列M中のデータをグループ分けして、クラスターのセットCjを生成する。クラスターのセットCj中の各クラスターは、点のセットを含む。
第III.D及びIII.E節において説明されるように、細胞状態の移行に影響を与える遺伝子及び転写因子の同定を可能にすることに加えて、第II及びIII節の方法はまた、細胞状態の移行に影響を与える、小分子などの摂動の同定も可能にする。まず、細胞を、移行の特定の軌道を辿るように誘導する摂動を同定するために、移行の可能な軌道が同定される。
この実施例の実験は、ニューロン及び/又は前駆細胞を促進するための方法を実証した。本明細書に記載される実験において、線維芽細胞(すなわち、初代マウス線維芽細胞)の開始集団を、Ascl1過剰発現レンチウイルスを含む組成物に曝した。48時間後、化合物(例えば、フォルスコリン、グレサチニブ、PD-0325901)、又はビヒクル(すなわち、DMSO又はエタノール)を、組成物に加えた。ニューロンの総数を、陽性Tuj1/Map2シグナル及び神経形態に基づいて手動で計数した。各実験について、各処理条件についてのニューロンの総数を、その実験に対する、DMSOで処理されたウェル中のニューロンの数によって正規化した。図10A及び図10Bに示されるように、線維芽細胞の開始集団から発生したニューロンの存在が、これらの実験において検出された。ニューロンの総数及びニューロンのパーセンテージの両方の倍率変化は、組成物に加えられる化合物に応じて、増加するか、減少するか、又は同じままであった。これらの実験は、本発明の方法が、線維芽細胞を含む細胞の開始集団からニューロン及び/又は前駆細胞を促進するのに有用であることを実証した。
2回目の継代の初代マウス胚性線維芽細胞(MEF)を、DMEM中10%のFBS、1×Glutamax、1×MEMの非必須アミノ酸、1mMのピルビン酸ナトリウム、0.05U/mlのpen/strep、及び55μMのβ-メルカプトエタノールを含むMEF培養培地中で、20,000~45,000個/ウェル(ロットに応じて)で、24ウェルプレート上で平板培養した。培養物中で24時間後、スピンフェクション(spinfection)(90分間にわたって32℃で、2000rpmでプレートを回転させる)によって、8μg/mlのポリブレンを含有するMEF培養培地中で、MEFを、Ascl1過剰発現レンチウイルスに感染させた。レンチウイルス生成については以下を参照されたい。48時間後、培地を、DMEM/F12、1%のN2、2%のB27 1:50、1×Glutamax、25μg/mlのインスリン、化合物又はビヒクル(DMSO又はエタノール)を含有する0.05U/mlのpen/strepを含む神経細胞用培地に交換した。化合物及びそれらの濃度を、以下から選択した:BI-2536(200nM)、シロスタゾール(1000nM)、ダブラフェニブ(2500nM)、エストラジオールシピオナート(2000nM)、EX-527(5000nM)、フェドラチニブ(1000nM)、フォレチニブ(200nM)、フォルスコリン(5000nM)、グレサチニブ(2500nM)、インジルビン-3-オキシム(2000nM)、KI20227(250nM)、KU 0060648(200nM)、m-3M3FBS(1000nM)、マヌマイシン(800nM)、PD-0325901(5000nM)、PHA-665752(1000nM)、キナクリン(200nM)、ロットレリン(1000nM)、セルメチニブ(100nM)、トログリタゾン(5000nM)、及びベムラフェニブ(5000nM)。半分培地交換を、補充される化合物とともに2~3日ごとに行った。
Ascl1感染の12日後の時点で、細胞を、4%のパラホルムアルデヒドで固定し、透過処理し(0.2%のTriton X100)、5%の血清(ロバ、子牛、ヤギ血清混合物)中でブロックし、4℃で一晩、又は室温で2時間にわたってウサギ抗Tuj1(1:1000)及びマウス抗Map2(1:500)抗体で染色した後、二次抗体及びDAPI染色を行った。
イメージングを、Molecular Devices ImageXpress Microにおいて行い;ウェルごとに36の画像を、10倍対物レンズで撮影した。ニューロンの総数を、陽性Tuj1/Map2シグナル及び神経形態に基づいて手動で計数した。各実験について、各処理条件についてのニューロンの総数を、その実験について、DMSOで処理されたウェル中のニューロンの数によって正規化した。
パッケージングプラスミド(Packaging plasmid)(SystemsBio、LV510A-1)又は同様のもの、及びAscl1過剰発現プラスミド(Origeneレンチウイルス発現ベクトルへとクローニングされたAscl1 cDNA cat# PS100064)とともに、Mirus TransIT Lenti Tranfection Reagent(Mirus、MIR 6603)によって293T細胞をトランスフェクトすることによって、レンチウイルスをパッケージングし、16,500RPMで1.5時間にわたってBeckmanCoulter超遠心機において濃縮した。48時間の時点でウサギ抗Ascl1(1:200;Abcam、ab74065-100UG)免疫蛍光染色によって判定した際に、90%以上の細胞のレンチウイルス感染を有する実験のみを追跡した。
実施形態1.摂動が細胞の移行に影響を与えるかどうかを予測するための方法であって、方法が、メモリ及び1つ以上のプロセッサを含むコンピュータシステムにおいて:第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーに電子的形態でアクセスする工程であって、変化した細胞状態は、第1の細胞状態から変化した細胞状態への細胞の移行を経て生じ、単一細胞の移行シグネチャーは、複数の細胞成分の同定及び、複数の細胞成分中の各それぞれの細胞成分について、それぞれの細胞成分の発現の変化と、第1の細胞状態と変化した細胞状態との間の細胞状態の変化との間の関連性を定量する対応する第1の重要度スコアを含む、工程と;複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーに電子的形態でアクセスする工程であって、摂動シグネチャーは、複数の細胞成分の全て又は一部の同定及び、複数の細胞成分の全て又は一部の中の各それぞれの細胞成分について、(i)複数の非摂動細胞と複数の摂動細胞との間のそれぞれの細胞成分の発現の変化と、(ii)複数の非摂動細胞と複数の摂動細胞との間の細胞状態の変化との間の関連性を定量する対応する第2の重要度スコアを含む、工程と;単一細胞の移行シグネチャー及び摂動シグネチャーを比較し、それによって、摂動が細胞の移行に影響を与えるかどうかを決定する工程とを含む、方法。
本明細書に引用される参考文献は、各個々の刊行物又は特許若しくは特許出願が、あらゆる目的のために全体が参照により援用されることが具体的に且つ個別に示されているのと同程度に、あらゆる目的のために全体が参照により本明細書に援用される。
本開示は以下の態様を含む。
項目1: 摂動が細胞の移行に影響を与えるかどうかを予測するための方法であって、前記方法が、
メモリ及び1つ以上のプロセッサを含むコンピュータシステムにおいて:
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーに電子的形態でアクセスする工程であって、前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、前記単一細胞の移行シグネチャーは、複数の細胞成分の同定及び、前記複数の細胞成分中の各それぞれの細胞成分について、前記それぞれの細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化との間の関連性を定量する対応する第1の重要度スコアを含む、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーに電子的形態でアクセスする工程であって、前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定及び、前記複数の細胞成分の全て又は一部の中の各それぞれの細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間のそれぞれの細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化との間の関連性を定量する対応する第2の重要度スコアを含む、工程と;
前記単一細胞の移行シグネチャー及び前記摂動シグネチャーを比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程とを含む、方法。
項目2: 前記単一細胞の移行シグネチャーにアクセスする工程が、
(i)第1の複数の第1の単一細胞細胞成分の発現データセット、及び(ii)第2の複数の第2の単一細胞細胞成分の発現データセットに基づいて、前記単一細胞の移行シグネチャーを決定する工程を含み、ここで:
前記第1の複数の第1の単一細胞細胞成分の発現データセット中の各それぞれの第1の単一細胞細胞成分の発現データセットが、前記第1の細胞状態の第1の複数の細胞の対応する単一細胞から得られ、
前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各それぞれの第2の単一細胞細胞成分の発現データセットが、前記変化した細胞状態の第2の複数の細胞の対応する単一細胞から得られる、項目1に記載の方法。
項目3: 前記第1の複数の単一細胞細胞成分の発現データセットの各それぞれのデータセットが、第1の複数の細胞成分ベクトル中の、対応する細胞成分ベクトルを含み、
前記第2の複数の単一細胞細胞成分の発現データセットの各それぞれのデータセットが、第2の複数の細胞成分ベクトル中の、対応する細胞成分ベクトルを含み、
前記第1及び第2の複数の細胞成分ベクトル中の各それぞれの細胞成分ベクトルが、複数の要素を含み、前記それぞれの細胞成分ベクトル中の各それぞれの要素が、前記複数の細胞成分中の対応する細胞成分に関連しており、前記第1及び第2の複数の単一細胞細胞成分の発現データセットの前記それぞれのデータセットによって表される前記対応する単一細胞についての前記対応する細胞成分の量を表す対応する値を含む、項目2に記載の方法。
項目4: 前記第1及び/又は前記第2の複数の単一細胞細胞成分の発現データセットに対して次元削減を行って、複数の次元削減成分を生成する工程と;
前記第1及び第2の複数の細胞成分ベクトル中の各それぞれの細胞成分ベクトルについて、前記複数の次元削減成分を、前記それぞれの細胞成分ベクトルに適用して、前記複数の次元削減成分中の各それぞれの次元削減成分についての次元削減成分値を含む対応する次元削減ベクトルを形成し、それによって、対応する第1及び第2の複数の次元削減ベクトルを形成する工程と;
クラスタリングを行って、クラスターのセットC j を生成する工程であって、各クラスターが、前記第1及び第2の複数の次元削減ベクトルのサブセットに対応する複数の点を含む、工程と;
前記クラスターのセットC j の第1のクラスターから前記第1の複数の細胞を同定する工程と;
前記クラスターのセットC j の第2のクラスターから前記第2の複数の細胞を同定する工程とをさらに含み、
前記方法が、任意に、前記対応する第1及び第2の複数の次元削減ベクトルを用いて多様体学習を行って、前記第1及び第2の複数の細胞における各細胞の他の各細胞に対する相対的な細胞状態を同定する工程をさらに含む、項目3に記載の方法。
項目5: 前記複数の非摂動細胞が、摂動に曝されていない対照細胞であり、又は前記非摂動細胞が、摂動に曝された非関連摂動細胞の平均である、項目1に記載の方法。
項目6: 前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをプルーニングして、前記複数の細胞成分を、任意にRNAレベルで測定される転写因子に制限する工程をさらに含む、項目1に記載の方法。
項目7: 前記単一細胞の移行シグネチャーを決定する工程が、
母平均の差の検定、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び一般化線形モデルのうちの1つを用いて、(i)前記第1の複数の第1の単一細胞細胞成分の発現データセットと、前記第2の複数の第2の単一細胞細胞成分の発現データセットとの間の前記複数の細胞成分にわたる細胞成分量の差を決定する工程を含む、項目2に記載の方法。
項目8: 細胞成分の発現の差の尺度が、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び一般化線形モデルのうちの1つを用いて、(i)第3の複数の第3の単一細胞細胞成分の発現データセットと、(ii)第4の複数の第4の単一細胞細胞成分の発現データセットとの間の細胞成分量の差を定量し、ここで:
前記第3の複数の第3の単一細胞細胞成分の発現データセット中の各それぞれの第3の単一細胞細胞成分の発現データセットが、前記複数の非摂動細胞中の対応する単一細胞から得られ、
前記第4の複数の第4の単一細胞細胞成分の発現データセット中の各それぞれの第4の単一細胞細胞成分の発現データセットが、摂動に曝された前記複数の摂動細胞中の第4の複数の細胞の対応する単一細胞から得られる、項目1に記載の方法。
項目9: 前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをフィルタリングして、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーに含まれる細胞成分の数を減少させる工程をさらに含み、
任意に、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをフィルタリングする工程が、閾値p値にしたがって、又は細胞成分の閾値数にしたがって、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーに含まれる細胞成分の数を減少させる工程を含む、項目1に記載の方法。
項目10: それぞれの細胞成分についての前記対応する第2の重要度スコアを決定する工程が、
前記複数の細胞成分中の各それぞれの細胞成分について、前記それぞれの細胞成分についての前記重要度スコアを、前記それぞれの細胞成分についての対応するマッチングスコアで置き換える工程と;
前記複数の細胞成分についての前記マッチングスコアを組み合わせて、前記摂動についてのマッチングスコアを生成する工程と;
前記摂動が、前記それぞれの摂動についての前記マッチングスコアに基づいて、前記第1の細胞状態と前記変化した細胞状態との間の細胞の移行に関連しているかどうかを決定する工程とを含み、
任意に、前記対応するマッチングスコアが、別個の又は連続したスコアを含む、項目1に記載の方法。
項目11: 前記重要度スコアを置き換える工程が、
前記それぞれの細胞成分についての前記単一細胞の移行シグネチャーからの前記細胞成分量及び前記それぞれの細胞成分についての前記摂動シグネチャーからの前記細胞成分量が両方とも上方制御される場合、前記重要度スコアを第1のスコアで置き換える工程と;
前記それぞれの細胞成分についての前記単一細胞の移行シグネチャーからの前記細胞成分量が上方制御され、前記それぞれの細胞成分についての前記摂動シグネチャーからの前記細胞成分量が下方制御される場合、前記重要度スコアを第2のスコアで置き換える工程と;
前記それぞれの細胞成分についての前記摂動シグネチャーからの前記細胞成分量が、有意に上方制御又は下方制御されない場合、前記重要度スコアを第3のスコアで置き換える工程とを含む、項目10に記載の方法。
項目12: 前記重要度スコアを置き換える工程が、
前記それぞれの細胞成分についての前記単一細胞の移行シグネチャーからの前記細胞成分量及び前記細胞成分についての前記摂動シグネチャーからの前記細胞成分量が両方とも下方制御される場合、前記重要度スコアを第1のスコアで置き換える工程と;
前記それぞれの細胞成分についての前記単一細胞の移行シグネチャーからの前記細胞成分量が下方制御され、前記細胞成分についての前記摂動シグネチャーからの前記細胞成分量が上方制御される場合、前記重要度スコアを第2のスコアで置き換える工程と;
前記細胞成分についての前記摂動シグネチャーからの前記細胞成分量が、有意に上方制御又は下方制御されない場合、前記重要度スコアを第3のスコアで置き換える工程とを含む、項目10に記載の方法。
項目13: 前記複数の細胞成分が、任意にRNAレベルで測定される複数の遺伝子を含む、項目1に記載の方法。
項目14: 前記第1の複数の第1の単一細胞細胞成分の発現データセット及び前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各単一細胞細胞成分の発現データセットが、単一細胞リボ核酸(RNA)シーケンシング(scRNA-seq)、scTag-seq、シーケンシングを用いたトランスポゼースでアクセス可能なクロマチンの単一細胞アッセイ(scATAC-seq)、CyTOF/SCoP、E-MS/Abseq、miRNA-seq、CITE-seq、及びそれらの任意の組合せ、並びに単一細胞細胞成分の発現データセットにおける活性化された経路を表す線形組合せなどの組合せを含む、上記のものの集合からなる群から選択される方法を用いて生成される、項目2に記載の方法。
項目15: 前記比較に基づいて、前記変化した細胞状態を促進するものとして前記摂動を同定する工程、又は
前記比較に基づいて、前記変化した細胞状態を阻害するものとして前記摂動を同定する工程をさらに含む、項目1に記載の方法。
項目16: 前記細胞の移行シグネチャー及び前記摂動シグネチャーが、異なるタイプの細胞成分を用いて生成される、項目1に記載の方法。
項目17: 前記細胞の移行シグネチャー及び前記摂動シグネチャーが、同じタイプの細胞成分を用いて生成される、項目1に記載の方法。
項目18: 前記電子的形態でアクセスする工程が、複数の摂動における各それぞれのシグネチャーについて行われ、それによって、複数の摂動シグネチャーを得て、
前記比較する工程が、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーを、複数の摂動シグネチャーにおける各それぞれのシグネチャーと比較し、それによって、前記第1の細胞状態と前記変化した細胞状態との間の細胞の移行に関連する前記複数の摂動のサブセットを決定する、項目1に記載の方法。
項目19: 1つ以上のプロセッサ及びメモリを含むコンピュータシステムであって、前記メモリが、摂動が細胞の移行に影響を与えるかどうかを予測するための方法を実行するための命令を保存し、前記方法が、
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーに電子的形態でアクセスする工程であって、前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、前記単一細胞の移行シグネチャーは、複数の細胞成分の同定及び、前記複数の細胞成分中の各それぞれの細胞成分について、前記それぞれの細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化との間の関連性を定量する対応する第1の重要度スコアを含む、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーに電子的形態でアクセスする工程であって、前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定及び、前記複数の細胞成分の全て又は一部の中の各それぞれの細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間のそれぞれの細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化との間の関連性を定量する対応する第2の重要度スコアを含む、工程と;
前記単一細胞の移行シグネチャー及び前記摂動シグネチャーを比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程とを含む、コンピュータシステム。
項目20: 摂動が細胞の移行に影響を与えるかどうかを予測するための、コンピュータによって実行可能な1つ以上のコンピュータプログラムを保存する非一時的なコンピュータ可読媒体であって、前記コンピュータが、1つ以上のプロセッサ及びメモリを含み、前記1つ以上のコンピュータプログラムが、
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーに電子的形態でアクセスする工程であって、前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、前記単一細胞の移行シグネチャーは、複数の細胞成分の同定及び、前記複数の細胞成分中の各それぞれの細胞成分について、前記それぞれの細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化との間の関連性を定量する対応する第1の重要度スコアを含む、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーに電子的形態でアクセスする工程であって、前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定及び、前記複数の細胞成分の全て又は一部の中の各それぞれの細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間のそれぞれの細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化との間の関連性を定量する対応する第2の重要度スコアを含む、工程と;
前記単一細胞の移行シグネチャー及び前記摂動シグネチャーを比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程と
を含む方法を実行するコンピュータで実行可能な命令を集合的に符号化する、非一時的なコンピュータ可読媒体。
Claims (13)
- 摂動が細胞の移行に影響を与えるかどうかを予測するための方法であって、前記方法が、
メモリ及び1つ以上のプロセッサを含むコンピュータシステムにおいて:
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーにアクセスする工程であって、
前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、
前記単一細胞の移行シグネチャーは、複数の細胞成分の同定を含み、
前記単一細胞の移行シグネチャーは、前記複数の細胞成分中の各細胞成分について、前記細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性を定量化する第1の重要度スコアをさらに含み、
前記単一細胞の移行シグネチャーにアクセスする工程が、(i)第1の複数の第1の単一細胞細胞成分の発現データセット、及び(ii)第2の複数の第2の単一細胞細胞成分の発現データセットを得る工程を含み、ここで:
前記第1の複数の第1の単一細胞細胞成分の発現データセット中の各それぞれの第1の単一細胞細胞成分の発現データセットが、前記第1の細胞状態の複数の第1の細胞のうちの単一の第1の細胞から得られ、
前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各それぞれの第2の単一細胞細胞成分の発現データセットが、前記変化した細胞状態の複数の第2の細胞のうちの単一の第2の細胞から得られ、
前記細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性は、前記第1の複数の第1の単一細胞細胞成分の発現データセットと前記第2の複数の第2の単一細胞細胞成分の発現データセットのそれぞれの細胞成分の発現の差である、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーにアクセスする工程であって、
前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定を含み、
前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の中の各細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化と、の間の関連性を定量化する第2の重要度スコアをさらに含み、
細胞成分の発現の差の尺度が、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び一般化線形モデルのうちの1つを用いて、(i)第3の複数の第3の単一細胞細胞成分の発現データセットと、(ii)第4の複数の第4の単一細胞細胞成分の発現データセットとの間の細胞成分量の差を定量化し、ここで:
前記第3の複数の第3の単一細胞細胞成分の発現データセット中の各それぞれの第3の単一細胞細胞成分の発現データセットが、前記複数の非摂動細胞中の単一細胞から得られ、
前記第4の複数の第4の単一細胞細胞成分の発現データセット中の各それぞれの第4の単一細胞細胞成分の発現データセットが、摂動に曝された前記複数の摂動細胞中の単一細胞から得られる、工程と;
前記単一細胞の移行シグネチャーの第1の重要度スコアを、前記摂動シグネチャーの第2の重要度スコアと比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程と、
を含む、方法。 - 前記第1の複数の第1の単一細胞細胞成分の発現データセットの各それぞれのデータセットが、第1の複数の細胞成分ベクトル中の細胞成分ベクトルを含み、
前記第2の複数の第2の単一細胞細胞成分の発現データセットの各それぞれのデータセットが、第2の複数の細胞成分ベクトル中の細胞成分ベクトルを含み、
前記第1及び第2の複数の細胞成分ベクトル中の各細胞成分ベクトルが、複数の要素を含み、前記細胞成分ベクトル中の各要素が、前記複数の細胞成分中の細胞成分に関連しており、前記第1及び第2の複数の単一細胞細胞成分の発現データセットの前記それぞれのデータセットによって表される前記単一細胞についての対応する細胞成分の量を表す対応する値を含む、請求項1に記載の方法。 - 前記第1の複数の第1の単一細胞細胞成分の発現データセット及び/又は前記第2の複数の第2の単一細胞細胞成分の発現データセットに対して次元削減を行って、複数の次元削減成分を生成する工程と;
前記第1及び第2の複数の細胞成分ベクトル中の各細胞成分ベクトルについて、前記複数の次元削減成分を、前記細胞成分ベクトルに適用して、前記複数の次元削減成分中の各次元削減成分についての次元削減成分値を含む次元削減ベクトルを形成し、それによって、第1及び第2の複数の次元削減ベクトルを形成する工程と;
クラスタリングを行って、クラスターのセットCjを生成する工程であって、各クラスターが、前記第1及び第2の複数の次元削減ベクトルのサブセットに対応する複数の点を含む、工程と;
前記クラスターのセットCjの第1のクラスターから前記第1の複数の細胞を同定する工程と;
前記クラスターのセットCjの第2のクラスターから前記第2の複数の細胞を同定する工程とをさらに含む、請求項2に記載の方法。 - 前記複数の非摂動細胞が、摂動に曝されていない対照細胞である、請求項1に記載の方法。
- 前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをプルーニングして、前記複数の細胞成分を、RNAレベルで測定される転写因子に制限する工程をさらに含む、請求項1に記載の方法。
- 前記単一細胞の移行シグネチャーを決定する工程が、
母平均の差の検定、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び
一般化線形モデルのうちの1つを用いて、(i)前記第1の複数の第1の単一細胞細胞成分の発現データセットと、(ii)前記第2の複数の第2の単一細胞細胞成分の発現データセットと、の間の前記複数の細胞成分にわたる細胞成分量の差を決定する工程を含む、請求項1に記載の方法。 - 前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをフィルタリングして、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーに含まれる細胞成分の数を減少させる工程をさらに含み、
前記単一細胞の移行シグネチャー及び前記摂動シグネチャーをフィルタリングする工程が、閾値p値にしたがって、又は細胞成分の閾値数にしたがって、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーに含まれる細胞成分の数を減少させる工程を含む、請求項1に記載の方法。 - 前記複数の細胞成分が、RNAレベルで測定される複数の遺伝子を含む、請求項1に記
載の方法。 - 前記第1の複数の第1の単一細胞細胞成分の発現データセット及び前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各単一細胞細胞成分の発現データセットが、単一細胞リボ核酸(RNA)シーケンシング(scRNA-seq)、scTag-seq、シーケンシングを用いたトランスポゼースでアクセス可能なクロマチンの単一細胞アッセイ(scATAC-seq)、CyTOF、SCoPE-MS、Abseq、miRNA-seq、CITE-seq、及びそれらの任意の組合せ、並びに単一細胞細胞成分の発現データセットにおける活性化された経路を表す線形組合せを含む、上記のものの集合からなる群から選択される方法を用いて生成される、請求項1に記載の方法。
- 前記比較に基づいて、前記変化した細胞状態を促進するものとして前記摂動を同定する工程、又は
前記比較に基づいて、前記変化した細胞状態を阻害するものとして前記摂動を同定する工程をさらに含む、請求項1に記載の方法。 - 前記アクセスする工程が、複数の摂動における各それぞれのシグネチャーについて行われ、それによって、複数の摂動シグネチャーを得て、
前記比較する工程が、前記単一細胞の移行シグネチャー及び前記摂動シグネチャーを、複数の摂動シグネチャーにおける各それぞれのシグネチャーと比較し、それによって、前記第1の細胞状態と前記変化した細胞状態との間の細胞の移行に関連する前記複数の摂動のサブセットを決定する、請求項1に記載の方法。 - 1つ以上のプロセッサ及びメモリを含むコンピュータシステムであって、前記メモリが、摂動が細胞の移行に影響を与えるかどうかを予測するための方法を実行するための命令を保存し、前記方法が、
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーにアクセスする工程であって、
前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、
前記単一細胞の移行シグネチャーは、複数の細胞成分の同定を含み、
前記単一細胞の移行シグネチャーは、前記複数の細胞成分中の各細胞成分について、前記細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性を定量化する第1の重要度スコアをさらに含み、
前記単一細胞の移行シグネチャーにアクセスする工程が、(i)第1の複数の第1の単一細胞細胞成分の発現データセット、及び(ii)第2の複数の第2の単一細胞細胞成分の発現データセットを得る工程を含み、ここで:
前記第1の複数の第1の単一細胞細胞成分の発現データセット中の各それぞれの第1の単一細胞細胞成分の発現データセットが、前記第1の細胞状態の複数の第1の細胞のうちの単一の第1の細胞から得られ、
前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各それぞれの第2の単一細胞細胞成分の発現データセットが、前記変化した細胞状態の複数の第2の細胞のうちの単一の第2の細胞から得られ、
前記それぞれの細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性は、前記第1の複数の第1の単一細胞細胞成分の発現データセットと前記第2の複数の第2の単一細胞細胞成分の発現データセットのそれぞれの細胞成分の発現の差である、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーにアクセスする工程であって、前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定を含み、
前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の中の各細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化と、の間の関連性を定量化する第2の重要度スコアをさらに含み、
細胞成分の発現の差の尺度が、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び一般化線形モデルのうちの1つを用いて、(i)第3の複数の第3の単一細胞細胞成分の発現データセットと、(ii)第4の複数の第4の単一細胞細胞成分の発現データセットとの間の細胞成分量の差を定量化し、ここで:
前記第3の複数の第3の単一細胞細胞成分の発現データセット中の各それぞれの第3の単一細胞細胞成分の発現データセットが、前記複数の非摂動細胞中の単一細胞から得られ、
前記第4の複数の第4の単一細胞細胞成分の発現データセット中の各それぞれの第4の単一細胞細胞成分の発現データセットが、摂動に曝された前記複数の摂動細胞中の単一細胞から得られる、工程と;
前記単一細胞の移行シグネチャーの第1の重要度スコアを、前記摂動シグネチャーの第2の重要度スコアと比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程と、
を含む、コンピュータシステム。 - 摂動が細胞の移行に影響を与えるかどうかを予測するための、コンピュータによって実行可能な1つ以上のコンピュータプログラムを保存する非一時的なコンピュータ可読媒体であって、前記コンピュータが、1つ以上のプロセッサ及びメモリを含み、前記1つ以上のコンピュータプログラムが、
第1の細胞状態と変化した細胞状態との間の細胞成分の発現の差の尺度を表す単一細胞の移行シグネチャーにアクセスする工程であって、
前記変化した細胞状態は、前記第1の細胞状態から前記変化した細胞状態への細胞の移行を経て生じ、
前記単一細胞の移行シグネチャーは、複数の細胞成分の同定を含み、
前記単一細胞の移行シグネチャーは、前記複数の細胞成分中の各細胞成分について、前記細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性を定量化する第1の重要度スコアをさらに含み、
前記単一細胞の移行シグネチャーにアクセスする工程が、(i)第1の複数の第1の単一細胞細胞成分の発現データセット、及び(ii)第2の複数の第2の単一細胞細胞成分の発現データセットを得る工程を含み、ここで:
前記第1の複数の第1の単一細胞細胞成分の発現データセット中の各それぞれの第1の単一細胞細胞成分の発現データセットが、前記第1の細胞状態の複数の第1の細胞のうちの単一の第1の細胞から得られ、
前記第2の複数の第2の単一細胞細胞成分の発現データセット中の各それぞれの第2の単一細胞細胞成分の発現データセットが、前記変化した細胞状態の複数の第2の細胞のうちの単一の第2の細胞から得られ、
前記それぞれの細胞成分の発現の変化と、前記第1の細胞状態と前記変化した細胞状態との間の細胞状態の変化と、の間の関連性は、前記第1の複数の第1の単一細胞細胞成分の発現データセットと前記第2の複数の第2の単一細胞細胞成分の発現データセットのそれぞれの細胞成分の発現の差である、工程と;
複数の非摂動細胞と摂動に曝された複数の摂動細胞との間の細胞成分の発現の差の尺度を表す摂動シグネチャーにアクセスする工程であって、前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の同定を含み、
前記摂動シグネチャーは、前記複数の細胞成分の全て又は一部の中の各細胞成分について、(i)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞成分の発現の変化と、(ii)前記複数の非摂動細胞と前記複数の摂動細胞との間の細胞状態の変化と、の間の関連性を定量化する第2の重要度スコアをさらに含み、
細胞成分の発現の差の尺度が、ウィルコクソンの順位和検定、t検定、ロジスティック回帰、及び一般化線形モデルのうちの1つを用いて、(i)第3の複数の第3の単一細胞細胞成分の発現データセットと、(ii)第4の複数の第4の単一細胞細胞成分の発現データセットとの間の細胞成分量の差を定量化し、ここで:
前記第3の複数の第3の単一細胞細胞成分の発現データセット中の各それぞれの第3の単一細胞細胞成分の発現データセットが、前記複数の非摂動細胞中の単一細胞から得られ、
前記第4の複数の第4の単一細胞細胞成分の発現データセット中の各それぞれの第4の単一細胞細胞成分の発現データセットが、摂動に曝された前記複数の摂動細胞中の単一細胞から得られる、工程と;
前記単一細胞の移行シグネチャーの第1の重要度スコアを、前記摂動シグネチャーの第2の重要度スコアと比較し、それによって、前記摂動が前記細胞の移行に影響を与えるかどうかを決定する工程と、
を含む方法を実行するコンピュータで実行可能な命令を集合的に符号化する、非一時的なコンピュータ可読媒体。
Applications Claiming Priority (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201862698701P | 2018-07-16 | 2018-07-16 | |
| US62/698,701 | 2018-07-16 | ||
| US201962805884P | 2019-02-14 | 2019-02-14 | |
| US201962805888P | 2019-02-14 | 2019-02-14 | |
| US62/805,888 | 2019-02-14 | ||
| US62/805,884 | 2019-02-14 | ||
| PCT/US2019/041976 WO2020018519A1 (en) | 2018-07-16 | 2019-07-16 | Methods of analyzing cells |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022501011A JP2022501011A (ja) | 2022-01-06 |
| JP7548672B2 true JP7548672B2 (ja) | 2024-09-10 |
Family
ID=69138815
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021501001A Active JP7548672B2 (ja) | 2018-07-16 | 2019-07-16 | 細胞を分析する方法 |
Country Status (7)
| Country | Link |
|---|---|
| US (2) | US20200020419A1 (ja) |
| EP (1) | EP3824080A4 (ja) |
| JP (1) | JP7548672B2 (ja) |
| KR (1) | KR20210031708A (ja) |
| CN (1) | CN112424866B (ja) |
| CA (1) | CA3103677A1 (ja) |
| WO (1) | WO2020018519A1 (ja) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11715551B2 (en) * | 2019-06-19 | 2023-08-01 | Recursion Pharmaceuticals, Inc. | Systems and methods for evaluating query perturbations |
| US20210071256A1 (en) * | 2019-09-11 | 2021-03-11 | Recursion Pharmaceuticals, Inc. | Systems and methods for pairwise inference of drug-gene interaction networks |
| EP3901956A1 (en) * | 2020-04-21 | 2021-10-27 | ETH Zürich | Methods of determining correspondences between biological properties of cells |
| CN111551706B (zh) * | 2020-04-29 | 2023-04-07 | 成都微康生物科技有限公司 | 一种一次加样多项目联检的碟式免疫检测方法 |
| US12412100B2 (en) * | 2021-01-22 | 2025-09-09 | International Business Machines Corporation | Cell state transition features from single cell data |
| CN112967755B (zh) * | 2021-03-04 | 2023-06-13 | 深圳大学 | 一种面向单细胞rna测序数据的细胞类型识别方法 |
| WO2022266256A1 (en) | 2021-06-15 | 2022-12-22 | Flagship Pioneering Innovations Vi, Llc | Methods and systems for associating cellular constituents with a cellular process of interest |
| CA3222355A1 (en) | 2021-06-15 | 2022-12-22 | Fabian Alexander Wolf | Systems and methods for associating compounds with physiological conditions using fingerprint analysis |
| WO2022266259A1 (en) * | 2021-06-15 | 2022-12-22 | Flagship Pioneering Innovations Vi, Llc | Systems and methods for associating compounds with physiological conditions using fingerprint analysis |
| US20240361301A1 (en) * | 2021-09-09 | 2024-10-31 | Flagship Pioneering Innovations Vi, Llc | Methods and compositions for modulating goblet cells and for muco-obstructive diseases |
| EP4399278A1 (en) * | 2021-09-09 | 2024-07-17 | Flagship Pioneering Innovations VI, LLC | Methods and compositions for modulating enteroendocrine cells |
| WO2024025950A1 (en) * | 2022-07-26 | 2024-02-01 | The Regents Of The University Of California | Predicting gene expression alterations caused by perturbations of transcription factors |
| WO2024129916A1 (en) * | 2022-12-13 | 2024-06-20 | Cellarity, Inc. | Systems and methods for predicting compounds associated with transcriptional signatures |
| WO2025064862A1 (en) * | 2023-09-22 | 2025-03-27 | Verily Life Sciences Llc | Modules for single-cell data conversion to transformed bulk-cell data and methods of using the same |
| WO2025122843A1 (en) * | 2023-12-06 | 2025-06-12 | Cellino Biotech, Inc. | Systems and methods for analyzing and managing cell culture processes |
| US12119090B1 (en) * | 2023-12-19 | 2024-10-15 | Recursion Pharmaceuticals, Inc. | Utilizing masked autoencoder generative models to extract microscopy representation autoencoder embeddings |
| US12119091B1 (en) | 2023-12-19 | 2024-10-15 | Recursion Pharmaceuticals, Inc. | Utilizing masked autoencoder generative models to extract microscopy representation autoencoder embeddings |
| CN119936411A (zh) * | 2025-03-04 | 2025-05-06 | 南京傲因生物科技有限公司 | 基于单细胞蛋白修饰的标志物及其检测分析方法和应用 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002037102A2 (en) | 2000-10-20 | 2002-05-10 | Children's Medical Center Corporation | Methods for analyzing dynamic changes in cellular informatics |
| JP2002518021A (ja) | 1998-06-19 | 2002-06-25 | ロゼッタ・インファーマティクス・インコーポレーテッド | 生物学的ネットワークモデルのテスト方法 |
| WO2017075294A1 (en) | 2015-10-28 | 2017-05-04 | The Board Institute Inc. | Assays for massively combinatorial perturbation profiling and cellular circuit reconstruction |
Family Cites Families (70)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040259073A1 (en) | 2003-03-07 | 2004-12-23 | Ariang Hassibi | Charge perturbation signature methods and devices for membrane analysis |
| WO2005107412A2 (en) | 2004-04-30 | 2005-11-17 | Rosetta Inpharmatics Llc | Systems and methods for reconstruction gene networks in segregating populations |
| US20070238094A1 (en) | 2005-12-09 | 2007-10-11 | Baylor Research Institute | Diagnosis, prognosis and monitoring of disease progression of systemic lupus erythematosus through blood leukocyte microarray analysis |
| EP2044431B1 (en) * | 2006-07-17 | 2019-04-24 | H. Lee Moffitt Cancer Center & Research Institute, Inc. | Computer systems and methods for selecting subjects for clinical trials |
| US9458495B2 (en) | 2008-07-01 | 2016-10-04 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and systems for assessment of clinical infertility |
| ES2560209T3 (es) * | 2009-01-20 | 2016-02-17 | The Board Of Trustees Of The Leland Stanford Junior University | Expresión genética en células individuales para el diagnóstico, pronóstico e identificación de dianas farmacológicas |
| US9315857B2 (en) | 2009-12-15 | 2016-04-19 | Cellular Research, Inc. | Digital counting of individual molecules by stochastic attachment of diverse label-tags |
| WO2012061510A2 (en) | 2010-11-03 | 2012-05-10 | Merck Sharp & Dohme Corp. | Methods of predicting cancer cell response to therapeutic agents |
| EP2675819B1 (en) | 2011-02-18 | 2020-04-08 | Bio-Rad Laboratories, Inc. | Compositions and methods for molecular labeling |
| WO2013009705A2 (en) | 2011-07-09 | 2013-01-17 | The Trustees Of Columbia University In The City Of New York | Biomarkers, methods, and compositions for inhibiting a multi-cancer mesenchymal transition mechanism |
| US9953129B2 (en) | 2011-09-23 | 2018-04-24 | Agency For Science, Technology And Research | Patient stratification and determining clinical outcome for cancer patients |
| US20150154352A1 (en) | 2012-06-21 | 2015-06-04 | Gigagen, Inc. | System and Methods for Genetic Analysis of Mixed Cell Populations |
| CA2879426A1 (en) * | 2013-01-14 | 2014-07-17 | Cellecta, Inc. | Methods and compositions for single cell expression profiling |
| US20160303133A1 (en) | 2013-04-05 | 2016-10-20 | Numedii, Inc. | Treatment of gastrointestinal and other disorders |
| JP6600302B2 (ja) * | 2013-06-12 | 2019-10-30 | ザ ジェネラル ホスピタル コーポレイション | 標的分子のマルチプレックス検出のための方法、キット、およびシステム、ならびにそれらの使用 |
| US10184933B2 (en) | 2013-10-01 | 2019-01-22 | The J. David Gladstone Industries | Compositions, systems and methods for gene expression noise drug screening and uses thereof |
| US10801070B2 (en) | 2013-11-25 | 2020-10-13 | The Broad Institute, Inc. | Compositions and methods for diagnosing, evaluating and treating cancer |
| JP6306370B2 (ja) | 2014-02-25 | 2018-04-04 | 千代田化工建設株式会社 | 芳香族化合物の水素化システムおよび水素化方法 |
| WO2015136516A1 (en) | 2014-03-13 | 2015-09-17 | Pathway Pharmaceuticals Ltd | System, method and software for analysis of intracellular signaling pathway activation using transcriptomic data |
| CN107873054B (zh) | 2014-09-09 | 2022-07-12 | 博德研究所 | 用于复合单细胞核酸分析的基于微滴的方法和设备 |
| US20170362654A1 (en) | 2014-12-09 | 2017-12-21 | The Trustees Of Princeton University | Biomarkers of oocyte quality |
| WO2016115537A2 (en) | 2015-01-15 | 2016-07-21 | Massachusetts Institute Of Technology | Systems, methods, and apparatus for in vitro single-cell identification and recovery |
| US10453555B2 (en) | 2015-01-16 | 2019-10-22 | University Of Virginia Patent Foundation | Parameterizing cell-to-cell regulatory heterogeneities via stochastic transcriptional profiles |
| WO2016191533A1 (en) | 2015-05-26 | 2016-12-01 | The Trustees Of Columbia University In The City Of New York | Rna printing and sequencing devices, methods, and systems |
| US10745741B2 (en) | 2015-06-26 | 2020-08-18 | European Molecular Biology Laboratory | Cell barcoding in microfluidics |
| EP3314020A1 (en) | 2015-06-29 | 2018-05-02 | The Broad Institute Inc. | Tumor and microenvironment gene expression, compositions of matter and methods of use thereof |
| WO2017040315A1 (en) | 2015-08-28 | 2017-03-09 | The Trustees Of Columbia University In The City Of New York | Virtual inference of protein activity by regulon enrichment analysis |
| CN108350510B (zh) | 2015-09-09 | 2022-06-03 | 普梭梅根公司 | 用于胃肠健康相关病症的源自微生物群系的诊断及治疗方法和系统 |
| US11996168B2 (en) | 2015-10-28 | 2024-05-28 | The Broad Institute, Inc. | Systems and methods for determining relative abundances of biomolecules |
| WO2017083817A1 (en) | 2015-11-13 | 2017-05-18 | Arizona Board Of Regents On Behalf Of Arizona State University | A device and method for high-throughput multiparameter measurements in one or more live and fixed cells |
| US20170154163A1 (en) | 2015-12-01 | 2017-06-01 | Ramot At Tel-Aviv University Ltd. | Clinically relevant synthetic lethality based method and system for cancer prognosis and therapy |
| US20170260584A1 (en) | 2016-02-11 | 2017-09-14 | 10X Genomics, Inc. | Cell population analysis using single nucleotide polymorphisms from single cell transcriptomes |
| US11127484B2 (en) | 2016-03-17 | 2021-09-21 | Regents Of The University Of Minnesota | Analysis of single cell transcriptomics |
| US11427861B2 (en) | 2016-03-17 | 2022-08-30 | Massachusetts Institute Of Technology | Methods for identifying and modulating co-occurant cellular phenotypes |
| US11262358B2 (en) | 2016-04-06 | 2022-03-01 | Technion Research & Development Foundation Limited | Infiltrating immune cell proportions predict anti-TNF response in colon biopsies |
| WO2017178345A1 (en) | 2016-04-15 | 2017-10-19 | Koninklijke Philips N.V. | Sub-population detection and quantization of receptor-ligand states for characterizing inter-cellular communication and intratumoral heterogeneity |
| CN109074430B (zh) | 2016-05-26 | 2022-03-29 | 贝克顿迪金森公司 | 分子标记计数调整方法 |
| US11913017B2 (en) | 2016-06-29 | 2024-02-27 | The Regents Of The University Of California | Efficient genetic screening method |
| US11225689B2 (en) | 2016-08-17 | 2022-01-18 | The Broad Institute, Inc. | Method for determination and identification of cell signatures and cell markers |
| AU2017331459B2 (en) | 2016-09-26 | 2023-04-13 | Becton, Dickinson And Company | Measurement of protein expression using reagents with barcoded oligonucleotide sequences |
| EP3519554A4 (en) | 2016-09-30 | 2020-05-20 | The Governing Council Of The University Of Toronto | System for identifying and targeting individual cells within a heterogeneous population for selective extraction of cellular content |
| CN106777870A (zh) | 2016-11-18 | 2017-05-31 | 邹欣 | 一种针对单细胞转录组数据的降噪声算法 |
| US12300357B2 (en) | 2016-12-14 | 2025-05-13 | FlowJo, LLC | Applied computer technology for management, synthesis, visualization, and exploration of parameters in large multi-parameter data sets |
| WO2018132635A1 (en) | 2017-01-12 | 2018-07-19 | Massachusetts Institute Of Technology | Methods for analyzing t cell receptors and b cell receptors |
| AU2018210695B2 (en) | 2017-01-20 | 2024-07-18 | The University Of British Columbia | Molecular subtyping, prognosis, and treatment of bladder cancer |
| WO2018140302A1 (en) | 2017-01-24 | 2018-08-02 | The Regents Of The University Of Michigan | Systems and methods for whole cell analysis |
| US10347365B2 (en) | 2017-02-08 | 2019-07-09 | 10X Genomics, Inc. | Systems and methods for visualizing a pattern in a dataset |
| GB201702847D0 (en) | 2017-02-22 | 2017-04-05 | Cancer Res Tech Ltd | Cell labelling, tracking and retrieval |
| US11963966B2 (en) | 2017-03-31 | 2024-04-23 | Dana-Farber Cancer Institute, Inc. | Compositions and methods for treating ovarian tumors |
| US20200071773A1 (en) | 2017-04-12 | 2020-03-05 | Massachusetts Eye And Ear Infirmary | Tumor signature for metastasis, compositions of matter methods of use thereof |
| WO2018204854A1 (en) | 2017-05-05 | 2018-11-08 | The United States Of America, As Represented By The Secretary, Department Of Health And Human Services | Methods of preparing a re-usable single cell and methods for analyzing the epigenome, transcriptome, and genome of a single cell |
| WO2018209324A2 (en) | 2017-05-11 | 2018-11-15 | The Broad Institute, Inc. | Methods and compositions of use of cd8+ tumor infiltrating lymphocyte subtypes and gene signatures thereof |
| CA3062985A1 (en) | 2017-05-16 | 2018-11-22 | The Chinese University Of Hong Kong | Integrative single-cell and cell-free plasma rna analysis |
| CN110870018B (zh) | 2017-05-19 | 2024-11-22 | 10X基因组学有限公司 | 用于分析数据集的系统和方法 |
| EP3631417B1 (en) | 2017-05-25 | 2024-02-14 | Flowjo, LLC | Visualization, comparative analysis, and automated difference detection for large multi-parameter data sets |
| WO2018232195A1 (en) | 2017-06-14 | 2018-12-20 | The Broad Institute, Inc. | Compositions and methods targeting complement component 3 for inhibiting tumor growth |
| WO2018236852A1 (en) | 2017-06-19 | 2018-12-27 | Jungla Inc. | Interpretation of genetic and genomic variants via an integrated computational and experimental deep mutational learning framework |
| US20200131564A1 (en) | 2017-07-07 | 2020-04-30 | Board Of Regents, The University Of Texas System | High-coverage and ultra-accurate immune repertoire sequencing using molecular identifiers |
| US12049643B2 (en) | 2017-07-14 | 2024-07-30 | The Broad Institute, Inc. | Methods and compositions for modulating cytotoxic lymphocyte activity |
| EP3655961A4 (en) | 2017-07-17 | 2021-09-01 | Massachusetts Institute of Technology | HEALTHY AND SICK BARRIER TISSUE CELL ATLAS |
| JP7545891B2 (ja) | 2017-07-21 | 2024-09-05 | ザ ボード オブ トラスティーズ オブ ザ レランド スタンフォード ジュニア ユニバーシティー | 混合細胞集団を分析するためのシステムおよび方法 |
| CN107451424A (zh) | 2017-07-31 | 2017-12-08 | 浙江绍兴千寻生物科技有限公司 | 大批量单细胞RNA‑seq数据质量控制和分析方法 |
| US11981957B2 (en) | 2017-08-21 | 2024-05-14 | The Johns Hopkins University | Methodology to identify biomarkers relevant to neurons in the brain by using non-invasive nasal biopsy |
| WO2019079362A1 (en) | 2017-10-16 | 2019-04-25 | Massachusetts Institute Of Technology | HOST-PATHOGENIC INTERACTION OF MYCOBACTERIUM TUBERCULOSIS |
| WO2019084046A1 (en) | 2017-10-23 | 2019-05-02 | The Broad Institute, Inc. | ENRICHING SINGLE CELL CELLULAR CONSTITUENT OF BARCODE SEQUENCING LIBRARY |
| US20190142722A1 (en) | 2017-11-10 | 2019-05-16 | The Trustees Of Columbia University In The City Of New York | Methods and compositions for promoting or inducing hair growth |
| US11994512B2 (en) | 2018-01-04 | 2024-05-28 | Massachusetts Institute Of Technology | Single-cell genomic methods to generate ex vivo cell systems that recapitulate in vivo biology with improved fidelity |
| CN109033743B (zh) | 2018-07-25 | 2021-01-01 | 上海交通大学 | 一种降低单细胞转录组数据中技术噪声的方法 |
| CN109486750A (zh) | 2018-10-16 | 2019-03-19 | 中国农业科学院北京畜牧兽医研究所 | 用于单细胞rna测序的羊卵母细胞和卵丘颗粒细胞采集方法 |
| CN109628555A (zh) | 2019-01-30 | 2019-04-16 | 杭州瑞普基因科技有限公司 | 适用于单细胞测序的人冰冻肿瘤组织细胞核分离方法 |
-
2019
- 2019-07-15 US US16/511,691 patent/US20200020419A1/en not_active Abandoned
- 2019-07-16 JP JP2021501001A patent/JP7548672B2/ja active Active
- 2019-07-16 WO PCT/US2019/041976 patent/WO2020018519A1/en not_active Ceased
- 2019-07-16 CN CN201980047446.0A patent/CN112424866B/zh active Active
- 2019-07-16 CA CA3103677A patent/CA3103677A1/en active Pending
- 2019-07-16 EP EP19838246.7A patent/EP3824080A4/en not_active Withdrawn
- 2019-07-16 KR KR1020217003121A patent/KR20210031708A/ko not_active Ceased
-
2023
- 2023-10-11 US US18/485,119 patent/US12282499B1/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002518021A (ja) | 1998-06-19 | 2002-06-25 | ロゼッタ・インファーマティクス・インコーポレーテッド | 生物学的ネットワークモデルのテスト方法 |
| WO2002037102A2 (en) | 2000-10-20 | 2002-05-10 | Children's Medical Center Corporation | Methods for analyzing dynamic changes in cellular informatics |
| WO2017075294A1 (en) | 2015-10-28 | 2017-05-04 | The Board Institute Inc. | Assays for massively combinatorial perturbation profiling and cellular circuit reconstruction |
Non-Patent Citations (3)
| Title |
|---|
| Current Opinion in Systems Biology,2017年12月28日,7, pp.36-46 |
| Proceedings of the National Academy of Sciences of the United States of America,2017年,114(9),2271-2276, Supporting Information |
| Science,2006年,313(5795),1929-1935 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20200020419A1 (en) | 2020-01-16 |
| EP3824080A4 (en) | 2022-05-04 |
| CN112424866A (zh) | 2021-02-26 |
| EP3824080A1 (en) | 2021-05-26 |
| JP2022501011A (ja) | 2022-01-06 |
| CN112424866B (zh) | 2025-01-07 |
| US20250124057A1 (en) | 2025-04-17 |
| WO2020018519A1 (en) | 2020-01-23 |
| CA3103677A1 (en) | 2020-01-23 |
| KR20210031708A (ko) | 2021-03-22 |
| US12282499B1 (en) | 2025-04-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112424866B (zh) | 分析细胞的方法 | |
| Bakken et al. | Comparative cellular analysis of motor cortex in human, marmoset and mouse | |
| Klein et al. | Mapping cells through time and space with moscot | |
| US8442772B2 (en) | Compositions and methods for defining cells | |
| Stanfield et al. | Myometrial transcriptional signatures of human parturition | |
| Ellwanger et al. | Transcriptional dynamics of hair-bundle morphogenesis revealed with CellTrails | |
| Alli Shaik et al. | Functional mapping of the zebrafish early embryo proteome and transcriptome | |
| KR20180086438A (ko) | 다양한 화학 구조를 가지는 생체이물질의 세포 유형 특이적 독성을 예측하기 위한 고처리량 영상화 기반 방법 | |
| JP2024525155A (ja) | フィンガープリント分析を使用して化合物を生理学的状態と関連付けるためのシステム及び方法 | |
| WO2020257501A1 (en) | Systems and methods for evaluating query perturbations | |
| EP4029019A1 (en) | Systems and methods for pairwise inference of drug-gene interaction networks | |
| Ansari et al. | A novel pathway analysis approach based on the unexplained disregulation of genes | |
| Que et al. | Transcriptomic homogeneity and an age-dependent onset of hemoglobin expression characterize morphological PV types | |
| WO2022266257A1 (en) | Systems and methods for associating compounds with properties using clique analysis of cell-based data | |
| US12374430B2 (en) | Systems and methods for associating compounds with properties using clique analysis of cell-based data | |
| He et al. | A statistical method for quantifying progenitor population reveals incipient cell fate commitments | |
| Ramirez et al. | Modeling a gene regulatory network of EMT hybrid states for mouse embryonic skin cells | |
| Yang et al. | A Manifold-Based Measure of Transcriptional Entropy for Quantifying Aging in Single Cells | |
| Patili et al. | imAgeScore, a Cell Painting-Based Predictor of Cellular Age for High-throughput Drug Screening Applications | |
| Crow et al. | Addressing the looming identity crisis in single cell RNA-seq | |
| Baek et al. | Comparison of transcriptomic and phenomic profiles for the prediction of drug mechanism | |
| Hans et al. | High-dimensional regression in cancer genomics | |
| Shpak et al. | Building an atlas of mechanobiology: high-throughput contractility screen of 2418 kinase inhibitors in five primary human cell types reveals selective divergent responses among related cell types | |
| Kunz et al. | DBSOMA: A Machine Learning Method that Identifies Chemical Modulators of Transcriptional States Uncovers Effectors of Beta-Cell Maturation | |
| Bondhus | A Computational Approach to Exploring the Role of Chromatin Modifiers in Development and Disease |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20211126 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220708 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230616 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230627 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20230926 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20231124 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231227 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240206 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240507 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240730 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240827 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7548672 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |