Multi-Omics Integration Reveals Three Potential Therapeutic Targets for Spina Bifida: NPC2, PCDH12 and IL10RB

기술 요약: 멀티오믹스 통합을 통한 척수 이분증의 세 가지 잠재적 치료 표적 규명

문제 제기
신경관 폐쇄 실패로 인해 발생하는 심각한 신경관 결손(NTD)인 척수 이분증(Spina bifida, SB)은 전 세계적으로 약 1,000명당 1명의 발생률을 보이는 중요한 글로벌 보건 부담 과제이다. 산전 진단(초음파, 모체 혈청 AFP) 및 외과적 중재의 발전에도 불구하고, 현재의 치료법으로는 신경학적 손상을 되돌리거나 근본적인 병태생리학적 문제를 완전히 해결할 수 없다. 그 원인은 엽산 대사, 평면 세포 극성(PCP) 신호 전달, 환경적 요인 등을 포함하여 매우 복합적이다. 신뢰할 수 있는 바이오마커와 새로운 치료 표적을 식별하는 데 있어 결정적인 공백이 존재하며, 표현형에서 단백질로 이어지는 관찰만으로는 인과 관계를 확립하거나 특정 질병 기전을 밝혀내는 데 한계가 있다.

연구 방법론
본 연구는 척수 이분증의 인과적 바이오마커와 약물 투여 가능 표적을 식별하기 위해 단백질 전 범위 멘델 무작위 분석(MR)과 전사체 데이터 분석을 결합한 멀티오믹스 통합 전략을 채택하였다.

1. 데이터 소스:

   • 단백질체학(Proteomics): 세 가지 대규모 단백질 양적 형질 위치(pQTL) 데이터셋을 활용하였다: Decode (SomaScan, 35,559명), UKB_PPP (Olink, 54,219명), 그리고 Prot_a (SomaScan, 3,301명).
   • 유전학: SB 유전적 연관성 데이터는 210명의 사례와 408,394명의 대조군으로 구성된 Pan-UK Biobank (PAN_UKB)에서 가져왔다.
   • 전사체학(Transcriptomics): SB 사례 및 대조군의 태반, 태아, 뇌, 척수 조직 샘플을 포함하는 두 개의 Gene Expression Omnibus (GEO) 데이터셋(GSE252552 및 GSE33111)을 분석하였다.

2. 분석 파이프라인:

   • 멘델 무작위 분석 (MR): 도구 변수(SNP)는 게놈 전 범위 유의성($p < 5 \times 10^{-8}$) 및 연관 불평형 프루닝($R^2 < 0.1$)을 기준으로 선택되었다. 인과 효과는 역분산 가중치(IVW), Wald ratio, 가중 중앙값(weighted median), MR-Egger 방법을 사용하여 추정되었다. 민감도 분석에는 이질성을 위한 Cochran's Q 검정, 다표현성(pleiotropy)을 위한 MR-PRESSO, Steiger 방향성 검정이 포함되었다.
   • 전사체 분석: 차등 발현 유전자(DEG) 분석은 limma 패키지를 사용하여 수행되었다. 유전자는 $|log2FC| > 0.25$및$p < 0.05$에 의해 필터링되었다. 기능적 농축(GO/KEGG)은 clusterProfiler를 사용하여 수행되었다.
   • 유전자 우선순위 지정: 후보 유전자는 MR 유의적 단백질(Bonferroni 교정 $p < 0.05$)과 일치하는 발현 방향(하향 조절된 유전자, OR < 1)을 보이는 차등 발현 유전자(DEG)의 교집합을 통해 식별되었다.
   • 검증 및 기전: 일반화된 요약 데이터 기반 멘델 무작위 분석(GSMR)을 사용하여 다표현성을 조정하였다. 단일 유전자 유전자 세트 농축 분석(GSEA)을 통해 생물학적 경로를 탐색하였다. 표현형 전 범위 연관 연구(PheWAS)를 통해 잠재적 부작용을 평가하였다.
   • 신약 개발: Enrichr 플랫폼(DSigDB)을 사용하여 잠재적 치료 화합물을 스크리닝하였다. 후보 화합물과 표적 단백질 간의 결합 친화도 및 안정성을 평가하기 위해 분자 도킹(SYBYL-X의 Surflex-Dock) 및 분자 동역학(MD) 시뮬레이션(GROMACS 2022.2)을 수행하였다.

