Differentially-expressed genes in pig Longissimus muscles with contrasting levels of fat, as identified by combined transcriptomic, reverse transcription PCR, and proteomic analyses

Comparative proteomic analyses to investigate premature browning in high‑oxygen modified atmosphere packaged beef patties

This study compared the proteomics of beef patties under high‑oxygen modified atmosphere packaging (HiOx-MAP) and vacuum packaging (VP) during heating. The color and oxidation stability of fresh patties, and myoglobin denaturation of cooked patties were also measured. The results suggested that HiOx-MAP patties contained more oxymyoglobin in fresh meat and had higher myoglobin denaturation during heating than VP patties, resulting in premature browning (PMB) during cooking. Proteomic analysis found that the overabundance of proteasome subunit beta type-2 (PSMB2) and peroxiredoxin-2 (PRDX2) in HiOx-55 °C, which can remove the damaged proteins and inhibit oxidation respectively, are of benefit to meat color stability during storage, however, this was still insufficient to inhibit the occurrence of PMB during cooking. The high abundance of lamin B1 (LMNB1) in VP-55 °C can maintain the stability of meat color. This research provides greater understanding, based on proteomic perspectives, of the molecular mechanism of PMB.

Identification of several lncRNA-mRNA pairs associated with marbling trait between Nanyang and Angus cattle

BACKGROUND

The marbling trait of cattle muscles, being a key indicator, played an important role in evaluating beef quality. Two breeds of cattle, namely a high-marbling (Angus) and a low-marbling (Nanyang) one, with their cattle muscles selected as our samples for transcriptome sequencing, were aimed to identify differentially expressed long non-coding RNAs (lncRNAs) and their targets associated with the marbling trait.

RESULTS

Transcriptome sequencing identified 487 and 283 differentially expressed mRNAs and lncRNAs respectively between the high-marbling (Angus) and low-marbling (Nanyang) cattle muscles. Twenty-seven pairs of differentially expressed lncRNAs-mRNAs, including eighteen lncRNAs and eleven target genes, were found to be involved in fat deposition and lipid metabolism. We established a positive correlation between fourteen up-regulated (NONBTAT000849.2, MSTRG.9591.1, NONBTAT031089.1, MSTRG.3720.1, NONBTAT029718.1, NONBTAT004228.2, NONBTAT007494.2, NONBTAT011094.2, NONBTAT015080.2, NONBTAT030943.1, NONBTAT021005.2, NONBTAT021004.2, NONBTAT025985.2, and NONBTAT023845.2) and four down-regulated (NONBTAT000850.2, MSTRG.22188.3, MSTRG.22188.4, and MSTRG.22188.5) lncRNAs and eleven genes related to adiponectin family protein (ADIPOQ), cytochrome P450 family (CYP4V2), 3-hydroxyacyl-CoA dehydratase family (HACD4), kinesin family (KIF5C), lipin family (LPIN2), perilipin family (PLIN1), prostaglandin family (PTGIS), solute carrier family (SLC16A7, SLC2213, and SLCO4C1), and containing a transmembrane domain protein family (VSTM1).

CONCLUSIONS

These candidate genes and lncRNAs can be regarded as being responsible for regulating the marbling trait of cattle. lncRNAs along with the variations in intramuscular fat marbling established a foundation for elucidating the genetic basis of high marbling in cattle.

Genome-Wide Association Study of Maternal Genetic Effects on Intramuscular Fat and Fatty Acid Composition in Rabbits

Maternal genetic effects (MGE) could affect meat quality traits such as intramuscular fat (IMF) and its fatty acid composition. However, it has been scarcely studied, especially in rabbits. The objectives of the present study were, first, to assess the importance of MGE on intramuscular fat and fatty acid composition by applying a Bayesian maternal animal model in two rabbit lines divergently selected for IMF. The second objective was to identify genomic regions and candidate genes of MGE that are associated with the traits of these offspring, using Bayesian methods in a Genome Wide Association Study (GWAS). Quantitative analyses were performed using data from 1982 rabbits, and 349 animals from the 9th generation and 76 dams of the 8th generation with 88,512 SNPs were used for the GWAS. The studied traits were IMF, saturated fatty acids (total SFA, C14:0; myristic acid, C16:0; palmitic acid and C18:0; stearic acid), monounsaturated fatty acids (total MUFA, C16:1n-7; palmitoleic acid and C18:1n-9; oleic acid), polyunsaturated fatty acids (total PUFA, C18:2n-6; linoleic acid, C18:3n-3; α-linolenic acid and C20:4n-6; arachidonic acid), MUFA/SFA and PUFA/SFA. The proportion of phenotypic variance explained by the maternal genetic effect ranged from 8 to 22% for IMF, depending on the model. For fatty acid composition, the proportion of phenotypic variance explained by maternal genetic effects varied from 10% (C18:0) to 46% (MUFA) in a model including both direct and additive maternal genetic effects, together with the common litter effect as a random variable. In particular, there were significant direct maternal genetic correlations for C16:0, C18:1n9, C18:2n6, SFA, MUFA, and PUFA with values ranging from -0.53 to -0.89. Relevant associated genomic regions were located on the rabbit chromosomes (OCU) OCU1, OCU5 and OCU19 containing some relevant candidates (TANC2, ACE, MAP3K3, TEX2, PRKCA, SH3GL2, CNTLN, RPGRIP1L and FTO) related to lipid metabolism, binding, and obesity. These regions explained about 1.2 to 13.9% of the total genomic variance of the traits studied. Our results showed an important maternal genetic effect on IMF and its fatty acid composition in rabbits and identified promising candidate genes associated with these traits.

Comparative Transcriptomic Analysis of mRNAs, miRNAs and lncRNAs in the Longissimus dorsi Muscles between Fat-Type and Lean-Type Pigs

