圆口铜鱼亲本与子代肝脏转录组特征及对比分析

通过Illumina HiseqTM 2000测序平台首次开展了圆口铜鱼(Coreius guichenoti)亲本与子代肝脏转录组测序并对比分析其测序结果。对转录组进行拼接和组装, 共获得80688个Unigene, Unigene的长度主要分布在401—600 bp, 占总Unigene的43.56%。与Nr、GO、COG、KEGG等公共数据库比对并进行功能注释, 经分析圆口铜鱼亲本与子代差异表达基因共2701个, 其中上调1240个, 下调1461个。差异基因表达模式聚类热图显示, 同一样本不同重复基因表达相似。将差异基因映射到代谢通路KEGG数据库进行富集分析, 并对37个KEGG显著富集的代谢通路(P<0.05)构建其基因共表达网络, 结果显示丙酸盐代谢、丙酮酸代谢、原核生物固碳途径、乙醛和二羧酸代谢、脂肪酸代谢、萜类骨架生物合成是KEGG差异表达基因的核心代谢途径, 且这些基因在圆口铜鱼子代中表达量均显著上调。该结果丰富了圆口铜鱼基因组数据, 并首次从分子水平比较了圆口铜鱼亲本和子代的代谢差异, 同时解释了在同一循环水系统中子代不易患病的现象。     

禽致病性大肠杆菌pagP基因缺失株相关毒力基因的RNA-Seq分析

本研究为筛选禽致病性大肠杆菌pagP基因缺失后表达量变化的相关毒力基因,对禽致病性大肠杆菌AE17及pagP基因缺失株进行RNA-Seq测序,从中筛选出与毒力相关的差异基因并对其表达量的变化进行GO和KEGG分析。以|Fold-change|≥2且FDR≤0.05为条件从测序结果中筛选得到372个差异表达基因,其中339个基因表达上调,33个基因表达下调。GO功能分类结果显示差异基因主要集中在物质跨膜运输、转Class录调控以及生物膜的形成等功能注释;KEGG数据库显示,差异表达基因主要富集在鞭毛装配、双组分系统、脂多糖的生物合成以及ABC转运蛋白等通路。pagP基因缺失会引起相关毒力基因表达量的变化,根据RNA-Seq测序结果筛选出7个对APEC具有重要作用的毒力基因,为深入研究APEC的致病机理提供依据。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

基于RNA-Seq分析MpigE在红曲霉色素生物合成中的转录调控

【目的】分析MpigE的缺失在转录水平对红曲色素产生的影响。【方法】对实验室保存的野生型紫色红曲霉(Monascus purpureus)Mp-21和△MpigE菌株进行高通量转录组测序、注释、表达差异基因功能富集分析和基因通路富集分析,在转录水平揭示MpigE缺失后红曲霉色素产量变化的原因。【结果】通过RNA-Seq测序,每个样品获得7.5–8.5 Gb的原始数据,经过拼接后得到7219个转录本(Unigenes),其中成功注释的为5692个。差异基因表达富集分析发现基因缺失菌株△MpigE相较于野生型菌株Mp-21上调差异基因达到199个,下调差异基因为293个。【结论】MpigE的缺失能够促进红曲霉中中央碳代谢和乙酰辅酶A代谢相关基因的表达以此影响色素生物合成。     

基于转录组分析铜绿假单胞菌DN1降解荧蒽特性

[背景]铜绿假单胞菌DN1是一株从石油污染土壤中分离筛选到的具有广谱降解功能的菌株。[目的]深入了解荧蒽胁迫条件下铜绿假单胞菌DN1降解污染物过程中重要的降解相关基因信息。[方法]通过高通量测序技术对铜绿假单胞菌DN1进行转录组测序,对其所有的转录本进行KEGG (kyoto encyclopedia of genes and genomes)分类和Pathway注释、GO (gene ontology)分类和富集分析。[结果]转录组测序显示:与对照组相比,荧蒽诱导组检测到6 189个基因,其中1 919个基因上调表达,1 603个基因下调表达。KEGG注释分析显示差异上调表达基因匹配到了112个KEGG代谢途径,注释到"代谢途径"的1 408个基因(约占总差异基因的73.4%)中有317个基因参与了碳氢化合物代谢及含有苯环结构的异源生物质的生物降解,占"代谢途径"的16.53%,暗示了菌株DN1降解荧蒽可能与这些途径有密切关系。另外,主要代谢途径中的差异表达基因主要集中在ABC转运系统、氨基酸生物合成、双组分系统及碳代谢,这些途径大多数参与了底物的识别转运、信号转导及基因表达调控。[结论]进一步拓展了铜绿假单胞菌DN1在荧蒽胁迫条件下的代谢途径和逆境反应,也为微生物修复环境污染物研究夯实了理论基础。     

