Genetic diversity and population structure of six Chinese indigenous pig breeds in the Taihu Lake region revealed by sequencing data

The Chinese indigenous pig breeds in the Taihu Lake region are the most prolific pig breeds in the world. In this study, we investigated the genetic diversity and population structure of six breeds, including Meishan, Erhualian, Mi, Fengjing, Shawutou and Jiaxing Black, in this region using whole‐genome SNP data. A high SNP with proportions of polymorphic markers ranging from 0.925 to 0.995 was exhibited by the Chinese indigenous pigs in the Taihu Lake region. The allelic richness and expected heterozygosity also were calculated and indicated that the genetic diversity of the Meishan breed was the greatest, whereas that of the Fengjing breed was the lowest. The genetic differentiation, as indicated by the fixation index, exhibited an overall mean of 0.149. Both neighbor‐joining tree and principal components analysis were able to distinguish the breeds from each other, but structure analysis indicated that the Mi and Erhualian breeds exhibited similar major signals of admixture. With this genome‐wide comprehensive survey of the genetic diversity and population structure of the indigenous Chinese pigs in the Taihu Lake region, we confirmed the rationality of the current breed classification of the pigs in this region.     

Population structure and genetic diversity of four Henan pig populations

To investigate the population structure and genetic diversity of Henan indigenous pig breeds, samples from a total of 78 pigs of 11 breeds were collected, including four pig populations from Henan Province, three Western commercial breeds, three Chinese native pig breeds from other provinces and one Asian wild boar. The genotyping datasets were obtained by genotyping‐by‐sequencing technology. We found a high degree of polymorphism and rapid linkage disequilibrium decay in Henan pigs. A neighbor‐joining tree, principal component analysis and structure analysis revealed that the Huainan and Erhualian pigs were clustered together and that the Queshan black pigs were clearly grouped together but that the Nanyang and Yuxi pigs were extensively admixed with Western pigs. In addition, heterozygosity values might indicate that Henan indigenous pigs, especially the Queshan black and Huainan pigs, were subjected to little selection during domestication. The results presented here indicate that Henan pig breeds were admixed from Western breeds, especially Nanyang and Yuxi pigs. Therefore, establishment of purification and rejuvenation systems to implement conservation strategies is urgent. In addition, it is also necessary to accelerate genetic resources improvement and utilization using modern breeding technologies, such as genomic selection and genome‐wide association studies.     

Genome‐wide genetic variation discovery in Chinese Taihu pig breeds using next generation sequencing

The Chinese Taihu pig breeds are an invaluable component of the world's pig genetic resources, and they are the most prolific breeds of swine in the world. In this study, the genomes of 252 pigs of the six indigenous breeds in the Taihu Lake region were sequenced using the genotyping by genome reducing and sequencing approach. A total of 950 million good reads were obtained using an Illumina Hiseq2000 at an average depth of 13× (for SNP calling) and an average coverage of 2.3%. In total, 122 632 indels, 31 444 insertions, 44 056 deletions and 455 CNVs (copy number variants) were identified in the genomes of the pigs. Approximately 2.3% of these genetic markers were mapped to gene exon regions, and 25% were in QTL regions related to economically important traits. The KEGG pathway or GO enrichment analyses revealed that genetic variants assumed to be large‐effect mutations were significantly overrepresented in 22 SNP, 56 indel, 26 insertion, 28 deletion and three CNV gene sets. A total of 343 breed‐specific SNPs were also identified in the six Chinese indigenous pigs. The findings from this study can contribute to future investigations of the genetic diversity, population structure, positive selection signals and molecular evolutionary history of these pigs at the genome level and can serve as a valuable reference for improving the breeding and cultivation of these pigs.     

Genetic diversity,population structure and phylogenetic relationships of three indigenous pig breeds from Jiangxi Province,China, in a worldwide panel of pigs

