華南植物園錐栗F睲SAP采樣策略及遺傳多樣性分析

劉孟　歐陽(yáng)學(xué)軍　陳潔　全慧敏　葉萬(wàn)輝　劉衛(wèi)

摘要： 基因組甲基化修飾受環(huán)境因素的影響。在以甲基化為代表的表觀遺傳學(xué)研究中，如何減少保存環(huán)境對(duì)異地采后樣品的影響，提高整個(gè)實(shí)驗(yàn)的準(zhǔn)確性和科學(xué)性，目前尚未有系統(tǒng)的認(rèn)知。該研究選取5種常用的采后樣品保存方式（液氮冷凍、-20 ℃冷凍、變色硅膠干燥、密封袋密封、75%酒精浸泡），分別用Wilcoxon signed ranks tests統(tǒng)計(jì)分析和UPGMA聚類(lèi)分析方法，對(duì)華南植物園錐栗進(jìn)行FMSAP研究，以期找出最佳保存方式。同時(shí)，利用正交試驗(yàn)法對(duì)FMSAP體系進(jìn)行優(yōu)化，篩選出9對(duì)引物（E3H/M2；E5H/M2；E6H/M1；E6H/M5；E8H/M1；E8H/M5；E9H/M2；E11H/M5；E14H/M1），并對(duì)不同發(fā)育時(shí)期的錐栗甲基化水平及遺傳多樣性進(jìn)行了論述。結(jié)果表明：在錐栗FMSAP的研究中，Willcoxon signed ranks tests統(tǒng)計(jì)分析和UPGMA聚類(lèi)分析結(jié)論一致，密封袋保存效果最佳；成熟葉半甲基化率（27.83%）和總甲基化率（51.13%）高于幼葉（21.35%，45.90%），全甲基化率（23.30%）低于幼葉（24.55%），平均多態(tài)位點(diǎn)百分?jǐn)?shù)39.60%，香農(nóng)信息指數(shù)0.207±0.002，表現(xiàn)出較高的甲基化水平和遺傳多樣性。

關(guān)鍵詞： 錐栗， 甲基化， FMSAP， 引物篩選， 樣品保存

中圖分類(lèi)號(hào)： Q943

文獻(xiàn)標(biāo)識(shí)碼： A

文章編號(hào)： 10003142（2017）01001508

Abstract： Genomic methylation characteristics were influenced by environmental factors. In the study of epigenetics， represented by DNA methylation， the preservation environmental of postharvest samples， especially the samples in the remote place， would have an influence on the later experiment. In this reason， standard sampling and preserve operation were crucial in the field of epigenetic. But the related research was few. Therefore， we studied Castanopsis chinensis in South China Botanical Garden with FMSAP sampling strategy using wilcoxon signed rank tests and unweighted pairgroup method with arithmetic means（UPGMA） methods， in order to find out the best way of sample preservation in five different environments （Liquid nitrogen， -20 ℃， allochroic silicagel， hermetic bag and 75% alcohol）. Using the orthogonal experiment method， we optimized the FMSAP system and screened nine pairs of primers （E3H/M2， E5H/M2， E6H/M1， E6H/M5， E8H/M1， E8H/M5， E9H/M2， E11H/M5 and E14H/M1）. Additionally， we briefly described different periods of methylation level and genetic diversity. The results showed that hermetic bag preservation was the most appropriate way. And the half and total methylation rate of mature leaves was 27.83% and 51.13%， higher than that of young leaves. The full methylation rate of mature leaves was 23.30%， lower than that of young leaves. The average percentage of polymorphic loci was 39.60%. Shannon information index was 0.207 ± 0.002， which showed a higher level of methylation and genetic diversity.

Key words： Castanopsis chinensis， methylation， FMSAP， primer selection， sample preservation

目前已有研究證據(jù)證明DNA甲基化這種重要的表觀遺傳方式，在植物生長(zhǎng)發(fā)育及進(jìn)化過(guò)程中起著重要的調(diào)節(jié)作用， 并顯著影響所在DNA區(qū)域的遺傳變異速率（Johannes et al， 2009；Paun et al， 2010；Angers et al， 2010）。近年來(lái)，DNA甲基化已成為生物學(xué)研究的熱點(diǎn)之一，特別在植物學(xué)研究上，由于植物固著生長(zhǎng)的特性，使之發(fā)展出以表型可塑性為主的適應(yīng)性機(jī)制，因此研究DNA甲基化對(duì)植物轉(zhuǎn)錄表達(dá)以及表型適應(yīng)性變化的影響已經(jīng)成為植物生理生態(tài)適應(yīng)機(jī)制和進(jìn)化的主要手段之一（高樂(lè)旋等，2008）。但是，相對(duì)于基因組序列特征，DNA甲基化的選擇性可遺傳特征對(duì)實(shí)驗(yàn)研究方法的影響卻極少被注意和研究。

