dna 5hmc 测序(化学法)

  • 简介
  • 实验流程
  • 优势
  • 结果展示
  •         在人类基因组中,5mc占据总胞嘧啶的比例为2-8%,并且5mc修饰在很多生物学过程中,如基因表达、基因组稳定性维持、父母印记、x染色体失活、发育调控、衰老以及癌症的发生发展等生物学过程中发挥着不可替代的作用。5hmc作为5mc的氧化产物,主要存在于脑组织、神经细胞和胚胎干细胞中,在这些样本中的5hmc的含量在0.3-0.7%。5hmc是通过tet氧化5mc而得到,这种氧化作用依赖于亚铁离子和a-酮戊二酸。5hmc在胚胎干细胞的维持、正常血细胞生成以及恶性肿瘤和受精卵的发育等生物学过程中起着重要的作用,同时也在dna的去甲基化过程中发挥着承上启下的作用。

              为了更好更深入地了解5hmc的生物学功能,需要特异性强、方法成熟的dna 5hmc富集技术。然而,传统的基于重亚硫酸盐的方法无法区分5-hmc和5-mc。 康成生物丨数谱生物dna 5hmc测序(化学法)能特异性富集dna 5hmc片段,同时有很高的富集效率,完全满足对dna 5hmc特异性检测的需要。该方法是基于beta-gt能够特异性地将葡萄糖结合在5hmc这一原理而开展的。康成生物丨数谱生物dna 5hmc测序(化学法)及数据分析,可以获得全基因组羟甲基化分布图,帮助研究者从一个新的角度解析胚胎发育、神经细胞分化以及疾病发生的分子机制。


    aksomics(原康成生物)为您提供一站式dna 5hmc测序(化学法)技术服务,您只需要提供保存完好的组织或细胞标本,我司就可以为您完成从dna提取、5hmc dna富集、文库构建、高通量测序到数据分析的一整套服务,并提供完整的实验报告。

  • 图释. dna 5hmc测序(化学法)测序的流程示意图

  • 高度特异性: 本方法可以从dna修饰(c,5mc,5hmc)中特异性识别5hmc


    figure 1. 分别对加入了含有3种不同修饰的spike in(分别为c、5mc和5 hmc)的dna样品(10 ug)进行富集试验,富集完成后,分别针对spike-in的p1,p2,p3 不同位置进行pcr检测,结果显示该方法具有良好的5hmc富集特异性。


    无密度偏好性: 该方法富集5hmc,不依赖于dna的修饰水平

    稳定性高:实验操作稳定性高,减少实验操作带来的误差。

    灵活度高:能够直接对有基因组信息的任意物种的进行dna 5hmc测序

    精确度高:能够在实际结合位点50个碱基范围内精确定位。

    提供数据可视化与文章发表级别的图谱


  • 1. peak识别及注释

    我司使用macs2进行peak识别,并根据ensembl数据库和ehancer数据库(enhanceratlas)的注释信息,对peak进行详细注释。

    (1)peak分类

    根据peak顶点与转录本的相对位置关系对peak进行分类,具体分类方法如下。



    figure1. peak分类


    table1. 各样本peaks区域


    (2)peak统计

    根据以上分类,对不同区域进行统计绘制饼状图。

    figure2. peak统计图(tfbs: 转录因子结合位点)


    对peaks和fragments在tss(转录起始位点)及转录本周围的分布情况进行统计并绘图。

     

    (a)peaks密度在tss周围的分布   (b) peaks密度在转录本周围的分布

    (c)标准化后的reads在tss周围的分布   (d) 标准化后的reads在转录本周围的分布

    figure3. peak分布情况


    修饰密度(peaks/c)用来展示修饰区域在整个染色体上的分布情况;将修饰水平分为低(0%~30%,绿色),中(30%~70%,蓝色),高(70%~100%,红色)三类,在circos图中进行展示(bin=1mbp)。

    figure4. 修饰密度在不同染色体上的分布 

    figure5. 不同修饰水平peaks在染色体上的分布


    2. 差异peak分析

    康成生物丨数谱生物使用diffbind进行组间或样本间比较的差异peak分析,并根据ensembl数据库和enhancer数据库(enhanceratlas)进行差异peak注释。

    table2. 差异peak区域


    根据差异peak注释到的不同区域,对区域数量进行统计绘制条形图,并以10kb为单位对peak在染色体上的分布进行统计并绘图。


    figure6. 差异peak统计展示  (左):差异peak分布统计图; (右):差异peak染色体分布图



    通过各样本在差异peak富集程度的系统聚类分析,可以了解各组样本间的关系。火山图通过-log10p_value和log2fold_change两个条件,对比较数据进行筛选,可以直观地展示两组数据之间富集差异的倍数变化及显著性的关系。


    figure 7. 差异peak聚类图和火山图 (左):聚类图; (右): 火山图



    3.差异peak的go & kegg富集分析

    (1)差异peak的go富集分析

    geneontology (go)是一种基因功能分类条目,有三个子条目,分别描述每个基因的分子功能(mf: molecular function)、细胞组成(cc: cellular component)、参与的生物学过程(bp: biological process)。通常的go富集分析就是利用统计学算法来找出一组差异表达基因和哪些具体功能条目联系最大,每个go条目都对应一个统计值p-value来表示显著性,p-value越小表示该go条目和输入的差异表达基因联系越大,即该组差异表达基因大部分具有该go条目对应的描述功能。

    table3. go富集分析结果



    figure 8. go富集结果图 (左):各go条目中出现的基因数;(右)mf、cc、bp前十富集条目柱状图



    (2)差异peak的pathway富集分析

    根据kegg database (kyoto encyclopedia of genes and genomes)中的生物学通路分类条目,来找出一组差异表达基因和哪些具体的功能条目联系最大。每个pathway条目都对应一个统计值p-value来表示显著性,p-value越小表示该go条目和输入的差异表达基因联系越大,即该组差异表达基因p-value越小表示该pathway条目和输入的差异表达基因联系越大。

    table 4. kegg富集分析结果


    figure 9. kegg富集结果图



    4. motif分析

    在遗传学中,motif是指一段短的且普遍存在的核苷酸或氨基酸序列。被认为有特定的生物学意义,比如dna上的蛋白结合位点。当motif出现在基因外显子区域时,其对蛋白结构具有重要作用。我司使用meme-chip(v4.9.1)中的dreme方法(适合于查找短的motif)来检测motif。将得到的motif与转录因子结合位点motif数据库(jasparcore2018)内的motif进行比较,从而得到相关motif的更多信息。

    table 5. motif比较结果


    figure 10. motif结果图 (左):motif结果图;(右)左图中序列反向互补图



    5. 比对结果可视化

    我司将测序数据储存为bedgraph文件,可以在igv或ucsc基因组浏览器中查看。

    figure 11. 比对结果可视化