張曉東+李彩霞+王連春+王元忠

摘要：根據(jù)滇龍膽轉(zhuǎn)錄組乙酰CoA?；D(zhuǎn)移酶基因GrAACT序列設(shè)計(jì)引物，使用逆轉(zhuǎn)錄PCR（RT-PCR）技術(shù)從滇龍膽幼葉擴(kuò)增該基因，并進(jìn)行序列分析、原核表達(dá)和組織特異性分析。結(jié)果表明：GrAACT基因（登錄號(hào)KJ917163）開(kāi)放閱讀框長(zhǎng) 1 215 bp，編碼404個(gè)氨基酸；GrAACT蛋白相對(duì)分子質(zhì)量41.41 ku，pI值為6.25，屬于AACT蛋白家族成員，無(wú)信號(hào)肽，為疏水穩(wěn)定蛋白，主要由α-螺旋、無(wú)規(guī)則卷曲構(gòu)成；GrAACT蛋白具有AACT蛋白保守結(jié)構(gòu)域：硫解酶N端結(jié)構(gòu)域（第12～272位氨基酸）、硫解酶C端結(jié)構(gòu)域（第282～402位氨基酸）、類硫解酶結(jié)構(gòu)域（第13～403位氨基酸）；在硫解酶C端結(jié)構(gòu)域中，包含硫解酶保守位點(diǎn)（第350～366位氨基酸）、硫解酶活性位點(diǎn)（第385～398位氨基酸）；滇龍膽GrAACT蛋白與長(zhǎng)春花CrAACT蛋白親緣關(guān)系最近；GrAACT基因在大腸桿菌中表達(dá)的重組蛋白相對(duì)分子質(zhì)量約為6741 ku；組織特異性表達(dá)分析結(jié)果表明，GrAACT基因主要在莖、根中表達(dá)。

關(guān)鍵詞：滇龍膽；乙酰CoA酰基轉(zhuǎn)移酶；基因克?。簧镄畔W(xué)分析；表達(dá)分析

中圖分類號(hào)： S184；R282.71；Q785文獻(xiàn)標(biāo)志碼： A文章編號(hào)：1002-1302（2016）06-0094-05

收稿日期：2015-05-27

基金項(xiàng)目：國(guó)家自然科學(xué)基金（編號(hào)：81260608）；國(guó)家科技支撐計(jì)劃（編號(hào)：2011BAI13B02-04）；云南省教育廳科學(xué)研究基金重點(diǎn)項(xiàng)目（編號(hào)：2015Z171）。

作者簡(jiǎn)介：張曉東（1980—），男，河南新鄭人，博士，副教授，從事植物代謝基因工程研究。E-mail：zxd95@126.com。

通信作者：王元忠，碩士，助理研究員，從事藥用植物資源評(píng)價(jià)與利用的研究。E-mail：boletus@126.com。滇龍膽是傳統(tǒng)中藥龍膽的來(lái)源植物之一，也是云南省重要道地藥材，主要分布于中國(guó)西南部，尤其是云南省山區(qū)地帶；其根入藥，用于治療各種炎癥、肝炎、風(fēng)濕病和膽囊炎等[1]。目前，醫(yī)用龍膽需求量每年遞增約10%，而現(xiàn)有滇龍膽產(chǎn)量低、品質(zhì)不高，且野生滇龍膽資源遭到人為大肆采挖[2]。為選育優(yōu)質(zhì)、高產(chǎn)、高抗新品種，2013年6月云南省已啟動(dòng)龍膽草航天育種工程[3]。從中醫(yī)藥現(xiàn)代化角度來(lái)看，滇龍膽主要藥效成分為龍膽苦苷，要從根本上解決龍膽草藥源問(wèn)題并保護(hù)野生滇龍膽，除篩選新品種之外，還必須弄清龍膽苦苷的生物合成途徑及其調(diào)控機(jī)理，為通過(guò)基因工程手段生產(chǎn)龍膽苦苷奠定基礎(chǔ)。

