배경

골격근은 인체에서 포도당 대사와 인슐린 반응을 주도하는 가장 큰 대사 기관이다. 제2형 당뇨병이나 이상지질혈증 같은 심장대사 질환 위험을 높이는 유전 변이의 상당수 역시 골격근의 대사 경로와 깊이 맞물려 있다. 전장 유전체 연관 분석(GWAS)은 지금까지 수천 개에 달하는 질환 위험 유전자좌를 발굴해 냈으나, 비부호화 영역에 위치한 비번역 변이가 질환으로 이어지는 분자적 경로는 여전히 미궁에 머물러 있었다.

종전의 전사체 분석이나 유전체 연구는 조직 파쇄액을 통째로 갈아 분석하는 대량 분석(bulk) 방식에 의존해 왔다. 근섬유를 비롯해 근위성세포, 섬유지방 전구세포, 내피세포, 면역세포 등이 복합적으로 뒤섞인 근육 조직에서 각 세포 유형이나 특정 생리학적 상태에 따른 미세한 유전자 조절 이상을 감지하기란 불가능에 가까웠다. 근섬유 자체가 다핵 세포라는 구조적 특수성 때문에 기존 단일세포 분리 기법을 온전히 적용하기 어렵다는 기술적 난관도 장벽이었다. 질환 연관 변이가 어떤 세포핵 상태에서 비정상적 조절을 유발하는지 밝혀낼 고해상도 분자 지도가 절실했던 이유다.

핵심 발견

연구진은 인간 골격근 조직에서 단리한 단일핵을 바탕으로 염색질 접근성을 측정하는 단일핵 염색질 접근성 분석(snATAC-seq)과 단일핵 전사체 분석(snRNA-seq)을 동시에 수행해 통합 다중오믹스 참조 지도를 구축했다. 근섬유의 유형별 아형(Type I, Type IIa, Type IIx)과 미세환경을 구성하는 기질 세포핵을 정밀하게 분리한 뒤, 염색질의 열림 상태와 전사체 활성을 1세포핵 수준에서 대조했다.

단일핵 수준의 발현 정량 형질 위치(eQTL)와 염색질 접근성 정량 형질 위치(caQTL)를 체계적으로 매핑한 결과, 특정 세포 상태와 생리학적 맥락에서만 작동하는 조절 변이군이 뚜렷하게 식별됐다. 조직 전체를 갈아 분석했을 때는 평균 신호에 묻혀 보이지 않던 수천 개의 세포 특이적 QTL이 수면 위로 드러난 셈이다. 대사 질환 GWAS 통계 신호와 이들 QTL 신호의 동시 국소화(colocalization) 분석을 결합하자, 질병 감수성을 좌우하는 인과적 유전자 조절 회로가 드러났다. 인슐린 신호 전달 및 포도당 수송체(GLUT4) 이동 경로에 관여하는 원위 조절 인자 부위에서 특정 위험 대립유전자가 염색질 접근성을 선택적으로 차단하고 유전자 발현을 왜곡하는 기전이 직접 확인됐다.

의미와 전망

이번 연구는 심장대사 질환을 유발하는 비부호화 유전 변이가 골격근 내부의 특정 세포 유형과 기능적 상태에서만 병리적 전사 장애를 일으킨다는 사실을 명확히 입증했다. 유전자 수준에 머물렀던 표적 발굴의 지평을 세포 상태 특이적 조절 부위로 좁힘으로써 신약 후보물질 선별의 실패율을 크게 낮출 토대가 마련됐다.

다만 정적 조직 생검 시료를 중심으로 지도가 작성된 만큼, 운동이나 급성 식이 변화 같은 동적 대사 자극 상황에서 염색질 환경이 어떻게 재구성되는지 추적하는 후속 연구가 뒤따라야 한다. 골격근 생검 검체의 인종적 다양성을 넓히고 단일핵 다중오믹스 파이프라인의 분석 비용을 낮추는 표준화 작업도 향후 해결해야 할 과제다.