JP5327773B2

JP5327773B2 - メタボリック症候群の検査方法

Info

Publication number: JP5327773B2
Application number: JP2007196100A
Authority: JP
Inventors: 紀久子堀田
Original assignee: RIKEN Institute of Physical and Chemical Research
Current assignee: RIKEN Institute of Physical and Chemical Research
Priority date: 2007-07-27
Filing date: 2007-07-27
Publication date: 2013-10-30
Anticipated expiration: 2027-07-27
Also published as: JP2009027986A

Description

本発明は遺伝子解析によりメタボリック症候群を検査する方法に関する。

メタボリック症候群は肥満、糖代謝異常、脂質異常、高血圧などの複数の代謝異常を呈する病態で（非特許文献１）、近年、食生活やライフスタイルの変化に伴って増加している。メタボリック症候群は進行すると糖尿病や動脈硬化などの疾患の発症につながる危険性があるため、早めにメタボリック症候群発症を予測し、食生活の改善等に努めることが重要である。
メタボリック症候群と遺伝的要因の関連に関しても幾つか報告があり（非特許文献２、３）、非特許文献２ではメタボリック症候群と関連する染色体上の部位について報告されている。しかしながら、この部位とメタボリック症候群との関連は必ずしも明確ではなく、より正確にメタボリック症候群発症を予測できる遺伝子解析法の開発が求められていた。

