DAP-seq与RNA-seq联合分析助力揭示AvERF73通过参与低氧应答和甲羟戊酸途径正调控猕猴桃的耐涝性

涝害(Waterlogging是长期阴雨或暴雨后,由于地表积水,地面径流不能及时排除,积水深度过大,时间过长,造成作物根部氧气不足,并产生乙醇等有毒有害物质,从而影响作物生长,甚至造成作物死亡,造成农业减产的灾害。

猕猴桃(Actinidia chinensis Planch)也称奇异果,原产于中国。猕猴桃因其独特的风味,富含维生素C,被誉为水果之王。目前已经跻身于世界主流消费水果之列。猕猴桃可以预防老年骨质疏松,抑制胆固醇的沉积,防治动脉硬化,改善心肌功能,防治心脏病,也能对抗癌起作用。

2024年2月2日,中国农业科学院郑州果树研究所果树生长发育与质量控制重点实验室的研究成果,发表在Horticultural Plant Journal期刊上(IF=5.7),文章题目为“AvERF73 positively regulates waterlogging tolerance in kiwifruit by participating in hypoxia response and mevalonate pathway” 。该研究使用DNA亲和测序(DAP-seq)技术和转录组测序(RNA-seq)技术揭示了AvERF73通过参与低氧应答和甲羟戊酸途径,正调控猕猴桃的耐涝性。

1.对AvERF73进行序列分析

作者在A. valvata KR5中鉴定出一个受洪涝诱导的ERF-VII基因,命名为AvERF73。随后,作者检测了AvERF73在不同组织中的表达情况,发现其在根系中的含量最高。系统发育树分析显示,AvERF73与AcERF74的蛋白序列最接近。蛋白序列分析显示,AvERF73含有一个AP2保守结构域,这个结构域含有Ala14和Asp19保守位点,说明AvERF73属于ERF转录家族之后作者将AvERF73与GFP融合表达,并在细胞核中检测到GFP荧光,证明AvERF73定位于细胞核中。

2.在拟南芥中过表达AvERF73提高了拟南芥的耐涝害能力

为研究AvERF73基因的功能,作者在拟南芥中过表达AvERF73基因,进行涝胁迫,发现过表达植株比WT植株更耐涝害,证实了AvERF73的过表达增强了拟南芥的耐涝害能力。

3.过表达AvERF73增强了猕猴桃的耐涝害能力

为进一步证实AvERF73在猕猴桃中的作用,作者利用“红阳”猕猴桃叶条外植体转化获得了两个过表达的猕猴桃品系OE53和OE65。涝胁迫后, 与过表达品系相比,WT植株出现了明显的叶片枯萎症状,叶绿素水平降低,根部H2O2和MDA含量升高,而过表达品系生长良好。表明AvERF73在应对涝胁迫方面发挥着有益的作用。

4.转录组分析鉴定AvERF73调控的基因

为了更好地探究AvERF73在涝胁迫下的分子机制,作者对OE53和WT植株进行转录组测序(RNA-seq)。qRT-PCR检测和相关性分析均表明转录组数据准确可靠。结果分析发现,在涝胁迫下,与WT植株相比,OE53植株的上调基因在“cellular response to hypoxia”途径显著富集,这些低氧应答基因(HRGs )包含很多转录因子,例如LOB,WRKY,锌指蛋白和NAC。,此外,上调基因在萜类化合物和甾醇生物合成过程也显著富集,并且编码MVA途径关键酶的基因表达上调。由此推测AvERF73可能调控HRGs的表达,同时激活MVA通路积累萜类化合物和甾醇来应对涝胁迫。

5.全基因组范围内鉴定AvERF73的靶基因

为了鉴定AvERF73调控的靶基因,作者以OE53和WT植物的根部为材料进行了DAP-seq两次重复共鉴定到3744个基因,同时鉴定到AvERF73结合的GCC-box (‘GCCGCC’) motif。随后作者将DAP-seqRNA-seq数据联合分析,找到了16个上调的AvERF73 可能直接靶向基因(DTGs)并进行重点关注。

6.AvERF73直接参与缺氧应答和MVA通路提高猕猴桃的耐涝渍能力 

为了进一步揭示AvERF73的潜在DTG。Venn图分析发现,在"cellular response to hypoxia"和MVA途径中,AcNAC022AcHMGS1基因为AvERF73的潜在DTGs。Dual-Luc实验发现AcNAC022AcHMGS1与AvERF73蛋白呈正相关关系。据此推测,AvERF73调控了AcNAC022AcHMGS1的表达。

小结

结果表明,AvERF73是猕猴桃耐涝性的重要调节因子。一方面,过表达AvERF73可以激活AcNAC022的表达,以应对涝胁迫引起的缺氧。另一方面,过表达AvERF73可能通过直接激活AcHMGS1表达来调节合成萜类化合物和甾醇的MVA途径。然而,AvERF73与AcNAC022AcHMGS1的相互作用调控猕猴桃耐涝能力的机制有待进一步探索。本研究结果为研究猕猴桃耐涝调控网络提供了候选基因和新思路。

  • 23
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
1. 使用MEME进行motif分析 MEME是一种常用的motif分析工具,可以用于从序列数据中识别出潜在的motif。利用dap-seq输出的peak序列,可以使用MEME进行motif分析。 首先,需要将dap-seq输出的peak序列转换为fasta格式的文件。可以使用bedtools将peak序列提取出来,并将其转换为fasta格式: ``` bedtools getfasta -fi genome.fa -bed peaks.bed -fo peaks.fa ``` 其中,genome.fa为参考基因组序列文件,peaks.bed为dap-seq输出的peak文件,peaks.fa为输出的fasta格式的peak序列文件。 然后,可以使用MEME对peaks.fa进行motif分析: ``` meme peaks.fa -oc meme_output -nmotifs 10 -minw 6 -maxw 20 ``` 其中,meme_output为输出结果的文件夹,-nmotifs指定需要识别的motif数量,-minw和-maxw分别指定motif的最小和最大长度。 2. 使用Homer进行motif分析 Homer是另一个常用的motif分析工具,也可以用于从序列数据中识别出潜在的motif。类似地,利用dap-seq输出的peak序列,可以使用Homer进行motif分析。 首先,需要将dap-seq输出的peak序列转换为bed格式的文件。可以使用bedtools将peak序列提取出来,并将其转换为bed格式: ``` bedtools sort -i peaks.bed > sorted_peaks.bed bedtools merge -i sorted_peaks.bed > merged_peaks.bed awk 'BEGIN{OFS="\t"}{print $1,$2,$3,"peak_"NR,".",$6}' merged_peaks.bed > peaks_homer.bed ``` 其中,peaks.bed为dap-seq输出的peak文件,sorted_peaks.bed和merged_peaks.bed为中间文件,peaks_homer.bed为转换后的bed格式的peak文件。 然后,可以使用Homer对peaks_homer.bed进行motif分析: ``` findMotifsGenome.pl peaks_homer.bed genome_dir homer_output -size 200 -p 8 ``` 其中,genome_dir为参考基因组序列文件夹,homer_output为输出结果的文件夹,-size指定motif的长度,-p指定使用的线程数。 3. 对比分析motif 对于使用不同的工具进行motif分析得到的结果,可以使用Tomtom进行对比分析。Tomtom是一个用于motif比对和聚类的工具,可以帮助用户在已知的motif数据库中搜索相似的motif,并将它们聚类为同一个motif家族。 首先,需要将使用不同工具得到的motif结果转换为meme格式的文件,并将其放入同一个文件夹中,如motif_dir。 然后,可以使用Tomtom进行对比分析: ``` tomtom -o tomtom_output -verbosity 1 -thresh 0.1 -eps -text -min-overlap 5 -dist pearson -no-ssc motif_dir/motif1.meme motif_dir/motif2.meme ``` 其中,tomtom_output为输出结果的文件夹,-thresh指定使用的阈值,-dist指定使用的距离度量方式,-no-ssc表示不使用自身比对,motif1.meme和motif2.meme为需要比对的motif文件。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值