JP2024034434A - 高感度かつ定量的な遺伝子検査方法 - Google Patents

Abstract

【課題】より定性検出能が向上し、定量範囲が拡張され、かつ定量性が向上した遺伝子検査方法を提供すること。【解決手段】ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査する方法であって、前記異なる遺伝子配列を２組以上のプライマーセットを用い、前記２組以上のプライマーセットの量比を変えて同時に増幅し、ＰＣＲのサイクル毎にＨＲＭ解析を行い、当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算し、スタンダードの量的数値とＣｉとの間で回帰式を算出し、サンプルのＣｉ値を代入することでサンプルの標的遺伝子配列の量的数値を求めることを特徴とする、サンプル中の異なる遺伝子配列の存在を定量的に検査する方法。【選択図】なし

Description

本発明は、ＰＣＲを用いて遺伝子増幅する際、混在する異なる配列のＤＮＡを同時に高感度に検出するのみならず、定量的に測定するＡｌｌｅｌｅ Ｓｐｅｃｉｆｉｃ ＰＣＲに関する。

がんの早期発見・早期診断には、現状ではＰＥＴ等の画像診断が最も普及している手段の一つである。しかし１ｃｍ程度の大きさに成長しないと発見できないことや高額である等、早期発見には不十分であり、拘束時間が長く集団検診には適していない。それらの短所を解決しうる手法の一つにリキッドバイオプシーがあるが、早期発見を想定する場合、がん細胞が画像診断で検出できない程度の微小な段階で液性検体から検出する必要がある。液性検体の代表的検体である血漿や尿においては正常細胞由来のせん断されたゲノムＤＮＡが存在するが、がん患者にはがん細胞由来のせん断されたＤＮＡ（ｃｉｒｃｕｌａｔｉｎｇ ｔｕｍｏｒ ＤＮＡ：ｃｔＤＮＡ）が混在し、その量は、がんのステージに左右される。そのため極初期のがんを検出するには大量の検体からＤＮＡを抽出して検出する必要があり、そこには正常細胞由来ＤＮＡが大過剰に存在することになる。そのため早期発見や微小残存病変の検出にはｃｔＤＮＡを、０．０１％を超える超高感度で解析する必要があるとされる（非特許文献１）。一般の遺伝子検査室に普及している機器を用い、蛍光プローブ等の高額な試薬を使用せずにそのような高感度で定量的な手法としては、本発明者が報告した特許文献１記載のＨｉｇｈｌｙ Ｑｕａｎｔｉｔａｔｉｖｅ Ａｌｌｅｌｅ Ｓｐｅｃｉｆｉｃ ＰＣＲ（ＨｉＱＡＳＰ）がある。

なお、従来のリキッドバイオプシーの検体として最も一般的な血漿に全血中のすべてのｃｔＤＮＡが存在するとは限らない。白血球にｃｔＤＮＡが濃縮されているという報告（非特許文献２）や細胞外小胞がｃｔＤＮＡ局在に重要であるとの報告（非特許文献３）がある。つまり血漿分離後に廃棄していた血球画分にｃｔＤＮＡがかなり残存している可能性が高い。ｃｔＤＮＡを漏れなく解析するには血漿以外の画分をすべて含んだ全血検体からＤＮＡを抽出して解析する必要があるが、そこにはｃｔＤＮＡ解析を妨害する正常ＤＮＡが大量に混在することになり、さらに超高感度な解析がより重要となる。

国際公開第２０１９／０４５０１６号明細書

Ｐａｎａｄｏｐｏｕｌｏｓ Ｎ．，Ｔｕｍｏｒ Ｌｉｑｕｉｄ Ｂｉｏｐｓｙ，２０２０ ＰＭＩＤ： ２９９７１９１７ ＰＭＩＤ： ３４５９５８４２

特許文献１記載のＨｉＱＡＳＰは、融解曲線解析を利用した従来の高解像度融解曲線解析＝Ｈｉｇｈ Ｒｅｓｏｌｕｔｉｏｎ Ｍｅｌｔ（ＨＲＭ）解析に対して、（１）各アレルの増幅プライマーにＴｍ値を上昇させるタグを付加してピーク分離を向上させることで定量性を改善し、（２）各アレル検出用プライマーセットの比を変化させることで検出・定量範囲を調整でき、ｃｔＤＮＡなどの微量の変異遺伝子を高感度で検出・定量できる方法である。しかし、一度の解析での定量性を示す範囲に制限がある点や、従来の方法と同様に極微量域での詳細な定量性が困難であること、ゲノムＤＮＡ量を増やすと中間ピークを生じる例が増加し、定量性が困難になる現象があった。またＨＲＭ解析一般的に、混在するＰＣＲ産物の量比はサイクル毎に変動しうるが、従来のＨＲＭ解析では解析したＰＣＲサイクル以外のサイクル時点の情報を得ることが不可能であった。
従って、本発明の課題は、より定量範囲が拡張され、かつ定量性が向上した超高感度遺伝子検査方法を提供することにある。

そこで、本発明者は、前記ＨｉＱＡＳＰと一般的なＨＲＭ技術をさらに改良すべく種々検討した結果、ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査するにあたり、前記異なる遺伝子配列を２組以上のプライマーセットを用い、前記２組以上のプライマーセットの量比を変えて同時に増幅し、ＰＣＲのサイクル毎にＨＲＭ解析を行うことで当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を算出し、それがスタンダードの量的数値と良好な相関を示すことを発見し、そこで算出した近似式を回帰式としてサンプルのＣｉ値を代入することで検体の標的遺伝子配列の量的数値を求めるという新たなアルゴリズムを開発し、より定量範囲が拡張され、かつ定量性が向上し、全血から抽出されたＤＮＡを用いた場合であっても、ｃｔＤＮＡなどの微量の遺伝子を超高感度で検出・定量できることを見出し、本発明を完成した。なお、メチル化ＤＮＡと非メチル化ＤＮＡを検出という時は、検体中のＤＮＡをバイサルファイト処理した後に変換されたＤＮＡ配列をＰＣＲ増幅可能なプライマーを用いて検出することをいう。

