Genome-wide association study identifies markers associated with carcass and meat quality traits in Italian Large White pigs

Genetic Structure and Genome-Wide Association Analysis of Growth and Reproductive Traits in Fengjing Pigs

The Fengjing pig is one of the local pig breed resources in China and has many excellent germplasm characteristics. However, research on its genome is lacking. To explore the degree of genetic diversity of the Fengjing pig and to deeply explore its excellent traits, this study took Fengjing pigs as the research object and used the Beadchip Array Infinium iSelect-96|XT KPS_PorcineBreedingChipV2 for genotyping. We analyzed the genetic diversity, relatedness, inbreeding coefficient, and population structure within the Fengjing pig population. Our findings revealed that the proportion of polymorphic markers (PN) was 0.469, and the effective population size was 6.8. The observed and expected heterozygosity were 0.301 and 0.287, respectively. The G-matrix results indicated moderate relatedness within the population, with certain individuals exhibiting closer genetic relationships. The NJ evolutionary tree classified Fengjing boars into five family lines. The average inbreeding coefficient based on ROH was 0.318, indicating a high level of inbreeding. GWAS identified twenty SNPs significantly associated with growth traits (WW, 2W, and 4W) and reproductive traits (TNB and AWB). Notably, WNT8B, RAD21, and HAO1 emerged as candidate genes influencing 2W, 4W, and TNB, respectively. Genes such as WNT8B were verified by querying the PigBiobank database. In conclusion, this study provides a foundational reference for the conservation and utilization of Fengjing pig germplasm resources and offers insights for future molecular breeding efforts in Fengjing pigs.

Deep Learning-Based Automated Approach for Determination of Pig Carcass Traits

Pig carcass traits are among the most economically significant characteristics and are crucial for genetic selection in breeding and enhancing the economic efficiency. Standardized and automated carcass phenotyping can greatly enhance the measurement efficiency and accuracy, thereby facilitating the selection and breeding of superior pig carcasses. In this study, we utilized phenotypic images and data from 3912 pigs to propose a deep learning-based approach for the automated determination of pig carcass phenotypic traits. Using the YOLOv8 algorithm, our carcass length determination model achieves an average accuracy of 99% on the test set. Additionally, our backfat segmentation model, YOLOV8n-seg, demonstrates robust segmentation performance, with a Mean IoU of 89.10. An analysis of the data distribution comparing manual and model-derived measurements revealed that differences in the carcass straight length are primarily concentrated between -2 cm and 4 cm, while differences in the carcass diagonal length are concentrated between -3 cm and 2 cm. To validate the method, we compared model measurements with manually obtained data, achieving coefficients of determination (R2) of 0.9164 for the carcass straight length, 0.9325 for the carcass diagonal length, and 0.7137 for the backfat thickness, indicating high reliability. Our findings provide valuable insights into automating carcass phenotype determination and grading in pig production.

Genetic diversity and breed-informative SNPs identification in domestic pig populations using coding SNPs

