Weighted gene co-expression network analysis reveals potential candidate genes affecting drip loss in pork

Drip loss is an essential evaluation indicator for pork quality. It is closely related to other meat quality indicators, including water-holding capacity, water loss rate and pH value at 45 min (pH₁) and 24 h post-mortem (pH₂), and is influenced by environmental and genetic factors and their interactions. We previously conducted differentially expressed gene analysis to identify candidate genes affecting drip loss using eight individuals with extremely high- and low-drip loss selected from 28 purebred Duroc pigs. Using 28 identical samples, in the present study, we performed weighted gene co-expression network analysis with drip loss and drip loss-related traits, including water-holding capacity, water loss rate, pH₁ and pH₂. A total of 25 modules were identified, and five of them correlated with at least two drip loss or drip loss-related traits. After functional enrichment analysis of genes in the five modules, three modules were found to be critical, as their genes were significantly involved in amino acid metabolism, immune response and apoptosis, which have potential relationships with drip loss. Furthermore, we identified five candidate genes affecting drip loss in one critical module, AASS, BCKDHB, ALDH6A1, MUT and MCCC1, as they overlapped with differentially expressed genes detected in our previous study, exhibited protein–protein interactions and had potential biological functions in affecting drip loss according to the literature. The outcomes of the present study enhance our understanding of the molecular mechanisms underlying drip loss and will aid in improving the pork quality.     

Associations of functional candidate genes derived from gene-expression profiles of prenatal porcine muscle tissue with meat quality and muscle deposition

Ten genes (ANK1, bR10D1, CA3, EPOR, HMGA2, MYPN, NME1, PDGFRA, ERC1, TTN), whose candidacy for meat-quality and carcass traits arises from their differential expression in prenatal muscle development, were examined for association in 1700 performance-tested fattening pigs of commercial purebred and crossbred herds of Duroc, Pietrain, Pietrain x (Landrace x Large White), Duroc x (Landrace x Large White) as well as in an experimental F(2) population based on a reciprocal cross of Duroc and Pietrain. Comparative sequencing revealed polymorphic sites segregating across commercial breeds. Genetic mapping results corresponded to pre-existing assignments to porcine chromosomes or current human-porcine comparative maps. Nine of these genes showed association with meat-quality and carcass traits at a nominal P-value of < or = 0.05; PDGFRA revealed no association reaching the P < or = 0.05 threshold. In particular, HMGA2, CA3, EPOR, NME1 and TTN were associated with meat colour, pH and conductivity of loin 24 h postmortem; CA3 and MYPN exhibited association with ham weight and lean content (FOM) respectively at P-values of < 0.003 that correspond to false discovery rates of < 0.05. However, none of the genes showed significant associations for a particular trait across all populations. The study revealed statistical-genetic evidence for association of the functional candidate genes with traits related to meat quality and muscle deposition. The polymorphisms detected are not likely causal, but markers were identified that are in linkage disequilibrium with causal genetic variation within particular populations.     

Sire breed comparisons for meat and eating quality traits in Australian pig populations

Sire breed comparisons for carcase (n = 1169), meat and eating quality (n = 686) traits were obtained using data recorded on progeny of known pedigree sired by Duroc (DU), Large White (LW), Landrace (LR) and Duroc Synthetic (DS) boars from crossbred (LW/LR) sows. Animals were reared in eco-shelters in large single-sex contemporary groups and slaughtered on an age constant basis at 22 weeks of age. Compared to progeny from other sire groups, animals sired by purebred Duroc boars tended to have improved eating quality traits (higher intramuscular fat and lower shear force) at the expense of poorer carcase characteristics (higher subcutaneous and belly fat). Animals sired by DS boars tended to be heavier and leaner than those sired by LW, LR or DU boars; intramuscular fat but not belly fat was correspondingly lower, while tenderness was generally consistent with that of the DU progeny. Significant variability of sire progeny groups within sire breed suggests that sire breed selection, potentially used for improving traits such as meat quality in commercial progeny, will be less accurate in the absence of sire-specific information, which is typically poorly recorded in this class of traits.     

