Integrated transcriptomic and proteomic analysis reveals molecular targets of vasopressin in GT1-7 GnRH neurons

技術要約：GT1-7 GnRHニューロンにおけるバソプレシンの統合的なトランスクリプトームおよびプロテオーム解析

問題提起
アルギニン・バソプレシン（AVP）は、視床下部ー下垂体ー性腺（HPG）軸の重要な構成要素であるゴナドトロピン放出ホルモン（GnRH）ニューロンの分泌活性を調節することが知られている。視交叉上核（SCN）および前視索におけるキスペプチン–GnRHネットワークを介したAVPシグナル伝達は十分に解明されているが、AVPがGnRHニューロン細胞に影響を与える直接的な分子メカニズムについては、依然として十分に特徴付けられていない。具体的には、AVP刺激に対するGnRHニューロンモデルの包括的なトランスクリプトームおよびプロテオームの応答に関する包括的なデータが不足している。さらに、広く用いられている不死化GnRHニューロンモデルであるGT1-7細胞株において、標準的なAVP受容体（AVPR1A、AVPR1B、AVPR2）が強固に定義されておらず、観察されたAVP誘発性応答の受容体基盤に関する不確実性が存在する。

方法論
本研究では、GT1-7細胞株を用いた統合的なマルチオミクス・アプローチを用いて、AVP誘発性の分子変化をマッピングした。

• 実験デザイン： GT1-7細胞を(Arg8)-バソプレシン 200 pg/mlで6時間処理し、ビークル処理コントロールと比較した。各群に対して3つの独立した生物学的レプリケートを用いた。
• トランスクリプトミクス： 全RNAを抽出し、Illumina HiSeq X Tenプラットフォームを用いたRNA-seqを実施した。発現変動解析にはDESeq2を用いた（閾値：|logFC| > 1, P < 0.05）。
• プロテオミクス： タンパク質を抽出し、消化後、iTRAQ 8-plexキットを用いてラベル化した。定量的解析はLC–MS/MSにより行った。データ処理にはProteome DiscovererおよびMascotを用い、発現変動はフォールド変化 > 1.2 または < 0.8 かつ P < 0.05 と定義した。
• バイオインフォマティクス： 機能濃縮解析（GO/KEKE）はDAVIDおよびGSEAを用いて実施した。タンパク質–タンパク質相互作用（PPI）ネットワークはSTRINGを用いて構築し、Cytoscapeで可視化した。ハブ遺伝子／タンパク質は、5つの中心性アルゴリズム（MCC、MNC、degree、closeness、stress）を用いて特定した。
• 統合： mRNAとタンパク質の両レベルで有意に変化した分子を特定するために、クロスレベル比較を行った。
• 検証： 標準的なAVP受容体の検出が低かったことを踏まえ、潜在的なシグナル伝達メカニズムを探るため、本研究ではGPR54/KISS1RアンタゴニストであるKiss234を利用した。これは、AVPに応答する細胞生存率、細胞内Ca²⁺、ATPレベル、GnRH分泌、およびプロモーター活性へのGPR54関連シグナルの寄与を評価するために使用された。また、ラットの弓状核を含む視床部切片を用いて、近接結紮アッセイ（PLA）を行い、AVP–GPR54の近接性を調べた。

主な結果

• トランスクリプトームの変化： AVP処理により1,094個の変動発現遺伝子（DEG）（アップレギュレート549個、ダウンレギュレート545個）が誘導された。主要なアップレギュレート遺伝子にはNupr1、Slc7a11、Trib3が含まれ、Rgs4、Tmem119、Ccng1はダウンレギュレートされた。
  • パスウェイ： 濃縮解析により、TNF、PI3K–Akt、MAPK、およびcAMPシグナル伝達経路が強調された。GSEAはさらに、軸索ガイダンスおよびWntシグナル伝達経路を特定した。
  • ハブ遺伝子： ネットワーク解析により、Cxcl10、Cxcl1、Cxcl12、Itpkb、Icam1、Ptgs2、Mmp9が中心的なハブ候補として特定された。
• プロテオームの変化： 本研究では2,648個の変動発現タンパク質（DEP）（増加1,362個、減少1,286個）を特定した。
  • パスウェイ： DEPは代謝経路、スプライソソーム、および小胞体におけるタンパク質プロセシングに濃縮されていた。
  • ハブタンパク質： Polr2iおよびPolr2g（RNAポリメラーゼIIのサブユニット）が、すべての中心性アルゴリズムにおいて一貫してハブタンパク質として特定された。
• 統合解析： mRNAとタンパク質の両レベルで有意に変化した分子は、わずか73個であった。このうち、35個が一致してアップレギュレートされ、21個が一致してダウンレギュレートされていた。
  • 機能的焦点： 重複する分子は、主にアラニン、アスパラギン酸、グルタミン酸、グリシン、セリン、トレオニン、システイン、およびメチオニン代謝を含むアミノ酸代謝経路に濃縮されていた。
• 受容体およびシグナル伝達の検証： 標準的なAVP受容体（Avpr1a、Avpr1b、Avpr2）はGT1-7細胞において強固には検出されなかった。しかし、GPR54アンタゴニストであるKiss234の前処理は、GnRH分泌、細胞内Ca²⁺、ATPレベル、およびGnrh1プロモーター活性におけるAVP誘発性の変化を減弱させた。PLAは、弓状核を含むラット視床部切片において、AVP–GPR54の近接性を確認した。

意義と主張
著者らは、AVP刺激下におけるGnRHニューロンの初の包括的なマルチオミクスプロファイルを紹介している。本研究は、将来の機能的および生理学的調査における候補分子ターゲットの貴重なリソースを提供すると主張している。

主な主張は以下の通りである：

1. メカニズムの洞察： AVPは、協調的な転写、転写後、および代謝調節を通じてGT1-7細胞の生物学に影響を与える。トランスクリプトームとプロテオームのデータの部分的な重複（73分子）は、プロテオームの結果を形成する上で、転写後および翻訳後のメカニズムが主要な役割を果たしていることを示唆している。
2. ハブ制御因子： 特定のハブ遺伝子（例：ケモカインであるCxcl10、Cxcl12）およびハブタンパク質（RNAポリメラーゼIIサブユニット）の特定は、AVPがいかにしてニューロンの興奮性、炎症、および全般的な転写能を調節するかを理解するための具体的な候補を提供する。
3. シグナル伝達の文脈： 本研究はGT1-7細胞における標準的なAVP受容体シグナル伝達を確定したものではないが、GPR54関連のシグナル伝達が、このモデルにおけるAVP応答性の変化に寄与している可能性を示す支持的な証拠を提供している。著者らは、これらの知見は、古典的なAVP受容体の特異性に対する証明ではなく、GPR54関連経路の可能性を示す証拠として解釈されるべきであることを明示しており、正確な受容体基盤には一次ニューロンと標的アプローチを用いたさらなる研究が必要であると述べている。
4. 代謝調節： 本研究は、アミノ酸代謝が転写および翻訳の両レベルにおいてAVPによって調節される主要なプロセスであることを強調しており、代謝状態とGnRHニューロン機能との間の潜在的な関連性を示唆している。

本論文は、AVPがGT1-7細胞において多層的な調節応答を引き起こすと結論付けており、これらの分子応答の全容を捉えるためには、統合的なトランスクリプトームおよびプロテオームのアプローチが不可欠であるとしている。