A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Effects of hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 polymorphisms on fat androstenone level and gene expression in Duroc pigs

A high level of androstenone in porcine adipose tissue is a major factor contributing to boar taint. Porcine hydroxy‐delta‐5‐steroid dehydrogenase, 3 beta‐ and steroid delta‐isomerase 1 (3β‐HSD, also known as HSD3B1) plays a key role in the hepatic metabolism that catalyzes androstenone to β‐androstenol. Therefore, 3β‐HSD is a candidate gene for boar taint. This study aimed to investigate functional 3β‐HSD polymorphisms in Duroc pigs. We found eight single nucleotide polymorphisms (SNPs) in the full‐length porcine 3β‐HSD. Four of the SNPs had restriction enzyme sites, and we genotyped them in 147 uncastrated male Duroc pigs using a polymerase chain reaction–restriction fragment length polymorphism method. Pigs with the GG genotype at the g.165262G>A locus (SNP5) had significantly lower androstenone levels than did those with other genotypes (P = 0.030). SNP5 also was associated with differences in 3β‐HSD mRNA levels: pigs with the GG genotype had higher levels than those with other genotypes (P = 0.019). The SNP5 polymorphism could affect the hepatic catabolism of androstenone and consequently impact androstenone accumulation in the adipose tissue. Therefore, SNP5 in the 3β‐HSD of Duroc pigs could be a useful selective marker for decreasing boar taint.     

