Comprehensive Analysis of Transcriptome and Metabolome Reveals Regulatory Mechanism of Intramuscular Fat Content in Beef Cattle

The intramuscular fat (IMF) content of beef determined the meat quality, and the market value of beef varies with different breeds. To provide some new approaches for improving meat quality and cattle breed improvement, 24-month-old Qinchuan cattle (Q, n = 6), Nanyang cattle (N, n = 6), and Japanese black cattle (J, n = 6) were selected. IMF content of the J group (16.92 ± 1.08%) is remarkably higher than that of indigenous Chinese cattle (Q, 13.38 ± 1.08%, and N, 12.35 ± 1.22%). Monounsaturated fatty acids and polyunsaturated fatty acids in the J group are higher than the Q and creatine, lysine, and glutamine are the three most abundant amino acids in beef, which contribute to the flavor formation. Similarly, IMF content-related genes were enriched in four vital KEGG pathways, including fatty acid metabolism, biosynthesis of unsaturated fatty acids, fatty acid elongation, and insulin resistance. Moreover, weighted genes coexpression network analysis (WGCNA) revealed that ITGB1 is the critical gene associated with the IMF content. This study compares transcriptome and metabolome of local and high-IMF cattle breeds, providing data for native cattle breeding and improvement of beef quality.

Genome-wide population structure, homozygosity, and heterozygosity patterns of Nero Siciliano pig in the framework of Italian and cosmopolitan breeds

Analysis of genomic data is becoming more and more common for the effective management of livestock breeding programmes, even in the case of local populations. In this work, the genome-wide data of Nero Siciliano pig breed were compared to that of wild boar, Italian local and cosmopolitan breeds to investigate its genetic structure, and runs of homozygosity (ROH) and heterozygosity patterns. The Nero Siciliano has been reported to have the highest rate of genetic diversity among the Italian breeds, and a genetic variability comparable to that of the cosmopolitan breeds. Analyses of genomic structure and relationships underlined its proximity to wild boar, and an internal substructure probably linked to different family lines. The breed showed a low value of inbreeding estimated from ROH, and the highest diversity index among the Italian breeds, even if lower than that of the cosmopolitans. Four ROH islands in three chromosomes (SSC8, SSC11, and SSC14) and one heterozygosity-rich region (SSC1) were identified in Nero Siciliano, highlighting genomic regions related to productive QTL. Across breeds, SSC8 and SSC14 were the chromosomes with most ROH islands, with Mora Romagnola and wild boar showing the highest level of autozygosity. Chromosomes SSC2, SSC6, SSC8 and SSC13 showed the majority of runs of heterozygosity regions, mainly found in the cosmopolitan pig breeds, which reported several genes associated with health-related QTL. The outlined results can help to better identify the genomic profile of this local breed in order to plan matings, maintain adequate internal diversity and exploit the production system.

Transcriptome Analysis Reveals the Age-Related Developmental Dynamics Pattern of the Longissimus Dorsi Muscle in Ningxiang Pigs

The growth and development of the Longissimus Dorsi muscle are complex, playing an important role in the determination of pork quality. The study of the Longissimus Dorsi muscle at the mRNA level is particularly crucial for finding molecular approaches to improving meat quality in pig breeding. The current study utilized transcriptome technology to explore the regulatory mechanisms of muscle growth and intramuscular fat (IMF) deposition in the Longissimus Dorsi muscle at three core developmental stages (natal stage on day 1, growing stage on day 60, and finishing stage on day 210) in Ningxiang pigs. Our results revealed 441 differentially expressed genes (DEGs) in common for day 1 vs. day 60 and day 60 vs. day 210, and GO (Gene Ontology) analysis showed that candidate genes RIPOR2, MEGF10, KLHL40, PLEC, TBX3, FBP2, and HOMER1 may be closely related to muscle growth and development, while KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis showed that DEGs (UBC, SLC27A5, RXRG, PRKCQ, PRKAG2, PPARGC1A, PLIN5, PLIN4, IRS2, and CPT1B) involved the PPAR (Peroxisome Proliferator-Activated Receptor) signaling pathway and adipocytokine signaling pathway, which might play a pivotal role in the regulation of IMF deposition. PPI (Protein-Protein Interaction Networks) analysis found that the STAT1 gene was the top hub gene. Taken together, our results provide evidence for the molecular mechanisms of growth and development and IMF deposition in Longissimus Dorsi muscle to optimize carcass mass.

Genome-wide association study reveals genetic loci and candidate genes for meat quality traits in a four-way crossbred pig population

Meat quality traits (MQTs) have gained more attention from breeders due to their increasing economic value in the commercial pig industry. In this genome-wide association study (GWAS), 223 four-way intercross pigs were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) and phenotyped for PH at 45 min post mortem (PH45), meat color score (MC), marbling score (MA), water loss rate (WL), drip loss (DL) in the longissimus muscle, and cooking loss (CL) in the psoas major muscle. A total of 227, 921 filtered single nucleotide polymorphisms (SNPs) evenly distributed across the entire genome were detected to perform GWAS. A total of 64 SNPs were identified for six meat quality traits using the mixed linear model (MLM), of which 24 SNPs were located in previously reported QTL regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43% to 16.32%. The genomic heritability estimates based on SNP for six meat-quality traits were low to moderate (0.07-0.47) being the lowest for CL and the highest for DL. A total of 30 genes located within 10 kb upstream or downstream of these significant SNPs were found. Furthermore, several candidate genes for MQTs were detected, including pH45 (GRM8), MC (ANKRD6), MA (MACROD2 and ABCG1), WL (TMEM50A), CL (PIP4K2A) and DL (CDYL2, CHL1, ABCA4, ZAG and SLC1A2). This study provided substantial new evidence for several candidate genes to participate in different pork quality traits. The identification of these SNPs and candidate genes provided a basis for molecular marker-assisted breeding and improvement of pork quality traits.

Combining genome-wide association study based on low-coverage whole genome sequencing and transcriptome analysis to reveal the key candidate genes affecting meat color in pigs

Meat color is an attractive trait that influences consumers' purchase decisions at the point of sale. To decipher the genetic basis of meat color traits, we performed a genome-wide association study based on low-coverage whole-genome sequencing. In total, 669 (Pietrain × Duroc) × (Landrace × Yorkshire) pigs were genotyped using low-coverage whole-genome sequencing. Single nucleotide polymorphism (SNP) calling and genotype imputation were performed using the BaseVar + STITCH channel. Six individuals with an average depth of 12.05× whole-genome resequencing were randomly selected to assess the accuracy of imputation. Heritability evaluation and genome-wide association study for meat color traits were conducted. Functional enrichment analysis of the candidate genes from genome-wide association study and integration analysis with our previous transcriptome data were conducted. The imputation accuracy parameters, allele frequency R² , concordance rate, and dosage R² were 0.959, 0.952, and 0.933, respectively. The heritability values of a*_45 min , b*_45 min , L*_45 min , C*, and H⁰ were 0.19, 0.11, 0.06, 0.16, and 0.26, respectively. In total, 3884 significant SNPs and 15 QTL, corresponding to 382 genes, were associated with meat color traits. Functional enrichment analysis revealed that 10 genes were the potential candidates for regulating meat color. Moreover, integration analysis revealed that DMRT2, EFNA5, FGF10, and COL11A2 were the most promising candidates affecting meat color. In summary, this study provides new insights into the molecular basis of meat color traits, and provides a new theoretical basis for the molecular breeding of meat color traits in pigs.

Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems

Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.

The comprehensive detection of miRNA and circRNA in the regulation of intramuscular and subcutaneous adipose tissue of Laiwu pig

