張圓圓　劉文敬　張斌斌　馬瑞娟　俞明亮

摘要： 環(huán)氧化物水解酶（EH）因其重要的生物學功能而在哺乳動物以及諸多植物中被廣泛關(guān)注，更是果實典型“桃香”氣味——內(nèi)酯芳香物質(zhì)生物合成的一個重要酶，但在桃等果實中的研究較少且鮮有該家族成員的系統(tǒng)報道或生物學功能的解析。為鑒別桃果實中與內(nèi)酯芳香物質(zhì)合成相關(guān)的環(huán)氧化物水解酶家族成員，本研究使用了同源序列比對和關(guān)鍵詞搜索等方法，在桃中共篩選獲得7個環(huán)氧化物水解酶家族成員。序列比對分析結(jié)果表明，這7個成員均有典型的α/β水解酶折疊結(jié)構(gòu)和環(huán)氧化物水解酶保守的序列片段。進化樹分析結(jié)果顯示，桃的環(huán)氧化物水解酶成員與擬南芥、煙草等其他物種中已被鑒別的環(huán)氧化物水解酶成員的親緣關(guān)系很近。基因表達分析結(jié)果顯示，在桃的發(fā)育成熟進程中7個EH基因均在果實的中果皮表達，在果實發(fā)育前期均呈現(xiàn)較高的表達水平，轉(zhuǎn)錄模式包括3種類型。綜合已報道的內(nèi)酯芳香物質(zhì)的生物合成通路以及前人指出的內(nèi)酯芳香物質(zhì)含量在果實成熟期開始顯著增加的變化規(guī)律，推測桃環(huán)氧化物水解酶成員的表達量可能與內(nèi)酯芳香物質(zhì)的積累負相關(guān)。本研究結(jié)果為后續(xù)深入挖掘桃或者其他果實環(huán)氧化物水解酶家族成員的生物學功能尤其是參與內(nèi)酯芳香物質(zhì)合成的分子機理提供了參考。

關(guān)鍵詞： 桃；內(nèi)酯芳香物質(zhì)；環(huán)氧化物水解酶

中圖分類號： S662.1 文獻標識碼： A 文章編號： 1000-4440（2023）01-0178-09

Identification of candidate epoxide hydrolase genes involved in the biosynthesis of lactone volatile compounds in peach （Prunus persica L.）

ZHANG Yuan-yuan¹， LIU Wen-jing^1，2， ZHANG Bin-bin¹， MA Rui-juan¹， YU Ming-liang¹

（1.Institute of Pomology， Jiangsu Academy of Agricultural Sciences/Jiangsu Key Laboratory for Horticultural Crop Genetic Improvement， Nanjing 210014， China；2.School of Food and Biological Engineering， Jiangsu University， Zhenjiang 212013， China）

Abstract： Epoxide hydrolases （EHs） were extensively characterized in mammals and numerous plant species for their multiple biological functions， and acted as one important enzyme in the biosynthesis of lactones——a group of volatile compounds which endowed fruit a typical peach-like aroma. However， there were few reports on the family members or biological functions of epoxide hydrolases in peach or other fruits. To identify the EH members involved in the biosynthesis of lactone volatile compounds in peach fruit， homologous sequence alignment and keyword search were used in this study. A total of seven EH family members were screened in peach， and sequence alignment analysis showed that they all had the typical α/β-hydrolase fold structure and conserved sequences. Phylogenetic tree analysis showed that EH members in peach were closely related to the identified EH members in Arabidopsis， tobacco and other plant species. The results of gene expression analysis indicated that these EH genes were expressed in the mesocarp of peach fruit during the development and maturation， and showed a high expression level in the early stages of fruit development. The transcription patterns included three types. Based on the reported biosynthetic pathways of lactones and the fact that the content of lactones began to increase significantly at fruit ripening stages， it was speculated that the expression levels of the epoxide hydrolase genes may be negatively correlated with the contents of lactones volatile compounds. Overall， these results lay a foundation for further exploration of the biological functions of EH family members in peach or other fruits， especially the molecular mechanism involved in the biosynthesis of lactone volatile compounds.

Key words： peach；lactones volatile compounds；epoxide hydrolases

