獼猴桃FAD基因家族鑒定及其在采后成熟過程中的表達(dá)分析

摘" " 要：【目的】分離并鑒定獼猴桃FAD基因家族成員，探究FAD基因在獼猴桃采后成熟過程中的表達(dá)模式及其與香氣物質(zhì)合成前體不飽和脂肪酸變化的關(guān)系。【方法】利用生物信息學(xué)方法鑒定并分析獼猴桃FAD基因家族成員，采用質(zhì)構(gòu)儀和氣相色譜儀測(cè)定硬度和脂肪酸含量，借助實(shí)時(shí)熒光定量PCR驗(yàn)證FAD在采后成熟過程中的表達(dá)特性?！窘Y(jié)果】在紅陽獼猴桃基因組中共鑒定出了26個(gè)FAD基因，分為6個(gè)亞族；該家族大多為堿性蛋白，大部分定位于內(nèi)質(zhì)網(wǎng)中；這些成員不均勻地分布在19條染色體上，種內(nèi)共有9對(duì)串聯(lián)重復(fù)基因和22對(duì)片段重復(fù)基因；在其啟動(dòng)子區(qū)域上發(fā)現(xiàn)大量的光響應(yīng)元件、植物激素響應(yīng)元件、逆境脅迫響應(yīng)元件和生長(zhǎng)發(fā)育相關(guān)元件；表達(dá)模式分析和qPCR驗(yàn)證實(shí)驗(yàn)表明，AcFAD2.2表達(dá)量均隨著成熟而不斷顯著上調(diào)；采后成熟過程中單不飽和脂肪酸（油酸）含量不斷顯著下降，雙不飽和脂肪酸（亞油酸）含量則不斷顯著提高，多不飽和脂肪酸（亞麻酸）含量在成熟早期和中期無顯著差別，而在后期顯著降低；獼猴桃硬度采后出現(xiàn)快速下降，存在明顯的后熟過程，此過程會(huì)形成特征香氣酯類物質(zhì)，而亞油酸又是酯類香氣物質(zhì)合成的主要前體物質(zhì)?！窘Y(jié)論】共鑒定了26個(gè)獼猴桃FAD基因成員，并篩選出了1個(gè)誘導(dǎo)獼猴桃采后成熟過程中亞油酸合成和積累的關(guān)鍵酶基因AcFAD2.2，亞油酸含量的增加伴隨著獼猴桃后熟散發(fā)特征性香氣，為進(jìn)一步探究FAD基因參與獼猴桃采后成熟過程不飽和脂肪酸轉(zhuǎn)化和香氣合成的生物學(xué)功能提供了參考依據(jù)。

關(guān)鍵詞：獼猴桃；FAD基因家族；不飽和脂肪酸；香氣合成；基因表達(dá)

中圖分類號(hào)：S663.4 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1009-9980（2024）11-2195-19

Identification of the FAD gene family in kiwifruit and its expression analysis during postharvest ripening

YANG Caining， ZHANG Yupei， YANG Congcong， CHEN Jinyin*， GAN Zengyu*

（Jiangxi Provincial Key Laboratory for Postharvest Storage and Preservation of Fruits amp; Vegetables/Jiangxi Agricultural University， Nanchang 330045， Jiangxi， China）

Abstract： 【Objective】 Unsaturated fatty acids play many physiological roles in plants， including the formation of triglycerides to provide energy for plant life activities; maintaining cell membrane homeostasis as a key component of cell membranes; participating in hormone regulation and signaling during biotic stress; playing an important role in responding to adversity such as low temperature and drought; and participating in synthesizing aroma substances as precursors during fruit ripening via the LOX pathway， α-oxidation pathway， or β-oxidation pathway. During fruit ripening， the unsaturated fatty acids are used as precursors to synthesize aroma substances through the LOX pathway， α-oxidation pathway， or β-oxidation pathway. At the same time， the unsaturated fatty acids also help the human body in numerous ways， including lowering lipids and blood pressure， enhancing fat metabolism， decreasing thrombosis， enhancing autoimmunity， and having anti-tumor properties. However， the human body is unable to synthesize unsaturated fatty acids， such as linoleic acid and linolenic acid， which can only be ingested from the diet. The animal body contains a high concentration of saturated fatty acids， whereas the plant body contains primarily unsaturated fatty acids. The fatty acid desaturases （FADs） are a group of enzymes that carry out the desaturation process， transforming saturated fatty acids into unsaturated fatty acids. To date， the FAD gene family has been identified in several species; however， it has not been reported in kiwifruit. In addition， most of the current studies on FAD genes have focused on the antistress function of unsaturated fatty acids， with little focus on how these fatty acids affect aroma volatiles. Therefore， the aim of this study was to isolate and characterize the kiwifruit FAD gene family and to clarify its expression pattern in kiwifruit during postharvest ripening and its relationship with the changes in unsaturated fatty acids， so as to lay a theoretical foundation for analyzing the formation of characteristic aroma in kiwifruit during postharvest ripening.【Methods】 Based on the Hong Yang v3 genomic data of Actinidia chinensis， we downloaded the Hidden Markov Model （HMM） file corresponding to the structural domains of the fatty acid desaturase from the Pfam protein family database and used Simple HMM Search in TBtools （v2.102） to preliminarily search the AcFAD genes with the kiwifruit protein data， and then used the SMART database to verify the structural domain information of the candidate protein sequences to finalize the kiwifruit FAD gene family members. Then we used the SMART database to validate the structural domain information of the initial screened candidate protein sequences to identify the kiwifruit FAD gene family members; we used the protein molecular weight calculation （SMS2） Nanjing DETA Bio-mirror website to analyze the physicochemical properties of the kiwifruit FAD gene family members; and we used the Cell-PLoc 2.0 website to predict the subcellular localization of the proteins; A phylogenetic evolutionary tree of FAD protein sequences of kiwifruit， Arabidopsis thaliana and cucumber was constructed using MEGA （version 11.0） software; the chromosomal localization was mapped using Advanced Circos in TBtools， and the replication events of kiwifruit and Arabidopsis thaliana FAD genes were analyzed by One Step MCScanX-Super Fast. The covariance analysis was performed; the exon-intron gene structures were analyzed by GSDS， a gene structure display server; the conserved motifs of AcFAD proteins were analyzed using the MEME online website; the promoter cis-acting elements of each gene were predicted using PlantCARE; and the gene expression heatmaps were generated using TBtools by normalizing the FPKM values of the transcriptomic data obtained. The expression pattern of the FAD genes was explored; the hardness and fatty acid content were determined using mass spectrometry and gas chromatography; and the expression characteristics of the FAD genes were verified with the help of the real-time fluorescence quantitative PCR during post-harvest ripening. 【Results】 In this study， a total of 26 FAD gene family members were identified from the whole Hong Yang v3 genome of Actinidia chinensis. These 26 kiwifruit FAD genes were named based on how closely they resembled the counterparts identified in Arabidopsis thaliana and cucumber. These genes were divided into six subfamilies： FAD3/FAD7/FAD8 （ω-3/Δ-15）， FAD2/FAD6 （ω-6/Δ-12）， FAB2 （Δ-9）， FAD4 （Δ-3）， DES/SLD， and FAD5/ADS （Δ-7）， each with a different number of subfamily members. All subfamilies had kiwifruit FAD family members， suggesting that the AcFAD proteins might be functionally diversified. A handful of them were acidic， but the majority were basic. The subcellular localization prediction revealed that kiwifruit FAD proteins were relatively dispersed in their localization and are distributed in all plant cell structures， with the highest number in the endoplasmic reticulum， followed by chloroplasts and cell membranes; the chromosome localization showed that the kiwifruit Hong Yang v3 genome had a total of 29 chromosomes， and 25 kiwifruit FAD genes were distributed on 19 different chromosomes （LG）. The majority of the FAD genes were localized in the anterior-middle end of the chromosomes. The segmental duplication events contributed to the diversity and evolution of the AcFADs， as evidenced by the 22 pairs of segmental duplicates and 9 pairs of tandem duplicates found in kiwifruit according to the intraspecific covariance mapping; the covariance mapping between kiwifruit and Arabidopsis thaliana interspecies covariance mapping showed the 32 pairs of covariance between the 23 kiwifruit FAD genes and the 17 Arabidopsis thaliana FAD genes， suggesting that kiwifruit and Arabidopsis thaliana would have more FAD homologous genes. In the same branch， most of the members had similar length and the same number of gene structures. The conserved motif analysis revealed that most kiwifruit FAD proteins contained motif 2， motif 4， and motif 14. Moreover， motif 2 and motif 4 were basically located at the C-terminal end， indicating that these motifs were strongly conserved and were typical FAD structural domains， which might perform similar functions. A large number of the light-responsive， phytohormone-responsive （methyl jasmonate， abscisic acid， growth hormone， salicylic acid， gibberellin）， stress-responsive （anaerobic， low temperature， drought）， and development-related components were found in the FAD promoter region of kiwifruit， and the heatmap expression pattern analysis and qPCR validation experiments showed that the expression of the AcFAD7.2 and the AcFAD2.2 were significantly up-regulated with maturation. During the postharvest ripening process， the content of monounsaturated fatty acids （oleic acid） decreased significantly， while the content of diunsaturated fatty acids （linoleic acid） increased significantly， and the content of polyunsaturated fatty acids （linolenic acid） did not significantly vary with the early and middle stages of ripening， while the content of polyunsaturated fatty acids （linolenic acid） decreased significantly during the late stage; the hardness of kiwifruit rapidly decreased after harvest， and a distinct post-ripening process took place that resulted in the formation of the distinctive aromatic ester substances. The characteristic fragrance esters were produced as a result of this process， and the primary precursor for the synthesis of the ester aroma compounds was linoleic acid. 【Conclusion】 A total of 26 kiwifruit FAD gene members were identified. A significant decrease in the oleic acid content was accompanied by a significant increase in the linoleic acid content during postharvest ripening of kiwifruit. The variations in the oleic and linoleic acid were highly and positively linked with the expression of the ω-6 fatty acid desaturase gene AcFAD2.2. This expression was consistently and considerably up-regulated throughout ripening. Although the expression of the ω-3 fatty acid desaturase gene AcFAD7.2 was up-regulated， the linolenic acid content remained essentially unchanged or even decreased. Therefore， the AcFAD2.2 would be a key enzyme gene involved in the synthesis and accumulation of unsaturated fatty acids in kiwifruit during postharvest ripening， and the increase in the linoleic acid content was accompanied by the characteristic aroma emitted by kiwifruit during postharvest ripening. This paper would provide a reference basis for further investigation of the biological functions of the FAD genes involved in the transformation of unsaturated fatty acids and aroma synthesis during the course of the postharvest ripening in kiwifruit.

