神农香菊自然居群遗传变异评价研究及核心种质筛选＊

目的 本文对神农香菊（Chrysanthemum indicum var. aromaticum）自然居群的遗传多样性、遗传结构及特征代谢产物开展研究，并以此为依据筛选核心种质，为神农香菊的资源保护和应用提供理论依据。方法 利用21对分别来源于基因组和转录组数据的简单重复序列引物（SSR）对我国神农架林区的151份神农香菊样本、41份疑似野菊样本和武汉市洪山区11份野菊对照样本进行居群遗传多样性和遗传结构分析。基于遗传信息，采用高效液相色谱技术（HPLC）对不同遗传分组的代表性纯合样本开展9种代谢物的差异分析，利用最小距离逐步取样策略（LDSS）筛选神农香菊的核心种质。结果 21对引物在所有研究样本中共扩增出202个等位基因，每个位点平均等位基因数为9.619个。神农架5个居群的平均有效基因数（N_e）为2.518，香农多样性指数（I）为1.004，多态性信息含量（PIC）为0.526，杂合度（H）为0.246。居群间具有较低的遗传分化（F_st<0.12）和较高的基因流水平（N_m>1）。遗传结构分析显示，神农架的所有样品可分为三个遗传组。异绿原酸C、蒙花苷、木犀草苷等物质在三个组的样本间有明显的含量差异。基于遗传信息，最终筛选出了来自神农架3个地方居群的45份核心种质。结论 本研究展示了神农香菊自然居群整体存在的中高水平的遗传多样性，同时神农香菊样本与邻（混）生的疑似野菊样本具有清晰的形态差异，但两者间可能存在基因流或遗传混杂。这些信息为研究神农香菊的起源进化和对其资源的开发应用提供了理论和数据支撑。     

干旱胁迫下忍冬全基因组DNA甲基化和转录组分析

目的 分析干旱胁迫下忍冬Lonicera japonica全基因组DNA甲基化水平变化、模式以及与基因表达的关系,为进一步探究忍冬响应干旱胁迫的分子机制奠定基础。方法 采用全基因组重亚硫酸盐甲基化测序法（whole-genome bisulfite sequencing,WGBS）,测定干旱胁迫下忍冬全基因组DNA甲基化水平,并联合转录组测序结果对差异甲基化基因相关的差异表达基因进行分析。结果 干旱胁迫下忍冬整体甲基化水平普遍上升;其中CG类型甲基化水平明显高于CHG或CHH甲基化类型;对不同干旱胁迫处理下的忍冬进行差异甲基化区域（differentially methylated regions,DMR）分析,发现1935、2158和1750个差异甲基化相关基因,基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）富集分析表明,显著富集的通路主要包括亚油酸代谢,氮代谢通路和次生代谢物的生物合成通路等。通过转录组测序技术,发现差异甲基化基因相关的差异表达基因主要富集于次生代谢产物合成相关通路,同时筛选出4个在干旱胁迫下基因表达量上升显著的忍冬MYB（Lonicera japonica MYB,LjMYB）LjMYB转录因子。结论 干旱胁迫下忍冬甲基化水平上升,推测DNA甲基化调控基因表达是干旱影响金银花有效成分生物合成的作用机制之一。     

不同生长年限巴戟天蒽醌类含量差异比较及转录组学挖掘蒽醌类生物合成相关基因

目的 研究不同生长年限巴戟天Morinda officinalis根部蒽醌类化合物的含量差异及生物合成相关基因。方法 采用HPLC测定1年生和6年生巴戟天根中蒽醌类化合物成分,并进行比较转录组分析,采用qRT-PCR验证基因表达结果。结果 HPLC检测结果显示在6年生根中蒽醌类物质显著增加,通过比较转录组分析,从1年生和6年生巴戟天的根中筛选到8 619条差异表达基因,KEGG分析表明差异基因主要富集在各种次生代谢物的生物合成、苯丙素生物合成、泛醌和其他萜类醌生物合成等通路中。此外,共鉴定到13条蒽醌生物合成途径相关的差异表达Unigenes,其中ICS、SK、CS、DHNA-CoA synthase均在6年生根中上调表达,而DXS、DXR、DAHPS、IDH均下调表达,DHQD/SDH既有上调表达的又有下调表达的Unigenes。进一步选择6个候选基因进行qRT-PCR验证,验证结果与转录组数据一致。结论 为鉴别参与巴戟天蒽醌生物合成的相关基因提供参考,为后期进一步深入研究蒽醌类生物合成的累积规律分子机制提供基础数据。     