Dongxiang Spotted, Pingxiang Two‐End‐Black and Yushan Black pigs are three indigenous breeds in Jiangxi Province, China, that have been listed in the national conservation program for Chinese indigenous livestock germplasm. Here, we investigated the genetic diversity and population structure of the nucleus populations of these three breeds in a worldwide context of European and Chinese pigs using Illumina Porcine 60K chips. Our data indicate that Dongxiang Spotted and Yushan Black pigs have recently experienced severe inbreeding. The two breeds show large runs‐of‐homozygosity values, long‐range extents of linkage disequilibrium and reduced observed heterozygosity. In contrast, Pingxiang Two‐End‐Black pigs are rich in genetic diversity and have few inbred individuals. Both phylogenetic and admixture analyses illustrate that Dongxiang Spotted and Yushan Black pigs are genetically close to their geographical neighbors in East China and that Pingxiang Two‐End‐Black pigs have a close relationship with three other Two‐End‐Black breeds in central China. Finally, we reconstructed the family structures of the three breeds and propose a reliable breeding strategy to better conserve these breeds.     

Genome‐wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.     

Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas of domestication

Microsatellite diversity in European and Chinese pigs was assessed using a pooled sampling method on 52 European and 46 Chinese pig populations. A Neighbor Joining analysis on genetic distances revealed that European breeds were grouped together and showed little evidence for geographic structure, although a southern European and English group could tentatively be assigned. Populations from international breeds formed breed specific clusters. The Chinese breeds formed a second major group, with the Sino-European synthetic Tia Meslan in-between the two large clusters. Within Chinese breeds, in contrast to the European pigs, a large degree of geographic structure was noted, in line with previous classification schemes for Chinese pigs that were based on morphology and geography. The Northern Chinese breeds were most similar to the European breeds. Although some overlap exists, Chinese breeds showed a higher average degree of heterozygosity and genetic distance compared to European ones. Between breed diversity was even more pronounced and was the highest in the Central Chinese pigs, reflecting the geographically central position in China. Comparing correlations between genetic distance and heterozygosity revealed that China and Europe represent different domestication or breed formation processes. A likely cause is a more diverse wild boar population in Asia, but various other possible contributing factors are discussed.     

Genetic relationships among Vietnamese local pigs investigated using genome‐wide SNP markers

Vietnam is one of the most important countries for pig domestication, and a total of 26 local breeds have been reported. In the present study, genetic relationships among the various pig breeds were investigated using 90 samples collected from local pigs (15 breeds) in 15 distantly separated, distinct areas of the country and six samples from Landrace pigs in Hanoi as an out‐group of a common Western breed. All samples were genotyped using the Illumina Porcine SNP60 v2 Genotyping BeadChip. We used 15 160–15 217 SNPs that showed a high degree of polymorphism in the Vietnamese breeds for identifying genetic relationships among the Vietnamese breeds. Principal components analysis showed that most pigs indigenous to Vietnam formed clusters correlated with their original geographic locations. Some Vietnamese breeds formed a cluster that was genetically related to the Western breed Landrace, suggesting the possibility of crossbreeding. These findings will be useful for the conservation and management of Vietnamese local pig breeds.     

Multilocus genotypic data reveal high genetic diversity and low population genetic structure of Iranian indigenous sheep

Iranian livestock diversity is still largely unexplored, in spite of the interest in the populations historically reared in this country located near the Fertile Crescent, a major livestock domestication centre. In this investigation, the genetic diversity and differentiation of 10 Iranian indigenous fat‐tailed sheep breeds were investigated using 18 microsatellite markers. Iranian breeds were found to host a high level of diversity. This conclusion is substantiated by the large number of alleles observed across loci (average 13.83, range 7–22) and by the high within‐breed expected heterozygosity (average 0.75, range 0.72–0.76). Iranian sheep have a low level of genetic differentiation, as indicated by the analysis of molecular variance, which allocated a very small proportion (1.67%) of total variation to the between‐population component, and by the small fixation index (F_ST = 0.02). Both Bayesian clustering and principal coordinates analysis revealed the absence of a detectable genetic structure. Also, no isolation by distance was observed through comparison of genetic and geographical distances. In spite of high within‐breed variation, signatures of inbreeding were detected by the F_IS indices, which were positive in all and statistically significant in three breeds. Possible factors explaining the patterns observed, such as considerable gene flow and inbreeding probably due to anthropogenic activities in the light of population management and conservation programmes, are discussed.     

Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel

Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47 401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair‐wise genetic differentiation (F_ST) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within‐breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.     