circRNAs, as miRNA sponges, participate in many important biological processes. However, it remains unclear whether circRNAs can regulate lipid metabolism. This study aimed to explore the competing endogenouse RNA (ceRNA) regulatory network that affects the difference between intramuscular fat (IMF) and subcutaneous fat (SCF) deposition, and to screen key circRNAs and their regulatory genes. In this experiment, we identified 265 differentially expressed circRNAs, of which 187 up-regulated circRNA and 78 down-regulated circRNA in IMF. Subsequently, we annotated the function of DEcircRNA's host genes, and found that DEcircRNA's host genes were mainly involved in GO terms (including cellular response to fatty acids, lysophosphatidic acid acyltransferase activity, R-SMAD binding, etc.) and signaling pathways (fatty acid biosynthesis, Citrate cycle, TGF- β Signal pathway) related to adipogenesis, differentiation and lipid metabolism. By constructing a circRNA-miRNA network, we screened out DEcircRNA that can competitively bind to more miRNAs as key circRNAs (circRNA_06424 and circRNA_08840). Through the functional annotation of indirect target genes and protein network analysis, we found that circRNA_06424 affects the expression of PPARD, MMP9, UBA7 and other indirect target genes by competitively binding to miRNAs such as ssc-miR-339-5p, ssc-miR-744 and ssc-miR-328, and participates in PPAR signaling pathway, Wnt signaling pathway, unsaturated fatty acid and other signaling pathways, resulting in the difference of fat deposition between IMF and SCF. This study provide a theoretical basis for further research investigating the differences of lipid metabolism in different adipose tissues, providing potential therapeutic targets for ectopic fat deposition and lipid metabolism diseases.

A polymorphism in porcine miR-22 is associated with pork color

MicroRNAs (miRNAs) are posttranscriptional regulators that play key roles in meat color regulation. Changes in miRNA expression affect their target mRNAs, leading to multifunctional effects on biological processes and phenotypes. In this study, a G > A mutation site located upstream of the precursor miR-22 sequence in Suhuai pigs was significantly correlated with the meat color parameter a^*(redness) of the porcine longissimus dorsi (LD) muscle. AA genotype individuals had the highest average meat color a^* value and the lowest miR-22 level. When G > A mutation was performed in the miR-22 overexpression vector, miR-22 expression significantly decreased. Considering that Ca²⁺ homeostasis is closely related to pig meat color, our results further demonstrated that ELOVL6 is a direct target of miR-22 in pigs. The effects of miR-22 on skeletal muscle intracellular Ca²⁺ were partially caused by the suppression of ELOVL6 expression.

Cross-population selection signatures in Canchim composite beef cattle

Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.

Effects of Neck-Arm Restraint Suspension of Beef Carcasses on Meat Quality and Proteome of Different Muscles During Post-mortem Aging

Beef quality is the first deciding factor for consumers to consider before purchasing. The aim of this study was to evaluate the effects of suspension and aging time on beef quality. We compared the differences in pH, drip loss, cooking loss, color, shear force, myofibril fragmentation index (MFI), and electron microscope of three muscle tissues between Achilles tendon (AT) and neck-arm restraint (NR) suspensions during seven aging periods (days 0, 1, 2, 3, 7, 14, and 21) after slaughter using the carcasses of six Xinjiang brown cattle. We found that NR suspension could significantly increase the water loss rate and MFI, as well as reduce the shear force compared to AT suspension. The muscle fiber structure with NR suspension was more severely damaged. The proteomics of longissimus dorsi was checked for the post-mortem days 1, 7, and 14. We detected 50, 26, and 29 differentially expressed proteins between NR and AT suspension at post-mortem days 1, 7, and 14, respectively. These proteins were involved in metabolic and muscle structure associated pathways and contributed to a comprehensive understanding of suspension-dependent meat quality regulation by proteins in beef cattle. To conclude, NR suspension can accelerate the aging time of beef carcasses, which will reduce the cost of carcass suspension and bring more benefits in the beef industry.

Whole-Genome Selective Scans Detect Genes Associated With Important Phenotypic Traits in Sheep (Ovis aries)

Sheep (Ovis aries) is one of the important livestock with diverse phenotypic traits. However, little is known about the molecular mechanism of diverse phenotypic traits in domestic sheep. Using the genome-wide high-density SNP data (600K) in 253 samples from 13 populations, we conducted the tests of selective sweeps (i.e., pairwise F_ST and XP-CLR) associated with several important phenotypic traits (e.g., tail types, horn morphology, prolificacy, coat pigmentation, ear size, milk production, meat production, body size and wool fineness). We identified strong selective signatures in previously reported (e.g., T, RXFP2, BMPR1B, TYRP1, MSRB3, TF, CEBPA, GPR21 and HOXC8) and novel genes associated with the traits, such as CERS6, BTG1, RYR3, SLC6A4, NNAT and OGT for fat deposition in the tails, FOXO4 for fertility, PTCH1 and EMX2 for ear size, and RMI1 and SCD5 for body size. Further gene annotation analysis showed that these genes were identified to be the most probable genes accounting for the diverse phenotypic traits. Our results provide novel insights into the genetic mechanisms underlying the traits and also new genetic markers for genetic improvement in sheep and other livestock.

An Integrative Analysis of Transcriptome and GWAS Data to Identify Potential Candidate Genes Influencing Meat Quality Traits in Pigs

Understanding the genetic factors behind meat quality traits is of great significance to animal breeding and production. We previously conducted a genome-wide association study (GWAS) for meat quality traits in a White Duroc × Erhualian F2 pig population using Illumina porcine 60K SNP data. Here, we further investigate the functional candidate genes and their network modules associated with meat quality traits by integrating transcriptomics and GWAS information. Quantitative trait transcript (QTT) analysis, gene expression QTL (eQTL) mapping, and weighted gene co-expression network analysis (WGCNA) were performed using the digital gene expression (DGE) data from 493 F2 pig's muscle and liver samples. Among the quantified 20,108 liver and 23,728 muscle transcripts, 535 liver and 1,014 muscle QTTs corresponding to 416 and 721 genes, respectively, were found to be significantly (p < 5 × 10^-4) correlated with 22 meat quality traits measured on longissiums dorsi muscle (LM) or semimembranosus muscle (SM). Transcripts associated with muscle glycolytic potential (GP) and pH values were enriched for genes involved in metabolic process. There were 42 QTTs (for 32 genes) shared by liver and muscle tissues, of which 10 QTTs represent GP- and/or pH-related genes, such as JUNB, ATF3, and PPP1R3B. Furthermore, a genome-wide eQTL mapping revealed a total of 3,054 eQTLs for all annotated transcripts in muscle (p < 2.08 × 10^-5), including 1,283 cis-eQTLs and 1771 trans-eQTLs. In addition, WGCNA identified five modules relevant to glycogen metabolism pathway and highlighted the connections between variations in meat quality traits and genes involved in energy process. Integrative analysis of GWAS loci, eQTL, and QTT demonstrated GALNT15/GALNTL2 and HTATIP2 as strong candidate genes for drip loss and pH drop from postmortem 45 min to 24 h, respectively. Our findings provide valuable insights into the genetic basis of meat quality traits and greatly expand the number of candidate genes that may be valuable for future functional analysis and genetic improvement of meat quality.

Detection of whole genome selection signatures of Pakistani Teddy goat

BACKGROUND

Natural and artificial selection tend to cause variability that contributes to shape the genome of livestock in a way that differentiates them among the animal kingdom. The particular aim here is to identify positive selection signatures with whole genome pooled-sequence data of Pakistani Teddy goat.

METHODS AND RESULTS