すなわち、本発明は、次の発明［１］～［１６］を提供するものである。
［１］ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査する方法であって、前記異なる遺伝子配列を２組以上のプライマーセットを用い、前記２組以上のプライマーセットについて超高感度を得られるよう量比を変えて同時に増幅し、ＰＣＲのサイクル毎にＨＲＭ解析を行い、当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算し、スタンダードの量的数値とＣｉとの間で回帰式を算出し、そこにサンプルのＣｉを代入することでサンプルの標的遺伝子配列の量的数値を求めることを特徴とする、サンプル中の異なる遺伝子配列の存在を定量的に検査する方法。
［２］前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算が、サイクル毎のＨＲＭ解析におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出する［１］記載の検査方法。
［３］前記スタンダードのＨＲＭ解析におけるスタンダードの量的数値とＣｉとの散布図から近似式と決定係数を算出し、最も決定係数の高い近似式を採用してサンプル定量用の回帰式として使用するものである［１］又は［２］記載の検査方法。
［４］前記２組以上のプライマーセットの１組のプライマーセットの一方又は両方にタグを付与する［１］～［３］のいずれかに記載の検査方法。
［５］前記タグが、ＤＮＡ、ＲＮＡ、ＰＮＡ又はＬＮＡである［４］記載の検査方法。
［６］前記タグが、ＰＣＲ産物のＴｍ値の差を大きくするために付与するものである［４］又は［５］記載の検査方法。
［７］前記タグが、ＰＣＲ産物のＴｍ値を少なくとも２℃上昇させるタグである［６］記載の検査方法。
［８］検査対象サンプルが、全血、血漿、血清、尿、大便、唾液、手術組織、生検組織、鼻粘膜液及び脳脊髄液から選ばれる体液サンプル、並びに食品、環境水から選ばれるサンプルである［１］～［７］のいずれかに記載の検査方法。
［９］検査対象遺伝子が、正常遺伝子と異常遺伝子が混在するサンプル中の少量の異常遺伝子である［１］～［８］のいずれかに記載の検査方法。
［１０］異常遺伝子が、がん遺伝子、がん細胞特有の遺伝子異常配列、ウイルス由来遺伝子、エピゲノム異常、法医学上の重要配列、及び食品又は環境水中の特定の遺伝子配列から選ばれる遺伝子である［９］記載の検査方法。
［１１］前記プライマーセットが、メチル化ＤＮＡと非メチル化ＤＮＡを検出するプライマーセット、又は野生型遺伝子と変異型遺伝子を検出するプライマーセットである［１］～［１０］のいずれかに記載の検査方法。
［１２］前記プライマーセットがＤＮＡメチル化異常を検出するプライマーセットであり、ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１、ＢＲＣＡ１、ＳＥＰＴ９、ＭＹＯ１Ｇ及びＴＳＰＹＬ５から選ばれる遺伝子のＤＮＡメチル化異常を定量的に検出するものである［１］～［１１］のいずれかに記載の検査方法。
［１３］サーマルサイクラー、分光蛍光光度計及び融解曲線解析ソフトを具備するリアルタイムＰＣＲ装置であって、増幅サイクル毎に融解プロファイルが作動して増幅サイクル毎のＨＲＭ解析を可能とし、得られる融解曲線にピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算する計算部を有することを特徴とする、リアルタイムＰＣＲ装置。
［１４］前記計算部には、スタンダードの量的数値と前記Ｃｉとの間で回帰式を算出し、そこにサンプルのＣｉを代入することでサンプルの標的遺伝子配列の量的数値を求める計算部が含まれている［１３］記載のリアルタイムＰＣＲ装置。
［１５］前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算が、サイクル毎の融解曲線におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出する手段である［１３］又は［１４］記載のリアルタイムＰＣＲ装置。
［１６］ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１、ＢＲＣＡ１、ＳＥＰＴ９、ＭＹＯ１Ｇ及びＴＳＰＹＬ５から選ばれる遺伝子のＤＮＡメチル化異常を定量的に検出するためのプライマーセットを含有するＤＮＡメチル化異常検出用キットであって、全血、尿、手術組織及び生検組織から選ばれる検体を用いてＰＣＲにより超高感度定量可能なキット。

本発明方法によれば、遺伝子配列の違いやエピゲノム異常をより高感度で定量的に解析することが可能となる。特に、がん組織やリキッドバイオプシーサンプルに、正常な遺伝子と共に微量に混在する遺伝子変異やエピゲノム異常をより定量的に高感度に検出することができる。

ＳＥＰＴ９のＲＴ－ＨＲＭ解析用プライマー配列の位置を示す。ＦＲＭの配列上にないグレーの線はターゲット遺伝子配列とは相補性のないタグを示す。 ＭＹＯ１ＧのＲＴ－ＨＲＭ解析用プライマー配列の位置を示す。ＦＲＭの配列上にないグレー線はターゲット遺伝子配列とは相補性のないタグを示す。 ＳＥＰＴ９及びＭＹＯ１ＧのＲＴ－ＨＲＭ解析の実施例の各サイクル毎のＨＲＭ解析のｗｅｌｌ毎のパネルを示す。中央の数字がサイクル数。 ＳＥＰＴ９のＲＴ－ＨＲＭ解析の実施例のサイクル毎の％Ｍの集計表を示す。最右列はＣｉ値。 ＭＹＯ１ＧのＲＴ－ＨＲＭ解析の実施例のサイクル毎の％Ｍの集計表を示す。最右列はＣｉ値。 ＳＥＰＴ９のメチル化ＤＮＡ定量の実施例の結果を示す。左上：スタンダードのｗｅｌｌ位置、Ｓ３Ｒ Ｕｎｉｔと前頁のＣｉ値、同一サンプルのＣｉ平均値（Ａｖ Ｃｉ）、そしてＳ３Ｒ ｕｎｉｔとＣｉから得られた検量線。下の５０サイクル目のＨＲＭ解析パネル（Ｍｅｌｔｉｎｇ Ｃｕｒｖｅ Ｐａｎｅｌ）とＬａｙｏｕｔにおける濃いグレー位置は、ＨＲＭ解析でＴｍ値の異なる非特異反応が認められ、解析不可としたＷｅｌｌ位置。薄いグレーはＭアレルピークの存在が確認できず、Ｃｉ値算出不可のＷｅｌｌ位置。右上の検体定量は左上の検量線の内、平均値で得られた回帰式に検体のＣｉ値を代入してＳ３Ｒ Ｕｎｉｔを算出したもの。ＣＥＡ列は、各症例のＣＥＡ腫瘍マーカーの数値。Ｎ／Ａはデータがないもの。 ＭＹＯ１Ｇのメチル化ＤＮＡ定量の実施例の結果を示す。左上：スタンダードのｗｅｌｌ位置、Ｍ２ Ｕｎｉｔと前頁のＣｉ値、同一サンプルのＣｉ平均値（Ａｖ Ｃｉ）、そしてＭ２ ｕｎｉｔとＣｉから得られた検量線。下の５０サイクル目のＨＲＭ解析パネル（Ｍｅｌｔｉｎｇ Ｃｕｒｖｅ Ｐａｎｅｌ）とＬａｙｏｕｔにおける薄いグレーはＭアレルピークの存在が確認できず、Ｃｉ値算出不可のＷｅｌｌ位置。右上の検体定量は左上の検量線の内、平均値で得られた回帰式に検体のＣｉ値を代入してＭ２ Ｕｎｉｔを算出したもの。ＣＥＡ列は、各症例のＣＥＡ腫瘍マーカーの数値。Ｎ／Ａはデータがないもの。

本発明の一態様は、ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査する方法であって、前記の異なる遺伝子配列を２組以上のプライマーセットを用い、前記２組以上のプライマーセットの量比を変えて同時に増幅し、ＰＣＲのサイクル毎にＨＲＭ解析を行い、当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算し、スタンダードの量的数値とＣｉとの間で回帰式を算出し、サンプルのＣｉ値を代入することでサンプルの標的遺伝子配列の量的数値を求めることを特徴とする、サンプル中の異なる遺伝子配列の存在を定量的に検査する方法である。

本発明方法は、ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査する方法であり、当該検出対象となる１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列としては、下記ＤＮＡメチル化異常の外、がん細胞における塩基の挿入欠失変異（ｉｎｄｅｌ）、点突然変異、ＳＮＰｓなどの一塩基の遺伝子多型なども挙げられる。また、種々のＤＮＡ鑑定、農業分野での遺伝子改変作物やアレルゲンの検査など、配列の異なる遺伝子が混在しているサンプルにおいて、定量的かつ高感度に遺伝子の検出を行うことができる。なお、ＤＮＡメチル化異常とは、通常のメチル化状態のＤＮＡが異常を生じた結果、メチル化（ハイパーメチレーション）もしくは非メチル化（ハイポメチレーション）した状態をいう。本発明でそれらを検出するには、対象のＤＮＡをバイサルファイト処理して塩基配列の変換を行い、その変換された配列を増幅するプライマーセットを用いる。
例えば、正常遺伝子と異常遺伝子が混在するサンプル中の少量の異常遺伝子を検出する場合が挙げられ、異常遺伝子としては、がん遺伝子、がん細胞由来の遺伝子異常配列、ウイルス由来遺伝子、エピゲノム異常、法医学上の重要配列、及び食品中の特定の遺伝子配列から選ばれる遺伝子などが挙げられる。
具体的な異常遺伝子としては、遺伝子組換え食品検出について、Ｐ－３５Ｓ、Ｔ－ＮＯＳ、Ｂｔ Ｃｒｙ１Ａｂ／Ａｃなどが挙げられる。イネの品種について、Ｍｉｃｒｏｓａｔｅｌｌｉｔｅ ｍａｒｋｅｒ （ＲＭ１、ＲＭ７２、ＲＭ２０６、ＲＭ２４１、ＲＭ２８７、ＢＡＤＥＸ７－５）などが挙げられる。リステリア菌（食中毒菌）の同定について、ｓｓｒＡ ｇｅｎｅ、ｒＲＮＡなどが挙げられる。法医学でＨＲＭ解析例のある遺伝子リストとしては、精液同定；ＤＡＣＴ１、毛髪などの解析難検体からのＤＮＡ型判定（個人鑑別）；ＡＢＯ、年齢推定； ＥＬＯＶＬ２、ＦＨＬ２などが挙げられる。ＤＮＡ高分解検体で使用されるＤＮＡ型判定（個人鑑別）；ｍｉｎｉＳＴＲ （Ｄ１Ｓ１１７１、Ｄ２Ｓ１２４２、Ｄ３Ｓ１５４５、Ｄ４Ｓ２３６６、Ｄ１２Ｓ３９１、Ｄ１６Ｓ３２５３、Ｄ２０Ｓ１６１、Ｄ２１Ｓ１４３７）などが挙げられる。

