利用WGCNA鉴定玉米非生物胁迫相关基因共表达网络

加权基因共表达网络分析(weighted gene co-expression network analysis, WGCNA)是一种经典的系统生物学分析方法,可用来鉴定协同表达的基因模块,探索模块与目标性状的生物学相关性,并挖掘模块网络中的核心基因。本文收集了低温胁迫、高温胁迫、干旱胁迫和盐胁迫处理下玉米(Zea mays L.)根、茎、叶等组织的58份转录组数据,利用WGCNA方法鉴定玉米非生物胁迫的基因共表达网络模块。过滤转录组数据中12,552个低表达基因后,利用余下27,204个高表达基因构建共表达网络,分析得到25个模块。根据玉米中已报道非生物胁迫相关基因与差异表达基因在模块中的分布,筛选出与低温胁迫、高温胁迫、干旱胁迫和盐胁迫最相关的mediumpurple4、ivory、coral2、darkseagreen4模块和响应多种胁迫的green模块。随后对这5个模块的基因进行GO富集分析,发现在这些模块内与非生物胁迫相关基因的在非生物胁迫调控相关功能显著富集,如胁迫响应、过氧化物酶活性等。对这5个模块作相关性分析,预测出Zm00001eb072870、Zm00001eb32...     

棉属D基因组3个野生棉种Bet v1基因鉴定及功能分析

【目的】对棉花D基因组中Bet v 1基因家族进行分析,比较其在不同抗性棉种间的表达模式差异,为深入研究Bet v 1基因在棉花抗黄萎病中的作用提供理论依据。【方法】通过生物信息学方法对D基因组中Bet v 1基因进行鉴定。通过雷蒙德氏棉(Gossypium raimondii,D5)、三裂棉(G. trilobum,D8)和瑟伯氏棉(G.thurberi,D1)转录组和实时荧光定量聚合酶链式反应(qRT-PCR)分析Bet v 1基因在黄萎病菌处理下的表达模式。利用病毒诱导的基因沉默技术(VIGS)对Bet v 1基因进行功能鉴定。【结果】[棉花D基因组中包含59个成员,其中57个基因带有内含子,分布于8条染色体上,多数为亲水性蛋白并定位于细胞质。]黄萎病菌胁迫条件下3个野生棉种的Bet v 1基因表达量与其抗病水平一致。将不同表达水平的基因分为3组,其中第3组基因响应黄萎病菌侵染并在抗病棉种瑟伯氏棉中高表达,表明该组基因可能与黄萎病胁迫应答反应有关。从中筛选出高水平表达的Bet v 1基因,在陆地棉中沉默相应的同源基因,棉株感病加重,揭示该基因在棉花抵御黄萎病菌侵染过程中起正调控作用。【结论】响应黄萎病菌胁迫的Bet v 1基因在棉花抗黄萎病复杂的生物过程中至关重要。本研究结果为棉花Bet v 1家族基因的深入研究提供依据,为进一步解析棉花Bet v 1基因的功能奠定基础。     

棉花抗黄萎病相关基因的序列和表达分析

黄萎病是棉花生产的主要病害,克隆抗病基因是培育棉花抗黄萎病品种的关键。本文根据前期的定位结果,结合棉花基因组测序信息,提取定位区段内基因组序列,预测获得63个基因。Gene Ontology分析表明63个基因参加多种生物进程,其中6个基因参与植物抗逆进程。根据Gene Ontology分析和前人研究结果,选择15个基因进行序列分析。启动子序列分析表明启动子区域包含各种抗逆的调控元件,其中6个基因包含了W-box元件。对海7124和苏棉8号棉花幼苗进行黄萎病菌处理,取病菌处理后不同时期的根,选择14个基因进行表达分析,结果表明在黄萎病菌处理后,有8个基因表达出现变化,其中在海7124和苏棉8号之间表达差异最大的基因是跨膜蛋白(Transmembrane protein 214-a isoform 1)、细胞色素P450(Cytochrome p450)和Udp-糖基转移酶UGT89A2(Udp-glycosyl transferase 89a2-like)基因。本研究为抗病基因克隆提供了候选基因。     

陆地棉TLP基因家族的全基因组鉴定及表达分析

