Genetic diversity of Bangladeshi native chickens based on complete sequence of mitochondrial DNA D-loop region

ABSTRACT

1. The aim of this study was to explore genetic diversity and possible origin of Bangladeshi (BD) native chickens. The complete mtDNA D-loop region was sequenced in 60 chickens representing five populations; naked neck, full feathered, Aseel, Hilly and autosomal dwarf. The 61 reference sequences representing different domestic chicken clades in China, India, Laos, Indonesia, Myanmar, and other Eurasian regions were included. The mtDNA D-loop sequence polymorphism and maternal origin of five BD populations were analysed.

2. A total of 35 polymorphic sites, and 21 haplotypes were detected in 60 mtDNA D-loop sequences. The haplotype and nucleotide diversity of the five populations were 0.921 ± 0.018 and 0.0061 ± 0.0019, respectively. Both mtDNA network and phylogenetic analysis indicated four clades (four haplogroups) in BD populations (21 haplotypes) along with 61 reference haplotypes. Clade E contained the most individuals (20) and haplotypes (11) of BD chickens, followed by clade D (17, 6), clade C (12, 2) and clade F (11, 2), respectively.

3. The higher number of unique haplotypes found in Yunnan, China, suggested that the origin of BD chickens was in this region. The haplotypes from different haplogroups were introduced in Bangladeshi chickens from India, China and Myanmar. The phylogenetic tree showed a close relationship of BD chickens with the clusters from India, China, Myanmar and Laos, and indicated the dispersion of BD chickens from these sources. The phylogenetic information revealed high genetic diversity of BD chickens because of their origin from different lineages with high genetic variation and distance, which was determined from four cluster and neighbour-joining trees.

4. In conclusion, BD populations had high genetic diversity. The mtDNA network profiles and phylogenetic trees showed multiple maternal origins of BD chickens from India, China, Myanmar and Laos.     

基于线粒体控制区的鸡不同杂交组合遗传多样性研究

旨在探讨鸡不同杂交组合线粒体控制区（mtDNA D-loop区）的遗传多样性和单倍型特性。选取固始鸡和隐性白羽鸡及其正、反交F1代、藏鸡以及F2代等6个群体共387个个体的mtDNA D-loop区进行测序,分析其遗传规律和单倍型特性,并与不同红色原鸡亚种进行聚类,分析其母系起源。结果显示,6个群体D-loop区全序列大小为1 231 bp,共检测到28个多态位点和1个C碱基缺失,共构成19种单倍型,分为A、B、C和E 4个单倍型群,其中,固始鸡和反交F1代主要为A、C单倍型,固始鸡A、C单倍型比例分别为53.42%和46.58%,反交F1代A、C单倍型比例分别为50.75%和49.25%;隐性白羽鸡、正交F1代和F2代优势单倍型均为E单倍型,占比分别为48.89%、48.84%和50.00%。6个鸡群体单倍型多样度（Hd）在0.496~0.729之间,核苷酸多样度（Pi）在0.003 40~0.005 41之间,Hd值和Pi值最大的均为正交F1代,其次为隐性白羽鸡和F2代,固始鸡和反交F1代群体遗传多样性接近。聚类分析显示,A、B单倍型群与滇南亚种交叉聚为一枝;E单倍型群与印度亚种交叉聚为一枝;C单倍型群与印度亚种、指名亚种、印尼亚种以及滇南亚种聚为一枝。结果提示,mtDNA D-loop区遵循严格的母系遗传,后代的遗传多样性和单倍型比例与其母本基本一致;我国家鸡群体具有多个红色原鸡母系起源,且主要起源于原鸡滇南亚种。     

Genetic Diversity of Lueyang Black-bone Chicken Based on Mitochondrial DNA D-loop Region Sequence

This study was aimed to assess mitochondrial DNA D-loop sequence diversity and origin of Lueyang Black-bone chicken.The mtDNA D-loop sequence of 30 individuals from Lueyang Black-bone chicken were amplified by PCR and subsequently sequenced.The mtDNA D-loop sequence of other chicken were collected from GenBank and used as reference sequences to analyze the diversity and origin of Lueyang Black-bone chicken.The results revealed that the average values of base composition of A,C,G and T in mtDNA D-loop the sequence of Lueyang Black-bone chicken were 26.6%,26.6%,13.4% and 33.4%,respectively.26 nucleotide polymorphic sites were transition.The average nucleotide diversity (Pi) of the sites and haplotype diversity (Hd) were 0.00705 and 1.000,and the value of Tajima's D was -0.47272.Phylogenetic tree showed that samples were clusted in 4 clades. In the research, it could be concluded that the genetic diversity was relatively rich and wealthy and there were 4 maternal origins to Lueyang Black-bone chicken population.     

