致病性与无致病性茄科罗尔斯通氏菌转录组测序与分析

为了挖掘茄科罗尔斯通氏菌（Ralstonia solanacearum）致病相关基因,分别提取致病性（RS+）和无致病性（RS-）茄科罗尔斯通氏菌的总RNA进行转录组测序。两组测序数据经组装后得到5 111条Contigs,RS-和RS+对比发现779条差异表达基因（DEGs）,其中上调表达376条,下调表达403条。功能注释显示,这些DEGs主要集中在氨基酸代谢、膜转运、细胞运动和翻译等途径。信号通路富集显示,茄科罗尔斯通氏菌对寄主的响应由众多信号通路共同参与,其作用方式可能包括增强细菌趋化性、强化鞭毛组装能力、提高细菌分泌系统、上调双组分信号转导系统、调整MAPK信号通路和增强氨基酸代谢等。转录组测序还发现许多毒力基因,其中茄科罗尔斯通氏菌Ⅲ型分泌系统的hrp基因簇中的23个基因,有16个在无致病性茄科罗尔斯通氏菌中下调表达。随机挑选4个DEGs进行qPCR分析,荧光定量结果与测序数据基本一致,说明基于转录组分析DEGs的表达较为可靠。研究为进一步深入探究茄科罗尔斯通氏菌致病分子机理提供参考。     

基于转录组学和代谢组学联合分析草珊瑚萜类化合物生物合成的组织特异性分布

本研究旨在探究草珊瑚叶和根中萜类化合物的组织特异性分布差异,解析其药效品质差异形成的分子机制。采用液相色谱-质谱联用技术(liquid chromatography-mass spectrometry, LC-MS)和Illumina HiSeq^TM高通量测序技术获得草珊瑚[Sarcandra glabra (thunb) naka]叶和根的代谢组学和转录组学数据。代谢组学结果表明参与叶和根中萜类合成的差异代谢物有50个,包括法尼西基半胱氨酸、甘油醛-3-磷酸、甲羟戊酸-5-磷酸等。转录学结果表明差异代谢酶基因有57条,包括ACTC、HMGCR、MVK、DXS与KS等,并预测了MYB、C2H2、AP2/ERF-ERF等7个转录因子参与调控草珊瑚不同组织部位中萜类的合成和积累差异。实时荧光定量PCR (quantitative real-time polymerase chain reaction, qRT-PCR)结果显示,随机选取的8个参与萜类合成的酶基因在草珊瑚不同组织部位中的表达趋势与转录组学测序结果一致。本研究有助于阐明草珊瑚叶和根临床疗效差异形成的分子机制,同时为草珊瑚的资源开发利用奠定基础。     

岷江百合WRKY转录因子基因的克隆与功能分析

WRKY转录因子家族在植物应对非生物和生物胁迫的防卫反应中起重要作用。岷江百合(Lilium regale Wilson)是高抗枯萎病的野生百合,该研究基于前期转录组测序分析,采用RT PCR方法从岷江百合中克隆得到WRKY转录因子基因LrWRKY4并分析其功能,以探讨岷江百合对枯萎病菌尖孢镰刀菌(Fusarium oxysporum)侵染的转录调控机制,为进一步研究岷江百合WRKY基因家族的功能奠定基础。结果显示：(1)LrWRKY4开放阅读框为993 bp,编码330个氨基酸,LrWRKY4含有一个高度保守的‘WRKYGQK’七肽序列和一个C₂H₂锌指基序,属于Ⅱc类WRKY转录因子。(2)成功构建GFP LrWRKY4融合载体并通过根癌农杆菌介导转化洋葱表皮细胞,激光共聚焦显微观察发现,GFP LrWRKY4融合蛋白表达的绿色荧光特异性地分布在洋葱表皮细胞的细胞核中。(3)成功构建了过表达载体pCAMBIA2300s LrWRKY4并转化烟草,得到11个T₂代转基因烟草株系,且转LrWRKY4基因烟草与野生型烟草(WT)在表型上无明显差异;尖孢镰刀菌接种根部和叶片的实验结果显示,转LrWRKY4基因烟草对尖孢镰刀菌的抗性较WT明显增强;qRT PCR分析显示,岷江百合LrWRKY4基因在11个转基因烟草株系中均有表达,且转基因烟草中JA/SA信号途径相关基因的表达上调,并诱导了部分病程蛋白相关基因以及抗氧化相关基因的表达上调。(4)岷江百合鳞片浸染LrWRKY4 RNAi载体后的腐烂程度和病变面积均远大于RNAi空载转化的鳞片;瞬时表达LrWRKY4 RNAi载体的岷江百合鳞片中LrWRKY4基因的表达水平较对照下降了约45.7%,接种尖孢镰刀菌72 h后表达水平下降了93.8%;瞬时表达LrWRKY4 RNAi后,一些JA/SA信号途径相关基因的表达水平明显下降。研究表明,岷江百合LrWRKY4编码一个定位于植物细胞核的Ⅱc类WRKY转录因子;LrWRKY4基因能够在转基因烟草中稳定表达,且过表达LrWRKY4基因提高了烟草对尖孢镰刀菌的抗性;但瞬时表达LrWRKY4 RNAi降低了岷江百合JA/SA信号途径相关基因的表达并增强了对尖孢镰刀菌的敏感性。推测LrWRKY4基因是岷江百合抗尖孢镰刀菌防卫反应中的正调控因子,可能通过参与JA/SA介导的信号传导途径,诱导防卫相关基因的表达从而调节其对尖孢镰刀菌的抗性。     

