基于线粒体 ATP6基因全序列分析中国绵羊痒螨(兔亚种)的遗传多样性
尹潇潇1, 李燕方1, 古江1, 廖艳2, 谢跃1, 杨光友1, 古小彬1,*
1.四川农业大学 动物医学院,四川 成都 611130
2.眉山市畜牧站,四川 眉山 620000
*通信作者,古小彬,E-mail: guxiaobin198225@126.com

作者简介:尹潇潇(1993—),女,四川成都人,硕士研究生,研究方向为动物寄生虫病学。E-mail: yinxiaoxiao1993@icloud.com

摘要

为进一步了解我国绵羊痒螨(兔亚种)种群的遗传变异情况和分类情况,采用PCR技术扩增绵羊痒螨(兔亚种)虫株的线粒体 ATP6基因全序列,并分析所得序列,旨在探讨中国华北、华东、华中、西北、西南地区绵羊痒螨(兔亚种)种群的遗传多样性和种群结构。本研究成功获得88条序列,长度均为672 bp,包含41个单倍型,表现出较高的遗传多样性( π=0.014 94)和单倍型多样性( Hd=0.925 81),且以种群内部的遗传变异为主(97.92%)。进一步分析发现,5个种群间遗传分化程度较弱( Fst=0.020 81),基因交流频繁( Nm=11.763 5),中性检验值(Tajima’s D=0.937 98,Fu’s Fs=0.522 06)为不显著的正值,结合错配分布曲线呈现多峰,表明我国绵羊痒螨(兔亚种)种群在进化过程中比较稳定,无种群扩张事件。从单倍型网络图和NJ树可知,中国绵羊痒螨(兔亚种)没有形成基于兔品种、温度带或地理分布的遗传结构。中国绵羊痒螨(兔亚种)种群遗传多样性高,但种群间无明显分化,未形成基于兔品种、温度带或地理分布的遗传结构。

关键词: 绵羊痒螨(兔亚种); 线粒体 ATP6基因; 遗传多样性
中图分类号:S858.26 文献标志码:A 文章编号:104-1524(2019)08-1231-08
Genetic diversity analysis of Psoroptes ovis var. cuniculi by full-length of mitochondrial ATP6 gene in China
YIN Xiaoxiao1, LI Yanfang1, GU Jiang1, LIAO Yan2, XIE Yue1, YANG Guangyou1, GU Xiaobin1,*
1.College of Veterinary Medicine, Sichuan Agricultural University, Chengdu 611130, China
2. Meishan Animal husbandry Station, Meishan 620000, China
Abstract

To explore the genetic variation and genetic taxonomy of Psoroptes ovis var. cuniculi in China, the complete sequence of mitochondrial ATP6 gene of P. ovis var. cuniculi isolated from 5 geographical regions (Huabei, Huadong, Huazhong, Xibei, and Xinan) were amplified by PCR sequencing techniques, and evaluated the genetic diversity and population structure of P. ovis var. cuniculi population. Eighty-eight sequences of mitochondrial ATP6 were obtained which were 672 bp, including 41 haplotypes. The whole population of P.ovis var. cuniculi in China showed high genetic diversity ( π=0.014 94) and haplotype diversity ( Hd=0.925 81), and the genetic variation within the population held the dominant positon(97.92%). The genetic differentiation index among the five populations was low ( Fst=0.02081), and there was no obvious differentiation and frequent gene exchange ( Nm=11.763 5). In P. ovis var. cuniculi population, the average values of Tajima’s D and Fu’s Fs were positive (Tajima’s D=0.937 98, Fu’s Fs=0.522 06) and there was no significant differentiation( P>0.05). The mismatch distribution curve showed multiple peaks, indicating that the population of P. ovis var. cuniculi was relatively stable in the evolutionary process and had no population expansion events. From NJ phylogenetic tree and haplotype network map, it can be seen that there was no genetic structure based on rabbit breed, temperature zone or geographical distribution. The results revealed that there was high level of genetic diversity, little population genetic differentiation, but population structure has not formed according to rabbit breed, temperature zone or geographical distribution.