Identification of genes related to intramuscular fat independent of backfat thickness in Duroc pigs using single‐step genome‐wide association

Intramuscular fat (IMF) is an important meat‐quality trait of pigs, which influences pork’s shearing force, hydraulics, tenderness and juicy flavor. However, to achieve a higher percentage of lean meat, pigs with lower backfat thickness (BF) are intensively selected for, which may lead to a reduction in pork quality. Therefore, the objective of this study was to locate loci that affect IMF without changing BF. A single‐step GWAS was performed on 950 Duroc pigs genotyped by a 50K SNP chip in order to detect genomic variants relevant to IMF and BF. The significant SNPs detected were afterwards divided into a BF subset (seven SNPs), an IMF subset (11 SNPs) and a subset of both traits (12 SNPs), according to their P‐value and LD. After SNP and QTL annotation, our results indicated that SSC1: 167938652, 166363826, 164829874 and 167171587 might be associated with IMF without changing BF. In the subset of both traits, we found that the combined effect of ALGA0006602 (SSC1: 159538854) and 12784636 (SSC1: 160773437) might improve the IMF without changing BF. Our gene annotation result showed that TLE3, ITGA11, SMAD6, PAQR5 and [RNF152^A/G × MC4R^A/A] genes might affect IMF independently of BF. We believe that the SNPs and genes identified in this study will be valuable for the future molecular breeding of IMF in Duroc pigs.     

SNP detection and genetic mapping of porcine genes encoding enzymes in hepatic metabolic pathways and evaluation of linkage with carcass traits

We have previously identified and mapped porcine expressed sequence tags (ESTs) derived from genes that are preferentially expressed in liver. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) in porcine genes encoding enzymes in hepatic metabolic pathways and use the SNPs for mapping. Furthermore, these genes, which are involved in utilization and partitioning of nutrients, were examined for their effects on carcass and meat quality traits by linkage analyses. In total, 100 ESTs were screened for SNPs by single strand conformation polymorphism analyses across a diverse panel of animals with a 36% success rate. Twelve of 36 polymorphic loci segregated in a three-generation Duroc x Berlin Miniature Pig (F2) resource population, the DUMI resource population, and were genetically mapped. Interval mapping of the corresponding chromosomes was performed to verify mapping of the genes within quantitative trait loci (QTL) regions detected in this resource population. QTL with genome-wide significance were detected in the vicinity of GNMT, ESTL147 and HGD. These loci therefore are positional candidate genes.     

Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2

Principal component analysis of traits related to carcass and meat properties were combined with microarray expression data for the identification of functional networks of genes and biological processes taking place during the conversion of muscle to meat. Principal components (PCs) with high loadings of meat quality traits were derived from phenotypic data of 572 animals of a porcine crossbreed population. Microarray data of 74 M. longissimus dorsi samples were correlated with PC datasets. Lists of significantly correlated genes were analyzed for enrichment of functional annotation groups as defined in the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases as well as the Ingenuity Pathways Analysis library. Ubiquitination, phosphorylation, mitochondrion dysfunction, actin, integrin, platelet-derived growth factor, epidermal growth factor, vascular endothelial growth factor, and Ca signaling pathways are correlated with meat quality. Among the significantly trait-associated genes, CAPZB, ANKRD1, and CTBP2 are promoted as candidate genes for meat quality that provide a link between the highlighted pathways. Knowledge of the relevant biological processes and the differential expression of members of the pathway will provide tools that are predictive for traits related to meat quality and that may also be diagnostic for many muscle defects or damages including muscle atrophy, dystrophy, and hypoxia.     

Evidence for new alleles in the protein kinase adenosine monophosphate-activated gamma(3)-subunit gene associated with low glycogen content in pig skeletal muscle and improved meat quality.

