1.内蒙古自治区地方病防治研究中心 呼和浩特 010031
【摘 要】[目的]为了探究应用DNA 条形码技术对鼠疫宿主动物进行分子生物学鉴定的可行性。[方法]本研究检测了内蒙古蒙古旱獭疫源地6 种啮齿动物的75 条线粒体COⅠ基因657bp 的片段序列，用MEGA5.0 软件计算GC 含量，经过双参数法计算遗传距离，采用邻接法（NJ）构建分析系统发育树。[结果]所有序列富含AT，所有物种G 含量及第3 位的G 含量均较少，各鼠种碱基含量随着物种以及密码子位置的不同而不同。种内平均遗传距离为0.1%，不同鼠种间平均遗传距离为22.63%，NJ 进化树能清楚的区分6 个不同的类群。[结论] 内蒙古蒙古旱獭疫源地的鼠疫宿主动物通过COⅠ基因应用DNA 条码技术可以进行鉴别，同时认为近年蒙古旱獭疫源地达乌尔黄鼠为优势种。
【关键词】DNA 条形码；COⅠ基因；鼠疫
【中图分类号】Q342+.2【文献标识码】A【文章编号】2096-0867（2015）-10-006-02
蒙古旱獭疫源地处于呼伦贝尔高原北部草原地带，东经115°50′～121°18′，北纬48°20′～49°50′之间，其范围涉及满洲里市和陈巴尔虎旗、新巴尔虎左旗、新巴尔虎右旗、鄂温克旗及牙克石市的一些地区。根据70 年代后期动物病调查队的调查和80 年代以来系统的鼠疫监测结果，在该疫源地共发现各种啮齿动物21 种[1]。但近年来在蒙古旱獭疫源地共获鼠 14 种[2]，其中常见种有达乌尔黄鼠、长爪沙鼠、蒙古旱獭、黑线毛足鼠、五趾跳鼠等，蒙古旱獭为优势种，在该疫源地2000-2011 年见共有 2 种啮齿动物（达乌尔黄鼠和五趾跳鼠）自然感染鼠疫，为了准确了解疫源地内啮齿动物的种类和数量，现利用DNA 条形码方法结合形态鉴定对蒙古旱獭疫源地内的鼠疫宿主动物进行鉴定研究。
1 材料和方法
采集内蒙古自治区呼伦贝尔高原蒙古旱獭鼠疫自然疫源地内的啮齿动物，现场初步进行鼠种形态鉴定，同时进行性别、年龄鉴定和体尺测定后，带回实验室解剖取肝脏组织，-80℃冷藏保存备用，实验样品信息见表1。
啮齿动物肝脏组织基因组DNA 提取方法按照TaKaRa 试剂盒指导书（D305A）的操作步骤进行。扩增时所用引物及使用条件参考Robins 等[3]（2007）的鼠类通用引物和Ivanova 等[4]（2007）的鸡尾酒引物，并由上海生工生物技术有限公司合成。PCR 产物经1.0%琼脂糖凝胶电泳检验样品的完整性，检测合格的PCR 产物直接送上海生工双向测序。同一份样品进行3 次PCR 试验，测序时应用自带的鸡尾酒引物。测序成功的序列首先用Chromos 软件观察测序峰图质量，如果峰图质量很差，不能准确判断碱基，则重新进行扩增和测序。再用ClustalX 软件进行多重序列比对，并进行适当的手工调整。然后将测定的序列在NCBI 上运行BLAST 程序进行序列同源性比较，以确保所获得的序列是目标序列。最后运用Mega5 软件分析各基因序列的碱基组成和变异， 基于K2P（Kimura-2-parameter）模型计算各物种的遗传距离。采用邻接法构建全部COⅠ基因序列的NJ（Neighbor-Joining）系统树，对NJ 树进行内部分支检验与1000 次Bootstrap 检验分析，来确定各支系的置信度。
表1 实验样品信息.png)
2 结果与分析
本文共获得75 条COⅠ基因序列片段（每条序列长度657bp），对所有序列的阅读框架进行确认，避免核基因组中的假线粒体基因序列干扰。75 个样本属于2 目5 科6 属6 种。
2.1 COⅠ基因序列分析
2.1.1 变异位点分析
对6 种啮齿动物75 条序列COⅠ进行比对，发现有变异位点234 个。所有DNA 条码都是物种所独有，不存在物种之间的共享；对所涉及的6 个物种，每个物种都有各自的鉴别位点。
2.1.2 GC 含量分析
各种碱基含量随着物种和密码子位置的不同而不同。达乌尔鼠兔和五趾跳鼠核苷酸含量C>T>A>G，而其他4 个鼠种的核苷酸含量T>A>C>G，其差别主要来自于C 碱基含量的差别。这些序列大部分是富含AT，GC 含量差别较大。五趾跳鼠和达乌尔鼠兔的GC 含量较高，而达乌尔黄鼠的GC 含量则较低（见表2）。该差别来源于密码子第1、3 位的GC 含量的差别，主要是第3 位的差别所致（见表3），第2 位的GC 含量几乎无差别。所有物种第3 位的G 含量均较少。
表2 各鼠种碱基含量表（AV±SE）.png)
表3 各鼠种不同密码子位置GC%含量表（AV±SE）.png)
2.2 K2P 遗传距离分析
2.3 进化树分析
基于COⅠ序列，应用邻位相连法（NJ）构建的进化树图中可看出，各鼠种被分成各自不同的类群，能清楚的区分6个不同的品种。
3.讨论
Hebert 等[5-8]选取COⅠ基因对动物不同种群进行种类鉴定，发现同属各种COⅠ序列的平均差异程度是11.3%，而种内的COⅠ序列差异程度基本小于2%，绝大部分低于1%。Borisenko[9]、Nicolas V]对哺乳动物（包括鼠科）的研究中，种内平均遗传距离基本小于2%，均与Hebert 的研究结果一致。本研究中我们通过脊椎动物COⅠ通用鸡尾酒引物，对蒙古旱獭鼠疫自然疫源地内的啮齿动物6 种75 个样本的COⅠ基因进行序列分析，发现种内平均遗传距离为0.1%，种内最大遗传距离为0.6%，略低于别人的研究结果。
图1 基于K2P 模型构建的NJ 进化树
Hebert 研究发现93.8%的脊椎动物用CO.png)
Ⅰ同样的基因区段，同属种间的遗传距离为4 到32%，平均为9.6%。2009 年Eaton MJ 对中非和南美的脊椎动物研究发现同属种间的平均遗传距离为9.8%，Borisenko[9]和Nicolas V 研究的哺乳动物其种间平均遗传距离分别为10.1%和11.2%。本研究发现种间平均遗传距离的为22.63%，大于种内平均遗传距离37 倍之多，其种间遗传距离略大于其他人的研究结果，造成该现象的原因可能是由于本文中所选的6 个品种分属于5 个不同的属，种间距离其实大部分是代表不同科间的遗传距离。
本研究测定了蒙古旱獭疫源地常见6 种啮齿动物的COⅠ基因657bp 的多变区，基本包括该疫源地常见鼠疫宿主动物，每种动物都有自己单独的条形码，这为品种鉴定提供了有力的证据，不同位点上突变频率不同，每种啮齿动物有其特异位点，而且6 个品种被分为各自不同的类群，能清楚的区分6 个不同的品种，因此，应用DNA 条形码基本可以鉴定蒙古旱獭疫源地的啮齿动物，同时可为形态鉴别提供佐证数据。同时可以发现蒙古旱獭疫源地啮齿动物的分布情况，近年黄鼠为优势种。但本实验涉及的品种以及个别品种测序的样本数量还相对较小，更可靠的结果有待于扩大实验品种范围及数量进行分析。
参考文献：
[1].方喜业.中国鼠疫自然疫源地[M].人民卫生出版社，1990：106-108.
[2].范蒙光，尉瑞平，李建云，等.内蒙古鼠疫自然疫源地啮齿动物现况调查[J].医学动物防制.2013，29（1）：5-9.
[3].Judith H.Robins，Melanie H，Elizabeth Matisoo-smith，et al.Identifying Rattus species using mitochondrial DNA[J].MolEcol Notes，2007，7（5）：717–729.
[4].Ivanova NV，Zemlak TS，Hanner R，et al.Universal primer cocktailsfor fish DNA barcoding[J].Mol Ecol Notes，2007，7（4）：544–548.
[5].Hebert PDN，Cywinska A，Ball SL，et al.Biologicalidentifications through DNA barcodes[J].Proc Biol Sci，2003，270（1512）：313-321.
[6].Hebert PDN，Penton EH，Burns JM，et al.Ten species inone：DNA barcoding reveals cryptic species in the neotropicalskipper butterfly Astraptes fulgerator[J].Proc Natl Acad Sci U SA.2004，101（41）：14812-14817.
[7].Hebert PDN ， Ratnasingham S ， deWaard JR ， etal.Barcoding animal life cytochrome c oxidase subunit 1divergences among closely related species[J].Proc Biol Sci.2003Aug 7；270 Suppl 1：S96-99.
[8].Hebert PDN，Stoeckle MY，Zemlak TS，et al.Identificationof birds through DNA barcodes[J].PLoS Biol，2004，2（10）：e312.
[9].Borisenko AV，Lim BK，Ivanova NV，et al.DNA barcodingin surveys of small mammal communities ： a field study inSuriname[J].Mol Ecol Resour，2008，8（3）：471-479.
[10].Nicolas V，Schaeffer B，Missoup AD，et al.Assessment of threemitochondrial genes（16S，Cytb，CO1）for identifying species in thePraomyini tribe（Rodentia Muridae）[J].PLoS One.2012，7（5）：e36586.
[11].Eaton MJ，Meyers GL，Kolokotronis SO，et al，Barcodingbushmeat：molecular identification of Central African and South Americanharvested vertebrates[J].Conserv Genet，2010，11（4）：1389–1404.
基金项目： 2012 年卫生公益性行业科研专项（201202021），内蒙古自治区卫生和计划生育委员会2013 年医疗卫生科研计划项目（201302039）。
作者简介：常子丽（1980-），女，主管技师，博士，从事地方病防治工作。
论文作者:常子丽 张忠兵 刘 芳 李振林 赵成祥 杨晓娟 范蒙光
论文发表刊物:《系统医学》2015年第1卷第10期
论文发表时间:2016/2/23
标签：源地论文; 旱獭论文; 蒙古论文; 序列论文; 鼠疫论文; 动物论文; 含量论文; 《系统医学》2015年第1卷第10期论文;