許明　楊志堅(jiān)　黃學(xué)敏　鄭金貴

摘要：【目的】分析藤茶高通量轉(zhuǎn)錄組序列，從中挖掘出黃酮類化合物合成相關(guān)基因，為進(jìn)一步揭示藤茶黃酮類化合物生物合成調(diào)控機(jī)制提供理論參考。【方法】分別采集藤茶的幼葉和成熟葉，提取其總RNA構(gòu)建cDNA文庫(kù)，采用Illumina HiSeqTM 4000高通量測(cè)序平臺(tái)對(duì)藤茶葉片進(jìn)行轉(zhuǎn)錄組測(cè)序，經(jīng)過(guò)濾處理后運(yùn)用Trinity組裝，將獲得的Unigene與Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG 7個(gè)數(shù)據(jù)庫(kù)進(jìn)行比對(duì)注釋，并預(yù)測(cè)Unigenes的編碼區(qū)序列（CDS）;基于KEGG信號(hào)通路富集分析，發(fā)掘藤茶黃酮類化合物合成相關(guān)基因?！窘Y(jié)果】藤茶葉片轉(zhuǎn)錄組測(cè)序獲得82126236條原始測(cè)序序列（Raw reads），過(guò)濾處理后得到80156972條高質(zhì)量序列（Clean reads），進(jìn)一步組裝拼接得到92472條Unigenes，平均長(zhǎng)度為1208 bp，N50長(zhǎng)度為1780 bp，其中，至少在1個(gè)數(shù)據(jù)庫(kù)注釋的Unigenes有84217條，占Unigenes總數(shù)的91.07%，有8944條Unigenes在7個(gè)數(shù)據(jù)庫(kù)均被注釋，占Unigenes總數(shù)的9.67%。在GO數(shù)據(jù)庫(kù)成功注釋的41116條Unigenes可分為生物學(xué)過(guò)程、細(xì)胞組分和分子功能三大類，共56個(gè)小類;在KOG數(shù)據(jù)庫(kù)注釋的14553條Unigenes可分成25類，其中，一般功能預(yù)測(cè)注釋成功的Unigenes最多（1946條）;其次是翻譯后修飾、蛋白質(zhì)翻轉(zhuǎn)、分子伴侶（1776條），參與次生代謝物質(zhì)的生物合成、轉(zhuǎn)運(yùn)和降解的Unigenes較少，僅有319條;KEGG信號(hào)通路富集分析發(fā)現(xiàn)，共有15262條Unigenes注釋到128條KEGG信號(hào)通路，以注釋為代謝的Unigenes最多，為8694條，其中篩選獲得有98個(gè)黃酮類化合物合成相關(guān)基因，分別編碼苯丙烷代謝通路的3種關(guān)鍵酶和類黃酮代謝通路的14種關(guān)鍵酶。藤茶葉片轉(zhuǎn)錄組Unigenes與Swiss-Prot和Nr數(shù)據(jù)庫(kù)比對(duì)，獲得52582條CDS序列，ESTScan 3.0.3預(yù)測(cè)獲得35535條CDS序列。【結(jié)論】藤茶在細(xì)胞過(guò)程、代謝過(guò)程、單有機(jī)體過(guò)程、細(xì)胞和細(xì)胞部分、結(jié)合和催化活性能力分布的基因較豐富，在一般功能、翻譯、翻譯后修飾、蛋白質(zhì)翻轉(zhuǎn)及分子伴侶的基因表達(dá)量較高，具有較強(qiáng)的碳水化合物代謝能力。多種關(guān)鍵酶基因參與藤茶黃酮類化合物的生物合成，推測(cè)其生物合成途徑存在多條分支，調(diào)控機(jī)制也較復(fù)雜。

關(guān)鍵詞： 藤茶;轉(zhuǎn)錄組;黃酮類化合物;基因挖掘;高通量測(cè)序

中圖分類號(hào)： S571.1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2020）08-1797-09

Transcriptome analysis of Ampelopsis grossedentata （Hand. Mazz.） W. T. Wang and mining of putative genes involved in flavonoid biosynthesis

XU Ming1，2， YANG Zhi-jian1，2， HUANG Xue-min3， ZHENG Jin-gui1，2*

（1Key Laboratory of Crop Biotechnology， Fujian Province University， Fuzhou? 350002， China; 2College of Agriculture，F(xiàn)ujian Agriculture and Forestry University， Fuzhou? 350002， China; 3Forestry Bureau of Yongchun County，

Fujian Province， Yongchun， Fujian? 362600）