略阳乌鸡线粒体DNA控制区(D-loop)的遗传多样性分析

为了研究略阳乌鸡线粒体DNA控制区（mtDNA D-loop）的遗传多样性和起源,本研究对30只略阳乌鸡样品的mtDNA D-loop全序列进行了PCR扩增和测序,结合GenBank中公布的其他鸡的D-loop区序列,分析略阳乌鸡线粒体多态性及其起源。结果表明,略阳乌鸡mtDNA D-loop区全序列中,A、C、G、T平均含量分别为26.6%、26.6%、13.4%和33.4%,26个核苷酸多态位点均为转换位点,核苷酸多样度（Pi）为0.00705,单倍型变异度（Hd）为1.000,中性检验Tajima's D值为-0.47272。通过群体构建的系统进化树发现,略阳乌鸡样品在系统进化树上聚为4大分支。研究结果表明,略阳乌鸡群体内个体序列变异程度较大,遗传多样性丰富,揭示略阳乌鸡在遗传组成上具有4个母系来源。     

建昌马mtDNA D-loop区遗传多样性及系统进化分析

试验旨在以线粒体DNA(mitochondrial DNA,mtDNA)为切入点,研究建昌马的母系遗传多样性与系统进化。从建昌马(n=39)血液中提取基因组DNA,用PCR方法扩增mtDNA D-loop区并直接测序,分析其高变区247 bp序列信息,统计mtDNA D-loop区的单倍型及变异位点,计算单倍型多样性(haplotype diversity,Hd)、核苷酸多样性(nucleotide diversity,Pi)和平均核苷酸变异数(average number of nucleotide differences,K)。构建包括建昌马在内的19个品种马的NJ系统进化树,计算各品种间的遗传距离。结果显示,试验获得了清晰的PCR扩增产物,并通过直接测序方法获得了约1200 bp的序列。39匹建昌马mtDNA D-loop区247 bp序列(其中1个样品缺失1 bp)的AT碱基含量为61.45%,属AT碱基对富集区,检测到33个多态性位点,共显示26种单倍型,其中4种为共享单倍型,且Hap7和Hap1为优势单倍型,单倍型多样性为0.947,核苷酸多样性为0.02399,平均核苷酸变异数为5.901,显示丰富的母系遗传多样性;NJ系统进化树显示,建昌马分布在A、C、D、E、F、G共6个支系中,约50%的样品分布在A支系,显示出复杂的母系起源;建昌马与关中马的遗传距离最小(0.021),其次是三河马、文山马、韩国车巨马(0.024),与韩国济州岛马遗传距离最大(0.032)。本研究结果表明,建昌马的mtDNA D-loop高变区遗传多样性丰富,具有多个母系起源,且A支系占有明显优势,与关中马、文山马可能有共同的母系起源。     

Genetic Diversity and Evolution of Danzhou Chicken Assessed with Mitochondrial Complete DNA D-loop Region Sequencing

The objective of this study was to determine the genetic diversity and evolution of Danzhou chicken.The complete mitochondrial DNA (mtDNA) D-loop regions of 36 Danzhou chickens were amplified,sequenced and analyzed.The sequencing reads were compared with the complete mtDNA D-loop sequence of several relative strains of chicken annotated in GenBank,and analyzed by bioinformatics methods.The genetic diversity and its evolutionary relationship in Danzhou chicken were analyzed.The results showed that the lengths of PCR products at the D-loop region were 1 210 bp,with 59.9% being A+T and 40.1% as C+G.The variable regions were 167-1 215 bp,and the high variable regions were mainly 167-367 bp.A total of 20 variable sites that defined 6 haplotypes were identified.The average haplotype diversity (Hd) and average number of nucleotide difference (k) were 0.571 and 6.449,respectively,the nucleotide diversity (Pi) was 0.00537,and the Tajima's D value of neutrality test was 1.61643.6 haplotypes could be grouped to 3 haplogroups (A,B and C) as determined by phylogenic analysis,with B clade,as the most abundant population.It concluded that the genetic diversity and haplotype diversity of Danzhou chicken were relatively low.Phylogenetic tree showed that the genetic composition of Danzhou chicken came from 3 maternal ancestors,Gallus gallus spadiceus,Gallus gallus bankiva and Gallus gallus jabouillei were potential ancestors.There was few influence of exotic lineage detected,which indicated that Danzhou chicken was a relatively conserved breed.     

