Genome-wide association analyses identify known and novel loci for teat number in Duroc pigs using single-locus and multi-locus models

Effects of dietary yeast mannan-rich fraction supplementation on growth performance, intestinal morphology, and lymphoid tissue characteristics in weaned piglets challenged with Escherichia Coli F4

We evaluated the effects of supplementing yeast mannan-reach-fraction on growth performance, jejunal morphology and lymphoid tissue characteristics in weaned piglets challenged with E. Coli F4. A total of 20 crossbred piglets were used. At weaning, piglets were assigned at random to one of four groups: piglets challenged and fed the basal diet supplemented with yeast mannan-rich fraction (C-MRF, n = 5); piglets challenged and fed the basal diet (C-BD, n = 5); piglets not challenged and fed the basal diet supplemented with yeast mannan-rich fraction (NC-MRF, n = 5), and piglets not challenged and fed the basal diet (NC-BD). Each dietary treatment had five replicates. On days 4, 5 and 10, piglets were orally challenged with 108 CFU/mL of E. Coli F4. C-MRF piglets had higher BW (p = 0.002; interactive effect) than C-BD piglets. C-MRF piglets had higher (p = 0.02; interactive effect) ADG in comparison with C-BD piglets. C-MRF piglets had higher (p = 0.04; interactive effect) ADFI than C-BD piglets. The diameter of lymphoid follicles was larger (p = 0.010; interactive effect) in the tonsils of C-MRF piglets than C-BD piglets. Lymphoid cells proliferation was greater in the mesenteric lymphnodes and ileum (p = 0.04 and p = 0.03, respectively) of C-MRF piglets. A reduction (p > 0.05) in E. Coli adherence in the ileum of piglets fed MRF was observed. In conclusion, the results of the present study demonstrate that dietary yeast mannan-rich fraction supplementation was effective in protecting weaned piglets against E. Coli F4 challenge.

ABCD4 is associated with mammary gland development in mammals

BACKGROUND

Mammary gland development is a critical process in mammals, crucial for their reproductive success and offspring nourishment. However, the functional roles of key candidate genes associated with teat number, including ABCD4, VRTN, PROX2, and DLST, in this developmental process remain elusive. To address this gap in knowledge, we conducted an in-depth investigation into the dynamic expression patterns, functional implications, and regulatory networks of these candidate genes during mouse mammary gland development.

RESULTS

In this study, the spatial and temporal patterns of key genes were characterized in mammary gland development. Using time-series single-cell data, we uncovered differences in the expression of A bcd4, Vrtn, Prox2, and Dlst in cell population of the mammary gland during embryonic and adult stages, while Vrtn was not detected in any cells. We found that only overexpression and knockdown of Abcd4 could inhibit proliferation and promote apoptosis of HC11 mammary epithelial cells, whereas Prox2 and Dlst had no significant effect on these cells. Using RNA-seq and qPCR, further analysis revealed that Abcd4 can induce widespread changes in the expression levels of genes involved in mammary gland development, such as Igfbp3, Ccl5, Tlr2, and Prlr, which were primarily associated with the MAPK, JAK-STAT, and PI3K-AKT pathways by functional enrichment.

CONCLUSIONS

These findings revealed ABCD4 as a candidate gene pivotal for regulating mammary gland development and lactation during pregnancy by influencing PRLR expression.

Genome-wide association studies for loin muscle area, loin muscle depth and backfat thickness in DLY pigs

This study aimed at identifying genes associated with loin muscle area (LMA), loin muscle depth (LMD) and backfat thickness (BFT). We performed single-trait and multi-trait genome-wide association studies (GWASs) after genotyping 685 Duroc × (Landrace × Yorkshire) (DLY) pigs using the Geneseek Porcine 50K SNP chip. In the single-trait GWASs, we identified two, eight and two significant SNPs associated with LMA, LMD and BFT, respectively, and searched genes within the 1 Mb region near the significant SNPs with relevant functions as candidate genes. Consequently, we identified one (DOCK5), three (PID1, PITX2, ELOVL6) and three (CCR1, PARP14, CASR) promising candidate genes for LMA, LMD and BFT, respectively. Moreover, the multi-trait GWAS identified four significant SNPs associated with the three traits. In conclusion, the GWAS analysis of LMA, LMD and BFT in a DLY pig population identified several associated SNPs and candidate genes, further deepening our understanding of the genetic basis of these traits, and they may be useful for marker-assisted selection to improve the three traits in DLY pigs.

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 for Somatic Skeletal Traits in Duroc × (Landrace × Yorkshire) Pigs

The pig bone weight trait holds significant economic importance in southern China. To expedite the selection of the pig bone weight trait in pig breeding, we conducted molecular genetic research on these specific traits. These traits encompass the bone weight of the scapula (SW), front leg bone weight (including humerus and ulna) (FLBW), hind leg bone weight (including femur and tibia) (HLBW), and spine bone weight (SBW). Up until now, the genetic structure related to these traits has not been thoroughly explored, primarily due to challenges associated with obtaining the phenotype data. In this study, we utilized genome-wide association studies (GWAS) to discern single nucleotide polymorphisms (SNPs) and genes associated with four bone weight traits within a population comprising 571 Duroc × (Landrace × Yorkshire) hybrid pigs (DLY). In the analyses, we employed a mixed linear model, and for the correction of multiple tests, both the false discovery rate and Bonferroni methods were utilized. Following functional annotation, candidate genes were identified based on their proximity to the candidate sites and their association with the bone weight traits. This study represents the inaugural application of GWAS for the identification of SNPs associated with individual bone weight in DLY pigs. Our analysis unveiled 26 SNPs and identified 12 promising candidate genes (OPRM1, SLC44A5, WASHC4, NOPCHAP1, RHOT1, GLP1R, TGFB3, PLCB1, TLR4, KCNJ2, ABCA6, and ABCA9) associated with the four bone weight traits. Furthermore, our findings on the genetic mechanisms influencing pig bone weight offer valuable insights as a reference for the genetic enhancement of pig bone traits.

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.