龍膽苦苷為裂環(huán)烯醚萜類化合物，在植物中主要通過(guò)質(zhì)體2-C-甲基-D-赤蘚糖醇-4-磷酸（MEP）、胞質(zhì)甲羥戊酸（MVA）途徑合成[4]。乙酰CoA?；D(zhuǎn)移酶（又稱乙酰乙酰CoA硫解酶，AACT，EC 2.3.1.9）是MVA途徑中的第1個(gè)限速酶，能夠催化2分子乙酰CoA形成1分子乙酰乙酰CoA[5]。目前，乙酰CoA?；D(zhuǎn)移酶基因AACT已從擬南芥[5]、丹參[6]、胡蘿卜[7]、白木香[8]、油茶[9]等許多植物中分離。AACT基因表達(dá)具有組織特異性，并被生物、非生物因素誘導(dǎo)。擬南芥中存在2個(gè)平行進(jìn)化的同源基因，其中AtAACT1主要在維管系統(tǒng)中表達(dá)，AtAACT2則在根尖、幼葉、莖上部和花藥中大量表達(dá)，二者存在功能冗余；AtAACT2的作用是為胞質(zhì)定位、MVA起源的異戊烯萜類生物合成途徑合成大量乙酰乙酰CoA前體[5]。在生物誘導(dǎo)劑（100 mg/mL酵母提取物）、非生物誘導(dǎo)劑（30 mmol/L Ag+）共同誘導(dǎo)下，丹參SmAACT基因在誘導(dǎo)后12 h的表達(dá)量達(dá)到最高值[10]。

目前，由于龍膽基因組資源極度匱乏[11]，導(dǎo)致龍膽苦苷生物合成途徑并不清楚[3]。本研究根據(jù)筆者實(shí)驗(yàn)室滇龍膽轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中的GrAACT基因序列，設(shè)計(jì)特異性引物，通過(guò)逆轉(zhuǎn)錄PCR（RT-PCR）技術(shù)成功從滇龍膽幼葉中擴(kuò)增到GrAACT基因，并進(jìn)行序列分析、原核表達(dá)和組織表達(dá)特異性分析，以期為闡明滇龍膽龍膽苦苷生物合成途徑及其調(diào)控機(jī)理奠定基礎(chǔ)。

1材料與方法

1.1材料

滇龍膽組培苗、盆栽苗均采自玉溪師范學(xué)院分子生物學(xué)實(shí)驗(yàn)室，經(jīng)云南省農(nóng)業(yè)科學(xué)院藥用植物研究所金航研究員鑒定為滇龍膽（Gentiana rigescens Franch. ex Hemsl）?；蚩寺∷貌牧蠟榈猃埬憻o(wú)菌苗幼葉，基因組織特異性表達(dá)分析所用材料為盆栽3年生滇龍膽的根、莖、葉，采樣日期為2014年5月17日。

1.2方法

按照多糖植物組織提取試劑RNAiso（TaKaRa，大連）說(shuō)明書提取滇龍膽幼葉總RNA；按照逆轉(zhuǎn)錄試劑盒（TaKaRa，大連）說(shuō)明書合成cDNA。根據(jù)原核表達(dá)載體pGEX-4T-1多克隆酶切位點(diǎn)和筆者所在實(shí)驗(yàn)室前期測(cè)序的滇龍膽轉(zhuǎn)錄組GrAACT基因序列，設(shè)計(jì)1對(duì)特異引物GrAACT BamHⅠ-F、GrAACT XhoⅠ-R（表1）。以cDNA為模板進(jìn)行PCR擴(kuò)增，反應(yīng)體系為：0.5 μL TransTaq Taq DNA Polymerase High Fidelity（2.5 U/μL，北京全式金生物技術(shù)有限公司），5 μL 10×TransTaq HiFi Buffer Ⅱ，4 μL dNTP（2.5 mmol/L，TaKaRa，大連），2 μL模板，各1 μL正反向引物（10 μmol/L），加dd H2O補(bǔ)足至50 μL。PCR反應(yīng)條件為：94 ℃ 3 min；94 ℃ 30 s，55.5 ℃ 30 s，72 ℃ 75 s，30個(gè)循環(huán)；72 ℃ 10 min。PCR產(chǎn)物經(jīng)1.0%瓊脂糖凝膠電泳分離，然后割膠，使用膠回收試劑盒（Qiagen，德國(guó)）對(duì)目的片段進(jìn)行回收，將其連接到pMD19-T載體（TaKaRa，大連）。轉(zhuǎn)化大腸桿菌DH5α（TaKaRa，大連）后進(jìn)行藍(lán)白斑篩選，挑取12個(gè)白斑搖菌；使用試劑盒提取質(zhì)粒（北京百泰克生物技術(shù)有限公司），經(jīng)酶切檢測(cè)正確后進(jìn)行測(cè)序[生工生物工程（上海）股份有限公司]，獲得重組質(zhì)粒pMD19-GrAACT。