在桃果實上百種揮發(fā)性芳香物質(zhì)中，內(nèi)酯是賦予其典型“桃香”氣味的重要貢獻組分，亦是桃果實芳香品質(zhì)研究和改良的重點。目前已知的內(nèi)酯合成途徑包括經(jīng)由飽和脂肪酸或不飽和脂肪酸的2條支路^[1]，其中不飽和脂肪酸支路通過一個重要酶——環(huán)氧化物水解酶（EH）的作用產(chǎn)生羥基脂肪酸，最終在自身環(huán)化或醇酰基轉(zhuǎn)移酶的作用下產(chǎn)生內(nèi)酯。Schttler和Boland^[2]指出環(huán)氧化物水解酶在成熟桃和草莓果實的內(nèi)酯芳香組分的形成中發(fā)揮關(guān)鍵作用，該酶在油桃果實中具有催化活性且催化產(chǎn)物通過后續(xù)代謝最終形成了內(nèi)酯物質(zhì)，由此推測不飽和脂肪酸如亞油酸和亞麻酸的環(huán)氧化作用可能是果實中發(fā)生飽和脂肪酸衍生物氧化的一種普遍途徑，并最終產(chǎn)生了γ-內(nèi)酯物質(zhì)和δ-內(nèi)酯物質(zhì)。

環(huán)氧化物水解酶存在于所有生物中，是一類能立體選擇地將水分子加成到環(huán)氧底物上生成相應的1，2-二醇類的酶，具有α/β水解酶折疊蛋白質(zhì)結(jié)構(gòu)，其活性位點包括核心結(jié)構(gòu)（通常由2個天冬氨酸殘基和1個組氨酸殘基構(gòu)成）和帽子結(jié)構(gòu)^[3-5]。由于諸多被EH催化的脂類基質(zhì)生物活性強，因此EH在植物和動物中具有重要而多樣的生物學作用，如參與宿主防御和發(fā)育控制等，對生理系統(tǒng)有著深遠的影響。根據(jù)細胞定位的不同和底物的特異性，高等生物中不同的環(huán)氧化物水解酶被分為7種亞型，包括可溶性環(huán)氧化物水解酶（sEH）、微粒體環(huán)氧化物水解酶（mEH）、保幼激素環(huán)氧化物水解酶（JhEH）、膽固醇環(huán)氧化物水解酶（ChEH）、羥環(huán)氧烯酸水解酶、白三烯A4環(huán)氧化物水解酶（LAH）和檸檬烯環(huán)氧化物水解酶（LEH），且不同類型和來源的環(huán)氧化物水解酶的蛋白質(zhì)相對分子質(zhì)量大小明顯不同^[6-8]。

環(huán)氧化物水解酶廣泛分布于發(fā)芽的種子、根、果實、塊莖和葉等組織中，并在大豆、綠豆、擬南芥、馬鈴薯、普通煙草等多種植物中被廣泛研究^[9-18]。在植物環(huán)氧化物水解酶中，有關(guān)sEH的相關(guān)研究目前報道較多。sEH定位于細胞質(zhì)基質(zhì)和乙醛酸循環(huán)體，以羥基酸和環(huán)氧脂肪酸等為催化底物，蛋白質(zhì)相對分子質(zhì)量約35 000，一般以單體或二聚體蛋白質(zhì)的形式存在^[4，6]。在模式植物擬南芥中，目前已有EH成員AtSEH和AtEH1等獲得報道，為其他植物EH的研究提供了借鑒。Kiyosue等^[10]從干旱脅迫處理的擬南芥中分離獲得AtSEH基因的cDNA序列，研究發(fā)現(xiàn)該基因的表達受到生長素和干旱脅迫的誘導并呈現(xiàn)組織特異性表達模式，其編碼的可溶性環(huán)氧化物水解酶的蛋白質(zhì)相對分子質(zhì)量預測為36 423。Pineau等^[18]的研究結(jié)果表明，AtEH1酶蛋白定位于細胞質(zhì)，具有環(huán)氧化物水解酶保守的核心催化位點，參與了角質(zhì)的形成。除擬南芥外，馬鈴薯也是植物中EH研究較為深入的物種。Mowbray等^[13]通過X射線衍射的方法揭示了馬鈴薯環(huán)氧化物水解酶StEH1的蛋白質(zhì)結(jié)構(gòu)，這是植物環(huán)氧化物水解酶結(jié)構(gòu)首次被報道，并促進了其后馬鈴薯EH及其變種蛋白質(zhì)結(jié)構(gòu)研究的逐步深入^{[14，17]}。此外，本氏煙草中過氧化物酶體靶向的EH成員NbEH2.1和NbEH2.2能夠特異地參與細菌等病原體的防御^[15]；環(huán)氧化物水解酶NbEH1.1和NbEH2.1來自2個進化分支，生化功能分析結(jié)果表明，這2種酶具有不同的底物特異性，EH1可能通過在脅迫環(huán)境下產(chǎn)生信號物質(zhì)激活植物的應激反應來發(fā)揮功能，EH2則可能參與角質(zhì)單體的形成^[16]。