侧孢短芽孢杆菌S2-31拮抗下细辛叶枯病菌转录组差异表达分析

通过分析侧孢短芽孢杆菌拮抗作用下叶枯病菌的转录组学特征,研究差异表达基因(DEGs)和代谢通路的富集情况,初步探索侧孢短芽孢杆菌拮抗叶枯病菌的分子机制。首先利用S2-31与叶枯病菌的对峙培养观察其拮抗作用,然后利用转录组测序探究侧孢短芽孢杆菌拮抗下和正常生长下的叶枯病菌的基因表达水平差异,并进行RT-qPCR验证,最后,利用相关数据库对DEGs进行注释和富集分析。对峙培养后,叶枯病菌菌丝出现皱缩、扭曲等现象。测序后共得到3681个DEGs,其中2224个相对上调,1457个相对下调。GO功能注释分析表明,上调和下调表达的DEGs均注释到生物过程8个亚群、细胞组分8个亚群和分子功能4个亚群。KEGG富集分析中,共有732个unigenes定位到115条生物学通路中,其中富集程度相对较高的通路有氨基糖和核苷糖代谢、糖基磷脂酰肌醇(GPI)锚定生物合成和过氧化物酶体;涉及差异基因较多的通路有淀粉和蔗糖代谢、剪接体和胞吞作用;上调表达的DEGs主要富集在代谢途径,下调表达的DEGs主要富集在遗传信息处理和细胞过程途径。从显著富集的通路中筛选出与菌丝生长发育相关的差异表达基因,结果表明侧孢短芽孢杆菌主要抑制病原菌菌丝过氧化物酶体的形成、胞吞过程和GPI锚定的合成等过程。差异表达基因的RT-qPCR验证结果与转录组测序的结果一致。在侧孢短芽孢杆菌S2-31的拮抗下抑制了叶枯病菌的生长,其转录组特征也发生显著变化,差异基因主要涉及代谢、细胞过程以及遗传信息处理等途径的相关通路。     

基于RNA-seq的无量山乌骨鸡(Gallus gallus)肝脏脂代谢相关miRNAs的筛选

为研究无量山乌骨鸡(Gallus gallus)肝组织脂代谢相关miRNA (microRNA)在不同发育阶段的表达特征,本研究采集出壳当日(D1)和168日龄(D168)母鸡肝组织样品作为试验材料,利用DNBSEQ平台进行测序,采用DEGseq筛选差异表达的miRNA及其靶基因,随机选取9个差异表达miRNA进行RT-qPCR验证,KEGG通路分类筛选出脂代谢相关miRNA并进行聚类分析,预测脂代谢相关miRNA靶基因并进行GO和KEGG通路功能富集,构建脂代谢相关miRNA和靶基因关联网络。分析结果表明,筛选出106个差异表达miRNAs,包括54个上调miRNA和52个下调miRNA;聚类得到41个脂代谢相关的miRNAs;预测到38个靶基因,对靶基因的功能注释确定主要富集于甘油磷脂代谢、脂肪酸代谢和鞘脂代谢等脂质代谢相关通路,novel-gga-miR2311-5p-DGKZ、 novel-gga-miR2047-3p-ACACA、 novel-gga-miR866-3p-DGKH是脂代谢相关候选miRNA-mRNA关系对。研究提示无量山乌骨鸡肝组织miRNA在不同发育阶段的表...     