Several quantitative trait loci (QTL) affecting muscle glycogen content and related traits were mapped to pig chromosome 15 using a three-generation intercross between Berkshire x Yorkshire pigs. On the basis of the QTL location the PRKAG3 (protein kinase, AMP-activated, gamma(3)-subunit) gene was considered to be a good candidate for the observed effects. Differences in the PRKAG3 gene sequences of the founder animals of the intercross were analyzed. The RN(-) mutation previously reported was not present in the cross but three missense substitutions and a polymorphic short interspersed element (SINE) were identified. To confirm the hypothesis that at least one of these mutations was associated with differences in meat quality, >1800 animals from several unrelated commercial lines were genotyped for the candidate substitutions and an association study was performed. The results demonstrate the presence of new economically important alleles of the PRKAG3 gene affecting the glycogen content in the muscle and the resulting meat quality. Haplotype analysis was shown to resolve the effects of PRKAG3 more clearly than analysis of individual polymorphisms. Because of their prevalence in the more common commercial breeds, the potential implications for the pig industry and consumers are considerably greater than the original discovery of the RN(-) mutation. Furthermore, these results illustrate that additional alleles of genes involved in major mutations may play a significant role in quantitative trait variation.     

Investigation of candidate genes for meat quality in dry-cured ham production: the porcine cathepsin B (CTSB) and cystatin B (CSTB) genes

Excessive softness is a serious defect of dry cured hams which seems related to high activity of lysosomal cysteine proteinases, such as cathepsin B, in fresh pork muscles a few days after slaughtering. As it has been shown that cathepsin B activity has a moderate heritability in Italian Large White pigs we started a candidate gene approach to identify the gene(s) that affect(s) this parameter. Here, we studied two candidate genes: cathepsin B (CTSB) and cystatin B (CSTB). We amplified and sequenced porcine DNA fragments for these two genes that were used to identify polymorphisms by SSCP and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis. Four and two alleles were detected at the CTSB and CSTB loci, respectively. Sequencing of the CSTB alleles showed a missense mutation that changes a codon for aspartic acid into a codon for asparagine in exon 3 of the gene. Allele frequencies for the two loci differed among the pig breeds studied (Large White, Landrace, Duroc, Belgian Landrace, Hampshire, Piétrain, Meishan, Cinta Senese, Casertana, Calabrese and Nero di Sicilia). Linkage, somatic cell hybrid panel and radiation hybrid panel analyses assigned CTSB to porcine chromosome (Sscr) 14 and CSTB to Sscr 13. The markers identified at the CTSB and CSTB loci were used in association studies with several traits of economic importance including parameters that may indicate the suitability of pig meat to produce dry-cured hams. Significant associations were observed between CTSB and back-fat thickness and between CSTB and average daily gain. In this study, cathepsin B activity was not associated with the polymorphisms identified at the CTSB and CSTB loci.     

Genome-wide association studies for seven production traits highlight genomic regions useful to dissect dry-cured ham quality and production traits in Duroc heavy pigs

Protected designation of origin dry-cured hams are obtained from heavy pigs (slaughtered at about 160 kg of live weight). A specific breeding program designed to improve meat quality for this production has included as key traits the level of intermuscular fat between the leg muscles and ham weight loss during the seasoning period together with a balance between fat and lean cuts. In this study we carried out genome-wide association studies for seven traits used in the genetic merit of Italian Duroc heavy pigs, five related to meat and carcass quality traits (visible intermuscular fat, ham weight loss at first salting, backfat thickness, ham weight and lean cuts), and two related to performance and efficiency traits (average daily gain and feed : gain ratio). A total of 573 performance-tested pigs were genotyped with the Illumina PorcineSNP60 BeadChip and genome-wide association analyses were carried out using the Bayes B approach with the 1 Mb window option of GenSel and random residuals for each of the seven traits. Detected windows were supported by independent single nucleotide polymorphism analyses with a linear mixed model (LMM) approach on the same animals for the same traits. A total of 30 windows identifying different quantitative trait loci (QTL) were detected and among those, 27 were confirmed by LMM in one of these traits. Among the confirmed windows, three QTL were reported for visible intermuscular fat, seven for ham weight loss at first salting and five and four for backfat thickness and lean cut, respectively. A total of eight QTL were detected for the other production traits. No overlapping QTL were reported except for one window on porcine chromosome 10 between lean cuts and ham weight that contained the CACNB2 gene that has been already associated with loin marbling score in other Duroc pigs. Several regions contained genes that have been already associated with production traits in other pig breeds, including Duroc lines, related to fat deposition or muscle structure. This work reports, for the first time, genome-wide association study results for several traits in Italian Duroc heavy pigs. These results will be useful to dissect the genetic basis for dry-cured ham production traits that determine the total genetic merit index of Italian Duroc pigs.     

Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle

Background

Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation.

Results

Following adjustment for false discovery (q-value < 0.05), 479 quantitative trait loci (QTL) were associated with at least one of the four carcass traits using a single SNP regression approach. Using a Bayesian approach, 46 QTL were associated (posterior probability > 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway.

Conclusions

A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-837) contains supplementary material, which is available to authorized users.     

Genome-wide association study reveals novel loci associated with fat-deposition and meat-quality traits in Pekin ducks

Meat-quality traits play an essential role in meat poultry production. To determine the genetic mechanisms of meat quality in Pekin ducks, we performed a large-scale GWAS to identify quantitative trait loci affecting meat quality in Pekin ducks. We measured 10 traits in 542 Pekin ducks and genotyped each duck using genotyping-by-sequencing. The genetic parameters (genomic heritability, genetic correlation) for 10 meat-quality related traits were evaluated. Based on the large genotype–phenotype dataset, we performed GWASs for all of these traits. A total of 33 significant QTL (P < 3.03 × 10⁻⁵) across 13 chromosomes were identified by loci-based analysis. Some newly identified candidate genes were discovered for fat-deposition and meat-quality traits, including PAG1 for body weight and eviscerated weight, INTU and NUP35 for abdominal fat weight and ratio, NUP3 and ARHGDIB for skin fat weight and ratio, GOLGA5 for breast muscle toughness and breast tenderness, and CTDSPL and PKP1 for breast muscle thickness. The current study is the first systematic report regarding duck meat quality.     

Polymorphisms on SSC15q21-q26 Containing QTL for reproduction in Swine and its association with litter size

Several quantitative trait loci (QTL) for important reproductive traits (ovulation rate) have been identified on the porcine chromosome 15 (SSC15). To assist in the selection of positional candidate swine genes for these QTL on SSC15, twenty-one genes had already been assigned to SSC15 in a previous study in our lab, by using the radiation hybrid panel IMpRH. Further polymorphism studies were carried out on these positional candidate genes with four breeds of pigs (Duroc, Erhualian, Dahuabai and Landrace) harboring significant differences in reproduction traits. A total of nineteen polymorphisms were found in 21 genes. Among these, seven in six genes were used for association studies, whereby NRP2 polymorphism was found to be significantly (p < 0.05) associated with litter-size traits. NRP2 might be a candidate gene for pig-litter size based on its chromosome location (Du et al., 2006), significant association with litter-size traits and relationships with Sema and the VEGF super families.     

Identification of endogenous microRNA references in porcine serum for quantitative real-time PCR normalization

MicroRNAs (miRNAs) are evolutionarily conserved small non-coding RNAs that regulate the expression of genes, and they affect important biological and physiological states. Circulating miRNAs in blood are useful markers of metabolism and economic traits. Expression levels of circulating miRNAs have been estimated using quantitative real-time PCR (qPCR). Proper normalization is critical for accurate miRNA expression analysis. However, there is no study which systematically presented endogenous reference genes for evaluating circulating miRNA expression in pigs. In this study, ten porcine miRNAs (let-7a, miR-16, miR-17, miR-23a, miR-26a, miR-93, miR-103, miR-107, miR-127 and miR-191), based on the literature, were chosen as candidate reference miRNAs in serum. We evaluated the expression stability value of these miRNAs in Berkshire, Duroc, Landrace and Yorkshire pigs using geNorm and NormFinder. We determined the optimal combination of reference miRNAs for qPCR experiments: miR-127 and miR-17 in Berkshire pigs; miR-127 and miR-93 in Duroc and Landrace pigs; miR-127 and miR-16 in Yorkshire pigs. miR-127 was the best reference gene in pigs, regardless of the breed. Our study is crucial for the discovery of novel biomarkers in pigs. The reference miRNAs presented in this study could be used as appropriate reference genes for the measurement of circulating miRNA levels in studies of physiological blood metabolites.