Lysosomal polygenic risk score in Parkinson's disease across populations.

기술 요약: 다양한 인구 집단 간 파킨슨병에 대한 리소좀 다유전성 위험 점수(Lysosomal Polygenic Risk Score)

문제 제기
파킨슨병(PD)은 유전적, 환경적, 노화 요인이 복합적으로 작용하는 임상적으로 이질적인 신경퇴행성 질환이다. 전전체 연관 분석(GWAS)을 통해 수많은 위험 변이가 식별되었으나, 이는 유전력의 극히 일부만을 설명한다. 리소좀 기능 장애는 특히 GBA1 유전자와 관련하여 파킨슨병의 잘 알려진 병리 기전이다. 그러나 특정 병원성 GBA1 변이를 제외한 더 넓은 범위의 리소좀 경로 내 공통 유전적 변이가 파킨슨병 위험 및 임상 표현형에 기여하는 정도는, 특히 다양한 조상 집단에 걸쳐 아직 충분히 탐구되지 않았다. 기존의 다유전성 위험 점수(PRS)는 조상 간 전이성이 부족한 경우가 많으며, 리소좀 기능에 대한 경로 특이적 PRS 모델 또한 대규모 다인종 코호트에서 종합적으로 평가된 바 없다.

연구 방법론
저자들은 리소좀 다유전성 위험 점수(lyso-PRS)를 구축하고 이를 Global Parkinson's Genetics Program (GP2) Release 11 데이터를 사용하여 평가하였다.

• 코호트: 연구는 10개의 뚜렷한 조상 집단(유럽인, 아프리카계 미국인, 아프리카인, 아슈케나지 유대인, 라틴계/아메리카 원주민, 중앙아시아인, 동아시아인, 중동인, 남아시아인, 복합 혼혈 이력)에 걸친 32,054명의 특발성 파킨슨병(iPD) 사례와 18,188명의 건강한 대조군(HC)을 분석하였다.
• 유전 데이터: 유전형은 Neurobooster 어레이에서 유래하였으며, Genools를 사용하여 임퓨테이션(imputation) 및 품질 관리를 수행하였다.
• 표현형 정의:
  • iPD: 확립된 파킨슨병 유전자(GBA1, LRRK2, SNCA 등 포함)의 병원성 변이가 있는 개인을 제외하여 정의함.
  • GBA1-PD: 병원성/가능성 있는 병원성 GBA1 변이를 가진 개인.
  • 비발현 보인자(NMC): GBA1 변이는 있으나 파킨슨병 진단을 받지 않은 개인.
• PRS 구축:
  • 리소좀 유전자는 Molecular Signatures Database(MSigDB)를 사용하여 정의하였으며, 총 519개의 고유 유전자를 포함하였다. GBA1과 LRRK2는 나머지 리소좀 경로의 다유전적 부담을 분리하기 위해 유전자 세트에서 제외하였다.
  • 대립유전자 가중치는 2025년 유럽 파킨슨병 GWAS 요약 통계량으로부터 도출되었다.
  • 미세 대립유전자 빈도(MAF)가 0.01 미만인 변이는 제외하였다. 필터링 후 217,053개의 SNP가 유지되었다.
  • 점수는 LD 클럼핑($r^2 = 0.1$, 250kb)을 적용한 PRSice-2를 사용하여 계산되었다.
• 통계 분석:
  • 성능은 AUC, 오즈비(Odds Ratio, OR), 그리고 균형 정확도(balanced accuracy)를 통해 평가되었다.
  • 임상 결과(발병 연령, UPDRS 점수, 인지 기능)와의 연관성은 성별, 연령, 주성분(PC)을 조정한 회귀 모델을 통해 테스트되었다.
  • 유럽 코호트는 훈련 세트(European-GP2)와 독립 검증 세트(European-TUEPAC)로 나누었다. 비유럽 코호트는 제한된 표본 크기로 인해 기본 p-값 임계값을 사용하였다.

주요 결과