4个引进山羊品种mtDNA控制区序列变异和系统发生关系研究

本研究测定了四川4个引进山羊品种24个个体的线粒体控制区全序列,并从GenBank获得山羊属2个野山羊种的2每控制区序列。利用MEGA2．0软件构建分子系统发育无根树。序列分析表明:山羊控制区线粒体控制全序列长度为1 212bp或1 213 bp,A+T含量占59．9％,其中64个核苷酸位点存在变异(约占5．28％),核苷酸多样度为 1．731％,这些差异共定义了16种单倍型,单倍型多样性为0．913±0．048。安哥拉山羊、波尔山羊都有自己独特的单倍型,与其他品种间都没有共享类型。使用NJ法构建了系统发育树,结果表明:2个野山羊种中,角(?)羊与家养山羊的关系相对较近,4个家养山羊品种有2个母系来源,在4个家养山羊品种内部,安哥拉山羊的分化要比其他3个品种早。     

唐古拉山牦牛群体的母系遗传多样性及遗传结构分析

[目的]从分子水平上探究青海省唐古拉山牦牛群体的母系遗传多样性、群体遗传结构及其遗传背景。[方法] 对52头唐古拉山牦牛个体mtDNA D-loop区序列进行测定后,使用生物信息学软件分析确定其核苷酸变异位点和单倍型数目,计算单倍型多样度和核苷酸多样度大小,并进行系统发育分析。[结果] 在619 bp唐古拉山牦牛D-loop区序列分析中,排除2处插入（缺失）后共检测到31处多态位点,包括单一多态位点5处和简约信息位点26处。根据序列间核苷酸变异共确定了13种单倍型,单倍型多样度和核苷酸多样度分别为0.821±0.043和0.007±0.004。与我国其他18个家牦牛品种和野牦牛相比,唐古拉山牦牛群体单倍型多样度和核苷酸多样度值均较低,表明该群体遗传变异较为贫乏,母系遗传多样性水平较低。以美洲野牛为外群,邻接法（即NJ法）构建的系统发育树结果显示：唐古拉山牦牛群体13种单倍型分布在A、B、C、D和E五种单倍型组中,且聚为2个大的母系分支（即I和II）,支系Ⅰ占比为77%,提示唐古拉山牦牛由2个母系支系组成,拥有2个母系起源且以支系Ⅰ为主。 [结论] 唐古拉山牦牛母系遗传多样性水平较低,由2个母系支系组成,以支系Ⅰ为主,推测其有2个母系起源。     

Genetic diversity of mtDNA D-loop sequences in four native Chinese chicken breeds

1. To explore the genetic diversity of Chinese indigenous chicken breeds, a 585 bp fragment of the mitochondrial DNA (mtDNA) region was sequenced in 102 birds from the Xichuan black-bone chicken, Yunyang black-bone chicken and Lushi chicken. In addition, 30 mtDNA D-loop sequences of Silkie fowls were downloaded from NCBI. The mtDNA D-loop sequence polymorphism and maternal origin of 4 chicken breeds were analysed in this study.

2. The results showed that a total of 33 mutation sites and 28 haplotypes were detected in the 4 chicken breeds. The haplotype diversity and nucleotide diversity of these 4 native breeds were 0.916 ± 0.014 and 0.012 ± 0.002, respectively. Three clusters were formed in 4 Chinese native chickens and 12 reference breeds. Both the Xichuan black-bone chicken and Yunyang black-bone chicken were grouped into one cluster. Four haplogroups (A, B, C and E) emerged in the median-joining network in these breeds.

3. It was concluded that these 4 Chinese chicken breeds had high genetic diversity. The phylogenetic tree and median network profiles showed that Chinese native chickens and its neighbouring countries had at least two maternal origins, one from Yunnan, China and another from Southeast Asia or its surrounding area.     

Analyses of Genetic Diversity and Phylogeny of mtDNA D-loop from Yak in Karakoram-Pamir Area