Background: The use of breed-informative genetic markers, specifically coding Single Nucleotide Polymorphisms (SNPs), is crucial for breed traceability, authentication of meat and dairy products, and the preservation and improvement of pig breeds. By identifying breed informative markers, we aimed to gain insights into the genetic mechanisms that influence production traits, enabling informed decisions in animal management and promoting sustainable pig production to meet the growing demand for animal products. Methods: Our dataset consists of 300 coding SNPs genotyped from three Italian commercial pig populations: Landrace, Yorkshire, and Duroc. Firstly, we analyzed the genetic diversity among the populations. Then, we applied a discriminant analysis of principal components to identify the most informative SNPs for discriminating between these populations. Lastly, we conducted a functional enrichment analysis to identify the most enriched pathways related to the genetic variation observed in the pig populations. Results: The alpha diversity indexes revealed a high genetic diversity within the three breeds. The higher proportion of observed heterozygosity than expected revealed an excess of heterozygotes in the populations that was supported by negative values of the fixation index (FIS) and deviations from the Hardy-Weinberg equilibrium. The Euclidean distance, the pairwise FST, and the pairwise Nei's GST genetic distances revealed that Yorkshire and Landrace breeds are genetically the closest, with distance values of 2.242, 0.029, and 0.033, respectively. Conversely, Landrace and Duroc breeds showed the highest genetic divergence, with distance values of 2.815, 0.048, and 0.052, respectively. We identified 28 significant SNPs that are related to phenotypic traits and these SNPs were able to differentiate between the pig breeds with high accuracy. The Functional Enrichment Analysis of the informative SNPs highlighted biological functions related to DNA packaging, chromatin integrity, and the preparation of DNA into higher-order structures. Conclusion: Our study sheds light on the genetic underpinnings of phenotypic variation among three Italian pig breeds, offering potential insights into the mechanisms driving breed differentiation. By prioritizing breed-specific coding SNPs, our approach enables a more focused analysis of specific genomic regions relevant to the research question compared to analyzing the entire genome.

A Comprehensive Genomic Analysis of Chinese Indigenous Ningxiang Pigs: Genomic Breed Compositions, Runs of Homozygosity, and Beyond

Ningxiang pigs are a renowned indigenous pig breed in China, known for their meat quality, disease resistance, and environmental adaptability. In recent decades, consumer demand for meats from indigenous breeds has grown significantly, fueling the selection and crossbreeding of Ningxiang pigs (NXP). The latter has raised concerns about the conservation and sustainable use of Ningxiang pigs as an important genetic resource. To address these concerns, we conducted a comprehensive genomic study using 2242 geographically identified Ningxiang pigs. The estimated genomic breed composition (GBC) suggested 2077 pigs as purebred Ningxiang pigs based on a ≥94% NXP-GBC cut-off. The remaining 165 pigs were claimed to be crosses, including those between Duroc and Ningxiang pigs and between Ningxiang and Shaziling pigs, and non-Ningxiang pigs. Runs of homozygosity (ROH) were identified in the 2077 purebred Ningxiang pigs. The number and length of ROH varied between individuals, with an average of 32.14 ROH per animal and an average total length of 202.4 Mb per animal. Short ROH (1-5 Mb) was the most abundant, representing 66.5% of all ROH and 32.6% of total ROH coverage. The genomic inbreeding estimate was low (0.089) in purebred Ningxiang pigs compared to imported western pig breeds. Nine ROH islands were identified, pinpointing candidate genes and QTLs associated with economic traits of interest, such as reproduction, carcass and growth traits, lipid metabolism, and fat deposition. Further investigation of these ROH islands and candidate genes is anticipated to better understand the genomics of Ningxiang pigs.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Meta-analysis of genome-wide association studies uncovers shared candidate genes across breeds for pig fatness trait

Background

Average backfat thickness (BFT) is a critical complex trait in pig and an important indicator for fat deposition and lean rate. Usually, genome-wide association study (GWAS) was used to discover quantitative trait loci (QTLs) of BFT in a single population. However, the power of GWAS is limited by sample size in a single population. Alternatively, meta-analysis of GWAS (metaGWAS) is an attractive method to increase the statistical power by integrating data from multiple breeds and populations. The aim of this study is to identify shared genetic characterization of BFT across breeds in pigs via metaGWAS. 

Results

In this study, we performed metaGWAS on BFT using 15,353 pigs (5,143 Duroc, 7,275 Yorkshire, and 2,935 Landrace) from 19 populations. We detected 40 genome-wide significant SNPs (Bonferroni corrected P < 0.05) and defined five breed-shared QTLs in across-breed metaGWAS. Markers within the five QTL regions explained 7 ~ 9% additive genetic variance and showed strong heritability enrichment. Furthermore, by integrating information from multiple bioinformatics databases, we annotated 46 candidate genes located in the five QTLs. Among them, three important (MC4R, PPARD, and SLC27A1) and seven suggestive candidate genes (PHLPP1, NUDT3, ILRUN, RELCH, KCNQ5, ITPR3, and U3) were identified.