Genetic diversity present within the near-complete mtDNA genome of 17 breeds of indigenous Chinese pigs

The genetic diversity present within the near-complete mitochondrial genome (15,982 bp) was determined from 17 indigenous Chinese pig breeds and 3 European breeds. Animals were selected from 17 Chinese breeds that reflect the large phenotypic diversity of Chinese pigs and represent each of the six breed types, which are grouped based on morphological characteristics. Analysis of nucleotide diversity confirmed a high level of divergence between animals of European versus Asian origin; however, much more limited variation was observed between the 17 indigenous Chinese breeds. Each had a unique haplotype, but the lowest pairwise sequence divergence was only 0.01 +/- 0.01%, observed between the Tongcheng and Yushan Black. Comparison of control region sequence diversity revealed the 17 Chinese breeds contain a lower average pairwise distance (0.61 +/- 0.19%) than a group of European commercial breeds (0.91 +/- 0.21%). The dendrogram constructed from the near-complete mtDNA sequences showed the Chinese sequences loosely clustering into two groups. Although some correspondence with geographic origin was present, notable differences between the dendrogram and the traditional pig breed grouping system were observed.     

Genetic evidence for the introgression of Western NR6A1 haplotype into Chinese Licha breed associated with increased vertebral number

There is evidence that NR6A1 is a strong candidate for being a causal gene underlying vertebral number in pigs. The Licha Black is one of the leanest Chinese indigenous pig breeds, having an average vertebral number of 21.5. The introgression of Western germplasm into Licha Black, resulting in increased vertebral number, has been assumed but is not confirmed. This study detected allele frequencies of the NR6A1 causative mutation (c.575T>C) in 519 pigs from three Western and seven Chinese breeds including Licha Black, and evaluated the genetic variation in a 650-kb region containing NR6A1 in the 10 breeds. Allele T for increased vertebral number was fixed in Western breeds. In contrast, this allele was very rare in most of the Chinese native breeds. Notably, the T allele was present in the Licha Black at a rather higher frequency (0.585) and in the Laiwu at lower frequency (0.250). As expected, selection pressure has wiped out the genetic variability in the 650 kb region in Western breeds. Conversely, Chinese indigenous breeds showed a high degree of genetic variability in this region. However, the Licha Black displayed dramatically reduced heterozygosity at the loci proximal to the causative mutation. Moreover, a high proportion (45.9%) of Licha Black pigs and a small number (21%) of Laiwu pigs had the Western NR6A1 haplotype, and the two breeds showed closer relationships with Western commercial breeds than other Chinese breeds in the phylogenic tree. When the results are taken together, this study supports the assumption that the Western NR6A1 haplotypes were introduced into Licha Black and possibly Laiwu and are associated with increased vertebral number.     

Nucleotide variability and linkage disequilibrium patterns in the porcine MUC4 gene

ABSTRACT: BACKGROUND: MUC4 is a type of membrane anchored glycoprotein and serves as the major constituent of mucus that covers epithelial surfaces of many tissues such as trachea, colon and cervix. MUC4 plays important roles in the lubrication and protection of the surface epithelium, cell proliferation and differentiation, immune response, cell adhesion and cancer development. To gain insights into the evolution of the porcine MUC4 gene, we surveyed the nucleotide variability and linkage disequilibrium (LD) within this gene in Chinese indigenous breeds and Western commercial breeds. RESULTS: A total of 53 SNPs covering the MUC4 gene were genotyped on 5 wild boars and 307 domestic pigs representing 11 Chinese breeds and 3 Western breeds. The nucleotide variability, haplotype phylogeny and LD extent of MUC4 were analyzed in these breeds. Both Chinese and Western breeds had considerable nucleotide diversity at the MUC4 locus. Western pig breeds like Duroc and Large White have comparable nucleotide diversity as many of Chinese breeds, thus artificial selection for lean pork production have not reduced the genetic variability of MUC4 in Western commercial breeds. Haplotype phylogeny analyses indicated that MUC4 had evolved divergently in Chinese and Western pigs. The dendrogram of genetic differentiation between breeds generally reflected demographic history and geographical distribution of these breeds. LD patterns were unexpectedly similar between Chinese and Western breeds, in which LD usually extended less than 20 kb. This is different from the presumed high LD extent (more than 100 kb) in Western commercial breeds. The significant positive Tajima'D, and Fu and Li's D statistics in a few Chinese and Western breeds implied that MUC4 might undergo balancing selection in domestic breeds. Nevertheless, we cautioned that the significant statistics could be upward biased by SNP ascertainment process. CONCLUSIONS: Chinese and Western breeds have similar nucleotide diversity but evolve divergently in the MUC4 region. Western breeds exhibited unusual low LD extent at the MUC4 locus, reflecting the complexity of nucleotide variability of pig genome. The finding suggests that high density (e.g. 1SNP/10 kb) markers are required to capture the underlying causal variants at such regions.     