Keyword: Psoroptes ovis var. cuniculi; mitochondrial ATP6 gene; genetic diversity

痒螨病是由痒螨科(Psoroptidae)、痒螨属(Psoroptes)的绵羊痒螨(P. ovis)寄生于绵羊、山羊、兔、水牛等家养动物及麋鹿等野生动物体表[1, 2], 引起的以结痂、脱毛、瘙痒为主要特征的慢性传染性皮肤病[3, 4, 5], 可造成养殖业重大的经济损失, 带来严重的动物福利问题, 同时威胁野生动物的生存[6]

中国养兔业发展迅速, 家兔饲养量达7亿只以上, 占全球饲养总量的40%以上[7]。痒螨病是危害我国养兔业的一种常见外寄生虫病, 在我国23个省中流行, 死亡率高达70%[8]。目前, 对我国各地绵羊痒螨(兔亚种)种群的种群变异、种群结构域基因交流信息的了解较少, 而这些信息对于研究绵羊痒螨(兔亚种)的动态历史及对养兔业中痒螨病的控制与消除是非常关键的。线粒体基因因其遗传稳定性好、进化速率快, 可用于痒螨遗传分化、遗传结构和系统发生关系研究[9]。本研究基于线粒体ATP6基因全序列分析, 旨在探讨我国绵羊痒螨(兔亚种)的遗传变异程度及遗传多样性。

1 材料与方法
1.1 绵羊痒螨(兔亚种)样品来源

从我国各地自然感染痒螨病的獭兔和肉兔耳部采集结痂, 收集在洁净的平皿中, 于37 ℃培养箱放置30 min, 除去痂皮, 收集虫体并鉴定[10], 最终获得88个痒螨虫株(表1), -20 ℃保存备用。

表1 我国绵羊痒螨(兔亚种)不同种群样品采集情况表 Table 1 Sampling sites of Psoroptes ovis var. cuniculi populations in China
1.2 引物设计

从绵羊痒螨(兔亚种)的线粒体基因组(GenBank No. KJ957822.1)获得ATP6基因全序列, 利用Primer Premier 5设计引物:ATP6-F: 5'-AGAGTGTAAAAACTAACTGTTCCAA-3'; ATP6-R: 5'-TAATCACCCCCAAAGCAAAA-3'。引物由生工生物工程(上海)股份有限公司合成。

1.3 单个虫体DNA的提取

取出-20 ℃保存的88份痒螨样品, 参考文献[11]的方法, 采用OMEGA公司Mollusc/Arthropod基因组DNA抽提试剂盒进行单个虫体基因组DNA提取, -20 ℃保存备用。

1.4 PCR扩增及测序

PCR扩增体系(25 μ L):Mixture 12.5 μ L, 上下游引物各1.0 μ L, 模板DNA 2.0 μ L, 双蒸水8.5 μ L。PCR反应条件:94 ℃ 5 min; 94 ℃ 45 s, 55 ℃ 45s, 72 ℃ 45 s, 35个循环; 72 ℃ 10 min。

扩增结束后, 将PCR产物进行1.0%琼脂糖凝胶电泳检测, 将条带大小符合预期的胶块切下, 经胶回收(天根普通琼脂糖凝胶DNA回收试剂盒, 编号:P4418)、胶回收产物连pMD19-T载体、转化DH5α 感受态细胞、蓝白斑筛选出阳性克隆, 交生工生物工程(上海)股份有限公司进行正反向测序。

1.5 序列分析

1.5.1 遗传多样性分析

测序序列经DNAMAN 5.2软件拼接、比对, 经DNASTAR 5.1软件剪切[12], 通过MEGA5.0软件计算核酸分歧度(nucleotide divergence, Dxy)、平均核苷酸差异指数(average index of nucleotide difference, K)和遗传距离[13], 运用DnaSP 5.1软件计算核苷酸多样性(nucleotide diversity, π )、单倍型多样性(haplotype diversity, Hd)和单倍型数(number of haplotypes, H)[14], 运用Arlequin 3.1软件作分子方差分析(AMOVA)计算出种群间遗传差异指数(Fst), 基因流(Nm)经下列公式计算:Nm=1/4(1/Fst-1)[15]

1.5.2 遗传结构分析

利用Arlequin 3.5软件计算中性检验值(neutrality tests)结果(Tajima’s D test, Fu’s FS test)[15], DnaSP 5.0软件分析错配分布曲线[14]。以欧洲屋尘螨(Dermatophagoides pteronyssinus, GenBank No. GQ469891.1)为外群, MEGA 5.1软件构建NJ系统发育树[13], 为检验聚类树置信度bootstraps 1 000次, 同时结合Network 5.0软件作单倍型网络图[16]

2 结果与分析
2.1 绵羊痒螨(兔亚种)ATP6基因全序列的碱基组成及变异

88个绵羊痒螨(兔亚种)样品全部扩增成功, 且测序获得ATP6基因全序列, 88条ATP6序列长度均为672 bp, 碱基平均组成由高到低分别是:T(45.53%)> A(22.92%)> G(17.82%)> C(13.73%), 表现出明显的AT偏倚性(68.45%)。88条ATP6序列共有60个变异位点, 其中简约信息位点28个, 单变异位点32个。

2.2 种群遗传多样性

表2所示, 我国5大地区绵羊痒螨(兔亚种)种群的单倍型数(H)为41, 单倍型多样性(Hd)范围为0.828 57(华中种群)~0.964 91(西北种群), 种群整体单倍型多样性(Hd)为0.925 81, 表明种群单倍型多样性丰富。5个种群中, 核苷酸多样性(π )最低的是华东种群(0.005 82), 西北种群最高(0.016 96), 种群整体核苷酸多样性表现较高水平(0.014 94)。我国5大地区种群的平均核苷酸差异指数(K)范围为8.166 67(华北种群)~11.397 66(西北种群)(表2)。

表2 绵羊痒螨(兔亚种)种群遗传多样性参数 Table 2 Genetic diversity indexes of P.ovis var. cuniculi populations in China
2.3 种群分化和遗传结构

表3所示, 我国西北种群与西南种群遗传距离最大0.016 47, 且差异极显著(P< 0.001); 华北种群与华中种群遗传距离最小0.012 28, 且差异极显著(P< 0.001); 核酸分歧度范围为0.014 49(华东种群与华北种群)~0.017 63(西北种群与西南种群), 整体表现与遗传距离的大小相一致。AMOVA分子变异分析显示, 我国绵羊痒螨(兔亚种)种群变异来源分别为种群内部(97.92%)与种群间(2.08%), 可见种群内分子遗传变异主要来自于种群内部。

表3 我国绵羊痒螨(兔亚种)各地理种群间的核酸分歧度(上三角)与遗传距离(下三角) Table 3 Pairwise comparisons of the nucleotide divergence (above triangle) and genetic distance (below triangle) between geographical populations of P. ovis var. cuniculi in China

表4所示, 我国5个地理种群间的Fst值为-0.044 52~0.106 12, 以华北与西南种群间的遗传分化最高(Fst=0.106 12), 西南、华北种群间基因交流程度最低(Nm=2.105 83)。我国绵羊痒螨(兔亚种)种群的总Fst值为0.020 81, 总基因流11.763 5, 表明我国种群内存在频繁的基因交流。

表4 我国绵羊痒螨(兔亚种)各地理种群间基因流(上三角)与遗传分化度(下三角) Table 4 Gene flow (above triangle) and genetic variation index values (below triangle) among 5 populations of P.ovis var. cuniculi in China
2.4 群体历史动态分析

中性检验结果见表2, 所有种群的Tajima’s D值均为正数。华中、西北种群的Fu’s Fs值为负值(-1.138 47和-2.714 84), 且差异不显著(P> 0.05), 其余3个种群Fu’s Fs值均为正值, 说明5个种群在历史进化过程中均没有经历过种群扩张。种群整体的Tajima’s D值和Fu’s Fs值均为正值(0.937 98和0.522 06), 且差异都不显著(P> 0.05)。核苷酸不配对差异分布曲线图具有多峰现象(图1), 说明我国5个绵羊痒螨(兔亚种)种群在进化过程中稳定, 无种群扩张事件。

图1 基于线粒体ATP6基因的绵羊痒螨(兔亚种)种群核苷酸不配对差异分布曲线图Fig.1 The distribution curve of nucleotide mismatch in P.ovis var. cuniculi population based on mitochondrial ATP6 gene

2.5 单倍型网络图及进化树分析

从构建的NJ树可见(图2), 外群欧洲屋尘螨(Dermatophagoides pteronyssinus, GenBank No. GQ469891.1)明显区分于绵羊痒螨(兔亚种)种群单独成一支, 而绵羊痒螨(兔亚种)虫株被分为2个大支, 且5个种群的单倍型在2个大支中都有分布。

图2 基于线粒体ATP6基因单倍型构建的系统发育树
CC, 华中地区; NWC, 西北地区; NC, 华北地区; EC, 华东地区; SWC, 西南地区。
Fig.2 Neighbor Joining (NJ) tree developed by haplotypes based on mitochondrial ATP6 gene
CC, Central China; NWC, Northwest China; NC, North China; EC, East China; SWC, Southwest China.

单倍型拓扑结构表明, 5个种群未出现基于地理位置的聚集, 即绵羊痒螨(兔亚种)没有形成基于地理分布的遗传变化(图2)。由单倍型网络图可知(图3), 大部分种群单倍型并未对应地理位置, 所有单倍型以星状网络结构进行分布共享, 大部分单倍型以5个单倍型(T1、T6、T7、T16、T24)为中心辐射分布, 其中T6可能是最为原始的单倍型, 并进化形成其它单倍型。整体单倍型分布聚类情况与NJ进化树类似, 说明我国绵羊痒螨(兔亚种)种群间未形成基于兔品种、温度带或地域来源的种群结构。

图3 基于线粒体ATP6基因构建的单倍型网络图
圆圈的大小与单倍型的分布频率成正比。
Fig.3 Haplotype network map based on mitochondrial ATP6 gene
The circle sizes were proportional to the distribution frequency of haplotype.

3 讨论

寄生虫的生活史、宿主种类、生活习性、外环境、人类活动等因素均会影响其种群变异与种群结构[17, 18, 19]。绵羊痒螨作为广宿主寄生的全球性流行外寄生虫, 其种群变异与种群结构研究非常有限。mtDNA以其严格的母系遗传和较快的进化速率[20], 能够全面地反映种群间和种群内的遗传差异[21, 22]。母系遗传可促使有利的突变迅速地形成、分离和表达, 并迅速在该群体中得到富集, 有利于探查自然选择在mtDNA留下的印记[23]。另外, ATP6、Cyt b基因对自然选择的作用较为敏感[24], 已见有学者将Cyt b基因用于绵羊痒螨(兔亚种)种群的相关研究[25]。因此, 本实验基于绵羊痒螨线粒体ATP6基因全序列探究我国绵羊痒螨(兔亚种)种群的遗传多样性, 以期为痒螨病的合理防控提供基础资料。

截止2018年12月, GenBank中未见绵羊痒螨(兔亚种)的ATP6基因序列, 故而本研究仅对我们扩增得到的88条ATP6序列进行分析。研究发现, 绵羊痒螨(兔亚种)具有较高的单倍型多样性, 且种群整体的遗传多样性水平较高。88个痒螨样本来自于我国7个省的10个采样点, 而这10个采样点当地的平均气候、降雨量、光照等自然因素均有差异, 由此可见, 高单倍型多样性和高核苷酸多样性可能有利于各地的绵羊痒螨(兔亚种)发生进化, 从而更好地适应当地环境而存活下来[26]

分析种群遗传距离可知, 种群间的遗传距离和地理距离未表现正相关关系。采用固定指数法度量种群间的遗传分化程度[27, 28]时, 我国绵羊痒螨(兔亚种)种群存在广泛的基因交流, 且遗传分化较弱。本研究的各地种群间地理距离较大, 且各种群间有山川、河流等天然屏障隔离, 单靠痒螨虫体的移动来完成不同地理种群间的基因交流较难。绵羊痒螨个体较小, 其离开宿主体后最长能存活18.8 d[29], 痒螨离开宿主后, 5 min内的最长运动距离不超过3 cm[30]

