JP2020536530A - 標的遺伝子発現の数学的モデル化を使用する、Notch細胞シグナル伝達経路活性の評価 - Google Patents

標的遺伝子発現の数学的モデル化を使用する、Notch細胞シグナル伝達経路活性の評価 Download PDF

Info

Publication number
JP2020536530A
JP2020536530A JP2020518608A JP2020518608A JP2020536530A JP 2020536530 A JP2020536530 A JP 2020536530A JP 2020518608 A JP2020518608 A JP 2020518608A JP 2020518608 A JP2020518608 A JP 2020518608A JP 2020536530 A JP2020536530 A JP 2020536530A
Authority
JP
Japan
Prior art keywords
notch
activity
target genes
cell signaling
signaling pathway
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.)
Withdrawn
Application number
JP2020518608A
Other languages
English (en)
Inventor
デ ストルプ アニア ファン
デ ストルプ アニア ファン
ローレンティウス ヘンリクス フランシスクス マリア ホルツァー
ローレンティウス ヘンリクス フランシスクス マリア ホルツァー
ウィルヘルムス フランシスクス ヨハネス フェルハーフ
ウィルヘルムス フランシスクス ヨハネス フェルハーフ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2020536530A publication Critical patent/JP2020536530A/ja
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6809Methods for determination or identification of nucleic acids involving differential detection
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/686Polymerase chain reaction [PCR]
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • C12Q1/6886Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/68Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B5/00ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/158Expression markers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/50Determining the risk of developing a disease
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/52Predicting or monitoring the response to treatment, e.g. for selection of therapy based on assay results in personalised medicine; Prognosis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/56Staging of a disease; Further complications associated with the disease
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Physics & Mathematics (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biotechnology (AREA)
  • Genetics & Genomics (AREA)
  • Biophysics (AREA)
  • Biochemistry (AREA)
  • Microbiology (AREA)
  • General Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Hematology (AREA)
  • Urology & Nephrology (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Data Mining & Analysis (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Cell Biology (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)

Abstract

本発明は、対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルに基づく、対象におけるNotch細胞シグナル伝達経路の活性を推測するための、コンピュータにより実施される方法に関する。本発明は、対象におけるNotch細胞シグナル伝達経路の活性を推測するための装置、非一時的記憶媒体、及びコンピュータプログラムにさらに関する。本発明は、対象の試料中のNotch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを測定するためのキット、対象におけるNotch細胞シグナル伝達経路の活性を推測するためのキット、及び方法の実施におけるこのようなキットの使用にさらに関する。

Description

本発明は、一般に、バイオインフォマティクス、ゲノムデータ処理、プロテオームデータ処理、及び、これらの関連技術の分野に関する。より詳細には、本発明は、デジタル処理デバイスにより実施される、対象におけるNotch細胞シグナル伝達経路の活性を推測するための、コンピュータにより実施される方法であって、コンピュータにより実施される方法は、推測することが、対象の試料中で測定されるNotch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルに基づく、コンピュータにより実施される方法に関する。本発明は、対象におけるNotch細胞シグナル伝達経路の活性を推測するための装置であって、装置は、方法を実施するように構成されたデジタルプロセッサを含む装置、対象におけるNotch細胞シグナル伝達経路の活性を推測するための非一時的記憶媒体であって、非一時的記憶媒体は、上記方法を実施するための、デジタル処理デバイスにより実行可能な命令を保存する非一時的記憶媒体、及び対象におけるNotch細胞シグナル伝達経路の活性を推測するためのコンピュータプログラムであって、コンピュータプログラムは、コンピュータプログラムがデジタル処理デバイス上で実行される場合に、デジタル処理デバイスに、方法を実施させるためのプログラムコード手段を含むコンピュータプログラムにさらに関する。本発明は、対象の試料中のNotch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを測定するためのキット、対象におけるNotch細胞シグナル伝達経路の活性を推測するためのキット、及び方法の実施におけるキットの使用にさらに関する。
ゲノム解析及びプロテオーム解析は、多様ながんが、がんの増殖及び進化、例えば、細胞の増殖及び転移において役割を果たす、特異的な遺伝子について、ゲノム突然変異/変異の特異的な組合せ、及び/又は高度若しくは低度の発現レベルと関連することが知られている腫瘍学など医療分野における臨床適用のためにかなり実現され、潜在的な可能性を有する。
Notchとは、胚の発生、免疫応答、及びがんに関与する多くの遺伝子の発現を調節する、誘導型転写因子である。がん(例えば、乳がん又は卵巣がん)などの病理学的障害に関して、異常なNotch経路活性は、重要な役割を果たす(Aster J.C.ら、「The varied roles of Notch in cancer」、Annual Review of Pathology、12巻、1号、2016年12月、245〜275頁を参照されたい)。Notch細胞シグナル伝達経路は、Notchファミリーに由来するタンパク質受容体と(細胞結合型)リガンドのファミリー(DSLファミリー)とからなり、リガンドのファミリーは結合型受容体の切断を誘導し、このとき、切断された細胞内断片が、核へと移行し、ここで、他のタンパク質と併せて、活性の転写因子複合体を形成し、これが、標的遺伝子の十分に規定されたセットに結合し、転写活性化させる。(また、Guruharsha K.G.ら、「The Notch signaling system:recent insights into the complexity of a conserved pathway」、Nature Reviews Genetics、13巻、2012年9月、654〜666頁に基づく、図1も参照されたい)。
ターゲティングされた薬による処置の正しい選択を可能とするためには、例えば、がんにおける、Notchシグナル伝達に関して、異常なNotchシグナル伝達活性を検出することが可能であることは重要である。現在、抗Notch療法が開発されている(Espinoza I.及びMiele L.、「Notch inhibitors for cancer treatment」、Pharmacology & Therapeutics、 139巻、2号、2013年8月、95〜110頁を参照されたい)。しかし現在、その活性の状態において、例えば、その不活性状態と比較して、腫瘍促進性である可能性が高いことを指し示す、Notch細胞シグナル伝達経路の機能状態及び活性のそれぞれについて評価するのに利用可能な臨床アッセイは存在しない。したがって、がん、例えば、乳がん、子宮頸がん、子宮内膜がん、卵巣がん、膵臓がん、若しくは前立腺がん、又は免疫障害などの疾患であって、異常なNotch細胞シグナル伝達経路の活性により、少なくとも部分的に駆動され、したがって、Notch細胞シグナル伝達経路の阻害剤に応答する可能性が高い疾患を有する患者を特徴付ける可能性を改善することが可能であることが所望される。
本発明の主要な態様に従い、上記の問題は、デジタル処理デバイスにより実施される、対象におけるNotch細胞シグナル伝達経路の活性を推測するための、コンピュータにより実施される方法であって、コンピュータにより実施される方法は、推測することが、
対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ以上の、例えば、3つ、4つ、5つ、6つ、7つ、8つ、9つ、10、11、12、以上の標的遺伝子の発現レベルを受信するステップと、
対象の試料中の3つ以上のNotch標的遺伝子の転写を制御する、Notch転写因子(TF)エレメントの活性レベルを決定するステップであって、決定するステップは、3つ以上のNotch標的遺伝子の発現レベルを、Notch TFエレメントの活性レベルと関連付ける較正された数学的経路モデルを査定することに基づく、決定するステップと、
対象におけるNotch細胞シグナル伝達経路の活性を、対象の試料中のNotch TFエレメントの決定された活性レベルに基づき推測するステップと
を有し、
3つ以上のNotch標的遺伝子が、CD28、CD44、DLGAP5、DTX1、EPHB3、FABP7、GFAP、GIMAP5、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、KLF5、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、PTCRA、SOX9、及びTNCからなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、NRARP、及びPTCRAからなる群から選択され、1つ又は複数の、例えば、2、3、4以上の、Notch標的遺伝子が、CD28、CD44、DLGAP5、EPHB3、FABP7、GFAP、GIMAP5、HES7、HEY1、HEYL、KLF5、NFKB2、NOX1、PBX1、PIN1、PLXND1、SOX9、及びTNCからなる群から選択される、
コンピュータにより実施される方法により解決される。
本明細書に、TFエレメントの「活性レベル」とは、その標的遺伝子の転写に関する、TFエレメントの活性のレベルを表示する。
本発明は、Notch細胞シグナル伝達経路内で生じる作用を同定する、適切な方式は、(とりわけ)Notch細胞シグナル伝達経路により制御される、Notch転写因子(TF)エレメントにより制御される、標的遺伝子の転写である、Notch細胞シグナル伝達経路シグナル伝達出力の測定に基づきうるという、本発明者らによる、新しい考え方に基づく。本発明者らによる、この新しい考え方は、TF活性レベルが、試料では、準定常状態にあり、これは、(とりわけ)Notch標的遺伝子の発現値により検出されうることを仮定する。本明細書で標的とされるNotch細胞シグナル伝達経路は、増殖、分化、及び創傷の治癒など、ヒトの、多くの細胞型において、多くの機能を制御することが知られている。がん(例えば、乳がん、子宮頸がん、子宮内膜がん、卵巣がん、膵臓がん、又は前立腺がん)などの病理学的障害に関して、異常なNotch細胞シグナル伝達活性は、重要な役割を果しており、これは標的遺伝子の発現プロファイル内で検出可能であり、したがって、較正された数学的経路モデルの手段として利用される。
本発明は、(i)対象の試料中の、Notch TFエレメントの活性レベルを決定するステップであって、決定するステップは、その転写が、Notch TFエレメントにより制御される、Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを、Notch TFエレメントの活性レベルと関連付ける、較正された数学的モデルを査定することに基づく、決定するステップと、(ii)対象におけるNotch細胞シグナル伝達経路の活性を、対象の試料中のNotch TFエレメントの決定された活性レベルに基づき推測するステップとにより、対象におけるNotch細胞シグナル伝達経路の活性を決定することを可能とする。これは、がん、例えば、乳がん、子宮頸がん、子宮内膜がん、卵巣がん、膵臓がん、又は前立腺がんなどの疾患であって、異常なNotch細胞シグナル伝達経路の活性により、少なくとも部分的に駆動され、したがって、Notch細胞シグナル伝達経路の阻害剤に応答する可能性が高い疾患を有する患者を特徴付ける可能性の改善を可能とすることが好ましい。特定の実施形態では、処置の決定は、特異的なNotch細胞シグナル伝達経路活性に基づきうる。特定の実施形態では、Notch細胞シグナル伝達状態は、Notch細胞シグナル伝達経路が活性であるオッズのカットオフ値、例えば、10:1、5:1、4:1、2:1、1:1、1:2、1:4、1:5、又は1:10に設定することができる。
本明細書では、「Notch転写因子エレメント」又は「Notch TFエレメント」又は「TFエレメント」という用語は、Notchタンパク質のうちの1つの、少なくとも細胞内ドメイン(対応する細胞内ドメインである、N1ICD、N2ICD、N3ICD、及びN4ICDを伴う、Notch1、Notch2、Notch3、及びNotch4)を、特異的なDNA配列への結合が可能な、DNA結合性転写因子CSL(CBF1/RBP−Jκ、Su(H)、及びLAG−1)、及び、好ましくは、転写を活性化させ、これにより、標的遺伝子の転写を制御するのに要求される、mastermind−like(MAML)ファミリー(MAML1、MAML2、及びMAML3)に由来する、1つの共活性化因子タンパク質などの補因子と共含有するタンパク質複合体であると規定される。好ましくは、用語は、Notch細胞内ドメイン(N1ICD、N2ICD、N3ICD、及びN4ICD)を結果としてもたらす、Notchタンパク質(Notch1、Notch2、Notch3、及びNotch4)のうちの1つの切断により誘発される、タンパク質又はタンパク質複合体である転写因子を指す。例えば、近傍の細胞上で発現する、DSLリガンド(DLL1、DLL3、DLL4、Jagged1、及びJagged2)は、Notchタンパク質/受容体の細胞外ドメインに結合し、細胞内のNotchシグナル伝達経路を誘発し、このNotch細胞内ドメインは、発現を制御するNotchシグナル伝達カスケードに参与することが知られている。
較正された数学的経路モデルは、Notch TFエレメントの活性レベルと、3つ若しくはこれを超えるNotch標的遺伝子の発現レベルとを関連付ける条件付き確率に基づく確率モデル、好ましくは、ベイジアンネットワークモデルであってもよく、又は3つ若しくはこれを超えるNotch標的遺伝子の発現レベルの1つ若しくは複数の線形結合に基づいてもよい。特に、Notch細胞シグナル伝達経路の活性の推測は、その内容が、それらの全体において本明細書に組み込まれる、国際特許出願公開第WO2013/011479A2号(「標的遺伝子発現の確率モデル化を使用する、細胞シグナル伝達経路活性の評価」)において開示される通りに、又は国際特許出願公開第WO2014/102668A2号(「標的遺伝子発現の線形結合を使用する、細胞シグナル伝達経路活性の評価」)において記載される通りに実施される。標的遺伝子発現の数学的モデル化を使用する、細胞シグナル伝達経路活性の推測に関するさらなる詳細は、Verhaegh W.ら、「Selection of personalized patient therapy through the use of knowledge−based computational models that identify tumor−driving signal transduction pathways」、Cancer Research、74巻、11号、2014年、2936〜2945頁において見出されうる。
本明細書で使用される「対象」という用語は、任意の生物を指す。一部の実施形態では、対象は、動物、好ましくは、哺乳動物である。ある特定の実施形態では、対象は、ヒト対象、好ましくは、医学的対象である。さらに他の実施形態では、対象は、細胞系である。
本明細書で使用される「標的遺伝子」という用語は、その転写が、Notch転写因子エレメントにより、直接的に、又は間接的に制御される遺伝子を意味する。「標的遺伝子」は、「直接的標的遺伝子」及び/又は「間接的標的遺伝子」(本明細書で記載される)である。さらに、「標的遺伝子」は、「直接的標的遺伝子」及び/又は「間接的標的遺伝子」(本明細書で記載される)である。
3つ以上のNotch標的遺伝子が、CD44、DTX1、EPHB3、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、CD44、EPHB3、HES7、HEY1、HEYL、NFKB2、NOX1、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択される方法(本明細書で記載される)が特に好ましい。
3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択される方法(本明細書で記載される)が特に好ましい。
特に適するNotch標的遺伝子については、以下の本文の節のほか、下記の例(例えば、下記の表1〜3を参照されたい)においても記載される。
したがって、好ましい実施形態に従い、Notch標的遺伝子は、下記の表1、表2、又は表3に列挙されるNotch標的遺伝子からなる群から選択される。
本発明者らにより、連鎖的な短いリストの中のNotch標的遺伝子が、ますます、Notch細胞シグナル伝達経路の活性を決定するための証拠となることが見出されている。
本発明の別の態様は、
Notch細胞シグナル伝達経路が対象において異常に作動しているのかどうかを、対象におけるNotch細胞シグナル伝達経路の推測された活性に基づき決定するステップ
をさらに有する方法(本明細書で記載される)に関する。
本発明はまた、
Notch細胞シグナル伝達経路の異常な作動を是正する薬物を、対象に処方することを推奨するステップ
をさらに有し、
推奨するステップは、Notch細胞シグナル伝達経路が、対象において異常に作動すると、Notch細胞シグナル伝達経路の推測された活性に基づき決定される場合に実施される、方法(本明細書で記載される)にも関する。
「細胞シグナル伝達経路は、作動が異常である」という語句は、経路の「活性」が、予期される通りではない場合を指し、この場合、「活性」という用語は、標的遺伝子を発現へと駆動するときの、転写因子複合体の活性、すなわち、標的遺伝子が転写される速度を指す。「正常」とは、不活性であることが予期される組織内で不活性であり、活性であることが予期される組織内で活性である場合である。さらに、「正常」と考えられる、ある特定のレベルの活性が存在し、これより高いか、又は低い、何らかのレベルは、「異常」であると考えられる。
本発明はまた、Notch細胞シグナル伝達経路の異常な作動が、対象において腫瘍プロモーターとしてNotch細胞シグナル伝達経路が作動するような作動である方法(本明細書で記載される)にも関する。
本発明に従い使用される試料は、抽出された試料、すなわち、対象から抽出された試料でありうる。試料の例は、対象の組織、細胞、血液、及び/又は体液を含むがこれらに限定されない。対象が、がんを有するか、又はがんを有する可能性がある医学的対象である場合、試料は、例えば、好ましくは、生検手順又は他の試料抽出手順を介して、がん病変、又はがんが疑われる病変、又は転移性腫瘍、又は体液が存在する体腔であって、がん細胞で汚染された体腔(例えば、胸腔又は腹腔又は膀胱腔)、又はがん細胞を含有する他の体液などから得られた試料でありうる。試料が抽出される細胞はまた、血液悪性腫瘍(白血病又はリンパ腫など)に由来する腫瘍細胞でもある。場合によって、細胞試料はまた、循環腫瘍細胞、すなわち、血流に入り、適切な単離技法、例えば、アフェレーシス又は従来の静脈の血採取法を使用して抽出される腫瘍細胞でもある。血液以外に、試料が抽出される体液は、尿、消化管内容物、又は滲出物である。本明細書で使用される「試料」という用語はまた、例えば、対象の組織及び/又は細胞及び/又は体液が、対象から採取され、例えば、顕微鏡のスライドガラス上に置かれ、特許請求される方法を実施するために、この試料の一部が、例えば、レーザーキャプチャーマイクロダイセクション(LCM)により、又はスライドガラスから、関心のある細胞を削ぎ落とすことにより、又は蛍光活性化細胞分取技法により抽出される場合も包摂する。加えて、本明細書で使用される「試料」という用語はまた、例えば、対象の組織及び/又は細胞及び/又は体液が、対象から採取され、 顕微鏡のスライドガラス上に置かれ、特許請求される方法が、スライドガラス上で実施される場合も包摂する。加えて、本明細書で使用される「試料」という用語はまた、例えば、細胞系及び/又は細胞培養物が、対象から採取された細胞/組織/体液に基づき作出される場合も包摂する。
開示される別の態様に従い、対象におけるNotch細胞シグナル伝達経路の活性を推測するための装置は、本明細書で記載される、本発明の方法を実施するように構成されたデジタルプロセッサを含む。
開示される別の態様に従い、対象におけるNotch細胞シグナル伝達経路の活性を推測するための非一時的記憶媒体は、本明細書で記載される、本発明の方法を実施するための、デジタル処理デバイスにより実行可能な命令を保存する。非一時的記憶媒体は、ハードドライブ若しくは他の磁気記憶媒体、光学ディスク若しくは他の光学記憶媒体、ランダムアクセスメモリ(RAM)、リードオンリーメモリ(ROM)、フラッシュメモリなどのコンピュータ−読取り型記憶媒体、又は他の電子記憶媒体、ネットワークサーバーなどである。デジタル処理デバイスは、携帯型デバイス(例えば、パーソナルデータアシスタント又はスマートフォン)、ノート型コンピュータ、デスクトップ型コンピュータ、タブレット型コンピュータ又はタブレット型デバイス、リモートネットワークサーバーなどである。
開示される別の態様に従い、対象におけるNotch細胞シグナル伝達経路の活性を推測するためのコンピュータプログラムは、コンピュータプログラムが、デジタル処理デバイス上で実行される場合に、デジタル処理デバイスに、本明細書で記載される、本発明の方法を実施させるためのプログラムコード手段を含む。デジタル処理デバイスは、携帯型デバイス(例えば、パーソナルデータアシスタント又はスマートフォン)、ノート型コンピュータ、デスクトップ型コンピュータ、タブレット型コンピュータ又はタブレット型デバイス、リモートネットワークサーバーなどでありえる。
開示される別の態様に従い、対象の試料中のNotch細胞シグナル伝達経路の3つ以上の、例えば、3つ、4つ、5つ、6つ、7つ、8つ、9つ、10、11、12、以上の標的遺伝子の発現レベルを測定するためのキットは、
対象の試料中の3つ以上のNotch標的遺伝子の発現レベルを決定するための1つ又は複数の構成要素
を含み、
3つ以上のNotch標的遺伝子は、CD28、CD44、DLGAP5、DTX1、EPHB3、FABP7、GFAP、GIMAP5、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、KLF5、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、PTCRA、SOX9、及びTNCからなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子は、DTX1、HES1、HES4、HES5、HEY2、MYC、NRARP、及びPTCRAからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子は、CD28、CD44、DLGAP5、EPHB3、FABP7、GFAP、GIMAP5、HES7、HEY1、HEYL、KLF5、NFKB2、NOX1、PBX1、PIN1、PLXND1、SOX9、及びTNCからなる群から選択される。
3つ以上のNotch標的遺伝子が、CD44、DTX1、EPHB3、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、CD44、EPHB3、HES7、HEY1、HEYL、NFKB2、NOX1、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択されるキット(本明細書で記載される)が特に好ましい。
3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択されるキット(本明細書で記載される)が特に好ましい。
3つ以上のNotch標的遺伝子の発現レベルを測定するための1つ又は複数の構成要素又は手段は、DNAアレイチップ、オリゴヌクレオチドアレイチップ、タンパク質アレイチップ、抗体、複数のプローブ、例えば、標識されたプローブ、RNA逆転写酵素シーケンシング構成要素のセット、及び/又はcDNA、増幅プライマーを含む、RNA若しくはDNAからなる群から選択されうる。ある実施形態では、キットは、本明細書で記載される、3つ以上のNotch標的遺伝子の、mRNA配列又はcDNA配列の部分を対象とする標識されたプローブのセットを含む。ある実施形態では、キットは、3つ以上のNotch標的遺伝子の、mRNA配列又はcDNA配列の部分を対象とするプライマー及びプローブのセットを含む。ある実施形態では、標識されたプローブは、標準化された96ウェルプレート内に含有される。ある実施形態では、キットは、基準遺伝子のセットを対象とするプライマー又はプローブをさらに含む。このような基準遺伝子は、例えば、本明細書で記載される、標的遺伝子発現レベルの正規化又は標準化に有用な、構成的に発現させた遺伝子でありうる。
ある実施形態では、対象の試料中のNotch細胞シグナル伝達経路の3つ以上の、例えば、3つ、4つ、5つ、6つ、7つ、8つ、9つ、10、11、12、以上の標的遺伝子の発現レベルを測定するためのキットは、
3つ以上のNotch標的遺伝子を対象とするポリメラーゼ連鎖反応プライマーと、
3つ以上のNotch標的遺伝子を対象とするプローブと
を含み、
3つ以上のNotch標的遺伝子は、CD28、CD44、DLGAP5、DTX1、EPHB3、FABP7、GFAP、GIMAP5、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、KLF5、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、PTCRA、SOX9、及びTNCからなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子は、DTX1、HES1、HES4、HES5、HEY2、MYC、NRARP、及びPTCRAからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子は、CD28、CD44、DLGAP5、EPHB3、FABP7、GFAP、GIMAP5、HES7、HEY1、HEYL、KLF5、NFKB2、NOX1、PBX1、PIN1、PLXND1、SOX9、及びTNCからなる群から選択される。
3つ以上のNotch標的遺伝子が、CD44、DTX1、EPHB3、HES1、HES4、HES5、HES7、HEY1、HEY2、HEYL、MYC、NFKB2、NOX1、NRARP、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、CD44、EPHB3、HES7、HEY1、HEYL、NFKB2、NOX1、PBX1、PIN1、PLXND1、及びSOX9からなる群から選択されるキット(本明細書で記載される)が特に好ましい。
3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、好ましくは、ここで、2つ以上の、例えば、3つ、4つ、5つ、6つ、以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数の、例えば、3つ、4つ、以上のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択されるキット(本明細書で記載される)が特に好ましい。
開示される別の態様に従い、対象におけるNotch細胞シグナル伝達経路の活性を推測するためのキットは、
本明細書で記載される、本発明のキットと、
本明細書で記載される、本発明の装置、本明細書で記載される、本発明の非一時的記憶媒体、又は本明細書で記載される、本発明のコンピュータプログラムと
を含む。
開示される別の態様に従い、本明細書で記載される、本発明のキットは、本明細書で記載される、本発明の方法の実施において使用される。
本明細書で記載される本発明は、例えば、以下のアクティビティ:
対象におけるNotch細胞シグナル伝達経路の推測された活性に基づく診断;
対象におけるNotch細胞シグナル伝達経路の推測された活性に基づく予後診断;
対象におけるNotch細胞シグナル伝達経路の推測された活性に基づく薬物処方;
対象におけるNotch細胞シグナル伝達経路の、論された活性に基づく薬物有効性の予測;
対象におけるNotch細胞シグナル伝達経路の推測された活性に基づく有害作用の予測;
薬物有効性のモニタリング;
薬物の開発;
アッセイの開発;
経路の探索;
がんの病期分類;
対象における、Notch細胞シグナル伝達経路の推測された活性に基づく臨床試験における対象の登録;
実施される後続の試験の選択;及び
コンパニオン診断試験の選択
のうちの少なくとも1つにおいてもまた、有利に使用されうる。
さらなる利点は、添付の図面、以下の記載を読み、これらを理解すれば、特に、本明細書の下記に提示される詳細な例を読めば、当業者に明らかであろう。
請求項1に記載の方法、請求項9に記載の装置、請求項10に記載の非一時的記憶媒体、請求項11に記載のコンピュータプログラム、請求項12から14のいずれか一項に記載のキット、及び請求項15に記載のキットの使用は、特に、従属請求項において規定される通り、類似の好ましい実施形態及び/又は同一の好ましい実施形態を有することを理解されたい。
本発明の好ましい実施形態はまた、従属請求項又は上記の実施形態の、それぞれの独立請求項との任意の組合せでもありうることを理解されたい。
本発明のこれらの態様及び他の態様は、本明細書の後出で記載される実施形態から明らかであり、これらを参照することにより解明されるであろう。
Notch細胞シグナル伝達経路を、概略的、かつ、例示的に示す図である。経路は、Notch細胞外ドメインが、DSLリガンドに結合すると、活性化する。受容体の切断の後、Notch細胞内ドメインは、核へと移行し、他のタンパク質と併せて、活性の転写因子複合体を形成する(Guruharsha K.G.ら、「The Notch signaling system:recent insights into the complexity of a conserved pathway」、Nature Reviews Genetics、13巻、2012年9月、654〜666頁を参照されたい;「TS」=転写スイッチ;「TG」=標的遺伝子)。 本明細書では、ベイジアンネットワークモデルである、数学的モデルであって、Notch細胞シグナル伝達経路の転写プログラムをモデル化するのに使用される数学的モデルを、概略的、かつ、例示的に示す図である。 対象の試料中で測定される、Notch細胞シグナル伝達経路の標的遺伝子の発現レベルに基づき、対象におけるNotch細胞シグナル伝達経路の活性を推測するための工程を例示的に図解するためのフローチャートを示す図である。 本明細書で記載される、較正された数学的経路モデルを得るための工程を例示的に図解するためのフローチャートを示す図である。 本明細書で記載される、対象の試料中の、Notch転写因子(TF)エレメントの活性レベルを決定するための工程を例示的に図解するためのフローチャートを示す図である。 離散可観測量を使用して、対象におけるNotch細胞シグナル伝達経路の活性を推測するための工程を例示的に図解するためのフローチャートを示す図である。 連続可観測量を使用して、対象におけるNotch細胞シグナル伝達経路の活性を推測するための工程を例示的に図解するためのフローチャートを示す図である。 Notch細胞シグナル伝達経路の標的遺伝子についてのRT−qPCR解析から、Cq値を決定するための工程を例示的に図解するためのフローチャートを示す図である。 11例の正常卵巣(群1)試料、及び20例の高悪性度漿液性乳頭状卵巣癌(群2)試料について公開されている発現データセット(データセット:GSE2109、GSE9891、GSE7307、GSE18520、GSE29450、GSE36668から採取された試料のサブセット)を使用する、表2による、18の標的遺伝子についてのショートリスト、及び本明細書で記載される方法に基づく、ベイジアンネットワークモデルの較正結果を示す図である。 11例の正常卵巣(群1)試料、及び20例の高悪性度漿液性乳頭状卵巣癌(群2)試料について公開されている発現データセット(データセット:GSE2109、GSE9891、GSE7307、GSE18520、GSE29450、GSE36668から採取された試料のサブセット)を使用する、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)、及び本明細書で記載される方法に基づく、ベイジアンネットワークモデルの較正結果を示す図である。 データセット:GSE6495に由来する、MOLT4細胞系の、3つの独立の培養物についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 データセット:GSE6495に由来する、MOLT4細胞系の、3つの独立の培養物についての、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 誘導型Notch3細胞内構築物をトランスフェクトされたIMR32細胞についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 漸増用量のプラスチック固定化Notchリガンドデルタ1 ext−IgGを伴い、72時間にわたり培養されたCD34+CD45RA−Lin−HPC(データセット:GSE29524)についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 高Notch活性を有することが知られているCUTLL1細胞についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 高Notch活性を有することが知られているCUTLL1細胞についての、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 COUP−TFII siRNAをトランスフェクトされたHUVEC細胞(データセット:GSE33301)についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 GSE6532、GSE9195、GSE12276、GSE20685、GSE21653、及びEMTAB365に由来する試料中の乳がん亜群についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 プラスチック固定化Notchリガンドデルタ1 ext−IgGの段階的投与を伴い、72時間にわたり培養されたCD34+CD45RA−Lin−HPC(データセット:GSE29524)についての、表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 高Notch活性を有することが知られているCUTLL1細胞についての、表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)、及び表3による、12の標的遺伝子についてのショートリストのそれぞれを使用してトレーニングされた、例示的なベイジアンネットワークモデルの間の相関を示す図である。 10のNotch標的遺伝子についてのリストと対比した7つのNotch標的遺伝子についてのリストを使用する、Notch細胞シグナル伝達経路活性の予測についての比較を示す図である。 12のNotch標的遺伝子についてのリストと対比した8つのNotch標的遺伝子についてのリストを使用する、Notch細胞シグナル伝達経路活性の予測についての比較を示す図である。 2例の対照胚性幹細胞(ESc)、2例の対照中胚葉前駆細胞(MPc)、OHT処理されていない、タモキシフェン誘導型NERT構築物(NotchIC)を含有する、2例のESc試料、OHT処理された、タモキシフェン誘導型NERT構築物(NotchIC)を含有する、2例のESc試料、OHT処理されていない、タモキシフェン誘導型NERT構築物(NotchIC)を含有する、4例のMPc試料、及びOHT処理された、タモキシフェン誘導型NERT構築物(NotchIC)を含有する、4例のMPc試料を含有する、公開されている発現データセットであるGSE15268を使用する、表6による、10の標的遺伝子についてのマウスリスト、及び本明細書で記載される方法に基づく、ベイジアンモデルの較正結果を示す図である。 構成的に活性の誘導型Notch1細胞内ドメイン(NICD1)を伴うマウス乳腺(データセット:GSE51628)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 Notch1を活性化させる、条件付きトランスジェニック系を伴うマウス卵黄嚢組織、及びRBPJ(Notch転写因子複合体の一部)機能喪失を伴うトランスジェニックマウスに由来するマウス卵黄嚢組織(データセット:GSE22418)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。 Notch2受容体の、条件付き機能獲得対立遺伝子を伴うマウス骨髄細胞(成体骨髄赤血球前駆細胞)(データセット:GSE46724)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す図である。
以下の例は、特に、好ましい方法、及びこれと関連する、選択された態様を例示するだけのものである。例においてもたらされる教示は、例えば、Notch細胞シグナル伝達経路の異常な活性を検出、予測、及び/又は診断するためのいくつかの試験及び/又はキットの構築に使用される。さらに、本明細書で記載される方法を使用すると、薬物処方を、有利に誘導することができ、薬物応答の予測及び薬物の有効性(及び/又は有害作用)のモニタリングを行うことができ、薬物耐性を予測及びモニタリングして、例えば、(コンパニオン診断試験のような)実施される後続の試験を選択することができる。以下の例は、本発明の範囲を限定するものとは見なさないものとする。
[実施例1]
数学的モデルの構築
国際特許出願公開第WO2013/011479A2号(「標的遺伝子発現の確率モデル化を使用する、細胞シグナル伝達経路活性の評価」)において、詳細に記載されている通り、確率モデル、例えば、ベイジアンネットワークモデルを構築し、本明細書では、Notch細胞シグナル伝達経路である細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルと、本明細書では、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御するTFエレメントであるNotch TFエレメントである転写因子(TF)エレメントの活性レベルとの間に、条件付き確率関係を組み込むことにより、このようなモデルを使用して、細胞シグナル伝達経路の活性を、高精度で決定することができる。さらに、確率モデルは、条件付き確率を補正し、かつ/又はさらなる情報源を表すように、新たなノードを、モデルへと付加することにより、その後の臨床研究により得られる、さらなる知見を組み込むように、たやすく更新することができる。このようにして、確率モデルは、最新の医学的知見を具体化するために、適宜更新されうる。
国際特許出願公開第WO2014/102668A2号(「標的遺伝子発現の線形結合を使用する、細胞シグナル伝達経路活性の評価」)において、詳細に記載されている、理解及び解釈が容易な、別の手法では、本明細書では、Notch細胞シグナル伝達経路である、細胞シグナル伝達経路の活性は、細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルと、本明細書では、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御するTFエレメントである、Notch TFエレメントである転写因子(TF)エレメントのレベルとの間の関係を組み込むことで、線形モデル又は(疑似)線形モデルを構築及び査定することにより、決定することができ、この線形モデルは、3つ以上の標的遺伝子の発現レベルの1つ又は複数の線形結合に基づく。
いずれの手法でも、3つ以上の標的遺伝子の発現レベルは、好ましくは、例えば、標的遺伝子のmRNA配列と関連するプローブを使用する、(RT)PCR技法及びマイクロアレイ技法、並びにRNAシーケンシングの結果でありうる、mRNAレベルの測定値である。別の実施形態では、3つ以上の標的遺伝子の発現レベルは、タンパク質レベル、例えば、標的遺伝子によりコードされるタンパク質の濃度及び/又は活性により測定することもできる。
前述の発現レベルは、任意選択で、適用によく適合する場合もあり、よく適合しない場合もある、多くの方式で転換される。例えば、発現レベル、例えば、マイクロアレイベースのmRNAレベルの、4つの異なる変換は、
「連続データ」、すなわち、MAS5.0及びfRMAなど、よく知られているアルゴリズムを使用して、マイクロアレイの前処理の後で得られる発現レベル、
「zスコア」、すなわち、全ての試料にわたる平均が0であり、標準偏差が1であるようにスケーリングされた、連続発現レベル、
「離散」、すなわち、ある特定の閾値を上回る、あらゆる発現を1とし、これを下回る値を0とする(例えば、プローブセットについての閾値は、多数の陽性臨床試料、及び同じ数の陰性臨床試料のセット内のその値の(重み付け)中央値として選び出される)、
「ファジー」、すなわち、以下のフォーマット:1/(1+exp((thr−expr)/se))のシグモイド関数を使用して、連続発現レベルが、0〜1の間の値へと転換されるが、式中、exprは、連続発現レベルであり、thrは、前述の閾値であり、seは、0〜1の間の差違に影響を及ぼすソフトニングパラメータである、連続発現レベル
である。
構築されうる、最も簡単な線形モデルのうちの1つは、本明細書では、Notch TFエレメントである、転写因子(TF)エレメントを表すノードを、第1の層に有し、例えば、マイクロアレイ実験又は(q)PCR実験において、特定の標的遺伝子と、特に、高度に相関する、例えば、1つのプローブセットによる、標的遺伝子発現レベルの直接的測定を表す、重み付けノードを、第2の層に有するモデルである。重みは、トレーニングデータセットに由来する計算に基づくか、計算専門家の知見に基づくかのどちらかである。標的遺伝子1つ当たり、おそらく、複数の発現レベルが測定される場合に(例えば、1つの標的遺伝子が、複数のプローブセットにより測定されうる、マイクロアレイ実験の場合に)、標的遺伝子1つ当たり、1つの発現レベルだけを使用する、この手法は、特に、簡単である。特定の標的遺伝子に使用される、1つの発現レベルを選択する、具体的な方式は、トレーニングデータセットの活性試料と、不活性試料とを、最も良好に分離することが可能なプローブセットに由来する発現レベルを使用することである。このプローブセットを決定する1つの方法は、統計学的検定、例えば、t検定を実施し、p値が最小となるプローブセットを選択することである。p値が最小となるプローブセットの、トレーニングデータセットの発現レベルとは、定義上、(知られている)活性試料の発現レベルと、及び不活性試料の発現レベルとが重複する確率が最小となるプローブセットの発現レベルである。別の選択方法は、オッズ比に基づく。このようなモデルでは、1つ又は複数の発現レベルは、3つ以上の標的遺伝子の各々についてもたらされ、1つ又は複数の線形結合は、3つ以上の標的遺伝子の各々について、重み付け項を含む線形結合を含み、各重み付け項は、それぞれの標的遺伝子についてもたらされる、1つ又は複数の発現レベルのうちの、1つの発現レベルだけに基づく。上記で記載した標的遺伝子1つ当たり、1つの発現レベルだけが選び出される場合、モデルを、「最も弁別力のあるプローブセット」モデルと呼ぶ。
「最も弁別力のあるプローブセット」モデルに対する代替法では、標的遺伝子1つ当たり、おそらく、複数の発現レベルが測定される場合に、標的遺伝子1つ当たりにもたらされる、全ての発現レベルを使用することも可能である。このようなモデルでは、1つ又は複数の発現レベルは、3つ以上の標的遺伝子の各々についてもたらされ、1つ又は複数の線形結合は、3つ以上の標的遺伝子についてもたらされる、1つ又は複数の発現レベルの、全ての発現レベルの線形結合を含む。言い換えれば、3つ以上の標的遺伝子の各々について、それぞれの標的遺伝子についてもたらされる、1つ又は複数の発現レベルの各々は、線形結合において、それ自身の(個別の)重みにより重み付けされる。この変化形は、「全てのプローブセット」モデルと呼ばれる。「全てのプローブセット」は、もたらされる全ての発現レベルを使用しながら、比較的簡単であることの利点を有する。
上記で記載したいずれのモデルも、それらが、TFエレメントの活性レベルが、3つ以上の標的遺伝子についての、1つ又は複数のプローブセットの発現レベルの線形結合に基づき計算される「単一層」モデルと考えられるものであるという点を共有する。
それぞれのモデルを査定することにより、本明細書では、Notch TFエレメントである、TFエレメントの活性レベルを決定した後で、本明細書では、Notch細胞シグナル伝達経路である、細胞シグナル伝達経路の活性を推測するために、決定されたTFエレメントの活性レベルを閾値化することができる。このような適切な閾値を計算するのに好ましい方法は、不活性の細胞シグナル伝達経路を有することが知られているトレーニング試料の、決定されたTFエレメントの活性レベルwlc(重み付け線形結合)と、活性の細胞シグナル伝達経路を伴うトレーニング試料のwlcとを比較することによる方法である。このようにする方法であって、また、これらの群における分散も考慮する方法は、閾値
Figure 2020536530
を使用することにより与えられ、式中、σ及びμは、トレーニング試料についての、決定された、TFエレメントの活性レベルwlcの標準偏差及び平均値である。活性のトレーニング試料中及び/又は不活性のトレーニング試料中で、少数の試料だけが利用可能な場合、2つの群の分散の平均に基づき、シュードカウントを計算された分散に加算する:
Figure 2020536530
式中、νは、群について決定されたTFエレメントの活性レベルwlcの分散であり、xは、正のシュードカウント、例えば、1又は10であり、nact及びnpasは、活性試料及び不活性試料のそれぞれの数である。次に、分散νの平方根を取ることにより、標準偏差σを得ることができる。
解釈を容易とするために、閾値を、決定されたTFエレメントの活性レベルであるwlcから差し引く結果として、負の値が、不活性の細胞シグナル伝達経路に対応し、正の値が、活性の細胞シグナル伝達経路に対応する、細胞シグナル伝達経路の活性スコアをもたらすことができる。
上記で記載された「単一層」モデルに対する代替法として、例ではまた、「二層」も使用される。このようなモデルでは、合計値は、その関連するプローブセットの、測定された強度に基づく線形結合を使用して、あらゆる標的遺伝子について計算される(「第1(下)層」)。その後、さらなる線形結合を使用して、計算された合計値を、細胞シグナル伝達経路の、他の標的遺伝子についての合計値と組み合わせる(「第2(上)層」)。ここでもまた、重みは、トレーニングデータセットから学習されるか、専門家の知見に基づくかのどちらか、又はこれらの組合せもある。言い換えれば、「二層」モデルでは、1つ又は複数の発現レベルは、3つ以上の標的遺伝子の各々についてもたらされ、1つ又は複数の線形結合は、3つ以上の標的遺伝子の各々について、それぞれの標的遺伝子についてもたらされる、1つ又は複数の発現レベルの、全ての発現レベルのうちの、第1の線形結合を含む。モデルは、3つ以上の標的遺伝子の各々について、重み付け項を含む、さらなる線形結合にさらに基づき、各重み付け項は、それぞれの標的遺伝子についての、第1の線形結合に基づく(「第2(上)層」)。
「二層」モデルの好ましい変化形では、合計値の計算は、トレーニングデータを使用して、各標的遺伝子について、閾値を規定することと、計算された線形結合から、閾値を差し引くことと、標的遺伝子の合計を求めることとを含む。この場合、閾値は、負の標的遺伝子合計値が、下方調節される標的遺伝子に対応し、正の標的遺伝子合計値が、上方調節される標的遺伝子に対応するように選び出される。また、例えば、上記で記載された変換(ファジー変換、離散変換など)のうちの1つを使用して、標的遺伝子の合計値を、「第2(上)層」において組み合わせる前に変換することも可能である。
「二層」モデルを査定することにより、TFエレメントの活性レベルを決定した後で、上記で記載した通りに、細胞シグナル伝達経路の活性を推測するために、決定されたTFエレメントの活性レベルを、閾値化することができる。
以下では、上記で記載したモデルを、「(疑似)線形」モデルとして、総体的に表示する。トレーニング及び確率モデル、例えば、ベイジアンネットワークモデルの使用についてのより詳細な記載は、下記の実施例3に提示される。
[実施例2]
標的遺伝子の選択
転写因子(TF)とは、特異的DNA配列への結合により、標的遺伝子からの転写を調節し、これにより、DNAから、mRNAへの遺伝情報の転写を制御することことが可能な、タンパク質複合体(すなわち、特異的構造内で一体に結合したタンパク質の組合せ)又はタンパク質である。本明細書では、TF複合体のこの作用により直接的に産生されるmRNAを、「直接的標的遺伝子」(転写因子の)と称する。細胞内シグナル伝達経路の活性化はまた、「間接的標的遺伝子」と称する、より副次的な遺伝子転写も結果としてもたらす。以下では、細胞シグナル伝達経路活性と、mRNAレベルとの直接的連関としての、直接的標的遺伝子を含むか、又はこれからなる、(疑似)線形モデル、又はベイジアンネットワークモデル(例示的な数学的モデルとしての)が好ましいが、直接的標的遺伝子と、間接的標的遺伝子との区別は、常に明白であるわけではない。本明細書では、利用可能な学術文献データに基づくスコアリング関数を使用して、直接的標的遺伝子を選択する方法が提示される。しかしながら、限定的な情報のほか、生物学的ばらつき、及び不確実性に起因して、間接的標的遺伝子の、意想外の選択を除外することはできない。標的遺伝子を選択するために、「www.ncbi.nlm.nih.gov/pubmed」においてアクセス可能であり、本明細書ではさらに、「Pubmed」とも称する、米国国立保健研究所のMEDLINEデータベースを援用して、標的遺伝子のリストを作成した。さらに、それらの発現の証拠となる性質に基づき、標的遺伝子の、2つのさらなるリストを選択した。
2016年の第4四半期〜2017年の第1四半期の期間において、「“Notch”AND“target gene”」などのクエリーを使用して、推定Notch標的遺伝子を含有する刊行物を検索した。Notch経路は、胚系統に応じて異なる(が、重複する)標的遺伝子プロファイルを活性化させる、胚性経路である(Meier−Stiegen F.ら、「Activated Notch1 target genes during embryonic cell differentiation depend on the cellular context and include lineage determinants and inhibitors」、PLoS One、5巻、7号、2010年7月を参照されたい)。検索は、外胚葉及び内胚葉由来の臓器/組織/細胞内で発現する標的遺伝子に、特段の重点を置いて、3つの異なる胚系統(外胚葉、内胚葉、中胚葉)からの、細胞型/組織/臓器派生物の間で、示差的に発現する、標的遺伝子のセットに焦点を絞った。結果として得られる刊行物を、下記で、より詳細に記載される方法に従い、手作業で、さらに解析した。
証拠が蓄積される科学的実験の種類に応じて、特異的な標的遺伝子についての科学的証拠に評定を与える、ランク付けシステムを使用することにより、特異的な細胞シグナル伝達経路のmRNA標的遺伝子を、学術文献から選択した。一部の実験の証拠が、例えば、Notch細胞シグナル伝達経路が活性であることが知られている細胞系のマイクロアレイ上の、プローブセットの強度の増大により検出される、mRNAの増大のように、遺伝子が、直接的標的遺伝子であることを示唆するに過ぎないのに対し、他の証拠は、同定されたNotch細胞シグナル伝達経路のTF結合部位と、細胞内の特異的細胞シグナル伝達経路の刺激の後の、クロマチン免疫沈降(ChIP)アッセイにおける、この部位の回収と、細胞系内の細胞シグナル伝達経路の、特異的な刺激の後における、mRNAの増大との組合せのように、極めて強力でありうる。
学術文献では、特異的な細胞シグナル伝達経路の標的遺伝子を見出すための、いくつかの種類の実験を同定することができる:
1.関心のある細胞シグナル伝達経路のTFの、ゲノム上のその結合部位への、直接的結合を示すChIP実験。例:例えば、Notchリガンドを伴う刺激、又はNICDのトランスフェクションによる、Notch細胞シグナル伝達経路の能動的誘導を伴うクロマチン免疫沈降(ChIP)技術、及びこの能動的誘導を伴わないChIP技術を使用することにより、細胞系のDNA内の、推定機能性Notch TF結合部位を、ヌクレオチド配列に純粋に基づいて認識される結合部位のサブセットとして同定した。推定機能性は、TFが、DNA結合部位に結合することを見出す、ChIP由来の証拠として同定した。
2.in vitroにおける、TFの、結合配列を含有するDNAの断片への結合を示す、電気泳動移動度シフトアッセイ(EMSA)。ChIPベースの証拠と比較して、EMSAベースの証拠は、in vivo状況へと置きかえることができないので、それほど強力ではない。
3.細胞シグナル伝達経路の刺激、及びマイクロアレイ、RNAシーケンシング、定量的PCR、又は他の技法を使用する、mRNAの発現の測定であり、Notch細胞シグナル伝達経路誘導型細胞系を使用し、少なくとも1つであるが、好ましくは、誘導(タンパク質への翻訳を阻害するので、誘導されたmRNAを、直接的標的遺伝子であると仮定する、シクロヘキシミドの存在下における)の後の、いくつかの時点において測定される、mRNAプロファイルを測定する。
4.3と同様であるが、代替的に、さらに下流における、ウェスタンブロットなど、タンパク質の存在度の測定を伴いながら、mRNAの発現を測定する。
5.Notch阻害剤、例えば、ガンマ−セクレターゼ阻害剤(GSI)を使用する、細胞シグナル伝達経路の阻害、及びマイクロアレイ、RNAシーケンシング、定量的PCR、又は他の技法を使用する、mRNAの発現の測定であり、Notch細胞シグナル伝達経路活性細胞系を使用し、少なくとも1つであるが、好ましくは、阻害の後の、いくつかの時点において測定される、mRNAプロファイルを測定する。
6.5と同様であるが、代替的に、さらに下流における、ウェスタンブロットなど、タンパク質の存在度の測定を伴いながら、mRNAの発現を測定する。
7.バイオインフォマティクス手法を使用する、ゲノム内のTF結合部位の同定。Notch TFエレメントについての例:CSL/RBP−J結合モチーフである5’−CGTGGGAA−3’を使用して、ソフトウェアプログラムを、ヒトゲノム配列上で実行し、潜在的結合部位を、遺伝子プロモーター領域内、及び他のゲノム領域内の両方において同定した。
8.シクロヘキシミドの非存在下であることを除き、3と同様である。
9.シクロヘキシミドの非存在下であることを除き、4と同様である。
最も簡単な形態では、遺伝子を、転写因子のNotchファミリーの標的遺伝子であると同定した、これらの実験の手法の各々について、あらゆる潜在的な遺伝子に、1ポイントを与えることができる。この相対ランク付け戦略を使用して、最も信頼できる標的遺伝子のリストを作成することができる。
代替的に、in vivoにおける直接的標的遺伝子の、最も大きな証拠をもたらす技術に、大きな数のポイントを与えることによる、別の方式によるランク付けを使用して、直接的標的遺伝子である可能性が最も高い標的遺伝子を同定することができる。上記のリストでは、これは、実験の手法1)については、9ポイント、2)については、8ポイントを意味し、実験の手法9)については、1ポイントへと低下することを意味する。このようなリストは、「標的遺伝子についての一般リスト」と呼ばれる。
生物学的変動及び不確実性にもかかわらず、本発明者らは、直接的標的遺伝子が、組織非依存的に誘導される可能性が最も高いことを仮定した。これらの標的遺伝子のリストは、「標的遺伝子についての、証拠をキュレートしたリスト」と呼ばれる。このような、標的遺伝子についての、証拠をキュレートしたリストを使用して、異なる組織供給源に由来する試料へと適用されうる、Notch細胞シグナル伝達経路の算出モデルを構築した。
以下は、Notch細胞シグナル伝達経路について、証拠をキュレートした標的遺伝子リストの選択を、どのように具体的に構築したのかを、例示的に説明する。
ChIP、EMSA、示差的発現、ノックダウン/ノックアウト、ルシフェラーゼ遺伝子レポーターアッセイ、刊行物において報告された配列解析など、各種の実験の証拠にポイントを与える、スコアリング関数を導入した。さらなる解析を実施して、1又は2種類の実験の証拠、例えば、示差的発現だけではなく、多様な種類の実験の証拠を有する遺伝子だけを許容した。2種類を超える、利用可能な実験の証拠を有する遺伝子(表1において示される)を選択した。
本発明者らは、標的遺伝子についての、証拠をキュレートしたリスト(表2に列挙される、「18の標的遺伝子についてのショートリスト」)のさらなる選択を行った。この選択は、証拠をキュレートしたリストの標的遺伝子であって、比較的少数の証拠、例えば、見出された1つの論考だけにおいて見出され、かつ/又は例えば、血液又は脳組織に高度に特異的な証拠を有する標的遺伝子を除外することにより行った。トレーニング試料から、Notchシグナル伝達経路の活性を決定するときに、強力な証拠となる「18の標的遺伝子についてのショートリスト」の標的遺伝子を、「12の標的遺伝子についてのショートリスト」(表3に列挙される、「12の標的遺伝子についてのショートリスト」)のために選択した。本明細書では、遺伝子発現オムニバス(GEO、www.ncbi.nlm.nih.gov/geo/;最終アクセス日:2016年12月3日)による高悪性度漿液性乳頭状卵巣がん患者(GSE2109及びGSE9891から取り出された、Notch活性サブセット)のセットと、正常卵巣組織試料の対応するセット(GSE7307、GSE18520、GSE29450及びGSE36668から取り出された、Notch不活性サブセット)のそれぞれに由来する患者試料との間で、最大のオッズ比(下記を参照されたい)を有し、かつ/又は証拠のランク付けで、極めて高いスコアを示した、12の標的遺伝子を選択した。
表1: Notch細胞シグナル伝達経路モデルにおいて使用される、Notch細胞シグナル伝達経路の「標的遺伝子についての、証拠をキュレートしたリスト」(26の標的遺伝子リスト)、及び標的遺伝子のmRNAの発現レベルを測定するのに使用される、関連するプローブセット
Figure 2020536530
表2: Notch標的遺伝子についての、証拠をキュレートしたリストに基づく、Notch標的遺伝子の「18の標的遺伝子についてのショートリスト」(関連するプローブセットは、表1の場合と同じである
Figure 2020536530
表3: Notch標的遺伝子についての、証拠をキュレートしたリストに基づく、Notch標的遺伝子の「12の標的遺伝子についてのショートリスト」(関連するプローブセットは、表1の場合と同じである)
Figure 2020536530
[実施例3]
数学的モデルのトレーニング及び使用
数学的モデルを使用して、対象における、本明細書では、Notch細胞シグナル伝達経路である、細胞シグナル伝達経路の活性を推測しうる前に、モデルを、適切にトレーニングしなければならない。
数学的経路モデルが、Notch TFエレメントの活性レベルと、対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルとを関連付ける条件付き確率に基づく確率モデル、例えば、ベイジアンネットワークモデルである場合、トレーニングは、好ましくは、国際特許出願公開第WO2013/011479A2号(「標的遺伝子発現の確率モデル化を使用する、細胞シグナル伝達経路活性の評価」)において、詳細に記載されている通りに実施される。
数学的経路モデルが、対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルの1つ又は複数の線形結合に基づく場合、トレーニングは、好ましくは、国際特許出願公開第WO2014/102668A2号(「標的遺伝子発現の線形結合を使用する、細胞シグナル伝達経路活性の評価」)において、詳細に記載されている通りに実施される。
本明細書では、図2に示される、例示的なベイジアンネットワークモデルを使用して、Notch細胞シグナル伝達経路の転写プログラムを、簡単にモデル化した。モデルは、3種類のノード:(a)第1の層である1における、転写因子(TF)エレメント(「存在しない」状態及び「存在する」状態を伴う);(b)第2の層である2における、標的遺伝子TG、TG、TG(「低下」状態及び「上昇」状態を伴う)、及び;(c)第3の層である3における、標的遺伝子の発現レベルへと連関させた測定ノードからなる。これらは、好ましくは、本明細書で使用されるマイクロアレイプローブセットである、PS1,1、PS1,2、PS1,3、PS2,1、PSn,1、PSn,m(「低度」状態及び「高度」状態を伴う)でありうるが、また、RNAseq又はRT−qPCRなど、他の遺伝子発現測定値でもありうる。
本明細書では、例示的なベイジアンネットワークモデルである数学的モデルの適切な実施は、マイクロアレイデータに基づく。モデルは、(i)標的遺伝子の発現レベルが、TFエレメントの活性化に、どのくらい依存するのか、及び(ii)同様に、プローブセットの強度が、それぞれの標的遺伝子の発現レベルに、どのくらい依存するのかについて記載する。(ii)のために、プローブセットの強度は、遺伝子発現オムニバス(GEO、www.ncbi.nlm.nih.gov/geo)及びArrayExpress(www.ebi.ac.uk/arrayexpress)から、広く入手可能な、fRMAにより前処理された、Affymetrix HG−U133Plus2.0マイクロアレイから取る。
例示的なベイジアンネットワークモデルは、本明細書では、Notch細胞シグナル伝達経路である、細胞シグナル伝達経路についての生物学の簡略化であり、生物学的測定は、典型的に、ノイズに富むので、確率的手法を選び取る、すなわち、(i)TFエレメントと、標的遺伝子との関係、及び(ii)標的遺伝子と、それらのそれぞれのプローブセットとの関係を、確率的に記載する。さらに、腫瘍の増殖を駆動する、発がん性の細胞シグナル伝達経路の活性は、一過性、かつ、劇的に変更されるのではなく、長期的又はさらに不可逆的に変更されることを仮定した。したがって、静的な細胞状態を解釈するための、例示的なベイジアンネットワークモデルを開発した。この理由で、複雑な動的細胞シグナル伝達経路の特徴は、モデルへと組み込まなかった。
例示的なベイジアンネットワークモデルを、構築し、較正すると(下記を参照されたい)、第3層である3における観測値としてのプローブセットの測定値を入力し、TFエレメントが「存在する」ためには、確率はいくつでなければならなかったのかを、モデルにおいて後ろ向きに推測することにより、モデルを、新たな試料のマイクロアレイデータに対して使用することができる。この場合、「存在する」とは、TFエレメントが、DNAに結合し、細胞シグナル伝達経路の標的遺伝子の転写を制御している現象であると考えられ、「存在しない」とは、TFエレメントが、転写を制御していない現象であると考えられる。よって、この確率は、本明細書では、Notch細胞シグナル伝達経路である、細胞シグナル伝達経路の活性を指し示すのに使用される、一次リードアウトである。これは、次に、細胞シグナル伝達経路が活性である確率に対する細胞シグナル伝達経路が不活性である確率との比(すなわち、オッズは、p/(1−p)により与えられ、ここで、pは、細胞シグナル伝達経路が活性であることの予測確率である)を取ることにより、細胞シグナル伝達経路が活性であるオッズへと置きかえることができる。
例示的なベイジアンネットワークモデルでは、定量的な確率的推測を可能とするために、確率的関係は定量化された。組織型にわたる一般化挙動を改善するために、(i)TFエレメントと、標的遺伝子との確率的関係について記載するパラメータを、注意深く厳選した。TFエレメントが「存在しない」場合、標的遺伝子は「低下」する可能性が極めて高く、よって、0.95の確率を、このために選び出し、0.05の確率を、標的遺伝子が「上昇」することのために選び出す。0.05(非ゼロ)の確率は、標的遺伝子が、他の因子により調節されるか、又は意想外に「上昇」することが観察される(例えば、測定ノイズのために)、(稀少な)確率に相応する。TFエレメントが「存在する」場合、0.70の確率で、標的遺伝子は「上昇」すると考えられ、0.30の確率で、標的遺伝子は「低下」すると考えられる。0.3の値は、標的遺伝子は、TFエレメントが存在する場合であっても、高度に発現しない原因がいくつかあるために、例えば、遺伝子のプロモーター領域がメチル化されているため、このように選び出される。標的遺伝子が、TFエレメントにより上方調節されず、下方調節される場合も、確率は、同様に選び出されるが、TFエレメントの存在時における下方調節を反映する。(ii)標的遺伝子と、それらのそれぞれのプローブセットとの関係について記載するパラメータを、実験データ上で較正した。後者について、この例では、実験データのために、活性のNotch細胞シグナル伝達経路を有することが知られている患者試料に由来するマイクロアレイデータを使用する一方で、異なるデータセットに由来する正常な健常試料を、不活性のNotch細胞シグナル伝達経路試料として使用したが、これはまた、細胞系による実験、又は既知の細胞シグナル伝達経路活性状態を伴う、他の患者試料を使用しても実施されうる。結果として得られる条件付き確率表を、以下により示す:
Figure 2020536530
Figure 2020536530
これらの表では、変数である、ALi,j、AHi,j、PLi,j、and PHi,jは、それぞれ、「低度」(L)又は「高度」(H)のプローブセットの強度を有する、転写複合体が「存在しない」(A)か、又は「存在する」(P)較正試料の数を指し示す。0〜1の極端な確率を回避するために、ダミーカウントを加算した。
観察されたプローブセットの強度を離散化するために、各プローブセットPSi,jについて、それを下回る観測値が「低度」と呼ばれ、それを上回る観測値が「高度」と呼ばれる、閾値ti,jを使用した。この閾値は、使用される較正データセット内のプローブセットの、(重み付け)中央値強度となるように選び出された。マイクロアレイデータのノイズ性のために、観察されたプローブセットの強度を、その閾値と比較する場合に、報告された強度の近傍において、標準偏差を0.25とする(log2スケール上で)正規分布を仮定し、閾値を下回る確率質量及び閾値を上回る確率質量を決定することにより、ファジー法を使用した。
上記で記載した、例示的なベイジアンネットワークの代わりに、上記の実施例1で記載された(疑似)線形モデルを援用した場合、ノードと閾値との相関の符号及び大きさを指し示す重みであって、ノードが「存在しない」のか、「存在する」のかを判定する重みを、モデルを使用して、被験試料中の、細胞シグナル伝達経路活性を推測しうる前に決定する必要がある。専門家の知見を使用して、重み及び閾値を、事前に充填しうるが、典型的には、好ましくは、そのグラウンドトゥルースが知られている、トレーニング試料の代表的セット、例えば、既知の「存在する」転写因子複合体(=活性の細胞シグナル伝達経路)、又は「存在しない」転写因子複合体(=不活性の細胞シグナル伝達経路)を伴う、試料中のプローブセットの発現データを使用して、モデルをトレーニングする。
本分野では、モデルの出力、ここでは、重み付け線形スコアを最適化するように、モデルトポロジーを考慮に入れ、モデルパラメータ、ここでは、重み及び閾値を変化させる、多数のトレーニングアルゴリズム(例えば、回帰)が知られている。代替的に、また、最適化アルゴリズムを必要とせずに、観察された発現レベルから直接、重みを計算することも可能である。
本明細書で「ブラックアンドホワイト」法と称する第1の方法は、各重みが、集合{−1、0、1}の要素である、三元システムへと縮約される。これを生物学的コンテキストで述べれば、−1及び1は、細胞シグナル伝達経路活性の場合、それぞれ、下方調節及び上方調節された、標的遺伝子又はプローブセットに対応する。プローブセット又は標的遺伝子が、上方調節又は下方調節されることが、統計学的に実証できない場合、プローブセット又は標的遺伝子は、0の重みを受け取る。一例では、活性の細胞シグナル伝達経路試料の発現レベルを、不活性の細胞シグナル伝達経路を伴う試料の発現レベルと対比した、左側及び右側の、2つのt検定を使用して、使用されるトレーニングデータを与えたときに、プローブ又は遺伝子が、上方調節されるのか、下方調節されるのかを決定することができる。活性試料の平均が、不活性試料の平均より、統計学的に大きい、すなわち、p値が、ある特定の閾値、例えば、0.3を下回る場合、標的遺伝子又はプローブセットは、上方調節されることが決定される。逆に、活性試料の平均が、不活性試料の平均より、統計学的に小さい場合、標的遺伝子又はプローブセットは、細胞シグナル伝達経路の活性化時に、下方調節されることが決定される。万一、最小のp値(左側又は右側)が、前述の閾値を超える場合、標的遺伝子又はプローブセットの重みは、0であると規定することができる。
本明細書で「対数オッズ」法と称する第2の方法は、対数(例えば、底をeとする)オッズ比に基づく。各標的遺伝子又は各プローブセットについてのオッズ比は、プローブセット/標的遺伝子のレベルが、対応する閾値、例えば、全てのトレーニング試料の(重み付け)中央値を上回る場合、及びこれを下回る場合の、陽性のトレーニング試料の数及び陰性のトレーニング試料の数に基づき計算する。ゼロによる除算を回避するために、疑似カウントを加算することができる。さらなる精緻化は、プローブセット/標的遺伝子レベルが、例えば、ある特定の指定された標準偏差(例えば、2logスケールで、0.25)、その観察された値の近傍に正規分布することを仮定し、閾値を上回る確率質量、確率質量これを下回る確率質量をカウントすることにより、いくぶんか、より確率的に、閾値を上回る/下回る試料をカウントすることである。本明細書に、決定論的測定値の代わりに、疑似カウントと組み合わせ、確率質量を使用して計算されたオッズ比を、「ソフトな」オッズ比と呼ぶ。
標的遺伝子発現の数学的モデル化を使用する、細胞シグナル伝達経路活性の推測に関するさらなる詳細は、Verhaegh W.ら、「Selection of personalized patient therapy through the use of knowledge−based computational models that identify tumor−driving signal transduction pathways」、Cancer Research、74巻、11号、2014年、2936〜2945頁において見出されうる。
本明細書では、遺伝子発現オムニバス(GEO、www.ncbi.nlm.nih.gov/geo/;最終アクセス日:2016年12月3日)による高悪性度漿液性乳頭状卵巣がん患者(データセット:GSE2109及びGSE9891)のセット、及び正常卵巣組織試料の対応するセット(データセット:GSE7307、GSE18520、GSE29450及びGSE36668)のそれぞれに由来する患者試料の発現について公開されているデータを使用した。高悪性度漿液性卵巣がんは、大半の症例において、活性のNotch細胞シグナル伝達経路を有することが知られているのに対し、正常卵巣組織試料は、不活性のNotch細胞シグナル伝達経路を有する。較正試料を選択する前に、データセットに対して、品質管理を実施して、試料が信頼できることを確認した。較正の目的で、Notchが最も活性である卵巣がん試料であって、全ての標的遺伝子、各個別の試料について、AffymetrixmによるRNA発現値を加算し、その後、合計値に従い、試料をランク付けすることにより決定される卵巣がん試料を、利用可能なセットから選択した。ランク付けが最高位である、20例の試料を、Notch活性であると仮定した。品質管理に合格した、12例の正常卵巣試料から、試料番号:GSM176237、GSM729048、GSM462651、GSM729050、GSM729051、GSM175789、GSM462652、GSM176131、GSM176318、GSM898306、GSM898307である、11例の試料を、Notch不活性の較正試料(1例の正常卵巣試料は、Notch活性であることが見出された)として選び出した(データセット:GSE42259に由来する試料もまた、Notch不活性の較正試料と考えられたが、品質管理の後、これらの試料のいずれも、残らなかった)。これらを使用して、Notchの活性及び不活性のそれぞれについて、モデルを較正した。較正されたモデルを、GEOデータベースに由来する、多数の公開データセットであって、Notch活性に関するグラウンドトゥルース、すなわち、Notch活性が、誘導又は阻害された細胞系(例えば、ガンマ−セクレターゼなどのNotch阻害剤で処理されるか、又はNotch3の細胞内発現を誘導する可能性を有する細胞系)を含有する公開データセット上で査定した。適用例として、モデルを、生存データが知られている乳がん試料のデータセット上で実行した。
図9は、11例の正常卵巣(群1)試料、及び20例の高悪性度漿液性乳頭状卵巣癌(群2)試料について公開されている発現データセット(データセット:GSE2109、GSE9891、GSE7307、GSE18520、GSE29450、GSE36668から採取された試料のサブセット)を使用する、表2による、18の標的遺伝子についてのショートリスト、及び本明細書で記載される方法に基づく、ベイジアンネットワークモデルの較正結果を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経路のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。モデルは、不活性の較正試料を、活性の較正試料から、明確に分離することが可能であった。
図10は、11例の正常卵巣(群1)試料、及び20例の高悪性度漿液性乳頭状卵巣癌(群2)試料について公開されている発現データセット(データセット:GSE2109、GSE9891、GSE7307、GSE18520、GSE29450、GSE36668から採取された試料のサブセット)を使用する、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)、及び本明細書で記載される方法に基づく、ベイジアンネットワークモデルの較正結果を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。ここでもまた、モデルは、不活性の較正試料を、活性の較正試料から、明確に分離することが可能であった。
以下では、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)、及び18の標的遺伝子についてのショートリストのそれぞれを使用してトレーニングされた、例示的なベイジアンネットワークモデルについてのバリデーション結果を、図11〜18に示す。
図11は、データセット:GSE6495に由来する、MOLT4細胞系の、3つの独立の培養物についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。MOLT4細胞系は、高度なNotchシグナル伝達を示すことが知られており、モデルは、これを、正確に予測した(群1)。細胞を、ガンマ−セクレターゼ阻害剤(GSI)である、5μMのDAPTで、48時間にわたり処理した(群2)。GSIは、Notchシグナル伝達を阻害することが知られており、モデルは、この群における、Notch活性の低下を正確に検出した(Dohda T.ら、「Notch signaling induces SKP2 expression and promotes reduction of p27Kip1 in T−cell acute lymphoblastic leukemia cell」、Experimental Cell Research、313巻、14号、2007年8月、3141〜3152頁を参照されたい)。
図12は、データセット:GSE6495に由来する、MOLT4細胞系の、3つの独立の培養物についての、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。MOLT4細胞系は、高度なNotchシグナル伝達を示すことが知られており、モデルは、これを、正確に予測した(群1)。細胞を、ガンマ−セクレターゼ阻害剤(GSI)である、5μMのDAPTで、48時間にわたり処理した(群2)。GSIは、Notchシグナル伝達を阻害することが知られており、モデルは、この群における、Notch活性の低下を正確に検出した(Dohda T.ら、「Notch signaling induces SKP2 expression and promotes reduction of p27Kip1 in T−cell acute lymphoblastic leukemia cell」、Experimental Cell Research、313巻、14号、2007年8月、3141〜3152頁を参照されたい)。
図13は、誘導型Notch3細胞内構築物をトランスフェクトされたIMR32細胞についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。50ng/mLのドキシサイクリンの存在下で、Notch3の細胞内発現を駆動する、2つの独立の単一細胞由来のクローン(c6、c8)が示される。t=0時間において、両方のクローンについて、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、低度のNotch活性を検出する。Notch3の細胞内発現の誘導の後、Notch活性は、両方のクローンにおいて上昇し、t=24時間において安定化することを正確に観察する。(データセット:GSE16477;van Nes J.ら、「A NOTCH3 Transcriptional Module Induces Cell Motility in Neuroblastoma」、Clinical Cancer Research、19巻、13号、2013年7月、3485〜3494頁)。
図14は、漸増用量のプラスチック固定化Notchリガンドデルタ1 ext−IgGを伴い、72時間にわたり培養されたCD34+CD45RA−Lin−HPC(データセット:GSE29524)についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、対照(群1)と比較した、デルタ1ext−IgG(群2)上で培養された細胞内の、高度のNotch活性を正確に予測する。
図15は、高Notch活性を有することが知られているCUTLL1細胞についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。ガンマ−セクレターゼ阻害剤(GSI)による処理は、Notchシグナル伝達を阻害する。データセット:GSE29544では、GSIのウォッシュアウトの2時間後、Notch活性は、高度であることが観察された。この図では、いずれの場合にも、Notch活性は、高度であることが予期されるので、非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞に由来するデータをプールする。6つの群を識別することができる:1)非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞の群である。この場合、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、この群における、高度なNotch活性を正確に予測する。2)GSIで処理されたCUTLL1細胞であり、これらについて、モデルは、低度のNotch活性を正確に予測する。3+4)Notchシグナル伝達に影響を及ぼすことが予期されない、空のMigRIレトロウイルスで処理されたCUTLL1細胞である。この場合、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、GSIのウォッシュアウト後における細胞(群3)、及びGSIで処理された細胞(群4)について、高度なNotch活性を正確に予測する。5+6)MigRI−ドミナントネガティブのMAML1ウイルスを形質導入されたCUTLL細胞である。DNMAML1は、Notchアンタゴニストであり、Notchシグナル伝達は、これらの細胞内で、低度であることが予期される。モデルは、GSIのウォッシュアウト後における細胞(群5)、及びGSIで処理された細胞(群6)のいずれについても、低度のNotch活性を正確に予測する(Wang H.ら、「Genome−wide analysis reveals conserved and divergent features of Notch1/RBPJ binding in human and murine T−lymphoblastic leukemia cells」、Proceedings of the National Academy of Sciences of the USA、108巻、36号、2011年、14908〜14913頁を参照されたい)。
図16は、高Notch活性を有することが知られているCUTLL1細胞についての、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。ガンマ−セクレターゼ阻害剤(GSI)による処理は、Notchシグナル伝達を阻害する。データセット:GSE29544では、GSIのウォッシュアウトの2時間後、Notch活性は高度であることが観察された。この図では、いずれの場合にも、Notch活性は、高度であることが予期されるので、非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞に由来するデータをプールする。6つの群を識別することができる:1)非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞の群である。この場合、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、この群における、高度なNotch活性を正確に予測する。2)GSIで処理されたCUTLL1細胞であり、これらについて、モデルは、低度のNotch活性を正確に予測する。3+4)Notchシグナル伝達に影響を及ぼすことが予期されない、空のMigRIレトロウイルスで処理されたCUTLL1細胞である。この場合、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルは、GSIのウォッシュアウト後における細胞(群3)、及びGSIで処理された細胞(群4)について、高度なNotch活性を正確に予測する。5+6)MigRI−ドミナントネガティブのMAML1ウイルスを形質導入されたCUTLL細胞である。DNMAML1は、Notchアンタゴニストであり、Notchシグナル伝達は、これらの細胞内で、低度であることが予期される。モデルは、GSIのウォッシュアウト後における細胞(群5)、及びGSIで処理された細胞(群6)のいずれについても、低度のNotch活性を正確に予測する(Wang H.ら、「Genome−wide analysis reveals conserved and divergent features of Notch1/RBPJ binding in human and murine T−lymphoblastic leukemia cells」、Proceedings of the National Academy of Sciences of the USA、108巻、36号、2011年、14908〜14913頁を参照されたい)。
図17は、COUP−TFII siRNAをトランスフェクトされたHUVEC細胞(データセット:GSE33301)についての、表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。COUP−TFIIは、Notchシグナル伝達を抑制することが知られている(You L.R.ら、「Suppression of Notch signaling by the COUP−TFII transcription factor regulates vein identity」、435巻、7038号、2005年5月、98〜104頁を参照されたい)。表2による、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、COUP−TFII siRNAをトランスフェクトされた細胞(群2)において、対照細胞(群1)と比較して高度のNotch活性を正確に検出する(Chen X.ら、「COUP−TFII is a major regulator of cell cycle and Notch signaling pathways」、26巻、8号、2012年8月、1268〜1277頁を参照されたい)。
図18は、GSE6532、GSE9195、GSE12276、GSE20685、GSE21653、及びEMTAB365に由来する試料中の乳がん亜群についての、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示す。Notch活性は、これらのデータセット内の、全ての乳がん試料中で高度であることが観察される。一元ANOVAに続いて、ゲームズ−ハウウェル事後検定を行う結果は、Basalと対比したNormL及びHER2と対比したLumAを除き、ほぼ全ての群が有意差を有することを示す。表4(亜群:Basal、HER2、LumA=Luminal A、LumB=Luminal B、NormL=正常様)を参照されたい。
表4:図18に示される、乳がん試料の異なる亜群を比較する、ゲームズ−ハウウェル事後検定の結果。<0.05のp値を、有意であると考える。
Figure 2020536530
表5は、18の標的遺伝子についてのショートリストを使用して、図18で使用されるデータセット上でトレーニングされた、例示的なベイジアンネットワークモデルのための、Notch活性についてのCox回帰の結果を示す。全ての試料を併せて、かつ、より具体的にLuminal A及びLuminal Bについては、本発明者らのモデルにより予測されるNotch活性の増大と共に、予後の大幅な増悪が見られる。これは、Notch1について試験結果が陽性である患者は、無病生存が短いことを見出した、近年の刊行物により裏付けられる(Zhong Y.ら、「NOTCH1 is a Poor Prognostic Factor for Breast Cancer and Is Associated With Breast Cancer Stem Cells」、Oncotargets and Therapy、9巻、2016年11月、6865〜6871頁を参照されたい)。
表5:表2による、18の標的遺伝子についてのショートリストを使用して、図18で使用されるデータセット上でトレーニングされた、例示的なベイジアンネットワークモデルのための、Notch活性についてのCox回帰の結果
Figure 2020536530
図19は、プラスチック固定化Notchリガンドデルタ1 ext−IgGの段階的投与を伴い、72時間にわたり培養されたCD34+CD45RA−Lin−HPC(データセット:GSE29524)についての、表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、対照(群1)と比較した、デルタ1ext−IgG(群2)上で培養された細胞内の、高度のNotch活性を正確に予測する。
図20は、高Notch活性を有することが知られているCUTLL1細胞についての、表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。略図では、垂直軸は、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれを指し示し、ここで、水平軸を上回る値は、TFエレメントが「存在する」/活性であるオッズが大きいことに対応し、水平軸を下回る値は、TFエレメントが「存在しない」/不活性であるオッズが、「存在する」/活性であるオッズより大きいことを指し示す。ガンマ−セクレターゼ阻害剤(GSI)による処理は、Notchシグナル伝達を阻害する。データセット:GSE29544では、GSIのウォッシュアウトの2時間後、Notch活性は高度であることが観察された。この図では、いずれの場合にも、Notch活性は、高度であることが予期されるので、非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞に由来するデータをプールする。6つの群を識別することができる:1)非処理CUTLL1細胞、及びGSIのウォッシュアウト後におけるCUTLL1細胞の群である。この場合、18の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、この群における、高度なNotch活性を正確に予測する。2)GSIで処理されたCUTLL1細胞であり、これらについて、モデルは、低度のNotch活性を正確に予測する。3+4)Notchシグナル伝達に影響を及ぼすことが予期されない、空のMigRIレトロウイルスで処理されたCUTLL1細胞である。この場合、表3による、12の標的遺伝子についてのショートリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルは、GSIのウォッシュアウト後における細胞(群3)、及びGSIで処理された細胞(群4)について、高度なNotch活性を正確に予測する。5+6)MigRI−ドミナントネガティブのMAML1ウイルスを形質導入されたCUTLL細胞である。DNMAML1は、Notchアンタゴニストであり、Notchシグナル伝達は、これらの細胞内で、低度であることが予期される。モデルは、GSIのウォッシュアウト後における細胞(群5)、及びGSIで処理された細胞(群6)のいずれについても、低度のNotch活性を正確に予測する(Wang H.ら、「Genome−wide analysis reveals conserved and divergent features of Notch1/RBPJ binding in human and murine T−lymphoblastic leukemia cells」、Proceedings of the National Academy of Sciences of the USA、108巻、36号、2011年、14908〜14913頁を参照されたい)。
図21は、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)、及び表3による、12の標的遺伝子についてのショートリストのそれぞれを使用してトレーニングされた、例示的なベイジアンネットワークモデルの間の相関を示す。略図では、水平軸は、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルにより予測される、活性のNotch細胞シグナル伝達経路及び不活性のNotch細胞シグナル伝達経のそれぞれに対応する、TFエレメントが「存在する」オッズ及び「存在しない」オッズのそれぞれ(log2のスケールにおける)を指し示す。垂直軸は、表3による、12の標的遺伝子についてのショートリスト(データセット:GSE5682、GSE5716、GSE6495、GSE9339、GSE14995、GSE15947、GSE16477、GSE16906、GSE18198、GSE20011、GSE20285、GSE20667、GSE24199、GSE27424、GSE29524、GSE29544、GSE29850、GSE29959、GSE32375、GSE33301、GSE33562、GSE34602、GSE35340、GSE36176、GSE37645、GSE39223、GSE42259、GSE46909、GSE49673、GSE53537、GSE54378、GSE57022、GSE61827、GSE74996、GSE81156、GSE82298)を使用してトレーニングされた、例示的なベイジアンネットワークモデルにより予測される、同じ情報を指し示す。2つのモデルは、2.2×10−16のp値、及び0.929の相関係数により、有意に相関する。
図22及び図23は、(i)7つのNotch標的遺伝子についてのリスト(DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARP)、及び10のNotch標的遺伝子についてのリスト(7つのNotch標的遺伝子に、EPHB3、SOX9、及びNFKB2を加えたリスト)、並びに(ii)8つのNotch標的遺伝子についてのリスト(DTX1、HES1、HES4、HES5、HEY2、MYC、NRARP、及びPTCRA)、及び12のNotch標的遺伝子についてのリスト(8つのNotch標的遺伝子に、HEYL、HEY1、PLXND1、及びGFAPを加えたリスト)を使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測についてのさらなる比較を示す。7つのNotch標的遺伝子は、表1〜3による標的遺伝子のリストの各々に含まれ、8つのNotch標的遺伝子は、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)だけ含まれる、さらなる標的遺伝子(PTCRA)を含む。10のNotch標的遺伝子についてのリストのうちの、3つのさらなる標的遺伝子は、表3による、12の標的遺伝子についてのショートリストから取り出され、12のNotch標的遺伝子についてのリストのうちの、4つのさらなる標的遺伝子であって、3つのさらなる標的遺伝子と異なる標的遺伝子は、表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)から取り出された。比較は、例示的に、表1〜3による標的遺伝子のリストの各々のサブセットである、7つのNotch標的遺伝子についてのリスト、及び表1による、標的遺伝子についての、証拠をキュレートしたリスト(26の標的遺伝子のリスト)のサブセットである、8つのNotch標的遺伝子ついてのリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測が、それぞれのリストに由来する、さらなる標的遺伝子を付加することにより、さらに改善されうることを示す。
図22は、10のNotch標的遺伝子についてのリストと対比した7つのNotch標的遺伝子についてのリストを使用する、Notch細胞シグナル伝達経路活性の予測についての比較を示す。モデルを、誘導型Notch3細胞内構築物をトランスフェクトされたIMR32細胞に由来する試料に対して実行した。略図では、水平軸は、時間単位の時間を指し示し、垂直軸は、相対Notch細胞シグナル伝達経路活性(log2のオッズスケールにおける)を指し示す。いずれのモデルも、Notch3細胞内構築物の誘導の後における、Notch活性の予期される増大を正確に示す。しかし、10の標的遺伝子によるモデル(点線)は、7つの標的遺伝子によるモデル(実線)と比較して、より大きな活性の増大を示す。Notch活性は、t=0時間において、0とし、比較を容易とした(データセット:GSE16477;また、van Nes J.ら、「A NOTCH3 Transcriptional Module Induces Cell Motility in Neuroblastoma」、Clinical Cancer Research、19巻、13号、2013年7月、3485〜3494頁も参照されたい)。
図23は、12のNotch標的遺伝子についてのリストと対比した8つのNotch標的遺伝子についてのリストを使用する、Notch細胞シグナル伝達経路活性の予測についての比較を示す。モデルを、Jag1レトロウイルスを感染させた子宮内膜間質細胞に由来する試料(データセット:GSE16906)に対して実行した。Jag1とは、結合すると、Notch受容体の切断を誘導し、これにより、最終的に、Notch標的遺伝子の転写を誘導するNotchリガンドである。12の標的遺伝子によるモデル(グラフの右側)は、対照細胞(図中の「C」)と、Jag1感染させた細胞(図中の「Jag1 INF」)との間の、Notch活性(log2によるオッズとして、垂直軸に示される)の、8つの標的遺伝子によるモデル(グラフの左側)と比較した、良好な分離を示す(また、Mikhailik A.ら、「Notch ligand−dependent gene expression in human endometrial stromal cells」、Biochemical and Biophysical Research Communications、388巻、3号、2009年10月、479〜482頁も参照されたい)。
較正された数学的モデル、例えば、例示的なベイジアンネットワークモデルを、マイクロアレイ又はRNAシーケンシングに由来する、入力されたmRNAデータに適用する代わりに、臨床適用では、例えば、qPCRを使用して、標的遺伝子のmRNAレベルを決定する、統合型プラットフォーム上で、試料測定を実施するための、専用アッセイを開発することが有益である。次いで、開示された標的遺伝子のRNA/DNA配列を使用して、どのプライマー及びプローブを、このようなプラットフォーム上で選択するのかを決定する。
このような専用アッセイのバリデーションは、マイクロアレイベースの数学的モデルを、基準モデルとして使用し、開発されたアッセイが、バリデーション試料のセットに対して、同様な結果をもたらすのかどうかを検証することにより行うことができる。専用アッセイに次いで、このバリデーションもまた、RNAシーケンシングデータを、入力された測定値として使用する、同様な数学的モデルを構築し、較正するために行うことができる。
較正された数学的モデル、例えば、例示的なベイジアンネットワークモデルを使用する、マイクロアレイ/RNAシーケンシングベースの探索に基づき、特異的な細胞シグナル伝達経路活性を、最もよく指し示すことが見出された、標的遺伝子のセット、例えば、表1〜3は、対象の試料及び/又はコンピュータ上で実施される、マルチプレックス定量的PCRアッセイに移行し、発現測定値を解釈し、かつ/若しくはNotch細胞シグナル伝達経路の活性を推測することができる。細胞シグナル伝達経路活性についての、このような試験(例えば、セントラルサービスラボにおける、FDAの承認若しくはCLIAの認証を受けた試験、又は探索使用だけのための、検査室により開発された試験)を開発するためには、規制機関による承認を得るのに、臨床試験における臨床の場でのバリデーションを必要とする、標準化された試験キットの開発が要求される。
本発明は、デジタル処理デバイスにより実施される、対象におけるNotch細胞シグナル伝達経路の活性を推測するための、コンピュータにより実施される方法であって、コンピュータにより実施される方法は、推測することが、対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルに基づく、コンピュータにより実施される方法に関する。本発明は、対象におけるNotch細胞シグナル伝達経路の活性を推測するための装置であって、装置は、方法を実施するように構成されたデジタルプロセッサを含む装置、対象におけるNotch細胞シグナル伝達経路の活性を推測するための非一時的記憶媒体であって、非一時的記憶媒体は、方法を実施するための、デジタル処理デバイスにより実行可能な命令を保存する非一時的記憶媒体、及び対象におけるNotch細胞シグナル伝達経路の活性を推測するためのコンピュータプログラムであって、コンピュータプログラムは、コンピュータプログラムが、デジタル処理デバイス上で実行される場合に、デジタル処理デバイスに、方法を実施させるためのプログラムコード手段を含むコンピュータプログラムにさらに関する。
方法は、例えば、Notch細胞シグナル伝達経路の(異常な)活性の診断、Notch細胞シグナル伝達経路の推測された活性に基づく予後診断、Notch細胞シグナル伝達経路の推測された活性に基づく臨床試験における対象の登録、実施される後続の試験の選択、コンパニオン診断試験の選択、臨床決定支援システムなどにおいて使用される。この点で、国際特許出願公開第WO2013/011479A2号(「標的遺伝子発現の確率モデル化を使用する、細胞シグナル伝達経路活性の評価」)、国際特許出願公開第WO2014/102668A2号(「標的遺伝子発現の線形結合を使用する、細胞シグナル伝達経路活性の評価」)、及びこれらの適用について、より詳細に記載する、Verhaegh W.ら、「Selection of personalized patient therapy through the use of knowledge−based computational models that identify tumor−driving signal transduction pathways」、Cancer Research、74巻、11号、2014年、2936〜2945頁が参照される。
[実施例4]
マウス組織についてのさらなる結果
シグナル伝達経路は、異なる種にわたり保存されていることが多く、同様な機能、及び同様な直接的標的遺伝子を有する。しかし、異なる種の間で、直接的標的遺伝子は、完全に同じではなく、遺伝子のDNA/mRNA配列も、一般に異なる。種間の遺伝子配列類似性(相同性)は、これらの種の間の、進化距離に依存し、例えば、マウスとヒトとの差違は、ヒトとトカゲとの差違より小さい。
これらの種間の類似性のために、(臓器/組織)の発生、細胞分裂、及び疾患などの生物学的過程について研究するために、動物モデルを使用することが多い。マウスは、そのヒトとの遺伝学的近接性のために、一般的なモデル生物である。例は、てんかん及びアルツハイマー病などの神経障害について研究するための、マウスモデルの使用である。このような障害については、ヒト組織を得ることが侵襲的であり(どのみち腫瘍の生検を採取することが多い、がんとは対照的に)、障害を模倣するマウスモデルが開発されている。
マウスモデルにおいて、シグナル伝達経路活性について評価することが可能であることは、それが、抽出された組織内の細胞の機能状態について、何ごとかを語りかけるので、極めて有用である。疾患マウスモデルの場合、これらのマウスモデルは通例、ヒト疾患を、可能な最良の形で反映するように作出されるので、シグナル伝達経路活性は、疾患のヒト変化形についての情報をもたらしうる。
Notch細胞シグナル伝達経路モデルは、元来、ヒト組織のために開発された、すなわち、表1〜3において選択された標的遺伝子は、ヒトにおける直接的標的遺伝子であり、モデルのための入力は、ヒトmRNA(例えば、マイクロアレイ実験、qPCR実験、又はRNAseq実験)の発現レベルであり、較正は、ヒト試料に由来する発現データ上で行われる。
本明細書ではまた、マウスにおける使用のための、Notch細胞シグナル伝達経路モデルも示される。マウスにおける直接的Notch細胞シグナル伝達経路の標的遺伝子を選択し、適切な較正試料(公開データベースに由来するAffymetrixマイクロアレイデータ)を使用することにより、マウスmRNA発現レベルを入力として使用し、Notch細胞シグナル伝達経路の活性をこの入力から推測するモデルを創出した。次いで、独立試料(公開データベースに由来するAffymetrixマイクロアレイデータ)を使用して、このモデルが、マウスにおけるNotch細胞シグナル伝達経路の活性を正確に測定することを示すため、このモデルをバリデーションした。
マウスNotch細胞シグナル伝達経路モデルのための直接的標的遺伝子の選択は、既に記載したのと同様な形で行った。ヒトNotchモデルのために使用した、26の遺伝子リストを、出発点として使用した。このリストを、証拠スコア(これは、既に記載した通りに計算した)によりランク付けし、「“mouse”AND“direct target gene”」などの検索キーワード、及びヒト直接的標的遺伝子について、既に見出された文献に由来する基準を使用して、ランク付けが首位の遺伝子について、文献検索を実施した。
まず、遺伝子が、マウスにおいて、実際に存在することを確認し、次いで、遺伝子がまた、マウスにおける、直接的Notch標的遺伝子でもあることを確認した。これは、ヒト標的遺伝子について使用されたのと同様な証拠(すなわち、転写因子複合体結合部位の存在、ChIP、ルシフェラーゼアッセイ、示差的発現、GSI処理など、実験の証拠)を使用して行った。証拠の複数の典拠が見出された場合は、遺伝子を、マウスNotchについての直接的標的遺伝子として受け入れた。このようにして、表6において示される通り、Notchマウスモデルについて、10の直接的標的遺伝子の選択を行った。
表6: Notch標的遺伝子についての、証拠をキュレートしたリストに基づく、Notch標的遺伝子の「10の標的遺伝子についてのマウスリスト」(Affymetrixマウスゲノム430 2.0アレイに由来する)
Figure 2020536530
Notchマウスモデルを、GEO(遺伝子発現オムニバス)データベースに由来する公開データセットである、データセットGSE15268に由来する試料上で較正した。このデータセットは、NotchIC(Notch細胞内ドメイン)誘導型構築物(ヒドロキシタモキシフェン(OHT)の添加により誘導される)を伴うマウス胚性幹細胞に由来するAffymetrixマイクロアレイデータを含有する。このデータセットから、NotchICを誘導しなかった4例の試料を、Notch不活性試料(GSM381312、GSM381313、GSM381317、GSM381316)として使用し、OHTを添加することにより、Notchを誘導した4例の試料を、Notch活性試料(GSM381324、GSM381325、GSM381320、GSM381321)として使用した。
次いで、較正されたNotchマウスモデルを、いくつかのデータセット:較正セット上、及びいくつかの独立バリデーションセット上で実行して、モデルが、Notch活性試料を、Notch不活性試料から、完全に区別するのに成功しうることを示した。これらの結果を、図24〜図27に示す。
図24は、2例の対照胚性幹細胞(図中のC ESc)、2例の対照中胚葉前駆細胞(図中のC MPc)、OHT処理されていないタモキシフェン誘導型NERT構築物(NotchIC)を含有する2例のESc試料(図中の「NERT ESc、OHTなし」)、OHT処理されたタモキシフェン誘導型NERT構築物(NotchIC)を含有する2例のESc試料(図中の「NERT ESc、OHT」)、OHT処理されていないタモキシフェン誘導型NERT構築物(NotchIC)を含有する4例のMPc試料(図中の「NERT MPc、OHTなし」)、及びOHT処理されたタモキシフェン誘導型NERT構築物(NotchIC)を含有する4例のMPc試料(図中の「NERT MPc、OHT」)を含有する、公開されている発現データセットであるGSE15268を使用する表6による10の標的遺伝子についてのマウスリスト、及び本明細書で記載される方法に基づく、ベイジアンモデルの較正結果を示す。モデルは、不活性の較正試料(対照ES及び対照MPc)を、活性の較正試料(NERT MPc、OHT)から、明確に分離することが可能であった。データセット内の他の試料もまた、正確に分離された(また、Meier−Stiegen F.ら、「Activated Notch1 Target Genes during Embryonic Cell Differentiation Depend on the Cellular Context and Include Lineage Determinants and Inhibitors」、PLoS One、5巻、7号、2010年7月も参照されたい)。
図25は、構成的に活性の誘導型Notch1細胞内ドメイン(NICD1)を伴うマウス乳腺(データセット:GSE51628)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。NICD1が誘導されない乳腺試料(図中の「Mg」)について、Notchマウスモデル(10の標的遺伝子)は、低度のNotch活性を検出する。予期される通り、ドキシサイクリンを使用して、NICD1が誘導される乳腺試料(図中の「Mg、NICD1 a」)は、大幅に高度のNotch活性を正確に示す。この図では、48時間後及び96時間後の時点が組み合わされている(また、Abravanel D.L.ら、「Notch promotes recurrence of dormant tumor cells following HER2/neu−targeted therapy」、Journal of Clinical Investigation、 125巻、6号、2015年6月、2484〜2496頁も参照されたい)。
図26は、Notch1を活性化させる、条件付きトランスジェニック系を伴うマウス卵黄嚢組織、及びRBPJ(Notch転写因子複合体の一部)機能喪失を伴うトランスジェニックマウスに由来するマウス卵黄嚢組織(データセット:GSE22418)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。野生型試料(図中の「Wt」)及びRBPJ機能喪失試料(図中の「RBPJ l−o−f」)のいずれも、低度のNotch活性を示し、Notch1が活性化させられる、卵黄嚢組織に由来する試料(図中の「Notch1 a」)は、予期される通り、Notch活性の上昇を示す(また、Copeland J.N.ら、「Notch signaling regulates remodeling and vessel diameter in the extraembryonic yolk sac」、BMC Developmental Biology、2011年2月も参照されたい)。
図27は、Notch2受容体の、条件付き機能獲得対立遺伝子を伴うマウス骨髄細胞(成体骨髄赤血球前駆細胞)(データセット:GSE46724)についての、表6による、10の標的遺伝子についてのマウスリストを使用してトレーニングされた、例示的なベイジアンネットワークモデルによる、Notch細胞シグナル伝達経路活性の予測を示す。マウスNotchモデル(10の標的遺伝子)は、ICN2陽性(細胞内Notch2)試料(図中の「ICN2 p」)について、ICN2陰性試料(図中の「ICN2 n」)と比較して、高度のNotch活性を正確に計算する(また、Oh P.ら、「In vivo mapping of notch pathway activity in normal and stress hematopoiesis」、Cell Stem Cell、13巻、1号、2013年8月、190〜204頁も参照されたい)。
[実施例5]
本発明を例示するためのさらなる情報
(1)遺伝子発現のレベルの測定
本明細書で記載される方法を使用して、Notch細胞シグナル伝達経路の活性を推測するのに、本明細書で記載される標的遺伝子の固有のセットから導出されるデータを、さらに活用する。
抽出された試料中の遺伝子発現レベルを解析するための方法は、一般に知られている。例えば、方法など、ノーザンブロット法、PCR、nestedPCR、定量的リアルタイムPCR(qPCR)、RNA−seq、又はマイクロアレイの使用は、いずれも、遺伝子発現レベルのデータを導出するのに使用することができる。本明細書では、標的遺伝子の遺伝子発現を解析するための、当技術分野で知られている全ての方法が想定される。
PCRベースの方法を使用して、遺伝子の発現産物を決定する方法は、特に、有用である。PCRを使用して、遺伝子発現のレベルを定量するためには、各増幅サイクルの後における、リアルタイムのPCR産物の蓄積を測定するのに、従来の定量的リアルタイムPCR(qPCR)を使用して、関心のある各PCR産物の量を推定することが典型的である。これは、典型的に、挿入色素、副溝結合色素、又は光の適用が、蛍光を発するように、レポーターを励起する、蛍光発生プローブなど、検出可能なレポーターを活用し、結果として得られる蛍光は、参照により本明細書に組み込まれる、米国特許第6,713,297号において開示されているCCDカメラ又は光電子増倍検出システムなどの、CCDカメラ又は光電子増倍検出システムを使用して検出されることが典型的である。
一部の実施形態では、定量的リアルタイムPCR(qPCR)アッセイにおいて、PCR産物の検出に使用されるプローブは、蛍光マーカーを含みうる。多数の蛍光マーカーが市販されている。例えば、Molecular Probes,Inc.(Eugene、Oreg.)は、多種多様な蛍光色素を販売している。非限定例は、Cy5、Cy3、TAMRA、R6G、R110、ROX、JOE、FAM、Texas Red(商標)、及びOregon Green(商標)を含む。さらなる蛍光マーカーは、qPCRアッセイにおける、従来の5’加水分解型プローブを伴う、IDT ZEN Double−Quenched Probesを含みうる。これらのプローブは、例えば、3’ TAMRA Quencher、3’ Black Hole Quencher(BHQ、Biosearch Technologies)、又はインターナル型のZEN Quencher及び3’ Iowa Black Fluorescent Quencher(IBFQ)を伴う5’FAM色素を含有しうる。
当技術分野でよく知られている方法を使用して、本発明に従う、有用な蛍光色素を、オリゴヌクレオチドプライマーへと接合させることができる。例えば、蛍光標識を、オリゴヌクレオチドへと付加するための、一般的な方式は、色素のN−ヒドロキシスクシンイミド(NHS)エステルを、標的上の反応性アミノ基と反応させることである。ヌクレオチドは、例えば、ヌクレオ塩基上の、アリルアミン基の組入れにより、反応性アミノ基を保有するように修飾することができる。アリルアミンを介する標識付けについては、例えば、参照により本明細書に組み込まれる、米国特許第5,476,928号及び同第5,958,691号において記載されている。ヌクレオチド、オリゴヌクレオチド、及びポリヌクレオチドを、蛍光標識付けする他の手段は、当業者によく知られている。
他の蛍光発生手法は、参照により本明細書に組み込まれる、米国特許第5,436,134号及び同第5,658,751号において開示される通り、任意の遺伝子発現産物から増幅されたDNAと共に挿入されると、蛍光を発する、SYBR−green色素など、一般的な検出システムの使用を含む。
標的遺伝子発現レベルを決定するための、別の有用な方法は、異なる整理学的条件の間の、遺伝子発現レベルの差違、又は発生時に、若しくは疾患進行の経過にわたり生じる変化を含む、トランスクリプトーム解析のために使用される、強力な解析ツールである、RNA−seqを含む。
遺伝子発現レベルを決定するための、別の手法は、マイクロアレイ、例えば、当技術分野でよく知られている、RNAマイクロアレイ及びDNAマイクロアレイの使用を含む。マイクロアレイを使用して、多数の遺伝子の発現を、同時に定量することができる。
(2)Notch細胞シグナル伝達の活性を決定するための、一般的ワークフロー
Notch細胞シグナル伝達の活性を、対象から単離された試料から推測するための工程を例示的に図解するフローチャートを、図3に示す。第1に、試料からmRNAを単離する(11)。第2に、当技術分野で知られている、遺伝子発現を測定するための方法を使用して、本明細書で記載される、少なくとも3つ以上のNotch標的遺伝子の固有のセットのmRNA発現レベルを測定する(12)。次に、3つ以上のNotch標的遺伝子の発現レベルを、Notch TFエレメントの活性レベルと関連付ける、較正された数学的経路モデル(14)を使用して、Notch転写因子(TF)エレメントの活性レベル(13)を決定する。最後に、対象の試料中の、Notch TFエレメントの、決定された活性レベルに基づき、対象におけるNotch細胞シグナル伝達経路の活性を推測する(15)。例えば、Notch細胞シグナル伝達経路を、活性が、ある特定の閾値を上回れば、活性であると決定し、活性が、ある特定の閾値を下回れば、不活性であると類別することができる。
(3)較正された数学的経路モデル
本明細書で想定される通り、本明細書で記載される、3つ以上のNotch標的遺伝子の固有のセットの発現レベルは、本明細書でさらに記載される、較正された数学的経路モデルを使用して、Notch TFエレメントの活性レベルを決定するのに使用される。較正された数学的経路モデルは、3つ以上のNotch標的遺伝子の発現レベルを、Notch TFエレメントの活性レベルへと関連付ける。
本明細書で想定される通り、較正された数学的経路モデルは、数学的経路モデルの適用に基づく。例えば、較正された数学的経路モデルは、確率モデル、例えば、ベイジアンネットワークモデル、又は線形モデル若しくは疑似線形モデルに基づきうる。
ある実施形態では、較正された数学的経路モデルは、Notch TFエレメントと、3つ以上のNotch標的遺伝子の発現レベルとを関連付ける、条件付き確率関係を組み込む確率モデルである。ある実施形態では、確率モデルは、ベイジアンネットワークモデルである。
代替的実施形態では、較正された数学的経路モデルは、線形モデルの場合もあり、疑似線形モデルの場合もある。ある実施形態では、線形モデル又は疑似線形モデルは、本明細書でさらに記載される、線形結合モデル又は疑似線形結合モデルである。
較正された数学的経路モデルを作出するための工程を、図4を例示的に図解するためのフローチャートに示す。最初のステップとして、mRNA発現レベルのトレーニングデータを、収集し、正規化する。データは、例えば、当技術分野で知られている、マイクロアレイプローブセットの強度(101)、リアルタイムPCRのCq値(102)、RNAseqの生リード(103)、又は代替的測定モダリティー(104)を使用して収集することができる。次いで、発現レベルの生データを、各方法に応じて、それぞれ、正規化アルゴリズム、例えば、フローズンロバストマルチアレイ解析(fRMA)又はMAS5.0を使用する正規化(111)、基準遺伝子の平均Cqに照らした正規化(112)、リードの、マッピングされたリード100万当たりの転写物1キロベース当たりのリード/フラグメント量(RPKM/FPKM)への正規化(113)、又は基準遺伝子/タンパク質に照らした正規化(114)を使用する正規化により正規化することができる。この正規化手順は、各方法に応じて、それぞれ、正規化されたプローブセットの強度(121)、正規化されたCq値(122)、正規化されたRPKM/FPKM(123)、又は正規化された測定値(124)をもたらし、これらは、トレーニング試料中の標的遺伝子発現レベルを指し示す。
トレーニングデータを正規化すると、1つ又は複数のトレーニング試料ID(131)を得、これらの特異的試料のトレーニングデータを、遺伝子発現を決定するための方法のうちの1つから得る(132)。トレーニング試料からの、最終的な遺伝子発現結果を、トレーニングデータ(133)として出力する。多様なトレーニング試料に由来するデータ(例えば、閾値、例えば、確率的ネットワーク又はベイジアンネットワークの場合のCPT、例えば、線形モデル又は疑似線形モデルの場合の重みなどを含む)の全てを組み込んで、モデルを較正する(144)。加えて、経路の標的遺伝子及び測定ノード(141)を使用して、例えば、図2に記載されるモデル構造を作出する(142)。次いで、結果として得られる、経路のモデル構造(143)を、トレーニングデータ(133)と共に組み込んで、モデルを較正するが(144)、ここで、標的遺伝子の遺伝子発現レベルは、転写因子エレメントの活性を示す。トレーニング試料中の、TFエレメントの決定の結果として、較正された経路モデル(145)を作出し、これにより、トレーニング試料中の標的遺伝子発現レベルに基づき、その後検討される関心のある試料、例えば、がんを伴う対象に由来する試料のNotch細胞シグナル伝達経路活性に割り当てる。
(4)TFエレメントの決定
TFエレメントの活性レベルを決定するための工程を例示的に図解するフローチャートを、図5に示す。対象から抽出された試料に由来する発現レベルのデータ(被験データ)(163)を、較正された数学的経路モデル(145)へと入力する。数学的経路モデルは、確率モデル、例えば、ベイジアンネットワークモデル、線形モデル、又は疑似線形モデルである。
数学的経路モデルは、Notch TFエレメントと、対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ若しくはこれを超える標的遺伝子の発現レベルとを関連付ける条件付き確率に基づく確率モデル、例えば、ベイジアンネットワークモデルであるか、又は対象の試料中で測定される、Notch細胞シグナル伝達経路の3つ若しくはこれを超える標的遺伝子の発現レベルの1つ若しくは複数の線形結合に基づく。特に、Notch細胞シグナル伝達経路の活性の決定は、その内容が、それらの全体において本明細書に組み込まれる、国際特許出願公開第WO2013/011479A2号(「標的遺伝子発現の確率モデル化を使用する、細胞シグナル伝達経路活性の評価」)において開示される通りに実施される。略述すると、データを、ベイジアンネットワーク(BN)推測エンジン判定(例えば、BNTツールボックス)へと入力する(154)。これは、BN内の全てのノードの、計算された周辺BN確率についての値のセット(155)をもたらす。これらの確率から、転写因子(TF)ノードの確率を決定し(156)、TFエレメントの活性レベル(157)を確立する。
代替的に、数学的モデルは、線形モデルである。例えば、線形モデルは、その内容が、それらの全体において本明細書に組み込まれる、国際特許出願公開第WO2014/102668A2号(「標的遺伝子発現の線形結合を使用する、細胞シグナル伝達経路活性の評価」)において記載される通りに使用されうる。標的遺伝子発現の数学的モデル化を使用して、細胞シグナル伝達経路活性を計算/決定することに関する、さらなる詳細はまた、Verhaegh W.ら、「Selection of personalized patient therapy through the use of knowledge−based computational models that identify tumor−driving signal transduction pathways」、Cancer Research、74巻、11号、2014年、2936〜2945頁においても見出されうる。略述すると、データを、計算された重み付け線形結合スコア(w/c)へと入力する(151)。これは、計算された重み付け線形結合スコアのための値のセット(152)をもたらす。これらの重み付け線形結合スコアから、転写因子(TF)ノードの重み付け線形結合スコアを決定し(153)、TFエレメントの活性レベル(157)を確立する。
(5)離散可観測量のための手順
対象におけるNotch細胞シグナル伝達経路の活性を、離散可観測量として推測するための工程を、図6を例示的に図解するためのフローチャートに示す。まず、被験試料を抽出し、被験試料ID(161)を与える。次に、mRNAの発現レベルについての被験データを収集し、正規化する(162)。被験データは、図5で、トレーニング試料について論じたのと同じ方法を使用して、マイクロアレイプローブセットの強度(101)、リアルタイムPCRのCq値(102)、RNAseqの生リード(103)、又は代替的測定モダリティー(104)を使用して収集することができる。次いで、発現レベルの生データを、各方法に応じて、それぞれ、アルゴリズム、例えば、fRMA又はMAS5.0を使用する正規化(111)、基準遺伝子の平均Cqに照らした正規化(112)、リードの、RPKM/FPKMへの正規化(113)、及び基準遺伝子/タンパク質に照らした正規化(114)を使用する正規化により正規化することができる。この正規化手順は、各方法に応じて、それぞれ、正規化されたプローブセットの強度(121)、正規化されたCq値(122)、正規化されたRPKM/FPKM(123)、又は正規化された測定値(124)をもたらす。
被験データを正規化すると、結果として得られる被験データ(163)を、較正された数学的経路モデル(145)に基づく閾値化ステップ(164)において解析する結果として、閾値化された被験データ(165)をもたらす。1つの非限定的な例では、離散可観測量を使用すると、ある特定の閾値を上回るあらゆる発現には、例えば、1の値を与え、閾値を下回る値には、0の値を与えるか、又は、代替的実施形態では、本明細書で記載される閾値を上回る確率質量を、閾値化された値として使用する。較正された数学的経路モデルに基づくと、この値は、TFエレメントの活性レベル(157)を表し、次いで、これを使用して、細胞シグナル伝達経路の活性を計算する(171)。最終的な出力は、対象における、細胞シグナル伝達経路の活性(172)をもたらす。
(6)連続可観測量のための手順
対象におけるNotch細胞シグナル伝達経路の活性を、連続可観測量として推測するための工程を、図7を例示的に図解するためのフローチャートに示す。まず、被験試料を抽出し、被験試料ID(161)を与える。次に、mRNAの発現レベルについての被験データを収集し、正規化する(162)。被験データは、図5で、トレーニング試料について論じたのと同じ方法を使用して、マイクロアレイプローブセットの強度(101)、リアルタイムPCRのCq値(102)、RNAseqの生リード(103)、又は代替的測定モダリティー(104)を使用して収集することができる。次いで、発現レベルの生データを、各方法に応じて、それぞれ、アルゴリズム、例えば、fRMA(111)、基準遺伝子の平均Cqに照らした正規化(112)、リードの、RPKM/FPKMへの正規化(113)、又は基準遺伝子/タンパク質に照らした正規化(114)を使用する正規化により正規化することができる。この正規化手順は、各方法に応じて、それぞれ、正規化されたプローブセットの強度(121)、正規化されたCq値(122)、正規化されたRPKM/FPKM(123)、又は正規化された測定値(124)をもたらす。
被験データを正規化すると、結果として得られる被験データ(163)を、較正された数学的経路モデル(145)により解析する。1つの非限定的な例として、連続可観測量を使用すると、本明細書で、さらに詳細に記載される通り、シグモイド関数を使用して、発現レベルを、0〜1の間の値へと転換する。本明細書で記載されるTFエレメントの決定を使用して、被験データを、較正された数学的経路モデルと組み合わせて解釈すると、結果として得られる値は、TFエレメントの活性レベル(157)を表し、次いで、これを使用して、細胞シグナル伝達経路の活性(171)を計算する。最終的な出力は、対象における、細胞シグナル伝達経路の活性(172)をもたらす。
(7)標的遺伝子発現レベルの決定手順
標的遺伝子発現レベルを、対象から抽出された試料から導出するための工程を、例示的に図解するフローチャートを、図8に示す。例示的な実施形態では、試料を、検査室で受領し、登録する。試料は、例えば、ホルマリン固定パラフィン包埋(FFPE)試料(181)又は新鮮凍結(FF)試料(180)を含みうる。FF試料は、直接溶解させることができる(183)。FFPE試料については、プロテイナーゼK添加時における、加熱インキュベーションステップにより、パラフィンを除去することができる(182)。次いで、細胞を溶解させ(183)、これにより、細胞膜及び核膜を破壊し、これにより、核酸(NA)を、さらなる加工のために利用可能とする。核酸を、例えば、ビーズ又はフィルターでありうる固相へと結合させる(184)。次いで、核酸を、洗浄緩衝液で洗浄して、溶解の後で存在する、全ての細胞破砕物を除去する(185)。次いで、清浄な核酸を、溶離緩衝液により、固相から引き離す(186)。DNAを、DNアーゼ処理により除去して、RNAだけが、試料中に存在することを確保する(187)。次いで、核酸試料を、RT−qPCR試料ミックス中で、直接使用することができる(188)。RT−qPCR試料ミックスは、RNA試料、cDNAを、RNA試料から調製するRT酵素、及びcDNAを増幅するPCR酵素、酵素の機能を確認する緩衝液を含有し、潜在的に、一定の濃縮物容量を設定するための分子グレードの水を含有しうる。次いで、試料ミックスを、乾燥RT−qPCRアッセイを含有する、マルチウェルプレート(すなわち、96ウェル又は384ウェルプレート)へと添加することができる(189)。次いで、指定されたプロトコールに従い、PCRマシン内で、RT−qPCRを作動させることができる(190)。例 PCRプロトコールは、i)50℃で、30分間;ii)95℃で、5分間;iii)95℃で、15秒間;iv)60℃で、45秒間;v)ステップiii及びivを反復する50サイクルを含む。次いで、生データに対して、二次導関数法を使用することにより、Cq値を決定する(191)。解析のために、Cq値をエクスポートする(192)。
(8)Notchに媒介される疾患及び障害、並びに処置の方法
本明細書で想定される通り、本発明の方法及び装置を活用して、対象、例えば、疾患若しくは障害を有することが疑われるか、又はこれを有する対象における、Notch細胞シグナル伝達経路活性について評価することができ、ここで、Notchシグナル伝達経路の状態は、疾患の存在又は進行の、全体的又は部分的な証拠となる。ある実施形態では、本明細書では、対象を処置する方法であって、処置する方法は、本明細書で記載される方法を使用して、対象から抽出された試料から導出される、Notch細胞シグナル伝達経路の活性状態に関する情報を受信するステップと、Notch細胞シグナル伝達経路の活性に関する情報が、活性のNotchシグナル伝達経路を示す場合、対象へと、Notch阻害剤を投与するステップとを含む、処置する方法が提示される。特定の実施形態では、Notch細胞シグナル伝達経路活性の指数は、Notch細胞シグナル伝達経路が活性であるオッズについての、10:1、5:1、4:1、2:1、1:1、1:2、1:4、1:5、1:10のカットオフ値に設定される。
本発明で使用されるNotch阻害剤は、よく知られている。Notch阻害剤の例は、DAPT、PF−03084014、MK−0752、RO−4929097、LY450139、BMS−708163、LY3039478、IMR−1、ジベンザゼピン、LY411575、FLI−06を含むがこれらに限定されない。
Notch経路は、多数の疾患において、役割を果たし、とりわけ、異なる種類の新生物、例えば、癌腫、肉腫、及び血液悪性腫瘍、免疫媒介性疾患、変性疾患、炎症性疾患、感染性疾患において、役割を果たす。これらの疾患は、胚系統由来の臓器又は組織に従い類別することができ、これらの臓器又は組織、例えば、脳、乳腺、皮膚、食道、消化管、血液(血液学的)、卵巣などにおいて、主に生じる。
特定の実施形態では、対象は、乳がん、肺がん、結腸がん、膵臓がん、脳がん、血液学的がん、卵巣がんを患っているか、又はこれらを患っていることが疑われる。特定の実施形態では、対象は、乳がんを患っているか、又はこれを患っていることが疑われる。
別の特定の実施形態では、対象は、脳がん、又はより好ましくは、神経芽細胞腫性がんを患っているか、又はこれを患っていることが疑われる。別の特定の実施形態では、対象は、血液がん、又はより好ましくは、T細胞リンパ芽球性白血病を患っているか、又はこれを患っていることが疑われる。
本出願は、いくつかの好ましい実施形態について記載する。前出の詳細な記載を読み、理解すれば、他の研究者らは、改変及び変更に想到する。本出願は、それらが、付属の特許請求の範囲又はそれらの同等物の範囲内に収まる限りにおいて、全てのこのような改変及び変更を含むものとみなされる。
開示される実施形態に対する、他の変動も、図面、本開示、及び付属の特許請求の範囲の検討から、当業者により理解され、特許請求された本発明の実施においてなされうる。
特許請求の範囲では、「〜を含むこと」という語は、他の要素又はステップを除外せず、単数形は、複数形を除外しない。
単一のユニット又はデバイスは、特許請求の範囲において列挙された、いくつかの項目の機能を果たす。ある特定の手段が、互いに異なる従属請求項内に列挙されるという事実だけで、これらの手段の組合せを使用しても有利でありえないことが指し示されるわけではない。
1つ又はいくつかのユニット又はデバイスにより実施される、危険性スコアの決定などの計算は、他の任意の数のユニット又はデバイスにより実施することができる。
コンピュータプログラムは、光学記憶媒体又は半導体媒体など、適切な媒体上に保存/分配され、他のハードウェアと併せて、又はこの一部として供給されるが、また、インターネット又は他の有線若しくは無線の通信システムなど、他の形態でも供給される。
[実施例6]
本出願において使用される配列表
配列表:
配列番号 遺伝子:
配列番号1 CD28
配列番号2 CD44
配列番号3 DLGAP5
配列番号4 DTX1
配列番号5 EPHB3
配列番号6 FABP7
配列番号7 GFAP
配列番号8 GIMAP5
配列番号9 HES1
配列番号10 HES4
配列番号11 HES5
配列番号12 HES7
配列番号13 HEY1
配列番号14 HEY2
配列番号15 HEYL
配列番号16 KLF5
配列番号17 MYC
配列番号18 NFKB2
配列番号19 NOX1
配列番号20 NRARP
配列番号21 PBX1
配列番号22 PIN1
配列番号23 PLXND1
配列番号24 PTCRA
配列番号25 SOX9
配列番号26 TNC


