Revealing New Candidate Genes for Teat Number Relevant Traits in Duroc Pigs Using Genome-Wide Association Studies

Genome-Wide Association Study and Phenotype Prediction of Reproductive Traits in Large White Pigs

In a study involving 385 Large White pigs, a genome-wide association study (GWAS) was conducted to investigate reproductive traits, specifically the number of healthy litters (NHs) and the number of weaned litters (NWs). Several SNP loci, including ALGA0098819, ALGA0037969, and H3GA0032302, were significantly associated with these traits. In the combined-parity analysis, candidate genes, such as BLVRA, STK17A, PSMA2, and C7orf25, were identified. GO and KEGG pathway enrichment analyses revealed that these genes are involved in key biological processes, including organic synthesis, the regulation of sperm activity, spermatogenesis, and meiosis. In the by-parity analysis, the PLCXD3 gene was significantly associated with the NW trait in the second and fourth parities, while RNASEH1, PYM1, and SEPTIN9 were linked to cell proliferation, DNA repair, and metabolism, suggesting their potential role in regulating reproductive traits. These findings provide new molecular markers for the genetic study of reproductive traits in Large White pigs. For the phenotypic prediction of NH and NW traits, several machine learning models (GBDT, RF, LightGBM, and Adaboost.R2), as well as traditional models (GBLUP, BRR, and BL), were evaluated using SNP data in varying proportions. After PCA processing, the GBDT model achieved the highest PCC for NH (0.141), while LightGBM reached the highest PCC for NW (0.146). The MAE, MSE, and RMSE results showed that the traditional models exhibited stable error rates, while the machine learning models performed comparatively better across the different SNP ratios. Overall, PCA processing provided some improvement in the predictive performance of all of the models, though the overall increase in accuracy was limited.

Identification of Goat Supernumerary Teat Phenotype Using Wide-Genomic Copy Number Variants

Supernumerary teats (SNTs) or nipples often emerge around the mammary line. This study performed a genome-wide selective sweep analysis (GWS) at the copy number variant (CNV) level using two selected signal calculation methods (VST and FST) to identify candidate genes associated with SNTs in goats. A total of 12,310 CNVs were identified from 37 animals and 123 CNVs, with the top 1% VST values including 84 candidate genes (CDGs). Of these CDGs, minichromosome maintenance complex component 3, ectodysplasin A receptor associated via death domain, and cullin 5 demonstrated functions closely related to mammary gland development. In addition, 123 CNVs with the top 1% FST values were annotated to 97 CDGs. 5-Hydroxytryptamine receptor 2A, CCAAT/enhancer-binding protein alpha, and the polymeric immunoglobulin receptor affect colostrum secretion through multiple signaling pathways. Two genes, namely, RNA-binding motif protein 46 and β-1,3-galactosyltransferase 5, showed a close relation to mammary gland development. Six CNVs were identified and annotated to five genes by intersecting the top 1% of candidate CNVs with both parameters. These genes include LOC102185621, LOC102190481, and UDP-glucose pyrophosphorylase 2, which potentially affect the occurrence of BC through multiple biological processes, such as cell detoxification, glycogen synthesis, and phospholipid metabolism. In conclusion, we discovered numerous genes related to mammary development and breast cancer (BC) through a GWS, which suggests the mechanism of SNTs in goats and a certain association between mammary cancer and SNTs.

Genome-wide detection of multiple variants associated with teat number in French Yorkshire pigs

BACKGROUND

Teat number is a vital reproductive trait in sows, crucial for providing immunity and nutrition to piglets during lactation. However, "missing heritability" in Single Nucleotide Polymorphism (SNP)-based Genome-Wide Association Studies (GWAS) has led to an increasing focus on structural variations in the genetic analysis of complex biological traits.

RESULTS

In this study, we generated a comprehensive CNV map in a population of French Yorkshire pigs (n = 644) and identified 429 CNVRs. Notably, 44% (189 CNVRs) of these were detected for the first time. Subsequently, we conducted GWAS for teat number in the French Yorkshire pig population using both 80K chip and its imputed data, as well as a GWAS analysis based on CNV regions (CNVR). Interestingly, 80K chip GWAS identified two SNPs located on Sus scrofa chromosome 5 (SSC5) that were simultaneously associated with Total Teat Number (TTN), Left Teat Number (LTN), and Right Teat Number (RTN). The leading SNP (WU_10.2_5_76130558) explained 3.33%, 2.69%, and 2.67% of the phenotypic variance for TTN, LTN, and RTN, respectively. Moreover, through imputed GWAS, we successfully identified 30 genetic variants associated with TTN located within the 73.22 -73.30 Mb region on SSC5. The two SNPs identified in the 80K chip GWAS were also located in this region. In addition, CNVR-based GWAS revealed three significant CNVRs associated with TTN. Finally, through gene annotation, we pinpointed two candidate genes, TRIM66 and PRICKLE1, which are related to diverse processes such as breast cancer and abnormal vertebral development.

