基于RNA-seq的番茄响应根结线虫MiPDCD6蛋白的SA信号通路基因鉴定与表达分析

为了揭示南方根结线虫MiPDCD6蛋白抑制番茄PTI免疫的分子机理，以番茄品种新金丰1号MiPDCD6超表达苗为试验材料，以番茄品种新金丰1号组培苗为对照，通过转录组测序技术分别对番茄MiPDCD6超表达苗和对照苗进行转录组测序。以番茄栽培品种Heinz 1706基因组作为参考基因组进行对比，利用FPKM法计算基因表达量，设定参数(|log₂ FC|>1且P<0.05)筛选差异表达基因。利用Gene Ontology(GO)数据库进行差异表达基因GO功能富集分析，统计每个GO term中的差异表达基因数目，计算出基因富集的显著性，找出富集显著的功能条目。利用KEGG数据库对差异表达基因进行Pathway富集分析，超几何分布检验方法计算每个Pathway中差异表达基因富集的显著性。以FDR和基因数量来衡量KEGG的富集程度。结合基因差异表达分析和功能富集分析，研究MiPDCD6蛋白对番茄PTI免疫相关通路基因的影响。结果发现：在超表达MiPDCD6番茄植株中，与野生型番茄相比，MiPDCD6过表达番茄植株中有2 366个差异表达基因(DEGs),其中1 ...     

基于转录组测序筛选鸡皮肤毛囊性状相关基因和关键通路

旨在通过对崇仁麻鸡背部皮肤毛囊组织的转录组学分析，筛选关键的差异表达基因和信号通路，探讨鸡皮肤毛囊性状的分子调控机制。采用TRIzol法提取皮肤毛囊组织总RNA,通过Illumina平台进行转录组高通量测序(RNA-Seq)。最后采用实时荧光定量的方法对转录组测序进行验证。研究结果如下：共筛选出了3 856个差异表达基因，在公鸡(BGB vs YGB)比较组合中筛选出了971个差异表达基因，其中274个表达上调，697个表达下调。在母鸡(BMB vs YMB)比较组合中筛选出了3 529个差异表达基因，其中1 477个上调，2 052个下调。在公鸡、母鸡共同的差异表达基因中筛选出了3个与皮肤毛囊性状有关的高表达基因KRT75、KRT6A和KRT14。在GO功能显著富集中，BGB vs YGB比较组合富集到了516个GO term, BMB vs YMB比较组合富集到了1 020个GO term,基因主要集中参与细胞黏附等活动。显著富集了多条KEGG通路，在BGB vs YGB比较组合中富集了8条显著的KEGG通路，在BMB vs YMB比较组合中富集了20条显著的KEGG通路，筛选出了...     

基于转录组测序技术挖掘影响布莱凯特黑牛睾丸大小的功能基因

旨在通过对布莱凯特黑牛大小睾丸转录组数据进行分析,挖掘影响睾丸生殖功能的关键候选基因,探究睾丸大小对牛生殖功能维持的分子机制。采集12月龄布莱凯特黑牛睾丸,根据质量、长径和短径将其分为2组(大小睾丸组各3头)进行高通量转录组测序;以|log₂(Fold Change)|≥1且P-value≤0.01为阈值筛选差异基因,对差异基因进行GO和KEGG功能富集分析,筛选出影响睾丸生殖功能的关键基因;对其进行CDS区克隆测序,并运用qRT-PCR、Western Blot、免疫组化等方法进行进一步研究。结果显示,差异表达基因共353个,其中114个基因在小睾丸中表达上调,239个基因表达下调。GO功能注释表明,差异基因参与代谢、发育、生殖等过程;KEGG富集分析显示,差异基因主要富集在Wnt信号传导、PI3K-Akt信号通路等与睾丸生殖功能相关的途径。最终筛选出睾丸生殖功能基因CCN4进行进一步研究。qRT-PCR和Western Blot结果显示,CCN4基因和蛋白在大睾丸中表达量显著低于小睾丸。免疫组化结果显示,CCN4蛋白主要在精细胞外的区域(特别是支持细胞)中表达...     

