Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP6963009B2 - 臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法 - Google Patents
[go: Go Back, main page]

JP6963009B2 - 臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法 - Google Patents

臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法 Download PDF

Info

Publication number
JP6963009B2
JP6963009B2 JP2019512963A JP2019512963A JP6963009B2 JP 6963009 B2 JP6963009 B2 JP 6963009B2 JP 2019512963 A JP2019512963 A JP 2019512963A JP 2019512963 A JP2019512963 A JP 2019512963A JP 6963009 B2 JP6963009 B2 JP 6963009B2
Authority
JP
Japan
Prior art keywords
cancer
human
scar
genes
gene
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
JP2019512963A
Other languages
English (en)
Other versions
JP2019520851A5 (ja
JP2019520851A (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.)
Oncoscar LLC
Original Assignee
Oncoscar LLC
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 Oncoscar LLC filed Critical Oncoscar LLC
Publication of JP2019520851A publication Critical patent/JP2019520851A/ja
Publication of JP2019520851A5 publication Critical patent/JP2019520851A5/ja
Application granted granted Critical
Publication of JP6963009B2 publication Critical patent/JP6963009B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61KPREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
    • A61K31/00Medicinal preparations containing organic active ingredients
    • A61K31/33Heterocyclic compounds
    • A61K31/395Heterocyclic compounds having nitrogen as a ring hetero atom, e.g. guanethidine or rifamycins
    • A61K31/495Heterocyclic compounds having nitrogen as a ring hetero atom, e.g. guanethidine or rifamycins having six-membered rings with two or more nitrogen atoms as the only ring heteroatoms, e.g. piperazine or tetrazines
    • A61K31/505Pyrimidines; Hydrogenated pyrimidines, e.g. trimethoprim
    • A61K31/513Pyrimidines; Hydrogenated pyrimidines, e.g. trimethoprim having oxo groups directly attached to the heterocyclic ring, e.g. cytosine
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61KPREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
    • A61K31/00Medicinal preparations containing organic active ingredients
    • A61K31/70Carbohydrates; Sugars; Derivatives thereof
    • A61K31/7088Compounds having three or more nucleosides or nucleotides
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61PSPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
    • A61P35/00Antineoplastic agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61PSPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
    • A61P43/00Drugs for specific purposes, not provided for in groups A61P1/00-A61P41/00
    • 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/70Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving virus or bacteriophage
    • C12Q1/701Specific hybridization probes
    • C12Q1/702Specific hybridization probes for retroviruses
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • 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
    • 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/106Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
    • 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/156Polymorphic or mutational 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/158Expression markers
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Genetics & Genomics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Molecular Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Immunology (AREA)
  • Public Health (AREA)
  • Analytical Chemistry (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Pharmacology & Pharmacy (AREA)
  • Medical Informatics (AREA)
  • Medicinal Chemistry (AREA)
  • Epidemiology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Evolutionary Biology (AREA)
  • Pathology (AREA)
  • Microbiology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Virology (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Chemical & Material Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Databases & Information Systems (AREA)
  • Bioethics (AREA)

Description

関連出願への相互参照
本出願は、米国仮出願第62/339007号(2016年5月19日出願)の利益を主張し、この出願はその全体が参照により本明細書に組み込まれる。
以下の関連出願の全てもまた、それらの全体が参照により組み込まれる:米国仮出願第60/875,061号(2006年12月15日出願);米国仮出願第60/823,577号(2006年8月25日出願);米国仮出願第60/822,705号(2006年8月17日出願);及び米国仮出願第60/787,818号(2006年3月31日出願)。
概要
一態様において、本開示は、とりわけ、患者内のがんの存在を診断するための、がんを有する被験体が様々な異なるタイプの処置レジメンの影響を受けやすいかどうかを決定するための、患者内のがんの処置を監視するための、新規の方法及びキットを対象とし、個別化され標的化されたがん治療を含むがん治療を提供する新規の方法を提供する。試験、監視及び処置されるべきがんは、前立腺がん、乳がん、肺がん、胃がん、卵巣がん、膀胱がん、リンパ腫、中皮腫、脳がん、肝臓がん、上記のうちいずれかの転移、及び血液がん(ALL、AML、及びCCLを含むがこれらに限定されない)を含むが、これらに限定されない。処置レジメンの初期に治療抵抗性である可能性の高い患者の同定は、より成功した成績を達成するための治療の変更につながる可能性がある。
一態様において、本開示は、とりわけ、がんを診断する又はがん治療成績を予測するための方法を対象とし、この方法は、同じ細胞において同時に、細胞の集団において、又はリキッドバイオプシー標本において、複数のマーカーの配列及び/又は発現レベルを検出する工程、及びそれらの配列及び/又は発現を、ある特定の閾値に関して定性的に異なる又は定量的に異なる(上又は下)ものであるとスコア付けする工程(ここで、マーカーはがんに関連する特定の経路に由来し、スコアはがん診断又はがん治療失敗の予後を示す)によるものである。この方法は、様々ながんについて、がんを診断する又はがん治療成績を予測するために使用することができる。一実施形態において、本方法は、遺伝的シグネチャー情報及びタンパク質シグネチャー情報を使用することによって、個体がSCARのネットワーク活性化を経験しているかどうかを決定する工程を含む。
一態様において、本開示は、とりわけ、がん前駆病変の初期の検出を含む、内在性ヒト幹細胞関連レトロウイルス(SCAR)のゲノム及びプロテオミクスシグネチャーの同定及び監視に基づく、臨床的に難治性の悪性腫瘍の診断及び治療標的化の新規の方法を対象とする。マーカーは、具体的にはSCARの経路及び「幹性」経路を含む、がんの制御に関与している任意の経路に由来することができる。マーカーは、mRNA、RNA、DNA、タンパク質、又はペプチドであることができる。一態様において、本開示は、とりわけ、内在性ヒト幹細胞関連レトロウイルス(SCAR)のゲノム及びプロテオミクスシグネチャーに基づいて、臨床的に難治性の悪性腫瘍のための処置を設計及び使用する新規の方法を対象とする。核酸、DNA、RNAなどを、一塩基ミスマッチ特異性で検出するための技術及び方法論の非限定的な例は、J.S. Gootenberg et al., "Nucleic acid detection with CRISPR-Cas13a/C2c2," Science, doi:10.1126/science.aam9321, 2017(これは、その全体が参照により本明細書に組み込まれる)において説明されるものを含む。
一態様において、本開示は、とりわけ、患者内のがんの存在を診断するための、がんを有する被験体が様々な異なるタイプの処置レジメンの影響を受けやすいかどうかを決定するための、患者内のがんの処置を監視するための、方法及びキットを対象としており、個別化され標的化されたがん治療を含むがん治療を提供する新規の方法を提供する。試験、監視及び処置されるべきがんは、前立腺がん、乳がん、肺がん、胃がん、卵巣がん、膀胱がん、リンパ腫、中皮腫、脳がん、肝臓がん、上記のうちいずれかの転移、及び血液がん(ALL、AML、及びCCLを含むがこれらに限定されない)を含むが、これらに限定されない。全体として、これらの方法の潜在的な実用的有用性が、ヒトがんの29の異なるタイプについて証明された。
一実施形態において、方法は、複数のマーカーの配列、複数のマーカーの発現レベルを、同じ細胞において同時に、細胞の集団において、又はリキッドバイオプシー標本において、同時に又は順次に検出する工程、及びそれらの配列及び/又は発現を異常であるとスコア付けする工程(ここで、マーカーはがんに関連する特定の経路に由来し、スコアはがん診断又はがん治療失敗の可能性の予後を示す)を含む。この方法は、様々ながんについて、がんを診断する又はがん治療成績を予測するために使用することができる。同じ細胞、細胞の集団、又は被験体からのリキッドバイオプシー標本における、少なくとも1つであるが好ましくは2つ以上のマーカーの同時共発現は、がんの診断指標及び被験体が標準的ながん治療に対して抵抗性であることの予測因子である。マーカーは、具体的にはSCARの経路、PcG経路及び「幹性」経路を含む、がんの制御に関与している任意の経路に由来することができる。マーカーは、mRNA、RNA、DNA、タンパク質、又はペプチドであることができる。
一態様において、本開示は、とりわけ、SCARの経路からの複数のマーカーの、閾値レベルを超える、同じ細胞における同時の発現(ここで、マーカーはがんに関連する経路内に見出される)が、がんを診断するための及びがんを有すると既に診断された患者が治療反応性又は治療抵抗性となるかどうかを予測するためのアッセイとして使用することができる、という新規の発見を対象とする。アッセイの要素は、少なくとも1つであるが好ましくは2つ以上のマーカーが、同じ細胞内で、細胞の集団内で、又はリキッドバイオプシー標本において、同時に検出されることである。マーカー検出は、免疫蛍光によるバーコーディング及び次世代シーケンシングを含む様々な検出手段により行うことができる。検出されるマーカーは、mRNA、RNA、DNA、タンパク質、及びペプチドを含む様々な産物であることができる。mRNA、RNA、及びDNAベースのマーカーについては、次世代シーケンシング及び/又はPCRを検出手段として使用することができる。加えて、核酸配列、タンパク質配列、タンパク質産物又は遺伝子コピー数は、本分野において既知の検出手段により同定することができる。検出されるマーカーは、がんに関連する様々な経路に由来することができる。マーカーの適切な経路は、発がん及び転移に関連する任意の経路を含み、より具体的には、SCARの経路、Polycomb群(PcG)クロマチンサイレンシング経路及び「幹性」経路を含む。
一態様において、本開示は、とりわけ、生体被験体においてがんを診断する又はがん治療成績を予測するための方法を対象とする。
一実施形態において、本方法は、生体サンプル(例えば、組織、細胞、体液の標本、生体液、バイオマーカー組成物など)を被験体から得る工程を含む。
一実施形態において、本方法は、マーカーをがんに関連する経路から選択する工程を含み、
一実施形態において、本方法は、少なくとも1つであるが好ましくは2つ以上のマーカーの同時の異常な配列及び/又は発現レベルについてスクリーニングする工程を含み、
一実施形態において、本方法は、それらの配列を、配列の質(配列全体又はその断片内の塩基の位置の規定される配列)が参照配列と比較して異なる場合に、異常であるとスコア付けする工程を含み、
一実施形態において、本方法は、検出された発現レベルがある特定の閾値を上回る場合に、それらの発現レベルを異常であるとスコア付けする工程を含む。
一実施形態において、本方法は、少なくとも1つであるが好ましくは2つ以上のそのようなマーカーの異常な配列及び/又は異常な発現レベルの存在が、被験体におけるがん診断又はがん治療失敗の予後を示すことを含む。
一実施形態において、マーカーの異常な配列及び/又は共発現レベルは、被験体におけるがんの存在を示すことができる、又は被験体におけるがん治療失敗を予測することができる。マーカーは、任意の適切ながん経路から選択することができ、好ましい実施形態においてSCARの又は「幹性」経路からのマーカーが含まれる。異常な配列検出のために、これらのマーカーは、ELF3;PCDH15;MALAT1;PTPN11;RB1;CHST6;NF1;VEZF1;TP53;SMAD4;KEAP1;STK11;PRX;ZNF28;IDH1;FEZ2;DPPA2;LPHN3;KIAA1244;EPHA7;EGFR;TLR4;DAB2IP;NOTCH1;GLUD2;DMD;KDM6A;KRAS;CDKN2A;DNMT3A;FLT3;NFE2L2;NPM1;MIR142;FOXL2;H3F3A;H3F3B;KMT2D;RNF43;TERT;ERBB2;PLCG1からなる群より選択される遺伝子であることができる。異常な発現検出のために、これらのマーカーは、PLCXD1、HKR1、ZNF283、ADA、AMACR+p63、ANK3、BCL2L1、BIRC5、BMI−1、BUB1、CCNB1、CCND1、CES1、CHAF1A、CRIP1、CRYAB、ESM1、EZH2、FGFR2、FOS、Gbx2、HCFC1、IER3、ITPR1、JUNB、KLF6、KI67、KNTC2、MGC5466、Phc1、RNF2、Suz12、TCF2、TRAP100、USP22、Wnt5A及びZFP36からなる群より選択される遺伝子であることができる。好ましい実施形態において、マーカーは、SCARの経路の制御及び下流遺伝要素、転写因子、及びメチル化パターンからなる群より選択される。一つの好ましい実施形態において検出される異常な配列、及び別の好ましい実施形態において検出される異常な共発現レベルは、SCARの経路の制御及び下流遺伝要素、転写因子、及びメチル化パターンのものである。検出されるマーカーは、mRNA、RNA、DNA、タンパク質、又はペプチドの形態である。
一実施形態において、少なくとも1つであるが好ましくは2つ以上のマーカーの異常な発現レベルは、本分野において既知の任意の検出手段によって検出することができ、検出手段には、次世代シーケンシング、多色定量的免疫蛍光共局在分析、蛍光in situハイブリダイゼーション、及び定量RT−PCR分析からなる群より選択される分析に細胞を供することが含まれるが、これに限定されない。
一態様において、本開示は、とりわけ、単細胞、細胞の集団、又はリキッドバイオプシーサンプルにおいて、少なくとも1つであるが好ましくは2つ以上のマーカーの異常な配列及び/又は共発現レベルを同時に検出するための方法を対象とする。一実施形態において、組織、細胞、又は体液の標本のサンプルを得る。一実施形態において、経路によって規定されるマーカーを選択する。一実施形態において、少なくとも1つであるが好ましくは2つ以上のマーカーの同時の異常な配列及び/又は発現レベルについてスクリーニングする。一実施形態において、それらの配列を、配列の質(配列全体又はその断片内の塩基の位置の配列)が参照配列と比較して異なる場合に、異常であるとスコア付けする。一実施形態において、検出された発現レベルがある特定の閾値を上回る場合に、それらの発現レベルを異常であるとスコア付けする。
一態様において、本開示は、とりわけ、単細胞、細胞の集団、又はリキッドバイオプシーサンプルにおいて、少なくとも1つであるが好ましくは2つ以上のマーカーの異常な配列及び/又は共発現レベルのうち少なくとも1つを検出するための方法を対象とする。一実施形態において、組織、細胞、又は体液の標本のサンプルを得る。一実施形態において、経路によって規定されるマーカーを選択する。一実施形態において、少なくとも1つであるが好ましくは2つ以上のマーカーの同時の異常な配列及び/又は発現レベルについてスクリーニングする。一実施形態において、それらの配列を、配列の質(配列全体又はその断片内の塩基の位置の配列)が参照配列と比較して異なる場合に、異常であるとスコア付けする。一実施形態において、検出された発現レベルがある特定の閾値を上回る場合に、それらの発現レベルを異常であるとスコア付けする。
一態様において、本開示は、とりわけ、単細胞、細胞の集団、又はリキッドバイオプシーサンプルにおいて、2つ以上のマーカーの異常な配列又は共発現レベルを同時に検出するのに有用なキットを対象とする。一態様において、本開示は、とりわけ、単細胞、細胞の集団、又はリキッドバイオプシーサンプルにおいて、2つ以上のマーカーの異常な配列又は共発現レベルのうち少なくとも1つを検出するのに有用なキットを対象とする。
一態様において、本開示は、とりわけ、任意の適切ながん経路から選択される分子マーカー(これは、好ましい実施形態において、SCARの又は「幹性」経路からのマーカーを含む)を保持する悪性腫瘍の標的治療の方法を対象とする。前記悪性腫瘍の治療標的化は、mRNA、RNA、DNA、タンパク質、又はペプチドの形態の検出されるマーカーによって導かれる。好ましい実施形態において、治療様式は、制御性SCAR遺伝子座及びSCARの経路の下流の遺伝要素からなる群より選択される分子標的に対して設計される。
本開示は、複数のマーカーについて配列、発現レベル、遺伝子レベル、転写レベルなどを検出することによって、がんを診断する、がん治療成績を予測する、がんを有する被験体が様々な異なるタイプの処置レジメンの影響を受けやすいかどうかを決定する、がん処置の有効性を監視する、がん診断又はがん治療失敗の予後を決定する、などのための、1又は複数の方法論又は技術を詳述する。
一実施形態において、処置不可能ながん(例えば、活性化された内在性ヒト幹細胞関連レトロウイルス(SCAR)ネットワークを有するもの)を診断するための1又は複数の方法論又は技術は、42の遺伝子(表S1に列挙される)の配列の突然変異の検出、特定のSCAR配列の転写レベルの分析、タンパク質配列のレベルの分析、シグネチャーにおける発現レベルの分析、遺伝子発現レベルの決定、及びデータセットS1、データセットS2、及びデータセットS3の遺伝子コピー数の決定のうち1又は複数を含む。
例えば、一実施形態において、方法論又は技術は、表3.3に列挙される1又は複数の遺伝子座の発現レベルに関連する異常な配列又は異常な発現レベルを示す1又は複数の入力を受け取ることに応答して、ユーザ固有のがん治療プロトコル、又はユーザ固有のがん診断を生成する工程を含む。ゲノムシグネチャー経路、シグネチャー評価方法などの非限定的な例は、米国特許第8,349,555号及び7,890,267号公報(これらのそれぞれは、その全体が参照により本明細書に組み込まれる)に見出すことができる。
一実施形態において、方法論又は技術は、補足表(Supplemental Table)S3に列挙される1又は複数のペプチドの発現レベルに関連する異常な発現レベルを示す1又は複数の入力を受け取ることに応答して、疾患状態(例えば、がん治療有効性状態、がん治療進行、がん予後、がん診断、治療失敗、再燃、再発など)の予測因子を生成する工程を含む。
一実施形態において、方法論又は技術は、補足表S4に列挙される遺伝子のSCARの経路活性化シグネチャーを示す1又は複数の入力を受け取ることに応答して、疾患状態(例えば、がん治療有効性状態、がん治療進行、がん予後、がん診断、治療失敗、再燃、再発など)の予測因子を生成する工程を含む。
一実施形態において、方法論又は技術は、補足表S5に列挙される1又は複数の遺伝子座の発現レベルに関連する異常な発現レベルを示す1又は複数の入力を受け取ることに応答して、SCAR活性化状態を生成する工程を含む。
一実施形態において、方法論又は技術は、補足表S6に列挙される1又は複数の遺伝子座の発現レベルに関連する異常な発現レベルを示す1又は複数の入力を受け取ることに応答して、疾患状態(例えば、がん治療有効性状態、がん治療進行、がん予後、がん診断、治療失敗、再燃、再発など)の予測因子を生成する工程を含む。
一実施形態において、方法論又は技術は、データセットS1に列挙される1又は複数の遺伝子座の発現レベル又はコピー数に関連する異常な発現レベル又は遺伝子コピー数を示す1又は複数の入力を受け取ることに応答して、疾患状態(例えば、がん治療有効性状態、がん治療進行、がん予後、がん診断、治療失敗、再燃、再発など)の予測因子を生成する工程を含む。
一実施形態において、方法論又は技術は、データセットS2に列挙される1又は複数の配列の発現レベルに関連する異常な発現レベルを示す1又は複数の入力を受け取ることに応答して、疾患状態(例えば、がん治療有効性状態、がん治療進行、がん予後、がん診断、治療失敗、再燃、再発など)の予測因子を生成する工程を含む。
一態様において、本開示は、とりわけ、SCAR配列に由来するゲノム遺伝子座によってコードされる共通のペプチド配列の同定の方法を対象とする。一実施形態において、本方法は、複数の異なるゲノム座標に位置するSCAR由来のゲノム遺伝子座の核酸配列を検索する工程;及び前記核酸配列内の全てのオープンリーディングフレーム(ORF)を同定する工程を含む。一実施形態において、本方法は、前記核酸配列によってコードされ、前記核酸配列から潜在的に転写される、全てのペプチド配列を同定する工程;及び複数の異なるゲノム座標に位置する複数の異なるSCAR由来のゲノム遺伝子座に共通のペプチド配列を同定する工程を、さらに含む。
一実施形態において、方法論又は技術は、遺伝的シグネチャー情報及びタンパク質シグネチャー情報を使用してSCARのネットワーク活性化を決定する工程を含む。一実施形態において、SCARのネットワーク活性化情報は、がん成績予後を生成するために使用される。例えば、活性化されたSCARのネットワークは、不良ながん治療成績又は予後不良を示す。
一実施形態において、方法論又は技術は、SCARネットワークマーカーの発現レベルを示す1又は複数の入力及び異常な配列を示す1又は複数の入力に基づいてがん関連成績を生成する工程を含む。
SCARのネットワークの非限定的な例は、i)対応するRNA分子の発現の検出に基づいて規定される、転写活性のあるSCARの遺伝子座;及びii)タンパク質コード遺伝子、非コードRNA分子、マイクロRNA、及びSCAR活性によって影響を受ける他の制御及び構造分子をコードする遺伝子を含む、下流のSCAR制御コード遺伝子の発現シグネチャーの、ゲノムワイドな概略を含む。
SCAR経路の非限定的な例は、単細胞及び細胞の集団を含む、特定の細胞及び/又は特定の生体サンプルにおいて転写活性のあるSCARの遺伝子座のサブセットを含む。
SCARの経路:単細胞及び細胞の集団を含む、特定の細胞及び/又は特定の生体サンプルにおけるゲノムワイドなSCARのネットワーク分析によって規定されるゲノム遺伝子座のサブセット。
シグネチャーの非限定的な例は、74遺伝子シグネチャー(例えば表S4を参照)、55遺伝子シグネチャー(例えば表S4を参照)、ヒト卵母細胞の単細胞分析によって規定されるSCARの経路シグネチャーを含み、これらの遺伝子の発現変化はHERV−H由来レトロウイルス配列の活性化された転写に関連しているようである。遺伝子記号は第1列に列挙されている。これらはコード遺伝子であり、その発現は、HERV−Hにコードされた制御転写物を標的化するshRNA干渉プロトコルを使用して、特定の方法で変化する(上方制御及び下方制御される)(対数変換された発現変化倍率は第2列に列挙されている)。ヒト卵母細胞におけるこれらの遺伝子の発現変化(対数変換された発現変化倍率は第3列に列挙されている)は、HERV−H経路活性化と一致する(r=−0.74043)、すなわち、発現がshHERVH干渉に続いて上方制御される遺伝子は、卵母細胞において下方制御されるようであり、反対に、発現がshHERVH干渉に続いて下方制御される遺伝子は、卵母細胞において上方制御されるようである。これらのシグネチャーの有用性は、前立腺がんサンプル及び前立腺上皮内新生物サンプルを含む、正常ヒト前立腺及び病的ヒト前立腺のサンプルの分析によって証明された(図1C及び2D)。表S4に列挙される個々の遺伝子のそれぞれの発現変化倍率は、当業者に既知の技術及び方法を使用して決定されるだろう。対応する遺伝子についての値は、それが第3列に列挙される卵母細胞の値について示されているとおりに、表S4において規定される順序で列挙されることとなる。次に、第2列及び第3列に列挙される値について、相関係数を計算する。相関係数の負の値は、SCARの経路活性化を示すものとして解釈されるべきものである。相関係数の正の値は、SCARの経路活性化の証拠がないことを示すものである。
一実施形態において、遺伝的シグネチャー及びタンパク質シグネチャーは、独立して疾患状態の予測因子として使用される。一実施形態において、現在のシグネチャーに列挙されるいくつかの特定の遺伝子/タンパク質標的は、がんに関連する可能性が高い。一実施形態において、現在のシグネチャーに列挙されるいくつかの特定の遺伝子/タンパク質標的は、SCARの経路及びネットワーク活性化を検出するために利用される。
図1A−1K(集合的に「図1」と呼ばれ、図1A−1Dは時折「A」と呼ばれ、図1E−1Hは時折「B」と呼ばれ、図1I−1Kは時折「C」と呼ばれる)。1細胞期対8細胞期の正倍数体及び異数体ヒト胚(図1A−1D)、発生的に生存可能な接合子対非生存可能な接合子(A、図1D)、及びインビボの成熟したヒト卵母細胞(図1E−1H)におけるHERVH制御遺伝子の異なる発現パターン。
A(図1A−1D):ヒト接合子対8細胞ヒト胚において異なって発現される合計36の統計的に有意な遺伝子は、hESCにおいてHERVH/LBP9によって制御される。これらの遺伝子のうち14個の発現は、正倍数体ヒト胚対異数体ヒト胚において顕著に異なり(図1A及び1C)、一方でこれらの遺伝子のうち22個の発現は、正倍数体ヒト胚対異数体ヒト胚において顕著に異ならない(図1B)。同様に、174のHERVH制御遺伝子の発現シグネチャーは、発生的に生存可能なヒト接合子及び非生存可能なヒト接合子において異なる(q<0.0005;A、図1D)。発生的に非生存可能なヒト接合子において上方制御された遺伝子がハイライトされている。
B(図1E−1H):マイクロアレイ分析は、成熟ヒト卵母細胞におけるHERVH制御遺伝子の遺伝子発現シグネチャーを同定する。
図2A−2M(集合的に「図2」と呼ばれ、図2A−2Dは時折「A」と呼ばれ、図2E−2Hは時折「B」と呼ばれ、図2I及び2Jは時折「C」と呼ばれ、図2K−2Mは時折「D」と呼ばれる)。ヒト着床前胚発生の様々な段階での(A−C)、及び正常前立腺上皮、正常前立腺間質、良性前立腺肥大、前立腺の萎縮性病変、前立腺上皮内新生物(PIN)の推定上の前立腺がん前駆病変、前立腺がん病変に隣接する形態学的に正常な前立腺上皮、限局性前立腺がん、及び転移性前立腺がんの臨床サンプルにおける(D)、個々のSCAR遺伝子座(A、B)、LncRNA HPAT3のSCAR−制御配列(C)、及びSCAR制御タンパク質コード遺伝子(D)の単細胞次世代シーケンシング(A−C)及びマイクロアレイ遺伝子発現解析(D)。
A−C、ヒト着床前胚の単細胞次世代RNA配列解析は、ヒト卵母細胞及び接合子における選択されたHERVH及びpHERVK遺伝子座の発現の活性化を明らかにする。ヒト着床前胚の各段階での個々のHERV遺伝子座の発現パターンが示されている。プロットされた発現値は、卵母細胞における発現レベルに対して正規化された平均発現値(A)又は胚発生の対応する段階の全ての個々の細胞における実際の測定値(B、C)のいずれかによって規定した。
D、正常前立腺上皮、正常前立腺間質、良性前立腺肥大、前立腺の萎縮性病変、前立腺上皮内新生物(PIN)の推定上の前立腺がん前駆病変、前立腺がん病変に隣接する形態学的に正常な前立腺上皮、限局性前立腺がん、及び転移性前立腺がんから切除された細胞で構成される、正常前立腺上皮から転移性前立腺腫瘍への悪性進行の仮想的順序の重要な段階を表す、臨床サンプルのマイクロアレイ遺伝子発現プロファイリング。
図3A−3D(集合的に「図3」と呼ばれ、図3A−3Dは以下の説明においてそれぞれA−Dとも呼ばれる)。SCAR標的化タンパク質コード遺伝子の遺伝子発現及び遺伝子コピー数の変化は、がん患者の長期生存率との顕著な関連を示す。宿主/ウイルスキメラ転写物の構造成分を含むタンパク質コード遺伝子の遺伝子コピー数及びmRNA発現レベルを、12のTCGAコホート(12の異なるがんタイプのPANCAN12研究)にわたる5,158の臨床サンプル及び全てのTCGAコホートにわたる12,093の臨床サンプルを含むTCGA汎がんデータベースにおけるカプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率との関連について評価した。遺伝子発現変化(A−C)及び遺伝子コピー数変化(D)の、TCGA PANCAN12研究のがん患者の長期生存率との顕著な関連を示すSCAR標的化遺伝子の例を示す(A;C;D)。がんの3つの個別のタイプ[前立腺がん(n=568)、乳がん(n=1,241)、及び直腸がん(n=187)]のTCGAコホートについてのこれらの関連の代表的な例を(B)に示す。遺伝子発現ヒートマップ及び対応するカプラン・マイヤー生存曲線を(A)に示す。遺伝子発現(左の画像)及びコピー数(右の画像)のヒートマップ及び関連するカプラン・マイヤー生存曲線を(D)に示す。縦の破線は、10年生存データ点を示す。対応するp値を補足データセットS1において報告する。
図4A−4D(集合的に「図4」と呼ばれ、図4A−4Dは以下の説明においてそれぞれA−Dとも呼ばれる)。ヒト特異的ウイルス/宿主キメラ転写物の翻訳されたアミノ酸配列のタンパク質アラインメントは、異なるSCAR遺伝子座によってコードされる保存タンパク質ドメインの異なるパターンを同定する。材料及び方法において説明されているように、ヒト特異的キメラ転写物のヌクレオチド配列をアミノ酸配列に翻訳し、BLASTタンパク質アラインメント分析に供した。ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内で最も頻繁に表される保存タンパク質ドメインは、GVQWアミノ酸配列であることに留意されたい(A;C;D)。
図5A−5D(集合的に「図5」と呼ばれ、図5A−5Dは以下の説明においてそれぞれA−Dとも呼ばれる)。chrX:278899-284216及びchrY:278899-284216のSCARのヒト特異的キメラウイルス/宿主転写物の同一の核酸配列に由来するGVQW保存タンパク質ドメインのヒト特異的拡大の進化的追跡。GVQW保存ドメインをコードするヌクレオチド配列を、数個の隣接するアミノ酸を含むように拡大した。これは、SCARの遺伝子座特異的ヌクレオチド配列を得るのに十分であった。GVQWをコードする配列のゲノム起源を、所与のゲノム配列及び対応する遺伝子座特異的SCAR由来の配列の100%ヌクレオチド配列同一性に基づいて、推測した。BLATアルゴリズムを利用して、特定のSCARの遺伝子座によってコードされるキメラウイルス/宿主転写物の配列と100%同一である、ヒト及び非ヒト霊長類のゲノム中のGVQWをコードするヌクレオチド配列の数を決定した。マウス及びラットゲノムにおいて、GVQW保存タンパク質ドメインをコードする配列は検出されなかったことに留意されたい。chrX:278899-284216及び/又はchrY:278899-284216のSCAR転写物に由来するGVQWをコードする配列のみが、ヒトゲノムにおいて顕著に拡大しているようであり(図3Cの赤のバー)、この拡大は、他の大型類人猿と比較したときのヒトプロテオームにおける、保存されたGVQWドメインを保持するタンパク質の数の著しい富化に関連している(図3D)。
図6A−6B(集合的に「図6」と呼ばれ、図6A−6Bは以下の説明においてそれぞれA−Bとも呼ばれる)。GVQW保存タンパク質ドメインを保持する21のジンクフィンガータンパク質の遺伝子レベルのコピー数の変化は、29の異なるタイプの悪性腫瘍を有すると診断されたがん患者の長期生存率との顕著な関連を示す。29のがんタイプを表す全てのTCGAコホートにわたる12,093の臨床サンプルを含むTCGA汎がんデータベースのカプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率との関連について、GVQW保存タンパク質ドメインを保持する現在同定された全てのジンクフィンガータンパク質の遺伝子コピー数を評価した。遺伝子コピー数変化(A)及び関連するカプラン・マイヤー生存曲線(B)のヒートマップを示す。カプラン・マイヤー生存率分析の結果が、GVQW保存タンパク質ドメイン及び3のSCAR標的化ジンクフィンガータンパク質(ZNF443;ZNF587;ZNF814)を保持する21のジンクフィンガータンパク質について示されている。報告されたp値は、Xenaがんゲノムブラウザデータ可視化ツール(http://xena.ucsc.edu/)によって生成したカプラン・マイヤー生存曲線からのものである。
図7A−7D(集合的に「図7」と呼ばれ、図7A−7Dは以下の説明においてそれぞれA−Dとも呼ばれる)。低い生存率及び高いがんによる死亡の可能性によって規定される、悪性腫瘍の臨床的難治性の体細胞非サイレント突然変異のシグネチャー。
A、がんによる死亡の表現型の体細胞非サイレント突然変異シグネチャーを保持する18の遺伝子の同定。18のトップスコアのヒト遺伝子を同定した。ここで、体細胞非サイレント突然変異(SNM)の最大の数が、全てのTCGAコホートにわたる12,093の腫瘍サンプルにおいて検出された(ただし、腫瘍におけるこれらの突然変異の存在が、カプラン・マイヤー生存率分析によって規定される顕著に高いがんによる死亡の可能性に関連している、という要件が満たされている)。上のパネルは、TP53遺伝子のSNMのプロファイルに整列させた患者の腫瘍サンプルの中の18の遺伝子のSNMの分布を示す。18の遺伝子のそれぞれのSNMを有するがん患者の数は、事象のパーセントとして報告される。陰影の付いたエリアは、SNMのないがん患者の相対的な数をハイライトしている。これらの18の遺伝子のそれぞれについてのカプラン・マイヤー生存曲線が、生存確率が顕著に低くがんによる死亡の可能性が高い患者を同定することに留意されたい。それ故に、腫瘍サンプルから単離されたこれらの18の遺伝子のそれぞれにおけるSNMの検出は、腫瘍にこれらの遺伝子のSNMを有しない患者と比較して、がん患者の不良な長期予後に関連している(図5A)。下線を付した遺伝子記号は、その発現が、hESCにおいてSCARによって制御される遺伝子を同定している。赤色の遺伝子記号はSCAR標的化遺伝子を表し、一方で黒色の遺伝子記号は既報の候補がんドライバー遺伝子を同定する。
B、TP53遺伝子のみについての、腫瘍にSNMを有する及び有しない7,509人のがん患者のカプラン・マイヤー生存率分析の比較(図7A、下記の左上の図);18遺伝子のSNMのシグネチャー(図7B、下記の右上の図);TP53のない26遺伝子のSNMのシグネチャー(図7C、下記の左下の図);TP53遺伝子を含む27遺伝子のSNMのシグネチャー(図7D、下記の右下の図)。
C、D、28(C)及び19(D)のがんタイプを有すると診断された患者における、悪性腫瘍の臨床的難治性の線形回帰分析。C、28のがんタイプのTCGA汎がんコホートからのがん患者の生存データを利用して、各がんタイプについての死亡事象のパーセントを計算した。得られた値を、がん患者の対応する群におけるSNMのがんによる死亡シグネチャーを有する患者のパーセントと整列させて、線形回帰分析に供した。D、19がんタイプについての米国における年齢調整されたがん発生率及び死亡率(100,000人あたり)を、疾病管理予防センター(CDC)の米国がん統計(USCS)の報告書から得た。各がんタイプについての推定死亡率は、対応する発生率の値及びSNMのがんによる死亡シグネチャーを有する患者のパーセントの値を掛け合わせることによって計算した。得られた値を対応するがんタイプについての実際の死亡率と整列させて、回帰分析に供した。
図8A−8B(集合的に「図8」と呼ばれ、図8A−8Bは以下の説明においてそれぞれA−Bとも呼ばれる)。SCAR幹性ネットワークの遺伝子のタンパク質発現変化は、低い長期生存率及び高いがんによる死亡の可能性と統計的に有意な関連を示す。
12のTCGAコホートにわたる5,158の臨床サンプルを含むTCGA汎がんデータベースにおける、カプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率との関連について、38のSCAR幹性ネットワークの遺伝子のタンパク質発現変化を評価した。全体として、23のSCAR制御遺伝子(60.5%)のタンパク質発現レベルの変化が、がん患者の長期生存確率と顕著な関連を示した(補足データセットS1)。タンパク質発現のヒートマップ及び関連するカプラン・マイヤー生存曲線を示す。対応するp値を補足データセットS1に報告する。
図9。転写活性のあるLTR7/HERVH SCARは、代替の非相同末端結合(NHEJ)DNA修復経路を利用することによって、宿主DNA(青線)の二本鎖切断(雷の形)の修復に寄与する。部分的な相同性領域を有するSCAR RNA(黒破線)の宿主DNAへの逆転写は、宿主DNAの二本鎖切断の部位でギャップを埋めるDNA分子(黒実線)を生成する。二本鎖DNA切断のSCAR関連修復のこのメカニズムの特徴は、ヒトゲノム中のLTR7/HERVH配列の挿入部位での祖先DNAセグメント(赤実線)の欠失の証拠である(詳細については表3及びテキストを参照されたい)。このプロセスは、SCARのヒト特異的組込み部位を生成し、宿主/ウイルスキメラ転写物(青/黒破線)の生成を容易にする可能性がある。DSB、二本鎖切断;NHEJ、非相同末端結合;RT、逆転写;SCAR、幹細胞関連レトロウイルス。
図10。血液、血清、及び血漿サンプル;循環腫瘍細胞;原発性及び転移性腫瘍サンプルにおける、SCARのネットワーク活性状態を監視するための、継続的な連続サンプリングに基づく、がん患者の臨床管理における意思決定プロセスのフローチャート。
この配列の任意の段階で活性化されたSCARのネットワークの遺伝的及び/又は分子的証拠の同定は、治療抵抗性の臨床的に致死的な疾患の表現型の診断に有利に働き、以下の治療選択肢の即座の検討のために必要なものをもたらすであろう:「次の候補の」積極的な処置プロトコル;SCARの経路を特異的に標的化する新規の治療及び/又はSCARのネットワークの活性状態を示す悪性腫瘍を有する患者に適切だと考えられる治療的介入。CTC、循環腫瘍細胞;FFPE、ホルマリン固定パラフィン包埋。Glinsky, GV. 2008. "Stemness" genomics law governs clinical behavior of human cancer: Implications for decision making in disease management. Journal of Clinical Oncology, 26: 2846-53から採用した。
補足図(Supplemental Figure)S1A−S1K(集合的に「補足図S1」と呼ばれる)。(図4に関連する)。hESCの内在性ヒトSCARによってコードされた宿主/ウイルスキメラ転写物の翻訳されたアミノ酸配列内の、保存タンパク質ドメイン発現の異なる及び共通の複数のパターンの追加の例。タンパク質BLASTアルゴリズム(http://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&BLAST_PROGRAMS=blastp&PAGE_TYPE=BlastSearch&SHOW_DEFAULTS=on&LINK_LOC=blasthome)及び保存タンパク質ドメインの同定及び可視化のための関連のウェブベースツール(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi?RID=3HZ5BMES01R&mode=all)(これらは、他の場所で詳細に説明された[80、81])を使用して、ヒト特異的キメラ転写物のヌクレオチド配列をアミノ酸配列に翻訳し、タンパク質アラインメント分析に供した。
ヒト特異的ウイルス/宿主キメラ転写物の翻訳されたアミノ酸配列のタンパク質アラインメントは、異なるSCAR遺伝子座によってコードされる保存タンパク質ドメインの異なるパターンを同定する。材料及び方法において説明されているように、ヒト特異的キメラ転写物のヌクレオチド配列をアミノ酸配列に翻訳し、BLASTタンパク質アラインメント分析に供した。ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内で最も頻繁に表される保存タンパク質ドメインは、GVQWアミノ酸配列であることに留意されたい。
補足図S2A−S2E(集合的に「補足図S2」と呼ばれ、補足図S2A−S2Dは、以下の説明においてそれぞれA−Dとも呼ばれ、補足図S2Eは補足図S2A−S2Eの編集物である)。(図6に関連する)。GVQW保存タンパク質ドメインを保持するジンクフィンガータンパク質の遺伝子発現及び遺伝子コピー数の変化は、がん患者の長期生存率との顕著な関連を示す。GVQW保存タンパク質ドメインを保持するジンクフィンガータンパク質の遺伝子コピー数(D)及びmRNA発現レベル(A−C)を、前立腺がん(n=568);乳がん(n=1,241);結腸がん(n=550);直腸がん(n=187);膵臓がん(n=196)を有すると診断されたがん患者;及び12のTCGAコホート(12の異なるがんタイプのPANCAN12研究)にわたる5,158の臨床サンプルを含むTCGA汎がんデータベースの、カプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率との関連について、評価した。がんの5つの個別のタイプ[前立腺がん(A);乳がん(B;C、左下のパネル);結腸がん(C;左上のパネル);直腸がん(C;右上のパネル);及び膵臓がん(C、右下のパネル)]のTCGAコホートにおける、遺伝子発現変化(A−C)の顕著な関連を示す、GVQW保存タンパク質ドメインを有するジンクフィンガータンパク質の代表的な例を示す。遺伝子コピー数変化のTCGA PANCAN12研究のがん患者の長期生存率との顕著な関連を示す、GVQW保存タンパク質ドメインを有するジンクフィンガータンパク質の例を、図4Dに示す。遺伝子発現ヒートマップ及び対応するカプラン・マイヤー生存曲線を(A−C)に示す。遺伝子発現(左の画像)及びエキソン発現(右の画像)のヒートマップ及び関連するカプラン・マイヤー生存曲線を(C)に示す。遺伝子発現(左の画像)及びコピー数(右の画像)のヒートマップ及び関連するカプラン・マイヤー生存曲線を(D)に示す。対応するp値を補足データセットS1に報告する。
補足図S3A−S3B(集合的に「補足図S3」と呼ばれ、補足図S3A−S3Bは、以下の説明においてそれぞれA−Bとも呼ばれる)。(図7に関連する)。フォローアップ生存データの完全な臨床記録を有する患者のみを含む、SNM遺伝子の分類性能の追加のカプラン・マイヤー生存率分析。
A、腫瘍にSNMを有する及び有しない7,258人のがん患者(左上及び左下の図)及び生存期間の昇順にソートした後同一のサイズ(n=2,419)のサブグループに階層化したがん患者(左上及び左下の図)のカプラン・マイヤー生存率分析の比較。この分析において、フォローアップ生存データの完全な臨床記録を有する患者のみを含めた。
B、がんによる死亡の表現型のSNMシグネチャーを保持する遺伝子の突然変異のフィンガープリントの可視化。臨床腫瘍サンプルから単離されたこれらの遺伝子には、突然変異が「散らばっている」ように見え、その圧倒的多数がSNMによって表されることに留意されたい。
補足図S4A−S4E(集合的に「補足図S4」と呼ばれ、補足図S4A−S4Dは、以下の説明においてそれぞれA−Dとも呼ばれ、補足図S4Eは補足図S4A−S4Eの編集物である)。HESCにおけるSCAR関連幹性ネットワークのマスター転写制御因子の遺伝子レベルのコピー数の変化(枠で囲んだKLF4;LBP9;NANOG;及びPOU5F1遺伝子のカプラン・マイヤープロット)及びSNMのがんによる死亡シグネチャーの遺伝子は、低い長期生存率及び高いがんによる死亡の可能性と統計的に有意な関連を示す。示されたタンパク質コード遺伝子の遺伝子レベルのコピー数変化を、12のTCGAコホートにわたる5,158の臨床サンプル(A;C)及び29のTCGAコホートにわたる12,093の臨床サンプル(B;D)を含む2つのTCGA汎がんデータベースにおいて、カプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率との関連について、独立して評価した。著しく類似した結果が、BMI1(C及びDの左下の図)及びEZH2(C及びDの右下の図)遺伝子のコピー数変化について観察され、その遺伝子のコピー数変化と、がんによる死亡の可能性が高いがん患者からの腫瘍における、Polycombクロマチンサイレンシング経路の活性化及び幹性遺伝子発現シグネチャーとの関連は、既に実証されたものである(37−51)ことに留意されたい。対応するp値を補足データセットS1に報告する。
補足図S5。11遺伝子のがんによる死亡シグネチャー及び3つのSCARネットワーク遺伝子(PLCXD1、HKR1、ZNF283)の発現シグネチャーの発現プロファイルに基づいて、異なるサブグループに階層化された前立腺がん患者における、治療成績のカプラン・マイヤー生存率分析。
表S1A(「表S1」とも呼ばれる)。複数の汎がんデータベースにおける治療失敗及びがんによる死亡の高い可能性との顕著な関連に基づいて同定された、体細胞非サイレント突然変異の分析のための42の遺伝子のパネル。
表S2A−S2C(集合的に「表S2」と呼ばれる)。
表S2A:両側p値:0.00090474;p=0.0009;図7Cに関連する。
表S2B:両側p値;図7Dに関連する。
表S2C:図7に関連する。
表S3A−S3B(集合的に「表S3」と呼ばれる)。
表S3A。ChrY_ChrX
表S3B。chr3_chr11
表S4A−S4B(集合的に「表S4」と呼ばれる)。
表S4A。74の遺伝子。
表S4B。55の遺伝子。
表S5A−S5C(集合的に「表S5」と呼ばれる)。
表S5A。ヒト胚形成の接合子段階で最も顕著な活性化を示すHERVH−遺伝子座。図2に関連する。
表S5B。ヒト胚形成の接合子段階で最も顕著な活性化を示す、HERVK−;HERVH−;及び他のSCAR遺伝子座。図2に関連する。
表S5C(部分的に図において示される)。ヒト胚形成及びヒト被験体における病理学的状態の発症に関係するSCAR配列。
表S6A−S6C(集合的に「表S6」と呼ばれる)。
表S6A。64のHERV1ヒト特異的キメラ転写物(ボノボ及びチンパンジーアラインメント失敗)。
表S6B。
表S6C。
詳細な説明
近年、新規の方法の、個人に合わせた、標的の調整されたがん治療を含む、幅広いがん処置プロトコルが開発されている。しばしば、非常に積極的ながん治療は、そのような治療によって生じる望ましくない副作用に起因して、末期がんのために保留される。しかしながら、そのような積極的な治療でさえ、そのような末期では一般に失敗する。最も積極的な治療にのみ応答するがんをより初期段階で同定できれば、そのようながんを有する患者の予後を大幅に向上できる可能性がある。
近年、そのような成績を予測する潜在的に有用なマーカーが同定されている。Glinsky, G.V. et al., J. Clin. Invest. 113: 913-923 (2004)は、遺伝子発現プロファイリングが前立腺がんの臨床成績を予測することを教示している。Van 't Veer et al., Nature 415: 530-536 (2002)は、遺伝子発現プロファイリングが乳がんの臨床成績を予測することを教示している。Glinsky et al., J. Clin. Invest. 115: 1503-1521 (2005)は、BMI1癌遺伝子の発現の変化が、正常及び白血病幹細胞の自己再生状態、並びに複数のタイプのがんを有する患者における治療失敗を予測する11遺伝子のがんによる死亡のシグネチャーの予後不良プロファイルに、機能的にリンクしていることを教示している。これらの研究は、マイクロアレイ遺伝子発現解析アプローチを利用した。
それ故に、がんの早期診断のための高度に正確な方法及び臨床現場に容易に適合可能ながん治療のための予後アッセイが、継続してますます必要とされている。そのような方法は、臨床検査室において容易に実施することができる状態の最新技術を利用するはずであり、標準的な治療レジメンに適用される様々ながんの抵抗性の可能性を正確に予測するはずである。
処置を発見、規定、及び設計し、処置を開発するために、そして転移性及び難治性のがんを処置するために、非常にたくさんの試みがなされている。それらは主に、迅速な細胞成長の基本的メカニズム又は異常ながん細胞代謝経路のいずれかを攻撃することによるが、ほとんど成功していない。最近、腫瘍及び微小転移に対するその攻撃において、免疫システムを有効化又は再有効化するいくつかの方法が、試験及び商業的使用においてかなり有望なデータを示しているが、転移性及び難治性の疾患を有する患者の大多数が、これらの免疫調節治療に対しても難治性であることが証明されている。それ故に、単独の治療剤として、又は他の様式、特に免疫調節と組み合わせて使用され、転移表現型を可能にする細胞メカニズムを根本的に攻撃するように設計された、新規のがん治療が必要とされている。そのような新規の治療は、がん細胞の転移及び生存に関与している重要遺伝子シグネチャーの理解から誘導されるはずである。
体細胞突然変異及び染色体不安定性は、がん細胞におけるゲノム異常の特徴である。異数体は、ヒト胚及び悪性固形腫瘍においてしばしば観察される、染色体不安定性の一般的な発現を表す。ヒト内在性レトロウイルス(HERV)由来遺伝子座の活性化は、着床前ヒト胚、hESC、及び複数のタイプのヒト悪性腫瘍において実証されている。HERV活性化が、ヒト胚発生の初期段階における染色体不安定性の頻繁な発生及びがんにおけるゲノム異常の出現に寄与する、共通の分子経路を浮き彫りにすることができるかどうかは不明のままである。
ヒト着床前胚の単細胞RNA配列解析は、卵母細胞から接合子への移行中の特定のLTR7/HERVH遺伝子座の活性化を明らかにし、ヒト胚における異数性に関連するHERVHネットワークシグネチャーを同定する。相関パターンの分析は、推定上の前駆病変(前立腺上皮内新生物)から限局性及び転移性の前立腺がんへのがん進行経路の存在を支持する、前立腺腫瘍の臨床サンプルの遺伝子発現プロファイルを有するインビボの成熟したヒト卵母細胞のHERVHネットワーク活性化のトランスクリプトームシグネチャーとリンクしている。既知の長期治療成績を有するがん患者の腫瘍サンプルにおいて、HERVHネットワークの活性化のシグネチャーを追跡することにより、治療失敗及びがんによる死亡の著しく異なる可能性を有するサブグループへの患者の階層化が可能になった。
29のがんタイプにわたる12,093の臨床腫瘍サンプルにおける、幹細胞関連レトロウイルス(SCAR)によって制御されるネットワークのヒト特異的遺伝要素のゲノムワイドな分析は、SCAR制御宿主遺伝子の転写物の及びタンパク質の発現、体細胞非サイレント突然変異(SNM)の存在、遺伝子レベルのコピー数変化によって規定される、臨床的に致死的な治療抵抗性疾患の汎がんゲノムシグネチャーを明らかにした。全てのがん死亡の73%超が、腫瘍にSNMのシグネチャーのある患者において発生した。米国国民集団におけるがん難治性の線形回帰分析は、臓器特異的がん死亡率が、腫瘍にSNMのシグネチャーのある患者のパーセンテージと直接相関していることを証明した。
SCARにコードされたRNA分子は、保存タンパク質ドメイン(CPD)として規定されるアミノ酸配列を含む、内在性タンパク質コード能力を保有する。SCARにコードされたCPDのマッピングは、ゲノムワイドに散在する、CPDの数千の遺伝子座特異的フィンガープリントを明らかにした。特定のCPDをコードするSCARの配列の進化的拡大は、CPDが見出される独特のタンパク質配列のヒトプロテオームにおける著しい富化をもたらした。これらの結果は、高い発現レベルのSCAR RNAを有する罹患細胞が、免疫療法介入のための魅力的で高度に特異的な分子標的を提供する、著しく多くの量のSCAR RNAにコードされたペプチドを有する可能性が高いことを示している。
ヒト特異的ウイルス/宿主キメラ転写物の分子構造の系統的解析は、ヒトゲノムにおけるSCARの組込みの特徴が、祖先DNAの複数種生物の欠失パターンであることを証明している。ヒト特異的挿入及び欠失を有するSCARの遺伝子座の異種間の追跡は、宿主細胞の即座の生存及び適応を増強し得るSCARの推定上の生物学的機能を浮き彫りにする、二本鎖DNA切断の修復における潜在的役割を示唆している。進化的規模では、数千のヒト特異的制御配列を蒔くことに加えて、SCARの活性は、DNA修復及びヒトゲノム全体にわたる特定のCPDの配列の拡散に関与しているようである。
本明細書に示す実施例は、分化細胞におけるSCAR制御幹性ネットワークの覚醒が、複数のタイプの臨床的に致死的な悪性腫瘍においてその後容易に検出可能であり、治療抵抗性表現型の出現に寄与している可能性が高い、多様な範囲にわたるゲノム異常の発生に関連していることを証明している。
キーワード:ヒト内在性幹細胞関連レトロウイルス(SCAR);ヒト特異的制御配列;ヒトESC;ヒト胚;多能性状態制御因子;NANOG;POU5F1(OCT4);CTCF;LTR7 RNA;長い末端反復、LTR;LTR7/HERVH;LTR5HS/HERVK;治療抵抗性がん;がん幹細胞
略語リスト
HERV、ヒト内在性レトロウイルス
hESC、ヒト胚性幹細胞
LINE、長い散在反復配列(long interspersed nuclear element)
lncRNA、長い非コードRNA
lincRNA、長い遺伝子間非コードRNA
LTR、長い末端反復
NANOG、Nanog ホメオボックス
POU5F1、POUクラス5ホメオボックス1
SCAR、幹細胞関連レトロウイルス
TCGA、がんゲノムアトラス(The Cancer Genome Atlas)
TE、転位因子
TF、転写因子
TFBS、転写因子結合部位
sncRNA、小さい非コードRNA
幹細胞関連レトロウイルス(SCAR)
内在性レトロウイルスの活性は、ヒト細胞において抑制されて、機能的ゲノム完全性に対する突然変異の潜在的に有害な効果を制限し、且つゲノム安定性の維持を保証する。これに関して、ヒト胚性幹細胞(hESC)及び初期段階のヒト胚は、著しく異なっているようである。ヒト内在性レトロウイルス(HERV)、具体的には、HERVH及びHERVKサブファミリーの発現は、hESCにおいて著しく活性化される[1−3]。ヒトゲノムにおけるLTR7/HERVH配列の挿入の割合の増強は、ヒト特異的転写結合部位[3]、及び長い非コードRNA[5]を含む、多能性コア転写因子の結合部位に関連しているようである[1;3;4]。hESCにおける転写因子結合部位の分析は、HERVHの発現が多能性制御回路によって制御されることを示唆している。なぜなら、50の最も高度に発現されたHERVH遺伝子座の長い末端反復(LTR)の80%が、NANOG及びPOU5F1を含む多能性コア転写因子によって占められているからである[1]。さらに、転位因子(TE)由来配列、特にLTR7/HERVH、LTR5_Hs/HERVK、及びL1HSは、hESCのゲノムに推定上の転写因子結合部位(TFBS)を有する候補ヒト特異的制御配列(HSRS)の99.8%を保持している[3]。HERVのこれらの特定のファミリーの一般的な機能的特徴(これらは、ヒト胚及びhESCにおいて、それらの活性発現によって媒介される)に基づいて[6−9]、それらは、内在性ヒト幹細胞関連レトロウイルス(SCAR)として指定された。
最近の研究は、ヒト着床前胚及び胚性幹細胞におけるSCARの活性化のメカニズム及び推定上の生物学的機能を浮き彫りにした。LTR7/HERVHサブファミリーは、ヒト胚の胚盤胞において迅速に脱メチル化及び上方制御され、hESCにおいて高度に発現された状態を維持する[10]。下流の完全長HERVH−int要素のためのプロモーターを典型的に保持する、LTR7、LTR7B、及びLTR7Yの配列が、最も高いレベルで発現されたことが見出され、ヒトESC及び誘導多能性幹細胞、iPSCにおいて、最も統計的に有意に上方制御されたレトロトランスポゾンであった[11]。HERVHサブファミリーのLTR、具体的には、LTR7がhESCにおいてエンハンサーとして機能し、HERVH配列がhESCの多能性及び同一性の維持に必要とされる核非コードRNAをコードしていることが証明されている[12]。分化したヒト細胞の誘導多能性幹細胞(iPSC)への再プログラミング、多能性の維持及び分化能力の再確立には、HERVHの一時的な空間時間的に制御された過活性化が必要である[13]。LTR7/HERVH活性の制御及びサイレンシングに失敗すると、神経系列において分化に欠陥のある表現型が生じる[13、14]。L1レトロトランスポゾンの活性化はまた、これらのプロセスに寄与し得る。なぜならば、ヒト及び他の大型類人猿のiPSCにおいて、L1転写及び転移の両方の顕著な活性が最近報告されたためである[15]。ヒト着床前胚及び胚性幹細胞の単細胞RNAシーケンシング[16、17]により、特定のレトロウイルス要素の著しい活性化によって規定される、初期のヒト胚性幹細胞の特定の異なる集団の同定が可能となった[18]。
内在性ヒトSCARの発見及びヒト胚形成におけるそれらの必須の役割の説得力のある証拠は、いくつかの即座の実用的意義を有する可能性がある。ヒトESC及びiPSCの異種集団は、最も幅広く頑強な多系列発生能力を有するナイーブ状態幹細胞を含み、それ故に、再生医療における多数の人命を救う治療用途に、かなり有望である。ナイーブ様hESCの特徴としての高いLTR7/HERVH発現の定義と一致して、LTR7/HERVH発現が著しく上昇したhESC及びヒト誘導多能性幹細胞(hiPSC)の亜集団は、ナイーブ様多能性幹細胞の重要な特性を示す[19]。さらに、ヒトナイーブ様多能性幹細胞は、LTR7/HERVHの転写の上昇のマーカーに基づいて、遺伝的にタグ付けされ、首尾よく単離され、インビトロで維持され得る[19]。胚性幹細胞特異的転写因子NANOG、POU5F1、KLF4、及びLBP9は、ヒト多能性幹細胞においてLTR7/HERVH転写を駆動する[19]。HERVH活性及びHERVH由来転写物での標的化された干渉は、ヒト多能性幹細胞の自己再生機能を大幅に損なう[19]。
LTR7/HERVHサブファミリーと同様に、HERVKレトロウイルスの最も進化的に新しいゲノム組込み部位を表す、低メチル化LTRでの、多能性マスター転写因子POU5F1(OCT4)によるLTR5_Hs/HERVKのトランス活性化は、正常ヒト胚形成中にHERVK発現を誘導する[20]。それは、胚性ゲノム活性化が、8細胞期で起こり、着床前胚盤胞におけるエピブラスト細胞の段階にわたって続き、胚盤胞生成物からのhESC誘導の間に停止することと一致する[20]。ヒト胚形成中のHERVK活性化の明白な実験的証拠が、Growらによって報告されている[20]。彼らは、ヒト胚盤胞におけるHERVKウイルス様粒子及びGagタンパク質の存在を証明した。これは、内在性ヒトレトロウイルスが、初期のヒト胚発生中に活性があり機能的であるという考えを支持するものである。この仮説と一致して、多能性細胞におけるHERVKウイルス−アクセアリータンパク質Recの過剰発現は、宿主タンパク質IFITM1レベルを増加させてウイルス感染を阻害するのに十分であった[20]。これは、ヒトの初期段階の胚におけるこの抗ウイルス防御メカニズムが、HERVK活性化によって誘発され得ることを示唆している。レトロトランスポゾンの活性化が胚性幹細胞における種特異的遺伝子発現をどのように調整するかについての詳細な分析が、最近のレビューで示されており[21]、これは、ヒト細胞における特定のレトロトランスポゾンの活性化と抑制の間で、進化の間に確立された精密な制御バランスを浮き彫りにしている。
最近の実験は、hESCにおいてSCARの重要な生物学的活性を媒介する重要なエフェクター分子を同定した。SCAR由来の長い非コードRNAは、hESCの多能性、機能的同一性、及び完全性を維持するための必須の制御分子として記載されている[12]。集合的に、これらの実験は、hESC及びiPSCを含むヒト多能性幹細胞の多能性維持及び機能的同一性のための、特定の内在性レトロウイルスの、持続的ではあるがしっかりと空間時間的に制御された活性の必須の役割を決定的に確立した。SCARの覚醒が、がん細胞における幹性ゲノムネットワークの活性化及び複数のタイプの悪性腫瘍を有すると診断された患者における臨床的に致死的ながんによる死亡の表現型の出現に関連している可能性があるという仮説が立てられている[6−9]。
要約すると、新たな合意見解は、ヒト着床前胚中の内在性幹細胞関連レトロウイルス(SCAR)、具体的にはLTR7/HERVH及びLTR5_Hs/HERVKサブファミリーの、空間時間的に制御された活性化が、ナイーブ状態ESCの多能性維持、機能的同一性及び完全性、並びに初期段階のヒト胚の抗ウイルス抵抗性のために必要である、ということである。SCARの発現は、分化したヒト細胞においてエピジェネティックにサイレンシングされ、SCAR活性の制御及び効率的なサイレンシングに失敗すれば、分化に欠陥のある表現型を生じる。がん細胞におけるSCAR遺伝子座のエピジェネティックなサイレンシングの逆転は、複数のタイプのヒト腫瘍におけるSCAR発現の活性化に関連しているようである(9及びその参考文献において概説されている)。
この寄稿では、ヒト着床前胚の単細胞RNA配列解析は、卵母細胞から接合子への移行中の特定のLTR7/HERVH遺伝子座の活性化を明らかにし、ヒト胚における異数性に関連するHERVHネットワークシグネチャーを同定する。相関パターンの分析は、前立腺上皮内新生物から限局性及び転移性の前立腺がんへのがん進行経路の存在を支持する、前立腺腫瘍の臨床サンプルの遺伝子発現プロファイルを有するインビボの成熟したヒト卵母細胞のHERVHネットワーク活性化のトランスクリプトームシグネチャーとリンクしている。臨床的に致死的な疾患を有するがん患者からの悪性腫瘍における、多様な範囲にわたるゲノム異常の発現は、がん細胞におけるSCARネットワークの活性化に関連している。29のがんタイプを表す、全てのTCGAコホートにわたる、12,093の臨床サンプルにおける、SCARネットワークのがんゲノムアトラス(TCGA)にガイドされた分析は、SCAR関連タンパク質コード遺伝子及び非コードRNA遺伝子座の体細胞非サイレント突然変異、遺伝子発現、遺伝子レベルのコピー数変化、タンパク質発現によって規定される、臨床的に致死的な治療抵抗性疾患の汎がんゲノムシグネチャーを明らかにした。
実験例の説明
単細胞トランスクリプトーム分析は、異数性をよく生成し発生的に非生存可能なヒト接合子におけるLTR7/HERVH制御遺伝子の発現の変化及び選択されたLTR7/HERVH遺伝子座からの活性転写を明らかにする。
染色体不安定性は、初期段階のヒト胚発生において一般的であり、異数体は卵割期ヒト胚の50−80%で観察される[Vanneste E, Voet T, Le Caignec C, Ampe M, Konings P, Melotte C, Debrock S, Amyere M, Vikkula M, Schuit F, Fryns JP, Verbeke G, D'Hooghe T, Moreau Y, Vermeesch JR. Chromosome instability is common in human cleavage-stage embryos. Nat Med. 2009; 15 :577-83; Johnson DS, Gemelos G, Baner J, Ryan A, Cinnioglu C, Banjevic M, Ross R, Alper M, Barrett B, Frederick J, Potter D, Behr B, Rabinowitz M. Preclinical validation of a microarray method for full molecular karyotyping of blastomeres in a 24-h protocol. Hum Reprod. 2010; 25: 1066-75; Chavez SL, Loewke KE, Han J, Moussavi F, Colls P, Munne S, Behr B, Reijo Pera RA. Dynamic blastomere behaviour reflects human embryo ploidy by the four-cell stage. Nat Commun. 2012; 3:1251; Vera-Rodriguez M, Chavez SL, Rubio C, Reijo Pera RA, Simon C. Prediction model for aneuploidy in early human embryo development revealed by single-cell analysis. Nat Commun. 2015; 6: 7601; Yanez LZ, Han J, Behr BB, Pera RA, Camarillo DB. Human oocyte developmental potential is predicted by mechanical properties within hours after fertilization. Nat Commun. 2016; 7: 10809]。
ヒト胚における異数体は、細胞周期停止、細胞生存能力の喪失、及び発生不全につながる、適切な発生を損なう。単細胞トランスクリプトーム分析は、接合子の遺伝子発現シグネチャーが、正倍数体及び異数体ヒト胚の発生を確実に予測し、発生的に生存可能な接合子と非生存可能な接合子とを区別することができたことを証明した[Vera-Rodriguez M, Chavez SL, Rubio C, Reijo Pera RA, Simon C. Prediction model for aneuploidy in early human embryo development revealed by single-cell analysis. Nat Commun. 2015; 6: 7601; Yanez LZ, Han J, Behr BB, Pera RA, Camarillo DB. Human oocyte developmental potential is predicted by mechanical properties within hours after fertilization. Nat Commun. 2016; 7: 10809]。
特定のLTR7/HERVH遺伝子座の活性化がヒト胚における異数体の発生に関連しているという仮説の妥当性試験は、これらの実験的パラダイムに適合し、以下の仮定に従わなければならない:
・高いLTR7/HERVH発現は、ヒト接合子において容易に検出可能であるはずである;
・接合子段階でLTR7/HERVH遺伝子座が活性化された細胞は、ヒト胚形成の後続の段階の間、持続しないはずである;及び
・異数性をよく生成するヒト胚の遺伝子発現シグネチャーは、顕著な数のLTR7/HERVH制御遺伝子を保持するはずである。
ヒト胚発生関連遺伝子の分析は、LTR7/HERVH制御遺伝子の数が、正倍数体胚と比較して、異数体胚において異なって発現される遺伝子の間で顕著に富化されていることを証明している(表1A)。対照的に、ヒト胚発生関連遺伝子の6の異なる遺伝子発現カテゴリーを表す他の遺伝子セットにおいて、LTR7/HERVH制御遺伝子の顕著な富化は実証されなかった(表1A)。LTR7/HERVH遺伝子座の活性化がヒト胚における異数体の発生に関連しているという仮説と一致して、shHERVH処置hESCの遺伝子発現シグネチャーと異数体胚対正倍数体胚で異なって発現される遺伝子で構成される接合子対8細胞胚の遺伝子発現プロファイルとの間で、顕著な相関が観察された(図1)。対照的に、shHERVH処置hESCの発現シグネチャーと異数体胚対正倍数体胚間で異なって発現されない遺伝子で構成される接合子対8細胞期胚の遺伝子発現プロファイルとの間では、顕著な相関は実証されなかった(図1)。HERVH制御遺伝子の発現が異なる発生能力を有するヒト接合子を区別するという考えと一致して、発生的に生存可能な接合子対非生存可能な接合子において異なって発現される全ての遺伝子の50パーセントが、hESCにおいてLBP9/HERVHによって制御される遺伝子で構成されたことが観察されている(図1)。
次に、LTR7/HERVH発現の活性化が、卵母細胞の受精後の胚形成の初期に起こり、それ故に、ヒト着床前胚の単細胞トランスクリプトーム分析中にヒト接合子において容易に観察され得る、という予測の妥当性を試験した。この考えと一致して、受精ヒト卵母細胞の接合子への移行中に、いくつかの規定されたLT7/HERVH遺伝子座の顕著な活性化が観察された(図2)。特に、接合子における高いLTR7/HERVH発現は、限定された数の特定のLTR7/HERVH遺伝子座のみに制限され、8細胞期を超えて持続しなかった(図2)。予想通り、ほとんどのLTR7/HERVH遺伝子座は、初期胚形成期中にはサイレントのままであり、後期胚盤胞期中に、エピブラスト形成期中に、及びhESC生成の開始時に、膨大な量の活性化を受ける[1−14;16−21]。仮説と一致して、接合子においてLTR7/HERVH遺伝子座の活性化された細胞の圧倒的多数は、ヒト胚形成の後続の段階の間持続しなかった(図2)が、パターン4細胞は、胚形成のエピブラスト期及びhESC生成段階でLTR7/HERVH発現の著しい増加を示した。ヒト胚形成中に発現プロファイルのパターン4を示すLTR7/HERVH遺伝子座の活性化は、基底状態多能性状態及びナイーブhESCの生成に関連している可能性が高い。この仮説は、hESCの多能性制御及び維持ネットワークにおいて重要な役割を果たすHPAT3 lincRNAのLTR7/HERVH配列の発現プロファイルの単細胞トランスクリプトーム分析によってさらに裏付けられる(図2)。
ヒト卵母細胞におけるLTR7/HERVHネットワーク活性化の遺伝子発現シグネチャーは、前立腺がん前駆病変、限局性及び転移性の前立腺がんを、正常前立腺上皮及び良性前立腺肥大と区別する。
胚形成の間、胚性ゲノム活性化の前に転写は起こらない。これは、胚形成の初期段階が、卵母細胞から排他的に受け継がれる母方の遺伝情報によって排他的に制御されることを示している。胚性ゲノムの転写活性化の主要な波が、ヒト胚形成の4細胞期から8細胞期で観察された[Dobson AT, Raja R, Abeyta MJ, Taylor T, Shen S, Haqq C, Pera RA. The unique transcriptome through day 3 of human preimplantation development. Hum. Mol. Genet. 2004; 13: 1461-1470]。これらの考察は、ヒト接合子において観察される高いHERVH遺伝子座の発現が、卵母細胞におけるそれらの活性転写状態に関連している可能性があることを示唆している。この考えと一致して、卵巣からの除去後数分以内に得られたヒト分裂中期II卵母細胞のトランスクリプトームの分析[Kocabas AM, Crosby J, Ross PJ, Otu HH, Beyhan Z, Can H, Tam WL, Rosa GJ, Halgren RG, Lim B, Fernandez E, Cibelli JB. The transcriptome of human oocytes. Proc Natl Acad Sci U S A. 2006; 103: 14027-32]は、異なって発現されたHERVH制御遺伝子の大きいセットを同定した(図1)。さらに、ヒト着床前胚の単細胞トランスクリプトーム分析は、ヒト卵母細胞における選択されたLTR7/HERVH遺伝子座の発現の直接の実験的証拠を明らかにした[図2]。ヒト卵母細胞におけるLTR7/HERVHネットワーク活性化の遺伝子発現シグネチャーの同定は、この遺伝子シグネチャーが悪性腫瘍の臨床サンプルにおけるLTR7/HERVHトランスクリプトーム活性化の検出に有用であり得るか否かを決定する機会を提供する。注目すべきことに、この分析は、ヒト卵母細胞におけるLTR7/HERVHネットワーク活性化の遺伝子発現シグネチャーが、前立腺がん前駆病変、限局性及び転移性の前立腺がんを、正常前立腺上皮、間質、及び良性前立腺肥大の臨床サンプルから区別するようであるということを明らかにする(図3)。
これらの観察は、LTR7/HERVHトランスクリプトームの活性化が、前立腺がん前駆病変(サンプルの31−46%)、限局性前立腺腺がん(サンプルの22−28%)、及び転移性前立腺がん(サンプルの45−60%)を構成する、前立腺上皮内新生物の臨床サンプルの大きいサブセットにおいて起こることを強く示唆している。集合的に、これらの結果は、LTR7/HERVH制御ネットワークの活性化が、臨床的に重要な前立腺がんの発症の初期に起こり、推定上の前駆病変(前立腺上皮内新生物)から限局性及び転移性の前立腺がんへの前立腺がん進行中に持続を示す、と主張する。
ヒト特異的キメラ宿主/ウイルス転写物の発現の差は、がん患者を、著しく異なる長期生存確率を有する複数のサブグループに分離する。
SCARの覚醒が、がん細胞における幹性ゲノムネットワークの活性化及び複数のタイプの悪性腫瘍を有すると診断された患者における臨床的に致死的ながんによる死亡の表現型の出現に関連しているという仮説が立てられている[6−9]。hESCゲノムの規定された領域へのSCARの挿入は、代替プロモーター、エキソン化、及び選択的スプライシングを生成することによって、宿主遺伝子及びキメラ宿主/ウイルス転写物の発現に、顕著に影響を及ぼすようである(18−20)。これらのデータは、SCARネットワークの活性化のゲノムシグネチャーが、SCAR由来の転写物、SCAR制御タンパク質コード遺伝子、キメラ宿主/ウイルス転写物、及び非コードRNAを含む、遺伝要素の様々な異なるクラスから成る可能性があることを示唆している。興味深いことに、完全長LTR7/HERVH遺伝子座の約75%が、ヒト及び非ヒト霊長類において高度に保存されているように見えるが(表1)、300超の遺伝子座が候補ヒト特異的調節エレメントを表しており、したがって、保存された霊長類特異的且つヒトに特有の制御SCAR由来の配列の生物学的役割の調査の必要性をはっきりと示している。注目すべきは、完全長ヒト特異的LTR7/HERVH配列が、不活性のLTR7/HERVH遺伝子座と比較して、転写活性のある遺伝子座の間で顕著に富化されていることである(表1)。それ故に、宿主/ウイルスキメラ転写物の構造成分を含むタンパク質コード遺伝子のmRNA発現プロファイルは、ヒト腫瘍における遺伝子座特異的SCAR活性化の潜在的な臨床的関連性の評価に有用であり得る。
SCAR活性化の潜在的な臨床的関連性を評価するために、カプラン・マイヤー生存率分析によって規定されるがん患者の長期生存確率に関連する、宿主/ウイルスキメラ転写物の構造成分を含むタンパク質コード遺伝子のmRNA発現レベルの変化のパターンを評価した(図1)。この分析の主眼は、ヒト特異的SCAR挿入を保持し、それ故に、候補ヒト特異的制御配列として規定された、宿主/ウイルスキメラ転写物であった(表1−3)。
12のTCGAコホート(12の異なるがんタイプのPANCAN12研究)にわたる5,158の臨床サンプル及び全てのTCGAコホートにわたる12,093の臨床サンプルを含む、2つのTCGA汎がんデータベースの問い合わせ(https://genomecancer.soe.ucsc.edu/proj/site/xena/datapages/)は、SCAR標的化タンパク質コード遺伝子の遺伝子発現及び遺伝子コピー数の変化が、がん患者の長期生存率と2つの異なる関連パターンを示すことを証明する(図1)。
関連パターンの1つは、SCAR標的化遺伝子の高い遺伝子発現レベルが、低いがん患者の生存の可能性に関連しているようである、という観察によって規定される。このパターンは、PLCXD1及びCCL26遺伝子について観察された(図1)。対照的に、第2の関連パターンは、SCAR標的化遺伝子の低い遺伝子発現レベルが低いがん患者の生存確率に関連しているという証拠によって示される。このパターンは、ZNF443、LRBA、TPT1、ABHD12B、及びLIN7A mRNAについて観察された(図1)。
TCGA乳がんコホート(1,241臨床サンプル);TCGA前立腺がんコホート(568臨床サンプル);及びTCGA直腸がんコホート(187臨床サンプル)を含む、がんタイプ特異的な患者の生存プロファイルの分析中に、TCGA汎がんデータセットと同様の関連パターンが観察された(図1B)。特に、前立腺及び直腸がんを有すると診断された患者の間で、診断及び治療後10年を超えて約100%の生存確率を有する個体で構成される患者の予後良好のサブグループを同定することが可能であるように思われる(図1及び補足図S2)。それ故に、宿主/ウイルスキメラ転写物の構造的要素を含むヒト特異的レトロウイルス挿入を有するSCAR標的化タンパク質コード遺伝子のmRNA発現レベル及び遺伝子コピー数の変化は、悪性腫瘍において観察される様々な異なるSCARの活性化パターンが、がん患者における臨床的に異なる成績に関連している、という仮説と一致しているようである。
高いがんによる死亡の可能性に関連する体細胞非サイレント突然変異のフィンガープリント。がん患者の効率的なエビデンスベースの個別化された管理、及び新規の診断、予後、及び治療の用途の開発のために、治療失敗、疾患再発、転移性進行、及び最終的にはがんによる死亡の高い確率によって規定される、悪性腫瘍の臨床的難治性の体細胞非サイレント突然変異の遺伝的シグネチャーを同定することは、特に有用であろう。この目的のために、SCARのゲノムネットワーク及びがんドライバー遺伝子を、体細胞非サイレント突然変異を獲得した遺伝子(腫瘍サンプルにおけるその検出は、高いがんによる死亡の可能性に関連する)について系統的に探索した。このタイプの関連の、複数の統計的に有意な事例が観察された:すなわち、SCAR関連ゲノムネットワークの遺伝子は、悪性腫瘍において体細胞非サイレント突然変異(SNM)を獲得し、これらの突然変異を有する腫瘍を有するがん患者は、顕著に低い長期生存確率及び高いがんによる死亡の可能性を示した(図5)。これらの観察は、臨床的に致死的ながんによる死亡の表現型の遺伝的ドライバーとして機能し得る遺伝子がSCAR関連ゲノムネットワーク内に存在することを示唆した。反対に、がんドライバーとして既に規定されている遺伝子のいくつかが、候補SCAR制御遺伝子のカテゴリーを構成し得ると期待することは妥当であった。
この仮説は、いくつの既報の候補がんドライバー遺伝子が、独立した実験においても、候補SCAR制御遺伝子(これはshRNAアプローチを使用して最近発見された[19])として同定されたかを決定することによって、試験した。高信頼性のがんドライバー遺伝子として報告された[22]291遺伝子のうち合計183遺伝子(63%)が、hESCの候補HERVH/LBP9制御遺伝子として同定された。同様に、ヒト腫瘍において顕著に変異した遺伝子として既に同定されている[23]127の遺伝子のうち75遺伝子(59%)が、候補HERVH/LBP9制御遺伝子であると報告された。最後に、がん遺伝子センサス(Cancer Gene Census)(http://cancer.sanger.ac.uk/census)の最新のリリースの572遺伝子のうち325遺伝子(57%)が、hESCの候補HERVH/LBP9制御遺伝子として同定された。集合的に、これらの観察は、腫瘍サンプルの複数のコホートにわたりポジティブセレクションのシグナルを示し、候補がんドライバー遺伝子として規定された遺伝子の大多数が、hESCにおいてHERVH/LBP9幹性経路によって制御されるようである、ということを示している。
これらの考察に基づいて、低い生存確率及び高いがんによる死亡の可能性を有する患者を分離する、18遺伝子のがんによる死亡のSNMのシグネチャーを同定した(図5)。腫瘍サンプルから単離されたこれらの18の遺伝子のそれぞれにおける体細胞非サイレント突然変異の検出は、腫瘍にこれらの遺伝子の体細胞非サイレント突然変異を有しない患者と比較して、がん患者の不良な長期予後に関連しているようである(図5)。意義深いことに、全てのがん死亡事象の約70%が、18遺伝子のがんによる死亡突然変異のシグネチャーによって規定される予後不良患者のサブグループにおいて起こった一方で、TP53突然変異シグネチャー単独が死亡事象の50%未満を占めたことが観察された(図5)。がんによる死亡のSNMのシグネチャーを含む18の遺伝子は、ヒト遺伝子であって、体細胞非サイレント突然変異の存在が、全てのTCGAコホートにわたる7,509の腫瘍サンプルの単一の汎がんデータセットにおいて検出され、1,934から8,272までの腫瘍サンプルに及ぶ9の汎がんデータセットのフォローアップ分析中に確認された、ヒト遺伝子を表す(ただし、腫瘍におけるこれらの突然変異の存在が、カプラン・マイヤー生存率分析によって規定される顕著に高いがんによる死亡の可能性に関連している、という要件が満たされている(下記参照))。特に、追加の9の顕著なSNM遺伝子がカプラン・マイヤー生存率分析に含まれていた場合に、SNMシグネチャーの分類力は、ほんのわずかに増加するようである(図5)。
SNM遺伝子のがん生存可能性分類性能を、いくつかの追加の分析を使用して確認した(補足図S3)。これらの分析において、フォローアップ生存データの完全な臨床記録を有する患者のみを含めた。腫瘍にSNMを有する及び有しない7,258人のがん患者のカプラン・マイヤー生存率分析の比較は、2つのSNM遺伝子を有する患者(生存期間中央値1,725日)又は1つのみのSNM遺伝子を有する患者(生存期間中央値1,944日)と比較して、腫瘍に少なくとも3つのSNM遺伝子を保持するがん患者が、最も短い生存期間中央値(1,438日)を示したことを証明する。腫瘍にSNM遺伝子を有しないがん患者は、最も長い生存期間中央値(4,068日)を有した。7,258人のがん患者を、生存期間の昇順にソートした後3つの同一のサイズ(n=2,419)のサブグループに階層化した場合、360日の生存期間中央値を有する患者の63.4%が、腫瘍にSNM遺伝子を有した一方で、それぞれ869日及び4,222日の生存期間中央値を有するがん患者の58.5%及び51.8%が、腫瘍にSNM遺伝子を有した(補足図S3A)。がんによる死亡の表現型のSNMシグネチャーを保持する遺伝子の突然変異のフィンガープリントの可視化は、臨床腫瘍サンプルから単離されたこれらの遺伝子には、突然変異が「散らばっている」ように見え、その圧倒的多数がSNMによって表されることを明らかにした(補足図S3B)。
興味深いことに、18のうち11(61%)のがんによる死亡のSNMのシグネチャー遺伝子が、15のヒト特異的NANOG結合部位の近くに位置している[3]。これは、これらの遺伝子が、hESCにおけるNANOG制御ネットワークの遺伝要素を表している可能性があることを示唆している。11のがんによる死亡のSNMのシグネチャー遺伝子の近くの、15のヒト特異的NANOG結合部位の配置は、偶然にのみ期待され得るよりも有意に高い(p=9.95E−05;超幾何分布検定)。これは、がんによる死亡のSNMのシグネチャー遺伝子の近くの顕著な配置富化を示さない、他のヒト特異的転写因子結合部位(CTCF;POU5F1;RNAPII)とは対照的である(データは示さず)。特に、これらの18の遺伝子全ての遺伝子コピー数の変化は、不良ながん患者の長期生存率に関連しているようであり(補足図S4)、したがって、遺伝子特異的な遺伝的変化の検出のための独立した分析エンドポイントを使用して、この遺伝子パネルの潜在的な診断値及び予後値を確認している。
次に、高いがんによる死亡の可能性に関連するSNMの遺伝子検出の探索を、複数の汎がんデータセット(下記参照)を採用して、ヒトがんにおいて顕著に変異した127の遺伝子[23]及びがんの体細胞突然変異のカタログ(catalogue of somatic mutations in cancer)、COSMIC(http://cancer.sanger.ac.uk/cosmic/census)に列挙された177の遺伝子を問い合わせて、実施した。全体として、42の遺伝子を同定した。これらの遺伝子は、悪性腫瘍の臨床サンプルにおいて体細胞非サイレント突然変異を獲得した。そして、これらの突然変異の存在は、不良な治療成績及びがんによる死亡の顕著に高い可能性に関連している(補足データセット3)。特に、がんによる死亡の表現型の突然変異のフィンガープリントを保持する、42の遺伝子のうち33の遺伝子(78.6%)が、SCAR関連ゲノムネットワークのメンバーを構成する(補足表S1及び補足データセットS3)。
高いがんによる死亡の可能性に関連するSNMのシグネチャーの妥当性確認解析。ゲノムワイドな高スループット実験における体細胞非サイレント突然変異(SNM)の検出は、重要な実験的且つ分析的課題を突きつける。SNMのコールは、同じDNAサンプルの処理中でさえ、多数の要因の影響を受ける。ライブラリー調製及びシーケンシングプラットフォームなどの技術的要因に加えて、シーケンシングリード及びコーリングアルゴリズムのマッピング、参照ゲノムデータベースの選択、ゲノムアノテーション、及び標的選択領域などの、分析及び計算の方法論の差異は全て、SNMの同定に寄与する。最終的に、遺伝子及びサンプルのブラックリストなどのアドホックプレ/ポストデータ処理の差異は、交絡要因となる可能性がある。これらの潜在的な変動原因を説明するために、がん患者の生存率とSNMコールとの間の関連の有意性を調べた。その際には、UCSC Xenaブラウザで利用可能な汎がんデータセットのために様々な異なる研究チームによって報告された体細胞非サイレント突然変異コールのデータベースを使用した。全体として、1,934から8,272までの腫瘍サンプルを含む10の汎がんデータセットを、この分析において評価した(補足データセットS3)。SNMのがんによる死亡の表現型シグネチャーの18の遺伝子全て(図5)を、少なくとも2つの汎がんデータセットにおいて統計的に有意な遺伝子としてスコア付けした(補足データセットS3)。18のSNMのシグネチャー遺伝子のうち17(94.4%)が、少なくとも3つのデータセットにおいて、統計的に有意な遺伝子として同定された。その中のSNMの突然変異は、カプラン・マイヤー分析によって規定される高いがんによる死亡の可能性に関連していた(補足データセットS3)。同様に、42の遺伝子のうち39の遺伝子(92.9%)におけるSNMの検出は、少なくとも2つの汎がんデータセットにおいて、顕著に高いがんによる死亡の可能性に関連していた(補足データセットS3)。まとめると、これらの観察は、ここで同定された遺伝子が、決定的な突然変異標的部位特異的妥当性確認実験及びフォローアップ構造−機能及び機構的研究の正当性を証明するのに十分に頑強な、有望な候補遺伝マーカーを表していることを、主張しているようである。
複数のタイプの悪性腫瘍を有すると診断された患者における悪性腫瘍の臨床的難治性の線形回帰分析は、がんによる死亡の可能性、がんタイプ、及び腫瘍におけるSNMのがんによる死亡シグネチャーの存在の間の関連の、顕著な証拠を明らかにした(図5)。1つの分析において、28のがんタイプのTCGA汎がんコホートからのがん患者の生存データを利用して、各がんタイプについての死亡事象のパーセントを計算した。得られた値を、がん患者の対応する群におけるSNMのがんによる死亡シグネチャーを有する患者のパーセントと整列させて、線形回帰分析に供した(図5C)。別の分析において、19がんタイプについての米国における年齢調整されたがん発生率及び死亡率(100,000人あたり)を、疾病管理予防センター(CDC)米国がん統計(USCS)報告から得た。各がんタイプについての推定死亡率は、対応する発生率の値及びSNMのがんによる死亡シグネチャーを有する患者のパーセントの値を掛け合わせることによって計算した。推定死亡率の値を、対応するがんタイプについての実際の死亡率と整列させて、回帰分析に供した(図5D)。いずれの場合においても、ひときわ顕著な相関が観察された。これは、腫瘍におけるSNMのシグネチャーの存在が、臨床的に致死的な疾患を発症する可能性が高いという分子シグナルを表し得るという仮説を、強く支持している。
集合的に、本分析は、臨床腫瘍サンプルにおけるSCAR関連ゲノムネットワークの規定された遺伝要素の活性化の分子的証拠が、悪性腫瘍の診断及び治療後の不良ながん患者の長期生存率によって規定される、臨床的に致死的ながんによる死亡の表現型の発現の高い可能性とリンクしているようであることを示す。不良ながん患者の生存率と、SCAR活性及びhESCにおける幹性ネットワークの維持のマスター転写制御因子を構成する遺伝子、すなわちKLF4、LBP9、POU5F1、及びNANOG、のコピー数変化との間に観察された顕著な相関は、この仮説を強く支持する(補足図S4)。これらのデータは、がん細胞におけるSCAR関連ゲノムネットワークの活性化が、選択的成長及び/又は生存上の利点を提供し、悪性進行中にポジティブセレクションの遺伝的シグナルを表し得ることを示唆している。
この結論は、TCGA PANCAN12コホートの臨床サンプルにおける、SCAR制御遺伝子によってコードされるタンパク質の発現の分析によって、さらに支持される(図6)。カプラン・マイヤー生存曲線に関連する全ての利用可能なタンパク質発現データを、38のHERVH/LBP9制御遺伝子について評価した。特に、23のSCAR制御遺伝子(60.5%)のタンパク質発現レベルの変化が、がん患者の長期生存確率と顕著な関連を示した(補足データセット1)。これらの高度に顕著な関連の例を図6に示しており、これは、SCAR関連幹性ゲノムネットワークの機能的変化が、がん患者における臨床的に致死的な疾患進行において役割を果たす可能性がある、という仮説を確認するものである。
本分析の結果に基づいて、ヒトがんの29の異なるタイプを表す全てのTCGAコホートにわたる12,093の臨床サンプルにおける、SCARのネットワークのTCGAにガイドされた調査が、体細胞非サイレント突然変異(SNM)の存在、遺伝子レベルのコピー数変化、SCAR制御宿主遺伝子の転写物の及びタンパク質の発現によって規定される、臨床的に致死的な治療抵抗性疾患の汎がんゲノムシグネチャーを明らかにした、という結論に達した。この通信文で報告される遺伝子は、決定的な突然変異標的部位特異的妥当性確認実験及びフォローアップ構造−機能及び機構的研究の正当性を証明するのに十分に頑強な、臨床的に致死的な形のヒトがんの有望な候補遺伝マーカーを表す。
複数の異なるSCARの遺伝子座の規定された遺伝的シグネチャーのゲノムワイドなマッピングは、ヒト特異的キメラ転写物によってコードされる保存タンパク質ドメインの、ヒトゲノムにおける著しい拡大を明らかにした。
ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内の保存タンパク質ドメインの分析は、複数の異なるSCARの遺伝子座が、保存タンパク質ドメインのコンビナトリアルパターンによって規定される異なるタンパク質コードシグネチャーを示すことを証明する(図2及び補足図S1)。個々のSCARの配列の系統的BLAST分析は、ウイルス配列の突然変異が、機能性ウイルスタンパク質の完全なコード能力を低下させたこと、及びある特定の保存タンパク質ドメインの残存構造のみが保存された状態を維持することを証明する(図2及び補足図S1)。特に、ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内で最も頻繁に表される保存タンパク質ドメインの1つは、GVQWアミノ酸配列である(図2及び3)。GVQWアミノ酸配列をコードする複数の異なるSCARの遺伝子座のヌクレオチド配列が容易に区別可能であるため、複数の異なるSCAR遺伝子座が蒔かれたヒトゲノム中のGVQWをコードする配列の数を特定することが可能であった。この分析が、ヒトゲノム全体にわたるGVQWドメインの拡散に対する、複数の異なるSCAR遺伝子座の拡大の相対的影響の評価に有用であり得る、という仮説が立てられている。
保存されたGVQWタンパク質ドメインをコードする複数の異なるSCARの遺伝子座の特定の遺伝的シグネチャーのゲノムワイドなマッピングは、ヒトゲノム全体にわたり散乱した数千の遺伝子座特異的遺伝的フィンガープリントを同定した。これは、親のSCARの配列と比較して、ギャップ又は挿入なしで100%配列同一性を有するヌクレオチド配列として規定された(図3)。注目すべきことに、この分析は、ヒトゲノム中のGVQW保存タンパク質ドメイン配列をコードするDNA配列の大多数が、ChrY:278899 -284215及びchrX:278899 -284215上のDNA配列に由来するヒト特異的キメラ転写物に由来しているようであることを明らかにした(図3)。特定のSCAR由来のヌクレオチド配列のこの拡大は、他の大型類人猿と比較したときの、ヒトプロテオーム内のGVQW保存タンパク質ドメインの著しい富化に寄与した可能性がある(図3)。
さらなる分析は、ジンクフィンガータンパク質が、GVQWドメインを保持するヒトゲノム中の最大のタンパク質ファミリーの1つを表すことを明らかにした。それ故に、GVQWドメインを保持するジンクフィンガータンパク質の発現が、治療後に異なる長期生存率を有するがん患者の悪性腫瘍において変化しているかどうかを決定することが、関心の対象となった。注目すべきことに、この分析は、GVQWドメインを保持するジンクフィンガータンパク質のmRNA発現レベル及び遺伝子コピー数の変化が、著しく異なる処置成績のサブグループにがん患者を分離するようであることを証明する(補足図S2)。遺伝子発現及び遺伝子コピー数の変化の観察されたパターンは、前立腺、乳、結腸、直腸、及び膵臓がんを有すると診断された患者の間での、治療失敗及びがんによる死亡の高い可能性を有する個体の同定に有用であるようである(補足図S2)。GVQWドメインの機能が何であるか、及びこのドメインの特定のタンパク質配列への挿入が、宿主タンパク質の構造−機能特性にどのようにして影響を及ぼすのかを、実験的に決定することが、関心の対象となるだろう。
注目すべきことに、GVQW保存タンパク質ドメインを有する21の全てのジンクフィンガータンパク質及び3のSCARネットワークジンクフィンガータンパク質遺伝子(ZNF443;ZNF587;ZNF814)の遺伝子レベルのコピー数変化は、TCGA汎がんコホートを含む12,093の臨床サンプルの、予後不良及びカプラン・マイヤー生存率分析によって規定される高いがんによる死亡の可能性と、高度に顕著な関連を示す(図4)。これらのデータは、がん患者の臨床管理のための、及び治療失敗及び疾患進行のリスクが高い個体の同定のための、保存されたGVQWドメインを含むジンクフィンガータンパク質の潜在的な診断的及び予後的価値に関する結論を強化する。
内在性ヒトSCARによってコードされるヒト特異的制御配列の生成における、DNA修復経路の推定上の役割。
哺乳動物細胞は、損傷の近傍で利用可能な、ほぼ全てのDNA分子で、DNA二本鎖切断(DSB)にパッチを当てることができる、高度に有効なDNA修復経路を効率的に利用するように進化してきた[24、25]。転位因子(TE)由来DNA配列(DNAトランスポゾン並びにLTR及び非LTR両方のレトロトランスポゾンを含む)の、DNA損傷部位での挿入は、DSBを修復するために真核細胞によって利用されているようである[26−31]。TE由来DNA捕捉の代替モデル、DNA DSB修復部位でのエンドヌクレアーゼ非依存性L1挿入メカニズムが提案されている[27、28、30]。この経路は、DNA修復欠損げっ歯類細胞株において最初に観察された[27]。後続の報告は、このメカニズムがヒトゲノムにおいても機能する可能性が高いことを示した[28、30−32]。TE挿入の非古典的メカニズムが、Aluエレメント[31]及びHERV−Kレトロウイルス[32]によって媒介されるDSB修復に関連している可能性があることが示唆されている。SCAR活性が、ヒト細胞におけるDNA修復に寄与しているかどうかを特定することが、関心の対象となった。
レトロトランスポゾンの様々なクラスについて観察される非古典的TE挿入メカニズムのコンセンサスシグネチャー特徴は、TE由来配列の挿入部位内での祖先DNA配列の欠失である。TE媒介性DSBに関連するヒト特異的欠失はしばしば、祖先DNA配列の数千の塩基対にまで拡張される[31、32]。SCARがDSB修復経路に寄与しているかどうかを特定するために、内在性ヒトSCARによってコードされる候補ヒト特異的制御配列(HSRS)を同定し、制御DNAのヒト特異的獲得(挿入)及び喪失(欠失)の存在について分析した(表1、2)。予想通り、ヒト多能性幹細胞HSRSにおいて転写活性のあるものの大多数(75.0%−79.5%)は、ヒト特異的挿入を含む(表2)。注目すべきことに、ヒトの2013年12月のUCSCゲノムブラウザ(GRCh38/hg38)アセンブリ(UCSC Genome Browser on Human Dec. 2013 (GRCh38/hg38) Assembly)(http://genome.ucsc.edu/cgi-bin/hgTracks?db=hg38&position=chr1%3A90820922-90821071&hgsid=441235989_eelAivpkubSY2AxzLhSXKL5ut7TN)の20の哺乳動物(17の霊長類)のLiftOverアルゴリズム及びMultizアラインメントを使用するDNA配列保存分析は、SCARにコードされたHSRSの74.4%−88.6%が、チンパンジーゲノム及びボノボゲノムとの比較によって規定される祖先DNA配列の欠失を含むことを明らかにした(表2)。特に、SCARにコードされたHSRSの40.0%−59.1%は、長さが1,000bpを超えるDNAセグメントの大きい連続的ヒト特異的喪失を含む。最も極端な例のいくつかは、チンパンジーのゲノムと比較して27,843bpのヒト特異的欠失(hg38座標:chr4:132,117,632-132,124,853)及びボノボのゲノムと比較して81,108bpのヒト特異的欠失(hg38座標:chr4:3,927,445-3,933,080)を含む。同様に、75,171bp(chr12:8,279,022-8,294,090)、35,326bp(chr4:3,927,445-3,933,080)、及び71,036bp(chr1:112,809,666-112,826,054)の大きいヒト特異的欠失が、それぞれゴリラゲノム、オランウータンゲノム及びテナガザルゲノムと比較して、SCARの挿入の異なる遺伝子座で検出された。
本分析は、霊長類の進化の間に異なるヒト特異的DNA喪失の複数の独立した事象を経験した、ヒト多能性幹細胞において転写活性のある、SCARにコードされたヒト特異的制御遺伝子座を101同定した(表2)。ヒト特異的欠失のカスケードパターンを示すこれらの101の遺伝子座のゲノム座標を、20の哺乳動物(17の霊長類)のMultizアラインメントのUCSCゲノムブラウザトラックを使用して、ヒトDNA配列を、非ヒト霊長類のオルソロガスな配列と比較することによって同定した。この分析において、ヒトゲノム中の連続的ヒト特異的DNA配列が、チンパンジー、ボノボ、ゴリラ、オランウータン、及びテナガザルから構成される群から選択される、非ヒト霊長類の少なくとも2つの異なる種のゲノムと比較して、ヒト特異的欠失の少なくとも2つの異なる事象を示す場合に、HSRSを、ヒト特異的欠失のカスケードパターンを有するゲノム遺伝子座として規定した。それ故に、ヒト特異的欠失のカスケードパターンを示すゲノム遺伝子座は、霊長類の進化の間の長い期間にわたって、複数の異なる連続的DNAセグメントの反復喪失を経験するようであり、これは、これらのゲノム位置でのSCAR配列の挿入によって媒介される、DSBの発生及びDNA分子の修復の反復サイクルのメカニズムと一致するだろう。
ヒト特異的SCARの組込み部位のこれらの独特の構造的特徴は、SCAR関連DSB修復の分子メカニズムが、代替非相同末端結合(Alt NHEJ)として知られるバックアップDNA修復経路に類似している可能性があることを示唆している。なぜならば、Alt NHEJ経路によって構築される修復接合部の顕著な特徴は、大きいDNA欠失、挿入、及びマイクロホモロジーのトラクトであるからである[33、34]。集合的に、これらのデータは、DSB修復のAlt NHEJ経路が、特定のゲノム位置でのSCARの挿入に寄与し、これがヒト多能性幹細胞において転写活性のあるHSRSの生成をもたらした可能性がある、という仮説を支持する(図7)。
潜在的な、生物学的、病態生理学的、診断的、及び治療的意義の説明
リキッドバイオプシー用途のための意義
がん細胞のアポトーシス死及び/又は壊死の結果として、悪性腫瘍が、DNAの無細胞断片を血流中に放出するという観察は、がん細胞に由来する循環無細胞(cfDNA)の分析に基づくリキッドバイオプシーの概念の、実験的及び臨床的がん研究への開示及び迅速な導入の道を開く。がん細胞に由来するcfDNAの負荷は、腫瘍進行度分類及び予後と相関しているようである、という合意見解が出てきた[Diaz LA Jr, Bardelli A. Liquid Biopsies: Genotyping Circulating Tumor DNA. J Clin Oncol. 2014;32: 579-86; Haber, D. A. & Velculescu, V. E. Blood-Based Analyses of Cancer: Circulating Tumor Cells and Circulating Tumor DNA. Cancer Discov. 2014; 4: 650-661; Bettegowda, C. et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci. Transl. Med. 2014; 6: 224ra24; Newman AM, Bratman SV, To J, Wynne JF, Eclov NC, Modlin LA, Liu CL, Neal JW, Wakelee HA, Merritt RE, Shrager JB, Loo BW Jr, Alizadeh AA, Diehn M. An ultrasensitive method for quantitating circulating tumor DNA with broad patient coverage. Nat. Med. Nat Med. 2014; 20: 548-54; Dawson SJ, Tsui DW, Murtaza M, Biggs H, Rueda OM, Chin SF, Dunning MJ, Gale D, Forshew T, Mahler-Araujo B, Rajan S, Humphray S, Becq J, Halsall D, Wallis M, Bentley D, Caldas C, Rosenfeld N. Analysis of circulating tumor DNA to monitor metastatic breast cancer. N. Engl. J. Med. 2013; 368: 1199-209; Garcia-Murillas I, Schiavon G, Weigelt B, Ng C, Hrebien S, Cutts RJ, Cheang M, Osin P, Nerurkar A, Kozarewa I, Garrido JA, Dowsett M, Reis-Filho JS, Smith IE, Turner NC. Mutation tracking in circulating tumor DNA predicts relapse in early breast cancer. Sci Transl Med. 2015;7: 302ra133]。次世代シーケンシング技術の最新の進歩は、腫瘍由来DNAの分析の感度、特異性、及び正確さを顕著に向上させた。原理的には、最先端の次世代シーケンシング技術の状態は、従来はがん細胞の直接分析によってのみ実証可能であった、体細胞ゲノム改変のための腫瘍由来cfDNAの遺伝子型判定を可能にしてきた。cfDNAベースのアッセイを使用して、個々の腫瘍における非常に不均一な一連の突然変異を容易に検出し確実に定量化する能力は、がん患者における個人に合わせた疾患管理の概念の臨床的実行を容易にする、様々な翻訳アプリケーションに使用することができる、腫瘍進化のダイナミクスのリアルタイムでの追跡において、非常に効率的であることが証明されている。
複数の翻訳アプリケーションに大きな期待が持たれているが、現在の形のリキッドバイオプシー技術には、大きな制限がある。これらの制限は、初期段階の固形腫瘍の診断又はがんドライバー遺伝子の治療的に実施可能な突然変異の予測同定のための、リキッドバイオプシーの意図される使用が、注意深く検討されるときに、特に明らかである。その現在の形において、リキッドバイオプシーは主として、予め選択されたがんドライバー遺伝子のセットにおいて体細胞突然変異を確実に検出しようという主な意図で、血液サンプル(血漿又は血清)から抽出されたcfDNAの徹底的高分解能シーケンシングに利用される。がん細胞由来cfDNAが血液中に出現するためには、腫瘍の血管新生が必要とされるだろうと予測することは、妥当だと思われる。しかしながら、がん患者における本質的に全ての固形腫瘍の発生の初期段階は、血管新生の必要がないことによって特徴付けられ、そして実際、長年にわたり、拡散によって十分な栄養物の供給を受ける、腫瘍発生及び進行の無血管段階を表すことは、十分に確立されている。この文脈において、血液中の腫瘍由来cfDNAの出現は、腫瘍の血管新生の証拠及び転移性疾患への悪性進行の高い可能性の分子シグナルとみなされるべきである。この論理の道筋と一致して、腫瘍由来cfDNAが、進行した固形腫瘍を有するがん患者の>90%の血液中に、確実に且つ再現性よく検出される一方で、初期段階のがんを有すると診断された患者の血液中での検出率は、約50%(又はそれ以下)に低下する。重要なことに、次世代シーケンシング技術の分析性能のさらなる向上がこれらの現実を劇的には変えないだろうということは、ほぼ確実である。
がん細胞由来cfDNAの主な起源が、アポトーシス死及び/又は壊死を受けている腫瘍細胞に由来するということが合意見解であるようである。血液サンプルから抽出された腫瘍由来cfDNAの起源が、生存可能で活発に分裂しているがん細胞又は腫瘍成長を維持するがん細胞の少数の亜集団、例えば元のがんの細胞、腫瘍開始細胞、又はがん幹細胞など、に由来することを一貫して証明する信頼できる証拠はない。それ故に、がん患者の血液から抽出された腫瘍由来cfDNAの突然変異シグネチャーが、腫瘍進化の過去の歴史を表しており、死んだがん細胞から抽出された遺伝情報に基づいて、リアルタイム突然変異状態を識別するための、又は腫瘍進化の将来を予測するための、信頼できる方法はない、と考えることは妥当である。
単一核分解能での、腫瘍進化中のゲノムワイドな突然変異ダイナミクスの最新の分析は、体細胞点突然変異は、異数体とは対照的に、徐々に進化し広範なクローン多様性を生じたことを明らかにした[Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi X, Chen K, Scheet P, Vattathil S, Liang H, Multani A, Zhang H, Zhao R, Michor F, Meric-Bernstam F, Navin NE. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014; 512: 155-160]。標的化された単一分子シーケンシングは、決定的に、腫瘍において検出される多様な点突然変異の多くが、腫瘍細胞集団の<10%の頻度で起こることを証明した。それとはかなり対照的に、異数体再配列は、腫瘍進化の初期に現れ、クローン増殖中に非常に安定したままであった[Wang, Y., et al. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014; 512: 155-160]。この寄稿は、異数体の発生とSCARネットワークの異常な活性とを結びつけ、活性化されたSCARの経路の遺伝子発現シグネチャーを、がん前駆病変、限局性腫瘍、及び転移性がんの臨床サンプルにおいて検出することができることを証明する。集合的に、これらの観察は、SCARネットワークの活性化及び関連するゲノム異常ががん前駆細胞において起こりやすく、転移性疾患への腫瘍進化及び進行を通してずっと持続することを、強く主張する。それ故に、ここで同定されたSCAR配列、SCAR/宿主遺伝子ハイブリッド配列、SCAR制御タンパク質コード遺伝子及び非コードRNA配列の検出は、リキッドバイオプシー技術を利用する、診断、予後、治療選択、及び疾患管理用途のための、目覚ましい機会を切り開くだろう。
核酸及びタンパク質を含む無細胞巨大分子は、エキソソームと呼ばれるナノスケールサイズの粒子にしばしば存在する。エキソソーム中のDNA分子及びRNA分子のパッケージングは、細胞外ヌクレアーゼによる分解からそれらを保護するようであり、microRNA及びlincRNAなどの生物学的に活性な核酸分子は、安定した状態を維持するようである。それ故に、リキッドバイオプシー分析のためのサンプル調製プロトコルは、エキソソーム富化及び精製ステップを含むことにより恩恵を受ける可能性が高い。
DNA修復におけるSCARの配列の推定上の役割及び転移性がん細胞の生存の増加
本分析は、DNA修復におけるSCARの妥当な生物学的役割を示唆する。その役割とは、不死がん細胞にとって特に有益である、即座の生存及び適応の利点を宿主細胞に提供することによって、レトロトランスポゾン駆動突然変異の潜在的に有害な効果を無効にすることができる、というものである。DNA修復経路の比較的高い活性にもかかわらず、hESCは、放射線誘導DNA損傷及びアポトーシスに対する高い感受性を示す[35、36]。hESCのアポトーシスに対する高い感受性は、DNA損傷に応える低いアポトーシス閾値に起因することが示唆されている[36]。それとはかなり対照的に、治療抵抗性悪性腫瘍、高度に転移性のがん細胞、及び循環腫瘍細胞における幹性経路の活性化の既報の実験的及び臨床的証拠は、様々な生体関連微小環境の変化及び様々な異なる化学的摂動によって誘導されるアポトーシスに対する、著しく高い抵抗性の発現との、遺伝的及び表現型の関連を、一貫して証明した[37−51]。正常ヒト多能性幹細胞と、活性化された幹性遺伝子ネットワークを有する腫瘍細胞の高度に悪性の集団との間の、ゲノム構造の根本的な差異によって規定される、これらの重要な生物学的差異は、がん幹細胞の腫瘍開始能力、容赦ない成長、自己再生、及び生存に関与している可能性が高い。腫瘍細胞におけるSCARの継続的な転写活性は、機能的ウイルスゲノムをコードすることが明らかに構造的にできないにもかかわらず、不変の潜在的に致命的な脅威である可能性がある。ヒトゲノムに組み込まれたSCARの配列の何千もの変異体が存在する。これは、SCARの遺伝子の多くの突然変異が、SCARの配列の内在性コピーでの組換えによって修復できることを示唆している。この仮説と一致して、必須のウイルス遺伝子に致死的欠失を有する変異型レトロウイルスの導入は、内在性DNA配列との相同組換えによって突然変異を修復した復帰変異ウイルスの拡散をもたらす可能性があることが証明されている[52]。
幹細胞関連レトロウイルスのゲノムネットワークは、臨床的に難治性の悪性腫瘍のシグネチャーを保持する
SCAR及び関連する幹性ゲノムネットワークの本分析は、ヒト特異的制御ネットワーク及び経路の発生に寄与した可能性のあるヒト特異的挿入及び/又は欠失を保持する遺伝子座に焦点を当てた。この戦略の選択のための主な論理の道筋の1つは、広く実証されている、ヒトと非ヒト霊長類との間での、がん発生率の明らかな主要な差異に基づく。非ヒト霊長類において、前立腺がんは本質的に存在せず、肺がんは非常に稀である(53−58)。全体的に、乳がん、前立腺がん、肺がん、結腸がん、卵巣がん、膵臓がん、及び胃がんを含む、一般的ながんの発生率は、約2%〜4%の範囲にあると推定される(53−57)。幹性ゲノムネットワークの回路内で作動する、ヒト特異的制御遺伝子座及び経路のヒトに特有の表現型の効果は、がん発生率におけるこれらの劇的な種特異的差異に寄与している可能性がある。
この考えに基づいて、最初の分析は、ヒト特異的SCAR挿入を保持する宿主/ウイルスキメラ転写物に焦点を当てた(表1−3;図1)。宿主/ウイルスキメラ転写物の構造的要素を含むヒト特異的レトロウイルス挿入を有するSCAR標的化タンパク質コード遺伝子のmRNA発現レベル及び遺伝子コピー数の観察された変化は、様々な異なるSCARの活性化パターンが、著しく異なるがん患者の長期生存率に関連しているという仮説を支持する。
次に、ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内の保存タンパク質ドメインの分析を行った。この分析は、複数の異なるSCARの遺伝子座が、保存タンパク質ドメインのコンビナトリアルパターンによって規定される異なるタンパク質コードシグネチャーを示すことを証明する(図2及び補足図S1)。ヒト特異的SCAR由来の宿主/ウイルスキメラ転写物によってコードされる翻訳されたアミノ酸配列内で最も頻繁に表される保存タンパク質ドメインの1つが、GVQWアミノ酸配列であることが観察された(図2及び3)。GVQWドメインをコードするヌクレオチド配列の規定されたSCAR遺伝子座特異的シグネチャーを使用して、ヒトゲノム中のGVQWアミノ酸配列をコードするDNA配列の大多数の起源が、ChrY:278899 -284215及びchrX:278899 -284215上のDNA配列によってコードされるヒト特異的キメラ転写物に由来することが、決定された(図3)。SCAR由来のヌクレオチド配列の拡散は、特定のGVQWをコードするDNA配列の著しい拡大、及び他の大型類人猿と比較したときの、ヒトプロテオーム内のGVQW保存タンパク質ドメインの約10倍の富化をもたらすようである(図3)。これらのデータは、継続的なSCAR活性の生物学的に重要な結果の1つが、ヒトゲノム全体にわたり、特定の保存タンパク質ドメインをコードするヌクレオチド配列を蒔いたことであると、強く主張している。
注目すべきことに、後続の分析は、GVQWドメインを保持するジンクフィンガータンパク質のmRNA発現レベル及び遺伝子コピー数の変化が、著しく異なる処置成績のサブグループにがん患者を分離することを証明する(図4及び補足図2)。遺伝子発現及びコピー数の変化の観察されたパターンは、前立腺、乳、結腸、直腸、及び膵臓がんを有すると診断された患者の間で、治療失敗及びがんによる死亡の高い可能性を有する個体を分離するようである(補足図2)。前立腺及び直腸がんを有すると診断された患者の間で、診断及び治療後10年を超えて約100%の生存確率を有する個体で構成される患者の予後良好のサブグループを同定することが可能であるように思われる(補足図2)。これは、個別化された、エビデンスベースの疾患管理意思決定プロセスにとって、非常に重要な臨床的意義を有する可能性がある。
SCAR活性の遺伝的シグネチャーが、診断用途及び予後用途に潜在的に有用であるかどうかを決定するために、SCARのゲノムネットワークを、体細胞非サイレント突然変異を獲得した遺伝子(腫瘍サンプルにおけるその検出は、高いがんによる死亡の可能性に関連する)について系統的に探索した。全てのTCGAコホートにわたる臨床腫瘍サンプルにおける獲得された体細胞非サイレント突然変異及び悪性腫瘍におけるこれらの突然変異の存在が、顕著に高いがんによる死亡の可能性に関連しているようであるという、この寄稿において、合計42のヒト遺伝子が同定された(図5;補足表S1;補足データセットS3)。がんによる死亡の表現型の突然変異のフィンガープリントを保持する遺伝子のかなりの大多数(42のうち33;78.6%)が、SCAR関連ゲノムネットワークのメンバーを構成する(補足表S1及び補足データセットS3)。したがってこれは、臨床腫瘍サンプルにおけるSCAR関連幹性ゲノムネットワークの規定される遺伝要素の活性化の分子的証拠が、カプラン・マイヤー生存率分析によって規定される臨床的に致死的ながんによる死亡の表現型の発現の高い可能性とリンクしているようである、ということを確認するものである。意義深いことに、全てのがん死亡事象の70%超が、がんによる死亡のSNMのシグネチャーによって規定される予後不良患者のサブグループにおいて起こったことが観察された(図5)。
この寄稿において報告した重要な結論の1つは、臨床腫瘍サンプルにおけるSCAR関連幹性ゲノムネットワークの規定される遺伝要素の活性変化の分子的証拠の検出が、悪性腫瘍の診断及び治療後の不良ながん患者の長期生存率によって規定される、疾患進行の臨床発現の高い可能性に関連しているようである、という観察に基づく。腫瘍におけるSCARネットワーク内の特定の遺伝子の関与の観察は、体細胞非サイレント突然変異及び遺伝子コピー数の変化の検出に基づいており、これは、がん細胞におけるSCAR関連ゲノムネットワークの活性変化が、選択的成長及び/又は生存上の利点を提供し、悪性進行中にポジティブセレクションの遺伝的シグナルを表し得ることを示唆している。意義深いことに、28のがんタイプを有すると診断された患者の長期生存率に基づいて特定された、悪性疾患の臨床的難治性は、腫瘍に体細胞非サイレント突然変異のシグネチャーを保持するがん患者のパーセンテージと直接相関している。それ故に、ここで報告されたがんによる死亡の表現型の遺伝的相関は、難治性のヒト悪性腫瘍に対する新規の診断、予後、及び治療用途の開発にとって、非常に魅力的な標的を表す可能性がある。
SCARのゲノムネットワークのヒト特異的構造−機能的特徴が、H. sapiensの生理及び病理の両方において独特の役割を果たし得るという考えと一致して、HERV−Hトランスクリプトームが、定方向選択の影響下でヒトにおいて最近進化してきたものであり、チンパンジーとヒトとの分離以来、宿主に対して検出可能な適応効果を発揮する可能性が高いことが報告されている(59)。病理学的及び生理学的状態におけるSCARの生物学的に重要な機能の調査は、感染性ウイルス粒子の検出及び単離に排他的に焦点を当てるべきではない。多くの他のHERVファミリーと同様に、SCARの配列の大多数は、進化中に複数の突然変異及び欠失を蓄積した。そして、HERV配列が自己複製性及び感染性であることは示されていない。
ヒトゲノムにおいて、HERV−Kファミリーは、レトロウイルスタンパク質の完全な又は部分的なコード能力を有する91のプロウイルス及び944のソロLTRを含む(60)。集合的に、HERV−Kプロウイルスは、感染力に必要な全てのレトロウイルス遺伝子のオープンリーディングフレームを維持し、たった3つのHERV−Kプロウイルスの間での潜在的組換えは、感染性レトロウイルスの産生を容易にする可能性がある(61)。しかしながら、ヒトにおけるがんの発生に対する、SCAR由来のレトロウイルス配列の重要な影響の新規の決定的な証拠は、感染性ウイルスの単離及びウイルス感染とがん発生との間の相関の確立を、必ずしも必要としない可能性がある。レトロウイルス配列の病理学的に重要な効果は、それらの生物学的活性の多くの異なるメカニズムに起因する可能性があり、以下の実験的証拠として証明することができる(62):
レトロウイルスの新規のがん特異的組込み部位の存在;
多くの異なる腫瘍における1個又は数個の宿主遺伝子の一貫した制御標的化;
レトロウイルス遺伝子(env;rec;np9)のタンパク質産物の発がん作用;
新規のスプライスドナー又はアクセプター部位、代替プロモーター、及び転写制御部位の寄与に起因する宿主遺伝子の発現に対する標的化された制御効果。
加えて、同じ及び/又は異なる染色体上の複数のSCARの配列の存在は、許容状態のクロマチンという文脈の中で、ゲノム遺伝子座間の組換え事象に起因する染色体再配列を容易にする可能性が高い。
本分析は、分化細胞におけるサイレンシングされたSCARの遺伝子座のエピジェネティックな活性化が、幹性制御ネットワークに関与することによって、細胞においてがん感受性状態を確立し得ることを示唆している。後続の突然変異誘発及びがんドライバー遺伝子の選択が、SCAR活性化幹性ネットワークを有する細胞において起こると主張することは妥当だと思われ、このことは、なぜ、高信頼性がんドライバー及びCOSMIC遺伝子のほぼ3分の2が、hESCにおいてSCARによって制御されるようであるのかを説明するだろう(上記参照)。この仮説の中心的な仮定は、がん幹細胞の前駆体としての機能を果たし得る、SCAR活性化幹性ネットワークを有する前がん状態の分化細胞の存在を予測し、その出現は、臨床的に難治性の悪性腫瘍の発生、腫瘍成長、がん進行、及び転移を、続いて促進するだろう。
材料及び方法
データソース及び分析プロトコル
この分析のために、もっぱら公的に利用可能なデータセット及びリソースを使用した。また、方法論的アプローチ及び計算パイプラインが、霊長類特異的遺伝子及びヒト特異的制御遺伝子座の発見について妥当性確認された[3;63−68]。shRNA実験を使用してhESCにおいて同定されたHERVH/LBP9制御遺伝子[19]を含む、SCAR関連幹性ゲノムネットワークを含む個々の遺伝要素を、転写活性のあるSCAR遺伝子座[12;16−20]、宿主/ウイルスキメラ転写物[18−20]、及びSCARがhESCゲノムに蒔かれたヒト特異的転写因子結合部位(TFBS)[3]を報告している、最近公表された寄稿から得た。
がんゲノムアトラス(TCGA)プロジェクトのウェブベースツールの最新のベータ版、UCSC Xena(http://xena.ucsc.edu/)、関連する臨床データ、及び数千の腫瘍サンプルにおいて同定された複数の機能的がんゲノミクスのエンドポイントを利用して、SCAR関連幹性ゲノムネットワークの個々の遺伝要素の遺伝子発現、体細胞非サイレント突然変異、及び遺伝子コピー数の臨床的に関連するパターンを、調査、分析、及び可視化した。これは、1万2千を超える注釈付きの臨床腫瘍サンプルの包括的な機能的がんゲノミクスデータセット(https://genomecancer.soe.ucsc.edu/proj/site/xena/datapages/)を問い合わせることによって行った。12のTCGAコホート(12の異なるがんタイプのPANCAN12研究)にわたる5,158の臨床サンプル及び全てのTCGAコホートにわたる12,088臨床サンプルを含む、2つのTCGA汎がんデータベース(https://genomecancer.soe.ucsc.edu/proj/site/xena/datapages/)を問い合わせることによって、高いがんによる死亡の可能性に関連する遺伝子発現、体細胞非サイレント突然変異、及びコピー数変化の汎がんシグネチャーを同定した。
配列保存分析は、カリフォルニア大学サンタクルーズ校(UCSC)のLiftOverアルゴリズムに基づいている。このLiftOverアルゴリズムは、0.95のMinMatch及びデフォルト設定の他の検索パラメーターで予め計算された全ゲノムBLASTZアラインメントの連鎖ファイルを使用して、ヒトブロックの座標を対応する非ヒトゲノムに変換するためのものである(http://genome.ucsc.edu/cgi-bin/hgLiftOver)。ヒトクエリーについての、LiftOverアルゴリズムによるBLASTZアラインメントの抽出は、LiftOverアウトプット「新たに欠失した(Deleted in new)」を生成する。このアウトプットは、ヒト配列が、所与の非ヒトゲノム中のいずれの鎖とも交差しないことを示す。これは、対象ゲノム中にクエリー配列が存在しないことを示し、非ヒト参照ゲノム中のヒト配列の存在又は非存在を推測するために使用された。BLASTアルゴリズム及び対応する参照ゲノムデータベースの最新のリリースを使用して、2013年4月から2015年10月までの期間、ヒト特異的制御配列を、それらの同一性及びゲノム特徴の妥当性確認をするために、手動でキュレートした。
宿主遺伝子に対するSCARの推定上の機能的に重要な制御効果の考察は、部分的に、対応する候補調節エレメント及び標的遺伝子のゲノムワイドな近接配置分析の結果に基づいた。近接配置分析中の近接の定量限界は、いくつかの測定基準に基づいて規定された。測定基準の1つは、hESCにおいて分析された調節タンパク質の50%の最も近い標的に対して実験的に規定された距離よりも、推定上の標的タンパク質コード遺伝子又はlncRNA遺伝子に近く、ヒト特異的制御配列を配置するゲノム座標を使用して規定された[69]。関心の対象となる各遺伝子について、hESCにおいてタンパク質コードTFによって制御される最も近い標的遺伝子までの距離の平均値より小さい、HSGRLと推定上の標的遺伝子との間のゲノム距離で、特定のHSGRLを同定し表にした。Guttmanらによって報告された[69]hESCにおけるTFの最も近い標的遺伝子までの距離に基づいて、タンパク質コード標的遺伝子及びlncRNA標的遺伝子についての対応する平均値を計算した。加えて、同じトポロジカル関連ドメイン(TAD)の境界内での共局在及び251の新皮質/前頭前皮質関連遺伝子の近くに位置するヒト特異的NANOG結合部位(HSNBS)の配置富化パターンに基づいて、近接配置測定基準を規定した[70]。HSNBSの配置富化分析は、1.5Mb未満のゲノム距離で最も顕著な富化を同定し、1.5Mbのゲノム距離で富化p値の鋭いピークを示した[70]。
この研究では、hESCゲノムにおいて同定された個々の調節エレメント及びクロマチン制御ドメインの包括的なデータベースを検討した。hESC中の3,127のトポロジカル関連ドメイン(TAD);6,823のhESC富化エンハンサー;hESC中の6,322の従来のエンハンサー及び684のスーパーエンハンサー(SE);mESC中の231のSE及び197のスーパーエンハンサードメイン(SED)のゲノム座標が、既に公表されている寄稿において報告された[2;71−74]。hESC中のNANOG−、POU5F1−、及びCTCF−結合部位及びヒト特異的TFBSの種特異的データセットが既に報告されており[3;4]、公的に利用可能である。規定されるゲノム領域の細胞タイプ特異的転写活性の可視化及び分析のために、RNA−SeqデータセットをUCSCデータリポジトリサイトから取り出した(http://genome.ucsc.edu/;[75])。単一塩基分解能でのヒトメチロームのゲノムワイドなマップが既に報告されており[76;77]、公的に利用可能である(http://neomorph.salk.edu/human_methylome)。可視化及び分析のためのヒストン修飾及び転写因子クロマチン免疫沈降シーケンス(ChIP−Seq)データセットを、UCSCデータリポジトリサイト(http://genome.ucsc.edu/;[78])から得た。ペアエンドタグシーケンシングによるクロマチン組込み解析(ChIA−PET)法によって決定された、RNAポリメラーゼII(PII)結合部位のゲノム座標を、MCF7及びK562ヒト細胞株について構築した飽和ライブラリーから得た[79]。選択されたヒト染色体の1−Mb及び1−kbの連続するセグメントあたりのタンパク質特異的結合事象の数を定量化し、得られた結合部位密度分布を可視化のためにプロットすることによって、染色体の所与のセグメントに対するTF結合の密度を推定した。WebLogoアルゴリズム(http://weblogo.berkeley.edu/logo.cgi)を使用して、複数の配列アラインメントの可視化を実施した。コンセンサスTF結合部位モチーフロゴは既報である[4;80;81]。
3人のネアンデルタール人、12人の新人、及び41,000年前のデニソワ人ゲノムの個々のゲノムにおけるHSGRLの保存の評価[82;83]を、個々のゲノム及びヒトゲノム参照データベース(http://genome.ucsc.edu/Neandertal/)から取り出した対応する配列の直接比較によって実施した。
タンパク質BLASTアルゴリズム(http://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&BLAST_PROGRAMS=blastp&PAGE_TYPE=BlastSearch&SHOW_DEFAULTS=on&LINK_LOC=blasthome)及び保存タンパク質ドメインの同定及び可視化のための関連のウェブベースツール(http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi?RID=3HZ5BMES01R&mode=all)(これらは、他の場所で詳細に説明された[84、85])を使用して、ヒト特異的キメラ転写物のヌクレオチド配列をアミノ酸配列に翻訳し、タンパク質アラインメント分析に供した。
米国における年齢調整されたがん発生率及び死亡率を、疾病管理予防センター(CDC)米国がん統計(USCS)報告から得た:
米国がん統計ワーキンググループ。米国がん統計:1999−2012年の発生率及び死亡率のウェブベースの報告。アトランタ:米国保健福祉省、疾患管理予防センター及び国立がん研究所;2015年。www.cdc.gov/uscsで利用可能。
公的に利用可能なデータセットの統計的分析
エラー率推定、バックグラウンド及び技術的ノイズの測定及びフィルタリング、特徴ピークコール、特徴選択、参照ヒトゲノムの対応する構築物へのゲノム座標の割り当て、及びデータ可視化を含む、公的に利用可能なゲノムデータセットの全ての統計的分析は、対応するデータ可視化トラック(http://genome.ucsc.edu/及びhttp://xena.ucsc.edu/)にリンクされている元の刊行物及び関連する参考文献において報告されているとおりに正確に実施した。統計的分析の任意の変更又は新規の要素は、結果の対応するセクションにおいて説明している。ピアソン相関係数の統計的有意性は、GraphPad Prismバージョン6.00ソフトウェアを使用して決定した。群間の事象の数の差異の有意性は、両側のフィッシャー直接検定及びカイ二乗検定を使用して計算し、事象間の重複の有意性は、超幾何分布検定を使用して決定した[86]。
参考文献
1. Santoni, F.A., Guerra, J., and Luban, J. HERV-H RNA is abundant in human embryonic stem cells and a precise marker for pluripotency. Retrovirology 2012; 9: 111.
2. Xie W, Schultz MD, Lister R, Hou Z, Rajagopal N, Ray P, Whitaker JW, Tian S, Hawkins RD, Leung D, Yang H, Wang T, Lee AY, Swanson SA, Zhang J, Zhu Y, Kim A, Nery JR, Urich MA, Kuan S, Yen CA, Klugman S, Yu P, Suknuntha K, Propson NE, Chen H, Edsall LE, Wagner U, Li Y, Ye Z, Kulkarni A, Xuan Z, Chung WY, Chi NC, Antosiewicz-Bourget JE, Slukvin I, Stewart R, Zhang MQ, Wang W, Thomson JA, Ecker JR, Ren B. Epigenomic analysis of multilineage differentiation of human embryonic stem cells. Cell 2013. 153: 1134-1148.
3. Glinsky, GV. Transposable Elements and DNA Methylation Create in Embryonic Stem Cells Human-Specific Regulatory Sequences Associated with Distal Enhancers and Noncoding RNAs. Genome Biol Evol. 2015; 7: 1432-54.
4. Kunarso, G, Chia, NY, Jeyakani, J, Hwang, C, Lu, ., Chan, YS, Ng, HH, and Bourque, G. Transposable elements have rewired the core regulatory network of human embryonic stem cells. Nat Genet. 2010; 42: 631-634.
5. Kelley, D, and Rinn, J. Transposable elements reveal a stem cell-specific class of long noncoding RNAs. Genome Biol. 2012; 13: R107.
6. Glinsky GV. Endogenous human stem cell-associated retroviruses. BioRxiv 2015; doi: http://dx.doi.org/10.1101/024273
7. Glinsky GV. SCARs: endogenous human stem cell-associated retroviruses and therapy-resistant malignant tumors. arXiv preprint 2015; arXiv:1508.02022 http://arxiv.org/abs/1508.02022
8. Glinsky GV. Viruses, stemness, embryogenesis, and cancer: a miracle leap toward molecular definition of novel oncotargets for therapy-resistant malignant tumors? Oncoscience 2015; 2: 751-754.
9. Glinsky GV. Activation of endogenous human Stem Cell-Associated Retroviruses and therapy-resistant phenotypes of malignant tumors. 2016. In revision.
10. Smith ZD, Chan MM, Humm KC, Karnik R, Mekhoubad S, Regev A, Eggan K, Meissner A. DNA methylation dynamics of the human preimplantation embryo. Nature 2014; 511: 611-615.
11. Fort A, Hashimoto K, Yamada D, Salimullah M, Keya CA, Saxena A, Bonetti A, Voineagu I, Bertin N, Kratz A, Noro Y, Wong CH, de Hoon M, Andersson R, Sandelin A, Suzuki H, Wei CL, Koseki H; FANTOM Consortium, Hasegawa Y, Forrest AR, Carninci P. Deep transcriptome profiling of mammalian stem cells supports a regulatory role for retrotransposons in pluripotency maintenance. Nature Genet. 2-14; 46: 558-566.
12. Lu X, Sachs F, Ramsay L, Jacques PE, Goke J, Bourque G, Ng HH. The retrovirus HERVH is a long noncoding RNA required for human embryonic stem cell identity. Nat Struct Mol Biol. 2014; 21:423-425.
13. Ohnuki M, Tanabe K1, Sutou K, Teramoto I, Sawamura Y, Narita M, Nakamura M, Tokunaga Y, Nakamura M, Watanabe A, Yamanaka S, Takahashi K. Dynamic regulation of human endogenous retroviruses mediates factor-induced reprogramming and differentiation potential. Proc Natl Acad Sci USA. 2014. 111:12426-31.
14. Koyanagi-Aoi M, Ohnuki M, Takahashi K, Okita K, Noma H, Sawamura Y, Teramoto I, Narita M, Sato Y, Ichisaka T, Amano N, Watanabe A, Morizane A, Yamada Y, Sato T, Takahashi J, Yamanaka S. Differentiation-defective phenotypes revealed by large-scale analyses of human pluripotent stem cells. Proc Natl Acad Sci USA. 2013; 110: 20569-74.
15. Marchetto MC, Narvaiza I, Denli AM, Benner C, Lazzarini TA, Nathanson JL, Paquola AC, Desai KN, Herai RH, Weitzman MD, Yeo GW, Muotri AR, Gage FH. (2013). Differential LINE-1 regulation in pluripotent stem cells of humans and other great apes. Nature 503: 525-529.
16. Xue Z, Huang K, Cai C, Cai L, Jiang CY, Feng Y, Liu Z, Zeng Q, Cheng L, Sun YE, Liu JY, Horvath S, Fan G. Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature 2013; 500: 593-597.
17. Yan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J, Huang J, Li M, Wu X, Wen L, Lao K, Li R, Qiao J, Tang F. Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol 2013; 20: 1131-1139.
18. Goke J, Lu X, Chan YS, Ng HH, Ly LH, Sachs F, Szczerbinska I. Dynamic transcription of distinct classes of endogenous retroviral elements marks specific populations of early human embryonic cells. Cell Stem Cell 2015; 16: 135-141.
19. Wang J, Xie G, Singh M, Ghanbarian AT, Rasko T, Szvetnik A, Cai H, Besser D, Prigione A, Fuchs NV, Schumann GG, Chen W, Lorincz MC, Ivics Z, Hurst LD, Izsvak Z. Primate-specific endogenous retrovirus-driven transcription defines naive-like stem cells. Nature 2014; 516: 405-9.
20. Grow EJ, Flynn RA, Chavez SL, Bayless NL, Wossidlo M, Wesche DJ, Martin L, Ware CB, Blish CA, Chang HY, Pera RA, Wysocka J. Intrinsic retroviral reactivation in human preimplantation embryos and pluripotent cells. Nature 2015; 522: 221-5.
21. Robbez-Masson L, Rowe HM. Retrotransposons shape species-specific embryonic stem cell gene expression. Retrovirology 2015; 12: 45.
22. Tamborero D1, Gonzalez-Perez A, Perez-Llamas C, Deu-Pons J, Kandoth C, Reimand J, Lawrence MS, Getz G, Bader GD, Ding L, Lopez-Bigas N. Comprehensive identification of mutational cancer driver genes across 12 tumor types. Sci Rep. 2013; 3: 2650.
23. Hoadley KA, Yau C, Wolf DM, Cherniack AD, Tamborero D, Ng S, Leiserson MD, Niu B, McLellan MD, Uzunangelov V, Zhang J, Kandoth C, Akbani R, Shen H, Omberg L, Chu A, Margolin AA, Van't Veer LJ, Lopez-Bigas N, Laird PW, Raphael BJ, Ding L, Robertson AG, Byers LA, Mills GB, Weinstein JN, Van Waes C, Chen Z, Collisson EA; Cancer Genome Atlas Research Network, Benz CC, Perou CM, Stuart JM. Multiplatform analysis of 12 cancer types reveals molecular classification within and across tissues of origin. Cell 2014; 158: 929-44.
24. Yu, X. and Gabriel, A. Patching broken chromosomes with extranuclear cellular DNA. Mol. Cell 1999; 4: 873-881.
25. Lin, Y. and Waldman, A.S. Promiscuous patching of broken chromosomes in mammalian cells with extrachromosomal DNA. Nucleic Acids Res. 2001; 29: 3975-3981.
26. Teng, S.C., Kim, B. and Gabriel, A. Retrotransposon reverse transcriptase-mediated repair of chromosomal breaks. Nature 1996; 383: 641-644.
27. Morrish, T.A., Gilbert, N., Myers, J.S., Vincent, B.J., Stamato, T.D., Taccioli, G.E., Batzer, M.A. and Moran, J.V. DNA repair mediated by endonuclease-independent LINE-1 retrotransposition. Nat. Genet. 2002; 31: 159-165.
28. Morrish TA, Garcia-Perez JL, Stamato TD, Taccioli GE, Sekiguchi J, Moran JV. Endonuclease-independent LINE-1 retrotransposition at mammalian telomeres. Nature. 2007; 446: 208-12.
29. Ichiyanagi, K., Nakajima, R., Kajikawa, M. and Okada, N. (2007) Novel retrotransposon analysis reveals multiple mobility pathways dictated by hosts. Genome Res. 2007; 17: 33-41.
30. Sen, S.K., Huang, C.T., Han, K., Batzer, M.A. Endonuclease-independent insertion provides an alternative pathway for L1 retrotransposition in the human genome. Nucleic Acids Res. 2007; 35: 3741-3751.
31. Srikanta D, Sen SK, Huang CT, Conlin EM, Rhodes RM, et al. An alternative pathway for Alu 63 retrotransposition suggests a role in DNA double strand break repair. Genomics 2009; 93: 205-212.
32. Shin W, Lee J, Son S-Y, Ahn K, Kim H-S, Han, K. Human-specific HERVK insertion causes genomic variations in the human genome. PLoS ONE 2013; 8: e60605.
33. Nussenzweig A, Nussenzweig MC. A backup DNA repair pathway moves to the forefront. Cell. 2007; 131: 223-225.
34. Iliakis G. Backup pathways of NHEJ in cells of higher eukaryotes: cell cycle dependence. Radiother Oncol. 2009; 92: 310-315.
35. Bogomazova AN, Lagarkova MA, Tskhovrebova LV, Shutova MV, Kiselev SL. Error-prone nonhomologous end joining repair operates in human pluripotent stem cells during late G2. Aging (Albany NY). 2011; 3: 584-96.
36. Fan J, Robert C, Jang YY, Liu H, Sharkis S, Baylin SB, Rassool FV. Human induced pluripotent cells resemble embryonic stem cells demonstrating enhanced levels of DNA repair and efficacy of nonhomologous end-joining. Mutat Res. 2011; 713: 8-17.
37. Glinsky GV, Glinskii AB, Berezovskaya O. Microarray analysis identifies a death-from-cancer signature predicting therapy failure in patients with multiple types of cancer. Journal of Clinical Investigation 2005; 115: 1503 - 21.
38. Glinsky GV. Death-from-cancer signatures and stem cell contribution to metastatic cancer. Cell Cycle 2005; 4: 1171 - 5.
39. Glinsky, GV. Genomic models of metastatic cancer: Functional analysis of death-from-cancer signature genes reveals aneuploid, anoikis-resistant, metastasis-enabling phenotype with altered cell cycle control and activated Polycomb Group (PcG) protein chromatin silencing pathway. Cell Cycle, 2006; 5: 1208-1216.
40. Berezovska, OP, Glinskii, AB, Yang, Z, Li, X-M, Hoffman, RM, Glinsky, GV. Essential role of the Polycomb Group (PcG) protein chromatin silencing pathway in metastatic prostate cancer. Cell Cycle, 2006; 5: 1886-1901.
41. Glinskii AB, Smith BA, Jiang P, Li XM, Yang M, Hoffman RM, Glinsky GV. Viable circulating metastatic cells produced in orthotopic but not ectopic prostate cancer models. Cancer Res. 2003; 63: 4239-43.
42. Berezovskaya O, Schimmer AD, Glinskii AB, Pinilla C, Hoffman RM, Reed JC, Glinsky GV. Increased expression of apoptosis inhibitor protein XIAP contributes to anoikis resistance of circulating human prostate cancer metastasis precursor cells. Cancer Res. 2005; 65: 2378-86.
43. Glinsky GV, Glinskii AB, Berezovskaya O, Smith BA, Jiang P, Li XM, Yang M, Hoffman RM. Dual-color-coded imaging of viable circulating prostate carcinoma cells reveals genetic exchange between tumor cells in vivo, contributing to highly metastatic phenotypes. Cell Cycle. 2006; 5: 191-7.
44. Holt, S., Glinsky, V.V., Ivanova, A.B., Glinsky, G.V. Resistance to apoptosis in human cells conferred by telomerase function and telomere stability. Molecular Carcinogenesis 1999; 25: 241-248.
45. Glinsky, G.V., Glinsky, V.V., Ivanova, A.B., Hueser, C.N. Apoptosis and metastasis: Increased apoptosis resistance of metastatic cancer cells is associated with the profound deficiency of apoptosis execution mechanisms. Cancer Letters 1997; 115: 185-193.
46. Glinsky, G.V. Apoptosis in metastatic cancer cells. Crit. Rev. Oncol/Hemat. 1997; 25: 175-186.
47. Glinsky, GV, Glinsky, VV. Apoptosis and metastasis: A superior resistance of metastatic cancer cells to programmed cell death. Cancer Letters 1996; 101: 43-51.
48. Glinsky GV. Stem cell origin of death-from-cancer phenotypes of human prostate and breast cancers. Stem Cells Reviews 2007; 3: 79-93.
49. Glinsky GV. "Stemness" genomics law governs clinical behavior of human cancer: Implications for decision making in disease management. Journal of Clinical Oncology 2008; 26:2 846-53.
50. Glinsky GV, Berezovska O, Glinskii A. Genetic signatures of regulatory circuitry of embryonic stem cells (ESC) identify therapy-resistant phenotypes in cancer patients diagnosed with multiple types of epithelial malignancies. Cancer Research 2007; 67 (9 Supplement):1272.
51. Glinskii A, Berezovskaya O, Sidorenko A, Glinsky G. Stemness pathways define therapy-resistant phenotypes of human cancers. Clinical Cancer Research 2008; 14 (15 Supplement):B38.
52. Schwartzberg P, Colicelli J, Goff SP. Recombination between a defective retrovirus and homologous sequences in host DNA: reversion by patch repair. J Virol. 1985; 53: 719-26.
53. McClure HM. Tumors in nonhuman primates: observations during a six-year period in the Yerkes primate center colony. Am J Phys Anthropol. 1973; 38:425-429.
54. Seibold HR, Wolf RH. Neoplasms and proliferative lesions in 1065 nonhuman primate necropsies. Lab Anim Sci. 1973; 23:533-539.
55. Beniashvili DS. An overview of the world literature on spontaneous tumors in nonhuman primates. J Med Primatol. 1989; 18:423-437.
56. Scott, G.B.D. 1992. Comparative primate pathology. Oxford University Press, New York, NY.
57. Waters DJ, Sakr WA, Hayden DW, Lang CM, McKinney L, Murphy GP, Radinsky R, Ramoner R, Richardson RC, Tindall DJ. Workgroup 4: spontaneous prostate carcinoma in dogs and nonhuman primates. Prostate. 1998; 36: 64-67.
58. Simmons HA, Mattison JA. The incidence of spontaneous neoplasia in two populations of captive rhesus macaques (Macaca mulatta). Antioxid Redox Signal. 2011; 14: 221-7.
59. Gemmell, P., Hein, J., Katzourakis, A. Orthologous endogenous retroviruses exhibit directional selection since the chimp-human split. Retrovirology 2015; 12: 52.
60. Subramanian, R.P., Wildschutte, J.H., Russo, C., Coffin, J.M. Identification, characterization, and comparative genomic distribution of the HERV-K (HML-2) group of human endogenous retroviruses. Retrovirology 2011; 8: 90.
61. Hohn, O., Hanke, K., Bannert, N. HERV-K(HML-2), the best preserved family of HERVs: Endogenization, expression, and implications in health and disease. Front Oncol 2013; 3: 246.
62. Bhardwaj, N., Coffin, J.M. Endogenous Retroviruses and Human Cancer: Is There Anything to the Rumors? Cell Host & Microbes 2014; 15: 255-250.
63. Kent, WJ. BLAT - the BLAST-like alignment tool. Genome Res. 2002; 12: 656-664.
64. Schwartz, S., Kent, W.J., Smit, A., Zhang, Z., Baertsch, R., Hardison, R.C., Haussler, D., and Miller, W. Human-mouse alignments with BLASTZ. Genome Res. 2003; 13: 103-107.
65. Tay, S.K., Blythe, J., and Lipovich, L. Global discovery of primate-specific genes in the human genome. Proc. Natl. Acad. Sci. USA 2009; 106: 12019-12024.
66. Capra, J.A., Erwin, G.D., McKinsey, G., Rubenstein, J.L., Pollard, K.S. Many human accelerated regions are developmental enhancers. Philos Trans R Soc Lond B Biol Sci. 2013; 368 (1632): 20130025.
67. Marnetto D, Molineris I, Grassi E, Provero P. Genome-wide identification and characterization of fixed human-specific regulatory regions. Am J Hum Genet 2014; 95: 39-48.
68. Gittelman RM, Hun E, Ay F, Madeoy J, Pennacchio L, Noble WS, Hawkins RD, Akey JM. 2015. Comprehensive identification and analysis of human accelerated regulatory DNA. Genome Res 2015; 25: 1245-55.
69. Guttman, M., Donaghey, J., Carey, B.W., Garber, M., Grenier, J.K., Munson, G., Young, G., Lucas, A.B., Ach, R., Bruhn, L., Yang, X., Amit, I., Meissner, A., Regev, A., Rinn, J.L., Root, D.E., and Lander, E.S. lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature 2011; 477: 295-300.
70. Glinsky, GV. Rapidly evolving in humans topologically associating domains. 2015. arXiv:1507.05368 .
71. Dixon, J.R., Selvaraj, S., Yue, F., Kim, A., Li, Y., Shen, Y., Hu, M., Liu, J.S., and Ren, B. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature 2012; 485: 376-380.
72. Dowen J.M., Fan Z.P., Hnisz D., Ren G., Abraham B.J., Zhang L.N., Weintraub A.S., Schuijers J., Lee T.I., Zhao K., Young RA. Control of cell identity genes occurs in insulated neighborhoods in mammalian chromosomes. Cell 2014; 159: 374-387.
73. Hnisz, D., Abraham, B.J., Lee, T.I., Lau, A., Saint-Andre´, V., Sigova, A.A., Hoke, H.A., and Young, RA. Super-enhancers in the control of cell identity and disease. Cell 2013; 155: 934-947.
74. Whyte, W.A., Orlando, D.A., Hnisz, D., Abraham, B.J., Lin, C.Y., Kagey, M.H., Rahl, P.B., Lee, T.I., and Young, RA. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 2013; 153: 307-319.
75. Meyer, L.R., Zweig, A.S., Hinrichs, A.S., Karolchik, D., Kuhn, R.M., Wong, M., Sloan, C.A., Rosenbloom, K.R., Roe, G., Rhead, B., Raney, B.J., Pohl, A., Malladi, V.S., Li, C.H., Lee, B.T., Learned, K., Kirkup, V., Hsu, F., Heitner, S., Harte, R.A., Haeussler, M., Guruvadoo, L., Goldman, M., Giardine, B.M., Fujita, P.A., Dreszer, T.R., Diekhans, M., Cline, M.S., Clawson, H., Barber, G.P., Haussler, D., and Kent, W.J. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013; 41: D64-69.
76. Lister, R., Pelizzola, M., Dowen, R.H., Hawkins, R.D., Hon, G., Tonti-Filippini, J., Nery, J.R., Lee, L., Ye, Z., Ngo, Q.M., Edsall, L., Antosiewicz-Bourget, J., Stewart, R., Ruotti, V., Millar, A.H., Thomson, J.A., Ren, B., and Ecker, JR. Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 2009; 462: 315-322.
77. Lister R, Mukamel EA, Nery JR, Urich M, Puddifoot CA, Johnson ND, Lucero J, Huang Y, Dwork AJ, Schultz MD, Yu M, Tonti-Filippini J, Heyn H, Hu S, Wu JC, Rao A, Esteller M, He C, Haghighi FG, Sejnowski TJ, Behrens MM, Ecker JR. Global epigenomic reconfiguration during mammalian brain development. Science 2013; 341: 1237905.
78. Rosenbloom, K.R., Sloan, C.A., Malladi, V.S., Dreszer, T.R., Learned, K., Kirkup, V.M., Wong, M.C., Maddren, M., Fang, R., Heitner, S.G., Lee, B.T., Barber, G.P., Harte, R.A., Diekhans, M., Long, J.C., Wilder, S.P., Zweig, A.S., Karolchik, D., Kuhn, R.M., Haussler, D., and Kent, WJ. ENCODE data in the UCSC Genome Browser: year 5 update. Nucleic Acids Res 2013; 41: D56-63.
79. Li, G., Ruan, X., Auerbach, R.K., Sandhu, K.S., Zheng, M., Wang, P., Poh, H.M., Goh, Y., Lim, J., Zhang, J., Sim, H.S., Peh, S.Q., Mulawadi, F.H., Ong, C.T., Orlov, Y.L., Hong, S., Zhang, Z., Landt, S., Raha, D., Euskirchen, G., Wei, C.L., Ge, W., Wang, H., Davis, C., Fisher-Aylor, K.I., Mortazavi, A., Gerstein, M., Gingeras, T., Wold, B., Sun, Y., Fullwood, M.J., Cheung, E., Liu, E., Sung, W.K., Snyder, M., and Ruan, Y. Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation. Cell 2012; 148: 84-98.
80. Wang, J., Zhuang, J., Iyer, S., Lin, X., Whitfield, T.W., Greven, M.C., Pierce, B.G., Dong, X., Kundaje, A., Cheng, Y., Rando, O.J., Birney, E., Myers, R.M., Noble, W.S., Snyder, M., and Weng, Z. Sequence features and chromatin structure around the genomic regions bound by 119 human transcription factors. Genome Res. 2012; 22: 1798-1812.
81. Ernst, J., and Kellis, M. 2013. Interplay between chromatin state, regulator binding, and regulatory motifs in six human cell types. Genome Res. 2013; 23: 1142-1154.
82. Reich, D., Green, R.E., Kircher, M., Krause, J., Patterson, N., Durand, E.Y., Viola, B., Briggs, A.W., Stenzel, U., Johnson, P.L., Maricic, T., Good, J.M., Marques-Bonet, T., Alkan, C., Fu, Q., Mallick, S., Li, H., Meyer, M., Eichler, E.E., Stoneking, M., Richards, M., Talamo, S., Shunkov, M.V., Derevianko, A.P., Hublin, J.J., Kelso, J., Slatkin, M., Paabo, S. Genetic history of an archaic hominin group from Denisova Cave in Siberia. Nature 2010; 468: 053-1060.
83. Meyer, M., Kircher, M., Gansauge, M.T., Li, H., Racimo, F., Mallick, S., Schraiber, J.G., Jay, F., Prufer, K., de Filippo, C., Sudmant, P.H., Alkan, C., Fu, Q., Do, R., Rohland, N., Tandon, A., Siebauer, M., Green, R.E., Bryc, K., Briggs, A.W., Stenzel, U., Dabney, J., Shendure, J., Kitzman, J., Hammer, M.F., Shunkov, M.V., Derevianko, A.P., Patterson, N., Andres, A.M., Eichler, E.E., Slatkin, M., Reich, D., Kelso, J., Paabo, S. A high-coverage genome sequence from an archaic Denisovan individual. Science 2012; 338: 222-226.
84. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011; 39: D225-9.
85. Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S2, Chitsaz F, Geer LY, Geer RC, He J, Gwadz M, Hurwitz DI, Lanczycki CJ, Lu F, Marchler GH, Song JS, Thanki N, Wang Z, Yamashita RA, Zhang D, Zheng C, Bryant SH. CDD: NCBI's conserved domain database. Nucleic Acids Res. 2015; 43: D222-6.
86. Tavazoie, S., Hughes, J.D., Campbell, M.J., Cho, R.J., and Church, GM. 1999. Systematic determination of genetic network architecture. Nat. Genet.1999; 22: 281-285.
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009

Figure 0006963009

Figure 0006963009
表4.図1のデータ
Figure 0006963009
表5.図4のデータ
Figure 0006963009
Figure 0006963009
表6.図5のデータ
Figure 0006963009
表7.図6のデータ
Figure 0006963009
表8.図7のデータ
Figure 0006963009
表9.図8のデータ(タンパク質P値)
Figure 0006963009
表10−14(データセットS2)は、チンパンジー、ボノボ、ゴリラ、オランウータン、テナガザル、及びアカゲザルのゲノムを含む、17種の霊長類のゲノムの配列との、ヒトゲノム配列の比較中の直接及び相互配列アラインメント変換失敗に基づいて規定される、ヒト特異的SCAR遺伝子座の説明を含む。表10−Xはまた、各SCAR遺伝子座について、17種の霊長類のゲノムへの配列アラインメントによって規定される、祖先DNAのヒト特異的欠失のサイズを示す。
表10.251b.c.失敗
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
表11.チンパンジー及びボノボのゲノムからヒトゲノムへの相互アラインメントの失敗に基づいて規定される、ヒト特異的SCAR
Figure 0006963009
Figure 0006963009
Figure 0006963009
表12.128hervhから。
Figure 0006963009
Figure 0006963009
表13.HERVH由来lincRNA。
Figure 0006963009
Figure 0006963009
表14.43のヒト特異的組込み部位
Figure 0006963009
表15.SNMの10のデータセット。
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
Figure 0006963009
表16.SNM P値。
Figure 0006963009
Figure 0006963009
表17.SNMのパーセント。
Figure 0006963009
注:表4−9は「データセットS1」であり、表10−14は「データセットS2」であり、表15−17は「データセットS3」である。
パラグラフ1:被験体においてがんを診断する又はがん治療成績を予測するための方法であって、内在性ヒト幹細胞関連レトロウイルス(SCAR)のプロテオミクスシグネチャー経路を示す1又は複数の入力及びゲノムシグネチャー経路を示す1又は複数の入力に応答して標的マーカー情報を生成する工程;及び標的マーカー情報を有する生体サンプルの検出される発現レベル及び配列情報の比較に応答して異常オブジェクト情報を生成する工程を含む、方法。
一実施形態において、異常オブジェクト情報を生成する工程は、異常オブジェクト情報をクライアント装置、ユーザインターフェイスなどに表示する工程を含む。一実施形態において、異常オブジェクト情報を生成する工程は、異常オブジェクト情報をリモートネットワークで交換する工程を含む。異常オブジェクト情報の非限定的な例は、異常配列情報、異常発現レベル情報、発現レベルが標的閾値情報を上回っていること、複数の塩基の検出された位置付け、配列異常スコアなどを含む。
さらなる異常オブジェクト情報の非限定的な例は、生体被験体から得られたサンプルに由来する参照情報を比較することによって得られた閾値レベルを示す情報;標的マーカー情報を有する生体サンプルの配列を示す少なくとも1つの入力と発現レベルを示す少なくとも1つの入力の比較を示す情報;などを含む。
パラグラフ2:標的マーカー情報を生成する工程が、SCAR経路を示す1又は複数の入力に応答して標的マーカー情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ3:標的マーカー情報を生成する工程が、SCAR経路標的遺伝子を示す1又は複数の入力に応答して標的マーカー情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ4:標的マーカー情報を生成する工程が、ELF3;PCDH15;MALAT1;PTPN11;RB1;CHST6;NF1;VEZF1;TP53;SMAD4;KEAP1;STK11;PRX;ZNF28;IDH1;FEZ2;DPPA2;LPHN3;KIAA1244;EPHA7;EGFR;TLR4;DAB2IP;NOTCH1;GLUD2;DMD;KDM6A;KRAS;CDKN2A;DNMT3A;FLT3;NFE2L2;NPM1;MIR142;FOXL2;H3F3A;H3F3B;KMT2D;RNF43;TERT;ERBB2;PLCG1のうち1又は複数に関連する標的マーカー情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ5:標的マーカー情報を生成する工程が、mRNA、RNA、DNA、ペプチド又はタンパク質のうち1又は複数に関連する標的マーカー情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ6:標的マーカー情報を生成する工程が、PLCXD1、HKR1、ZNF283、ADA、AMACR+p63、ANK3、BCL2L1、BIRC5、BMI−1、BUB1、CCNB1、CCND1、CES1、CHAF1A、CRIP1、CRYAB、ESM1、EZH2、FGFR2、FOS、Gbx2、HCFC1、IER3、ITPR1、JUNB、KLF6、KI67、KNTC2、MGC5466、Phc1、RNF2、Suz12、TCF2、TRAP100、USP22、Wnt5A及びZFP36のうち1又は複数に関連する標的マーカー情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ7:異常オブジェクト情報を生成する工程が、生体サンプルに関連する配列の質が1又は複数の参照配列と比較して異なる場合に異常配列情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ8:異常オブジェクト情報を生成する工程が、1又は複数の参照配列と比較して、生体サンプルに関連する配列全体内の複数の塩基の異なる位置付けを示す1又は複数の入力に応答して異常配列情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ9:異常オブジェクト情報を生成する工程が、1又は複数の参照配列と比較して、生体サンプルに関連する配列の異なる断片を示す1又は複数の入力に応答して異常配列情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ10:異常オブジェクト情報を生成する工程が、発現レベルが標的閾値を超えた場合を示す1又は複数の入力に応答して異常発現レベル情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ11:異常オブジェクト情報を生成する工程が、検出された発現レベルが標的閾値を上回る場合に発現レベル異常スコアを決定する工程を含む、パラグラフ1に記載の方法。
パラグラフ12:異常オブジェクト情報を生成する工程が、生体サンプルに関連する複数の塩基の検出された位置付けが1又は複数の参照配列と比較して異なる場合に、配列異常スコアを決定する工程を含む、パラグラフ1に記載の方法。
パラグラフ13:異常オブジェクト情報を生成する工程が、次世代シーケンシング、多色定量的免疫蛍光共局在分析、蛍光in situハイブリダイゼーション、及び定量RT−PCR分析からの1又は複数の入力に応答して配列異常スコアを決定する工程を含む、パラグラフ1に記載の方法。
パラグラフ14:異常オブジェクト情報を生成する工程が、既知の診断又は既知の治療後の臨床成績を有する生体被験体から得られたサンプルに由来する参照情報を比較することによって閾値レベルを決定する工程を含む、パラグラフ1に記載の方法。
パラグラフ15:パラグラフ14に記載の方法であって、少なくとも2つのマーカーの標的閾値係数を上回る異常な発現及び発現レベルを示す1又は複数の入力に応答して、がん治療有効性状態、がん治療進行、がん予後、がん診断を生成する工程をさらに含む、方法。
パラグラフ16:異常オブジェクト情報を生成する工程が、異常配列情報及びマーカー共発現レベル情報を生成する工程を含む、パラグラフ1に記載の方法。
パラグラフ17:パラグラフ1に記載の方法であって、
異常な配列及び閾値マーカー共発現レベルを示す1又は複数の入力に応答してがん治療有効性状態を生成する工程
をさらに含む、方法。
パラグラフ18:パラグラフ1に記載の方法であって、異常な配列及び閾値マーカー共発現レベルを示す1又は複数の入力に応答して、生体被験体におけるがんの存在又は非存在を示す情報を生成する工程をさらに含む、方法。
パラグラフ19:被験体においてがんを診断するための又はがん治療成績を予測するためのシステムであって、内在性ヒト幹細胞関連レトロウイルス(SCAR)のプロテオミクスシグネチャー経路を示す1又は複数の入力及びゲノムシグネチャー経路を示す1又は複数の入力に応答して、標的マーカー情報を生成するように構成される回路;及び標的マーカー情報を有する生体サンプルの配列を示す少なくとも1つの入力と発現レベルを示す少なくとも1つの入力の比較に応答して、異常オブジェクト情報を生成するように構成される回路を含む、システム。
パラグラフ20:パラグラフ19に記載のシステムであって、異常な配列及び閾値マーカー共発現レベルを示す1又は複数の入力に応答して、生体被験体におけるがんの存在又は非存在を示す情報を生成するように構成される回路をさらに含む、システム。
パラグラフ21:パラグラフ19に記載のシステムであって、少なくとも2つのマーカーの標的閾値係数を上回る異常な発現及び発現レベルを示す1又は複数の入力に応答して、がん治療有効性状態、がん治療進行、がん予後、がん診断を生成するように構成される回路をさらに含む、システム。
パラグラフ22:パラグラフ19に記載のシステムであって、異常な配列及び閾値マーカー共発現レベルを示す1又は複数の入力に応答して、がん治療有効性状態を生成するように構成される回路をさらに含む、システム。
パラグラフ23:がんを処置するためのシステムであって、がんを有すると診断された被験体における幹細胞関連レトロウイルス(SCAR)経路活性化に関連する情報を取得するように構成される回路;及びがんを有すると診断された被験体における幹細胞関連レトロウイルス(SCAR)経路活性化に関連する1又は複数の入力に応答して、単一の治療剤又は複数の治療剤の組み合わせを同定し、ユーザ固有の処置プロトコルを生成するように構成される回路を含む、システム。
パラグラフ24:被験体においてがんを診断する又はがん治療成績を予測するための方法であって、内在性ヒト幹細胞関連レトロウイルス(SCAR)のゲノム及びプロテオミクスシグネチャーに関与する経路に関連する1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在について生体サンプルを同時にスクリーニングする工程;生体サンプルに関連する配列を、配列の質が参照配列と比較して異なる場合に、異常であるとスコア付けする工程;及び検出された発現レベルが標的閾値係数を上回る場合に、生体サンプルに関連する発現レベルを異常であるとスコア付けする工程を含む、方法。一実施形態において、被験体においてがんを診断する又はがん治療成績を予測するための方法であって、内在性ヒト幹細胞関連レトロウイルス(SCAR)のゲノム及びプロテオミクスシグネチャーに関与する経路に関連する1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在の少なくとも1つについて生体サンプルをスクリーニングする工程;生体サンプルに関連する配列を、配列の質が参照配列と比較して異なる場合に、異常であるとスコア付けする工程;及び検出された発現レベルが標的閾値係数を上回る場合に、生体サンプルに関連する発現レベルを異常であるとスコア付けする工程を含む、方法。
パラグラフ25:内在性SCARのゲノム及びプロテオミクスシグネチャーに関与する経路に関連する1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在について生体サンプルを同時にスクリーニングする工程が、生体被験体におけるがん診断又はがん治療失敗の予後を示す1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在について、生体サンプルを同時にスクリーニングする工程を含む、パラグラフ24に記載の方法。
パラグラフ26:パラグラフ25に記載の方法であって、生体被験体におけるがん診断又はがん治療失敗の予後に関連する異常な配列又は異常な発現レベルを示す1又は複数の入力に応答してユーザ固有のがん治療プロトコルを生成する工程をさらに含む、方法。
パラグラフ27:内在性SCARのゲノム及びプロテオミクスシグネチャーに関与する経路に関連する1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在について生体サンプルを同時にスクリーニングする工程が、生体被験体におけるがん治療の進行を示す1又は複数の標的マーカーの異常な配列及び異常な発現レベルの存在について、生体サンプルを同時にスクリーニングする工程を含む、パラグラフ24に記載の方法。
パラグラフ28:パラグラフ27に記載の方法であって、生体被験体におけるがん治療の進行に関連する異常な配列又は異常な発現レベルを示す1又は複数の入力に応答して、ユーザ固有のがん治療プロトコルを生成する工程をさらに含む、方法。
パラグラフ29:検出閾値が、既知の診断又は既知の治療後の臨床成績を有する被験体から得られたサンプルの参照データベースにおける値と比較することによって決定されており、試験サンプルにおける、少なくとも1つであるが好ましくは2つ以上のマーカーの異常な発現レベルの存在及び2つ以上のそのようなマーカーの異常な発現の存在が、がん診断又はがん治療失敗の予後を示す、又は被験体におけるがん治療の進行を示す、パラグラフ24に記載の方法。
パラグラフ30:試験された各患者の成績データをサンプルの参照データベースに追加すること、及び自動化された及び/又は再帰的な方法で、手動で又はコンピュータによる方法を使用して、ローカルに、リモートサーバに、又はクラウドに保存されたデータを使用して、診断、予後、又は将来のがん治療の明細の正確さを連続して向上させることによって、検出閾値が連続して改良される、パラグラフ24に記載の方法。
パラグラフ31:前記サンプル表現型が、がん、非がん、再発、非再発、再燃、非再燃、浸潤性、非浸潤性、転移性、非転移性、限局性、腫瘍サイズ、腫瘍悪性度、Gleasonスコア、生存予後、リンパ節の状態、腫瘍の病期、分化の程度、年齢、ホルモンレセプターの状態、腫瘍抗原レベル(PSAレベル、PSMAレベル、サバイビンレベル、がん胎児性タンパク質レベル、精巣抗原レベルを含むが、これらに限定されない)、組織学的タイプ、免疫細胞のレベル、表現型及び遺伝子型、及び活性化状態、及び無病生存率からなる群より選択される、パラグラフ24に記載の方法。
パラグラフ32:前記閾値係数が、絶対値≧0.5を有する、パラグラフ24に記載の方法。
パラグラフ33:前記閾値係数が、絶対値≧0.6を有する、パラグラフ24に記載の方法。
パラグラフ34:前記閾値係数が、絶対値≧0.7を有する、パラグラフ24に記載の方法。
パラグラフ35:前記閾値係数が、絶対値≧0.8を有する、パラグラフ24に記載の方法。
パラグラフ36:前記閾値係数が、絶対値≧0.9を有する、パラグラフ24に記載の方法。
パラグラフ37:前記閾値係数が、絶対値≧0.95を有する、パラグラフ24に記載の方法。
パラグラフ38:前記閾値係数が、絶対値≧0.99を有する、パラグラフ24に記載の方法。
パラグラフ39:前記閾値係数が、絶対値≧0.995を有する、パラグラフ24に記載の方法。
パラグラフ40:前記閾値係数が、絶対値≧0.999を有する、パラグラフ24に記載の方法。
パラグラフ41:サンプル表現型を分類するための検出閾値を決定する方法であって、マーカーのサブセットを同定し、パラグラフ24に記載の方法にしたがって細胞におけるマーカー発現をスコア付けする工程;及び既知の表現型を有する被験体からのサンプルの参照データベースを使用して、複数の異なる検出閾値でサンプル分類の正確さを決定する工程を含む、方法。
パラグラフ42:自動化された及び/又は再帰的な方法で、手動で又はコンピュータによる方法を使用して、ローカルに、リモートサーバに、又はクラウドに保存されたデータを使用して、サンプル分類の正確さを決定する工程を含む、パラグラフ41に記載の方法。
パラグラフ43:前記検出閾値の最高性能の大きさを決定し、前記大きさを使用して、サンプル表現型を分類する際の前記確立された検出閾値の信頼性を評価する工程をさらに含む、パラグラフ41に記載の方法。
パラグラフ44:自動化された及び/又は再帰的な方法で、手動で又はコンピュータによる方法を使用して、ローカルに、リモートサーバに、又はクラウドに保存されたデータを使用して、前記検出閾値の最高性能の大きさを決定し、前記大きさを使用して、サンプル表現型を分類する際の前記確立された検出閾値の信頼性を評価する工程をさらに含む、パラグラフ41に記載の方法。
パラグラフ45:前記検出閾値の最高性能の大きさを使用して分類されていないサンプルをスコア付けし、サンプル表現型を前記サンプルに割り当てる工程をさらに含む、パラグラフ41に記載の方法。
パラグラフ46:手動で又はコンピュータによる方法を使用して、ローカルに、リモートサーバに、又はクラウドに保存されたデータを使用して、前記検出閾値の最高性能の大きさを使用して分類されていないサンプルをスコア付けし、サンプル表現型を前記サンプルに割り当てる工程をさらに含む、パラグラフ41に記載の方法。
パラグラフ47:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS1、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列から本質的に成る、パラグラフ41に記載の方法。
パラグラフ48:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS1、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列の90%から本質的に成る、パラグラフ41に記載の方法。
パラグラフ49:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS1、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列の80%から本質的に成る、パラグラフ41に記載の方法。
パラグラフ50:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS1、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列の70%から本質的に成る、パラグラフ41に記載の方法。
パラグラフ51:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS1、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列の60%から本質的に成る、パラグラフ41に記載の方法。
パラグラフ52:マーカーの前記サブセットが、表1A、表1、表2、表3、表S1、表S3、表S4、表S5、表S6、データセットS2、データセットS3において同定される遺伝子、遺伝子座、及び配列の50%から本質的に成る、パラグラフ41に記載の方法。
パラグラフ53:がんを処置する方法であって、がんを有すると診断された被験体におけるSCARの経路活性化の分子シグナルを検出する工程;SCARの経路活性化の前記分子シグナルの検出に基づいて、活性化されたSCARの遺伝子座及び/又は下流のSCAR制御遺伝子座を標的としたユーザ固有の治療的処置を生成する工程を含む、方法。
パラグラフ54:活性化されたSCAR経路の規定されたゲノム要素をサイレンシングするために、ユーザ固有の治療的処置が、ゲノム編集に基づき、このゲノム編集が、CRISPR/Cas9複合体媒介ゲノム編集を含むがこれに限定されない、パラグラフ53に記載の方法。
パラグラフ55:活性化されたSCAR経路の規定されたゲノム要素を活性化するために、ユーザ固有の治療的処置が、ゲノム編集に基づき、このゲノム編集が、CRISPR/Cas9複合体媒介ゲノム編集を含むがこれに限定されない、パラグラフ53に記載の方法。
パラグラフ56:ユーザ固有の治療的処置が、高活性抗レトロウイルス療法(HAART)の適用に基づく、パラグラフ53に記載の方法。
パラグラフ57:ユーザ固有の治療的処置が、抗レトロウイルス薬、Raltegravir(RAL、Isentress、旧MK−0518)の投与に基づく、パラグラフ53に記載の方法。
パラグラフ58:ユーザ固有の治療的処置が、転写活性のあるSCARの遺伝子座及び/又は活性化されたSCAR経路の規定されたゲノム要素に対するアンチセンス療法の適用に基づく、パラグラフ53に記載の方法。
パラグラフ59:ユーザ固有の治療的処置が、アンタゴニスト抗体又はその断片、アゴニスト抗体又はその断片、自家細胞、同種異系細胞、ペプチド、小分子、シグナル伝達タンパク質又はその断片、又は上記のうち2つ以上を含む組成物及び単一分子又は細胞療法において上記のうち2つ以上のすべて又は一部を含む組成物を含むがこれらに限定されず、活性化されたSCAR配列によってコードされるタンパク質及び/又はペプチドに対する、標的化された免疫療法の適用に基づく、パラグラフ53に記載の方法。
パラグラフ60:単独の機能として又は免疫システムの抗がんモジュレーターの活性を増大させるために、パラグラフ39−45に記載の方法が使用されて、処置された被験体の腫瘍において腫瘍浸潤リンパ球を増強する、がんを処置する方法。

Claims (2)

  1. 被験体において、がんの予後不良の指標を提供するための方法であって、
    (i)図TS4A1、図TS4A2、及び図TS4A3に規定される74遺伝子のSCARの経路シグネチャーの遺伝子;および/または(ii)図TS4B1及び図TS4B2に規定される55遺伝子のSCARの経路シグネチャーの遺伝子の発現レベルを検出し、そして(i)および/または(ii)の遺伝子セットのそれぞれの遺伝子の発現を非悪性の体細胞組織におけるそれぞれの遺伝子の発現量である参照遺伝子の発現量と比較し、がんと非悪性の体細胞組織における当該遺伝子の発現の相関係数を決定することを含む方法によって当該がんにおけるSCAR経路の活性化を評価し、
    ここで、正の相関係数はSCAR経路が活性化していないことを示し、負の相関係数はSCAR経路が活性化していることとがんの予後が不良であることを示す指標である、前記方法。
  2. 前記がんは前立腺がんである、請求項1に記載の方法。
JP2019512963A 2016-05-19 2017-05-19 臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法 Active JP6963009B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201662339007P 2016-05-19 2016-05-19
US62/339,007 2016-05-19
PCT/US2017/033678 WO2017201497A1 (en) 2016-05-19 2017-05-19 Methods of diagnosis and therapeutic targeting of clinically intractable malignant tumors

Publications (3)

Publication Number Publication Date
JP2019520851A JP2019520851A (ja) 2019-07-25
JP2019520851A5 JP2019520851A5 (ja) 2020-07-16
JP6963009B2 true JP6963009B2 (ja) 2021-11-05

Family

ID=59363200

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019512963A Active JP6963009B2 (ja) 2016-05-19 2017-05-19 臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法

Country Status (3)

Country Link
US (2) US20180057890A1 (ja)
JP (1) JP6963009B2 (ja)
WO (1) WO2017201497A1 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102933367B1 (ko) * 2018-02-27 2026-03-03 코넬 유니버시티 게놈-와이드 통합을 통한 순환 종양 dna의 초민감 검출
CN108220446B (zh) * 2018-03-29 2020-06-30 青岛泱深生物医药有限公司 Linc01356作为分子标志物在胃癌中的应用
CN110527682B (zh) * 2018-05-25 2023-09-01 中国科学院深圳先进技术研究院 一种长链非编码rna及其应用
KR20220030945A (ko) 2019-05-23 2022-03-11 크리스티아나 케어 헬스 서비시즈 인코포레이티드 암의 치료를 위한 nrf2의 유전자 녹아웃
JP7585237B2 (ja) 2019-05-23 2024-11-18 クリスティアナ ケア ジーン エディティング インスティテュート,インコーポレーテッド がんを処置するためのバリアントnrf2の遺伝子ノックアウト
CN110556158B (zh) * 2019-08-30 2022-02-15 山西农业大学 抗心肌纤维化药物的筛选方法
CN110874502B (zh) * 2019-11-11 2020-09-01 中国人民解放军国防科技大学 基于多阶段试验数据折合的航天产品可靠性评估方法
JP2023512449A (ja) * 2020-01-13 2023-03-27 エイチ リー モフィット キャンサー センター アンド リサーチ インスティテュート インコーポレイテッド Tap63制御された発がん性長鎖非コードrna
CN111893188B (zh) * 2020-08-21 2023-01-31 河北医科大学第二医院 生物标志物linc01356在宫颈癌诊断和治疗中的应用
WO2023152664A1 (en) * 2022-02-09 2023-08-17 B.Y. Quantitative Medicine Limited Analytic platform using npm1-associated genes interaction network for identifying genetic traits
WO2024163478A1 (en) * 2023-01-31 2024-08-08 Guardant Health, Inc. Non-invasive monitoring of genomic alterations induced by gene-editing therapies
CN116240243A (zh) * 2023-03-16 2023-06-09 中国人民解放军军事科学院军事医学研究院 基于dCas9-VP64的转录激活系统构建的A549稳定表达细胞系及其应用
CN121506510A (zh) * 2026-01-14 2026-02-10 奥明星程(杭州)生物科技有限公司 一种乳腺癌淋巴结转移预测的方法及其应用

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004535765A (ja) * 2000-12-07 2004-12-02 カイロン コーポレイション 前立腺癌においてアップレギュレートされた内因性レトロウイルス
EP1365034A3 (en) * 2002-05-21 2004-02-18 Bayer HealthCare AG Methods and compositions for the prediction, diagnosis, prognosis, prevention and treatment of malignant neoplasia
AU2006234897A1 (en) 2005-03-16 2006-10-19 Sidney Kimmel Cancer Center Methods and compositions for predicting death from cancer and prostate cancer survival using gene expression signatures
CA2648021A1 (en) * 2006-03-31 2007-10-11 Ordway Research Institute, Inc. Prognostic and diagnostic method for cancer therapy
FR2984360B1 (fr) * 2011-12-20 2017-11-24 Biomerieux Sa Procede pour le diagnostic ou le pronostic in vitro du cancer du colon
KR101874390B1 (ko) * 2013-09-26 2018-07-04 파이브3 제노믹스, 엘엘씨 바이러스-연관 종양을 위한 시스템, 방법, 및 조성물
WO2016011558A1 (en) * 2014-07-23 2016-01-28 Ontario Institute For Cancer Research Systems, devices and methods for constructing and using a biomarker
US9373070B2 (en) * 2014-08-04 2016-06-21 Avery Dennison Corporation Use of RFID chip as assembly facilitator

Also Published As

Publication number Publication date
WO2017201497A1 (en) 2017-11-23
US20180057890A1 (en) 2018-03-01
JP2019520851A (ja) 2019-07-25
US20230056481A1 (en) 2023-02-23

Similar Documents

Publication Publication Date Title
JP6963009B2 (ja) 臨床的に難治性の悪性腫瘍の診断及び治療標的化の方法
US12454725B2 (en) Method of detecting cancer through generalized loss of stability of epigenetic domains and compositions thereof
CN110603329B (zh) 用于诊断肝细胞癌和肺癌的甲基化标志物
KR102028375B1 (ko) 희귀 돌연변이 및 카피수 변이를 검출하기 위한 시스템 및 방법
Wang et al. Whole-exome sequencing of human pancreatic cancers and characterization of genomic instability caused by MLH1 haploinsufficiency and complete deficiency
Spencer et al. Comparison of clinical targeted next-generation sequence data from formalin-fixed and fresh-frozen tissue specimens
Birch et al. Chromosome 3 anomalies investigated by genome wide SNP analysis of benign, low malignant potential and low grade ovarian serous tumours
US20210292844A1 (en) Prostate cancer aggressiveness biomarkers
EP3704267A1 (en) Cancer detection, classification, prognostication, therapy prediction and therapy monitoring using methylome analysis
CN103797120B (zh) 前列腺癌的生物学标志物、治疗靶点及其用途
AU2019351522B2 (en) Second generation sequencing-based method for detecting microsatellite stability and genome changes by means of plasma
US20190309352A1 (en) Multimodal assay for detecting nucleic acid aberrations
CN108779487A (zh) 用于检测甲基化状态的核酸和方法
Xie et al. Genomic and transcriptomic landscape of human gastrointestinal stromal tumors
US20190249257A1 (en) Long intergenic non-coding rna as pancancer biomarker
CN104846070B (zh) 前列腺癌的生物学标志物、治疗靶点及其用途
US20250297320A1 (en) Methylation signatures in cell-free dna for tumor classification and early detection
Zawistowski et al. Unifying genomics and transcriptomics in single cells with ResolveOME amplification chemistry to illuminate oncogenic and drug resistance mechanisms
WO2014190927A1 (zh) 胰腺神经内分泌肿瘤易感基因位点及检测方法和试剂盒
van Dessel et al. The genomic landscape of metastatic castration-resistant prostate cancers using whole genome sequencing reveals multiple distinct genotypes with potential clinical impact
KR20130065058A (ko) 자발적 dna 이중 나선 손상 부위의 확인 방법 및 이의 이용
Sottoriva et al. The co-evolution of the genome and epigenome in colorectal cancer
Einstein Prognostic Biomarkers in Active Surveillance: Parsing Risk in Early-Stage Prostate Cancer
CN108753959A (zh) 一种位于disc1fp1基因的与放疗引起的放射性脑损伤相关的snp标志物及其应用
Dalley et al. Evaluation of Endobronchial Ultrasound-Guided Transbronchial Needle Aspiration (EBUS-TBNA) Samples from Advanced Non-Small Cell Lung Cancer for Whole Genome, Whole Exome and Comprehensive Panel Sequencing

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200513

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200602

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20201218

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210105

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20210127

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20210127

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20210127

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20210331

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210624

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210921

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211014

R150 Certificate of patent or registration of utility model

Ref document number: 6963009

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250