CONCLUSIONS

Our research provides an in-depth analysis of the complex genetic structure underlying teat number, contributing to the genetic enhancement of sows with improved reproductive performance and, ultimately, bolstering the economic benefits of swine production enterprises.

Weighted single-step genome-wide association study to reveal new candidate genes for productive traits of Landrace pig in Korea

The objective of this study was to identify genomic regions and candidate genes associated with productive traits using a total of 37,099 productive records and 6,683 single nucleotide polymorphism (SNP) data obtained from five Great-Grand-Parents (GGP) farms in Landrace. The estimated of heritabilities for days to 105 kg (AGE), average daily gain (ADG), backfat thickness (BF), and eye muscle area (EMA) were 0.49, 0.49, 0.56, and 0.23, respectively. We identified a genetic window that explained 2.05%-2.34% for each trait of the total genetic variance. We observed a clear partitioning of the four traits into two groups, and the most significant genomic region for AGE and ADG were located on the Sus scrofa chromosome (SSC) 1, while BF and EMA were located on SSC 2. We conducted Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), which revealed results in three biological processes, four cellular component, three molecular function, and six KEGG pathway. Significant SNPs can be used as markers for quantitative trait loci (QTL) investigation and genomic selection (GS) for productive traits in Landrace pig.

Enhancement of selectivity, 25-hydroxyvitamin D3 level, alkaline phosphatase activity and reproductive performance in gilts and primiparous sows using 14-epimer of 25-hydroxyvitamin D3

Selecting breed-worthy gilts as sow replacements is essential for continuity of pig production cycle. Though vitamin D3 (VD3) is known to enhance reproductive performance of multiparous sows, there is still a knowledge gap on its impact in developing gilts and primiparous sows. This study was aimed to quantify plasma 25-hydroxyvitamin D3 (25(OH)D3), serum alkaline phosphatase (ALP), and examine the reproductive performance of primiparous sows fed diets supplemented with regular VD3, and its 25(OH)D3 epimers. The study sample comprised 10-week-old replacement gilts (50 % Landrace x 50 % Yorkshire, N = 180) assigned in a randomized complete block design to three treatments [2,000 IU/kg of VD3 (T1), 25 µg/kg of 14‑epi-25(OH)D3, half dose (T2), and 50 µg/kg of 25(OH)D3 (T3)] equilibrated to 2,000 IU/kg in base diets. Selections occurred at 22, 27 and 35 weeks of age, respectively. Plasma 25(OH)D3, serum alkaline phosphatase (ALP), bone structure and reproductive performance were analyzed. Dietary treatments influenced carpus (P = 0.023), fore view stance (P = 0.017), infantile vulva (P = 0.014), inverted (P = 0.048), and prominent teat (P < 0.001). Post-partum 25(OH)D3 concentration and ALP activity were elevated by day 25 (P < 0.001). Treatment diets also influenced total born (P < 0.001), born alive (P = 0.048), and still born (P = 0.049). Two factors affect circulating 25(OH)D3 and ALP activity: physiological changes in sows during lactation, and dietary 25(OH)D3 intake. 14‑epi-25(OH)D3 is a potent metabolite for improving maturation of reproductive organs in developing gilts. It also reduces still birth in primiparous sows.

Genome-wide association analysis revealed new QTL and candidate genes affecting the teat number in Dutch Large White pigs

Teat number (TNUM) is an important reproductive trait of sows, which affects the weaning survival rate of piglets. In this study, 1166 Dutch Large White pigs with TNUM phenotype were used as the research object. These pigs were genotyped by 50K SNP chip and the chip data were further imputed to the resequencing level. The estimated heritabilities of left teat number (LTN), right teat number (RTN) and total teat number (TTN) were 0.21, 0.19 and 0.3, respectively. Based on chip data, significant SNPs for RTN on SSC2, SSC5, SSC9 and SSC13 were identified using genome-wide association analysis (GWAS). Significant SNPs for TTN were identified on SSC2, SSC5 and SSC7. Based on imputed data, the GWAS identified a significant SNP (rs329158522) for LTN on SSC17, two significant SNPs (rs342855242 and rs80813115) for RTN on SSC2 and SSC9, and two significant SNPs (rs327003548 and rs326943811) for TTN on SSC5 and SSC6. Among them, four novel QTL were discovered. The Bayesian fine-mapping method was used to fine map the QTL identified in the GWAS of the imputed data, and the confidence intervals of QTL affecting LTN (SSC17: 45.22-46.20 Mb), RTN (SSC9: 122.18-122.80 Mb) and TTN (SSC5: 14.01-15.91 Mb, SSC6: 120.06-121.25 Mb) were detected. A total of 52 candidate genes were obtained. Furthermore, we identified five candidate genes, WNT10B, AQP5, FMNL3, NUAK1 and CKAP4, for the first time, which involved in breast development and other related functions by gene annotation. Overall, this study provides new molecular markers for the breeding of teat number in pigs.