【目的】类甜蛋白(Thaumatin-like proteins,TLP)参与多种植物病原真菌侵染后的防御反应。本研究通过在全基因组分析TLP基因,为深入探究其在介导棉花抗黄萎病通路中的分子机理提供基础。【方法】分别利用生物信息学方法和荧光定量聚合酶链式反应技术对陆地棉中的TLP基因进行全基因组鉴定和响应棉花黄萎病菌侵染表达分析。【结果】共鉴定出88个TLP基因;通过系统发育分析和序列特征分析可将棉花TLP家族分为10类。TLP基因家族外显子数量介于1～5之间,内含子介于0～4之间,同组成员具有相似的基因结构。大部分TLP含有5个保守的motif,具有相同的基序组织模式,均为Motif 5_4_2_3_1。染色体定位显示87个TLP基因定位于陆地棉20条染色体上,其中A亚组和D亚组各分布有42和45个,其在染色体上的分布存在串联重复现象。响应病原菌侵染表达分析显示,黄萎病菌侵染可以诱导6个候选基因的表达,不同抗性品种相对表达量达到峰值的时间不同,且耐病品种(GZ-1)上述基因表达量比感病品种(86-1)有增高的趋势。【结论】该研究结果有助于了解棉花TLP基因家族的进化与功能。     

棉花三元杂种(HBT)衍生系应答黄萎病菌侵染反应

棉花黄萎病是危害棉花维管束的严重病害,从陆地棉、海岛棉及近缘野生种杂交后代中筛选抗黄萎病种质能有效拓宽陆地棉抗源基础。本研究通过人工接种强致病力、落叶型大丽轮枝菌(Verticillium dahliae Kleb.)菌株Vd991,对包括海岛棉(G.barbadense)、陆地棉(G.hirsutum)和瑟伯氏棉(G.thurberi Todro)三元杂种衍生系在内的42份棉花种质的抗病性进行了鉴定,并以棉花Gh ACTIN14(Gen Bank:AY305733)为内标,设计抗病相关基因Gb DIR1、Gb DIR2、β-1,3葡聚糖酶基因的q RT-PCR(quantitative Real-time PCR)特异引物,使用SYBR Green PCR试剂盒标记反应产物,在Applied Biosystems 7500 Real-Time PCR System上采用2-ΔΔCt法分析高抗黄萎病的三元杂种衍生系海陆野96-5受黄萎病菌侵染的调控表达特征。结果显示:筛选出的16份抗黄萎病种质中有2份高抗种质和12份抗病种质均为陆海瑟三元杂交种衍生系或衍生品种;q RT-PCR分析表明高抗病品系海陆野96-5受黄萎病菌侵染1 h时,Gb DIR1、Gb DIR2和β-1,3葡聚糖酶即开始上调表达,Gb DIR1到8 h时上调最高,达到处理前20倍,Gb DIR2则在24 h和96 h表达量较高,分别上调15倍和29倍,β-1,3葡聚糖酶在8 h和48 h时上调表达2倍以上。本研究结果为高抗黄萎病资源的利用奠定了基础。     

陆地棉PPO基因全基因组鉴定及对黄萎病菌的响应分析

【目的】多酚氧化酶(Polyphenol oxidases,PPO)广泛参与植物抵抗病虫害等生物逆境胁迫的过程。本研究旨在研究棉花中的PPO基因及其表达模式,验证多酚氧化酶与棉花抗黄萎病的关系。【方法】分析了黄萎病菌V991侵染下异源四倍体陆地棉TM-1根系中PPO酶活力变化,并利用TM-1的基因组数据库,鉴定相关基因,并进行生物信息学和表达分析。【结果】共筛选到13个候选GhPPO基因。这些基因都不含内含子,所编码的蛋白包含2个铜离子结合位点和PPO的保守序列(酪氨酸酶、DWL和KFDV保守结构域)。进化分析显示,陆地棉PPO基因的种内相似性大于种间相似性,并且相互之间形成了多对重复基因。GhPPO的表达量差异较大,少数基因具有组织特异表达的特性,多数基因在棉花不同组织中表达量都很低。实时荧光定量聚合酶链式反应分析结果表明:GhPPO6D在棉花根系中表达量较高,且在黄萎病菌V991侵染后上调。【结论】GhPPO6D可能参与了棉花与黄萎病菌的互作过程,与V991侵染后棉花根系PPO酶活力上升相关。     

基于WGCNA的马铃薯根系抗旱相关共表达模块鉴定和核心基因发掘

