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

A review of porcine skeletal muscle plasticity and implications for genetic improvement of carcass and meat quality

Skeletal muscle is characterized by a remarkable plasticity to adapt to stimuli such as contractile activity, loading conditions, substrate supply or environmental factors. The existing knowledge of muscle plasticity along with developed genetic and genomic technologies, have enabled creating animal breeding strategies and allowed for implementing agriculturally successful porcine genetic improvement programs. The primary focus of this review paper is on pig skeletal muscle plasticity as it relates to genetic improvement of desirable carcass composition and pork quality traits. Biological constraints between practically realized breeding objectives, pig skeletal muscle biology, and pork quality are also discussed. Future applications of genetic and genomic technologies and plausible focus on new breeding objectives enhancing pork production sustainability are proposed as well.

Revelation of genetic diversity and genomic footprints of adaptation in Indian pig breeds

In the present study, the genetic diversity measures among four Indian domestic breeds of pig namely Agonda Goan, Ghurrah, Ghungroo, and Nicobari, of different agro-climatic regions of country were explored and compared with European commercial breeds, European wild boar and Chinese domestic breeds. The double digest restriction site-associated DNA sequencing (ddRADseq) data of Indian pigs (102) and Landrace (10 animals) were generated and whole genome sequencing data of exotic pigs (60 animals) from public data repository were used in the study. The principal component analysis (PCA), admixture analysis and phylogenetic analysis revealed that Indian breeds were closer in ancestry to Chinese breeds than European breeds. European breeds exhibited highest genetic diversity measures among all the considered breeds. Among Indian breeds, Agonda Goan and Ghurrah were found to be more genetically diverse than Nicobari and Ghungroo. The selection signature regions in Indian pigs were explored using iHS and XP-EHH, and during iHS analysis, it was observed that genes related to growth, reproduction, health, meat quality, sensory perception and behavior were found to be under selection pressure in Indian pig breeds. Strong selection signatures were recorded in 24.25-25.25 Mb region of SSC18, 123.25-124 Mb region of SSC15 and 118.75-119.5 Mb region of SSC2 in most of the Indian breeds upon pairwise comparison with European commercial breeds using XP-EHH. These regions were harboring some important genes such as EPHA4 for thermotolerance, TAS2R16, FEZF1, CADPS2 and PTPRZ1 for adaptability to scavenging system of rearing, TRIM36 and PGGT1B for disease resistance and CCDC112, PIAS1, FEM1B and ITGA11 for reproduction.

Genetic lines influenced the texture, collagen and intramuscular fat of pork longissimus and semimembranosus

The effects of genetic lines (PM-LR - Pure maternal, Landrace-type; PM-LW - Pure maternal, Large White-type; PM-D - Pure maternal, Duroc-type; SynT-LWLR - Synthetic terminal, Large white and Landrace-type; PT-D - Pure terminal, Duroc-type; and PT-LW - Pure terminal, Large White-type) on carcass traits and pork quality, including cooking loss, Warner-Bratzler shear force (WBSF), texture profile analysis (TPA), collagen content and solubility and intramuscular fat (IMF) content of Longissimus thoracis et lumborum (LTL) and Semimembranosus (SM), were investigated. The results showed that all terminal lines had lower backfat thickness than maternal lines, but muscles from line SynT-LWLR had the highest hardness and cohesiveness and the lowest IMF content. The SM showed higher hardness, adhesiveness, springiness, chewiness, collagen, and IMF content but lower WBSF and collagen solubility than LTL. Collagen and IMF content affected WBSF and TPA parameters across both muscles and in LTL alone. Significant relationships between collagen content and texture parameters were found in line PT-D. Genetic lines and muscles affect the contribution of collagen and IMF to pork texture.

Comparative muscle transcriptome of Mali and Hampshire breeds of pigs: a preliminary study

Muscle development is an important priority of pig breeding programs. There is a considerable variation in muscularity between the breeds, but the regulation mechanisms of genes underlying myogenesis are still unclear. Transcriptome data from two breeds of pigs with divergent muscularity (Mali and Hampshire) were integrated with histology, immunofluorescence and meat yield to identify differences in myogenesis during the early growth phase. The muscle transcriptomics analysis revealed 17,721 common, 1413 and 1115 unique transcripts to Hampshire and Mali, respectively. This study identified 908 differentially expressed genes (p < 0.05; log2FC > ±1) in the muscle samples, of which 550 were upregulated and 358 were downregulated in Hampshire pigs, indicating differences in physiological process related to muscle function and development. Expression of genes related to myoblast fusion (MYMK), skeletal muscle satellite cell proliferation (ANGPT1, CDON) and growth factors (HGF, IGF1, IGF2) were higher in Hampshire than Mali, even though transcript levels of several other myogenesis-related genes (MYF6, MYOG, MSTN) were similar. The number of fibers per fascicle and the expression of myogenic marker proteins (MYOD1, MYOG and PAX7) were more in Hampshire as compared to Mali breed of pig, supporting results of transcriptome studies. The results suggest that differences in muscularity between breeds could be related to the regulation of myoblast fusion and myogenic activities. The present study will help to identify genes that could be explored for their utility in the selection of animals with different muscularities.

