海马齿根系响应盐胁迫的转录组分析

【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用Illumina测序技术对0 mmol/L NaCl （对照组）和400 mmol/L NaCl胁迫处理（盐胁迫处理组）下的海马齿根系进行转录组测序分析,从中筛选出差异表达基因,选取13个基因进行实时荧光定量PCR （qRTPCR）检测,以验证转录组数据的可靠性。【结果】在海马齿根系转录组中共鉴定出305145个转录本,平均长度为622 bp,其中,对照组有146177个长度>300 bp的转录本,盐胁迫处理组有72173个长度>300 bp的转录本;共有65535条Unigenes在Nr、GO、Swiss-Prot、COG和KEGG五大数据库注释成功,占Unigenes总数的52.36%。对照组和盐胁迫处理组共有65535个差异Unigenes,其中,有182个热休克蛋白基因。对照组和盐胁迫处理组间共有24042个差异表达基因,从中选取13个基因进行qRT-PCR检测,结果显示,9个基因表达上调,其余4个基因表达下调,与转录组测序结果一致。24042个差异表达基因中,共有10106个显著差异基因富集到129条代谢通路,其中富集程度排名前10的代谢途径为核糖体、次级代谢生物合成、RNA转运、内吞作用、剪接体、甘油磷脂代谢、内质网加工、吞噬、醚脂类代谢和植物-病原体相互作用,参与盐胁迫相关的硫代谢、脯氨酸积累、活性氧（ROS）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。     

截形叶螨雌成螨应对高温胁迫的转录组分析

为明确截形叶螨Tetranychus truncatus雌成螨应对热胁迫的差异表达基因，采用Illumina HiSeq 4000进行转录组测序。结果表明，截形叶螨雌成螨热胁迫后共检测到17 577个差异基因，其中12 831个上调，4 746个下调。GO功能富集发现，雌成螨热胁迫后的差异基因主要富集于细胞过程、催化活性和细胞，KEGG通路富集分析发现，差异基因主要富集在代谢途径和溶酶体等。从转录组测序数据中筛选出与高温存在关联的抗氧化酶和热激蛋白等基因。采用RT-qPCR技术检测这些基因在截形叶螨雌成螨上的表达特征，发现17条基因的上下调趋势与转录组测序结果一致。这为后续深入研究截形叶螨与热胁迫相关基因的功能奠定了基础。     

基于转录组学的巨菌草内生链霉菌拮抗玉蜀黍平脐蠕孢的分子机理

通过Illumina Hiseq4000高通量测序平台,采用PE150测序策略研究了玉蜀黍平脐蠕孢在链霉菌发酵液胁迫下(S组)的转录组变化,共得到128 613 786个clean reads,筛选出4 559个差异基因.与无发酵液的对照组相比,S组中有2 283个基因表达上调,2 276个基因表达下调;从中挑选出表达量差异最大的22个基因作为内生链霉菌调控拮抗的拟关键基因,这些基因多涉及氨基酸合成和代谢以及信号转导途径,说明链霉菌发酵液对玉蜀黍平脐蠕孢的氨基酸生物合成和MAPK信号通路相关基因的影响显著;采用实时荧光定量PCR验证其中5个基因的表达量,其趋势与转录组分析结果基本一致.     

生物炭添加下枯草芽孢杆菌应对镉胁迫的分子机制

为了探究生物炭添加前后枯草芽孢杆菌应对镉胁迫的分子机制，本研究基于转录组测序，对比分析了枯草芽孢杆菌在4种体系（正常体系、镉胁迫体系、生物炭体系、生物炭-镉胁迫体系）中的基因表达差异。结果表明：在镉胁迫压力下，枯草芽孢杆菌共产生1 932个差异表达基因，其中上调基因967个，下调基因965个。枯草芽孢杆菌的代谢活性受到明显抑制，但是其鞭毛组装和细菌趋化性通路明显增强，离子转运和阳离子跨膜转运蛋白功能条目的基因表达也呈现上调，说明枯草芽孢杆菌通过上调细胞运动及离子转运相关基因的表达来规避镉毒性压力。相比于镉胁迫体系，在生物炭-镉胁迫体系中，细菌共产生1 541个差异表达基因，其中上调基因691个，下调基因850个。生物炭的存在增强了枯草芽孢杆菌生物过程、细胞组分和分子功能基因的表达，枯草芽孢杆菌物质和能源代谢相关通路的基因表达均呈现上调，并且细菌体内镉抗性基因czcD、cadA和膜功能相关差异基因的表达也以上调为主，说明生物炭可以缓解镉胁迫对枯草芽孢杆菌的毒性效应，并可刺激枯草芽孢杆菌的金属抗性能力。     

