姚雪倩，岳川，楊國一，陳丹，張冬桃，陳桂信，葉乃興

?

茶樹牻牛兒基牻牛兒基焦磷酸合成酶基因的克隆及表達分析

姚雪倩，岳川，楊國一，陳丹，張冬桃，陳桂信*，葉乃興*

福建農(nóng)林大學園藝學院中國烏龍茶協(xié)同創(chuàng)新中心，福建福州 350002

以鐵觀音芽葉為材料，驗證從茶樹轉(zhuǎn)錄組數(shù)據(jù)庫中篩選到的1條牻牛兒基牻牛兒基焦磷酸合成酶（GGDPS）的編碼全長序列cDNA。該cDNA全長1?661?bp，含有1個1?137?bp完整的開放閱讀框，命名為。該基因編碼378個氨基酸，氨基酸序列具有類異戊二烯合成酶家族的5個保守域和2個特征功能域；序列分析顯示，該cDNA與其他植物GGDPS高度保守，與三七（）的親緣關(guān)系最近。CsGGDPS屬于不穩(wěn)定、親水蛋白，可能定位到葉綠體中，不存在跨膜結(jié)構(gòu)，無信號肽，發(fā)生磷酸化的位點可能有20個；二級結(jié)構(gòu)主要由α-螺旋構(gòu)成，三級結(jié)構(gòu)與擬南芥GGPPS11匹配度最高。實時熒光定量PCR結(jié)果表明，在茶樹發(fā)育過程和不同葉位的表達量隨芽葉成熟度的增加呈上升趨勢，隨著做青過程的進行，的表達量逐漸升高；在父本黃旦、母本鐵觀音和子一代金觀音中均有表達，但表達量存在差異。

茶樹；牻牛兒基牻牛兒基焦磷酸合成酶；茶葉香氣；基因表達

茶樹（(L.) O. Kuntze）是多年生常綠木本植物。茶樹品種、茶園生態(tài)及茶葉加工工藝等是影響茶葉香氣的重要因素。萜烯類化合物作為茶葉香氣的重要組分之一，其含量占春季茶鮮葉總揮發(fā)油的51.26%，主要包含芳樟醇、香葉醇、橙花叔醇、法尼烯、萜品烯、杜松烯等物質(zhì)[1-2]。植物中萜烯類化合物合成途徑復雜，由多種酶參與其合成調(diào)控。其中，牻牛兒基牻牛兒基焦磷酸合成酶（Geranylgeranyl diphosphate synthase, GGDPS），是植物萜類物質(zhì)合成途徑中的關(guān)鍵酶之一[3-4]。植物萜類化合物合成途徑的第一階段經(jīng)甲羥戊酸途徑（Mevalonic acid, MVA）和甲基赤鮮糖醇磷酸途徑（Methyl- erythritol-phosph-ate, MEP）生成的終產(chǎn)物——異戊烯基焦磷酸（Isopenteny pyrophosphate, IPP）和其雙鍵異構(gòu)體二甲烯丙基焦磷酸酯（Dimethylallyl diphosphate, DMAPP），是合成牻牛兒基牻牛兒基焦磷酸（Geranylgeranyl diphosphate, GGDP）的關(guān)鍵底物。第二階段中，GGDPS能夠以3分子IPP與1分子DMAPP、1分子IPP和1分子法尼基焦磷酸（Famesyl pyrophosphate, FPP）或者2分子IPP和1分子牻牛兒基焦磷酸（Geranyl diphosphate, GPP）為底物，分別催化縮合生成GGDP[5-6]。GGDP是一種類異戊二烯化合物，該通用前體物在多步不同酶促反應下參與葉綠素、類胡蘿卜素、赤霉素、脫落酸、獨腳金內(nèi)酯、輔酶Q、維生素K、維生素E、多酚和多萜醇等的生物合成，同時還能調(diào)控植物代謝中“碳流”的方向和催化蛋白的異戊二烯化等[7-10]。因此，GGDPS被認為是植物體內(nèi)合成二萜、四萜和多萜類化合物等重要初級和次級代謝途徑中的關(guān)鍵酶，在植物生長、發(fā)育、繁殖和防御中扮演重要角色。盡管茶樹中萜烯類合成途徑中的多種關(guān)鍵酶基因已有報道，如乙酰輔酶A?；D(zhuǎn)移酶（Acetyl-CoA C-acetyltransferase, AACT）[11]、l-脫氧-D-木酮糖-5-磷酸還原異構(gòu)酶（1-deoxy- D- xylulose-5- phosphate reductoisomerase, DXR）和3-羥基-3-甲基戊二酰輔酶A還原酶（3-hydroxy-3- methylglutaryl coenzyme A reductase, HMGR）[12]等，但GGDPS基因及其在茶葉香氣形成中的作用未見報道。