近年來，環(huán)氧化物水解酶作為桃果實內(nèi)酯芳香物質(zhì)合成的參與酶在多項研究中獲得關(guān)注。Vecchietti等^[19]創(chuàng)建了多基因型多發(fā)育階段的桃外果皮和中果皮的表達序列標簽數(shù)據(jù)庫，并結(jié)合公共數(shù)據(jù)庫的桃表達序列標簽信息，通過綜合分析芳香物質(zhì)的合成和轉(zhuǎn)錄模式篩選出了與芳香相關(guān)的重要基因，其中包含1個內(nèi)酯合成的關(guān)鍵基因——環(huán)氧化物水解酶編碼基因，可從果實材料中克隆獲得該環(huán)氧化物水解酶編碼基因1 350 bp的全長序列?；贐olero桃與OroA桃的芳香差異極大的特點，Pirona等^[20]對這2種基因型的桃構(gòu)建了微陣列數(shù)據(jù)以分析桃果實在3個發(fā)育成熟階段中基因的表達量變化，發(fā)現(xiàn)2個環(huán)氧化物水解酶基因EPH2（ppa008854m）和EPH3（ppa009153m）在2種桃的發(fā)育成熟階段呈現(xiàn)差異表達模式。Li等^[21]以溶質(zhì)桃為試驗材料通過基因表達分析鑒定與芳香相關(guān)的基因，發(fā)現(xiàn)ppa007004m、ppa008918m、ppa008854m、ppa008756m、ppa013585m等環(huán)氧化物水解酶基因在果實的不同發(fā)育階段以及4個品種之間呈現(xiàn)差異表達。

目前，有關(guān)桃果實內(nèi)酯合成相關(guān)環(huán)氧化物水解酶基因的報道呈零星分布，多由差異表達基因分析獲得，且不同研究涉及的環(huán)氧化物水解酶基因成員之間存在重疊，鮮有文獻系統(tǒng)報道桃EH基因家族成員的構(gòu)成和序列特征或深入解析其具體的生物學功能等。本研究篩選獲得了桃EH家族候選成員，對候選成員編碼的氨基酸序列和蛋白質(zhì)結(jié)構(gòu)進行了分析，并分析了果實發(fā)育成熟過程中各成員基因的表達量變化，鑒別桃果實中內(nèi)酯合成相關(guān)的候選環(huán)氧化物水解酶成員，以期為深入研究果實環(huán)氧化物水解酶的生物學功能和果實內(nèi)酯芳香物質(zhì)合成的分子機理提供借鑒。

1 材料與方法

1.1 試驗材料

本研究以2020年取自國家果樹種質(zhì)南京桃資源圃的白花水蜜（Prunus persica L.）的桃果實為試驗材料。對S1～S4發(fā)育成熟階段（分別對應盛花后43 d，77 d，115 d，144 d）的桃果實依次進行采收，選取大小和色澤一致且無機械傷的果實用于研究。每個生物學重復選取6個果實，設置3組生物學重復。中果皮部位切碎后速凍于液氮，樣品存放于-80 ℃冰箱備用。

1.2 總RNA提取與cDNA合成

使用RNAprep Pure試劑盒[天根生化科技（北京）有限公司產(chǎn)品]提取S1～S4發(fā)育成熟階段桃果實樣品的總RNA。凝膠電泳檢測RNA的完整性。

使用PrimeScript反轉(zhuǎn)錄試劑盒[寶生物工程（大連）有限公司產(chǎn)品]去除基因組DNA殘留，并對1 μg的RNA進行cDNA合成。合成的cDNA經(jīng)過DEPC（焦碳酸二乙酯）稀釋后作為后續(xù)實時熒光定量PCR的模板。

1.3 桃EH家族候選成員搜索與蛋白質(zhì)保守結(jié)構(gòu)域分析

基于JGI（https：//phytozome-next.jgi.doe.gov/）和NCBI（https：//www.ncbi.nlm.nih.gov/）數(shù)據(jù)庫，參考已報道的桃（Prunus persica L.）和擬南芥（Arabidopsis thaliana L.）的EH家族成員^{[10，18，20-21]}的蛋白質(zhì)氨基酸序列，通過同源序列比對、α/β水解酶結(jié)構(gòu)家族成員搜索、可溶性環(huán)氧化物水解酶的關(guān)鍵詞搜索等方法獲得桃EH家族的候選成員；使用NCBI數(shù)據(jù)庫的CCD在線工具（https：//www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi）分析桃EH家族候選成員的保守結(jié)構(gòu)域（表1）。

1.4 氨基酸序列與蛋白質(zhì)結(jié)構(gòu)分析

利用ISPript3.0在線工具（https：//espript.ibcp.fr/ESPript/ESPript/index.php）^[22]進行桃EH候選成員與其他物種已報道的EH成員的氨基酸序列的比對與蛋白質(zhì)結(jié)構(gòu)的分析。