党参不同组织化学成分分布及转录组学分析

目的 分析党参不同组织转录组特征，解析党参根中活性成分积累的遗传基础，为党参的高品质生产与种植提供理论依据。方法 选取盛花期党参根、茎、叶、花不同组织为实验材料，采用高效液相色谱法（HPLC）检测党参苷Ⅰ、党参炔苷、苍术内酯Ⅲ含量，利用RNA-Seq对不同组织进行转录组测序，通过基因本体（GO）、京都基因和基因组百科全书（KEGG）数据库富集分析筛选分析差异表达基因，探讨党参不同组织化学成分分布特征及转录谱特征。结果 党参根中党参多糖及党参苷Ⅰ含量显著高于其他各组织。党参根转录谱与茎、叶、花存在显著差异，差异基因表达聚类分析、GO、KEGG富集分析发现差异基因主要富集在黄酮类、苯丙烷类生物合成、蔗糖淀粉代谢、植物激素信号传导、植物病原菌互作、MAPK级联信号传递、三磷酸腺苷结合转运盒（ABC）转运蛋白等途径。苯丙烷类化合物生物合成相关基因在根、花中表达量显著上调，ABC转运蛋白则多在花中高表达，党参多糖合成关键酶基因蔗糖：蔗糖1-果糖基转移酶（1-SST）、果聚糖1-外切水解酶（1-Feh）及大量与次生代谢产物积累密切相关的胁迫响应基因等在党参根中显著上调表达。结论 党参根中活性成分含量与转录谱特征与茎、叶、花等组织存在显著差异，表现出组织特异性，同时胁迫响应相关基因及活性成分菊粉型果聚糖、苯丙烷类化合物合成相关基因在根中表达上调。     

菊花脑芳樟醇合酶基因CnTPS1的克隆与原核表达分析＊

目的 克隆菊花脑芳樟醇合酶基因CnTPS1的全长编码序列，利用原核系统表达融合蛋白，为进一步研究该基因在菊花脑萜类合成中的功能提供理论依据。方法 以菊花脑基因组数据为基础，设计特异性引物，PCR扩增CnTPS1的全长编码序列，利用生物信息学分析软件分析序列特征。利用qRT-PCR技术分析CnTPS1基因在不同组织中的基因表达量。构建原核表达载体，体外诱导目的蛋白表达。结果 CnTPS1编码序列全长1749bp，编码582个氨基酸；基因表达模式分析表明该基因在茎和管状花中表达量较高；原核表达系统能诱导出67.58kDa大小的蛋白。结论 首次从菊花脑中克隆得到一个芳樟醇合酶基因CnTPS1，运用生物信息学方法对其编码蛋白的理化性质、结构特征等进行了分析预测，分析了该基因的组织表达模式，并在原核表达系统中成功诱导表达出目标蛋白，这些结果将为菊花脑萜类合酶基因的功能以及萜类物质生物合成途径的解析提供理论依据。     

响应ABA的甘草bZIP转录因子的基因表达及调控研究

目的 研究响应脱落酸（ABA）的甘草碱性亮氨酸拉链（GubZIPs）转录因子的基因表达差异,克隆并分析甘草关键酶基因的启动子序列,为解析bZIP转录因子调控甘草活性成分生物合成提供参考。方法 采用实时荧光定量聚合酶链式反应（qRT-PCR）分析甘草毛状根中GubZIPs转录因子基因在外源ABA处理下的表达模式及在3年生甘草中的空间表达差异;利用TBtools软件从甘草基因组中提取关键酶基因的启动子片段并克隆,利用生物信息学分析启动子上潜在的关键调控元件。结果 甘草GubZIP1、GubZIP5、GubZIP33和GubZIP56基因在不同的ABA处理条件下表达模式不同,其中GubZIP1和GubZIP33的响应最为显著。分析GubZIPs基因在甘草不同组织的表达差异发现,GubZIP1和GubZIP5在根部高水平表达,GubZIP33和GubZIP56在叶中高水平表达。基于甘草基因组克隆获得甘草生物合成途径中关键酶查耳酮异构酶（CHI）、查耳酮合成酶（CHS）、β-香树脂醇合酶（β-AS）基因的启动子序列,分析结果显示,CHI基因的启动子区域有2个G-box和1个A-box顺式作用元件,CHS基因的启动子区域有2个ABRE和1个G-box顺式作用元件,β-AS基因的启动子区域有1个ABRE和1个G-box顺式作用元件。结论 甘草毛状根中GubZIP1、GubZIP5、GubZIP33和GubZIP56基因对ABA均有显著响应,采用50 mg·L^–1的ABA处理3 h是研究甘草中ABA相关通路较为适宜的方案。此外甘草的关键酶基因启动子上存在可与bZIP转录因子结合和响应ABA等激素的顺式作用元件;研究结果可为深入解析响应ABA的bZIP转录因子调控甘草活性成分生物合成的分子机制提供参考。     

