国際バイオバンク連携によるヒト疾患リスク遺伝子アトラスを構築

日本主導型の国際メタアナリシスでゲノム情報に基づく疾患の精密分類を提案

日本主導型の国際メタアナリシスでゲノム情報に基づく疾患の精密分類を提案

2021-10-01大阪大学,東京都健康長寿医療センター,東北大学東北メディカル・メガバンク機構,東京大学,理化学研究所,日本医療研究開発機構

研究成果のポイント

バイオバンク・ジャパン18万人のゲノムデータを基に過去最大220の健康・医療データ(多因子疾患・希少疾患・バイオマーカー^※1・服薬データ)に対する網羅的ゲノムワイド関連解析(GWAS)^※2を実施。
国際バイオバンク連携によりイギリス・フィンランドのバイオバンク^※3と計63万人のメタアナリシスを実施し、5,000以上の新規遺伝的リスク関連領域を発見。GWAS結果を公開するデータベースPheWeb.jpを構築し、データシェアリングを推進。
GWAS要約統計量^※4の特異値分解^※5と、東北大学東北メディカル・メガバンク機構のメタボローム^※6データ等、他のオミクスデータへのプロジェクション解析^※7を実施することで、ゲノム情報を基にした従来の疾患分類の精密・層別化方法を提唱。

概要

全世界のヒトゲノムデータの蓄積に伴ってこの20年ほどで実施されてきたGWASにより、ヒトの遺伝的変異と将来の疾患発症リスクとの関係が網羅的に判明しつつあります。得られたゲノムデータを基に、長い医学の歴史から構築されてきた多彩なヒトの病気の分類について、客観的に振り返ってみたらどうなるでしょうか?この疑問に適切に答えるためには、従来のゲノム研究の課題である、①:ゲノム研究データが欧米人集団に偏重して構築されていること、②:研究対象となる疾患の網羅性が低かったこと、③:得られた大量のGWAS結果を医学的・生物学的に解釈する方法論が確立されていないこと、の3点を解決する必要がありました。

大阪大学大学院医学系研究科の坂上沙央里助教(研究当時、現ハーバード大学医学部博士研究員)、金井仁弘特別研究学生(ハーバード大学医学部　博士課程)、岡田随象教授(遺伝統計学/理化学研究所生命医科学研究センター自己免疫疾患研究チーム　客員主管研究員)、東京大学大学院新領域創成科学研究科　松田浩一教授らの研究グループは、国際バイオバンク連携を通じて、バイオバンク・ジャパン(日本)・UKバイオバンク(英国)・FinnGen(フィンランド)の計63万人のゲノムデータと健康・医療データの網羅的な解析を実施しました。多因子疾患・希少疾患・バイオマーカー・服薬データまでを網羅する過去最大220のヒト形質のゲノムワイド関連解析(GWAS)により、ヒト疾患に関わる5000以上の新規遺伝的リスク関連領域が発見されました。研究グループは、各疾患の遺伝的リスク構造が遺伝的集団^※8を超え共有されていることを示すとともに、得られたゲノム解析アトラスを医療に役立てる方法として、GWAS要約統計量を特異値分解し遺伝学を基にした疾患群の再分類を試みました(図1)。東北大学東北メディカル・メガバンク機構が構築したメタボロームデータ等、他のオミクスデータへのプロジェクション解析、エピゲノム情報やパスウェイ等との多分野融合的な解析を実施し、GWAS要約統計量行列を数理学的に分解する工程を経て、より複雑な疾患(例:心筋梗塞)のリスクを単純な形質の足し合わせとして解釈すること(例:コレステロールと血圧)や、ゲノム変異を基にした疾患の精密分類・層別化方法の提唱(例:I型アレルギーとIV型アレルギー疾患の分類)に成功しました。さらに、本研究結果を広く利用促進することを目的にデータシェアリングサイトPheWeb.jp を構築しデータを無償・非制限公開し、Polygenic risk score(PRS)^※9をはじめとするゲノム個別化医療の社会実装基盤として広く利用されることが期待されます。

図1日本主導の国際バイオバンク連携により、220のヒト形質のGWASと横断的メタアナリシスを実施して結果を公開するとともに、医療に役立てる方法論を提唱。

本研究成果は、米国科学誌「Nature Genetics」に、2021年10月1日(金)午前0時(日本時間)に公開されました。

研究の背景