1.5 系統(tǒng)進化樹構(gòu)建

采用Clustalx軟件進行桃EH候選成員的蛋白質(zhì)氨基酸序列和其他物種已報道的EH成員的比對，采用Figtree （version 1.3.1）軟件進行系統(tǒng)進化樹的構(gòu)建。

1.6 RT-qPCR分析

采用實時熒光定量PCR方法，利用CFX96實時熒光定量PCR儀[伯樂生命醫(yī)學產(chǎn)品（上海）有限公司產(chǎn)品]檢測桃EH家族候選基因在S1～S4發(fā)育成熟階段的表達量變化。反應體系包括SYBR PCR 預混液[伯樂生命醫(yī)學產(chǎn)品（上海）有限公司產(chǎn)品]10.0 μl，RT-qPCR上下游引物各1.0 μl，滅菌DEPC水6.0 μl，稀釋的cDNA模板2.0 μl。根據(jù)EH基因成員的編碼區(qū)堿基序列，利用NCBI/Primer-BLAST在線設計各基因相應的RT-qPCR引物，經(jīng)過熔解曲線分析和回收產(chǎn)物測序分析等方法檢驗引物的擴增特異性^[23]。以PpTEF2基因為內(nèi)參基因^[24]，設定S1階段的各基因相對表達量為1，使用2^-△△Ct方法計算相對表達量。RT-qPCR引物序列見表2。

1.7 統(tǒng)計分析與圖表制作

使用Microsoft excel 2019計算標準差，使用Origin8.0進行相關(guān)圖表的制作。

2 結(jié)果與分析

2.1 桃EH家族候選成員的篩選和氨基酸序列比對

通過蛋白質(zhì)保守結(jié)構(gòu)域分析，首先從相關(guān)序列中篩選出屬于α/β水解酶結(jié)構(gòu)超級家族的成員，它們均具有保守的MhpC（Pimeloyl-ACP methyl ester carboxylesterase）結(jié)構(gòu)域。隨后參考馬鈴薯環(huán)氧化物水解酶StEH1的蛋白質(zhì)結(jié)構(gòu)（蛋白質(zhì)PDB編號：2CJP），通過氨基酸序列比對和蛋白質(zhì)結(jié)構(gòu)分析鑒別出7個蛋白質(zhì)具有環(huán)氧化物水解酶保守的α/β水解酶折疊結(jié)構(gòu)域，將其命名為PpEH1～PpEH7，其編碼的基因為桃EH家族的候選基因（表1）。桃EH家族候選成員所編碼的蛋白質(zhì)中，蛋白質(zhì)氨基酸序列最短為313個氨基酸（PpEH7），最長為339個氨基酸（PpEH1）。這些成員的MhpC結(jié)構(gòu)域的位置存在差異。

將桃與已報道的其他物種EH成員的氨基酸序列進行比對分析，發(fā)現(xiàn)它們的α/β水解酶結(jié)構(gòu)呈現(xiàn)較大的相似性，氨基酸序列的平均相似度達到60.35%。桃的EH成員均具有標準的α/β水解酶折疊結(jié)構(gòu)，8條β片層被多個α螺旋結(jié)構(gòu)包圍，而高度保守的三聯(lián)體組成了催化殘基，包括位于β5之后的親核氨基酸（天冬氨酸），位于最后一個β片層之后的負責活化水分子的組氨酸以及輔佐組氨酸的位于β7片層之后的酸性催化殘基（天冬氨酸）（圖1）。三位一體的結(jié)構(gòu)之外，2個氧陰離子洞位置的氨基酸在桃EH成員之間非常保守，包括β3附近的色氨酸和β5親核天冬氨酸旁的色氨酸（PpEH3）或苯丙氨酸。同時，所有比對的EH成員在β4結(jié)束的位置上有保守的RGYG（D/L）（S/T）序列片段。

2.2 桃EH家族候選成員的蛋白質(zhì)進化樹分析

利用Figtree1.3.1軟件構(gòu)建桃和其他植物EH成員的蛋白質(zhì)系統(tǒng)進化樹（圖2）。在進化樹的3個分支中，擬南芥的2個EH成員（AtEH4和AtEH5）位于第一分支，馬鈴薯、粗皮檸檬、大豆和油菜的4個EH蛋白與擬南芥的2個EH成員（AtEH2和AtEH3）共同位于第二分支，桃的7個EH成員、擬南芥的2個EH成員（AtEH1和AtEH6）和煙草的EH成員共同聚類在第三分支。

在第二分支，馬鈴薯和粗皮檸檬的EH成員聚類在一起，大豆、油菜和擬南芥的2個EH成員聚類在一起，具有較近的親緣關(guān)系。在第三分支，PpEH1、PpEH4、PpEH5、PpEH6、PpEH7和煙草EH成員聚類在一起，PpEH2、PpEH3和擬南芥AtEH6聚類在一起，且上述EH成員均與AtEH1親緣關(guān)系很近。