In pigs, meat quality and production are two important traits affecting the pig industry and human health. Compared to lean-type pigs, fat-type pigs contain higher intramuscular fat (IMF) contents, better taste and nutritional value. To uncover genetic factors controlling differences related to IMF in pig muscle, we performed RNA-seq analysis on the transcriptomes of the Longissimus dorsi (LD) muscle of Laiwu pigs (LW, fat-type pigs) and commercial Duroc × Landrace × Yorkshire pigs (DLY, lean-type pigs) at 150 d to compare the expression profiles of mRNA, miRNA and lncRNA. A total of 225 mRNAs, 12 miRNAs and 57 lncRNAs were found to be differentially expressed at the criteria of |log2(foldchange)| > 1 and q < 0.05. The mRNA expression of LDHB was significantly higher in the LD muscle of LW compared to DLY pigs with log2(foldchange) being 9.66. Using protein interaction prediction method, we identified more interactions of estrogen-related receptor alpha (ESRRA) associated with upregulated mRNAs, whereas versican (VCAN) and proenkephalin (PENK) were associated with downregulated mRNAs in LW pigs. Integrated analysis on differentially expressed (DE) mRNAs and miRNAs in the LD muscle between LW and DLY pigs revealed two network modules: between five upregulated mRNA genes (GALNT15, FKBP5, PPARGC1A, LOC110258214 and LOC110258215) and six downregulated miRNA genes (ssc-let-7a, ssc-miR190-3p, ssc-miR356-5p, ssc-miR573-5p, ssc-miR204-5p and ssc-miR-10383), and between three downregulated DE mRNA genes (IFRD1, LOC110258600 and LOC102158401) and six upregulated DE miRNA genes (ssc-miR1379-3p, ssc-miR1379-5p, ssc-miR397-5p, ssc-miR1358-5p, ssc-miR299-5p and ssc-miR1156-5p) in LW pigs. Based on the mRNA and ncRNA binding site targeting database, we constructed a regulatory network with miRNA as the center and mRNA and lncRNA as the target genes, including GALNT15/ssc-let-7a/LOC100523888, IFRD1/ssc-miR1379-5p/CD99, etc., forming a ceRNA network in the LD muscles that are differentially expressed between LW and DLY pigs. Collectively, these data may provide resources for further investigation of molecular mechanisms underlying differences in meat traits between lean- and fat-type pigs.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in Muscle Development of Wu'an Goat Longissimus Dorsi Tissues

The growth and development of skeletal muscle is a physiological process regulated by a variety of genes and signaling pathways. As a posttranscriptional regulatory factor, circRNA plays a certain regulatory role in the development of animal skeletal muscle in the form of a miRNA sponge. However, the role of circRNAs in muscle development and growth in goats is still unclear. In our study, apparent differences in muscle fibers in Wu'an goats of different ages was firstly detected by hematoxylin-eosin (HE) staining, the circRNA expression profiles of longissimus dorsi muscles from 1-month-old (mon1) and 9-month-old (mon9) goats were screened by RNA-seq and verified by RT-qPCR. The host genes of differentially expressed (DE) circRNAs were predicted, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses (KEGG) of host genes with DE circRNAs were performed to explore the functions of circRNAs. The circRNA-miRNA-mRNA networks were then constructed using Cytoscape software. Ten significantly differentially expressed circRNAs were also verified in the mon1 and mon9 groups by RT-qPCR. Luciferase Reporter Assay was used to verify the binding site between circRNA and its targeted miRNA. The results showed that a total of 686 DE circRNAs were identified between the mon9 and mon1 groups, of which 357 were upregulated and 329 were downregulated. Subsequently, the 467 host genes of DE circRNAs were predicted using Find_circ and CIRI software. The circRNA-miRNA-mRNA network contained 201 circRNAs, 85 miRNAs, and 581 mRNAs; the host mRNAs were associated with "muscle fiber development" and "AMPK signaling pathway" and were enriched in the FoxO signaling pathway. Competing endogenous RNA (ceRNA) network analysis showed that novel_circ_0005314, novel_circ_0005319, novel_circ_0009256, novel_circ_0009845, novel_circ_0005934 and novel_circ_0000134 may play important roles in skeletal muscle growth and development between the mon9 and mon1 groups. Luciferase Reporter Assay confirmed the combination between novel_circ_0005319 and chi-miR-199a-5p, novel_circ_0005934 and chi-miR-450-3p and novel_circ_0000134 and chi-miR-655. Our results provide specific information related to goat muscle development and a reference for the goat circRNA profile.

Post-synaptic scaffold protein TANC2 in psychiatric and somatic disease risk

Understanding the shared genetic aetiology of psychiatric and medical comorbidity in neurodevelopmental disorders (NDDs) could improve patient diagnosis, stratification and treatment options. Rare tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (TANC2)-disrupting variants were disease causing in NDD patients. The post-synaptic scaffold protein TANC2 is essential for dendrite formation in synaptic plasticity and plays an unclarified but critical role in development. We here report a novel homozygous-viable Tanc2-disrupted function model in which mutant mice were hyperactive and had impaired sensorimotor gating consistent with NDD patient psychiatric endophenotypes. Yet, a multi-systemic analysis revealed the pleiotropic effects of Tanc2 outside the brain, such as growth failure and hepatocellular damage. This was associated with aberrant liver function including altered hepatocellular metabolism. Integrative analysis indicates that these disrupted Tanc2 systemic effects relate to interaction with Hippo developmental signalling pathway proteins and will increase the risk for comorbid somatic disease. This highlights how NDD gene pleiotropy can augment medical comorbidity susceptibility, underscoring the benefit of holistic NDD patient diagnosis and treatment for which large-scale preclinical functional genomics can provide complementary pleiotropic gene function information.

Summary: Disruption of mouse Tanc2 causes brain and liver abnormality, increasing psychiatric and somatic disease risk long term, highlighting the benefit of holistic diagnosis and treatment approaches for human neurodevelopmental disorder.

Trait correlated expression combined with eQTL and ASE analyses identified novel candidate genes affecting intramuscular fat

BACKGROUND

Intramuscular fat (IMF) content is a determining factor for meat taste. The Luchuan pig is a fat-type local breed in southern China that is famous for its desirable meat quality due to high IMF, however, the crossbred offspring of Luchuan sows and Duroc boars displayed within-population variation on meat quality, and the reason remains unknown.