Conclusion

QTLs and candidate genes underlying BFT across breeds were identified via metaGWAS from multiple populations. Our findings contribute to the understanding of the genetic architecture of BFT and the regulating mechanism underlying fat deposition in pigs.

Weighted Single-Step Genomic Best Linear Unbiased Prediction Method Application for Assessing Pigs on Meat Productivity and Reproduction Traits

Changes in the accuracy of the genomic estimates obtained by the ssGBLUP and wssGBLUP methods were evaluated using different reference groups. The weighting procedure’s reasonableness of application Pwas considered to improve the accuracy of genomic predictions for meat, fattening and reproduction traits in pigs. Six reference groups were formed to assess the genomic data quantity impact on the accuracy of predicted values (groups of genotyped animals). The datasets included 62,927 records of meat and fattening productivity (fat thickness over 6–7 ribs (BF1, mm)), muscle depth (MD, mm) and precocity up to 100 kg (age, days) and 16,070 observations of reproductive qualities (the number of all born piglets (TNB) and the number of live-born piglets (NBA), according to the results of the first farrowing). The wssGBLUP method has an advantage over ssGBLUP in terms of estimation reliability. When using a small reference group, the difference in the accuracy of ssGBLUP over BLUP AM is from −1.9 to +7.3 percent points, while for wssGBLUP, the change in accuracy varies from +18.2 to +87.3 percent points. Furthermore, the superiority of the wssGBLUP is also maintained for the largest group of genotyped animals: from +4.7 to +15.9 percent points for ssGBLUP and from +21.1 to +90.5 percent points for wssGBLUP. However, for all analyzed traits, the number of markers explaining 5% of genetic variability varied from 71 to 108, and the number of such SNPs varied depending on the size of the reference group (79–88 for BF1, 72–81 for MD, 71–108 for age). The results of the genetic variation distribution have the greatest similarity between groups of about 1000 and about 1500 individuals. Thus, the size of the reference group of more than 1000 individuals gives more stable results for the estimation based on the wssGBLUP method, while using the reference group of 500 individuals can lead to distorted results of GEBV.

Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis

Backfat thickness (BFT) is complex and economically important traits in the pig industry, since it reflects fat deposition and can be used to measure the carcass lean meat percentage in pigs. In this study, all 6,550 pigs were genotyped using the Geneseek Porcine 50K SNP Chip to identify SNPs related to BFT and to search for candidate genes through genome-wide association analysis in two Duroc populations. In total, 80 SNPs, including 39 significant and 41 suggestive SNPs, and 6 QTLs were identified significantly associated with the BFT. In addition, 9 candidate genes, including a proven major gene MC4R, 3 important candidate genes (RYR1, HMGA1, and NUDT3) which were previously described as related to BFT, and 5 novel candidate genes (SIRT2, NKAIN2, AMH, SORCS1, and SORCS3) were found based on their potential functional roles in BFT. The functions of candidate genes and gene set enrichment analysis indicate that most important pathways are related to energy homeostasis and adipogenesis. Finally, our data suggest that most of the candidate genes can be directly used for genetic improvement through molecular markers, except that the MC4R gene has an antagonistic effect on growth rate and carcass lean meat percentage in breeding. Our results will advance our understanding of the complex genetic architecture of BFT traits and laid the foundation for additional genetic studies to increase carcass lean meat percentage of pig through marker-assisted selection and/or genomic selection.

Single-Step Genome Wide Association Study Identifies QTL Signals for Untrimmed and Trimmed Thigh Weight in Italian Crossbred Pigs for Dry-Cured Ham Production

Simple Summary