古来より、ヒト疾患の診断分類は、患者の多彩な症候を観察し、似た症状を分類整理し、それに基づき最適な治療を選ぶために経験的に構築されてきました。先端的な検査方法や画像診断が発達した現在においても、従来の診断分類は徐々に形を変えながら踏襲されていますが、ゲノムデータを基にこの分類を客観的に振り返ってみたらどうなるでしょうか?ゲノムに基づく疾患分類によって従来の分類の正しさを確かめることができるかもしれませんし、あるいは似通った疾患群によりよい精密な分類を与えることが可能になるかもしれません。数十万人規模のゲノムデータと多彩な臨床データを蓄えたバイオバンクが各国に存在する現在、国際共同研究を通じてこの疑問に初めて答えを出すことができるようになってきました。

本研究の成果

研究グループは、国際バイオバンク連携を通じてバイオバンク・ジャパン(日本)・UKバイオバンク(英国)・FinnGen(フィンランド)の計63万人のゲノムデータと、多因子疾患・希少疾患・バイオマーカー・服薬データまでを網羅するヘルスデータをテキストマイニングを用いて網羅的に解析し、過去最大220のヒト形質のゲノムワイド関連解析(GWAS)を実施しました。本研究では、従来のゲノム研究の課題であった欧米人集団へのゲノムデータ偏重に対して、アジア人集団として最大規模を誇るバイオバンク・ジャパンのデータを基底に、UKバイオバンク・FinnGenと協調した幅広い形質を対象としたことに意義があります。その結果、集団特異的変異を含む14,814の疾患感受性領域の発見に至り、このうち5,343の領域が新規報告となりました。今回初めて網羅的に多数の疾患の遺伝的構造を日本人集団と欧米人集団とで比較したところ、各疾患の遺伝的リスク構造は遺伝的集団を超えて共通しており、ゲノムによって定義された類似疾患グループ(例:心筋梗塞・不安定狭心症などの心血管疾患)も遺伝的集団間で共有されていることがわかりました。これは、地域や環境要因を超えて歴史的に形成されてきた診断分類の妥当性を、大規模ゲノム解析によって再確認できた知見と言えるかもしれません。

現在のゲノム研究の最大の課題は、GWAS解析から得られた大量の要約統計量が複雑すぎて、それらを生物学的に解釈し、直接的に医療に役立てるのが困難であることでした。そこで研究グループは、全ゲノムに渡る遺伝的変異ごとの疾患リスク統計量行列を、特異値分解(SVD)という数理統計学手法を用いて行列分解しました。さらに分解の結果得られたゲノムに基づく疾患リスク要素を、バイオマーカーGWASや東北大学東北メディカル・メガバンク機構で実施したメタボロームGWAS結果へ投影(プロジェクション)しました(図2)。これにより、心筋梗塞のリスク変異がコレステロールと血圧を説明するリスク変異の組み合わせで説明されるなど、複雑な疾患のリスクの背景をより単純な形質のリスク要素の和で説明したり、それらの要素を説明する生物学的パスウェイと重要な細胞組織を特定したりすることが可能になりました。この分解要素を基にして、似通った複数のアレルギー疾患を、喘息やアレルギー性鼻炎など古典的にI型アレルギーと分類されてきた疾患群と、金属アレルギーや接触性皮膚炎など古典的にIV型アレルギーと分類されてきた疾患群との2群に分類することに成功しました。この結果は、遺伝学を基にした疾患診断の再分類の成功例として意義深いものと考えています。

図2要約統計量行列に特異値分解を行い、多彩な疾患群のオミクス情報を基に再整理し、異なる疾患に共通する要素を抽出して、その背景にある生物学的パスウェイを特定することに成功。

本研究成果が社会に与える影響(本研究成果の意義)

本研究の成果はわが国発のゲノム解析データ基盤として今後のさらなる国際連携のために重要であり、GWAS解析結果をデータシェアリングサイトPheWeb.jp にて無償・非制限公開を実施しています。このデータを基にした新規ゲノム解析手法の開発や、Polygenic risk score(PRS)をはじめとするゲノム個別化医療の社会実装基盤として広く利用されることが期待されます。

特記事項

本研究成果は、2021年10月1日(金)午前0時(日本時間)〔9月30日(木)午後4時(英国時間)〕に米国科学誌「Nature Genetics」(オンライン)に掲載されました。

タイトル