Genome-Wide Association Studies and Runs of Homozygosity to Identify Reproduction-Related Genes in Yorkshire Pig Population

Reproductive traits hold considerable economic importance in pig breeding and production. However, candidate genes underpinning the reproductive traits are still poorly identified. In the present study, we executed a genome-wide association study (GWAS) and runs of homozygosity (ROH) analysis using the PorcineSNP50 BeadChip array for 585 Yorkshire pigs. Results from the GWAS identified two genome-wide significant and eighteen suggestive significant single nucleotide polymorphisms (SNPs) associated with seven reproductive traits. Furthermore, we identified candidate genes, including ELMO1, AOAH, INSIG2, NUP205, LYPLAL1, RPL34, LIPH, RNF7, GRK7, ETV5, FYN, and SLC30A5, which were chosen due to adjoining significant SNPs and their functions in immunity, fertilization, embryonic development, and sperm quality. Several genes were found in ROH islands associated with spermatozoa, development of the fetus, mature eggs, and litter size, including INSL6, TAF4B, E2F7, RTL1, CDKN1C, and GDF9. This study will provide insight into the genetic basis for pig reproductive traits, facilitating reproduction improvement using the marker-based selection methods.

Genome-wide association analysis unveils candidate genes and loci associated with aplasia cutis congenita in pigs

BACKGROUND

Aplasia cutis congenita (ACC) is a rare genetic disorder characterized by the localized or widespread absence of skin in humans and animals. Individuals with ACC may experience developmental abnormalities in the skeletal and muscular systems, as well as potential complications. Localized and isolated cases of ACC can be treated through surgical and medical interventions, while extensive cases of ACC may result in neonatal mortality. The presence of ACC in pigs has implications for animal welfare. It contributes to an elevated mortality rate among piglets at birth, leading to substantial economic losses in the pig farming industry. In order to elucidate candidate genetic loci associated with ACC, we performed a Genome-Wide Association Study analysis on 216 Duroc pigs. The primary goal of this study was to identify candidate genes that associated with ACC.

RESULTS

This study identified nine significant SNPs associated with ACC. Further analysis revealed the presence of two quantitative trait loci, 483 kb (5:18,196,971-18,680,098) on SSC 5 and 159 kb (13:20,713,440-207294431 bp) on SSC13. By annotating candidate genes within a 1 Mb region surrounding the significant SNPs, a total of 11 candidate genes were identified on SSC5 and SSC13, including KRT71, KRT1, KRT4, ITGB7, CSAD, RARG, SP7, PFKL, TRPM2, SUMO3, and TSPEAR.

CONCLUSIONS

The results of this study further elucidate the potential mechanisms underlying and genetic architecture of ACC and identify reliable candidate genes. These results lay the foundation for treating and understanding ACC in humans.

Identification of Candidate Genes for Economically Important Carcass Cutting in Commercial Pigs through GWAS

During the process of pork production, the carcasses of pigs are divided and sold, which provides better economic benefits and market competitiveness for pork production than selling the carcass as a whole. Due to the significant cost of post-slaughter phenotypic measurement, the genetic architecture of tenderloin weight (TLNW) and rib weight (RIBW)-important components of pig carcass economic value-remain unknown. In this study, we conducted genome-wide association studies (GWAS) for TLNW and RIBW traits in a population of 431 Duroc × Landrace × Yorkshire (DLY) pigs. In our study, the most significant single nucleotide polymorphism (SNP) associated with TLNW was identified as ASGA0085853 (3.28 Mb) on Sus scrofa chromosome 12 (SSC12), while for RIBW, it was Affx-1115046258 (172.45 Mb) on SSC13. Through haplotype block analysis, we discovered a novel quantitative trait locus (QTL) associated with TLNW, spanning a 5 kb region on SSC12, and a novel RIBW-associated QTL spanning 1.42 Mb on SSC13. Furthermore, we hypothesized that three candidate genes, TIMP2 and EML1, and SMN1, are associated with TLNW and RIBW, respectively. Our research not only addresses the knowledge gap regarding TLNW, but also serves as a valuable reference for studying RIBW. The identified SNP loci strongly associated with TLNW and RIBW may prove useful for marker-assisted selection in pig breeding programs.

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 Body Conformation Traits in a Three-Way Crossbred Commercial Pig Population

Body conformation is the most direct production index, which can fully reflect pig growth status and is closely related to critical economic traits. In this study, we conducted a genome-wide association study (GWAS) on body conformation traits in a population of 1518 Duroc × (Landrace × Yorkshire) commercial pigs. These traits included body length (BL), body height (BH), chest circumference (CC), abdominal circumference (AC), and waist circumference (WC). Both the mixed linear model (MLM) and fixed and random model circulating probability unification (FarmCPU) approaches were employed for the analysis. Our findings revealed 60 significant single nucleotide polymorphisms (SNPs) associated with these body conformation traits in the crossbred pig population. Specifically, sixteen SNPs were significantly associated with BL, three SNPs with BH, thirteen SNPs with CC, twelve SNPs with AC, and sixteen SNPs with WC. Moreover, we identified several promising candidate genes located within the genomic regions associated with body conformation traits. These candidate genes include INTS10, KIRREL3, SOX21, BMP2, MAP4K3, SOD3, FAM160B1, ATL2, SPRED2, SEC16B, and RASAL2. Furthermore, our analysis revealed a novel significant quantitative trait locus (QTL) on SSC7 specifically associated with waist circumference, spanning an 84 kb interval. Overall, the identification of these significant SNPs and potential candidate genes in crossbred commercial pigs enhances our understanding of the genetic basis underlying body conformation traits. Additionally, these findings provide valuable genetic resources for pig breeding programs.

