"第10染色体長腕24領域に存在する一塩基多型を分析し、該分析結果に基づいて骨関節疾患を検査することを特徴とする、骨関節疾患の発症リスクおよび/または発症の有無の判定方法。前記骨関節疾患が、側弯症である、請求項1に記載の方法。前記一塩基多型が、配列番号1~3から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、請求項1または2に記載の方法。前記連鎖不平衡の関係にある塩基が、配列番号4~11から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基である、請求項1~3のいずれか1項に記載の方法。配列番号1~11から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列を有する骨関節疾患検査用プローブ。配列番号1~11から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域を増幅することのできる骨関節疾患検査用プライマー。本発明は、側弯症等の骨関節疾患の発症リスクや発症の有無を正確に検査する方法、及び該方法に用いられる検査試薬を提供することを課題とする。本発明者らは上記課題の解決のために鋭意検討した結果、第10染色体長腕24領域10q24領域に存在する一塩基多型SNPが思春期特発性側弯症AISと関連することを同定した。そして、これらの多型を調べることにより側弯症等の骨関節疾患の発症リスクや発症の推定を正確に実施できることを見出し、本発明を完成するに至った。すなわち、本発明は以下の通りである。1第10染色体長腕24領域に存在する一塩基多型を分析し、該分析結果に基づいて骨関節疾患を検査することを特徴とする、骨関節疾患の発症リスクおよび/または発症の有無の判定方法。2前記骨関節疾患が、側弯症である、1に記載の方法。3前記一塩基多型が、配列番号1~3から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、1または2に記載の方法。4前記連鎖不平衡の関係にある塩基が、配列番号4~11から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基である、1~3のいずれかに記載の方法。5配列番号1~11から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列を有する骨関節疾患検査用プローブ。6配列番号1~11から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域を増幅することのできる骨関節疾患検査用プライマー。本発明によれば、これまで予測が困難であった側弯症等の骨関節疾患の発症リスク罹患リスクを正確かつ簡便に予測することができる。また、側弯症等の骨関節疾患の発症を正確かつ簡便に判定することができる。したがって、本発明は側弯症等の骨関節疾患の予防や早期治療に貢献するものである。<1>本発明の方法本発明の方法は、ヒトの第10染色体長腕24領域10q24領域に存在するSNPを分析し、該分析結果に基づいて骨関節疾患を検査することを特徴とする、骨関節疾患の発症リスクおよび/または発症の有無の判定方法である。すなわち、本発明において、検査とは骨関節疾患の発症リスクの検査及び骨関節疾患の発症の有無の検査を含む。本発明の方法においては、SNPの分析結果を、骨関節疾患の発症リスクおよび/または発症の有無と関連付ける。骨関節疾患としては、特に限定されないが、例えば、脊椎の疾患、具体的には、側弯症が挙げられる。本発明の方法は、特に、従来病因が特定されていない特発性側弯症の検査に好適に用いられる。側弯症は、先天性、若年性、思春期性、成人性等、いずれの時期に発症するものであってもよいが、思春期性側弯症であるのが好ましく、思春期特発性側弯症AISであるのがより好ましい。本発明の方法は、いずれの人種の被検者に対しても用いることができるが、特に、日本人等のアジア人の被検者に好適に用いることができる。また、本発明の方法は、いずれの性別の被検者に対しても用いることができる。10q24領域に存在する具体的なSNPとしては、ヒトrs11190870、rs625039、およびrs11598564を挙げることができる。ここで、rs番号はNationalCenterforBiotechnologyInformationのdbSNPデータベースhttpwwwncbinlmnihgovprojectsSNPの登録番号を示す。これら3つのSNPsは、10q2431~q2432領域のLBX1遺伝子および仮想遺伝子FLJ41350が含まれる連鎖不平衡ブロックに位置する。よって、特に、当該連鎖不平衡ブロックに存在するSNPを解析することによって骨関節疾患を検査することができる。LBX1遺伝子としては、具体的には、GenBankAccessionNoNC00001010の102986733~102988717の領域の相補配列が挙げられる。仮想遺伝子FLJ41350としては、具体的には、GenBankAccessionNoNC00001010の102989351~102998616の領域が挙げられる。rs11190870はGenBankAccessionNoNT03005912の21727733番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は骨関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で骨関節疾患の可能性または発症リスクが高い。rs625039はGenBankAccessionNoNT03005912の21742175番目の塩基におけるグアニンG/アデニンAの多型を意味し、この塩基がGである場合は骨関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GA>AAの順で骨関節疾患の可能性または発症リスクが高い。rs11598564はGenBankAccessionNoNT03005912の21713130番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がGである場合は骨関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GA>AAの順で骨関節疾患の可能性または発症リスクが高い。なお、上記3個のSNPについて、SNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、配列番号1~3に示した。61番目の塩基が多型を有する。本発明においては、上記塩基に相当する塩基を解析する。上記塩基に相当する塩基とは、上記領域における該当塩基を意味する。すなわち、上記塩基に相当する塩基を解析することには、仮に人種の違いなどによって上記配列がSNP以外の位置で若干変化したとしても、上記領域における該当塩基を解析することが含まれる。また、本発明において解析する塩基は上記のものに限定されず、上記の塩基と連鎖不平衡にある塩基の多型を分析してもよい。ここで上記の塩基と連鎖不平衡にある塩基とは、上記の塩基とr2>0.5、好ましくはr2>0.8、さらに好ましくはr2>0.9の関係を満たす塩基をいう。r2は連鎖不平衡係数である。また上記の塩基と連鎖不平衡にある塩基は、例えば、HapMapデータベース等を用いて同定することができる。もしくは、複数人通常は20-40人程度から採取したDNAをシークエンサーにて配列解析し、連鎖不平衡にあるSNPを探索することにより同定することもできる。上記の塩基と連鎖不平衡にある塩基は、遺伝子型を考慮して解析した場合は、リスクアレルのホモ接合体>リスクアレルと非リスクアレルのへテロ接合体>非リスクアレルのホモ接合体の順で骨関節疾患の可能性または発症リスクが高い。rs11190870とr2>0.5で連鎖不平衡にある塩基としては、上記のrs625039およびrs11598564に加え、rs12771674、rs11598177、rs1407409、rs1322331、rs594791、rs679206、rs678741、およびrs10883597が挙げられる。これら各塩基について、rs11190870に対する連鎖不平衡係数r2、アレルの組み合わせ、およびリスクアレル等を表1に示す。表中、アレルおよびリスクアレルの表記はNCBIが公開している参照配列Build36におけるforwardstrandに対応している。また、rs12771674、rs11598177、rs1407409、rs1322331、rs594791、rs679206、rs678741、およびrs10883597について、SNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、それぞれ配列番号4~11に示した。61番目の塩基が多型を有する。Allele1がメジャーアレルである。Allele12は、humanreferencesequenceのstrandに基づいて示した。上記SNPの塩基の種類を調べ、得られた結果を上記のような基準に基づいて骨関節疾患と関連付けることにより、骨関節疾患を検査することができる。上記SNPは単独で解析されてもよいし、上記SNPの少なくとも1つを含む複数のSNPsをまとめて解析ハプロタイプ解析してもよい。例えば、上記SNPの複数をまとめて解析してもよいし、上記SNPの少なくとも1つと、骨関節疾患と関連する既知のSNPs非特許文献12等や当該既知のSNPsと連鎖不平衡にあるSNPsとを組み合わせて解析してもよい。骨関節疾患と関連する複数のSNPsをまとめて解析すれば、骨関節疾患の検査の精度が向上する。なお、いずれのSNPも、二本鎖DNAのどちらの鎖を解析してもよい。例えば、LBX1遺伝子の配列はセンス鎖を解析してもよいし、アンチセンス鎖を解析してもよい。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~3から選ばれる塩基配列の61番目の塩基を含む配列、又はその相補配列を有する10塩基以上の長さのプローブや、当該塩基と連鎖不平衡の関係にある塩基を含む配列、又はその相補配列を有する10塩基以上の長さのプローブが挙げられる。なお、当該塩基と連鎖不平衡の関係にある塩基及びその前後の領域の塩基配列は、例えば、NationalCenterforBiotechnologyInformationのdbSNPデータベースhttpwwwncbinlmnihgovprojectsSNPから取得できる。プローブの長さは好ましくは、15~35塩基であり、より好ましくは20~35塩基である。また、プライマーとしては、上記SNP部位を増幅するためのPCRに用いることのできるプライマー、又は上記SNP部位を配列解析シークエンシングするために用いることのできるプライマーが挙げられる。具体的には、配列番号1~3から選ばれる塩基配列の61番目の塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーや、当該塩基と連鎖不平衡の関係にある塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーが挙げられる。このようなプライマーの長さは10~50塩基が好ましく、15~35塩基がより好ましく、20~35塩基がさらに好ましい。上記SNP部位をシークエンシングするためのプライマーとしては、上記塩基の5’側領域、好ましくは30~100塩基上流の配列を有するプライマーや、上記塩基の3’側領域、好ましくは30~100塩基下流の領域に相補的な配列を有するプライマーが例示される。PCRによる増幅の有無で多型を判定するために用いるプライマーとしては、上記塩基を含む配列を有し、上記塩基を3’側に含むプライマーや、上記塩基を含む配列の相補配列を有し、上記塩基の相補塩基を3’側に含むプライマーなどが例示される。なお、本発明の検査用試薬はこれらのプライマーやプローブに加えて、PCR用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。以下、本発明を実施例によりさらに具体的に説明する。但し、本発明はこれらの実施例に限定されない。思春期特発性側弯症AISと関連するSNPsの同定AIS感受性を決定する遺伝的変異を同定するために、日本人被検者を用いてゲノムワイド関連解析GWASを行った。GWASとは、疾患等の表現型に関わる遺伝的変異を探索する遺伝統計学的手法である。例えば、ヒトゲノム全体を網羅するような数十万~100万ヶ所のSNPsを用いて、ある疾患の患者ケースとその疾患にかかっていない被験者コントロールとの間で、多型の頻度に差があるかどうかを統計的に検定することで、疾患と関連する遺伝的変異を見出すことができる。<被検者>本実施例に用いた被検者の特徴を表2に示す。a:ORC大阪御堂筋ロータリークラブから採用した対照被検者は健康なボランティアからなる。b:BBJから採用したGWASの対照被検者はAIS以外の13種の疾患を有する患者からなる。c:BBJから採用した再現試験の対照被検者はAIS以外の42種の疾患を有する患者からなる。AISは女性に多く発症することから、被検者としてはいずれも日本人女性を採用した。GWASに用いたAIS被検者ケースcase1050名は、8ヶ所の病院で募集した。再現試験replicationstudyに用いたAIS被検者ケースcase326名は、8ヶ所の病院で募集した。全AIS被検者は、臨床検査および放射線検査を受け、側弯症の専門医によりAISと診断された。なお、側弯の指標であるコブ角Cobbangleは全例で、20°を超えていた。GWASに用いた対照被検者コントロールcontrol1474名は、BBJに登録されたAIS以外の疾患の患者、および大阪御堂筋ロータリークラブで募集した健康なボランティアからなる。再現試験に用いた対照被検者9823名は、BBJに登録されたAIS以外の疾患の患者から採用した。本研究は東京大学医科学研究所および理化学研究所横浜研究所の倫理委員会によって承認され、全ての被検者および一部の被検者の両親からインフォームドコンセントを得た。<GWAS>AIS被検者1050名の遺伝子型を、IlluminaHuman610GenotypingBeadChipIllumina社USAを用いて解析した。対照被検者326名の遺伝子型を、IlluminaHumanHap550v3GenotypingBeadChipIllumina社USAを用いて解析した。クオリティコントロールとして、コール率callrateが0.99未満のSNPs、HardyWeinberg平衡検定でのP値がカットオフ値以下P≦106のSNPs、常染色体上にないSNPsnonautosomalSNPs、単型のSNPsmonomorphicSNPs、AIS被検者と対照被検者とで共有されていないSNPs、あいまいなコールambiguouscallのSNPs、identifybystateIBS≧1.7の近縁ペアの内の低コール率の被検者AIS被検者16名、対照被検者1名、および主成分分析PCAにより外れ値outlierであると判定された被検者AIS被検者1名を除外した。クオリティコントロールを通過したAIS被検者1033名および対照被検者1473名の455121個のSNPsについて、コクランアーミテージ傾向検定CochranArmitagetrendtestによりAISとの関連を評価した。なお、genomicinflationfactorは109であり、集団の偏りによる偽陽性の関連が得られる可能性は低いと考えられた。GWASの結果、10q24に存在する3個のSNPsrs11190870、rs625039、およびrs11598564が、ゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<1.10×107=0.05/455,121を満たし、AISとの有意な関連を示した表3のGWAS欄。当該閾値は、455,121回のテストに基づくボンフェローニ補正後のP<0.05に相当する。Allele1がメジャーアレルである。Allele12は、humanreferencesequenceのstrandに基づいて示した。Risk;リスクアレルRiskAlleleRAF;リスクアレル頻度RiskAlleleFrequencyCI;信頼区間ConfidenceIntervala:コクランアーミテージ傾向検定のP値。b:95CIでのアレルのオッズ比OddsRatio。c:Breslow-Day検定のP値。d:統合解析combinedanalysisの結果は、MantelHaenszel法により算出した。<再現試験>これら3個のSNPsとAISの関連を確認するため、GWASで解析した被検者とは独立の被検者を用いて再現試験replicationstudyを行った。AIS被検者326名の遺伝子型は、multiplexPCRinvaderassayThirdWaveTechnologiesにより解析した。対照被検者9823名の遺伝子型は、IlluminaHumanHap550v3GenotypingBeadChipIllumina社USAを用いて解析した。再現試験の結果、3個のSNPs全てが、ボンフェローニ補正後の有意性の閾値であるP<1.67×102=0.05/3を満たし、AISとの有意な関連を示した表3のReplication欄。また、Breslow-Day検定のP値はいずれも0.05以下であり表3のPhet欄、GWASの被検者と再現試験の被検者とでは有意な不均一性significantheterogeneityは認められなかった。GWASおよび再現試験のいずれにおいても、AISとの最も強い関連が認められたSNPsはrs11190870であった。MantelHaenszel法を用いてGWASと再現試験の統合解析combinedanalysisを行った結果、rs11190870はP=1.24×1019オッズ比=1.56でAISと関連した表3のCombined欄。10q24領域の連鎖不平衡マップおよび関連解析の結果を図1に示す。rs11190870、rs625039、およびrs11598564はいずれも、10q2431~q2432領域のLBX1遺伝子および仮想遺伝子FLJ41350が含まれる約80kbの連鎖不平衡ブロックに位置しており、より具体的には、rs11190870はLBX1遺伝子の3’フランキング領域に位置していた図1。なお、これら3個のSNPsによるハプロタイプ解析を行ったところ、rs11190870より強くAISと関連するハプロタイプは見出されなかった。これまでにゲノムワイド連鎖解析によりAISの病因となる遺伝子座が見出されているが、10q24領域は報告されていなかった非特許文献1、2。また、TDT法に基づくGWASにおいては、上記3個のSNPsrs11190870、rs625039、およびrs11598564の内、rs11598564がAISとの関連の上位100位に含まれていたが、TDT法におけるP=8.19×105全人種またはP=1.03×103非ヒスパニックであり有意な関連は認められなかった非特許文献12。いずれの知見においてもゲノムワイド水準でAISとの関連が認められた遺伝子や一塩基多型SNPsは報告されておらず、すなわち本願は、ゲノムワイド水準で有意にAISと関連するSNPsを初めて見出したものである。LBX1遺伝子は、ホメオボックス遺伝子としてクローニングされ、脊椎動物においては、脊髄後索、後脳、心臓神経堤の一部、筋原細胞等での発現が報告されている。本発明者も脊髄および骨格筋におけるLBX1遺伝子の発現を確認した。マウスにおいて、Lbx1タンパク質は、脊髄後索ニューロンや、後脳における体性感覚ニューロンのidentityの決定因子として機能することが知られている。脊髄後索は体性感覚情報の伝達に関与し、AIS等の側弯症患者ではしばしば体性感覚機能の異常が認められるGuoXetalSpine31E437E440。また、痙攣性両側麻痺SpasticDiplegiaの患者においては、後根神経切除dorsalrhizotomy後に側弯症の罹患率が増大することが知られているSteinbokPetalJNeurosurg102363373。以上より、LBX1遺伝子は、体性感覚機能の異常を介して側弯症の病因となっている可能性がある。<補完imputation>また、10q24領域とAISの関連をさらに解析するため、rs11190870が含まれる上記連鎖不平衡ブロックに存在する他のSNPsGWASでは解析されなかったものの遺伝子型をMACHhttpwwwsphumicheducsgylimachindexhtmlにより解析し、AISとの関連を評価した。補完imputationの結果、新たに3個のSNPsが、ゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<1.10×107を満たし、AISと有意に関連することが見出された図1の下パネルの白い四角。これら3個のSNPsはLBX1遺伝子の3’フランキング領域に位置し、いずれもrs11190870に対しr2>0.7で連鎖不平衡の関係にあった。なお、これら3個のSNPsとAISとの関連は、rs11190870とAISとの関連よりも弱いものであった。<男性被検者を用いた関連解析>女性被検者でAISとの強い関連が認められた3個のSNPsrs11190870、rs625039、およびrs11598564について、男性被検者におけるAISとの関連を検証した。94名の男性AIS被検者、および、1849名の男性対照被検者は、女性被検者と同様の基準で採用した。AIS被検者の遺伝子型は、multiplexPCRinvaderassayThirdWaveTechnologiesにより解析した。対照被検者の遺伝子型は、IlluminaHumanHap550v3GenotypingBeadChipIllumina社USAを用いて解析した。上記3個のSNPsについて、コクランアーミテージ傾向検定CochranArmitagetrendtestによりAISとの関連を評価した。その結果、3個のSNPs全てについて、AISとの関連が再現された表4。オッズ比でみると、rs11190870とAISとの関連は女性の集団での関連よりも強く、また、rs625039およびrs11598564とAISとの関連は女性の集団での関連と同等であった。以上より、rs11190870、rs625039、およびrs11598564は、男性においてもAIS感受性多型であり、男性におけるAISの検査にも用いることができると明らかになった。RAF;リスクアレル頻度RiskAlleleFrequencyCI;信頼区間ConfidenceIntervala:コクランアーミテージ傾向検定のP値。b:95CIでのアレルのオッズ比OddsRatio。"