The chimerical genome of Isla del Coco feral pigs (Costa Rica), an isolated population since 1793 but with remarkable levels of diversity

The history of domestic species and of their wild ancestors is not a simple one, and feral processes can clarify key aspects of this history, including the adaptive processes triggered by new environments. Here, we provide a comprehensive genomic study of Isla del Coco (Costa Rica) feral pigs, a unique population that was allegedly founded by two individuals and has remained isolated since 1793. Using SNP arrays and genome sequencing, we show that Cocos pigs are hybrids between Asian and European pigs, as are modern international pig breeds. This conclusively shows that, as early as the 18th century, British vessels were loading crossbred pigs in Great Britain and transporting them overseas. We find that the Y chromosome has Asian origin, which has not been reported in any international pig breed. Chinese haplotypes seem to have been transmitted independently between Cocos and other pig breeds, suggesting independent introgression events and a complex pattern of admixing. Although data are compatible with a founder population of N = 2, variability levels are as high in Cocos pigs as in international pig breeds (~1.9 SNPs/kb) and higher than in European wild boars or local breeds (~1.7 SNPs/kb). Nevertheless, we also report a 10‐Mb region with a marked decrease in variability across all samples that contains four genes (CPE, H3F3C, SC4MOL and KHL2) previously identified as highly differentiated between wild and domestic pigs. This work therefore illustrates how feral population genomic studies can help to resolve the history of domestic species and associated admixture events.     

Genetic diversity and population structure of 10 Chinese indigenous egg-type duck breeds assessed by microsatellite polymorphism

The genetic structure and diversity of 10 Chinese indigenous egg-type duck breeds were investigated using 29 microsatellite markers. The total number of animals examined were 569, on average 57 animals per breed were selected. The microsatellite marker set analysed provided 177 alleles (mean 6.1 alleles per locus, ranging from 3 to 10). All populations showed high levels of heterozygosity with the lowest estimate of 0.539 for the Jinding ducks, and the highest 0.609 observed for Jingjiang partridge ducks. The global heterozygote deficit across all populations (F _IT) amounted to −0.363. About 10% of the total genetic variability originated from differences among breeds, with all loci contributing significantly. An unrooted consensus tree was constructed using the NeighborNet tree based on the Reynold’s genetic distance. The structure software was used to assess genetic clustering of these egg-type duck breeds. Clustering analysis provided an accurate representation of the current genetic relations among the breeds. An integrated analysis was undertaken to obtain information on the population dynamics in Chinese indigenous egg-type duck breeds, and to better determine the conservation priorities.     

Whole‐genome SNP data unravel population structure and signatures of selection for black plumage of indigenous chicken breeds from Jiangxi province,China