Along with the traditional traits, swine breeding programs for Italian dry-cured ham production have recently aimed to include novel phenotypes. The identification of the genomic regions underlying such new traits helps to untangle their genetic architecture and may provide useful information to be integrated in genetic selection. With this aim, we estimated genetic parameters and conducted a single step genome wide association studies (GWAS) on untrimmed and trimmed thigh weight considering two pig crossbred lines approved for Italian Protected Designation of Origin ham production. Quantitative trait loci (QTLs) were characterized based on the variance of 10-SNP sliding windows genomic estimated breeding values. In particular, we identified interesting QTL signals on several chromosomes, notably on chromosome 4, 6, 7 and 15. A high heritability and genetic correlation were observed for the two traits under investigation and although independent studies including other pig populations are required to disentangle the possible effects of specific linkage disequilibrium in our population, our findings suggest that such QTL could be investigated in future pig breeding programs to improve the reliability of genomic estimated breeding values for the dry-cured ham production.

Abstract

Protected Designation of Origin (PDO) dry-cured ham is the most important product in the Italian pig breeding industry, mainly oriented to produce heavy pig carcasses to obtain hams of the right weight and maturity. Recently, along with the traditional traits swine breeding programs have aimed to include novel carcass traits. The identification at the genome level of quantitative trait loci (QTLs) affecting such new traits helps to reveal their genetic determinism and may provide information to be integrated in prediction models in order to improve prediction accuracy as well as to identify candidate genes underlying such traits. This study aimed to estimate genetic parameters and perform a single step genome wide association studies (ssGWAS) on novel carcass traits such as untrimmed (UTW) and trimmed thigh weight (TTW) in two pig crossbred lines approved for the ham production of the Italian PDO. With this purpose, phenotypes were collected from ~1800 animals and 240 pigs were genotyped with Illumina PorcineSNP60 Beadchip. The single-step genomic BLUP procedure was used for the heritability estimation and to implement the ssGWAS. QTL were characterized based on the variance of 10-SNP sliding window genomic estimated breeding values. Moderate heritabilities were detected and QTL signals were identified on chromosome 1, 4, 6, 7, 11 and 15 for both traits. As expected, the genetic correlation among the two traits was very high (~0.99). The QTL regions encompassed a total of 249 unique candidate genes, some of which were already reported in association with growth, carcass or ham weight traits in pigs. Although independent studies are required to further verify our findings and disentangle the possible effects of specific linkage disequilibrium in our population, our results support the potential use of such new QTL information in future breeding programs to improve the reliability of genomic prediction.

How Does Protein Nutrition Affect the Epigenetic Changes in Pig? A Review

Simple Summary

Epigenetic mechanisms regulate gene expression and depend of nutrition. In farm animals, and concretely, in pigs, some papers on protein nutrition have been realized to improve several productive traits. Changes in protein diet influence on epigenetic mechanisms that could affect productive and reproductive traits in individuals and their offspring. The purpose of this review was to update the current knowledge about the effects of these nutritional changes on epigenetic mechanisms in pigs.

Abstract

Epigenetic changes regulate gene expression and depend of external factors, such as environment and nutrition. In pigs, several studies on protein nutrition have been performed to improve productive and reproductive traits. Indeed, these studies aimed not only to determine broad protein requirements but also pigs’ essential amino acids requirements. Moreover, recent studies tried to determine these nutritional requirements for each individual, which is known as protein precision nutrition. However, nutritional changes could affect different epigenetic mechanisms, modifying metabolic pathways both in a given individual and its offspring. Modifications in protein nutrition, such as change in the amino acid profile, increase or decrease in protein levels, or the addition of metabolites that condition protein requirements, could affect the regulation of some genes, such as myostatin, insulin growth factor, or genes controlling cholesterol and glucose metabolism pathways. This review summarizes the impact of most common protein nutritional strategies on epigenetic changes and describes their effects on regulation of gene expression in pigs. In a context where animal nutrition is shifting towards precision protein nutrition (PPN), further studies evaluating the effects of PPN on animal epigenetic are necessary.