运用MSAP技术测定谷子基因组DNA胞嘧啶甲基化水平

张雨欣，阮亚男，赵宸，薛敏敏，李博，王晶晶，刘洋，王凯玺，王红艳



运用MSAP技术测定谷子基因组DNA胞嘧啶甲基化水平

张雨欣1*，阮亚男1*，赵宸1，薛敏敏1，李博1，王晶晶1，刘洋1，王凯玺2，王红艳1

1 辽宁大学 生命科学院 植物表观遗传与进化实验室，辽宁 沈阳 110036 2 辽宁省水土保持研究所，辽宁 朝阳 122000

张雨欣, 阮亚男, 赵宸, 等. 运用MSAP技术测定谷子基因组DNA胞嘧啶甲基化水平. 生物工程学报, 2019, 35(2): 263–269.ZhangYX, Ruan YN, Zhao C, et al. Analysis of genomic DNA methylation level in foxtail millet by Methylation Sensitive Amplified Polymorphism. Chin J Biotech, 2019, 35(2): 263–269.

DNA甲基化是真核生物一种重要的表观修饰形式。为了探讨谷子基因组DNA胞嘧啶甲基化的水平和模式，以谷子的两个品种朝谷58号和豫谷1号为实验材料，利用RⅠ和Ⅱ/Ⅰ双酶切建立适合于谷子基因组的甲基化敏感扩增多态性 (MSAP) 分析体系。结果表明，从100对MSAP选扩引物中，筛选出32对MSAP引物组合，在朝谷58号和豫谷1号中分别扩增产生1 615、1 482条清晰可辨且可重复的DNA条带，其中包括3种类型的甲基化条带，朝谷58号和豫谷1号的基因组中CCGG序列胞嘧啶甲基化水平分别为6.93%和8.77%。这种谷子不同品种间甲基化水平和分布位点的差异为从表观遗传学的角度培育新品种提供了初步的理论依据和参考。

DNA甲基化，谷子，甲基化敏感扩增多态性 (MSAP)

DNA甲基化是一种重要的表观遗传学修饰方式，与许多重要的生物过程相关，其中包括染色质结构维持、基因印记、细胞分化等[1]。DNA甲基化是指在DNA甲基转移酶 (DNA methyltransferase)的催化作用下，将S-腺苷甲硫氨酸 (SAM) 上的甲基转移到胞嘧啶的第5位碳原子上，形成5-甲基胞嘧啶 (5-mC)[2]。植物细胞中，DNA甲基化主要发生在CG、CHG和CHH位点上 (H指A、C、T三种碱基)，由不同的DNA甲基转移酶调 控[3]。DNA甲基化主要分布于着丝粒区和重复序列中，其中转座子常常被高度甲基化，使得植物的基因组保持相对稳定[4]。植物的甲基化水平会因为物种、器官、生长时期不同而不同，并以此来调节植物生长发育、器官分化甚至影响植物进化[1,4-7]。随着对DNA甲基化研究的不断深入，已经发展了很多方法来检测DNA甲基化的水平，其中MSAP技术以其多态性高、信息量大、操作简单等优点被广泛应用于植物DNA甲基化检测[7-10]。

甲基化敏感扩增多态性，即MSAP技术 (Methylation sensitive amplified polymorphism) 是由AFLP (Amplified fragment length polymorphism) 技术发展而来的[5]。其原理是运用限制性内切酶Ⅱ和Ⅰ识别CCGG序列，限制性内切酶RⅠ识别GAATTC序列，将Ⅱ和Ⅰ分别和RⅠ组合对基因组进行酶切，再加上相应的限制性内切酶的接头，根据接头序列设计引物，对CCGG位点甲基化进行特异性扩增[6-7]。虽然Ⅱ和Ⅰ两种限制性内切酶识别的位点相同，但两者对DNA甲基化敏感程度存在很大的差异，利用其敏感性的不同可以检测出总甲基化水平，并区分DNA的全甲基化和半甲基化等甲基化状态。已有众多研究者采用这一方法取得了一系列的研究成果，王聪聪等利用MSAP技术分析不同单株叶籽银杏DNA 甲基化水平与模式，统计发现不同单株DNA甲基化水平差异极显著，正常银杏DNA甲基化水平低于叶籽银杏DNA甲基化水平[8]。Zheng等利用MSAP技术评估水稻两个品种在幼苗期处于正常条件和干旱胁迫下的原代和第6代的DNA甲基化变化，发现29%干旱诱导的DNA甲基化变异位点在胁迫解除后仍然可以保持，并且在干旱条件下DNA甲基化的变异具有组织特异性，其中分蘖期比抽穗期和孕穗期的甲基化水平显著下降[9]。李照令等利用MSAP技术研究不同浓度Cd处理对拟南芥幼苗基因组DNA甲基化水平的影响，发现3个处理组的DNA甲基化水平均高于对照，并且随着处理浓度的增加，其甲基化水平随之降低[10]。

