JPWO2020145351A1 - 薬物動態関連遺伝子の網羅的配列解析法とそれに使用されるプライマーセット - Google Patents
薬物動態関連遺伝子の網羅的配列解析法とそれに使用されるプライマーセット Download PDFInfo
- Publication number
- JPWO2020145351A1 JPWO2020145351A1 JP2020565207A JP2020565207A JPWO2020145351A1 JP WO2020145351 A1 JPWO2020145351 A1 JP WO2020145351A1 JP 2020565207 A JP2020565207 A JP 2020565207A JP 2020565207 A JP2020565207 A JP 2020565207A JP WO2020145351 A1 JPWO2020145351 A1 JP WO2020145351A1
- Authority
- JP
- Japan
- Prior art keywords
- pharmacokinetic
- types
- related genes
- genes
- sequences
- 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.)
- Pending
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012300 Sequence Analysis Methods 0.000 title claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 45
- 230000003321 amplification Effects 0.000 claims abstract description 16
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 16
- 230000004544 DNA amplification Effects 0.000 claims abstract description 5
- 239000003153 chemical reaction reagent Substances 0.000 claims description 3
- 108700026220 vif Genes Proteins 0.000 claims description 2
- 230000035772 mutation Effects 0.000 description 18
- 239000003814 drug Substances 0.000 description 13
- 229940079593 drug Drugs 0.000 description 13
- 108010074918 Cytochrome P-450 CYP1A1 Proteins 0.000 description 9
- 108010074922 Cytochrome P-450 CYP1A2 Proteins 0.000 description 9
- 102100031476 Cytochrome P450 1A1 Human genes 0.000 description 9
- 102100026533 Cytochrome P450 1A2 Human genes 0.000 description 9
- 230000000694 effects Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 7
- 238000007482 whole exome sequencing Methods 0.000 description 7
- 238000003505 heat denaturation Methods 0.000 description 6
- 238000000137 annealing Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- GUAHPAJOXVYFON-ZETCQYMHSA-N (8S)-8-amino-7-oxononanoic acid zwitterion Chemical compound C[C@H](N)C(=O)CCCCCC(O)=O GUAHPAJOXVYFON-ZETCQYMHSA-N 0.000 description 3
- 102000004190 Enzymes Human genes 0.000 description 3
- 108090000790 Enzymes Proteins 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000036267 drug metabolism Effects 0.000 description 3
- 238000007403 mPCR Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012163 sequencing technique Methods 0.000 description 3
- 238000012070 whole genome sequencing analysis Methods 0.000 description 3
- 108010078791 Carrier Proteins Proteins 0.000 description 2
- 238000011529 RT qPCR Methods 0.000 description 2
- 239000000090 biomarker Substances 0.000 description 2
- 230000001351 cycling effect Effects 0.000 description 2
- 210000004185 liver Anatomy 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 102000054765 polymorphisms of proteins Human genes 0.000 description 2
- 238000000746 purification Methods 0.000 description 2
- 238000003753 real-time PCR Methods 0.000 description 2
- 238000007480 sanger sequencing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 108010020070 Cytochrome P-450 CYP2B6 Proteins 0.000 description 1
- 108010026925 Cytochrome P-450 CYP2C19 Proteins 0.000 description 1
- 108010000543 Cytochrome P-450 CYP2C9 Proteins 0.000 description 1
- 108010001237 Cytochrome P-450 CYP2D6 Proteins 0.000 description 1
- 108010081668 Cytochrome P-450 CYP3A Proteins 0.000 description 1
- 102000002004 Cytochrome P-450 Enzyme System Human genes 0.000 description 1
- 108010015742 Cytochrome P-450 Enzyme System Proteins 0.000 description 1
- 102100038739 Cytochrome P450 2B6 Human genes 0.000 description 1
- 102100029363 Cytochrome P450 2C19 Human genes 0.000 description 1
- 102100029358 Cytochrome P450 2C9 Human genes 0.000 description 1
- 102100021704 Cytochrome P450 2D6 Human genes 0.000 description 1
- 102100039208 Cytochrome P450 3A5 Human genes 0.000 description 1
- 102100024902 Cytochrome P450 4F2 Human genes 0.000 description 1
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 1
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 1
- 102100022334 Dihydropyrimidine dehydrogenase [NADP(+)] Human genes 0.000 description 1
- 208000030453 Drug-Related Side Effects and Adverse reaction Diseases 0.000 description 1
- 101000909122 Homo sapiens Cytochrome P450 4F2 Proteins 0.000 description 1
- 101000902632 Homo sapiens Dihydropyrimidine dehydrogenase [NADP(+)] Proteins 0.000 description 1
- 101000991945 Homo sapiens Nucleotide triphosphate diphosphatase NUDT15 Proteins 0.000 description 1
- 101000744394 Homo sapiens Oxidized purine nucleoside triphosphate hydrolase Proteins 0.000 description 1
- 101000799388 Homo sapiens Thiopurine S-methyltransferase Proteins 0.000 description 1
- 101000621945 Homo sapiens Vitamin K epoxide reductase complex subunit 1 Proteins 0.000 description 1
- 208000029523 Interstitial Lung disease Diseases 0.000 description 1
- 206010067125 Liver injury Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 102100030661 Nucleotide triphosphate diphosphatase NUDT15 Human genes 0.000 description 1
- 102100039792 Oxidized purine nucleoside triphosphate hydrolase Human genes 0.000 description 1
- 108091008109 Pseudogenes Proteins 0.000 description 1
- 102000057361 Pseudogenes Human genes 0.000 description 1
- 108091006731 SLCO1B1 Proteins 0.000 description 1
- 102100027233 Solute carrier organic anion transporter family member 1B1 Human genes 0.000 description 1
- 102100034162 Thiopurine S-methyltransferase Human genes 0.000 description 1
- 102100029152 UDP-glucuronosyltransferase 1A1 Human genes 0.000 description 1
- 101710205316 UDP-glucuronosyltransferase 1A1 Proteins 0.000 description 1
- 102100023485 Vitamin K epoxide reductase complex subunit 1 Human genes 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000000941 bile Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 101150055214 cyp1a1 gene Proteins 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003623 enhancer Substances 0.000 description 1
- 230000029142 excretion Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 210000001035 gastrointestinal tract Anatomy 0.000 description 1
- 108091008053 gene clusters Proteins 0.000 description 1
- 238000003205 genotyping method Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 231100000234 hepatic damage Toxicity 0.000 description 1
- 230000008818 liver damage Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 210000002700 urine Anatomy 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
- C12N15/09—Recombinant DNA-technology
- C12N15/11—DNA or RNA fragments; Modified forms thereof; Non-coding nucleic acids having a biological activity
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/686—Polymerase chain reaction [PCR]
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Genetics & Genomics (AREA)
- Biotechnology (AREA)
- General Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biochemistry (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Microbiology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Biomedical Technology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Plant Pathology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法であって、前記複数種類の薬物動態関連遺伝子の配列を相互に区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含み、前記プライマーセットはそれぞれ260〜490bpの領域を増幅するように設計されていることを特徴とする方法。
Description
本発明は遺伝子解析に関し、より詳細には遺伝子解析を利用した研究、診断、医療等に関する。
消化管からの吸収、肝臓での薬物代謝、胆汁・尿への排泄などの薬物体内動態における個人差により、薬の血中濃度に個人間変動が生じる結果、薬の効果や副作用の発現リスクに個人差が見られることが知られている。このような薬物動態の個人差は薬物代謝及び薬物の生体膜輸送に関連する遺伝子群の塩基配列の違いによることが多く、このような遺伝子配列の個人差を調べることは、薬の効果や副作用を予測するためのゲノムバイオマーカーの同定につながる。
しかしながら、薬物動態関連遺伝子は複数のメンバーを構成因子とするファミリーとして存在することが多く、互いに非常に類似した配列を有するため、PCRでは各構成因子を一緒に増幅してしまうことがあり、また、配列解析では、配列の違いが実際の配列の違いによるものか、配列解析における誤読（疑陽性）によるものかが区別できないこともあった。
このように、薬物動態関連遺伝子の多くは複雑なゲノム領域に存在するため、従来のPCRや次世代シークエンサーを用いる方法では高精度かつ網羅的な遺伝子型判定は困難であった。
このように、薬物動態関連遺伝子の多くは複雑なゲノム領域に存在するため、従来のPCRや次世代シークエンサーを用いる方法では高精度かつ網羅的な遺伝子型判定は困難であった。
Gordon et al., Hum Mol Genet, 2014, Apr. 15, Vol. 23（8）:1957-63
Twist et al., NPJ Genom Med, 2016, Jan. 13, 1:15007
上記のように薬物動態関連遺伝子を網羅的に配列解析する方法としては、次世代シークエンサーを用いた全ゲノムシークエンシング（WGS）や全エクソームシークエンシング（WES）などがある。しかしながら、WGSやWESは時間とコストがかかるうえ、薬物動態関連遺伝子が複雑なゲノム領域に存在し、高い配列相同性を有することから、薬物動態関連遺伝子の変異を精度よく解析することは困難であった。そこで、本発明は薬物動態関連遺伝子の配列を効率よく、かつ精度よく網羅的に解析するための方法を提供することを課題とする。
本発明者らは上記課題を解決するために鋭意検討を行った。その結果、薬物動態関連遺伝子の配列を相互に区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットであって、それぞれ260〜490bpの領域を増幅するように設計されたプライマーセットを用い、複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および第一工程で得られた遺伝子増幅産物を配列解析する第二工程、を行うことにより、被験者由来の複数種類の薬物動態関連遺伝子の配列を効率よく、かつ高精度で網羅的に解析することができることを見出した。
より具体的には、薬物代謝に関与することが報告されている60種類の薬物代謝酵素、37種類の薬物トランスポーターを含む薬物動態関連100遺伝子をリストアップした。それらの遺伝子の翻訳領域をターゲットとし、特定のゲノム領域のみを増幅可能なマルチプレックスPCR法でターゲット領域を濃縮することにより、100遺伝子の機能に影響する遺伝子多型及び極めて稀な変異 (レアバリアント) を網羅的かつ高精度に検出することができるターゲット・リシークエンス解析パネルを開発することに成功した。以上に基づき、本発明を完成させた。
より具体的には、薬物代謝に関与することが報告されている60種類の薬物代謝酵素、37種類の薬物トランスポーターを含む薬物動態関連100遺伝子をリストアップした。それらの遺伝子の翻訳領域をターゲットとし、特定のゲノム領域のみを増幅可能なマルチプレックスPCR法でターゲット領域を濃縮することにより、100遺伝子の機能に影響する遺伝子多型及び極めて稀な変異 (レアバリアント) を網羅的かつ高精度に検出することができるターゲット・リシークエンス解析パネルを開発することに成功した。以上に基づき、本発明を完成させた。
本発明の要旨は以下のとおりである。
[１]被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法であって、
前記複数種類の薬物動態関連遺伝子の配列を相互に区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および
第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含み、
前記プライマーセットはそれぞれ260〜490bpの領域を増幅するように設計されていることを特徴とする方法。
[２]前記複数種類の薬物動態関連遺伝子それぞれについて、複数の領域の増幅反応を行う、[１]に記載の方法。
[３]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される25種類以上の遺伝子である、[１]または[２]に記載の方法。
[４]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される50種類以上の遺伝子である、[１]または[２]に記載の方法。
[５]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類全ての遺伝子である、[１]または[２]に記載の方法。
[６]前記複数の反応系は、６〜１０の反応系である、[１]〜[５]のいずれかに記載の方法。
[７]表１に記載されたプライマーセットを使用する、[１]〜[６]のいずれかに記載の方法。
[８]被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法に使用される試薬であって、表２に記載のプライマーセットから選択される複数のプライマーセットを含む試薬。
[１]被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法であって、
前記複数種類の薬物動態関連遺伝子の配列を相互に区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および
第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含み、
前記プライマーセットはそれぞれ260〜490bpの領域を増幅するように設計されていることを特徴とする方法。
[２]前記複数種類の薬物動態関連遺伝子それぞれについて、複数の領域の増幅反応を行う、[１]に記載の方法。
[３]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される25種類以上の遺伝子である、[１]または[２]に記載の方法。
[４]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される50種類以上の遺伝子である、[１]または[２]に記載の方法。
[５]解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類全ての遺伝子である、[１]または[２]に記載の方法。
[６]前記複数の反応系は、６〜１０の反応系である、[１]〜[５]のいずれかに記載の方法。
[７]表１に記載されたプライマーセットを使用する、[１]〜[６]のいずれかに記載の方法。
[８]被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法に使用される試薬であって、表２に記載のプライマーセットから選択される複数のプライマーセットを含む試薬。
本発明によれば、薬物動態関連遺伝子の配列を効率よく、低コストかつ精度よく網羅的に解析することができる。例えば、本発明の解析手法を用いて、がん分子標的薬などによる肝障害や間質性肺炎などの重篤な副作用が起こった患者由来のDNAサンプルを解析することにより、医学的に有用なゲノムバイオマーカーとなる機能的多型を同定することができる。
さらに、本発明によれば、GWASでは検出できない稀な変異 (rare variant) を解析可能である。また、重篤な副作用に関する変異の有無を調べることにより、薬剤の投与前に副作用予測ができ、治療方針の決定に有用である。また、本発明によれば、多くの対象のサンプルを高スループットで解析できる。
本発明の方法は、被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法であって、
前記複数種類の薬物動態関連遺伝子の配列を区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および
第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含む。
前記複数種類の薬物動態関連遺伝子の配列を区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および
第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含む。
＜薬物動態関連遺伝子＞
薬物動態関連遺伝子としては、薬物代謝酵素または薬物トランスポーターをコードする遺伝子が含まれるが、その変異が薬物の副作用に関与する可能性が有る限りにおいて、これらに分類されない遺伝子でもよい。例えば、表１に記載のNUDT1、NUDT15、VKORC1なども薬物動態関連遺伝子に含まれる。薬物動態関連遺伝子は肝臓で発現する遺伝子であることが好ましい。
薬物動態関連遺伝子としては、薬物代謝酵素または薬物トランスポーターをコードする遺伝子が含まれるが、その変異が薬物の副作用に関与する可能性が有る限りにおいて、これらに分類されない遺伝子でもよい。例えば、表１に記載のNUDT1、NUDT15、VKORC1なども薬物動態関連遺伝子に含まれる。薬物動態関連遺伝子は肝臓で発現する遺伝子であることが好ましい。
薬物動態関連遺伝子としてより具体的には、表１に記載されたものが挙げられる。
薬物動態関連遺伝子は２５種類以上解析することが好ましく、５０種類以上解析することが好ましく、１００種類全て解析することがさらに好ましい。
薬物動態関連遺伝子の中では、その変異が薬剤の副作用に関与することが知られているCYP2D6、CYP2C9、CYP2C19、CYP4F2、CYP2B6、DPYD、TPMT、SLCO1B1、UGT1A1、CYP3A5などがあり、少なくともこれらを解析することが好ましい。
薬物動態関連遺伝子は２５種類以上解析することが好ましく、５０種類以上解析することが好ましく、１００種類全て解析することがさらに好ましい。
薬物動態関連遺伝子の中では、その変異が薬剤の副作用に関与することが知られているCYP2D6、CYP2C9、CYP2C19、CYP4F2、CYP2B6、DPYD、TPMT、SLCO1B1、UGT1A1、CYP3A5などがあり、少なくともこれらを解析することが好ましい。
以下、各工程について説明する。
＜第一工程＞
第一工程では薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う。
第一工程では薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う。
＜プライマーセット＞
プライマーセットは複数種類の薬物動態関連遺伝子の配列を区別して増幅するために薬物動態関連遺伝子ごとに設計される。
例えば、A,B,C,D,E,F,G,H,I,Jの１０種類の遺伝子を解析する場合、これらA〜Jの１０種類の遺伝子に対し、それぞれを特異的に増幅するためのプライマーセット１０組を用意する。各遺伝子に対するプライマーセットは、公知の遺伝子配列に基づき、他の遺伝子とは相同性の少ない領域に設定することができる。
プライマーセットは複数種類の薬物動態関連遺伝子の配列を区別して増幅するために薬物動態関連遺伝子ごとに設計される。
例えば、A,B,C,D,E,F,G,H,I,Jの１０種類の遺伝子を解析する場合、これらA〜Jの１０種類の遺伝子に対し、それぞれを特異的に増幅するためのプライマーセット１０組を用意する。各遺伝子に対するプライマーセットは、公知の遺伝子配列に基づき、他の遺伝子とは相同性の少ない領域に設定することができる。
特に、同じファミリーに属する遺伝子同士を区別して増幅するためには、両遺伝子の配列を比較し、配列が異なる２ヵ所にプライマーを設計することが好ましい。
また、チトクロームP450等の薬物動態関連遺伝子は偽遺伝子や類似性の高い配列が多い。そのため汎用されている次世代シークエンサー (NGS) の２ｘ１００bp程度のリード長では各遺伝子の配列を区別することができなくなるため、前記プライマーセットはそれぞれ260〜490bpの領域を増幅するように設計される。この範囲の長さの領域を増幅し、解析することで、増幅効率を低下させることなく配列を増幅でき、さらに、相同性の高い配列同士を正確に区別して精度の高い配列解析ができる。
なお、１つの薬物動態関連遺伝子につき、複数の領域を増幅して配列解析してもよい。これにより、より多くの変異を検出でき、遺伝子同士の区別もより正確に行うことができる。
増幅する領域はコード領域（エクソン）であることが好ましく、複数の領域を増幅する場合、各エクソンについてそれぞれ増幅することが好ましい。なお、1つのエクソンについて複数領域を増幅してもよい。
増幅する領域はコード領域（エクソン）であることが好ましく、複数の領域を増幅する場合、各エクソンについてそれぞれ増幅することが好ましい。なお、1つのエクソンについて複数領域を増幅してもよい。
図２に、CYP1A1とCYP1A2を例として、プライマーセットの設計の例を示す。
CYP1A1とCYP1A2は同じファミリー（CYP1Aファミリー）に属する遺伝子であり、配列が互いに類似している。
これらの配列を比較し、エクソン１に３ヵ所、エクソン２に１ヵ所、エクソン３に１ヵ所、エクソン４に１ヵ所、エクソン５に１ヵ所、エクソン６に２ヵ所、合計９カ所、両者で配列が異なる領域を増幅する。
CYP1A1のプライマーF1（forward primer 1）、R1（reverse primer 1）、およびCYP1A2のプライマーF1、R1はそれぞれCYP1A1、CYP1A2の類似領域（約300bp）を増幅するように設計されているが、前者はCYP1A1のみ、後者はCYP1A2のみにハイブリダイズするように設計されている。
CYP1A1とCYP1A2は同じファミリー（CYP1Aファミリー）に属する遺伝子であり、配列が互いに類似している。
これらの配列を比較し、エクソン１に３ヵ所、エクソン２に１ヵ所、エクソン３に１ヵ所、エクソン４に１ヵ所、エクソン５に１ヵ所、エクソン６に２ヵ所、合計９カ所、両者で配列が異なる領域を増幅する。
CYP1A1のプライマーF1（forward primer 1）、R1（reverse primer 1）、およびCYP1A2のプライマーF1、R1はそれぞれCYP1A1、CYP1A2の類似領域（約300bp）を増幅するように設計されているが、前者はCYP1A1のみ、後者はCYP1A2のみにハイブリダイズするように設計されている。
各遺伝子に対するプライマーセットの例を表２に示す。解析対象遺伝子の種類に応じ、これらプライマーセットから適宜選択して使用することができる。例えば、表１の100遺伝子のうち、25種類の遺伝子を解析する場合、選択された25種類の遺伝子につき、表２に記載のプライマーセットを使用すればよい。
＜増幅反応＞
前記複数種類の薬物動態関連遺伝子それぞれについて、増幅反応を行う
増幅反応は複数の反応系に分けて行う。ここで複数の反応系とは、好ましくは、４以上、１５以下、より好ましくは６以上、１０以下、特に好ましくは８である。
前記複数種類の薬物動態関連遺伝子それぞれについて、増幅反応を行う
増幅反応は複数の反応系に分けて行う。ここで複数の反応系とは、好ましくは、４以上、１５以下、より好ましくは６以上、１０以下、特に好ましくは８である。
例えば、１００種類の遺伝子につき、それぞれ８ヵ所ずつ増幅する場合、合計８００の反応を行うことになるが、これらを複数の反応系に分ける。例えば、１反応系当たり、１００プライマーセットを用いた１００の反応を行い、合計８つの反応系で反応を行うことができる。一反応系において、例えば、５０〜２００のプライマーセットを用いた５０〜２００の増幅反応を同時に行うことができる。
なお、各反応系の反応数は必ずしも同じである必要はない。
なお、各反応系の反応数は必ずしも同じである必要はない。
このように複数の反応系を分ける場合、特に同一ファミリーに存在する遺伝子同士を区別して増幅する場合、これら遺伝子同士の対応領域を増幅するための反応は別々の反応系で行うようにすることが好ましい。
例えば、図２に示すように、CYP1A1とCYP1A2が解析対象に含まれる場合は、これらを確実に分けて増幅するために、CYP1A1とCYP1A2で類似する領域を増幅するプライマーセット同士、ここではCYP1A1のF1-R1とCYP1A2のF1-R1は別々の反応系に加え、増幅反応を別々に行う。これにより、類似遺伝子同士を分けて増幅することができるので、配列解析における疑陽性を確実に避けることができる。
例えば、図２に示すように、CYP1A1とCYP1A2が解析対象に含まれる場合は、これらを確実に分けて増幅するために、CYP1A1とCYP1A2で類似する領域を増幅するプライマーセット同士、ここではCYP1A1のF1-R1とCYP1A2のF1-R1は別々の反応系に加え、増幅反応を別々に行う。これにより、類似遺伝子同士を分けて増幅することができるので、配列解析における疑陽性を確実に避けることができる。
反応は通常のPCRの条件に準じて行うことができる。
例えば、熱変性、アニーリング、増幅のサイクルを２０〜３０サイクル行う。これらの各ステップの温度はプライマーのTm値や使用する酵素などに応じて適宜設定できる。
例えば、熱変性、アニーリング、増幅のサイクルを２０〜３０サイクル行う。これらの各ステップの温度はプライマーのTm値や使用する酵素などに応じて適宜設定できる。
＜第二工程＞
第二工程では、第一工程で得られた各薬物動態関連遺伝子の増幅産物について配列解析を行う。
なお、配列解析の前に、配列解析用の配列を付加するために、配列解析用の配列を含むプライマーを用いて増幅してもよい
配列解析の方法は特に制限されないが、次世代シークエンサーなどを用いることが好ましい。
次世代シークエンサーなどを用いて網羅的に解析し、各遺伝子を野生型配列と比較することで変異の有無などを調べることができる。
第二工程では、第一工程で得られた各薬物動態関連遺伝子の増幅産物について配列解析を行う。
なお、配列解析の前に、配列解析用の配列を付加するために、配列解析用の配列を含むプライマーを用いて増幅してもよい
配列解析の方法は特に制限されないが、次世代シークエンサーなどを用いることが好ましい。
次世代シークエンサーなどを用いて網羅的に解析し、各遺伝子を野生型配列と比較することで変異の有無などを調べることができる。
本発明の配列解析法により、薬物動態関連遺伝子に変異が見られた場合、その変異のタイプに応じ、薬物の投与を調節することができる。例えば、特定の薬剤に対する副作用に対して高リスク型変異を有する場合、その薬剤の投与を中止したり、投与量を減らしたりすることができる。
以下、実施例を挙げて本発明を具体的に説明する。ただし、本発明は以下の実施態様には限定されない。
本発明の方法を用いたターゲット・リシークエンシングのワークフローを図１に示す。
表１に記載の薬物動態関連100遺伝子における翻訳（コード）領域を特異的に増幅することができる1st PCR用のプライマー (1,102セット) を設計した（表２）。プライマーセットを8つのグループに分け、それぞれのグループでマルチプレックスPCRを実施する。得られたPCR産物に、イルミナ社MiSeqの解析に必要なアダプター配列や多検体識別用のインデックス配列を付加するために2ndPCRを実施し、得られたライブラリーをイルミナ社のMiSeqでシークエンシングする (2 × 250 bp)。
具体的手順は以下の通り。
表１に記載の薬物動態関連100遺伝子における翻訳（コード）領域を特異的に増幅することができる1st PCR用のプライマー (1,102セット) を設計した（表２）。プライマーセットを8つのグループに分け、それぞれのグループでマルチプレックスPCRを実施する。得られたPCR産物に、イルミナ社MiSeqの解析に必要なアダプター配列や多検体識別用のインデックス配列を付加するために2ndPCRを実施し、得られたライブラリーをイルミナ社のMiSeqでシークエンシングする (2 × 250 bp)。
具体的手順は以下の通り。
（１）プライマーの設計
薬物動態関連100遺伝子の翻訳領域は合計162,216 bpであった。これらの翻訳領域を特異的に増幅するために、260-490 bpのPCR産物が得られるように、15-30 bpのプライマーを設計した。最終的に、1,102セットのプライマー (表２) を8つのグループに分けた後、1st PCRを実施し、159,347 bpの領域を増幅した。
薬物動態関連100遺伝子の翻訳領域は合計162,216 bpであった。これらの翻訳領域を特異的に増幅するために、260-490 bpのPCR産物が得られるように、15-30 bpのプライマーを設計した。最終的に、1,102セットのプライマー (表２) を8つのグループに分けた後、1st PCRを実施し、159,347 bpの領域を増幅した。
（２）PCR
1stPCRは、サーマルサイクラー(GeneAmp PCR System 9700、Thermo Fisher Scientific) により、5 ng/μLテンプレートDNA 2 μL、2×Platinum Multiplex PCR Master Mix (Thermo Fisher Scientific) 5 μL、GCエンハンサー 1 μL、0.1 nMプライマー溶液2 μL、反応溶液総量10 μLについて、初期熱変性を95℃、2分間、サイクリングは、熱変性を95℃、30秒間、アニーリングを60℃、90秒間、伸長反応を72℃、3分間で25サイクルとし、最終伸長を72℃、10分間行った。
1stPCRは、サーマルサイクラー(GeneAmp PCR System 9700、Thermo Fisher Scientific) により、5 ng/μLテンプレートDNA 2 μL、2×Platinum Multiplex PCR Master Mix (Thermo Fisher Scientific) 5 μL、GCエンハンサー 1 μL、0.1 nMプライマー溶液2 μL、反応溶液総量10 μLについて、初期熱変性を95℃、2分間、サイクリングは、熱変性を95℃、30秒間、アニーリングを60℃、90秒間、伸長反応を72℃、3分間で25サイクルとし、最終伸長を72℃、10分間行った。
2ndPCRは、サーマルサイクラー(GeneAmp PCR System 9700、Thermo Fisher Scientific) により、5×KAPA HiFi Fidelity Buffer (Kapa Biosystems) 2 μL、10 mM dNTP 0.3 μL、0.2 U KAPA HiFi HotStart DNA Polymerase (Kapa Biosystems)、プライマー溶液 (5 μM equimolar mix of index 1 and index 2) 0.05 μL、1st PCR産物 2 μL、反応溶液総量5 μLについて、初期熱変性を98℃、45秒間、サイクリングは、熱変性を99℃、15秒間、アニーリングを65℃、30秒間、伸長反応を72℃、30秒間で4サイクルとし、最終伸長を72℃、1分間行った。
qPCRをリアルタイムPCRシステム (Applied Biosystems 7900HT Fast Real-Time PCR System、Thermo Fisher Scientific) により、2×KAPA SYBR FAST qPCR Master Mix (Kapa Biosystems) 5 μL、2ndPCR産物の精製後産物 2 μL、反応溶液総量10 μLについて、初期熱変性を95℃、5分間、サイクリングは、熱変性を95℃、30秒間、アニーリング及び伸長反応を60℃、45秒間で40サイクルの条件下で行った。
（３）ライブラリー精製とシークエンスラン2ndPCR産物より、AMPure XP (Beckman Coulter) を用いて目的サイズのDNA断片を回収し、ライブラリーを精製した。各グループのライブラリーを4 nMに希釈した後、等量混合し、Illumina MiSeqシステム及びMiSeq Reagent Kit v2を用いてシークエンシングを行った(2 × 250 bp、500サイクル)。
（４）変異の検出
変異の検出を、シークエンシングによって得られたFastqファイルのヒト参照ゲノム(hg19)へのマッピングにＢＷＡ(Li H et al、Bioinformatics、２００９、２５：１７５４）を用いて、変異の検出及びフィルタリングにGATK（McKenna A et al、Genome Res、２０１０、２０（９）：１２９７−３０３）を用いて、行った。また、コピー数変異(CNV)の検出を、BED Tools（Quinlan AR、２０１４、４７：11、12、１−３４）を用いて、行った。
変異の検出を、シークエンシングによって得られたFastqファイルのヒト参照ゲノム(hg19)へのマッピングにＢＷＡ(Li H et al、Bioinformatics、２００９、２５：１７５４）を用いて、変異の検出及びフィルタリングにGATK（McKenna A et al、Genome Res、２０１０、２０（９）：１２９７−３０３）を用いて、行った。また、コピー数変異(CNV)の検出を、BED Tools（Quinlan AR、２０１４、４７：11、12、１−３４）を用いて、行った。
その結果、図３に示すように、一人のゲノムから１６３の変異を検出した。同じゲノムにWESを用いたところ１６２の変異を検出した。変異検出のゴールドスタンダードであるサンガーシークエンシングの結果と比較したところ、本発明の感度及び特異度は共に１００％を示した。一方WESの感度及び特異度はそれぞれ９１．４％及び９９．９９２％を示した。WESとサンガーシークエンシングの結果の不一致の内訳は偽陽性が１３変異及びリード数の不足による未検出が１４変異であった。
Claims (8)
- 被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法であって、
前記複数種類の薬物動態関連遺伝子の配列を相互に区別して増幅するために薬物動態関連遺伝子ごとに設計されたプライマーセットを用い、前記複数種類の薬物動態関連遺伝子それぞれについての増幅反応を複数の反応系に分けて行う第一工程、および
第一工程で得られた遺伝子増幅産物を配列解析する第二工程を含み、
前記プライマーセットはそれぞれ260〜490bpの領域を増幅するように設計されていることを特徴とする方法。 - 前記複数種類の薬物動態関連遺伝子それぞれについて、複数の領域の増幅反応を行う、請求項１に記載の方法。
- 解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される25種類以上の遺伝子である、請求項１または２に記載の方法。
- 解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類から選択される50種類以上の遺伝子である、請求項１または２に記載の方法。
- 解析対象の前記複数種類の薬物動態関連遺伝子は、表１の100種類全ての遺伝子である、請求項１または２に記載の方法。
- 前記複数の反応系は、６〜１０の反応系である、請求項１〜５のいずれか一項に記載の方法。
- 表１に記載されたプライマーセットを使用する、請求項１〜６のいずれか一項に記載の方法。
- 被験者由来の複数種類の薬物動態関連遺伝子の配列を網羅的に解析する方法に使用される試薬であって、表２に記載のプライマーセットから選択される複数のプライマーセットを含む試薬。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019001821 | 2019-01-09 | ||
| JP2019001821 | 2019-01-09 | ||
| PCT/JP2020/000485 WO2020145351A1 (ja) | 2019-01-09 | 2020-01-09 | 薬物動態関連遺伝子の網羅的配列解析法とそれに使用されるプライマーセット |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPWO2020145351A1 true JPWO2020145351A1 (ja) | 2021-11-25 |
Family
ID=71520270
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020565207A Pending JPWO2020145351A1 (ja) | 2019-01-09 | 2020-01-09 | 薬物動態関連遺伝子の網羅的配列解析法とそれに使用されるプライマーセット |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JPWO2020145351A1 (ja) |
| WO (1) | WO2020145351A1 (ja) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003144176A (ja) * | 2000-12-27 | 2003-05-20 | Inst Of Physical & Chemical Res | 遺伝子多型の検出方法 |
| JP2004000039A (ja) * | 2002-05-30 | 2004-01-08 | Inst Of Physical & Chemical Res | 遺伝子多型の検出方法 |
| JP2004000004A (ja) * | 2002-04-23 | 2004-01-08 | Kokuritsu Iyakuhin Shokuhin Eisei Kenkyusho | 薬剤代謝へ影響を及ぼすｃｙｐ３ａ４遺伝子多型、およびその利用 |
| WO2007055255A1 (ja) * | 2005-11-08 | 2007-05-18 | Olympus Corporation | 複数の識別用核酸配列を増幅する方法 |
| JP2018506307A (ja) * | 2015-02-11 | 2018-03-08 | パラゴン ゲノミクス, インク.Paragon Genomics, Inc. | 非特異的増幅産物を減少させる為の方法及び組成物 |
-
2020
- 2020-01-09 JP JP2020565207A patent/JPWO2020145351A1/ja active Pending
- 2020-01-09 WO PCT/JP2020/000485 patent/WO2020145351A1/ja active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003144176A (ja) * | 2000-12-27 | 2003-05-20 | Inst Of Physical & Chemical Res | 遺伝子多型の検出方法 |
| JP2004000004A (ja) * | 2002-04-23 | 2004-01-08 | Kokuritsu Iyakuhin Shokuhin Eisei Kenkyusho | 薬剤代謝へ影響を及ぼすｃｙｐ３ａ４遺伝子多型、およびその利用 |
| JP2004000039A (ja) * | 2002-05-30 | 2004-01-08 | Inst Of Physical & Chemical Res | 遺伝子多型の検出方法 |
| WO2007055255A1 (ja) * | 2005-11-08 | 2007-05-18 | Olympus Corporation | 複数の識別用核酸配列を増幅する方法 |
| JP2018506307A (ja) * | 2015-02-11 | 2018-03-08 | パラゴン ゲノミクス, インク.Paragon Genomics, Inc. | 非特異的増幅産物を減少させる為の方法及び組成物 |
Non-Patent Citations (2)
| Title |
|---|
| BMC BIOINFORMATICS, vol. 18:14, JPN7020001018, 2017, ISSN: 0005179918 * |
| ONCOTARGET, vol. 9, no. 51, JPN7020001019, 2018, pages 29789 - 29800, ISSN: 0005179917 * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020145351A1 (ja) | 2020-07-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210222245A1 (en) | Autism associated genetic markers | |
| Lam et al. | Noninvasive prenatal diagnosis of monogenic diseases by targeted massively parallel sequencing of maternal plasma: application to β-thalassemia | |
| Xu et al. | Long-molecule sequencing: a new approach for identification of clinically significant DNA variants in α-thalassemia and β-thalassemia carriers | |
| KR102028375B1 (ko) | 희귀 돌연변이 및 카피수 변이를 검출하기 위한 시스템 및 방법 | |
| Kote-Jarai et al. | Seven prostate cancer susceptibility loci identified by a multi-stage genome-wide association study | |
| JP2023017006A (ja) | Ｄｎａのサイズに基づく解析 | |
| Qin et al. | Ultra deep sequencing detects a low rate of mosaic mutations in tuberous sclerosis complex | |
| Milbury et al. | COLD-PCR enrichment of rare cancer mutations prior to targeted amplicon resequencing | |
| JP7198199B2 (ja) | マイクロサテライト不安定性を同定する方法 | |
| Betensky et al. | Clonal evolution and clinical significance of copy number neutral loss of heterozygosity of chromosome arm 6p in acquired aplastic anemia | |
| US20100184044A1 (en) | Detecting Genetic Abnormalities | |
| Song et al. | Selective multiplexed enrichment for the detection and quantitation of low-fraction DNA variants via low-depth sequencing | |
| Thompson et al. | Single-step capture and sequencing of natural DNA for detection of BRCA1 mutations | |
| Jiang et al. | Identification of thalassemia gene cluster deletion by long‐read whole‐genome sequencing (LR‐WGS) | |
| Kristensen et al. | Targeted ultradeep next‐generation sequencing as a method for KIT D 816 V mutation analysis in mastocytosis | |
| Watson et al. | Long‐read nanopore DNA sequencing can resolve complex intragenic duplication/deletion variants, providing information to enable preimplantation genetic diagnosis | |
| Verma et al. | A nanopore sequencing-based pharmacogenomic panel to personalize tuberculosis drug dosing | |
| Said et al. | Genomics in premature infants: a non-invasive strategy to obtain high-quality DNA | |
| WO2020145351A1 (ja) | 薬物動態関連遺伝子の網羅的配列解析法とそれに使用されるプライマーセット | |
| Jadhav et al. | A phenome-wide association study of methylated GC-rich repeats identifies a GCC repeat expansion in AFF3 as a significant cause of intellectual disability | |
| EP4414457A1 (en) | Method for analyzing structural polymorphism, primer pair set, and method for designing primer pair set | |
| JP2023041117A (ja) | 薬物の副作用関連遺伝子の配列解析方法とそれに使用されるプライマーセット | |
| US11028442B2 (en) | Simultaneous detection of multiple nucleic acid templates using modified primers | |
| WO2022239485A1 (ja) | 急性骨髄性白血病の遺伝子パネル検査のためのアンプリコンｄｎａライブラリー、キット及びそれらの使用 | |
| Calabrese et al. | Methods for the identification of mitochondrial DNA variants |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20221206 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231024 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20231222 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240216 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20240604 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240904 |