本研究以茶樹芽葉為材料，克隆并驗證從茶樹轉(zhuǎn)錄組中篩選到的全長cDNA序列。對該基因的生物信息學特征進行了全面的分析；采用實時熒光定量PCR技術(shù)研究該基因在鐵觀音茶樹芽葉發(fā)育過程、不同葉位，以及做青過程不同階段的動態(tài)表達模式，并分析了該基因在父母本和子一代間的表達量變化，為揭示烏龍茶萜烯類香氣形成的分子機制和香氣相關(guān)基因在茶樹親子代間的差異表達模式提供參考。

1 材料與方法

1.1 試驗材料及試劑

供試材料：2016年3月中旬至5月中旬于福建農(nóng)林大學南區(qū)教學茶場采集以下4種樣品：（1）從3月中旬開始，對鐵觀音茶樹進行掛牌，每隔1?d對鐵觀音春梢芽葉發(fā)育過程進行觀察，摘取各生育期長勢一致的芽頭、第一葉、第二葉、第三葉和第四葉，用于研究在茶樹芽葉發(fā)育過程不同階段的表達模式變化；（2）以鐵觀音（母本）、黃旦（父本）和金觀音（雜交子一代）品種的一芽二葉為供試材料，用于研究在父母本和子一代中的表達量變化；（3）采摘鐵觀音春梢中小開面時期的第一葉、第二葉、第三葉和第四葉，研究在不同葉位的表達模式；（4）以鐵觀音中小開面二三葉為鮮葉原料，按閩南烏龍茶做青工藝進行做青[13]，并對鮮葉、曬青葉以及一搖后、二搖后、三搖后、殺青前各階段進行取樣。采集的樣品用錫箔紙包裹并標記，液氮速凍后移入—80℃冰箱保存?zhèn)溆?，每個樣品3次重復。

試劑：RNAprep Pure Plant Kit（Polysaccharides & Polyphenolics-rich）和Universal DNA Purification Kit購自天根生化科技有限公司；SuperScriptTMII Reverse Transcriptase購自Invitrogen公司；KlenTaq LA Polymerase Mix和DL2000 DNA Marker購自Clontech公司；DNA Polymerase High Fidelity (HiFi)、-T1 Simple Cloning Kit和-T1 Phage Resistant Chemically Competent Cell、One-Step gDNA Removal and cDNA Synthesis SuperMix購自北京全式金公司；2×PCR StarMix with Loading Dye購自普洛斯生物科技有限公司；瓊脂糖Agarose為西班牙進口試劑；50×TAE電泳緩沖液和DEPC處理水購自北京索萊寶公司；引物合成和測序由上海鉑尚生物技術(shù)有限公司完成。

1.2 試驗方法

1.2.1 供試材料總RNA的提取

取適量供試樣品，液氮研磨后，按照天根植物多糖多酚提取試劑盒的操作步驟提取茶樹總RNA。提取的總RNA用1%的瓊脂糖凝膠電泳檢測其完整性，用NanoDrop 2000超微量分光光度計測定RNA的濃度和純度，合格的RNA樣品于—80℃保存。

1.2.2 cDNA的合成

以鐵觀音一芽二葉為材料，以提取的總RNA為模板，參照SMARTTMcDNA Library Construction Kit說明書，采用SMART兩步逆轉(zhuǎn)錄法合成cDNA第一鏈，采用抑制LA-PCR合成cDNA第二鏈，即制備好普通RT-PCR模板。按照全式金的One-Step gDNA Removal and cDNA Synthesis SuperMix試劑盒方法合成的cDNA作為熒光定量PCR的模板。

1.2.3的全長驗證

