茄科植物中HCT基因家族的鑒定及進(jìn)化和表達(dá)分析

李倩　蓋江濤　白蓓蓓　葉秀旭　王鵬

摘要：莽草酸/奎寧酸羥基肉桂酰轉(zhuǎn)移酶（shikimic acid/quinic acid hydroxy cinnamyl transferase，簡(jiǎn)稱HCT）是辣椒素苯丙烷代謝途徑中關(guān)鍵的限速酶，其催化的產(chǎn)物對(duì)茄科植物具有重要的生理學(xué)意義。為了探清HCT在茄科植物基因組中的狀況，重點(diǎn)考察了其在茄科植物基因組中的分布，并利用生物信息學(xué)分析的方法，對(duì)HCT基因進(jìn)行全面的鑒定，共獲得擬南芥、煙草、番茄和辣椒的13條基因序列。其次，構(gòu)建了HCT在茄科植物中的系統(tǒng)發(fā)育樹(shù)，并對(duì)其進(jìn)行相關(guān)的結(jié)構(gòu)、功能和表達(dá)的分析。結(jié)果揭示了茄科植物中HCT基因的親緣關(guān)系;HCT基因編碼的蛋白結(jié)構(gòu)域高度保守，該保守性體現(xiàn)在植物的結(jié)構(gòu)和功能上。本研究為進(jìn)一步研究茄科植物中的HCT基因的酶學(xué)及生理功能提供了基礎(chǔ)。

關(guān)鍵詞：茄科植物;莽草酸/奎寧酸羥基肉桂酰轉(zhuǎn)移酶;HCT基因;辣椒素;進(jìn)化;表達(dá)

中圖分類號(hào)： S641.301文獻(xiàn)標(biāo)志碼： A

文章編號(hào)：1002-1302（2019）19-0065-03

收稿日期：2018-07-03

基金項(xiàng)目：中國(guó)熱帶農(nóng)業(yè)科學(xué)院基本科研業(yè)務(wù)專項(xiàng) （編號(hào)：1630032018023）。

作者簡(jiǎn)介：李 倩（1994—），女，安徽池州人，碩士研究生，研究方向?yàn)橹参锘蚪M學(xué)。E-mail：a827680831@163.com。

通信作者：王 鵬，博士，副研究員，主要從事植物天然產(chǎn)物生物技術(shù)、植物基因組學(xué)。E-mail：pwang521@163.com。

辣椒為茄科（Solanaceae）辣椒屬（Capsicum）一年或有限多年生草本植物，是人們常食用的一種風(fēng)味蔬菜。研究發(fā)現(xiàn)，辣椒的果實(shí)中含有豐富的維生素、礦物質(zhì)以及一些微量元素，可以維持人體正常的生理活動(dòng)，提高身體免疫力，增強(qiáng)人體機(jī)能。辣椒素（capsaicin）類物質(zhì)是辣椒中辣味的物質(zhì)基礎(chǔ)，早在1876年，辣椒素就由Thresh分離出來(lái)并命名，Nelson最先解析了辣椒素的結(jié)構(gòu)，表明它是一種含香草酰胺的生物堿[1]，主要的物質(zhì)是辣椒堿和二氫辣椒堿[2]。辣椒素的生物合成是個(gè)復(fù)雜的途徑，其間會(huì)涉及到多種酶和代謝物的參與，然而本質(zhì)上，辣椒素主要由苯丙氨酸途徑合成的香草基胺和由支鏈脂肪酸途徑合成的8-甲基-6-癸烯酰經(jīng)辣椒素合成酶（CS）催化縮合而成[3-4]。研究者在辣椒素的合成途徑中，發(fā)掘到了PAL、C4H、HCT、C3H、COMT、KAS、pAMT、ACL和FAT等多個(gè)參與辣椒素生物合成的基因[3]，查閱相關(guān)文獻(xiàn)發(fā)現(xiàn)，其中莽草酸/奎寧酸羥基肉桂酰轉(zhuǎn)移酶（HCT）是在木質(zhì)素生物合成路徑中發(fā)現(xiàn)較晚的一種酶。

目前普遍認(rèn)同的木質(zhì)素生物合成途徑大致可以分為3個(gè)步驟：（1）莽草酸途徑，由植物光合作用同化產(chǎn)物到芳香氨基酸——苯丙氨酸、絡(luò)氨酸和色氨酸的合成過(guò)程;（2）苯丙烷途徑，苯丙氨酸到羥基肉桂酸及其輔酶A酯類;（3）木質(zhì)素合成的特異途徑，羥基肉桂酸輔酶酯類到合成木質(zhì)素單體及其聚合物的過(guò)程 [5-7] 。其中研究的重點(diǎn)大都是放在苯丙烷途徑和木質(zhì)素合成特異途徑這2個(gè)階段。