This experiment was conducted to clarify the genetic diversity,genetic differentiation and phylogenetic status of yak in Karakoram-Pamir area.The mtDNA D-loop region sequence was selected as a molecular marker,and the sequence and genetic diversity of the mtDNA D-loop region of yak in Karakoram-Pamir area were analyzed by PCR direct sequencing and bioinformatics methods.The yak sequence in GenBank was used.The maximum likelihood method was used to construct the phylogenetic tree and the intermediary network relationship.The results showed that the mtDNA D-loop sequence of yak in Karakoram-Pamir area was rich in A and T bases,with AT content of 61.2%,and there were 63 polymorphic loci,accounting for 7.04% of the total number of nucleotides.The results indicated that A and T bases were rich in the mtDNA D-loop sequences at 61.2%.There were 63 mutation sites,accounting for 7.04% of all nucleotides,The average haplotype diversity (Hd) was 0.806,the average nucleotide diversity (π) was 0.01528,and the average nucleotide difference (K) was 13.509,indicating that the yak was rich in genetic diversity in Karakoram-Pamir area;Through phylogenetic analysis,there were two branches in yak in China,forming two branches and six small clades.The yak in Karakoram-Pamir area involved in this study had two different maternal origins.Additionally,yak in the Karakoram-Pamir area was less shared with other breeds of yak haplotypes.In the branch C,the yak group in the Karakoram-Pamir area accounts for a large proportion and was shared with wild yak.The yak population in Karakoram-Pamir area had a unique genetic background,which might be the result of early domestication of wild yaks.It was suggested to increase the identification of yak breeds and the formulation of breed standards in this area,and strengthen the protection of yak genetic resources in this area.According to the current situation of the population,wild blood yaks were introduced for purification and rejuvenation to prevent breed degeneration and decrease of genetic diversity.The introduction of foreign yak breeds and disorderly hybridization were reduced to ensure the characteristics of this breed of high-quality yak breed resources.     

不同类型肉鸡线粒体单倍型与遗传起源研究

旨在探讨鸡线粒体D-loop区单倍型特性以及与生长速度的相关性,并分析其遗传起源。选取不同类型肉鸡品种（配套系）6个,测定其生产性能,并对6个群体共计314个个体线粒体D-loop区全长进行测序,分析其单倍型特性;同时与不同红色原鸡亚种进行聚类,分析其母系起源。结果显示,6个群体D-loop区全序列共检测到37个突变位点,构成40种单倍型,分为A、B、C和E 4个单倍型群,其中AA肉鸡、罗斯308、禽雁麻鸡和裕禾1号肉鸡主要为E单倍型,占比分别为85.92%、50.00%、100.00%和70.21%;园丰麻鸡2号和港丰瑶黑麻鸡E单倍型占比相对较低,分别为15.00%和22.39%,园丰麻鸡2号B单倍型含量最高,占比66.67%,港丰瑶黑麻鸡C单倍型含量最高,占比40.30%。相关性分析显示,初生重与E单倍型比例之间呈极显著正相关（P<0.01）,与A、B和C单倍型比例之间均呈负相关;E单倍型比例与公母平均体重约1.8 kg时日龄之间呈显著负相关（P<0.05）,而与公母平均体重约1.8 kg时饲料转化比之间呈极显著负相关（P<0.01）。聚类分析显示,A、B单倍型所有个体均与原鸡滇南亚种聚为一枝;E单倍型所有个体均与原鸡印度亚种聚为一枝;C单倍型所有个体与原鸡印度亚种、滇南亚种、指名亚种以及印尼亚种交叉聚为一类。结果提示,线粒体单倍型与肉鸡生长速度之间相关显著,E单倍型与肉鸡生长速度具有较强的正相关;我国家鸡群体母系起源丰富,但主要起源于原鸡滇南亚种。     

中国地方猪种线粒体基因组遗传多样性及其系统发育分析

为从分子水平上探究中国地方猪种遗传多样性和分类地位,本试验采用生物信息学方法比较了6个类型共22个中国地方猪种的线粒体基因组全序列,分析了其多态性,并构建了6个类型猪种线粒体D-loop区单倍型的网络中介图以及基于线粒体D-loop序列、Cytb基因、完整编码区序列的系统进化树。结果表明,6个类型22个猪种中共检测到了144个多态位点,22种单倍型,说明地方猪种具有丰富的遗传多样性;地方猪线粒体基因组序列中核苷酸变异以转换为主,且Ti/Tv大于转换/颠换比临界值（2.0）,变异位点均符合中性突变。6个类型猪种间遗传距离均较小,且有共享单倍型。系统进化树结果表明,6种类型地方猪种主要聚为两个支系。表明线粒体D-loop序列及Cytb基因均可作为研究种内系统发育、起源进化的分子标记。     

