Causality of eosinophil cells and diabetic nephropathy: combined Single-cell RNA sequencing and Mendelian randomization analysis

기술 요약: 호산구 세포와 당뇨병성 신병증의 인과관계

문제 제기
당뇨병성 신병증(DN)은 지속적인 알부민뇨와 신장 기능 저하를 특징으로 하는 심각한 당뇨병성 미세혈관 합병증으로, 높은 이환율과 사망률을 초도한다. 염증 과정이 DN의 주요 동인이라는 점은 이미 알려져 있으나, 혈액 백혈구의 일종인 호산구와 DN 사이의 구체적인 인과관계는 아직 확인되지 않았다. 기존의 관찰 연구들은 호산구 수치와 DN 중증도 사이의 상관관계를 시사했으나, 인과관계를 명확히 규명하거나 단일 세포 수준에서의 근본적인 분자 기전을 밝혀내지는 못했다. 본 연구는 호산구와 DN 사이의 인과적 연결 고리를 조사하고, 관련된 특정 유전자 및 조절 기전을 식별함으로써 이러한 공백을 메우고자 한다.

연구 방법론
본 연구는 멘델 무작위 분석(Mendelian Randomization, MR)과 단일 세포 RNA 시퀀싱(scRNA-seq) 분석을 결합한 이중 접근 전략을 채택하였다.

