低温胁迫对冰菜转录组水平响应分析

利用二代高通量测序技术对低温胁迫处理的冰菜进行测序,构建冰菜转录组数据库.分别得到24.13 Gb有效数据和24045条Unigene的注释,得到DEGs 1902个(T0 vs T1)和2134个(T0 vs T2).T0 vs T1组和T0 vs T2组分别有40和41个功能小类化归GO数据库;分别有20和24个功能分类注释到KOG数据库;155和272条基因注释到KEGG数据库,并分别富集在74和105条代谢通路.T0 vs T1组DEGs主要注释到植物信号转导等4个代谢通路;T0 vs T2组DEGs注释到苯丙醇类生物合成等11个代谢通路,其中正向影响代谢途径:丙酮醇类生物合成、嘌呤代谢、谷胱甘肽代谢、脂肪酸代谢、类黄酮代谢、氨基酸的生物合成代谢途径等;负向影响代谢途径:植物-病原互作、植物激素信号转导、淀粉和蔗糖代谢等途径.通过对淀粉和蔗糖代谢途径关键基因分析表明:低温胁迫1 h (T1),海藻糖6-磷酸合酶、海藻糖-6-磷酸酯酶、β-淀粉酶、葡萄糖-1-磷酸腺苷酰转移酶、糖原磷酸化酶等5个关键基因表现为上调表达,未见下调表达基因;低温胁迫36 h (T2),海藻糖-6-磷酸合酶、己糖激酶、β-淀粉酶等3个关键基因上调表达,葡萄糖内酯-1,3-β-葡萄糖苷酶基因下调表达.选取淀粉和蔗糖代谢途径中8个DEGs,经RT-qPCR分析,8个DEGs的相对表达量与转录组表达水平相符.     

狭叶薰衣草叶片转录组耐寒相关差异表达基因分析

以生长在哈尔滨地区的狭叶薰衣草(Lavandula angustifolia Mill.)为试材,采用高通量转录组测序的方法,研究了田间不同温度(26.67℃, F; 32.07℃, S;-0.76℃, T)对薰衣草叶片基因差异表达的影响,以期挖掘与耐寒相关的功能基因,探讨薰衣草耐寒的分子机制。结果表明,温度对薰衣草差异表达基因具有显著影响,温差越大导致基因差异表达变化越大,差异表达基因与环境因子CCA分析表明温度是诱导基因差异表达的首要因素。在KEGG富集通路中显著富集的有植物与病原体互作通路、MAPK信号通路和昼夜节律通路,在这些通路中,基因Pti5、SUGT1、WRKY25、GI、FKF1等均响应植物低温胁迫并显著表达,且昼夜节律调控对薰衣草耐寒起到重要作用。利用STRING蛋白质互作数据库筛选出KEGG通路中的抗寒相关枢纽基因有FKF1、GI、CRY2等。从转录因子家族中鉴定出1 137个转录因子,隶属于51个转录因子家族,参照拟南芥数据库筛选出的枢纽基因AT4G30825、AT4G17020、AT1G55750均来自bZIP家族,且在KEGG通路中响应低温后表达量普遍升高。通...     

RNA-seq揭示碱性盐(NaHCO3)对细叶百合鳞茎基因表达的影响

本研究旨在揭示细叶百合在盐碱逆境中的基因表达表达情况,为合理利用盐碱地和大面积种植百合提供理论基础。以一年生的细叶百合鳞茎为材料,经过20 mmol/L NaHCO₃处理24 h后,用Illumina HiSeqTM2000测序平台进行转录组测序。共测得56828个mRNA的注释信息,其中,55433个Unigene被注释到Nr数据库,26973条Unigene被注释到KOG数据库,23610条Unigene被注释到GO数据库,13142条Unigene被注释到KEGG数据库的五大类中。共鉴定出390个差异表达基因。选取了9个与盐碱胁迫密切相关的基因进行qPCR验证,9个基因的表达结果与转录组的结果基本趋于一致。通过细叶百合在碳酸盐(NaHCO₃)逆境下的转录组对比数据,提供了许多的基因表达通路和表达量的差异,为合理利用盐碱地和大面积种植百合提供理论基础。     

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

为了获得甘薯耐盐转录组序列信息,挖掘差异表达基因及其相关代谢途径,本文以盐胁迫处理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代谢通路。以部分差异表达基因的转录组数据为基础,...     

