Increasing feed efficiency is a major goal of breeders as it can reduce production cost and energy consumption. However, the genetic architecture of feeding behavior and feed efficiency traits remains elusive. To investigate the genetic architecture of feed efficiency in pigs, three feeding behavior traits (daily feed intake, number of daily visits to feeder, and duration of each visit) and two feed efficiency traits (feed conversion ratio and residual feed intake) were considered. We performed genome-wide association studies (GWASs) of the five traits using a population of 1,008 Duroc pigs genotyped with an Illumina Porcine SNP50K BeadChip. A total of 9 genome-wide (P < 1.54E-06) and 35 suggestive (P < 3.08E-05) single nucleotide polymorphisms (SNPs) were detected. Two pleiotropic quantitative trait loci (QTLs) on SSC 1 and SSC 7 were found to affect more than one trait. Markers WU_10.2_7_18377044 and DRGA0001676 are two key SNPs for these two pleiotropic QTLs. Marker WU_10.2_7_18377044 on SSC 7 contributed 2.16 and 2.37% of the observed phenotypic variance for DFI and RFI, respectively. The other SNP DRGA0001676 on SSC 1 explained 3.22 and 5.46% of the observed phenotypic variance for FCR and RFI, respectively. Finally, functions of candidate genes and gene set enrichment analysis indicate that most of the significant pathways are associated with hormonal and digestive gland secretion during feeding. This study advances our understanding of the genetic mechanisms of feeding behavior and feed efficiency traits and provide an opportunity for increasing feeding efficiency using marker-assisted selection or genomic selection in pigs.

Gut microbiota has indispensable roles in nutrient digestion and energy harvesting, especially in processing the indigestible components of dietary polysaccharides. Searching for the microbial taxa and functional capacity of the gut microbiome associated with feed efficiency (FE) can provide important knowledge to increase profitability and sustainability of the swine industry. In the current study, we performed a comparative analysis of the fecal microbiota in 50 commercial Duroc × (Landrace × Yorkshire) (DLY) pigs with polarizing FE using 16S rRNA gene sequencing and shotgun metagenomic sequencing. There was a different microbial community structure in the fecal microbiota of pigs with different FE. Random forest analysis identified 24 operational taxonomic units (OTUs) as potential biomarkers to improve swine FE. Multiple comparison analysis detected 8 OTUs with a significant difference or tendency toward a difference between high- and low-FE pigs (P < 0.01, q < 0.1). The high-FE pigs had a greater abundance of OTUs that were from the Lachnospiraceae and Prevotellaceae families and the Escherichia-Shigella and Streptococcus genera than low-FE pigs. A sub-species Streptococcus gallolyticus subsp. gallolyticus could be an important candidate for improving FE. The functional capacity analysis found 18 KEGG pathways and CAZy EC activities that were different between high- and low-FE pigs. The fecal microbiota in high FE pigs have greater functional capacity to degrade dietary cellulose, polysaccharides, and protein and may have a greater abundance of microbes that can promote intestinal health. These results provided insights for improving porcine FE through modulating the gut microbiome.

Background More teats are necessary for sows to nurse larger litters to provide immunity and nutrient for piglets prior to weaning. Previous studies have reported the strong effect of an insertion mutation in the Vertebrae Development Associated (VRTN) gene on Sus scrofa chromosome 7 (SSC7) that increased the number of thoracic vertebrae and teat number in pigs. We used genome-wide association studies (GWAS) to map genetic markers and genes associated with teat number in two Duroc pig populations with different genetic backgrounds. A single marker method and several multi-locus methods were utilized. A meta-analysis that combined the effects and P-values of 34,681 single nucleotide polymorphisms (SNPs) that were common in the results of single marker GWAS of American and Canadian Duroc pigs was conducted. We also performed association tests between the VRTN insertion and teat number in the same populations. Results A total of 97 SNPs were found to be associated with teat number. Among these, six, eight and seven SNPs were consistently detected with two, three and four multi-locus methods, respectively. Seven SNPs were concordantly identified between single marker and multi-locus methods. Moreover, 26 SNPs were newly found by multi-locus methods to be associated with teat number. Notably, we detected one consistent quantitative trait locus (QTL) on SSC7 for teat number using single-locus and meta-analysis of GWAS and the top SNP (rs692640845) explained 8.68% phenotypic variance of teat number in the Canadian Duroc pigs. The associations between the VRTN insertion and teat number in two Duroc pig populations were substantially weaker. Further analysis revealed that the effect of VRTN on teat number may be mediated by its LD with the true causal mutation. Conclusions Our study suggested that VRTN insertion may not be a strong or the only candidate causal mutation for the QTL on SSC7 for teat number in the analyzed Duroc pig populations. The combination of single-locus and multi-locus GWAS detected additional SNPs that were absent using only one model. The identified SNPs will be useful for the genetic improvement of teat number in pigs by assigning higher weights to associated SNPs in genomic selection.