在茶樹高通量全轉(zhuǎn)錄組（SRX1676258）數(shù)據(jù)庫中找到1個GGDPS基因全長序列，利用該序列設計全長引物-fl-F、-fl-R（表1），以鐵觀音一芽二葉總RNA的反轉(zhuǎn)錄第二鏈cDNA為模板，進行RT-PCR擴增。PCR反應體系：10×HiFi Buffer Ⅱ2.5?μL，2.5?mmol·L-1dNTPs 2?μL，-fl-F和-fl-R（10?μmol·L-1）各0.5?μL，HiFi Polymerase 0.125?μL，Template DNA 1?μL，DEPC處理水補至25?μL。程序：95℃預變性5?min，95℃變性30?s，55℃退火90?s，72℃延伸2?min，35個循環(huán)；72℃延伸10?min，最后保持4℃。擴增產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測，在紫外燈下切割目的片段大小的條帶，并用DNA凝膠回收試劑盒進行回收純化。新鮮PCR產(chǎn)物與-T1載體連接，將重組質(zhì)粒轉(zhuǎn)入-T1化學感受態(tài)細胞，37℃培養(yǎng)過夜。通過藍白斑篩選與菌落PCR鑒定得到陽性克隆，挑選3個陽性菌液送至上海鉑尚生物技術(shù)有限公司測序。

1.2.4 CsGGDPS的生物信息學特征分析

對所得全長cDNA序列，利用NCBI網(wǎng)站在線查找ORF和同源性檢索等；利用protparam、protscale、TMHMM-2.0、NetPhos、SignalP、TargetP軟件分別對保守區(qū)所編碼的蛋白質(zhì)進行理化性質(zhì)、親疏水性、跨膜結(jié)構(gòu)、磷酸化位點、信號肽預測分析和亞細胞定位；利用predictprotein和swissmodel兩個軟件預測分析蛋白質(zhì)的二級結(jié)構(gòu)和三級結(jié)構(gòu)；利用DNAMAN軟件進行CsGGDPS氨基酸多序列比對；利用MEGA5.0軟件構(gòu)建茶樹GGDPS與其他物種GGDPS的系統(tǒng)進化樹。

1.2.5 CsGGDPS實時熒光定量PCR表達分析

根據(jù)全長cDNA序列設計熒光定量PCR引物-QPCR-F、-QPCR-R（表1）。以茶樹為內(nèi)參基因，檢測在鐵觀音芽葉發(fā)育過程、不同葉位、做青過程，以及親本與子代間的相對表達量變化。數(shù)據(jù)結(jié)果采用Excel 2003和2–ΔΔT法進行定量分析。

表1 引物序列

2 結(jié)果與分析

2.1 茶樹CsGGDPS全長cDNA序列的驗證及序列分析

從茶樹轉(zhuǎn)錄組中獲得的1?597?bp的序列，NCBI BLASTx比對顯示它具有完整的開放閱讀框（ORF）。設計引物擴增驗證其全長ORF，PCR結(jié)果顯示擴增出的條帶大小與預期目的片段大小基本一致，約1?500?bp（圖1）。測序結(jié)果拼接后，獲得全長序列。BLASTn分析表明，它與公布的茶樹基因（登錄號：KU892076）的同源性達到99%。在鐵觀音中擴增出GGDPS基因，命名為。該序列全長1?661?bp，含有200?bp 5′非編碼區(qū)（UTR）、1?137?bp ORF（201~1?337?bp），324?bp 3′非編碼區(qū)（UTR）。

圖1 鐵觀音CsGGDPS擴增產(chǎn)物電泳圖

2.2 茶樹CsGGDPS的生物信息學分析

2.2.1 蛋白質(zhì)理化性質(zhì)預測

在線軟件預測該蛋白質(zhì)相對分子質(zhì)量為41.12?kDa，理論等電點為6.37，由378個氨基酸組成，半衰期為30?h。該蛋白不穩(wěn)定系數(shù)為41.26，整個多肽鏈中分值最低（－3.333）的氨基酸親水性最強；分值最高（2.433）的氨基酸疏水性最強（圖2），總平均親水性為－0.072，推測該蛋白為不穩(wěn)定、親水性蛋白。

圖2 CsGGDPS蛋白親/疏水性分析

2.2.2 跨膜結(jié)構(gòu)、信號肽和亞細胞定位預測

