本研究通过冰山雪莲(Paeonia rockii’Bing Shan Xue Lian’)和牡丹主栽品种凤丹(Paeonia ostii’Phoenix White’)根组织转录组测序进行差异表达基因分析,结果显示,冰山雪莲和凤丹基因表达量相近,差异表达基因487条,其中228条上调,259条下调;冰山雪莲特有表达基因数113 558,凤丹58 114,共有表达基因数55491;两品种4 228条低表达量基因聚集成2个子集,6 259条高表达量基因聚集成2个子集;GO富集注释差异基因1 487条,其中上调基因数1 137,下调基因数1 200,上调差异基因的主要功能与氧化还原和营养库活性有关;差异基因KEGG显著富集的代谢通路73条,富集前20的代谢通路中参与核糖体通路的基因最多,富集因子最高的代谢通路为硫胺素代谢、其次为C5-分支二元酸代谢和二萜类生物合成,参与植物病原体互作、谷胱甘肽代谢和RNA降解的上调差异基因数目相对较多,其次为异喹啉、托烷、哌啶、吡啶生物碱和萜类骨架合成、苯丙氨酸、酪氨酸、甘氨酸、丝氨酸和苏氨酸、肌醇磷酸和戊糖磷酸代谢、DNA复制和核苷酸剪切修复。通过转录组测序全面分析了冰山雪莲和凤丹根组织差异基因表达状况,为牡丹功能基因的研究提供了理论依据。