低温胁迫下枇杷幼果转录组的De novo组装和功能注释

【目的】分析低温胁迫下枇杷幼果的转录组,以探讨幼果冻害的分子基础,为深入鉴定枇杷果实生长发育、抗寒性基因和抗寒蛋白提供资料。【方法】以花后75d的"大五星"枇杷幼果为试材,研究低温胁迫(3℃条件下处理12h)下幼果转录组的De novo组装和功能注释,并对其进行生物学信息分析,探索基因表达差异。【结果】低温胁迫下枇杷幼果通过De novo组装产生116 723条contigs,并应用Illumina高能量测序技术得到64 814条unigenes,其中有42 830条被nr、KEGG、COG、GO注释。此外,组装和注释了4 017个枇杷幼果在低温胁迫下表达的差异基因。【结论】发现在枇杷低温耐性中起重要作用的油菜素甾醇生物合成和磷脂酰肌醇信号系统的基因表达上调。     

盐胁迫培养下季也蒙毕赤酵母的转录组学差异分析

在本实验室从海南省西沙群岛海沙里分离出1株季也蒙毕赤酵母(Meyerozyma guilliermondii),并证实该酵母菌株有较好的耐盐性,在12%NaCl胁迫培养条件下仍可生长的基础上,本文通过Illumina HiSeqTM高通量测序技术对经盐胁迫培养与非盐胁迫培养24 h的M. guilliermondii进行转录组测序比较。结果表明:2组样品共有1 027个显著性差异表达基因,其中458个基因表达上调,569个基因表达下调。通过基因本体(gene ontology, GO)功能注释发现,经盐胁迫处理后M. guilliermondii的差异表达基因主要富集在生物过程分类中,其中核苷酸代谢、糖代谢及辅酶代谢基因产生较大差异。对差异基因进行京都基因与基因组百科全书数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)通路富集分析,发现大部分富集通路都与细胞分裂和代谢有关,与GO富集结果相对应。上述结果可以为盐胁迫培养对M. guilliermondii的生长影响及进一步的生物学探究提供科学依据。     

镉胁迫下美洲黑杨无性系‘中菏1号’转录组分析

  目的  通过转录组测序和分析,发掘美洲黑杨无性系‘中菏1号’耐镉关键功能基因,揭示美洲黑杨在镉胁迫下的响应机制。  方法  以‘中菏1号’为试验材料,通过沙培试验法对处在不同质量浓度（0 mg/L（ZH1C0）,10 mg/L（ZH1C1）,20 mg/L（ZH1C2））Cd2+ 环境下的无性系进行胁迫处理,而后采集叶片,提取RNA,并使用Illumina HiSeqTM2500系统对转录组进行高通量测序。  结果  主成分分析确定PC1作为区分不同Cd2+ 浓度样本的标准;经过差异表达分析,在ZH1C0和ZH1C2之间共检测到4 109条差异基因,其中2 741条基因显著上调表达,1 368条基因显著下调表达;在ZH1C1和ZH1C2之间检测到3 710条差异基因,其中1 969条基因显著上调表达,1 741条基因显著下调表达;不同数据库之间的注释表明Cd2+ 不仅会影响叶绿体、线粒体等细胞器的微结构,而且会干扰光合作用和呼吸作用中的电子传递。转录因子分析发现,WRKY、MYB、ZIP、ERF、bHLH在镉胁迫响应和耐受中具有重要作用。  结论  在‘中菏1号’中,异戊二烯化植物蛋白在阻止Cd2+ 进入细胞内部发挥了重要作用;谷胱甘肽不仅可以促进细胞内Cd2+ 的螯合,还可以缓解细胞内的氧化胁迫;ABC家族蛋白和MATE家族蛋白是‘中菏1号’体内转移Cd2+ 螯合蛋白的关键载体。      

盐胁迫下转 DcCIPK24 拟南芥的转录组比较分析