Genome-Wide Association Study Identifies the Crucial Candidate Genes for Teat Number in Crossbred Commercial Pigs

The number of teats is a crucial reproductive trait with significant economic implications on maternal capacity and litter size. Consequently, improving this trait is essential to facilitate genetic selection for increased litter size. In this study, we performed a genome-wide association study (GWAS) of the number of teats in a three-way crossbred commercial Duroc × (Landrace × Yorkshire) (DLY) pig population comprising 1518 animals genotyped with the 50K BeadChip. Our analysis identified crucial quantitative trait loci (QTL) for the number of teats, containing the ABCD4 and VRTN genes on porcine chromosome 7. Our results establish SNP variants of ABCD4 and VRTN as new molecular markers for improving the number of teats in DLY pigs. Furthermore, the most significant noteworthy single nucleotide polymorphism (SNP) (7_97568284) was identified within the ABCD4 gene, exhibiting a significant association with the total teat number traits. This SNP accounted for a substantial proportion of the genetic variance, explaining 6.64% of the observed variation. These findings reveal a novel gene on SSC7 for the number of teats and provide a deeper understanding of the genetic mechanisms underlying reproductive traits.

Comparison of two multi-trait association testing methods and sequence-based fine mapping of six additive QTL in Swiss Large White pigs

BACKGROUND

Genetic correlations between complex traits suggest that pleiotropic variants contribute to trait variation. Genome-wide association studies (GWAS) aim to uncover the genetic underpinnings of traits. Multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS enable detecting variants associated with multiple phenotypes. In this study, we used array-derived genotypes and phenotypes for 24 reproduction, production, and conformation traits to explore differences between the two methods and used imputed sequence variant genotypes to fine-map six quantitative trait loci (QTL).

RESULTS

We considered genotypes at 44,733 SNPs for 5,753 pigs from the Swiss Large White breed that had deregressed breeding values for 24 traits. Single-trait association analyses revealed eleven QTL that affected 15 traits. Multi-trait association testing and the meta-analysis of the single-trait GWAS revealed between 3 and 6 QTL, respectively, in three groups of traits. The multi-trait methods revealed three loci that were not detected in the single-trait GWAS. Four QTL that were identified in the single-trait GWAS, remained undetected in the multi-trait analyses. To pinpoint candidate causal variants for the QTL, we imputed the array-derived genotypes to the sequence level using a sequenced reference panel consisting of 421 pigs. This approach provided genotypes at 16 million imputed sequence variants with a mean accuracy of imputation of 0.94. The fine-mapping of six QTL with imputed sequence variant genotypes revealed four previously proposed causal mutations among the top variants.

CONCLUSIONS

Our findings in a medium-size cohort of pigs suggest that multivariate association testing and the meta-analysis of summary statistics from single-trait GWAS provide very similar results. Although multi-trait association methods provide a useful overview of pleiotropic loci segregating in mapping populations, the investigation of single-trait association studies is still advised, as multi-trait methods may miss QTL that are uncovered in single-trait GWAS.

Integrated Single-Trait and Multi-Trait GWASs Reveal the Genetic Architecture of Internal Organ Weight in Pigs

Internal organ weight is an essential indicator of growth status as it reflects the level of growth and development in pigs. However, the associated genetic architecture has not been well explored because phenotypes are difficult to obtain. Herein, we performed single-trait and multi-trait genome-wide association studies (GWASs) to map the genetic markers and genes associated with six internal organ weight traits (including heart weight, liver weight, spleen weight, lung weight, kidney weight, and stomach weight) in 1518 three-way crossbred commercial pigs. In summation, single-trait GWASs identified a total of 24 significant single- nucleotide polymorphisms (SNPs) and 5 promising candidate genes, namely, TPK1, POU6F2, PBX3, UNC5C, and BMPR1B, as being associated with the six internal organ weight traits analyzed. Multi-trait GWAS identified four SNPs with polymorphisms localized on the APK1, ANO6, and UNC5C genes and improved the statistical efficacy of single-trait GWASs. Furthermore, our study was the first to use GWASs to identify SNPs associated with stomach weight in pigs. In conclusion, our exploration of the genetic architecture of internal organ weights helps us better understand growth traits, and the key SNPs identified could play a potential role in animal breeding programs.

Analysis of Copy Number Variation in the Whole Genome of Normal-Haired and Long-Haired Tianzhu White Yaks

Long-haired individuals in the Tianzhu white yak population are a unique genetic resource, and have important landscape value. Copy number variation (CNV) is an important source of phenotypic variation in mammals. In this study, we used resequencing technology to detect the whole genome of 10 long-haired Tianzhu white yaks (LTWY) and 10 normal-haired Tianzhu white yaks (NTWY), and analyzed the differences of CNV in the genome of LTWYs and NTWYs. A total of 110268 CNVs were identified, 2006 CNVRs were defined, and the distribution map of these CNVRs on chromosomes was constructed. The comparison of LTWYs and NTWYs identified 80 differential CNVR-harbored genes, which were enriched in lipid metabolism, cell migration and other functions. Notably, some differential genes were identified as associated with hair growth and hair-follicle development (e.g., ASTN2, ATM, COL22A1, GK5, SLIT3, PM20D1, and SGCZ). In general, we present the first genome-wide analysis of CNV in LTWYs and NTWYs. Our results can provide new insights into the phenotypic variation of different hair lengths in Tianzhu white yaks.

