摘要[目的]为获得刀豆转录组信息特征,丰富其基因数据库,并为后续遗传图谱构建、品种选育以及解析刀豆黄酮类化合物次生代谢途径和调控机制研究提供基础数据.[方法]基于Illumina 6000平台对刀豆各组织部位进行转录组测序,组装所有转录本并进行去冗余.对获得的去冗余基因(unigenes)在NR、Swiss-Prot、KEGG、KOG、GO五大数据库中进行功能注释.随后进行基因表达量分析,并对刀豆各组织部位及3个不同时期的果荚进行基因差异表达分析.对差异表达基因进行京都基因与基因组百科全书(KEGG)富集分析,并对类黄酮生物合成途径相关基因进行注释和分类.[结果](1)测序组装后共获得137 278条去冗余的unigenes.(2)在五大数据库中总计获得功能注释76 793条(占unigenes总量的55.95％),具体为:KOG数据库40 965条(占注释总量的29.84％)、KEGG数据库32 776条(占注释总量的23.88％)、NR数据库67 373条(占注释总量的49.08％)、Swiss-Prot数据库56 397条(占注释总量的41.08％)、GO数据库31 399条(占注释总量的22.87％).(3)通过比较3个不同时期的刀豆果荚,分别鉴定到1015、6 445和7207条差异表达基因.(4)KEGG富集分析表明,差异表达基因显著富集于植物激素信号转导、苯丙烷代谢、类黄酮生物合成等通路.(5)在类黄酮生物合成途径中共注释到57个基因,这些基因可能在刀豆黄酮类化合物生物合成中发挥重要作用.[结论]成功获取刀豆的转录组数据,组装获得大量unigenes并进行功能注释,揭示不同组织部位及果荚发育时期的基因表达差异及其富集通路,特别是发现57个与类黄酮生物合成相关的潜在关键基因.这些结果显著丰富刀豆基因数据库,为刀豆的遗传图谱构建、分子育种以及深入解析其黄酮类化合物次生代谢途径和调控机制奠定数据基础.