がんにおいてＤＮＡメチル化異常が報告されている遺伝子としては、ＭｅｔｈＨＣ（Ａ ｄａｔａｂａｓｅ ｏｆ Ｍｅｔｈｙｌａｔｉｏｎ ａｎｄ ｇｅｎｅ ｅｘｐｒｅｓｓｉｏｎ ｉｎ Ｈｕｍａｎ Ｃａｎｃｅｒ、ｈｔｔｐ：／／ｍｅｔｈｈｃ．ｍｂｃ．ｎｃｔｕ．ｅｄｕ．ｔｗ／ｐｈｐ／ｉｎｄｅｘ．ｐｈｐ）に挙げられている遺伝子など、多数の遺伝子が報告されている。
さらに、ＳＥＰＴ９は、大腸がんに特異的にメチレーションが報告されており、血液中のＳＥＰＴ９ＤＮＡのメチレーションを検出する診断薬が販売されている。また、ＴＡＣ１、ＥＹＡ４、ＴＭＥＦＦ２、ＮＧＦＲ、ＡＤＡＭＴＳ２、ＡＤＡＭＴＳ１６、ＴＳＰＹＬ５、ＭＹＯ１Ｇなどの遺伝子も、大腸がんでメチル化されているとの報告がある。これらの遺伝子でも、メチレーションを検出することにより、大腸がんを特異的に検出することが可能となる。
また、ＥＧＦＲのチロシンキナーゼ阻害薬（ＥＧＦＲ－ＴＫＩ）など分子標的薬によって治療に抵抗性を有する変異が生じた場合の検出にも応用することができる。
さらに、ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１及びＢＲＣＡ１等から選ばれる遺伝子のメチル化ＤＮＡの検出・定量にも応用することができる。

本発明の検査は、リアルタイムＰＣＲ装置を用いて行われる。リアルタイムＰＣＲ法は、ＰＣＲの増幅量をリアルタイムでモニターし解析する方法である。この方法にはサーマルサイクラーと分光蛍光光度計を一体化したリアルタイムＰＣＲ装置が用いられる。一般的なリアルタイムＰＣＲでは、段階希釈した既知量のＤＮＡをスタンダードとしてＰＣＲを行い、これをもとに、増幅が指数関数的に起こる領域で一定の増幅産物量になるサイクル数（ｔｈｒｅｓｈｏｌｄ ｃｙｃｌｅ；Ｃｔ値）を横軸に、初発のＤＮＡ量を縦軸にプロットし、検量線を作成する。未知濃度のサンプルについても、同じ条件下で反応を行い、Ｃｔ値を求める。この値と検量線から、サンプル中の目的のＤＮＡ量を測定する。通常、リアルタイムＰＣＲのモニターには蛍光試薬が用いられるが、本発明においては、蛍光試薬としてインターカレーター試薬を用いるのが好ましい。

本発明に用いられるプライマーセットとしては、一組のＰＣＲプライマーセットの片方又は両方のプライマーにタグを付与したプライマーセットが好ましい。
塩基組成や鎖長が異なるＰＣＲ産物は、Ｔｍ値が異なる産物として増幅し、ＨＲＭ解析により定量することができる。このＴｍ値の差をより広げるために、各一組のＰＣＲプライマーセットの両方のプライマーにタグを付与するのが好ましい。

Ｔｍ値（Ｍｅｌｔｉｎｇ Ｔｅｍｐｅｒａｔｅｕｒｅ）とは、二重鎖核酸が加熱により融解する際、５０％が乖離する温度として定義される。増幅したＰＣＲ産物のＴｍ値が異なればＨＲＭ解析によって、ＰＣＲ産物を定量的に解析することができる。したがって、ＨＲＭ解析を行うことのできるＰＣＲ機器であればどのような機器を用いてもよい。例えば、リアルタイムＰＣＲ機器に付属している融解曲線解析ソフトを利用し、そこから得られたデータよりピーク高や面積を算出することで増幅されたＰＣＲ産物を定量することができる。リアルタイムＰＣＲ機器に付属している融解曲線解析ソフトであればどのようなソフトを使用してもよい。本発明方法では、得られたＰＣＲ産物を複数回ＨＲＭ解析するだけでよく、その後に電気泳動、シークエンスなどによる解析を必要としないため、簡便であるだけではなく、コンタミネーションのリスクがないという長所がある。

前記のように、本発明において、各一組のプライマーセットの片方又は両方のプライマーにタグを付与するのは、増幅する２組のＰＣＲ産物のＴｍ値の差をより広げるためである。すなわちタグは、Ｆｏｒｗａｒｄプライマー及びＲｅｖｅｒｓｅプライマーの片方又は両方に付与する。
付与するタグは、Ｔｍ値に差が生じればよく、標的核酸とのハイブリダイズを阻害しなければどのようなものを用いてもよい。具体的には、ＤＮＡ、ＲＮＡ、ＰＮＡ（ペプチド核酸）、ヌクレオチドアナログを使用したＬＮＡ（Ｌｏｃｋｅｄ ｎｕｃｌｅｉｃ ａｃｉｄｓ）などを用いることができる。また、増幅される２種のＰＣＲ産物の配列のＴｍ値が大きく異なり、ＨＲＭ解析によって分離することが可能であれば、タグを付与せずにプライマーを設計することもできる。
増幅するＰＣＲ産物のＴｍ値は、ＨＲＭ解析においてピーク分離が可能であることが必須であるが、少なくとも２℃以上、好ましくは３℃以上、より好ましくは５℃以上の差があることが好ましい。また、ＨＲＭ解析でのピーク形状や必要な感度によっては、６℃以上、より好ましくは７℃以上の差があることが好ましい。

本発明方法では、異なる遺伝子配列を２組以上の当該両側タグ付きプライマーセットを用いて同時に増幅させる。この時、２組以上のプライマーセットのプライマー比を変えてＰＣＲを行うことにより、至適な定量範囲を調節することが可能となる。したがって、異なる遺伝子配列間に極端な量的な差があっても検出・定量することができる。ダイナミックレンジが広く、大きな違いを広い範囲で定量できるリアルタイム定量ＰＣＲ法と比較すると小さな差であっても検出・定量することが可能となる。また、従来の競合ＰＣＲ法と比較すると、高感度領域で定量性をもたらすことができる。なお、この２組以上のプライマーセットの比を変えてＰＣＲを行う手段は、スタンダードにおいても、実際の検査を行う場合も、採用される。