Genome-Wide Association Studies for Flesh Color and Intramuscular Fat in (Duroc × Landrace × Large White) Crossbred Commercial Pigs

The intuitive impression of pork is extremely important in terms of whether consumers are enthusiastic about purchasing it. Flesh color and intramuscular fat (IMF) are indispensable indicators in meat quality assessment. In this study, we determined the flesh color and intramuscular fat at 45 min and 12 h after slaughter (45 mFC, 45 mIMF, 12 hFC, and 12 hIMF) of 1518 commercial Duroc × Landrace × Large White (DLY) pigs. We performed a single nucleotide polymorphism (SNP) genome-wide association study (GWAS) analysis with 28,066 SNPs. This experiment found that the correlation between 45 mFC and 12 hFC was 0.343. The correlation between 45 mIMF and 12 hIMF was 0.238. The heritability of the traits 45 mFC, 12 hFC, 45 mIMF, and 12 hIMF was 0.112, 0.217, 0.139, and 0.178, respectively, and we identified seven SNPs for flesh color and three SNPs for IMF. Finally, several candidate genes regulating these four traits were identified. Three candidate genes related to flesh color were provided: SNCAIP and PRR16 on SSC2, ST3GAL4 on SSC5, and GALR1 on SSC1. A total of three candidate genes related to intramuscular fat were found, including ABLIM3 on SSC2, DPH5 on SSC4, and DOCK10 on SSC15. Furthermore, GO and KEGG analysis revealed that these genes are involved in the regulation of apoptosis and are implicated in functions such as pigmentation and skeletal muscle metabolism. This study applied GWAS to analyze the scoring results of flesh color and IMF in different time periods, and it further revealed the genetic structure of flesh color and IMF traits, which may provide important genetic loci for the subsequent improvement of pig meat quality traits.

A composite strategy of genome-wide association study and copy number variation analysis for carcass traits in a Duroc pig population

BACKGROUND

Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs.

RESULTS

Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function.

CONCLUSIONS

We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

A genome-wide association study of important reproduction traits in large white pigs

Augmenting the reproductive efficiency of sows remains the predominant challenge in the swine industry. This work was aimed at scrutinizing vital genetic markers for reproductive traits in this animal. This entailed probing of the records of vital attributes of Large White pigs (n = 695) inclusive of the total number of born (TNB), number of born alive (NBA), number of weaned pigs (NWP), number of healthy births (NHS), total litter weight of piglets born alive (BALWT), weaning litter weight (WNWT), and corrected litter weight at 21 days (W21). A genome-wide association study (GWAS) for the four litter traits and three traits of litter weight in the Denmark Large White population then ensued. We discovered seven significantly related SNPs and eleven potential candidate genes (e.g., TUSC3, THRB for TNB; STT3B for NBA). The subsequent functional enrichment analysis of these genes showed that the significant gene were associated with steroid hormone receptor activity. Our findings indicated that the genes TUSC3, THRB and STT3B probably contribute to litter traits in this population. This work reveals genetic mechanisms of reproduction traits and also supports ensuing genetic improvement employing marker-assisted selection in Large White pigs.

Identification of Homozygous Regions With Adverse Effects on the Five Economic Traits of Duroc Pigs

Runs of homozygosity (ROH) are widely used to estimate genomic inbreeding, which is linked to inbreeding depression on phenotypes. However, the adverse effects of specific homozygous regions on phenotypic characteristics are rarely studied in livestock. In this study, the 50 K SNP data of 3,770 S21 Duroc (American origin) and 2,096 S22 Duroc (Canadian origin) pigs were used to investigate the harmful ROH regions on five economic traits. The results showed that the two Duroc lines had different numbers and distributions of unfavorable ROHs, which may be related to the different selection directions and intensities between the two lines. A total of 114 and 58 ROH segments were found with significant adverse effects on the economic traits of S21 and S22 pigs, respectively. Serval pleiotropic ROHs were detected to reduce two or multiple phenotypic performances in two Duroc populations. Candidate genes in these shared regions were mainly related to growth, fertility, immunity, and fat deposition. We also observed that some ROH genotypes may cause opposite effects on different traits. This study not only enhances our understanding of the adverse effects of ROH on phenotypes, but also indicates that ROH information could be incorporated into breeding programs to estimate and control the detrimental effects of homozygous regions.

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.

Runs of Homozygosity Uncover Potential Functional-Altering Mutation Associated With Body Weight and Length in Two Duroc Pig Lines