2.3 桃EH家族候選基因在果實發(fā)育成熟階段的表達模式

從果實的第一次快速膨大期至成熟期，7個EH基因均有表達，在發(fā)育成熟階段前期呈現(xiàn)了較高的表達水平，在整個發(fā)育成熟進程中呈現(xiàn)整體逐漸下調(diào)的表達趨勢（圖3）。

桃EH家族候選基因表達變化分為3種模式。模式Ⅰ表現(xiàn)為3個基因（PpEH1、PpEH3、PpEH6）的表達量伴隨果實發(fā)育成熟進程逐步下調(diào)，模式Ⅱ表現(xiàn)為1個基因（PpEH2）的表達量在發(fā)育成熟階段的前期上調(diào)而后期下調(diào)，模式Ⅲ表現(xiàn)為3個基因（PpEH4、PpEH5、PpEH7）的表達量在發(fā)育成熟階段的前期有下調(diào)和上調(diào)趨勢而隨后下調(diào)，整體呈現(xiàn)下調(diào)趨勢。

表達模式Ⅰ的3個基因在S1階段的表達量最高，而在S4階段的表達量下調(diào)至S1階段的4.5%以下，是整個發(fā)育成熟階段表達量的最低值。模式Ⅱ的EH基因PpEH2存在一個上調(diào)表達過程，最高表達量出現(xiàn)在S2階段，最低表達量出現(xiàn)在S4階段，S4階段的基因表達量下調(diào)至S2表達量的2.3%。模式Ⅲ中，PpEH4在S4階段的基因表達量下調(diào)至最高表達量（S3階段）的37.2%，PpEH5在S4階段的基因表達量下調(diào)至最高表達量（S3階段）的31.9%；PpEH7在S4階段的基因表達量下調(diào)至最高表達量（S1階段）的8.5%，下調(diào)幅度最大。綜合上述分析，桃的EH家族基因在S4階段之前均具有較高的表達水平，在成熟期S4階段的基因表達水平相對前期呈現(xiàn)顯著的下調(diào)；結(jié)合前人研究所報道的內(nèi)酯芳香物質(zhì)的含量在成熟后顯著增加的規(guī)律，可知這些EH基因表達量可能與內(nèi)酯芳香物質(zhì)的積累負相關(guān)。

3 討論

3.1 環(huán)氧化物水解酶的結(jié)構(gòu)特征

環(huán)氧化物水解酶具有α/β水解酶折疊結(jié)構(gòu)，屬于包含有酯酶、蛋白酶、脂肪酶、脫鹵酶、裂解酶等功能的α/β水解酶超級家族中的一類^[25-26]。標準的α/β水解酶折疊是由幾乎平行的8條β片層（其中只有第二個β片層與其他片層反向平行）和包圍在兩邊的多個α螺旋共同組成的，且α/β水解酶的催化殘基都是由1個高度保守的三聯(lián)體組成。根據(jù)賈佳^[27]統(tǒng)計的EH家族的核心催化位點的特征，EHs除了上述三聯(lián)體結(jié)構(gòu)之外，在質(zhì)子供體、氧陰離子洞位置的氨基酸的選擇也具有高度保守性，且在β4結(jié)束的位置存在保守片段GXGXS（甘氨酸-可變氨基酸-甘氨酸-可變氨基酸-絲氨酸）。本研究中，我們鑒別的7個桃的EH家族候選成員均具有標準的α/β水解酶折疊結(jié)構(gòu)，與馬鈴薯、柑橘、煙草、擬南芥等其他植物中所報道的EH成員結(jié)構(gòu)的相似度高，且催化殘基包括2個天冬氨酸和1個組氨酸。在質(zhì)子供體的氨基酸選擇上，除了PpEH3之外，分析的其他EH成員位于α5螺旋的質(zhì)子供體均為酪氨酸。本研究所涉及的EH成員的氧陰離子洞位置的氨基酸的選擇均很保守，且具有GXGX（甘氨酸-可變氨基酸-甘氨酸-可變氨基酸）的保守片段，與賈佳^[27]的結(jié)論一致。

3.2 植物環(huán)氧化物水解酶的生物學功能