基于转录组测序筛选黄牛低氧适应性相关差异基因

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总RNA,建立测序文库后利用Illumina HiSeqTM 4000进行测序,经生物信息学分析筛选出差异表达转录本,并对差异表达基因进行GO富集和KEGG通路分析,最后采用荧光定量PCR(RT-qPCR)验证测序数据的准确性。结果表明:藏黄牛相比于三江黄牛得到显著差异表达基因608个,其中上调基因467个,下调基因141个,推测表达量较高的B2M、COX7C、NFATC1、ECE1、CAST基因是与低氧适应性相关的基因;差异转录本GO功能富集分析可知,在分子功能、细胞组分、生物过程分别富集到327,253,2 191个条目,其中富集程度最高的条目分别是细胞发育过程、结合与转录因子活性、细胞前缘。KEGG分析结果表明,差异转录本显著富集在8条通路,筛选低氧适应性候选基因NFATC1是MP-PKG信号通路与T细胞受体信号通路成员,推测其通过cGMP-PKG通路与T细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

转录组学鉴定高粱营养器官间功能差异的决定基因

本研究采用Illumina HiSeq 4000高通量测序技术,对高粱根、茎和叶的转录组进行测序以及生物信息学分析,以期鉴定出决定营养器官功能异质性的核心基因。结果表明,转录组测序获得165 413 338条原始测序序列(Raw reads),过滤后获得160195618条序列(Clean reads);共鉴定出33204个基因及230685个可变剪接事件;三种营养器官基因表达模式的相似度较低,共筛选出10 083个差异表达基因,其中在不同营养器官均有表达差异的决定基因1320个。对决定基因进行GO富集分析,得到2946个GO功能注释(p0.05),分布于光合作用、结构分子活性和细胞组成相关的GO项中;KEGG富集分析表明差异基因在光合作用、黄酮和黄酮醇生物合成以及叶酸参与的"一碳单位"代谢途径富集程度显著。本研究为深入探讨高粱营养器官间功能差异性机制及器官特异性启动子克隆等工作奠定了数据基础。     

2种黑籽南瓜响应枯萎病菌侵染的转录组学研究

为了筛选白皮黑籽南瓜和绿皮花纹黑籽南瓜抗枯萎病的差异表达基因,挖掘参与抗性相关的代谢通路。以尖孢镰刀菌侵染2 d后和未侵染的2种幼苗叶片作为材料,利用Illumina Hiseq 2500测序技术进行转录组测序分析。以未侵染的幼苗作为对照,白皮和绿皮花纹黑籽南瓜差异表达基因分别有2082和1116个。其中141个基因在2种黑籽南瓜中均呈现出差异表达,62个基因具有相同表达趋势。对差异表达基因进行GO功能注释,结果显示2种黑籽南瓜差异表达基因主要富集在生物学过程中。KEGG代谢通路分析发现白皮有22个差异表达基因显著富集在植物激素信号通路,绿皮花纹有30个差异表达基因富集在苯丙烷生物合成通路。本研究挖掘出黑籽南瓜抗枯萎病差异表达基因,预测了抗性通路,为瓜类抗病机制的深入研究奠定了理论基础。     

盐胁迫下不同甘薯品种的转录组测序分析

为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理0 d、3 d和6d的耐盐甘薯品种‘榕薯819’以及不耐盐甘薯品种‘榕薯910’的叶片为材料,借助高通量测序技术进行转录组测序分析。结果表明, 2个品种共获得157,252条Unigenes,平均组装长度为576 bp。其中有83,264条Unigenes在七大数据库中得到注释,占总数的52.95%。NR注释分类结果显示,在牵牛花(Ipomoea nil)中比对到同源序列的Unigenes最多,共43,620条,占总数的57.05%。Unigenes在KOG数据库中的注释主要富集在普通功能预测(8752个)、信号转导机制(5067个)以及翻译后修饰、蛋白转换、分子伴侣(4471个)中。差异表达分析显示,在‘榕薯819’中,盐处理3 d和6d的样品差异表达基因数分别为323个和3752个,共参与了33个GO功能分类项和302条KEGG代谢通路。在‘榕薯910’中,差异表达基因数则分别为5554个和7395个,共参与了50个GO功能分类项,涉及了329条KEGG代谢通路。以部分差异表达基因的转录组数据为基础,...     

金花葵开花前后差异基因及代谢物分析

