Interactome modules underlie shared spatiotemporal patterns in specific psychiatric and neurodegenerative disorders

기술 요약: 특정 정신 질환 및 신경 퇴행성 질환에서 공유된 시공간적 패턴을 뒷받받하는 인터랙톰 모듈

문제 제기
정신 의학적 및 신경학적 임상 증후군은 증상 측면에서 상당한 중첩(예: 강박, 정동, 정신병 및 운동 클러스터)을 보이며, 종종 발병 연령 또한 겹치는 양상을 보인다. 유전적 상관관계 연구를 통해 이러한 중첩을 정량화할 수는 있었으나, 위험 유전자, 특히 그들의 1차 상호작는(first-order interactors)의 전사체 및 인터랙톰 패턴에 대한 상세하고 가설 생성적인 연구는 부족한 실정이다. 기존의 분류법은 기저의 생물학적 복잡성을 반영하지 못하는 경우가 많으며, 진단적 경계(예: 강박장애(OCD)와 파킨슨병(PD)과 같은 정신 질환과 신경 퇴행성 질환 사이)를 가로지르는 위험 유전자의 메커니즘은 여전히 불분명하다.

연구 방법론
본 연구는 임상적 중첩이 존재하는 다섯 가지 증후군인 강박장애(OCD), 헌팅턴병(HD), 파킨슨병(PD), 조현병(SZ), 양극성 장애(BD)를 분석하기 위해 전사체 및 인터랙톰 프레임워크를 통합하였다.

1. 유전자 컴파일: GWAS 카탈로그에서 위험 유전자를 수집하였다. 증후군 간의 유전자 수 차이를 해결하기 위해, 저자들은 PD, SZ, BD에 대해 높은 유전적 부하(예: 가족력, 조기 발병, 치료 저항성)를 반영하는 로커스(loci)를 포함하도록 필터링하였다. 이를 통해 선별된 173개의 질환 관련 유전자 세트(HD 59개, PD 29개, OCD 34개, SZ 31개, BD 20개)를 도출하였다.
2. 네트워크 구축: 첫 번째 단계의 단백질-단백질 상호작용(PPI)을 BioGRID 및 HPRD에서 검색하여, 위험 유전자와 그들의 즉각적인 상호작용 유전자를 포함하는 증후군 특이적 네트워크를 구축하였다. 이를 통해 1,862개의 고유 유전자로 구성된 결합 네트워크를 생성하였다.
3. 연결성 분석: 차수 보존 몬테카를로 시뮬레이션(1,000회 반복)을 사용하여, 질환 유전자와 그 상호작용 유전자들이 무작위 기대치보다 유의미하게 더 응집된 서브 네트워크를 형성하는지 테스트하였다. 견고성 검증을 위해 상위 1% 및 5%의 허브 노드를 제거하며 실험을 진행하였다.
4. 시공간적 발현 프로파일링: BrainSpan Atlas(524개 샘플, 26개 뇌 구조, 31개 시간 지점)와 GTEx의 발현 데이터를 사용하였다. 계층적 클러스터링을 통해 공간적 클러스터를 식별하였다. 모듈 식별은 Cytoscape 내의 MCODE 알고리즘을 사용하여 수행되었다.
5. 기능 및 지역적 농축: 유전자 온톨로지(GO), KEGG, Reactome 분석을 통해 생물학적 과정을 식별하였다. 하이퍼기하 검정(Hypergeometric tests)을 통해 특정 뇌 영역 및 발달 단계에서의 모듈 농축도를 결정하였다.

주요 결과