跨膜區(qū)域顯示（圖3）整條肽鏈位于膜外，故CsGGDPS不存在跨膜結(jié)構(gòu)，非膜蛋白，可能是在游離核糖體上合成蛋白多肽。TargetP1.1 Server在線分析含有信號肽（Signal peptide, SP）的可能性分值為0.007，說明CsGGDPS不具有信號肽，屬于非分泌蛋白，該結(jié)果與SignalP-HMM預測結(jié)果一致。亞細胞定位預測顯示，全氨基酸序列中葉綠體運輸肽（Chloroplast transit peptide, cTP）的可能性分值最高為0.764，線粒體靶向肽（Mitochondrial targeting peptide, mTP）的可能性分值為0.247，其他可能性分值為0.051，可靠等級值為3，由此推斷CsGGDPS蛋白定位到葉綠體上。在細胞質(zhì)中合成GGDPS，再在葉綠體轉(zhuǎn)運肽的協(xié)助下進入葉綠體，于葉綠體基質(zhì)中與底物結(jié)合催化合成GGDP。

圖3 CsGGDPS蛋白跨膜區(qū)域分析

2.2.3 CsGGDPS蛋白的二級結(jié)構(gòu)和三級結(jié)構(gòu)預測

CsGGDPS蛋白的二級結(jié)構(gòu)主要由α-螺旋構(gòu)成，占47.62%，其次為不規(guī)則卷曲，占30.42%，N端含有較多不規(guī)則卷曲，可能存在活性位點，延伸鏈和β-轉(zhuǎn)角同樣分布于整個肽鏈中，分別占15.08%和6.88%。CsGGDPS的三級結(jié)構(gòu)建模與二級結(jié)構(gòu)預測類似（圖4），預測結(jié)果評估顯示GMQE分值為0.71，說明構(gòu)建的CsGGDPS蛋白三級結(jié)構(gòu)可信度高。建模結(jié)果顯示，CsGGDPS為同源二聚體蛋白，與擬南芥的GGPPS11（5e8l.1.B）匹配度最高（79.39%），因此借助模式植物擬南芥該基因的表達情況可為深入研究奠定基礎。

圖4 CsGGDPS蛋白三級結(jié)構(gòu)分析

2.2.4 磷酸化位點預測

通過NetPhos 2.0 Serve對CsGGDPS進行磷酸化結(jié)構(gòu)的預測。結(jié)果表明，絲氨酸的磷酸化位點可能有14個，位于多肽鏈的第23、25、60、97、107、215、216、229、244、249、282、313、340、356位氨基酸，其中第107、215、216、313、356位氨基酸的預測值均達到0.980以上，說明這些位點的絲氨酸發(fā)生磷酸化的可能性較高。蘇氨酸可能的磷酸化位點有6個，分別為第26、48、66、76、211、264位氨基酸。

2.2.5 同源性比對和分子系統(tǒng)進化樹分析

將CsGGDPS的氨基酸序列進行BLASTp比對，結(jié)果顯示（圖5），該氨基酸序列屬于類異戊二烯家族，含有底物結(jié)合袋（Substrate binding pocket）、底物-Mg2+結(jié)合位點（Substrate-Mg2+binding site）、鏈長決定區(qū)（Chain length determination region）、催化殘基的活性位點，以及3個特異結(jié)構(gòu)域“Trans_ IPPS_HT”、“IspA”和“Polypreny_synt”等。CsGGDPS與多個物種中的GGDPS具有較高的相似性，相似性均在70%以上，其中與三七（）、長春花（）、芝麻（）、米團花（）、葡萄（）的相似度分別為77%、75%、74%、73%、71%；將它與8個物種GGDPS氨基酸序列進行多重序列比對，結(jié)合功能結(jié)構(gòu)域分析，結(jié)果顯示（圖6），這些植物的GGDPS基因都包含聚異戊二烯合成酶的5個特征性功能結(jié)構(gòu)域（ⅠⅤ），同時含有2個幾乎所有萜類合成酶共有的保守區(qū)天冬氨酸富集基序DDXXXXD和DDXXD（X可為任意氨基酸）。CsGGDPS在5′端前86個氨基酸位點保守性較低，從第87位開始到之后的位點氨基酸保守性較高。

將19個物種的GGDPS氨基酸序列采用鄰近連接法（Neighbor-Joining, NJ）構(gòu)建系統(tǒng)進化樹（圖7）。GGDPS在分子進化樹上被分為兩大類，其中CsGGDPS與三七（）的聚在一起，親緣關(guān)系最近，同時與長春花（）和龍膽（）的親緣關(guān)系也較近；而與木本植物，如蘋果、白梨、胡楊等的關(guān)系較遠。