谷子又称为粟，是起源于我国的传统优势作物之一，一直在北方大面积种植。由于谷子具有抗旱耐瘠、水分利用率高、适应性强、营养价值高等优点，近年来受到研究人员的重视[11]。谷子的测序工作在2012年完成，基因组大小约为515 Mb，基因组序列相对较小，富含重复序列，并且有一些转录组测序也已完成，为功能基因的研究提供了良好条件[12-15]。但在表观遗传学的研究领域，关于谷子的研究还处于刚刚起步的阶段，谷子在正常生长条件下全基因组甲基化的水平尚未深入研究。因此，本研究以两个优势种——朝谷58号和豫谷1号为例，对谷子基因组DNA胞嘧啶甲基化水平进行评估，以了解谷子在正常生长条件下全基因组甲基化水平差异和模式特点，以期为进一步从表观遗传学的角度培育谷子新品种提供初步的理论依据和参考。

1 材料与方法

1.1 材料

朝谷58号和豫谷1号，分别由辽宁省水土保持研究所和中国农业科学院作物科学研究所提供。其中，豫谷1号是目前谷子遗传研究中应用最广泛的品种，也是基因组测序的品种[15-16]；朝谷58是适应辽西干旱半干旱地区自然条件的优秀品种，抗逆性强，且丰产稳产，是谷子抗逆性状基因挖掘的好材料[17]。两个品种均播种于灭菌土壤中，并在自然光照条件下生长，8周后选取叶片完全展开且长势一致的倒三叶、倒四叶叶片进行采集，液氮冷冻后在−80 ℃冰箱中保存。

1.2 方法

1.2.1 DNA提取

采用改良的CTAB方法提取谷子叶片的基因组DNA[18]。所用试剂购自北京鼎国昌盛生物技术有限责任公司。

1.2.2 DNA甲基化敏感多态性分析

分别用RⅠ/Ⅱ、RⅠ/Ⅰ组合对基因组DNA进行酶切，连接接头 (表1)，实验所用酶购买自NEB公司，所用引物序列由生工生物工程 (上海) 股份有限公司合成 (表2)。每对引物进行5%的聚丙烯酰胺凝胶电泳、银染后对清晰可辨、可重复的带型进行记录，并对其位点进行MSAP分析。实验根据泳道内有无条带，将清晰可重复的条带记为“1”，同一位置没有条带记为“0”。统计每对引物扩增出的总条带数和其中的甲基化条带数，每个引物组合做3次重复。

表1 MSAP接头

表2 MSAP扩增引物

Note:andrepresentRⅠandⅡ/Ⅰ, respectively.

酶切连接体系、预扩、选扩体系参照Wang等反应体系进行[19]。

PCR产物电泳及银染：PCR扩增后的产物经5%聚丙烯酰胺凝胶电泳及硝酸银染色后，统计凝胶上各种模式的片段[19]。

2 结果与分析

2.1 引物筛选和DNA甲基化图谱构建

本研究是通过使用RⅠ和Ⅱ/Ⅰ的引物组合扩增PCR片段，从而得到基因组甲基化“CCGG”位点的变异模式，引物组合由10个RⅠ引物和10个Ⅱ /Ⅰ引物自由组合而成，共100对组合。从中筛选出32对可以扩增出条带分布均匀、带多清晰的引物组合完成选择性扩增实验，以朝谷58号为例，选用的32对引物组合及其扩增条带的统计情况见表3。

MSAP中Ⅱ和Ⅰ是两种同裂酶，它们的酶切位点均为CCGG，但二者对酶切位点甲基化修饰的敏感性不同，其中Ⅱ可以识别并切开外侧胞嘧啶半甲基化位点，Ⅰ可以识别并切开内侧胞嘧啶全甲基化位点[20]。因此DNA样品经RⅠ/Ⅱ () 和RⅠ/Ⅰ () 酶切扩增后的产物在PAGE胶上会出现如下3种甲基化类型：Ⅰ型：+和+都有带，表明CCGG位点未发生甲基化；Ⅱ型：+有带，无带，表明CCGG位点处于半甲基化状态，即DNA甲基化只发生在一条链上，其互补链上没有；Ⅲ型：+无带，+有带，表明CCGG序列中的内侧CG位点处于全甲基化状态，即DNA甲基化同时发生在两条链的内侧胞嘧啶上。图1展示了朝谷58号和豫谷1号不同引物的选择性扩增图谱，表4统计了朝谷58号和豫谷1号的胞嘧啶甲基化水平。

