基于病例对照组的检验差异甲基化基因功能区域的方法

摘要

DNA甲基化是一种重要的表观遗传学修饰,能够参与基因表达调控从而影响生理活动.在病例-对照组研究中,利用甲基化水平定位与疾病相关联的基因功能区域,有助于了解甲基化调控基因表达的机制,为后续研究提供疾病关联基因的线索.目前检验甲基化差异区域的方法存在一定不足:一是缺乏定义差异甲基化区域的统一标准;二是检测结果存在横跨多个基因功能区域的现象,难以进行生物学解释.本文基于得分检验(Zscore test)提出了以基因功能区域为单位的检验差异甲基化水平的方法cZST.该方法的特点是依靠位点物理距离和染色体全局信息修正协方差矩阵.模拟研究表明cZST在亚硫酸氢盐测序数据样本量小、读数低、扰动位点多的情形下具有高功效,且优于现有方法.将cZST运用到乳腺导管癌真实数据中,检测出与疾病显著相关的661个基因功能区和81个单位点.基因注释分析说明了cZST方法的有效性,并发现Dbl、Pleckstrin基因家族及Rho调控基因的甲基化与乳腺导管癌紧密相关,为后续研究提供指引.