大绒鼠Eothenomys miletus隶属于仓鼠科Cricetidae绒鼠属，在云南省分布广泛(罗泽珣，2000)，是云南野鼠鼠疫自然疫源地的主要储存和流行宿主动物之一。云南省特殊的地理环境和气候条件，可能导致了大绒鼠形态特征不稳定、同种异型和近缘种形态特征相似等问题：王应祥(2003)将大绒鼠分为2个地理亚种：指名亚种E. m. miletus和怒江亚种E. m. confinii；整个绒鼠属的系统关系及分类地位尚存争议(刘少英，刘洋，2005)。因此，云南不同地理位置和气候环境中分布的大绒鼠的遗传多样性与遗传分化值得研究和探讨。

动物线粒体DNA(mtDNA)分子小而稳定、母系遗传、进化速度快等特征(Birky et al., 1989)，已成为遗传多样性研究的一个强有力的工具，它包含的细胞色素b(cyt b)基因是研究物种地理变异问题最可信的线粒体DNA标记之一(李楠楠等，2012)。一些学者对大绒鼠的遗传结构进行了初步研究(沐远等，2015；Mu et al., 2019)，但尚未对其他区域分布的大绒鼠的遗传信息进行比较分析，本文对云南7个地区的大绒鼠cyt b基因进行检测分析，为大绒鼠乃至绒鼠属的系统研究提供分子数据。

本研究从玉龙县、剑川县、祥云县、大理市、巧家县、禄劝县、腾冲市的7个采样点共采集到大绒鼠140只(表 1)。利用鼠夹从采样点的农田及灌木丛中捕获大绒鼠，捕获个体按形态学特征，并剥制头骨标本进行鉴定，取肌肉保存于95%乙醇中，-20 ℃保存。

用肌肉标本按照QIAGEN试剂盒产品说明书步骤提取基因组DNA。cyt b的PCR扩增使用哺乳动物通用引物L14724：5'-CGAAGCTTGATATGAAAAACCATCGTTG-3'，H15915R：5'-GGAATTCATCTCTCCGGTTTACAAGAC-3'；扩增体系为25 μL：模板2 μL，10 μmol·L^-1引物各0.5 μL，Mix 12.5 μL，ddH₂O 9.5 μL。扩增条件：95 ℃ 5 min；95 ℃ 45 s，55 ℃ 1 min，72 ℃ 1 min，35个循环；72 ℃ 10 min。实验过程以水作为空白对照。扩增产物经琼脂糖凝胶电泳检测后送上海生工生物工程技术服务有限公司测序。

利用Clustal X(Thompson et al., 1997)对测序获得的有效样本序列进行比对。利用MEGA 5(Tamura et al., 2011)基于Kimura双参数距离模型计算碱基组成、转换与颠换值；以齐氏姬鼠Apodemus chevrieri(GenBank登录号：AB096821)为外群，基于类平均法(UPGMA)构建分子系统发育树，各分支的自举检验值由1 000次重复检验所得。利用Dnasp 5.0(Rozas et al., 2003)计算单倍型数目、变异位点数、简约信息位点数、单倍型多样性、核苷酸多样性和基因流(N_m)。其中，N_m < 1表明群体间可能由于遗传漂变而发生分化；N_m>1表明群体间的基因流水平较高，群体间遗传分化较小。利用PopART(Leigh & Bryant，2015)基于中连法(Median joining method)构建大绒鼠单倍型间的网络关系图。运用Arlequin 3.1(Excoffier et al., 1992)计算遗传分化系数(F_st)，进行分子变异分析(AMOVA)。通过IBD 1.53(Bohonak，2002)检测地理距离与遗传距离的相关性，用F_st值表示遗传距离，地理距离采用采样点间的地球表面直线距离。

扩增结果经琼脂糖凝胶电泳检测未发现非特异性条带，共获得140只大绒鼠的扩增条带，部分样本扩增结果见图 1。

图 1 大绒鼠cyt b基因部分序列扩增结果 Fig. 1 Amplification of partial cyt b gene of Eothenomys miletus 1~4. 样本, M. Marker 1-4. sample, M. Marker

图选项

共获得140只大绒鼠样本cyt b基因部分序列(GenBank登录号：MW567503~MW567642)，在1 027 bp的序列中，共有核苷酸变异位点128个，包含单变异位点109个和简约信息位点19个，转换和颠换(TS/TV)比为2.0，A、T、G、C平均碱基组成分别为：29%、27%、30.4%、13.6%，A+T含量显著高于G+C含量。109个变异位点共定义51种单倍型，分别命名为Hap1~Hap51。

7个采样点大绒鼠的单倍型多样性为0.711~ 0.907，核苷酸多样性为0.002 25~0.002 98，整体的单倍型多样性和核苷酸多样性分别为0.977± 0.004和0.018 45±0.001 78。所有种群均为高单倍型多样性和低核苷酸多样性的模式。Hap20为大理、祥云、剑川种群共享，Hap21、Hap22、Hap24为大理、祥云种群共享，其余47个单倍型为各个种群特有(表 2)。玉龙和剑川种群的单倍型多样性相对较高，巧家种群的核苷酸多样性最高(表 3)。

