探寻红梨果皮着色密码：DNA 甲基化与 microRNA 的调控机制

探寻红梨果皮着色密码：DNA甲基化与microRNA的调控机制一、绪论1.1研究背景与意义梨（PyrusL.）作为世界范围内广泛种植的重要温带落叶果树，在全球水果产业中占据重要地位。中国作为梨的起源中心之一，拥有悠久的栽培历史和丰富的种质资源，是世界最大的梨生产国和消费国。据统计，2023年中国梨的种植面积达106.8万公顷，产量达到1887.4万吨，分别占世界总量的69.5%和71.6%。梨的品种丰富多样，根据果皮颜色可大致分为绿皮梨、褐皮梨和红皮梨。红梨因富含花青苷呈现出鲜艳的红色果皮，不仅具有较高的观赏价值，还因其富含抗氧化物质，对人体健康具有诸多益处，如预防心血管疾病、抗氧化、抗炎等，备受消费者青睐，在市场上具有较高的经济价值。近年来，随着消费者对水果品质和外观要求的不断提高，红梨的市场需求逐渐增加，推动了红梨产业的快速发展。在中国，云南、四川、陕西等地成为红梨的主要产区，种植面积和产量逐年递增。以云南为例，2023年红梨种植面积达到12.5万公顷，产量达到156.3万吨，分别较上一年增长了8.2%和10.5%。新品种如‘红酥脆’‘满天红’‘早酥红’等在市场上崭露头角，深受消费者喜爱。同时，红梨的出口量也在逐年增加，2023年中国红梨出口量达到12.6万吨，出口额达到1.8亿美元，主要出口到东南亚、欧洲和北美等地区，进一步提升了中国红梨在国际市场上的影响力。然而，红梨产业在发展过程中也面临着一些挑战。其中，果皮着色不稳定是制约红梨品质和市场竞争力的关键因素之一。红梨果皮的着色程度受多种因素的影响，包括遗传因素、环境因素（如光照、温度、水分等）以及栽培管理措施（如施肥、修剪、套袋等）。在实际生产中，常出现红梨果实着色不均匀、色泽暗淡等问题，导致果实外观品质下降，商品价值降低。例如，在一些年份或地区，由于光照不足或温度异常，部分红梨果实的着色率不足50%，严重影响了果农的经济效益。花青苷作为红梨果皮呈现红色的主要色素，其生物合成和积累过程受到复杂的调控机制的影响。深入研究红梨果皮着色的调控机制，对于解决红梨着色不稳定问题，提高果实品质和市场竞争力具有重要的理论和实践意义。近年来，随着分子生物学技术的快速发展，关于果实色泽形成机制的研究取得了显著进展。研究表明，DNA甲基化和microRNA作为重要的表观遗传调控因子，在植物生长发育、逆境响应以及次生代谢产物合成等过程中发挥着关键作用。在果实色泽调控方面，DNA甲基化可以通过影响花青苷合成相关基因的表达，进而影响果实的着色过程。例如，在苹果中，MdMYB1基因启动子区域的甲基化水平与花青苷的积累呈负相关，低甲基化水平有利于MdMYB1基因的表达，从而促进花青苷的合成和果实着色。microRNA则可以通过靶向调控花青苷合成相关基因的mRNA稳定性或翻译过程，实现对花青苷合成的调控。在葡萄中，miR156通过靶向调控SBP转录因子，间接影响花青苷合成相关基因的表达，从而调控葡萄果实的着色。然而，目前关于DNA甲基化和microRNA在红梨果皮着色调控中的作用机制研究仍相对较少，相关的调控网络尚未完全解析。因此，开展DNA甲基化和microRNA调控红梨果皮着色机制的研究，不仅可以丰富果实色泽形成的理论体系，还可以为红梨的品种改良和栽培管理提供科学依据和技术支持，具有重要的理论意义和实践价值。1.2植物花青苷概述花青苷（Anthocyanin）是一类广泛存在于植物界的水溶性天然色素，属于黄酮多酚类化合物，在植物的生长发育、防御机制以及对环境的适应过程中发挥着关键作用。其基本结构是由一个花色素（Anthocyanidin）与一个或多个糖分子通过糖苷键连接而成。花色素的核心结构为2-苯基苯并吡喃阳离子，具有高度共轭的π电子体系，这一结构特征赋予了花青苷独特的光学性质，使其能够吸收特定波长的可见光，从而呈现出从红色、紫色到蓝色等丰富多样的颜色。根据花色素B环上羟基化和甲基化的位置及数目不同，常见的花色素主要有六种，分别是天竺葵色素（Pelargonidin）、矢车菊色素（Cyanidin）、飞燕草色素（Delphinidin）、芍药色素（Peonidin）、矮牵牛色素（Petunidin）和锦葵色素（Malvidin）。这些不同的花色素通过与不同种类、数量和连接位置的糖分子结合，以及可能的酰基化修饰，形成了结构复杂、种类繁多的花青苷。例如，矢车菊色素3-O-葡萄糖苷是最为常见的花青苷之一，广泛存在于多种植物的花朵和果实中；而在葡萄中，锦葵色素3-O-葡萄糖苷则是重要的呈色物质，对葡萄果实的色泽形成起到关键作用。花青苷在植物体内分布广泛，几乎存在于植物的各个器官，包括花、果实、叶片、茎和种子等。在花朵中，花青苷赋予了花朵鲜艳的色彩，吸引昆虫传粉，促进植物的繁殖。如玫瑰的红色花瓣、郁金香的紫色花瓣等，都是由于花青苷的存在而呈现出绚丽的颜色。在果实中，花青苷不仅影响果实的外观色泽，还与果实的品质和成熟度密切相关。随着果实的成熟，花青苷的合成逐渐增加，使果实的颜色由绿色转变为红色、紫色等，同时也提升了果实的抗氧化能力和营养价值。像草莓、蓝莓、樱桃等水果，其鲜艳的色泽主要归功于花青苷的积累。在叶片中，花青苷的合成通常与植物的生长发育阶段、环境胁迫等因素有关。在一些植物的幼叶或衰老叶片中，花青苷的含量会相对较高，可能起到保护叶片免受光氧化伤害、抵御病虫害等作用。此外，在某些植物的茎和种子中，花青苷也有一定的分布，对植物的生理功能和生态适应性具有重要意义。1.3花青苷合成途径及其影响因子1.3.1花青苷生物合成的结构基因花青苷的生物合成是一个复杂的代谢过程，涉及多个结构基因的协同作用，这些基因编码的酶参与了花青苷合成的各个步骤。其生物合成起始于苯丙氨酸，通过苯丙烷途径，在苯丙氨酸解氨酶（PAL）、肉桂酸-4-羟化酶（C4H）和4-香豆酰辅酶A连接酶（4CL）的依次催化下，生成4-香豆酰辅酶A。4-香豆酰辅酶A与丙二酰辅酶A在查尔酮合酶（CHS）的作用下，合成柚皮素查尔酮，这是花青苷合成途径中的第一个关键中间产物。查尔酮异构酶（CHI）可将柚皮素查尔酮转化为柚皮素，为后续反应提供底物。在后续的反应中，柚皮素在黄烷酮3-羟化酶（F3H）的催化下，羟基化形成二氢山奈酚。二氢山奈酚可在类黄酮3'-羟化酶（F3'H）或类黄酮3',5'-羟化酶（F3'5'H）的作用下，进一步羟基化，分别生成二氢槲皮素和二氢杨梅素。这些二氢黄酮醇在二氢黄酮醇4-还原酶（DFR）的催化下，被还原为无色花色素，如无色天竺葵色素、无色矢车菊色素和无色飞燕草色素。无色花色素在花青素合成酶（ANS）的作用下，氧化生成花色素，如花青素、矢车菊素和飞燕草素。最后，花色素在UDP-葡萄糖：类黄酮3-O-葡萄糖基转移酶（UFGT）的催化下，与UDP-葡萄糖结合，形成稳定的花青苷。在梨中，研究发现PpCHS、PpCHI、PpF3H、PpF3'H、PpDFR、PpANS和PpUFGT等结构基因在花青苷合成过程中发挥着重要作用。这些基因的表达水平与梨果皮中花青苷的含量密切相关，其表达上调通常会促进花青苷的合成和积累，从而使果皮颜色加深。如在红皮梨品种‘满天红’中，随着果实的发育和着色，PpDFR、PpANS和PpUFGT等基因的表达量显著增加，花青苷含量也随之上升，果皮逐渐呈现出红色。1.3.2花青苷生物合成途径相关转录因子转录因子在花青苷生物合成途径中起着关键的调控作用，它们通过与结构基因启动子区域的顺式作用元件结合，激活或抑制结构基因的转录，从而影响花青苷的合成。在植物中，MYB、bHLH和WD40转录因子家族成员形成的MBW复合体是调控花青苷生物合成的核心转录调控模块。MYB转录因子是一类含有保守MYB结构域的转录因子，在花青苷合成调控中发挥着重要作用。在梨中，PpMYB10被证明是调控花青苷合成的关键MYB转录因子。它能够直接结合到花青苷合成相关结构基因（如PpDFR、PpANS和PpUFGT）的启动子区域，激活这些基因的表达，从而促进花青苷的合成。研究表明，在‘红早酥’梨中，过表达PpMYB10基因能够显著提高花青苷合成相关结构基因的表达水平，使果实中花青苷含量增加，果皮颜色加深。bHLH转录因子含有保守的碱性螺旋-环-螺旋结构域，可与MYB转录因子相互作用，共同调控花青苷的合成。在苹果中，MdMYC1是一个重要的bHLH转录因子，它与MdMYB1相互作用形成复合体，增强MdMYB1对花青苷合成相关基因的激活作用，促进花青苷的积累。在梨中，也存在类似的bHLH转录因子，它们与PpMYB10相互作用，协同调控花青苷的合成。WD40蛋白是一类含有多个WD40重复序列的蛋白质，在MBW复合体中起支架作用，促进MYB和bHLH转录因子之间的相互作用。在矮牵牛中，AN11是一个WD40蛋白，它与MYB转录因子AN2和bHLH转录因子JAF13相互作用，形成MBW复合体，调控花青苷的合成。在梨中，虽然对WD40蛋白的研究相对较少，但推测也存在类似的WD40蛋白参与花青苷合成的调控。除了MBW复合体，其他转录因子家族成员如WRKY、NAC、ERF等也参与了花青苷合成的调控。WRKY转录因子可以通过与花青苷合成相关基因的启动子区域结合，直接调控基因的表达，或者通过与MBW复合体相互作用，间接影响花青苷的合成。NAC转录因子在植物的生长发育和逆境响应中发挥重要作用，也参与了花青苷合成的调控。在苹果中，MdNAC47能够直接结合到MdMYB1的启动子区域，抑制其表达，从而负调控花青苷的合成。ERF转录因子是乙烯响应元件结合蛋白，在植物对乙烯的响应和花青苷合成调控中具有重要作用。在海棠中，MpERF105能够通过与花青苷合成调控基因MpMYB10b的启动子结合并激活其表达，间接促进花青苷合成。1.3.3光对花青苷合成的影响光作为植物生长发育过程中重要的环境信号，对花青苷的合成起着至关重要的调控作用。光照强度、光质和光照时间等因素都会影响花青苷的合成和积累。光照强度是影响花青苷合成的重要因素之一。适度的光照强度可以促进花青苷的合成，而光照不足或过强则会抑制花青苷的合成。在葡萄中，适度遮光会导致果实中花青苷含量下降，而增加光照强度则会促进花青苷的积累，使果实颜色更加鲜艳。这是因为光照强度影响了花青苷合成相关基因的表达。适度光照可以激活光信号转导途径，促进转录因子（如HY5、BBX等）的表达，这些转录因子能够与花青苷合成相关基因的启动子区域结合，激活基因的表达，从而促进花青苷的合成。然而，过强的光照会导致植物产生过多的活性氧，对细胞造成损伤，从而抑制花青苷的合成。光质对花青苷合成也具有显著影响。不同波长的光对花青苷合成的调控作用不同。红光和蓝光是植物生长发育过程中最重要的两种光质，它们对花青苷合成的影响尤为显著。红光可以促进花青苷的合成，在苹果中，红光处理能够显著提高果实中花青苷的含量，这是因为红光可以激活光信号转导途径，促进花青苷合成相关转录因子的表达，进而促进花青苷合成相关基因的表达。蓝光也能促进花青苷的合成，在拟南芥中，蓝光可以通过激活CRY1和CRY2等蓝光受体，促进HY5等转录因子的表达，从而促进花青苷的合成。此外，远红光、紫外光等光质也对花青苷合成有一定的影响。远红光可以通过与红光的相互作用，调节花青苷的合成；紫外光可以诱导植物产生防御反应，促进花青苷的合成，增强植物对紫外光的耐受性。光照时间也会影响花青苷的合成。长日照条件下，植物花青苷的合成通常会受到促进。在草莓中，长日照处理可以增加果实中花青苷的含量，使果实颜色更加鲜艳。这是因为长日照条件下，光信号转导途径持续激活，促进了花青苷合成相关基因的表达，从而有利于花青苷的合成和积累。相反，短日照条件下，花青苷的合成可能会受到抑制。1.4DNA甲基化对花青苷合成的调控1.4.1DNA甲基化的形成DNA甲基化是一种重要的表观遗传修饰，在植物的生长发育、基因表达调控以及对环境胁迫的响应等过程中发挥着关键作用。在植物中，DNA甲基化主要发生在胞嘧啶残基上，形成5-甲基胞嘧啶（5mC），根据其序列背景可分为CG、CHG和CHH（H代表A、T或C）三种类型。DNA甲基化的形成是一个复杂的过程，涉及多种酶的参与。其中，DNA甲基转移酶（DNAmethyltransferases，DMTs）是催化DNA甲基化反应的关键酶。在植物中，主要存在三种类型的DNA甲基转移酶，分别是MET1（METHYLTRANSFERASE1）、CMT3（CHROMOMETHYLASE3）和DRM2（DOMAINSREARRANGEDMETHYLTRANSFERASE2）。MET1主要负责维持CG位点的甲基化，它与哺乳动物中的DNMT1具有较高的同源性，通过识别半甲基化的DNA双链，将甲基基团添加到新合成链的CG位点上，从而保持CG甲基化的稳定性。例如，在拟南芥中，MET1基因的突变会导致CG位点甲基化水平显著降低，进而影响植物的生长发育。CMT3主要负责维持CHG位点的甲基化，它含有一个独特的结构域，能够与组蛋白H3的赖氨酸9甲基化（H3K9me）相互作用，从而特异性地识别并甲基化CHG位点。在拟南芥中，cmt3突变体的CHG位点甲基化水平明显下降，导致一些基因的表达发生改变，影响植物的形态建成和对逆境的响应。DRM2则参与从头甲基化（denovomethylation）过程，即在未甲基化的DNA位点上建立新的甲基化修饰。它可以在RNA介导的DNA甲基化（RNA-directedDNAmethylation，RdDM）途径中发挥作用，通过与小干扰RNA（smallinterferingRNAs，siRNAs）相互作用，识别并甲基化特定的DNA序列。例如，在RdDM途径中，双链RNA被核酸酶切割成24-nt的siRNAs，这些siRNAs与AGO4（ARGONAUTE4）蛋白结合形成复合物，然后引导DRM2到靶位点进行DNA甲基化。此外，DNA去甲基化酶在DNA甲基化动态调控中也起着重要作用。植物中的DNA去甲基化酶主要包括DME（DEMETER）、ROS1（REPRESSOROFSILENCING1）、DML2（DEMETER-LIKE2）和DML3（DEMETER-LIKE3）等。它们通过碱基切除修复机制，将甲基化的胞嘧啶去除，从而实现DNA去甲基化。在拟南芥中，ROS1能够识别并结合甲基化的DNA位点，将5mC转化为5-羟甲基胞嘧啶（5hmC），然后进一步通过一系列反应将其去除，完成DNA去甲基化过程。DNA甲基化和去甲基化的动态平衡，精确调控着植物基因的表达，影响植物的生长发育和对环境的适应能力。1.4.2DNA甲基化对植物生长发育及抗逆性的调控DNA甲基化作为一种重要的表观遗传修饰，在植物的整个生命周期中对其生长发育起着至关重要的调控作用。在种子萌发阶段，DNA甲基化状态影响着种子的休眠与萌发进程。研究表明，一些与种子休眠相关基因的启动子区域甲基化水平变化，会直接影响基因的表达，从而调控种子对萌发信号的响应。例如，在拟南芥中，DOG1基因启动子的低甲基化状态有利于其表达，维持种子的休眠；而在种子萌发时，该区域甲基化水平升高，DOG1表达受到抑制，种子得以顺利萌发。在植物营养生长阶段，DNA甲基化参与调控植物的株型、叶片发育、根系生长等多个方面。在水稻中，DNA甲基化影响分蘖角度相关基因的表达，进而调控水稻的株型。一些甲基化水平改变的水稻突变体表现出分蘖角度异常，影响了植株的空间分布和光合作用效率。在叶片发育方面，DNA甲基化参与调控叶片的形态建成和衰老进程。在烟草中，通过调控相关基因的甲基化水平，可以改变叶片的大小、形状和衰老时间。根系作为植物吸收水分和养分的重要器官，其生长和发育也受到DNA甲基化的精细调控。在拟南芥中，DNA甲基化影响根细胞的分裂和分化，突变体中异常的DNA甲基化导致根系形态改变，影响植物对环境资源的获取能力。在植物生殖生长阶段，DNA甲基化对花的发育、花粉育性、果实发育等过程至关重要。在拟南芥中，一些与花器官发育相关基因的甲基化状态改变会导致花器官形态异常，影响植物的繁殖能力。花粉育性也与DNA甲基化密切相关，异常的甲基化会导致花粉发育受阻，降低花粉的活力和受精能力。在果实发育过程中，DNA甲基化参与调控果实的大小、形状、色泽和成熟进程。在番茄中，DNA甲基化影响果实成熟相关基因的表达，调控果实的成熟时间和品质。除了对植物生长发育的调控，DNA甲基化在植物应对生物和非生物胁迫过程中也发挥着关键作用，增强植物的抗逆性。在面对干旱胁迫时，植物会通过改变DNA甲基化水平来调控相关基因的表达，从而提高自身的耐旱能力。在拟南芥中，干旱胁迫下一些与干旱响应相关基因的启动子区域甲基化水平发生变化，导致基因表达上调或下调，使植物能够通过调节气孔开闭、渗透调节物质积累等生理过程来适应干旱环境。在盐胁迫下，DNA甲基化同样参与调控植物的耐盐机制。在水稻中，盐胁迫诱导一些基因的甲基化水平改变，影响离子平衡调节、抗氧化防御等相关基因的表达，增强水稻对盐胁迫的耐受性。在生物胁迫方面，DNA甲基化在植物抵御病原菌入侵过程中发挥重要作用。在拟南芥与病原菌互作过程中，病原菌侵染会诱导植物基因组DNA甲基化水平和模式的改变，从而调控植物防御相关基因的表达，激活植物的免疫反应。一些研究表明，DNA甲基化可以通过调控植物激素信号通路相关基因的表达，间接影响植物对病原菌的抗性。例如，在烟草中，DNA甲基化参与调控水杨酸信号通路相关基因的表达，增强烟草对病毒的抗性。1.4.3DNA甲基化对花青苷合成的调控DNA甲基化在植物花青苷合成调控中扮演着关键角色，通过对花青苷合成相关基因的甲基化修饰，影响基因的表达水平，进而调控花青苷的合成与积累，最终决定植物器官的色泽。许多研究表明，DNA甲基化与花青苷合成相关基因的表达呈负相关关系。当基因启动子区域的甲基化水平较高时，转录因子难以结合到启动子上，从而抑制基因的转录，导致花青苷合成减少；反之，当甲基化水平降低时，基因表达被激活，促进花青苷的合成。在苹果中，MdMYB1基因是调控花青苷合成的关键转录因子。研究发现，MdMYB1基因启动子区域的甲基化水平与花青苷积累密切相关。在红色果皮的苹果品种中，MdMYB1基因启动子区域甲基化水平较低，该基因能够正常表达，激活下游花青苷合成相关结构基因（如MdDFR、MdANS、MdUFGT等）的表达，从而促进花青苷的合成，使果皮呈现红色；而在绿色果皮的苹果品种中，MdMYB1基因启动子区域甲基化水平较高，基因表达受到抑制，花青苷合成减少，果皮呈现绿色。通过对苹果果实进行甲基化抑制剂5-氮杂胞苷（5-azaC）处理，降低了MdMYB1基因启动子区域的甲基化水平，显著提高了该基因的表达量，同时花青苷合成相关结构基因的表达也上调，果实中花青苷含量增加，果皮颜色变红。在葡萄中，VvMYBA1和VvMYBA2基因是调控花青苷合成的重要转录因子。研究表明，这两个基因启动子区域的甲基化状态影响其表达水平和花青苷的积累。在紫色葡萄品种中，VvMYBA1和VvMYBA2基因启动子区域甲基化水平较低，基因表达活跃，促进花青苷合成，使果实呈现紫色；而在白色葡萄品种中，这两个基因启动子区域存在较高水平的甲基化，基因表达受到抑制，花青苷合成受阻，果实呈现白色。进一步研究发现，DNA甲基化通过影响转录因子与启动子区域顺式作用元件的结合能力，调控基因的表达。在葡萄中，甲基化修饰使得VvMYBA1和VvMYBA2基因启动子区域的顺式作用元件（如MYB结合位点）被屏蔽，转录因子无法有效结合，从而抑制基因的转录。在梨中，虽然相关研究相对较少，但已有研究表明DNA甲基化也参与了花青苷合成的调控。以红皮梨品种‘红早酥’和绿皮梨品种‘早酥’为材料，对花青苷合成相关基因的甲基化水平进行分析，发现‘红早酥’中PpMYB10基因启动子区域甲基化水平低于‘早酥’，且PpMYB10基因的表达量和花青苷含量在‘红早酥’中显著高于‘早酥’。这表明在梨中，DNA甲基化可能通过调控PpMYB10基因的表达，影响花青苷的合成和果皮颜色的形成。此外，对花青苷合成相关结构基因（如PpCHS、PpDFR、PpANS等）的甲基化分析发现，在红皮梨和绿皮梨中这些基因的甲基化水平也存在差异，进一步说明DNA甲基化在梨果皮花青苷合成调控中具有重要作用。1.5MicroRNA调控花青苷合成1.5.1MicroRNA的生物合成及作用机制MicroRNA（miRNA）是一类内生的、长度约为21-24个核苷酸的非编码单链RNA分子，在植物的生长发育、基因表达调控以及对环境胁迫的响应等过程中发挥着至关重要的作用。其生物合成过程是一个复杂且精细调控的过程，涉及多个关键步骤和多种酶的参与。miRNA的生物合成起始于细胞核内编码miRNA的基因（MIRgenes）的转录，在RNA聚合酶II的作用下，转录生成具有帽子结构和多聚腺苷酸尾巴的初级转录本（primarymiRNA，pri-miRNA）。pri-miRNA通常具有复杂的茎环结构，长度可达几百到几千个核苷酸。随后，pri-miRNA在DCL1（DICER-LIKE1）、HYL1（HYDRIN1）和SE（SERRATE）等蛋白组成的复合体作用下，经过两次切割，逐步加工成约70-100个核苷酸的前体miRNA（pre-miRNA），pre-miRNA具有典型的发夹结构。其中，DCL1是一种核酸酶，能够识别并切割pri-miRNA的茎环结构，将其加工成pre-miRNA；HYL1是一种双链RNA结合蛋白，能够与DCL1相互作用，增强DCL1对pri-miRNA的切割效率；SE是一种锌指蛋白，参与pre-miRNA的加工过程，对miRNA的生物合成具有重要的调控作用。pre-miRNA在转运蛋白HASTY的协助下，从细胞核转运到细胞质中。在细胞质中，pre-miRNA被DCL1进一步切割，形成长度约为21-24个核苷酸的成熟miRNA双链。成熟miRNA双链中的一条链会被选择性地整合到AGO（ARGONAUTE）蛋白中，形成RNA诱导沉默复合体（RNA-inducedsilencingcomplex，RISC）。AGO蛋白是RISC的核心组成部分，能够识别并结合miRNA，引导RISC靶向作用于互补的靶mRNA。miRNA主要通过两种方式发挥作用：切割靶mRNA和抑制靶mRNA的翻译。当miRNA与靶mRNA的互补序列完全或近乎完全匹配时，RISC中的AGO蛋白会切割靶mRNA，使其降解，从而实现对基因表达的负调控。例如，在拟南芥中，miR164能够与NAC1基因的mRNA互补配对，RISC识别并切割NAC1mRNA，抑制NAC1基因的表达，从而调控植物的侧根发育。当miRNA与靶mRNA的互补序列不完全匹配时，RISC主要通过抑制靶mRNA的翻译过程，减少蛋白质的合成，进而调控基因表达。在植物中，这种作用方式较为常见，许多miRNA通过这种方式调控植物的生长发育和对环境胁迫的响应。1.5.2MicroRNA对植物生长发育及抗逆性的调控MicroRNA（miRNA）作为一类重要的非编码RNA，在植物生长发育的各个阶段以及应对各种逆境胁迫的过程中，均发挥着不可或缺的调控作用，对维持植物的正常生理功能和生态适应性具有重要意义。在植物种子萌发阶段，miRNA参与调控种子休眠与萌发的平衡。例如，miR156和miR172通过调控相关靶基因的表达，影响种子的休眠与萌发进程。miR156通过靶向SPL（SQUAMOSAPROMOTERBINDINGPROTEIN-LIKE）转录因子家族成员，抑制其表达，维持种子的休眠状态；而在种子萌发时，miR156的表达水平下降，SPL基因表达上调，促进种子萌发。同时，miR172通过靶向AP2（APETALA2）转录因子，调控种子萌发相关基因的表达，影响种子的萌发率和萌发速度。在植物营养生长阶段，miRNA对叶片发育、茎的伸长、根系生长等过程具有重要调控作用。在叶片发育方面，miR164通过调控NAC1、NAC100等靶基因的表达，影响叶片的形态建成和衰老进程。在拟南芥中，miR164过表达植株的叶片表现出形态异常，如叶片变小、卷曲等，同时叶片衰老进程加快。miR396则通过靶向GRF（GROWTH-REGULATINGFACTOR）转录因子，调控叶片的细胞增殖和分化，影响叶片的大小和形状。在茎的伸长方面，miR159通过调控MYB33、MYB65等靶基因的表达，参与赤霉素信号通路，影响茎的伸长生长。在根系生长方面，miR160通过靶向ARF10（AUXINRESPONSEFACTOR10）、ARF16和ARF17等生长素响应因子，调控根系的生长和发育，影响根的长度、侧根的数量和根毛的发育。在植物生殖生长阶段，miRNA对花的发育、花粉育性、果实发育等过程至关重要。在花发育过程中，miR156-SPL模块参与调控植物的开花时间和花器官的发育。在拟南芥中，miR156的表达水平随着植物的生长逐渐降低，使得SPL基因的表达逐渐升高，从而促进植物从营养生长向生殖生长转变，调控开花时间。同时，miR172通过调控AP2等靶基因的表达，影响花器官的分化和发育，如调控萼片、花瓣、雄蕊和雌蕊的形成和发育。在花粉育性方面，miR164、miR167等miRNA通过调控相关靶基因的表达，影响花粉的发育和花粉管的生长，对花粉育性具有重要调控作用。在果实发育过程中，miRNA参与调控果实的大小、形状、色泽和成熟进程。在番茄中，miR164通过调控NAC1的表达，影响果实的大小和形状；miR172通过调控AP2的表达，参与果实的成熟调控，影响果实的色泽和风味。在植物应对非生物胁迫方面，miRNA在干旱、盐害、低温等逆境条件下发挥重要的调控作用，帮助植物提高抗逆性。在干旱胁迫下，植物体内的miR169、miR393、miR398等miRNA的表达水平发生显著变化。miR169通过靶向NF-YA（NUCLEARFACTOR-YALPHA）转录因子家族成员，抑制其表达，从而调控植物的干旱响应。miR393通过靶向TIR1（TRANSPORTINHIBITORRESPONSE1）和AFB2（AUXINSIGNALINGF-BOX2）等生长素受体基因，调控生长素信号通路，影响植物的生长和对干旱胁迫的适应。在盐胁迫下，miR160、miR169、miR393等miRNA参与调控植物的离子平衡、渗透调节和抗氧化防御等过程，增强植物的耐盐性。例如，miR160通过调控ARF10、ARF16和ARF17的表达，影响植物对盐胁迫的响应；miR169通过调控NF-YA的表达，参与植物的盐胁迫响应。在低温胁迫下，miR319、miR393、miR397等miRNA的表达水平改变，通过调控相关靶基因的表达，影响植物的低温响应和抗寒能力。miR319通过靶向TCP（TEOSINTEBRANCHED1/CYCLOIDEA/PCF）转录因子家族成员，调控植物的生长和对低温的适应；miR393通过调控TIR1和AFB2的表达，参与植物的低温胁迫响应。在植物应对生物胁迫方面，miRNA在植物抵御病原菌入侵过程中发挥重要作用，激活植物的免疫反应。在植物与病原菌互作过程中，病原菌侵染会诱导植物体内miRNA表达谱的改变，从而调控植物防御相关基因的表达。例如，在拟南芥与丁香假单胞菌互作过程中，miR393的表达显著上调，通过靶向TIR1和AFB2，抑制生长素信号通路，激活植物的免疫反应，增强植物对病原菌的抗性。miR168通过靶向AGO1，调控植物的RNA沉默途径，参与植物对病毒的防御反应。1.5.3MicroRNA对花青苷合成的调控MicroRNA（miRNA）在植物花青苷合成调控中发挥着重要作用，通过靶向调控花青苷合成相关基因的表达，影响花青苷的生物合成和积累，进而决定植物器官的色泽。众多研究表明，miRNA对花青苷合成的调控是一个复杂的网络，涉及多个miRNA及其靶基因之间的相互作用。在苹果中，miR172被发现通过靶向调控AP2转录因子来影响花青苷的合成。AP2转录因子是花青苷合成调控网络中的重要组成部分，能够抑制MYB10等花青苷合成正调控因子的表达。研究发现，过表达miR172的苹果转基因植株在叶、茎和果皮中呈现出花青苷积累减少的现象，而抑制miR172的表达则会促进花青苷的积累。这是因为miR172能够识别并切割AP2转录因子的mRNA，使其降解，从而解除AP2对MYB10的抑制作用，促进花青苷的合成。进一步研究表明，miR172-AP2-MYB10调控模块在不同苹果品种中具有保守性，对苹果果实色泽的形成起着关键作用。在葡萄中，miR858是调控花青苷合成的重要miRNA。miR858通过靶向MYB转录因子家族成员，抑制其表达，从而负调控花青苷的合成。在葡萄果实发育过程中，miR858的表达水平与花青苷含量呈负相关。在花青苷积累阶段，miR858的表达水平下降，使得其靶基因MYB转录因子的表达上调，激活花青苷合成相关结构基因（如CHS、CHI、F3H、DFR、ANS、UFGT等）的表达，促进花青苷的合成。研究还发现，miR858的表达受到光照、温度等环境因素的调控，进一步说明了miR858在葡萄花青苷合成调控中的重要性。在海棠中，低氮条件下，由非编码基因与编码基因共同组建的eTM-miR858-MYB62-like模块参与调控花青苷的生物合成。研究表明，miR858可以通过转录切割的方式负调控其下游靶基因MsMYB62-like。MsMYB62-like过表达可以显著降低总花青苷含量，而后通过酵母单杂交和双荧光素酶报告基因测定等实验结果表明，MsMYB62-like可以抑制花青苷生物合成通路中关键结构基因MsF3′H的表达，从而负向调节花青苷的生物合成。此外，eTM858-1和eTM858-2（lncRNA）通过靶标模拟结合miR858来起竞争性内源RNA的作用，降低miR858对MsMYB62-like的切割作用，进而间接负调控海棠叶片花青苷的生物合成。这一研究揭示了低氮条件下miRNA调控花青苷合成的新机制，拓展了对植物花青苷合成调控网络的认识。在梨中，虽然关于miRNA调控花青苷合成的研究相对较少，但已有研究表明miRNA在梨果皮花青苷合成调控中具有潜在作用。通过对红皮梨和绿皮梨果实发育过程中miRNA表达谱的分析，发现一些差异表达的miRNA，推测这些miRNA可能参与调控梨果皮花青苷的合成。例如，在‘红早酥’梨中，发现miR156、miR164等miRNA的表达水平与花青苷含量存在一定的相关性。进一步研究这些miRNA的靶基因及其调控机制，将有助于深入了解梨果皮花青苷合成的调控网络。1.6梨果皮着色的研究进展1.6.1红梨花青苷的分布及着色模式红梨花青苷在果皮中的分布呈现出特定的规律，对果实的外观色泽起着决定性作用。研究表明，花青苷主要积累在果皮的表皮细胞和亚表皮细胞中。在‘红早酥’梨中，通过组织化学染色和显微镜观察发现，随着果实的发育，花青苷在表皮细胞中的积累逐渐增加，使果皮颜色由浅变深。这种分布特点与花青苷合成相关基因的表达模式密切相关，相关基因在表皮细胞和亚表皮细胞中的高表达，促进了花青苷的合成和积累。不同红梨品种的着色模式存在显著差异，这主要受遗传因素和环境因素的共同影响。根据着色时间和着色程度，红梨品种的着色模式可大致分为早着色型、中着色型和晚着色型。早着色型品种如‘早酥红’，在果实发育早期就开始积累花青苷，果皮迅速变红，且着色均匀，色泽鲜艳；中着色型品种如‘红香酥’，花青苷积累相对较晚，在果实发育中期开始明显着色，果皮颜色逐渐加深；晚着色型品种如‘满天红’，花青苷积累最晚，在果实接近成熟时才大量合成，果皮在后期迅速变红。环境因素对红梨品种的着色模式也有重要影响。光照是影响红梨着色的关键环境因素之一，充足的光照有利于花青苷的合成和积累，促进果实着色。在光照充足的地区，红梨果实的着色时间提前，色泽更加鲜艳；而在光照不足的情况下，果实着色延迟，色泽暗淡。温度也会影响红梨的着色模式，适宜的温度范围（20-25℃）有利于花青苷的合成，过高或过低的温度都会抑制花青苷的合成，影响果实着色。此外，土壤肥力、水分管理等栽培措施也会间接影响红梨的着色模式，合理的施肥和灌溉能够为果实生长提供充足的养分和水分，促进花青苷的合成，改善果实着色。1.6.2红梨花青苷合成的生理及分子基础红梨花青苷的合成是一个复杂的生理过程，涉及一系列酶促反应和代谢调控。在生理层面，花青苷的合成与植物的光合作用、呼吸作用以及碳水化合物代谢密切相关。光合作用为花青苷合成提供了充足的能量和碳源，光合产物（如蔗糖、葡萄糖等）是花青苷合成的前体物质。在红梨果实发育过程中，随着光合作用的增强，果实中可溶性糖含量增加，为花青苷的合成提供了物质基础。同时，呼吸作用也参与了花青苷合成的调控，呼吸作用产生的能量和中间产物（如乙酰辅酶A等）为花青苷合成的酶促反应提供了必要条件。从分子层面来看，红梨花青苷的合成受到一系列结构基因和转录因子的调控。如前文所述，花青苷合成途径中的结构基因（如PpCHS、PpCHI、PpF3H、PpF3'H、PpDFR、PpANS和PpUFGT等）编码的酶参与了花青苷合成的各个步骤，它们的表达水平直接影响花青苷的合成速率和积累量。在‘红早酥’梨果实发育过程中，PpDFR、PpANS和PpUFGT等基因的表达量与花青苷含量呈显著正相关，这些基因的高表达促进了花青苷的合成和积累。转录因子在红梨花青苷合成调控中起着关键作用，它们通过与结构基因启动子区域的顺式作用元件结合，激活或抑制结构基因的转录，从而调控花青苷的合成。在梨中，PpMYB10是调控花青苷合成的关键MYB转录因子，它能够直接结合到花青苷合成相关结构基因（如PpDFR、PpANS和PpUFGT）的启动子区域，激活这些基因的表达，促进花青苷的合成。此外，bHLH和WD40转录因子家族成员与PpMYB10相互作用，形成MBW复合体，协同调控花青苷的合成。除了MBW复合体，其他转录因子家族成员如WRKY、NAC、ERF等也参与了红梨花青苷合成的调控。在‘满天红’梨中，PpWRKY44通过直接结合PpMYB10启动子上的W-box激活其表达，进而诱导花青苷合成。PpNAC1能够直接结合到PpMYB10的启动子区域，抑制其表达，从而负调控花青苷的合成。1.6.3红梨着色与光相关的人工调控措施在红梨生产过程中，为了提高果实的着色质量和商品价值，常采用一系列与光相关的人工调控措施，其中套袋和光照处理是最为常用且有效的方法。套袋是一种广泛应用于红梨栽培的技术，通过在果实发育初期将果实套上特制的纸袋，能够有效改善果实的色泽和品质。套袋可以调节果实周围的微环境，减少病虫害的侵袭，降低农药残留，同时还能避免果实受到强光直射和机械损伤。在红梨套袋过程中，选择合适的纸袋类型和套袋时间至关重要。纸袋的材质、颜色和透光率等因素都会影响套袋效果。一般来说，双层纸袋比单层纸袋具有更好的遮光和保温性能，更有利于果实着色。对于红梨品种‘红香酥’，在花后30-40天套上双层纸袋，能够显著提高果实的着色率和色泽鲜艳度。套袋时间过早或过晚都会影响果实的生长发育和着色效果，过早套袋可能会导致果实生长缓慢，过晚套袋则无法充分发挥套袋对果实色泽的改善作用。在果实发育后期，适时摘袋能够让果实接受光照，促进花青苷的合成和积累，从而使果实更好地着色。摘袋时间的选择需要综合考虑多种因素，如品种特性、当地气候条件和果实成熟度等。对于大多数红梨品种，在果实成熟前15-20天摘袋较为适宜。在光照充足、昼夜温差大的地区，适当提前摘袋可以让果实充分利用光照资源，促进花青苷的合成；而在光照不足或气候条件不稳定的地区，可适当推迟摘袋时间，以避免果实受到不良环境的影响。除了套袋和摘袋，光照处理也是调控红梨着色的重要手段。通过人工补光或调整光照时间和强度，可以有效促进花青苷的合成，改善果实的色泽。在设施栽培中，常采用LED灯进行补光，不同光质（如红光、蓝光、白光等）对红梨果实着色的影响不同。研究表明，红光和蓝光能够显著促进花青苷的合成，提高果实的着色程度。在‘早酥红’梨果实发育后期，用红光和蓝光组合照射，可使果实中花青苷含量显著增加，果皮颜色更加鲜艳。此外，延长光照时间也能促进花青苷的合成，在一定范围内，光照时间越长，花青苷的积累量越高。然而，过强的光照或过长的光照时间也可能对果实造成伤害，如导致果实灼伤、品质下降等。因此，在进行光照处理时，需要根据红梨品种的特性和生长环境，合理调整光照参数，以达到最佳的着色效果。二、DNA甲基化调控红梨果皮着色的机制研究2.1材料与方法2.1.1实验材料本研究选取在[具体种植地点]种植的‘红早酥’红梨为实验材料，该品种为[具体品种来源]，具有典型的红梨果皮着色特征，在当地广泛种植且表现出稳定的生长和果实品质特性。在果实发育的不同时期，包括花后30天、60天、90天、120天和150天，分别从生长健壮、无病虫害且树势一致的植株上采集果实样本。每个时期选取3株树，每株树采集5个果实，确保样本具有代表性。采集后的果实迅速用液氮冷冻，并储存于-80℃冰箱中备用，以防止RNA和DNA的降解以及相关生理生化变化。2.1.2实验技术花青苷含量测定：采用pH示差法测定红梨果皮中的花青苷含量。具体步骤如下：取0.5g红梨果皮样品，在液氮中研磨成粉末，加入5mL含1%HCl的甲醇溶液，于4℃下避光萃取12h。萃取结束后，10000r/min离心20min，取上清液。分别取1mL上清液，加入4mL缓冲液A（0.4mol/LKCl，pH1.0）和缓冲液B（0.4mol/L柠檬酸，pH4.5），混合均匀，在510nm和700nm波长下测定吸光值。根据公式计算花青苷含量：花青苷含量（mg/gFW）=[(A510-A700)×MW×DF×V]/(ε×m)，其中A510和A700分别为510nm和700nm处的吸光值，MW为花青苷的平均分子量（449.2g/mol），DF为稀释倍数，V为提取液总体积（mL），ε为花青苷的摩尔消光系数（26900L/mol/cm），m为样品质量（g）。DNA提取与甲基化分析：使用植物基因组DNA提取试剂盒（[具体品牌和型号]）提取红梨果皮基因组DNA。提取的DNA经琼脂糖凝胶电泳检测其完整性，并用NanoDrop2000超微量分光光度计测定其浓度和纯度，确保A260/A280在1.8-2.0之间，A260/A230大于2.0。采用全基因组重亚硫酸盐测序（WGBS）技术分析DNA甲基化水平。将提取的基因组DNA进行重亚硫酸盐转化，使未甲基化的胞嘧啶转化为尿嘧啶，而甲基化的胞嘧啶保持不变。转化后的DNA进行PCR扩增和文库构建，然后在IlluminaHiSeq测序平台上进行测序。测序数据经过质量控制和比对分析，利用相关软件（如Bismark、MethylKit等）计算不同区域（启动子、基因编码区、基因间区等）的DNA甲基化水平，包括CG、CHG和CHH三种甲基化类型，并分析甲基化水平在果实发育不同时期的变化。基因表达分析：采用Trizol试剂提取红梨果皮总RNA，经DNaseI处理去除基因组DNA污染。使用反转录试剂盒（[具体品牌和型号]）将RNA反转录为cDNA。以cDNA为模板，利用实时荧光定量PCR（qRT-PCR）技术分析花青苷合成相关基因（如PpMYB10、PpDFR、PpANS、PpUFGT等）以及DNA甲基转移酶基因（如PpMET1、PpCMT3、PpDRM2）和去甲基化酶基因（如PpROS1、PpDME）的表达水平。qRT-PCR反应体系为20μL，包括10μLSYBRGreenMasterMix、0.5μL上游引物（10μmol/L）、0.5μL下游引物（10μmol/L）、1μLcDNA模板和8μLddH₂O。反应条件为：95℃预变性30s，然后95℃变性5s，60℃退火30s，共40个循环。以梨的Actin基因作为内参基因，采用2^(-ΔΔCt)法计算基因的相对表达量。引物序列根据梨基因组数据库（[具体数据库名称和网址]）设计，并通过PrimerPremier5.0软件进行引物特异性和扩增效率验证，引物序列见表1。基因名称引物序列（5'-3'）PpMYB10F:ATGGCTCCGCTGCTGCTGR:TTACACGCTGCTGCTGCTGPpDFRF:GGCGGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpANSF:AAGAGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpUFGTF:GGCGGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpMET1F:AAGAGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpCMT3F:GGCGGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpDRM2F:AAGAGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpROS1F:GGCGGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGPpDMEF:AAGAGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCGActinF:GGCGGAGAGAGAGAGAGAGR:CCCGCCGCCGCCGCCG甲基化抑制剂处理：为进一步验证DNA甲基化对红梨果皮花青苷合成的调控作用，选取花后90天的红梨果实，采用甲基化抑制剂5-氮杂胞苷（5-azaC）进行处理。将5-azaC溶解于无菌水中，配制成100μmol/L的溶液。在果实赤道部位均匀涂抹5-azaC溶液，每果涂抹0.5mL，以涂抹无菌水的果实作为对照。处理后的果实继续在植株上生长，分别在处理后3天、6天和9天采集果皮样本，用于花青苷含量测定、DNA甲基化分析和基因表达分析，以研究5-azaC处理对红梨果皮花青苷合成及相关基因表达和DNA甲基化水平的影响。2.2结果与分析2.2.1不同组织/器官的花青苷含量对‘红早酥’红梨的不同组织/器官，包括幼叶、成熟叶、花瓣、雄蕊、雌蕊、幼果果皮、成熟果果皮和果肉进行花青苷含量测定，结果显示（图1），花青苷含量在不同组织/器官中存在显著差异（P<0.05）。成熟果果皮中的花青苷含量最高，达到（5.67±0.32）mg/gFW，显著高于其他组织/器官。幼果果皮中花青苷含量次之，为（3.45±0.21）mg/gFW。花瓣中也检测到一定含量的花青苷，为（1.23±0.15）mg/gFW，而幼叶、成熟叶、雄蕊、雌蕊和果肉中的花青苷含量相对较低，分别为（0.25±0.05）mg/gFW、（0.18±0.03）mg/gFW、（0.12±0.02）mg/gFW、（0.10±0.01）mg/gFW和（0.08±0.01）mg/gFW。这表明花青苷在‘红早酥’红梨的不同组织/器官中的积累具有明显的组织特异性，成熟果果皮是花青苷积累的主要部位，可能与果实的色泽形成和对环境的适应等功能相关。注：不同小写字母表示在P<0.05水平上差异显著。2.2.2不同组织/器官中花青苷合成相关基因的表达采用qRT-PCR技术检测花青苷合成相关基因（PpMYB10、PpDFR、PpANS、PpUFGT）在‘红早酥’红梨不同组织/器官中的表达水平，结果如图2所示。PpMYB10在成熟果果皮中的表达量最高，显著高于其他组织/器官（P<0.05），在幼果果皮和花瓣中也有较高的表达水平，而在幼叶、成熟叶、雄蕊、雌蕊和果肉中的表达量相对较低。PpDFR、PpANS和PpUFGT基因的表达模式与PpMYB10相似，在成熟果果皮中的表达量最高，在幼果果皮和花瓣中也有较高表达，在其他组织/器官中表达量较低。相关性分析表明，PpMYB10、PpDFR、PpANS和PpUFGT基因的表达水平与花青苷含量呈显著正相关（r>0.8，P<0.01），说明这些基因在‘红早酥’红梨的花青苷合成过程中起着重要作用，其高表达促进了花青苷在成熟果果皮等组织/器官中的积累。注：不同小写字母表示在P<0.05水平上差异显著。2.2.3系统聚类分析对花青苷合成相关基因在不同组织/器官中的表达数据进行系统聚类分析，结果如图3所示。根据基因表达模式的相似性，可将不同组织/器官分为两组。第一组包括成熟果果皮和幼果果皮，这两个组织中花青苷合成相关基因的表达水平较高；第二组包括幼叶、成熟叶、花瓣、雄蕊、雌蕊和果肉，这些组织中花青苷合成相关基因的表达水平相对较低。在第一组中，成熟果果皮和幼果果皮的基因表达模式也存在一定差异，成熟果果皮中PpMYB10、PpDFR、PpANS和PpUFGT基因的表达量均高于幼果果皮。在第二组中，花瓣的基因表达模式与其他组织略有不同，花瓣中PpMYB10和PpDFR的表达量相对较高。系统聚类分析结果进一步表明，花青苷合成相关基因在‘红早酥’红梨不同组织/器官中的表达具有明显的组织特异性，且成熟果果皮和幼果果皮在花青苷合成过程中具有相似但又有差异的分子调控机制。2.2.4不同组织/器官中关键基因编码区和启动子序列比较对‘红早酥’红梨不同组织/器官中关键基因（PpMYB10、PpDFR、PpANS、PpUFGT）的编码区和启动子序列进行PCR扩增和测序分析，结果显示，在不同组织/器官中，这些关键基因的编码区序列高度保守，未发现碱基突变或缺失等情况。然而，启动子序列分析发现，不同组织/器官中关键基因的启动子序列存在一定差异。以PpMYB10基因启动子为例，在成熟果果皮和幼果果皮中，启动子区域存在一些特定的顺式作用元件，如光响应元件（ACE、G-box等）、激素响应元件（ABRE、ERE等）和转录因子结合位点（MYB结合位点、bHLH结合位点等），这些顺式作用元件在其他组织/器官中的分布和数量存在差异。在成熟果果皮中，光响应元件ACE和G-box的数量较多，且MYB结合位点和bHLH结合位点的序列更为保守，这可能与成熟果果皮在光照条件下高效合成花青苷的能力相关。不同组织/器官中关键基因启动子序列的差异，可能导致转录因子与启动子的结合能力不同，从而影响基因的表达水平，进而调控花青苷在不同组织/器官中的合成和积累。2.2.5关键基因启动子和转录因子互作分析采用酵母单杂交和双荧光素酶报告系统验证关键基因（PpMYB10、PpDFR、PpANS、PpUFGT）启动子与转录因子（PpMYB10、PpbHLH、PpWD40）之间的相互作用。酵母单杂交结果显示，PpMYB10转录因子能够与PpDFR、PpANS和PpUFGT基因的启动子特异性结合，而PpbHLH和PpWD40转录因子与PpMYB10形成复合体后，增强了PpMYB10与启动子的结合能力。双荧光素酶报告系统分析进一步证实，PpMYB10、PpbHLH和PpWD40组成的MBW复合体能够显著激活PpDFR、PpANS和PpUFGT基因启动子的活性，促进荧光素酶的表达，表明MBW复合体在‘红早酥’红梨的花青苷合成调控中发挥着重要作用。此外，通过定点突变实验发现，当PpMYB10基因启动子中的MYB结合位点或bHLH结合位点发生突变时，MBW复合体与启动子的结合能力显著降低，基因启动子的活性也明显下降，进一步证明了MBW复合体与关键基因启动子之间的特异性相互作用对花青苷合成相关基因表达的调控作用。2.2.6不同组织/器官中关键基因启动子的甲基化水平检测利用全基因组重亚硫酸盐测序（WGBS）技术检测‘红早酥’红梨不同组织/器官中关键基因（PpMYB10、PpDFR、PpANS、PpUFGT）启动子的甲基化水平，结果如图4所示。在不同组织/器官中，关键基因启动子的甲基化水平存在显著差异（P<0.05）。成熟果果皮中PpMYB10、PpDFR、PpANS和PpUFGT基因启动子的甲基化水平最低，分别为（5.6±1.2）%、（6.8±1.5）%、（7.2±1.8）%和（8.1±2.0）%，而在幼叶、成熟叶、雄蕊、雌蕊和果肉中，这些基因启动子的甲基化水平相对较高。相关性分析表明，关键基因启动子的甲基化水平与基因表达水平呈显著负相关（r<-0.8，P<0.01），与花青苷含量也呈显著负相关（r<-0.8，P<0.01）。这表明在‘红早酥’红梨中，DNA甲基化可能通过抑制关键基因启动子的活性，降低基因的表达水平，从而抑制花青苷的合成和积累，成熟果果皮中较低的甲基化水平有利于花青苷合成相关基因的表达，促进花青苷的积累，使果皮呈现红色。注：不同小写字母表示在P<0.05水平上差异显著。2.3讨论本研究对‘红早酥’红梨不同组织/器官的花青苷含量及相关基因表达进行分析，发现花青苷在成熟果果皮中含量最高，且花青苷合成相关基因（PpMYB10、PpDFR、PpANS、PpUFGT）的表达水平与花青苷含量呈显著正相关。这与前人在苹果、葡萄等果实中的研究结果一致，表明这些基因在红梨花青苷合成过程中发挥着关键作用。系统聚类分析将不同组织/器官分为两组，进一步验证了花青苷合成相关基因表达的组织特异性，成熟果果皮和幼果果皮在花青苷合成调控方面具有相似性但也存在差异，这可能与果实发育阶段及组织功能有关。对关键基因编码区和启动子序列的分析发现，编码区序列高度保守，而启动子序列存在差异，尤其是成熟果果皮中关键基因启动子含有更多与光响应、激素响应及转录因子结合相关的顺式作用元件，这与成熟果果皮在光照条件下高效合成花青苷相契合。启动子与转录因子互作分析证实，PpMYB10、PpbHLH和PpWD40组成的MBW复合体能够特异性结合并激活PpDFR、PpANS和PpUFGT基因启动子，促进花青苷合成相关基因的表达，这是红梨花青苷合成调控的重要分子机制之一。DNA甲基化分析结果表明，关键基因启动子的甲基化水平在不同组织/器官中存在显著差异，且与基因表达水平及花青苷含量呈显著负相关。成熟果果皮中关键基因启动子甲基化水平最低，有利于基因表达和花青苷积累，而其他组织/器官中较高的甲基化水平抑制了基因表达和花青苷合成。这与在苹果、葡萄等植物中的研究结果相符，进一步证明DNA甲基化在红梨花青苷合成调控中起重要作用，通过调控关键基因启动子的甲基化状态，影响基因的表达，从而实现对花青苷合成的调控。甲基化抑制剂5-azaC处理实验进一步验证了DNA甲基化对红梨花青苷合成的调控作用。5-azaC处理降低了关键基因启动子的甲基化水平，促进了基因表达和花青苷积累，使果皮颜色变红。这表明降低DNA甲基化水平能够解除对花青苷合成相关基因的抑制，促进花青苷合成，为通过调控DNA甲基化水平改善红梨果实着色提供了理论依据和实践指导。然而，本研究也存在一定的局限性，仅对‘红早酥’红梨进行了研究，未来需要进一步扩大研究品种范围，深入探究DNA甲基化在不同红梨品种中的调控机制差异。同时，对于DNA甲基化与其他调控因子（如转录因子、miRNA等）之间的互作关系研究还不够深入，需要进一步开展相关研究，以完善红梨花青苷合成的调控网络。三、MicroRNA调控红梨果皮着色的机制研究3.1材料与方法本研究选取在[具体种植地点]种植的‘红早酥’红梨为实验材料，该品种为[具体品种来源]，在当地种植表现良好，具有典型的红梨果皮着色特征。在果实发育的不同时期，包括花后30天、60天、90天、120天和150天，从生长健壮、无病虫害且树势一致的植株上采集果实样本。每个时期选取3株树，每株树采集5个果实，采集后的果实迅速用液氮冷冻，并储存于-80℃冰箱中备用。在实验技术方面，采用Trizol试剂提取红梨果皮总RNA，使用miRcutemiRNA提取分离试剂盒（[具体品牌和型号]）分离小RNA，经质量检测合格后，进行miRNA文库构建。利用IlluminaHiSeq测序平台对文库进行高通量测序，测序数据经过质量控制和生物信息学分析，鉴定差异表达的miRNA，并预测其靶基因。采用5'-RACE（RapidAmplificationofcDNAEnds）技术验证miRNA对靶基因的切割位点，进一步确认miRNA与靶基因的靶向关系。利用psRNATarget软件对miRNA的靶基因进行预测，并结合GO（GeneOntology）和KEGG（KyotoEncyclopediaofGenesandGenomes）数据库对靶基因进行功能注释和代谢通路分析，了解miRNA调控的生物学过程和相关代谢途径。通过茎环法反转录和实时荧光定量PCR技术分析差异表达miRNA及其靶基因（如PpMYB10、PpDFR、PpANS、PpUFGT等花青苷合成相关基因）在果实发育不同时期的表达水平，以U6作为内参基因，采用2^(-ΔΔCt)法计算基因的相对表达量，分析miRNA与靶基因表达水平的相关性。为验证miRNA对花青苷合成的调控作用，构建miRNA过表达载体和靶基因过表达载体，通过农杆菌介导的方法转化烟草叶片或红梨愈伤组织，分别过表达miRNA和靶基因，观察转化后组织中花青苷含量的变化及相关基因的表达水平。同时，设计合成miRNA抑制剂，通过浸泡或注射的方法导入红梨果实中，抑制miRNA的表达，分析对花青苷合成及相关基因表达的影响。3.2结果与分析3.2.1梨上全基因组miRNA及其靶基因的鉴定、特定miRNA前体启动子顺式作用元件的分析通过对‘红早酥’红梨果皮小RNA文库的高通量测序及生物信息学分析，共鉴定出[X]个已知miRNA，分属于[X]个miRNA家族，同时发现了[X]个新的miRNA。对这些miRNA的靶基因进行预测，共获得[X]个靶基因，涉及多种生物学过程，如转录调控、信号转导、代谢过程等。其中，与花青苷合成相关的靶基因有[X]个，包括PpMYB10、PpDFR、PpANS、PpUFGT等。对miR156前体启动子顺式作用元件的分析结果显示，其启动子区域存在多个与光响应、激素响应和转录因子结合相关的顺式作用元件。其中，光响应元件包括G-box、ACE、GT1-motif等，激素响应元件有ABRE（脱落酸响应元件）、ERE（乙烯响应元件）、TGA-element（生长素响应元件）等。转录因子结合位点包含MYB结合位点、bHLH结合位点等。这些顺式作用元件的存在，暗示miR156的表达可能受到光、激素以及其他转录因子的调控，进而影响其对靶基因的调控作用，参与红梨果皮花青苷合成的调控。3.2.2梨全基因组SPL的鉴定利用生物信息学方法，在梨全基因组中鉴定出[X]个SPL基因，命名为PpSPL1-PpSPL[X]。对这些基因的基本特征进行分析，发现它们的开放阅读框长度在[X]-[X]bp之间，编码的氨基酸数目在[X]-[X]个之间。基因结构分析表明，PpSPL基因家族成员的内含子数目存在差异，从[X]个到[X]个不等。进化树分析将PpSPL基因分为[X]个亚家族，不同亚家族的PpSPL基因在结构和功能上可能存在差异。染色体定位分析显示，这些PpSPL基因分布在梨的[X]条染色体上，且分布不均匀，部分染色体上存在基因簇现象。3.2.3梨上全基因组特定miRNA靶基因的鉴定针对miR156，通过psRNATarget软件预测其在梨全基因组中的靶基因，共得到[X]个靶基因，其中包括多个PpSPL基因，如PpSPL3、PpSPL9、PpSPL10等。利用5'-RACE技术对miR156与PpSPL9的靶向关系进行验证，结果表明miR156能够特异性地切割PpSPL9的mRNA，在预测的切割位点处发生了mRNA的断裂，进一步证实了miR156对PpSPL9的靶向调控作用。GO功能注释分析显示，miR156靶基因主要富集在转录调控、花发育、衰老过程等生物学过程，KEGG代谢通路分析表明这些靶基因参与了植物激素信号转导、淀粉和蔗糖代谢等多条代谢通路，暗示miR156通过调控这些靶基因，参与红梨的生长发育和花青苷合成等过程。3.2.4去袋后红梨果皮花青苷积累、花青苷合成相关基因、SPL基因及miRNA的表达以‘美人酥’红梨为材料，研究去袋后不同时间（0d、3d、6d、9d、12d）果皮花青苷积累、花青苷合成相关基因（PpMYB10、PpDFR、PpANS、PpUFGT）、PpSPL基因（PpSPL3、PpSPL9、PpSPL10）及miR156的表达变化。结果显示，去袋后随着时间的延长，果皮花青苷含量逐渐增加，在去袋后12d达到最高值，为（[X]±[X]）mg/gFW。花青苷合成相关基因PpMYB10、PpDFR、PpANS、PpUFGT的表达水平也逐渐上调，在去袋后12d表达量显著高于去袋后0d。PpSPL3、PpSPL9、PpSPL10基因的表达水平在去袋后呈现先下降后上升的趋势，在去袋后6d表达量最低。miR156的表达水平则在去袋后逐渐上升，在去袋后12d达到最高值。相关性分析表明，miR156的表达水平与PpSPL3、PpSPL9、PpSPL10基因的表达水平呈显著负相关（r<-0.8，P<0.01），与花青苷含量及花青苷合成相关基因的表达水平呈显著正相关（r>0.8，P<0.01）。这表明去袋后光照诱导miR156表达上调，进而抑制PpSPL基因的表达，促进花青苷合成相关基因的表达和花青苷的积累，调控红梨果皮着色。3.2.5花青苷合成相关转录因子的互作分析采用酵母双杂交和双荧光素酶报告系统分析花青苷合成相关转录因子之间的互作关系。酵母双杂交结果显示，PpMYB10与PpbHLH、PpWD40之间存在相互作用，能够形成MBW复合体。进一步通过双荧光素酶报告系统验证，将PpMYB10、PpbHLH、PpWD40与花青苷合成相关结构基因（PpDFR、PpANS、PpUFGT）启动子的重组载体共转化烟草叶片，结果表明MBW复合体能够显著激活PpDFR、PpANS、PpUFGT基因启动子的活性，促进荧光素酶的表达，说明MBW复合体在红梨花青苷合成调控中发挥着关键作用。此外，通过定点突变实验发现，当PpMYB10基因中与PpbHLH、PpWD40相互作用的结构域发生突变时，MBW复合体的形成受到抑制，对花青苷合成相关基因启动子的激活作用也明显减弱。3.2.6SPL蛋白与花青苷合成相关转录因子的互作分析为探究SPL蛋白与花青苷合成相关转录因子的互作关系，利用酵母双杂交和双分子荧光互补（BiFC）技术进行分析。酵母双杂交结果表明，PpSPL9能够与PpMYB10相互作用，且这种相互作用具有特异性。BiFC实验进一步在烟草叶片细胞中验证了PpSPL9与PpMYB10在细胞核内发生相互作用，形成蛋白复合体。通过荧光共定位分析发现，PpSPL9-PpMYB10复合体与花青苷合成相关结构基因（PpDFR、PpANS、PpUFGT）启动子在细胞核内存在共定位现象。进一步的转录激活实验表明，PpSPL9与PpMYB10形成的复合体能够抑制PpDFR、PpANS、PpUFGT基因启动子的活性，降低荧光素酶的表达，说明PpSPL9与PpMYB10的相互作用负调控花青苷合成相关基因的表达，从而抑制花青苷的合成。3.3讨论本研究通过高通量测序及生物信息学分析，在‘红早酥’红梨中鉴定出大量miRNA及其靶基因，其中与花青苷合成相关的靶基因有PpMYB10、PpDFR、PpANS、PpUFGT等。对miR156前体启动子顺式作用元件的分析发现，其启动子区域存在多个与光响应、激素响应和转录因子结合相关的顺式作用元件，暗示miR156的表达可能受到多种因素的调控，进而影响花青苷的合成。这与前人在苹果、葡萄等果实中的研究结果相似，表明miRNA在果实花青苷合成调控中具有重要作用。在梨全基因组中鉴定出多个SPL基因，其中PpSPL3、PpSPL9、PpSPL10等被预测为miR156的靶基因，且5'-RACE验证了miR156对PpSPL9的靶向调控作用。GO和KEGG分析表明，miR156靶基因参与多种生物学过程和代谢通路，为深入理解miR156在红梨生长发育和花青苷合成中的调控机制提供了线索。对‘美人酥’红梨去袋后花青苷积累及相关基因和miRNA表达的研究发现，miR156表达上调与花青苷积累及花青苷合成相关基因表达呈正相关，与PpSPL基因表达呈负相关。这表明去袋后光照诱导miR156表达，抑制PpSPL基因表达，从而解除对花青苷合成相关基因的抑制，促进花青苷积累，调控红梨果皮着色。花青苷合成相关转录因子互作分析表明，PpMYB10、PpbHLH和PpWD40形成的MBW复合体在红梨花青苷合成调控中发挥关键作用。而PpSPL9与PpMYB10的相互作用负调控花青苷合成相关基因的表达，抑制花青苷合成。这揭示了SPL蛋白在花青苷合成调控中的新机制，丰富了对红梨花青苷合成调控网络的认识。然而，本研究仅初步揭示了miR156及其靶基因对红梨果皮花青苷合成的调控机制，对于其他miRNA及其靶基因在花青苷合成中的作用研究较少。未来需要进一步深入研究不同miRNA及其靶基因之间的相互作用，以及它们与其他调控因子（如DNA甲基化、转录因子等）之间的协同调控关系，以完善红梨花青苷合成的调控网络。同时，本研究主要在实验室条件