为深入研究金花葵花朵开花前后基因和代谢物的变化,利用转录组学和代谢组学技术相结合的方法对金花葵花蕾和花朵进行检测。结果显示,通过转录组分析鉴定了206 636个Unigenes,筛选出42 618个差异表达Unigenes,其中包括63个差异表达转录因子家族,24个转录调节因子家族。GO分析结果显示,差异表达基因主要富集于囊泡介导的逆行运输等生物学过程。KOG功能注释显示,差异表达基因功能以通用功能预测居多,其次是信号转导机制、翻译后修饰蛋白周转以及碳水化合物的转运和代谢。差异表达基因KEGG富集分析表明,差异基因主要富集于代谢、植物激素和信号转导、淀粉和蔗糖代谢等代谢通路。通过代谢组学检测,筛选到差异显著代谢物135个,主要包括脂类、氨基酸及衍生物和黄酮类等,KEGG富集分析表明,差异代谢物主要富集于甘油磷脂代谢、糖基磷脂酰肌醇(GPI)-锚定生物合成和次生代谢产物的生物合成-未分类等过程,其中,缬氨酸,亮氨酸和异亮氨酸的降解代谢通路同时出现在了基因和代谢物KEGG代谢通路富集Top 20中。     

玉米自交系响应花粒期高温胁迫差异表达基因的分析

为了探明玉米重要自交系花粒期高温胁迫下转录组基因的表达差异,发掘玉米响应高温胁迫的关键基因和蛋白质,明确高温胁迫引起玉米减产的机理,从而利用分子生物学技术手段提高玉米耐高温胁迫性。首先,利用Ampha~(TM)Z30花粉活性分析仪检测高温胁迫和正常生长条件下玉米自交系昌7-2、郑58的花粉活性。然后收集花粉提取总RNA,利用Illumina Hiseq2000~(TM)高通量测序技术分别对昌7-2、郑58花粒期高温胁迫材料和正常生长条件下的对照材料进行全转录组测序,序列结果分析后获得差异显著表达基因。通过差异基因维恩图(VENN图)重叠区域分析、GO功能分类和富集统计、KEGG pathway通路分类和富集分析以及转录因子分析,获得玉米花粒期高温胁迫相关的重要差异表达基因和蛋白质。花粉活性检测结果显示,高温胁迫后和正常生长条件下,昌7-2花粉活性分别为24. 37%,42. 89%,郑58花粉活性分别为35. 57%,64. 83%。高温胁迫后在昌7-2花粉转录组中共检测到4 176个显著差异表达基因,郑58花粉转录组中共检测到5 487个显著差异表达基因。KEGG pathway通路分类和富集分析结果显示,高温胁迫后郑58的花粉转录组中有2 062个差异表达基因富集到399个相关通路上,昌7-2的花粉转录组中有1 943个差异表达基因富集到352个相关通路上。转录因子分析结果表明,共获得55个转录因子家族,其中,有45个转录因子包含10个以上差异表达基因。研究表明,通过对玉米自交系花粒期高温胁迫后花粉转录组测序获得了大量的差异显著表达基因信息,昌7-2、郑58对高温胁迫响应机制存在明显差异,热激蛋白(Hsps)、热激转录因子(Hsfs)、生长素(IAA)基因和磷脂酰基醇-4,5-二磷酸基因家族(PIP2)在响应玉米花粒期高温胁迫过程中起着重要作用。     

基于转录组筛选肉牛骨骼肌差异基因

旨在基于转录组筛选肉牛骨骼肌差异基因,探究差异基因对骨骼肌转录调控的影响,鉴定其与肌肉生长发育的关系,为研究骨骼肌生长发育的分子机制提供理论依据。采取12月龄的布莱凯特黑牛和鲁西黄牛背最长肌,利用Illumina 2500进行转录组测序,通过综合生物信息学分析(差异基因筛选、GO分类、KEGG富集分析和蛋白互作网络等)筛选骨骼肌差异的调控基因。布莱凯特黑牛和鲁西黄牛2个文库分别存在13 242,13 350个基因,共计13 758个基因。5 500个基因表达量相同,约占总参考基因的39. 98%,8 258个基因表达量有所不同,占比60. 02%。获得554个差异显著的基因作为差异表达基因,其中296个基因表达上调,258个基因表达下调。GO功能注释结果显示,生物学过程、细胞组分和分子功能分别包括22 641,7 720和4 414个基因,其中29个基因与肌细胞的发育和分化有关; KEGG富集到Wnt、CMC、TGF-β、DCM、ROAC、VSMC 6条肌肉生长发育相关通路,筛选骨骼肌候选基因MYL2、MYL3和MYLPF。MYL2、MYL3和MYLPF是Wnt、DMC和TGF-β信号通路中的成员,参与了肌细胞和肌肉组织的形成的过程,推测MYL2、MYL3和MYLPF在骨骼肌生长发育过程中发挥重要作用。     