通常、ＨＲＭ解析は、必要なＰＣＲサイクルの増幅終了後に１回のみ行うが、本発明では、複数のＰＣＲのサイクル毎にＨＲＭ解析を行うことにより、より詳細な定性的データを得られるのみならず、Ｃｉ値を求めて標準曲線を作成して定量を行うことで、ダイナミックレンジの広い定量解析を行うことができる。
ＨＲＭ解析のためのＰＣＲ機器には、ここではＱｕａｎｔ Ｓｔｕｄｉｏ ３（商標）（サーモフィッシャーサイエンティフィック株式会社）を使用しているが、同様の機器であればどのような機器を用いても解析することができる。そのような機器として、例えば、ＳｔｅｐＯｎｅ＆ＳｔｅｐＯｎｅＰｌｕｓ（商標）リアルタイムＰＣＲシステム（サーモフィッシャーサイエンティフィック株式会社）、ＬｉｇｈｔＣｙｃｌｅｒ（商標）４８０システム（ロシュ・ダイアグノスティックス株式会社）、Ｔｈｅｒｍａｌ Ｃｙｃｌｅｒ Ｄｉｃｅ（商標）Ｒｅａｌ Ｔｉｍｅ Ｓｙｓｔｅｍ ＩＩＩ（タカラバイオ株式会社）、ＣＦＸ９６ Ｔｏｕｃｈ（商標）Ｒｅａｌ－Ｔｉｍｅ ＰＣＲ Ｄｅｔｅｃｔｉｏｎ Ｓｙｓｔｅｍ（バイオ・ラッド ラボラトリーズ株式会社）がある。すでに広く普及している機器を使用することができることは、本検査方法のメリットである。

ＰＣＲ増幅に用いる試薬としては、一般に使用されるＰＣＲ試薬の他に、ＭｅｌｔＤｏｃｔｏｒ（商標）ＨＲＭ Ｍａｓｔｅｒ Ｍｉｘ、Ｆａｓｔ ＳＹＢＲ（商標）Ｇｒｅｅｎ Ｍａｓｔｅｒ Ｍｉｘ、ＰｏｗｅｒＵｐＳＹＢＲ（商標）Ｇｒｅｅｎ Ｍａｓｔｅｒ Ｍｉｘ、Ｌｕｍｉｎａｒｉｓ Ｃｏｌｏｒ ＨＲＭ Ｍａｓｔｅｒ Ｍｉｘ（以上、サーモフィッシャーサイエンティフィック株式会社）、ＫＡＰＡ ＨＲＭ ＦＡＳＴ ＰＣＲ Ｋｉｔ（ＫＡＰＡＢＩＯＳＹＳＴＥＭＳ社）、ＳｓｏＦａｓｔ ＥｖａＧｒｅｅｎ Ｓｕｐｅｒｍｉｘ（バイオ・ラッド ラボラトリーズ株式会社）等を使用することができる。ＨＲＭとしてはＥｖａＧｒｅｅｎ、ＳＹＴＯ ９、ＬＣ ＧｒｅｅｎやＭｅｌｔＤｏｃｔｏｒなどのＳａｔｕｒａｔｉｏｎ Ｄｙｅインターカレーター試薬の使用が必要となる。

本発明においては、ＰＣＲのサイクル毎にＨＲＭ解析を行い、当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算し、スタンダードの量的数値とＣｉとの間で回帰式を算出し、サンプルのＣｉ値を代入することでサンプルの標的遺伝子配列の量的数値を求める。
ここで、前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算は、例えばサイクル毎の融解曲線におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出することができる。
ＩＮＴＥＲＣＥＰＴ関数とは、単純線形回帰における回帰直線の切片を取得する方法であり、表計算ソフトエクセルに付属しているものを使用することができる。
前記スタンダードの融解曲線におけるＣｉとサンプルの融解曲線におけるＣｉとの間での回帰式の算出は、ＨＲＭ解析におけるピーク高又はピーク比率とサイクル数の散布図から決定係数を算出し、最も決定係数の高い近似式を算出することにより行うことができる。ここで、回帰式は、エクセルの散布図作成機能を使用し、最も決定係数（相関係数の２乗）の高い式を使用すればよい。

一般的なリアルタイムＰＣＲ（ｑＰＣＲ）の場合、スタンダードを用いた定量にはＣｔ値を算出して、スタンダートの数値とＣｔ値が相関することを利用して回帰曲線＝近似曲線の関数を元に算出することが一般的である（Ｘ軸にサイクル数、Ｙ軸に蛍光値）。本発明においては、前記のＣｉ値を用いた前記の回帰式を採用することにより、検体由来の混在ＰＣＲ産物の定量値を算出することが可能となった。これにより症例の経時的数値変動を追跡することが可能となった。例えば、ターゲットを腫瘍特異的ＤＮＡマーカーにしてリキッドバイオプシーを行えば、現在使用されているＣＥＡやＣＡ１９－９等の腫瘍マーカーを、感度、特異度で上回ることができる。

また、インターナルプローブを用いたｑＰＣＲ法やデジタルＰＣＲ法は、定量的かつ高感度な解析方法ではあるものの、同じ反応液中に目的外の非特異ＰＣＲ産物が発生しても検知することができない。超高感度なＰＣＲ条件では非特異反応が発生し易いが、目的のＰＣＲ産物の増幅に悪影響を及ぼしやすく、間違った結論を導かないよう除外することが望ましいにも関わらず、それらの解析方法では不可能である。本発明では、Ｔｍ値の異なる非特異ＰＣＲ産物の発生を検知してデータを除外することが可能であるため、それらのインターナルプローブを用いた解析方法より精度の高い結果を得ることが可能となる。また、高価なインターナルプローブを用いる必要がないことから、付加的にコストの必要がないという利点がある。

通常のリアルタイムＰＣＲ装置は、前述のようにサーマルサイクラーと分光蛍光光度計を一体化したリアルタイムＰＣＲ装置であり、本発明ではこれに融解曲線解析ソフトを具備するリアルタイムＰＣＲ装置が用いられる。しかし、このようなＰＣＲ装置において融解曲線解析は、必要な回数の増幅終了後に実施できるように設定されており、増幅サイクル毎に実施できるようには設定されていない。本発明方法を実施するには、増幅サイクル毎に融解曲線解析機能を作動でき、前記Ｃｉを算出できるリアルタイムＰＣＲ装置を用いるのが好ましい。
従って、本発明の別の一態様は、サーマルサイクラー、分光蛍光光度計及び融解曲線解析ソフトを具備するリアルタイムＰＣＲ装置であって、増幅サイクル毎に融解プロファイルが作動して増幅サイクル毎のＨＲＭ解析を可能とし、得られる融解曲線にピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算する計算部を有することを特徴とする、リアルタイムＰＣＲ装置である。
前記計算部には、スタンダードの量的数値と前記Ｃｉとの間で回帰式を算出し、そこにサンプルのＣｉを代入することでサンプルの標的遺伝子配列の量的数値を求める計算部が含まれているのが好ましい。
前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算は、サイクル毎の融解曲線におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出する手段が好ましい。

本発明方法に供するサンプルは、血液、組織などの臨床検体、農産物等あらゆるものから抽出した核酸を対象とすることができる。特にサンプルが、全血、血漿、血清、尿、大便、唾液、手術組織、生検組織、鼻粘膜液及び脳脊髄液から選ばれる体液サンプルである場合に好ましく用いられる。
以下に、血液サンプルによるリキッドバイオプシーを想定したＤＮＡ抽出法を記載するが、生検組織サンプル、他の試料であってもこれに準じて精製し解析を行うことができる。
新鮮なＥＤＴＡ採血検体をＱＩＡａｍｐ ＤＮＡ Ｂｌｏｏｄ Ｍａｘｉ Ｋｉｔ（キアゲン社）を使用して、末梢血ＤＮＡを抽出・精製した。なおＮｕｃｌｅｏＳｎａｐ（商標）ＤＮＡ Ｐｌａｓｍａ キット（マッハライ・ナーゲル社）、ＥＺ１ ｃｃｆＤＮＡ キット（キアゲン株式会社）、ＭａｇＭＡＸ Ｃｅｌｌ－Ｆｒｅｅ ＤＮＡ Ｉｓｏｌａｔｉｏｎ キット（サーモフィッシャーサイエンティフィック株式会社）を使用して、末梢血循環ＤＮＡ（Ｃｅｌｌ ｆｒｅｅ ＤＮＡ；ｃｆＤＮＡ）を抽出・精製してもよい。なお、上記キットは例示であり、同様のキットであればどのようなキットを使用してもよい。