從目前的研究看，一方面如低溫誘導(dǎo)的金魚(yú)草（Antirrhinum majus）的轉(zhuǎn)座子Tam 3序列的甲基化水平下降，在一個(gè)世代內(nèi)可逆地響應(yīng)溫度變化（Hashinda et al， 2003）；干早脅迫誘導(dǎo)水稻DNA甲基化水平發(fā)生變化，復(fù)水后有70%回復(fù)到初始狀態(tài)（Wang et al， 2010）；煙草異染色質(zhì)DNA受到脅迫后甲基化水平上升，當(dāng)脅迫移除后部分可逆地去甲基化（Kovar et al，1997）。另一方面，許多研究表明，表觀遺傳變異與遺傳變異有緊密的聯(lián)系，且可遺傳的表觀遺傳變異具有適應(yīng)性，能夠影響適應(yīng)的速度和結(jié)果，有助于微進(jìn)化（Kaeppler & Phillips， 1993； Stokes et al， 2002； Klironomos et al， 2013；Herrera et al， 2013）。許多由環(huán)境改變或脅迫引發(fā)的植物DNA甲基化，通過(guò)有絲分裂和減數(shù)分裂，可以在一些植物中遺傳數(shù)代（Calarco et al， 2012；Klose & Bird， 2006；Heard & Martienssen， 2014；Morgan et al， 2005； Paszkowski & Grossniklaus， 2011）。

甲基化的選擇性可遺傳性對(duì)具體的試驗(yàn)研究帶來(lái)了特殊的要求?；蚪M內(nèi)，許多位點(diǎn)上甲基化特征的穩(wěn)定性較低，容易在短時(shí)間內(nèi)因環(huán)境因素而變化。尤其在木本植物的生態(tài)適應(yīng)機(jī)制的研究中，因采樣地點(diǎn)偏遠(yuǎn)，采樣量大等現(xiàn)實(shí)問(wèn)題，無(wú)法立即對(duì)DNA進(jìn)行提取。而在以甲基化為代表的表觀遺傳研究中，樣品采后保存環(huán)境對(duì)后期實(shí)驗(yàn)具有影響。因此，研究如何進(jìn)行規(guī)范的采樣和樣品保存操作具有至關(guān)重要的現(xiàn)實(shí)意義。

錐栗作為森林生態(tài)系統(tǒng)的優(yōu)勢(shì)種，在生物多樣化格局形成中有著舉足輕重的地位（王崢?lè)宓龋?001）。尤其是從微觀角度研究錐栗遺傳多樣性分布格局和生理適應(yīng)性調(diào)節(jié)機(jī)制已然成為了研究熱點(diǎn)（王崢?lè)宓龋?004）。但在研究初期，錐栗樣品采集后，尤其是異地采樣后，如何科學(xué)完好的保存樣品，最大可能的減少外界環(huán)境因素對(duì)DNA甲基化的影響，則鮮有研究報(bào)道。因此，本文重點(diǎn)就華南植物園錐栗葉片不同保存方式對(duì)FMSAP的影響，不同生長(zhǎng)階段甲基化水平及遺傳多樣性信息進(jìn)行了系統(tǒng)研究，以期為研究錐栗甲基化修飾制定科學(xué)合理的采樣策略和遺傳多樣性分析提供理論指導(dǎo)。

1材料與方法

1.1 樣品采集

以華南植物園內(nèi)25棵錐栗作為研究對(duì)象，分別進(jìn)行樣品保存試驗(yàn)和老幼葉試驗(yàn)。

樣品保存試驗(yàn)：每棵錐栗采集5片健康成熟葉片，共125片。葉片沿中脈一分為二，一半為實(shí)驗(yàn)組，一半為對(duì)照組。老幼葉試驗(yàn)：任選3棵胸徑相近個(gè)體，每棵采3片成熟葉和3片幼葉，共18片。