Claims (12)

  1. デジタル処理デバイスにより実施される、対象におけるNotch細胞シグナル伝達経路の活性を推測するための、コンピュータにより実施される方法であって、コンピュータにより実施される方法は、前記推測することが、
    前記対象の試料中で測定される、前記Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを受信するステップと、
    前記対象の前記試料中の3つ以上のNotch標的遺伝子の転写を制御する、Notch転写因子(TF)エレメントの活性レベルを決定するステップであって、前記決定は、前記3つ以上のNotch標的遺伝子の前記発現レベルを、前記Notch TFエレメントの前記活性レベルと関連付ける、較正された数学的経路モデルを査定することに基づく、前記決定するステップと、
    前記対象における前記Notch細胞シグナル伝達経路の前記活性を、前記対象の前記試料中の前記Notch TFエレメントの決定された活性レベルに基づき推測するステップと
    を有し、
    前記3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、ここで、2つ以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択される、
    コンピュータにより実施される方法。
  2. 前記Notch細胞シグナル伝達経路が前記対象において異常に作動しているのかどうかを、前記対象における前記Notch細胞シグナル伝達経路の推測された活性に基づき決定するステップ
    をさらに有する、請求項1に記載のコンピュータにより実施される方法。
  3. 前記Notch細胞シグナル伝達経路の異常な作動を是正する薬物を、前記対象に処方することを推奨するステップ
    をさらに有し、
    前記推奨は、前記Notch細胞シグナル伝達経路が、前記対象において異常に作動すると、前記Notch細胞シグナル伝達経路の推測された前記活性に基づき決定される場合に実施される、請求項2に記載のコンピュータにより実施される方法。
  4. 前記Notch細胞シグナル伝達経路の前記異常な作動は、前記Notch細胞シグナル伝達経路が前記対象において腫瘍プロモーターとして作動する作動である、請求項2又は3に記載のコンピュータにより実施される方法。
  5. 以下のアクティビティ:
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく診断;
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく予後診断;
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく薬物処方;
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく薬物有効性の予測;
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく有害作用の予測;
    薬物有効性のモニタリング;
    薬物の開発;
    アッセイの開発;
    経路の探索;
    がんの病期分類;
    前記対象における前記Notch細胞シグナル伝達経路の推測された前記活性に基づく臨床試験における前記対象の登録;
    実施される後続の試験の選択;及び
    コンパニオン診断試験の選択
    のうちの少なくとも1つにおいて使用される、請求項1から4のいずれか一項に記載のコンピュータにより実施される方法。
  6. 較正された前記数学的経路モデルが、前記Notch TFエレメントの前記活性レベルと、前記3つ若しくはこれを超えるNotch標的遺伝子の前記発現レベルとを関連付ける条件付き確率に基づく確率モデル、好ましくは、ベイジアンネットワークモデルであるか、又は前記数学的経路モデルが、前記3つ若しくはこれを超えるNotch標的遺伝子の前記発現レベルの1つ若しくは複数の線形結合に基づく、請求項1から5のいずれか一項に記載のコンピュータにより実施される方法。
  7. 対象におけるNotch細胞シグナル伝達経路の活性を推測するための装置であって、前記装置は、請求項1から6のいずれか一項に記載のコンピュータにより実施される方法を実施するように構成されたデジタルプロセッサを含む、装置。
  8. 対象におけるNotch細胞シグナル伝達経路の活性を推測するための非一時的記憶媒体であって、前記非一時的記憶媒体は、請求項1から6のいずれか一項に記載のコンピュータにより実施される方法を実施するための、デジタル処理デバイスにより実行可能な命令を保存する、非一時的記憶媒体。
  9. 対象におけるNotch細胞シグナル伝達経路の活性を推測するためのコンピュータプログラムであって、前記コンピュータプログラムは、前記コンピュータプログラムがデジタル処理デバイス上で実行される場合に、前記デジタル処理デバイスに、請求項1から6のいずれか一項に記載のコンピュータにより実施される方法を実施させるためのプログラムコード手段を含む、コンピュータプログラム。
  10. 対象の試料中の前記Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを測定するためのキットであって、前記キットは、
    前記3つ以上のNotch標的遺伝子を対象とするポリメラーゼ連鎖反応用のプライマーと、
    前記3つ以上のNotch標的遺伝子を対象とするプローブと、
    請求項7に記載の装置、請求項8に記載の非一時的記憶媒体、又は請求項9に記載のコンピュータプログラムと
    を含み、
    前記3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、ここで、2つ以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択される、
    キット。
  11. 対象におけるNotch細胞シグナル伝達経路の活性を推測するためのキットであって、前記キットは、
    前記対象の試料中の前記Notch細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを決定するための1つ又は複数の構成要素と、
    請求項7に記載の装置、請求項8に記載の非一時的記憶媒体、又は請求項9に記載のコンピュータプログラムと
    を含み、
    前記1つ又は複数の構成要素が、好ましくは、DNAアレイチップ、オリゴヌクレオチドアレイチップ、タンパク質アレイチップ、抗体、複数のプローブ、例えば、標識されたプローブ、RNA逆転写酵素シーケンシング構成要素のセット、及び/又はcDNA、増幅プライマーを含む、RNA若しくはDNAからなる群から選択され、
    ここで、前記3つ以上のNotch標的遺伝子が、DTX1、EPHB3、HES1、HES4、HES5、HEY2、MYC、NFKB2、NRARP、PIN1、PLXND1、及びSOX9からなる群から選択され、ここで、2つ以上のNotch標的遺伝子が、DTX1、HES1、HES4、HES5、HEY2、MYC、及びNRARPからなる群から選択され、1つ又は複数のNotch標的遺伝子が、EPHB3、NFKB2、PIN1、PLXND1、及びSOX9からなる群から選択される、
    キット。
  12. 請求項1から6のいずれか一項に記載のコンピュータにより実施される方法の実施における、請求項10又は11に記載のキットの使用。

