Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Whole-genome sequencing revealed genetic diversity, structure and patterns of selection in Guizhou indigenous chickens

BACKGROUND

The eight phenotypically distinguishable indigenous chicken breeds in Guizhou province of China are great resources for high-quality development of the poultry industry in China. However, their full value and potential have yet to be understood in depth. To illustrate the genetic diversity, the relationship and population structure, and the genetic variation patterns shaped by selection in Guizhou indigenous chickens, we performed a genome-wide analysis of 240 chickens from 8 phenotypically and geographically representative Guizhou chicken breeds and 60 chickens from 2 commercial chicken breeds (one broiler and one layer), together with 10 red jungle fowls (RJF) genomes available from previous studies.

RESULTS

The results obtained in this present study showed that Guizhou chicken breed populations harbored higher genetic diversity as compared to commercial chicken breeds, however unequal polymorphisms were present within Guizhou indigenous chicken breeds. The results from the population structure analysis markedly reflected the breeding history and the geographical distribution of Guizhou indigenous chickens, whereas, some breeds with complex genetic structure were ungrouped into one cluster. In addition, we confirmed mutual introgression within Guizhou indigenous chicken breeds and from commercial chicken breeds. Furthermore, selective sweep analysis revealed candidate genes which were associated with specific and common phenotypic characteristics evolved rapidly after domestication of Guizhou local chicken breeds and economic traits such as egg production performance, growth performance, and body size.

CONCLUSION

Taken together, the results obtained from the comprehensive analysis of the genetic diversity, genetic relationships and population structures in this study showed that Guizhou indigenous chicken breeds harbor great potential for commercial utilization, however effective conservation measures are currently needed. Additionally, the present study drew a genome-wide selection signature draft for eight Guizhou indigenous chicken breeds and two commercial breeds, as well as established a resource that can be exploited in chicken breeding programs to manipulate the genes associated with desired phenotypes. Therefore, this study will provide an essential genetic basis for further research, conservation, and breeding of Guizhou indigenous chickens.

Genetic dissection and genomic prediction for pork cuts and carcass morphology traits in pig

BACKGROUND

As pre-cut and pre-packaged chilled meat becomes increasingly popular, integrating the carcass-cutting process into the pig industry chain has become a trend. Identifying quantitative trait loci (QTLs) of pork cuts would facilitate the selection of pigs with a higher overall value. However, previous studies solely focused on evaluating the phenotypic and genetic parameters of pork cuts, neglecting the investigation of QTLs influencing these traits. This study involved 17 pork cuts and 12 morphology traits from 2,012 pigs across four populations genotyped using CC1 PorcineSNP50 BeadChips. Our aim was to identify QTLs and evaluate the accuracy of genomic estimated breed values (GEBVs) for pork cuts.

RESULTS

We identified 14 QTLs and 112 QTLs for 17 pork cuts by GWAS using haplotype and imputation genotypes, respectively. Specifically, we found that HMGA1, VRTN and BMP2 were associated with body length and weight. Subsequent analysis revealed that HMGA1 primarily affects the size of fore leg bones, VRTN primarily affects the number of vertebrates, and BMP2 primarily affects the length of vertebrae and the size of hind leg bones. The prediction accuracy was defined as the correlation between the adjusted phenotype and GEBVs in the validation population, divided by the square root of the trait's heritability. The prediction accuracy of GEBVs for pork cuts varied from 0.342 to 0.693. Notably, ribs, boneless picnic shoulder, tenderloin, hind leg bones, and scapula bones exhibited prediction accuracies exceeding 0.600. Employing better models, increasing marker density through genotype imputation, and pre-selecting markers significantly improved the prediction accuracy of GEBVs.

CONCLUSIONS

We performed the first study to dissect the genetic mechanism of pork cuts and identified a large number of significant QTLs and potential candidate genes. These findings carry significant implications for the breeding of pork cuts through marker-assisted and genomic selection. Additionally, we have constructed the first reference populations for genomic selection of pork cuts in pigs.

Genetic dissection of 26 meat cut, meat quality and carcass traits in four pig populations

BACKGROUND

Currently, meat cut traits are integrated in pig breeding objectives to gain extra profit. However, little is known about the heritability of meat cut proportions (MCP) and their correlations with other traits. The aims of this study were to assess the heritability and genetic correlation of MCP with carcass and meat quality traits using single nucleotide polymorphism chips and conduct a genome-wide association study (GWAS) to identify candidate genes for MCP.

RESULTS