次にＤＮＡのメチル化異常を検出する本発明方法の一実施態様について説明する。
ＤＮＡのメチル化異常をＰＣＲで検出するためには、ＤＮＡをバイサルファイト（亜硫酸水素ナトリウム）処理して塩基置換を利用する方法を用いるのがよい。バイサルファイト処理を行うと、ＤＮＡでは、一般にシトシン塩基（Ｃ）がウラシル塩基（Ｕ）に変換されるが、５位の炭素のメチル化修飾されたシトシン塩基は変換されないため、正常と異常のメチル化ＤＮＡでは遺伝子配列に違いを生じる。
バイサルファイト処理は、精製後のＤＮＡを、ＥｐｉＴｅｃｔ Ｂｉｓｕｌｆｉｔｅ キット（キアゲン株式会社）、ＥＺ ＤＮＡ Ｍｅｔｈｙｌａｔｉｏｎキット（ＺＹＭＯ ＲＥＳＥＡＲＣＨ 社）、ＭｅｔｈｙｌＥａｓｙ（商標）Ｘｃｅｅｄ Ｒａｐｉｄ ＤＮＡ Ｂｉｓｕｌｐｈｉｔｅ Ｍｏｄｉｆｉｃａｔｉｏｎ キット（Ｈｕｍａｎ Ｇｅｎｅｔｉｃ Ｓｉｇｎａｔｕｒｅｓ社）等を使用してＣのバイサルファイト変換を行う。ＤＮＡ配列内のＣがメチル化されていなければ脱アミノ化によりウラシル（Ｕ）塩基に変換され、Ｃの５位がメチル化されていれば、Ｕへの変換が生じない。この後のＰＣＲにより、Ｕの位置はチミン塩基（Ｔ）に置き換わるため、ＣからＴへの変換が観察される。
ＤＮＡメチル化異常はＣｐＧ配列の存在する広範な領域で生じることが多い。したがって、通常周辺のＣｐＧ配列が同時に異常をきたすため、Ｆｏｒｗａｒｄ、Ｒｅｖｅｒｓｅの双方のプライマーをＤＮＡメチル化異常が生じている配列上に設定することができる。その結果、０．０１％以下で混在する異常メチル化ＤＮＡを検出することが可能である。したがって、ＤＮＡメチル化異常は点突然変異などの遺伝子変異と比較して、非常に感度良く定量することができる。一般に、プライマーの３'側から１～数塩基の位置に、塩基配列の異なる部分が含まれるようにプライマーを設計すると感度良く異なる配列を検出することができる。

本発明では、特許文献１のＨｉＱＡＳＰで使用していたＴｍ値を上昇させるためのＰｒｉｍｅｒに付加するＴａｇは通常、ｃｔＤＮＡ（メチル化アレル）用プライマーセットの片方のみでよいが両方のＰｒｉｍｅｒに付加することで、付加しない正常ＤＮＡ（非メチル化アレル）由来ＰＣＲ産物とのＴｍ値の差が一層増大した。想定できうるこの現象に加え、ＰＣＲ産物間の競合性が薄れている現象が観察された。さらに、ＰＣＲ産物のＨＲＭ解析をＰＣＲのサイクル毎に行うと、スタンダードのピークレシオがサイクル毎に変化する現象が観察された。ここで前記Ｃｉによる回帰式を用いて解析する新アルゴリズムを適用して極微量の標的配列の定量解析が可能となった。この現象を利用すると１セットのスタンダードを１ ｒｕｎで解析するだけで、マイナーピークの現れる範囲（メジャーピークとの合計値に対する％が算出可能）を広くすることが可能となった。

一般にＰＣＲ実験系を新たに構築する際に、至適なサイクル数は経験的に何度か試験することが多い。特に高感度な解析を行いたい場合はサイクル数を最大にしたいところ、サイクル数が多すぎると時間がかかるのみならず、非特異増幅が出易くなり間違った結果を得やすくなる。そこで非特異増幅が出ていないか定性的に確認する一般的な手法にＨＲＭ解析がある。これに対し、本発明方法であれば、一般的なＰＣＲにおいても至適なサイクル数を１度の実験で幅広いサイクル数で確認できるため、ＰＣＲプログラムに組み込まれれば多くの研究者に有用な手法となる。
また、Ｃｏｖｉｄ－１９等の下水中の感染症原因ウイルスのＶａｒｉａｎｔ解析においては、多数のＶａｒｉａｎｔが予想されるところ、後期のサイクルのみでは多数あるが増幅効率の悪いＶａｒｉｎａｔが増幅の良い少数のＶａｒｉａｎｔに隠されて見られなくなる現象が想定される。これに対し、本発明方法であれば、早期サイクルでそれが確認できると考えられる。

現在より早期のがんの診断や微小残存病変（Ｍｉｎｉｍａｌ Ｒｅｓｉｄｕａｌ Ｄｅｓｅａｓｅ：ＭＲＤ）の診断が可能になれば、これまで以上に早期の治療が可能となり、治療成績の向上が期待できる。ＣＯＶＩＤ－１９パンデミックで世界に普及した安価なリアルタイムＰＣＲ装置を利用してのがん検診が可能となるため、がん患者の早期治療によるがん死亡率の低減が世界的にもたらされることが期待できる。

本発明方法によれば、サンプルとして全血を用いた場合であっても、前記メチル化ＤＮＡのような極微量の遺伝子を高感度で検出できるようになった。従って、本発明の別の一態様は、ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１、ＢＲＣＡ１、ＳＥＰＴ９、ＭＹＯ１Ｇ及びＴＳＰＹＬ５等から選ばれる遺伝子のメチル化ＤＮＡを定量的に検出するためのプライマーセットを含有するＲＴ－ＨＲＭ解析用キットであって、全血、尿、手術組織及び生検組織から選ばれる検体を用いてＰＣＲにより超高感度定量可能なキットである。

本発明のキットには、全血を用いてメチル化ＤＮＡが検出できることを記載したプロトコールが含まれることが好ましい。また、プロトコールには、本発明検査方法で検査することが記載されているのがより好ましい。

次に実施例を挙げて、本発明をさらに詳細に説明するが、本発明はこれら実施例に限定されるものではない。

［実施例１］大腸がん症例の末梢血ＤＮＡ中のがん由来ＤＮＡのメチル化ＤＮＡマーカー（ＳＥＰＴ９、ＭＹＯ１Ｇ）による定量
本発明方法により、大腸がん患者の全血ＤＮＡを用い、ＳＥＰＴ９、ＭＹＯ１Ｇのメチレーションを検出することができるか検討した。用いた試料は大腸がん症例（Ｊ１９， Ｊ２２， Ｊ２３）及び比較として市販の健常人末梢血由来ＤＮＡ（ＰｒｏＭ）を用いた。各症例の採血時期は、Ｂ１が手術直前、Ｂ２が術後１週間、Ｂ３が術後１か月である。
下記に使用したオリゴＤＮＡ配列を示す。なお、小文字配列はターゲット遺伝子配列とは相補性のないタグを示す。なお、ここでのメチル化アレル（Ｍアレル）とはＰＣＲ増幅領域におけるＣｐＧサイトのＣ（シトシン）の５'がメチル化されているエピアレルとし、非メチル化アレル（Ｕアレル）とはそれがメチル化されていないものとする。