JP2020518608A 2017-10-02 2018-09-28 標的遺伝子発現の数学的モデル化を使用する、Notch細胞シグナル伝達経路活性の評価 Withdrawn JP2020536530A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP17194288.1 2017-10-02
EP17194288.1A EP3462349A1 (en) 2017-10-02 2017-10-02 Assessment of notch cellular signaling pathway activity using mathematical modelling of target gene expression
PCT/EP2018/076488 WO2019068585A1 (en) 2017-10-02 2018-09-28 EVALUATION OF THE NOTCH CELL SIGNALING PATH ACTIVITY USING THE MATHEMATICAL MODELING OF THE TARGET GENE EXPRESSION

Publications (1)

Publication Number Publication Date
JP2020536530A true JP2020536530A (ja) 2020-12-17

Family

ID=60119805

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020518608A Withdrawn JP2020536530A (ja) 2017-10-02 2018-09-28 標的遺伝子発現の数学的モデル化を使用する、Notch細胞シグナル伝達経路活性の評価

Country Status (8)

Country Link
US (1) US20190100790A1 (ja)
EP (2) EP3462349A1 (ja)
JP (1) JP2020536530A (ja)
CN (1) CN111183233A (ja)
AU (1) AU2018344383A1 (ja)
BR (1) BR112020006484A2 (ja)
CA (1) CA3076635A1 (ja)
WO (1) WO2019068585A1 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3066107A1 (en) 2016-06-03 2017-12-07 Hough Ear Institute Combination therapies for inner ear sensory hair cell regeneration/replacement
EP3835433A1 (en) * 2019-12-12 2021-06-16 Koninklijke Philips N.V. Notch signaling pathway activity as prognostic marker in bladder cancer
EP3882363A1 (en) 2020-03-17 2021-09-22 Koninklijke Philips N.V. Prognostic pathways for high risk sepsis patients
EP3978628A1 (en) 2020-10-01 2022-04-06 Koninklijke Philips N.V. Prognostic pathways for viral infections
EP4136261A1 (en) 2020-04-16 2023-02-22 InnoSIGN B.V. Prognostic pathways for viral infections
EP3940704A1 (en) * 2020-07-14 2022-01-19 Koninklijke Philips N.V. Method for determining the differentiation state of a stem cell
EP3960875A1 (en) 2020-08-28 2022-03-02 Koninklijke Philips N.V. Pcr method and kit for determining pathway activity
EP3974540A1 (en) 2020-09-25 2022-03-30 Koninklijke Philips N.V. Method for predicting immunotherapy resistance
EP4015651A1 (en) 2020-12-17 2022-06-22 Koninklijke Philips N.V. Treatment prediction and effectiveness of anti-tnf alpha treatment in ibd patients
CN112582030B (zh) * 2020-12-18 2023-08-15 广州大学 一种基于dna存储介质的文本存储方法
EP4039825A1 (en) 2021-02-09 2022-08-10 Koninklijke Philips N.V. Comparison and standardization of cell and tissue culture
CN117295824A (zh) 2021-03-11 2023-12-26 皇家飞利浦有限公司 用于高风险脓毒症患者的预后通路
CN113177359B (zh) * 2021-04-30 2023-04-18 上海电机学院 一种基于假人模型的身体组织状态预测方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4711955A (en) 1981-04-17 1987-12-08 Yale University Modified nucleotides and methods of preparing and using same
US5660985A (en) 1990-06-11 1997-08-26 Nexstar Pharmaceuticals, Inc. High affinity nucleic acid ligands containing modified nucleotides
US5658751A (en) 1993-04-13 1997-08-19 Molecular Probes, Inc. Substituted unsymmetrical cyanine dyes with selected permeability
US5436134A (en) 1993-04-13 1995-07-25 Molecular Probes, Inc. Cyclic-substituted unsymmetrical cyanine dyes
US6783934B1 (en) 2000-05-01 2004-08-31 Cepheid, Inc. Methods for quantitative analysis of nucleic acid amplification reaction
US20110166028A1 (en) * 2007-08-28 2011-07-07 Donald Bergstrom Methods for predicting treatment response based on the expression profiles of biomarker genes in notch mediated cancers
US7544476B1 (en) * 2008-07-11 2009-06-09 Aveo Pharmaceuticals, Inc. Identifying cancers sensitive to treatment with inhibitors of notch signaling
EP2549399A1 (en) 2011-07-19 2013-01-23 Koninklijke Philips Electronics N.V. Assessment of Wnt pathway activity using probabilistic modeling of target gene expression
EP2606884A1 (en) * 2011-12-21 2013-06-26 Ecole Polytechnique Fédérale de Lausanne (EPFL) Inhibitors of notch signaling pathway and use thereof in treatment of cancers
KR102279844B1 (ko) 2012-12-26 2021-07-21 코닌클리케 필립스 엔.브이. 타깃 유전자 발현의 선형 조합(들)을 사용한 세포 시그널링 경로 활성의 평가
AU2015334840B2 (en) * 2014-10-24 2021-10-21 Innosign B.V. Assessment of TGF-beta cellular signaling pathway activity using mathematical modelling of target gene expression
EP3210144B1 (en) * 2014-10-24 2020-10-21 Koninklijke Philips N.V. Medical prognosis and prediction of treatment response using multiple cellular signaling pathway activities
US20160125127A1 (en) * 2014-10-29 2016-05-05 Council Of Scientific & Industrial Research Identification of minimal combinations of oncoproteins in notch pathway to suppress human glioblastoma
JP7028763B2 (ja) * 2015-08-14 2022-03-02 コーニンクレッカ フィリップス エヌ ヴェ 標的遺伝子発現の数学的モデリングを使用するNFkB細胞シグナリング経路活性の評価
KR20160086775A (ko) * 2016-02-11 2016-07-20 성균관대학교산학협력단 Pin1 저해제를 유효성분으로 포함하는 뇌졸중 예방 또는 치료용 약학적 조성물