【 目 的】 探 究 铁 皮 石 斛（Dendrobium catenatum Lindl.）DcCIPK24 的 下 游 响 应 基 因， 为 研 究DcCIPK24 提高耐盐性的分子调控途径提供指导。【方法】利用 Illumina NovaSeq 6000 测序平台对 200 mmol/L NaCl 处理和对照组的转 DcCIPK24 拟南芥、野生型拟南芥进行测序，分析差异表达基因的富集通路，并对通路中的差异基因进行荧光定量 PCR 验证。【结果】盐胁迫下，从 2 个转基因拟南芥株系与野生型拟南芥的比较组合中共鉴定出 96 个差异表达基因；GO 注释分析发现 96 个差异基因主要被富集在分子功能相关通路；KEGG 富集分析发现差异基因主要被注释在氨基糖和核苷酸糖代谢、植物 MAPK 信号通路和甘油酯代谢通路中，选取上调表达的差异基因 AtGPAT5、AtSIRK、AtMEE25、AtUGE5 和下调表达基因 AtAMT1-2、AtATTI3、AtCYP71A25、AtLHCB4.3 进行实时荧光定量 PCR 验证，结果表明盐胁迫下 8 个差异基因表达趋势与转录组数据基本一致。【结论】盐胁迫下，DcCIPK24 主要通过氨基糖和核苷酸糖代谢途径、植物 MAPK 信号途径和甘油酯代谢通路响应耐盐。     

外源茉莉酸对燕麦在不同干旱胁迫下转录组的影响

为探究外源茉莉酸(Jasmonic acid,JA)对不同干旱胁迫下燕麦转录组的影响,本研究测定轻度(0.017 mol/L PEG-6000处理)和重度干旱胁迫(0.034 mol/L PEG-6000处理)处理及再添加JA处理的‘蒙燕1号’转录组,并对表达显著差异基因(DEGs)进行GO和KEGG功能富集分析。结果表明:GO功能富集分析显示,不同干旱胁迫处理下,外源JA诱导的DEGs均主要富集在膜、催化酶、转移酶和对刺激反应等代谢途径。KEGG功能富集分析显示,在未进行干旱胁迫条件下(CK),外源JA诱导1 955个基因表达上调,4 666个下调;DEGs主要富集的通路为代谢途径和次生代谢物的合成,其中苯丙素的生物合成通路基因表达发生显著变化。轻度干旱胁迫下,外源施用JA诱导1 574个基因表达上调,933个基因下调;DEGs主要富集的通路为次生代谢的生物合成、植物激素信号转导、精氨酸和脯氨酸代谢通路,其中脱落酸和水杨酸信号转导、精氨酸和脯氨酸代谢通路中基因表达均上调;重度干旱胁迫下,外源JA诱导223个基因表达上调,456个下调;DEGs显著富集的只有鞘脂代谢通路。与未施用JA相比,不同干旱胁迫下,施用JA均可诱导细胞分裂素(CTK)代谢通路中CTK的N-糖基化基因(UGT73C5)上调和氧化酶/脱氢酶基因(CKX11)下调,以及脱落酸(ABA)信号转导相关基因(THI1.3)和未知功能基因(AT1G71250)表达发生显著变化。综上,轻度和重度干旱胁迫下,外源JA能诱导燕麦响应干旱胁迫的CTK和ABA相关基因表达发生显著变化。     

荞麦根的转录组学分析及黄酮合成基因的鉴定

荞麦基因资源相对匮乏。以甜荞根和金荞根为材料,利用Illumina Hi SeqTM2000对其进行转录组测序,经过de novo从头组装得到64 618条Unigene,其中37 546条基因得到了有效注释。对甜荞根和金荞根中的差异基因进行筛选,共获得了4 709个差异基因,其中2 460个上调表达,2 249个下调表达。这些差异基因的GO注释集中在10个细胞组分(cellular component)、12个分子功能(molecular function)和22个生物学过程(biological process)中。对差异基因参与的代谢途径进行富集,上调基因中共有40个代谢途径显著富集,下调基因中共有18个代谢途径显著富集,这些显著富集的代谢途径参与甜荞根和金荞根的生物学调控。最后对转录组中注释到黄酮合成途径中的关键基因进行了分析,共有21个基因注释到黄酮合成途径中的关键酶。不仅丰富了荞麦的基因资源,为荞麦分子生物学研究提供数据基础,也为筛选甜荞根和金荞根中相关基因的表达趋势提供参考依据。     

