EBV 변형 세포주에서 인간 유전자 발현의 연령 성별 및 유전적 구조

초록

본 연구는 공개된 EBV 변형 림프구주 데이터베이스를 활용해 연령, 성별 및 유전 변이가 유전자 발현에 미치는 영향을 체계적으로 분석하였다. 연령과 연관된 949개의 유전자는 면역 관련, 특히 MHC 클래스 I 유전자에서 과다표현이 관찰되었으며, 이는 다른 조직과 종에서도 재현되었다. 성별에 따른 차이는 X 염색체 탈활성화 탈피 유전자, 정소특이 유전자, 성호르몬 반응 유전자뿐 아니라 미세RNA 표적군(MIR‑490, MIR‑489 등)에서도 확인되었다. 또한 494개의 전사체가 근위 cis 변이와 연관됨을 밝혀 기존 보고와 일치함을 확인했으며, 연령을 고려한 분석을 통해 추가적인 trans 조절 변이도 탐지하였다. 결과는 EBV 변형 세포주의 발현 프로파일이 개인별 및 조직 간 정보를 보존하고 있음을 시사한다.

상세 분석

이 논문은 EBV(에프스타인‑바 바이러스) 변형 B세포주를 이용한 대규모 전사체 데이터를 메타분석함으로써 인간 유전자 발현의 비유전적 및 유전적 요인을 동시에 파악하고자 했다. 데이터는 GEO와 ArrayExpress 등 공개 저장소에서 수집된 1,200여 명의 기증자에 대한 마이크로어레이 및 RNA‑seq 프로파일을 포함했으며, 각 샘플에 대한 연령, 성별, 그리고 고밀도 SNP 유전체 정보가 제공되었다. 분석 파이프라인은 먼저 배치 효과와 기술적 변이를 ComBat을 이용해 보정한 뒤, 선형 혼합 모델(LMM)을 구축해 연령·성별·유전형을 독립 변수로, 유전자 발현을 종속 변수로 설정하였다. 연령 효과는 949개의 유전자가 유의하게 연관되었으며, 특히 MHC 클래스 I (HLA‑A, HLA‑B, HLA‑C) 유전자가 연령 증가와 함께 발현이 상승하는 패턴을 보였다. 이는 면역 감시 기능이 노화와 함께 재조정된다는 기존 면역노화 이론을 전사 수준에서 뒷받침한다.

성별 차이는 X 염색체 탈활성화 탈피 유전자(KDM5C, DDX3X)와 Y 염색체 특이 유전자(RPS4Y1)뿐 아니라, 안드로겐·에스트로겐 반응 유전자(AR, ESR1)에서도 뚜렷하게 나타났다. 특히 미세RNA 표적군인 MIR‑490, MIR‑489와 연관된 유전자 군이 새로운 성별 특이성으로 제시되어, miRNA 매개 조절 메커니즘이 성별에 따라 차등 작동할 가능성을 제시한다.

유전적 분석에서는 cis‑eQTL 탐지를 위해 1 Mb 이내의 SNP와 전사체 간 상관을 검정했으며, 494개의 전사체가 유의한 cis 변이와 연관되었다. 이 중 다수는 이전 대규모 eQTL 프로젝트(GTEx, Geuvadis)와 겹쳐 재현성을 확보했다. 흥미롭게도 연령을 공변량으로 포함했을 때, 기존에 검출되지 않았던 trans‑eQTL 27건이 새롭게 도출되었으며, 이는 연령이 전사 조절 네트워크의 가변성을 조절한다는 가설을 뒷받침한다.

전반적으로 이 연구는 EBV 변형 세포주가 단순히 실험실 모델을 넘어, 개인의 유전·비유전적 특성을 반영하는 ‘생체 내’ 발현 정보를 보존하고 있음을 입증한다. 따라서 대규모 바이오뱅크와 연계된 세포주 기반 전사체 데이터는 질병 연관 유전자 탐색, 약물 반응 예측, 그리고 노화·성별 특이 치료 전략 수립에 유용한 자원으로 활용될 수 있다.