采后樣品均立即帶回實(shí)驗(yàn)室，進(jìn)行下一步研究。

1.2 DNA提取

樣品保存試驗(yàn)：5個(gè)試驗(yàn)組（液氮冷凍處理組、-20 ℃冷凍處理組、變色硅膠干燥處理組、密封袋密封處理組、75%酒精浸泡組），處理24 h后，用改良的CTAB法提取DNA（陳昆松等，2004），對(duì)照組及老幼葉試驗(yàn)葉片，采后立即提取DNA。

Nano Drop 2000/2000c 超微量紫外可見(jiàn)分光光度計(jì)和1%瓊脂糖凝膠電泳對(duì)DNA進(jìn)行篩選檢測(cè)，-20 ℃冷凍保存。

1.3 FMSAP

熒光標(biāo)記甲基化敏感擴(kuò)增多態(tài)性（fluorescent labeled methylation sensitive amplified polymorphism， FMSAP）是在甲基化敏感擴(kuò)增多態(tài)性（methylation sensitive amplification polymorphism， MSAP）的基礎(chǔ)上對(duì)選擇性引物進(jìn)行了熒光標(biāo)記。FMSAP包括4個(gè)反應(yīng)體系，酶切體系、連接體系、預(yù)擴(kuò)增體系和選擇性擴(kuò)增體系（褚會(huì)娟等，2011）。通過(guò)梯度設(shè)置參數(shù)，優(yōu)化體系中的酶切時(shí)間、連接時(shí)間、退火溫度及時(shí)間等因素，結(jié)合正交試驗(yàn)法，對(duì)預(yù)擴(kuò)增體系和選擇性擴(kuò)增體系進(jìn)行優(yōu)化，得到錐栗FMSAP技術(shù)的最優(yōu)反應(yīng)體系。接頭序列見(jiàn)表1。

酶切體系：50 μL（DNA 500 ng，10×Buffer 5 μL，EcoR Ⅰ 10 U，Hpa Ⅱ/Msp Ⅰ 10 U，10 mg·mL1 100×BSA 0.2 μL），37 ℃恒溫孵育6 h，65 ℃變性10 min。

連接體系：50 μL（酶切液12.5 μL，5 pmol·μl1 EcoR Ⅰ接頭1 μL，50 pmol·μl1 Hpa Ⅱ/Msp Ⅰ接頭1 μL， 5 U·μL1 T4連接酶0.5 μL，10×T4 Buffer 5 μL），16 ℃連接過(guò)夜，65 ℃變性8 min（Vos et al，1995）。

預(yù)擴(kuò)增體系：25 μL（5 U·μL1 Taq聚合酶0.4 μL，10×PCR Buffer 2.5 μL，模板5 μL，2.5 pmol·μL1 dNTPs 1.6 μL，10 pmol·μL1引物各1 μL），94 ℃預(yù)變性4 min；94 ℃變性30 s，56 ℃退火30 s，72 ℃延伸1 min，25個(gè)循環(huán)；72 ℃延伸10 min。

選擇性擴(kuò)增體系：25 μL（5U·μL1 Taq聚合酶0.2 μL，2.5 pmol·μL1 dNTPs 1.6 μL，10×PCR Buffer 2.5 μL，10 pmol·μL1引物各3 μL），94 ℃預(yù)變性4 min；94 ℃變性30 s，65 ℃退火30 s，每個(gè)循環(huán)退火溫度降低0.7 ℃，72 ℃延伸1 min，12個(gè)循環(huán)；94 ℃變性30 s，56 ℃退火30 s，72 ℃延伸1 min，23個(gè)循環(huán)；72 ℃延伸10 min。

1.4 引物篩選

在EcoR Ⅰ、Hap Ⅱ/Msp Ⅰ通用引物的基礎(chǔ)上，通過(guò)添加隨機(jī)堿基，組合出128個(gè)引物對(duì)（Xiong et al，1999）。用優(yōu)化的FMSAP體系進(jìn)行多態(tài)性擴(kuò)增，產(chǎn)物先用1%瓊脂糖凝膠電泳初步篩選，再用6%變性聚丙烯酰胺凝膠電泳篩選，最終選出9對(duì)可擴(kuò)增出完整清晰且多態(tài)性高條帶的引物對(duì)（表2，圖1）。

1.5 數(shù)據(jù)分析