权重基因共表达网络分析(Weighted Gene Co-expression Network Analysis, WGCNA)是系统生物学皀一种研究方法,在多样本转录组数据中挖掘与目标性状相兲皀基因模块有较广泛皀应用。为了深入探究马铃薯应对干旱胁迫皀分子机制,本研究以国际马铃薯中心引迚栻培种C16 (CIP 397077.16)和C119 (CIP 398098.119)为试验材料,将其无菌组培苗利用甘露醇模拟干旱胁迫,处理0 h、2 h、6 h、12 h和24 h,取其根系迚行转录组测序,每样设3个生物学重复,共30个样本。基于以上转录组数据,利用WGCNA构建与抗逆生理性状相兲联皀权重基因共表达网络,得到15个与根系抗旱密切相兲皀基因共表达模块,幵从4个与目标性状兲联度最高皀模块中収掘到数个核心基因,功能注释表明其中大部分参与干旱胁迫调控通路。这些结果为迚一步研究马铃薯根系抗旱皀分子遗传机制提供了线索。     

棉花WRKY22基因分离及其对黄萎病的抗性分析

为了克隆棉花抗黄萎病相关基因,利用生物信息学、病毒诱导基因沉默技术(VIGS)、实时定量PCR(q PCR)和接菌分析来研究棉花WRKY基因对大丽轮枝菌的诱导响应和抗性。结果从棉花中筛选出8个与拟南芥直系同源的棉花抗病相关WRKY基因,这些直系同源WRKY基因表达受到大丽轮枝菌浸染诱导响应。其中有6个GhWRKY基因的表达均受到大丽轮枝菌诱导上调,而GhWRKY70和GhWRKY48的表达量随着不同的时间点呈现上调或下调。GhWRKY22等5个基因表达量在24,72 h分别表现为上调,呈现出双峰曲线。对GhWRKY22表达特征进一步分析,结果显示,GhWRKY22基因在茎里优势表达,同时表达受到大丽轮枝菌、水杨酸和茉莉酸激素的诱导。GhWRKY22基因沉默植株对大丽轮枝菌的敏感性增加,抗病标志基因(PR1、PR3、PR4、PR5、PAL和PDF1. 2)表达量也显著降低,揭示GhWRKY22是通过SA和JA信号途径参与棉花对大丽轮枝菌的响应过程。总之,棉花中8个WRKY基因参与棉花对黄萎病抗性调控,其中GhWRKY22基因正调控棉花的抗性,可作为棉花抗病育种的候选基因。     

黄萎病菌侵染对感病棉花品种叶片蛋白质组的影响

为进一步探索黄萎病菌与棉花的相互作用,以感黄萎病棉花品种冀棉11为实验材料,在接种大丽轮枝菌V991 3 d后提取叶片蛋白,运用双向电泳和质谱分析技术研究黄萎病菌侵染后叶片蛋白质组学的变化。分析发现,与对照相比差异表达量在2倍以上的蛋白点有22个,其中上调表达的蛋白点10个,下调表达的蛋白点12个。质谱鉴定发现,上调表达的蛋白包括ATP合成酶β亚基,Ru Bis CO活化酶1,Ru Bis CO活化酶α2以及果糖-1,6-二磷酸醛缩酶,下调表达的蛋白包括质体叶绿素AB结合蛋白,核酮糖二磷酸羧化酶大链以及核酮糖二磷酸羧化酶小链。分别对下调表达的蛋白rubisco小亚基基因(RBCS)和叶绿素AB结合蛋白基因(cab)、上调表达的蛋白Ru Bis CO活化酶1基因(RCA1)进行荧光定量PCR,发现接种黄萎病菌后3个基因在m RNA水平与蛋白水平上的变化是一致的。通过分析推测,黄萎病菌通过与参与光合能量代谢的蛋白相互作用,破坏植物的光合作用,近而引起棉花叶片的萎蔫与凋零。     

棉花GbSTK基因调控开花和黄萎病抗性的功能研究