转录组鉴定‘三白’石榴品种籽粒发育相关差异表达基因

为了研究石榴籽粒发育遗传机制及其在育种中的应用,以‘三白’石榴品种开花后60 d和120 d籽粒为试材,进行转录组测序分析。对表达差异基因进行GO注释及KEGG富集分析。结果分析表明,与‘三白’石榴开花后60 d籽粒相比,在开花后120 d籽粒中总共检测到8 409个显著差异表达基因;KEGG富集分析表明差异表达基因显著富集在激素信号转导及其它代谢产物的合成和代谢等18条代谢通路上;5类可变剪接类型中鉴定出了76个靶基因,其中17个为差异表达基因。定位到的差异表达基因可以为解析石榴籽粒发育以及基因克隆提供理论指导。     

灵武长枣嫁接和根蘖繁殖植株果实转录组差异分析

为解析嫁接繁殖和根蘖繁殖灵武长枣果实糖和酸转化差异,本研究以两种繁殖方式的成熟期果实为材料,利用高通量测序平台Illumina/Nova-Seq 2500进行了转录组测序,对转录组数据质量评估获得两样品的序列长度分别为24.48 Gb (嫁接)和28.64 Gb (根蘖),错误率为0.03%,GC含量44.11%～44.56%,Q30大于90.54%,转录组数据与冬枣基因组序列比对,所有的转录组Reads均分布于冬枣12条染色体上,基因定位的Reads数量嫁接繁殖比根蘖繁殖多;差异表达基因分析筛选到27个显著差异表达基因,其中9个显著上调,18个显著下调;GO富集分析表明15个差异基因显著富集到GO的三大功能类别分子功能、细胞组分、生物过程中;KEGG富集分析表明10个差异表达基因显著富集到12条代谢通路中,其中包括3个上调基因和7个下调基因;进一步筛选到调控糖及有机酸代谢关键差异基因BAM1、PCO和HAOX1在淀粉和蔗糖代谢、牛胆素和亚牛磺酸代谢、乙醛酸和二羧酸代谢中显著表达,并在糖和有机酸的相互转化中起着关键调控作用。本研究表明嫁接繁殖和根蘖繁殖灵武长枣果实在糖和有机酸代谢存在差异,为灵武长枣优质果品繁育提供理论基础。     

加工番茄耐盐突变体耐盐相关基因的转录组分析

为筛选鉴定番茄叶片响应盐胁迫应答基因,本研究通过甲基磺酸乙酯(methyl ethyl sulfonate, EMS)诱变处理加工番茄‘JW9’获得的16株耐盐突变植株。经单独收种后,利用(RNA sequencing, RNA-Sep)技术对获得的耐盐突变体加工番茄幼苗叶片进行转录组测序分析。提取样品总RNA进行转录组cDNA文库制备后进行Illumina HiSeq/MiSeq测序。对转录组测序结果进行分析共获得了102 Gb大约810 974 186个高质量的Clean阅读序列,Cufflinks软件重构转录本后得到38 190个已知转录本,1 647个未知转录本。在测序样本中共筛选出4 367个基因在样品间显著差异表达,其中上调表达基因2 606个,下调表达基因1 761个。经qRT-PCR检测表明,10个随机选择的差异表达基因(differentially expressed genes, DEGs)表达量变化与转录组测序结果一致。通过GO和KEGG通路富集分析,从具有耐盐功能的丝氨酸/苏氨酸蛋白激酶和ATPase酶合成相关的基因中,找到5个显著差异表达基因Solyc08g066270.1和Solyc03g083420.2以及Solyc10g018530.1、Solyc09g082890.1和Solyc02g014190.2可能为耐盐相关基因。本研究从耐盐突变体与原始亲本的差异表达基因中筛选出耐盐候选基因,为进一步鉴定和克隆出番茄的耐盐基因培育耐盐品种番茄提供理论参考。     

