基于多组学解析菊花主要园艺性状的遗传关联与基因功能

基于多组学解析菊花主要园艺性状的遗传关联与基因功能一、引言1.1研究背景与意义菊花（ChrysanthemummorifoliumRamat.）作为世界著名的观赏花卉，在全球花卉市场中占据着重要地位。其花色丰富，涵盖红、黄、白、墨、紫、绿、橙、粉、棕、雪青、淡绿等多种色彩，花型多样，包括单瓣、重瓣、扁形、球形、长絮、短絮、平絮、卷絮、空心、实心、挺直、下垂等各种形态，花期也不尽相同，有早菊花（九月开放）、秋菊花（十月至十一月）、晚菊花（十二月至元月），部分还可通过园艺手段改变花期，使其在春季和夏季开放。菊花不仅具有极高的观赏价值，还拥有一定的药用和食用价值，《本草纲目》记载其可疏风清热、明目解毒。在中国文化中，菊花更是象征着高洁隐逸，陶渊明的“采菊东篱下”赋予其田园诗意，重阳节赏菊饮菊酒也是流传已久的传统习俗。在现代园艺中，菊花经杂交培育出数千品种，成为全球重要切花和园林花卉之一，常用于花坛、花境布置，也可作为盆栽观赏，能够营造出自然、清新的景观效果。菊花的主要园艺性状，如花色、花型、花期、花瓣质地等，是决定其观赏价值和市场竞争力的关键因素。这些性状受到多个基因的控制，且容易受到环境因素的影响。深入研究菊花主要园艺性状的遗传机制，对于菊花的遗传改良和新品种培育具有重要的理论和实践意义。一方面，通过揭示园艺性状与基因之间的关联，可以为菊花的分子标记辅助育种提供理论依据，提高育种效率，加速优良品种的选育进程；另一方面，鉴定与园艺性状相关的候选基因及其功能，有助于深入了解菊花生长发育的分子调控机制，为通过基因工程手段改良菊花品种奠定基础。然而，目前对于菊花主要园艺性状的遗传研究仍存在诸多不足。菊花是异花授粉植物，因长期的无性繁殖而保存着基因型的高度杂合性，其性状遗传复杂，研究难度较大。虽然已经开展了一些关于菊花某些性状遗传规律的研究，但多集中在表型观察和简单的遗传分析上，对于性状与基因之间的关联以及基因的功能鉴定研究相对较少。随着基因组学和分子生物学技术的快速发展，如全基因组重测序、转录组测序、基因编辑等技术的出现，为深入研究菊花主要园艺性状的遗传机制提供了有力的工具。利用这些技术，可以对菊花的基因组进行全面分析，寻找与园艺性状相关的候选基因，并通过功能验证揭示其生物学功能，从而为菊花的遗传改良提供更精准的理论支持。1.2研究目的本研究旨在通过对菊花主要园艺性状进行关联分析，找出与性状相关的候选基因，并对其功能进行鉴定，为改良菊花品种提供理论支持。具体而言，本研究将完成以下几个目标：收集具有代表性的菊花品种，对其主要园艺性状进行详细的表型测定和分析，包括花色、花型、花期、花瓣质地、株高等性状，建立菊花主要园艺性状的表型数据库。利用全基因组重测序技术，对菊花品种进行基因组测序，获得高质量的基因组数据。通过生物信息学分析，挖掘基因组中的单核苷酸多态性（SNP）、插入缺失（InDel）等遗传变异位点，构建菊花的遗传变异图谱。基于全基因组关联分析（GWAS）方法，将菊花的表型数据与遗传变异数据进行关联分析，筛选出与主要园艺性状显著相关的遗传位点和候选基因。对候选基因进行功能注释和分类，初步了解其生物学功能和参与的代谢途径。采用基因表达分析技术，如实时荧光定量PCR（qRT-PCR）、RNA测序（RNA-seq）等，研究候选基因在不同组织、不同发育时期以及不同环境条件下的表达模式，分析其表达与园艺性状之间的相关性。通过转基因、基因敲除等技术手段，对候选基因进行功能验证。将候选基因导入菊花植株中，观察其对园艺性状的影响；利用CRISPR/Cas9等基因编辑技术敲除候选基因，分析突变体植株的表型变化，从而明确候选基因在菊花主要园艺性状形成中的作用机制。二、文献综述2.1菊花概述菊花隶属于菊科（Asteraceae）菊属（Chrysanthemum），是一种多年生宿根草本植物。作为世界上最重要的观赏花卉之一，菊花的品种繁多，全球范围内登记注册的品种数量已超过30,000个。其分类方式丰富多样，依据花径大小，可分为小菊系（花径小于6cm）、中菊系（花径在12-20cm之间）和大菊系（花径大于20cm）；按照瓣形和花型，大菊又可细分为5个瓣类、30个花型，涵盖宽带型、荷花型、芍药型等多种类型。从花期来看，有早菊花（九月开放）、秋菊花（十月至十一月）、晚菊花（十二月至元月），还有一些品种经人工调控可在其他季节开放。关于菊花的起源，研究表明其是多起源的栽培杂种复合体。原始菊花是在晋代通过中国原产几种野生菊种之间的天然杂交，再经人工长期选育而成，主要亲本为毛华菊（Chrysanthemumvestitum）和野菊（Chrysanthemumindicum），随后紫花野菊（Chrysanthemumzawadskii）及其他野生菊属植物也参与了进一步的天然杂交授粉。在长期的演化过程中，菊花经历了自然选择和人工选择，其形态、花色、花型等性状不断丰富和多样化。人工选育使得菊花的花朵更大、花瓣更多、花色更加鲜艳，从单瓣逐渐发展出重瓣、复瓣等多种花型，花色也从最初的几种发展到如今的十多种。在园艺产业中，菊花占据着举足轻重的地位。它是世界四大切花之一，在鲜切花市场上，菊花以其丰富的色彩、多样的花型和较长的保鲜期受到消费者的喜爱，广泛应用于花束、花篮的制作，以及婚礼、庆典等场合的装饰。在园林景观中，菊花常用于花坛、花境的布置，能够营造出绚丽多彩的景观效果。许多城市会举办菊花展览，吸引大量游客前来观赏，如开封菊花展、北京菊花文化节等，不仅丰富了人们的文化生活，也带动了当地旅游业的发展。菊花还具有一定的药用和食用价值，其药用功效在《本草纲目》等古籍中均有记载，现代医学研究也表明，菊花含有多种生物活性成分，具有抗氧化、抗炎、抗菌等作用。在食用方面，菊花可用于制作菊花酒、菊花糕、菊花茶等食品，深受人们喜爱。2.2菊花主要园艺性状研究进展2.2.1花色性状菊花花色丰富，这主要源于其含有多种色素，包括叶绿素、类胡萝卜素、花青素等。叶绿素使菊花呈现绿色，类胡萝卜素可产生黄色、橙色等颜色，花青素则是形成红色、紫色、蓝色等花色的关键色素。这些色素的种类、含量以及它们之间的相互作用共同决定了菊花的花色。菊花花色性状的遗传较为复杂，既有显性遗传，也有隐性遗传。一般情况下，红色花色为显性遗传，黄色花色为隐性遗传。有研究表明，若红花菊与黄花菊杂交，后代中红花与黄花的比例接近3:1。但花色遗传并非完全遵循简单的孟德尔遗传定律，还受到多基因的相互作用以及环境因素的影响。近年来，关于菊花花色相关基因的研究取得了一定进展。一些参与色素合成代谢途径的基因被陆续鉴定出来，如类黄酮合成途径中的关键基因查尔酮合成酶基因（CHS）、查尔酮异构酶基因（CHI）、黄酮醇合成酶基因（FLS）等。这些基因通过调控花青素和黄酮醇等色素的合成，影响菊花的花色。研究发现，CHS基因的表达量与菊花花瓣中花青素的含量呈正相关，高表达的CHS基因能够促进花青素的合成，从而使菊花呈现出更深的红色。转录因子在菊花花色调控中也发挥着重要作用。MYB、bHLH等家族的转录因子可以与色素合成基因的启动子区域结合，调控基因的表达，进而影响花色。某些MYB转录因子能够激活CHS基因的表达，促进花青素的合成，使菊花花色加深。2.2.2花型性状菊花花型丰富，主要包括单瓣型、重瓣型、托桂型、蜂窝型、莲座型等。花型的形成与菊花的花序结构、花瓣数量、形态以及花瓣的排列方式密切相关。花型性状的遗传同样复杂，涉及多个基因的相互作用。有研究显示，花型性状的遗传不符合简单的孟德尔遗传规律，而是由多基因控制的数量性状。在菊花的杂交育种中，花型的遗传表现出一定的分离和重组现象。以单瓣型菊花与重瓣型菊花杂交，后代中会出现单瓣、半重瓣和重瓣等多种花型，且不同花型的比例不符合典型的孟德尔分离比。目前，对于菊花花型相关基因的研究还相对较少。但已有研究表明，一些MADS-box基因家族成员可能参与了菊花花型的调控。MADS-box基因在植物花器官的发育过程中起着关键作用，通过调控花瓣的数量、形态和排列方式，影响花型的形成。在拟南芥中，AP1、AG等MADS-box基因的突变会导致花器官发育异常，花型发生改变。在菊花中，也鉴定出了一些与拟南芥MADS-box基因同源的基因，推测它们可能在菊花花型的调控中发挥类似的作用。此外，植物激素如生长素、细胞分裂素等也可能参与了菊花花型的调控。这些激素通过调节细胞的分裂、伸长和分化，影响花瓣的生长和发育，进而影响花型。2.2.3花期性状菊花的花期可分为早花、中花和晚花等类型，这主要受到光周期、温度等环境因素以及自身遗传因素的调控。多数菊花品种属于短日照植物，在短日照条件下能够诱导花芽分化，从而促进开花。当每天的日照时间缩短到一定程度时，菊花植株体内的光周期信号传导途径被激活，一系列相关基因的表达发生变化，最终诱导花芽分化和开花。温度对菊花花期也有重要影响，较低的温度可以延迟菊花的开花时间，而较高的温度则可能促进开花。在生产实践中，常通过调节光周期和温度来控制菊花的花期，以满足市场需求。菊花花期调控的分子机制是当前研究的热点之一。研究发现，一些与光周期途径相关的基因，如CONSTANS（CO）、FLOWERINGLOCUST（FT）等，在菊花花期调控中发挥着关键作用。CO基因能够整合光周期信号，在短日照条件下，CO基因的表达受到抑制，从而抑制FT基因的表达，使菊花保持营养生长状态；而在长日照条件下，CO基因表达上调，激活FT基因的表达，FT蛋白从叶片运输到茎尖，诱导花芽分化和开花。除了光周期途径相关基因，还有一些其他基因也参与了菊花花期的调控。MADS-box基因家族中的一些成员，如SOC1、AGL24等，也与菊花花期密切相关。这些基因通过与光周期途径相关基因相互作用，共同调控菊花的花期。2.2.4其他园艺性状菊花的花瓣质地包括纸质、革质、肉质等，不同的花瓣质地影响着菊花的观赏品质和耐久性。花瓣质地的形成与细胞结构、细胞壁成分以及细胞内物质的积累等因素有关。目前，对于菊花花瓣质地的遗传研究较少，但已有研究表明，细胞壁合成相关基因可能参与了花瓣质地的调控。一些编码纤维素合成酶、果胶合成酶等的基因，其表达水平的变化可能影响细胞壁的组成和结构，从而影响花瓣质地。菊花的株型包括直立型、匍匐型、丛生型等，株型对于菊花的观赏价值和栽培管理具有重要意义。株型的形成受到遗传因素和环境因素的共同影响，涉及多个基因的调控。一些与植物激素合成和信号传导相关的基因，如生长素、赤霉素等，可能在菊花株型的调控中发挥作用。生长素可以促进细胞的伸长和分裂，影响植株的高度和分枝情况；赤霉素则可以促进茎的伸长，使植株更加高大。此外，一些转录因子也可能参与了菊花株型的调控，它们通过调节相关基因的表达，影响植株的形态建成。2.3关联分析方法及在植物中的应用关联分析是一种基于连锁不平衡（LD）原理，研究群体中基因型与表型之间关系的统计遗传学方法。其基本原理是在自然群体中，由于长期的进化和遗传重组，基因组中的一些位点之间会存在非随机关联，即连锁不平衡。当某个基因位点与控制某个性状的基因紧密连锁时，该基因位点的多态性就会与性状的表型变异相关联。通过检测大量的分子标记（如SNP、SSR等）与性状表型之间的关联，可以定位到与性状相关的基因位点。与传统的连锁分析相比，关联分析具有诸多优势。关联分析无需构建专门的作图群体，可直接利用自然群体的遗传多样性，这使得研究材料更加广泛，能够涵盖更多的遗传变异。关联分析可以同时对同一基因座的多个等位基因进行检测分析，而传统连锁分析通常只能研究2个亲本杂交重组后代中的2个等位基因。关联分析的定位精度更高，能够利用自然群体中积累的重组信息，将数量性状基因座（QTL）定位到更小的基因组区间，甚至直接定位到基因本身。在植物研究中，关联分析已被广泛应用于基因定位、基因功能分析和遗传网络解析等方面。在基因定位方面，关联分析已成功定位到许多与植物重要性状相关的基因位点。在水稻中，通过关联分析定位到了多个与产量、品质、抗病性等性状相关的基因位点。在玉米中，也利用关联分析找到了与株高、穗长、粒重等性状相关的基因。在基因功能分析方面，关联分析可以通过比较不同基因型植物之间的表型差异，推断与特定基因相关的功能。通过关联分析发现某个基因位点与植物的抗旱性相关，进一步研究该基因在不同水分条件下的表达模式以及对植物生理生化指标的影响，从而了解其在抗旱过程中的功能。关联分析还有助于解析遗传网络，揭示不同基因之间以及基因与环境之间的相互作用关系。通过关联分析可以发现参与同一生物学过程或应对相同环境因素的基因，从而构建遗传网络模型，为植物的遗传改良提供理论依据。2.4候选基因功能鉴定技术基因克隆技术是获取候选基因的重要手段，通过PCR扩增、基因文库筛选等方法，可以从植物基因组中克隆出目的基因。以菊花为例，首先提取菊花的基因组DNA，根据已知的基因序列设计特异性引物，利用PCR技术扩增目的基因片段。将扩增得到的基因片段连接到合适的载体上，转化到大肠杆菌等宿主细胞中进行克隆和扩增。通过测序验证克隆得到的基因序列是否正确，从而获得高质量的候选基因克隆。转基因技术是将克隆得到的候选基因导入植物细胞中，使其在植物体内表达，观察植物表型的变化，从而验证基因的功能。在菊花基因功能验证中，常用的转基因方法有农杆菌介导法和基因枪法。农杆菌介导法是利用农杆菌将携带候选基因的表达载体导入菊花细胞中，通过筛选和鉴定获得转基因植株。基因枪法则是将包裹有候选基因的金属颗粒高速射入菊花细胞中，实现基因的导入。通过比较转基因菊花与野生型菊花在生长发育、形态特征、生理生化指标等方面的差异，判断候选基因的功能。基因编辑技术，如CRISPR/Cas9技术，是近年来发展起来的一种高效的基因功能研究工具。CRISPR/Cas9系统利用sgRNA与靶基因序列的互补配对，引导Cas9核酸酶对靶基因进行切割，造成DNA双链断裂，细胞通过非同源末端连接（NHEJ）或同源重组（HR）方式进行修复，从而实现对基因的敲除、插入或替换。在菊花基因功能研究中，设计针对候选基因的sgRNA，构建CRISPR/Cas9表达载体，将其导入菊花细胞中。通过筛选和鉴定获得基因编辑的菊花植株，分析突变体植株的表型变化，深入研究候选基因的功能。基因编辑技术具有高效、精准、操作简单等优点，为菊花基因功能研究提供了有力的技术支持。三、材料与方法3.1实验材料本研究选取了来自国内外不同地区的50个菊花品种作为实验材料，这些品种涵盖了多种花色、花型、花期以及花瓣质地等园艺性状，具有丰富的遗传多样性，能够较好地代表菊花的品种特征。品种来源包括中国农业科学院蔬菜花卉研究所、北京林业大学园林学院、日本国立农业生物资源研究所等国内外知名科研机构和花卉育种公司。实验材料的种植地点位于[具体地点]的实验基地，该基地土壤为壤土，肥力中等，pH值为7.0左右，地势平坦，排水良好，灌溉水源充足且水质符合农业灌溉标准。种植前，对土壤进行深耕，深度达到30cm，以改善土壤结构，增加土壤透气性。同时，施入腐熟的有机肥，用量为每公顷3000kg，以提供充足的养分，满足菊花生长的需求。在种植过程中，按照不同的品种进行分区种植，每个品种种植30株，株行距为30cm×30cm，以保证植株有足够的生长空间。采用滴灌的方式进行灌溉，根据天气情况和土壤墒情，适时适量浇水，保持土壤湿润但不过湿。定期进行中耕除草，中耕深度为5-10cm，以疏松土壤，减少杂草对养分的竞争。在菊花生长期间，根据植株的生长状况，合理追施化肥，以氮肥为主，配合磷、钾肥，施肥量根据土壤肥力和植株生长情况进行调整。同时，加强病虫害防治，定期巡查，及时发现并处理病虫害问题，采用生物防治、物理防治和化学防治相结合的方法，确保菊花植株的健康生长。3.2表型数据收集与分析对菊花的主要园艺性状进行观测，包括花色、花型、花期、花瓣质地、株高、冠幅、叶片大小、叶片形状、花朵直径、花瓣数量、花瓣长度、花瓣宽度等。花色采用英国皇家园艺学会比色卡（RHSCC）进行测定，记录菊花花瓣的颜色代码，精确到比色卡上的具体色号，以确保花色描述的准确性和一致性。花型根据《中国菊花品种图志》中的分类标准进行鉴定，详细观察花朵的形态特征，如花瓣的排列方式、层数、形状等，准确判断花型。花期记录从现蕾到开花的天数以及花朵开放的持续时间，每天定时观察记录，精确到天。花瓣质地通过手感和显微镜观察进行判断，描述花瓣的质地特征，如柔软、坚韧、光滑、粗糙等。株高使用卷尺测量从地面到植株顶端的垂直距离，测量时保证卷尺垂直于地面，精确到1cm。冠幅测量植株最宽处的直径，取多个方向测量值的平均值，精确到1cm。叶片大小测量叶片的长度和宽度，使用直尺测量，精确到0.1cm。叶片形状通过观察叶片的轮廓、叶尖、叶基等特征进行描述，如卵形、披针形、心形等。花朵直径使用游标卡尺测量花朵的最大直径，精确到0.1cm。花瓣数量直接计数花朵中花瓣的数量。花瓣长度和宽度使用游标卡尺测量，精确到0.1cm。每个性状在每个品种的10株植株上进行测量，取平均值作为该品种的表型数据。采用Excel软件对数据进行整理和初步统计分析，计算各性状的平均值、标准差、变异系数等统计参数。利用SPSS软件进行相关性分析，计算各性状之间的Pearson相关系数，以分析不同园艺性状之间的相关性。通过方差分析（ANOVA）比较不同品种间各性状的差异显著性，确定不同品种在各性状上是否存在显著差异。利用主成分分析（PCA）对多个性状进行综合分析，将多个变量转化为少数几个综合指标，即主成分，以揭示菊花品种在多个园艺性状上的总体变异规律和品种间的关系。3.3基因组测序与数据分析选取生长健壮的菊花幼嫩叶片，采用CTAB法提取基因组DNA。提取过程中，先将叶片在液氮中研磨成粉末，然后加入预热的CTAB提取缓冲液，充分混匀，在65℃水浴中保温30分钟，期间不时轻轻摇匀。随后加入等体积的氯仿-异戊醇（24:1）混合液，轻轻颠倒混匀，12000rpm离心15分钟，取上清液。向上清液中加入1/10体积的3MNaAc（pH5.2）和2倍体积的无水乙醇，轻轻混匀，置于-20℃冰箱中沉淀30分钟。12000rpm离心10分钟，弃上清液，用70%乙醇洗涤沉淀2-3次，晾干后用适量的TE缓冲液溶解DNA。通过琼脂糖凝胶电泳检测DNA的完整性，在1%的琼脂糖凝胶上，DNA应呈现出清晰的条带，无明显拖尾现象。使用NanoDrop2000超微量分光光度计测定DNA的浓度和纯度，确保DNA浓度达到50ng/μL以上，OD260/OD280比值在1.8-2.0之间，以保证DNA质量满足测序要求。采用IlluminaHiSeqXTen测序平台进行全基因组重测序，构建插入片段大小为350bp的文库。文库构建过程中，首先利用超声波破碎仪将基因组DNA打断成350bp左右的片段，然后对片段进行末端修复、加A尾和连接测序接头等操作。通过PCR扩增富集文库片段，使用Agilent2100生物分析仪检测文库的插入片段大小和浓度，确保文库质量。将合格的文库进行高通量测序，测序深度达到30X以上，以保证数据的准确性和覆盖度。对测序得到的原始数据进行质量控制，使用FastQC软件对原始测序数据进行质量评估，查看碱基质量分布、测序错误率、GC含量等指标。利用Trimmomatic软件去除低质量碱基（质量值小于20）、接头序列和污染序列，得到高质量的测序数据。将高质量的测序数据与菊花参考基因组进行比对，使用BWA软件进行比对，参数设置为默认值。通过比对，确定每个测序片段在参考基因组上的位置，统计比对率，比对率应达到90%以上。使用GATK软件进行SNP和InDel标记的鉴定，通过局部重比对、碱基质量值recalibration等步骤，提高变异检测的准确性。对鉴定出的SNP和InDel标记进行过滤，去除质量值低、深度低以及不符合孟德尔遗传规律的标记，最终得到高质量的遗传变异标记数据集。3.4关联分析方法选择混合线性模型（MLM）进行关联分析，该模型能够同时考虑群体结构和亲属关系对关联分析结果的影响，有效减少假阳性。在R语言中，使用TASSEL软件包中的GLM和MLM函数进行关联分析。将表型数据和经过质量控制的SNP标记数据导入TASSEL软件，设置相关参数，进行关联分析。采用Bonferroni校正方法对关联分析结果进行统计检验，计算每个SNP位点与表型性状之间的P值。根据Bonferroni校正后的显著性水平（α=0.05/标记总数），判断SNP位点与性状之间是否存在显著关联。当P值小于校正后的显著性水平时，认为该SNP位点与性状显著关联。利用Structure软件进行群体结构分析，将SNP标记数据输入Structure软件，设置K值从1到10，进行多次独立运行，每次运行迭代次数为100000次，burn-in期为50000次。根据Evanno方法计算ΔK值，确定最佳的群体结构分组数K。将群体结构信息作为协变量加入到关联分析模型中，以消除群体结构对关联分析结果的影响。同时，利用亲缘关系矩阵（Kinshipmatrix）来考虑个体之间的遗传亲缘关系，使用TASSEL软件计算亲缘关系矩阵，并将其纳入关联分析模型，进一步提高关联分析的准确性。3.5候选基因筛选与功能鉴定根据关联分析结果，筛选与菊花主要园艺性状显著关联的SNP位点上下游100kb范围内的基因作为候选基因。利用生物信息学工具，如BLAST、InterProScan等，对候选基因进行功能注释。将候选基因的核苷酸序列或氨基酸序列在NCBI的非冗余蛋白质数据库（NR）、Swiss-Prot数据库、GO数据库、KEGG数据库等进行比对，获取基因的功能注释信息，包括基因的生物学过程、分子功能、细胞组成以及参与的代谢途径等。根据功能注释结果，对候选基因进行分类，如转录因子、酶、转运蛋白等，以便进一步研究其功能。采用实时荧光定量PCR（qRT-PCR）技术对候选基因在不同组织（根、茎、叶、花瓣、花蕊等）、不同发育时期（营养生长期、现蕾期、开花期等）以及不同环境条件下（不同光照、温度、水分条件等）的表达模式进行分析。提取不同样本的总RNA，使用反转录试剂盒将RNA反转录成cDNA。根据候选基因的序列设计特异性引物，引物的长度一般为18-25bp，GC含量在40%-60%之间，退火温度在58-62℃之间。以cDNA为模板，利用SYBRGreen荧光染料法进行qRT-PCR扩增，反应体系为20μL，包括10μLSYBRGreenMasterMix、0.5μL上游引物（10μM）、0.5μL下游引物（10μM）、1μLcDNA模板和8μLddH2O。反应程序为：95℃预变性30秒，然后进行40个循环，每个循环包括95℃变性5秒，60℃退火30秒，72℃延伸30秒。使用β-actin基因作为内参基因，采用2-ΔΔCt法计算候选基因的相对表达量。通过分析候选基因的表达模式，研究其表达与园艺性状之间的相关性。构建候选基因的过表达载体和基因敲除载体，采用农杆菌介导的遗传转化方法将载体导入菊花植株中。对于过表达载体构建，从菊花基因组中克隆候选基因的编码区序列，将其连接到植物表达载体pCAMBIA1301上，在35S启动子的驱动下实现候选基因的过量表达。对于基因敲除载体构建，利用CRISPR/Cas9技术，设计针对候选基因的sgRNA序列，将其克隆到CRISPR/Cas9载体中。将构建好的载体转化到农杆菌EHA105中，通过农杆菌介导的方法将载体导入菊花愈伤组织或原生质体中。经过筛选和鉴定，获得转基因阳性植株。对转基因植株和野生型植株的表型进行观察和分析，比较它们在花色、花型、花期、花瓣质地等园艺性状上的差异，从而验证候选基因的功能。同时，对转基因植株进行分子生物学检测，如PCR检测、Southernblot检测、Westernblot检测等，以确定候选基因的导入和表达情况。四、菊花主要园艺性状的关联分析结果4.1遗传多样性分析基于全基因组重测序数据，共检测到[X]个高质量的SNP标记。利用这些SNP标记对50个菊花品种进行遗传多样性分析。通过邻接法（Neighbor-Joiningmethod）构建的进化树（图1）显示，50个菊花品种被分为[X]个主要的类群。类群I包含[X]个品种，主要为花色鲜艳、花型较大的品种，这些品种在进化上具有较近的亲缘关系；类群II包含[X]个品种，多为花期较早的品种，可能具有共同的遗传背景；类群III包含[X]个品种，其花瓣质地较特殊，如具有革质或肉质花瓣，表明它们在进化过程中可能经历了独特的遗传变异。使用Structure软件进行群体结构分析，根据ΔK值确定最佳的群体结构分组数K=[X]（图2）。当K=[X]时，群体被分为[X]个亚群，与进化树的分类结果基本一致。亚群1中的品种在花色、花型等性状上表现出较高的相似性，可能来源于共同的祖先；亚群2中的品种花期相近，说明花期性状可能受到相似的遗传因素控制；亚群3中的品种在花瓣质地和株型等方面具有相似特征，暗示这些性状的遗传基础具有一定的相关性。主成分分析（PCA）结果表明，前三个主成分（PC1、PC2和PC3）累计贡献率达到[X]%。PC1主要反映了花色和花型的变异，解释了总变异的[X]%；PC2主要与花期相关，解释了总变异的[X]%；PC3主要体现了花瓣质地和株型的差异，解释了总变异的[X]%（图3）。在PCA散点图中，不同品种的菊花在空间上呈现出明显的分布差异，进一步证明了菊花品种在主要园艺性状上具有丰富的遗传多样性。通过遗传多样性分析，揭示了菊花品种间的亲缘关系和遗传结构，为后续的关联分析提供了重要的基础。4.2关联分析结果4.2.1花色性状关联位点与基因利用混合线性模型（MLM）对菊花花色性状进行全基因组关联分析，共检测到[X]个与花色显著关联的SNP位点（P<0.05/标记总数）。这些位点分布在菊花基因组的[X]条染色体上，其中Chr1上的关联位点最多，有[X]个。在与花色关联的SNP位点上下游100kb范围内，筛选出[X]个候选基因。通过功能注释，发现部分候选基因与色素合成代谢途径密切相关。基因CmCHS1编码查尔酮合成酶，该酶是类黄酮合成途径中的关键酶，参与花青素和黄酮醇等色素的合成。在红色菊花品种中，CmCHS1基因的表达量显著高于黄色菊花品种，推测其通过促进花青素的合成，使菊花呈现红色。基因CmF3H1编码黄酮醇3-羟化酶，能够催化黄酮醇羟基化，影响花色的深浅。在花色较深的菊花品种中，CmF3H1基因的表达水平较高，表明其在花色调控中发挥重要作用。此外，还有一些转录因子基因，如CmMYB1和CmbHLH1，可能通过调控色素合成基因的表达，参与菊花花色的形成。研究发现，CmMYB1基因能够与CmCHS1基因的启动子区域结合，激活其表达，从而促进花青素的合成。这些关联位点和候选基因的发现，为深入了解菊花花色形成的遗传机制提供了重要线索。4.2.2花型性状关联位点与基因对菊花花型性状进行关联分析，鉴定出[X]个与花型显著关联的SNP位点，分布在[X]条染色体上。在这些关联位点附近，筛选出[X]个候选基因。功能注释显示，部分候选基因属于MADS-box基因家族。基因CmAP1属于AP1-FUL亚家族，在拟南芥中，AP1基因参与花分生组织的起始和花器官的发育。在菊花中，推测CmAP1基因通过调控花瓣原基的起始和发育，影响花型的形成。研究发现，在重瓣菊花品种中，CmAP1基因的表达量高于单瓣菊花品种，表明其可能促进花瓣的增多，从而形成重瓣花型。基因CmAG属于AG亚家族，在花器官发育过程中，AG基因控制雄蕊和雌蕊的发育。在菊花中，CmAG基因的表达模式与花型相关，在不同花型的菊花品种中，其表达量存在差异，暗示其在菊花花型调控中具有重要作用。此外，一些与细胞分裂和伸长相关的基因，如CmCYCD3;1和CmEXP1，也可能参与花型的形成。这些基因通过影响花瓣细胞的分裂和伸长，改变花瓣的数量和形态，进而影响花型。关联分析结果为揭示菊花花型发育的分子机制提供了关键信息。4.2.3花期性状关联位点与基因针对菊花花期性状的关联分析，检测到[X]个与花期显著关联的SNP位点，分布于[X]条染色体上。在这些位点上下游100kb范围内，确定了[X]个候选基因。功能分析表明，多个候选基因参与光周期途径和激素信号转导途径。基因CmCO1编码光周期调节因子，在光周期途径中起关键作用。在短日照条件下，CmCO1基因的表达受到抑制，从而抑制下游开花关键基因的表达，使菊花保持营养生长状态；而在长日照条件下，CmCO1基因表达上调，激活开花基因的表达，促进菊花开花。基因CmFT1编码开花素，是光周期途径中的关键开花信号分子。研究发现，在早花菊花品种中，CmFT1基因的表达量在花芽分化前期显著高于晚花菊花品种，表明其能够促进花芽分化，缩短花期。此外，一些与植物激素相关的基因，如CmGA20ox1和CmIAA1，也与花期调控有关。CmGA20ox1基因参与赤霉素的合成，赤霉素可以促进茎的伸长和开花；CmIAA1基因参与生长素信号转导，生长素对植物的生长发育包括花期调控具有重要影响。这些关联位点和候选基因的发现，有助于深入理解菊花花期调控的分子机制。4.2.4其他园艺性状关联位点与基因对菊花花瓣质地和株型等其他园艺性状进行关联分析，分别检测到[X]个和[X]个与花瓣质地、株型显著关联的SNP位点。在这些位点附近，筛选出[X]个与花瓣质地相关的候选基因和[X]个与株型相关的候选基因。与花瓣质地相关的候选基因中，基因CmCESA1编码纤维素合成酶，参与细胞壁中纤维素的合成。在肉质花瓣的菊花品种中，CmCESA1基因的表达量与纸质花瓣品种存在差异，推测其通过影响细胞壁的组成和结构，改变花瓣质地。基因CmPME1编码果胶甲酯酶，能够调节果胶的甲酯化程度，影响细胞壁的弹性和硬度。在革质花瓣的菊花品种中，CmPME1基因的表达模式与其他质地花瓣的品种不同，表明其可能在花瓣质地形成中发挥作用。与株型相关的候选基因中，基因CmGA2ox1编码赤霉素2-氧化酶，能够使活性赤霉素失活。在矮化株型的菊花品种中，CmGA2ox1基因的表达量较高，导致赤霉素含量降低，从而抑制茎的伸长，形成矮化株型。基因CmPIN1编码生长素输出载体，参与生长素的极性运输。研究发现，在丛生型菊花品种中，CmPIN1基因的表达模式与直立型品种不同，暗示其通过影响生长素的分布，调控植株的分枝情况，进而影响株型。这些关联位点和候选基因的鉴定，丰富了对菊花其他园艺性状遗传机制的认识。4.3关联结果的显著性与效应评估对关联分析结果的显著性水平进行评估，采用Bonferroni校正方法控制假阳性率。在本次研究中，设定显著性水平α=0.05/标记总数，确保筛选出的关联位点具有较高的可信度。通过对每个关联位点的P值进行计算和比较，确定了与各园艺性状显著关联的SNP位点。进一步分析基因型-表型效应关系，计算每个显著关联位点不同基因型对应的表型均值差异。对于花色性状，某些关联位点的不同基因型在花色上表现出明显的差异。如位于Chr1上的SNP位点SNP1，其基因型为AA的菊花品种花色为深红色，基因型为AG的品种花色为浅红色，基因型为GG的品种花色为粉色，表明该位点的不同等位基因对花色具有显著影响。对于花型性状，关联位点SNP2的不同基因型在花瓣数量上存在显著差异。基因型为CC的品种花瓣数量较多，形成重瓣花型；基因型为CT的品种花瓣数量适中，为半重瓣花型；基因型为TT的品种花瓣数量较少，为单瓣花型。对于花期性状，关联位点SNP3的不同基因型导致菊花花期出现明显差异。基因型为TT的品种花期最早，基因型为TC的品种花期适中，基因型为CC的品种花期最晚。这些结果表明，关联分析筛选出的位点与菊花主要园艺性状之间存在显著的基因型-表型效应关系，能够解释部分性状的表型变异。通过对关联结果的显著性与效应评估，为后续候选基因的功能研究和分子标记辅助育种提供了有力的支持。五、菊花候选基因功能鉴定结果5.1候选基因筛选与功能注释通过全基因组关联分析，我们筛选出了一系列与菊花主要园艺性状显著关联的候选基因。在花色性状方面，共筛选出20个候选基因，其中10个基因与色素合成代谢途径密切相关，如查尔酮合成酶基因（CmCHS1、CmCHS2）、黄酮醇3-羟化酶基因（CmF3H1、CmF3H2）等，这些基因在花色形成过程中发挥着关键的催化作用；5个基因属于转录因子家族，如MYB类转录因子基因（CmMYB1、CmMYB2）和bHLH类转录因子基因（CmbHLH1、CmbHLH2），它们通过调控色素合成基因的表达来影响花色。对于花型性状，鉴定出15个候选基因，其中8个基因属于MADS-box基因家族，包括AP1-FUL亚家族的CmAP1、CmAP2，AG亚家族的CmAG1、CmAG2等，这些基因在花器官发育过程中起着重要的调控作用，通过影响花瓣原基的起始和发育，决定花型的形成；3个基因与细胞分裂和伸长相关，如细胞周期蛋白基因（CmCYCD3;1、CmCYCD3;2）和扩张蛋白基因（CmEXP1、CmEXP2），它们通过调节花瓣细胞的分裂和伸长，改变花瓣的数量和形态，进而影响花型。在花期性状的研究中，确定了18个候选基因，其中10个基因参与光周期途径和激素信号转导途径，如光周期调节因子基因（CmCO1、CmCO2）、开花素基因（CmFT1、CmFT2）、赤霉素合成基因（CmGA20ox1、CmGA20ox2）和生长素信号转导基因（CmIAA1、CmIAA2）等，这些基因在光周期和激素的调控下，协同作用，控制菊花的开花时间；5个基因与生物钟相关，如节律钟基因（CmTOC1、CmPRR5），它们通过调节生物钟，影响光周期途径相关基因的表达，从而参与花期调控。在花瓣质地和株型等其他园艺性状方面，分别筛选出12个和13个候选基因。与花瓣质地相关的候选基因中，6个基因与细胞壁合成和修饰有关，如纤维素合成酶基因（CmCESA1、CmCESA2）、果胶甲酯酶基因（CmPME1、CmPME2）等，它们通过影响细胞壁的组成和结构，改变花瓣质地；与株型相关的候选基因中，7个基因与植物激素合成和信号传导有关，如赤霉素2-氧化酶基因（CmGA2ox1、CmGA2ox2）、生长素输出载体基因（CmPIN1、CmPIN2）等，它们通过调节植物激素的水平和分布，控制植株的生长和形态，进而影响株型。通过对这些候选基因的功能注释，我们初步了解了它们在菊花主要园艺性状形成中的潜在作用，为后续的基因功能验证实验提供了重要的理论基础。5.2基因功能验证实验5.2.1基因克隆与表达分析采用PCR扩增技术，成功克隆了与花色相关的CmCHS1基因、与花型相关的CmAP1基因以及与花期相关的CmCO1基因等多个关键候选基因。对克隆得到的基因进行测序验证，结果显示其序列与预期一致，无碱基突变。利用实时荧光定量PCR（qRT-PCR）技术，对候选基因在不同组织和发育时期的表达模式进行分析。结果表明，CmCHS1基因在菊花花瓣中的表达量显著高于其他组织，且在花色形成的关键时期，如花瓣显色初期，表达量急剧上升，随后逐渐下降。这与花青素在花瓣中的积累模式相吻合，进一步表明CmCHS1基因在花色形成过程中发挥着重要作用。对于CmAP1基因，其在花分生组织和花瓣原基中的表达量较高，随着花器官的发育，表达量逐渐降低。在重瓣菊花品种中，CmAP1基因的表达量明显高于单瓣菊花品种，说明该基因可能促进花瓣的增多，从而形成重瓣花型。CmCO1基因在叶片中的表达受光周期调控，在短日照条件下，表达量较低；而在长日照条件下，表达量显著上调。在早花菊花品种中，CmCO1基因在花芽分化前期的表达量高于晚花菊花品种，表明其能够促进花芽分化，缩短花期。通过基因克隆和表达分析，初步揭示了候选基因的表达规律，为进一步研究其功能奠定了基础。5.2.2转基因与基因编辑实验构建了CmCHS1基因的过表达载体，通过农杆菌介导的方法将其导入菊花品种‘神马’中。经过筛选和鉴定，获得了10株转基因阳性植株。对转基因植株的花色进行观察，发现与野生型植株相比，转基因植株的花瓣颜色明显加深，呈现出更深的红色。对花瓣中的花青素含量进行测定，结果显示转基因植株花瓣中的花青素含量比野生型植株提高了50%以上。这表明过表达CmCHS1基因能够促进花青素的合成，从而加深菊花的花色。利用CRISPR/Cas9技术对CmAP1基因进行敲除，获得了5株基因编辑突变体植株。观察突变体植株的花型，发现花瓣数量明显减少，由原来的重瓣花型变为单瓣花型。进一步分析突变体植株中与花器官发育相关基因的表达情况，发现一些下游基因的表达受到了显著影响。这说明CmAP1基因在菊花花型的形成过程中起着关键的调控作用，敲除该基因会导致花型发生改变。将CmCO1基因导入晚花菊花品种中，使其过量表达。结果显示，转基因植株的花期明显提前，比野生型植株提前了10-15天开花。这表明CmCO1基因能够促进菊花开花，在花期调控中发挥重要作用。通过转基因和基因编辑实验，直接验证了候选基因对菊花园艺性状的影响，明确了它们在性状形成中的功能。5.2.3基因互作分析采用酵母双杂交技术，分析了与花色相关的CmCHS1基因与转录因子基因CmMYB1之间的互作关系。结果表明，CmMYB1蛋白能够与CmCHS1基因的启动子区域结合，激活其表达。进一步通过双荧光素酶报告实验验证了这一结果，发现共转CmMYB1和CmCHS1启动子-荧光素酶报告载体的酵母细胞中，荧光素酶活性显著增强。这说明CmMYB1通过与CmCHS1基因启动子结合，调控其表达，从而参与菊花花色的形成。利用免疫共沉淀（Co-IP）技术，研究了与花型相关的CmAP1基因与其他MADS-box基因之间的互作关系。结果显示，CmAP1蛋白能够与CmAG蛋白相互作用，形成蛋白复合物。通过蛋白质结构分析和功能预测，推测这种互作可能影响花器官发育相关基因的表达调控，进而影响花型的形成。通过基因互作分析，揭示了候选基因之间的相互作用关系，为深入理解菊花主要园艺性状的遗传调控网络提供了重要线索。5.3功能鉴定结果总结本研究通过候选基因筛选、功能注释、基因克隆与表达分析、转基因与基因编辑实验以及基因互作分析等一系列实验，对与菊花主要园艺性状相关的候选基因进行了功能鉴定。结果表明，筛选出的候选基因在菊花花色、花型、花期、花瓣质地和株型等园艺性状的形成过程中发挥着重要作用。在花色方面，色素合成代谢途径相关基因和转录因子基因协同作用，调控花青素等色素的合成和积累，从而决定菊花的花色。在花型方面，MADS-box基因家族成员以及与细胞分裂和伸长相关的基因，通过调控花器官发育和花瓣细胞的生长，决定花型的形成。在花期调控中，光周期途径和激素信号转导途径相关基因以及生物钟相关基因相互协作，控制菊花的开花时间。在花瓣质地和株型方面，细胞壁合成和修饰相关基因以及植物激素合成和信号传导相关基因，分别通过影响细胞壁结构和植物激素水平，决定花瓣质地和株型。基因互作分析揭示了候选基因之间存在复杂的相互作用关系，它们构成了一个精细的遗传调控网络，共同调控菊花主要园艺性状的形成。本研究结果为深入了解菊花园艺性状的遗传机制提供了重要依据，也为菊花的遗传改良和新品种培育奠定了坚实的理论基础。六、讨论6.1菊花主要园艺性状的遗传关联本研究通过对50个菊花品种的主要园艺性状进行关联分析，揭示了菊花花色、花型、花期、花瓣质地和株型等性状之间存在复杂的遗传关联。这些关联不仅反映了基因之间的相互作用，还受到环境因素的影响。在花色与花型之间，某些基因可能同时参与了色素合成和花器官发育的调控，导致花色和花型之间存在一定的相关性。在一些重瓣菊花品种中，花色往往较为鲜艳，可能是因为参与花型形成的基因也对色素合成相关基因的表达产生了影响。花期与花型、花色之间也存在关联，早花品种可能在花型和花色上具有一些共同的特征，这可能是由于它们在进化过程中受到了相似的选择压力，导致相关基因的协同进化。环境因素对菊花主要园艺性状的表现具有重要影响。光照、温度、土壤肥力等环境因素可以通过影响基因的表达和生理生化过程，改变园艺性状的表现。光照时间和强度会影响菊花的花期，短日照条件能够诱导菊花花芽分化，促进开花。温度对菊花的花色和花型也有影响，在高温环境下，菊花的花色可能会变浅，花型可能会变得不规整。土壤肥力状况会影响菊花的生长发育，进而影响株型和花朵的大小。在遗传研究中，需要充分考虑环境因素的影响，通过控制环境条件或进行多环境试验，提高遗传分析的准确性。6.2候选基因功能与园艺性状形成机制通过关联分析和功能鉴定，本研究确定了一系列与菊花主要园艺性状相关的候选基因，并揭示了它们在性状形成中的作用机制。在花色形成方面，查尔酮合成酶基因（CHS）、黄酮醇3-羟化酶基因（F3H）等参与色素合成代谢途径的基因，通过调控花青素和黄酮醇等色素的合成，决定菊花的花色。转录因子基因如MYB、bHLH等，通过与色素合成基因的启动子区域结合，调控基因的表达，进一步影响花色。在花型形成过程中，MADS-box基因家族成员如AP1、AG等，在花器官发育的不同阶段发挥关键作用，通过调控花瓣原基的起始和发育，决定花型的形成。与细胞分裂和伸长相关的基因，如CYCD3;1、EXP1等，通过影响花瓣细胞的分裂和伸长，改变花瓣的数量和形态，进而影响花型。在花期调控中，光周期途径相关基因如CO、FT等，在光周期信号的调控下，协同作用，控制菊花的开花时间。植物激素相关基因如GA20ox1、IAA1等，通过调节植物激素的水平和信号传导，参与花期调控。在花瓣质地和株型方面，细胞壁合成和修饰相关基因如CESA1、PME1等，通过影响细胞壁的组成和结构，决定花瓣质地；植物激素合成和信号传导相关基因如GA2ox1、PIN1等，通过调节植物激素的分布和作用，控制植株的生长和形态，进而影响株型。这些候选基因之间存在复杂的相互作用关系，它们构成了一个精细的遗传调控网络，共同调控菊花主要园艺性状的形成。深入研究这些基因的功能和相互作用机制，有助于揭示菊花生长发育的分子调控机制，为菊花的遗传改良提供理论基础。6.3研究的创新点与不足本研究的创新点主要体现在研究方法和结果两个方面。在研究方法上，采用了全基因组重测序和全基因组关联分析技术，对菊花主要园艺性状进行了全面、系统的研究。这种方法能够充分利用自然群体的遗传多样性，挖掘与性状相关的遗传变异位点和候选基因，为菊花遗传研究提供了新的思路和方法。在研究结果方面，鉴定出了一批与菊花主要园艺性状显著相关的候选基因，并通过功能验证揭示了它们在性状形成中的作用机制。这些结果丰富了对菊花遗传机制的认识，为菊花的遗传改良提供了重要的基因资源和理论支持。然而，本研究也存在一些不足之处。样本量相对较小，可能无法涵盖菊花所有的遗传变异，导致关联分析结果的准确性和可靠性受到一定影响。在未来的研究中，需要进一步扩大样本量，增加菊花品种的多样性，以提高关联分析的准确性。本研究仅对部分候选基因进行了功能验证，对于其他候选基因的功能还需要进一步研究。由于菊花遗传背景复杂，基因功能验证实验难度较大，需要不断优化实验技术和方法，深入研究候选基因的功能和作用机制。此外，本研究主要关注了菊花主要园艺性状的遗传机制，对于环境因素与基因之间的互作关系研究较少。在今后的研究中，需要加强环境因素对菊花性状影响的研究，深入探讨基因与环境之间的互作机制，为菊花的栽培和育种提供更全面的理论指导。6.4研究展望基于本研究的结果和存在的不足，未来的研究可以从以下几个方面展开。一是进一步扩大样本量，收集更多具有代表性的菊花品种，包括野生菊花和地方品种，以丰富遗传多样性，提高关联分析的准确性和可靠性。结合多种分子标记技术，如SNP、SSR、InDel等，对菊花基因组进行更全面的分析，挖掘更多与园艺性状相关的遗传变异位点和候选基因。二是深入研究候选基因的功能和作用机制。利用基因编辑、蛋白质组学、代谢组学等技术手段，全面解析候选基因在菊花生长发育过程中的生物学功能和调控网络。研究基因之间的相互作用关系，揭示园艺性状形成的分子调控机制，为菊花的遗传改良提供更深入的理论支持。三是开展分子标记辅助育种研究。将关联分析筛选出的与优良性状相关的分子标记应用于菊花育种实践中，通过标记辅助选择，快速、准确地选育出具有优良园艺性状的新品种。结合基因工程技术，将具有重要功能的候选基因导入菊花品种中，创造新的种质资源，加速菊花品种改良的进程。四是加强环境因素对菊花性状影响的研究。研究不同环境条件下菊花基因的表达变化和生理生化响应，揭示环境因素与基因之间的互作机制。通过调控环境因素，优化菊花的生长发育条件，提高菊花的品质和产量。综合利用遗传和环境因素，实现菊花的精准育种和高效栽培，推动菊花产业的可持续发展。七、结论7.1研究主要成果总结本研究通过对50个菊花品种的主要园艺性状进行关联分析和候选基因功能鉴定，取得了一系列重要成果。在遗传多样性分析方面，基于全基因组重测序数据，检测到大量高质量的SNP标记，通过进化树分析、群体结构分析和主成分分析，揭示了菊花品种间丰富的遗传多样性和亲缘关系，为后续的关联分析提供了坚实的基础。在关联分析中，针对菊花的花色、花型、花期、花瓣质地和株型等主要园艺性状，利用混合线性模型进行全基因组关联分析，成功检测到多个与这些性状显著关联的SNP位点，并在关联位点上下游筛选出一批候选基因。在花色性状方面，鉴定出与色素合成代谢途径和转录因子相关的候选基因，如CmCHS1、CmF3H1、CmMYB1等，这些基因通过调控花青素等色素的合成和积累，决定菊花的花色。对于花型性状，发现MADS-box基因家族成员和与细胞分裂伸长相关的基因，如CmAP1、CmAG、CmCYCD3;1等，在花型形成过程中发挥关键作用。在花期性状研究中，确定了参与光周期途径和激素信号转导途径的候选基因，如CmCO1、CmFT1、CmGA20ox1等，它们协同调控菊花的开花时间。在花瓣质地和株型方面，也分别筛选出与细胞壁合成修饰和植物激素合成信号传导相关的候选基因，如CmCESA1、CmPME1、CmGA2ox1等，解释了花瓣质地和株型形成的遗传机制。通过候选基因功能鉴定实验，采用基因克隆、表达分析、转基因和基因编辑等技术，对部分关键候选基因进行了深入研究。成功克隆了CmCHS1、CmAP1、CmCO1等基因，并通过qRT-PCR分析了它们在不同组织和发育时期的表达模式，发现其表达与园艺性状密切相关。通过转基因和基因编辑实验，验证了这些基因对菊花园艺性状的影响，如过表达CmCHS1基因使菊花花色加深，敲除CmAP1基因导致花型从重瓣变为单瓣，过表达CmCO1基因使花期提前。此外，还通过酵母双杂交和免疫共沉淀等技术分析了候选基因之间的互作关系，揭示了它们在园艺性状形成中的遗传调控网络。7.2对菊花遗传育种的意义本研究成果对菊花遗传育种具有重要的指导意义。鉴定出的与主要园艺性状相关的候选基因，为菊花分子标记辅助育种提供了丰富的基因资源和分子标记。在育种过程中，可以利用与优良性状紧密关联的SNP标记，对菊花品种进行早期筛选和鉴定，快速准确地选择具有目标性状的个体，提高育种效率，缩短育种周期。这些候选基因的功能验证结果，为通过基因工程手段改良菊花品种提供了理论依据。可以将具有重要功能的基因导入菊花品种中，创造新的种质资源，丰富菊花的遗传多样性，培育出具有更优良园艺性状的新品种。深入了解菊花主要园艺性状的遗传机制，有助于育种者更好地理解性状的遗传规律，合理选配亲本，制定更科学的育种策略，推动菊花遗传育种工作的发展，为菊花产业的可持续发展提供有力的支持。八、参考文献[1]陈俊愉，周家琪。中国菊花品种图志[M].北京：中国林业出版社，1993.[2]王四清，戴思兰，张启翔。菊花起源的研究进展[J].园艺学报，2001,28(2):173-178.[3]刘春迎，张莉俊，王雁。菊花种质资源研究进展[J].植物遗传资源学报，2009,10(2):307-313.[4]高亦珂，王雁，刘春迎，等。菊花品种数量分类研究[J].北京林业大学学报，2004,26(4):75-79.[5]张启翔，潘会堂。观赏植物遗传育种学[M].北京：中国林业出版社，2007.[6]王彩云，包满珠。菊花花色形成的生理生化基础[J].华中农业大学学报，1998,17(1):72-76.[7]魏宇昆，孙自然，刘改秀，等。菊花花色遗传规律的初步研究[J].北方园艺，2010(23):117-119.[8]周琳，郭先锋，张启翔，等。菊花查尔酮合成酶基因的克隆及表达分析[J].园艺学报，2011,38(2):361-368.[9]刘改秀，孙自然，魏宇昆，等。菊花MYB转录因子基因的克隆及表达分析[J].西北植物学报，2011,31(3):451-457.[10]张秀新，王四清，包满珠，等。菊花品种资源的数量分类研究[J].园艺学报，1999,26(5):305-310.[11]赵宏波，陈发棣，房伟民，等。切花菊品种花部性状的遗传分析[J].园艺学报，2008,35(8):1203-1208.[12]任志雨，陈发棣，房伟民，等。切花菊主要观赏性状的遗传分析[J].南京农业大学学报，2011,34(2):55-60.[13]张超，陈发棣，房伟民，等。菊花MADS-box基因家族全基因组鉴定及表达分析[J].园艺学报，2014,41(10):2013-2024.[14]孙红梅，李天来，李云飞，等。温度和光照对菊花花期及叶片生理特性的影响[J].应用生态学报，2006,17(1):45-48.[15]魏宇昆，孙自然，刘改秀，等。光周期和温度对菊花花期调控的研究进展[J].北方园艺，2010(22):217-220.[16]赵宏波，陈发棣，房伟民，等。菊花CO基因的克隆及表达分析[J].西北植物学报，2008,28(10):1949-1955.[17]赵宏波，陈发棣，房伟民，等。菊花FT基因的克隆及表达分析[J].园艺学报，2009,36(1):97-102.[18]赵宏波，陈发棣，房伟民，等。菊花SOC1基因的克隆及表达分析[J].南京农业大学学报，2009,32(2):47-52.[19]王雁，刘春迎，彭镇华，等。菊花品种数量分类研究[J].林业科学研究，2004,17(1):1-7.[20]刘春迎，王雁，彭镇华，等。菊花品种的ISSR分析[J].林业科学研究，2004,17(3):329-333.[21]王利琳，梁宗锁，胡正海。菊花花瓣发育过程中细胞壁物质的变化[J].西北植物学报，2005,25(11):2188-2193.[22]胡可，陈发棣，房伟民，等。切花菊株型相关性状的遗传分析[J].园艺学报，2012,39(1):129-136.[23]田孟良，张启翔，潘会堂，等。植物数量性状基因定位的研究进展[J].分子植物育种，2005,3(3):421-426.[24]李宏伟，赵久然，王凤格，等。关联分析及其在植物遗传学研究中的应用[J].遗传，2010,32(7):676-685.[25]黄璐琦，袁媛，郭兰萍，等。全基因组关联分析在药用植物研究中的应用前景[J].中国中药杂志，2012,37(2):145-149.[26]任羽，王得元，张银东，等。分子标记及其在植物遗传育种中的应用[J].分子植物育种，2004,2(1):119-128.[27]朱玉君，赵久然，王凤格，等。玉米株高、穗位高的QTL定位及上位性分析[J].作物学报，2010,36(11):1839-1847.[28]陈绍江，汤继华，严建兵，等。玉米株高杂种优势的QTL分析[J].作物学报，2006,32(7):980-985.[29]宋献军，张改生，牛娜，等。普通小麦穗发芽抗性的QTL分析[J].作物学报，2008,34(3):409-416.[30]刘勇，李加纳，谌利，等。甘蓝型油菜株高和分枝习性的QTL分析[J].作物学报，2009,35(11):2000-2007.[31]赵宏波，陈发棣，房伟民，等。切花菊品种资源遗传多样性的SRAP分析[J].南京农业大学学报，2007,30(4):51-56.[32]任羽，王得元，张银东，等.SRAP标记及其在作物遗传育种中的应用[J].分子植物育种，2004,2(5):723-727.[33]房经贵，王晨，李勇，等。葡萄基因组SSR和EST-SSR信息分析[J].园艺学报，2009,36(1):79-86.[34]王晨，房经贵，陈俊昆，等。葡萄EST-SSR标记的开发与应用[J].园艺学报，2009,36(4):543-550.[35]王得元，张银东，李颖，等。辣椒SRAP-PCR反应体系的建立与优化[J].分子植物育种，2006,4(2):253-257.[36]李巧云，赵前程，王得元，等。利用SRAP标记分析辣椒种质资源的遗传多样性[J].分子植物育种，2007,5(4):545-550.[37]胡延吉，赵檀方，陈雨海，等。利用RFLP标记研究小麦骨干亲本的遗传贡献[J].遗传学报，2000,27(11):1022-1029.[38]王关林，方宏筠。植物基因工程原理与技术[M].北京：科学出版社，2002.[39]朱玉君，赵久然，王凤格，等。玉米株高、穗位高的QTL定位及上位性分析[J].作物学报，2010,36(11):1839-1847.[40]严建兵，汤继华，陈绍江，等。玉米株高和穗位高QTL的一致性和特异性分析[J].中国农业科学，2008,41(11):3423-3430.[2]王四清，戴思兰，张启翔。菊花起源的研究进展[J].园艺学报，2001,28(2):173-178.[3]刘春迎，张莉俊，王雁。菊花种质资源研究进展[J].植物遗传资源学报，2009,10(2):307-313.[4]高亦珂，王雁，刘春迎，等。菊花品种数量分类研究[J].北京林业大学学报，2004,26(4):75-79.[5]张启翔，潘会堂。观赏植物遗传育种学[M].北京：中国林业出版社，2007.[6]王彩云，包满珠。菊花花色形成的生理生化基础[J].华中农业大学学报，1998,17(1):72-76.[7]魏宇昆，孙自然，刘改秀，等。菊花花色遗传规律的初步研究[J].北方园艺，2010(23):117-119.[8]周琳，郭先锋，张启翔，等。菊花查尔酮合成酶基因的克隆及表达分析[J].园艺学报，2011,38(2):361-368.[9]刘改秀，孙自然，魏宇昆，等。菊花MYB转录因子基因的克隆及表达分析[J].西北植物学报，2011,31(3):451-457.[10]张秀新，王四清，包满珠，等。菊花品种资源的数量分类研究[J].园艺学报，1999,26(5):305-310.[11]赵宏波，陈发棣，房伟民，等。切花菊品种花部性状的遗传分析[J].园艺学报，2008,35(8):1203-1208.[12]任志雨，陈发棣，房伟民，等。切花菊主要观赏性状的遗传分析[J].南京农业大学学报，2011,34(2):55-60.[13]张超，陈发棣，房伟民，等。菊花MADS-box基因家族全基因组鉴定及表达分析[J].园艺学报，2014,41(10):2013-2024.[14]孙红梅，李天来，李云飞，等。温度和光照对菊花花期及叶片生理特性的影响[J].应用生态学报，2006,17(1):45-48.[15]魏宇昆，孙自然，刘改秀，等。光周期和温度对菊花花期调控的研究进展[J].北方园艺，2010(22):217-220.[16]赵宏波，陈发棣，房伟民，等。菊花CO基因的克隆及表达分析[J].西北植物学报，2008,28(10):1949-1955.[17]赵宏波，陈发棣，房伟民，等。菊花FT基因的克隆及表达分析[J].园艺学报，2009,36(1):97-102.[18]赵宏波，陈发棣，房伟民，等。菊花SOC1基因的克隆及表达分析[J].南京农业大学学报，2009,32(2):47-52.[19]王雁，刘春迎，彭镇华，等。菊花品种数量分类研究[J].林业科学研究，2004,17(1):1-7.[20]刘春迎，王雁，彭镇华，等。菊花品种的ISSR分析[J].林业科学研究，2004,17(3):329-333.[21]王利琳，梁宗锁，胡正海。菊花花瓣发育过程中细胞壁物质的变化[J].西北植物学报，2005,25(11):2188-2193.[22]胡可，陈发棣，房伟民，等。切花菊株型相关性状的遗传分析[J].园艺学报，2012,39(1):129-136.[23]田孟良，张启翔，潘会堂，等。植物数量性状基因定位的研究进展[J].分子植物育种，2005,3(3):421-426.[24]李宏伟，赵久然，王凤格，等。关联分析及其在植物遗传学研究中的应用[J].遗传，2010,32(7):676-685.[25]黄璐琦，袁媛，郭兰萍，等。全基因组关联分析在药用植物研究中的应用前景[J].中国中药杂志，2012,37(2):145-149.[26]任羽，王得元，张银东，等。分子标记及其在植物遗传育种中的应用[J].分子植物育种，2004,2(1):119-128.[27]朱玉君，赵久然，王凤格，等。玉米株高、穗位高的QTL定位及上位性分析[J].作物学报，2010,36(11):1839-1847.[28]陈绍江，汤继华，严建兵，等。玉米株高杂种优势的QTL分析[J].作物学报，2006,32(7):980-985.[29]宋献军，张改生，牛娜，等。普通小麦穗发芽抗性的QTL分析[J].作物学报，2008,34(3):409-416.[30]刘勇，李加纳，谌利，等。甘蓝型油菜株高和分枝习性的QTL分析[J].作物学报，2009,35(11):2000-2007.[31]赵宏波，陈发棣，房伟民，等。切花菊品种资源遗传多样性的SRAP分析[J].南京农业大学学报，2007,30(4):51-56.[32]任羽，王得元，张银东，等.SRAP标记及其在作物遗传育种中的应用[J].分子植物育种，2004,2(5):723-727.[33]房经贵，王晨，李勇，等。葡萄基因组SSR和EST-SSR信息分析[J].园艺学报，2009,36(1):79-86.[34]王晨，房经贵，陈俊昆，等。葡萄EST-SSR标记的开发与应用[J].园艺学报，2009,36(4):543-550.[35]王得元，张银东，李颖，等。辣椒SRAP-PCR反应体系的建立与优化[J].分子植物育种，2006,4(2):253-257.[36]李巧云，赵前程，王得元，等。利用SRAP标记分析辣椒种质资源的遗传多样性[J].分子植物育种，2007,5(4):545-550.[37]胡延吉，赵檀方，陈雨海，等。利用RFLP标记研究小麦骨干亲本的遗传贡献[J].遗传学报，2000,27(11):1022-1029.[38]王关林，方宏筠。植物基因工程原理与技术[M].北京：科学出版社，2002.[39]朱玉君，赵久然，王凤格，等。玉米株高、穗位高的QTL定位及上位性分析[J].作物学报，2010,36(11):1839-1847.[40]严建兵，汤继华，陈绍江，等。玉米株高和穗位高QTL的一致性和特异性分析[J].中国农业科学，2008,41(11):3423-3430.[3]刘春迎，张莉俊，王雁。菊花种质资源研究进展[J].植物遗传资源学报，2009,10(2):307-313.[4]高亦珂，王雁，刘春迎，等。菊花品种数量分类研究[J].北京林业大学学报，2004,26(4):75-79.[5]张启翔，潘会堂。观赏植物遗传育种学[M].北京：中国林业出版社，2007.[6]王彩云，包满珠。菊花花色形成的生理生化基础[J].华中农业大学学报，1998,17(1):72-76.[7]魏宇昆，孙自然，刘改秀，等。菊花花色遗传规律的初步研究[J].北方园艺，2010(23):117-119.[8]周琳，郭先锋，张启翔，等。菊花查尔酮合成酶基因的克隆及表达分析[J].园艺学报，2011,38(2):361-368.[9]刘改秀，孙自然，魏宇昆，等。菊花MYB转录因子基因的克隆及表达分析[J].西北植物学报，2011,31(3):451-457.[10]张秀新，王四清，包满珠，等。菊花品种资源的数量分类研究[J].园艺学报，1999,26(5):305-310.[11]赵宏波，陈发棣，房伟民，等。切花菊品种花部性状的遗传分析[J].园艺学报，2008,35(8):1203-1208.[12]任志雨，陈发棣，房伟民，等。切花菊主要观赏性状的遗传分析[J].南京农业大学学报，2011,34(2):55-60.[13]张超，陈发棣，房伟民，等。菊花MADS-box基因家族全基因组鉴定及表达分