Abstract：【Objective】The study aimed to obtain the high throughput transcriptome sequence? information and explore the genes related to the flavonoid biosynthesis in Ampelopsis grossedentata（Hand. Mazz.） W. T. Wang. It provided the theoretical basis for further revealing the biosynthetic pathway and regulation mechanism of flavonoids in A. grosse-dentata. 【Method】The young and mature leaves of A. grossedentata were collected respectively for the extraction of total RNA which was used to construct cDNA library. The transcriptome of A. grossedentata leaves was sequenced using an Illumina HiSeqTM 4000 high throughput platform. The sequences were assembled with Trinity after data filtering. The obtained Unigenes were compared and annotated with seven databases， including Nr， Nt， Pfam， Swiss-Prot， GO，KO and KOG， and the coding region sequences（CDS） were predicted by Blast and ESTScan software. KEGG signal pathway enrichment analysis was used to identify genes related to flavonoid synthesis in A. grossedentata. 【Result】A total of 82126236 raw reads were obtained with the transcriptome sequencing. After reads filtering， a total of 80156972 clean reads were obtained， which were used to assemble and merge into 92472 Unigenes with an average length of 1208 bp and a N50 length of 1780 bp.Of these unigenes，? 84217（91.07%） were successfully annotated in at least one of the Nr， Nt， Swiss-Prot， KEGG， COG and GO databases， and 8944 unigenes（9.67%） were annotated in all databases. Among them，41116 Unigenes annotated in GO database were matched to 56 sub-classes of biological function， cell component and molecular function，and 14553 unigenes annotated in KOG database were divided into 25 sub-classes. The number of the unigenes annotated in the general function was the greatest（1946）; the post translational modification， protein turn over and chaperones was the second（1776）; the unigenes involved in the biosynthesis of secondary metabolites， transportation and degradation were less（only 319）. The KEGG signal pathway enrichment analysis results showed that a total of 15262 unigenes were participated in 128 KEGG pathways. The unigenes annotated with metabolism were the greatest in number（8694）， 98 of which participated in the flavonoid synthesis pathway of A. grossedentata. They encoded three key enzymes of phenylpropane metabolic pathway and 14 key enzymes of flavonoid metabolic pathway. By comparing with Swiss-prot and Nr database， 52582 CDS were obtained. A total of 35535 CDS sequences were obtained via ESTScan 3.0.3 prediction. 【Conclusion】A. grossedentata has abundant genes related to cell process， metabolic process， single organism process， cell and cell part， binding and catalytic activities. It shows a high level of gene expression in general function， translation， post translational modification， protein turn over， chaperones.And it has a strong carbohydrate metabolic ability as well. Many key enzymes are involved in the biosynthesis of flavonoids in A. grossedentata， suggesting that there might be multiple branches in the biosynthetic pathway， and the regulation mechanisms might also be complex.

Key words： Ampelopsis grossedentata（Hand. Mazz.） W. T. Wang; transcriptome; flavonoids; gene mining; high-throughput sequencing

Foundation item： National Science and Technology Support Program（2013BAD01B05）;Science and Technology Innovation Project of Fujian Agriculture and Forestry University（KFA17424A）

0 引言

【研究意義】藤茶[Ampelopsis grossedentata（Hand.Mazz.） W. T. Wang]是我國(guó)栽培歷史悠久的一種茶藥兩用植物，其主要藥效成分為黃酮類化合物，具有良好的開發(fā)利用前景（何桂霞等，2007;冉京燕等，2016）。但由于藤茶基因組及轉(zhuǎn)錄組數(shù)據(jù)缺乏，遺傳背景尚不清楚，嚴(yán)重制約了藤茶種質(zhì)資源鑒定和保護(hù)、遺傳多樣性分析及分子育種等工作的開展。此外，產(chǎn)地、栽培條件等諸多因素均會(huì)導(dǎo)致藤茶黃酮含量存在明顯差異，從而影響藤茶產(chǎn)品的質(zhì)量和療效（王家勝等，2014;冉京燕等，2016）。因此，開展藤茶轉(zhuǎn)錄組及其黃酮類化合物的代謝途徑研究，了解藤茶遺傳信息及其藥效成分合成的分子調(diào)控機(jī)制，對(duì)藤茶種質(zhì)資源的持續(xù)高效利用及藥用品質(zhì)提升具有重要意義?！厩叭搜芯窟M(jìn)展】迄今，國(guó)內(nèi)外針對(duì)藤茶的研究主要集中在化學(xué)成分的分離鑒定及藥理作用等方面（Zhang et al.，2006;譚沙等，2015），其基因水平的相關(guān)研究?jī)H有少量報(bào)道。黃海波（2008）從生理生化層面分析了藤茶不同季節(jié)苯丙氨酸解氨酶（PAL）活性與藤茶黃酮含量變化的相關(guān)性;付明等（2013）利用RT-PCR克隆藤茶查爾酮合成酶（CHS）基因;許明等（2017）篩選了適用于藤茶不同組織實(shí)時(shí)熒光定量PCR分析的內(nèi)參基因。轉(zhuǎn)錄組測(cè)序是指利用第二代高通量測(cè)序技術(shù)進(jìn)行cDNA測(cè)序，能在無(wú)參考基因組信息的情況下，對(duì)某一物種的轉(zhuǎn)錄組進(jìn)行全面分析，現(xiàn)已成為揭示特定生物學(xué)過(guò)程分子機(jī)制的一種有效手段（Wilhelm and Landry，2009;嚴(yán)志祥等，2019）。目前，利用高通量測(cè)序技術(shù)已完成較多藥用植物的轉(zhuǎn)錄組分析，并鑒定出大量與黃酮類化合物合成相關(guān)的基因。馬婧等（2016）對(duì)草麻黃進(jìn)行轉(zhuǎn)錄組測(cè)序分析，鑒定出黃酮類化合物代謝相關(guān)基因70個(gè);劉偉等（2018）對(duì)銀杏不同生長(zhǎng)時(shí)期葉片進(jìn)行轉(zhuǎn)錄組測(cè)序分析，篩選出50個(gè)銀杏類黃酮生物合成相關(guān)的候選基因;鄒福賢等（2019）對(duì)3個(gè)生長(zhǎng)時(shí)期的金線蓮進(jìn)行轉(zhuǎn)錄組測(cè)序，找出22個(gè)黃酮類代謝通路上的差異表達(dá)基因，并結(jié)合其黃酮類成分，繪制出相應(yīng)的生物合成途徑;Yang等（2019）對(duì)不同季節(jié)的大葉蛇葡萄葉片進(jìn)行轉(zhuǎn)錄組分析，鑒定出31個(gè)類黃酮生物合成相關(guān)基因和50個(gè)黃酮轉(zhuǎn)運(yùn)和生物修飾相關(guān)基因?！颈狙芯壳腥朦c(diǎn)】至今，鮮見(jiàn)有關(guān)藤茶轉(zhuǎn)錄組測(cè)序及其黃酮類化合物合成相關(guān)基因挖掘的研究報(bào)道?！緮M解決的關(guān)鍵問(wèn)題】利用Illumina HiseqTM 4000高通量測(cè)序平臺(tái)對(duì)藤茶葉片進(jìn)行轉(zhuǎn)錄組測(cè)序，將獲得的序列信息進(jìn)行拼接和組裝，并利用生物大分子公共數(shù)據(jù)庫(kù)對(duì)所得基因進(jìn)行功能注釋和分類，挖掘出藤茶黃酮類化合物生物合成相關(guān)基因，為進(jìn)一步揭示黃酮類化合物生物合成調(diào)控機(jī)制提供理論參考。

