Next Article in Journal
The Study of Amorphous Kaempferol Dispersions Involving FT-IR Spectroscopy
Previous Article in Journal
Effective Antiviral Application of Antisense in Plants by Exploiting Accessible Sites in the Target RNA
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 24
10.3390/ijms242417154
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Genetic Diversity and Population Structure of Four Native Pig Populations in Gansu Province

by
Zunqiang Yan
1,
Kelin Song
1,
Pengfei Wang
1,
Shuangbao Gun
1,* and
Xi Long
2,*
1
College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730070, China
2
Chongqing Academy of Animal Sciences, Chongqing 402460, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(24), 17154; https://doi.org/10.3390/ijms242417154
Submission received: 31 October 2023 / Revised: 30 November 2023 / Accepted: 4 December 2023 / Published: 5 December 2023
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Indigenous pig populations, including Bamei pigs (BM), Hezuo pigs (HZ), Huixian Qingni Black pigs (HX), and Minxian Black pigs (MX) in Gansu Province, live in a particular climate and a relatively closed geographical environment. These local pig breeds are characterized by excellent characteristics (e.g., cold tolerance, robust disease resistance, and superior meat quality). In the past few years, pig populations in Gansu Province have decreased significantly because of their poor lean meat percentage, high fat content, and slow growth rate. Maintaining the diversity of these four breeds can act as a source of new alleles to be incorporated into commercial breeds which are more susceptible to disease and less adaptable to changing conditions because of inbreeding. Genomic data analysis is adequate for determining the genetic diversity and livestock breeding population structure, even in local pig populations. However, the genetic diversity and population structure of the four native pig populations in Gansu Province are still unknown. Thus, we used “Zhongxin-I” porcine chip for the SNP detection of 102 individuals living on four pig conservation farms. A total of 57,466 SNPs were identified among the four pig breeds. The linkage disequilibrium (LD) plot showed that MX had the highest level of LD, followed by BM, HZ, and HX. The observed heterozygosity (Ho) in all four populations was higher than the expected heterozygosity (He). A principal component analysis (PCA) demonstrated that the four local pig populations were isolated. The identity displayed by the state matrix and G matrix heat map results indicated that small numbers of individuals among the four pig breeds had a high genetic distance and weak genetic relationships. The results of the population genetic structure of BM, HZ, HX, and MX pigs showed a slight genetic diversity loss. Our findings enabled us to better understand the genome characteristics of these four indigenous pig populations, which will provide novel insights for the future germplasm conservation and utilization of these indigenous pig populations.

1. Introduction

Domesticated pigs (Sus scrofa), a species originating from wild boar, have played indispensable roles in productive human life [1,2]. Chinese farmers have contributed to forming various pig breeds (e.g., Taihu pig with a high litter size and Tibetan pig with strong adaptability) in different climate environments (e.g., tropical climates, cold climates, and low-oxygen climates). This long history of pig breeding in China has led the country to bear the largest number of pig breeders worldwide [3,4,5]. The data show that China has a rich variety of pig breeds, reaching 102 breeds (including 83 local breeds, 13 cultivated breeds, and six introduced breeds), with local breeds accounting for 81.37% of the total in China. Local pig breeds in China can be classified into six major types based on their geographical distribution, namely the North China type, South China type, Central China type, Southwest China type, Plateau type, and Jianghai type, which exhibit excellent characteristics (e.g., good meat quality, robust disease resistance, and strong vitality) [6,7,8]. Recently, market forces have been pushing out local breeds all over the world, because local pig breeds often have long growth cycles and poor lean meat percentages [6]. China is now a major pig-raising nation due to it having the highest production and inventory of pigs worldwide. Many pig farms are more likely to produce European pig breeds (e.g., Landrace, Yorkshire, Duroc, Pietran, and Hampshire) because of their high lean meat production, reasonable growth rate, and excellent feed conversion efficiency [9]. These commercial pigs have gradually become dominant in the Chinese pig industry [4]. As a result, many local pig breeds have been replaced or crossed with imported pig breeds, leading to the mixing of bloodlines or even the extinction of some valuable breeds [10,11]. Livestock genetic resources are national strategic resources related to seedstock industry revitalization. Due to their essential role in agriculture in China, pig genetic resources must be taken seriously [12]. Strengthening the conservation of pig genetic resources is essential, because conserving the diversity of indigenous pig breeds has cultural value and conserves candidate genes for disease resistance for commercial breeds’ genetic improvement. It is thus important for the sustainable development of Chinese animal husbandry [13,14].
Gansu Province, located in northwest China, has a unique climate characterized by cold temperatures and less rainfall. The special climatic conditions and the relatively closed geographical environment have generated four well adapted pig breeds: Bamei (BM), Hezuo (HZ), Huixian Qingni (HX), and Minxian (MX). These germplasm resources display many attractive traits (e.g., stress resistance, strong disease resistance, and high-quality meat). The Hezuo pig is also a plateau miniature pig commonly treated as a model organism for investigating plateau adaptability [15,16].
There is no denying that commercial pig breeds threaten indigenous pigs due to their excellent performance [1]. Commercial pigs have resulted in an apparent sharp decline in the indigenous pig breed population in Gansu Province. The Chinese government has taken many measures and efforts to conserve valuable indigenous pig breeds, such as establishing provincial conservation farms and performing national conservation programs. Additionally, the population structure of these four indigenous pig breeds in Gansu Province suffered considerable reductions due to African swine fevers outbreaks. However, the genetic structure and diversity of these four pig breeds based on genome level remain unknown. Thus, it is necessary to thoroughly investigate the genetic structure and diversity of existing pig breeds in Gansu Province. This contributes to managing effective conservation measures for farmers and maintaining the diversity of these pig breeds. In this paper, based on the data of the “Zhongxin-I” porcine chip, we selected BM pigs, HZ pigs, HX pigs, and MX pigs from conserved populations as the research objects to conduct a comprehensive analysis of the phylogenetic relationships and genetic diversity among four pig breeds in Gansu Province. Our findings will provide new insights into scientific conservation and management strategies for these valuable genetic resources in Gansu Province.