Key words： Kiwifruit; FAD gene family; Unsaturated fatty acids; Aroma synthesis; Gene expression

脂肪酸是植物細(xì)胞的重要組成成分，它形成甘油三酯為生命活動(dòng)提供能量[1]；作為細(xì)胞膜的關(guān)鍵成分維持細(xì)胞膜的穩(wěn)態(tài)[2]；在生物脅迫中參與激素調(diào)節(jié)與信號(hào)轉(zhuǎn)導(dǎo)[3]；在響應(yīng)低溫和干旱等逆境脅迫中發(fā)揮重要作用[4-5]；在果實(shí)成熟過程中通過LOX途徑、α-氧化途徑或β-氧化途徑參與合成香氣物質(zhì)[6]。植物中以不飽和脂肪酸為主，飽和脂肪酸轉(zhuǎn)化成不飽和脂肪酸的過程稱為去飽和，是由一系列脂肪酸去飽和酶（fatty acid desaturases，F(xiàn)ADs）來完成的[7]。FADs根據(jù)溶解度可分為可溶性去飽和酶和膜結(jié)合去飽和酶。硬脂酰ACP去飽和酶（FAB2/SAD）是質(zhì)體基質(zhì)中唯一已知的可溶性FAD，它在Δ-9位置（從羧基端開始的第9號(hào)碳與10號(hào)碳之間）引入雙鍵，催化硬脂酸轉(zhuǎn)化為油酸，剩余的FAD均是膜結(jié)合的[8]。根據(jù)功能不同，膜結(jié)合FADs又進(jìn)一步分為FAD2/FAD6（ω-6/Δ-12）、FAD3/FAD7/FAD8 （ω-3/Δ-15）、FAD4（Δ-3）、FAD5/ADS（Δ-7）和DES/SLD五大亞族。其中，F(xiàn)AD2和FAD6是ω-6/Δ-12去飽和酶，即FAD2在內(nèi)質(zhì)網(wǎng)、FAD6在質(zhì)體中于油酸的Δ-12位置，也就是ω-6位置（甲基端開始的第6號(hào)碳與7號(hào)碳之間）插入雙鍵生成亞油酸[9]；FAD3在內(nèi)質(zhì)網(wǎng)、FAD7和FAD8在質(zhì)體中于亞油酸的Δ-15位置（ω-3位置）插入雙鍵生成亞麻酸[10]；FAD4和FAD5/ADS都在質(zhì)體中分別作用于磷脂酰甘油和單半乳糖二?；视椭械淖貦八岷铣勺貦坝退醄11]；DES和SLD都是在內(nèi)質(zhì)網(wǎng)中參與脂肪酸衍生物鞘脂代謝的酶，SLD催化t18∶0鞘脂長(zhǎng)鏈堿基（LCB）C8位去飽和生成t18∶1鞘脂[12]，DES催化二氫鞘氨醇（d18：0）生成鞘氨醇-1-磷酸（d18：1），DES通常不會(huì)單獨(dú)起作用，而是和SLD共同參與鞘脂代謝[13]。

1992年有學(xué)者先后從擬南芥突變體中克隆得到FAD3和FAD2 [14-15]，1993年和1994年從擬南芥中先后克隆出FAD7和FAD8 [16]，1993年從藍(lán)藻中克隆出FAD6 [17]，SLD最先于1998年從向日葵中分離得到[12]。目前研究比較多的是ω-6和ω-3型FAD基因，已經(jīng)在香蕉[18]、茄子[19]、番茄[20]等中分離出其家族成員。當(dāng)前對(duì)FAD基因研究熱點(diǎn)主要集中在抗逆功能，尤其是低溫冷害脅迫，大量研究證實(shí)了脂肪酸去飽和酶能夠調(diào)整膜脂不飽和脂肪酸的比例和成分，從而維持細(xì)胞膜流動(dòng)性，達(dá)到增強(qiáng)對(duì)低溫抵抗力的目的[21]。然而，關(guān)于FAD基因在果實(shí)成熟中的功能鮮有報(bào)道。

果實(shí)中的揮發(fā)性芳香物質(zhì)主要通過脂肪酸代謝途徑、萜類代謝途徑和氨基酸代謝途徑合成。其中，直鏈脂肪族醇、醛、酮和酯類物質(zhì)主要來源于脂肪酸的代謝，因此脂肪酸是形成果實(shí)香氣物質(zhì)的主要前體物質(zhì)[22]。脂肪酸代謝途徑中對(duì)香氣物質(zhì)合成貢獻(xiàn)最大的是不飽和脂肪酸的代謝，這是通過LOX酶催化的第一步和隨后的反應(yīng)，通常被稱為脂氧合酶途徑。脂肪酸去飽和酶FAD在脂肪酸碳?xì)滏溕弦腚p鍵，產(chǎn)生不飽和脂肪酸，在將油酸（18∶1）轉(zhuǎn)化為亞油酸（18∶2）和亞麻酸（18∶3）過程中發(fā)揮著關(guān)鍵的作用[23]。

獼猴桃果實(shí)中因富含有機(jī)物質(zhì)以及人體所必需的多種維生素而備受消費(fèi)者青睞。獼猴桃在發(fā)育階段（未成熟階段）主要產(chǎn)生C6醇和醛等“青香型”芳香物質(zhì)，但在采收后熟過程中酯類等“果香型”芳香物質(zhì)增多，因此常會(huì)散發(fā)出成熟階段特有的果香味，也影響著消費(fèi)者的偏好[24]。目前關(guān)于獼猴桃采后成熟過程中的香氣研究較多，主要集中在貯藏處理方式、成熟度和品種等因素對(duì)其香氣物質(zhì)化學(xué)成分種類和含量比例的影響方面[25-26]，香氣脂肪酸代謝合成途徑的相關(guān)酶基因（如LOX、AAT、ADH等）也已被分離鑒定出來[27-29]。然而對(duì)其合成前體物質(zhì)不飽和脂肪酸的研究鮮有報(bào)道，目前發(fā)現(xiàn)低溫貯藏延緩了獼猴桃果實(shí)中亞油酸和亞麻酸的分解，抑制了脂肪酸LOX代謝途徑中相關(guān)酶活性及其基因的表達(dá)，從而造成酯類香氣物質(zhì)的種類減少和相對(duì)含量降低，保持了較多醛酮類物質(zhì)種類和較高的相對(duì)含量，維持了果實(shí)的特征風(fēng)味[30]，其他相關(guān)的研究更多集中在獼猴桃籽和果肉中脂肪酸組分方面，發(fā)現(xiàn)其富含多種不飽和脂肪酸[31-32]，獼猴桃FAD基因家族也未進(jìn)行分離鑒定。紅陽獼猴桃屬于中華獼猴桃，因其果心具放射狀紅色條紋而得名，味甜，維生素C含量極高，香氣濃郁[33]。因此，筆者以紅陽獼猴桃為研究試材，采用生物信息學(xué)的方法對(duì)獼猴桃FAD基因家族成員進(jìn)行分離和鑒定，然后再基于轉(zhuǎn)錄組數(shù)據(jù)對(duì)各成員在果實(shí)采后成熟過程中的表達(dá)模式進(jìn)行分析，最后通過實(shí)時(shí)熒光定量PCR驗(yàn)證FADs基因的表達(dá)特性并篩選出差異表達(dá)的候選基因，為解析FAD基因在獼猴桃果實(shí)采后成熟過程中不飽和脂肪酸代謝和香氣物質(zhì)形成的生物學(xué)功能提供一定的理論基礎(chǔ)。

1 材料和方法

1.1 植物材料

材料為紅陽獼猴桃（Actinidia chinensis ‘Hongyang’）果實(shí)。于2023年8月中旬（盛花后130 d左右）從江西省奉新縣果業(yè)局獼猴桃種質(zhì)資源圃采摘試驗(yàn)果。選取樹齡相同、長(zhǎng)勢(shì)良好的母樹，采摘大小均勻、果形一致、無病蟲害和機(jī)械損傷的果實(shí)，采摘當(dāng)天立即運(yùn)回實(shí)驗(yàn)室，經(jīng)人工挑選剔除癍痂、傷果后置于陰涼通風(fēng)處發(fā)汗到次日。經(jīng)發(fā)汗散除田間熱的果實(shí)用塑料袋密封貯藏于室溫（20±1）℃中，于貯藏后0、4、10 d取樣，每次隨機(jī)挑選18個(gè)果實(shí)，6個(gè)果實(shí)為1次重復(fù)，設(shè)置3次生物學(xué)重復(fù)，經(jīng)測(cè)量果肉硬度等成熟指標(biāo)后去掉頭部和尾部，剝除果皮、中柱和籽，取果肉部分混合均勻經(jīng)液氮速凍后保存于-80 ℃?zhèn)溆谩?/p>

1.2 獼猴桃FAD基因家族成員鑒定及其蛋白特性分析

