背景

骨格筋は人体において、グルコース代謝とインスリン反応を主導する最大の代謝器官である。第2型糖尿病や脂質異常症などの心血管代謝疾患のリスクを高める遺伝子変異の多くも、骨格筋の代謝経路と深く関わっている。ゲノムワイド関連解析（GWAS）は、これまで数千に及ぶ疾患リスク遺伝子座を発見してきたが、非コード領域に位置する非翻訳変異が疾患へとつながる分子経路は、依然として謎に包まれていた。

従来のトランスクリプトーム解析やゲノム研究は、組織破砕液を丸ごと解析するバルク（bulk）解析に依存してきた。筋線維をはじめ、筋衛星細胞、脂肪前駆細胞、内皮細胞、免疫細胞などが複雑に混ざり合う筋肉組織において、各細胞型や特定の生理学的状態による微細な遺伝子調節異常を感知することは、不可能に近いことだった。筋線維自体が多核細胞という構造的特性を持つため、従来の単一細胞分離技術を完全に適用することが難しいという技術的困難も障壁となっていた。疾患関連変異がどのような細胞核の状態において異常な調節を引き起こすのかを明らかにする、高解像度の分子マップが切実に求められていたのである。

主な発見

研究チームは、ヒト骨格筋組織から単離した単一核に基づき、クロマチン接近性を測定する単一核クロマチン接近性解析（snATAC-seq）と単一核トランスクリプトーム解析（snRNA-seq）を同時に実施し、統合マルチオミクス参照マップを構築した。筋線維の型別サブタイプ（Type I, Type IIa, Type IIx）および微小環境を構成する基質細胞核を精密に分離した後、クロマチンの開裂状態とトランスクリプトーム活性を単一核レベルで対照した。

単一核レベルの発現定量形質遺伝子座（eQTL）およびクロマチン接近性定量形質遺伝子座（caQTL）を体系的にマッピングした結果、特定の細胞状態や生理学的文脈においてのみ機能する調節変異群が明確に識別された。組織全体を解析した場合には平均的な信号に埋もれて見えなかった数千の細胞特異的QTLが、表面化したのである。代謝疾患のGWAS統計信号とこれらのQTL信号の共局在（colocalization）解析を組み合わせることで、疾患感受性を左右する因果的な遺伝子調節回路が明らかになった。インスリンシグナル伝達およびグルコース輸送体（GLUT4）の移動経路に関与する遠位調節因子部位において、特定の疾患リスクアレルがクロマチン接近性を選択的に遮断し、遺伝子発現を歪曲させるメカニズムが直接確認された。

意義と展望

本研究は、心血管代謝疾患を誘発する非コード遺伝子変異が、骨格筋内部の特定の細胞型や機能的状態においてのみ病理的な転写障害を引き起こす事実を明確に立証した。遺伝子レベルに留まっていた標的探索の地平を、細胞状態特異的な調節部位へと絞り込むことで、新薬候補物質の選定における失敗率を大幅に低減させる基盤が整った。

ただし、静的な組織生検試料を中心にマップが作成されたため、運動や急激な食事変化などの動的な代謝刺激下において、クロマチン環境がどのように再構成されるかを追跡する後続研究が必要である。骨格筋生検検体の人種的多様性を拡大し、単一核マルチオミクスパイプラインの解析コストを低減させる標準化作業も、今後の課題である。