为了揭示棉花腺体发育的分子机制,本研究利用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个差异表达转录因子,可能在腺体发育和棉酚合成过程中发挥关键作用。