拟南芥不同生态型中MET1依赖性CG甲基化对全基因组转录活性的沟道化作用

Canalization of genome-wide transcriptional activity in Arabidopsis thaliana accessions by MET1-dependent CG methylation

作者信息Thanvi Srikant, Wei Yuan, Kenneth Wayne Berendzen, Adrián Contreras-Garrido, Hajk-Georg Drost, Rebecca Schwab, Detlef Weigel
PMID36539836
发布时间2022-12
DOI10.1186/s13059-022-02833-5

实验完整度

高

包含18个生态型的多组学分析(BS-seq、ATAC-seq、RNA-seq)、功能验证(表型分析、分离比统计)以及机制探讨(DMR与DEG关联分析)。

主要模型

18个拟南芥早花生态型的met1突变体 Col-0、Ler-1等野生型拟南芥生态型

重点核对

使用CRISPR-Cas9敲除MET1,产生移码突变(18个生态型) 获取纯合met1突变体,部分另有杂合子及第二代纯合子 同一植株莲座叶配对进行BS-seq、ATAC-seq和RNA-seq,每个基因型三个生物学重复 BS-seq数据分析:Bismark比对TAIR10,MethylScore识别DMR,CG-DMR共2388个 ATAC-seq数据分析:MACS2 callpeak,DiffBind识别差异可及区域(dACRs),获得9505个HV-dACRs

摘要

背景:尽管CG甲基化在基因表达和转座子沉默中具有保守作用,但野生拟南芥不同生态型之间的全基因组CG甲基化差异显著。结果:为了验证我们关于CG甲基化整体减少会降低拟南芥生态型中表观基因组、转录组和表型多样性的假设,我们在18个早花生态型中敲除了MET1(CG甲基化所必需)。所有生态型的纯合met1突变体都表现出共同的发育缺陷,如矮化、开花延迟,此外还有生态型特异性的莲座叶结构、角果形态和育性异常。对全基因组甲基化、染色质可及性和转录组的整合分析证实,MET1失活大幅降低了CG甲基化,并在数千个位点改变了染色质可及性。虽然转座子激活的影响在所有生态型中同样剧烈,但对非转座子基因的定量影响差异很大。在去除全基因组CG甲基化后,生态型的整体表达谱彼此之间的差异变得更大,尽管少数在野生型生态型间表达谱多样的基因在突变体中倾向于变得更相似。大多数差异表达基因并未在顺式位置表现出染色质可及性或CG甲基化的改变,这表明MET1缺失可能对基因表达产生深远的间接效应,且这些效应在不同生态型间差异很大。结论:对不同拟南芥生态型中MET1需求性的系统分析揭示了CG甲基化的双重作用:对于许多基因,CG甲基化似乎使表达水平沟道化,甲基化掩盖了调控差异;然而,对于一小部分基因,CG甲基化增加了超越遗传编码差异的表达多样性。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
MET1依赖性CG甲基化如何影响拟南芥不同生态型间基因表达多样性和表型多样性?
核心机制
MET1介导的CG甲基化对大多数基因的表达水平起沟道化作用,掩盖了潜在的调控差异;但对少数基因,CG甲基化可增加表达多样性。这种效应依赖于野生型中的初始甲基化状态和遗传背景。
主要证据
在18个生态型中敲除MET1,结合BS-seq、ATAC-seq和RNA-seq分析,发现CG甲基化缺失后,多数基因表达变化与局部甲基化或染色质可及性改变不直接相关,表明大量间接效应;同时发现了少数基因的差异表达与顺式CG-DMR相关。
研究意义
该研究揭示了CG甲基化在不同遗传背景下对基因表达的调控具有上下文依赖性,强调了研究参考生态型Col-0之外的生态型的重要性,为理解表观遗传变异在自然群体中的作用提供了资源。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

CRISPR-Cas9敲除MET1

在18个早花拟南芥生态型中创建MET1功能缺失突变体,以研究CG甲基化缺失的效应。

利用CRISPR-Cas9系统,设计靶向MET1外显子7的gRNA,通过农杆菌介导的花序浸蘸法转化18个生态型。筛选T1代,并在T2代分离出不含转基因的杂合突变体,进一步繁殖获得纯合突变体。

2

多组学数据收集

在同一植株的莲座叶中配对获取甲基化、染色质可及性和转录组数据,以关联表观遗传状态与基因表达。

对野生型和met1突变体植株的莲座叶分别进行BS-seq、ATAC-seq和RNA-seq文库制备和测序。每个基因型设置三个生物学重复。

3

DNA甲基化和染色质可及性分析

评估MET1缺失对全基因组CG甲基化和染色质可及性的影响,并识别差异甲基化区域和差异可及区域。

使用Bismark比对BS-seq数据,MethylScore鉴定DMR;使用MACS2对ATAC-seq数据进行peak calling,DiffBind鉴定差异可及区域(dACRs)。

4

转录组分析

鉴定met1突变体与野生型之间的差异表达基因(DEGs),并分析其与TE和非TE基因的关系。

使用RSEM和bowtie2比对RNA-seq reads,DESeq2分析差异表达,并区分TE-DEGs和非TE-DEGs。

5

表型分析

评估met1突变体在不同生态型中的发育表型,包括莲座叶形态、角果形态和育性。

观察并记录野生型和met1突变体的莲座叶、角果和育性表型,并结合分离比统计评估纯合突变体的传递率。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
DNeasy植物基因组DNA提取试剂盒Qiagen--
TruSeq Nano文库制备试剂盒Illumina--
Illumina HiSeq3000测序仪Illumina--
MoFlo XDP流式分选仪Beckman Coulter--
Tn5转座酶Homemade--
DNA清洁与浓缩试剂盒Zymo--
SPRI磁珠----
CytoFlex流式细胞仪Beckman Coulter--
飞利浦GreenPower TLED模块Philips--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
CRISPR-Cas9敲除MET1
gRNA序列、靶向外显子、转化方法(花序浸蘸)、筛选标记(mCherry)、热激处理条件
阅读提示:Methods 部分“CRISPR/Cas9 knockout of MET1 in 18 A. thaliana accessions”
多组学数据收集
植物生长条件(光照、温度、湿度)、取样时间(发芽后25天)、取样组织(莲座叶)、生物学重复数(3)
阅读提示:Methods 部分“Plant growth conditions and tissue collection for large-scale sequencing”
BS-seq分析
比对软件及参数、DMR calling参数(最小覆盖度、甲基化差异阈值)、参考基因组(TAIR10)
阅读提示:Methods 部分“Processing of Bisulfite-seq data and DMR calling”
ATAC-seq分析
核提取方法、FACS分选参数、Tn5转座酶量、PCR循环数、peak calling参数
阅读提示:Methods 部分“Nuclei isolation for ATAC-seq”、“Fluorescence-activated cell sorting (FACS) for ATAC-seq”、“Processing of ATAC-seq libraries and peak calling”
RNA-seq分析
RNA提取方法、文库制备协议、测序深度、比对和定量软件(RSEM)、差异表达分析模型
阅读提示:Methods 部分“RNA extraction and RNA-seq library prep”、“Mapping and identification of DEGs”