促生解淀粉芽胞杆菌SQR9定殖对黄瓜根系基因表达的转录组学研究

[目的]本文旨在阐明黄瓜植株接种促生解淀粉芽胞杆菌SQR9后其根系基因的表达谱变化特征,从植物响应微生物的角度揭示植物根际促生细菌(plant growth-promoting rhizobacteria, PGPR)的作用机制。[方法]利用Illumina高通量测序技术研究接种菌株SQR9对黄瓜根系基因表达谱影响的转录组特征,对显著差异表达基因进行GO功能和KEGG富集分析。[结果]与未接种的对照黄瓜植株比较,菌株SQR9接种72 h后黄瓜根系有484个基因的表达发生显著变化,包括300个上调基因和184个下调基因。Real-time PCR验证表明基因表达差异趋势与转录组测序结果一致,证明测序结果可信。黄瓜根系响应菌株SQR9的差异基因主要参与碳水化合物代谢、氨基酸代谢、次级化合物代谢和信号转导途径等。在植物促生方面,氨基酸代谢途径中ASP1基因(编码天冬氨酸转氨酶)表达上调约4倍;糖酵解途径中PDC基因(编码丙酮酸脱羧酶)和卡尔文循环中FBA6基因(编码果糖二磷酸醛缩酶)表达均上调3～5倍;生长素响应途径中编码生长素抑制子的IAA8基因表达显著下调,而编码响应蛋白的基因SAUR77和SAUR21表达均上调3倍左右,表明菌株SQR9产生的外源吲哚乙酸(indole-3-acetic acid, IAA)激活了植物内源的IAA信号途径。在植物免疫方面,免疫信号途径转录因子的编码基因WRKY29和ERF1B表达上调,激活植物系统抗性并有助于提高抗病能力;编码抗胁迫因子的C_4H和ADH1基因表达也显著上调,提高植物抗逆能力。[结论]菌株SQR9可通过调控黄瓜根系代谢相关基因表达和激活免疫途径相关基因的方式促进黄瓜植株生长并提高其抗病能力。     

高温胁迫下半夏叶片的转录组分析

【目的】通过分析高温（38℃）胁迫下半夏（Pinellia ternate）叶片转录组数据,挖掘其响应高温胁迫相关的基因,探究高温胁迫影响半夏的分子机制。【方法】通过Illumina HiSeq测序对高温胁迫处理后的半夏叶片进行转录组分析,并用qRT-PCR技术对DEGs进行验证。【结果】转录组测序共获得38.48 Gb Clean data,与对照组（25℃）相比,高温处理后1 138个DEGS上调表达,578个DEGs下调表达。GO富集分析表明,DEGs主要集中在转录因子复合体、内肽酶活性和血红素结合;KEGG富集分析表明,内质网蛋白质加工代谢通路的DEGs最多。基因表达水平分析表明,植物激素信号转导相关基因（PtARP1、PtGRP1、PtABF）、7个谷胱甘肽S转移酶基因PtGST和23个热激蛋白基因PtHSP在高温胁迫下表达量上调。qRT-PCR分析结果表明,9个DEGs在高温胁迫下的差异表达与转录组分析结果一致。【结论】通过对高温胁迫下半夏叶片进行转录组分析,挖掘高温胁迫相关基因,为进一步研究半夏高温胁迫响应机制提供理论参考。     

混合盐碱胁迫下不同抗性荷花品种比较转录组分析

为研究荷花的耐盐碱机制,以混合盐碱胁迫处理前后的抗性荷花品种黄帅TO及敏感品种台湾磨盘莲SE为试验材料,进行转录组测序,并运用qRT-PCR法验证测序结果的准确性.2个品种在混合盐碱胁迫下共检测到4 916个差异表达基因,其中,上调表达基因3 101个,下调表达基因1 815个,对抗性品种的49个特有差异基因进行分析发...     

干旱胁迫棉花转录组DNA损伤修复相关基因的分析

【目的】研究干旱胁迫下,棉花DNA损伤修复基因的表达情况,从整体水平探讨棉花DNA损伤修复相关基因表达与抗旱的相关性。【方法】用2.5%PEG6000处理棉花幼苗,利用RNA-Seq技术对干旱胁迫下棉花幼苗的转录组进行测序。【结果】从棉花干旱胁迫响应的转录组中,筛选获得差异表达的DNA损伤修复相关基因共51个,其中差异表达的上调基因23个,差异表达的下调基因28个,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。选取4个差异基因进行生物信息学分析及qRT-PCR验证,HMGB1、rec A1、UDGs和GMP synthase基因的表达变幅有一定的差异,但基因的表达趋势一致,棉花转录组测序结果通过qRTPCR验证是可靠的。【结论】DNA损伤修复相关基因可能与干旱胁迫有一定的相关性,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。     