利用GeneMarker2.2.0軟件讀取熒光標(biāo)記引物擴(kuò)增產(chǎn)物的高效毛細(xì)管電泳（High performance capillary electrophoresis，HPCE）結(jié)果（錢(qián)曉偉等，2014；李金龍等，2014），得到表示擴(kuò)增片段大小的熒光峰值圖，選取30～500 bp的片段進(jìn)行下一步的研究。

利用半?yún)^(qū)間放置法將表示片段大小的數(shù)據(jù)轉(zhuǎn)換成0/1二元數(shù)據(jù)矩陣，即某一特定位點(diǎn)，有條帶記為“1”，無(wú)條帶記為“0”（陳秋博等，2014；夏寒冰等，2009）。根據(jù)Hpa Ⅱ/EcoR Ⅰ和Msp Ⅰ/EcoR Ⅰ雙酶切位點(diǎn)帶型差異，將基因組甲基化模式分為3種：Ⅰ型，記作“11”，雙鏈內(nèi)胞嘧啶全甲基化；Ⅱ型，記作“10”，單鏈外胞嘧啶半甲基化；Ⅲ型，記作“01”，胞嘧啶未甲基化（Salmon et al，2008）。

根據(jù)Gervera et al（2002）的方法，把MSAP二元矩陣分為甲基化敏感多態(tài)矩陣（Methylation Sensitive Polymorphisms，MSP），即將原矩陣中“01”“10”的位點(diǎn)記為“1”，其它記為“0”，以及甲基化不敏感多態(tài)矩陣（Methylation Insensitive Polymorphisms，MISP），即將原矩陣中“11”位點(diǎn)記為“1”，其它記為“0”。

接著，將原始數(shù)據(jù)進(jìn)行l(wèi)og函數(shù)歸一化標(biāo)準(zhǔn)處理，KolmogorovSmirnova檢驗(yàn)，分析實(shí)驗(yàn)組與對(duì)照組擴(kuò)增片段數(shù)是否符合正態(tài)分布規(guī)律。如果符合，實(shí)驗(yàn)組與對(duì)照組兩兩配對(duì)，用參數(shù)檢驗(yàn)法進(jìn)一步分析，否則，選擇非參數(shù)檢驗(yàn)法（郭東星等，2000）。

2結(jié)果與分析

2.1 不同樣品保存方式對(duì)錐栗甲基化影響的統(tǒng)計(jì)分析

從篩選出的9對(duì)引物中選取3對(duì)（E3H/M2，E5H/M2，E9H/M2）分別對(duì)5種樣品保存處理組的DNA進(jìn)行FMSAP擴(kuò)增?；诨蚪M甲基化帶型分類(lèi)，統(tǒng)計(jì)實(shí)驗(yàn)組與對(duì)照組各帶型的擴(kuò)增條帶數(shù)，然后通過(guò)兩兩配對(duì)的統(tǒng)計(jì)分析方法， 找出實(shí)驗(yàn)組與對(duì)

照組差異最小的處理組，選出最優(yōu)樣品保存方式。KolmogorovSmirnova檢驗(yàn)結(jié)果（表3）顯示，30組數(shù)據(jù)中有18組不符合正態(tài)分布規(guī)律。因此，實(shí)驗(yàn)組與對(duì)照組配對(duì)差異性分析選用非參數(shù)檢驗(yàn)法，Wilcoxon signed ranks tests，檢驗(yàn)結(jié)果見(jiàn)表4。

Wilcoxon signed ranks tests 檢驗(yàn)結(jié)果顯示，與自身對(duì)照組相比，用75%酒精保存的錐栗葉片，3種甲基化模式擴(kuò)增條帶數(shù)均差異極顯著；用-20 ℃冷凍保存的，半甲基化、全甲基化擴(kuò)增條帶數(shù)差異極顯著，未甲基化無(wú)顯著差異；用變色硅膠和液氮保存的，半甲基化和全甲基化擴(kuò)增條帶數(shù)差異不顯著，未甲基化差異極顯著；用密封袋保存的，3種甲基化模式均無(wú)差異。

5種樣品保存方式對(duì)DNA甲基化修飾的影響程度由小到大依次為密封袋保存、液氮保存和變色硅膠保存、-20 ℃冷凍保存、75%酒精保存。初步判定最佳樣品保存方式是密封袋保存。

2.2 不同樣品保存方式對(duì)錐栗甲基化影響的聚類(lèi)分析

