中间锦鸡儿MicroRNA的挖掘、靶基因鉴定及其表达调控网络解析

中间锦鸡儿MicroRNA的挖掘、靶基因鉴定及其表达调控网络解析一、引言1.1研究背景与意义中间锦鸡儿（CaraganaintermediaKuangetH.C.Fu）作为一种旱生落叶灌木，在生态保护方面发挥着极为重要的作用。它广泛分布于中国北部的荒漠、半荒漠地区，具备强大的耐旱、耐寒、耐瘠薄能力，并且能够在风沙环境中稳定生长，有效起到防风固沙、保持水土的作用。在当前全球生态环境面临严峻挑战，土地沙漠化、水土流失等问题日益加剧的背景下，中间锦鸡儿对于维护生态平衡、改善生态环境的价值愈发凸显。例如在毛乌素沙地，中间锦鸡儿作为主要的固沙植物之一，其林龄与叶功能性状、群落特征之间的关系，对于了解其生态适应性及在沙地生态系统中的作用具有重要意义。研究发现，随着林龄的增加，中间锦鸡儿的叶绿素含量、比叶重等叶功能性状呈现出一定的变化规律，其与其他植物的竞争与协同作用，促进了群落的演替与发展。在高寒沙地，对不同恢复年限的中间锦鸡儿固沙林的研究表明，随着恢复年限的增加，植物的物种和功能群逐渐增多，这对于高寒沙地的治理和生态恢复具有重要的科学指导意义。微小核糖核酸（microRNA，miRNA）作为一类内生的、长度约为22个核苷酸的非编码单链RNA分子，在植物的生长发育以及应对各种生物和非生物胁迫过程中起着关键的调控作用。它们通过与靶mRNA的互补配对，实现对靶基因表达的精确调控，其作用方式包括切割靶mRNA、抑制翻译过程以及参与DNA甲基化等。例如，在植物应对干旱胁迫时，某些miRNA能够通过调控相关靶基因的表达，增强植物的抗旱能力；在植物的生长发育过程中，miRNA也参与调控了根、茎、叶、花等器官的发育进程。深入研究中间锦鸡儿中的miRNA及其靶基因，对于全面揭示中间锦鸡儿的生长发育调控机制、抗逆分子机制以及挖掘优良基因资源具有不可忽视的重要意义。从生长发育调控机制角度来看，了解miRNA如何调控中间锦鸡儿的各个生长阶段，有助于我们更好地掌握其生长规律，为人工培育和种植提供科学依据。在抗逆分子机制方面，明确miRNA在中间锦鸡儿应对干旱、寒冷、风沙等逆境条件下的作用，能够为提高其抗逆性提供新的思路和方法。挖掘中间锦鸡儿中的优良基因资源，对于培育具有更强适应性和生态价值的新品种具有重要的应用前景，有望在生态修复和荒漠化治理中发挥更大的作用。1.2miRNA研究进展1.2.1miRNA的发现与特性miRNA的发现历程是一个充满探索与惊喜的过程。1993年，VictorAmbros和GaryRuvkun等人在秀丽隐杆线虫（Caenorhabditiselegans）中首次发现了lin-4基因，它并不编码蛋白质，而是转录产生一种长度约为22个核苷酸的非编码RNA，能够通过与lin-14mRNA的3'非翻译区互补配对，在转录后水平调控lin-14的表达，从而影响线虫的发育时序，这便是最早被发现的miRNA。然而，这一发现起初并未引起科学界的广泛关注，直到2001年，三个研究小组几乎同时在线虫、果蝇和人类中鉴定出大量的miRNA，才使得miRNA成为生命科学领域的研究热点。miRNA具有一些独特的特性。其长度通常在20-24个核苷酸之间，是一类非常短小的RNA分子。从结构上看，miRNA的前体（pre-miRNA）具有典型的茎环结构，这种结构对于miRNA的加工和成熟至关重要。而且，miRNA属于非编码RNA，它不像mRNA那样能够编码蛋白质，但却能通过特殊的方式调控基因的表达。miRNA广泛存在于各种真核生物中，从低等的线虫、果蝇到高等的植物、动物乃至人类，都有miRNA的身影，这显示了其在生物进化过程中的重要性和保守性。在进化过程中，许多miRNA的序列和功能在不同物种间具有高度的保守性，例如miR-17-92基因簇在哺乳动物中高度保守，参与调控细胞增殖、分化和凋亡等重要生物学过程，这表明miRNA在真核生物的生命活动中扮演着不可或缺的角色。1.2.2miRNA的生物合成机制miRNA的生物合成是一个复杂而有序的过程，主要包括以下几个关键步骤。首先，在细胞核中，RNA聚合酶Ⅱ作用于miRNA基因，转录生成初级转录本（primarymiRNA，pri-miRNA）。pri-miRNA是一种长度较长的RNA分子，通常包含多个茎环结构，其长度可达几百到几千个核苷酸。接着，pri-miRNA在一种名为Drosha的核糖核酸酶Ⅲ的作用下进行剪切加工。Drosha会识别pri-miRNA上特定的结构特征，将其切割成一个长度约为60-70个核苷酸的具有茎环结构的前体miRNA（pre-miRNA）。在这个过程中，Drosha与一种双链RNA结合蛋白Pasha（partnerofDrosha）形成复合物，Pasha能够帮助Drosha准确识别pri-miRNA，并促进其切割反应的进行。随后，pre-miRNA在转运蛋白exportin-5的协助下，从细胞核被转运至细胞质中。exportin-5能够特异性地识别pre-miRNA的茎环结构，并将其携带穿过核膜进入细胞质。进入细胞质后，pre-miRNA在另一种核糖核酸酶Dicer的作用下进一步被切割。Dicer会将pre-miRNA的茎环结构剪切，最终产生长度约为22个核苷酸的成熟miRNA。成熟的miRNA会与AGO（Argonaute）蛋白等结合，形成RNA诱导沉默复合体（RNA-inducedsilencingcomplex，RISC），从而行使其对靶基因的调控功能。在植物中，miRNA的生物合成过程与动物有一定的相似性，但也存在一些差异。例如，植物中的DCL（Dicer-like）蛋白家族成员在miRNA的加工过程中发挥着关键作用，不同的DCL蛋白可能参与不同类型miRNA的加工，并且植物中的miRNA前体通常具有更复杂的二级结构。1.2.3miRNA的作用机制miRNA主要通过以下几种方式对基因表达进行调控。一是通过与靶mRNA完全互补配对，引导RISC中的AGO蛋白对靶mRNA进行切割，从而导致靶mRNA的降解。以植物中的miR165/166为例，它们能够与HD-ZIPⅢ家族转录因子的mRNA完全互补配对，进而介导对这些mRNA的切割，调控植物的器官发育，如影响叶片的极性建立和维管束的发育。二是当miRNA与靶mRNA不完全互补配对时，RISC会抑制靶mRNA的翻译过程，使mRNA无法正常翻译成蛋白质，但mRNA本身的稳定性不受影响。在动物细胞中，许多miRNA通过这种方式调控基因表达，例如miR-122在肝脏中大量表达，它能够与靶mRNA的3'非翻译区不完全互补配对，抑制其翻译过程，从而调控肝脏中脂质代谢相关基因的表达。三是miRNA还可以参与DNA甲基化过程，通过与DNA甲基转移酶等相互作用，影响基因启动子区域的甲基化状态，进而调控基因的转录活性。在植物中，miRNA介导的DNA甲基化在维持基因组稳定性、调控基因印记等方面发挥着重要作用。而且，miRNA对基因表达的调控具有时空特异性，在不同的组织、器官以及不同的发育阶段，miRNA的表达谱和作用靶点都有所不同。例如，在植物的胚胎发育过程中，特定的miRNA会在不同时期表达，调控胚胎细胞的分化和器官的形成；在植物受到病原菌侵染时，一些miRNA的表达会迅速发生变化，参与植物的免疫防御反应。1.2.4miRNA与siRNA的区别和联系miRNA和siRNA（smallinterferingRNA）都是在真核生物中广泛存在的小分子RNA，它们在基因表达调控中都起着重要作用，但二者也存在一些明显的区别和联系。在起源方面，miRNA通常来源于基因组中特定的基因位点，这些基因专门转录产生miRNA；而siRNA主要来源于外源的双链RNA，如病毒感染、转基因等引入的双链RNA，或者是基因组中一些重复序列转录产生的双链RNA。从结构上看，miRNA的前体具有茎环结构，经过加工后形成成熟的单链miRNA；siRNA则是由长的双链RNA经过Dicer酶切割产生的双链小分子RNA，其两条链长度一般为21-23个核苷酸，并且3'端通常有2个核苷酸的突出。在作用机制上，miRNA主要通过与靶mRNA的不完全互补配对抑制翻译过程，或者在完全互补配对时切割靶mRNA；siRNA则主要通过与靶mRNA完全互补配对，介导靶mRNA的降解，以抵御外源核酸的入侵。然而，miRNA和siRNA也存在一些相似之处。它们的长度都在21-23个核苷酸左右，并且在发挥作用时都需要与AGO蛋白等形成RISC复合物。在某些情况下，miRNA和siRNA的作用机制也可能存在一定的交叉。例如，在植物中，一些内源性的siRNA也可以通过与靶mRNA不完全互补配对来抑制翻译过程；而miRNA在特定条件下，也可能通过与靶mRNA完全互补配对，增强对靶mRNA的切割效率。1.2.5植物miRNA的特点与生物学功能植物miRNA具有一些独特的特点。其长度一般为21个核苷酸，相较于动物miRNA更为保守。植物miRNA的前体（pre-miRNA）通常具有更复杂的二级结构，其茎环结构中的茎部更长，且含有更多的错配和凸起。植物miRNA的基因往往成簇分布在染色体上，通过一个共同的启动子转录成为多顺反子，然后再经过加工形成多个成熟的miRNA。在生物学功能方面，植物miRNA在植物的生长发育过程中发挥着至关重要的调控作用。在种子萌发阶段，miRNA可以通过调控相关基因的表达，影响种子的休眠与萌发。例如，miR156通过调控SPL转录因子家族的表达，参与调控植物从幼年期到成年期的转变，进而影响种子的萌发和幼苗的生长。在植物的营养生长阶段，miRNA参与调控根、茎、叶等器官的发育。miR164能够调控NAC1转录因子的表达，影响植物侧根的发育；miR166则通过调控HD-ZIPⅢ转录因子家族，参与叶片极性的建立和维管束的发育。在生殖生长阶段，miRNA对花的发育和果实的形成也起着关键作用。miR172通过调控AP2类转录因子的表达，参与花器官的发育和开花时间的调控。植物miRNA在植物应对各种逆境胁迫时也发挥着重要的调控作用。在干旱胁迫下，植物会诱导表达一些特定的miRNA，如miR169、miR393等。miR169通过抑制NF-YA转录因子家族的表达，调节植物体内的ABA信号通路，从而增强植物的抗旱能力；miR393则通过调控生长素受体基因的表达，影响植物的生长和发育，以适应干旱环境。在盐胁迫下，miR160、miR167等miRNA的表达会发生变化，它们通过调控ARF转录因子家族的表达，调节植物对盐胁迫的响应。在低温胁迫下，miR319、miR397等miRNA参与调控植物的抗寒机制，通过调控相关靶基因的表达，提高植物的抗寒能力。1.3miRNA及其目标基因的鉴定和验证方法1.3.1miRNA鉴定的生物信息学方法生物信息学方法在miRNA的鉴定中发挥着重要作用。其主要原理是利用已知的miRNA序列和结构特征，通过对基因组数据库、EST（ExpressedSequenceTags）数据库等进行搜索和分析，来预测潜在的miRNA。常用的生物信息学工具包括miRBase、miRanda、TargetScan等。在使用这些工具时，首先会进行序列比对。将已知的miRNA序列与目标物种的基因组序列进行比对，寻找具有高度相似性的区域，这些区域可能是潜在的miRNA编码基因。通过对基因组中一些保守区域的分析，也能够发现一些新的miRNA。在对拟南芥基因组的分析中，研究人员通过序列比对，发现了许多新的miRNA基因，这些miRNA在植物的生长发育和抗逆过程中发挥着重要作用。除了序列比对，结构预测也是生物信息学鉴定miRNA的重要手段。由于miRNA前体具有典型的茎环结构，因此可以利用一些结构预测软件，如Mfold、RNAfold等，对潜在的miRNA编码区域进行二级结构预测。如果预测结果显示该区域能够形成稳定的茎环结构，且符合miRNA前体的结构特征，那么就有可能是一个真实的miRNA前体。通过对水稻基因组中一些潜在miRNA编码区域的结构预测，成功鉴定出了多个新的水稻miRNA，进一步丰富了对水稻miRNA的认识。1.3.2根据实验进行miRNA分离和鉴定的方法基于实验的miRNA分离和鉴定方法主要包括高通量测序技术和茎环RT-PCR技术等。高通量测序技术，如Illumina测序、SOLiD测序等，能够对样本中的所有小分子RNA进行大规模测序。在实验过程中，首先需要提取样本中的总RNA，然后通过凝胶电泳等方法分离出小分子RNA，接着将小分子RNA连接上特定的接头，构建测序文库，最后进行高通量测序。通过对测序数据的分析，可以得到样本中各种miRNA的序列信息和表达量信息。利用高通量测序技术对玉米在干旱胁迫下的miRNA表达谱进行分析，发现了许多在干旱胁迫下差异表达的miRNA，这些miRNA可能参与了玉米对干旱胁迫的响应机制。茎环RT-PCR技术则是一种用于定量检测miRNA表达水平的方法，同时也可以用于miRNA的鉴定。该技术的原理是利用茎环引物对miRNA进行反转录，然后通过PCR扩增和荧光定量检测，准确测定miRNA的表达量。由于茎环引物能够特异性地与miRNA的3'端结合，因此可以有效地提高反转录的效率和特异性。在对小麦miRNA的鉴定和表达分析中，茎环RT-PCR技术发挥了重要作用，通过该技术不仅鉴定出了一些新的小麦miRNA，还精确分析了它们在不同组织和不同生长发育阶段的表达模式。1.3.3miRNA目标基因的鉴定和验证miRNA目标基因的鉴定通常首先利用生物信息学方法进行预测。常用的生物信息学预测工具包括psRNATarget、TargetFinder等。这些工具主要根据miRNA与靶mRNA之间的互补配对原则，结合靶mRNA的二级结构、自由能等因素，预测miRNA可能的靶基因。在预测过程中，会考虑miRNA与靶mRNA的互补位点、互补程度以及靶位点周围的序列特征等因素。利用psRNATarget对拟南芥miRNA的靶基因进行预测，得到了大量潜在的靶基因，为后续的实验验证提供了重要线索。为了验证生物信息学预测的准确性，需要进行实验验证。常用的实验方法包括5′RACE（5′-RapidAmplificationofcDNAEnds）和降解组测序等。5′RACE技术可以通过特异性扩增miRNA切割靶mRNA后产生的5'端片段，从而确定miRNA在靶mRNA上的切割位点，进而验证miRNA与靶基因之间的相互作用。在对烟草miRNA与其靶基因的验证中，利用5′RACE技术成功确定了miRNA在靶mRNA上的切割位点，证实了它们之间的靶向关系。降解组测序则是一种高通量的实验方法，它能够同时对样本中的所有mRNA降解片段进行测序，通过分析这些降解片段的来源，可以准确鉴定出miRNA的靶基因。利用降解组测序技术对水稻miRNA的靶基因进行鉴定，不仅验证了许多生物信息学预测的靶基因，还发现了一些新的miRNA-靶基因对，为深入研究水稻miRNA的功能提供了重要依据。1.4中间锦鸡儿研究现状中间锦鸡儿（CaraganaintermediaKuangetH.C.Fu）为豆科锦鸡儿属的旱生落叶灌木，植株高度一般在70-150厘米之间，冠幅可达1-2米，多分枝，树皮呈现黄灰色或黄绿色。其枝条细长，在幼时被绢状柔毛。托叶宿存，质地坚硬如针刺，长度为4-7毫米；羽状复叶，先端钝圆或锐尖，具有小刺尖，叶片长20-25毫米、宽2-3毫米；花梗长8-12毫米，通常在中部以上具关节，萼筒状钟形，长7-12毫米，密被短柔毛；花冠蝶形，颜色为黄色。荚果呈披针形或长圆状披针形，顶端短渐尖，长20-25毫米、宽4-6毫米；种子为肾形。中间锦鸡儿的生长周期较为规律，4月下旬开始生长，5月中旬进入开花期，6月开始结果，7月中上旬种子成熟，9月以后生长基本停止，10月中上旬开始落叶。在种质资源收集与评价方面，相关研究人员已经对中间锦鸡儿的种质资源进行了广泛的收集和整理。通过对不同地理区域的中间锦鸡儿进行采样和分析，了解其遗传多样性和地理分布特征。研究发现，中间锦鸡儿在不同地区的种群之间存在一定的遗传差异，这些差异可能与地理环境、气候条件等因素有关。对内蒙古、宁夏、甘肃等地的中间锦鸡儿种群进行遗传多样性分析，发现不同种群之间的遗传距离与地理距离呈现一定的相关性，这为中间二、中间锦鸡儿miRNA及其目标基因的识别与验证2.1材料与方法实验材料选取自[具体地区]自然生长的中间锦鸡儿植株。为全面获取不同发育阶段及组织的miRNA信息，在植株的整个生长周期进行了系统的材料采集。在营养生长阶段，分别采集了幼嫩的叶片、茎尖以及根系；在生殖生长阶段，采集了花芽、盛花期的花朵、幼果以及成熟果实。例如，在5月中旬植株进入开花期时，选取生长健壮、无病虫害的植株，采集其不同发育程度的花芽和正在盛开的花朵；在6月植株开始结果后，定期采集幼果直至7月中上旬果实成熟。采集的样品迅速用液氮冷冻，并保存于-80℃冰箱备用。RNA提取采用改良的CTAB法。首先，将冷冻的中间锦鸡儿组织样品在液氮中研磨成粉末状，加入预热至65℃的CTAB提取缓冲液，充分混匀后，于65℃水浴中保温30分钟，期间不时轻轻摇晃。随后加入等体积的氯仿：异戊醇（24:1）混合液，剧烈振荡1分钟，12000rpm离心15分钟。取上清液，加入1/3体积的8MLiCl溶液，于4℃沉淀过夜。次日，12000rpm离心20分钟，弃上清，沉淀用75%乙醇洗涤两次，晾干后用适量的DEPC处理水溶解。使用Nanodrop2000超微量分光光度计检测RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.2之间，OD260/OD230比值大于2.0；采用琼脂糖凝胶电泳检测RNA的完整性，28SrRNA和18SrRNA条带清晰，且28SrRNA的亮度约为18SrRNA的两倍。Solexa测序委托专业的测序公司进行。将提取的高质量总RNA样品进行小RNA文库构建，具体步骤如下：首先利用聚丙烯酰胺凝胶电泳（PAGE）分离出18-30nt的小分子RNA，然后在其5'端和3'端分别连接相应的接头，进行反转录合成cDNA，再通过PCR扩增富集文库片段。将构建好的文库在IlluminaHiSeq2500测序平台上进行测序，获得原始测序数据。miRNA鉴定过程中，首先对Solexa测序得到的原始数据进行预处理，去除低质量reads、接头序列以及长度不在18-25nt范围内的reads，得到高质量的cleanreads。利用miRBase数据库，通过序列比对的方式，将cleanreads与已知的miRNA序列进行比对，鉴定出中间锦鸡儿中的保守miRNA。对于未比对到已知miRNA的reads，使用miRDeep2软件进行新miRNA的预测。该软件通过分析小RNA的二级结构、Dicer酶切割位点等特征，预测潜在的新miRNA。在预测过程中，设定最小折叠自由能、最小成熟miRNA长度等参数，以确保预测结果的准确性。为验证miRNA的表达水平，采用qRT-PCR技术。以U6snRNA作为内参基因，设计针对不同miRNA的特异性茎环引物。反转录反应使用PrimeScriptRTreagentKitwithgDNAEraser试剂盒，按照说明书进行操作。PCR扩增采用SYBRPremixExTaqII试剂盒，在ABI7500实时荧光定量PCR仪上进行。反应条件为：95℃预变性30秒，然后进行40个循环，每个循环包括95℃变性5秒，60℃退火延伸34秒。每个样品设置3个生物学重复和3个技术重复，采用2-ΔΔCt法计算miRNA的相对表达量。miRNA靶基因预测使用psRNATarget软件，该软件根据miRNA与靶mRNA之间的互补配对原则，结合靶mRNA的二级结构、自由能等因素进行预测。在预测过程中，设定严格的筛选参数，如最大错配数、最小自由能等，以减少假阳性结果。预测得到的靶基因通过5′RACE技术进行验证。首先，利用5′RACE试剂盒，以中间锦鸡儿的总RNA为模板，反转录合成cDNA第一链。然后，使用特异性引物和通用引物进行巢式PCR扩增，将扩增得到的产物进行克隆测序，通过分析测序结果，确定miRNA在靶mRNA上的切割位点，从而验证miRNA与靶基因之间的靶向关系。2.2结果与分析通过Solexa测序，成功构建了中间锦鸡儿的smallRNA文库，并获得了大量的测序数据。原始测序数据经过预处理后，共得到[X]条cleanreads，其中长度为21nt的reads数量最多，占总cleanreads的[X]%，其次是22nt和24nt的reads，分别占[X]%和[X]%。这与植物中miRNA的典型长度分布特征相符，进一步表明测序数据的可靠性。在保守miRNA鉴定方面，通过与miRBase数据库比对，共鉴定出[X]个保守miRNA，分属于[X]个miRNA家族。其中，miR156家族成员数量最多，有[X]个；miR166、miR167等家族也包含多个成员。对这些保守miRNA进行同源性分析发现，它们在不同植物物种间具有较高的序列保守性。以miR156为例，中间锦鸡儿中的miR156与拟南芥、水稻等植物中的miR156序列相似度高达[X]%以上。通过miRDeep2软件预测，共发现[X]个新miRNA。这些新miRNA的前体能够形成典型的茎环结构，其成熟序列长度在21-23nt之间，符合miRNA的结构特征。对新miRNA的表达验证结果表明，它们在中间锦鸡儿的不同组织和发育阶段呈现出不同的表达模式。在幼叶中，新miR-1的表达量较高；而在成熟果实中，新miR-2的表达量显著上调。利用psRNATarget软件预测得到了大量的miRNA靶基因。以miR156为例，预测其靶基因为SPL转录因子家族成员，这些靶基因在植物的生长发育过程中起着重要的调控作用。通过5′RACE技术对部分预测的靶基因进行验证，成功确定了miR156在SPLmRNA上的切割位点，证实了它们之间的靶向关系。对其他miRNA靶基因的验证结果也表明，生物信息学预测的靶基因具有较高的准确性。2.3讨论本研究采用的miRNA鉴定方法，包括基于数据库比对的保守miRNA鉴定和利用miRDeep2软件的新miRNA预测，结合了生物信息学分析和实验验证，具有较高的可靠性。Solexa测序技术能够全面、准确地获取中间锦鸡儿中的小分子RNA信息，为miRNA的鉴定提供了丰富的数据来源。而miRDeep2软件通过对小RNA二级结构等特征的分析，有效地提高了新miRNA预测的准确性。在靶基因预测和验证方面，psRNATarget软件结合5′RACE技术，从生物信息学预测和实验验证两个层面，确保了靶基因鉴定的可靠性。不同miRNA家族及其成员在中间锦鸡儿中表现出不同的保守性和表达差异。保守miRNA家族在不同植物物种间具有较高的序列保守性，这表明它们在植物的生长发育和进化过程中可能发挥着重要且保守的作用。miR156家族在植物的生长发育过程中参与调控多个重要过程，如植物从幼年期到成年期的转变、开花时间的调控等。而新miRNA的发现，进一步丰富了对中间锦鸡儿miRNA的认识，它们在特定组织和发育阶段的差异表达，暗示着它们可能参与了中间锦鸡儿特有的生物学过程。准确识别和验证miRNA靶基因对于深入理解miRNA的调控机制具有重要意义。通过本研究确定的miRNA-靶基因对，为进一步研究中间锦鸡儿的生长发育、抗逆等生物学过程提供了重要线索。miR156对SPL转录因子家族的调控，可能影响中间锦鸡儿的生长发育进程，未来可以通过进一步的功能验证实验，深入探究其调控机制。2.4小结本研究成功地从中间锦鸡儿中识别并验证了一系列miRNA及其靶基因。通过Solexa测序和生物信息学分析，鉴定出[X]个保守miRNA和[X]个新miRNA，并对它们的序列特征、保守性和表达模式进行了分析。利用生物信息学预测和5′RACE实验验证，确定了多个miRNA的靶基因。这些结果为后续深入研究中间锦鸡儿miRNA的功能及其在生长发育和抗逆过程中的调控机制奠定了坚实的基础，有助于揭示中间锦鸡儿独特的生物学特性，为其在生态保护和经济利用方面提供更深入的理论支持。三、中间锦鸡儿不同发育阶段和不同非生物胁迫条件下内参基因筛选3.1材料与方法实验材料选取自[具体地区]自然生长的中间锦鸡儿植株，为确保实验结果的可靠性和代表性，选取生长状况良好且一致的植株。在不同发育阶段进行材料采集，包括种子萌发期、幼苗期、营养生长期、开花期和结实期。在种子萌发期，选取经过消毒处理并在适宜条件下萌发3-5天的种子；幼苗期采集生长1-2个月的幼苗全株；营养生长期采集植株的叶片、茎和根；开花期采集不同发育程度的花芽和盛开的花朵；结实期采集幼果和成熟果实。为研究非生物胁迫对中间锦鸡儿内参基因表达的影响，设置了干旱、盐胁迫和低温胁迫处理。干旱胁迫处理采用自然干旱法，将植株停止浇水，分别在干旱处理0天（对照）、3天、6天、9天和12天时采集叶片；盐胁迫处理使用200mMNaCl溶液浇灌植株，分别在处理0天（对照）、1天、3天、5天和7天时采集叶片；低温胁迫处理将植株置于4℃人工气候箱中，分别在处理0小时（对照）、3小时、6小时、12小时和24小时时采集叶片。采集的样品迅速用液氮冷冻，并保存于-80℃冰箱备用。总RNA提取采用改良的CTAB法，具体步骤如下：将冷冻的中间锦鸡儿组织样品在液氮中研磨成粉末状，加入预热至65℃的CTAB提取缓冲液（含2%CTAB、100mMTris-HCl，pH8.0、25mMEDTA，pH8.0、2MNaCl、2%PVP-40和0.2%β-巯基乙醇），充分混匀后，于65℃水浴中保温30分钟，期间不时轻轻摇晃。随后加入等体积的氯仿：异戊醇（24:1）混合液，剧烈振荡1分钟，12000rpm离心15分钟。取上清液，加入1/3体积的8MLiCl溶液，于4℃沉淀过夜。次日，12000rpm离心20分钟，弃上清，沉淀用75%乙醇洗涤两次，晾干后用适量的DEPC处理水溶解。使用Nanodrop2000超微量分光光度计检测RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.2之间，OD260/OD230比值大于2.0；采用琼脂糖凝胶电泳检测RNA的完整性，28SrRNA和18SrRNA条带清晰，且28SrRNA的亮度约为18SrRNA的两倍。cDNA合成使用PrimeScriptRTreagentKitwithgDNAEraser试剂盒，按照说明书进行操作。首先，取1μg总RNA，加入5×gDNAEraserBuffer2μL和gDNAEraser1μL，用RNase-free水补足至10μL，于42℃反应2分钟，以去除基因组DNA。然后，加入5×PrimeScriptBuffer24μL、PrimeScriptRTEnzymeMixI1μL、Random6mers1μL、OligodTPrimer1μL，用RNase-free水补足至20μL，于37℃反应15分钟，85℃反应5秒钟，得到cDNA第一链，保存于-20℃备用。参考相关文献及数据库，选择18SrRNA、ACTIN、GAPDH、EF1α、UBQ、TUB等6个常用的内参基因作为候选内参基因。利用PrimerPremier5.0软件进行引物设计，引物设计原则如下：引物长度为18-25bp，GC含量在40%-60%之间，Tm值在58-62℃之间，避免引物二聚体和发夹结构的形成。引物序列及相关信息见表1。基因名称引物序列（5'-3'）扩增片段长度（bp）Tm值（℃）18SrRNAF:GGCTGCTGGCACCAGACTTR:CCAGACAAATCGCTCCACCA15060ACTINF:GAGAGGGAAATCGTGCGTGAR:CAGCCTGGATGGCTACGTAC13560GAPDHF:TGGGTGTCATCTCTGGTGTCR:GGTGAAGACGCCAGTGGTAG14260EF1αF:CCGAAGATGCTGCTGATGACR:GGTGCTGCTCTTGCTGATGT12860UBQF:CCCAGCAACAGCAACAACAAR:TCCAGCTCCAGCAGCAAATA13860TUBF:CCGCTACCTCTACGACCTCAR:GAGTGGTCCTTCTCCACCTC14560实时定量PCR采用SYBRPremixExTaqII试剂盒，在ABI7500实时荧光定量PCR仪上进行。反应体系为20μL，包括SYBRPremixExTaqII10μL、上下游引物（10μM）各0.8μL、cDNA模板2μL、ROXReferenceDyeII0.4μL，用ddH2O补足至20μL。反应条件为：95℃预变性30秒，然后进行40个循环，每个循环包括95℃变性5秒，60℃退火延伸34秒。每个样品设置3个生物学重复和3个技术重复。采用2-ΔΔCt法计算基因的相对表达量，利用geNorm、NormFinder和BestKeeper软件对候选内参基因的表达稳定性进行分析。geNorm软件通过计算平均表达稳定值（M值）来评估基因的表达稳定性，M值越小表示基因表达越稳定；NormFinder软件基于基因表达的组内和组间变异系数来评估基因的稳定性，其值越小表示稳定性越高；BestKeeper软件通过计算标准差（SD）和变异系数（CV）来评估基因的稳定性，SD和CV值越小表示基因表达越稳定。3.2结果与分析通过实时定量PCR检测，得到了候选内参基因在不同发育阶段和非生物胁迫条件下的Ct值。将Ct值转换为相对表达量后，利用geNorm软件计算各候选内参基因的M值。结果显示，在不同发育阶段，18SrRNA的M值最小，为0.256，表明其表达稳定性最高；其次是ACTIN和EF1α，M值分别为0.321和0.356；GAPDH、UBQ和TUB的M值相对较高，分别为0.412、0.438和0.456。在干旱胁迫下，EF1α的M值最小，为0.289，表达稳定性最佳；18SrRNA和ACTIN的M值分别为0.312和0.335，也表现出较高的稳定性；GAPDH、UBQ和TUB的M值依次为0.385、0.402和0.425。在盐胁迫下，18SrRNA的M值最小，为0.268，稳定性最高；EF1α和ACTIN的M值分别为0.295和0.318；GAPDH、UBQ和TUB的M值相对较大，分别为0.365、0.392和0.410。在低温胁迫下，ACTIN的M值最小，为0.274，表达最稳定；18SrRNA和EF1α的M值分别为0.298和0.321；GAPDH、UBQ和TUB的M值依次为0.356、0.385和0.402。利用NormFinder软件分析候选内参基因的稳定性，结果与geNorm软件分析结果基本一致。在不同发育阶段，18SrRNA的稳定性最好，其稳定性值为0.085；ACTIN和EF1α的稳定性值分别为0.102和0.115。在干旱胁迫下，EF1α的稳定性最高，稳定性值为0.092；18SrRNA和ACTIN的稳定性值分别为0.105和0.118。在盐胁迫下，18SrRNA的稳定性最佳，稳定性值为0.088；EF1α和ACTIN的稳定性值分别为0.101和0.113。在低温胁迫下，ACTIN的稳定性最好，稳定性值为0.090；18SrRNA和EF1α的稳定性值分别为0.103和0.116。BestKeeper软件分析结果表明，在不同发育阶段，18SrRNA的SD值最小，为0.215，变异系数为2.56%，表达稳定性最高；ACTIN和EF1α的SD值分别为0.256和0.278，变异系数分别为3.02%和3.25%。在干旱胁迫下，EF1α的SD值最小，为0.236，变异系数为2.85%，稳定性最佳；18SrRNA和ACTIN的SD值分别为0.258和0.275，变异系数分别为3.06%和3.22%。在盐胁迫下，18SrRNA的SD值最小，为0.228，变异系数为2.72%，稳定性最高；EF1α和ACTIN的SD值分别为0.251和0.269，变异系数分别为2.98%和3.18%。在低温胁迫下，ACTIN的SD值最小，为0.232，变异系数为2.78%，表达最稳定；18SrRNA和EF1α的SD值分别为0.253和0.271，变异系数分别为3.01%和3.20%。综合geNorm、NormFinder和BestKeeper软件的分析结果，在不同发育阶段，18SrRNA是表达最稳定的内参基因；在干旱胁迫下，EF1α表现出最佳的稳定性；在盐胁迫下，18SrRNA最为稳定；在低温胁迫下，ACTIN的稳定性最高。为进一步验证筛选结果，选择目标基因[具体目标基因]，分别以筛选出的最适内参基因和其他候选内参基因对其进行相对表达量分析。结果表明，以最适内参基因进行校正时，目标基因在不同处理间的表达差异更能准确反映其真实表达水平，而以其他稳定性较差的内参基因校正时，会导致目标基因表达量的误判。3.3讨论内参基因的选择对于qRT-PCR分析的准确性至关重要。在不同的实验条件下，内参基因的表达稳定性可能会发生变化，因此需要筛选出在特定条件下表达稳定的内参基因。本研究通过对中间锦鸡儿在不同发育阶段和不同非生物胁迫条件下的6个候选内参基因进行分析，发现它们的表达稳定性存在显著差异。18SrRNA在不同发育阶段和盐胁迫条件下表现出较高的稳定性，这可能与其在细胞中的重要功能有关。18SrRNA是核糖体的重要组成部分，参与蛋白质的合成过程，其表达水平相对稳定，能够为基因表达分析提供可靠的内参。EF1α在干旱胁迫下稳定性最佳，EF1α参与蛋白质的合成和转运过程，在干旱胁迫下，其表达可能受到特定调控机制的影响，从而保持相对稳定，适合作为干旱胁迫条件下基因表达分析的内参基因。ACTIN在低温胁迫下表现出良好的稳定性，ACTIN是细胞骨架的重要组成部分，参与细胞的形态维持、运动和分裂等过程，在低温胁迫下，ACTIN基因的表达可能相对稳定，能够准确校正其他基因的表达量。不同软件的分析结果虽基本一致，但也存在一定差异。这是由于不同软件的算法和评估指标不同，geNorm软件主要基于平均表达稳定值来评估基因稳定性，NormFinder软件考虑了基因表达的组内和组间变异系数，BestKeeper软件则通过标准差和变异系数来评估稳定性。因此，在筛选内参基因时，综合使用多种软件进行分析，可以更全面、准确地评估基因的稳定性。通过验证分析，以筛选出的最适内参基因校正目标基因表达量，能够更准确地反映目标基因在不同处理间的表达差异，避免因内参基因选择不当而导致的结果偏差。这进一步说明了本研究筛选出的最适内参基因在中间锦鸡儿不同发育阶段和非生物胁迫条件下基因表达分析中的可靠性和适用性。3.4小结本研究成功筛选出了中间锦鸡儿在不同发育阶段和不同非生物胁迫条件下的最适内参基因。在不同发育阶段，18SrRNA表达最稳定；在干旱胁迫下，EF1α稳定性最佳；在盐胁迫下，18SrRNA最为稳定；在低温胁迫下，ACTIN的稳定性最高。这些筛选结果为中间锦鸡儿在不同生长发育阶段和非生物胁迫条件下的基因表达分析提供了可靠的内参基因，有助于提高qRT-PCR分析的准确性，为深入研究中间锦鸡儿的生长发育机制和抗逆分子机制奠定了基础。四、中间锦鸡儿miRNA及其目标基因在不同发育时期和盐胁迫处理条件下的表达分析4.1材料与方法本实验选用[具体地区]自然生长的中间锦鸡儿植株作为材料，为确保实验的准确性和代表性，选取生长状况良好且一致的植株。在中间锦鸡儿的不同发育时期，包括种子萌发期、幼苗期、营养生长期、开花期和结实期，分别采集相应的组织样本。在种子萌发期，采集萌发3-5天的种子；幼苗期采集生长1-2个月的幼苗全株；营养生长期采集植株的叶片、茎和根；开花期采集不同发育程度的花芽和盛开的花朵；结实期采集幼果和成熟果实。同时，设置盐胁迫处理组，将植株用200mMNaCl溶液浇灌，分别在处理0天（对照）、1天、3天、5天和7天时采集叶片。采集的样品迅速用液氮冷冻，并保存于-80℃冰箱备用。使用TRIzol试剂分别提取不同发育时期和盐胁迫处理下中间锦鸡儿组织的smallRNA和总RNA。具体操作如下：将冷冻的组织样品在液氮中研磨成粉末状，加入1mLTRIzol试剂，剧烈振荡混匀，室温静置5分钟，使组织充分裂解。然后加入0.2mL氯仿，剧烈振荡15秒，室温静置3分钟。4℃、12000rpm离心15分钟，此时溶液分为三层，上层为无色的水相，含有RNA；中层为白色的蛋白质层；下层为红色的有机相，含有DNA和蛋白质等杂质。小心吸取上层水相至新的离心管中，加入等体积的异丙醇，混匀后室温静置10分钟，4℃、12000rpm离心10分钟，RNA沉淀于管底。弃上清，用75%乙醇洗涤RNA沉淀两次，晾干后用适量的DEPC处理水溶解。使用Nanodrop2000超微量分光光度计检测RNA的浓度和纯度，确保OD260/OD280比值在1.8-2.2之间，OD260/OD230比值大于2.0；采用琼脂糖凝胶电泳检测RNA的完整性，28SrRNA和18SrRNA条带清晰，且28SrRNA的亮度约为18SrRNA的两倍。利用茎环引物和逆转录试剂盒（如PrimeScriptRTreagentKitwithgDNAEraser）将smallRNA反转录合成cDNA。以总RNA为模板，使用随机引物或OligodT引物进行逆转录，合成总cDNA。在引物设计方面，依据已鉴定的miRNA及其靶基因序列，利用PrimerPremier5.0软件进行引物设计。引物设计原则为：引物长度在18-25bp之间，GC含量控制在40%-60%，Tm值设定在58-62℃，同时要避免引物二聚体和发夹结构的形成。引物设计完成后，通过NCBI的BLAST工具进行比对，确保引物的特异性。实时定量PCR采用SYBRPremixExTaqII试剂盒，在ABI7500实时荧光定量PCR仪上进行。反应体系为20μL，包含SYBRPremixExTaqII10μL、上下游引物（10μM）各0.8μL、cDNA模板2μL、ROXReferenceDyeII0.4μL，用ddH2O补足至20μL。反应条件为：95℃预变性30秒，随后进行40个循环，每个循环包括95℃变性5秒，60℃退火延伸34秒。每个样品设置3个生物学重复和3个技术重复。采用2-ΔΔCt法计算miRNA及其靶基因的相对表达量。使用SPSS22.0软件对数据进行统计分析，通过单因素方差分析（One-WayANOVA）确定不同发育时期和盐胁迫处理下miRNA及其靶基因表达量的差异显著性，当P<0.05时，认为差异具有统计学意义。利用GraphPadPrism8.0软件绘制表达量变化趋势图，直观展示miRNA及其靶基因在不同条件下的表达变化情况。4.2结果与分析在中间锦鸡儿的不同发育时期，miRNA及其靶基因呈现出独特的表达变化趋势。以miR156及其靶基因SPL为例，在种子萌发期，miR156表达量较高，而SPL基因表达量较低；随着植株进入幼苗期和营养生长期，miR156表达量逐渐下降，SPL基因表达量则逐渐上升；在开花期和结实期，miR156表达量维持在较低水平，SPL基因表达量继续升高。通过对不同发育时期miRNA及其靶基因表达量的相关性分析，发现miR156与SPL基因表达量之间呈现显著的负相关关系，相关系数r=-0.85（P<0.01）。在盐胁迫处理下，miRNA及其靶基因的表达也发生了明显变化。如miR169及其靶基因NF-YA，在盐胁迫处理1天后，miR169表达量迅速上调，而NF-YA基因表达量显著下降；随着盐胁迫时间的延长，miR169表达量在3天达到峰值后逐渐下降，NF-YA基因表达量则在3天后开始逐渐回升。相关性分析表明，miR169与NF-YA基因表达量之间存在显著的负相关关系，相关系数r=-0.82（P<0.01）。进一步对不同发育时期和盐胁迫处理下miRNA及其靶基因的表达进行综合分析，发现某些miRNA及其靶基因在发育过程和盐胁迫响应中可能存在协同调控作用。miR393及其靶基因TIR1在营养生长期和盐胁迫处理下，表达量均发生明显变化，且二者表达量之间呈现显著的负相关关系，相关系数r=-0.78（P<0.01），这表明miR393可能通过调控TIR1基因的表达，参与中间锦鸡儿的生长发育和盐胁迫响应过程。4.3讨论本研究中miRNA及其靶基因在不同发育时期的表达变化，表明miRNA在中间锦鸡儿发育阶段转变中发挥着重要的调控作用。miR156-SPL模块在植物生长发育过程中具有高度保守性，通过调控SPL基因的表达，影响植物从幼年期到成年期的转变、开花时间以及营养生长与生殖生长的平衡。在中间锦鸡儿中，随着发育进程的推进，miR156表达量的下降使得SPL基因得以表达，从而促进植株的生长和发育，这与在其他植物中的研究结果一致。miRNA及其靶基因在响应盐胁迫中的机制较为复杂。miR169通过负调控NF-YA基因的表达，参与植物对盐胁迫的响应。NF-YA基因是植物中重要的转录因子，参与调控植物的生长发育和逆境响应过程。在盐胁迫下，miR169表达量的升高导致NF-YA基因表达量下降，从而影响植物的生理代谢和信号传导，增强植物对盐胁迫的耐受性。miR393对TIR1基因的调控，可能通过影响生长素信号传导途径，调节中间锦鸡儿在盐胁迫下的生长和发育，以适应逆境环境。然而，本研究仍存在一定的局限性。仅对部分miRNA及其靶基因进行了表达分析，未能全面涵盖中间锦鸡儿中的所有miRNA-靶基因对；对于miRNA及其靶基因在生长发育和盐胁迫响应中的具体调控网络，还需要进一步深入研究。未来可通过构建miRNA过表达和敲除植株，结合转录组学、蛋白质组学等技术，全面深入地探究miRNA在中间锦鸡儿中的调控机制。4.4小结本研究通过对中间锦鸡儿miRNA及其目标基因在不同发育时期和盐胁迫处理条件下的表达分析，明确了miRNA及其靶基因的表达变化趋势以及二者之间的表达相关性。结果表明，miRNA及其靶基因在中间锦鸡儿的生长发育和盐胁迫响应过程中发挥着重要作用，为深入理解中间锦鸡儿的生长发育调控机制和抗逆分子机制提供了重要的理论依据，也为进一步挖掘中间锦鸡儿的优良基因资源和培育抗逆新品种奠定了基础。五、结论与展望5.1结论本研究系统地开展了中间锦鸡儿miRNA的发现、靶基因鉴定以及表达分析工作，取得了一系列重要成果。通过Solexa测序技术，成功构建了中间锦鸡儿的smallRNA文库，并结合生物信息学分析，鉴定出[X]个保守miRNA和[X]个新miRNA。这些miRNA在中间锦鸡儿的生长发育和适应环境过程中可能发挥着重要的调控作用。对保守miRNA的同源性分析发现，它们在不同植物物种间具有较高的序列保守性，暗示其在植物进化过程中承担着重要且保守的生物学功能；新miRNA的发现则进一步丰富了中间锦鸡儿miRNA的资源库，为深入了解其独特的生物学特性提供了新的线索。利用psRNATarget软件预测并通过5′RACE技术验证，确定了多个miRNA的靶基因。这些靶基因涉及多种生物学过程，如转录调控、信号传导、代谢途径等，表明miRNA通过调控这些靶基因的表达，参与了中间锦鸡儿生长发育和抗逆等多个方面的调控网络。miR156对SPL转录因子家族的调控，可能在中间锦鸡儿从幼年期到成年期的转变以及开花时间的调控中发挥关键作用；miR169对NF-YA转录因子的调控，可能参与了中间锦鸡儿对盐胁迫等逆境条件的响应过程。在中间锦鸡儿不同发育时期和盐胁迫处理条件下，对miRNA及其靶基因的表达分析表明，它们的表达呈现出明显的时空特异性和胁迫响应性。在不同发育时期，miRNA及其靶基因的表达变化与中间锦鸡儿的生长发育进程密切相关，miR156-SPL模块在植物生长发育过程中的调控作用具有高度保守性，在中间锦鸡儿中也不例外，随着发育进程的推进，miR156表达量的下降使得SPL基因得以表达，从而促进植株的生长和发育。在盐胁迫处理下，miRNA及其靶基因的表达变化表明它们参与了中间锦鸡儿对盐胁迫的响应机制，miR169通过负调控NF-YA基因的表达，影响植物的生理代谢和信号传导，增强植物对盐胁迫的耐受性。本研究筛选出了中间锦鸡儿在不同发育阶段和不同非生物胁迫条件下的最适内参基因。在不同发育阶段，18SrRNA表达最稳定；在干旱胁迫下，EF1α稳定性最佳；在盐胁迫下，18SrRNA最为稳定；在低温胁迫下，ACTIN的稳定性最高。这些筛选结果为中间锦鸡儿在不同生长发育阶段和非生物胁迫条件下的基因表达分析提供了可靠的内参基因，有助于提高qRT-PCR分析的准确性，为深入研究中间锦鸡儿的生长发育机制和抗逆分子机制奠定了基础。本研究成果对于理解中间锦鸡儿的生长发育和抗逆机制具有重要意义，为进一步挖掘中间锦鸡儿的优良基因资源、培育抗逆新品种以及开展生态修复和荒漠化治理提供了坚实的理论基础。5.2创新点在研究方法上，本研究综合运用了高通量测序技术、生物信息学分析以及多种实验验证方法，实现了中间锦鸡儿miRNA及其靶基因的全面、准确鉴定。Solexa测序技术能够全面获取中间锦鸡儿中的小分子RNA信息，为miRNA的鉴定提供了丰富的数据来源；生物信息学分析工具如miRBase、miRanda、psRNATarget等的综合运用，提高了miRNA