氮胁迫下小麦叶片转录组分析

氮素是小麦生长发育过程中重要的营养元素,利用转录组技术鉴定参与小麦低氮胁迫分子调控网络的基因,对揭示小麦耐低氮分子机理、开展耐低氮育种具有重要参考意义。采用Illumina HiseqTM 2500高通量测序技术,对小麦品种晋麦47正常生长和低氮胁迫下的叶片进行转录组测序,筛选差异表达基因并在GO、KEGG 数据库中进行比对注释,分析小麦低氮胁迫响应相关基因。结果显示,对照组与低氮组分别获得高质量序列52 383 726和52 192 061条,检测到差异表达基因1 267个,其中上调表达基因179个,下调表达基因1 088个。差异基因GO功能注释到3个大类的44个功能组。差异表达基因被注释到178个途径上,主要富集于氨基酸代谢、碳水化合物代谢、脂类代谢和信号传导等途径。转录因子分析发现,在低氮条件下变化明显的转录因子家族包括WRKY、MYB和NAC等。     

氮胁迫下小麦叶片转录组分析

氮素是小麦生长发育过程中重要的营养元素,利用转录组技术鉴定参与小麦低氮胁迫分子调控网络的基因,对揭示小麦耐低氮分子机理、开展耐低氮育种具有重要参考意义。采用Illumina HiseqTM 2500高通量测序技术,对小麦品种晋麦47正常生长和低氮胁迫下的叶片进行转录组测序,筛选差异表达基因并在GO、KEGG 数据库中进行比对注释,分析小麦低氮胁迫响应相关基因。结果显示,对照组与低氮组分别获得高质量序列52 383 726和52 192 061条,检测到差异表达基因1 267个,其中上调表达基因179个,下调表达基因1 088个。差异基因GO功能注释到3个大类的44个功能组。差异表达基因被注释到178个途径上,主要富集于氨基酸代谢、碳水化合物代谢、脂类代谢和信号传导等途径。转录因子分析发现,在低氮条件下变化明显的转录因子家族包括WRKY、MYB和NAC等。     

RNA-seq技术鉴定Bt水稻对拟环纹豹蛛肠道基因表达的影响

为了探究Bt水稻对拟环纹豹蛛体内基因表达水平的影响,利用Illumina 2000测序平台对蜘蛛肠道进行转录组测序并分析其基因表达情况,从分子角度揭示Bt水稻的生物安全性。对肠道测序共获得65 028 186条clean reads,经Trinity软件拼接得到135 829条长度大于200bp的Unigene用于整个转录组分析。用FPKM法计算2个样本中基因的表达量,并以负二项分布检验法分析其显著性。结果显示,富集Bt的拟环纹豹蛛体内相对于正常蜘蛛显著差异表达的基因有142个,包括1个下调141个上调,并且其差异表达倍数都在8倍以上(FDR0.05,Log2Ration2)。利用Blast程序将差异基因与COG,GO和KEGG数据库进行比对,对差异基因进行功能注释。注释结果描述的酶、生物学过程和代谢途径与细胞内产能反应相关,如细胞色素c氧化酶,氧化磷酸化等。推测Cry1Ab可能对拟环纹豹蛛体内能量代谢相关基因的表达有一定的影响。     

匍匐翦股颖接种立枯丝核菌后基因表达变化的转录组学分析

