基于微卫星DNA的甘肃马鹿遗传多样性解析及与青海马鹿亲缘关系探究

基于微卫星DNA的甘肃马鹿遗传多样性解析及与青海马鹿亲缘关系探究一、引言1.1研究背景马鹿（Cervuselaphus）作为国家二级重点保护野生动物，在生态系统中扮演着举足轻重的角色。作为大型林栖哺乳动物，马鹿是北方森林生态系统的重要组成部分，对维持生态平衡意义重大。其觅食行为能够影响植物的分布与生长，进而调节植被结构；同时，马鹿也是众多捕食者的重要猎物来源，在食物链中处于关键位置，其种群数量的变化会对整个生态系统的能量流动和物质循环产生深远影响，是衡量生态系统健康状态的重要指标。甘肃马鹿和青海马鹿是马鹿在中国的重要类群，然而，二者如今均面临着严峻的濒危困境。由于人类活动的强烈干扰，包括大规模的森林砍伐、草原开垦以及基础设施建设，导致甘肃马鹿和青海马鹿的栖息地遭到严重破坏与碎片化，适宜生存的空间大幅缩减。同时，非法捕猎现象屡禁不止，进一步加剧了其种群数量的减少。在过去的几十年间，甘肃马鹿和青海马鹿的种群数量均呈现出急剧下降的趋势，已被世界自然保护联盟（IUCN）列为濒危物种，生存状况岌岌可危。保护生物学领域高度重视遗传多样性，它是物种适应环境变化、维持种群生存与进化潜力的关键因素。对于濒危物种而言，深入了解其遗传多样性状况，能够为制定科学有效的保护策略提供坚实依据。微卫星DNA，又被称为短串联重复序列（Shorttandemrepeat，STR），以1-6个核苷酸为核心序列，首尾相连组成串联重复序列，重复数通常为10-20次。因其具有多态信息含量高、共显性遗传等诸多优点，微卫星DNA已成为评估物种遗传多样性和分析亲缘关系的重要分子标记，在品种资源分类及遗传多样性的评估和保存中发挥着关键作用。1.2研究目的与意义本研究旨在运用微卫星DNA技术，精准剖析甘肃马鹿的遗传多样性，深入探究其与青海马鹿的亲缘关系。通过对甘肃马鹿多个微卫星位点的细致分析，精确评估其遗传多样性水平，包括基因杂合度、多态信息含量等关键指标，全面揭示甘肃马鹿种群的遗传结构和遗传变异情况。同时，通过与青海马鹿的对比分析，明确二者之间的遗传距离和遗传分化程度，追溯它们的演化历史和遗传交流情况。从理论层面来看，深入了解甘肃马鹿和青海马鹿的遗传多样性，能够为马鹿属物种的分类学研究提供更为精确的遗传数据，有力推动物种进化理论的发展，加深对生物多样性形成和维持机制的认识。从实践意义而言，掌握甘肃马鹿的遗传多样性状况，对于制定科学有效的保护策略至关重要。了解种群的遗传结构，能够明确不同种群间的亲缘关系，避免近亲繁殖，维持种群的遗传健康，提高其生存和适应能力。而揭示甘肃马鹿与青海马鹿的亲缘关系，有助于确定马鹿类物种的系统发育关系，为物种的保护和繁殖提供关键的遗传学依据，对保护生物多样性和维护生态平衡具有重要意义。1.3国内外研究现状在国外，微卫星DNA技术在马鹿遗传多样性研究方面已取得了丰富成果。早期，欧洲学者通过对欧洲马鹿多个微卫星位点的分析，揭示了其种群内存在着显著的遗传分化，不同地理区域的马鹿种群在基因频率上存在明显差异，这种差异与地理隔离和生态环境的差异密切相关。北美地区对加拿大马鹿的研究也表明，微卫星标记能够有效区分不同亚种的马鹿，为其种群管理和保护提供了有力的遗传学依据。国内对于马鹿遗传多样性的研究起步相对较晚，但近年来发展迅速。新疆地区利用微卫星DNA技术对塔里木马鹿的研究发现，该种群具有较高的遗传变异度，遗传多样性相对丰富。黑龙江省完达山地区对东北马鹿的研究筛选出了8对多态性高、扩增效果良好的微卫星位点，成功进行了个体识别和性别鉴定，结果表明该地区马鹿具有较高的遗传多样性水平。在亲缘关系研究方面，国内外学者主要通过计算遗传距离、构建系统发育树等方法来推断马鹿不同种群或亚种之间的亲缘关系。国外研究通过对欧洲马鹿和亚洲马鹿部分种群的微卫星数据分析，初步揭示了它们之间的遗传关系，发现二者在演化过程中存在一定程度的基因交流。国内研究则针对中国境内不同地区的马鹿种群展开，通过对比不同种群马鹿的微卫星位点信息，探讨了它们之间的亲缘远近，为马鹿的分类和保护提供了重要参考。然而，针对甘肃马鹿遗传多样性的研究仍较为匮乏，对其与青海马鹿亲缘关系的研究更是鲜见报道。现有的研究大多集中在马鹿的某一特定种群或亚种，缺乏对甘肃马鹿和青海马鹿这两个重要类群的系统研究。同时，在研究方法上，虽然微卫星DNA技术已广泛应用，但如何进一步优化实验流程、提高数据准确性，以及如何将微卫星数据与其他分子标记或生态数据相结合，仍有待深入探索。此外，在保护实践中，如何将遗传多样性研究成果有效转化为保护策略，也是当前亟待解决的问题。二、研究区域与方法2.1研究区域概况2.1.1甘肃马鹿栖息地甘肃马鹿主要栖息于甘肃省祁连山林区，该区域地理位置独特，处于青藏高原、蒙新高原和黄土高原的交汇地带。其地理坐标大致为东经97°25′-103°36′，北纬36°63′-39°36′，东西跨度达1200千米，南北宽度在50-120千米之间，是甘肃马鹿重要的生存家园。祁连山林区地形复杂，海拔高度在3500-5000米之间，属于高山森林地带，具备多样化的生态环境，为甘肃马鹿提供了丰富的食物资源和适宜的栖息场所。该地区气候呈现出明显的垂直变化特征，年平均气温约为3.6℃，无霜期在90-120天左右。年降雨量处于350-500毫米的范围，年蒸发量则为1000-1200毫米，相对湿度维持在50%-70%之间。复杂的地形地貌和多样的气候条件，造就了祁连山林区丰富的植被类型。从低海拔到高海拔，依次分布着荒漠草原、山地草原、山地森林、高山灌丛和高山草甸等植被类型。甘肃马鹿主要以草本植物、灌木的嫩枝和树叶等为食，林区内丰富的植被资源为其提供了充足的食物来源。同时，茂密的森林和复杂的地形为甘肃马鹿提供了良好的隐蔽场所，有助于它们躲避天敌，保障自身安全。2.1.2青海马鹿栖息地青海马鹿主要分布于青海省的祁连山地区以及柴达木盆地周边区域。祁连山地区与甘肃马鹿的栖息地存在部分重叠，其地形同样以高山和山地为主，海拔高度多在3000米以上，地势起伏较大，山峰林立，峡谷纵横。柴达木盆地周边则呈现出独特的荒漠和半荒漠景观，地形相对较为平坦，但气候干旱，植被稀疏。青海马鹿栖息地的气候具有典型的高原大陆性特征，寒冷干燥，气温年较差和日较差都较大。年平均气温在-6℃-2℃之间，冬季漫长而寒冷，夏季短暂且凉爽。年降水量较少，大部分地区年降水量不足200毫米，且降水分布极不均匀，主要集中在夏季的几个月中。由于气候干旱，该地区植被以耐旱的草本植物、灌木和荒漠植物为主，如针茅、沙棘、红砂等。青海马鹿能够适应这种干旱的环境，以这些植物为食，维持自身的生存和繁衍。在冬季，当高山地区食物匮乏时，青海马鹿会向低海拔的山谷和向阳缓坡地带迁移，以寻找更多的食物和较为温暖的栖息环境。2.2实验材料2.2.1样本采集在甘肃马鹿的主要栖息地甘肃省祁连山林区，以及青海马鹿的分布区域青海省祁连山地区和柴达木盆地周边，进行样本采集工作。为确保样本的代表性，在不同的季节、不同的地理区域设置多个采样点。在祁连山林区，选择海拔高度不同、植被类型各异的区域，如高山森林、山地草原等，以涵盖甘肃马鹿的多样生存环境。在青海马鹿的栖息地，同样兼顾祁连山地区的高山环境和柴达木盆地周边的荒漠、半荒漠环境。对于马鹿样本的采集，优先选择血液样本。使用经过严格消毒的一次性注射器，在马鹿的颈静脉部位进行采血操作。为了降低对马鹿的惊扰和伤害，采血过程中尽量保持操作的轻柔与迅速。每只马鹿采集约5-10毫升的血液，将采集到的血液缓慢注入含有抗凝剂（如EDTA-K2）的无菌离心管中，轻轻颠倒混匀，防止血液凝固。在无法获取血液样本的情况下，采集马鹿的组织样本作为替代。使用经过灭菌处理的手术刀和镊子，在马鹿自然脱落的毛发根部、耳部皮肤或其他不易造成严重伤害的部位，采集少量的组织样本。例如，对于自然脱落的毛发，尽量选取带有毛囊的部分，因为毛囊中含有丰富的细胞，可提供足够的DNA。采集的组织样本立即放入无菌的离心管中，避免污染。在整个采样过程中，严格遵循野生动物保护法律法规，确保采样活动不会对马鹿的生存和繁衍造成不良影响。同时，详细记录每一个样本的采集信息，包括采集地点的经纬度、海拔高度、采集时间、马鹿的性别、年龄（如有可能判断）等，为后续的数据分析提供全面的背景资料。2.2.2样本保存与运输样本采集完成后，立即进行妥善的保存处理，以确保样本的质量不受影响。血液样本在加入抗凝剂并混匀后，迅速放置于冰盒中进行低温保存。冰盒中的冰块能够维持样本周围的低温环境，减缓DNA的降解速度。在野外采样点，尽量缩短样本在常温下的暴露时间，尽快将样本转移至冰盒中。组织样本同样放置于冰盒中保存。对于组织样本，为了进一步防止组织的腐败和DNA的降解，可以在样本中加入适量的DNA保护液，如RNAlater溶液。这种保护液能够稳定细胞内的核酸，抑制核酸酶的活性，从而有效保护DNA的完整性。在样本运输环节，采用专门的生物样本运输箱。运输箱内部配备有足量的冰袋，确保在运输过程中样本始终处于低温状态。冰袋的数量和放置位置经过精心设计，以保证运输箱内的温度均匀分布，避免局部温度过高导致样本受损。同时，运输箱具有良好的密封性和抗震性能，能够防止样本受到外界环境的污染和震动影响。将装有样本的离心管整齐放置在运输箱内的样本固定架上，确保样本在运输过程中不会发生晃动和碰撞。在运输过程中，密切关注运输箱内的温度变化，使用温度记录仪实时监测温度。一旦发现温度异常，及时采取措施进行调整，如补充冰袋或更换损坏的保温材料。如果样本需要进行长途运输，选择专业的冷链物流服务。这些物流公司具备完善的冷链运输体系，能够确保样本在整个运输过程中始终处于适宜的低温环境。在运输前，与物流公司进行充分沟通，明确样本的特殊运输要求，如温度范围、运输时间限制等。同时，为样本购买足额的运输保险，以降低运输过程中的风险。通过严格的样本保存和运输措施，确保采集到的马鹿样本能够以最佳状态送达实验室，为后续的实验分析提供可靠的材料基础。2.3实验方法2.3.1DNA提取采用磁珠法DNA提取试剂盒对采集的马鹿样本进行DNA提取。将血液样本或组织样本从冰盒中取出，放置在超净工作台中。对于血液样本，取200微升全血加入到含有细胞破碎缓冲液的离心管中，充分混匀，使血细胞完全裂解。对于组织样本，用无菌剪刀将其剪切成约1-2毫米的小块，放入含有细胞破碎缓冲液的离心管中，利用组织研磨仪进行充分研磨，确保细胞破碎完全。将裂解后的样本在离心机中以12000转/分钟的速度离心5分钟，使细胞碎片沉淀到离心管底部。小心吸取上清液，转移至新的离心管中。按照磁珠法DNA提取试剂盒的说明书，向上清液中加入适量的磁珠悬浮液，充分混匀，使磁珠与DNA分子特异性结合，形成DNA-磁珠复合物。将离心管放置在磁珠分离架上，静置2-3分钟，使磁珠复合物沉降到管底。此时，利用磁力架或磁力装置将磁珠吸附在管壁上，小心倒掉上清液，避免磁珠丢失。向含有磁珠的离心管中加入洗涤缓冲液，轻轻摇晃离心管，使磁珠复合物充分悬浮，以去除杂质。再次将离心管放置在磁珠分离架上，静置2-3分钟，倒掉上清液。重复洗涤步骤1-2次，确保磁珠表面的杂质被彻底清除。向离心管中加入适量的洗脱缓冲液，轻轻摇晃离心管，使磁珠复合物洗脱出DNA分子。将离心管放置在磁珠分离架上，静置2-3分钟，使磁珠沉降到管底。小心吸取上清液，转移至新的离心管中，即得到提取的DNA溶液。使用分光光度计检测提取的DNA的浓度和纯度。将DNA溶液稀释一定倍数后，加入到分光光度计的比色皿中，设置波长为260纳米和280纳米，分别测定吸光度值。根据公式计算DNA的浓度：DNA浓度（μg/μL）=A260×稀释倍数×50/1000。同时，通过A260/A280的比值来评估DNA的纯度，理想的比值范围在1.8-2.0之间。若比值偏离此范围，说明DNA可能存在蛋白质或RNA污染，需要进一步纯化。此外，还可以采用琼脂糖凝胶电泳对DNA的完整性进行检测。将DNA样品与适量的上样缓冲液混合后，加入到琼脂糖凝胶的加样孔中。在电泳槽中加入适量的电泳缓冲液，接通电源，设置电压和时间，进行电泳。电泳结束后，将凝胶放入含有核酸染料的溶液中染色15-20分钟，然后在凝胶成像系统下观察DNA条带。若DNA条带清晰、完整，无明显拖尾现象，说明提取的DNA质量良好，可用于后续实验。2.3.2微卫星位点选择通过查阅相关的马鹿遗传学研究文献，筛选出10个在马鹿遗传多样性研究中广泛应用且多态性较高的微卫星位点，分别为BM1824、BM1818、INRA005、INRA037、ETH225、ETH10、TPOX、CSRD247、SRCRSP1和RT1。这些位点在马鹿的不同染色体上分布较为均匀，能够全面反映马鹿的遗传信息。例如，BM1824位点位于马鹿的第1号染色体上，已被证实能够有效区分马鹿的不同种群；INRA005位点位于第5号染色体，在多个马鹿种群的遗传分析中表现出较高的多态性。选择这10个微卫星位点的依据主要有以下几点：首先，它们在马鹿属物种中具有较高的保守性，能够在不同地区的马鹿种群中稳定扩增，保证实验结果的可靠性。其次，已有研究表明这些位点的多态信息含量（PIC）较高，通常PIC值大于0.5，这意味着它们能够提供丰富的遗传信息，有助于准确评估遗传多样性。此外，这些位点在之前的马鹿遗传研究中得到了广泛应用，积累了大量的参考数据，便于与本研究结果进行对比和分析。2.3.3PCR扩增PCR扩增反应在25微升的反应体系中进行。其中，10×PCR缓冲液2.5微升，提供PCR反应所需的缓冲环境，维持反应体系的pH值稳定；2.5mmol/L的dNTP混合物2微升，作为DNA合成的原料，包含四种脱氧核苷酸（dATP、dTTP、dCTP、dGTP），为新合成的DNA链提供碱基；25mmol/L的MgCl21.5微升，Mg2+是TaqDNA聚合酶的激活剂，能够增强酶的活性，促进DNA合成；10μmol/L的上下游引物各0.5微升，引物是根据所选微卫星位点的侧翼序列设计合成的，能够特异性地结合到目标DNA区域，引导DNA的扩增；5U/μL的TaqDNA聚合酶0.2微升，负责催化DNA的合成反应，以模板DNA为指导，将dNTP逐个添加到引物的3'端，延伸DNA链；模板DNA1微升，含有待扩增的目标DNA序列；最后用ddH2O补足至25微升，确保反应体系的总体积符合要求。PCR扩增程序如下：首先进行预变性，将反应体系在95℃下加热5分钟，使模板DNA的双链完全解开，为后续的引物结合和DNA合成提供单链模板。然后进行35个循环的扩增，每个循环包括94℃变性30秒，使DNA双链再次解链；55℃退火30秒，引物与变性后的单链模板DNA特异性结合；72℃延伸30秒，在TaqDNA聚合酶的作用下，以dNTP为原料，从引物的3'端开始合成新的DNA链。循环结束后，再进行72℃延伸10分钟，确保所有新合成的DNA链都得到充分延伸，使扩增产物更加完整。在进行PCR扩增操作时，首先将PCR反应所需的各种试剂从冰箱中取出，放置在冰盒上融化，避免试剂温度过高导致活性降低。按照上述反应体系，依次准确吸取各种试剂，加入到无菌的PCR薄壁管中。在加样过程中，使用移液器时要确保量程准确，避免产生误差。每加完一种试剂，都要轻轻吹打混匀，确保试剂充分混合。加样完成后，将PCR薄壁管放入PCR仪中，按照设定好的扩增程序进行反应。在PCR仪运行过程中，要密切关注仪器的运行状态，如温度、时间等参数是否正常，确保扩增反应顺利进行。2.3.4琼脂糖凝胶电泳制备1.5%的琼脂糖凝胶。称取1.5克琼脂糖粉末，加入到100毫升1×TAE电泳缓冲液中，充分搅拌均匀。将混合液放入微波炉中加热，每隔15-20秒取出振荡一次，直至琼脂糖完全溶解，溶液变得澄清透明。待溶液冷却至50-60℃左右时，加入5微升核酸染料（如GoldView），轻轻摇匀，避免产生气泡。将制胶模具放置在水平台上，插入合适的梳子，然后将冷却后的琼脂糖溶液缓慢倒入模具中，厚度约为5-8毫米。待凝胶完全凝固后，小心拔出梳子，将凝胶放入电泳槽中，加入适量的1×TAE电泳缓冲液，使缓冲液刚好没过凝胶表面。取5微升PCR扩增产物与1微升6×上样缓冲液混合均匀，用移液器将混合液缓慢加入到凝胶的加样孔中。在加样过程中，要注意避免加样枪头触碰凝胶孔壁，防止样品泄漏。同时，在第一个加样孔中加入适量的DNA分子量标准（Marker），用于判断扩增产物的大小。接通电源，设置电压为120V，电泳时间为30-40分钟。在电泳过程中，DNA分子会在电场的作用下向正极移动，不同大小的DNA分子由于迁移速率不同，会在凝胶上形成不同的条带。电泳结束后，将凝胶从电泳槽中取出，放入凝胶成像系统中。打开凝胶成像系统的光源和相机，调整焦距和曝光时间，使凝胶条带清晰可见。观察凝胶上的条带分布情况，根据DNA分子量标准判断扩增产物的大小。若扩增产物的条带清晰、单一，且大小与预期相符，说明PCR扩增成功；若出现多条杂带或条带模糊不清，则可能存在非特异性扩增或扩增失败的情况，需要进一步分析原因并优化实验条件。2.3.5数据分析方法利用Popgene3.2软件对实验数据进行遗传多样性分析。通过该软件计算甘肃马鹿和青海马鹿在各个微卫星位点上的等位基因数（Na），它反映了种群中基因的丰富程度，等位基因数越多，说明基因的多样性越高。同时计算有效等位基因数（Ne），有效等位基因数考虑了等位基因的频率分布，能够更准确地反映基因的有效变异程度。观测杂合度（Ho）和期望杂合度（He）也是重要的分析指标，观测杂合度是指实际观察到的杂合子在种群中的比例，期望杂合度则是基于哈迪-温伯格平衡定律计算得到的理论杂合度，两者的差异可以反映种群中是否存在遗传偏离。多态信息含量（PIC）则综合考虑了等位基因数和基因频率，用于评估微卫星位点的多态性高低，PIC值越大，表明该位点提供的遗传信息越丰富。使用Arlequin3.5软件进行遗传分化分析。通过计算F-统计量中的Fst值来评估甘肃马鹿和青海马鹿种群之间的遗传分化程度。Fst值的范围在0-1之间，0表示两个种群之间没有遗传分化，完全同质；1表示两个种群之间完全分化，没有基因交流。一般认为，Fst值小于0.05表示种群间遗传分化程度很低，0.05-0.15表示中等程度的遗传分化，大于0.15则表示遗传分化程度较高。通过分析Fst值，可以明确甘肃马鹿和青海马鹿之间的遗传差异大小，为进一步探讨它们的亲缘关系提供数据支持。运用MEGA7.0软件构建系统发育树。采用邻接法（Neighbor-Joiningmethod），基于遗传距离矩阵构建甘肃马鹿和青海马鹿的系统发育树。在构建过程中，通过多次重复计算（如1000次bootstrap检验）来评估分支的可靠性。系统发育树能够直观地展示甘肃马鹿和青海马鹿个体或种群之间的亲缘关系远近，从进化的角度揭示它们的遗传联系。通过分析系统发育树的拓扑结构和分支长度，可以推断两个种群的演化历史和遗传分化时间，为深入了解马鹿的进化历程提供重要线索。三、甘肃马鹿遗传多样性分析3.1遗传多样性指标统计3.1.1等位基因数通过对甘肃马鹿样本在10个微卫星位点上的PCR扩增产物进行分析，统计出每个位点的等位基因数（Na）。结果显示，10个微卫星位点共检测到[X]个等位基因，每个位点的等位基因数范围为[最小值]-[最大值]。其中，位点BM1824的等位基因数最多，达到了[X]个；而位点TPOX的等位基因数最少，为[X]个。平均等位基因数为[X]个，表明甘肃马鹿在这些微卫星位点上具有一定程度的遗传变异。等位基因数是衡量遗传多样性的基础指标，较多的等位基因意味着种群具有更丰富的遗传信息，能够为种群的适应性和进化提供更多的遗传基础。例如，在其他鹿科动物的研究中，如梅花鹿的某些种群，平均等位基因数在[X]个左右，相比之下，甘肃马鹿的等位基因数处于[较高/较低/相当]水平。这可能与甘肃马鹿的种群历史、栖息地环境以及基因交流情况等多种因素有关。高海拔且复杂多样的祁连山林区为甘肃马鹿提供了多样的生态位，在长期的生存过程中，不同的环境压力可能促使马鹿种群在基因层面发生适应性变化，从而积累了较多的等位基因。同时，历史上甘肃马鹿可能与周边其他马鹿种群存在一定程度的基因交流，也增加了其等位基因的丰富度。3.1.2有效等位基因数有效等位基因数（Ne）能够更准确地反映基因的有效变异程度，它考虑了等位基因的频率分布情况。利用Popgene3.2软件计算得出，甘肃马鹿群体的有效等位基因数在各个位点上的范围为[最小值]-[最大值]，平均有效等位基因数为[X]个。与等位基因数相比，有效等位基因数普遍较低，这表明部分等位基因在种群中的频率较低，对遗传多样性的实际贡献相对较小。以位点BM1824为例，其等位基因数为[X]个，但有效等位基因数仅为[X]个，说明该位点上存在一些稀有等位基因，它们在种群中的频率较低，虽然增加了等位基因的数量，但对遗传多样性的实际影响相对有限。有效等位基因数的多少反映了种群中真正起作用的遗传变异水平，较低的有效等位基因数可能意味着种群在某些方面的遗传稳定性较高，但同时也可能降低了种群对环境变化的适应能力。在甘肃马鹿的保护和管理中，需要关注有效等位基因数的变化，避免因遗传漂变等因素导致稀有等位基因的丢失，从而维持种群的遗传多样性。3.1.3基因多样性基因多样性（GD）是衡量种群遗传多样性的重要指标之一，它反映了种群中基因的杂合程度和变异水平。经计算，甘肃马鹿在10个微卫星位点上的基因多样性范围为[最小值]-[最大值]，平均基因多样性为[X]。较高的基因多样性表明甘肃马鹿种群具有较为丰富的遗传变异，这有利于种群在面对环境变化时保持适应性。例如，当栖息地的气候发生变化或出现新的疾病时，丰富的遗传多样性使得种群中可能存在具有适应这些变化能力的个体，从而增加了种群生存和繁衍的机会。与其他地区马鹿种群的基因多样性相比，甘肃马鹿的基因多样性处于[较高/较低/相当]水平。如新疆地区的塔里木马鹿，其基因多样性为[X]，与甘肃马鹿的基因多样性存在一定差异。这种差异可能是由于地理隔离、生态环境差异以及人类活动影响程度不同等多种因素造成的。甘肃马鹿栖息地的生态环境相对较为复杂，可能促进了基因的多样化；而塔里木马鹿可能由于长期处于相对孤立的地理环境中，基因交流相对较少，导致基因多样性水平有所不同。3.1.4杂合度杂合度分为观察杂合度（Ho）和期望杂合度（He）。观察杂合度是指实际观察到的杂合子在种群中的比例，期望杂合度则是基于哈迪-温伯格平衡定律计算得到的理论杂合度。统计结果显示，甘肃马鹿的观察杂合度范围为[最小值]-[最大值]，平均观察杂合度为[X]；期望杂合度范围为[最小值]-[最大值]，平均期望杂合度为[X]。在大多数位点上，观察杂合度低于期望杂合度，这表明甘肃马鹿种群中存在一定程度的近交现象或杂合子缺失。以位点INRA005为例，其观察杂合度为[X]，期望杂合度为[X]，观察杂合度明显低于期望杂合度。近交可能导致有害隐性基因的纯合表达，增加种群中遗传疾病的发生概率，降低种群的生存能力和适应性。造成甘肃马鹿近交现象的原因可能与栖息地破碎化有关，人类活动导致马鹿栖息地被分割成小块，使得种群间的基因交流受到限制，从而增加了近亲繁殖的机会。此外，种群数量的减少也可能使得马鹿在选择配偶时的范围缩小，进一步加剧了近交现象。在甘肃马鹿的保护过程中，需要采取措施促进种群间的基因交流，扩大种群数量，以降低近交程度，提高种群的遗传健康水平。3.2遗传结构分析3.2.1主成分分析（PCA）利用主成分分析（PCA）方法对甘肃马鹿种群的遗传结构进行分析。将甘肃马鹿在10个微卫星位点上的基因频率数据导入到R语言软件中，运用相关的R包（如adegenet）进行PCA分析。通过PCA分析，将高维的遗传数据降维到低维空间，以直观地展示甘肃马鹿个体之间的遗传关系。在PCA分析结果的散点图中，横坐标和纵坐标分别代表第一主成分（PC1）和第二主成分（PC2），每个点表示一个甘肃马鹿个体。结果显示，第一主成分和第二主成分累计贡献率达到[X]%，能够较好地反映甘肃马鹿种群的遗传变异情况。从散点图的分布来看，甘肃马鹿个体并没有呈现出明显的聚类分布，而是相对较为分散地分布在图中。这表明甘肃马鹿种群内部的遗传结构相对较为复杂，个体之间的遗传差异较小，不存在明显的遗传亚结构。例如，在某些研究中，当种群存在明显的地理隔离或遗传分化时，个体在PCA散点图中会呈现出明显的聚类现象，不同的聚类代表不同的遗传亚群。而甘肃马鹿的这种分散分布情况，可能是由于其栖息地相对连续，种群间的基因交流较为频繁，没有形成明显的地理隔离，从而使得遗传差异在个体间均匀分布。3.2.2结构分析（STRUCTURE）借助STRUCTURE软件，对甘肃马鹿的遗传结构进行深入推断。将甘肃马鹿样本的微卫星数据整理成符合STRUCTURE软件输入格式的文件，设置相关参数进行分析。在参数设置中，将聚类数K的范围设定为1-10，运行多次独立模拟，每次模拟的迭代次数设置为100000次，燃烧期（burn-inperiod）设置为50000次。通过多次模拟，确保结果的可靠性和稳定性。运行STRUCTURE软件后，利用StructureHarvester软件对结果进行处理，计算不同K值下的似然值（LnP(D)）和ΔK值。似然值反映了数据与假设聚类数K的拟合程度，ΔK值则用于确定最佳的聚类数。结果表明，当K=[X]时，ΔK值达到峰值，似然值也相对较高，因此确定甘肃马鹿的最佳聚类数为[X]。这意味着甘肃马鹿种群可以划分为[X]个遗传簇。从STRUCTURE软件输出的结果图中可以直观地看到，不同颜色代表不同的遗传簇，每个个体被分配到不同遗传簇的比例用柱状图表示。在甘肃马鹿种群中，部分个体的遗传组成较为单一，主要来自于某一个遗传簇；而另一些个体则表现出混合的遗传组成，由多个遗传簇的基因组成。这种遗传组成的差异反映了甘肃马鹿种群内部存在一定程度的遗传多样性和遗传分化。例如，在某些区域的甘肃马鹿个体中，可能由于长期的地理隔离或生态环境差异，与其他区域的马鹿基因交流较少，导致其遗传组成相对单一；而在一些基因交流频繁的区域，马鹿个体则呈现出混合的遗传组成。通过STRUCTURE分析，能够更深入地了解甘肃马鹿种群的遗传结构，为进一步研究其遗传多样性和保护策略提供重要依据。3.3遗传多样性结果讨论3.3.1遗传多样性水平评估甘肃马鹿在10个微卫星位点上展现出一定程度的遗传多样性。从等位基因数来看，共检测到[X]个等位基因，平均等位基因数为[X]个，这一数值相较于一些遗传多样性较为丰富的马鹿种群，如塔里木马鹿平均等位基因数达到[X]个，处于相对较低水平；但与部分受到严重威胁、遗传多样性急剧下降的马鹿种群相比，仍保持着一定的遗传变异基础。较多的等位基因意味着种群拥有更丰富的遗传信息，这些遗传信息是种群在长期进化过程中积累下来的宝贵资源，能够为种群适应环境变化提供更多的可能性。当面对气候变化、食物资源波动或新的疾病威胁时，丰富的等位基因可能使种群中存在能够适应这些变化的个体，从而保证种群的生存和繁衍。有效等位基因数能够更准确地反映基因的有效变异程度，甘肃马鹿群体的平均有效等位基因数为[X]个，低于等位基因数。这表明在甘肃马鹿种群中，部分等位基因的频率较低，它们虽然增加了等位基因的数量，但对遗传多样性的实际贡献相对较小。例如，某些稀有等位基因可能只在极少数个体中存在，在种群的遗传结构和适应性中所起的作用有限。然而，这些稀有等位基因仍然具有潜在的价值，它们可能携带了一些特殊的遗传信息，在未来环境发生剧烈变化时，这些稀有等位基因可能会成为种群适应新环境的关键。基因多样性反映了种群中基因的杂合程度和变异水平，甘肃马鹿的平均基因多样性为[X]，显示出其具有较为丰富的遗传变异。较高的基因多样性使得甘肃马鹿种群在面对复杂多变的环境时，具有更强的适应能力。例如，当栖息地的生态环境发生改变时，具有丰富基因多样性的种群更有可能拥有适应新环境的基因组合，从而提高种群的生存概率。与其他地区马鹿种群相比，甘肃马鹿的基因多样性处于[较高/较低/相当]水平。这种差异可能与它们的进化历史、栖息地环境以及人类活动的影响程度密切相关。杂合度分析结果显示，甘肃马鹿的观察杂合度平均为[X]，期望杂合度平均为[X]，观察杂合度低于期望杂合度。这一现象表明甘肃马鹿种群中存在一定程度的近交现象或杂合子缺失。近交可能导致有害隐性基因的纯合表达，增加种群中遗传疾病的发生概率，降低种群的生存能力和适应性。例如，在一些近亲繁殖较为严重的动物种群中，出现了生长发育迟缓、免疫力下降等问题。造成甘肃马鹿近交现象的原因可能是多方面的，其中栖息地破碎化是一个重要因素。人类活动导致马鹿栖息地被分割成小块，使得种群间的基因交流受到限制，马鹿在选择配偶时的范围缩小，从而增加了近亲繁殖的机会。此外，种群数量的减少也进一步加剧了近交现象。3.3.2影响遗传多样性的因素栖息地破碎化是影响甘肃马鹿遗传多样性的重要因素之一。随着人类活动的不断扩张，如大规模的森林砍伐、道路建设和农业开垦，甘肃马鹿的栖息地被分割成众多小块，形成了一个个孤立的“岛屿”。这种破碎化的栖息地使得马鹿种群之间的交流变得困难，基因流动受到严重阻碍。不同种群的马鹿难以进行有效的基因交换，导致遗传漂变作用增强，一些稀有等位基因可能在小种群中逐渐丢失，从而降低了整个种群的遗传多样性。例如，在祁连山林区，由于道路的修建，将原本连续的马鹿栖息地一分为二，使得两侧的马鹿种群难以相互接触，基因交流几乎中断。长期以往，这两个种群在遗传上的差异逐渐增大，遗传多样性也随之降低。人类活动的干扰对甘肃马鹿的遗传多样性产生了负面影响。非法捕猎行为使得甘肃马鹿的种群数量急剧减少，导致遗传瓶颈效应的出现。当种群数量大幅下降时，许多等位基因可能会随着个体的消失而丢失，种群的遗传多样性也会随之降低。同时，人类活动还可能改变马鹿的栖息地环境，影响其食物资源和生存空间。例如，过度放牧导致草原退化，马鹿的食物来源减少；工业污染可能破坏马鹿的栖息地生态平衡，影响其繁殖和生存。这些因素都间接影响了甘肃马鹿的遗传多样性。种群数量的变化与遗传多样性密切相关。当种群数量充足时，个体之间的基因交流频繁，能够有效地维持遗传多样性。然而，由于栖息地破坏和人类捕猎等原因，甘肃马鹿的种群数量不断减少。小种群更容易受到遗传漂变的影响，基因频率的随机波动可能导致某些等位基因的固定或丢失，从而降低遗传多样性。此外，小种群中近亲繁殖的概率增加，进一步加剧了遗传多样性的下降。例如，在一些甘肃马鹿种群数量稀少的地区，已经观察到了遗传多样性明显降低的现象，表现为等位基因数减少、杂合度降低等。自然选择和遗传漂变在甘肃马鹿的遗传多样性形成和维持中也发挥着重要作用。在自然选择的作用下，适应环境的基因得以保留和传播，不适应环境的基因逐渐被淘汰。例如，在祁连山林区，气候寒冷，冬季食物资源匮乏，那些具有更强耐寒能力和适应冬季食物的基因在马鹿种群中逐渐占据优势。而遗传漂变则是由于种群数量有限，基因频率在世代传递过程中发生随机变化。在小种群中，遗传漂变的作用更为明显，可能导致一些对种群生存和繁衍没有明显优势或劣势的基因随机固定或丢失，从而影响遗传多样性。四、甘肃马鹿与青海马鹿亲缘关系分析4.1遗传距离计算运用Popgene3.2软件，基于Nei's遗传距离模型计算甘肃马鹿与青海马鹿之间的遗传距离。Nei's遗传距离能够准确衡量两个种群在基因频率上的差异程度，通过比较等位基因频率在不同种群间的分布情况，反映种群间的遗传分化水平。在计算过程中，软件将甘肃马鹿和青海马鹿在10个微卫星位点上的等位基因频率数据进行分析处理。结果显示，甘肃马鹿与青海马鹿之间的遗传距离为[X]。这一数值表明，二者在遗传物质上存在一定程度的差异，但差异并不十分显著。与其他相关研究中不同马鹿种群间的遗传距离进行对比。例如，在对塔里木马鹿与天山马鹿的研究中，它们之间的遗传距离为[X]，明显大于甘肃马鹿与青海马鹿之间的遗传距离。这说明甘肃马鹿与青海马鹿在亲缘关系上相对更近，在演化历程中可能有着更为紧密的联系，基因交流更为频繁。而塔里木马鹿与天山马鹿由于地理隔离等因素，在长期的进化过程中积累了更多的遗传差异，导致遗传距离较大。遗传距离的大小受到多种因素的影响。地理分布是一个重要因素，甘肃马鹿与青海马鹿的栖息地存在部分重叠，特别是在祁连山地区，这使得它们在历史上有更多机会进行基因交流，从而减小了遗传距离。相比之下，地理隔离明显的马鹿种群，由于缺乏基因交流，遗传差异逐渐积累，遗传距离增大。此外，生态环境的相似性也会对遗传距离产生影响。如果两个种群所处的生态环境相似，它们面临的自然选择压力相近，基因频率的变化趋势可能较为一致，进而遗传距离相对较小。甘肃马鹿和青海马鹿在一定程度上共享相似的生态环境，这也可能是导致它们遗传距离相对较小的原因之一。4.2系统发育树构建4.2.1邻接法（NJ）运用MEGA7.0软件，基于邻接法构建甘肃马鹿与青海马鹿的系统发育树。邻接法是一种基于距离矩阵的聚类算法，其核心思想是通过逐步合并距离最近的分类单元，构建系统发育树。在构建过程中，首先利用Popgene3.2软件计算得到的Nei's遗传距离矩阵作为输入数据。该遗传距离矩阵包含了甘肃马鹿和青海马鹿各个个体之间的遗传距离信息，是邻接法构建系统发育树的基础。在MEGA7.0软件中，进行如下参数设置：选择“Neighbor-Joining”作为建树方法，bootstrap检验的重复次数设置为1000次。bootstrap检验是一种用于评估系统发育树分支可靠性的方法，通过多次重复抽样和建树，统计每个分支在重复构建的树中出现的频率，以此来判断分支的可信度。当bootstrap值较高（如大于70%）时，表明该分支在多次重复构建中较为稳定，可信度较高。构建完成的系统发育树结果显示，甘肃马鹿和青海马鹿的个体并没有完全按照种群各自聚类。部分甘肃马鹿个体与青海马鹿个体在系统发育树上相互交错分布，这直观地表明甘肃马鹿与青海马鹿之间存在着较为紧密的亲缘关系。它们在演化过程中可能存在频繁的基因交流，导致部分个体的遗传特征相似，无法依据种群进行明显的区分。例如，在系统发育树的某一分支上，同时出现了来自甘肃马鹿种群和青海马鹿种群的个体，且该分支的bootstrap值达到了[X]%，说明这一分支的可信度较高，进一步证实了两者之间存在紧密的遗传联系。4.2.2最大似然法（ML）采用最大似然法在MEGA7.0软件中构建系统发育树，以与邻接法的结果进行比较。最大似然法是基于概率模型的建树方法，它假设在给定的进化模型下，某棵系统发育树能够最合理地解释所观察到的遗传数据，那么这棵树就是最有可能的系统发育树。在运用最大似然法时，首先需要选择合适的核苷酸替代模型。通过ModelTest软件进行模型选择，根据赤池信息准则（AIC）等指标，确定最适合本研究数据的核苷酸替代模型为[具体模型名称]。该模型能够准确地描述DNA序列在进化过程中的碱基替换规律，为最大似然法构建系统发育树提供可靠的基础。在MEGA7.0软件中，选择“MaximumLikelihood”作为建树方法，并根据ModelTest软件的结果设置核苷酸替代模型参数。同样进行1000次bootstrap检验，以评估分支的可靠性。最大似然法构建的系统发育树结果与邻接法具有一定的相似性。在最大似然法构建的系统发育树中，甘肃马鹿和青海马鹿的个体同样存在相互交错的分布情况，这进一步验证了两者之间存在密切的亲缘关系。然而，两种方法构建的系统发育树在一些细节上也存在差异。例如，某些分支的长度和拓扑结构在两种方法构建的树中略有不同。这些差异可能是由于两种建树方法的原理和假设不同导致的。邻接法基于遗传距离进行聚类，而最大似然法基于概率模型进行推断。尽管存在这些差异，但两种方法都一致表明甘肃马鹿与青海马鹿之间存在紧密的亲缘关系，从不同角度为研究两者的演化关系提供了有力的证据。4.3亲缘关系结果讨论4.3.1亲缘关系程度分析甘肃马鹿与青海马鹿之间的遗传距离为[X]，此数值表明二者在遗传物质上存在一定程度的差异，但差异并不十分显著。从遗传距离的角度来看，相较于其他一些马鹿种群间较大的遗传距离，如塔里木马鹿与天山马鹿之间的遗传距离达到[X]，甘肃马鹿与青海马鹿的遗传距离相对较小，这直观地反映出二者在亲缘关系上较为紧密。在系统发育树的构建结果中，无论是邻接法还是最大似然法，都显示甘肃马鹿和青海马鹿的个体存在相互交错分布的情况。部分甘肃马鹿个体与青海马鹿个体在系统发育树上无法依据种群进行明显区分，这进一步证实了它们之间存在紧密的亲缘关系。这种亲缘关系的紧密性暗示着两个种群在演化历程中可能有着频繁的基因交流。在漫长的历史进程中，由于它们的栖息地存在部分重叠，特别是在祁连山地区，为两个种群的马鹿提供了相遇和交配的机会，使得基因能够在两个种群间流动，从而减小了遗传差异，拉近了亲缘关系。4.3.2分化原因探讨尽管甘肃马鹿与青海马鹿亲缘关系较近，但仍然存在一定程度的遗传分化，这可能是由多种因素共同作用导致的。地理隔离是一个重要因素，虽然二者栖息地有重叠部分，但在一些边缘区域或特定的地理条件下，仍然可能存在一定程度的隔离。例如，青海马鹿分布的柴达木盆地周边区域，与甘肃马鹿主要栖息的祁连山林区在地理环境上存在较大差异，盆地周边的荒漠和半荒漠景观形成了相对隔离的生态环境，使得青海马鹿在这一区域逐渐形成了独特的遗传特征。随着时间的推移，不同地理区域的马鹿种群在基因交流受限的情况下，各自积累了不同的遗传变异，从而导致遗传分化。生态环境的差异也对二者的分化起到了推动作用。甘肃马鹿栖息地的祁连山林区气候相对湿润，植被类型丰富，以高山森林和山地草原为主；而青海马鹿栖息地的部分区域，如柴达木盆地周边，气候干旱，植被以耐旱的荒漠植物为主。不同的生态环境对马鹿的生存和繁殖产生了不同的选择压力。在湿润的祁连山林区，马鹿可能更适应以草本植物和树叶为食，对水源的依赖程度较高；而在干旱的柴达木盆地周边，马鹿则需要适应耐旱植物作为食物，并且具备更强的水分利用能力。在长期的自然选择过程中，适应不同生态环境的基因在两个种群中逐渐固定下来，导致了遗传分化。此外，人类活动也可能对甘肃马鹿和青海马鹿的遗传分化产生了影响。随着人类活动范围的不断扩大，如道路建设、农业开垦和工业发展，破坏了马鹿的栖息地，导致栖息地破碎化。这使得原本连续的马鹿种群被分割成多个小种群，限制了种群间的基因交流。对于甘肃马鹿和青海马鹿来说，栖息地破碎化可能加剧了它们之间的遗传分化。小种群更容易受到遗传漂变的影响，基因频率的随机波动可能导致某些等位基因在小种群中逐渐固定或丢失，从而增加了种群间的遗传差异。同时，人类的捕猎活动也可能改变了马鹿的种群结构和遗传组成，进一步影响了它们的遗传分化进程。五、保护建议5.1基于遗传多样性的保护策略根据本研究对甘肃马鹿遗传多样性的分析结果，甘肃马鹿在10个微卫星位点上呈现出一定程度的遗传多样性，但也存在近交现象和遗传多样性下降的潜在风险。为了有效保护甘肃马鹿，应采取以下基于遗传多样性的保护策略。鉴于甘肃马鹿部分种群存在近交现象，导致有害隐性基因纯合表达的风险增加，进而降低种群的生存能力和适应性，有必要建立遗传管理档案。详细记录每只马鹿的谱系信息，包括其父母、祖父母等亲缘关系，以及在各个微卫星位点上的基因信息。通过遗传管理档案，能够准确了解马鹿个体之间的亲缘关系，避免近亲繁殖。在进行人工繁育或种群迁移时，根据遗传管理档案，选择亲缘关系较远的个体进行配对或迁移，促进基因交流，增加遗传多样性。例如，在某马鹿繁育中心，通过建立遗传管理档案，对马鹿的亲缘关系进行严格把控，成功避免了近亲繁殖，使得繁育出的马鹿种群遗传健康状况得到显著改善。为了促进甘肃马鹿种群间的基因交流，可采取种群迁移和再引入等措施。对于一些由于栖息地破碎化而导致基因交流受阻的小种群，可以选择合适的时机，将部分个体迁移到其他种群中，增加种群间的基因流动。在迁移过程中，要充分考虑迁移个体的健康状况、年龄结构以及目标种群的容纳能力等因素。例如，在经过充分的科学论证和前期准备后，将甘肃马鹿某小种群中的部分年轻、健康个体迁移到与之亲缘关系较近但遗传差异相对较大的另一种群中。迁移后，通过定期监测种群的遗传多样性指标，发现新种群的遗传多样性得到了有效提升，种群的适应性和生存能力也有所增强。同时，对于已经灭绝或数量极少的甘肃马鹿种群分布区域，可以进行再引入工作。选择遗传多样性丰富、健康状况良好的个体，经过驯化和适应期后，放归到适宜的栖息地中，重建种群。在再引入过程中，要加强对放归个体的跟踪监测，确保其能够适应新环境并成功繁衍后代。5.2基于亲缘关系的保护措施由于甘肃马鹿与青海马鹿存在紧密的亲缘关系，且在演化历程中可能有着频繁的基因交流，在制定保护策略时，应将二者视为一个整体进行综合保护。建立跨区域的联合保护机制，加强甘肃和青海两地在马鹿保护方面的合作与交流。通过共享科研数据、技术经验和管理模式，共同制定科学合理的保护规划，实现资源的优化配置和高效利用。例如，在祁连山地区，甘肃和青海可以共同划定马鹿的核心保护区，加强对该区域的生态保护和监测，确保马鹿栖息地的完整性和生态功能的稳定性。在人工繁育过程中，充分考虑甘肃马鹿和青海马鹿的亲缘关系，避免近亲繁殖。制定科学的繁育计划，合理选择配对个体，促进基因交流，以维持和提高种群的遗传多样性。可以建立马鹿遗传资源库，保存马鹿的精液、胚胎等遗传材料，为人工繁育提供丰富的遗传资源。同时，利用现代生物技术，如胚胎移植、人工授精等，提高繁育效率和质量。例如，通过人工授精技术，将具有优良遗传特性的甘肃马鹿或青海马鹿的精液引入到其他种群中，增加种群的遗传多样性。加强对马鹿栖息地的生态修复工作，针对由于人类活动导致的栖息地破碎化、退化等问题，采取植树造林、草原恢复、湿地保护等措施，改善马鹿的生存环境。在生态修复过程中，要充分考虑甘肃马鹿和青海马鹿的生态需求，营造适宜它们生存和繁衍的生态环境。例如，在祁连山林区和柴达木盆地周边，通过种植马鹿喜爱的食物植物，增加食物资源的丰富度；恢复湿地生态系统，为马鹿提供水源和栖息场所。同时，减少人类活动对马鹿栖息地的干扰，限制在马鹿栖息地内的开发建设活动，确保马鹿有足够的生存空间。5.3保护工作的展望未来，甘肃马鹿和青海马鹿的保护工作任重道远，需要多方面的持续努力。在科研领域，应进一步深化对马鹿遗传多样性和亲缘关系的研究。随着分子生物学技术的不断发展，可运用全基因组测序等更先进的技术手段，全面深入地揭示马鹿的遗传信息。这将有助于发现更多与马鹿生存、繁殖和适应环境相关的基因，为保护策略的制定提供更精准的遗传依据。例如，通过全基因组测序，能够准确识别出在马鹿适应高海拔环境或应对疾病过程中起关键作用的基因，从而有针对性地开展保护工作。在保护行动方面，要加大对马鹿栖息地的保护力度。建立更多的自然保护区和生态廊道，确保马鹿栖息地的连通性，减少栖息地破碎化对马鹿种群的影响。同时，加强对保护区的管理和监测，运用现代信息技术，如卫星遥感、无人机监测等，实时掌握马鹿的种群动态和栖息地变化情况。在内蒙古阿鲁科尔沁旗境内的高格斯台罕乌拉国家级自然保护区，通过实行全区域围封管理，设置多个管护站和瞭望站，24小时守护巡山护林，有效保护了野生马鹿的繁衍生息。未来可借鉴类似的成功经验，加强对甘肃马鹿和青海马鹿栖息地的保护管理。公众教育也是马鹿保护工作的重要环节。通过开展科普宣传活动，提高公众对马鹿保护重要性的认识，增强公众的保护意识。鼓励公众参与到马鹿保护行动中来，形成全社会共同保护的良好氛围。例如，组织志愿者参与马鹿栖息地的巡查和保护工作，举办马鹿保护主题的科普讲座和展览等。此外，还可以利用社交媒体、网络平台等传播渠道，广泛宣传马鹿的生态价值和保护现状，吸引更多人关注马鹿保护事业。六、结论与展望6.1研究结论总结本研究运用微卫星DNA技术，对甘肃马鹿的遗传多样性进行了深入分析，并探究了其与青海马鹿的亲缘关系。在遗传多样性分析方面，10个微卫星位点共检测到[X]个等位基因，平均等位基因数为[X]个，显示出甘肃马鹿具有一定的遗传变异，但相较于部分遗传多样性丰富的马鹿种群，等位基因数处于相对较低水平。平均有效等位基因数为[X]个，低于等位基因数，表明部分等位基因频率较低，对遗传多样性实际贡献有限。平均基因多样性为[X]，说明甘肃马鹿种群遗传变异较为丰富，具备一定的适应环境变化的潜力。然而，观察杂合度平均为[X]，期望杂合度平均为[X]，观察杂合度低于期望杂合度，揭示出甘肃马鹿种群存在近交现象或杂合子缺失，这可能对种群的生存和繁衍产生不利影响。通过主成分分析（PCA），发现甘肃马鹿个体在散点图上相对分散，种群内部遗传结构复杂，个体间遗传差异较小，无明显遗传亚结构，这可能与栖息地相对连续、基因交流频繁有关。利用STRUCTURE软件分析，确定甘肃马鹿最佳聚类数为[X]，种群可划分为[X]个遗传簇，部分个体遗传组成单一，部分个体呈现混合遗传组成，反映出种群内部存在一定程度的遗传多样性和遗传分化。在亲缘关系分析中，甘肃马鹿与青海马鹿的遗传距离为[X]，表明二者遗传物质存在一定差异，但并不显著，亲缘关系较为紧密。基于邻接法和最大似然法构建的系统发育树均显示，甘肃马鹿和青海马鹿个体相互交错分布，进一步证实了它们之间存在频繁的基因交流和紧密的亲缘关系。然而，二者仍存在一定遗传分化，地理隔离、生态环境差异以及人类活动干扰是导致分化的主要因素。青海马鹿分布的柴达木盆地周边区域与甘肃马鹿主要栖息的祁连山林区地理环境不同，