王靜宇 陳曉慧 賴鐘雄

摘? 要：表觀遺傳調(diào)控是指基于非DNA序列的改變所致基因功能的變化，最終導(dǎo)致生物多效性表型。表觀遺傳調(diào)控因子通過多種途徑調(diào)控基因表達(dá)，參與植物幾乎所有的生長發(fā)育過程。本文介紹了DNA甲基化、組蛋白修飾、RNA甲基化、染色質(zhì)重塑、非編碼RNA、基因印記和母性效應(yīng)等植物表觀遺傳修飾主要類型的分子機(jī)制，以及它們?cè)诩?xì)胞增殖、細(xì)胞分化、生物與非生物脅迫、開花、果實(shí)成熟和胚胎發(fā)育中的作用機(jī)制和生物學(xué)功能，并對(duì)其在農(nóng)業(yè)上的應(yīng)用作了展望。

關(guān)鍵詞：表觀遺傳;植物發(fā)育;DNA甲基化;組蛋白修飾;染色質(zhì)重塑

中圖分類號(hào)：Q943.2? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

Abstract： Epigenetic regulation refers to the functions of genes based on the change of non-DNA sequence， which finally leads to the biological pleiotropic phenotype. Epigenetic regulators are involved in almost all plant growth and development by directing gene expression at multiple levels. In this paper， the main types and molecular mechanisms of plant epigenetic modifications， such as DNA methylation， histone modification， RNA methylation， chromatin remodeling， non-coding RNA， gene imprinting and maternal effects， as well as their action mechanisms and biological functions in cell proliferation， cell differentiation， biotic and abiotic stresses， flowering， fruit maturation and embryo development were reviewed while their applications in agriculture were finally prospected.

Keywords： epigenetics; plant development; DNA methylation; histone modification; chromatin remodeling

DOI： 10.3969/j.issn.1000-2561.2020.10.016

在1957年，由Waddington提出表觀遺傳學(xué)（epigenetics）這一概念之前，經(jīng)典遺傳學(xué)一直支配著人們對(duì)生物遺傳理論的認(rèn)知。隨著對(duì)生命科學(xué)的深入研究，出現(xiàn)了一些無法用經(jīng)典遺傳學(xué)解釋的生命現(xiàn)象，例如，基因組完全一致的同卵雙生即使成長于相同的環(huán)境，他們的性格、身體素質(zhì)等也會(huì)存在明顯的不同;某些表型只由其中一個(gè)親本的基因決定，另一親本的基因卻保持沉默;馬、驢正反交的后代會(huì)有較大的差異等。如今，作為遺傳學(xué)中的前沿分支學(xué)科，表觀遺傳學(xué)的發(fā)展則為這些現(xiàn)象的理解提供了可能。表觀遺傳調(diào)控是指基于非DNA序列的改變所致基因功能的變化，這些穩(wěn)定的、可遺傳的變化最終導(dǎo)致不同的表型。它與孟德爾的核內(nèi)遺傳規(guī)律并不相符，由此認(rèn)為基因組中可能主要蘊(yùn)含2種類型的遺傳信息：一類是由DNA序列提供的經(jīng)典遺傳信息;另一類則是指導(dǎo)基因表達(dá)模式分化的表觀遺傳信息。如今隨著分子生物技術(shù)的迅速發(fā)展和完善，表觀遺傳調(diào)控在人類醫(yī)藥學(xué)、微生物代謝和植物生長發(fā)育過程中的研究和應(yīng)用取得了相當(dāng)?shù)倪M(jìn)展[1]。近幾年，表觀遺傳修飾在擬南芥（Arabidopsis thaliana）、水稻（Oryza sativa）等模式植物中的開花、響應(yīng)逆境脅迫、衰老等重要發(fā)育過程中的調(diào)控機(jī)制有了較為系統(tǒng)的認(rèn)知。本文就植物表觀遺傳修飾類型及其在植物重要生物學(xué)過程中的分子機(jī)制進(jìn)行系統(tǒng)表述。

1? 植物表觀遺傳調(diào)控類型及其分子機(jī)制

植物表觀遺傳主要涉及DNA甲基化、組蛋白修飾、RNA甲基化、染色質(zhì)重塑、非編碼RNA、基因組印跡和母性效應(yīng)等典型的調(diào)控類型。這些修飾類型彼此關(guān)聯(lián)，通過指導(dǎo)基因的時(shí)空表達(dá)參與生物學(xué)表型的重塑，對(duì)維持基因組穩(wěn)定和植物正常的生長發(fā)育至關(guān)重要。

1.1? DNA甲基化

DNA甲基化可以通過改變?nèi)旧|(zhì)結(jié)構(gòu)、DNA構(gòu)象、組蛋白修飾以及DNA與蛋白質(zhì)的相互作用方式調(diào)控基因的表達(dá);也可以通過抑制轉(zhuǎn)座子和外源DNA的轉(zhuǎn)錄，維持基因組的穩(wěn)定性，并在植物響應(yīng)環(huán)境脅迫、花期調(diào)控、體細(xì)胞無性系發(fā)育等過程中執(zhí)行重要功能。真核生物中的DNA甲基化主要以5-mC和N6-mA形式存在，其中5-mC主要圍繞增強(qiáng)子和啟動(dòng)子位點(diǎn)富集，是基因轉(zhuǎn)錄受抑制的表現(xiàn)，N6-mA則在轉(zhuǎn)錄起始位點(diǎn)有較高的豐度，是基因轉(zhuǎn)錄活性較高的標(biāo)志[2]。

目前研究顯示植物DNA甲基化主要有2種模式：存在于對(duì)稱CG和CHG序列中的維持甲基化和非對(duì)稱CHH（H指A、T和G）序列中的從頭甲基化。植物中的從頭甲基化被稱為RdDM途徑，主要由一系列RNA（包括siRNA、scaffold RNA）和蛋白質(zhì)介導(dǎo)（圖1）。在經(jīng)典RdDM途徑中，首先轉(zhuǎn)座子和高度重復(fù)序列區(qū)域在RNA聚合酶Ⅳ的作用下合成ssRNA（single-stranded RNA），隨后RNA聚合酶RDR2以之為模板合成dsRNA（double-stranded RNA），接著由DCL3（dicer-like 3）對(duì)其切割產(chǎn)生siRNA（short interfering RNA），經(jīng)HEN1（hua-enhancer 1）加工修飾后的siRNA招募DRM2（domains rearranged methylase 2），并與scaffold RNA和AGO蛋白（argonaute protein）結(jié)合配對(duì)甲基化同源DNA序列[3]。近期研究表明，與RNA聚合酶Ⅳ同源的RNA聚合酶Ⅱ也可誘導(dǎo)非經(jīng)典RdDM途徑。在此途徑中，RNA聚合酶Ⅱ不僅可以轉(zhuǎn)錄合成siRNA和scaffold RNA，而且可以通過招募RNA聚合酶Ⅳ、Ⅴ在RdDM靶標(biāo)區(qū)域合成siRNA。CG和CHG維持甲基化分別由甲基轉(zhuǎn)移酶MET1和CMT3介導(dǎo)。CMT3和H3K9me2組蛋白甲基轉(zhuǎn)移酶SUVH4的缺失會(huì)降低DNA甲基化水平，但它們?cè)贑HG甲基化中的互作機(jī)制仍有待了解[4]。同時(shí)，DNA甲基化可由DNA糖基化酶/裂解酶ROS1介導(dǎo)的去甲基化途徑逆轉(zhuǎn)，存在于基因組中的主動(dòng)或被動(dòng)去甲基化可激活沉默狀態(tài)的基因。被動(dòng)去甲基化是指基因組中DNA甲基轉(zhuǎn)移酶濃度或活性降低時(shí)導(dǎo)致未發(fā)生甲基化的胞嘧啶取代已甲基化的胞嘧啶或缺乏甲基供體致使DNA甲基化水平下降。主動(dòng)去甲基化是ROS1基因家族切除5-mC后由堿基切除修復(fù)觸發(fā)的DNA去甲基化機(jī)制。DNA糖基化酶還包DME（demeter）、DML2（demeter- like protein 2）和DML3。

1.3.2? m5C? 研究表明，mRNA和tRNA m5C修飾位點(diǎn)分布于真核生物和古細(xì)菌中，但不存在于真細(xì)菌中，mRNA m5C水平在擬南芥不同組織細(xì)胞中呈現(xiàn)空間變化，提示m5C修飾的物種保守性和組織特異性[14]。m5C甲基化主要由NSUN（NOL1/NOP2/sun）家族蛋白介導(dǎo)，該家族中的NSUN2（也稱為TRM4）成員為mRNA m5C特異的甲基轉(zhuǎn)移酶，主要依賴2個(gè)半胱氨酸位點(diǎn)發(fā)揮催化活性，它的缺失可導(dǎo)致m5C甲基化水平明顯下降。近期，從擬南芥中鑒定出的TRM4B對(duì)植株根的發(fā)育及響應(yīng)氧化應(yīng)激等過程具有重要作用，且TRM4B介導(dǎo)的m5C修飾可能對(duì)隨后的翻譯過程起負(fù)面調(diào)控作用[14]。此外，DNMT2（DNA methyltrans-ferase 2）是RNA m5C和DNA甲基化的通用甲基轉(zhuǎn)移酶，人體DNMT2介導(dǎo)tRNAAsp的C38和miRNA甲基化，但不介導(dǎo)DNA甲基化[15]。m5C位點(diǎn)發(fā)生甲基化修飾后，可與ALYREF特異結(jié)合調(diào)控mRNA的出核效率。早期研究發(fā)現(xiàn)RNA m5C修飾可能通過DNA甲基轉(zhuǎn)移酶TET去除，但近期實(shí)驗(yàn)顯示TET對(duì)RNA中m5C去甲基化起始位點(diǎn)5-mrC的酶解速率遠(yuǎn)低于對(duì)應(yīng)DNA鏈中的5-mdC，且目前尚未完全解析TET-RNA復(fù)合物的組分，因此TET介導(dǎo)m5C去甲基化機(jī)制仍有待明確[16]。