從獼猴桃基因組網(wǎng)站（https：//kiwifruitgenome.org/）下載Hong Yang v3基因組數(shù)據(jù)，包括基因序列數(shù)據(jù)文件、蛋白序列數(shù)據(jù)文件、GFF3注釋文件等。從Pfam蛋白家族數(shù)據(jù)庫（http：//pfam.xfam.org/）下載FA_去飽和酶（PF00487）、FA_去飽和酶2（PF03405）和TMEM189 （PF10520）結(jié)構(gòu)域?qū)?yīng)的隱馬爾科夫模型（HMM）文件，使用TBtools（v2.102）中的Simple HMM Search，以e值≤1e-5為判據(jù)，與獼猴桃蛋白數(shù)據(jù)進(jìn)行比對(duì)，初步篩選出AcFAD基因。然后使用SMART數(shù)據(jù)庫（https：//smart.embl-heidelberg.de/）對(duì)初篩得到的候選蛋白序列進(jìn)行結(jié)構(gòu)域信息驗(yàn)證，最終確定了26個(gè)獼猴桃FAD基因家族成員。在獼猴桃基因組網(wǎng)站中獲得這26個(gè)AcFADs編碼區(qū)長(zhǎng)度。利用蛋白分子質(zhì)量計(jì)算-SMS2南京德泰生物鏡像網(wǎng)站（https：//www.detaibio.com/sms2/protein_mw.html）對(duì)26個(gè)獼猴桃FAD基因家族成員進(jìn)行蛋白質(zhì)理化性質(zhì)分析，獲得氨基酸數(shù)量、蛋白分子質(zhì)量（MW）、等電點(diǎn)（pI）等。通過Cell-PLoc 2.0網(wǎng)站（http：//www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/）中的Plant-mPLoc對(duì)AcFAD蛋白進(jìn)行亞細(xì)胞定位預(yù)測(cè)。

1.3 獼猴桃FAD基因家族多序列比對(duì)與系統(tǒng)發(fā)育分析

擬南芥和黃瓜的FAD蛋白質(zhì)序列來自Cheng等[18]和Dong等[34]的報(bào)道，利用MEGA（version 11.0）軟件中的ClustalW對(duì)獼猴桃、擬南芥和黃瓜的FAD蛋白序列進(jìn)行多序列比對(duì)。系統(tǒng)發(fā)育分析采用鄰接法（neighbor-joining，NJ），在p-distance、部分刪除50%和1000次Bootstrap重復(fù)參數(shù)下構(gòu)建系統(tǒng)發(fā)育樹，并使用Evolview v2（https：//evolgenius.info//evolview-v2/）對(duì)進(jìn)化樹進(jìn)行美化。

1.4 獼猴桃FAD基因家族染色體定位和共線性分析

根據(jù)獼猴桃和擬南芥的全基因組序列文件和GFF3注釋文件，使用TBtools中的Fasta Stats提取獼猴桃全基因組染色體長(zhǎng)度并輸出染色體骨架文件，利用One Step MCScanX獲得全基因組共線性文件并據(jù)此文件使用Text Block Extract and Filter提取獼猴桃FAD家族基因共線性文件，使用File Transformat for MicroSynteny Viewer提取配對(duì)基因?qū)θ旧w信息，最后使用Advanced Circos繪制獼猴桃FAD染色體定位圖。使用One Step MCScanX-Super Fast對(duì)獼猴桃和擬南芥FAD基因的復(fù)制事件進(jìn)行共線性分析，并使用Dual Systeny Plot作可視化圖。

1.5 獼猴桃FAD基因結(jié)構(gòu)與蛋白保守基序分析

通過基因結(jié)構(gòu)顯示服務(wù)器GSDS（http：//gsds.gao-lab.org/）分析獼猴桃FAD基因的外顯子-內(nèi)含子基因結(jié)構(gòu)。使用MEME在線網(wǎng)站（https：//meme-suite.org/meme/tools/meme）分析AcFAD蛋白的保守基序，最大基序數(shù)設(shè)置為20，并利用TBtools中的Gene Structure View繪制蛋白保守基序。

1.6 獼猴桃FAD基因家族啟動(dòng)子順式作用元件分析

使用TBtools從獼猴桃基因組數(shù)據(jù)庫中提取每個(gè)AcFAD基因起始密碼子（ATG）上游1500 bp啟動(dòng)子序列，利用PlantCARE（http：//bioinformatics.psb.ugent.be/webtools/plantcare-/html/）預(yù)測(cè)各個(gè)基因的啟動(dòng)子順式作用元件，再利用TBtools對(duì)獲得的順式作用元件進(jìn)行可視化繪圖。

1.7 AcFADs在果實(shí)采后成熟過程中的基因表達(dá)分析

將常溫貯藏0、4、10 d獼猴桃果實(shí)冷凍樣品分裝3管設(shè)置3次生物學(xué)重復(fù)，在上海派森諾生物科技有限公司完成轉(zhuǎn)錄組測(cè)序。將獲得的轉(zhuǎn)錄組數(shù)據(jù)提取上述26個(gè)AcFAD基因的FPKM值，將log2（FPKM+1）歸一化處理后得到的數(shù)據(jù)使用TBtools生成基因表達(dá)熱圖。

1.8 果實(shí)硬度的測(cè)量

使用質(zhì)構(gòu)儀（TMS-Touch，美國(guó)）測(cè)定果實(shí)的硬度，每個(gè)取樣點(diǎn)隨機(jī)選取18個(gè)果實(shí)，重復(fù)測(cè)量3次，單位：N[35]。

1.9 果實(shí)脂肪酸含量的測(cè)定

將獼猴桃果肉冷凍樣品粉末2 g與正己烷∶異丙醇（3∶2，體積比）15 mL和6.7% Na2SO4 7.5 mL混合，10 000g離心10 min（4 ℃），上清液用氮?dú)獯蹈?。加入甲醇∶甲苯∶H2SO4（88∶10∶2，體積比）制備脂肪酸甲酯（FAMEs）。冷卻后，加入1 mL庚烷和0.5 g無水Na2SO4進(jìn)行FAME提取。為了檢測(cè)脂肪酸含量，使用安捷倫6890 N氣相色譜儀，配備火焰電離檢測(cè)器和DB-Wax色譜柱（0.25 mm，30 m，0.25 μm；Jamp;W Scientific）。注入器和檢測(cè)器溫度為230 ℃。初始烤箱溫度為50 ℃，以25 ℃·min-1速率增加到200 ℃，然后以3 ℃·min -1速率增加到230 ℃。載氣為氮?dú)猓魉贋? mL·min-1。加入外源十七烷酸（C17∶0）作為內(nèi)標(biāo)定量計(jì)算脂肪酸含量。

1.10 RNA提取、cDNA合成及實(shí)時(shí)熒光定量PCR分析

使用多糖多酚植物總RNA提取試劑盒（天根，北京）提取獼猴桃果肉總RNA，每個(gè)樣品設(shè)置3次技術(shù)重復(fù)，然后通過凝膠電泳和超微量分光光度計(jì)（Bio-Rad，美國(guó)）測(cè)定提取的R N A的質(zhì)量和濃度。將提取好的RNA使用三代逆轉(zhuǎn)錄預(yù)混液MonScriptTM RTIII All-in-One Mix with dsDNase（莫納生物，蘇州）反轉(zhuǎn)錄成第一鏈cDNA，并用ddH2O稀釋3倍，最后采用半定量PCR檢測(cè)cDNA的質(zhì)量。通過WoLF PSORT網(wǎng)站（https：//www.genscript.com/wolf-psort.html）設(shè)計(jì)實(shí)時(shí)熒光定量引物（表1）。實(shí)時(shí)熒光定量PCR（qPCR）分析參照之前的研究方法[36]。

1.11 數(shù)據(jù)統(tǒng)計(jì)與分析

所有數(shù)據(jù)結(jié)果以均數(shù)±標(biāo)準(zhǔn)誤差（SD）表示。采用SPSS Statistics 25軟件中的單因素方差檢驗(yàn)進(jìn)行顯著性差異分析并使用Duncan新復(fù)極差字母標(biāo)記法標(biāo)注顯著性，最后使用 GraphPad Prism 8.0繪制圖表。

2 結(jié)果與分析

2.1 獼猴桃FAD基因家族成員鑒定及其蛋白特性分析

筆者從Hong Yang v3全基因組中共鑒定出26個(gè)FAD基因家族成員（表2），根據(jù)其與擬南芥和黃瓜中FAD基因親緣關(guān)系的遠(yuǎn)近命名，包括2個(gè)AcFAD2s，2個(gè)AcFAD6s，3個(gè)AcFAD3s，3個(gè)AcFAD7s，1個(gè)AcFAB2，1個(gè)AcFAD4，3個(gè)AcFAD5s，2個(gè)AcADSs，7個(gè)AcSLDs和2個(gè)AcDESs。它們的編碼序列（CDS）長(zhǎng)度差異較大，范圍為321 bp（AcFAD5.1）～4392 bp（AcFAD7.3），編碼106～1463個(gè)氨基酸。蛋白質(zhì)理論分子質(zhì)量（MW）為11.51～161.8 ku，等電點(diǎn)為4.93（AcFAD5.1）～10.09（AcADS1.1），其中只有AcFAD7.3、AcFAD5.1和AcSLD1.4這3個(gè)蛋白的pI小于7，說明獼猴桃FAD蛋白大多數(shù)為堿性，少數(shù)為酸性。亞細(xì)胞定位預(yù)測(cè)結(jié)果表明，獼猴桃FAD蛋白定位比較分散，在植物細(xì)胞各結(jié)構(gòu)中均有分布。其中AcFAD2蛋白都只定位在內(nèi)質(zhì)網(wǎng)中；AcFAD6蛋白只定位在葉綠體中；AcFAD3、AcFAD7和AcADS蛋白在葉綠體和內(nèi)質(zhì)網(wǎng)中都有定位；AcFAB2定位在細(xì)胞膜和細(xì)胞核中；AcFAD4定位在細(xì)胞膜、葉綠體、線粒體和細(xì)胞核中；AcFAD5.1定位在細(xì)胞膜、葉綠體、高爾基體和細(xì)胞核中，另2個(gè)AcFAD只定位在內(nèi)質(zhì)網(wǎng)中；AcSLD1.1定位在細(xì)胞膜和內(nèi)質(zhì)網(wǎng)中，其他AcSLD只定位在內(nèi)質(zhì)網(wǎng)中；AcDES1.2定位在細(xì)胞膜和內(nèi)質(zhì)網(wǎng)中，AcDES1.1只定位在內(nèi)質(zhì)網(wǎng)中。