ＳＥＰＴ９遺伝子用（Ｓ３Ｒ）：
Ｕアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ “Ｓ３Ｒ＿Ｕｆ１” ５'－ＡＧＴＴＴＡＧＴＡＴＴＴＡＴＴＴＴＴＧＡＡＧＴＴＴ（配列番号１）
Ｕアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ “Ｓ３Ｒ＿Ｕｒ３” ５'－ＣＣＡＴＣＡＴＡＴＣＡＡＡＣＣＣＣＡＣＡ（配列番号２）
Ｍアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ “Ｓ３Ｒ＿Ｓ１５Ｍｆ１” ５'－ｇｃｃｇｃｃｇｃｃｇｃｃｇｃｃＴＴＴＡＧＴＡＴＴＴＡＴＴＴＴＣＧＡＡＧＴＴＣ（配列番号３）
Ｍアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ “Ｓ３Ｒ＿Ｓ１５Ｍｒ７” ５'－ｇｃｔｇｇｔｃｔｇｇｃｔｇｇｃＣＡＴＡＴＣＧＡＡＣＣＣＣＧＣＧ（配列番号４）
Ｕアレル用スタンダードオリゴＤＮＡ “Ｓ３ＲＳｔｄ＿Ｕｎｍｅｔ” ５'－ＡＡＧＣＧＴＧＴＣＡＡＴＴＣＴＴＡＧＣＧＧＴＴＡＧＧＴＴＴＡＴＴＴＧＴＡＧＧＧＴＴＴＴＴＴＴＴＡＧＴＡＴＧＴＴＴＧＴＧＧＴＴＧＴＡＧＴＡＧＴＴＡＧＴＴＴＡＧＴＡＴＴＴＡＴＴＴＴＴＧＡＡＧＴＴＴＧＡＡＡＴＧＡＴＴＴＴＡＴＴＴＡＧＴＴＧＴＧＴＧＴＴＧＡＴＴＧＴＧＧＧＧＴＴＴＧＡＴＡＴＧＡＴＧＧＴＴＧＧＴＧＧＧＴＡＧＴＧＧＧＴＴＧＴＧＴＧＧＡＧＧＧＴＡＧＴＧＧＴＧＡＧＧＡＴＧＴＴＡＡＧＧＴＴＴＡＧＡＴＧＴＴＴＴＴＧＴＧＴＡＧＧＡＧＧＧＡＴＧＡＴＧＡＴＴＴＴＴＴＴＴＡＴＧＴＴＴＴＴＧＴＧＧＴＴＴＴＡＡＴＴＴＧＧＴＧＴＴＴＴＧＴＴＡＴＴＴＴＴＧＡＴＴＴＧＧＴＧＡＡＴＡＴＡＴＴＴＴＡＧＡＧＡＡＧＴＴＡＡＡＡＴＧＧＴＴＧＴＴＡＴＧＡＡＧＡＧＧＴＴＴＴＴＴＴＧＧＴＡＧＧＣＣＣＡＣＡＡＧＡＣＴＴＴＴ（配列番号５）
Ｍアレル用スタンダードオリゴＤＮＡ “Ｓ３ＲＳｔｄ＿Ｍｅｔ” ５'－ＧＴＴＡＧＧＴＴＴＡＴＴＣＧＴＡＧＧＧＴＴＴＴＴＴＴＴＡＧＴＡＣＧＴＴＣＧＣＧＧＴＣＧＴＡＧＴＡＧＴＴＡＧＴＴＴＡＧＴＡＴＴＴＡＴＴＴＴＣＧＡＡＧＴＴＣＧＡＡＡＴＧＡＴＴＴＴＡＴＴＴＡＧＴＴＧＣＧＣＧＴＴＧＡＴＣＧＣＧＧＧＧＴＴＣＧＡＴＡＴＧＡＴＧＧＴＴＧＧＴＧＧＧＴＡＧＣＧＧＧＴＣＧＣＧＣＧＧＡＧＧＧＴＡＧＣＧＧＣＧＡＧＧＡＣＧＴＴＡＡＧＧＴＴＴＡＧＡＣＧＴＴＴＴＣＧＴＧＴＡＧＧＡＧＧＧＡＣＧＡＣＧＡＴＴＴＴＴＴＴＴＡＣＧＴＴＴＴＣＧＴＧＧＴＴＴＴＡＡＴＴＣＧＧＣＧＴＴＴＴＧＴＴＡＴＴＴＴＴＧＡＴＴＣＧＧＴＧＡＡＴＡＴＡＴＴＴＴＡＧＡＧＡＡＧＴＴＡＡＡＡＴＧＧＴＣＧＴＴＡＣＧＡＡＧＡＧＧＴＴＴＴＴＴＴ（配列番号６）

ＭＹＯ１Ｇ遺伝子用（Ｍ２）：
Ｕアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ “Ｍ２＿Ｕｆ３” ５'－ＴＴＴＴＴＡＡＡＧＴＧＧＧＧＡＧＧＡＴＴＴＡＧＴＡＴＴＴ（配列番号７）
Ｕアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ “Ｍ２＿Ｕｒ１”
５'－ＴＡＡＣＣＣＡＴＴＴＡＣＡＣＡＣＣＡＡＡＡＡＣＡ（配列番号８）
Ｍアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ “Ｍ２＿Ｓ１５Ｍｆ３” ５'－ｃｃｇｃｃｃｃｃｃｃｃｃｃｃｃＡＡＡＧＣＧＧＧＧＡＧＧＡＴＴＴＡＧＴＡＴＴＣ（配列番号９）
Ｍアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ “Ｍ２＿Ｓ１５Ｍｒ１” ５'－ｃｇｇｃｃｇｇｃｇｃｇｃｃｔｃＣＣＧＴＴＴＡＣＡＣＧＣＣＡＡＡＡＡＣＧ（配列番号１０）
Ｕアレル用スタンダードオリゴＤＮＡ “Ｍ２Ｓｔｄ＿Ｕｎｍｅｔ” ５'－ＴＧＴＧＴＴＧＴＴＡＧＧＧＴＴＧＧＡＡＧＴＴＡＡＴＴＴＡＴＧＧＡＧＴＴＴＴＧＴＴＴＴＴＴＧＴＴＴＧＴＧＧＧＴＧＧＡＧＧＧＧＡＧＴＴＴＴＴＧＡＡＧＧＡＧＧＴＴＡＧＴＧＡＡＧＡＡＴＴＴＴＴＴＡＡＡＧＴＧＧＧＧＡＧＧＡＴＴＴＡＧＴＡＴＴＴＧＧＴＴＴＴＴＴＧＴＴＴＴＴＧＧＴＧＴＧＴＡＡＡＴＧＧＧＴＴＡＧＡＴＴＧＴＴＡＧＧＴＡＴＡＧＴＴＴＴＴＴＴＧＴＧＧＧＧＴＧＴＡＧＧＧＡＴＡＧＧＡＡＧＴＡＧＧＧＴＴＴＡＧ（配列番号１１）
Ｍアレル用スタンダードオリゴＤＮＡ “Ｍ２Ｓｔｄ＿Ｍｅｔ” ５'－ＣＧＴＧＴＴＧＴＴＡＧＧＧＴＴＧＧＡＡＧＴＴＡＡＴＴＴＡＣＧＧＡＧＴＴＴＴＧＴＴＴＴＴＴＧＴＴＴＧＴＧＧＧＴＧＧＡＧＧＧＧＡＧＴＴＴＴＣＧＡＡＧＧＡＧＧＴＴＡＧＣＧＡＡＧＡＡＴＴＴＴＴＴＡＡＡＧＣＧＧＧＧＡＧＧＡＴＴＴＡＧＴＡＴＴＣＧＧＴＴＴＴＴＣＧＴＴＴＴＴＧＧＣＧＴＧＴＡＡＡＣＧＧＧＴＴＡＧＡＴＣＧＴＴＡＧＧＴＡＴＡＧＴＴＴＴＴＴＴＧＴＧＧＧＧＣＧＴＡＧＧＧＡＴＡＧＧＡＡＧＴＡＧＧＧＴＴＴＡＧ（配列番号１２）

Ｕアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ＋Ｕアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ＝各５μＭ：ＦＲ－Ｕ
Ｍアレル用Ｆｏｒｗａｒｄ Ｐｒｉｍｅｒ＋Ｍアレル用Ｒｅｖｅｒｓｅ Ｐｒｉｍｅｒ＝各５μＭ：ＦＲ－Ｍ
予めＦＲ－ＵとＦＲ－Ｍを調整しておき、ＰＣＲ時にＦＲ－Ｕ：ＦＲ－Ｍ比をＳＥＰＴ９では１．７５：１、ＭＹＯ１Ｇでは１：２．５で使用した（今回のデータを得る際に使用した比率であり、必ずしも固定したものではない。通常は各アッセイ構築の際に各種比率での検討を行い、高感度定量に関して至適な比率を使用する）。

