JP7065609B2 - 複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測 - Google Patents

複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測 Download PDF

Info

Publication number
JP7065609B2
JP7065609B2 JP2017522163A JP2017522163A JP7065609B2 JP 7065609 B2 JP7065609 B2 JP 7065609B2 JP 2017522163 A JP2017522163 A JP 2017522163A JP 2017522163 A JP2017522163 A JP 2017522163A JP 7065609 B2 JP7065609 B2 JP 7065609B2
Authority
JP
Japan
Prior art keywords
pathway
risk
activity
cell signaling
target genes
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
Application number
JP2017522163A
Other languages
English (en)
Other versions
JP2018502553A (ja
JP7065609B6 (ja
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 JP2018502553A publication Critical patent/JP2018502553A/ja
Application granted granted Critical
Publication of JP7065609B2 publication Critical patent/JP7065609B2/ja
Publication of JP7065609B6 publication Critical patent/JP7065609B6/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • 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
    • G16B5/20Probabilistic models
    • 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
    • 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
    • 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
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • 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
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/20Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
    • 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
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/30Detection of binding sites or motifs
    • 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
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/112Disease subtyping, staging or classification
    • 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/118Prognosis of disease development
    • 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
    • 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/16Primer sets for multiplex assays

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Genetics & Genomics (AREA)
  • Biophysics (AREA)
  • Biotechnology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Evolutionary Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Pathology (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Public Health (AREA)
  • Immunology (AREA)
  • Biomedical Technology (AREA)
  • Microbiology (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Physiology (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)
  • Probability & Statistics with Applications (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Description

本明細書に記載の主題は、主にバイオ・インフォマティクス、ゲノム処理技術分野、プロテオミクス処理技術分野、及び関連技術分野に関する。より詳細には、本発明は、対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための、デジタル処理デバイスによって実行される、コンピュータ実装方法に関連する。このリスク・スコアは、対象における2つ以上の細胞シグナル伝達経路についての推測活性の組合せに基づいて決定される。本発明はさらに、対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための、上記の方法を実行するように構成された、デジタル・プロセッサを備える装置と、対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための、上記の方法を実行するためにデジタル処理デバイスによって実行可能である命令を記憶する非一時的記憶媒体と、対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための、デジタル処理デバイスで実行されたときにデジタル処理デバイスに方法を実行させるプログラム・コード手段を含むコンピュータ・プログラムとに関する。本発明はさらに、対象の試料中の2つ以上の細胞シグナル伝達経路それぞれの3つ以上の標的遺伝子の発現レベルを測定するためのキットと、患者が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのキットと、上記の方法を実行する際の、このキットの使用とに関する。
ゲノム及びプロテオミクスの分析は、腫瘍学などの医療分野における臨床応用に対し、実質的に実現されており、また潜在的な有望性がある。様々な癌が、ゲノムの突然変異/バリエーション/異常メチル化のパターンの特定の組合せ、及び/又は特定の遺伝子の高発現レベル若しくは低発現レベルと関連することが知られており、これらは、癌の成長及び進化、たとえば細胞増殖及び転移においての役割を果たす。たとえば、Wntシグナル伝達経路は、細胞増殖の調節に影響を及ぼし、高度に調節されている。調節が損なわれることによる高いWnt経路活性は、癌、中でも悪性大腸腫瘍と相関付けられてきた。どのような特定の動作原理にも限定されないが、悪性大腸細胞のWnt経路の調節解除が高いWnt経路活性につながり、ひいては、悪性大腸細胞の細胞増殖、すなわち大腸癌の転移をもたらすと考えられている。一方で、異常に低い経路活性もまた、たとえば骨粗しょう症の場合において注目される。健康及び疾患に関して、細胞の分裂、機能及び/又は分化における同様の役割を果たす他の経路としては、(たとえば、ER、PR、AR、PPAR、GR、VitD、TGF-β、ノッチ、ヘッジホッグ、FGF、NFkB、VEGF、及びPDGF等の)細胞シグナル伝達経路がある。
ゲノム及びプロテオミクスのデータを取得する技術が、臨床環境において容易に利用できるようになっている。たとえば、マイクロアレイによる測定が、遺伝子発現レベル、タンパク質レベル、メチル化などを評価するのに日常的に使用されている。自動化された遺伝子配列決定により、DNA及びmRNAにおける遺伝子のバリエーション/突然変異/異常メチル化パターン識別が、高い費用対効果で可能になっている。遺伝子配列決定時の定量的なmRNAレベルの評価が、遺伝子発現レベルを評価するための臨床ツールとして有望である。
治療専門家、たとえば腫瘍専門家にとっての主な難題の1つは、患者の予後について経験や知識に基づく推測をすることである。その理由は、この情報が治療の選択に影響を及ぼすからである。個々の患者の癌組織試料に基づくゲノムミクス、トランスクリプトミクス及びプロテオミクス(並びに他の「オミクス」)解析により、患者の予後評価に潜在的に寄与できる情報が得られる。
しかし、関連する臨床情報を引き出すためにこれらの複雑なデータを解釈するのは難題であることが分かっており、大部分が依然として未解決である。患者の予後は、たとえば、「(疾患の)再発までの時間」、「(疾患の)進行までの時間」、「(疾患の)発生の時点」、又は「死亡(疾患)までの時間」のように、いくつかの方法で定量的に示すことができる。
本発明の主態様によれば、上記の問題は、対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのコンピュータ実装方法によって解決されるか、又は少なくとも軽減され、その決定するステップは、
対象における2つ以上の細胞シグナル伝達経路のそれぞれの活性を、対象の試料において測定されたそれぞれの細胞シグナル伝達経路の3個以上の、たとえば、3個、4個、5個、6個、7個、8個、9個、10個、11個、又は12個以上の、標的遺伝子の発現レベルに基づいて推測するステップと、
リスク・スコアを、推測された活性を組み合わせたものに基づいて決定するステップと、
を含み、
臨床事象は、疾患再発、疾患進行、疾患発症、及び疾患が原因の死亡のうちの1つであり、疾患は癌、好ましくは乳癌であり、
細胞シグナル伝達経路は、ホスファチジルイノシチド3-キナーゼ(PI3K)経路と、Wnt経路、エストロゲン受容体(ER)経路、及びヘッジホッグ(HH)経路のうちの1つ以上の経路とを含む。
本発明は、好ましくは、疾患と関連する臨床事象を規定期間内に、たとえば3カ月以内、6カ月以内、1年以内、18カ月以内、2年以内、30カ月以内、3年以内、42カ月以内、4年以内、5年以内、6年以内、7年以内、8年以内、9年以内、又は10年以上内に経験するリスクがある対象の識別を可能にする。
本明細書で「対象」という用語は、任意の生物を指す。いくつかの実施形態では、対象は動物、好ましくは哺乳動物である。いくつかの実施形態では、対象はヒトであり、好ましくは医療対象者である。
本明細書で「転写因子要素」という用語は、好ましくは、活性転写因子の中間体若しくは前駆体のタンパク質若しくはタンパク質複合体、又は指定の標的遺伝子発現を制御する活性転写因子タンパク質若しくはタンパク質複合体を指す。
本明細書で「標的遺伝子」という用語は、転写がそれぞれの転写因子要素によって直接又は間接に制御される遺伝子を意味する。「標的遺伝子」は、「直接標的遺伝子」及び/又は「間接標的遺伝子」(本明細書に記載)であり得る。
対象における細胞シグナル伝達経路の活性を推測することは、たとえば、特に(i)対象の試料において測定された細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを含む入力の組に対する細胞シグナル伝達経路を表す較正済み確率的経路モデルの、好ましくはベイジアン・ネットワークの、一部分を評価することによって、また(ii)対象における転写因子(TF)要素の活性レベルを推定することによって行うことができ、このTF要素は、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御し、推定することは、TF要素の活性レベルと、対象の試料において測定された細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルとを関連付ける条件付き確率に基づいており、また(iii)細胞シグナル伝達経路の活性を対象の試料におけるTF要素の推定活性レベルに基づいて推測することによって行うことができる。これについては、公開されている国際特許出願WO2013/011479A2(「Assessment of cellular signaling pathway activity using probabilistic modeling of target gene expression」)に詳細に記載されており、同出願の内容はその全体が本明細書に組み込まれる。
例示的な一代替形態では、対象における細胞シグナル伝達経路の活性を推測することは、たとえば、特に(i)対象の試料における転写因子(TF)要素の活性レベルを決定することによって行うことができ、このTF要素は、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御し、決定することは、細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルをTF要素の活性レベルに関係付ける較正済み数学経路モデルを評価することに基づいており、数学経路モデルは、3つ以上の標的遺伝子の発現レベルの1つ以上の一次結合に基づいており、また(ii)対象における細胞シグナル伝達経路の活性を対象の試料におけるTF要素の決定された活性レベルに基づいて推定することによって行うことができる。これについては、公開されている国際特許出願WO2014/102668A2(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)に詳細に記載されている。
好ましい実施形態によれば、細胞シグナル伝達経路は、Wnt経路、ER経路、及びHH経路を含む。
PI3K経路、Wnt経路、ER経路、及びHH経路のそれぞれは、好ましくは、その経路と関連する転写因子(TF)複合体の転写活性に最終的につながる細胞シグナル伝達経路として定義される。好ましくは、これらは、それぞれ少なくともFOXOファミリー・メンバー、β-カテニン/TCF4、ERα二量体、及びGLIファミリー・メンバーから成る。
本発明は、PI3K経路と、活性がPI3K経路の活性とは実質的に負の相関関係にあるFOXO TFファミリーとに狙いを定めている。すなわちFOXOの活性は、PI3K経路の不活性と実質的に相関関係があるのに対し、FOXOの不活性は、PI3K経路の活性と実質的に相関関係がある。
細胞シグナル伝達経路は、Wnt経路及び/又はER経路及び/又はHH経路を含むことが好ましく、リスク・スコアは、示されたリスクがWnt経路の推測活性の増加、及び/又はHH経路の推測活性の増加と共に単調増加し、かつ/又はER経路の推測活性の増加と共に単調減少するように規定される。
リスク・スコアは、示されたリスクがPI3K経路の推測活性の増加と共に単調増加するように規定されることもまた好ましい。
MPSは、項w・Pと、項w・P、w・Pe、及びw・Pのうちの1つ以上とを含む合計から成り、ここで、P、P、P、及びPは、それぞれPI3K経路、Wnt経路、ER経路、及びHH経路の推測活性を表し、w、w、w、及びwは、対象が臨床事象を既定期間内に経験するリスクと、PI3K経路、Wnt経路、ER経路、及びHH経路それぞれの活性との間の相関関係を表す一定重み付け係数である。
一定重み付け係数w、w、w、及びwは、それぞれの細胞シグナル伝達経路のコックス比例ハザード・モデルを臨床データに適合させることにより得られるコックスの係数の値に基づいてそれぞれ決定されるか、又は決定されている。
本発明は、本明細書に記載の細胞シグナル伝達経路に生じる効果を識別する適切な方法が、細胞シグナル伝達経路のシグナル伝達出力の測定に基づくことができるという本発明者らの新考案に基づいており、このシグナル伝達出力は、とりわけ、シグナル伝達経路が制御する転写因子(TF)要素によって制御される、本明細書に記載の固有の標的遺伝子の転写である。本発明者らによるこの新考案では、TF活性レベルが、試料において、とりわけ一意的に識別された標的遺伝子の発現値によって検出できる準安定状態にあると仮定している。
特に、発現レベルが数学経路モデルで解析される細胞シグナル伝達経路標的遺伝子の固有の組が特定された。数学経路モデルで使用するために、評価される各細胞シグナル伝達経路から3個以上の、たとえば、3個、4個、5個、6個、7個、8個、9個、10個、11個、又は12個以上の標的遺伝子を解析してリスク・スコアを策定することができる。
好ましいのは、
3つ以上のPI3K標的遺伝子が、ATP8A1、BCL2L11、BNIP3、BTG1、C10orf10、CAT、CBLB、CCND1、CCND2、CDKN1B、DDB1、DYRK2、ERBB3、EREG、ESR1、EXT1、FASLG、FGFR2、GADD45A、IGF1R、IGFBP1、IGFBP3、INSR、LGMN、MXI1、PPM1D、SEMA3C、SEPP1、SESN1、SLC5A3、SMAD4、SOD2、TLE4、及びTNFSF10から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、ADRA2C、ASCL2、AXIN2、BMP7、CCND1、CD44、COL18A1、DEFA6、DKK1、EPHB2、EPHB3、FAT1、FZD7、GLUL、HNF1A、CXCL8(以前にはIL8として知られていた)、CEMIP(以前にはKIAA1199として知られていた)、KLF6、LECT2、LEF1、LGR5、MYC、NKD1、OAT、PPARG、REGIB、RNF43、SLC1A2、SOX9、SP5、TBX3、TCF7L2、TDGF1、及びZNRF3から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、AP1B1、ATP5J、COL18A1、COX7A2L、CTSD、DSCAM、EBAG9、ESR1、HSPB1、KRT19、NDUFV3、NRIP1、PGR、PISD、PRDM15、PTMA、RARA、SOD1、TFF1、TRIM25、XBP1、GREB1、IGFBP4、MYC、SGK3、WISP2、ERBB2、CA12、CDH26、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、HHIP、SPP1、TSC22D1、CCND2、H19、IGFBP6、TOM1、JUP、FOXA2、MYCN、NKX2-2、NKX2-8、RAB34、MIF、GLI3、FST、BCL2、CTSL1、TCEA2、MYLK、FYN、PITRM1、CFLAR、IL1R2、S100A7、S100A9、CCND1、JAG2、FOXM1、FOXF1、及びFOXL1から成る群から選択されることである。
さらに好ましいのは、
3つ以上のPI3K標的遺伝子が、AGRP、BCL2L11、BCL6、BNIP3、BTG1、CAT、CAV1、CCND1、CCND2、CCNG2、CDKN1A、CDKN1B、ESR1、FASLG、FBXO32、GADD45A、INSR、MXI1、NOS3、PCK1、POMC、PPARGC1A、PRDX3、RBL2、SOD2、TNFSF10から成る群から選択され、好ましくはFBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、CEMIP、AXIN2、CD44、RNF43、MYC、TBX3、TDGF1、SOX9、ASCL2、CXCL8、SP5、ZNRF3、EPHB2、LGR5、EPHB3、KLF6、CCND1、DEFA6、及びFZD7から成る群から選択され、好ましくはAXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、CDH26、SGK3、PGR、GREB1、CA12、XBP1、CELSR2、WISP2、DSCAM、ERBB2、CTSD、TFF1、PDZK1、IGFBP4、ESR1、SOD1、AP1B1、及びNRIP1から成る群から選択され、好ましくはTFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、IGFBP6、SPP1、CCND2、FST、FOXL1、CFLAR、TSC22D1、RAB34、S100A9、S100A7、MYCN、FOXM1、GLI3、TCEA2、FYN、及びCTSL1から成る群から選択され、好ましくはGLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9から成る群から選択されることである。
特に好ましいのは、
3つ以上のPI3K標的遺伝子が、FBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1であり、かつ/又は
3つ以上のWnt標的遺伝子が、AXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6であり、かつ/又は
3つ以上のER標的遺伝子が、TFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2であり、かつ/又は、
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9であることである。
本発明の別の態様、この方法はさらに、
対象を、対象が臨床事象を規定期間内に経験することの別々に示されたリスクと関連付けられた複数のリスク群のうちの少なくとも1つに割り当てるステップ、
及び/又は
対象に推奨される治療を、対象が臨床事象を規定期間内に経験することの示されたリスクに基づいて決定するステップをさらに含む。
本発明により使用されるべき試料は、抽出試料、すなわち対象から取り出された試料とすることができる。試料の例には、それだけには限らないが、対象の組織、細胞、血液及び/又は体液が含まれる。これは、たとえば、癌病変部から、又は癌の疑いがある病変部から、又は転移腫瘍から、又は癌細胞で汚染されている流体が存在する体腔(たとえば、胸膜若しくは腹部の腔、又は膀胱腔)から、又は癌細胞を含有する他の体液などから、好ましくは組織診手順又は他の試料抽出手順を介して得られた試料とすることができる。試料が取り出される細胞はまた、悪性血液疾患(白血病又はリンパ腫など)からの腫瘍細胞とすることができる。場合によっては、細胞試料はまた、循環腫瘍細胞、すなわち血流に入った腫瘍細胞であることもあり、適切な分離技法、たとえばアフェレーシス又は従来の静脈血採血を用いて取り出すことができる。血液は別として、試料が取り出される体液は、尿、胃腸内容物、又は血管外湧出物であり得る。本明細書で「試料」という用語はまた、たとえば、対象の組織及び/又は細胞及び/又は体液が対象から採取され、たとえば顕微鏡スライド若しくは固定剤に載せられている場合と、特許請求された方法を実施するために、この試料の一部分が、たとえば、レーザーキャプチャ法(LCM)によって、又は穿孔によって、又は対象の細胞をスライドからかき取ることによって、又は蛍光活性化細胞選別技法によって取り出される場合とを包含する。加えて、本明細書で「試料」という用語はまた、たとえば、対象の組織及び/又は細胞及び/又は体液が対象から採取され、顕微鏡スライドに載せられ、また特許請求された方法がスライド上で実施される場合も含む。
この方法はさらに、リスク・スコア及び/又は推測活性のうちの少なくとも1つを、1つ又は複数の追加予後検査から得られた1つ又は複数の追加リスク・スコアと組み合わせて結合リスク・スコアを得ることを含むのが好ましく、この結合リスク・スコアは、対象が臨床事象を規定期間内に経験するリスクを示す。1つ又は複数の追加予後検査には特に、Oncotype DX(登録商標)乳癌検査、Mammostrat(登録商標)乳癌検査、MammaPrint(登録商標)乳癌検査、EndoPredict(登録商標)乳癌検査、BluePrint(商標)乳癌検査、CompanDx(登録商標)乳癌検査、Breast Cancer Index(役務標章)(HOXB13/IL17BR)、OncotypeDX(登録商標)結腸癌検査、及び/又は、遺伝子/タンパク質Ki67の発現を測定することによって行われる増殖検査が含まれ得る。
上記のように、臨床事象は、疾患再発、疾患進行、疾患発症、及び疾患が原因の死亡のうちの1つであり、疾患が癌であり、好ましくは乳癌である。この場合、臨床事象が既定期間内に生じるリスクは優先的に、所与の治療後の(「癌治療応答予測」とも呼ばれる)、又は全く治療なしでの(「癌予後」とも呼ばれる)癌の再発の、すなわちぶり返しの、リスクになる。再発は、局部的であること(すなわち、元の腫瘍のそば)も、離れていること(すなわち、元の側を超えた転移)もある。別の代替形態では、臨床事象が既定期間内に生じるリスクは、癌の進行のリスク、癌の発生のリスク、又は癌が原因の死亡のリスクである。
別の開示態様によれば、対象が疾患と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための装置は、本明細書に記載の本明細書の方法を実施するように構成されたデジタル・プロセッサを備える。
別の開示態様によれば、対象が疾患と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための非一時的記憶媒体は、本明細書に記載の本発明の方法を実施するデジタル処理デバイスによって実行可能である命令を記憶する。非一時的記憶媒体は、ハード・ドライブ若しくは他の磁気記憶媒体、光ディスク若しくは他の光記憶媒体、ランダム・アクセス・メモリ(RAM)、読み出し専用メモリ(ROM)、フラッシュ・メモリ、又は他の電子記憶媒体、ネットワーク・サーバなどの、コンピュータ可読記憶媒体とすることができる。デジタル処理デバイスは、手持ち型デバイス(たとえば、携帯情報端末又はスマートフォン)、ノート型コンピュータ、デスクトップ・コンピュータ、タブレット・コンピュータ若しくはデバイス、リモート・ネットワーク・サーバなどであり得る。
別の開示態様によれば、対象が疾患と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのコンピュータ・プログラムは、コンピュータ・プログラムがデジタル処理デバイス上で実行されたとき、デジタル処理デバイスに本明細書に記載の本発明の方法を実施させるプログラム・コード手段を含む。デジタル処理デバイスは、手持ち型デバイス(たとえば、携帯情報端末又はスマートフォン)、ノート型コンピュータ、デスクトップ・コンピュータ、タブレット・コンピュータ若しくはデバイス、リモート・ネットワーク・サーバなどであり得る。
別の開示態様によれば、対象の試料における2つ以上の細胞シグナル伝達経路それぞれの、3個以上の、たとえば、3個、4個、5個、6個、7個、8個、9個、10個、11個、又は12個以上の、標的遺伝子の発現レベルを測定するためのキットは、
対象の試料における、それぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを決定するための1つ又は複数の構成要素を含み、
細胞シグナル伝達経路は、TGF-β経路と、PI3K経路、Wnt経路、ER経路、及びHH経路のうちの1つ以上とを含む。
それぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを測定するための1つ又は複数の構成要素又は手段は、DNAアレイ・チップ、オリゴヌクレオチド・アレイ・チップ、タンパク質アレイ・チップ、抗体、複数のプローブ、たとえば標識されたプローブ、RNAリバーサ転写酵素配列決定構成要素の組、及び/又は、RNA又はcDNAを含むDNAの増幅プライマーから成る群から選択することができる。一実施形態では、キットは、本明細書に記載の3つ以上の標的遺伝子のmRNA又はcDNA配列の一部分を対象とする標識されたプローブの組を含む。一実施形態では、キットは、さらに以下に記載された3つ以上の標的遺伝子のmRNA又はcDNA配列の一部分を対象とするプライマーとプローブの組、たとえば、表16から表19の配列から選択された特定のプライマーとプローブの組を含む。一実施形態では、標識されたプローブは、標準化96ウェル・プレートに含まれる。一実施形態では、キットはさらに、たとえば表20に表された参照遺伝子の組を対象とするプライマー又はプローブを含む。このような参照遺伝子は、たとえば、本明細書に記載の標的遺伝子発現レベルの発現レベルを正規化又は標準化するのに有用な、構造的に発現された遺伝子とすることができる。
一実施形態では、対象の試料における2つ以上の細胞シグナル伝達経路それぞれの、3つ以上の標的遺伝子の発現レベルを測定するためのキットは、
それぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子を対象とするポリメラーゼ連鎖反応プライマーと、
それぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子を対象とするプローブとを含み、
細胞シグナル伝達経路がPI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上とを含む。
好ましいのは、
3つ以上のPI3K標的遺伝子が、ATP8A1、BCL2L11、BNIP3、BTG1、C10orf10、CAT、CBLB、CCND1、CCND2、CDKN1B、DDB1、DYRK2、ERBB3、EREG、ESR1、EXT1、FASLG、FGFR2、GADD45A、IGF1R、IGFBP1、IGFBP3、INSR、LGMN、MXI1、PPM1D、SEMA3C、SEPP1、SESN1、SLC5A3、SMAD4、SOD2、TLE4、及びTNFSF10から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、ADRA2C、ASCL2、AXIN2、BMP7、CCND1、CD44、COL18A1、DEFA6、DKK1、EPHB2、EPHB3、FAT1、FZD7、GLUL、HNF1A、CXCL8、CEMIP、KLF6、LECT2、LEF1、LGR5、MYC、NKD1、OAT、PPARG、REG1B、RNF43、SLC1A2、SOX9、SP5、TBX3、TCF7L2、TDGF1、及びZNRF3から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、AP1B1、ATP5J、COL18A1、COX7A2L、CTSD、DSCAM、EBAG9、ESR1、HSPB1、KRT19、NDUFV3、NRIP1、PGR、PISD、PRDM15、PTMA、RARA、SOD1、TFF1、TRIM25、XBP1、GREB1、IGFBP4、MYC、SGK3、WISP2、ERBB2、CA12、CDH26、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、HHIP、SPP1、TSC22D1、CCND2、H19、IGFBP6、TOM1、JUP、FOXA2、MYCN、NKX2-2、NKX2-8、RAB34、MIF、GLI3、FST、BCL2、CTSL1、TCEA2、MYLK、FYN、PITRM1、CFLAR、IL1R2、S100A7、S100A9、CCND1、JAG2、FOXM1、FOXF1、及びFOXL1から成る群から選択されることである。
さらに好ましいのは、
3つ以上のPI3K標的遺伝子が、AGRP、BCL2L11、BCL6、BNIP3、BTG1、CAT、CAV1、CCND1、CCND2、CCNG2、CDKN1A、CDKN1B、ESR1、FASLG、FBXO32、GADD45A、INSR、MXI1、NOS3、PCK1、POMC、PPARGC1A、PRDX3、RBL2、SOD2、TNFSF10から成る群から選択され、好ましくはFBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、CEMIP、AXIN2、CD44、RNF43、MYC、TBX3、TDGF1、SOX9、ASCL2、CXCL8、SP5、ZNRF3、EPHB2、LGR5、EPHB3、KLF6、CCND1、DEFA6、及びFZD7から成る群から選択され、好ましくはAXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、CDH26、SGK3、PGR、GREB1、CA12、XBP1、CELSR2、WISP2、DSCAM、ERBB2、CTSD、TFF1、PDZK1、IGFBP4、ESR1、SOD1、AP1B1、及びNRIP1から成る群から選択され、好ましくはTFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、IGFBP6、SPP1、CCND2、FST、FOXL1、CFLAR、TSC22D1、RAB34、S100A9、S100A7、MYCN、FOXM1、GLI3、TCEA2、FYN、及びCTSL1、から成る群から選択され、好ましくはGLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9から成る群から選択されることである。
特に好ましいのは、
3つ以上のPI3K標的遺伝子が、FBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1であり、かつ/又は
3つ以上のWnt標的遺伝子が、AXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6であり、かつ/又は
3つ以上のER標的遺伝子が、TFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2であり、かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9であることである。
別の開示態様によれば、対象が疾患と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのキットは、
本明細書に記載の本発明のキットと、
本明細書に記載の本発明の装置、本明細書に記載の本発明の非一時的記憶媒体、又は本明細書に記載の本発明のコンピュータ・プログラムとを含む。
別の開示態様によれば、対象の試料において2つ以上の細胞シグナル伝達経路それぞれの3個以上の標的遺伝子の発現レベルを測定するための、及び、任意で、対象が疾患と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのキットは、
対象の試料において2つ以上の細胞シグナル伝達経路それぞれの、3個以上の、たとえば3個、4個、5個、6個、7個、8個、9個、10個、11個、又は12個以上の、標的遺伝子の発現レベルを決定するための1つ又は複数の構成要素を含み、
1つ又は複数の構成要素は、好ましくは、(例えばDNAアレイ・チップ、オリゴヌクレオチド・アレイ・チップ、タンパク質アレイ・チップ等の)マイクロアレイ・チップ、抗体、複数のプローブ、たとえば標識されたプローブ、RNAリバーサ転写酵素配列決定構成要素の組、及び/又は、RNA又はcDNAを含むDNAの増幅プライマーから成る群から選択され、
細胞シグナル伝達経路は、PI3Kシグナル伝達経路と、Wntシグナル伝達経路、ERシグナル伝達経路、及びHHシグナル伝達経路のうちの1つ以上の経路とを含み、
任意で、当該キットは、本明細書に記載の本発明の装置、本明細書に記載の本発明の非一時的記憶媒体、又は本明細書に記載の本発明のコンピュータ・プログラムを含む。
好ましいのは、
3つ以上のPI3K標的遺伝子が、ATP8A1、BCL2L11、BNIP3、BTG1、C1Oorf1O、CAT、CB1B、CCND1、CCND2、CDKN1B、DDB1、DYRK2、ERBB3、EREG、ESR1、EXT1、FASLG、FGFR2、GADD45A、IGF1R、IGFBP1、IGFBP3、INSR、LGMN、MXI1、PPM1D、SEMA3C、SEPP1、SESN1、SLC5A3、SMAD4、SOD2、TLE4、及びTNFSF10から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、ADRA2C、ASCL2、AXIN2、BMP7、CCND1、CD44、COL18A1、DEFA6、DKK1、EPHB2、EPHB3、FAT1、FZD7、GLUL、HNF1A、CXCL8、CEMIP、KLF6、LECT2、LEF1、LGR5、MYC、NKD1、OAT、PPARG、REG1B、RNF43、SLC1A2、SOX9、SP5、TBX3、TCF7L2、TDGF1、及びZNRF3から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、AP1B1、ATP5J、COL18A1、COX7A2L、CTSD、DSCAM、EBAG9、ESR1、HSPB1、KRT19、NDUFV3、NRIP1、PGR、PISD、PRDM15、PTMA、RARA、SOD1、TFF1、TRIM25、XBP1、GREB1、IGFBP4、MYC、SGK3、WISP2、ERBB2、CA12、CDH26、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、HHIP、SPP1、TSC22D1、CCND2、H19、IGFBP6、TOM1、JUP、FOXA2、MYCN、NKX2-2、NKX2-8、RAB34、MIF、GLI3、FST、BCL2、CTSL1、TCEA2、MYLK、FYN、PITRM1、CFLAR、IL1R2、S100A7、S100A9、CCND1、JAG2、FOXM1、FOXF1、及びFOXL1から成る群から選択されることである。
さらに好ましいのは、
3つ以上のPI3K標的遺伝子が、AGRP、BCL2L11、BCL6、BNIP3、BTG1、CAT、CAV1、CCND1、CCND2、CCNG2、CDKN1A、CDKN1B、ESR1、FASLG、FBXO32、GADD45A、INSR、MXI1、NOS3、PCK1、POMC、PPARGC1A、PRDX3、RBL2、SOD2、TNFSF10、から成る群から選択され、好ましくは、FBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1から成る群から選択され、
かつ/又は
3つ以上のWnt標的遺伝子が、CEMIP、AXIN2、CD44、RNF43、MYC、TBX3、TDGF1、SOX9、ASCL2、CXCL8、SP5、ZNRF3、EPHB2、LGR5、EPHB3、KLF6、CCND1、DEFA6、及びFZD7から成る群から選択され、好ましくは、AXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6から成る群から選択され、
かつ/又は
3つ以上のER標的遺伝子が、CDH26、SGK3、PGR、GREB1、CA12、XBP1、CELSR2、WISP2、DSCAM、ERBB2、CTSD、TFF1、PDZK1、IGFBP4、ESR1、SOD1、AP1B1、及びNRIP1、から成る群から選択され、好ましくは、TFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2から成る群から選択され、
かつ/又は
3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、IGFBP6、SPP1、CCND2、FST、FOXL1、CFLAR、TSC22D1、RAB34、S100A9、S100A7、MYCN、FOXM1、GLI3、TCEA2、FYN、及びCTSL1、から成る群から選択され、好ましくは、GLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9から成る群から選択されることである。
特に好ましいのは、
3つ以上のPI3K標的遺伝子は、FBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1であり、かつ/又は
3つ以上のWnt標的遺伝子は、AXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6であり、かつ/又は
3つ以上のER標的遺伝子は、TFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2あり、かつ/又は
3つ以上のHH標的遺伝子は、GLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9であることである。
別の開示態様によれば、本明細書に記載の本発明のキットは、本明細書に記載の本発明の方法を実施する際に使用される。
1つの利点は、たとえば、PI3K経路、Wnt経路、ER経路、及びHH経路の確率論的又は他の数学経路モデルを使用することなどによる2つ以上の細胞シグナル伝達経路の解析に基づいて、特に、対象が疾患、たとえば癌、特に乳癌と関連した臨床事象を既定期間内に経験するリスクであって、細胞シグナル伝達経路の推測活性を組み合わせたものに基づいて決定されるリスク・スコアによって示されたリスクに基づいて、対象に対する治療を決定することによって臨床推奨を提供するように適合されている臨床判断支援(CDS)システムにある。
別の利点は、対象が疾患、たとえば癌、特に乳癌と関連した臨床事象を既定期間内に経験する別々のリスクであって、2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づいて決定されるリスク・スコアによって示された別々のリスクと関連付けられた、複数のリスク群のうちの少なくとも1つに対象を割り当てるように適合されているCDSシステムにある。
別の利点は、対象が疾患、たとえば癌、特に乳癌と関連した臨床事象を既定期間内に経験する、2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づいて決定されるリスク・スコアを、1つ又は複数の追加予後検査から得られた1つ又は複数のリスク・スコアと組み合わせることにある。
本明細書に記載の本発明はまた、たとえば、2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づく予後及び/又は予測に、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づく、たとえば化学療法及び/又はホルモン治療の薬効の予測に、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づく薬効の監視に、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づいて、監視の頻度、又はより具体的には、療法応答監視の頻度を決定することに、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づく薬物開発に、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づくアッセイ開発に、及び/又は
2つ以上の細胞シグナル伝達経路の推測活性を組み合わせたものに基づく癌進行度分類に、関連して有利に使用することもでき、
それぞれの場合で、細胞シグナル伝達経路には、PI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上とが含まれる。
さらなる利点は、添付の図、及び以下の説明を読み理解すれば、特に本明細書で以下に提示される詳細な実施例を読めば、当業者には明らかになろう。
請求項1に記載の方法、請求項に記載の装置、請求項10に記載の非一時的記憶媒体、請求項11に記載のコンピュータ・プログラム、請求項12から14に記載のキット、及び請求項15に記載のキットの使用には、同様の、及び/又は同一の好ましい、特に、従属請求項に定義された実施形態があることを理解されたい。
本発明の好ましい実施形態はまた、従属請求項又は上記の実施形態を、それぞれの独立請求項と任意に組み合わせたものにもできることを理解されたい。
本発明の上記その他の態様は、以下に記載の実施形態から明らかになり、またそれを参照して解明されよう。
本発明では、標的遺伝子の固有の組の発現レベルの解析を利用する。特に適切な標的遺伝子は、以下の文字経路、並びに以下の実施例(たとえば、以下の表1から表13参照)に記載されている。
すなわち、一実施形態では標的遺伝子が、表1、表2、表3、表4、表5、表6、表7、表8、表9、表10、表11、表12、又は表13に列挙された標的遺伝子から成る群から選択される。
PI3K経路、Wnt経路、ER経路、及びHH経路それぞれの転写プログラムをモデル化するために使用される、本明細書ではベイジアン・ネットワーク・モデルである、数学モデルを概略的及び例示的に示す図である。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路とWnt経路の各推測活性を合わせたものである、MPSpwリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=6.1e-4)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路及びER経路の各推測活性を合わせたものである、MPSpeリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=9.2e-8)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路及びHH経路の各推測活性を合わせたものである、MPSphリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=1.4e-8)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路、Wnt経路及びER経路の各推測活性を合わせたものである、MPSpweリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=1.7e-9)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路、Wnt経路及びHH経路の各推測活性を合わせたものである、MPSpwhリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=5.4e-6)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路、ER経路及びHH経路の各推測活性を合わせたものである、MPSpehリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=2.4e-9)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路、Wnt経路、ER経路及びHH経路の各推測活性を合わせたものである、MPSpwehリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=1.8e-9)。 E-MTAB-365、GSE20685及びGSE21653の乳癌患者の無病生存期間のカプラン・マイヤー・グラフである。3つの患者群は、PI3K経路、Wnt経路、ER経路及びHH経路の選択標的遺伝子と関連付けられた各プローブセットを合わせたものである、MPSprobesetsリスク・スコアの三分位値に基づいて分けられている。高リスク患者と低リスク患者の生存曲線間の差は明らかに有意である(対数順位検定は、p=3.0e-6)。 5年(下方、実線)及び10年(上方、点線)の無病生存期間の可能性を、例として非密封MPSpwehを使用して示すグラフである。 対象が臨床事象を特定の期間内に経験するリスクを示すリスク・スコアを決定するように構成された、本明細書に開示の臨床判断支援(CDS)システムを示す図である。 リスク・スコアをPI3K経路及び追加の細胞シグナル伝達経路の標的遺伝子の発現レベルの測定に基づいて決定するプロセスを例示的に示す流れ図である。 生存データを用いた複数経路スコア(MPS)モデルを較正するプロセスを例示的に示す流れ図である。 較正された複数経路スコア(MPS)モデルからリスク・スコアを計算するプロセスを例示的に示す流れ図である。 細胞シグナル伝達経路の標的遺伝子のRT-qPCR解析からCq値を決定するプロセスを例示的に示す流れ図である。
以下の例は単に、特に好ましい方法、及びそれに関連する選択された態様を例示するにすぎない。以下で提供される教示を用いて、いくつかの試験及び/又はキットを構築することができる。以下の例は、本発明の範囲を限定するものと解釈されるべきではない。
実施例1 2つ以上の細胞シグナル伝達経路の活性の推測
公開されている国際特許出願WO2013/011479A2(「Assessment of cellular signaling pathway activity using probabilistic modeling of target gene expression」)に詳細に記載されているように、確率モデル、たとえばベイジアン・モデルを構築すると共に、いくつかの異なる標的遺伝子の発現レベルと細胞シグナル伝達経路の活性との間の条件付き確率関係を組み込むことによって、このようなモデルを使用して、細胞シグナル伝達経路の活性を高い精度で決定することができる。さらに、確率モデルは、容易に更新して、より最近の臨床研究によって得られた付加的な知識を組み込むことが、条件付き確率を調整することによって、及び/又は新しいノードをモデルに追加して付加的な情報源を表すことによって可能である。このようにして、確率モデルは、最近の医学的知識を具現化するのに適切なように更新することができる。
この手法を使用する場合、Wnt標的遺伝子、ER標的遺伝子、及びHH標的遺伝子は、好ましくはWO2013/011479A2のセクション「Example 3: Selection of target genes」及び「Example 4: Comparison of evidence curated list and broad literature list」に記載された方法に従って選択され、確率モデルは、好ましくはWO2013/011479A2の「Example 5: Training and using the Bayesian network」に記載された方法に従って訓練される。Wnt経路、ER経路、及びAR経路の活性を決定するために使用される標的遺伝子の適切な選択は、添付の特許請求の範囲に定義されている。
公開されている国際特許出願WO2014/102668A2(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)に詳細に記載されている、理解及び解釈するのが容易な別の手法では、特定の細胞シグナル伝達経路の活性は、細胞シグナル伝達経路の1つ又は複数の標的遺伝子の発現レベルと転写因子(TF)要素の活性レベルとの間の関係を組み込む数学モデル(たとえば、線形又は(擬似)線形モデル)を構築することによって決定される。このTF要素は、細胞シグナル伝達経路の1つ又は複数の標的遺伝子の転写を制御し、このモデルは、1つ又は複数の標的遺伝子の発現レベルの1つ又は複数の一次結合に基づいている。
この手法を使用する場合、Wnt標的遺伝子、ER標的遺伝子、及びHH標的遺伝子は、好ましくはWO2014/102668A2のセクション「Example 2: Selection of target genes」及び「Example 3: Comparison of evidence curated list and broad literature list」に記載された方法に従って選択され、確率モデルは、好ましくはWO2014/102668A2の「Example 4: Training and using the mathematical model」に記載された方法に従って訓練される。添付の特許請求の範囲に定義された標的遺伝子の選択もまた、Wnt経路、ER経路、及びHH経路の活性をこの後者の手法によって決定するのに有効である。
2つの異なる手法に関して、1つ又は複数の標的遺伝子の発現レベルは、好ましくはmRNAのレベルの測定値とすることができ、この測定値は、たとえば、標的遺伝子mRNA配列と関連付けられたプローブを使用する(RT)-PCR及びマイクロアレイ技法、並びにRNA配列決定の結果とすることができる。別の実施形態では、1つ又は複数の標的遺伝子の発現レベルは、タンパク質レベルによって、たとえば標的遺伝子によってコードされたタンパク質の濃度によって、測定することができる。
上記の発現レベルは、任意選択で、この用途によりよく適切であることもないこともある多くの方法で変換することができる。たとえば、4つの異なる発現レベル、たとえばマイクロアレイによるmRNAレベルの変換は、
「連続データ」にすること、すなわち、MAS5.0及びfRMAなどのよく知られているアルゴリズムを使用してマイクロアレイの処理後に得られる発現レベルとすること、
「zスコア」にすること、すなわち、全試料の平均が0になり標準偏差が1になるように基準化された連続発現レベルとすること、
「離散的」にすること、すなわち、ある特定の閾値を超えるすべての発現が1に、それ未満が0に設定される(たとえば、プローブセットの閾値は、いくつかの陽性の臨床試料と同じ数の陰性の臨床試料の組におけるその値の中央値に選ぶことができる)ものとすること、
「ファジー」にすること、すなわち、連続発現レベルが0と1の間の値に、次式
1/(1+exp((thr-expr)/se))
のシグモイド関数を使用して変換されるものとすることができ、ここで、exprは連続発現レベルであり、thrは前述の閾値であり、seは0と1の間の差に影響を及ぼす軟化パラメータである。
構築することができる最も簡単なモデルの1つは、第1の層の転写因子(TF)を表すノードと、たとえばマイクロアレイ又は(q)PCR実験において、第2の層の特定の標的遺伝子と特に高い相関関係がある1つのプローブセットによる、標的遺伝子発現強度レベルの直接測定値を表す重み付けノードとを有するモデルである。重みは、訓練データ・セットからの計算に基づくこと、又は専門的知識に基づくことができる。複数の発現レベルが標的遺伝子ごとに測定される可能性がある場合に(たとえば、1つの標的遺伝子が複数のプローブセットを用いて測定され得るマイクロアレイ実験の場合に)、標的遺伝子ごとに1つの発現レベルしか使用しないというこの手法は、特に単純である。特定の標的遺伝子に使用される1つの発現レベルを選択する具体的な方法は、訓練データ・セットの活性試料と不活性試料を最善に分離できるプローブセットによる発現レベルを使用することである。このプローブセットを決定する1つの方法は、統計的検定、たとえばt検定を行い、最小のp値を持つプローブセットを選択することである。最小のp値を持つ、訓練データ・セットの発現レベルのプローブは、定義により、(既知の)活性及び不活性試料の発現レベルが重なり合う確率が最小のプローブである。もう1つの選択方法は、オッズ比に基づく。このようなモデルでは、1つ又は複数の発現レベルが1つ又は複数の標的遺伝子のそれぞれに与えられ、1つ又は複数の一次結合が、1つ又は複数の標的遺伝子のそれぞれに重み付け項を含む一次結合を含み、各重み付け項は、それぞれの標的遺伝子に与えられた1つ又は複数の発現レベルのうちのただ1つの発現レベルに基づいている。ただ1つの発現レベルが標的遺伝子ごとに上述のように選ばれる場合、そのモデルは「最も判別的なプローブセット」モデルと呼ぶことができる。
「最も判別的なプローブセット」モデルの代替形態では、複数の発現レベルが標的遺伝子ごとに測定される可能性がある場合に、標的遺伝子ごとに与えられるすべての発現レベルを利用することが可能である。このようなモデルでは、1つ又は複数の発現レベルが1つ又は複数の標的遺伝子のそれぞれに与えられ、1つ又は複数の一次結合が、1つ又は複数の標的遺伝子に与えられた1つ又は複数の発現レベルのうちの全発現レベルの一次結合を含む。言い換えると、1つ又は複数の標的遺伝子それぞれについて、それぞれの標的遺伝子に与えられた1つ又は複数の発現レベルのそれぞれに、一次結合においてそれ自体の(個々の)重みによって重み付けすることができる。この変形形態は、「全プローブセット」モデルと呼ぶことができる。これには、与えられた発現レベルすべてを利用しながらも比較的単純であるという利点がある。
上述の両モデルには、「単層」モデルとみなせるものであるという共通点があり、TF要素の活性レベルは、発現レベルの一次結合に基づいて計算される。
TF要素の活性レベルがそれぞれのモデルを評価することによって決定された後、決定されたTF要素活性レベルは、細胞シグナル伝達経路の活性を推測するための閾値とすることができる。このような適切な閾値を計算する方法は、不活性経路を有することが分かっている訓練試料の決定されたTF要素活性レベルwlcと、活性経路を持つ訓練試料とを比較することである。そのようにする、またこれらの群において分散も考慮に入れる方法は、閾値
Figure 0007065609000001
を使用して与えられ、
ここで、σ及びμは訓練試料の標準偏差及び平均値である。少数の試料しか活性及び/又は不活性訓練試料において入手できない場合には、擬似数を計算分散に、2つの群の分散の平均に基づいて加えることができる。
Figure 0007065609000002
ここで、νは、群の決定されたTF要素活性レベルwlcの分散であり、xは正の擬似数、たとえば1又は10であり、nact及びnpasはそれぞれ活性試料及び不活性試料の数である。標準偏差σは、分散νの平方根をとることによって得ることができる。
閾値は、解釈を容易にするために、TF要素wlcの決定された活性レベルから差し引いて、細胞シグナル伝達経路の活性スコアを得ることができ、それにより、負の値が不活性細胞シグナル伝達経路に対応し、正の値が活性細胞シグナル伝達経路に対応するようになる。
説明した「単層」モデルの代替形態として、ある経路の活性信号伝達を実験的に決定することを表す「2層」モデルを使用することができる。すべての標的遺伝子について要約レベルが、一次結合を使用して、その関連するプローブセットの測定強度に基づいて計算される(「第1の(下部)層」)。計算要約値は引き続き、別の一次結合を使用して、経路の他の標的遺伝子の要約値と結合される(「第2の(上部)層」)。重みは、訓練データ・セットから、若しくは専門的知識に基づいて、又はこれらを合わせたものに基づいて、知ることができる。別の言い方をすると、「2層」モデルにおいて、1つ又は複数の発現レベルが1つ又は複数の標的遺伝子のそれぞれに与えられ、1つ又は複数の「一次結合」は、1つ又は複数の標的遺伝子のそれぞれに対し、それぞれの標的遺伝子に与えられた1つ又は複数の発現レベルのうちの全発現レベルの第1の一次結合を含む(第1の(下部)層)。このモデルはさらに、1つ又は複数の標的遺伝子のそれぞれに重み付け項を含む別の一次結合に基づいており、各重み付け項は、それぞれの標的遺伝子の第1の一次結合に基づいている(第2の(上部)層」)。
要約値の計算は、「2層」モデルの好ましいバージョンでは、訓練データを使用して各標的遺伝子の閾値を定義すること、及び計算された一次結合から閾値を差し引いて遺伝子要約を得ることを含む。ここで、閾値は、負の遺伝子要約レベルが下方調節標的遺伝子に対応し、正の遺伝子要約レベルが上方調節標的遺伝子に対応するように選ぶことができる。また、遺伝子要約値が、たとえば上述の変換(ファジー、離散など)の1つを使用して、それらが「第2の(上部)層」)で結合される前に変換されるということも可能性がある。
TF要素の活性レベルが、「2層」モデルを評価することによって決定された後、決定されたTF要素活性レベルは、上述のように、細胞シグナル伝達経路の活性を推測するための閾値とすることができる。
本明細書では、WO2014/102668A2に関連して上述したモデルが「(擬似)線形モデル」と総称される。
数学モデル構築に関する上記の説明は、PI3K経路の活性の推測、標的遺伝子の選択及び訓練にも当てはまり、また数学モデルの使用は、Wnt経路、ER経路、及びHH経路と比較して、PI3K経路に対してある程度修正された。したがって、これらのステップについては、PI3K経路について以下でより詳細に説明する。
(A)PI3K経路
(i)標的遺伝子の選択
転写因子(TF)は、タンパク質複合体(すなわち、特定の構造で結合したタンパク質の複合化したもの)又はタンパク質であり、特定のDNA配列に結合することによって標的遺伝子からの転写を調節し、それによってDNAからmRNAへの遺伝情報の転写を制御することができる。この転写複合体の作用により直接産生されたmRNAは、本明細書では(転写因子の)「直接標的遺伝子」と呼ばれる。細胞シグナル伝達経路活性化によりまた、「間接標的遺伝子」と呼ばれる、さらなる第2の遺伝子転写が生じ得る。以下では、細胞シグナル伝達経路活性とmRNAレベルの間の直接連結として直接標的遺伝子を含む、又は直接標的遺伝子から成るベイジアン・ネットワーク・モデルが(例示的数学モデルとして)選ばれているが、直接標的遺伝子と間接標的遺伝子の差異がいつも明らかであるとは限らない。本明細書では、スコアリング関数を使用して、利用可能な科学文献データに基づき直接標的遺伝子を選択する方法が提示される。それにもかかわらず、限られた情報並びに生物学的バリエーション及び不確実なことがらにより、間接標的遺伝子の偶発的な選択を除外することはできない。標的遺伝子を選択するために、現在利用可能な科学文献の2つのリポジトリが、標的遺伝子の2つのリストを生成するために使用された。
標的遺伝子の第1のリストは、「www.ncbi.nlm.nih.gov/pubmed」でアクセスすることができ、本明細書では別に「Pubmed」と呼ばれる、国立保健研究所のMEDLINEデータベースから検索された科学文献に基づいて生成された。推定FOXO標的遺伝子を含む刊行物は、2013年の第1四半期の期間に、(FOXO AND 「target gene」)などの照会文を使用することによって探索された。得られた刊行物はさらに、以下で詳細に説明する方法論に従って、手作業で解析された。
特定の細胞シグナル伝達経路mRNA標的遺伝子が科学文献から、ランキング・システムを使用することによって選択され、特定の標的遺伝子の科学的証拠が、証拠が累積された科学実験のタイプに応じて格付けされた。一部の実験的証拠は、たとえば、PI3K細胞シグナル伝達軸が活性であることが分かっている細胞株のマイクロアレイ上で増加するmRNAのように、遺伝子が標的遺伝子であることを単に示唆するにすぎないが、別の証拠は、同定された細胞シグナル伝達経路TF結合部位と、細胞内の特定の細胞シグナル伝達経路の刺激後のクロマチン免疫沈降(ChIP)アッセイにおけるこの部位の検索と、細胞株における細胞シグナル伝達経路の特定の刺激後のmRNAの増加とを合わせたもののように、非常に強力になり得る。
特定の細胞シグナル伝達経路標的遺伝子を見つける実験のいくつかのタイプを科学文献で特定することができる。
1.細胞シグナル伝達経路-TFからゲノム上のその結合部位までの直接結合が示される、ChIP実験。例として、クロマチン免疫沈降(ChIP)技術を引き続き使用することによって、PI3K経路の活性誘導がある細胞株とない細胞株のDNAの推定機能FOXO TF結合部位が、純粋にヌクレオチド配列に基づいて認識された結合部位のサブセットとして同定された。推定機能が、TFがDNA結合部位に結合することが見出されたというChIP由来の証拠として同定された。
2.結合配列を含むDNAの断片へのTFの結合を生体外で示す電気泳動移動度シフト(EMSA)アッセイ。ChIPによる証拠と比較して、EMSAによる証拠は、それを生体外状況に変換することができないので、あまり強力ではない。
3.細胞シグナル伝達経路を刺激すること、及びマイクロアレイ上でmRNAプロファイルを測定すること、或いはRNA配列決定を使用すること、細胞シグナル伝達経路誘導性細胞株を使用すること、及び誘導後のいくつかの時点で測定されたmRNAプロファイルを、タンパク質への翻訳を阻害するシクロヘキシミドの存在下で測定すること。したがって誘導されたmRNAは直接標的遺伝子であると考えられる。
4.3と同様であるが、mRNAの量を測定するために定量PCRを使用すること。
5.バイオ・インフォマティクス手法を使用してゲノム内のTF結合部位を同定すること。FOXO TF要素の例として、保存されたFOXO結合モチーフ5’-TTGTTTAC-3’を使用して、ヒトゲノム配列に対しソフトウェア・プログラムが実行され、潜在的な結合部位が、遺伝子プロモータ領域でも他の遺伝子領域でも同定された。
6.3と同様であるが、ただシクロヘキシミドがない状態で。
7.4と同様であるが、ただシクロヘキシミドがない状態で。
8.細胞シグナル伝達経路が活性であることが分かっている特定の組織又は細胞試料のmRNA発現プロファイリングであるが、適切なネガティブ・コントロール条件がない状態で。
最も単純な形で、標的mRNAが同定されたこれらの実験的手法のそれぞれについて、すべての潜在的標的mRNAに1ポイントを与えることができる。
或いは、ポイントを漸増的に与えて、1つの技術が1ポイントになるようにすることができ、第2の技術が第2のポイントを付加し、以下同様とすることもできる。この比較的単純なランキング方策を使用して、最も信頼できる標的遺伝子のリストを作ることができる。
或いは別の方法のランキングを使用して、生体外直接標的遺伝子に最も多くの証拠をもたらす技術により多くのポイント数を与えることによって、直接標的遺伝子である可能性が最も高い標的遺伝子を同定することができる。上記のリストでは、これは、実験手法1)では8ポイント、2)では7ポイントになり、実験手法8)では1ポイントまで下がっていくということになる。このようなリストは、「一般的標的遺伝子リスト」と呼ぶことができる。
生物学的バリエーション及び不確実性にもかかわらず、本発明者らは、直接標的遺伝子が、組織に依存せずに誘導される可能性が最も高いと仮定した。これらの標的遺伝子のリストは「標的遺伝子の証拠精選リスト」と呼ぶことができる。このような標的遺伝子の証拠精選リストを使用して、様々な組織源から来る試料に適用できるPI3K経路の計算モデルを構築した。
以下では、どのようにして証拠精選標的遺伝子リストの選択がPI3K経路について具体的に構築されたかを例示的に示す。
「モデル」の入力として使用されるPI3K標的遺伝子を選択する目的のために、以下の3つの基準が使用された。
1.遺伝子プロモータ/エンハンサ領域がFOXO結合モチーフを含む。
a.FOXO結合モチーフは、PI3K経路の活性に反応することが、たとえば、特定のFOXOモチーフがレポーター遺伝子に連結される一過性トランスフェクション・アッセイによって証明されなければならない。
b.FOXOモチーフの存在が、たとえば、遺伝子プロモータ/エンハンサ領域の豊富なモチーフ解析によって確認されなければならない。
2.FOXOは、問題となっている遺伝子のプロモータ/エンハンサ領域に生体外で(差次的に)結合し、これはたとえば、ChIP/CHIP実験又は別のクロマチン免疫沈降技法によって実証される。
a.FOXOは、PI3K経路が活性ではないとき、遺伝子のプロモータ/エンハンサ領域に結合することが証明され、また
b.PI3K経路が活性であるとき、(好ましくは)遺伝子のプロモータ/エンハンサ領域に結合しない(又は弱く結合する)。
3.遺伝子は、PI3K経路の活性が変化するとき差次的に転写され、これはたとえば、
a.リアルタイムPCR若しくはマイクロアレイ実験による、問題となっている遺伝子のmRNAの折りたたみ強化、又は
b.免疫沈降アッセイによる、RNA Pol IIが遺伝子のプロモータ領域に結合することの実証、
によって実証される。
選択は、十分かつ適切に記述された実験的証拠が集められて、上述の3つすべての基準が満たされたことを証明された遺伝子を、PI3K標的遺伝子として定義することによって行われた。PI3Kの差次的結合の証拠を集めるための適切な実験は、たとえば、タモキシフェンに曝露されたとき又は曝露されないときに、タモキシフェンに反応してPI3K経路の活性を発現する癌細胞株(たとえば、FOXO.A3.ERなどの、タモキシフェン誘導性FOXO構造物を形質移入された細胞株)におけるChIP-Seq実験の結果を比較することである。mRNA転写の証拠を集めることについても同様である。
上記では、上述の手法を使用して見つかった証拠に基づいていくつかの標的遺伝子を選択するために使用された、標的遺伝子選択手順の一般的手法及びより具体的な例について論じている。PI3K経路についてベイジアン・ネットワーク・モデルに使用された標的遺伝子のリストは、表1に示されている。
Figure 0007065609000003
Figure 0007065609000004
標的遺伝子の第2のリストは、Thomson-ReutersのMetacore(最終アクセス2013年5月14日)の中に提供されている、手作業で精選された科学文献のデータベースを使用して生成された。このデータベースは、ヒトFOXO転写因子(すなわち、FOXOl、FOX03A、FOX04及びFOX06)のファミリーの下流で直接転写調節される遺伝子に関して照会された。この照会により、336個の推定FOXO標的遺伝子が得られ、これらはさらに以下のように解析された。まず、支持する刊行物が1つしかない推定FOXO標的遺伝子がすべて取り除かれた。次に、刊行物で報告された、ChIP、EMSA、差次的発現、ノックダウン/アウト、ルシフェラーゼ遺伝子レポーターアッセイ、配列解析などの実験的証拠のタイプごとにポイントを与えるスコアリング関数が導入された。同一の実験的証拠が場合により複数の刊行物に記載されていて、その結果、対応するポイント数が得られ、たとえば、2つの刊行物が1つのChIP結果について述べていると、単一のChIP結果に与えられるスコアの2倍になる。さらなる解析が行われて、1つのタイプの実験的証拠しかない、たとえば差次的発現しかないのではなく、多様なタイプの実験的証拠がある遺伝子だけが認められた。最後に、全推定FOXO標的遺伝子について証拠スコアが計算され、6つ以上の証拠スコアがある推定FOXO標的遺伝子がすべて選択された(表2に示す)。6つというカットオフレベルは、経路活性を決定するにはおおよそ30個の標的遺伝子でほぼ十分であることが事前に明らかになっているので、ヒューリスティックに選ばれた。
これらの標的遺伝子のリストは「標的遺伝子のデータベースによるリスト」と呼ぶことができる。このような精選標的遺伝子リストを使用して、異なる組織源から来る試料に適用できる計算モデルを構築した。
Figure 0007065609000005
標的遺伝子の第3のリストは前述の2つのリスト、すなわち証拠精選リスト(表1参照)、及びデータベースによるリスト(表2参照)に基づいて生成された。これら2つのリストから遺伝子をさらに選択するために、3つの基準が使用された。第1の基準は、標的遺伝子に起因する機能に関連する。遺伝子に起因する機能は科学文献に見出すことができるが、NIHのOMIMデータベース(「http://www.ncbi.nlm.nih.gov/omim」によってアクセス可能)などの公共データベースで入手できることが多い。アポトーシス、細胞周期、腫瘍抑制/進行、DNA修復、分化などの、癌にとって極めて重要なプロセスに関与することに帰することが分かった、表1の証拠精選リスト及び表2のデータベースによるリストからの標的遺伝子が、第3のリストに選択された。最後に、既知の低PI3K/高FOXO活性に対して既知の高PI3K/低FOXO活性を用いた細胞株実験で、高い差次的発現があることが分かった標的遺伝子が選択された。本明細書では、複数試料について平均されたFOXO転写の「オン」と「オフ」の状態間に(本明細書では、プローブセット・レベルで)20.5の最小発現差がある標的遺伝子が第3のリストに含められた。第3の基準は、最も判別可能な標的遺伝子を選択することが特に狙いとされた。既知の高PI3K/低FOXO活性を用いた複数の試料の細胞株実験、及び既知の低PI3K/高FOXO活性を用いた複数の試料の細胞株実験における発現レベルに基づいて、オッズ比(OR)が計算された。本明細書では、オッズ比は、測定の不確実性をあらわすカットオフ及び軟境界として中央値を使用して、プローブセットごとに計算された。証拠精選リスト及びデータベースによるリストの標的遺伝子は、「軟」オッズ比に従って順位付けされ、最高順位(OR>2)及び最低順位(OR<1/2、すなわち負に調節された標的遺伝子)の標的遺伝子が、標的遺伝子の第3のリストに選択された。
遺伝子の機能、「オン」対「オフ」シグナル伝達における差次的発現、及び高オッズ比を考慮に入れると、PI3Kシグナル伝達経路の活性の決定の際により証拠となると考えられた標的遺伝子のセットが見つかった(表3に示す)。このような標的遺伝子のリストは、「標的遺伝子のショートリスト」と呼ぶことができる。したがって、表3に報告された標的遺伝子は、本発明によれば特に好ましい。それにもかかわらず、マイクロアレイなどの取得技術により遺伝子の大きい組の発現レベルを比較的容易に取得できることをかんがみて、表3の標的遺伝子の一部又は全部を利用すること、並びに任意選択で、表1及び表2の残りの標的遺伝子のうちの1つ、2つ、一部又は全部を付加的に使用することが企図されている。
Figure 0007065609000006
標的遺伝子の証拠精選リスト(表1参照)の別の選択が、本発明者らによって行われた。訓練試料からPI3K経路の活性を決定する上でより証拠となることが判明した証拠精選リストの標的遺伝子が選択された。標的遺伝子が、「軟」オッズ比に基づいて選択され、上位12個の遺伝子が「12標的遺伝子ショートリスト」(表4参照)に選択された。
Figure 0007065609000007
(ii)数学経路モデルの訓練及び使用
数学経路モデルを使用して対象の細胞シグナル伝達経路の、ここではPI3K経路の、活性を推測することができるようにするには、モデルが適切に訓練されなければならない。
数学経路モデルが、FOXO TF要素に関連する条件付き確率と、対象の試料で測定されたPI3K経路の3つ以上の標的遺伝子の発現レベルとに基づく、確率的なモデル、たとえばベイジアン・ネットワーク・モデルである場合、訓練は、好ましくは、公開されている国際特許出願WO2013/011479A2(「Assessment of cellular signaling pathway activity using probabilistic modeling of target gene expression」)に詳細に記載されているように実施することができる。
数学経路モデルが、対象の試料で測定されたPI3K経路の3つ以上の標的遺伝子の発現レベルの1つ又は複数の一次結合に基づいている場合、訓練は、未公開米国特許仮出願第61/745839号(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)に詳細に記載されているように実施することができる。
本明細書では、図1に示された例示的なベイジアン・ネットワーク・モデルが使用されて、PI3K経路の転写プログラムが簡単にモデル化された。モデルは、3つのタイプのノードから成る。すなわち、(a)第1の層1としての転写因子(TF)要素(状態「不在」及び「存在」で)、(b)第2の層2としての標的遺伝子TG、TG、TG(状態「上」及び「下」で)、並びに第3の層3としての(c)標的遺伝子の発現レベルと関連付けられた測定ノード、である。これらは、本明細書で以前に用いられたように、マイクロアレイプローブセットPS1,1、PS1,2、PS1,3、PS2,1、PSn,1、PSn,m(状態「低」及び「高」で)とすることができるが、RNAseq又はRT-qPCRなどの他の遺伝子発現測定値とすることもできる。
数学経路モデルの、本明細書では例示的ベイジアン・ネットワーク・モデルの、適切な実施態様はマイクロアレイ・データに基づいている。このモデルは、(i)標的遺伝子の発現レベルがどのようにTF要素の活性化に依存するか、及び(ii)プローブセット強度が、ひいては、どのようにそれぞれの標的遺伝子の発現レベルに依存するか、について記述する。後者については、プローブセット強度は、Gene Expression Omnibus(GEO, www.ncbi.nlm.nih.gov/geo)及びArrayExpress(www. ebi.ac.uk/arrayexpress)から広範に入手可能である、fRMA事前処理アフィメトリクスHG-U133Plus2.0マイクロアレイにより取得することができる。
例示的なベイジアン・ネットワーク・モデルは細胞シグナル伝達経路の、本明細書ではPI3K経路の、生物学の単純化であるので、また生物学的測定には一般にノイズが伴うので、確率論的アプローチが選択される。すなわち、(i)TF要素と標的遺伝子、及び(ii)標的遺伝子とそれらの各プローブセット、の間の関係が確率論的用語で記述される。さらに、腫瘍増殖を促進する癌化細胞シグナル伝達経路の活性は、一時的及び動的に変化するのではなく、長期に、さらには不可逆的にも変化すると仮定された。したがって、例示的なベイジアン・ネットワーク・モデルは、静的細胞状態の解釈のために開発された。この理由のために、複雑な動的細胞シグナル伝達経路特性はモデルに組み込まれなかった。
例示的なベイジアン・ネットワーク・モデルが構築され較正されると(下記参照)、モデルは、第3の層3における観察結果としてのプローブセット測定値を入力することによって、またTF要素が「存在」する確率がどうなるはずであったかをモデルで逆に推測することによって、新しい試料のマイクロアレイ・データに対して使用することができる。ここで、「存在」とは、TF要素がDNAに結合し、細胞シグナル伝達経路の標的遺伝子の転写を制御している現象であると考えられ、「不在」とはTFエレメントが転写を制御していない場合であると考えられる。したがって、この後者の確率は、細胞シグナル伝達経路の、本明細書ではPI3K細胞シグナル伝達経路の、活性を示すために使用できる一次読み出しであり、これは次に、活性であることと不活性であることの確率の比をとることによって、細胞シグナル伝達経路が活性であるオッズに変換することができる(すなわち、pを細胞シグナル伝達経路が活性であることの予測確率とすれば、オッズはp/(1-p)で与えられる)。
例示的なベイジアン・ネットワーク・モデルでは、確率論的関係が、定量的な確率的推論を可能にするために定量的になっている。組織型全体にわたって一般化挙動を改善するために、(i)TF要素と標的遺伝子、の間の確率論的関係を記述するパラメータが慎重に精選された。TF要素が「不在」である場合、標的遺伝子は「下」になる可能性が最も高く、したがって、これには0.95の確率が選ばれ、標的遺伝子が「上」になることには0.05の確率が選ばれる。後者(非ゼロ)の確率は、標的遺伝子が他の因子によって調節される、又は(たとえば測定ノイズの故に)偶発的に「上」と観察される(まれな)可能性に相当するものである。TF要素が「存在」する場合、0.70の確率で標的遺伝子は「上」と考えられ、0.30の確率で標的遺伝子は「下」と考えられる。後者の値は、TF要素が存在しても標的遺伝子が高度に発現しないいくつかの理由があり得るので、たとえば、遺伝子のプロモータ領域がメチル化されるので、このように選ばれる。標的遺伝子がTF要素によって上には調節されないが下には調節される場合には、確率は同じようにして選ばれるが、TF要素が存在することに対し下方調節を反映する。(ii)標的遺伝子とそれらの各プローブセット、の間の関係を記述するパラメータは、実験データに対して較正された。後者に関して、この例では、定義された活性及び不活性経路設定を用いた細胞株実験からのマイクロアレイ・データが使用されたが、これは、既知の細胞シグナル伝達経路活性状態の患者試料を使用して行うこともできる。得られた条件付き確率表を以下に示す。
Figure 0007065609000008
これらの表では、変数ALi,j、AHi,j、PLi,j、及びPHi,jは、「低」(L)又は「高」(H)プローブセット強度をそれぞれ有する「不在」(A)又は「存在」(P)転写複合体がある較正試料の数を示す。0及び1の極端な確率を回避するために、ダミー数が加えられた。
観測されたプローブセット強度を離散化するために、各プローブセットPSi,jに対し閾値ti,jが使用され、これに満たない観察値は「低」と呼ばれ、これを超える観察値は「高」と呼ばれる。この閾値は、使用された較正データセット中のプローブセットの(重み付け)中心強度になるように選ばれた。マイクロアレイ・データにノイズが伴うことにより、ファジー法が、観察プローブセット強度をその閾値と比較するときに、報告強度まわりで0.25(log2スケールで)の標準偏差による正規分布を仮定することによって、また閾値の下及び上の確率質量を決定することによって、使用された。
上述の例示的なベイジアン・ネットワークの代わりに上述の(擬似)線形モデルが使用された場合には、モデルを使用して試験試料の細胞シグナル伝達経路活性を推測できるようにするには、ノードと、ノードが「不在」であるか、それとも「存在」すると呼ぶための閾値との間の相関性の符号及び大きさを示す重みが決定される必要がある。専門的知識を用いてこれらの重み及び閾値を先験的に埋めることができるが、一般には、モデルは訓練試料の代表的な組を使用して訓練され、その好ましくはグランド・トルースが、たとえば、既知の「存在」転写因子複合体(=活性細胞シグナル伝達経路)又は「不在」転写因子複合体(=不活性細胞シグナル伝達経路)を有する試料中のプローブセットの発現データが、知られている。
当分野では、モデル・トポロジを考慮に入れる、またモデル・パラメータ、ここでは重み及び閾値を、モデル出力、ここでは重み付け線形スコアが最適になるように変更する、多数の訓練アルゴリズム(たとえば、回帰)が知られている。或いは、観察された発現レベルから直接重みを計算することも、最適化アルゴリズムを必要とせずに、可能である。
第1の方法、本明細書では「白黒」法という名称、は煎じ詰めると3進法であり、各重みは組{-1,0,1}のうちの1つの要素になる。これが生物学的文脈において言われる場合には、-1及び1が、細胞シグナル伝達経路活性の場合にそれぞれ下方及び上方調節される標的遺伝子又はプローブセットに対応する。プローブセット又は標的遺伝子が上方調節されるか、又は下方調節されるか統計的に証明され得ない場合には、プローブセット又は標的遺伝子は0の重みを受け取る。1つの例では、左側及び右側の、活性細胞シグナル伝達経路試料の発現レベル対不活性細胞シグナル伝達経路を有する試料の発現レベル、の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,Vol.74,No.11,2014年、2936~2945頁、に見出すことができる。
本明細書では、40HTにより刺激すると誘導可能であるFOXOコンストラクトの安定したトランスフェクションを用いたHUVEC細胞の発現、についての公的に入手可能なデータ(Gene Expression Omnibusから入手可能なGSE16573、最終アクセス2014年10月6日)が、PI3K経路モデルを訓練する一例として使用された。40HTにより12時間刺激された、誘導可能なFOXOコンストラクトを持つ細胞株はFOXO活性試料(n=3)とみなされたのに対し、不活性FOXO試料は、40HT刺激なしのコンストラクトを持つ細胞株であった(n=3)。
(B)Wnt経路
Wnt標的遺伝子の選択については以前に、WO2013/011479A2及びWO2014/102668A2に記載されている。Wnt経路の「標的遺伝子の証拠精選リスト」(表5参照)は、Wnt経路の「標的遺伝子のショートリスト」(表6参照)及びWnt標的遺伝子の「12標的遺伝子ショートリスト」(表7参照)を生成するために、上記のこの実施例でPI3K標的遺伝子について説明されているように使用された。
Figure 0007065609000009
Figure 0007065609000010
Figure 0007065609000011
(C)ER経路
WO2013/011479A2及びWO2014/102668A2と比較して、本明細書では、新しい文献証拠が追加されたので、ER標的遺伝子の順位がわずかに変えられていることに留意されたい。ER標的遺伝子は、WO2014/102668A2の実施例3に記載されたのと同じようにして選択され順位付けされた。遺伝子は、文献証拠スコアと各遺伝子の個々の機能とを組み合わせてモデル内の活性経路と不活性経路とを区別することによって、順位付けされた。この順位付けは、エストロゲンを枯渇させ、引き続き枯渇させたままにされた、又は1nMエストロゲンに24時間曝露された(GSE35428)MCF7細胞株試料の訓練セットを用いてモデルを訓練したときに、また、モデルをその訓練セットと、MCF7細胞がエストロゲンを枯渇させ、引き続き枯渇させたままにされた、又は10nM若しくは25nMエストロゲン(それぞれ、GSE11352及びGSE8597)に曝露された別の2つの訓練セットとを用いて試験したときに遺伝子ごとに得られた、重み付けされた偽正比率と偽負比率の一次結合に基づいた。
(重み付けされた偽正数と偽負数の組合せ(オッズ比の代わり)は、様々なセットに使用された異なる実験条件に相当するように使用されたことに留意されたい。異なる重みは、偽正数(負数)はモデルの結果であり、異なる実験条件に試料が置かれた結果ではなかったという発明者の確信に従って、設定された。たとえば、すべての実験において、MCF7細胞株試料は、まずエストロゲンをある期間枯渇させた後にエストロゲンに曝露されるか、又はさらにもう24時間枯渇させた。枯渇時間が短いと、エストロゲン枯渇にもかかわらず経路が活性のままになっている可能性があり、この場合には、試験試料と訓練試料の両方が同じ時間枯渇されたときよりも、偽正数は重みが小さい。)
以下でさらに詳細に論じる、さらなる文献調査と、活性試料と不活性試料の間の差次的発現の大きさの検査とに基づいて、PDZK1がER経路の直接標的遺伝子として選択された。この(PI3Kの)実施例で説明した類似の方法論を用いて、推定ER標的遺伝子の実験的証拠に関する追加の科学論文を手作業で評価した後、いくつかの追加の推定ER標的遺伝子が同定された。
推定ER標的遺伝子が、エストロゲン応答要素(ERE)モチーフを含む遺伝子プロモータ/エンハンサ領域の存在に関して、解析された。EREモチーフは、たとえば、特定のEREモチーフがレポーター遺伝子と関連付けられる一過性トランスフェクション・アッセイによって、エストロゲンに応答することが証明されなければならない。EREモチーフの存在は、たとえば、遺伝子プロモータ/エンハンサ領域の強化モチーフ解析によって確認されなければならない。加えて、ERは、問題となっている遺伝子のプロモータ/エンハンサ領域に生体外で(差次的に)結合し、これはたとえば、ChIP/CHIP実験又は別のクロマチン免疫沈降アッセイによって実証される。たとえば、ERは、ER経路が活性であるとき、遺伝子のプロモータ/エンハンサ領域に結合することが証明されるはずであり、また、たとえば、ER経路が活性でない場合は、遺伝子のプロモータ/エンハンサ領域に結合しない(又は弱くしか結合しない)。最後に、遺伝子は、ER経路が活性であるとき差次的に転写され、これはたとえば、リアルタイムPCR若しくはマイクロアレイ実験による、問題となっている遺伝子のmRNAの折りたたみ強化によって、又は、免疫沈降アッセイによる、RNA Pol IIが遺伝子のプロモータ領域に結合することの実証によって、実証される。
選択は、十分かつ適切に記述された実験的証拠が文献から集められ、上述の3つすべての基準が満たされたことを証明された遺伝子をER標的遺伝子として定義することによって行われた。ERの差次的結合の証拠を集めるための適切な実験は、たとえば、エストロゲンに曝露されたとき又は曝露されないときに、エストロゲンに反応する癌細胞株(たとえば、MCF-7細胞株)についての、たとえばChIP/CHIP実験の結果を比較することである。すべての追加の科学論文を評価した後、すべての推定標的遺伝子の新しい順位付けが、文献中に見出された実験的証拠の強さに基づいた。その結果、ER細胞シグナル伝達経路の1つの推定標的遺伝子PDZK1が、設定閾値を超える実験的証拠スコアに達した。したがって、PDZK1がER経路の真の直接標的遺伝子であるとみなされた。
ER標的遺伝子の当初の選択では、「軟」オッズ比を用いて計算された、不活性試料に対し活性試料を区別する能力だけが考慮された。今の解析では、差次的発現の大きさも評価に含まれた。差次的発現シグナルの大きさは、適切に設計されたアッセイの重要な特徴として「軟」オッズ比に次ぐので、この新しい選択方法は、当初の基準に比べて改善になるものと期待される。差次的遺伝子発現の大きさは、アフィメトリクスHG1133Plus2データセット、すなわちGSE35427、GSE11352、GSE21618、GSE8597と、エストラジオール(E2)又はコントロールで刺激された複数の乳癌細胞株を含む2つの社内生成データベースとから選択したものについて、ER活性(オン)試料とER不活性(オフ)試料の間の平均遺伝子発現の差を平均することによって推定された。平均遺伝子発現は、それぞれのデータセットの遺伝子に関連したアフィメトリクス・プローブセットごとに別々に計算された。顕著に差次的に発現されたプローブセットだけが、平均に関して考慮された。PDZK1の、エストラジオールで刺激された試料すなわちER活性試料と、コントロール/未刺激試料すなわちER不活性試料との間の平均の差次的発現は2.08であった。この差次的発現は例外的に高く(上方調節された遺伝子全体としての平均は0.88)、差次的発現が最も高い標的遺伝子、たとえば平均の差次的発現が2.14であるPGRに匹敵する。加えて、PDZK1の「軟」オッズ比(平均26.6)も、平均(19.03)より高い。
以下の例において、我々は、元の13ER標的遺伝子リスト(GREBl、PGR、XBP1、CA12、SOD1、CTSD、IGFBP4、TFF1、SGK3、NRIP1、CELSR2、WISP2、及びAP1B1)モデル(以下ではショート・リスト・モデルと呼ばれる)を、PDZK1及び元の13ER標的遺伝子リストを用いて構築した新しい14ER標的遺伝子モデル(以下ではショート・リスト+PDZK1モデルと呼ばれる)と比較した。両方のベイジアン・ネットワーク・モデルが、違いはER標的遺伝子のリストのみで、全く同じようにして(アフィメトリクスHGU133Plus2 GSE8597データセットを使用して)訓練された。
例1では、ER経路活性が、アフィメトリクスHGU133Plus2データセットから選択したものについて計算された。このデータセットは、256個のER陽性乳癌試料、195個のER陰性乳癌試料、27個の正常な乳房組織試料、及び94個の未知のER状態乳癌試料を含む、典型的な乳癌及び正常な乳房の組織試料(公開データセットGSE12276、GSE10870、及びGSE21653)を例証するものである。ER経路は、ER陰性乳癌及び正常乳房では不活性であると予想されるのに対し、ER陽性乳癌の約50~70%は、ホルモン療法への反応に基づき、活性であると予想される。ショート・リスト・モデルにより(74%が)、またショート・リスト+PDZK1により(73%が)活性であると予測されるER陽性乳癌試料の比率は、ホルモン療法に反応するER陽性癌患者の比率に匹敵し、また類似している。さらに、ショート・リスト+PDZK1(平均log2オッズ比が2.73)リスト・モデルによって、ER陽性試料全体として計算されたER活性化の確率の平均は、ショート・リスト・モデルによって予測された活性化の平均確率よりもわずかに高くなって(平均log2オッズ比が2.70であり、log2オッズ比スケールで0.03の差がある)、この種の試料に匹敵するものになる。PDZK1を含めることの予想外の有益な技法的効果が、ER陰性乳癌試料及び正常組織試料を解析するときに生じる。すなわち、ショート・リスト+PDZK1リスト・モデルによって計算されたER活性化の確率の平均(平均log2オッズ比が-7.3)は、ショート・リスト・モデルによって予測された活性化の平均確率よりもかなり低くなって(平均log2オッズ比が6.8で、log2オッズ比スケールで0.5の差があり、ウィルコクソン順位検定2面pv=0.02)、この状況ではショート・リスト+PDZK1モデルがショート・モデルよりも技法的に優れたものになる。さらに、この改善は、1つ又は複数の標的遺伝子がモデルに追加される場合に予期される予測経路活性の些細なスケーリングよりも大きく、したがって、PDZK1の追加は、予想外の有利な技法的効果をもたらす。
例2では、ER経路活性が、公開アフィメトリクスHGU133Plus2データセットGSE8597、GSE35428、GSE11352について計算された。このデータセットは、エストロゲン感受性乳房細胞株(この場合ではMCF7)がエストロゲン刺激に曝露される又はエストロゲン刺激が与えられない実験を例証するものである。エストロゲンに曝露されるとMCF7細胞株のER経路が活性化し、エストロゲンが与えられないとMCF7細胞株のER経路が遮断されることがよく知られている。また、この場合、ショート・リスト+PDZK1モデルはショート・リスト・モデルよりも、MCF7細胞株がエストロゲンに曝露される場合でも、ショート・リスト+PDZK1モデルによって計算された予測活性(平均log2オッズ比が14.7)がショート・リスト・モデルによって計算された予測活性(平均log2オッズ比が14.0であり、log2オッズ比スケールで0.7の差)よりも高い場合でも、技法的に優れているように見える。エストロゲン刺激が与えられないすべての試料についてショート・リスト+PDZK1モデルによって計算された予測活性(平均log2オッズ比が-7.7)は、エストロゲンが与えられなかった27個の試料の85%についてショート・リスト・モデルによって計算された予測活性(平均log2オッズ比が-7.3であり、log2オッズ比スケールで0.4の差)よりも低い。また、この改善は、1つ又は複数の標的遺伝子がモデルに追加される場合に予期される予測経路活性の些細なスケーリングよりも大きく、したがって、PDZK1の追加は、予想外の有利な技法的効果をもたらす。
PCRアッセイにおいて新しい遺伝子の効果を究明するために、以下の例で我々は、11ER標的遺伝子リスト(GREB1、PGR、XBP1、CA12、SOD1、CTSD、IGFBP4、TFF1、SGK3、NRIP1、CELSR2、ERBB2、and ESR1)モデル(以下ではPCRリスト・モデルと呼ばれる)を、PDZK1及び上述の11ER標的遺伝子リストを用いて構築した新しい12ER標的遺伝子モデル(以下ではPCRリスト+PDZK1モデルと呼ばれる)と比較した。両方のベイジアン・ネットワーク・モデルが、違いはPCRリスト+PDZK1モデルにPDZK1 ER標的遺伝子を追加することのみで、全く同じようにして(MCF7細胞株についての社内のエストロゲン欠乏/刺激実験から、RT-qPCRによって生成された遺伝子発現データを使用して)訓練された。ER経路活性は、6個がエストロゲンを与えられず、6個がエストロゲンで刺激された合計12個の試料について計算された。ここでもまた、PDZK1を含むモデル(PCRリスト+PDZK1モデル)は、PDZK1がないモデル(PCRリスト・モデル)よりも、エストロゲンに曝露される場合でも、PCRリスト+PDZK1モデルによって計算された予測活性(平均log2オッズ比が4.7)がPCRリスト・モデルによって計算された予測活性(平均log2オッズ比が3.9であり、log2オッズ比スケールで0.8の差)よりも高い場合でも、技法的に優れているように見える。エストロゲンが与えられない試料についてPCRリスト+PDZK1モデルによって計算された予測活性(平均log2オッズ比が-5.1)は、ショート・リスト・モデルによって計算された予測活性(平均log2オッズ比が-4.5であり、log2オッズ比スケールで0.6の差)よりも低い。この差は、少量の「プローブ」を使用して試料ER標的遺伝子プロファイルを測定するモデルにおいては、プローブが有する識別力が通常は小さいので(低い平均予測活性に留意されたい)、非常に重要である。結論として、この改善は、1つ又は複数の標的遺伝子がモデルに追加される場合に予期される予測経路活性の些細なスケーリングよりも大きく、したがって、PDZK1の追加は、予想外の有利な技法的効果をもたらす。
上で論じたように、ER標的遺伝子の選択については、以前にWO2013/011479A2及びWO2014/102668A2に記載されている。HH経路の「標的遺伝子の証拠精選リスト」は、ER経路の「標的遺伝子のショートリスト」及びER標的遺伝子の「12標的遺伝子ショートリスト」を生成するために、さらなる文献調査及びPDZK1標的遺伝子を含めることに基づいて、上記のこの実施例でPI3K標的遺伝子について説明されているように使用された。
Figure 0007065609000012
Figure 0007065609000013
Figure 0007065609000014
(D)HH経路
HH標的遺伝子の選択については、以前にWO2013/011479A2及びWO2014/102668A2に記載されている。HH経路の「標的遺伝子の証拠精選リスト」は、HH経路の「標的遺伝子のショートリスト」及びHH標的遺伝子の「12標的遺伝子ショートリスト」を生成するために、上記のこの実施例でPI3K標的遺伝子について説明されているように使用された。
Figure 0007065609000015
Figure 0007065609000016
Figure 0007065609000017
実施例2 リスク・スコアの決定
一般に、対象が臨床事象を特定の期間内に経験するリスクを示す、また対象における2つ以上の細胞シグナル伝達経路の推測活性を合わせたものに基づく、リスク・スコアを決定するための多くの異なる式を考案することができる。すなわち、
MPS=F(P)+X、ここでi=1...N (3)
ここで、MPSはリスク・スコアを表し(「MPS」という用語は、本明細書では、リスク・スコアが2つ以上の細胞シグナル伝達経路の推測活性の影響を受けることを表すための「複数経路スコア」の略語として用いられる)、Pは細胞シグナル伝達経路iの活性を表し、Nは、リスク・スコアを計算するために使用される細胞シグナル伝達経路の総数を表し、Xは、式に入る可能性がある、あり得る別の要素及び/又はパラメータのためのスペースホルダである。このような式は、より具体的には、所与の変数に関するある次数の多項式にも、その変数の一次結合にもなり得る。このような多項式における重み付け係数及び冪は専門知識に基づいて設定することができるが、一般には、既知のグランド・トルースを用いて設定された訓練データ、たとえばサバイバル・データを使用して式(3)の重み付け係数及び冪の推定値を得る。推測活性は、式(3)を使用して集約され、続いてMPSを生成する。次に、スコアリング関数の重み付け係数及び冪は、高いMPSが、患者が臨床事象を経験する高い確率と相関するように、また逆も同様に、最適化される。サバイバル・データとのスコアリング関数の相関性の最適化は、多数の解析技法を使用して、たとえば、コックス比例ハザード検定(本明細書で以前に用いた)、対数順位検定、傾斜下降又は手動適応などの標準的な最適化技法と併せたカプラン-マイヤー推定量、などを使用して行うことができる。
これらの実験において、本発明者らは、細胞シグナル伝達経路の活性と再発リスクの間に冪法則応答を予期する理由を見出さなかった。したがって、式(3)は次のように簡略化することができる。
MPS=w・P+w・P+...+w・P+X (4)
ここで、w,...,wは重み付け係数を表す。
この実施例では、臨床事象は癌、特に乳癌であり、PI3K経路、Wnt経路、ER経路、HH経路の各推測活性が、本明細書、並びに公開国際特許出願WO2013/011479A2(「Assessment of cellular signaling pathway activity using probabilistic modeling of target gene expression」)及び/又は公開国際特許出願WO2014/102668A2(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)で詳細に論じられているように、考察されている。
本明細書で好ましくは使用される式は、PI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上との活性を考慮に入れている。これらの式は、癌生物学研究から、並びに、公的に利用可能なデータセットにおいて生存と、PI3K経路、Wnt経路、ER経路及びHH経路の各活性との間に本発明者らによって発見された相関関係から導出された、本発明者らの観察結果に基づいている。Wnt経路及びHH経路のような初期発生の経路は、癌幹細胞と呼ばれる、表現型のようにさらに幹細胞に戻った癌細胞によって引き起こされる転移において役割を果たすと考えられる。実のところ、本発明者らは、転移癌細胞がシーディング位置において別の器官又は組織へと分割を開始できるようにする、癌転移における役割を果たすWnt経路などの初期発生の経路の十分な指標が得られると考える。転移は、より不良な予後に関連していると共に癌再発の一形態を表し、したがって、Wnt経路及びHH経路などの癌細胞における初期発生の経路の活性が、より悪い予後を予測するものになると本発明者らに期待された。癌進行及び転移におけるWnt経路及びHH経路の推定された役割は、前臨床研究に基づいており、その活性を測定する方法が得られていなかったので、研究対象として提示されてこなかった。加えて、本発明者らは、生存のための(比較的)防御的な機構であるER経路の活性と、より悪い予後と相互に関連があるPI3K経路の活性との間の相関関係を示す十分な指標を、公的に利用可能なデータセット中に発見した。したがって、ER経路の不活性とPI3K経路の活性とが、乳癌患者の不良な転帰と相互に関連があることが本発明者らによって見出された。
生物学研究からの本発明者らのこれらの観察結果、並びにPI3K経路、Wnt経路及びHH経路の活性が癌再発及び生存全体に役割を果たし得ることと、ER経路の活性が良好な臨床成績につながるように見えることとの臨床相関性が、本明細書では以下の好ましい式として組み合わされており、これは式(4)の特別な場合である。
MPS=w・P+w・P+w・P+w・P+X (5)
ここで、P、P、P、及びPは、それぞれPI3K経路、Wnt経路、ER経路及びHH経路の推測活性を示し(たとえば、0と1の間の範囲で)、wは正の一定重み付け係数であり、w及びwは非負の一定重み付け係数であり、wは非正の一定重み付け係数である。この式によると、対象が臨床事象を特定の期間内に経験する示されたリスクは、合計値の増加と共に単調に増加する。
以下の例では、本発明者らは例示的に、表3に示された標的遺伝子のショートリスト及び本明細書で論じた訓練を使用するPI3K経路と、WO2013/011479A2の表1に示された標的遺伝子の証拠精選リスト及び本明細書で論じた訓練を使用するWnt経路と、WO2013/011479A2の表2に示された標的遺伝子の証拠精選リスト及び本明細書で論じた訓練を使用するER経路と、WO2013/011479A2の表3に示された標的遺伝子の証拠精選リスト及び本明細書で論じた訓練を使用するHH経路とからなるベイジアン・ネットワークからの推測活性を使用した。別法として、経路活性は、本明細書及びより詳細にWO2014/102668A2で論じられた、(擬似)線形モデルを使用するなどの代替方法の手段によって推測することができ、或いは、本明細書で例示的に使用された標的遺伝子のリストを、推測経路活性に関して比較可能な結果を得ることが証明されたその証拠となる性質に基づいて、証拠精選リストから標的遺伝子を別に選択したものと置き換えることもできる。代替リストは、本明細書ではPI3K経路(表1及び表2参照)について論じられ、WO2013/011479A2では、Wnt経路(WO2013/011479A2の表6参照)、ER経路(WO2013/011479A2の表7参照)、及びHH経路(WO2013/011479A2の表8参照)について論じられている。
本明細書で、我々は、コックスの比例ハザード・モデルを使用して重み付け係数w、w、w、及びwの適切な値を推測する好ましい方法について説明する。コックスの比例ハザード・モデルは、適切な数(好ましくは癌種類の多様なサブセットを表す、好ましくは>100)の、推測活性がP、P、P及びPである試料、及びサバイバル・データすなわち生存時間から成る訓練セットと、検閲情報とに、たとえばMATLAB(MATLAB R2014a、The MathWorks Inc.,Natick,MA)又はR(v3.0.3,R Core Team(2014)。Rは統計計算のための言語及び環境。R Foundation for Statistical Computing,Vienna,Austria)を使用して適合される。例示的に、公的に得られる乳癌の、Guyの病院に由来するGSE6532からの試料(n=87)、及びhttp://www.ncbi.nlm.nih.gov/geo/、最終アクセス2014年7月20日、でアクセス可能なGSE9195からの試料(n=77)が、訓練データセットとして使用された。コックスの比例ハザード回帰モデルが、すべての経路の活性に適合され、その結果、経路活性ごとのコックスの係数、その関連する係数推定値の標準誤り(SE)、コックスの係数の指数であるハザード比(HR)、ハザード比の95%信頼区間、及びコックスの係数及び標準誤りから導出されたp値が、表14に見ることができるように得られた。係数推定値の符号は、経路活性が、負係数の場合に臨床事象に対して防御的であるかどうかを示し、又は正係数の場合に悪い予後を予測する。係数の絶対値は、予後に関するリスク・スコアの強さを示す。
Figure 0007065609000018
本発明者らによって、訓練データ・セットについてそれぞれの細胞シグナル伝達経路の活性に適合されたコックスの係数が、たとえば表14に示されるように、リスク・スコアの線形重み付け係数として使用するのに良い値であることが見出された。したがって、これらのコックスの係数は、式(5)の重み付け係数として使用されることが好ましい。リスク・スコアの決定に使用するためのその適性が、以下に説明するように、非常に詳細に評価された。
まず、PI3K経路の活性が、Wnt経路、ER経路、及びHH経路それぞれの活性と組み合わされて、以下の式が得られた。
MPSpw=0.80(±0.41)・P+1.30(±0.85)・P (6)
MPSpe=0.80(±0.41)・P+(-1.02(±0.52))・P (7)
MPSph=0.80(±0.41)・Pp+0.83(±0.54)・P (8)
次に、PI3K経路の活性が、Wnt経路、ER経路、及びHH経路から成る群の異なる2つの経路の活性と組み合わされて、以下の式が得られた。
MPSpwe=0.80(±0.41)・P+1.30(±0.85)・P+(-1.02(±0.52))・P (9)
MPSpwh=0.80(±0.41)・P+1.30(±0.85)・P+0.83(±0.54)・P (10)
MPSpeh=0.80(±0.41)・P+(-1.02(±0.52))・P+0.83(±0.54)・P (11)
コックスの係数が、式(5)に列記されたPI3K経路、Wnt経路、ER経路、及びHH経路の活性の一次結合をパラメータ化するのに用いられることが特に好ましく、それにより以下の式が得られる。
MPSpweh=0.80(±0.41)・P+1.30(±0.85)・Pw+(-1.02(±0.52))・P+0.83(±0.54)・P (12)
ここで、係数の標準誤りが括弧間に列記されている。
別法として、(擬似)線形モデルを、本明細書及びより詳細にWO2014/102668A2(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)で説明されているように使用して経路活性を推測することができ、またこれらの推測活性を、上記で論じられたのと同様に、確率的なモデルによって推測された経路活性と共に使用することができる。経路活性のこれらの線形モデルを式(6)から(12)に挿入すると、最終的に、加算の展開後に、たった1回の加算を伴う式に一般化できる一次結合になる。
MPSprobesets=Σwij・Eij (13)
ここで、Σは全j個の経路のうちの全i個のプローブセットの合計であり、ここで経路は例示的にPI3K、Wnt、ER及びHHであり、wijはプローブセットに関連した重みであり、これは、「単層」線形モデルついては経路に関連した重みとプローブセット又は経路との積に等しく、「2層」線形モデルについては標的遺伝子とプローブセットとの積に等しい。本明細書では、重みwijは例示的に、j番目の経路のi番目のプローブセットの、訓練データ・セットから推定されたコックスの係数に等しく選択され、Eijはj番目の経路のi番目のプローブセットである。当業者であれば、この式を(RT-q)PCR、配列決定、mRNA fish法、並びに、本明細書で例示的に使用されたアフィメトリクスHG-U133Plus2.0に由来するプローブセットの代わりに標的遺伝子の発現レベルを検出する他の適切な方法など、他の測定プラットフォームに適合させることができよう。
次に、本明細書に記載のリスク・スコアが他の3つのデータセットの組合せについて検査された。GSE20685及びGSE21653は、http://www.ncbi.nlm.nih.gov/geo/、最終アクセス2014年7月20日、でアクセス可能な遺伝子発現オムニバスにおいて得られるのに対し、E-MTAB-365は、http://www.ebi.ac.uk/arrayexpress/ experiments/、最終アクセス2014年7月20日、でアクセス可能なArrayExpressにおいて得られる。3つのデータセットでは、合計1005人の乳癌患者の様々な組を、完全な生存時間及び検閲データと組み合わせる。これらの患者のリスク・スコアが式(6)から(13)に従って計算され、次に、リスク・スコアの予後値が、このような予後値を量子化する2つの方法を使用して調査された。これらは、コックスの比例ハザード回帰モデル、及び対数順位統計検定と合わせたカプラン-マイヤー・グラフである。
第1の方法は、コックスの比例ハザード・モデルを、1つ又は複数の共変量によって生存データに適合させる。手短に言えば、このようなハザード・モデルは、共変量の(数)値に基づいて母集団内の生存(臨床事象)に関する変動を説明する。適合の結果として、含まれる各共変量に、コックスの係数の指数であるハザード比(HR)が割り当てられ、このハザード比は、臨床事象の関連リスクを共変量の値に基づいて定量化し、たとえば2のHRは、共変量の値に1の増加がある患者に対して、対象となる臨床事象の2倍高いリスクに相当する。詳細には、HR=1という値は、この共変量が生存に影響を及ぼさないことを意味するのに対し、HR<1では、共変量数の増加がより低いリスクを表し、共変量数の減少がより高いリスクを表し、またHR>1では、共変量数の増加がより高いリスクを表し、共変量数の減少がより低いリスクを表す。ハザード比と共に、95%信頼区間及びp値が報告される(すなわち、ハザード比が1よりも著しく小さい、又は大きいという片側確率)。全リスク・スコアが、ハザード比の直接比較を簡単明瞭にするために、リスク・スコアのスケール(最小値から最大値まで)が1になるように基準化される。
後者の方法には、臨床事象の生存確率を時間の関数として表すカプラン-マイヤー曲線のグラフ化が含まれる。たとえば、母集団内の異なるリスク群について、例示的な予後検査に基づいてカプラン-マイヤー曲線をグラフ化することによって、例示的な臨床事象のリスク分別の特性を可視化することができる。すなわち、より分岐するリスク群は、リスク・スコアが、リスクのある患者の層別化ではより良いことを示す。この特性は、2つの生存曲線が、完全な追跡調査期間を考慮に入れて等しくなる確率(p値)を計算する対数順位検定によって、さらに定量化することができる。
PI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上との推測活性を少なくとも使用するリスク・スコアの結果は、本明細書で提示されるように、個々の推測活性P、P、P、及びP、すなわち、PI3K経路、Wnt経路、ER経路、及びHH経路それぞれの推測活性と比較して、また、本明細書で提示されるように、P、P、Pの非一次結合、及びGenomic Healthの乳癌Oncotype DX(登録商標)検定と比較して、評価された。P、P、及びPの非一次結合は、以下のように計算される。
MPSewh=-P+max(P,P) (14)
MPSewhは、乳癌患者における再発の適切な予知因子であることが示された。これは式(14)を用いて計算され、患者は、低リスク、中リスク及び高リスクの患者に、式に示されたMPSewhの閾値を使用して、すなわちそれぞれ-0.1及び0.1において、層別化された。Oncotype DX(登録商標)検定は、ER陽性乳癌患者に関する再発の適切な予知因子であることが示された。Oncotype DX(登録商標)検定は、ここではハザード比の直接比較のために0と1の間に基準化される0と100の間のリスク・スコア又は再発スコア(RS)を返し、このRSは、遺伝子のパネルについて測定された発現レベルの組合せに基づいて計算される。(S.Paikらの「A multi-gene assay to predict recurrence of Tamoxifen-treated,node-negative breast cancer」、The New England Journal of Medicine,Vol.351,No.27,2004、2817~2826頁、C.Fanらの「Concordance among gene-expression-based predictors for breast cancer」、The New England Journal of Medicine,Vol.355,No.6,2006,560~569頁参照)。RSは、ER陽性、HER2陰性(タンパク質染色又はFISH)、ノード陰性乳癌患者の10年生存に関して最適化される。RSは、上記のデータセットに報告されたマイクロアレイ発現データを使用して、Fanらにより報告された手順に従って計算され(C.Fanら(2006)参照)、患者は続いて、低リスク患者、中リスク患者、及び高リスク患者に、Oncotype DX(登録商標)リスク層別化アルゴリズム(Paikら(2004)参照)に従って分けられた。
まず、コックスの比例ハザード回帰が、基準化リスク・スコアについてE-MTAB365、GSE20685及びGSE21653の乳癌患者を用いて行われた。計算された一変量のコックスの係数、その標準誤り、ハザード比、関連する95%信頼区間、及びp値が表15に示されている。印象的なことに、PI3K経路の活性を他の細胞シグナル伝達経路のうちの1つの活性と組み合わせるすべてのリスク・スコアが、コックスの係数の高い絶対値で表されるように、個々の経路活性よりも適切に機能し、これは、他の細胞シグナル伝達経路の活性と一緒にPI3K経路の活性を組み合わせたものが、臨床事象の予後、この場合には無病生存、に関して個々の経路の活性よりも適切に機能したことを示す。加えて、2つの細胞シグナル伝達経路の組合せ活性のp値もまた、この優位性を示す。その理由は、これらのp値が一般に、PI3K経路の活性を別の細胞シグナル伝達経路の活性と組み合わせたものについては、個々の経路活性のp値よりも小さいからである。PI3K経路の活性を他の2つの細胞シグナル伝達経路の活性と組み合わせることによりまた、コックスの係数(及びp値)は、2つの経路活性に基づくリスク・スコアと比較して改善した。PI3K経路、Wnt経路、ER経路、及びHH経路の活性を、式(12)に示されるように組み合わせるMPSpwehリスク・スコアは、最適には、個々の経路活性よりも、並びにPI3K経路の活性を2つ又は3つの細胞シグナル伝達経路の活性と組み合わせたものよりも、係数、標準誤り、HR及びp値に関して明白であるように、適切に機能する。加えて、MPSpwehスコアに使用されるものと同じプローブセットを含むMPSprobesetsリスク・スコアは、PI3K経路の活性、及び他の1つの細胞シグナル伝達経路の活性を含むリスク・スコアよりも、コックスの回帰結果から明らかなように、優れている。それにもかかわらず、MPSprobesetsの性能は、MPSpeh及びMPSpwehよりも悪く、これは、多数の適合係数により(MPSpeh及びMPSpwehそれぞれで3個及び4個の係数に対し、MPSprobesetsでは276個の係数)、リスク・スコアを訓練データに「過度に適合させること」の結果である可能性がある。PI3K経路の活性と他の1つ又は複数の経路の活性とを組み合わせるすべてのリスク・スコアが、MPSewh及びRSリスク・スコアよりも、それぞれのコックスの係数から明らかなように、適切に機能した。
Figure 0007065609000019
次に、対象のリスク・スコアの予後層別化が、対数順位検定と併せてカプラン-マイヤー・グラフを使用して解析された。本明細書に記載の新しいリスク・スコアのための単純化アルゴリズムを例示的に使用して、患者をそのリスク・スコアに応じて層別化した。1005人の患者が、増加するリスク・スコアからなる等しいサイズの3つの群(n=335)に分けられる。すなわち、カットオフが全患者のそれぞれのリスク・スコアの三分位値のところにある。リスク層別化の、前述の方法に対する他の変形形態は、当業者によって理解され、既知の最適化技法を使用して行われ得る。たとえば、ヨーデンのJ統計を使用してリスク閾値を推測することができる。比較のために含められた他のリスク・スコアのリスク層別化は、その発明者によって記述されたように行われる。すなわち、PI3K経路、Wnt経路、ER経路、及びHH経路の活性を使用して、それぞれの経路が活性であるか、すなわち0から1のスケールで0.5を超える活性であるか、それとも不活性であるか、すなわち0から1のスケールで0.5以下であるかどうかに応じて患者を層別化する。MPSewhが-0.1以下の患者は低リスクであると考えられ、MPSewhが0.1以上の患者は高リスクであると考えられ、残りのすべての患者(MPSewhが-0.1と0.1の間)は中リスクであると考えられる。一方で、RSが18未満の患者は低リスクであると考えられ、RSが31以上の患者は高リスクであると考えられ、残りのすべての患者(RSが18から31の間)は中リスクであると考えられる。本明細書に記載の新しいリスク・スコアのカプラン-メイヤー・グラフ、すなわち、MPSpw(図2参照)、MPSpe(図3参照)、MPSph(図4参照)、MPSpwe(図5参照)、MPSpwh(図6参照)、MPSpeh(図7参照)、MPSpweh(図8参照)及びMPSprobesets(図9参照)が、図2から図9に提示されている。これらのグラフにおいて、垂直軸は、患者群に対する無再発生存を示し、水平軸が年単位の時間を示す。低、中、及び高リスク群(それぞれ335人の患者を含む)は、それぞれ実線(特徴的に上の線)、点線(特徴的に中間の線)、及び破線(特徴的に下の線)で示されている。これらのグラフは、対象が臨床事象を特定の期間内に経験する可能性のあるリスクの、異なる群間の明確な区別を示す。リスク層別化におけるこの差異は、対数順位検定によって量子化することができる。ここでは、最も高いリスク群対最も低いリスク群(個々の経路活性の場合には、これは活性対不活性になる)のカプラン-マイヤー曲線を比較することが選択された。対数順位p値は、表15の最後の欄に表されている。カプラン-マイヤー・グラフ、及び関連する対数順位統計はさらに、PI3K経路の活性と別の1つの細胞シグナル伝達経路の活性とを含むリスク・スコアの利点を例証するものである。その理由は、このリスク・スコアを使用して疾患再発の低リスク又は高リスクの患者を層別化できるからである。
図10は、5年(実線)及び10年(点線)の無病生存期間の可能性を、例として非密封MPSpwehを使用して示す。区分的曲線は、-0.5を超える値に対して疾患再発の可能性/リスクの強度の(単調)増加を示し、この値未満では、リスクは横ばいになるように見える。したがって、カットオフをこの値の近くに置くことは理に適っている。さらに、ユーザが使いやすいように、複数経路スコアは、負の値を含む範囲をカバーする代わりに再基準化して、ゼロから始めてある特定の正の数まで、たとえば0から15、又は0から100の間に広がることもできる。たとえば、これらの閾値を含む再基準化MPSpwehは、以下のように見え得る。
Figure 0007065609000020
GSE6532及びGSE9195の乳癌患者の初期訓練セットについて訓練されたMPSpw、MPSpe、MPSph、MPSpweh及びMPSprobesetsリスク・スコアは、乳癌試料の他のデータセットを適切に一般化することが示された。別法として、リスク・スコアは、以前に記載したデータセット、すなわち、GSE6532、GSE9195、E-MTAB-365、GSE20685及びGSE21653について同時に(生存データがある合計1169人の患者)、以前に論じた推定コックスの係数を用いて訓練することもできる。これにより、以下のリスク・スコアが得られる。
MPSpw=0.70(±0.17)・P+0.38(±0.26)・P (16)
MPSpe=0.70(±0.17)・P+(-0.87(±0.18))・P (17)
MPSph=0.70(±0.17)・P+0.90(±0.20)・P (18)
MPSpwe=0.70(±0.17)・P+0.38(±0.26)・P+(-0.87(±0.18))P (19)
MPSpwh=0.70(±0.17)・P+0.38(±0.26)・P+0.90(±0.20)・P (20)
MPSpeh=0.70(±0.17)・P+(-0.87(±0.18))・P+0.90(±0.20)・P (21)
MPSpweh=0.70(±0.17)・P+0.38(±0.26)・P+(-0.87(±0.18))・P+0.90(±0.20)・P (22)
別法として、リスク・スコアの係数は、データセットついて別々に推定されたコックスの係数を組み合わせることによって決定することもできる。別々に決定されたコックスの係数は、その標準誤りと共に、最大の可能性のある推定を用いて各経路の活性の真の係数を推定するのに使用される。Guyの病院からの、患者の両データセット、GSE6532及びGSE9195は、そのサイズが小さいので一緒にして1つの訓練データセットにされた。最も可能性のある係数の値は、別々に決定された係数推定値を、データセットに含まれる患者数を係数推定の標準誤りで割ったもので重み付けすることによって決定された。
Figure 0007065609000021
ここで、nはデータセットiに含まれる患者の数であり、b^は真の係数値の推定量であり、bはデータセットiのコックスの係数であり、σはデータセットiから推定されたコックスの係数の標準誤りである。この最小化は、PI3K経路、Wnt経路、ER経路、及びHH経路の活性それぞれについて行われた。真の係数推定値の分散は、フィッシャー情報行列を使用して決定された。これらの値を使用して前述の経路活性の一次結合をパラメータ化することにより、以下のリスク・スコアが得られる。
MPSpw0.65(±0.09)・P+0.12(±0.09)・P (24)
MPSpe=0.65(±0.09)・P+(-0.85(±0.06))・P (25)
MPSph=0.65(±0.09)・P+0.68(±0.16)・P (26)
MPSpwe=0.65(±0.09)・P+0.12(±0.09)・P+(-0.85(±0.06))・P (27)
MPSpwh=0.65(±0.09)・P+0.12(±0.09)・P+0.68(±0.16)・P (28)
MPSpwh=0.65(±0.09)・P+(-0.85(±0.06))・Pe+0.68(±0.16)・P (29)
MPSpweh=0.65(±0.09)・P+0.12(±0.09)・P+(-0.85(±0.06))・P+0.68(±0.16)・P (30)
実施例3 CDSアプリケーション
図11(本明細書に開示されているように、対象が臨床事象を特定の期間内に経験するリスクを示すリスク・スコアを決定するように構成された臨床判断支援(CDS)システムを図示する)を参照すると、臨床判断支援(CDS)システム10が、適切に構成されたコンピュータ12として実現される。コンピュータ12は、ハード・ドライブ若しくは他の磁気記憶媒体、光ディスク若しくは他の光記憶媒体、ランダム・アクセス・メモリ(RAM)、読み出し専用メモリ(ROM)、フラッシュ・メモリ、又は他の電子記憶媒体、ネットワーク・サーバなどの非一時的記憶媒体(図示せず)に記憶されている適切なソフトウェア、ファームウェア、又は他の命令を実行することによってCDSシステム10として動作するように構成することができる。説明的なCDSシステム10は、説明的なコンピュータ12によって具現化されているが、より一般的にCDSシステムは、本明細書で論述されている臨床判断支援方法を実施するように構成されたデジタル・プロセッサを備える、デジタル処理デバイス又は装置によって具現化することができる。たとえば、デジタル処理デバイスは、手持ち型デバイス(たとえば、CDSアプリケーションを実行する携帯情報端末又はスマートフォン)、ノート型コンピュータ、デスクトップ・コンピュータ、タブレット・コンピュータ若しくはデバイス、リモート・ネットワーク・サーバなどであり得る。コンピュータ12又は他のデジタル処理デバイスは通常、表示装置14を含むか、又はそれに動作可能に接続され、この表示装置を介して臨床判断支援推奨を含む情報が医療関係者に表示される。コンピュータ12又は他のデジタル処理デバイスはまた通常、説明的なキーボード16、若しくはマウス、トラックボール、トラックパッド、タッチ・センサ式画面(場合により表示装置14と一体化)、又は他のポインタ・ベースのユーザ入力デバイスなどの、1つ又は複数のユーザ入力デバイスを含むか、又はそれに動作可能に接続され、このユーザ入力デバイスを介して医療関係者が、CDSシステム10を制御するための操作コマンド、CDSシステム10で使用するためのデータなどの情報を入力することができる。
CDSシステム10は入力として、対象(たとえば、癌専門医、外科医若しくは他の医療関係者の治療を受けている病院患者若しくは外来患者、或いは結腸癌、乳癌若しくは肝臓癌などの特定の種類の癌があることが分かっている、又は疑われている、癌のスクリーニング若しくは他の診察を受けている人)に関する情報を受け取る。CDSシステム10は、表示装置14を介して(又は、人間が知覚できる出力が得られる音声合成装置若しくは他のデバイスを介して)医療関係者に提示される臨床判断支援推奨を生成するために、この入力情報に様々なデータ解析アルゴリズムを適用する。いくつかの実施形態では、これらのアルゴリズムは、臨床指針を患者に適用することを含み得る。臨床指針は、通常では医療専門家委員会の推奨に基づいて構築され、任意選択で、臨床指針に目を通しやすくするための臨床「流れ図」の形にフォーマットされる、標準又は「基準」の治療推奨の記憶セットである。様々な実施形態では、CDS10のデータ処理アルゴリズムは追加として、又は別法として、様々な診断アルゴリズム又は臨床試験アルゴリズムを含むことができ、これらのアルゴリズムは、本明細書で開示されている機械学習法などの臨床判断推奨を抽出するために、入力情報に対して実施される。
本明細書に開示されている説明的なCDSシステム(たとえば、CDSシステム10)では、CDSデータ解析アルゴリズムは、1つ又は複数の診断アルゴリズム又は臨床試験アルゴリズムを含み、これは、1つ又は複数の医療検査所18で取得された入力ゲノム及び/又はプロテオミクスの情報に対して実施される。これらの検査所は「現地に」、すなわち、対象が医療検査及び/又は治療を受けている病院又は他の場所に、又は「現地外に」様々に設置することができる、たとえば、対象の試料を(郵便又は他の配達サービスによって)受け取る専門化及び集中化した検査所であり、試料は、対象から取り出されている(たとえば、癌病変部から、又は癌の疑いがある病変部から、又は転移腫瘍から、又は癌細胞で汚染されている流体が存在する体腔(たとえば、胸膜若しくは腹部の腔、又は膀胱腔)から、又は癌細胞を含有する他の体液などから、好ましくは組織診手順又は他の試料抽出手順を介して得られた試料)。試料が取り出された細胞はまた、悪性血液疾患(白血病又はリンパ腫など)からの腫瘍細胞でもあり得る。場合によっては、細胞試料はまた、循環腫瘍細胞、すなわち血流に入った腫瘍細胞であることもあり、適切な分離技法、たとえばアフェレーシス又は従来の静脈血採血を用いて取り出すことができる。血液は別として、試料が抽出される体液は、尿、胃腸内容物、又は血管外湧出物であり得る。
試料は検査所で処理されて、ゲノム又はプロテオミクスの情報が生成される。たとえば、試料は、マイクロアレイ(当技術分野ではまた様々に、遺伝子チップ、DNAチップ、バイオ・チップなどとも呼ばれている)を使用して、又は定量的ポリメラーゼ連鎖反応(qPCR)処理によって、処理して、対象の遺伝子の発現レベルなどの証拠となるゲノム又はプロテオミクスの情報を、たとえば、遺伝子から転写されるメッセンジャーリボ核酸(mRNA)のレベル、又は遺伝子から転写されたmRNAから翻訳されるたんぱく質のレベルの形で測定する。別の例として、試料は、遺伝子配列決定検査所で処理して、デオキシリボ核酸(DNA)の配列を生成することも、RNA配列、コピー数変化、メチル化などを生成することもできる。他の企図される測定手法には、病理スライド上で行われる免疫組織化学的検査(IHC)、細胞診断、蛍光in situハイブリダイゼーション(FISH)、近接ライゲーション・アッセイなどが含まれる。マイクロアレイ処理、質量分析、遺伝子配列決定、又は他の検査所技法によって生成できる他の情報には、メチル化情報が含まれる。このようなゲノム及び/又はプロテオミクスの測定を様々に組み合わせることもまた実施することができる。
いくつかの実施形態では、医療検査所18は、対象の試料に対していくつかの標準化データ取得を実施して、大量のゲノム及び/又はプロテオミクスのデータを生成する。たとえば、標準化データ取得技法により、1つ又は複数の染色体又は染色体部分の、又はゲノム全体の、(任意選択で配列された)DNA配列を生成することができる。標準マイクロアレイを適用すると、多数の遺伝子の発現レベル、様々なメチル化データなど、数千又は数万のデータ項目を生成することができる。同様に、PCRによる測定を用いて、遺伝子のうちの選択されたものの発現レベルを測定することができる。この大量のゲノム及び/又はプロテオミクスのデータ、又はそのうちの選択された部分がCDSシステム10に入力されて処理され、それによって、臨床判断支援推奨を策定するための臨床的に有用な情報が創出される。
開示されたCDSシステム及び関連する方法は、ゲノム及び/又はプロテオミクスのデータを処理して様々な細胞シグナル伝達経路の活性を評価すること、並びに対象が臨床事象(たとえば、癌)を特定の期間内に経験するリスクを示すリスク・スコアを決定することに関する。しかし、開示されたCDSシステム(たとえば、CDSシステム10)は任意選択で、バイタルサイン監視データ、患者病歴データ、患者人口統計データ(たとえば、性別、年齢など)、患者医療画像データなどの様々な患者データに基づく記憶された臨床指針に従って、臨床判断支援推奨を生成することなど、多様な付加機能をさらに含むことができることを理解されたい。別法として、いくつかの実施形態では、CDSシステム10の機能は、本明細書に開示のように、ゲノム及び/又はプロテオミクスのデータ解析だけを行って細胞シグナル伝達経路の活性を評価すること、及び対象が臨床事象(たとえば、癌)を特定の期間内に経験するリスクを示すリスク・スコアを決定することに限定することもできる。
例示的な図11を引き続き参照すると、CDSシステム10は、対象における1つ又は複数の細胞シグナル伝達経路の、ここではPI3K経路、並びにWnt経路、ER経路、及びHH経路のうちの1つ以上の、活性22を、それだけには限らないが、対象の試料において測定された細胞シグナル伝達経路の1つ又は複数の標的遺伝子の発現レベル20に基づいて、推測する。PI3K経路、Wnt経路、ER経路、及びHH経路は、これらの経路の調節が失われることが癌の増殖の原因になり得るので、腫瘍学の様々な分野で注目されている。約10個から15個の関連するシグナル伝達経路があり、各癌は、調節解除されている少なくとも1つの主要経路によって駆動される。何か特定の動作原理に制限されなければ、これらの経路は細胞増殖を調節し、その結果、癌細胞におけるこれらの経路の調節が失われることにより経路が「常にオン」になり、それによって癌細胞の増殖が加速される可能性があり、これがひいては、癌の成長、浸潤又は転移(拡散)として顕在化する。
細胞シグナル伝達経路を形成するタンパク質カスケードの一部である中間タンパク質などの、細胞シグナル伝達経路の調節タンパク質をコードする遺伝子のmRNA発現レベルの測定は、調節タンパク質発現レベルの間接的測定であり、実際の調節タンパク質発現レベルとの強い(細胞シグナル伝達経路の活性全体とはずっと弱い)相関関係があることもないこともある。細胞シグナル伝達経路は、標的遺伝子の転写を直接調節し、したがって、標的遺伝子から転写されたmRNAの発現レベルは、この調節活性の直接の結果になる。それゆえに、CDSシステム10は、1つ又は複数の細胞シグナル伝達経路(ここでは、PI3K経路、並びにWnt経路、ER経路、及びHH経路のうちの1つ以上)の活性を、細胞シグナル伝達経路の1つ又は複数の標的遺伝子の発現レベル(代用測定としてmRNA又はタンパク質レベル)に基づいて推測する。これにより、CDSシステム10が経路の活性を、測定された標的遺伝子の発現レベルによって得られた直接の情報に基づいて推測することが確実になる。
推測活性、この例ではP、P、P及びP、すなわちPI3K経路、Wnt経路、ER経路、及びHH経路の推測活性を次に使用して、本明細書に詳細に記載のように、対象が臨床事象、この例では癌、特に乳癌を特定の期間内に経験するリスクを示すリスク・スコアを24で決定する。リスク・スコアは、推測活性を組み合わせたものに基づく。たとえば、リスク・スコアは、式(4)又は(5)に関して詳細に説明されたように計算された「複数経路スコア」(MPS)とすることができる。
決定されたMPSに基づいて、この例ではCDSシステム10は、26で対象を、対象が臨床事象を特定の期間内に経験することの別々に示されたリスクと関連付けられた複数のリスク群のうちの少なくとも1つに割り当て、かつ/又は28で、対象に推奨される治療を、対象が臨床事象を特定の期間内に経験することの示されたリスクに基づいて決定する。
特定の患者についてのMPS及び/又はリスク分類をCDSシステムによって決定すること、又は本明細書に記載のMPS及びリスク分類の独立した実施により、患者の診断又は治療又は監視/経過観察に関わる腫瘍医、外科医、又は他の医療関係者が治療を、患者が長期生存の最善の機会を得られるように、同時に、望ましくない副作用、特に積極的な化学療法及び/又は標的治療及び/又は免疫療法及び/又は放射線療法及び/又は手術の副作用が最小限になるように、適応させることが可能になる。したがって、たとえば、癌再発のリスクが低い患者、すなわちMPSが低い患者、及び/又は本明細書に記載のリスク層別化アルゴリズムに基づいて低リスクとして分類された患者は現在、通常はホルモン療法だけで、又はホルモン療法を組み合わせたもの、たとえば抗エストロゲン剤及び/又はアロマターゼ阻害薬と低毒性化学療法剤とで治療される。一方で、癌再発のリスクが中間又は高い患者、すなわちMPSが中から高の患者、及び/又は本明細書に記載のリスク層別化アルゴリズムに基づいて中リスク又は高リスクとして分類された患者は現在、通常はアントラサイクリン及び/又はタキサンがベースの治療養生法などの、より積極的な化学療法で治療されることになる。加えて、MPSは、場合により他の患者検査結果及び/又は他の予後又は予測的な(たとえば、コンパニオン診断)検査の結果と組み合わせて、患者を、タモキシフェン、トラスツズマブ、ベバシズマブなどの標的薬剤で、及び/又は、現在は患者の特定の癌に対する主系統治療プロトコルの一部ではない他の治療薬剤(たとえば、免疫療法)で、及び/又は放射線療法、たとえば密封小線源療法で、及び/又は、たとえば一次治療の前及び/又は後の、異なる治療のタイミングなどの他の治療選択肢で、治療する決定を行うことができる。
対象が臨床事象(たとえば、癌)を特定の期間内に経験するリスクを示すものとして決定リスク・スコア(MPS)を直接使用する代わりに、CDSシステム10が、リスク・スコアを1つ又は複数の追加予後検査から得られた1つ又は複数の追加リスク・スコアと組み合わせて結合リスク・スコアが得られるように構成されることが可能であり、この結合リスク・スコアが、対象が臨床事象を特定の期間内に経験するリスクを示すことに留意されたい。1つ又は複数の追加予後検査には特に、Oncotype DX(登録商標)乳癌検査、Mammostrat(登録商標)乳癌検査、MammaPrint(登録商標)乳癌検査、EndoPredict(登録商標)乳癌検査、BluePrint(商標)乳癌検査、CompanDx(登録商標)乳癌検査、Breast Cancer Index(役務標章)(HOXB13/IL17BR)、OncotypeDX(登録商標)結腸癌検査、及び/又は、遺伝子/タンパク質Ki67の発現を測定することによって行われる増殖検査が含まれ得る。
実施例4 本発明を説明するためのさらなる情報
(1)遺伝子発現のレベルの測定
本発明に記載の標的遺伝子の固有の組から導出されたデータが、本明細書に記載の方法を用いて細胞シグナル伝達経路の活性を推測するために、さらに利用される。
取り出された試料における遺伝子発現レベルを解析する方法は、広く知られている。たとえば、ノーザン・ブロット法、PCRの使用、ネステッドPCR、定量的リアルタイムPCR(qPCR)、RNA-seq、又はマイクロアレイなどの方法はすべて、遺伝子発現レベル・データを導出するのに使用することができる。当技術分野で知られている、標的遺伝子の遺伝子発現を解析するためのすべての方法が本明細書では企図されている。
遺伝子の発現産物をPCRベースの方法を用いて決定する方法が特に使用され得る。遺伝子発現のレベルをPCRの使用により定量化するために、対象の各PCR産物の量が通常は、従来の定量的リアルタイムPCR(qPCR)を使用して推定されて、PCR産物の蓄積が、増幅の各サイクル後にリアルタイムで測定される。これには通常、挿入色素、副溝結合色素、又は蛍光発生プローブなどの検出可能レポーターを利用し、それによって、光を当てるとレポーターが励起されて蛍光を発生させ、発生させた蛍光が通常は、参照により本明細書に組み込まれる米国特許第6,713,297号に開示されているものなどの、CCDカメラ又は光増倍管検出システムを使用して検出される。
いくつかの実施形態では、定量的リアルタイムPCR(qPCR)アッセイにおいてPCR産物の検出に使用されるプローブは、蛍光マーカーを含み得る。多数の蛍光マーカーが市販されている。たとえば、Molecular Probes,Inc.(Eugene,Oreg.)は、多種多様な蛍光色素を販売している。非限定的な例には、Cy5、Cy3、TAMRA、R6G、R110、ROX、JOE、FAM、Texas Red(商標)、及びOregon Green(商標)が含まれる。付加的な蛍光マーカーには、IDT ZEN Double-Quenched Probeが、qPCRアッセイにおける従来の5’加水分解プローブと共に含まれ得る。これらのプローブは、たとえば、5’FAM色素を3’TAMRA Quencher、3’Black Hole Quencher(BHQ,Biosearch Technologies)、又は内部ZEN Quencher及び3’Iowa Black Fluorescent Quencher(IBFQ)と共に含有し得る。
本発明によると有用な蛍光色素は、当技術分野でよく知られている方法を用いて、オリゴヌクレオチド・プライマーに付着させることができる。たとえば、蛍光標識をオリゴヌクレオチドに付加するための1つの一般的な方法は、色素のN-ヒドロキシスクシンイミドエステル(NHS)エステルを標的上の反応性アミノ基と反応させるものである。ヌクレオチドは、たとえば核酸塩基上にアリルアミン基を含むことによって、反応性アミノ基を保有するように修飾することができる。アリルアミンによる標識は、たとえば、参照により本明細書に組み込まれる米国特許第5,476,928号及び第5,958,691号に記載されている。蛍光でヌクレオチド、オリゴヌクレオチド及びポリヌクレオチドを標識する他の手段は、当業者によく知られている。
他の蛍光発生の手法には、SYBR-green色素などの遺伝子検出系の使用が含まれ、これは、参照により本明細書に組み込まれる米国特許第5,436,134号及び第5,658,751号に開示されているように、任意の遺伝子発現産物からの増幅DNAとインターカレートされると蛍光を発生させる。
標的遺伝子発現レベルを決定するための別の有効な方法には、異なる生理学的条件間の遺伝子発現レベル差、又は発生中若しくは疾患進行の間に起こる変化を含めてトランスクリプトーム解析に使用される強力な解析ツールである、RNA-seqが含まれる。
遺伝子発現レベルを決定するための別の手法には、マイクロアレイ、たとえばRNA及びDNAマイクロアレイを使用することが含まれ、これは当技術分野でよく知られている。マイクロアレイは、多数の遺伝子の発現を同時に定量化するのに使用することができる。
(2)PI3K、Wnt、ER、及びHH細胞シグナル伝達の活性を決定するための一般化されたワークフロー
本発明は、特定の臨床事象を経験する対象のリスク・スコアを計算するためにPI3K、Wnt、ER、及びHHの経路の機能状態又は活性を評価する、本明細書に開示された新規で改善された方法及び装置を提供する。
PI3K細胞シグナル伝達及び他の細胞シグナル伝達の活性を対象から取り出された試料により決定するプロセスを例示的に図示する流れ図が、図12に示されている。第一に、試料からmRNAが分離される(11)。第二に、本明細書に記載のように、少なくとも3つ以上のPI3K標的遺伝子からなる固有の組のmRNA発現レベルが、当技術分野で知られている遺伝子発現を測定する方法を使用して測定される(12)。次に、PI3K転写因子(TF)要素の活性レベルが、3つ以上の標的遺伝子の発現レベルをPI3K TF要素の活性レベルに関係付ける較正済み数学経路モデル(14)を使用して決定される(13)。次に、対象におけるPI3K経路の活性が、対象の試料におけるPI3K TF要素の決定された活性レベルに基づいて推測される(15)。
図12の右側に示されるように、PI3K TF要素の活性レベルを決定した後に、少なくとも1つの追加細胞シグナル伝達経路(すなわち、Wnt、ER、及びHHのうちの1つ以上)のTF要素の活性レベルが決定される。一例として、本明細書に記載のように、追加細胞シグナル伝達経路からの3つ以上の標的遺伝子からなる固有の組のmRNA発現レベルが、当技術分野で知られている遺伝子発現を測定する方法を使用して測定される(16)。次に、TF要素の活性レベルが、追加細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルをTF要素の活性レベルに関係付ける較正済み数学経路モデル(14)を使用して決定される(17)。次に、対象における追加細胞シグナル伝達経路の活性が、対象の試料におけるTF要素の決定された活性レベルに基づいて推測される(18)。次に、PI3K及び追加細胞シグナル伝達経路の活性が複数経路スコア(MPS)に変換され(19)、これは、疾患と関連する臨床事象を対象が規定期間内に経験するリスクを示す。MPSとしてリスク・スコアを決定することは、較正済み複数経路スコア(MPS)モデルを評価することとして理解することができ、モデルのパラメータは、たとえば、本明細書に記載のように、重み付け係数(たとえば、w、w、w、及びw)を含む。最後に、試料には、計算されたMPSに基づいた臨床事象を経験するリスク・スコアが割り当てられる(20)。
(3)複数経路スコア(MPS)モデルの較正及び複数経路スコア(MPS)の決定
本明細書で企図されているように、臨床事象が起こるリスクに対応するリスク・スコアは、以下でさらに説明するように、臨床事象と関連する細胞シグナル伝達経路の活性を含む較正済み複数経路スコア(MPS)モデルを使用して決定することができる。
本発明で使用される較正済み複数経路スコア(MPS)モデルは、対象の臨床事象及び推測経路活性についての容易に得られる臨床データを用いて較正することができる。MPSモデルを生存データで較正するためのプロセスを例示的に示す流れ図が、図13に示されている。最初のステップとして、較正済み数学経路モデルを使用して推測された関連経路活性が、経路活性データベース(201)から検索される。経路活性データベースは、PI3K経路活性(205)及び少なくとも1つの追加経路の経路活性を含む。たとえば、経路活性データベースは、ER経路活性(202)、Wnt経路活性(203)、HH経路活性(204)、及びPI3K経路活性(205)を含む。試料の特定の訓練セットのID(218)が次に使用されて、関連経路活性(219)と、たとえば、生存データ・データベース(221)から受け取られる生存データ(220)(生存が、解析される臨床事象である場合)とを受け取る。経路活性が次に、ER経路活性、Wnt経路活性、HH経路活性、及びPI3K経路活性の場合に、それぞれP、P、P、及びPの出力を用いて選択される(222)。生存データは、MPSが使用される所与の期間内の生存時間及び検閲データを反映する変数Surv及びcens(223)に変換される。経路活性及び生存データは次に、コックスの比例ハザード・モデル(224)に適合され、これにより、適合されたコックスの比例ハザード・モデル(225)が得られる。コックスの比例ハザード・モデルからコックスの係数が集められ(226)、次に、出力w、w、w、及びwの重みに割り当てられる(227)。MPS構造(228)と重みが一緒に取得されてMPSモデルを較正し(229)、較正済みMPSモデル(210)が出力される。
較正済みMPSモデルからリスク・スコアを決定するためのプロセスを例示的に示す流れ図が、図14に示されている。最初のステップとして、較正済み経路モデルを使用して推測された関連経路活性が、経路活性データベース(201)から検索される。経路活性データベースは、PI3K経路活性(205)及び少なくとも1つの追加経路の経路活性を含む。たとえば、経路活性データベースは、ER経路活性(202)、Wnt経路活性(203)、HH経路活性(204)、及びPI3K経路活性(205)を含む。患者試料は次に識別され(206)、初期経路活性が試料及びデータベースから、転写因子の測定値又は遺伝子発現レベルとして、関連経路について集められる(207)。関連経路それぞれの総合活性が次に、P、P、P、及びPの出力を用いて推測される(208)。これらの活性は次に、較正済みMPSモデル(210)を使用して、リスク・スコア(209)に変換される。この初期リスク・スコアはさらに、他の関連データを用いて、患者の最終リスク・スコア(211)を生成するようにさらに調整することができる。この最終スコアを使用して次に、表示(212)、割り当て(213)、又は治療についての決定(214)を行い、それによって、それぞれ表示リスク・スコア(215)、割り当てリスク・スコア(216)、又は決定治療(217)の結果が生成される。
対象における細胞シグナル伝達経路の活性を推測することは、たとえば、特に(i)対象の試料において測定された、細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを含む入力の組に対する細胞シグナル伝達経路を表す較正済み確率的経路モデル、好ましくはベイジアン・ネットワーク、の一部分を評価することによって、また(ii)対象における転写因子(TF)要素の活性レベルを推定することによって行うことができ、このTF要素は、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御し、推定することは、TF要素の活性レベルと、対象の試料において測定された細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルとに関連する条件付き確率とに基づいており、また(iii)細胞シグナル伝達経路の活性を対象の試料におけるTF要素の推定活性レベルに基づいて推測することによって行うことができる。これについては、公開されている国際特許出願WO2013/011479A2(「Assessment of cellular signaling pathway activity using probabilistic modeling of target gene expression」)に詳細に記載されており、同出願の内容はその全体が本明細書に組み込まれる。
例示的な一代替形態では、対象における細胞シグナル伝達経路の活性を推測することは、たとえば、特に(i)対象の試料における転写因子(TF)要素の活性レベルを決定することによって行うことができ、このTF要素は、細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御し、決定することは、細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルをTF要素の活性レベルに関係付ける較正済み数学経路モデルを評価することに基づいており、数学経路モデルは、3つ以上の標的遺伝子の発現レベルの1つ以上の一次結合に基づいており、また(ii)対象における細胞シグナル伝達経路の活性を対象の試料におけるTF要素の決定された活性レベルに基づいて推定することによって行うことができる。これについては、公開されている国際特許出願WO2014/102668A2(「Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions」)に詳細に記載されている。
一実施形態では、細胞シグナル伝達経路がWnt経路及び/又はER経路及び/又はHH経路を含み、リスク・スコアが、示されたリスクがWnt経路の推測活性の増加及び/又はHH経路の推測活性の増加と共に単調増加するように、かつ/又はER経路の推測活性の増加と共に単調減少するように規定される方法を提供する。
一実施形態では、リスク・スコアが、示されたリスクがPI3K経路の推測活性の増加と共に単調増加するように規定される方法が提供される。
一実施形態では、推測活性を組み合わせたものは、項w・Pと、項w・P、w・Pe、及びw・Pのうちの1つ以上とを含む合計を含み、ここで、P、P、P、及びPはそれぞれ、PI3K経路、Wnt経路、ER経路、及びHH経路の推測活性を表し、w、w、及びwは正の一定重み付け係数であり、wは負の一定重み付け係数であり、対象が臨床事象を特定の期間内に経験する示されたリスクは、合計値の増加と共に単調増加する。
いくつかの実施形態では、一定重み付け係数w、w、w、及びwは、それぞれの細胞シグナル伝達経路のコックス比例ハザード・モデルを臨床データに適合させることにより得られるコックスの係数の値に基づいて、それぞれ決定されるか、又は決定されている。たとえば、係数推定値の符号は、経路活性が、負係数の場合に臨床事象に対して防御的であるかどうかを示し、又は正係数の場合に良くない、若しくは悪い予後を予測する。係数の絶対値は、予後に関するリスク・スコアの強さを示す。
一実施形態では、臨床事象は癌転移であり、w、w、及びwは非負の一定重み付け係数であり、wは非正の一定重み付け係数である。これらの係数と共にMPSは、対象が臨床事象を特定の期間内に経験する示されたリスクが、合計値の増加と共に単調に増加することを示す。
(4)標的遺伝子発現レベル決定手順
対象から取り出された試料から標的遺伝子発現レベルを導出するためのプロセスを例示的に示す流れ図が、図15に示されている。例示的な実施形態では、試料は検査所で受け取られ登録される。試料は、たとえば、ホルマリン固定パラフィン包埋(FFPE)試料(181)又は新鮮凍結(FF)試料(180)を含み得る。FF試料は、すぐ溶解することができる(183)。FFPE試料では、パラフィンは、プロテイナーゼKを追加するときに加熱インキュベーション・ステップによって除去することができる(182)。次に細胞を溶解し(183)、それにより細胞及び核膜を破壊し、それにより、さらなる処理に利用可能な核酸(NA)が得られるようになる。核酸は、たとえばビーズ又はフィルタであり得る固相に結合される(184)。核酸は次に、洗浄バッファで洗浄されて、溶解後に存在する細胞デブリがすべて除去される(185)。清浄な核酸が次に固相から、溶出バッファを用いて引き離される(186)。DNAがDNAse処理によって取り除かれて、RNAだけが試料中に存在することが確実になる(187)。核酸試料は次に、RT-qPCR試料混合物中ですぐに使用することができる(188)。RT-qPCR試料混合物は、RNA試料、cDNAをRNA試料から調製するためのRT酵素、及びcDNAを増幅するためのPCR酵素、酵素の機能を確保するための緩衝液を含有し、また場合により、固定量の濃度を設定するために分子グレード水を含有し得る。試料混合物は次に、乾燥RT-qPCRアッセイを含むマルチウェル・プレート(すなわち、96ウェル又は384ウェル・プレート)に加えることができる(189)。RT-qPCRを次に、PCR機内で指定プロトコルに従って実行することができる(190)。例示的なPCRプロトコルは、i)50℃で30分、ii)95℃で5分、iii)95℃で15秒、iv)60℃で45秒、v)ステップiii及びivの50サイクル繰り返し、を含む。Cq値が次に、生データによって、第2の誘導法を用いて決定される(191)。Cq値が解析のためにエクスポートされる(192)。
(5)疾患、障害、及び治療方法
本明細書で企図されているように、本発明の方法及び装置は、対象、たとえば疾患若しくは障害があることが疑われる対象、又は疾患若しくは障害がある対象において、PI3K、Wnt、ER、及び/又はHH細胞シグナル伝達経路活性を評価するために利用することができ、シグナル伝達経路のうちの1つの状態は、疾患の存在又は進行の、全体的又は部分的な証拠となるものである。一実施形態では、対象を治療する方法が本明細書で提供され、この方法は、対象から本明細書に記載の方法を用いて分離された試料から導出されたPI3K、Wnt、ER、及び/又はHH細胞シグナル伝達経路の活性状態に関する情報を受け取ること、並びに、細胞シグナル伝達経路の活性に関する情報が、活性のPI3K、Wnt、ER、及び/又はHHシグナル伝達経路を示す場合に、対象にPI3K、Wnt、ER、及び/又はHHの阻害剤を投与することを含む。
本発明において使用され得るPI3K阻害剤はよく知られている。PI3Kの例に含まれるものには、それだけには限らないが、ワートマニン、デメトキシビリジン、ペリフォシン、イデラシブ、ピクチリシブ、パロミド529、ZSTK474、PWT33597、CUDC-907、及びAEZS-136、デュベリシブ、GS-9820、BKM120、GDC-0032(Taselisib)(2-[4-[2-(2-イソプロピル-5-メチル-1,2,4-トリアゾール-3-イル)-5,6-ジヒドロイミダゾ[1,2-d][1,4]ベンゾオキサゼピン-9-イル]ピラゾール-1-イル]-2-メチルプロパンアミド)、MLN-1117((2R)-l-フェノキシ-2ブタニルヒドロゲン(S)-メチルホスホネート;オルメチル(オキソ){[(2R)-l-フェノキシ-2-ブタニル]オキシ}ホスホニウム))、BYL-719((2S)-N1-[4-メチル-5-[2-(2,2,2-トリフルオロ-1、1-ジメチルエチル)-4-ピリジニル]-2-チアゾリル]-1,2-ピロリジンジカルボキサミド)、GSK2126458(2,4-ジフルオロ-N-{2-(メチルオキシ)-5-[4-(4-ピリダジニル)-6-キノリニル]-3-ピリジニル}ベンゼンスルホンアミド)(オミパリシブ)、TGX-221((±)-7-メチル-2-(モルホリン-4-イル)-9-(l-フェニルアミノエチル)-ピリド[l,2-a]-ピリミジン-4-オン)、GSK2636771(2-メチル-1-(2-メチル-3-(トリフルオロメチル)ベンジル)-6-モルホリノ-lH-ベンゾ[d]イミダゾール-4-カルボン酸ジヒドロクロリド)、KIN-193((R)-2-((l-(7-メチル-2-モルホリノ-4-オキソ-4H-ピリド[1,2-a]ピリミジン-9-イル)エチル)アミノ)安息香酸)、TGR-1202/RP5264、GS-9820((S)-l-(4-((2-(2-アミノピリミジン-5-イル)-7-メチル-4-モヒドロキシプロパン-1-オン)、GS-1101(5-フルオロ-3-フェニル-2-([S)]-l-[9H-プリン-6-イルアミノ]-プロピル)-3H-キナゾリン-4-オン)、AMG-319、GSK-2269557、SAR245409(N-(4-(N-(3-((3,5-ジメトキシフェニル)アミノ)キノキザリン-2-イル)スルファモイル)フェニル)-3-メトキシ-4メチルベンズアミド)、BAY80-6946(2-アミノ-N-(7-メトキシ-8-(3-モルホリノプロポキシ)-2,3-ジヒドロイミダゾ[l,2-c]キナズ)、AS252424(5-[l-[5-(4-フルオロ-2-ヒドロキシ-フェニル)-フラン-2-イル]-メタ-(Z)-イリデン]-チアゾリジン-2,4-ジオン)、CZ24832(5-(2-アミノ-8-フルオロ-[l,2,4]トリアゾロ[1,5-a]ピリジン-6-イル)-N-tert-ブチルピリジン-3-スルホンアミド)、Buparlisib(5-[2,6-ジ(4-モルホリニル)-4-ピリミジニル]-4-(トリフルオロメチル)-2-ピリジンアミン)、GDC-0941(2-(lH-インダゾール-4-イル)-6-[[4-(メチルスルホニル)-l-ピペラジニル]メチル]-4-(4-モルホリニル)チエノ[3,2-d]ピリミジン)、GDC-0980((S)-l-(4-((2-(2-アミノピリミジン-5-イル)-7-メチル-4-モルホリノチエノ[3,2-d]ピリミジン-6イル)メチル)ピペラジン-l-イル)-2-ヒドロキシプロパン-l-オン(RG7422としても知られている))、SF1126((8S,14S,17S)-14-(カルボキシメチル)-8-(3-グアニジノプロピル)-17-(ヒドロキシメチル)-3,6,9,12,15-ペンタオキソ-1-(4-(4-オキソ-8-フェニル-4H-クロメン-2-イル)モルホリノ-4-イウム)-2-オキサ-7,10,13,16-テトラアザオクタデカン-18-カルボン酸メチルエステル)、PF-05212384(N-[4-[[4-(ジメチルアミノ)-l-ピペリジニル]カルボニル]フェニル-N’-[4-(4,6-ジ-4-モルホリニル-l,3,5-トリアジン-2-イル)フェニル]尿素)(ゲダトリシブ)、LY3023414、BEZ235(2-メチル-2-{4-[3-メチル-2-オキソ-8-(キノリン-3-イル)-2,3-ジヒドロ-lH-イミダゾ[4,5-c]キノリン-l-イル]フェニル}プロパンニトリル)(ダクトリシブ)、XL-765(N-(3-(N-(3-(3,5-ジメトキシフェニルアミノ)キノキサリン-2-イル)スルファモイル)フェニル)-3-メトキシ-4-メチルベンズアミド)、及びGSK1059615(5-[[4-(4-ピリジニル)-6-キノリニル]メチレン]-2,4-チアゾリデンジオン)、PX886([(3aR,6E,9S,9aR,10R,llaS)-6-[[ビス(プロプ-2-エニル)アミノ]メチリデン]-5-ヒドロキシ-9-(メトキシメチル)-9a,11a-ジメチル-1,4,7-トリオキソ-2,3,3a,9,10,11-ヘキサヒドロインデノ[4,5h]イソクロメン-10-イル]アセテート(ソノリシブとしても知られている))、LY294002、AZD8186、PF-4989216、pilaralisib、GNE-317、PI-3065、PI-103、NU7441(KU-57788)、HS173、VS-5584(SB2343)、CZC24832、TG100-115、A66、YM201636、CAY10505、PIK-75、PIK-93、AS-605240、BGT226(NVP-BGT226)、AZD6482、voxtalisib、アルペリシブ、IC-87114、TGI100713、CH5132799、PKI-402、コパンリシブ(BAY80-6946)、XL 147、PIK-90、PIK-293、PIK-294、3-MA(3-メチルアデニン)、AS-252424、AS-604850、アピトリシブ(GDC-0980、RG7422)、並びにWO2014/071109に記載されている構造がある。或いは、PI3Kの下流のmTOR複合体の阻害剤は、異常なPI3K活性の有用な阻害剤である。mTOR阻害剤の例には、それだけには限らないが、エベロリムス、テムシロリムス、リダフォロリムスが含まれる。或いは、PI3Kの上流のHER2複合体の阻害剤は、異常なPI3K活性の有用な阻害剤である。HER2阻害剤の例には、それだけには限らないが、トラスツズマブ、ラパチニブ、ペルツズマブが含まれる。
内分泌療法をエストロゲン受容体陽性である乳癌に施すことができる。本発明において用いられ得る内分泌療法治療はよく知られている。内分泌療法は、i)通常にはゴナドトロピン放出ホルモン・アゴニスト(GnRHa)を使用して得られる、卵巣機能を抑制するもの、ii)選択的なエストロゲン受容体修飾因子若しくは下方制御因子(SERM又はSARD)、又はiii)アロマターゼ阻害剤(AI)、或いはこれらが組み合わされたものを投与することから成る。卵巣機能を抑制するものには、たとえば、ゴナドトロピン放出ホルモン・アゴニスト(GnRHa)が含まれる。ゴナドトロピン放出ホルモン・アゴニスト(GnRHa)の例には、ブセレリン、デスロレリン、ゴナドレリン、ゴセレリン、ヒストレリン、リュープロレリン、ナファレリン、及びトリプトレリンが含まれ得る。選択的なエストロゲン受容体修飾因子(SERM)には、たとえば、タモキシフェン、トレミフェン、ラロキシフェン、ラソフォキシフェン、バゼドキシフェン、クロミフェン、オルメロキシフェン、オスペミフェン、アフィモキシフェン、及びアルゾキシフェンが含まれる。選択的なエストロゲン受容体下方制御因子(SEDRD)には、たとえば、フルベストラント、SR16234、及びZK191703が含まれる。アロマターゼ阻害剤には、たとえば、アナストロゾール、レトロゾール、ボロゾール、エキセメスタン、アミノグルテチミド、テストラクトン、フォルメスタン、ファドロゾール、アンドロステンジオン、4-ヒドロキシアンドロステンジオン,1,4,6-アンドロスタトリエン-3,17-ジオン、又は4-アンドロステン-3,6,17-トリオンが含まれる。一実施形態では、アロマターゼ阻害剤は、非ステロイド性アロマターゼ阻害剤である。
Wnt阻害剤はよく知られており、それだけには限らないが、ピルビニウム、IWR-l-endo、IWP-2、FH535、WIKI4、IWP-L6、KY02111、LGK-974、Wnt-C59、XAV929、3289-8625、FJ9、NSC 668036、PFK115-584、CGP049090、iCRT3、iCRT5、iCRT14、ICG-001、デメトキシ・クルクミン、CCT036477、KY02111、PNU-74654、又はPRI-724が含まれる。
HH阻害剤はよく知られており、それだけには限らないが、シクロピアミン、SANT1-SANT4、CUR-61414、HhAntag-691、GDC-0449、MK4101、IPI-926、BMS-833923、ロボトニキニン、イトラコナゾール、エリベッジ、オドムゾ、カルシトリオール、コレカルシフェロール、IPI-906、RU-SKI39、又はKAAD-シクロパミン、NVP-LDE225、TAK-441、XL-139、LY2940680、NVP-LEQ506、イトラコナゾール、MRT-10、MRT83、PF-04449913、GANT-61、GANT-58、HPI-1、HPI-3、又はHPI-4が含まれる。
一実施形態では、疾患又は障害は、自己免疫性又は他の免疫性障害、癌、気管支喘息、心臓疾患、糖尿病、遺伝性出血性末梢血管拡張症、マルファン症候群、脈管エーラス・ダンロス症候群、ロイス・ディーツ症候群、パーキンソン病、慢性腎疾患、多発性硬化症、肝線維症、肺線維症又は腎線維症などの線維性疾患、デュピュイトラン病、又はアルツハイマー病のうちの1つである。
特定の実施形態では、対象は癌に罹患している、又はその疑いがあり、その癌は、たとえば、それだけには限らないが、原発性腫瘍又は転位腫瘍、固形腫瘍、たとえば、メラノーマ、肺癌(肺腺癌、基底細胞癌、扁平上皮癌、大細胞癌、細気管支肺胞癌、気管支原性肺癌、非小細胞癌、小細胞癌、中脾腫を含む)と、乳癌(腺管癌、小葉癌、炎症性乳癌、明細胞癌、粘液性癌、漿膜腔乳癌を含む)と、結腸直腸癌(結腸癌、直腸癌、結腸直腸腺癌)と、肛門癌と、膵臓癌(膵臓腺癌、膵島腺癌、神経内分泌腫瘍を含む)と、前立腺癌、前立腺腺癌と、卵巣癌(漿液性腫瘍を含む卵巣上皮癌又は表面上皮間質性腫瘍、子宮内膜性腫瘍及び粘液性嚢胞腺癌、性索間質性腫瘍)と、肝臓及び胆管癌(肝細胞癌、胆管細胞癌、血管腫を含む)と、食道癌(食道腺癌及び扁平上皮癌を含む)と、口腔及び中喉頭扁平上皮癌と、唾液腺嚢胞癌と、膀胱癌と、膀胱癌と、子宮の癌(子宮内膜腺癌、眼性、子宮乳頭漿液性癌、子宮明細胞癌、子宮肉腫及び平滑筋肉腫、混合ミュラー管腫瘍を含む)と、神経膠腫、神経膠芽腫、髄芽腫、及び脳の他の腫瘍と、腎臓癌(腎細胞癌、明細胞癌、ウィルムス腫瘍を含む)と、頭及び首の癌(扁平上皮癌を含む)と、胃の癌(胃癌、胃腺癌、消化管間質性腫瘍)と、精巣癌と、胚細胞性腫瘍と、神経内分泌腫瘍と、子宮頸癌と、消化管、乳房、及び他の器官のカルチノイドと、印環細胞癌と、肉腫を含む間葉腫、線維肉腫、血管腫、血管腫症、血管周皮細胞腫、偽血管腫間質性過形成、筋組織新生物、線維種症、炎症性筋線維芽細胞腫瘍、脂肪腫、血管脂肪腫、顆粒細胞腫、神経線維種、シュワン腫、血管肉腫、脂肪肉腫、横紋筋肉腫、骨肉腫、平滑筋腫、平滑筋肉腫、皮膚、メラノーマを含む、頸部、網膜芽細胞腫、頭頸部癌、膵臓、脳、甲状腺、精巣、腎臓、膀胱、軟組織、副腎、尿道、陰茎癌、粘液肉腫、軟骨肉腫、骨肉腫、脊索腫、悪性線維性組織球腫、琳派肉腫、中皮腫、扁平上皮癌と、類表皮癌、悪性皮膚付属器腫瘍、腺癌、環細胞腫、環細胞癌、腎細胞癌、副腎腫、胆管細胞癌、移行上皮癌、絨毛癌、精上皮癌、胚的細胞癌、未分化神経膠腫と、多形神経膠芽腫、神経芽細胞腫、髄芽腫、悪性髄膜腫、悪性シュワン腫、神経線維肉腫、副甲状腺癌、甲状腺の髄様癌、気管支カルチノイド、褐色細胞腫、膵島細胞癌、悪性カルチノイド、悪性傍神経節腫、メラノーマ、メルケル細胞新生物、葉状腫瘍、唾液腺癌、胸腺癌と、とりわけ膣癌である。
一実施形態では、本明細書に記載の方法は、リンパ腫又はリンパ球若しくは骨髄球増殖の障害又は異常に罹患している宿主を治療するのに有効である。たとえば、非ホジキンリンパ腫のホジキンリンパ腫に罹患している対象。たとえば、対象は非ホジキンリンパ腫に罹患している場合があり、これには、それだけには限らないが、AIDS関連リンパ腫、未分化大細胞型リンパ腫、血管免疫芽細胞性リンパ腫、芽細胞性NK細胞リンパ腫、バーキットリンパ腫、バーキット様リンパ腫(小型非分割細胞リンパ腫)、慢性リンパ性白血病/小型リンパ球性リンパ腫、皮膚T細胞性リンパ腫、びまん性大細胞型B細胞リンパ腫、腸症型T細胞リンパ腫、濾胞性リンパ腫、肝脾ガンマ・デルタT細胞リンパ腫、リンパ芽球性リンパ腫、マントル細胞リンパ腫、辺縁帯リンパ腫、鼻T細胞リンパ腫、小児リンパ腫、末梢T細胞リンパ腫、中枢神経系リンパ腫、T細胞白血病、形質転換リンパ腫、治療関連T細胞リンパ腫、又はワルデンシュトレーム・マクログロブリン血症などがある。
或いは、対象はホジキンリンパ腫に罹患していることがあり、これには、それだけには限らないが、結節硬化型古典的ホジキンリンパ腫(CHL)、混合細胞型CHL、リンパ球枯渇型CHL、リンパ球豊富型CHL、リンパ球優位型ホジキンリンパ腫、又は結節性リンパ球優位型HLがある。
一実施形態では、対象は、特異的T細胞、B細胞、又はNK細胞によるリンパ腫、増殖性障害、又は異常に罹患していることがある。たとえば、対象は特異的T細胞又はNK細胞リンパ腫をに罹患している場合があり、これには、たとえば、それだけには限らないが、末梢T細胞リンパ腫、たとえば、末梢T細胞リンパ腫、及び特にことわらなければ末梢T細胞リンパ腫(PTCL-NOS)、又は未分化大細胞型リンパ腫、たとえば未分化リンパ腫キナーゼ(ALK)陽性、ALK陰性未分化大細胞型リンパ腫、又は原発性皮膚未分化大細胞型リンパ腫、又は自己免疫芽細胞性リンパ腫、又は皮膚T細胞リンパ腫、たとえば菌状息肉症、セザリー症候群、原発性皮膚未分化大細胞型リンパ腫、原発性皮膚CD30+T細胞リンパ増殖性障害、又は原発性皮膚侵襲性表皮向性CD8+細胞障害性T細胞リンパ腫、又は原発性皮膚ガンマ・デルタT細胞リンパ腫、又は原発性皮膚小/中CD4+T細胞リンパ腫、及びリンパ腫様丘疹症、又は成人T細胞白血病/リンパ腫(ATLL)、又は芽細胞性NK細胞リンパ腫、又は腸症型T細胞リンパ腫、又は肝脾ガンマ・デルタT細胞リンパ腫、又はリンパ芽球性リンパ腫、又は鼻NK/T細胞リンパ腫、又は治療関連T細胞リンパ腫、又はたとえば、実質臓器又は骨髄移植の後に現われるリンパ腫、又はT細胞前リンパ球性白血病、又はT細胞大顆粒リンパ球性白血病、又はNK細胞の慢性リンパ増殖性障害、又は侵襲性NK細胞白血病、又は子供の全身性EBV+T細胞リンパ球増殖性疾患(慢性活動性EBV感染を伴う)、又は種痘様水泡症様リンパ腫、又は成人T細胞白血病/リンパ腫、又は腸症関連T細胞リンパ腫、又は肝脾T細胞リンパ腫、又は皮下脂肪組織炎様T細胞リンパ腫などがある。
或いは、対象は特異的B細胞リンパ腫又は細胞増殖性障害に罹患していることがあり、これには、それだけには限らないが、多発性骨髄腫、又はびまん性大細胞型B細胞リンパ腫、又は濾胞性リンパ腫、又は粘膜関連リンパ組織リンパ腫(MALT)、又は小細胞型リンパ球性リンパ腫、又はマントル細胞リンパ腫(MCL)、又はバーキットリンパ腫、又は縦隔大細胞型B細胞リンパ腫、又はワルデンシュトレーム・マクログロブリン血症、又は結節辺縁帯B細胞リンパ腫(NMZL)、又は脾臓辺縁帯リンパ腫(SMZL)、又は血管内大細胞型B細胞リンパ腫、又は原発性滲出リンパ腫若しくはリンパ腫様肉芽腫症、又は慢性リンパ性白血病/小リンパ球性リンパ腫、又はB細胞前リンパ球性白血病、又はヘアリー細胞白血病、又は脾リンパ腫/白血病、分類不能、又は脾びまん性赤髄小細胞型B細胞リンパ腫、又はヘアリー細胞白血病変種、又はリンパ形質細胞性リンパ腫、又は重鎖疾患、たとえば、α重鎖疾患、γ重鎖疾患、μ重鎖疾患、又は形質細胞性骨髄腫、又は骨の孤立性形質細胞腫、又は骨外性形質細胞腫、又は原発性皮膚濾胞中心リンパ腫、又はT細胞/組織球豊富大細胞型B細胞リンパ腫、又は慢性炎症関連DLBCL、又は高齢者のエプスタイン・バー・ウイルス(EBV)+DLBCL、又は原発性縦隔(胸腺)大細胞型B細胞リンパ腫、又は原発性皮膚DLBCL、脚型、又はALK+大細胞型B細胞リンパ腫、又は形質芽球性リンパ腫、又はHHV8関連多中心性に生じる大細胞型B細胞リンパ腫、又はキャッスルマン病、又はびまん性大細胞型B細胞リンパ腫とバーキットリンパ腫の中間の特徴を持つ、分類不可能なB細胞リンパ腫、又はびまん性大細胞型B細胞リンパ腫と古典的ホジキンリンパ腫の中間の特徴を持つ、分類不可能なB細胞リンパ腫、又は結節硬化型古典的ホジキンリンパ腫、又はリンパ球豊富古典的ホジキンリンパ腫、又は混合細胞性古典的ホジキンリンパ腫、又はリンパ球枯渇古典的ホジキンリンパ腫などがある。
一実施形態では、対象は白血病に罹患している。たとえば、対象は、リンパ球性又は骨髄性起源の急性又は慢性白血病に罹患していることがあり、これには、それだけには限らないが、急性リンパ性白血病(ALL)、又は急性骨髄性白血病(AML)、又は慢性リンパ性白血病(CLL)、又は慢性骨髄性白血病(CML)、又は若年性骨髄単球性白血病(JMML)、又はヘアリー細胞白血病(HCL)、又は急性前骨髄球性白血病(AMLのサブタイプ)、又はT細胞前リンパ性白血病(TPLL)、又は大顆粒リンパ球性白血病、或いは成人T細胞慢性白血病、又は大顆粒リンパ球性白血病(LGL)などがある。一実施形態では、対象は急性骨髄性白血病に罹患しており、これには、たとえば未分化型AML(M0)、又は骨髄芽球白血病(M1、又は最小限細胞成熟あり/なし)、又は骨髄芽球白血病(M2、又は細胞成熟あり)、又は前骨髄球性白血病(M3若しくはM3変種[M3V])、又は骨髄球性白血病(好酸球増加症を伴うM4若しくはM4変種[M4E])、又は単球性白血病(M5)、又は赤白血病(M6)、又は巨核芽球白血病(M7)がある。
特定の実施形態では、対象は、乳癌、肺癌、結腸癌、膵臓癌、又は脳癌に罹患しているか、又は罹患している疑いがある。特定の実施形態では、対象は乳癌に罹患しているか、又は罹患して疑いがある。
癌の特定の実施形態では、臨床事象を経験するリスクが高い患者は、化学療法又は標的治療を、それだけには限らないが、手術、放射線療法、(標的)薬物療法などの標準の治療様式に加えて受けることができる。或いは、臨床事象を経験するリスクが低い患者は、それだけには限らないが、手術、放射線療法、化学療法などの標準の治療様式を控えることができる。
一実施形態では、治療薬を投与するか、それとも治療薬を投与することを控えるかの決定は、閾値MPSスコア、たとえば患者を低リスク群に割り当てるために確立された閾値、又は患者を高リスク群に割り当てるために確立された閾値に基づくことができる。たとえば、一実施形態で、患者を低リスク群に割り当てるための閾値は、5、6、7、8、9、10年又はそれ以上で5%、10%、15%、20%以下の臨床事象のリスクに基づくことができるのに対し、患者を高リスク群に割り当てるための閾値は、5、6、7、8、9、10年又はそれ以上で20%、25%、30%、35%、40%、45%、50%又はそれ以上と等しいか、又はそれを超える臨床事象のリスクに基づくことができる。たとえば、上記の例証を用いると、MPStpwehの特定の場合では、これにより、低リスク患者群に対する閾値が-0.5、-0.4、-0.3、-0.2、-0.1、0になり、高リスク患者群に対する閾値が0、0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9、1、1.1、1.2になる。
本発明の一態様では、対象を低リスク群又は高リスク群に割り当て得ることに関する臨床事象には、癌再発、進行、転位、癌による死亡、又は本明細書の別のところに記載の臨床事象が含まれ得る。
特定の実施形態では、高リスク又は低リスクの割り当ては、乳癌がある対象に対してであり、ER+若しくはHR+腫瘍又は管腔A若しくは管腔Bサブタイプがあって臨床事象を経験するリスクが高い患者は、(術前)補助化学療法を、それだけには限らないが、タモキシフェン又はアロマターゼ阻害剤などのホルモン治療に加えて受けることができる。ER+腫瘍又は管腔A若しくは管腔Bサブタイプがあって臨床事象を経験するリスクが低い患者は、(術前)補助ホルモン治療を受ける(かつ化学療法を控える)ことができる。HER2+/HR-腫瘍又はHER2豊富サブタイプがあって臨床事象を経験するリスクが高い患者は、(術前)補助化学療法を、それだけには限らないが、トラスツズマブなどの抗HER2治療に加えて受けることができるのに対し、HER2+/HR-腫瘍又はHER2豊富サブタイプがあって臨床事象を経験するリスクが低い患者は、(術前)補助抗HER2治療を受ける(かつ化学療法を控える)ことができる。HER2+/HR+腫瘍があって臨床事象を経験するリスクが高い患者は、抗HER2治療と共に(術前)補助化学療法を、それだけには限らないが、タモキシフェン又はアロマターゼ阻害剤などのホルモン治療に加えて受けることができるのに対し、HER2+/HR+腫瘍があって臨床事象を経験するリスクが低い患者は、(術前)補助ホルモン治療を受ける(かつ化学療法及び/又は抗HER2治療を控える)ことができる。三重陰性(HER2-/ER-/PR-又はHER2-/HR-)腫瘍又は基底サブタイプがあって臨床事象を経験するリスクが高い患者は、(術前)補助化学療法を、それだけには限らないが本明細書に記載されているものなどの標的治療に加えて、受けることができるのに対し、三重陰性腫瘍又は基底サブタイプがあって臨床事象を経験するリスクが低い患者は、(術前)補助標的治療を受ける(かつ化学療法を控える)ことができる。
実施例5 リスク・スコアを決定するためのキット及び解析ツール
それぞれの細胞シグナル伝達経路の活性を最善に示すことが、たとえば、ベイジアン・モデル又は(擬似)線形モデルを使用するマイクロアレイ/RNA配列決定ベースの調査に基づいて分かっている標的遺伝子の組は、たとえば、対象の試料に対して実施されるべき定量的多重PCRアッセイ又は専用マイクロアレイ・バイオ・チップに翻訳することができる。本明細書に記載の遺伝子配列から選択されたものを使用して、たとえばRT-PCRのためのプライマー・プローブ・セット、又はマイクロアレイ開発のためのオリゴヌクレオチドを選択することができる。経路活性及びリスク・スコア決定のためのこのようなFDA承認検査を開発するには、標準化検査キットの開発が必要とされ、このキットは、規制上の承認を得るために臨床試験で臨床的に検証される必要がある。
本出願では、いくつかの好ましい実施形態を記載している。変更及び変形は、これまでの詳細な記述を読み理解することにより他の人にも想起されよう。本出願は、このような変更及び変形を、それが添付の特許請求の範囲、又はその等価物の範囲に入る限りにおいて含むと解釈されるものである。
開示された実施形態に対する他の変形形態は、特許請求された本発明を実践するにおいて、図面、開示、及び添付の特許請求の範囲を検討することにより、当業者によって理解され、もたらされ得る。
特許請求の範囲で、「備える」という語は他の要素又はステップを除外せず、また不定冠詞「a」又は「an」は複数を除外しない。
単一のユニット又はデバイスは、特許請求の範囲に列挙されたいくつかの物品の機能を実現することができる。いくつかの方策が互いに異なる独立請求項に列挙されているにすぎないことは、これらの方策の組合せを有利に使用できないことを示すものではない。
1つ又はいくつかのユニット又はデバイスによって実行されるリスク・スコアの決定のような計算は、他の任意の数のユニット又はデバイスによって実行することができる。
コンピュータ・プログラムは、光記憶媒体、又は他のハードウェアと一緒に、若しくはその一部として供給される固体媒体などの、適切な非一時的媒体に記憶/配布することができるが、インターネット、又は他の有線若しくは無線通信システムなどを介する他の形で配布することもできる。
特許請求の範囲におけるいかなる参照符号も、その範囲を限定するものと解釈されるべきではない。
実施例6 適用に使用される配列リスト
Seq.No. Gene:

Seq. 1 ADRA2C
Seq. 2 AGRP
Seq. 3 AP1B1
Seq. 4 ASCL2
Seq. 5 ATG14
Seq. 6 ATP5J
Seq. 7 ATP8A1
Seq. 8 AXIN2
Seq. 9 BCL2
Seq. 10 BCL2L11
Seq. 11 BCL6
Seq. 12 BIRC5
Seq. 13 BMP7
Seq. 14 BNIP3
Seq. 15 BTG1
Seq. 16 C10orf10
Seq. 17 CA12
Seq. 18 CAT
Seq. 19 CAV1
Seq. 20 CBLB
Seq. 21 CCND1
Seq. 22 CCND2
Seq. 23 CCNG2
Seq. 24 CD44
Seq. 25 CDH26
Seq. 26 CDKN1A
Seq. 27 CDKN1B
Seq. 28 CELSR2
Seq. 29 CFLAR
Seq. 30 COL18A1
Seq. 31 COX7A2L
Seq. 32 CTSD
Seq. 33 CTSL
Seq. 34 DDB1
Seq. 35 DEFA6
Seq. 36 DKK1
Seq. 37 DSCAM
Seq. 38 DYRK2
Seq. 39 EBAG9
Seq. 40 EPHB2
Seq. 41 EPHB3
Seq. 42 ERBB2
Seq. 43 ERBB3
Seq. 44 EREG
Seq. 45 ESR1
Seq. 46 EXT1
Seq. 47 FASLG
Seq. 48 FAT1
Seq. 49 FBXO32
Seq. 50 FGFR2
Seq. 51 FOXA2
Seq. 52 FOXF1
Seq. 53 FOXL1
Seq. 54 FOXM1
Seq. 55 FST
Seq. 56 FYN
Seq. 57 FZD7
Seq. 58 GADD45A
Seq. 59 GLI1
Seq. 60 GLI3
Seq. 61 GLUL
Seq. 62 GREB1
Seq. 63 H19
Seq. 64 HHIP
Seq. 65 HNF1A
Seq. 66 HSPB1
Seq. 67 IGF1R
Seq. 68 IGFBP1
Seq. 69 IGFBP3
Seq. 70 IGFBP4
Seq. 71 IGFBP6
Seq. 72 IL1R2
Seq. 73 CXCL8(previously known as IL8)
Seq. 74 INSR
Seq. 75 JAG2
Seq. 76 JUP
Seq. 77 CEMIP(previously known as KIAA1199)
Seq. 78 KLF2
Seq. 79 KLF4
Seq. 80 KLF6
Seq. 81 KRT19
Seq. 82 LECT2
Seq. 83 LEF1
Seq. 84 LGMN
Seq. 85 LGR5
Seq. 86 MIF
Seq. 87 MXI1
Seq. 88 MYC
Seq. 89 MYCN
Seq. 90 MYLK
Seq. 91 MYOD1
Seq. 92 NDUFV3
Seq. 93 NKD1
Seq. 94 NKX2-2
Seq. 95 NKX2-8
Seq. 96 NOS3
Seq. 97 NRIP1
Seq. 98 OAT
Seq. 99 PCK1
Seq. 100 PDK4
Seq. 101 PGR
Seq. 102 PISD
Seq. 103 PITRM1
Seq. 104 POMC
Seq. 105 PPARG
Seq. 106 PPARGC1A
Seq. 107 PPM1D
Seq. 108 PRDM15
Seq. 109 PRDX3
Seq. 110 PTCH1
Seq. 111 PTCH2
Seq. 112 PTMA
Seq. 113 RAB34
Seq. 114 RAG1
Seq. 115 RAG2
Seq. 116 RARA
Seq. 117 RBL2
Seq. 118 REG1B
Seq. 119 RNF43
Seq. 120 S100A7
Seq. 121 S100A9
Seq. 122 SEMA3C
Seq. 123 SEPP1
Seq. 124 SESN1
Seq. 125 SGK3
Seq. 126 SIRT1
Seq. 127 SLC1A2
Seq. 128 SLC5A3
Seq. 129 SMAD4
Seq. 130 SOD1
Seq. 131 SOD2
Seq. 132 SOX9
Seq. 133 SP5
Seq. 134 SPP1
Seq. 135 STK11
Seq. 136 TBX3
Seq. 137 TCEA2
Seq. 138 TCF7L2
Seq. 139 TDGF1
Seq. 140 TFF1
Seq. 141 TLE4
Seq. 142 TNFSF10
Seq. 143 TOM1
Seq. 144 TRIM25
Seq. 145 TSC22D1
Seq. 146 TXNIP
Seq. 147 WISP2
Seq. 148 XBP1
Seq. 149 ZNRF3
Seq. 150 SERPINE1
Seq. 151 PDZK1
Figure 0007065609000022
Figure 0007065609000023
Figure 0007065609000024
Figure 0007065609000025
Figure 0007065609000026

Claims (13)

  1. 対象が疾患に関連した臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための、デジタル処理デバイスによって実行される方法であって、
    前記デジタル処理デバイスが、前記対象における2つ以上の細胞シグナル伝達経路のそれぞれの活性を、前記対象の試料において測定されたそれぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子の発現レベルを、前記対象の試料における前記それぞれの細胞シグナル伝達経路のそれぞれの転写因子(TF)要素の活性レベルに関係付ける較正済み数学経路モデルを評価することに基づいて推測するステップであり、前記TF要素は、前記それぞれの細胞シグナル伝達経路の3つ以上の標的遺伝子の転写を制御する、ステップと、
    前記デジタル処理デバイスが、前記推測された活性を組み合わせたものに基づいて、前記リスク・スコアを決定するステップと、
    を含み、
    前記臨床事象は、疾患再発、疾患進行、及び疾患が原因の死亡のうちの1つであり、前記疾患は乳癌であり、
    前記細胞シグナル伝達経路がPI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上の経路とを含
    前記PI3K経路、前記Wnt経路、前記ER経路、及び前記HH経路は、前記それぞれのTF要素の転写活性に最終的につながる細胞シグナル伝達経路として定義され、
    PI3K TF要素はFOXOファミリーメンバーを含み、Wnt TF要素はβカテニン/TCF4を含み、ER TF要素はERα二量体を含み、HH TF要素はGLIファミリーメンバーを含む、方法。
  2. 前記細胞シグナル伝達経路が、前記Wnt経路及び/又は前記ER経路及び/又は前記HH経路を含み、前記リスク・スコアが、前記示されたリスクが前記Wnt経路の推測された活性の増加、及び/又は前記HH経路の推測された活性の増加と共に単調増加し、かつ/又は、前記ER経路の推測された活性の増加と共に単調減少するように規定される、請求項1に記載の方法。
  3. 前記リスク・スコアが、前記示されたリスクが前記PI3K細胞シグナル伝達経路の推測された活性の増加と共に単調増加するように規定される、請求項1又は2に記載の方法。
  4. 前記リスク・スコアが、項w・Pと、項w・P、w・Pe、及びw・Pのうちの1つ以上とを含む合計から成る複数経路スコア(MPS)に基づき、ここで、P、P、P、及びPが、それぞれ前記PI3K経路、前記Wnt経路、前記ER経路、及び前記HH経路の前記推測された活性を表し、w、w、w、及びwが、前記対象が前記臨床事象を前記規定期間内に経験する前記リスクと、前記PI3K経路、前記Wnt経路、前記ER経路、及び前記HH経路それぞれの前記活性との間の相関関係を表す一定重み付け係数である、請求項1から3のいずれか一項に記載の方法。
  5. 前記一定重み付け係数w、w、w、及びwが、前記それぞれの細胞シグナル伝達経路のコックス比例ハザード・モデルを臨床データに適合させることにより得られるコックスの係数の値に基づいてそれぞれ決定されるか、又は決定されている、請求項4に記載の方法。
  6. 前記3つ以上のPI3K標的遺伝子が、AGRP、BCL2L11、BCL6、BNIP3、BTG1、CAT、CAV1、CCND1、CCND2、CCNG2、CDKN1A、CDKN1B、ESR1、FASLG、FBXO32、GADD45A、INSR、MXI1、NOS3、PCK1、POMC、PPARGC1A、PRDX3、RBL2、SOD2、TNFSF10から成る群から選択されるか、又は、FBXO32、BCL2L11、SOD2、TNFSF10、BCL6、BTG1、CCNG2、CDKN1B、BNIP3、GADD45A、INSR、及びMXI1から成る群から選択され、
    かつ/又は
    前記3つ以上のWnt標的遺伝子が、CEMIP、AXIN2、CD44、RNF43、MYC、TBX3、TDGF1、SOX9、ASCL2、CXCL8、SP5、ZNRF3、EPHB2、LGR5、EPHB3、KLF6、CCND1、DEFA6、及びFZD7から成る群から選択されるか、又は、AXIN2、CD44、LGR5、CEMIP、MYC、CXCL8、SOX9、EPHB3、RNF43、TDGF1、ZNRF3、及びDEFA6から成る群から選択され、
    かつ/又は
    前記3つ以上のER標的遺伝子が、CDH26、SGK3、PGR、GREB1、CA12、XBP1、CELSR2、WISP2、DSCAM、ERBB2、CTSD、TFF1、PDZK1、IGFBP4、ESR1、SOD1、AP1B1、及びNRIP1から成る群から選択されるか、又は、TFF1、GREB1、PGR、SGK3、PDZK1、IGFBP4、NRIP1、CA12、XBP1、ERBB2、ESR1、及びCELSR2から成る群から選択され、
    かつ/又は
    前記3つ以上のHH標的遺伝子が、GLI1、PTCH1、PTCH2、IGFBP6、SPP1、CCND2、FST、FOXL1、CFLAR、TSC22D1、RAB34、S100A9、S100A7、MYCN、FOXM1、GLI3、TCEA2、FYN、及びCTSL1から成る群から選択されるか、又は、GLI1、PTCH1、PTCH2、CCND2、IGFBP6、MYCN、FST、RAB34、GLI3、CFLAR、S100A7、及びS100A9から成る群から選択される、請求項1から5のいずれか一項に記載の方法。
  7. 前記デジタル処理デバイスが、前記対象を、前記対象が臨床事象を規定期間内に経験することの別々に示されたリスクと関連付けられた複数のリスク群のうちの少なくとも1つに割り当てるステップ、
    及び/又は
    前記デジタル処理デバイスが、前記対象に推奨される治療を、前記対象が臨床事象を規定期間内に経験することの前記示されたリスクに基づいて決定するステップをさらに含む、請求項1からのいずれか一項に記載の方法。
  8. 対象が乳癌と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための装置であって、請求項1からのいずれか一項に記載の方法を実施するように構成されたデジタル・プロセッサを備える装置。
  9. 対象が乳癌と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するための非一時的記憶媒体であって、請求項1からのいずれか一項に記載の方法を実施するデジタル処理デバイスによって実行可能である命令を記憶する非一時的記憶媒体。
  10. 対象が乳癌と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのコンピュータ・プログラムであって、前記コンピュータ・プログラムがデジタル処理デバイス上で実行されたとき、前記デジタル処理デバイスに請求項1からのいずれか一項に記載の方法を実施させるプログラム・コード手段を含む、コンピュータ・プログラム。
  11. 対象が乳癌と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのキットであって、
    対象の試料における2つ以上の細胞シグナル伝達経路それぞれの、3つ以上の標的遺伝子の発現レベルを測定するためのキットであり、前記それぞれの細胞シグナル伝達経路の前記3つ以上の標的遺伝子を対象とするポリメラーゼ連鎖反応プライマーと、前記それぞれの細胞シグナル伝達経路の前記3つ以上の標的遺伝子を対象とするプローブと、を含み、前記細胞シグナル伝達経路がPI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上の経路とを含む、キットと、
    請求項に記載の装置、請求項に記載の非一時的記憶媒体、又は請求項10に記載のコンピュータ・プログラムと、
    を含む、キット。
  12. 対象の試料において2つ以上の細胞シグナル伝達経路それぞれの3つ以上の標的遺伝子の発現レベルを測定するための、及び、対象が乳癌と関連する臨床事象を規定期間内に経験するリスクを示すリスク・スコアを決定するためのキットであって、
    前記対象の試料において2つ以上の細胞シグナル伝達経路それぞれの、3つ以上の標的遺伝子の発現レベルを決定するための1つ又は複数の構成要素を含み、
    前記1つ又は複数の構成要素が、マイクロアレイ・チップ、抗体、複数のプローブ、たとえば標識されたプローブ、RNA逆転写酵素を用いた配列決定に使用される構成要素の組、及び/又は、RNA又はcDNAを含むDNAの増幅プライマーから成る群から選択され、
    前記細胞シグナル伝達経路がPI3K経路と、Wnt経路、ER経路、及びHH経路のうちの1つ以上の経路とを含み
    求項に記載の装置、請求項に記載の非一時的記憶媒体、又は請求項10に記載のコンピュータ・プログラムを含むキット。
  13. 請求項1からのいずれか一項に記載の方法を実施する際の、請求項11又は12に記載のキットの使用。
JP2017522163A 2014-10-24 2015-10-26 複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測 Active JP7065609B6 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP14190273 2014-10-24
EP14190273.4 2014-10-24
PCT/EP2015/074701 WO2016062892A1 (en) 2014-10-24 2015-10-26 Medical prognosis and prediction of treatment response using multiple cellular signaling pathway activities

Publications (3)

Publication Number Publication Date
JP2018502553A JP2018502553A (ja) 2018-02-01
JP7065609B2 true JP7065609B2 (ja) 2022-05-12
JP7065609B6 JP7065609B6 (ja) 2022-06-06

Family

ID=51846475

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017522163A Active JP7065609B6 (ja) 2014-10-24 2015-10-26 複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測

Country Status (8)

Country Link
US (1) US11640845B2 (ja)
EP (1) EP3210143B1 (ja)
JP (1) JP7065609B6 (ja)
CN (1) CN108064311B (ja)
AU (1) AU2015334841B2 (ja)
CA (1) CA2965217A1 (ja)
ES (1) ES2857953T3 (ja)
WO (1) WO2016062892A1 (ja)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11045352B2 (en) * 2014-05-12 2021-06-29 Gholam A. Peyman Methods for treatment of dry eye and other acute or chronic inflammatory processes
US11648261B2 (en) 2014-05-12 2023-05-16 Gholam A. Peyman Method of treating, reducing, or alleviating a medical condition in a patient
US10925889B2 (en) 2014-05-12 2021-02-23 Gholam A. Peyman Method of treating, reducing, or alleviating a medical condition in a patient
US20190128899A1 (en) * 2016-06-13 2019-05-02 Koninklijke Philips N.V. Method for inferring activity of a transcription factor of a signal transduction pathway in a subject
US10621236B2 (en) * 2016-09-16 2020-04-14 At&T Intellectual Property I, L.P. Concept based querying of graph databases
WO2018096076A1 (en) * 2016-11-25 2018-05-31 Koninklijke Philips N.V. Method to distinguish tumor suppressive foxo activity from oxidative stress
EP3431582A1 (en) 2017-07-18 2019-01-23 Koninklijke Philips N.V. Cell culturing materials
WO2019033013A1 (en) * 2017-08-11 2019-02-14 Eric Meittunen MANAGEMENT OF PATIENT PANELS OF PRIMARY CARE
EP3462348A1 (en) * 2017-09-28 2019-04-03 Koninklijke Philips N.V. Bayesian inference
EP3692170A1 (en) 2017-10-02 2020-08-12 Koninklijke Philips N.V. Determining functional status of immune cells types and immune response
EP3503118A1 (en) * 2017-12-21 2019-06-26 Koninklijke Philips N.V. Methods and apparatus for reducing risk to a subject undergoing radiotherapy-based treatment
CN110557808A (zh) * 2018-06-01 2019-12-10 北京嘀嘀无限科技发展有限公司 Tbox控制方法、装置及设备
CN109243535A (zh) * 2018-09-13 2019-01-18 河南财经政法大学 同步编程模型Hama BSP下基于蝶形网络的双聚类挖掘方法
CN109439753B (zh) * 2018-11-28 2022-05-06 四川大学华西医院 检测基因表达水平的试剂的应用以及乳腺癌患者nac疗效预测模型的构建方法
WO2020168290A1 (en) * 2019-02-14 2020-08-20 University Of Kentucky Research Foundation N-aryl benzenesulfonamides for use in treating cancers, bacterial diseases, metabolic diseases, and traumatic brain injury
CN110157792A (zh) * 2019-04-22 2019-08-23 中山大学孙逸仙纪念医院 血清外泌体has_circ_0004771在制备酒精依赖综合征诊断试剂中的应用
US11707518B2 (en) 2019-04-28 2023-07-25 Gholam A. Peyman Method of treating, reducing, or alleviating a medical condition in a patient
US11674185B2 (en) * 2019-05-03 2023-06-13 Koninklijke Philips N.V. Methods of prognosis in high-grade serous ovarian cancer
EP3812474A1 (en) 2019-10-22 2021-04-28 Koninklijke Philips N.V. Methods of prognosis in high-grade serous ovarian cancer
EP3739588A1 (en) 2019-05-13 2020-11-18 Koninklijke Philips N.V. Assessment of multiple signaling pathway activity score in airway epithelial cells to predict airway epithelial abnormality and airway cancer risk
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
WO2021209567A1 (en) 2020-04-16 2021-10-21 Koninklijke Philips N.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
EP3965119A1 (en) 2020-09-04 2022-03-09 Koninklijke Philips N.V. Methods for estimating heterogeneity of a tumour based on values for two or more genome mutation and/or gene expression related parameter, as well as corresponding devices
EP3974540A1 (en) 2020-09-25 2022-03-30 Koninklijke Philips N.V. Method for predicting immunotherapy resistance
EP4255578A1 (en) * 2020-12-07 2023-10-11 Foundation Medicine, Inc. Comprehensive genomic profiling (cgp) of metastatic invasive lobular carcinomas reveals heterogeneity
EP4015651A1 (en) 2020-12-17 2022-06-22 Koninklijke Philips N.V. Treatment prediction and effectiveness of anti-tnf alpha treatment in ibd patients
CN112646864A (zh) * 2020-12-30 2021-04-13 杭州联川基因诊断技术有限公司 一种检测esr1基因表达的引物、探针、试剂盒和检测方法
EP4039825A1 (en) 2021-02-09 2022-08-10 Koninklijke Philips N.V. Comparison and standardization of cell and tissue culture
WO2022189530A1 (en) 2021-03-11 2022-09-15 Koninklijke Philips N.V. Prognostic pathways for high risk sepsis patients

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007532113A (ja) 2004-04-09 2007-11-15 ジェノミック ヘルス, インコーポレイテッド 化学療法剤に対する応答を予測するための遺伝子発現マーカー
US20090186024A1 (en) 2005-05-13 2009-07-23 Nevins Joseph R Gene Expression Signatures for Oncogenic Pathway Deregulation
JP2012525159A (ja) 2009-05-01 2012-10-22 ジェノミック ヘルス, インコーポレイテッド 結腸直腸癌の再発および化学療法に対する応答の可能性における遺伝子発現プロファイルアルゴリズムおよび試験
JP2014520566A (ja) 2011-07-19 2014-08-25 コーニンクレッカ フィリップス エヌ ヴェ 標的遺伝子発現の確率モデルを使用した細胞シグナル経路の活性の評価

Family Cites Families (48)

* 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
US6004761A (en) 1986-11-19 1999-12-21 Sanofi Method for detecting cancer using monoclonal antibodies to new mucin epitopes
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
US5545531A (en) 1995-06-07 1996-08-13 Affymax Technologies N.V. Methods for making a device for concurrently processing multiple biological chip assays
US6720149B1 (en) 1995-06-07 2004-04-13 Affymetrix, Inc. Methods for concurrently processing multiple biological chip assays
US6146897A (en) 1995-11-13 2000-11-14 Bio-Rad Laboratories Method for the detection of cellular abnormalities using Fourier transform infrared spectroscopy
US6391550B1 (en) 1996-09-19 2002-05-21 Affymetrix, Inc. Identification of molecular sequence signatures and methods involving the same
NZ516848A (en) 1997-06-20 2004-03-26 Ciphergen Biosystems Inc Retentate chromatography apparatus with applications in biology and medicine
JP2001515234A (ja) 1997-07-25 2001-09-18 アフィメトリックス インコーポレイテッド 多型性データベースを提供するためのシステム
US6953662B2 (en) 1997-08-29 2005-10-11 Human Genome Sciences, Inc. Follistatin-3
US6020135A (en) 1998-03-27 2000-02-01 Affymetrix, Inc. P53-regulated genes
US6884578B2 (en) 2000-03-31 2005-04-26 Affymetrix, Inc. Genes differentially expressed in secretory versus proliferative endometrium
JP2004522980A (ja) 2000-11-16 2004-07-29 シファーゲン バイオシステムズ, インコーポレイテッド 質量スペクトルを分析する方法
CA2438391A1 (en) 2001-02-16 2002-10-10 William E. Rich Method for correlating gene expression profiles with protein expression profiles
EP2799555B1 (en) 2002-03-13 2017-02-22 Genomic Health, Inc. Gene expression profiling in biopsied tumor tissues
US7097976B2 (en) 2002-06-17 2006-08-29 Affymetrix, Inc. Methods of analysis of allelic imbalance
JP4606879B2 (ja) 2002-11-15 2011-01-05 ジェノミック ヘルス, インコーポレイテッド Egfr陽性癌の遺伝子発現プロファイリング
US20040231909A1 (en) 2003-01-15 2004-11-25 Tai-Yang Luh Motorized vehicle having forward and backward differential structure
ES2609234T3 (es) 2003-06-24 2017-04-19 Genomic Health, Inc. Predicción de la probabilidad de recidiva de cáncer
CA2531967C (en) 2003-07-10 2013-07-16 Genomic Health, Inc. Expression profile algorithm and test for cancer prognosis
WO2006052862A1 (en) 2004-11-05 2006-05-18 Genomic Health, Inc. Predicting response to chemotherapy using gene expression markers
CA2602951A1 (en) 2005-03-23 2006-09-28 Bio-Rad Laboratories, Inc. Diverse chemical libraries bound to small particles with paramagnetic properties
US20060234911A1 (en) 2005-03-24 2006-10-19 Hoffmann F M Method of reversing epithelial mesenchymal transition
KR100806274B1 (ko) 2005-12-06 2008-02-22 한국전자통신연구원 멀티 쓰레디드 프로세서 기반의 병렬 시스템을 위한 적응형실행 방법
NZ593226A (en) 2006-01-11 2012-10-26 Genomic Health Inc Gene expression markers (efnb2) for colorectal cancer prognosis
WO2007123772A2 (en) 2006-03-31 2007-11-01 Genomic Health, Inc. Genes involved in estrogen metabolism
US8518639B2 (en) 2006-04-11 2013-08-27 Bio-Rad Innovations HPV detection and quantification by real-time multiplex amplification
WO2008115419A2 (en) * 2007-03-15 2008-09-25 Genomic Health, Inc. Gene expression markers for prediction of patient response to chemotherapy
US7816084B2 (en) 2007-11-30 2010-10-19 Applied Genomics, Inc. TLE3 as a marker for chemotherapy
WO2009114836A1 (en) 2008-03-14 2009-09-17 Genomic Health, Inc. Gene expression markers for prediction of patient response to chemotherapy
WO2009124251A1 (en) 2008-04-03 2009-10-08 Sloan-Kettering Institute For Cancer Research Gene signatures for the prognosis of cancer
US20120009581A1 (en) 2008-07-08 2012-01-12 Bankaitis-Davis Danute M Gene Expression Profiling for Predicting the Survivability of Prostate Cancer Subjects
US8541170B2 (en) 2008-11-17 2013-09-24 Veracyte, Inc. Methods and compositions of molecular profiling for disease diagnostics
US8765383B2 (en) 2009-04-07 2014-07-01 Genomic Health, Inc. Methods of predicting cancer risk using gene expression in premalignant tissue
EP2460005A4 (en) 2009-07-31 2012-11-21 Translational Genomics Res Inst METHOD FOR ASSESSING A CANCER PROGRESSION RISK
US8451450B2 (en) 2009-09-14 2013-05-28 Bio-Rad Laboratories, Inc. Near real time optical phase conjugation
CN102859528A (zh) 2010-05-19 2013-01-02 加利福尼亚大学董事会 使用生物网络识别药物靶点的系统和方法
US8703736B2 (en) 2011-04-04 2014-04-22 The Translational Genomics Research Institute Therapeutic target for pancreatic cancer cells
US9970057B2 (en) 2011-05-06 2018-05-15 Albert Einstein College Of Medicine, Inc. Human invasion signature for prognosis of metastatic risk
WO2013003384A1 (en) * 2011-06-27 2013-01-03 Dana-Farber Cancer Institute, Inc. Signatures and determinants associated with prostate cancer progression and methods of use thereof
US20140303133A1 (en) 2011-11-18 2014-10-09 Vanderbilt University Markers of Triple-Negative Breast Cancer And Uses Thereof
WO2013084144A1 (en) 2011-12-08 2013-06-13 Koninklijke Philips Electronics N.V. Semiconductor light emitting device with thick metal layers
US8725426B2 (en) 2012-01-31 2014-05-13 Genomic Health, Inc. Gene expression profile algorithm and test for determining prognosis of prostate cancer
AU2013368945B2 (en) 2012-12-26 2020-01-23 Innosign B.V. Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions
CA2909991A1 (en) 2013-04-26 2014-10-30 Koninklijke Philips N.V. Medical prognosis and prediction of treatment response using multiple cellular signalling pathway activities
EP3013986B1 (en) 2014-01-03 2016-11-16 Koninklijke Philips N.V. Assessment of the pi3k cellular signaling pathway activity using mathematical modelling of target gene expression

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007532113A (ja) 2004-04-09 2007-11-15 ジェノミック ヘルス, インコーポレイテッド 化学療法剤に対する応答を予測するための遺伝子発現マーカー
US20090186024A1 (en) 2005-05-13 2009-07-23 Nevins Joseph R Gene Expression Signatures for Oncogenic Pathway Deregulation
JP2012525159A (ja) 2009-05-01 2012-10-22 ジェノミック ヘルス, インコーポレイテッド 結腸直腸癌の再発および化学療法に対する応答の可能性における遺伝子発現プロファイルアルゴリズムおよび試験
JP2014520566A (ja) 2011-07-19 2014-08-25 コーニンクレッカ フィリップス エヌ ヴェ 標的遺伝子発現の確率モデルを使用した細胞シグナル経路の活性の評価

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BMC Cancer,2010年,Vol.10,No.604,pp.1-20

Also Published As

Publication number Publication date
CA2965217A1 (en) 2016-04-28
US20160117443A1 (en) 2016-04-28
WO2016062892A1 (en) 2016-04-28
CN108064311A (zh) 2018-05-22
EP3210143B1 (en) 2020-12-09
JP2018502553A (ja) 2018-02-01
ES2857953T3 (es) 2021-09-29
JP7065609B6 (ja) 2022-06-06
CN108064311B (zh) 2022-10-28
AU2015334841A1 (en) 2017-06-15
US11640845B2 (en) 2023-05-02
AU2015334841B2 (en) 2022-02-03
EP3210143A1 (en) 2017-08-30

Similar Documents

Publication Publication Date Title
JP7065609B2 (ja) 複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測
JP7065610B2 (ja) 複数の細胞シグナル伝達経路活性を用いる治療応答の医学的予後及び予測
JP6721665B2 (ja) 目的遺伝子発現の線形結合を用いた細胞信号伝達経路活性の評価
JP6603208B2 (ja) 複数の細胞内シグナル伝達経路活性を用いた治療反応の医学的予後及び予測
JP6077178B2 (ja) 標的遺伝子発現の数学的モデリングを用いるpi3k細胞内シグナル伝達経路活性の評価
JP6204354B2 (ja) 標的遺伝子発現の確率モデルを使用した細胞シグナル経路の活性の評価

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170622

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181024

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190924

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20191216

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200319

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20200623

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201022

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20201022

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20201029

C21 Notice of transfer of a case for reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C21

Effective date: 20201104

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20201204

C211 Notice of termination of reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C211

Effective date: 20201208

C22 Notice of designation (change) of administrative judge

Free format text: JAPANESE INTERMEDIATE CODE: C22

Effective date: 20210706

C22 Notice of designation (change) of administrative judge

Free format text: JAPANESE INTERMEDIATE CODE: C22

Effective date: 20220118

C23 Notice of termination of proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C23

Effective date: 20220301

C03 Trial/appeal decision taken

Free format text: JAPANESE INTERMEDIATE CODE: C03

Effective date: 20220405

C30A Notification sent

Free format text: JAPANESE INTERMEDIATE CODE: C3012

Effective date: 20220405

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220426

R150 Certificate of patent or registration of utility model

Ref document number: 7065609

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150