层出镰刀菌乌头酸酶家族的功能研究

【目的】层出镰刀菌是引起苜蓿根腐病的主要病原菌之一,探究层出镰刀菌中乌头酸酶家族蛋白的功能特性,为深入认识层出镰刀菌基础生理代谢的分子机制提供依据。【方法】利用hmmsearch工具,对真菌中含有乌头酸酶结构域的蛋白进行检索,并进行系统进化分析;通过实时荧光定量PCR及SWISS MODEL建模技术分别分析FpACO基因的表达模式与蛋白结构;利用同源重组双交换方法构建层出镰刀菌乌头酸酶基因敲除突变体;分析ΔFpACO3、ΔFpACO4-1、ΔFpACO4-2敲除突变体的生长、产孢、孢子形态、环境胁迫响应及致病力等表型变化;进一步测定敲除突变体中线粒体代谢相关生理生化指标的变化情况。【结果】FpACO4-1与FpACO4-2在产孢及孢子形态发生中发挥作用;FpACO3、FpACO4-1、FpACO4-2参与调控层出镰刀菌对细胞壁胁迫及金属离子胁迫的敏感性;FpACO3、FpACO4-1、FpACO4-2影响线粒体代谢,包括总乌头酸酶活性、ATP含量、过氧化氢含量及三羧酸循环关键基因表达等。【结论】乌头酸酶家族参与调控层出镰刀菌产孢、孢子形态分化、细胞壁胁迫及金属离子胁迫响应和线粒体代谢等过程。     

不同抗逆性菠萝品种的差异基因和差异代谢物分析

该研究对2个不同抗逆性的菠萝品种PZ2(抗逆性较差)和PZ3(抗逆性强)进行转录组和代谢组比较分析,以探索品种间抗逆性差异的分子生理机制,为菠萝抗性育种提供理论依据。(1)转录组结果显示,以PZ2为对照,在2个品种间共筛选到1 667个差异表达基因,上调和下调的差异表达基因分别为770个和897个;其中筛选到20个与抗逆性相关的差异表达基因,包括WRKY和MYB转录因子基因以及ASR3、SOD1、POD48、HSP20、GSTF1、GSTU17等基因。(2)代谢组分析表明,以PZ2为对照,在2个品种间共筛选到208个差异代谢物,含量上调和下调的差异代谢物分别为98个和110个,其中22种差异代谢物与抗逆性相关,包括氨基酸及其衍生物、脂类、黄酮类、糖类和糖苷等。(3)qRT PCR分析验证表明,所选的8个差异表达基因在菠萝品种PZ2与PZ3中的差异表达趋势与转录组测序中差异基因表达水平变化趋势基本一致。(4)关联分析表明,差异表达基因和差异代谢物共同富集的通路有类黄酮生物合成,苯丙素生物合成,莨宕烷、哌啶和吡啶丙酸盐生物合成,半胱氨酸与蛋氨酸代谢等4条。研究发现,抗逆性强的菠萝品种的非生物逆境反应相关基因、抗氧化酶基因、逆境响应结构基因和调控基因的表达水平高于抗逆性较差的品种,且黄酮类、氨基酸衍生物类及木脂素、脂质等代谢物的含量也高于抗逆性较差的品种。     