RESULTS

In the present study, we identified 212 IMF-correlated genes (FDR ≤ 0.01) using correlation analysis between gene expression level and the value of IMF content. The IMF-correlated genes were significantly enriched in the processes of lipid metabolism and mitochondrial energy metabolism, as well as the AMPK/PPAR signaling pathway. From the IMF-correlated genes, we identified 99 genes associated with expression quantitative trait locus (eQTL) or allele-specific expression (ASE) signals, including 21 genes identified by both cis-eQTL and ASE analyses and 12 genes identified by trans-eQTL analysis. Genome-wide association study (GWAS) of IMF identified a significant QTL on SSC14 (p-value = 2.51E^-7), and the nearest IMF-correlated gene SFXN4 (r = 0.28, FDR = 4.00E^-4) was proposed as the candidate gene. Furthermore, we highlighted another three novel IMF candidate genes, namely AGT, EMG1, and PCTP, by integrated analysis of GWAS, eQTL, and IMF-gene correlation analysis.

CONCLUSIONS

The AMPK/PPAR signaling pathway together with the processes of lipid and mitochondrial energy metabolism plays a vital role in regulating porcine IMF content. Trait correlated expression combined with eQTL and ASE analysis highlighted a priority list of genes, which compensated for the shortcoming of GWAS, thereby accelerating the mining of causal genes of IMF.

Transcriptomic Profiling of Skeletal Muscle Reveals Candidate Genes Influencing Muscle Growth and Associated Lipid Composition in Portuguese Local Pig Breeds

Simple Summary

Screening and interpretation of differentially expressed genes and associated biological pathways was conducted among experimental groups with divergent phenotypes providing valuable information about the metabolic events occurring and identification of candidate genes with major regulation roles. This comparative transcriptomic analysis includes the first RNA-seq analysis of the Longissimus lumborum muscle tissue from two Portuguese autochthonous pig breeds with different genetic backgrounds, Alentejano and Bísaro. Moreover, a complementary candidate gene approach was employed to analyse, by real time qPCR, the expression profile of relevant genes involved in lipid metabolism, and therefore with potential impacts on meat composition. This study contributes to explaining the biological basis of phenotypical differences occurring between breeds, particularly the ones related to meat quality traits that affect consumer interest.

Abstract

Gene expression is one of the main factors to influence meat quality by modulating fatty acid metabolism, composition, and deposition rates in muscle tissue. This study aimed to explore the transcriptomics of the Longissimus lumborum muscle in two local pig breeds with distinct genetic background using next-generation sequencing technology and Real-Time qPCR. RNA-seq yielded 49 differentially expressed genes between breeds, 34 overexpressed in the Alentejano (AL) and 15 in the Bísaro (BI) breed. Specific slow type myosin heavy chain components were associated with AL (MYH7) and BI (MYH3) pigs, while an overexpression of MAP3K14 in AL may be associated with their lower loin proportion, induced insulin resistance, and increased inflammatory response via NFkB activation. Overexpression of RUFY1 in AL pigs may explain the higher intramuscular (IMF) content via higher GLUT4 recruitment and consequently higher glucose uptake that can be stored as fat. Several candidate genes for lipid metabolism, excluded in the RNA-seq analysis due to low counts, such as ACLY, ADIPOQ, ELOVL6, LEP and ME1 were identified by qPCR as main gene factors defining the processes that influence meat composition and quality. These results agree with the fatter profile of the AL pig breed and adiponectin resistance can be postulated as responsible for the overexpression of MAP3K14′s coding product NIK, failing to restore insulin sensitivity.

Effects of Intensive Fattening With Total Mixed Rations on Carcass Characteristics, Meat Quality, and Meat Chemical Composition of Yak and Mechanism Based on Serum and Transcriptomic Profiles

The objective of this study was to investigate the effects of intensive fattening with total mixed rations (TMR) on carcass characteristics, meat quality, and chemical composition of the yak meat. Theoretical data has been provided for evaluating the quality of yak meat during natural grazing and short-term fattening. Based on the analysis, we found that in fattening yak, the carcass weight (CWT) was increased by 106.43%, whereas the cooking loss, tenderness, and drop loss were significantly improved due to higher intramuscular fat content and lower moisture (P < 0.05). Protein, fat, calcium, and amino acids were also much higher (P < 0.01) in fattening yak compared with the grazing yak. The levels of albumin (ALB), lactate dehydrogenase (LDH), triglyceride (TRIG), and amylase (AMYL) in serum indicated better nutritional status for fattening yaks. The transcriptomics analysis showed that the high expression of ACSL1 and ACACB genes improved the synthesis and deposition of fat in fattening yak, whereas the regulation of SLC7A8, ATP1A4, ATP1A1, SLC3A2, and CPA3 gene expression weakened the proteolysis. These results indicated that fattening with TMR improves the yield and quality of the yak meat.

Label-free quantitative proteomic analysis of M. longissimus dorsi from cattle during dietary restriction and subsequent compensatory growth

Compensatory growth (CG) is a naturally occurring physiological process whereby an animal has the ability to undergo enhanced growth following a period of restricted feeding. This studies objective was to identify key proteins involved in the expression of CG. Forty Holstein Friesian bulls were equally assigned to one of four groups. R1 and R2 groups were subjected to restricted feed allowance for 125 days (Period 1). A1 and A2 animals had ad libitum access to feed in Period 1. Following Period 1, all animals from R1 and A1 were slaughtered. Remaining animals (R2 and A2) were slaughtered following ad libitum access to feed for successive 55 days (Period 2). M. longissimus dorsi samples were collected at slaughter from all animals. Proteins were isolated from samples and subjected to label-free mass spectrometry proteomic quantification. Proteins which were differentially abundant during CG (n = 39) were involved in cellular binding processes, oxidative phosphorylation and mitochondrial function. There was also evidence for up regulation of three pathways involved in nucleotide biosynthesis. Genetic variants in or regulating genes pertaining to proteins identified in this study may hold potential for use as DNA based biomarkers for genomic selection of animals with a greater ability to undergo CG.