Genome-wide association study to reveal new candidate genes using single-step approaches for productive traits of Yorkshire pig in Korea

OBJECTIVE

The objective is to identify genomic regions and candidate genes associated with age to 105 kg (AGE), average daily gain (ADG), backfat thickness (BF), and eye muscle area (EMA) in Yorkshire pig.

METHODS

This study used a total of 104,380 records and 11,854 single nucleotide polymorphism (SNP) data obtained from Illumina porcine 60K chip. The estimated genomic breeding values (GEBVs) and SNP effects were estimated by single-step genomic best linear unbiased prediction (ssGBLUP).

RESULTS

The heritabilities of AGE, ADG, BF, and EMA were 0.50, 0.49, 0.49, and 0.23, respectively. We identified significant SNP markers surpassing the Bonferroni correction threshold (1.68×10-6), with a total of 9 markers associated with both AGE and ADG, and 4 markers associated with BF and EMA. Genome-wide association study (GWAS) analyses revealed notable chromosomal regions linked to AGE and ADG on Sus scrofa chromosome (SSC) 1, 6, 8, and 16; BF on SSC 2, 5, and 8; and EMA on SSC 1. Additionally, we observed strong linkage disequilibrium on SSC 1. Finally, we performed enrichment analyses using gene ontology and Kyoto encyclopedia of genes and genomes (KEGG), which revealed significant enrichments in eight biological processes, one cellular component, one molecular function, and one KEGG pathway.

CONCLUSION

The identified SNP markers for productive traits are expected to provide valuable information for genetic improvement as an understanding of their expression.

Genome-Wide Association Study Meta-Analysis Elucidates Genetic Structure and Identifies Candidate Genes of Teat Number Traits in Pigs

The teat number is a pivotal reproductive trait that significantly influences the survival rate of piglets. A meta-analysis is a robust instrument, enhancing the universality of research findings and improving statistical power by increasing the sample size. This study aimed to identify universal candidate genes associated with teat number traits using a genome-wide association study (GWAS) meta-analysis with three breeds. We identified 21 chromosome threshold significant single-nucleotide polymorphisms (SNPs) associated with five teat number traits in single-breed and cross-breed meta-GWAS analyses. Using a co-localization analysis of expression quantitative trait loci and GWAS loci, we detected four unique genes that were co-localized with cross-breed GWAS loci associated with teat number traits. Through a meta-analysis and integrative analysis, we identified more reliable candidate genes associated with multiple-breed teat number traits. Our research provides new information for exploring the genetic mechanism affecting pig teat number for breeding selection and improvement.

Genome-wide association study using a single-step approach for teat number in Duroc, Landrace and Yorkshire pigs in Korea

We investigated the genetic basis of teat number in sows, which is an important factor in their reproductive performance. We collected genotyping data from 20 353 pigs of three breeds (Duroc, Landrace and Yorkshire) using the Porcine SNP60K Bead Chip, and analyzed phenotypic data from 240 603 pigs. The heritability values of total teat number were 0.33 ± 0.02, 0.51 ± 0.01 and 0.50 ± 0.01 in Duroc, Landrace and Yorkshire pigs, respectively. A genome-wide association study was used to identify significant chromosomal regions associated with teat number in SSC7 and SSC9 in Duroc pig, SSC3, SSC7 and SSC18 in Landrace pig, and SSC7, SSC8 and SSC10 in Yorkshire pig. Among the markers, MARC0038565, located between the vertnin (VRTN) and synapse differentiation-inducing 1-like (SYNDIG1L) genes, showed the strongest association in the Duroc pig and was significant in all breeds. In Landrace and Yorkshire pigs, the most significant markers were located within the apoptosis resistant E3 ubiquitin protein ligase 1 (AREL1) and latent transforming growth factor beta-binding protein 2 (LTBP2) genes in SSC7, respectively. VRTN is a candidate gene regulating the teat number. Most markers were located in SSC7, indicating their significance in determining teat number and their potential as valuable genomic selection targets for improving this trait. Extensive linkage disequilibrium blocks were identified in SSC7, supporting their use in genomic selection strategies. Our study provides valuable insights into the genetic architecture of teat numbers in pigs, and helps identify candidate genes and genomic regions that may contribute to this economically important trait.

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 of Parasite Resistance to Gastrointestinal Nematodes in Corriedale Sheep

