Genome-wide association studies identify genetic determinants of synucleinopathy biomarkers. 

기술 요약: 전전체 게놈 연관 분석을 통한 시뉴클레인 병증 바이오마커의 유전적 결정 요인 식별

문제 제기
파킨슨병(PD), 루이소체 치매(DLB), 다계통 위축증(MSA)을 포함한 $\alpha$-시뉴클레인 병증은 임상적 및 생물학적 이질성으로 인해 진단 및 치료에 상당한 어려움을 초래한다. $\alpha$-시뉴클레인 시드 증폭 분석법(SAA), 아밀로이드-$\beta$(A$\beta$), 타우(tau), 그리고 리소좀 지질과 같은 바이오마커는 조기 진단, 환자 층화 및 임상 시험 농축에 매우 중요하지만, 유전적 수식 인자(genetic modifiers)로 인해 그 활용도가 복잡해진다. 유전적 변이는 바이오마커 수치에 영향을 미칠 수 있으며, 이는 연구 및 임상 환경에서 해석의 편향을 일으킬 가능성이 있다. 예를 들어, 특정 유전자 보유자(예: PRKN 또는 LRRK2)는 질병 상태와 무관하게 변화된 바이오마커 프로파일을 보일 수 있다. 따라서 연구의 정확한 해석을 보장하고 중개 연구의 관련성을 높이기 위해 바이오마커 수치를 조절하는 공통 유전 변이를 체계적으로 식별할 필요가 있다.

연구 방법론
본 연구는 산발성 파킨슨병 사례, 건강한 대조군, SWEDD(도파민 결핍 증거가 없는 스캔 대상자), 그리고 전구기(prodromal) 개인을 포함하여 유럽 혈통의 581명으로 구성된 파킨슨병 진행 마커 이니셔티브(PPMI) 코호트 데이터를 활용하였다.

• 바이오마커 분석: 뇌척수액(CSF), 혈장, 소변에서 63개의 후보 바이오마커를 분석하였다. 여기에는 $\alpha$-시뉴클레인, A$\beta$, 타우 이소형, 신경필라멘트 경쇄(NfL), 그리고 다양한 지질 대사체(bis(monoacylglycerol)phosphate [BMP] 이소형, 세라마이드, 라크토실세라마이드, 스핑고미엘린)가 포함되었다.
• 유전 분석: 품질 관리 및 람다 통계량이 과잉된 마커를 제외한 후, 59개 바이오마커에 대한 전전체 게놈 연관 분석(GWAS)을 수행하였다. 분석에는 TOPMed 참조 패널을 사용하여 임퓨테이션(imputation)된 두 개의 유전형 데이터셋(NeuroX 및 OmniExpress Exome+)이 사용되었다.
• 통계적 접근법: 선형 회귀를 사용하여 공통 변이와 바이오마커 수치 사이의 연관성을 테스트하였으며, 연령, 성별, 질병 상태 및 상위 5개 주성분을 보정하였다. 두 유전형 어레이 모두에서 측정된 바이오마커에 대해 메타 분석을 수행하였다. 독립적인 신호를 식별하기 위해 조건부 및 결합 분석(GCTA-COJO)이 사용되었다. 또한, PD 상태에 대한 바이오마커의 예측력을 평가하기 위해 로지스틱 회귀 및 AUC(Area Under the Curve) 분석을 수행하였다.
• 재현: 선택된 결과는 컬럼비아 대학교의 외부 코호트에서 검증되었다.

주요 기여 및 결과
본 연구는 확립된 것뿐만 아니라 새로운 시뉴클레인 병증 바이오마커의 유전적 조절 인자를 식별하였다:

