Functional Genomic Analyses Reveal transfer RNA Modification Enzymes as Risk Genes for Bipolar I Disorder and Schizophrenia

기술 요약: 양극성 I형 장애 및 조현병에서의 RNA 변형 효소에 대한 기능적 유전체 분석

문제 제기
전장 유전체 연관 분석(GWAS)을 통해 조현병(SCZ) 및 양극성 I형 장애(BIPI)와 같은 신경정신질량 질환과 관련된 수천 개의 유전적 변이가 식별되었으나, 이러한 광범위한 유전적 신호를 구체적이고 치료 가능한 생물학적 기제로 번역하는 것은 여전히 큰 과제로 남아 있다. 후성유전적 조절은 알려진 취약성 축이지만, "에피트랜스크립톰(epitranscriptome)"—특히 RNA 변형 단백질(RMP)에 의해 매개되는 RNA 변형—의 역할은 신경정신질환 위험의 맥락에서 체계적으로 평가되지 않았다. 100개 이상의 RMP가 확인되었음에도 불구하고, 이들이 정신질환의 유전적 구조에 기여하는 바는 표준 경로 분석을 위한 기능적 유전자 세트로 분류되지 않았다는 점 때문에 여전히 미개척 영역으로 남아 있다.

방법론
저자들은 16가지 신경정신질환의 위험을 결정하는 RMP 유전자를 체계적으로 우선순위화하기 위해 다단계 기능적 유전체 파이프라인을 사용하였다. 본 연구는 6가지 RNA 바이오타입에 걸친 31가지 RNA 변형 유형에 해당하는 123개의 인간 RMP 큐레이션 세트를 활용하였다.

분석 워크플로우는 다음 접근 방식들을 통합하였다:

1. 유전자 기반 연관성 테스트: MAGMA 및 hMAGMA(Hi-C 염색질 상호작용 데이터를 포함)를 사용하여 RMP와 PGC(Psychiatric Genomics Consortium)의 GWAS 요약 통계 간의 유전자 세트 연관성을 테스트하였다.
2. 유전력 추정: GCTB–SBayesR를 적용하여 개별 RMP 유전자 본체(gene body)가 질환 유전력에 기여하는 정도를 추정하였다.
3. 미세 매핑(Fine-Mapping): PolyFun 및 SuSiE를 사용하여 위험 로커스 내의 인과적 변이를 식별하였으며, 이때 SNP당 유전력 추정치를 반영하였다.
4. QTL 통합 분석:
   • SMR (Summary-data-based Mendelian Randomization): GWAS 데이터와 뇌 피질 eQTL 및 sQTL 데이터(BrainMeta)를 통합하여 공동 연관성을 테스트하였고, 이어서 다표현형(pleiotropy)과 연관 불균형(linkage)을 구분하기 위해 HEIDI 테스트를 수행하였다.
   • TWAS (Transcriptome-Wide Association Studies): 유전적으로 예측된 유전자 발현과 질환 위험 사이의 연관성을 테스트하기 위해 GTEx v8 뇌 피질 모델을 이용한 FUSION을 사용하였다.
   • Two-Sample MR: IVW, MR-Egger 및 기타 방법들을 사용하여 방향성 일관성을 검증하였다.
   • 공동 국소화(Colocalization): GWAS와 QTL 신호가 공통의 인과적 변이를 공유하는지 결정하기 위해 베이지안 공동 국소화(coloc)를 적용하였다 (PP.H4 > 0.6).
5. 검증: 차등 발현 분석 및 공변량을 조정한 로지스틱 회귀를 통해 BrainSeq Phase II 코호트(해마, 배외측 전전두엽 피질, 미상핵)의 벌크 RNA-seq 데이터를 사용하여 우선순위가 지정된 후보들을 검증하였다.

주요 기여 및 결과
본 연구는 RNA 사이토신-5 메틸전이효소인 NSUN2를 BIPI와 SCZ 모두에서 재현 가능한 공유 위험 로커스로 식별하였다.