携带Paenibacillus polymyxa WLY78固氮基因簇(nif)的重组大肠杆菌Escherichia coli78-7转录组分析

[目的]来自Paenibacillus polymyxa WLY78的固氮基因簇（nifBHDKEfNXhesAnifV）可以转化入Escherichia coli中表达并使重组大肠杆菌合成有固氮活性的固氮酶。本文拟通过对重组大肠杆菌E.coli 78-7的转录组分析以提高其固氮能力。[方法]对固氮条件（无氧无NH₄⁺）和非固氮条件（空气和100 mmol/L NH₄⁺）培养的重组大肠杆菌E.coli 78-7进行转录组分析。[结果]nif基因在两种培养条件下显著表达,说明在重组大肠杆菌中可规避原菌中氧气和NH₄⁺对nif基因的负调控。对于固氮过程必需的非nif基因,如参与钼、硫、铁元素转运的mod、cys和feoAB,这些基因在两种培养条件下表达水平有差异。而参与铁硫簇合成的suf和isc基因簇在两条件下表达水平差异巨大。此外,参与氮代谢的基因在固氮条件下显著上调。[结论]重组大肠杆菌中与固氮相关的非nif基因在该菌的固氮过程中具有较大影响,本文对在异源宿主中调高固氮酶活性研究具有重要意义。     

miRNA响应印度梨形孢定殖调控大麦生长发育的作用机制

【背景】内生真菌印度梨形孢(Piriformospora indica)定殖植物可以显著促进植物生长发育。miRNA已被证实在植物体的生长发育中具有调控作用。【目的】揭示印度梨形孢定殖大麦促进大麦生长发育过程中miRNA对印度梨形孢定殖的响应及对大麦生长发育的调控作用。【方法】提取大麦总RNA,实施转录组测序并进行序列比对与数据挖掘;使用高效液相色谱检测大麦生长素等激素水平变化。【结果】印度梨形孢对大麦有显著促生作用;全转录组测序结果显示：印度梨形孢侵染3 d较空白对照有18个差异表达的miRNA,其中11个miRNA上调、7个miRNA下调;侵染7 d与空白对照相比24个差异表达的miRNA,其中11个miRNA上调、13个miRNA下调;侵染3 d与侵染7 d相比有3个miRNA上调、6个miRNA下调。GO功能富集分析与KEGG通路分析显示,差异表达miRNA的靶基因主要参与转录、细胞分裂、生长素信号的感知和转导、光合作用和激素刺激响应。靶基因所参与的途径与大麦生长发育密切相关,暗示miRNA对印度梨形孢定殖过程做出了积极响应。代谢产物分析表明miRNA参与的调控路径的代谢产物发生改变。【结论】本研究以miRNA为入手点,探究了miRNA对大麦生长发育的调控机制,为揭示印度梨形孢的促生机制提供了新的研究方向。     

转录组分析罗伯茨绿僵菌精胺合成酶MrSPS介导真菌血腔定殖功能

