存在することを明らかにしました。発症機構や新たな診断法・治療法の開発につながると期待できます。炎症性腸疾患は大腸や小腸などに慢性の炎症や潰瘍を起こす難病の一つで、日本においても患者数は年々増加しており、新薬開発が望まれています。これまでゲノムワイド関連解析(GWAS)[2]によって、200以上の疾患発症に関わる因子である遺伝子バリアント[3]が明らかになりました。グループは、これまでのゲノム解析に加え、独自に構築したeQTL[4](遺伝子発現量の個人差に関与するゲノム領域)データを新たな手法で統合解析しました。その結果、約100の新たな原因遺伝子候補を同定するとともに、予想に反して近傍の複数遺伝子の発現量が同時に変化することで発症につながっている可能性が示されました。今回開発した統合解析手法は他の疾患にも応用できると考えられます。Communications』(6月21日付け)に掲載されます。大腸や小腸などに慢性の炎症や潰瘍を起こす疾患の総称で、クローン病[5]と潰瘍性大腸炎[6]に分類されます。10代~20代の若者が発症しやすい傾向にありますが完治させることは難しく、継続的な経過観察および治療が必要となる難病の一つです。欧米諸国において患者数が多いことが知られていますが、日本でも患者数は年々増加しています。発症には衛生環境や食生活などに加えて、遺伝的な要因も影響を強く及ぼしていることが知られています。ジェネティクス・コンソーシアム[7]」が中心となり、これまでゲノムワイド関連解析(GWAS)を用いて、疾患発症に関わるゲノム領域が200以上明らかになりました注1)。それらのゲノム領域に存在すると想定される疾患発症の原因遺伝子同定は、疾患発症機構の解明や新しい治療薬の開発につながります。これまで、低頻度の遺伝子バリアントの解析が可能なイムノチップ解析[8]注2)、低頻度の遺伝子バリアントの集積を解明するシークエンス解析注3)、高度な統計的解析手法注4)などを用いて解析が行われましたが、原因遺伝子を同定できたのは一部のゲノム領域のみでした。国際共同研究グループはまず、GWASにより明らかになった遺伝子バリアントの多くは遺伝子の翻訳領域外に存在することから、「それらの遺伝子バリアントが疾患発症に重要な細胞内で、近傍に存在する遺伝子の発現量を変化させることで発症に至る」という仮説を立てました。この仮説に基づいた解析を行うためには、ゲノム全体にわたり、遺伝子バリアントと遺伝子発現量の関係を明らかにしたeQTL(遺伝子発現量の個人差に関与するゲノム領域)データが必要ですが、そのようなeQTLデータは細胞の種類によって異なるため、炎症性腸疾患の発症に重要だと考えられる免疫細胞や腸管の上皮細胞におけるeQTLデータが不可欠だと考えられました。そこで、リエージュ大学付属病院において、健康な約350人の血液から6種類の免疫細胞と、腸管の3カ所の上皮細胞をそれぞれ収集しました。続いて、各遺伝子の発現量を測定するとともに、ゲノム全体の遺伝子バリアントを解析し、合計23,650個の遺伝子バリアントと遺伝子発現量の関係を示したeQTLデータを構築しました。次に、疾患発症の原因遺伝子の同定を目指し、遺伝子バリアントと疾患の関係を示すGWASデータとeQTLデータを統合解析するための新たな手法を開発しました。図1のように、遺伝子バリアントが近傍遺伝子の発現量の変化を引き起こすことが発症につながる場合は、最も疾患への影響が大きい遺伝子バリアントは遺伝子の発現量も最も大きく変え、疾患への影響が弱いバリアントになるにつれ発現量の変化も小さくなると予想されることから、GWASデータとeQTLデータの関連(P値:観察した事象が偶然に起こる確率)のパターンが類似していると考えられます。図1左には、GWASデータとeQTLデータに共通して、同じゲノム位置に遺伝子バリアントの関連の強さを示すピークが見られます。図1右のように、それぞれの遺伝子バリアントの関係を発現量の増減も考慮した散布図で表わすと、強い負の相関があることが分かります。このことから、遺伝子バリアントは対象遺伝子の発現量を低下させることで、疾患発症に関与すると考えられます。逆に、図2のように、遺伝子バリアントが近傍遺伝子の発現量変化に関わっているが、疾患発症には関与していないと考えられる場合は、異なるゲノム位置に遺伝子バリアントの関連の強さを示すピークが存在し、かつ互いに相関が認められません。ジェネティクス・コンソーシアムから入手したGWASデータと、200以上の疾患発症に関わるゲノム領域について統合解析を行った結果、約100個の遺伝子を原因遺伝子の有力な候補として同定しました。この中には、炎症性腸疾患の原因遺伝子としてすでに知られるATG16L1、CARD9、FUT2などや、免疫関連遺伝子のIL18R1、IL6ST、THEMISに加えて、APEH、ANKRD55、CISD1、CPEB4、DOCK7、ERAP2、GNA12、GPX1、GSDMB、ORMDL3、SKAP2、UBE2L3、ZMIZ1など、現状では疾患発症との関連が不明な遺伝子も含まれていました。今回、示した各ゲノム領域についてのGWASデータとeQTLデータの関係については、本研究のウェブサイト注5)において確認することができます。各ゲノム領域からGWASデータとeQTLデータの両方に類似のパターンを示す遺伝子は一つだけと予想しましたが、図3に示すように、一つのゲノム領域から複数の遺伝子のeQTLとGWASデータが同じパターンを示しているものが25カ所も同定されました。これは、これまでの想定と異なり、複数の近傍遺伝子が同時に疾患発症に寄与するという複雑な制御機構を持つ可能性を示しています。最後に、得られた原因遺伝子候補の中に真の原因遺伝子が実際に含まれているかを検証するために、低頻度の遺伝子バリアントの集積という概念を用いました。目的の遺伝子について疾患群と対照群のシークエンス解析を行って低頻度の遺伝子バリアントを同定し、それが疾患群または対照群の一方に偏りがあれば、それは原因遺伝子としてのさらなる証拠になるという考え方です。チームで開発したターゲットシークエンス手法注6)を用いて解析しました。その結果、疾患群と対照群の間に強い偏り(P<6.9×10-4)が認められ、今回同定した遺伝子群の中に多く真の原因遺伝子が含まれていることが示されました。今後、同定された100個の各遺伝子について炎症性腸疾患にどのようにつながるか機能や動物モデルの解析を行うことで、炎症性腸疾患の発症機構や、新たな診断法・治療法の開発につながると期待できます。今回開発したeQTLデータとGWASデータの統合解析手法は、他の疾患にも応用できると考えられます。E-MTAB-6666)および遺伝子の発現量のデータ(E-MTAB-6667)は、欧州バイオインフォマティクス研究所(EBI)で公開しています。炎症や潰瘍を起こす難病。クローン病と潰瘍性大腸炎に分類される。10~20歳代の若者が発症しやすく、下痢や血便、体重減少が見られる。完治が難しく生涯にわたって経過観察・治療が必要なため、新薬開発が望まれている。ゲノムワイド関連解析(GWAS)疾患の関連遺伝子を見つける代表的な方法。数が多いほど検出力が上がるため、近年はコンソーシアムを結成し、大規模検体を用いた解析が行われている。遺伝子バリアントヒトのDNA配列は30億の塩基対からなるが、その配列の個人間の違いを遺伝子バリアントと呼ぶ。eQTL遺伝子の発現量の個人差に関与するゲノム領域。