Runs of homozygosity (ROH) are widely used to investigate genetic diversity, demographic history, and positive selection signatures of livestock. Commercial breeds provide excellent materials to reveal the landscape of ROH shaped during the intense selection process. Here, we used the GeneSeek Porcine 50K single-nucleotide polymorphism (SNP) Chip data of 3,770 American Duroc (AD) and 2,096 Canadian Duroc (CD) pigs to analyze the genome-wide ROH. First, we showed that AD had a moderate genetic differentiation with CD pigs, and AD had more abundant genetic diversity and significantly lower level of inbreeding than CD pigs. In addition, sows had larger levels of homozygosity than boars in AD pigs. These differences may be caused by differences in the selective intensity. Next, ROH hotspots revealed that many candidate genes are putatively under selection for growth, sperm, and muscle development in two lines. Population-specific ROHs inferred that AD pigs may have a special selection for female reproduction, while CD pigs may have a special selection for immunity. Moreover, in the overlapping ROH hotspots of two Duroc populations, we observed a missense mutation (rs81216249) located in the growth and fat deposition-related supergene (ARSB-DMGDH-BHMT) region. The derived allele of this variant originated from European pigs and was nearly fixed in Duroc pigs. Further selective sweep and association analyses indicated that this supergene was subjected to strong selection and probably contributed to the improvement of body weight and length in Duroc pigs. These findings will enhance our understanding of ROH patterns in different Duroc lines and provide promising trait-related genes and a functional-altering marker that can be used for genetic improvement of pigs.

Assessment of Heterozygosity and Genome-Wide Analysis of Heterozygosity Regions in Two Duroc Pig Populations

Heterozygosity can effectively reflect the diverse models of population structure and demographic history. However, the genomic distribution of heterozygotes and the correlation between regions of heterozygosity (runs of heterozygosity, ROHet) and phenotypes are largely understudied in livestock. The objective of this study was to identify ROHet in the Duroc pig genome, and investigate the relationships between ROHet and eight important economic traits. Here, we genotyped 3,770 American Duroc (S21) and 2,096 Canadian Duroc (S22) pigs using 50 K single nucleotide polymorphism array to analyze heterozygosity. A total of 145,010 and 84,396 ROHets were characterized for S21 and S22 populations, respectively. ROHet segments were mostly enriched in 1–2 Mb length classification (75.48% in S21 and 72.25% in S22). The average genome length covered by ROHet was 66.53 ± 12.20 Mb in S21 and 73.32 ± 13.77 Mb in S22 pigs. Additionally, we detected 20 and 13 ROHet islands in S21 and S22 pigs. Genes in these genomic regions were mainly involved in the biological processes of immunity and reproduction. Finally, the genome-wide ROHet-phenotypes association analysis revealed that 130 ROHets of S21 and 84 ROHets of S22 were significantly associated with eight economic traits. Among the candidate genes in the significant ROHet regions, 16 genes related to growth, metabolism, and meat quality were considered as candidate genes for important economic traits of pigs. This work preliminarily explores the effect of heterozygosity-rich regions in the pig genome on production performance and provides new insights for subsequent research on pig genetic improvement.

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.

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Genome-Wide Detection of Genetic Loci and Candidate Genes for Body Conformation Traits in Duroc × Landrace × Yorkshire Crossbred Pigs

The Duroc × (Landrace × Yorkshire) hybrid pigs (DLY) are the most popular commercial pigs, providing consumers with the largest source of pork. In order to gain more insights into the genetic architecture of economically important traits in pigs, we performed a genome-wide association study (GWAS) using the GeneSeek Porcine 50 K SNP Chip to map the genetic markers and genes associated with body conformation traits (BCT) in 311 DLY pigs. The quantitative traits analyzed included body weight (BW), carcass length (CL), body length (BL), body height (BH), and body mass index (BMI). BMI was defined as BMI_CL, BMI_BL, and BMI_BH, respectively, based on CL, BL, and BH phenotypic data. We identified 82 SNPs for the seven traits by GEMMA-based and FarmCPU-based GWASs. Both methods detected two quantitative trait loci (QTL) on SSC8 and SSC17 for body conformation traits. Several candidate genes (such as TNFAIP3, KDM4C, HSPG2, BMP2, PLCB4, and GRM5) were found to be associated with body weight and body conformation traits in pigs. Notably, the BMP2 gene had pleiotropic effects on CL, BL, BH, BMI_CL, and BMI_BL and is proposed as a strong candidate gene for body size due to its involvement in growth and bone development. Furthermore, gene set enrichment analysis indicated that most of the pathway terms are associated with regulation of cell growth, negative regulation of cell population proliferation, and chondrocyte differentiation. We anticipate that these results further advance our understanding of the genetic architecture of body conformation traits in the popular commercial DLY pigs and provide new insights into the genetic architecture of BMI in pigs.

Genomic Analyses Revealed the Genetic Difference and Potential Selection Genes of Growth Traits in Two Duroc Lines

Duroc pigs are famous for their high growth rate, feed conversion efficiency, and lean meat percentage. Given that they have been subjected to artificial selection and breeding in multiple countries, various lines with obvious differences in production performance have formed. In this study, we genotyped 3,770 American Duroc (AD) pigs and 2,098 Canadian Duroc (CD) pigs using the GeneSeek Porcine SNP50 Beadchip to dissect the genetic differences and potential selection genes of growth traits in these two Duroc pig lines. Population structure detection showed that there were significant genetic differences between the two Duroc pig lines. Hence, we performed F _ST and cross-population extended haplotype homozygosity (XP-EHH) analyses between the two lines. As a result, we identified 38 annotated genes that were significantly enriched in the gland development pathway in the AD line, and 61 annotated genes that were significantly enriched in the immune-related pathway in the CD line. For three growth traits including backfat thickness (BFT), loin muscle depth (LMD), and loin muscle area (LMA), we then performed selection signature detection at 5 and 10% levels within the line and identified different selected regions and a series of candidate genes that are involved in lipid metabolism and skeletal muscle development or repair, such as IRX3, EBF2, WNT10B, TLR2, PITX3, and SGCD. The differences in selected regions and genes between the two lines may be the cause of the differences in growth traits. Our study suggests significant genetic differences between the AD and CD lines, which provide a theoretical basis for selecting different Duroc lines as sires for different needs.