表3 朝谷58号32对引物组合及其扩增条带统计表

Note:andrepresentRⅠandⅡ/Ⅰ,respectively.

图1 朝谷58号 (A) 和豫谷1号 (B) 不同引物组合选择性扩增图

2.2 甲基化修饰水平分析

2.2.1 朝谷58号的甲基化水平

从100对MSAP 选扩引物中，选出32对适合朝谷58号的MSAP引物组合，共扩增产生1 615条清晰可辨且可重复的DNA条带。从表4可以看出，辽宁省优势品种朝谷58在全基因组水平上CCGG位点胞嘧啶甲基化比例为6.93%，其中内侧全甲基化比例为4.27%，单链半甲基化比例为2.66%。从结果看出朝谷58号在CCGG位点上的胞嘧啶甲基化主要发生在内侧胞嘧啶上。

2.2.2 豫谷1号的甲基化水平

从100对MSAP 选扩引物中，选出32对适合于豫谷1号的MSAP 引物组合，共扩增产生1 482条清晰可辨且可重复的DNA条带。在全基因组水平上CCGG位点的甲基化比例为8.77%，根据条带类型统计，其中内侧胞嘧啶全甲基化比例为6.41%，半甲基化比例为2.36%，豫谷1号在CCGG位点上的胞嘧啶甲基化也主要发生在内侧胞嘧啶上。

2.2.3 朝谷58号与豫谷1号全基因组甲基化水平的比较分析

在计算全基因组CCGG位点上的胞嘧啶甲基化水平后 (表4)，发现朝谷58号的甲基化水平低于豫谷1号，这种甲基化水平的差异是由于两个品种胞嘧啶内侧全甲基化位点差异造成的。同时，两个品种的胞嘧啶内侧全甲基化水平均高于单链半甲基化水平，说明在谷子基因组中这种对称性的胞嘧啶内侧全甲基化位点是甲基化修饰的主要形式。此外，单链胞嘧啶半甲基化水平在两个品种中没有明显差异。仪治本等[21]在运用MSAP技术研究玉米杂交种及其亲本基因组DNA胞嘧啶甲基化水平中，其选用的第一组亲本Mo 17与Suwan 5的全甲基化水平分别为17.1%和24.4%，第二组亲本Suwan 1与太3221的全甲基化水平分别为24.6%和21.6%，与本实验的研究结果相同，说明同种作物的不同品种之间的全甲基化水平也不同。

表4 朝谷58号和豫谷1号的胞嘧啶甲基化水平

Note: total methylated rate (%) = [(Ⅱ+Ⅲ)/(Ⅰ+Ⅱ+Ⅲ)]×100; fully methylated ratio(%) = [(Ⅲ)/(Ⅰ+Ⅱ+Ⅲ)]×100; hemi methylated ratio(%) = [(Ⅱ)/(Ⅰ+Ⅱ+Ⅲ)]×100; non methylated ratio(%) = [(Ⅰ)/(Ⅰ+Ⅱ+Ⅲ)]×100.

3 讨论

DNA甲基化在植物的生长发育中具有重要意义，甲基化的变化会影响农作物的表型、花期和产量[22-23]。另外，Zhong等研究表明在番茄果实成熟过程中DNA甲基化的变化发挥重要作用[24]。本研究通过100对引物组合筛选获得32对有效的引物组合，检测两个谷子品种在CCGG位点DNA甲基化水平。朝谷58号和豫谷1号分别检测到1 615条和1 482条条带，但由于MSAP技术使用的限制性内切酶只能识别和切割CCGG序列，对其他位点的甲基化不能检测，其检测位点具有局限性，所以对于基因组中实际甲基化水平评估偏低[20]，即实验数据得到朝谷58号和豫谷1号全基因组CCGG序列胞嘧啶甲基化所占的比例为6.93%和8.77%，要低于谷子两个品种各自的甲基化水平。Xiong等在研究水稻杂交种及其亲本的甲基化模式中发现，其中亲本Zhenshan 97和Minghui 63的胞嘧啶甲基化水平基本相同为16.3%，而杂交后代的Shanyou 63的胞嘧啶甲基化水平低于两种亲本，说明其杂交后代和亲本之间的甲基化水平也会存在差异[25]。谷子不同品种间的甲基化差异可能是其表型和性状差异原因，这种差异可以为培育新型的谷子品种提供表观遗传学方向的思路。