2.3 CsGGDPS相對熒光定量PCR表達分析

采用實時熒光定量PCR技術(shù)，檢測在鐵觀音芽葉發(fā)育過程不同階段、不同葉位和做青過程的表達量變化。結(jié)果顯示，在發(fā)育過程中，的表達量總體呈上升趨勢，以芽頭為對照，一葉時表達量約為芽頭的2倍，從二葉到四葉過程中，的表達量呈明顯上升趨勢，四葉期的表達量最高，約為芽頭的4倍（圖8-A）。在不同葉位中，從第一葉到第四葉，的表達量呈上升趨勢，其中第三葉的表達量達到最高，約為第一葉的2倍，第四葉中仍保持較高的表達量水平（圖8-B）。在做青過程中，從鮮葉到曬青葉，的表達量較低且增長緩慢；一搖后表達量增加，與鮮葉相比，二者相差近6倍；隨著做青過程的進行，從二搖后到殺青前，的表達量顯著升高，殺青前的表達量達到最高，約為鮮葉的20倍（圖8-C）。

研究檢測了在鐵觀音母本（♀）、黃旦父本（♂）及它們的雜交子一代金觀音3個品種間的表達水平。結(jié)果顯示（圖9），在親本黃旦、鐵觀音和子一代金觀音中均有表達，屬于基因差異表達模式中的共同表達型。然而，在鐵觀音中的表達量顯著高于黃旦，它們的雜交后代金觀音中的表達水平介于父母本之間，呈現(xiàn)偏低親的雜種表達量現(xiàn)象，這可能與3個品種間香氣差異有關(guān)。

圖9 茶樹親本與子一代間CsGGDPS表達分析

3 討論

茶樹中芳香物質(zhì)代謝包括不飽和脂肪族醇、芳香族醇和萜烯類化合物3個方面的生物合成與轉(zhuǎn)化[14]。本研究克隆了鐵觀音中CsGGDPS基因，該基因編碼蛋白是植物萜烯類化合物合成途徑中的關(guān)鍵酶之一。CsGGDPS基因cDNA全長1?661?bp，與其他植物的具有較高的相似性，氨基酸多重序列比對顯示該序列含5個保守結(jié)構(gòu)域以及GGDPS的特異功能域——FARM（the first aspartate-rich motif）和SARM（the second aspartate-rich motif），這2個基序是GGDPS關(guān)鍵的催化活性位點，對Mg2+橋與底物分子二磷酸基團的結(jié)合至關(guān)重要[15-16]。