Paired-end alignment of 635,357,043 reads of Teddy goat with (ARS1) reference genome assembly was carried out. Pooled-Heterozygosity (Hp) and Tajima's D (TD) are applied for validation and getting better hits of selection signals, while pairwise F_ST statistics is conducted on Teddy vs. Bezoar (wild goat ancestor) for genomic differentiation, moreover annotation of regions under positive selection was also performed. Hp score with - ZHp > 5 detected six windows having highest hits on Chr. 29, 9, 25, 15 and 14 that harbor HRASLS5, LACE1 and AXIN1 genes which are candidate for embryonic development, lactation and body height. Secondly, - ZTD value of > 3.3 showed 4 windows with very strong hits on Chr.5 & 9 which harbor STIM1 and ADM genes related to body mass and weight. Lastly, - ZF_ST < - 5 generated four strong signals on Chr.5 & 12 harbor LOC102183233 gene. Other significant selection signatures encompass genes associated with wool production, prolificacy and coat colors traits in this breed.

CONCLUSIONS

In brief, this study identified the genes under selection in Pakistani Teddy goat that will be helpful to refining the marker-assisted breeding policies and converging required production traits within and across other goat breeds and to explore full genetic potential of this valued species of livestock.

Meat tenderness: advances in biology, biochemistry, molecular mechanisms and new technologies

Meat tenderness is an important quality trait critical to consumer acceptance, and determines satisfaction, repeat purchase and willingness-to-pay premium prices. Recent advances in tenderness research from a variety of perspectives are presented. Our understanding of molecular factors influencing tenderization are discussed in relation to glycolysis, calcium release, protease activation, apoptosis and heat shock proteins, the use of proteomic analysis for monitoring changes, proteomic biomarkers and oxidative/nitrosative stress. Each of these structural, metabolic and molecular determinants of meat tenderness are then discussed in greater detail in relation to animal variation, postmortem influences, and changes during cooking, with a focus on recent advances. Innovations in postmortem technologies and enzymes for meat tenderization are discussed including their potential commercial application. Continued success of the meat industry relies on ongoing advances in our understanding, and in industry innovation. The recent advances in fundamental and applied research on meat tenderness in relation to the various sectors of the supply chain will enable such innovation.

Accurate prediction and genome-wide association analysis of digital intramuscular fat content in longissimus muscle of pigs

Intramuscular fat (IMF) content is a critical indicator of pork quality that affects directly the purchasing desire of consumers. However, to measure IMF content is both laborious and costly, preventing our understanding of its genetic determinants and improvement. In the present study, we constructed an accurate and fast image acquisition and analysis system, to extract and calculate the digital IMF content, the proportion of fat areas in the image (PFAI) of the longissimus muscle of 1709 animals from multiple pig populations. PFAI was highly significantly correlated with marbling scores (MS; 0.95, r²  = 0.90), and also with IMF contents chemically defined for 80 samples (0.79, r²  = 0.63; more accurate than direct analysis between IMF contents and MS). The processing time for one image is only 2.31 s. Genome-wide association analysis on PFAI for all 1709 animals identified 14 suggestive significant SNPs and 1 genome-wide significant SNP. On MS, we identified nine suggestive significant SNPs, and seven of them were also identified in PFAI. Furthermore, the significance (-log P) values of the seven common SNPs are higher in PFAI than in MS. Novel candidate genes of biological importance for IMF content were also discovered. Our imaging systems developed for prediction of digital IMF content is closer to IMF measured by Soxhlet extraction and slightly more accurate than MS. It can achieve fast and high-throughput IMF phenotype, which can be used in improvement of pork quality.

Genome-Wide Association Study of Meat Quality Traits in a Three-Way Crossbred Commercial Pig Population

Meat quality is an important trait for pig-breeding programs aiming to meet consumers' demands. Geneticists must improve meat quality based on their understanding of the underlying genetic mechanisms. Previous studies showed that most meat-quality indicators were low-to-moderate heritability traits; therefore, improving meat quality using conventional techniques remains a challenge. Here, we performed a genome-wide association study of meat-quality traits using the GeneSeek Porcine SNP50K BeadChip in 582 crossbred Duroc × (Landrace × Yorkshire) commercial pigs (249 males and 333 females). Meat conductivity, marbling score, moisture, meat color, pH, and intramuscular fat (IMF) content were investigated. The genome-wide association study was performed using both fixed and random model Circulating Probability Unification (FarmCPU) and a mixed linear model (MLM) with the rMVP software. The genomic heritability of the studied traits ranged from 0.13 ± 0.07 to 0.55 ± 0.08 for conductivity and meat color, respectively. Thirty-two single-nucleotide polymorphisms (SNPs) were identified for meat quality in the crossbred pigs using both FarmCPU and MLM. Among the detected SNPs, five, nine, seven, four, six, and five were significantly associated with conductivity, IMF, marbling score, meat color, moisture, and pH, respectively. Several candidate genes for meat quality were identified in the detected genomic regions. These findings will contribute to the ongoing improvement of meat quality, meeting consumer demands and improving the economic outlook for the swine industry.

Genetic parameter estimation and gene association analyses for meat quality traits in open-air free-range Iberian pigs

Meat quality of Iberian pigs is defined by the combination of their genetic characteristics and the particular production system. To carry out a genetic analysis of the main meat quality traits, we estimated their heritabilities, genetic correlations and the association effects of 32 selected SNPs of 12 candidate genes. A total of ten traits were measured in longissimus dorsi samples from 1,199 Iberian pigs fattened in the traditional free-range system: water holding capacity (thawing, cooking and centrifuge force water losses), instrumental colour (lightness L*, redness a* and yellowness b*), myoglobin content, shear force on cooked meat, and shear force and maximum compression force on dry-cured loin. Estimated heritability values were low to moderate (0.01 to 0.43) being the lowest for L* and the highest for cooking loss. Strong genetic correlations between water holding capacity traits (0.93 to 0.96) and between myoglobin content and a* (0.94) were observed. The association analyses revealed 19 SNPs significantly associated with different traits. Consistent and strong effects were observed between PRKAG3 SNPs (rs319678464G > C and rs330427832C > T), MYH3_rs81437544T > C, CASP3_rs319658214G > T and CTSL_rs332171512A > G and water losses. Also for CAPN1_rs81358667G > A and CASP3_rs319658214G > T and shear force. The SNPs mapping on PRKAG3 showed the highest effects on Minolta colour traits. Genotyping of these SNPs could be useful for the selection of Iberian young boars with similar estimated breeding values for productive traits.

Recent Updates in Meat Color Research: Integrating Traditional and High-Throughput Approaches

Deviation from a bright cherry-red color of fresh meat results in less consumer acceptance and either discounted or discarded products in the value chain. Tissue homeostasis changes immediately after exsanguination, leading to acidification of muscle. Any alteration in pH drop can influence both muscle structure and enzymatic activity related to oxygen consumption and the redox state of myoglobin. This review focuses on both fundamental and applied approaches to under-stand the effects of pH on biochemical changes, oxygen diffusion, and its impact on meat color. Recent updates utilizing high-throughput “omics” approaches to elucidate the biochemical changes associated with high-pH meat are also dis-cussed. The fundamental aspects affecting fresh meat color are complex and highly interrelated with factors ranging from live animal production to preharvest environmental issues, muscle to meat conversion, and numerous facets along the merchandising chain of marketing meat to consumers.

Whole-genome sequencing association analysis reveals the genetic architecture of meat quality traits in Chinese Qingyu pigs