1 材料與方法

1. 1 試驗(yàn)材料

供試藤茶引自福建尤溪縣下村林場(chǎng)。Qubit RNA Assay Kit購(gòu)自美國(guó)Invitrogen公司。主要儀器設(shè)備：NanoDrop ND2000超微量分光光度計(jì)（Thermo，美國(guó)）、Agilent 2100生物分析儀（Agilent，美國(guó)）。

1. 2 總RNA提取

分別采集藤茶的幼葉和成熟葉，參照許明等（2015）的方法提取其總RNA，并采用1.0%瓊脂糖凝膠電泳進(jìn)行檢測(cè)。利用Qubit RNA Assay Kit測(cè)定RNA濃度。

1. 3 cDNA文庫(kù)構(gòu)建及測(cè)序

將上述2種葉片RNA等比例混合成1 mg樣品，以干冰保存送至北京諾禾致源生物信息科技有限公司進(jìn)行cDNA文庫(kù)構(gòu)建，并采用Illumina HiSeqTM 4000高通量測(cè)序平臺(tái)進(jìn)行轉(zhuǎn)錄組測(cè)序。

1. 4 轉(zhuǎn)錄組組裝

對(duì)測(cè)序獲得的原始測(cè)序序列（Raw reads）進(jìn)行過(guò)濾處理，得到高質(zhì)量序列（Clean reads），然后采用Trinity對(duì)其進(jìn)行拼接，將得到的轉(zhuǎn)錄本（Transcript）作為參考序列，并取每條基因中最長(zhǎng)的轉(zhuǎn)錄本作為Unigene，用于后續(xù)轉(zhuǎn)錄組序列分析。

1. 5 基因功能注釋和分類

為獲得全面的基因功能信息，將拼接得到的Unigene分別在Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG等7個(gè)數(shù)據(jù)庫(kù)中進(jìn)行功能注釋和分類，以及進(jìn)行KEGG信息通路富集分析。

1. 6 編碼區(qū)序列（CDS）預(yù)測(cè)及分析

按照Swiss-Prot和Nr數(shù)據(jù)庫(kù)的優(yōu)先級(jí)順序，將藤茶葉片轉(zhuǎn)錄組Unigene進(jìn)行BLAST比對(duì)。從比對(duì)上的序列中提取出Transcript開放閱讀框（ORF），并將CDS序列翻譯成氨基酸序列。未比對(duì)上或比對(duì)上未預(yù)測(cè)出結(jié)果的序列則采用ESTScan 3.0.3預(yù)測(cè)其ORF，根據(jù)CDS序列翻譯獲得氨基酸序列。

2 結(jié)果與分析

2. 1 轉(zhuǎn)錄組測(cè)序及組裝結(jié)果

藤茶葉片轉(zhuǎn)錄組測(cè)序獲得82126236條Raw reads，過(guò)濾處理后得到80156972條Clean reads，有效測(cè)序長(zhǎng)度達(dá)12.02 Gb，堿基錯(cuò)誤率為0.02%，Q20和Q30分別為97.30%和92.97%，GC含量為43.36%，說(shuō)明測(cè)序質(zhì)量好，可用于后續(xù)分析。

采用Trinity將所有Clean reads拼接得到153528條轉(zhuǎn)錄本，平均長(zhǎng)度為841 bp，N50為1541 bp，其中，以長(zhǎng)度200～500 bp的序列數(shù)量最多，為83963條，占54.69%;以長(zhǎng)度>2000 bp的序列數(shù)量最少，為16368條，占10.66%。將轉(zhuǎn)錄本進(jìn)一步組裝拼接得到92472條Unigenes，平均長(zhǎng)度為1208 bp，N50為1780 bp，其中，以長(zhǎng)度500～1000 bp的序列數(shù)量最多，為26916條，占29.11%;以長(zhǎng)度>2000 bp的序列數(shù)量最少，為16368條，占17.70%（表1）。

2. 2 Unigene功能注釋結(jié)果