2. Results

2.1. Genetic Variation and SNP Sites Discovery

One hundred and two samples were initially collected for genotyping with the “Zhongxin-I” porcine chip. All samples passed the quality control filter. A total of 57,466 SNPs were detected across four pig populations. Among these SNPs, there were 529 markers with a minor allele frequency (MAF) of less than 0.01. There were 721 markers with an SNP detection rate of less than 0.90. A total of 4252 markers located on the X chromosome were found. In addition, the number of SNPs with insertions and deletions was six. After quality control, 49,648 markers were retained for subsequent analysis. The distribution of SNPs on the chromosomes was visualized based on the SNP distribution plot. The results show that SNPs were distributed on all chromosomes. The results indicate that chromosome 1 contained the most SNPs (Figure 1).

2.2. Genetic Diversity and Linkage Disequilibrium (LD) Analysis

Six statistics were calculated to evaluate the genetic diversity of the four pig breeds. We noticed that the MAF of the four indigenous pig breeds was 0.21–0.28. The MAF of the the HX population was the highest at 0.28, followed by the MX, BM, and HZ populations (Table 1, Figure 2A). Heterozygosity is often used to calculate the genetic diversity of populations. This paper used the observed heterozygosity (Ho) and expected heterozygosity (He) parameters to measure the genetic diversity of four pig populations. The Ho in all four populations was higher than He. The HZ pig population was found to have the lowest Ho and He values at 0.35 and 0.32, respectively. At the same time, the MX pig population showed the highest Ho value at 0.43. The HX pig population had the highest He value at 0.37 (Table 1, Figure 2B,C). Regarding nucleotide diversity (Pi) across four indigenous pig populations, the HX pig population displayed the highest Pi value at 0.39. In contrast, the HZ pig population had the lowest Pi value at 0.33 (Table 1, Figure 2D). We also evaluated the polymorphic marker ratio (Pn) and effective population content (Ne) based on the autosomal genome information of the HX, MX, BM, and HZ conserved pig populations. The HZ pig population had the lowest Pn value at 0.82. In contrast, the HX pig population had the highest Pn value at 0.92 (Table 1). The Ne values of the HX, MX, BM, and HZ pig populations were 2.2, 2.5, 1.4, and 1.7, respectively (Table 1).
Among the four populations, the LD plot showed that the MX pig population had the highest level of LD, followed by the BM, HZ, and HX pig populations. The LD coefficient (r2) decreased as the distance between loci increased in all four populations. The decay rate of LD was observed to be the fastest in the HX pig population, which is consistent with the LD plot results (Figure 2E). We used the FST value for statistical analysis to determine the genetic dissimilarity between each pair of four pig populations. The genetic differentiation among the HX, MX, BM, and HZ pig populations ranged from 0.12 to 0.19. The genetic differentiation in MX vs. HX and HZ vs. BM was 0.14 and 0.13, respectively. BM vs. HX had the lowest genetic differentiation value, at 0.12. The highest genetic differentiation value was 0.19 in the HZ vs. MX populations (Figure 2F).
A neighbor-joining (NJ) tree was created to investigate the genetic distances among the four indigenous pig breeds, providing a clear phylogenetic relationship among these pigs. The NJ tree for the four pig populations showed that the HX, MX, BM, and HZ pig populations have a monophyletic clade (Figure 3A). Principal component analysis (PCA) can identify the principal components representing the population structure based on individual genetic correlations. The PCA analysis displayed four local pig populations, indicating significant genetic differences between these groups. The first and second PCA explained 12.07% and 4.64% of the total variation, respectively (Figure 3B). We observed two to six subgroups (K-values) based on the research population and performed cross-validation on the clustering results. Based on the valley value of the cross-validation error rate, we determined that the optimal number of subgroups was K = 4 (Figure 3C). According to the population structure analysis, four populations were distinguished from each other when K was equal to 4 (Figure 3D).

2.3. Genetic Distance of the Conserved Population