基于转录组测序揭示光质对刺五加实生苗基因表达的调控作用

目的 对刺五加Acanthopanax senticosus实生苗进行转录组测序,分析其响应不同光质调控的分子机制。方法 以生长60 d的刺五加实生苗为材料,设置不同光质处理组（LED-1、LED-2）和对照组[高压钠灯（high pressure sodium,HPS）],应用Illumina HiSeq^TM对其叶片进行转录组测序,对测序后的结果进行基因功能注释、差异表达基因（differentially expressed genes,DEGs）筛选和共表达网络分析。结果 转录组测序共获得126 252条Unigene。其中92 681条被注释,3个处理组中筛选出DEGs 7664个。基因本体（gene ontology,GO）富集结果表明,DEGs在代谢过程、刺激响应、细胞结构、催化活性等功能中得到显著富集。京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）途径富集分析表明,DEGs显著富集在植物激素传导、糖代谢、苯丙烷类生物合成、类黄酮生物合成。适度光质胁迫可以促进刺五加实生苗中苯丙素类化合物生物合成途径关键酶基因的表达,是促进刺五加有效成分累积的基础。结论 通过高通量转录组测序,揭示了不同光质对刺五加实生苗基因表达的调控特征,可为深入研究刺五加药效成分的生物合成机制及栽培中的光环境设置提供科学依据。     

艾活性成分生物合成调控相关的WRKY转录因子筛选及表达分析＊

艾（Artemisia argyi）是常用的中草药和灸用原料植物，具有重要的药用和经济价值，然而其药用活性成分的合成调控研究鲜有报道。WRKY转录因子家族是植物中最大的转录因子家族之一，在植物生长发育、抗逆以及次生代谢合成调控等方面发挥重要作用。本研究基于转录组数据系统鉴定艾的WRKY家族转录因子，分析其蛋白理化性质、进化关系和保守基序等特性，并研究该家族基因在茉莉酸胁迫下的表达模式。实验结果表明，从转录组水平共鉴定到37个AarWRKY转录因子，其氨基酸数目为119-547个，分子量为13.4-59.8 kDa，理论等电点为4.86-10.03。根据拟南芥WRKY家族分类，34个AarWRKY蛋白归类在I、II、III三组，而3个蛋白未归类。AarWRKY蛋白均具有保守的WRKY结构域，且同一组蛋白具有相似的保守基序。基因表达谱分析以及利用qRT-PCR方法检测基因表达差异发现，外源茉莉酸甲酯处理后，6个AarWRKY基因的表达发生了显著差异，其中3个基因上调表达，3个基因下调表达。进一步蛋白互作网络预测发现这些差异表达的AarWRKY基因可能参与JA响应和萜类化合物的生物合成调控等途径。本研究为艾WRKY家族基因的功能研究及其在艾叶活性成分累积中的调控作用解析奠定基础。     

剑叶龙血树4CL基因的克隆及其功能预测

