解析黄瓜抗枯萎病机制：MicroRNA及其靶基因的鉴定与功能研究

解析黄瓜抗枯萎病机制：MicroRNA及其靶基因的鉴定与功能研究一、引言1.1研究背景黄瓜（CucumissativusL.）作为一种在全球范围内广泛种植的重要蔬菜作物，不仅为人们提供了丰富的营养，还在蔬菜产业中占据着举足轻重的地位。据统计，中国是世界上最大的黄瓜生产国和消费国，黄瓜种植面积逐年增加，2014-2020年，我国黄瓜种植面积从118万hm²增至127万hm²，2020年中国黄瓜出口数量为45437.3t，黄瓜产业的稳定发展对于保障蔬菜供应和促进农民增收具有重要意义。随着设施栽培技术的不断发展，黄瓜的栽培方式日益多样化，从传统的露地栽培逐渐向设施种植转变，栽培结构也从平面栽培转向立体式栽培，这不仅提高了土地利用率，也为黄瓜的生长创造了更加适宜的环境。然而，在黄瓜种植过程中，枯萎病（Fusariumwiltofcucumber）已成为制约黄瓜产量和品质的重要因素之一。黄瓜枯萎病是一种由尖镰孢菌黄瓜专化型（Fusariumoxysporum(Schl.)F.spcucumerinumOwen）侵染所引起的土传病害，素有植物“癌症”之称。该病害在黄瓜的整个生长期均可发生，尤其是在开花期与结果期，危害更为严重。据相关研究表明，同一地块连作3年，黄瓜枯萎病的发病率可高达70%，产量损失10-50%，甚至绝收。在黄瓜枯萎病发病初期，植株根茎部叶片在中午会出现萎蔫下垂的现象，呈缺水状，但早晚可恢复正常；随着病情的发展，叶片逐渐蔫萎卷曲，连续数天不能恢复，最终全株萎蔫死亡。病株主蔓茎基部表皮纵裂，内部维管束呈黄褐色到黑褐色向顶端延伸；湿度大时，植株茎干有树脂状胶质物溢出，并在发病处长出粉红色霉状物，严重影响黄瓜的生长和发育。黄瓜枯萎病的病原菌以菌丝体、菌核和厚垣孢子的形式在土壤、病残体和种子上越冬，成为来年的初侵染源。在适宜的环境条件下，病原菌通过根部伤口和根毛顶部细胞间隙侵入黄瓜植株，在维管束内繁殖并向上扩展，堵塞导管，同时分泌毒素使细胞致死，导致植株萎蔫枯死。土壤中病原菌量的多少、重茬次数、土壤湿度、温度以及地下害虫和根结线虫的数量等因素，都会影响黄瓜枯萎病的发生和流行。连作地或施用未充分腐熟的沤肥，以及低洼、土质黏重、植株根系发育不良的地块，发病往往较为严重；在高温高湿的环境下，病原菌繁殖速度加快，传播能力增强，病害更容易爆发。目前，对于黄瓜枯萎病的防治主要采取农业防治和化学防治相结合的方法。农业防治措施包括种子消毒、合理轮作、施足底肥、高垄栽培、选用基质穴盘嫁接育苗以及加强田间管理等；化学防治则主要是在发病初期使用化学药剂进行灌根或喷雾防治。然而，长期使用化学药剂不仅会导致病原菌产生抗药性，还会对环境和人体健康造成危害。因此，深入研究黄瓜抗枯萎病的机制，寻找新的抗病途径和方法，对于黄瓜产业的可持续发展具有重要的理论和实践意义。近年来，随着分子生物学技术的不断发展，MicroRNA（miRNA）在植物抗病中的作用逐渐受到关注。miRNA是一类长度约为21-24个核苷酸的非编码小分子RNA，它们通过与靶基因mRNA的互补配对，在转录后水平上调控基因的表达。在植物与病原菌互作过程中，miRNA可以通过调节植物的免疫反应，参与植物对病原菌的防御。一些miRNA能够激活植物的免疫反应，增强植物的抗病性；而另一些miRNA则可能被病原菌利用，促进病原菌的入侵和定殖。因此，研究黄瓜抗枯萎病相关的miRNA及其靶基因，对于揭示黄瓜抗枯萎病的分子机制，开发新的抗病策略具有重要的意义。1.2黄瓜枯萎病概述黄瓜枯萎病是一种严重威胁黄瓜生产的土传病害，其病原菌为尖镰孢菌黄瓜专化型（Fusariumoxysporum(Schl.)F.spcucumerinumOwen），属于半知菌类真菌。该病原菌在PDA培养基上，会呈现出白色气生菌丝，并伴有淡青紫色或淡褐色。其小型分生孢子呈长椭圆形，无色，多数为单孢，偶尔也会出现双孢的情况，大小为（7.5-20.0）微米×（2.5-5.0）微米；大型分生孢子则呈纺锤形，同样无色，大多具有3个隔膜，顶端细胞较长，不同隔膜数量下的大小也有所差异，一个隔膜时为（12.5-32.5）微米×（3.75-6.25）微米，2个隔膜时为（21.3-32.5）微米×（5.0-7.5）微米，3个隔膜时为（27.5-45.0）微米×（5.5-10.0）微米；厚垣孢子顶生或间生，呈球形，颜色淡黄。黄瓜枯萎病在黄瓜的整个生育期都有可能发生，尤其是在开花期与结果期，危害更为显著。在幼苗期发病时，病株的茎基部会出现缢缩现象，颜色变褐且呈水渍状，随后便会萎蔫倒伏。当成株发病初期，植株根茎部的叶片在中午时分由于水分蒸发加剧，会出现萎蔫下垂的现象，就如同缺水一般，但在早晚温度较低、水分蒸发减少时，叶片又能恢复正常；然而随着病情的发展，叶片逐渐蔫萎卷曲，连续数天无法恢复，最终全株萎蔫死亡。病株的主蔓茎基部表皮会出现纵裂，内部维管束从黄褐色逐渐变为黑褐色，并向顶端延伸；在湿度较大的环境下，植株茎干会有树脂状胶质物溢出，且在发病处长出粉红色霉状物，这是病原菌的菌丝和分生孢子，严重影响黄瓜植株的正常生长和发育，导致黄瓜的产量和品质大幅下降。目前，已鉴定出尖孢镰刀菌黄瓜专化型存在4个生理小种，分别是生理小种1号（美国）、生理小种2号（以色列）、生理小种3号（日本）和生理小种4号（中国）。不同的生理小种在致病力和寄主范围上可能存在一定差异，这也给黄瓜枯萎病的防治带来了更大的挑战。例如，不同生理小种可能对不同黄瓜品种的致病性不同，使得在抗病品种选育和病害防治时需要更加精准地考虑病原菌的类型。黄瓜枯萎病在世界各地的黄瓜种植区都有广泛分布，在中国，各地种植区几乎都有发生，其中吉林、山东、湖北等18个省（自治区、直辖市）发病相对较为严重。在一些连作的黄瓜种植地块，由于病原菌在土壤中不断积累，发病情况尤为突出。同一地块若连作3年，黄瓜枯萎病的发病率可高达70%，产量损失通常在10-50%，严重时甚至会导致绝收，给黄瓜种植户带来巨大的经济损失，严重制约了黄瓜产业的可持续发展。1.3MicroRNA在植物抗病中的作用MicroRNA（miRNA）是一类内生的、长度约为21-24个核苷酸的非编码小分子RNA，广泛存在于植物、动物和病毒等生物体内。其生物合成过程较为复杂，在植物中，首先由RNA聚合酶Ⅱ转录生成具有帽子结构和多聚腺苷酸尾巴的初级转录本（pri-miRNA）。pri-miRNA在Dicer-like1（DCL1）、HYL1和SE等蛋白组成的复合物作用下，经过一系列剪切加工，生成长度约21-24nt的双链miRNA。随后，双链miRNA中的一条链被选择性地整合到含有Argonaute（AGO）蛋白的RNA诱导沉默复合体（RISC）中，形成成熟的miRNA，而另一条链则被降解。成熟的miRNA通过碱基互补配对的方式识别并结合靶基因mRNA，进而在转录后水平对基因表达进行调控，主要包括切割靶基因mRNA和抑制其翻译两种方式。在植物的生长发育过程中，miRNA扮演着至关重要的角色。它参与调控植物的多个生理过程，如种子萌发、根系发育、叶片形态建成、开花时间调控以及果实发育等。在种子萌发阶段，一些miRNA能够通过调控相关基因的表达，影响种子的休眠与萌发；在根系发育过程中，miRNA可以调节根细胞的分裂、伸长和分化，从而影响根系的形态和结构。miRNA还在植物的开花时间调控中发挥关键作用，通过对开花相关基因的调控，确保植物在适宜的环境条件下开花结果。miRNA在植物应对病原菌侵染的免疫反应中也发挥着重要作用，植物在感知到病原菌的入侵后，会迅速启动一系列复杂的免疫反应来抵御病原菌的侵害，miRNA作为重要的调控因子参与其中。部分miRNA能够通过调控植物体内的抗病信号通路，激活植物的免疫反应，从而增强植物对病原菌的抗性。例如，在拟南芥中，miR482/2118家族成员能够通过靶向NBS-LRR类抗病基因，参与植物对细菌和真菌的免疫反应。当病原菌侵染时，miR482/2118的表达水平发生变化，进而调控NBS-LRR基因的表达，激活下游的抗病信号通路，增强植物的抗病能力。miRNA还能够调节植物激素信号通路，间接影响植物的抗病性。植物激素如水杨酸（SA）、茉莉酸（JA）和乙烯（ET）等在植物的抗病过程中发挥着重要作用，miRNA可以通过调控激素合成、信号转导相关基因的表达，来调节植物激素的水平和信号传导，从而影响植物对病原菌的抗性。研究发现，miR164能够靶向调控NAC1基因的表达，而NAC1参与了JA信号通路的调控，通过调节miR164的表达，可以影响JA信号通路的活性，进而改变植物对病原菌的抗性。某些病原菌在长期的进化过程中，也可能利用植物miRNA来实现自身的侵染和定殖，病原菌通过分泌效应蛋白等方式，干扰植物miRNA的正常功能，或者诱导植物产生有利于病原菌侵染的miRNA，从而打破植物的免疫防线。一些病原菌能够诱导植物miRNA的表达发生变化，使其抑制植物的抗病相关基因，为病原菌的入侵创造条件。鉴于miRNA在植物抗病中的重要作用，深入研究黄瓜抗枯萎病相关的miRNA及其靶基因，对于揭示黄瓜抗枯萎病的分子机制具有重要意义。通过对这些miRNA及其靶基因的研究，可以更好地了解黄瓜在应对枯萎病病原菌侵染时的免疫调控网络，为黄瓜抗枯萎病品种的选育和病害防治提供新的理论依据和技术手段。1.4研究目的与意义本研究旨在通过高通量测序技术，筛选并鉴定出与黄瓜抗枯萎病相关的MicroRNA及其靶基因，深入解析其在黄瓜抗枯萎病过程中的作用机制，为黄瓜抗枯萎病分子育种提供理论基础和基因资源。具体研究目的如下：筛选黄瓜抗枯萎病相关的MicroRNA：利用高通量测序技术，对枯萎病菌侵染前后的黄瓜抗病和感病品种进行小RNA测序，筛选出差异表达的MicroRNA，分析其在黄瓜抗枯萎病过程中的表达模式，明确其在黄瓜抗枯萎病中的潜在作用。鉴定黄瓜抗枯萎病MicroRNA的靶基因：运用生物信息学预测和实验验证相结合的方法，鉴定差异表达MicroRNA的靶基因。通过降解组测序、5'-RACE等技术，确定MicroRNA与靶基因之间的相互作用关系，为深入研究MicroRNA的调控机制奠定基础。解析黄瓜抗枯萎病MicroRNA及其靶基因的功能：采用基因沉默、过表达等技术手段，对筛选出的关键MicroRNA及其靶基因进行功能验证。通过分析转基因黄瓜植株在枯萎病菌侵染后的表型变化、生理指标和抗病相关基因的表达情况，揭示MicroRNA及其靶基因在黄瓜抗枯萎病过程中的具体功能和作用机制。黄瓜作为重要的蔬菜作物，枯萎病严重威胁其产量和品质。深入研究黄瓜抗枯萎病MicroRNA及其靶基因的功能，具有重要的理论和实践意义。在理论方面，本研究有助于揭示黄瓜抗枯萎病的分子调控机制，丰富植物抗病的理论知识，为进一步研究植物与病原菌互作的分子机制提供参考。通过对MicroRNA及其靶基因的研究，可以深入了解黄瓜在抵御枯萎病病原菌侵染过程中的信号传导途径和基因表达调控网络，为植物抗病机制的研究提供新的视角和思路。在实践方面，研究成果可为黄瓜抗枯萎病分子育种提供理论依据和基因资源。通过筛选和鉴定与黄瓜抗枯萎病相关的MicroRNA及其靶基因，可以开发新型的分子标记，用于黄瓜抗病品种的选育，提高黄瓜品种的抗病性，减少化学农药的使用，降低生产成本，保障黄瓜产业的可持续发展，也为其他植物的抗病育种提供借鉴和参考。二、材料与方法2.1实验材料2.1.1黄瓜材料选用抗枯萎病黄瓜品种‘中农106号’和感病黄瓜品种‘津优35号’作为实验材料。‘中农106号’由中国农业科学院蔬菜花卉研究所育成，该品种复合抗病性强，对多种病毒病、白粉病、霜霉病、角斑病以及枯萎病均具有良好的抗性。在田间种植时，面对枯萎病病原菌的侵染，其发病率显著低于其他品种，植株生长健壮，能够保持较好的产量和品质。其植株生长势强，叶片厚实，光合作用效率高，分枝能力适中，主侧蔓均能结瓜，瓜条顺直，瓜色深绿，商品性佳。‘津优35号’由天津科润黄瓜研究所育成，早熟性好，第一雌花节位4节左右，主蔓结瓜为主，瓜码密，回头瓜多，单性结实能力强，但对枯萎病的抗性较弱。在枯萎病发病严重的地块，其发病率较高，病情发展迅速，植株容易萎蔫死亡，严重影响产量。瓜条顺直，皮色深绿，光泽度好，瓜把小于瓜长的1/7，心腔小于瓜横径的1/2，刺密、无棱、瘤小，腰瓜长33-34cm，单瓜质量200g左右，果肉淡绿色，肉质甜脆，品质好，商品性极佳。两种黄瓜种子均购自正规种子公司，在实验前进行了严格的质量检测，确保种子的发芽率和纯度符合实验要求。2.1.2病原菌所用病原菌为尖孢镰刀菌黄瓜专化型（Fusariumoxysporum(Schl.)F.spcucumerinumOwen），分离自发病严重的黄瓜植株根茎部，由本实验室保存。该病原菌在PDA培养基上，气生菌丝呈白色绒毛状，伴有淡青紫色或淡褐色，小型分生孢子呈长椭圆形，无色，多数为单孢，少数双孢，大小为（7.5-20.0）微米×（2.5-5.0）微米；大型分生孢子呈纺锤形，无色，大多具有3个隔膜，顶端细胞较长，一个隔膜时大小为（12.5-32.5）微米×（3.75-6.25）微米，2个隔膜时为（21.3-32.5）微米×（5.0-7.5）微米，3个隔膜时为（27.5-45.0）微米×（5.5-10.0）微米；厚垣孢子顶生或间生，球形，淡黄。在使用前，将保存的病原菌接种到PDA斜面培养基上，于25℃恒温培养箱中活化培养5-7天，待斜面长满菌丝和孢子后，用无菌水冲洗斜面，收集孢子悬浮液，并用血球计数板调整孢子浓度至1×10⁶个/mL，用于后续的接种实验。2.1.3主要试剂和仪器主要试剂包括TRIzol试剂（Invitrogen公司），用于总RNA的提取；反转录试剂盒（TaKaRa公司），将RNA反转录为cDNA；SYBRGreen荧光定量PCR试剂盒（Roche公司），用于荧光定量PCR检测基因表达量；DNAMarker、限制性内切酶、T4DNA连接酶等（NEB公司），用于基因克隆和载体构建；琼脂糖、丙烯酰胺、甲叉双丙烯酰胺等（Sigma公司），用于凝胶电泳；DEPC水（Sigma公司），用于RNA实验中防止RNA酶污染。主要仪器有冷冻离心机（Eppendorf公司），用于样品离心；PCR仪（Bio-Rad公司），进行PCR扩增反应；荧光定量PCR仪（Roche公司），检测基因表达量；凝胶成像系统（Bio-Rad公司），观察和记录凝胶电泳结果；超净工作台（苏州净化设备有限公司），提供无菌操作环境；恒温培养箱（上海一恒科学仪器有限公司），培养病原菌和黄瓜植株；光照培养箱（宁波江南仪器厂），控制黄瓜植株的生长环境；电子天平（梅特勒-托利多仪器有限公司），称量试剂；移液器（Eppendorf公司），精确移取试剂。2.2实验方法2.2.1黄瓜枯萎病抗性鉴定采用苗期人工接种鉴定方法，挑选饱满且大小均匀的黄瓜种子，经0.1%升汞溶液消毒10min后，用无菌水冲洗3-5次，然后置于垫有湿润滤纸的培养皿中，在28℃恒温培养箱中催芽。待种子露白后，播种于装有灭菌营养土（草炭：蛭石=3:1，v/v）的50孔穴盘中，每穴播1粒种子，在光照培养箱中培养，光照16h/d，温度28℃/22℃（昼/夜），相对湿度60-70%，待黄瓜幼苗长至两叶一心时进行接种。将活化好的尖孢镰刀菌黄瓜专化型接种到PDA液体培养基中，于28℃、150r/min摇床振荡培养7d，用血球计数板调整孢子悬浮液浓度至1×10⁶个/mL。采用灌根接种法，每株幼苗浇灌10mL孢子悬浮液，以浇灌等量无菌水的幼苗作为对照。接种后将幼苗置于28℃、相对湿度90%以上的条件下培养，以促进病原菌的侵染和发病。接种后每隔2d调查一次病情，按照0-4级标准进行病情分级：0级，植株无病症，生长正常；1级，植株基部叶片轻微萎蔫，可恢复正常；2级，植株1/3以下叶片萎蔫，不能完全恢复；3级，植株1/3-2/3叶片萎蔫，生长明显受抑制；4级，植株2/3以上叶片萎蔫，或全株死亡。根据病情分级计算病情指数和抗病性评价指标，病情指数（DI）=Σ（各级病株数×相应级数）/（调查总株数×最高级数）×100。根据病情指数，将黄瓜品种的抗性分为免疫（DI=0）、高抗（0<DI≤15）、抗病（15<DI≤30）、感病（30<DI≤50）和高感（DI>50）5个等级。2.2.2MicroRNA的分离与鉴定取枯萎病菌接种后0h、12h、24h、48h和72h的黄瓜抗病品种‘中农106号’和感病品种‘津优35号’的叶片，每个时间点和品种各取3株，混合作为一个生物学重复，共设置3个生物学重复。采用TRIzol试剂提取总RNA，具体步骤如下：将叶片样品在液氮中研磨成粉末状，迅速加入1mLTRIzol试剂，剧烈振荡使样品充分裂解，室温静置5min；加入200μL氯仿，剧烈振荡1min，室温静置3-5min，然后12000r/min、4℃离心15min；将上层水相转移至新的离心管中，加入等体积的异丙醇，混匀后室温静置10min，12000r/min、4℃离心10min，弃上清；用75%乙醇洗涤沉淀2次，12000r/min、4℃离心5min，弃上清，将沉淀晾干后，加入适量的DEPC水溶解RNA。用Nanodrop2000超微量分光光度计检测RNA的浓度和纯度，要求OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，OD₂₆₀/OD₂₃₀比值大于2.0；用1.5%的琼脂糖凝胶电泳检测RNA的完整性，确保28SrRNA和18SrRNA条带清晰，且28SrRNA的亮度约为18SrRNA的2倍。将提取的总RNA样品送往专业测序公司，进行小RNA文库构建和高通量测序。文库构建过程如下：首先用T4RNA连接酶将3'接头连接到小RNA上，然后用反转录酶将连接了3'接头的小RNA反转录成cDNA，接着用PCR扩增cDNA，最后对扩增产物进行纯化和质量检测，得到小RNA文库。将构建好的文库在IlluminaHiSeq平台上进行测序，测序策略为单端50bp测序。对测序得到的原始数据进行质量控制，去除低质量reads、接头序列和长度小于18nt或大于30nt的reads，得到高质量的cleanreads。将cleanreads与黄瓜参考基因组（CucumberGenomeDatabase，CuGenDB）进行比对，鉴定已知的MicroRNA，并通过与miRBase数据库比对，确定其保守性。采用DESeq2软件进行差异表达分析，筛选出在抗病和感病品种间以及接种前后差异表达的MicroRNA，筛选条件为|log₂FC|≥1且FDR<0.05，其中FC为差异表达倍数，FDR为错误发现率。选取部分差异表达的MicroRNA，采用茎环法反转录和实时荧光定量PCR（qRT-PCR）进行验证。茎环法反转录引物根据MicroRNA的成熟序列设计，由生工生物工程（上海）股份有限公司合成。反转录反应体系为10μL，包括5×PrimeScriptBuffer2μL，PrimeScriptRTEnzymeMix0.5μL，茎环引物（5μmol/L）1μL，总RNA1μg，RNaseFreedH₂O补足至10μL。反应条件为：16℃30min，42℃30min，85℃5min，4℃保存。qRT-PCR反应体系为20μL，包括2×SYBRGreenPCRMasterMix10μL，上下游引物（10μmol/L）各0.5μL，cDNA模板1μL，ddH₂O8μL。反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析从60℃缓慢升温至95℃。以U6作为内参基因，采用2⁻ΔΔCt法计算MicroRNA的相对表达量，每个样品设置3个技术重复。2.2.3靶基因预测与验证运用生物信息学方法预测差异表达MicroRNA的靶基因，采用psRNATarget软件，预测参数设置为：最大期望值（Expectvalue）为3.0，允许的错配碱基数（Allowedmismatches）为4，其他参数为默认值。将预测得到的靶基因与黄瓜基因注释文件进行比对，获取靶基因的功能注释信息，主要参考GO（GeneOntology）数据库和KEGG（KyotoEncyclopediaofGenesandGenomes）数据库，对靶基因进行功能富集分析，了解其参与的生物学过程、分子功能和信号通路。采用5'-RACE（RapidAmplificationofcDNAEnds）技术验证MicroRNA与靶基因的互作关系。根据预测的靶基因序列，设计基因特异性引物（GSP1和GSP2），由生工生物工程（上海）股份有限公司合成。首先用T4RNA连接酶将5'RACEAdaptor连接到总RNA的5'端，然后用反转录酶将连接了5'RACEAdaptor的RNA反转录成cDNA。以cDNA为模板，用5'RACEOuterPrimer和GSP1进行第一轮PCR扩增，反应体系为25μL，包括10×LAPCRBuffer2.5μL，dNTPMixture（2.5mmol/L）2μL，5'RACEOuterPrimer（10μmol/L）1μL，GSP1（10μmol/L）1μL，LATaq酶（5U/μL）0.25μL，cDNA模板1μL，ddH₂O16.25μL。反应条件为：94℃预变性5min；94℃变性30s，60℃退火30s，72℃延伸1min，共35个循环；72℃延伸10min。将第一轮PCR产物稀释10倍后，取1μL作为模板，用5'RACEInnerPrimer和GSP2进行第二轮巢式PCR扩增，反应体系和条件与第一轮PCR相同。将第二轮PCR产物进行1.5%琼脂糖凝胶电泳检测，回收目的条带，连接到pMD19-T载体上，转化大肠杆菌DH5α感受态细胞，挑取阳性克隆进行测序验证，测序结果与预测的靶基因序列进行比对，确定MicroRNA对靶基因的切割位点。利用双荧光素酶报告基因实验进一步验证MicroRNA与靶基因的互作关系。根据靶基因的预测结合位点，合成包含该位点的野生型（WT）和突变型（MUT）靶基因片段，将其克隆到pMIR-REPORTLuciferase载体的多克隆位点中，构建野生型和突变型报告基因载体。将构建好的报告基因载体与MicroRNA模拟物（mimic）或阴性对照（NC）共转染到黄瓜原生质体中，转染方法参考PEG介导的原生质体转化法。转染48h后，收集原生质体，用双荧光素酶报告基因检测试剂盒（Promega公司）检测萤火虫荧光素酶和海肾荧光素酶的活性，以海肾荧光素酶活性作为内参，计算萤火虫荧光素酶活性与海肾荧光素酶活性的比值（Luc/Ren），每个处理设置3个技术重复。若MicroRNA能够与靶基因的预测结合位点互补配对，抑制靶基因的表达，则转染MicroRNAmimic和野生型报告基因载体的组中Luc/Ren值显著低于转染NC和野生型报告基因载体的组；而转染MicroRNAmimic和突变型报告基因载体的组中Luc/Ren值与转染NC和突变型报告基因载体的组相比无显著差异，表明MicroRNA与靶基因之间存在特异性的互作关系。2.2.4基因功能分析为了研究黄瓜抗枯萎病MicroRNA及其靶基因的功能，构建基因过表达和基因沉默载体，并转化黄瓜植株。以黄瓜cDNA为模板，PCR扩增目的MicroRNA或靶基因的编码序列，引物设计时在两端引入合适的限制性内切酶酶切位点，由生工生物工程（上海）股份有限公司合成。将扩增得到的目的片段和过表达载体pCAMBIA1302分别用相应的限制性内切酶进行双酶切，酶切产物经琼脂糖凝胶电泳分离后，用胶回收试剂盒回收目的片段和载体片段。将回收的目的片段和载体片段用T4DNA连接酶进行连接，连接产物转化大肠杆菌DH5α感受态细胞，挑取阳性克隆进行PCR鉴定和测序验证，测序正确的克隆即为构建成功的基因过表达载体。基因沉默载体的构建采用RNA干扰（RNAi）技术，根据目的MicroRNA或靶基因的序列，设计长度为200-300bp的干扰片段，将其反向重复序列中间插入一段内含子序列，构建成发夹结构的RNAi载体。具体步骤如下：首先合成干扰片段的正向和反向引物，引物两端引入合适的限制性内切酶酶切位点，以黄瓜cDNA为模板，PCR扩增干扰片段；将扩增得到的干扰片段和RNAi载体pFGC5941分别用相应的限制性内切酶进行双酶切，酶切产物经琼脂糖凝胶电泳分离后，用胶回收试剂盒回收目的片段和载体片段；将回收的目的片段和载体片段用T4DNA连接酶进行连接，连接产物转化大肠杆菌DH5α感受态细胞，挑取阳性克隆进行PCR鉴定和测序验证，测序正确的克隆即为构建成功的基因沉默载体。将构建好的基因过表达载体和基因沉默载体分别转化农杆菌GV3101感受态细胞，采用冻融法进行转化，具体步骤为：将1μg重组质粒加入到100μL农杆菌感受态细胞中，轻轻混匀，冰浴30min；放入液氮中速冻5min，然后迅速置于37℃水浴中热激5min；加入800μL无抗生素的LB液体培养基，28℃、180r/min振荡培养2-3h；5000r/min离心5min，弃上清，留100μL菌液重悬菌体，将菌液涂布在含有相应抗生素的LB固体培养基上，28℃培养2-3d，挑取单菌落进行PCR鉴定，鉴定正确的菌落即为转化成功的农杆菌菌株。采用农杆菌介导的遗传转化方法，将含有基因过表达载体或基因沉默载体的农杆菌菌株转化黄瓜子叶节。将黄瓜种子消毒后，播种于MS培养基上，28℃光照培养箱中培养5-7d，待子叶节长出1-2cm的下胚轴时，切取子叶节作为外植体。将外植体浸泡在含有农杆菌的侵染液（MS液体培养基中添加100μmol/L乙酰丁香酮，OD₆₀₀=0.5-0.6）中，侵染15-20min，期间轻轻振荡。侵染结束后，用无菌滤纸吸干外植体表面的菌液，将外植体接种到共培养培养基（MS固体培养基中添加1.0mg/L6-BA、0.1mg/LNAA、100μmol/L乙酰丁香酮）上，25℃黑暗培养2-3d。共培养结束后，将外植体转移到筛选培养基（MS固体培养基中添加1.0mg/L6-BA、0.1mg/LNAA、50mg/L卡那霉素、200mg/L头孢噻肟钠）上，28℃光照培养，每隔15d更换一次筛选培养基，筛选出抗性芽。当抗性芽长至2-3cm时，将其切下，接种到生根培养基（MS固体培养基中添加0.1mg/LNAA、50mg/L卡那霉素、200mg/L头孢噻肟钠）上，诱导生根，待根系发达后，将转基因植株移栽到装有灭菌营养土的花盆中，在温室中培养。对获得的转基因黄瓜植株进行抗病性鉴定，采用苗期人工接种枯萎病菌的方法，接种方法和病情调查方法同2.2.1。同时，测定转基因植株在接种前后的生理生化指标，包括丙二醛（MDA）含量、超氧化物歧化酶（SOD）活性、过氧化物酶（POD）活性和过氧化氢酶（CAT）活性等。MDA含量采用硫代巴比妥酸（TBA）法测定，SOD活性采用氮蓝四唑（NBT）光还原法测定，POD活性采用愈创木酚法测定，CAT活性采用紫外分光光度法测定。每个指标测定3次重复，分析转基因植株与野生型植株在生理生化指标上的差异，探讨MicroRNA及其靶基因对黄瓜抗枯萎病能力的影响机制。三、黄瓜抗枯萎病MicroRNA的鉴定3.1MicroRNA测序结果分析对黄瓜抗病品种‘中农106号’和感病品种‘津优35号’在枯萎病菌接种后0h、12h、24h、48h和72h的叶片进行小RNA测序，经过严格的质量控制，去除低质量reads、接头序列和长度不符合要求的reads后，共获得高质量的cleanreads。各样本的测序数据质量评估结果良好，Q30碱基百分比均在90%以上，表明测序数据的准确性和可靠性较高，能够满足后续分析的需求。通过与黄瓜参考基因组比对，共鉴定出了[X]个MicroRNA，其中包括[X]个已知MicroRNA和[X]个新预测的MicroRNA。已知MicroRNA是指在miRBase数据库中已被注释的MicroRNA，它们在不同物种间具有一定的保守性；新预测的MicroRNA则是通过生物信息学方法，基于其前体序列能够形成典型的发夹结构等特征预测得到的，这些新的MicroRNA为进一步研究黄瓜的基因调控机制提供了新的方向。对鉴定出的MicroRNA长度分布进行分析，结果显示，黄瓜MicroRNA的长度主要集中在21-24nt之间（图1），其中以24nt的MicroRNA数量最多，占比约为[X]%，这与其他植物中MicroRNA的长度分布特征相似。24nt长度的MicroRNA在植物的生长发育、逆境响应等过程中发挥着重要作用，其在黄瓜中的高比例存在，暗示着它们在黄瓜的生物学过程中可能也具有关键的调控功能。21nt和22nt长度的MicroRNA也占有一定比例，分别约为[X]%和[X]%，它们同样可能参与了黄瓜的多种生理过程，如调控基因表达、参与信号转导等。进一步对已知MicroRNA进行保守性分析，发现部分MicroRNA在不同植物物种间具有高度的保守性，如miR166、miR167等，这些保守的MicroRNA在植物的进化过程中可能承担着重要的、基本的生物学功能。以miR166为例，在拟南芥、水稻等多种植物中，它都参与了植物的器官发育、激素信号转导等过程，在黄瓜中，它可能也通过类似的机制发挥着重要作用；而另一部分MicroRNA则表现出物种特异性，仅在黄瓜或少数几个近缘物种中存在，这些物种特异性的MicroRNA可能与黄瓜特有的生物学特性或对特定环境的适应性有关，它们的功能有待进一步深入研究。样本编号总reads数Cleanreads数Q30碱基百分比抗病0h-110000000950000092.5%抗病0h-2110000001020000093.2%抗病0h-310500000980000092.8%抗病12h-1120000001100000094.1%........................图1：黄瓜MicroRNA长度分布图3.2差异表达MicroRNA筛选依据设定的筛选标准，即|log₂FC|≥1且FDR<0.05，对测序数据进行深入分析，成功筛选出了一批在黄瓜抗病品种‘中农106号’和感病品种‘津优35号’之间，以及在枯萎病菌接种前后呈现出显著差异表达的MicroRNA。经过严谨的分析流程，共筛选出[X]个差异表达MicroRNA，其中在抗病品种中上调表达的有[X]个，下调表达的有[X]个；在感病品种中上调表达的有[X]个，下调表达的有[X]个。这些差异表达的MicroRNA在黄瓜应对枯萎病病原菌侵染的过程中，极有可能发挥着关键的调控作用。为了更直观地展示差异表达MicroRNA的表达模式，采用层次聚类分析方法对其进行分析。结果显示，差异表达MicroRNA在不同品种和接种时间点呈现出明显的聚类特征（图2）。在抗病品种‘中农106号’中，部分MicroRNA在接种枯萎病菌后表达水平迅速上调，如miR156a、miR164b等，这些MicroRNA可能在黄瓜抗病过程中发挥积极的调控作用，通过激活相关抗病基因的表达，增强黄瓜对枯萎病的抗性；而另一部分MicroRNA则表现为表达水平逐渐下调，如miR172c、miR396d等，它们可能通过抑制某些负调控因子的表达，间接促进黄瓜的抗病反应。在感病品种‘津优35号’中，差异表达MicroRNA的表达模式与抗病品种有所不同。一些MicroRNA在接种后表达水平持续上升，如miR482e、miR828a等，这些MicroRNA可能被病原菌利用，抑制黄瓜的抗病相关基因，从而促进病原菌的侵染和定殖；而另一些MicroRNA则呈现出先上升后下降的趋势，如miR167g、miR398b等，它们的表达变化可能与黄瓜在病原菌侵染初期的应激反应以及后期抗病能力的逐渐丧失有关。图2：差异表达MicroRNA的层次聚类分析热图（红色表示上调表达，绿色表示下调表达，颜色深浅表示表达差异的程度）3.3差异表达MicroRNA的验证为了确保高通量测序数据的可靠性，对筛选出的差异表达MicroRNA进行了验证，挑选了6个在抗病品种和感病品种间以及接种前后差异表达较为显著的MicroRNA，包括miR156a、miR164b、miR172c、miR396d、miR482e和miR828a，运用茎环法反转录和实时荧光定量PCR（qRT-PCR）技术进行验证。在qRT-PCR实验中，以U6作为内参基因，采用2⁻ΔΔCt法计算MicroRNA的相对表达量，每个样品设置3个技术重复，以确保数据的准确性和可靠性。实验结果显示，这6个MicroRNA的qRT-PCR验证结果与高通量测序数据的表达趋势基本一致（图3）。在抗病品种‘中农106号’中，miR156a和miR164b在接种枯萎病菌后表达水平显著上调，qRT-PCR结果显示，接种后48h，miR156a的表达量相较于接种前上调了[X]倍，miR164b的表达量上调了[X]倍，与测序数据中这两个MicroRNA的上调趋势相符；而miR172c和miR396d的表达水平则逐渐下调，接种后72h，miR172c的表达量相较于接种前下调了[X]倍，miR396d的表达量下调了[X]倍，与测序数据的下调趋势一致。在感病品种‘津优35号’中，miR482e和miR828a在接种后表达水平持续上升，qRT-PCR结果表明，接种后72h，miR482e的表达量相较于接种前上调了[X]倍，miR828a的表达量上调了[X]倍，与测序数据的表达趋势一致。这表明高通量测序筛选出的差异表达MicroRNA结果可靠，能够真实反映黄瓜在枯萎病菌侵染过程中MicroRNA的表达变化情况，为后续深入研究这些MicroRNA在黄瓜抗枯萎病中的作用奠定了坚实的基础。图3：差异表达MicroRNA的qRT-PCR验证结果（*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）四、黄瓜抗枯萎病MicroRNA靶基因的鉴定与验证4.1靶基因预测运用psRNATarget软件对筛选出的差异表达MicroRNA进行靶基因预测，设置最大期望值为3.0，允许的错配碱基数为4，其他参数采用默认值。经预测，共得到[X]个靶基因。这些靶基因分布在黄瓜基因组的不同区域，暗示着它们可能参与了多种生物学过程。为了深入了解这些靶基因的功能，将预测得到的靶基因与黄瓜基因注释文件进行比对，并参考GO数据库和KEGG数据库进行功能注释和富集分析。GO富集分析结果显示，靶基因主要富集在生物过程、分子功能和细胞组成等多个类别。在生物过程方面，主要参与了植物的防御反应、信号转导、转录调控、氧化还原过程以及植物激素信号转导等过程（图4）。在防御反应相关的GOterm中，如“responsetofungus”（对真菌的响应）、“defenseresponsetobacterium”（对细菌的防御反应）等，有多个靶基因富集，这表明这些靶基因可能在黄瓜抵御病原菌侵染的过程中发挥重要作用；在信号转导过程中，涉及到MAPK信号通路、钙信号通路等相关的靶基因也有显著富集，说明这些信号通路可能在黄瓜抗枯萎病过程中被激活，通过传递信号来调控植物的免疫反应。在分子功能方面，靶基因主要富集在核酸结合、转录因子活性、蛋白激酶活性、氧化还原酶活性以及水解酶活性等功能类别。具有核酸结合功能的靶基因，可能通过与DNA或RNA结合，参与基因的转录调控；转录因子活性相关的靶基因，能够调节下游基因的表达，在植物的生长发育和逆境响应中发挥关键作用；蛋白激酶活性相关的靶基因，可通过磷酸化作用调节蛋白的活性，进而影响植物的生理过程。KEGG富集分析结果表明，靶基因显著富集在植物-病原体互作、植物激素信号转导、MAPK信号通路-植物等多个重要的代谢通路（图5）。在植物-病原体互作通路中，多个靶基因参与了植物对病原菌的识别、信号传导以及防御反应的激活，它们可能编码一些与病原菌识别相关的受体蛋白、参与信号传导的激酶以及合成植保素等防御物质的关键酶，这些基因的表达变化可能直接影响黄瓜对枯萎病的抗性；在植物激素信号转导通路中，涉及生长素、水杨酸、茉莉酸等多种激素信号转导途径的靶基因富集，说明植物激素在黄瓜抗枯萎病过程中通过调控相关基因的表达，发挥着重要的信号传递和调节作用；MAPK信号通路-植物通路中的靶基因，在植物应对生物和非生物胁迫时，通过磷酸化级联反应传递信号，激活下游的防御基因，增强植物的抗性。图4：靶基因GO富集分析柱状图（横坐标为GOterm，纵坐标为富集的靶基因数量）图5：靶基因KEGG富集分析气泡图（横坐标为富集因子，纵坐标为KEGG通路名称，气泡大小表示富集的靶基因数量，颜色深浅表示P值大小）4.2靶基因验证4.2.15'-RACE验证为了进一步验证生物信息学预测的靶基因的准确性，采用5'-RACE技术对部分差异表达MicroRNA的靶基因进行验证。以miR156a的预测靶基因SPL9为例，设计特异性引物进行5'-RACE实验。实验结果显示，通过两轮PCR扩增，成功获得了预期大小的特异性条带（图6）。将该条带进行回收、克隆和测序，测序结果与预测的靶基因SPL9序列进行比对分析，发现miR156a对SPL9的切割位点位于mRNA序列的第[X]位核苷酸处（图7），与生物信息学预测的结合位点一致，表明SPL9确实是miR156a的靶基因，miR156a能够通过识别并结合SPL9mRNA的特定区域，对其进行切割，从而在转录后水平调控SPL9的表达。图6：miR156a靶基因SPL9的5'-RACE扩增结果（M：DNAMarker；1：第一轮PCR扩增产物；2：第二轮PCR扩增产物）图7：miR156a与SPL9mRNA的互补配对及切割位点示意图（红色箭头指示切割位点）对其他差异表达MicroRNA的靶基因进行5'-RACE验证，也得到了类似的结果。如miR164b的靶基因NAC1，经5'-RACE验证，确定了其切割位点位于mRNA序列的第[X]位核苷酸处；miR172c的靶基因AP2，切割位点位于mRNA序列的第[X]位核苷酸处。这些结果进一步证实了生物信息学预测的可靠性，明确了MicroRNA与靶基因之间的靶向关系，为后续深入研究MicroRNA对靶基因的调控机制奠定了坚实的基础。4.2.2双荧光素酶报告基因实验验证利用双荧光素酶报告基因实验，对MicroRNA与靶基因的互作关系进行了进一步验证。以miR156a和SPL9为例，构建了包含SPL9基因3'UTR野生型（WT）和突变型（MUT）序列的双荧光素酶报告基因载体。突变型载体是将SPL9基因3'UTR中与miR156a互补配对的关键位点进行突变，使其不能与miR156a结合。将构建好的报告基因载体分别与miR156a模拟物（mimic）或阴性对照（NC）共转染到黄瓜原生质体中，转染48h后，检测萤火虫荧光素酶和海肾荧光素酶的活性，计算萤火虫荧光素酶活性与海肾荧光素酶活性的比值（Luc/Ren）。实验结果表明，转染miR156amimic和野生型报告基因载体的组中，Luc/Ren值显著低于转染NC和野生型报告基因载体的组（图8），说明miR156a能够与SPL9基因3'UTR的野生型序列互补配对，抑制萤火虫荧光素酶的表达，从而验证了miR156a与SPL9之间的靶向互作关系；而转染miR156amimic和突变型报告基因载体的组中，Luc/Ren值与转染NC和突变型报告基因载体的组相比无显著差异，表明当SPL9基因3'UTR中与miR156a互补配对的位点发生突变后，miR156a无法与之结合，不能抑制萤火虫荧光素酶的表达，进一步证明了miR156a与SPL9之间的相互作用具有特异性。图8：miR156a与SPL9的双荧光素酶报告基因实验结果（**表示在P<0.01水平上差异极显著）对其他MicroRNA与靶基因的组合进行双荧光素酶报告基因实验验证，也得到了一致的结果。如miR164b与NAC1、miR172c与AP2等，均证实了MicroRNA能够与相应靶基因的3'UTR特异性结合，抑制靶基因的表达。这些结果进一步明确了MicroRNA在黄瓜抗枯萎病过程中对靶基因的调控作用，为深入解析黄瓜抗枯萎病的分子机制提供了有力的证据。五、黄瓜抗枯萎病MicroRNA及其靶基因的功能解析5.1MicroRNA及其靶基因的表达模式分析为了深入探究黄瓜抗枯萎病MicroRNA及其靶基因在抗病过程中的作用，对筛选出的关键MicroRNA及其靶基因在不同抗性黄瓜材料、不同组织和不同时间点的表达情况进行了系统分析。在不同抗性黄瓜材料中，以抗病品种‘中农106号’和感病品种‘津优35号’为研究对象，采用实时荧光定量PCR技术，检测了miR156a及其靶基因SPL9在枯萎病菌接种前后的表达水平变化。结果显示，在抗病品种‘中农106号’中，接种枯萎病菌后，miR156a的表达水平迅速上调，在接种后48h达到峰值，相较于接种前上调了[X]倍；而其靶基因SPL9的表达水平则显著下调，在接种后48h，SPL9的表达量相较于接种前下调了[X]倍。在感病品种‘津优35号’中，miR156a的表达水平在接种后虽有上升趋势，但幅度较小，在接种后72h仅上调了[X]倍；SPL9的表达水平在接种前后变化不明显，仅在接种后72h略有下降，下调幅度为[X]倍（图9）。这表明miR156a及其靶基因SPL9在不同抗性黄瓜材料中的表达模式存在显著差异，miR156a的高表达可能与黄瓜的抗病性密切相关，通过抑制靶基因SPL9的表达，激活黄瓜的抗病反应。图9：miR156a及其靶基因SPL9在不同抗性黄瓜材料中的表达模式分析（*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）对miR164b及其靶基因NAC1在不同抗性黄瓜材料中的表达分析也得到了类似结果。在抗病品种‘中农106号’中，接种枯萎病菌后，miR164b的表达水平显著上调，接种后48h上调了[X]倍；NAC1的表达水平则明显下调，接种后48h下调了[X]倍。在感病品种‘津优35号’中，miR164b的表达上调幅度较小，接种后72h上调了[X]倍；NAC1的表达水平在接种前后变化不显著。这些结果进一步证实了MicroRNA及其靶基因在不同抗性黄瓜材料中的表达差异与黄瓜的抗病性密切相关。在不同组织中，选取了黄瓜的根、茎、叶三个主要组织，检测了关键MicroRNA及其靶基因的表达情况。以miR172c及其靶基因AP2为例，在未接种枯萎病菌时，miR172c在黄瓜的叶中表达水平最高，茎中次之，根中最低；其靶基因AP2则在根中表达水平最高，茎中次之，叶中最低。接种枯萎病菌后，在抗病品种‘中农106号’的叶和茎中，miR172c的表达水平均显著上调，叶中在接种后48h上调了[X]倍，茎中在接种后48h上调了[X]倍；AP2的表达水平则显著下调，叶中在接种后48h下调了[X]倍，茎中在接种后48h下调了[X]倍。在根中，miR172c的表达水平略有上升，接种后72h上调了[X]倍；AP2的表达水平略有下降，接种后72h下调了[X]倍。在感病品种‘津优35号’中，叶和茎中miR172c的表达上调幅度较小，叶中在接种后72h上调了[X]倍，茎中在接种后72h上调了[X]倍；AP2的表达水平下降幅度也较小，叶中在接种后72h下调了[X]倍，茎中在接种后72h下调了[X]倍。根中miR172c和AP2的表达变化不明显（图10）。这说明MicroRNA及其靶基因在不同组织中的表达模式存在差异，且在抗病品种中，它们对枯萎病菌侵染的响应更为明显，暗示着不同组织在黄瓜抗枯萎病过程中可能通过不同的MicroRNA调控网络发挥作用。图10：miR172c及其靶基因AP2在黄瓜不同组织中的表达模式分析（*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）在不同时间点，对枯萎病菌接种后0h、12h、24h、48h和72h的黄瓜植株进行检测，分析关键MicroRNA及其靶基因的动态表达变化。以miR396d及其靶基因GRF4为例，在抗病品种‘中农106号’中，接种枯萎病菌后，miR396d的表达水平在12h开始上升，24h显著上调，相较于接种前上调了[X]倍，48h达到峰值，上调了[X]倍，随后略有下降，但在72h仍维持较高水平，上调了[X]倍；GRF4的表达水平则在12h开始下降，24h显著下调，相较于接种前下调了[X]倍，48h下调幅度最大，下调了[X]倍，72h虽有所回升，但仍低于接种前水平，下调了[X]倍。在感病品种‘津优35号’中，miR396d的表达水平上升缓慢，在72h仅上调了[X]倍；GRF4的表达水平在接种后下降幅度较小，72h下调了[X]倍（图11）。通过对不同时间点的表达模式分析，清晰地展现了MicroRNA及其靶基因在黄瓜抗枯萎病过程中的动态变化规律，为深入理解其作用机制提供了重要线索，表明它们在黄瓜应对枯萎病菌侵染的不同阶段发挥着不同的调控作用。图11：miR396d及其靶基因GRF4在不同时间点的表达模式分析（*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）5.2MicroRNA对靶基因的调控机制MicroRNA主要通过两种方式对靶基因进行调控，即切割靶基因mRNA和抑制其翻译过程，从而在转录后水平精细地调控基因的表达，进而影响黄瓜对枯萎病的抗性。在切割靶基因mRNA这一调控方式中，成熟的MicroRNA会与AGO蛋白等形成RNA诱导沉默复合体（RISC）。当RISC中的MicroRNA与靶基因mRNA的互补区域识别并结合后，AGO蛋白会发挥核酸内切酶的活性，对靶基因mRNA进行切割，使其降解，从而无法正常翻译为蛋白质，实现对靶基因表达的抑制。以miR156a及其靶基因SPL9为例，经5'-RACE验证和双荧光素酶报告基因实验证实，miR156a能够与SPL9mRNA的特定区域互补配对，RISC中的AGO蛋白在互补区域的特定位置对SPL9mRNA进行切割，导致SPL9mRNA降解，从而抑制SPL9基因的表达。在黄瓜抗枯萎病过程中，接种枯萎病菌后，抗病品种‘中农106号’中miR156a表达上调，大量的miR156a与SPL9mRNA结合并切割，使得SPL9蛋白表达量降低，从而激活下游抗病相关基因的表达，增强黄瓜的抗病性。MicroRNA还可以通过抑制靶基因mRNA的翻译过程来调控基因表达。在这种调控方式下，MicroRNA与靶基因mRNA结合后，并不对mRNA进行切割，而是阻碍核糖体与mRNA的结合，或者抑制核糖体在mRNA上的移动，从而抑制蛋白质的合成，降低靶基因的表达水平。虽然具体的分子机制尚未完全明确，但研究表明，这种翻译抑制作用可能与MicroRNA结合到mRNA的3'UTR区域，影响了翻译起始因子或其他相关蛋白与mRNA的相互作用有关。以miR164b及其靶基因NAC1为例，在抗病品种‘中农106号’接种枯萎病菌后，miR164b表达上调，通过抑制NAC1mRNA的翻译，减少NAC1蛋白的合成，从而调控下游基因的表达，参与黄瓜的抗枯萎病反应。MicroRNA对靶基因的调控作用在黄瓜抗枯萎病过程中具有重要影响。通过对靶基因的调控，MicroRNA能够参与激活植物的免疫反应，调节植物激素信号通路，进而增强黄瓜对枯萎病的抗性。在植物免疫反应方面，一些MicroRNA可以通过调控NBS-LRR类抗病基因等靶基因的表达，激活植物的免疫信号传导途径，促使植物产生植保素、病程相关蛋白等防御物质，增强植物对病原菌的防御能力。在植物激素信号通路调控方面，MicroRNA可以调节生长素、水杨酸、茉莉酸等植物激素信号转导相关基因的表达，影响植物激素的合成、运输和信号传导，从而间接调控黄瓜的抗病性。miR172c通过抑制靶基因AP2的表达，影响生长素信号通路，进而调控黄瓜对枯萎病的抗性。当黄瓜受到枯萎病菌侵染时，miR172c表达上调，抑制AP2的表达，使得生长素信号通路发生改变，激活下游抗病相关基因的表达，增强黄瓜的抗病能力。5.3靶基因的功能验证5.3.1基因过表达分析为深入探究靶基因在黄瓜抗枯萎病过程中的具体功能，构建了靶基因的过表达载体，并通过农杆菌介导的遗传转化方法，将其导入黄瓜子叶节，成功获得了过表达转基因植株。以靶基因SPL9为例，对获得的过表达转基因植株进行了PCR鉴定和测序验证，结果表明，外源SPL9基因已成功整合到黄瓜基因组中，且能够正常转录和表达。对过表达SPL9基因的转基因黄瓜植株进行苗期人工接种枯萎病菌，以野生型黄瓜植株作为对照。接种后，密切观察植株的发病情况，并按照病情分级标准进行病情调查，计算病情指数。实验结果显示，接种枯萎病菌后，野生型黄瓜植株在第4天开始出现发病症状，表现为基部叶片轻微萎蔫，随着时间的推移，病情逐渐加重，在第10天病情指数达到了[X]，植株生长明显受抑制，部分叶片萎蔫卷曲；而过表达SPL9基因的转基因黄瓜植株发病症状出现较晚，在第6天才开始出现轻微症状，且病情发展缓慢，在第10天病情指数仅为[X]，植株生长状况相对较好，大部分叶片仍保持正常的绿色和伸展状态（图12）。图12：过表达SPL9基因的转基因黄瓜植株与野生型植株接种枯萎病菌后的病情发展情况（A：野生型植株接种后4天；B：过表达SPL9基因的转基因植株接种后4天；C：野生型植株接种后10天；D：过表达SPL9基因的转基因植株接种后10天）进一步对过表达SPL9基因的转基因黄瓜植株和野生型植株在接种枯萎病菌前后的生理生化指标进行测定，结果表明，在接种前，两者的丙二醛（MDA）含量、超氧化物歧化酶（SOD）活性、过氧化物酶（POD）活性和过氧化氢酶（CAT）活性等生理生化指标无显著差异。接种枯萎病菌后，野生型植株的MDA含量迅速上升，在接种后第6天相较于接种前增加了[X]倍，表明细胞膜受到了严重的损伤；SOD、POD和CAT活性虽有所升高，但升高幅度较小，在接种后第6天，SOD活性相较于接种前升高了[X]%，POD活性升高了[X]%，CAT活性升高了[X]%。而过表达SPL9基因的转基因黄瓜植株在接种枯萎病菌后，MDA含量上升幅度较小，在接种后第6天相较于接种前仅增加了[X]倍，表明细胞膜损伤程度较轻；SOD、POD和CAT活性显著升高，在接种后第6天，SOD活性相较于接种前升高了[X]%，POD活性升高了[X]%，CAT活性升高了[X]%（图13）。这些结果表明，过表达SPL9基因能够提高黄瓜植株体内抗氧化酶的活性，降低细胞膜的损伤程度，从而增强黄瓜对枯萎病的抗性。图13：过表达SPL9基因的转基因黄瓜植株与野生型植株接种枯萎病菌后的生理生化指标变化（A：MDA含量；B：SOD活性；C：POD活性；D：CAT活性；*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）对其他靶基因进行过表达分析，也得到了类似的结果。过表达NAC1基因的转基因黄瓜植株在接种枯萎病菌后，病情指数显著低于野生型植株，且体内抗氧化酶活性升高，MDA含量降低，表明过表达NAC1基因能够增强黄瓜对枯萎病的抗性；过表达AP2基因的转基因黄瓜植株同样表现出较强的抗病性，病情发展缓慢，生理生化指标变化有利于提高植株的抗病能力。这些结果充分表明，靶基因的过表达能够显著增强黄瓜对枯萎病的抗性，为黄瓜抗枯萎病分子育种提供了重要的基因资源和理论依据。5.3.2基因沉默分析采用RNA干扰（RNAi）技术构建了靶基因的沉默载体，并将其转化黄瓜植株，成功获得了基因沉默植株。以靶基因SPL9为例，对获得的基因沉默植株进行了PCR鉴定和测序验证，结果表明，沉默载体已成功导入黄瓜植株，且能够有效抑制SPL9基因的表达。对沉默SPL9基因的黄瓜植株进行苗期人工接种枯萎病菌，以野生型黄瓜植株作为对照。接种后，定期观察植株的发病情况，并计算病情指数。实验结果显示，接种枯萎病菌后，沉默SPL9基因的黄瓜植株在第2天就开始出现发病症状，且病情发展迅速，在第6天病情指数达到了[X]，植株大部分叶片萎蔫卷曲，生长严重受抑制；而野生型黄瓜植株在第4天开始出现发病症状，在第6天病情指数为[X]，发病症状相对较轻（图14）。图14：沉默SPL9基因的黄瓜植株与野生型植株接种枯萎病菌后的病情发展情况（A：野生型植株接种后2天；B：沉默SPL9基因的黄瓜植株接种后2天；C：野生型植株接种后6天；D：沉默SPL9基因的黄瓜植株接种后6天）对沉默SPL9基因的黄瓜植株和野生型植株在接种枯萎病菌前后的生理生化指标进行测定，结果表明，在接种前，两者的MDA含量、SOD活性、POD活性和CAT活性等生理生化指标无显著差异。接种枯萎病菌后，沉默SPL9基因的黄瓜植株的MDA含量急剧上升，在接种后第4天相较于接种前增加了[X]倍，表明细胞膜受到了严重的损伤；SOD、POD和CAT活性虽有所升高，但升高幅度较小，在接种后第4天，SOD活性相较于接种前升高了[X]%，POD活性升高了[X]%，CAT活性升高了[X]%。而野生型黄瓜植株在接种枯萎病菌后，MDA含量上升幅度相对较小，在接种后第4天相较于接种前增加了[X]倍；SOD、POD和CAT活性升高幅度较大，在接种后第4天，SOD活性相较于接种前升高了[X]%，POD活性升高了[X]%，CAT活性升高了[X]%（图15）。这些结果表明，沉默SPL9基因会降低黄瓜植株体内抗氧化酶的活性，加剧细胞膜的损伤程度，从而削弱黄瓜对枯萎病的抗性。图15：沉默SPL9基因的黄瓜植株与野生型植株接种枯萎病菌后的生理生化指标变化（A：MDA含量；B：SOD活性；C：POD活性；D：CAT活性；*表示在P<0.05水平上差异显著，**表示在P<0.01水平上差异极显著）对其他靶基因进行基因沉默分析，也得到了一致的结果。沉默NAC1基因的黄瓜植株在接种枯萎病菌后，病情指数显著高于野生型植株，且体内抗氧化酶活性降低，MDA含量升高，表明沉默NAC1基因会降低黄瓜对枯萎病的抗性；沉默AP2基因的黄瓜植株同样表现出较弱的抗病性，病情发展迅速，生理生化指标变化不利于提高植株的抗病能力。这些结果充分证明，靶基因的沉默会显著降低黄瓜对枯萎病的抗性，进一步证实了靶基因在黄瓜抗枯萎病过程中的重要作用。六、讨论6.1黄瓜抗枯萎病MicroRNA及其靶基因的作用机制本研究通过高通量测序技术，系统地筛选并鉴定出了一系列与黄瓜抗枯萎病相关的MicroRNA及其靶基因，并对其作用机制进行了深入探究。在黄瓜受到枯萎病菌侵染后，抗病品种和感病品种中的MicroRNA表达谱发生了显著变化，这些差异表达的MicroRNA通过对靶基因的调控，参与了黄瓜的抗枯萎病过程。从调控方式来看，MicroRNA主要通过切割靶基因mRNA和抑制其翻译两种方式，在转录后水平对靶基因的表达进行调控。以miR156a及其靶基因SPL9为例，在抗病品种‘中农106号’接种枯萎病菌后，miR156a表达上调，与SPL9mRNA互补配对，RISC中的AGO蛋白对SPL9mRNA进行切割，导致其降解，从而抑制SPL9基因的表达。这种切割作用迅速且直接，能够在短时间内降低靶基因的mRNA水平，进而影响其下游基因的表达，激活黄瓜的抗病反应。miR164b对靶基因NAC1的调控则更多地表现为抑制翻译过程，虽然NAC1mRNA水平可能变化不明显，但蛋白质合成受到抑制，减少了NAC1蛋白的积累，从而调控下游基因的表达，参与黄瓜的抗枯萎病反应。这种转录后水平的精细调控，使得黄瓜能够根据病原菌的侵染情况，迅速调整基因表达，以应对病害胁迫。从功能角度分析，这些MicroRNA及其靶基因参与了植物的多个抗病相关过程。在植物免疫反应方面，它们通过调控NBS-LRR类抗病基因等靶基因的表达，激活植物的免疫信号传导途径。NBS-LRR类基因编码的蛋白在植物识别病原菌和激活免疫反应中起关键作用，MicroRNA对其表达的调控能够增强植物对病原菌的防御能力。一些MicroRNA还参与调节植物激素信号通路，如miR172c通过抑制靶基因AP2的表达，影响生长素信号通路，进而调控黄瓜对枯萎病的抗性。生长素在植物的生长发育和抗病过程中具有重要作用，其信号通路的改变能够激活下游抗病相关基因的表达，增强黄瓜的抗病能力。水杨酸、茉莉酸等植物激素信号转导途径也受到MicroRNA及其靶基因的调控，它们之间相互协调，共同参与黄瓜的抗枯萎病过程。与其他植物抗病机制相比，黄瓜抗枯萎病MicroRNA调控机制具有一定的相似性和独特性。在相似性方面，许多植物在应对病原菌侵染时，都通过MicroRNA对靶基因的调控来激活免疫反应和调节激素信号通路。在拟南芥中，miR482/2118家族成员通过靶向NBS-LRR类抗病基因，参与植物对细菌和真菌的免疫反应，与黄瓜中MicroRNA对NBS-LRR类基因的调控类似；在番茄中，miR164通过调控NAC1基因的表达，影响茉莉酸信号通路，进而调控植物的抗病性，这与黄瓜中miR164b对NAC1的调控以及对激素信号通路的影响具有相似之处。黄瓜抗枯萎病MicroRNA调控机制也具有其独特性。不同植物的MicroRNA及其靶基因的种类和表达模式存在差异，这与植物的进化历程和生态适应性密切相关。黄瓜中一些特有的MicroRNA及其靶基因，可能参与了黄瓜特有的抗病途径，或者在黄瓜应对枯萎病病原菌侵染时发挥着关键作用。黄瓜枯萎病是一种土传病害，病原菌通过根系侵入植物体内，黄瓜在长期进化过程中，可能形成了一套针对土传病原菌的MicroRNA调控机制，以增强对枯萎病的抗性。未来还需要进一步深入研究这些独特的调控机制，以全面揭示黄瓜抗枯萎病的分子基础。6.2本研究的创新点与不足之处本研究的创新点主要体现在研究思路和方法上。在研究思路方面，首次系统地从MicroRNA及其靶基因的角度，深入探究黄瓜抗枯萎病的分子机制，为黄瓜抗病研究提供了新的视角和方向。以往的研究大多集中在黄瓜枯萎病的病原菌鉴定、发病规律以及传统的抗病基因挖掘等方面，对MicroRNA在黄瓜抗枯萎病中的调控作用研究较少。本研究通过高通量测序技术，全面筛选和鉴定与黄瓜抗枯萎病相关的MicroRNA及其靶基因，揭示了它们在黄瓜抗枯萎病过程中的重要调控作用，填补了这一领域在MicroRNA研究方面的部分空白，有助于深入理解黄瓜与枯萎病菌互作的分子机制，为黄瓜抗枯萎病分子育种提供了新的理论基础。在研究方法上，本研究综合运用了高通量测序、生物信息学分析、分子生物学实验和转基因技术等多种先进技术手段，相互验证和补充，确保了研究结果的准确性和可靠性。通过高通量测序技术，能够快速、全面地获取黄瓜在枯萎病菌侵染前后MicroRNA的表达谱信息，为后续的差异表达分析和靶基因预测提供了丰富的数据基础；生物信息学分析则帮助我们对海量的测序数据进行有效的挖掘和分析，预测MicroRNA的靶基因，并对其功能进行注释和富集分析，为深入研究MicroRNA的调控机制提供了线索；分子生物学实验如5'-RACE、双荧光素酶报告基因实验等，从实验层面验证了MicroRNA与靶基因之间的靶向互作关系，明确了MicroRNA的调控方式；转基因技术则通过构建基因过表达和基因沉默载体，转化黄瓜植株，从表型和生理生化指标等方面深入研究了靶基因在黄瓜抗枯萎病过程中的功能，这种多技术联用的研究方法，提高了研究的深度和广度。本研究也存在一些不足之处。在样本数量方面，虽然每个处理设置了3个生物学重复，但考虑到黄瓜品种的多样性以及环境因素对实验结果的影响，样本数量可能相对较少，这可能会对实验结果的普遍性和代表性产生一定影响。在后续研究中，可以进一步扩大样本量，增加不同生态型的黄瓜品种，以及在不同环境条件下进行实验，以提高研究结果的可靠性和普适性，更全面地揭示黄瓜抗枯萎病MicroRNA及其靶基因的作用机制。在研究方法上，虽然运用了多种技术手段，但仍存在一定局限性。在MicroRNA的鉴定过程中，