1.2.1GrAACT基因的生物信息學(xué)分析用NCBI網(wǎng)站的BLAST程序進(jìn)行序列比對(duì)，用Genetyx 6.1.8進(jìn)行翻譯，用DNAMAN 7進(jìn)行多序列比對(duì)，用Clustal X2.1進(jìn)行比對(duì)，然后用MEGA 6.0軟件內(nèi)置的NJ法構(gòu)建系統(tǒng)進(jìn)化樹(shù)，設(shè)置Bootstrap=1 000；利用在線數(shù)據(jù)庫(kù)（http：//molbiol. edu.ru/eng/scripts/01_11.html）進(jìn)行稀有密碼子分析。用ChloroP服務(wù)器v1.1進(jìn)行葉綠體轉(zhuǎn)運(yùn)肽預(yù)測(cè)；用Interpro軟件進(jìn)行保守結(jié)構(gòu)域預(yù)測(cè)；用ProtScalee軟件進(jìn)行疏水性分析；用PredictProtein對(duì)二級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)；用Phyre2（http：//www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi？id=index）對(duì)三級(jí)結(jié)構(gòu)進(jìn)行預(yù)測(cè)；用Expasy中的TMHMM工具預(yù)測(cè)蛋白的跨膜螺旋區(qū)；用在線工具WOLF PSORT預(yù)測(cè)蛋白的亞細(xì)胞定位情況。

1.2.2GrAACT基因的熒光實(shí)時(shí)定量分析分別取3年生滇龍膽的根、莖、葉，提取總RNA，用DNase I處理除去基因組DNA。用反轉(zhuǎn)錄試劑盒（TaKaRa）合成第1鏈cDNA。以轉(zhuǎn)錄組中GrGAPDH基因（GenBank登錄號(hào)KM061807）作為內(nèi)參設(shè)計(jì)特異性引物（表1），PCR條件：95 ℃ 3 min，95 ℃ 15 s，60 ℃ 31 s。根據(jù)GrAACT基因開(kāi)放閱讀柜（ORF）序列設(shè)計(jì)特異性引物（表1），使用SuperReal PreMix Plus試劑盒[天根生化科技（北京）有限公司]進(jìn)行qPCR，擴(kuò)增條件為：95 ℃ 3 min，95 ℃ 15 s，60 ℃ 31 s；每個(gè)反應(yīng)重復(fù)3次。反應(yīng)在ABI7000熒光定量PCR儀（Applied Biosystems）上進(jìn)行擴(kuò)增，擴(kuò)增曲線、溶解曲線、標(biāo)準(zhǔn)曲線由定量PCR儀軟件自動(dòng)生成。使用內(nèi)參基因GrGAPDH表達(dá)校準(zhǔn)后，計(jì)算根、莖、葉中GrAACT基因的相對(duì)表達(dá)量。采用比較CT值的“2-ΔΔCT”的方法進(jìn)行定量數(shù)據(jù)的分析處理。

2結(jié)果與分析

2.1滇龍膽GrAACT基因的克隆

以滇龍膽幼葉cDNA為模板，使用表1中的基因特異性引物GrAACTBamHI-F、GrAACTBamHI-R擴(kuò)增出約 1 200 bp 的片段，詳見(jiàn)圖1。通過(guò)TA克隆獲得重組質(zhì)粒 pMD19-GrAACT。

2.2GrAACT基因的生物信息學(xué)分析