Selection of genetically resistant animals is one alternative to reduce the negative impact of gastrointestinal nematodes (GIN) on sheep production. The aim of this study was to identify genomic regions associated with GIN resistance in Corriedale sheep by single-step genome-wide association studies (ssGWAS) using 170, 507 and 50K single nucleotide polymorphisms (SNPs). Analysis included 19,547 lambs with faecal egg counts (FEC) records, a pedigree file of 40,056 animals and 454, 711 and 383 genotypes from 170, 507 and 50K SNPs, respectively. Genomic estimated breeding values (GEBV) were obtained with single-step genomic BLUP methodology (ssGBLUP), using a univariate animal model, which included contemporary group, type of birth and age of dam as class fixed effects and age at FEC recording as covariate. The SNP effects as wells as p-values were estimated with POSTGSF90 program. Significance level was defined by a chromosome-wise False Discovery Rate of 5%. Significant genomic regions were identified in chromosomes 1, 3, 12 and 19 with the 170 SNP set, in chromosomes 7, 12 and 24 using the 507 SNP chip and only in chromosome 7 with the 50K SNP chip. Candidate genes located in these regions, using Oar_v4.0 as reference genome, were TIMP3, TLR5, LEPR and TLR9 (170 SNPs), SYNDIG1L and MGRN1 (507 SNP chip) and INO80, TLN2, TSHR and EIF2AK4 (50K SNP chip). These results validate genomic regions associated with FEC previously identified in Corriedale and other breeds and report new candidate regions for further investigation.

Genome-Wide Association Studies, Runs of Homozygosity Analysis, and Copy Number Variation Detection to Identify Reproduction-Related Genes in Bama Xiang Pigs

Litter size and teat number are economically important traits in the porcine industry. However, the genetic mechanisms influencing these traits remain unknown. In this study, we analyzed the genetic basis of litter size and teat number in Bama Xiang pigs and evaluated the genomic inbreeding coefficients of this breed. We conducted a genome-wide association study to identify runs of homozygosity (ROH), and copy number variation (CNV) using the novel Illumina PorcineSNP50 BeadChip array in Bama Xiang pigs and annotated the related genes in significant single nucleotide polymorphisms and common copy number variation region (CCNVR). We calculated the ROH-based genomic inbreeding coefficients (FROH) and the Spearman coefficient between FROH and reproduction traits. We completed a mixed linear model association analysis to identify the effect of high-frequency copy number variation (HCNVR; over 5%) on Bama Xiang pig reproductive traits using TASSEL software. Across eight chromosomes, we identified 29 significant single nucleotide polymorphisms, and 12 genes were considered important candidates for litter-size traits based on their vital roles in sperm structure, spermatogenesis, sperm function, ovarian or follicular function, and male/female infertility. We identified 9,322 ROHs; the litter-size traits had a significant negative correlation to FROH. A total of 3,317 CNVs, 24 CCNVR, and 50 HCNVR were identified using cnvPartition and PennCNV. Eleven genes related to reproduction were identified in CCNVRs, including seven genes related to the testis and sperm function in CCNVR1 (chr1 from 311585283 to 315307620). Two candidate genes (NEURL1 and SH3PXD2A) related to reproduction traits were identified in HCNVR34. The result suggests that these genes may improve the litter size of Bama Xiang by marker-assisted selection. However, attention should be paid to deter inbreeding in Bama Xiang pigs to conserve their genetic diversity.

A Whole Genome Sequencing-Based Genome-Wide Association Study Reveals the Potential Associations of Teat Number in Qingping Pigs

Teat number plays an important role in the reproductive performance of sows and the growth of piglets. However, the quantitative trait loci (QTLs) and candidate genes for the teat number-related traits in Qingping pigs remain unknown. In this study, we performed GWAS based on whole-genome single-nucleotide polymorphisms (SNPs) and insertions/deletions (Indels) for the total number of teats and five other related traits in 100 Qingping pigs. SNPs and Indels of all 100 pigs were genotyped using 10× whole genome resequencing. GWAS using General Linear Models (GLM) detected a total of 28 SNPs and 45 Indels as peak markers for these six traits. We also performed GWAS for the absolute difference between left and right teat number (ADIFF) using Fixed and random model Circulating Probability Unification (FarmCPU). The most strongly associated SNP and Indel with a distance of 562,788 bp were significantly associated with ADIFF in both GLM and FarmCPU models. In the 1-Mb regions of the most strongly associated SNP and Indel, there were five annotated genes, including TRIML1, TRIML2, ZFP42, FAT1 and MTNR1A. We also highlighted TBX3 as an interesting candidate gene for SSC14. Enrichment analysis of candidate genes suggested the Wnt signaling pathway may contribute to teat number-related traits. This study expanded significant marker-trait associations for teat number and provided useful molecular markers and candidate genes for teat number improvement in the breeding of sows.