评估CpG甲基化状态的保持程度
量化小鼠原始生殖细胞和植入前胚胎全局去甲基化前后CpG甲基化状态的保持比例。
比较E6.5外胚层和精子中的BS-seq数据，计算高甲基化（>80%）和低甲基化（<20%）CpG在两者中的一致性，并据此估算保持比例。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
原始生殖细胞和植入前胚胎中转录因子对DNA再甲基化的保护可解释跨代表观遗传
Protection from DNA re-methylation by transcription factors in primordial germ cells and pre-implantation embryos can explain trans-generational epigenetic inheritance
实验完整度
高整合BS-seq、DNase-seq、ATAC-seq、RNA-seq等公共数据进行多阶段分析，并辅以ChIP-seq和TRAP亲和力预测，形成了从相关性到机制推断的多层级证据。
主要模型
小鼠原始生殖细胞 小鼠植入前胚胎 胚胎干细胞
重点核对
E14.5m PGC DNase-Hi位点的定义（FPKM > 20）及对照DNase-Lo位点（FPKM = 0） E14.5m PGC DNase-Hi位点在ESC中的ATAC-seq信号（亚核小体片段<115bp） CpG甲基化水平计算时仅纳入BS-seq覆盖度≥10的CpG TRAP亲和力评分使用TEPIC计算，并取区域内所有转录因子的最大亲和力 DNase-Hi和DNase-Lo位点均来自DNase-seq峰内fimo预测的转录因子结合序列
摘要
背景：越来越多的证据表明，某些表观表型可以通过哺乳动物的雄性和雌性生殖系跨代传递。由于在原始生殖细胞发育过程中和受精后，亲本染色体中存在的大多数已知表观遗传标记会全局丢失，这些观察结果一直难以解释。结果：通过整合已发表的在原始生殖细胞和植入前发育多个阶段采集的BS-seq、DNase-seq、ATAC-seq和RNA-seq数据，我们发现，尽管全局CpG甲基化在原始生殖细胞中降至10%，在囊胚内细胞团中降至20%，但全基因组大部分CpG的甲基化状态在全局去甲基化后得以恢复。我们估计这类甲基化状态得以保留的CpG比例为78%。进一步，我们发现，在生殖系和胚胎发育的全局再甲基化阶段，转录因子结合位点的CpG在所有观察到的发育阶段均保持低甲基化。相反，在全局再甲基化阶段未被转录因子结合的位点的CpG，在全局去甲基化之前具有高甲基化水平，在全局去甲基化过程中被去甲基化，然后被再甲基化。结论：结果表明，转录因子可以通过确保CpG甲基化状态得以维持，在生殖细胞和植入前发育过程中充当表观遗传信息的载体。这些发现为哺乳动物表观遗传信息跨代传递的机制描述提供了基础。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 转录因子是否能在原始生殖细胞和植入前胚胎发育的全局DNA去甲基化过程中充当表观遗传信息的载体，从而维持CpG甲基化状态并解释跨代表观遗传？
- 核心机制
- 转录因子在全局再甲基化阶段的结合保护其结合位点附近的CpG免受从头甲基化，从而维持低甲基化状态；而未被结合的位点则遵循全局甲基化动态，在去甲基化后被再甲基化。
- 主要证据
- 整合BS-seq、DNase-seq、ATAC-seq和RNA-seq数据，发现E14.5m PGC中DNase-Hi位点在PGC和植入前胚胎中保持低甲基化，而DNase-Lo位点则经历去甲基化和再甲基化；ESC中ATAC-seq信号与PGC DNase-Hi位点高度重叠；ChIP-seq显示Mll2、Ash2l、Tet1富集及H3K4me3修饰。
- 研究意义
- 为哺乳动物跨代表观遗传的机制提供了描述基础，并提出了可检验的预测，即转录因子结合可维持CpG甲基化状态，从而介导跨代继承。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
分析PGC发育中转录因子结合与甲基化的关联
检验在PGC再甲基化阶段，转录因子结合位点是否保持低甲基化。
使用E9.5-E16.5 PGC的DNase-seq数据，将E14.5m高DNase信号（DNase-Hi）与无信号（DNase-Lo）位点比较，并分析其甲基化水平变化。
分析ESC和植入前胚胎中的对应模式
检验PGC中观察到的转录因子结合与低甲基化模式是否在ESC和植入前胚胎中重现。
利用ESC中的ATAC-seq数据（亚核小体片段<115bp）评估E14.5m PGC DNase-Hi位点的可及性，并分析这些位点在植入前胚胎中的甲基化水平。
分析转录因子结合位点上的共调节因子
探索转录因子结合是否与抑制从头甲基化的共因子相关。
检查ESC中ChIP-seq数据，比较DNase-Hi和DNase-Lo位点上Mll2、Ash2l、Tet1、H3K4me3、H3K27me3和Dnmt3a1的富集情况。
识别可重编程的转录因子结合位点
寻找在PGC发育中甲基化状态发生改变的位点及其潜在作用转录因子。
定义E9.5时低DNase信号而E14.5m高DNase信号的位点（E9.5-trace, E14.5m-Hi），以及相反模式的位点，并比较其甲基化水平变化，进行基序富集分析。
评估胚胎与成体组织间的CpG重编程
量化从胚胎到成体组织发育过程中CpG甲基化状态的重编程程度。
分析E7.5胚胎与成体肠道及其他组织的BS-seq数据，比较DNase-Hi和DNase-Lo位点的甲基化变化，并进行基序富集。
分析高亲和力结合位点的保护效应
检验转录因子结合位点的亲和力是否影响其保护低甲基化的能力。
使用TRAP计算转录因子结合位点的亲和力，比较DNase-Hi和DNase-Lo位点的亲和力分布，并分析高亲和力位点的甲基化水平。
分析IAP元件的甲基化保持与转录因子亲和力
探究IAP转座子对去甲基化的抵抗力是否与其转录因子亲和力有关。
对IAP LTR进行基序富集和TRAP亲和力分析，比较可重编程转录因子与非可重编程转录因子的亲和力差异，并检查DNase可及的IAP的甲基化水平。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
分析DNA甲基化
分析基因表达
分析组蛋白修饰和转录因子结合
预测转录因子结合亲和力
路径富集分析
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 数据分析 | Trimmomatic | -- | 0.38 |
| bowtie | -- | -- | |
| MarkDuplicates (picard-tools) | picard-tools | 2.1.0 | |
| macs2 | -- | -- | |
| bedtools | -- | -- | |
| pyadapter_trim.py | -- | -- | |
| bowtie2 | -- | -- | |
| liftover | -- | -- | |
| bismark | -- | v0.19.0 | |
| deduplicate_bismark | -- | -- | |
| bismark_methylation_extractor | -- | -- | |
| tophat2 | -- | -- | |
| cuffdiff | -- | -- | |
| 序列比对和注释 | Dfam数据库 | -- | -- |
| 数据分析 | R函数pheatmap | -- | -- |
| R函数stringdist() | -- | -- | |
| TEPIC | -- | -- | |
| Enrichr | -- | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据来源与处理 | 所有分析基于公开数据集，需核对数据GEO/ENA编号及版本（mm9/mm10） 阅读提示：Methods中的数据可用性部分 |
| DNase-Hi/Lo位点定义 | DNase-Hi: E14.5m PGC中RPKM>20；DNase-Lo: RPKM=0；需使用fimo扫描的转录因子结合序列 阅读提示：Results章节'TF binding at CpGs in E14.5m PGCs'及Methods |
| 甲基化水平计算 | 仅使用覆盖度≥10的CpG，加权平均按读段数加权 阅读提示：Methods中'BS-seq analysis'部分 |
| ATAC-seq分析 | 仅使用片段长度50-115bp的读段代表转录因子结合（TF ATAC） 阅读提示：Methods中'ATAC-seq analysis'部分 |
| 转录因子亲和力计算 | 使用TEPIC和TRAP计算，取区域内最大亲和力；需指定转录因子集合（如E16.5m PGC中表达） 阅读提示：Methods中'Estimation of TF binding affinities'部分 |
| IAP分析 | 使用mm9参考基因组中的IAP注释，计算最大亲和力时需区分可重编程和不可重编程转录因子 阅读提示：Results中'IAPs possess...'部分及Methods |