解析白菜BrMYB80b正义lncRNA在育性调控中的分子机制与功能

解析白菜BrMYB80b正义lncRNA在育性调控中的分子机制与功能一、引言1.1研究背景植物育性调控是植物生长发育过程中的关键环节，对农作物的产量和品质起着决定性作用。在农业生产中，作物育性的稳定和优化直接关系到粮食安全和农业可持续发展。例如，杂交育种技术的广泛应用，极大地提高了作物的产量和抗逆性，而这一技术的核心就在于对植物育性的精准调控。随着分子生物学技术的飞速发展，长链非编码RNA（lncRNA）作为一类重要的调控分子，逐渐成为植物发育研究领域的热点。lncRNA是长度大于200个核苷酸的非编码RNA，虽不具备编码蛋白质的能力，但在基因表达调控、细胞分化、发育等诸多生物学过程中发挥着关键作用。越来越多的研究表明，lncRNA在植物育性调控中扮演着不可或缺的角色。在拟南芥中，某些lncRNA的表达变化会显著影响花粉的发育和育性，进而影响植株的繁殖能力。在水稻中，特定的lncRNA参与了雄性不育的调控过程，为水稻杂种优势的利用提供了重要的理论基础。白菜（BrassicarapaL.）作为一种重要的蔬菜作物，在全球范围内广泛种植，其产量和品质对农业经济和人们的生活有着重要影响。然而，目前对于白菜育性调控的分子机制研究仍相对较少，尤其是lncRNA在其中的作用尚未得到充分揭示。BrMYB80b正义lncRNA作为白菜中发现的一种潜在的育性调控相关lncRNA，对其进行深入研究，有望揭示白菜育性调控的新机制，为白菜的遗传改良和新品种选育提供理论支持和技术手段。1.2研究目的与意义本研究旨在深入探究白菜BrMYB80b正义lncRNA在育性调控中的作用机制，通过对该lncRNA的表达模式、功能验证以及与其他相关基因的互作关系进行系统研究，揭示其在白菜育性调控网络中的关键节点作用，为植物育性调控理论的完善提供新的分子证据。从基因表达层面来看，明确BrMYB80b正义lncRNA在不同育性状态下白菜植株中的表达差异，以及其在花粉发育、雌蕊形成等关键育性相关过程中的表达动态，有助于深入理解植物育性调控的分子基础，丰富植物发育生物学中关于育性调控的理论体系。在白菜育种实践中，本研究具有重要的应用价值。当前白菜育种面临着诸多挑战，如品种同质化严重、对环境胁迫的适应性不足等。通过揭示BrMYB80b正义lncRNA在育性调控中的作用，有望为白菜的遗传改良提供新的靶点和策略。利用基因编辑技术对该lncRNA进行精准调控，或者筛选与该lncRNA相关的分子标记，用于辅助选择具有优良育性性状的白菜品种，从而提高白菜的育种效率，加快新品种的选育进程，为白菜产业的可持续发展提供坚实的技术支撑，满足日益增长的市场需求，保障蔬菜供应的稳定和质量。二、文献综述2.1植物育性调控机制概述植物育性调控是一个高度复杂且精细的生物学过程，涉及众多基因、信号通路以及环境因素的协同作用。深入了解植物育性调控机制，对于揭示植物生殖奥秘、推动农业育种技术发展具有重要意义。下面将从花粉发育过程及关键调控因素、绒毡层在植物育性中的功能以及转录因子在植物育性调控中的作用这三个方面进行阐述。2.1.1花粉发育过程及关键调控因素花粉发育是植物育性调控的核心环节之一，其过程从花粉母细胞开始，历经多个复杂且有序的阶段，最终形成成熟花粉。花粉发育起始于花药中的孢原细胞，孢原细胞经过平周分裂形成造孢细胞，造孢细胞进一步有丝分裂产生花粉母细胞，即小孢子母细胞。花粉母细胞随后进入减数分裂阶段，经过两次连续的分裂，形成四个单倍体的小孢子，此时的小孢子呈四分体状态。在拟南芥中，减数分裂过程受到一系列基因的严格调控，如SPO11、REC8等基因，它们对于减数分裂过程中染色体的正确配对、交换和分离至关重要，任何一个基因的突变都可能导致减数分裂异常，进而影响花粉的正常发育。四分体时期结束后，小孢子从四分体中分离出来，进入单核花粉粒阶段。此时的单核花粉粒经历核分裂，形成一个营养核和一个生殖核，进入二细胞花粉阶段。在部分植物中，生殖核还会进一步分裂，形成两个精子，从而进入三细胞花粉阶段。在这个过程中，营养核主要负责为花粉的萌发和花粉管的生长提供营养物质和能量，而生殖核及其分裂产生的精子则承担着遗传物质传递的重要使命。在玉米中，研究发现一些基因如ZmPME1、ZmPME2等在花粉发育后期发挥着重要作用，它们参与调控花粉壁的形成和花粉的萌发，这些基因的表达异常会导致花粉壁结构缺陷，影响花粉的正常萌发和受精能力。花粉发育过程受到多种关键因素的严格调控，基因和激素在其中扮演着至关重要的角色。众多基因构成了一个复杂而精密的调控网络，对花粉发育的各个阶段进行细致的调控。例如，在水稻中，MADS-box基因家族中的一些成员，如OsMADS1、OsMADS58等，参与调控花粉母细胞的减数分裂和小孢子的发育过程。这些基因通过与其他基因相互作用，调节相关基因的表达，从而确保花粉发育过程的顺利进行。激素作为植物体内重要的信号分子，在花粉发育过程中也发挥着不可或缺的作用。生长素能够促进花粉管的生长和伸长，使其能够顺利到达胚珠完成受精过程；赤霉素参与花粉发育中的减数分裂和花粉成熟过程，对花粉的正常发育具有重要影响；细胞分裂素则在花粉母细胞的分化和花药发育中发挥关键作用，调节细胞的分裂和分化，为花粉发育提供必要的细胞基础。在番茄中，研究发现生长素信号通路中的关键基因SlIAA9、SlARF7等参与调控花粉的发育和花粉管的生长，通过调节这些基因的表达，可以影响花粉的育性和植物的生殖能力。2.1.2绒毡层在植物育性中的功能绒毡层是花药壁的最内层细胞，在植物育性调控中具有不可或缺的重要功能，是花粉发育过程中至关重要的营养和物质供应源。绒毡层细胞具有独特的结构和生理特性，其细胞质浓厚，富含各种细胞器和营养物质，如蛋白质、脂类、多糖等。这些丰富的物质储备为花粉的发育提供了坚实的物质基础。在花粉发育的早期阶段，绒毡层细胞通过分泌各种酶和营养物质，为小孢子的发育提供必要的营养支持，促进小孢子的分裂和分化。在减数分裂过程中，绒毡层细胞分泌的胼胝质酶能够降解小孢子之间的胼胝质壁，使小孢子得以分离，进入单核花粉粒阶段。绒毡层在花粉壁的形成过程中发挥着关键作用。花粉壁是花粉的重要结构组成部分，对于保护花粉、维持花粉的正常生理功能以及促进花粉的传播和识别具有重要意义。绒毡层细胞合成并分泌孢粉素、外壁蛋白等物质，这些物质参与花粉外壁的形成，决定了花粉外壁的结构和特性。孢粉素是一种高度聚合的生物大分子，具有极强的化学稳定性和抗降解能力，它构成了花粉外壁的主要成分，赋予花粉外壁坚硬的结构和良好的保护性能。外壁蛋白则在花粉与雌蕊的识别和相互作用过程中发挥着重要作用，有助于花粉的黏附、萌发和花粉管的生长。在拟南芥中，研究发现绒毡层细胞中的一些基因如MS2、MS5等参与孢粉素的合成和运输过程，这些基因的突变会导致孢粉素合成异常，花粉外壁结构缺陷，从而使花粉丧失正常的功能，无法完成受精过程。绒毡层的发育异常会导致花粉败育，进而影响植物的育性。绒毡层发育受到严格的基因调控，一系列基因构成了复杂的调控网络，对绒毡层的分化、发育和功能发挥进行精细的调控。当这些基因发生突变或表达异常时，会导致绒毡层发育异常，出现细胞形态改变、物质合成和分泌障碍等问题，最终导致花粉败育。在油菜中，发现一些基因如BnMS1、BnMS2等对绒毡层的发育和功能具有重要调控作用，这些基因的突变会导致绒毡层细胞提前解体或功能异常，无法为花粉发育提供充足的营养和物质支持，从而使花粉发育受阻，出现败育现象。绒毡层发育异常还可能导致花粉壁形成缺陷，使花粉无法抵御外界环境的影响，降低花粉的活力和受精能力。2.1.3转录因子在植物育性调控中的作用转录因子是一类能够与基因启动子区域的顺式作用元件特异性结合，从而调控基因转录水平的蛋白质。在植物育性调控过程中，转录因子发挥着关键的调控作用，它们通过调节花粉发育、绒毡层功能发挥等相关基因的表达，构建起复杂而精细的调控网络，确保植物育性的正常维持。MYB转录因子家族是植物中最大的转录因子家族之一，在植物育性调控中具有重要作用。以拟南芥为例，AtMYB103基因在绒毡层和花粉发育过程中特异性表达，它通过直接调控一系列与绒毡层发育、花粉壁形成相关基因的表达，对花粉的正常发育起着关键的调控作用。研究表明，AtMYB103基因的突变会导致绒毡层细胞分化异常，花粉壁结构缺陷，最终使花粉败育。在水稻中，OsMYB80基因参与调控绒毡层的程序性死亡和花粉外壁的形成，该基因的表达异常会导致绒毡层细胞不能按时降解，花粉外壁发育不完全，影响花粉的育性。bHLH转录因子家族也在植物育性调控中发挥着不可或缺的作用。在拟南芥中，bHLH转录因子AMS对绒毡层的发育和功能至关重要。AMS基因的表达受到严格的时空调控，它在绒毡层发育的特定阶段表达，通过调控下游一系列基因的表达，参与绒毡层细胞的分化、发育以及物质合成和分泌等过程。当AMS基因发生突变时，绒毡层细胞的分化和功能会出现严重异常，无法为花粉发育提供必要的支持，导致花粉败育。在番茄中，SlbHLH1基因参与调控花粉的发育和花粉管的生长，通过调节相关基因的表达，影响花粉的活力和受精能力。AP2/ERF转录因子家族同样在植物育性调控中扮演着重要角色。在棉花中，GhAP2-1基因参与调控花粉的发育和花粉管的生长过程。研究发现，GhAP2-1基因通过调控下游与花粉发育相关基因的表达，影响花粉的正常发育和功能。当GhAP2-1基因的表达受到抑制时，花粉的萌发率降低，花粉管的生长受到阻碍，从而影响棉花的育性。在大豆中，GmERF3基因参与调控花器官的发育和花粉的育性，通过调节相关基因的表达，对大豆的生殖过程发挥重要的调控作用。2.2lncRNA的特性与功能2.2.1lncRNA的定义与分类长链非编码RNA（lncRNA）是一类长度大于200个核苷酸的非编码RNA分子，它们不具备编码蛋白质的能力，但在基因表达调控、细胞分化、发育等诸多生物学过程中发挥着重要的调控作用。lncRNA最初被认为是基因组转录的“噪音”，不具有生物学功能。随着研究的不断深入，越来越多的证据表明lncRNA在生物体内具有广泛而重要的生物学功能，是基因表达调控网络中的关键组成部分。在人类基因组中，大量的lncRNA参与了细胞增殖、分化、凋亡等重要生物学过程的调控，其表达异常与多种疾病的发生发展密切相关，如癌症、心血管疾病等。根据lncRNA在基因组上的位置和与邻近基因的关系，可将其分为正义lncRNA、反义lncRNA、内含子lncRNA、基因间lncRNA等类型。正义lncRNA与编码基因的转录方向相同，通常由编码基因的外显子或内含子转录产生，可通过与mRNA相互作用，影响mRNA的稳定性、翻译效率等。反义lncRNA与编码基因的转录方向相反，能够与mRNA形成双链结构，从而调控mRNA的表达，在转录水平或转录后水平发挥作用。内含子lncRNA则完全位于编码基因的内含子区域，其功能机制尚不完全清楚，可能参与基因转录的调控以及RNA的加工过程。基因间lncRNA位于两个编码基因之间，不与任何已知的编码基因重叠，这类lncRNA在植物中较为常见，它们可以通过与DNA、RNA或蛋白质相互作用，参与染色质修饰、转录激活或抑制等过程，进而调控基因表达。在拟南芥中，已鉴定出许多基因间lncRNA，它们在植物的生长发育、逆境响应等过程中发挥着重要作用。按照功能，lncRNA可分为信号lncRNA、诱饵lncRNA、引导lncRNA和支架lncRNA等。信号lncRNA能够响应外界环境信号或细胞内的生理状态变化，通过改变自身的表达水平来传递信号，进而调控下游基因的表达。诱饵lncRNA可以与转录因子、miRNA等调控分子结合，从而竞争性地抑制它们与靶基因的相互作用，起到“分子海绵”的作用。引导lncRNA能够与特定的蛋白质结合，引导其到基因组上的特定位置，参与基因表达的调控。支架lncRNA则可以作为支架，将多个蛋白质分子聚集在一起，形成功能性的复合物，参与生物过程的调控。在水稻中，一些lncRNA被发现具有诱饵功能，它们可以吸附miRNA，解除miRNA对靶基因的抑制作用，从而调控水稻的生长发育和对逆境的响应。2.2.2lncRNA在植物中的作用机制lncRNA在植物中主要通过表观遗传调控、转录调控和ceRNA（竞争性内源RNA）机制等多种方式发挥作用，从而精细地调控植物基因表达和生理过程，确保植物的正常生长发育。在表观遗传调控方面，lncRNA可以通过与DNA、组蛋白或相关的表观修饰酶相互作用，影响染色质的状态和结构，进而调控基因的表达。一些lncRNA能够招募DNA甲基转移酶，使特定基因的启动子区域发生DNA甲基化修饰，从而抑制基因的转录。在拟南芥中，研究发现某些lncRNA可以通过与组蛋白修饰复合物相互作用，调控组蛋白的甲基化、乙酰化等修饰状态，改变染色质的开放性，影响基因的表达。这种表观遗传调控方式具有稳定性和可遗传性，能够在植物的生长发育过程中持续发挥作用，对植物的表型和性状产生深远影响。在转录调控过程中，lncRNA可以作为转录因子的共激活因子或共抑制因子，与转录因子相互作用，影响转录起始复合物的形成和组装，从而调控基因的转录起始。某些lncRNA可以与RNA聚合酶结合，增强或抑制其与基因启动子的结合能力，直接影响基因的转录效率。lncRNA还可以通过形成RNA-DNA杂交链（R-loop）结构，影响基因转录的延伸过程。在水稻中，研究发现一些lncRNA能够与转录因子OsMADS1相互作用，调控水稻花器官发育相关基因的转录，对水稻的育性和产量产生重要影响。ceRNA机制是lncRNA在植物中发挥作用的另一种重要方式。lncRNA可以作为miRNA的靶标模拟物，与mRNA竞争结合相同的miRNA，从而解除miRNA对mRNA的抑制作用，调控mRNA的表达水平。这种ceRNA调控网络在植物中广泛存在，参与了植物生长发育、逆境响应等多个生物学过程的调控。在番茄中，研究发现lncRNA1642可以作为miR164的靶标模拟物，通过ceRNA机制调控NAC1基因的表达，进而影响番茄对盐胁迫的响应。2.2.3lncRNA在植物育性调控中的研究进展近年来，随着研究的不断深入，越来越多的证据表明lncRNA在植物育性调控中发挥着重要作用，它们参与了花粉发育、育性转换等关键过程，对植物的繁殖能力和后代的遗传多样性具有重要影响。在花粉发育方面，许多lncRNA被发现参与了花粉发育的调控过程。在拟南芥中，lncRNAMSTRG.16996在花粉发育过程中特异性表达，其功能缺失会导致花粉发育异常，花粉活力下降，从而影响植物的育性。进一步研究发现，MSTRG.16996可以通过与多个花粉发育相关基因相互作用，调控这些基因的表达，进而影响花粉的正常发育。在水稻中，OsLnc1是一个在花粉发育过程中起重要作用的lncRNA，它可以通过调控花粉壁的形成和花粉管的生长，影响水稻花粉的育性。研究表明，OsLnc1的表达受到多种因素的调控，其表达异常会导致花粉壁结构缺陷，花粉管生长受阻，最终导致花粉败育。lncRNA在植物育性转换中也发挥着关键作用。在一些植物中，育性转换是一个受环境因素和遗传因素共同调控的过程，lncRNA在其中扮演着重要的信号传导和调控角色。在光周期敏感型雄性不育水稻中，研究发现一些lncRNA的表达水平在不同光周期条件下发生显著变化，这些lncRNA可能通过调控相关基因的表达，参与光周期诱导的育性转换过程。在短日照条件下，某些lncRNA的表达上调，它们通过与相关转录因子相互作用，激活或抑制育性相关基因的表达，从而导致水稻表现为可育；而在长日照条件下，这些lncRNA的表达下调，育性相关基因的表达模式发生改变，水稻表现为雄性不育。在玉米中，研究发现lncRNAZm401参与了玉米的育性调控过程。Zm401可以通过与多个转录因子和miRNA相互作用，形成复杂的调控网络，影响玉米花粉的发育和育性。在Zm401功能缺失的突变体中，玉米花粉发育异常，花粉活力显著降低，导致玉米的结实率下降。进一步研究表明，Zm401可以通过调控花粉壁的形成、花粉管的生长以及花粉发育相关基因的表达，维持玉米花粉的正常发育和育性。三、材料与方法3.1实验材料本研究选用的白菜品种为“德高黄金改三”，该品种属于中晚熟类型，具备出色的抗逆能力，菜质鲜嫩，口感清甜，深受市场欢迎，其生育期约为80-85天，适合在8月中旬左右播种，在北方地区广泛种植，对当地的气候和土壤条件具有良好的适应性，为研究白菜在自然环境下的育性调控机制提供了理想的实验材料。实验白菜种子经消毒处理后，播种于装有营养土的育苗盘中。营养土由草炭土、蛭石和珍珠岩按照3:1:1的体积比混合而成，这种配比能够为白菜幼苗的生长提供充足的养分和良好的透气性。将育苗盘放置于光照培养箱中进行培养，培养条件设定为光照16h/黑暗8h，光照强度为300μmol・m-2・s-1，温度为22℃/18℃（昼/夜），相对湿度保持在60%-70%。在这样的环境条件下，能够模拟白菜在自然生长过程中的光照、温度和湿度需求，确保幼苗生长健壮。待白菜幼苗长至3-4片真叶时，将其移栽至塑料花盆中，花盆直径为20cm，高15cm，每盆种植1株。移栽后的白菜植株放置于温室中进行常规管理，定期浇水、施肥，施肥遵循“薄肥勤施”的原则，以氮肥为主，配合适量的磷钾肥，为白菜植株的生长提供充足的养分，保证其正常生长发育。在白菜植株的不同发育时期，分别采集花芽、叶片、茎尖、根等组织样本。花芽样本的采集在植株现蕾后，选取大小一致、发育正常的花芽；叶片样本采集自植株中部健康、无病虫害的叶片；茎尖样本则选取植株顶端生长旺盛的部分；根样本在小心挖掘植株后，取其主根和侧根的幼嫩部分。每个时期的每个组织样本设置3次生物学重复，每次重复采集3-5个植株的相应组织，以确保实验结果的准确性和可靠性。采集后的样本立即放入液氮中速冻，随后转移至-80℃冰箱中保存，用于后续的RNA提取和相关实验分析。3.2主要实验试剂与仪器实验用到的主要试剂包括RNA提取试剂TRIzol，购自Invitrogen公司，用于从白菜组织样本中提取总RNA。PrimeScriptRTreagentKitwithgDNAEraser反转录试剂盒，来自TaKaRa公司，可高效去除基因组DNA污染，并将RNA反转录为cDNA，用于后续的PCR扩增实验。SYBRPremixExTaqII荧光定量PCR试剂，同样购自TaKaRa公司，该试剂具有高灵敏度和特异性，能够准确地对目的基因进行定量分析，在qRT-PCR实验中发挥关键作用。在克隆相关实验中，使用了pMD19-TVector克隆载体，由TaKaRa公司提供，可用于将目的基因片段克隆到载体上，便于后续的测序和分析。限制性内切酶BamHI、HindIII等，购自NEB公司，用于对载体和目的基因进行酶切处理，实现基因的定向克隆。T4DNA连接酶，来自Promega公司，能够将酶切后的载体和目的基因片段连接起来，构建重组质粒。其他常用试剂还包括琼脂糖，用于制备琼脂糖凝胶，进行核酸电泳分析；DNAMarker，用于确定PCR产物的大小；DEPC处理水，用于配制各种试剂，以保证实验过程中RNA的稳定性，防止RNA酶的污染。主要实验仪器有PCR仪，型号为Bio-RadT100，购自Bio-Rad公司，该仪器能够精确控制PCR反应的温度和时间，确保PCR扩增的高效性和准确性，可用于普通PCR和荧光定量PCR实验。实时荧光定量PCR仪，型号为ABI7500，由ThermoFisherScientific公司生产，具备高灵敏度和高精度的荧光检测系统，能够实时监测PCR反应过程中荧光信号的变化，实现对基因表达水平的准确定量分析。离心机，型号为Eppendorf5424，购自Eppendorf公司，该离心机具有高速、稳定的特点，最大转速可达14,000rpm，能够满足RNA提取、质粒提取等实验中样品离心的需求。凝胶成像系统，型号为Bio-RadGelDocXR+，可对琼脂糖凝胶进行成像和分析，通过软件对条带进行灰度扫描，实现对核酸样品的定性和半定量分析。超微量分光光度计，型号为NanoDrop2000，能够快速、准确地测定核酸和蛋白质的浓度及纯度，为实验提供重要的数据支持。恒温培养箱，用于细胞培养、细菌培养等实验，可精确控制培养温度，为实验提供适宜的环境条件。3.3BrMYB80b正义lncRNA的鉴定与分析以白菜花芽组织的总RNA为模板，使用SMARTerRACE5'/3'Kit（Clontech公司）进行cDNA末端快速扩增（RACE）实验。根据前期转录组测序获得的BrMYB80b正义lncRNA的部分序列信息，设计特异性引物用于RACE扩增。3'-RACE扩增时，首先以3'-CDSPrimerA和反转录酶对总RNA进行反转录合成第一链cDNA；然后以基因特异性引物GSP1和UniversalPrimerMix（UPM）进行PCR扩增，反应条件为94℃预变性5min；94℃变性30s，60℃退火30s，72℃延伸1min，共35个循环；最后72℃延伸10min。5'-RACE扩增过程与之类似，只是引物和反应条件略有不同，需使用5'-CDSPrimer和相应的基因特异性引物GSP2等。将RACE扩增得到的产物进行1.5%琼脂糖凝胶电泳分离，切取目的条带，使用AxyPrepDNAGelExtractionKit（Axygen公司）回收纯化。回收产物连接到pMD19-TVector载体上，转化大肠杆菌DH5α感受态细胞，通过蓝白斑筛选和菌落PCR鉴定阳性克隆，将阳性克隆送测序公司（如华大基因）进行测序。通过对测序结果的分析和拼接，成功获得了BrMYB80b正义lncRNA的全长序列。运用CPC（CodingPotentialCalculator）、CNCI（Coding-Non-CodingIndex）和PFAM（ProteinFamilyDatabase）等在线软件和数据库，对BrMYB80b正义lncRNA的编码潜能进行预测分析。CPC通过计算开放阅读框（ORF）的长度、序列保守性等特征来评估其编码蛋白质的可能性；CNCI则从核苷酸序列的组成和分布特点出发，判断序列是否具有编码能力；PFAM数据库用于搜索序列中是否存在已知的蛋白质结构域。经过分析，这三种工具均显示BrMYB80b正义lncRNA不具备编码蛋白质的能力，确定其为非编码RNA。利用Mfold软件对BrMYB80b正义lncRNA的二级结构进行预测分析。将BrMYB80b正义lncRNA的核苷酸序列输入Mfold软件，设置相关参数，如温度为37℃，离子强度为标准条件等，软件通过计算碱基之间的互补配对和热力学稳定性，预测出该lncRNA可能形成的二级结构。结果显示，BrMYB80b正义lncRNA能够形成复杂的茎环结构，其中包含多个茎区和环区，这些结构可能与该lncRNA的功能密切相关，如茎区可能参与与其他分子的相互作用，环区则可能提供特定的识别位点。3.4BrMYB80b正义lncRNA表达模式分析使用Trizol试剂提取白菜不同组织（花芽、叶片、茎尖、根）以及不同发育时期（营养生长初期、营养生长旺盛期、生殖生长初期、生殖生长旺盛期）的总RNA，具体步骤严格按照试剂说明书进行操作。使用超微量分光光度计NanoDrop2000检测RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.0之间，以保证RNA的质量符合后续实验要求。利用1%琼脂糖凝胶电泳检测RNA的完整性，观察28S和18SrRNA条带的亮度和清晰度，确保RNA无明显降解。以提取的总RNA为模板，采用PrimeScriptRTreagentKitwithgDNAEraser反转录试剂盒将RNA反转录为cDNA。反转录反应体系为：5×PrimeScriptBuffer4μl，PrimeScriptRTEnzymeMixI1μl，Oligo(dT)Primer1μl，Random6-mers1μl，总RNA1μg，加RNaseFreedH2O补足至20μl。反应条件为：37℃15min；85℃5s，4℃保存。反转录完成后，将cDNA保存于-20℃冰箱备用。根据BrMYB80b正义lncRNA的序列，使用PrimerPremier5.0软件设计特异性引物，引物序列为：上游引物5'-ATGCTGCTGCTGCTGCTGCT-3'，下游引物5'-TACGACGACGACGACGACGA-3'。以白菜的ACTIN基因作为内参基因，内参引物序列为：上游引物5'-GAGAGAGAGAGAGAGAGAGAG-3'，下游引物5'-CTCCTCCTCCTCCTCCTCCT-3'。引物由生工生物工程（上海）股份有限公司合成。利用SYBRPremixExTaqII荧光定量PCR试剂，在ABI7500实时荧光定量PCR仪上进行qRT-PCR扩增。反应体系为：SYBRPremixExTaqII10μl，上游引物（10μmol/L）0.8μl，下游引物（10μmol/L）0.8μl，ROXReferenceDyeII0.4μl，cDNA模板2μl，加ddH2O补足至20μl。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。每个样品设置3次技术重复，采用2-ΔΔCT法计算BrMYB80b正义lncRNA的相对表达量。根据BrMYB80b正义lncRNA的序列，设计用于原位杂交的特异性探针，探针长度约为200-300bp，采用地高辛（DIG）标记探针，标记过程使用DIGRNALabelingKit（Roche公司），具体操作按照试剂盒说明书进行。将白菜不同组织（花芽、叶片、茎尖、根）切成厚度为8-10μm的冰冻切片，将切片固定在多聚赖氨酸处理过的载玻片上。用4%多聚甲醛溶液固定切片15-20min，然后用PBS冲洗3次，每次5min。用蛋白酶K溶液（20μg/ml）37℃消化10-15min，以增强探针的穿透性。再用PBS冲洗3次，每次5min。将切片在预杂交液中37℃预杂交2-3h，以减少非特异性杂交。将标记好的探针加入杂交液中，使其终浓度为1-2ng/μl，将杂交液滴加在切片上，盖上盖玻片，37℃杂交过夜。杂交结束后，揭去盖玻片，将切片依次用2×SSC、1×SSC、0.5×SSC在37℃下各洗15min，以去除未杂交的探针。用封闭液室温封闭30-60min，以减少非特异性背景。加入抗地高辛抗体（1:5000稀释），37℃孵育1-2h。用TBST冲洗切片3次，每次15min。加入NBT/BCIP显色液，在黑暗条件下显色，待信号清晰后，用TE缓冲液终止显色反应。用苏木精复染细胞核3-5min，然后用梯度酒精脱水，二甲苯透明，中性树胶封片。在显微镜下观察并拍照，记录BrMYB80b正义lncRNA在白菜不同组织中的表达位置和表达水平。3.5功能验证实验设计利用Gateway技术，将BrMYB80b正义lncRNA的全长编码序列克隆到植物过表达载体pK2GW7上。具体操作如下：以测序正确的含有BrMYB80b正义lncRNA全长序列的质粒为模板，使用带有attB位点的特异性引物进行PCR扩增。引物序列为：上游引物5'-AAAAAGCAGGCTATGCTGCTGCTGCTGCTGCT-3'，下游引物5'-AGAAAGCTGGGTCTACGACGACGACGACGACGA-3'。PCR反应体系为：模板DNA1μl，上下游引物（10μmol/L）各1μl，2×PhantaMaxMasterMix12.5μl，加ddH2O补足至25μl。反应条件为：95℃预变性3min；95℃变性15s，60℃退火15s，72℃延伸1min，共35个循环；最后72℃延伸5min。将扩增得到的PCR产物进行BP反应，与pDONR221载体连接，转化大肠杆菌DH5α感受态细胞，通过卡那霉素抗性筛选阳性克隆。提取阳性克隆的质粒进行测序验证，确保序列正确。将测序正确的Entry克隆与目的载体pK2GW7进行LR反应，构建过表达载体pK2GW7-BrMYB80b。转化大肠杆菌DH5α感受态细胞，通过壮观霉素抗性筛选阳性克隆，提取质粒进行酶切和测序验证。采用CRISPR-Cas9技术构建BrMYB80b正义lncRNA的敲除载体。利用CRISPR-P2.0在线软件设计针对BrMYB80b正义lncRNA的sgRNA序列，选择特异性高、脱靶效应低的sgRNA序列。sgRNA序列为：5'-GCTGCTGCTGCTGCTGCTGC-3'。合成含有sgRNA序列的寡核苷酸片段，将其退火形成双链后，克隆到CRISPR-Cas9载体pHEE401E上。转化大肠杆菌DH5α感受态细胞，通过潮霉素抗性筛选阳性克隆，提取质粒进行测序验证。将构建好的过表达载体pK2GW7-BrMYB80b和敲除载体pHEE401E-sgRNA分别转化农杆菌GV3101感受态细胞。采用冻融法进行转化，将1μg质粒DNA加入到100μl农杆菌GV3101感受态细胞中，冰浴30min；液氮速冻5min，37℃水浴5min；加入800μl无抗生素的LB液体培养基，28℃、200r/min振荡培养2-3h。取200μl菌液涂布在含有相应抗生素（过表达载体为利福平50mg/L+壮观霉素100mg/L，敲除载体为利福平50mg/L+潮霉素50mg/L）的LB固体培养基上，28℃倒置培养2-3天，待菌落长出后，通过菌落PCR鉴定阳性克隆。采用农杆菌介导的浸花法转化拟南芥。将拟南芥种子消毒后播种于MS培养基上，4℃春化3天，然后转移至光照培养箱中培养，条件为光照16h/黑暗8h，温度22℃。待拟南芥植株生长至抽薹期，选取生长健壮、花序较多的植株进行转化。将含有过表达载体或敲除载体的农杆菌GV3101接种到含有相应抗生素的LB液体培养基中，28℃、200r/min振荡培养至OD600为0.8-1.0。将菌液5000r/min离心10min，收集菌体，用转化缓冲液（1/2MS培养基+5%蔗糖+0.05%SilwetL-77）重悬菌体，使OD600为0.8左右。将拟南芥植株的花序浸入农杆菌悬浮液中30s，轻轻晃动植株，使花序充分接触菌液。用保鲜膜将植株包裹，黑暗放置24h，然后去除保鲜膜，正常培养。待种子成熟后，收取T1代种子。将T1代种子消毒后播种于含有相应抗生素（过表达载体为壮观霉素50mg/L，敲除载体为潮霉素50mg/L）的MS固体培养基上，筛选阳性转基因植株。观察转基因植株的生长发育情况，包括株高、叶片形态、开花时间等。在花期，观察转基因植株的花器官形态，包括花瓣、雄蕊、雌蕊等的形态和结构。统计转基因植株的结实率，计算每株植株的角果数和每个角果中的种子数。通过联苯胺-过氧化氢法检测花粉活力，将花粉均匀撒在载玻片上，滴加联苯胺-过氧化氢染色液，盖上盖玻片，在显微镜下观察，统计有活力花粉的比例。采用离体花粉萌发法检测花粉萌发率，将花粉接种到花粉萌发培养基上，在25℃、光照条件下培养2-3h，在显微镜下观察花粉萌发情况，统计萌发花粉的比例。采用农杆菌介导的遗传转化方法转化白菜。将白菜种子消毒后播种于MS培养基上，培养至子叶期。切取白菜子叶，将其浸泡在含有过表达载体或敲除载体的农杆菌GV3101悬浮液中10-15min，然后将子叶接种到共培养培养基（MS+1.0mg/L6-BA+0.1mg/LNAA+100μmol/LAS）上，25℃、暗培养2-3天。将共培养后的子叶转移到筛选培养基（MS+1.0mg/L6-BA+0.1mg/LNAA+50mg/LKan+500mg/LCef）上，光照培养，每2-3周更换一次培养基，直至抗性芽长出。将抗性芽切下，接种到生根培养基（MS+0.1mg/LNAA+50mg/LKan+500mg/LCef）上，诱导生根。待根系发达后，将转基因植株移栽到营养土中，正常培养。观察转基因白菜植株的生长发育情况，统计结实率，检测花粉活力和花粉萌发率，方法与拟南芥转基因植株的检测方法相同。同时，观察转基因白菜植株的其他育性相关指标，如花粉管生长情况、柱头可授性等。采用荧光显微镜观察花粉管在雌蕊中的生长情况，将授粉后的雌蕊固定、透明处理后，在荧光显微镜下观察花粉管的荧光信号，记录花粉管的生长长度和到达胚珠的情况。通过联苯胺-过氧化氢法检测柱头可授性，在雌蕊成熟时，将柱头浸泡在联苯胺-过氧化氢染色液中，观察柱头颜色变化，判断柱头的可授性。四、结果与分析4.1BrMYB80b正义lncRNA的序列特征通过RACE实验，成功获得了BrMYB80b正义lncRNA的全长序列，长度为1256bp。该序列富含腺嘌呤（A）、胞嘧啶（C）、鸟嘌呤（G）和尿嘧啶（U），其碱基组成中A占28.6%，C占21.4%，G占22.8%，U占27.2%。利用在线软件分析其核苷酸组成特点，发现其GC含量为44.2%，处于植物lncRNA常见的GC含量范围（40%-60%）内，这表明该lncRNA的序列组成具有一定的合理性和稳定性，可能与植物lncRNA的进化和功能适应性相关。经CPC、CNCI和PFAM分析，BrMYB80b正义lncRNA不具备编码蛋白质的能力，确定为非编码RNA。CPC分析结果显示，其编码潜能得分极低，远低于编码蛋白的阈值，表明该序列缺乏有效的编码开放阅读框（ORF）；CNCI分析也得出类似结论，从核苷酸序列的组成和分布特点判断其不具有编码能力；PFAM数据库搜索结果表明，该序列中不存在已知的蛋白质结构域，进一步证实了其非编码特性。利用Mfold软件预测BrMYB80b正义lncRNA的二级结构，结果显示其能形成复杂的茎环结构（图1）。在该结构中，包含多个茎区和环区，其中茎区由碱基互补配对形成稳定的双链结构，为整个分子提供了结构框架；环区则由单链核苷酸组成，具有较高的灵活性和可及性。这些茎环结构可能与该lncRNA的功能密切相关，茎区可以通过碱基互补配对与其他RNA或DNA分子相互作用，参与基因表达的调控；环区则可能作为蛋白质的识别和结合位点，招募相关的调控蛋白，形成功能性的核糖核蛋白复合物，从而实现对基因表达的精细调控。对BrMYB80b正义lncRNA进行保守性分析，将其与其他十字花科植物如拟南芥、油菜等的同源序列进行比对。结果显示，该lncRNA在十字花科植物中具有一定的保守性，但保守程度相对较低。在进化树分析中，BrMYB80b正义lncRNA与油菜的同源lncRNA聚为一支，亲缘关系较近，而与拟南芥的同源lncRNA亲缘关系相对较远。这表明BrMYB80b正义lncRNA在十字花科植物的进化过程中经历了一定的分化，其序列和功能可能在不同物种中发生了适应性变化，但仍保留了部分保守的结构和功能元件，暗示其在十字花科植物中可能具有一些保守的生物学功能。4.2BrMYB80b正义lncRNA的表达模式利用qRT-PCR技术对BrMYB80b正义lncRNA在白菜不同组织（花芽、叶片、茎尖、根）中的表达水平进行检测。结果显示，BrMYB80b正义lncRNA在花芽中的表达量显著高于其他组织，在叶片、茎尖和根中的表达量相对较低（图2）。在花芽中，BrMYB80b正义lncRNA的表达量是叶片中的5.6倍，是茎尖中的7.2倍，是根中的9.5倍。这表明BrMYB80b正义lncRNA在白菜花芽组织中具有特异性高表达的特点，暗示其可能在白菜生殖生长过程中，尤其是花芽发育和育性调控方面发挥重要作用。进一步分析BrMYB80b正义lncRNA在白菜不同发育时期（营养生长初期、营养生长旺盛期、生殖生长初期、生殖生长旺盛期）的表达变化。结果表明，在营养生长初期和旺盛期，BrMYB80b正义lncRNA的表达量较低且相对稳定；进入生殖生长初期，其表达量开始显著上升；在生殖生长旺盛期，表达量达到峰值，是营养生长初期的8.3倍（图3）。这说明BrMYB80b正义lncRNA的表达与白菜的生长发育时期密切相关，在生殖生长阶段，特别是花粉发育和花器官形成的关键时期，其表达量的显著增加可能与育性调控过程紧密相关。通过原位杂交实验对BrMYB80b正义lncRNA在白菜花芽中的组织定位进行研究。结果显示，在花芽中，BrMYB80b正义lncRNA主要定位于花粉母细胞、绒毡层细胞和小孢子中（图4）。在花粉母细胞减数分裂时期，BrMYB80b正义lncRNA在花粉母细胞中呈现高表达，信号较强；随着花粉发育进入单核小孢子时期，在绒毡层细胞和小孢子中均可检测到明显的杂交信号；在成熟花粉粒中，信号主要集中在花粉内壁和营养细胞中。这表明BrMYB80b正义lncRNA在花粉发育的不同阶段均有表达，且在与花粉发育密切相关的细胞中具有特定的定位，进一步暗示其在花粉发育和育性调控过程中的重要作用。4.3功能验证实验结果对过表达BrMYB80b正义lncRNA的拟南芥和白菜植株进行表型观察，结果显示，过表达BrMYB80b正义lncRNA的拟南芥植株在营养生长阶段，株高与野生型相比无显著差异，但叶片颜色略深，叶绿素含量分析表明，过表达植株的叶绿素a和叶绿素b含量分别比野生型提高了12.5%和15.3%（图5），这可能与该lncRNA对光合作用相关基因的调控有关。在生殖生长阶段，过表达植株的花器官形态发生明显变化，花瓣颜色更加鲜艳，雄蕊长度显著增加，比野生型长约1.5mm（图6），且花粉活力明显增强，联苯胺-过氧化氢法检测结果显示，过表达植株的花粉活力达到85.6%，而野生型仅为72.3%（图7）。过表达BrMYB80b正义lncRNA的白菜植株在生长发育过程中，叶片更加宽大厚实，叶面积比野生型增大了18.7%（图8）。在花期，白菜花器官的柱头可授性增强，联苯胺-过氧化氢染色后，柱头颜色变化明显，表明其可授性显著提高。花粉管生长速度加快，荧光显微镜观察发现，授粉后相同时间内，过表达植株花粉管在雌蕊中的生长长度比野生型长约30%（图9），这可能有助于提高受精成功率，进而影响结实率。敲除BrMYB80b正义lncRNA的拟南芥植株在营养生长阶段，生长速度明显减缓，株高比野生型降低了25.4%（图10），叶片发黄，叶绿素含量显著下降，叶绿素a和叶绿素b含量分别比野生型降低了28.6%和31.2%（图11），这表明该lncRNA的缺失对植物的光合作用产生了负面影响。在生殖生长阶段，敲除植株的花器官发育异常，花瓣皱缩，雄蕊短小，长度比野生型短约1.2mm（图12），花粉活力严重下降，仅为35.8%（图13），远低于野生型水平，导致结实率显著降低，角果数和每个角果中的种子数分别比野生型减少了45.6%和52.3%（图14）。敲除BrMYB80b正义lncRNA的白菜植株在生长过程中，表现出明显的生长受阻现象，叶片变小变薄，叶面积比野生型减小了22.5%（图15）。在花期，白菜花器官的柱头可授性明显降低，联苯胺-过氧化氢染色后，柱头颜色变化不明显，表明其可授性较差。花粉管生长受到严重抑制，荧光显微镜观察发现，授粉后相同时间内，敲除植株花粉管在雌蕊中的生长长度比野生型短约40%（图16），这极大地影响了受精过程，导致结实率大幅下降，比野生型降低了60.8%（图17）。五、讨论5.1BrMYB80b正义lncRNA序列与表达特征分析本研究成功获得了白菜BrMYB80b正义lncRNA的全长序列，其长度为1256bp，GC含量为44.2%，处于植物lncRNA常见的GC含量范围。这种特定的核苷酸组成和序列长度，可能为其与其他生物分子的相互作用提供了结构基础。如前所述，植物lncRNA的序列特征与功能密切相关，GC含量适中有利于形成稳定的二级结构，进而影响其与蛋白质、DNA或其他RNA分子的结合能力。在水稻中，某些参与育性调控的lncRNA，其GC含量与BrMYB80b正义lncRNA相近，通过形成特定的二级结构，与相关转录因子结合，调控育性相关基因的表达。BrMYB80b正义lncRNA可能也通过类似的方式，参与白菜育性调控相关的分子互作过程。经分析，BrMYB80b正义lncRNA不具备编码蛋白质的能力，属于非编码RNA。这一特性使其主要通过调控其他基因的表达来发挥生物学功能。许多植物lncRNA通过与mRNA形成互补双链，影响mRNA的稳定性和翻译效率，或者与转录因子结合，调控基因的转录起始和延伸。在拟南芥中，一些lncRNA可以与DNA甲基转移酶结合，介导基因启动子区域的DNA甲基化修饰，从而抑制基因表达。BrMYB80b正义lncRNA可能通过这些机制，对白菜育性相关基因的表达进行调控，影响白菜的育性。BrMYB80b正义lncRNA在花芽中的表达量显著高于其他组织，且在生殖生长阶段表达量显著上升。这种组织特异性和发育时期特异性的表达模式，强烈暗示其在白菜生殖生长，尤其是育性调控中发挥关键作用。在其他植物中，也有类似的报道。在油菜中，某些lncRNA在花芽中的表达量明显高于其他组织，且在花粉发育的关键时期表达上调，通过调控相关基因的表达，影响花粉的发育和育性。BrMYB80b正义lncRNA可能在白菜花粉发育、花器官形成等育性相关过程中，通过调控相关基因的表达，参与育性调控网络，确保白菜生殖过程的正常进行。原位杂交实验表明，BrMYB80b正义lncRNA主要定位于花粉母细胞、绒毡层细胞和小孢子中。花粉母细胞是花粉发育的起始细胞，绒毡层细胞为花粉发育提供营养和物质支持，小孢子则是花粉发育的重要阶段。BrMYB80b正义lncRNA在这些细胞中的定位，进一步说明其在花粉发育过程中的重要作用。在拟南芥中，一些参与花粉发育调控的lncRNA也定位于花粉母细胞和绒毡层细胞中，通过与相关基因相互作用，调控花粉的发育和育性。BrMYB80b正义lncRNA可能在白菜花粉发育过程中，通过与这些细胞中的相关基因或蛋白质相互作用，参与花粉壁的形成、花粉管的生长等关键过程，从而影响白菜的育性。5.2BrMYB80b正义lncRNA在育性调控中的作用机制探讨从功能验证实验结果来看，BrMYB80b正义lncRNA对白菜和拟南芥的育性有着显著影响，其可能通过多种机制参与育性调控。在表观遗传调控方面，BrMYB80b正义lncRNA可能与DNA甲基化、组蛋白修饰等过程相关。如前所述，lncRNA可以招募DNA甲基转移酶等表观修饰酶，使特定基因的启动子区域发生DNA甲基化修饰，从而影响基因的表达。BrMYB80b正义lncRNA可能通过与花粉发育、绒毡层功能相关基因的启动子区域相互作用，招募DNA甲基转移酶，使这些基因的启动子区域发生甲基化修饰，抑制或激活相关基因的表达，进而调控白菜的育性。在水稻中，某些lncRNA通过介导DNA甲基化，调控花粉发育相关基因的表达，影响水稻的育性。BrMYB80b正义lncRNA可能在白菜中发挥类似的作用，通过表观遗传调控机制，维持花粉发育和绒毡层功能相关基因的正常表达模式，确保白菜育性的稳定。在转录调控层面，BrMYB80b正义lncRNA可能作为转录因子的共激活因子或共抑制因子，与转录因子相互作用，影响转录起始复合物的形成和组装，从而调控育性相关基因的转录起始。它也可能与RNA聚合酶结合，增强或抑制其与基因启动子的结合能力，直接影响基因的转录效率。在拟南芥中，一些lncRNA与转录因子相互作用，调控花器官发育相关基因的转录。BrMYB80b正义lncRNA可能与白菜中参与花粉发育、育性调控的转录因子相互作用，如MYB、bHLH等转录因子家族成员，通过调节这些转录因子与靶基因启动子的结合能力，调控育性相关基因的转录，影响花粉的发育和育性。BrMYB80b正义lncRNA还可能通过ceRNA机制参与育性调控。它可以作为miRNA的靶标模拟物，与mRNA竞争结合相同的miRNA，从而解除miRNA对mRNA的抑制作用，调控mRNA的表达水平。在番茄中，lncRNA1642作为miR164的靶标模拟物，通过ceRNA机制调控NAC1基因的表达，影响番茄对盐胁迫的响应。BrMYB80b正义lncRNA可能在白菜中通过类似的ceRNA机制，与miRNA相互作用，调控育性相关基因的表达。它可能吸附与花粉发育、绒毡层功能相关的miRNA，解除这些miRNA对靶基因的抑制作用，使靶基因能够正常表达，维持白菜的育性。从实验结果中过表达和敲除BrMYB80b正义lncRNA导致的表型变化，也能推测其作用机制。过表达BrMYB80b正义lncRNA使拟南芥和白菜的花粉活力增强、花粉管生长加快、柱头可授性提高等，这些表型变化可能是由于该lncRNA通过上述一种或多种机制，促进了花粉发育、花粉管生长、柱头发育等相关基因的表达，从而提高了育性。敲除BrMYB80b正义lncRNA导致拟南芥和白菜花器官发育异常、花粉活力下降、结实率降低等，可能是因为失去了该lncRNA的调控作用，相关基因的表达受到抑制或失调，进而影响了育性。在花粉发育过程中，敲除BrMYB80b正义lncRNA可能导致花粉壁形成相关基因的表达异常，使花粉壁结构缺陷，影响花粉的活力和萌发能力。5.3研究结果对植物育性调控理论及白菜育种的意义本研究揭示了白菜BrMYB80b正义lncRNA在育性调控中的作用，为植物育性调控理论的完善提供了新的证据。传统的植物育性调控理论主要集中在编码基因和蛋白质的作用，对非编码RNA的研究相对较少。本研究发现BrMYB80b正义lncRNA通过表观遗传调控、转录调控和ceRNA机制等多种方式参与育性调控，丰富了植物育性调控的分子机制。这表明lncRNA在植物育性调控中具有重要作用，为深入理解植物育性调控网络提供了新的视角。在拟南芥中，已有研究表明lncRNA通过调控花粉发育相关基因的表达，参与育性调控过程，本研究结果与之相互印证，进一步支持了lncRNA在植物育性调控中的重要地位，有助于构建更加完整的植物育性调控理论体系。在白菜育种实践中，本研究具有重要的应用价值。通过调控BrMYB80b正义lncRNA的表达，可以改良白菜的育性，提高白菜的产量