Genome‐wide association study for the level of serum electrolytes in Italian Large White pigs

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome‐wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance‐tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome‐wide efficient mixed‐model association. For the level of Ca²⁺, eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (P_{nominal value} < 5.00E‐05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.     

Molecular Characterization of Five Porcine Candidate Genes for Drip Loss in Pork

Drip loss is the loss of fluid from a piece of meat without mechanical force and represents an important meat quality trait. Previous work revealed a quantitative trait locus (QTL) for drip loss in pork in an experimental Duroc x Piétrain (DUPI) F2 family on SSC 5. Based on functional data indicating their possible involvement in water holding capacity and their expression in skeletal muscle, we selected five positional candidates (ACO2, ADSL, CBY1, KCNJ4, PLA2AG6) out of 130 predicted genes in the QTL interval for further analysis. We performed a mutation analysis of all coding exons and discovered 204 polymorphisms. We genotyped 39 single nucleotide polymorphisms (SNPs) in 192 Piétrain pigs with extreme drip loss phenotypes and detected a possible association with drip loss for one non-coding SNP in the ADSL gene (ss107793818, p_raw = 0.021). Correspondingly, ADSL diplotypes were associated with drip loss and pH1 of M. longissimus dorsi. However, after correction for multiple testing, none of the tested SNPs were significantly associated with drip loss. One possible explanation for these results is that one of the QTL-alleles from the experimental DUPI family may be fixed or nearly fixed in the tested Piétrain population.     

Genome-Wide Association Study Singles Out SCD and LEPR as the Two Main Loci Influencing Intramuscular Fat Content and Fatty Acid Composition in Duroc Pigs

Intramuscular fat (IMF) content and fatty acid composition affect the organoleptic quality and nutritional value of pork. A genome-wide association study was performed on 138 Duroc pigs genotyped with a 60k SNP chip to detect biologically relevant genomic variants influencing fat content and composition. Despite the limited sample size, the genome-wide association study was powerful enough to detect the association between fatty acid composition and a known haplotypic variant in SCD (SSC14) and to reveal an association of IMF and fatty acid composition in the LEPR region (SSC6). The association of LEPR was later validated with an independent set of 853 pigs using a candidate quantitative trait nucleotide. The SCD gene is responsible for the biosynthesis of oleic acid (C18:1) from stearic acid. This locus affected the stearic to oleic desaturation index (C18:1/C18:0), C18:1, and saturated (SFA) and monounsaturated (MUFA) fatty acids content. These effects were consistently detected in gluteus medius, longissimus dorsi, and subcutaneous fat. The association of LEPR with fatty acid composition was detected only in muscle and was, at least in part, a consequence of its effect on IMF content, with increased IMF resulting in more SFA, less polyunsaturated fatty acids (PUFA), and greater SFA/PUFA ratio. Marker substitution effects estimated with a subset of 65 animals were used to predict the genomic estimated breeding values of 70 animals born 7 years later. Although predictions with the whole SNP chip information were in relatively high correlation with observed SFA, MUFA, and C18:1/C18:0 (0.48–0.60), IMF content and composition were in general better predicted by using only SNPs at the SCD and LEPR loci, in which case the correlation between predicted and observed values was in the range of 0.36 to 0.54 for all traits. Results indicate that markers in the SCD and LEPR genes can be useful to select for optimum fatty acid profiles of pork.     

Genome‐wide association analysis for growth,muscularity and meat quality in Piétrain pigs

Improvement in growth and meat quality is one of the main objectives in sire line pig breeding programmes. Mapping quantitative trait loci for these traits using experimental crosses and a linkage‐based approach has been performed frequently in the past. The Piétrain breed often was involved as a founder breed to establish the experimental crosses. This breed was selected for muscularity and leanness but shows relatively poor meat quality. It is frequently used as a sire line breed. With the advent of genome‐wide and dense SNP chips in pig genomic research, it is possible to also conduct genome‐wide association studies within the Piétrain breed. In this study, around 500 progeny‐tested sires were genotyped with 60k SNPs. Data filtering showed that around 48k SNPs were useable in this sample. These SNPs were used to conduct a genome‐wide association study for growth, muscularity and meat quality traits. Because it is known that a mutation in the RYR1 gene located on chromosome 6 shows a major effect on meat quality, this mutation was included in the models. Single‐marker and multimarker association analyses were performed. The results revealed between zero and eight significant associations per trait with P < 5 × 10^?5. Of special interest are SNPs located on SSC6, SSC10 and SSC15.     

Evaluation of quantitative trait loci affecting intramuscular fat and reproductive traits in pigs using marker‐assisted introgression

We investigated the effects of previously identified quantitative trait loci (QTL) in an experimental backcross (BC) between Chinese Meishan pigs and commercial Duroc pigs. We performed marker‐assisted introgression of two QTL for intramuscular fat (IMF) content (IMF population) and three QTL for reproductive traits (reproduction population) from a donor Meishan pig into a recipient Duroc pig. At the fourth BC generation of the IMF population and third BC generation of the reproduction population, carrier animals were selected for the production of animals homozygous for the QTL. Our previous studies have shown that the presence of a Meishan allele on the IMF QTL is associated with low IMF values, and the Meishan allele on the reproductive QTL is associated with large litters. In this study, the presence of a Duroc allele at the IMF QTL on SSC9 resulted in a 0.27% increase in IMF (additive effect = 0.27 ± 0.08), whereas the presence of a Meishan allele at the IMF QTL on SSC7 resulted in a 0.34% increase in IMF (additive effect = ?0.34 ± 0.09). The presence of the Meishan allele at the IMF QTL on SSC7 thus had the opposite effect to our previous studies, that is, increased IMF. In the reproduction population, we observed no differences between the genotypes of the three QTL in regard to number of corpora lutea or litter size. Marker‐assisted introgression at these QTL is thus unlikely to result in an associated increase in litter size. These results show that it is possible to introgress alleles from other breeds into a selection population using molecular markers; any unexpected results might be associated with the genetic background.     

Genome‐wide association study reveals regions associated with gestation length in two pig populations

Reproduction traits, such as gestation length (GLE), play an important role in dam line breeding in pigs. The objective of our study was to identify single nucleotide polymorphisms (SNPs) that are associated with GLE in two pig populations. Genotypes and deregressed breeding values were available for 2081 Dutch Landrace‐based (DL) and 2301 Large White‐based (LW) pigs. We identified two QTL regions for GLE, one in each population. For DL, three associated SNPs were detected in one QTL region spanning 0.52 Mbp on Sus scrofa chromosome (SSC) 2. For LW, four associated SNPs were detected in one region of 0.14 Mbp on SSC5. The region on SSC2 contains the heparin‐binding EGF‐like growth factor (HBEGF) gene, which promotes embryo implantation and has been described to be involved in embryo survival throughout gestation. The associated SNP can be used for marker‐assisted selection in the studied populations, and further studies of the HBEGF gene are warranted to investigate its role in GLE.     

A genome‐wide association study for a proxy of intermuscular fat level in the Italian Large White breed identifies genomic regions affecting an important quality parameter for dry‐cured hams

Intermuscular fat content in protected designations of origin dry‐cured hams is a very important meat quality trait that affects the acceptability of the product by the consumers. An excess in intermuscular fat (defined as the level of fat deposition between leg muscles) is a defect that depreciates the final product. In this study we carried out a genome‐wide association study for visible intermuscular fat (VIF) of hams in the Italian Large White pig breed. This trait was evaluated on the exposed muscles of green legs in 1122 performance‐tested gilts by trained personnel, according to a classification scale useful for routine and cheap evaluation. All animals were genotyped with the Illumina PorcineSNP60 BeadChip. The genome‐wide association study identified three QTL regions on porcine chromosome 1 (SSC1; accounting for ~79% of the SNPs below the 5.0E?04 threshold) and SSC2, two on SSC7 and one each on SSC3, SSC6, SSC9, SSC11, SSC13, SSC15, SSC16 and SSC17. The most significant SNP (ALGA0004143 on SSC1 at 77.3 Mb; P_FDR < 0.05), included in the largest QTL region which spanned about 6.8 Mb on SSC1, is located within the glutamate ionotropic receptor kainate type subunit 2 (GRIK2) gene. Functional annotation of all genes included in QTL regions for VIF suggested that intermuscular fat in the Italian Large White breed is a complex trait apparently influenced by complex biological mechanisms also involving obesity‐related processes. These QTL target mainly chromosome regions different from those affecting subcutaneous and intramuscular fat deposition.     

Non-Synonymous SNP in the ApoR Gene Associated with Pork Meat Quality

Two novel non-synonymous SNPs in the 2^nd and 3^rd exons of the porcine ApoR gene are reported. One was identified as a novel SNP significantly associated with multiple traits of pork meat quality. The data can provide a useful resource for developing a marker in the genetic improvement of pigs.     

A genome‐wide association study of copy number variations with umbilical hernia in swine

Umbilical hernia (UH) is one of the most common congenital defects in pigs, leading to considerable economic loss and serious animal welfare problems. To test whether copy number variations (CNVs) contribute to pig UH, we performed a case–control genome‐wide CNV association study on 905 pigs from the Duroc, Landrace and Yorkshire breeds using the Porcine SNP60 BeadChip and penncnv algorithm. We first constructed a genomic map comprising 6193 CNVs that pertain to 737 CNV regions. Then, we identified eight CNVs significantly associated with the risk for UH in the three pig breeds. Six of seven significantly associated CNVs were validated using quantitative real‐time PCR. Notably, a rare CNV (CNV14:13030843–13059455) encompassing the NUGGC gene was strongly associated with UH (permutation‐corrected P = 0.0015) in Duroc pigs. This CNV occurred exclusively in seven Duroc UH‐affected individuals. SNPs surrounding the CNV did not show association signals, indicating that rare CNVs may play an important role in complex pig diseases such as UH. The NUGGC gene has been implicated in human omphalocele and inguinal hernia. Our finding supports that CNVs, including the NUGGC CNV, contribute to the pathogenesis of pig UH.     

Genomic regions affecting backfat thickness and cannon bone circumference identified by genome‐wide association study in a Duroc pig population

We performed a genome‐wide association study using the porcine 60K SNP array to detect QTL regions for nine traits in a three‐generational Duroc samples (n = 651), viz. generations 1, 2 and 3 from a population selected over five generations using a closed nucleus breeding scheme. We applied a linear mixed model for association mapping to detect SNP effects, adjusting for fixed effects (sex and season) and random polygenic effects (reflecting genetic relatedness), and derived a likelihood ratio statistic for each SNP using the efficient mixed‐model association method. We detected a region on SSC6 for backfat thickness (BFT) and on SSC7 for cannon bone circumference (CANNON), with a genome‐wide significance of P < 0.01 after Bonferroni correction. These regions had been detected previously in other pig populations. Six genes are located in the BFT‐associated region, while the CANNON‐associated region includes 66 genes. In the future, significantly associated SNPs, derived by sequencing the coding regions of the six genes in the BFT region, can be used in marker‐assisted selection of BFT, whereas haplotypes constructed from the SSC7 region with strong LD can be used to select for the CANNON trait in our resource family.     

Integrating genome-wide association study with RNA-seq revealed DBI as a good candidate gene for intramuscular fat content in Beijing black pigs

Increasing intramuscular fat (IMF) content can enhance the sensory quality of meat, including tenderness, juiciness, flavor, and color. Genome-wide association study and RNA-sequencing (RNA-seq) analysis were used to identify candidate IMF genes in Beijing Black pigs, a popular species among consumers in northern China. Two and three single nucleotide polymorphisms were significantly associated with IMF in SSC13 and SSC15 respectively. Solute carrier family 4 member 7 (SLC4A7) on SSC13 and insulin induced gene 2 (INSIG2), coiled-coil domain containing 93 (CCDC93), and diazepam binding inhibitor acyl-CoA binding protein (DBI) on SSC15 are good candidate genes in this population. Furthermore, RNA-seq analysis was performed between high and low IMF groups, and 534 differentially expressed genes were identified. In addition, based on differentially expressed genes, Kyoto Encyclopedia of Genes and Genomes analysis revealed that peroxisome proliferator-activated receptors and FoxO signaling pathway pathways might contribute to IMF. Moreover, the DBI gene was identified as a candidate for IMF both by genome-wide association study and RNA-seq analysis, suggesting that it might be a crucial candidate gene for influencing IMF in Beijing Black pigs.     

Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population

Pork quality is an economically important trait and one of the main selection criteria for breeding in the swine industry. In this genome-wide association study (GWAS), 455 pigs from a porcine Large White × Minzhu intercross population were genotyped using the Illumina PorcineSNP60K Beadchip, and phenotyped for intramuscular fat content (IMF), marbling, moisture, color L*, color a*, color b* and color score in the longissimus muscle (LM). Association tests between each trait and the SNPs were performed via the Genome Wide Rapid Association using the Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approach. From the Ensembl porcine database, SNP annotation was implemented using Sus scrofa Build 9. A total of 45 SNPs showed significant association with one or multiple meat quality traits. Of the 45 SNPs, 36 were located on SSC12. These significantly associated SNPs aligned to or were in close approximation to previously reported quantitative trait loci (QTL) and some were located within introns of previously reported candidate genes. Two haplotype blocks ASGA0100525-ASGA0055225-ALGA0067099-MARC0004712-DIAS0000861, and ASGA0085522-H3GA0056170 were detected in the significant region. The first block contained the genes MYH1, MYH2 and MYH4. A SNP (ASGA0094812) within an intron of the USP43 gene was significantly associated with five meat quality traits. The present results effectively narrowed down the associated regions compared to previous QTL studies and revealed haplotypes and candidate genes on SSC12 for meat quality traits in pigs.     

Genome‐wide association analysis identifies quantitative trait loci for growth in a Landrace purebred population

Growth‐related traits are complex and economically important in the livestock industry. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting growth in pigs. A genome‐wide association study (GWAS) was performed using the porcine single‐nucleotide polymorphism (SNP) 60K bead chip. A mixed‐effects model and linear regression approach were used for the GWAS. The data used in the study included 490 purebred Landrace pigs. All experimental animals were genotyped with 39 438 SNPs located throughout the pig autosomes. We identified a strong association between a SNP marker on chromosome 16 and body weight at 71 days of age (ALGA0092396, P = 5.35 × 10^?9, Bonferroni adjusted P < 0.05). The SNP marker was located near the genomic region containing IRX4, which encodes iroquois homeobox 4. This SNP marker could be useful in the selective breeding program after validating its effect on other populations.     

Genome-wide identification of quantitative trait loci for carcass composition and meat quality in a large-scale White Duroc × Chinese Erhualian resource population

Carcass and meat quality traits are economically important in pigs. In this study, 17 carcass composition traits and 23 meat quality traits were recorded in 1028 F₂ animals from a White Duroc × Erhualian resource population. All pigs in this experimental population were genotyped for 194 informative markers covering the entire porcine genome. Seventy-seven genome-wide significant quantitative trait loci (QTL) for carcass traits and 68 for meat quality were mapped to 34 genomic regions. These results not only confirmed many previously reported QTL but also revealed novel regions associated with the measured traits. For carcass traits, the most prominent QTL was identified for carcass length and head weight at 57 cM on SSC7, which explained up to 50% of the phenotypic variance and had a 95% confidence interval of only 3 cM. Moreover, QTL for kidney and spleen weight and lengths of cervical vertebrae were reported for the first time in pigs. For meat quality traits, two significant QTL on SSC5 and X were identified for both intramuscular fat content and marbling score in the longissimus muscle, while three significant QTL on SSC1 and SSC9 were found exclusively for IMF. Both LM and the semimembranous muscle showed common QTL for colour score on SSC4, 5, 7, 8, 13 and X and discordant QTL on other chromosomes. White Duroc alleles at a majority of QTL detected were favourable for carcass composition, while favourable QTL alleles for meat quality originated from both White Duroc and Erhualian.     

Evaluation of the causality of the low‐density lipoprotein receptor gene (LDLR) for serum lipids in pigs

A significant quantitative trait locus (QTL) for low‐density lipoprotein cholesterol (LDL‐C) and total cholesterol (TC) was identified around the LDLR gene on chromosome 2 (SSC2) in a White Duroc × Erhualian F₂ resource population and Sutai pigs in our previous study. However, in previous reports, the causality of LDLR with serum lipids is controversial in pigs. To systematically assess the causality of LDLR with serum lipids, association analyses were successively performed in three populations: Sutai pigs, a White Duroc × Erhualian F₂ resource population and a Duroc × (Landrace × Large White) population. We first performed a haplotype‐based association study with 60K SNP genotyping data and evidenced the significant association with LDL‐C and TC around the LDLR gene region. We also found that there is more than one QTL for LDL‐C and TC on SSC2. Then, we evaluated the causalities of two missense mutations, c.1812C>T and c.1520A>G, with LDL‐C and TC. We revealed that the c.1812C>T SNP showed the strongest association with LDL‐C (P = 5.40 × 10^?11) and TC (P = 3.64 × 10^?8) and explained all the QTL effect in Sutai pigs. Haplotype analysis found that two missense SNPs locate within a 1.93‐Mb haplotype block. One major haplotype showed the strongest significant association with LDL‐C (P = 4.62 × 10^?18) and TC (P = 1.06 × 10^?9). However, the c.1812C>T SNP was not identified in the White Duroc × Erhualian intercross, and the association of c.1520A>G with both LDL‐C and TC did not achieve significance in this F₂ population, suggesting population heterogeneity. Both missense mutations were identified in the Duroc × (Landrace × Large White) population and showed significant associations with LDL‐C and TC. Our data give evidence that the LDLR gene should be a candidate causative gene for LDL‐C and TC in pigs, but heterogeneity exists in different populations.