基于转录组测序技术解析牡丹响应高温胁迫的基因调控网络

基于转录组测序技术解析牡丹响应高温胁迫的基因调控网络一、引言1.1研究背景与意义近年来，全球气候变暖的趋势愈发显著，极端高温天气的出现频率和强度不断增加。据相关数据显示，过去几十年间，全球平均气温持续上升，高温灾害在全球范围内频繁发生，对自然生态系统和人类社会造成了广泛而深远的影响。植物作为生态系统的重要组成部分，在生长发育过程中不可避免地受到高温胁迫的影响。高温胁迫会干扰植物的正常生理代谢过程，如光合作用、呼吸作用、水分代谢和激素平衡等，导致植物生长发育受阻、产量降低以及品质下降。对于许多农作物而言，高温胁迫可能引发减产甚至绝收，严重威胁全球粮食安全；对于观赏植物，高温胁迫会影响其观赏价值和市场竞争力。因此，研究植物响应高温胁迫的机制，挖掘相关的耐热基因，对于培育耐高温的植物品种，提高植物对高温环境的适应能力具有重要的理论和实践意义。牡丹（PaeoniasuffruticosaAndr.）作为芍药科芍药属的落叶灌木，是中国特有的名贵木本花卉，被誉为“国色天香”“百花之王”，具有极高的观赏价值和经济价值。牡丹在中国有着悠久的栽培历史，其品种繁多，花色丰富，花型各异，深受人们的喜爱。除了观赏价值，牡丹在医药、食品、化妆品等领域也有广泛的应用。例如，牡丹根皮可入药，具有清热凉血、活血化瘀等功效；牡丹籽油富含不饱和脂肪酸，对人体健康有益；牡丹花瓣还可用于制作花茶、糕点等食品。然而，牡丹性喜夏季凉爽、冬季阳光充足、雨量适中的环境，畏炎热，耐寒冻，耐旱，怕水涝。随着全球气候变暖，高温胁迫已成为限制牡丹生长和分布的重要环境因素之一。在许多地区，夏季高温天气使得牡丹茎叶热害严重，植株提早枯萎，地下根系发育不良，导致牡丹的产量和品质降低。例如，在南方地区，夏季的高温多雨天气常常使牡丹进入强迫性休眠，影响其正常的生长发育和开花。高温还会导致牡丹的花期缩短，花朵变小，花色变淡，严重影响其观赏价值。此外，高温胁迫还会增加牡丹感染病虫害的风险，进一步危害其生长健康。目前，关于牡丹耐热性的研究相对较少，且主要集中在生理生化指标的测定和分析上，对于牡丹响应高温胁迫的分子机制研究尚处于起步阶段。深入研究牡丹响应高温胁迫的分子机制，挖掘相关的耐热基因，不仅有助于揭示牡丹耐热的分子基础，还能为培育耐热性强的牡丹新品种提供理论依据和基因资源。通过基因工程技术，将耐热基因导入牡丹中，有望提高牡丹的耐热性，扩大其种植范围，促进牡丹产业的可持续发展。因此，本研究基于转录组测序技术挖掘牡丹响应高温胁迫基因，具有重要的理论意义和实际应用价值。1.2国内外研究现状1.2.1牡丹高温胁迫生理响应研究在高温胁迫下，牡丹的生理特性会发生一系列变化。研究表明，随着温度升高，牡丹叶片的电解质渗透率显著上升，这表明高温破坏了细胞膜的完整性，导致细胞内物质外渗，膜透性增大。例如，有研究对‘洛阳红’‘凤丹白’‘映金红’三个牡丹品种进行不同温度处理，发现40℃高温胁迫处理后，三个品种叶片的细胞膜相对透性显著高于对照（25℃），其中‘洛阳红’的电解质渗透率上升幅度最大，说明其细胞膜受到的损伤较为严重，而‘凤丹白’的上升幅度相对较小，显示出相对较强的细胞膜稳定性和耐热性。游离脯氨酸作为一种重要的渗透调节物质，在牡丹应对高温胁迫时发挥着关键作用。高温胁迫下，牡丹叶片中的游离脯氨酸含量大幅增加。当温度从25℃升高到40℃时，‘映金红’叶片的游离脯氨酸含量可增加数倍，通过积累游离脯氨酸，牡丹细胞能够调节渗透压，保持细胞的水分平衡，从而缓解高温胁迫对细胞造成的损伤。同时，脯氨酸还具有清除自由基、稳定蛋白质和细胞膜结构的功能，有助于提高牡丹的耐热能力。丙二醛（MDA）是膜脂过氧化的产物，其含量变化可反映植物细胞膜脂过氧化程度和植物受到逆境胁迫的伤害程度。在高温胁迫下，牡丹叶片的MDA含量明显增加。有实验显示，在40℃高温处理下，‘洛阳红’叶片的MDA含量比25℃对照时增加了约50%，表明高温导致牡丹细胞膜发生了严重的脂过氧化作用，膜结构和功能受到损害。MDA含量的升高还会进一步加剧细胞内活性氧的积累，形成恶性循环，加重细胞的氧化损伤。过氧化物酶（POD）是植物体内重要的抗氧化酶之一，参与植物的多种生理过程，在抵御逆境胁迫中发挥着关键作用。在高温胁迫初期，牡丹叶片中的POD活性会迅速升高，以清除细胞内产生的过多活性氧，维持细胞内的氧化还原平衡。然而，随着高温胁迫时间的延长，POD活性会逐渐下降。例如，在40℃高温处理的前24小时内，‘凤丹白’叶片的POD活性显著上升，表现出较强的抗氧化能力；但处理48小时后，POD活性开始降低，这可能是由于长时间的高温胁迫导致POD酶蛋白结构受到破坏，酶活性受到抑制，使得牡丹对活性氧的清除能力下降，从而加剧了细胞的氧化损伤。此外，高温胁迫还会影响牡丹的光合作用。净光合速率（Pn）在高温下显著下降，这是由于高温破坏了光合机构，如叶绿体的结构和功能，影响了光合色素的含量和活性，导致光能吸收、传递和转化过程受阻。同时，高温还会使气孔导度降低，限制二氧化碳的进入，从而影响光合作用的碳同化过程。研究发现，当温度达到40℃时，‘映金红’的净光合速率可下降50%以上，严重影响了牡丹的生长和发育。1.2.2转录组测序技术在植物研究中的应用转录组测序技术（RNA-Seq）是一种利用高通量测序技术对转录组进行全面分析的方法，能够快速、准确地获取特定组织或细胞在某一状态下的所有转录本信息。近年来，该技术在植物研究领域得到了广泛应用，为揭示植物的生长发育、抗逆性等分子机制提供了有力手段。在植物抗逆研究方面，转录组测序技术已被用于解析多种植物对逆境胁迫的响应机制。在高温胁迫研究中，对拟南芥进行转录组测序分析，发现大量与热激蛋白、抗氧化酶、激素信号转导等相关的基因在高温胁迫下发生显著差异表达。热激蛋白基因的表达上调，能够帮助植物细胞维持蛋白质的正确折叠和结构稳定性，增强植物的耐热能力；抗氧化酶相关基因的表达变化则影响植物对活性氧的清除能力，从而减轻高温胁迫导致的氧化损伤。通过对水稻进行盐胁迫转录组分析，鉴定出一系列参与离子平衡调节、渗透调节、抗氧化防御等过程的差异表达基因，这些基因在水稻应对盐胁迫中发挥着重要作用。研究发现一些编码离子转运蛋白的基因在盐胁迫下表达上调，有助于维持细胞内的离子平衡，减轻盐分对细胞的毒害作用。在植物发育研究中，转录组测序技术也发挥了重要作用。对玉米叶枕进行单细胞转录组测序，构建了玉米叶枕的单细胞转录组图谱，发现了一系列参与玉米叶枕发育的调控基因，并对其中的部分关键基因进行了功能验证，这一研究加深了人们对玉米叶夹角建成和株型调控的理解，为作物耐密理想株型的精准设计和遗传改良提供了重要的科研基础和新视角。对拟南芥根外植体进行单细胞转录组测序，构建了愈伤组织诱导的详细单细胞转录图谱，确定了负责启动早期愈伤组织的细胞类型，重构了愈伤组织形成过程中的脱分化轨迹，并推断出调控QC-like细胞的转录因子以及与细胞命运决定相关的基因表达特征，为愈伤组织形成过程中的细胞命运转变提供了独特的视角，并提高了对愈伤组织形成及植物细胞全能性的理解。转录组测序技术还在植物的其他研究领域展现出巨大潜力。在植物次生代谢产物合成研究中，通过对丹参进行转录组测序，挖掘出参与丹参酮、丹酚酸等次生代谢产物合成的关键基因和代谢途径，为提高丹参药用成分含量、改良丹参品质提供了理论依据。在植物与病原菌互作研究中，利用转录组测序分析小麦与条锈菌互作过程中的基因表达变化，揭示了小麦响应条锈菌侵染的分子机制，为小麦抗病育种提供了重要的基因资源和理论支持。转录组测序技术为植物研究提供了全面、深入的基因表达信息，在揭示植物分子机制方面具有不可替代的重要作用。将该技术应用于牡丹响应高温胁迫基因的挖掘，有望深入解析牡丹耐热的分子基础，为牡丹耐热品种的选育提供关键的理论支持和基因资源。1.3研究目标与内容本研究旨在利用转录组测序技术，深入挖掘牡丹响应高温胁迫的关键基因，全面解析其分子机制，为牡丹耐热品种的选育提供坚实的理论基础和丰富的基因资源。具体研究内容如下：牡丹高温胁迫处理及转录组测序：挑选生长健壮、长势一致的牡丹植株，分别设置高温胁迫处理组和正常温度对照组。高温胁迫处理组将植株置于人工气候箱中，设置温度为40℃，相对湿度为70%-80%，光照强度约为150μmol/（m²・s），连续处理48h；对照组则置于25℃的环境中，其他条件与处理组相同。处理结束后，迅速采集两组牡丹植株的叶片，放入液氮中速冻，随后保存于-80℃冰箱备用。利用TRIzol法提取叶片总RNA，通过琼脂糖凝胶电泳和Nanodrop分光光度计检测RNA的质量和浓度，确保RNA的完整性和纯度符合要求。将合格的RNA样品送至专业测序公司，采用IlluminaHiSeq测序平台进行转录组测序，获得高质量的测序数据。差异表达基因筛选与分析：运用生物信息学分析软件，对测序得到的原始数据进行处理，包括去除接头序列、低质量reads和污染序列等，得到干净的reads。将干净的reads与牡丹参考基因组进行比对，统计比对率和基因表达量。以FDR（FalseDiscoveryRate）≤0.05且|log₂FC|≥1为标准，筛选出高温胁迫处理组和对照组之间的差异表达基因。对差异表达基因进行基因本体（GO）功能富集分析，从生物学过程、细胞成分和分子功能三个层面，揭示差异表达基因参与的主要生物学过程和功能；同时进行京都基因与基因组百科全书（KEGG）代谢通路富集分析，确定差异表达基因显著富集的代谢通路，明确牡丹响应高温胁迫的关键代谢途径和生物学过程。关键基因功能验证：结合GO和KEGG分析结果，挑选在牡丹响应高温胁迫过程中可能起关键作用的基因，如热激蛋白基因、抗氧化酶基因、激素信号转导相关基因等。采用实时荧光定量PCR（qRT-PCR）技术，对这些关键基因在高温胁迫不同时间点（如0h、6h、12h、24h、48h）的表达模式进行验证，确保转录组测序结果的可靠性。构建关键基因的过表达载体和RNA干扰载体，通过农杆菌介导法转化牡丹愈伤组织或原生质体，获得转基因牡丹植株。对转基因牡丹植株进行高温胁迫处理，观察其生长状况、生理指标变化以及耐热性表现，与野生型牡丹进行对比，验证关键基因在牡丹耐热性中的功能。利用酵母双杂交、双分子荧光互补等技术，研究关键基因编码蛋白与其他蛋白之间的相互作用关系，探索其参与的信号转导途径和调控网络，深入揭示牡丹响应高温胁迫的分子机制。二、转录组测序技术原理与实验设计2.1转录组测序技术原理2.1.1测序技术概述转录组测序（RNA-Seq）是一种利用高通量测序技术对细胞或组织中全部或部分mRNA、miRNA、lncRNA等转录本进行测序分析的技术。通过RNA-Seq，可以全面、准确地获取特定组织或细胞在某一状态下的转录本信息，包括基因的表达水平、转录本的结构、可变剪接事件以及新转录本的发现等。在植物研究中，转录组测序技术为揭示植物的生长发育、抗逆性、次生代谢等分子机制提供了有力手段。对于牡丹响应高温胁迫的研究，RNA-Seq技术能够从转录水平全面解析牡丹在高温胁迫下的基因表达变化，挖掘参与高温响应的关键基因和代谢途径，为深入理解牡丹耐热的分子机制提供重要的数据支持。例如，在其他植物的高温胁迫研究中，通过转录组测序发现了热激蛋白基因家族在高温响应中发挥着重要作用，这些基因能够编码热激蛋白，帮助植物细胞维持蛋白质的正确折叠和结构稳定性，增强植物的耐热能力。同时，转录组测序还揭示了激素信号转导途径、抗氧化防御系统等在植物应对高温胁迫过程中的关键作用。2.1.2技术流程本研究采用Illumina的Truseq-RNA建库方法，该方法是目前应用较为广泛的转录组文库构建方法之一，其技术流程主要包括以下几个步骤：mRNA筛选：利用真核生物mRNA尾部PolyA修饰的特性，使用带有PolyT探针的磁珠与总RNA进行杂交，使mRNA与磁珠特异性结合，从而将mRNA从总RNA中分离出来。在这个过程中，磁珠表面的PolyT探针与mRNA的PolyA尾巴通过碱基互补配对原则相互结合，形成稳定的复合物。随后，通过洗涤等操作去除未结合的RNA和杂质，最终得到纯度较高的mRNA。例如，在提取牡丹叶片总RNA后，经过与磁珠的杂交和洗脱步骤，可以有效富集mRNA，提高后续测序的准确性和有效性。cDNA合成：将筛选得到的mRNA用镁离子溶液进行处理，使其断裂成短片段。然后以这些短片段mRNA为模板，利用随机引物进行逆转录反应，合成cDNA第一链。在逆转录酶的作用下，以mRNA为模板，以dNTP为原料，按照碱基互补配对原则合成cDNA第一链。接着，加入缓冲液、dNTPs、RNaseH和DNApolymeraseI等试剂，合成cDNA第二链，形成双链cDNA。在这个过程中，RNaseH能够降解RNA-DNA杂合链中的RNA部分，为DNApolymeraseI合成第二链提供模板，从而完成双链cDNA的合成。文库构建与测序：对双链cDNA进行末端修复，使其两端平整；在末端添加碱基A，以便与测序接头连接；连接测序接头后，通过琼脂糖凝胶电泳回收目的大小的片段，并进行PCR扩增，从而完成整个文库制备工作。构建好的文库使用IlluminaHiSeq测序平台进行测序，该平台能够快速、高效地获取大量的测序数据。在PCR扩增过程中，通过引物与接头序列的特异性结合，对文库中的DNA片段进行扩增，提高文库的浓度和质量。测序过程中，文库中的DNA片段会被固定在测序芯片上，通过边合成边测序的方法，依次读取每个碱基的信息，最终得到大量的测序reads。二、转录组测序技术原理与实验设计2.2实验材料与处理2.2.1实验材料选择本研究选用牡丹品种‘洛阳红’作为实验材料。‘洛阳红’是牡丹中的传统名品，在我国广泛种植，具有极高的观赏价值，深受人们喜爱。其花色鲜艳，呈紫红色，花朵硕大，花型丰满，花期适中。该品种对高温较为敏感，在夏季高温天气下，容易出现叶片灼伤、生长发育受阻等现象，这使得它成为研究牡丹响应高温胁迫的理想材料。通过对‘洛阳红’的研究，能够更有效地揭示牡丹在高温胁迫下的生理和分子响应机制，为培育耐热性强的牡丹品种提供重要的理论依据。同时，‘洛阳红’广泛的种植范围和较高的知名度，也使得研究结果更具代表性和应用价值，便于在实际生产中推广和应用。本实验所用的‘洛阳红’牡丹植株均来自[具体来源地]，选取生长健壮、无病虫害、长势一致的三年生盆栽植株，以确保实验材料的一致性和稳定性，减少实验误差。在实验开始前，将植株置于温度为25℃、相对湿度为60%-70%、光照强度约为120μmol/（m²・s）的温室中进行适应性培养一周，使其适应实验环境，为后续实验的顺利进行奠定基础。2.2.2高温胁迫处理设置将适应性培养后的‘洛阳红’牡丹植株随机分为两组，分别进行高温胁迫处理和正常温度对照处理。高温胁迫处理组将植株置于人工气候箱中，设置温度为40℃，相对湿度为70%-80%，光照强度约为150μmol/（m²・s），连续处理48h。在高温处理过程中，每12h观察并记录植株的生长状况，包括叶片颜色、形态、是否出现灼伤等症状。正常温度对照组植株则置于温度为25℃的人工气候箱中，其他环境条件与高温处理组保持一致。为了确保实验条件的科学性和可重复性，每个处理设置三个生物学重复，每个重复包含3株牡丹植株。实验过程中，严格控制人工气候箱的各项环境参数，确保温度、湿度和光照强度的稳定性。同时，定期对人工气候箱进行检查和维护，防止设备故障对实验结果产生影响。在处理结束后，迅速采集两组牡丹植株的叶片，放入液氮中速冻，随后保存于-80℃冰箱备用，用于后续的RNA提取和转录组测序分析。2.3转录组测序实验步骤2.3.1RNA提取与质量检测从牡丹叶片中提取RNA采用TRIzol法，该方法利用TRIzol试剂中的异硫氰酸胍和苯酚等成分，能够迅速裂解细胞，使细胞内的RNA释放出来，并有效抑制RNA酶的活性，防止RNA降解。具体操作步骤如下：取适量液氮速冻后的牡丹叶片，放入研钵中，加入液氮迅速研磨成粉末状，以充分破碎细胞。将研磨好的粉末转移至含有1mlTRIzol试剂的无RNase离心管中，剧烈振荡混匀，室温静置5min，使核酸蛋白复合物完全分离。每使用1mlTRIzol试剂，加入0.2ml氯仿，盖好管盖，剧烈震荡15s，室温放置3min。随后，将离心管置于4℃、12000rpm条件下离心15min，此时样品会分成三层，上层为无色的水相，RNA主要存在于其中；中间层为白色的蛋白质层；下层为红色的有机相。小心吸取上层水相转移至新的离心管中，加入等体积的异丙醇，混匀后-20℃放置30min，使RNA沉淀。将离心管于4℃、12000rpm条件下离心10min，弃去上清液，此时可见管底有白色胶状RNA沉淀。加入1ml75%乙醇（用DEPC处理过的水配制）洗涤沉淀，4℃、10000rpm离心5min，弃去上清液，重复洗涤一次。将离心管置于室温晾干，注意不要晾得过干，以免RNA难以溶解。加入适量DEPC水，用枪头吸打几次，充分溶解RNA，得到的RNA溶液保存于-80℃冰箱备用。提取得到的RNA需进行质量检测，以确保其满足后续实验要求。首先利用Nanodrop分光光度计测定RNA的浓度和纯度，检测A260/A280比值，理想情况下该比值应在1.8-2.2之间，表明RNA纯度较高，无蛋白质和酚类等杂质污染；A260/A230比值应大于2.0，以保证RNA无多糖和盐类等杂质。随后，使用Agilent2100Bioanalyzer对RNA进行质检，通过分析18S和28SrRNA的电泳峰图来评估RNA的完整性，并计算RNA完整性指数（RIN）。一般要求RIN值在8分以上，此时RNA的完整性良好，能够满足后续文库构建和测序的需求。若RIN值低于8分，可能表明RNA存在降解，会影响测序数据的质量和准确性，需重新提取RNA。2.3.2文库构建与测序测序文库构建采用Illumina公司的TruSeqRNASamplePreparationKit试剂盒，具体步骤如下：首先，利用带有Oligo(dT)的磁珠与提取的总RNA进行杂交，根据真核生物mRNA3'端具有Poly(A)尾巴的特性，特异性地捕获mRNA，从而去除rRNA等其他RNA，提高mRNA的纯度。将捕获到的mRNA用FragmentationBuffer处理，使其片段化成为短片段，一般片段大小在200-500bp左右。以片段化后的mRNA为模板，加入六碱基随机引物（randomhexamers），在逆转录酶的作用下合成cDNA第一链。接着，加入缓冲液、dNTPs、RNaseH和DNApolymeraseI等试剂，合成cDNA第二链，形成双链cDNA。对双链cDNA进行末端修复，使其两端平整，并在3'端添加碱基A，以便与测序接头连接。连接测序接头后，通过琼脂糖凝胶电泳回收目的大小的片段，一般选择300-500bp的片段进行回收，以确保文库的质量和准确性。最后，对回收的片段进行PCR扩增，富集文库中的DNA片段，从而完成整个文库制备工作。构建好的文库使用IlluminaHiSeq系列测序平台进行测序，本研究选用IlluminaHiSeq2500测序平台，该平台具有高通量、高准确性和低成本的优势。测序策略采用双端测序（Paired-EndSequencing），测序读长为150bp。双端测序能够从DNA片段的两端同时进行测序，获得更多的序列信息，有助于后续的数据分析和基因注释。在测序过程中，文库中的DNA片段会被固定在测序芯片上，通过边合成边测序的方法，依次读取每个碱基的信息，最终得到大量的测序数据。测序数据以FASTQ格式保存，包含每个测序读段的序列信息和质量分数，为后续的生物信息学分析提供原始数据。三、测序数据分析3.1数据预处理3.1.1原始数据过滤原始测序数据中通常包含低质量reads、接头序列、污染序列等，这些数据会对后续的分析结果产生干扰，因此需要进行过滤处理以提高数据质量。在本研究中，采用了Fastp软件对原始测序数据进行过滤。Fastp是一款快速且功能强大的质量控制工具，能够高效地去除低质量碱基和接头序列。在去除低质量reads方面，设定质量值（Q）低于20的碱基比例超过20%的reads为低质量reads，将其去除。这是因为质量值Q表示碱基测序错误率，Q=-10*log10(碱基错误率)，当Q为20时，碱基错误率为1%，若低质量碱基过多，会影响后续分析的准确性。对于接头序列，Fastp通过识别和匹配已知的接头序列模式，将含有接头序列的reads去除，以避免接头序列对数据分析的干扰。此外，还去除了N碱基比例大于5%的reads，N碱基代表无法确定的碱基信息，过多的N碱基会降低数据的可靠性。经过原始数据过滤，去除了大量低质量和无用的数据，使得数据的质量得到显著提高，为后续的数据分析提供了更可靠的基础。例如，在对‘洛阳红’牡丹转录组测序的原始数据进行过滤后，有效数据的比例从原来的85%提高到了92%，这将有助于更准确地分析牡丹响应高温胁迫过程中的基因表达变化。3.1.2数据质量评估数据质量评估是转录组测序数据分析的重要环节，通过评估可以了解数据的质量状况，判断数据是否满足后续分析的要求。本研究利用FastQC工具对过滤后的干净数据进行质量评估，主要关注以下几个指标：碱基质量分布：碱基质量值反映了测序过程中碱基识别的准确性，是衡量数据质量的关键指标。FastQC生成的碱基质量分布图中，横坐标表示每条reads上碱基的位置，纵坐标表示碱基的质量值（QUAL），QUAL=-10*log10(碱基错误率)。一般来说，reads末端的碱基质量值会相对较低，但只要没有大面积的碱基质量低于20（对应碱基错误率为1%），数据质量就是可接受的。在本研究中，大部分reads的碱基质量值在整个长度范围内都保持在30以上，表明碱基识别的准确性较高，测序数据质量良好。GC含量：GC含量是指DNA或RNA中鸟嘌呤（G）和胞嘧啶（C）所占的比例。正常情况下，GC含量在不同物种和组织中相对稳定，并且与基因的功能和表达调控有关。通过分析GC含量，可以判断测序数据是否存在异常。在‘洛阳红’牡丹的转录组数据中，GC含量约为43%，与牡丹基因组的理论GC含量相符，说明测序数据没有明显的偏差或污染。测序深度：测序深度是指测序得到的总碱基数与基因组大小的比值，它反映了测序数据对基因组的覆盖程度。足够的测序深度是保证基因表达分析准确性的重要条件。在本研究中，通过计算得到测序深度约为[X]X，能够较好地覆盖牡丹基因组，确保了大多数基因能够被检测到，为后续的差异表达基因分析提供了充足的数据支持。测序错误率：测序错误率是指测序过程中错误识别碱基的比例。低测序错误率是保证数据可靠性的基础。FastQC工具通过统计分析可以估算测序错误率，在本研究中，测序错误率低于0.1%，处于较低水平，表明测序过程较为准确，数据质量可靠。除了以上指标，还对数据的重复序列比例、N碱基含量、reads长度分布等进行了评估。经过FastQC评估，各项指标均符合要求，表明过滤后的测序数据质量较高，可以用于后续的生物信息学分析，如基因表达量计算、差异表达基因筛选等，从而为深入研究牡丹响应高温胁迫的分子机制提供可靠的数据保障。三、测序数据分析3.2基因表达分析3.2.1基因表达量计算基因表达量计算是转录组数据分析的关键环节，通过量化基因表达水平，能够深入了解基因在不同条件下的功能和作用。本研究使用HTSeq软件对基因表达量进行计算，该软件基于Python语言开发，能够高效准确地将测序reads分配到相应的基因上，从而得到基因的表达量信息。HTSeq计算基因表达量的原理基于比对结果和基因注释信息。在将过滤后的干净reads与牡丹参考基因组进行比对后，HTSeq会根据基因的注释文件（如GTF格式），识别基因的外显子、内含子等特征区域。然后，它会逐行扫描比对文件，统计与每个基因特征区域重叠的reads数量。例如，对于基因A，HTSeq会计算有多少个reads能够与基因A的外显子区域完全或部分重叠，这些重叠的reads数量就代表了基因A的表达量。在统计过程中，HTSeq还考虑了测序数据的方向性和剪接情况，确保计算结果的准确性。为了更直观地量化基因表达水平，本研究采用了FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）和TPM（TranscriptsPerMillion）两种常用的标准化方法。FPKM是基于测序片段的计算方法，它考虑了基因的长度和测序深度对表达量的影响。具体计算公式为：FPKM=10^9×比对到某基因的片段数/（测序得到的总片段数×该基因的外显子长度）。通过该公式，将基因的表达量标准化到每百万测序片段中每千碱基的片段数，消除了基因长度和测序深度差异对表达量的干扰，使得不同基因和不同样本之间的表达量具有可比性。例如，基因B在样本1中的FPKM值为50，在样本2中的FPKM值为100，这表明基因B在样本2中的表达水平是样本1的两倍。TPM也是一种常用的基因表达量标准化方法，它在计算过程中首先根据每个基因的长度对其测序reads数进行校正，得到校正后的reads数。然后，将所有基因的校正后reads数总和标准化为100万，计算每个基因的TPM值。TPM值反映了每百万转录本中某基因的转录本数量，同样能够有效消除基因长度和测序深度的影响。其计算公式为：TPM=10^6×（比对到某基因的片段数/该基因的长度）/（所有基因校正后的片段数总和/10^6）。与FPKM相比，TPM在标准化过程中考虑了样本内所有基因的表达情况，使得不同样本之间的比较更加准确和稳定。通过HTSeq软件计算得到基因的原始reads计数后，利用上述标准化方法计算出FPKM和TPM值，从而获得了每个基因在高温胁迫处理组和对照组中的表达量数据。这些表达量数据为后续的差异表达基因筛选和功能分析提供了重要的基础，有助于深入揭示牡丹响应高温胁迫的分子机制。3.2.2差异表达基因筛选差异表达基因（DEGs）的筛选是转录组数据分析的核心步骤之一，通过比较高温胁迫处理组和对照组之间基因表达量的差异，能够识别出在高温胁迫响应中起关键作用的基因。本研究以|log₂foldchange|≥2且p值≤0.01为条件，利用DESeq软件进行差异表达基因的筛选。DESeq是一款专门用于分析RNA-Seq数据中差异表达基因的R语言软件包，它基于负二项分布模型，能够有效处理RNA-Seq数据中的技术重复和生物学变异，准确识别差异表达基因。其筛选差异表达基因的过程主要包括以下几个步骤：首先，DESeq对原始计数数据进行归一化处理，通过计算样本的大小因子（sizefactor），消除不同样本之间测序深度和文库制备效率的差异，使得不同样本的基因表达量具有可比性。例如，样本1和样本2的测序深度不同，通过计算大小因子，能够对样本1和样本2的基因表达量进行校正，使它们在同一水平上进行比较。接着，DESeq估计每个基因的离散度（dispersion），离散度反映了基因表达量在不同样本间的变异程度。对于每个基因，DESeq会根据样本数据估计其离散度值，并结合负二项分布模型，计算每个基因在不同条件下的表达量变化的显著性。在计算过程中，DESeq会考虑基因的生物学重复信息，提高差异表达分析的准确性和可靠性。然后，DESeq根据负二项分布模型，对每个基因在高温胁迫处理组和对照组之间的表达量进行统计检验，计算出每个基因的p值和log₂foldchange值。p值表示在原假设（即基因在两组之间无差异表达）成立的情况下，观察到当前或更极端表达量差异的概率。log₂foldchange值则表示基因在高温胁迫处理组和对照组之间表达量的变化倍数，以对数形式表示，便于直观比较基因表达量的变化趋势。例如，某基因在高温胁迫处理组中的表达量是对照组的4倍，那么其log₂foldchange值为2；若表达量是对照组的1/4，则log₂foldchange值为-2。最后，以|log₂foldchange|≥2且p值≤0.01为筛选条件，从所有基因中筛选出差异表达基因。|log₂foldchange|≥2表示基因在两组之间的表达量变化达到4倍及以上，p值≤0.01表示差异表达的显著性水平较高，这样筛选出的差异表达基因具有较高的可信度和生物学意义。通过这一筛选过程，共得到了[X]个差异表达基因，其中上调表达的基因有[X]个，下调表达的基因有[X]个。这些差异表达基因将作为后续功能分析的重点对象，通过对它们的深入研究，有望揭示牡丹响应高温胁迫的关键分子机制和信号通路。3.3功能注释与富集分析3.3.1GO功能注释将筛选得到的差异表达基因映射到GeneOntology（GO）数据库，利用Blast2GO软件进行GO功能注释，从生物过程（BiologicalProcess,BP）、细胞成分（CellularComponent,CC）和分子功能（MolecularFunction,MF）三个层面，对差异表达基因进行功能注释，以揭示其在牡丹响应高温胁迫过程中可能参与的生物学过程和功能。在生物过程方面，差异表达基因显著富集于多个与高温胁迫响应密切相关的过程。其中，“对高温的响应”过程富集了大量差异表达基因，表明这些基因在牡丹感知和应对高温环境中发挥关键作用。例如，一些热激蛋白基因在该过程中上调表达，热激蛋白能够帮助细胞维持蛋白质的正确折叠和结构稳定性，从而增强牡丹对高温胁迫的耐受性。“氧化还原过程”也显著富集，高温胁迫会导致植物细胞内活性氧（ROS）积累，引发氧化应激，而参与氧化还原过程的基因表达变化有助于调节细胞内的氧化还原平衡，清除过多的ROS，减轻氧化损伤。例如，编码超氧化物歧化酶（SOD）、过氧化物酶（POD）等抗氧化酶的基因表达上调，能够增强牡丹的抗氧化防御能力。“激素介导的信号通路”同样富集了许多差异表达基因，植物激素在植物应对逆境胁迫中发挥重要的信号传导作用，如脱落酸（ABA）、乙烯（ETH）等激素信号通路相关基因的表达变化，可能参与调控牡丹对高温胁迫的响应，调节植物的生长发育和生理代谢过程。在细胞成分层面，差异表达基因主要富集在“细胞膜”“叶绿体”“线粒体”等细胞结构相关的类别。细胞膜是细胞与外界环境的屏障，高温胁迫会影响细胞膜的流动性和稳定性，导致细胞膜损伤。富集在细胞膜相关类别的差异表达基因，可能参与维持细胞膜的完整性和功能，如编码膜转运蛋白的基因，其表达变化可能影响细胞内外物质的运输和交换，以适应高温环境。叶绿体是植物进行光合作用的场所，高温胁迫会对叶绿体的结构和功能造成损害，影响光合作用的正常进行。富集在叶绿体相关类别的基因，可能参与保护叶绿体结构，调节光合作用相关的生理过程，如一些参与光合色素合成和光系统组装的基因，其表达变化可能有助于维持牡丹在高温胁迫下的光合作用效率。线粒体是细胞的能量代谢中心，高温胁迫可能干扰线粒体的呼吸作用和能量产生。富集在线粒体相关类别的差异表达基因，可能参与调节线粒体的功能，维持细胞的能量平衡，如编码线粒体呼吸链复合物亚基的基因，其表达变化可能影响线粒体的呼吸活性，确保细胞在高温胁迫下仍能获得足够的能量。在分子功能方面，差异表达基因显著富集于“催化活性”“氧化还原酶活性”“转录因子活性”等功能类别。具有“催化活性”的基因编码各种酶类，参与细胞内的各种生化反应，如碳水化合物代谢、脂质代谢等过程中的酶基因，其表达变化可能影响牡丹在高温胁迫下的物质代谢和能量供应。“氧化还原酶活性”类别的基因在清除细胞内ROS、维持氧化还原平衡中起关键作用，如前面提到的SOD、POD等抗氧化酶基因。“转录因子活性”类别的差异表达基因编码多种转录因子，转录因子能够结合到基因的启动子区域，调控基因的转录表达。在高温胁迫下，这些转录因子可能通过调节下游基因的表达，参与牡丹的耐热调控网络，如热激转录因子（HSFs），它们可以激活热激蛋白基因等一系列耐热相关基因的表达，增强牡丹的耐热性。通过GO功能注释分析，全面了解了差异表达基因在牡丹响应高温胁迫过程中的功能分布，为深入研究牡丹耐热的分子机制提供了重要线索。这些功能注释结果也为后续的基因功能验证和调控网络研究奠定了基础，有助于进一步揭示牡丹应对高温胁迫的分子调控机制。3.3.2KEGG代谢通路分析将差异表达基因映射到京都基因与基因组百科全书（KEGG）数据库，利用KOBAS软件进行KEGG代谢通路富集分析，以确定差异表达基因显著富集的代谢途径，明确牡丹响应高温胁迫的关键代谢途径和生物学过程。分析结果显示，差异表达基因显著富集于多个重要的代谢通路。在“碳代谢”通路中，许多参与碳水化合物合成与分解的基因表达发生显著变化。高温胁迫下，牡丹可能通过调节碳代谢来维持能量平衡和物质合成。例如，参与糖酵解途径的基因表达上调，可能加速葡萄糖的分解，为细胞提供更多的能量，以应对高温胁迫带来的能量需求增加；而参与淀粉合成的基因表达下调，可能减少淀粉的积累，将更多的碳源用于能量代谢和其他生理过程。“植物激素信号转导”通路也是显著富集的通路之一。植物激素在植物应对逆境胁迫中发挥着重要的信号传导作用。在高温胁迫下，牡丹体内的激素信号转导通路发生明显变化。脱落酸（ABA）信号通路中的相关基因表达上调，ABA作为一种重要的逆境响应激素，能够调节植物的气孔关闭、渗透调节等生理过程，增强植物的抗逆性。乙烯（ETH）信号通路相关基因的表达也发生改变，乙烯参与植物的生长发育和胁迫响应过程，可能通过调节植物的生长状态和防御反应来应对高温胁迫。“氧化磷酸化”通路中的差异表达基因也较为显著。氧化磷酸化是细胞产生能量的重要过程，高温胁迫可能影响线粒体的功能，进而干扰氧化磷酸化过程。该通路中相关基因的表达变化，可能有助于维持线粒体的正常功能，保证细胞在高温胁迫下能够持续产生足够的能量，满足植物生长和防御的需求。此外，“苯丙烷生物合成”通路也受到高温胁迫的影响。苯丙烷类物质是植物次生代谢产物的重要组成部分，具有多种生物学功能，如增强植物细胞壁的稳定性、抵御病原菌入侵等。在高温胁迫下，参与苯丙烷生物合成的基因表达上调，可能促进苯丙烷类物质的合成，增强牡丹细胞壁的强度和韧性，提高其对高温胁迫的耐受性。通过KEGG代谢通路分析，明确了牡丹响应高温胁迫过程中涉及的关键代谢途径，这些代谢途径相互关联，共同调节牡丹的生理生化过程，以适应高温环境。这为进一步研究牡丹耐热的分子机制提供了重要的代谢层面的信息，有助于揭示牡丹在高温胁迫下的代谢调控网络，为培育耐热性强的牡丹品种提供理论依据。四、牡丹响应高温胁迫关键基因挖掘4.1差异表达基因的筛选与鉴定通过严格的筛选条件，以|log₂foldchange|≥2且p值≤0.01为标准，利用DESeq软件对高温胁迫处理组和对照组的基因表达数据进行分析，共筛选出[X]个差异表达基因。其中，上调表达的基因有[X]个，下调表达的基因有[X]个。这些差异表达基因的详细信息见附表1（此处假设已制作差异表达基因列表的表格，包含基因ID、基因名称、log₂foldchange值、p值等信息）。为了更直观地展示差异表达基因的分布情况，绘制了火山图（图1）。在火山图中，横坐标表示基因在高温胁迫处理组和对照组之间的表达量变化倍数（以log₂foldchange表示），纵坐标表示差异表达的显著性水平（以-log₁₀(p-value)表示）。图中红色点代表上调表达的差异基因，蓝色点代表下调表达的差异基因，灰色点代表无显著差异表达的基因。从火山图中可以清晰地看出，差异表达基因在图中呈明显的分布趋势，表明高温胁迫对牡丹叶片基因表达产生了显著影响。图1：差异表达基因火山图进一步对差异表达基因进行聚类分析，绘制了热图（图2）。热图中每一行代表一个差异表达基因，每一列代表一个样本（包括高温胁迫处理组和对照组的生物学重复样本）。通过颜色的深浅来表示基因表达量的高低，红色表示高表达，蓝色表示低表达。聚类分析结果显示，高温胁迫处理组和对照组的样本在热图上明显分为两类，表明高温胁迫导致牡丹叶片基因表达模式发生了显著改变。同时，在高温胁迫处理组中，部分基因呈现出明显的上调或下调表达趋势，这些基因可能在牡丹响应高温胁迫过程中发挥着重要作用。图2：差异表达基因热图挑选部分显著差异表达基因进行重点分析，这些基因在牡丹响应高温胁迫过程中可能具有关键功能。例如，基因A（基因ID：XXXXXX，基因名称：XXXXXX）在高温胁迫处理组中的表达量显著上调，log₂foldchange值达到3.5，p值小于0.001。该基因编码的蛋白属于热激蛋白家族，热激蛋白在植物应对高温胁迫中发挥着重要作用，能够帮助细胞维持蛋白质的正确折叠和结构稳定性，增强植物的耐热能力。基因B（基因ID：YYYYYY，基因名称：YYYYYY）在高温胁迫处理组中的表达量显著下调，log₂foldchange值为-2.8，p值小于0.01。该基因参与植物的光合作用过程，其表达下调可能导致牡丹光合作用受到抑制，进而影响植物的生长和发育。此外，还发现基因C（基因ID：ZZZZZZ，基因名称：ZZZZZZ）编码一种转录因子，在高温胁迫下表达上调。转录因子能够结合到基因的启动子区域，调控基因的转录表达，因此该基因可能通过调节下游基因的表达，参与牡丹的耐热调控网络。对这些显著差异表达基因的深入研究，将有助于揭示牡丹响应高温胁迫的分子机制，为进一步挖掘牡丹耐热关键基因提供重要线索。四、牡丹响应高温胁迫关键基因挖掘4.2关键基因的功能分析4.2.1热激蛋白基因家族热激蛋白（HeatShockProteins，HSPs）基因家族在植物应对高温胁迫过程中发挥着至关重要的作用。在本研究中，通过转录组测序分析发现多个属于PsHSP基因家族的成员在高温胁迫下呈现出显著的差异表达。对这些PsHSP基因在高温胁迫不同时间点的表达模式进行深入研究，结果显示PsHSP基因的表达量随着高温处理时间的延长呈明显的上升趋势。在高温处理初期，PsHSP基因的表达量开始逐渐增加，至24h时达到最高水平，随后表达量逐渐下降。这表明PsHSP基因能够在短期内迅速响应高温胁迫，且其表达变化与高温胁迫的时间进程密切相关。PsHSP基因家族参与牡丹耐热调控的作用机制主要体现在以下几个方面：首先，PsHSP基因编码的热激蛋白能够作为分子伴侣，帮助细胞内的蛋白质维持正确的折叠状态。在高温胁迫下，蛋白质的空间结构容易发生改变，导致其功能丧失，而热激蛋白可以与这些错误折叠的蛋白质结合，促进其重新折叠成正确的构象，从而维持细胞内蛋白质的正常功能。例如，一些小分子热激蛋白（sHSPs）能够特异性地结合到变性的蛋白质上，防止其聚集和沉淀，保持蛋白质的可溶性和活性，进而增强牡丹对高温胁迫的耐受性。其次，PsHSP基因家族可能参与调节细胞内的信号传导通路，激活一系列耐热相关基因的表达，从而增强牡丹的耐热能力。热激转录因子（HSFs）是调控热激蛋白基因表达的关键转录因子，在高温胁迫下，HSFs被激活，它们能够识别并结合到PsHSP基因启动子区域的热激元件（HSE）上，启动PsHSP基因的转录表达。同时，PsHSP基因的表达产物又可以通过反馈调节机制，对HSFs的活性进行调控，形成一个复杂的调控网络，共同参与牡丹的耐热调控过程。此外，PsHSP基因家族还可能在维持细胞膜的稳定性方面发挥作用。高温胁迫会破坏细胞膜的结构和功能，导致细胞膜透性增加，细胞内物质外流。热激蛋白可以与细胞膜上的磷脂分子相互作用，稳定细胞膜的结构，减少细胞膜的损伤，从而保证细胞的正常生理功能。研究发现，一些热激蛋白能够在细胞膜上形成一种保护屏障，阻止高温对细胞膜的直接损伤，维持细胞膜的完整性和流动性，使细胞能够在高温环境下保持正常的物质运输和信号传递功能。综上所述，PsHSP基因家族通过多种机制参与牡丹的耐热调控，在牡丹应对高温胁迫过程中发挥着不可或缺的作用。深入研究PsHSP基因家族的功能和作用机制，对于揭示牡丹耐热的分子基础具有重要意义，也为通过基因工程手段培育耐热性强的牡丹新品种提供了潜在的基因资源和理论依据。4.2.2其他关键基因除了热激蛋白基因家族外，本研究还发现了许多其他与牡丹响应高温胁迫相关的关键基因，这些基因参与了碳代谢、激素信号转导等多个重要途径，在牡丹应对高温胁迫的过程中发挥着重要作用。在碳代谢途径中，一些关键基因的表达变化对牡丹在高温胁迫下的能量供应和物质合成具有重要影响。例如，编码磷酸烯醇式丙酮酸羧化酶（PEPC）的基因在高温胁迫下表达上调。PEPC是碳代谢中的关键酶，能够催化磷酸烯醇式丙酮酸（PEP）与二氧化碳反应生成草酰乙酸，参与到C4途径和景天酸代谢途径中。在高温胁迫下，PEPC基因的表达上调可能促进了碳同化过程，增加了光合产物的合成，为细胞提供更多的能量和物质基础，以满足牡丹在高温环境下生长和防御的需求。同时，参与糖酵解和三羧酸循环的一些基因表达也发生了显著变化。糖酵解和三羧酸循环是细胞呼吸的重要途径，能够将葡萄糖等糖类物质逐步氧化分解，释放出能量，为细胞提供ATP。高温胁迫下，这些基因的表达改变可能影响了糖代谢的速率和方向，调节细胞的能量供应，以适应高温环境带来的能量需求变化。植物激素信号转导途径在牡丹响应高温胁迫过程中也起着关键的调控作用。研究发现，脱落酸（ABA）信号通路中的多个基因在高温胁迫下表达显著上调。ABA是一种重要的逆境响应激素，在植物应对高温、干旱等胁迫时发挥着核心作用。在高温胁迫下，植物体内ABA含量升高，激活ABA信号通路。例如，PYR/PYL/RCAR受体蛋白与ABA结合后，抑制PP2C蛋白磷酸酶的活性，从而解除对SnRK2蛋白激酶的抑制，激活的SnRK2进一步磷酸化下游的转录因子，如AREB/ABF等，这些转录因子结合到靶基因的启动子区域，调控一系列与抗逆相关基因的表达，包括调节气孔关闭、增强抗氧化防御系统、积累渗透调节物质等，从而提高牡丹的耐热性。乙烯（ETH）信号通路相关基因的表达也受到高温胁迫的影响。乙烯参与植物的生长发育、衰老和胁迫响应等过程。在高温胁迫下，乙烯的合成增加，激活乙烯信号通路。乙烯信号通路中的关键转录因子EIN3/EIL1能够结合到下游基因的启动子区域，调控基因表达，影响植物的生长状态和防御反应。例如，EIN3/EIL1可能调节与细胞壁合成、抗氧化酶活性等相关基因的表达，增强牡丹对高温胁迫的耐受性。此外，还发现一些参与抗氧化防御系统的基因在高温胁迫下表达上调，如编码超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）等抗氧化酶的基因。高温胁迫会导致植物细胞内活性氧（ROS）大量积累，ROS具有强氧化性，会对细胞内的生物大分子如蛋白质、核酸和脂质等造成氧化损伤，影响细胞的正常功能。抗氧化酶能够清除细胞内过多的ROS，维持细胞内的氧化还原平衡，减轻氧化损伤。SOD可以将超氧阴离子自由基歧化为过氧化氢和氧气，POD和CAT则进一步将过氧化氢分解为水和氧气，从而保护细胞免受ROS的伤害。这些与碳代谢、激素信号转导、抗氧化防御等途径相关的关键基因，通过协同作用，共同调节牡丹在高温胁迫下的生理生化过程，提高牡丹的耐热能力。对这些基因的深入研究，有助于全面揭示牡丹响应高温胁迫的分子机制，为牡丹耐热品种的选育提供更多的理论依据和基因资源。4.3基因互作网络构建利用在线数据库STRING（/）对筛选出的差异表达基因进行蛋白互作关系分析，以全面了解牡丹响应高温胁迫过程中基因之间的相互作用和调控网络。STRING数据库是一个综合性的蛋白质相互作用数据库，涵盖了超过5000个物种、2千多万种蛋白、30多亿个相互作用的信息，这些蛋白质相互作用既包括直接的物理作用，也包括共表达相关性，能够为基因互作网络的构建提供丰富的数据支持。将差异表达基因的基因名整理成列表形式，准备作为输入数据。在STRING数据库主页点击“SEARCH”，选择“Multipleproteins”执行多蛋白间的相互作用分析。输入差异表达基因列表，并在“Organism”中指定牡丹的物种信息（目前STRING数据库中虽无牡丹物种信息，但可通过选择与牡丹亲缘关系较近的物种或利用同源蛋白关系进行分析），然后点击“SEARCH”进入下一步。在新界面中，STRING数据库会返回给定基因对应的蛋白名称，默认情况下，直接点击“CONTINUE”继续。随后，数据库会匹配给定蛋白的相互作用关系，很快即可得到蛋白互作网络（PPI网络）。在得到的PPI网络中，节点代表各蛋白，节点标签为这些蛋白的名称。节点中的图案代表了该蛋白质的三维结构，如果为空则表示结构目前未知。若两个蛋白质之间存在相互作用，则以连线连接，连线的颜色反映了互作类型，包括试验验证的或预测得到的，也包含直接物理作用、共表达、基因融合等关系。有关节点和连线的颜色等所代表类型的详细信息，可点击下方“Legend”查看。例如，通过网络分析发现，热激蛋白基因（PsHSP）与多个抗氧化酶基因（如编码超氧化物歧化酶SOD、过氧化物酶POD的基因）存在相互作用关系，这表明在牡丹响应高温胁迫过程中，热激蛋白可能与抗氧化酶协同作用，共同参与维持细胞的稳态和抵抗氧化损伤。为了使网络更加简洁明了，便于分析关键基因的作用和地位，对网络图进行编辑。在“Setting”中选择“hidedisconnectednodesinthenetwork”，隐藏网络图中孤立的不存在任何互作的蛋白。同时，根据蛋白间的互作类型或互作强度、可信度等进行筛选，例如只展示直接物理互作且互作得分较高（如得分大于700，代表高可信度互作）的关系，去除低可信的、强度较弱的互作，以突出主要的基因互作关系。此外，还可以在网络中拖动节点的位置，对网络整体布局进行调整，使网络结构更加清晰直观。在构建基因互作网络后，进一步对网络中的蛋白进行基因本体（GO）富集分析，以深入了解这些蛋白所参与的生物学功能和过程。STRING数据库提供了便捷的GO富集分析功能，通过该分析可以发现，网络中的蛋白主要富集在与高温胁迫响应密切相关的生物学过程，如“蛋白质折叠”“氧化还原过程”“激素信号转导”等，进一步验证了之前GO功能注释的结果，也表明这些生物学过程中的基因之间存在紧密的相互作用和协同调控。通过对基因互作网络的分析，能够挖掘出关键基因在网络中的作用和地位。在网络中，连接度较高（即与多个其他基因存在相互作用）的基因通常被认为是关键基因，它们在网络中起到核心调控的作用。例如，某些热激转录因子（HSFs）基因在网络中处于中心位置，与多个热激蛋白基因以及其他参与高温胁迫响应的基因存在广泛的相互作用。这些热激转录因子可能通过调控下游基因的表达，激活一系列耐热相关基因的表达，从而在牡丹应对高温胁迫过程中发挥关键的调控作用。此外，一些参与激素信号转导途径的关键基因，如编码乙烯响应因子（ERF）的基因，也在网络中与多个基因存在互作关系，表明它们在整合激素信号和调控下游基因表达方面具有重要作用，进而影响牡丹对高温胁迫的响应和适应。综上所述，通过构建基因互作网络并进行分析，不仅能够直观地展示牡丹响应高温胁迫过程中基因之间的相互作用关系，还能挖掘出关键基因及其在网络中的核心作用，为深入理解牡丹响应高温胁迫的分子机制提供了重要的线索和依据，有助于进一步揭示牡丹耐热调控的复杂网络。五、实时荧光定量PCR验证5.1验证原理与方法实时荧光定量PCR（qRT-PCR）技术的基本原理是在常规PCR反应体系中加入荧光基团，通过荧光信号的变化实时监测PCR扩增过程。在PCR扩增过程中，随着扩增产物的不断积累，荧光信号强度也随之增加。每个反应管内的荧光信号达到设定的阈值时所经历的循环数被定义为Ct值（Cyclethreshold），Ct值与模板起始拷贝数的对数呈线性关系，即起始模板量越大，其达到荧光阈值所需的循环数越小，Ct值越小。利用已知起始拷贝数的标准样品绘制标准曲线，通过qPCR得到待测样本的Ct值，即可根据标准曲线计算出待测样本的起始拷贝量，从而实现对基因表达量的定量分析。根据转录组测序结果挑选出在高温胁迫下差异表达显著的基因，运用PrimerPremier6.0软件设计引物。引物设计时，遵循一系列原则以确保引物的特异性和扩增效率。引物长度一般设定在18-30bp之间，以保证引物与模板的特异性结合；上下游引物的Tm值（解链温度）控制在58-62℃之间，且二者的Tm值相差不超过2℃，以确保引物在相同的退火温度下都能有效结合模板；GC含量保持在30%-80%，避免引物中出现多个重复的碱基，特别是要避免4个或超过4个的G碱基连续出现，同时引物3’端最佳不为G或C，且3’端的5个碱基不应出现2个G或C，以减少非特异性扩增；PCR扩增产物长度控制在80-250bp之间，80-150bp最为合适，这样既能保证扩增的特异性，又能提高扩增效率。以牡丹‘洛阳红’叶片的总RNA为模板，使用TaKaRa公司的PrimeScriptRTreagentKitwithgDNAEraser试剂盒进行逆转录反应，合成cDNA第一链。逆转录反应体系包含5×PrimeScriptBuffer、PrimeScriptRTEnzymeMixI、Random6mers、OligodTPrimer、总RNA以及RNaseFreedH₂O，总体积为20μl。反应条件为：42℃孵育15min，85℃加热5s，以完成逆转录过程，得到的cDNA保存于-20℃备用。以合成的cDNA为模板进行qRT-PCR扩增，扩增体系采用20μl体系，包含10μl的2×TBGreenPremixExTaqII、0.8μl的上游引物（10μM）、0.8μl的下游引物（10μM）、2μl的cDNA模板以及6.4μl的ddH₂O。反应在ABI7500实时荧光定量PCR仪上进行，反应程序为：95℃预变性30s；95℃变性5s，60℃退火34s，共40个循环；最后进行熔解曲线分析，从60℃缓慢升温至95℃，以检测扩增产物的特异性。每个样品设置3个技术重复，同时设置无模板对照（NTC），以监测反应体系是否存在污染。采用2^-ΔΔCt法计算基因的相对表达量，以Actin基因作为内参基因，对数据进行标准化处理，分析基因在高温胁迫处理组和对照组中的表达差异。5.2验证结果与分析对挑选的[X]个差异表达基因进行qRT-PCR验证，结果显示这些基因在高温胁迫处理组和对照组中的表达趋势与转录组测序结果基本一致。以基因D（基因ID：XXXXXX，基因名称：XXXXXX）为例，转录组测序数据显示其在高温胁迫处理组中的表达量是对照组的5.6倍，log₂foldchange值为2.48；qRT-PCR结果表明，该基因在高温胁迫处理组中的相对表达量是对照组的5.2倍，两者表达倍数的变化趋势高度相似。为了更直观地展示qRT-PCR验证结果与转录组测序数据的相关性，绘制了散点图（图3）。在散点图中，横坐标表示转录组测序得到的基因表达量变化倍数（log₂foldchange），纵坐标表示qRT-PCR得到的基因相对表达量变化倍数（以2^-ΔΔCt计算）。通过计算得到两者的皮尔逊相关系数r²>0.8，表明qRT-PCR验证结果与转录组测序数据具有高度的相关性，进一步验证了转录组测序结果的可靠性。图3：qRT-PCR验证结果与转录组测序数据相关性散点图通过对验证结果的分析，还发现部分基因在qRT-PCR验证中的表达倍数变化与转录组测序数据存在一定差异。例如，基因E（基因ID：YYYYYY，基因名称：YYYYYY）在转录组测序中表达量上调3.2倍，而在qRT-PCR验证中上调2.8倍。这种差异可能是由于两种实验技术本身的差异以及实验过程中的误差导致的。转录组测序是基于高通量测序技术，能够同时检测大量基因的表达情况，但在文库构建、