Single-cell multiomic mapping of genetic predisposition to childhood B-cell acute lymphoblastic leukemia

기술 요약: 소아 B세포 급성 림프구 성 백혈병의 유전적 소인에 대한 단일 세포 멀티오믹스 매핑

문제 제기
소아 B세포 급성 림프구 성 백혈병(B-ALL)은 가장 흔한 소아 악성 종양이며, 유전된 생식선 유전적 변이가 질병 감수성에 상당한 역할을 한다. 전전체 연관 분석(GWAS)을 통해 최소 16개의 재현 가능한 위험 로커스(loci)가 식별되었으나, 그 근저에 있는 메커니즘은 여전히 잘 알려져 있지 않다. 주요 과제는 GWAS 신호가 연관 불평형(linkage disequilibrium, LD)으로 인해 수많은 연결된 변이들을 포함하고 있어, 인과적 변이와 비인과적 이웃 변이를 구분하기 어렵다는 점이다. 또한, 계산적 미세 매핑(fine-mapping)을 통해 후보 단일 염기 다형성(SNP)을 우선순위화할 수 있지만, 이러한 예측은 확률적인 수준에 머문다. 결정적으로, 위험 변이가 조혈 발달 과정 중 효과를 발휘하는 정밀한 세포 맥락, 즉 희귀하고 일시적인 B세포 전구세포 내에서 이를 기능적으로 해결할 체계적인 접근법이 부족하다.

방법론
저자들은 통계 유전학을 고해 resolution 단일 세포 멀티오믹스와 결적으로 결합하여 위험 변이를 그 기능적 메커니즘에 매핑하는 체계적이고 통합적인 프레임워크를 사용하였다:

1. 통계적 미세 매핑(Statistical Fine-Mapping): 연구진은 대규모 인종 혼합 B-ALL 메타 GWAS 데이터에 세 가지 상보적인 알고리즘(FINEMAP, SuSiE, CARMA)을 적용하여 16개 위험 로커스에 대한 95% 신뢰 집합(credible sets)을 정의했다. 이 집합들의 합집합은 1,054개의 유력한 인과적 변이를 생성했으며, 이들 대부분은 비부호화(non-coding) 영역에 위치했다.
2. 단일 세포 멀티오믹 매핑(scMultiome): 저자들은 190,309개의 세포와 34개의 구별된 세포 유형을 사용하여 인간 조혈 작용의 고해상도 참조 지도를 구축했다(염색질 접근성(ATAC-seq) 및 RNA 발현(scRNA-seq) 데이터 통합). 이 지도는 SCAVENGE 알고리즘을 사용하여 특정 세포 상태에서의 위험 변이 농축도를 평가하는 데 사용되었으며, 이 알고리즘은 세포 간 네트워크를 통해 유전적 연관 신호를 전파하여 형질 관련성 점수(trait-relevance scores, TRS)를 계산한다.
3. 대립유전자 특이적 염색질 접근성(caQTL) 분석: 156명의 B-ALL 환자로부터 얻은 벌크 ATAC-seq 데이터셋을 활용하여, 저자들은 염색질 접근성 양적 형질 위치(caQTL)를 식별했다. 이들을 미세 매핑된 위험 변이와 중첩시켜, 조절 요소의 활성을 변화시키는 "caSNP"(염색질 접근성 SNP)를 탐지했다.
4. 전사 인자(TF) 모티프 파괴: 연구진은 caQTL 효과 크기와 예측된 TF 결합 친화도 변화(motifbreakR 사용)를 통합하여, 특정 위험 대립유전자가 어떤 TF 결의 결합 부위를 파괴하거나 생성하는지 식별했다. 이를 통해 어떤 TF가 특정 위험 변이의 조절 영향을 매개하는지 추론할 수 있었다.
5. 변이-대상 유전자 매핑: 조절 변이를 대상 유전자에 연결하기 위해, 저자들은 두 가지 직교하는 접근 방식을 결합했다: (i) scMultiome 지도 내에서의 염색질 접근성과 유전자 발현 간의 상관관계, (ii) B-ALL 블라스트(blasts)의 Hi-C 및 ATAC-seq 데이터에 적용된 ABC(Activity-By-Contact) 모델.
6. 기능적 검증: 연구는 ELK3 로커스와 rs2268508 변이에 집중되었다. 검증 방법은 다음과 같다:
   • 리포터 분석(Reporter Assays): REH 및 Nalm-6 세포에서 렌티바이러스 리포터 구조체를 사용하여 위험 대립유전자와 비위험 대립유전자의 인핸서 활성을 테스트했다.
   • CRISPR 활성화(CRISPRa): 일차 배양된 CD34+ 조혈 모세포(HSPC)에서 내인성 인핸서를 활성화하여 ELK3 상향 조절을 측정했다.
   • 과발현 및 결손(Overexpression and Knockout): ELK3 과발현 세포의 RNA-seq 및 HSPC 유래 B세포의 scRNA-seq를 수행하여, ELK3 결손이 전사적 결과 및 분화 결함에 미치는 영향을 평가했다.