棉花是重要的经济作物,实现既抗病又早熟是棉花重要的育种目标,但棉花抗病与早熟基因之间的关联性研究报道甚少。本课题组前期鉴定到一个响应黄萎病菌诱导的丝氨酸/苏氨酸蛋白激酶基因GbSTK,可以提高转基因拟南芥的黄萎病抗性。本研究发现多个抗病棉花品种中STK基因的表达量显著高于感病对照品种H208。在抗病品种农大601中沉默STK基因,沉默植株在接菌后20d,其黄萎病抗性显著降低,病情指数由27.5(耐病)升高到63.2(感病),表明该基因正向调控棉花黄萎病抗性。对转GbSTK基因拟南芥表型鉴定发现,其开花时莲座叶片数平均为7.5个,显著低于空载体对照植株(11.9个),转基因植株开花时间较对照平均提早5~7 d。进一步对拟南芥开花途径标志基因FT、SOC1和LFY的转录水平检测发现,转基因不同株系中FT和SOC1的表达水平显著高于对照植株,而LFY基因的表达受影响较小,表明STK基因可以影响开花关键基因FT和SOC1的表达,进而调控植株提早开花。基于酵母双杂交互作蛋白筛选试验,初步获得30个与GbSTK具有互作关系的候选蛋白,为进一步揭示GbSTK的分子调控网络奠定了基础。本研究首次鉴定出可同时提高黄萎病抗性和提早开花的一个功能基因,为棉花抗病早熟遗传改良提供了参考依据。     

陆地棉漆酶基因家族鉴定及在黄萎病菌胁迫下的表达分析

黄萎病是降低棉花产量和品质较为严重的维管束病害。棉花在抵抗病原菌的过程中,抗病基因起着尤为重要的作用。漆酶是一种多功能氧化酶,在植物木质素合成和提高植株抗逆性等方面发挥着重要作用。高质量版本的基因组是提升基因家族分析准确性的必要条件。本研究以目前最新的陆地棉TM-1高质量版本基因组为参考,通过生物信息学分析鉴定了陆地棉基因组中的Laccase(GhirLAC)基因家族,并对其进行了理化性质、基因结构、染色体定位以及在黄萎病胁迫下的表达模式分析。结果表明,陆地棉TM-1基因组中共包含83个GhirLAC家族成员,分布在24条染色体上,所有GhirLAC蛋白均定位在胞外,且具有相同/相似的保守基序。系统发育树分析显示,GhirLAC基因家族成员可分为7个亚组。结合黄萎病胁迫下棉花转录组数据,将GhirLAC基因家族各成员的表达分为3种模式,其中第1类和第2类GhirLAC基因在黄萎病菌胁迫下分别呈现有规律的表达量升高和降低,推测这些基因在棉花抗黄萎病反应中发挥重要功能。荧光定量PCR结果表明,GhirLAC02(GhLAC4)、GhirLAC38(GhLAC11)和GhirLAC20(GhLAC12)3个候选基因均受黄萎病菌诱导表达,其表达趋势与转录组数据相吻合。本研究为以后深入解析棉花GhirLAC基因的抗病功能及分子机制奠定基础。     

用全基因组关联作图和共表达网络分析鉴定油菜种子硫苷含量的候选基因

油菜籽饼粕是畜禽养殖中重要的蛋白原料, 但饼粕中的硫苷是一种抗营养物质, 食用过多会对禽畜产生毒害, 因此挖掘油菜籽粒硫苷含量的候选基因对油菜种子低硫苷育种具有重要现实意义。本研究连续4年种植1个含157份材料的油菜自然群体, 结合重测序数据对种子硫苷含量进行全基因组关联分析(GWAS), 并对15份低硫苷和15份高硫苷材料进行种子发育早期的转录组测序, 通过权重基因共表达网络分析(WGCNA)鉴定种子硫苷含量的候选基因。用GWAS共检测到45个与种子硫苷含量显著相关的SNP, 单个位点解释的表型变异为13.5%~23.3%, 主要分布在A09、C02和C09染色体的3个区间中, 覆盖5个已知的硫苷代谢基因。用WGCNA分析发现高、低硫苷材料之间的2275个差异表达基因, 可分为12个基因模块, 其中1个模块的基因显著富集在已知的硫苷生物合成途径, 对该模块内163个基因的权重分析得到13个候选基因。经检测, GWAS和WGCNA共得到的18个候选基因中, 有14个候选基因的表达量与种子硫苷含量显著相关(r = 0.376~0.638, P<0.05)。用两种方法鉴定到1个共同的候选基因BnaC02g41790D (基因名MAM1), 与该基因连锁的5个SNP构成5种单体型, 等位基因效应分析发现, 自然群体中63%的材料(99/157)为Hap 5, 平均硫苷含量为50.79 μmol g^-1, 与另外4种单体型(95.04~110.28 μmol g^-1)存在极显著差异(P<0.01)。本研究结合GWAS和WGCNA两种方法鉴定了油菜种子硫苷含量的候选基因, 可为复杂性状候选基因的筛选提供参考。     