説明と同意を得た大腸がん患者末梢血検体は、リキッドバイオプシー用採血管「ＰＡＸｇｅｎｅ^TM ｃｃｆＤＮＡ採血管１０ｍＬ」（ＢＤ社）にて、採血して、「ＱＩＡａｍｐ ＤＮＡ Ｂｌｏｏｄ Ｍａｘｉ Ｋｉｔ」（キアゲン社）にてＤＮＡを抽出精製した。その後「ＥＺ ＤＮＡ Ｍｅｔｈｙｌａｔｉｏｎ－Ｌｉｇｈｔｎｉｎｇ ｋｉｔ」（ＺｙｍｏＲｅｓｅａｒｃｈ社）によりバイサルファイト処理し、解析試料とした。
１ウェルあたりの鋳型ＤＮＡ量は、検体の解析には元のゲノムＤＮＡ量約１００ｎｇを出発材料としたバイサルファイト処理試料を用い、スタンダードオリゴＤＮＡについてはＵアレル用スタンダードオリゴＤＮＡがすべて３０，３００コピー相当（ヒトゲノム１００ｎｇ当たりのコピー数）、Ｍアレル用スタンダードオリゴＤＮＡは各々のｗｅｌｌレイアウト記載の数値を使用した（三重測定）。ＳＥＰＴ９とＭＹＯＧ１のＰＣＲ増幅領域は正常ＤＮＡではＵアレル、大腸がんではＭアレルであるため全血由来ＤＮＡから微量ｃｔＤＮＡを定量するために、Ｍアレル用スタンダードオリゴは３２０コピー相当を最大量として２倍毎の希釈段階を使用した（単位は、それぞれのオリゴＤＮＡ濃度をコピー数換算したもので、ＳＥＰＴ９についてはＳ３Ｒ Ｕｎｉｔ、ＭＹＯ１ＧについてはＭ２ Ｕｎｉｔと表記）。スタンダードオリゴの希釈には、１μｇ／ｍｌ ｏｆ Ｃａｒｒｉｅｒ ＲＮＡ（Ｑｉａｇｅｎ社）ｉｎ ＩＤＴ ｂｕｆｆｅｒ （ＩＤＴ社）を使用した。市販の健常人末梢血由来ＤＮＡはプロメガ社のＨｕｍａｎ Ｇｅｎｏｍｉｃ ＤＮＡ （Ｇ１４７１）を使用した。

ＰＣＲ反応液の組成を以下に記載する。ただし本発明方法（Ｒｅａｌ―Ｔｉｍｅ－ＨＲＭ＝ＲＴ－ＨＲＭ）を実施するにあたりＰＣＲ機器とＰＣＲ試薬等、濃度や量比、温度は下記に限定されるものではなく、一般的なＰＣＲに準じる。濃度はＰＣＲ反応液中の最終濃度を示す。なお以下の記載以外はメーカーの取り扱い説明書に従った。
ＰＣＲ試薬；ＥｐｉＴａｑ^TM ＨＳ（Ｔａｋａｒａ－ｂｉｏ社）。
ＭｇＣｌ₂濃度：２．５μＭ
ｄｄＮＴＰ濃度：２００μＭ ｅａｃｈ
各Ｐｒｉｍｅｒの合計濃度：１．４μＭ
蛍光試薬；２０ｘＥｖａＧｒｅｅｎ^TM（Ｂｉｏｔｉｕｍ， Ｉｎｃ）濃度：１ｘ
蛍光リアルタイムＰＣＲ装置にはＱｕａｎｔ Ｓｔｕｄｉｏ ３を、解析ソフトには機器付属のＱｕａｎｔＳｔｕｄｉｏ^TM Ｄｅｓｉｇｎ＆Ａｎａｌｙｓｉｓ Ｓｏｆｔｗａｒｅ ｖ．１．４．３を使用した。ＰＣＲは、最初の３０サイクルは通常のＨＲＭ解析とし、その後のサイクルでサイクル毎のＨＲＭデータ取得を行い、取得したデータを集計して、本発明方法によるＲＴ－ＨＲＭ解析を行った。以下にＰＣＲとＲＴ－ＨＲＭ解析条件を示す。

Ｄｅｎａｔｕｒｅ：９５℃ ｆｏｒ ２ ｍｉｎ －＞ＰＣＲ：（９５℃ ｆｏｒ ５ ｓｅｃ－Ｔ_AE ｆｏｒ １５ ｓｅｃ）３０ ｃｙｃｌｅｓ－＞Ｍｃ：６８℃－９２℃ ａｔ ０．１℃／ｓｅｃ－Ｔ_AE ｆｏｒ ２０ ｓｅｃ－＞Ｄａｔａ ｅｘｐｏｒｔ ｆｏｒ ３０ ｃｙｃｌｅｓ

次に、本発明方法（ＲＴ－ＨＲＭ）の特徴であるサイクル毎のＨＲＭデータ取得は以下を行い、データのエキスポートをマニュアルでサイクル毎に繰り返し実施した。
Ｍｃ：６８℃－９２℃ ａｔ ０．１℃／ｓｅｃ－Ｔ_AE ｆｏｒ ２０ ｓｅｃ－＞Ｄａｔａ ｅｘｐｏｒｔ
Ｔ_AE：Ａｎｅａｌｉｎｇ／Ｅｘｔｅｎｓｉｏｎ ｔｅｍｐｅｒａｔｕｒｅ：５８℃ ｆｏｒ ＳＥＰＴ９，６４℃ ｆｏｒ ＭＹＯ１Ｇ
なお、複数サイクル毎の本発明方法（ＲＴ－ＨＲＭ）データ取得は最初のＸサイクルのＰＣＲ、Ｍｃ後、以下のサイクル数＝Ｙを調整して、データのエキスポートを繰り返すことで実施可能である。
ＰＣＲ：（９５℃ ｆｏｒ ５ ｓｅｃ－Ｔ_AE ｆｏｒ １５ ｓｅｃ） Ｙ ｃｙｃｌｅｓ－＞Ｍｃ：６８℃－９２℃ ａｔ ０．１℃／ｓｅｃ－Ｔ_AE ｆｏｒ ２０ ｓｅｃ－＞Ｄａｔａ ｅｘｐｏｒｔ：この場合Ｘ＋１＋Ｙサイクルのデータが取得される。この工程は必要回数実施可能で、Ｙ部は用途に応じて可変。また、Ｍｃ部の０．１℃／ｓｅｃも可変である。
マニュアルでのデータ取得と後述するデータ解析については、機器内のソフトをプログラミングして組み込めば手間のかからない解析となる。
図１に、ＳＥＰＴ９ メチル化ＤＮＡ（大腸がんにおけるｃｔＤＮＡ）ＲＴ－ＨＲＭ解析用プライマー配列の位置を示す。図２に、ＭＹＯ１Ｇ メチル化ＤＮＡ（大腸がんにおけるｃｔＤＮＡ）ＲＴ－ＨＲＭ解析用プライマー配列の位置を示す。

図３はＳＥＰＴ９及びＭＹＯ１Ｇについてサイクル毎のＨＲＭグラフを示す。中央の数字がサイクル数を示す。１マスが１ ＷｅｌｌのＨＲＭパターンであり、ｘ軸が温度、ｙ軸は－ｄＲＦＵ／ｄｔ（温度上昇に対する蛍光シグナル減少変化の負の微分値）。サイクル毎の各アレルの存在比が変化する様子が確認でき、Ｅｎｄｐｏｉｎｔサイクルのみの従来のＨＲＭでは試料内の正確な定性的、定量的情報を得られていないことが理解できる。またサイクル毎の変化が定量的であることは、段階的な融解曲線のみならず、次に示す新たに考案したＣｉ（Ｉｎｔｅｒｃｅｐｔ Ｃｙｃｌｅ）値とスタンダードによる検量線で明確に確認できる。

