摘要：【目的】以自然突变的楸树雄性不育花芽为研究材料,分析与楸树雄性不育相关基因的表达模式,从基因的表达水平上揭示楸树雄性不育的分子机制,为楸树及木本植物的雄性不育研究提供有价值的参考。【方法】采用高通量测序技术对自然状态下楸树的雄性不育的花芽以及可育的花芽进行转录组测序,通过生物信息学对可育花芽与不育花芽的转录本进行比较分析,预测和筛选出与楸树雄性不育有关的基因。【结果】转录组测序共产生27.18 Gb数据,组装并去冗余后得到86 076个Unigene,然后将Unigene比对到7大功能数据库（NR,NT,GO,COG,KEGG,Swissprot,Interpro）进行功能注释,最终被7大数据库中任意一个数据库注释上的Unigene总数为64 600（75.05%）。将试验组（雄性不育花芽,SL）与对照组（可育花芽,FL）所测得的Unigene进行表达量的分析后,筛选出表达量有差异且可信度高的差异表达基因。在噪音分布中（差异表达倍数在2以上,可信度在0.8以上）共筛选出试验组中有6 915个基因上调表达,3 504个基因下调表达。在泊松分布中（差异表达倍数在2以上,错误发生率在0.001以下）,SL-1 vs FL-1、SL-2 vs FL-2、SL-3 vs FL-3三个生物学重复中得出差异上调表达基因分别有13 979,13 513,13 055个,差异下调表达基因分别有12 170,13 807,10 411个。差异基因的GO功能分析表明,在生物过程中显著性富集的条目集群频率较高的是生殖过程、生殖发育过程、生殖系统发育、生殖结构发育,在分子功能中只有生长素流出跨膜转运蛋白活性显著性富集。差异基因KEGG通路分析中,差异基因映射到127个不同的生物途径,显著富集代谢通路中差异基因的生物途径主要为：代谢途径、次生代谢产物的生物合成、剪接体、RNA转运以及甘油磷脂及淀粉和蔗糖的代谢。在用差异基因与同源基因比对分析中,共比对出246个同源性高�

注：因版权方要求，不能公开全文，如需全文，请咨询杂志社