1. 조상별 성능:

   • lyso-PRS는 유럽 인구 집단에서 가장 강력한 예측 성능을 보였다. 유럽-GP2 훈련 코호트에서 AUC는 0.602 ($p = 2.02 \times 10^{-172}$)였고, 유럽-TUEPAC 검증 코호트에서는 0.632 ($p = 5.34 \times 10^{-20}$)였다.
   • 비유럽 인구 집단에서는 성능이 현저히 낮았다 (AUC 범위: 0.51–0.59). 가장 낮은 성능은 아프리카계 인구에서 관찰되었다 (AUC = 0.509).
   • 비유럽 집단에서 lyso-PRS 점수는 일반적으로 건강한 대조군보다 iPD에서 높게 나타났으나, 아프리카계 미국인 및 아프리카인 코호트에서는 유의성에 도달하지 못했다.

2. 유럽인에서의 PD 아형 구분:

   • lyso-PRS는 GBA1-PD를 건강한 대조군(AUC = 0.621) 및 비발현 보인자(AUC = 0.621) 모두로부터 성공적으로 구분하였다.
   • 그러나 이 점수는 GBA1-PD와 특발성 파킨슨병(iPD) 사이에서 약한 판별력을 보였다 (AUC = 0.515). 이는 해당 점수가 GBA1 변이 상태를 구체적으로 포착하기보다는 일반적인 파킨슨병 위험을 포착함을 시사한다.

3. 임상적 연관성 (유럽 코호트):

   • 발병 연령: 높은 lyso-PRS는 유럽-GP2 코호트에서 더 이른 발병 연령과 관련이 있었으나 ($\beta = -0.21$, $p = 0.022$), 효과 방향은 일치했음에도 불구하고 유럽-TUEPAC 코호트에서는 재현되지 않았다 ($p = 0.14$).
   • 운동 증상: 높은 lyso-PRS는 훈련 코호트에서 일상생활 운동 경험(UPDRS-II) 및 운동 심각도(UPDRS-III)의 악화와 관련이 있었으나, 검증 코호트에서는 이러한 결과가 재현되지 않았다.
   • 운동 합병증: 높은 lyso-PRS와 UPDRS-IV(운동 합병증) 사이의 일관된 양의 상관관계가 훈련 코호트 ($\beta = 0.32$, $p = 0.0088$)와 검증 코호트 ($\beta = 0.22$, $p = 0.045$) 모두에서 발견되었다.
   • 인지 기능: 검증 코호트에서 높은 lyso-PRS는 더 나쁜 인지 수행 능력과 관련이 있었으나 ($\beta = -0.17$, $p = 3.58 \times 10^{-3}$), 훈련 코호트에서는 유의한 연관성이 발견되지 않았다.
   • 성별: 두 코호트 모두에서 lyso-PRS와 성별 사이의 유의한 연관성은 발견되지 않았다.

의의 및 주장
저자들은 본 연구가 10개 인구 집단에 대한 가장 광범위하고 첫 번째인 리소좀 다유전성 위험의 종합적 평가라고 주장한다. 주요 의의는 다음과 같다:

• 경로 특이적 위험: (GBA1 변이와 독립적인) 리소좀 경로 내의 공통 유전적 변이는 파킨슨병 감수성에 기여하며, GBA1-PD를 비발현 보인자와 구분하여 위험을 계층화할 수 있다.
• 가산적 위험 모델: 이러한 결과는 GBA1 변이와 더 넓은 범위의 리소좀 다유전적 부담이 독립적으로 파킨슨병 위험을 증가시킨다는 가산적 효과를 시사한다.
• 임상적 관련성: 리소좀 다유전적 부담은 특정 임상 특징, 특히 운동 합병증(UPDRS-IV)에 영향을 미칠 수 있으며, 이는 리소좀 기능 장애를 치료 관련 운동 문제와 연결할 가능성을 보여준다.
• 조상 관련 한계: 본 연구는 PRS 성능의 조상 간 격차를 강조하며, 이를 유럽 GWAS 요약 통계량에 대한 의존성과 연관 불평형(LD) 패턴의 차이로 돌린다. 저자들은 전 세계적으로 다유전적 위험 예측 정확도를 높이기 위해 인구 특이적 GWAS 데이터가 시급히 필요함을 강조한다.

결론적으로, lyso-PRS는 생물학적 계층화 및 리소좀 표적 치료를 위한 고위험 환자 식별에 있어 유망성을 보이지만, 현재 그 효용성은 유럽 인구 집단에 국한되어 있으며, 생물학적 타당성을 검증하고 조상 간 전이성을 개선하기 위해 추가적인 대규모의 다양한 연구가 필요하다.