Ten indigenous chicken breeds were originally distributed in Jiangxi Province, China, and they define a critical component of Chinese chicken genetic resources. We have investigated the population genetics of seven Jiangxi chicken breeds using 600K chicken BeadChip SNP data. To provide a genome‐wide perspective for the population structure of all 10 Jiangxi chicken breeds, we herein genotyped 78 additional individuals from the seven breeds and 63 chickens from three uninvestigated breeds—Yugan Black (YG), Nancheng Black (NC) and Wanzai Yellow using 55K chicken SNP arrays. We then explored merged data of 17 101 SNPs from 235 individuals to infer the population structure of the 10 breeds. We showed that NC and YG are two regional populations of the same breed, as individuals from the two populations clustered together to form a branch separate from the other breeds in the neighbor‐joining tree, they always grouped together in multidimensional principal component analyses and they displayed an identical pattern of ancestral lineage composition. Hence, NC and YG should be considered a single breed in the state‐supported conservation scheme. Moreover, we conducted a genome scan for signatures of selection for black plumage. bayescan and hapflk analyses of two contrasting groups (three black‐feathered breeds vs. six non‐black‐feathered breeds) consistently detected 25 putative regions under selection. Nine pigmentation‐ associated genes (DCT, SLC24A5, SLC30A4, MYO5A, CYP19A1, NADK2, SLC45A2, GNAQ and DCP2) reside within these regions, and these genes are interesting candidates for black plumage and provide a starting point for further identification of causative mutations for black feathers in chicken.     

Genetic structure,relationships and admixture with wild relatives in native pig breeds from Iberia and its islands

Background

Native pig breeds in the Iberian Peninsula are broadly classified as belonging to either the Celtic or the Mediterranean breed groups, but there are other local populations that do not fit into any of these groups. Most of the native pig breeds in Iberia are in danger of extinction, and the assessment of their genetic diversity and population structure, relationships and possible admixture between breeds, and the appraisal of conservation alternatives are crucial to adopt appropriate management strategies.

Methods

A panel of 24 microsatellite markers was used to genotype 844 animals representing the 17 most important native swine breeds and wild populations existing in Portugal and Spain and various statistical tools were applied to analyze the results.

Results

Genetic diversity was high in the breeds studied, with an overall mean of 13.6 alleles per locus and an average expected heterozygosity of 0.80. Signs of genetic bottlenecks were observed in breeds with a small census size, and population substructure was present in some of the breeds with larger census sizes. Variability among breeds accounted for about 20% of the total genetic diversity, and was explained mostly by differences among the Celtic, Mediterranean and Basque breed groups, rather than by differences between domestic and wild pigs. Breeds clustered closely according to group, and proximity was detected between wild pigs and the Mediterranean cluster of breeds. Most breeds had their own structure and identity, with very little evidence of admixture, except for the Retinto and Entrepelado varieties of the Mediterranean group, which are very similar. Genetic influence of the identified breed clusters extends beyond the specific geographical areas across borders throughout the Iberian Peninsula, with a very sharp transition from one breed group to another. Analysis of conservation priorities confirms that the ranking of a breed for conservation depends on the emphasis placed on its contribution to the between- and within-breed components of genetic diversity.

Conclusions

Native pig breeds in Iberia reveal high levels of genetic diversity, a solid breed structure and a clear organization in well-defined clusters.     

Genomic diversity and population structure of Mexican and Spanish bovine Lidia breed

The Lidia bovine breed is distinguished for its low genetic exchangeability given its selection on aggressive behavior, its management uniqueness and its subdivided structure. In this study, we present a comprehensive genome‐wide analysis of genetic diversity, population structure and admixture of 468 animals from Mexican and Spanish Lidia breed populations and 64 samples belonging to 10 Spanish native and American‐creole breeds using 37 148 single nucleotide polymorphisms. We found similar average inbreeding values in the Lidia breed, with different distributions within groups; variability of inbreeding values among Spanish lineages was significant and no differences were found among the Mexican sub‐populations. Together, the high F_IS of the lineages and the behavior of the runs of homozygosity are consequences of the lineage's small effective population sizes, contributing to their inbreeding increase. Population admixture analysis discarded any influence on the genetic structure of the Lidia populations from the Spanish native and American‐creole breeds. In addition, both Lidia populations depicted different genetic origins, with the exception of some Mexican individuals whose origins traced back to recent Spanish importations.     

Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation

After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual‐based approaches (principal component analysis and population structure) and population‐based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.