Predicting the Quality of Meat: Myth or Reality?

This review is aimed at providing an overview of recent advances made in the field of meat quality prediction, particularly in Europe. The different methods used in research labs or by the production sectors for the development of equations and tools based on different types of biological (genomic or phenotypic) or physical (spectroscopy) markers are discussed. Through the various examples, it appears that although biological markers have been identified, quality parameters go through a complex determinism process. This makes the development of generic molecular tests even more difficult. However, in recent years, progress in the development of predictive tools has benefited from technological breakthroughs in genomics, proteomics, and metabolomics. Concerning spectroscopy, the most significant progress was achieved using near-infrared spectroscopy (NIRS) to predict the composition and nutritional value of meats. However, predicting the functional properties of meats using this method-mainly, the sensorial quality-is more difficult. Finally, the example of the MSA (Meat Standards Australia) phenotypic model, which predicts the eating quality of beef based on a combination of upstream and downstream data, is described. Its benefit for the beef industry has been extensively demonstrated in Australia, and its generic performance has already been proven in several countries.

Effects of long-term feeding of rapeseed meal on skeletal muscle transcriptome, production efficiency and meat quality traits in Norwegian Landrace growing-finishing pigs

This study was performed to investigate the effects of dietary inclusion of 20% rapeseed meal (RSM) as an alternative to soybean meal (SBM) in a three-month feeding experiment with growing finishing pigs. Dietary alteration affected growth performance, several carcass traits and transcriptional responses in the skeletal muscle, but did not affect measured meat quality traits. In general, pigs fed the RSM test diet exhibited reduced growth performance compared to pigs on SBM control diet. Significant transcriptional changes in the skeletal muscle of growing pigs fed RSM diet were likely the consequence of an increased amount of fiber and higher polyunsaturated fatty acids, and presence of bioactive phytochemicals, such as glucosinolates. RNAseq pipeline using Tophat2-Cuffdiff identified 57 upregulated and 63 downregulated genes in RSM compared to SBM pigs. Significantly enriched among downregulated pathways was p53-mediated signalling involved in cellular proliferation, while activation of negative growth regulators (IER5, KLF10, BTG2, KLF11, RETREG1, PRUNE2) in RSM fed pigs provided further evidence for reduced proliferation and increased cellular death, in accordance with the observed reduction in performance traits. Upregulation of well-known metabolic controllers (PDK4, UCP3, ESRRG and ESRRB), involved in energy homeostasis (glucose and lipid metabolism, and mitochondrial function), suggested less available energy and nutrients in RSM pigs. Furthermore, several genes supported more pronounced proteolysis (ABTB1, OTUD1, PADI2, SPP1) and reduced protein synthesis (THBS1, HSF4, AP1S2) in RSM muscle tissue. In parallel, higher levels of NR4A3, PDK4 and FGF21, and a drop in adropin, ELOVL6 and CIDEC/FSP27 indicated increased lipolysis and fatty acid oxidation, reflective of lower dressing percentage. Finally, pigs exposed to RSM showed greater expression level of genes responsive to oxidative stress, indicated by upregulation of GPX1, GPX2, and TXNIP.

Data from the Farmgate-to-Meat Continuum Including Omics-Based Biomarkers to Better Understand the Variability of Beef Tenderness: An Integromics Approach

This study is based on an integromic approach of 71 young bulls' data from the farmgate-to-meat continuum including omics-based biomarkers, to understand beef tenderness variability in two muscle cuts that differ by their contractile and metabolic properties. By the means of chemometrics using partial least-squares (PLS) and principal component regressions (PCR), important variables from a list of 49 that characterize four levels of the continuum (rearing factors-carcass-muscle-meat) were identified to explain tenderness of Longissimus thoracis (LT) and Semitendinosus (ST) muscles evaluated by a sensory panel and instrumental Warner-Bratzler shear force (WBSF). The PLS and PCR analyses validated 16 and 15 variables for LT and 12 and 14 for ST from the whole continuum to explain sensory tenderness and WBSF, respectively. Among the explanatory variables in the four models and in line with the role of apoptosis in tenderness determinism, HSP70-1A/B (a heat shock protein) was retained to explain beef tenderness irrespective of muscle and evaluation method. Similarly, dressing percentage from the carcass level was another robust predictor but in a muscle-dependent direction manner. HSP20, ENO3, and MyHC-I as three muscle protein biomarkers and dry matter intake (DMI) as a rearing factor were involved in three models to explain beef tenderness. This study highlighted also that several variables were muscle-specific irrespective of the evaluation method of tenderness. For LT muscle, six variables including three carcass traits (fatness score, fat carcass %, and muscle carcass %), two muscle biomarkers (HSP70-8 and MyHC-IIx/b), and one meat quality trait (pH_3h) were found. For ST muscle, five variables were validated from two rearing factors (average daily gain and feed efficiency) and three structural protein biomarkers (α-actin, MyBP-H, and CapZ-β). Finally, for WBSF only, lactate dehydrogenase chain B (LDH-B) was retained positively for LT and negatively for ST muscles. Overall, this trial showed that tenderness of LT and ST muscle cuts is influenced by variables belonging to the whole continuum with relationships that depend on both the muscle type and the evaluation method. It further highlighted the potential of integromic/chemometric approaches on the farmgate-to-meat continuum data to better understand the sophisticated biological processes that orchestrate the conversion of muscle into meat and tenderness determinism.

Genome-wide analysis of circular RNAs in prenatal and postnatal muscle of sheep