[1] Richards EJ. DNA methylation and plant development. Trends Genetics, 1997, 13(8): 319–323.

[2] Feng SH, Jacobsen SE, Reik W. Epigenetic reprogramming in plant and animal development. Science, 2010, 330(6004): 623–627.

[3] Finnegan EJ, Kovac KA. Plant DNA methyltransferases. Plant Mol Biol, 2000, 43(2/3): 189–201.

[4] Singer T, Yordan C, Martienssen RA. Robertson’s mutator transposons in.are regulated by the chromatin-remodeling gene decrease in DNA methylation (DDM1). Genes Dev, 2001, 15(5): 591–602.

[5] Wang ZC, Li ZA, Li SP. The technology of MSAP and its applications in plant. Biotechnol Bull, 2006(S1): 195–196 (in Chinese). 王子成, 李忠爱, 李锁平. MSAP技术及其在植物上的应用. 生物技术通报, 2006(S1): 195–196.

[6] Han YN, Zhao XJ, Cai L. Application of methylation-sensitive amplified polymorphism technology in plant stress resistance. Biotechnol Bull, 2010(6): 71–74, 79 (in Chinese). 韩雅楠, 赵秀娟, 蔡禄. MSAP技术在植物抗逆性方面的应用. 生物技术通报, 2010(6): 71–74, 79.

[7] Li WG, Chang TJ, Gong HM. Methylation-sensitive amplified polymorphism and its application to plant genetics. Biotechnology, 2008, 18(1): 83–87 (in Chinese).李卫国, 常天俊, 龚红梅. MSAP技术及其在植物遗传学研究中的应用. 生物技术, 2008, 18(1): 83–87.

[8] Wang CC, Xing SY, Li JH, et al. DNA methlation levels and patterns in different ginkgoL. var. Epiphylla Mak. by methlation-sensitive amplification polymorphism (MSAP). J Nucl Agric Sci, 2013, 27(4): 399–407 (in Chinese). 王聪聪, 邢世岩, 李际红, 等. 不同单株叶籽银杏DNA甲基化水平与模式的MSAP分析. 核农学报, 2013, 27(4): 399–407.

[9] Zheng XG, Chen L, Li MS, et al. Transgenerational variations in DNA methylation induced by drought stress in two rice varieties with distinguished difference to drought resistance. PLoS ONE, 2013, 8(11): e80253.

[10] Li ZL, Wang HT, Chen HJ, et al. Studying genomic methylation of arabidopsis thaliana seedlings under cadmium stress using MSAP. J Agro-Environ Sci, 2014, 33(1): 28–36 (in Chinese). 李照令, 王鹤潼, 陈瑞娟, 等. 运用MSAP研究镉胁迫对拟南芥幼苗基因甲基化的影响. 农业环境科学学报, 2014, 33(1): 28–36.

[11] Xu ZG, Fu LB, Zhong FN. Analysis on the regional comparative advantage of grain production in China. Chin Agric Resour Region Plann, 2001, 22(1): 45–48 (in Chinese). 徐志刚, 傅龙波, 钟甫宁. 中国粮食生产的区域比较优势分析. 中国农业资源与区划, 2001, 22(1): 45–48.

[12] Zhang GY, Liu X, Quan ZW, et al. Genome sequence of foxtail millet () provides insights into grass evolution and biofuel potential. Nat Biotechnol, 2012, 30(6): 549–554.

[13] Jo Y, Lian S, Cho JK, et al.transcriptome assembly ofvariety Taejin. Genomics Data, 2016, 8: 121–122.

[14] Tang S, Li L, Wang YQ, et al. Genotype-specific physiological and transcriptomic responses to drought stress in(an emerging model for Panicoideae grasses). Sci Rep, 2017, 7: 10009.

[15] Bennetzen JL, Schmutz J, Wang H, et al. Reference genome sequence of the model plant Setaria. Nature Biotechnol, 2012, 30(6): 555–561.

[16] Li W, Zhi H, Wang YF, et al. Analysis on the conditions of EMS treatment on Foxtail Millet. J Hebei Agric Sci, 2010, 14(11): 77–79 (in Chinese). 李伟, 智慧, 王永芳, 等. 谷子EMS诱变的处理条件分析. 河北农业科学, 2010, 14(11): 77–79.