Based on the identity by state (IBS) genetic distance matrix of the population structure, among the four populations, we found that small numbers of individuals had high genetic distances (Figure 4). Across the four indigenous pig populations, the HX pig population had the highest average genetic distance value at 0.29, ranging from 0.14 to 0.34 (Figure 4A). The MX pig population had the lowest average genetic distance value at 0.26, ranging from 0.14 to 0.32 (Figure 4B). We calculated the IBS distance as varying from 0.16 to 0.37 for the conserved population of BM pigs (Figure 4C). The genetic distance ranged from 0.15 to 0.33 in the HZ conserved pig population (Figure 4D).
We also calculated the genetic relationship based on the G matrix. The visualization results of the G matrix of the four conserved populations showed that most individuals have moderate genetic relationships (Figure 5). We found that the genetic relationship of the MX population was the weakest, with a value of 0.03 (Figure 5B). Additionally, the HZ pig population produced the highest value for its genetic relationship at 0.05 (Figure 5D). The genetic relationship values of the BM and HX pig populations were 0.04 and 0.04, respectively (Figure 5A,C).

3. Discussion

Domestic animal (including livestock and poultry) genetic diversity is closely related to human production, representing valuable genetic resources for cultivating new and excellent breeds of domestic animals and improving the core competitiveness of livestock products [17]. Human beings must protect domestic animal genetic diversity, because livestock genetic diversity can drive the high-quality development of animal husbandry, accelerate the progress of human civilization and promote economic development [18,19,20]. Based on data from the Domestic Animal Diversity Information System and Food and Agricultural Organization in 2004, one-third of the pig breeds in the world are cultivated in China [2,20,21]. However, some commercial breeds account for the highest proportion of genetics utilized for production and dominate the Chinese pig industry. Consequently, the number of Chinese indigenous breeds has decreased considerably over the past few decades [18].
The HX, MX, BM, and HZ pig breeds, which are essential Chinese genetic resources, are mainly distributed in the Gansu Province of China. These pigs are distinguished by their moderate intramuscular fat content, good flavor, tender meat quality, apparent cold tolerance and robust disease resistance. Western commercial pigs threaten to replace these pig breeds, leading to the erosion of valuable and particular genetically adapted resources in Gansu Province. Genetic diversity and population structure analyses are crucial for protecting these indigenous pig populations [21]. This study aimed to determine the population structure and genetic diversity to provide rational scientific conservation strategies for HX, MX, BM, and HZ pig populations. More and more studies have investigated the genetic diversity and population structure of various species of livestock (e.g., Min pig, Tongcheng pig, Tibetan sheep, Hainan cattle, and BoHuai goat) using whole-genome sequencing or SNP chips [14,22,23,24,25,26,27]. For instance, Meng et al. used a KPS porcine breeding chip v1 50K for SNP to study the genetic diversity of the Min pig-conserved population. They found that the genetic diversity of the conserved Min pig population was low, indicating a certain degree of inbreeding in the Min pig [28]. Yang et al. performed whole-genome sequencing to detect the population genetic structure and diversity of Beijing Black pigs. They demonstrated that the core population size was small in the conserved farm. However, existing breeding schemes effectively avoid inbreeding to a certain extent [11]. These studies indicated that using whole-genome sequencing or SNP chips to explore genetic diversity and population structure provides the necessary data to implement conservation programs.
This study used high-density SNP chip data to evaluate the conservation state of HX, MX, BM, and HZ pig populations for the first time. The expected heterozygosity (He) values of HX, MX, BM, and HZ conserved pig populations were lower than those in Jinhua pigs (He = 0.43) [29], Laiwu pigs (He = 0.38) [30], Erhualian pigs (He = 0.38) [5], and Meishan pigs (He = 0.38) [5]. The He values of HX, MX, BM, and HZ conserved pig populations were higher than those of Min pigs (He = 0.33) [28], Liangshan pigs (He = 0.35) [20], Wuzhishan pigs (He = 0.31) [31], and Licha Black pigs (He = 0.35) [12]. The moderate He of HX, MX, BM, and HZ conserved pig populations demonstrated that many SNP loci have been subjected to inbreeding. Additionally, the observed heterozygosity (Ho) values were higher than those of He in this study, indicating that they have been influenced by gene flow (such as the introduction of a foreign lineage). Unlike Western commercial pigs, indigenous pigs in China do not undergo intense artificial selection pressure [31]. As a consequence, Chinese indigenous pigs display abundant genetic diversity. The values of the Pi parameter of the four pig populations in this study were 0.38, 0.36, 0.37, and 0.33, respectively, which were higher than those of commercial pig breeds, such as Landrace pig (Pi = 0.0019) [31] and Large White pig (Pi = 0.0018) [31]. Effective population size (Ne) is an important parameter in pig genetic breeding. We speculated that some excellent commercial pig breeds have been introduced, leading to threats to indigenous pig breeds in Gansu Province. Although the conserved population of the four pig breeds were established, the Ne remained low because of the limited number of unique families available. It is a fact that many indigenous domestic animals in China are facing the same situation. Thus, it is vital to immediately minimize consanguinity in indigenous pig breed populations by paying more attention to breeding plans.
LD describes the nonrandom linear association of alleles at different loci generated by migration, selection, and genetic drift [32]. LD analysis can reveal the intensity of selection that breeds have been subjected to during their evolutionary process and provide information on the evolution of the population [33,34]. When selection drives LD, the slower the LD decay, the greater the intensity of selection, thereby inferring the differences in the selection intensity in various populations. This study found that the LD coefficients of the MX pig populations were highest among the four populations, followed by BM, HZ, and HX. These results suggest that the four pig populations posed weaker artificial selection intensity during domestication. Herein, we analyzed the genetic structure among four pig populations using a phylogenetic tree, PCA, and structure methods. Considering the significance of the pig population in the breeding process, the phylogenetic tree of pig populations was drawn using the NJ method. The phylogenetic tree has monophyletic clades for each breed, including BM, HZ, HX and MX. When K = 4 in the structure analysis, a clear division among the four pigs was observed, consistent with the FST value calculated between each other. These results suggest that the four local pig breeds displayed obvious genetic differences because of their unique geographic isolation.
Currently, most of the local pig breeds in China are conserved by closed breeding, which leads to the genomic kinship of the population being easily affected by selection. Therefore, understanding the kinship of the conserved populations is crucial for conserving the populations. The genomic kinship matrix established using a whole-genome SNP marker panel reflects the individual kinship when the marker density is suitable. In this study, IBS distance matrix and G matrix analyses showed that most individuals in the HX, MX, BM, and HZ pig populations were distinct, and few individuals were close to each other. Meanwhile, we calculated the genetic distance among the four pig populations as 0.29, 0.26, 0.29, and 0.26, respectively, which is approximately equal to the genetic distance of Liangshan pigs (0.28) [20].