Circular RNAs (circRNAs), a type of non-coding RNA with circular structure, were generated by back splicing and widely expressed in animals and plants. Recent studies have shown that circRNAs extensively participate in cell proliferation, cell differentiation, cell autophagy and other biological processes. However, the role and expression of circRNAs in the development and growth of muscle have not been studied in sheep. In our study, we first used RNA-seq to study the circRNAs in prenatal and postnatal longissimus dorsi muscle of sheep. A total of 6113 circRNAs were detected from the RNA-seq data. Several circRNAs were identified using reverse transcription PCR, DNA sequencing and RNase R digestion experiments. The expression levels of several circRNAs in prenatal and postnatal muscle were confirmed by Real-Time RT-PCR. The gene ontology (GO) and KEGG enrichment analysis of the host gene of the circRNAs showed that these circRNAs were mainly involved in the growth and development of muscle related signaling pathways. These circRNAs might sponge microRNAs (miRNAs) in predicted circRNA-miRNA-mRNA networks. The circRNAs expression profiles in muscle provided an important reference for the study of circRNAs in sheep.

Analysis of extended haplotype in Korean cattle (Hanwoo) population

Korean cattle (Hanwoo) are categorized into three breeds based on color: brown, brindle, and black. Among these breeds, brown Hanwoo has been subjected to intensive selection to improve meat traits. To identify genetic traces driven by recent selection in brown Hanwoo, we scanned the genomes of brown and brindle Hanwoo using a bovine SNP chip. We identified 17 candidate selection signatures in brown Hanwoo and sequenced four candidate regions from 10 individuals each of brown and brindle Hanwoo. In particular, non-synonymous SNPs in the ADSL gene (K88M, L189H, and R302Q) might have had mutational effects on protein structure as a result of altering the purine pathway during nucleotide breakdown. The ADSL gene was previously reported to affect meat quality and yield in livestock. Meat quality and yield are main breeding goals for brown Hanwoo, and our results support a potential causal influence of non-synonymous SNPs in the ADSL gene. [BMB Reports 2016; 49(9): 514-519]

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio)

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F₂ resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

Integrative Analysis of Metabolomic, Proteomic and Genomic Data to Reveal Functional Pathways and Candidate Genes for Drip Loss in Pigs

The aim of this study was to integrate multi omics data to characterize underlying functional pathways and candidate genes for drip loss in pigs. The consideration of different omics levels allows elucidating the black box of phenotype expression. Metabolite and protein profiling was applied in Musculus longissimus dorsi samples of 97 Duroc × Pietrain pigs. In total, 126 and 35 annotated metabolites and proteins were quantified, respectively. In addition, all animals were genotyped with the porcine 60 k Illumina beadchip. An enrichment analysis resulted in 10 pathways, amongst others, sphingolipid metabolism and glycolysis/gluconeogenesis, with significant influence on drip loss. Drip loss and 22 metabolic components were analyzed as intermediate phenotypes within a genome-wide association study (GWAS). We detected significantly associated genetic markers and candidate genes for drip loss and for most of the metabolic components. On chromosome 18, a region with promising candidate genes was identified based on SNPs associated with drip loss, the protein "phosphoglycerate mutase 2" and the metabolite glycine. We hypothesize that association studies based on intermediate phenotypes are able to provide comprehensive insights in the genetic variation of genes directly involved in the metabolism of performance traits. In this way, the analyses contribute to identify reliable candidate genes.

Integrative analysis of transcriptomics and proteomics of skeletal muscles of the Chinese indigenous Shaziling pig compared with the Yorkshire breed

Background

The Shaziling pig (Sus scrofa) is a well-known indigenous breed in China. One of its main advantages over European breeds is its high meat quality. However, little genetic information is available for the Shaziling pig. To screen for differentially expressed genes and proteins that might be responsible for the meat quality, the longissimus dorsi muscles from Shaziling and Yorkshire pig breeds were investigated using an integrative analysis of transcriptomics and proteomics, involving high-throughput sequencing, the two-dimensional gel electrophoresis, and mass spectrometry.

Results

Sequencing produced 79,320 unigenes by de novo assembly, and 488 differentially expressed genes in the longissimus dorsi muscle of Shaziling pig compared with the Yorkshire breed were identified. Gene Ontology term enrichment of biological functions and Kyoto Encyclopedia of Genes and Genomes analysis showed that the gene products were mainly involved in metabolism, protein binding, and regulation of skeletal muscle development. At the protein level, 23 differentially expressed proteins were identified, which were potentially associated with fatty acid metabolism, the glycolytic pathway, and skeletal muscle growth. Eight differentially expressed genes were confirmed by real-time PCR. These results give an insight into the mechanisms underlying the formation of skeletal muscle in the Shaziling pig.

Conclusions

Certain differentially expressed genes and proteins are involved in fatty acid metabolism, intramuscular fat deposition, and skeletal muscle growth in the Shaziling pig. These results provide candidate genes for improving meat quality and will promote further transcriptomic research in Shaziling pigs.

Electronic supplementary material

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

Beef quality with different intramuscular fat content and proteomic analysis using isobaric tag for relative and absolute quantitation of differentially expressed proteins

Intramuscular fat (IMF) is an important trait for beef eating quality. The mechanism of how IMF is deposited in beef cattle muscle is not clear at the molecular level. The muscle (M. longissimus lumborum: LL) of a group of Xiangxi yellow×Angus cattle with high fat levels (HF), was compared to the muscle of a low fat group (LF). The meat quality and the expressed protein patterns were compared. It was shown that LL from the HF animals had a greater fat content (P<0.05) and lower moisture content (P<0.05) than LL from LF animals. Forty seven sarcoplasmic proteins were differentially expressed and identified between the two groups. These proteins are involved in 6 molecular functions and 16 biological processes, and affect the Mitogen-activated protein kinases pathway, insulin pathway and c-Jun N-terminal kinases leading to greater IMF deposition. Cattle in the HF group had greater oxidative capacity and lower glycolytic levels suggesting a greater energetic efficiency.

Different Statistical Approaches to Investigate Porcine Muscle Metabolome Profiles to Highlight New Biomarkers for Pork Quality Assessment