1. 멘델 무작위 분석 (MR):

   • 데이터 소스: 호산구 세포 수(ieu-b-33, n=563,946)와 DN(finn-b-DM_NEPHROPATHY, n=213,746)에 대한 전전체 연관성 연구(GWAS) 데이터를 IEU Open GWAS 및 GWAS Catalog 데이터베이스에서 추출하였다.
   • 도구 변수 (IVs): 호산구 수와 연관된 440개의 SNP를 IV로 선택하였으며(P < 5×10⁻⁸), 연관 불평형(r² < 0.001)을 기준으로 필터링하고 F-통계량 > 10을 통해 검증하였다.
   • 분석: 역분산 가중법(IVW), MR Egger, 가중 중앙값(weighted median), 단순 중앙값(simple median), 최대 우도법(maximum likelihood)을 사용하여 단변량 MR(UVMR)을 수행하였다.
   • 민감도 분석: 결과의 견고성을 확보하기 위해 이질성(Cochran's Q), 수평적 다표현형성(MR Egger intercept), 그리고 leave-one-out(LOO) 분석을 실시하였다. Steiger 방향성 검정을 통해 인과 방향을 확인하였다.

2. 단일 세포 RNA 시퀀싱 (scRNA-seq) 통합:

   • 데이터셋: GEO 데이터베이스에서 벌크 RNA-seq 데이터(GSE320122)와 scRNA-seq 데이터(GSE131882)를 확보하였다.
   • 데이터 처리: Seurat 패키지를 사용하여 표준 품질 관리, 정규화 및 배치 효과 보정을 수행하였다. 세포 유형은 마커 유전자를 기반으로 주석을 달았다.
   • 유전자 식별: 벌크 데이터(GSE30122)에서 차등 발현 유전자(DEG)를 식별하였고, 단일 세포 데이터(GSE131882)에서 호산구 특이적 sc-DEG를 도출하였다.
   • 교집합 분석: 121개의 벌크 DEG와 153개의 호산구 sc-DEG를 교차 분석하여 5개의 핵심 유전자를 식별하였다.
   • 기능 분석: 유전자 온톨로지(GO), KEGG 경로 농축, 유전자 세트 변이 분석(GSVA), 의사시간 궤적 분석(Monocle2)을 수행하였다. 세포 간 통신은 CellChat을 사용하여 분석하였다.
   • 조절 네트워크 구축: MicroCosm, TargetScan, starBase, JASPAR 데이터베이스를 사용하여 경쟁적 내인성 RNA(ceRNA) 네트워크(lncRNA-miRNA-mRNA)와 전사 인자(TF) 상호작용을 예측하고, Cytoscape을 통해 시각화하였다.

주요 결과

• 인과적 추론: MR 분석 결과, 호산구 세포 수가 DN의 인과적 위험 요인임을 확인하였다 (IVW: P = 0.025, OR = 1.153, 95% CI: 1.018–1.306). 민감도 분석(이질성, 다표현형성, LOO)은 이 결과의 신뢰성을 뒷받침하였으며, Steiger 검정은 올바른 인과 방향을 확인하였다.
• 유전자 식별:
  • 벌크 DN 샘플에서 121개의 DEG가 식별되었고, 호산구에서 153개의 sc-DEG가 발견되었다.
  • 5가지 핵심 유전자: 교차 분석을 통해 IL7R, CD53, BPTF, FNBP4, ANKRD36B가 호산구와 DN을 연결하는 중심 유전자로 식별되었다.
  • 발현 패턴: IL7R과 CD53은 호산구에서 높게 발현되었고, BPTF는 상피 세포에서, FNBP4와 ANKRD36B는 사구체 메산지움 세포에서 높게 발현되었다.
• 기능적 농축:
  • 핵심 유전자들은 케모카인 신호 전달, 산화적 인산화, 세포 부착 및 엔도시토시스(endocytosis)를 포함하는 경로에 농축되어 있었다.
  • GSVA 분석 결과, DN 그룹의 호산구는 활성 산소 종(ROS), 혈관 신생 및 상피-간엽 이행(EMT)과 관련된 하향 조절된 경로들이 농축되어 있었다.
• 궤적 및 역동성: 의사시간(pseudotime) 분석 결과, 모든 발달 단계에서 DN 그룹의 호산구 수가 일관되게 높게 유지됨을 확인하였다. 핵심 유전자의 발현은 초기에는 감소하다가, 호산구 발달의 중기 및 후기 단계에서 DN 그룹에서 유의하게 증가하였다.
• 조절 네트워크:
  • miRNA: 핵심 유전자를 표적으로 하는 13개의 miRNA가 예측되었다 (예: FNBP4와 ANKRD36B를 표적으로 하는 hsa-miR-17-5p).
  • 전사 인자: 핵심 유전자를 조절하는 여러 TF가 예측되었다. 여기에는 IL7R 및 CD53을 조절하는 GATA2와 STAT3, FNBP4 및 BPTF를 조절하는 PPARG, 그리고 ANKRD3B, BPTF, CD53을 조절하는 FOXC1/HOXA5가 포함된다.

의의 및 주장
저자들은 본 연구가 관찰적 상관관계를 넘어 호산구 수와 당뇨병성 신병증 발생 사이의 인과적 관계를 입증한 첫 번째 증거를 제공한다고 주장한다. MR과 scRNA-seq을 통합함으로써, 본 연구는 IL7R, CD53, BPTF, FNBP4, ANKRD36B를 이 관계를 매개하는 잠재적 핵심 유전자로 식별하였다.

본 논문은 이 유전자들이 호산구의 증식, 활성화, 부착, 이동 및 산화적 인산화에 영향을 미침으로써 "유전자-호산구-신장 손상" 조절 축을 형성할 수 있다고 가정한다. 이 축은 신장 염증성 손상과 섬유화를 악화시켜 DN 진행을 유도하는 것으로 가설을 세웠다. 이러한 발견은 특히 질병의 면역병리학적 기전에서 호산구의 역할을 강조하며, DN의 진단 및 치료를 위한 새로운 후보 타겟을 제시한다.

한계점
저자들은 결론이 주로 생물정보학적 분석과 연관성 검증에 기반하고 있음을 인정하였다. 이들은 선정된 유전자들의 예측 능력을 검증하고 근본적인 분자 기전을 완전히 규명하기 위해서는 전향적 코호트 연구와 기능적 실험이 필요하다고 언급하였다. 또한, 유전적 집계(genetic aggregation)를 사용하는 통계적 접근 방식이 개인 간 차이에 대한 분석 범위를 제한한다는 점도 명시하였다.