(57)【特許請求の範囲】
【請求項１】
関節リウマチ患者から得た生体試料における、下記（１）〜（６）から選択されるいずれか１つの一塩基多型(SNP)の組み合わせを解析する工程と、
解析結果に基づいて、前記患者におけるメトトレキサート(MTX)の有効性を判定する工程とを含む、
関節リウマチ患者におけるメトトレキサートの有効性の診断を補助する方法：
（１）FPGS c.64A>G、FPGS c.*192A>G、GGH c.452C>T、GGH c.-401C>T、SLC19A1 c.696T>C、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A；
（２）FPGS c.64A>G、FPGS c.*192A>G、GGH c.452C>T、GGH c.-401C>T、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A；
（３）FPGS c.64A>G、FPGS c.*192A>G、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A；
（４）FPGS c.64A>G、FPGS g.130562725G>A、GGH c.452C>T、GGH c.-401C>T、SLC19A1 c.696T>C、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A；
（５）FPGS c.64A>G、FPGS g.130562725G>A、GGH c.452C>T、GGH c.-401C>T、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A；及び
（６）FPGS c.64A>G、FPGS g.130562725G>A、PPARG c.*1411G>C、FMO2 c.585A>G、EPHX1 c.357G>A、及びSLC28A3 c.267G>A。
【請求項２】
判定工程において、前記解析結果からMTXの有効性に関する数値指標を取得し、取得した数値指標と、所定の閾値とを比較した結果に基づいて、前記患者におけるMTXの有効性を判定する請求項１に記載の方法。
【請求項３】
MTXの有効性に関する数値指標が、多重ロジスティックモデルにより算出される予測値である請求項２に記載の方法。
【請求項４】
判定工程において、前記予測値が前記所定の閾値以上である場合、前記患者においてMTXの有効性が高いと判定し、前記予測値が前記所定の閾値より小さい場合、前記患者におけるMTXの有効性が低いと判定する請求項３に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、関節リウマチ患者におけるメトトレキサートの有効性の診断を補助する方法に関する。
【背景技術】
【０００２】
現在、関節リウマチの治療では、葉酸代謝拮抗薬として知られるメトトレキサート(MTX)が第一選択薬とされている。MTXの強力な免疫抑制作用によって関節破壊が阻止されるので、MTXは関節リウマチの症状のコントロールに大きな役割を果たしている。一般に、MTXによる治療効果は、その使用量に依存して増強するとされている。しかし、実際は、MTXの有効性には大きな個人差がある。個々の関節リウマチ患者について、治療効果の高い最適な量のMTXで治療することが望まれるが、患者におけるMTXの有効性を予測することは非常に難しい。
【０００３】
近年、MTXの代謝産物であるポリグルタミル化メトトレキサート(MTX-PG)の血中濃度が、MTXの有効性と関連していることが注目されている。例えば、非特許文献１及び非特許文献２には、MTXからMTX-PGへの代謝において、SLC19A1(RFC-1とも呼ばれる)、GGH及びFPGSのタンパク質分子が関与することが報告されている。具体的には、SLC19A1は、MTXの細胞内へのトランスポーターであり、GGHは、MTXをポリグルタミル化してMTX-PGに変換する酵素であり、FPGSは、MTX-PGを脱ポリグルタミル化してMTXに変換する酵素である。また、非特許文献３には、GGH遺伝子及びSLC19A1遺伝子の遺伝子多型がMTXの有効性の判定に利用できることが報告されている。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Davila L.及びRanganathan P., Nat. Rev. Rheumatol, vol.7, p.537-550 (2011)
【非特許文献２】Wessels J.A.M.ら, Arthritis Rheum, vol.56, p.1765-1775 (2007)
【非特許文献３】Hayashi H.ら, J. Clin. Pharmacy Therapeutics, vol.34, p.355-361 (2009)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
個々の患者についてMTXの有効性を判定できれば、その判定結果は、関節リウマチの治療方針を決定するための有用な指標となる。本発明者らは、上記の報告に鑑みて、MTXの代謝に関連するSLC19A1、GGH及びFPGSの各遺伝子のSNPを解析し、それらの解析結果に基づいてMTXの有効性を判定することを試みた。しかし、予想に反して、有効性について精度の高い判定結果を得ることができなかった。そのため、より高い精度でMTXの有効性を判定することを可能にするSNPマーカーセットのさらなる開発が求められる。
【課題を解決するための手段】
【０００６】
よって、本発明は、関節リウマチ患者から得た生体試料における、関節リウマチ患者から得た生体試料における、FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子の一塩基多型(SNP)の少なくとも１つと、FMO2遺伝子のSNPの少なくとも１つと、SLC28A3遺伝子のSNPの少なくとも１つと、EPHX1遺伝子のSNPの少なくとも１つと、PPARG遺伝子のSNPの少なくとも１つとを解析する工程と、解析結果に基づいて、上記の患者におけるMTXの有効性を判定する工程とを含む、関節リウマチ患者におけるMTXの有効性の診断を補助する方法を提供する。
【発明の効果】
【０００７】
本発明によれば、関節リウマチ患者について、MTXが有効であるか否かの診断を補助することを可能にする。
【図面の簡単な説明】
【０００８】
【図１】MTXの有効性の判定装置の一例を示した概略図である。
【図２】図１に示される判定装置のハードウェア構成を示すブロック図である。
【図３】図１に示される判定装置を用いたMTXの有効性の判定のフローチャートである。
【図４】判別式１を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図５】判別式２を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図６】判別式３を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図７】判別式４を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図８】判別式５を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図９】判別式６を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図１０】判別式７を用いて、MTXの有効性の判定を行ったときのROC曲線である。
【図１１】試薬キットの一例を示した概略図である。
【発明を実施するための形態】
【０００９】
本実施形態に係る関節リウマチ患者におけるメトトレキサートの有効性を判定する方法(以下、単に「方法」ともいう)では、まず、関節リウマチ患者から得た生体試料における、FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子のSNPの少なくとも１つと、FMO2遺伝子のSNPの少なくとも１つと、SLC28A3遺伝子のSNPの少なくとも１つと、EPHX1遺伝子のSNPの少なくとも１つと、PPARG遺伝子のSNPの少なくとも１つとを解析する。
【００１０】
SNPは、遺伝子多型の一種であり、ゲノムDNAの塩基配列が１カ所だけ変異した状態又はその部位を指す。SNPは、その変異が所定の母集団において１％以上の頻度で見られる点で、点突然変異とは異なるものである。本明細書では、SNPを「遺伝子多型部位における変異」ともいう場合がある。
【００１１】
関節リウマチ患者(以下、「RA患者」とも呼ぶ)は、関節リウマチであると診断された患者であれば特に限定されない。本実施形態に係る方法は、MTXの有効性の判定に関連することから、RA患者としては、MTXによる治療をまだ受けていない患者、及びMTXを投与されているが有効性がまだ判明していない患者が、本実施形態に係る方法の被験者として好ましい。したがって、本実施形態に係る方法は、「RA患者におけるMTXの有効性の予測する方法」と解釈することもできる。
【００１２】
本実施形態において、生体試料は、RA患者のゲノムDNAを含む生体由来の試料であれば特に制限されない。そのような生体試料としては、体液、尿、手術又は生検により採取した組織及び細胞などが挙げられる。体液としては、血液、血漿、血清、滑液、リンパ液、腹水、骨髄液、乳頭分泌液などが挙げられる。採取の簡便性の観点から、生体試料としては末梢血が好ましい。
【００１３】
本実施形態においては、生体試料からゲノムDNAを抽出することにより、測定用試料を調製することが好ましい。生体試料からのゲノムDNAの抽出方法は、当該技術において公知である。例えば、生体試料と、界面活性剤(例えばコール酸ナトリウム、ドデシル硫酸ナトリウムなど)を含む処理液とを混合し、得られた混合液に物理的処理(撹拌、ホモジナイズ、超音波破砕など)を施して、生体試料に含まれるゲノムDNAを該混合液中に遊離させることによって、ゲノムDNAを抽出することができる。この場合、混合液を遠心分離して細胞破片を沈殿させ、遊離したゲノムDNAを含む上清を後述の検出方法に用いることが好ましい。また、得られた上清を当該技術において公知の方法により精製してもよい。なお、生体試料からのゲノムDNAの抽出及び精製は、市販のキットを用いて行うこともできる。
【００１４】
本実施形態においては、検出対象の遺伝子多型部位が存在する領域を含む限り、RA患者のゲノムDNAから合成されるcDNA又はcRNAを測定用試料として用いてもよい。なお、ゲノムDNAからcDNA又はcRNAを合成する方法自体は当該技術において公知である。以下、RA患者の生体試料から得られるゲノムDNA、ならびにこのゲノムDNAから得られるcDNA及びcRNAを総称して、「RA患者の生体試料由来の核酸」ともいう。
【００１５】
上記の各遺伝子の名称を、以下の表１に示す。なお、SLC19A1遺伝子は、RFC-1(reduced folate carrier-1)とも呼ばれる。これらの遺伝子については、解析対象のSNPを含むゲノムDNAの塩基配列自体は公知である。これらの塩基配列は、例えば米国国立医学図書館の国立生物情報センター(National Center for Biotechnology Information：NCBI)により提供されるデータベース(http://www.ncbi.nlm.nih.gov/)などの公知のデータベースから知ることができる。
【００１６】
【表１】
0006542037.tif 000002
【００１７】
本実施形態において、解析対象となるSNPを有する遺伝子の組み合わせは、FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子と、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子の４遺伝子とを含む、少なくとも５つの遺伝子の組み合わせであれば特に限定されない。FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子は、特に限定されないが、FPGS遺伝子を含むことが好ましい。
【００１８】
本実施形態では、上記の各遺伝子について解析すべきSNPの位置及び数は特に限定されず、上記の各遺伝子について公知のSNPから適宜選択できる。好ましくは、上記の各遺伝子について、次に示される部位のSNPを解析する。FPGS遺伝子については、c.64位、c.*192位及びg.130562725位から選択される少なくとも１つの部位、好ましくはc.64位である。GGH遺伝子については、c.452位、c.-401位及びc.16位から選択される少なくとも１つの部位である。SLC19A1遺伝子については、c.696位及びc.80位である。FMO2遺伝子については、c.585位である。SLC28A3遺伝子については、c.267位である。EPHX1遺伝子については、c.357位である。PPARG遺伝子については、c.*1411位である。
【００１９】
上記の各遺伝子におけるSNPが存在する部位の表記について、以下に説明する。なお、この表記法は、Human Genome Variation Society (http://www.hgvs.org/)により定められたガイドラインに従っている。
【００２０】
FPGS遺伝子のc.64位は、FPGS遺伝子のコーディング配列の5'末端から下流側へ数えて64番目の塩基の位置を指す。ここで、「遺伝子のコーディング配列の5'末端から下流側へ数える」とは、ゲノムDNA上の遺伝子の塩基配列において、該遺伝子のコーディング配列の5'末端にある開始コドン(ATG)のアデニンを１番目として、下流側(3'側)へ塩基を数えていくことを意図する。
【００２１】
FPGS c.*192位は、FPGS遺伝子のコーディング配列の3'末端から下流側へ数えて192番目の塩基の位置を指す。ここで、「遺伝子のコーディング配列の3'末端から下流側へ数える」とは、ゲノムDNA上の遺伝子の塩基配列において、該遺伝子のコーディング配列の3'末端の塩基の3'側に隣接する塩基を１番目として、下流側(3'側)へ塩基を数えていくことを意図する。
【００２２】
FPGS g.130562725位は、ゲノムDNA上のFPGS遺伝子の塩基配列の１番目の塩基から下流側へ数えて130562725番目の塩基の位置を指す。
【００２３】
GGH c.452位は、GGH遺伝子のコーディング配列の5'末端から下流側へ数えて452番目の塩基の位置を指す。
【００２４】
GGH c.-401位は、GGH遺伝子のコーディング配列の5'末端から上流側へ数えて401番目の塩基の位置である。ここで、「遺伝子のコーディング配列の5'末端から上流側へ数える」とは、ゲノムDNA上の遺伝子の塩基配列において、該遺伝子のコーディング配列の5'末端にある開始コドン(ATG)のアデニンの5'側に隣接する塩基を１番目として、上流側(5'側)へ塩基を数えていくことを意図する。
【００２５】
GGH c.16位は、GGH遺伝子のコーディング配列の5'末端から下流側へ数えて16番目の塩基の位置を指す。
【００２６】
SLC19A1 c.696位は、SLC19A1遺伝子のコーディング配列の5'末端から下流側へ数えて696番目の塩基の位置を指す。
【００２７】
SLC19A1 c.80位は、SLC19A1遺伝子のコーディング配列の5'末端から下流側へ数えて80番目の塩基の位置を指す。
【００２８】
FMO2遺伝子のc.585位は、FMO2遺伝子のコーディング配列の5'末端から下流側へ数えて585番目の塩基の位置を指す。
【００２９】
SLC28A3遺伝子のc.267位は、SLC28A3遺伝子のコーディング配列の5'末端から下流側へ数えて267番目の塩基の位置を指す。
【００３０】
EPHX1遺伝子のc.357位は、EPHX1遺伝子のコーディング配列の5'末端から下流側へ数えて357番目の塩基の位置を指す。
【００３１】
PPARG遺伝子のc.*1411位とは、PPARG遺伝子のコーディング配列の3'末端から下流側へ数えて1411番目の塩基の位置を指す。
【００３２】
本実施形態においては、FPGS遺伝子のSNPは、FPGS c.64A>G、FPGS c.*192A>G及びFPGS g.130562725G>Aから選択される少なくとも１つが好ましい。特に、FPGS c.64A>Gを解析することが好ましい。また、FPGS g.130562725G>AとFPGS c.*192A>Gとの間には多重共線性が認められるので、これらのうち、いずれか一方のみを解析することがより好ましい。FPGS g.130562725G>A及びFPGS c.*192A>Gのうち、いずれを選択するかは特に限定されない。
【００３３】
本実施形態においては、GGH遺伝子のSNPは、GGH c.452C>T、GGH c.-401C>T及びGGH c.16T>Cから選択される少なくとも１つが好ましい。GGH c.16T>CとGGH c.-401C>Tとの間には多重共線性が認められるので、これらのうち、いずれか一方のみを解析することがより好ましい。GGH c.16T>C及びGGH c.-401C>Tのうち、いずれを選択するかは特に限定されない。
【００３４】
本実施形態においては、SLC19A1遺伝子のSNPは、SLC19A1 c.696T>C及びSLC19A1 c.80G>Aから選択される少なくとも１つが好ましい。SLC19A1 c.696T>CとSLC19A1 c.80G>Aとの間には多重共線性が認められるので、これらのうち、いずれか一方のみを解析することがより好ましい。SLC19A1 c.696T>C及びSLC19A1 c.80G>Aのうち、いずれを選択するかは特に限定されない。
【００３５】
本実施形態においては、FMO2遺伝子のSNPは、FMO2 c.585A>Gが好ましい。SLC28A3遺伝子のSNPは、SLC28A3 c.267G>Aが好ましい。EPHX1遺伝子のSNPは、EPHX1 c.357G>Aが好ましい。PPARG遺伝子のSNPは、PPARG c.*1411G>Cが好ましい。
【００３６】
上記のSNPの表記では、「>」の左側には変異前の塩基が示され、右側には変異後の塩基が示される。例えば、「FPGS c.64A>G」は、FPGS遺伝子のc.64位におけるアデニンからグアニンへの変異を示す。
【００３７】
上記のSNPのゲノムDNA上の位置は、NCBIにより提供されるデータベース(http://www.ncbi.nlm.nih.gov/)などの公知のデータベースから知ることができる。各SNPのRS番号（Reference SNP ID number）を表２に示す。表２中の配列番号１〜１２は、各SNPの前の500塩基及び後の500塩基を含む計1001塩基の配列を表している。すなわち、配列番号１〜１２で表される塩基配列において、501番目の塩基の位置がSNPの位置である。なお、配列番号１〜１２で表される塩基配列では、SNPの位置にある塩基(501番目の塩基)は、変異前の塩基を示している。
【００３８】
【表２】
0006542037.tif 000003
【００３９】
本実施形態において、SNPの解析は、SNPの有無、すなわち、遺伝子多型部位において変異が存在するか否かを決定することであってもよい。また、SNPの解析は、対立遺伝子のSNPの数を取得することであってもよいし、SNPの有無と相関する値を取得してもよい。SNPの有無と相関する値としては、例えば、後述のSNPの解析手段により取得される測定値そのものであってもよいし、該測定値に基づいて取得される値であってもよい。測定値に基づいて取得される値は、特に限定されず、例えば、和、差、積、比、対数、統計学的代表値(例えば平均値、中央値など)などが挙げられる。また、測定値に基づいて取得される値は、所定の式などによって、上記の測定値から算出される値であってもよい。
【００４０】
本実施形態では、SNPの解析手段は特に限定されず、当該技術において公知のSNP検出方法及びSNPタイピング方法から適宜選択すればよい。そのような方法としては、例えば、ダイレクトシーケンス法、クエンチング・プローブ(QProbe)法(特許第4950414号参照)、マイクロアレイ法などが挙げられる。また、TaqMan(商標) SNP Genotyping Assaysのような市販のSNPタイピング用キットを用いてもよい。
【００４１】
ダイレクトシーケンス法では、ゲノムDNA中のSNPを含む領域を、所定のプライマーセットを用いるPCR法によって増幅し、得られた増幅産物の塩基配列を解析することによって、SNPの存否を検出する。表１に示される遺伝子はいずれも、解析対象のSNPを含むゲノムDNAの塩基配列自体は公知であるので、その塩基配列に基づいて増幅用プライマー及びシーケンス用プライマーを設計すればよい。なお、プライマーは、当該技術において公知の核酸合成法により合成することができる。
【００４２】
QProbe法でSNPを検出する場合は、ゲノムDNA中のSNPを含む領域を、所定のプライマーセットを用いるPCR法によって増幅し、得られた増幅産物の塩基配列に相補的な配列を有するQProbeとハイブリダイズさせる。ここで、増幅産物にSNPが存在する場合、該増幅産物とQProbeとが解離する温度が、増幅産物にSNPが存在しない場合の温度と異なるので、QProbe法では、この差を利用してSNPを検出する。なお、QProbe自体は当該技術において公知であり、一般に入手可能である。
【００４３】
マイクロアレイ法でSNPを検出する場合、解析用マイクロアレイは、上記のSNPが存在する領域の塩基配列に相補的な核酸プローブを、基板上に固定して作製できる。なお、このようなマイクロアレイは、当該技術において公知の方法により作製できる。マイクロアレイによる解析では、生体試料由来の核酸は、当該技術において公知の標識物質により標識されていることが好ましい。なお、標識物質としては、蛍光物質、ビオチンなどのハプテン、放射性物質などが挙げられる。また、蛍光物質としては、Cy3、Cy5、FITC、Alexa Fluor(商標)などが挙げられる。このようにDNAを標識することにより、マイクロアレイ上のプローブからのシグナルの測定が容易になる。なお、DNAをこれらの標識物質で標識する方法は、当該技術において公知である。
【００４４】
上記のシグナルは、マイクロアレイの種類に応じて適切なシグナルであり得る。例えば、シグナルは、マイクロアレイの各プローブとハイブリダイズしたDNA断片が存在する場合に発生する電気的シグナルであってもよいし、生体試料由来の核酸が標識されている場合は、標識物質から生じる蛍光、発光などのシグナルであってもよい。シグナルの検出は、通常のマイクロアレイ測定装置に備えられたスキャナーにより行うことができる。スキャナーとしては、例えば、GeneChip(登録商標) Scanner3000 7G(Affymetrix社)などが挙げられる。
【００４５】
本実施形態に係る方法では、上記の解析工程で得られた解析結果に基づいて、RA患者におけるMTXの有効性を判定する。
【００４６】
本実施形態においては、解析結果に基づいて、RA患者について、MTXの有効性が高いか否かを判定してもよい。ここで、MTXの有効性の判定は、RA患者に対するMTXの治療効果の予測、又は、RA患者におけるMTXでの治療反応性の予測と解釈することもできる。例えば、MTXの有効性が高いと判定された場合、患者はMTXによる治療効果が高いと予測できる。一方、MTXの有効性が低いと判定された場合、患者はMTXによる治療効果が低いと予測できる。
【００４７】
本実施形態では、解析結果からMTXの有効性に関する数値指標を取得し、取得した数値指標と、所定の閾値とを比較した結果に基づいて、RA患者におけるMTXの有効性を判定することが好ましい。MTXの有効性に関する数値指標としては、例えば、上記の各遺伝子についてのSNPの有無、SNPの数、又はSNPの有無と相関する値を多変量解析して取得される数値指標が挙げられる。そのような多変量解析により取得される数値指標としては、例えば、多重ロジスティックモデルにより算出される予測値が挙げられる。なお、多変量解析自体は、当該技術において公知であり、市販のソフトウェアを用いてコンピュータに実行させることができる。
【００４８】
判定は、多重ロジスティックモデルに基づいて構築された回帰式を用いて行うことができる。そのような回帰式としては、例えば、下記の式(I)の回帰式を用いることができる。この式(I)の回帰式は、連続尺度による回帰式である。判定は、解析した各SNPに関する情報を式(I)に代入して得られるＰ値と、所定の閾値とを比較することにより行われる。ここで、SNPに関する情報として、遺伝子座において対立遺伝子の遺伝子多型部位のいずれにも変異がない場合は「０」を、対立遺伝子の遺伝子多型部位のいずれか一方に変異がある場合は「１」を、対立遺伝子の遺伝子多型部位の両方に変異がある場合は「２」を用いることができる。これらの値を、下記の式(I)の*SNPnに代入する。b0は切片の係数、b1〜bnは各SNPの係数であり、判別に使用するSNPの種類によって適宜設定できる。
【００４９】
【数１】
0006542037.tif 000004
【００５０】
あるいは、回帰式としては、下記の式(II)の回帰式を用いてもよい。この式(II)の回帰式は、名義尺度による回帰式である。判定は、解析した各SNPに関する情報を用いて式(II)から得られるＰ値と、所定の閾値とを比較することにより行われる。式(II)では、SNPに関する情報は、遺伝子座において対立遺伝子の遺伝子多型部位のいずれにも変異がない場合は「WT」であり、対立遺伝子の遺伝子多型部位のいずれか一方に変異がある場合は「HT」であり、対立遺伝子の遺伝子多型部位の両方に変異がある場合は「MT」である。式(II)においては、SNPに関する情報、すなわち、判別に使用するSNPの種類及びその遺伝子型(WT、HT又はMT)によって、各SNPの係数であるb1〜bnに代入する値が異なる。b1〜bnに代入する値は、当該技術において公知の方法により算出できる。式(II)において、b0は切片の係数である。
【００５１】
【数２】
0006542037.tif 000005
【００５２】
上記の式(I)又は(II)により算出されるＰ値は、MTXの有効性の予測値(以下、「有効性予測値」ともいう)を示す。本実施形態においては、所定の閾値を例えば0.5と設定したとき、有効性予測値が0.5以上であれば、RA患者においてMTXの有効性が高い(又は、MTXによる治療効果が高い)と判定できる。また、有効性予測値が0.5より小さければ、RA患者においてMTXの有効性が低い(又は、MTXによる治療効果が低い)と判定できる。なお、所定の閾値は、当業者が任意に設定できる。また、上記の回帰式(I)及び(II)を、本実施形態に係る方法における有効性予測値の判別式としてもよい。
【００５３】
本発明には、RA患者におけるMTXの有効性の診断を補助する方法も含まれる。この方法によれば、MTXの有効性についての判定結果が医師等に提供され、医師等によるMTXの有効性の診断を補助することができる。すなわち、医師等は、MTXの有効性の判定結果に基づいて、RA患者についてMTXの有効性を診断し、MTXを投与するか否かを判断できる。また、被験者がMTXをすでに投与されているRA患者の場合は、医師等は、MTXの投与を継続するか否か、又はMTXの投与量を変更するか否かを判断できる。
【００５４】
本発明には、RA患者におけるMTXの有効性の判定装置も含まれる。そのような装置としては、例えば、次のとおりである。
【００５５】
プロセッサ及び該プロセッサの制御下にあるメモリを含むコンピュータを備え、
該メモリには、下記のステップ：
関節リウマチ患者から得た生体試料における、FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子の一塩基多型(SNP)の少なくとも１つと、FMO2遺伝子のSNPの少なくとも１つと、SLC28A3遺伝子のSNPの少なくとも１つと、EPHX1遺伝子のSNPの少なくとも１つと、PPARG遺伝子のSNPの少なくとも１つとを解析するステップと、
取得した解析結果に基づいて、関節リウマチ患者におけるメトトレキサート(MTX)の有効性が高いか否かを判定するステップと
を該コンピュータに実行させるためのコンピュータプログラムが記録されている、
関節リウマチ患者におけるMTXの有効性の判定装置。
【００５６】
また、本発明には、RA患者におけるMTXの有効性の判定をコンピュータに実行させるためのコンピュータプログラムも含まれる。そのようなコンピュータプログラムとしては、例えば、次のとおりである。
【００５７】
コンピュータに読み取り可能な媒体に記録されているコンピュータプログラムであって、
下記のステップ：
関節リウマチ患者から得た生体試料における、FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子から選択される少なくとも１つの遺伝子の一塩基多型(SNP)の少なくとも１つと、FMO2遺伝子のSNPの少なくとも１つと、SLC28A3遺伝子のSNPの少なくとも１つと、EPHX1遺伝子のSNPの少なくとも１つと、PPARG遺伝子のSNPの少なくとも１つとを解析するステップと、
取得した解析結果に基づいて、関節リウマチ患者におけるメトトレキサート(MTX)の有効性が高いか否かを判定するステップと
を該コンピュータに実行させる、
RA患者におけるMTXの有効性の判定をコンピュータに実行させるためのコンピュータプログラム。
【００５８】
上記の媒体は、上記のコンピュータプログラムが非一時的に記録され、且つコンピュータに読取可能な媒体であってもよい。
【００５９】
以下に、本発明の方法を実施するのに好適な装置の一例を、図面を参照して説明する。しかし、本実施形態は、この例のみに限定されるものではない。図１は、MTXの有効性の判定装置の一例を示した概略図である。図１に示された判定装置１０は、測定装置２０と、該測定装置２０と接続されたコンピュータシステム３０とを含んでいる。
【００６０】
本実施形態において、測定装置２０は、マイクロアレイ上のプローブと結合した核酸に基づくシグナルを検出するマイクロアレイスキャナーである。本実施形態において、シグナルは、蛍光シグナルなどの光学的情報である。測定用試料と接触させたマイクロアレイを測定装置２０にセットすると、測定装置２０は、マイクロアレイ上のプローブに結合した、RA患者の生体試料由来の核酸に基づく光学的情報を取得し、得られた光学的情報をコンピュータシステム３０に送信する。
【００６１】
マイクロアレイスキャナーは、マイクロアレイ上のプローブに結合した核酸に基づくシグナルの検出が可能であれば特に限定されない。シグナルは、RA患者の生体試料由来の核酸の標識に用いられた標識物質によって異なることから、マイクロアレイスキャナーは、標識物質の種類に応じて適宜選択することができる。例えば、標識物質が放射性物質である場合、測定装置２０として、該放射性物質から生じる放射線を検出可能なマイクロアレイスキャナーが用いられる。
【００６２】
なお、SNPをダイレクトシーケンス法により検出する場合、測定装置２０は、核酸増幅装置及びシーケンス解析装置からなる装置であってもよい。この場合、測定用試料、核酸増幅用の酵素及びプライマーなどを含む反応液を測定装置２０にセットし、核酸増幅法によって反応液中の核酸を増幅させる。そして、測定装置２０は、増幅産物の塩基配列を解析して配列情報を取得し、得られた配列情報をコンピュータシステム３０に送信する。
【００６３】
コンピュータシステム３０は、コンピュータ本体３００と、入力部３０１と、検体情報や判定結果などを表示する表示部３０２とを含む。コンピュータシステム３０は、測定装置２０から光学的情報を受信する。そして、コンピュータシステム３０のプロセッサは、光学的情報に基づいて、RA患者におけるMTXの有効性を判定するプログラムを実行する。なお、コンピュータシステム３０は、図１に示されるように測定装置２０とは別個の機器であってもよいし、測定装置２０を内包する機器であってもよい。後者の場合、コンピュータシステム３０は、それ自体で判定装置１０となってもよい。
【００６４】
コンピュータ本体３００は、図２に示されるように、ＣＰＵ(Central Processing Unit)３１０と、ＲＯＭ(Read Only Memory)３１１と、ＲＡＭ(Random Access Memory)３１２と、ハードディスク３１３と、入出力インターフェイス３１４と、読出装置３１５と、通信インターフェイス３１６と、画像出力インターフェイス３１７とを備えている。ＣＰＵ３１０、ＲＯＭ３１１、ＲＡＭ３１２、ハードディスク３１３、入出力インターフェイス３１４、読出装置３１５、通信インターフェイス３１６及び画像出力インターフェイス３１７は、バス３１８によってデータ通信可能に接続されている。また、測定装置２０は、通信インターフェイス３１６により、コンピュータシステム３０と通信可能に接続されている。
【００６５】
ＣＰＵ３１０は、ＲＯＭ３１１に記憶されているプログラム及びＲＡＭ３１２にロードされたプログラムを実行することが可能である。ＣＰＵ３１０は、有効性予測値を算出し、ＲＯＭ３１１に格納されている判別式を読み出し、有効性を判定する。ＣＰＵ３１０は、判定結果を出力して表示部３０２に表示させる。
【００６６】
ＲＯＭ３１１は、マスクＲＯＭ、ＰＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭなどによって構成されている。ＲＯＭ３１１には、前述のようにＣＰＵ３１０によって実行されるプログラム及びこれに用いるデータが記録されている。ＲＯＭ３１１には、所定の閾値、多重ロジスティックモデルの判別式が記録されていてもよい。
【００６７】
ＲＡＭ３１２は、ＳＲＡＭ、ＤＲＡＭなどによって構成されている。ＲＡＭ３１２は、ＲＯＭ３１１及びハードディスク３１３に記録されているプログラムの読み出しに用いられる。ＲＡＭ３１２はまた、これらのプログラムを実行するときに、ＣＰＵ３１０の作業領域として利用される。
【００６８】
ハードディスク３１３は、ＣＰＵ３１０に実行させるためのオペレーティングシステム、アプリケーションプログラム(MTXの有効性の判定のためのコンピュータプログラム)などのコンピュータプログラム及び当該コンピュータプログラムの実行に用いるデータがインストールされている。ハードディスク３１３には、所定の閾値、多重ロジスティックモデルの判別式が記録されていてもよい。
【００６９】
読出装置３１５は、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤ−ＲＯＭドライブなどによって構成されている。読出装置３１５は、可搬型記録媒体４０に記録されたプログラム又はデータを読み出すことができる。
【００７０】
入出力インターフェイス３１４は、例えば、ＵＳＢ、ＩＥＥＥ１３９４、ＲＳ−２３２Ｃなどのシリアルインターフェイスと、ＳＣＳＩ、ＩＤＥ、ＩＥＥＥ１２８４などのパラレルインターフェイスと、Ｄ／Ａ変換器、Ａ／Ｄ変換器などからなるアナログインターフェイスとから構成されている。入出力インターフェイス３１４には、キーボード、マウスなどの入力部３０１が接続されている。操作者は、当該入力部３０１により、コンピュータ本体３００に各種の指令を入力することが可能である。
【００７１】
通信インターフェイス３１６は、例えば、Ethernet(登録商標)インターフェイスなどである。コンピュータ本体３００は、通信インターフェイス３１６により、プリンタなどへの印刷データの送信も可能である。
【００７２】
画像出力インターフェイス３１７は、ＬＣＤ、ＣＲＴなどで構成される表示部３０２に接続されている。これにより、表示部３０２は、ＣＰＵ３１０から与えられた画像データに応じた映像信号を出力できる。表示部３０２は、入力された映像信号にしたがって画像(画面)を表示する。
【００７３】
次に、判定装置１０による、MTXの有効性の高低を判定する処理手順を説明する。ここでは、マイクロアレイ上のプローブに結合した、RA患者の生体試料由来の核酸に基づく蛍光情報から、多重ロジスティックモデルにより有効性予測値を取得し、得られた有効性予測値を用いて判定を行なう場合を例として説明する。しかし、本実施形態は、この例のみに限定されるものではない。
【００７４】
図３を参照して、ステップＳ１０１において、判定装置１０のＣＰＵ３１０は、測定装置２０から蛍光情報を取得する。次に、ステップＳ１０２において、ＣＰＵ３１０は、取得した蛍光情報から蛍光強度を算出し、ＲＡＭ３１２に記憶する。そして、ステップＳ１０３において、ＣＰＵ３１０は、ＲＡＭ３１２に記憶された前記蛍光強度から各遺伝子多型部位における変異の有無を決定し、ＲＯＭ３１１又はハードディスク３１３に記憶された判別式にしたがって、多重ロジスティックモデルに基づく有効性予測値Ｐを算出する。
【００７５】
その後、ステップＳ１０４において、ＣＰＵ３１０は、算出された有効性予測値と、ＲＯＭ３１１又はハードディスク３１３に記憶された所定の閾値とを用いて、RA患者におけるMTXの有効性の高低を判定する。ここで、有効性予測値Ｐが所定の閾値よりも小さいとき、処理は、ステップＳ１０５に進行し、ＣＰＵ３１０は、RA患者におけるMTXの有効性は低いことを示す判定結果をＲＡＭ３１２に記憶する。一方、有効性予測値Ｐが所定の閾値よりも低くないとき(すなわち、有効性予測値Ｐが閾値以上であるとき)、処理は、ステップＳ１０６に進行し、ＣＰＵ３１０は、RA患者におけるMTXの有効性が高いことを示す判定結果をＲＡＭ３１２に記憶する。
【００７６】
そして、ステップＳ１０７において、ＣＰＵ３１０は、判定結果を出力し、表示部３０２に表示させたり、プリンタに印刷させたりする。これにより、RA患者においてMTXの有効性が高いか否かの判定を補助する情報を医師などに提供することができる。
【００７７】
上述のSNPを検出するための各種試薬は、図１１に示されるような試薬キットとしてユーザに提供され得る。この試薬キットは、上述のSNPの検出方法に応じた試薬を含む。例えば、この試薬キットは、プライマーやプローブなどのSNP検出用オリゴヌクレオチドを収容した試薬容器と、蛍光物質などの標識物質を収容した試薬容器と、DNAポリメラーゼなどの核酸増幅酵素を収容した試薬容器とを含む。さらに、検出方法に応じて、マイクロアレイやTaqMan(商標) Probeなどの試薬を含んでいてもよい。また、この試薬キットは、試薬の添付文書を同梱していてもよい。この添付文書には、たとえば、試薬キットの構成や検出プロトコールなどが記載される。図１１において、５０は、試薬キットを示し、５１〜５３は、試薬容器を示し、５４は、添付文書を示し、５５は、梱包箱を示す。
【００７８】
以下に、本発明を実施例により詳細に説明するが、本発明はこれらの実施例に限定されるものではない。
【実施例】
【００７９】
比較例１
MTXの代謝に関連するFPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子のSNPを解析し、その解析結果に基づいて、RA患者におけるMTXの有効性を判定することを試みた。
【００８０】
(１)被験者
まず、関節リウマチと診断され、MTXで治療中の患者の診療情報を調査し、以下の２群の患者(合計100名)を選出した。
1. MTXの投与量が８mg/week以下であって、且つ、血液検査によるALT値又はAST値が40 IU/L以上を示した患者(44名)。
2. MTXの投与量が８mg/week以上であって、且つ、血液検査によるALT値又はAST値が40 IU/L以下の患者(56名)。
【００８１】
上記の100名の患者のうち、MTX投与後に、CRP(Ｃ反応性タンパク)の数値に基づくDAS(disease activity score)28(以下、「DAS28(CRP)」ともいう)及び肝機能を経時的に測定した49名を、被験者として選出した。これらの患者について、MTX投与の24週後に、DAS28(CRP)により疾患活動性を評価し、EULAR(欧州リウマチ学会)治療反応性基準によりMTXの有効性を判定した。その結果、49名のうち、有効群(good)が22名であり、無効群(moderate及びnone)が27名であった。
【００８２】
(２)生体試料
生体試料として、上記の被験者(49名)の末梢血を採取した。QIAamp DNA blood mini kit(QIAGEN社)又はQuickGene-810(Kurabo社)を用いて、末梢血(200μL)から血液細胞のゲノムDNAを抽出した。抽出したゲノムDNAは、後述の遺伝子多型解析に使用するまで、−30℃にて保管した。
【００８３】
(３)遺伝子多型解析
各被験者のゲノムDNAについて、上記の３遺伝子に存在する８つのSNPを、TaqMan(商標) SNP Genotyping Assays(ライフテクノロジーズ社)を用いて解析した。上記の８つのSNP及びそれらの解析に用いたキットのAssay IDを、表３に示す。
【００８４】
【表３】
0006542037.tif 000006
【００８５】
なお、FPGS c.64A>G及びSLC19A1 c.80G>Aの解析には、ライフテクノロジーズ社に委託して製造したTaqMan(商標) SNP Genotyping Assaysのカスタム製品を用いた。この製品に含まれるプライマー及びTaqMan Probeの塩基配列を、以下に示す。
【００８６】
FPGS c.64A>G解析用プライマー及びプローブ
・フォワードプライマー：5'- CGCGCCGCTCTATTCCT -3'  (配列番号１３)
・リバースプライマー：5'- GGCCACGCGCTCAAG -3'  (配列番号１４)
・TaqMan Probe Reporter 1：5'- CGCGGCATAACGA -3'  (配列番号１５)
・TaqMan Probe Reporter 2：5'- CGCGGCGTAACGA -3'  (配列番号１６)
【００８７】
SLC19A1 c.80G>A解析用プライマー及びプローブ
・フォワードプライマー：5'- GCCTGACCCCGAGCTC -3'  (配列番号１７)
・リバースプライマー：5'- CATGAAGCCGTAGAAGCAAAGGTA -3'  (配列番号１８)
・TaqMan Probe Reporter 1：5'- ACACGAGGTGCCGCC -3'  (配列番号１９)
・TaqMan Probe Reporter 2：5'- ACACGAGGCGCCGCC -3'  (配列番号２０)
【００８８】
各被験者のゲノムDNA(10 ng)について、上記のTaqMan(商標) SNP Genotyping AssaysとTaqMan(商標) GTXpress Master Mix(ライフテクノロジーズ社)を用いて、Real-Time PCR法にて測定を行った。なお、具体的な操作は、ライフテクノロジーズ社が提供するプロトコールに従って行った。
【００８９】
49名の被験者についての遺伝子多型解析の結果、表３に示す８つのSNPのうち、GGH c.16T>CとGGH c.-401C>Tとの間、及びSLC19A1 c.80G>AとSLC19A1 c.696T>Cとの間、FPGS g.130562725G>AとFPGS c.*192A>Gとの間には、多重共線性を認めた。表４に、上記の被験者における各組み合わせの挙動を示す。また、表５に、これらのSNPの組み合わせについてのSpearmanの相関係数を示す。
【００９０】
【表４】
0006542037.tif 000007
【００９１】
【表５】
0006542037.tif 000008
【００９２】
表４に示されるように、GGH c.16T>CとGGH c.-401C>Tは、上記の被験者において全く同じ挙動を示した。また、SLC19A1 c.80G>AとSLC19A1 c.696T>Cも同様であった。FPGS g.130562725G>AとFPGS c.*192A>Gについては、挙動は同じではないが、非常に近いことが分かった。また、表５に示されるように、多重共線性の関係にある２つのSNP(SNP1及びSNP2)には、強い相関性があった。よって、後述の統計解析では、多重共線性の関係にある２つのSNPのうち、いずれか一方のみを用いればよい。本比較例１及び後述の実施例１〜３では、FPGS c.*192A>G、GGH c.-401C>T及びSLC19A1 c.696T>Cを選択した。
【００９３】
(４)統計解析
FPGS c.64A>G、FPGS c.*192A>G、GGH c.452C>T、GGH c.-401C>T及びSLC19A1 c.696T>Cの５つのSNPについて、被験者49名のデータを用いて多重ロジスティック解析を行い、上記の回帰式(II)に基づいて判別式を構築した(以下、構築した判別式を「判別式１」と呼ぶ)。多重ロジスティック回帰分析は、統計計算用ソフトウェア「Ｒ」によって行った。「Ｒ」はhttp://www.r-project.org/からダウンロードできる。この判別式１では、変異アレル数を名義尺度(WT、HTまたはMT)として扱った。判別式１の係数を、表６に示す。これらの係数は、統計解析ソフトウェア「JMP10」(SAS社製)を用いて計算した。表６中、FPGS c.64A>GについてWTの係数が空欄になっているが、これは、被験者49名においてWTに該当する者がいなかったからである。なお、FPGS遺伝子のc.64位について、日本人のGenotype頻度をHapMap (http://hapmap.ncbi.nlm.hih.gov/index.html.ja)で検索したところ、WTは０％、HTは9.8％、MTは90.2％であった。
【００９４】
【表６】
0006542037.tif 000009
【００９５】
判別式１を用いて、ROC解析を行った。具体的には、判別式１により上記(１)で述べた有効群及び無効群について、MTXの有効性の判定を行った。判別式１のb1〜b5に、各SNPについて対立遺伝子のいずれにも変異がない場合は表６の「WT」欄の係数を代入し、対立遺伝子のいずれか一方に変異がある場合は「HT」欄の係数を代入し、対立遺伝子の遺伝子多型部位の両方に変異がある場合は「MT」欄の係数を代入し、有効性予測値を算出した。有効性予測値が0.5以上となった場合はMTXが有効であると判定し、0.5未満となった場合はMTXが有効でないと判定した。また、Ｘ軸を「１−特異度」の値(すなわち、偽陽性率)、Ｙ軸を「感度」(すなわち、陽性率)として、ROC曲線を作成し、ROC曲線下面積(AUC)を算出した。なお、ROC曲線は、JMP10(SAS社製)により作成した。ROC曲線を図４に示す。判定結果の感度及び特異度、並びに図４のROC曲線のAUCを以下の表７に示す。
【００９６】
【表７】
0006542037.tif 000010
【００９７】
(５)判定結果
表７に示されるように、判別式１による判定は、感度が55％、特異度が89％、AUCが0.75であり、十分に精度の高い判定とはいえない。MTX代謝に関連する上記の３遺伝子のSNPだけでは、MTXの有効性を精度よく判定することは難しいことが示唆された。なお、GGH c.-401C>Tは、上記の被験者における挙動がGGH c.16T>Cと完全に一致しているので、GGH c.-401C>TをGGH c.16T>Cに代えても、同じ判定結果が得られると考えられる。同様に、SLC19A1 c.696T>CをSLC19A1 c.80G>Aに代えても、同じ判定結果が得られると考えられる。また、FPGS c.*192A>Gは、上記の被験者における挙動がFPGS g.130562725C>Tと非常に近いので、FPGS c.*192A>GをFPGS g.130562725G>Aに代えても、同様の判定結果が得られると考えられる。
【００９８】
実施例１
FPGS遺伝子、GGH遺伝子及びSLC19A1遺伝子のSNPに加えて、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPをさらに解析し、その解析結果に基づいて、RA患者におけるMTXの有効性を判定することを試みた。
【００９９】
(１)被験者及び生体試料
実施例１では、比較例１で取得した被験者(有効群22名及び無効群27名)のゲノムDNAを、SNPの解析に用いた。
【０１００】
(２)遺伝子多型解析
各被験者のゲノムDNAについて、上記の７遺伝子に存在する９つのSNPを、TaqMan(商標) SNP Genotyping Assays(ライフテクノロジーズ社)及びダイレクトシーケンス法を用いて解析した。上記の９つのSNP及びそれらの解析方法を、表８に示す。FPGS c.64A>Gの解析には、上記のTaqMan(商標) SNP Genotyping Assaysのカスタム製品を用いた。
【０１０１】
【表８】
0006542037.tif 000011
【０１０２】
(2-1)TaqMan(商標) SNP Genotyping Assaysによる解析
各被験者のゲノムDNA(10 ng)について、上記のTaqMan(商標) SNP Genotyping AssaysとTaqMan(商標) GTXpress Master Mix(ライフテクノロジーズ社)を用いて、Real-Time PCR法にて測定を行った。なお、具体的な操作は、ライフテクノロジーズ社が提供するプロトコールに従って行った。
【０１０３】
(2-2)ダイレクトシーケンス法による解析
各被験者のゲノムDNAを鋳型核酸として用いて、遺伝子多型部位を含むFMO2遺伝子をPCR反応により増幅した。増幅に用いたフォワードプライマーおよびリバースプライマーの配列を以下に示す。
【０１０４】
＜FMO2 c.585A>G＞
・フォワードプライマー：5'- TCCAGAAAGGAAAAGCTGGCAATG -3'  (配列番号２１)
・リバースプライマー：5'- GAGCCATCTCCCAGAGAAGTGAA -3'  (配列番号２２)
【０１０５】
PCRの反応系には、ゲノムDNA(２ng／反応)、各プライマーおよびポリメラーゼとしてTaKaRa Ex-TaqまたはTaKaRa LA-Taq(タカラバイオ株式会社)を添加した。ポリメラーゼとしてTaKaRa LA-Taqを用いた反応液は、該ポリメラーゼに添付のプロトコールに従って調製した。ポリメラーゼとしてEx-Taqを用いた反応液は、該ポリメラーゼに添付のプロトコールに記載された組成の反応液にDMSOを終濃度５％(V/V)となるように添加して調製した。反応終了後、反応液をAmpure XP(Beckman Coulter社)を用いて精製した。精製した反応液の一部をアガロース電気泳動に付して、増幅産物の存在を確認した。そして、以下に示すシーケンス用プライマーを用いて、増幅産物の塩基配列を解析した。
【０１０６】
＜FMO2 c.585A>G＞
・シーケンス用プライマー：5'- TGACATAGTTGCTCTGGAGC -3'  (配列番号２３)
【０１０７】
(３)統計解析
表８に示す９つのSNPについて、被験者49名のデータを用いて多重ロジスティック解析を行い、上記の回帰式(II)に基づいて判別式を構築した(以下、構築した判別式を「判別式２」と呼ぶ)。多重ロジスティック回帰分析は、比較例１と同様にして行った。判別式２の係数を、表９に示す。これらの係数は、JMP10(SAS社製)を用いて計算した。
【０１０８】
【表９】
0006542037.tif 000012
【０１０９】
判別式２を用いて、ROC解析を行った。具体的には、判別式２により、有効群及び無効群について、比較例１と同様にしてMTXの有効性の判定を行った。有効性予測値が0.5以上となった場合はMTXが有効であると判定し、0.5未満となった場合はMTXが有効でないと判定した。また、ROC曲線を作成してAUCを算出した。なお、ROC曲線は、JMP10(SAS社製)により作成した。ROC曲線を図５に示す。判定結果の感度及び特異度、並びに図５のROC曲線のAUCを以下の表１０に示す。
【０１１０】
【表１０】
0006542037.tif 000013
【０１１１】
(５)判定結果
表１０に示されるように、判別式２による判定は、感度が100％、特異度が85％、AUCが0.98であり、判定の精度が極めて高いことがわかった。MTX代謝に関連する３遺伝子のSNPに加えて、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPをさらに解析することにより、MTXの有効性の判定精度が顕著に上昇することが示された。なお、GGH c.-401C>Tは、上記の被験者における挙動がGGH c.16T>Cと完全に一致しているので、GGH c.-401C>TをGGH c.16T>Cに代えても、同じ判定結果が得られると考えられる。同様に、SLC19A1 c.696T>CをSLC19A1 c.80G>Aに代えても、同じ判定結果が得られると考えられる。
【０１１２】
実施例２
MTXの代謝に関連する３遺伝子のうちFPGS遺伝子及びGGH遺伝子のSNPに加えて、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPを解析し、その解析結果に基づいて、RA患者におけるMTXの有効性を判定することを試みた。
【０１１３】
(１)被験者及び生体試料
実施例２では、比較例１で取得した被験者(有効群22名及び無効群27名)のゲノムDNAを、SNPの解析に用いた。
【０１１４】
(２)遺伝子多型解析
各被験者のゲノムDNAについて、上記の６遺伝子に存在する８つのSNPを、実施例１と同様にして解析した。上記の８つのSNPは、FPGS c.64A>G、FPGS c.*192A>G、GGH c.452C>T、GGH c.-401C>T、FMO2 c.585A>G、SLC28A3 c.267G>A、EPHX1 c.357G>A及びPPARG c.*1411G>Cである。
【０１１５】
(３)統計解析
上記の８つのSNPについて、被験者49名のデータを用いて多重ロジスティック解析を行い、上記の回帰式(II)に基づいて判別式を構築した(以下、構築した判別式を「判別式３」と呼ぶ)。多重ロジスティック回帰分析は、比較例１と同様にして行った。判別式３の係数を、表１１に示す。これらの係数は、JMP10(SAS社製)を用いて計算した。
【０１１６】
【表１１】
0006542037.tif 000014
【０１１７】
判別式３を用いて、ROC解析を行った。具体的には、判別式３により、有効群及び無効群について、比較例１と同様にしてMTXの有効性の判定を行った。有効性予測値が0.5以上となった場合はMTXが有効であると判定し、0.5未満となった場合はMTXが有効でないと判定した。また、ROC曲線を作成してAUCを算出した。なお、ROC曲線は、JMP10(SAS社製)により作成した。ROC曲線を図６に示す。判定結果の感度及び特異度、並びに図６のROC曲線のAUCを以下の表１２に示す。
【０１１８】
【表１２】
0006542037.tif 000015
【０１１９】
(５)判定結果
表１２に示されるように、判別式３による判定は、感度が95％、特異度が89％、AUCが0.97であり、判定の精度が極めて高いことがわかった。MTX代謝に関連する３遺伝子のうちFPGS遺伝子及びGGH遺伝子のSNPと、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPとを解析することにより、精度の高い判定が可能であることが示唆された。なお、GGH c.-401C>Tは、上記の被験者における挙動がGGH c.16T>Cと完全に一致しているので、GGH c.-401C>TをGGH c.16T>Cに代えても、同じ判定結果が得られると考えられる。
【０１２０】
実施例３
MTXの代謝に関連する３遺伝子のうちFPGS遺伝子のSNPに加えて、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPを解析し、その解析結果に基づいて、RA患者におけるMTXの有効性を判定することを試みた。
【０１２１】
(１)被験者及び生体試料
実施例３では、比較例１で取得した被験者(有効群22名及び無効群27名)のゲノムDNAを、SNPの解析に用いた。
【０１２２】
(２)遺伝子多型解析
各被験者のゲノムDNAについて、上記の５遺伝子に存在する６つのSNPを、実施例１と同様にして解析した。上記の６つのSNPは、FPGS c.64A>G、FPGS c.*192A>G、FMO2 c.585A>G、SLC28A3 c.267G>A、EPHX1 c.357G>A及びPPARG c.*1411G>Cである。
【０１２３】
(３)統計解析
上記の６つのSNPについて、被験者49名のデータを用いて多重ロジスティック解析を行い、上記の回帰式(II)に基づいて判別式を構築した(以下、構築した判別式を「判別式４」と呼ぶ)。多重ロジスティック回帰分析は、比較例１と同様にして行った。判別式４の係数を、表１３に示す。これらの係数は、JMP10(SAS社製)を用いて計算した。
【０１２４】
【表１３】
0006542037.tif 000016
【０１２５】
判別式４を用いて、ROC解析を行った。具体的には、判別式４により、有効群及び無効群について、比較例１と同様にしてMTXの有効性の判定を行った。有効性予測値が0.5以上となった場合はMTXが有効であると判定し、0.5未満となった場合はMTXが有効でないと判定した。また、ROC曲線を作成してAUCを算出した。なお、ROC曲線は、JMP10(SAS社製)により作成した。ROC曲線を図７に示す。判定結果の感度及び特異度、並びに図７のROC曲線のAUCを以下の表１４に示す。
【０１２６】
【表１４】
0006542037.tif 000017
【０１２７】
(５)判定結果
表１４に示されるように、判別式４による判定は、感度が91％、特異度が89％、AUCが0.95であり、判定の精度が極めて高いことがわかった。MTX代謝に関連する３遺伝子のうちFPGS遺伝子のSNPと、FMO2遺伝子、SLC28A3遺伝子、EPHX1遺伝子及びPPARG遺伝子のSNPとを解析することにより、精度の高い判定が可能であることが示唆された。
【０１２８】
実施例４
実施例１〜３のそれぞれで用いたSNPのうち、FPGS c.*192A>Gに代えて、FPGS g.130562725G>Aを用いて解析を行い、その解析結果に基づいて、RA患者におけるMTXの有効性を判定した。
【０１２９】
(１)被験者及び生体試料
実施例４では、比較例１で取得した被験者(有効群22名及び無効群27名)のゲノムDNAを、SNPの解析に用いた。
【０１３０】
(２)遺伝子多型解析
各被験者のゲノムDNAについて、FPGS g.130562725G>Aを比較例１と同様にして解析し、FPGS c.64A>G、GGH c.452C>T、GGH c.-401C>T、FMO2 c.585A>G、SLC28A3 c.267G>A、EPHX1 c.357G>A及びPPARG c.*1411G>Cを、実施例１と同様にして解析した。
【０１３１】
(３)統計解析
実施例４では、以下の表１５に示されるSNPの組み合わせを用いた。表１５中、「●」が付されたSNPの解析データが、後述の判別式に用いられたことを示す。
【０１３２】
【表１５】
0006542037.tif 000018
【０１３３】
表１５に示されるSNPについて、被験者49名のデータを用いて多重ロジスティック解析を行い、上記の回帰式(II)に基づいて判別式を構築した(以下、構築した判別式を、それぞれ「判別式５」、「判別式６」及び「判別式７」と呼ぶ)。多重ロジスティック回帰分析は、比較例１と同様にして行った。各判別式の係数を、表１６〜１８に示す。これらの係数は、JMP10(SAS社製)を用いて計算した。
【０１３４】
【表１６】
0006542037.tif 000019
【０１３５】
【表１７】
0006542037.tif 000020
【０１３６】
【表１８】
0006542037.tif 000021
【０１３７】
判別式５〜７のそれぞれを用いて、ROC解析を行った。具体的には、各判別式により、有効群及び無効群について、比較例１と同様にしてMTXの有効性の判定を行った。有効性予測値が0.5以上となった場合はMTXが有効であると判定し、0.5未満となった場合はMTXが有効でないと判定した。また、ROC曲線を作成してAUCを算出した。なお、ROC曲線は、JMP10(SAS社製)により作成した。判別式５〜７のROC曲線をそれぞれ図８〜１０に示す。判定結果の感度及び特異度、並びにROC曲線のAUCを以下の表１９に示す。
【０１３８】
【表１９】
0006542037.tif 000022
【０１３９】
表１９に示されるように、判別式５〜７のいずれを用いても、判定の精度が極めて高いことがわかった。よって、実施例１〜３のそれぞれで用いたSNPのうち、FPGS c.*192A>Gに代えて、FPGS g.130562725G>Aを用いて解析を行っても、RA患者におけるMTXの有効性を精度よく判定できることが示された。
【符号の説明】
【０１４０】
１０      判定装置
２０      測定装置
３０      コンピュータシステム
４０      記録媒体
５０      試薬キット
５１〜５３    試薬容器
５４      添付文書
５５      梱包箱
３００    コンピュータ本体
３０１    入力部
３０２    表示部
３１０    ＣＰＵ
３１１    ＲＯＭ
３１２    ＲＡＭ
３１３    ハードディスク
３１４    入出力インターフェイス
３１５    読出装置
３１６    通信インターフェイス
３１７    画像出力インターフェイス
３１８    バス
【図１】
0006542037.tif 000023
【図２】
0006542037.tif 000024
【図３】
0006542037.tif 000025
【図４】
0006542037.tif 000026
【図５】
0006542037.tif 000027
【図６】
0006542037.tif 000028
【図７】
0006542037.tif 000029
【図８】
0006542037.tif 000030
【図９】
0006542037.tif 000031
【図１０】
0006542037.tif 000032
【図１１】
0006542037.tif 000033
【配列表】
0006542037000001.app