劉堰珺，馬　靜，王廣龍，黃　蔚，王　楓，熊愛生

?

胡蘿卜肉桂醇脫氫酶基因的克隆及其對(duì)非生物脅迫的響應(yīng)

劉堰珺，馬靜，王廣龍，黃蔚，王楓，熊愛生*

(南京農(nóng)業(yè)大學(xué)作物遺傳與種質(zhì)創(chuàng)新國(guó)家重點(diǎn)實(shí)驗(yàn)室，園藝學(xué)院，農(nóng)業(yè)部華東地區(qū)園藝作物生物學(xué)與種質(zhì)創(chuàng)制重點(diǎn)實(shí)驗(yàn)室，南京210095)

該實(shí)驗(yàn)采用RT-PCR技術(shù)，從胡蘿卜‘黑田五寸’中克隆獲得了編碼肉桂醇脫氫酶(Cinnamyl alcohol dehydrogenase, CAD)的基因DcCAD。DcCAD序列長(zhǎng)1 074 bp，編碼357個(gè)氨基酸。親水/疏水性分析表明，DcCAD屬于親水性蛋白。系統(tǒng)進(jìn)化分析顯示，DcCAD與番茄CAD(XP_010314515.1)親緣關(guān)系最近。熒光定量PCR結(jié)果顯示，DcCAD基因在胡蘿卜根、葉片和葉柄中的表達(dá)量差異顯著，其相對(duì)表達(dá)量為葉片>根>葉柄。DcCAD基因?qū)Ω邷?38 ℃)、低溫(4 ℃)、干旱(20% PEG)和鹽(0.2 mol·L-1NaCl)脅迫均有響應(yīng)，尤其對(duì)高溫脅迫和低溫脅迫響應(yīng)明顯，而且高溫處理后1 h和低溫處理后2 h表達(dá)量最高。研究推測(cè),DcCAD基因?qū)}卜抗逆性具有一定的作用。

肉桂醇脫氫酶；同源分析；非生物脅迫；表達(dá)分析；胡蘿卜

木質(zhì)素(lignin)主要分布于木質(zhì)化的次生細(xì)胞壁中，它為植物細(xì)胞壁提供了必要的強(qiáng)度、疏水性以及對(duì)外界惡劣環(huán)境的抗性[1-3]。木質(zhì)素的生物合成是一個(gè)由多種酶參與的復(fù)雜調(diào)控網(wǎng)絡(luò)，肉桂醇脫氫酶(cinnamyl alcohol dehydrogenase, CAD)作為其中的關(guān)鍵酶，人們?cè)诤茉缰熬烷_始對(duì)其進(jìn)行研究[4]。國(guó)內(nèi)外已對(duì)植物中肉桂醇脫氫酶的結(jié)構(gòu)性質(zhì)、相應(yīng)的功能以及轉(zhuǎn)基因方面進(jìn)行了一些研究[5]。

由于肉桂醇脫氫酶對(duì)木質(zhì)素單體的生物合成起著重要的作用，目前，已有多個(gè)編碼肉桂醇脫氫酶的基因從不同植物中被成功克隆，如火炬松(Pinustaeda)[6]、九眼獨(dú)活(Araliacordata)[7]、紫花苜蓿(Medicagosativa)[8-9]、多年生黑麥草(Loliumperenne)[10]和高粱(Sorghumbicolor)[11]等。研究也證實(shí)，調(diào)控木質(zhì)素合成途徑中相關(guān)基因的表達(dá)，可以控制木質(zhì)素的生物合成。通過(guò)分析楊樹(Populuslaurifolia)[12]、小麥(Triticumaestivum)[13]、二穗短柄草(Brachypodiumdistachyon)[14]中肉桂醇脫氫酶基因的表達(dá)模式發(fā)現(xiàn)，肉桂醇脫氫酶基因在植物中以多基因家族的形式存在。植物在受到各種生物及非生物脅迫時(shí)，其木質(zhì)素代謝相關(guān)基因的表達(dá)會(huì)發(fā)生變化，其木質(zhì)素含量也會(huì)相應(yīng)改變。

胡蘿卜(DaucuscarotaL.)是傘形科胡蘿卜屬的一種重要蔬菜，在世界各地廣泛栽培。胡蘿卜以肉質(zhì)根為食用部位，含有豐富的營(yíng)養(yǎng)物質(zhì)和水溶性纖維[15]。作為根用型蔬菜，胡蘿卜肉質(zhì)根中木質(zhì)素含量的增加會(huì)影響胡蘿卜的品質(zhì)。目前，有關(guān)傘形科蔬菜作物中木質(zhì)素代謝相關(guān)研究報(bào)道很少。