將92472條Unigenes分別在7個(gè)數(shù)據(jù)庫(kù)進(jìn)行比對(duì)注釋，每個(gè)數(shù)據(jù)庫(kù)中注釋成功的Unigenes數(shù)量及其占Unigenes總數(shù)的比例如表2所示。在7個(gè)數(shù)據(jù)庫(kù)中，至少在1個(gè)數(shù)據(jù)庫(kù)注釋的Unigenes有84217條，占Unigenes總數(shù)的91.07%，其中，在Nt數(shù)據(jù)庫(kù)中注釋的Unigenes最多（80266條），占總數(shù)的86.80%，其次是Nr數(shù)據(jù)庫(kù)（58148條），占Unigenes總數(shù)的62.88%，KOG數(shù)據(jù)庫(kù)中注釋的Unigenes數(shù)量最少（14553條），占Unigenes總數(shù)的15.73%;有8944條Unigenes在7個(gè)數(shù)據(jù)庫(kù)均被注釋，占Unigenes總數(shù)的9.67%。

基于Nr數(shù)據(jù)庫(kù)比對(duì)注釋的結(jié)果，統(tǒng)計(jì)并繪制比對(duì)上的物種分布情況（圖1），發(fā)現(xiàn)注釋到葡萄的Unigenes數(shù)量最多，所占比例高達(dá)88.17%，遠(yuǎn)高于可可（1.17%）、蓮（1.11%）、麻風(fēng)樹（0.69%）、甜橙（0.60%）及其他物種（8.27%），究其原因可能是藤茶與葡萄同屬于葡萄科植物，親緣關(guān)系更近。

2. 3 Unigene的GO功能分類結(jié)果

從圖2可知，在GO數(shù)據(jù)庫(kù)成功注釋的41116條Unigenes可分為生物學(xué)過(guò)程、細(xì)胞組分和分子功能三大類，共56個(gè)小類。由于存在同一個(gè)Unigene注釋到不同功能的情況，有106044條Unigenes注釋歸屬于生物學(xué)過(guò)程，其中，以細(xì)胞過(guò)程、代謝過(guò)程和單有機(jī)體過(guò)程最多;有66390條Unigenes注釋歸屬于細(xì)胞組分，其中以細(xì)胞和細(xì)胞部分最多;有51272條Unigenes注釋歸屬于分子功能，其中，主要以結(jié)合和催化活性為主。

2. 4 Unigene的KOG功能分類結(jié)果

通過(guò)KOG注釋可將14553條Unigenes分成25類（圖3），其中，一般功能預(yù)測(cè)注釋成功的Unigenes最多（1946條）;其次是翻譯后修飾、蛋白質(zhì)翻轉(zhuǎn)、分子伴侶（1776條）;翻譯、核糖體結(jié)構(gòu)和生物起源注釋成功的Unigenes有1250條;參與次生代謝物質(zhì)的生物合成、轉(zhuǎn)運(yùn)和降解的Unigenes較少，僅有319條。此外，有988條未知功能的Unigenes有待進(jìn)一步研究。

2. 5 Unigene的KEGG信號(hào)通路富集分析結(jié)果

利用KEGG數(shù)據(jù)庫(kù)中的通路數(shù)據(jù)庫(kù)（Pathway databases）對(duì)Unigene的功能進(jìn)行富集分析（表3），結(jié)果發(fā)現(xiàn)共有15262條Unigenes能注釋到128條KEGG信號(hào)通路。這些通路可分成細(xì)胞進(jìn)程、環(huán)境信號(hào)處理、遺傳信息處理、代謝和有機(jī)系統(tǒng)五大類。其中，注釋為代謝的Unigenes最多，為8694條，占注釋Unigenes總數(shù)的56.97%，其次是注釋為遺傳信息處理的Unigenes，為4120條，占注釋Unigenes總數(shù)的26.99%。

KEGG信號(hào)通路富集分析結(jié)果（表3）表明，在五大類通路中，與代謝有關(guān)的通路最多，有98條，占所有通路的76.56%，其中，Unigene注釋最多為碳水化合物（1671條），其次是氨基酸代謝（1176條）。有512條Unigenes參與苯丙素類、生物堿、黃酮、花青素等生物合成相關(guān)的11個(gè)其他次生代謝標(biāo)準(zhǔn)通路。

2. 6 Unigene的CDS序列預(yù)測(cè)結(jié)果

按照Swiss-Prot和Nr數(shù)據(jù)庫(kù)的優(yōu)先級(jí)順序，將藤茶葉片轉(zhuǎn)錄組Unigene在這2個(gè)蛋白數(shù)據(jù)庫(kù)進(jìn)行BLAST比對(duì)，獲得52582條CDS序列;將未比對(duì)上或比對(duì)上未預(yù)測(cè)出結(jié)果的序列采用ESTScan 3.0.3進(jìn)行CDS預(yù)測(cè)，結(jié)果獲得35535條CDS序列。在獲得的88117條CDS序列中，以長(zhǎng)度0～200 nt的CDS序列最多，有52828條;其次是長(zhǎng)度201～400和401～600 nt的CDS序列，分別為20564和8796條，其余長(zhǎng)度范圍的CDS序列均不足4000條（圖4）。

2. 7 黃酮類物質(zhì)生物合成相關(guān)基因挖掘及分析結(jié)果

