转录因子结合区域的DNA甲基化模式表征其功能和进化背景

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

作者信息Martina Rimoldi, Ning Wang, Jilin Zhang, Diego Villar, Duncan T Odom, Jussi Taipale, Paul Flicek, Maša Roller
PMID38844976
发布时间2024-06-06
DOI10.1186/s13059-024-03218-6

实验完整度

高

包含WGBS和ChIP-seq实验数据的生成与分析,并整合系统发育比较、聚类和统计检验,形成多个独立证据层级。

主要模型

人类肝脏组织 猕猴肝脏组织 小鼠肝脏组织 大鼠肝脏组织 狗肝脏组织

重点核对

WGBS数据的测序覆盖度和深度 ChIP-seq数据的峰调用和保守性 物种间直系同源区域的比对 DNA甲基化模式的聚类参数选择

摘要

背景:DNA甲基化是一种重要的表观遗传修饰,在调节基因组功能中具有多种作用。其水平在基因组上空间相关,通常在抑制区域高,而在转录因子(TF)结合位点和活性调控区域低。然而,建立全基因组和TF结合位点甲基化模式的机制仍不清楚。结果:在这里,我们使用比较方法研究DNA甲基化与哺乳动物TF结合进化的关联。具体来说,我们通过实验描绘DNA甲基化,并将其与已发表的五种不同TF(CTCF,CEBPA,HNF4A,ONECUT1,FOXA1)在五种哺乳动物(人,猕猴,小鼠,大鼠,狗)肝脏中的占据谱相结合。TF结合位点甲基化水平低,但它们也经常具有中等甲基化水平。此外,结合位点受到其更宽结合区域中CpG甲基化状态的影响,即使核心结合基序中不存在CpG。采用分类和聚类方法,我们提取了TF结合区域DNA甲基化水平的独特且物种保守的模式。CEBPA,HNF4A,ONECUT1和FOXA1共享相同的甲基化模式,而CTCF的模式不同。这些模式表征了TF结合区域的不同功能和染色质景观。利用我们的系统发育框架,我们发现TF占据进化丢失后DNA甲基化增加,表明协调进化。此外,每种甲基化模式都有其自身的进化轨迹,反映其基因组背景。结论:我们的表观基因组分析表明DNA甲基化在跨物种TF结合变化中的作用,包括特定DNA甲基化谱表征TF结合并与其调控活性、染色质背景和进化轨迹相关。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
DNA甲基化在哺乳动物转录因子结合进化中的作用及其与结合位点功能的关系。
核心机制
DNA甲基化与TF结合协同进化:结合位点丢失后甲基化增加,甲基化模式与染色质背景和调控功能相关。
主要证据
五种哺乳动物肝脏WGBS和ChIP-seq数据,发现TF结合区域普遍低甲基化或中间甲基化,且不同甲基化模式与启动子、增强子等调控元件相关。
研究意义
表明DNA甲基化是TF结合跨物种变化的重要关联因素,为理解基因调控进化提供表观遗传视角。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

生成WGBS和收集ChIP-seq数据

获取五种哺乳动物肝脏的甲基化和TF占据数据

对五个物种肝脏样本进行WGBS测序,并从公共数据库获取五种TF的ChIP-seq数据。

2

定义TF结合区域和结合位点

确定TF结合区域及结合位点的基因组坐标

使用MACS2调用ChIP-seq峰,归一化峰长,并利用MEME进行基序发现,定义TF结合区域和结合位点。

3

分析TF结合区域的甲基化水平

描述TF结合区域的甲基化水平和分布模式

计算TFBR和TFBS的CpG数量和平均甲基化水平,分析其分布特征。

4

聚类甲基化谱

识别TF结合区域的不同甲基化模式

使用BPRMeth包对扩展的TFBR进行甲基化谱建模和聚类。

5

关联功能基因组注释

探究甲基化模式与调控元件和染色质状态的关系

将聚类结果与启动子、增强子、UMR/LMR/FMR等注释进行重叠分析。

6

进化保守性分析

评估TF结合和甲基化在物种间的保守性及共进化

利用EPO比对将TFBR映射到其他物种,定义保守结合、丢失、获得等类别,并比较甲基化水平。

7

统计关联分析

检验甲基化模式与进化保守性的关联

使用卡方检验分析甲基化模式与保守类别的独立性。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

转录因子结合检测
峰调用与基序发现
甲基化谱建模与聚类
功能注释

产品清单

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
WGBS实验
肝脏样本来源、重复数、测序平台、数据覆盖率
阅读提示:Methods: Whole-genome bisulfite sequencing (WGBS) protocol; 补充表1
ChIP-seq数据处理
数据获取编号、峰调用参数、重复性
阅读提示:Methods: ChIP-seq data processing
TFBR定义
峰长度归一化方法、CpG覆盖度阈值
阅读提示:Methods: Defining transcription factor binding regions (TFBRs) and their methylation coverage
甲基化谱聚类
扩展长度、聚类数、RBF数
阅读提示:Methods: 5mC and CpG profiles
进化保守性分析
物种比对版本、保守性定义
阅读提示:Methods: Evolutionary conservation of TF binding regions