【背景】精胺在植物应对逆境胁迫、动物抵抗疲劳和衰老、真菌生长代谢等过程中发挥重要作用,但目前在昆虫病原真菌中的研究未见报道。【目的】在分子水平上探究罗伯茨绿僵菌精胺合成关键酶——精胺合成酶在昆虫血腔定殖中的作用机制。【方法】显微注射法测定Mrsps敲除株ΔMrsps的致病力变化,并观察血腔中ΔMrsps生长状态;收集ΔMrsps和野生型WT注射侵染30 h后的大蜡螟血淋巴进行转录组测序,分别与罗伯茨绿僵菌和大蜡螟参考基因组进行比对分析,并结合定量PCR进行验证。【结果】与WT和回补株ΔMrsps-cp相比较,ΔMrsps致病力显著下降,而且随着注射浓度的降低,ΔMrsps致病力下降越显著。侵染36 h后WT和ΔMrsps孢子都能正常萌发且开始以类酵母状态生长,60 h后,相较于WT,ΔMrsps的生长繁殖数量较少。转录组共检测到3 202个罗伯茨绿僵菌基因,其中1 769个基因在ΔMrsps中表达上调,922个基因表达下调;差异表达基因涉及碳水化合物代谢、运输、分解代谢、翻译和氨基酸代谢等多条途径;筛选出28个血腔致病相关基因全部在ΔMrsps中表达下调;定量PCR检测发现在整个血腔定殖阶段免疫逃避蛋白Mcl1基因和血腔定殖Colonization of hemocoel 1基因在WT和ΔMrsps-cp中的表达量高于ΔMrsps。共检测到13 249个大蜡螟基因,其中4 026个差异表达基因;KEGG注释分析显示大量差异表达基因富集到内分泌系统和免疫系统等途径;深入分析发现22个差异表达基因归属于Toll和Imd信号通路,其中18个基因在ΔMrsps侵染的大蜡螟中表达上调,表明ΔMrsps侵染大蜡螟过程中更易引起免疫系统的激活。【结论】揭示了Mrsps在罗伯茨绿僵菌血腔定殖阶段作用的分子机制,为进一步揭示精胺在真菌中的作用机理提供了理论基础。     

嗜水气单胞菌中转录因子AHA_1581对细菌生理功能调控机制的研究

[目的] LuxR家族转录因子能够抑制或刺激不同功能类型基因的表达,来维持细胞功能的稳定性。嗜水气单胞菌是水产养殖中重要的致病菌之一,目前对该菌中的LuxR家族转录因子功能的研究还较少。[方法] 本研究利用含有sacB标记的自杀载体pRE112和同源重组技术敲除LuxR家族转录因子AHA_1581基因。[结果] 生理表型测定结果发现,ΔAHA_1581的运动与胞外蛋白的酶活增强、生物被膜形成能力降低,且耐受低温、卡那霉素、庆大霉素胁迫,但是对K₂Cr₂O₇更加敏感。进一步对野生型Ah和ΔAHA_1581的定量蛋白质组学分析,共鉴定到2654个蛋白,其中59个蛋白下调表达,142个蛋白上调表达。生物信息学分析表明AHA_1581参与调控双组分调节系统、丙酮酸代谢、碳代谢、TCA循环等细菌重要生理过程,以及细菌耐药基因和毒力因子的差异表达。[结论] 了解AHA_1581基因在调控细菌毒力以及生物过程中所起的重要作用,对预防和控制嗜水气单胞菌引起疾病的发生和传播可能具有重要的科学意义。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

三峡库区消落带野生狗牙根对水淹适生性的转录组分析

为探究狗牙根对水淹环境的适生性机制,取长势一致的狗牙根分为未水淹对照组(A1)和30 cm沉水处理组(A2),A2水淹15 d后两组同时采样,提取样本总RNA作转录组测序。A1和A2测序数据经de novo组装分别获得128031和83363条高质量Unigene;A1和A2对比发现28844条差异表达基因(differentially expressed genes, DEGs),其中11858条DEGs上调表达,16986条DEGs下调表达;功能注释显示,这些DEGs主要集中在转录、翻译、碳水化合物代谢和环境适应等途径;信号通路富集显示,包括植物对病原体的防御、物质代谢、光合作用、形态发育和抗氧化等过程的103个信号通路与狗牙根对水淹生境的适生性有关;随机挑选4个DEGs作qRT-PCR分析,荧光定量结果与测序数据基本一致,说明基于转录组分析DEGs的表达情况较为可靠。该工作或可为进一步挖掘与狗牙根水淹适生性相关的调控基因和功能基因提供基础。     

紫外辐照恢复早期耐辐射奇球菌的转录组学分析（英文）

