五指毛桃为桑科榕属植物粗叶榕（Ficus hirta Vahl）,其掌状叶形如五根手指,果实状如毛桃,因此被称为五指毛桃,又因其根皮具有牛奶的清香,所以又有“五指牛奶”的美称,有文献报道同属植物极简榕（Ficus simplissima Lour.）的干燥根也作为五指毛桃,是岭南地区习用中草药[1]。五指毛桃为地道的民族用药,广东、广西、福建、湖南、云南等省（区）14个少数民族均有用药记载,目前产区主要为广西、广东等地,以野生供应市场需求,其用药部位为根茎。中医认为,五指毛桃性平、味甘、辛,对脾虚浮肿、肺痨咳嗽、食少无力、风湿痹痛、盗汗、产后无乳等症有很好的疗效[2]。目前已从五指毛桃中鉴定出多种化学成分,如香豆素类、黄酮类、甾醇类等[3,4,5],具有抗氧化、抗炎、抗菌、抗病毒和抗肿瘤等作用[6,7],同时还能提高机体的免疫能力[8,9];通过对其挥发油类化合物进行测定,发现十六酸、油酸、亚油酸等为其主要含有的挥发油类成分[10]。黄酮类化合物是一种重要的天然有机化合物,具有抗氧化、抗癌症、抗肿瘤以及提高机体免疫机能等方面的功效,其为五指毛桃活性物质中数量最多的一类,如芹菜素、橙皮苷、甲氧基黄酮[11,12]、牡荆苷[13]、柚皮素[3]、槲皮素、金合欢素[4]等,具有较高的研究和药用价值。
随着高通量转录组测序技术的发展,不仅在模式植物及重要的农作物上具有普遍的运用,近几年也在药用植物上得到了广泛的运用[14]。该技术可在缺少参考基因组情况下全面分析某一物种的转录水平,从而确定转录本及其核酸序列,为理解物种的生物学及生物化学等方面提供技术手段[15]。目前,国内外在西洋参[16]、甘草[17]、丹参[18]、金线莲[19]、半枝莲[20]、穿心莲[21]等药用植物开展转录组研究,从中挖掘和分离出了大量与次生代谢产物相关的基因。目前,关于五指毛桃的研究主要集中在生物学特征、活性成分的测定、药理作用及栽培种植等方面,然而,五指毛桃的转录组和功能基因组等分子水平上的研究较少报道。本研究利用Illumina HiSeq X Ten测序平台对五指毛桃的根、茎和叶进行了大规模的转录组学分析,阐述了五指毛桃不同组织部位的生物学过程和代谢调节相关途径及基因表达的变化,重点分析了黄酮类化合物生物合成途径中关键基因的表达特性,以期为五指毛桃活性成分的生物合成与调控等分子机制的研究及次生代谢工程的发展提供理论依据。
1 材料与方法
1.1 材料
野外采集五指毛桃种子,植株特征为：灌木,嫩枝中空,小枝、叶、果实均被金黄色长硬毛,有乳汁。叶片通常掌状分裂;瘦果椭圆形。植株经漳州卫生职业学院李珍讲师鉴定为粗叶榕（Ficus hirta Vahl）。种子委托广西新众农业有限公司育苗,得到2月龄五指毛桃小苗。将植株的根、茎、叶分开采集后迅速置于液氮中冻存,然后保存于‒80 ℃冰箱用于后续转录组测序及相关试验。
1.2 方法
1.2.1 RNA提取 采用Trizol（Invitrogen, CA, USA）试剂盒进行样品总RNA的提取,并用DNase I（Takara Bio, Japan）消化基因组中的DNA。总RNA的浓度和质量分别用安捷伦2100型（Agilent Technologies, USA）生物分析仪测定,并用琼脂糖凝胶电泳评估总RNA的完整性。选取质量、浓度和完整性符合要求的总RNA用于后续cDNA文库的构建。
1.2.2 cDNA文库的构建和高通量测序 首先将带有Oligo dT（Qiagen, China）磁珠与mRNA完全结合,纯化mRNA,然后加入打断剂（Frag/Rrime buffer）,用移液枪反复吸打混匀后置于PCR仪中进行mRNA片段化。以该片段为模板,用随机六聚体引物和SuperScript III（Qiagen, China）逆转录酶进行第一链cDNA的合成,然后以第一链cDNA为模板,RNase H和DNA聚合酶I（Qiagen, China）进行第二链cDNA的合成。用QIAquick PCR提取试剂盒（Qiagen, China）将双链cDNA纯化及末端修复/dA尾添加,接着用T4连接酶进行接头连接,然后将连接产物纯化、片段大小分选和文库扩增。用凝胶电泳、安捷伦2100型（Agilent Technologies, USA）生物分析仪和定量PCR仪（T100™ Thermal Cycler, BIO-RAD）检测所构建的文库质量,最后,用Illumina HiSeq X Ten测序平台对合格的文库进行RNA-seq,以获得原始数据（Raw reads）,该高通量测序由生工生物工程（上海）股份有限公司完成。
1.2.3 质量控制与De novo组装 将测序获得的原始数据（Raw reads）进行质量控制,去除接头序列、低质碱基（Q值<20）、含量超过10%的N碱基的序列,获得高质量的Clean reads。使用Trinity程序将Clean数据De novo组装成转录本,对拼接获得的转录本去除冗余,取每个转录本聚类最长的转录本作为unigene,用于后续的分析。
1.2.4 Unigene功能注释 将获得的unigene用NCBI Blastx后比对到蛋白质数据库Nr、KOG、Swiss-Prot、TrEMBL、CDD、PFAM和KEGG（E-value≤0.00001）,获得与序列相似度最高的蛋白质进行五指毛桃unigene的功能注释。用topGO进行GO注释和分类[22],并将unigene用NCBI Blastx后与NR数据库进行比对,获得与五指毛桃相近的物种以及同源序列的功能信息。
1.2.5 差异表达基因的筛选和GO和KEGG富集分析 利用RSEM（RNA-seq by Expectation Maximization）法计算基因的表达量,并用TPM（Transcripts Per Million）方法对相关基因表达水平进行归一化处理。再采用DEseq进行差异分析,并对统计得到P值进行FDR矫正获得q值,为了得到显著差异的基因,将筛选条件设为q<0.005且差异倍数|FoldChange|>2。分别使用topGO（http://www.bioconductor.org/packages/release/bioc/html/topGO.html）和clusterProfiler（http://www. bioconductor.org/packages/release/bioc/html/clusterProfiler.html）进行差异表达基因的GO和KEGG富集分析,且认为P<0.05为显著富集。
1.2.6 五指毛桃unigene SSR和SNP/InDel分析 使用MISA（https://webblast.ipk-gatersleben.de/ misa/）对五指毛桃根、茎和叶基因转录本进行SSR检测,找出SSR标记之后,使用Primer3完成SSR引物设计,参数均设为默认。使用BCFtools（http://www.htslib.org/doc/bcftools.html）进行SNP/InDel calling,抽提各样本中的SNP/InDel。并对其进行质量值大于20、覆盖度大于8的过滤。
1.2.7 黄酮类化合物生物合成关键基因的RT- qPCR验证 为了进行黄酮类化合物生物合成关键基因的RT-qPCR验证,从五指毛桃根、茎、叶的每个样品中提取总RNA 500 ng,利用SYBR Ex Script™试剂盒（TaKaRa, China）,用随机引物和Oligo dT引物逆转录成cDNA。RT-qPCR引物参照引物设计原则,利用DNAMAN 6.0软件进行设计（表1）,引物送北京六合华大科技有限公司合成。RT-qPCR操作方法按照试剂盒（TaKaRa, China）操作说明进行,参照同属植物桑树的相关研究以ACT3为内参基因[23],通过罗氏Light Cyclers荧光定量PCR仪进行,采用RT-qPCR检测9个关键基因的相对表达量,并用2‒ΔΔCT法进行计算,利用Graphpad Prism软件作图。
表1 实时荧光定量PCR所用引物Tab. 1 Sequences of specific primers used for quantitative real-time PCR |
| 基因ID Gene ID | 基因名称 Gene name | 引物序列（5°‒3°） Prime sequence (5°‒3°) | 产物大小 Product size/bp | 退火温度 Tm/℃ |
|---|---|---|---|---|
| TRINITY_DN34248_c1_g1 | DFR | F: CTCCTAGAGCGCGGTTACAC | 149 | 60 |
| R: CTTCACTAAAGCTTCCCTCCAC | ||||
| TRINITY_DN33113_c1_g1 | CHS | F: GACTGTCAAGGTGGGTGTCA | 101 | 60 |
| R: GGACCTCGCTGAGAACAACA | ||||
| TRINITY_DN33700_c3_g1 | CHI | F: ACCTACCATTTCCCTGCAGC | 146 | 60 |
| R: GACAGAGCTGAAGACCCACC | ||||
| TRINITY_DN30694_c0_g4 | C4H | F: GCCATCGTGATCTCAAAGC | 161 | 59 |
| R: ATGCGGAGGAGGAATACGT | ||||
| TRINITY_DN34559_c2_g3 | HST | F: TGATGCTGACTTTGGATGG | 107 | 59 |
| R: ACAAACTCCCGTCGTTGG | ||||
| TRINITY_DN29797_c0_g2 | CCOMT | F: TTGGAAGGACTGGGTGGTG | 118 | 60 |
| R: ACCCGTCTCATCTGCTTCTG | ||||
| TRINITY_DN34838_c3_g2 | F3'H | F: CGGTCGAGTTCTTCTTGGA | 167 | 59 |
| R: AATGCTGATGGTGAGGGTG | ||||
| TRINITY_DN27691_c0_g1 | LDOX | F: TTCTCCTTGTTGACGAGCC | 179 | 58 |
| R: GCCCTCACCTTCATCCTAC | ||||
| TRINITY_DN31152_c2_g1 | AAT | F: TGAAACTGCGAACCAAACAG | 133 | 60 |
| R: TTGAGGACGGAAGACATGG | ||||
| TRINITY_DN34055_c1_g1 | ACT3 | F: TTAGACGTTGTGTAGCCAGCA | 174 | 60 |
| R: CTACTTGCGTTGTGATCCGG |
2 结果与分析
2.1 转录组测序质量分析和De novo组装
利用RNA-seq技术对五指毛桃的根、茎、叶进行转录组测序,共生成138 712 412个高质量reads,Q20碱基百分比（测序错误率小于1%）分布介于98.58%~98.63%之间,Q30碱基百分比分布（测序错误率小于0.1%）介于94.92%~95.05%之间,GC碱基含量介于47.89%~49.17%之间（表2）。质量评估结果表明五指毛桃的转录组测序数据质量较高,可用于后续分析。使用Trinity将高质量数据进行De novo组装成转录本,重新组装成104 335条unigenes,平均长度为578 bp,其中N50长度为1025 bp,unigene序列长度分布见图1,这些获得的转录本可作为五指毛桃参考的转录本。
表2 五指毛桃转录组数据统计和质量评估Tab. 2 Statistics and raw reads and clean reads in transcriptome of F. hirta Vahl |
| 样品 Sample | 原始数据Raw data | 高质量数据Clean data | ||||
|---|---|---|---|---|---|---|
| 读取总数 Total reads count | 碱基总数 Total bases count/bp | Q30碱基比率 Q30 bases ratio/% | 读取总数 Total reads count | 碱基总数 Total bases count/bp | Q30碱基比率 Q30 bases ratio/% | |
| 根 | 54 938 740 | 8 240 811 000 | 90.00 | 47 635 582 | 6 586 283 353 | 94.92 |
| 茎 | 49 382 530 | 7 407 379 500 | 91.62 | 45 976 160 | 6 337 511 886 | 95.03 |
| 叶 | 48 584 930 | 7 287 739 500 | 91.73 | 45 100 670 | 6 203 066 849 | 95.05 |
2.2 五指毛桃unigene的注释分析
将获得的unigene用NCBI Blastx后比对到蛋白质数据库CDD、KOG、Nr、Swiss-Prot和KEGG等9个数据库,根据比对蛋白数据库的序列同源性结果可知（表3）,至少在1个数据库注释成功的unigene有62142条,占59.56%,其中与GO数据库匹配的unigene数目最多为51552条,占总数的49.41%;其次为Swiss-Prot数据库（49 172个）,占47.13%;在所比对数据库中都有注释的unigene数目为2397条,占2.21%。为了研究五指毛桃转录本序列与相近物种的近似情况和同源序列的功能信息,将序列与Nr数据库进行比对,分析unigene数据的物种分布,得到各序列之间的最佳匹配结果如图2所示。不同物种在不同水平上与五指毛桃unigene序列存在相似,其中与川桑（Morus notabilis）表现出明显的相似性,达到16.23%;其次是小麦（Hordeum vulgare subsp. Vulgare）和枣（Ziziphus jujuba）,相似性分别为2.67%和2.44%,与软克里藻（Klebsormidium flaccidum）和胶球藻C-169（Coccomyxa subellipsoidea C-169）相似性分别为1.53%和0.93%。
表3 五指毛桃unigenes在公共数据库中的注释数量和比例Tab. 3 Percentage and numbers of F. hirta Vahl annotated in public databases |
| 注释的数据库 Annotated in database | 注释到的转录本数量 Number of genes | 比例 Percentage/% |
|---|---|---|
| 保守结构域数据库（CDD） | 41548 | 39.82 |
| 真核生物蛋白质同源簇数据库（KOG） | 32530 | 31.18 |
| 美国国立生物技术信息中心非冗余蛋白库（Nr） | 46498 | 44.57 |
| 美国国立生物技术信息中心核苷酸序列数据库（Nt） | 36979 | 35.44 |
| 蛋白质家族数据库（Pfam） | 29100 | 27.89 |
| 瑞士蛋白质序列及注解数据（Swiss-Prot） | 49172 | 47.13 |
| 欧洲分子生物学实验室蛋白质数据库（TrEMBL） | 46164 | 44.25 |
| 基因本体数据库（GO） | 51552 | 49.41 |
| 京都基因与基因组百科全书（KEGG） | 3618 | 3.47 |
| 以上数据库至少一个注释成功 | 62142 | 59.56 |
| 以上数据库均注释成功 | 2307 | 2.21 |
对五指毛桃转录本的潜在功能进行预测和分类,将所有转录本与KOG数据库进行比对。结果如图3所示,35661条转录本被分成26个KOG群,其中翻译-核糖体结构和生物发生基因数目最多,为4840个,占4.64%;其次是翻译后修饰-蛋白质转换-分子伴侣为3669个,占3.52%;一般功能预测类有3633个,占3.48%;信号转导机制有3428个,占3.29%。同时,有1262个（1.21%）转录本参与次生代谢物的生物合成、运输和分解代谢。
图3 五指毛桃unigene序列KOG功能分类横轴上每个字母代表一个KOG的功能分类,见图注右侧;纵轴为注释到该分类下的基因数目。 Fig. 3 KOG functional classification of unigenes from F. hirta Vahl Each letter on the horizontal axis represents a functional classification of KOG, specific annotation comments are on the right figure. The vertical axis shows the number of genes annotated to this classification. |
基于GO分析,对所有unigene的潜在功能进行分类,如图4所示,共获得404727条GO注释信息,其中,在生物过程中,可将五指毛桃unigene分成26类,参与细胞过程的数量最多,为34056个,其次是代谢过程为30357个;在分子功能方面,五指毛桃的unigene被分成19类,参与结合和剪切活性的个数远远多于其他,分别为29343和27867个;在细胞组分上,五指毛桃的unigene可被分成22类,参与细胞的数目最多,为38314个。
图4 五指毛桃unigene序列GO功能分类横轴为GO的二级分类,纵轴为该分类基因个数(右)及其占被注释上基因总数的百分比(左)。不同颜色代表不同的GO分类。 Fig. 4 GO functional classification of unigenes from F. hirta Vahl The horizontal axis represents the secondary classification of GO, and the vertical axis indicates the number of genes in the classification (right) and their percentage of the total number of genes annotated (left). Different colors represent different orthologs |
利用KEGG数据库预测分析五指毛桃unigene可能参与的生物化学途径,结果如图5所示。总共有7744条unigene被注释到285个KEGG途径中,又可细分为细胞过程、环境信息处理、基因信息处理、代谢、生物系统,其中运输和分解代谢（289个）、信号转导（501个）、翻译（799个）、碳水化合物代谢（436个）、内分泌系统（187个）是每个分类中参与基因个数最多的。
图5 五指毛桃unigene序列KEGG功能分类横轴为KEGG二级分类,纵轴为该分类基因个数。不同颜色代表不同的分类。 Fig. 5 KEGG functional classification of unigenes from F. hirta Vahl The horizontal axis represents the secondary classification of KEGG, and the vertical axis indicates the number of genes in the classification. Different colors represent different orthologs |
2.3 五指毛桃根、茎、叶差异表达基因的分析
通过筛选,共获得27382个差异表达的基因,占总基因数的26.24%。用DESeq法对五指毛桃根、茎、叶3个样品进行两两比较,发现Leaf_vs_Root之间差异表达的基因个数最多为11755个,9550个基因上调表达,2205个基因下调表达;而Root_vs_Stem之间差异表达基因数目最少为5559个（图6）。
将差异表达的基因进行KEGG富集分析,结果发现在Leaf_vs_Root、Root_vs_Stem、Stem_vs_ Leaf分别有1363、624和1006个差异表达的基因在KEGG途径中显著富集,其中Leaf_vs_Root相比富集差异表达的基因个数较多的是核糖体（135个）、碳代谢（81个）、植物激素信号转导（56个）;Root_vs_Stem相比,核糖体（73个）、碳代谢（36个）是富集差异表达基因个数较多的通路;Stem_vs_Leaf相比,富集差异表达基因个数较多的是碳代谢（60个）、植物激素信号转导（53个）、淀粉和蔗糖代谢（40个）。分析各组分间具有显著差异的通路,发现Leaf_vs_Root、Root_vs_Stem、Stem_vs_Leaf分别有21、21和22条具有显著富集的通路,其中植物信号转导（plant hormone signal transduction）、植物MAPK信号通路（MAPK signaling pathway-plant）、黄酮类生物合成（flavonoid biosynthesis）为3组两两比较均共有显著富集的通路,而黄酮类生物合成代谢通路中差异基因的改变最为明显。
黄酮类生物合成通路是直接与次生代谢产物黄酮类化合物具有直接关系,分析后共获得与该代谢通路中相关基因个数为16个基因显著性差异表达,涉及到10种酶,分别为二氢黄酮醇4-还原酶（dihydroflavonol-4-reductase, DFR）、查尔酮合成酶（chalcone synthase, CHS）、反式肉桂酸4-单加氧酶（trans-cinnamate 4-monooxygenase, C4H）、查尔酮/查尔酮-黄烷酮异构酶（chalcone isomerase、Chalcone-flavanone isomerase, CHI）、莽草酸羟基肉桂酰转移酶（shikimate hydroxycinnamoyltransferase, HST）、咖啡酰辅酶A-O-甲基转移酶（caffeoyl-CoA O-methyltransferase, CCOMT）、黄酮类3°-单加氧酶（flavonoid 3°-monooxygenase, F3°H）、百花色素双加氧酶（leucoanthocyanidin dioxygenase, LDOX）、乙醇乙酰基转移酶（alcohol acetyltransferase, AAT）、细胞色素P450（cytochrome P450 family protein, CYP）,各类酶对应基因在五指毛桃根、茎、叶中的TPM值及相关详细信息见表4。通过对这些参与黄酮类生物合成途径中相关基因在五指毛桃根、茎、叶的TPM值的分析可知,8个基因在叶片中有最高的TPM值,7个基因在茎中有最高的TPM值,只有1个基因在根中有最高的TPM值。
表4 五指毛桃黄酮类生物合成通路中差异表达基因Tab. 4 Differently expression genes related to flavonoid biosynthesis of F. hirta Vahl |
| 酶 Enzyme | KO途径 KO pathway | 基因ID Gene ID | TPM值 TPM value | 长度 Length/bp | ||
|---|---|---|---|---|---|---|
| 根Root | 茎Stem | 叶Leaf | ||||
| 二氢黄酮醇4-还原酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN34248_c1_g1 | 13 | 167 | 4146 | 1369 |
| ko00941 Flavonoid biosynthesis | TRINITY_DN34248_c1_g4 | 30 | 219 | 363 | 1507 | |
| 查尔酮合成酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN33113_c1_g1 | 131 | 576 | 2721 | 1595 |
| ko00941 Flavonoid biosynthesis | TRINITY_DN34707_c1_g1 | 19 | 210 | 1187 | 1152 | |
| 查尔酮异构酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN33700_c3_g1 | 324 | 420 | 2427 | 1275 |
| 查尔酮-黄烷酮异构酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN31035_c3_g1 | 1 | 23 | 134 | 1071 |
| ko00941 Flavonoid biosynthesis | TRINITY_DN31035_c3_g2 | 40 | 161 | 23 | 972 | |
| 反式肉桂酸4-单加 氧酶 | ko00940 Phenylpropanoid biosynthesis | TRINITY_DN30694_c0_g4 | 6013 | 15421 | 12224 | 1952 |
| ko00940 Phenylpropanoid biosynthesis | TRINITY_DN26685_c0_g1 | 28 | 299 | 51 | 1879 | |
| 莽草酸羟基肉桂酰 转移酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN34559_c2_g3 | 1566 | 2822 | 221 | 1291 |
| 咖啡酰辅酶A O-甲基转移酶 | ko00945 Stilbenoid, diarylheptanoid and gingerol biosynthesis | TRINITY_DN29797_c0_g2 | 2046 | 1273 | 580 | 1081 |
| 黄酮类3’-单加氧酶 | ko00944 Flavone and flavonol biosynthesis | TRINITY_DN34838_c3_g2 | 450 | 244 | 883 | 1810 |
| 百花色素双加氧酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN27691_c0_g1 | 49 | 480 | 0 | 1292 |
| 乙醇乙酰基转移酶 | ko00941 Flavonoid biosynthesis | TRINITY_DN31328_c2_g2 | 125 | 1338 | 10 | 1575 |
| ko00941 Flavonoid biosynthesis | TRINITY_DN31152_c2_g1 | 2277 | 3452 | 7936 | 2147 | |
| 细胞色素P450 | ko00941 Flavonoid biosynthesis | TRINITY_DN33319_c0_g2 | 5575 | 10620 | 1494 | 1915 |
2.4 五指毛桃根、茎、叶基因转录本的SSR和SNP/InDel分析
通过五指毛桃根、茎、叶的基因转录本进行SSR分析,结果表明共有26 520条unigene有SSR位点分布。单碱基重复和双碱基重复是SSR的主要类型（图7）,分别占50.29%和31.46%,A/T和AT/AG分别为单核苷酸和二核苷酸的主要重复类型。
通过五指毛桃根、茎、叶的基因转录本SNP/InDel分析发现,其根、茎叶基因转录本SNP分别为47 051、60 922和70 272个,发生InDel事件分别为8665、10 046和10 123个,其中SNP的主要类型为C->T/G->A/A->G转换（图7）。
2.5 五指毛桃黄酮类化合物生物合成关键基因RT-qPCR分析
对黄酮类化合物生物合成中9个关键基因进行RT-qPCR,以根为参照样品,茎和叶中各基因的相对表达量见图8。这9个关键基因在根、茎和叶中的相对表达量与RNA-seq数据中TPM值的变化趋势基本一致,其中DFR、CHS、CHI、F3'H和AAT在五指毛桃的叶中有最高的相对表达量,C4H、HST和LDOX在茎中有最高的相对表达量,而CCOMT在根中有最高的相对表达量。
图8 五指毛桃黄酮类化合物生物合成关键基因RT-qPCR验证DFR：二氢黄酮醇4-还原酶;CHS：查尔酮合成酶;CHI：查尔酮异构酶;C4H：反式肉桂酸4-单加氧酶;HST：莽草酸羟基肉桂酰转移酶;CCOMT：咖啡酰辅酶A-O-甲基转移酶;F3°H：黄酮类3°-单加氧酶;LDOX：百花色素双加氧酶;AAT：乙醇乙酰基转移酶。 Fig. 8 RT-qPCR verification for key genes involved in flavonoid biosynthesis pathway in F. hirta Vahl DFR: dihydroflavonol-4-reductase; CHS: chalcone synthase; CHI: chalcone isomerase; C4H: trans-cinnamate 4-monooxygenase; HST: shikimate hydroxycinnamoyltransferase; CCOMT: caffeoyl-CoA O-methyltransferase; F3°H: flavonoid 3'-monooxygenase; LDOX: leucoanthocyanidin dioxygenase; AAT: alcohol acetyltransferase. |
3 讨论
本研究首次利用高通量测序技术对五指毛桃的根、茎和叶进行转录组测序,总共获得104335条unigenes,这些unigenes平均长度为578 bp,N50平均长度为1025 bp,Q20、Q30碱基百分比均大于94%,在GO、Siss-Prot、TrEMBL、Nr等4个数据库中均有大量unigenes被注释,说明本次转录组测序数据质量高,为五指毛桃有效成分生物合成及育种等相关研究提供理论基础。
黄酮类化合物是植物体内一类重要的次生代谢产物,大多与糖类结合形成糖苷。黄酮类具有多种多样的生物学功能,能调节果实风味、改变花果色泽等,还可作为一种重要的抗逆性物质,能够提高植物抗性[24,25,26]。同时,黄酮类化合物对人体健康还具有重要的生理功能,其可清除人体自由基,在抑菌、降血脂、消炎、抗病毒、抗癌等方面具有重要的作用[27,28]。药用植物黄酮类化合物含量与其种植环境和种植时间均具有密切的关系,对金线莲根、茎、叶中黄酮含量进行测定,发现叶的总黄酮含量占全草的总黄酮含量80%以上,而茎和根约各占不到10%,对组织培养、种植1个月和种植4个月的金线莲新鲜叶片的黄酮含量进行测定,发现黄酮类成分含量随着种植时间的延长而增加,并对黄酮类生物合成通路中的6个关键基因进行RT-qPCR验证,发现6个基因的相对表达量在种植4个月的金线莲叶片中最高[19];不同光质对植物中黄酮类化合物含量也有影响,蓝光能够促进大豆芽苗菜、青钱柳叶（Cyclocarya paliurus）中黄酮的合成[29,30]。
CHS是黄酮类化合物合成途径中的第1个关键酶,在该酶的作用下形成查尔酮,查尔酮能够通过进一步衍生形成各类黄酮类化合物[31],研究表明CHS的表达与否直接与植物黄酮类化合物含量多少有关[32]。CHI是黄酮类化合物代谢途径中的第2个关键酶,该酶是植物中查尔酮转化成黄烷酮必须的酶[33];刘长英等[34]克隆获得了桑树（Morus alba L.）查耳酮异构酶基因（MaCHI）全长序列,采用RT-qPCR法分析MaCHI在不同组织中的表达水平,发现该基因在叶片和果中的表达量较高,在根中表达量较低。F3°H属于细胞色素P450类,是黄酮类化合物代谢途径上的关键酶,能调控黄酮类化合物合成途径中的代谢流向,它催化生成二氢黄酮醇,而二氢黄酮醇是黄酮类化合物生物合成中生成黄酮、花色素和异黄酮等黄酮类化合物的重要中间产物[35]。二氢黄酮醇4-还原酶DFR是黄酮类生物合成途径中的下游关键酶,它能催化黄烷酮醇形成无色花色素,是花青素合成途径中的关键酶,无色花色素可作为花青素和缩合单宁的前体物质[32]。LDOX是一种2-酮戊二酸依赖性酶,通过抑制LDOX基因发现石竹目中花青素含量较少[36]。通过RNA-seq数据和RT-qPCR结果可知,黄酮类化合物合成途径中关键基因主要在五指毛桃叶和茎中高表达,因此,这几个基因的表达情况与黄酮类化合物的含有具有直接的关系。
本研究以2月龄的五指毛桃根、茎、叶为材料进行转录组测序,3个样品两两比较后发现存在大量差异表达的基因,对差异表达基因进行KEGG通路富集分析,结果发现主要涉及植物信号转导、植物MAPK信号通路、黄酮类生物合成这3条通路。对黄酮类生物合成通路中大部分关键酶基因的TPM值进行分析,发现大部分基因在叶片和茎中高表达,而仅有1个基因在根中高表达,同时,对其中9个关键基因进行RT-qPCR验证,发现这9个基因在五指毛桃根、茎和叶中的相对表达量变化趋势与RNA-seq数据一致,大部分基因都在茎和叶中高表达,因此,我们推测幼苗五指毛桃的黄酮类化合物主要积累在茎和叶片中,这与民间习惯用药部位根茎具有一定的差异,可能原因是本次转录组所用的材料为幼苗,而民间药用一般为种植多年的植株。五指毛桃根、茎和叶中黄酮类等次生代谢产物随着种植时间的延长,由叶片往根茎运输,并在根中大量积累,该研究成果为扩大用药资源具有重要的意义。因此,黄酮类生物合成的关键基因在植物生长过程中具有复杂的调控作用,其网络调控有待于进一步的研究。