2.2.1GrAACT基因的ORF分析和多序列比對(duì)分析利用Genetyx軟件對(duì)GrAACT基因的ORF序列進(jìn)行分析，結(jié)果顯示：該基因（GenBank登錄號(hào)KJ91716）長(zhǎng)1 215 bp，編碼404個(gè)氨基酸。利用GenBank數(shù)據(jù)庫(kù)的BLASTp程序?qū)rAACT蛋白進(jìn)行比對(duì)分析，結(jié)果表明：滇龍膽GrAACT蛋白與長(zhǎng)春花CrAACT蛋白的序列相似性最高（88.61%），與單子葉植物玉米ZmAACT蛋白的相似性（78.55%）稍低。利用DNAMAN 7將GrAACT蛋白序列與NCBI中相似性較高的部分序列進(jìn)行多序列比對(duì)分析，結(jié)果表明：GrAACT蛋白與已知蛋白序列相比保守性較高（圖2）。利用Mega 6.0將GrAACT氨基酸序列與NCBI中相似性較高的部分序列進(jìn)行系統(tǒng)發(fā)育分析，結(jié)果顯示：滇龍膽GrAACT蛋白與長(zhǎng)春花CrAACT蛋白處于同一進(jìn)化支（圖3），表明二者親緣關(guān)系較近。

2.2.2GrAACT蛋白的理化特性分析使用ExPASy ProtParam tool對(duì)GrAACT蛋白進(jìn)行分析，結(jié)果表明：GrAACT蛋白單體相對(duì)分子質(zhì)量為41 410 u，pI值為6.25，這與丹參SmAACT蛋白[6]、油茶CoAACT蛋白[9]、擬南芥ATAACT2蛋白[5]類似；帶正電氨基酸殘基（Arg+Lys）數(shù)為36個(gè)，帶負(fù)電氨基酸殘基（Asp+Glu）數(shù)為38個(gè)，化學(xué)式為：C1812H2959N509O562S17；不穩(wěn)定指數(shù)為27.92，屬穩(wěn)定蛋白；脂肪指數(shù)為9856，總平均疏水性（GRAVY）為0.178，為疏水蛋白（圖4）。GrAACT蛋白含20種基本氨基酸，其中丙氨酸含量最高，為13.9%；其次是甘氨酸、纈氨酸，含量分別為11.9%、9.9%；色氨酸含量最低，為0.5%。

2.2.3GrAACT蛋白的二級(jí)結(jié)構(gòu)、三級(jí)結(jié)構(gòu)分析利用SSpro方法對(duì)GrAACT蛋白進(jìn)行二級(jí)結(jié)構(gòu)分析，結(jié)果表明：該蛋白二級(jí)結(jié)構(gòu)中α-螺旋（H）占39.11%，無(wú)規(guī)則卷曲（C）占40.10%，延伸帶（E）占20.79%。利用Swiss-Model Workspace的自動(dòng)模式預(yù)測(cè)GrAACT蛋白的三級(jí)結(jié)構(gòu)，結(jié)果見(jiàn)圖5，該模型以人線粒體乙酰CoA?；D(zhuǎn)移酶[2f2s.1.C]為模板，在第13—404位氨基酸處建模，序列相似度為52.73%。

2.2.4GrAACT蛋白保守結(jié)構(gòu)域分析使用InterPro在線工具對(duì)GrAACT蛋白的保守結(jié)構(gòu)域進(jìn)行分析，結(jié)果表明：GrAACT

蛋白包含3個(gè)保守結(jié)構(gòu)域，它們分別為硫解酶N端結(jié)構(gòu)域（IPR020616，第12—272位氨基酸）、硫解酶C端結(jié)構(gòu)域（IPR020617，第282—402位氨基酸）、類硫解酶結(jié)構(gòu)域（IPR016039，第13—403位氨基酸）（圖6）。在硫解酶C端結(jié)構(gòu)域中，包含硫解酶保守位點(diǎn)（THIOLASE_2，第350—366位氨基酸）、硫解酶活性位點(diǎn)（THIOLASE_3，第385—398位氨基酸）。