ミオシンIXA（MYO9A）遺伝子は当初、Bardet-Biedl 症候群 (BBS)の原因遺伝子BBS4として同定されたが（非特許文献４）、その後、BBS4は別の遺伝子であったことが報告されたため（非特許文献５）、MYO9A遺伝子の機能は不明である。また、MYO9A遺伝子のメタボリック症候群との関わりについては全く知られていなかった。
Diabetes 53:2087-2094 (2004) Nat Clin Pract Cardiovasc Med 3:482-489 (2006) Obesity Research 13:381-490 (2005) Genomics 59:150-160 (1999) Nat Genet 28:188-191 (2001

本発明は、メタボリック症候群の発症や進行を正確に予測するための検査方法を提供することを課題とする。

本発明者らは、上記課題を解決するために鋭意検討を重ねた。その結果、MYO9A遺伝子上に存在する塩基の一塩基多型または該塩基と連鎖不平衡にある塩基の一塩基多型がメタボリック症候群に関連することを見出し、本発明を完成するに至った。

すなわち、本発明は以下の通りである。
（１）ミオシンIXA（MYO9A）遺伝子上に存在する塩基の一塩基多型または該塩基と連鎖不平衡にある塩基の一塩基多型を分析し、該分析結果に基づいてメタボリック症候群を検査する方法。
（２）前記一塩基多型が、配列番号１〜９のいずれかの塩基配列の塩基番号６１番目の塩基に相当する塩基における多型である、（１）の方法。
（３）配列番号１〜９のいずれかの塩基配列において、塩基番号６１番目の塩基を含む１０塩基以上の配列、又はその相補配列を有するメタボリック症候群検査用プローブ。
（４）配列番号１〜９のいずれかの塩基配列において、塩基番号６１番目の塩基を含む領域を増幅することのできるメタボリック症候群検査用プライマー。
（５）（３）のプローブまたは（４）のプライマーを含む、メタボリック症候群検査用キット。

本発明の検査方法により、メタボリック症候群の発症リスクや進行を正確に予測することができるため、メタボリック症候群を予防したり、糖尿病や動脈硬化への発展を抑えたりすることができる。

＜１＞本発明の検査方法
本発明の検査方法は、ミオシンIXA（MYO9A）遺伝子上に存在する塩基の一塩基多型（SNP；single nucleotide polymorphism）または該塩基と連鎖不平衡にある塩基の一塩基多型を分析し、該分析結果に基づいてメタボリック症候群を検査する方法である。メタボリック症候群としては、例えば、J Atheroscler Thromb 13:202-208, 2006に記載の基準に基づいて定められる病態が例示される。なお、本発明において、「検査」とは現在メタボリック症候群を発症しているどうかの検査、および将来メタボリック症候群になるかどうかを予測するための検査を含む。

MYO9A遺伝子としては、ヒトMYO9A遺伝子が好ましく、例えば、GenBank Accession No. NM_006901.1に登録された配列を有する遺伝子を挙げることができる。ただし、該遺伝子は人種の違いなどによって１又は複数の塩基に置換や欠失等が存在する可能性があるため、上記配列の遺伝子に限定されない。
MYO9A遺伝子上に存在する塩基の一塩基多型としては、MYO9A遺伝子のエクソンやイントロンに存在する塩基の一塩基多型だけでなく、MYO9A遺伝子の５’−隣接領域や３’−隣接領域に存在する塩基の一塩基多型も含まれる。

メタボリック症候群に関連するMYO9A遺伝子の一塩基多型は特に制限されないが、例えば、以下に示すSNP-1〜SNP-9が挙げられる。

表１において、例えば、SNP-5はGenBank Accession No. NT_010194.16の42961005番目の塩基におけるアデニン（A）／グアニン（G）の多型を意味し、この塩基がAである場合はメタボリック症候群になる確率が高い。また、アレルを考慮して解析した場合は、SNP-5がAA＞AG＞GGの順でメタボリック症候群になる確率が高い。SNP-5はMYO9A遺伝子のイントロン内に存在するSNPである。dbSNPはNational Center for Biotechnology InformationのdbSNPデータベース（//www.ncbi.nlm.nih.gov/projects/SNP/）の登録番号を示す。
その他のＳＮＰについても同様であるが、いずれについても左側の塩基がメタボリック症候群になりやすい塩基である（例えば、SNP-4の場合、A/GのうちAの場合がメタボリック症候群になりやすく、SNP-8の場合、T/GのうちTの場合がメタボリック症候群になりやすい）。
なお、SNP-1〜SNP-９について、ＳＮＰ塩基及びその前後６０ｂｐの領域を含む合計１２１ｂｐの長さの配列を、それぞれ配列番号１〜９に示した。それぞれ６１番目の塩基が多型を有する。
これらの塩基に相当する塩基を本発明においては解析する。ここで、「相当する」とは、ヒトMYO9A遺伝子上の上記配列を有する領域中の該当塩基を意味し、仮に、人種の違いなどによって上記配列がＳＮＰ以外の位置で若干変化したとしても、その中の該当塩基を解析することも含む。

上記ＳＮＰの塩基の種類を調べることによって、メタボリック症候群を検査することができる。検査するＳＮＰの数は、一種類でもよいし、複数（ハプロタイプ解析）でもよい。なお、SNP-1がGであれば他のSNPｓも左側の塩基（SNP-2＝C、SNP-3＝G、SNP-4＝A、SNP-5＝A、SNP-6＝C、SNP-7＝C、SNP-8＝T、SNP-9＝G）であるというように、上記SNPのタイプは互いに相関するため、一箇所のみを解析することによっても十分正確なメタボリック症候群の予測を行うことができる。なお、MYO9A遺伝子の配列はセンス鎖を解析してもよいし、アンチセンス鎖を解析してもよい。

また、本発明において解析する塩基は上記のものに限定されず、上記の塩基と連鎖不平衡にある塩基の多型を分析してもよい。ここで「上記の塩基と連鎖不平衡にある塩基」とは、上記の塩基とｒ²＞０．５の関係を満たす塩基をいう。

MYO9A遺伝子の一塩基多型の解析に用いる試料としては、染色体ＤＮＡを含む試料であれば特に制限されないが、例えば、血液、尿等の体液サンプル、細胞、毛髪等の体毛、爪などが挙げられる。一塩基多型の解析にはこれらの試料を直接使用することもできるが、これらの試料から染色体ＤＮＡを常法により単離し、これを用いて解析することが好ましい。

MYO9A遺伝子の一塩基多型の解析は、通常の一塩基多型解析方法によって行うことができる。例えば、シークエンス解析、ＰＣＲ、ハイブリダイゼーションなどが挙げられるが、これらに限定されない。

シークエンスは通常の方法により行うことができる。具体的には、多型を示す塩基の５’側数十塩基の位置に設定したプライマーを使用してシークエンス反応を行い、その解析結果から、該当する位置がどの種類の塩基であるかを決定することができる。なお、シークエンスを行う場合、あらかじめ多型を含む断片をＰＣＲなどによって増幅しておくことが好ましい。

また、ＰＣＲによる増幅の有無を調べることによっても解析することができる。例えば、多型を示す塩基を含む領域に対応する配列を有し、かつ、各多型に対応するプライマーをそれぞれ用意する。それぞれのプライマーを使用してＰＣＲを行い、増幅産物の有無によってどのタイプの多型であるかを決定することができる。

また、多型を含むDNA断片を増幅し、増幅産物の電気泳動における移動度の違いによってどのタイプの多型であるかを決定することもできる。このような方法としては、例えば、PCR-SSCP（ｓｉｎｇｌｅ−ｓｔｒａｎｄｃｏｎｆｏｒｍａｔｉｏｎｐｏｌｙｍｏｒｐｈｉｓｍ）法（Ｇｅｎｏｍｉｃｓ．１９９２Ｊａｎ１；１２（１）：１３９−１４６．）が挙げられる。具体的には、まず、MYO9A遺伝子の多型部位を含むＤＮＡを
増幅し、増幅したＤＮＡを一本鎖ＤＮＡに解離させる。次いで、解離させた一本鎖ＤＮＡを非変性ゲル上で分離し、分離した一本鎖ＤＮＡのゲル上での移動度の違いによってどのタイプの多型であるかを決定することができる。

さらに、多型を示す塩基が制限酵素認識配列に含まれる場合は、制限酵素による切断の有無によって解析することもできる（ＲＦＬＰ法）。この場合、まず、ＤＮＡ試料をＰＣＲで増幅し、それを制限酵素により切断する。次いで、ＤＮＡ断片を分離し、検出されたＤＮＡ断片の大きさによってどのタイプの多型であるかを決定することができる。

ハイブリダイゼーションの有無を調べることによって多型の種類を解析することも可能である。すなわち、各多型に対応するプローブを用意し、いずれのプローブにハイブリダイズするかを調べることによってSNPがいずれの塩基であるかを調べることもできる。

＜２＞本発明の検査用試薬
本発明はまた、メタボリック症候群を検査するためのプライマーやプローブなどの検査試薬を提供する。このようなプローブとしては、MYO9A遺伝子における上記多型部位を含み、ハイブリダイズの有無によって多型部位の塩基の種類を判定できるプローブが挙げられる。具体的には、配列番号１〜９のいずれかの塩基配列の６１番目の塩基を含む配列、又はその相補配列を有する１０塩基以上の長さのプローブが挙げられる。プローブの長さはより好ましくは１５〜３５塩基である。

また、プライマーとしては、MYO9A遺伝子における上記多型部位を増幅するためのPCRに用いることのできるプライマー、又は上記多型部位を配列解析（シークエンシング）するために用いることのできるプライマーが挙げられる。具体的には、配列番号１〜９のいずれかの塩基配列の６１番目の塩基含む領域を増幅したりシークエンシングしたりすることのできるプライマーが挙げられる。このようなプライマーの長さは１０〜５０塩基が好ましく、１５〜３５塩基がより好ましい。
上記多型部位をシークエンシングするためのプライマーとしては、上記塩基の５’側領域、好ましくは30〜100塩基上流の配列を有するプライマーや、上記塩基の３’側領域、好ましくは30〜100塩基下流の領域に相補的な配列を有するプライマーが例示される。PCRによる増幅の有無で多型を判定するために用いるプライマーとしては、上記塩基を含む配列を有し、上記塩基を３’側に含むプライマーや、上記塩基を含む配列の相補配列を有し、上記塩基の相補塩基を３’側に含むプライマーなどが例示される。
なお、本発明の検査用試薬はこれらのプライマーやプローブに加えて、PCR用のポリメラーゼやバッファー、ハイブリダイゼーション用試薬などを含むものであってもよい。

以下、本発明を実施例によりさらに具体的に説明する。但し、本発明はこれらの実施例に限定されない。

日本各地の医療機関から集められたメタボリック症候群患者および健常者ボランティアを被検者とし、以下のようにして、メタボリック症候群に関連するＳＮＰｓの探索を行った。

下記の基準（J Atheroscler Thromb 13:202-208, 2006）でメタボリック症候群と診断された日本人７２９名（男性３３７名、女性３９２名：平均年齢５４±１４歳）を症例とした。
＜メタボリック症候群＞
１．トリグリセリド値１５０ｍｇ／ｄＬ以上及び／又はＨＤＬ−コレステロール値４０ｍｇ／ｄＬ未満、あるいはこれらの脂質異常を治療中
２．収縮期血圧１３０ｍｍＨｇ以上及び／又は拡張期血圧８５ｍｍＨｇ以上、あるいは高血圧の治療中
３．空腹時血糖１１０ｍｇ／ｄＬ以上、または糖尿病の治療中
上記１〜３のうちの２項目以上に該当し、ＢＭＩ２５ｋｇ／ｍ²以上の者をメタボリック症候群と診断した。（なお、一般のメタボリック症候群の診断にはＢＭＩの代わりにウエスト周囲が使用される。すなわち、ウエスト周囲が男性で８５ｃｍ以上、女性で９０ｃｍ以上、かつ、上記１〜３のうちの２項目以上に該当すればメタボリック症候群と診断される。ウエスト周囲径、男性８５ｃｍ、女性９０ｃｍは臍の高さで測定したCTによる内臓脂肪面積が１００ｃｍ²に相当する。）

＜コントロール１＞
上記基準によりメタボリック症候群に該当しなかった日本人１５２４名（男性５４３名、女性９８１名：平均年齢５３±１６歳）をコントロール１とした。
＜コントロール２＞
コントロール２はコントロール１の中から、ＢＭＩ２５ｋｇ／ｍ²未満で上記代謝異常のない者を選んだ（４４１名（男性１１２名、女性３２９名：平均年齢４７±１６歳）。
症例、コントロール１、コントール２の性質を表２にまとめた。
なお、いずれの被検者からも書面のインフォームドコンセントを得たうえで、各病院施設の倫理委員会及び理化学研究所の承認を得たプロトコルに基づいて解析を行った。

＜ＤＮＡの調製とＳＮＰジェノタイピング＞
試料として、各被検者の血液サンプルから染色体ＤＮＡを、通常の手順にて調製した。
Obesity Research 13:381-490., 2005に挙げられた遺伝子群におけるSNPｓからメタボリック症候群関連SNPｓを探索した。まずは、その遺伝子群におけるSNPｓをＪＳＮＰデータベース（IMS-JST: Institute of Medical Science-Japan Science and Technology Agency Japanese SNP database：Nucleic. Acids. Res., 30, 158-162）より探索したところ、１２２遺伝子について合計７５５個のSNPｓが見つかった。そのうち、マイナーアレル頻度が０．２より大きく、予想されるアレル頻度がＨＷＥ（ｐ＞０．００１）から大きく逸れない３３６個のSNPｓを本試験では採用した。なお、HWEテストはHardy-Weinberg equilibrium test（Nealsen, D.M., Ehm, M.G. & Weir, B.S. Am. J. Hum. Gent. 63, 153-1
540 (1998)）を示す。
これらのSNPsに対してインベーダープローブ（Third Wave Technologies）を設計し、症例群とコントロール群について、マルチＰＣＲとインベーダーアッセイ（Ohnishi, Y. et al. J. Hum. Genet. 46, 471-477 (2001)）の組み合わせにより各SNPｓをジェノタイピングした。

まず、症例（Case）群とコントロール１群を用いて、ケース−コントロール相関スタディを行った。該スタディにおいて最もP値が小さかった１０個のSNPｓを表３に示す。また、これら１０個のSNPｓに関して症例群とコントロール２群を用いてケース−コントロール相関スタディを行った結果を表４に示す。

なお、表４におけるアレル１，２は表３と同じ塩基を示している。表３、４において、NPYはNeuro peptide Y遺伝子、MKKSはMcKusick-Kaufman syndrome遺伝子を示す。

表３，４において、症例群とコントロール群の間でのジェノタイプとアレルの頻度はFisher's exact testによって解析した。相関テストは３種類のモードで行った。（１）アレル頻度モード：２×２の分割表を用いてケースとコントロール間でアレル頻度を比較した、（２）優性モード：２×２の分割表を用いてアレル１のホモザイゴートとヘテロザイゴートを他（アレル２のホモザイゴート）と比較した、（３）劣性モード：２×２の分割表を用いてアレル１のホモザイゴートの頻度をその他と比較した。危険率（Odds ratio）及び95%信頼区間（c.i.）はWoolfの方法に従った。さらに、各SNP解析において並べ替え
検定を行い、各SNPｓについて最小P値の経験分布を作成し、その分布に基づいて得られた経験P値を補正P値とした。

次に、症例群とコントロール１群でリニアディスクリミナント分析（Annals of Eugenics 7:179-188., 1936）を行った。ここでは、年齢（歳）、性別（男性／女性、0/1）、BMI（kg/m²）、空腹時血糖（mg/dL）、トリグリセリド（mg/dL）、HDL-コレステロール（mg/dL）、収縮期血圧（mmHg）、拡張期血圧（mmHg）、糖尿病・脂質異常・高血圧の治療（治療無／有、0/1）を変数に用いた。
結果を図１に示す。
症例群とコントロール群はある程度区別できたが、中間スコアの者が多かった。そして、ディスクリミナントスコアの相関比は０．５３であり、症例群とコントロール群を完全には区別できていないと考えられた。

そこで、上記10個のSNPｓとｒｓ２２７８９８７（MYO9A遺伝子）のディスクリミナントスコアに対する影響を症例群とコントロール１群で調べた。結果を表５に示す。その結果、MYO9A遺伝子のｒｓ２３０６５７６の劣性ホモ（AA）のディスクリミナントスコアがAG+GGに対して有意に小さかった。（P=0.000077）また、MYO9A遺伝子のｒｓ２３０６４８９、ｒｓ３７８４３１３およびｒｓ２２７８９８７の劣性ホモも低いP値を示した。

Rs番号はＪＳＮＰデータベースの番号である。

以上によって、MYO9A遺伝子のSNPｓがメタボリック症候群に相関することがわかった。
MYO9A遺伝子には上記４個のSNPｓ以外にもSNPｓが存在する。そこで、マイナーアレル頻度が５％以上のSNPｓを抽出し、各SNPｓに関してLD（連鎖不平衡）係数（Genomics 29:311-322, 1995）を算出したところ、表６のようになり、これらのSNPｓが連鎖不平衡にあることがわかった。したがって、上記４個のSNPｓ以外にもｒｓ２６２５５２９、ｒｓ２２９１２８０、ｒｓ２３０６５７５、ｒｓ２３０６４８８、ｒｓ３８０９５０５のSNPｓもメタボリック症候群に相関することが示唆された。

メタボリック症候群は加齢とともに増加し、男性のほうが女性より６−７倍多いことが報告されている。そこで、性別と年齢の影響を調べるため、logistic回帰分析を行った。Objective変数にはメタボリック症候群か否かの２値変数を使用し、Explanatory変数は年齢、性別、SNP（ｒｓ２３０６５７６）を使用した。統計解析にはsoftware R (http://www.r-project.org/)を使用した。
結果を表７に示す。それによると、年齢及び性別で補正した後もMYO9A遺伝子のSNPｓはメタボリック症候群と有意に相関することがわかった。すなわち、MYO9A遺伝子のSNPｓは性別、年齢などの他の変数の影響を受けることなく、メタボリック症候群に相関すること
が示唆された。

CTスキャンによって測定された内臓脂肪面積が１００ｃｍ²以上であり、上記１〜３のうちの２項目以上を満たす者をケースとし、内臓脂肪面積が１００ｃｍ²未満で上記項目のいずれにも該当しない者をコントロールとして、MYO9A遺伝子のSNP（ｒｓ２３０６５７６）との相関を調べた。その結果、表８に示されるように、MYO9A遺伝子のSNPは、内臓脂肪面積が１００ｃｍ²以上で上記１〜３のうちの２項目以上を満たすという指標によって診断されるメタボリック症候群に対しても有意に相関することがわかった。

なお、表８におけるアレル１，２は表３と同じ塩基を示している。

ディスクリミナントスコアのプロット。横軸がサンプルを示し、縦軸がディスクリミナントスコアを示す。黒丸がコントロール１を示し、×印が症例を示す。

Claims

ミオシンIXA（MYO9A）遺伝子上に存在する、配列番号１，４，５，８のいずれかの塩基配列の塩基番号６１番目の塩基の一塩基多型を分析し、該分析結果に基づいて、配列番号１の６１番目の塩基に相当する塩基がＧの場合、メタボリック症候群になる確率が高い；配列番号４の６１番目の塩基に相当する塩基がＡの場合、メタボリック症候群になる確率が高い；配列番号５の６１番目の塩基に相当する塩基がＡの場合、メタボリック症候群になる確率が高い；配列番号８の６１番目の塩基に相当する塩基がＴの場合、メタボリック症候群になる確率が高いと判定される、日本人におけるメタボリック症候群を検査する方法。
配列番号１，４，５，８のいずれかの塩基配列において、塩基番号６１番目の塩基を含む１５塩基以上の配列、又はその相補配列からなるプローブを含む、日本人におけるメタボリック症候群検査用試薬。
配列番号１，４，５，８のいずれかの塩基配列において、塩基番号６１番目の塩基を含む領域を増幅することのできるプライマーであって、該塩基の３０〜１００塩基上流の１５塩基以上の配列、又は該塩基の３０〜１００塩基下流の領域に相補的な１５塩基以上の配列からなるプライマーを含む、日本人におけるメタボリック症候群検査用試薬。
請求項２に記載のプローブまたは請求項３に記載のプライマーを含む、日本人におけるメタボリック症候群検査用キット。