Germline genomic and methylomic dynamics following three generations of early-life metabolic challenges

기술 요약: 3세대에 걸친 초기 생애 대사적 도전(metabolic challenges)에 따른 생식선 게놈 및 메틸로믹 역학

문제 제기
과영양과 같은 초기 생애의 대사적 도전과 같은 환경적 및 식이적 요인은 건강과 발달에 다세대적 영향을 미치는 것으로 알려져 있다. 이전 연구들은 대사적 교란이 후성유전적 기제(예: miRNA, DNA 메틸화)를 통해 부계 생식선을 통해 전달될 수 있음을 입증해 왔으나, 다세대적(연속적인) 노출에 따른 장기적인 게놈 영향에 대해서는 여전히 이해가 부족하다. 특히, 세 세대에 걸친 지속적인 초기 생애 대사 스트레스가 생식선 게놈(SNP 및 CNV)에 어떻게 영향을 미치는지, 그리고 이러한 게놈 변화가 에피게놈 변화(DNA 메틸화)에 의해 유도되거나 이와 상관관계가 있는지 여부는 불분명하다. 기존의 대부분의 연구는 연속적인 다세대 노출(환경적 비만 유도 물질의 더 현실적인 시나리오임에도 불구하고)보다는 단일 발달 노출(세대 간 전이)에 초점을 맞추고 있다.

연구 방법
본 연구는 초기 생애 대사적 도전을 모사하기 위해 리터 크기(litter size) 감소를 이용한 아웃브레드(outbred) ICR (CD1) 계통의 마우스 모델을 사용하였다.

• 실험 설계: F0–F2까지 총 3세대를 추적하였다. 과영양(Overnutrition, ON) 그룹은 어미당 새끼 수를 4마리로 줄인 리터를 가졌으며, 대조군(Control, CT)은 8마리로 설정하였다. 이 모델은 신생아 비만 및 그에 따른 대사 증후군을 유도하는 데 확립된 모델이다. 남성 개체들은 부계 혈통 효과를 분리하기 위해 각 세대에서 외부의 순수(naïve) 암컷과 교배되었다.
• 표현형 분석: 생후 4개월령에 생리적 및 대사적 파라미터(체중, 포도당, 인슐린, 간 및 지방 조직 무게)를 평가하였다.
• 시퀀싱 전략: 정량적 유전 변이 분석을 위한 GBS(Genotyping-by-Sequencing)와 메틸로믹 프로파일링을 위한 GBS-MeDIP(Methyl-Immunoprecipitation)를 병행하여 정자 DNA를 분석하였다.
• 게놈 분석:
  • SNPs: GATK 베스트 프랙티스를 사용하여 식별하였다. 유전적 유사성과 처리 효과를 평가하기 위해 주성분 분석(PCA) 및 계층적 클러스터링(Identity by State)을 사용하였다.
  • CNVs: GATK의 GermlineCNVCaller를 사용하여 생식선 복제수 변이(CNV)를 호출하였다. 위양성을 제거하기 위해 무작위 서브샘플링 대조군을 포함한 엄격한 필터링 과정을 적용하였다.
  • 주석 달기(Annotation): 전이 인자(TEs) 및 반복 요소(REs)에 대한 CNV 주석을 달기 위해 순열 검정(permutation tests)을 사용하여 농축도를 결정하였다.
• 에피게놈 분석:
  • DMRs: Mann-Whitney 검정을 사용하여 세대 내(CT vs. ON) 및 세대 간(F1 vs. F2) 비교를 통해 차별적 메틸화 영역(DMR)을 식별하였다.
  • 경로 분석: 메틸화된 윈도우와 관련된 유전자들에 대해 KEGG 및 PathBIX를 사용하여 경로 농축 분석을 수행하였다.
• 통합: 부계 메틸화 수준과 자손에서 나타나는 새로운 SNP 또는 CNV의 출현 사이의 관계를 조사하기 위해 Spearman 순위 상관관계를 계산하였다.

주요 결과