The Chinese Qingyu pig breed is an invaluable indigenous genetic resource. However, few studies have investigated the genetic architecture of meat quality traits in Qingyu pigs. Here, 30 purebred Qingyu pigs were subjected to whole-genome sequencing. After quality control, 18 436 759 SNPs were retained. Genome-wide association studies (GWAS) were then performed for meat pH and color at three postmortem time points (45 min, 24 h, and 48 h) using single-marker regression analysis. In total, 11 and 69 SNPs were associated with meat pH and color of the longissimus thoracis muscle (LTM), respectively, while 54 and 29 SNPs were associated with meat pH and color of the semimembranosus muscle (SM), respectively. Seven SNPs associated with pork pH were shared by all three postmortem time points. Several candidate genes for meat traits were identified, including four genes (CXXC5, RYR3, BNIP3, and MYCT1) related to skeletal muscle development, regulation of Ca²⁺ release in the muscle, and anaerobic respiration, which are promising candidates for selecting superior meat quality traits in Qingyu pigs. To our knowledge, this is the first study investigating the postmortem genetic architecture of pork pH and color in Qingyu pigs. Our findings further the current understanding of the genetic factors influencing meat quality.

Genomics of Heat Tolerance in Reproductive Performance Investigated in Four Independent Maternal Lines of Pigs

Improving swine climatic resilience through genomic selection has the potential to minimize welfare issues and increase the industry profitability. The main objective of this study was to investigate the genetic and genomic determinism of tolerance to heat stress in four independent purebred populations of swine. Three female reproductive traits were investigated: total number of piglets born (TNB), number of piglets born alive (NBA) and average birth weight (ABW). More than 80,000 phenotypic and 12,000 genotyped individuals were included in this study. Genomic random-regression models were fitted regressing the phenotypes of interest on a set of 95 environmental covariates extracted from public weather station records. The models yielded estimates of (genomic) reactions norms for individual pigs, as indicator of heat tolerance. Heat tolerance is a heritable trait, although the heritabilities are larger under comfortable than heat-stress conditions (larger than 0.05 vs. 0.02 for TNB; 0.10 vs. 0.05 for NBA; larger than 0.20 vs. 0.10 for ABW). TNB showed the lowest genetic correlation (-38%) between divergent climatic conditions, being the trait with the strongest impact of genotype by environment interaction, while NBA and ABW showed values slightly negative or equal to zero reporting a milder impact of the genotype by environment interaction. After estimating genetic parameters, a genome-wide association study was performed based on the single-step GBLUP method. Heat tolerance was observed to be a highly polygenic trait. Multiple and non-overlapping genomic regions were identified for each trait based on the genomic breeding values for reproductive performance under comfortable or heat stress conditions. Relevant regions were found on chromosomes (SSC) 1, 3, 5, 6, 9, 11, and 12, although there were important regions across all autosomal chromosomes. The genomic region located on SSC9 appears to be of particular interest since it was identified for two traits (TNB and NBA) and in two independent populations. Heat tolerance based on reproductive performance indicators is a heritable trait and genetic progress for heat tolerance can be achieved through genetic or genomic selection. Various genomic regions and candidate genes with important biological functions were identified, which will be of great value for future functional genomic studies.

Muscle glycogen level and occurrence of acid meat in commercial hybrid pigs are regulated by two low-frequency causal variants with large effects and multiple common variants with small effects

Background

Meat production from the commercial crossbred Duroc × (Landrace × Yorkshire) (DLY) pig is predominant in the pork industry, but its meat quality is often impaired by low ultimate pH (pHu). Muscle glycogen level at slaughter is closely associated with pHu and meat technological quality, but its genetic basis remains elusive. The aim of this study was to identify genes and/or causative mutations associated with muscle glycogen level and other meat quality traits by performing a genome-wide association study (GWAS) and additional analyses in a population of 610 DLY pigs.

Results

Our initial GWAS identified a genome-wide significant (P = 2.54e−11) quantitative trait locus (QTL) on SSC15 (SSC for Sus scrofa chromosome) for the level of residual glycogen and glucose (RG) in the longissimus muscle at 45 min post-mortem. Then, we demonstrated that a low-frequency (minor allele frequency = 0.014) R200Q missense mutation in the PRKAG3 (RN) gene caused this major QTL effect on RG. Moreover, we showed that the 200Q (RN^–) allele was introgressed from the Hampshire breed into more than one of the parental breeds of the DLY pigs. After conditioning on R200Q, re-association analysis revealed three additional QTL for RG on SSC3 and 4, and on an unmapped scaffold (AEMK02000452.1). The SSC3 QTL was most likely caused by a splice mutation (g.8283C>A) in the PHKG1 gene that we had previously identified. Based on functional annotation, the genes TMCO1 on SSC4 and CKB on the scaffold represent promising candidate genes for the other two QTL. There were significant interaction effects of the GWAS tag SNPs at those two loci with PRKAG3 R200Q on RG. In addition, a number of common variants with potentially smaller effects on RG (P < 10⁻⁴) were uncovered by a second conditional GWAS after adjusting for the two causal SNPs, R200Q and g.8283C>A.

Conclusions

We found that the RN^– allele segregates in the parental lines of our DLY population and strongly influences its meat quality. Our findings also indicate that the genetic basis of RG in DLY can be mainly attributed to two major genes (PRKAG3 and PHKG1), along with many minor genes.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0488-0) contains supplementary material, which is available to authorized users.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

Data Mining and Validation of AMPK Pathway as a Novel Candidate Role Affecting Intramuscular Fat Content in Pigs

Simple Summary

Intramuscular fat (IMF) is increasingly being recognized as a key meat trait in the modern pork industry. The aims of this research were to identify potential signaling pathways associated with IMF content in the longissimus dorsi (LD) muscle of different pig breeds and investigate the gene expression levels in the screened signaling pathways. Our results indicated that the AMPK signaling pathway may be related to IMF deposition in the LD muscle of pigs. The results of qRT-PCR analysis showed that the expression of ten key hub genes (AMPK, ADIPOR1, ADIPOR2, LKB1, CAMKKβ, CPT1A, CPT1B, PGC-1α, CD36, and ACC1) differed between the LD muscle of Min and Large White pigs. The protein expression levels of AMPK, LKB1, CaMKK2, CPT1A, and ACC1 were similar to the genes expression patterns in the LD muscle of Large White pigs. The results of this study provide novel insights into the regulatory function of the AMPK signaling pathway in relation to the IMF content in the LD muscle of different pigs.

Abstract

Intramuscular fat (IMF) is an important economic trait for pork quality and a complex quantitative trait regulated by multiple genes. The objective of this work was to investigate the novel transcriptional effects of a multigene pathway on IMF deposition in the longissimus dorsi (LD) muscles of pigs. Potential signaling pathways were screened by mining data from three gene expression profiles in the Gene Expression Omnibus (GEO) database. We designed quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) arrays for the candidate signaling pathways to verify the results in the LD muscles of two pig breeds with different IMF contents (Large White and Min). Western blot analysis was used to detect the expression levels of several candidate proteins. Our results showed that the AMPK signaling pathway was screened via bioinformatics analysis. Ten key hub genes of this signaling pathway (AMPK, ADIPOR1, ADIPOR2, LKB1, CAMKKβ, CPT1A, CPT1B, PGC-1α, CD36, and ACC1) were differentially expressed between the Large White and Min pigs. Western blot analysis further confirmed that LKB1/CaMKK2-AMPK-ACC1-CPT1A axis dominates the activity of AMPK signaling pathway. Statistical analyses revealed that AMPK signaling pathway activity clearly varied among the two pig breeds. Based on these results, we concluded that the activation of the AMPK signaling pathway plays a positive role in reducing IMF deposition in pigs.

Comprehensive Analysis of Differentially Expressed mRNA, lncRNA and circRNA and Their ceRNA Networks in the Longissimus Dorsi Muscle of Two Different Pig Breeds

