(57)【特許請求の範囲】
【請求項１】
第１０染色体長腕２４領域に存在する、配列番号１〜３，６，７，９〜１１から選ばれる塩基配列の塩基番号６１番目の塩基に相当する一塩基多型を分析し、該分析結果に基づいて、
配列番号１の６１番目の塩基に相当する塩基がＴの場合；
配列番号２の６１番目の塩基に相当する塩基がＧの場合；  
配列番号３の６１番目の塩基に相当する塩基がＧの場合；
配列番号６の６１番目の塩基に相当する塩基がＴの場合；
配列番号７の６１番目の塩基に相当する塩基がＡの場合；
配列番号９の６１番目の塩基に相当する塩基がＴの場合；
配列番号１０の６１番目の塩基に相当する塩基がＡの場合；または、
配列番号１１の６１番目の塩基に相当する塩基がＣの場合、
側弯症の発症リスクが高い、または、側弯症の可能性が高いと判定されることにより側弯症を検査することを特徴とする、側弯症の発症リスクおよび／または発症の有無の判定方法。
【請求項２】
配列番号１〜３，６，７，９〜１１から選ばれる塩基配列において、塩基番号６１番目の塩基を含む１５塩基以上の配列、又はその相補配列を有する側弯症検査用プローブ。
【請求項３】
配列番号１〜３，６，７，９〜１１から選ばれる塩基配列において、塩基番号６１番目の塩基を含む領域を増幅することのできる側弯症検査用プライマー。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は側弯症等の骨・関節疾患の発症リスクや発症の有無を判定するための検査方法及び該検査方法に用いられる試薬に関する。
【背景技術】
【０００２】
側弯症は、脊椎が左右方向に弯曲した病態をいう。側弯症の内、その病因が明らかでないものを特発性側弯症と呼び、側弯症の８０〜９０％を占める。
【０００３】
側弯症は、通学年代の児童に多く発症し、特に、女子に多く発症する。１０歳から骨格が成熟するまでの児童において、側弯の指標であるコブ角（cobb angle）で少なくとも１０°の弯曲を示し、且つ、その病因が明らかでないものを思春期特発性側弯症（Adolescent idiopathic scoliosis；ＡＩＳ）と定義する。ＡＩＳの発症頻度は、通学年代の児童
の内、日本で２％、世界では２〜３％であり、日本では毎年１万人程度が発症している。
【０００４】
側弯症の診断は主にＸ線検査により行われるが、Ｘ線検査は発症前や発症初期における側弯症の診断には有効でない。また、特発性側弯症の治療は、対症療法により行われるのみで、病因が明らかでない以上、原因療法は行われていない。よって、側弯症の発症前診断（リスク診断）や早期診断を可能にし、また、原因治療を可能にするため、側弯症と関連する遺伝子や一塩基多型（ＳＮＰｓ）の同定が望まれている。
【０００５】
ゲノムワイド連鎖解析（genome-wide linkage analysis）によりＡＩＳの病因となる多くの遺伝子座が見出されており、ＡＩＳは複雑な遺伝的素因に基づくと考えられている（非特許文献１、２）。また、候補遺伝子解析によりＡＩＳ感受性遺伝子が報告されているが（非特許文献３〜９）、いずれの遺伝子についても別の人種を被検者とした場合にＡＩＳとの関連の再現性が得られていない（非特許文献１０、１１）。
【０００６】
さらに、近年、伝達不平衡解析（Transmission Disequilibrium Test；ＴＤＴ）法に基づくゲノムワイド関連解析（Genome-wide association study；ＧＷＡＳ）によって、Ａ
ＩＳと関連する遺伝子の候補が報告された（非特許文献１２）。しかしながら、それらは、多重検定補正後のケース−コントロール関連解析においてはＡＩＳとの関連の再現性が得られていない。
【０００７】
以上の通り、側弯症感受性遺伝子に関する知見は存在するものの、ゲノムワイド水準で側弯症との関連が認められた遺伝子や一塩基多型（ＳＮＰｓ）は報告されていない。
【０００８】
第１０染色体長腕２４領域（10q24）に存在するLBX1遺伝子は、ホメオボックス遺伝子
としてクローニングされており、脊髄後索、後脳、心臓神経堤の一部、筋原細胞等での発現が報告されている。しかしながら、LBX1遺伝子と側弯症との関連は知られていない。
【先行技術文献】
【非特許文献】
【０００９】
【非特許文献１】Wise, C.A. et al. Curr. Genomics 9, 51-59 (2008).
【非特許文献２】Raggio, C. L. et al. J. Orthop. Res. 27, 1366-1372 (2009).
【非特許文献３】Wu, J. et al. Spine 31, 1131-1136 (2006).
【非特許文献４】Zhang, H.Q. et al. Spine 34, 760-764 (2009).
【非特許文献５】Chen, Z. et al. Eur. J. Hum. Genet. 17, 525-532 (2009).
【非特許文献６】Qiu, X. S. et al. Spine 32, 1748-1753 (2007).
【非特許文献７】Wang, H. et al. Spine 33, 2199-2203 (2008).
【非特許文献８】Inoue, M. et al. Stud. Health Technol. Inform. 91, 90-96 (2002).
【非特許文献９】Yeung, H.Y. et al. Stud. Health Technol. Inform. 123, 18-24 (2006).
【非特許文献１０】Takahashi, Y. et al. J. Orthop. Res. 29, 834-837(2011).
【非特許文献１１】Takahashi, Y. et al. J. Orthop. Res. (2011).
【非特許文献１２】Sharma, S. et al. Hum Mol Genet. 20, 1456-1466 (2011).
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本発明は、側弯症等の骨・関節疾患の発症リスクや発症の有無を正確に検査する方法、及び該方法に用いられる検査試薬を提供することを課題とする。
【課題を解決するための手段】
【００１１】
本発明者らは上記課題の解決のために鋭意検討した結果、第１０染色体長腕２４領域（10q24領域）に存在する一塩基多型（ＳＮＰ）が思春期特発性側弯症（ＡＩＳ）と関連す
ることを同定した。そして、これらの多型を調べることにより側弯症等の骨・関節疾患の発症リスクや発症の推定を正確に実施できることを見出し、本発明を完成するに至った。
【００１２】
すなわち、本発明は以下の通りである。
[１]
第１０染色体長腕２４領域に存在する一塩基多型を分析し、該分析結果に基づいて骨・関節疾患を検査することを特徴とする、骨・関節疾患の発症リスクおよび／または発症の有無の判定方法。
[２]
前記骨・関節疾患が、側弯症である、[１]に記載の方法。
[３]
前記一塩基多型が、配列番号１〜３から選ばれる塩基配列の塩基番号６１番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、[１]または[２]に記載の方法。
[４]
前記連鎖不平衡の関係にある塩基が、配列番号４〜１１から選ばれる塩基配列の塩基番号６１番目の塩基に相当する塩基である、[１]〜[３]のいずれかに記載の方法。
[５]
配列番号１〜１１から選ばれる塩基配列において、塩基番号６１番目の塩基を含む１０塩基以上の配列、又はその相補配列を有する骨・関節疾患検査用プローブ。
[６]
配列番号１〜１１から選ばれる塩基配列において、塩基番号６１番目の塩基を含む領域を増幅することのできる骨・関節疾患検査用プライマー。
【発明の効果】
【００１３】
本発明によれば、これまで予測が困難であった側弯症等の骨・関節疾患の発症リスク（罹患リスク）を正確かつ簡便に予測することができる。また、側弯症等の骨・関節疾患の発症を正確かつ簡便に判定することができる。したがって、本発明は側弯症等の骨・関節疾患の予防や早期治療に貢献するものである。
【図面の簡単な説明】
【００１４】
【図１】10q24領域の連鎖不平衡マップ（上パネル）およびＡＩＳとの関連解析の結果（下パネル）を示す図。連鎖不平衡マップは、Phase II HapMap (release24) JPTに基づいて作成し、ＧＷＡＳで遺伝子型を解析したＳＮＰｓをボックスで囲んだ。関連解析の結果は、コクラン・アーミテージ傾向検定のＰ値の-log10として示す。関連解析の結果において、黒い菱形はＧＷＡＳの結果を、白い四角は補完（imputation）の結果を示す。
【発明を実施するための形態】
【００１５】
＜１＞本発明の方法
本発明の方法は、ヒトの第１０染色体長腕２４領域（10q24領域）に存在するＳＮＰを
分析し、該分析結果に基づいて骨・関節疾患を検査することを特徴とする、骨・関節疾患の発症リスクおよび／または発症の有無の判定方法である。すなわち、本発明において、「検査」とは骨・関節疾患の発症リスクの検査及び骨・関節疾患の発症の有無の検査を含む。本発明の方法においては、ＳＮＰの分析結果を、骨・関節疾患の発症リスクおよび／または発症の有無と関連付ける。
【００１６】
骨・関節疾患としては、特に限定されないが、例えば、脊椎の疾患、具体的には、側弯症が挙げられる。本発明の方法は、特に、従来病因が特定されていない特発性側弯症の検査に好適に用いられる。側弯症は、先天性、若年性、思春期性、成人性等、いずれの時期に発症するものであってもよいが、思春期性側弯症であるのが好ましく、思春期特発性側弯症（ＡＩＳ）であるのがより好ましい。
【００１７】
本発明の方法は、いずれの人種の被検者に対しても用いることができるが、特に、日本人等のアジア人の被検者に好適に用いることができる。また、本発明の方法は、いずれの性別の被検者に対しても用いることができる。
【００１８】
10q24領域に存在する具体的なＳＮＰとしては、ヒトrs11190870、rs625039、およびrs11598564を挙げることができる。ここで、rs番号はNational Center for Biotechnology InformationのdbSNPデータベース（http//www.ncbi.nlm.nih.gov/projects/SNP/）の登録
番号を示す。これら３つのＳＮＰｓは、10q24.31〜q24.32領域のLBX1遺伝子および仮想遺伝子FLJ41350が含まれる連鎖不平衡ブロックに位置する。よって、特に、当該連鎖不平衡ブロックに存在するＳＮＰを解析することによって骨・関節疾患を検査することができる。LBX1遺伝子としては、具体的には、GenBank Accession No. NC_000010.10の102986733
〜102988717の領域の相補配列が挙げられる。仮想遺伝子FLJ41350としては、具体的には
、GenBank Accession No. NC_000010.10の102989351〜102998616の領域が挙げられる。
【００１９】
rs11190870はGenBank Accession No. NT_030059.12の21727733番目の塩基におけるチミン（Ｔ）／シトシン（Ｃ）の多型を意味し、この塩基がＴである場合は骨・関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、ＴＴ＞ＴＣ＞ＣＣの順で骨・関節疾患の可能性または発症リスクが高い。
【００２０】
rs625039はGenBank Accession No. NT_030059.12の21742175番目の塩基におけるグアニン（Ｇ）／アデニン（Ａ）の多型を意味し、この塩基がＧである場合は骨・関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、ＧＧ＞ＧＡ＞ＡＡの順で骨・関節疾患の可能性または発症リスクが高い。
【００２１】
rs11598564はGenBank Accession No. NT_030059.12の21713130番目の塩基におけるアデニン（Ａ）／グアニン（Ｇ）の多型を意味し、この塩基がＧである場合は骨・関節疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、ＧＧ＞ＧＡ＞ＡＡの順で骨・関節疾患の可能性または発症リスクが高い。
【００２２】
なお、上記３個のＳＮＰについて、ＳＮＰ塩基及びその前後６０ｂｐの領域を含む合計１２１ｂｐの長さの配列を、配列番号１〜３に示した。６１番目の塩基が多型を有する。
【００２３】
本発明においては、上記塩基に相当する塩基を解析する。「上記塩基に相当する塩基」とは、上記領域における該当塩基を意味する。すなわち、「上記塩基に相当する塩基を解析する」ことには、仮に人種の違いなどによって上記配列がＳＮＰ以外の位置で若干変化したとしても、上記領域における該当塩基を解析することが含まれる。
【００２４】
また、本発明において解析する塩基は上記のものに限定されず、上記の塩基と連鎖不平衡にある塩基の多型を分析してもよい。ここで「上記の塩基と連鎖不平衡にある塩基」とは、上記の塩基とｒ2＞０．５、好ましくはｒ2＞０．８、さらに好ましくはｒ2＞０．９
の関係を満たす塩基をいう。ｒ2は連鎖不平衡係数である。また上記の塩基と連鎖不平衡
にある塩基は、例えば、HapMapデータベース(http://www.hapmap.org/index.html.ja)等
を用いて同定することができる。もしくは、複数人（通常は２０−４０人程度）から採取したＤＮＡをシークエンサーにて配列解析し、連鎖不平衡にあるＳＮＰを探索することにより同定することもできる。上記の塩基と連鎖不平衡にある塩基は、遺伝子型を考慮して解析した場合は、リスクアレルのホモ接合体  ＞  リスクアレルと非リスクアレルのへテロ接合体  ＞  非リスクアレルのホモ接合体の順で骨・関節疾患の可能性または発症リスクが高い。
【００２５】
rs11190870とｒ2＞０．５で連鎖不平衡にある塩基としては、上記のrs625039およびrs11598564に加え、rs12771674、rs11598177、rs1407409、rs1322331、rs594791、rs679206
、rs678741、およびrs10883597が挙げられる。これら各塩基について、rs11190870に対する連鎖不平衡係数（ｒ2）、アレルの組み合わせ、およびリスクアレル等を表１に示す。
表中、アレルおよびリスクアレルの表記はＮＣＢＩが公開している参照配列（Build 36）におけるforward strandに対応している。また、rs12771674、rs11598177、rs1407409、rs1322331、rs594791、rs679206、rs678741、およびrs10883597について、ＳＮＰ塩基及びその前後６０ｂｐの領域を含む合計１２１ｂｐの長さの配列を、それぞれ配列番号４〜１１に示した。６１番目の塩基が多型を有する。
【００２６】
【表１】
0005895317.tif 000002
Allele 1がメジャーアレルである。
Allele 1/2は、human reference sequenceの(+) strandに基づいて示した。
【００２７】
上記ＳＮＰの塩基の種類を調べ、得られた結果を上記のような基準に基づいて骨・関節疾患と関連付けることにより、骨・関節疾患を検査することができる。上記ＳＮＰは単独で解析されてもよいし、上記ＳＮＰの少なくとも１つを含む複数のＳＮＰｓをまとめて解析（ハプロタイプ解析）してもよい。例えば、上記ＳＮＰの複数をまとめて解析してもよいし、上記ＳＮＰの少なくとも１つと、骨・関節疾患と関連する既知のＳＮＰｓ（非特許文献１２等）や当該既知のＳＮＰｓと連鎖不平衡にあるＳＮＰｓとを組み合わせて解析してもよい。骨・関節疾患と関連する複数のＳＮＰｓをまとめて解析すれば、骨・関節疾患の検査の精度が向上する。なお、いずれのＳＮＰも、二本鎖ＤＮＡのどちらの鎖を解析してもよい。例えば、ＬＢＸ１遺伝子の配列はセンス鎖を解析してもよいし、アンチセンス鎖を解析してもよい。
【００２８】
ＳＮＰの解析に用いる試料としては、染色体ＤＮＡを含む試料であれば特に制限されないが、例えば、血液、尿等の体液サンプル、口腔粘膜などの細胞、毛髪等の体毛などが挙げられる。ＳＮＰの解析にはこれらの試料を直接使用することもできるが、これらの試料から染色体ＤＮＡを常法により単離し、これを用いて解析することが好ましい。
【００２９】
ＳＮＰの解析は、通常の遺伝子多型解析方法によって行うことができる。例えば、シークエンス解析、ＰＣＲ、ハイブリダイゼーション、インベーダー法などが挙げられるが、これらに限定されない。
【００３０】
シークエンス解析は通常の方法により行うことができる。具体的には、多型を示す塩基の５’側  数十塩基の位置に設定したプライマーを使用してシークエンス反応を行い、その解析結果から、該当する位置がどの種類の塩基であるかを決定することができる。なお、シークエンス反応の前に、あらかじめＳＮＰ部位を含む断片をＰＣＲなどによって増幅しておくことが好ましい。
【００３１】
また、ＳＮＰの解析は、ＰＣＲによる増幅の有無を調べることによって行うことができる。例えば、多型を示す塩基を含む領域に対応する配列を有し、かつ、３’末端が各多型に対応するプライマーをそれぞれ用意する。それぞれのプライマーを使用してＰＣＲを行い、増幅産物の有無によってどのタイプの多型であるかを決定することができる。また、ＬＡＭＰ法（特許第３３１３３５８号明細書）、ＮＡＳＢＡ法（Nucleic Acid Sequence-Based Amplification；特許２８４３５８６号明細書）、ＩＣＡＮ法（特開２００２−２
３３３７９号公報）などによって増幅の有無を調べることもできる。その他、単鎖増幅法を用いてもよい。
【００３２】
また、ＳＮＰ部位を含むＤＮＡ断片を増幅し、増幅産物の電気泳動における移動度の違いによってどのタイプの多型であるかを決定することもできる。このような方法としては、例えば、ＰＣＲ−ＳＳＣＰ（single-strand conformation polymorphism）法（Genomics. 1992 Jan 1; 12(1): 139-146.）が挙げられる。具体的には、まず、目的のＳＮＰを含むＤＮＡを増幅し、増幅したＤＮＡを一本鎖ＤＮＡに解離させる。次いで、解離させた一本鎖ＤＮＡを非変性ゲル上で分離し、分離した一本鎖ＤＮＡのゲル上での移動度の違いによってどのタイプの多型であるかを決定することができる。
【００３３】
さらに、多型を示す塩基が制限酵素認識配列に含まれる場合は、制限酵素による切断の有無によって解析することもできる（ＲＦＬＰ法）。この場合、まず、ＤＮＡ試料を制限酵素により切断する。次いで、ＤＮＡ断片を分離し、検出されたＤＮＡ断片の大きさによってどのタイプの多型であるかを決定することができる。
【００３４】
また、ハイブリダイゼーションの有無を調べることによって多型の種類を解析することも可能である。すなわち、各塩基に対応するプローブを用意し、いずれのプローブにハイブリダイズするかを調べることによってＳＮＰがいずれの塩基であるかを調べることもできる。
【００３５】
このようにしてＳＮＰがいずれの塩基であるかを決定することで、骨・関節疾患を検査するためのデータを得ることができる。
【００３６】
＜２＞本発明の検査用試薬
本発明はまた、側弯症等の骨・関節疾患を検査するためのプライマーやプローブなどの検査試薬を提供する。このようなプローブとしては、上記ＳＮＰ部位を含み、ハイブリダイズの有無によってＳＮＰ部位の塩基の種類を判定できるプローブが挙げられる。具体的には、配列番号１〜３から選ばれる塩基配列の６１番目の塩基を含む配列、又はその相補配列を有する１０塩基以上の長さのプローブや、当該塩基と連鎖不平衡の関係にある塩基を含む配列、又はその相補配列を有する１０塩基以上の長さのプローブが挙げられる。なお、「当該塩基と連鎖不平衡の関係にある塩基」及びその前後の領域の塩基配列は、例えば、National Center for Biotechnology InformationのdbSNPデータベース（http//www.ncbi.nlm.nih.gov/projects/SNP/）から取得できる。プローブの長さは好ましくは、１５〜３５塩基であり、より好ましくは２０〜３５塩基である。
【００３７】
また、プライマーとしては、上記ＳＮＰ部位を増幅するためのＰＣＲに用いることのできるプライマー、又は上記ＳＮＰ部位を配列解析（シークエンシング）するために用いることのできるプライマーが挙げられる。具体的には、配列番号１〜３から選ばれる塩基配列の６１番目の塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーや、当該塩基と連鎖不平衡の関係にある塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーが挙げられる。このようなプライマーの長さは１０〜５０塩基が好ましく、１５〜３５塩基がより好ましく、２０〜３５塩基がさらに好ましい。
【００３８】
上記ＳＮＰ部位をシークエンシングするためのプライマーとしては、上記塩基の５’側領域、好ましくは３０〜１００塩基上流の配列を有するプライマーや、上記塩基の３’側領域、好ましくは３０〜１００塩基下流の領域に相補的な配列を有するプライマーが例示される。ＰＣＲによる増幅の有無で多型を判定するために用いるプライマーとしては、上記塩基を含む配列を有し、上記塩基を３’側に含むプライマーや、上記塩基を含む配列の相補配列を有し、上記塩基の相補塩基を３’側に含むプライマーなどが例示される。
【００３９】
なお、本発明の検査用試薬はこれらのプライマーやプローブに加えて、ＰＣＲ用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。
【実施例】
【００４０】
以下、本発明を実施例によりさらに具体的に説明する。但し、本発明はこれらの実施例に限定されない。
【００４１】
[思春期特発性側弯症（ＡＩＳ）と関連するＳＮＰｓの同定]
ＡＩＳ感受性を決定する遺伝的変異を同定するために、日本人被検者を用いてゲノムワイド関連解析（ＧＷＡＳ）を行った。ＧＷＡＳとは、疾患等の表現型に関わる遺伝的変異を探索する遺伝統計学的手法である。例えば、ヒトゲノム全体を網羅するような数十万〜１００万ヶ所のＳＮＰｓを用いて、ある疾患の患者（ケース）とその疾患にかかっていない被験者（コントロール）との間で、多型の頻度に差があるかどうかを統計的に検定する
ことで、疾患と関連する遺伝的変異を見出すことができる。
【００４２】
＜被検者＞
本実施例に用いた被検者の特徴を表２に示す。
【００４３】
【表２】
0005895317.tif 000003
ａ：ＯＲＣ（大阪御堂筋ロータリークラブ）から採用した対照被検者は健康なボランティアからなる。
ｂ：ＢＢＪから採用したＧＷＡＳの対照被検者はＡＩＳ以外の１３種の疾患を有する患者からなる。
ｃ：ＢＢＪから採用した再現試験の対照被検者はＡＩＳ以外の４２種の疾患を有する患者からなる。
【００４４】
ＡＩＳは女性に多く発症することから、被検者としてはいずれも日本人女性を採用した。ＧＷＡＳに用いたＡＩＳ被検者（ケース（case））1,050名は、８ヶ所の病院で募集し
た。再現試験（replication study）に用いたＡＩＳ被検者（ケース（case））326名は、８ヶ所の病院で募集した。全ＡＩＳ被検者は、臨床検査および放射線検査を受け、側弯症の専門医によりＡＩＳと診断された。なお、側弯の指標であるコブ角（Cobb angle）は全例で、２０°を超えていた。ＧＷＡＳに用いた対照被検者（コントロール（control））1,474名は、BBJに登録されたＡＩＳ以外の疾患の患者、および大阪御堂筋ロータリークラ
ブで募集した健康なボランティアからなる。再現試験に用いた対照被検者9,823名は、BBJに登録されたＡＩＳ以外の疾患の患者から採用した。
【００４５】
本研究は東京大学医科学研究所および理化学研究所横浜研究所の倫理委員会によって承認され、全ての被検者および一部の被検者の両親からインフォームドコンセントを得た。
【００４６】
＜ＧＷＡＳ＞
ＡＩＳ被検者1,050名の遺伝子型を、Illumina Human610 Genotyping BeadChip（Illumina社（USA））を用いて解析した。対照被検者326名の遺伝子型を、Illumina HumanHap550v3 Genotyping BeadChip（Illumina社（USA））を用いて解析した。クオリティコントロ
ールとして、コール率（call rate）が０．９９未満のＳＮＰｓ、Hardy-Weinberg平衡検
定でのＰ値がカットオフ値以下（Ｐ≦１０-6）のＳＮＰｓ、常染色体上にないＳＮＰｓ（non-autosomal SNPs）、単型のＳＮＰｓ（monomorphic SNPs）、ＡＩＳ被検者と対照被検者とで共有されていないＳＮＰｓ、あいまいなコール（ambiguous call）のＳＮＰｓ、identify-by-state（ＩＢＳ）≧１．７の近縁ペアの内の低コール率の被検者（ＡＩＳ被検
者１６名、対照被検者１名）、および主成分分析（ＰＣＡ）により外れ値（outlier）で
あると判定された被検者（ＡＩＳ被検者１名）を除外した。クオリティコントロールを通過したＡＩＳ被検者1,033名および対照被検者1,473名の455,121個のＳＮＰｓについて、
コクラン・アーミテージ傾向検定（Cochran-Armitage trend test）によりＡＩＳとの関
連を評価した。なお、genomic inflation factorは1.09であり、集団の偏りによる偽陽性の関連が得られる可能性は低いと考えられた。
【００４７】
ＧＷＡＳの結果、10q24に存在する３個のＳＮＰｓ（rs11190870、rs625039、およびrs11598564）が、ゲノムワイドな有意性（genome-wide significance）の閾値として設定さ
れたＰ＜１．１０×１０-7（＝０．０５／４５５，１２１）を満たし、ＡＩＳとの有意な関連を示した（表３の「ＧＷＡＳ」欄）。当該閾値は、４５５，１２１回のテストに基づくボンフェローニ補正後のＰ＜０．０５に相当する。
【００４８】
【表３】
0005895317.tif 000004
Allele 1がメジャーアレルである。
Allele 1/2は、human reference sequenceの(+) strandに基づいて示した。
Ｒｉｓｋ；リスクアレル（Risk Allele）
ＲＡＦ；リスクアレル頻度（Risk Allele Frequency）
ＣＩ；信頼区間（Confidence Interval）
ａ：コクラン・アーミテージ傾向検定のＰ値。
ｂ：９５％ＣＩでのアレルのオッズ比（Odds Ratio）。
ｃ：Ｂｒｅｓｌｏｗ−Ｄａｙ検定のＰ値。
ｄ：統合解析（combined analysis）の結果は、Mantel-Haenszel法により算出した。
【００４９】
＜再現試験＞
これら３個のＳＮＰｓとＡＩＳの関連を確認するため、ＧＷＡＳで解析した被検者とは独立の被検者を用いて再現試験（replication study）を行った。ＡＩＳ被検者326名の遺伝子型は、multiplex-PCR invader assay（Third Wave Technologies）により解析した。対照被検者9,823名の遺伝子型は、Illumina HumanHap550v3 Genotyping BeadChip（Illumina社（USA））を用いて解析した。
【００５０】
再現試験の結果、３個のＳＮＰｓ全てが、ボンフェローニ補正後の有意性の閾値であるＰ＜１．６７×１０-2（＝０．０５／３）を満たし、ＡＩＳとの有意な関連を示した（表３の「Replication」欄）。また、Ｂｒｅｓｌｏｗ−Ｄａｙ検定のＰ値はいずれも０．０
５以下であり（表３の「Phet」欄）、ＧＷＡＳの被検者と再現試験の被検者とでは有意な
不均一性（significant heterogeneity）は認められなかった。
【００５１】
ＧＷＡＳおよび再現試験のいずれにおいても、ＡＩＳとの最も強い関連が認められたＳＮＰｓはrs11190870であった。Mantel-Haenszel法を用いてＧＷＡＳと再現試験の統合解
析（combined analysis）を行った結果、rs11190870はＰ＝１．２４×１０-19（オッズ比＝１．５６）でＡＩＳと関連した（表３の「Combined」欄）。
【００５２】
10q24領域の連鎖不平衡マップおよび関連解析の結果を図１に示す。rs11190870、rs625039、およびrs11598564はいずれも、10q24.31〜q24.32領域のLBX1遺伝子および仮想遺伝
子FLJ41350が含まれる約８０ｋｂの連鎖不平衡ブロックに位置しており、より具体的には、rs11190870はLBX1遺伝子の３’フランキング領域に位置していた（図１）。なお、これら３個のＳＮＰｓによるハプロタイプ解析を行ったところ、rs11190870より強くＡＩＳと関連するハプロタイプは見出されなかった。
【００５３】
これまでにゲノムワイド連鎖解析によりＡＩＳの病因となる遺伝子座が見出されているが、10q24領域は報告されていなかった（非特許文献１、２）。また、ＴＤＴ法に基づく
ＧＷＡＳにおいては、上記３個のＳＮＰｓ（rs11190870、rs625039、およびrs11598564）の内、rs11598564がＡＩＳとの関連の上位１００位に含まれていたが、ＴＤＴ法におけるＰ＝８．１９×１０-5（全人種）またはＰ＝１．０３×１０-3（非ヒスパニック）であり有意な関連は認められなかった（非特許文献１２）。いずれの知見においてもゲノムワイド水準でＡＩＳとの関連が認められた遺伝子や一塩基多型（ＳＮＰｓ）は報告されておらず、すなわち本願は、ゲノムワイド水準で有意にＡＩＳと関連するＳＮＰｓを初めて見出したものである。
【００５４】
LBX1遺伝子は、ホメオボックス遺伝子としてクローニングされ、脊椎動物においては、脊髄後索、後脳、心臓神経堤の一部、筋原細胞等での発現が報告されている。本発明者も脊髄および骨格筋におけるLBX1遺伝子の発現を確認した。マウスにおいて、Lbx1タンパク質は、脊髄後索ニューロンや、後脳における体性感覚ニューロンのidentityの決定因子として機能することが知られている。脊髄後索は体性感覚情報の伝達に関与し、ＡＩＳ等の側弯症患者ではしばしば体性感覚機能の異常が認められる（Guo, X. et al. Spine 31, E437-E440 (2006)）。また、痙攣性両側麻痺（Spastic Diplegia）の患者においては、後
根神経切除（dorsal rhizotomy）後に側弯症の罹患率が増大することが知られている（Steinbok, P. et al. J. Neurosurg. 102, 363-373. (2005)）。以上より、LBX1遺伝子は、体性感覚機能の異常を介して側弯症の病因となっている可能性がある。
【００５５】
＜補完（imputation）＞
また、10q24領域とＡＩＳの関連をさらに解析するため、rs11190870が含まれる上記連
鎖不平衡ブロックに存在する他のＳＮＰｓ（ＧＷＡＳでは解析されなかったもの）の遺伝子型をＭＡＣＨ（http://www.sph.umich.edu/csg/yli/mach/index.html）により解析し、ＡＩＳとの関連を評価した。補完（imputation）の結果、新たに３個のＳＮＰｓが、ゲノムワイドな有意性（genome-wide significance）の閾値として設定されたＰ＜１．１０×１０-7を満たし、ＡＩＳと有意に関連することが見出された（図１の下パネルの白い四角）。これら３個のＳＮＰｓはLBX1遺伝子の３’フランキング領域に位置し、いずれもrs11190870に対しｒ2＞０．７で連鎖不平衡の関係にあった。なお、これら３個のＳＮＰｓと
ＡＩＳとの関連は、rs11190870とＡＩＳとの関連よりも弱いものであった。
【００５６】
＜男性被検者を用いた関連解析＞
女性被検者でＡＩＳとの強い関連が認められた３個のＳＮＰｓ（rs11190870、rs625039、およびrs11598564）について、男性被検者におけるＡＩＳとの関連を検証した。９４名の男性ＡＩＳ被検者、および、１８４９名の男性対照被検者は、女性被検者と同様の基準
で採用した。ＡＩＳ被検者の遺伝子型は、multiplex-PCR invader assay（Third Wave Technologies）により解析した。対照被検者の遺伝子型は、Illumina HumanHap550v3 Genotyping BeadChip（Illumina社（USA））を用いて解析した。上記３個のＳＮＰｓについて
、コクラン・アーミテージ傾向検定（Cochran-Armitage trend test）によりＡＩＳとの
関連を評価した。その結果、３個のＳＮＰｓ全てについて、ＡＩＳとの関連が再現された（表４）。オッズ比でみると、rs11190870とＡＩＳとの関連は女性の集団での関連よりも強く、また、rs625039およびrs11598564とＡＩＳとの関連は女性の集団での関連と同等であった。以上より、rs11190870、rs625039、およびrs11598564は、男性においてもＡＩＳ感受性多型であり、男性におけるＡＩＳの検査にも用いることができると明らかになった。
【００５７】
【表４】
0005895317.tif 000005
ＲＡＦ；リスクアレル頻度（Risk Allele Frequency）
ＣＩ；信頼区間（Confidence Interval）
ａ：コクラン・アーミテージ傾向検定のＰ値。
ｂ：９５％ＣＩでのアレルのオッズ比（Odds Ratio）。
【産業上の利用可能性】
【００５８】
以上の通り、全ゲノムレベルでのケース−コントロール（Case-Control）関連解析により、ゲノムワイド水準を満たしてＡＩＳと関連する領域が見出された。当該領域に存在するＳＮＰｓは側弯症等の骨・関節疾患の検査に有用であり、側弯症等の骨・関節疾患の予防および／または治療に貢献するものである。
【図１】
0005895317.tif 000006
【配列表】
0005895317000001.app