Population structure, genetic diversity and phylogenetic analysis of different rural and commercial chickens of Pakistan using complete sequence of mtDNA D-loop

Chicken is the most important poultry bird in Pakistan that not only provides nutrition but also contributes to country's economy. The Pakistani chicken and its germplasm resources are not genetically characterized and explored. Here, we focused at rural and commercial chickens of Pakistan and explored genetic diversity, population structure and phylogeny. We first collected feather samples from Rural and Broiler populations of Pakistani chickens, isolated DNA and sequenced complete D-loop of mtDNA. The length of complete D-loop ranged from 1231 to 1234 bp in Pakistani chickens. The GC content was 39%. Hotspots of mutations were three hypervariable sites (HVS). Most of the variations (77%) were in HVS1. In a total, 26 polymorphic sites defined 12 haplotypes and all major haplogroups (A-I) in genetic structure of Pakistani chickens. Genetic diversity remained relatively very low in Broiler (Pi = 0.00212 ± 0.00136). There was a low sharing of matrilineages between the two populations (Fst = 0.170). With high Hd value (0.825 ± 0.051) and presence of all nine major haplogroups the rural chicken population showed relatively rich genepool. Finally we did molecular phylogenetic analysis and inferred phylogeny. Presence of subcontinent specific haplogroups E3 and I and clustering of Indian red junglefowl closely with Pakistani chickens in Bayesian inference tree, provide further evidence for an independent domestication event of chicken in subcontinent.