2.2 獼猴桃FAD基因的系統(tǒng)進(jìn)化樹分析

為了闡明AcFAD蛋白的功能和進(jìn)化關(guān)系，將26個(gè)獼猴桃、19個(gè)擬南芥和22個(gè)黃瓜的FAD蛋白序列構(gòu)建了系統(tǒng)發(fā)育樹。這些FAD可分為6個(gè)亞族，包括FAD3/FAD7/FAD8（ω-3/Δ-15）、FAD2/FAD6 （ω-6/Δ-12）、FAB2（Δ-9）、FAD4（Δ-3）、DES/SLD和FAD5/ADS（Δ-7），每個(gè)亞族分別有6、4、1、1、9和5個(gè)AcFAD成員（圖1）。各亞族均有獼猴桃FAD家族成員的分布暗示著AcFAD蛋白可能具有功能上的多樣性。在大多數(shù)亞族中獼猴桃與擬南芥親緣關(guān)系較近，與黃瓜的親緣關(guān)系較遠(yuǎn)。

2.3 獼猴桃FAD基因的染色體定位和共線性分析

為了研究獼猴桃FAD基因家族的遺傳差異，分析了AcFADs的染色體定位。結(jié)果（圖2）顯示獼猴桃Hong Yang v3基因組共有29條染色體，25個(gè)獼猴桃FAD基因分布在19條不同的染色體（LG）上，1個(gè)獼猴桃FAD基因（AcSLD1.1）分布在Contig（CG）01622重疊群上。大多數(shù)FAD基因定位于染色體的前中端。0和20號(hào)染色體上都分布著3個(gè)FAD基因，1和15號(hào)染色體有2個(gè)FAD基因，剩下的15條染色體上（3、4、5、6、8、9、11、13、14、16、18、21、22、23和27號(hào)）都只有1個(gè)FAD基因。由此可見，獼猴桃FAD基因和染色體存在一定的關(guān)聯(lián)性。

為了進(jìn)一步了解FAD基因家族的進(jìn)化關(guān)系，進(jìn)行了獼猴桃物種內(nèi)以及獼猴桃與擬南芥物種間的共線性分析。種內(nèi)共線性圖譜顯示有9對(duì)串聯(lián)重復(fù)基因，分別是LG0中的AcDES1.2、AcFAD7.2和AcSLD2.1；LG1的AcFAD6.2和AcADS1.1；LG3的AcFAD7.1；LG4的AcFAD5.2；LG5的AcFAD4和AcSLD1.3。此外還鑒定出22個(gè)獼猴桃FAD基因產(chǎn)生了22對(duì)片段重復(fù)基因，其中AcDES1.2（LG0）與AcDES1.1（LG20）存在共線性；AcFAD7.2（LG0）分別與AcFAD7.1（LG3）、AcFAD7.3（LG22）、AcFAD3.1（LG20）、AcFAD3.2（LG9）存在共線性；AcSLD2.1（LG0）與AcSLD2.3（LG23）、AcSLD1.1（CG01622）、AcSLD2.2（LG16）、AcSLD1.2（LG14）存在共線性；AcFAD6.2（LG1）與AcFAD6.1（LG15）、AcFAD7.1（LG3）、AcFAD3.2（LG9）存在共線性；AcADS1.1（LG1）與AcADS1.2（LG27）存在共線性；AcFAD7.1（LG3）與AcFAD3.1（LG20）、AcFAD3.2（LG9）、AcFAD7.3（LG22）、AcFAD3.3（LG15）存在共線性；AcFAD5.2（LG4）與AcFAD5.3（LG11）、AcFAD5.1（LG21）存在共線性；AcSLD1.3（LG6）與AcSLD2.1（LG0）、AcSLD2.3（LG23）、AcSLD1.1（CG01622）存在共線性（圖2）。由此可見，獼猴桃FAD基因家族的串聯(lián)重復(fù)事件伴隨有片段重復(fù)事件，所以基因重復(fù)事件（串聯(lián)重復(fù)和片段重復(fù)）對(duì)AcFADs的多樣性和進(jìn)化具有促進(jìn)作用。獼猴桃與擬南芥物種間共線性圖譜顯示，23個(gè)獼猴桃FAD基因與17個(gè)擬南芥FAD基因之間存在32對(duì)共線性關(guān)系（圖3），說明獼猴桃和擬南芥的FAD同源基因較多。

2.4 獼猴桃FAD基因結(jié)構(gòu)與蛋白保守基序分析

為了深入了解AcFADs結(jié)構(gòu)的異同，提取獼猴桃FAD基因家族成員的CDS序列及其基因組序列繪制了外顯子-內(nèi)含子基因結(jié)構(gòu)圖（圖4），結(jié)果顯示內(nèi)含子數(shù)量在0～33個(gè)之間，外顯子數(shù)量在1～34個(gè)之間（一般比對(duì)應(yīng)內(nèi)含子數(shù)目多1個(gè)），其中FAD3/FAD7家族中除了AcFAD7.3內(nèi)含子數(shù)量高達(dá)33個(gè)外，其他基因成員的內(nèi)含子數(shù)量都穩(wěn)定在7～8個(gè)；FAD2.1沒有內(nèi)含子，F(xiàn)AD2.2有1個(gè)內(nèi)含子，F(xiàn)AD6成員各有7個(gè)和5個(gè)內(nèi)含子，ADS成員均只有2個(gè)內(nèi)含子，DES成員均只有1個(gè)內(nèi)含子，F(xiàn)AD4、SLD、FAD5成員包含1～2個(gè)內(nèi)含子，F(xiàn)AB2有4個(gè)內(nèi)含子。AcSLD1.2和AcSLD2.2沒有5’UTR，AcFAB2沒有3’UTR，AcFAD7.3、AcADS1.2、AcFAD4以及SLD其余的5個(gè)成員沒有5’UTR和3’UTR。由此可見，在同一分支中，大多數(shù)成員具有相似的長(zhǎng)度和相同數(shù)量的結(jié)構(gòu)分布。

通過在線網(wǎng)站MEME預(yù)測(cè)了FAD蛋白序列中的20個(gè)保守基序（圖5），結(jié)果發(fā)現(xiàn)大多數(shù)獼猴桃FAD蛋白均含有motif 2、motif 4和motif 14，且motif 2和motif 4基本都位于C端，說明這些基序保守性較強(qiáng)，是典型的FAD結(jié)構(gòu)域，可能行使相似的功能。motif 15是FAD5蛋白特有的基序，motif 18和motif 19是DES蛋白特有的基序，F(xiàn)AD5與DES蛋白共有motif 17，F(xiàn)AD4和ADS1.2蛋白沒有保守基序，F(xiàn)AD3、FAD2和DES家族所有蛋白成員的基序種類、數(shù)量和排列方式各自完全一致。同一亞族內(nèi)的蛋白表現(xiàn)出相似的基序組成，說明同一亞族成員高度保守，具有較近的親緣關(guān)系。但亞族間也存在明顯差異，說明一些基因在進(jìn)化過程中發(fā)生了功能分化。

2.5 獼猴桃FAD基因啟動(dòng)子順式作用元件分析

為探索獼猴桃FAD基因家族順式作用元件的功能，對(duì)該家族上游1500 bp的啟動(dòng)子序列進(jìn)行了順式作用元件預(yù)測(cè)（圖6）。結(jié)果預(yù)測(cè)出47種順式作用元件，其中光響應(yīng)元件數(shù)量最多，共計(jì)298個(gè)（ACE 2個(gè)，AE-box 8個(gè)，AT1-motif 4個(gè)，ATC-motif 5個(gè)，ATCT-motif 7個(gè)，Box 4 68個(gè)，Box Ⅱ 1個(gè)，CAG-motif 1個(gè)，chs-CMA1a 5個(gè)，chs-CMA2a 4個(gè)，GA-motif 4個(gè)，GATA-motif 20個(gè)，G-Box 62個(gè)，GT1-motif 28個(gè)，I-box 9個(gè)，LAMP-element 4個(gè)，MRE 22個(gè)，Sp1 12個(gè)，TCCC-motif 10個(gè)，TCT-motif 21個(gè)，TGA-element 1個(gè)）；其次是植物激素響應(yīng)元件（196個(gè)），包括茉莉酸甲酯響應(yīng)元件72個(gè)（CGTCA-motif 36個(gè)，TGACG-motif 36個(gè)）、脫落酸響應(yīng)元件59個(gè)（ABRE）、生長(zhǎng)素響應(yīng)元件24個(gè)（AAAC-motif 1個(gè)，AuxRE 1個(gè)，AuxRR-core 4個(gè)，TGA-box 2個(gè)，TGA-element 16個(gè)）、水楊酸響應(yīng)元件21個(gè)（TCA-element）和赤霉素響應(yīng)元件20個(gè)（GARE-motif 8個(gè)，P-box 7個(gè)，TATC-box 5個(gè)）；再者是逆境脅迫響應(yīng)元件（110個(gè)），包括厭氧誘導(dǎo)所必需的元件69個(gè)（ARE 60個(gè)，GC-motif 9個(gè)）、低溫響應(yīng)元件23個(gè)（LTR）、參與防御和應(yīng)激反應(yīng)元件11個(gè)（TC-rich repeats）和參與干旱誘導(dǎo)的MYB結(jié)合位點(diǎn)元件7個(gè)（MBS）；最后是植物生長(zhǎng)發(fā)育相關(guān)元件（54個(gè)），包括MYBHv1結(jié)合位點(diǎn)元件13個(gè)（CCAAT-box）、與分生組織表達(dá)有關(guān)的元件12個(gè)（CAT-box）、參與玉米蛋白代謝調(diào)控的元件8個(gè)（O2-site）、參與胚乳表達(dá)的元件7個(gè)（AACA_motif 1個(gè)，GCN4_motif 6個(gè)）、參富含AT的DNA結(jié)合蛋白（ATBP-1）的結(jié)合位點(diǎn)元件6個(gè)（AT-rich element）、與晝夜節(jié)律控制有關(guān)的元件3個(gè)（circadian）、參與細(xì)胞周期調(diào)控的元件2個(gè)（MSA-like）、參與葉綠體中胚層細(xì)胞分化的元件2個(gè)（HD-Zip 1）和參與黃酮類生物合成基因調(diào)控的MYB結(jié)合位點(diǎn)元件1個(gè)（MBSI）。此外，還有3個(gè)最大誘導(dǎo)劑介導(dǎo)的激活元件（2個(gè)）（AT-rich sequence）。除了AcADS1.1外，其他AcFADs啟動(dòng)子序列上均有激素響應(yīng)元件，因此這些基因極有可能參與激素調(diào)控的生理過程；所有獼猴桃FAD基因家族成員都含有逆境脅迫響應(yīng)元件，所以獼猴桃FAD基因在逆境脅迫應(yīng)答過程中發(fā)揮著重要的作用；AcFAD6.2、AcFAD3.2、AcADS1.1和AcADS1.2沒有調(diào)控生長(zhǎng)發(fā)育的順式作用元件。

2.6 獼猴桃FAD基因家族表達(dá)模式分析

為了了解獼猴桃FAD基因在采后成熟過程中的表達(dá)模式，對(duì)紅陽獼猴桃采后0、4、10 d果肉樣品轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行了熱圖繪制以表征其表達(dá)豐度。結(jié)果（圖7）顯示，AcFAD3.1、AcFAD3.3、AcSLD1.4和AcDES1.1這4個(gè)基因成員在獼猴桃整個(gè)采后成熟過程中均不表達(dá)；AcFAD6.1、AcFAD7.1、AcSLD1.2和AcSLD2.1在采后成熟早期（常溫貯藏0 d）表現(xiàn)出最高的表達(dá)水平，隨著后熟時(shí)間的增加其表達(dá)水平也隨之下降，在成熟后期（常溫貯10 d）達(dá)到最低的表達(dá)水平；AcFAD3.2、AcFAD5.3、AcADS1.1、AcADS1.2和AcFAD6.2同樣在采后成熟早期表現(xiàn)出最高的表達(dá)水平，但在采后成熟中期（常溫貯4 d）表現(xiàn)出最低的表達(dá)水平，后期其表達(dá)水平又有所提高；AcFAD7.3、AcFAD5.1在成熟早期和后期的表達(dá)水平一致，中期略有下降；AcFAD4僅在成熟早期和中期表達(dá)且表達(dá)水平較低；AcSLD2.3在成熟中期和后期表達(dá)水平基本一致，早期表達(dá)水平最低；AcSLD1.1在成熟后期表達(dá)水平最高，中期表達(dá)水平最低；AcFAD2.1只在成熟后期表達(dá)且表達(dá)水平極低；AcDES1.2、AcSLD1.3、AcFAD5.2、AcFAD7.2、AcSLD2.2、AcFAB2和AcFAD2.2的表達(dá)水平均隨著成熟而不斷提高，且值得注意的是AcFAD2.2的表達(dá)水平極高且各樣點(diǎn)間表達(dá)水平增加1.4倍左右，后期相比于早期增加了2倍。因此，獼猴桃FAD基因家族各成員在獼猴桃采后成熟過程中有不同的表達(dá)模式，暗示著其在成熟過程中的各階段發(fā)揮著不同的作用。

2.7 獼猴桃果實(shí)成熟過程中硬度和不飽和脂肪酸含量變化

在獼猴桃采后成熟過程中硬度呈下降趨勢(shì)，并于采后10 d達(dá)到了可食用成熟度（圖8-A）。不飽和脂肪酸是獼猴桃香氣合成的重要前體物質(zhì)，對(duì)采后獼猴桃果實(shí)中不飽和脂肪酸含量進(jìn)行測(cè)定，發(fā)現(xiàn)油酸（OA，18∶1）含量隨著獼猴桃采后成熟過程不斷下降，亞油酸（LA，18∶2）則表現(xiàn)出與之相反的趨勢(shì)，亞麻酸（LeA，18∶3）含量在成熟早期和中期無顯著差別，而在后期含量顯著降低（圖8-B～D）。

2.8 獼猴桃FAD基因家族熒光定量分析

從獼猴桃轉(zhuǎn)錄組熱圖中挑選出在成熟后期差異表達(dá)上調(diào)的8個(gè)基因（AcSLD1.1、AcSLD1.3、AcSLD2.2、AcSLD2.3、AcFAD5.2、AcFAD7.2、AcFAB2和AcFAD2.2）進(jìn)行了qPCR驗(yàn)證。結(jié)果（圖9）顯示，在獼猴桃整個(gè)采后成熟過程中AcSLD1.3表達(dá)量無顯著變化；AcSLD2.2表達(dá)量先上調(diào)后下調(diào)但變化幅度不大；AcSLD1.1、AcSLD2.3和AcFAB2表達(dá)量在成熟早期和中期無顯著變化，在后期顯著上調(diào)，但表達(dá)量較低；與早期相比，AcFAD5.2在后期表達(dá)量顯著上調(diào)，中期表達(dá)量相對(duì)于早期和后期無顯著變化；AcFAD7.2和AcFAD2.2表達(dá)量都是隨著成熟而不斷顯著上調(diào)，但AcFAD7.2表達(dá)量并不是很高，中期表達(dá)量相比于早期上調(diào)2.42倍，后期上調(diào)量是早期的13.09倍，而AcFAD2.2中期表達(dá)量比早期上調(diào)8.44倍，后期表達(dá)量比中期上調(diào)9.05倍，后期上調(diào)量是早期的76.43倍，這與轉(zhuǎn)錄組數(shù)據(jù)的變化模式一致，因此推測(cè)AcFAD2.2極有可能是參與編碼獼猴桃采后成熟過程中w-6脂肪酸去飽和酶的關(guān)鍵候選基因，使得單不飽和脂肪酸轉(zhuǎn)化為雙不飽和脂肪酸，從而導(dǎo)致了油酸的減少和亞油酸的積累。

3 討 論

脂肪酸去飽和酶對(duì)植物生長(zhǎng)發(fā)育、抵抗生物與非生物脅迫和品質(zhì)形成等方面具有正向的積極作用[1-4]。近年來，已經(jīng)在多種植物中對(duì)FAD基因家族進(jìn)行了鑒定，如亞麻中有43個(gè)[37]、茄子有38個(gè)[19]、核桃中有24個(gè)[38]、桃中有6個(gè)FAD基因家族成員[39]，而本研究中筆者首次在獼猴桃中鑒定出了26個(gè)FAD基因。由此可見，不同物種間FAD成員數(shù)量存在一定差異，這可能是其基因組大小不同和進(jìn)化過程中某些基因的復(fù)制、分化或缺失造成的。對(duì)獼猴桃FAD家族蛋白序列進(jìn)行理化性質(zhì)分析發(fā)現(xiàn)，其分子質(zhì)量與氨基酸數(shù)量成正比，而在茄子[19]、番茄[20]和亞麻[37]等中也有同樣的規(guī)律。獼猴桃FAD的亞細(xì)胞定位結(jié)果發(fā)現(xiàn)FAD2s和多數(shù)SLDs/DESs都定位在內(nèi)質(zhì)網(wǎng)中，F(xiàn)AD6s定位在葉綠體中，F(xiàn)AD3s/FAD7s定位在葉綠體和內(nèi)質(zhì)網(wǎng)中，說明不同成員在這些不同的特定部位中發(fā)揮作用，與前人的研究結(jié)果一致[9-10，12]。獼猴桃與19個(gè)擬南芥FAD同樣都可以分為6大亞族，在FAB2和FAD4亞族中這兩個(gè)物種都只有一個(gè)FAD成員，而在其他亞族中成員數(shù)量分布存在較大的差異（尤其是SLD和FAD3/FAD7/FAD8亞族），該現(xiàn)象與番茄FAD家族相似[20]，暗示著這些亞族基因在進(jìn)化過程中可能發(fā)生了特異性擴(kuò)增。本研究在獼猴桃FAD家族中鑒定出9對(duì)串聯(lián)重復(fù)基因和22對(duì)片段重復(fù)基因，同樣地在甘藍(lán)型油菜FAD家族中存在3對(duì)串聯(lián)重復(fù)和25對(duì)片段重復(fù)基因[8]、在小麥中形成26對(duì)串聯(lián)重復(fù)和126對(duì)片段重復(fù)基因[40]，因此推測(cè)FAD基因家族的擴(kuò)增主要是基因片段引起的。具有相似的外顯子-內(nèi)含子基因結(jié)構(gòu)的AcFAD成員通常聚為同一分支，且其編碼的蛋白具有相似的基序組成，但也有部分成員（AcFAD4和AcADS1.2）丟失了保守基序，在茄子FAD相應(yīng)亞族中也發(fā)現(xiàn)此現(xiàn)象[19]，說明FAD家族在物種間具有相對(duì)保守的進(jìn)化趨勢(shì)。對(duì)獼猴桃FAD啟動(dòng)子序列進(jìn)行分析，發(fā)現(xiàn)大量的光響應(yīng)元件、植物激素響應(yīng)元件（茉莉酸甲酯，脫落酸，生長(zhǎng)素，水楊酸，赤霉素）、逆境脅迫響應(yīng)元件（厭氧，低溫，干旱）和生長(zhǎng)發(fā)育相關(guān)元件，尤其是逆境脅迫響應(yīng)元件遍布在每個(gè)成員中，說明獼猴桃FAD的表達(dá)受到光的調(diào)控，并廣泛參與了植物生長(zhǎng)發(fā)育、成熟及抗逆過程，在其他呼吸躍變型果實(shí)中（如番茄和香蕉等）也發(fā)現(xiàn)了相同的順式作用元件[18，20]。