1.3.3? m1A? 目前m1A已在不同tRNA的第9、14、22、57和58核苷酸位點(diǎn)鑒定。其中m1A14和m1A22的修飾高度保守，僅存在于哺乳動(dòng)物中的細(xì)胞質(zhì)（cyt）tRNAPhe和細(xì)菌tRNA中，m1A57則僅作為水解脫氨基生成的1-甲基肌苷中間體短暫存在，因此，tRNA的第9和第58位為m1A修飾最佳研究位點(diǎn)。在真核生物中，（cyt）tRNA m1A58特異甲基轉(zhuǎn)移酶是由Trm6和Trm61這2個(gè)亞基組成的異四聚體，其中，催化亞基Trm61與輔因子SAM的活性甲基結(jié)合，Trm6由于缺少有關(guān)輔因子調(diào)節(jié)的保守基序而無法與之結(jié)合。每個(gè)異四聚體將2個(gè)tRNA分子結(jié)合到蛋白質(zhì)復(fù)合物2個(gè)遠(yuǎn)端的‘L形表面上，在此過程中tRNA中的D環(huán)和T環(huán)分離，包含待修飾堿基A58的T環(huán)通過C56堿基和糖磷酸骨架形成的大量氫鍵與活性位點(diǎn)穩(wěn)定結(jié)合;真核m1A9甲基化由Trm10介導(dǎo)，但Trm10與HSD10（17β-hydroxysteroid dehydrogenase type 10）復(fù)合才具有甲基轉(zhuǎn)移酶活性，Trm10-HSD10復(fù)合物以SAM為甲基供體催化m1A9位點(diǎn)的甲基化[17]。研究表明，m1A9和m1A58均可通過誘導(dǎo)起始tRNA的正確折疊提高其穩(wěn)定性，其中m1A58去甲基化酶ALKBH1已被鑒定，它的缺失會(huì)導(dǎo)致起始tRNA的積累，這與m1A58在起始tRNA結(jié)構(gòu)穩(wěn)定性的作用一致[18]。此外，m1A富集于mRNA的5'-UTR區(qū)，并影響其定位、剪切、翻譯和穩(wěn)定性等，且葡萄糖和氨基酸的含量與mRNA m1A甲基化水平正相關(guān)。目前已將ALKB家族中的ALKBH3作為逆轉(zhuǎn)mRNA m1A甲基化并參與其動(dòng)態(tài)調(diào)控的候選分子。

1.4? 染色質(zhì)重塑

組蛋白共價(jià)修飾和依賴ATP的染色質(zhì)重塑復(fù)合物CRCs（ATP-dependent chromatin remodeling complex）可以通過改變核小體的組成和定位局部或整體重塑染色質(zhì)結(jié)構(gòu)，從而改變DNA的“可及性”以調(diào)節(jié)組蛋白和DNA的互作。CRCs依賴其核心亞基SWI2/SNF2 ATPase中共有的DEXDc和HELICc結(jié)構(gòu)域水解ATP供能，在調(diào)控基因轉(zhuǎn)錄、DNA修復(fù)和轉(zhuǎn)錄重組等核過程發(fā)揮重要作用?；贏TPase結(jié)構(gòu)的差異，CRCs被分為SWI/SNF、ISWI、CHD和INO80/SWR1這4個(gè)亞家族，它們具有各自特異的結(jié)構(gòu)域和識(shí)別作用（圖2）。

1.4.1? SWI/SNF? SWI/SNF復(fù)合體是植物中最重要的一類保守CRCs，大多數(shù)SWI/SNF復(fù)合體可通過改變DNA與組蛋白的結(jié)合狀態(tài)體外滑動(dòng)或分解核小體引起染色質(zhì)重塑。由1個(gè)SWI2/SNF2 ATPase、1個(gè)SNF5和2個(gè)SWI3亞基構(gòu)成的SWI/ SNF核心亞復(fù)合體即具有基本的染色質(zhì)重塑活性，這些核心亞基可與一系列轉(zhuǎn)錄因子和核蛋白互作調(diào)控下游基因的表達(dá)，在植物開花、葉片發(fā)育、胚胎發(fā)生和應(yīng)激反應(yīng)等過程發(fā)揮作用。擬南芥基因組編碼4個(gè)SWI2/SNF2（BRM、SYD、CHR12/MINU1、CHR 23/MINU 2）、1個(gè)SNF5（BSH）和4個(gè)SWI3家族成員（SWI3A， B， C， D）[20]。其中只有BRM（brahma）具有SWI/SNF-ATPase典型的C端Bromo結(jié)構(gòu)域，并被作為組蛋白乙?；慕Y(jié)合蛋白，Bromo結(jié)構(gòu)域通過與C末端的DNA結(jié)合區(qū)域構(gòu)成核小體結(jié)合模塊介導(dǎo)BRM和其他蛋白的互作，由此推測(cè)BRM最有可能是植物SWI/SNF復(fù)合體的催化亞基。擬南芥BRM是一個(gè)組成龐大、功能多樣的蛋白質(zhì)復(fù)合物，由SUMO連接酶MMS21對(duì)BRM進(jìn)行翻譯后修飾以維持其穩(wěn)定性[21]。在擬南芥基因組中，BRM與近1000個(gè)基因的啟動(dòng)子和終止子區(qū)域互作調(diào)節(jié)靶基因的轉(zhuǎn)錄，表明BRM介導(dǎo)的轉(zhuǎn)錄可能受到其他亞基、互作蛋白或特定染色質(zhì)區(qū)域中的不同核小體靶標(biāo)的影響。此外，核小體的穩(wěn)定性與組蛋白H2A和其變體H2A.Z的替換相關(guān)，研究表明，SWI/SNF和H2A.Z定位于擬南芥基因組中數(shù)千個(gè)位點(diǎn)，它們?cè)谶@些位點(diǎn)中既相互作用又彼此拮抗，共同調(diào)控相關(guān)核小體動(dòng)力學(xué)[22]。

1.4.2? ISWI? 在真核生物中，ISWI蛋白通過與DDT結(jié)構(gòu)域蛋白（DNA-binding homeobox and different transcription factor）形成保守CRCs改變核小體的位置;在此過程中，ISWI-ATPase C端的SANT、SLIDE結(jié)構(gòu)域通過與核小體連接DNA互作將ISWI蛋白錨定在核小體上，并促進(jìn)核小體延DNA鏈滑動(dòng)。擬南芥基因組編碼的2種ISWI蛋白CHR11和CHR17負(fù)責(zé)調(diào)控核小體在基因區(qū)的分布“間距”，并通過與多種DDT結(jié)構(gòu)域蛋白互作影響植物營養(yǎng)器官和花器官的發(fā)育[23]。

1.4.3? CHD? 植物CHD蛋白分為CHD1和CHD3兩類。CHD1蛋白的C端具有DNA結(jié)合區(qū)，并且可以通過獨(dú)立或協(xié)作DNA復(fù)制后染色質(zhì)的組裝指導(dǎo)基因轉(zhuǎn)錄所需要的核小體定位;CHD3缺少DNA結(jié)合區(qū)，但含有與H3K4me3特異結(jié)合并促進(jìn)下游基因轉(zhuǎn)錄激活的PHD結(jié)構(gòu)域。擬南芥編碼1個(gè)CHD1同源蛋白CHR5和3個(gè)與CHD3同源的PKL（pickle）、CHR4、CHR7蛋白。其中，CHR5是SNC1（suppressor of npr1， constitutive 1）在植物免疫反應(yīng)中誘導(dǎo)表達(dá)的正調(diào)控因子，PKL則通過調(diào)節(jié)H3K27me3的表達(dá)水平參與植物的生長發(fā)育[24-25]。

1.4.4? INO80/SWR1? 與其他亞家族相比，INO80/SWR1家族重塑因子的ATPase結(jié)構(gòu)存在一個(gè)相當(dāng)長的“插入?yún)^(qū)”，用于招募Rvb1/2解旋酶。該家族包含重塑因子INO80和SWR1，INO80復(fù)合體催化核小體H2A-H2B取代H2A.Z-H2B，SWR1則反向催化二聚體的替換。在真核生物中兩者高度保守，并共享RVB1和RVB2等多個(gè)不可或缺的亞基。但也具有各自特異的組成單位，其中ARP5和ARP6分別為INO80和SWR1的高

度特異性亞基。擬南芥INO80同源蛋白AtINO80與ARP5互作，對(duì)復(fù)制脅迫條件下染色體穩(wěn)定性的維持具有重要作用[26]。然而，AtINO80向局部染色質(zhì)的募集會(huì)受到ARP6介導(dǎo)的H2A.Z定位的影響。ARP6作為SWR1的特異亞基，最近被證實(shí)與ISWI的作用蛋白MBD9（methyl-CpG- binding domain 9）存在組合互作[27];另外，ARP6和ARP5可聯(lián)合SWR1的其它核心亞基包括PIE1、SWC6參與調(diào)節(jié)DNA DSB（double-strand breakage）修復(fù)途徑同源重組的頻率?？傊?，INO80/SWR1復(fù)合體在DSB修復(fù)、染色體分離等多種核過程中存在亞基互作和組蛋白串?dāng)_，然而對(duì)其部分亞基包括RVB1、RVB2A/2B和精確的重塑機(jī)制如復(fù)合體被招募到其靶位點(diǎn)的途徑仍是未知。