黃酮類化合物的主要合成途徑：首先通過(guò)苯丙烷途徑將苯丙氨酸轉(zhuǎn)化為香豆酰-CoA，然后香豆酰-CoA進(jìn)入黃酮合成途徑，與3分子丙二酰CoA結(jié)合生成查爾酮，再經(jīng)分子內(nèi)的環(huán)化反應(yīng)生成二氫黃酮類化合物，最后通過(guò)分支途徑合成黃酮、異黃酮、黃酮醇、黃烷醇和花色素等黃酮類化合物（鄒麗秋等，2016）?；贙EGG信號(hào)通路富集分析結(jié)果，共篩選獲得98個(gè)藤茶黃酮類化合物合成相關(guān)基因（表4），其中18個(gè)映射到苯丙烷代謝通路上，分別編碼苯丙氨酸向香豆酰-CoA轉(zhuǎn)化所需的3種關(guān)鍵酶：苯丙氨酸解氨酶（PAL）、肉桂酸-4-羥化酶（C4H）和對(duì)香豆酸—輔酶A連接酶（4CL）;有80個(gè)Unigenes映射到類黃酮類化合物代謝通路上，編碼該通路中的14種關(guān)鍵酶，包括從香豆酰-CoA向二氫黃酮轉(zhuǎn)化所需的查爾酮合成酶（CHS）和查耳酮異構(gòu)酶（CHI），以及從二氫黃酮到黃酮醇和黃烷醇合成支路所需的黃烷酮-3-羥化酶（F3H）、類黃酮-3-羥化酶（F3H）、類黃酮-3，5-羥化酶（F35H）、黃酮醇合成酶（FLS）、二氫黃酮醇還原酶（DFR）、花青素還原酶（ANR）、無(wú)色花色素還原酶（LAR）、無(wú)色花色素雙加氧酶（LDOX）、香豆酰酯3羥化酶（C3H）、咖啡酰輔酶A-O甲基轉(zhuǎn)移酶（CCOAOMT）、黃酮醇3-甲基轉(zhuǎn)移酶（F3OMT）和莽草酸羥基肉桂轉(zhuǎn)移酶（HCT）共12種。

3 討論

本研究利用第二代高通量測(cè)序技術(shù)對(duì)藤茶葉片轉(zhuǎn)錄組進(jìn)行測(cè)序分析，共獲得92472條Unigenes。與其他藥用植物比較，藤茶轉(zhuǎn)錄組測(cè)序得到的Unigenes多于銀杏（劉偉等，2018）、大葉蛇葡萄（Yang et al.，2019）和金線蓮（鄒福賢等，2019）等，與山茱萸（朱畇昊等，2017）和掌葉大黃（李歡等，2018）相當(dāng)，究其原因可能與物種基因型和取材等因素有關(guān)。從Unigene長(zhǎng)度來(lái)看，本研究得到Unigene的平均長(zhǎng)度為1208 bp，N50為1780 bp，二者均明顯高于連翹（861 bp，1497 bp）（王興春等，2015）、草麻黃（453 bp，525 bp）（馬婧等，2016）及短絲木犀（697 bp，1200 bp）（陳林等，2016）等，與櫻椒樹（1116 bp，1768 bp）（Zhou et al.，2015）和掌葉大黃（1108 bp，1794 bp）（李歡等，2018）等相近。N50是評(píng)價(jià)組裝序列完整性的重要指標(biāo)，其數(shù)值越大表示組裝得到的長(zhǎng)片段越多，當(dāng)N50≥800 bp說(shuō)明組裝得到序列完整性較好（胡俊杰等，2017）。本研究中N50>800 bp，表明藤茶轉(zhuǎn)錄組組裝后得到的序列完整性較好，能滿足轉(zhuǎn)錄組分析的要求，從而保證功能注釋能得到較理想的比對(duì)結(jié)果。

本研究將獲得的92472條Unigenes在7個(gè)Nr數(shù)據(jù)庫(kù)中進(jìn)行比對(duì)注釋，結(jié)果顯示，有84217條Unigenes（91.07%）至少在1個(gè)數(shù)據(jù)庫(kù)被注釋，有8944條Unigenes（9.67%）在7個(gè)數(shù)據(jù)庫(kù)均被注釋，且有20956條（22.66%）、41116條（44.46%）和14553條（15.73%）分別在KO、GO和KOG數(shù)據(jù)庫(kù)被注釋。從本研究Nr數(shù)據(jù)庫(kù)同源比對(duì)匹配的物種分布結(jié)果可知，注釋到葡萄的Unigenes最多，所占比例高達(dá)88.17%。與其他藥用植物相比，藤茶轉(zhuǎn)錄組Unigenes的注釋率（91.07%）明顯高于草麻黃（81.81%）（馬婧等，2016）、銀杏（81.67%）（劉偉等，2018）和夏枯草（53.13%）（朱畇昊等，2019）等，可能是因?yàn)樘俨铻槠咸训慕壩锓N，在現(xiàn)有數(shù)據(jù)庫(kù)葡萄的蛋白序列信息較豐富，導(dǎo)致藤茶Unigenes注釋率也較高。