2.2.5GrAACT蛋白信號(hào)肽分析利用ExPASy SignalP 4.1服務(wù)器分析GrAACT蛋白，未發(fā)現(xiàn)信號(hào)肽，表明該蛋白為非分泌型蛋白。利用TMHMM工具預(yù)測(cè)GrAACT蛋白的跨膜螺旋區(qū)，結(jié)果表明：GrAACT蛋白不含有跨膜區(qū)域（圖7）。

2.2.6GrAACT蛋白亞細(xì)胞定位分析使用WoLF PSORT軟件進(jìn)行亞細(xì)胞定位分析，結(jié)果表明，GrAACT蛋白在細(xì)胞質(zhì)、葉綠體、過(guò)氧化物酶體的定位系數(shù)分別為6.0、4.0、3.0。

2.2.7GrAACT基因稀有密碼子分析使用在線軟件對(duì)GrAACT基因進(jìn)行稀有密碼子分析，結(jié)果表明：GrAACT基因中稀有密碼子占1.48%，且無(wú)二聯(lián)或三聯(lián)稀有密碼子連續(xù)出現(xiàn)的情況，因此可選用大腸桿菌表達(dá)菌BL21或Rosetta（DE3）進(jìn)行原核表達(dá)。

2.3GrAACT基因原核表達(dá)載體的構(gòu)建

使用酶切檢測(cè)重組質(zhì)粒pGEX-4T-1-GrAACT，結(jié)果表明：BamHⅠ、XhoⅠ雙酶切可切出目的基因和載體，且二者長(zhǎng)度之和等于XhoⅠ單酶切產(chǎn)物長(zhǎng)度（圖8），表明GrAACT基因已成功插入載體pGEX-4T-1。

2.4GrAACT基因的原核表達(dá)

將重組質(zhì)粒pGEX-4T-1-GrAACT轉(zhuǎn)化大腸桿菌Rosetta（DE3）后，在37 ℃、1 mmol/L IPTG（終濃度）下，分別

誘導(dǎo)表達(dá)0、2、4、6 h后，提取細(xì)菌總蛋白進(jìn)行SDS-PAGE分析。結(jié)果表明：與對(duì)照（第1～3泳道）相比，pGEX-4T-1-GrAACT轉(zhuǎn)化菌經(jīng)IPTG誘導(dǎo)后，在相對(duì)分子質(zhì)量67 410 u（含GST蛋白26 000 u）左右有1條蛋白條帶，并且隨誘導(dǎo)時(shí)間增加，其蛋白含量逐漸升高（圖9），表明重組質(zhì)粒pGEX-4T-1-GrAACT在大腸桿菌Rosetta（DE3）中成功誘導(dǎo)表達(dá)了GrAACT蛋白。當(dāng)溫度為37 ℃、誘導(dǎo)時(shí)間為6 h時(shí)，蛋白表達(dá)量最大（圖9），細(xì)胞破碎后可直接用于蛋白純化。

2.5GrAACT基因的組織表達(dá)分析

取3年生滇龍膽的根、莖、葉，提取總RNA，反轉(zhuǎn)錄成cDNA后，通過(guò)實(shí)時(shí)熒光定量RT-PCR分析GrAACT基因在根、莖、葉中的表達(dá)情況。圖10結(jié)果表明：GrAACT基因在莖中表達(dá)量最高，分別約為葉、根的13.17、1.36倍，表達(dá)量最低的是葉。

3討論

滇龍膽藥效成分龍膽苦苷是通過(guò)MVA、MEP途徑合成的，在前期筆者已經(jīng)克隆了MEP途徑關(guān)鍵酶1-脫氧-D-木酮糖-5-磷酸還原異構(gòu)酶基因GrDXR[12]。為闡明龍膽苦苷生物合成途徑及其調(diào)控機(jī)理，本研究克隆了MVA途徑中的GrAACT基因。AACT是MVA途徑的第1個(gè)催化酶，因此滇龍膽中GrAACT基因的表達(dá)情況可以直接或間接影響龍膽苦苷的生物合成。研究表明，植物中如丹參[6]、白木香[8]、油茶[9]AACT蛋白均包含硫解酶N端結(jié)構(gòu)域、C端結(jié)構(gòu)域和保守位點(diǎn)。本研究從滇龍膽幼葉中克隆獲得GrAACT基因，其編碼蛋白序列比對(duì)分析結(jié)果表明，GrAACT蛋白屬于AACT家族蛋白，具有硫解酶N端、C端結(jié)構(gòu)域，并且在其C端具有保守位點(diǎn)：NVHGGAVALGHPLGCSG（第350—366位氨基酸）。這些結(jié)果表明，所克隆基因的確為滇龍膽乙酰CoA?；D(zhuǎn)移酶基因。

