Proteome-wide QTL mapping enables gene-protein-phenotype metabolic network construction in a genetically diverse MASLD mouse model

技術要約：遺伝的多様性を持つMASLDマウスモデルにおけるプロテオームワイドQTLマッピング

問題提起
代謝異常関連脂肪性肝疾患（MASLD）およびその重症形態である代謝異常関連脂肪肝炎（MASH）は、世界中で約38%の成人に影響を及ぼしている。高い遺伝率（30〜60%）が推定されているにもかかわらず、疾患の病態生理を駆動する複雑な分子ネットワークの理解が不完全であるため、治療の選択肢は限られている。既存の動物モデルは、単一の遺伝的背景や食事誘発性の表現型に依存していることが多く、肝細胞のバルーニング（風船様変化）や線維化といった疾患進行の典型的な症状を十分に捉えられていない。さらに、トランスクリプトミクスは広く利用されてきたが、タンパク質の存在量は、細胞プロセスにおける機能的なエフェクターであり、mRNAレベルからは予測が困難な場合が多い。特に、代謝フラックスの中心となる翻訳後調節を受けるタンパク質においては顕著である。したがって、遺伝的変異、ディープ・プロテオミクス・プロファイリング、および表現型データを統合し、遺伝的多様性のある文脈において「遺伝子ータンパク質ー表現型」のネットワークをマッピングするリソースが切実に求められている。

方法論
著者らは、MASLDへの感受性をモデル化するために、遺伝的多様性を持つマウスコホートを構築した。これには、疾患感受性が異なる4つの創始系統（C57BL/6J、129S1/SvImJ、CAST/EiJ、およびPWK/PhJ）を用いた4元交雑が含まれ、444匹のF2マウスが生成された。コホートはサーモニョートラリティ（熱中立条件）で飼育され、肝損傷を誘発するために高脂肪のウェスタンダイエットが与えられた。

プロテオミクスのスループットにおける歴史的な制限を克服するため、本研究では次世代質量分析（Orbitrap Astral）を用いたデータ非依存的取得（DIA）を活用した。これにより、444個すべての肝臓サンプルから、1ヶ月足らずでディープ・プロテオミクス・データを取得することが可能となり、1サンプルあたり平均9,790個のタンパク質が検出された。研究では以下のカスタムパイプラインを用いた：

1. ジェノタイプ・マップの生成： 創始系統およびF2マウスのRNA-seqデータからバリアントをコールし、4つの創始者ハプロタイプを区別する高解像度のジェノタイプ・マップを構築した。
2. pQTLマッピングの実行： タンパク質存在量を制御するゲノム領域（pQTL）を特定するため、雄および雌の集団に対して個別に定量的形質遺伝子座（QTL）解析を実施し、cis-作用（遺伝子の±1 Mbp以内）とtrans-作用（遠位）のメカニズムを区別した。
3. 表現型の統合： 性別を共変量として考慮し、ケンドールのタウを用いてタンパク質存在量と組織学的脂肪変性スコアを相関させた。
4. 媒介分析（Mediation Analysis）： 候補となる媒介因子を特定するため、特定のタンパク質存在量（例：Ubr1）が、ゲノム領域の脂肪変性スコアまたは遠位のタンパク質存在量への影響を媒介しているかどうかを判断する媒介分析を行った。
5. ユビキチン・プロファイリング： ジグリシン（K-GG）残基濃縮を用い、候補タンパク質の直接的なユビキチン化標的を同定した。
6. ヒトでの検証： 同様のプロテオミクス手法を用いて、減量手術を受けた患者のヒト肝臓サンプル38検体を解析し、種を越えた一致性を評価した。

主な結果