Circular RNA (circRNA) and long non-coding RNA (lncRNA) are known to participate in adipogenesis and myogenic differentiation, but their impact on porcine muscle traits is not well understood. We compared their expressional profiles in the longissimus dorsi muscle of Chinese Huainan pigs (HN, the fat type) and Western commercial Duroc×(Landrace×Yorkshire) (DLY, the thin type) pigs, and 854 mRNAs, 233 lncRNAs, and 66 circRNAs (p < 0.05 and ｜log₂FoldChange｜>1) were found to be differentially expressed. The differentially expressed mRNA and circRNA parental genes were enriched in the Wnt signaling pathway (adipogenesis), the transition between fast and slow fibers (myogenic differentiation), and alanine, aspartate and glutamate metabolism (pork flavor). The potential lncRNAs/circRNAs-miRNAs-mRNAs regulatory networks shared MYOD1, PPARD, miR-423-5p and miR-874, which were associated with skeletal muscle muscular proliferation, differentiation/regeneration and adipogenesis. Taken together, these differentially expressed non-coding RNAs may be involved in the molecular basis of muscle traits, acting as the competing endogenous RNA (ceRNA) for miRNAs.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs

Background

Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs.

Results

Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat.

Conclusion

Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5386-2) contains supplementary material, which is available to authorized users.

Genome-wide association study for intramuscular fat content in Chinese Lulai black pigs

Objective

Intramuscular fat (IMF) content plays an important role in meat quality. Identification of single nucleotide polymorphisms (SNPs) and genes related to pig IMF, especially using pig populations with high IMF content variation, can help to establish novel molecular breeding tools for optimizing IMF in pork and unveil the mechanisms that underlie fat metabolism.

Methods

We collected muscle samples of 453 Chinese Lulai black pigs, measured IMF content by Soxhlet petroleum-ether extraction method, and genotyped genome-wide SNPs using GeneSeek Genomic Profiler Porcine HD BeadChip. Then a genome-wide association study was performed using a linear mixed model implemented in the GEMMA software.

Results

A total of 43 SNPs were identified to be significantly associated with IMF content by the cutoff p<0.001. Among these significant SNPs, the greatest number of SNPs (n = 19) were detected on Chr.9, and two linkage disequilibrium blocks were formed among them. Additionally, 17 significant SNPs are mapped to previously reported quantitative trait loci (QTLs) of IMF and confirmed previous QTLs studies. Forty-two annotated genes centering these significant SNPs were obtained from Ensembl database. Overrepresentation test of pathways and gene ontology (GO) terms revealed some enriched reactome pathways and GO terms, which mainly involved regulation of basic material transport, energy metabolic process and signaling pathway.

Conclusion

These findings improve our understanding of the genetic architecture of IMF content in pork and facilitate the follow-up study of fine-mapping genes that influence fat deposition in muscle.

Gene Co-expression Analysis Indicates Potential Pathways and Regulators of Beef Tenderness in Nellore Cattle

Beef tenderness, a complex trait affected by many factors, is economically important to beef quality, industry, and consumer’s palatability. In this study, RNA-Seq was used in network analysis to better understand the biological processes that lead to differences in beef tenderness. Skeletal muscle transcriptional profiles from 24 Nellore steers, selected by extreme estimated breeding values (EBVs) for shear force after 14 days of aging, were analyzed and 22 differentially expressed transcripts were identified. Among these were genes encoding ribosomal proteins, glutathione transporter ATP-binding cassette, sub-family C (CFTR/MRP), member 4 (ABCC4), and synaptotagmin IV (SYT4). Complementary co-expression analyses using Partial Correlation with Information Theory (PCIT), Phenotypic Impact Factor (PIF) and the Regulatory Impact Factor (RIF) methods identified candidate regulators and related pathways. The PCIT analysis identified ubiquitin specific peptidase 2 (USP2), growth factor receptor-bound protein 10 (GBR10), anoctamin 1 (ANO1), and transmembrane BAX inhibitor motif containing 4 (TMBIM4) as the most differentially hubbed (DH) transcripts. The transcripts that had a significant correlation with USP2, GBR10, ANO1, and TMBIM4 enriched for proteasome KEGG pathway. RIF analysis identified microRNAs as candidate regulators of variation in tenderness, including bta-mir-133a-2 and bta-mir-22. Both microRNAs have target genes present in the calcium signaling pathway and apoptosis. PIF analysis identified myoglobin (MB), enolase 3 (ENO3), and carbonic anhydrase 3 (CA3) as potentially having fundamental roles in tenderness. Pathways identified in our study impacted in beef tenderness included: calcium signaling, apoptosis, and proteolysis. These findings underscore some of the complex molecular mechanisms that control beef tenderness in Nellore cattle.

A whole-genome sequence based association study on pork eating quality traits and cooking loss in a specially designed heterogeneous F6 pig population

To determine the genetic basis of pork eating quality traits and cooking loss, we herein performed a genome-wide association study (GWAS) for tenderness, juiciness, oiliness, umami, overall liking and cooking loss by using whole genome sequences of heterogeneous stock F6 pigs which were generated by crossing 4 typical western pig breeds (Duroc, Landrace, Large White and Pietrain) and 4 typical Asian pig breeds (Erhualian, Laiwu, Bamaxiang and Tibetan). We identified 50 associated loci (QTLs) and most of them are novel. Seven loci also showed pleiotropic associations with different traits. In addition, we identified multiple promising candidate genes for these traits, including PAK1 and AQP11 for cooking loss, EP300 for tenderness, SDK1 for juiciness, FITM2 and 5-linked MYH genes for oiliness, and TNNI2 and TNNT3 for overall liking. Our results provide not only a better understanding of the genetic basis for meat quality, but also a potential application in future breeding for these complex traits.

Genome-Wide Association of Myoglobin Concentrations in Pork Loins

Lean color is a major focus for identifying pork loins for export markets, and myoglobin is the primary pigment driving pork color. Thus, increasing myoglobin concentration should increase redness of pork products and the number of loins acceptable for exportation. Therefore, understanding genetic variation and parameters affecting myoglobin concentration is critical for improving pork color. The objective of this study was to identify genetic markers associated with myoglobin concentration in pork loin muscle. Ultimate pH and myoglobin concentrations were measured in longissimus thoracis et lumborum samples of pigs (n = 599) from two different commercial finishing swine facilities. A Bayes-C model implemented in GenSel identified regions within 7 chromosomes that explained greater than 63% of the genetic variance in myoglobin concentration. Chromosome 7 had 1 significant region which accounted for 37% of the genetic variance, while chromosome 14 had 4 significant regions accounting for 9.8% of the genetic variance. Candidate genes in the region on chromosome 7 were involved in iron homeostasis, and genes in the significant regions on chromosome 14 were involved in calcium regulation. Genes identified in this study represent potential biomarkers that could be used to select for higher myoglobin concentrations in pork, which may improve lean meat color.

Genome-wide association studies, meta-analyses and derived gene network for meat quality and carcass traits in pigs

A large number of quantitative trait loci (QTL) for meat quality and carcass traits has been reported in pigs over the past 20 years. However, few QTL have been validated and the biological meaning of the genes associated to these QTL has been underexploited. In this context, a meta-analysis was performed to compare the significant markers with meta-QTL previously reported in literature. Genome association studies were performed for 12 traits, from which 144 SNPs were found out to be significant (P < 0.05). They were validated in the meta-analysis and used to build the Association Weight Matrix, a matrix framework employed to investigate co-association of pairwise SNP across phenotypes enabling to derive a gene network. A total of 45 genes were selected from the Association Weight Matrix analysis, from which 25 significant transcription factors were identified and used to construct the networks associated to meat quality and carcass traits. These networks allowed the identification of key transcription factors, such as SOX5 and NKX2–5, gene–gene interactions (e.g. ATP5A1, JPH1, DPT and NEDD4) and pathways related to the regulation of adipose tissue metabolism and skeletal muscle development. Validated SNPs and knowledge of key genes driving these important industry traits might assist future strategies in pig breeding.