The aim of this study was to elucidate the underlying biochemical processes to identify potential key molecules of meat quality traits drip loss, pH of meat 1 h post-mortem (pH1), pH in meat 24 h post-mortem (pH24) and meat color. An untargeted metabolomics approach detected the profiles of 393 annotated and 1,600 unknown metabolites in 97 Duroc × Pietrain pigs. Despite obvious differences regarding the statistical approaches, the four applied methods, namely correlation analysis, principal component analysis, weighted network analysis (WNA) and random forest regression (RFR), revealed mainly concordant results. Our findings lead to the conclusion that meat quality traits pH1, pH24 and color are strongly influenced by processes of post-mortem energy metabolism like glycolysis and pentose phosphate pathway, whereas drip loss is significantly associated with metabolites of lipid metabolism. In case of drip loss, RFR was the most suitable method to identify reliable biomarkers and to predict the phenotype based on metabolites. On the other hand, WNA provides the best parameters to investigate the metabolite interactions and to clarify the complex molecular background of meat quality traits. In summary, it was possible to attain findings on the interaction of meat quality traits and their underlying biochemical processes. The detected key metabolites might be better indicators of meat quality especially of drip loss than the measured phenotype itself and potentially might be used as bio indicators.

Recent advances in omic technologies for meat quality management

The knowledge of the molecular organization of living organisms evolved considerably during the last years. The methodologies associated also progressed with the development of the high-throughput sequencing (SNP array, RNAseq, etc.) and of genomic tools allowing the simultaneous analysis of hundreds or thousands of genes, proteins or metabolites. In farm animals, some proteins, mRNAs or metabolites whose abundance has been associated with meat quality traits have been detected in pig, cattle, chicken. They constitute biomarkers for the assessment and prediction of qualities of interest in each species, with potential biomarkers across species. The ongoing development of rapid methods will allow their use for decision-making and management tools in slaughterhouses, to better allocate carcasses or cuts to the appropriate markets. Besides, their application on living animals will help to improve genetic selection and to adapt a breeding system to fulfill expected quality level. The ultimate goal is to propose effective molecular tools for the management of product quality in meat production chains.

An Optimized Trichloroacetic Acid/Acetone Precipitation Method for Two-Dimensional Gel Electrophoresis Analysis of Qinchuan Cattle Longissimus Dorsi Muscle Containing High Proportion of Marbling

Longissimus dorsi muscle (LD) proteomics provides a novel opportunity to reveal the molecular mechanism behind intramuscular fat deposition. Unfortunately, the vast amounts of lipids and nucleic acids in this tissue hampered LD proteomics analysis. Trichloroacetic acid (TCA)/acetone precipitation is a widely used method to remove contaminants from protein samples. However, the high speed centrifugation employed in this method produces hard precipitates, which restrict contaminant elimination and protein re-dissolution. To address the problem, the centrifugation precipitates were first grinded with a glass tissue grinder and then washed with 90% acetone (TCA/acetone-G-W) in the present study. According to our result, the treatment for solid precipitate facilitated non-protein contaminant removal and protein re-dissolution, ultimately improving two-dimensional gel electrophoresis (2-DE) analysis. Additionally, we also evaluated the effect of sample drying on 2-DE profile as well as protein yield. It was found that 30 min air-drying did not result in significant protein loss, but reduced horizontal streaking and smearing on 2-DE gel compared to 10 min. In summary, we developed an optimized TCA/acetone precipitation method for protein extraction of LD, in which the modifications improved the effectiveness of TCA/acetone method.

The Longissimus and Semimembranosus muscles display marked differences in their gene expression profiles in pig

BACKGROUND

Meat quality depends on skeletal muscle structure and metabolic properties. While most studies carried on pigs focus on the Longissimus muscle (LM) for fresh meat consumption, Semimembranosus (SM) is also of interest because of its importance for cooked ham production. Even if both muscles are classified as glycolytic muscles, they exhibit dissimilar myofiber composition and metabolic characteristics. The comparison of LM and SM transcriptome profiles undertaken in this study may thus clarify the biological events underlying their phenotypic differences which might influence several meat quality traits.

METHODOLOGY/PRINCIPAL FINDINGS

Muscular transcriptome analyses were performed using a custom pig muscle microarray: the 15 K Genmascqchip. A total of 3823 genes were differentially expressed between the two muscles (Benjamini-Hochberg adjusted P value ≤0.05), out of which 1690 and 2133 were overrepresented in LM and SM respectively. The microarray data were validated using the expression level of seven differentially expressed genes quantified by real-time RT-PCR. A set of 1047 differentially expressed genes with a muscle fold change ratio above 1.5 was used for functional characterization. Functional annotation emphasized five main clusters associated to transcriptome muscle differences. These five clusters were related to energy metabolism, cell cycle, gene expression, anatomical structure development and signal transduction/immune response.

CONCLUSIONS/SIGNIFICANCE

This study revealed strong transcriptome differences between LM and SM. These results suggest that skeletal muscle discrepancies might arise essentially from different post-natal myogenic activities.

Transcriptome analysis of porcine M. semimembranosus divergent in intramuscular fat as a consequence of dietary protein restriction

Background

Intramuscular fat (IMF) content is positively correlated with aspects of pork palatability, including flavour, juiciness and overall acceptability. The ratio of energy to protein in the finishing diet of growing pigs can impact on IMF content with consequences for pork quality. The objective of this study was to compare gene expression profiles of Musculus semimembranosus (SM) of animals divergent for IMF as a consequence of protein dietary restriction in an isocaloric diet. The animal model was derived through the imposition of low or high protein diets during the finisher stage in Duroc gilts. RNA was extracted from post mortem SM tissue, processed and hybridised to Affymetrix porcine GeneChip® arrays.

Results

IMF content of SM muscle was increased on the low protein diet (3.60 ± 0.38% versus 1.92 ± 0.35%). Backfat depth was also greater in animals on the low protein diet, and average daily gain and feed conversion ratio were lower, but muscle depth, protein content and moisture content were not affected. A total of 542 annotated genes were differentially expressed (DE) between animals on low and high protein diets, with 351 down-regulated and 191 up-regulated on the low protein diet. Transcript differences were validated for a subset of DE genes by qPCR. Alterations in functions related to cell cycle, muscle growth, extracellular matrix organisation, collagen development, lipogenesis and lipolysis, were observed. Expression of adipokines including LEP, TNFα and HIF1α were increased and the hypoxic stress response was induced. Many of the identified transcriptomic responses have also been observed in genetic and fetal programming models of differential IMF accumulation, indicating they may be robust biological indicators of IMF content.

