Welcome to Chinese Journal of Tropical Crops,
Omics & Biotechnology

Identification and Expression Analysis of the DNA Methyltransferase and Demethylase Gene Families in Capsicum annuum L.

  • ZHANG Ying 1 ,
  • CAI Xiaotao 1 ,
  • XIE Bingchun 1 ,
  • WEI Lili 1 ,
  • XU Xiaowan 2 ,
  • ZHANG Bipei 1 ,
  • WU Zhiming , 1, *
Expand
  • 1. College of Horticulture and Landscape Architecture, Zhongkai University of Agriculture and Engineering, Guangzhou, Guangdong 510225, China
  • 2. Vegetable Research Institute, Guangdong Academy of Agricultural Sciences / Guangdong Key Laboratory for New Technology Research of Vegetables, Guangzhou, Guangdong 510640, China
* WU Zhiming,E-mail: .

Received date: 2021-07-20

  Request revised date: 2021-08-16

  Online published: 2022-02-22

Copyright

Copyright reserved © 2022.

Abstract

DNA methylation or demethylation is a highly conserved epigenetic modification involved in the regulation of numerous biological processes, including plant growth and development, secondary metabolism, and response to abiotic stresses, which controlled by DNA methyltransferase and demethylase genes. To fully understand the characteristics of DNA methylation modifying enzyme genes in the pepper genome, a total of 10 putative DNA methyltransferase and 4 demethylase genes were identified in pepper in the present study by using bioinformatics methods. The amino acids encoded by the genes were between 294-2037 AA. The genes were located on 10 different chromosomes, the number of exons contained in the 14 genes is between 1 and 21. The conserved motif compositions and exon-intron structures were systematically analyzed, and the results strongly supported the classification. Transcriptome data analysis proved that the expression levels of family member genes in different organs of pepper and different stages of fruit development were different. Among them, CaCMT1 and CaMET1-like showed the lowest expression in all organs, while CaDRM1-like2 and CaROS1-like2 showed the highest. The expression pattern of genes under high temperature and salt stress was analyzed by qPCR. Compared with the whole period of high temperature treatment, it was found that the expression of methylation modifying enzyme gene changed most obviously in materials treated for 3 hours, and the genes with the highest up-regulation were CaROS1-like2 and CaROS1-like1, respectively. The methylation modifying enzyme gene of the materials induced by salt stress was the most sensitive when treated for 12 hours, and the expression of 9 genes reached the peak, and the genes with the highest up-regulation were CaROS1-like2 and CaDRM1-like1. This study would provide a theoretical reference for further revealing the regulatory role of pepper DNA methylation modifying enzyme gene family in epigenetics.

Cite this article

ZHANG Ying , CAI Xiaotao , XIE Bingchun , WEI Lili , XU Xiaowan , ZHANG Bipei , WU Zhiming . Identification and Expression Analysis of the DNA Methyltransferase and Demethylase Gene Families in Capsicum annuum L.[J]. Chinese Journal of Tropical Crops, 2022 , 43(2) : 251 -261 . DOI: 10.3969/j.issn.1000-2561.2022.02.004

表观遗传(epigenetics)是指不依赖DNA序列的改变所产生的可遗传的变异[1],在动植物界均普遍存在。DNA甲基化(DNA methylation)与去甲基化(demethylation)是表观遗传修饰中最早发现且非常保守的分子机制之一,通过改变染色质结构、DNA构象、DNA稳定性或DNA与蛋白质的相互作用方式调控基因的表达[2] 。已有研究表明,植物DNA甲基化参与许多生物学过程的调控,在维持基因组稳定[3]、组织器官的生长发育[4,5,6]及响应胁迫应答[7,8]等过程中均发挥了重要作用[9]
植物DNA甲基化修饰过程通常由DNA甲基转移酶(DNA methyltransferase)和DNA去甲基化酶(DNA demethylase)共同作用,前者主要包括甲基转移酶MET1(Methyltransferase 1)、染色质甲基化酶CMT3(Chromomethylase 3)和结构域重排甲基转移酶DRMs(Domain-rearranged methyltransferases)等,分别负责CG、CHG和CHH(H=A,C或T)的甲基化修饰[10]。DNA去甲基化酶基因通常都含有保守的DNA糖苷酶结构域,属于DNA葡萄糖基酶(DNA glycosylases)家族,如DME(DEMETER)、DML2和DML3(DEMETER-like proteins 2 and 3)和沉默抑制子ROS1(Repressor of Silencing 1)等[11,12] 。基于基因组水平的生物信息学分析,目前研究者已分别从拟南芥[11]、水稻[13]、玉米[14]、花生[15]、丹参[16]、番茄[17]、茶[18]和石斛[19]等多种植物中鉴定出8~14个不等的DNA甲基化修饰酶基因,对这些基因功能的研究主要集中在模式植物如拟南芥[20,21,22]、水稻[23,24]和番茄[4, 25-26]中,基因突变后可引起拟南芥、水稻整个或部分基因组DNA甲基化水平发生改变。
辣椒(Capsicum spp.)是重要的茄科蔬菜作物,原产于墨西哥和中南美洲,约在明朝末年经丝绸之路和东南亚海路传入我国。我国辣椒年播种面积超过200万hm2,是种植面积最大的蔬菜种类之一[27]。本研究团队前期运用甲基化敏感扩增多态性技术(methylation sensitive amplified polymorphism, MSAP)分析了一年生辣椒(Capsicum annuum L.)和中华辣椒(Capsicum chinense Jacq.)种间与种内DNA甲基化多样性,发现辣椒表观遗传十分丰富[28]。同时分析了高温高湿胁迫下辣椒基因组DNA甲基化水平及状态的变化,推测DNA去甲基化可能是辣椒耐高温高湿的机制之一。最近也有研究表明,DNA甲基化修饰酶基因与植物激素相互作用,参与辣椒果实成熟的调控[6]。其他有关辣椒DNA甲基化表观遗传调控的研究甚少。
2014年,韩国[29]和中国[30]研究团队连续完成了对辣椒基因组的测序,为后续辣椒应用基础研究提供了重要支撑。本研究通过生物信息学手段从辣椒基因组中鉴定DNA甲基转移酶和DNA去甲基化酶家族成员,对其结构、染色体定位、系统进化、蛋白质保守序列进行分析,并比较该基因家族在辣椒不同组织及在逆境处理下的表达模式,为进一步挖掘辣椒甲基化修饰酶基因的功能提供参考依据。