4. Materials and Methods

4.1. Animal Care

All animal work was conducted following the guidelines on the care and use of experimental animals established by the Ministry of Agriculture in China. The Committee on Animal Ethics of Gansu Agricultural University approved this study under Permit No. 2006-398.

4.2. Animal and Sample Collection

One hundred and two ear tissue samples, including 20 BM pigs, 20 HZ pigs, 28 HX pigs, and 34 MX pigs, were randomly collected from conservation farms and core areas of origin (Table 2).

4.3. DNA Extraction

Genomic DNA was extracted from the ear tissues using a standard phenol/chloroform extraction method. DNA integrity and purity were detected using a Nanodrop 2000 nucleic acid protein analyzer (Thermo, Shanghai, China) and agarose gel electrophoresis (BIO RAD, Wuhan, China). The extracted DNA samples were stored at −20 °C for further analysis.

4.4. Single Nucleotide Polymorphism Genotyping

The qualified 102 DNA samples were genotyped with the “Zhongxin-I” porcine chip containing 51,315 SNPs (Beijing Compass Biotechnology Co., Ltd., Beijing, China), which provided the highest throughput array format based on the data. Nine populations of Chinese local pigs (e.g., Bamaxiang pig, Sutai pig, and Laiwu pig) provided reference genome information [28]. Therefore, the “Zhongxin-I” porcine chip was selected for this study.

4.5. Quality Control

Strict quality control criteria were implemented to enhance the accuracy of the data processing. The genotype data underwent quality evaluation using Plink (v1.90) software [35]. The panel eliminated unreliable genotypes that did not meet the following four criteria [28]: (i) individuals with an SNP call rate greater than 90%; (ii) SNP MAF greater than 0.01; (iii) Hardy–Weinberg equilibrium (HWE) test p-value less than 0.000001; and (iv) SNPs of autosomal loci were used.

4.6. Population Genetic Diversity Analysis