주요 결과

• 위험의 세포 맥락: SCAVENGE 분석 결과, B-ALL 위험 변이는 초기 B세포 발달 단계, 특히 V(D)J pro-B 및 small pre-B 세포에서 유의하게 농축됨을 밝혔다. 이러한 농축은 다른 조혈 계통에서는 관찰되지 않았으며, 이는 유전된 감수성이 주로 활발한 V(D)J 재조합 시기의 조절 프로그램을 교란함으로써 작용함을 시사한다.
• TF 네트워크 교란: 연구는 모티프 활성 또는 발현이 위험 변이 농축과 상관관계를 보이는 74개의 전사 인자를 식별했다. 이 TF들(EBF1, PAX5, RUNX1, IRF4 포함)은 B세포 계통 분화에 필수적이다. caQTL 분석 결과, 미세 매핑된 조절 변이의 거의 절반이 이러한 핵심 TF(예: RUNX1, SPI-B, PU.1)의 결합에 강력한 변화를 일으키는 것으로 나타났다.
• 고신뢰도 감수성 유전자: CRE-유전자 상관관계와 ABC 모델 접촉을 통합함으로써, 저자들은 알려진 드라이버(예: ARID5B, IKZF1, GATA3)와 새로운 후보를 포함한 34개의 고신뢰도 B-ALL 감수성 유전자를 식별했다.
• ELK3 및 rs2268508의 메커니즘:
  • ELK3의 인트론 인핸서에 위치한 rs2268508의 위험 대립유전자(C)는 더 강한 PU.1(SPI1) 결합 모티프를 생성하고 B세포 전구세포에서 염색질 접근성을 증가시킨다.
  • 기능 분석을 통해 위험 대립유전자가 강력한 인핸서로 작용하여 ELK3 발현을 선택적으로 상향 조절함을 확인했다.
  • ELK3 과발현은 NF-κB 신호 전달 및 림프구 활성화 경로를 억제하는 동시에 핵심 백혈병 줄기세포 시그니처를 촉진하며, 이는 ERG 과발현에 의해 유도되는 전사적 변화를 재현한다.
  • 반대로, 일차 배양 HSPC에서 ELK3를 결손시키면 B세포 성숙이 저해되어, 미성숙 및 성숙 B세포 상태로 전환되는 세포의 비율이 감소하고 B세포 전구체 드라이버들이 하향 조절되었다.

의의 및 주장
본 논문은 유전적 소인의 근저에 있는 조절 교란을 해독하는 종합적인 "변이-기능(variant-to-function)" 지도를 구축했다고 주장한다. 통계적 미세 매핑을 단일 세포 멀티오믹스 및 기능적 검증과 통합함으로써, 저자들은 흔한 생식선 변이가 어떻게 B세포 전구체 조절 네트워크를 재구성하는지에 대한 기전적 통찰을 제공한다.

구체적으로, 본 연구는 ELK3를 B세포 발달 및 백혈병 발생의 이전에 알려지지 않은 조절자로 규명했다. 저자들은 rs2268508 위험 대립유전자가 ELK3 발현을 강화함으로써 B세포 성숙을 방해(NF-κB 억제를 통해)하는 동시에 백혈병 줄기세포 상태를 촉진하여 B-ALL 감수성을 높인다고 상정한다. 이러한 이중적 역할은 또 다른 ETS 패밀리 TF인 ERG의 거동과 유사하다.

저자들은 이 통합적 멀티오믹스 프레임워크가 GWAS 발견을 기전적 이해로 전환하는 확장 가능한 플랫폼을 제공하며, 향로적으로 유전적 위험이 악성 종양을 개시하는 데 수렴되는 특정 조절 회로와 발달 창(developmental windows)을 식별함으로써 맞춤형 예방 접근법을 마련할 수 있다고 결론짓는다.