1 材料与方法

1.1 材料

供试材料为一年生辣椒(Capsicum annuum L.)栽培种‘遵辣1号'(Zunla-1),由遵义市农业科学院辣椒研究所提供。种子经浸泡过夜,28℃催芽,出芽后播种于装有专用育苗基质的50孔穴盘,待植株长至5~6片真叶,一批露地种植于仲恺农业工程学院教学科研基地,按常规栽培管理;另一批用于逆境处理试验。逆境试验选择长势良好且长势一致的植株,分别经42℃高温和200 mmol/L NaCl胁迫处理,以未经处理正常生长的材料作为对照(CK),处理后0、3、6、12 h采取辣椒叶片,每个处理10株,设3个生物学重复,液氮速冻后立即放入-80℃超低温冰箱保存备用。

1.2 方法

1.2.1 辣椒DNA甲基化修饰酶基因家族成员的鉴定 在拟南芥基因组数据库(https://www. arabidopsis.org/)[31]搜索获得拟南芥中DNA甲基化修饰酶基因:AtMET1(AT5G49160.1)、AtCMT1(AT1G80740.1)、AtCMT2(AT4G19020.1)、AtCMT3(AT1G69770.1)、AtDNMT2(AT5G25480.1)、AtDRM1(AT5G15380.1)、AtDME(AT5G04560.2)、AtDML2(AT3G10010.1)、AtDML3(AT4G34060.1)和AtROS1(AT2G36490.1),下载基因相应的蛋白质序列,在茄科植物基因组数据库(https:// solgenomics.net/marker/SGN-M8338/details)[32]选择辣椒基因组蛋白质序列库(Zunla v2.0)进行BLASTP比对。利用ProtParam( http://web.expasy.org/protparam/)对甲基化修饰酶基因进行理化性质分析,预测等电点。
1.2.2 基因定位及结构分析 从辣椒基因组网站获取DNA甲基化基因所在的染色体位置及染色体大小,利用MapChart 2.2制作染色体定位图。同时从数据库下载基因的DNA和cDNA序列,利用在线工具GSDS 2.0( http://gsds.cbi.pku.edu.cn/)对甲基化修饰酶基因的结构进行展示[33]。利用MEME( https://meme-suite.org/meme/)在线分析软件分析蛋白质中的保守结构域。
1.2.3 构建系统进化树 为研究DNA甲基化修饰酶基因的进化关系,从NCBI分别下载马铃薯、番茄等物种中相应基因的蛋白质序列,采用MEGA-X软件内置的Clustal W软件对蛋白序列进行多重比对,默认软件本身参数设置。使用邻接法(neighbor-joining method,NJ)构建系统进化树,并对构建的系统进化树进行自检,bootstrap值设为1000[34]
1.2.4 组织特异性表达特征分析 利用辣椒基因组测序中获得的转录组数据集[35],分析辣椒甲基化修饰酶基因的表达特征,其中包含‘遵辣1号'辣椒不同组织(根、茎、叶、花和不同时期的果实等)所有注释基因的表达数据。转录本丰度用 FPKM(fragments per kilo bases per million reads)值表示。家族成员在不同组织中的FPKM值,取对数(log2FPKM)进行转换,用TBtools绘制表达量热图[36]
1.2.5 RNA提取与qPCR表达分析 采用试剂盒(华越洋)提取组织RNA。用1.5%普通琼脂糖凝胶电泳检测RNA的完整性,同时用核酸蛋白仪测定RNA的浓度。反转录合成cDNA并以此为模板,在CFX96 Connect实时定量PCR分析仪(Bio-Rad)上进行基因表达分析。扩增体系为25 μL:其中含100 ng/μL cDNA 2 μL,上下游引物(0.2 μmol/L)各0.5 μL,TB Green preminx Ex Taq(TaKaRa公司)12.5 μL,ddH2O 9.5 μL,扩增程序为:95℃预变性30 s;95℃变性5 s,60℃退火30 s,40个循环。用2-ΔΔCT法计算基因相对表达量,所得数据运用Origin软件作图。本研究所用qPCR引物见表1
表1 实时荧光定量PCR引物序列

Tab. 1 Primer sequence of qPCR

基因名称
Gene name
正向引物序列(5°-3°)
Forward primer sequence (5°-3°)
反向引物序列(5°-3°)
Reverse primer sequence (5°-3°)
CaMET1 TGCCAAAGAATCCCGAAGC TGCCATCCCAATCCAACCT
CaMET1-like TCAGAGCACCTGGAGTAAA TCTCAAGAAGGGAAGCAAC
CaCMT1 GGGTCTCACTCCCACAGA ATTTCATCCCTTCGCTCA
CaCMT2 AGGCGAGGACGATTACTT CCTGGGCACGATGATTAG
CaCMT3 CGAAATGCCTTACACTGAC CAACTGAAGAGGACGATGA
CaDRM1 CTTGGCAACTCATTCCAGGTT TCTGATTTGTTTGCTCCCACC
CaDRM1-like1 ACTTGGAACCCTGCATTA GCTTCTGGACCACATCTC
CaDRM1-like2 AAGTTGCTCTTCACCGTCTCG TGATTAGTTTGCTCCCACCAG
CaDRM1-like3 GGTGACAGACAGCGAGAT AGCAATGTTTCAGCGATA
CaDME CCTCAACAGACTGGTAAGGGAT GCCGCACGCATTCTACAC
CaDML3 AACTTGAGCACGACATTCTT ACATGACATCTAGCCACTTC
CaROS1-like1 TGGCACGGATGCGTCTTGTTC CCACGGCACTCACTTCGGGTT
CaROS1-like2 GGGGAAACTAACTATGTGC GATCTCCTCCTTAACCTCT
CaActin GGTGACGAGGCTCAATCCAA CTCTGGAGCCACACGAAGTT

2 结果与分析

2.1 辣椒DNA甲基化修饰酶基因的鉴定与分布

通过序列比对并去除重复,在‘遵辣1号'基因组中共鉴定获得14个DNA甲基化修饰酶基因(表2),包括10个DNA甲基转移酶基因和4个DNA去甲基化酶基因,分别是2个MET基因(CaMET1CaMET1-like),3个CMT基因(CaCMT1CaCMT2CaCMT3),4个DRM基因(CaDRM1CaDRM1-like1CaDRM1-like2CaDRM1-like3)和1个DNMT2基因(CaDNMT2),2个DME和DML基因(CaDMECaDML3)和2个ROS1基因(CaROS1-like1CaROS1- like2)。这些基因编码蛋白质的氨基酸长度范围为294~2037 aa,CDS序列长度为885~5835 bp,等电点在4.77~9.41之间。从基因在染色体上的分布看,除6号和11号染色体外,其余10条染色体均有分布(图1)。14个基因中有3个分布在12号染色体中,2个分布在1号染色体上,2个分布在10号染色体上,其余均单独分布在2号、3号、4号、5号、7号、8号和9号染色体上。
表2 辣椒基因组中的DNA甲基化修饰酶基因

Tab. 2 DNA methylation modifying enzyme genes in pepper genome

基因名称
Gene name
Zunla-1中对应的
基因名称
Gene name in Zunla-1
编码序列
Coding
sequence/bp
氨基酸序列特征
Amino acid sequence characteristics
拟南芥中同源基因名称
Name of homologous
genes in Arabidopsis
长度
Length/aa
等电点
pI
分子量
Molecular weight/bp
CaMET1 Capana04g000012 4659 1552 5.74 174 644.61 AT5G49160.1
CaMET1-like Capana12g001109 4662 1553 6.29 174 633.99 AT5G49160.1
CaCMT1 Capana01g001654 2262 753 4.88 84 873.54 AT1G80740.1
CaCMT2 Capana01g004297 3399 1132 6.29 174 633.99 AT4G19020.1
CaCMT3 Capana12g000016 2244 747 5.08 84 224.55 AT1G69770.1
CaDRM1 Capana10g002486 1851 616 4.98 69 452.50 AT5G15380.1
CaDRM1-like1 Capana07g000549 1944 647 4.81 72 631.62 AT5G15380.1
CaDRM1-like2 Capana02g000460 1329 442 5.29 50 091.82 AT5G15380.1
CaDRM1-like3 Capana05g002234 2037 678 5.38 76 366.68 AT5G15380.1
CaDNMT2 Capana08g000655 885 294 4.77 33 142.29 AT5G25480.1
CaDME Capana12g002335 5835 1944 9.41 184 423.43 AT5G04560.2
CaDML3 Capana03g000092 5666 2037 6.40 180 852.54 AT4G34060.1
CaROS1-like1 Capana09g002341 5496 1831 5.58 126 389.21 AT3G10010.1
CaROS1-like2 Capana10g001947 5451 1816 5.87 124 735.91 AT2G36490.1
图1 辣椒DNA甲基化修饰酶基因在染色体上的分布

Fig. 1 Distribution of DNA methylation modifying enzyme gene in pepper on chromosome

2.2 辣椒DNA甲基化修饰酶基因的结构分析

利用GSDS 2.0和MEME结构域在线分析工具,得到辣椒甲基化修饰酶基因的结构和蛋白质Motif分析的分布图(图2图3)。图2展示了基因的内含子与外显子,其中CaDML3含有21个外显子,而CaDRM1-like3仅有1个外显子。通过分析预测辣椒DNA甲基化修饰酶基因编码蛋白中的保守结构域可知,关系较近的基因拥有的保守域基本一致。如CaMET1CaMET1-like同时包含motif2、moti4、motif5和motif9;CaCMT1、CaCMT2CaCMT3共同拥有motif2、moti4和motif9;而motif7和motif8仅在DRMs中出现,motif1、motif3、motif6和motif10仅在DNA去甲基化酶基因(CaDML3、CaDMECaROS1-like)中出现(图3)。保守motif的分布预示着辣椒的这些甲基化修饰酶基因进化过程相对保守。
图2 辣椒DNA甲基化修饰酶基因的结构

Fig. 2 Structure of DNA methylation modification enzyme genes in pepper

图3 辣椒DNA甲基化修饰酶基因的蛋白质Motif分析

Fig. 3 Protein motif analysis of DNA methylation modifying enzyme gene in pepper

2.3 辣椒DNA甲基化修饰酶基因的系统进化分析

为探究辣椒和其他物种DNA甲基化修饰酶基因的同源进化关系,以拟南芥、辣椒、番茄和马铃薯中DNA甲基化修饰酶基因编码的蛋白质序列为材料,利用MEGA-X构建系统进化树(图4)。结果显示,辣椒基因组中的DNA甲基化修饰相关酶基因与拟南芥中该基因家族的分类结果一致。首先可将DNA甲基化修饰相关酶基因分为两大类,即I类DNA甲基转移酶基因和II类DNA去甲基化酶基因。DNA甲基转移酶基因家族可分为a、b、c和d四个亚家族,分别为DRMs、DNMT2、MET1和CMT3。DNMT2亚家族中的成员最少,仅包含1个基因(CaDNMT2),CMT3亚家族中包含3个辣椒CMT基因。DRMs亚家族中成员最丰富,包含4个辣椒DRM基因。在同科的马铃薯和番茄中,几乎都能找到辣椒DNA甲基化修饰相关酶基因对应的同源基因。
图4 辣椒、拟南芥、马铃薯和番茄DNA甲基化修饰酶基因系统进化树

:辣椒甲基化酶修饰基因; :拟南芥甲基化酶修饰基因; :马铃薯甲基化酶修饰基因; :番茄甲基化酶修饰基因。

Fig. 4 Phylogenetic tree of DNA methylation modifying enzyme genes in pepper, arabidopsis, potato and tomato

: Pepper methylase modified gene; : Arabidopsis methylase modified gene; : Potato methylase modified gene; : Tomato methylase modified gene.

2.4 辣椒DNA甲基化修饰酶基因在不同组织的时空表达模式分析

利用辣椒不同组织(包括根、茎、叶、花蕾、花、不同发育时期的果实)RNA-Seq数据库中辣椒DNA甲基化修饰酶基因对应转录本的FPKM值,绘制基因表达热图(图5)。结果显示,CaCMT1CaMET1-like在辣椒所有组织的表达量均很低,特别是CaCMT1,几乎在辣椒所有组织中均不表达;而CaDRM1-like2CaROS1-like2在所有组织中表达量均较高,特别是CaDRM1-like2,在根、茎和叶中的表达量均是最高的。在花蕾中,CaROS1-like1CaROS1-like2表达量最高,而在开放的花中,CaROS1-like2CaDRM1-like2CaDME表达量居前三。
图5 辣椒DNA甲基化修饰酶基因在不同组织及果实发育中的表达情况

Dev1~Dev5分别表示0~1、1~3、3~4、4~5 cm和成熟青果,Dev6表示转色期的果实,Dev7~Dev9分别表示转色后3、5、7 d的果实。红色越深,表示基因表达量越高;反之,蓝色越深,表示基因表达量越低。

Fig. 5 Expression abundance of DNA methyltransferase and demethylase genes in different tissues and developing fruits

Dev1-dev5 represent 0-1 cm, 1-3 cm, 3-4 cm, 4-5 cm and mature green fruits respectively, Dev6 represents fruits at color transformation stage, and Dev7-Dev9 represents fruits at 3 d, 5 d and 7 d after color transformation, respectively. The darker the red,the higher the gene expression; Conversely, the darker the blue, the lower the gene expression.

本研究详细分析了相关基因在辣椒果实发育与成熟过程中的表达情况。结果显示,除CaCMT1几乎不表达外,其他基因呈现不同的表达模式。在果实膨大过程中(F-Dev1~F-Dev5),CaROS1- like2CaDRM1-like2表达量一直较高,随着果实的发育,CaMET1CaMET1-likeCaROS1-like1的表达量逐步升高,在绿熟期前后达到最大值;与之相反,随着果实的发育,CaDRM1CaDMECaDML3的表达量逐步降低,在绿熟期前后达到最小值。在果实成熟(变色转红,F-Dev5~F- Dev9)过程中,CaMET1CaMET1-likeCaCMT2CaROS1-like1均表现出逐步下调的趋势,在完全成熟的果实中表达量达到最低值。

2.5 辣椒DNA甲基化修饰酶基因在高温胁迫下的时空表达模式分析

对高温胁迫处理下基因的表达模式进行qPCR验证(图6),结果显示,在高温胁迫处理后,基因表达呈现不同的变化趋势。其中,CaCMT2、CaCMT3CaMET1受高温胁迫表达量显著降低,分别在3 h和6 h达到显著差异。有9个基因呈现先上升后降低再上升的规律。处理后3 h,除CaCMT3、CaDRM1、CaDME、CaROS1- like1、CaMET1外,另外8个基因相比于其他时期表达量最突出。其中表达量变化最显著的是CaROS1-like2,与对照相比其表达量上调8倍,变化量最不明显的是CaCMT2,上调不到1倍。CaDRM1CaROS1-like1在处理12 h后,表达量出现峰值,分别上调11倍和25倍。DNA甲基化修饰酶基因的表达量随高温胁迫呈现不同的变化趋势,推测这些基因在辣椒响应高温胁迫过程中具有一定的调控作用。
图6 高温胁迫下辣椒叶片DNA甲基化修饰酶基因的表达情况

*表示差异显著(P<0.05),**表示差异极显著(P<0.01)。

Fig. 6 Expression analysis of DNA methylation modifying enzyme genes in leaves of pepper under high temperature treatment

* indicates significant difference (P<0.05), ** indicates extremely significant difference (P<0.01).

2.6 辣椒DNA甲基化修饰酶基因在盐胁迫下的时空表达模式分析

对盐胁迫下的13个基因进行qPCR验证,结果见图7。13个基因在盐胁迫后12 h内的表达量发生显著变化。CaCMT1、CaDRM1-like1、CaDML3CaROS1-like1受盐胁迫诱导表达,呈现逐步升高的趋势,在12 h表达量达到最高,其中表达量差异变化最大的是CaDRM1-like1,与对照相比上调9倍,变化最小的CaROS1-like1,与对照相比上调1.5倍。与之相反,CaCMT2、CaCMT3CaDRM1-like2CaDME在盐胁迫处理后,表达量显著降低。CaDRM1-like2、CaDRM1-like3、CaDMECaROS1-like2呈现先抑制表达后诱导表达的趋势。辣椒DNA甲基化修饰酶基因在盐胁迫下呈现不同的表达模式,推测这些基因参与辣椒盐胁迫响应诱导的机制可能不同。
图7 盐胁迫下辣椒叶片DNA甲基化修饰酶基因的表达情况

**表示差异极显著(P<0.01)。

Fig. 7 Expression analysis of DNA methylation modifying enzyme genes in leaves of pepper under NaCl stress

** indicates extremely significant difference (P<0.01).

3 讨论

DNA甲基化是由甲基化修饰酶基因催化S-腺苷甲硫氨酸(SAM)上的甲基与基因组DNA上的胞嘧啶第5位碳原子进行共价结合,将其修饰为5-甲基胞嘧啶(5-mC)的过程。DNA甲基化修饰酶(包括DNA甲基转移酶和DNA去甲基化酶)在调控甲基化与否和甲基化水平高低中起着非常重要的作用。DNA甲基化修饰酶家族的全基因组学鉴定已在拟南芥[10]、水稻[13]、玉米[14]和番茄[15]等多种植物中有过报道,部分基因如MET1[20, 37]ROS1[24, 38]在模式植物拟南芥和水稻中的功能也有一些深入研究。然而关于辣椒DNA甲基化修饰酶基因的研究还鲜见报道。因此,本研究基于已测序的辣椒基因组数据库[37,38,39]鉴定出DNA甲基转移酶和DNA去甲基化酶基因家族成员,并通过全面的生物信息学分析和表达分析,为辣椒甲基化修饰酶基因的功能挖掘提供参考依据。
辣椒基因组中共包含10个DNA甲基转移酶基因和4个DNA去甲基化酶基因,与同科的番茄、马铃薯中报道的基因数量一致[17, 25],比花生基因组中报道的少,这一差异的发生可能因为不同物种间基因组大小和复制方式存在差异。一般来说,在系统进化树上聚集在同一亚族的基因可能具有类似的功能[40,41]。本研究对拟南芥、番茄、马铃薯和辣椒的甲基化修饰酶基因进行系统进化关系分析,辣椒DNA甲基化修饰酶基因家族可分为a、b、c和d四个亚家族,分别为DRMs、DNMT2、MET1和CMT3,分类结果与前人在拟南芥[10]中的报道一致。说明辣椒DNA甲基化转移酶基因家族进化相对保守。
基因的组织特性表达特征一定程度上与基因的功能密切相关。已有研究结果表明,DNA甲基化修饰酶基因在植物不同组织和不同发育时期的果实表现出特异性表达。如柑橘CsMETl在幼苗和嫩叶中偏好表达[42],番茄中的SlMETI在幼嫩组织中表达量高[43]。番茄中4个DNA去甲基化酶基因(SlDML1~SlDML4),SlDML4在所有组织中的表达量极低,SlDML1SlDML3在叶片、花和幼果中表达量高,并均随器官发育成熟,表达量降低;随果实发育与成熟,从转色期开始SlDML2的表达量极显著升高,在橙色果实中表达量最高[25]。本研究发现,CaMET1在辣椒茎和根组织中的相对表达量最高,CaROS1-like1CaROS1-like2在辣椒花蕾组织中相对表达量高于其他基因,推测其在花蕾发育中发挥促进作用。辣椒果实发育过程中,随着果实的发育,CaMET1CaMET1-likeCaROS1-like1的表达量逐步升高,在绿熟期前后达到最大值,反方,随着果实的发育,CaDRM1CaDMECaDML3 的表达量逐步降低,在绿熟期前后达到最小值,推测他们参与了辣椒果实发育与成熟过程的调控。另外,对柑橘果实发育与成熟过程的甲基化组学分析表明,柑橘果实发育与成熟过程伴随DNA甲基化水平的升高,DNA去甲基化酶基因表达量降低,最终激活ABA信号基因的表达触发了柑橘果实的成熟[44]。今后也可对辣椒果实发育过程中的甲基化组进行分析,有利于更深入地解析DNA甲基化在辣椒果实发育中的调控作用。
当植物遭受非生物胁迫时,其体内的DNA甲基化会迅速发生动态变化,从而调控相关胁迫应答基因的表达,实现植物对非生物胁迫的快速响应[4, 45-46]。而DNA甲基化修饰酶基因也会响应胁迫应答,基因的表达量发生变化。如玉米中的ZmMET1aZmMET1bZmMET2aZmMET2b在干旱和盐胁迫下表达量显著下调,而ZmMET3aZmMET3bZmMET3cZmMET4的表达不受干旱和盐胁迫诱导[14]。在低温胁迫下,番茄SlMET1SlDRM5的表达明显受到抑制;而在盐胁迫下,SlCMT2的表达量显著升高,在4 h达到峰值,SlCMT3在12 h后表达量增加了13倍,而SlMET1SlDRM8的表达几乎不受盐胁迫影响[17]。茶树中6个DNA甲基转移酶基因,包括CsMET1CsDRM1CsDRM2aCsDRM2bCsCMT3CsDNMT2在冷胁迫下均显著下调,干旱胁迫下,CsMET1CsDRM1CsDRM2aCsDNMT2在12 h表达量显著降低,而CsROS1CsDML3CsMET1a等基因显著上调表达[18]。石斛基因组中DNA甲基转移酶和DNA去甲基化酶基因的表达也受冷胁迫和干旱胁迫诱导或抑制,呈现不同的表达特性[19]。本研究分析了辣椒DNA甲基化修饰酶基因在高温(42℃处理)和盐(200 mmol/L NaCl)胁迫下基因的表达情况,发现CaCMT2、CaCMT3CaMET1受高温胁迫表达量显著降低,同时CaDRM1CaROS1-like1CaROS1-like2等表达量显著上调,研究结果与在茶树中[18]的研究一致。在盐胁迫下,有的基因如CaCMT1、CaDRM1-like1、CaDML3CaROS1- like1的表达量逐步升高,有的基因如CaCMT2、CaCMT3CaDRM1-like2CaDME的表达量显著降低。也有一些基因如CaDRM1-like2、CaDRM1-like3、CaDMECaROS1-like2呈现先抑制表达后诱导表达的趋势。辣椒DNA甲基化修饰酶基因在盐胁迫下呈现不同的表达模式,推测这些基因参与辣椒盐胁迫响应诱导的机制可能不同。本研究结果为进一步揭示辣椒甲基化修饰酶基因的功能提供了参考依据。
[1]
ROBERTO B, TU S J, DANNY R. Molecular signals of epigenetic states[J]. Science, 2010, 330(6004): 612-616.

DOI

[2]
王静宇, 陈晓慧, 赖钟雄. 植物表观遗传修饰的分子机制及其生物学功能[J]. 热带作物学报, 2020, 41(10): 2099-2112.

WANG J Y, CHEN X H, LAI Z X. Molecular mechanism and biological function of plant epigenetic modification[J]. Journal of Tropical Crops, 2020, 41(10): 2099-2112. (in Chinese)

[3]
ZEMACH A, KIM M Y, SILVA P. Local DNA hypomethylation activates genes in rice endosperm[J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(43): 18729-18734.

[4]
MA Y, MIN L, WANG M J, WANG C Z, ZHAO Y L, LI Y Y, FANG Q D, WU Y L, XIE S, DING Y H, SU X J, HU Q, ZHANG Q H, LI X Y, ZHANG X L. Disrupted genome methylation in response to high temperature has distinct effects on microspore abortion and anther indehiscence[J]. Plant Cell, 2018, 30(7): 1387-1403.

DOI

[5]
ZHOU L L, TIAN S P, QIN G Z. RNA methylomes reveal the m6A-mediated regulation of DNA demethylase gene SlDML2 in tomato fruit ripening[J]. Genome Biology, 2019, 20(1): 156-163.

DOI

[6]
XIAO K, CHEN J, HE Q M, WANG Y X, SHEN H L, SUN L. DNA methylation is involved in the regulation of pepper fruit ripening and interacts with phytohormones[J]. Journal of Experimental Botany, 2020, 71(6): 1928-1942.

DOI

[7]
DEEPSHIKHA A, SANJAY K, MEENU K. Physcomitrella patens DNA methyltransferase 2 is required for recovery from salt and osmotic stress[J]. FEBS Journal, 2016, 283(3): 556-570.

DOI

[8]
FLÁVIA T, ADRIANA S H, PAULO C G. A role for epigenetic regulation in the adaptation and stress responses of non-model plants[J]. Frontiers in Plant Science, 2019, 10: 246-255.

DOI

[9]
ASSAF Z Y K, PEDRO S, JESSICA A R, BRADLEY D, MATTHEW D B, DANIEL Z, ROBERT L F. Local DNA hypomethylation activates genes in rice endosperm[J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(43): 18729-18734.

[10]
SIMON W C, IAN R H, STEVEN E J. Gardening the genome: DNA methylation in Arabidopsis thaliana[J]. Nature Reviews Genetics, 2005, 6(7): 351-360.

DOI

[11]
DELERIS A, HALTER T, NAVARRO L. DNA methylation and demethylation in plant immunity[J]. Annual Review of Phytopathology, 2016, 5(4): 579-603.

[12]
GALLUSCI P, HODGMAN C, TEYSSIER E. DNA methylation and chromatin regulation during fleshy fruit development and ripening[J]. Frontiers in Plant Science, 2016, 7(8): 807-820.

[13]
AHMAD F, HUANG X, LAN H X, HUMA T, BAO Y M, HUANG J, ZHANG H S. Comprehensive gene expression analysis of the DNA (cytosine-5) methyltransferase family in rice (Oryza sativa L.)[J]. Genetics and Molecular Research, 2014, 13(3): 5159-5172.

DOI

[14]
QIAN Y X, XI Y L, CHENG B J, ZHU S W. Genome-wide identification and expression profiling of DNA methyltransferase gene family in maize[J]. Plant Cell Reports, 2014, 33(10): 1661-1672.

DOI

[15]
WANG P F, GAO C, BIAN X T. Genome-wide identification and comparative analysis of cytosine-5 DNA methyltransferase and demethylase families in wild and cultivated peanut[J]. Frontiers in Plant Science, 2016, 7(1): 7-11.

[16]
LI J, LI C L, LU S F. Identification and characterization of the cytosine-5 DNA methyltransferase gene family in Salvia miltiorrhiza[J]. PeerJ, 2018, 6(6): 4461-4476.

[17]
GUO X H, XIE Q, LI B Y, SU H Z. Molecular characterization and transcription analysis of DNA methyltransferase genes in tomato (Solanum lycopersicum)[J]. Genetics and Molecular Biology, 2020, 43(1): 295-310.

[18]
ZHU C, ZHANG S, ZHOU C Z, CHEN L, FU H F, LI X Z, LIN Y L, LAI Z X, GUO YU Q. Genome-wide investigation and transcriptional analysis of cytosine-5 DNA methyltransferase and DNA demethylase gene families in tea plant (Camellia sinensis ) under abiotic stress and withering processing[J]. PeerJ, 2020, 8(8): 8432-8440.

[19]
YU Z M, ZHANG G H, TEIXEIRA S J, LI M Z, ZHAO C H, HE C M, SI C, ZHANG M Z, DUAN J. Genome-wide identification and analysis of DNA methyltransferase and demethylase gene families in Dendrobium officinale reveal their potential functions in polysaccharide accumulation[J]. BMC Plant Biology, 2021, 21(1): 2811-2838.

[20]
ELENA Z, MATT G, ANDREA K, PETER M. De novo DNA methylation activity of METHYLTRANSFERASE 1 (MET1) partially restores body methylation in Arabidopsis thaliana[J]. The Plant Journal, 2012, 71(6): 1029-1037.

DOI

[21]
PARK J, FROST J M, PARK K, OHR H, PARK G T, KIM S, EOM H, LEE I, BROOKS J S, FISCHER R L, CHOI Y. Control of DEMETER DNA demethylase gene transcription in male and female gamete companion cells in Arabidopsis thaliana[J]. Proceedings of the National Academy of Sciences of the United States of America, 2017, 114(8): 2078-2083.

[22]
JIANG J J, LIU J, SANDERS D, QIAN S M, REN W D, SONG J, LIU F Q, ZHONG X H. UVR8 interacts with de novo DNA methyltransferase and suppresses DNA methylation in Arabidopsis[J]. Nature Plants, 2021, 7(2): 184-197.

DOI

[23]
TAKAKI Y, YASUYO J, RIE T, IKUO N, SHIGERU I. The MET1b gene encoding a maintenance DNA methyltransferase is indispensable for normal development in rice[J]. Plant Molecular Biology, 2014, 85(3): 219-232.

DOI

[24]
LIU J X, WU X B, YAO X F, YU R, LARKIN P, LIU C M. Mutations in the DNA demethylase OsROS1 result in a thickened aleurone and improved nutritional value in rice grains[J]. Proceedings of the National Academy of Sciences of the United States of America, 2018, 115(44): 11327-11332.

[25]
LIU R, HOW-KIT A, STAMMITTI L, TEYSSIER E, ROLIN D, MORTAIN-BERTRAND A, HALLE S, LIU M C, KONG J H, WU C Q, CHARLOTTE D, CHAPMAN N H, MAUCOURT M, HODGMAN T C, TOST J, BOUZAYEN M, HONG Y G, SEYMOUR G, GIOVANNONI J, GALLUSCI P. A DEMETER-like DNA demethylase governs tomato fruit ripening[J]. Proceedings of the National Academy of Sciences of the United States of America, 2015, 112(34): 10804-10809.

[26]
YU Y, KAI T, TATSIANA U D, LIU W S, LV S H, LANG Z B, WANG X G, GAO J H, WANG W, NI W FE, CHU Z Q, ZHANG H, AVTAR K, ZHU J K, ZHANG H M. Critical function of DNA methyltransferase 1 in tomato development and regulation of the DNA methylome and transcriptome[J]. Journal of Integrative Plant Biology, 2019, 61(12): 1224-1242.

DOI

[27]
王立浩, 马艳青, 张宝玺. 我国辣椒品种市场需求与育种趋势[J]. 中国蔬菜, 2019(8): 1-4.

WANG L H, MA Y Q, ZHANG B X. Market demand and breeding trend of pepper varieties in China[J]. Chinese Vegetables, 2019(8): 1-4. (in Chinese)

[28]
李涛, 徐小万, 李颖, 李明珠, 王恒明, 徐晓美, 罗少波. 一年生辣椒(Capsicum annuum L.)与中华辣椒(Capsicum chinense Jacquin)DNA甲基化多样性分析[J]. 分子植物育种, 2014, 12(2): 306-315.

LI T, XU X W, LI Y, LI M Z, WANG H M, XU X M, LUO S B. DNA methylation diversity analysis of annual pepper (Capicum annuum L.) and Chinese pepper (Capicum chinense Jacquin)[J]. Molecular Plant Breeding, 2014, 12(2): 306-315. (in Chinese)

[29]
HU J H, CHEN X J, ZHANG H Y, DING Y. Genome-wide analysis of DNA methylation in photoperiod- and thermo- sensitive male sterile rice Peiai 64S[J]. BMC Genomics, 2015, 16(1): 102-134.

DOI

[30]
KUMAR G, RATTAN U K, SINGH A K. Chilling-mediated DNA methylation changes during dormancy and its release reveal the importance of epigenetic regulation during winter dormancy in apple (Malus × domestica Borkh.)[J]. PLoS One, 2016, 11(2): 934-966.

[31]
Reiser L, RHEE S Y. Using the Arabidopsis information resource (TAIR) to find information about Arabidopsis genes[J]. Current Protocols in Bioinformatics, 2005, 9(1): 1-51.

[32]
FERNANDEZ P N, MENDA N, EDWARDS J D, SAHA S, TECLE I Y, STRICKLER S R, BOMBARELY A, FISHER Y T, PUJAR A, FOERSTER H, YAN A, MUELLER L A. The sol genomics network (SGN)--from genotype to phenotype to breeding[J]. Nucleic Acids Research, 2015, 4(3): 36-41.

[33]
HU B, JIN J P, GUO A Y, ZHANG H, LUO J C, GAO G. GSDS 2.0: an upgraded gene feature visualization server[J]. Narnia, 2015, 31(8): 1296-1297.

[34]
KUMAR S, STECHER G, LI M, KNYAZ C, TAMURA K. MEGA X: molecular evolutionary genetics analysis across computing platforms[J]. Molecular Biology and Evolution, 2018, 35(6): 1547-1549.

DOI

[35]
QIN C, YU C H, SHEN Y, FANG X D, CHEN L. Whole- genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization[J]. Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(14): 5135-5140.

[36]
CHEN C J, CHEN H, ZHANG Y, THORMA R H, MARGARET H F, HE Y H, XIA R. TBtools: An integrative toolkit developed for interactive analyses of big biological data[J]. Molecular Plant, 2020, 13(8): 1194-1202.

DOI

[37]
SCHMIDT A, WÖHRMANN H, RAISSIG M. The polycomb group protein MEDEA and the DNA methyltransferase MET1 interact to repress autonomous endosperm development in Arabidopsis[J]. The Plant Journal, 2013, 7(3): 776-787.

[38]
JUNE S K, SATOSHI K, KAZUO S, KAZUKO Y. DNA demethylase ROS1 prevents inheritable DREB1A/CBF3 repression by transgene-induced promoter methylation in the Arabidopsis ice1-1 mutant[J]. Plant Molecular Biology, 2020, 9(30): 575-582.

[39]
刘卓毅, 于文菲, 蔡文丽, 刘子珊, 张雨, 袁媛, 伍炳华, 吕美玲. 辣椒CTL基因家族的全基因组鉴定和表达特征分析[J]. 热带作物学报, 2021, 5(25): 1-14.

LIU Z Y, YU W F, CAI W L, LIU Z S, ZHANG Y, YUAN Y, WU B H, LV M L. Genome wide identification and expression characterization of CTL gene family in pepper[J]. Journal of Tropical Crops, 2021, 5(25): 1-14. (in Chinese)

[40]
KIM S, PARK M, YEOM S. Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species[J]. Nature Genetics, 2014, 4(6): 270-278.

[41]
吕萌, 倪志勇, 王娟, 李波, 罗淑萍, 范玲. 棉花甲基化酶基因COMTCCoAOMT的组织特异性表达分析[J]. 核农学报, 2010, 24(4): 713-719.

LU M, NI Z Y, WANG J, LI B, LUO S P, FAN L. Tissue specific expression analysis of cotton methylase genes COMT and CCoAOMT[J]. Journal of Nuclear Agriculture, 2010, 24(4): 713-719. (in Chinese)

[42]
徐记迪. 柑橘全基因组DNA甲基化分析及调控作用研究[D]. 武汉: 华中农业大学, 2015.

XU J D. Analysis and regulation of genome-wide DNA methylation in Citrus[D]. Wuhan: Huazhong Agricultural University, 2015. (in Chinese)

[43]
TEYSSIER E, BERNACCHIA G, MAURY S, KIT A, STAMMITTI-BER L, ROLIN D, GALLUSCI P. Tissue dependent variations of DNA methylation and endoreduplication levels during tomato fruit development and ripening[J]. Planta, 2008, 228(3): 391-399.

DOI

[44]
HUANG H, LIU R, NIU Q F, TANG K, ZHANG B, ZHANG H, CHEN K S, ZHU J K, LANG Z B. Global increase in DNA methylation during orange fruit development and ripening[J]. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(4): 1430-1436.

[45]
XU J D, ZHOU S H, GONG X Q, SONG Y, WAN N S, MA F W, GUAN Q M. Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple[J]. Plant Biotechnology Journal, 2018, 16(2): 672-687.

DOI

[46]
ZHANG B, DENIES M, JIAO C, XU Y M, CHEN K M, FEI Z J, JAMES J, HARRY J K. Chilling-induced tomato flavor loss is associated with altered volatile synthesis and transient changes in DNA methylation[J]. Proceedings of the National Academy of Sciences, 2016, 113(49): 12580-12585.

DOI

Outlines

/