環(huán)氧化物水解酶在生物體系中發(fā)揮著重要的作用，包括外源化合物代謝、信號調(diào)節(jié)、細胞保護等，在哺乳動物中的研究相對深入，在植物中的研究略為滯后。Neuteboom等^[28]鑒別了在菠蘿根部呈現(xiàn)組織特異表達的基因，其中PFE258是長為1 199 bp的環(huán)氧化物酶基因。大戟屬植物（Euphorbia lagascae）的種子胚乳中含有高水平的環(huán)氧化脂肪酸-羥基乙酸，而環(huán)氧化物水解酶參與了種子萌發(fā)過程中羥基乙酸的氧化過程，Edqvist和Farbos^[29]從該植物中分離到了這種酶，并發(fā)現(xiàn)其編碼基因的表達受到種子萌發(fā)過程的誘導，而且外源激素如生長素和乙稀等的刺激也可增加該基因的轉(zhuǎn)錄。鄭柳城和朱宏波^[30]通過在水稻蛋白質(zhì)組數(shù)據(jù)庫中進行BLAST搜索，獲得了10個水溶性環(huán)氧化物水解酶，并對編碼這些酶的基因開展了表達模式分析、蛋白質(zhì)氨基酸序列聯(lián)配、三級結(jié)構(gòu)預測等生物信息學分析，預測了它們可能具有與抗逆相關(guān)的重要功能。在擬南芥中，近期有研究關(guān)注了一些EH的生物學功能，并證明EH參與了宿主防御的重要功能。AtEH1參與角質(zhì)的形成^[18]，催化水解C18脂肪酸環(huán)氧化物生成鄰二醇；在敲除了AtEH1基因的擬南芥突變體中，葉片和種子角質(zhì)中C18脂肪酸環(huán)氧化物含量會升高，并伴隨鄰二醇產(chǎn)物含量的降低。現(xiàn)有研究結(jié)果強調(diào)了植物EH在植物抗逆性方面的重要功能。

3.3 桃果實環(huán)氧化物水解酶基因的表達分析

目前涉及果實EH鑒別和生物學功能的研究較少，其中桃果實EH參與內(nèi)酯芳香物質(zhì)合成的相關(guān)研究備受關(guān)注。Pirona等^[20]報道了2個可能的EH成員，其中EPH3（編號：ppa009153m）的保守結(jié)構(gòu)域分析結(jié)果顯示其屬于HAD水解酶超家族，而EPH2（編號：ppa008854m）對應于本研究所鑒別的PpEH2。Pirona等^[20]對更為芳香的溶質(zhì)品種Bolero桃和香氣較弱的不溶質(zhì)品種OroA桃開展研究，發(fā)現(xiàn)EPH2在Bolero桃果實發(fā)育前期上調(diào)表達，在Bolero和OroA桃果實的發(fā)育后期均呈現(xiàn)下調(diào)表達。本研究所用品種白花水蜜為溶質(zhì)桃而且香氣濃郁，在該品種果實發(fā)育的S2階段PpEH2有明顯的上調(diào)表達，且在整個成熟期呈現(xiàn)下調(diào)表達，類似于Bolero桃的表達趨勢，與Pirona等^[20]的研究結(jié)果一致。Li等^[21]報道了若干可能的EH成員，其中EPH2（編號：ppa008854m）對應于本研究所鑒別的PpEH2，且與Pirona等^[20]所報道的EPH2是同一個基因；EPH1（編號：ppa008918m）對應于本研究所鑒別的PpEH6；但所報道的基因編號為ppa007004m和ppa013585m基因編碼的不屬于α/β水解酶結(jié)構(gòu)家族，另一個基因編號為ppa008756m基因編碼的不具有EH家族的保守蛋白質(zhì)結(jié)構(gòu)。本研究結(jié)果表明PpEH6和PpEH2在成熟期S4階段相比于S3階段表達量呈現(xiàn)下調(diào)趨勢，而Li等^[21]更精細地分析了這2個基因在HJ1階段（成熟早期）、HJ2階段（商業(yè)采收成熟期）、HJ3階段（完全成熟期）表達量的變化：EPH1（對應本研究的PpEH6）在湖景蜜露果實HJ2階段相比HJ1階段表達量呈現(xiàn)下調(diào)的趨勢，EPH2（對應本研究的PpEH2）在HJ2階段表達下調(diào)后在HJ3階段又呈現(xiàn)上調(diào)趨勢。