黃酮類化合物的種類較多、分布較廣，其生物合成途徑目前已基本清楚，各步驟催化反應(yīng)的所需酶及其基因均有大量研究（Hichri et al.，2011;鄒麗秋等，2016）。藤茶是富含黃酮類化合物的一種天然植物，其嫩莖葉中的總黃酮含量高達(dá)43.4%～45.52%，遠(yuǎn)高于其他植物，且藤茶中的黃酮類化合物種類多樣，除二氫楊梅素外，還含有其他多種黃酮類成分，如槲皮素、花旗松素、山奈酚、藤茶素、楊梅苷、橙皮素、洋芹素、大黃素和阿福豆素等（Zhang and Yang，2001;何桂霞等，2007;付明等，2015）。本研究對(duì)藤茶轉(zhuǎn)錄組Unigenes進(jìn)行KEGG信號(hào)通路富集分析，結(jié)果發(fā)現(xiàn)在五大類通路中，有98條與代謝有關(guān)的通路，占所有通路的76.56%，通過(guò)進(jìn)一步篩選，共獲得98個(gè)黃酮類化合物合成相關(guān)的基因，分別編碼17種關(guān)鍵酶，包括3個(gè)苯丙烷代謝途徑的所需酶（PAL、4CL和C4H）和14個(gè)類黃酮代謝途徑的所需酶（CHS、CHI、F3H、F3H、F35H、FLS、DFR、ANR、LDOX、LAR、F3OMT、HCT、C3H和CCOAOMT），為深入研究藤茶黃酮生物合成途徑及其相關(guān)關(guān)鍵基因的篩選、功能鑒定及調(diào)控機(jī)制打下了基礎(chǔ)。今后還將通過(guò)實(shí)時(shí)熒光定量PCR等技術(shù)檢測(cè)上述關(guān)鍵酶基因的差異表達(dá)情況，同時(shí)結(jié)合藥效成分含量分析，解析上述基因在藤茶藥效成分生物合成過(guò)程的功能。

4 結(jié)論

藤茶在細(xì)胞過(guò)程、代謝過(guò)程、單有機(jī)體過(guò)程、細(xì)胞和細(xì)胞部分、結(jié)合和催化活性能力分布的基因較豐富，在一般功能、翻譯、翻譯后修飾、蛋白質(zhì)翻轉(zhuǎn)及分子伴侶的基因表達(dá)量較高，具有較強(qiáng)的碳水化合物代謝能力。多種關(guān)鍵酶基因參與藤茶黃酮類化合物的生物合成，推測(cè)其生物合成途徑存在多條分支，調(diào)控機(jī)制也較復(fù)雜。

參考文獻(xiàn)：

陳林，李龍娜，戴亞平，楊國(guó)棟. 2016. 短絲木犀轉(zhuǎn)錄組測(cè)序及類胡蘿卜素生物合成相關(guān)基因表達(dá)分析[J]. 南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），40（5）：21-28. [Chen L，Li L N，Dai Y P，Yang G D. 2016. De novo transcriptome sequencing and analysis of carotenoids biosynthesis related gene expression in Osmanthus serrulatus[J]. Journal of Nanjing Forestry University（Natural Sciences Edition），40（5）：21-28.]

付明，黎曉英，王登宇，郭孟齊. 2015. 顯齒蛇葡萄葉中黃酮類化合物的研究[J]. 中國(guó)藥學(xué)雜志，50（7）：574-578. [Fu M，Li X Y，Wang D Y，Guo M Q. 2015. Flavonoid constituents of leaves of Ampelopsis grossedentata（Hand-Mazz） W. T. wang[J]. Chinese Pharmaceutical Journal，50（7）：574-578.]

付明，魏麟，余娟，余小林. 2013. 顯齒蛇葡萄查耳酮合成酶基因cDNA克隆及蛋白質(zhì)序列分析[J]. 中草藥，44（1）：85-89. [Fu M，Wei L，Yu J，Yu X L. 2013. cDNA cloning and protein sequence analysis of chalcone synthase gene in leaves of Ampelopsis grossedentata[J]. Chinese Traditional and Herbal Drugs，44（1）：85-89.]

何桂霞，裴剛杜，方麓，歐陽(yáng)文，李斌. 2007. 藤茶化學(xué)成分的研究[J]. 中國(guó)現(xiàn)代中藥，9（12）：11-13. [He G X，Pei G D，F(xiàn)ang L，Ouyang W，Li B. 2007. Studies on chemical constituents of Ampelopsis grossdentata[J]. Modern Chinese Medicine，9（12）：11-13.]

胡俊杰，孟翔，周佳濱，楊陸興，劉善海，李潤(rùn)釗. 2017. 扶桑綿粉蚧轉(zhuǎn)錄組分析[J]. 昆蟲學(xué)報(bào)，60（1）：9-17. [Hu J J，Meng X，Zhou J B，Yang L X，Liu S H，Li R Z. 2017. Transcriptome analysis of the cotton mealybug，Phenacoccus solenopsis（Hemiptera：Pseudococcidae）[J]. Acta Entomologica Sinica，60（1）：9-17.]

黃海波. 2008. 藤茶PAL特性及其與黃酮和DMY含量變化相關(guān)性的研究[D]. 武漢：華中農(nóng)業(yè)大學(xué). [Huang H B. 2008. Studies on the characters of PAL and its correlation between flavones and DMY on contents diversification in Tengcha[D]. Wuhan：Huazhong Agricultural University.]

李歡，張娜，李依民，黑小斌，李元敏，鄧翀，顏永剛，劉蒙蒙，張崗. 2018. 利用轉(zhuǎn)錄組測(cè)序挖掘掌葉大黃蒽醌類生物合成相關(guān)基因[J]. 藥學(xué)學(xué)報(bào)，53（11）：1908-1917. [Li H，Zhang N，Li Y M，Hei X B，Li Y M，Deng C，Yan Y G，Liu M M，Zhang G. 2018. High-throughput transcripto-mic sequencing of Rheum palmatum L. seedlings and elucidation of genes in anthraquinone biosynthesis[J]. Acta Pharmaceutica Sinica，53（11）：1908-1917.]