1. 게놈 영향 (SNPs): SNP 데이터의 PCA 결과, 서로 관련 없는 ON 가계들이 CT 가계가 친족 관계에 따라 엄격하게 클러스터링되는 것과 달리 함께 클러스터링됨을 확인하였다. 이는 대사적 도전이 공유된 게놈 시그니처를 유도했음을 시사한다. 처리는 유전적 분산의 약 4.6%를 설명하였다. 또한, ON 그룹은 CT 그룹에 비해 자손에서 새로운 SNP의 출현이 감소하였는데, 이는 유전적 다양성에 대한 제약을 나타낸다.
2. 게놈 불안정성 (CNVs): ON 그룹은 특정 CNV 이벤트(F1에서 F2로 유지된 7개의 추정 이벤트)의 출현을 보였다. 이러한 이벤트는 Long Interspersed Nuclear Elements (LINEs) 및 Long Terminal Repeats (LTRs)에 유의미하게 농축되어 있었다. 결실(deletions)은 세대를 거쳐 안정적으로 유지된 반면, 중복(duplications)은 확장되는 경향을 보였다.
3. 에피게놈 영향: 메틸로옴의 전역 PCA 결과, ON과 CT 그룹 사이에 명확한 분리가 나타나지 않았으며, 이는 광범위한 전역적 메틸화 변화가 없음을 의미한다. 그러나 경로 농축 분석 결과, ON 그룹은 CT 그룹에 비해 다당류 대사, 지방 분해, P53, mTOR 및 hedgehog 신호 전달과 관련된 경로에서 농축 신호가 약하거나 결여되어 있었다.
4. 국소적 메틸화 변화: DMR 분석 결과, F2 세대에서 특정 변화가 확인되었으며, 여기에는 LINEs 및 SINEs에서의 저메틸화와 Fhod3 및 Etl4 유전자의 고메틸화가 포함된다. 세대 간 비교를 통해 Fhod3 및 gamma-satellite 반복 서열에서의 DMR도 확인되었다.
5. 게놈-에피게놈 역학:
   • SNPs: CT 그룹에서는 높은 부계 메틸화가 자손의 높은 SNP 출현과 상관관계를 보였다(메틸화된 사이토신의 알려진 돌연변이 유발 가능성과 일치함). 반대로, ON 그룹에서는 역상관관계가 관찰되었는데, 이는 대사 스트레스 하에서 메틸화된 CpGs가 산화적 손상에 대한 보호 기제 덕분에 돌연변이에 덜 취약할 수 있음을 시사한다.
   • CNVs: CNV에서의 LINEs/LTRs 농축은 RE 메틸화의 조절 장애와 결합되어, RE 메틸화의 불균형이 CNV 형성을 유도할 수 있음을 시사한다.

의의 및 주장
저자들은 본 연구가 다세대 대사적 도전이 생식선에서 게놈 및 에피게놈 변화를 모두 유도할 수 있다는 증거를 제공한다고 주장한다. 주요 주장은 다음과 같다:

• 제약 vs. 불안정성: 대사적 도전은 SNP 수준에서는 유전적 다양성을 제약(새로운 돌연변이 감소)하는 동시에, 구조적 수준(CNV)에서는 게놈 불안정성을 유도하는 것으로 보인다.
• 기제: 관찰된 게놈 불안정성(CNVs)은 전역적인 DNA 메틸화 변화보다는 전이 인자(특히 LINEs 및 LTRs)의 조절 장애에 의해 매개될 가능성이 높다.
• 현실성: 3세대에 걸친 연속적인 노출을 추적함으로써, 본 연구는 단일 노출 세대 간 전이 모델에 비해 환경적 비만 유도 물질이 진화적 잠재력과 게놈 구조에 미치는 영향을 더 현실적으로 보여준다.
• 독창성: 본 연구는 특히 다세대 노출 시나리오에서 정자 DNA 메틸화 변화가 게놈에 어떻게 영향을 미치는지 다룬 최초의 통제된 연구이며, 에피게놈의 교란이 어떻게 구조적 게놈 변이로 이어질 수 있는지에 대한 복잡한 상호작용을 강조한다.

본 논문은 초기 생애의 대사적 도전이 전이 인자 활성에 의해 매개될 수 있는 장기적인 게놈 구조 변화를 초래하며, 게놈과 에피게놈 역학을 현실적인 노출 시나리오 하에서 함께 연구하는 것이 필수적임을 결론짓는다.