Plink (v1.90) software [35] was employed to compute the expected heterozygosity (He), observed heterozygosity (Ho), and polymorphic loci ratio (Pn). Nucleotide diversity (Pi) and the population differentiation index (FST) were calculated using VCFtools (v0.1.17) software [36]. Additionally, population effective size (Ne) was determined using SNeP (v1.1) software [37] to analyze the genetic diversity of the BM, HZ, HX, and MX pig populations in Gansu Province. Additionally, raincloud plots for the MAF, Ho, He, and Pi statistical indicators in the four pig populations were obtained using OmicShare tools, a free online platform for data analyses (http://www.omicshare.com/tools (accessed on 6 October 2023)).

4.7. Population Genetic Structure Analysis

PCA of the four pig populations was performed using GCTA (v1.94) software [38]. An evolutionary tree of the four pig populations was constructed using MEGA (v11) software [39], employing the NJ method. Admixture analysis was performed at various levels (K = 2 to 6) for the four populations using Admixture (v1.30) software [40].

4.8. LD Analysis

LD refers to the nonrandom combination of alleles at different population loci [41]. It indicates that the probability of the simultaneous occurrence of two alleles on the same chromosome is higher than the probability expected based on random distribution within the population [42,43]. This condition is called LD and is commonly represented using D’ and r2 values. This study analyzed the LD levels of four pig populations using PopLDdecay (v1.90) software [44].

4.9. Population Differentiation Index Analysis

FST, the population genetic differentiation index, measures population differentiation and genetic distance. Its values range from 0 to 1. A value of 0 indicates random mating and complete genetic similarity between two populations, while 1 implies complete isolation and no genetic similarity between two populations. FST analysis was performed using VCFtools (v0.1.17) software [36] in this study.

4.10. Relatedness Analysis

Plink (v1.90) software [35] was used to calculate the genetic distance based on identity by state (IBS) for the four pig populations, and an IBS matrix was constructed. This matrix was then used to analyze the genetic distance between individuals within the population. GCTA (v1.94) software [38] was utilized to build the G matrix, which was used to analyze relatedness among individuals within the population of four pig populations. The results of the analysis were visualized using R programming language.

5. Conclusions

This study analyzed the genetic diversity and structure of HX, MX, BM, and HZ pig populations at the molecular level based on 50K SNP chip data. We found that the genetic diversity and the effective content of the conserved HX and MX populations were high. These results demonstrate that conserved pig populations were unstratified purebred populations with rich genetic diversity, indicating that the HX and MX pig populations have been effectively conserved. However, the genetic diversity level of the BM and HZ pig populations was inadequate. Thus, we need to take some effective measures to avoid inbreeding to ensure the long-term preservation of genetic resources via establishing preservation measures and formulating a scientific mating plan. This study will provide valuable information for the conservation, protection, and utilization of Gansu-native pigs in China.

Author Contributions

Conceptualization, Z.Y.; methodology, P.W.; software, P.W. and K.S.; validation, Z.Y., P.W. and K.S.; formal analysis, K.S.; investigation, Z.Y.; resources, S.G.; data curation, P.W.; writing—original draft preparation, Z.Y.; writing—review and editing, S.G. and X.L.; visualization, P.W.; supervision, S.G.; project administration, Z.Y.; funding acquisition, S.G. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key Research and Development Program (2021YFD1200303) from Ministry of Science and Technology in China; Gansu New Pig Breeds Resource Mining Project (GSCQ-2022-01) from Animal Husbandry Technology Extension Station of Gansu Province; Protection and Quality Improvement of Gansu Local Pig Germplasm Resources Program (GSLK-2021-13) from Animal Husbandry Technology Extension Station of Gansu Province; Animal Husbandry Pig Industry Technology Innovation Team Project of Gansu Agricultural University (GAU-XKTD-2022-25) from Gansu Agricultural University; Artificial Insemination Technology and System in Small and Medium-Sized Pig Farms Program (GSLK-2022-10) from Animal Husbandry Technology Extension Station of Gansu Province.

Institutional Review Board Statement

All animal experiments were performed based on animal ethical standards and regulations approved by the Committee for Animal Ethics of the College of Animal Science and Technology, Gansu Agricultural University (approval number 2006–398).

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data used in this paper are publicly available and can be obtained upon reasonable request from the corresponding author.

Acknowledgments

We thank the College of Animal Science and Technology, Gansu Agricultural University, for allowing this research to be performed in the Key Laboratory of Animal Genetics, Breeding, and Reproduction. We thank Beijing Compass Biotechnology Co., Ltd. for providing technical support. Additionally, we thank Qiaoli Yang, Xiaoyu Huang, Xiaoli Gao, and Jiaojiao Yang, who helped to revise the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huang, M.; Yang, B.; Chen, H.; Zhang, H.; Wu, Z.; Ai, H.; Ren, J.; Huang, L. The fine-scale genetic structure and selection signals of Chinese indigenous pigs. Evol. Appl. 2020, 13, 458–475. [Google Scholar] [CrossRef]
  2. McGlone, J.J.; Archer, C.; Henderson, M. Interpretive review: Semiochemicals in domestic pigs and dogs. Front. Vet. Sci. 2022, 9, 967980. [Google Scholar] [CrossRef] [PubMed]
  3. Fang, Y.; Hao, X.; Xu, Z.; Sun, H.; Zhao, Q.; Cao, R.; Zhang, Z.; Ma, P.; Sun, Y.; Qi, Z.; et al. Genome-wide detection of runs of homozygosity in Laiwu pigs revealed by sequencing data. Front. Genet. 2021, 12, 629966. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, C.; Wang, H.Y.; Zhang, Y.; Tang, Z.L.; Li, K.; Liu, B. Genome-wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs. Mol. Ecol. Resour. 2015, 15, 414–424. [Google Scholar] [CrossRef]
  5. Wang, Z.; Chen, Q.; Yang, Y.; Liao, R.; Zhao, J.; Zhang, Z.; Chen, Z.; Zhang, X.; Xue, M.; Yang, H.; et al. Genetic diversity and population structure of six Chinese indigenous pig breeds in the Taihu lake region revealed by sequencing data. Anim. Genet. 2015, 46, 697–701. [Google Scholar] [CrossRef] [PubMed]
  6. Pugliese, C.; Sirtori, F. Quality of meat and meat products produced from Southern European pig breeds. Meat Sci. 2012, 90, 511–518. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, C.; Luo, J.Q.; Yu, B.; Zheng, P.; Huang, Z.Q.; Mao, X.B.; He, J.; Yu, J.; Chen, J.L.; Chen, D.W. Dietary resveratrol supplementation improves meat quality of finishing pigs through changing muscle fiber characteristics and antioxidative status. Meat Sci. 2015, 102, 15–21. [Google Scholar] [CrossRef] [PubMed]
  8. Zhao, X.; Mo, D.L.; Li, A.N.; Gong, W.; Xiao, S.Q.; Zhang, Y.; Qin, L.M.; Niu, Y.N.; Guo, Y.X.; Liu, X.H.; et al. Comparative analyses by sequencing of transcriptomes during skeletal muscle development between pig breeds differing in muscle growth rate and fatness. PLoS ONE 2011, 6, e19774. [Google Scholar]
  9. Yang, S.L.; Wang, Z.G.; Liu, B.; Zhang, G.X.; Zhao, S.H.; Yu, M.; Fan, B.; Li, M.H.; Xiong, T.A.; Li, K. Genetic variation and relationships of eighteen Chinese indigenous pig breeds. Genet. Sel. Evol. GSE 2003, 35, 657–671. [Google Scholar] [CrossRef]
  10. Xu, Z.; Sun, H.; Zhang, Z.; Zhao, Q.B.; Olasege, B.S.; Li, Q.M.; Yue, Y.; Ma, P.P.; Zhang, X.Z.; Wang, Q.S.; et al. Assessment of autozygosity derived from runs of homozygosity in Jinhua pigs disclosed by sequencing data. Front. Genet. 2019, 10, 274. [Google Scholar] [CrossRef] [PubMed]
  11. Yang, W.; Liu, Z.; Zhao, Q.; Du, H.; Yu, J.; Wang, H.; Liu, X.; Liu, H.; Jing, X.; Yang, H.; et al. Population genetic structure and selection signature analysis of Beijing black pig. Front. Genet. 2022, 13, 860669. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Y.; Dong, R.; Li, X.; Cui, C.; Yu, G. Analysis of the genetic diversity and family structure of the Licha black pig population on jiaodong peninsula, Shandong province, China. Animals 2022, 12, 1045. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, X.Y.; D’Alessandro, E.; Chi, C.L.; Moawad, A.S.; Zong, W.C.; Chen, C.; Song, C.Y. Genetic evaluation and population structure of jiangsu native pigs in China revealed by sine insertion polymorphisms. Animals 2022, 12, 1345. [Google Scholar] [CrossRef]
  14. Yuan, J.; Zhou, X.; Xu, G.; Xu, S.; Liu, B. Genetic diversity and population structure of Tongcheng pigs in china using whole-genome SNP chip. Front. Genet. 2022, 13, 910521. [Google Scholar] [CrossRef] [PubMed]
  15. Tang, Y.R.; Zhang, B.; Shi, H.X.; Yan, Z.Q.; Wang, P.F.; Yang, Q.L.; Huang, X.Y.; Li, J.Y.; Wang, Z.K.; Gun, S. Cloning, expression analysis and localization of DAZl gene implicated in germ cell development of male Hezuo pig. Anim. Biotechnol. 2023. [Google Scholar] [CrossRef]
  16. Yan, J.H.; Wang, P.F.; Yan, Z.Q.; Yang, Q.L.; Huang, X.Y.; Gao, X.L.; Li, J.Y.; Wang, Z.K.; Gao, Y.; Gun, S. Cloning of STC-1 and analysis of its differential expression in Hezuo pig. Anim. Biotechnol. 2023. [Google Scholar] [CrossRef] [PubMed]
  17. Cai, Y.; Quan, J.; Gao, C.; Ge, Q.; Jiao, T.; Guo, Y.; Zheng, W.; Zhao, S. Multiple domestication centers revealed by the geographical distribution of Chinese native pigs. Animals 2019, 9, 709. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, Z.; Xiao, Q.; Zhang, Q.Q.; Sun, H.; Chen, J.C.; Li, Z.C.; Xue, M.; Ma, P.P.; Yang, H.J.; Xu, N.Y.; et al. Genomic analysis reveals genes affecting distinct phenotypes among different Chinese and Western pig breeds. Sci. Rep. 2018, 8, 13352. [Google Scholar] [CrossRef] [PubMed]
  19. Cao, R.; Feng, J.; Xu, Y.J.; Fang, Y.F.; Zhao, W.; Zhang, Z.Y.; Zhang, Z.; Li, M.; Wang, Q.S.; Pan, Y.C. Genomic signatures reveal breeding effects of Lulai pigs. Genes 2022, 13, 1969. [Google Scholar] [CrossRef]
  20. Liu, B.; Shen, L.; Guo, Z.; Gan, M.; Chen, Y.; Yang, R.; Niu, L.; Jiang, D.; Zhong, Z.; Li, X.; et al. Single nucleotide polymorphism-based analysis of the genetic structure of Liangshan pig population. Anim. Biosci. 2021, 34, 1105–1115. [Google Scholar] [CrossRef] [PubMed]
  21. Halimani, T.E.; Mapiye, O.; Marandure, T.; Januarie, D.; Imbayarwo-Chikosi, V.E.; Dzama, K. Domestic free-range pig genetic resources in Southern africa: Progress and prospects. Diversity 2020, 12, 68. [Google Scholar] [CrossRef]
  22. Chen, S.; Chang, C.C.; Cui, K.; Yang, W.J.; Li, B.L.; Ni, S.H.; Zhang, W.C.; Li, S.Y.; Li, X.B.; Wu, G.S.; et al. Whole-genome analyses reveal the genomic diversity and selection signatures of Hainan cattle. Livest. Sci. 2023, 275, 105311. [Google Scholar] [CrossRef]
  23. Schmid, M.; Gurtler, J.; Schneider, H.; Bennewitz, J. Population structure and genomic analyses of carcass back length, shoulder width and leg width in purebred and crossbred german Merino sheep. Small Rumin. Res. 2023, 226, 107052. [Google Scholar] [CrossRef]
  24. Tian, D.H.; Han, B.Y.; Li, X.; Liu, D.H.; Zhou, B.C.; Zhao, C.C.; Zhang, N.; Wang, L.; Pei, Q.B.; Zhao, K. Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing. Anim. Biosci. 2023, 36, 991–1002. [Google Scholar] [CrossRef]
  25. Vijh, R.K.; Sharma, U.; Kapoor, P.; Raheja, M.; Arora, R.; Ahlawat, S.; Dureja, V. Design and validation of high-density SNP array of goats and population stratification of indian goat breeds. Gene 2023, 885, 147691. [Google Scholar] [CrossRef] [PubMed]
  26. Yao, Z.; Zhang, S.J.; Wang, X.W.; Guo, Y.W.; Xin, X.L.; Zhang, Z.J.; Xu, Z.J.; Wang, E.Y.; Jiang, Y.; Huang, Y.Z. Genetic diversity and signatures of selection in Bohuai goat revealed by whole-genome sequencing. BMC Genom. 2023, 24, 116. [Google Scholar] [CrossRef] [PubMed]
  27. Zhong, Z.Q.; Li, R.; Wang, Z.; Tian, S.S.; Xie, X.F.; Wang, Z.Y.; Na, W.; Wang, Q.S.; Pan, Y.; Xiao, Q. Genome-wide scans for selection signatures in indigenous pigs revealed candidate genes relating to heat tolerance. Animal 2023, 17, 100882. [Google Scholar] [CrossRef] [PubMed]
  28. Meng, F.B.; Cai, J.C.; Wang, C.A.; Fu, D.C.; Di, S.W.; Wang, X.B.; Chang, Y.; Xu, C.Z. Single nucleotide polymorphism-based analysis of the genetic structure of the Min pig conserved population. Anim. Biosci. 2022, 35, 1839–1849. [Google Scholar] [CrossRef] [PubMed]
  29. Li, S.J.; Yang, S.L.; Zhao, S.H.; Fan, B.; Yu, M.; Wang, H.S.; Li, M.H.; Liu, B.; Xiong, T.A.; Li, K. Genetic diversity analyses of 10 indigenous Chinese pig populations based on 20 microsatellites. J. Anim. Sci. 2004, 82, 368–374. [Google Scholar] [CrossRef]
  30. Wang, X.P.; Zhang, H.; Huang, M.; Tang, J.H.; Yang, L.J.; Yu, Z.Q.; Li, D.S.; Li, G.X.; Jiang, Y.C.; Sun, Y.X.; et al. Whole-genome SNP markers reveal conservation status, signatures of selection, and introgression in Chinese Laiwu pigs. Evol. Appl. 2021, 14, 383–398. [Google Scholar] [CrossRef] [PubMed]
  31. Zhong, Z.; Wang, Z.; Xie, X.; Tian, S.; Wang, F.; Wang, Q.; Ni, S.; Pan, Y.; Xiao, Q. Evaluation of the genetic diversity, population structure and selection signatures of three native Chinese pig populations. Animals 2023, 13, 2010. [Google Scholar] [CrossRef] [PubMed]
  32. Gibson, J.; Tapper, W.; Ennis, S.; Collins, A. Exome-based linkage disequilibrium maps of individual genes: Functional clustering and relationship to disease. Hum. Genet. 2013, 132, 233–243. [Google Scholar] [CrossRef] [PubMed]
  33. Hedrick, P.W. Assortative mating and linkage disequilibrium. G3-Genes Genomes Genet. 2017, 7, 55–62. [Google Scholar] [CrossRef] [PubMed]
  34. Smith, R.D. The nonlinear structure of linkage disequilibrium. Theor. Popul. Biol. 2020, 134, 160–170. [Google Scholar] [CrossRef] [PubMed]
  35. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. Plink: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [PubMed]
  36. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef]
  37. Barbato, M.; Orozco-TerWengel, P.; Tapio, M.; Bruford, M.W. SNeP: A tool to estimate trends in recent effective population size trajectories using genome-wide SNP data. Front. Genet. 2015, 6, 109. [Google Scholar] [CrossRef]
  38. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef]
  39. Tamura, K.; Stecher, G.; Kumar, S. MEGA 11 molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  40. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef]
  41. Qanbari, S. On the extent of linkage disequilibrium in the genome of farm animals. Front. Genet. 2020, 10, 1304. [Google Scholar] [CrossRef] [PubMed]
  42. Charles, B.A.; Shriner, D.; Rotimi, C.N. Accounting for linkage disequilibrium in association analysis of diverse populations. Genet. Epidemiol. 2014, 38, 265–273. [Google Scholar] [CrossRef] [PubMed]
  43. Clark, A.G.; Wang, X.; Matise, T. Contrasting methods of quantifying fine structure of human recombination. Annu. Rev. Genom. Hum. Genet. 2010, 11, 45–64. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, C.; Dong, S.S.; Xu, J.Y.; He, W.M.; Yang, T.L. Poplddecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 2019, 35, 1786–1788. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 17154 g001