HCT是控制木質(zhì)素中的H-單體和G/S-單體相互轉(zhuǎn)化和聚合的關(guān)鍵酶，在苯丙烷3-羥基化步驟上游和下游的苯丙氨酸途徑中發(fā)揮著關(guān)鍵的雙重調(diào)節(jié)作用[8-9]，其中H-單體主要呈現(xiàn)在禾本科植物中，而S-單體主要在被子植物中存在[6]。Hoffmann等從煙梗提取物中分離純化出一種可以水解羥基肉桂酰基-CoA酯的蛋白質(zhì)，從它的序列分析結(jié)果揭示出該煙草基因?qū)儆谝环N植物酰基轉(zhuǎn)移酶的基因家族，其成員也具有各種功能[10]。最佳的酰基供體和酰基轉(zhuǎn)移酶是香豆酰輔酶A和咖啡酰輔酶A，該酶能夠催化莽草酸和奎酸酯的合成，也催化其逆反應(yīng)，從而對(duì)苯丙烷途徑中的木質(zhì)素合成有了明確認(rèn)識(shí)。

HCT在茄科植物辣椒的辣椒素生物合成通路中可催化對(duì)-香豆酰輔酶A，使其形成C3H底物的同時(shí)將C3H的產(chǎn)物咖啡酰莽草酸/奎寧酸進(jìn)一步催化生成咖啡酰輔酶A。在木質(zhì)素單體的合成通路中，不同的木質(zhì)素單體形成的分水嶺是C3H/HCT，其決定了植物體內(nèi)的碳源流向，是苯丙烷途徑中的限速酶，也是控制木質(zhì)素H-單體與G/S-單體的一個(gè)關(guān)鍵點(diǎn)，木質(zhì)素不同單體間的比例差異對(duì)它的含量也會(huì)有所影響[8]。木質(zhì)素自身具有疏水的化學(xué)性質(zhì)，在植物的生長(zhǎng)發(fā)育過(guò)程中，既可以保持植物細(xì)胞內(nèi)的水分，促進(jìn)植物抵御干旱等不良生態(tài)環(huán)境，還可以提高植物抗病、抗擊外來(lái)干擾，維持正常生長(zhǎng)的能力。因此推論，對(duì)于植物中HCT基因的研究具有重要的生物學(xué)意義和功能。目前，已經(jīng)從煙草[10]、擬南芥[11]、番茄[12]、桉樹(shù)[13]、歐美楊[14]中克隆出來(lái)HCT基因，然而，關(guān)于HCT基因在其他茄科植物中的克隆和鑒定，國(guó)內(nèi)鮮有相關(guān)的系統(tǒng)生物信息學(xué)研究報(bào)道。

本研究以擬南芥、煙草和番茄中的HCT基因序列作為參考，以茄科植物中的辣椒作為研究對(duì)象，通過(guò)BLAST序列比對(duì)的方法獲得這些植物HCT基因家族的全部序列，利用生物信息學(xué)分析方法中的系統(tǒng)發(fā)育分析、聚類分析等方法，來(lái)預(yù)測(cè)和分析茄科植物中的HCT基因的相關(guān)性質(zhì)和功能，為后續(xù)HCT基因的研究和利用提供理論基礎(chǔ)。

1 材料與方法

1.1 保守結(jié)構(gòu)域鑒定

把HCT基因的蛋白序列作為query訪問(wèn)，輸入PfamA數(shù)據(jù)庫(kù)開(kāi)始檢索，獲取相應(yīng)的結(jié)構(gòu)域信息[15]。

1.2 數(shù)據(jù)庫(kù)搜索

搜索擬南芥、煙草、番茄、辣椒的HCT基因組，各基因組數(shù)據(jù)從Phytozome中下載得到[16]。通過(guò)搜索該蛋白質(zhì)的模型數(shù)據(jù)庫(kù)，再?gòu)慕Y(jié)構(gòu)域進(jìn)行選擇，獲取HCT家族的序列信息。

1.3 序列比對(duì)與進(jìn)化樹(shù)構(gòu)建

使用probcons比對(duì)之前已經(jīng)搜索到的蛋白質(zhì)序列，再利用Phy ML構(gòu)建與之相對(duì)應(yīng)的系統(tǒng)發(fā)育樹(shù)，之后進(jìn)行1 000次bootstrap分析[17-18]，最后采用Figtree來(lái)顯示構(gòu)建的進(jìn)化樹(shù)結(jié)果。

1.4 相對(duì)表達(dá)量測(cè)定