干旱胁迫下枇杷叶片的转录组分析

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行de novo拼接、功能注释和ORF预测,将差异表达基因在COG、GO和KEGG数据库中进行比对注释。测序结果表明,获得转录本共88 530个,平均长度为740.64 bp,ORF41 748条。COG、GO和KEGG数据库将转录本分别划分为24,54个功能类别及291条代谢通路中。25 197个差异表达的基因在30条代谢通路中显著富集,与酶活性、激素合成代谢和信号转导等相关的差异基因积极响应枇杷干旱,其中双萜类、油菜素类固醇和类胡萝卜素3条生物合成途径中的差异基因呈现出较为一致的表达。采用实时荧光定量(qRT-RCR)对选取的差异基因进行验证,其中过氧化物酶、蛋白激酶byr2、丝氨酸苏氨酸蛋白激酶、脱落酸8'-羟化酶、吲哚-3-乙酰乙酸合酶相关基因上调表达,细胞色素P450 734A1基因下调表达。为枇杷提供了较为全面的基因信息和代谢途径数据,为干旱胁迫下枇杷分子调控机制的深入研究奠定了基础。     

不同硬度西瓜果皮的转录组测序及相关基因表达分析

为解析西瓜果皮硬度不同的相关分子机制,挖掘西瓜果皮硬度相关的关键基因提供理论基础,以生育期相似但果皮硬度差异显著的高硬度（901）4-1-1-M和低硬度BSH为试验材料,选取果皮硬度差异最大时期即西瓜授粉30 d后的果皮,利用Illumina HiSeq^TM测序平台进行转录组测序分析,采用实时荧光定量的PCR技术对测序结果进行验证。结果显示,转录组测序共筛选获得1 085个差异表达基因,其中,上调表达基因555个,下调表达基因530个。GO富集分析表明,1 085个差异表达基因显著富集在细胞组分、分子功能和生物过程中的细胞壁、细胞外围、外部封装结构、胞外区、四吡咯骨架、氧化还原酶活性、转移酶活性、果胶酯酶活性等功能类别。KEGG富集分析表明,Cla97C04G075800、Cla97C02G044950、Cla97C09G165820、Cla97C10G195660、Cla97C01G025380等19个差异基因被注释富集程度最高的苯丙烷代谢通路,该路径最终形成木质素、5羟基愈创木酚木质素、愈创木酚木质素和对羟基苯酚木质素代谢物,qRT-PCR和转录组测序数据相关...     

棉花无腺体近等基因系差异表达基因分析

为了揭示棉花腺体发育的分子机制,本研究利用RNA-Seq对棉花无腺体性状近等基因系Z12/Z12YW的幼嫩叶片进行转录组测序,寻找差异表达基因。转录组测序结果表明,与有腺体棉Z12相比,无腺体棉Z12YW中共有306个基因表达发生变化,其中282个基因下调,24个基因上调。差异表达基因GO(Gene ontology)功能注释显示,细胞组分、胞外基质组分的功能变化以及细胞杀伤生物学过程在腺体形成过程中发挥重要作用;COG(Cluster of orthologous groups)功能和KEGG(Kyoto encyclopedia of genes and genomes)通路显著性富集发现,棉酚等萜烯物质的合成主要受甲羟戊酸途径上游基因的表达调控。此外,本研究共筛选到13个差异表达转录因子,可能在腺体发育和棉酚合成过程中发挥关键作用。     

基于转录组测序技术挖掘大豆蛋白质合成相关基因

通过转录组测序技术挖掘蛋白质合成相关基因,为解析蛋白质合成机制奠定基础。以遗传背景相近但蛋白含量差异较大的大豆高蛋白品系冀HJ117(蛋白质含量52. 99%)及其回交亲本冀豆12(蛋白质含量46. 48%)为研究材料,利用转录组测序技术和生物信息学分析,以期挖掘大豆籽粒蛋白质合成相关基因。通过对转录组数据中冀HJ117和冀豆12差异表达基因的分析筛选,共得到336个差异表达基因,其中冀HJ117较冀豆12有195个上调表达基因,141个下调表达基因。通过GO功能富集分析,发现在分子功能类型中注释的差异基因主要与催化和连接等功能相关; KEGG显著富集分析发现差异表达基因主要富集在蛋白质内质网合成途径中,该途径共筛选到34个差异表达基因,其中33个基因在冀HJ117中表达量较高。从该途径中筛选出9个候选基因,采用荧光定量PCR方法检测其表达量,表达趋势与RNA-Seq检测结果基本一致,证实了RNA-Seq数据的可靠性。RNA-Seq测序结果获得了与蛋白质合成相关基因的基本信息,蛋白质内质网合成途径可能是冀HJ117与冀豆12蛋白质含量产生差异的重要通路。     