金线莲应答高温胁迫的转录组学特征分析

为了研究高温胁迫下金线莲转录组表达特征，分析其热激响应机制。以正常培养和高温胁迫(45^oC)的金线莲(NYJ2)为材料，利用Illumina HiSeq^TM2000平台进行测序，通过Trinity软件进行De novo组装。结果共获得75688个unigene，有28323个unigene在Nr、KOG、KEGG和Swisspro数据库中得到功能注释，注释率为37.42%。其中，17532个unigene在KOG数据库中可归类到25个功能家族；10108个unigene被KEGG数据库注释到128个代谢途径中；17485个unigene被GO数据库的生物过程、分子功能和细胞组成三大功能注释到47个条目中。777个unigene在高温胁迫前后差异表达显著，其中胁迫后上调表达为362个，下调表达为415个，获得响应热胁迫的基因24个，包括热激蛋白的基因1个，PSⅠ和PSⅡ相关的基因15个、叶绿体rbcL基因3个、PLD基因2个，CAT编码基因、GAPDH基因、CYP编码基因各1个。本研究从转录组水平分析了金线莲对热胁迫的响应，这些数据将为功能基因的鉴定奠定基础，为培育耐高温金线莲品系提供了参考依据。     

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

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

基于转录组测序分析的黄精种子休眠解除相关差异基因研究

旨在获得黄精转录组数据库并挖掘参与其种子发育和休眠解除相关基因,以休眠解除前后的黄精种胚为试材,利用新一代高通量测序手段对供试样品进行转录组测序,并进行系统的生物信息学分析。黄精种子休眠解除前后样品中共得到79716个差异表达基因,上调的表达基因有60074个,下调的表达基因有19642个。休眠解除前后的黄精种胚中共有130284个差异表达基因被GO功能注释到生物进程、分子功能和细胞组分3个大类56个亚类,注释的差异表达基因与代谢过程、生物调控、细胞组分合成和酶催化活性等密切相关。KEGG代谢通路结果表明,共有65038个差异表达基因,涉及138个代谢通路,主要参与碳代谢、次生代谢产物的生物合成和多糖的代谢。基于KEGG数据库中注释结果,共发现15条与黄精种胚休眠解除相关的代谢通路。黄精种子发育与休眠解除过程,大量的种胚形态建成、多糖分解及蛋白质合成差异基因参与表达,并涉及到多个代谢途径的相互作用,构成复杂的休眠解除调控网络。     

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

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行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基因下调表达。为枇杷提供了较为全面的基因信息和代谢途径数据,为干旱胁迫下枇杷分子调控机制的深入研究奠定了基础。     

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

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

梁山慈竹高通量转录组测序及差异表达基因分析

分析梁山慈竹及其体细胞突变体的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为梁山慈竹遗传改良提供理论依据。利用RNA-Seq技术进行转录组测序,对测序结果进行de novo拼接和功能注释;并对差异表达基因进行筛选及COG、GO、KEGG数据库中进行比对注释,此外,基于Swiss-Prot功能注释结果,分析纤维素和木质素相关功能基因的表达量差异。测序结果表明,共获得86 575 631条reads,de novo组装得到84 741条unigenes,共有49 829条被Nr、COG、GO、KEGG、Swiss-Prot注释。从梁山慈竹实生植株(对照)和体细胞突变体No.30这2个测序样本中,筛选出3 572条差异表达unigenes,757条差异表达unigenes在COG分类体系中具有详细的蛋白功能释义,2 213条差异表达unigenes在GO数据库具有功能定义,385条unigenes被注释到94条KEGG Pathways中。纤维素合成相关纤维素合酶、过氧化物酶、泛素连接酶和热休克蛋白在梁山慈竹体细胞突变体No.30中表达量升高,木质素合成相关MYB4、4-香豆酸CoA连接酶、肉桂醇脱氢酶、肉桂酰-CoA还原酶和漆酶在突变体中表达量降低。提供了全面的梁山慈竹转录组信息,获得了一批在梁山慈竹纤维素和木质素生物合成过程中有重要功能的基因序列。     

菠萝蜜叶绿素缺失突变体茎次生生长转录组分析