棉花GbSTK基因调控开花和黄萎病抗性的功能研究

Cotton is an important cash crop, and achieving disease resistance as well as early maturity is an important breeding goal. However, the molecular relationship between disease resistance and early maturity of cotton is still a research gap. In previous study, we identified a GbSTK gene response to V. dahliae that could enhance Verticillium wilt resistance in transgenic Arabidopsis. In this study, we found that GbSTK was more highly expressed in resistance and/or tolerant varieties than in susceptible varieties. Knocking down of GbSTK expression in Nongda 601 showed that the silenced plants displayed serious disease symptoms, with disease index from 27.5 (tolerant) to 63.2 (susceptible), suggesting that GbSTK positively regulates cotton Verticillium wilt resistance. More interesting, we found that GbSTK could promote the development of Arabidopsis. At the flowering time point, the average number of rosette leaves in transgenic lines was 7.5, significantly less than the mock plants, and the transgenic plants bloomed five to seven days earlier than the control. The FT and SOC1 gene, regarded as marker genes in Arabidopsis flowering signaling pathway, maintained higher expression level in transgenic lines than in the mock, indicating that GbSTK could accelerate flowering via regulating FT and SOC1 gene expression. We also obtained 30 candidate proteins that could interacted with GbSTK via yeast double hybrid test, which helped to further revealing the molecular regulatory network of GbSTK. Taken together, we first identified a functional gene that could improve disease resistance and early flowering, which would provide a reference for genetic improvement of cotton disease resistant and early maturity.     

利用WGCNA鉴定玉米株高和穗位高基因共表达模块

株高和穗位高是玉米株型的重要影响因子,与产量性状紧密相关。加权基因共表达网络分析(weightedgene co-expression network analysis, WGCNA)是探索基因网络与特定性状间关联关系的重要方法,为株高和穗位高相关基因的挖掘提供新途径。本研究利用郑58、掖478、昌7-2和黄早四及其组配的杂交种郑单958、安玉5号、郑58/黄早四和掖478/黄早四,结合其在45,000株hm^–2和67,500株hm^–2条件下的转录组数据,采用WGCNA构建了2种密度条件下的共表达网络,分别得到24个和21个共表达模块,并鉴定到与株高和穗位高显著且高度相关(相关系数的绝对值>0.50)的共表达模块15个,其中两性状相同的模块共6个。基因功能富集分析结果表明,株高和穗位高共表达模块主要参与生长发育、光合作用、响应光刺激、植物激素、碳水化合物合成/代谢等重要活动。根据模块内基因的连接度,发现乙烯响应因子EREB14、硫胺素酶TENA2、磷酸甘油酸激酶PGK、谷胱甘肽转移酶GST2和琥珀酸脱氢酶SUDH7等是模块内的核心基因。...     

整合GWAS和WGCNA分析挖掘甘蓝型油菜黄籽微效作用位点

甘蓝型油菜是世界上最重要的油料作物之一, 黄籽是提高品质的重要育种目标。本研究以520份具有代表性的甘蓝型油菜品种(系)为材料, 结合种子发育过程中8个时期的转录组数据, 采取整合全基因组关联分析(GWAS)和权重基因共表达网络分析(WGCNA)的策略, 挖掘油菜黄籽性状微效作用位点, 2年共检测到199个SNP位点, 在SNP位点附近共挖掘出1826个名义候选基因。利用R语言中的WGCNA软件包构建了8个共表达模块, 基因功能富集分析显示, turquoise模块和blue模块与黄籽表型相关。苯丙烷代谢途径、类黄酮途径的关键基因BnATCAD4、BnF3H以及BnANS为turquoise模块的枢纽基因(hub gene)。通过已知的黄籽相关基因, 挖掘出了一部分黄籽微效作用基因, 这些基因多参与苯丙烷、类黄酮以及原花青素代谢途径。本研究挖掘的这些位点和候选基因可作为影响油菜黄籽形成的重要候选区域和基因, 有助于探究甘蓝型油菜黄籽基因资源信息、揭示油菜黄籽性状的遗传基础和分子机制、丰富分子育种理论以及提高油菜品质。     