• 인터랙톰 연결성: 질환 유전자 세트 자체는 유의미한 내부 집합 연결성을 보이지 않았다. 그러나 1차 상호작용 유전자를 포함했을 때, 네트워크는 무작위 기대치를 크게 상회하는 광범위한 연결성(34,702개 엣지)을 나타냈다($p < 0.001$). 이러한 농축은 상위 허브를 제거한 후에도 지속되었으며, 이는 연결성이 글로벌 네트워크 구조가 아닌 분산된 질환 관련 상호작용에 의해 구동됨을 시사한다.
• 공간적 클러스터링: 가용한 발현 데이터가 있는 1,826개 네트워크 유전자를 계층적 클러스터링한 결과, 두 개의 공간적 클러스터로 분리되었다:
  • 공간적 클러스터 1: 태아 구조(종창부, 횡문 입술 등)에서 발현되는 유전자에 농축되었다. 이 클러스터는 OCD 네트워크에 대한 유의미한 농축을 보였다.
  • 공간적 클러스터 2: 성인 뇌 구조(해마, 선조체, 운동 피질)에서 발현되는 유전자에 농축되었다. 이 클러스터는 HD 네트워크에 대한 유의미한 농축을 보였다.
  • 주목할 점은, 늦은 발병 증후군임에도 불구하고 PD 네트워크가 OCD 네트워크의 태아기 시그니처와 유사한 농축 프로파일을 보였다는 것이다.
• 모듈 수준 분석: 다섯 가지 증후군 네트워크 전반에서 44개의 위상적 모듈을 식별하였다. 특정 모듈들이 관찰된 공간적 농축을 주도하였다:
  • OCD-1, PD-3, PD-4는 공간적 클러스터 1(태아기)에 농축되었다.
  • HD-1, HD-3, SZ-5는 공간적 클러스터 2(성인기)에 농축되었다.
• 공유 메커니즘 (PD-OCD 축): OCD-1, PD-3, PD-4 모듈은 인터랙톰 내에서 근접한 위치에 있음이 확인되었다. 이들은 리보솜 관련 과정(예: 리보핵단백질 복합체 생합성, RNA 대사)에 대한 유의미한 농축을 공유하였다. 또한 이 모듈들은 초기 발달 단계 동안 발달 중인 배외측 전전두엽 피질에서 발현되는 유전자들에 대한 공동 농축을 보였다.
• 검증: 몬테카를로 분석을 통해 결합된 OCD-1, PD-3, PD-4 모듈 내 리보솜 관련 유전자의 출현이 비무작위적임을 확인하였다($p = 0.002$). 또한, 이 모듈들에 속한 24개의 리보솜 관련 유전자 중 19개가 OCD와 PD 모두에 대해 공개된 사후 뇌 조직 데이터셋에서 차별 발현되었다.

의의 및 주장
본 논문은 유전적 위험이 뇌 영역 및 발달 시간대를 가로질러 조직된, 분산되어 있으면서도 일관된 모듈을 통해 작용한다는 네트워크 기반의 뇌 질환 모델을 지지한다.

• 메커니즘적 통찰: 본 연구는 질환 간의 중첩이 전체 증후군 네트워크 전체에서 균일하게 일어나는 것이 아니라, 고도로 연결된 특정 단백질 모듈에 의해 구동된다는 점을 시사한다.
• PD-OCD 연결 고리: 본 연구의 주요 발견은 파킨슨병(PD)과 강박장애(OCD) 사이의 잠재적인 메커니즘적 연결 고리를 식별한 것이다. 발병 연령의 차이에도 불구하고, 두 질환은 전전두엽 피질에서의 태아기 발현 시그니처와 공통된 리보솜 관련 과정의 관여를 공유한다. 저자들은 이를 검증 가능한 가설로 제시한다: 즉, 신경 발달 과정 중 리보솜 유전자의 교란이 두 질환 모두에 대한 감수성을 변화시킬 수 있다는 것이다.
• 방법론적 기여: 본 연구는 PPI 네트워크와 그 구성 모듈을 (전체 네트워크 분석이 아닌) 지역적 발현 프로파일과 연결하는 것이, 진단적 경계에 의해 가려져 있던 특정 생물학적 테마를 드러낼 수 있음을 보여준다.

저자들은 이러한 발견이 공유된 테마를 포착하고 검증 가능한 가설을 생성하는 데 도움이 된다고 결론지었으나, 인터랙톰 모델의 정적인 특성과 임상적 공존성을 확증하기 위한 대규모 환자 데이터의 필요성 등 한계점도 인정하였다.