本研究以胡蘿卜‘黑田五寸’為實(shí)驗(yàn)材料，基于轉(zhuǎn)錄組和基因組數(shù)據(jù)，通過(guò)RT- PCR方法克隆出木質(zhì)素合成關(guān)鍵酶—肉桂醇脫氫酶的基因(DcCAD)，并對(duì)其進(jìn)行了較為詳盡的分析。通過(guò)實(shí)時(shí)定量PCR技術(shù)，研究DcCAD基因在胡蘿卜不同組織及不同逆境處理下的表達(dá)情況，為進(jìn)一步深入研究胡蘿卜中肉桂醇脫氫酶的功能及應(yīng)用奠定基礎(chǔ)。

1　材料和方法

1.1材料

胡蘿卜‘黑田五寸’于2015年10月份種植于南京農(nóng)業(yè)大學(xué)作物遺傳與種質(zhì)創(chuàng)新國(guó)家重點(diǎn)實(shí)驗(yàn)室人工氣候室。胡蘿卜植株種植在蛭石和有機(jī)質(zhì)1∶1的混合基質(zhì)中，光照和溫度設(shè)定為：14 h 25 ℃光照以及10 h 18 ℃黑暗，相對(duì)濕度60%～70%以及240 μmol·m-2·s-1光照強(qiáng)度。

大腸桿菌菌株DH5α由本實(shí)驗(yàn)室保存；質(zhì)粒載體pMD19-T vector simple、ExTaq酶、反轉(zhuǎn)錄酶MluI、DL marker 2000和dNTP熒光定量染料SYBR Green I等均購(gòu)自大連TaKaRa公司。

1.2方法

1.2.1取樣對(duì)‘黑田五寸’2月齡植株分別進(jìn)行高溫(38 ℃)、低溫(4 ℃)、干旱(20% PEG)和鹽(0.2 mol·L-1NaCl)處理，分別取處理 0、1、2、4、8和24 h的成熟健康葉片，同時(shí)取未處理植株葉片為對(duì)照。取2月齡‘黑田五寸’植株，剪取其生長(zhǎng)健康的葉片用于葉片總RNA提取，反轉(zhuǎn)錄成cDNA，用于克隆DcCAD基因。取2月齡‘黑田五寸’植株的葉片、葉柄、根為材料進(jìn)行組織表達(dá)分析。每次取樣3個(gè)重復(fù)，取樣后立即用液氮速凍，并保存于-80 ℃冰箱。

1.2.2RNA提取及cDNA合成用RNA simple Total Kit(Tiangen公司，中國(guó)北京)提取總RNA，RNA濃度和質(zhì)量利用Nanodrop ND-100分光光度計(jì)(NanoDrop公司，美國(guó))檢測(cè)。利用Prime Script RT reagent Ki(TaKaRa公司，中國(guó)大連)將提取的總RNA反轉(zhuǎn)錄成cDNA。所有cDNA樣品稀釋10倍用于PCR。

1.2.3胡蘿卜DcCAD基因序列的獲得與克隆基于本實(shí)驗(yàn)室測(cè)定的胡蘿卜轉(zhuǎn)錄組數(shù)據(jù)和基因組數(shù)據(jù)[16]，拼接得到DcCAD基因的序列，并設(shè)計(jì)正向引物DcCAD-F(5′-ATGGCTCAAACAACTCCCAATCACA-3′)和反向引物DcCAD-R(5′-CTATGGTTTTCCAGCAATATCAAT-3′)。以2月齡‘黑田五寸’新鮮主根cDNA第一鏈為模板進(jìn)行擴(kuò)增，反應(yīng)條件為94 ℃預(yù)變性4 min；94 ℃變性30 s；54 ℃退火30 s；72 ℃延伸1 min；35個(gè)循環(huán)；72 ℃延伸10 min。利用1.2%瓊脂糖凝膠電泳分離擴(kuò)增產(chǎn)物，回收后連接到pMD19-T simple載體上，并轉(zhuǎn)化到大腸桿菌DH5α。隨后提取質(zhì)粒進(jìn)行PCR鑒定，委托南京金斯瑞生物科技有限公司進(jìn)行DNA測(cè)序。

1.2.4序列分析不同植物相關(guān)基因和氨基酸序列均來(lái)自NCBI數(shù)據(jù)庫(kù)；利用BioXM 2.6軟件對(duì)DcCAD基因的可閱讀框長(zhǎng)度進(jìn)行預(yù)測(cè)，并翻譯成對(duì)應(yīng)的氨基酸序列。通過(guò)DNAMAN進(jìn)行氨基酸序列比對(duì)及親水性/疏水性分析，然后用Mega 5構(gòu)建進(jìn)化樹和編輯并生成報(bào)告圖形[17]。蛋白質(zhì)基本性質(zhì)分析使用ExPASy相關(guān)軟件完成[18]。通過(guò)SWISS-MODEL(http://swiss-model.expase.org)預(yù)測(cè)DcCAD蛋白的三維結(jié)構(gòu)。

1.2.5實(shí)時(shí)定量PCR反應(yīng)實(shí)時(shí)定量PCR采用SYBR PremixExTaq試劑盒(大連TaKaRa公司)，在ABI 7300熒光定量PCR儀上完成，采用iQTM5 Real-time PCR System完成，按照操作說(shuō)明進(jìn)行。相對(duì)定量方法參照2-ΔΔCt法。以胡蘿卜Tubulin為參照基因[19]，引物為DcA15-F(5′-GAGTGGAGTTACCTGCTGCCTTC-3′)和DcA15-R(5′-ATGTAGACGAGGGAACGGAATCCAG-3′)。DcCAD基因表達(dá)分析引物為DcCAD-F(5′-CGGATAAGCCTCTGGACCTACCT-3′)和DcCAD-R(5′-TTGACAACTTCGATACCACACCTGAT-3′)，目標(biāo)基因與Tubulin一起擴(kuò)增。實(shí)驗(yàn)進(jìn)行3次生物學(xué)重復(fù)，對(duì)數(shù)值進(jìn)行了差異顯著水平分析。

1.2.6數(shù)據(jù)分析采用 Excel 2003和SPSS 16.0進(jìn)行數(shù)據(jù)整理和基因表達(dá)分析，采用Duncan測(cè)驗(yàn)分析數(shù)據(jù)差異顯著性。

2　結(jié)果與分析

2.1胡蘿卜DcCAD基因克隆

PCR結(jié)果(圖1)顯示，以胡蘿卜‘黑田五寸’的cDNA為模版擴(kuò)增得到1 000 bp左右產(chǎn)物，其長(zhǎng)度與

預(yù)期相一致。測(cè)序結(jié)果顯示，DcCAD基因長(zhǎng)1 074 bp, 編碼357個(gè)氨基酸(圖2)，預(yù)測(cè)其蛋白質(zhì)相對(duì)分子質(zhì)量為39.09 kD，等電點(diǎn)6.7。

2.2氨基酸序列比對(duì)和理化性質(zhì)分析

對(duì)DcCAD的氨基酸序列進(jìn)行BLASTp同源性檢索，結(jié)果(圖)表明：DcCAD的氨基酸序列與擬南芥(Arabidopsisthaliana，NP_195643.1)、黃蒿(Artemisiaannua，ACB54931.1)、茶樹(Camelliasinensis，ACV74414.1)、辣椒(Capsicumannuum，ACF17645.1)、黃瓜(Cucumissativus，XP_004136373.1)、黑麥草(Loliumperenne，AAL99536.1)、蘋果(Malusdomestica，XP_008359757.1)、煙草(Nicotianaattenuata，AFP43763.1)、馬鈴薯(Solanumtuberosum，XP_006356896.1)、小麥(Triticumaestivum，CDM85988.1)、玉米(Zeamays，NP_001136687.1)、水稻(Oryzasativa，NP_001052290.1)、高粱(Sorghumbicolor，XP_002446076.1)、菠菜(Spinaciaoleracea，KNA09004.1)和蓖麻(Ricinuscommunis，XP_011083448.1)等植物的氨基酸序列具有較高的相似度，相似性均在60%以上(圖3)。

圖1　胡蘿卜DcCAD基因克隆圖譜Fig. 1　PCR amplification of DcCAD gene from carrot

圖2　胡蘿卜DcCAD基因及其編碼氨基酸序列Fig. 2　The nucleotide acid and deduced amino acid sequences of DcCAD gene from carrot

星號(hào)表示3個(gè)催化鋅結(jié)合位點(diǎn)；箭頭表示3個(gè)鋅結(jié)合位點(diǎn)； CAD蛋白登錄號(hào)依次為：蓖麻(XP_011083448.1)、菠菜(KNA09004.1)、茶樹(ACV74414.1)、高粱(XP_002446076.1)、黑麥草(AAL99536.1)、胡蘿卜(ALN38479.1)、黃瓜(XP_004136373.1)、黃蒿(ACB54931.1)、辣椒(ACF17645.1)、馬鈴薯(XP_006356896.1)、擬南芥(NP_195643.1)、蘋果(XP_008359757.1)、水稻(NP_001052290.1)、小麥(CDM85988.1)、煙草(AFP43763.1)和玉米(NP_001136687.1)圖3　胡蘿卜DcCAD的氨基酸序列與其他物種類似蛋白質(zhì)氨基酸序列的多重比對(duì)The asterisk marks the 3 catalytic Zn binding sites; The arrow marks the 3 structural Zn binding sites; The accession numbers of CAD from different plants were: Ricinus communis (XP_011083448.1), Spinacia oleracea (KNA09004.1), Camellia sinensis(ACV74414.1), Sorghum bicolor (XP_002446076.1), Lolium perenne (AAL99536.1), Daucus carota(ALN38479.1), Cucumis sativus (XP_004136373.1), Artemisia annua (ACB54931.1), Capsicum annuum (ACF17645.1), Solanum tuberosum (XP_006356896.1), Arabidopsis thaliana (NP_195643.1), Malus domestica (XP_008359757.1), Oryza sativa (NP_001052290.1), Triticum aestivum (CDM85988.1), Nicotiana attenuata (AFP43763.1), and Zea mays (NP_001136687.1)Fig. 3　Alignment of amino acid sequences of carrot DcCAD and other different plants

植物Plant登錄號(hào)Accessionnumber氨基酸殘基數(shù)Numberofaminoacids理論相對(duì)分子質(zhì)量Theoreticalrelativemolecularmass/Da理論等電點(diǎn)TheoreticalpI酸性氨基酸比例Ratioofpositiveaminoacids/%堿性氨基酸比例Ratioofnegativeaminoacids/%脂肪族氨基酸指數(shù)Aliphaticindex/%平均疏水性Averageofhydrophobicity胡蘿卜D.carotaALN38479.135739088.06.730111287.09-0.133小麥T.aestivumCDM85988.135438503.66.384121586.98-0.093擬南芥A.thalianaNP_195643.136043058.56.668121590.51-0.077黃蒿A.annuaACB54931.136147931.55.629121696.23-0.003茶樹C.sinensisACV74414.136052691.57.095131592.29-0.144辣椒C.annuumACF17645.135753961.26.179131694.14-0.143黃瓜C.sativusXP_004136373.136259078.97.497131590.59-0.210黑麥草L.perenneAAL99536.137072474.26.696131697.15-0.122蘋果M.domesticaXP_008359757.136678411.98.122131594.79-0.193煙草N.attenuataAFP43763.135782483.45.413131697.14-0.189馬鈴薯S.tuberosumXP_006356896.136087503.26.567131697.33-0.213玉米Z.maysNP_001136687.136138711.37.269131692.33-0.143水稻O.sativaNP_001052290.136043535.17.230131692.57-0.145高粱S.bicolorXP_002446076.136247805.76.575131693.91-0.143菠菜S.oleraceaKNA09004.135957354.57.286131689.81-0.195蓖麻R.communisXP_011083448.136252929.36.868131694.25-0.199

對(duì)上述植物進(jìn)行氨基酸組成成分及理化性質(zhì)分析，結(jié)果(表1)發(fā)現(xiàn)，各植物氨基酸殘基數(shù)在354～370之間不等；理論相對(duì)分子質(zhì)量的差異較大，理論等電點(diǎn)在5.413～8.122之間；堿性氨基酸所占比例略高于酸性氨基酸，脂肪族氨基酸的比例在86.98%～97.33%之間；蛋白質(zhì)可溶性預(yù)測(cè)中平均疏水性為-0.213～-0.003，其中黃蒿的疏水性最強(qiáng)，馬鈴薯的疏水性最弱。

2.3氨基酸親水性/疏水性分析及結(jié)構(gòu)域分析

對(duì)DcCAD氨基酸序列進(jìn)行疏水性/親水性分析，結(jié)果(圖4)表明，DcCAD疏水區(qū)域中，第187～202位疏水性較強(qiáng)，第165位的丙氨酸(Ala)和第189位、276位的纈氨酸(Val)疏水性最強(qiáng)。親水區(qū)域中182～187位親水性較強(qiáng)，第33位的天冬酰胺(Asn)、112位的谷氨酸(Glu)和223位的賴氨酸(Lys)疏水性最強(qiáng)。結(jié)果表明，DcCAD屬于親水性蛋白。對(duì)胡蘿卜DcCAD編碼的氨基酸序列進(jìn)行Blastp同源性檢索，結(jié)果顯示胡蘿卜DcCAD屬于MDR超家族，歸為CAD1家族，具有多個(gè)Zn結(jié)合位點(diǎn)，屬于多結(jié)構(gòu)域蛋白序列。

圖4　胡蘿卜DcCAD氨基酸序列的親水性和疏水性Fig. 4　Predicted hydrophilicity and hydrophobicty of deduced amino acid sequence of DcCAD from carrot

2.4蛋白三級(jí)結(jié)構(gòu)預(yù)測(cè)和進(jìn)化分析

胡蘿卜DcCAD進(jìn)行蛋白質(zhì)三級(jí)結(jié)構(gòu)同源建模分析的結(jié)果如圖5所示。胡蘿卜DcCAD蛋白有9個(gè)α螺旋和16個(gè)β折疊，與擬南芥的CAD蛋白有較高的相似性。通過(guò)Blastp檢索，將DcCAD的氨基酸序列與GenBank已報(bào)道的其他植物的CAD蛋白序列進(jìn)行比對(duì)，用Mega 5構(gòu)建系統(tǒng)進(jìn)化樹(圖6)。

圖6　胡蘿卜DcCAD與其他植物CAD的系統(tǒng)進(jìn)化樹Fig. 6　Phylogenetic tree of DcCAD from carrot and CAD from other species

圖5　胡蘿卜DcCAD三級(jí)結(jié)構(gòu)Fig. 5　The three-dimension structural model of DcCAD from carrot

圖中小寫字母表示0.05水平差異顯著性圖7　胡蘿卜DcCAD基因在不同組織中的表達(dá)情況The different letters mean significant difference at 0.05 levelFig. 7　Expression analysis of the DcCAD gene in different tissues of carrot

豎條上方不同字母代表同一處理在0.05水平上的差異顯著性圖8　胡蘿卜DcCAD基因在不同逆境以及不同時(shí)間段處理的表達(dá)Different letters on the vertical bars indicate significant difference in the same treatment at 0.05 levelFig. 8　Expression profiles of the DcCAD gene under different abiotic stress treatments in carrot

結(jié)果表明：胡蘿卜與茄科的番茄親緣關(guān)系最近，兩者氨基酸相似性達(dá)78.89%。與菠菜、茶樹、葡萄CAD進(jìn)化關(guān)系較近，與薔薇科的蘋果、胡麻科的芝麻、山茶科的蓖麻以及楊柳科的楊樹和毛果楊CAD進(jìn)化關(guān)系較遠(yuǎn)，與葫蘆科的黃瓜、甜瓜以及桉樹CAD進(jìn)化關(guān)系最遠(yuǎn)。

2.5胡蘿卜DcCAD基因在不同組織及不同非生物脅迫條件下的表達(dá)

以胡蘿卜‘黑田五寸’的根、葉柄和葉片為材料，通過(guò)熒光定量PCR檢測(cè)DcCAD基因在胡蘿卜不同組織中的表達(dá)情況。結(jié)果(圖7)顯示，DcCAD基因在胡蘿卜不同組織中的表達(dá)量不同，且差異較大。葉片中表達(dá)量高于葉柄和根，葉片中表達(dá)量分別是根的10倍、葉柄的17.7倍。

用熒光定量PCR方法測(cè)定胡蘿卜‘黑田五寸’在高溫(38 ℃)、低溫(4 ℃)、干旱(20% PEG)和鹽(0.2 mol·L-1NaCl)4種脅迫處理后的表達(dá)情況。結(jié)果(圖8)顯示，38 ℃處理下，1、2、4、8和24 h表達(dá)量均較高，分別是未處理對(duì)照的65.8、7.28、8.79、39.66和21.47倍。4 ℃處理下隨著處理時(shí)間的延長(zhǎng)，DcCAD基因的表達(dá)量呈現(xiàn)低-高-低的趨勢(shì)，且處理2 h后表達(dá)量最高，是對(duì)照的18.75倍。20% PEG處理下，除2 h和24 h外，其余處理時(shí)間的表達(dá)低于對(duì)照，2 h表達(dá)量達(dá)到峰值，是對(duì)照的3.95倍。鹽處理下，DcCAD基因總體表現(xiàn)為先上升，后下降，再上升，并于24 h達(dá)到峰值。

3　討　論

肉桂醇脫氫酶在植物木質(zhì)素生物合成中發(fā)揮著重要的作用，其活性的高低影響著植物木質(zhì)素的含量和木質(zhì)素單體的組成[20-21]。目前已經(jīng)從多種植物中分離編碼肉桂醇脫氫酶基因，并通過(guò)轉(zhuǎn)基因方法對(duì)其功能進(jìn)行了驗(yàn)證。本研究以‘黑田五寸’胡蘿卜為研究對(duì)象，克隆獲得編碼肉桂醇脫氫酶基因DcCAD。氨基酸理化性質(zhì)分析結(jié)果顯示，DcCAD屬親水性蛋白，但不同物種間的氨基酸理化性質(zhì)存在一定差異，這可能與其參與不同的逆境調(diào)控有關(guān)[22]。通過(guò)對(duì)胡蘿卜DcCAD基因及編碼蛋白的同源性分析發(fā)現(xiàn)，與DcCAD相似性最高的是番茄CAD。DcCAD與葡萄科、山茶科、藜科的CAD相似性較高，與單子葉禾本科、葫蘆科、楊柳科的CAD親緣關(guān)系較遠(yuǎn)。植株中的CAD在不同科之間保守性并不高，DcCAD與其他物種同源基因的相似性只有60 %左右。本研究通過(guò)對(duì)不同物種CAD進(jìn)化關(guān)系進(jìn)行分析，為探討植物中CAD起源和進(jìn)化提供一定的借鑒。

植物對(duì)外界環(huán)境變化的調(diào)控網(wǎng)絡(luò)十分復(fù)雜，當(dāng)處于逆境條件下時(shí)，各組織器官會(huì)行使其特定的生理生化功能，產(chǎn)生相應(yīng)的應(yīng)對(duì)機(jī)制[23-25]。木質(zhì)化過(guò)程也是植物響應(yīng)逆境脅迫的一種方式，它可能被臭氧、重金屬、干旱、紫外線、機(jī)械損傷、病蟲害等逆境所誘導(dǎo)[26-29]。肉桂醇脫氫酶通過(guò)將醛類轉(zhuǎn)化為相應(yīng)的醇類來(lái)參與木質(zhì)素的生物合成，同時(shí)它還與植物的逆境應(yīng)答有關(guān)[26,30-31]。當(dāng)植株受到外界環(huán)境的傷害時(shí)，會(huì)合成木質(zhì)素來(lái)響應(yīng)逆境脅迫。本研究中對(duì)胡蘿卜植株進(jìn)行高溫、低溫、干旱、鹽脅迫處理，發(fā)現(xiàn)DcCAD基因?qū)Ψ巧锩{迫響應(yīng)明顯。此外DcCAD基因在胡蘿卜不同組織中的表達(dá)量不同，在葉片中的表達(dá)量最高，葉柄中的表達(dá)量最低。這些結(jié)果說(shuō)明胡蘿卜DcCAD基因?qū)δ婢趁{迫有響應(yīng)，由此可推測(cè)它對(duì)胡蘿卜的抗逆性起著一定的作用。

[1]PETER G, NEALE D. Molecular basis for the evolution of xylem lignification[J].CurrentOpinioninPlantBiology, 2004, 7(6): 737-742.

[2]GUO D M, RAN J H, WANG X Q. Evolution of the cinnamyl/sinapyl alcohol dehydrogenase (CAD/SAD) gene family: The emergence of real lignin is associated with the origin of bona fide CAD[J].JournalofMolecularEvolution, 2010, 71(3): 202-218.

[3]WENG J K, CHAPPLE C. The origin and evolution of lignin biosynthesis[J].NewPhytologist, 2010, 187(2):273-285.

[4]張魯斌, 谷會(huì), 弓德強(qiáng), 等. 植物肉桂醇脫氫酶及其基因研究進(jìn)展[J].西北植物學(xué)報(bào), 2011, 31(1):204-211.

ZHANG L B, GU H, GONG D Q，etal. Research progress of cinnamyl alcohol dehydrogenase and its gene[J].ActaBotanicaBoreali-OccidentaliaSinica, 2011, 31(1):204-211.

[5]曹佳強(qiáng), 李波, 楊洋, 等. 木質(zhì)素生物合成中肉桂醇脫氫酶基因(CAD)的研究進(jìn)展[J].分子植物育種, 2014,10(5):1 034-1 043.

CAO J Q, LI B, YANG Y,etal. Advanced progress in the cinnamyl alcohol dehydrogenase gene involved in lignin biosynthesis[J].MolecularPlantBreeding, 2014, 10(5):1 034-1 043.

[6]O’MALLEY D M, PORTER S, SEDEROFF R R. Purification, characterization, and cloning of cinnamyl alcohol dehydrogenase in loblolly pine (PinustaedaL.)[J].PlantPhysiology, 1992, 98(4), 1 364-1 371.

[7]HIBINO T, SHIBATE D, CHEN J Q,etal. Cinnamyl alcohol dehydrogenase fromAraliacordata: cloning of the cDNA and expression of the gene in lignified tissues[J].PlantCellPhysiology, 1993, 34(5): 659-665.

[8]BRILL E M, ABRAHAMS S, HAYES C M,etal. Molecular characterization and expression of a wound-inducible cDNA encoding a novel cinnamyl-alcohol dehydrogenase enzyme in lucerne (MedicagosativaL.)[J].PlantMolecularBiology, 1999, 41(2): 279-291.

[9]VAN DOORSSELAERE J, BAUCHER M, FEUILLET C,etal. Isolation of cinnamyl alcohol dehydrogenase cDNAs from two important economic species: alfalfa and poplar. Demonstration of a high homology of the gene within angiosperms[J].PlantPhysiologyandBiochemistry, 1995, 33(1): 105-109.

[10]LYNCH D, LIDGETT A, MCINNES R,etal. Isolation and characterisation of three cinnamyl alcohol dehydrogenase homologue cDNAs from perennial ryegrass (LoliumperenneL.)[J].JournalofPlantPhysiology, 2002, 159(6):653-660.

[11]TSURUTA S, EBINA M, NAKAGAWA H,etal. Isolation and characterization of cDNA encoding cinnamyl alcohol dehydrogenase (CAD) in sorghum (Sorghumbicolor(L.) Moench)[J].JapaneseSocietyofGrasslandScience, 2007, 53(2): 103-109.

[12]EUDES A, GEORGE A, MUKERJEE P,etal. Biosynthesis and incorporation of side-chain-truncated lignin monomers to reduce lignin polymerization and enhance saccharification[J].PlantBiotechnologyJournal, 2010, 10(5): 609-620.

[13]MA Q H. Functional analysis of a cinnamyl alcohol dehydrogenase involved in lignin biosynthesis in wheat[J].JournalofExperimentalBotany, 2012, 61(10): 2 735-2 744.

[14]BUKH C, NORD-LARSEN P H, RASMUSSEN S K. Phylogeny and structure of the cinnamyl alcohol dehydrogenase gene family inBrachypodiumdistachyon[J].JournalofExperimentalBotany, 2012, 63(17): 6 223-6 236.

[15]莊飛云, 歐承剛, 趙志偉, 等. 胡蘿卜育種回顧及展望[J].中國(guó)蔬菜, 2008, 15(3): 41-44.

ZHUANG F Y, OU C G, ZHAO Z W,etal. Review and prospect of carrot breeding[J].ChinaVegetables, 2008, 15(3):41-44.

[16]XU Z S, TAN H W, WANG F,etal. CarrotDB: a genomic and transcriptomic database for carrot[J].DatabasetheJournalofBiologicalDatabases&Curation, 2014, 40 (5):1-8.

[17]TAMURA K, PETERSON D, PETERSON N,etal. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J].MolecularBiologyandEvolution, 2011, 28(10): 2 731-2 739.

[18]GASTEIGER E, GATTIKER A, HOOGLAND C,etal. ExPASy: the proteomics server for in-depth protein knowledge and analysis[J].NucleicAcidsResearch, 2003, 31(13): 3 784-3 788.

[19]TIAN C, JIANG Q, WANG F,etal. Selection of suitable reference genes for qPCR normalization under abiotic stresses and hormone stimuli in carrot leaves[J].PLOSOne, 2015, 10(2):128-136.

[20]PILATE G, GUINEY E, HOLT K. Field and pulping performances of transgenic trees with altered lignification[J].NatureBiotechnology,2002, 20(4):607-612.

[21]CHABANNES M, BARAKATE A, LAPIERRE C,etal. Strong decrease in lignin content without significant alteration of plant development is induced by simulataneous downregulation of cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD) in tobacco plants[J].PlantJournal, 2001, 28(3):257-270.

[22]YAMAGUCHI-SHINOZAKI K, SHINOZAKI K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses[J].AnnualReviewofPlantBiology, 2006, 57(6): 781-803.

[23]HU Y, LI W C, XU Y Q,etal. Differential expression of candidate genes for lignin biosynthesis under drought stress in maize leaves[J].JournalofAppliedGenetics, 2009, 50(3): 213-223.

[24]ZELLER G, HENZ S R, WIDMER C K,etal. Stress-induced changes in theArabidopsisthalianatranscriptome analyzed using whole-genome tiling arrays[J].PlantJournal, 2009, 58(6): 1 068-1 082.

[25]WALLEY J W, DEHESH K. Molecular mechanisms regulating rapid stress signaling networks inArabidopsis[J].JournalofIntegrativePlantBiology, 2010, 52(4): 354-359.

[26]RICHET N, TOZO K, ATIF D,etal. The response to daylight or continuous ozone of phenylpropanoid and lignin biosynthesis pathways in poplar differs between leaves and wood[J].Planta, 2012, 236(2): 727-737.

[27]BETZ G, KNAPPE C, LAPIERRE C,etal. Ozone affects shikimate pathway transcripts and monomeric lignin composition in European beech (FagussylvaticaL.)[J].EuropeanJournalofForestResearch, 2009, 128(2):109-116.

[28]TRONCHET M, BALAGUé C, KROJ T,etal. Cinnamyl alcohol dehydrogenases-C and D, key enzymes in lignin biosynthesis, play an essential role in disease resistance inArabidopsis[J].MolecularPlantPathology, 2010, 11(1):83-92.

[29]ZHAO Q, DIXON R A. Transcriptional networks for lignin biosynthesis: more complex than we thought[J].TrendsinPlantScience, 2011, 16(4): 227-233.

[30]BI C, CHEN F, JACKSON L,etal. Expression of lignin biosynthetic genes in wheat during development and upon infection by fungal pathogens[J].PlantMolecularBiologyReporter, 2011, 29(1): 149-161.