【目的】确定匍匐翦股颖(Agrostis stolonifera)接种立枯丝核菌(Rhizoctonia solani)后基因种类和表达量在转录水平的变化规律,明确草坪草病原菌侵染响应的关键基因。【方法】匍匐翦股颖生长14 d后采用麦粒培养物接种立枯丝核菌,接种3 d后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,用Trinity组装匍匐翦股颖转录组,以组装子为参考,以|log2(fold change)|1,q-value0.005为阈值选取感病和健康匍匐翦股颖叶片转录组的差异表达基因,并用i TAK软件分析其转录因子家族及表达变化,与植物R基因库进行blast分析R蛋白分类、利用Mapman软件分析生物胁迫信号通路相关基因的表达变化。【结果】高通量测序得到125 253 092条高质量待分析reads。经Trinity从头组装后,得到466 761条转录本。过半数的转录本长度为700 bp以上,组装结果 N50=1 100 bp。使用CD-HIT选择334 212条转录本(所有转录本的71.60%)作为Unigene,平均长度573 bp,N50=791 bp。接种后植物比接种前植物基因有7 937个上调表达,1 570个下调表达。上调基因中296个,下调基因有142个都可被定义为转录因子,分布在58个转录因子家族中,其中锌指蛋白包含转录因子C2H2最多,有54个,C3H次之,为22个。差异基因中451个可定义为植物R蛋白表达基因,可分为33类,其中包含NBS-LRR结构域的抗病蛋白、LRR受体蛋白激酶、ABC-2类型的转运蛋白、U-box结构域蛋白激酶和热激蛋白这5类基因变化最显著。差异基因中大量上调表达基因可富集在病原识别、活性氧消除、信号传导、细胞凋亡、病程相关蛋白等生物胁迫相关基因类别,下调基因显著富集在植物生长发育相关类别和通路。q RT-PCR验证了随机挑选差异基因的表达量变化,均与RNA-seq分析结果一致,其中包括12个C2H2转录因子基因,10个C3H转录因子基因,以及12个R蛋白基因。【结论】病原菌侵染后,引起匍匐翦股颖大量基因表达变化,其中转录因子、R蛋白以及抗性相关基因多为上调表达,作用为抑制病原菌扩散,而生长发育相关基因表达下调,这些进程共同使匍匐翦股颖产生对立枯丝核菌的先天基础抗性。     

基于高通量测序技术分析生物质炭可溶性组分处理不结球白菜叶片的转录组学分析

[目的]本文旨在探究生物质炭表面可溶性组分对作物生长的潜在作用机制。[方法]采用新一代高通量测序手段——RNA-Seq技术对添加生物质炭可溶性组分处理和未添加可溶性组分的不结球白菜叶片进行转录组分析,研究其差异基因的功能,并利用实时荧光定量PCR技术(RT-qPCR)对关键基因进行检测验证。[结果]添加生物质炭可溶性组分的不结球白菜叶片中共检测到499个差异基因,其中上调表达基因152个,下调表达基因347个。RT-qPCR验证分析与转录组测定结果一致。基因功能注释显示差异基因多集中于细胞壁的多糖代谢和碳水化合物代谢等生物过程,主要调控果胶裂解酶、果胶甲酯酶、葡糖基转移酶等活性。[结论]果胶裂解酶和果胶甲酯酶同源基因的下调表达、葡糖基转移酶活性同源基因的上调表达证实生物质炭可溶性组分可以在增加植物细胞壁稳定性的同时提高作物对逆境条件的抵抗能力,从基因表达的角度揭示了生物质炭可溶性组分对不结球白菜的影响。     

腾冲红花油茶蜜腺2个发育时期转录组差异性分析

【目的】植物蜜腺基因能够调控其分泌产物花蜜中的化学组成,从而影响植物对传粉者的吸引力,本研究主要揭示腾冲红花油茶蜜腺发育与分泌过程中其基因表达的变化情况。【方法】采用高通量Illumina Hiseq测序手段对花蕾期蜜腺与泌蜜高峰期蜜腺进行转录组测序,并对测序数据进行相关分析研究,探讨蜜腺发育与分泌过程中的调控基因。【结果】在蜜腺发育与分泌过程中共检测到58 488个差异基因,其中表达上调的基因有1 602个,表达下调的基因有693个。GO分析结果表明:被注释生物学过程的差异表达基因有3 096个,其中富集差异基因数量较多的生物学过程为胞内过程、生物调节过程和新陈代谢过程;富集于分子功能的差异表达基因有1 280个,其中富集差异基因数量较多的分子功能为连接功能、催化活性和载体活性。Pathway显著性富集分析发现:有2 147个差异表达基因参与了292条代谢途径,其中富集差异基因数量最多的途径是代谢途径(941个,占43.83%),主要是参与糖代谢、氨基酸代谢等过程。【结论】蜜腺发育及花蜜分泌过程中大多数差异表达基因与花蜜化学组成的调控过程有关。