Genetic variation and relationships in the mitochondrial DNA D-loop region of Qinghai indigenous and commercial pig breeds
Cellular & Molecular Biology Letters volume 23, Article number: 31 (2018)
Bamei pigs are special germplasm resources in Qinghai. In the specific environment of the Qinghai plateau, local breeds, through long-term natural and artificial selection, have acquired a relatively stable population structure and genetic diversity. This study revealed Bamei pigs’ genetic diversity at the molecular level and the degree of introgression of commercial breeds into Bamei pigs.
The mitochondrial DNA D-loop region was amplified and sequenced using the ABI 3130 DNA sequencer. Sequences were collected, edited and aligned using the MEGA 5.0, DNASTAR and ClustalX 1.81 software. Haplotypes of all sequences and genetic diversity were analyzed by DnaSP 5.0 software.
Results and conclusions
We identified a total of 73 polymorphic sites in our 489 novel sequences, including 1 single variable site and 72 parsimony informative sites.
Genetic diversity was analyzed in this study. The results showed that haplotype diversity, nucleotide diversity and the average number of nucleotide differences of Bamei pigs were lower than those of commercial pigs. Synthetic evaluation of genetic diversity through principal component analysis indicated that Bamei pigs also showed low genetic diversity. In summary, the results of sequence analysis showed that Bamei pigs were low in genetic diversity, and the introgression of commercial pigs into Bamei pigs was serious.
With the high intensification of pig production, highly productive breeds became predominant in the pig production system, and the quantity of indigenous breeds was decreased, so breed resource protection become more important . The Bamei pig is the indigenous pig breed in Qinghai Province of China, where there is a plateau continental climate. Under the influence of the ecological environment and natural and artificial selection, the Bamei pig was developed, which has the characteristics of good adaptability, high resistance, large litter size, good maternal qualities, high fat deposition ability, good meat quality and flavor, genetic stability, strong inbreeding resistance, adaptation to cold climate conditions and extensive feeding and management. However, they have some shortcomings, for example, slow growth, low lactation quantity, and low fattening benefit . In the past 50 years, Qinghai Province has introduced some pig breeds such as Berkshire, white, Landrace, Duroc, Hampshire, England, Yorkshire and other foreign breeds, which were crossed with the Bamei pig through an economic hybrid and three-way hybrid. The offspring were fed with high nutritional diets. Although the production performance was increased, the Bamei pig breed is at risk of extinction .
mtDNA sequence analysis has enabled important investigation of the origin and diversification of animal populations [4,5,6,7,8]. mtDNA contains the displacement (D)-loop, containing regulatory sequences controlling both replication and transcription of mtDNA [1, 9]. In the present study, we examined sequence variation in mtDNA from Bamei pig populations and commercial pig populations. We also assessed the relative impact of commercial pigs on Bamei pig populations. This study aimed to acquire information on Bamei pig genetic diversity. The findings will be helpful for conservation and sustainable use of Bamei pig resources.
Sampling and sequencing
A total of 489 samples, including 489 individuals (4 pig breeds) distributed in Qinghai Province of China, were collected. (Information on collected samples is provided in Table 1.) Only ear tissues were collected into microcentrifuge tubes containing 75% ethanol, and preserved at − 80 °C. Animals were released immediately following treatment of the wounds with antiseptic.
Genomic DNA was extracted using the phenol-chloroform extraction method . A fragment of the D-loop region was amplified using the primers 5′-CCAAAAACAAAGCAGAGTGTAC-3′ and 5′-CGTTATGAGCTACCGTTATA-3′. The PCR reaction mixture consisted of 25 μL, containing 12.5 μL 2 × Eco Taq PCR Supermix containing 1 U Taq polymerase, 500 mΜ dNTPs, and 10×Taq buffer (Beijing TransGen Biotech Co., Ltd., China), 0.1 μg of template DNA, 0.4 μL of 10 pmol/mL of each primer and 11.6 μL of ddH2O. The cycling conditions were initial denaturation at 94 °C for 5 min, followed by 33 cycles of 94 °C for 30 s, 56 °C for 30 s and 72 °C for 30 s, and a final extension for 5 min at 72 °C. Amplified DNA fragments were purified following agarose gel electrophoresis and sequenced using the ABI 3130 DNA sequencer (Applied Biosystems, Foster City, CA, USA).
A total of 489 mtDNA D-loop sequences were obtained and analyzed in this study, including 115 sequences from indigenous pigs (Bamei pigs) and 374 sequences from commercial pigs (Duroc, Landrace and Yorkshire).
Original sequence data were obtained using the ABI PRISM DNA sequencer software. Sequences were edited using the DNASTAR software and aligned using ClustalX 1.81 . In all sequences the length of alignment was reduced to 435 bp, and these were used to perform additional analyses. In addition, the haplotypes of all sequences and genetic diversity were analyzed using DnaSP 5.0 software . Correlation analysis and principal component analysis (PCA) were investigated by SPSS 19.0.
Genetic diversity analysis
The 435 bp control region of mtDNA was used to analyze single nucleotide polymorphisms (SNPs) for all 489 sequences. No insertion/deletions (indels) were detected in our novel sequences. We identified a total of 73 polymorphic sites, including 1 single variable site and 72 parsimony informative sites. The four types of nucleotide mutations identified were transitions, transversions, insertions and deletions. The transition: transversion ratio R (Ts/Tv) was 13.29:1, indicating a strong transitional bias that is common in mammalian mitochondrial evolution .
The genetic diversity of Bamei pigs and commercial pigs was calculated (Table 2). Haplotype diversity (Hd) of pig populations was between 0.491 and 0.856. All commercial pigs had high haplotype diversities of over 0.7. Yorkshire had the highest haplotype diversity (0.856 ± 0.018). Bamei pigs had the lowest haplotype diversity (0.491 ± 0.055). The result showed that haplotype diversity of Bamei pigs was lower than that of commercial pigs. Nucleotide diversities (Pis) were in the range 0.00264–0.01559. Landrace had the highest nucleotide diversity (0.01564), while Bamei pigs had the lowest (0.00264). Average nucleotide diversity of commercial pigs was higher than that of Bamei pigs. For all varieties of pigs the average number of nucleotide differences (Ks) was between 6.815 and 1.15; for Landrace it was the highest (6.815), and for Bamei pigs the lowest (1.15) – lower than for all commercial pigs.
The correlation of Pi, Hd and K was analyzed using SPSS 19.0 (Table 3). The results showed that Hd, Pi and K were positively correlated with each other, which indicated that all three indexes influenced the abundant degree of genetic diversity. Principal component analysis (PCA) was used for synthetic assessment of genetic diversity. The results are shown in Table 4. We obtained a synthesized assessment score (Fz). The Fz score indicated the highest genetic diversity in Landrace and the lowest for Bamei pigs.
Haplotype analysis of sequences
The 435 bp control region of mtDNA was used to calculate haplotypes for all 489 sequences. In total, 43 haplotypes were identified according to the distribution of variable sites (Table 5). Distribution frequencies of haplotypes indicated no equilibrium. The highest frequency haplotype was Hap34, which was shared by 97 sequences. The lowest frequency haplotypes were 26, and each haplotype harbored a sequence. Another 17 haplotypes were shared by two or more sequences. Hap31 and Hap34 were advantageous haplotypes which were present in more than 90 sequences. The Bamei porcine population shared 11 haplotypes. The commercial porcine population shared 36 haplotypes.
Introgression of commercial pigs into Qinghai indigenous pigs
Forty-three haplotypes in 489 individuals from indigenous and commercial breeds were identified. Eleven haplotypes were identified in indigenous pigs and 36 haplotypes were identified in commercial pigs. Seventeen shared haplotypes were identified and distributed in 96 indigenous and 123 commercial pigs. The ratio of the number of indigenous pigs with shared haplotypes and the total of indigenous pigs (Sc/S) showed the degree to which indigenous pigs were affected by commercial pigs. Sc/S was 83.48% (Table 6). Our data showed that Bamei breed was impacted by commercial breeds (Landrace, Yorkshire and Duroc).
Genetic diversity of Bamei and commercial porcine populations
Hd, Pi and K were the basic parameters which were used to assess genetic diversity. Hd is a measure of the uniqueness of a particular haplotype in a given population , which reflects haplotype abundance in a population. Pi and K measure the degree of intrapopulation haplotype mutation . Our study examined genetic diversity of Bamei and commercial porcine populations, and the results showed that the genetic diversities of both Bamei and commercial porcine populations were low, which is consistent with the trend of the genetic diversity of global livestock populations declining . Bamei pigs had the lowest haplotype diversity due to commercial hybridization. Bamei pigs have a low growth rate, so commercial pigs were used to cross with Bamei pigs to improve their performance. The Bamei population is, therefore, becoming smaller, and it is even becoming a critically endangered breed. Principal component analysis (PCA) is a statistical procedure to reduce the dimensionality of a data set by transformation to a new set of variables (the principal components) to summarize the features of the data . The diversity parameters, including Hd, Pi and K, were analyzed using PCA. Fz ranged from − 2.345 to 0.405, and its rank indicated that Bamei pigs had the lowest level of diversity, while Landrace had the highest diversity. This is consistent with the analysis of other studies which showed that indigenous pigs have a lower level of genetic diversity [17,18,19]. Genetic diversity is essential for continued breeding. This is especially true in the situation where future breeding goals differ from current goals . The Bamei pig as an indigenous porcine breed could be a useful resource for porcine production, so the protection of the Bamei pig from extinction will be important.
Introgression of commercial breeds into Qinghai indigenous breeds
Chinese indigenous breeds are known for their low growth rate and poor feed conversion efficiency, while commercial breeds are characterized by high meat yield, fast growth rate and excellent feed efficiency, so commercial breeds have been extensively applied in the swine industry . Increasing use of commercial lines threatens indigenous breeds and decreases genetic diversity. The Bamei pig similarly faces the threat of extinction due to the introgression of commercial pigs. Shared haplotypes were identified in 83.48% of Bamei pigs. Ninety-six haplotypes that were shared with commercial pigs were identified in 115 Bamei pigs. The reason for a high shared haplotype frequency in the Bamei population was crossbreeding with commercial lines. Current commercial lines were introduced and crossed with indigenous breeds, increasing the lineage of commercial lines. Despite its potentially economically unfavorable characteristics, the Bamei pig was recognized as an important genetic resource of the indigenous pig population due to specific traits (indigenous adaptation, strong adversity resistance, good maternal qualities). Additionally, the Bamei pig population was becoming small. The breeders are trying to maintain the Bamei pig breed and increase its number, so the founder effect might occur, which could result in a reduction in genetic diversity.
Our results showed that diversity parameters, such as haplotype diversity, nucleotide diversity and the average number of nucleotide differences, of Bamei pigs were lower than those of commercial pigs. Synthetic evaluation of genetic diversity through principal component analysis indicated that Bamei pigs also have low genetic diversity. The introgression of commercial pigs into Bamei pigs was serious due to the influence of commercial pigs.
Gao J, Ren J, Ken S, et al. The population genetic relationships among 47 Chinese and western pig populations inferred by AFLP markers. Chin J Anim Vet Sci. 2009;40(3):285–90.
The editorial committee of Chinese livestock and poultry resources. Chinese livestock and poultry resources. Beijing: Chinese Agriculture Press; 2004. p. 16.
Jin YC. Protection and utilization of Qinghai Bamei pig resource. Pig Prod Sci. 2006;11:70–4.
Kim KI, Lee JH, Li K, et al. Phylogenetic relationships of Asian and European pig breeds determined by mitochondrial DNA D-loop sequence polymorphism. Anim Genet. 2002;33:19–25.
Giuffra E, Kijas JM, Amarger V, Carlborg O, Jeon JT, et al. The origin of the domestic pig: independent domestication and subsequent introgression. Genetics. 2000;154:1785–91.
Okumura N, Kurosawa Y, Kobayashi E. Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs. Anim Genet. 2001;32:139–47.
Larson G, Dobney K, Albarella U, et al. Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science. 2005;307(5715):1618–21.
Ciesielski GL, Oliveira MT. Kaguni1 LS. Animal mitochondrial DNA replication. Enzymes. 2016;39:255–92.
Toro MA, Rodrigáñez JJ, Silió L. Rodr’ıguez MC. Genealogical analysis of a closed herd of black hairless Iberian pigs, cons. Biol. 2000;14:1843–51.
Sambrook J, Frietsch EF, Maniatis T. Molecular cloning: a laboratory manual. 3rd ed. Cold Spring Harbor, NY, USA: Cold Spring Harbor Laboratory; 2001.
Thompson JD, Gibson TJ, Plewniak F, et al. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876–82.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.
Kim KI, Yang YH, Lee SS. Phylogenetic relationships of Cheju horses to other breeds as determined by mtDNA D-loop sequence polymorphism. Anim Genet. 1989;30:102–8.
Nei M, Tajima F. DNA polymorphism detectable by restriction endonucleases. Genetics. 1981;97:145.
FAO. In: Rischkowsky B, Pilling D, editors. The state of the world’s animal genetic resources for food and agriculture. Rome: FAO; 2007. Available at: http://www.fao.org/docrep/010/a1250e/a1250e00.htm.
Jolliffe IT. Principal Component Analysis. New York: Springer; 1986.
Kim TH, Kim KS, Choi BH. Genetic structure of pig breeds from Korea and China using microsatellite loci analysis. J Anim Sci. 2005;83:2255–63.
Kim KS, Yeo JS, Kim JW. Assessment of genetic diversity of Korean native pig (Sus scrofa) using AFLP markers. Gen Genet Syst. 2002b;77:361–8.
Li KY, Chen C, Moran BF, Zhao S, Peng Z. Analysis of diversity and genetic relationships between four Chinese indigenous pig breeds and one Australian commercial pig breed. Anim Genet. 2000;31:322–5.
China national commission of animal genetic resources. Animal genetic resources in China pigs. Beijing: Chinese Agriculture Press; 2011.
Ai HS, Huang LS, Ren J. Genetic diversity, linkage disequilibrium and selection signatures in Chinese and western pigs revealed by genome-wide SNP markers. PLoS One. 2013;8(2):e56001.
We thank Junshan Xing for help with experiments. We are grateful to Wenzhi Li and Caibao for assistance with sample collections. We thank Victoria and one of the native English speaking editors for helping in proofreading and editing the English of the final manuscript.
JXZ and BCY conceived and designed the experiments. JXZ and GQS performed the experiments. JXZ and WXC analyzed the data. JXZ and WXC contributed reagents/materials/analysis tools. JXZ SGQ wrote the paper. All authors read and approved the final manuscript.
This work was supported by the Young Foundation of Qinghai University (No. 2014-QNY-5).
Ethics approval and consent to participate
Studies involving animals: The animal experiments were approved by the Animal Care and Use Committee of Qinghai University, Xining, China.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, J., Yang, B., Wen, X. et al. Genetic variation and relationships in the mitochondrial DNA D-loop region of Qinghai indigenous and commercial pig breeds. Cell Mol Biol Lett 23, 31 (2018). https://doi.org/10.1186/s11658-018-0097-x