根據(jù)番茄HCT基因的CDs序列設(shè)計(jì)引物，提取番茄樣品的總RNA，反轉(zhuǎn)成cDNA，篩選引物，之后做熒光定量測(cè)定相對(duì)表達(dá)量。

2 結(jié)果與分析

2.1 結(jié)構(gòu)域鑒定

通過(guò)與Pfam A數(shù)據(jù)庫(kù)匹配，從已知擬南芥的HCT基因的蛋白質(zhì)序列中，鑒定得到了1個(gè)保守結(jié)構(gòu)域，編號(hào)為PF02458.11，顯示該蛋白結(jié)構(gòu)域高度保守，基本占據(jù)了整個(gè)蛋白質(zhì)序列的98%，存在很少的非保守序列區(qū)域（圖1）。

2.2 含PF02458.11結(jié)構(gòu)域的基因鑒定

為全面鑒定含有編號(hào)為PF02458.11的結(jié)構(gòu)域的基因序列，選取了擬南芥（Arabidopsis thaliana）和茄科植物中煙草（Nicotiana tabacum）、番茄（Solanum lycopersicum）、辣椒（Capsicum annuum）為目標(biāo)物種，通過(guò)BLASTP搜索目標(biāo)物種基因組的蛋白質(zhì)模型文件，在BLAST結(jié)果中篩選出PF02458.11結(jié)構(gòu)域的序列，最終得到13個(gè)HCT家族成員（表1）。

2.3 進(jìn)化樹(shù)構(gòu)建

為闡述植物中HCT家族各成員的進(jìn)化關(guān)系，通過(guò)最大似然法構(gòu)建了HCT家族進(jìn)化樹(shù)，結(jié)果（圖2）顯示，此家族可分為2條明顯的獨(dú)立分支（分別命名為Ⅰ和Ⅱ），共有13條氨基酸序列，其中擬南芥1條、番茄2條、煙草4條、辣椒6條。煙草的2個(gè)HCT基因、番茄的1個(gè)HCT基因、辣椒的2個(gè)HCT基因聚為一支，說(shuō)明這些基因的親緣關(guān)系較近;煙草的2個(gè)HCT基因、番茄的1個(gè)HCT基因、辣椒的4個(gè)HCT基因、擬南芥的1個(gè)HCT基因聚為另一支，其中擬南芥單獨(dú)分支，說(shuō)明擬南芥與這些基因的親緣關(guān)系較遠(yuǎn)（圖3）。

2.4 HCT基因相對(duì)表達(dá)量分析

通過(guò)建立熱圖來(lái)進(jìn)行茄科植物中番茄的HCT基因相對(duì)表達(dá)量的差異分析。測(cè)定結(jié)果（圖4）表明，共有野生型和突變型2個(gè)基因型，包含3個(gè)發(fā)育時(shí)期：幼果期、轉(zhuǎn)色期和成熟期，整體來(lái)看它們?cè)诨蛑械谋磉_(dá)量較符合生長(zhǎng)發(fā)育的規(guī)律。其中野生型內(nèi)的Solyc07g043670.1.1、Solyc12g006330.1.1、Solyc09g014280.1.1基因和突變型內(nèi)的Solyc09g092270.2.1基因在幼果期表達(dá)量低，在成熟期表達(dá)量高。

3 討論與結(jié)論

木質(zhì)素的生物合成是植物體內(nèi)碳素次級(jí)代謝的一個(gè)環(huán)節(jié)，HCT是木質(zhì)素合成途徑中研究較晚的一種酶，同時(shí)在辣椒素苯丙氨酸合成通路的苯丙烷C3羥基化的代謝途徑中充當(dāng)關(guān)鍵酶和限速酶的重要角色，其代謝產(chǎn)物在植物生長(zhǎng)發(fā)育的生命活動(dòng)過(guò)程中肩負(fù)著不可或缺的作用，這些物質(zhì)的相對(duì)含量也總是與HCT基因的活性保持著關(guān)聯(lián)性，因而HCT基因?qū)τ谥参锏纳{(diào)控具有重大的生物學(xué)意義。

本研究對(duì)擬南芥、煙草、番茄和辣椒這4種植物的13條HCT基因家族蛋白構(gòu)建系統(tǒng)發(fā)育分析，結(jié)果顯示了茄科植物之間HCT基因的親緣進(jìn)化關(guān)系，還鑒定到另一個(gè)分支的7條序列與擬南芥的HCT基因分別在2個(gè)分支，親緣關(guān)系較遠(yuǎn)，推測(cè)這些基因可能會(huì)在茄科植物中發(fā)揮重要的生理功能，同時(shí)隨著研究的深入，不斷增加的HCT基因完整序列將會(huì)構(gòu)建出更加豐富的系統(tǒng)發(fā)育樹(shù)。相對(duì)表達(dá)量的測(cè)定結(jié)果顯示，本試驗(yàn)設(shè)定的參考對(duì)象番茄符合植物生長(zhǎng)發(fā)育的規(guī)律，在幼果期表達(dá)量低，成熟期表達(dá)量高。

目前，人們對(duì)于苯丙烷代謝途徑中的基因研究較多，雖有一些關(guān)于植物中HCT基因的研究，但鮮有研究從茄科植物的HCT基因家族層面來(lái)分析它的系統(tǒng)進(jìn)化，包括在辣椒素的相關(guān)合成通路中的生長(zhǎng)調(diào)控影響，還有該基因的表達(dá)與辣椒辣味之間的聯(lián)系，以及比較茄科植物中HCT基因的性質(zhì)和功能等等一系列的問(wèn)題亟待解決。因此本研究以茄科植物為研究對(duì)象，對(duì)它們的序列結(jié)構(gòu)、性質(zhì)、進(jìn)化進(jìn)行系統(tǒng)分析，為后續(xù)研究茄科植物中的HCT基因提供參考依據(jù)和理論基礎(chǔ)，促進(jìn)茄科植物在育種工作中的進(jìn)一步發(fā)展。

參考文獻(xiàn)：

[1]Nelson E K. The constitution of capsaicin—The pungent principle of Capsicum[J]. Journal of the American Chemical Society，1919，188（2）：269.

[2]Sarpras M，Gaur R，Sharma V，et al. Comparative analysis of fruit metabolites and pungency candidate genes expression between Bhut Jolokia and other Capsicum species[J]. PLoS One，2016，11（12）：e0167791.

[3]成善漢，吳艷閣，賀申魁. 辣椒素生物合成路徑及調(diào)控基因研究的最新進(jìn)展[J]. 分子植物育種，2008，6（2）：335-340.

[4]Bennett DJ，Kirby[KG*2]G[KG*2]W. Constitution and biosynthesis of capsaicin[J]. Journal of the Chemical Society C Organic，1968，4：442-446.

[5]王雪霞. HCT基因的克隆及其鑒定[D]. 大連：大連工業(yè)大學(xué)，2009.

[6]張士花. 不同4CL基因家族成員對(duì)植物苯丙烷衍生物代謝調(diào)控的研究[D]. 濟(jì)南：山東師范大學(xué)，2011.

[7]Boerjan W，Ralph J，Baucher M. Lignin biosynthesis[J]. Annual Review of Plant Biology，2003，54（1）：519-546.

[8]王雪霞，薛永常，趙文超. 木質(zhì)素生物合成中C3H/HCT的研究進(jìn)展[J]. 生命的化學(xué)，2008，28（5）：650-653.

[9]崔 雪，于月華，申麗婕，等. 海島棉GbHCT基因克隆及生物信息學(xué)分析[J]. 生物技術(shù)通報(bào)，2016，32（4）：80-86.

[10]Hoffmann L，Maury S，Martz F，et al. Purification，cloning，and properties of an acyltransferase controlling shikimate and quinate ester intermediates in phenylpropanoid metabolism[J]. The Journal of Biological Chemistry，2003，278（1）：95-103.

[11]Niggeweg R，Michael A J，Martin C. Engineering plants with increased levels of the antioxidant chlorogenic acid[J]. Nature Biotechnology，2004，22（6）：746-754.

[12]Besseau S，Hoffmann L，Geoffroy P A，et al. Flavonoid accumulation in Arabidopsis repressed in lignin synthesis affects auxin transport and plant growth[J]. Plant Cell，2007，19（1）：148-162.

[13]Harakava R. Genes encoding enzymes of the lignin biosynthesis pathway in Eucalyptus[J]. Genetics & Molecular Biology，2005，28（3）：601-607.

[14]王雪霞，曹 方，薛永常. 楊樹(shù)木質(zhì)素合成酶hct基因的克隆及核苷酸序列分析[J]. 遼寧林業(yè)科技，2009（6）：4-7.

[15]Finn R D，Coggill P，Eberhardt R Y，et al. The Pfam protein families database：towards a more sustainable future[J]. Nucleic Acids Research，2016，44（1）：D279-D285.

[16]Goodstein D M，Shu S，Howson R，et al. Phytozome：a comparative platform for green plant genomics[J]. Nucleic Acids Research，2012，40：1178-1186.

[17]Do C B，Mahabhashyam M S，Brudno M，et al. ProbCons：probabilistic consistency-based multiple sequence alignment[J]. Genome Research，2005，15（2）：330-340.

[18]Guindon S，Delsuc F，Dufayard J F，et al. Estimating maximum likelihood phylogenies with PhyML[J]. Methods in Molecular Biology，2009，537：113-137.