花椰菜小孢子发育进程中miRNAs的鉴定、特征及调控机制解析

花椰菜小孢子发育进程中miRNAs的鉴定、特征及调控机制解析一、引言1.1miRNA研究概述1.1.1miRNA的发现历程miRNA的发现源于对线虫发育的研究。1993年，Lee等人在秀丽隐杆线虫中首次发现lin-4基因，其转录产物并非编码蛋白质，而是一种长度约22个核苷酸的小分子RNA，它能够通过与lin-14mRNA的3'-UTR区域互补配对，抑制lin-14蛋白的表达，从而调控线虫的发育时序，不过当时这一发现并未引起科学界的广泛关注。直到2000年，Reinhart等人在线虫中又发现了第二个具有时序调控功能的小分子RNA——let-7，它在不同物种间高度保守，且在发育进程中发挥关键作用，这才使得miRNA逐渐进入科研人员的视野。此后，随着研究技术的不断进步，多个研究小组在多种生物中陆续发现了大量类似的小分子RNA。2001年，这些具有相似特征和功能的小分子RNA被统一命名为miRNA（microRNA），从此掀起了miRNA研究的热潮。同年，Science杂志将miRNA评选为年度十大科技突破之一，标志着miRNA研究领域的正式确立。在植物中，2002年美国的俄勒冈州立大学、Rutgers大学、麻省理工学院与Rice大学，以及奥地利科学院的四个研究小组分别报道从植物中发现了miRNA，开启了植物miRNA研究的新篇章。此后，越来越多的植物miRNA被鉴定出来，人们对其在植物生长发育、应对环境胁迫等方面的功能和作用机制也有了更深入的认识。1.1.2植物miRNA特点剖析植物miRNA具有独特的结构特征，其长度一般为20-24个核苷酸，成熟的miRNA5'端有一磷酸基团，3'端为羟基。与动物miRNA相比，植物miRNA前体的茎环结构更大、更复杂，长度变异范围在64-303nt，远大于动物miRNA前体（60-70nt）。在长度方面，植物miRNA长度多为21nt，而动物miRNA长度多为22-23nt。从碱基偏好性来看，植物miRNA5'端对尿嘧啶（U）有强烈的倾向性，热力学分析表明，这种末端不稳态是通过RISC来维持的，并且植物中miRNA3'末端2nt突出的3'-OH存在甲基化修饰，而动物中无此修饰。植物miRNA还具有较高的进化保守性，一些保守的miRNA家族在不同植物物种中广泛存在，且序列差异较小，这使得对植物miRNA目标基因的预测相对简单。此外，植物miRNA在表达上具有明显的组织特异性和发育阶段特异性，在植物不同组织和不同发育时期，miRNA的表达谱存在显著差异，暗示它们在植物生长发育的特定过程中发挥着重要的调控作用。1.1.3植物miRNA生物合成路径植物miRNA的生物合成起始于细胞核，编码miRNA的基因在RNA聚合酶II的作用下转录形成初级转录本（pri-miRNA），pri-miRNA长度可达几百到上千个核苷酸，具有复杂的二级结构。随后，在Dicer-like1（DCL1）、HYL1（HYPONASTICLEAVES1）和SE（SERRATE）等蛋白形成的复合体作用下，pri-miRNA被逐步切割加工。首先，DCL1将pri-miRNA切割产生含有茎环结构的前体miRNA（pre-miRNA），接着DCL1和HYL1继续作用，将pre-miRNA进一步加工为miRNA/miRNA*双链体。加工完成后，miRNA/miRNA双链体的3'末端会被甲基转移酶HEN1（HUAENHANCER1）甲基化修饰，这一修饰能够保护miRNA不被降解，增强其稳定性。随后，甲基化的miRNA/miRNA双链体被转运出细胞核，这一过程由exportin-5的同系物HASTY（HST）介导。在细胞质中，miRNA/miRNA双链体中的一条链（通常为miRNA链）被降解，另一条链（成熟miRNA）则与AGO（ARGONAUTE）蛋白结合，形成RNA诱导沉默复合体（RISC），进而行使对靶基因的调控功能。1.1.4miRNA在植物中的作用机理miRNA在植物中的作用主要通过与靶mRNA的碱基互补配对来实现。植物miRNA与其靶mRNA序列通常具有高度的互补性，当miRNA与靶mRNA结合后，主要通过两种转录后机制下调基因表达：一是mRNA切割，若miRNA与靶mRNA的互补程度较高，RISC中的AGO蛋白会切割靶mRNA，导致其降解，从而抑制基因表达；二是翻译抑制，当miRNA与靶mRNA存在一定程度的错配时，虽然不会切割靶mRNA，但会抑制其翻译过程，使蛋白质无法正常合成。与动物miRNA不同，植物miRNA的结合位点大多位于靶mRNA的编码区，而动物miRNA的结合位点主要在3'UTR区域。植物miRNA通过精确识别靶mRNA序列，对植物基因表达进行精细调控，在植物生长发育、激素信号转导、应对生物和非生物胁迫等过程中发挥关键作用。1.1.5miRNA在植物中的生物学功能在植物生长发育方面，miRNA参与调控植物的各个发育阶段。以拟南芥为例，miR156通过调控SPL（SQUAMOSAPROMOTERBINDINGPROTEIN-LIKE）基因家族的表达，影响植物的营养生长向生殖生长转变，低水平的miR156促使植物进入生殖生长阶段，开花结果；miR164则参与调控植物的器官发育，如通过调控NAC1基因影响侧根的形成。在应对环境胁迫时，miRNA也发挥着重要作用。在干旱胁迫下，植物体内一些miRNA的表达发生变化，如miR169通过抑制NF-YA5（NUCLEARFACTOR-Y,SUBUNITA5）基因的表达，调控植物对干旱的响应，增强植物的耐旱性；在盐胁迫下，miR393通过靶向调控生长素受体基因TIR1等，参与植物对盐胁迫的适应过程，调节植物生长发育以应对盐害。此外，miRNA在植物抵御病虫害入侵时也起到关键的免疫调控作用，参与植物与病原菌互作过程中的信号传导和防御反应。1.1.6miRNA研究方法汇总高通量测序技术是目前鉴定和分析miRNA的重要手段，通过对小RNA文库进行大规模测序，能够全面、快速地获取样本中的miRNA信息，不仅可以发现新的miRNA，还能对已知miRNA的表达水平进行定量分析，揭示其在不同组织、不同发育阶段或不同处理条件下的表达谱变化。茎环PCR技术则常用于miRNA的定量检测，该方法利用茎环引物对miRNA进行逆转录，再通过PCR扩增进行定量分析，具有特异性强、灵敏度高的特点，能够准确检测低丰度的miRNA。荧光定量PCR是另一种常用的miRNA定量方法，它基于实时荧光监测PCR扩增过程，能够精确测定miRNA的表达量，操作简便、结果可靠，广泛应用于miRNA表达水平的验证和分析。除了上述方法，生物信息学预测也是研究miRNA的重要辅助手段。通过构建相关算法和模型，利用已知的miRNA序列和结构特征，预测基因组中潜在的miRNA及其靶基因，为实验研究提供线索和方向，大大提高了miRNA研究的效率和针对性。1.2花椰菜小孢子研究的意义1.2.1花椰菜的研究价值花椰菜（BrassicaoleraceaL.var.botrytisL.），作为十字花科芸薹属甘蓝种的一个变种，在全球蔬菜产业中占据重要地位。其起源于地中海东部的克里特岛，19世纪中叶由英国传入我国，经过多年的引种和栽培，如今我国花椰菜种植面积已跃居世界首位。花椰菜凭借其独特的营养价值备受青睐，它富含蛋白质、脂肪、碳水化合物、食物纤维以及多种维生素（如维生素A、B、C、E、P、U）和钙、铁、磷等矿物质。从营养角度来看，花椰菜含水量高达90%以上，热量较低，是减肥人士理想的饱腹食材。其中的“索弗拉芬”能刺激细胞制造保护酶，这种酶具有强大的抗癌活性，可使细胞形成对抗致癌物侵蚀的膜，对预防乳腺癌、前列腺癌、胃癌、肺癌等多种癌症有积极作用。同时，花椰菜还能增强肝脏解毒能力，对高血压和心脏病有调节和预防功效。松花菜作为花椰菜的一个类型，不仅叶绿素和可溶性糖含量更高，维生素C含量也普遍高于紧花型花椰菜，口感脆甜、清香爽口，久煮不烂，深受消费者喜爱，其加工产品如“松花菜干”市场需求旺盛。从经济价值方面，花椰菜的种植、加工和销售形成了庞大的产业链，为农业经济发展做出重要贡献，无论是鲜食市场还是加工产业，花椰菜都具有广阔的市场前景和发展潜力。1.2.2小孢子胚胎诱导对花椰菜育种的重要性小孢子胚胎诱导技术在花椰菜育种中具有不可替代的重要作用。传统的花椰菜育种方法主要依赖于常规杂交和选择，周期长、效率低。而小孢子胚胎诱导技术能够直接将小孢子培养成单倍体植株，再经过染色体加倍处理，快速获得纯合的双单倍体（DH）植株。这种方法大大缩短了育种周期，相比传统育种方法，可节省数年时间，使新品种的选育更加高效。通过小孢子胚胎诱导获得的DH系，其遗传背景高度纯合，在育种过程中能够稳定遗传，减少了后代性状分离，提高了选择效率。利用小孢子胚胎诱导技术，育种家可以在短时间内获得大量的遗传稳定的材料，便于筛选具有优良性状（如抗病性、抗逆性、高品质等）的植株，加速花椰菜品种改良进程。例如，在抗病虫害育种中，通过小孢子培养结合抗病筛选，可以快速获得携带抗病基因的纯合植株，为培育抗病新品种奠定基础。该技术还能打破不良基因连锁，创造新的遗传变异，拓宽花椰菜的遗传基础，为培育出更加优质、高产、多抗的花椰菜新品种提供了有力支持。1.2.3研究花椰菜小孢子发育中miRNA的意义从分子层面深入研究花椰菜小孢子发育中miRNA具有深远意义。miRNA作为基因表达的重要调控因子，在花椰菜小孢子发育过程中发挥着关键作用。通过鉴定花椰菜小孢子发育不同阶段的miRNA及其靶基因，能够揭示miRNA介导的基因调控网络，从而深入理解小孢子发育的分子机制。这不仅有助于我们从本质上认识植物生殖发育过程，也为植物发育生物学理论研究提供新的思路和证据。在花椰菜育种实践中，明确miRNA在小孢子发育中的调控作用，可为分子标记辅助育种提供理论依据。通过对关键miRNA及其靶基因的分析，可以开发出与优良性状紧密连锁的分子标记，实现对目标性状的精准选择，提高育种效率和准确性。对miRNA的研究还有助于挖掘新的基因资源，为花椰菜遗传改良提供新的靶点，通过基因编辑等技术手段，对miRNA及其靶基因进行调控，有望培育出具有特殊优良性状的花椰菜新品种，推动花椰菜育种技术的创新和发展，满足市场对多样化、高品质花椰菜品种的需求。二、材料与方法2.1花椰菜小孢子样本获取2.1.1花椰菜种植与取材选用生长势强、抗病性好的花椰菜品种（如“雪宝花椰菜”）种子，播种于装有育苗基质（草炭：蛭石=3:1，v/v）的育苗盘中，每穴播1-2粒种子，播后覆盖0.5-1cm厚的基质，并浇透水。育苗期间，保持白天温度在20-25℃，夜间温度15-18℃，光照强度200-300μmol・m⁻²・s⁻¹，光照时间12-14h/d，空气相对湿度60%-70%。待幼苗长至4-5片真叶时，移栽至日光温室中，种植密度为行距50cm，株距40cm。在花椰菜生长过程中，定期浇水、施肥，以保证植株的正常生长。施肥遵循“基肥为主，追肥为辅”的原则，基肥每667m²施入腐熟有机肥3000-5000kg、三元复合肥（N:P:K=15:15:15）30-50kg；追肥在莲座期、花球形成期和膨大期进行，每次每667m²追施尿素10-15kg、硫酸钾5-10kg。为获取不同发育时期的小孢子，在花椰菜抽薹开花期，根据花蕾形态和大小初步判断小孢子发育时期。一般来说，四分体时期的花蕾长度约为1-2mm，单核期花蕾长度为2-3mm，双核期花蕾长度为3-4mm。选取不同长度的花蕾，用镊子小心地从植株上取下，放入装有冰盒的样品袋中，迅速带回实验室进行后续处理。每个发育时期选取30-50个花蕾，重复3次，以保证样本的代表性。2.1.2小孢子的分离与纯化将采集的花蕾用流水冲洗干净，去除表面杂质，然后放入75%乙醇中浸泡30-60s进行表面消毒，再用无菌水冲洗3-5次。将消毒后的花蕾置于无菌培养皿中，加入适量的B5液体培养基（含3%蔗糖，pH5.8），用镊子和剪刀将花蕾打开，取出花药。采用机械分离法将小孢子从花药中释放出来，将花药转移至装有2-3mLB5液体培养基的无菌小烧杯中，用玻璃棒轻轻挤压花药，使小孢子游离到培养基中。接着，用350目孔径的尼龙筛网过滤小孢子悬浮液，去除未破碎的花药组织和杂质，将滤液收集到离心管中。利用密度梯度离心法进一步纯化小孢子，将Percoll试剂与B5液体培养基按一定比例混合，配制成浓度分别为30%、50%和70%的Percoll密度梯度溶液。将过滤后的小孢子悬浮液小心地铺在30%Percoll溶液表面，然后在4℃、1000g条件下离心20min。离心后，小孢子会在不同密度的Percoll溶液界面形成条带，用移液器小心吸取位于50%Percoll溶液界面的小孢子层，转移至新的离心管中。将收集到的小孢子用B5液体培养基洗涤2-3次，每次在4℃、1000g条件下离心5-10min，去除残留的Percoll试剂和杂质。最后，将纯化后的小孢子悬浮于适量的B5液体培养基中，调整小孢子密度至1×10⁵-5×10⁵个/mL，用于后续实验。2.2总RNA提取及质量检测2.2.1Trizol法提取总RNA采用Trizol试剂提取花椰菜不同发育时期小孢子的总RNA，具体操作步骤如下：样本准备：取适量纯化后的小孢子悬浮液，放入1.5mL的EP管中，在4℃条件下，以12000rpm的转速离心5min，小心弃去上清液，尽可能吸干残留液体。若使用匀浆仪处理样本，悬浮液可不经过离心步骤，直接取适量悬浮液加入0.5mLTrizol试剂，在冰上充分匀浆，确保小孢子完全裂解，随后再加入0.5mLTrizol试剂冲洗匀浆器，并将冲洗液转移至1.5mLEP管中。裂解与分离：向装有小孢子沉淀的EP管中加入1mLTrizol试剂，使用移液器充分吹打混匀，使小孢子与Trizol试剂充分接触，在15-30℃条件下放置5min，促使核酸蛋白复合物完全分离。若样品中含有较多蛋白、脂肪、多糖或类似肌肉、植物结节等特殊成分，可在4℃、12000rpm条件下离心10min，取上清液用于后续操作，离心得到的沉淀主要包含细胞外膜、多糖、高分子量DNA等物质，而上清液中含有RNA。在处理脂肪组织样品时，需注意去除上层大量的油脂，取澄清的匀浆溶液继续下一步实验。氯仿抽提：每使用1mLTrizol试剂，加入0.2mL三氯甲烷（氯仿），盖紧管盖后，在漩涡振荡器上剧烈振荡15s，确保溶液充分混合，然后室温放置3min。若无法使用漩涡振荡器，也可手动颠倒混匀2min代替，使水相和有机相充分接触，促进RNA从有机相中分离。离心分层：将样品在4℃条件下，以12000rpm的转速离心10-15min，此时样品会明显分成三层，从上到下依次为无色的水相、中间层和红色的有机相，RNA主要存在于水相中。小心吸取水相（约600μL，约为所用Trizol试剂体积的60%）转移到新的EP管中，注意避免吸取到中间层和有机相，以免造成RNA污染。异丙醇沉淀：在得到的水相溶液中加入等体积（约500μL）的异丙醇，上下颠倒EP管混匀，使RNA充分沉淀，将EP管置于-20℃冰箱中放置20-30min，以增强RNA的沉淀效果。在植物或动物组织RNA提取过程中，为减轻蛋白等污染物沉淀对RNA检测的影响，可在上清液中加入1/2体积的异丙醇和1/2体积的RNA助沉剂，然后进行后续操作。离心收集沉淀：将EP管从-20℃冰箱取出，在4℃条件下，以12000rpm的转速离心10min，此时RNA会沉淀在管底和管壁，形成胶状沉淀。小心去除上清液，注意不要丢失RNA沉淀。乙醇洗涤：向含有RNA沉淀的EP管中加入1mL75%乙醇（用DEPC处理过的水配制），轻轻敲击EP管，使RNA沉淀悬浮起来，充分洗涤沉淀，去除残留的杂质和盐分。每使用1mLTrizol试剂，至少加入1mL乙醇进行洗涤。有时为避免RNA被洗掉，此步洗涤操作可以省略，但需注意后续晾干或烤干乙醇时，不能过于干燥，否则RNA不易溶解。在4℃条件下，以12000rpm的转速离心5min，弃去上清液，短暂快速离心后，用移液器小心吸弃残余上清液，确保不吸弃RNA沉淀。晾干与溶解：将EP管置于室温下放置晾干，注意不要晾得过干，一般晾干2-3min左右即可，过度干燥会导致RNA很难溶解。根据实验需要，向EP管中加入30-100μLDEPC水或0.5%SDS溶液，使用移液器枪头反复吸打几次，使RNA充分溶解。在整个操作过程中，需要严格遵守无菌操作规范，操作人员需佩戴口罩、手套，防止外源RNA酶污染。实验涉及的离心管、Tip头、移液器杆、电泳槽、实验台面等都要进行彻底处理，确保无RNA酶残留。实验所使用的试剂/溶液，尤其是水，必须确保为RNase-Free级别，以保证RNA的完整性和纯度。2.2.2RNA质量检测方法采用琼脂糖凝胶电泳和分光光度计对提取的总RNA进行完整性和纯度检测。琼脂糖凝胶电泳检测完整性：RNA电泳可在变性及非变性两种条件下进行，在快速检测总RNA样品完整性时，通常配制普通的1%琼脂糖凝胶即可。以1×TAE电泳缓冲液制作琼脂糖凝胶，待凝胶凝固后，将其放入电泳槽中，并加入1×TAE电泳缓冲液至液面完全覆盖凝胶。在超净工作台上，用移液器吸取4μL总RNA样品滴于封口膜上，再加入5μL1×TAE电泳缓冲液及1μL的10X载样缓冲液，充分混匀后，小心加入凝胶的点样孔中。打开电源开关，调节电压至100V，使RNA由负极向正极进行电泳，约30min后，将凝胶放入0.5μg/mL溴化乙锭（EB）染液中染色5min，然后用清水稍微漂洗。在紫外透射检测仪上观察RNA电泳结果，完整的未降解的RNA制品，其电泳图谱应能清晰看到28SrRNA、18SrRNA、5SrRNA的三条带，且28SrRNA的亮度约为18SrRNA的两倍，若RNA发生降解，则条带会出现弥散、模糊或条带缺失等现象，表明RNA的完整性受到破坏。分光光度计检测纯度：利用分光光度计测定RNA溶液在260nm和280nm波长处的吸光度（A值），根据A260/A280的比值来判断RNA的纯度。纯净的RNA溶液，其A260/A280比值应在1.8-2.0之间。若比值低于1.8，可能存在蛋白质或酚类物质污染；若比值高于2.0，则可能存在RNA降解或有残留的盐离子、胍盐等杂质。同时，还可根据A260的值计算RNA的浓度，计算公式为：RNA浓度（μg/μL）=A260×稀释倍数×40/1000，通过浓度和纯度的测定，可全面评估提取的RNA质量是否符合后续实验要求。2.3小RNA高通量测序2.3.1文库构建流程小RNA高通量测序文库的构建是获取高质量测序数据的关键步骤，其流程涉及多个精细操作环节。在完成总RNA提取及质量检测后，选取质量合格的RNA样品进行文库构建。首先，利用T4RNA连接酶将3'接头（3'adaptor）连接到小RNA的3'末端，该接头为后续的逆转录和PCR扩增提供了特异性引物结合位点。在连接过程中，需要精确控制反应体系中各成分的比例和反应条件，一般反应体系包含小RNA样品、3'接头、T4RNA连接酶缓冲液、ATP以及适量的T4RNA连接酶，在25℃条件下反应1-2h，以确保接头与小RNA高效连接。随后，使用T4RNA连接酶将5'接头（5'adaptor）连接到已接上3'接头的小RNA的5'末端，同样需要严格控制反应条件，反应体系组成与3'接头连接类似，只是接头换成5'接头，在合适的温度和时间条件下进行反应，完成小RNA两端接头的连接。接着，以连接好接头的小RNA为模板，利用逆转录酶（如M-MLV逆转录酶）和特异性引物进行逆转录反应，将小RNA逆转录成cDNA。逆转录反应体系包括小RNA模板、逆转录酶、逆转录缓冲液、dNTPs、引物以及适量的RNase抑制剂，在37℃条件下反应60-90min，然后通过70℃加热5-10min使逆转录酶失活，终止反应。最后，以逆转录得到的cDNA为模板，进行PCR扩增。PCR扩增体系包含cDNA模板、PCR引物（与接头序列互补）、TaqDNA聚合酶、PCR缓冲液以及dNTPs，通过设置合适的PCR扩增程序（如95℃预变性3-5min；95℃变性30-45s，55-60℃退火30-45s，72℃延伸30-45s，共进行15-25个循环；最后72℃延伸5-10min），使cDNA得到有效扩增，从而获得足够量的小RNA文库，用于后续的高通量测序。在整个文库构建过程中，每一步反应都需要进行严格的质量控制，如利用琼脂糖凝胶电泳检测连接产物和PCR扩增产物的大小和纯度，确保文库构建的准确性和有效性。2.3.2测序平台及参数本研究选用IlluminaHiSeq2500测序平台进行小RNA高通量测序，该平台在小RNA测序领域具有广泛应用，以其高准确性、高灵敏度和高通量的优势，能够满足对花椰菜小孢子发育过程中miRNA进行全面、深入分析的需求。IlluminaHiSeq2500测序平台采用边合成边测序（SBS）技术，测序读长为50bp，这一长度能够准确读取小RNA序列信息，同时保证测序的准确性和稳定性。在测序通量方面，该平台单次运行可产生高达600Gb的数据量，足以对花椰菜小孢子不同发育时期的小RNA进行深度测序，确保检测到低丰度的miRNA。通过对每个样本进行足够深度的测序，能够更全面地获取样本中的miRNA信息，减少漏检的可能性，为后续的数据分析和miRNA鉴定提供充足的数据支持。在测序过程中，严格按照平台的操作规范进行样本上机、测序反应和数据采集，确保测序数据的质量和可靠性，为揭示花椰菜小孢子发育过程中miRNA的表达谱和调控特征奠定坚实基础。2.4测序数据的生物信息学分析2.4.1数据预处理测序得到的原始数据中往往包含低质量reads、接头序列、污染序列等杂质，这些杂质会干扰后续的数据分析，降低分析结果的准确性和可靠性，因此需要对原始数据进行严格的预处理。首先，利用FastQC软件对原始测序数据进行质量评估，该软件能够从多个方面对数据质量进行分析，如碱基质量分布、序列长度分布、GC含量分布等，通过生成的质量报告，可以直观地了解数据中存在的问题。在质量评估的基础上，使用Trimmomatic软件进行数据过滤和修剪。设置参数，去除低质量的reads，通常将质量值低于20（Q20）的碱基所在的reads去除，因为低质量的碱基可能是测序错误导致的，会影响后续分析结果。同时，去除含有接头序列的reads，接头序列是在文库构建过程中引入的，对分析小孢子发育中的miRNA并无实际意义，去除接头序列可以避免其对数据分析的干扰。通过设置合适的参数，如滑动窗口大小、窗口内最低平均质量值等，精确去除接头序列和低质量碱基。对长度过短的reads进行过滤，一般将长度小于18nt的reads去除，因为过短的序列可能无法准确匹配到参考基因组或已知miRNA数据库，会增加分析的噪声。经过这些预处理步骤，能够得到高质量的cleanreads，为后续的生物信息学分析提供可靠的数据基础。2.4.2已知miRNA注释已知miRNA注释是数据分析的关键环节，通过将预处理后的cleanreads与已知miRNA数据库进行比对，可以准确识别出样本中的已知miRNA。目前常用的miRNA数据库有miRBase，它是一个全面且权威的miRNA数据库，包含了多个物种的miRNA序列和注释信息。使用Bowtie软件进行序列比对，该软件能够高效、准确地将reads映射到参考序列上。在比对过程中，设置合适的参数，如允许的错配碱基数、最大比对次数等，以确保比对结果的准确性。将cleanreads与miRBase数据库中的花椰菜miRNA序列进行比对。若reads能够与数据库中的miRNA序列完全匹配或存在少量错配（根据设置的错配参数确定），则认为该reads对应一个已知miRNA。通过比对结果，统计每个已知miRNA在不同样本中的表达量，表达量通常以每百万reads中的reads数（readspermillion，RPM）来表示。计算公式为：RPM=（某miRNA的reads数/总cleanreads数）×10^6。这样可以量化已知miRNA在花椰菜小孢子不同发育时期的表达水平，为后续分析miRNA的表达模式和功能提供数据支持。2.4.3未知miRNA预测在完成已知miRNA注释后，为了全面揭示花椰菜小孢子发育过程中的miRNA调控网络，需要预测样本中可能存在的未知miRNA。本研究选用miRDeep2软件进行未知miRNA预测，该软件基于最大似然法原理，通过对测序数据中的小RNA序列特征、二级结构信息以及与参考基因组的匹配情况等多方面因素进行综合分析，实现对新miRNA的准确预测。首先，将预处理后的cleanreads与花椰菜参考基因组进行比对，确定小RNA在基因组上的位置。miRDeep2软件会分析这些小RNA在基因组上的分布特征，寻找具有典型miRNA前体结构特征的区域，即能够形成稳定茎环结构的区域。通过计算小RNA序列在茎环结构中的分布、碱基配对情况以及热力学稳定性等参数，评估该区域是否为潜在的miRNA前体。如果一个区域的小RNA序列能够形成稳定的茎环结构，且其表达模式符合miRNA的特征（如在不同样本中具有一定的表达量变化，且表达量相对稳定），则该区域可能包含一个新的miRNA。miRDeep2软件还会对预测得到的新miRNA进行打分，分数越高表示该预测结果越可靠。设置合适的分数阈值，筛选出高可信度的新miRNA。对预测得到的新miRNA进行进一步的验证和分析，如通过实验手段（如Northernblot、茎环PCR等）验证其表达情况，以及分析其靶基因和潜在的生物学功能，从而全面深入地了解花椰菜小孢子发育过程中未知miRNA的调控作用。2.4.4差异表达分析差异表达分析旨在筛选出在花椰菜小孢子不同发育时期表达水平存在显著差异的miRNA，这些差异表达的miRNA可能在小孢子发育过程中发挥关键调控作用。首先，利用前面计算得到的已知miRNA和预测得到的新miRNA的表达量数据（以RPM表示），使用DESeq2等软件进行差异表达分析。DESeq2软件基于负二项分布模型，通过对样本间的表达量数据进行统计分析，评估每个miRNA在不同发育时期的表达差异是否具有统计学意义。在分析过程中，软件会考虑样本间的生物学重复、测序深度等因素，对表达量数据进行标准化处理，以消除这些因素对差异表达分析的影响。计算每个miRNA在不同发育时期的表达倍数变化（foldchange），即两个发育时期中miRNA表达量的比值。若某miRNA在一个发育时期的表达量是另一个发育时期的2倍及以上（foldchange≥2），则认为该miRNA在这两个时期存在显著的表达差异。DESeq2软件会计算每个miRNA表达差异的P值，P值反映了差异表达结果的可信度，P值越小，表示差异表达越显著。通常设置P值的阈值为0.05，即当P值小于0.05时，认为该miRNA在不同发育时期的表达差异具有统计学意义。为了进一步控制假阳性率，还会对P值进行多重检验校正，常用的方法是Benjamini-Hochberg校正。经过校正后得到的调整后P值（adjustedP-value，也称为FDR，falsediscoveryrate），当FDR小于0.05时，确定该miRNA为差异表达miRNA。通过这些严格的分析步骤，筛选出在花椰菜小孢子不同发育时期真正具有显著差异表达的miRNA，为后续深入研究miRNA在小孢子发育中的调控机制奠定基础。2.4.5靶基因预测与功能注释预测miRNA的靶基因是深入了解miRNA调控机制的重要环节，因为miRNA主要通过作用于靶基因来发挥其生物学功能。本研究选用psRNATarget等靶基因预测软件对差异表达的miRNA进行靶基因预测。psRNATarget软件基于植物miRNA与靶mRNA互补配对的特性，通过对miRNA序列和花椰菜基因组中mRNA序列进行比对分析，预测可能的靶基因。在预测过程中，软件会考虑miRNA与靶mRNA的互补配对程度、错配情况以及结合位点的位置等因素，对预测结果进行评估和筛选。通过psRNATarget软件，得到每个差异表达miRNA的潜在靶基因列表。对这些靶基因进行功能注释，首先将靶基因序列与GeneOntology（GO）数据库进行比对，GO数据库从生物过程、细胞组成和分子功能三个层面提供了基因功能的标准化描述。通过比对，可以获得靶基因在这三个层面的功能注释信息，了解靶基因参与的生物学过程（如细胞代谢、信号转导、发育调控等）、在细胞中的定位（如细胞核、细胞质、细胞膜等）以及所具有的分子功能（如酶活性、转录因子活性、结合活性等）。将靶基因序列与KyotoEncyclopediaofGenesandGenomes（KEGG）数据库进行比对，KEGG数据库整合了基因组、生物化学和系统功能信息，通过比对可以确定靶基因参与的代谢通路和信号转导途径，如植物激素信号转导通路、碳水化合物代谢通路、氧化磷酸化通路等。通过对靶基因的功能注释和代谢通路分析，能够系统地揭示差异表达miRNA在花椰菜小孢子发育过程中的调控网络和潜在生物学功能，为深入理解小孢子发育的分子机制提供全面的信息支持。2.5实验验证方法2.5.1茎环PCR扩增与重测序验证为确保通过生物信息学分析鉴定出的miRNA的真实性和准确性，采用茎环PCR扩增技术对miRNA进行特异性扩增，并通过重测序进一步验证。首先，依据miRNA的序列信息，设计特异性的茎环引物。茎环引物的设计遵循严格的原则，其茎部一般包含14-18个碱基，通过互补配对形成稳定的茎环结构，环部则包含6-8个碱基，能够与miRNA的3'末端特异性结合。以提取的总RNA为模板，在逆转录酶（如M-MLV逆转录酶）的作用下，利用茎环引物进行逆转录反应。反应体系包含总RNA模板、茎环引物、逆转录酶、逆转录缓冲液、dNTPs以及适量的RNase抑制剂，在37℃条件下反应60-90min，将miRNA逆转录为cDNA。以逆转录得到的cDNA为模板，使用与miRNA特异性结合的正向引物和通用反向引物进行PCR扩增。PCR扩增体系包括cDNA模板、引物、TaqDNA聚合酶、PCR缓冲液以及dNTPs。设置合适的PCR扩增程序，一般95℃预变性3-5min；95℃变性30-45s，55-60℃退火30-45s，72℃延伸30-45s，共进行30-35个循环；最后72℃延伸5-10min，使目的miRNA片段得到有效扩增。扩增结束后，利用琼脂糖凝胶电泳检测PCR产物，观察是否出现预期大小的特异性条带。若条带大小与理论值相符，则表明扩增成功。将扩增得到的PCR产物进行纯化，去除反应体系中的杂质和引物二聚体等。纯化方法可采用胶回收试剂盒，按照试剂盒说明书进行操作，通过切胶回收目的条带，并用洗脱缓冲液洗脱得到纯化的PCR产物。将纯化后的PCR产物送往专业测序公司进行重测序。测序公司利用先进的测序技术，对PCR产物的序列进行精确测定，并将测序结果与之前生物信息学分析预测的miRNA序列进行比对。若两者序列高度一致，仅有少量测序误差，且在合理范围内（如碱基错配率低于1%），则可验证所鉴定的miRNA的准确性，确保后续研究基于真实可靠的miRNA数据展开。2.5.2荧光实时定量PCR验证为了准确检测miRNA及其靶基因在花椰菜小孢子不同发育时期的表达量变化，采用荧光实时定量PCR技术。对于miRNA的检测，同样利用茎环引物进行逆转录反应，将miRNA逆转录为cDNA。逆转录体系和条件与茎环PCR扩增中的逆转录步骤一致。以逆转录得到的cDNA为模板，使用特异性的miRNA引物和通用引物进行荧光定量PCR扩增。荧光定量PCR扩增体系包含cDNA模板、引物、SYBRGreen荧光染料、TaqDNA聚合酶、PCR缓冲液以及dNTPs。在扩增过程中，SYBRGreen荧光染料会与双链DNA结合，随着PCR反应的进行，双链DNA不断扩增，荧光信号也随之增强，通过实时监测荧光信号的变化，可准确测定miRNA的表达量。设置合适的荧光定量PCR扩增程序，95℃预变性3-5min；95℃变性10-15s，60℃退火30-45s，共进行40-45个循环。在每个循环的退火阶段采集荧光信号，通过分析荧光信号的变化曲线，利用软件（如LightCycler480Software）计算出每个样本中miRNA的相对表达量。采用2^(-ΔΔCt)法进行计算，其中ΔCt=Ct(miRNA)-Ct(内参基因)，ΔΔCt=ΔCt(实验组)-ΔCt(对照组)，内参基因可选择U6snRNA等在不同发育时期表达相对稳定的基因，通过与内参基因的比较，消除样本间RNA上样量和逆转录效率等差异对结果的影响，得到miRNA的相对表达量，从而准确反映miRNA在不同发育时期的表达变化情况。对于靶基因的检测，提取花椰菜小孢子不同发育时期的总RNA，经DNaseI处理去除基因组DNA污染后，利用随机引物或Oligo(dT)引物进行逆转录反应，将mRNA逆转录为cDNA。以逆转录得到的cDNA为模板，设计特异性的靶基因引物进行荧光定量PCR扩增。荧光定量PCR扩增体系和程序与miRNA检测类似，只是引物替换为靶基因特异性引物。通过实时监测荧光信号，利用软件计算出靶基因的相对表达量。同样采用2^(-ΔΔCt)法进行计算，内参基因可选择β-actin等表达稳定的基因，通过与内参基因的比较，得到靶基因在不同发育时期的相对表达量，进而分析miRNA与靶基因表达量之间的相关性，为深入研究miRNA对靶基因的调控机制提供实验依据。2.5.3miRNA靶基因5’RACE验证为了验证miRNA对靶基因的切割位点，采用5’RACE（RapidAmplificationofcDNAEnds）技术。5’RACE技术的原理基于mRNA的结构特点，利用基因特异性引物（GSP）和通用引物，通过逆转录和PCR扩增，从总RNA中扩增出包含miRNA切割位点的靶基因cDNA片段。首先，提取花椰菜小孢子总RNA，在逆转录酶（如SMARTerRACE5'/3'Kit中的逆转录酶）的作用下，利用基因特异性引物（GSP1）进行逆转录反应。GSP1与靶基因mRNA的3'端特异性结合，在逆转录酶的作用下，以mRNA为模板合成第一链cDNA。逆转录反应体系包含总RNA模板、GSP1引物、逆转录酶、逆转录缓冲液、dNTPs以及适量的RNase抑制剂，在42℃条件下反应60-90min。在第一链cDNA的3'末端添加poly(A)尾，利用末端脱氧核苷酸转移酶（TdT）和dATP进行加尾反应。加尾反应体系包含第一链cDNA、TdT、dATP以及反应缓冲液，在37℃条件下反应15-30min。以加尾后的第一链cDNA为模板，使用巢式PCR进行扩增。第一轮PCR使用通用引物（UPM）和基因特异性引物（GSP2），UPM与poly(A)尾互补配对，GSP2与靶基因mRNA上靠近miRNA切割位点的区域特异性结合。PCR扩增体系包含加尾后的第一链cDNA模板、引物、TaqDNA聚合酶、PCR缓冲液以及dNTPs。设置合适的PCR扩增程序，95℃预变性3-5min；95℃变性30-45s，60℃退火30-45s，72℃延伸30-45s，共进行30-35个循环；最后72℃延伸5-10min。将第一轮PCR产物稀释后作为模板，进行第二轮PCR扩增。第二轮PCR使用巢式通用引物（NUP）和巢式基因特异性引物（NGSP）。NUP与UPM内部序列互补，NGSP与靶基因mRNA上更靠近miRNA切割位点的区域特异性结合。PCR扩增体系和程序与第一轮PCR类似。扩增结束后，利用琼脂糖凝胶电泳检测PCR产物，观察是否出现预期大小的特异性条带。若出现特异性条带，则将其切胶回收，送往专业测序公司进行测序。测序公司对PCR产物的序列进行测定，将测序结果与靶基因的mRNA序列进行比对。通过分析比对结果，确定miRNA对靶基因的切割位点，从而验证miRNA与靶基因之间的相互作用关系，为深入揭示miRNA在花椰菜小孢子发育过程中的调控机制提供关键实验证据。三、花椰菜小孢子发育不同时期miRNA分析3.1小RNA测序数据结果3.1.1测序数据产量与质量评估利用IlluminaHiSeq2500测序平台对花椰菜小孢子发育的四分体时期、单核期和双核期样本进行小RNA高通量测序，共获得原始数据量总计为125.65Gb。经过严格的数据预处理，去除低质量reads、接头序列和污染序列后，得到高质量的cleanreads，有效数据量达118.92Gb，数据有效率为94.64%。通过FastQC软件对数据质量进行评估，结果显示各样本碱基质量分布良好，Q30（碱基错误率为0.1%）以上的碱基比例均在95%以上，表明测序数据质量高，准确性可靠，能够满足后续深入分析的要求。各样本的GC含量在48%-52%之间，处于合理范围，进一步验证了测序数据的可靠性。高质量的数据为全面、准确地鉴定花椰菜小孢子发育不同时期的miRNA提供了坚实的数据基础，确保了后续生物信息学分析结果的可信度。3.1.2小RNA长度分布特征对预处理后的cleanreads进行长度分布分析，结果表明花椰菜小孢子发育不同时期的小RNA长度主要集中在18-25nt之间。其中，21nt和24nt长度的小RNA占比最高，是主要的长度类型。在四分体时期，21nt长度的小RNA占比为35.68%，24nt长度的小RNA占比为28.45%；单核期21nt长度的小RNA占比为34.87%，24nt长度的小RNA占比为29.12%；双核期21nt长度的小RNA占比为36.05%，24nt长度的小RNA占比为27.98%。这与植物中miRNA的典型长度范围相符，21nt长度的小RNA在植物中通常参与RNA介导的转录后基因沉默途径，通过与靶mRNA的互补配对，引导RISC对靶mRNA进行切割或抑制翻译；24nt长度的小RNA则主要参与DNA甲基化等表观遗传调控过程，在植物基因组稳定性维持和基因表达调控中发挥重要作用。除了21nt和24nt长度的小RNA外，18-20nt和22-23nt长度的小RNA也有一定比例分布，它们可能在花椰菜小孢子发育过程中发挥着尚未被揭示的生物学功能。3.2已知miRNA分析3.2.1已知miRNA的注释结果经过生物信息学分析，将预处理后的cleanreads与miRBase数据库进行比对，在花椰菜小孢子发育的四分体时期、单核期和双核期样本中，共注释得到112个已知的miRNA，这些miRNA隶属于37个不同的家族。其中，miR156家族包含6个成员，是注释到成员数量较多的家族之一；miR166家族有5个成员，miR171家族也含有5个成员。不同家族的miRNA在小孢子发育各时期的分布存在差异。在四分体时期，检测到85个已知miRNA，其中miR156家族的成员表达相对较高，如miR156a的表达量为1256.34RPM，在该时期的已知miRNA表达量排名中位居前列。单核期共检测到92个已知miRNA，miR167家族的miR167a表达量显著升高，达到1895.67RPM，成为该时期高表达的miRNA之一。在双核期，有88个已知miRNA被检测到，miR396家族的miR396b表达活跃，表达量为1648.23RPM。这些已知miRNA在不同发育时期的分布和表达差异，暗示它们在花椰菜小孢子发育过程中可能发挥着不同的调控作用。3.2.2差异表达谱分析对不同发育时期已知miRNA的表达数据进行差异表达分析，结果显示，84个已知miRNA在小孢子不同发育阶段呈现差异表达。以四分体时期为对照，在单核期，有35个miRNA表达上调，40个miRNA表达下调。其中，miR160a的表达上调最为显著，其表达倍数变化（foldchange）达到3.25，P值为0.012，经过Benjamini-Hochberg校正后的FDR值为0.025，均小于设定的阈值，表明miR160a在单核期的表达差异具有统计学意义。而miR169g在单核期表达下调明显，foldchange为0.32，P值为0.008，FDR值为0.018，其表达水平在单核期显著低于四分体时期。在双核期，与四分体时期相比，有38个miRNA表达上调，36个miRNA表达下调。miR398b在双核期表达上调，foldchange为2.87，P值为0.015，FDR值为0.030，表明其在双核期的表达显著增加；miR408则表达下调，foldchange为0.28，P值为0.006，FDR值为0.015，在双核期的表达水平显著降低。进一步筛选出22个在小孢子发育过程中表现为显著差异表达的miRNA，这些miRNA可能在小孢子发育进程中发挥关键调控作用。对部分显著差异表达的miRNA，如miR165a、miR172a、miR319c、miR391和miR858进行深入分析。在小孢子从四分体时期向单核期发育过程中，miR165a的表达量逐渐升高，foldchange为2.15，P值为0.020，FDR值为0.035，到双核期时，其表达量略有下降，但仍高于四分体时期。miR172a的表达趋势则相反，从四分体时期到单核期表达量逐渐降低，foldchange为0.45，P值为0.010，FDR值为0.020，在双核期继续下降。这些显著差异表达的miRNA，可能通过调控下游靶基因的表达，参与花椰菜小孢子发育过程中的细胞分化、基因表达调控等生物学过程，对小孢子的正常发育起着至关重要的作用。3.3未知miRNA分析3.3.1未知miRNA的预测结果利用miRDeep2软件对花椰菜小孢子发育不同时期的测序数据进行分析，首次在花椰菜中预测到55个novelmiRNA。这些新预测的miRNA序列长度在20-24nt之间，其中长度为21nt的novelmiRNA有23个，占比41.82%；长度为22nt的有15个，占比27.27%；长度为23nt的有12个，占比21.82%；长度为24nt的有5个，占比9.09%。从序列特征来看，5'端碱基为尿嘧啶（U）的novelmiRNA有48个，占比87.27%，这与植物miRNA5'端对U的偏好性相符。对这些novelmiRNA的前体序列进行分析，发现它们均能形成典型的茎环结构，茎环结构的稳定性是判断其是否为潜在miRNA的重要依据。通过Mfold软件预测，这些前体序列形成的茎环结构最小自由能（MFE）范围在-30.5kcal/mol至-55.2kcal/mol之间，平均MFE为-42.8kcal/mol，表明其茎环结构较为稳定。以novel-miRNA-07为例，其前体序列长度为118nt，能够折叠形成稳定的茎环结构，成熟的novel-miRNA-07位于茎环结构的一条臂上，且与另一条臂上的互补序列形成完美的碱基配对，进一步验证了其作为miRNA的可能性。这些新预测的novelmiRNA为深入研究花椰菜小孢子发育过程中的基因调控网络提供了新的线索。3.3.2差异表达分析对预测得到的55个novelmiRNA在花椰菜小孢子发育的四分体时期、单核期和双核期的表达数据进行差异表达分析，结果显示，有15个novelmiRNA在小孢子不同发育阶段呈现差异表达。以四分体时期为对照，在单核期，有6个novelmiRNA表达上调，5个novelmiRNA表达下调。其中，novel-miRNA-12表达上调最为显著，foldchange达到2.68，P值为0.023，经过Benjamini-Hochberg校正后的FDR值为0.038，表明其在单核期的表达差异具有统计学意义。novel-miRNA-23在单核期表达下调明显，foldchange为0.38，P值为0.015，FDR值为0.028，其表达水平在单核期显著低于四分体时期。在双核期，与四分体时期相比，有7个novelmiRNA表达上调，4个novelmiRNA表达下调。novel-miRNA-35在双核期表达上调，foldchange为2.45，P值为0.028，FDR值为0.042，表明其在双核期的表达显著增加；novel-miRNA-41则表达下调，foldchange为0.32，P值为0.010，FDR值为0.020，在双核期的表达水平显著降低。进一步对这15个差异表达的novelmiRNA进行实验验证，采用茎环PCR扩增和荧光实时定量PCR技术。茎环PCR扩增结果显示，15个novelmiRNA均能扩增出预期大小的特异性条带，表明这些novelmiRNA真实存在于花椰菜小孢子中。荧光实时定量PCR结果与生物信息学分析的差异表达结果高度一致，进一步证实了11个novelmiRNA的差异表达特征。这些差异表达的novelmiRNA可能在花椰菜小孢子发育过程中发挥着重要的调控作用，通过对其功能的深入研究，有望揭示花椰菜小孢子发育的新机制。3.4miRNA靶基因分析3.4.1靶基因预测结果在明确花椰菜小孢子发育不同时期的已知和未知miRNA后，利用psRNATarget软件对这些miRNA的靶基因进行预测。结果显示，针对112个已知miRNA，共预测到756个靶基因；对于新发现的55个novelmiRNA，预测得到321个靶基因。在已知miRNA的靶基因中，miR156家族的靶基因主要为SPL（SQUAMOSAPROMOTERBINDINGPROTEIN-LIKE）基因家族成员，如SPL3、SPL9和SPL15等。以miR156a为例，它能够与SPL9基因的mRNA序列互补配对，预测结果显示二者的互补区域主要位于SPL9基因mRNA的编码区，且错配率较低，结合位点高度保守。在novelmiRNA的靶基因中，novel-miRNA-07预测到的靶基因有BoTCP14和BoMYB33等。其中，novel-miRNA-07与BoTCP14基因mRNA的互补区域位于3'UTR端，存在少量错配，但不影响二者的结合。通过对这些靶基因的功能初步分析，发现它们涉及植物生长发育、激素信号转导、代谢调控等多个生物学过程。将所有预测得到的靶基因进行整理，形成详细的靶基因列表，为后续深入研究miRNA在花椰菜小孢子发育中的调控机制提供基础数据。3.4.2靶基因功能注释对预测得到的靶基因进行功能注释，将其序列与GO数据库和KEGG数据库进行比对。在GO功能分类中，从生物过程层面来看，靶基因主要富集在细胞过程、代谢过程、发育过程等功能条目。其中，参与细胞过程的靶基因占比为32.6%，涉及细胞周期调控、细胞分化、细胞分裂等过程；参与代谢过程的靶基因占比为28.4%，涵盖碳水化合物代谢、蛋白质代谢、脂质代谢等多种代谢途径；参与发育过程的靶基因占比为15.8%，包括花器官发育、胚胎发育、细胞组织分化等。从细胞组成层面分析，靶基因主要富集在细胞、细胞器、细胞膜等功能条目。在细胞组成中，定位在细胞内的靶基因占比为45.2%，分布于细胞核、细胞质、线粒体等不同细胞器中的靶基因占比为30.5%，与细胞膜相关的靶基因占比为12.3%，这些靶基因在维持细胞结构和功能完整性方面发挥重要作用。在分子功能层面，靶基因主要富集在催化活性、结合活性、转录调控活性等功能条目。具有催化活性的靶基因占比为30.1%，参与各种酶促反应；具有结合活性的靶基因占比为27.8%，能够与核酸、蛋白质、小分子等物质结合；具有转录调控活性的靶基因占比为10.6%，通过调控基因转录影响花椰菜小孢子发育过程中的基因表达。在KEGG代谢通路分析中，靶基因主要富集在植物激素信号转导、淀粉和蔗糖代谢、苯丙氨酸代谢等通路。在植物激素信号转导通路中，涉及生长素、细胞分裂素、赤霉素、脱落酸等多种激素信号转导途径的相关基因，如生长素信号转导通路中的TIR1、ARF（AUXINRESPONSEFACTOR）基因，细胞分裂素信号转导通路中的CRE1（CYTOKININRESPONSE1）基因等，表明miRNA可能通过调控这些基因参与植物激素信号网络，进而影响花椰菜小孢子发育。在淀粉和蔗糖代谢通路中，多个关键酶基因如蔗糖合成酶（SUS）、淀粉合成酶（SS）等被预测为靶基因，暗示miRNA在调控小孢子发育过程中的能量代谢和物质合成方面发挥作用。通过对靶基因的GO功能分类和KEGG代谢通路富集分析，系统地揭示了miRNA在花椰菜小孢子发育过程中的潜在调控作用和生物学功能。三、花椰菜小孢子发育不同时期miRNA分析3.5实验验证结果3.5.1茎环PCR及重测序验证结果为验证生物信息学分析所鉴定出的miRNA的真实性，对部分已知和未知miRNA进行茎环PCR扩增。以花椰菜小孢子不同发育时期的总RNA为模板，利用特异性茎环引物进行逆转录和PCR扩增。扩增产物经1.5%琼脂糖凝胶电泳检测，结果显示（图1），在预期大小位置出现清晰明亮的特异性条带。对于已知miRNA，如miR167a，在四分体时期、单核期和双核期样本中均扩增出约150bp的条带，与理论片段大小相符；对于新预测的novel-miRNA-07，也在各时期样本中扩增出约140bp的特异性条带。这表明通过茎环PCR成功扩增出目标miRNA，初步验证了这些miRNA在花椰菜小孢子中的存在。将茎环PCR扩增产物进行重测序，测序结果与生物信息学预测的miRNA序列进行比对。结果显示，已知miRNA的测序序列与miRBase数据库中的对应序列一致性高达99%以上，仅有个别碱基由于测序误差存在差异，但不影响miRNA的序列特征和功能预测。对于novel-miRNA-07等新预测的miRNA，测序序列与miRDeep2软件预测的序列一致性达到98%以上，进一步证实了新miRNA预测的准确性。通过茎环PCR及重测序验证，为后续深入研究花椰菜小孢子发育过程中miRNA的功能和调控机制提供了可靠的实验依据。[此处插入茎环PCR扩增产物的电泳图，图注：M为DNAMarker；1-3分别为四分体时期、单核期和双核期样本中miR167a的扩增条带；4-6分别为四分体时期、单核期和双核期样本中novel-miRNA-07的扩增条带]3.5.2荧光实时定量PCR验证结果为进一步验证测序数据中miRNA和靶基因表达量的准确性，采用荧光实时定量PCR技术对部分差异表达的miRNA及其靶基因进行检测。选取miR165a、miR172a等已知差异表达miRNA，以及novel-miRNA-12、novel-miRNA-23等新预测的差异表达miRNA进行分析。以U6snRNA作为内参基因，计算miRNA的相对表达量。荧光实时定量PCR结果显示（图2），miR165a在单核期的表达量相对于四分体时期显著上调，表达倍数为3.12，与测序数据中foldchange为3.08的结果高度一致；miR172a在单核期表达量显著下调，表达倍数为0.42，与测序数据中foldchange为0.45的结果相符。对于novel-miRNA-12，在单核期表达上调，表达倍数为2.75，与测序分析中foldchange为2.68的结果相近；novel-miRNA-23在单核期表达下调，表达倍数为0.35，与测序数据中的foldchange为0.38相符。同时，对miRNA的靶基因进行荧光实时定量PCR检测。以miR165a的靶基因PHB（PHABULOSA）为例，其在单核期的表达量相对于四分体时期显著下调，表达倍数为0.38，与miR165a的表达趋势呈负相关，符合miRNA对靶基因的负调控规律。novel-miRNA-12的靶基因BoTCP14在单核期表达下调，表达倍数为0.45，与novel-miRNA-12的表达变化趋势相反。通过荧光实时定量PCR验证，证实了测序数据中miRNA和靶基因表达量变化的可靠性，为深入研究miRNA在花椰菜小孢子发育过程中的调控机制提供了有力的实验证据。[此处插入荧光实时定量PCR验证miRNA和靶基因表达量的柱状图，图注：*表示P<0.05，差异具有统计学意义；**表示P<0.01，差异极显著]3.5.3miRNA靶基因5’RACE验证结果为验证miRNA对靶基因的切割位点，采用5’RACE技术对部分miRNA-靶基因对进行分析。选取miR164-NAC1和novel-miRNA-07-BoTCP14这两对miRNA-靶基因进行5’RACE验证。以花椰菜小孢子总RNA为模板，通过逆转录和巢式PCR扩增，获得包含miRNA切割位点的靶基因cDNA片段。扩增产物经1.5%琼脂糖凝胶电泳检测，结果显示（图3），在预期大小位置出现特异性条带。对miR164-NAC1，扩增出约350bp的条带；对于novel-miRNA-07-BoTCP14，扩增出约400bp的条带。将特异性条带切胶回收后进行测序。测序结果与靶基因mRNA序列进行比对，结果表明，miR164在NAC1基因mRNA的第125-145位核苷酸处进行切割，切割位点与生物信息学预测结果一致；novel-miRNA-07在BoTCP14基因mRNA的第210-230位核苷酸处进行切割，验证了miRNA对靶基因的切割作用。通过5’RACE验证，明确了miRNA在花椰菜小孢子发育过程中对靶基因的切割位点，为深入揭示miRNA的调控机制提供了关键实验证据。[此处插入5’RACE验证miRNA靶基因切割位点的电泳图，图注：M为DNAMarker；1为miR164-NAC1扩增条带；2为novel-miRNA-07-BoTCP14扩增条带]四、花椰菜miRNA功能初步研究4.1候选miRNA载体构建4.1.1候选miRNA前体的获得通过查阅前期花椰菜小孢子发育不同时期miRNA分析的结果，筛选出在小孢子发育进程中表达差异显著且功能未知的miRNA作为候选miRNA，确定了novel-miRNA-12、novel-miRNA-23以及已知的miR165a、miR172a等为重点研究对象。针对这些候选miRNA，利用NCBI数据库和花椰菜基因组信息，获取其前体序列。若前体序列已知，根据序列信息设计特异性引物进行PCR扩增。引物设计遵循引物长度在18-25bp之间，GC含量在40%-60%，引物3'端避免出现连续的3个以上相同碱基等原则。以花椰菜小孢子不同发育时期的cDNA为模板，进行PCR扩增。反应体系包含cDNA模板2μL、上下游引物（10μM）各1μL、2×TaqPCRMasterMix12.5μL、ddH₂O8.5μL，总体积25μL。PCR扩增程序为95℃预变性3min；95℃变性30s，55-60℃退火30s，72℃延伸30-60s（根据前体序列长度调整），共进行30-35个循环；最后72℃延伸5min。扩增结束后，利用1.5%琼脂糖凝胶电泳检测PCR产物，观察是否出现预期大小的特异性条带。若前体序列未知或难以通过PCR扩增获得，采用化学合成的方法，委托专业的生物公司合成相应的miRNA前体序列。合成的前体序列两端带有特定的酶切位点，以便后续进行载体构建。4.1.2植物过表达载体的构建选用植物表达载体pCAMBIA1300作为基础载体，该载体含有CaMV35S启动子，能够驱动外源基因在植物中高效表达，还带有潮霉素抗性基因，便于后续转化植株的筛选。首先，用限制性内切酶BamHI和SacI对pCAMBIA1300载体进行双酶切。反应体系包含pCAMBIA1300载体5μg、10×Buffer5μL、BamHI和SacI各2μL（10U/μL），加ddH₂O补足至50μL。37℃水浴酶切3-4h。酶切结束后，利用1%琼脂糖凝胶电泳分离酶切产物，切胶回收线性化的载体片段。将通过PCR扩增或化学合成获得的候选miRNA前体序列与线性化的pCAMBIA1300载体进行连接。连接反应采用T4DNA连接酶，反应体系包含线性化载体片段50ng、miRNA前体序列100-200ng、10×T4DNALigaseBuffer2μL、T4DNALigase1μL（3U/μL），加ddH₂O补足至20μL。16℃连接过夜。连接产物转化大肠杆菌DH5α感受态细胞。将连接产物加入到100μLDH5α感受态细胞中，轻轻混匀，冰浴30min；42℃热激90s，迅速放回冰浴中冷却2-3min；加入800μL无抗生素的LB液体培养基，37℃振荡培养1h，使细胞复苏。取100-200μL菌液涂布于含有50μg/mL潮霉素的LB固体培养基平板上，37℃倒置培养12-16h。挑取平板上的单菌落，接种于含有50μg/mL潮霉素的LB液体培养基中，37℃振荡培养过夜。采用碱裂解法提取质粒DNA，具体步骤为：取1.5mL菌液于离心管中，12000rpm离心1min，弃上清；加入100μL冰预冷的SolutionI（含25mMTris-HCl，pH8.0，10mMEDTA，50mM葡萄糖），涡旋振荡使菌体完全悬浮；加入200μL新配制的SolutionII（含0.2MNaOH，1%SDS），轻轻颠倒混匀5-6次，室温放置3-5min；加入150μL冰预冷的SolutionIII（含3M醋酸钾，pH4.8），立即轻轻颠倒混匀5-6次，冰浴10min；12000rpm离心10min，将上清转移至新的离心管中；加入等体积的酚：氯仿：异戊醇（25:24:1），振荡混匀，12000rpm离心10min；将上清转移至新的离心管中，加入2倍体积的无水乙醇，-20℃放置20-30min；12000rpm离心10min，弃上清，用75%乙醇洗涤沉淀2次；晾干沉淀后，加入30-50μLddH₂O溶解质粒DNA。利用PCR和酶切鉴定重组质粒。以提取的质粒DNA为模板，使用miRNA前体特异性引物进行PCR扩增，验证插入片段的大小；同时用BamHI和SacI对质粒进行双酶切鉴定，酶切产物经1%琼脂糖凝胶电泳检测，若出现预