本研究結(jié)果表明，在鐵觀音芽葉發(fā)育過程和不同葉位的相對表達量總體隨芽葉成熟度的增加而呈上升的趨勢，與劉晶晶等[17]研究萜類化合物在茶樹發(fā)育階段的合成與釋放大體表現(xiàn)為幼葉中多、老葉中少的結(jié)果一致。張冬桃[18]分析了萜類代謝途徑的關(guān)鍵基因，如等在鐵觀音發(fā)育過程和不同葉位的表達量情況，結(jié)果表明這些基因的表達均隨芽葉成熟進程呈上升趨勢。從CsGGDPS的分子系統(tǒng)進化樹可以看出，它與三七聚為一類。閔丹丹[4]研究發(fā)現(xiàn)，三七可能與三七葉綠素、類胡蘿卜素等二萜物質(zhì)的合成代謝及葉綠體的發(fā)育有關(guān)，由此可推測可能與三七的功能相近，伴隨芽葉成熟度的增加，茶樹體內(nèi)不斷合成葉綠素和類胡蘿卜素等供茶樹生長發(fā)育必需的物質(zhì)，萜類前體物質(zhì)也逐步積累[19]，最終導致定位于葉綠體中的萜類合成酶的表達量上升。這與烏龍茶生產(chǎn)加工需采摘一定成熟度芽葉的理論相吻合。雖然前人有研究表明，未受損傷的茶樹鮮葉基本無味或味輕微[20]，但以中小開面為標準采摘的烏龍茶鮮葉中能檢測到更多微量的游離態(tài)萜類香氣物質(zhì)，為烏龍茶制茶過程中香氣品質(zhì)的形成與轉(zhuǎn)化提供物質(zhì)基礎。做青是烏龍茶獨特香氣形成的關(guān)鍵工序，其中的搖青工序使茶鮮葉破損，細胞組織機械損傷，加速萜烯糖苷的水解。β-葡萄糖苷酶在水解萜烯類香氣前驅(qū)體中發(fā)揮重要作用[21]。張秀云等[22]研究發(fā)現(xiàn)，β-葡萄糖苷酶在做青過程中呈峰型變化，做青結(jié)束前2?h活性達到最高，糖苷鍵合態(tài)香氣物質(zhì)水解為游離態(tài)的能力最強，促使己烯酯類、芳樟醇氧化物、倍半萜烯類、順-茉莉酮、茉莉內(nèi)酯和苯乙醛等芳香物質(zhì)大量形成[23]。做青過程中茶鮮葉內(nèi)部氧氣、水分、膜透性發(fā)生變化，搖青使脂肪酸氧化裂解生成的青葉醇和青葉醛大部分揮發(fā)，留下少量具有清香的反式青葉醇[24]。烏龍茶做青過程中還伴隨胡蘿卜素類的氧化降解產(chǎn)物生成，如β-紫羅酮、二氫海葵內(nèi)酯、茶螺烯酮等香氣物質(zhì)[25]。烏龍茶品質(zhì)特征主要形成于三搖或其后階段，橙花叔醇、香葉醇、芳樟醇及其氧化物、吲哚、α-法尼烯、順-茉莉酮等香氣組分在一定搖青強度范圍內(nèi)不斷積累[26-29]。本研究發(fā)現(xiàn)，在鐵觀音做青過程中表達量顯著增加，殺青前其表達量最高可達到鮮葉的20倍。由前面的研究結(jié)果分析，3次搖青靜置后的做青葉，具有典型的蘭花香，對烏龍茶品質(zhì)具有重要的影響。因此，由在鐵觀音芽葉發(fā)育過程、不同葉位和做青過程的表達結(jié)果推測，該基因與烏龍茶萜烯類香氣品質(zhì)的形成密切相關(guān)。

在母本鐵觀音、父本黃旦和子一代金觀音中的表達差異較大。金觀音中的表達量接近父本而顯著低于母本，屬于基因表達量差異中的偏低親表達模式，也稱作非加性表達[30]。這與Li等[31]提出的關(guān)于雙親中不同表達量的基因比相同表達量的基因更易呈現(xiàn)非加性表達模式的觀點相近，如苜蓿中這一比值為3∶1。郭吉春等[32]對金觀音的選育研究認為，金觀音的遺傳性狀趨向母本，該觀點與本實驗在金觀音和鐵觀音中表達量有明顯差異的結(jié)果不一致。因此推測，雙親和子一代中表達的差異既有可能受順式作用元件、反式作用元件或二者共同作用的影響，也有可能與它們的基因表達調(diào)控、表觀遺傳機制有關(guān)。鐵觀音和黃旦作為重要的烏龍茶育種材料，通過雜交選育了一批優(yōu)良的茶樹品種，豐富了烏龍茶品種資源，對烏龍茶的發(fā)展產(chǎn)生了深遠的影響。但由于茶樹是異花授粉的遺傳雜合體，親本與F1代間在一定程度上存在遺傳差異[33]，所以鐵觀音和黃旦雜交后基因的遺傳規(guī)律還有待后續(xù)研究。茶葉香氣這一品質(zhì)性狀屬于數(shù)量性狀，數(shù)量性狀從微效多基因控制過渡到主基因-多基因混合控制形成植物雜種優(yōu)勢[34]。因此，萜類合成酶相關(guān)基因在雜種和親本中的遺傳趨勢有待進一步研究驗證，以期為雜種的早期選擇、雜種優(yōu)勢形成理論和烏龍茶加工過程中萜烯類香氣形成等的解析提供理論依據(jù)。

[1] 張正竹, 施兆鵬, 宛曉春. 萜類物質(zhì)與茶葉香氣（綜述）[J]. 安徽農(nóng)業(yè)大學學報, 2000, 27 (1): 51-54.

[2] 吳勇. 萜烯類化合物與茶葉香氣[J]. 化學工程與裝備, 2009, 32(11): 123-125.

[3] 魏攀, 孟利軍, 陳千思, 等. 煙草牻牛兒基牻牛兒基焦磷酸合成酶基因的克隆和功能分析[J]. 煙草科技, 2016, 49(4): 8-15.

