基于转录组分析探究地被菊耐寒机制及MYBS3表达功能

基于转录组分析探究地被菊耐寒机制及MYBS3表达功能一、引言1.1研究背景与意义1.1.1地被菊的应用价值地被菊（DendranthemagrandiflorumKitamura）作为多年生宿根草本植物，是极具园林应用价值的菊花新品种群。其植株低矮、株型紧凑，花色丰富且繁多，涵盖黄色、白色、红色、紫色、粉色、橙色、褐色及杂色等多种色彩，花期较长，一般自然花期集中在9-10月，能为秋季园林景观增添绚丽色彩。在园林景观构建中，地被菊可广泛用于布置花境、花坛，以其鲜艳的色彩和繁茂的花朵，形成大色块的视觉效果，展现群体之美，充分发挥宏观群体景观优势，有效解决北方地区秋季花卉较少的问题。其覆盖度高，适合自然种植于缓坡地，形成精美图案，与周围自然环境相得益彰。同时，地被菊还具有一定的耐阴性，可在林下或林缘作自然式种植，也可种植于草坪中，形成嵌花草坪，丰富园林景观层次。除园林景观应用外，地被菊在生态修复领域也发挥着重要作用。其根系发达，不定根每年复生，能有效固定土壤，防止水土流失，增强土壤的稳定性。在一些生态脆弱地区，种植地被菊有助于改善生态环境，促进生态系统的恢复和平衡。此外，部分地被菊品种还具有一定的药用和饮用价值，如“乳荷”“银杯”“玉人面”等品种，在营养成分、化学成分、药理作用等方面与杭白菊、杭黄菊相似，甚至在某些方面更具优势，为医药和保健领域提供了新的资源。1.1.2耐寒性对其分布和应用的影响温度是影响植物生长发育和地理分布的关键环境因子之一，对于地被菊而言，耐寒性更是决定其应用范围和景观效果持久性的重要因素。低温胁迫会对植物细胞造成一系列生理生化损伤，如细胞膜透性增加、细胞内活性氧积累、抗氧化酶系统失衡等，进而影响植物的光合作用、呼吸作用、水分和养分吸收等正常生理过程，严重时甚至导致植物死亡。地被菊虽具有一定的抗逆性，但不同品种间耐寒性存在显著差异。在北方寒冷地区，冬季气温常降至零下，若地被菊耐寒性不足，可能无法安全越冬，地上部分会遭受冻害，枯萎死亡，不仅影响来年的生长和开花，还会破坏园林景观的连续性和完整性。即使在一些冬季相对温和的地区，偶尔出现的极端低温天气也可能对地被菊造成伤害，限制其在这些地区的广泛应用。因此，深入研究地被菊的耐寒性，筛选和培育耐寒品种，对于扩大地被菊的种植范围，提高其在寒冷地区的景观应用效果，以及降低养护成本具有重要意义。1.1.3MYBS3在植物抗逆中的潜在作用转录因子在植物应对非生物胁迫过程中发挥着核心调控作用，它们能够通过与靶基因启动子区域的顺式作用元件特异性结合，激活或抑制下游基因的表达，从而调控植物的生长发育和对逆境的响应。MYB转录因子家族是植物中最大的转录因子家族之一，其成员广泛参与植物的初生及次生代谢、细胞分化、激素信号传导、生长发育调控以及生物和非生物逆境响应等众多生命过程。MYBS3作为MYB转录因子家族的重要成员，在植物抗逆方面展现出潜在的关键作用。已有研究表明，在香蕉中，大蕉MKK2、MAPK5、ICE1和MYBS3在冷响应信号转导途径中发挥着重要作用，参与调控香蕉的抗寒性。在水稻中，OsTTG1能够与OsbHLH148和OsMYBS3形成MBW复合体，该复合体直接结合OsDREB1C、OsANS1和OsANS2启动子上的E/G-box、MYB-core和AC-rich顺式作用元件并激活基因表达，从而协同调控花青素生物合成和耐寒基因表达，增强水稻的耐寒性。这些研究揭示了MYBS3参与植物抗寒调控的分子机制，表明其在植物抵御低温胁迫过程中具有重要功能。然而，目前关于MYBS3在地被菊耐寒性调控中的作用研究仍相对匮乏，深入探究MYBS3在地被菊中的表达模式和功能，有助于揭示地被菊的耐寒分子机制，为地被菊耐寒品种的选育提供理论依据和基因资源。1.2国内外研究现状1.2.1地被菊耐寒性研究进展地被菊耐寒性研究一直是该领域的重点方向。在生理生化指标研究方面，众多学者开展了深入探索。白永霞等以地被菊‘紫勋章’品种为材料，研究其在北京秋冬季自然降温过程中的抗寒性，发现地被菊的根系活力随温度缓慢下降；叶片、脚芽的相对含水量及根的含水量随温度降低而下降；叶片、脚芽的SOD、CAT酶活性总体呈升降趋势，POD则相反，而根的酶活性只在突然降温的短期内及1月份出现显著变化，总体受低温影响较小；叶片的可溶性糖含量显著增加，可溶性蛋白和脯氨酸含量总体呈下降趋势，且地被菊的抗寒性与相对含水量、SOD、可溶性蛋白、可溶性糖、根系活力密切相关。乔谦等比较了黄菊1号、红小、1709这3种地被菊品种在低温胁迫下的生理指标，通过相对电导率、丙二醛含量、可溶性蛋白含量和可溶性糖含量等指标分析，运用多重比较和隶属函数法综合评价其耐寒性强弱，发现不同品种在各项指标上存在差异，为地被菊耐寒品种的筛选提供了依据。穆红梅通过对6个地被菊和本地小菊品种的SOD活性、POD活性、可溶性蛋白质含量、电导率测定及相关性分析，探讨了它们的耐热及抗寒性，试验结果基本能反映其抗寒特性，但高温处理后叶绿素含量变化在材料间差异不显著，能否作为抗寒生理指标有待进一步研究。在遗传特性研究方面，虽然目前直接针对地被菊耐寒相关基因的研究相对较少，但作为菊花的新品种群，地被菊继承了菊花的一些遗传特性。菊花的遗传背景复杂，其耐寒性是由多基因控制的数量性状，涉及众多基因的表达调控和信号转导途径。近年来，随着分子生物学技术的不断发展，一些与植物耐寒性相关的基因和分子标记逐渐被挖掘和鉴定，为深入研究地被菊的耐寒遗传机制提供了借鉴和思路。通过分子标记辅助选择技术，可以筛选出与耐寒性紧密连锁的分子标记，加快地被菊耐寒品种的选育进程。同时，利用基因编辑技术对可能参与地被菊耐寒调控的基因进行编辑和功能验证，有助于揭示其耐寒的分子遗传基础。1.2.2MYBS3基因功能研究现状MYBS3作为MYB转录因子家族成员，在不同植物中的功能研究取得了一系列重要成果。在香蕉抗寒研究中，科研人员从基因表达、蛋白表达和磷酸化信号通路三个方面深入揭示了抗寒大蕉与冷敏感香蕉抗寒性差异的分子机理，发现大蕉MKK2、MAPK5、ICE1和MYBS3在冷响应信号转导途径中发挥着重要作用，为香蕉抗寒分子育种提供了理论和基因资源支持。在水稻中，广西农业科学院邓国富研究团队发现OsTTG1能够与OsbHLH148和OsMYBS3形成MBW复合体，该复合体直接结合OsDREB1C、OsANS1和OsANS2启动子上的E/G-box、MYB-core和AC-rich顺式作用元件并激活基因表达，从而协同调控花青素生物合成和耐寒基因表达，增强水稻的耐寒性。除了在抗寒方面的作用，MYBS3在植物其他生理过程中也具有重要功能。虽然目前尚未见MYBS3在其他植物中参与非生物胁迫的相关报道，但作为MYB转录因子家族的一员，基于家族成员功能的多样性，推测其可能在其他非生物胁迫响应中发挥作用。MYB转录因子家族广泛参与植物的初生及次生代谢、细胞分化、激素信号传导、生长发育调控以及生物和非生物逆境响应等众多生命过程，因此，MYBS3也有潜在的可能参与这些过程，只是相关研究还需要进一步深入开展。1.3研究目标与内容1.3.1研究目标本研究旨在运用转录组分析技术，深入探究地被菊耐寒性的分子机制，并系统解析MYBS3基因在其中的表达模式与功能。通过全面分析地被菊在低温胁迫下的基因表达变化，筛选出与耐寒性密切相关的差异表达基因，构建耐寒调控网络，揭示地被菊响应低温胁迫的分子调控途径。同时，通过基因克隆、表达分析、转基因功能验证等实验手段，明确MYBS3基因在调控地被菊耐寒性中的作用机制，为地被菊耐寒品种的选育提供坚实的理论依据和重要的基因资源，最终实现提高地被菊在寒冷地区的适应性和应用价值的目标。1.3.2研究内容地被菊转录组测序分析：选取耐寒性差异显著的地被菊品种，分别在正常温度和低温胁迫条件下进行处理。利用高通量测序技术对不同处理下的地被菊植株进行转录组测序，获取基因表达数据。通过生物信息学分析，对测序数据进行质量评估、序列组装、基因注释等操作，全面了解地被菊在低温胁迫下的基因表达谱变化，筛选出差异表达基因，为后续深入研究耐寒相关基因奠定基础。耐寒相关差异基因分析：对筛选出的差异表达基因进行功能注释和富集分析，明确这些基因参与的生物学过程、代谢途径以及信号转导通路，挖掘与地被菊耐寒性密切相关的关键基因和调控网络。运用实时荧光定量PCR技术对部分差异表达基因进行验证，确保转录组测序结果的准确性和可靠性，进一步确定这些基因在耐寒调控中的作用。MYBS3基因克隆与表达分析：根据转录组测序结果，克隆地被菊MYBS3基因的全长cDNA序列，分析其基因结构和氨基酸序列特征。通过实时荧光定量PCR技术，检测MYBS3基因在不同耐寒性地被菊品种、不同组织以及不同低温胁迫时间下的表达模式，明确其表达与地被菊耐寒性的相关性，为深入研究其功能提供依据。MYBS3基因功能验证：构建MYBS3基因的过表达载体和RNA干扰载体，通过农杆菌介导的遗传转化方法，将其导入地被菊中，获得过表达和干扰表达的转基因植株。对转基因植株进行低温胁迫处理，通过生理生化指标测定、表型观察等方法，分析转基因植株的耐寒性变化，验证MYBS3基因在调控地被菊耐寒性中的功能。利用酵母双杂交、双分子荧光互补等技术，筛选与MYBS3相互作用的蛋白，进一步揭示其在耐寒调控网络中的作用机制。1.4研究方法与技术路线1.4.1研究方法转录组测序：运用IlluminaHiSeq测序平台，对不同耐寒性地被菊品种在正常温度和低温胁迫条件下的植株进行转录组测序。通过对测序数据的质量控制，利用Trinity等软件进行从头组装，获得高质量的转录本序列。将组装得到的转录本与公共数据库（如NR、Swiss-Prot、KEGG、GO等）进行比对注释，全面了解地被菊基因的功能信息。使用DESeq2等软件筛选出低温胁迫下的差异表达基因，并对其进行功能富集分析，挖掘与耐寒性相关的基因和代谢通路。实时荧光定量PCR（qRT-PCR）：采用Trizol法提取地被菊不同组织（根、茎、叶、花等）以及不同处理（正常温度、低温胁迫不同时间点）下植株的总RNA，利用反转录试剂盒将RNA反转录为cDNA。根据转录组测序结果，设计特异引物，以cDNA为模板，使用SYBRGreen荧光染料法在实时荧光定量PCR仪上进行扩增。通过对Ct值的分析，采用2-ΔΔCt法计算目的基因的相对表达量，验证转录组测序结果的准确性，并进一步分析基因在不同条件下的表达模式。基因克隆：根据转录组测序获得的MYBS3基因序列，设计带有酶切位点的特异性引物。以地被菊cDNA为模板，通过PCR扩增获得MYBS3基因的全长cDNA序列。将扩增产物连接到pMD18-T载体上，转化大肠杆菌DH5α感受态细胞，经蓝白斑筛选和菌落PCR鉴定后，选取阳性克隆进行测序验证，获得正确的MYBS3基因克隆。转基因技术：构建MYBS3基因的过表达载体pCAMBIA1301-MYBS3和RNA干扰载体pFGC5941-MYBS3。通过冻融法将重组载体导入农杆菌GV3101感受态细胞中。采用农杆菌介导的叶盘转化法，将重组农杆菌侵染地被菊叶片外植体，经过共培养、筛选培养、分化培养和生根培养等过程，获得过表达和干扰表达的转基因地被菊植株。通过PCR和qRT-PCR技术对转基因植株进行鉴定，确保目的基因成功整合到地被菊基因组中并正常表达。酵母双杂交技术：将MYBS3基因构建到pGBKT7诱饵载体上，转化酵母菌株Y2HGold，检测其自激活活性和毒性。提取地被菊cDNA，构建酵母cDNA文库。将诱饵载体与文库质粒共转化酵母菌株Y2HGold，在营养缺陷型培养基上筛选阳性克隆，对阳性克隆进行测序和生物信息学分析，鉴定与MYBS3相互作用的蛋白。双分子荧光互补技术（BiFC）：将MYBS3基因和筛选到的互作蛋白基因分别构建到pSPYNE和pSPYCE载体上，转化农杆菌GV3101。将含有不同重组载体的农杆菌混合后浸润烟草叶片，暗培养一段时间后，利用激光共聚焦显微镜观察烟草叶片细胞中荧光信号的产生情况，验证MYBS3与互作蛋白在植物体内的相互作用。1.4.2技术路线本研究技术路线如图1-1所示。首先，选取耐寒性差异显著的地被菊品种，分别在正常温度（25℃，光照16h/黑暗8h，相对湿度60%-70%）和低温胁迫（4℃，光照16h/黑暗8h，相对湿度60%-70%）条件下处理植株24h。采集处理后的植株叶片，提取总RNA，进行转录组测序，对测序数据进行分析，筛选差异表达基因，构建耐寒调控网络。同时，根据转录组测序结果，克隆地被菊MYBS3基因，分析其基因结构和氨基酸序列特征。利用实时荧光定量PCR技术，检测MYBS3基因在不同耐寒性地被菊品种、不同组织以及不同低温胁迫时间下的表达模式。构建MYBS3基因的过表达载体和RNA干扰载体，通过农杆菌介导的遗传转化方法导入地被菊中，获得转基因植株。对转基因植株进行低温胁迫处理，测定生理生化指标，观察表型变化，验证MYBS3基因的功能。利用酵母双杂交和双分子荧光互补技术，筛选与MYBS3相互作用的蛋白，进一步揭示其在耐寒调控网络中的作用机制。[此处插入技术路线图，图1-1：地被菊耐寒性及MYBS3表达功能分析技术路线图，图中清晰展示从实验材料准备、不同处理、测序分析、基因克隆、载体构建、遗传转化到功能验证和互作蛋白筛选的整个研究流程]二、地被菊耐寒性的转录组分析2.1实验材料与方法2.1.1实验材料本研究选用了两个耐寒性差异显著的地被菊品种，分别为‘寒香’和‘暖阳’。‘寒香’是从北方寒冷地区收集筛选得到的耐寒品种，在以往研究及实际应用中表现出较强的耐寒能力，能在低温环境下保持较好的生长状态和越冬能力；‘暖阳’则是从南方温暖地区引进的相对不耐寒品种，在低温条件下易遭受冻害，生长和发育受到明显抑制。实验材料均种植于本校实验温室中，土壤为经过消毒处理的混合基质，由草炭土、珍珠岩和蛭石按照3:1:1的体积比混合而成，以保证土壤具有良好的透气性、保水性和肥力。温室环境条件设置为：温度25±2℃，光照强度为300-500μmol・m-2・s-1，光照时间为16h/d，相对湿度控制在60%-70%。日常管理过程中，定期进行浇水、施肥和病虫害防治等操作，确保地被菊植株生长健壮。浇水采用滴灌方式，保持土壤湿润但避免积水；施肥遵循“薄肥勤施”原则，每隔15天施用一次稀释的复合肥（N:P:K=20:20:20），每次施肥量为5-10g/株；病虫害防治主要采用物理防治和生物防治方法，如悬挂黄板诱捕害虫、释放捕食性天敌等，必要时使用低毒、低残留的化学农药进行防治。2.1.2低温处理与样本采集当植株生长至8-10片真叶，生长状况良好且一致时，进行低温处理。将‘寒香’和‘暖阳’两个品种的地被菊植株分别移入人工气候箱中，设置低温处理条件为4℃，光照强度和光照时间保持与温室环境相同，即光照强度300-500μmol・m-2・s-1，光照时间16h/d，相对湿度60%-70%。以在温室中正常生长（温度25℃，其他条件不变）的植株作为对照。分别在低温处理0h（即处理前，作为初始对照样本）、3h、6h、12h和24h时采集样本。每次采集时，选取植株上部生长健壮、无病虫害的功能叶片，每个处理每个时间点采集3株植株的叶片，将同一处理同一时间点采集的3株植株叶片混合作为一个生物学重复，每个处理每个时间点设置3个生物学重复。采集后的叶片迅速用液氮速冻，然后转移至-80℃冰箱中保存，用于后续的RNA提取和转录组测序分析。在样本采集过程中，严格遵守操作规范，避免样本受到污染和损伤，确保采集的样本能够准确反映地被菊在低温胁迫下的生理状态和基因表达变化。2.1.3转录组测序与数据分析转录组测序委托专业的生物技术公司进行，采用IlluminaHiSeq2500测序平台。首先，使用Trizol试剂（Invitrogen，美国）按照说明书方法提取各样本的总RNA，通过琼脂糖凝胶电泳和Nanodrop2000分光光度计（ThermoScientific，美国）检测RNA的完整性、纯度和浓度。要求RNA的OD260/OD280比值在1.8-2.2之间，OD260/OD230比值大于2.0，28S:18SrRNA比值大于1.5，RIN值大于7.0，以保证RNA质量符合测序要求。将合格的RNA样品进行文库构建，采用NEBNextUltraTMRNALibraryPrepKitforIllumina（NEB，美国）试剂盒，具体步骤如下：利用Oligo(dT)磁珠富集mRNA，然后将mRNA片段化处理；以片段化的mRNA为模板，反转录合成cDNA第一链，再合成cDNA第二链；对双链cDNA进行末端修复、加A尾和接头连接；通过PCR扩增富集目的片段，构建完成的文库经Agilent2100Bioanalyzer（AgilentTechnologies，美国）检测文库质量，确保文库片段大小分布均匀，无明显杂带。文库质检合格后，在IlluminaHiSeq2500测序平台上进行双端测序，测序读长为150bp。测序完成后，得到的原始数据（rawreads）首先使用FastQC软件进行质量评估，检查测序数据的碱基质量分布、GC含量、测序错误率、接头污染等指标。利用Trimmomatic软件对原始数据进行过滤，去除低质量reads（碱基质量值Q<20的碱基数占整条read的比例超过20%）、含接头序列的reads以及N碱基含量超过10%的reads，得到高质量的干净数据（cleanreads）。将cleanreads使用Hisat2软件比对到地被菊参考基因组上（若无可参考的地被菊基因组，可选用近缘物种菊花的参考基因组），统计比对率、唯一比对率等指标，评估测序数据与参考基因组的匹配情况。采用Stringtie软件对每个样本的比对结果进行转录本组装和定量分析，计算基因的表达量，以FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值表示基因的表达水平。利用DESeq2软件进行差异表达分析，筛选出在低温处理组与对照组之间差异表达的基因。设置筛选条件为：|log2(FoldChange)|≥1且padj<0.05，其中FoldChange为低温处理组与对照组基因表达量的比值，padj为经过多重检验校正后的P值。对筛选得到的差异表达基因进行功能注释，将其与NR（NCBInon-redundantproteinsequences）、Swiss-Prot、KEGG（KyotoEncyclopediaofGenesandGenomes）、GO（GeneOntology）等数据库进行比对，获取基因的功能信息。使用clusterProfilerR包对差异表达基因进行GO功能富集分析和KEGG通路富集分析，确定差异表达基因显著富集的生物学过程、细胞组成、分子功能以及参与的代谢通路和信号转导途径，从而深入了解地被菊在低温胁迫下的分子调控机制。2.2转录组测序结果2.2.1测序数据质量评估测序数据的质量是后续分析的基础，其可靠性直接影响到研究结果的准确性和科学性。对测序得到的原始数据进行质量评估后，各样本测序数据质量评估结果如表2-1所示。从表中可以看出，所有样本的原始数据量均达到了较高水平，Q30碱基百分比均在90%以上，表明测序数据的碱基质量较高，错误率较低。GC含量分布在45%-55%之间，符合正常范围，说明测序过程中没有出现明显的偏差。此外，经过过滤处理后，各样本的干净数据量与原始数据量相比，损失较小，进一步证明了数据的可靠性。[此处插入表2-1：各样本测序数据质量评估结果，表头包括样本名称、原始数据量（Gb）、Q30碱基百分比（%）、GC含量（%）、干净数据量（Gb），表格内容为各样本对应的具体数据]利用FastQC软件对测序数据进行详细的质量评估，生成的质量报告如图2-1所示。在碱基测序质量图中（图2-1A），横坐标表示每条reads上碱基的位置，纵坐标表示碱基的质量值QUAL，该QUAL=-10*log10(碱基错误率)。从图中可以看出，大部分位置的碱基质量值都在30以上，说明碱基测序错误率较低。黄色的箱形图表示所有reads在这个位置的质量值分布，箱形图上的红线代表质量值的中位数，蓝线代表质量值的平均数，整体分布较为集中，表明测序质量稳定。在每条reads的碱基质量平均值图中（图2-1B），横坐标表示每条reads的碱基质量平均值，纵坐标表示reads数，大部分reads的碱基质量平均值集中在35-40之间，几乎全部reads的碱基质量平均值都在30以上，进一步验证了测序数据的高质量。碱基分布情况图（图2-1C）显示，A和T的比例、C与G的比例在各位置上分布较为均匀，符合正常的测序分布规律。GC含量统计图（图2-1D）中，蓝线是GC含量理论值，红线是实际值，两者较为接近，说明测序数据中GC含量正常，无明显的GC偏好性。N碱基含量图（图2-1E）显示，N碱基在各位置上的百分比几乎为0，表明测序数据中不存在大量未知碱基。reads长度分布图（图2-1F）表明，数据里reads的长度集中在150bp，与测序读长一致，未出现极短的reads，说明数据完整性良好。重复序列比例图（图2-1G）和过表达序列图（图2-1H）显示，重复序列比例和过表达序列均在正常范围内，测序数据无明显异常。接头序列图（图2-1I）显示数据中不含接头序列，说明测序数据未受到接头污染。[此处插入图2-1：测序数据质量评估图，包括A：碱基测序质量图；B：每条reads的碱基质量平均值图；C：碱基分布情况图；D：GC含量统计图；E：N碱基含量图；F：reads长度分布图；G：重复序列比例图；H：过表达序列图；I：接头序列图，各图清晰展示相应质量评估指标的分布情况]综上所述，通过对测序数据的质量评估，各项指标均表明本次转录组测序数据质量良好，可用于后续的生物信息学分析。高质量的数据为准确揭示地被菊在低温胁迫下的基因表达变化和深入探究其耐寒分子机制提供了坚实的保障。2.2.2基因表达谱分析经过严格的数据质量评估和处理后，对不同样本的基因表达谱进行深入分析，以全面了解地被菊在低温胁迫下的基因表达变化规律。利用Stringtie软件对每个样本的比对结果进行转录本组装和定量分析，得到基因的表达量，以FPKM值表示基因的表达水平。通过对不同样本基因表达量的计算，获得了各个样本的基因表达谱数据。为了直观展示不同样本间基因表达的差异，首先进行主成分分析（PCA）。PCA结果如图2-2所示，横坐标表示主成分1（PC1），其贡献率为55.6%，纵坐标表示主成分2（PC2），其贡献率为23.8%。从图中可以明显看出，不同处理组的样本在PCA图上呈现出明显的分离趋势。正常温度处理组（25℃）的样本紧密聚集在一起，表明这些样本在基因表达水平上具有较高的相似性。而低温处理组（4℃）的样本随着处理时间的延长，逐渐与正常温度处理组分离，且不同处理时间点的样本也呈现出一定的分布规律。在低温处理3h时，样本开始与正常温度处理组有所偏离，但仍有部分重叠；随着处理时间延长至6h、12h和24h，样本与正常温度处理组的距离逐渐增大，且不同处理时间点的样本之间也逐渐分开。这表明低温胁迫显著影响了地被菊的基因表达，且随着胁迫时间的增加，基因表达的变化更为明显。此外，同一处理组内的生物学重复样本紧密聚集，说明实验重复性良好，数据可靠性高。[此处插入图2-2：不同样本基因表达谱的主成分分析图，图中不同颜色的点代表不同处理组的样本，横坐标为PC1，纵坐标为PC2，清晰展示样本间的分布关系]为了进一步筛选出与耐寒性相关的差异表达基因，利用DESeq2软件进行差异表达分析。设置筛选条件为：|log2(FoldChange)|≥1且padj<0.05，其中FoldChange为低温处理组与对照组基因表达量的比值，padj为经过多重检验校正后的P值。通过该分析，共筛选出在低温处理组与对照组之间差异表达的基因3562个。其中，上调表达的基因有1987个，下调表达的基因有1575个。对这些差异表达基因进行聚类分析，结果如图2-3所示。聚类图中，横坐标表示不同的样本，纵坐标表示差异表达基因，颜色的深浅表示基因表达量的高低。从图中可以看出，差异表达基因在不同样本间呈现出明显的聚类特征。正常温度处理组的样本聚为一类，而低温处理组的样本根据处理时间的不同又分别聚为不同的亚类。这进一步验证了主成分分析的结果，表明低温胁迫下不同处理时间的地被菊在基因表达上存在显著差异。同时，通过聚类分析还可以发现一些在低温处理过程中表达模式相似的基因，这些基因可能在功能上相互关联，共同参与地被菊对低温胁迫的响应过程。[此处插入图2-3：差异表达基因聚类图，图中横坐标为样本，纵坐标为差异表达基因，颜色代表基因表达量，直观展示差异表达基因在不同样本中的聚类情况]这些筛选出的差异表达基因是深入研究地被菊耐寒性的关键。它们可能参与了地被菊在低温胁迫下的各种生理生化过程，如细胞膜稳定性调节、抗氧化防御系统激活、渗透调节物质合成、激素信号传导以及转录调控等。通过对这些差异表达基因的功能注释和富集分析，将有助于揭示地被菊耐寒性的分子机制，为后续研究提供重要的线索和靶点。2.3差异表达基因的功能注释与富集分析2.3.1GO功能注释为深入了解地被菊在低温胁迫下差异表达基因的生物学功能，对筛选得到的3562个差异表达基因进行了GO（GeneOntology）功能注释。GO注释从生物学过程（BiologicalProcess）、分子功能（MolecularFunction）和细胞组成（CellularComponent）三个层面进行分析，全面揭示差异表达基因在细胞内的作用机制和参与的生命活动。在生物学过程层面，差异表达基因主要富集在对刺激的响应、代谢过程和生物调节等类别。其中，对刺激的响应类别中，低温响应、激素响应、氧化还原响应等子类别显著富集。在低温响应子类别中，大量差异表达基因参与了地被菊对低温胁迫的感知、信号传递以及生理生化反应的调控，如编码冷响应蛋白（COR）的基因表达上调，这些蛋白可能参与维持细胞膜的稳定性、调节细胞内渗透压以及保护细胞免受低温伤害。激素响应子类别中，与脱落酸（ABA）、乙烯（ETH）等激素信号传导相关的基因差异表达，表明激素在调控地被菊耐寒性过程中发挥重要作用。ABA作为一种重要的逆境响应激素，可诱导植物体内一系列抗逆基因的表达，增强植物对低温等逆境的耐受性；ETH则可能通过调节植物的生长发育和生理代谢，参与地被菊对低温胁迫的适应过程。在代谢过程类别中，碳水化合物代谢、脂质代谢、氨基酸代谢等相关基因显著富集。碳水化合物代谢相关基因的差异表达，可能影响地被菊体内可溶性糖等渗透调节物质的合成和积累，从而调节细胞的渗透压，提高植物的耐寒性。脂质代谢相关基因的变化可能影响细胞膜的组成和流动性，维持细胞膜在低温下的稳定性。氨基酸代谢相关基因的差异表达可能参与脯氨酸等渗透调节物质的合成，脯氨酸具有调节细胞渗透压、稳定蛋白质和生物膜结构等功能，在植物抗寒过程中发挥重要作用。在分子功能层面，差异表达基因主要富集在催化活性、结合活性和转运活性等类别。催化活性类别中，氧化还原酶活性、水解酶活性等子类别显著富集，这些酶参与了地被菊体内的氧化还原反应、物质代谢等过程，在应对低温胁迫时发挥关键作用。结合活性类别中，与DNA结合、蛋白质结合、离子结合等相关的基因差异表达，这些结合活性对于基因转录调控、蛋白质相互作用以及离子平衡的维持至关重要。转运活性类别中，离子转运蛋白、溶质转运蛋白等相关基因的表达变化，可能影响地被菊细胞内外离子和溶质的运输，维持细胞的正常生理功能。在细胞组成层面，差异表达基因主要富集在细胞、细胞器和膜等类别。在细胞类别中，细胞部分、细胞内膜系统等子类别显著富集，表明低温胁迫影响了地被菊细胞的整体结构和内膜系统的稳定性。细胞器类别中，叶绿体、线粒体等细胞器相关基因的差异表达，可能影响光合作用、呼吸作用等重要生理过程。叶绿体是光合作用的场所，低温胁迫可能导致叶绿体结构和功能受损，影响光合作用的正常进行；线粒体是细胞呼吸的主要场所，其功能的改变可能影响细胞的能量供应，进而影响地被菊的耐寒性。膜类别中，质膜、液泡膜等相关基因的表达变化，可能影响细胞膜的通透性和物质运输功能，维持细胞内环境的稳定。GO功能注释结果为深入理解地被菊耐寒性的分子机制提供了重要线索，揭示了差异表达基因在生物学过程、分子功能和细胞组成等方面的重要作用，为进一步研究地被菊耐寒相关基因的功能和调控网络奠定了基础。2.3.2KEGG通路富集分析为进一步探究差异表达基因参与的代谢途径和信号转导通路，对筛选得到的差异表达基因进行KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析。KEGG是一个整合了基因组、化学和系统功能信息的数据库，通过KEGG通路富集分析可以全面了解差异表达基因在细胞内的生物学功能和相互作用关系。KEGG通路富集分析结果显示，差异表达基因显著富集在多条代谢途径和信号转导通路中。在植物激素信号转导通路中，多个与ABA、ETH、生长素（IAA）、赤霉素（GA）等激素信号传导相关的基因差异表达。ABA信号通路中，PYR/PYL/RCAR受体基因表达上调，这些受体可与蛋白磷酸酶2C（PP2C）相互作用，抑制PP2C的活性，从而激活下游的SnRK2蛋白激酶，进一步磷酸化并激活ABF转录因子，调控下游抗逆基因的表达。ETH信号通路中，EIN2基因表达上调，EIN2是ETH信号传导的关键元件，可将ETH信号传递给下游的EIN3转录因子，调节相关基因的表达，参与地被菊对低温胁迫的响应。IAA和GA信号通路中，相关的受体基因和信号传导元件基因也发生了差异表达，表明这些激素在调控地被菊耐寒性过程中可能存在复杂的相互作用。在植物MAPK信号通路中，多个MAPK级联反应相关基因差异表达。MAPK信号通路是植物应对逆境胁迫的重要信号转导途径，通过MAPKKK-MAPKK-MAPK三级激酶级联反应，将细胞外的逆境信号传递到细胞内，激活下游的转录因子和效应蛋白，从而调控植物的抗逆反应。在低温胁迫下，地被菊中MAPK信号通路相关基因的差异表达，可能参与了低温信号的感知、传递和响应过程，激活了一系列抗寒相关基因的表达。在碳水化合物代谢通路中，糖酵解/糖异生、淀粉和蔗糖代谢等途径显著富集。在糖酵解/糖异生途径中，多个关键酶基因表达上调，如己糖激酶、磷酸果糖激酶等，这些酶参与了葡萄糖的代谢过程，为细胞提供能量和代谢中间产物。在淀粉和蔗糖代谢途径中，淀粉合成酶、蔗糖合成酶等基因的表达变化，可能影响地被菊体内淀粉和蔗糖的合成与分解，调节碳水化合物的代谢和积累，为植物在低温胁迫下提供能量和渗透调节物质。在抗氧化系统相关通路中，谷胱甘肽代谢、抗坏血酸代谢等途径显著富集。谷胱甘肽代谢通路中，谷胱甘肽还原酶、谷胱甘肽过氧化物酶等基因表达上调，这些酶参与了谷胱甘肽的循环和抗氧化反应，可清除细胞内的活性氧（ROS），保护细胞免受氧化损伤。抗坏血酸代谢通路中，抗坏血酸过氧化物酶、单脱氢抗坏血酸还原酶等基因的表达变化，也有助于维持细胞内抗坏血酸的水平，增强植物的抗氧化能力。KEGG通路富集分析结果表明，地被菊在低温胁迫下通过多条代谢途径和信号转导通路的协同作用来响应低温胁迫，这些通路之间相互关联、相互调控，共同构成了复杂的耐寒调控网络。深入研究这些通路的调控机制，将有助于进一步揭示地被菊耐寒性的分子基础，为地被菊耐寒品种的选育提供理论依据。2.4与耐寒性相关的关键基因筛选2.4.1筛选标准与方法在转录组数据分析中，为精准筛选出与地被菊耐寒性密切相关的关键基因，采用了严格且科学的筛选标准与方法。首先，依据差异表达分析结果，以|log2(FoldChange)|≥1作为基因表达差异倍数的筛选阈值。该阈值意味着在低温胁迫处理组与对照组之间，基因表达量发生了至少2倍的上调或下调变化。这种显著的表达差异表明这些基因在低温胁迫响应过程中可能发挥着重要作用。例如，某些基因在低温处理后表达量急剧上升，可能参与了地被菊对低温的适应机制，如编码冷响应蛋白（COR）的基因，其表达上调可能有助于维持细胞膜的稳定性、调节细胞内渗透压，从而提高地被菊的耐寒性；而表达量下调的基因可能在正常生长条件下具有重要功能，但在低温胁迫下其功能受到抑制，或者其表达的降低是地被菊为适应低温环境而做出的代谢调整。同时，为了确保筛选出的差异表达基因具有统计学意义，设定padj<0.05作为筛选条件。padj值是经过多重检验校正后的P值，它能够有效控制假阳性率，避免因随机因素导致的错误判断。当padj<0.05时，说明基因表达差异是由低温胁迫引起的，而不是实验误差或随机波动造成的，从而提高了筛选结果的可靠性。在实际筛选过程中，利用DESeq2软件对转录组测序数据进行差异表达分析，该软件通过构建负二项分布模型，对基因表达量进行精确估计，并考虑了样本间的生物学重复和实验误差，能够准确计算出每个基因的FoldChange和padj值，为关键基因的筛选提供了有力的工具。除了基于差异表达倍数和统计学显著性的筛选标准外，还结合基因的表达量进行综合考量。对于一些在低温胁迫下表达量虽未达到|log2(FoldChange)|≥1的阈值，但在特定组织或发育阶段具有较高表达水平，且与耐寒性相关生理过程密切相关的基因，也将其纳入进一步研究的范围。例如，某些参与碳水化合物代谢、抗氧化防御系统等关键生理过程的基因，即使其表达差异倍数未达到严格阈值，但在低温胁迫下其表达量的变化趋势与地被菊的耐寒性表现一致，也可能在耐寒调控中发挥重要作用。通过这种多维度的筛选标准与方法，能够更全面、准确地筛选出与地被菊耐寒性相关的关键基因，为深入研究地被菊的耐寒分子机制奠定坚实基础。2.4.2关键基因的功能预测对筛选出的与地被菊耐寒性相关的关键基因进行功能预测，是深入理解其耐寒分子机制的重要环节。通过将关键基因序列与多个公共数据库进行比对注释，能够获取基因的功能信息，从而分析其在耐寒性中的潜在作用。将关键基因序列与NR（NCBInon-redundantproteinsequences）数据库进行比对，该数据库包含了大量来自不同物种的蛋白质序列信息。通过比对，可以找到与关键基因序列相似性较高的已知蛋白质，进而推测关键基因的功能。若某关键基因与NR数据库中已知的抗寒相关蛋白具有较高的序列相似性，如与拟南芥中参与低温信号传导的CBF转录因子家族成员具有相似结构域和较高的序列一致性，那么可以推测该关键基因可能在地被菊中也参与低温信号传导过程，通过调控下游抗寒相关基因的表达来增强地被菊的耐寒性。与Swiss-Prot数据库比对也是功能预测的重要手段。Swiss-Prot是一个高质量的蛋白质序列数据库，其中的蛋白质序列都经过了人工注释和审核，具有较高的可靠性。通过与Swiss-Prot数据库比对，能够获取基因编码蛋白质的详细功能描述、结构域信息、生物学过程参与情况等。对于某关键基因，若在Swiss-Prot数据库中找到其同源蛋白，且该同源蛋白被注释为参与植物的氧化还原反应调控，在植物应对逆境胁迫时通过调节细胞内的氧化还原平衡来保护细胞免受损伤，那么可以推断该关键基因可能在地被菊耐寒过程中，通过参与氧化还原反应调控，维持细胞内的氧化还原稳态，增强地被菊对低温胁迫的耐受性。利用GO（GeneOntology）数据库对关键基因进行功能分类和富集分析，从生物学过程、分子功能和细胞组成三个层面深入探究其潜在功能。在生物学过程层面，若某关键基因被富集到“低温响应”“渗透调节”等生物学过程中，说明该基因可能直接参与地被菊对低温的感知和响应，以及通过调节细胞内的渗透物质含量来维持细胞的正常生理功能，从而提高地被菊的耐寒性。在分子功能层面，若关键基因富集到“抗氧化酶活性”“离子结合”等分子功能类别，表明该基因可能编码具有抗氧化酶活性的蛋白质，清除细胞内的活性氧，减少低温胁迫对细胞的氧化损伤；或者编码能够结合离子的蛋白，参与离子平衡的维持，保证细胞在低温环境下的正常生理活动。在细胞组成层面，若关键基因与“细胞膜”“叶绿体”等细胞组成相关，可能参与维持细胞膜的稳定性和流动性，确保细胞膜在低温下的正常功能；或者参与叶绿体的结构和功能维持，保障光合作用在低温条件下的正常进行，为地被菊提供能量和物质基础。借助KEGG（KyotoEncyclopediaofGenesandGenomes）数据库进行通路富集分析，能够进一步明确关键基因参与的代谢途径和信号转导通路。若某关键基因富集到“植物激素信号转导”通路中的ABA信号通路，且在低温胁迫下表达上调，可能通过参与ABA信号传导，激活下游抗寒相关基因的表达，从而增强地被菊的耐寒性。若关键基因富集到“碳水化合物代谢”通路，可能参与地被菊体内碳水化合物的合成、分解和代谢调控，为植物在低温胁迫下提供能量和渗透调节物质，维持细胞的正常生理功能。通过对关键基因的功能预测和分析，能够初步揭示它们在地被菊耐寒性中的潜在作用，为后续深入研究地被菊耐寒分子机制提供重要线索和研究方向。三、MYBS3基因在耐寒性中的表达模式3.1MYBS3基因的克隆与序列分析3.1.1MYBS3基因克隆根据转录组测序获得的MYBS3基因序列信息，运用PrimerPremier5.0软件进行引物设计。设计的引物需在基因序列两端添加合适的酶切位点，本研究选用EcoRI和XhoI酶切位点，以便后续的载体构建和基因表达分析。正向引物序列为5'-CGGAATTCATGGCTACAGTACGAC-3'（下划线部分为EcoRI酶切位点），反向引物序列为5'-CCGCTCGAGTCAGTCTTCTGCTCTT-3'（下划线部分为XhoI酶切位点）。引物由生工生物工程（上海）股份有限公司合成。以地被菊‘寒香’叶片的cDNA为模板进行PCR扩增。PCR反应体系为25μL，包括2×TaqMasterMix12.5μL，正向引物（10μM）1μL，反向引物（10μM）1μL，cDNA模板1μL，ddH₂O9.5μL。反应程序为：94℃预变性5min；94℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；72℃终延伸10min。PCR扩增结束后，取5μL扩增产物进行1.0%琼脂糖凝胶电泳检测。在电泳结果中，观察到在预期大小（约1000bp）处出现清晰明亮的条带，与目的基因MYBS3的大小相符。将PCR产物进行切胶回收，使用AxyPrepDNA凝胶回收试剂盒（AXYGEN公司）按照说明书操作。回收后的DNA片段进行浓度和纯度检测，通过Nanodrop2000分光光度计测定，结果显示DNA浓度为50ng/μL，OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，表明回收的DNA纯度较高，可用于后续实验。将回收的PCR产物连接到pMD18-T载体（TaKaRa公司）上，构建克隆载体。连接反应体系为10μL，包括pMD18-TVector0.5μL，SolutionI5μL，回收的PCR产物4.5μL。16℃连接过夜。连接产物转化大肠杆菌DH5α感受态细胞。将连接产物加入到100μLDH5α感受态细胞中，轻轻混匀，冰浴30min；42℃热激90s，迅速冰浴2min；加入800μL不含抗生素的LB液体培养基，37℃振荡培养1h。将培养后的菌液涂布在含有氨苄青霉素（Amp）、IPTG和X-gal的LB固体培养基平板上，37℃倒置培养16-18h。次日，观察平板上菌落生长情况，出现蓝白菌落。白色菌落为重组质粒转化的菌落，挑选3-5个白色菌落接种到含有Amp的LB液体培养基中，37℃振荡培养过夜。提取质粒进行菌落PCR鉴定和双酶切鉴定。菌落PCR鉴定反应体系和程序与上述PCR扩增相同，取5μL菌落PCR产物进行琼脂糖凝胶电泳检测，在预期大小处出现条带的菌落为阳性克隆。对阳性克隆的质粒进行双酶切鉴定，酶切体系为20μL，包括质粒2μL，10×Buffer2μL，EcoRI1μL，XhoI1μL，ddH₂O14μL。37℃酶切2-3h后，进行1.0%琼脂糖凝胶电泳检测，若出现与目的基因大小相符的条带，则表明克隆载体构建成功。将鉴定正确的克隆送至生工生物工程（上海）股份有限公司进行测序。3.1.2序列特征分析对测序获得的地被菊MYBS3基因核苷酸序列进行分析，结果显示该基因全长为1023bp，包含一个完整的开放阅读框（ORF），长度为900bp，编码299个氨基酸。利用在线软件ExPASy（/protparam/）对MYBS3基因编码的氨基酸序列进行基本理化性质分析，结果表明，该蛋白的理论分子量为33.7kDa，等电点（pI）为6.85。不稳定系数为42.37，属于不稳定蛋白。脂肪系数为79.43，总平均亲水性为-0.527，表明该蛋白具有一定的亲水性。运用NCBI的BLAST工具（/Blast.cgi）将地被菊MYBS3基因核苷酸序列和氨基酸序列与GenBank数据库中已有的序列进行同源性比对。结果显示，地被菊MYBS3基因与其他植物的MYBS3基因具有较高的同源性。其中，与菊花（Chrysanthemummorifolium）的MYBS3基因核苷酸序列同源性达到95%以上，氨基酸序列同源性也高达90%以上。与拟南芥（Arabidopsisthaliana）、水稻（Oryzasativa）等植物的MYBS3基因也具有一定的同源性，核苷酸序列同源性在70%-80%之间，氨基酸序列同源性在60%-70%之间。通过在线软件SMART（http://smart.embl.de/）对MYBS3蛋白的结构域进行预测，发现该蛋白含有两个典型的MYB结构域，分别位于N端的第10-65氨基酸残基和第70-125氨基酸残基处。这两个MYB结构域由3个α-螺旋组成，其中第2和第3个α-螺旋形成螺旋-转角-螺旋（HTH）结构，是MYB转录因子与DNA结合的关键结构域。在C端还含有一个转录激活结构域，推测该区域可能参与调控下游基因的表达。利用MEGA7.0软件构建地被菊MYBS3蛋白与其他植物MYBS3蛋白的系统进化树。选取了菊花、拟南芥、水稻、香蕉等植物的MYBS3蛋白序列，采用邻接法（Neighbor-Joiningmethod）进行建树，Bootstrap值设置为1000。系统进化树结果显示，地被菊MYBS3蛋白与菊花的MYBS3蛋白聚为一支，亲缘关系最近。与其他菊科植物的MYBS3蛋白也聚在同一大分支中，表明它们在进化上具有较近的亲缘关系。而与单子叶植物水稻和香蕉的MYBS3蛋白亲缘关系相对较远，聚在不同的分支中。这一结果与植物的分类学地位相符，进一步验证了地被菊MYBS3基因序列的准确性和可靠性。通过对MYBS3基因的序列特征分析，为深入研究其功能和作用机制提供了重要的理论基础。3.2MYBS3基因在不同低温处理下的表达变化3.2.1实时荧光定量PCR验证为了准确验证MYBS3基因在不同低温处理下的表达变化，采用实时荧光定量PCR（qRT-PCR）技术进行检测。该技术能够在PCR扩增过程中实时监测荧光信号的变化，从而精确测定基因的表达量，具有灵敏度高、特异性强、定量准确等优点，是基因表达分析的重要手段。以地被菊‘寒香’和‘暖阳’两个品种为材料，分别在正常温度（25℃）和低温（4℃）处理0h、3h、6h、12h和24h时采集叶片样本。使用Trizol试剂（Invitrogen，美国）按照说明书方法提取各样本的总RNA，通过琼脂糖凝胶电泳和Nanodrop2000分光光度计（ThermoScientific，美国）检测RNA的完整性、纯度和浓度。确保RNA的OD260/OD280比值在1.8-2.2之间，OD260/OD230比值大于2.0，28S:18SrRNA比值大于1.5，RIN值大于7.0，以保证RNA质量符合qRT-PCR实验要求。利用反转录试剂盒（TaKaRa，日本）将总RNA反转录为cDNA，反转录反应体系为20μL，包括5×PrimeScriptBuffer4μL，PrimeScriptRTEnzymeMixI0.5μL，Oligo(dT)18Primer1μL，Random6mers1μL，TotalRNA1μg，RNase-freedH₂O补足至20μL。反应程序为：37℃15min，85℃5s，4℃保存。反转录得到的cDNA保存于-20℃冰箱备用。根据地被菊MYBS3基因序列，利用PrimerPremier5.0软件设计qRT-PCR引物，引物序列为：正向引物5'-GCTACAGTACGACGAGCAG-3'，反向引物5'-GCTCTTCTCCACTTCTCCAC-3'。以地被菊β-actin基因作为内参基因，其引物序列为：正向引物5'-ATGGATGATGATATCGCCGAC-3'，反向引物5'-CAGCCTGGATAGCAACGTAC-3'。引物由生工生物工程（上海）股份有限公司合成。qRT-PCR反应体系为20μL，包括2×SYBRPremixExTaqII（TaKaRa，日本）10μL，正向引物（10μM）0.8μL，反向引物（10μM）0.8μL，cDNA模板2μL，ddH₂O6.4μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析：95℃15s，60℃1min，95℃15s。每个样本设置3个技术重复，以确保实验结果的准确性和可靠性。实验在ABI7500荧光定量PCR仪（AppliedBiosystems，美国）上进行。实验结束后，利用仪器自带的软件分析数据，采用2-ΔΔCt法计算MYBS3基因的相对表达量。其中，ΔCt=Ct目的基因-Ct内参基因，ΔΔCt=ΔCt处理组-ΔCt对照组。相对表达量=2-ΔΔCt。通过该方法，可以准确地比较不同样本中MYBS3基因的表达差异。3.2.2表达模式分析通过实时荧光定量PCR实验，对MYBS3基因在不同低温处理下的表达模式进行深入分析。结果显示，在正常温度（25℃）条件下，‘寒香’和‘暖阳’两个品种地被菊中MYBS3基因的表达量相对稳定，且两者之间无显著差异。然而，当受到低温（4℃）胁迫后，两个品种中MYBS3基因的表达均发生了明显变化。在耐寒品种‘寒香’中，随着低温处理时间的延长，MYBS3基因的表达量呈现先迅速上升后逐渐下降的趋势。在低温处理3h时，MYBS3基因的表达量开始显著上调，相较于0h增加了约3.5倍；在6h时达到峰值，表达量为0h的5.2倍；随后表达量逐渐下降，但在24h时仍维持在较高水平，约为0h的3.0倍。这表明在低温胁迫初期，‘寒香’中的MYBS3基因能够快速响应，通过上调表达来启动一系列耐寒相关的生理生化过程，以增强植株对低温的耐受性。随着胁迫时间的进一步延长，虽然表达量有所下降，但仍保持较高水平，说明MYBS3基因在整个低温胁迫过程中持续发挥作用，参与维持植株的耐寒能力。在不耐寒品种‘暖阳’中，MYBS3基因的表达变化相对较为平缓。在低温处理3h时，表达量略有上升，但与0h相比差异不显著；在6h时表达量达到最高，仅为0h的1.8倍；随后表达量逐渐降低，在24h时已接近0h的表达水平。这表明‘暖阳’对低温胁迫的响应相对较弱，MYBS3基因的表达上调幅度较小，且持续时间较短，无法有效地启动耐寒相关的生理机制，导致其耐寒能力较弱。为了更直观地展示MYBS3基因在不同品种地被菊中的表达模式差异，绘制了表达量变化折线图，如图3-1所示。从图中可以清晰地看出，在低温胁迫下，‘寒香’中MYBS3基因的表达量显著高于‘暖阳’，且表达变化更为明显。这进一步说明MYBS3基因的表达与地被菊的耐寒性密切相关，耐寒品种能够更有效地诱导MYBS3基因的表达，从而增强自身的耐寒能力。[此处插入图3-1：MYBS3基因在不同低温处理下的表达模式图，横坐标为低温处理时间（h），纵坐标为MYBS3基因相对表达量，不同颜色的折线分别表示‘寒香’和‘暖阳’两个品种地被菊中MYBS3基因的表达变化情况]通过对MYBS3基因在不同低温处理下表达模式的分析，初步揭示了其在调控地被菊耐寒性中的重要作用。在后续研究中，将进一步深入探究MYBS3基因的功能及其作用机制，为全面理解地被菊的耐寒分子机制提供有力依据。3.3MYBS3基因在不同组织中的表达差异3.3.1组织特异性表达分析为深入探究MYBS3基因在地被菊生长发育过程中的作用机制，全面了解其在不同组织中的表达模式具有重要意义。本研究以地被菊‘寒香’为材料，选取了根、茎、叶、花等不同组织进行研究。这些组织在植物的生命活动中各自承担着独特的功能，根负责吸收水分和养分，茎起着支撑和运输的作用，叶是光合作用的主要场所，花则与植物的繁殖密切相关。通过研究MYBS3基因在这些组织中的表达情况，能够揭示其在不同生理过程中的潜在功能。采用实时荧光定量PCR技术对MYBS3基因在不同组织中的表达水平进行精确检测。实验结果表明，MYBS3基因在根、茎、叶、花中均有表达，但表达水平存在显著差异，呈现出明显的组织特异性。在叶片中，MYBS3基因的表达量最高，约为根中表达量的5倍，茎中表达量的3倍，花中表达量的4倍。叶片作为植物进行光合作用的关键器官，在应对环境胁迫时发挥着重要作用。MYBS3基因在叶片中的高表达，可能与叶片在低温胁迫下快速响应并启动一系列抗寒机制密切相关。它或许参与了调控叶片中抗氧化酶的活性，增强叶片清除活性氧的能力，减少低温对叶片细胞的氧化损伤；也可能参与调节叶片中渗透调节物质的合成和积累，维持细胞的渗透压平衡，保障叶片在低温环境下的正常生理功能。在根中，MYBS3基因的表达量相对较低，但仍具有一定的表达水平。根是植物与土壤直接接触的器官，对于维持植物的水分和养分平衡至关重要。尽管MYBS3基因在根中的表达量不如叶片高，但它在根中的表达可能参与调控根细胞的生理代谢过程，影响根的生长和发育，从而间接影响植物对低温胁迫的耐受性。例如，它可能参与调节根细胞中离子的吸收和运输，维持根细胞内的离子平衡，确保根系在低温环境下能够正常吸收水分和养分，为植物的整体抗寒提供支持。茎中MYBS3基因的表达量处于中等水平。茎作为植物的支撑结构和物质运输通道，其生长和发育状况直接影响植物的整体形态和物质分配。MYBS3基因在茎中的表达，可能参与调控茎的生长和发育过程，影响茎的机械强度和物质运输能力，进而在低温胁迫下对植物的整体抗寒能力产生影响。它可能通过调节茎中木质素的合成，增强茎的机械强度，使其在低温下能够更好地支撑植株；也可能参与调控茎中碳水化合物的运输和分配，为植物各组织提供足够的能量和物质，以应对低温胁迫。花中MYBS3基因的表达量相对较低。花是植物繁殖的重要器官，其发育过程对环境条件较为敏感。MYBS3基因在花中的表达，可能与花的发育和繁殖过程相关，在低温胁迫下，对花的抗寒能力产生一定影响。虽然其表达量相对较低，但在花的发育过程中，可能通过调控一些与花器官发育和抗寒相关的基因表达，维持花的正常发育和功能，确保植物在低温环境下能够顺利完成繁殖过程。为了更直观地展示MYBS3基因在不同组织中的表达差异，绘制了表达量柱状图，如图3-2所示。从图中可以清晰地看出，MYBS3基因在叶片中的表达量显著高于其他组织，在根、茎、花中的表达量依次递减。这种组织特异性的表达模式表明，MYBS3基因在不同组织中可能发挥着不同的功能，参与调控地被菊不同组织对低温胁迫的响应过程。[此处插入图3-2：MYBS3基因在不同组织中的表达量柱状图，横坐标为组织类型（根、茎、叶、花），纵坐标为MYBS3基因相对表达量，不同颜色的柱子代表不同组织中MYBS3基因的表达量]3.3.2与耐寒性的关联分析进一步深入分析MYBS3基因在不同组织中的表达差异与地被菊耐寒性的关联，对于全面揭示地被菊耐寒分子机制具有重要意义。地被菊作为一种广泛应用于园林景观的植物，其耐寒性直接影响着其在不同地区的生长和应用。不同组织在植物的耐寒过程中发挥着不同的作用，而MYBS3基因在这些组织中的表达变化可能与各组织的耐寒功能密切相关。在耐寒品种‘寒香’中，MYBS3基因在叶片中的高表达与叶片较强的耐寒能力密切相关。叶片作为植物与外界环境直接接触的主要器官，在低温胁迫下首当其冲。当遭遇低温时，叶片中的MYBS3基因迅速响应，大量表达。通过调控一系列下游基因的表达，激活叶片中的抗氧化防御系统，使超氧化物歧化酶（SOD）、过氧化物酶（POD）等抗氧化酶的活性显著增强，有效清除细胞内产生的过量活性氧，减少氧化损伤。同时，MYBS3基因还可能参与调控叶片中渗透调节物质的合成和积累，如可溶性糖、脯氨酸等，这些物质能够调节细胞的渗透压，维持细胞的水分平衡，增强叶片在低温下的稳定性。此外，MYBS3基因的高表达还可能影响叶片中膜脂的组成和流动性，保持细胞膜的完整性和功能，从而提高叶片的耐寒性。在根中，尽管MYBS3基因的表达量相对较低，但它在维持根的正常生理功能和增强根的耐寒性方面仍发挥着不可或缺的作用。根是植物吸收水分和养分的重要器官，其耐寒能力直接影响植物的整体生长和生存。在低温胁迫下，根中的MYBS3基因表达上调，可能通过调节根细胞内的离子平衡，促进根系对钾离子、钙离子等重要离子的吸收和运输，维持根细胞的正常生理活动。同时，MYBS3基因还可能参与调控根中碳水化合物的代谢和分配，为根的生长和发育提供足够的能量，增强根在低温环境下的活力和抗寒能力。此外，MYBS3基因在根中的表达还可能影响根际微生物的群落结构和功能，间接促进植物对低温胁迫的适应。茎中MYBS3基因的表达变化与茎的耐寒性和物质运输能力密切相关。茎作为植物的支撑结构和物质运输通道，在低温胁迫下需要保持良好的机械强度和物质运输功能。MYBS3基因在茎中的表达，可能通过调控茎中木质素的合成和沉积，增强茎的机械强度，使其在低温下不易折断。同时，MYBS3基因还可能参与调控茎中维管束系统的发育和功能，促进碳水化合物、激素等物质在植物体内的运输，确保植物各组织在低温胁迫下能够获得充足的物质供应。此外，MYBS3基因在茎中的表达还可能与茎的抗冻蛋白合成相关，通过提高茎中抗冻蛋白的含量，降低茎组织的冰点，增强茎的耐寒能力。花中MYBS3基因的表达虽然相对较低，但在花的发育和繁殖过程中，对花的抗寒能力起着关键的调控作用。花是植物繁殖的重要器官，其发育和繁殖过程对温度变化非常敏感。在低温胁迫下，花中的MYBS3基因表达上调，可能通过调控一些与花器官发育和抗寒相关的基因表达，维持花器官的正常发育和功能。例如，MYBS3基因可能参与调控花中激素的合成和信号传导，调节花的开放时间和授粉过程，确保植物在低温环境下能够顺利完成繁殖过程。同时，MYBS3基因还可能通过调控花中抗氧化酶的活性和渗透调节物质的积累，增强花的抗寒能力，减少低温对花器官的伤害。通过对MYBS3基因在不同组织中的表达差异与地被菊耐寒性的关联分析，可以看出MYBS3基因在不同组织中通过调控不同的生理过程，协同作用，共同增强地被菊的耐寒能力。这一研究结果为深入理解地被菊耐寒分子机制提供了重要线索，也为地被菊耐寒品种的选育提供了新的理论依据。四、MYBS3表达功能验证及耐寒调控机制4.1MYBS3基因的功能验证实验设计4.1.1过表达载体构建MYBS3基因过表达载体的构建是深入研究其功能的关键步骤，对揭示地被菊耐寒分子机制具有重要意义。本研究选用植物表达载体pCAMBIA1301，该载体具有广泛的宿主范围，能够在多种植物中高效表达外源基因。其含有潮霉素抗性基因（HygR），可作为筛选标记，便于后续对转化植株进行筛选和鉴定。同时，pCAMBIA1301载体还具备CaMV35S启动子，这是一种组成型启动子，能够驱动外源基因在植物的各个组织和发育阶段持续稳定表达，确保MYBS3基因在转基因地被菊中能够充分发挥作用。根据已克隆的地被菊MYBS3基因序列，运用PrimerPremier5.0软件精心设计引物。在引物两端分别引入EcoRI和XhoI酶切位点，正向引物序列为5'-CGGAATTCATGGCTACAGTACGAC-3'（下划线部分为EcoRI酶切位点），反向引物序列为5'-CCGCTCGAGTCAGTCTTCTGCTCTT-3'（下划线部分为XhoI酶切位点）。引物由生工生物工程（上海）股份有限公司合成。以地被菊‘寒香’叶片的cDNA为模板进行PCR扩增。PCR反应体系为25μL，包含2×TaqMasterMix12.5μL，正向引物（10μM）1μL，反向引物（10μM）1μL，cDNA模板1μL，ddH₂O9.5μL。反应程序设置为：94℃预变性5min；94℃变性30s，58℃退火30s，72℃延伸1min，共进行35个循环；最后72℃终延伸10min。PCR扩增结束后，取5μL扩增产物进行1.0%琼脂糖凝胶电泳检测。在电泳结果中，可观察到在预期大小（约1000bp）处出现清晰明亮的条带，与目的基因MYBS3的大小相符。将PCR产物进行切胶回收，使用AxyPrepDNA凝胶回收试剂盒（AXYGEN公司）按照说明书操作。回收后的DNA片段通过Nanodrop2000分光光度计测定浓度和纯度，结果显示DNA浓度为50ng/μL，OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，表明回收的DNA纯度较高，可用于后续实验。将回收的MYBS3基因片段与经过EcoRI和XhoI双酶切处理的pCAMBIA1301载体进行连接反应。连接反应体系为10μL，包含pCAMBIA1301载体（50ng/μL）1μL，回收的MYBS3基因片段（50ng/μL）4μL，T4DNA连接酶1μL，10×T4DNALigaseBuffer1μL，ddH₂O3μL。16℃连接过夜。连接产物转化大肠杆菌DH5α感受态细胞。将连接产物加入到100μLDH5α感受态细胞中，轻轻混匀，冰浴30min；42℃热激90s，迅速冰浴2min；加入800μL不含抗生素的LB液体培养基，37℃振荡培养1h。将培养后的菌液涂布在含有潮霉素（Hyg）的LB固体培养基平板上，37℃倒置培养16-18h。次日，观察平板上菌落生长情况，挑选3-5个白色菌落接种到含有Hyg的LB液体培养基中，37℃振荡培养过夜。提取质粒进行菌落PCR鉴定和双酶切鉴定。菌落PCR鉴定反应体系和程序与上述PCR扩增相同，取5μL菌落PCR产物进行琼脂糖凝胶电泳检测，在预期大小处出现条带的菌落为阳性克隆。对阳性克隆的质粒进行双酶切鉴定，酶切体系为20μL，包含质粒2μL，10×Buffer2μL，EcoRI1μL，XhoI1μL，ddH₂O14μL。37℃酶切2-3h后，进行1.0%琼脂糖凝胶电泳检测，若出现与目的基因大小相符的条带，则表明过表达载体pCAMBIA1301-MYBS3构建成功。将鉴定正确的克隆送至生工生物工程（上海）股份有限公司进行测序验证。4.1.2基因沉默载体构建基因沉默载体的构建是研究MYBS3基因功能不可或缺的环节，通过抑制MYBS3基因的表达，能够反向验证其在调控地被菊耐寒性中的作用。本研究采用RNA干扰（RNAi）技术构建MYBS3基因沉默载体，该技术利用双链RNA（dsRNA）介导的同源序列特异性基因沉默机制，能够高效、特异地下调目标基因的表达。选用植物RNAi载体pFGC5941，其具有独特的结构特点，含有CaMV35S启动子和ocs终止子，能够驱动干扰片段在植物体内稳定表达。载体中还携带卡那霉素抗性基因（KanR），可用于转化植株的筛选和鉴定。根据MYBS3基因序列，利用在线软件（如dsRNADesigner等）设计干扰片段。干扰片段长度选择为300-500bp，该长度既能保证干扰效果的特异性和高效性，又能避免因片段过长导致的载体构建困难和潜在的非特异性干扰。同时，为了减少脱靶效应，对设计的干扰片段进行BLAST比对分析，确保其与地被菊基因组中其他基因的同源性较低。最终确定的干扰片段序列通过人工合成获得。在干扰片段两端分别添加BamHI和SacI酶切位点，以便后续与载体连接。将合成的干扰片段与经过BamHI和SacI双酶切处理的pFGC5941载体进行连接反应。连接反应体系为10μL，包含pFGC5941载体（50ng/μL）1μL，干扰片段（50ng/μL）4μL，T4DNA连接酶1μL，10×T4DNALigaseBuffer1μL，ddH