利用非參數(shù)檢驗(yàn)法，統(tǒng)計(jì)分析不同樣品保存方式下擴(kuò)增條帶數(shù)差異，可以初步判定不同處理方式對(duì)DNA甲基化修飾影響的大小程度，但無(wú)法排除其中個(gè)體差異和位點(diǎn)差異。利用個(gè)體位點(diǎn)信息進(jìn)行聚類(lèi)分析，將相似組歸為一類(lèi)，相異組進(jìn)行區(qū)分，可以解決這個(gè)問(wèn)題，驗(yàn)證上述結(jié)論。

將5種處理組的實(shí)驗(yàn)組與對(duì)照組各視為一類(lèi)，根據(jù)基因頻率得到 MSP矩陣的Nei氏遺傳距離，用MEGA軟件進(jìn)行非加權(quán)組（Unweighted pairgroup method with arithmetic means， UPGMA）聚類(lèi)分析（圖2）（Zhang et al，2015）。結(jié)果顯示，隨著遺傳距離的增加，每種處理的實(shí)驗(yàn)組與對(duì)照組先后聚在一起，依據(jù)聚在一起的先后順序，判定不同樣品保存方式對(duì)甲基化影響。影響程度由小到大依次為密封袋保存、變色硅膠保存、液氮保存，-20 ℃冷凍保存、75%酒精保存。與上述Wilcoxon signed ranks tests結(jié)果基本一致。最終得出結(jié)論，利用FMSAP技術(shù)對(duì)錐栗甲基化修飾進(jìn)行研究時(shí)，密封袋保存采后葉片樣品，能最大程度的減少對(duì)后期實(shí)驗(yàn)的影響，為最佳保存方式。

2.3 錐栗不同發(fā)育階段甲基化特異性分析

錐栗成熟葉和幼葉分別進(jìn)行FMSAP擴(kuò)增，成熟葉擴(kuò)增出2 828條，幼葉3 377條。成熟葉半甲基化率（27.83%）和總甲基化率（51.13%）均高于幼葉（21.35%，45.90%），而全甲基化率（23.30%）低于幼葉（24.55%）（表5）。同一個(gè)體不同發(fā)育階段甲基化水平的差異，表明錐栗在生長(zhǎng)發(fā)育過(guò)程中，生長(zhǎng)環(huán)境的變化可能誘導(dǎo)了甲基化調(diào)控，出現(xiàn)了一定程度的去甲基化過(guò)程（孫慧敏，2013；郭廣平，2011）。

2.4 錐栗甲基化模式與遺傳多樣性分析

用樣品保存試驗(yàn)對(duì)照組數(shù)據(jù)，對(duì)華南植物園錐栗甲基化模式與遺傳多樣性進(jìn)行分析。擴(kuò)增產(chǎn)物中共檢測(cè)出471個(gè)位點(diǎn)，總擴(kuò)增條帶數(shù)56 488條，其中全甲基化條帶數(shù)15 920條，半甲基化條帶數(shù)13 340條，總甲基化率51.80%，平均多態(tài)位點(diǎn)百分?jǐn)?shù)為39.60%，香農(nóng)信息指數(shù)0.207±0.002，說(shuō)明華南植物園錐栗具有較高的甲基化水平和遺傳多樣性（表6）。

GenALEx軟件得出錐栗個(gè)體遺傳距離的范圍為71～175，通過(guò)設(shè)置重復(fù)個(gè)體，得到Nei Genetic Distance的范圍為 0.031～0.205（高寶嘉等，2010）。將MISP矩陣與MSP矩陣進(jìn)行 Mantel檢驗(yàn)，r=0.806，P=0.001，基因組遺傳變異與表觀遺傳變異顯著相關(guān)（覃光蓮和譚勁英，2014）。

利用Pearson相關(guān)分析將全甲基化率、半甲基化率、未甲基化率、總甲基化率分別與胸徑進(jìn)行相關(guān)檢驗(yàn)，相關(guān)系數(shù)分別0.215、-0.012、-0.14、0.131，即全甲基化率和總甲基化率與胸徑呈正相關(guān)，未甲基化率、半甲基化率與胸徑呈負(fù)相關(guān)。

3討論