[31]KAUR H, SHAKER K, HEINZEL N,etal. Environmental stresses of field growth allow cinnamyl alcohol dehydrogenase-deficientNicotianaattenuataplants to compensate for their structural deficiencies[J].PlantPhysiology, 2012, 159(4): 1 545-1 570.

(編輯:宋亞珍)

Cloning and Expression Profile Analysis of the Gene Encoding Cinnamyl Alcohol Dehydrogenase under Abiotic Stress in Carrot

LIU Yanjun, MA Jing, WANG Guanglong, HUANG Wei, WANG Feng, XIONG Aisheng*

(State Key Laboratory of Crop Genetics and Germplasm Enhancement, Ministry of Agriculture Key Laboratory of Biology and Germplasm Enhancement of Horticultural Crops in East China, College of Horticulture, Nanjing Agricultural University, Nanjing 210095, China)

Here, the gene encoding DcCAD were cloned from carrot cultivar ‘Kurodagosun’ with RT-PCR method. The length ofDcCADwas 1 074 bp, and encoding 357 amino acids. The DcCAD may belong to hydrophilic protein. Phylogenetic analysis showed that DcCAD closed to CAD from tomato(XP_010314515.1). Quantitative real-time PCR analysis showed that the expression profiles of theDcCADgene are significantly different in different tissues (root, leaf blade, and petiole) of carrot. Relative expression level of theDcCADgene is leaf blade>root>petiole. TheDcCADgene is responsive to various abiotic stresses, such as heat (38 ℃), cold (4 ℃), drought (20% PEG) and salt (0.2 mol·L-1NaCl). The expression level was significantly increased after high temperature or low temperature treatments, respectively. The expression levels peaked at 1 h heat treatment and 2 h cold treatment, respectively. These results suggested thatDcCADgene participate in the response process of carrot under the abiotic stress.

cinnamyl alcohol dehydrogenase; homologous analysis; abiotic stress; expression analysis;Daucuscarota

1000-4025(2016)07-1294-08

10.7606/j.issn.1000-4025.2016.07.1294

2016-04-20;修改稿收到日期:2016-06-13

江蘇省自然科學(xué)基金杰出青年基金(BK20130027)；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)項(xiàng)目(PAPD)

劉堰珺(1993-)，女,碩士研究生,主要從事蔬菜分子生物學(xué)研究。E-mail: 2015104076@njau.edu.cn

熊愛生,教授、博士生導(dǎo)師,主要從事蔬菜遺傳育種與分子生物學(xué)研究。E-mail: xiongaisheng@njau.edu.cn

Q785; Q786; S631.2

A