基于RNA-seq分析美洲大蠊提取物对结直肠癌细胞信号通路的影响

美洲大蠊Periplaneta americana提取物能够抑制多种肿瘤细胞生长,甚至能使肿瘤细胞发生凋亡,但是其对结直肠癌细胞信号通路的影响尚不清楚。本文基于RNA-seq技术初步分析了美洲大蠊提取物联合顺铂处理结直肠癌细胞后,细胞内信号通路的变化。结果表明:30μM顺铂和0.6%康复新液联合处理组与对照组相比,共有1 901个差异基因,其中1 555个基因表达上调,346个基因表达下调;30μM顺铂和0.8%精粉酵解液联合处理组与对照组比,共有2 587个差异基因,其中2 183个基因表达上调,404个基因表达下调;30μM顺铂和100μg·m L-1精粉联合处理组与对照组相比,共有1 488个差异基因,其中1 164个基因表达上调,324个基因表达下调。用WEGO和DAVID进行差异基因的GO和KEGG富集分析发现,美洲大蠊提取物联合顺铂处理组与对照组相比,表达上调的差异基因主要富集在p53、细胞粘附分子、MAPK、细胞凋亡等信号通路;表达下调的差异基因主要富集在氨基酰-tRNA生物合成、氨基酸生物合成、抗生素生物合成、一碳代谢等信号通路。     

基于高通量测序技术分析家兔感染猪瘟兔化弱毒株后脾脏的转录组学变化

本研究基于高通量测序技术对感染猪瘟兔化弱毒株(C株)后0h,12h,24h,30h,36h,48h的家兔脾脏组织进行测序,以家兔基因组作为参考,筛选感染后各时间点和对照组(0h)之间的差异基因。利用Uniprot和NCBI数据库查找各组中变化水平显著的前10个差异基因的生物学功能。结果发现,感染后12h,24h,30h的前10个差异基因只有少数与炎症、免疫、凋亡等相关;感染后36h,48h的前10个差异基因完全相同,其中B2M、RLA-DR-ALPHA、CD74和IGJ参与抗病毒过程。GO功能注释结果显示:感染后各时间点的差异基因都和免疫、代谢、调控途径紧密相关。KEGG通路富集分析发现,感染后24h的差异基因富集到RIG-I样受体信号通路,感染后30h的差异基因富集到黏着斑和细胞外基质-受体相互作用通路,感染后36h和48h的差异基因共同富集到的通路有20个,其中蛋白酶体,溶酶体,核糖体,趋化因子信号通路,B细胞受体信号通路,抗原加工提呈通路和FcγR介导的吞噬通路与抗病毒相关。本研究建立的家兔感染C株后脾脏组织的差异表达基因数据库,将为进一步阐述家兔适应C株的分子机制奠定基础。     

多年宿根甘蔗的转录组分析及差异基因挖掘

甘蔗(Saccharum hybirds)宿根性直接关系到甘蔗生产成本和种植效益,品种、环境和栽培措施均会影响甘蔗宿根能力,但品种的种性是影响宿根性最关键因素。国内外从分子生物学层面解释甘蔗宿根性差异的文献鲜见报道,从转录水平分析不同宿根年限GR2号和ROC22号基因表达的差异。结果表明：转录组测序共得到100 558条转录本和25 582条Unigene。获得53 790条Unigene的注释结果,分别在NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam数据库进行比对。GO功能注释共分成三大类,51小类,6年宿根蔗注释到1 029个差异基因,3年宿根蔗注释到3 391个差异基因。主要KEGG代谢通路有8条,分别涉及植物激素信号转导,淀粉和蔗糖代谢,甘氨酸、丝氨酸和苏氨酸代谢,苯丙素生物合成,同源重组,DNA复制,错配修复及植物病原体相互作用。筛选出脱落酸(ABA)相关差异基因6个,分别为脱落酸不敏感蛋白2基因(ABI2)、bZIP转录因子超家族蛋白、碱性亮氨酸拉链型转录因子基因(ABI5)、aba响应元件结合因子1基因(ABF1)、G-box结合因子基因(GBFs)...     