Genome-Wide Association Study for Body Length, Body Height, and Total Teat Number in Large White Pigs

Body length, body height, and total teat number are economically important traits in pig breeding, as these traits are usually associated with the growth, reproductivity, and longevity potential of piglets. Here, we report a genetic analysis of these traits using a population comprising 2,068 Large White pigs. A genotyping-by-sequencing (GBS) approach was used to provide high-density genome-wide SNP discovery and genotyping. Univariate and bivariate animal models were used to estimate heritability and genetic correlations. The results showed that heritability estimates for body length, body height, and total teat number were 0.25 ± 0.04, 0.11 ± 0.03, and 0.22 ± 0.04, respectively. The genetic correlation between body length and body height exhibited a strongly positive correlation (0.63 ± 0.15), while a positive but low genetic correlation was observed between total teat number and body length. Furthermore, we used two different genome-wide association study (GWAS) approaches: single-locus GWAS and weighted single-step GWAS (WssGWAS), to identify candidate genes for these traits. Single-locus GWAS detected 76, 13, and 29 significant single-nucleotide polymorphisms (SNPs) associated with body length, body height, and total teat number. Notably, the most significant SNP (S17_15781294), which is located 20 kb downstream of the BMP2 gene, explained 9.09% of the genetic variance for body length traits, and it also explained 9.57% of the genetic variance for body height traits. In addition, another significant SNP (S7_97595973), which is located in the ABCD4 gene, explained 8.92% of the genetic variance for total teat number traits. GWAS results for these traits identified some candidate genomic regions, such as SSC6: 14.96–15.02 Mb, SSC7: 97.18–98.18 Mb, SSC14: 128.29–131.15 Mb, SSC17: 15.39–17.27 Mb, and SSC17: 22.04–24.15 Mb, providing a starting point for further inheritance research. Most quantitative trait loci were detected by single-locus GWAS and WssGWAS. These findings reveal the complexity of the genetic mechanism of the three traits and provide guidance for subsequent genetic improvement through genome selection.

Genome-wide association studies reveals polygenic genetic architecture of litter traits in Duroc pigs

With continuous improvement of sow litter size, breeders are gradually paying more attention to the quality of litter traits that directly impact the production efficiency of pig companies, such as the rate of piglets born alive (RBA) and the rate of healthy births (RHB). The objectives of this study are to dissect the genetic basis of litter traits in pig and to identify valuable genes and genetic markers, especially pleiotropic, for pig breeding. Herein, 1140 Duroc pigs and 2046 reproduction records, 5 litter traits, including the number of healthy births (NHB), number of deformed fetuses (NDF), number of stillborn (NSB), RBA, and RHB, were used in this study. Subsequently, a genome-wide association study (GWAS) was performed for the five litter traits in the first two parities from two Duroc populations. A total of 76 significantly related SNPs and 10 potential candidate genes (CAV1, DAB2, FGF12, FHOD3, DYNC2H1, GRHL1, TCTN3, PYROXD2, MMP8, MMP13, and PGR) were detected, including 13 pleiotropic SNPs that affected more than one litter trait. Finally, the functional enrichment analysis of functional genes that were closest to these significant SNPs indicated that most of the significant pathways were associated with hormone secretion and embryo and organ development. This study advances our understanding of the genetic mechanisms of litter traits, especially the survival rate of piglets born, and provides an opportunity to increase the quality of litter using marker-assisted selection or genomic selection in pigs.

Accelerated deciphering of the genetic architecture of agricultural economic traits in pigs using a low-coverage whole-genome sequencing strategy

BACKGROUND

Uncovering the genetic architecture of economic traits in pigs is important for agricultural breeding. However, high-density haplotype reference panels are unavailable in most agricultural species, limiting accurate genotype imputation in large populations. Moreover, the infinitesimal model of quantitative traits implies that weak association signals tend to be spread across most of the genome, further complicating the genetic analysis. Hence, there is a need to develop new methods for sequencing large cohorts without large reference panels.

RESULTS

We describe a Tn5-based highly accurate, cost- and time-efficient, low-coverage sequencing method to obtain 11.3 million whole-genome single-nucleotide polymorphisms in 2,869 Duroc boars at a mean depth of 0.73×. On the basis of these single-nucleotide polymorphisms, a genome-wide association study was performed, resulting in 14 quantitative trait loci (QTLs) for 7 of 21 important agricultural traits in pigs. These QTLs harbour genes, such as ABCD4 for total teat number and HMGA1 for back fat thickness, and provided a starting point for further investigation. The inheritance models of the different traits varied greatly. Most follow the minor-polygene model, but this can be attributed to different reasons, such as the shaping of genetic architecture by artificial selection for this population and sufficiently interconnected minor gene regulatory networks.

CONCLUSIONS

Genome-wide association study results for 21 important agricultural traits identified 14 QTLs/genes and showed their genetic architectures, providing guidance for genetic improvement harnessing genomic features. The Tn5-based low-coverage sequencing method can be applied to large-scale genome studies for any species without a good reference panel and can be used for agricultural breeding.

A large-scale genome-wide association analysis reveals QTL and candidate genes for intramuscular fat content in Duroc pigs

