A cross-platform Mendelian randomization study of mitochondrial DNA copy number across psychiatric disorders

기술 요약: 정신질환에 걸친 미토콘드리아 DNA 복제수(mtDNA-CN)에 관한 교차 플랫폼 멘델 무작위 분석 연구

문제 제기
미토콘드리아 DNA 복제수(mtDNA-CN)는 미토콘드리아 생물학을 나타내는 접근 가능한 말초 대리 지표 역할을 하지만, 정신적 취약성과의 관계를 해석하기는 매우 어렵다. 말초 혈액 mtDNA-CN은 분석 플랫폼의 차이, 조직 맥락, 그리고 혈액 세포 구성에 의해 혼란을 겪는다. 관찰 연구들은 자폐 스펙트럼 장애(ASD), 주의력 결핍 과잉 행동 장애(ADHD), 조현병, 양극성 장애와 같은 질환에서 변동된 mtDNA-CN을 보고해 왔으나, 그 결과는 이질적이다. 이러한 연관성이 상위 유전적 취약성을 반영하는 것인지, 아니면 질병의 경과, 치료, 염증 또는 세포 구성의 하위 상관관계인지는 여전히 불분명하다. 또한, 이전의 멘델 무작위 분석(MR) 연구들은 노출 패널의 차이와 여러 GWAS 간의 상당한 샘플 중복(특히 UK Biobank 내에서의 중복)을 고려하지 못함으로 인해 일관되지 않은 결과를 도출하였다.

방법론
본 연구는 두 샘플 멘델 무작위 분석 프레임워크를 사용하여 광범한 정신적 표현형에 걸쳐 유전적으로 대리된 mtDNA-CN을 평가하였다. 분석 파이프라인은 다음과 같은 핵심 요소들을 포함한다:

• 노출 데이터: 세 가지 분석 접근 방식을 아우르는 네 개의 부분적으로 중복된 유럽 혈통 GWAS 노출 패널이 활용되었다:
  1. qPCR 정량화 (Longchamps et al., N = 465,809).
  2. SNP 어레이 강도 추정 (Chong et al., N = 383,476).
  3. 전장 유전체 시퀀싱(WGS) 유래 mtDNA-CN (조정되지 않은 버전 및 혈액 세포 구성 조정 버전 각각 N = 155,998, Gupta et al.).
• 결과 데이터: 본 연구는 다섯 가지 주요 질환(ASD, ADHD, 양극성 장애, 주요 우울 장애, 조현병), 아홉 가지 이차적 결과(불안, 강박 장애, PTSD, 항우울제 반응 포함), 그리고 14개 질환의 게놈 구조 방정식 모델로부터 도출된 여섯 가지 교차 질환 요인을 평가하였다. 신장과 BMI는 음성 대조군 및 생물학적 대조물로 사용되었다.
• 도구 선택: 게놈 전체 유의 수준($p < 5 \times 10^{-8}$)의 변이들이 도구로 선택되었으며, 방향성을 보장하기 위해 LD 클럼핑($r^2 < 0.001$)과 Steiger 필터링이 적용되었다.
• 통계적 합성: 본 연구의 핵심적인 방법론적 혁신은 네 개의 노출 패널을 독립적인 반복이 아닌 상관관계가 있는 것으로 취급한 것이다. 이는 공유된 UK Biobank 참가자들로 인한 중복 때문이다. 패널 수준의 추정치는 경험적으로 관찰된 평균 패널 간 상관관계($\rho = 0.70$)를 기반으로 한 설계 효과(design-effect) 보정을 통해 합성되었으며, 공유 샘플 의존성을 고려하여 표준 오차를 팽창시켰다.
• 민감도 및 기전 분석: 역방향 MR, 샘플 중복 보정을 위한 MRlap, 유전자 세트 농축(MitoCarta3.0), 그리고 13개 GTEx 뇌 조직에 대한 탐색적 뇌 eQTL 공국소화(colocalization)를 포함하였다.

주요 기여