주요 결과

• 인과적 단백질 식별: 단백질 전 범위 MR을 통해 SB와 통계적으로 유의한 인과 관계를 가진 34개의 단백질을 식별하였다.
• 후보 유전자 우선순위 지정: MR 결과와 전사체 데이터를 교차함으로써 세 가지 하향 조절된 후보 유전자를 우선순위화하였다:
  • NPC2: OR = 0.753 (95% CI: 0.646–0.879); 콜레스테롤 대사에 관여함.
  • PCDH12: OR = 0.636 (95% CI: 0.507–0.798); 세포 접착에 관여함.
  • IL10RB: OR = 0.665 (95% CI: 0.549–0.805); 면역 조절에 관여함.
• 경로 농축: 단일 유전자 GSEA는 JAK–STAT, 사이토카인-사이토카인 수용체 상호작용(CCRI), 세포 접착 분자(CAMs) 경로에서 유의미한 농축을 나타냈다.
• GSMR 검증: GSMR 분석은 세 유전자 모두에 대해 유의미한 인과 관계($p < 0.05$)를 확인하였으며, 이는 SB 병태생리에 대한 이들의 역할을 강화한다.
• 약물 표적 예측:
  • Tamibarotene이 NPC2와 IL10RB 모두에 대해 유의미한 예측 상호작용을 보이는 최상위 후보 화합물로 식별되었다.
  • Decitabine 또한 두 유전자에 대한 중첩을 보였다.
  • **콜레스테롤(Cholesterol)**은 NPC2에 대한 내인성 양성 대조군 역할을 하였다.
• 구조적 검증:
  • 분자 도킹: NPC2는 콜레스테롤과 강한 결합(총 점수: 7.262)을 보여 계산 접근법을 검증하였다. Tamibarotene은 IL10RB와 유리한 결합(총 점수: 6.019)을 보였다.
  • MD 시뮬레이션: NPC2-콜레스테롤 복합체는 빠른 안정성에 도달하였다. IL10RB-tamibarotene 복합체는 50 ns 시뮬레이션의 대부분 동안 안정적이었으나, 후기 단계(34–50 ns)에서 구조적 변동을 보여 내인성 리간드에 비해 낮은 안정성을 나타냈다.

의의 및 주장
저자들은 본 연구가 엄격한 멀티오믹스 통합 전략을 통해 NPC2, PCDH12, IL10RB를 하향 조절된 새로운 보호 인자이자 척수 이분증의 잠재적 치료 표적으로 성공적으로 식별했다고 주장한다.

• 기전적 통찰: 본 연구 결과는 콜레스테롤 수송(NPC2), 신경관 폐쇄 중의 세포 접착(PCDH12), 그리고 JAK-STAT 및 CCRI 경로를 통한 면역 조절(IL10RB)을 포함한 SB의 특정 생물학적 기전을 시사한다.
• 방법론적 검증: 연구는 알려진 NPC2-콜레스테롤 상호작용을 양성 대조군으로 확인함으로써 계산 워크플로우를 검증하였다.
• 치료적 잠재력: 논문은 tamibarotene이 IL10RB(및 잠재적으로 NPC2)를 표적으로 하는 유망한 후보 화합물임을 강조하며, 약물 개발의 새로운 경로를 제시한다.
• 주의 깊은 전망: 저자들은 본 연구 결과가 주로 생물정보학적임을 명시하였다. 그들은 tamibarotene의 임상적 적용 가능성이 여전히 추측 단계이며, 약리학적 활성과 기전적 역할을 확인하기 위해 추가적인 in vitro 및 in vivo 검증이 필요함을 강조하였다. 또한, MR 및 GEO 데이터셋의 적은 사례 수(통계적 검정력 문제)와 비유럽인 인구 데이터의 부족이 결과의 일반화에 영향을 미칠 수 있음을 인정하였다.