有研究結(jié)果表明，內(nèi)酯芳香物質(zhì)伴隨桃果實的成熟逐漸積累。Zhang等^[31]分析了內(nèi)酯物質(zhì)在湖景蜜露和奉化玉露桃果實的不同發(fā)育成熟階段（S0：未成熟階段；S1：早期成熟階段；S2：商業(yè)采收成熟階段；S3：完全成熟階段）的含量變化，發(fā)現(xiàn)γ-癸內(nèi)酯是其中含量最高的內(nèi)酯物質(zhì)，且該內(nèi)酯在2個品種桃果實中的含量均在S1階段呈現(xiàn)顯著增加，在S3階段達到最高。Peng等^[32]對盛花后114～132 d的Fenghuayulu桃果實每隔2 d進行1次取樣，并對芳香物質(zhì)進行檢測，發(fā)現(xiàn)從花后124 d開始可以在果實中檢測到γ-癸內(nèi)酯，且在隨后的階段中含量逐漸上升，在最后一個取樣點γ-癸內(nèi)酯含量達到最高。結(jié)合本研究的基因表達分析結(jié)果，在桃果實中表達的7個EH基因在成熟期的S3階段均有一定的表達量，可能與內(nèi)酯的合成相關(guān)聯(lián)。由于內(nèi)酯的合成通路尚未完全闡明，而環(huán)氧化物水解酶在內(nèi)酯合成通路的中游發(fā)揮作用，因此EH成員催化特征的鑒別、各EH成員催化產(chǎn)物的明確、EH產(chǎn)物是否以及如何作為前體物質(zhì)參與內(nèi)酯的合成等將是后續(xù)相關(guān)研究需要解決的關(guān)鍵問題。同時需要注意桃的EH是否具有多重生物學功能，如在其他物種中所報道的參與逆境防御的功能等。

綜上所述，本研究系統(tǒng)鑒別出7個桃的EH家族成員，序列特征分析結(jié)果表明它們具有典型的α/β水解酶折疊結(jié)構(gòu)，屬于環(huán)氧化物水解酶家族。進化樹分析結(jié)果顯示，桃的EH成員與擬南芥、煙草、油菜、大豆、粗皮檸檬等物種的EH的聚類很近?；虮磉_分析結(jié)果顯示，它們均在果實的S1～S4發(fā)育成熟階段表達，表達模式分為3種，基本呈現(xiàn)成熟前期表達量較高而在整個成熟期整體表達量下調(diào)的總趨勢，與內(nèi)酯芳香物質(zhì)的積累負相關(guān)。本研究對于桃環(huán)氧化物水解酶家族成員的構(gòu)成、進化關(guān)系和基因表達模式進行了系統(tǒng)分析，為后續(xù)桃乃至其他果實中內(nèi)酯物質(zhì)合成相關(guān)的EH成員的深入研究提供參考。

參考文獻：

[1] 王貴章，王貴禧，梁麗松，等. 桃果實芳香揮發(fā)物及其生物合成研究進展 [J].食品科學，2014，35（17）：278-284.

[2] SCH?TTLER M，BOLAND W. Biosynthesis of dodecano-4-lactone in ripening fruits： crucial role of an epoxide-hydrolase in enantioselective generation of aroma components of the nectarine （Prunus persica var. nucipersica） and the strawberry （Fragaria ananassa） [J]. Helvetica Chimica Acta，1996，79：1488-1496.

[3] ZOU J，HALLBERG B M，BERGFORS T，et al. Structure of Aspergillus niger epoxide hydrolase at 1.8 A resolution： implications for the structure and function of the mammalian microsomal class of epoxide hydrolases [J]. Structure，2000，8（2）：111-122.

[4] 盛艷旻. 環(huán)氧化物水解酶產(chǎn)生菌的篩選及其發(fā)酵條件和對苯基縮水甘油醚類化合物拆分條件的研究 [D]. 杭州：浙江工業(yè)大學，2010.

[5] 陳文靜. 綠豆環(huán)氧化物水解酶催化環(huán)氧苯乙烯不對稱水解反應 [D]. 廣州：華南理工大學，2012.

[6] NEWMAN J W，MORISSEAU C，HAMMOCK B D. Epoxide hydrolases： their roles and interactions with lipid metabolism [J]. Progress in Lipid Research，2005，44（1） ：1-51.

[7] MORISSEAU C. Role of epoxide hydrolase in lipid metabolism [J]. Biochimie，2013，95（1）：91-95.

[8] 婁文勇，趙 瑩，彭 飛，等. 環(huán)氧化物水解酶的研究進展 [J]. 華南師范大學學報（自然科學版），2017，49（6）：1-6.

[9] BLE E，SCHUBER F. Occurrence of fatty acid epoxide hydrolases in soybean （Glycine max）. Purification and characterization of the soluble form [J]. Biochemical Journal，1992，282（3）：711-714.

[10]KIYOSUE T，BEETHAM J K，PINOT F，et al. Characterization of an Arabidopsis cDNA for a soluble epoxide hydrolase gene that is inducible by auxin and water stress [J]. The Plant Journal，1994，6（2）：259-269.

[11]MORISSEAU C，BEETHAM J K，PINOT F，et al. Cress and potato soluble epoxide hydrolases： purification， biochemical characterization， and comparison to mammalian enzymes [J]. Archives of Biochemistry and Biophysics，2000，378（2）：321-332.