劉偉，王俊燚，李萌，董金金，何崇單，汪貴斌，郁萬(wàn)文，王義強(qiáng). 2018. 基于轉(zhuǎn)錄組測(cè)序的銀杏類黃酮生物合成關(guān)鍵基因表達(dá)分析[J]. 中草藥，49（23）：5633-5639. [Liu W，Wang J Y，Li M，Dong J J，He C D，Wang G B，Yu W W，Wang Y Q. 2018. Transcriptome sequencing analysis of gene expression of flavonoid biosynthesis in Ginkgo biloba[J]. Chinese Traditional and Herbal Drugs，49（23）：5633-5639.]

馬婧，成鐵龍，孫燦岳，鄧楠，史勝青，江澤平. 2016. 草麻黃高通量轉(zhuǎn)錄組分析及黃酮類代謝途徑相關(guān)基因的鑒定[J]. 浙江農(nóng)業(yè)學(xué)報(bào)，28（4）：609-617. [Ma J，Cheng T L，Sun C Y，Deng N，Shi S Q，Jiang Z P. 2016. Characterization of transcriptome reveals pathway of flavonoids in Ephedra sinica Stapf[J]. Acta Agriculturae Zhejiangensis，28（4）：609-617.]

冉京燕，方建國(guó)，謝雪佳，熊微，王文清. 2016. 藤茶的本草資源學(xué)研究概況[J]. 中草藥，47（20）：3728-3735. [Ran J Y，F(xiàn)ang J G，Xie X J，Xiong W，Wang W Q. 2016. Scientific research on herbal resource of vine tea[J]. Chinese Traditional and Herbal Drugs，47（20）：3728-3735.]

譚沙，羅靜，李建新，宋珊珊. 2015. 藤茶研究進(jìn)展[J]. 安徽農(nóng)業(yè)科學(xué)，43（31）：71-75. [Tan S，Luo J，Li J X，Song S S. 2015. Research progress of vine tea[J]. Journal of Anhui Agricultural Sciences，43（31）：71-75.]

王家勝，何磊磊，張妮，孔琪，玉秋萍，余正文. 2014. 不同產(chǎn)地顯齒蛇葡萄中二氫楊梅素測(cè)定[J]. 中成藥，36（1）：145-147. [Wang J S，He L L，Zhang N，Kong Q，Yu Q P，Yu Z W. 2014. Dihydromyricetin in Ampelosis grossedentata collected from different habitats[J]. Chinese Traditional Patent Medicine，36（1）：145-147.]

王興春，譚河林，陳釗，孟令芝，王文斌，范圣此. 2015. 基于RNA-Seq技術(shù)的連翹轉(zhuǎn)錄組組裝與分析及SSR分子標(biāo)記的開發(fā)[J]. 中國(guó)科學(xué)（生命科學(xué)），45（3）：301-310. [Wang X C，Tan H L，Chen Z，Meng L Z，Wang W B，F(xiàn)an S C. 2015. Assembly and characterization of the transcriptome and development of SSR markers in Forsythia suspensa based on RNA-Seq technology[J]. Scientia Sinica（Vitae），45（3）：301-310.]

許明，楊志堅(jiān)，鄭金貴. 2015. 藤茶葉片總RNA的提取及苯丙氨酸解氨酶基因（pal）片段的克隆[J]. 生物學(xué)雜志，32（2）：96-99. [Xu M，Yang Z J，Zheng J G. 2015. Extraction of total RNA from leaves of Ampelopsis grossedentata and cloning of phenylalanine ammonia lyase gene fragment[J]. Journal of Biology，32（2）：96-99.]

許明，伊恒杰，趙帥，張玉文，楊志堅(jiān)，鄭金貴. 2017. 顯齒蛇葡萄實(shí)時(shí)熒光定量PCR內(nèi)參基因的篩選與驗(yàn)證[J]. 中草藥，6（3）：1192-1198. [Xu M，Yi H J，Zhao S，Zhang Y W，Yang Z J，Zheng J G. 2017. Selection and validation of reference genes for quantitative RT-PCR analysis in Ampelopsis grossedentata[J]. Chinese Traditional and Her-bal Drugs，6（3）：1192-1198.]

嚴(yán)志祥，羅茜，張翼冠，趙軍寧. 2019. 藥用植物轉(zhuǎn)錄組研究現(xiàn)狀與展望[J]. 中國(guó)藥學(xué)雜志，54（7）：513-520. [Yan Z X，Luo X，Zhang Y G，Zhao J N. 2019. A brief review on medicinal plants transcriptome research[J]. Chinese Pharmaceutical Journal，54（7）：513-520.]

朱畇昊，董誠(chéng)明，鄭曉珂，馮衛(wèi)生，劉孟奇，趙樂(lè). 2017. 基于轉(zhuǎn)錄組測(cè)序的山茱萸次生代謝生物合成相關(guān)基因的挖掘[J]. 中國(guó)中藥雜志，42（2）：213-219.[Zhu Y H，Dong C M，Zheng X K，F(xiàn)eng W S，Liu M Q，Zhao L. 2017. Transcriptome analysis reveals genes involved in biosynthesis of secondary metabolism in Cornus officinalis[J]. China Journal of Chinese Materia Medica，42（2）：213-219.]

朱畇昊，張夢(mèng)佳，李璐，趙樂(lè)，董誠(chéng)明. 2019. 夏枯草的轉(zhuǎn)錄組測(cè)序與次生代謝產(chǎn)物生物合成相關(guān)基因的挖掘[J]. 中草藥，50（5）：1220-1226. [Zhu Y H，Zhang M J，Li L，Zhao L，Dong C M. 2019. Transcriptome analysis of Prunella vulgaris and identification of putative genes involved in secondary metabolism biosynthesis[J]. Chinese Traditional and Herbal Drugs，50（5）：1220-1226.]