1.5? 非編碼RNA

非編碼RNA（non-coding RNAs， ncRNAs）是一類由生物體基因組轉(zhuǎn)錄產(chǎn)生的含有遺傳物質(zhì)但不編碼蛋白質(zhì)的分子。全基因組和轉(zhuǎn)錄組測(cè)序結(jié)果表明，超過50%的擬南芥基因組可轉(zhuǎn)錄成RNA，其中絕大多數(shù)是缺乏蛋白編碼能力但在真核生物的生長發(fā)育和應(yīng)激反應(yīng)中發(fā)揮重要作用的ncRNAs。根據(jù)ncRNAs的功能差異，將其分為組成型ncRNAs（constitutive ncRNAs）和調(diào)控型ncRNAs（regulatory ncRNAs）兩大類，前者包括tRNA、rRNA、snoRNA（small nucleolar RNA）、snRNA（small nuclear RNA），后者主要包括miRNA（microRNA）、lncRNA（long ncRNA）和siRNA等。其中，調(diào)控型ncRNAs中研究最為透徹的miRNA是一類長為20～24 nt且高度保守的小分子RNA，通過與靶位點(diǎn)的堿基互補(bǔ)配對(duì)實(shí)現(xiàn)轉(zhuǎn)錄本的切割、降解和翻譯抑制，從而參與調(diào)節(jié)植物的生長發(fā)育。與miRNA不同，lncRNA是一類序列長度超過200 nt、作用機(jī)制復(fù)雜多樣的低保守性轉(zhuǎn)錄物，主要與大分子DNA、RNA和蛋白質(zhì)等組合互作，在轉(zhuǎn)錄、轉(zhuǎn)錄后和表觀遺傳水平調(diào)控RNA代謝、蛋白質(zhì)修飾、染色質(zhì)重塑等，同時(shí)lncRNA的轉(zhuǎn)錄、剪切也受到DNA甲基化、RNA和組蛋白修飾的調(diào)節(jié)。研究表明，lncRNA主要通過參與信號(hào)通路、募集RNA結(jié)合蛋白、招募核糖核蛋白或染色質(zhì)修飾酶至特定位點(diǎn)以及作為核糖核蛋白復(fù)合物的組合支架靶向蛋白質(zhì)編碼基因[28]。

1.6? 基因印記

基因印記（genic imprinting）是指在配子或合子發(fā)生期間，來自不同親本的等位基因由于受到DNA甲基化和FISPcG介導(dǎo)的組蛋白差異修飾而表現(xiàn)偏好親源的選擇性表達(dá)。其中母源特異表達(dá)而父源被印記或沉默的基因稱為母本表達(dá)的印記基因（mat ernally expressed imprinted gene， MEG），反之稱為父本表達(dá)的印記基因（paternally expressed imprinted gene， PEG）。近幾年相繼在擬南芥、高粱（Sorghum bicolor）和水稻中發(fā)現(xiàn)多個(gè)印記基因，其中只有4個(gè)被鑒定為PEG（PHE1、PEG1、HDG3和At5G62110），而有20個(gè)表現(xiàn)為MEG印記。這些印記基因絕大多數(shù)在胚乳中表達(dá)，它們可通過調(diào)控早期胚乳的細(xì)胞化決定種子發(fā)育大小，并在種子花青素的合成以及營養(yǎng)物質(zhì)的運(yùn)輸與分配中發(fā)揮重要作用。

1.6.1? MEG印記機(jī)制? 植物基因組中由差異DNA甲基化形成MEG等位基因表達(dá)和沉默的機(jī)制主要與DNA糖基化酶DME的特異表達(dá)相關(guān)。DME具有5-mC切除活性，并在雌配子體的中央細(xì)胞中特異表達(dá)，同時(shí)DNA甲基轉(zhuǎn)移酶MET1的表達(dá)受到成視網(wǎng)膜細(xì)胞瘤途徑（retinoblastoma pathway）的抑制，此途徑涉及RBR1（retinobl astoma related1）及其伴侶分子MS1（multicopy supressor of ira1）的作用;但MET1的單獨(dú)抑制并不足以移除MEG母源等位基因的甲基化印記，因此，MEG在卵細(xì)胞中缺乏DME和MET1被抑制的情況下仍維持甲基化狀態(tài)。雄配子體中的MET1雖正常表達(dá)，但不存在DME的活性轉(zhuǎn)錄，因此最終只有MEG母源等位基因在胚乳中表達(dá)，而在胚胎中保持沉默[29]。

1.6.2? PEG 印記機(jī)制? 目前有關(guān)PEG的印記機(jī)制主要在印記基因PHE1的研究中較為完善，PHE1等位基因的表達(dá)與多梳抑制復(fù)合體FISPcG2相關(guān)。在雌配子體的中央細(xì)胞中，DME的特異表達(dá)導(dǎo)致PHE1 3'端差異甲基化區(qū)（diffe r entially methylated region， DMR）的去甲基化，隨后FISPcG2復(fù)合體與PHE1啟動(dòng)子區(qū)的多梳響應(yīng)元件（polycomb response element， PRE）和去甲基化的DMR形成穩(wěn)定的染色質(zhì)內(nèi)環(huán)，同時(shí)FISPcG2復(fù)合體通過催化H3K27me3的形成抑制PHE1的表達(dá);而在雄配子體中由于DMR的甲基化修飾和FISPcG復(fù)合體的缺失，PHE1的啟動(dòng)子區(qū)域不含或含有少量的甲基化位點(diǎn)，因此，PHE1的父源等位基因在胚乳中活躍表達(dá)[30]。

1.7? 母性效應(yīng)

母性效應(yīng)是指在不改變遺傳因素的條件下，由母本核基因影響或決定其子代某一性狀的表現(xiàn)型，以緩沖子代在異質(zhì)環(huán)境中的生存壓力增強(qiáng)其適應(yīng)性。早期胚胎的發(fā)育即涉及典型的母性效應(yīng)。自受精開始，特異儲(chǔ)存在卵母細(xì)胞中的一些基因啟動(dòng)表達(dá)，其轉(zhuǎn)錄產(chǎn)物mRNA或編碼蛋白質(zhì)在卵母細(xì)胞向胚胎的過渡中發(fā)揮著關(guān)鍵的作用，而這些基因（母性效應(yīng)基因）一旦有所缺失，即使相應(yīng)的父源基因正常表達(dá)也會(huì)導(dǎo)致胚胎發(fā)育缺陷或停滯。目前，植物中有關(guān)母性效應(yīng)基因的研究主要集中于擬南芥，擬南芥母性效應(yīng)胚胎滯育基因MEE/MIS1是FISPcG2的組成亞基，負(fù)責(zé)調(diào)控印記基因的表達(dá)與沉默[31];MEE/MIS1在擬南芥響應(yīng)干旱脅迫、適時(shí)開花和胚珠發(fā)育等過程中發(fā)揮著重要的作用。我們從龍眼中鑒定出的擬南芥MEE/MIS1同源基因MEE70-1a/b主要定位于過氧化物酶體和細(xì)胞核，生物信息學(xué)分析顯示它們可能協(xié)同調(diào)控種子萌發(fā)和胚胎發(fā)育等過程，與擬南芥相比可能具有更復(fù)雜的代謝途徑[32];但迄今為止，母性效應(yīng)在植物中的研究較為落后，其在早期胚胎發(fā)育中的詳細(xì)調(diào)控機(jī)制仍需進(jìn)一步研究。

2? 植物表觀遺傳的生物學(xué)功能

2.1? 細(xì)胞增殖

表觀遺傳修飾廣泛參與了細(xì)胞對(duì)基因表達(dá)的控制，在細(xì)胞生長、分化、增殖和疾病狀態(tài)中起關(guān)鍵性的作用。Pathania等[33]研究發(fā)現(xiàn)DNA甲基轉(zhuǎn)移酶DNMT1對(duì)于乳腺干細(xì)胞和癌癥干細(xì)胞的維持是必不可少的，在乳腺發(fā)育和腫瘤發(fā)展中起到了重要的作用;EZH2（enhancer of zeste homolog 2）主要對(duì)組蛋白H3K27進(jìn)行甲基化，由此沉默下游基因，在細(xì)胞增殖、分化和腫瘤形成方面起到了重要的作用。Yin等[34]通過研究表明，EZH2通過組蛋白甲基轉(zhuǎn)移酶活性調(diào)控了NK（natural killer）細(xì)胞的分化和功能，去除EZH2或抑制其活性可以促使造血干細(xì)胞及祖細(xì)胞生成更多更強(qiáng)的NK細(xì)胞。此外，表觀遺傳調(diào)控異常會(huì)干擾干細(xì)胞功能和組織內(nèi)穩(wěn)態(tài)，其嚴(yán)重程度取決于出現(xiàn)問題的組織類型和表觀遺傳學(xué)調(diào)控子。

愈傷組織的形成是高等植物組織培養(yǎng)中器官從頭再生過程的一個(gè)關(guān)鍵步驟。近年的研究發(fā)現(xiàn)Polycomb和ISWI表觀遺傳通路對(duì)于葉片形成愈傷組織至關(guān)重要。Polycomb通過組蛋白甲基化途徑控制基因組表觀信息在細(xì)胞命運(yùn)轉(zhuǎn)變中的重排，主要負(fù)責(zé)沉默葉特征基因的表達(dá)。ISWI是ATP依賴型染色質(zhì)重塑因子，在擬南芥中由CHR11與CHR17兩個(gè)基因編碼，chr11-1、chr17-1雙突變體可導(dǎo)致其葉片外植體無法再生愈傷組織[35]。