根據(jù)底物鏈長(zhǎng)特異性，硫解酶被分為2類：第I類是降解型硫解酶（EC 2.3.1.16），定位于線粒體、過(guò)氧化物酶體，在脂肪酸氧化中催化3-酮脂酰CoA斷裂產(chǎn)生乙酰乙酰CoA、縮短的?；鵆oA；第Ⅱ類是合成型硫解酶（EC 2.3.1.9），定位于細(xì)胞質(zhì)，能夠催化2分子乙酰CoA形成1分子乙酰乙酰CoA[13]。生物信息學(xué)分析結(jié)果表明，GrAACT蛋白定位于細(xì)胞質(zhì)，屬于Ⅱ型硫解酶，參與胞質(zhì)中的MVA途徑。在植物中，AACT蛋白的定位還存在爭(zhēng)議。在擬南芥中，參與MVA生物合成途徑的AtAACT2蛋白定位于細(xì)胞質(zhì)[14]；而對(duì)擬南芥過(guò)氧化物酶體蛋白質(zhì)組研究結(jié)果表明，該蛋白也存在于過(guò)氧化物酶體[15]。在向日葵（Helianthus annuus）乙醛酸循環(huán)體中，HaAACT能夠催化朝向乙酰乙酰CoA方向的反應(yīng)[16]，說(shuō)明過(guò)氧化物酶體中也存在AACT。GrAACT蛋白三維模型預(yù)測(cè)結(jié)果表明，該蛋白的活性形式為四聚體。在本研究中，GrAACT單體相對(duì)分子質(zhì)量為41 410 u，其活性形式可能是單體，也可能是四聚體，這需要通過(guò)非變性SDS-PAGE進(jìn)行進(jìn)一步檢測(cè)。

在擬南芥中，參與MVA途徑的AtAACT2基因在根尖、幼葉、莖上部和花藥中大量表達(dá)[5]。組織表達(dá)特異性檢測(cè)結(jié)果表明，GrAACT基因在根中的表達(dá)量較高，因此推測(cè)：滇龍膽根部為龍膽苦苷生物合成的主要器官。滇龍膽用藥部分為根，因?yàn)榕c葉、莖相比，根中龍膽苦苷含量最高[17-18]，這也為上述假說(shuō)提供了證據(jù)。在擬南芥中，參與MVA途徑的AtAACT2基因是多效基因。首先，AtAACT2對(duì)于胚胎形成和正常的雄性配子傳遞是必需的，AtAACT2突變體不能存活；其次，AtAACT2基因的RNAi株系表現(xiàn)出多效的表型，包括減弱的頂端優(yōu)勢(shì)、延長(zhǎng)的生活周期和開(kāi)花時(shí)間、雄性不育、矮化、種子產(chǎn)量降低、根長(zhǎng)變短和植物甾醇類質(zhì)和量的改變等；再次，當(dāng)外源添加甲羥戊酸后，AtAACT2基因的RNAi株系以上的表型和生化改變都可逆轉(zhuǎn)[5]。本研究中的GrAACT蛋白與AtAACT2蛋白相似性高達(dá)83.87%，GrAACT基因是否也為多效基因，需要進(jìn)一步試驗(yàn)驗(yàn)證。

本研究為滇龍膽的分子生物學(xué)研究提供基因資源，并為GrAACT基因功能的解析奠定基礎(chǔ)。下一步將對(duì)GrAACT蛋白進(jìn)行純化和酶活分析，通過(guò)互補(bǔ)試驗(yàn)或過(guò)表達(dá)研究GrAACT基因的功能，從而為龍膽苦苷生物合成途徑的闡明奠定基礎(chǔ)。