“A cross-population atlas of genetic associations for 220 human phenotypes”

著者名

Saori Sakaue^{1-5, 46,}*, Masahiro Kanai^{1, 5-9, 46}, Yosuke Tanigawa¹⁰, Juha Karjalainen^5-7,9, Mitja Kurki^5-7,9, Seizo Koshiba^11,12, Akira Narita¹¹, Takahiro Konuma¹, Kenichi Yamamoto^1,13,14, Masato Akiyama^2,15, Kazuyoshi Ishigaki^2-5, Akari Suzuki¹⁶, Ken Suzuki¹, Wataru Obara¹⁷, Ken Yamaji¹⁸, Kazuhisa Takahashi¹⁹, Satoshi Asai^20,21, Yasuo Takahashi²¹, Takao　Suzuki²², Nobuaki Shinozaki²², Hiroki Yamaguchi²³, Shiro Minami²⁴, Shigeo Murayama²⁵, Kozo Yoshimori²⁶, Satoshi Nagayama²⁷, Daisuke Obata²⁸, Masahiko Higashiyama²⁹, Akihide Masumoto³⁰, Yukihiro Koretsune³¹, FinnGen, Kaoru Ito³², Chikashi Terao², Toshimasa Yamauchi³³, Issei Komuro³⁴, Takashi Kadowaki^33,35, Gen Tamiya^11,12,36,37, Masayuki Yamamoto^11,12,36, Yusuke Nakamura^38,39, Michiaki Kubo⁴⁰, Yoshinori Murakami⁴¹, Kazuhiko Yamamoto¹⁶, Yoichiro Kamatani^2,42, Aarno Palotie^5,9,43, Manuel A. Rivas¹⁰, Mark J. Daly^5-7,9, Koichi Matsuda^44,*, Yukinori Okada^{1,2,14,42,45,}*
(*責任著者)

所属

大阪大学大学院医学系研究科　遺伝統計学
理化学研究所生命医科学研究センター　ゲノム解析応用研究チーム
ハーバード大学医学部　Center for Data Sciences
ブリガム&ウィメンズ病院　Divisions of Genetics and Rheumatology
ブロード研究所　Program in Medical and Population Genetics
マサチューセッツ総合病院　Analytic and Translational Genetics Unit
ブロード研究所　Stanley Center for Psychiatric Research
ハーバード大学医学部　Department of Biomedical Informatics
ヘルシンキ大学　Institute for Molecular Medicine Finland (FIMM)
スタンフォード大学　Department of Biomedical Data Science
東北大学　東北メディカル・メガバンク機構
東北大学　未来型医療創成センター(INGEM)
大阪大学大学院医学系研究科　小児科学
大阪大学免疫学フロンティア研究センター　免疫統計学
九州大学医学部　眼科
理化学研究所生命医科学研究センター　自己免疫疾患研究チーム
岩手医科大学　泌尿器科学講座
順天堂大学大学院　膠原病・リウマチ内科学
順天堂大学大学院　呼吸器内科学
日本大学医学部　生体機能医学系薬理学分野
日本大学医学部　臨床試験研究センター
徳洲会グループ
日本医科大学　血液内科
日本医科大学　生体機能制御学分野
東京都健康長寿医療センター
結核予防会　複十字病院
がん研究会　がん研究所
滋賀医科大学　医学部附属病院臨床研究開発センター
大阪国際がんセンター　呼吸器外科
麻生飯塚病院
国立病院機構大阪医療センター
理化学研究所生命医科学研究センター　循環器ゲノミクス・インフォマティクス研究チーム
東京大学医学部附属病院　糖尿病・代謝内科
東京大学医学部附属病院　循環器内科
虎の門病院
東北大学大学院医学系研究科
理化学研究所　革新知能統合研究センター
東京大学医科学研究所　ヒトゲノム解析センター
がんプレシジョン医療研究センター
理化学研究所生命医科学研究センター
東京大学医科学研究所　癌・細胞増殖部門　人癌病因遺伝子分野
東京大学大学院新領域創成科学研究科　メディカル情報生命専攻　複雑形質ゲノム解析分野
マサチューセッツ総合病院　Psychiatric & Neurodevelopmental Genetics Unit
東京大学大学院新領域創成科学研究科　メディカル情報生命専攻　クリニカルシークエンス分野
大阪大学先導的学際研究機構　生命医科学融合フロンティア研究部門
共同筆頭著者