目的 耐辐射奇球菌是一种对紫外线、电离、干燥和化学试剂具有较强抗性的极端微生物。然而，该菌在紫外辐照后恢复早期的分子响应还不完全清楚。本文的目的是揭示耐辐射奇球菌在这一阶段的转录组响应。方法 本研究采用RNA-seq技术，测定了正常和紫外辐照培养条件下耐辐射奇球菌的转录组。为确定关键的差异表达基因及其调控关系，进行了功能富集分析。选取部分关键差异表达基因，进行实时定量PCR实验验证。利用以往研究中的转录组数据，寻找紫外辐照、电离辐射和干燥胁迫条件下公共的差异表达基因。构建了蛋白质-蛋白质相互作用网络；对蛋白质互作网络中的枢纽基因和主要模块进行了鉴定；对这些枢纽基因和模块进行了功能富集分析。结果 紫外辐照后的恢复早期，上调基因数量是下调基因数量的2倍以上，且多数与应激反应和DNA修复有关。恢复早期的修复途径主要有单链退火（SSA）途径（涉及基因：ddr A-D）、非同源端连接（NHEJ）途径（涉及基因：lig B、ppr A）和核苷酸切除修复（NER）途径（涉及基因：uvr A-C），前两种途径为同源重组（HR）做准备，而NER途径去除紫外线照射带来的嘧啶二聚体。通过比较紫外辐照、电离辐...     

NaCl胁迫下哈茨木霉ACCC32524的转录组和代谢组分析

【目的】通过分析NaCl胁迫下哈茨木霉(Trichoderma harzianum)ACCC32524转录组和代谢组数据,研究差异表达基因及次级代谢产物的变化情况,初步探索响应NaCl胁迫的分子机制。【方法】利用Illumina HiSeq XTen高通量测序平台完成0、0.4、0.6 mol/L NaCl浓度胁迫培养下哈茨木霉ACCC32524的转录组测序,GC-TOF-MS技术完成对0mol/L和0.6mol/LNaCl胁迫培养下的差异次级代谢产物检测,利用相关软件及数据库对差异表达基因(DEGs)和次级代谢产物的注释、筛选和分类,并进行RT-qPCR验证。【结果】本研究分别得到0.4 mol/L和0.6 mol/L NaCl胁迫下417和733条差异表达基因;GO富集分析显示,分别有318和582条差异表达基因注释到生物学过程、分子功能和细胞组分3个一级分类和40个二级分类;COG分类结果表明分别有232和414条转录本为20个类别,涉及差异表达基因最多的分别为氨基酸的转运和代谢、一般功能预测、碳水化合物的转运和代谢;KEGG代谢途径分析结果表明,分别有75和96条基因归到25个代谢通路中(P≤0.05),其中涉及差异基因最多的是氨基酸的生物合成和2-氧代羧酸代谢通路。从转录组数据中共筛选出与渗透调节、离子转运、活性氧清除等22个耐盐相关基因。0 mol/L和0.6 mol/L NaCl胁迫下的代谢组数据中共筛选出101个差异次级代谢产物,包括8种积累量上调和93种下调物质,其中36个得到定性,分属于糖类、有机酸和氨基酸等9个分类中。RT-qPCR验证挑选的差异表达基因的表达量变化,均与RNA-seq分析结果一致。【结论】NaCl胁迫下引起哈茨木霉ACCC32524基因及次级代谢产物发生明显变化,细胞代谢途径发生明显偏移,这些进程共同作用减少NaCl对细胞的毒害作用,为木霉菌的耐盐机理研究提供重要信息。     

乙烯诱导水仙成花相关代谢物和基因的筛选与分析

为了解乙烯诱导水仙(Narcissus tazetta var.chinensis)成花的生理和分子机制,利用代谢组和转录组测序技术,筛选乙烯诱导水仙成花的差异表达代谢物和基因.结果表明,乙烯处理的侧芽检测到12个差异表达代谢物(DEM),包括7个上调,5个下调,其中,(±)7-表茉莉酸、多巴胺、亚精胺可能与乙烯诱导水...     

基于比较转录组学分析NaCl胁迫影响德尔卑沙门氏菌的耐渗机制