Conclusion

An extensive perturbation of overall energy metabolism in muscle occurs in response to protein restriction. A low protein diet can modulate IMF content of the SM by altering gene pathways involved in lipid biosynthesis and degradation; however this nutritional challenge negatively impacts protein synthesis pathways, with potential consequences for growth.

Transcriptional analysis of intramuscular fatty acid composition in the longissimus thoracis muscle of Iberian × Landrace back-crossed pigs

This study aimed at identifying differential gene expression conditional on the fatty acid profile of the longissimus thoracis (Lt) muscle, a prime cut of economic relevance for fresh and cured pork production. A population of 110 Iberian (25%) × Landrace (75%) back-crossed pigs was used, because these two breeds exhibit extreme profiles of intramuscular saturated fatty acid, monounsaturated fatty acid (MUFA) and polyunsaturated fatty acid (PUFA) contents. Total RNA from Lt muscle was individually hybridized to GeneChip Porcine Genome arrays (Affymetrix). A principal component analysis was performed with data from the 110 animals to select 40 extreme animals based on the total fatty acid profile and the MUFA composition (MAP). Comparison of global transcription levels between extreme fatty acid profile pigs (n = 40) resulted in 219 differentially expressed probes (false discovery rate <0.10). Gene ontology, pathway and network analysis indicated that animals with higher percentages of PUFA exhibit a shift toward a more oxidative muscular metabolism state, with a raise in mitochondria function (PPARGC1A, ATF2), fatty acid uptake and oxidation (FABP5, MGLL). On the other hand, 87 probes were differentially expressed between MUFA composition groups (n = 40; false discovery rate <0.10). In particular, muscles rich in n-7 MUFA expressed higher levels of genes involved in lipid metabolism (GLUL, CRAT, PLA2G15) and lower levels of fatty acid elongation genes (ELOVL5). Moreover, the chromosomal position of FABP5, PAQR3, MGLL, PPARGC1A, GLUL and ELOVL5 co-localized with very relevant QTL for fat deposition and composition described in the same resource population. This study represents a complementary approach to identifying genes underlying these QTL effects.

Associations between muscle gene expression pattern and technological and sensory meat traits highlight new biomarkers for pork quality assessment

Meat quality (MQ) results from complex phenomenon and despite improved knowledge on MQ development, its variability remains high. The identification of biomarkers and the further development of rapid tests would thus be helpful to evaluate MQ in pork industries. Using transcriptomics, the present study aimed at identifying biomarkers of eight pork quality traits: ultimate pH, drip loss, lightness, redness, hue angle, intramuscular fat, shear force and tenderness, based on an experimental design inducing a high variability in MQ. Associations between microarray gene expression and pork traits (n=50 pigs) highlighted numerous potential biomarkers of MQ. Using quantitative RT-PCR, 113 transcript-trait correlations including 40 of these genes were confirmed (P<0.05, |r|≤0.73), out of which 60 were validated (P<0.05, |r|≤0.68) on complementary experimental data (n=50). Multiple regression models including 3 to 5 genes explained up to 59% of MQ trait variability. Moreover, functional analysis of correlated-trait genes provided information on the biological phenomena underlying MQ.

Differential effects of reduced protein diets on fatty acid composition and gene expression in muscle and subcutaneous adipose tissue of Alentejana purebred and Large White × Landrace × Pietrain crossbred pigs

The present study assessed the effect of pig genotype (fatty v. lean) and dietary protein and lysine (Lys) levels (normal v. reduced) on intramuscular fat (IMF) content, subcutaneous adipose tissue (SAT) deposition, fatty acid composition and mRNA levels of genes controlling lipid metabolism. The experiment was conducted on sixty intact male pigs (thirty Alentejana purebred and thirty Large White × Landrace × Pietrain crossbred), from 60 to 93 kg of live weight. Animals were divided into three groups fed with the following diets: control diet equilibrated for Lys (17·5 % crude protein (CP) and 0·7 % Lys), reduced protein diet (RPD) equilibrated for Lys (13·2 % CP and 0·6 % Lys) and RPD not equilibrated for Lys (13·1 % CP and 0·4 % Lys). It was shown that the RPD increased fat deposition in the longissimus lumborum muscle in the lean but not in the fatty pig genotype. It is strongly suggested that the effect of RPD on the longissimus lumborum muscle of crossbred pigs is mediated via Lys restriction. The increase in IMF content under the RPD was accompanied by increased stearoyl-CoA desaturase (SCD) and PPARG mRNA levels. RPD did not alter backfat thickness, but increased the total fatty acid content in both lean and fatty pig genotype. The higher amount of SAT in fatty pigs, when compared with the lean ones, was associated with the higher expression levels of ACACA, CEBPA, FASN and SCD genes. Taken together, the data indicate that the mechanisms regulating fat deposition in pigs are genotype and tissue specific, and are associated with the expression regulation of the key lipogenic genes.

Expression levels of candidate genes for intramuscular fat deposition in two Banna mini-pig inbred lines divergently selected for fatness traits

Intramuscular fat (IMF) content plays an important role in meat quality. Many genes involved in lipid and energy metabolism were identified as candidate genes for IMF deposition, since genetic polymorphisms within these genes were associated with IMF content. However, there is less information on the expression levels of these genes in the muscle tissue. This study aimed at investigating the expression levels of sterol regulating element binding protein-1c (SREBP-1c), diacylglycerol acyltransferase (DGAT-1), heart-fatty acids binding protein (H-FABP), leptin receptor (LEPR) and melanocortin 4 receptor (MC4R) genes and proteins in two divergent Banna mini-pig inbred lines (BMIL). A similar growth performance was found in both the fat and the lean BMIL. The fat meat and IMF content in the fat BMIL were significantly higher than in the lean BMIL, but the lean meat content was lower. The serum triacylglycerol (TAG) and free fatty acid (FFA) contents were significantly higher in the fat than in the lean BMIL. The expression levels of SREBP-1c, DGAT-1 and H-FABP genes and proteins in fat BMIL were increased compared to the lean BMIL. However, the expression levels of LEPR and MC4R genes and proteins were lower.

