以下の1~5のいずれかの領域に存在する一塩基多型を分析し、該分析結果に基づいて免疫疾患または閉塞性肺疾患を検査することを特徴とする、免疫疾患または閉塞性肺疾患の発症リスクおよび/または発症の有無の判定方法。1第4染色体長腕31領域2第12染色体長腕13領域3第10染色体短腕14領域4第5染色体長腕22領域5第6染色体短腕21領域前記免疫疾患または閉塞性肺疾患が、気管支喘息である、請求項1に記載の方法。前記一塩基多型が、配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、請求項1または2に記載の方法。下記1または2の配列を有する免疫疾患または閉塞性肺疾患検査用プローブ。1配列番号1~24から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列。2配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基と連鎖不平衡の関係にある塩基を含む10塩基以上の配列、又はその相補配列。下記1または2の領域を増幅することのできる免疫疾患または閉塞性肺疾患検査用プライマー。1配列番号1~24から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域。2配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基と連鎖不平衡の関係にある塩基を含む領域。本発明は、気管支喘息等の免疫疾患または閉塞性肺疾患の発症リスクや発症の有無を正確に検査する方法、及び該方法に用いられる検査試薬を提供することを課題とする。本発明者らは上記課題の解決のために鋭意検討した結果、第4染色体長腕31領域4q31領域、第5染色体長腕22領域5q22領域、第6染色体短腕21領域6p21領域、第10染色体短腕14領域10p14領域、または第12染色体長腕13領域12q13領域に存在する一塩基多型SNPが気管支喘息と関連することを同定した。そして、これらの多型を調べることにより気管支喘息等の免疫疾患または閉塞性肺疾患の発症リスクや発症の推定を正確に実施できることを見出し、本発明を完成するに至った。すなわち、本発明は以下の通りである。1以下の1~5のいずれかの領域に存在する一塩基多型を分析し、該分析結果に基づいて免疫疾患または閉塞性肺疾患を検査することを特徴とする、免疫疾患または閉塞性肺疾患の発症リスクおよび/または発症の有無の判定方法。1第4染色体長腕31領域2第12染色体長腕13領域3第10染色体短腕14領域4第5染色体長腕22領域5第6染色体短腕21領域2前記免疫疾患または閉塞性肺疾患が、気管支喘息である、1に記載の方法。3前記一塩基多型が、配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基、または該塩基と連鎖不平衡の関係にある塩基における一塩基多型である、1または2に記載の方法。4下記1または2の配列を有する免疫疾患または閉塞性肺疾患検査用プローブ。1配列番号1~24から選ばれる塩基配列において、塩基番号61番目の塩基を含む10塩基以上の配列、又はその相補配列。2配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基と連鎖不平衡の関係にある塩基を含む10塩基以上の配列、又はその相補配列。5下記1または2の領域を増幅することのできる免疫疾患または閉塞性肺疾患検査用プライマー。1配列番号1~24から選ばれる塩基配列において、塩基番号61番目の塩基を含む領域。2配列番号1~24から選ばれる塩基配列の塩基番号61番目の塩基に相当する塩基と連鎖不平衡の関係にある塩基を含む領域。本発明によれば、これまで予測が困難であった気管支喘息等の免疫疾患または閉塞性肺疾患の発症リスク罹患リスクを正確かつ簡便に予測することができる。また、気管支喘息等の免疫疾患または閉塞性肺疾患の発症を正確かつ簡便に判定することができる。したがって、本発明は気管支喘息等の免疫疾患または閉塞性肺疾患の予防や早期治療に貢献するものである。<1>本発明の方法本発明の方法は、ヒトの第4染色体長腕31領域4q31領域、第5染色体長腕22領域5q22領域、第6染色体短腕21領域6p21領域、第10染色体短腕14領域10p14領域、または第12染色体長腕13領域12q13領域に存在するSNPを分析し、該分析結果に基づいて免疫疾患または閉塞性肺疾患を検査することを特徴とする、免疫疾患または閉塞性肺疾患の発症リスクおよび/または発症の有無の判定方法である。すなわち、本発明において、検査とは免疫疾患または閉塞性肺疾患の発症リスクの検査及び免疫疾患または閉塞性肺疾患の発症の有無の検査を含む。本発明の方法においては、SNPの分析結果を、免疫疾患または閉塞性肺疾患の発症リスクおよび/または発症の有無と関連付ける。免疫疾患としては、特に限定されないが、例えば、気道アレルギー疾患やアトピー性皮膚炎などのアレルギー疾患が好ましく、気道アレルギー疾患がより好ましい。気道アレルギー疾患としては、例えば、気管支喘息が好ましい。閉塞性肺疾患としては、特に限定されないが、例えば、COPD慢性閉塞性肺疾患、Chronicobstructivepulmonarydisease、気管支喘息、慢性気管支炎、肺気腫、びまん性汎細気管支炎などが好ましく、気管支喘息がより好ましい。気管支喘息には成人気管支喘息、小児気管支喘息、アトピー性気管支喘息などが含まれるが、成人気管支喘息が好ましい。4q31領域に存在する具体的なSNPとしては、ヒトrs7686660およびrs3805236を挙げることができる。ここで、rs番号はNationalCenterforBiotechnologyInformationのdbSNPデータベースhttpwwwncbinlmnihgovprojectsSNPの登録番号を示す。これらSNPsは、4q31領域のUSP38遺伝子およびGAB1遺伝子が含まれる連鎖不平衡ブロックに位置する。よって、特に、当該連鎖不平衡ブロックに存在するSNPを解析することによって免疫疾患または閉塞性肺疾患を検査することができる。USP38遺伝子としては、具体的には、GenBankAccessionNoNC00000411の144106070~144143141の領域が挙げられる。GAB1遺伝子としては、具体的には、GenBankAccessionNoNC00000411の144257983~144395718の領域が挙げられる。rs7686660はGenBankAccessionNoNT01635418の68551306番目の塩基におけるチミンT/グアニンGの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3805236はGenBankAccessionNoNT01635418の68905884番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GA>AAの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。5q22領域に存在する具体的なSNPとしては、ヒトrs1837253を挙げることができる。このSNPは、5q22領域のTSLP遺伝子およびWDR36遺伝子が含まれる連鎖不平衡ブロックに位置する。よって、特に、当該連鎖不平衡ブロックに存在するSNPを解析することによって免疫疾患または閉塞性肺疾患を検査することができる。TSLP遺伝子としては、具体的には、GenBankAccessionNoNC0000059の110407390~110413722の領域が挙げられる。WDR36遺伝子としては、具体的には、GenBankAccessionNoNC0000059の110427870~110466200の領域が挙げられる。また、特に、TSLP遺伝子はこの領域における気管支喘息関連遺伝子である可能性が高い。rs1837253はGenBankAccessionNoNT0347725の12816885番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がCである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。6p21領域に存在する具体的なSNPとしては、ヒトrs204995、rs204993、rs404860、rs443198、s3130299、rs412657、rs570963、rs10484560、rs3129943、rs3117098、rs10947262、rs3806156、rs3129890、rs7775228、rs6457617、rs2858308、rs9275698、rs9500927を挙げることができる。これらSNPsは、6p21領域の主要組織適合遺伝子複合体領域MHCregionに位置する。よって、特に、MHC領域に存在するSNPを解析することによって免疫疾患または閉塞性肺疾患を検査することができる。rs204995は、GenBankAccessionNoNT00759214の23012535番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs204993は、GenBankAccessionNoNT00759214の23013831番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs404860は、GenBankAccessionNoNT00759214の23042595番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs443198は、GenBankAccessionNoNT00759214の23048656番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3130299は、GenBankAccessionNoNT00759214の23061787番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GA>AAの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs412657は、GenBankAccessionNoNT00759214の23069335番目の塩基におけるアデニンA/シトシンCの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs570963は、GenBankAccessionNoNT00759214の23147844番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs10484560は、GenBankAccessionNoNT00759214の23156387番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3129943は、GenBankAccessionNoNT00759214の23196945番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3117098は、GenBankAccessionNoNT00759214の23216763番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GA>AAの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs10947262は、GenBankAccessionNoNT00759214の23231562番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がCである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3806156は、GenBankAccessionNoNT00759214の23231948番目の塩基におけるチミンT/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs3129890は、GenBankAccessionNoNT00759214の23272523番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs7775228は、GenBankAccessionNoNT00759214の23516329番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs6457617は、GenBankAccessionNoNT00759214の23522101番目の塩基におけるアデニンA/グアニンGの多型を意味し、この塩基がAである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、AA>AG>GGの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs2858308は、GenBankAccessionNoNT00759214の23528250番目の塩基におけるチミンT/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs9275698は、GenBankAccessionNoNT00759214の23546223番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs9500927は、GenBankAccessionNoNT00759214の23819611番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がTである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、TT>TC>CCの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。10p14領域に存在する具体的なSNPとしては、ヒトrs10508372を挙げることができる。rs10508372はGenBankAccessionNoNT0775692の3334914番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がCである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。12q13領域に存在する具体的なSNPとしては、ヒトrs2069408およびrs1701704を挙げることができる。これらSNPsは、12q13領域のIKZF4遺伝子Eos遺伝子等が含まれる連鎖不平衡ブロックに位置する。よって、特に、当該連鎖不平衡ブロックに存在するSNPを解析することによって免疫疾患または閉塞性肺疾患を検査することができる。IKZF4遺伝子Eos遺伝子としては、具体的には、GenBankAccessionNoNC00001211の56414689~56432219の領域が挙げられる。rs2069408はGenBankAccessionNoNT02941911の18507627番目の塩基におけるチミンT/シトシンCの多型を意味し、この塩基がCである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、CC>CT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。rs1701704はGenBankAccessionNoNT02941911の18555793番目の塩基におけるチミンT/グアニンGの多型を意味し、この塩基がGである場合は免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。また、遺伝子型を考慮して解析した場合は、GG>GT>TTの順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。なお、上記24個のSNPについて、SNP塩基及びその前後60bpの領域を含む合計121bpの長さの配列を、配列番号1~24に示した。61番目の塩基が多型を有する。また、上記24個のSNPの詳細については、後述の表2にも示す。本発明においては、上記塩基に相当する塩基を解析する。上記塩基に相当する塩基とは、上記領域における該当塩基を意味する。すなわち、上記塩基に相当する塩基を解析することには、仮に人種の違いなどによって上記配列がSNP以外の位置で若干変化したとしても、上記領域における該当塩基を解析することが含まれる。また、本発明において解析する塩基は上記のものに限定されず、上記の塩基と連鎖不平衡にある塩基の多型を分析してもよい。ここで上記の塩基と連鎖不平衡にある塩基とは、上記の塩基とr2>0.5、好ましくはr2>0.8、さらに好ましくはr2>0.9の関係を満たす塩基をいう。r2は連鎖不平衡係数である。また上記の塩基と連鎖不平衡にある塩基は、例えば、HapMapデータベース等を用いて同定することができる。もしくは、複数人通常は20-40人程度から採取したDNAをシークエンサーにて配列解析し、連鎖不平衡にあるSNPを探索することにより同定することもできる。上記の塩基と連鎖不平衡にある塩基は、遺伝子型を考慮して解析した場合は、リスクアレルのホモ接合体>リスクアレルと非リスクアレルのへテロ接合体>非リスクアレルのホモ接合体の順で免疫疾患または閉塞性肺疾患の可能性または発症リスクが高い。上記SNPの塩基の種類を調べ、得られた結果を上記のような基準に基づいて免疫疾患または閉塞性肺疾患と関連付けることにより、免疫疾患または閉塞性肺疾患を検査することができる。上記SNPは単独で解析されてもよいし、上記SNPの少なくとも1つを含む複数のSNPsをまとめて解析ハプロタイプ解析してもよい。例えば、上記SNPの複数をまとめて解析してもよいし、上記SNPの少なくとも1つと、免疫疾患または閉塞性肺疾患と関連する既知のSNPs非特許文献1~9や当該既知のSNPsと連鎖不平衡にあるSNPsとを組み合わせて解析してもよい。具体的には、例えば、rs7686660、rs1837253、rs404860、rs10508372、およびrs1701704からなる群から選択される少なくとも2つ以上のSNPsを組み合わせて解析してもよい。免疫疾患または閉塞性肺疾患と関連する複数のSNPsをまとめて解析すれば、免疫疾患または閉塞性肺疾患の検査の精度が向上する。なお、いずれのSNPも、二本鎖DNAのどちらの鎖を解析してもよい。例えば、MHC領域の遺伝子やTSLP遺伝子の配列はセンス鎖を解析してもよいし、アンチセンス鎖を解析してもよい。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~24から選ばれる塩基配列の61番目の塩基を含む配列、又はその相補配列を有する10塩基以上の長さのプローブや、当該塩基と連鎖不平衡の関係にある塩基を含む配列、又はその相補配列を有する10塩基以上の長さのプローブが挙げられる。なお、当該塩基と連鎖不平衡の関係にある塩基及びその前後の領域の塩基配列は、例えば、NationalCenterforBiotechnologyInformationのdbSNPデータベースhttpwwwncbinlmnihgovprojectsSNPから取得できる。プローブの長さは好ましくは、15~35塩基であり、より好ましくは20~35塩基である。また、プライマーとしては、上記SNP部位を増幅するためのPCRに用いることのできるプライマー、又は上記SNP部位を配列解析シークエンシングするために用いることのできるプライマーが挙げられる。具体的には、配列番号1~24から選ばれる塩基配列の61番目の塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーや、当該塩基と連鎖不平衡の関係にある塩基を含む領域を増幅したりシークエンシングしたりすることのできるプライマーが挙げられる。このようなプライマーの長さは10~50塩基が好ましく、15~35塩基がより好ましく、20~35塩基がさらに好ましい。上記SNP部位をシークエンシングするためのプライマーとしては、上記塩基の5’側領域、好ましくは30~100塩基上流の配列を有するプライマーや、上記塩基の3’側領域、好ましくは30~100塩基下流の領域に相補的な配列を有するプライマーが例示される。PCRによる増幅の有無で多型を判定するために用いるプライマーとしては、上記塩基を含む配列を有し、上記塩基を3’側に含むプライマーや、上記塩基を含む配列の相補配列を有し、上記塩基の相補塩基を3’側に含むプライマーなどが例示される。なお、本発明の検査用試薬はこれらのプライマーやプローブに加えて、PCR用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。以下、本発明を実施例によりさらに具体的に説明する。但し、本発明はこれらの実施例に限定されない。1気管支喘息と関連するSNPsの同定気管支喘息感受性を決定する遺伝的変異を同定するために、日本人被検者を用いてゲノムワイド関連解析GWASを行った。GWASとは、疾患等の表現型に関わる遺伝的変異を探索する遺伝統計学的手法である。例えば、ヒトゲノム全体を網羅するような数十万~100万ヶ所のSNPsを用いて、ある疾患の患者ケースとその疾患にかかっていない被験者コントロールとの間で、多型の頻度に差があるかどうかを統計的に検定することで、疾患と関連する遺伝的変異を見出すことができる。<被検者>本実施例に用いた被検者の特徴を表1に示す。GWASに用いた気管支喘息の被検者ケースcase1532名は、複数の病院で募集した。再現試験replicationstudyに用いた気管支喘息の被検者5639名は、東京大学医科学研究所のBioBankJapanに登録された気管支喘息の患者から採用した。気管支喘息の被検者は、いずれも日本人であり、米国胸部疾患協会AmericanThoracicSociety;ATSの基準に基づき医師により気管支喘息と診断された。GWASに用いた対照被検者コントロールcontrol3304名は、BBJに登録された2403名の患者、および大阪御堂筋ロータリークラブで募集した901名の健康なボランティアからなる。再現試験に用いた対照被検者24608名は、BBJに登録された患者から採用した。なお、対照被検者として用いたBBJに登録された患者とは、気管支喘息、アレルギー性鼻炎allergicrhinitis、アトピー性皮膚炎atopicdermatitis、肺線維症lungfibrosis、慢性閉塞性肺疾患ChronicObstructivePulmonaryDisease以外の疾患の患者である。対照被検者は、いずれも日本人である。本研究は東京大学医科学研究所および理化学研究所横浜研究所の倫理委員会によって承認され、全ての参加者からインフォームドコンセントを得た。<GWAS>気管支喘息の被検者1532名の遺伝子型を、IlluminaHumanHap610QuadGenotypingBeadChipIllumina社USAを用いて解析した。対照被検者3304名の遺伝子型を、IlluminaHumanHap550GenotypingBeadChipIllumina社USAを用いて解析した。クオリティコントロールとして、主成分分析PCAにより外れ値outlierであると判定された1名の被検者を除外し、さらに、マイナーアレル頻度MAFsが0.01未満のSNPs、コール率callrateが0.99未満のSNPs、HardyWeinberg平衡検定でのP値がカットオフ値未満コントロールでP<106のSNPsを除外した。クオリティコントロールを通過した合計458847個のSNPsについて、コクランアーミテージ傾向検定CochranArmitagetrendtestにより気管支喘息との関連を評価した。なお、genomicinflationfactorGCは107であり、サンプルサイズで補正したλ1000は103であり、集団の偏りによる偽陽性の関連が得られる可能性は低いと考えられた。GWASの結果、6p21のMHCclassIII領域に存在するrs8192565に気管支喘息との最も強い関連が認められた。また、その他、NOTCH4遺伝子座とC6orf10遺伝子座の2つのSNPsが、ゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<5.0×108を満たし、気管支喘息との有意な関連を示した。さらに、ここで得られた日本人のGWASの結果を用いて、気管支喘息と関連することが既に報告されている領域について再現性の確認を行なった。例えば、最近行なわれたヨーロッパのコンソーシアムによる大規模メタ解析では、喘息とゲノムワイドレベルの強い関連をしめす領域が10カ所報告されている非特許文献7。日本人のGWASでは、10カ所の内、6p21のrs1063355、15q22のrs11071559、15q22のrs744910の3ヶ所で気管支喘息との関連が再現された。<再現試験>成人気管支喘息に関連するゲノム領域をさらに同定するため、GWASで解析した被検者とは独立の被検者を用いて、計102個のSNPsについて再現試験replicationstudyを行った。当該102個のSNPsは、GWASのコクランアーミテージ傾向検定でP<1×104を満たしたSNPsの内、当該102個のSNPsのいずれかとr2>0.9で連鎖不平衡の関係にある31個のSNPsを除いたものである。気管支喘息の被検者5639名の遺伝子型は、TaqManassayAppliedBiosystemsまたはmultiplexPCRinvaderassayThirdWaveTechnologiesにより解析した。対照被検者24608名の遺伝子型は、IlluminaHumanHap610QuadGenotypingBeadChipIllumina社USAを用いて解析した。再現試験の結果、102個のSNPsの内、5つの領域に存在する計24個のSNPsが、ボンフェローニ補正後の有意性の閾値であるP<4.9×104=0.05/102を満たし、気管支喘息との有意な関連を示した表2のreplication欄。これら24個のSNPsの詳細を表2に示す。さらに、これら24個のSNPsについて、MantelHaenszel法を用いてGWASと再現試験の統合解析combinedanalysisを行った結果、24個全てがゲノムワイドな有意性genomewidesignificanceの閾値として設定されたP<5.0×108を満たし、気管支喘息との有意な関連を示した表2のcombined欄。なお、統合解析の結果、各領域で気管支喘息との最も強い関連が認められたSNPsは以下のとおりである。4q315q226p2110p1412q13Risk;リスクアレルRiskAlleleRAF;リスクアレル頻度RiskAlleleFrequencya:GWASおよび再現試験replicationstudyのP値は、コクランアーミテージ傾向検定により算出した。b:GWASのP値は、λ値で補正したPGCを示す。c:OR;オッズ比OddsRatioおよびCI;信頼区間ConfidenceIntervalは、非リスクアレルを対照として算出した。d:統合解析combinedanalysisの結果は、MantelHaenszel法により算出した。これら5つの領域の内、5q22のTSLP遺伝子の転写開始点の57kb上流に存在するrs1837253は、カナダ人集団を被検者とする候補遺伝子解析により気管支喘息と気道過敏性との関連が報告されているが非特許文献8、ゲノムワイド水準でTSLP遺伝子座と気管支喘息とが関連するかは知られていなかった。また、6p21のMHC領域HLA遺伝子座に位置するMHCclassIIのHLADQ領域は、白人集団を被検者とする研究により気管支喘息および気道過敏性との関連が報告されている非特許文献7、9が、HLADQ領域以外のMHC領域と気管支喘息との関連は報告されていなかった。これら2つの領域はいずれも、アジア人集団における気管支喘息との関連は報告されておらず、これらの領域のSNPsが日本人等のアジア人集団における気管支喘息等の免疫疾患または閉塞性肺疾患の検査に有用であることは本実施例により初めて見出されたものである。また、4q31、10p14、12q13の3領域については、最近行なわれたヨーロッパのコンソーシアムによる大規模メタ解析非特許文献7でも気管支喘息との関連は報告されておらず、本実施例において初めて気管支喘息との関連が見出された。なお、当該ヨーロッパの大規模メタ解析に用いられた気管支喘息の被検者はその654が小児期発症喘息の被検者であったのに対し、本実施例では小児期発症喘息の被検者は196に過ぎず、成人喘息の被検者が中心であった。気管支喘息は多様性が高い疾患であり、年齢に特徴的な病態も認められることから、本実施例で解析した集団とヨーロッパの大規模メタ解析で用いられた集団とでは喘息の形質が異なる可能性がある。統合解析で気管支喘息との有意な関連を示した24個のSNPsの内、18ヶ所のSNPsは6p21の主要組織適合遺伝子複合体領域MHCregion、HLA遺伝子座ともいうに位置している。その内、最も強い関連はMHCclassIII領域のrs404860統合解析でP407×1023、OR121で認められた。GWASのデータに基づく6p21領域の関連解析の結果と連鎖不平衡マップを図1に示す。最近行なわれたヨーロッパのコンソーシアムによる大規模メタ解析では、MHCclassIIのHLADQ領域に位置するrs9273349と気管支喘息との関連が報告されており、当該SNPは、小児期発症childhoodonsetよりも成人発症lateronsetに強く影響するとされている非特許文献7。そこで、rs9273349と、本実施例のGWASで気管支喘息との最も強い関連を示したrs8192565統合解析ではP<5.0×108を満たさなかったため表2には記載せずとの関係を、GWASのデータを用いたロジスティック回帰分析により検討した。なお、本実施例のGWASにはrs9273349のデータが含まれていないため、代わりに、HapMapJPTおよびCEUデータベースでrs9273349とr2=1で完全連鎖するrs1063355のデータを用いた。その結果、rs1063355で補正しても、MHCclassIII領域のSNPsは気管支喘息と強く関連することが明らかとなった。すなわち、本実施例において見出されたMHCclassIII領域と気管支喘息との関連は、既に報告されていたHLADQ領域と気管支喘息との関連とは独立であると考えられる。また、6p21のAGER遺伝子座に存在するrs2070600は、1秒率FEV1FVCに関連する事が報告されている。1秒率とは呼吸機能の指標の1つであり、1秒率の低下は気管支喘息や慢性閉塞性肺疾患等の閉塞性肺疾患の特徴である。本実施例においてrs2070600は気管支喘息との関連を示さなかったが、統合解析で気管支喘息との最も強い関連を示したrs404860はrs2070600の近傍329kb下流に位置し、これらSNPsは弱い連鎖不平衡を示すr20115、D0805。このことから気管支喘息と呼吸機能低下とで共通の遺伝要因が存在する可能性も考えられる。さらに、MHC領域以外の4つの関連領域において、気管支喘息と関連する遺伝子を同定するため、GWASのデータを用いて関連の強さのマッピングを行なった。マッピングには、日本人のHapMapphaseIIIIIrelease27の情報をもとに、マイナーアレル頻度minorallelefrequency;MAF≧0.05でr2<0.8のSNPsをタグSNPstagSNPsとして選出して用いた。各領域の関連解析の結果と連鎖不平衡マップを図2~5に示す。13個のタグSNPsを用いた5q22領域の詳細なマッピングの結果、TSLPとWDR36の2つの遺伝子を含む約88kbの連鎖不平衡ブロックが気管支喘息の関連領域として同定され、最も強い関連はrs1837253で認められた図2。TSLP遺伝子は、気道ウイルス感染、複数のプロテアーゼ活性を持つアレルゲン、炎症性サイトカイン、およびタバコにより気道上皮細胞において強く誘導される遺伝子であり、環境要因の感知や、環境要因の暴露後のTh2細胞応答の誘導において重要な役割を果たしている。よって、TSLP遺伝子は、この領域における気管支喘息関連遺伝子である可能性が高い。22個のタグSNPsを用いた10p14領域の詳細なマッピングの結果、rs10508372が存在する約112kbの連鎖不平衡ブロックが気管支喘息の関連領域として同定された図3。この領域はいわゆる遺伝子砂漠と呼ばれる既知の遺伝子が存在しない領域であるが、1Mb上流にはTh2細胞の分化における主要な調節因子をコードするGATA3遺伝子が存在する。エンハンサーはゲノム上Mb単位で離れた遺伝子の発現を調節することも可能であり、この気管支喘息の関連領域上に存在する遺伝子多型がGATA3遺伝子発現をシスとして調節する活性があるのか等、詳細な検討が待たれる。12個のタグSNPsを用いた12q13領域の詳細なマッピングの結果、rs1701704が存在する約201kbの連鎖不平衡ブロックが気管支喘息の関連領域として同定された図4。rs1701704は、これまでにI型糖尿病や円形脱毛症との関連が報告されている。この連鎖不平衡ブロックには13個の遺伝子が含まれており、rs1701704はIKZF4遺伝子Eos遺伝子の2kb上流に位置していた。IKZF4遺伝子は、FoxP3依存性遺伝子サイレンシングFoxP3dependentgenesilencingの共調節因子として調節性T細胞の分化に関与することが知られている。調節性T細胞は肺のホメオスタシスを保ち、本来は無害である吸入抗原に対する有害な免疫応答を抑制する機能がある。11個のタグSNPsを用いた4q31領域の詳細なマッピングの結果、rs1397527が存在する約497kbの連鎖不平衡ブロックが気管支喘息の関連領域として同定された図5。この連鎖不平衡ブロックには、USP38とGAB1の2つの遺伝子が含まれている。USP38遺伝子は機能未知のubiquitinspecificpeptidase38をコードする。GAB1遺伝子はscaffoldingadopterproteinをコードし、IL3、IL6、IFNα、およびIFNγ等のサイトカイン受容体、並びにB細胞およびT細胞受容体を介して活性化されるシグナル伝達経路において重要な役割を果たす。さらに、本実施例で同定された5つの気管支喘息の関連領域について、日本人以外の人種においても気管支喘息との関連が認められるか検討した。499名のヒスパニック系でない白人集団からなる成人喘息の被検者、およびIlluminasiControldatabaseから取得した639名の対照被検者を用いて、5つの領域の5SNPs統合解析の結果、各領域で気管支喘息との最も強い関連が認められたrs7686660、rs1837253、rs404860、rs10508372、rs1701704について関連解析を行った。結果、5q22に存在するrs1837253で気管支喘息との有意な関連が認められた。このことから、ゲノム上のTSLPを含む領域に日本人と白人との共通の成人喘息の遺伝要因が存在する可能性が示唆される。