植物器官的發(fā)育取決于精細(xì)的生長調(diào)控。組蛋白去乙?；窰DACs通過形成具有不同調(diào)節(jié)亞基的多種共抑制因子復(fù)合物來實(shí)現(xiàn)其在特異的組織或細(xì)胞中的功能，調(diào)節(jié)特異靶標(biāo)基因的表達(dá)。其中，HDC1（histone deacetylase complex 1）作為核心與HDAC復(fù)合物的多個(gè)組成蛋白相互作用，調(diào)節(jié)植物的器官大小。最近，Zhang等[36]報(bào)道了黃瓜（Cucumis sativus）果實(shí)細(xì)胞增殖的表觀遺傳調(diào)控機(jī)制，該研究鑒定了一個(gè)細(xì)胞增殖受到抑制的短瓜突變體sf2。SF2編碼了一個(gè)擬南芥HDC1的同源蛋白，主要在分生組織中表達(dá)，這與其通過HDAC復(fù)合體調(diào)節(jié)細(xì)胞增殖的作用相一致。在sf2突變體中，SF2第515位的甘氨酸（G）突變?yōu)楣劝彼幔‥），導(dǎo)致SF2在全基因組水平的結(jié)合能力減弱，組蛋白乙?；缴摺Ｈ蚪M水平比對(duì)分析表明，SF2直接靶向并促進(jìn)多種植物激素途徑和細(xì)胞周期調(diào)控關(guān)鍵基因的組蛋白去乙酰化，并在細(xì)胞增殖過程中抑制或激活它們的表達(dá)。進(jìn)一步研究表明，SF2通過直接靶向細(xì)胞分裂素和多胺的生物合成與代謝來控制果實(shí)細(xì)胞增殖。

2.2? 細(xì)胞分化

細(xì)胞分化是指多能性細(xì)胞特化為不同類型細(xì)胞的過程，這一過程依賴于表觀遺傳修飾對(duì)發(fā)育調(diào)控因子時(shí)空表達(dá)的有序調(diào)節(jié)。研究表明，大麥（Hordeum vulgare）小孢子重塑胚胎發(fā)育的過程中伴隨著細(xì)胞命運(yùn)和發(fā)育程序變化的DNA甲基化重編程，DNA甲基化在胚細(xì)胞第一次發(fā)生分裂的過程中維持低水平狀態(tài)，但隨著后期胚胎的分化逐漸升高[37];玉米（Zea mays）葉片組織細(xì)胞從增殖到分化的過程中伴隨著DNA甲基化水平的顯著變化，其中參與細(xì)胞周期、染色質(zhì)重塑等過程的基因發(fā)生了差異甲基化，但基因區(qū)內(nèi)的差異甲基化與轉(zhuǎn)錄水平無相關(guān)性，只有位于基因上游的差異甲基化與轉(zhuǎn)錄負(fù)相關(guān)。近期研究發(fā)現(xiàn)，組蛋白修飾復(fù)合物PcG-PRC2突變體中H3K 27 me3的缺失導(dǎo)致多個(gè)階段的發(fā)育異常，包括愈傷組織再生、胚胎到幼苗的過渡[38]。此外，PRC2組件CLF突變體中根分生組織調(diào)節(jié)因子WOX5（wuschel-related homeobox 5）、PLT（plethoral）和AGLs（agamous-like）被顯著上調(diào)，最終導(dǎo)致根干細(xì)胞活性增強(qiáng)和分生區(qū)擴(kuò)大，而染色質(zhì)重塑因子PKL突變體則表現(xiàn)出了相反的表型，提示PcG可能通過拮抗PKL維持根分生組織活性[39]。最近研究表明，PKL參與植物基因組H3K27me3的動(dòng)態(tài)調(diào)節(jié)也支持了這一推斷[25]。另外，組蛋白乙酰轉(zhuǎn)移酶GCN5和重塑子BRM可通過調(diào)控PLT的表達(dá)參與根尖干細(xì)胞微環(huán)境的調(diào)節(jié)，然而GCN5是否通過調(diào)節(jié)PLT的乙酰化水平影響其表達(dá)還有待證明。在水稻中，轉(zhuǎn)錄因子WOX11通過募集GCN5及其互作蛋白ADA2激活下游基因的轉(zhuǎn)錄，從而促進(jìn)冠根的伸長;而組蛋白去乙?；窰DT1/2可通過靶向赤霉素氧化酶GAox2的啟動(dòng)子區(qū)抑制其在分生區(qū)的表達(dá)，對(duì)調(diào)節(jié)根尖靜止中心的細(xì)胞分裂發(fā)揮了重要作用[40]。在莖尖分生組織中，AG（agmous）在WUS（wuschel）位點(diǎn)通過招募CLF、EMF2等PcG因子抑制其表達(dá)，最終在干細(xì)胞身份無法維持的情況下促進(jìn)花器官的形成[41]。在胚胎發(fā)生的過程中，KNOX（knotted1-like homeobox）負(fù)責(zé)維持周邊區(qū)內(nèi)部細(xì)胞的增殖和未分化狀態(tài)，而PcG組件LHP1、CLF、VRN2、FIE、RING1A/B和EMF4參與了KNOX的轉(zhuǎn)錄調(diào)節(jié)[42]。總之，PcG介導(dǎo)的轉(zhuǎn)錄抑制在平衡細(xì)胞增殖和分化中發(fā)揮了關(guān)鍵的作用。

2.3? 生物脅迫

植物通過病原體相關(guān)分子模式（PAMP）觸發(fā)的免疫（PTI）和效應(yīng)器觸發(fā)的免疫（ETI）保護(hù)自己免受病原體的攻擊。PTI和ETI途徑的激活依賴于多種調(diào)節(jié)機(jī)制嚴(yán)格調(diào)控的轉(zhuǎn)錄重編程。研究發(fā)現(xiàn)PolⅤ依賴型DNA甲基化的缺失增強(qiáng)了擬南芥對(duì)丁香假單胞菌番茄致病變種DC3000的免疫[43]。在小麥（Triticum aestivum）中，接種白粉菌（Blumeria graminis f. sp.， Bgt）通過降低AGO4a相關(guān)聯(lián)的siRNA水平下調(diào)CHH甲基化，從而激活A(yù)eGlu/PR2-like等防御基因以增強(qiáng)對(duì)Bgt的抗性，提示DNA甲基化可能作為細(xì)菌和活體營養(yǎng)性真菌抗性的負(fù)調(diào)控因子[44]，這與早期DNA甲基化抑制劑5-氮雜胞苷（5-azadeoxycytidine）處理增強(qiáng)水稻白葉枯病細(xì)菌抗性的實(shí)驗(yàn)結(jié)果一致。然而，RdDM通過促進(jìn)茉莉酸信號(hào)正調(diào)節(jié)植物對(duì)灰葡萄孢（Botrytis cinerea）和黃瓜枯萎菌（Plecto sph aerella cucumerina）等死體營養(yǎng)性真菌的抗性，提示植物的免疫應(yīng)答途徑中存在不同的DNA甲基化調(diào)節(jié)機(jī)制[45]。與此類似，染色質(zhì)重塑復(fù)合體SWR1和H2A.Z也可誘導(dǎo)擬南芥對(duì)活體和死體營養(yǎng)病原物不同的防御機(jī)制。在水稻中，SNF2類重塑因子BRHIS1通過與單泛素化組蛋白變體包括H2B.7和H2A.Xa/Xb/3互作促進(jìn)不依賴水楊酸信號(hào)防御基因的表達(dá)，從而誘導(dǎo)對(duì)稻瘟病菌的先天免疫[46]。植物免疫也需要SDG組蛋白甲基轉(zhuǎn)移酶依賴性轉(zhuǎn)錄組的2個(gè)基因CCR2和CER3。它們分別編碼類胡蘿卜素和角質(zhì)層臘質(zhì)生物合成酶。研究表明SDG8/25突變體中H2Bub1靶向的抗病基因SNC1以及CCR2、CER3位點(diǎn)的H2B泛素化水平均顯著降低，提示SDG8/25可能通過影響H3K4、H3K36甲基化水平直接或通過H2B泛素化間接調(diào)節(jié)免疫基因的表達(dá)[47]。

2.4? 非生物脅迫

2.4.1? 低溫脅迫? CBF（C repeat binding factor）- COR（cold responsive）是植物低溫響應(yīng)機(jī)制的核心組分。研究發(fā)現(xiàn)染色質(zhì)重塑因子PKL介導(dǎo)了CBF-COR低溫應(yīng)答途徑，PKL參與RdDM途徑的調(diào)節(jié)，并通過與組蛋白甲基轉(zhuǎn)移酶CLF和SWR1類染色質(zhì)重塑因子PIE1協(xié)作調(diào)控COR基因中的H3K27me3水平，因此，PKL可能通過H3K27me3的轉(zhuǎn)錄抑制效應(yīng)調(diào)節(jié)COR基因的染色質(zhì)狀態(tài)以響應(yīng)低溫脅迫[48]。此外，低溫脅迫轉(zhuǎn)錄因子的表達(dá)還受組蛋白乙酰化酶HAT和去乙?；窰DAC的動(dòng)態(tài)調(diào)節(jié)。正常條件下，擬南芥去乙酰化酶HD2C被活性轉(zhuǎn)錄的E3泛素連接酶HOS15募集至COR的啟動(dòng)子區(qū)，并通過去乙酰化作用抑制其表達(dá)。低溫處理下，HOS15介導(dǎo)的HD2C泛素化降解可提高COR的乙?；剑瑫r(shí)轉(zhuǎn)錄因子CBF被招募至下游COR的啟動(dòng)子區(qū)以提高其轉(zhuǎn)錄活性[49]。在血橙（Citrus sinensis）中，低溫脅迫可以特異誘導(dǎo)果實(shí)中反轉(zhuǎn)座子的轉(zhuǎn)錄激活，從而促進(jìn)下游Ruby基因的表達(dá)以影響隨后的花色苷合成[50];但在不同品種中花色苷合成所需的Ruby表達(dá)水平有較大差異，提示合理利用低溫脅迫可以形成理想的農(nóng)藝性狀。