果實(shí)風(fēng)味作為水果重要的品質(zhì)性狀，是在成熟過程中形成的品質(zhì)之一，也是影響其銷量以及消費(fèi)者偏好的重要因素，果實(shí)的風(fēng)味取決于糖酸以及揮發(fā)性芳香化合物的種類和含量。獼猴桃作為典型的呼吸躍變型果實(shí)，在成熟過程中，存在明顯的呼吸躍變并促發(fā)內(nèi)源乙烯大量合成，從而加速各種揮發(fā)性化合物的合成，逐漸產(chǎn)生獨(dú)特的果香。在LOX途徑中，亞油酸和亞麻酸首先被氧化為脂肪酸氫過氧化物，隨后被氫過氧化物裂解酶裂解形成己醛和己烯醛。然后，C6醛通過醇脫氫酶還原為相應(yīng)的C6醇，隨后通過醇?；D(zhuǎn)移酶轉(zhuǎn)化為酯[30]。哈密瓜亞麻酸含量在貯藏過程中均有所增加，貯藏前期油酸含量急劇增多而中后期則呈下降趨勢(shì)，亞油酸含量與之相反，呈上升趨勢(shì)[41]，這與本研究中發(fā)現(xiàn)獼猴桃在采后成熟過程中油酸含量隨成熟進(jìn)程而顯著下降，亞油酸表現(xiàn)出與之負(fù)相關(guān)的趨勢(shì)，亞麻酸含量在早中期無顯著變化，而在后期有所下降的現(xiàn)象一致。同時(shí)，也有研究表明油酸在油棕果實(shí)采后成熟階段含量均減少，而亞油酸只在成熟早期階段減少，卻在后期階段顯著增加[42]。番茄果實(shí)成熟階段亞油酸含量最高，油酸含量下降，亞油酸含量逐漸升高，亞麻酸含量無顯著變化[43]。以上研究表明大多數(shù)果實(shí)成熟過程往往伴有亞油酸含量的增加和油酸含量的下降。但也有研究發(fā)現(xiàn)了與本研究不一致的現(xiàn)象，如杧果采后成熟過程中亞油酸含量逐漸下降，亞麻酸含量逐漸升高[44]，在香蕉果實(shí)的成熟進(jìn)程中亞麻酸增加，呼吸躍變前期亞麻酸含量增加緩慢，至呼吸躍變上升期迅速增加，而亞油酸含量則有所下降[45]。這可能是熱帶水果和亞熱帶水果差異及其芳香物質(zhì)種類不同所致。

基于FAD家族在獼猴桃成熟各階段的表達(dá)模式，發(fā)現(xiàn)了AcFAD7.2和AcFAD2.2這兩個(gè)表達(dá)量均隨著成熟而不斷顯著上調(diào)的家族成員。前人研究顯示，在油茶和擬南芥中分別同源和異源過表達(dá)ω-3脂肪酸去飽和酶CoFAD7后，亞麻酸的含量顯著增加并且提高了亞麻酸/亞油酸的比值[46]。此外，在煙草中分別過量表達(dá)桃 PpFAD3-1和PpFAD3-2后發(fā)現(xiàn)亞油酸含量大幅下降，亞麻酸含量顯著上升[47]。但在本研究中并未發(fā)現(xiàn)亞麻酸含量增多，反而出現(xiàn)下降，說明獼猴桃成熟過程中亞麻酸的合成速率可能低于其作為前體物質(zhì)參與芳香物質(zhì)合成的速率，造成亞麻酸積累量降低，因此編碼ω-3脂肪酸去飽和酶基因AcFAD7.2在獼猴桃采后成熟不飽和脂肪酸合成和積累中發(fā)揮的作用并不是很顯著。對(duì)大腸桿菌中異源表達(dá)桃的PpFAD2-1和PpFAD2-2及其酶活性分析，發(fā)現(xiàn)這兩個(gè)ω-6脂肪酸去飽和酶基因都能促進(jìn)油酸向亞油酸的去飽和[48]；對(duì)橄欖進(jìn)行OeFAD2-2和OeFAD2-5超量表達(dá)，發(fā)現(xiàn)其具有將油酸轉(zhuǎn)化為亞油酸并提高亞油酸含量的能力[49]；在番茄中進(jìn)行異源酵母表達(dá)SlFAD2-1和SlFAD2-2能將油酸轉(zhuǎn)化生成亞油酸并顯著增加了亞油酸含量[50]。在本研究中也同樣發(fā)現(xiàn)了油酸含量減少和亞油酸含量增多的現(xiàn)象，作為與番茄SlFAD2-1和SlFAD2-2高度同源的ω-6脂肪酸去飽和酶基因，AcFAD2.2的表達(dá)水平也隨著成熟而不斷攀高，其表達(dá)模式與其發(fā)揮作用的兩種不飽和脂肪酸（油酸和亞油酸）含量變化存在明顯的正相關(guān)，說明AcFAD2.2是參與獼猴桃采后成熟過程中脂肪酸代謝的關(guān)鍵基因，在不飽和脂肪酸合成與積累過程中發(fā)揮著極為重要的作用。

4 結(jié) 論

筆者在獼猴桃中鑒定出了26個(gè)FAD基因家族成員，根據(jù)其發(fā)揮作用的部位被定位到細(xì)胞的不同部位，因此絕大多數(shù)成員定位于內(nèi)質(zhì)網(wǎng)或葉綠體中，這些成員隨機(jī)分布在19條不同染色體上，片段重復(fù)事件是獼猴桃FAD擴(kuò)增和進(jìn)化的主要?jiǎng)恿Γ讷J猴桃FAD啟動(dòng)子序列上發(fā)現(xiàn)大量的光響應(yīng)元件、植物激素響應(yīng)元件（茉莉酸甲酯、脫落酸、生長(zhǎng)素、水楊酸、赤霉素）、逆境脅迫響應(yīng)元件（厭氧、低溫、干旱）和生長(zhǎng)發(fā)育相關(guān)元件。獼猴桃采后成熟過程中油酸含量下降、亞油酸含量增多，并且AcFAD2.2表達(dá)水平隨著成熟而不斷顯著上調(diào)，由此可以推斷出AcFAD2.2是油酸去飽和生成亞油酸的關(guān)鍵酶基因，并且可能會(huì)通過光、激素等響應(yīng)方式在獼猴桃采后成熟過程中表達(dá)來調(diào)控不飽和脂肪酸的合成與積累。因此，今后可借助分子生物學(xué)技術(shù)對(duì)此關(guān)鍵候選基因進(jìn)行克隆、功能驗(yàn)證和上游調(diào)控試驗(yàn)，深入解析獼猴桃FAD基因在采后成熟過程中不飽和脂肪酸積累和香氣物質(zhì)合成中的功能機(jī)制和調(diào)控網(wǎng)絡(luò)。

參考文獻(xiàn) References：

[1] HE M，QIN C X，WANG X，DING N Z. Plant unsaturated fatty acids：Biosynthesis and regulation[J]. Frontiers in Plant Science，2020，11：390.

[2] FAN R S，LI L，CAI G，YE J，LIU M H，WANG S H，LI Z Q. Molecular cloning and function analysis of FAD2 gene in Idesia polycarpa[J]. Phytochemistry，2019，168：112114.

[3] 王利民，符真珠，高杰，董曉宇，張晶，袁欣，蔣卉，王慧娟，李艷敏，師曼，張和臣. 植物不飽和脂肪酸的生物合成及調(diào)控[J]. 基因組學(xué)與應(yīng)用生物學(xué)，2020，39（1）：254-258.

WANG Limin，F(xiàn)U Zhenzhu，GAO Jie，DONG Xiaoyu，ZHANG Jing，YUAN Xin，JIANG Hui，WANG Huijuan，LI Yanmin，SHI Man，ZHANG Hechen. Molecular mechanism of unsaturated fatty acids synthesis and regulation in plant[J]. Genomics and Applied Biology，2020，39（1）：254-258.

[4] 薛曉夢(mèng)，李建國(guó)，白冬梅，晏立英，萬麗云，康彥平，淮東欣，雷永，廖伯壽. 花生FAD2基因家族表達(dá)分析及其對(duì)低溫脅迫的響應(yīng)[J]. 作物學(xué)報(bào)，2019，45（10）：1586-1594.

XUE Xiaomeng，LI Jianguo，BAI Dongmei，YAN Liying，WAN Liyun，KANG Yanping，HUAI Dongxin，LEI Yong，LIAO Boshou. Expression profiles of FAD2 genes and their responses to cold stress in peanut[J]. Acta Agronomica Sinica，2019，45（10）：1586-1594.

[5] KUGLER A，ZORIN B，DIDI-COHEN S，SIBIRYAK M，GORELOVA O，ISMAGULOVA T，KOKABI K，KUMARI P，LUKYANOV A，BOUSSIBA S，SOLOVCHENKO A，KHOZIN-GOLDBERG I. Long-chain polyunsaturated fatty acids in the green microalga Lobosphaera incisa contribute to tolerance to abiotic stresses[J]. Plant and Cell Physiology，2019，60（6）：1205-1223.