A genome-wide association study for loin depth and muscle pH in pigs from intensely selected purebred lines

BACKGROUND

Genome-wide association studies (GWAS) aim at identifying genomic regions involved in phenotype expression, but identifying causative variants is difficult. Pig Combined Annotation Dependent Depletion (pCADD) scores provide a measure of the predicted consequences of genetic variants. Incorporating pCADD into the GWAS pipeline may help their identification. Our objective was to identify genomic regions associated with loin depth and muscle pH, and identify regions of interest for fine-mapping and further experimental work. Genotypes for ~ 40,000 single nucleotide morphisms (SNPs) were used to perform GWAS for these two traits, using de-regressed breeding values (dEBV) for 329,964 pigs from four commercial lines. Imputed sequence data was used to identify SNPs in strong ([Formula: see text] 0.80) linkage disequilibrium with lead GWAS SNPs with the highest pCADD scores.

RESULTS

Fifteen distinct regions were associated with loin depth and one with loin pH at genome-wide significance. Regions on chromosomes 1, 2, 5, 7, and 16, explained between 0.06 and 3.55% of the additive genetic variance and were strongly associated with loin depth. Only a small part of the additive genetic variance in muscle pH was attributed to SNPs. The results of our pCADD analysis suggests that high-scoring pCADD variants are enriched for missense mutations. Two close but distinct regions on SSC1 were associated with loin depth, and pCADD identified the previously identified missense variant within the MC4R gene for one of the lines. For loin pH, pCADD identified a synonymous variant in the RNF25 gene (SSC15) as the most likely candidate for the muscle pH association. The missense mutation in the PRKAG3 gene known to affect glycogen content was not prioritised by pCADD for loin pH.

CONCLUSIONS

For loin depth, we identified several strong candidate regions for further statistical fine-mapping that are supported in the literature, and two novel regions. For loin muscle pH, we identified one previously identified associated region. We found mixed evidence for the utility of pCADD as an extension of heuristic fine-mapping. The next step is to perform more sophisticated fine-mapping and expression quantitative trait loci (eQTL) analysis, and then interrogate candidate variants in vitro by perturbation-CRISPR assays.

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.

Selection and validation of differentially expressed metabolic and immune genes in weaned Ghurrah versus crossbred piglets

In the present investigation, differentially expressed genes (DEGs) were studied using RNA sequencing (RNA-seq) technique in porcine peripheral blood mononuclear cells (PBMC) of weaned Ghurrah and crossbred piglets at 3-month age. Transcriptomic analysis was done using three different packages, namely, EBSeq, DESeq2, and edgeR, to identify the DEGs between Ghurrah and crossbred piglets. Total 7717 DEGs were commonly identified by all three packages, out of which 4151 genes found to be up-regulated, and 3566 genes were down-regulated. Functional annotation of these DEGs indicated metabolism as the most commonly enriched category followed by the immune response. Genes related to metabolism and growth were up-regulated in crossbred piglets as compared with Ghurrah piglets, whereas immunity-related genes were up-regulated in Ghurrah piglets elucidating the disease resistance nature of this indigenous breed over crossbred counterparts. Further, eight DEGs, namely, LRP-1, ADCY4, ERRFI1, LDHD, ARG1, OASL, MGARP, and S100A8, were validated by qRT-PCR in a separate set of biological samples and found to be in concordance with RNA-seq results. Finding in the present study provides insight into genes and their molecular mechanisms governing difference in growth performance between Ghurrah and crossbred pigs.

Assessment of carcass and meat quality of longissimus thoracis and semimembranosus muscles of Hanwoo with Korean beef grading standards