1. 소변 BMP 이소형 및 LRRK2/GBA1:

   • 소변 BMP 이소형(특히 2,2'-di-22:6, total di-18:1, total di-22:6)과 LRRK2 p.G2019S 변이(양의 상관관계) 및 GBA1 p.N370S 변이(음의 상관관계) 사이의 강력한 연관성을 확인하였다.
   • 이러한 연관성은 외부 코호트에서도 재현되었다.
   • 결정적으로, BMP 수치가 이러한 유전형에 의해 구동됨에도 불구하고, 본 연구에서는 특정 유전자 보유자 그룹 내에서 PD 사례와 대조군 간의 BMP 수치에 유의미한 차이가 없음을 발견하였다. 이는 BMP가 PD의 단독 진단 바이오마커라기보다는 LRRK2 경로의 지표로서 더 잘 기능함을 시사한다.

2. CSF A$\beta$ 및 APOE:

   • CSF A$\beta$ 수치와 APOE $\epsilon$4 대립유전자(rs429358) 사이에 강력한 음의 연관성이 확인되었으며, 이는 아밀로이드 조절에서 APOE의 핵심적인 역할을 뒷받침한다. 이 신호는 A$\beta$/tau 비율에서 관찰된 연관성도 주도하였다.

3. CSF 세라마이드 및 새로운 로커스(Loci):

   • CSF 세라마이드 이소형(C22, C23, C24)과 MCF2L2 로커스 사이의 새로운 연관성이 탐지되었다. 주목할 점은 MCF2L2가 알려진 PD 위험 로커스(MCCC1;LAMP3;MCF2L2)의 일부이지만, 본 연구의 주요 변이는 일차적인 PD GWAS 신호와 강한 연관 불평형(linkage disequilibrium) 상태에 있지 않았다는 것이다.
   • GMNN 로커스 또한 CSF 세라마이드 수치와 연관되어 있었으며, 이는 $\alpha$-시뉴클레인 병증 바이오마커와의 이전에 보고되지 않은 연결을 나타낸다.

4. CSF 타우 및 TP63:

   • CSF 총 타우(t-tau) 수치와 TP63 근처의 인터제닉(intergenic) 로커스 사이에 유의미한 음의 연관성이 발견되었다. 인산화된 타우(p-tau)에서도 유사한 경향이 관찰되었다. 이 영역은 이전에 타우 생물학에 관여하는 것으로 알려진 바가 없었다.

5. PD 및 DLB 위험 로커스와 리소좀 지질:

   • 여러 PD 위험 로커스가 리소좀 지질 바이오마커와 연관되어 있어, 경로 수준의 수렴을 보여주었다:
     • MAPT 로커스: 혈장 라크토실세라마이드(GL2) 수치 증가와 연관됨.
     • SIPA1L2: 혈장 GL2 C23 증가와 연관됨.
     • MCCC1/LAMP3: 혈장 스핑고미엘린(SM) 수치와 연관됨.
     • RAB29: 혈장 세라마이드 수치와 연관됨.
   • DLB 분석에서는 APOE 변이가 CSF A$\beta$ 및 A$\beta/tau$ 비율의 감소와 연관되어 있었다.

의의 및 주장
본 논문은 유전적 배경이 잠재적인 $\alpha$-시뉴클레인 병증 바이오마커의 수치를 유의미하게 형성한다고 결론짓는다. 저자들은 데이터의 오해석을 방지하고 환자 층화를 개선하기 위해 바이오마커 기반 연구, 임상 시험 및 치료 개발 시 이러한 유전적 영향이 반드시 고려되어야 한다고 주장한다.

구체적으로, 본 연구는 다음을 강조한다:

• BMP와 같은 바이오마커는 LRRK2 및 GBA1 유전형에 의해 견고하게 구동되지만, 질병 상태의 차이를 반영하지 않을 수 있으므로, 이를 활용하는 방식(진단 도구가 아닌 경로 지표로서의 활용)에 대한 전환이 필요하다.
• TP63 및 GMNN과 같은 새로운 로커스는 신경 퇴행과 관련된 이전에 인식되지 않은 생물학적 경로를 나타낼 수 있다.
• 정밀 의료 접근법이 PD 및 DLB 분야에서 발전함에 따라, 신뢰할 수 있고 실행 가능한 바이오마커를 개발하기 위해서는 유전적 맥락의 통합이 필수적이다.

저자들은 유럽 혈통에 국한된 연구 대상, 횡단적 데이터의 특성(인과적 추론 불가), 희귀 변이 연관성에 대한 재현을 위한 제한된 샘령 규모 등의 한계를 인정하였다. 이들은 자신의 연구 결과를 검증하고 정교화하기 위해 포괄적인 유전 및 바이오마커 데이터를 갖춘 더 크고 다양한 코호트가 필요함을 강조하였다.