[6] 蔡洪芳. 1-MCP/NO/MeJA對(duì)采后桃果實(shí)脂肪酸途徑香氣物質(zhì)的調(diào)控研究[D]. 南京：南京農(nóng)業(yè)大學(xué)，2020.

CAI Hongfang. Study on the regulation of 1-MCP/NO/MeJA on aroma volatiles in fatty acid pathway of postharvest peach fruit[D]. Nanjing：Nanjing Agricultural University，2020.

[7] 趙訓(xùn)超，魏玉磊，丁冬，劉夢(mèng)，蓋勝男，張今杰，邵文靜，李嘉欣，徐晶宇. 甜蕎麥脂肪酸脫氫酶基因（FeFAD）家族的鑒定與分析[J]. 東北農(nóng)業(yè)科學(xué)，2021，46（1）：36-41.

ZHAO Xunchao，WEI Yulei，DING Dong，LIU Meng，GAI Shengnan，ZHANG Jinjie，SHAO Wenjing，LI Jiaxin，XU Jingyu. Genome-wide identification and bioinformatics analysis of fatty acid desaturase gene （FeFAD） family in common buckwheat[J]. Journal of Northeast Agricultural Sciences，2021，46（1）：36-41.

[8] XUE Y F，CHEN B J，WANG R，WIN A N，LI J N，CHAI Y R. Genome-wide survey and characterization of fatty acid desaturase gene family in Brassica napus and its parental species[J]. Applied Biochemistry and Biotechnology，2018，184（2）：582-598.

[9] DU C，CHEN Y Y，WANG K，YANG Z，ZHAO C Z，JIA Q L，TAYLOR D C，ZHANG M. Strong co-suppression impedes an increase in polyunsaturated fatty acids in seeds overexpressing FAD2[J]. Journal of Experimental Botany，2019，70（3）：985-994.

[10] PENG Z Y，RUAN J，TIAN H Y，SHAN L，MENG J J，GUO F，ZHANG Z M，DING H，WAN S B，LI X G. The family of peanut fatty acid desaturase genes and a functional analysis of four ω-3 AhFAD3 members[J]. Plant Molecular Biology Reporter，2020，38（2）：209-221.

[11] SHANKLIN J，CAHOON E B. Desaturation and related modifications of fatty acids[J]. Annual Review of Plant Physiology and Plant Molecular Biology，1998，49：611-641.

[12] 李昊遠(yuǎn)，郝翠翠，陳明娜，陳娜，王冕，潘麗娟，王通，禹山林，侯艷華，遲曉元. 花生鞘脂Δ8去飽和酶基因（AhSLD2）的克隆與表達(dá)分析[J]. 花生學(xué)報(bào)，2018，47（2）：24-29.

LI Haoyuan，HAO Cuicui，CHEN Mingna，CHEN Na，WANG Mian，PAN Lijuan，WANG Tong，YU Shanlin，HOU Yanhua，CHI Xiaoyuan. Cloning and expression analysis of sphingolipid Δ8 desaturase （AhSLD2） gene in peanut[J]. Journal of Peanut Science，2018，47（2）：24-29.

[13] MICHAELSON L V，Z?UNER S，MARKHAM J E，HASLAM R P，DESIKAN R，MUGFORD S，ALBRECHT S，WARNECKE D，SPERLING P，HEINZ E，NAPIER J A. Functional characterization of a higher plant sphingolipid Delta4-desaturase：Defining the role of sphingosine and sphingosine-1-phosphate in Arabidopsis[J]. Plant Physiology，2009，149（1）：487-498.

[14] ARONDEL V，LEMIEUX B，HWANG I，GIBSON S，GOODMAN H M，SOMERVILLE C R. Map-based cloning of a gene controlling omega-3 fatty acid desaturation in Arabidopsis[J]. Science，1992，258（5086）：1353-1355.

[15] OKULEY J，LIGHTNER J，F(xiàn)ELDMANN K，YADAV N，LARK E，BROWSE J. Arabidopsis FAD2 gene encodes the enzyme that is essential for polyunsaturated lipid synthesis[J]. The Plant Cell，1994，6（1）：147-158.

[16] IBA K，GIBSON S，NISHIUCHI T，F(xiàn)USE T，NISHIMURA M，ARONDEL V，HUGLY S，SOMERVILLE C. A gene encoding a chloroplast ω-3 fatty acid desaturase complements alterations in fatty acid desaturation and chloroplast copy number of the fad7 mutant of Arabidopsis thaliana[J]. Journal of Biological Chemistry，1993，268（32）：24099-24105.

[17] REDDY A S，NUCCIO M L，GROSS L M，THOMAS T L. Isolation of a Δ6-desaturase gene from the Cyanobacterium synechocystis sp. strain PCC 6803 by gain-of-function expression in Anabaena sp. strain PCC 7120[J]. Plant Molecular Biology，1993，22（2）：293-300.

[18] CHENG C Z，LIU F，SUN X L，WANG B，LIU J P，NI X T，HU C H，DENG G M，TONG Z，ZHANG Y Y，Lü P T. Genome-wide identification of FAD gene family and their contributions to the temperature stresses and mutualistic and parasitic fungi colonization responses in banana[J]. International Journal of Biological Macromolecules，2022，204：661-676.

[19] 朱宗文，張愛冬，吳雪霞，查丁石. 生物信息學(xué)鑒定分析茄子脂肪酸去飽和酶（FAD）基因家族[J]. 分子植物育種，2023，21（8）：2453-2463.

ZHU Zongwen，ZHANG Aidong，WU Xuexia，ZHA Dingshi. Identification and bioinformatics analysis of fatty acid desaturase （FAD） gene family in eggplant （Solanum melongena L.）[J]. Molecular Plant Breeding，2023，21（8）：2453-2463.

[20] 張明亞，龐勝群，劉玉東，蘇永峰，牛博文，韓瓊瓊. 番茄FAD基因家族的鑒定與表達(dá)分析[J]. 生物技術(shù)通報(bào)，2024，40（7）：150-162.

ZHANG Mingya，PANG Shengqun，LIU Yudong，SU Yongfeng，NIU Bowen，HAN Qiongqiong. Identification and expression analysis of FAD gene family in Solanum lycopersicum[J]. Biotechnology Bulletin，2024，40（7）：150-162.

[21] 高謝旺，譚安琪，胡信暢，祝孟洋，阮穎，劉春林. 利用CRISPR/Cas9技術(shù)創(chuàng)制高油酸甘藍(lán)型油菜新種質(zhì)[J]. 植物遺傳資源學(xué)報(bào)，2020，21（4）：1002-1008.

GAO Xiewang，TAN Anqi，HU Xinchang，ZHU Mengyang，RUAN Ying，LIU Chunlin. Creation of new germplasm of high-oleic rapeseed using CRISPR/Cas9[J]. Journal of Plant Genetic Resources，2020，21（4）：1002-1008.

[22] 王慶華，王磊，吳文江，郭家選，沈元月，吳國(guó)良. 果實(shí)香氣物質(zhì)的合成及其激素調(diào)控研究進(jìn)展[J/OL]. 分子植物育種，2021：1-11（2021-11-02）[2024-09-07]. https：//kns.cnki.net/kcms/detail/46.1068.S.20211029.1846.006.html.

WANG Qinghua，WANG Lei，WU Wenjiang，GUO Jiaxuan，SHEN Yuanyue，WU Guoliang. Advances in aroma compounds biosynthesis and hormone regulation of fruit[J/OL]. Molecular Plant Breeding，2021：1-11（2021-11-02）[2024-09-07]. https：//kns.cnki.net/kcms/detail/46.1068.S.20211029.1846.006.html.

[23] 蔡璨，白玉，韓藝，郭佳欣，沙偉，張梅娟，彭疑芳，馬天意. 植物多不飽和脂肪酸的研究進(jìn)展[J]. 高師理科學(xué)刊，2023，43（9）：64-69.

CAI Can，BAI Yu，HAN Yi，GUOJIA Xin，SHA Wei，ZHANG Meijuan，PENG Yifang，MA Tianyi. Research progress of plant polyunsaturated fatty acids[J]. Journal of Science of Teachers’ College and University，2023，43（9）：64-69.

[24] 陳成，王依，楊勇，閻永齊. 采收成熟度對(duì)‘金艷’獼猴桃果實(shí)品質(zhì)及香氣成分的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2020，36（31）：28-36.

CHEN Cheng，WANG Yi，YANG Yong，YAN Yongqi. Effects of maturity stage on fruit quality and aroma components of ‘Jinyan’ kiwifruit[J]. Chinese Agricultural Science Bulletin，2020，36（31）：28-36.

[25] 陳義挺，賴瑞聯(lián)，馮新，程春振，鐘春水，高敏霞，吳如健. 不同貯藏條件下獼猴桃香氣成分的變化規(guī)律研究[J]. 熱帶作物學(xué)報(bào)，2020，41（6）：1251-1256.

CHEN Yiting，LAI Ruilian，F(xiàn)ENG Xin，CHENG Chunzhen，ZHONG Chunshui，GAO Minxia，WU Rujian. Change of aroma components in different storage conditions of kiwifruit[J]. Chinese Journal of Tropical Crops，2020，41（6）：1251-1256.

[26] COZZOLINO R，DE GIULIO B，PETRICCIONE M，MARTIGNETTI A，MALORNI L，ZAMPELLA L，LAURINO C，PELLICANO M P. Comparative analysis of volatile metabolites，quality and sensory attributes of Actinidia chinensis fruit[J]. Food Chemistry，2020，316：126340.

[27] 張波，李鮮，陳昆松. 基于EST庫的獼猴桃脂氧合酶基因家族成員的克隆[J]. 園藝學(xué)報(bào)，2008，35（3）：337-342.

ZHANG Bo，LI Xian，CHEN Kunsong. Molecular cloning of lipoxygenase gene family members in kiwifruit based on EST database[J]. Acta Horticulturae Sinica，2008，35（3）：337-342.