Genome-wide association study in an F2 Duroc x Pietrain resource population for economically important meat quality and carcass traits

Meat quality is essential for consumer acceptance, it ultimately impacts pork production profitability and it is subject to genetic control. The objective of this study was to map genomic regions associated with economically important meat quality and carcass traits. We performed a genome-wide association (GWA) analysis to map regions associated with 38 meat quality and carcass traits recorded for 948 F2 pigs from the Michigan State University Duroc × Pietrain resource population. The F0, F1, and 336 F2 pigs were genotyped with the Illumina Porcine SNP60 BeadChip, while the remaining F2 pigs were genotyped with the GeneSeek Genomic Profiler for Porcine Low Desnisty (LD) chip, and imputed with high accuracy ( = 0.97). Altogether the genomic dataset comprised 1,019 animals and 44,911 SNP. A Gaussian linear mixed model was fitted to estimate the breeding values and the variance components. A linear transformation was performed to estimate the marker effects and variances. Type I error rate was controlled at a False Discovery Rate of 5%. Seven putative QTL found in this study were previously reported in other studies. Two novel QTL associated with tenderness (TEN) were located on SSC3 [135.6:137.5Mb; False Discovery rate (FDR) < 0.03] and SSC5 (67.3:69.1Mb; FDR < 0.02). The QTL region identified on SSC15 includes Protein Kinase AMP-activated ɣ 3-subunit gene (), which has been associated with 24-h pH (pH24), drip loss (DL) and cook yield (CY). Also, novel candidate genes were identified for TEN in the region on SSC5 [A Kinase (PRKA) Anchor Protein 3 (], and for tenth rib backfat thickness (BF10) [Carnitine O-Acetyltransferase ()] on SSC1. The association of gene polymorphisms with pork quality traits has been reported for several pig populations. However, there are no SNP for this gene on the chip used, thus we genotyped the animals for 2 non-synonymous variants ( and ). We then performed a GWA conditioning on the genotype of both SNP and was associated with pH24, DL, protein content (PRO) and CY ( < 0.004) and T30N with Juiciness, TEN, shear force, pH24, PRO, and CY < 0.04). Finally, we performed a GWA conditioning on the genotype of the SNP peak detected in this study, and T30N remained associated only with PRO ( < 0.02). Therefore, in this study we identified 2 novel QTL regions, suggest 2 novel candidate genes, and conclude that other SNP in PRKAG3 or nearby gene(s) explain the observed associations on SSC15 in this population.

Enhancement of porcine intramuscular fat content by overexpression of the cytosolic form of phosphoenolpyruvate carboxykinase in skeletal muscle

Intramuscular fat (IMF) content has been generally recognized as a desirable trait in pork meat because of its positive effect on eating quality. An effective approach to enhance IMF content in pork is the generation of transgenic pigs. In this study, we used somatic cell nuclear transfer (SCNT) to generate cloned pigs exhibiting ectopic expression of phosphoenolpyruvate carboxykinase (PEPCK-C) driven by an α-skeletal-actin gene promoter, which was specifically expressed in skeletal muscle. Using qRT-PCR and Western blot analysis, we demonstrated that PEPCK-C was functionally expressed and had a significant effect on total fatty acid content in the skeletal muscle of the transgenic pigs, while the n-6/n-3 polyunsaturated fatty acid (PUFA) ratio showed no difference between transgenic and control pigs. Thus, genetically engineered PEPCK-C^mus pigs may be an effective solution for the production of IMF-enriched pork.

Whole-genome association analysis of pork meat pH revealed three significant regions and several potential genes in Finnish Yorkshire pigs

Background

One of the most commonly used quality measurements of pork is pH measured 24 h after slaughter. The most probable mode of inheritance for this trait is oligogenic with several known major genes, such as PRKAG3. In this study, we used whole-genome SNP genotypes of over 700 AI boars; after a quality check, 42,385 SNPs remained for association analysis. All the boars were purebred Finnish Yorkshire. To account for relatedness of the animals, a pedigree-based relationship matrix was used in a mixed linear model to test the effect of SNPs on pH measured from loin. A bioinformatics analysis was performed to identify the most promising genes in the significant regions related to meat quality.

Results

Genome-wide association study (GWAS) revealed three significant chromosomal regions: one on chromosome 3 (39.9 Mb–40.1 Mb) and two on chromosome 15 (58.5 Mb–60.5 Mb and 132 Mb–135 Mb including PRKAG3). A conditional analysis with a significant SNP in the PRKAG3 region, MARC0083357, as a covariate in the model retained the significant SNPs on chromosome 3. Even though linkage disequilibrium was relatively high over a long distance between MARC0083357 and other significant SNPs on chromosome 15, some SNPs retained their significance in the conditional analysis, even in the vicinity of PRKAG3. The significant regions harbored several genes, including two genes involved in cyclic AMP (cAMP) signaling: ADCY9 and CREBBP. Based on functional and transcription factor-gene networks, the most promising candidate genes for meat pH are ADCY9, CREBBP, TRAP1, NRG1, PRKAG3, VIL1, TNS1, and IGFBP5, and the key transcription factors related to these genes are HNF4A, PPARG, and Nkx2-5.

Conclusions

Based on SNP association, pathway, and transcription factor analysis, we were able to identify several genes with potential to control muscle cell homeostasis and meat quality. The associated SNPs can be used in selection for better pork. We also showed that post-GWAS analysis reveals important information about the genes’ potential role on meat quality. The gained information can be used in later functional studies.

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH₂₄), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH₂₄ SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Identification of candidate genes associated with porcine meat color traits by genome-wide transcriptome analysis

Meat color is considered to be the most important indicator of meat quality, however, the molecular mechanisms underlying traits related to meat color remain mostly unknown. In this study, to elucidate the molecular basis of meat color, we constructed six cDNA libraries from biceps femoris (Bf) and soleus (Sol), which exhibit obvious differences in meat color, and analyzed the whole-transcriptome differences between Bf (white muscle) and Sol (red muscle) using high-throughput sequencing technology. Using DEseq2 method, we identified 138 differentially expressed genes (DEGs) between Bf and Sol. Using DEGseq method, we identified 770, 810, and 476 DEGs in comparisons between Bf and Sol in three separate animals. Of these DEGs, 52 were overlapping DEGs. Using these data, we determined the enriched GO terms, metabolic pathways and candidate genes associated with meat color traits. Additionally, we mapped 114 non-redundant DEGs to the meat color QTLs via a comparative analysis with the porcine quantitative trait loci (QTL) database. Overall, our data serve as a valuable resource for identifying genes whose functions are critical for meat color traits and can accelerate studies of the molecular mechanisms of meat color formation.

Identification and characterization of novel single nucleotide polymorphism markers for fat deposition in muscle tissue of pigs using amplified fragment length polymorphism

Objective

This study was conducted to identify and evaluate the effective single nucleotide polymorphism (SNP) markers for fat deposition in the longissimus dorsi muscles of pigs using the amplified fragment length polymorphism (AFLP) approach.

Methods

Sixty-four selective primer combinations were used to identify the AFLP markers in the 20 highest- and 20 lowest-intramuscular fat (IMF) content phenotypes. Five AFLP fragments were converted into simple codominant SNP markers. These SNP markers were tested in terms of their association with IMF content and fatty acid (FA) composition traits in 620 commercially crossbred pigs.

