基于生物信息学的玉米抗纹枯病相关miRNA预测与功能解析

基于生物信息学的玉米抗纹枯病相关miRNA预测与功能解析一、引言1.1研究背景玉米（ZeamaysL.）作为全球重要的粮食、饲料及工业原料作物，在农业生产与国民经济中占据关键地位。根据联合国粮食及农业组织（FAO）的数据，2023年全球玉米总产量达到12.1亿吨，种植面积约为1.97亿公顷，广泛分布于美洲、亚洲、欧洲和非洲等地区。在中国，玉米是三大主粮之一，2023年国内玉米产量达2.88亿吨，种植面积达4307万公顷，其产量与种植面积对保障国家粮食安全及促进农业经济发展意义重大。玉米不仅为人类提供碳水化合物、蛋白质、维生素和矿物质等营养成分，也是畜牧业重要的饲料来源，支撑着肉类、蛋类和奶制品的生产。同时，玉米在工业领域用途广泛，可加工成淀粉、糖浆、玉米油和生物燃料等多种产品，其在食品、造纸、纺织和能源等行业发挥着不可或缺的作用。玉米纹枯病（Cornsheathblight）是由立枯丝核菌（RhizoctoniasolaniKühn）引起的一种世界性真菌病害，对玉米生产构成严重威胁。该病菌寄主范围广泛，可侵染水稻、小麦、高粱等多种禾本科作物。在玉米上，纹枯病主要侵害叶鞘、果穗和茎秆。发病初期，叶鞘基部出现淡褐色水渍状小斑，随着病情发展，病斑逐渐扩大并融合成云纹状灰白色大斑，导致叶鞘腐败、叶片枯死，并向上蔓延。湿度大时，发病部位可见茂盛的菌丝体，后期形成白色小绒球，逐渐变成褐色的菌核。严重时，病株茎秆松软，组织解体，果穗干缩、腐败，造成玉米严重减产甚至绝收。近年来，随着玉米种植面积的扩大、密植程度的增加以及多年连作等因素的影响，玉米纹枯病的发生呈加重趋势。据统计，在中国主要玉米产区，纹枯病的发病率已超过70%，减产幅度达10%-20%，严重时可高达35%。在全球范围内，每年因玉米纹枯病造成的经济损失高达数十亿美元。例如，在东南亚地区，由于高温高湿的气候条件适宜纹枯病的发生，该地区的玉米产量受到严重影响。此外，随着气候变化，极端天气事件增多，纹枯病的发生范围和危害程度可能进一步扩大。纹枯病的发生与多种因素密切相关。病原菌以菌核在土壤中越冬，成为翌年初次侵染源。在适宜的温湿度条件下，菌核萌发长出菌丝，侵染玉米基部叶鞘，并逐渐向上蔓延。温度在25-30℃、相对湿度90%以上是玉米纹枯病发生发展的适宜气候条件。在病害发生期，特别是梅雨季节，雨期长、雨日多，湿度大，病害发生重。此外，土壤肥沃程度、施肥合理性、种植密度、地势和排水条件等也会影响纹枯病的发生。土壤肥沃发病较轻，反之则发病较重；偏施氮肥、种植密度过大、植株间通风透光条件差、地势低洼、排水不良的田块发病重；连作田比轮作田发病重。防治玉米纹枯病对于保障玉米产业的可持续发展具有重要意义。目前，针对玉米纹枯病的防治主要采取农业防治、化学防治和选育抗病品种等措施。农业防治包括加强清洁田园，将玉米收获后的秸秆远离农田堆放，田间残茬进行烧毁或深埋等无害化处理；控制种植密度，根据玉米品种、土壤和光照条件确定合理密度；加强田间管理，及时剥除下部病叶鞘及叶片；科学施肥，控制氮肥施入量，增施磷肥、钾肥和锌肥。化学防治主要选用5%的井冈霉素水剂1500倍液、40%核净可粉剂1000溶液或50%甲基硫菌灵可湿性粉剂500倍液等进行喷雾防治。然而，农业防治措施实施较为繁琐，且效果受环境因素影响较大；化学防治虽能在短期内有效控制病害，但长期使用易导致病原菌产生抗药性，同时造成环境污染和食品安全问题；而选育抗病品种是最经济、有效的防治手段，但目前优质高抗玉米品种相对匮乏，且传统育种周期长、效率低。因此，深入研究玉米抗纹枯病的分子机制，挖掘新的抗纹枯病基因资源，对于培育抗病品种、实现玉米纹枯病的绿色防控具有重要的理论和实践意义。1.2miRNA的研究概述微小核糖核酸（microRNA，miRNA）是一类内生的、长度约为21-24个核苷酸的非编码单链RNA分子，在真核生物中广泛存在。1993年，美国科学家VictorAmbros和GaryRuvkun在秀丽隐杆线虫（Caenorhabditiselegans）中首次发现了lin-4基因，它并不编码蛋白质，而是产生一段长度约22个核苷酸的非编码RNA，通过与lin-14基因的mRNA的3'非翻译区互补配对，抑制lin-14mRNA的翻译过程，从而调控线虫的发育时序，这便是最早被发现的miRNA。然而，这一发现起初并未引起科学界的广泛关注。直到2001年，3个研究小组几乎同时在线虫、果蝇和人类中鉴定出大量类似的小分子RNA，才掀起了miRNA的研究热潮。此后，随着高通量测序技术和生物信息学的飞速发展，越来越多的miRNA被发现和鉴定，其在生物体内的重要调控功能也逐渐被揭示。miRNA具有一些独特的特点。首先，它长度较短，一般为21-24个核苷酸，这种短小的结构使其能够高效地与靶mRNA相互作用。其次，miRNA在生物进化过程中具有高度的保守性，一些miRNA家族在不同物种中具有相似的序列和功能，这表明它们在生物的生存和繁衍中发挥着至关重要的作用。此外，miRNA的表达具有时空特异性，在不同的组织、器官以及发育阶段，miRNA的表达谱存在显著差异，这也决定了其调控功能的特异性。miRNA的来源主要是基因组DNA，其生物合成过程较为复杂。在细胞核中，由RNA聚合酶Ⅱ转录生成初始miRNA（primarymiRNA，pri-miRNA），pri-miRNA通常具有较长的核苷酸序列，且包含多个茎环结构。随后，pri-miRNA在Drosha酶和DGCR8蛋白组成的复合体作用下，被切割成约70-100个核苷酸长度的前体miRNA（precursormiRNA，pre-miRNA），pre-miRNA呈发夹状结构。接着，pre-miRNA在Exportin-5转运蛋白的协助下，从细胞核转运到细胞质中。在细胞质中，pre-miRNA被Dicer酶进一步切割，形成成熟的miRNA双链。最后，miRNA双链中的一条链被降解，另一条链则与AGO（Argonaute）蛋白等结合形成RNA诱导沉默复合体（RNA-inducedsilencingcomplex，RISC），发挥其生物学功能。在植物生长发育过程中，miRNA扮演着不可或缺的角色。例如，在植物的胚胎发育过程中，miRNA参与调控胚胎的极性建立、细胞分化和器官形成。研究发现，miR160和miR167通过调控生长素响应因子（ARFs）的表达，影响拟南芥胚胎的发育。在植物的营养生长阶段，miRNA参与调控叶片的形态建成、茎的伸长和根系的发育。miR164通过靶向NAC1基因，调控侧根的发育；miR396通过调控生长调节因子（GRFs）的表达，影响叶片的大小和形状。在植物的生殖生长阶段，miRNA参与调控花的发育、花粉的形成和果实的成熟。miR172通过调控AP2-like转录因子的表达，影响花器官的发育和开花时间。miRNA在植物抗逆过程中也发挥着关键作用，其作用机制主要通过调控相关基因的表达，使植物能够更好地适应逆境环境。在植物应对生物胁迫时，例如受到病原菌侵染，miRNA可参与植物的免疫反应。研究表明，在水稻受到稻瘟病菌侵染时，miR168通过靶向AGO1基因，调控水稻对稻瘟病的抗性；miR393通过靶向生长素受体基因TIR1和AFB2，激活植物的防御反应，增强水稻对稻瘟病和白叶枯病的抗性。在玉米抗纹枯病过程中，也可能存在类似的miRNA调控机制，通过识别并结合纹枯病菌侵染相关的靶mRNA，影响相关基因的表达，从而参与玉米对纹枯病的防御反应。在非生物胁迫方面，如干旱、盐渍、低温等逆境条件下，miRNA同样发挥重要作用。在干旱胁迫下，植物体内的miR169通过靶向NF-YA转录因子基因，调控植物的抗旱性；miR398通过靶向铜锌超氧化物歧化酶（Cu/Zn-SOD）基因，参与植物对氧化胁迫的响应。在盐胁迫下，miR171通过调控SCL转录因子基因的表达，影响植物的耐盐性。1.3生物信息学在miRNA研究中的应用随着生命科学研究进入大数据时代，生物信息学作为一门由生命科学与信息科学等多学科交叉融合的新兴学科，在miRNA研究领域发挥着不可或缺的作用。植物miRNA研究从小规模实验逐步向大规模计算分析方向迈进，生物信息学的应用已成为该领域的研究热点。在植物miRNA的研究进程中，生物信息学具有至关重要的意义。传统的实验方法，如克隆测序、芯片技术和深度测序等，虽能直接检测miRNA，但存在成本高、效率低、通量有限等弊端。例如，克隆测序需构建文库，操作繁琐，且易遗漏低丰度miRNA；芯片技术虽能实现高通量检测，但受限于探针设计，对新miRNA的发现能力不足；深度测序虽能获取大量miRNA信息，但后续数据分析复杂。而生物信息学方法能有效弥补这些不足，它可利用计算机强大的运算能力和数据分析算法，对海量的基因组数据、转录组数据进行高效处理和分析，从庞大的生物数据中快速筛选出潜在的miRNA及其靶基因，为实验研究提供精准的候选目标，极大地提高了研究效率，降低了研究成本，推动植物miRNA研究从局部、单一的研究模式向系统、全面的方向发展。在miRNA预测方面，常用的生物信息学方法主要基于miRNA的结构特征、序列保守性以及与已知miRNA的同源性等。从结构特征角度，miRNA通常由具有发夹结构的前体转录本加工而来，其前体具有特定的二级结构，如茎环结构，且茎部区域碱基互补配对程度较高，形成稳定的双链结构。基于此，可通过RNA二级结构预测算法，如Mfold、RNAfold等软件，对基因组序列进行扫描，预测可能形成miRNA前体发夹结构的区域。这些软件利用热力学原理，计算RNA序列形成不同二级结构的自由能，从而筛选出具有最低自由能、最稳定的发夹结构，作为潜在的miRNA前体。例如，在拟南芥miRNA研究中，通过Mfold软件对其基因组进行分析，成功预测出多个新的miRNA前体，后续实验验证了这些预测结果的可靠性。从序列保守性角度，由于miRNA在物种进化过程中具有一定的保守性，一些miRNA家族在不同物种中具有相似的序列，可通过将待预测物种的基因组序列与已知物种的miRNA序列进行比对，寻找保守的序列区域，进而预测潜在的miRNA。常用的比对工具如BLAST（BasicLocalAlignmentSearchTool），能快速在不同物种的序列数据库中搜索相似序列。以水稻和玉米为例，通过BLAST比对，发现了一些在两者中保守存在的miRNA序列，为进一步研究这些miRNA在禾本科植物中的功能提供了线索。从同源性角度，若一个新的miRNA与已知的miRNA在序列和结构上具有较高的同源性，那么它很可能也是一个功能性miRNA。通过构建miRNA序列数据库，将新发现的RNA序列与数据库中的已知miRNA进行同源性分析，可判断其是否为新的miRNA。例如，在番茄miRNA研究中，通过同源性分析，从新测序数据中鉴定出多个与已知植物miRNA具有高度同源性的新miRNA。预测miRNA靶基因同样依赖多种生物信息学方法。基于序列互补配对原理，miRNA通过与靶基因mRNA的特定区域互补配对，实现对靶基因表达的调控，因此可利用算法寻找与miRNA序列互补配对的mRNA区域，确定潜在的靶基因。常用的靶基因预测软件有TargetFinder、psRNATarget等。其中，TargetFinder主要针对动物miRNA靶基因预测，它考虑了miRNA与mRNA互补配对的自由能、种子序列的匹配程度等因素；psRNATarget则专门用于植物miRNA靶基因预测，它结合了植物miRNA与靶基因互补配对的特点，如植物miRNA与靶基因的互补配对更为严格，通常要求几乎完全互补，能更准确地预测植物miRNA的靶基因。在实际应用中，以玉米抗纹枯病相关miRNA研究为例，利用psRNATarget软件对预测出的miRNA进行靶基因预测，发现了多个与植物免疫反应、信号传导相关的潜在靶基因，为后续研究miRNA在玉米抗纹枯病中的作用机制提供了重要线索。此外，还可结合基因表达数据、蛋白质-蛋白质相互作用数据等多组学数据，综合分析确定miRNA的靶基因。通过整合不同层面的数据信息，能更全面、准确地揭示miRNA与靶基因之间的调控关系，避免单一方法预测带来的假阳性或假阴性结果。1.4研究目的与意义本研究旨在借助生物信息学手段，深入预测玉米抗纹枯病相关的miRNA，并全面分析其功能，具体目标如下：利用生物信息学方法，从玉米基因组数据中精准预测与抗纹枯病相关的miRNA，构建玉米抗纹枯病相关miRNA的数据集；运用生物信息学工具和算法，预测上述miRNA的靶基因，探究miRNA与靶基因之间的调控关系；从生物学功能、信号通路等方面，深入分析预测得到的miRNA及其靶基因在玉米抗纹枯病过程中的潜在作用机制；通过实验验证部分预测结果，为后续深入研究玉米抗纹枯病的分子机制提供坚实的理论基础和可靠的实验依据。本研究在理论和实践方面均具有重要意义。在理论层面，深入揭示玉米抗纹枯病相关miRNA及其靶基因的作用机制，将极大地丰富植物抗病分子生物学的理论体系，有助于深入理解植物与病原菌互作的分子基础，为进一步探究植物抗病的复杂调控网络提供新的思路和方向。在实践应用方面，研究成果可为玉米抗病育种提供关键的基因资源和分子标记，加速抗病品种的选育进程，提高玉米的抗病能力，减少纹枯病对玉米产量和品质的影响，从而保障粮食安全；此外，为玉米纹枯病的绿色防控提供新的理论依据和技术支持，有助于减少化学农药的使用，降低环境污染，促进农业的可持续发展。二、材料与方法2.1试验材料玉米供试材料选用对纹枯病具有不同抗性的玉米自交系，包括高抗纹枯病的自交系R15，该自交系经多年田间接种鉴定及抗性评价，在多种环境条件下对纹枯病表现出稳定的高抗性；感病自交系黄早4，其对纹枯病较为敏感，在常规种植条件下易感染纹枯病且发病程度较重。将这些玉米自交系种植于试验田中，按照常规的玉米种植管理方法进行栽培，包括适时播种、合理密植、科学施肥与灌溉等，以确保玉米植株的正常生长发育，为后续实验提供充足且生长状态良好的材料。纹枯病接菌菌种为立枯丝核菌（RhizoctoniasolaniKühn）AG-1-IA融合群，该融合群是我国玉米纹枯病的优势致病融合群，具有较强的致病性。菌种从玉米纹枯病典型病株上分离获得，经过形态学观察及分子生物学鉴定确定其种类。将分离得到的菌种在马铃薯葡萄糖琼脂（PDA）培养基上进行培养，PDA培养基配方为：马铃薯200g、葡萄糖20g、琼脂15-20g、水1000mL，自然pH值。培养条件为28℃恒温培养箱中培养5-7天，待菌落长满平板且菌丝生长旺盛时，用于后续的接种实验。克隆载体选用pMD18-TVector，购自TaKaRa公司。该载体具有氨苄青霉素抗性基因，便于后续的阳性克隆筛选；其多克隆位点丰富，可方便地插入目的基因片段。宿主菌为大肠杆菌DH5α，购自TransGenBiotech公司。大肠杆菌DH5α具有易于转化、生长迅速等特点，适合作为克隆载体的宿主菌进行基因克隆和扩增实验。在实验前，将大肠杆菌DH5α接种于含有LB培养基的试管中，LB培养基配方为：胰蛋白胨10g、酵母提取物5g、氯化钠10g、琼脂15-20g（固体培养基添加）、水1000mL，pH7.0-7.2。于37℃、200rpm摇床振荡培养过夜，使其处于对数生长期，用于感受态细胞的制备及后续的转化实验。在生物信息学分析过程中，使用的数据库包括miRBase数据库（/），该数据库是国际上权威的miRNA序列数据库，收录了大量已知的miRNA序列及其相关信息，涵盖了多种物种，为miRNA的序列比对、同源性分析等提供了重要的数据支持；NCBI（NationalCenterforBiotechnologyInformation）数据库（/），包含丰富的核酸、蛋白质等生物分子序列数据，以及基因表达、生物医学文献等多种类型的数据资源，在基因注释、功能分析等方面发挥着关键作用；玉米基因组数据库（MaizeGDB，/），专门针对玉米基因组数据进行整合和管理，提供玉米基因的结构、功能、表达等详细信息，是玉米基因研究的重要数据平台。使用的分析软件有Mfold软件，用于预测RNA的二级结构，通过计算RNA序列形成不同二级结构的自由能，寻找最稳定的发夹结构，以此预测潜在的miRNA前体；RNAfold软件，同样基于热力学原理预测RNA二级结构，在预测miRNA前体的发夹结构方面具有较高的准确性和可靠性；BLAST软件，进行序列比对分析，可快速在不同物种的序列数据库中搜索相似序列，用于判断预测的miRNA与已知miRNA的同源性，以及预测miRNA的靶基因；TargetFinder软件，用于动物miRNA靶基因预测，考虑了miRNA与mRNA互补配对的自由能、种子序列的匹配程度等因素，在分析miRNA与靶基因的相互作用关系中发挥重要作用；psRNATarget软件，专门用于植物miRNA靶基因预测，结合了植物miRNA与靶基因互补配对的特点，能更准确地预测植物miRNA的靶基因，为研究植物miRNA的调控机制提供有力工具。主要试剂包括RNA提取试剂TRIzol，购自Invitrogen公司，用于从玉米组织中提取总RNA，其能够有效地裂解细胞，保持RNA的完整性，并通过分层、沉淀等步骤实现RNA的纯化；反转录试剂盒，采用TaKaRa公司的PrimeScriptRTreagentKitwithgDNAEraser，可将提取的总RNA反转录为cDNA，用于后续的PCR扩增及基因表达分析；PCR扩增试剂，选用TransGenBiotech公司的TransTaqHiFiPCRSuperMix，具有高保真度、扩增效率高等优点，能准确地扩增目的基因片段；DNAMarker，如DL2000DNAMarker，购自TaKaRa公司，用于在琼脂糖凝胶电泳中确定PCR产物的大小；限制性内切酶，如EcoRI、HindIII等，购自NewEnglandBiolabs公司，用于切割DNA片段，构建重组质粒；T4DNA连接酶，购自TaKaRa公司，可将切割后的DNA片段与载体进行连接，形成重组DNA分子。主要试剂盒有质粒提取试剂盒，采用TIANGEN公司的PlasmidMiniKit，用于从大肠杆菌中提取重组质粒，操作简便、快捷，能高效地获得高质量的质粒DNA；胶回收试剂盒，选用OMEGA公司的GelExtractionKit，可从琼脂糖凝胶中回收目的DNA片段，去除杂质和引物等，保证回收DNA的纯度和完整性。酶类及相关试剂包括RNase抑制剂，购自TaKaRa公司，可抑制RNA酶的活性，防止RNA在提取和后续操作过程中被降解；dNTPs，由四种脱氧核糖核苷酸组成，购自TaKaRa公司，为PCR扩增提供原料；MgCl₂，用于调节PCR反应体系的离子强度，影响TaqDNA聚合酶的活性，购自Sigma公司。实验主要仪器设备有高速冷冻离心机，型号为Eppendorf5424R，用于细胞裂解液的离心分离、RNA和DNA的沉淀等操作，可在低温条件下高速离心，保证生物分子的活性和完整性；PCR仪，型号为Bio-RadT100，用于进行PCR扩增反应，可精确控制反应温度和时间，实现目的基因的大量扩增；凝胶成像系统，型号为Bio-RadGelDocXR+，用于观察和分析琼脂糖凝胶电泳后的DNA条带，通过成像和分析软件，可对条带的亮度、大小等进行定量分析；核酸蛋白测定仪，型号为NanoDrop2000，用于测定RNA、DNA和蛋白质的浓度和纯度，通过紫外分光光度法，快速准确地获取样品的相关信息；恒温培养箱，型号为ThermoScientificHeratherm，用于培养大肠杆菌、玉米纹枯病菌等微生物，提供稳定的温度环境；超净工作台，型号为苏净安泰SW-CJ-2FD，为实验操作提供无菌环境，防止微生物污染；移液器，包括EppendorfResearchplus系列移液器，用于准确移取各种试剂和样品，量程涵盖0.1-1000μL，满足不同实验需求。2.2实验方法2.2.1miRNA的生物信息学分析从miRBase数据库中下载已知的植物miRNA序列，将下载的序列与玉米基因组数据库中的序列进行比对，通过BLAST软件进行序列相似性搜索，设定E值阈值为1e-5，筛选出与玉米基因组匹配且相似度较高的序列作为候选玉米miRNA。运用Mfold和RNAfold软件对候选miRNA序列进行二级结构预测，根据miRNA前体具有典型发夹结构的特征，筛选出具有稳定发夹结构、最小自由能较低的序列作为最终预测的玉米miRNA。利用psRNATarget软件对预测得到的玉米miRNA进行靶基因预测。将预测的miRNA序列输入到psRNATarget软件中，选择玉米基因组数据库作为参考基因组，设定错配数不超过4个，允许的最大凸起数为2个，严格匹配区错配、缺失的惩罚分数为1，G:U惩罚为0.5，非严格匹配区的惩罚分数减半等参数，运行软件进行靶基因预测。预测完成后，对得到的靶基因进行功能注释，通过与NCBI数据库中的基因功能信息进行比对，分析靶基因在生物过程、分子功能和细胞组成等方面的功能，筛选出与植物抗病、信号传导、代谢途径等相关的靶基因进行深入研究。2.2.2半定量检测预测得到的玉米microRNA前体选取生长至三叶期的玉米高抗纹枯病自交系R15和感病自交系黄早4植株，每个自交系选取3株生长状况一致的植株，分别在每株植株的第二片完全展开叶上进行纹枯病菌接种处理。将培养好的立枯丝核菌菌丝块（直径约5mm）用无菌镊子小心放置在叶片中部，并用湿润的脱脂棉覆盖，以保持接种部位的湿度，促进病菌侵染。同时设置未接种的健康植株作为对照，每个处理重复3次。接种后，将植株置于温度为28℃、相对湿度90%以上的人工气候箱中培养，分别在接种后0h、12h、24h、48h和72h采集叶片样品，迅速放入液氮中速冻，然后转移至-80℃冰箱保存备用。采用TRIzol法提取各处理和对照样品的总RNA。取约100mg冷冻的叶片组织，在液氮中充分研磨成粉末状，加入1mLTRIzol试剂，剧烈振荡混匀，室温静置5min，使组织充分裂解。然后加入0.2mL氯仿，剧烈振荡15s，室温静置2min，使溶液分层。4℃、12000g离心15min，将上层水相转移至新的离心管中，加入等体积的异丙醇，轻轻混匀，室温静置10min，使RNA沉淀。4℃、12000g离心10min，弃上清，RNA沉淀用1mL75%乙醇洗涤2次，4℃、7500g离心5min，弃上清，将RNA沉淀在室温下晾干，加入适量的DEPC处理水溶解RNA。用核酸蛋白测定仪测定RNA的浓度和纯度，要求OD260/OD280比值在1.8-2.0之间，OD260/OD230比值大于2.0，以确保RNA的质量。使用TaKaRa公司的PrimeScriptRTreagentKitwithgDNAEraser进行RNA反转录为cDNA。在冰上配制20μL反转录反应体系，包括5×PrimeScriptBuffer4μL、PrimeScriptRTEnzymeMixI1μL、OligodTPrimer（50μM）1μL、Random6mers（100μM）1μL、总RNA1μg，用DEPC处理水补足至20μL。将反应体系轻轻混匀，短暂离心后，置于PCR仪中进行反转录反应，反应条件为：37℃15min，85℃5s，4℃保存。反应结束后，得到的cDNA可直接用于后续实验或保存于-20℃冰箱备用。以玉米看家基因18SrRNA作为内参基因，设计其引物序列。18SrRNA上游引物：5'-CGGCTACCACATCCAAGGAA-3'，下游引物：5'-GCTGGAATTACCGCGGCT-3'。同时设计预测得到的玉米microRNA前体的特异性引物，上游引物根据前体序列的5'端设计，下游引物根据前体序列的3'端设计，引物长度一般为18-25bp，Tm值在55-65℃之间。以反转录得到的cDNA为模板，分别进行看家基因和microRNA前体的PCR扩增。PCR反应体系为25μL，包括2×TransTaqHiFiPCRSuperMix12.5μL、上游引物（10μM）1μL、下游引物（10μM）1μL、cDNA模板1μL，用ddH2O补足至25μL。PCR反应条件为：95℃预变性5min；95℃变性30s，58℃退火30s，72℃延伸30s，共进行28-32个循环（根据目的基因的表达量调整循环数，确保扩增处于线性阶段）；72℃终延伸10min。扩增结束后，取5μLPCR产物进行1.5%琼脂糖凝胶电泳检测，电泳缓冲液为1×TAE，电压为120V，电泳时间约30min。在凝胶成像系统下观察电泳条带，根据看家基因条带的亮度调整各样品cDNA的浓度，使各样本看家基因的表达量一致，以消除模板量差异对后续结果的影响。根据预测得到的玉米microRNA前体序列，设计特异性引物进行PCR扩增。PCR反应体系和条件同上述调整cDNA浓度时的PCR反应。扩增结束后，取10μLPCR产物进行2%琼脂糖凝胶电泳检测，电泳条件同上。在凝胶成像系统下观察电泳条带，用干净的手术刀将目的条带从凝胶上切下，放入1.5mL离心管中，使用OMEGA公司的GelExtractionKit进行胶回收。按照试剂盒说明书的步骤，加入适量的BindingBuffer，65℃水浴加热10min，使凝胶完全溶解。将溶解后的溶液转移至吸附柱中，12000g离心1min，弃滤液。用WashBuffer洗涤吸附柱2次，每次12000g离心1min，弃滤液。将吸附柱放入新的离心管中，加入适量的ElutionBuffer，室温静置2min，12000g离心1min，收集离心管中的溶液，即为回收的目的片段。将回收的目的片段连接到pMD18-TVector上，构建重组质粒。连接反应体系为10μL，包括pMD18-TVector1μL、回收的目的片段4μL、SolutionI5μL，轻轻混匀，16℃连接过夜。将连接产物转化到大肠杆菌DH5α感受态细胞中。从-80℃冰箱中取出DH5α感受态细胞，置于冰上解冻，加入10μL连接产物，轻轻混匀，冰浴30min。然后将离心管放入42℃水浴中热激90s，迅速放回冰浴中冷却2min。加入800μL不含抗生素的LB液体培养基，37℃、200rpm振荡培养1h，使菌体复苏。将培养后的菌液均匀涂布在含有氨苄青霉素（Amp）的LB固体培养基平板上，37℃倒置培养过夜。次日，从平板上挑取白色单菌落，接种到含有Amp的LB液体培养基中，37℃、200rpm振荡培养过夜。使用TIANGEN公司的PlasmidMiniKit提取重组质粒，按照试剂盒说明书的步骤进行操作。提取得到的重组质粒送测序公司进行测序，将测序结果与预测的玉米microRNA前体序列进行比对，验证序列的准确性。2.2.3RNA加尾和引物延伸实时荧光定量PCR法实时定量检测microRNA采用mirVanamiRNAIsolationKit从上述接种纹枯病菌不同时间的玉米叶片样品和对照样品中分离小分子RNA。取约100mg冷冻的叶片组织，在液氮中研磨成粉末，加入适量的裂解液，剧烈振荡混匀，室温静置5min。然后加入氯仿，振荡混匀，室温静置2min，4℃、12000g离心15min，将上层水相转移至新的离心管中。加入等体积的无水乙醇，轻轻混匀，将混合液转移至过滤柱中，12000g离心1min，弃滤液。用WashSolution洗涤过滤柱2次，每次12000g离心1min，弃滤液。将过滤柱放入新的离心管中，加入适量的ElutionSolution，室温静置2min，12000g离心1min，收集离心管中的溶液，即为分离得到的小分子RNA。用核酸蛋白测定仪测定小分子RNA的浓度和纯度，确保其质量符合后续实验要求。使用TaKaRa公司的miRcute增强型miRNA荧光定量检测试剂盒对分离得到的小分子RNA进行实时荧光定量PCR验证及表达分析。首先对miRNA进行加尾反应，在冰上配制10μL加尾反应体系，包括5×miRNAPoly(A)PolymeraseBuffer2μL、ATP(10mM)1μL、miR-PAP1μL、小分子RNA5μL，用DEPC处理水补足至10μL。将反应体系轻轻混匀，短暂离心后，置于37℃恒温金属浴中反应60min。然后进行引物延伸反应，在冰上配制20μL引物延伸反应体系，包括5×miRNARTBuffer4μL、dNTPs(10mM)2μL、miRNARTPrimer1μL、M-MLVReverseTranscriptase1μL、加尾后的miRNA5μL，用DEPC处理水补足至20μL。将反应体系轻轻混匀，短暂离心后，置于PCR仪中进行反应，反应条件为：42℃60min，70℃15min，4℃保存。反应结束后，得到的cDNA可直接用于实时荧光定量PCR或保存于-20℃冰箱备用。以U6snRNA作为内参基因，设计其引物序列。U6上游引物：5'-CTCGCTTCGGCAGCACA-3'，下游引物：5'-AACGCTTCACGAATTTGCGT-3'。同时设计预测得到的玉米miRNA的特异性引物，引物设计遵循茎环引物设计原则，根据miRNA的成熟序列设计茎环引物的环部序列，使其与miRNA的3'端互补配对，茎部序列为通用序列。以引物延伸反应得到的cDNA为模板，分别进行U6和miRNA的实时荧光定量PCR扩增。实时荧光定量PCR反应体系为20μL，包括2×miRcutemiRNAPremix10μL、上游引物（10μM）0.4μL、下游引物（10μM）0.4μL、cDNA模板2μL，用ddH2O补足至20μL。反应在实时荧光定量PCR仪上进行，反应条件为：95℃预变性15min；95℃变性10s，60℃退火延伸30s，共进行40个循环。在每个循环的退火延伸阶段采集荧光信号，反应结束后，根据仪器自带的软件分析Ct值。采用2-ΔΔCt法计算miRNA的相对表达量，比较不同处理和对照样品中miRNA的表达差异，分析miRNA在玉米抗纹枯病过程中的表达模式。三、结果与分析3.1microRNA及microRNA靶基因的预测结果通过生物信息学分析，从玉米基因组数据中成功预测得到了一系列潜在的microRNA。经序列比对和二级结构分析，共筛选出56个具有典型miRNA特征的序列，这些序列的前体能够形成稳定的发夹结构，且最小自由能较低，符合miRNA前体的结构特点。部分预测得到的microRNA及其前体的序列信息如表1所示。microRNA名称前体序列成熟序列发夹结构最小自由能（kcal/mol）miR156agugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugcugc3.2预测得到的microRNA前体半定量结果对各处理和对照样品进行总RNA提取后，通过核酸蛋白测定仪检测RNA的浓度和纯度。结果显示，所有样品的OD260/OD280比值均在1.8-2.0之间，OD260/OD230比值均大于2.0，表明提取的总RNA纯度较高，无蛋白质和多糖等杂质污染，完整性良好，可满足后续反转录及PCR实验的要求。将提取的总RNA反转录为cDNA后，以玉米看家基因18SrRNA为内参，进行PCR扩增。扩增产物经1.5%琼脂糖凝胶电泳检测，结果如图1所示。从图中可以看出，各样本中18SrRNA条带亮度基本一致，表明通过以18SrRNA为内参调整各样品cDNA浓度的方法有效，消除了模板量差异对后续结果的影响，确保了后续实验结果的准确性和可靠性。[此处插入图1：各样本看家基因18SrRNAPCR扩增电泳图，图中M为DNAMarker，1-6分别为不同处理和对照样品的cDNA模板扩增产物][此处插入图1：各样本看家基因18SrRNAPCR扩增电泳图，图中M为DNAMarker，1-6分别为不同处理和对照样品的cDNA模板扩增产物]利用设计的特异性引物对预测得到的玉米microRNA前体进行PCR扩增，扩增产物经2%琼脂糖凝胶电泳检测。结果显示，在预期大小的位置出现了清晰的条带，与预测的microRNA前体大小相符，表明成功扩增出了目的片段。对扩增得到的目的片段进行回收，回收产物经1%琼脂糖凝胶电泳检测，条带单一、明亮，无明显杂带，说明回收的目的片段纯度较高，可用于后续的克隆实验。将回收的目的片段连接到pMD18-TVector上，转化大肠杆菌DH5α，挑取白色单菌落进行培养并提取重组质粒。重组质粒经EcoRI和HindIII双酶切鉴定，酶切产物经1%琼脂糖凝胶电泳检测，在预期大小的位置出现了两条清晰的条带，一条为载体片段，大小约为2692bp，另一条为插入的目的片段，大小与预测的microRNA前体大小一致，表明重组质粒构建成功。将重组质粒送测序公司进行测序，测序结果与预测的玉米microRNA前体序列进行比对，结果显示两者完全一致，进一步验证了预测得到的玉米microRNA前体序列的准确性。3.3实时荧光定量PCR法实时定量检测microRNA结果对分离得到的小分子RNA进行加尾和引物延伸反应，成功获得了用于实时荧光定量PCR的cDNA模板。以U6snRNA为内参基因，对预测得到的玉米miRNA进行实时荧光定量PCR扩增。扩增结果显示，各miRNA均能得到特异性扩增，熔解曲线分析表明，各miRNA的扩增产物熔解曲线均为单一峰，无杂峰出现，说明扩增特异性良好，未出现非特异性扩增和引物二聚体等问题。通过标准曲线法对各miRNA进行定量分析，以不同浓度梯度的已知拷贝数的miRNA标准品进行实时荧光定量PCR扩增，绘制标准曲线。结果显示，内参基因U6和各miRNA的标准曲线均具有良好的线性关系，相关系数R²均大于0.99，扩增效率在90%-110%之间，满足实时荧光定量PCR的实验要求，能够准确地对样品中的miRNA进行定量分析。[此处插入图2：内参基因U6和各miRNA的标准曲线和熔解曲线，图中横坐标为log10(起始拷贝数)，纵坐标为Ct值，熔解曲线横坐标为温度(℃)，纵坐标为荧光信号的一阶导数][此处插入图2：内参基因U6和各miRNA的标准曲线和熔解曲线，图中横坐标为log10(起始拷贝数)，纵坐标为Ct值，熔解曲线横坐标为温度(℃)，纵坐标为荧光信号的一阶导数]对高抗纹枯病自交系R15和感病自交系黄早4接菌后不同时间的叶片样品中各miRNA的相对表达量进行分析。结果表明，在R15中，miR164在接菌后12h表达量开始显著上调，在24h达到峰值，随后逐渐下降，但仍高于对照水平；miR167在接菌后24h表达量显著上调，48h时达到最高，之后略有下降；而在黄早4中，miR164和miR167的表达量在接菌后虽有变化，但幅度相对较小，且与对照相比，差异不显著。miR156在R15中接菌后表达量呈先下降后上升的趋势，在24h时表达量最低，随后逐渐升高；在黄早4中，其表达量变化不明显。这些结果表明，不同miRNA在玉米抗纹枯病过程中的表达模式存在差异，且在抗病自交系和感病自交系中的表达变化也有所不同，暗示它们在玉米抗纹枯病机制中可能发挥着不同的作用。[此处插入图3：各miRNA在R15和黄早4接菌后的不同时间的相对表达量变化图，图中横坐标为接菌后的时间(h)，纵坐标为miRNA的相对表达量，以未接菌的对照样品为1，误差线表示标准误，*表示差异显著(P<0.05)，**表示差异极显著(P<0.01)][此处插入图3：各miRNA在R15和黄早4接菌后的不同时间的相对表达量变化图，图中横坐标为接菌后的时间(h)，纵坐标为miRNA的相对表达量，以未接菌的对照样品为1，误差线表示标准误，*表示差异显著(P<0.05)，**表示差异极显著(P<0.01)]进一步对各miRNA的靶基因进行预测，利用psRNATarget软件，共预测得到了多个潜在的靶基因。对这些靶基因进行功能注释分析发现，miR164的靶基因主要包括一些NAC转录因子基因，这些转录因子在植物的生长发育、细胞分化以及逆境响应等过程中发挥重要作用；miR167的靶基因包括生长素响应因子ARF6和ARF8，参与生长素信号传导途径，调控植物的生长发育和对逆境的响应；miR156的靶基因主要为SBP-box转录因子基因，该类转录因子在植物的生长发育、开花调控以及对逆境胁迫的响应中具有重要功能。这些靶基因的功能与植物的抗病过程密切相关，进一步说明了预测得到的miRNA在玉米抗纹枯病机制中可能通过调控这些靶基因的表达来发挥作用。四、讨论4.1miRNA生物信息学鉴定的优缺点在本研究中，生物信息学方法在鉴定玉米抗纹枯病相关miRNA方面展现出显著优势。借助生物信息学工具，能够快速处理海量的基因组数据，从玉米庞大的基因组中高效筛选出潜在的miRNA序列。与传统的实验鉴定方法相比，大大缩短了研究周期，降低了研究成本。在预测miRNA靶基因时，psRNATarget等生物信息学软件能够根据miRNA与靶基因的互补配对原则，快速预测出潜在的靶基因，为后续研究miRNA的调控机制提供了重要线索。然而，生物信息学鉴定方法也存在一定的局限性。在miRNA预测过程中，基于序列相似性和结构特征的预测方法可能会产生假阳性结果。由于miRNA的结构和序列存在一定的灵活性，一些并非真正miRNA的序列可能因具有类似的结构特征而被误判为潜在的miRNA。例如，某些非编码RNA片段或mRNA降解产物，其二级结构可能与miRNA前体的发夹结构相似，从而干扰了miRNA的准确预测。在靶基因预测方面，生物信息学方法虽然能提供大量潜在的靶基因，但这些预测结果的准确性仍有待进一步验证。不同的靶基因预测软件由于算法和参数设置的差异，预测结果可能存在较大差异，导致难以确定真正的靶基因。并且，生物信息学预测往往忽略了基因表达的时空特异性以及miRNA与靶基因之间复杂的相互作用环境，这些因素都可能影响预测结果的可靠性。为提高miRNA预测的准确性，可采取多种策略。一方面，整合多组学数据进行分析是有效的途径。结合转录组数据，可了解基因的表达情况，优先筛选在纹枯病菌侵染后表达有显著变化的miRNA，提高预测的针对性；结合蛋白质组数据，能从蛋白质水平验证miRNA对靶基因的调控作用，进一步确认预测结果的可靠性。另一方面，优化生物信息学算法也是关键。研发更精准的算法，综合考虑miRNA的多种特征，如序列保守性、结构稳定性、与靶基因的互补配对模式以及在不同组织和发育阶段的表达特异性等，减少假阳性和假阴性结果的出现。在后续研究中，还需加强生物信息学预测结果与实验验证的结合，通过实验手段对预测得到的miRNA及其靶基因进行验证，进一步完善玉米抗纹枯病相关miRNA的研究。4.2miR168与miR168~*的真实性讨论在本研究预测得到的玉米miRNA中，miR168及其互补链miR168~*备受关注。从序列特征来看，miR168成熟序列长度为21个核苷酸，其5’端具有典型的磷酸基团，3’端为羟基，符合miRNA的基本结构特征。通过Mfold和RNAfold软件预测其前体的二级结构，发现miR168前体能够形成稳定的发夹结构，茎部区域碱基互补配对程度高，最小自由能较低，这进一步支持了其作为miRNA的可能性。在其他植物中，如拟南芥、水稻等，miR168也被广泛研究并证实具有重要的生物学功能。在拟南芥中，miR168通过靶向AGO1基因，调控AGO1蛋白的表达水平，进而影响植物的生长发育和对逆境的响应。在水稻中，miR168同样参与了对稻瘟病菌侵染的免疫反应，其表达量在病菌侵染后发生显著变化。在玉米抗纹枯病过程中，miR168的表达模式显示出与抗病性的密切关联。实时荧光定量PCR结果表明，在高抗纹枯病自交系R15中，接菌后miR168的表达量迅速上调，在24h时达到峰值，随后逐渐下降，但仍维持在较高水平；而在感病自交系黄早4中，miR168的表达量虽有变化，但幅度相对较小，且显著低于R15中的表达水平。这表明miR168的高表达可能与玉米对纹枯病的抗性增强相关。其作用机制可能是通过靶向相关靶基因，调控植物的免疫反应。根据psRNATarget软件预测，miR168的潜在靶基因包括一些参与植物信号传导和防御反应的基因，如丝裂原活化蛋白激酶（MAPK）基因家族成员。MAPK级联反应在植物应对生物胁迫过程中发挥关键作用，通过激活下游的转录因子和防御基因，增强植物的抗病能力。miR168可能通过调控这些靶基因的表达，参与玉米对纹枯病的防御反应。关于miR168~*，其与miR168来源于同一前体的不同臂，在生物合成过程中，两者的产生量和稳定性可能存在差异。在本研究中，通过深度测序数据分析发现，miR168~*的表达量相对较低，但其在纹枯病菌侵染后的表达变化趋势与miR168具有一定的相关性。在R15中，miR168~*的表达量在接菌后也呈现出先上升后下降的趋势，尽管变化幅度不如miR168明显，但这暗示着miR168~*可能在玉米抗纹枯病过程中也发挥着一定的作用。一些研究表明，miRNA的互补链在某些情况下也能与AGO蛋白结合形成RISC，参与基因表达的调控。虽然miR168~*的具体功能尚不完全明确，但它可能与miR168协同作用，共同调节玉米对纹枯病的抗性。例如，miR168~*可能通过与miR168竞争结合靶基因，或者与其他RNA分子相互作用，影响miR168的调控网络，从而在玉米抗纹枯病过程中发挥潜在的调节作用。4.3上调表达的miRNA与玉米纹枯病的关系在本研究中，通过实时荧光定量PCR分析发现，在玉米接种纹枯病菌后，部分miRNA呈现上调表达的趋势。其中，miR164在高抗纹枯病自交系R15中，接菌后12h表达量开始显著上调，在24h达到峰值，随后逐渐下降，但仍高于对照水平。miR164的靶基因主要为NAC转录因子基因。NAC转录因子在植物的生长发育、细胞分化以及逆境响应等过程中发挥着至关重要的作用。在植物应对病原菌侵染时，NAC转录因子可通过调控下游防御基因的表达，参与植物的免疫反应。研究表明，某些NAC转录因子能够激活植物细胞壁相关防御基因的表达，增强细胞壁的强度和稳定性，从而限制病原菌的入侵和扩展。miR164可能通过靶向调控NAC转录因子基因的表达，参与玉米对纹枯病的防御反应。在纹枯病菌侵染初期，miR164表达量上调，可能导致NAC转录因子基因的表达受到抑制，从而解除对下游防御基因的抑制作用，使防御基因得以表达，增强玉米的抗病性。但随着侵染时间的延长，miR164表达量下降，NAC转录因子基因的表达逐渐恢复，可能参与调控植物的生长发育及其他生理过程，以维持植物在逆境下的正常生长。miR167在R15中接菌后24h表达量显著上调，48h时达到最高，之后略有下降。其靶基因包括生长素响应因子ARF6和ARF8，参与生长素信号传导途径。生长素在植物的生长发育过程中起着关键作用，同时也参与植物对逆境的响应。在植物受到病原菌侵染时，生长素信号传导途径会发生改变，以调节植物的免疫反应。ARF6和ARF8作为生长素响应因子，可与生长素响应基因的启动子区域结合，调控基因的表达。miR167可能通过调控ARF6和ARF8的表达，影响生长素信号传导途径，进而参与玉米对纹枯病的抗性。在纹枯病菌侵染后，miR167表达量上调，可能抑制ARF6和ARF8的表达，从而调节生长素信号传导，激活植物的防御反应，增强玉米对纹枯病的抗性。同时，生长素信号的改变也可能影响植物的生长发育，如抑制植株的纵向生长，使更多的能量和物质分配到防御过程中，以应对纹枯病菌的侵染。这些上调表达的miRNA在玉米抗纹枯病过程中发挥着重要作用，它们通过调控各自靶基因的表达，参与植物的免疫反应、信号传导以及生长发育等过程，共同调节玉米对纹枯病的抗性。不同miRNA及其靶基因之间可能存在复杂的相互作用网络，协同调控玉米的抗病机制。在后续研究中，深入探究这些miRNA与靶基因之间的具体作用机制，以及它们在玉米抗纹枯病网络中的相互关系，将有助于进一步揭示玉米抗纹枯病的分子机制，为玉米抗病育种提供更坚实的理论基础。4.4下调表达的miRNA与玉米纹枯病的关系在玉米接种纹枯病菌的实验中，一些miRNA呈现出下调表达的趋势，这表明它们在玉米抗纹枯病过程中可能发挥着独特的作用。以miR156为例，在高抗纹枯病自交系R15中，接菌后其表达量呈先下降后上升的趋势，在24h时表达量最低，随后逐渐升高；而在感病自交系黄早4中，其表达量变化不明显。miR156的靶基因主要为SBP-box转录因子基因，该类转录因子在植物的生长发育、开花调控以及对逆境胁迫的响应中具有重要功能。在植物应对纹枯病菌侵染时，细胞壁是抵御病原菌入侵的第一道防线。研究表明，某些SBP-box转录因子可通过调控细胞壁合成相关基因的表达，增强细胞壁的强度和稳定性，从而限制纹枯病菌的侵入。在纹枯病菌侵染初期，R15中miR156表达量下调，可能导致SBP-box转录因子基因的表达上调，进而促进细胞壁合成相关基因的表达，增强细胞壁的防御能力，抵御纹枯病菌的侵染。随着侵染时间的延长，miR156表达量逐渐上升，对SBP-box转录因子基因的抑制作用逐渐恢复，可能是植物在平衡防御反应与正常生长发育之间的关系。蛋白质可逆磷酸化是植物细胞内重要的信号传导机制之一，参与植物对多种逆境胁迫的响应。有研究发现，SBP-box转录因子可通过与蛋白激酶或蛋白磷酸酶相互作用，参与蛋白质可逆磷酸化过程，进而调节植物的免疫反应。miR156通过调控SBP-box转录因子基因的表达，可能间接影响蛋白质可逆磷酸化信号通路，在玉米抗纹枯病过程中发挥调控作用。当纹枯病菌侵染时，miR156表达量的变化可能导致SBP-box转录因子的表达改变，进而影响相关蛋白激酶和蛋白磷酸酶的活性，调节下游防御基因的表达，增强玉米的抗病性。植物自身调控是一个复杂的网络，miR156及其靶基因在其中也扮演着重要角色。SBP-box转录因子不仅参与植物的生长发育调控，还在植物激素信号传导、氧化还原平衡调节等过程中发挥作用。在玉米抗纹枯病过程中，miR156表达量的变化可能通过调控SBP-box转录因子，影响植物自身的调控网络，使植物能够更好地应对纹枯病菌的侵染。miR156可能通过调节SBP-box转录因子，影响植物激素的合成、运输和信号传导，从而调节植物的生长发育和免疫反应。谷胱甘肽-S-转移酶（GST）是植物体内一类重要的解毒酶，参与植物对多种逆境胁迫的响应，包括病原菌侵染。研究发现，SBP-box转录因子可调控GST基因的表达，增强植物的抗氧化能力和解毒能力，从而提高植物的抗病性。在玉米抗纹枯病过程中，miR156表达量的下调可能导致SBP-box转录因子基因表达上调，进而促进GST基因的表达，增强玉米对纹枯病菌侵染所产生的氧化胁迫的抵抗能力，减轻病菌对植物细胞的损伤。激素在植物的生长发育和逆境响应中起着关键作用。SBP-box转录因子可参与生长素、赤霉素、脱落酸等多种激素的信号传导途径，调节植物的生长发育和对逆境的响应。在玉米抗纹枯病过程中，miR156可能通过调控SBP-box转录因子，影响激素信号传导，从而调节玉米的抗病性。在纹枯病菌侵染时，miR156表达量的变化可能导致SBP-box转录因子对激素信号传导途径的调节发生改变，使植物能够调整自身的生长发育和防御反应，以适应逆境环境。诱导相关蛋白是植物在受到病原菌侵染后产生