【目的】研究高渗胁迫条件下德尔卑沙门氏菌(Salmonella enterica subsp. enterica Derby, S. Derby)的转录组调控机制,分析差异表达基因(differentially expressed genes, DEGs)表达水平,探究在高渗胁迫影响下德尔卑沙门氏菌耐渗反应的相关代谢通路。【方法】通过高渗胁迫诱导德尔卑沙门氏菌的耐渗性,提取菌株的总RNA,去除rRNA,构建cDNA文库。利用转录组测序技术及生物学信息技术分析相关DEGs,并通过实时荧光定量PCR (real-time fluorescence quantitative PCR, qRT-PCR)进行验证。【结果】胁迫组德尔卑沙门氏菌通过转录组测序结果发现有3 950个DEGs,其中具有显著上调的基因21个,显著下调基因38个。涉及到细胞膜蛋白、氨基酸的代谢等相关基因上调,协助德尔卑沙门氏菌在高渗环境中存活。与此同时,胁迫组德尔卑沙门氏菌的糖转运系统(sugar transport system, PTS)、糖酵解过程以及抗氧化性相关基因表达显著下调,这是由于高渗环境菌体需要在体内储存大量糖类等物质,从而降低了糖原的消耗,进而导致细胞外膜的脂多糖合成受到抑制,降低了高渗胁迫下德尔卑沙门氏菌细胞膜表面的O抗原的合成。【结论】高渗环境诱导后显著提高了德尔卑沙门氏菌的耐渗性,其中Na⁺/H⁺逆向转运蛋以及谷氨酸的代谢通路发挥着重要的作用,为进一步了解以及更好地控制其在食品中的污染提供了理论依据。     

暗褐网柄牛肝菌菌核形成相关基因分析

【背景】暗褐网柄牛肝菌(Phlebopus portentosus)是第一个能够人工栽培的食用牛肝菌,人工栽培过程中,不同菌株会形成数量不等的菌核。【目的】探明不同菌株产核差异机制。【方法】采集多菌核(JH1)、寡菌核(JH2)菌株的成熟菌核及无菌核(JH3)菌株培养相同时间的菌丝体进行转录组测序,分析差异表达基因对菌核形成的作用和功能。【结果】KEGG富集分析显示,JH2 vs. JH1互比,苯丙氨酸、酪氨酸和色氨酸的生物合成,精氨酸和脯氨酸代谢,半胱氨酸及蛋氨酸代谢显著富集;JH3 vs. JH1互比,乙醛酸和二羧酸代谢显著富集;JH3 vs. JH2互比,谷胱甘肽、乙醛酸和二羧酸代谢显著富集。菌核形成相关基因分析显示,从JH2 vs. JH1、JH3 vs. JH1和JH3 vs. JH2的差异表达基因中分别筛选到69、118和82条与信号转导、感知刺激、防御、碳水化合物活性酶等有关的基因,其中碳水化合物活性酶基因数量最多。三个比较组共有的碳水化合物活性酶基因在JH1中的表达量高于JH2、JH3,表明JH1更能充分利用底物营养以形成更多菌核。【结论】本研究从转录组水平初步分析了暗...     

枫香叶片变色期全长转录组测序及分析

枫香因其树形优美,入秋后叶色红艳或橙黄,极具观赏价值,是优良的景观生态树种。为了解枫香叶片变色及其次级代谢过程的遗传基础,该文以枫香5个叶片变色期叶片混合样品为材料,利用单分子实时测序技术(PacBio平台)对其进行全长转录组测序。结果表明:(1)全长转录组测序共获得41.04 Gb的高质量数据,从中鉴定出全长非嵌合序列563 180条,通过聚类和去冗余,获得27 269条高质量全长转录本。在27 269条全长转录本中预测到2 035条长链非编码RNA(lncRNA),并检测出14 892个简单重复序列(SSR)位点和1 856个转录因子。(2)基因注释结果表明,NR、GO、COG、KEGG 等8个数据库共注释了24 857条转录本,KEGG数据库共获得了124个条代谢途径,主要有核糖体、碳代谢、氨基酸生物合成等,在类黄酮和叶绿素代谢途径中分别有49和71个转录本参与。上述结果初步揭示了枫香叶片变色期转录组信息以及功能特性,为后续研究枫香叶片变色分子机制、色素代谢合成途径和调控、相关功能基因克隆以及叶色改良提供基础数据。