目的 对剑叶龙血树4CL家族基因进行克隆及功能预测,筛选出可能参与黄酮类生物合成相关的4CL基因,为深入阐述龙血竭的形成机制奠定基础。方法 在前期转录组测序数据的基础上,初步筛选4CL家族基因,并对其进行基因克隆及功能预测,筛选可能与黄酮类及木质素类化合物合成相关的4CL基因;以剑叶龙血树组培苗为研究对象,对其进行过氧化氢、茉莉酸、水杨酸等胁迫处理,实时荧光定量聚合酶链式反应（qRT-PCR）分析不同胁迫处理条件下4CL基因的相对表达变化规律。结果 在剑叶龙血树转录组数据库中共注释得到6个4CL基因并成功完成其基因克隆;其中3个基因（Dc4CL1、Dc4CL2和Dc4CL3）在龙血竭形成过程中呈现差异表达,Dc4CL1基因和Dc4CL2基因的表达模式与木质素类化合物生物合成的关键酶基因表达变化趋势一致,而Dc4CL3基因的表达模式与黄酮类化合物生物合成的关键酶基因表达变化趋势一致;系统发育分析结果显示,Dc4CL1和Dc4CL2属于Ⅰ类4CL基因,而Dc4CL3基因属于Ⅱ类4CLs基因;胁迫处理结果显示,Dc4CL均有差异性表达变化,其中Dc4CL3基因的相对表达量显著高于其他Dc4CL基因。结论 对龙血竭形成过程中起关键催化作用的4CL基因进行初步的筛选,其中Dc4CL1基因和Dc4CL2基因（Ⅰ类）可能主要参与木质素类化合物的生物合成,Dc4CL3基因（Ⅱ类）可能是龙血竭形成过程中黄酮类化合物生物合成的关键催化酶基因。     

基于高通量测序技术探讨补骨脂宁治疗肺癌的潜在机制

目的 探索补骨脂宁治疗肺癌的潜在分子机制。方法 利用不同浓度的补骨脂宁干预A549细胞,采用噻唑蓝（MTT）法、高通量测序技术、趋势分析、基因本体（GO）功能注释和京都基因与基因组百科全书（KEGG）通路富集分析的方法分别检测A549细胞的增殖能力,筛选补骨脂宁抑制A549细胞增殖的关键基因以及关键通路,并利用实时荧光定量聚合酶链式反应（Real-time PCR）对测序的结果进行验证,探讨补骨脂宁治疗肺癌的潜在机制。结果 补骨脂宁可明显抑制A549细胞的增殖,并调节细胞内4 364个基因的表达。通过趋势分析发现这些差异基因分别聚类在20个基因表达模式中,其中前4个基因表达模式具有统计学意义,且均为显著下调趋势,说明补骨脂宁主要通过抑制某些基因的表达,发挥抗增殖的作用;通过组间差异基因的比较,笔者发现随着补骨脂宁浓度的增加,下调基因也随之增加,细胞增殖趋势越加减弱,与趋势分析结果一致,因此,笔者着重分析高浓度补骨脂宁干预组的278个差异基因,GO分子功能结果显示,补骨脂宁主要改变细胞和代谢的过程,KEGG通路富集结果显示,补骨脂宁主要通过影响类固醇生物合成、脂肪酸代谢、不饱和脂肪酸的生物合成、酮体的合成与降解和类固醇激素的合成等通路发挥作用,因此选取显著富集于脂质代谢通路并且P<0.01的差异表达基因LSS,硬脂酰基辅酶A脱氢酶（SCD）,3-羟基-3-甲基戊二酰辅酶A合酶基因（HMGCS1）,血管生成素样蛋白4（ANGPTL4）进行验证,Real-time PCR显示补骨脂宁可以抑制LSS,SCD,HMGCS1的mRNA表达,升高ANGPTL4的mRNA表达,与测序结果基本一致。结论 补骨脂宁主要通过靶向脂质代谢途径相关的代谢因子,抑制A549细胞增殖,从而达到治疗肺癌的作用。     

神农香菊挥发油GC-MS指纹图谱

目的:建立了神农香菊挥发油GC-MS指纹图谱分析测定方法,为神农香菊药材的质量控制提供参考。方法:采用水蒸气蒸馏法提取神农香菊挥发油,用GC-MS联用技术对神农香菊进行指纹图谱测定,并采用国家药品监督管理局推荐的中药指纹图谱计算软件进行计算,建立神农香菊挥发油的共有指纹图谱。结果:神农香菊挥发性成分中含有23个特征性指标成分,初步建立了以此23个共有峰为特征指纹信息的GC-MS指纹图谱。结论:该方法具有较好的重复性、精密度和稳定性(RSD均5%),建立的神农香菊挥发油指纹图谱为其质量控制奠定了基础。     

神农香菊茎叶总黄酮与总酚的提取纯化和抗氧化活性研究＊

