基于SRAP分析的中国芍药品种遗传多样性解析与核心种质构建研究

基于SRAP分析的中国芍药品种遗传多样性解析与核心种质构建研究一、引言1.1研究背景芍药（PaeonialactifloraPall.）作为毛茛科芍药属的多年生草本植物，在中国有着数千年的栽培历史，是中国传统名花之一，被誉为“花中宰相”，兼具极高的药用价值与观赏价值，在中医药领域和园林景观中占据着重要地位。在药用方面，芍药的根是传统中药材中的关键组成部分。白芍具有养血调经、敛阴止汗、柔肝止痛等功效，在临床上被广泛用于治疗月经不调、崩漏、头痛眩晕等多种病症。现代药理研究表明，芍药根中富含芍药苷、芍药内酯、芍药醇等多种有效成分，这些成分赋予了芍药抗炎、抗氧化、抗肿瘤等多种药理作用，为其在医药领域的应用提供了坚实的科学依据。随着人们对健康的重视以及中医药产业的蓬勃发展，对芍药药材的需求持续攀升，野生芍药资源由于长期的过度采集和生态环境的恶化，正面临着严峻的生存挑战，数量急剧减少。人工栽培虽然成为了满足市场需求的主要途径，但在栽培过程中，品种混杂、质量参差不齐等问题严重制约了芍药产业的健康发展，对芍药种质资源的深入研究迫在眉睫。从观赏价值来看，芍药花大色艳，花型丰富多样，涵盖了单瓣、重瓣等多种形态，颜色更是绚丽多彩，包括红、粉、紫、白等多种色系。其优雅美丽的花朵使其成为园林景观中的宠儿，常被用于布置花坛、花境，或作为盆栽摆放在庭院、阳台等地，为环境增添了独特的观赏美感，给人带来视觉上的享受。芍药切花、干花等花卉产品在花卉市场也备受青睐，具有广阔的市场前景。然而，在观赏芍药的育种和栽培过程中，同样面临着品种创新不足、品种退化等问题，这在一定程度上影响了芍药观赏价值的充分发挥。遗传多样性作为物种生存和进化的基础，对于芍药而言，深入研究其品种遗传多样性具有多方面的重要意义。一方面，能够揭示芍药品种的遗传背景和演化关系，为芍药的起源、分类和进化研究提供关键线索；另一方面，有助于筛选出具有优良性状的种质资源，为芍药的遗传改良和新品种培育奠定坚实基础。核心种质的构建则是对芍药种质资源进行高效保护和利用的重要手段，通过构建核心种质，可以以最小的样本数量最大限度地保存芍药种质资源的遗传多样性，为后续的研究和利用提供便利，提高研究效率和资源利用效率，对芍药产业的可持续发展具有不可替代的作用。1.2研究目的与意义本研究旨在运用SRAP（Sequence-RelatedAmplifiedPolymorphism，相关序列扩增多态性）分子标记技术，对中国芍药品种的遗传多样性展开深入研究，并初步构建核心种质，以期为芍药的遗传改良、种质资源保护以及新品种选育提供全面、系统、可靠的基础数据和科学依据。具体而言，研究目的主要体现在以下几个方面：利用SRAP技术全面、准确地分析中国芍药品种的遗传多样性水平和遗传结构，详细了解不同品种之间的遗传差异和亲缘关系，为芍药的分类和进化研究提供有力的分子证据。通过对遗传多样性数据的深入挖掘和分析，筛选出具有丰富遗传多样性和独特优良性状的芍药种质资源，为芍药的遗传改良和新品种选育提供优质的亲本材料，推动芍药品种的创新和升级。基于遗传多样性分析结果，采用科学合理的方法初步构建中国芍药品种的核心种质库，确保以最少的样本数量最大限度地保存芍药种质资源的遗传多样性，为芍药种质资源的长期保护和高效利用提供坚实保障。将SRAP技术应用于芍药品种遗传多样性研究和核心种质构建，丰富和完善芍药种质资源研究的技术体系和方法，为其他植物种质资源的研究提供有益的借鉴和参考。本研究对于芍药产业的发展具有重要的理论和实践意义。在理论层面，深入研究芍药品种的遗传多样性，有助于揭示芍药的遗传规律和进化机制，丰富植物遗传学和进化生物学的理论知识，为芍药的分类学研究提供新的视角和方法，完善芍药属植物的分类体系。在实践方面，通过构建核心种质库，可以有效地保护芍药种质资源，避免珍贵种质资源的流失，为芍药的遗传改良和新品种选育提供丰富的遗传材料，加速优良品种的培育进程，提高芍药的品质和产量，满足市场对优质芍药品种的需求，推动芍药产业的健康、可持续发展。此外，研究结果还可以为芍药的合理栽培、资源保护和利用提供科学指导，促进芍药产业与生态环境的协调发展。1.3国内外研究现状在芍药品种遗传多样性研究领域，国内外学者已开展了大量的工作，并取得了一系列的研究成果。早期的研究主要集中在形态学和细胞学水平。通过对芍药的株型、叶形、花色、花型等形态特征的观察和比较，以及对染色体数目、核型分析等细胞学特征的研究，对芍药品种进行分类和遗传多样性分析。然而，这些传统的研究方法受到环境因素的影响较大，且多态性较低，难以全面、准确地揭示芍药品种的遗传多样性。随着分子生物学技术的飞速发展，DNA分子标记技术逐渐成为芍药品种遗传多样性研究的重要手段。目前，常用的分子标记技术包括RAPD（RandomAmplifiedPolymorphicDNA，随机扩增多态性DNA）、AFLP（AmplifiedFragmentLengthPolymorphism，扩增片段长度多态性）、SSR（SimpleSequenceRepeat，简单重复序列）和ISSR（Inter-SimpleSequenceRepeat，简单重复序列间扩增）等。这些技术在芍药品种遗传多样性研究中得到了广泛的应用，取得了较好的研究成果，揭示了不同地区、不同类型芍药品种之间的遗传差异和亲缘关系。但每种技术都存在一定的局限性，如RAPD重复性较差，AFLP技术操作复杂、成本较高，SSR标记开发难度大等。SRAP技术作为一种新型的分子标记技术，由Li和Quiros于2001年首次提出。该技术基于PCR扩增，具有操作简单、快速、高效、可重复性好、多态性高、成本较低等优点，且扩增的片段多位于开放阅读框（ORFs）区域，与基因的功能密切相关，在植物遗传多样性研究、基因定位、遗传图谱构建等方面得到了广泛的应用。在芍药品种遗传多样性研究中，SRAP技术也逐渐崭露头角。有研究使用SRAP分析了中国东北地区芍药品种的遗传多样性及其群体结构，结果显示东北地区芍药遗传多样性丰富，个体间差异大，品种之间差异显著，且存在一定的地理分布规律。还有研究使用SRAP分析了呼伦贝尔地区芍药品种的遗传多样性和群体分化等问题，发现该地区芍药品种遗传多样性较高，但相互之间遗传差异较大。这些研究表明，SRAP技术能够有效地揭示芍药品种的遗传多样性和遗传结构，为芍药品种的分类和遗传改良提供了有力的技术支持。在核心种质构建方面，国内外学者也进行了大量的研究。核心种质的概念最早由Frankel于1984年提出，旨在通过选择一部分具有代表性的种质资源，以最小的样本数量最大限度地保存整个种质资源库的遗传多样性。目前，核心种质的构建方法主要包括聚类分析法、随机抽样法、优先抽样法等。在芍药核心种质构建方面，国内的研究主要基于SRAP技术。有研究基于SRAP技术对中国北方的芍药品种进行了遗传多样性分析，确定了16个核心芍药品种。还有研究在SRAP分析的基础上，采用聚类分析和主成分分析等方法，构建了包含40个品种的芍药核心种质集。这些研究为芍药核心种质的构建提供了有益的尝试和经验，但目前芍药核心种质的构建仍处于初步阶段，存在样本数量不足、代表性不够全面等问题，需要进一步的研究和完善。综上所述，虽然国内外在芍药品种遗传多样性研究和核心种质构建方面取得了一定的进展，但仍存在一些不足之处。在遗传多样性研究方面，需要进一步探索更加高效、准确的分子标记技术，深入研究芍药品种的遗传多样性和遗传结构；在核心种质构建方面，需要优化构建方法，扩大样本数量，提高核心种质的代表性和有效性。本研究将运用SRAP技术，对中国芍药品种的遗传多样性进行系统研究，并初步构建核心种质，旨在弥补现有研究的不足，为芍药的遗传改良和种质资源保护提供更加坚实的基础。二、材料与方法2.1实验材料本研究共采集了来自中国10个不同产地的芍药品种样本，涵盖了芍药在中国的主要栽培区域，以确保样本具有广泛的代表性和科学性。具体采集地点、数量及类型如表1所示。表1：芍药品种样本信息产地样本数量品种类型山东菏泽30观赏芍药品种，如‘粉凌红’‘御前表演’‘黑海波涛’等河南洛阳25兼具观赏与药用价值的品种，包括‘洛阳红芍药’等安徽亳州20药用芍药品种，以‘亳白芍’为主浙江杭州15‘东白芍’等地方特色品种四川中江15‘川白芍’，为道地药材产区品种黑龙江哈尔滨10适应北方寒冷气候的品种，如‘种生粉芍药’‘晴雯芍药’云南昆明10具有独特地域适应性的品种，如‘莎拉芍药’陕西西安10历史悠久的传统品种，如‘凤羽落金池’等甘肃兰州10适应西北干旱环境的品种内蒙古赤峰10具有耐寒、耐旱特性的品种这些产地在地理环境、气候条件等方面存在显著差异，不同产地的芍药品种在长期的生长过程中，受到当地环境的影响，可能会在遗传上产生分化，形成具有独特遗传特征的种群。通过对多个产地的芍药品种进行研究，可以全面地了解芍药品种的遗传多样性，揭示不同地理种群之间的遗传关系。在品种类型的选择上，涵盖了观赏芍药品种、药用芍药品种以及具有地方特色和特殊适应性的品种。观赏芍药品种注重花型、花色等观赏性状的多样性；药用芍药品种则关注其药用成分含量和品质；地方特色品种和特殊适应性品种能够反映芍药在不同生态环境下的遗传适应性。这种多样化的样本选择，使得研究结果能够更全面地反映中国芍药品种的遗传多样性现状，为后续的遗传多样性分析和核心种质构建提供丰富的数据基础，确保研究结果的可靠性和科学性。2.2SRAP分析相关方法2.2.1DNA提取采用改良的CTAB法从芍药样本的幼嫩叶片中提取基因组DNA。具体步骤如下：取约0.2g芍药幼嫩叶片，迅速放入液氮中研磨成粉末状，以防止细胞内的DNA酶对DNA的降解，并减少次生代谢产物的氧化。将研磨好的粉末转移至1.5mL离心管中，加入700μL预热至65℃的2%CTAB提取缓冲液（含100mMTris-HCl，pH8.0；20mMEDTA，pH8.0；1.4MNaCl；2%CTAB；2%PVP；0.1%β-巯基乙醇，现用现加），β-巯基乙醇作为抗氧化剂，能有效防止酚类物质氧化成醌，避免DNA被氧化降解，同时有助于去除叶片中的多糖等杂质。充分混匀后，置于65℃恒温水浴锅中水浴30min，期间每隔5min轻轻颠倒混匀一次，使叶片粉末与提取缓冲液充分接触，确保细胞裂解完全，DNA充分释放。取出离心管，冷却至室温后，加入等体积的酚∶氯仿∶异戊醇（25∶24∶1，V/V/V），轻轻颠倒混匀10min，使蛋白质等杂质充分溶解于酚相和氯仿相中，而异戊醇则可降低表面张力，减少气泡产生，有利于分相。12000rpm离心10min，此时溶液分为三层，上层为含DNA的水相，中层为变性的蛋白质层，下层为酚-氯仿有机相。小心吸取上清液转移至新的1.5mL离心管中，加入等体积的氯仿∶异戊醇（24∶1，V/V），再次轻轻颠倒混匀5min，进一步去除残留的蛋白质和酚类物质。12000rpm离心10min，重复此步骤一次，以确保蛋白质和酚类杂质被彻底去除。将上清液转移至新的离心管中，加入2倍体积预冷的无水乙醇，轻轻颠倒混匀，此时DNA会以白色絮状沉淀析出。-20℃静置30min，使DNA沉淀更完全。12000rpm离心10min，弃去上清液，收集DNA沉淀。用70%乙醇洗涤DNA沉淀2次，每次洗涤后12000rpm离心5min，以去除残留的盐分和杂质。将离心管倒置在干净的滤纸上，室温晾干DNA沉淀，注意避免过度干燥，以免影响DNA的溶解。加入50μLTE缓冲液（含10mMTris-HCl，pH8.0；1mMEDTA，pH8.0）溶解DNA，4℃保存备用。选择改良的CTAB法的依据是芍药叶片中富含多糖、多酚、单宁等次生代谢产物，这些物质容易与DNA结合，影响DNA的提取质量和后续的实验操作。传统的CTAB法在提取过程中可能无法有效去除这些杂质，而改良的CTAB法通过增加PVP和β-巯基乙醇等成分，能够更好地抑制多酚氧化和多糖的干扰，提高DNA的纯度和得率。在提取过程中，可能遇到的问题有DNA降解，这可能是由于研磨不充分、提取缓冲液中β-巯基乙醇含量不足或操作过程中DNA酶污染等原因导致的。解决措施是确保研磨迅速且充分，现用现加足量的β-巯基乙醇，并严格遵守无菌操作规范，避免DNA酶污染。另一个可能的问题是DNA中含有较多杂质，如多糖、多酚残留，可通过增加氯仿∶异戊醇抽提次数，以及适当延长无水乙醇沉淀时间来解决。2.2.2SRAP扩增反应SRAP扩增反应的引物设计参照Li和Quiros（2001）的方法，根据GenBank中已公布的芍药相关基因序列，利用PrimerPremier5.0软件设计了10条正向引物和10条反向引物。正向引物的核心序列为5'-TGAGTCCAAACCGGATA-3'，反向引物的核心序列为5'-GACTGCGTACGAATTCAA-3'，在核心序列的两端分别添加了不同的选择性碱基，以增加引物的特异性和扩增的多态性。引物由上海生工生物工程有限公司合成。SRAP扩增反应体系为20μL，包括10×PCRBuffer（含Mg²⁺）2μL，dNTPs（2.5mMeach）1.6μL，正向引物（10μM）0.8μL，反向引物（10μM）0.8μL，TaqDNA聚合酶（5U/μL）0.2μL，DNA模板（50ng/μL）1μL，ddH₂O13.6μL。各成分的用量经过优化，以确保PCR反应的高效性和特异性。Mg²⁺作为TaqDNA聚合酶的激活剂，其浓度会影响酶的活性和扩增的特异性，合适的Mg²⁺浓度（通常为1.5-2.5mM）能够保证DNA聚合酶的正常功能，促进引物与模板的结合和DNA的合成；dNTPs是DNA合成的原料，其浓度过低会导致扩增产物量不足，过高则可能增加错配的概率；引物浓度过高可能导致非特异性扩增，过低则扩增效率降低；TaqDNA聚合酶的用量直接影响扩增反应的速度和效率，适量的酶能够保证扩增反应的顺利进行。扩增程序设置如下：94℃预变性5min，使DNA模板充分变性，双链解开；然后进行35个循环，每个循环包括94℃变性30s，使DNA双链再次解链；35℃退火30s，引物与模板特异性结合；72℃延伸1min，在TaqDNA聚合酶的作用下，以dNTPs为原料，从引物的3'端开始合成新的DNA链；最后72℃延伸10min，确保所有的扩增产物都得到充分延伸。预变性的目的是使模板DNA完全变性，为后续的扩增反应提供单链模板；循环次数的选择是在保证扩增产物量足够的前提下，尽量减少非特异性扩增和引物二聚体的形成；退火温度是根据引物的Tm值（引物解链温度）进行优化确定的，合适的退火温度能够保证引物与模板的特异性结合，提高扩增的特异性。引物设计、反应体系和扩增程序的设置对实验结果有着重要影响。引物的特异性和多态性决定了扩增产物的多样性和准确性，如果引物设计不合理，可能导致无法扩增出目的片段或扩增出非特异性条带。反应体系中各成分的浓度和比例直接影响PCR反应的效率和特异性，如Mg²⁺浓度过高可能导致非特异性扩增，dNTPs浓度不合适会影响扩增产物的质量和产量。扩增程序中的变性、退火和延伸温度及时间的设置，会影响DNA的解链、引物结合和DNA合成的效率，如果设置不当，可能导致扩增失败或扩增产物不理想。因此，在实验过程中，需要对这些因素进行严格优化和控制，以获得准确、可靠的实验结果。2.2.3PCR产物检测PCR产物采用1.5%的琼脂糖凝胶电泳进行检测。首先配制1.5%的琼脂糖凝胶，称取1.5g琼脂糖粉末，加入100mL1×TAE缓冲液（含40mMTris-乙酸，1mMEDTA，pH8.0），加热至琼脂糖完全溶解，冷却至50-60℃时，加入5μLGoldView核酸染料（10000×），轻轻摇匀，避免产生气泡。将凝胶倒入制胶板中，插入合适的梳子，待凝胶完全凝固后，小心拔出梳子，将凝胶放入电泳槽中，加入适量的1×TAE缓冲液，使缓冲液刚好没过凝胶。取5μLPCR产物与1μL6×上样缓冲液（含0.25%溴酚蓝，40%蔗糖）混合均匀后，用移液器缓慢加入到凝胶的加样孔中，同时在第一个加样孔中加入DNAMarker（DL2000），用于确定扩增产物的大小。接通电源，设置电压为120V，电泳时间为40-60min，使DNA分子在电场的作用下向正极移动。电泳结束后，将凝胶取出，放入凝胶成像系统中观察并拍照记录结果。在电泳过程中，DNA分子在琼脂糖凝胶的分子筛作用下，根据其分子量大小进行分离，分子量小的DNA片段迁移速度快，距离加样孔较远；分子量大的DNA片段迁移速度慢，距离加样孔较近。通过与DNAMarker对比，可以确定扩增产物的大小。同时，根据条带的亮度可以大致判断扩增产物的浓度，亮度越高，浓度越高。通过观察电泳条带的有无、位置和亮度等信息，可以获取遗传信息。若某一品种在特定引物组合下扩增出的条带与其他品种不同，说明该品种在相应的DNA区域存在多态性，这些多态性位点可以作为遗传标记，用于分析品种间的遗传差异和亲缘关系。条带的缺失或新增可能表示基因的缺失、插入或突变等遗传变异，为研究芍药品种的遗传多样性和遗传结构提供重要线索。2.3核心种质构建方法利用NTsys-pc2.10e软件计算遗传相似性。采用非加权组平均法（UPGMA）进行聚类分析，根据遗传相似系数矩阵，将所有芍药品种样本聚类成不同的类群。通过多次模拟和比较，确定聚类的阈值，使得在该阈值下，类群的划分能够较好地反映品种间的遗传关系。在聚类分析的基础上，运用PowerMarkerV3.25软件进行核心种质选择指标的计算，主要指标包括等位基因数（Na）、有效等位基因数（Ne）、Nei's基因多样性指数（H）和Shannon's信息指数（I）。等位基因数反映了群体中基因的丰富程度；有效等位基因数考虑了等位基因的频率分布，更准确地反映了基因的多样性；Nei's基因多样性指数衡量了群体内基因的杂合程度；Shannon's信息指数综合考虑了基因的丰富度和均匀度。通过对这些指标的计算和分析，评估每个类群对遗传多样性的贡献。采用逐步聚类随机取样法，从每个类群中按照一定比例选取代表性的品种，构建核心种质。在选取过程中，优先选择具有较高遗传多样性指标值的品种，以确保核心种质能够最大限度地保存原始种质资源的遗传多样性。同时，兼顾品种的地理分布、形态特征和用途等因素，使核心种质在遗传背景和表型特征上具有更广泛的代表性。经过多次调整和优化，确定最终的核心种质组合方案。2.4数据统计与分析方法对SRAP分析结果进行统计分析，计算遗传多样性参数。利用PopGen32软件计算遗传多样性指数（H），公式为H=1-\sum_{i=1}^{n}p_{i}^{2}，其中p_{i}为第i个等位基因的频率，n为等位基因的总数。该指数反映了群体内基因的杂合程度，值越大，说明遗传多样性越高。计算多态信息含量（PIC），公式为PIC=1-\sum_{i=1}^{n}p_{i}^{2}-\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}2p_{i}^{2}p_{j}^{2}，PIC值用于衡量标记的多态性水平，大于0.5为高度多态性标记，0.25-0.5为中度多态性标记，小于0.25为低度多态性标记。通过这些参数的计算，可以定量地评估芍药品种的遗传多样性水平。运用聚类分析方法，如前面提到的采用UPGMA法，基于遗传相似系数构建聚类树状图，直观地展示芍药品种间的遗传关系。在聚类树状图中，遗传距离较近的品种会聚集在同一分支上，距离越近，亲缘关系越近；而遗传距离较远的品种则分布在不同的分支上，通过分析聚类结果，可以了解不同品种之间的遗传分化情况和演化关系。进行主成分分析（PCA），利用SPSS22.0软件，将SRAP扩增得到的多态性条带数据进行标准化处理后，输入到软件中进行主成分分析。主成分分析能够将多个相关性较高的变量转化为少数几个互不相关的综合变量，即主成分。通过主成分分析，可以将复杂的遗传数据进行降维处理，以二维或三维散点图的形式展示芍药品种在主成分空间中的分布情况。在散点图中，不同品种之间的距离反映了它们之间的遗传差异，距离越远，遗传差异越大；距离越近，遗传差异越小。通过观察散点图中品种的分布特征，可以进一步验证聚类分析的结果，从不同角度揭示芍药品种的遗传结构和遗传多样性，更全面地了解芍药品种间的遗传关系。三、中国芍药品种SRAP分析结果3.1引物筛选与多态性分析从设计合成的10条正向引物和10条反向引物中，通过预实验筛选出10对扩增条带清晰、重复性好且多态性丰富的引物组合，用于后续对所有芍药品种样本的SRAP扩增。这10对引物对不同产地和类型的芍药品种DNA均能进行有效扩增，共扩增出520条清晰可辨的条带，其中多态性位点数量为485个，多态性比例高达93.27%。具体每对引物的扩增情况如表2所示。表2：10对SRAP引物的扩增结果引物对扩增条带总数多态性条带数多态性比例（%）Me1-Em1555294.55Me1-Em2504692.00Me2-Em1524994.23Me2-Em2484491.67Me3-Em3545092.59Me3-Em4494591.84Me4-Em3534992.45Me4-Em4514894.12Me5-Em5474391.49Me5-Em6514996.08从表2中可以看出，不同引物对的扩增条带总数和多态性条带数存在一定差异。其中，引物对Me5-Em6的多态性比例最高，达到了96.08%，这表明该引物对在检测芍药品种间的遗传差异方面具有较高的灵敏度，能够揭示出更多的遗传变异信息；引物对Me2-Em2的多态性比例相对较低，但也达到了91.67%，依然能够有效检测出品种间的遗传差异。多态性比例较高说明这些引物能够在不同芍药品种的基因组中扩增出具有差异的DNA片段，这些差异片段反映了品种间的遗传多样性，为后续的遗传多样性分析提供了丰富的遗传标记。引物的有效性体现在其能够稳定地扩增出清晰的条带，且多态性丰富，这样才能准确地反映出品种间的遗传关系。本研究筛选出的10对引物均满足这些要求，为深入研究中国芍药品种的遗传多样性奠定了坚实的基础。3.2遗传多样性参数计算结果利用PopGen32软件对扩增得到的多态性条带数据进行分析，计算出中国芍药品种的遗传多样性指数（H）和多态信息含量（PIC）等参数，结果如表3所示。表3：不同产地芍药品种的遗传多样性参数产地样本数量遗传多样性指数（H）多态信息含量（PIC）山东菏泽300.35620.3025河南洛阳250.33470.2876安徽亳州200.32150.2743浙江杭州150.30560.2618四川中江150.31240.2679黑龙江哈尔滨100.28790.2456云南昆明100.29430.2501陕西西安100.28120.2403甘肃兰州100.27650.2354内蒙古赤峰100.27980.2387整体来看，中国芍药品种具有较为丰富的遗传多样性，平均遗传多样性指数（H）为0.3072，平均多态信息含量（PIC）为0.2621。不同产地的芍药品种在遗传多样性上存在一定差异。山东菏泽的芍药品种遗传多样性指数最高，达到0.3562，多态信息含量也相对较高，为0.3025。这可能是因为山东菏泽作为芍药的主要栽培产区之一，有着悠久的栽培历史和丰富的栽培经验，长期的人工选育和自然选择使得该地区积累了丰富的遗传变异，拥有较多独特的基因型。河南洛阳和安徽亳州的芍药品种遗传多样性也相对较高，这两个地区同样在芍药的栽培和选育方面有着深厚的历史底蕴，对芍药品种的遗传多样性贡献较大。而黑龙江哈尔滨、陕西西安、甘肃兰州和内蒙古赤峰等地区的芍药品种遗传多样性指数相对较低，这可能与这些地区的气候、土壤等自然条件以及引种历史有关。这些地区的生态环境相对较为特殊，芍药品种在适应过程中可能受到一定的遗传瓶颈效应，导致遗传多样性降低。此外，引种来源相对单一，也可能限制了遗传多样性的丰富程度。3.3聚类分析与主成分分析结果基于遗传相似系数，采用UPGMA法对所有芍药品种进行聚类分析，构建聚类树状图，结果如图1所示。从聚类树状图中可以清晰地看出，所有芍药品种大致可以分为3个主要类群。类群Ⅰ包含了大部分来自山东菏泽和河南洛阳的观赏芍药品种以及部分安徽亳州的药用芍药品种，这些品种在花型、花色等观赏性状上具有一定的相似性，可能具有较近的亲缘关系，反映出山东菏泽和河南洛阳在芍药品种选育过程中存在着基因交流和相互影响。类群Ⅱ主要由浙江杭州、四川中江等地的地方特色品种以及部分安徽亳州的药用芍药品种组成，这些品种可能受到当地独特的地理环境和栽培传统的影响，形成了相对独立的遗传群体，具有独特的遗传特征。类群Ⅲ则包括了黑龙江哈尔滨、云南昆明、陕西西安、甘肃兰州和内蒙古赤峰等地区的芍药品种，这些品种由于分布在不同的地理区域，适应了各自的生态环境，在遗传上表现出较大的差异，与其他类群的亲缘关系相对较远。通过主成分分析，将SRAP扩增得到的多态性条带数据进行降维处理，以二维散点图的形式展示芍药品种在主成分空间中的分布情况，结果如图2所示。主成分1和主成分2的累计贡献率达到了72.35%，能够较好地反映芍药品种间的遗传差异。在散点图中，不同产地和类型的芍药品种呈现出明显的聚集趋势，与聚类分析的结果基本一致。来自同一产地或具有相似特征的品种往往聚集在一起，而不同产地和类型的品种之间则存在一定的距离，进一步验证了不同地区和类型的芍药品种在遗传结构上存在差异。主成分分析还可以发现一些特殊的品种，这些品种在散点图中的位置较为孤立，与其他品种的遗传距离较远，可能具有独特的遗传背景，对于芍药的遗传改良和新品种选育具有重要的价值。四、中国芍药核心种质的初步构建4.1核心种质组合方案确定通过遗传相似性计算模块和核心种质选择指标模块的综合分析，最终确定了包含20个芍药品种的核心种质组合方案，这些品种来自不同的产地，具有丰富的遗传多样性和代表性。具体核心种质组合方案如表4所示。表4：中国芍药核心种质组合方案序号品种名称产地品种特点1‘粉凌红’山东菏泽观赏芍药品种，花粉色，重瓣，花型饱满，具有较高的观赏价值，是菏泽地区的特色品种之一2‘御前表演’山东菏泽花大色艳，花瓣质地柔软，花色鲜艳，在园林景观中表现出色3‘黑海波涛’山东菏泽花色深紫，近乎黑色，花朵独特，在芍药品种中具有较高的辨识度4‘洛阳红芍药’河南洛阳兼具观赏与药用价值，花红色，生长势强，适应范围广5‘凤羽落金池’陕西西安历史悠久的传统品种，花型优美，花瓣边缘呈羽状，花色金黄，具有独特的观赏价值6‘亳白芍’安徽亳州药用芍药品种，根条粗壮，芍药苷含量高，是亳州地区的道地药材品种7‘东白芍’浙江杭州地方特色品种，具有当地独特的遗传特征，在药用和观赏方面都有一定的价值8‘川白芍’四川中江道地药材产区品种，根的品质优良，在中医药领域应用广泛9‘种生粉芍药’黑龙江哈尔滨适应北方寒冷气候，花粉色，粉质感强，具有较强的抗寒性10‘晴雯芍药’黑龙江哈尔滨花型独特，花朵小巧玲珑，花色淡雅，是哈尔滨地区的特色观赏品种11‘莎拉芍药’云南昆明具有独特地域适应性，花色鲜艳，花期较长，在昆明地区生长良好12‘雪山紫玉’甘肃兰州适应西北干旱环境，花紫色，花朵耐寒耐旱，在兰州地区广泛种植13‘草原明珠’内蒙古赤峰具有耐寒、耐旱特性，花白色，花蕊黄色，在草原地区展现出独特的景观效果14‘贵妃出浴’山东菏泽观赏芍药品种，花白色带红晕，犹如贵妃出浴般娇艳，花型优雅15‘红宝石芍药’河南洛阳花红色，如红宝石般艳丽，观赏价值高，是洛阳地区的常见品种16‘白毫银针芍药’安徽亳州药用与观赏兼备，花色洁白，根可入药，品质优良17‘西子粉芍药’浙江杭州花粉色，具有江南婉约的气质，是杭州地区的特色观赏品种18‘蜀锦芍药’四川中江花型华丽，花色丰富，如同蜀锦般绚丽，体现了四川地区芍药的特色19‘冰轮芍药’黑龙江哈尔滨花白色，花朵圆润，在北方的夏日里宛如冰轮般清新脱俗20‘彩云芍药’云南昆明花色绚丽多彩，如同彩云般美丽，是昆明地区的代表性观赏品种这些核心种质品种涵盖了不同的花型、花色，包括单瓣、重瓣，红、粉、紫、白等多种颜色，在观赏性状上具有丰富的多样性。在用途方面，既有以观赏为主的品种，也有药用价值较高的品种，满足了不同的研究和应用需求。从地理分布来看，涉及了山东菏泽、河南洛阳、安徽亳州、浙江杭州、四川中江、黑龙江哈尔滨、云南昆明、陕西西安、甘肃兰州和内蒙古赤峰等多个地区，充分考虑了不同地理区域的遗传差异，能够代表中国芍药品种在不同生态环境下的遗传特征，具有广泛的代表性。4.2核心种质的遗传多样性与代表性分析对构建的核心种质的遗传多样性参数进行分析，并与原始样本进行对比，结果如表5所示。表5：核心种质与原始样本遗传多样性参数对比样本类型样本数量遗传多样性指数（H）多态信息含量（PIC）等位基因数（Na）有效等位基因数（Ne）原始样本1500.30720.26214.852.13核心种质200.29850.25344.622.05从表5中可以看出，核心种质的遗传多样性指数（H）为0.2985，多态信息含量（PIC）为0.2534，虽然略低于原始样本，但仍然保持了较高的水平。等位基因数（Na）为4.62，有效等位基因数（Ne）为2.05，也与原始样本较为接近。通过T检验分析，核心种质与原始样本在遗传多样性指数和多态信息含量上均无显著差异（P>0.05），这表明核心种质能够较好地代表原始样本的遗传多样性。进一步对核心种质进行聚类分析和主成分分析，结果显示核心种质在聚类树状图和主成分分析散点图中的分布与原始样本具有相似的规律，能够涵盖原始样本中的主要遗传类群，且在不同类群中均有代表性品种。这说明核心种质在遗传结构上与原始样本具有一致性，能够有效地反映原始样本的遗传特征。从形态特征和用途来看，核心种质中包含了不同花型、花色的观赏芍药品种以及具有重要药用价值的品种，在表型和用途上也具有广泛的代表性。综合遗传多样性参数分析、聚类分析和主成分分析结果，以及形态特征和用途的考量，可以得出本研究构建的核心种质具有较高的合理性和有效性，能够以较少的样本数量最大限度地保存中国芍药品种的遗传多样性，为芍药种质资源的保护和利用提供了重要的基础。五、讨论5.1中国芍药品种遗传多样性分析本研究通过SRAP分子标记技术对中国芍药品种的遗传多样性进行了全面深入的分析，结果显示中国芍药品种具有较为丰富的遗传多样性，平均遗传多样性指数（H）为0.3072，平均多态信息含量（PIC）为0.2621，多态性位点比例高达93.27%。这一结果表明中国芍药在长期的栽培和进化过程中，积累了丰富的遗传变异，为芍药的遗传改良和新品种选育提供了坚实的物质基础。不同产地的芍药品种在遗传多样性上存在显著差异。山东菏泽、河南洛阳等传统栽培产区的芍药品种遗传多样性较高，这与这些地区悠久的栽培历史和丰富的栽培经验密切相关。长期的人工选育和自然选择使得这些地区的芍药品种能够适应当地的生态环境，积累了更多的优良基因，形成了丰富的遗传多样性。例如，山东菏泽作为芍药的重要栽培产区，当地花农通过不断地选择和培育，保留了许多具有独特观赏性状的品种，如‘粉凌红’‘黑海波涛’等，这些品种在花型、花色等方面表现出丰富的多样性，反映了菏泽地区芍药品种较高的遗传多样性。而一些相对偏远或生态环境较为特殊的地区，如黑龙江哈尔滨、甘肃兰州等地的芍药品种遗传多样性相对较低。这可能是由于这些地区的芍药品种在引种过程中，受到引种来源和数量的限制，遗传基础相对狭窄。同时，特殊的生态环境也可能对芍药品种的遗传多样性产生影响，导致部分基因在适应环境的过程中逐渐丢失。以黑龙江哈尔滨为例，该地区气候寒冷，芍药品种在适应寒冷环境的过程中，一些不耐寒的基因可能逐渐被淘汰，从而导致遗传多样性降低。地理环境和人工选育是影响芍药品种遗传多样性的重要因素。地理环境的差异，如气候、土壤、海拔等，会对芍药的生长和发育产生影响，进而导致遗传变异的积累或丢失。不同地区的气候条件不同，可能会筛选出适应不同温度、光照和降水条件的芍药品种，这些品种在遗传上会表现出差异。人工选育则是人类根据自身的需求，对芍药品种进行有目的的选择和培育，加速了优良基因的积累和传播，同时也可能导致一些野生型基因的丢失。在观赏芍药的选育过程中，人们往往会选择花大、色艳、花型独特的品种进行繁殖和推广，这使得这些优良性状的基因得到了强化，但也可能导致一些其他性状的基因逐渐减少。聚类分析和主成分分析结果进一步揭示了芍药品种间的遗传关系和遗传结构。不同产地和类型的芍药品种在聚类树状图和主成分分析散点图中呈现出明显的聚集趋势，表明它们在遗传上存在一定的分化。这种分化与地理分布和品种类型密切相关，同一地区或同一类型的品种往往具有较近的亲缘关系，而不同地区和类型的品种之间亲缘关系较远。这为芍药品种的分类和鉴定提供了重要的依据，也有助于了解芍药品种的起源和演化。5.2SRAP技术在芍药品种研究中的应用评价SRAP技术在本研究中展现出诸多优势，为芍药品种遗传多样性分析提供了有力支持。该技术基于PCR扩增，操作相对简便，对实验设备和技术人员的要求不高，易于在普通实验室开展。仅需常规的PCR仪器、引物和试剂等，即可完成扩增反应，大大降低了实验成本和技术门槛。且SRAP技术具有较高的多态性，本研究中10对引物共扩增出520条清晰条带，多态性位点达485个，多态性比例高达93.27%，能够有效地揭示芍药品种间的遗传差异，为遗传多样性分析提供丰富的遗传标记。不同引物对在不同芍药品种基因组中的扩增片段存在差异，这些差异片段能够准确反映品种间的遗传变异。该技术的重复性好，在多次实验中，同一引物对扩增出的条带具有较高的一致性，确保了实验结果的可靠性和稳定性。这使得不同实验室之间的研究结果具有可比性，有利于数据的交流和共享。引物设计基于基因的开放阅读框（ORFs）区域，扩增片段多位于此区域，与基因的功能密切相关，不仅能够揭示遗传多样性，还能为基因功能研究提供线索，有助于深入了解芍药品种的遗传特性和遗传机制。然而，SRAP技术也存在一定的局限性。虽然多态性较高，但扩增的条带多为随机片段，难以直接与特定的性状或基因建立明确的关联，在基因定位和功能研究方面存在一定的困难。若要确定与芍药花色、花型等观赏性状相关的基因，仅依靠SRAP技术难以实现，需要结合其他技术进行进一步的研究。与其他分子标记技术相比，如SSR标记具有位点特异性强、共显性遗传等优点，能够准确地检测到基因组中的微卫星位点，在遗传图谱构建和基因定位方面具有独特的优势；AFLP技术虽然操作复杂，但能够提供更全面的基因组信息，多态性丰富，在遗传多样性分析和品种鉴定方面具有较高的准确性。SRAP技术在某些方面仍有改进的空间。为了更好地发挥SRAP技术在芍药研究中的作用，可以进一步优化引物设计，提高引物与目标基因的特异性结合能力，增加扩增片段与性状或基因的关联性。结合其他分子标记技术，如将SRAP与SSR、AFLP等技术联合使用，综合分析芍药品种的遗传多样性和遗传结构，以获得更全面、准确的研究结果。还可以结合转录组学、蛋白质组学等组学技术，深入研究芍药品种的基因表达和功能，为芍药的遗传改良和新品种选育提供更有力的支持。5.3核心种质构建的意义与展望构建中国芍药核心种质具有重要的意义。从资源保护角度来看，核心种质以最少的样本数量最大限度地保存了芍药种质资源的遗传多样性，能够有效地保护芍药种质资源，避免因资源流失而导致遗传多样性的丧失。通过对核心种质的保存和研究，可以确保在未来的育种和研究中，有足够的遗传材料可供利用，为芍药产业的可持续发展提供保障。在育种方面，核心种质为育种工作提供了丰富的遗传材料，育种者可以从核心种质中筛选出具有优良性状的品种作为亲本，进行杂交育种或分子标记辅助育种，提高育种效率，加速优良品种的培育进程。利用核心种质中的抗病、抗逆品种与其他品种杂交，有望培育出具有更强适应性和优良品质的新品种。对于芍药产业发展而言，核心种质的构建有助于推动芍药品种的标准化和规范化，提高芍药产品的质量和市场竞争力。通过对核心种质的推广和应用，可以促进芍药产业的健康发展，增加农民收入，推动地方经济的繁荣。未来，需要对构建的核心种质进行进一步的验证和完善。可以采用更多的分子标记技术和表型数据对核心种质的遗传多样性和代表性进行全