• プロテオミクスの深度と多様性： 本研究では、コホート全体で12,411個のユニークなタンパク質が検出された。このうち、9,776個が50%以上のサンプルで検出され、QTL解析に使用された。
• pQTLの発見： 解析により2,652個のタンパク質に対するpQTLが同定され、4,000を超える異なるゲノム領域にマッピングされた。研究では、cis-pQTL（雌で2,418個、雄で2,269個）とtrans-pQTL（雌で3,924個、雄で4,560個）を区別した。既知のMASLD関連遺伝子であるPnpla3、Hsd17b13、Mboat7が確認されたほか、新規の領域も同定された。
• ゲノミック・ホットスポット： 第2染色体にtrans-pQTLの有意なクラスターが観察され、これは脂肪変性スコアの高いLODピークと重複していた。この領域には、微小管関連タンパク質およびMCM複合体（Mcm2-5、Mcm7）の制御が含まれており、これらは肝細胞癌（HCC）でアップレギュレートされるが、マウスにおいては線維化やHCCがまだ進行していない段階で脂肪変性と相関していることが判明した。
• 中心的な制御因子としてのUbr1の同定： 媒介分析により、主要な候補としてE3ユビキチンリガーゼであるUbr1が特定された。Ubr1の存在量は脂肪変性スコアと強い正の相関を示し（雄において領域の効果の約42%を媒介）、細胞内輸送、細胞骨格組織、および脂質代謝に関与する多数の遠位タンパク質の存在量を媒介していることが示された。
• メカニズムの洞察： ユビキチン部位プロファイリングにより、Ubr1の存在量がプロテオスタシス（例：リボソームタンパク質、シャペロニン）および脂質代謝（例：Slc27a2、Scd1、Hsd17b13）に関与するタンパク質の協調的なユビキチン化の変化と関連していることが明らかになった。特筆すべき点として、Ubr1はPlin2の分解を標的とすることで知られているが、本データは、代替的な要因によって駆動される可能性のある、増加したUbr1がPlin2の増加と相関するという複雑な調節関係を示唆していた。
• 種を越えた検証： ヒト肝臓コホートにおいて、UBR1の存在量はMASH患者で有意に上昇しており、肝脂肪量と相関していた。ヒトのUBR1レベルは、マウスの知見と同様に、ユビキチン化経路および代謝酵素の協調的な変化と相関していた。

意義および主張
本論文は、遺伝的変異、タンパク質存在量、および代謝表現型を結合する「マイニング可能なプロテオゲノミクス・リソース」を提供すると主張している。主な意義は、遺伝的多様性を持つ集団におけるQTLマッピングにおいて、ディープ・プロテオミクス・プロファイリングを統合することで、トランスクリプトミクス研究が見落とす可能性のある複雑な調節ネットワークを明らかにできることを示した点にある。

具体的には、著者らは以下の点を主張している：

1. 新規の調節ノード： プロテオスタシスと脂質代謝を繋ぐ中心的な制御因子としてUbr1を同定したことは、Ubr1の存在量に影響を与える遺伝的多型が、代謝機能不全への素因となる可能性を示唆している。
2. 手法の進歩： 本研究は、大規模なpQTLマッピングにおいて現代的なMS機器（Orbitrap Astral）およびDIAワークフローの有用性を検証し、数百のサンプルにおいてほぼ完全なプロテオーム・カバレッジを達成した。これにより、プロテオミクスが複雑な疾患構造を解明するためのゲノミクスやトランスクリプトミクスに匹敵する地位を確立した。
3. トランスレーショナルな関連性： 遺伝的多様性を持つマウスモデルとヒト肝臓コホートとの間の種を越えた一致性は、単一の背景を持つマウスモデルの歴史的な限界に対処し、ヒトの疾患への汎用性を持つ候補遺伝子を特定するための遺伝的参照集団の使用を支持している。
4. コミュニティ・リソース： 著者らは、pQTLマップの可視化、共調節パターンの探索、および種を越えた一致性の評価を、計算科学的な専門知識なしで行えるインタラクティブなウェブベースのツールであるMAPLE（MASLD Associated Proteomic Loci Explorer）を通じてデータを公開した。

本研究は、MASLDの感受性と進行を支配する遺伝子ータンパク質ネットワークを解明することにより、バイオマーカー探索および治療標的の設定に新たな洞察を提供すると結論付けている。