나노그 결핍 마우스 배아줄기세포에서 동적 히스톤 아세틸화와 유전자 발현 매핑

초록

본 연구는 나노그를 억제한 마우스 배아줄기세포(ESC)에서 5일간 진행되는 분화 과정 동안 히스톤 H3K9,14 아세틸화와 전사량 변화를 전역적으로 측정하였다. 시간에 따라 아세틸화 감소와 mRNA 발현 변화가 유의하게 상관함을 확인했으며, 특히 ESC 특이 유전자군은 높은 초기 아세틸화 신호와 급격한 감소를 보였지만 발현 변화는 일관되지 않았다.

상세 분석

이 논문은 나노그(Nanog) 억제에 의해 유도된 마우스 ESC의 초기 5일간 분화 과정을 모델로, 히스톤 H3K9,14 아세틸화(H3K9,14ac)와 전사체 변화를 동시에 정량화하였다. 데이터는 ChIP‑seq 기반 아세틸화 프로파일과 마이크로어레이(또는 RNA‑seq) 기반 mRNA 발현 프로파일을 사용했으며, 각각 24시간 간격으로 수집된 5개의 시점에 대해 전사체 전반에 걸쳐 매핑되었다.

첫 번째 분석 단계에서는 아세틸화 신호를 유전자 프로모터(±2 kb)와 전사체 전반에 걸쳐 구간화하고, 각 구간별 평균 신호를 정규화하였다. 이후 시간별 차이값을 계산해 Δacetylation을 도출하고, 동일하게 Δexpression을 구해 두 변수 간 Pearson 및 Spearman 상관관계를 평가했다. 흥미롭게도 전체 유전체 수준에서 상관계수는 초기(0 h)에는 약 0.25 수준이었으나, 120 h에 이르러 0.45 이상으로 증가하였다. 이는 분화 진행에 따라 히스톤 아세틸화와 전사 활성 사이의 연관성이 강화된다는 가설을 뒷받침한다.

다음으로, 감독학습 모델(랜덤 포레스트, 서포트 벡터 머신, 로지스틱 회귀)을 이용해 Δacetylation을 입력 변수로 Δexpression(상승/하강/불변)을 예측하였다. 데이터 불균형을 해소하기 위해 SMOTE 기반 오버샘플링과 언더샘플링을 적용했으며, 10‑fold 교차검증을 수행했다. 가장 높은 정확도(≈78 %)와 AUC(0.81)를 보인 모델은 랜덤 포레스트였으며, 변수 중요도 분석 결과는 프로모터 근처 아세틸화 변화가 가장 큰 기여를 함을 보여준다.

특히 ESC 특이 유전자군(PluriNet, RNAi 스크린에서 도출된 핵심 조절인자 등)을 별도로 분석했을 때, 이들 유전자는 초기 단계에서 평균 아세틸화 신호가 전유전체 평균보다 1.8배 높았으며, 5일 동안 평균 40 % 이상의 아세틸화 감소를 보였다. 그러나 발현 변화는 일관되지 않아, 일부는 강하게 다운레귤레이션되었지만, 다른 일부는 변동이 미미하거나 오히려 상승하였다. 이는 히스톤 아세틸화가 전사 조절의 하나의 층일 뿐이며, 메틸화, 크로마틴 재구성, 전사인자 결합 등 다중 오믹스 요인이 복합적으로 작용함을 시사한다.

마지막으로, 저자들은 아세틸화와 전사 사이의 시간 지연 효과를 탐색하기 위해 교차상관 분석을 수행했으며, 최적 지연시간이 12–24 h임을 확인했다. 이는 히스톤 변형이 전사 활성화/억제에 직접적인 원인이라기보다, 선행적인 ‘프리시그널’ 역할을 할 가능성을 제시한다. 전체적으로 이 연구는 동적 에피제네틱 변화를 정량화하고, 이를 통해 유전자 발현을 예측하는 프레임워크를 제공함으로써, 줄기세포 운명 결정 메커니즘을 이해하는 데 중요한 발판을 마련한다.