[12]SUMMERER S，HANANO A，UTSUMI S，et al. Stereochemical features of the hydrolysis of 9，10-epoxystearic acid catalysed by plant and mammalian epoxide hydrolases [J]. Biochemical Journal，2002，366（2）：471-480.

[13]MOWBRAY S L，ELFSTR?M L T，AHLGREN K M，et al. X-ray structure of potato epoxide hydrolase sheds light on substrate specificity in plant enzymes [J]. Protein Science，2006，15（7）：1628-1637.

[14]THOMAEUS A，NAWORYTA A，MOWBRAY S L，et al. Removal of distal protein-water hydrogen bonds in a plant epoxide hydrolase increases catalytic turnover but decreases thermostability [J]. Protein Science，2008，17：1275-1284.

[15]WIJEKOON C P，GOODWIN P H，VALLIANI M，et al. The role of a putative peroxisomal-targeted epoxide hydrolase of Nicotiana benthamiana in interactions with Colletotrichum destructivum， C. orbiculare or Pseudomonas syringae pv. tabaci [J]. Plant Science，2011，181（2）：177-187.

[16]HUANG F C，SCHWAB W. Molecular characterization of NbEH1 and NbEH2， two epoxide hydrolases from Nicotiana benthamiana [J]. Phytochemistry，2013，90：6-15.

[17]CARLSSONJ，BAUER P，DOBRITZSCH D，et al. Laboratory evolved enzymes provide snapshots of the development of enantioconvergence in enzyme-catalyzed epoxide hydrolysis [J]. Chembiochem，2016，17（18）：1693-1697.

[18]PINEAU E，XU L，RENAULT H，et al. Arabidopsis thaliana epoxide hydrolase1 （AtEH1） is a cytosolic epoxide hydrolase involved in the synthesis of poly-hydroxylated cutin monomers [J]. New Phytologist，2017，215：173-186.

[19]VECCHIETTI A，LAZZARI B，ORTUGNO C，et al. Comparative analysis of expressed sequence tags from tissues in ripening stages of peach （Prunus persica L. Batsch） [J]. Tree Genetics & Genomes，2009，5（3）：377-391.

[20]PIRONA R，VECCHIETTI A，LAZZARI B，et al. Expression profiling of genes involved in the formation of aroma in two peach genotypes [J]. Plant Biology，2013，15（3）：443-451.

[21]LI X W，JIANG J，ZHANG L P，et al. Identification of volatile and softening-related genes using digital gene expression profiles in melting peach [J]. Tree Genetics & Genomes，2015，11：1-15.

[22]ROBERT X，GOUET P. Deciphering key features in protein structures with the new ENDscript server [J]. Nucleic Acids Research，2014，42（1）：320-324.

[23]ZHANG Y Y，YIN X R，XIAO Y W，et al. An ethylene response factor-MYB transcription complex regulates furaneol biosynthesis by activating quinone oxidoreductase expression in strawberry [J]. Plant Physiology，2018，178（1）：189-201.

[24]TONG Z G，GAO Z H，WANG F，et al. Selection of reliable reference genes for gene expression studies in peach using real-time PCR [J]. BMC Molecular Biology，2009，10：1-13.

[25]OLLIS D L，CHEAH E，CYGLER M，et al. The α/β hydrolase fold [J]. Protein Engineering Design & Selection，1992，5（3）：197-211.

[26]徐麗詩. α/β水解酶亞家族共進化網(wǎng)絡的算法優(yōu)化及功能預測 [D]. 上海：上海交通大學，2016.

[27]賈 佳. 環(huán)氧化物水解酶序列識別和分類模型的構(gòu)建 [D]. 杭州：浙江大學，2006.

[28]NEUTEBOOM L W，KUNIMITSU W Y，WEBB D，et al. Characterization and tissue-regulated expression of genes involved in pineapple （Ananas comosus L.） root development [J]. Plant Science，2002，163（5）：1021-1035.

[29]EDQVIST J，F(xiàn)ARBOS I. A germination-specific epoxide hydrolase from Euphorbia lagascae [J]. Planta，2003，216：403-412.

[30]鄭柳城，朱宏波. 水稻水溶性環(huán)氧化合物水解酶的生物信息學分析 [J]. 生物信息學，2009，7（2）：108-112， 139.

[31]ZHANG L P，LI H Y，GAO L，et al. Acyl-CoA oxidase 1 is involved in γ-decalactone release from peach （Prunus persica） fruit [J]. Plant Cell Reports，2017，36（6）：829-842.

[32]PENG B，YU M L，ZHANG B B，et al. Differences in PpAAT1 Activity in high- and low-aroma peach varieties affect γ-decalactone production [J]. Plant Physiology，2020，182（4）：2065-2080.

（責任編輯：陳海霞）