鄒福賢，許文，黃澤豪，張勛，陳抒云，林羽，徐偉. 2019. 金線蓮轉(zhuǎn)錄組測(cè)序及其黃酮類合成相關(guān)基因分析[J]. 中國(guó)藥科大學(xué)學(xué)報(bào)，50（1）：66-74. [Zou F X，Xu W，Huang Z H，Zhang X，Chen S Y，Lin Y，Xu W. 2019. Analysis of transcriptome sequencing and related genes of flavonoid biosynthesis from Anoectochilus roxburghii[J]. Journal of China Pharmaceutical University，50（1）：66-74.]

鄒麗秋，王彩霞，匡雪君，李瀅，孫超. 2016. 黃酮類化合物合成途徑及合成生物學(xué)研究進(jìn)展[J]. 中國(guó)中藥雜志，41（22）：4124-4128. [Zou L Q，Wang C X，Kuang X J，Li Y，Sun C. 2016. Advance in flavonoids biosynthetic pathway and synthetic biology[J]. China Journal of Chinese Materia Medica，41（22）：4124-4128.]

Grabherr M G，Haas B J，Yassour M，Levin J Z，Thompson D A，Amit I，Adiconis X，F(xiàn)an L，Raychowdhury R，Zeng Q，Chen Z，Mauceli E，Hacohen N，Gnirke A，Rhind N，di Palma F，Birren B W，Nusbaum C，Lindblad-Toh K，F(xiàn)riedman N，Regev A. 2011. Full-length transcriptome assembly from RNA-Seq data without are ference genome[J]. Nature Biotechnology，29（7），644-652.

Hichri I，Barrieu F，Bogs J，Kappel C，Delrot S，Lauvergeat V. 2011. Recent advances in the transcriptional regulation of the flavonoid biosynthetic pathway[J]. Journal of Experimental Botany，62（8）：2465-2483.

Wilhelm B T，Landry J R. 2009. RNA-Seq-quantitative measurement of expression through massively parallel RNA-sequencing[J]. Methods，48（3）：249-257.

Yang M，Zhou P N，Gui C，Da G Z，Gong L，Zhang X Q. 2019. Comparative transcriptome analysis of Ampelopsis megalophylla for identifying genes involved inflavonoid biosynthesis and accumulation during different seasons[J]. Molecules，24（7）：1267-1282.

Zhang Y S，Yang W L. 2001. Study on functional constituents of Ampelopsis grossdentata[J]. Food Science，22（9）：75-77.

Zhang Y S，Zhang Q Y，Wang B，Li L Y，Zhao Y Y. 2006. Chemical constituents from Ampelopsis grossedentata[J]. Journal of Chinese Pharmaceutical Sciences，15（4）：211-214.

Zhou X J，Wang Y Y，Xu Y N，Yan R S，Zhao P，Liu W Z. 2015. De Novo characterization of flower bud transcriptomes and the development of EST-SSR markers for the endangered tree Tapiscia sinensis[J]. International Journal of Molecular Sciences，16（6）：12855-12870.

（責(zé)任編輯 陳 燕）

收稿時(shí)間：2019-08-15

基金項(xiàng)目：國(guó)家科技支撐計(jì)劃項(xiàng)目（2013BAD01B05）;福建農(nóng)林大學(xué)科技創(chuàng)新專項(xiàng)（KFA17424A）

作者簡(jiǎn)介：*為通訊作者，鄭金貴（1949-），教授，主要從事地方特色作物種質(zhì)資源挖掘與創(chuàng)新利用研究工作，E-mail：jgzheng@ fafu.edu.cn。許明（1978-），博士，主要從事藥用植物生物技術(shù)研究工作，E-mail：xmfau@163.com

優(yōu) 秀 青 年 學(xué) 者 論 壇

許明（1978-），博士，主要從事藥用植物生物技術(shù)研究工作。先后主持或作為主要成員參與國(guó)家科技支撐計(jì)劃項(xiàng)目“地方特色作物種質(zhì)資源發(fā)掘與創(chuàng)新利用”、國(guó)家轉(zhuǎn)基因重大專項(xiàng)“功能型轉(zhuǎn)基因水稻新品種的培育”、農(nóng)業(yè)農(nóng)村部“948”計(jì)劃項(xiàng)目“臺(tái)灣優(yōu)質(zhì)農(nóng)產(chǎn)品技術(shù)的引進(jìn)、創(chuàng)新與推廣”和福建省特色現(xiàn)代農(nóng)業(yè)項(xiàng)目“特種作物新品種試驗(yàn)篩選與功能品質(zhì)鑒定”等國(guó)家級(jí)、省部級(jí)科研項(xiàng)目10余項(xiàng);在《Molecular Genetics and Genomics》《International Journal of Molecular Sciences》《Iranian Journal of Biotechno-logy》《中國(guó)農(nóng)業(yè)科學(xué)》《南方農(nóng)業(yè)學(xué)報(bào)》《西北植物學(xué)報(bào)》《中草藥》《生物技術(shù)通報(bào)》等國(guó)內(nèi)外核心期刊上發(fā)表學(xué)術(shù)論文20余篇;申報(bào)國(guó)家發(fā)明專利4項(xiàng)，獲福建省科學(xué)技術(shù)獎(jiǎng)三等獎(jiǎng)1項(xiàng)。