| 查看: 1579 | 回复: 7 | |||
[交流]
手把手教你做全基因组DNA甲基化测序分析
|
|||
本期，我们讲讲全基因组DNA甲基化实验怎么做，从技术原理、建库测序流程、信息分析流程和研究套路等四方面详细介绍。 01全基因组甲基化测序技术原理 表观修饰不需要改变 DNA 序列便能实现对性状的改变，表观修饰的改变与基因功能乃至细胞状态、发育、衰老、疾病等存在重要的关联。在众多的表观遗传修饰中，最为重要且研究最为广泛的修饰之一是 DNA 甲基化，而全基因组甲基化测序（WGBS-seq）无疑是最有效的研究手段。 全基因组甲基化测序利用重亚硫酸盐能够将未甲基化的胞嘧啶（C）转化为胸腺嘧啶 （T）的特性，将基因组用重亚硫酸盐处理后测序，即可根据单个 C 位点上未转化为 C 未转化为 T 的 reads 数目与所有覆盖的 reads 数目的比例，计算得到甲基化率。该技术对于全面研究胚胎发育、衰老机制、疾病发生发展的表观遗传机制，以及筛选疾病相关的表观遗传学标记位点具有重要的应用价值。 02全基因组甲基化测序建库流程 （一）样品检测 对DNA样品的检测主要包括2种方法： 琼脂糖凝胶电泳分析DNA降解程度以及是否有污染，检测具有明显的主带，且条带清晰； Qubit 2.0对DNA浓度进行精确定量，DNA检测总量不低于1ug。 （二）文库构建 样本检测合格后，使用Bioruptor系统将1µg样品基因组DNA与未甲基化的lambda DNA混合，然后将其片段化，平均大小约为250bp。片段化后，纯化的随机片段化DNA随后用T4 DNA聚合酶，Klenow片段和T4多核苷酸激酶的混合物进行修复，钝化和磷酸化末端。随后使用Klenow片段（3'-5'exo-）对钝的DNA片段进行3'腺苷酸化，然后与连接5'-甲基胞嘧啶而不是使用T4 DNA连接酶的胞嘧啶连接的衔接子进行连接。完成每个步骤后，使用磁珠纯化DNA。之后，根据说明使用ZYMO EZ DNA甲基化金试剂盒将未甲基化的胞嘧啶转化为尿嘧啶。最后，用JumpStart Taq DNA聚合酶进行PCR扩增，再使用磁珠对PCR产物进行纯化获得最终文库。 （三）文库质检 文库构建完成后，先使用Qubit2.0进行初步定量，稀释文库至1ng/ul，随后使用Agilent 2100对文库的insert size进行检测，insert size符合预期后，使用qPCR方法对文库的有效浓度进行准确定量（文库有效浓度> 2nM），以保证文库质量。 （四）上机测序 文库检测合格后，把不同文库按照有效浓度及目标下机数据量的需求pooling后在illumina Nova平台测序，测序策略为PE150。 03全基因组甲基化测序信息分析流程 （一）原始下机数据质控 原始下机数据为FASTQ格式，是高通量测序的标准格式。FASTQ文件每四行为一个单位，包含一条测序序列（read）的信息。该单位第一行为read的ID，一般以@符号开头；第二行为测序的序列，也就是reads的序列；第三行一般是一个+号，或者与第一行的信息相同；第四行是碱基质量值，是对第二行序列的碱基的准确性的描述，一个碱基会对应一个碱基质量值，所以这一行和第二行的长度相同。以下为一条read信息的示例： 原始下机数据包含建库时引进的接头序列以及质量过低的碱基，这些因素会导致后续比对到基因组的reads较少，从而导致得到的信息较少，因此需要进行过滤。利用trim_galore软件对原始数据进行去除接头序列及低质量碱基等质控步骤。 （二）序列比对 经过质控的reads需要根据与参考基因组的序列相似度比对到参考基因组上。相比于常规基因组及转录组测序，WGBS测序方法产生的数据的特点决定其在比对时存在三大困难： DNA片段正链和负链经过重亚硫酸盐转化后将不再反向互补，再经过PCR，便会产生四条不同的序列，这将大大增加比对时的计算量。 经过重亚硫酸盐转化后，DNA序列大部分C碱基被转化成T碱基，因此序列含大量T而缺乏C；经过PCR后，产生的互补链则含有大量A而缺乏G。这样便导致序列的复杂度降低（即序列的组成特征更单一），从而增加比对的难度。 C和T的比对是不对称的。经过重亚硫酸盐转化后，序列中非甲基化的C碱基（占大部分）被转化为T，这将导致测序序列与参考基因组不匹配，T既可能应该比对到T上，有可能应该比对到C上；而C则只能比对到C上。这也增加了比对的难度。 利用BSMAP软件进行比对。BSMAP进行比对时，先以参考基因组上C碱基的位置作为指导，将reads中对应参考基因组C碱基位置的T标记为C，其他T保持不变，从而使reads可以直接比对到参考基因组。 （三）甲基化水平计算 甲基化水平可根据未转化为 T 的 C 与转化为 T 的 C 的 reads 的比例计算得到，即： Beta-value = C-reads / (C-reads + T-reads) * 100% 其中，Beta-value 即为该胞嘧啶的甲基化水平，C-reads 为覆盖该位点的支持甲基化的reads 数目（测得该位点为 C 的 reads），T-reads 为覆盖该位点的不支持甲基化的 reads 数目（测得该位点为 T 的 reads）。计算原理示意图如下： （四）差异甲基化区域（DMR）鉴定及统计 DMR检测使用权威期刊发表的metilene软件。该软件先将基因组进行预分段，以排除较长序列中不包含CG位点的片段。随后，利用二元分隔算法，递归缩小检测范围，以搜索得到组间累积平均甲基化差异最大的区域，作为可能的DMR；最后，结合双重统计学检验（MWU-test和2D KS-test），得到准确的DMR。检测原理如下图所示： 本分析检测DMR的标准如下： 区域平均甲基化差异不小于0.1； CpG位点数不少于5个； 区域长度不小于50 bp； 甲基化水平差异统计检验的校正P值小于0.05； 2D KS-test检验P值小于0.05。 04总结：全基因组甲基化研究思路 DNA甲基化组学研究的核心内容在于对DNA甲基化数据的挖掘。DNA甲基化一般遵循三个步骤进行数据挖掘。 首先，进行整体全基因组甲基化变化的分析，包括平均甲基化水平变化、甲基化水平分布变化、降维分析、聚类分析、相关性分析等。 其次，进行甲基化差异水平分析，筛选具体差异基因，包括DMC/DMR/DMG鉴定、DMC/DMR在基因组元件上的分布、DMC/DMR的TF结合分析、时序甲基化数据的分析策略、DMG的功能分析等。 最后，将甲基化组学&转录组学关联分析，包括Meta genes整体关联、DMG-DEG对应关联、网络关联等。 05全基因组甲基化研究案例 Whole-Genome Bisulfite Sequencing of Two Distinct Interconvertible DNA Methylomes of Mouse Embryonic Stem Cells. 两种状态的小鼠胚胎干细胞的甲基化组学研究 背景 小鼠胚胎干细胞一般生长在含有血清的基质中，被称作血清干细胞(serum ESCs)；加两种激酶抑制因子使胚胎干细胞在无血清的情况下更能保持多能性的基态，这种干细胞称为2i干细胞(2i ESCs)；这两种状态的胚胎干细胞可以互相转化。以前这方面的甲基化研究大多基于质谱，覆盖度和研究结果有限，尚缺乏2i胚胎干细胞的甲基化组学研究。 方法 利用全基因组重亚硫酸盐甲基化测序（WGBS），对这两种可互相转换的小鼠胚胎干细胞进行甲基化组学研究 结论 全面准确的检测了两种小鼠胚胎干细胞的DNA甲基化修饰并进行了系统的比较；同serum ESCs相比，雄性2iESCs全局低甲基化；在血清中，雌性ESCs跟雄性2i ESCs类似呈现全局低甲基化，而在2i ESCs状态下，甲基化水平会进一步降低。 参考文献：[1] Ashburner, M. and C. A. Ball, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet, 2000, 25 (1): 25-9. [2] Dirk Schübeler. Function and information content of DNA methylation. Nature, 2015, 517: 321–326. [3] Frank Jühling et al. metilene: Fast and sensitive calling of differentially methylated regions from bisulfite sequencing data. Genome Research, 2016, 26: 256-262. [4] Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research, 2000,28(1): 27-30. [5] Tadafumi Kato Kazuya Iwamoto. Comprehensive DNA methylation and hydroxymethylation analysis in the human brain and its implication in mental disorders. Neuropharmacology, 2014, 80: 133-139. [6] Xiaojing Yang et al. Gene Body Methylation Can Alter Gene Expression and Is a Therapeutic Target in Cancer. Cancer Cell 26, 577–590. [7] Yuanxin Xi et al. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics, 2009, 10:232. [8] Gao F, et al. De novo DNA methylation during monkey pre-implantation embryogenesis. Cell Res. 2017 Apr;27(4):526-539. pii: cr201725. 以上就是关于全基因组甲基化测序实验流程和分析思路的介绍，全面的DNA甲基化研究整体解决方案。 |
回复此楼» 猜你喜欢
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有7人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
三无人员面上一把中散金币
+2/1922
-
发点金币玩
+5/590
-
面上得中
+1/470
-
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/228
-
中国基，散金币！！！！！！
+5/159
-
北京化工大学何静教授团队诚聘青年英才和博士后
+1/95
-
东华大学材料科学与工程学院丁奇老师招收2026级硕士研究生
+1/89
-
欢迎报考中科院苏州纳米所的氢能、燃料电池方向的博士生！（长期有效）
+1/57
-
中国科学院苏州纳米技术与纳米仿生研究所 仿生自适应智能系统实验室招收联培生
+1/39
-
现在生物医药相关有什么兼职可以做吗？
+1/35
-
国家蛋白质科学中心（北京），招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/18
-
新审查指南堵死了"一案双申"的老路，实用新型和发明只能二选一了？
+1/13
-
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/7
-
项目助理招聘启事
+1/5
-
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/4
-
细胞治疗全流程关键检测技术与研究策略
+1/3
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/2
-
量子物理团队招收博士生与博士后
+1/1
-
EI会议征稿！第二届先进复合材料、聚合物和纳米技术国际学术会议(ACMPN 2027)
+1/1
-
收到青A意见了，这句评语给我看破防了，你是多少年远离正经科研圈了呀？？
+1/1
简单回复
jiantu2楼
2026-08-04 15:12
回复
chance℃(金币+1): 谢谢参与
tfang3楼
2026-08-04 15:12
回复
chance℃(金币+1): 谢谢参与
2026-08-04 15:48
回复
chance℃(金币+1): 谢谢参与
2026-08-04 16:46
回复
chance℃(金币+1): 谢谢参与
lymjjh6楼
2026-08-04 16:51
回复
chance℃(金币+1): 谢谢参与
庚辰辛丑7楼
2026-08-04 17:05
回复
ckxbear8楼
2026-08-04 17:45
回复