Seventeen MCP, 12 carcass, and seven meat quality traits were measured in 2012 pigs from four populations (Landrace; Yorkshire; Landrace and Yorkshire hybrid pigs; Duroc, and Landrace and Yorkshire hybrid pigs). Estimates of the heritability for MCP ranged from 0.10 to 0.55, with most estimates being moderate to high and highly consistent across populations. In the combined population, the heritability estimates for the proportions of scapula bone, loin, back fat, leg bones, and boneless picnic shoulder were 0.44 ± 0.04, 0.36 ± 0.04, 0.44 ± 0.04, 0.38 ± 0.04, and 0.39 ± 0.04, respectively. Proportion of middle cuts was genetically significantly positively correlated with intramuscular fat content and backfat depth. Proportion of ribs was genetically positively correlated with carcass oblique length and straight length (0.35 ± 0.08 to 0.45 ± 0.07) and negatively correlated with backfat depth (- 0.26 ± 0.10 to - 0.45 ± 0.10). However, weak or nonsignificant genetic correlations were observed between most MCP, indicating their independence. Twenty-eight quantitative trait loci (QTL) for MCP were detected by GWAS, and 24 new candidate genes related to MCP were identified, which are involved with growth, height, and skeletal development. Most importantly, we found that the development of the bones in different parts of the body may be regulated by different genes, among which HMGA1 may be the strongest candidate gene affecting forelimb bone development. Moreover, as previously shown, VRTN is a causal gene affecting vertebra number, and BMP2 may be the strongest candidate gene affecting hindlimb bone development.

CONCLUSIONS

Our results indicate that breeding programs for MCP have the potential to enhance carcass composition by increasing the proportion of expensive cuts and decreasing the proportion of inexpensive cuts. Since MCP are post-slaughter traits, the QTL and candidate genes related to these traits can be used for marker-assisted and genomic selection.

Genome-Wide Re-Sequencing Data Reveals the Population Structure and Selection Signatures of Tunchang Pigs in China

Tunchang pig is one population of Hainan pig in the Hainan Province of China, with the characteristics of delicious meat, strong adaptability, and high resistance to diseases. To explore the genetic diversity and population structure of Tunchang pigs and uncover their germplasm characteristics, 10 unrelated Tunchang pigs were re-sequenced using the Illumina NovaSeq 150 bp paired-end platform with an average depth of 10×. Sequencing data from 36 individuals of 7 other pig breeds (including 4 local Chinese pig breeds (5 Jinhua, 5 Meishan, 5 Rongchang, and 6 Wuzhishan), and 3 commonly used commercial pig breeds (5 Duorc, 5 Landrace, and 5 Large White)) were downloaded from the NCBI public database. After analysis of genetic diversity and population structure, it has been found that compared to commercial pigs, Tunchang pigs have higher genetic diversity and are genetically close to native Chinese breeds. Three methods, FST, θπ, and XP-EHH, were used to detect selection signals for three breeds of pigs: Tunchang, Duroc, and Landrace. A total of 2117 significantly selected regions and 201 candidate genes were screened. Gene enrichment analysis showed that candidate genes were mainly associated with good adaptability, disease resistance, and lipid metabolism traits. Finally, further screening was conducted to identify potential candidate genes related to phenotypic traits, including meat quality (SELENOV, CBR4, TNNT1, TNNT3, VPS13A, PLD3, SRFBP1, and SSPN), immune regulation (CD48, FBL, PTPRH, GNA14, LOX, SLAMF6, CALCOCO1, IRGC, and ZNF667), growth and development (SYT5, PRX, PPP1R12C, and SMG9), reproduction (LGALS13 and EPG5), vision (SLC9A8 and KCNV2), energy metabolism (ATP5G2), cell migration (EPS8L1), and olfaction (GRK3). In summary, our research results provide a genomic overview of the genetic variation, genetic diversity, and population structure of the Tunchang pig population, which will be valuable for breeding and conservation of Tunchang pigs in the future.

Identification of Candidate Genes for Min Pig Villi Hair Traits by Genome-Wide Association of Copy Number Variation

The Min pig is a famous native pig breed in northeast China, which has the special genetic character of villi hair growth in cold seasons. At present, little research has focused on the genetic mechanism of villi hair growth in Min pigs. Copy number variations (CNVs) are a type of variant that may influence many traits. In this study, we first investigated the phenotype of Large White × Min pigs' F2 pig villi hair in detail and then performed a CNV-based genome-wide association study (GWAS) between CNVs and pig villi hair appearance. Finally, a total number of 15 significant CNVRs were found to be associated with Min pig villi hair. The most significant CNVR was located on chromosome 1. Nearest gene annotation analysis indicated that the pig villi hair traits may be associated with the biological process of the G-protein-coupled receptor signaling pathway. QTL overlapping analysis found that among the CNVRs, 14 CNVRs could be co-located with known QTLs. Some genes such as MCHR2, LTBP2, and GFRA2 may be candidate genes for pig villi traits and are worth further study. Our study may provide a basic reference for the selection and breeding of cold-resistant pigs and outdoor breeding.