利用WGCNA鉴定非生物胁迫相关基因共表达网络

加权共表达网络分析(Weighted Gene Co-expression Network Analysis, WGCNA)是用来描述不同样品之间基因关联模式的系统生物学方法,可以用来鉴定高度协同变化的基因集。本研究利用正常水稻组织共47份转录组数据,通过冷胁迫、干旱胁迫、盐胁迫不同的处理方式,使用WGCNA方法,根据已克隆基因的报道与以上3种胁迫相关的关键基因,探究不同逆境下基因之间的调控关系。通过对低表达量基因的过滤,最终利用筛选的30,339个表达的基因来构建共表达矩阵,得到15个模块。分析发现已知的水稻3种相关基因在各个模块均有存在,于是对预测到的靶基因进行GO富集分析。对3种胁迫下处理的转录组数据进行差异表达基因分析,结合已报道与胁迫相关的基因,选取各胁迫相关的2个模块进行了基因调控网络的构建。鉴于3种胁迫相关基因在green模块中大量分布,通过对green模块下各自特有的基因和共有的基因的GO功能富集分析,并对共有的基因构建调控网络,挖掘到2599个与3种胁迫都相关的基因,并预测出25个抗逆相关的关键基因,为水稻的抗逆及综合抗逆能力等研究提供了新思路。     

利用WGCNA进行玉米花期基因共表达模块鉴定

权重基因共表达网络分析(weighted gene co-expression network analysis, WGCNA)是系统生物学的一种研究方法,在挖掘生物学数据与特定性状之间的生物学关系方面具有十分重要的作用。本研究利用玉米(Zea mays L.)自交系B73的14份不同发育阶段的转录组数据,筛选掉低表达丰度的基因,最终得到了22,426个高表达的基因用于创建基因表达矩阵;利用不同组织作为性状,创建表型矩阵。然后利用R软件中的WGCNA包建立了共表达网络,共得到20个模块。本研究将与组织相关性高于0.65的模块定义为组织特异性模块,最终鉴定到14个组织特异性模块。利用在线网站Agrigo对组织特异性模块中的基因进行GO (gene ontology)富集分析,发现14个模块中均可以得到富集种类。开花作为玉米生育周期中的一个重要生理过程,不仅代表着植物从营养生长到生殖生长的转变,也关系到产量、株高和抗逆性等农艺性状。本研究发现8个组织特异性模块中的基因可以富集到与开花调控的代谢通路。此外,有17个已经报道过的开花时间调控基因存在于共表达模块中,并且主要分布在Blue模块和Darkmagenta模块,因此本研究重点关注了这2个模块内部的基因调控网络。本研究通过计算不同组织中的基因表达丰度,并联合权重基因共表达网络分析的方法,鉴定到了具有生物学意义的共表达基因模块,挖掘到了数个开花相关的模块,有助于揭示玉米开花调控的遗传机制。     

甘蓝型油菜响应低氮胁迫的表达谱分析

随着人们对作物产量的需求不断提高,氮肥被过量施用,而作物的氮素利用率(NUE)却在不断降低。本研究从低氮胁迫下油菜的生理变化入手,结合高通量的数字基因表达谱测序技术,分析了油菜在低氮0、3、72h下的转录组差异响应,鉴定了氮的吸收﹑转运﹑分配和转录因子等方面的差异表达基因。结果表明,甘蓝型油菜在低氮处理后,氮优先分配到地上部,硝酸还原酶(NR)活性显著降低,而谷氨酞胺合成酶(GS)活性升高,油菜植株总氮浓度降低, NUE升高。基因基因本位论(GO)功能与京都基因和基因组百科全书(KEGG)代谢通路分析表明,地上部差异表达基因主要是参与代谢过程、蛋白结合、离子结合、阴离子结合等,根中差异表达基因主要是参与分子功能、初级代谢过程、离子结合、阴离子结合等。基因表达谱分析表明,低氮胁迫72 h后,根中BnaGLNs家族基因表达大部分升高;根中BnaWRKY33s和BnaWRKY70s的基因表达量降低;BnaMYB4s、BnaMYB44s和BnaMYB51s亚家族中的大部分基因的表达量降低;BnaNIAs家族中的大部分基因表达上调;在BnaNRT2.1s和BnaNRT3.1s亚家族中,根中Bna...