植物從種子萌發(fā)到死亡的整個(gè)生命過(guò)程中，生物或非生物因素的誘導(dǎo)均伴隨著DNA甲基化的特異性改變（Vanyushin & Ashapkin， 2011）。DNA甲基化通常在植物中扮演兩個(gè)角色。一是保護(hù)基因組抵御轉(zhuǎn)座子的干擾，進(jìn)而調(diào)控基因的表達(dá)（Chan et al， 2005）。在逆境脅迫下，通過(guò)DNA甲基轉(zhuǎn)移酶的作用（Goll & Bestor， 2005）和相關(guān)機(jī)制，如RNA干擾（RNAi）的介導(dǎo)（Lukens & Zhan， 2007），甲基基團(tuán)可以迅速、可逆地對(duì)植物基因組DNA進(jìn)行修飾，避免過(guò)渡的、不必要的基因重組與種群多樣性（Boyko & Kovalchuk， 2008）。二是在物種受到環(huán)境脅迫時(shí)，DNA甲基化模式對(duì)不同逆境脅迫表現(xiàn)出不同的水平響應(yīng)（Verhoeven et al， 2010），使植物在環(huán)境脅迫下作出快速應(yīng)答（Peng & Zhang， 2009； Fulnecˇek & Kovarˇík， 2014），從而根據(jù)發(fā)育階段、環(huán)境誘因的變化對(duì)基因表達(dá)模式進(jìn)行“可塑性”調(diào)節(jié)（Boyko & Kovalchuk， 2008； Grant-Downton & Dickinson， 2006； Rapp & Wendel； 2005）。

DNA甲基化所扮演的兩個(gè)角色使得其MSAP帶譜比AFLP蘊(yùn)含更加復(fù)雜的機(jī)理。本研究從擴(kuò)增產(chǎn)物位點(diǎn)差異和擴(kuò)增產(chǎn)物總量差異兩個(gè)水平，利用統(tǒng)計(jì)分析和聚類(lèi)分析的方法，分別分析了不同樣品保存方式之間FMSAP的差異。非參數(shù)檢驗(yàn)法Wilcoxon signed ranks tests和UPGMA聚類(lèi)結(jié)果表明，采后樣品保存方式顯著影響FMSAP的帶譜特征，24 h內(nèi)的常溫密封保存對(duì)樣品的影響最小，而其它保存條件顯著影響表觀遺傳修飾。

雖然目前有不少研究表明，甲基化遺傳變異標(biāo)記（如FMSAP）可用于區(qū)分種內(nèi)遺傳差異。但在實(shí)際分析FMSAP結(jié)果時(shí)，情況可能更加復(fù)雜，特別在多年生木本植物中，我們的研究表明不同發(fā)育階段錐栗甲基化修飾的特異性以及胸徑與不同甲基化模式之間存在相關(guān)性（結(jié)果尚未公布）。今后，在對(duì)研究結(jié)果分析時(shí)，這些均值得重點(diǎn)關(guān)注。

植物在不同生長(zhǎng)發(fā)育階段，通過(guò)表觀遺傳修飾介導(dǎo)基因表達(dá)，進(jìn)而適應(yīng)環(huán)境的不斷變化。通常認(rèn)為，DNA甲基化水平的升高通常抑制基因表達(dá)，而去甲基化可促進(jìn)基因的表達(dá)（Suzuki & Bird， 2008）。Sherman & Talbert（2002）研究發(fā)現(xiàn)，春化可以通過(guò)降低冬小麥的甲基化水平誘導(dǎo)冬小麥開(kāi)花。馬開(kāi)峰等（2012）研究發(fā)現(xiàn)，毛白楊子代和親代之間甲基化水平和模式存在較大的差異；本研究分別對(duì)錐栗成熟葉和幼葉甲基化水平和模式進(jìn)行分析，發(fā)現(xiàn)成熟葉擴(kuò)增條帶數(shù)（2 828條）低于幼葉（3 377條），半甲基化率（27.83%）和總甲基化率（51.13%）高于幼葉（21.35%，45.90%）。成熟葉和幼葉DNA甲基化表現(xiàn)出的發(fā)育特異性，與前人在其他植物的研究結(jié)果一致。因此，在錐栗甲基化修飾研究過(guò)程中，有效區(qū)分不同生長(zhǎng)發(fā)育階段的樣品，是制定合理有效采樣方案的前提。

參考文獻(xiàn)：

ANGERS B， CASTONGUAY E， MASSICOTTE R， 2010. Environmentally induced phenotypes and DNA methylation： how to deal with unpredictable conditions until the next generation and after [J]. Mol Ecol， 19（7）： 1283-1295.