1. 교차 플랫폼 합성: 본 연구는 네 가지 서로 다른 mtDNA-CN GWAS 패널(qPCR, SNP-array, WGS)을 체계적으로 통합하고, 중복된 샘플로 인해 발생하는 상관관계를 명시적으로 모델링한 최초의 연구이다.
2. 설계 효과 보정: 중복된 패널을 독립적인 것으로 취급하는 대신 설계 효과 보정을 적용함으로써, 본 연구는 mtDNA-CN의 인과적 효과에 대해 보다 보수적이고 통 статисти적으로 엄격한 추론을 제공한다.
3. 요인 수준 분석: 분석은 개별 질환을 넘어 교차 질환 요인(예: 신경발달, 조현-양극성, 일반 정신병리)으로 확장되어 범진단적 관점을 제공하였다.
4. 기전 주석 달기: 본 연구는 MR을 유전자 세트 농축 및 뇌 eQTL 공국소화와 결합하여 가설 생성용 후보 유전자(예: RECELL4, SUOX, MSRA)를 명명하였으나, 매개 관계를 확립하는 데 따르는 한계점 또한 인정하였다.

결과

• 주요 발견: 유전적으로 대리된 높은 mtDNA-CN은 자폐 스펙트럼 장애(ASD)에 대한 낮은 취약성과 견고하게 연관되어 있었다. 이 역의 연관성($\beta = -0.182$)은 설계 효과 보정 및 Benjamini-Hochberg 허위 발견율(BH-FDR) 보정 후에도 유의미하게 유지되었다 ($q = 0.030$).
• 신경발달 요인: ADHD와 ASD-ADHD 신경발달 요인(F3)은 방향적으로 일치하는 음의 추정치를 보였다. 그러나 이들은 엄격한 BH-FDR 보정을 통과하지 못했다 (ADHD 보정 $p = 0.173$; F3 보정 $q = 0.062$). 저자들은 F3가 ASD와 ADHD의 복합 로딩이므로, F3 결과를 ASD 결과의 독립적인 재현으로 해석해서는 안 된다고 언급하였다.
• 무효 결과: 조현병, 양극성 장애, 주요 우울 장애, 일반 정신병리(pFactor) 또는 기타 이차적 결과(불안, 강박 장애, PTSD 등)에서는 견고한 연관성이 관찰되지 않았다.
• 민감도 분석: 역방향 MR은 정신적 취약성이 mtDNA-CN에 영향을 미친다는 증거를 제공하지 않았다. 음성 대조군(신장)은 연관성을 보이지 않았으며, BMI는 예상된 이질성을 보였다. UK Biobank 중복 제거 결과와 MRlap을 사용한 민감도 분석은 주요 발견의 방향성을 확인하였다.
• 생물학적 주석: 유전자 세트 농축은 네 가지 노출 패널 중 세 곳에서 mtDNA-CN 도구의 미토콘드리아 관련성을 뒷받침하였다. 그러나 결과 GWAS(ASD, ADHD, F3)는 그에 상응하는 양의 미토콘드리아 농축을 보여주지 않았다. 공국소화 분석은 몇몇 후보 유전자를 지목했으나 매개를 확립하지는 못했다.

의의 및 주장
본 논문은 연구 결과가 유전적으로 대리된 mtDNA-CN과 신경발달 취약성 사이의 특정 ASD 중심 역의 연관성을 지지한다는 점을 주장한다. 저자들은 이 연관성이 분석 플랫폼 및 공유 샘플 상관관계에는 견고하지만, 혈액 세포 구성 조정에는 민감하다는 점을 강조한다. 즉, 이 연관성은 qPCR, SNP-array, 조정되지 않은 WGS 패널에서는 유의미했으나, 혈액 세포 구성이 조정된 WGS 패널에서는 유의성에 도달하지 못했다.

저자들은 이러한 결과들을 "확정적이라기보다 가설을 강화하는 것"으로 규정한다. 그들은 데이터가 신경발달 장애에서 미토콘드리아 경로의 잠재적 역할을 시사하지만, 더 넓은 기전적 해석(예: 직접적인 인과 경로 대 뉴클레오타이드 대사 또는 세포 증식을 통한 상관된 다표현형성)은 추가적인 검증이 필요하다고 주장한다. 본 연구는 광범위한 범정신병적 mtDNA-CN 메커니즘을 명시적으로 거부하는데, 이는 정신병적, 정동적 또는 일반 정신병리 요인에서 견고한 효과가 발견되지 않았기 때문이다. 본 연구는 미토콘드리아 경로가 신경발달 장애에 대한 기전적 조사를 위해 지목되었으나, 조직 특이적 검증과 다변량 MR이 인과관계와 기전을 확립하기 위해 필요하다는 결론을 내린다.