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

기술 요약: GT1-7 GnRH 뉴런 내 바소프레신에 대한 통합적 전사체 및 단백질체 분석

문제 제기
아르기닌 바소프레신(AVP)은 성선 자극 호르몬 방출 호르몬(GnRH) 뉴런의 분비 활동을 조절하는 것으로 알려져 있으며, 이는 시상하부-뇌하수체-성선(HPG) 축의 핵심 구성 요소이다. 시삭핵(SCN) 및 전시삭 영역의 키스펩틴-GnRH 네트워크를 통한 AVP 신호 전달은 잘 입증되어 있으나, AVP가 GnRH 뉴로날 세포에 영향을 미치는 직접적인 분자 기전은 여전히 명확히 규명되지 않았다. 구체적으로, AVP 자극에 대한 GnRH 뉴런 모델의 전반적인 전사체 및 단백질체 반응에 관한 포괄적인 데이터가 부족하다. 또한, 널리 사용되는 불멸화된 GnRH 뉴런 모델인 GT1-7 세포주에서 정통 AVP 수용체(AVPR1A, AVPR1B, AVPR2)가 강력하게 정의되지 않았으며, 이는 관찰된 AVP 유도 반응의 수용체 기반에 대한 불확실성을 야기한다.

방법론
본 연구는 AVP 유도 분자 변화를 매핑하기 위해 GT1-7 세포주를 이용한 통합적 멀티오믹스 접근법을 채택하였다.

• 실험 설계: GT1-7 세포를 (Arg8)-바소프레신 200 pg/ml 농도로 6시간 동안 처리하였으며, 이를 용매 처리된 대조군과 비교하였다. 각 그룹당 3개의 독립적인 생물학적 반복 실험을 수행하였다.
• 전사체 분석(Transcriptomics): 총 RNA를 추출하여 Illumina HiSeq X Ten 플랫폼을 통해 RNA-seq를 수행하였다. 차등 발현 분석은 DESeq2를 사용하여 수행하였다 (임계값: |logFC| > 1, P < 0.05).
• 단백질체 분석(Proteomics): 단백질을 추출, 소화시킨 후 iTRAQ 8-plex 키트로 라벨링하였다. 정량 분석은 LC–MS/MS를 통해 수행되었다. 데이터는 Proteome Discoverer와 Mascot으로 처리되었으며, 차등 발현은 fold change > 1.2 또는 < 0.8 및 P < 0.05로 정의하였다.
• 생물정보학: 기능적 농축 분석(GO/KEGG)은 DAVID 및 GSEA를 사용하여 수행되었다. 단백질-단백질 상호작용(PPI) 네트워크는 STRING을 사용하여 구축되었고 Cytoscape를 통해 시각화되었다. 허브 유전자/단백질은 5가지 중심성 알고리즘(MCC, MNC, degree, closeness, stress)을 사용하여 식별되었다.
• 통합 분석: mRNA와 단백질 수준 모두에서 유의하게 변화된 분자를 식별하기 위해 교차 수준 비교를 수행하였다.
• 검증: 정통 AVP 수용체의 낮은 검출률을 고려하여 잠재적인 신호 전달 기전을 탐구하기 위해, 본 연구는 GPR54/KISS1R 길항제인 Kiss234를 활용하였다. 이를 통해 세포 생존력, 세포 내 Ca²⁺, ATP 수준, GnRH 분비 및 프로모터 활성에 대한 GPR54 관련 신호 전달의 기여도를 평가하였다. AVP와 GPR54의 근접성을 조사하기 위해, 궁상핵을 포함하는 쥐의 시상하부 절편에서 근접 결찰 분석(PLA)을 수행하였다.

주요 결과

• 전사체 변화: AVP 처리는 1,094개의 차등 발현 유전자(DEGs)(상향 조절 549개, 하향 조절 545개)를 유도하였다. 주요 상향 조절 유전자는 Nupr1, Slc7a11, Trib3였으며, Rgs4, Tmem119, Ccng1은 하향 조절되었다.
  • 경로: 농축 분석 결과 TNF, PI3K–Akt, MAPK 및 cAMP 신호 전달 경로가 강조되었다. GSEA는 축삭 유도(axon guidance) 및 Wnt 신호 전달 경로를 추가로 식별하였다.
  • 허브 유전자: 네트워크 분석을 통해 Cxcl10, Cxcl1, Cxcl12, Itpkb, Icam1, Ptgs2, Mmp9가 중심 허브 후보임을 확인하였다.
• 단백질체 변화: 본 연구는 2,648개의 차등 발현 단백질(DEPs)(증가 1,362개, 감소 1,286개)을 식별하였다.
  • 경로: DEPs는 대사 경로, 스플라이싱오좀(spliceosome), 그리고 소포체 내 단백질 가공 과정에 농축되어 있었다.
  • 허브 단백질: 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 세포에서 다층적인 조절 반응을 유발하며, 통합적인 전사체 및 단백질체 접근법이 이러한 분자적 반응의 전체 스펙트럼을 포착하는 데 필수적임을 밝히고 있다.