Also Published As

Publication number Publication date
EP3692537A1 (en) 2020-08-12
US20190100790A1 (en) 2019-04-04
WO2019068585A1 (en) 2019-04-11
EP3462349A1 (en) 2019-04-03
CN111183233A (zh) 2020-05-19
AU2018344383A1 (en) 2020-05-21
BR112020006484A2 (pt) 2020-10-13
CA3076635A1 (en) 2019-04-11
EP3692537B1 (en) 2022-01-19

Similar Documents

Publication Publication Date Title
JP2020536530A (ja) 標的遺伝子発現の数学的モデル化を使用する、Notch細胞シグナル伝達経路活性の評価
JP6721665B2 (ja) 目的遺伝子発現の線形結合を用いた細胞信号伝達経路活性の評価
JP6415712B2 (ja) 標的遺伝子発現の数学的モデル化を用いるTGF−β細胞シグナル伝達経路活性の評価
JP7368483B2 (ja) 相同組換え欠損を推定するための統合された機械学習フレームワーク
Halloran et al. Molecular assessment of disease states in kidney transplant biopsy samples
JP6931125B2 (ja) 標的遺伝子発現の数学的モデル化を使用する、jak−stat1/2細胞シグナル伝達経路活性の評価
JP7028763B2 (ja) 標的遺伝子発現の数学的モデリングを使用するNFkB細胞シグナリング経路活性の評価
EP3729439B1 (en) Assessment of mapk-ap 1 cellular signaling pathway activity using mathematical modelling of target gene expression
KR102256141B1 (ko) 다수의 세포 신호전달 경로 활성을 이용하는 치료 반응의 의학적 예후 및 예측
JP6204354B2 (ja) 標的遺伝子発現の確率モデルを使用した細胞シグナル経路の活性の評価
JP2020535823A (ja) 標的遺伝子発現の数学的モデル化を使用する、jak−stat3細胞シグナル伝達経路活性の評価
JP2022544604A (ja) がん検体において細胞経路調節不全を検出するためのシステム及び方法
CA2953368C (en) Methods for diagnosing risk of renal allograft fibrosis and rejection
CN111164700A (zh) 使用下一代测序的免疫-肿瘤学应用
CN115701286A (zh) 使用无循环mRNA谱分析检测阿尔茨海默病风险的系统和方法
EP4058603A1 (en) Assessment of pr cellular signaling pathway activity using mathematical modelling of target gene expression
Tomoiaga et al. Single-Cell, Human Sperm Transcriptomes and Variants from Fathers of Autistic and Healthy Children
CA2965442C (en) Assessment of tgf-b cellular signaling pathway activity using mathematical modelling of target gene expression

Legal Events

Date Code Title Description
RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20200713

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210924

A761 Written withdrawal of application

Free format text: JAPANESE INTERMEDIATE CODE: A761

Effective date: 20220509