"第6染色体短腕22領域または第9染色体長腕21領域に存在する一塩基多型を分析し、該分析結果に基づいて肥満を検査することを特徴とする、肥満の発症リスクおよび/または発症の有無の判定方法。前記一塩基多型が、配列番号1または2の塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、請求項1に記載の方法。前記連鎖不平衡の関係にある塩基が、配列番号3~72から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基である、請求項2に記載の方法。配列番号1~72から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列を有する肥満検査用プローブ。配列番号1~72から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域を増幅することのできる肥満検査用プライマー。本発明は、肥満の発症リスクや発症の有無を正確に検査する方法、及び該方法に用いられる検査試薬を提供することを課題とする。本発明者らは上記課題の解決のために鋭意検討した結果、第6染色体短腕22領域または第9染色体長腕21領域に存在する一塩基多型SNPが肥満の指標であるBMIと関連することを同定した。そして、これらの多型を調べることにより肥満の発症リスクや発症の推定を正確に実施できることを見出し、本発明を完成するに至った。すなわち、本発明は以下の通りである。1第6染色体短腕22領域または第9染色体長腕21領域に存在する一塩基多型を分析し、該分析結果に基づいて肥満を検査することを特徴とする、肥満の発症リスクおよび/または発症の有無の判定方法。2前記一塩基多型が、配列番号1または2の塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、1に記載の方法。3前記連鎖不平衡の関係にある塩基が、配列番号3~72から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基である、2に記載の方法。4配列番号1~72から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列を有する肥満検査用プローブ。5配列番号1~72から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域を増幅することのできる肥満検査用プライマー。本発明によれば、これまで予測が困難であった肥満の発症リスク罹患リスクを正確かつ簡便に予測することができる。また、肥満の発症を正確かつ簡便に判定することができる。したがって、本発明は肥満の予防や早期治療に貢献するものである。<1>本発明の方法本発明の方法は、ヒトの第6染色体短腕22領域または第9染色体長腕21領域に存在するSNPを分析し、該分析結果に基づいて肥満を検査することを特徴とする、肥満の発症リスクおよび/または発症の有無の判定方法である。すなわち、本発明において、検査とは肥満の発症リスクの検査及び肥満の発症の有無の検査を含む。本発明の方法においては、SNPの分析結果を、肥満の発症リスクおよび/または発症の有無と関連付ける。ヒト第6染色体短腕22領域6p22領域に存在する具体的なSNPとしては、ヒトrs2206734を挙げることができる。ここで、rs番号はNationalCenterforBiotechnologyInformationNCBIのdbSNPデータベースhttpwwwncbinlmnihgovprojectsSNPの登録番号を示す。rs2206734はGenBankAccessionNoNT00759214の11553135番目の塩基におけるシトシンC/チミンTの多型を意味し、この塩基がCである場合は肥満の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CT>TTの順で肥満の可能性または発症リスクが高い。rs2206734は、CDKAL1遺伝子座に存在し、CDKAL1遺伝子座は肥満のリスク遺伝子座であると考えられる。よって、特に、CDKAL1遺伝子座に存在するSNPを解析することによって肥満を検査することができる。CDKAL1遺伝子としては、具体的には、GenBankAccessionNoNC00000611の20534688~21232635の領域が挙げられる。ヒト第9染色体長腕21領域9q21領域に存在する具体的なSNPとしては、ヒトrs11142387を挙げることができる。rs11142387はGenBankAccessionNoNT02393517の2162864番目の塩基におけるアデニンA/シトシンCの多型を意味し、この塩基がCである場合は肥満の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CA>AAの順で肥満の可能性または発症リスクが高い。rs11142387は、KLF9遺伝子のプロモーター領域に存在し、KLF9遺伝子座は肥満のリスク遺伝子座であると考えられる。よって、特に、KLF9遺伝子座に存在するSNPを解析することによって肥満を検査することができる。KLF9遺伝子としては、具体的には、GenBankAccessionNoNC00000911の72999513~73029573の領域の相補配列が挙げられる。なお、rs2206734およびrs11142387について、SNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、それぞれ配列番号1および2に示した。61番目の塩基が多型を有する。本発明においては、上記塩基に相当する塩基を解析する。上記塩基に相当する塩基とは、上記領域における該当塩基を意味する。すなわち、上記塩基に相当する塩基を解析することには、仮に人種の違いなどによって上記配列がSNP以外の位置で若干変化したとしても、上記領域における該当塩基を解析することが含まれる。また、本発明において解析する塩基は上記のものに限定されず、上記の塩基と連鎖不平衡にある塩基の多型を分析してもよい。ここで上記の塩基と連鎖不平衡にある塩基とは、上記の塩基とr2>0.5、好ましくはr2>0.8、さらに好ましくはr2>0.9の関係を満たす塩基をいう。r2は連鎖不平衡係数である。また上記の塩基と連鎖不平衡にある塩基は、例えば、HapMapデータベース等を用いて同定することができる。もしくは、複数人通常は20-40人程度から採取したDNAをシークエンサーにて配列解析し、連鎖不平衡にあるSNPを探索することにより同定することもできる。rs2206734とr2>0.5で連鎖不平衡にある塩基としては、rs9368219、rs7766070、rs9356744、rs9368222、rs9350271、rs10440833、rs6456364、rs2328545、rs742642、rs9295475、rs9348440、rs9358355、rs9368216、rs7752906、rs10946398、rs4710940、rs4712522、rs4712525、rs4712526、rs6456368、rs6456369、rs6906327、rs7748382、rs7752780、rs7754840、rs7774594、rs9460544、rs9460545、rs9460546、rs6456367、rs7772603、rs9358356、rs4712523、rs9295474、rs9465850、rs2328528、rs6940200、rs7751957、rs7768642、rs9465846、rs9465847、rs9356747、rs9356746、rs4515379、rs9465837、rs9465841、rs2328548、rs1040558、rs10946403、rs11753081、rs6935599、rs7747752、rs7767391、rs9295478、rs9356748、rs9358357、rs9358358、rs9368224、rs12111351、rs9460550、rs9368226、rs9465871、rs2328529、rs7756992、rs1569699、rs2328549、rs9350270、およびrs9356743が挙げられる。これらSNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、それぞれ配列番号3~70に示した。61番目の塩基が多型を有する。これらSNPは、リスクアレルのホモ接合体>リスクアレルと非リスクアレルのへテロ接合体>非リスクアレルのホモ接合体の順で肥満の可能性または発症リスクが高い。rs11142387とr2>0.5で連鎖不平衡にある塩基としては、rs10780928およびrs4745008が挙げられる。これらSNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、それぞれ配列番号71および72に示した。61番目の塩基が多型を有する。これらSNPは、リスクアレルのホモ接合体>リスクアレルと非リスクアレルのへテロ接合体>非リスクアレルのホモ接合体の順で肥満の可能性または発症リスクが高い。rs2206734またはrs11142387とr2>0.5で連鎖不平衡にあるこれら各塩基について、rs2206734またはrs11142387に対する連鎖不平衡係数r2、アレルの組み合わせ、およびリスクアレルを表1、2に示す。表1、2に示した連鎖不平衡係数r2は、HapMapデータベース上で公開されている、日本人および中国人90人のジェノタイプデータPhaseII、release24に基づき推定されたものである。リスクアレルの表記はNCBIが公開している参照配列Build36におけるforwardstrandに対応している。上記SNPの塩基の種類を調べ、得られた結果を上記のような基準に基づいて肥満と関連付けることにより、肥満を検査することができる。上記SNPは単独で解析されてもよいし、上記SNPの少なくとも1つを含む複数のSNPsをまとめて解析ハプロタイプ解析してもよい。例えば、上記SNPの複数をまとめて解析してもよいし、上記SNPの少なくとも1つと、肥満と関連する既知のSNPs非特許文献3~12や当該既知のSNPsと連鎖不平衡にあるSNPsとを組み合わせて解析してもよい。肥満と関連する複数のSNPsをまとめて解析すれば、肥満の検査の精度が向上する。なお、いずれのSNPも、二本鎖DNAのどちらの鎖を解析してもよい。例えば、CDKAL1遺伝子またはKLF9遺伝子の配列はセンス鎖を解析してもよいし、アンチセンス鎖を解析してもよい。SNPの解析に用いる試料としては、染色体DNAを含む試料であれば特に制限されないが、例えば、血液、尿等の体液サンプル、口腔粘膜などの細胞、毛髪等の体毛などが挙げられる。SNPの解析にはこれらの試料を直接使用することもできるが、これらの試料から染色体DNAを常法により単離し、これを用いて解析することが好ましい。SNPの解析は、通常の遺伝子多型解析方法によって行うことができる。例えば、シークエンス解析、PCR、ハイブリダイゼーション、インベーダー法などが挙げられるが、これらに限定されない。シークエンス解析は通常の方法により行うことができる。具体的には、多型を示す塩基の5’側数十塩基の位置に設定したプライマーを使用してシークエンス反応を行い、その解析結果から、該当する位置がどの種類の塩基であるかを決定することができる。なお、シークエンス反応の前に、あらかじめSNP部位を含む断片をPCRなどによって増幅しておくことが好ましい。また、SNPの解析は、PCRによる増幅の有無を調べることによって行うことができる。例えば、多型を示す塩基を含む領域に対応する配列を有し、かつ、3’末端が各多型に対応するプライマーをそれぞれ用意する。それぞれのプライマーを使用してPCRを行い、増幅産物の有無によってどのタイプの多型であるかを決定することができる。また、LAMP法特許第3313358号明細書、NASBA法NucleicAcidSequenceBasedAmplification;特許2843586号明細書、ICAN法特開2002-233379号公報などによって増幅の有無を調べることもできる。その他、単鎖増幅法を用いてもよい。また、SNP部位を含むDNA断片を増幅し、増幅産物の電気泳動における移動度の違いによってどのタイプの多型であるかを決定することもできる。このような方法としては、例えば、PCR-SSCPsinglestrandconformationpolymorphism法Genomics1992Jan112139146が挙げられる。具体的には、まず、目的のSNPを含むDNAを増幅し、増幅したDNAを一本鎖DNAに解離させる。次いで、解離させた一本鎖DNAを非変性ゲル上で分離し、分離した一本鎖DNAのゲル上での移動度の違いによってどのタイプの多型であるかを決定することができる。さらに、多型を示す塩基が制限酵素認識配列に含まれる場合は、制限酵素による切断の有無によって解析することもできるRFLP法。この場合、まず、DNA試料を制限酵素により切断する。次いで、DNA断片を分離し、検出されたDNA断片の大きさによってどのタイプの多型であるかを決定することができる。また、ハイブリダイゼーションの有無を調べることによって多型の種類を解析することも可能である。すなわち、各塩基に対応するプローブを用意し、いずれのプローブにハイブリダイズするかを調べることによってSNPがいずれの塩基であるかを調べることもできる。このようにしてSNPがいずれの塩基であるかを決定することで、肥満を検査するためのデータを得ることができる。<2>本発明の検査用試薬本発明はまた、肥満を検査するためのプライマーやプローブなどの検査試薬を提供する。このようなプローブとしては、上記SNP部位を含み、ハイブリダイズの有無によってSNP部位の塩基の種類を判定できるプローブが挙げられる。具体的には、配列番号1~72から選ばれる塩基配列の61番目の塩基を含む配列、又はその相補配列を有する10塩基以上の長さのプローブが挙げられる。プローブの長さは好ましくは、15~35塩基であり、より好ましくは20~35塩基である。また、プライマーとしては、上記SNP部位を増幅するためのPCRに用いることのできるプライマー、又は上記SNP部位を配列解析シークエンシングするために用いることのできるプライマーが挙げられる。具体的には、配列番号1~72から選ばれる塩基配列の61番目の塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーが挙げられる。このようなプライマーの長さは10~50塩基が好ましく、15~35塩基がより好ましく、20~35塩基がさらに好ましい。上記SNP部位をシークエンシングするためのプライマーとしては、上記塩基の5’側領域、好ましくは30~100塩基上流の配列を有するプライマーや、上記塩基の3’側領域、好ましくは30~100塩基下流の領域に相補的な配列を有するプライマーが例示される。PCRによる増幅の有無で多型を判定するために用いるプライマーとしては、上記塩基を含む配列を有し、上記塩基を3’側に含むプライマーや、上記塩基を含む配列の相補配列を有し、上記塩基の相補塩基を3’側に含むプライマーなどが例示される。なお、本発明の検査用試薬はこれらのプライマーやプローブに加えて、PCR用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。以下、本発明を実施例によりさらに具体的に説明する。但し、本発明はこれらの実施例に限定されない。1肥満と関連するSNPsの同定肥満感受性を決定する遺伝的変異を同定するために、日本人被検者を用いてゲノムワイド関連解析GWASを行った。GWASとは、疾患等の表現型に関わる遺伝的変異を探索する遺伝統計学的手法である。例えば、ヒトゲノム全体を網羅するような数十万~100万ヶ所のSNPsを用いて、ある形質の平均値の分布において、多型の組み合わせジェノタイプによる差があるかどうかを統計的に検定することで、疾患と関連する遺伝的変異を見出すことができる。本実施例において、肥満の指標としては、ボディマス指数BodyMassIndex;BMIを用いた。<被検者>GWASに用いた26620名の被検者、再現試験replicationstudy1に用いた3763名の被検者、および再現試験2に用いた4147名の被検者は、東京大学医科学研究所のBioBankJapanに登録された種々の疾患の日本人患者から採用した。これら被検者のうち、BMI≧27.5の被検者を、肥満の被検者ケースcaseとして用いた。それ以外の被検者を、対照被検者コントロールcontrolとして用いた。GWAS、再現試験1、および再現試験2に用いた被検者における肥満の被検者の比率は、それぞれ9.3、7.2、7.3であった。再現試験3には、27715名の東アジア人種の8集団からなる被検者を用いた。GWAS、および再現試験1~3に用いた被検者集団は、互いに独立である。本実施例に用いた被検者の特徴を表3に示す。a:BMI≧27.5の被検者を、肥満obesityとした。b:BioBankJapanに登録された種々の疾患の患者からなる被検者。c:東アジア人種の8集団KAREShanghaiSiMESSDCSSP2IMCJTaiwanCRCからなる被検者。本研究は東京大学医科学研究所および理化学研究所横浜研究所の倫理委員会によって承認され、全ての参加者からインフォームドコンセントを得た。<GWAS>被検者26620名の遺伝子型を、IlluminaHuman610QuadGenotypingBeadChipIllumina社USAを用いて解析した。なお、クオリティコントロールとして、コール率callrateが0.98未満の被検者、コール率callrateが0.99未満のSNPs、あいまいなコールambiguouscallのSNPs、identifybystateIBS>1.7の近縁ペアの内の低コール率の被検者、および主成分分析PCAにより外れ値outlierであると判定された被検者を除外した。クオリティコントロールを通過したマイナーアレル頻度MAFs≧0.01を満たす常染色体上の480,103個のSNPsについて、遺伝子型のデータを取得した。各SNPsについて、線形回帰モデルlinearregressionmodelを用いてBMIとの関連を評価した。GWASの結果、ゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<5.0×108を満たしてBMIと有意に関連する3ヶ所のSNPs9q21領域のrs11142387、11p14領域のrs2030323、および19q13領域のrs11671664が見出された図1a、表4。a:GWASおよび全3回の再現試験の複合解析でP<5.0×105を満たしたSNPsを示す。b:BMIを増大させるアレルをA1アレル1として、NCBIBuild36のforwardstrandに基づいて示した。c:アレル1の頻度。d:標準化したBMI平均値=0、標準偏差=1におけるアレル1の効果量。e:独立した3度の再現試験の総合結果。f:非特許文献12の結果に基づく。<再現試験>GWASで見出された関連性を検証するため、独立した3度の再現試験replicationstudyを行った。再現試験1では、GWASにおいてP<5.0×105を満たした36個のSNPsについて、BMIとの関連を検証した。被検者3763名の遺伝子型をIlluminaHumanHap550v3GenotypingBeadChipIllumina社USAを用いて解析し、GWASと同様のクオリティコントロールの下でBMIとの関連を検証した。GWASと再現試験1の複合解析combinedanalysisの結果、P<5.0×105を満たした11個のSNPsを再現試験2および3に供した。再現試験2では、被検者4147名の遺伝子型を、IlluminaHumanOmniExpressGenotypingBeadChipIllumina社USAを用いて解析し、GWASと同様のクオリティコントロールの下でBMIとの関連を検証した。再現試験3では、被検者27715名の遺伝子型をAffymetrixGenomeWideHumanSNPArray60Affymetrix社USA、AffymetrixGenomeWideHumanSNParray50Affymetrix社USA、Illumina1MduoGenotypingBeadChipIllumina社USA、IlluminaHumanHap610QuadGenotypingBeadChipIllumina社USA、およびIlluminaHumanHap550GenotypingBeadChipIllumina社USAを用いて解析した。GWASと全3回の再現試験の複合解析combinedanalysisの結果、P<5.0×105を満たした9個のSNPsを表4に示す。これらの内、7か所のSNPsが、ゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<5.0×108を満たし、BMIと有意に関連することが確認された。なお、これら7個のSNPsには、GWASにおいて有意性を満たした3個のSNPsも含まれている。これら7個のSNPsが位置する遺伝子座の内、1q25領域のSEC16B遺伝子座、11p14領域のBDNF遺伝子座、16q12領域のFTO遺伝子座、18p21領域のMC4R遺伝子座、および19q13領域のGIPR遺伝子座の5つについては、すでにヨーロッパ人においてBMIとの関連が見出されている。一方、残りのrs2206734およびrs11142387が位置する領域6p22および9q21については、主にヨーロッパ人被検者特に成人被検者を用いたこれまでの研究ではBMIとの関連が報告されていなかった。rs2206734の位置する領域におけるSNPsのプロットおよび遺伝子構造を図1bに示す。また、rs11142387の位置する領域におけるSNPsのプロットおよび遺伝子構造を図1cに示す。rs2206734は、CDKAL1遺伝子座のコード領域に位置している。CDKAL1遺伝子は、サイクリン依存性キナーゼ5制御サブユニット関連タンパク質1ライク1cyclindependentkinase5regulatorysubunitassociatedprotein1like1をコードする遺伝子である。CDKAL1遺伝子座は、ドイツ人集団を用いた研究で8歳児のBMIと関連することが報告されている非特許文献13。また、CDKAL1遺伝子座は、デンマーク人集団を用いた研究で出生時体重と関連することが報告されている非特許文献14。しかしながら、アジア人集団におけるCDKAL1遺伝子座とBMIとの関連は報告されておらず、また、成人におけるCDKAL1遺伝子座とBMIとの関連も報告されていなかった。よって、特に、CDKAL1遺伝子座のSNPsを解析することでアジア人種の肥満や成人の肥満を検査できることは、本発明者により初めて見出された知見である。また、CDKAL1遺伝子座は2型糖尿病のリスク遺伝子座であることが知られている非特許文献13~15。そこで、rs2206734およびrs11142387について、共変数を含むロジスティック回帰モデルlogisticregressionmodelを用いて肥満および2型糖尿病との関連を評価した。結果、rs2206734CDKAL1遺伝子座のTアレルは、BMIを低下させるのに対し、2型糖尿病のリスクを増大させることが明らかとなった表5。通常、BMIが増大すれば2型糖尿病のリスクも増大すると考えられる。よって、rs2206734のBMIに対する効果は2型糖尿病に対する効果から想定されるものと反対であり、rs2206734のCアレルがBMIのリスクアレルとして見出されたのは驚くべき結果であった。a:BMIを増大させるアレルをA1アレル1として、NCBIBuild36のforwardstrandに基づいて示した。b:GWASと全3回の再現試験の被検者が含まれる。c:肥満および2型糖尿病におけるアレル1の効果量。rs11142387は、KLF9遺伝子のプロモーター領域に位置している。rs11142387を含む連鎖不平衡LDブロックは複数の遺伝子を含むが、タグSNPであるrs11142387がKLF9遺伝子座に位置していることから、KLF9遺伝子が原因遺伝子の候補である。KLF9遺伝子は、ジンクフィンガー転写因子の1つをコードする遺伝子であり、種々の生理反応に関与する。ここで、脂肪細胞分化の重要な要素としてperoxisomeproliferatoractivatedreceptorPPARγが知られており、PPARγは肥満を含む多くの代謝特性の病態と関連しているPeiHetalCellDeathDiffer18315327。最近の研究では、KLF9遺伝子がコードするタンパク質は、PPARγのトランス活性化を介した脂肪生成促進性proadipogenicの転写因子であることが示唆されている非特許文献16。また、KLFとしては、KLF7遺伝子座が肥満と関連すること、およびKLF5遺伝子がコードするタンパク質が脂肪細胞分化を制御していることが知られている。しかしながら、KLF9遺伝子座とBMIとの関連は報告されていなかった。なお、rs2206734CDKAL1遺伝子座と異なり、rs11142387KLF9遺伝子座は2型糖尿病との関連は認められなかった表5。<民族間での比較>さらに、民族間での肥満の遺伝学的差異を評価するため、GWASと全3回の再現試験の複合解析combinedanalysisでP<5.0×105を満たした9個のSNPsについて、GIANTコンソーシアムの123,865名のヨーロッパ人被検者のメタ解析の結果非特許文献12を利用してBMIとの関連を評価した表4。9個のSNPsの内、rs2206734CDKAL1遺伝子座を含む5つが、ヨーロッパ人集団においてもBMIと有意P<0.028、falsediscoveryrate<0.05に関連した。一方、rs11142387KLF9遺伝子座については、ヨーロッパ人集団においてはBMIとの関連が認められなかった。よって、特にrs11142387KLF9遺伝子座および当該塩基と連鎖不平衡の関係にあるSNPsは、アジア人集団における肥満の検査に好適に用いることができる。次に、主にヨーロッパ人集団においてBMIとの関連が報告されている既知のSNPs非特許文献13~15、ZobelDPetalEurJEndocrinol160603609について、本実施例のGWASで得られたデータを用いて日本人集団におけるBMIとの関連を評価した表6。a:BMIとの関連が報告されている既知のSNPs。b:BMIを増大させるアレルをA1アレル1として、NCBIBuild36のforwardstrandに基づいて示した。c:GWAS被検者におけるアレル1の頻度。d:標準化したBMI平均値=0、標準偏差=1におけるアレル1の効果量。e:P<0.02多型遺伝子座の数に基づくfalsediscoveryrate<0.05を有意とした。f:統計的検出力statisticalpowerは、アレル頻度とGWASのサンプルサイズに基づいて算出し、効果量を0.035と仮定した条件での有意性の閾値を0.02とした。BMIとの関連が報告されている既知のSNPsの内、本実施例で既にBMIとの関連が再確認された5つの遺伝子座に加え、さらに10ヶ所の遺伝子座が日本人被検者においても同様にBMIに影響することが明らかとなったP<0.02falsediscoveryrate<0.05。"