[4] 閔丹丹, 唐美瓊, 李剛, 等. 三七香葉基香葉基焦磷酸合酶基因的克隆及表達分析[J]. 中國中藥雜志, 2015, 40(11): 2090-2095.

[5] 楊瀅. 植物萜類代謝途徑關(guān)鍵酶的比較及其在丹參中的表達分析[D]. 西安: 陜西師范大學, 2012: 1-2.

[6] Beck G, Coman D, Herren E, et al. Characterization of thesynthase gene family in[J]. Plant Molecular Biology, 2013, 82(4): 393-416.

[7] Yamamura Y, Mizuguchi Y, Taura F, et al. Transcriptional activation of a geranylgeranyl diphosphate synthase gene，，isolated fromby treatment with methyl jasmonate and yeast extract [J]. Journal of Natural Medicines, 2014, 68(4): 748-753.

[8] 李澤鋒, 魏攀, 夏玉珍, 等. 煙草牻牛兒基牻牛兒基焦磷酸合成酶基因家族的全基因組鑒定[J]. 煙草科技, 2015, 48(6): 1-8.

[9] 張萌, 蘇平, 劉雨佳, 等. 雷公藤牻牛兒基牻牛兒基焦磷酸合酶基因全長cDNA的獲得及生物信息學分析[J]. 中國中藥雜志, 2015, 40(6): 1066-1070.

[10] 韓立敏. 菘藍牦牛兒基牦牛兒基焦磷酸合成酶基因() 的克隆及其表達特性分析[J]. 基因組學與應用生物學, 2015, 34( 6): 1172-1178.

[11] Wei X, Yan X, Fu M W, et al. In silico analysis and feeding assays of some genes in the early steps of terpenoid biosynthetic pathway in[J]. Journal of Tea, 2013, 39(4): 191-198.

[12] 徐燕. 茶樹萜類合成途徑關(guān)鍵基因克隆及表達研究[D]. 杭州: 浙江大學, 2013: 12.

[13] 溫立香, 郭雅玲, 黃壽輝. 烏龍茶做青技術(shù)的研究進展[J]. 安徽農(nóng)業(yè)科學, 2015, 43(24): 215-217.

[14] 宛曉春. 茶葉生物化學[M]. 3版. 北京: 中國農(nóng)業(yè)出版社, 2003: 119-120.

[15] HEMMI H, NOIKE M, NAKAYAMA T, et al. An alternative mechanism of product chain-length determination in type III granylgeranyl diphosphate synthase [J]. Eur J Biochem, 2003, 270(10): 2186-2194.

[16] Song L, Poulter C D. Yeast farnesyl-diphosphate synthase: site-directed mutagenesis of residues in highly conserved prenyltransferase domainsⅠand Ⅱ[J]. Proc Natl Acad Sci USA, 1994, 91(8): 3044-3048.

[17] 劉晶晶, 王富民, 劉國峰, 等. 茶樹萜類香氣物質(zhì)代謝譜與相關(guān)基因表達譜時空變化的關(guān)系[J]. 園藝學報, 2014, 41(10): 2094-2106.

[18] 張冬桃. 烏龍茶做青過程香氣形成相關(guān)基因的分離與表達[D]. 福州: 福建農(nóng)林大學, 2016: 54.

[19] Dudareva N, Pichersky E. Biology of floral scent[M]. Boca Raton: CRC Press, 2006: 55-78.

[20] 許寧, 陳宗懋, 游小清. 引誘茶尺蠖天敵寄生蜂的茶樹揮發(fā)物的分離與鑒定[J]. 昆蟲學報, 1999, 42(2): 126-131.

[21] 吳穎, 戴永峰, 張凌云. 做青工藝對烏龍茶品質(zhì)影響研究進展[J]. 廣東茶業(yè), 2013, 31(5): 8-11.

[22] 張秀云, 方世輝, 夏濤. 烏龍茶萎凋做青中β-葡萄糖苷酶活性變化研究[J]. 安徽農(nóng)業(yè)大學學報, 2000, 27(2): 164-166.

[23] 楊意成, 梁月榮. 烏龍茶花香形成機理的研究[J]. 茶葉, 2008, 34(1): 10-14.

[24] 梁曉嵐, 陳春林. 烏龍茶香氣形成機理初探[J]. 廣東農(nóng)業(yè)科學, 1996, 23(4): 23-25.