• NSUN2 연관성: 개별 SNP는 전장 유전체 유의 수준에 도달하지 못했으나, NSUN2는 MAGMA 및 hMAGMA 테스트에서 두 질환 모두에 대해 유의미한 유전자 수준의 연관성을 보였다. SMR 및 TWAS 분석은 NSUN2 발현과 질환 위험 사이의 음의 상관관계를 밝혀냈으며, 이는 신경 보호 역할을 시사한다. 공동 국소화 분석은 NSUN2 eQTL과 BIPI 위험 사이의 공유 인과적 변이를 지지하였다.
• 스플라이싱 기제: 연구는 위험이 선택적 스플라이싱(alternative splicing)을 통해 매개된다는 점을 강조하였다. 비부호화 또는 절단된 이소형(isoform)을 생성하는 엑손 건너뛰기(exon skipping)를 유발하는 특정 sQTL은 질환 위험과 양의 상관관계가 있었던 반면, 단백질 코딩 이소형을 촉진하는 사건들은 음의 상관관계가 있었다. 벌크 RNA-seq 검증 결과, SCZ 환자의 뇌(특히 남성 해마 및 여성 배외측 전전두엽 피질)에서 코딩 NSUN2 이소형의 발현이 감소함을 확인하였다.
• 광범위한 tRNA 효소 관여: NSUN2 외에도 다른 tRNA 변형 효소들의 유의미한 연관성을 확인하였다:
  • NSUN6 및 TYW5는 SCZ에 대해 전장 유전체 유의 수준에 도달하였다.
  • TRMT61A 및 QTRT1은 BIPI 또는 SCZ와 유의미한 연관성을 보였다.
  • MRM1은 BIPI와 연관되었다.
  • (COLOCdb를 포함한) 통합 분석은 이러한 유전자들이 여러 QTL 양상(eQTL, sQTL, mQTL)에 걸쳐 연관성이 있음을 확인하였다.
• 우선순위 지정: 가중치 점수 시스템을 통해 NSUN2를 BIPI 및 SCZ에 대해 가장 높은 우선순위 유전자로 선정하였으며, 그 뒤를 SCZ에 대한 NSUN6 및 TYW5가 이었다.

의의 및 주장
본 논문은 BIPI 및 SCZ에 대한 공통 유전적 위험이 RNA 변형 기구, 특히 tRNA 변형의 조절 변이에 수렴된다고 주장하며, tRNA 변형의 조절 장애가 정신질환 감수성에 기여하는 분자 경로임을 지목한다.

주요 결론은 다음과 같다:

• 수렴적 기제: 연구 결과는 tRNA 변형 효소(m5C, m1A, yW, Q 변형의 작성자들)가 더 흔히 연구되는 mRNA m6A 경로와는 구별되는, BIPI 및 SCZ의 병리학적 근저에 있는 수렴적 분자 기제를 나타낸다는 것을 시사한다.
• 공유 위험 유전자로서의 NSUN2: NSUN2는 유전적 위험이 총 유전자 풍부도가 아닌 이소형 특이적 발현 및 스플라이싱의 조절 장애를 통해 작용하는 공유 위험 유전자로 식별되었다.
• 신경 보호 역할: 데이터는 온전한 NSUN2 기능이 신경 보호 효과가 있다는 모델을 지지하며, 그 조절 장애(코딩 이소형 감소를 통해)가 질환 감수성에 기여한다.
• 한계점: 저자들은 현재의 RMP 목록이 포괄적이지 않으며, 뇌 영역 및 세포 유형에 따른 RMP 발현의 이질성이 여전히 과제임을 겸허히 언급하였다. 이들은 본 연구가 후보들을 우선순위화한 것이며, 신경 발생 과정에서의 구체적인 역할을 정의하기 위해서는 동물 모델 및 뇌 오가노이드를 통한 추가 검증이 필요함을 강조하였다.

본 연구는 이 후보들의 역할을 신경 발생 과정에서의 역할을 확인하기 위한 동물 모델 검증의 필요성을 제외하고는 새로운 임상 적용이나 향다 실험적 제안을 제시하지 않는다.