ASSOCIATION ANALYSIS OF MITOCHONDRIAL HETEROPLASMIC VARIANTS AND CARDIOMETABOLIC TRAITS

기술 요약: 미토콘드리아 이형접합 변이와 심혈관 대사 형질의 연관성 분석

문제 제기
미토콘드리아 기능 장애는 심혈관 및 대사 질환의 병인에 관여하는 것으로 알려져 있으나, 흔한 복합 심혈관 대사 질환(CMD)의 맥락에서 미토콘드리아 이형접합성(heteroplasmy)—즉, 세포 내에 혼합된 미토콘드리아 DNA(mtDNA) 분자 집단이 공존하는 현상—의 구체적인 역할은 여전히 명확히 밝혀지지 않았다. 드문 단일 염기 변이가 심각한 미토콘드리아 질환을 유발한다는 것은 잘 알려져 있지만, 비만, 고혈압, 당뇨병, 고지혈증과 같은 흔한 질환에서 드문 이형접합 변이가 기여하는 바를 조사하는 것은 대규모 인간 코호트의 심층 시퀀싱 데이터 부족으로 인해 어려움이 있었다. 또한, 이형접합 변이의 98%는 매우 드물며(단 한 명 또는 소수의 개인에게만 존재) 낮은 변이 대립유전자 분율(VAF)을 가지므로, 이를 탐지하기 위한 특화된 통계적 프레임워크가 필요하다.

방법론
본 연구는 6개의 다인종 TOPMed 코호트(ARIC, CARDIA, CHS, FHS, JHS, MESA)에 참여한 16,882명의 심층 전전체 게놈 시퀀싱(WGS) 데이터를 활용하였다. 분석은 8가지 심혈관 대사 형질인 체질량 지수(BMI), 비만, 수축기 혈압(SBP), 고혈려, 혈당(BG), 당뇨병, 저밀도 지질단백질(LDL), 고지혈증에 초점을 맞추었다.

주요 방법론적 구성 요소는 다음과 같다:

• 이형접합성 식별: 수정된 캠브리지 참조 서열(rCRS)을 사용하여 MToolBox를 통해 이형접합 변이를 식별하였으며, 핵 미토콘드리아 DNA 단편(NUMTs)의 영향을 최소화하기 위해 엄격한 VAF 임계값인 5%에서 95%를 적용하였다. VAF가 95%를 초과하는 변이는 동형접합(homoplasmic)으로 분류하였다.
• 유전적 코딩: 이형접합 변이의 생물학적 측면을 포착하기 위해 세 가지 서로 다른 코딩 정의를 적용하였다:
  1. 이진 코딩(Binary Coding): 지표 함수(5% ≤ VAF ≤ 95%이면 1, 그렇지 않으면 0).
  2. VAF 가중 코딩(VAF-weighted Coding): 실제 변이 대립유전자 분율을 포함.
  3. Mito-score 코딩(Mito-score Coding): 국소적 치환 불내성에서 유도된 기능적 제약 점수(Mito-score)를 포함.
• 통계적 프레임워크: 16개의 mtDNA 유전자에 대해 네 가지 유전자 기반 연관성 테스트(Burden, SKAT, SKAT-O, ACAT-O)를 수행하였다. 분석은 아프리카계 미국인(AA)과 유럽계 미국인(EA) 조상 그룹에 대해 각각 수행된 후 메타 분석되었다. 공변량으로는 연령, 성별, BMI(비BMI 형질의 경우)가 포함되었으며, 이형접합 부담이 채혈 연도와 상관관계가 있는 특정 코호트(FHS, JHS)의 경우 방문 연도를 추가로 조정하였다.
• 유의성 임계값: 16개 mtDNA 유전자 영역에 대해 본페로니 교정(Bonferroni correction)을 적용하였다 ($P \le 0.00313$).

주요 결과
본 연구는 본페로니 교정 후 12개의 유의미한 유전자-형질 연관성을 식별하였으며, 세 가지 코딩 정의 모두에서 일관된 효과 방향이 관찰되었다.

• 가장 강력한 연관성: 가장 유의미한 발견은 EA 참여자에서 MT-CO1 유전자(Cytochrome C Oxidase I 인코딩) 내의 이형접합 변이와 고지혈증 사이의 관계였다. Mito-score 코딩(정의 3)과 Burden 테스트를 사용했을 때, 오즈비(OR)는 0.28 (95% CI: 0.17–0.46, $P = 3.4 \times 10^{-7}$)이었다.
• 기타 유의미한 연관성: 추가적인 유의미한 연관성은 다음과 같이 검출되었다:
  • 고지혈증: EA 참여자의 전전체 게놈 이형접합 부담 (OR = 0.79, $P = 0.00082$).
  • 약물 조정 LDL: EA 참여자의 MT-CO1 변이 (Burden 테스트 $\beta = -13.4$, $P = 6.0 \times 10^{-4}$).
  • 약물 조정 SBP: EA 참여자의 MT-CO3 변이 (SKAT $P = 0.0005 \text{ (또는 } 0.00051)$).
  • 비만: AA 참여자의 MT-CO2 변이 (Burden 테스트 OR = 0.57, $P = 0.0015$).
  • BMI, 혈당, 당뇨병에 대해서도 다양한 유전자 및 조상 그룹에서 연관성이 발견되었다.
• 조상 차이: EA 하위 집단의 더 작은 표본 크기와 질병 유병률의 차이로 인해, AA 참여자에 비해 EA 참여자에서 더 적은 수의 유전자-형질 연관성이 확인되었다. 이형접합성 관련 형질 간의 중첩은 AA와 EA 그룹 사이에서 발견되지 않았다.
• 테스트 성능: 결과는 서로 다른 유전자-형질 연관성이 서로 다른 분포의 인과 변이에 의존함을 시사했다. 즉, 일부는 일관된 효과를 가진 다수의 드문 변이에 의해 구동되었고(Burden 테스트에 의해 검출), 다른 일부는 더 적은 비율의 변이 또는 혼합된 효과 방향에 의해 구동되었다(SKAT에 의해 검출).

의의 및 주장
저자들은 본 연구가 대규모 다인종 WGS 데이터셋을 사용하여 전체 mtDNA 게놈에 걸친 드문 이형접합 변이를 체계적으로 평가한 첫 번째 연구라고 주장한다. 연구 결과는 특정 mtDNA 유전자 내의 이형접합 변이가 심혈관 대사 형질에 기여함을 강조하며, 이는 CMD의 병태생리에 대한 새로운 통찰을 제공한다.

본 연구는 여러 코딩 정의(이진, VAF 가중, Mito-score)를 사용하는 것이 유사한 결과를 낳았음을 강조하며, 이는 발견의 신뢰성을 강화한다. 저자들은 양(+)과 음(-)의 연관성이 모두 발견된 것은 mtDNA 게놈 내에 상반된 효과 방향을 가진 이형접합 변이가 존재함을 시사한다고 언급하였다. 본 연구는 이러한 통계적 연관성을 확립하였으나, 저자들은 생물학적 기전은 아직 연구 과제로 남아 있다고 겸허히 밝히며, 향-후 연구에서 이형접합성의 역할을 완전히 밝히기 위해 mtDNA-nDNA 상호작용 및 인과적 매개 경로를 탐구해야 한다고 제언하였다.