DOI番号

10.1038/s41588-021-00931-x

なお、本研究は、日本医療研究開発機構(AMED)ゲノム医療実現バイオバンク利活用プログラム:B-cureのうち、ゲノム研究バイオバンク(旧:疾患克服に向けたゲノム医療実現プロジェクト(オーダーメイド医療の実現プログラム))、ならびにゲノム医療実現推進プラットフォーム・先端ゲノム研究開発:GRIFIN「遺伝統計学に基づく日本人集団のゲノム個別化医療の実装」の一環として行われ、大阪大学大学院医学系研究科バイオインフォマティクス・イニシアティブの協力を得て行われました。本研究で使用したサンプルは、「オーダーメイド医療の実現プログラム」において収集されたものです。また、本研究で使用した一部のデータはAMEDによる「東北メディカル・メガバンク計画」において得られました。

用語説明

※1　バイオマーカー: ヒトの疾患病態や生物的特徴の説明に役立つ、測定指標のこと。体重、血圧などの身体測定値や、コレステロールや尿酸などの血液検査値が含まれる。
※2　ゲノムワイド関連解析研究(GWAS:Genome Wide Association Study): ヒトゲノム配列上に存在する数百～数千万カ所の遺伝的変異(とくに一塩基多型・SNP)とヒトの疾患発症リスクや個人差(形質)との関連を網羅的に検討する、遺伝統計解析手法。
※3　バイオバンク: 疾患疫学や病態研究などを目的に、多数のヒトのDNA、血清、尿、組織などの検体を収集、蓄積、管理する施設のこと。近年では国家レベルで数十万人を対象とするバイオバンクが構築され、個人の検体とともに電子カルテ上の臨床情報やその後の予後などの追跡情報も蓄積される例が多い。
※4　要約統計量(summary statistics): GWASの関連解析結果の主となる統計量をまとめた表のことで、一般的には全ゲノムに渡る遺伝的変異の名称や位置情報と、線形回帰やロジスティック回帰の結果のP値(=統計学的な有意性)・beta値(=リスクの方向性・大きさ)・標準誤差、解析サンプル数などの情報がまとめられたもののことを指す。元のゲノムデータや形質データと異なりプライバシーの問題がないため、GWAS結果を他機関に共有する際の基本的な情報として利用される。
※5　特異値分解(SVD:singular value decomposition): 線形代数において、正方行列に限らない任意の行列を複数の行列の積として分解する手法の一つ。
任意の実行列Aに対して得られる特異値分解A=UΣVでは、ΣはAの特異値を対角項にもつ矩形行列、UとVは直交行列となる。本論文ではtruncated SVDという応用手法により、指定した階数(< rank(A))までの分解を実施し、複雑な遺伝子変異―疾患間の遺伝的相関関係を、得られた成分の組み合わせとして表現した。
※6　メタボローム(metabolome): 生体内の低分子化合物質(代謝物:metabolite)の総体(-ome)のこと。質量分析法等により網羅的に代謝物を測定する技術が近年発達し、数千種類におよぶ代謝物の量を測定することができる。
※7　プロジェクション解析: 疾患形質に対する特異値分解によって得られた遺伝的変異と各要素の関係の大きさと符号を表す重み情報を用いて、新たなGWAS要約統計量を特異値分解の要素に投影し、それぞれの要素に対してどのくらいの貢献度があるか定量する解析のこと。
※8　遺伝的集団(population): 英語におけるpopulation、すなわちある遺伝学的背景を共有した任意の集団(データから科学的に定義された任意の集団)のことで、本研究においては日本語訳の一例として「遺伝的集団」を使用している。より狭義のancestry(遺伝学的背景を共有した集団かつ、歴史的にも遺伝的祖先を共有していると考えられる集団、日本語訳の例:遺伝的祖先)や、社会的に構築されたニュアンスを含むrace(日本語訳の例:人種)やethnicity(日本語訳の例:民族)とは異なる意味で使用している。
※9　Polygenic risk score(PRS): 大規模なゲノムワイド関連解析研究(GWAS)により疾患や形質との関連が示唆された数十～数千の遺伝的変異の重み付きの和を個人ごとに計算したスコア。このスコアは実際の疾患発症リスクと相関することが示されており、集団内でスコアの分布を調べることで、特にその疾患のリスクが高い個人を特定することができる。