[25] Yang Z Y, Baldermann S, Watanabe N. Recent studies of the volatile compounds in tea [J]. Food Research International, 2013, 53(2): 585-599.

[26] 黃福平, 陳榮冰, 梁月榮, 等. 烏龍茶做青過程中香氣組成的動態(tài)變化及其與品質(zhì)的關(guān)系[J]. 茶葉科學, 2003, 23(1): 31-37.

[27] 陳林, 陳鍵, 陳泉賓, 等. 做青工藝對烏龍茶香氣組成化學模式的影響[J]. 茶葉科學, 2014, 34(4): 387-395.

[28] 王日為, 張麗霞, 楊偉麗, 等. 烏龍茶做青過程中香氣的動態(tài)變化規(guī)律[J]. 湖南農(nóng)業(yè)大學學報, 1999, 25(3): 194-199.

[29] Lin J, Zhang P, Pan Z Q, et al. Discrimination of Oolong tea () varieties based on feature extraction and selection from aromatic profiles analysed by HS-SPME/ GC-MS [J]. Food Chemistry, 2013, 141(1): 259-265.

[30] 許晨璐, 孫曉梅, 張守攻. 基因差異表達與雜種優(yōu)勢形成機制探討[J]. 遺傳, 2013, 35(6): 714-726.

[31] Li X H, Wei Y L, Nettleton D, et al. Comparative gene expression profiles between heterotic and non-heterotic hybrids of tetraploid[J]. BMC Plant Biol, 2009, 9(1): 107-115.

[32] 郭吉春, 葉乃興, 楊如興, 等. 茶樹雜交種金觀音、黃觀音的選育與應用[J]. 福建茶葉, 2008, 31(1): 11-14.

[33] 王讓劍, 楊軍, 孔祥瑞, 等. 福建15個茶樹品種SSR遺傳差異分析與指紋圖譜建立[J]. 福建農(nóng)業(yè)學報, 2014, 29(10): 970-975.

[34] 蓋鈞鎰, 章元明, 王建康. 植物數(shù)量性狀遺傳體系[M]. 北京: 科學出版社, 2003: 145-156.

Cloning and Expression Analysis of Geranylgeranyl Diphosphate Synthase Genein Tea Plant ()

YAO Xueqian, YUE Chuan, YANG Guoyi, CHEN Dan, ZHANG Dongtao, CHEN Guixin*, YE Naixing*

College of Horticulture, Fujian Agriculture and Forestry University, Collaborative Innovation Center of Chinese Oolong Tea Industry, Fuzhou 350002, China

A full-length cDNA sequence encoding geranylgeranyl diphosphate synthase (GGDPS) was isolated from transcriptome database of tea plant, cloned fromTieguanyin and named as. The cDNA length ofwas 1?661?bp, with an open reading frame (ORF) of 1?137?bp and deduced protein of 378 amino acids. The protein was deduced to contain 5 conserved domains with 2 functional domains of Isoprenoid-Biosyn-C1 superfamily. The sequence analysis showed that CsGGDPS was highly conserved and had the closest genetic relationship with. CsGGDPS was an instability and hydrophilic protein, which was predicted to be located in chloroplast but with no transmembrane structure and signal peptide. There were 20 phosphorylation sites within the polypeptide chain. Alpha helix was predicted to be the major secondary structure of CsGGDPS. The three-dimension structure of CsGGDPS was highly similar to GGPPS11 from. The quantitative real-time PCR showed thatexpression was increased during the developmental process and increased with the age of tea leaves. Meanwhile, its expression was also enhanced during the Zuoqing procedure.was ubiquitously expressed in theHuangdan,Tieguanyin and their first filial generationJinguanyin, but with different expression levels.

, GGDPS, tea aroma, gene expression

S571.1；Q946.5

A

1000-369X(2017)01-086-11

2016-08-01

2016-08-23

國家自然科學基金（31600555）、福建省“2011協(xié)同創(chuàng)新中心”中國烏龍茶產(chǎn)業(yè)協(xié)同創(chuàng)新中心專項（閩教科〔2015〕75號）、福建茶產(chǎn)業(yè)農(nóng)技推廣服務試點建設（KNJ-151001）。

姚雪倩，女，碩士研究生，主要從事茶樹栽培育種與生物技術(shù)研究。