Population Genetic Analysis of Six Chinese Indigenous Pig Meta-Populations Based on Geographically Isolated Regions

The diversification of indigenous pig breeds in China has resulted from multiple climate, topographic, and human cultural influences. The numerous indigenous pig breeds can be geographically divided into six meta-populations; however, their genetic relationships, contributions to genetic diversity, and genetic signatures remain unclear. Whole-genome SNP data for 613 indigenous pigs from the six Chinese meta-populations were obtained and analyzed. Population genetic analyses confirmed significant genetic differentiation and a moderate mixture among the Chinese indigenous pig meta-populations. The North China (NC) meta-population had the largest contribution to genetic and allelic diversity. Evidence from selective sweep signatures revealed that genes related to fat deposition and heat stress response (EPAS1, NFE2L2, VPS13A, SPRY1, PLA2G4A, and UBE3D) were potentially involved in adaptations to cold and heat. These findings from population genetic analyses provide a better understanding of indigenous pig characteristics in different environments and a theoretical basis for future work on the conservation and breeding of Chinese indigenous pigs.

A web tool for the global identification of pig breeds

BACKGROUND

Natural and artificial selection for more than 9000 years have led to a variety of domestic pig breeds. Accurate identification of pig breeds is important for breed conservation, sustainable breeding, pork traceability, and local resource registration.

RESULTS

We evaluated the performance of four selectors and six classifiers for breed identification using a wide range of pig breeds (N = 91). The internal cross-validation and external independent testing showed that partial least squares regression (PLSR) was the most effective selector and partial least squares-discriminant analysis (PLS-DA) was the most powerful classifier for breed identification among many breeds. Five-fold cross-validation indicated that using PLSR as the selector and PLS-DA as the classifier to discriminate 91 pig breeds yielded 98.4% accuracy with only 3K single nucleotide polymorphisms (SNPs). We also constructed a reference dataset with 124 pig breeds and used it to develop the web tool iDIGs ( http://alphaindex.zju.edu.cn/iDIGs_en/ ) as a comprehensive application for global pig breed identification. iDIGs allows users to (1) identify pig breeds without a reference population and (2) design small panels to discriminate several specific pig breeds.

CONCLUSIONS

In this study, we proved that breed identification among a wide range of pig breeds is feasible and we developed a web tool for such pig breed identification.

Population differentiated copy number variation between Eurasian wild boar and domesticated pig populations

Sus scrofa is a globally distributed livestock species that still maintains two different ways of life: wild and domesticated. Herein, we detected copy number variation (CNV) of 328 animals using short read alignment on Sscrofa11.1. We compared CNV among five groups of porcine populations: Asian domesticated (AD), European domesticated (ED), Asian wild (AW), European wild (EW), and Near Eastern wild (NEW). In total, 21,673 genes were identified on 154,872 copy number variation region (CNVR). Differences in gene copy numbers between populations were measured by considering the variance-based value [Formula: see text] and the one-way ANOVA test followed by Scheffe test. As a result, 111 genes were suggested as copy number variable genes. Abnormally gained copy number on EEA1 in all populations was suggested the presence of minor CNV in the reference genome assembly, Sscrofa11.1. Copy number variable genes were related to meat quality, immune response, and reproduction traits. Hierarchical clustering of all individuals and mean pairwise [Formula: see text] in breed level were visualized genetic relationship of 328 individuals and 56 populations separately. Our findings have shown how the complex history of pig evolution appears in genome-wide CNV of various populations with different regions and lifestyles.

Transcriptome Analysis Reveals Genes Contributed to Min Pig Villi Hair Follicle in Different Seasons

The Min pig, a local pig breed in China, has a special trait which has intermittent villus and coat hair regeneration. However, the regulation and mechanism of villus in Min pigs have not yet been described. We observed and described the phenotype of Min pig dermal villi in detail and sequenced the mRNA transcriptome of Min pig hair follicles. A total of 1520 differentially expressed genes (DEG) were obtained.K-means hierarchical clustering showed that there was a significant expression pattern difference in winter compared with summer. Gene enrichment and network analysis results showed that the hair growth in Min pigs was closely related to the composition of desmosomes and regulated by an interaction network composed of eight core genes, namely DSP, DSC3, DSG4, PKP1, TGM1, KRT4, KRT15, and KRT84. Methylation analysis of promoters of target genes showed that the PKP1 gene was demethylated. Our study will help to supplement current knowledge of the growth mechanism of different types of hair.