Genetic Diversity Analysis of Shaanxi Moschus berezovskii Based on mtDNA Cytb Gene and D-loop Region Sequence

【Objective】 This study was aimed to explore the genetic diversity of Shaanxi Moschus berezovskii population,and understand the genetic information of Moschus berezovskii.【Method】 The hair of Moschus berezovskii was collected to extract DNA,the mitochondrial DNA(mtDNA) cytochrome b(Cytb) gene and D-loop sequences of 43 Moschus berezovskii individuals were determined,and the base composition was counted.All sequences were integrated and compared using ClustalX 2.0 software to obtain nucleotide polymorphic sites (SNPs) in the population.The nucleotide diversity (Pi),number of haplotype (H),haplotype diversity (Hd) and average number of nucleotide differences (K) were calculated by DNASP 5.10 software.The genetic distance among different haplotypes of Cytb gene and D-loop sequences was calculated by Mega 7.0 software,and Neighbor-Joining (NJ) phylogenetic tree was constructed.【Result】 The AT content of Cytb gene and D-loop region were higher than GC content,indicating there was bias in base composition.There were 241 and 383 SNPs of Cytb gene and D-loop region,respectively.The nucleotide diversity of Cytb gene and D-loop region were 0.28343 and 0.07707,and the haplotype diversity was 0.983 and 0.975,respectively,indicating that the population genetic diversity was rich.The genetic distances of 35 haplotypes of Cytb gene ranged from 0.002 to 0.831,and 29 haplotypes of D-loop region ranged from 0.006 to 1.342.The phylogenetic tree showed that there were two mitochondrial lineages,indicating that there were two mitochondrial maternal origins.The evolutionary analysis of D-loop region also supported this conclusion.【Conclusion】 The nucleotide diversity and haplotype diversity of Moschus berezovskii population were high,and the genetic diversity was rich.At the same time further supported the view of Moschus berezovskii and Moshus moschiferus belonged to a branch of the view.     

关中驴线粒体DNA D-loop多态性分析

本文对 6头关中驴线粒体DNAD loop区 399bp的核苷酸序列进行了分析。结果发现 ,关中驴的D loop区核苷酸变异只有转换 1种形式。 6头关中驴D loop区的核苷酸序列组成 3种单倍型 ,单倍型比例为 5 0 .0 0 % ,说明关中驴mtDNA遗传多样性正逐步丧失 ,需要加强驴种质资源保护。以欧洲驴D loop序列为对照 ,关中驴 3种单倍型的核苷酸变异率分别为 4 .2 1%、4 .5 1%和 0 .2 5 %。在获得的 399bpD loop碱基序列中 ,共检测出核苷酸多态位点18个 ,多态位点比例为 4 .5 1%。从D loop区核苷酸序列的 3种单倍型分析 ,发现关中驴可能有 2种不同的母系起源     

中国家鸡和红色原鸡mtDNA控制区遗传多态性及系统进化分析

通过线粒体DNA控制区的结构和多态性来研究中国家鸡和红色原鸡的遗传多态性与系统进化。测定14个中国地方鸡种和红色原鸡2个亚种的256个个体线粒体DNA控制区部分序列约560bp,结果表明,A、C、G、T这4种核苷酸的平均比例分别为25．∞％、37．40％、4．40％和33．20％。共发现44个变异位点,约占分析位点总数的7．86％,没有观测到插入／缺失,颠换和转换之比为0．13;共具有32种单倍型,9种为共享单倍型;16个群体内单倍型多样度从0到0．964,单倍型变异度总体为0．909&#177;0．014,整体的平均核苷酸差异数为7．276,核苷酸多样度为1．851％。群体间核苷酸分歧度（Dxy）在0．747％～3．125％之间变化,核苷酸净遗传距离（Da）为0．015％～2．633％。16个群体表现出较高水平的遗传多态性,群体间表现出显著的遗传分化。群体遗传多态性和亲缘关系分析表明,一些中国家鸡的群体（如固始鸡和仙居鸡）起源于泰国红色原鸡Gallus gallu sgallus亚种,一些中国家鸡的群体（如茶花鸡和藏鸡等）起源于中国红色原鸡Gallus gallus spadiceus亚种,在一些中国地方鸡种还同时具有这2种红色原鸡的遗传贡献;认为中国家鸡起源于泰国或单纯起源于中国的观点都是不全面的。     

Lineage classification of canine inheritable disorders using mitochondrial DNA haplotypes