功能雄性不育茄子差异基因及代谢途径分析

以茄子功能雄性不育系 (S13) 和可育系 (F142) 为材料,对花蕾不同发育时期的细胞学特征进行观察,发现S13花药中环形细胞簇的解体时期比F142延迟2 d,开花当天F142的裂口组织及内绒毡层的细胞解体,而S13中无此现象发生。转录组测序分析表明,F142和S13开花前8天、5天和开花当天的花药共有1 436个差异表达基因 (DEGs) (651个表达上调,785个表达下调)。GO显著性分析表明差异基因主要富集在生物学过程中涉及单细胞生物过程、代谢过程和细胞过程的单基因簇 (Unigene) 较多,分子功能中涉及较多的有催化活性和结合功能;KEGG注释发现,差异基因主要富集在次生代谢产物的生物合成、代谢途径、内质网蛋白质加工、氨基酸的生物合成、碳代谢和植物激素信号转导。选取16 465个基因进行加权基因共表达网络构建分析 (Weighted gene co-expression network analysis,WGCNA),共鉴定到15个基因共表达模块,其中3个为S13在3个花蕾发育时期高度相关特异性模块 (Plum2、Royalblue和Bisque4模块);KEGG富集表明,特异性模块可以富集到苯丙烷类生物合成、光合作用、卟啉与叶绿素代谢、α-亚麻酸代谢、多糖的生物合成和代谢、脂肪酸降解以及戊糖与葡萄糖醛酸的相互转化等相关通路,这些基因可能在S13花蕾发育过程中发挥重要作用。研究结果为进一步探索茄子花药开裂机制提供了一定的参考。     

草菇耐低温突变菌株Vtlt-1和原始菌株V23转录组差异的研究

为了培育草菇耐低温菌株与解析其耐低温的分子机制,采用紫外诱变的方法,选育出耐低温草菇菌株Vtlt-1,并利用表达谱芯片技术,比较突变菌株Vtlt-1与原始菌株V23的表达差异基因,筛选差异显著基因后利用GO（gene ontology）功能注释和KEGG（kyoto encyclopedia of genes and genomes）通路富集等方法分析,结果发现与V23相比Vtlt-1表达显著差异基因共有1 600个,其中704个基因上调表达,896个基因下调表达。针对差异基因的GO功能分类结果发现：生物学过程方面,差异基因主要分布在金属离子结合、氧化还原过程、碳水化合物的代谢与核酸的结合。细胞组分方面主要与细胞核相关;分子功能中,氧化还原活性以及解旋酶活性,依赖于ATP 的解旋酶活性、DNA指导的RNA聚合酶活性。KEGG注释结果发现差异表达基因主要富集在氨基酸与氮类物质的代谢、脂肪酸与生物碱的合成这两方面,此外还富集到核糖体的生物合成、细胞色素P450、RNA聚合酶、硫和氨基酸的代谢、核酸的修复等通路。这些结果为解析草菇耐低温机制提供分子依据和理论基础。     

RNA-seq用于获得共轭亚油酸对贵妃鸡肌内脂肪代谢调控的关键基因

通过转录物组测序获得在贵妃鸡基础日粮中添加共轭亚油酸(CLA)对肌内脂肪代谢的差异表达基因,经生物信息学分析获得相关的信号通路及可能发挥重要作用的候选基因,为CLA对肌内脂肪沉积的分子机制奠定基础。本研究选用55日龄健康的贵妃鸡为试验动物,在基础日粮中添加CLA 0%、1%和2%,预饲期1周,正饲期6周。屠宰采集胸肌组织进行转录物组测序,对测序数据进行差异表达分析,差异表达基因GO功能和差异表达基因KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qRT-PCR对差异表达基因进行验证。结果显示,共获得1 065个差异表达基因,其中上调基因703个,下调基因362个。GO富集结果显示,差异表达基因主要富集在生物过程的细胞过程、单一生物过程、生物调节和代谢过程。KEGG信号通路富集显示,差异表达基因显著富集在黏着斑、不饱和脂肪酸生物合成、脂肪酸生物合成和类固醇生物合成等信号通路中,发现11个主要与肌内脂肪代谢相关的候选基因,分别是FADS1、FADS2、ELOVL5、ACOX2、SLC27A1、FABP5、LPL、LOC107050163、ENSGALG00000030996、ENSGALG00000005043和ENSGALG00000048882。并随机选取6个基因进行qRT-PCR验证,其相对表达量变化趋势与测序结果一致。本研究筛选到CLA影响贵妃鸡胸肌脂类代谢相关的差异表达基因,并对11个主要参与脂肪代谢相关的基因进行分析,为揭示CLA调控肌内脂肪沉积的分子机制奠定基础。     

