调查H2A。Z沉积在植物染色质中,我们生成了H2A的高分辨率全基因组图。Z在拟南芥中的应用体内我们用于亲和纯化果蝇染色质的生物素化系统21.我们标记了拟南芥H2A。Z带有由大肠杆菌生物素连接酶BirA(生物素连合酶识别肽,BLRP),并创建共表达BLRP-H2A的转基因植物。Z和BirA。细胞定位显示BLRP-H2A。Z具有弥散的核分布,但被排除在异色色心之外(补充图1)与内源性H2A的模式相同。Z轴17.用微球菌核酸酶消化后,大部分为单核小体(补充图1),我们从根组织中纯化生物素化染色质,并将相关DNA与对照DNA在代表整个拟南芥基因组的贴片微阵列上进行联合杂交22为了确保我们的结果不受潜在标记伪影的影响,我们用内源性H2A抗体重复了实验。Z轴17我们还绘制了根中的DNA甲基化图(我们之前已经发布了来自气生组织的数据集22).
链霉亲和素下拉和免疫沉淀生成的图谱几乎相同(和补充图2). 最显著的特征是与DNA甲基化强烈的定量反相关(Pearson’s第页= −0.81;补充表1-2). H2A的明显峰。基因5′端附近的Z也很明显(). 为了更好地可视化H2A。Z分布,我们将所有拟南芥注释序列(包括基因、假基因和转座子)排列在它们的5′端,并将它们从1号染色体的顶部堆叠到5号染色体的底部(和补充图2). 该路线的一个明显特征是高H2A的垂直条带。大致对应于转录开始后第一个核小体的Z。这种H2A模式。Z沉积与酵母和人类中的沉积一致10–15表明这是真核基因的一个普遍特征。还有五条明显的低H2A水平条纹。Z公司。这些对应于五个拟南芥着丝粒周围富含转座子、高度甲基化的异染色质。这种结合模式正好与DNA甲基化模式相反(和补充图2).
拟南芥H2A高分辨率图谱。Z和DNA甲基化a、,H2A。拟南芥2号染色体的Z(绿色)和DNA甲基化(蓝色)图谱。每个竖线代表原木2测试样本信号除以输入控制信号的信号比(log2(测试/输入))。黑色圆圈表示着丝粒序列间隙的位置。b–c,常染色(位置547000–587000,b条)和一个异色(4407000–4463000,c(c))基因组区域。空中组织和根的DNA甲基化以蓝色显示;从两个独立的BLRP-H2A获得的HA2.Z剖面。Z转基因株系和内源性H2A的免疫沉淀。Z以绿色显示。线的上方和下方分别显示了顶部和底部链上的基因和转座子。H2A的5′峰。基因中的Z由中的方框强调b条H2A水平相对较高的非甲基化转座子。中的方框强调Zc(c).
H2A。Z和DNA甲基化相互排斥a–b,所有TAIR 7注释序列(31762)均在5′端对齐,并从1号染色体顶部堆叠到5号染色体底部。BLRP-H2A。Z在中显示为热图一; 根DNA甲基化显示在b条注意H2A之间的高度反相关。Z和甲基化。c–d,非甲基化转座元件(列于补充表3). BLRP-H2A。Z在中显示为热图c(c); 根DNA甲基化显示在d、e、,所有TAIR 7注释序列均为k个-表示聚集(k个=3)基于BLRP-H2A。Z图案,并显示为热图。为了进行比较,相同序列的根DNA甲基化在(f).
甲基化在基因组中分布不均匀。转座子甲基化程度高且均匀,而一些基因甲基化时间短,而大多数基因是非甲基化的22–26这三组序列显示了H2A的相应三相分布。Z信号:低H2A。Z水平存在于转座子、甲基化基因的中间水平和非甲基化基因中(补充图3). 一种可能性是H2A含量低。转座子中的Z是由内在序列偏好引起的,而不是DNA甲基化。为了测试这一点,我们检测了少量(49)未甲基化的拟南芥转座子(补充表3). 很明显,所有这些转座子都具有高H2A。Z水平,表明H2A较低。Z合并本身不是转座子的特征(和). 非甲基化转座子也缺乏任何可识别的H2A。Z峰值,表明这些是内源性基因的独特特征。无监督k个-表示基于H2A的拟南芥注释序列聚类。Z模式产生了三组与非甲基化基因、体甲基化基因和转座子密切相关的基因(,补充图4和补充表4). 同样是H2A。Z和DNA甲基化水平显示出显著的反相关(). DNA甲基化和H2A。因此,Z是相互排斥的染色质特征,我们的分析表明,这种关系独立于序列上下文、转录或转录潜能(补充数据和补充图5-11).
到目前为止,我们的结果表明甲基化和H2A之间存在强烈的反相关性。Z沉积,但我们无法区分哪个是因果关系。为了解决这个问题,我们利用了一个DNA甲基转移酶突变为零的株系MET1,金属1-64,27.中的突变MET1型导致总体DNA甲基化显著降低,但也导致其他甲基转移酶介导的显著超甲基化26我们推断,如果DNA甲基化影响H2A。Z沉积,DNA甲基化的变化应通过H2A的变化反映出来。Z分布。值得注意的是,因为满足1导致DNA甲基化的损失和获得,我们应该看到H2A的获得和损失。为了验证我们的假设,我们绘制了H2A图。Z、 以及DNA甲基化和转录金属1-6植物。
DNA甲基化的变化确实引起了H2A的变化。Z分布(和补充图12-13). 为了可视化这些变化,我们减去了野生型(WT)H2A。来自的Z数据集满足1H2A。Z数据集,因此高值表示H2A增加。Z并入满足1(补充图12). 信息位点示例如所示. The鱼类和野生动物管理局该基因通常具有5′甲基化并且缺少H2A。Z峰,失去启动子甲基化,获得5′H2A。Z英寸满足1(). 反转录转座子在5g13205在WT中严重甲基化,但失去甲基化并获得H2A。Z英寸满足1(). 基因在1g22000编码F-box蛋白,在满足1导致其5′H2A损失。Z峰值().
H2A。Z合并变更表1-6突变植物a–c,野生型(WT)根DNA甲基化(深蓝色),金属1-6根DNA甲基化(紫色),WT H2A。Z(抗体,绿色),WT H2A。从中减去Z轮廓金属1-6H2A。Z剖面(两组独立的配对实验,浅蓝色),以及金属1-6/WT转录(红色)鱼类和野生动物管理局在里面一,类copia转置元素在5g13205失去甲基化并获得H2A。Z英寸金属1-6在里面b条和F-box基因在1g22000它被高度甲基化并失去H2A。Z英寸金属1-6在里面c(c).的5′区域鱼类和野生动物管理局WT中的甲基化通过a.天,所有TAIR 7注释序列均在5′端对齐,并从1号染色体顶部到5号染色体底部堆叠。WT H2A。从中减去Z图案金属1-6H2A。Z图案显示为热图。之后的数据相同k个-意味着聚类(k个=3)如所示e(电子)为了进行比较,序列的根DNA甲基化排列如下e(电子)在中显示为热图f.g、,WT甲基化水平(左)和金属1-6探针的甲基化水平(右)表示H2A显著降低。Z英寸金属1-6(补充图12). 三个独立的甲基化数据集的直方图是累积的。背景中的灰色直方图显示了所有探头的信号分布。小时,核密度图,用于跟踪数据集中显示的所有探针的频率分布d日(黑色痕迹),转座因子上调金属1-6(红色迹线)和转座因子在金属1-6(蓝色痕迹)。
全面了解H2A。Z动态输入金属1-6,我们对齐并排列了所有带注释的拟南芥序列,如。在这个剖面上,同样明显的中心周围条纹也很明显(和补充图13)–H2A。转座元件中的Z水平升高,它们失去了大部分甲基化,并在满足122,23数据的无偏分类产生了三个簇,分别大致包含非甲基化基因、甲基化基因和转座子(,补充图13和补充表4,序列分类如下22). H2A的变化。Z与获得H2A的DNA甲基化序列密切相关。Z英寸满足1在WT中甲基化(和补充图13). 相反,H2A降低的位点。Z掺入在WT中未甲基化,但在金属1-6(). 总的来说,DNA甲基化的变化反映在H2A的变化上。Z以一种强烈认为甲基化抑制H2A的方式表达。Z公司。
因为一些转座子和基因在满足1植物22,我们有机会测试H2A。Z掺入受到甲基化的负面影响或受到转录的正面影响。在基因内,WT和H2A中的DNA甲基化之间存在着强大的相关性。Z中的更改金属1-6(平均皮尔逊第页= 0.51,补充表2),但转录与H2A之间没有相关性。Z变化(皮尔逊平均值第页= 0.05).鱼类和野生动物管理局,在满足1,降低了H2A水平。Z在基因体中,在WT中没有甲基化(). 类似地,在少数野生型中没有甲基化的转座子中,有两个(见4g10690和在5g35205)然而,在满足1(补充图14). 两者的H2A也较少。Z英寸满足1而在野生型中,与其他转座子相反。
因为只有大约一半的转座元件在满足1,我们可以问这些元素是否优先获得H2A。Z、 如H2A所预期。Z掺入与转录活性相关。为了确保数据集的大小和甲基化不是问题,我们比较了12500个代表激活转座子的探针和12500个表示沉默转座子并具有相同甲基化特征的探针。我们发现这两个转座子类在H2A中的富集程度相同。Z轴(和补充图15). 因此,DNA甲基化的改变,而不是转录,导致H2A的重新分布。我们观察到的Z满足1.
我们的结果表明,DNA甲基化排除了H2A。一个有趣的问题是H2A是否存在。Z也可以排除甲基化。我们的一些数据表明,情况确实如此。H2A最引人注目的特征。Z掺入,即5′基因峰,与DNA甲基化无关(和补充图4、6),但甲基化被强烈排除在这个区域之外22,23同样,H2A越高。转录较少的基因体内的Z水平(补充数据和补充图6-7)这可能解释了令人费解的观察结果,即基因甲基化的几率随着转录的增加而增加(高达约70第个百分位数)22,23.
为了解决这个问题,我们绘制了植物中DNA甲基化的图谱项目1(Swr1的保守催化成分)破坏H2A的正常沉积。Z轴17。甲基化的总体模式图1-5植物仍与WT相似(补充表5),但DNA甲基化有适度但持续的增加(补充图16). 可视化甲基化变化第1张我们从中减去匹配WT对照(F2同胞)的甲基化模式第1张并将结果数据显示为热图(和补充图16). 该分析揭示了基因体的全基因组超甲基化。使用ChIPOTle算法28,我们确定了1201个高甲基化区域(对应1172个基因)用于进一步分析(阈值p<10−7,补充表6).
H2A。Z保护DNA甲基化a、,所有TAIR 7注释序列均在5′端对齐,并从1号染色体顶部到5号染色体底部堆叠。WT甲基化模式从第1张甲基化模式显示为热图。b、,五个位点的亚硫酸氢盐测序结果。我们对每个基因型的12个克隆进行了测序,除了地址:1g69850(10个克隆第1张)和收件人4g38190(11个克隆第1张).c、,用TaqI消化来自亚硫酸氢转化基因组DNA的PCR产物,TaqI识别TCGA,只有在C未转化(因此甲基化)时才会切割。L=100 bp梯形图,Unc=未切割PCR产物,TaqI=用TaqI消化的PCR产物。请注意,在第1张与WT相比。日期:,所有基因在5′端对齐,每个100-bp间隔的平均分数从距离基因2 kb(负数)到基因3 kb(正数)绘制。使用5点滑动窗口对数据进行平滑处理。虚线表示对齐点。e、,基因根据转录水平分为百分位。红线追踪了在第1张在每个百分位内(左Y轴)。黑线表示DNA甲基化富集(所有基因),绿线表示H2A。非甲基化基因中的Z富集(右Y轴)。使用10点滑动窗口对数据进行平滑处理。右侧Y轴的刻度设置为从零开始,以便比较甲基化和H2A。Z.公司。
在植物中,任何胞嘧啶都可能发生DNA甲基化5大多数甲基化是在对称的CG位点中发现的,就像在动物中一样,是由MET1介导的,但在其他序列环境中也有大量的甲基化被其他甲基转移酶催化(因此在满足1)25,26。为了确定第1张突变在不同情况下影响DNA甲基化,我们使用亚硫酸氢盐测序分析了ChIPOTle评分为高甲基化的五个位点中单个胞嘧啶的甲基化:地址:1g69850(硝酸盐转运蛋白),地址:3G22340(类COPIA反转录转座子),收件人4g03480(含蛋白质的锚蛋白重复序列),收件人4g38190(纤维素合成酶)和地址:5g37450(一种蛋白激酶)。所有五个样本都显示CG甲基化适度但持续增加(),确认微阵列分析。WT或WT中的任何一个位点几乎没有非CG甲基化第1张(数据未显示)。有趣的是,所有基因座在WT中都有一些甲基化,所以我们在第1张可能主要是由正常轻度甲基化位点甲基化增加引起的,而不是从头开始先前未甲基化的基因座的甲基化。
鉴于第1页突变时,我们问过甲基化位点是否代表整个基因组。正如所料,第1张高甲基化基因含有高水平的H2A。WT中的Z(即通常在非甲基化基因中发现的Z;). 它们通常也富含低转录基因,在30岁左右含量最高第个转录百分位数(). 这种模式与正常的甲基化基因非常不同,甲基化基因在70岁左右最为普遍第个百分位数(),也不同于非甲基化基因,后者富含低表达和高表达基因22.第1张然而,高甲基化基因确实与H2A的总体分布密切相关。Z轴(). 这些位点还包括49个富含H2A的转座子中的17个。Z和未甲基化的WT(补充表3、6),10倍的过度代表(p=10−4,费希尔精确测试)。因此,DNA甲基化的首选靶序列(基因体和转座子)在第1张与WT中这些序列中存在的低水平DNA甲基化一致(). 高水平的H2A。在这些位点上发现的Z明显保护了它们不发生完全的DNA甲基化,这可能解释了所观察到的基因转录和DNA甲基化之间的关系22.
几十年来,甲基化如何沉默基因一直是一个棘手的问题。一个流行的模型是,与甲基化DNA结合的蛋白质通过招募组蛋白去乙酰化酶而产生沉默6然而,对小鼠进行的仔细的基因破坏研究表明,这些蛋白质不太可能完全解释甲基化诱导的抑制29,30以前的工作提供了强有力的证据证明H2A。Z有助于提升推广人员的能力16–19因此,排除H2A。Z代表了DNA甲基化导致基因沉默的一种新机制。H2A。反过来,Z掺入可能保护基因启动子免受DNA甲基化的影响,从而促进基因活性并防止沉默。假设DNA甲基化和H2A。Z都是古老的染色质成分,它们的相互作用可能在调节真核基因表达中起着重要的普遍作用。