Intramuscular fat (IMF) is an important quantitative trait of meat, which affects the associated sensory properties and nutritional value of pork. To gain a better understanding of the genetic determinants of IMF, we used a composite strategy, including single-locus and multi-locus association analyses to perform genome-wide association studies (GWAS) for IMF in 1,490 Duroc boars. We estimated the genomic heritability of IMF to be 0.23 ± 0.04. A total of 30 single nucleotide polymorphisms (SNPs) were found to be significantly associated with IMF. The single-locus mixed linear model (MLM) and multiple-locus methods multi-locus random-SNP-effect mixed linear model (mrMLM), fast multi-locus random-SNP-effect efficient mixed model association (FASTmrEMMA), and integrative sure independence screening expectation maximization Bayesian least absolute shrinkage and selection operator model (ISIS EM-BLASSO) analyses identified 5, 9, 8, and 21 significant SNPs, respectively. Interestingly, a novel quantitative trait locus (QTL) on SSC 7 was found to affect IMF. In addition, 10 candidate genes ( BDKRB2 , GTF2IRD1 , UTRN , TMEM138 , DPYD , CASQ2 , ZNF518B , S1PR1 , GPC6 , and GLI1 ) were found to be associated with IMF based on their potential functional roles in IMF. GO analysis showed that most of the genes were involved in muscle and organ development. A significantly enriched KEGG pathway, the sphingolipid signaling pathway, was reported to be associated with fat deposition and obesity. Identification of novel variants and functional genes will advance our understanding of the genetic mechanisms of IMF and provide specific opportunities for marker-assisted or genomic selection in pigs. In general, such a composite single-locus and multi-locus strategy for GWAS may be useful for understanding the genetic architecture of economic traits in livestock.

Loin muscle area (LMA) and loin muscle depth (LMD) are important traits influencing the production performance of breeding pigs. However, the genetic architecture of these two traits is still poorly understood. To discern the genetic architecture of LMA and LMD, a material consisting of 6043 Duroc pigs belonging to two populations with different genetic backgrounds was collected and applied in genome-wide association studies (GWAS) with a genome-wide distributed panel of 50K single nucleotide polymorphisms (SNPs). To improve the power of detection for common SNPs, we conducted a meta-analysis in these two pig populations and uncovered additional significant SNPs. As a result, we identified 75 significant SNPs for LMA and LMD on SSC6, 7, 12, 16, and 18. Among them, 25 common SNPs were associated with LMA and LMD. One pleiotropic quantitative trait locus (QTL), which was located on SSC7 with a 283 kb interval, was identified to affect LMA and LMD. Marker ALGA0040260 is a key SNP for this QTL, explained 1.77% and 2.48% of the phenotypic variance for LMA and LMD, respectively. Another genetic region on SSC16 (709 kb) was detected and displayed prominent association with LMA and the peak SNP, WU_10.2_16_35829257, contributed 1.83% of the phenotypic variance for LMA. Further bioinformatics analysis determined eight promising candidate genes ( GCLC , GPX8 , DAXX , FGF21 , TAF11 , SPDEF , NUDT3 , and PACSIN1 ) with functions in glutathione metabolism, adipose and muscle tissues development and lipid metabolism. This study provides the first GWAS for the LMA and LMD of Duroc breed to analyze the underlying genetic variants through a large sample size. The findings further advance our understanding and help elucidate the genetic architecture of LMA, LMD and growth-related traits in pigs.

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.

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 FST 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.

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.