[28] CROWHURST R N，GLEAVE A P，MACRAE E A，AMPOMAH-DWAMENA C，ATKINSON R G，BEUNING L L，BULLEY S M，CHAGNE D，MARSH K B，MATICH A J，MONTEFIORI M，NEWCOMB R D，SCHAFFER R J，USADEL B，ALLAN A C，BOLDINGH H L，BOWEN J H，DAVY M W，ECKLOFF R，F(xiàn)ERGUSON A R，F(xiàn)RASER L G，GERA E，HELLENS R P，JANSSEN B J，KLAGES K，LO K R，MACDIARMID R M，NAIN B，MCNEILAGE M A，RASSAM M，RICHARDSON A C，RIKKERINK E H，ROSS G S，SCHR?DER R，SNOWDEN K C，SOULEYRE E J F，TEMPLETON M D，WALTON E F，WANG D，WANG M Y，WANG Y Y，WOOD M，WU R M，YAUK Y K，LAING W A. Analysis of expressed sequence tags from Actinidia：Applications of a cross species EST database for gene discovery in the areas of flavor，health，color and ripening[J]. BMC Genomics，2008，9：351.

[29] ZHANG J Y，HUANG S N，CHEN Y H，WANG G，GUO Z R. Identification and characterization of two waterlogging responsive alcohol dehydrogenase genes （AdADH1 and AdADH2） in Actinidia deliciosa[J]. Molecular Breeding，2017，37（4）：52.

[30] 陶淑華，陳麗，蔣鎮(zhèn)燁，宋倩倩，宋亦超，姜天甲，鄭小林. 低溫貯藏對(duì)美味獼猴桃布魯諾果實(shí)主要揮發(fā)性物質(zhì)和脂肪酸代謝的影響[J]. 核農(nóng)學(xué)報(bào)，2020，34（2）：288-297.

TAO Shuhua，CHEN Li，JIANG Zhenye，SONG Qianqian，SONG Yichao，JIANG Tianjia，ZHENG Xiaolin. Effects of lower temperature on flavor components and fatty acid pathway in harvested kiwifruit （Actinidia deliciosa cv. Bruno）[J]. Journal of Nuclear Agricultural Sciences，2020，34（2）：288-297.

[31] 李可，袁懷瑜，朱永清，周艷，夏陳，趙楠，李華佳. 不同品種獼猴桃籽油脂肪酸組成的PCA分析[J]. 中國(guó)調(diào)味品，2021，46（2）：70-74.

LI Ke，YUAN Huaiyu，ZHU Yongqing，ZHOU Yan，XIA Chen，ZHAO Nan，LI Huajia. PCA analysis of fatty acid composition of kiwi seed oils with different varieties[J]. China Condiment，2021，46（2）：70-74.

[32] ZHANG B，YIN X R，LI X，YANG S L，F(xiàn)ERGUSON I B，CHEN K S. Lipoxygenase gene expression in ripening kiwifruit in relation to ethylene and aroma production[J]. Journal of Agricultural and Food Chemistry，2009，57（7）：2875-2881.

[33] 甘武. 獼猴桃果實(shí)品質(zhì)和香氣成分分析研究[D]. 南昌：江西農(nóng)業(yè)大學(xué)，2018.

GAN Wu. Analysis of fruit quality and aroma components of kiwifruit[D]. Nanchang：Jiangxi Agricultural University，2018.

[34] DONG C J，CAO N，ZHANG Z G，SHANG Q M. Characterization of the fatty acid desaturase genes in cucumber：Structure，phylogeny，and expression patterns[J]. PLoS One，2016，11（3）：e0149917.

[35] GAN Z Y，SHAN N，F(xiàn)EI L Y，WAN C P，CHEN J Y. Isolation of the 9-Cis-epoxycarotenoid dioxygenase （NCED） gene from kiwifruit and its effects on postharvest softening and ripening[J]. Scientia Horticulturae，2020，261：109020.

[36] 袁馨，徐云鶴，張雨培，單楠，陳楚英，萬春鵬，開文斌，翟夏琬，陳金印，甘增宇. 獼猴桃后熟過程中ABA響應(yīng)結(jié)合因子AcAREB1調(diào)控AcGH3.1的表達(dá)[J]. 園藝學(xué)報(bào)，2023，50（1）：53-64.

YUAN Xin，XU Yunhe，ZHANG Yupei，SHAN Nan，CHEN Chuying，WAN Chunpeng，KAI Wenbin，ZHAI Xiawan，CHEN Jinyin，GAN Zengyu. Studies on AcAREB1 regulating the expression of AcGH3.1 during postharvest ripening of kiwifruit[J]. Acta Horticulturae Sinica，2023，50（1）：53-64.

[37] 侯靜靜，趙利，王斌. 亞麻FAD基因家族的生物信息學(xué)鑒定分析[J]. 寒旱農(nóng)業(yè)科學(xué)，2023（3）：246-253.

HOU Jingjing，ZHAO Li，WANG Bin. Identification and bioinformatics analysis of FAD gene family in Linum usitatissimum L.[J]. Journal of Cold-Arid Agricultural Sciences，2023（3）：246-253.

[38] LIU K，ZHAO S G，WANG S，WANG H X，ZHANG Z H. Identification and analysis of the FAD gene family in walnuts （Juglans regia L.） based on transcriptome data[J]. BMC Genomics，2020，21（1）：299.

[39] 金正楠. 轉(zhuǎn)錄因子PpNAC1和表觀遺傳修飾通過調(diào)控PpFAD3-1表達(dá)參與桃果實(shí)芳香物質(zhì)合成[D]. 杭州：浙江大學(xué)，2022.

JIN Zhengnan. Transcription factor PpNAC1 and epigenetic modification are involved in the synthesis of volatiles in peach fruits by regulating the expression of PpFAD3-1[D]. Hangzhou：Zhejiang University，2022.

[40] HAJIAHMADI Z，ABEDI A，WEI H，SUN W B，RUAN H H，ZHUGE Q，MOVAHEDI A. Identification，evolution，expression，and docking studies of fatty acid desaturase genes in wheat （Triticum aestivum L.）[J]. BMC Genomics，2020，21（1）：778.

[41] 王靜，茅林春，楊璐，李學(xué)文，張輝，呂卓，劉彩虹，李乾，侯琛元. 草酸處理對(duì)采后哈密瓜果實(shí)膜脂代謝的影響[J]. 中國(guó)食品學(xué)報(bào)，2019，19（8）：189-198.

WANG Jing，MAO Linchun，YANG Lu，LI Xuewen，ZHANG Hui，Lü Zhuo，LIU Caihong，LI Qian，HOU Chenyuan. Effect of oxalic acid on reduction of membrane lipids metabolism of Hami melon’s fruit in postharvest[J]. Journal of Chinese Institute of Food Science and Technology，2019，19（8）：189-198.

[42] 吳秋妃，楊程，張淑巖，韋露，馮美利，李睿，周麗霞，曹紅星. 油棕果實(shí)發(fā)育和采后脂肪酸合成轉(zhuǎn)錄代謝差異分析[J]. 熱帶作物學(xué)報(bào)，2024，45（2）：234-246.

WU Qiufei，YANG Cheng，ZHANG Shuyan，WEI Lu，F(xiàn)ENG Meili，LI Rui，ZHOU Lixia，CAO Hongxing. Differential analysis of fatty acid synthesis，transcriptional metabolism during fruit development and postharvest in oil palm[J]. Chinese Journal of Tropical Crops，2024，45（2）：234-246.

[43] STAVECKIEN? J，KULAITIEN? J，LEVICKIEN? D，VAITKEVI?IEN? N. Changes in fatty acid content in Solanum spp. fruits during ripening[J]. Plants，2023，12（2）：268.

[44] DESHPANDE A B，CHIDLEY H G，OAK P S，PUJARI K H，GIRI A P，GUPTA V S. Data on changes in the fatty acid composition during fruit development and ripening of three mango cultivars （Alphonso，Pairi and Kent） varying in lactone content[J]. Data in Brief，2016，9：480-491.

[45] WADE N L. Membrane lipid composition and tissue leakage of pre- and early-climacteric banana fruit[J]. Postharvest Biology and Technology，1995，5（1/2）：139-147.

[46] 張嘉錫. 油茶亞麻酸合成關(guān)鍵基因CoFAD7的功能分析和優(yōu)異等位變異挖掘[D]. 長(zhǎng)沙：中南林業(yè)科技大學(xué)，2024.

ZHANG Jiaxi. Functional analysis of CoFAD7，a key gene for the linolenic acid synthesis in Camellia oleifera，and mining of excellent allelic variation[D]. Changsha：Central South University of Forestry and Technology， 2024.

[47] WANG J J，LIU H R，GAO J，HUANG Y J，ZHANG B，CHEN K S. Two ω-3 FADs are associated with peach fruit volatile formation[J]. International Journal of Molecular Sciences，2016，17（4）：464.

[48] PENG B，GU Z X，ZHOU Y F，NING Y Z，XU H Y，LI G，NI Y，SUN P P，XIE Z Q，SHI S P，DARK A，SONG Z Z. Potential role of fatty acid desaturase 2 in regulating peach aroma formation[J]. Postharvest Biology and Technology，2023，204：112473.

[49] HERNáNDEZ M L，SICARDO M D，BELAJ A，MARTíNEZ-RIVAS J M. The oleic/linoleic acid ratio in olive （Olea europaea L.） fruit mesocarp is mainly controlled by OeFAD2-2 and OeFAD2-5 genes together with the different specificity of extraplastidial acyltransferase enzymes[J]. Frontiers in Plant Science，2021，12：653997.

[50] LEE M W，PADILLA C S，GUPTA C，GALLA A，PEREIRA A，LI J M，GOGGIN F L. The FATTY ACID DESATURASE2 family in tomato contributes to primary metabolism and stress responses[J]. Plant Physiology，2020，182（2）：1083-1099.