金钱鱼下丘脑转录组分析及其对雌激素在体注射的响应

为探究XX雌性和XY雄性金钱鱼(Scatophagus argus)下丘脑中的差异表达基因及雌激素(17β-Estradiol,E2)在体注射对XY雄鱼下丘脑中基因表达的影响,开展了转录组测序分析,包括测序数据质控、基因功能注释,差异基因筛选、鉴定和差异基因功能富集分析等,并利用实时定量PCR(Real-time quantitative PCR, qPCR)检测了金钱鱼下丘脑中相关基因的表达。金钱鱼下丘脑转录组共测得Clean reads 275833710,测序质量30(Q30)大于95%, GC含量值大于48%。其中, Ctrl-XX-H vs. Ctrl-XY-H组,共筛选和鉴定了91个差异表达基因,包括36个上调基因和55个下调基因。在Ctrl-XY-H vs. E2-XY-H组,筛选和鉴定了28个差异表达基因,包括11个上调基因和17个下调基因。GO和KEGG富集分析发现,差异表达基因主要显著富集在细胞、单生物过程,膜、膜组分,结合等生物学功能及泛醌和其他萜类-醌生物合成及类固醇激素生物合成通路等。转录组和qPCR结果表明, XX和XY金钱鱼下丘脑中prl、ccna1和hs...     

黑果枸杞果实成熟发育过程表达谱差异分析

对黑果枸杞(Lycium ruthenicum Murr.)3个不同发育时期的果实进行转录组测序,分析对比黑果枸杞果实不同发育时期相关基因表达谱的变化。提取果实RNA进行Illumina高通量测序,利用GO、KEGG等公共数据库进行功能注释、分类,利用数字基因表达谱技术对比分析。结果显示,KEGG pathway富集分析表明亮氨酸生物合成途径在变色期对比青果期有7个差异表达基因,黑熟期对比变色期有35个差异表达基因,且全部为上调表达。在糖酵解途径中3个不可逆的反应步骤限速酶基因在黑熟期均为上调表达。对黑果枸杞亮氨酸合成和糖酵解涉及到的基因进行了功能和代谢通路的分析,为黑果枸杞种质资源探索提供理论基础。     

RNA-Seq探索低水平利奈唑胺耐药粪肠球菌机制

为了探索低水平利奈唑胺耐药粪肠球菌机制,本研究采用Illumina HiSeq 4000高通量转录组测序技术对低水平利奈唑胺耐药粪肠球菌(P10748 MIC:8 mg/L)和利奈唑胺敏感粪肠球菌(3138 MIC:2 mg/L)进行测序及生物信息学分析,以标准菌株ATCC29212作为质量评估参考。利用GO和KEGG数据库对差异表达基因进行注释与富集分析,并利用荧光定量PCR对测序中的差异表达基因进行验证。结果显示,测序共获得了3.57 Gb有效数据,通过De novo拼接获得1 920条unigenes,总长度为2 122 210 bp,平均长度为1 105 bp。差异基因模式聚类分析显示共有150个显著性差异表达基因(FDR≤0.001,|log2 Ratio|≥1),其中141个上调,9个下调。其中生物膜形成和外排泵相关基因esp、optrA、fexA在耐药菌中显著上调。GO功能注释结果显示,催化活性、代谢过程和细胞是注释最多的条目;Pathway富集分析发现,肽聚糖生物合成、缬氨酸亮氨酸和异亮氨酸的降解和硫代谢为显著性富集通路(p<0.05)。荧光定量PCR结果与测序结果基本一致。推测膜转运蛋白和生物膜形成可能在耐药机制中具有重要作用,为低水平利奈唑胺耐药肠球菌机制提供了可参考的耐药靶标。     