玉米发芽期响应盐胁迫的转录组分析

玉米种子发芽期对盐分较为敏感,为了进一步明确玉米发芽期对盐胁迫的应答机制,本研究利用高通量转录组测序技术,对玉米耐盐自交系昌7-2在150 mmol/L NaCl胁迫下的胚芽进行转录组测序和数据分析。结果表明：共有40 290个Unigenes获得功能注释,其中差异表达基因615个。在差异表达基因的GO富集分析中,对富集程度明显的前20个类别做出功能分析。有21个DEGs富集到苯丙烷类的生物合成途径,是KEGG富集分析最显著的通路。对差异表达基因进行COG分类,共得到25个不同的COG功能注释。通过对目标基因的挖掘共得到13个较重要的转录因子,涉及到MYB、WRKY、AP2、bZIP、bHLH和NAC等6个家族。本研究为挖掘玉米耐盐相关基因、解析玉米响应盐胁迫的分子调控机制提供科学依据。     

盐、低磷以及干旱胁迫下乌拉尔甘草根转录组分析

对乌拉尔甘草根进行转录组测序分析,挖掘响应盐、低磷和干旱胁迫的相关基因,探讨甘草抵御逆境胁迫的分子机制。以150 mmol/L NaCl模拟盐胁迫、10μmol/L KH₂PO₄模拟低磷胁迫和15%聚乙二醇(PEG-6000)模拟干旱胁迫处理下的乌拉尔甘草根作为研究材料,应用Illumina HiSeq 2000测序平台对根进行转录组测序,对测序结果进行基因功能注释分析和差异表达基因(DEGs)筛选。分别获得盐、低磷和干旱胁迫处理下甘草根中4 294、3 201和4 719个DEGs。GO富集结果表明,3种逆境胁迫下甘草根中的DEGs均显著富集到细胞过程、代谢过程、刺激响应、细胞、细胞器、催化活性等功能过程中;KEGG途径富集分析表明,3种胁迫处理下的DEGs显著富集的代谢途径主要涉及基因表达调控、生物合成及代谢、信号转导等。转录因子鉴定结果显示3种胁迫处理下差异表达数量最多的转录因子主要来自ERF、bHLH、MYB-related、WRKY、Dof、NAC等转录因子家族;基因表达水平分析显示,参与植物激素代谢及信号转导相关的DEGs在盐胁迫...     

银杏三个年龄阶段雄花芽的转录组分析及差异基因筛选

本研究筛选银杏三个年龄阶段开花调控的关键基因,揭示银杏开花调控年龄途径的分子机理,为促进银杏早花、分子育种和推广种植提供科学依据。本研究利用高通量测序技术及生物信息学工具,对三个年龄阶段的银杏雄花芽进行转录组测序,对测序数据进行分析,筛选出年龄途径关键差异表达基因。转录组测序转录组测序共产生57.45 Gb原始数据,注释到8大功能数据库(GO, COG, KEGG, KOG, NR, Pfam, SwissProt, egg NOG)上的Unigene总数为35 058个。通过KEGG通路注释和GO富集分析,将Unigene分别归于55个GO类别和126个代谢途径。差异表达基因分析显示,幼年雄花芽较开花一年雄花芽有37个基因上调,75个基因下调;幼年雄花芽较开花多年雄花芽有592个基因上调,871个基因下调;开花一年雄花芽较开花多年雄花芽有961个基因上调,1 203个基因下调。本研究发掘出大量的开花相关的基因,最终筛选出SPL (gene.Gb_23724, gene.Gb_03922)、AP2 (gene.Gb_00766)、MADS-box (gene.Gb_01886, gene.Gb_15398, gene.Gb_28337, Gingko_newGene_2213)、GA调节蛋白(gene.Gb_34467, gene.Gb_28606, gene.Gb_33214)和DELLAs蛋白(gene. Gb_34644)等11个年龄途径开花调控的关键基因。