This study aimed at identifying genomic regions and genes associated with intramuscular fat content (IMF) in Duroc pigs using a weighted single-step GWAS. Data from 3912 pigs, of which 3770 animals were genotyped with GeneSeek Porcine 50K Bead chip, were used for the association analysis. We identified 19 genomic regions that each explained >1% of the additive genetic variance associated with IMF. Notably, a consistent QTL on SSC7 (117.42-117.92 Mb) was confirmed, explaining 3.70% of the additive genetic variance, and two genes, BDKRB2 and ATG2B, were highlighted as promising candidates for IMF. Two QTL (SSC7, 94.19-94.64 Mb; SSC14, 123.25-123.75 Mb), which harbored MED6 and MAP3K9 genes and TCF7L2 gene respectively, were newly identified as associated with IMF. In conclusion, we identified a consistent QTL and additional genomic regions and genes that contributed to the genetic variance of IMF using a large-scale sample size of genotyped pigs and genealogical information.

Genome-wide detection of CNV regions and their potential association with growth and fatness traits in Duroc pigs

Background

In the process of pig breeding, the average daily gain (ADG), days to 100 kg (AGE), and backfat thickness (BFT) are directly related to growth rate and fatness. However, the genetic mechanisms involved are not well understood. Copy number variation (CNV), an important source of genetic diversity, can affect a variety of complex traits and diseases and has gradually been thrust into the limelight. In this study, we reported the genome-wide CNVs of Duroc pigs using SNP genotyping data from 6627 animals. We also performed a copy number variation region (CNVR)-based genome-wide association studies (GWAS) for growth and fatness traits in two Duroc populations.

Results

Our study identified 953 nonredundant CNVRs in U.S. and Canadian Duroc pigs, covering 246.89 Mb (~ 10.90%) of the pig autosomal genome. Of these, 802 CNVRs were in U.S. Duroc pigs with 499 CNVRs were in Canadian Duroc pigs, indicating 348 CNVRs were shared by the two populations. Experimentally, 77.8% of nine randomly selected CNVRs were validated through quantitative PCR (qPCR). We also identified 35 CNVRs with significant association with growth and fatness traits using CNVR-based GWAS. Ten of these CNVRs were associated with both ADG and AGE traits in U.S. Duroc pigs. Notably, four CNVRs showed significant associations with ADG, AGE, and BFT, indicating that these CNVRs may play a pleiotropic role in regulating pig growth and fat deposition. In Canadian Duroc pigs, nine CNVRs were significantly associated with both ADG and AGE traits. Further bioinformatic analysis identified a subset of potential candidate genes, including PDGFA, GPER1, PNPLA2 and BSCL2.

Conclusions

The present study provides a necessary supplement to the CNV map of the Duroc genome through large-scale population genotyping. In addition, the CNVR-based GWAS results provide a meaningful way to elucidate the genetic mechanisms underlying complex traits. The identified CNVRs can be used as molecular markers for genetic improvement in the molecular-guided breeding of modern commercial pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07654-7.

Genetic dissection of flowering time in flax (Linum usitatissimum L.) through single- and multi-locus genome-wide association studies

In a rapidly changing climate, flowering time (FL) adaptation is important to maximize seed yield in flax (Linum usitatissimum L.). However, our understanding of the genetic mechanism underlying FL in this multipurpose crop remains limited. With the aim of dissecting the genetic architecture of FL in flax, a genome-wide association study (GWAS) was performed on 200 accessions of the flax core collection evaluated in four environments. Two single-locus and six multi-locus models were applied using 70,935 curated single nucleotide polymorphism (SNP) markers. A total of 40 quantitative trait nucleotides (QTNs) associated with 27 quantitative trait loci (QTL) were identified in at least two environments. The number of QTL with positive-effect alleles in accessions was significantly correlated with FL (r = 0.77 to 0.82), indicating principally additive gene actions. Nine QTL were significant in at least three of the four environments accounting for 3.06-14.71% of FL variation. These stable QTL spanned regions that harbored 27 Arabidopsis thaliana and Oryza sativa FL-related orthologous genes including FLOWERING LOCUS T (Lus10013532), FLOWERING LOCUS D (Lus10028817), transcriptional regulator SUPERMAN (Lus10021215), and gibberellin 2-beta-dioxygenase 2 (Lus10037816). In silico gene expression analysis of the 27 FL candidate gene orthologous suggested that they might play roles in the transition from vegetative to reproductive phase, flower development and fertilization. Our results provide new insights into the QTL architecture of flowering time in flax, identify potential candidate genes for further studies, and demonstrate the effectiveness of combining different GWAS models for the genetic dissection of complex traits.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

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

Simple Summary

Number of teats is very important for lactating sows. We conducted genome-wide association studies (GWAS) and estimated the genetic parameters for traits related to teat number. Results showed that there were nine and 22 SNPs exceeding genome-wide significance and suggestive significance levels, respectively. Eighteen genes annotated near them were concentrated on chromosomes 7 and 10. Among them, three new candidate genes were located on the genomic regions around the significant SNPs. Our findings provide new insight into investigating the complex genetic mechanism of traits related to teat number in pigs.

Abstract

The number of teats is related to the nursing ability of sows. In the present study, we conducted genome-wide association studies (GWAS) for traits related to teat number in Duroc pig population. Two mixed models, one for counted and another for binary phenotypic traits, were employed to analyze seven traits: the right (RTN), left (LTN), and total (TTN) teat numbers; maximum teat number on a side (MAX); left minus right side teat number (LR); the absolute value of LR (ALR); and the presence of symmetry between left and right teat numbers (SLR). We identified 11, 1, 4, 13, and 9 significant SNPs associated with traits RTN, LTN, MAX, TTN, and SLR, respectively. One significant SNP (MARC0038565) was found to be simultaneous associated with RTN, LTN, MAX and TTN. Two annotated genes (VRTN and SYNDIG1L) were located in genomic region around this SNP. Three significant SNPs were shown to be associated with TTN, RTN and MAX traits. Seven significant SNPs were simultaneously detected in two traits of TTN and RTN. Other two SNPs were only identified in TTN. These 13 SNPs were clustered in the genomic region between 96.10—98.09 Mb on chromosome 7. Moreover, nine significant SNPs were shown to be significantly associated with SLR. In total, four and 22 SNPs surpassed genome-wide significance and suggestive significance levels, respectively. Among candidate genes annotated, eight genes have documented association with the teat number relevant traits. Out of them, DPF3 genes on Sus scrofa chromosome (SSC) 7 and the NRP1 gene on SSC 10 were new candidate genes identified in this study. Our findings demonstrate the genetic mechanism of teat number relevant traits and provide a reference to further improve reproductive performances in practical pig breeding programs.