BOYKO A， KOVALCHUK I， 2008. Epigenetic control of plant stress response [J]. Environ Mol Mutagen， 49（1）：61-72.

CALARCO JP， BORGES F， DONOGHUE MTA， et al， 2012. Reprogramming of DNA methylation in pollen guides epigenetic inheritance via small RNA [J]. Cell， 151（1）： 194-205.

CERVERA MT， RUIZGARCIA L， MARTINEZZAPATER J， 2002. Analysis of DNA methylation in Arabidopsis thaliana based on methylationsensitive AFLP markers [J]. Mol Genet Genom， 268（4）：543-552.

CHAN SWL， HENDERSON IR， JACOBSEN SE， 2005. Gardening the genome： DNA methylation in Arabidopsis thaliana [J]. Nat Rev Genet， 6（5）： 351-360.

CHEN KS， LI F， XU CJ， et al， 2004. An efficient macromethod of genomic DNA isolation from Actinidia chinensis leaves [J]. Hereditas， 26 （4）：529-531. [陳昆松，李方，徐昌杰，等，2004. 改良CTAB法用于多年生植物組織基因組DNA的大量提取 [J]. 遺傳，26（4）：529-531.]

CHEN QF， ZHOU ZY， LI JL， et al， 2014. Automatic processing of the initial data of methylation level detected by fluorescently labeled MSAP [J]. Lett Biotechnol， 25（4）：527-531. [陳秋博，鄒智元，李金龍，等，2014. MSAP熒光檢測(cè)初始數(shù)據(jù)的自動(dòng)分析研究 [J]. 生物技術(shù)通訊，25（4）：527-531.]

CHU HJ， ZHANG JJ， WU Y， et al， 2011. Optimization and application of MSAP molecular marker reaction system of Salvia miltiorrhiza by orthogonal design [J]. Acta Bot BorealOccident Sin， 31（5）：920-928. [褚會(huì)娟，張今今，王喆，等，2011. 丹參MSAP分子標(biāo)記技術(shù)體系的優(yōu)化及其應(yīng)用 [J]. 西北植物學(xué)報(bào)，31（5）：920-928.]

FULNECˇEK J， KOVARˇK A， 2014. How to interpret methylation sensitive amplified polymorphism （MSAP） profiles？ [J]. BMC Genet， 15（1）： 1.

GAO BJ， DU J， GAO SH， et al， 2010. Genetic diversity and differentiations of fall webworm （Hyphantria cunea） populations [J]. Sci Silv Sin， 46（8）：120-124. [高寶嘉，杜娟，高素紅，等，2010. 美國(guó)白蛾種群的遺傳多樣性與遺傳分化 [J]. 林業(yè)科學(xué)，46（8）：120-124.]

GAO LX， CHEN JK， YANG J， 2008. Phenotypic plasticity： EcoDevo and evolution [J]. J Syst Evol， 46 （4）：441-451. [高樂(lè)旋，陳家寬，楊繼，表型可塑性變異的生態(tài)-發(fā)育機(jī)制及其進(jìn)化意義 [J]. 植物分類(lèi)學(xué)報(bào) 46（4）：441-451.]

GOLL MG， BESTOR TH， 2005. Eukaryotic cytosine methyltransferases [J]. Annu Rev Biochem， 74： 481-514.

GRANTDOWNTON RT， DICKINSON H， 2006. Epigenetics and its implications for plant biology 2. The ‘epigenetic epiphany： epigenetics， evolution and beyond [J]. Ann Bot， 97（1）：11-27.

GUO DJ， HE DW， 2000. The rank sum test for two independent samples and for paired datathe simplified formula of rank sum test [J]. Chin J Health Stat， 17（2）：83-85. [郭東星，何大衛(wèi)，2000. 成組設(shè)計(jì)的兩樣本比較秩和檢驗(yàn)與配對(duì)設(shè)計(jì)差值的符合秩和檢驗(yàn) [J]. 中國(guó)衛(wèi)生統(tǒng)計(jì)，17（2）：83-85.]

GUO GP， 2011. Research on the features of DNA methylation in the process of growth and development in Bamboo [D]. Beijing： Chinese Academy of Forestry. [郭廣平，2011. 竹類(lèi)植物生長(zhǎng)發(fā)育過(guò)程中的 DNA 甲基化研究 [D]. 北京：中國(guó)林業(yè)科學(xué)研究院.]