Results

The SSC7 g.4937240C>G marker showed an association with IMF content (p<0.05). The SSC9 g.5496647_5496662insdel marker showed a significant association with IMF content and arachidonic levels (p<0.05). The SSC10 g.71225134G>A marker revealed an association with palmitoleic and ω9 FA levels (p<0.05), while the SSC17 g.61976696G>T marker showed a significant association with IMF content and FA levels of palmitoleic, eicosenoic, arachidonic, monounsaturated fatty acids, and ω9 FA levels. However, no significant association of SSC8 g.47338181G>A was observed with any IMF and FA levels in this study.

Conclusion

Four SNP markers (SSC7 g.4937240C>G, SSC9 g.5496647_5496662insdel, SSC10 g.71225134G>A, and SSC17 g.61976696G>T) were found to be associated with IMF and/or FA content traits in commercially crossbred pigs. These findings provide evidence of the novel SNP markers as being potentially useful for selecting pigs with the desirable IMF content and FA composition.

Implementing meta-analysis from genome-wide association studies for pork quality traits

Pork quality plays an important role in the meat processing industry. Thus, different methodologies have been implemented to elucidate the genetic architecture of traits affecting meat quality. One of the most common and widely used approaches is to perform genome-wide association (GWA) studies. However, a limitation of many GWA in animal breeding is the limited power due to small sample sizes in animal populations. One alternative is to implement a meta-analysis of GWA (MA-GWA) combining results from independent association studies. The objective of this study was to identify significant genomic regions associated with meat quality traits by performing MA-GWA for 8 different traits in 3 independent pig populations. Results from MA-GWA were used to search for genes possibly associated with the set of evaluated traits. Data from 3 pig data sets (U.S. Meat Animal Research Center, commercial, and Michigan State University Pig Resource Population) were used. A MA was implemented by combining -scores derived for each SNP in every population and then weighting them using the inverse of estimated variance of SNP effects. A search for annotated genes retrieved genes previously reported as candidates for shear force (calpain-1 catalytic subunit [] and calpastatin []), as well as for ultimate pH, purge loss, and cook loss (protein kinase, AMP-activated, γ 3 noncatalytic subunit []). In addition, novel candidate genes were identified for intramuscular fat and cook loss (acyl-CoA synthetase family member 3 mitochondrial []) and for the objective measure of muscle redness, CIE a* (glycogen synthase 1, muscle [] and ferritin, light polypeptide []). Thus, implementation of MA-GWA allowed integration of results for economically relevant traits and identified novel genes to be tested as candidates for meat quality traits in pig populations.

Identification of signatures of selection for intramuscular fat and backfat thickness in two Duroc populations

Intramuscular fat (IMF) content is an important trait affecting the quality of pork. Two Duroc populations, one under positive selection for IMF and the other selected for decreased backfat but under stabilizing selection for IMF, were used to identify signatures of selection associated with IMF using 60,000 single-nucleotide polymorphism data. The effects of selection were analyzed between 2 lines or groups representing selected and control animals within each population using a discriminant analysis of principal components and Wright’s fixation index (FST). Moreover, extended haplotype homozygosity-based approaches were used to examine the changes in haplotype frequency due to recent selection. Each statistical method identified 10–20 selection signatures. A few haplotype-based signatures of selection agreed with results from a genome-wide association study (GWAS), while FST measures showed a better agreement with GWAS results. Agreement of marker-trait associations and signatures of selection was limited, and further examination will be necessary to understand the effect of selection on IMF and why some regions identified by GWAS did not appear to respond to the selection practiced. The genes in 21 consensus selection signatures were examined. Several genes with an effect on overall fatness were identified, but further research is needed to assess whether or not some of them could have a specific effect on IMF.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

A gene and protein expression study on four porcine genes related to intramuscular fat deposition

Intramuscular fat (IMF) content has a prominent role in meat quality, affecting sensory attributes such as flavour and texture. In the present research, we studied in samples of porcine Semimembranosus muscle four genes related to lipid metabolism and whose gene expressions have been associated to IMF deposition: FASN, SCD, LIPE and LPL. We analysed both mRNA and protein expressions in two groups of Italian Large White pigs divergent for Semimembranosus IMF deposition, with the aim of comparing the levels of four genes and enzymes between the two groups and identifying possible coexpression links. The obtained results suggest a prominent role of LIPE enzyme in IMF hydrolysis, as the samples with low IMF deposition show a significantly higher amount of this lipase. Finally, a poorly known correlation was found between LIPE and FASN enzymes only in female individuals. These results provide new information for the understanding of IMF deposition.

Genome Wide Association Studies (GWAS) Identify QTL on SSC2 and SSC17 Affecting Loin Peak Shear Force in Crossbred Commercial Pigs

Of all the meat quality traits, tenderness is considered the most important with regard to eating quality and market value. In this study we have utilised genome wide association studies (GWAS) for peak shear force (PSF) of loin muscle as a measure of tenderness for 1,976 crossbred commercial pigs, genotyped for 42,721 informative SNPs using the Illumina PorcineSNP60 Beadchip. Four 1 Mb genomic regions, three on SSC2 (at 4 Mb, 5 Mb and 109 Mb) and one on SSC17 (at 20 Mb), were detected which collectively explained about 15.30% and 3.07% of the total genetic and phenotypic variance for PSF respectively. Markers ASGA0008566, ASGA0008695, DRGA0003285 and ASGA0075615 in the four regions were strongly associated with the effects. Analysis of the reference genome sequence in the region with the most important SNPs for SSC2_5 identified FRMD8, SLC25A45 and LTBP3 as potential candidate genes for meat tenderness on the basis of functional annotation of these genes. The region SSC2_109 was close to a previously reported candidate gene CAST; however, the very weak LD between DRGA0003285 (the best marker representing region SSC2_109) and CAST indicated the potential for additional genes which are distinct from, or interact with, CAST to affect meat tenderness. Limited information of known genes in regions SSC2_109 and SSC17_20 restricts further analysis. Re-sequencing of these regions for informative animals may help to resolve the molecular architecture and identify new candidate genes and causative mutations affecting this trait. These findings contribute significantly to our knowledge of the genomic regions affecting pork shear force and will potentially lead to new insights into the molecular mechanisms regulating meat tenderness.

Meta-analysis of genome-wide association from genomic prediction models

Genome-wide association (GWA) studies based on GBLUP models are a common practice in animal breeding. However, effect sizes of GWA tests are small, requiring larger sample sizes to enhance power of detection of rare variants. Because of difficulties in increasing sample size in animal populations, one alternative is to implement a meta-analysis (MA), combining information and results from independent GWA studies. Although this methodology has been used widely in human genetics, implementation in animal breeding has been limited. Thus, we present methods to implement a MA of GWA, describing the proper approach to compute weights derived from multiple genomic evaluations based on animal-centric GBLUP models. Application to real datasets shows that MA increases power of detection of associations in comparison with population-level GWA, allowing for population structure and heterogeneity of variance components across populations to be accounted for. Another advantage of MA is that it does not require access to genotype data that is required for a joint analysis. Scripts related to the implementation of this approach, which consider the strength of association as well as the sign, are distributed and thus account for heterogeneity in association phase between QTL and SNPs. Thus, MA of GWA is an attractive alternative to summarizing results from multiple genomic studies, avoiding restrictions with genotype data sharing, definition of fixed effects and different scales of measurement of evaluated traits.

Copy number variation-based genome wide association study reveals additional variants contributing to meat quality in Swine