2.4.2? 高溫脅迫? 表觀遺傳修飾在植物的高溫應(yīng)激反應(yīng)中主要具有3個(gè)方面的信號(hào)轉(zhuǎn)導(dǎo)途徑。其一，高溫應(yīng)答途徑中的核心轉(zhuǎn)錄因子HSFA1（heat shock transcription factor A1）負(fù)責(zé)激活下游HSFA2的表達(dá)，隨后HSFA2通過招募H3K9ac和H3K4me3促進(jìn)下游APX2、HSP18、HSP22等熱休克基因的轉(zhuǎn)錄[51];其二，高溫脅迫通過克服RdDM途徑對(duì)其靶標(biāo)SDC（suppressor of drm1 drm2 cmt3）的沉默效應(yīng)誘導(dǎo)下游熱休克基因的表達(dá)[52];其三，高溫脅迫通過誘導(dǎo)ONSEN（Ty1/ Copia-like）反轉(zhuǎn)座子的激活介導(dǎo)“跨代”脅迫記憶的形成，此過程中ONSEN的活性受RdDM途徑中siRNA的嚴(yán)格調(diào)控。早期研究發(fā)現(xiàn)HSPs（heat shock proteins）和HSFs的表達(dá)水平與水稻、大豆（Glycine max）、番茄（Solanum lycopersicum）等植物的耐熱能力呈正相關(guān)，如番茄HSP21和大豆HSFA1的過表達(dá)可分別增強(qiáng)它們?cè)诟邷孛{迫中的耐受性，但在無脅迫條件下也會(huì)對(duì)它們產(chǎn)生生長遲緩等負(fù)面作用[53]。

2.4.3? 鹽脅迫? 在植物中，HKT1（high-affinity k+ channel 1）與SOS（salt overly sensitive）途徑協(xié)作賦予植物在鹽脅迫下的耐受性。正常條件下，RdDM依賴型的DNA甲基化使擬南芥HKT1和鹽誘導(dǎo)轉(zhuǎn)錄因子MYB74維持低水平的轉(zhuǎn)錄，然而在鹽脅迫條件下隨著MYB74的上調(diào)表達(dá)，幾乎無法檢測(cè)到其甲基化和siRNA活性[54]。小麥HKT1也具有類似的DNA甲基化調(diào)節(jié)途徑，鹽脅迫特異誘導(dǎo)鹽敏感和耐鹽基因型品種芽和根中的m5C甲基化，從而抑制HKT2;1和HKT2;3的表達(dá)，然而根中HKT14的下調(diào)與m5C的甲基化作用無關(guān)[55]。此外，組蛋白修飾如甲基化、乙?；土姿峄矃⑴c植物的鹽脅迫應(yīng)答途徑。在正常處理下，2種組蛋白去乙?；窰D2C和HDA6協(xié)作抑制脫落酸響應(yīng)因子ABA1和ABI2的表達(dá)。其中，ABI2介導(dǎo)SOS2途徑的失活及其靶標(biāo)的去磷酸化，高鹽條件會(huì)抑制HD2C的活性，并誘導(dǎo)與染色質(zhì)密度相關(guān)的H3 Ser-10磷酸化，H3K9K14ac和H3K4me3等活性轉(zhuǎn)錄標(biāo)記被募集到鹽脅迫應(yīng)答因子中，而H3K9me2和H3K27me3等抑制轉(zhuǎn)錄標(biāo)記明顯下調(diào)。研究表明，高鹽處理下擬南芥HKT1的轉(zhuǎn)錄激活與其基因區(qū)H3K27me3的去除相關(guān)[56]。此外，H3、H4的乙?；虷3 Ser-10磷酸化水平在擬南芥和煙草（Nicotiana tabacum）的不同脅迫反應(yīng)中均明顯上調(diào)，但具有不同的調(diào)節(jié)機(jī)制。

2.4.4? 干旱脅迫? 干旱脅迫誘導(dǎo)ABA的生物合成，ABA則對(duì)植物的抗旱性起反饋調(diào)節(jié)作用。干旱脅迫可降低脅迫應(yīng)答因子RD29A和RD29B啟動(dòng)子區(qū)的核小體密度，該區(qū)域所含的ABRE（ABA-responsive element）可通過招募AREB（ABA-responsive element binding 1）和DREB（dehydration responsive element binding protein）促進(jìn)RD29A和RD29B的表達(dá)，且干旱持續(xù)時(shí)間與它們的表達(dá)量呈正相關(guān)。研究表明，SWI/SNF染色質(zhì)重塑因子的核心組分CHR12在干旱脅迫誘導(dǎo)的擬南芥瞬時(shí)生長停滯中發(fā)揮作用，另一組分SNF5在豌豆（Pisum sativum）的應(yīng)激反應(yīng)中誘導(dǎo)表達(dá)，提示SWI/SNF重塑復(fù)合物對(duì)不同植物干旱響應(yīng)途徑中染色體結(jié)構(gòu)的調(diào)節(jié)具有重要作用[57]。此外，干旱響應(yīng)基因的表達(dá)還受組蛋白修飾的動(dòng)態(tài)調(diào)控。ABA關(guān)鍵合成酶NCE3的活性與其基因區(qū)內(nèi)H3K4me3的水平相關(guān)，RD29A、RD29B、RD22和RAP2.4等干旱應(yīng)答因子的啟動(dòng)子區(qū)通過富集H3K4me3和H3K9ac實(shí)現(xiàn)活性轉(zhuǎn)錄，這些組蛋白活性標(biāo)記的豐度還隨著干旱程度的升高而升高。而在脅迫恢復(fù)的過程中，H3K9ac擁有比H3K4更快的降解速度[58]。最近研究發(fā)現(xiàn)，干旱脅迫應(yīng)答因子中H3K27me3水平的降低也是賦予植物抗旱性的因素之一，其結(jié)合蛋白PRC1復(fù)合物通過轉(zhuǎn)錄抑制干旱響應(yīng)的NAC轉(zhuǎn)錄因子ANAC019和ANAC055對(duì)ABA依賴的抗旱性起負(fù)調(diào)控作用。在毛果楊（Populus trichocarpa）中，多數(shù)NAC家族成員的啟動(dòng)子區(qū)含有ABRE元件，AREB可與之結(jié)合，并招募ADA2b-GCN5組蛋白乙酰轉(zhuǎn)移酶將H3K9ac富集于NAC中以促進(jìn)其表達(dá)，進(jìn)而增強(qiáng)植物的抗旱性[59]。對(duì)毛果楊的進(jìn)一步研究發(fā)現(xiàn)，干旱脅迫誘導(dǎo)的DNA甲基化指導(dǎo)相關(guān)基因的表達(dá)，但其調(diào)控機(jī)制仍處于未知狀態(tài)。已證實(shí)含有轉(zhuǎn)座元件MITE的玉米NAC111可通過H3K9me2和RdDM途徑抑制干旱脅迫應(yīng)答基因的表達(dá)[60]。

2.5? 開花

不同類型的花器官由各輪中同源異型A、B、C、D、E基因的特異組合所決定，大多數(shù)植物的花同源異型基因編碼MADS結(jié)構(gòu)域蛋白，其中A類基因編碼AP1（apetala 1）和AP2，B類編碼AP3和PI（pistillata），C類和E類則分別編碼AG和SEP（sepallata）轉(zhuǎn)錄因子。Smaczniak等[61]根據(jù)擬南芥花發(fā)育過程中的標(biāo)志性事件將其劃分為12個(gè)階段，在第1～2階段即花原基形成期間，花同源異型B和C類基因的表達(dá)被EMF2-PRC2復(fù)合物介導(dǎo)的H3K27me3沉默，SWI/SNF重塑因子SYD、BRM被成花關(guān)鍵基因LFY（leafy）、SEP3招募至AG和AP3位點(diǎn)進(jìn)行轉(zhuǎn)錄激活。與此同時(shí)，H3K4me3甲基轉(zhuǎn)移酶trxG蛋白ULT1通過與同源ATX1/SDG27互作激活A(yù)G和其他MADS家族基因的表達(dá)。從第3階段花萼原基出現(xiàn)開始，AG通過招募PcG抑制WUS，并通過驅(qū)除KNU（knuckles）基因座中的PcG誘導(dǎo)其轉(zhuǎn)錄，直至第6階段完全沉默WUS并終止花干細(xì)胞的分化。在此過程中，乙酰轉(zhuǎn)移酶AtGCN5可能參與調(diào)節(jié)AG和WUS的時(shí)空表達(dá)，同時(shí)H3K27me3去甲基化酶REF6通過與AG和SEP3互作促進(jìn)下游靶標(biāo)包括誘導(dǎo)雄蕊和心皮正常發(fā)育因子的轉(zhuǎn)錄[62];在第6階段，MADS通過招募染色質(zhì)重塑因子CHR4/11/17與SEPs、AP1、AG、PI和AP3等同源異性蛋白形成復(fù)合物調(diào)節(jié)長雄蕊房室的發(fā)育[63]。在第8～12階段，H3K4甲基轉(zhuǎn)移酶SDG2可能在雌雄配子體發(fā)育過程中介導(dǎo)基因的活性轉(zhuǎn)錄，而FIS2-PRC2復(fù)合物在配子發(fā)生和種子發(fā)育過程中的基因沉默中起作用。此外，H3K36甲基轉(zhuǎn)移酶SDG4在花粉管的形成中發(fā)揮重要作用。在花的轉(zhuǎn)變過程中，關(guān)鍵開花阻遏基因FLC（flowering locus C）需要被多梳響應(yīng)元件PRE和染色質(zhì)重塑因子INO80與H2A.Z的互作所沉默。此外，SWR1有利于FLC的轉(zhuǎn)錄激活，最新研究發(fā)現(xiàn)其組成亞基YAF9A在控制開花方面存在一條依賴H4和H2A.Z乙?；男滦驼{(diào)控通路[64]。

2.6? 果實(shí)成熟

近年研究發(fā)現(xiàn)m5C、m6A以及組蛋白甲基化在番茄果實(shí)成熟過程中建立了重要的分子聯(lián)系，其中非活性轉(zhuǎn)錄的DNA甲基化通過靶向m6A去甲基化酶ALKBH2調(diào)控m6A水平，ALKBH2則作用于ROS1高度同源的DML2對(duì)DNA甲基化進(jìn)行負(fù)反饋調(diào)節(jié)，同時(shí)DML2通過催化成熟誘導(dǎo)基因的去甲基化以激活它們的轉(zhuǎn)錄，包括參與細(xì)胞壁水解（PG2a、PL）、類胡蘿卜素（PSY1、Z-ISO、ZDS）以及乙烯生物合成和信號(hào)轉(zhuǎn)導(dǎo)的基因（ACS、ACO和ETR），從而促進(jìn)果實(shí)的成熟和軟化，ALKBH2是否同時(shí)靶向這些基因還有待進(jìn)一步研究[65-66];而組蛋白去甲基化酶JMJ6既可以通過降低DML2的H3K27me3水平促進(jìn)其表達(dá)，也可以直接靶向相關(guān)基因參與果實(shí)成熟的調(diào)節(jié)[67];此外，HDA1型組蛋白去乙?；窰DA1和HDA2的下調(diào)可促進(jìn)乙烯的生物合成和類胡蘿卜素的積累，進(jìn)而加速果實(shí)成熟，而HD2型去乙酰化酶HDT3的沉默卻發(fā)揮相反的作用，組蛋白乙酰化與甲基化是否存在聯(lián)系仍是未知[68-70]。草莓（Fragaria ananassa Duch.）果實(shí)的成熟也需要全基因組范圍的DNA低甲基化，但在此期間DML2并不被上調(diào)，DNA甲基轉(zhuǎn)移酶和RdDM活性的降低導(dǎo)致DNA去甲基化占據(jù)主導(dǎo)地位[71];香蕉（Musa nana Lour.）HDA1則通過與轉(zhuǎn)錄因子ERF11（ethylene response factor11）互作抑制成熟相關(guān)基因的表達(dá)[72]，提示誘導(dǎo)果實(shí)成熟的調(diào)控機(jī)制可能因其種類而異。

2.7? 胚胎發(fā)育

眾多研究表明，不同表觀遺傳調(diào)控類型介導(dǎo)了植物胚胎的發(fā)育，其中最主要的是DNA甲基化和組蛋白修飾。Bouyer等[73]研究發(fā)現(xiàn)成熟胚胎中CHH甲基化增加導(dǎo)致總DNA甲基化水平升高，CHH甲基化在早期和成熟胚胎中的差異在富含轉(zhuǎn)座子和基因缺乏的著絲粒周圍區(qū)域中最為明顯，提示胚胎發(fā)育過程中染色體存在不同的甲基化動(dòng)力學(xué)。對(duì)擬南芥DNA甲基化相關(guān)基因MET1和CMT3功能缺失突變體研究發(fā)現(xiàn)，其胚胎發(fā)育在細(xì)胞分裂平面和數(shù)量發(fā)生異常變化，且胚胎生存能力降低。兩種組蛋白去乙?；窰DAC6/ HDAC19也在胚胎發(fā)育中具有重要的調(diào)控作用，在擬南芥萌芽過程中，HDAC6/HDAC19通過抑制胚胎特異性基因功能直接或間接冗余的調(diào)控胚胎特性[74]。在種子發(fā)育期間，CHD類染色質(zhì)重塑因子CHR5和PKL通過拮抗作用于染色質(zhì)不同方面的修飾調(diào)節(jié)胚胎發(fā)育因子的轉(zhuǎn)錄，CHR5可通過克服核小體屏障與ABI3（abscisic acidinsensitive 3）和FUS3（fusca 3）的啟動(dòng)子區(qū)結(jié)合，并降低轉(zhuǎn)錄起始位點(diǎn)附近的核小體密度;PKL則通過促進(jìn)H3K27me3抑制相關(guān)基因的表達(dá)[75]。此外，種子發(fā)育的核心調(diào)控因子LEC1（leafy cotyledon 1）通過與不同轉(zhuǎn)錄因子互作參與調(diào)控胚胎形態(tài)發(fā)生和種子成熟等多個(gè)生物學(xué)過程。近期揭示了LEC1在早期胚胎中通過激活FLC重置春化狀態(tài)的分子機(jī)制。在胚胎發(fā)育早期，LEC1通過在FLC位點(diǎn)募集組蛋白甲基轉(zhuǎn)移酶EFS和染色質(zhì)重塑因子SWR1C促進(jìn)H3K36me3水平的升高和轉(zhuǎn)錄抑制標(biāo)記H3K27me3的消除，由此形成的FLC活性染色質(zhì)狀態(tài)持續(xù)到整個(gè)胚胎發(fā)育時(shí)期及幼苗期，隨后冬季低溫則通過招募PRC2復(fù)合體富集H3K27me3以沉默F(xiàn)LC活性狀態(tài)，并使這種沉默狀態(tài)維持到春季回暖期形成“低溫記憶”，這種記憶直至胚胎發(fā)育早期以LEC1激活FLC的形式被擦除[76]。

植物體胚發(fā)生過程中基因的時(shí)空表達(dá)同樣受到表觀遺傳修飾尤其是DNA甲基化的動(dòng)態(tài)調(diào)節(jié)。DNA甲基化模式的改變或修飾可參與植物體胚發(fā)生的啟動(dòng)和發(fā)育過程中多種候選標(biāo)記基因的調(diào)控。例如在輻射松（Pinus radiata D. Don）體細(xì)胞胚胎發(fā)生中，SERK1（somatic embryogenesis receptor-like kinase）、BBM（baby boom）和WOX2（WUS homeobox 2）基因的表達(dá)水平降低與DNA甲基化總體水平的增加相關(guān)[77]。此外，生長素可能以SAM和SAH（S-腺苷半胱氨酸）為媒介調(diào)控DNA甲基化。在培養(yǎng)基中去除2，4-D時(shí)可抑制乙烯的生物合成，引起SAM與SAH的濃度比值升高并最終導(dǎo)致DNA甲基化水平增加[78]。近期研究發(fā)現(xiàn)DNA和組蛋白甲基化在體胚誘導(dǎo)過程中存在不同的組合互作，其中RdDM和H3K9me2可能協(xié)同促進(jìn)陸地棉（Gossypium hirsutum）NEC（non-embryonic callus）到EC（embryonic callus）階段的CHH甲基化動(dòng)力學(xué)，而在EC至SE（somatic embryo）階段中伴隨著DNA甲基化和H3K9me2水平的降低;在中?？Х龋–offea canephora）中，DNA甲基化的減少和H3K9me2、H3K27me3水平的降低相關(guān)[79-80]。

總之，以DNA甲基化和組蛋白修飾為代表的表觀遺傳調(diào)控在植物胚胎發(fā)育過程中發(fā)揮著至關(guān)重要的作用。

3? 展望

表觀遺傳修飾廣泛地參與了植物的生長發(fā)育，重要發(fā)育因子中表觀遺傳標(biāo)記的動(dòng)態(tài)變化導(dǎo)致染色質(zhì)“可及”狀態(tài)的改變，進(jìn)而在轉(zhuǎn)錄或翻譯水平調(diào)節(jié)發(fā)育因子的表達(dá)。然而，表觀遺傳在植物響應(yīng)非生物及生物脅迫、胚胎發(fā)育、果實(shí)早期發(fā)育中的作用機(jī)制仍不清楚。例如，對(duì)非生物脅迫如何通過表觀遺傳調(diào)節(jié)染色質(zhì)狀態(tài)和隨后的轉(zhuǎn)錄重編程仍缺乏完整和清晰的認(rèn)知。總之，關(guān)于植物生長發(fā)育過程中的表觀遺傳密碼仍需要更深入的研究。在植物表觀遺傳的高通量測(cè)序和染色質(zhì)分析技術(shù)快速發(fā)展的背景下，越來越多的植物表觀基因組將會(huì)被破譯，這將大大促進(jìn)我們對(duì)模式植物及農(nóng)作物生長發(fā)育過程中表觀遺傳調(diào)控機(jī)制的研究、理解和應(yīng)用。植物表觀遺傳在農(nóng)業(yè)上應(yīng)用前景廣闊，研究表明，玉米、杉木（Cunninghamia lanceolata）和高粱的親本和雜種之間的特異基因組位點(diǎn)存在差異甲基化，水稻自交系之間也存在類似甲基化模式的差異，并且與F1代雜種某些DMRs的轉(zhuǎn)錄水平相關(guān)。這些發(fā)現(xiàn)支持了表觀遺傳調(diào)控與雜種優(yōu)勢(shì)之間的潛在聯(lián)系。最近，在龍眼上進(jìn)行了一系列體胚發(fā)生的關(guān)鍵基因包括RNA甲基化m6A修飾酶基因、組蛋白甲基化SDG基因家族的鑒定，以及DNA甲基化基因鑒定和甲基化水平的測(cè)定為利用表觀遺傳修飾的調(diào)控培育無核或焦核龍眼新品種的培育奠定了基礎(chǔ)。此外，龍眼母性效應(yīng)基因MEE70的鑒定也為穩(wěn)定無籽育種工作的開展提供了良好的理論基礎(chǔ)?？傊?，眾多關(guān)于表觀遺傳修飾在植物重要發(fā)育過程中的研究暗示了其在雜種優(yōu)勢(shì)利用、多倍體育種、抗性育種、開花與花形態(tài)調(diào)控、胚胎發(fā)育調(diào)控、果實(shí)成熟調(diào)控等領(lǐng)域具有巨大的應(yīng)用潛力，這將開辟新的育種思路，特別是結(jié)合基因編輯調(diào)控表觀遺傳的分子育種，將會(huì)顯示出巨大的應(yīng)用潛力。

參考文獻(xiàn)

Miro-Blanch J， Yanes O. Epigenetic regulation at the interplay between gut microbiota and host metabolism[J]. Frontiers in Genetics， 2019， 10： 638.

Greer E L， Blanco M A， Gu L， et al. DNA methylation on N6-adenine in C. elegans[J]. Cell， 2015， 161（4）： 868-878.

Ji L J， Chen X M. Regulation of small RNA stability： Methylation and beyond[J]. Cell research， 2012， 22（4）： 624-636.

Seymour D K， Becker C. The causes and consequences of DNA methylome variation in plants[J]. Current Opinion in Plant Biology， 2017， 36： 56-63.

袁? 超， 張少偉， 牛? 義， 等. 植物DNA甲基化作用機(jī)制的研究進(jìn)展[J]. 生物工程學(xué)報(bào)， 2020， 36（5）： 838-848.

Mozgova I， Hennig L. The polycomb group protein regulatory network[J]. Annual Review of Plant Biology， 2015， 66： 269-296.

Li S H， Ali S F， Duan X T， et al. JMJD1B demethylates H4R3me2s and H3K9me2 to facilitate gene expression for development of hematopoietic stem and progenitor cells[J]. Cell Reports， 2018， 23（2）： 389-403.

Roudier F， Ahmed I， Bérard C， et al. Integrative epigenomic mapping defines four main chromatin states in Arabidopsis[J]. The EMBO Journal， 2011， 30（10）： 1928-1938.

Xue H， Yao T H， Cao M， et al. Structural basis of nucleosome recognition and modification by MLL methyltransferases[J]. Nature， 2019， 573（7774）： 445-449.

Su Y H， Wang S L， Zhang F， et al. Phosphorylation of histone H2A at serine 95： a plant-specific mark involved in flowering time regulation and H2A. Z deposition[J]. The Plant Cell， 2017， 29（9）： 2197-2213.

陳榮珠. 龍眼體細(xì)胞胚胎發(fā)生早期AGO基因家族的全基因組鑒定、表達(dá)與功能分析[D]. 福州： 福建農(nóng)林大學(xué)， 2020.

Pendleton K E， Chen B B， Liu K Q， et al. The U6 snRNA m6A methyltransferase METTL16 regulates SAM synthetase intron retention[J]. Cell， 2017， 169（5）： 824-835.

王靜宇， 陳曉慧， 申? 序， 等. 龍眼體胚發(fā)生過程中RNA甲基化相關(guān)基因全基因組鑒定及表達(dá)分析[J]. 熱帶作物學(xué)報(bào)， 2019， 40（10）： 1902-1913.

Cui X A， Liang Z， Shen L S， et al. 5-Methylcytosine RNA methylation in Arabidopsis thaliana[J]. Molecular plant， 2017， 10（11）： 1387-1399.

Goll M G， Kirpekar F， Maggert K A， et al. Methylation of tRNAAsp by the DNA methyltransferase homolog Dnmt2[J]. Science， 2006， 311（5759）： 395-398.

Fu L J， Guerrero C R， Zhong N， et al. Tet-mediated formation of 5-hydroxymethylcytosine in RNA[J]. Journal of the American Chemical Society， 2014， 136（33）： 11582-11585.

Vilardo E， Nachbagauer C， Buzet A， et al. A subcomplex of human mitochondrial RNase P is a bifunctional methyltransferase-extensive moonlighting in mitochondrial tRNA biogenesis[J]. Nucleic Acids Research， 2012， 40（22）： 11583- 11593.

Kawarada L， Suzuki T， Ohira T， et al. ALKBH1 is an RNA dioxygenase responsible for cytoplasmic and mitochondrial tRNA modifications[J]. Nucleic Acids Research， 2017， 45（12）： 7401-7415.

Clapier C R， Cairns B R. The biology of chromatin remodeling complexes[J]. Annual Review of Biochemistry， 2009， 78： 273-304.

Han S K， Wu M F， Cui S J， et al. Roles and activities of chromatin remodeling ATP ases in plants[J]. The Plant Journal， 2015， 83（1）： 62-77.

Zhang J J， Lai J B， Wang F G， et al. A SUMO ligase AtMMS21 regulates the stability of the chromatin remodeler BRAHMA in root development[J]. Plant Physiology， 2017， 173（3）： 1574-1582.

Torres E S， Deal R B. The histone variant H2A. Z and chromatin remodeler BRAHMA act coordinately and antagonistically to regulate transcription and nucleosome dynamics in Arabidopsis[J]. The Plant Journal， 2019， 99（1）： 144-162.

Li G， Liu S J， Wang J W， et al. ISWI proteins participate in the genome-wide nucleosome distribution in Arabidopsis[J]. The Plant Journal， 2014， 78（4）： 706-714.

Zou B H， Sun Q， Zhang W L， et al. The Arabidopsis chromatin-remodeling factor CHR5 regulates plant immune responses and nucleosome occupancy[J]. Plant and Cell Physiology， 2017， 58（12）： 2202-2216.

Carter B， Bishop B， Ho K K， et al. The chromatin remodelers PKL and PIE1 act in an epigenetic pathway that determines H3K27me3 homeostasis in Arabidopsis[J]. The Plant Cell， 2018， 30（6）： 1337-1352.

Kang H J， Zhang C， An Z X， et al. AtINO80 and AtARP5 physically interact and play common as well as distinct roles in regulating plant growth and development[J]. New Phytologist， 2019， 223（1）： 336-353.

Potok M E， Wang Y F， Xu L H， et al. Arabidopsis SWR1-associated protein methyl-CpG-binding domain 9 is required for histone H2A. Z deposition[J]. Nature Communications， 2019， 10（4）： 301-311.

Wang K C， Chang H Y. Molecular mechanisms of long noncoding RNAs[J]. Molecular Cell， 2011， 43（6）： 904-914.

Jullien P E， Berger F. Gamete-specific epigenetic mechanisms shape genomic imprinting[J]. Current Opinion in Plant Biology， 2009， 12（5）： 637-642.

Makarevich G， Villar C B R， Erilova A， et al. Mechanism of PHERES1 imprinting in Arabidopsis[J]. Journal of Cell Science， 2008， 121（6）： 906-912.

Dumbliauskas E， Lechner E， Jaciubek M， et al. The Arabidopsis CUL4-DDB1 complex interacts with MSI1 and is required to maintain MEDEA parental imprinting[J]. The EMBO Journal， 2011， 30（4）： 731-743.

盧秉國， 練從龍， 賴鐘雄. 龍眼MEE70/MSI1基因DlMEE70-1a及其可變剪接體DlMEE70-1b的分離及其在體胚發(fā)生過程中的表達(dá)分析[J]. 熱帶作物學(xué)報(bào)， 2015， 36（1）： 59-67.

Pathania R， Ramachandran S， Elangovan S， et al. DNMT1 is essential for mammary and cancer stem cell maintenance and tumorigenesis[J]. Nature Communications， 2015， 6（1）： 6910.

Yin J， Leavenworth J W， Li Y， et al. Ezh2 regulates differentiation and function of natural killer cells through histone methyltransferase activity[J]. Proceedings of the National Academy of Sciences， 2015， 112（52）： 15988-15993

Sheng L H， Hu X M， Du Y J， et al. Non-canonical WOX11-mediated root branching contributes to plasticity in Arabidopsis root system architecture[J]. Development， 2017， 144（17）： 3126-3133.

Zhang Z， Wang B W， Wang S H， et al. Genome-wide target mapping shows histone deacetylase complex1 regulates cell proliferation in cucumber fruit[J]. Plant Physiology， 2020， 182（1）： 167-184.

El-Tantawy A A， Solís M T， Risue?o M C， et al. Changes in DNA methylation levels and nuclear distribution patterns after microspore reprogramming to embryogenesis in barley[J]. Cytogenetic and Genome Research， 2014， 143（1/3）： 200- 208.

Takatsuka H， Umeda M. Epigenetic control of cell division and cell differentiation in the root apex[J]. Frontiers in Plant Science， 2015， 6： 1178.

Aichinger E， Villar C B R， Di Mambro R， et al. The CHD3 chromatin remodeler PICKLE and polycomb group proteins antagonistically regulate meristem activity in the Arabidopsis root[J]. The Plant Cell， 2011， 23（3）： 1047-1060.

Li H C， Torres-Garcia J， Latrasse D， et al. Plant-specific histone deacetylases HDT1/2 regulate GIBBERELLIN 2-OXIDASE2 expression to control Arabidopsis root meristem cell number[J]. The Plant Cell， 2017， 29（9）： 2183-2196.

Liu X G， Kim Y J， Müller R， et al. AGAMOUS terminates floral stem cell maintenance in Arabidopsis by directly repressing WUSCHEL through recruitment of Polycomb Group proteins[J]. The Plant Cell， 2011， 23（10）： 3654-3670.

Aichinger E， Kornet N， Friedrich T， et al. Plant stem cell niches[J]. Annual Review of Plant Biology， 2012， 63： 615-636.

Yu A， Lepère G， Jay F， et al. Dynamics and biological relevance of DNA demethylation in Arabidopsis antibacterial defense[J]. Proceedings of the National Academy of Sciences， 2013， 110（6）： 2389-2394.

Geng S F， Kong X C， Song G Y， et al. DNA methylation dynamics during the interaction of wheat progenitor Aegilops tauschii with the obligate biotrophic fungus Blumeria graminis f. sp. tritici[J]. New Phytologist， 2019， 221（2）： 1023-1035.

López A， Ramírez V， García-Andrade J， et al. The RNA silencing enzyme RNA polymerase V is required for plant immunity[J]. PLoS Genetic， 2011， 7（12）： e1002434.

Li X Y， Jiang Y X， Ji Z C， et al. BRHIS 1 suppresses rice innate immunity through binding to monoubiquitinated H2A and H2B variants[J]. EMBO Reports， 2015， 16（9）： 1192- 1202.

Lee S H， Fu F Y， Xu S M， et al. Global regulation of plant immunity by histone lysine methyl transferases[J]. The Plant Cell， 2016， 28（7）： 1640-1661.

Yang R， Hong Y C， Ren Z Z， et al. A role for PICKLE in the regulation of cold and salt stress tolerance in Arabidopsis[J]. Frontiers in Plant Science， 2019， 10： 900.

Park J， Lim C J， Shen M Z， et al. Epigenetic switch from repressive to permissive chromatin in response to cold stress[J]. Proceedings of the National Academy of Sciences， 2018， 115（23）： E5400-E5409.

Butelli E， Licciardello C， Zhang Y， et al. Retrotransposons control fruit-specific， cold-dependent accumulation of anthocyanins in blood oranges[J]. The Plant Cell， 2012， 24（3）： 1242-1255.

L?mke J， Brzezinka K， Altmann S， et al. A hit-and-run heat shock factor governs sustained histone methylation and transcriptional stress memory[J]. The EMBO Journal， 2016， 35（2）： 162-175.

Sanchez D H， Paszkowski J. Heat-induced release of epigenetic silencing reveals the concealed role of an imprinted plant gene[J]. PLoS Genetic， 2014， 10（11）： e1004806.

Fragkostefanakis S， Roeth S， Schleiff E， et al. Prospects of engineering thermotolerance in crops through modulation of heat stress transcription factor and heat shock protein networks[J]. Plant， Cell & Environment， 2015， 38（9）： 1881- 1895.

Xu R， Wang Y H， Zheng H， et al. Salt-induced transcription factor MYB74 is regulated by the RNA-directed DNA methylation pathway in Arabidopsis[J]. Journal of Experimental Botany， 2015， 66（19）： 5997-6008.

Kumar S， Beena A S， Awana M， et al. Salt-induced tissue-specific cytosine methylation downregulates expression of HKT genes in contrasting wheat （Triticum aestivum L.） genotypes[J]. DNA and Cell Biology， 2017， 36（4）： 283-294.

Sani E， Herzyk P， Perrella G， et al. Hyperosmotic priming of Arabidopsis seedlings establishes a long-term somatic memory accompanied by specific changes of the epigenome[J]. Genome Biology， 2013， 14（6）： R59.

Mlynárová L， Nap J P， Bisseling T. The SWI/SNF chromatin-remodeling gene AtCHR12 mediates temporary growth arrest in Arabidopsis thaliana upon perceiving environmental stress[J]. The Plant Journal， 2007， 51（5）： 874-885.

Kim J M， To T K， Seki M. An epigenetic integrator： New insights into genome regulation， environmental stress responses and developmental controls by histone deacetylase 6[J]. Plant and Cell Physiology， 2012， 53（5）： 794-800.

Li S， Lin Y C， Wang P Y， et al. The AREB1 transcription factor influences histone acetylation to regulate drought responses and tolerance in Populus trichocarpa[J]. The Plant Cell， 2019， 31（3）： 663-686.

Mao H D， Wang H W， Liu S X， et al. A transposable element in a NAC gene is associated with drought tolerance in maize seedlings[J]. Nature Communications， 2015， 6（1）： 8326.

Smaczniak C， Immink R G H， Mui?o J M， et al. Characterization of MADS-domain transcription factor complexes in Arabidopsis flower development[J]. Proceedings of the National Academy of Sciences， 2012， 109（5）： 1560-1565.

Sun B， Looi L S， Guo S Y， et al. Timing mechanism dependent on cell division is invoked by polycomb eviction in plant stem cells[J]. Science， 2014， 343（6170）： 1248559.

Li G， Liu S J， Wang J W， et al. ISWI proteins participate in the genome-wide nucleosome distribution in Arabidopsis[J]. The Plant Journal， 2014， 78（4）： 706-714.

Crevillén P， Gómez-Zambrano ?， López J A， et al. Arabidopsis YAF9 histone readers modulate flowering time through NuA4-complex-dependent H4 and H2A. Z histone acetylation at FLC chromatin[J]. New Phytologist， 2019， 222（4）： 1893-1908.

Lang Z B， Wang Y H， Tang K， et al. Critical roles of DNA demethylation in the activation of ripening-induced genes and inhibition of ripening-repressed genes in tomato fruit[J]. Proceedings of the National Academy of Sciences， 2017， 114（22）： E4511-E4519.

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.

Li Z W， Jiang G X， Liu X C， et al. Histone demethylase SlJMJ6 promotes fruit ripening by removing H3K27 methylation of ripening-related genes in tomato[J]. New Phytologist， 2020， 227（4）： 1138-1156.

Guo J E， Hu Z L， Li F F， et al. Silencing of histone deacetylase SlHDT3 delays fruit ripening and suppresses carotenoid accumulation in tomato[J]. Plant Science， 2017， 265： 29-38.

Guo J E， Hu Z L， Yu X H， et al. A histone deacetylase gene， SlHDA3， acts as a negative regulator of fruit ripening and carotenoid accumulation[J]. Plant Cell Reports， 2018， 37（1）： 125-135.

Guo J E， Hu Z L， Zhu M K， et al. The tomato histone deacetylase SlHDA1 contributes to the repression of fruit ripening and carotenoid accumulation[J]. Scientific Reports， 2017， 7（1）： 7930.

Cheng J F， Niu Q F， Zhang B， et al. Downregulation of RdDM during strawberry fruit ripening[J]. Genome Biology， 2018， 19： 212.

Han Y C， Kuang J F， Chen J Y， et al. Banana transcription factor MaERF11 recruits histone deacetylase MaHDA1 and represses the expression of MaACO1 and expansins during fruit ripening[J]. Plant Physiology， 2016， 171（2）： 1070-1084.

Bouyer D， Kramdi A， Kassam M， et al. DNA methylation dynamics during early plant life[J]. Genome biology， 2017， 18（1）： 179.

Tanaka M， Kikuchi A， Kamada H. The Arabidopsis histone deacetylases HDA6 and HDA19 contribute to the repression of embryonic properties after germination[J]. Plant physiology， 2008， 146（1）： 149-161.

Shen Y， Devic M， Lepiniec L， et al. Chromodomain， Helicase and DNA-binding CHD1 protein， CHR5， are involved in establishing active chromatin state of seed maturation genes[J]. Plant Biotechnology Journal， 2015， 13（6）： 811-820.

Tao Z， Shen L S， Gu X F， et al. Embryonic epigenetic reprogramming by a pioneer transcription factor in plants[J]. Nature， 2017， 551（7678）： 124-128.

Karim R， Tan Y S， Singh P， et al. Expression and DNA methylation of SERK， BBM， LEC2 and WUS genes in in vitro cultures of Boesenbergia rotunda （L.） Mansf[J]. Physiology and Molecular Biology of Plants， 2018， 24（5）： 741-751.

Jo L， Dos Santos A L W， Bueno C A， et al. Proteomic analysis and polyamines， ethylene and reactive oxygen species levels of Araucaria angustifolia （Brazilian pine） embryogenic cultures with different embryogenic potential[J]. Tree Physiology， 2014， 34（1）： 94-104.

Li J Y， Wang M J， Li Y J， et al. Multi-omics analyses reveal epigenomics basis for cotton somatic embryogenesis through successive regeneration acclimation process[J]. Plant Biotechnology Journal， 2019， 17（2）： 435-450.

Nic-Can G I， Lopez-Torres A， Barredo-Pool F， et al. New insights into somatic embryogenesis： LEAFY COTYLEDON1， BABY BOOM1 and WUSCHEL-RELATED HOMEOBOX4 are epigenetically regulated in Coffea canephora[J]. PLoS One， 2013， 8（8）： e72160.