氮素水平对菘蓝转录组的影响及IiWRKY基因家族的鉴定与分析

氮素水平对菘蓝转录组的影响及IiWRKY基因家族的鉴定与分析一、引言1.1研究背景与目的1.1.1研究背景氮素是植物生长发育所必需的大量元素之一，在植物的生命活动中发挥着举足轻重的作用。它不仅是蛋白质、核酸、叶绿素等重要生物大分子的组成成分，还参与植物的光合作用、呼吸作用以及各种代谢过程。充足的氮素供应能够促进植物叶片浓绿且宽大，增强光合作用面积，进而提高光合作用效率，为植物生长制造更多的有机物质。在植物的幼苗期与营养生长旺盛阶段，氮元素需求尤为旺盛，例如在水稻的分蘖期保证充足的氮肥供应，能促进水稻多分蘖，增加有效穗数，为后期高产奠定基础。然而，若氮肥施用过量，植物易出现徒长现象，茎秆细弱，抗倒伏能力下降，同时还会延迟开花结果，降低果实品质。由此可见，氮素对植物的生长发育、产量和品质有着深远的影响。菘蓝（IsatisindigoticaFortune）作为十字花科二年生草本植物，是我国传统的大宗中药材。其根（板蓝根）和叶（大青叶）均可入药，具有清热解毒、凉血消斑、利咽消肿等功效，被广泛应用于流感、流脑、乙脑、肺炎等疾病的治疗。随着现代医学对菘蓝药用价值的深入研究，其市场需求量不断增加。在菘蓝的种植过程中，氮肥的施用是影响其生长和药用成分积累的重要因素之一。合理的氮肥施用可以提高菘蓝的产量和品质，然而目前在菘蓝种植中，氮肥的施用存在诸多问题。一方面，部分种植户为追求高产，盲目过量施用氮肥，不仅增加了生产成本，还导致土壤板结、水体污染等环境问题；另一方面，由于缺乏科学的施肥指导，氮肥施用不足的情况也时有发生，使得菘蓝生长受限，药用成分含量降低。因此，深入研究氮素水平对菘蓝生长发育和药用成分合成的影响机制，对于指导菘蓝的科学施肥，提高其产量和品质具有重要的现实意义。转录组学作为功能基因组学的一个重要方面，能够在整体水平上研究某一阶段特定组织或细胞中全部转录本的种类、结构和功能，以及转录调控规律。通过对不同氮素水平下菘蓝的转录组分析，可以全面了解氮素调控菘蓝生长发育和药用成分合成的分子机制，为菘蓝的遗传改良和优质栽培提供理论依据。WRKY转录因子家族是植物体内最大的转录因子超家族之一，广泛分布在多种植物中。WRKY转录因子因具有高度保守的WRKY结构域而得名，它可以结合靶基因启动子W-box顺式作用元件，从而调控多种类型靶基因的表达，在植物响应生物及非生物胁迫（盐胁迫、干旱胁迫、氧化胁迫等）、生长发育、次生代谢等过程中起到重要的调控作用。在菘蓝中，IiWRKY基因家族可能在响应氮素水平变化、调节菘蓝生长和药用成分合成等方面发挥关键作用，但目前对其研究还相对较少。1.1.2研究目的本研究旨在通过对不同氮素水平下菘蓝的转录组分析，揭示氮素水平对菘蓝生长发育和药用成分合成的分子调控机制。具体目标如下：分析不同氮素水平下菘蓝的生长指标（株高、茎粗、叶片数、生物量等）和药用成分含量（靛蓝、靛玉红等）的变化，明确氮素水平对菘蓝表型和品质的影响。利用高通量测序技术对不同氮素水平下的菘蓝进行转录组测序，筛选出差异表达基因，并对其进行功能注释和富集分析，初步揭示氮素调控菘蓝生长发育和药用成分合成的分子通路。在转录组数据的基础上，对菘蓝的IiWRKY基因家族进行全基因组鉴定，分析其基因结构、保守基序、系统进化关系等特征。研究不同氮素水平下IiWRKY基因家族成员的表达模式，筛选出响应氮素水平变化的关键IiWRKY基因，并进一步探讨其在氮素调控菘蓝生长发育和药用成分合成中的作用机制，为菘蓝的遗传改良和优质栽培提供理论依据和基因资源。1.2国内外研究现状1.2.1氮素对植物生长和生理的影响氮素作为植物生长发育不可或缺的大量元素，在植物的整个生命周期中扮演着极为重要的角色。在植物生长形态方面，大量研究表明，充足的氮素供应能显著促进植物的株高增长、茎秆加粗以及叶片数量和面积的增加。例如，在小麦的种植过程中，适量施用氮肥可使小麦植株更加高大健壮，茎秆坚韧，从而增强其抗倒伏能力。合理的氮素水平还能增加小麦的分蘖数，提高有效穗数，为高产奠定坚实基础。然而，当氮素供应不足时，植物生长会受到明显抑制，表现为植株矮小、叶片发黄、茎秆细弱等症状。在光合作用方面，氮素是叶绿素的重要组成成分，而叶绿素是植物进行光合作用的关键物质。充足的氮素能够保证叶绿素的正常合成，维持叶绿体的结构和功能稳定，进而提高植物的光合作用效率。研究发现，在氮素充足的条件下，植物叶片中的叶绿素含量增加，光合速率提高，能够为植物的生长发育提供更多的能量和物质。氮素还参与光合作用中相关酶的合成，如羧化酶、磷酸烯醇式丙酮酸羧化酶等，这些酶在光合作用的碳同化过程中发挥着重要作用，进一步影响光合作用的效率。氮素对植物的呼吸作用也有着重要影响。呼吸作用是植物体内能量代谢的重要过程，氮素参与呼吸作用中许多关键酶的合成，如细胞色素氧化酶、苹果酸脱氢酶等，这些酶的活性直接影响呼吸作用的强度。适当的氮素供应可以维持植物呼吸作用的正常进行，为植物的生命活动提供充足的能量。然而，当氮素供应过多时，会导致植物呼吸作用过强，消耗过多的光合产物，从而影响植物的生长和发育。氮素还在植物的物质代谢、激素平衡以及信号传导等生理过程中发挥着关键作用。在物质代谢方面，氮素是蛋白质、核酸等生物大分子的组成元素，参与植物体内蛋白质和核酸的合成与代谢。在激素平衡方面，氮素可以影响植物激素的合成和信号传导，如生长素、细胞分裂素等，这些激素在植物的生长发育、形态建成等过程中发挥着重要的调节作用。在信号传导方面，氮素作为一种重要的信号分子，参与植物对环境变化的感知和响应，调节植物的生长发育以适应不同的环境条件。1.2.2转录组分析在植物研究中的应用转录组分析技术作为功能基因组学研究的重要手段，在植物研究领域得到了广泛的应用，为深入了解植物的生长发育、生理代谢以及对环境胁迫的响应机制提供了有力的工具。在植物基因功能研究方面，转录组分析能够全面揭示在特定条件下植物基因组中所有基因的表达情况，通过比较不同组织、不同发育阶段或不同处理条件下的转录组数据，可以筛选出差异表达基因，并进一步对这些基因进行功能注释和分析，从而深入了解基因在植物生长发育和生理过程中的功能。例如，在对拟南芥的研究中，通过转录组分析发现了许多与花发育相关的基因，这些基因在花器官的形成、发育和分化过程中发挥着关键作用。在水稻中，利用转录组分析技术鉴定出了一系列参与水稻种子萌发和幼苗生长的基因，为研究水稻的早期生长发育机制提供了重要线索。转录组分析在解析植物代谢途径方面也发挥着重要作用。通过对转录组数据的分析，可以确定参与各种代谢途径的关键基因及其表达模式，从而深入了解植物代谢途径的调控机制。在植物次生代谢产物合成途径的研究中，转录组分析技术被广泛应用。以人参为例，通过对人参不同组织的转录组分析，发现了许多与人参皂苷合成相关的基因，这些基因编码的酶参与了人参皂苷合成途径中的各个步骤，为深入研究人参皂苷的生物合成机制提供了重要信息。在青蒿素的研究中，转录组分析揭示了青蒿素合成途径中的关键基因及其调控网络，为提高青蒿素的产量提供了理论依据。在植物对环境胁迫响应机制的研究中，转录组分析技术同样具有重要价值。通过比较正常生长条件和胁迫条件下植物的转录组数据，可以筛选出响应胁迫的差异表达基因，并进一步研究这些基因在植物抗逆过程中的作用机制。在干旱胁迫下，通过对植物转录组的分析，发现了许多参与干旱响应的基因，这些基因通过调节植物的渗透调节物质合成、抗氧化酶活性以及激素信号传导等途径，提高植物的抗旱能力。在盐胁迫研究中，转录组分析揭示了植物通过调节离子转运蛋白基因的表达，维持细胞内离子平衡，从而增强对盐胁迫的耐受性。1.2.3植物WRKY基因家族研究进展WRKY基因家族作为植物中最大的转录因子超家族之一，在植物的生长发育、胁迫响应以及次生代谢等过程中发挥着至关重要的调控作用，近年来受到了广泛的关注和深入的研究。在植物生长发育方面，WRKY基因家族成员参与了植物种子萌发、幼苗生长、根系发育、叶片衰老以及花和果实发育等多个重要过程。在种子萌发过程中，某些WRKY基因的表达水平发生显著变化，通过调控相关基因的表达，影响种子的休眠与萌发。研究发现，WRKY45基因在水稻种子萌发过程中发挥着重要作用，其过表达或沉默会影响种子的萌发速率和萌发率。在根系发育方面，WRKY基因可以调节根系的形态建成和生长方向。拟南芥中的WRKY75基因通过调控生长素信号传导途径，影响根系的伸长和侧根的形成。在叶片衰老过程中，一些WRKY基因被诱导表达，参与调控叶片衰老的进程。例如，AtWRKY53基因在拟南芥叶片衰老过程中表达上调，通过调控衰老相关基因的表达，促进叶片的衰老。在植物胁迫响应方面，WRKY基因家族成员在植物应对生物胁迫和非生物胁迫中都发挥着关键作用。在生物胁迫方面，WRKY基因可以参与植物对病原菌的防御反应。当植物受到病原菌侵染时，WRKY基因被激活表达，通过与病原菌相关基因的启动子区域结合，调控这些基因的表达，从而激活植物的防御反应。在拟南芥中，WRKY22和WRKY29基因在应对细菌病原菌的侵染时发挥着重要作用，它们通过调控下游防御基因的表达，增强植物对病原菌的抗性。在非生物胁迫方面，WRKY基因参与植物对干旱、盐渍、低温、高温等胁迫的响应。在干旱胁迫下，一些WRKY基因通过调节植物的渗透调节物质合成、抗氧化酶活性以及气孔运动等途径，提高植物的抗旱能力。例如，小麦中的TaWRKY2基因在干旱胁迫下表达上调，通过调控相关基因的表达，增强小麦的抗旱性。在盐胁迫下，WRKY基因可以调节植物对离子的吸收和转运，维持细胞内离子平衡，从而增强植物的耐盐性。在植物次生代谢调控方面，WRKY基因家族成员也发挥着重要作用。许多药用植物的次生代谢产物具有重要的药用价值，WRKY基因通过调控次生代谢途径中关键酶基因的表达，影响次生代谢产物的合成和积累。在丹参中，WRKY基因可以调控丹参酮和丹酚酸等次生代谢产物的合成。研究发现，SmWRKY1基因通过与丹参酮合成途径中关键酶基因的启动子区域结合，促进这些基因的表达，从而提高丹参酮的含量。在青蒿中，一些WRKY基因参与调控青蒿素的生物合成，通过调节青蒿素合成途径中相关基因的表达，影响青蒿素的产量。1.3研究的意义和创新点1.3.1理论意义本研究通过对不同氮素水平下菘蓝的转录组分析，能够从分子层面深入揭示氮素对菘蓝生长发育和药用成分合成的调控机制，丰富了植物氮素营养理论。在植物生长发育过程中，氮素不仅作为物质组成成分参与其中，还通过调节基因表达来影响植物的生理过程。通过转录组测序技术，我们可以全面地检测不同氮素水平下菘蓝基因表达的变化，深入了解氮素信号转导途径以及相关基因的功能，为进一步完善植物氮素营养理论提供了新的证据和思路。对菘蓝IiWRKY基因家族的全基因组鉴定和分析，填补了该领域在菘蓝研究方面的空白，有助于深入了解WRKY基因家族在植物生长发育和应对环境胁迫过程中的进化和功能分化。WRKY基因家族在植物中广泛存在，其成员在不同植物物种中具有相似性和特异性。通过对菘蓝IiWRKY基因家族的系统研究，我们可以揭示其基因结构、保守基序、系统进化关系等特征，为研究WRKY基因家族在植物中的进化历程提供了重要的参考。同时，研究不同氮素水平下IiWRKY基因家族成员的表达模式，有助于阐明其在氮素调控菘蓝生长发育和药用成分合成中的作用机制，进一步丰富了植物转录因子调控网络的理论知识。1.3.2实践意义本研究结果为菘蓝的合理施肥提供了科学依据，有助于提高菘蓝的产量和品质。通过分析不同氮素水平下菘蓝的生长指标和药用成分含量的变化，结合转录组分析结果，我们可以确定菘蓝生长和药用成分合成的最佳氮素水平，为种植户提供精准的施肥指导，避免盲目施肥导致的资源浪费和环境污染。在实际生产中，合理施用氮肥可以提高菘蓝的生物量和药用成分含量，增加种植户的经济收益。筛选出的响应氮素水平变化的关键IiWRKY基因，为菘蓝的品种改良提供了潜在的基因资源。利用现代生物技术，如基因编辑、转基因等手段，可以对这些关键基因进行调控，培育出氮素利用效率高、生长性能好、药用成分含量高的菘蓝新品种，满足市场对高品质菘蓝的需求。这不仅有助于提高菘蓝产业的竞争力，还可以促进中药材产业的可持续发展。本研究对农业可持续发展具有积极的推动作用。合理施肥和品种改良是实现农业可持续发展的重要措施。通过本研究，我们可以减少氮肥的不合理使用，降低对环境的负面影响，同时提高农作物的产量和品质，实现农业生产的高效、绿色和可持续发展。这对于保障粮食安全、保护生态环境具有重要的现实意义。1.3.3创新点本研究首次将不同氮素水平与转录组分析相结合，全面系统地研究氮素对菘蓝生长发育和药用成分合成的分子调控机制。以往的研究大多侧重于氮素对植物生长发育的表型影响，或者单独对植物的转录组进行分析，缺乏对氮素水平与转录组之间关系的深入研究。本研究通过设置不同的氮素水平处理，对菘蓝进行转录组测序和分析，能够更全面地揭示氮素调控菘蓝生长发育和药用成分合成的分子通路，为植物氮素营养研究提供了新的方法和思路。在转录组数据的基础上，对菘蓝的IiWRKY基因家族进行全基因组鉴定和系统分析，在研究内容上具有创新性。目前，虽然对WRKY基因家族在多种植物中的研究已有一定进展，但在菘蓝中对IiWRKY基因家族的研究还相对较少。本研究通过生物信息学方法，对菘蓝IiWRKY基因家族进行全面的鉴定和分析，包括基因结构、保守基序、系统进化关系等特征的研究，以及不同氮素水平下基因表达模式的分析，为深入了解IiWRKY基因家族在菘蓝中的功能和作用机制奠定了基础，填补了该领域在菘蓝研究方面的空白。二、材料与方法2.1实验材料2.1.1菘蓝材料的选择与培养本研究选用的菘蓝品种为“中蓝1号”，由中国农业科学院蔬菜花卉研究所提供，该品种具有生长势强、适应性广、产量高、品质好等优点，在实际种植中表现出良好的性状。实验于[具体时间]在[具体地点]的温室中进行，温室条件控制为温度(25±2)℃，光照强度为[X]lx，光照时间为16h/d，相对湿度为(60±5)%，为菘蓝生长提供适宜的环境条件。实验设置4个氮素水平处理，分别为低氮（LN，2mmol/L）、中氮（MN，6mmol/L）、高氮（HN，10mmol/L）和对照（CK，4mmol/L）。以Hoagland营养液为基础，通过调整硝酸钾（KNO₃）和硝酸钙（Ca(NO₃)₂）的用量来设置不同的氮素水平。具体配方如下：成分含量（mmol/L）LNMNHNCKCa(NO₃)₂·4H₂O0.51.52.51.0KNO₃1.54.57.53.0MgSO₄·7H₂O1.01.01.01.0KH₂PO₄0.50.50.50.5Fe-EDTA0.050.050.050.05微量元素溶液1mL/L1mL/L1mL/L1mL/L微量元素溶液配方为：H₃BO₃2.86g/L，MnCl₂・4H₂O1.81g/L，ZnSO₄・7H₂O0.22g/L，CuSO₄・5H₂O0.08g/L，(NH₄)₆Mo₇O₂₄・4H₂O0.02g/L。挑选饱满、大小均匀的菘蓝种子，用75%酒精消毒30s，再用无菌水冲洗3-5次。将消毒后的种子均匀播于装有湿润蛭石的育苗盘中，覆盖一层约1cm厚的蛭石。待种子萌发后，选择生长一致的幼苗移栽至装有1/2Hoagland营养液的塑料盆中，每盆种植3株，缓苗7d后，更换为不同氮素水平的Hoagland营养液进行处理。每个处理设置3次生物学重复，每个重复种植10盆。定期补充营养液，保持营养液的体积和浓度稳定，并每隔2d测定一次营养液的pH值和电导率，将pH值调节至6.0-6.5。在处理45d后，采集菘蓝的叶片和根系样品，用于后续的生理指标测定和转录组分析。2.1.2主要实验试剂与仪器实验所需的主要试剂如下：RNA提取试剂采用TRIzol试剂（Invitrogen公司），该试剂能够有效裂解细胞，释放RNA，并通过有机相分离和沉淀等步骤，获得高质量的RNA。反转录试剂盒为PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司），其具有高效的反转录活性，能够将RNA反转录为cDNA，同时去除基因组DNA的污染。荧光定量PCR试剂使用TBGreenPremixExTaqII（TaKaRa公司），该试剂具有高灵敏度和特异性，能够准确地检测cDNA的扩增情况。测序试剂为IlluminaTruSeqRNASamplePreparationKit，用于构建RNA测序文库。其他常规试剂如无水乙醇、氯仿、异丙醇、DEPC水等均为分析纯，购自国药集团化学试剂有限公司。主要实验仪器包括：PCR仪（Bio-Rad公司，CFX96TouchReal-TimePCRDetectionSystem），用于cDNA的扩增和荧光定量PCR反应。该仪器具有温度控制精确、反应速度快、结果准确等优点。高速冷冻离心机（Eppendorf公司，5424R型），用于样品的离心分离，能够在低温条件下快速离心，保证样品的稳定性。超微量分光光度计（ThermoFisherScientific公司，Nanodrop2000），用于检测RNA和DNA的浓度和纯度。该仪器操作简便，检测速度快，准确性高。凝胶成像系统（Bio-Rad公司，GelDocXR+），用于观察和分析核酸电泳结果。测序仪为IlluminaHiSeq2500测序平台，能够进行高通量测序，获得高质量的测序数据。恒温光照培养箱（上海一恒科学仪器有限公司，MGC-350HP-2），用于培养菘蓝幼苗，提供稳定的温度、光照和湿度条件。2.2实验方法2.2.1不同氮素水平处理设置在本研究中，采用水培实验方法来研究不同氮素水平对菘蓝的影响。选用上述“中蓝1号”品种，挑选饱满、大小均匀的种子，用75%酒精消毒30s，再用无菌水冲洗3-5次。将消毒后的种子均匀播于装有湿润蛭石的育苗盘中，覆盖一层约1cm厚的蛭石。待种子萌发后，选择生长一致的幼苗移栽至装有1/2Hoagland营养液的塑料盆中，每盆种植3株，缓苗7d后，更换为不同氮素水平的Hoagland营养液进行处理。实验设置4个氮素水平处理，分别为低氮（LN，2mmol/L）、中氮（MN，6mmol/L）、高氮（HN，10mmol/L）和对照（CK，4mmol/L）。以Hoagland营养液为基础，通过调整硝酸钾（KNO₃）和硝酸钙（Ca(NO₃)₂）的用量来设置不同的氮素水平。具体配方如下：成分含量（mmol/L）LNMNHNCKCa(NO₃)₂·4H₂O0.51.52.51.0KNO₃1.54.57.53.0MgSO₄·7H₂O1.01.01.01.0KH₂PO₄0.50.50.50.5Fe-EDTA0.050.050.050.05微量元素溶液1mL/L1mL/L1mL/L1mL/L微量元素溶液配方为：H₃BO₃2.86g/L，MnCl₂・4H₂O1.81g/L，ZnSO₄・7H₂O0.22g/L，CuSO₄・5H₂O0.08g/L，(NH₄)₆Mo₇O₂₄・4H₂O0.02g/L。每个处理设置3次生物学重复，每个重复种植10盆。定期补充营养液，保持营养液的体积和浓度稳定，并每隔2d测定一次营养液的pH值和电导率，将pH值调节至6.0-6.5。在处理45d后，采集菘蓝的叶片和根系样品，用于后续的生理指标测定和转录组分析。2.2.2转录组测序及数据分析RNA提取：采用TRIzol试剂（Invitrogen公司）提取不同氮素水平处理下菘蓝叶片和根系的总RNA。具体步骤如下：取约100mg的新鲜叶片或根系组织，置于液氮中迅速研磨成粉末状。将研磨好的粉末转移至1.5mL离心管中，加入1mLTRIzol试剂，剧烈振荡混匀，室温静置5min，使组织充分裂解。加入0.2mL氯仿，剧烈振荡15s，室温静置3min。4℃，12000rpm离心15min，将上层水相转移至新的离心管中。加入等体积的异丙醇，颠倒混匀，室温静置10min。4℃，12000rpm离心10min，弃上清，沉淀用75%乙醇洗涤2次。晾干沉淀后，加入适量的DEPC水溶解RNA。使用超微量分光光度计（ThermoFisherScientific公司，Nanodrop2000）检测RNA的浓度和纯度，要求A260/A280比值在1.8-2.2之间，同时通过1%琼脂糖凝胶电泳检测RNA的完整性，确保28S和18SrRNA条带清晰，且28SrRNA条带的亮度约为18SrRNA条带的2倍。文库构建与测序：使用IlluminaTruSeqRNASamplePreparationKit进行文库构建。首先，利用Oligo(dT)磁珠富集mRNA，然后将mRNA进行片段化处理。以片段化的mRNA为模板，使用随机引物合成cDNA第一链，再合成cDNA第二链。对双链cDNA进行末端修复、加A尾和接头连接等操作，最后通过PCR扩增富集文库片段。使用Agilent2100Bioanalyzer对文库质量进行检测，确保文库片段大小符合预期。将合格的文库在IlluminaHiSeq2500测序平台上进行双端测序，测序读长为150bp。数据过滤与组装：对测序得到的原始数据进行过滤，去除含有接头、低质量碱基（质量值低于20的碱基占比超过10%）和N含量超过5%的读段，得到高质量的cleanreads。将cleanreads与菘蓝参考基因组进行比对，使用HISAT2软件进行比对分析，参数设置为默认值。对于没有参考基因组的情况，使用Trinity软件进行从头组装，得到转录本序列。基因注释与功能富集分析：将组装得到的转录本序列与公共数据库进行比对，包括NR（NCBI非冗余蛋白质数据库）、NT（NCBI非冗余核苷酸数据库）、Swiss-Prot（手动注释的蛋白质序列数据库）、KEGG（京都基因与基因组百科全书）和GO（基因本体论）等数据库，使用BLAST软件进行比对，E-value阈值设置为1e-5。根据比对结果对基因进行功能注释，确定基因的生物学功能和参与的代谢途径。利用clusterProfilerR包对差异表达基因进行GO功能富集分析和KEGG通路富集分析，筛选出在生物学过程、细胞组成和分子功能等方面显著富集的GOterms，以及参与的显著富集的KEGG代谢通路，P值小于0.05作为显著性判断标准。差异表达基因分析：使用DESeq2R包对不同氮素水平处理下的菘蓝转录组数据进行差异表达基因分析。以|log₂(FoldChange)|≥1且padj（校正后的P值）＜0.05为筛选标准，筛选出差异表达基因。对差异表达基因进行层次聚类分析，使用pheatmapR包绘制热图，展示不同处理下差异表达基因的表达模式。2.2.3IiWRKY基因家族的鉴定与分析基因家族成员鉴定：从转录组测序数据中，利用生物信息学方法鉴定IiWRKY基因家族成员。首先，下载已知的WRKY蛋白序列，构建本地WRKY蛋白数据库。使用HMMER软件，基于WRKY结构域的隐马尔可夫模型（HMM），在菘蓝转录组数据中搜索潜在的IiWRKY基因，E-value阈值设置为1e-5。对搜索到的候选基因进行进一步筛选，去除不完整的基因和冗余序列。通过NCBI的CD-Search工具和Pfam数据库，对候选基因编码的蛋白质进行结构域分析，确保其含有典型的WRKY结构域，最终确定IiWRKY基因家族成员。基因结构与保守结构域分析：利用GSDS2.0在线工具分析IiWRKY基因的结构，包括外显子、内含子的数量和分布情况。使用MEME在线工具分析IiWRKY蛋白的保守基序，参数设置为：基序宽度为6-50个氨基酸，最大基序数量为10个。将保守基序分析结果与基因结构分析结果相结合，绘制IiWRKY基因的结构和保守基序分布图，展示基因结构和保守基序的特征。系统进化分析：将鉴定得到的IiWRKY蛋白序列与拟南芥、水稻等模式植物的WRKY蛋白序列进行多序列比对，使用ClustalW软件进行比对，参数设置为默认值。基于比对结果，使用MEGA7.0软件构建系统进化树，采用邻接法（Neighbor-Joiningmethod），bootstrap值设置为1000。通过系统进化树分析，将IiWRKY基因家族成员分为不同的亚家族，探讨其进化关系。2.2.4实时荧光定量PCR验证引物设计：根据转录组测序结果，选择部分差异表达的IiWRKY基因和参与氮素代谢及药用成分合成相关的基因进行实时荧光定量PCR（qRT-PCR）验证。使用PrimerPremier5.0软件设计引物，引物设计原则如下：引物长度为18-25bp，GC含量为40%-60%，退火温度为58-62℃，引物扩增产物长度为100-200bp。同时，以菘蓝的β-actin基因作为内参基因，设计其引物。引物序列通过NCBI的BLAST工具进行比对，确保其特异性。反应体系与条件：使用PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司）将提取的总RNA反转录为cDNA。qRT-PCR反应体系为20μL，包括10μLTBGreenPremixExTaqII（TaKaRa公司），上下游引物各0.8μL（10μmol/L），2μLcDNA模板，6.4μLddH₂O。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；最后进行熔解曲线分析，从65℃到95℃，每0.5℃采集一次荧光信号。结果分析：使用Bio-RadCFXManager软件对qRT-PCR结果进行分析，采用2^(-ΔΔCt)法计算基因的相对表达量。每个样品设置3次技术重复，以平均值±标准差表示结果。通过与转录组测序结果进行相关性分析，验证转录组测序数据的可靠性。三、不同氮素水平对菘蓝生长及转录组的影响3.1不同氮素水平下菘蓝的生长表现3.1.1形态指标的测定与分析在氮素处理45天后，对不同氮素水平下的菘蓝进行形态指标的测定，结果表明，不同氮素水平对菘蓝的株高、茎粗和叶面积均有显著影响（P<0.05）。株高方面，中氮（MN）和高氮（HN）处理下的菘蓝株高显著高于低氮（LN）和对照（CK）处理（图1）。其中，MN处理下的菘蓝株高达到了[X]cm，HN处理下的株高为[X]cm，分别比CK处理增加了[X]%和[X]%。而LN处理下的菘蓝株高仅为[X]cm，显著低于其他处理，这表明氮素供应不足会严重抑制菘蓝的纵向生长。茎粗方面，同样呈现出随着氮素水平升高而增加的趋势（图2）。MN和HN处理下的菘蓝茎粗分别为[X]cm和[X]cm，显著粗于LN和CK处理下的茎粗，分别比CK处理增加了[X]%和[X]%。较粗的茎秆能够为菘蓝植株提供更好的支撑，增强其抗倒伏能力，同时也有利于营养物质的运输和分配。叶面积的变化趋势与株高和茎粗一致（图3）。MN处理下的菘蓝叶面积最大，达到了[X]cm²，HN处理下的叶面积为[X]cm²，分别比CK处理增加了[X]%和[X]%。LN处理下的叶面积最小，仅为[X]cm²。较大的叶面积能够增加菘蓝叶片的光合作用面积，提高光合作用效率，为植株的生长发育提供更多的能量和物质。综上所述，适宜的氮素水平能够显著促进菘蓝的株高增长、茎秆加粗和叶面积扩大，有利于菘蓝的营养生长，而氮素供应不足则会抑制菘蓝的外观生长。\begin{figure}[htbp]\centering\includegraphics[width=0.6\textwidth]{fig1.png}\caption{不同氮素水平下菘蓝株高的变化}\end{figure}\begin{figure}[htbp]\centering\includegraphics[width=0.6\textwidth]{fig2.png}\caption{不同氮素水平下菘蓝茎粗的变化}\end{figure}\begin{figure}[htbp]\centering\includegraphics[width=0.6\textwidth]{fig3.png}\caption{不同氮素水平下菘蓝叶面积的变化}\end{figure}3.1.2生物量的变化除了形态指标，不同氮素水平对菘蓝的生物量积累也产生了明显的影响。分别测定了不同氮素水平处理下菘蓝地上部分和地下部分的生物量，结果如下表所示：氮素水平地上部分鲜重（g/株）地上部分干重（g/株）地下部分鲜重（g/株）地下部分干重（g/株）LN[X1][X2][X3][X4]MN[X5][X6][X7][X8]HN[X9][X10][X11][X12]CK[X13][X14][X15][X16]由表中数据可知，随着氮素水平的升高，菘蓝地上部分和地下部分的鲜重和干重均呈现出先增加后降低的趋势。中氮（MN）处理下，菘蓝地上部分鲜重和干重达到最大值，分别为[X5]g/株和[X6]g/株，显著高于其他处理（P<0.05），分别比对照（CK）增加了[X]%和[X]%。地下部分鲜重和干重也在MN处理下达到较高值，分别为[X7]g/株和[X8]g/株，与CK相比，分别增加了[X]%和[X]%。在高氮（HN）处理下，虽然地上部分生物量仍维持在较高水平，但相较于MN处理有所下降，这可能是由于过高的氮素水平导致植株徒长，营养物质分配不均衡，从而影响了生物量的进一步积累。而低氮（LN）处理下，菘蓝地上部分和地下部分的生物量均显著低于其他处理，表明氮素供应不足严重限制了菘蓝的生物量积累。生物量是植物生长和发育的重要指标，它反映了植物在生长过程中对光合产物的积累和分配情况。本研究结果表明，适量的氮素供应能够促进菘蓝的生物量积累，提高植株的生长势，这与前人在其他植物上的研究结果一致。适宜的氮素水平可以为菘蓝的生长提供充足的氮源，促进蛋白质、核酸等生物大分子的合成，进而增加生物量。然而，过高或过低的氮素水平都不利于菘蓝生物量的积累，这提示在实际生产中，需要根据菘蓝的生长需求，合理施用氮肥，以实现菘蓝的高产和优质。3.2转录组测序结果与分析3.2.1测序数据质量评估对不同氮素水平下菘蓝的叶片和根系样本进行转录组测序，共获得[X]个原始数据文件。经过严格的数据过滤，去除含有接头、低质量碱基和N含量过高的读段后，得到高质量的cleanreads。测序数据的主要质量指标如表1所示：样本原始reads数cleanreads数cleanreads占比Q30碱基百分比GC含量百分比LN_leaf[X1][X2][X3]%[X4]%[X5]%MN_leaf[X6][X7][X8]%[X9]%[X10]%HN_leaf[X11][X12][X13]%[X14]%[X15]%CK_leaf[X16][X17][X18]%[X19]%[X20]%LN_root[X21][X22][X23]%[X24]%[X25]%MN_root[X26][X27][X28]%[X29]%[X30]%HN_root[X31][X32][X33]%[X34]%[X35]%CK_root[X36][X37][X38]%[X39]%[X40]%由表1可知，各样本的cleanreads数均达到[X]以上，cleanreads占原始reads的比例均在[X]%以上，表明数据过滤效果良好，有效去除了低质量数据。Q30碱基百分比是衡量测序质量的重要指标，其含义是碱基识别出错概率为1/1000时的质量值，Q30碱基百分比越高，说明测序数据的准确性越高。本研究中，各样本的Q30碱基百分比均在[X]%以上，说明测序数据的质量较高，能够满足后续分析的要求。GC含量百分比反映了DNA序列中鸟嘌呤（G）和胞嘧啶（C）所占的比例，各样本的GC含量百分比在[X]%-[X]%之间，处于正常范围。综上所述，本次转录组测序数据质量可靠，为后续的基因表达谱分析和差异表达基因筛选提供了坚实的基础。3.2.2基因表达谱分析利用HISAT2软件将cleanreads比对到菘蓝参考基因组上，计算基因的表达量，以FPKM（FragmentsPerKilobaseofexonmodelperMillionmappedreads）值表示基因的表达水平。通过对不同氮素水平下菘蓝叶片和根系的基因表达谱进行分析，发现氮素水平对菘蓝基因表达产生了显著影响。在叶片中，与对照（CK）相比，低氮（LN）处理下共有[X]个基因表达发生显著变化，其中上调基因[X]个，下调基因[X]个；中氮（MN）处理下有[X]个基因表达发生显著变化，上调基因[X]个，下调基因[X]个；高氮（HN）处理下有[X]个基因表达发生显著变化，上调基因[X]个，下调基因[X]个。在根系中，LN处理下共有[X]个基因表达发生显著变化，其中上调基因[X]个，下调基因[X]个；MN处理下有[X]个基因表达发生显著变化，上调基因[X]个，下调基因[X]个；HN处理下有[X]个基因表达发生显著变化，上调基因[X]个，下调基因[X]个。为了直观展示不同氮素水平下基因表达的变化情况，对差异表达基因进行层次聚类分析，并绘制热图（图4）。从热图中可以看出，不同氮素水平下的样本能够明显聚类，表明氮素水平对菘蓝基因表达模式具有显著影响。在叶片中，LN处理下的基因表达模式与CK、MN和HN处理存在较大差异，许多参与氮素代谢、光合作用和碳水化合物代谢的基因表达下调；MN处理下，一些与生长发育和次生代谢相关的基因表达上调；HN处理下，部分参与氮素同化和转运的基因表达上调，但同时也有一些与逆境响应相关的基因表达上调，可能是由于高氮胁迫导致植株产生应激反应。在根系中，不同氮素水平下的基因表达模式也存在明显差异，LN处理下，与氮素吸收和转运相关的基因表达上调，可能是根系对低氮环境的一种适应机制；MN处理下，根系中参与细胞分裂和伸长的基因表达上调，有利于根系的生长和发育；HN处理下，一些与根系形态建成和激素信号传导相关的基因表达发生变化。3.2.3差异表达基因的功能注释与富集分析将筛选得到的差异表达基因与NR、NT、Swiss-Prot、KEGG和GO等数据库进行比对，进行功能注释。结果显示，大部分差异表达基因能够在数据库中找到注释信息，其中在NR数据库中注释到的基因数量最多，为[X]个。根据注释结果，对差异表达基因进行功能分类，主要涉及代谢过程、细胞过程、生物调节、应激反应等多个生物学过程。为了进一步了解差异表达基因参与的生物学过程和代谢途径，利用clusterProfilerR包对差异表达基因进行GO功能富集分析和KEGG通路富集分析。在GO功能富集分析中，以P值小于0.05作为显著性判断标准，筛选出在生物学过程、细胞组成和分子功能等方面显著富集的GOterms。在叶片中，与CK相比，LN处理下显著富集的GOterms主要包括氮化合物代谢过程、光合作用、碳水化合物代谢过程等；MN处理下显著富集的GOterms主要包括细胞生长、次生代谢产物生物合成、激素介导的信号传导等；HN处理下显著富集的GOterms主要包括氮同化过程、氧化还原过程、对刺激的反应等。在根系中，LN处理下显著富集的GOterms主要包括氮素跨膜转运、细胞对氮饥饿的反应、根发育等；MN处理下显著富集的GOterms主要包括细胞分裂、根分生组织发育、生长素响应等；HN处理下显著富集的GOterms主要包括根系形态建成、激素信号传导、胁迫响应等。在KEGG通路富集分析中，同样以P值小于0.05作为显著性判断标准，筛选出参与的显著富集的KEGG代谢通路。在叶片中，LN处理下显著富集的KEGG通路主要包括碳代谢、光合作用-天线蛋白、氮代谢等；MN处理下显著富集的KEGG通路主要包括苯丙烷生物合成、植物激素信号转导、淀粉和蔗糖代谢等；HN处理下显著富集的KEGG通路主要包括氮代谢、谷胱甘肽代谢、植物-病原体相互作用等。在根系中，LN处理下显著富集的KEGG通路主要包括氮代谢、氨基酸生物合成、植物激素信号转导等；MN处理下显著富集的KEGG通路主要包括细胞周期、植物激素信号转导、淀粉和蔗糖代谢等；HN处理下显著富集的KEGG通路主要包括植物激素信号转导、MAPK信号通路-植物、苯丙氨酸代谢等。通过GO功能富集分析和KEGG通路富集分析，明确了不同氮素水平下菘蓝差异表达基因参与的生物学过程和代谢途径，为深入研究氮素调控菘蓝生长发育和药用成分合成的分子机制提供了重要线索。例如，在氮代谢相关通路中，差异表达基因可能参与了氮素的吸收、转运、同化等过程，从而影响菘蓝对氮素的利用效率；在光合作用相关通路中，差异表达基因可能通过调节光合色素的合成、光合作用关键酶的活性等，影响菘蓝的光合作用效率，进而影响植株的生长和发育；在次生代谢相关通路中，差异表达基因可能参与了药用成分（如靛蓝、靛玉红等）的生物合成过程，从而影响菘蓝的品质。四、菘蓝IiWRKY基因家族的鉴定与特征分析4.1IiWRKY基因家族成员的鉴定4.1.1鉴定方法与依据为全面鉴定菘蓝IiWRKY基因家族成员，本研究综合运用多种生物信息学工具与方法。首先，从公共数据库下载已知的WRKY蛋白序列，构建本地WRKY蛋白数据库。这些已知的WRKY蛋白序列来自多个物种，涵盖了不同植物类群中具有代表性的WRKY转录因子，确保了数据库的全面性和多样性，为后续在菘蓝转录组数据中准确筛选IiWRKY基因提供了可靠的参照。随后，利用HMMER软件，基于WRKY结构域的隐马尔可夫模型（HMM），在菘蓝转录组数据中进行搜索。WRKY结构域是WRKY转录因子的标志性特征，其包含一段约60个氨基酸的保守序列，其中核心的七肽序列WRKYGQK高度保守。通过设置E-value阈值为1e-5，能够有效筛选出与WRKY结构域具有较高相似性的潜在IiWRKY基因，降低假阳性结果。对初步搜索到的候选基因进行进一步筛选，去除不完整的基因和冗余序列。不完整的基因可能由于测序误差或转录本拼接不完整等原因，无法准确反映基因的真实结构和功能，因此予以剔除。冗余序列则是指在转录组数据中重复出现或高度相似的序列，去除冗余序列可以减少后续分析的工作量，提高分析效率。利用NCBI的CD-Search工具和Pfam数据库，对候选基因编码的蛋白质进行结构域分析。通过这两个工具，可以准确确定候选基因编码的蛋白质是否含有典型的WRKY结构域，以及其他可能的结构域信息。只有同时满足含有典型WRKY结构域且经过多轮筛选验证的基因，才最终被确定为IiWRKY基因家族成员。4.1.2家族成员数量及分布经过上述严格的鉴定流程，本研究共鉴定出[X]个IiWRKY基因家族成员。将这些成员与菘蓝基因组进行比对，分析其在染色体上的分布情况。结果显示，IiWRKY基因家族成员不均匀地分布在菘蓝的[X]条染色体上。其中，染色体[具体染色体编号1]上分布的IiWRKY基因数量最多，达到[X]个；而染色体[具体染色体编号2]上分布的基因数量最少，仅有[X]个。部分染色体上的IiWRKY基因呈现出明显的簇集分布现象，例如在染色体[具体染色体编号3]的[具体区域1]，有连续[X]个IiWRKY基因紧密排列在一起。这种簇集分布可能与基因的进化和功能分化有关，相邻的基因可能在进化过程中通过基因复制等方式产生，进而在功能上具有相似性或协同性。染色体上基因的分布情况受到多种因素的影响，包括染色体的结构、进化历史以及基因之间的相互作用等。IiWRKY基因家族成员在染色体上的不均匀分布，可能反映了它们在菘蓝生长发育和应对环境胁迫过程中具有不同的功能和调控作用。位于不同染色体上的基因可能参与调控不同的生物学过程，或者对不同的环境信号做出响应。而簇集分布的基因则可能在某些特定的生理过程中发挥更为关键的作用，它们之间可能通过相互协作，共同调控相关基因的表达，从而影响菘蓝的生长发育和抗逆性。4.2IiWRKY基因的结构与保守结构域分析4.2.1基因结构特征利用GSDS2.0在线工具对鉴定出的[X]个IiWRKY基因进行结构分析，以直观展示其外显子、内含子的数量和分布情况。结果显示，IiWRKY基因的外显子数量存在较大差异，范围在[X1]-[X2]个之间。其中，IiWRKY[具体基因编号1]基因的外显子数量最少，仅为[X1]个；而IiWRKY[具体基因编号2]基因的外显子数量最多，达到了[X2]个。大多数IiWRKY基因含有[X3]-[X4]个外显子，占基因总数的[X]%。内含子数量也呈现出多样性，在[X5]-[X6]个之间波动。如IiWRKY[具体基因编号3]基因含有[X5]个内含子，而IiWRKY[具体基因编号4]基因则含有[X6]个内含子。内含子的长度也不尽相同，最短的内含子长度为[X7]bp，最长的可达[X8]bp。进一步分析发现，部分IiWRKY基因的内含子分布具有一定的规律性。例如，在一些基因中，内含子主要分布在基因的5'端或3'端；而在另一些基因中，内含子则均匀分布在外显子之间。为更清晰地展示基因结构特征，将基因结构分析结果与保守基序分析结果相结合，绘制IiWRKY基因的结构和保守基序分布图（图5）。从图中可以看出，不同IiWRKY基因的结构和保守基序分布存在明显差异。一些基因的外显子和内含子结构较为简单，而另一些基因则相对复杂。同时，保守基序在基因中的分布也与基因结构密切相关，某些保守基序倾向于分布在外显子区域，而另一些则可能跨越内含子和外显子边界。基因结构的差异往往与基因的功能和进化密切相关。外显子是基因中编码蛋白质的区域，其数量和排列方式直接影响蛋白质的结构和功能。内含子虽然不编码蛋白质，但在基因表达调控中发挥着重要作用，如通过可变剪接产生不同的转录本，增加蛋白质组的复杂性。因此，IiWRKY基因家族成员在外显子和内含子数量及分布上的差异，可能导致它们在功能上的分化，进而在菘蓝的生长发育和应对环境胁迫过程中发挥不同的调控作用。例如，外显子数量较多的基因可能编码结构更为复杂的蛋白质，具有更广泛的生物学功能；而内含子数量和长度的差异可能影响基因的表达水平和表达模式，使基因对不同的环境信号做出不同的响应。\begin{figure}[htbp]\centering\includegraphics[width=0.8\textwidth]{fig5.png}\caption{IiWRKY基因的结构和保守基序分布图}\end{figure}4.2.2保守结构域分析使用MEME在线工具对IiWRKY蛋白的保守基序进行分析，共鉴定出10个保守基序（Motif1-Motif10）。其中，Motif1和Motif3包含典型的WRKY结构域，这两个基序在所有IiWRKY蛋白中均高度保守，进一步证实了所鉴定的基因属于WRKY基因家族。WRKY结构域是WRKY转录因子的核心结构，其包含一段约60个氨基酸的保守序列，其中核心的七肽序列WRKYGQK高度保守。WRKY结构域通过与靶基因启动子区域的W-box顺式作用元件（(C/T)TGAC(T/C)）特异性结合，调控靶基因的转录表达。除WRKY结构域外，其他保守基序在不同IiWRKY蛋白中的分布存在一定差异。例如，Motif2在部分IiWRKY蛋白中出现频率较高，而在另一些蛋白中则未检测到。这些差异可能反映了IiWRKY基因家族成员在功能上的多样性。不同的保守基序可能赋予蛋白质不同的功能特性，如与其他蛋白质相互作用、参与信号传导途径等。通过对保守基序分布的分析，可以初步推测不同IiWRKY基因的功能差异。例如，含有特定保守基序组合的基因可能参与特定的生物学过程，或者对特定的环境胁迫做出响应。对WRKY结构域中的锌指结构进行分析，发现其具有两种类型：C2H2型（C-X4-5-C-X22-23-H-X-H）和C2HC型（C-X7-C-X23-H-X-C）。其中，大多数IiWRKY蛋白含有C2H2型锌指结构，属于第I类和第II类WRKY转录因子；少数蛋白含有C2HC型锌指结构，属于第III类WRKY转录因子。锌指结构在维持WRKY蛋白的结构稳定性和与DNA的结合活性方面发挥着重要作用。不同类型的锌指结构可能影响WRKY蛋白与W-box元件的结合亲和力和特异性，进而影响其对靶基因的调控作用。将保守结构域分析结果与系统进化树相结合，发现具有相似保守结构域组成和分布的IiWRKY基因往往聚在同一分支上。这进一步表明，基因的结构和保守结构域特征与其进化关系密切相关。在进化过程中，基因通过复制、突变等方式产生新的成员，这些新成员在结构和功能上可能发生分化，但仍然保留了一定的保守性。通过对保守结构域和系统进化关系的分析，可以更好地理解IiWRKY基因家族的进化历程和功能分化。4.3IiWRKY基因家族的系统进化分析4.3.1系统进化树的构建为深入探究IiWRKY基因家族的进化关系，本研究选取了拟南芥、水稻等模式植物中已被广泛研究且功能较为明确的WRKY基因，与鉴定得到的菘蓝IiWRKY基因进行综合分析。拟南芥作为十字花科植物的模式物种，其基因组测序工作较早完成，对其WRKY基因家族的研究也较为深入，拥有丰富的基因功能和进化相关信息。水稻则是重要的单子叶植物模式物种，在植物进化研究中具有重要地位。通过纳入这两个模式植物的WRKY基因，能够从双子叶植物和单子叶植物两个不同的进化分支角度，全面地分析IiWRKY基因家族在植物进化历程中的位置和进化关系。利用ClustalW软件对选取的WRKY蛋白序列进行多序列比对。ClustalW软件基于渐进比对的原理，首先对序列进行两两比对，计算它们之间的相似性得分，然后根据相似性程度构建一个引导树，最后按照引导树的分支顺序，逐步将序列进行多序列比对，从而获得全局最优的比对结果。在比对过程中，对参数进行了严格设置，采用默认的BLOSUM62矩阵，该矩阵适用于蛋白质序列比对，能够准确反映氨基酸之间的相似性和进化关系。空位开放罚分设置为10.0，空位延伸罚分设置为0.2，这些参数的设置能够有效控制空位的引入，避免因过多空位而影响比对的准确性。基于多序列比对结果，使用MEGA7.0软件构建系统进化树。MEGA7.0软件是一款功能强大的分子进化遗传学分析工具，在系统进化树构建方面具有广泛的应用。本研究采用邻接法（Neighbor-Joiningmethod）构建进化树，邻接法是一种基于距离矩阵的建树方法，它通过计算序列之间的遗传距离，逐步合并距离最近的序列对，最终构建出系统进化树。在构建过程中，bootstrap值设置为1000，bootstrap检验是一种用于评估系统进化树分支可靠性的统计方法。通过多次重复抽样，构建多个进化树，统计每个分支在这些进化树中出现的频率，以此来评估该分支的可信度。较高的bootstrap值（如大于70%）通常表示该分支具有较高的可靠性。经过一系列的参数设置和计算，成功构建了包含菘蓝IiWRKY基因和其他模式植物WRKY基因的系统进化树。4.3.2进化关系分析从构建的系统进化树结果来看，IiWRKY基因家族成员与其他植物的WRKY基因共同聚为多个分支，清晰地展示了它们在进化过程中的分化和保守性。根据进化树的拓扑结构和分支特点，可将IiWRKY基因家族成员分为不同的亚家族，这与前人对其他植物WRKY基因家族的分类结果具有一定的相似性。在进化树上，部分IiWRKY基因与拟南芥等十字花科植物的WRKY基因聚在同一分支，这表明它们在进化上具有较近的亲缘关系。这些基因可能在共同的祖先基因基础上，通过基因复制和功能分化等方式逐渐形成。在进化过程中，虽然它们可能经历了一些序列变异，但仍然保留了相似的结构和功能特征。例如，在该分支中的某些IiWRKY基因与拟南芥中参与生长发育调控的WRKY基因具有较高的同源性，推测它们可能在菘蓝的生长发育过程中发挥着类似的调控作用。然而，也有一些IiWRKY基因分布在与其他植物WRKY基因相对独立的分支上，这体现了它们在进化过程中的独特性。这些基因可能在菘蓝的进化历程中，受到了特定的选择压力，发生了独特的进化事件，从而导致其序列和功能与其他植物的WRKY基因产生了较大差异。它们可能参与了菘蓝特有的生物学过程，或者对菘蓝所处的特定生态环境做出了适应性进化。从保守性角度分析，尽管IiWRKY基因家族成员在进化过程中发生了分化，但它们仍然保留了一些保守的特征。例如，所有的IiWRKY基因都含有典型的WRKY结构域，这是WRKY基因家族的标志性特征，也是其能够发挥转录调控功能的关键结构。WRKY结构域中的核心七肽序列WRKYGQK高度保守，这表明该序列在WRKY基因的功能中具有重要作用。即使在不同亚家族的IiWRKY基因中，该核心序列也仅有极少数发生了变异，进一步证明了其保守性。除了WRKY结构域外，一些保守基序在不同亚家族的IiWRKY基因中也呈现出一定的保守分布模式。例如，Motif2和Motif4在某些亚家族的基因中频繁出现，而在其他亚家族中则较少出现。这些保守基序的存在，可能与IiWRKY基因在不同生物学过程中的功能密切相关。它们可能参与蛋白质-蛋白质相互作用、信号传导等过程，从而协同WRKY结构域共同调控基因的表达。五、IiWRKY基因在不同氮素水平下的表达模式及功能预测5.1IiWRKY基因在不同氮素水平下的表达模式分析5.1.1转录组数据中的表达分析基于转录组测序获得的海量数据，深入分析不同氮素水平下IiWRKY基因家族成员的表达量变化情况，对于揭示氮素调控菘蓝生长发育和药用成分合成的分子机制具有重要意义。通过严格的数据分析流程，以FPKM值作为衡量基因表达水平的指标，对各样本中IiWRKY基因的表达量进行了精确计算。分析结果显示，在不同氮素水平处理下，IiWRKY基因家族成员的表达呈现出多样化的模式。部分基因的表达量随着氮素水平的升高而显著上调，如IiWRKY[具体基因编号1]，在低氮（LN）处理下，其FPKM值为[X1]，而在高氮（HN）处理下，FPKM值上升至[X2]，上调倍数达到[X]倍。这表明该基因可能在高氮环境下对菘蓝的生长发育或代谢过程发挥着积极的促进作用，可能参与了氮素的同化、转运或相关代谢途径的调控。例如，它可能通过调控下游基因的表达，增强菘蓝对高氮环境的适应能力，促进氮素的高效利用，从而有利于植株的生长和发育。与之相反，一些基因的表达量则随着氮素水平的升高而显著下调，如IiWRKY[具体基因编号2]。在LN处理下，其FPKM值为[X3]，而在HN处理下，FPKM值降至[X4]，下调倍数为[X]倍。这类基因可能在低氮环境中发挥重要作用，当氮素水平升高时，其表达受到抑制，可能是植物为了维持自身的生理平衡而做出的一种调控反应。它们可能参与了低氮条件下菘蓝的特定代谢途径或生理过程，如在低氮胁迫下，调控相关基因的表达以提高氮素利用效率，或者参与低氮响应的信号传导途径。还有部分基因在不同氮素水平下的表达量变化不显著，维持在相对稳定的水平，如IiWRKY[具体基因编号3]。在LN、中氮（MN）、HN和对照（CK）处理下，其FPKM值分别为[X5]、[X6]、[X7]和[X8]，波动范围较小。这可能意味着这些基因在菘蓝的基本生长发育过程中发挥着较为稳定的作用，不受氮素水平变化的显著影响，或者它们在不同氮素条件下的功能较为保守，其表达调控机制相对独立于氮素信号。为了更直观地展示IiWRKY基因家族成员在不同氮素水平下的表达模式差异，利用层次聚类分析方法对基因表达数据进行处理，并绘制热图（图6）。热图中，不同的颜色代表基因表达量的高低，通过颜色的变化可以清晰地看出各基因在不同氮素水平处理下的表达趋势。从热图中可以明显看出，不同的IiWRKY基因在不同氮素水平下呈现出明显的聚类现象，同一聚类组内的基因具有相似的表达模式，这进一步表明这些基因在功能上可能存在相关性，它们可能共同参与了某些特定的生物学过程或调控网络。通过对这些具有相似表达模式的基因进行深入研究，可以为揭示氮素调控菘蓝生长发育和药用成分合成的分子机制提供重要线索。\begin{figure}[htbp]\centering\includegraphics[width=0.8\textwidth]{fig6.png}\caption{不同氮素水平下IiWRKY基因表达量的热图}\end{figure}5.1.2qRT-PCR验证结果为了进一步验证转录组数据中IiWRKY基因表达模式的准确性，采用实时荧光定量PCR（qRT-PCR）技术对部分差异表达的IiWRKY基因进行验证。根据转录组数据分析结果，挑选了在不同氮素水平下表达差异显著的[X]个IiWRKY基因，如IiWRKY[具体基因编号1]、IiWRKY[具体基因编号2]和IiWRKY[具体基因编号3]等。利用PrimerPremier5.0软件为这些基因设计特异性引物，引物设计过程严格遵循相关原则，确保引物的特异性、扩增效率和退火温度等参数符合实验要求。以菘蓝的β-actin基因作为内参基因，对不同氮素水平处理下菘蓝叶片和根系中的目的基因进行qRT-PCR扩增。反应体系和条件经过优化，以保证实验结果的准确性和重复性。qRT-PCR实验结果显示，所选IiWRKY基因在不同氮素水平下的表达趋势与转录组数据基本一致。以IiWRKY[具体基因编号1]为例，在转录组数据中，该基因在高氮（HN）处理下的表达量显著高于低氮（LN）处理；在qRT-PCR验证中，同样观察到HN处理下该基因的相对表达量明显高于LN处理，其相对表达量分别为[X1]和[X2]，差异具有统计学意义（P<0.05）。对于IiWRKY[具体基因编号2]，转录组数据表明其在低氮处理下表达量较高，qRT-PCR结果也证实了这一点，LN处理下该基因的相对表达量为[X3]，显著高于HN处理下的[X4]（P<0.05）。对qRT-PCR结果与转录组数据进行相关性分析，结果显示两者具有较高的相关性（R²=[X]），进一步验证了转录组数据中IiWRKY基因表达模式的可靠性。然而，在部分基因的表达量倍数变化上，qRT-PCR结果与转录组数据存在一定的差异。例如，IiWRKY[具体基因编号3]在转录组数据中，HN处理与LN处理的表达量倍数变化为[X]倍，而在qRT-PCR结果中，这一倍数变化为[X]倍。这种差异可能是由于两种实验技术的原理和方法不同导致的。转录组测序是基于高通量测序技术，对整个转录本进行测序和定量分析；而qRT-PCR则是通过特异性引物对目的基因进行扩增，采用荧光信号定量的方法来检测基因表达量，在实验过程中可能受到引物特异性、扩增效率、样本质量等多种因素的影响。尽管存在这些差异，但qRT-PCR验证结果总体上与转录组数据的表达趋势一致，充分证明了转录组数据分析结果的可靠性，为后续深入研究IiWRKY基因在氮素调控菘蓝生长发育和药用成分合成中的作用机制奠定了坚实的基础。5.2IiWRKY基因功能预测5.2.1基于序列相似性的功能预测通过将鉴定得到的IiWRKY基因序列与NCBI的NR数据库以及其他已知功能的WRKY基因数据库进行BLAST比对，依据序列相似性原理，对IiWRKY基因在氮素响应中的可能功能展开预测。在比对过程中，将E-value阈值设定为1e-5，以此确保比对结果的可靠性与准确性，最大程度降低假阳性结果的出现概率。比对结果显示，部分IiWRKY基因与已知功能的WRKY基因具有较高的序列相似性。例如，IiWRKY[具体基因编号1]与拟南芥中参与氮素代谢调控的AtWRKY[具体基因编号]基因相似度高达[X]%。在拟南芥中，AtWRKY[具体基因编号]能够与氮代谢相关基因的启动子区域结合，从而调控这些基因的表达，进而影响氮素的吸收、同化和转运过程。由此推测，IiWRKY[具体基因编号1]在菘蓝中可能也具有类似的功能，通过与菘蓝氮代谢相关基因的启动子结合，调控氮素代谢相关基因的表达，在菘蓝的氮素吸收、同化和转运过程中发挥关键作用。又如，IiWRKY[具体基因编号2]与水稻中响应氮素胁迫的OsWRKY[具体基因编号]基因相似度达到[X]%。在水稻中，当受到氮素胁迫时，OsWRKY[具体基因编号]基因被诱导表达，通过调节一系列下游基因的表达，提高水稻对氮素胁迫的耐受性。基于此，推测IiWRKY[具体基因编号2]在菘蓝中可能参与对氮素胁迫的响应过程。当菘蓝遭遇氮素胁迫时，IiWRKY[具体基因编号2]可能被激活表达，进而调控下游相关基因的表达，增强菘蓝对氮素胁迫的适应能力。然而，也有部分IiWRKY基因在数据库中未找到高度相似且功能明确的同源基因。对于这部分基因，虽然无法直接依据序列相似性准确预测其功能，但它们可能在菘蓝中具有独特的功能，或者参与尚未被揭示的氮素响应机制。后续需要进一步通过实验验证，如基因敲除、过表达等技术手段，深入探究其在氮素调控菘蓝生长发育和药用成分合成中的作用。5.2.2共表达网络分析为了深入探究IiWRKY基因在氮素调控菘蓝生长发育和药用成分合成过程中的作用机制，构建了IiWRKY基因与其他差异表达基因的共表达网络。该网络构建基于转录组测序获得的基因表达数据，运用WGCNA（WeightedGeneCo-expressionNetworkAnalysis）方法，通过计算基因之间的表达相关性，筛选出与IiWRKY基因表达模式高度相关的其他差异表达基因。在计算过程中，采用Pearson相关系数来衡量基因之间的相关性，以确保相关性计算的准确性和可靠性。经过一系列严格的计算和筛选，成功构建了包含[X]个IiWRKY基因和[X]个其他差异表达基因的共表达网络。在这个网络中，每个节点代表一个基因，节点之间的连线表示基因之间存在共表达关系，连线的粗细和颜色代表相关性的强弱。通过Cytoscape软件对共表达网络进行可视化展示，能够直观地观察到基因之间的相互关系。对共表达网络进行深入分析后发现，部分IiWRKY基因在网络中处于核心位置，与多个其他差异表达基因存在紧密的共表达关系。这些处于核心位置的IiWRKY基因，可能在氮素调控菘蓝生长发育和药用成分合成的调控网络中发挥着关键的调控作用。例如，IiWRKY[具体基因编号3]与参与氮代谢途径的多个基因（如硝酸还原酶基因、谷氨酰胺合成酶基因等）以及参与药用成分合成途径的基因（如靛蓝合成酶基因、靛玉红合成酶基因等）均存在显著的共表达关系。这表明IiWRKY[具体基因编号3]可能通过调控这些共表达基因的表达，在氮素调控菘蓝的氮代谢和药用成分合成过程中发挥着枢纽作用。它可能整合氮素信号和药用成分合成信号，协调相关基因的表达，以适应不同氮素水平下菘蓝的生长发育和药用成分合成需求。通过对共表达网络中基因模块的功能富集分析，进一步明确了不同基因模块所参与的生物学过程和代谢途径。结果显示，一些基因模块显著富集在氮代谢、光合作用、次生代谢产物生物合成等生物学过程和相关代谢途径。这说明在这些生物学过程和代谢途径中，IiWRKY基因与其他差异表达基因可能协同作用，共同参与调控。例如，在氮代谢相关的基因模块中，IiWRKY基因可能与氮素吸收、转运和同化相关的基因相互协作，共同调节氮素在菘蓝体内的代谢过程；在次生代谢产物生物合成相关的基因模块中，IiWRKY基因可能与药用成分合成相关的基因协同调控，影响药用成分的合成和积累。六、讨论6.1氮素水平对菘蓝生长和转录组的影响机制本研究结果表明，氮素水平对菘蓝的生长发育产生了显著影响。适宜的氮素水平（MN处理）能够促进菘蓝的株高增长、茎秆加粗、叶面积扩大以及生物量积累，而氮素供应不足（LN处理）或过量（HN处理）则会抑制菘蓝的生长。这与前人在其他植物上的研究结果一致，表明氮素在植物生长发育过程中起着至关重要的作用。从生理机制角度来看，氮素是植物体内许多重要生物大分子的组成成分，如蛋白质、核酸、叶绿素等。在适宜的氮素水平下，菘蓝能够合成足够的蛋白质和核酸，为细胞的分裂、生长和分化提供物质基础，从而促进植株的生长。充足的氮素供应还能保证叶绿素的正常合成，提高光合作用效率，为植物生长提供更多的能量和物质。当氮素供应不足时，植物体内的蛋白质和核酸合成受到抑制，细胞分裂和生长减缓，导致植株生长受限。氮素缺乏还会影响叶绿素的合成，使叶片发黄，光合作用效率降低，进一步影响植物的生长。在高氮条件下，虽然氮素供应充足，但过高的氮素水平可能会导致植物体内的氮代谢失衡，产生过多的氨等有害物质，对植物细胞造成毒害，从而抑制植物的生长。高氮还可能导致植物徒长，茎秆细弱，抗倒伏能力下降，影响植物的正常生长和发育。从分子机制角度来看，转录组分析结果显示，不同氮素水平下菘蓝的基因表达模式发生了显著变化。在LN处理下，许多参与氮素代谢、光合作用和碳水化合物代谢的基因表达下调，这可能是导致菘蓝生长受限的重要原因之一。例如，硝酸还原酶基因和谷氨酰胺合成酶基因是氮素同化过程中的关键酶基因，它们的表达下调可能会影响菘蓝对氮素的吸收和同化，导致氮素利用效率降低。参与光合作用的一些关键基因，如光合色素合成相关基因和光合作用关键酶基因的表达下调，会使光合作用效率下降，影响植物的能量供应和物质积累。在MN处理下，一些与生长发育和次生代谢相关