Pork quality is important both to the meat processing industry and consumers' purchasing attitude. Copy number variation (CNV) is a burgeoning kind of variants that may influence meat quality. In this study, a genome-wide association study (GWAS) was performed between CNVs and meat quality traits in swine. After false discovery rate (FDR) correction, a total of 8 CNVs on 6 chromosomes were identified to be significantly associated with at least one meat quality trait. All of the 8 CNVs were verified by next generation sequencing and six of them were verified by qPCR. Only the haplotype block containing CNV12 is adjacent to significant SNPs associated with meat quality, suggesting the effects of those CNVs were not likely captured by tag SNPs. The DNA dosage and EST expression of CNV12, which overlap with an obesity related gene Netrin-1 (Ntn1), were consistent with Ntn1 RNA expression, suggesting the CNV12 might be involved in the expression regulation of Ntn1 and finally influence meat quality. We concluded that CNVs may contribute to the genetic variations of meat quality beyond SNPs, and several candidate CNVs were worth further exploration.

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) commercial population

BACKGROUND

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F2 crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

METHODS

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F2 and Sutai pig populations.

RESULTS

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

CONCLUSIONS

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

Genome-wide association studies (GWAS) identify a QTL close to PRKAG3 affecting meat pH and colour in crossbred commercial pigs

Background

Improving meat quality is a high priority for the pork industry to satisfy consumers’ preferences. GWAS have become a state-of-the-art approach to genetically improve economically important traits. However, GWAS focused on pork quality are still relatively rare.

Results

Six genomic regions were shown to affect loin pH and Minolta colour a* and b* on both loin and ham through GWAS in 1943 crossbred commercial pigs. Five of them, located on Sus scrofa chromosome (SSC) 1, SSC5, SSC9, SSC16 and SSCX, were associated with meat colour. However, the most promising region was detected on SSC15 spanning 133–134 Mb which explained 3.51% - 17.06% of genetic variance for five measurements of pH and colour. Three SNPs (ASGA0070625, MARC0083357 and MARC0039273) in very strong LD were considered most likely to account for the effects in this region. ASGA0070625 is located in intron 2 of ZNF142, and the other two markers are close to PRKAG3, STK36, TTLL7 and CDK5R2. After fitting MARC0083357 (the closest SNP to PRKAG3) as a fixed factor, six SNPs still remained significant for at least one trait. Four of them are intragenic with ARPC2, TMBIM1, NRAMP1 and VIL1, while the remaining two are close to RUFY4 and CDK5R2. The gene network constructed demonstrated strong connections of these genes with two major hubs of PRKAG3 and UBC in the super-pathways of cell-to-cell signaling and interaction, cellular function and maintenance. All these pathways play important roles in maintaining the integral architecture and functionality of muscle cells facing the dramatic changes that occur after exsanguination, which is in agreement with the GWAS results found in this study.

Conclusions

There may be other markers and/or genes in this region besides PRKAG3 that have an important effect on pH and colour. The potential markers and their interactions with PRKAG3 require further investigation

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0192-1) contains supplementary material, which is available to authorized users.

Genome-wide association analysis reveals genetic loci and candidate genes for meat quality traits in Chinese Laiwu pigs

Meat quality traits have economically significant impacts on the pig industry, and can be improved using molecular approaches in pig breeding. Since 1994 when the first genome-wide scan for quantitative trait loci (QTLs) in pig was reported, over the past two decades, numerous QTLs have been identified for meat quality traits by family based linkage analyses. However, little is known about the genetic variants for meat quality traits in Chinese purebred or outbred populations. To unveil it, we performed a genome-wide association study for 10 meat quality traits in Chinese purebred Laiwu pigs. In total, 75 significant SNPs (P < 1.01 × 10(-6)) and 33 suggestive SNPs (P < 2.03 × 10(-5)) were identified. On SSC12, a region between 56.22 and 61.49 Mb harbored a cluster of SNPs that were associated with meat color parameters (L*, lightness; a*, redness; b*, yellowness) and moisture content of longissimus muscle (LM) and semimembranosus muscle at the genome-wide significance level. A region on SSC4 also has pleiotropic effects on moisture content and drip loss of LM. In addition, this study revealed at least five novel QTLs and several candidate genes including 4-linked MYH genes (MYH1, MYH2, MYH3, and MYH13), MAL2, LPAR1, and PRKAG3 at four significant loci. Except for the SSC12 QTL, other QTLs are likely tissue-specific. These results provide new insights into the genetic basis of meat quality traits in Chinese Laiwu pigs and some significant SNPs reported here could be incorporated into the selection programs involving this breed.

Identification of novel putative adipomyokines by a cross-species annotation of secretomes and expression profiles

Adipose tissue and skeletal muscle are organs that respond strongly to obesity and physical activity exhibiting high secretory activity. To identify novel putative adipomyokines, comparative expression studies of skeletal muscle and adipose tissue of lean (C57BL/6J) and obese (C57BL/6J on a high-fat diet and NZO) mice, of sedentary and endurance trained C57BL/6J mice and of cattle characterized by different amounts of intramuscular fat were combined with human secretome data and scored. In highly regulated transcripts, we identified 119 myokines, 79 adipokines and 22 adipomyokines. Network analysis of these candidates revealed remodelling of extracellular matrix and tissue fibrosis as relevant functions of several of these candidates. Given the pathophysiogical relevance of fibrosis for adipose-muscle-cross-talk in obesity and type 2 diabetes and its physiological role in exercise adaptation and meat quality of farm animals, they represent interesting candidates for further investigations in different research areas and species.

Expression of variant transcripts of the potassium channel tetramerization domain-containing 15 (KCTD15) gene and their association with fatness traits in chickens

The aim of this study was to characterize the structure, expression, and biological functions of potassium channel tetramerization domain containing 15 (KCTD15) in chickens. We compared the KCTD15 expression level in samples of hypothalamic, adipose, and liver tissue of Xinghua chickens that were maintained on different dietary status. An association analysis of KCTD15 gene variant transcripts with fatness traits in a F2 resource population of chickens was performed. Three KCTD15 transcripts were identified in which the complete transcript was predominantly expressed in adipose tissue and the hypothalamus. The chicken KCTD15 gene was regulated by both feeding and fasting and consumption of a high-fat diet. The expression level of KCTD15 gene was markedly decreased in hypothalamus and liver of fasted and refed chickens (P < 0.05) and significantly downregulated in adipose tissue by the high-fat diet (P < 0.05). Three single-nucleotide polymorphisms of the KCTD15 gene were significantly associated with a number of fatness traits in chicken (P < 0.05). These results suggest that KCTD15 have a potential role regulation of obesity and fat metabolism in chickens.

Genome-wide association study for wool production traits in a Chinese Merino sheep population

Genome-wide association studies (GWAS) provide a powerful approach for identifying quantitative trait loci without prior knowledge of location or function. To identify loci associated with wool production traits, we performed a genome-wide association study on a total of 765 Chinese Merino sheep (JunKen type) genotyped with 50 K single nucleotide polymorphisms (SNPs). In the present study, five wool production traits were examined: fiber diameter, fiber diameter coefficient of variation, fineness dispersion, staple length and crimp. We detected 28 genome-wide significant SNPs for fiber diameter, fiber diameter coefficient of variation, fineness dispersion, and crimp trait in the Chinese Merino sheep. About 43% of the significant SNP markers were located within known or predicted genes, including YWHAZ, KRTCAP3, TSPEAR, PIK3R4, KIF16B, PTPN3, GPRC5A, DDX47, TCF9, TPTE2, EPHA5 and NBEA genes. Our results not only confirm the results of previous reports, but also provide a suite of novel SNP markers and candidate genes associated with wool traits. Our findings will be useful for exploring the genetic control of wool traits in sheep.