本研究中, 长距离痒螨种群间的频繁基因交流可能与以下因素有关:(1)我国各省间兔的频繁调运。以四川省为例, 四川肉兔的消费占全国30%, 每年需要从山东、陕西、云南等地调入兔只, 尤以山东调入最多, 这种频繁的调运导致四川与山东两地兔体寄生的痒螨群体间进行基因交流的机会增大, 从而可能导致西南种群与华东种群的基因交流趋于无穷大。(2)我国养殖户的养兔模式和痒螨的生活史。我国养兔场多以兔自然交配的形式繁殖后代, 兔的交配增加了不同痒螨群体的接触机会, 从而增加群体间的基因交流机会; 同时, 痒螨属于终生寄生虫, 其整个生活史均在宿主体完成, 这种呈群居生活且同一群体内交配繁殖的特性, 可增加群体间的基因交流, 降低分化度, 这点在线虫中亦得到证实[17]。(3)绵羊痒螨宿主范围广泛。已有研究证实, 分离自兔的痒螨可感染牛, 分离自绵羊的痒螨可以感染牛和家兔[31, 32], 这种宿主的广谱性可增加基因流动与交流。

我国痒螨整体种群的中性检测Tajima’s D和Fu’s Fs检验结果均为不显著的正值, 歧点分布为明显的双峰分布, 表明中国地区的绵羊痒螨种群在长期的历史上一直处于较为稳定的状态, 种群早期可能并未经历过群体扩张[33]。NJ进化树显示, 我国绵羊痒螨各种群的单倍型并未基于地理来源聚类, 而是混杂分布, 未形成明显的系统地理格局, 这点在单倍型网络图中亦得到进一步的证实。由此可见, 我国绵羊痒螨(兔亚种)种群未形成基于兔品种、温度带或者地理来源的种群遗传结构, 与赵习彬等[25]研究的结论一致, 这可能与痒螨在我国养兔业中的高流行率及虫体高密度寄生有关[8]。痒螨在养兔业中广泛的流行率和虫体高密度的寄生方式导致种群遗传结构的缺乏, 这点在微细毛圆线虫(Trichostrongylus tenuis)的种群遗传研究中得到证实[34]

参考文献
[1] ZAMAN M A, REHMAN T U, ABBAS R Z, et al. Therapeutic potential of ivermectin, doramectin and trichlorophan against Psoroptes ovis in sheep and cattle of Cholistan[J]. Pakistan Veterinary Journal, 2017, 37(2): 233-235. [本文引用:1]
[2] 王保健, 朱俊杨, 杨光友, . 痒螨抗原研究进展[J]. 动物医学进展, 2015, 36(7): 99-103.
WANG B J, ZHU J Y, YANG G Y, et al. Advance in antigen of Psoroptes mites[J]. Progress in Veterinary Medicine, 2015, 36(7): 99-103. (in Chinese with English abstract) [本文引用:1]
[3] 贾小勇, 杨光友, 古小彬, . 兔痒螨和水牛痒螨第二转录间隔区(ITS-2)基因序列分析[J]. 昆虫学报, 2008, 51(5): 545-550.
JIA X Y, YANG G Y, GU X B, et al. Nucleotide variation in rDNA ITS-2 differentiates Psoroptes isolates from rabbits and buffalo[J]. Acta Entomologica Sinica, 2008, 51(5): 545-550. (in Chinese with English abstract) [本文引用:1]
[4] 付雪红, 张常锁, 李静, . 绵羊痒螨的生物学研究进展[J]. 中国畜牧兽医, 2007, 34(6): 82-85.
FU X H, ZHANG C S, LI J, et al. Biological progress on Psoroptes ovis[J]. China Animal Husband ry & Veterinary Medicine, 2007, 34(6): 82-85. (in Chinese) [本文引用:1]
[5] BURGESS S, BARTLEY K, MARR E, et al. Draft genome assembly of the sheep scab mite, Psoroptes ovis[J]. Genome Announcements, 2018, 6(16): e00265-18. [本文引用:1]
[6] 刘学波. 一例兔痒螨病的诊治报告[J]. 畜牧兽医科技信息, 2016(4): 127.
LIU X B. Diagnosis and treatment of in Psoroptes ovis in rabbits[J]. Chinese Journal of Animal Husband ry and Veterinary Medicine, 2016(4): 127. (in Chinese) [本文引用:1]
[7] ANDRIANTSOANIRINA V, ARIEY F, IZRI A, et al. Sarcoptes scabiei mites in humans are distributed into three genetically distinct clades[J]. Clinical Microbiology and Infection, 2015, 21(12): 1107-1114. [本文引用:1]
[8] ZHENG W, TANG Q, ZHANG R, et al. Vaccination with recombinant actin from scab mites and evaluation of its protective efficacy against Psoroptes cuniculiin fection[J]. Parasite Immunology, 2013, 35(2): 91-98. [本文引用:2]
[9] WANG S, GU X B, FU Y, et al. Molecular taxonomic relationships of Psoroptes and Chorioptes mites from China based on COI and 18S rDNA gene sequences[J]. Veterinary Parasitology, 2012, 184(2/4): 392-397. [本文引用:1]
[10] 陈志蓉, 连宏军, 徐雪平, . 绵羊痒螨发育分期观察及生物学特性研究[J]. 中国兽医科技, 2000, 30(7): 5-7.
CHEN Z R, LIAN H J, XU X P, et al. Biological features and observation of development of Psoroptes communis Ovis[J]. Chinese Journal of Veterinary Science and Technology, 2000, 30(7): 5-7. (in Chinese) [本文引用:1]
[11] 李燕方, 古江, 袁媛, . 基于线粒体12S基因全序列分析我国绵羊痒螨(兔亚种)的遗传多样性[J]. 畜牧兽医学报, 2018, 49(11): 2468-2476.
LI Y F, GU J, YUAN Y, et al. Genetic diversity of Psoroptes ovis var. cuniculi determined by the full-length of mitochondrial 12S gene in China[J]. Chinese Journal of Animal and Veterinary Sciences, 2018, 49(11): 2468-2476. (in Chinese with English abstract) [本文引用:1]
[12] BURLAND T G. DNASTAR’s lasergene sequence analysis software[J]. Methods in Molecular Biology, 2000, 132: 71-91. [本文引用:1]
[13] TAMURA K, PETERSON D, PETERSON N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods[J]. Molecular Biology and Evolution, 2011, 28(10): 2731-2739. [本文引用:2]
[14] LIBRADO P, ROZAS J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics, 2009, 25(11): 1451-1452. [本文引用:2]
[15] EXCOFFIER L, LISCHER H E L. Arlequin suite ver 3. 5: a new series of programs to perform population genetics analyses under Linux and Windows[J]. Molecular Ecology Resources, 2010, 10(3): 564-567. [本文引用:2]
[16] 宋森, 王小立, 周蓉, . 基于线粒体DNA控制区的石鸡华北亚种的种群历史动态研究[J]. 生态学报, 2015, 35(2): 280-289.
SONG S, WANG X L, ZHOU R, et al. Demographic history of Alectoris chukar pubescens based on mitochondrial DNA control region[J]. Acta Ecologica Sinica, 2015, 35(2): 280-289. (in Chinese with English abstract) [本文引用:1]
[17] COLE R, VINEY M. The population genetics of parasitic Nematodes of wild animals[J]. Parasites & Vectors, 2018, 11: 590. [本文引用:2]
[18] VAZQUEZ-PRIETO S, VILAS R, PANIAGUA E, et al. Influence of life history traits on the population genetic structure of parasitic helminths: a minireview[J]. Folia Parasitologica, 2015, 62. [本文引用:1]
[19] BARRETT L G, THRALL P H, BURDON J J, et al. Life history determines genetic structure and evolutionary potential of host-parasite interactions[J]. Trends in Ecology & Evolution, 2008, 23(12): 678-685. [本文引用:1]
[20] RUIZ-PESINI E, LOTT M T, PROCACCIO V, et al. An enhanced MITOMAP with a global mtDNA mutational phylogeny[J]. Nucleic Acids Research, 2007, 35: 823-828. [本文引用:1]
[21] LÓPEZ-LÓPEZ A, VOGLER A P. The mitogenome phylogeny of Adephaga( Coleoptera)[J]. Molecular Phylogenetics and Evolution, 2017, 114: 166-174. [本文引用:1]
[22] WANG A R, KIM M J, JEONG S Y, et al. Complete mitochondrial genome sequence of Cicindela anchoralis chevrolat, 1845 (Coleoptera: Carabidae)[J]. Mitochondrial DNA Part B, 2018, 3(1): 282-283. [本文引用:1]
[23] FAY J C, CHUNGI W. Hitchhiking under positive Darwinian selection[J]. Genetics, 2000, 155(3): 1405-1413. [本文引用:1]
[24] MISHMAR D, RUIZ-PESINI E, GOLIK P, et al. Natural selection shaped regional mtDNA variation in humans[J]. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(1): 171-176. [本文引用:1]
[25] 赵习彬, 王保健, 简克灵, . 基于线粒体细胞色素氧化酶b(cytb)全基因序列分析我国部分地区兔痒螨的遗传变异特征[J]. 畜牧兽医学报, 2017, 48(4): 714-721.
ZHAO X B, WANG B J, JIAN K L, et al. Analysis of the characteristics of genetic variability within Psoroptes cuniculi isolates in some regions of China, inferred by mitochondrial cytochrome oxidase b gene[J]. Chinese Journal of Animal and Veterinary Sciences, 2017, 48(4): 714-721. (in Chinese with English abstract) [本文引用:2]
[26] 任旭琴. 遗传多样性及其研究方法[J]. 淮阴工学院学报, 2002, 11(5): 6-8.
REN X Q. Genetic diversity and the methods[J]. Journal of Huaiyin Industry College, 2002, 11(5): 6-8. (in Chinese with English abstract) [本文引用:1]
[27] CHAKRABORTY R, DANKER-HOPFE H. Analysis of population structure: a comparative study of different estimators of wright’s fixation indices[J]. Hand book of Statistics, 1991, 8(8): 203-254. [本文引用:1]
[28] SLATKIN M. Gene flow and the geographic structure of natural populations[J]. Science, 1987, 236(4803): 787-792. [本文引用:1]
[29] THIND B B, FORD H L. Maintaining Psoroptes ovis(Acari: Psoroptidae) off-host in the laboratory[J]. Experimental and Applied Acarology, 2007, 43(3): 205-211. [本文引用:1]
[30] PEGLER K R, WALL R. Tactic responses of the parasitic mite, Psoroptes ovis to light and temperature[J]. Experimental and Applied Acarology, 2004, 33(1-2): 69-79. [本文引用:1]
[31] KIRKWOOD A C. Some observations on the biology and control of the sheep scab mite Psoroptes ovis(Hering) in Britain[J]. Veterinary Parasitology, 1985, 18(3): 269-279. [本文引用:1]
[32] MELENEY W P. Experimentally induced bovine psoroptic acariasis in a rabbit[J]. American Journal of Veterinary Research, 1967, 28(124): 892-894. [本文引用:1]
[33] 谢跃, 周璇, 孙韵, . 基于线粒体12S对中国两大山系大熊猫西氏贝蛔虫种群遗传多样性的分析[J]. 兽类学报, 2015, 35(3): 328-335.
XIE Y, ZHOU X, SUN Y, et al. Mitochondrial 12S-based analysis of genetic diversity of Baylisascaris schroederi in giant pand as from two mountain ranges in China[J]. Acta Theriologica Sinica, 2015, 35(3): 328-335. (in Chinese with English abstract) [本文引用:1]
[34] SEIVWRIGHT L J, REDPATH S M, MOUGEOT F, et al. Faecal egg counts provide a reliable measure of Trichostrongylus tenuis intensities in free-living red grouse Lagopus lagopus scoticus[J]. Journal of Helminthology, 2004, 78(1): 69-76. [本文引用:1]