To estimate the maternal effects of dog breeds using mitochondrial DNA(mtDNA) haplotypes in the dogs with several clinical disorders, 600 base pairs of mtDNA D-loop region were amplified from 365 dogs and were determined for mtDNA sequences. The diversity of the 600-bp sequences was classified into 64 haplotypes, including 46 newly discovered haplotypes, and the haplotypes were grouped into four clusters I to IV. Lineage analysis using the mtDNA haplotype indicated that each dog breed genetically comprises one or a few mtDNA haplotypes. When the relationship between genetic background and occurrences of clinical diseases was estimated, canine lineage analysis using mtDNA haplotype revealed that the disorders distributed in the dominant mtDNA haplotypes of each dog breed, but no disorder closely associated with mtDNA haplotypes was detected.     

4种鼢鼠线粒体DNA D-loop区全序列分析

采用PCR测序技术对采自甘肃境内的4种鼢鼠22个个体的线粒体DAND-loop区全序列进行了测定。结果表明：鼢鼠线粒体DNAD-loop序列长度为893、894、895、896或899bp,核苷酸位点突变类型有5种,即转换、颠换、插入、缺失及转换与颠换共存。碱基转换以T〈=〉C形式为主。其中T、C、A、G4种核苷酸的平均比例分别为34．1％、24．3％、30．1％、11．5％,A＋T含量（64．2％）明显高于G＋C含量（35．8％）。这4种鼢鼠22条线粒体DNAD-loop区发现20种单倍型,单倍型比例为81．82％,说明中国鼢鼠mtDNA遗传多态性很丰富。从系统发育树分析,4种鼢鼠明显聚为两类,推测鼢鼠可能有两个起源。     

贵州地方鸡种的遗传变异研究

为了研究贵州地方鸡种的遗传多样性与系统进化,试验通过PCR直接测序的方法,测定了6个贵州地方鸡种的90个个体线粒体DNA控制区部分序列,并进行了序列分析。结果表明,DNA控制区序列长度为780 bp,A、G、C、T 4种核苷酸的平均比例分别为26.14%、13.38%、25.98%、34.50%;共发现40个变异位点(不包含种内变异位点),占分析位点总数的5.13%;共确定30个单倍型,6个群体内单倍型变异度总体为0.914,变化范围为0.506～0.976。6个群体间表现出显著的遗传分化和较高水平的群体遗传多态性。     

浙江省本地及引进鹅种遗传多样性研究

本研究利用PCR和DNA克隆测序技术扩增并分析了浙江省5个地方鹅种和1个引进家鹅品种96个个体的线粒体DNA D环区(D-loop)全序列。结果表明:浙江省境内的家鹅线粒体DNA D-loop全序长为1 166～1 179 bp,共发现89个变异位点,定义了36种单倍型,群体总核苷酸多样度达到了0.007 25,单倍型多样度为0.862,说明鹅群的遗传资源较为丰富,其中磐石灰鹅的遗传多样性最丰富,朗德鹅的遗传多样性最匮乏,而且没有发现朗德鹅与其他本地鹅种之间的显著差异。同时,个体聚类图中同品种内的个体没有聚在一起,说明浙江省鹅品种的母系来源比较混杂,品种选育偏重于父系,应同时结合基因组水平对其遗传结构进行分析。     

深县猪mtDNA CytB基因的遗传多样性与母系起源研究

为探讨深县猪群体内的遗传多样性及母系起源分化,采用PCR直接测序法测定了40头深县猪的mtDNA CytB基因长1 140 bp片段的全序列,并结合GenBank已公布的15个猪种mtDNA CytB基因全序列,采用邻接法构建深县猪、部分其他地方猪种和引进猪种的系统发育树。结果显示,在1 140 bp长的序列中A、T、G、C的含量分别为25.93%、31.72%、28.90%及13.45%,其中A+T的含量(57.65%)高于G+C的含量(42.35%);共发现2个变异位点,无插入和缺失突变,全部为转换位点,其中简约信息位点2个,未发现单一信息位点。40条序列共定义了3种单倍型,单倍型多样性(HD)为0.312,核苷酸多样性(Pi)为0.00026,平均核苷酸差异数(K)为0.327。在深县猪mtDNA CytB基因全序列NJ进化树中,三种单倍型聚为同一分支,母系来源单一,与莱芜猪、大蒲莲猪等山东品种遗传距离较近。结果表明,深县猪群体内遗传多样性比较贫乏,与山东地区的华北型黑猪种群有更近的亲缘关系。