Functional analysis of muscle gene expression profiles associated with tenderness and intramuscular fat content in pork

Warner-Bratzler shear force (WBSF) and % intramuscular fat content (IMF) are objective meat quality measurements that are significantly correlated with aspects of palatability such as tenderness, flavour and juiciness. Using cDNA microarrays, Musculus longissimus transcriptomic profiles at slaughter were compared in samples displaying lower or higher IMF (n=8) and WBSF values on day 1 post mortem (n=8). 101 identified genes were differentially expressed in relation to WBSF, while 160 genes were associated with differences in IMF. Reduced expression of protein synthesis genes and enhanced expression of genes involved in protein degradation were associated with lower WBSF values on day 1. Pathways including oxidative phosphorylation and the citrate cycle were significantly associated with higher IMF. Many lipid oxidation and fatty acid metabolism pathway genes were down-regulated in high IMF tissue, suggesting a suppression of fatty acid turnover in muscle with higher fat content. Identified genes provide targets for the discovery of novel genetic variation influential on pork palatability.

Saliva proteomics as an emerging, non-invasive tool to study livestock physiology, nutrition and diseases

Saliva is an extraordinary fluid in terms of research and diagnostic possibilities. Its composition in electrolytes, hormones and especially its proteome contains information about feeding status, nutritional requirements and adaptations to diet and environment, and also about health status of animals. It is easy to collect on a non-invasive and routine basis without any need for special training. Therefore, the analysis of salivary proteomes is going to emerge into a field of high interest with the future goal to maintain and improve livestock productivity and welfare. Moreover, the comprehensive analysis and identification of salivary proteins and peptides in whole and glandular saliva is a necessary pre-requisite to identify animal disease biomarkers and a powerful tool to better understand animal physiology. This review focuses on the different approaches used to study the salivary proteomes of farm animals, in respect to the physiology of nutrition and food perception in relation to food choices. The potential of animal saliva as a source of disease biomarkers will also be pointed out. Special emphasis is laid on the 'ruminating triad' - cattle, goat and sheep - as well as swine as major species of animal production in Western and Southern Europe.

Functional genomics of the muscle response to restraint and transport in chickens

In the present study, we used global approaches (proteomics, transcriptomics, and metabolomics) to assess the molecular basis of the muscle response to stress in chickens. A restraint test, combined with transport for 2 h (RT test) was chosen as the potentially stressful situation. Chickens (6 wk old) were either nontreated (control chickens) or submitted to the RT test (treated chickens). The RT test induced a 6-fold increase in corticosterone concentrations, suggesting hypothalamic-pituitary-adrenal axis activation. The RT test decreased the relative abundance of several hexose phosphates [glucose-1-P (G1P), glucose-6-P (G6P), fructose-6-P (F6P), and mannose-6-P (M6P)] in thigh muscle. In addition, 55 transcripts, among which 39 corresponded to unique annotated genes, were significantly up- (12 genes) or downregulated (27 genes) by treatment. Similarly, 45 proteic spots, among which 29 corresponded to unique annotated proteins, were overexpressed (11 proteins), underexpressed (14 proteins), or only expressed in treated chickens. Integrative analysis of differentially expressed genes and proteins showed that most transcripts and proteins belong to 2 networks whose genes were mainly related with cytoskeleton structure or carbohydrate metabolism. Whereas the decrease in energetic metabolites suggested an activation of glycogenolysis and glycolysis in response to the RT test, the reduced expression of genes and proteins involved in these pathways suggested the opposite. We hypothesized that the prolonged RT test resulted in a repression of glycogenolysis and glycolysis in thigh muscle of chickens. The down-expression of genes and proteins involved in the formation of fiber stress after the RT test suggests a reinforcement of myofibrils in response to stress.

Use of a bovine genome array to identify new biological pathways for beef marbling in Hanwoo (Korean Cattle)

BACKGROUND

Marbling (intramuscular fat) is a valuable trait that impacts on meat quality and an important factor determining price of beef in the Korean beef market. Animals that are destined for this high marbling market are fed a high concentrate ration for approximately 30 months in the Korean finishing farms. However, this feeding strategy leads to inefficiencies and excessive fat production. This study aimed to identify candidate genes and pathways associated with intramuscular fat deposition on highly divergent marbling phenotypes in adult Hanwoo cattle.

RESULTS

Bovine genome array analysis was conducted to detect differentially expressed genes (DEGs) in m. longissimus with divergent marbling phenotype (marbling score 2 to 7). Three data-processing methods (MAS5.0, GCRMA and RMA) were used to test for differential expression (DE). Statistical analysis identified 21 significant transcripts from at least two data-processing methods (P < 0.01). All 21 differentially expressed genes were validated by real-time PCR. Results showed a high concordance in the gene expression fold change between the microarrays and the real time PCR data. Gene Ontology (GO) and pathway analysis demonstrated that some genes (ADAMTS4, CYP51A and SQLE) over expressed in high marbled animals are involved in a protein catabolic process and a cholesterol biosynthesis process. In addition, pathway analysis also revealed that ADAMTS4 is activated by three regulators (IL-17A, TNFα and TGFβ1). QRT-PCR was used to investigate gene expression of these regulators in muscle with divergent intramuscular fat contents. The results demonstrate that ADAMTS4 and TGFβ1 are associated with increasing marbling fat. An ADAMTS4/TGFβ1 pathway seems to be associated with the phenotypic differences between high and low marbled groups.

CONCLUSIONS

Marbling differences are possibly a function of complex signaling pathway interactions between muscle and fat. These results suggest that ADAMTS4, which is involved in connective tissue degradation, could play a role in an important biological pathway for building up marbling in cattle. Moreover, ADAMTS4 and TGFβ1could potentially be used as an early biological marker for marbling fat content in the early stages of growth.