[17] Wang SS, Zhang N, Wang KX, et al. Genetic diversity of foxtail millet cultivars in western Liaoning province revealed by SSR markers. Mol Plant Breed, 2015, 13(5): 1091–1097 (in Chinese). 王姗姗, 张宁, 王凯玺, 等. 中国辽西地区谷子品种遗传多样性的SSR分析. 分子植物育种, 2015, 13(5): 1091–1097.

[18] Kidwell KK, Osborn TC. Simple plant DNA isolation procedures//Beckmann JS, Osborn TC. Plant Genomes: Methods for Genetic and Physical Mapping. Netherlands: Springer 1992: 1–13.

[19] Wang QM, Wang YZ, Sun LL, et al. Direct and indirect organogenesis ofand assessment of DNA methylation changes in various regenerated plantlets. Plant Cell Rep, 2012, 31(7): 1283–1296.

[20] Dong ZY, Wang YM, Zhang ZJ, et al. Extent and pattern of DNA methylation alteration in rice lines derived from introgressive hybridization of rice andGriseb. Theoret Appl Genet, 2006, 113(2): 196–205.

[21] Yi ZB, Sun Y, Niu TT, et al. Genomic DNA cystoine methylations of corn hybrids and their parents. Acta Bot Bor-Occidental Sin, 2005, 25(12): 2420–2425 (in Chinese). 仪治本, 孙毅, 牛天堂, 等. 玉米杂交种及其亲本基因组胞嘧啶甲基化水平研究. 西北植物学报, 2005, 25(12): 2420–2425.

[22] Gehring M, Henikoff S. DNA methylation dynamics in plant genomes. Biochim Biophys Acta Gene Struct Expr, 2007, 1769(5/6): 276–286.

[23] Sung S, He YH, Eshoo TW, et al. Epigenetic maintenance of the vernalized state inrequires LIKE HETEROCHROMATIN PROTEIN 1. Nat Genetics, 2006, 38(6): 706–710.

[24] Zhong SL, Fei ZJ, Chen YR, et al. Single-base resolution methylomes of tomato fruit development reveal epigenome modifications associated with ripening. Nat Biotechnol, 2013, 31(2): 154–159.

[25] Xiong LZ, Xu CG, Maroof MAS, et al. Patterns of cytosine methylation in an elite rice hybrid and its parental lines, detected by a methylation-sensitive amplification polymorphism technique. Mol Gen Genet, 1999, 261(3): 439–446.

Analysis of genomic DNA methylation level in foxtail millet by Methylation Sensitive Amplified Polymorphism

Yuxin Zhang1*, Yanan Ruan1*, Chen Zhao1, Minmin Xue1, Bo Li1, Jingjing Wang1, Yang Liu1, Kaixi Wang2, and Hongyan Wang1

1 Laboratory of Plant Epigenetics and Evolution, School of Life Sciences, Liaoning University, Shenyang 110036, Liaoning, China 2 Liaoning Institute of Soil and Water Conservation, Chaoyang 122000, Liaoning, China

DNA methylation is an important type of epigenetic modification in eukaryotes. In order to research genome-wide methylation levels and patterns in foxtail millet (), the Methylation Sensitive Amplified Polymorphism (MSAP) analysis (employing double digestion withR I andIII) was established and applied in two foxtail millet cultivars (Chaogu 58 and Yugu 1). The results showed that 32 pairs of MSAP primers were selected from 100 MSAP primers, and 1 615 and 1 482 clearly distinguishable and reproducible bands were amplified from Chaogu 58 and Yugu 1 respectively, including 3 types of methylation patterns. Cytosine methylation levels of CCGG context in Chaogu 58 and Yugu 1 were characterized as 6.93% and 8.77% respectively. Such different genomic DNA methylation levels between two foxtail millet varieties may provide a preliminary reference for the cultivation of this crop from a novel epigenetic viewpoint.

DNA methylation, foxtail millet (), methylation-sensitive amplification polymorphism (MSAP)

May 25, 2018;

August 13, 2018

National Natural Science Foundation of China (No. 31100172), Program for Liaoning Excellent Talents in University (No. LJQ2013003), Guidance Program of Natural Science Fundation, Department of Science and Technology, Liaoning Province (No. 20180550686), and Project in Shenyang Science and Technology Bureau (No. F16-205-1-49).

Hongyan Wang. Tel: +86-24-62202232; E-mail: hongyan2003@126.com

10.13345/j.cjb.180220

*These authors contributed equally to this work.

国家自然科学基金 (No. 31100172)，辽宁省高等学校优秀人才支持计划 (No. LJQ2013003)，辽宁省科技厅自然基金指导计划 (No. 20180550686)，沈阳市科技局项目 (No. F16-205-1-49) 资助。

(本文责编 郝丽芳)