The carcass traits, meat quality traits, and fatty acid profile of two muscles, longissimus thoracis (LT) and semimembranosus (SM), in Hanwoo steers were investigated against the Korean beef quality grade (QG) and yield grade (YG) standards. Back fat thickness and carcass weight most affect the YGs, while the marbling score is the primary determinant of the QGs. The muscle type greatly influenced the meat qualities, sensory properties, and fatty acid profiles of Hanwoo meat. In terms of sensory aspects, the muscle type affected the QG for tenderness and juiciness, with the LT being more desirable than the SM for all the sensory attributes. SM meat is potentially beneficial for human health because it has a lower amount of unhealthy saturated fatty acids and a greater amount of healthier polyunsaturated fatty acids, especially essential fatty acids, than the LT. The present study indicates that carcass based common grading system does not reflect the real value of SM meat.

Differences in gene expression profiles for subcutaneous adipose, liver, and skeletal muscle tissues between Meishan and Landrace pigs with different backfat thicknesses

Backfat thickness is one of the most important traits of commercially raised pigs. Meishan pigs are renowned for having thicker backfat than Landrace pigs. To examine the genetic factors responsible for the differences, we first produced female crossbred pig lines by mating Landrace (L) × Large White (W) × Duroc (D) females (LWD) with Landrace (L) or Meishan (M) boars (i.e., LWD × L = LWDL for Landrace offspring and LWD × M = LWDM for the Meishan offspring). We confirmed that LWDM pigs indeed had a thicker backfat than LWDL pigs. Next, we performed gene expression microarray analysis in both genetic lines to examine differentially expressed genes (DEGs) in energy metabolism-related tissues, subcutaneous adipose (fat), liver, and longissimus dorsi muscle tissues. We analyzed the annotation of DEGs (2-fold cutoff) to functionally categorize them by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. The number of DEGs in muscle tissues of both lines was much less than that in fat and liver tissues, indicating that DEGs in muscle tissues may not contribute much to differences in backfat thickness. In contrast, several genes related to muscle (in fat tissue) and lipid metabolism (in liver tissue) were more upregulated in LWDM pigs than LWDL pigs, indicating that those DEGs might be responsible for differences in backfat thickness. The different genome-wide gene expression profiles in the fat, liver, and muscle tissues between genetic lines can provide useful information for pig breeders.

Invited review: Measurable biomarkers linked to meat quality from different pig production systems

Abstract. Biological processes underlie all livestock traits, including post-mortem meat quality traits. Biomarkers are molecular components of the biological processes showing differential expression associated with the phenotype of the trait. The phenotypes of the meat quality traits are determined by the animal's genotype interacting with the environment affecting the expression of the genome. The omics technologies enable measuring the expression of the genome at all levels: transcriptome, proteome, and metabolome. Associations between the phenotype of the traits and expressions measured with the omics techniques are a first step in developing biomarkers. Biomarkers enable the monitoring, diagnosis, and prediction of changes in meat quality related to external (environmental, e.g. feed and animal management conditions) stimuli and interactions with the genotype. In this paper we review the development of biomarkers for meat quality of pigs in diverse pig breeds, environments, and pork production chains.

Molecular motors and nuclear movements in muscle

Muscle fibers have the particularity of containing numerous nuclei evenly distributed and positioned next to the plasma membrane. This unique disposition is the result of sequential events of nuclear movements that start when myoblasts fuse together and end with the clustering of few nuclei under the neuromuscular junction. Nuclei are mispositioned in multiple muscle disorders therefore the mechanisms of nuclear positioning can be novel targets for muscle disorders therapies. The 2 first nuclear movements that occur upon myoblast fusion require different microtubule motors. We performed a siRNA screen against all the microtubules motors and quantified nuclei behavior after fusion and inside the myotube. The different motors we found to be involved in the nuclear behaviors and the analysis of motors expression suggest a competition between both movement mechanisms, which potentially relies on the discrepancy between myoblast and myotube microtubules stability.

Developmental Stage, Muscle and Genetic Type Modify Muscle Transcriptome in Pigs: Effects on Gene Expression and Regulatory Factors Involved in Growth and Metabolism

