(57)【特許請求の範囲】
【請求項１】
以下の工程（Ａ）〜（Ｅ）を備えたことを特徴とする子宮平滑筋における腫瘍を診断するためのデータを収集する方法。
（Ａ）採取された子宮平滑筋腫、子宮平滑筋肉腫、及び正常子宮平滑筋層の組織からゲノムＤＮＡを抽出する工程；
（Ｂ）工程（Ａ）で抽出した各組織のゲノムＤＮＡ中の、ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域におけるＣｐＧ配列のメチル化頻度を測定する工程；
（Ｃ）被検体から採取された子宮平滑筋における腫瘍組織からゲノムＤＮＡを抽出する工程；
（Ｄ）工程（Ｃ）で抽出したゲノムＤＮＡ中の、ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域をコードする領域におけるＣｐＧ配列のメチル化頻度を測定する工程；
（Ｅ）前記被検体における子宮平滑筋腫及び子宮平滑筋肉腫を診断するためのデータを収集するため、工程（Ｄ）で測定したＣｐＧ配列のメチル化頻度データを、基準となる工程（Ｂ）で測定したＣｐＧ配列のメチル化頻度データとともに階層的クラスタリングを行う工程；
【請求項２】
工程（Ｂ）及び（Ｄ）において、
ＡＬＸ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列であり、
ＣＢＬＮ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列であり、
ＣＯＲＩＮ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列であり、
ＦＯＸＰ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列であり、ＧＡＴＡ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列であり、
ＩＧＬＯＮ５遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列であり、
ＮＰＴＸ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列であり、
ＮＴＲＫ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列であり、ＰＲＬ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列であり、ＳＴＥＡＰ４遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列であることを特徴とする請求項１記載の方法。
【請求項３】
ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域におけるＣｐＧ配列のメチル化頻度を測定するためのプライマー若しくはプローブ、又はそれらの標識物を備えたことを特徴とする、子宮平滑筋における腫瘍の診断用キット。
【請求項４】
ＡＬＸ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列であり、
ＣＢＬＮ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列であり、
ＣＯＲＩＮ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列であり、
ＦＯＸＰ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列であり、ＧＡＴＡ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列であり、
ＩＧＬＯＮ５遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列であり、
ＮＰＴＸ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列であり、
ＮＴＲＫ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列であり、ＰＲＬ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列であり、ＳＴＥＡＰ４遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列であることを特徴とする請求項３記載のキット。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、子宮平滑筋における腫瘍を診断するためのデータを収集する方法や子宮平滑筋における腫瘍の診断用キットに関する。
【背景技術】
【０００２】
子宮肉腫は子宮平滑筋から発生する悪性腫瘍であり、比較的早期の症例でも肺や肝臓に転移を起こす極めて悪性度が高い難治性の腫瘍である。子宮平滑筋に由来する腫瘍には、良性である子宮平滑筋腫（子宮筋腫）と、悪性である子宮平滑筋肉腫（子宮肉腫）がある。子宮肉腫は、子宮筋腫と比較しその発症頻度は１％以下であるが、５年生存率が約３７％であり、極めて予後不良な腫瘍である。子宮肉腫に対して、既存の化学療法は殆ど延命効果が望めず、有効な治療法は現在、外科的手術以外に確立されていない。
【０００３】
一方、子宮筋腫は良性腫瘍であり、性成熟期女性の３０％以上に認められる最も発症頻度の高い腫瘍である。近年、晩婚化や妊娠高齢化といった女性のライフスタイルの変化により、子宮筋腫治療の方法として、妊娠可能な状態で子宮を温存する子宮筋腫摘出手術（子宮核出術）が選択される機会が増えている。
【０００４】
子宮肉腫は子宮筋腫に比較し発症頻度は低いが、両者は発生場所や形状が類似しているため、患者の子宮温存を決定する際には腫瘍の良性若しくは悪性の診断が非常に重要となる。しかしながら、現行の診断は病理標本における分裂細胞数、細胞異型、及び壊死の有無の３要因を元に総合的に判断するという形態学的基準に依っている。
形態学的基準は施設或いは診断した個人により判断が異なる場合があるため、客観性のある信頼性の高い診断法が望まれている。しかし、子宮平滑筋腫と子宮平滑筋肉腫を区別する有用なバイオマーカーは現状存在しない。
【０００５】
現在、最も信頼性の高い子宮筋腫バイオマーカーはＭＥＤ１２遺伝子の突然変異であり、約７０％の子宮筋腫に遺伝子変異が検出されるという報告がある（非特許文献１：Science. 2011 Oct 14;334(6053):252-5.）。また、子宮肉腫が発症することにより発現が低下することが見出されているＬＭＰ２タンパクの免疫染色による診断が報告されている（非特許文献２：Sci Rep. 2011;1:180.）。しかしこれらの判定方法はいずれも、実際の診断に用いるには客観性及び信頼性に著しく欠ける。
【０００６】
一方、近年、エピジェネティクス（クロマチンへの後天的な修飾により、遺伝子配列の変化を伴わずに遺伝子発現が制御されることに起因する遺伝学又は分子生物学の研究分野）への関心が高まり、がん等の疾患との関連についての研究が進められている。中でも、ゲノムＤＮＡ上の遺伝子のメチル化が注目されている。高等真核生物のゲノムＤＮＡ配列中に存在する５’−ＣＧ−３’ＤＮＡ（以下、ＣｐＧ）部位では、グアニン（Ｇ）の５’側に位置するシトシン（Ｃ）がメチル化修飾される現象が知られており、このＣｐＧのメチル化修飾は、遺伝子発現に影響を及ぼすと考えられている。特にＣｐＧ部位に富む領域（ＣｐＧアイランド）が遺伝子のプロモーター領域内に存在する場合には、遺伝子発現に対して重要な影響を及ぼす。通常、ゲノム上の多くの遺伝子はこれらのメチル化修飾から保護されているが、何らかの原因により、遺伝子のプロモーター領域に存在するＣｐＧアイランドがメチル化された場合、遺伝子の転写が抑制されることになる。このため、ＣｐＧアイランドのメチル化異常によって、例えば、ヒト生体内におけるがん抑制遺伝子の転写が不活性化された場合、細胞増殖の制御が効かなくなり、がんなどの細胞増殖性疾患が進行してしまうことになる。実際、いくつかのがんにおいて、特定の遺伝子におけるメチル化頻度が上昇していることが報告されており、最近では、特定の遺伝子のＣｐＧアイランドのメチル化頻度を、特定のがんの診断に利用する試みがいくつか行われている。例えば、ＢＡＳＰ１等の遺伝子のＣｐＧアイランドのメチル化の程度を検出することによって、肝臓がんを診断する方法が知られている（特許文献１：特開２００８−２４５６３５号公報）。また、子宮がん等の婦人科がんにおいては、ｈＭＬＨ−1、ＣＤＫＮ２Ａ／ｐ１６の各遺伝子におけるメチル化頻度が上昇しているとの報告がなされている（非特許文献３：Eur J Gynaecol Oncol. 2009; 30:267-270.）。更に、ＤＣＣ、ＧＨＳＲ、ＡＪＡＰ１、ＺＮＦ５６０、ＦＥＲＤ３Ｌ、ＦＬＪ２３５１４、ＢＸＤＣ１、ＨＫＲ１、ＤＲＤ４、ＰＥＮＫ、ＦＡＭ１２Ｂ、ＨＩＳＴ１Ｈ４Ｆ、ＮＥＦ３、ＳＫＩＰ、ＣＸ３６、ＥＯＭＥＳ、及び、ＳＯＲＣＳ３からなる遺伝子群から選ばれる１又は２種以上の遺伝子のＣｐＧ部位におけるメチル化の頻度を、子宮体がんの判定に用いる方法が報告されている（特許文献２：特開２０１１−１６０７１１号公報）。しかし、がんにおける遺伝子のメチル化頻度の上昇の程度は遺伝子等によって様々であり、実用的な婦人科がんの診断に耐え得る程度のメチル化頻度の上昇が生じる遺伝子はこれまで知られていない。また、本発明者らはこれまで、子宮筋腫と子宮平滑筋（正常筋層）のゲノムワイドなＤＮＡメチル化解析（ＤＮＡメチローム）を行い、ＤＮＡメチル化がｍＲＮＡ発現と比較し、より明確に正常筋層と子宮筋腫を区別できることを見出している（非特許文献４：PLoS One, 8(6): e66632, 2013）が、これは正常筋層と子宮筋腫とを区別するものであり、子宮平滑筋腫と子宮平滑筋肉腫とを区別するものではない。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００８−２４５６３５号公報
【特許文献２】特開２０１１−１６０７１１号公報
【非特許文献】
【０００８】
【非特許文献１】Science. 2011 Oct 14;334(6053):252-5.
【非特許文献２】Sci Rep. 2011;1:180.
【非特許文献３】Eur J Gynaecol Oncol. 2009; 30:267-270.
【非特許文献４】PLoS One, 8(6): e66632, 2013
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の課題は、子宮平滑筋に生じた腫瘍が良性腫瘍（子宮平滑筋腫）か、悪性腫瘍（子宮平滑筋肉腫）のいずれかであるかを精度よく診断できる診断マーカーを提供することにある。
【課題を解決するための手段】
【００１０】
発明者らはこれまで、子宮筋腫と子宮平滑筋（正常筋層）のゲノムワイドなＤＮＡメチル化解析（ＤＮＡメチローム）を行い、ＤＮＡメチル化がｍＲＮＡ発現と比較し、より明確に正常筋層と子宮筋腫を区別できることを見いだしている。これらＤＮＡメチル化が起こった遺伝子からまず、メチル化解析に基づき正常筋層と子宮筋腫とが判別可能かを検討したが、１遺伝子のメチル化レベル（度）を検出しただけでは判別することは難しいことが明らかとなった。そこで、判別精度を保ちつつメチル化解析の労力を最小限にできる遺伝子の組み合わせを検討したところ、ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４の１０遺伝子（以下、これらを総称して「本件バイオマーカー遺伝子」ということがある）をコードする領域又は前記各遺伝子の転写制御領域をコードする領域（以下、「本件バイオマーカー遺伝子領域」ということがある）におけるＣｐＧ配列のメチル化頻度を測定し、多変量解析を行い、ＣｐＧメチル化に関する本件バイオマーカー遺伝子パネルを作製することにより、正常筋層と子宮筋腫とが明瞭に判別できるのみならず、腫瘍が子宮平滑筋腫か子宮平滑筋肉腫のいずれかであるかを判定できることを見いだし、本発明を完成するに至った。
【００１１】
すなわち、本発明は、以下のとおりである。
［１］以下の工程（Ａ）〜（Ｅ）を備えたことを特徴とする子宮平滑筋における腫瘍を診断するためのデータを収集する方法。
（Ａ）採取された子宮平滑筋腫、子宮平滑筋肉腫、及び正常子宮平滑筋層の組織からゲノムＤＮＡを抽出する工程；
（Ｂ）工程（Ａ）で抽出した各組織のゲノムＤＮＡ中の、ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域におけるＣｐＧ配列のメチル化頻度を測定する工程；
（Ｃ）被検体から採取された腫瘍組織からゲノムＤＮＡを抽出する工程；
（Ｄ）工程（Ｃ）で抽出したゲノムＤＮＡ中の、ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域をコードする領域におけるＣｐＧ配列のメチル化頻度を測定する工程；
（Ｅ）工程（Ｂ）で測定したＣｐＧ配列のメチル化頻度、及び工程（Ｄ）で測定したＣｐＧ配列のメチル化頻度に基づき多変量解析を行い、前記被検体における子宮平滑筋腫及び子宮平滑筋肉腫を診断するためのデータを収集する工程；
［２］多変量解析が階層的クラスタリング解析であることを特徴とする上記［１］記載の方法。
［３］工程（Ｂ）及び（Ｄ）において、
ＡＬＸ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列であり、
ＣＢＬＮ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列であり、
ＣＯＲＩＮ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列であり、
ＦＯＸＰ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列であり、
ＧＡＴＡ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列であり、
ＩＧＬＯＮ５遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列であり、
ＮＰＴＸ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列であり、
ＮＴＲＫ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列であり、
ＰＲＬ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列であり、
ＳＴＥＡＰ４遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列であることを特徴とする上記［１］又は［２］記載の方法。
［４］ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４遺伝子をコードする領域又は前記各遺伝子の転写制御領域におけるＣｐＧ配列のメチル化頻度を測定するためのプライマー若しくはプローブ、又はそれらの標識物を備えたことを特徴とする、子宮平滑筋における腫瘍の診断用キット。
［５］ＡＬＸ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列であり、
ＣＢＬＮ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列であり、
ＣＯＲＩＮ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列であり、
ＦＯＸＰ１遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列であり、
ＧＡＴＡ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列であり、
ＩＧＬＯＮ５遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列であり、
ＮＰＴＸ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列であり、
ＮＴＲＫ２遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列であり、
ＰＲＬ遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列であり、
ＳＴＥＡＰ４遺伝子をコードする領域又は前記遺伝子の転写制御領域におけるＣｐＧ配列が、配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列であることを特徴とする上記［４］記載のキット。
【００１２】
上記本発明の方法は、医師による子宮平滑筋における腫瘍の診断を補助する方法であって、医師による診断行為を含まない。
【００１３】
また、本発明の実施の他の形態として、子宮平滑筋における腫瘍の診断方法を挙げることができる。
【発明の効果】
【００１４】
患者の子宮温存を決定する際には腫瘍の良性若しくは悪性の診断が非常に重要となるが、子宮平滑筋腫と子宮平滑筋肉腫を区別する有用なバイオマーカーは存在しなかった。ＣｐＧメチル化に関する本件バイオマーカー遺伝子パネルを用いることにより、平滑筋由来の腫瘍が生じた患者が子宮平滑筋腫か、子宮平滑筋肉腫のいずれかを罹患しているかについて精度よく判定することができる。
【図面の簡単な説明】
【００１５】
【図１】各１８症例の正常筋層検体、及び子宮筋腫検体を用い、ＡＬＸ１遺伝子のメチル化頻度を解析した結果である。
【図２】各１８症例の正常筋層検体、及び子宮筋腫検体について、本件バイオマーカー遺伝子領域のＣｐＧメチル化データを取得し、階層的クラスタリングを行った結果である。
【図３】１２症例の子宮肉腫検体、及び肉腫細胞株４サンプル（ヒト子宮肉腫細胞株ＳＫＮ、及びヒト卵巣肉腫細胞株ＲＫＮをそれぞれ２サンプルずつ）について、ＣＯＢＲＡ法により本件バイオマーカー遺伝子領域のＣｐＧメチル化データを取得し、図２にて構築した各１８症例の正常筋層検体、及び子宮筋腫検体における階層的クラスタリングの結果と共に示した結果である。
【図４】新規な２症例（検体Ｘ、及びＹ）について、ＣＯＢＲＡ法により本件バイオマーカー遺伝子領域のＣｐＧメチル化データを取得し、そのデータを図３にて示した肉腫、筋腫、及び筋層検体のメチル化データに追加して階層的クラスタリングを行った結果である。
【発明を実施するための形態】
【００１６】
本発明の子宮平滑筋における腫瘍を診断するためのデータを収集する方法（以下、「本件収集方法」ということがある）としては、子宮平滑筋腫及び子宮平滑筋肉腫、並びに正常子宮平滑筋層の組織からゲノムＤＮＡを抽出する工程（Ａ）；工程（Ａ）で抽出した各組織のゲノムＤＮＡ中の本件バイオマーカー遺伝子領域におけるＣｐＧ配列（部位）（以下、「本件ＣｐＧ部位」とも表示する。）がメチル化されている度合い（レベル又は頻度）を測定する工程（Ｂ）；被検体から採取された腫瘍組織からゲノムＤＮＡを抽出する工程（Ｃ）；工程（Ｃ）で抽出したゲノムＤＮＡ中の本件バイオマーカー遺伝子領域におけるＣｐＧ配列がメチル化されている頻度を測定する工程（Ｄ）；工程（Ｂ）で測定したメチル化頻度、及び工程（Ｄ）で測定したメチル化頻度に基づき多変量解析を行い、前記被検体における子宮平滑筋腫及び子宮平滑筋肉腫を診断するためのデータを収集する工程（Ｅ）；の工程（Ａ）〜（Ｅ）を備えた方法であれば特に制限はされない。子宮平滑筋腫の組織由来の本件ＣｐＧ部位は、正常子宮平滑筋層の組織由来の本件ＣｐＧ部位と比べメチル化頻度が変化し、また、子宮平滑筋肉腫組織由来の本件ＣｐＧ部位は、子宮平滑筋腫の組織由来の本件ＣｐＧ部位とは異なるメチル化頻度を示すため、上記工程（Ｅ）において多変量解析を行い、ＣｐＧメチル化に関する本件バイオマーカー遺伝子パネルを作成すると、メチル化頻度を指標とし腫瘍が子宮平滑筋腫か子宮平滑筋肉腫のいずれかであるかを判定することができる。なお、本明細書における本件バイオマーカー遺伝子領域におけるＣｐＧ部位とは、ゲノムＤＮＡ上で該ＣｐＧ部位に最も近い位置に存在する遺伝子が、その本件バイオマーカー遺伝子領域の遺伝子であることを意味する。また、本件ＣｐＧ部位が該遺伝子の近傍に複数存在し、ＣｐＧアイランドを形成している場合は、それら複数のＣｐＧ部位を子宮平滑筋腫と子宮平滑筋肉腫とを区別するバイオマーカーとして使用することが好ましいが、近傍のＣｐＧ部位のメチル化頻度は、相関が見られるとされているため、特定の１カ所のＣｐＧ部位のメチル化頻度が増加していれば、該ＣｐＧ部位が形成するＣｐＧアイランドにおけるメチル化頻度も増加していると評価することができる。
【００１７】
本発明において、「ＣｐＧ配列がメチル化されている」とは、シトシン（Ｃ）の次にグアニン（Ｇ）が現れる２塩基配列（ＣＧ又はＧＣ配列）において、シトシン塩基５位の炭素がメチル化された状態を意味する。
【００１８】
本件収集方法において被検体の生物種としては、ヒト、マウス、ラット、ハムスター、モルモット、サル、ウシ、ブタ、ウマ、ウサギ、ヒツジ、ヤギ、ネコ、イヌ等の哺乳動物を好適に例示することができ、中でもヒト、マウスをより好適に例示することができ、特にヒトを好適に例示することができる。
【００１９】
上記（Ａ）工程及び（Ｃ）工程においてゲノムＤＮＡを抽出する方法としては、組織からゲノムＤＮＡを抽出する方法である限り特に制限されず、例えば、MagNA PureLC DNA Isolation Kit I（ロシュ・ダイアグノスティックス製）やAllprep DNA/RNA mini kit、及びQIAamp DNA FFPE Tissue kit（QIAGEN社製）などの市販品を添付のプロトコールにしたがって用いる方法の他、後述の実施例に記載されているような、フェノール・クロロホルム処理、及びエタノールあるいはイソプロパノール沈殿等を利用した常法（Molecular Cloning第３版Volume1のプロトコール参照）を例示することができる。
【００２０】
上記（Ｂ）工程及び（Ｄ）工程中において、本件ＣｐＧ部位のメチル化頻度を測定する方法としては、本件ＣｐＧ部位のメチル化を検出・定量し得る方法である限り特に制限されず、本件ＣｐＧ部位のポリヌクレオチド配列に特異的にハイブリダイズし得るポリヌクレオチド（プローブ）の標識物を用いたサザンハイブリダイゼーション法や、メチル化感受性制限酵素を用いて解析する方法や、亜硫酸水素塩（バイサルファイト）処理を利用したバイサルファイト法を例示することができる。上記ポリヌクレオチド（プローブ）のヌクレオチド数としては、本件ＣｐＧ部位にハイブリダイズし得る限り、特に制限されず、例えば７個以上を例示することができ、特異性の観点からは、１５個以上、より好ましくは２５個以上を好適に例示することができる。また、上記ポリヌクレオチド（プローブ）の標識物における標識物質としては、ペルオキシダーゼ（例えば、horseradish peroxidase）、アルカリフォスファターゼ、β−Ｄ−ガラクトシダーゼ、グルコースオキシダーゼ、グルコ−ス−６−ホスフェートデヒドロゲナーゼ、アルコール脱水素酵素、リンゴ酸脱水素酵素、ペニシリナーゼ、カタラーゼ、アポグルコースオキシダーゼ、ウレアーゼ、ルシフェラーゼ若しくはアセチルコリンエステラーゼ等の酵素、フルオレスセインイソチオシアネート、フィコビリタンパク、希土類金属キレート、ダンシルクロライド若しくはテトラメチルローダミンイソチオシアネート等の蛍光物質、緑色蛍光タンパク質（Green Fluorescence Protein；ＧＦＰ）、シアン蛍光タンパク質（Cyan Fluorescence Protein；ＣＦＰ）、青色蛍光タンパク質（Blue Fluorescence Protein；ＢＦＰ）、黄色蛍光タンパク質（Yellow Fluorescence Protein；ＹＦＰ）、赤色蛍光タンパク質（Red Fluorescence Protein；ＲＦＰ）、ルシフェラーゼ（luciferase）等の蛍光タンパク質、３Ｈ 、１４Ｃ、１２５Ｉ若しくは１３１Ｉ等の放射性同位体、ビオチン、アビジン、又は化学発光物質を挙げることができる。
【００２１】
上記メチル化感受性制限酵素を用いて解析する方法は、制限酵素の認識配列中のシトシンがメチル化されることで、ＤＮＡを切断できなくなる現象を利用したものである。必要な反応時間は、制限酵素の種類により異なるが、１〜数時間程度である。また、これまでは、制限酵素を用いていることから検出可能な配列が限られていたが、現在では、入手可能なメチル化感受性制限酵素は１００種類以上あり、認識配列も多彩であるため、本件ＣｐＧ部位のほとんどを標的とすることが可能である。メチル化感受性制限酵素により処理後、本件ＣｐＧ部位の切断効率は、サザンハイブリダイゼーションやＰＣＲ（polymerase chain reaction）等の方法で解析することができる。
【００２２】
また、上記バイサルファイト処理を利用したバイサルファイト法とは、バイサルファイト処理によって、本件ＣｐＧ部位における非メチル化シトシンのみをウラシルに変換するバイサルファイト反応を利用した方法であり、現在最も一般的に用いられている方法である。具体的には、以下の（i）〜（iv）の４つの方法等を例示することができる。
【００２３】
（i）バイサルファイトシーケンシング法
サンプルであるゲノムＤＮＡに対してバイサルファイト処理を行うことによって、ゲノムＤＮＡにおける、メチル化修飾を受けたシトシンはそのまま変換せず、非メチル化シトシンのみをウラシルに変換させる。このゲノムＤＮＡについてシーケンス反応を行うと、ウラシルはチミンとして表現されることとなる。バイサルファイト処理前後のゲノムＤＮＡの配列データを比較し、バイサルファイト処理前後のいずれもシトシンである部位はメチル化シトシンであることが分かり、バイサルファイト処理前にシトシンでありバイサルファイト処理後にチミンとなった部位は非メチル化シトシンであることが分かる。
【００２４】
（ii）メチル化特異的ＰＣＲ（ＭＳＰ、Methylation-Specific PCR)法
サンプルであるゲノムＤＮＡに対してバイサルファイト処理を行うと、メチル化ＣｐＧ部位と非メチル化ＣｐＧ部位とは、異なった塩基配列が生じることとなる。このことを利用して、各々において異なる塩基配列の部位に特異的なＰＣＲプライマーを設計して、ＰＣＲ産物の有無でメチル化状態を検出する方法である。このＭＳＰ法では、バイサルファイト処理をしたＤＮＡをそのまま解析に使うことができるため、短時間で結果を確認することが可能である。
【００２５】
（iii）ＣＯＢＲＡ（COmbined Bisulfite Restriction Analysis）法
バイサルファイト処理によりＤＮＡメチル化状態依存的に特定の制限酵素認識配列の塩基配列が変化することを利用し、バイサルファイト−ＰＣＲ後に酵素処理したＰＣＲ産物を電気泳動することでメチル化頻度を測定する方法であり、後述の実施例で用いている方法でもある。特定の認識配列におけるメチル化頻度を、ＰＣＲ産物の切断の有無の量比にて定量的に解析できる。
【００２６】
（iv）HumanMethylation450BeadChip（イルミナ株式会社製）を用いた方法
この方法は、バイサルファイトシーケンシング法の１種であるバイサルファイトパイロシークエンシング法を基盤とした、ゲノムワイドなＤＮＡメチル化（ＤＮＡメチローム）解析法である。サンプルから抽出したゲノムＤＮＡについてバイサルファイト処理を行った後、制限酵素によってＤＮＡを断片化し、該ＤＮＡ断片上の本件ＣｐＧ部位におけるメチル化頻度（メチル化レベル）を、付属のBeadChipを用いて検出する。このBeadChipには、本件ＣｐＧ部位を含む約４５万個の各ＣｐＧ部位におけるシトシンがウラシルに変換したＤＮＡ配列に特異的にハイブリダイズし得るオリゴヌクレオチドプローブ（非メチル検出用プローブ）の標識物と、変換していないＤＮＡ配列に特異的にハイブリダイズし得るオリゴヌクレオチドプローブ（メチル検出用プローブ）の標識物がそれぞれ固定されている。それぞれのプローブにハイブリダイズしたＤＮＡ断片の量を蛍光で簡便に検出することによって、ゲノムＤＮＡ上の本件ＣｐＧ部位におけるメチル化頻度を測定することができる。本件ＣｐＧ部位におけるメチル化頻度の好適な指標として、β値を挙げることができる。β値とは、例えば、標識物質が蛍光物質や蛍光タンパク質の場合、測定によって得られた、各ＣｐＧ部位に対応する非メチル検出用プローブの蛍光値（signal A）、及び、メチル検出用プローブの蛍光値（signal B）について、以下の計算式により算出される値である。
β＝（本件ＣｐＧ部位におけるsignal Bの最大値）／（本件ＣｐＧ部位におけるsignal Aの最大値＋本件ＣｐＧ部位におけるsignal Bの最大値＋１００）；
この計算式によると、各ＣｐＧ部位のメチル化頻度が、０（完全非メチル化）〜１（完全メチル化）の範囲で算出されることとなる。
【００２７】
上記（Ｅ）工程中における、ＤＮＡメチル化頻度に基づき多変量解析を行う方法としては、抽出した各組織のゲノムＤＮＡのメチル化状態の類似度を統計的に解析できる手法であれば特に制限されないが、階層的クラスタリング、自己組織化マッピングあるいは主成分分析等を好適に例示することができ、なかでも一般に広く認識されている階層的クラスタリングを用いることが好ましい。
【００２８】
上記工程（Ｂ）及び（Ｄ）において、本件ＣｐＧ部位におけるメチル化の頻度を測定する領域は特に制限されないが、好ましくは、
配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列、
配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列、
配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列、
配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列、
配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列、
配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列、
配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列、
配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列、
配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列、
配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列、
のメチル化頻度を測定することが好ましい。
【００２９】
上記（Ｅ）工程は、「工程（Ｂ）において測定した本件ＣｐＧ部位におけるメチル化の頻度を、ＤＮＡメチル化頻度に基づき多変量解析を行い、バイオマーカー遺伝子パネルを作製する工程」である。ここで、上記工程（Ａ）において用いた検体組織と、上記工程（Ｃ）において用いた検体組織とは同種の組織であることが好ましい。また、該検体組織は、被検体と同じ生物種のもの（被検体がヒトである場合にはヒトのもの）を用いることが好ましい。
【００３０】
このように、検体組織から抽出したゲノムＤＮＡ中の、本件ＣｐＧ部位におけるメチル化頻度に基づいて、該検体組織が子宮平滑筋腫、子宮平滑筋肉腫、及び子宮平滑筋層のいずれにあたるかを判定することができる。より具体的には、工程（Ｄ）において測定した被検体組織におけるメチル化頻度データを、基準となる工程（Ｂ）で測定したメチル化頻度データと共に、工程（Ｅ）による多変量解析に供することで、前記被検体が子宮平滑筋腫、子宮平滑筋肉腫、及び子宮平滑筋層のいずれにあたるかを区別することができる。
【００３１】
本発明の子宮平滑筋における腫瘍の診断用キット（以下、「本件診断用キット」ということがある）としては、本件バイオマーカー遺伝子領域におけるＣｐＧ配列のメチル化頻度を測定するためのプライマー若しくはプローブ、又はそれらの標識物を備えた、子宮平滑筋における腫瘍を診断するためのキットであれば特に制限されず、ここで本件診断用キットは、子宮平滑筋における腫瘍を診断するためのキットに関する用途発明であり、このキットには、一般にこの種の診断キットに用いられる成分、例えば担体、ｐＨ緩衝剤、安定剤の他、取扱説明書、子宮平滑筋における腫瘍を診断するための説明書等の添付文書が通常含まれる。
【００３２】
本件診断用キットにおけるプライマーとしては、本件ＣｐＧ部位のメチル化頻度を測定することができ、かつ、ＰＣＲにより増幅可能なプライマーであればよく、例えば、本件ＣｐＧ部位の上流のヌクレオチド配列の相補配列にアニールするフォワードプライマーと、本件ＣｐＧ部位の下流のヌクレオチド配列にアニールするリバースプライマーとのプライマーセットや、本件ＣｐＧ部位の上流のヌクレオチド配列の相補配列にアニールし、かつ本件ＣｐＧ部位のシトシンで３’末端がアニールするフォワードプライマーと、本件ＣｐＧ部位の下流のヌクレオチド配列にアニールするリバースプライマーとのプライマーセットや、本件ＣｐＧ部位の上流のヌクレオチド配列の相補配列にアニールするフォワードプライマーと、本件ＣｐＧ部位のシトシンで３’末端がアニールするリバースプライマーとのプライマーセットや、本件ＣｐＧ部位の上流のヌクレオチド配列の相補配列にアニールし、かつ本件ＣｐＧ部位のシトシンがウラシルに変換された場合にかかるウラシルの部位で３’末端がアニールするフォワードプライマーと、本件ＣｐＧ部位の下流のヌクレオチド配列にアニールするリバースプライマーとのプライマーセットや、本件ＣｐＧ部位の上流のヌクレオチド配列の相補配列にアニールするフォワードプライマーと、本件ＣｐＧ部位のシトシンがウラシルに変換された場合にかかるウラシルの部位で３’末端がアニールするリバースプライマーとのプライマーセットを挙げることができる。また、プライマー配列の長さ、アニーリングする部位等は、ＤＮＡの増幅効率や特異性を考慮して適宜選択することができる。例えば、プライマー配列の長さとしては、１０〜３０塩基を選択することができる。なお、ＮＣＢＩのデータベースに登録されている本件バイオマーカー遺伝子のヌクレオチド配列の５’末端側を上流とし、３’末端側を下流とする。
【００３３】
本件診断用キットにおけるプローブとしては、本件ＣｐＧ部位を含むヌクレオチド配列にハイブリダイゼーションするプローブや、本件ＣｐＧ部位のシトシンがウラシルに変換された上記ヌクレオチド配列にハイブリダイゼーションするプローブを挙げることができ、プローブの長さ、ハイブリダイズする部位等は、ハイブリダイゼーションの効率や特異性を考慮して適宜選択することができる。
【００３４】
本件診断用キットの標識物における標識物質としては、前述したように、ペルオキシダーゼ（例えば、horseradish peroxidase）、アルカリフォスファターゼ、β−Ｄ−ガラクトシダーゼ、グルコースオキシダーゼ、グルコ−ス−６−ホスフェートデヒドロゲナーゼ、アルコール脱水素酵素、リンゴ酸脱水素酵素、ペニシリナーゼ、カタラーゼ、アポグルコースオキシダーゼ、ウレアーゼ、ルシフェラーゼ若しくはアセチルコリンエステラーゼ等の酵素、フルオレスセインイソチオシアネート、フィコビリタンパク、希土類金属キレート、ダンシルクロライド若しくはテトラメチルローダミンイソチオシアネート等の蛍光物質、緑色蛍光タンパク質（Green Fluorescence Protein；ＧＦＰ）、シアン蛍光タンパク質（Cyan Fluorescence Protein；ＣＦＰ）、青色蛍光タンパク質（Blue Fluorescence Protein；ＢＦＰ）、黄色蛍光タンパク質（Yellow Fluorescence Protein；ＹＦＰ）、赤色蛍光タンパク質（Red Fluorescence Protein；ＲＦＰ）、ルシフェラーゼ（luciferase）等の蛍光タンパク質、３Ｈ 、１４Ｃ、１２５Ｉ若しくは１３１Ｉ等の放射性同位体、ビオチン、アビジン、又は化学発光物質を挙げることができる。
【００３５】
以下、実施例により本発明をより具体的に説明するが、本発明の技術的範囲はこれらの例示に限定されるものではない。
【実施例１】
【００３６】
［サンプル組織の採取及びＤＮＡの抽出］
山口大学医学部、及び医学部附属病院遺伝子解析研究に関する倫理審査委員会の承認のもと、医学的適応により子宮全摘術を受けた患者より、組織サンプルを得た。用いた検体は、正常筋層１８検体、子宮筋腫１８検体、及び子宮肉腫１２検体だった。また、これとは別にヒト平滑筋肉腫由来細胞株２種（ヒト子宮肉腫細胞株ＳＫＮ、及びヒト卵巣肉腫細胞株ＲＫＮ）（独立行政法人医薬基盤研究所 ＪＣＲＢ細胞バンクより購入）を用いた。
組織検体からゲノムＤＮＡを抽出した。組織片１０〜３０ｍｇを、液体窒素で凍結後、粉砕し粉末状にした。ＤＮＡ抽出は一般的な方法で行った。粉末状組織片をProteinase K（QIAGEN社製）で処理した後、フェノール／クロロフォルムを用いて抽出を行い、エタノール沈殿で精製した。
【００３７】
［筋腫特異的バイオマーカー候補遺伝子の選定］
出願人はこれまでに、子宮筋腫と正常筋層それぞれ各３検体のＤＮＡメチロームデータをもとに、子宮筋腫特異的にメチル化変異が生じている１２０遺伝子を抽出している（PLoS One, 8(6): e66632, 2013）。ＤＮＡメチル化パターンは組織・細胞種ごとに特異的である。子宮筋腫と子宮肉腫はともに子宮平滑筋に由来するが本質的には異なる腫瘍であるため、出願人は、子宮筋腫特異的にメチル化変異が生じている遺伝子群が子宮筋腫と子宮肉腫とを明確に区別するバイオマーカーとなり得る可能性は高いと考えた。そこで、前記１２０遺伝子のうち、ＤＮＡメチル化解析用プライマーの構築に成功した３３遺伝子を子宮筋腫特異的メチル化変異遺伝子の候補として選出した。３３遺伝子を表１に示す。
【００３８】
【表１】
0006583817.tif 000002
【００３９】
［ＤＮＡメチル化頻度の測定］
ＤＮＡメチル化頻度の測定には、ＣＯＢＲＡ法を用いた。ゲノムＤＮＡをバイサルファイト処理した後にバイサルファイト−ＰＣＲを行うと、非メチル化シトシンはチミンに変換され、他方、メチル化シトシンはシトシンのままで変換されないという反応が生じる。この反応により、特定の制限酵素認識配列の塩基配列がメチル化状態に依存して変化することを利用し、酵素処理したＰＣＲ産物を電気泳動することでメチル化頻度を解析する。
バイサルファイト処理は、EpiTect(R) Bisulfite kit（Qiagen社製）を用い、キット付属の説明書に基づいて行った。１μｇのゲノムＤＮＡに対し、バイサルファイト反応を、９５℃  ５分間、６５℃  ８５分間、９５℃  ５分間、及び６５℃  １７５分間の条件で行った。バイサルファイト−ＰＣＲは、１／１００量のバイサルファイトＤＮＡをテンプレートとし、Biotaq HS DNA polymerase（Bioline社製）を用いて、９５℃  １０分間、［９５℃  ０．５分間、６０℃  ０．５分間、７２℃  ０．５分間］×３５サイクル、７２℃  ７分間の反応条件で行った。制限酵素は配列TCGAを認識するTaqI、及び配列ACGTを認識するHpyCH4IVを使用した。制限酵素処理には半量のＰＣＲ産物を使用し、TaqIは３７℃、HpyCH4IVは６５℃でそれぞれ２時間以上反応した。制限酵素処理したＰＣＲ産物はアガロースゲル電気泳動で分離し、メチル化率を算出した。
【００４０】
まずは、表１に記載の遺伝子において、ＣＯＢＲＡ法による１遺伝子のメチル化解析に基づき、正常筋層と子宮筋腫とを判別可能かを検討したが、１遺伝子のメチル化頻度だけでは多数の症例を明瞭に区別することは難しいことが明らかとなった。１例として各１８症例の正常筋層検体、及び子宮筋腫検体を用い、ＡＬＸ１遺伝子のメチル化頻度を解析した結果を図１に示す。
【実施例２】
【００４１】
［筋腫特異的バイオマーカーの選出］
実施例１を踏まえ、筋腫特異的バイオマーカー遺伝子パネルを構築するにあたり、判別精度を保ちつつメチル化解析の労力を最小限にできる遺伝子数を検討した。
子宮筋腫特異的メチル化変異遺伝子の候補として選出した上記３３遺伝子について、３症例の子宮筋腫検体を用いて１遺伝子のメチル化頻度をＣＯＢＲＡ法により解析した。その結果、全ての筋腫検体で高メチル化変異或いは低メチル化変異が生じている、１７個の遺伝子を筋腫特異的バイオマーカー遺伝子として同定した。更に、子宮筋腫検体を１１症例に増やし同様の解析を行い、７割以上の筋腫検体で共通して高メチル化変異あるいは低メチル化変異が生じた１０個の遺伝子（ＡＬＸ１、ＣＢＬＮ１、ＣＯＲＩＮ、ＦＯＸＰ１、ＧＡＴＡ２、ＩＧＬＯＮ５、ＮＰＴＸ２、ＮＴＲＫ２、ＰＲＬ、及びＳＴＥＡＰ４）を特定した。
特定した筋腫特異的バイオマーカー遺伝子について、バイサルファイト−ＰＣＲに用いたＰＣＲプライマーを表２に示す。表２に示したプライマー配列は、バイサルファイト処理により非メチル化シトシンがウラシル（チミン）に変換された配列を増幅できるように設計したものである。また、本件バイオマーカー遺伝子のバイサルファイト−ＰＣＲにて増幅されたゲノム領域のヌクレオチド配列を、配列番号１〜１０に示す（配列番号１：ＡＬＸ１、配列番号２：ＣＢＬＮ１、配列番号３：ＣＯＲＩＮ、配列番号４：ＦＯＸＰ１、配列番号５：ＧＡＴＡ２、配列番号６：ＩＧＬＯＮ５、配列番号７：ＮＰＴＸ２、配列番号８：ＮＴＲＫ２、配列番号９：ＰＲＬ、配列番号１０：ＳＴＥＡＰ４）。ＣＯＢＲＡ法にてメチル化頻度を解析したＣｐＧ部位は、以下のとおりである。
配列番号１に示されるヌクレオチド配列の１４４−１４５番目のヌクレオチド配列
配列番号２に示されるヌクレオチド配列の１４６−１４７番目のヌクレオチド配列
配列番号３に示されるヌクレオチド配列の１８１−１８２番目のヌクレオチド配列
配列番号４に示されるヌクレオチド配列の６７−６８番目のヌクレオチド配列
配列番号５に示されるヌクレオチド配列の１７１−１７２番目のヌクレオチド配列
配列番号６に示されるヌクレオチド配列の３７７−３７８番目のヌクレオチド配列
配列番号７に示されるヌクレオチド配列の１５９−１６０番目のヌクレオチド配列
配列番号８に示されるヌクレオチド配列の６５−６６番目のヌクレオチド配列
配列番号９に示されるヌクレオチド配列の２５９−２６０番目のヌクレオチド配列
配列番号１０に示されるヌクレオチド配列の１０２−１０３番目のヌクレオチド配列
【００４２】
【表２】
0006583817.tif 000003
【００４３】
［ＤＮＡメチル化変異子宮筋腫特異的バイオマーカー遺伝子パネルの構築］
各症例のメチル化状態の類似度を測る尺度としてピアソン相関を用い、クラスタリングのアルゴリズムとして平均連結法を採用したクラスタリングを行った。実際の階層的クラスタリングには、オープンリソースソフトであるMultiExperiment Viewer（ＭｅＶ）を使用し、症例方向のクラスタリングを行った。
各１８症例の正常筋層検体、及び子宮筋腫検体における階層的クラスタリングの結果を図２に示す。選出した前述の１０遺伝子により、１８症例の正常筋層と子宮筋腫がクラスタリングにより完全に異なるクラスターに区別された。一方で、筋腫検体においては樹形図（デンドログラム）の約半分の高さ（点線）で切り取った場合、メチル化パターンが若干異なる３つのサブグループ（Ｌ１、Ｌ２およびＬ３）に分類されることも示された。
【実施例３】
【００４４】
［クラスタリングによる子宮肉腫検体の分別］
１２症例の子宮肉腫検体、及び肉腫細胞株４サンプル（ヒト子宮肉腫細胞株ＳＫＮ、及びヒト卵巣肉腫細胞株ＲＫＮをそれぞれ２サンプルずつ）を、実施例２にて構築した各１８症例の正常筋層検体、及び子宮筋腫検体における階層的クラスタリングの結果と共に示した結果を図３に示す。
８症例の肉腫検体と４株の肉腫細胞株は、筋腫・筋層とは異なるクラスターに分別された。一方、４症例の肉腫検体は筋腫のサブグループＬ２に分別された。筋腫を含まないクラスターに分類された検体のうち、筋層のクラスターに含まれなかったものは肉腫と判断できる。筋腫を含むクラスターに分類された検体のうち、サブグループＬ１およびＬ３にクラスタリングされたものは筋腫と判断でき、サブグループＬ２にクラスタリングされたものは４０％の確率で肉種の可能性を持つ検体と判断できる。
【実施例４】
【００４５】
［階層的クラスタリングによる子宮肉腫の診断］
追加した新規の症例データが肉腫、筋腫、及び筋層検体を含むどのクラスターに分類されるかを調べることで、その症例が肉腫であるか否かを診断することができる。新規な２症例（検体Ｘ、及びＹ）について、ＣＯＢＲＡ法により前述の１０遺伝子領域のメチル化データを取得し、そのデータを実施例３にて示した肉腫、筋腫および筋層検体のメチル化データに追加して階層的クラスタリングを行った。検体Ｘ、及びＹは病理診断により、それぞれ子宮筋腫、及び子宮肉腫と診断されている。階層的クラスタリングの結果を図４に示す。
検体Ｘ、及びＹは階層的クラスタリングにより、それぞれ子宮筋腫、及び子宮肉腫に分別された。これより、本願筋腫特異的バイオマーカー遺伝子パネルは、子宮肉腫の診断に応用可能であることが示された。
【産業上の利用可能性】
【００４６】
本発明によると、子宮平滑筋における腫瘍が良性腫瘍（子宮平滑筋腫）か、悪性腫瘍（子宮平滑筋肉腫）のいずれかであるかを精度よく診断できるため、子宮平滑筋における腫瘍の悪性度をより正確に分類し、効果的な治療が可能となり、ＱＯＬ（Quality of Life）の向上や医療費削減、及び妊娠可能な子宮温存による少子化問題の解消等の効果が期待される。
【図１】
0006583817.tif 000004
【図２】
0006583817.tif 000005
【図３】
0006583817.tif 000006
【図４】
0006583817.tif 000007
【配列表】
0006583817000001.app