The Association of an SNP in the EXOC4 Gene and Reproductive Traits Suggests Its Use as a Breeding Marker in Pigs

In mammals, the exocyst complex component 4 (EXOC4) gene has often been reported to be involved in vesicle transport. The SNP rs81471943 (C/T) is located in the intron of porcine EXOC4, while six quantitative trait loci (QTL) within 5-10 Mb around EXOC4 are associated with ovary weight, teat number, total offspring born alive, and corpus luteum number. However, the molecular mechanisms between EXOC4 and the reproductive performance of pigs remains to be elucidated. In this study, rs81471943 was genotyped from a total of 994 Duroc sows, and the genotype and allele frequency of SNP rs81471943 (C/T) were statistically analyzed. Then, the associations between SNP rs81471943 and four reproductive traits, including number of piglets born alive (NBA), litter weight at birth (LWB), number of piglets weaned (NW), and litter weight at weaning (LWW), were determined. Sanger sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) were utilized to identify the rs81471943 genotype. We found that the genotype frequency of CC was significantly higher than that of CT and TT, and CC was the most frequent genotype for NBA, LWB, NW, and LWW. Moreover, 5'-deletion and luciferase assays identified a positive transcription regulatory element in the EXOC4 promoter. After exploring the EXOC4 promoter, SNP -1781G/A linked with SNP rs81471943 (C/T) were identified by analysis of the transcription activity of the haplotypes, and SNP -1781 G/A may influence the potential binding of P53, E26 transformation specific sequence -like 1 transcription factor (ELK1), and myeloid zinc finger 1 (MZF1). These findings provide useful information for identifying a molecular marker of EXOC4-assisted selection in pig breeding.

Weighted Single-Step GWAS Identified Candidate Genes Associated with Growth Traits in a Duroc Pig Population

Growth traits are important economic traits of pigs that are controlled by several major genes and multiple minor genes. To better understand the genetic architecture of growth traits, we performed a weighted single-step genome-wide association study (wssGWAS) to identify genomic regions and candidate genes that are associated with days to 100 kg (AGE), average daily gain (ADG), backfat thickness (BF) and lean meat percentage (LMP) in a Duroc pig population. In this study, 3945 individuals with phenotypic and genealogical information, of which 2084 pigs were genotyped with a 50 K single-nucleotide polymorphism (SNP) array, were used for association analyses. We found that the most significant regions explained 2.56–3.07% of genetic variance for four traits, and the detected significant regions (>1%) explained 17.07%, 18.59%, 23.87% and 21.94% for four traits. Finally, 21 genes that have been reported to be associated with metabolism, bone growth, and fat deposition were treated as candidate genes for growth traits in pigs. Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses implied that the identified genes took part in bone formation, the immune system, and digestion. In conclusion, such full use of phenotypic, genotypic, and genealogical information will accelerate the genetic improvement of growth traits in pigs.

A meta-analysis of genome-wide association studies for average daily gain and lean meat percentage in two Duroc pig populations

Background

Average daily gain (ADG) and lean meat percentage (LMP) are the main production performance indicators of pigs. Nevertheless, the genetic architecture of ADG and LMP is still elusive. Here, we conducted genome-wide association studies (GWAS) and meta-analysis for ADG and LMP in 3770 American and 2090 Canadian Duroc pigs.

Results

In the American Duroc pigs, one novel pleiotropic quantitative trait locus (QTL) on Sus scrofa chromosome 1 (SSC1) was identified to be associated with ADG and LMP, which spans 2.53 Mb (from 159.66 to 162.19 Mb). In the Canadian Duroc pigs, two novel QTLs on SSC1 were detected for LMP, which were situated in 3.86 Mb (from 157.99 to 161.85 Mb) and 555 kb (from 37.63 to 38.19 Mb) regions. The meta-analysis identified ten and 20 additional SNPs for ADG and LMP, respectively. Finally, four genes (PHLPP1, STC1, DYRK1B, and PIK3C2A) were detected to be associated with ADG and/or LMP. Further bioinformatics analysis showed that the candidate genes for ADG are mainly involved in bone growth and development, whereas the candidate genes for LMP mainly participated in adipose tissue and muscle tissue growth and development.

Conclusions

We performed GWAS and meta-analysis for ADG and LMP based on a large sample size consisting of two Duroc pig populations. One pleiotropic QTL that shared a 2.19 Mb haplotype block from 159.66 to 161.85 Mb on SSC1 was found to affect ADG and LMP in the two Duroc pig populations. Furthermore, the combination of single-population and meta-analysis of GWAS improved the efficiency of detecting additional SNPs for the analyzed traits. Our results provide new insights into the genetic architecture of ADG and LMP traits in pigs. Moreover, some significant SNPs associated with ADG and/or LMP in this study may be useful for marker-assisted selection in pig breeding.