HASHIDA S， KITAMURA K， MIKAMI T， et al， 2003. Temperature shift coordinately changes the activity and the methylation state of transposon Tam 3 in Antirrhinum majus [J]. Plant Physiol， 132（3）： 1207-1216.

HEARD E， MARTIENSSEN RA， 2014. Transgenerational epigenetic inheritance： myths and mechanisms [J]. Cell， 157（1）： 95-109.

HERRERA CM， BAZAGA P， 2013. Epigenetic correlates of plant phenotypic plasticity： DNA methylation differs between prickly and nonprickly leaves in heterophyllous Ilex aquifolium （Aquifoliaceae） trees [J]. Bot J Linn Soc， 171（3）： 441-452.

JOHANNES F， PORCHER E， TEIXEIRA FK， et al， 2009. Assessing the impact of transgenerational epigenetic variation on complex traits [J]. PLoS Genet， 5（6）： e1000530.

KAEPPLER HF， PEDERSEN JF， 1997. Evaluation of 41 elite and exotic inbred Sorghum genotypes for high quality callus production [J]. Plant Cell， Tissue Organ Cult， 48（1）： 71-75.

KLIRONOMOS FD， BERG J， COLLINS S， 2013. How epigenetic mutations can affect genetic evolution： model and mechanism [J]. Bioessays， 35（6）： 571-578.

KLOSE RJ， BIRD AP， 2006. Genomic DNA methylation： the mark and its mediators [J]. Trends Biochem Sci， 31（2）： 89-97.

KOVAR A， KOUKALOVA B， BEZDE M， et al， 1997. Hypermethylation of tobacco heterochromatic loci in response to osmotic stress [J]. Theor Appl Genet， 95（1-2）： 301-306.

LI JL， TANG SQ， ZHOU ZY， et al， 2014. Establishment of fluorescence labeling and capillary electrophoresis in MSAP for Beijing You chicken [J]. Hereditas， 36（5）：495-502. [李金龍，唐韶青，鄒智元，等，2014. 北京油雞MSAP毛細(xì)管電泳熒光檢測(cè)技術(shù)的建立. 遺傳，36（5）：495-502.]

LUKENS LN， ZHAN S， 2007. The plant genomes methylation status and response to stress： implications for plant improvement [J]. Curr Opin Plant Biol， 10（3）： 317-322.

MA KF， ZHANG ZY， WANG SQ， et al， 2012. Optimization of MSAP system and preliminary analysis of DNA methylation in Populus tomentosa [J]. J NE For Univ， 40（12）：1-7. [馬開(kāi)峰，張志毅，王斯琪，等，2012. 毛白楊MSAP體系優(yōu)化及DNA甲基化的初步分析 [J]. 東北林業(yè)大學(xué)學(xué)報(bào)，40（12）：1-7.]

MORGAN HD， SANTOS F， GREEN K， et al， 2005. Epigenetic reprogramming in mammals [J]. Human Mol Genet， 14（suppl 1）： R47-R58.

PASZKOWSKI J， GROSSNIKLAUS U， 2011. Selected aspects of transgenerational epigenetic inheritance and resetting in plants [J]. Curr Opin Plant Biol， 14（2）： 195-203.

PAUN O， BATEMAN RM， FAY MF， et al， 2010. Stable epigenetic effects impact adaptation in allopolyploid orchids （Dactylorhiza： Orchidaceae） [J]. Mol Biol Evol， 27（11）： 2465-2473.

PENG H， ZHANG J， 2009. Plant genomic DNA methylation in response to stresses： Potential applications and challenges in plant breeding [J]. Prog Nat Sci， 19（9）： 1037-1045.

QIAN XW， YUAN YL， RONG P， et al， 2014. Analysis of methylationsensitive amplified polymorphism in cotton genome under salt stress based on capillary electrophoresis [J]. Chin J Biochem Mol Biol， 30（3）：298-306. [錢(qián)曉偉，袁有祿，榮平，等，2014. 鹽脅迫下棉花基因組基于毛細(xì)管電泳的MSAP分析 [J]. 中國(guó)生物化學(xué)與分子生物學(xué)報(bào)，30（3）：298-306.]

QIN GL， TAN JY， 2014. Mantel test and its appications in exology [J]. J Biomathem， 29 （03）：507-512. [覃光蓮， 譚勁英， 2014.