Iberian pig production includes purebred (IB) and Duroc-crossbred (IBxDU) pigs, which show important differences in growth, fattening and tissue composition. This experiment was conducted to investigate the effects of genetic type and muscle (Longissimus dorsi (LD) vs Biceps femoris (BF)) on gene expression and transcriptional regulation at two developmental stages. Nine IB and 10 IBxDU piglets were slaughtered at birth, and seven IB and 10 IBxDU at four months of age (growing period). Carcass traits and LD intramuscular fat (IMF) content were measured. Muscle transcriptome was analyzed on LD samples with RNA-Seq technology. Carcasses were smaller in IB than in IBxDU neonates (p < 0.001), while growing IB pigs showed greater IMF content (p < 0.05). Gene expression was affected (p < 0.01 and Fold change > 1.5) by the developmental stage (5,812 genes), muscle type (135 genes), and genetic type (261 genes at birth and 113 at growth). Newborns transcriptome reflected a highly proliferative developmental stage, while older pigs showed upregulation of catabolic and muscle functioning processes. Regarding the genetic type effect, IBxDU newborns showed enrichment of gene pathways involved in muscle growth, in agreement with the higher prenatal growth observed in these pigs. However, IB growing pigs showed enrichment of pathways involved in protein deposition and cellular growth, supporting the compensatory gain experienced by IB pigs during this period. Moreover, newborn and growing IB pigs showed more active glucose and lipid metabolism than IBxDU pigs. Moreover, LD muscle seems to have more active muscular and cell growth, while BF points towards lipid metabolism and fat deposition. Several regulators controlling transcriptome changes in both genotypes were identified across muscles and ages (SIM1, PVALB, MEFs, TCF7L2 or FOXO1), being strong candidate genes to drive expression and thus, phenotypic differences between IB and IBxDU pigs. Many of the identified regulators were known to be involved in muscle and adipose tissues development, but others not previously associated with pig muscle growth were also identified, as PVALB, KLF1 or IRF2. The present study discloses potential molecular mechanisms underlying phenotypic differences observed between IB and IBxDU pigs and highlights candidate genes implicated in these molecular mechanisms.

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.

Comparative Analysis of Muscle Transcriptome between Pig Genotypes Identifies Genes and Regulatory Mechanisms Associated to Growth, Fatness and Metabolism

Iberian ham production includes both purebred (IB) and Duroc-crossbred (IBxDU) Iberian pigs, which show important differences in meat quality and production traits, such as muscle growth and fatness. This experiment was conducted to investigate gene expression differences, transcriptional regulation and genetic polymorphisms that could be associated with the observed phenotypic differences between IB and IBxDU pigs. Nine IB and 10 IBxDU pigs were slaughtered at birth. Morphometric measures and blood samples were obtained and samples from Biceps femoris muscle were employed for compositional and transcriptome analysis by RNA-Seq technology. Phenotypic differences were evident at this early age, including greater body size and weight in IBxDU and greater Biceps femoris intramuscular fat and plasma cholesterol content in IB newborns. We detected 149 differentially expressed genes between IB and IBxDU neonates (p < 0.01 and Fold-Change > 1. 5). Several were related to adipose and muscle tissues development (DLK1, FGF21 or UBC). The functional interpretation of the transcriptomic differences revealed enrichment of functions and pathways related to lipid metabolism in IB and to cellular and muscle growth in IBxDU pigs. Protein catabolism, cholesterol biosynthesis and immune system were functions enriched in both genotypes. We identified transcription factors potentially affecting the observed gene expression differences. Some of them have known functions on adipogenesis (CEBPA, EGRs), lipid metabolism (PPARGC1B) and myogenesis (FOXOs, MEF2D, MYOD1), which suggest a key role in the meat quality differences existing between IB and IBxDU hams. We also identified several polymorphisms showing differential segregation between IB and IBxDU pigs. Among them, non-synonymous variants were detected in several transcription factors as PPARGC1B and TRIM63 genes, which could be associated to altered gene function. Taken together, these results provide information about candidate genes, metabolic pathways and genetic polymorphisms potentially involved in phenotypic differences between IB and IBxDU pigs associated to meat quality and production traits.

Genome-wide study on intramuscular fat in Italian Large White pig breed using the PorcineSNP60 BeadChip

Genome-wide association study results are presented for intramuscular fat in Italian Large White pig breed. A total of 886 individuals were genotyped with PorcineSNP60 BeadChip. After quality control performed with plink software and in R environment, 49 208 markers remained for the association analysis. The genome-wide association studies was conducted using linear mixed model implemented in GenABEL. We detected seven new SNPs of genes till now not found associated to intramuscular fat (IMF). Three markers map in a wide intergenic region rich of QTL linked to fat traits, one map 388 kb upstream the gene SDK1, one map inside PPP3CA gene, one inside SCPEP1 gene and the last is not mapped in the porcine genome yet. Associations here presented indicate a moderate effect of these genes on IMF. In particular, PPP3CA, that is involved in the oxidative metabolism of skeletal muscle, could be considerated as an interesting candidate gene for IMF content in pigs. However, further studies are needed to clarify the role of these genes on the physiological processes involved in IMF regulation. These results may be useful to control this trait that is important in terms of nutritional, technological and organoleptic characteristics of fresh meat and processed products.

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.