目的 提高神农香菊茎叶总黄酮与总酚的产率和纯度并探究其抗氧化机制。方法 采用正交试验设计，优化神农香菊茎叶总黄酮和总酚的提取参数。使用大孔树脂吸附技术，对提取物中总黄酮和总酚进行纯化。通过UPLC-MS/MS分析纯化前后提取物中黄酮和酚酸类成分的含量变化。采用1，1-二苯基-2-苦肼基（DPPH）清除试验和分子对接评价提取物的抗氧化潜力。结果 最佳提取工艺为：料液比1：5 g·mL^-1，乙醇浓度70%，提取时间60 min。最佳纯化条件为：上样浓度0.12 g·mL^-1，上样体积3 BV，上样流速3.0 BV/h，洗脱溶剂80%乙醇，洗脱流速4.5 BV/h，洗脱体积4 BV。与纯化前相比，纯化后提取物中总黄酮和总酚纯度分别提高了2.85倍和4.03倍。UPLC-MS/MS分析证实，纯化后部分黄酮和酚类物质含量显著增加。此外，神农香菊茎叶提取物显示出较强的DPPH自由基清除活性，其黄酮成分与黄嘌呤氧化酶具有较强结合能力。结论 本文对神农香菊茎叶总黄酮与总酚的富集和纯化具有重要指导意义，并表明了其潜在抗氧化价值。     

款冬叶不同发育阶段的转录组学分析

目的比较不同生长时期款冬叶中苯丙素类成分生物合成相关基因的表达,推测苯丙素类成分的最佳合成积累时期,为款冬叶资源的合理利用提供科学依据。方法利用Illumina Hi Seq2500测序技术对不同时期的款冬叶进行转录组测序,获得转录组数据后进行基因功能注释等生物信息学分析,对苯丙素类生物合成相关基因在不同时期的表达进行比较。结果通过转录组测序技术获得46 793个Unigenes,平均长度为952.144 8 bp,其中有4 774个可以被NR、Swiss-Prot、eggNOG、GO、KEGG等公共数据库共同注释。对注释得到的Unigenes进行KEGG代谢通路分析,结果显示款冬叶转录组中有144条Unigenes参与萜类、黄酮类和苯丙素类生物合成,其中萜类65条,苯丙素类64条,黄酮类15条。进一步对参与苯丙素类生物合成相关基因进行动态比较,发现其关键酶PAL、4CL、HCT、CCoAOMT等在9月的表达量最高,即苯丙素类成分在9月的生物合成量可能最高。结论为苯丙素类成分的生物合成途径及调控分析奠定了基础,也为款冬叶资源的开发利用奠定了科学依据。     

赶黄草的转录组测序及分析

目的获得赶黄草Penthorum chinense转录基因参考序列和表达量信息,探索赶黄草药用活性成分生物合成的遗传基础,为赶黄草的重要功能基因研究提供参考。方法采用Illumina Hi Seq 2000高通量测序技术,对赶黄草进行转录组测序,对获得的数据进行过滤组装,对得到的unigene进行比对注释,同时对本研究获得赶黄草活性成分合成相关代谢通路基因进行分析。结果共获得了40 005 442条有效的短读序,通过de novo拼接得到42 306个unigene,开放阅读框(ORF)分析共计得到518个ORF,其中75个ORF具有转录因子结构域。通过KEGG分析,检测到33、32、59、68个unigene分别参与黄酮类生物合成、固醇类生物合成、萜类骨架生物合成和羧酸代谢。结论本研究发现的这些基因对赶黄草遗传改良和药用活性物质的生物合成路径相关基因研究有着重要意义。     

利用转录组分析款冬萜类化合物生物合成关键酶基因及表达特征

目的挖掘与款冬萜类化合物生物合成途径相关的关键酶基因。方法利用Illumina HiSeq2500测序平台对野生款冬花蕾和叶进行转录组测序,使用Trinity软件de novo从头组装,并通过基因功能注释、差异表达基因分析等生物信息学方法进行分析,挖掘款冬中萜类化合物生物合成途径相关的酶基因。结果转录组测序获得39 912 371条高质量reads(SRA号:SRR9113366,SRR9113367),组装获得91 118条Unigene,55 830条Unigene在NR、Swiss-Prot、GO、COG、KEGG等数据库得到注释。鉴定出参与款冬萜类化合物生物合成相关酶基因的129个Unigene,包括91个参与三萜骨架生物合成酶基因,32个萜类合成酶基因,6个细胞色素P450基因,其中差异表达基因25个,涉及上游MVA途径的4个差异基因在花蕾中表达量高于叶,MEP途径的5个差异基因在叶中表达高于花蕾,另外还有10个基因在叶中高表达,9个基因在花蕾中高表达。根据差异基因HMGR、TPS、AS、CYP450基因在花蕾中高表达,推测可能与花蕾中萜类物质含量高有关。结论从款冬转录组数据库中初步获得了可能参与萜类化合物生物合成途径的候选关键酶基因,为进一步阐明款冬萜类化合物生物合成的分子机制奠定基础。     