サイクル毎にＥｘｐｏｒｔしたＭｅｌｔ Ｃｕｒｖｅ Ｄｅｒｉｖａｔｉｖｅ データを集計して、ＵアレルＰＣＲ産物とＭアレルＰＣＲ産物のＴｍ値（ＰＣＲ産物の融解温度；－ｄＲＦＵ／ｄｔ値が最大値を取る温度）に対して、サイクル毎に－ｄＲＦＵ／ｄｔ値のピーク高の数値を算出した。その際、バックグラウンドは各ピークの中間域の温度の数値を使用して各ピーク数値から差し引いた。ピークの比率はピーク数値の和を１００％として求めた。Ｍアレルの比率＝％Ｍをサイクル毎にプロットして、％Ｍが最初に０％でなくなるサイクルの数値と次のサイクルの数値とでエクセルのＩＮＴＥＲＣＥＰＴ関数を使用し、サイクルをＹ、％ＭをＸとしてＹ切片を求めＣｉ値とした（図４及び図５）。このＣｉ値により検量線を作成したものが、図６及び図７であり、ＨＲＭ解析により非特異シグナルと判定したｗｅｌｌのデータを除いた検量線において、決定係数０．９以上の良好な相関を示した。ゆえに新たに開発したＣｉを用いたＲＴ－ＨＲＭ解析の定量アルゴリズムは良好に機能することを示した。同一ロットのスタンダードを使用して本手法を実施すれば検体採取時期の異なる検体の経時的な解析により、数値変動の追跡が可能となるため、新たな腫瘍マーカー解析手法として有用と考えられる。

ここで示すＳＥＰＴ９とＭＹＯ１Ｇのデータにおいて、上記のＴｍ値のピークは、Ｕアレル由来のピークとＭアレル由来のピークの２種類を対象とし、その２種類以外のピークは非特異ＰＣＲ産物として当該ｗｅｌｌのデータを除外対象とした（非特異シグナルを判別できるＨＲＭ解析ならではの利点）。図６のパネルにはＭアレル由来ピークに外れＴｍが認められるｗｅｌｌ（濃いグレーｗｅｌｌ）があり、また図６と図７でＭアレル由来ピークが認められないｗｅｌｌ（薄いグレーｗｅｌｌ）があり、それらのデータは上記のように定量解析から除外した。
なおＳＥＰＴ９のＵアレル由来ピークのＴｍ分布は７３．２℃－７４．３℃、Ｍアレル由来ピークのＴｍ分布は８４．０℃－８５．０℃で、ＭＹＯ１ＧのＵアレル由来ピークのＴｍ分布は７３．１℃－７４．２℃、Ｍアレル由来ピークのＴｍ分布は８６．６℃－８７．３℃であった。

Claims (16)

1. ＰＣＲによって１種又は２種以上の遺伝子領域に混在する異なる遺伝子配列を検査する方法であって、前記異なる遺伝子配列を２組以上のプライマーセットを用い、前記２組以上のプライマーセットの量比を変えて同時に増幅し、ＰＣＲのサイクル毎にＨＲＭ解析を行い、当該ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算し、スタンダードの量的数値とＣｉとの間で回帰式を算出し、サンプルのＣｉ値を代入することでサンプルの標的遺伝子配列の量的数値を求めることを特徴とする、サンプル中の異なる遺伝子配列の存在を定量的に検査する方法。
2. 前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算が、サイクル毎のＨＲＭ解析におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出する請求項１記載の検査方法。
3. 前記スタンダードのＨＲＭ解析におけるスタンダードの量的数値とＣｉとの散布図から近似式と決定係数を算出し、最も決定係数の高い近似式を採用してサンプル定量用の回帰式として使用するものである請求項１又は２記載の検査方法。
4. 前記２組以上のプライマーセットの１組のプライマーセットの一方又は両方にタグを付与する請求項１～３のいずれか１項記載の検査方法。
5. 前記タグが、ＤＮＡ、ＲＮＡ、ＰＮＡ又はＬＮＡである請求項４記載の検査方法。
6. 前記タグが、ＰＣＲ産物のＴｍ値の差を大きくするために付与するものである請求項４又は５記載の検査方法。
7. 前記タグが、ＰＣＲ産物のＴｍ値を少なくとも２℃上昇させるタグである請求項６記載の検査方法。
8. 検査対象サンプルが、全血、血漿、血清、尿、大便、唾液、手術組織、生検組織、鼻粘膜液及び脳脊髄液から選ばれる体液サンプルである請求項１～７のいずれか１項記載の検査方法。
9. 検査対象遺伝子が、正常遺伝子と異常遺伝子が混在するサンプル中の少量の異常遺伝子である請求項１～８のいずれか１項記載の検査方法。
10. 異常遺伝子が、がん遺伝子、がん細胞特有の遺伝子異常配列、ウイルス由来遺伝子、エピゲノム異常、法医学上の重要配列、及び食品又は環境水中の特定の遺伝子配列から選ばれる遺伝子である請求項９記載の検査方法。
11. 前記プライマーセットが、メチル化ＤＮＡと非メチル化ＤＮＡを検出するプライマーセット、又は野生型遺伝子と変異型遺伝子を検出するプライマーセットである請求項１～１０のいずれか１項記載の検査方法。
12. 前記プライマーセットがＤＮＡメチル化異常を検出するプライマーセットであり、ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１、ＢＲＣＡ１、ＳＥＰＴ９、ＭＹＯ１Ｇ及びＴＳＰＹＬ５から選ばれる遺伝子のＤＮＡメチル化異常を定量的に検出するものである請求項１～１１のいずれか１項記載の検査方法。
13. サーマルサイクラー、分光蛍光光度計及び融解曲線解析ソフトを具備するリアルタイムＰＣＲ装置であって、増幅サイクル毎に融解プロファイルが作動して増幅サイクル毎のＨＲＭ解析を可能とし、得られる融解曲線にピークが生じ始める直前のサイクル数換算値（Ｃｉ）を計算する計算部を有することを特徴とする、リアルタイムＰＣＲ装置。
14. 前記計算部には、スタンダードの量的数値と前記Ｃｉとの間で回帰式を算出し、そこにサンプルのＣｉを代入することでサンプルの標的遺伝子配列の量的数値を求める計算部が含まれている請求項１３記載のリアルタイムＰＣＲ装置。
15. 前記ＨＲＭ解析においてピークが生じ始める直前のサイクル数換算値（Ｃｉ）の計算が、サイクル毎の融解曲線におけるピーク比率（Ｙ）とサイクル数（Ｘ）から、ＴｈｒｅｓｈｏｌｄラインをＹ軸＝０上に置き、ＩＮＴＥＲＣＥＰＴ関数を用いてＹ切片（Ｙ値＝０）のＸ軸値（Ｃｉ）を算出する手段である請求項１３又は１４記載のリアルタイムＰＣＲ装置。
16. ＣＤＫＮ２Ａ、ＭＧＭＴ、ＴＦＰＩ２、ＮＤＲＧ４、ＳＤＣ２、ＢＭＰ３、ＰＲＩＭＡ１、ＳＦＲＰ２、ＶＩＭ、ＣＤＨ１、ＣＵＸ２、ＲＥＧ１Ａ、ＣＣＮＤ２、ＰＬＡＵ、ＳＯＣＳ１、ＴＨＢＳ、ＶＨＬ、ＲＡＳＳＦ１、ＺＮＦ１５４、ＢＮＣ１、ＳＳＴ１、ＭＬＨ１、ＢＲＣＡ１、ＳＥＰＴ９、ＭＹＯ１Ｇ及びＴＳＰＹＬ５から選ばれる遺伝子のＤＮＡメチル化異常を定量的に検出するためのプライマーセットを含有するＤＮＡメチル化異常検出用キットであって、全血、尿、手術組織及び生検組織から選ばれる検体を用いてＰＣＲにより超高感度定量可能なキット。