參考文獻(xiàn)：

[1]Suyama Y，Kurimoto S I，Kawazoe K，et al. Rigenolide A，a new secoiridoid glucoside with a cyclobutane skeleton，and three new acylated secoiridoid glucosides from Gentiana rigescens Franch[J]. Fitoterapia，2013，91：166-172.

[2]金航，張霽，張金渝，等. 滇龍膽[M]. 昆明：云南科技出版社，2013.

[3]Zhang X，Allan A C，Li C，et al. De novo assembly and characterization of the transcriptome of the Chinese medicinal herb，Gentiana rigescens[J]. International Journal of Molecular Sciences，2015，16（5）：11550-11573.

[4]Hua W，Zheng P，He Y，et al. An insight into the genes involved in secoiridoid biosynthesis in Gentiana macrophylla by RNA-seq[J]. Molecular Biology Reports，2014，41（7）：4817-4825.

[5]Jin H，Song Z，Nikolau B J. Reverse genetic characterization of two paralogous acetoacetyl CoA thiolase genes in Arabidopsis reveals their importance in plant growth and development[J]. The Plant Journal，2012，70（6）：1015-1032.

[6]崔光紅，王學(xué)勇，馮華，等. 丹參乙酰CoA?；D(zhuǎn)移酶基因全長(zhǎng)克隆和SNP分析[J]. 藥學(xué)學(xué)報(bào)，2010，45（6）：785-790.

[7]Bach T J，Boronat A，Campos N，et al. Mevalonate biosynthesis in plants[J]. Critical Reviews in Biochemistry & Molecular Biology，1999，34（2）：107-122.

[8]劉娟，徐艷紅，楊勇，等. 白木香乙酰乙?；o酶A硫解酶基因（AsAACT）的克隆與表達(dá)分析[J]. 中國(guó)中藥雜志，2014，39（6）：972-980.

[9]張琳，譚曉風(fēng)，胡姣，等. 油茶乙酰CoA?；D(zhuǎn)移酶基因 cDNA 克隆及序列特征分析[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)，2011，31（8）：108-112.

[10]Gao W，Sun H X，Xiao H，et al. Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis in Salvia miltiorrhiza[J]. BMC Genomics，2014，15（1）：1-14.

[11]Nakatsuka T，Yamada E，Saito M，et al. Construction of the first genetic linkage map of Japanese gentian （Gentianaceae）[J]. BMC Genomics，2012，13（1）：1-15.

[12]張曉東，趙靜，李彩霞，等. 滇龍膽1-脫氧-D-木酮糖-5-磷酸還原異構(gòu)酶基因GrDXR的克隆、序列分析與原核表達(dá)[J]. 中草藥，2014，45（16）：2378-2384.

[13]Meng Y，Li J. Cloning，expression and characterization of a thiolase gene from Clostridium pasteurianum[J]. Biotechnology Letters，2006，28（16）：1227-1232.

[14]Carrie C，Murcha M W，Millar A H，et al. Nine 3-ketoacyl-CoA thiolases （KATs） and acetoacetyl-CoA thiolases （ACATs） encoded by five genes in Arabidopsis thaliana are targeted either to peroxisomes or cytosol but not to mitochondria[J]. Plant Molecular Biology，2007，63（1）：97-108.

[15]Reumann S，Babujee L，Ma C，et al. Proteome analysis of Arabidopsis leaf peroxisomes reveals novel targeting peptides，metabolic pathways，and defense mechanisms[J]. The Plant Cell，2007，19（10）：3170-3193.

[16]Oeljeklaus S，F(xiàn)ischer K，Gerhardt B. Glyoxysomal acetoacetyl-CoA thiolase and 3-oxoacyl-CoA thiolase from sunflower cotyledons[J]. Planta，2002，214（4）：597-607.

[17]楊美權(quán)，張金渝，沈濤，等. 不同栽培模式對(duì)滇龍膽中龍膽苦苷含量的影響[J]. 江蘇農(nóng)業(yè)科學(xué)，2011（1）：287-289.

[18]王彩云，張曉東，沈濤，等. 龍膽苦苷生物合成途徑研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué)，2014，42（3）：4-10.