对于木本植物而言茎的次生生长是植株发育和木材形成的基础,本研究以菠萝蜜(Artocarpus heterophyllus)正常植株为对照(control check, CK),将苗龄40 d的叶绿素缺失突变体(chlorophyll deficient mutant,CDM)和CK次生生长茎作为材料进行转录组测序,测序数据利用Trinity软件进行de novo组装共得到Unigenes 265 709个,去除低表达量(≤0.5)的转录本后筛选得到差异表达基因(differentially expressed genes,DEGs) 43 881个。利用加权基因共表达网络分析(WGCNA)找到与表型相关度最高的模块。将此模块中上下调的基因分别注释到KEGG和GO数据库,KEGG分类结果显示DEGs主要注释到代谢条目,GO分类结果中DEGs主要注释到与膜、细胞以及催化活性相关的过程。同时上调的基因在对活性氧的反应、抗坏血酸和醛酸的代谢、谷胱甘肽代谢等与胁迫相关的通路富集。这些结果可能暗示CDM茎的次生生长受到代谢和细胞过程的调控,生物或非生物因素引起的胁迫可能也参与到影响茎次生生长的机制当中。木本植物茎的次生生长涉及各种代谢通路和信号传导等复杂过程。通过转录组分析本研究为进一步了解白化病对木本植物茎次生生长的影响提供基础。     

姜黄根茎的高通量转录组测序与生物信息学分析

为获得姜黄(Curcuma longa)的转录组特征信息,本研究采用Illumina HiSeqΧTen高通量测序平台对姜黄根茎进行高通量转录组测序并进行系统的生物信息学分析。共获得7.18Gb Clean数据,组装了50194条unigenes,平均长度961.3 bp,N50为1 339 bp。数据库比对显示,姜黄根茎转录组unigenes在NR、Swiss-Prot、KEGG、KOG、eggNOG、GO、Pfam数据库中分别注释到38 802条(77.30%)、27 869条(55.52%)、14 725条(29.34%)、22 225条(44.28%)、37 317条(74.35%)、25 863条(51.53%)、26 137条(52.07%)。注释结果显示,姜黄与野生型马来西亚蕉的同源序列最多,unigenes在GO数据库中注释到参与生物过程、细胞组分和分子功能3个大类50小类,KOG功能分类获得25个不同的功能群,涉及128个KEGG代谢通路,其中包括21个次生代谢通路。在植物抗性基因(PRG)数据库中分别注释到3 718条unigenes;借助MISA软件发现7 183...     

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

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总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细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

转录组解析外源ABA对玉米脱水速率的影响

为了探明迪卡517、郑单1002喷施外源ABA后脱水速率变化的分子机制,发掘影响脱水速率的关键差异表达基因及代谢通路,利用Illumina HiSeq~(TM)2500测序仪分别对喷施ABA前后迪卡517(脱水速率较快)和郑单1002(脱水速率较慢)的穗位叶进行高通量转录组测序。原始序列经过质量控制、基因组比对、测序饱和度、基因覆盖度及冗余序列分析后进行基因功能注释。结果显示,迪卡517外源ABA处理与正常生长条件对比组检测到1 732个差异表达基因,其中1 251个为上调表达基因,481个为下调表达基因;郑单1002外源ABA处理与正常生长条件对比组检测到52个差异表达基因,其中24个为上调表达基因;迪卡517外源ABA处理与郑单1002外源ABA处理对比组检测到3 867个差异表达基因,其中1 946个为上调表达基因。不同对比组中筛选到差异表达基因GO分类情况、COG基因产物分类、KEGG代谢通路分析不同。转录因子分析共获得多个重要的转录因子家族,主要有AGC蛋白激酶家族、AP2/ERF转录因子家族、ARID转录因子、AUX/IAA转录因子家族、Alfin-like转录因子、Aur转录因子家族、B3转录因子家族、BBR-BPC转录因子家族、BES1转录因子、BUB转录因子、C2C2-CO-like转录因子等。Zm00001d035000、Zm00001d042063、Zm00001d045314等共同差异表达基因均在迪卡517和郑单1002中表达。对所有检测到的差异表达基因GO分类分析获得73个差异表达基因与水分代谢相关,推测这些基因作为外源ABA调控的下游基因起作用。     

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

为了探明玉米重要自交系花粒期高温胁迫下转录组基因的表达差异,发掘玉米响应高温胁迫的关键基因和蛋白质,明确高温胁迫引起玉米减产的机理,从而利用分子生物学技术手段提高玉米耐高温胁迫性。首先,利用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)在响应玉米花粒期高温胁迫过程中起着重要作用。     