乌头萜类生物合成代谢的转录组学分析

目的:对乌头的6个器官组织,主根、侧根、茎、叶、花、果实进行转录组测序分析,以研究参与其萜类次级代谢生物合成相关基因的表达谱,挖掘其功能基因。方法:采用Illumina Hi Seq 2500平台测序、组装得unigenes,与公共数据库NR,Swiss-Prot,GO,COG,KOG,KEGG比对、注释以获得差异基因,使用实时荧光定量聚合酶链式反应(Real-time PCR)验证参与萜类次级代谢生物合成的关键基因。结果:本研究共获得156 967 635条Clean Reads(28 254 174 300 bp),经序列合并拼接后获得103 337个unigenes,通过核苷酸和蛋白序列等方面的同源性分析,表明其中37.31%(38 554个unigenes)与其他生物的已知基因具有不同程度的同源性,通过功能分类研究和代谢途径分析(KEGG)获得参与萜类合成158个unigenes(编码5个关键酶)。结论:这些发现的基因将为乌头萜类化合物的生物合成途径及分子机制提供基础数据。     

鱼腥草乙酰辅酶A酰基转移酶基因克隆、表达及生物信息学分析

目的通过克隆鱼腥草乙酰辅酶A酰基转移酶(acetyl-Co A C-acetyltransferase,AACT)基因c DNA全长,并展开生物信息学分析和表达分析,为解析鱼腥草萜类次生代谢产物的生物合成机制奠定基础。方法采用RT-PCR方法获得AACT基因c DNA序列并进行生物信息学分析;利用实时荧光定量PCR方法检测了AACT基因在鱼腥草的地下茎、地上茎、叶、花中的表达情况。结果克隆获得的AACT基因c DNA全长为1 218 bp,编码405个氨基酸。生物信息学预测AACT蛋白不含跨膜区,不含信号肽。AACT基因在鱼腥草地上茎中的表达量最高,达到1.49,其次是地下茎0.96,花和叶中的表达量相对较低,分别为0.20和0.10。结论首次从鱼腥草中克隆了AACT基因,为进一步阐明该基因在鱼腥草萜类化合物代谢途径中的重要作用奠定基础。     

苗药八爪金龙转录组测序与次生代谢产物合成相关基因的挖掘

目的对八爪金龙根进行转录组测序,挖掘次生代谢合成相关基因,探索八爪金龙次生代谢产物生物合成的分子基础。方法采用IlluminaHi Seq4000高通量测序技术,对八爪金龙根进行转录组测序。使用Trinity软件对获得的Unigenes数据进行过滤组装,运用KEGG数据库对Unigenes进行注释。结果共获得52249条Unigenes,其中31391条被公共数据库成功注释,1507条Unigenes被注释到次生代谢产物合成。通过对转录组数据深入挖掘发现,参与八爪金龙苯丙素生物合成的Unigenes共有126条,参与萜类化合物骨架生物合成的Unigenes数量为73条,参与黄酮类化合物生物合成Unigenes有58条,参与次生代谢后修饰的Unigenes共有253条。筛选出9条Unigenes编码4个与八爪金龙香豆素生物合成的关键酶,23条Unigenes编码8个与黄酮生物合成的关键酶,39条Unigenes编码19个与萜类生物合成的关键酶,140条CYP450基因和113条UGT基因可能参与次生代谢物的修饰。微卫星识别软件(microsatellite identification tool,MISA)分析发现八爪金龙转录组包含17 400个简单重复序列(simple sequence repeats,SSR)。结论通过高通量转录组测序,初步揭示了参与八爪金龙次生代谢产物合成相关的基因,为进一步研究八爪金龙次生代谢产物合成途径关键酶的功能及其调控机制奠定了基础。