濒危植物秦岭石蝴蝶花瓣数量变异机理研究

为了探明秦岭石蝴蝶花瓣数量变异原因,该研究采用Illumina HiSeq 2500高通量测序技术对秦岭石蝴蝶两种花型发育的早期和晚期进行转录组测序,挖掘参与其花发育相关的差异基因,并探讨花器变异的可能机制。结果显示：(1)与NR数据库进行比对,共有52 677个Unigene注释到NR库,占Unigene总数的46.25%,与旋蒴苣苔(Dorcoceras hygrometricum)的序列同源性最高(54.29%)。(2)GO富集分析结果显示,正常2 3型和变异2 4型差异基因富集最显著的GO条目为：细胞组分类的细胞膜、分子功能类的反转运蛋白活性和酶抑制剂活性、生物过程类的跨膜转运和催化活性的负调控等;KEGG富集分析结果显示,正常2 3型和变异2 4型差异基因富集最显著的KEGG通路为：植物激素信号转导、脂肪酸延长、戊糖和葡萄糖醛酸酯的相互转化、苯丙烷生物合成、硫代谢、类黄酮生物合成、玉米素的生物合成等途径。(3)对表达差异基因进行筛选并进一步对4个比较组进行交叉比对分析,确定了6个可能与秦岭石蝴蝶花器官发育相关的基因,分别为PqMIF2、PqMYB340、PqMYB305、PqGATA12、PqCCD4和PqZBED;qRT PCR验证发现,其表达趋势和转录组测序分析结果一致。该研究结果为秦岭石蝴蝶的花器官发育和系统进化及其濒危机制方面研究提供了重要的信息。     

低氧胁迫下唇?及杂交F1代肝脏转录组学分析

为探讨低氧胁迫对唇?Hemibarbus labeo和杂交?基因表达水平的影响，基于Illumina HiSeq高通量测序平台，对低氧条件下唇?和杂交?肝脏组织进行转录组测序并开展生物信息学分析。结果显示，唇?和杂交?分别产生了126 646 494个和130 613 321个clean reads，组装后共获得了48 027个unigenes，分别在NR、SwissProt、GO、COG、KOG、eggNOG4.5、KEGG数据库中获得32 253个注释；低氧胁迫时，分别有3 589个和3 177个差异表达基因（DEGs）上调，4 075个和5 473个下调；GO功能富集分析发现，DEGs与细胞结合、催化活性和代谢进程密切相关，KEGG通路富集分析发现，DEGs富集在一些与代谢和细胞生长增殖、分化、凋亡相关的信号通路中。研究结果为深入开展唇?和杂交?溶解氧胁迫的调控机制和养殖育种提供了参考。     

基于转录组测序的羽衣甘蓝叶色相关基因分析

叶色是羽衣甘蓝重要的观赏性状之一。本研究采用Illumina Hi-Seq2500高通量测序技术,对基因型纯合的紫叶和白叶羽衣甘蓝叶片进行转录组测序,筛选差异基因并与GO和KEGG数据库比对进行注释分析,分析羽衣甘蓝叶色形成相关基因。结果显示,获得高质量短读序共104 608 770条,筛选出紫叶相对白叶的差异表达基因1 993个,其中上调表达基因1 094个,下调表达基因899个。根据GO功能分类可分为生物过程、细胞组分和分子功能3大类64功能组。根据KEGG代谢通路分析可以分为171类,在叶色相关的类黄酮生物合成途径中黄酮醇合成酶(FLS)、二氢黄酮醇4-还原酶(DFR)上调以及类胡萝卜素生物合成途径中的类胡萝卜素β-环化酶上调与紫叶形成关系密切。本研究丰富了羽衣甘蓝的转录组信息,获得了一些差异表达基因,为进一步研究羽衣甘蓝叶色形成的遗传机制提供了有价值的信息。