种群间遗传差异分析结果显示，7个地理种群的F_st值为0.179 68~0.949 45，遗传差异均有统计学意义(P＜0.001)(表 4)。分子变异分析(AMOVA)结果表明，群体间的遗传变异为84.29%，群体内的为15.71%。种群间的基因流为0.026 62~2.282 72，除玉龙、祥云、剑川、大理4个地理距离较近种群之间的N_m值大于1之外，其他均小于1，且值均较低(表 5)。

采用UPGMA法构建的系统进化树显示，7个采样点的大绒鼠种群分化为4个地域分支(图 2)，第一分支包括玉龙、祥云、剑川、大理的单倍型，第二分支为禄劝的单倍型，第三分支为巧家单倍型，第四分支为腾冲单倍型。单倍型进化网络关系图显示出与进化树相同的结果(图 3)，大绒鼠已经出现了地理种群分化。

图 2 基于云南7个采样点的大绒鼠cyt b基因序列构建的单倍型UPGMA树 Fig. 2 UPGMA tree of haplotypes based on cyt b gene of Eothenomys miletus from 7 localities of Yunnan Province

图选项

图 3 基于cyt b基因序列建立的云南7个采样点的大绒鼠单倍型进化网络关系图 Fig. 3 Network relationships of haplotypes based on cyt b gene of Eothenomys miletus from 7 localities of Yunnan Province 每个实心圆代表 1个单倍型, 大小与单倍型出现的频率呈正比; 空心圆代表缺失单倍型; 每条线代表 1次突变 Each solid circlerepresents a haplotype, and its size is proportional to haplotype frequency; an open circle represents a missing haplotype; each line represents a mutation

图选项

通过IBD1.53对大绒鼠7个地理种群间的地理距离与遗传距离进行检验，结果显示，地理距离与遗传距离之间呈显著的正相关(r=0.699 7，P < 0.05，R²=0.49)。

遗传多样性的高低与物种的适应能力、生存能力密切相关(季维智等，1999；Frankham et al., 2010)，本研究中大绒鼠单倍型多样性为0.977±0.004，核苷酸多样性为0.018 45±0.001 78，多样性水平与其常同域分布的齐氏姬鼠相当(朱万龙等，2013)，玉龙和剑川大绒鼠种群单倍型多样性高于其他种群，大绒鼠在这些地区小兽中的构成比也较高，说明大绒鼠在这些地区的适应能力和生存能力较强。大绒鼠种群总体表现出高单倍型多样性、低核苷酸多样性的特点，Yuan等(2012)认为高单倍型多样性、低核苷酸多样性主要与单个碱基的变异、地理种群处于扩张有关，大绒鼠种群总体表现出高单倍型多样性、低核苷酸多样性的特点可能也与这2种原因有关。

基因分化系数和基因流是衡量遗传分化的重要指标(Nei，1987)，F_st值为0.15~0.25时说明遗传分化较大(Wright，1965)，基于线粒体cyt b序列的AMOVA分析结果显示，大绒鼠各地理种群间遗传分化指数为0.179 68~0.949 45，表明云南地区大绒鼠各地理种群间存在明显的遗传分化。大绒鼠7个地理种群间基因流为0.026 62~2.282 72，玉龙、祥云、剑川、大理4个地理距离较近种群组合之间的N_m值大于1之外，其他种群组合间均小于1，且值均较低，说明多数种群间基因流水平贫乏。在单倍型系统进化树中，大绒鼠种群分化为4个地域分支，玉龙、祥云、剑川、大理为一分支，其他种群分别独立为一支，在实际的地理位置中，玉龙、剑川、大理、祥云的地理距离较近，而腾冲、巧家、禄劝与其他几个地区距离较远。单倍型进化网络关系图显示出与进化树相同的结果，综合研究结果表明大绒鼠已经出现了明显的种群分化。

遗传分化的程度会随着地理距离的增加而增加(Floyd et al., 2005)，IBD检测结果显示，F_st与地理距离呈显著正相关，大绒鼠种群间遗传分化中49.0%的变异可以用地理隔离解释，这表明种群间地理隔离对遗传分化起重要作用，这可能与大绒鼠的迁移扩散能力较弱有关。云南气候多样，山川河流纵横，地形地势复杂，独特的生境可能造成大绒鼠种群间的地理隔离，是否因地理隔离形成独立的亚种还有待于进一步研究。

综上所述，云南地区大绒鼠高单倍型多样性、低核苷酸多样性，以及种群间高遗传分化的遗传结构特点，可能与其生物学特征和扩散能力等因素有关，同时地理距离隔离对其遗传分化起了重要作用。然而，由于本研究仅用了线粒体cyt b基因作为分子标记，可能存在局限性，欲全面了解大绒鼠的遗传多样性、遗传分化程度，还需要扩大大绒鼠采样范围，加入微卫星、核基因等标记的分析，进一步研究。