基于高通量测序的南酸枣转录组分析

以南酸枣雌株叶片和成熟果实为材料,利用Illumina Hi Seq TM 4 000测序平台对其转录组进行测序,共获得14.15 G有效数据,通过序列拼接组装得到46 936条Unigene,平均长度为1 287 bp。36 299条Unigene(77.33%)在7大数据库(NR,NT,KO,Swiss-Prot,PFAM,GO和KOG)中得到注释,其中与GO数据库比对上的26 014条Unigene可分为生物过程、细胞组分和分子功能3大类53分支,注释到KOG数据库中的8 459条Unigene依据功能可分为25类。与KEGG数据库比对,6 225条Unigene分属5大类246条代谢通路中。南酸枣叶片和果实转录组中共发现5 631个差异表达基因,包括1 930个上调基因和3 701个下调基因。本研究获得的转录组数据将有助于开展南酸枣功能基因挖掘与利用、分子辅助育种和其种质资源遗传改良等方面的研究。     

百合未授粉成熟雌蕊转录组测序研究

本研究旨在大量发掘百合雌蕊特异表达基因,为研究百合杂交障碍的分子机制提供依据。利用Illumina高通量测序平台,进行百合未授粉成熟雌蕊转录组测序和数据组装,并对获得的Unigene进行功能注释、分类和代谢通路分析。“干柱头”雌蕊LoP1和“湿柱头”雌蕊LoP2分别获得40792条和39708条Unigene。2个样品的Unigene通过序列拼接,去冗余处理,以及同源转录本聚类,得到44175条All-Unigene。基于NR、NT、Swiss-Prot、GO、KEGG和COG等6个数据库进行相似性比对（E值≤10-5）,依次分别有30343、21647、21281、19274、12964和22227个基因被注释。KEGG分析中,19274条AllUnigene被注释到128个代谢通路上。基因分析显示,百合未授粉成熟雌蕊已为授粉受精做好了物质、能量、信号转导以及抗病原等准备工作。     

基于高通量测序的露地菊(Chysanthemum×grandiflora)盐胁迫转录组分析

为揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,本研究以盐胁迫处理的露地菊及其对照为材料,使用Illumina Hiseq2500高通量测序平台对转录组进行测序,分别获得了60370448和71415448条Clean reads,通过序列拼接组装得到45591条Unigene,平均长度724 bp。有37675条Unigene在七大数据库(COG,GO,KEGG,KOG,Pfam,Swiss-Prot,NR)中得到注释。通过比对露地菊盐胁迫处理组和对照组样品间Unigene的表达量及在各数据库中的注释情况,统计得到:有4143条差异表达基因获得注释;有2441条差异表达基因在GO数据库中获得功能注释;注释到COG数据库中的2281条差异表达基因依据功能可分为25类;有1062条差异基因映射到KEGGPathway数据库中,涉及了199个代谢通路,包括核糖体途径、植物激素信号传导途径、淀粉蔗糖代谢、碳代谢、氨基酸的生物合成等。本研究获得的转录组数据将有助于揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,及相关抗性基因的挖掘和分子辅助育种等方面的研究。     

水稻无侧根突变体RM109及其野生型对照的转录组分析

侧根是水稻根系的重要组成部分,具有重要的意义。本研究通过Illumina Hiseq2000测序平台,对无侧根突变体RM109与其野生型对照根系进行转录组测序,进而比较分析差异表达基因,初步揭示影响水稻侧根发育的相关基因。结果表明无侧根突变体与野生型对照根系相比出现了2 327个差异基因,其中上调表达的有1 131个,下调表达的有1 196个,有36个基因表达差异倍数高达100倍以上(|log2Ratio|≥6.7)。通过Gene ontology (GO)数据库、KEGG pathway数据库比对分析注释差异表达基因的功能和参与的分子调控途径。GO分析结果表明被注释的差异基因划分成50个功能类别,KEGG分析共有352个显著差异基因可以详细注释到149条分类代谢途径中,共有36个表达差异基因富集于植物激素信号转导途径,36个差异表达基因中有2个是与植物根系生长发育相关的Aux/IAA基因,6个植物生长素响应SAUR基因。本研究为进一步为水稻侧根发育分子机理研究提供基础,为水稻根系育种提供理论参考。     

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

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