Figure 1. Distribution of SNPs on chromosomes.
Figure 1. Distribution of SNPs on chromosomes.
Ijms 24 17154 g001
Ijms 24 17154 g002
Figure 2. Analysis of the genetic diversity of four pig populations. Orange, blue, red and green color represent HX, MX, BM and HZ in figure (A–D), respectively. (A) Raincloud plot for minor allele frequency (MAF). (B) Observed heterozygosity (Ho) displayed by a raincloud plot. (C) A raincloud plot exhibiting the expected heterozygosity (He). (D) A raincloud plot describing the change in nucleotide diversity (Pi). (E) Drawn linkage disequilibrium (LD) decay plot describing the change in the degree of LD among the four pig populations. (F) Heat map of FST distance for four pig populations.
Figure 2. Analysis of the genetic diversity of four pig populations. Orange, blue, red and green color represent HX, MX, BM and HZ in figure (A–D), respectively. (A) Raincloud plot for minor allele frequency (MAF). (B) Observed heterozygosity (Ho) displayed by a raincloud plot. (C) A raincloud plot exhibiting the expected heterozygosity (He). (D) A raincloud plot describing the change in nucleotide diversity (Pi). (E) Drawn linkage disequilibrium (LD) decay plot describing the change in the degree of LD among the four pig populations. (F) Heat map of FST distance for four pig populations.
Ijms 24 17154 g002
Ijms 24 17154 g003
Figure 3. Analysis of the genetic structure of four pig populations. (A) A neighbor-joining tree for all individuals. (B) Plot of the first and second principal components resulting from a principal component analysis of all pig populations. (C) The cross-validation error rate. (D) The plot of the ancestry compositions of four breeds uncovered with the assumed number of ancestries (K = 2–6). Each column represents an individual, in which the length of different colors fragments represents the proportion of an ancestor in the individual genome. A common color represents one ancestral cluster and the abscissa represents the group name.
Figure 3. Analysis of the genetic structure of four pig populations. (A) A neighbor-joining tree for all individuals. (B) Plot of the first and second principal components resulting from a principal component analysis of all pig populations. (C) The cross-validation error rate. (D) The plot of the ancestry compositions of four breeds uncovered with the assumed number of ancestries (K = 2–6). Each column represents an individual, in which the length of different colors fragments represents the proportion of an ancestor in the individual genome. A common color represents one ancestral cluster and the abscissa represents the group name.
Ijms 24 17154 g003
Ijms 24 17154 g004
Figure 4. Identity by state (IBS) distance matrix of HX (A), MX (B), BM (C), and HZ (D) pig-conserved populations. Each small square in the figure represents the genetic distance value between the two pairs from the first individual to the last individual. The smaller the value (the closer it is to green), the smaller the genetic distance between two individuals, meaning they were extremely similar, and vice versa.
Figure 4. Identity by state (IBS) distance matrix of HX (A), MX (B), BM (C), and HZ (D) pig-conserved populations. Each small square in the figure represents the genetic distance value between the two pairs from the first individual to the last individual. The smaller the value (the closer it is to green), the smaller the genetic distance between two individuals, meaning they were extremely similar, and vice versa.
Ijms 24 17154 g004
Ijms 24 17154 g005
Figure 5. The G matrix heat map of the conserved population of HX (A), MX (B), BM (C), and HZ (D) pigs. Each small square exhibits the value of the relationship between two pairs, from the first individual to the last individual. The larger the value (the closer it is to purple), the closer the relationship between two individuals.
Figure 5. The G matrix heat map of the conserved population of HX (A), MX (B), BM (C), and HZ (D) pigs. Each small square exhibits the value of the relationship between two pairs, from the first individual to the last individual. The larger the value (the closer it is to purple), the closer the relationship between two individuals.
Ijms 24 17154 g005
Table 1. Genetic diversity indicators of four indigenous pig breeds in Gansu Province.
Table 1. Genetic diversity indicators of four indigenous pig breeds in Gansu Province.
BreedMAFHoHePiPnNe
HX0.280.390.370.380.922.2
MX0.240.430.360.360.822.5
BM0.240.370.360.370.851.4
HZ0.210.350.320.330.821.7
Table 2. Basic information on samples from four domestic pig breeds in Gansu Province.
Table 2. Basic information on samples from four domestic pig breeds in Gansu Province.
BreedAbbreviationSample TypeSample SizeSampling Location
Bamei pigBMEar tissue20Lingtai, Pingliang, Gansu
Hezuo pigHZEar tissue20Hezuo, Gannan, Gansu
Huixian Qingni Black pigHXEar tissue28Huixian, Longnan, Gansu
Minxian Black pigMXEar tissue34Minxian, Dingxi, Gansu
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, Z.; Song, K.; Wang, P.; Gun, S.; Long, X. Evaluation of the Genetic Diversity and Population Structure of Four Native Pig Populations in Gansu Province. Int. J. Mol. Sci. 2023, 24, 17154. https://doi.org/10.3390/ijms242417154

AMA Style

Yan Z, Song K, Wang P, Gun S, Long X. Evaluation of the Genetic Diversity and Population Structure of Four Native Pig Populations in Gansu Province. International Journal of Molecular Sciences. 2023; 24(24):17154. https://doi.org/10.3390/ijms242417154

Chicago/Turabian Style

Yan, Zunqiang, Kelin Song, Pengfei Wang, Shuangbao Gun, and Xi Long. 2023. "Evaluation of the Genetic Diversity and Population Structure of Four Native Pig Populations in Gansu Province" International Journal of Molecular Sciences 24, no. 24: 17154. https://doi.org/10.3390/ijms242417154

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop