A muscle transcriptome analysis identifies positional candidate genes for a complex trait in pig

Skeletal muscle transcriptome in healthy aging

Age-associated changes in gene expression in skeletal muscle of healthy individuals reflect accumulation of damage and compensatory adaptations to preserve tissue integrity. To characterize these changes, RNA was extracted and sequenced from muscle biopsies collected from 53 healthy individuals (22-83 years old) of the GESTALT study of the National Institute on Aging-NIH. Expression levels of 57,205 protein-coding and non-coding RNAs were studied as a function of aging by linear and negative binomial regression models. From both models, 1134 RNAs changed significantly with age. The most differentially abundant mRNAs encoded proteins implicated in several age-related processes, including cellular senescence, insulin signaling, and myogenesis. Specific mRNA isoforms that changed significantly with age in skeletal muscle were enriched for proteins involved in oxidative phosphorylation and adipogenesis. Our study establishes a detailed framework of the global transcriptome and mRNA isoforms that govern muscle damage and homeostasis with age.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

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.

Comparative transcriptomic analysis reveals beneficial effect of dietary mulberry leaves on the muscle quality of finishing pigs

Background

The aim of this study was to investigate the effect of dietary mulberry leaves on the transcriptome profiles of finishing pigs. RNA‐Seq was used to identify differentially expressed genes (DEGs) in the longissimus dorsi of 56 pigs fed either a traditional diet or diets supplemented with 3%, 6% or 9% mulberry leaf powder, and both gene ontology (GO) function and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses were performed. Furthermore, protein–protein interaction (PPI) network and the subnet module analysis were used to identify genes with beneficial potential, and quantitative real‐time polymerase chain reaction (qRT‐PCR) was used to validate the expression patterns revealed by RNA‐Seq.

Results

Pigs fed with the 6% mulberry diet exhibited greater average daily gain, lower water loss and lower shear force than the control group and yielded 531 DEGs, including 271 and 260 upregulated and downregulated genes, respectively. Function analysis revealed that the DEGs were significantly enriched in functions related to muscle growth and development. Furthermore, several genes (i.e. ACOT4, ECHS1, HACD1, NPR1, ADCY2, MGLL and IRS1) were enriched in a KEGG pathway that was associated with fatty acid metabolism, and in the PPI subnet module, four of eight node genes, namely TNNC1, MYL3, TCAP and TNNT1, were associated with muscle formation and development. The upregulation of these genes, including TNNC1, TNNT1 and MYL3, was confirmed by qRT‐PCR.

Conclusions

Dietary mulberry leaves (6%) may improve the muscle quality of pigs by modulating the expression of several key genes, such as TNNC1, MYL3 and TNNT1.

The study was aimed to explain the effect of the inclusion of mulberry in the diet of pigs on transcriptome profiling. The inclusion of mulberry in the diet might be helpful in muscle formation and development of pigs by modulating the expression levels of three genes including TNNC1, MYL3 and TNNT1

NK, S. S, Fairoze MN, Kaur M, Sharma A, Girdhar Y, M. SR, Devatkal SK, Ahlawat S, Vijh RK, S. MS. Transcriptome profiling of longissimus thoracis muscles identifies highly connected differentially expressed genes in meat type sheep of India

This study describes the muscle transcriptome profile of Bandur breed, a consumer favoured, meat type sheep of India. The transcriptome was compared to the less desirable, unregistered local sheep population, in order to understand the molecular factors related to muscle traits in Indian sheep breeds. Bandur sheep have tender muscles and higher backfat thickness than local sheep. The longissimus thoracis transcriptome profiles of Bandur and local sheep were obtained using RNA sequencing (RNA Seq). The animals were male, non-castrated, with uniform age and reared under similar environment, as well as management conditions. We could identify 568 significantly up-regulated and 538 significantly down-regulated genes in Bandur sheep (p≤0.05). Among these, 181 up-regulated and 142 down-regulated genes in Bandur sheep, with a fold change ≥1.5, were considered for further analysis. Significant Gene Ontology terms for the up-regulated dataset in Bandur sheep included transporter activity, substrate specific transmembrane, lipid and fatty acid binding. The down-regulated activities in Bandur sheep were mainly related to RNA degradation, regulation of ERK1 and ERK2 cascades and innate immune response. The MAPK signaling pathway, Adipocytokine signaling pathway and PPAR signaling pathway were enriched for Bandur sheep. The highly connected genes identified by network analysis were CNOT2, CNOT6, HSPB1, HSPA6, MAP3K14 and PPARD, which may be important regulators of energy metabolism, cellular stress and fatty acid metabolism in the skeletal muscles. These key genes affect the CCR4-NOT complex, PPAR and MAPK signaling pathways. The highly connected genes identified in this study, form interesting candidates for further research on muscle traits in Bandur sheep.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Differential expression and co-expression gene networks reveal candidate biomarkers of boar taint in non-castrated pigs

Boar taint (BT) is an offensive odour or taste observed in pork from a proportion of non-castrated male pigs. Surgical castration is effective in avoiding BT, but animal welfare issues have created an incentive for alternatives such as genomic selection. In order to find candidate biomarkers, gene expression profiles were analysed from tissues of non-castrated pigs grouped by their genetic merit of BT. Differential expression analysis revealed substantial changes with log-transformed fold changes of liver and testis from -3.39 to 2.96 and -7.51 to 3.53, respectively. Co-expression network analysis revealed one module with a correlation of -0.27 in liver and three modules with correlations of 0.31, -0.44 and -0.49 in testis. Differential expression and co-expression analysis revealed candidate biomarkers with varying biological functions: phase I (COQ3, COX6C, CYP2J2, CYP2B6, ACOX2) and phase II metabolism (GSTO1, GSR, FMO3) of skatole and androstenone in liver to steroidgenesis (HSD17B7, HSD17B8, CYP27A1), regulation of steroidgenesis (STARD10, CYB5R3) and GnRH signalling (MAPK3, MAP2K2, MAP3K2) in testis. Overrepresented pathways included "Ribosome", "Protein export" and "Oxidative phosphorylation" in liver and "Steroid hormone biosynthesis" and "Gap junction" in testis. Future work should evaluate the biomarkers in large populations to ensure their usefulness in genomic selection programs.

Transcriptome analysis reveals that constant heat stress modifies the metabolism and structure of the porcine longissimus dorsi skeletal muscle

Exposure to high ambient temperatures is detrimental to pig rearing and porcine meat quality. Deep molecular sequencing allows for genomic characterization of porcine skeletal muscles and helps understand how the genomic landscape may impact meat quality. To this end, we performed mRNA-seq to molecularly dissect the impact of heat stress on porcine skeletal muscles, longissimus dorsi. Sixteen castrated, male DLY pigs [which are crossbreeds between Duroc (D) boars and Landrace (L) × Yorkshire (Y) sows, 79.0 ± 1.5 kg BW] were evenly split into two groups that were subjected to either control (CON) (22 °C; 55 % humidity) or constant heat stress (H30; 30 °C; 55 % humidity) conditions for 21 days. Seventy-eight genes were found to be differentially expressed, of which 37 were up-regulated and 41 were down-regulated owing to constant heat stress. We predicted 5247 unknown genes and 6108 novel transcribed units attributed to alternative splicing (AS) events in the skeletal muscle. Furthermore, 30,761 and 31,360 AS events were observed in the CON and H30 RNA-seq libraries, respectively. The differentially expressed genes in the porcine skeletal muscles were involved in glycolysis, lactate metabolism, lipid metabolism, cellular defense, and stress responses. Additionally, the expression levels of these genes were associated with variations in meat quality between the CON and H30 groups, indicating that heat stress modulated genes crucial to skeletal muscle development and metabolism. Our transcriptomic analysis provides valuable information for understanding the molecular mechanisms governing porcine skeletal muscle development. Such insights may lead to innovative strategies to improve meat quality of pigs under heat stress.

Genomic Footprints in Selected and Unselected Beef Cattle Breeds in Korea

Korean Hanwoo cattle have been subjected to intensive artificial selection over the past four decades to improve meat production traits. Another three cattle varieties very closely related to Hanwoo reside in Korea (Jeju Black and Brindle) and in China (Yanbian). These breeds have not been part of a breeding scheme to improve production traits. Here, we compare the selected Hanwoo against these similar but presumed to be unselected populations to identify genomic regions that have been under recent selection pressure due to the breeding program. Rsb statistics were used to contrast the genomes of Hanwoo versus a pooled sample of the three unselected population (UN). We identified 37 significant SNPs (FDR corrected) in the HW/UN comparison and 21 known protein coding genes were within 1 MB to the identified SNPs. These genes were previously reported to affect traits important for meat production (14 genes), reproduction including mammary gland development (3 genes), coat color (2 genes), and genes affecting behavioral traits in a broader sense (2 genes). We subsequently sequenced (Illumina HiSeq 2000 platform) 10 individuals of the brown Hanwoo and the Chinese Yanbian to identify SNPs within the candidate genomic regions. Based on allele frequency differences, haplotype structures, and literature research, we singled out one non-synonymous SNP in the APP gene (APP: c.569C>T, Ala199Val) and predicted the mutational effect on the protein structure. We found that protein-protein interactions might be impaired due to increased exposed hydrophobic surfaces of the mutated protein. The APP gene has also been reported to affect meat tenderness in pigs and obesity in humans. Meat tenderness has been linked to intramuscular fat content, which is one of the main breeding goals for brown Hanwoo, potentially supporting a causal influence of the herein described nsSNP in the APP gene.

Polymorphisms in the regulatory region of the porcine MYLPF gene are related to meat quality traits in the Large White breed

The MYLPF gene encodes fast myosin regulatory light chain, and is a positional and functional candidate gene for meat quality. The aim of this study was to identify associations between SNPs in the promoter region of the porcine MYLPF gene and meat quality traits. A total of 22 SNPs were identified in a population of crossbred animals (n=86) and based on minor allele frequency and proximity to the transcription start site, five SNPs were genotyped in purebred; Large White (n=98), Duroc (n=99) and Pietrain (n=98) pigs. No associations were observed in the Pietrain breed, while the Duroc breed was almost monomorphic for all SNPs. In the Large White breed SNP g-1314A>G and linked SNPS g.-871T>G, g.-566T>C, g.-403C>G were associated with ultimate pH and driploss (P<0.05). This study identified associations between MYLPF and meat quality and highlights the importance of considering the genetic background within gene-assisted selection programmes.

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.

Growth performance, meat quality traits, and genetic mapping of quantitative trait loci in 3 generations of Japanese quail populations (Coturnix japonica)

The current research was conducted to compare growth, carcass traits, pH, intramuscular collagen (IMC) properties, and genetic bases of IMC and carcasses (breast-muscle weight) of different lines and generations of adult males and females of Japanese quail (Coturnix japonica). Forty-four quails (generation F0), 22 Pharaoh (F-33) meat-type males and 22 Standard (S-22) laying-type females, were crossed to produce the F1 hybrids generation. The F2 generation was created by mating one F1 male with one F1 female, full siblings. The birds, randomly chosen from F0 (22 males and 22 females), F1 (22 males and 22 females), and F2 (84 males and 152 females) were raised to 20 wk of age in collective cages. Quails were fed ad libitum commercial diets. At slaughter, all birds were individually weighed (after a fasting period of 12 h) and dressing yield (without giblets) was calculated. The carcasses were then dissected. Genomic DNA was extracted from all of the blood, and 30 microsatellite markers located on 2 quail chromosomes were genotyped. The F -: 33 quails had higher in vivo and postmortem performances and a higher abdominal fat percentage than those of the egg line. Meat from S -: 22 quails had a slower collagen maturation (hydroxylysylpyridinoline crosslink/collagen) and a higher ultimate pH. The F1 and F2 generations showed an evident sexual dimorphism, and an additional effect could be due to hybrid heterosis evident in F2. Meat from quails of F1 and F2 generations had a lower IMC amount with a higher degree of collagen maturation compared with parental lines. Two statistically significant QTL have been detected on quail chromosome 2 (CJA02): a QTL with an additive effect (0.50) for IMC in the marker bracket GUJ0037 and GUJ0093; a second QTL with additive (1.32) and dominant (1.91) effects for breast-muscle weight in the marker bracket GUJ0084 and GUJ0073. To our knowledge, this is the first report of a QTL associated with breast-muscle weight and IMC in quail and poultry species, respectively.

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.

A genetical genomics approach reveals new candidates and confirms known candidate genes for drip loss in a porcine resource population

In this study lean meat water-holding capacity (WHC) of a Duroc × Pietrain (DuPi) resource population with corresponding genotypes and transcriptomes was investigated using genetical genomics. WHC was characterized by drip loss measured in M. longissimus dorsi. The 60K Illumina SNP chips identified genotypes of 169 F2 DuPi animals. Whole-genome transcriptomes of muscle samples were available for 132 F2 animals using the Affymetrix 24K GeneChip® Porcine Genome Array. Performing genome-wide association studies of transcriptional profiles, which are correlated with phenotypes, allows elucidation of cis- and trans-regulation. Expression levels of 1,228 genes were significantly correlated with drip loss and were further analyzed for enrichment of functional annotation groups as defined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. A hypergeometric gene set enrichment test was performed and revealed glycolysis/glyconeogenesis, pentose phosphate pathway, and pyruvate metabolism as the most promising pathways. For 267 selected transcripts, expression quantitative trait loci (eQTL) analysis was performed and revealed a total of 1,541 significant associations. Because of positional accordance of the gene underlying transcript and the eQTL location, it was possible to identify eight eQTL that can be assumed to be cis-regulated. Comparing the results of gene set enrichment and the eQTL detection tests, molecular networks and potential candidate genes, which seemed to play key roles in the expression of WHC, were detected. The α-1-microglobulin/bikunin precursor (AMBP) gene was assumed to be cis-regulated and was part of the glycolysis pathway. This approach supports the identification of trait-associated SNPs and the further biological understanding of complex traits.

Identification of the novel candidate genes and variants in boar liver tissues with divergent skatole levels using RNA deep sequencing

Boar taint is the unpleasant odour of meat derived from non-castrated male pigs, caused by the accumulation of androstenone and skatole in fat. Skatole is a tryptophan metabolite produced by intestinal bacteria in gut and catabolised in liver. Since boar taint affects consumer's preference, the aim of this study was to perform transcriptome profiling in liver of boars with divergent skatole levels in backfat by using RNA-Seq. The total number of reads produced for each liver sample ranged from 11.8 to 39.0 million. Approximately 448 genes were differentially regulated (p-adjusted <0.05). Among them, 383 genes were up-regulated in higher skatole group and 65 were down-regulated (p<0.01, FC>1.5). Differentially regulated genes in the high skatole liver samples were enriched in metabolic processes such as small molecule biochemistry, protein synthesis, lipid and amino acid metabolism. Pathway analysis identified the remodeling of epithelial adherens junction and TCA cycle as the most dominant pathways which may play important roles in skatole metabolism. Differential gene expression analysis identified candidate genes in ATP synthesis, cytochrome P450, keratin, phosphoglucomutase, isocitrate dehydrogenase and solute carrier family. Additionally, polymorphism and association analysis revealed that mutations in ATP5B, KRT8, PGM1, SLC22A7 and IDH1 genes could be potential markers for skatole levels in boars. Furthermore, expression analysis of exon usage of three genes (ATP5B, KRT8 and PGM1) revealed significant differential expression of exons of these genes in different skatole levels. These polymorphisms and exon expression differences may have impacts on the gene activity ultimately leading to skatole variation and could be used as genetic marker for boar taint related traits. However, further validation is required to confirm the effect of these genetic markers in other pig populations in order to be used in genomic selection against boar taint in pig breeding programs.

The structure of a gene co-expression network reveals biological functions underlying eQTLs

What are the commonalities between genes, whose expression level is partially controlled by eQTL, especially with regard to biological functions? Moreover, how are these genes related to a phenotype of interest? These issues are particularly difficult to address when the genome annotation is incomplete, as is the case for mammalian species. Moreover, the direct link between gene expression and a phenotype of interest may be weak, and thus difficult to handle. In this framework, the use of a co-expression network has proven useful: it is a robust approach for modeling a complex system of genetic regulations, and to infer knowledge for yet unknown genes. In this article, a case study was conducted with a mammalian species. It showed that the use of a co-expression network based on partial correlation, combined with a relevant clustering of nodes, leads to an enrichment of biological functions of around 83%. Moreover, the use of a spatial statistics approach allowed us to superimpose additional information related to a phenotype; this lead to highlighting specific genes or gene clusters that are related to the network structure and the phenotype. Three main results are worth noting: first, key genes were highlighted as a potential focus for forthcoming biological experiments; second, a set of biological functions, which support a list of genes under partial eQTL control, was set up by an overview of the global structure of the gene expression network; third, pH was found correlated with gene clusters, and then with related biological functions, as a result of a spatial analysis of the network topology.

Transcriptional profiling and miRNA-dependent regulatory network analysis of longissimus dorsi muscle during prenatal and adult stages in two distinct pig breeds

MicroRNAs (miRNAs) and mRNAs establish a complex regulatory network influencing diverse biological pathways including muscle development and growth. Elucidating miRNA-dependent regulatory networks involved in muscle development could provide additional insights into muscle traits largely predefined during prenatal development. The present study aimed to determine differentially expressed transcripts and functional miRNA-mRNA relationships associated with different stages of skeletal muscle development in two pig breeds, German Landrace and Pietrain, distinct in muscle characteristics. A comparative transcriptional profiling of longissimus dorsi muscle tissues from fetuses at 35, 63 and 91 days post-conception as well as adult pigs (180 days postnatum) was performed using the Affymetrix GeneChip porcine genome microarray. Differential expression patterns were identified to be associated with muscularly developmental stages and breed types. The integration of miRNA expression data and ingenuity pathways analysis (ipa) pathway analysis revealed several miRNA-dependent regulatory networks related to muscle growth and development. The present results provide insights into muscle biology for further improvement of porcine meat quality.

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.

Molecular characterization and expression analysis of NDUFS4 gene in m. longissimus dorsi of Laiwu pig (Sus scrofa)

To study the molecular basis of intramuscular fat (IMF) deposition, suppression subtractive hybridization was used to investigate the differences in gene expression between m. longissimus dorsi (LD) of high IMF Laiwu pig group and low IMF Laiwu pig group. From two specific subtractive cDNA libraries, the expression-upregulated clone HL-27 was selected by reverse Northern high-density blot, and then identified to be pig mitochondrial NADH dehydrogenase (ubiquinone) Fe-S protein 4 (NDUFS4). Pig NDUFS4 full-length cDNA was cloned by RACE, and contains a 528 bp-open reading frame (ORF) encoding 175 amino acid residues. The derived amino acid sequence of NDUFS4 is well conserved compared with NDUFS4 of various species with higher degree of sequence similarity with other mammalian (86.3-92.6 %) than amphibian, aves, and fishes (70.2-81.1 %), and contains one N-linked glycosylation site, one O-linked glycosylation site, seven Ser phosphorylation sites and five Thr phosphorylation sites. A-G mutation was found at nt 122 site of ORF between Laiwu pig and Large White, which results in the K-R mutation at 41 site of protein sequence. Real-time PCR analysis indicated that the level of NDUFS4 mRNA expression was higher in high IMF Laiwu pig group than in low IMF Laiwu pig group, and in Laiwu pig than in Large White. The tissue expression of the pig NDUFS4 gene showed a tissue-specific pattern: highly expressed in LD muscle, spleen and kidney, but hardly expressed in lung, stomach and large intestine. The possible role of NDUFS4 and its relation to IMF deposition are discussed.

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.

Genetic variability of transcript abundance in pig peri-mortem skeletal muscle: eQTL localized genes involved in stress response, cell death, muscle disorders and metabolism

Background

The genetics of transcript-level variation is an exciting field that has recently given rise to many studies. Genetical genomics studies have mainly focused on cell lines, blood cells or adipose tissues, from human clinical samples or mice inbred lines. Few eQTL studies have focused on animal tissues sampled from outbred populations to reflect natural genetic variation of gene expression levels in animals. In this work, we analyzed gene expression in a whole tissue, pig skeletal muscle sampled from individuals from a half sib F2 family shortly after slaughtering.

Results

QTL detection on transcriptome measurements was performed on a family structured population. The analysis identified 335 eQTLs affecting the expression of 272 transcripts. The ontologic annotation of these eQTLs revealed an over-representation of genes encoding proteins involved in processes that are expected to be induced during muscle development and metabolism, cell morphology, assembly and organization and also in stress response and apoptosis. A gene functional network approach was used to evidence existing biological relationships between all the genes whose expression levels are influenced by eQTLs. eQTLs localization revealed a significant clustered organization of about half the genes located on segments of chromosome 1, 2, 10, 13, 16, and 18. Finally, the combined expression and genetic approaches pointed to putative cis-drivers of gene expression programs in skeletal muscle as COQ4 (SSC1), LOC100513192 (SSC18) where both the gene transcription unit and the eQTL affecting its expression level were shown to be localized in the same genomic region. This suggests cis-causing genetic polymorphims affecting gene expression levels, with (e.g. COQ4) or without (e.g. LOC100513192) potential pleiotropic effects that affect the expression of other genes (cluster of trans-eQTLs).

Conclusion

Genetic analysis of transcription levels revealed dependence among molecular phenotypes as being affected by variation at the same loci. We observed the genetic variation of molecular phenotypes in a specific situation of cellular stress thus contributing to a better description of muscle physiologic response. In turn, this suggests that large amounts of genetic variation, mediated through transcriptional networks, can drive transient cell response phenotypes and contribute to organismal adaptative potential.

Transcriptional profiling during foetal skeletal muscle development of Piau and Yorkshire-Landrace cross-bred pigs

Skeletal muscle development is a complex process involving the coordinated expression of thousands of genes. The aim of this study was to identify differentially expressed genes in longissimus dorsi (LD) muscle of pigs at 40 and 70 days (d) of gestation (developmental stages encompassing primary and secondary fibre formation) in Yorkshire-Landrace (YL) cross-bred pigs and Piau pigs (a naturalized Brazilian breed), which are two breed types that differ in muscularity. Foetuses were obtained from gilts at each gestational age (n = 3 YL; n = 4 Piau), and transcriptional profiling was performed using the Pigoligoarray microarray containing 20 400 oligonucleotides. A total of 486 oligonucleotides were differentially expressed (fold change (FC) ≥ 1.5; false discovery rate (FDR) ≤ 0.05) between 40 and 70 d gestation in either YL or Piau pigs, and a total of 1300 oligonucleotides were differentially expressed (FC ≥ 1.5; FDR ≤ 0.05) between YL and Piau pigs at either age. Gene ontology annotation and pathway analyses determined functional classifications for differentially expressed genes and revealed breed type-specific developmental expression patterns. Thirteen genes were selected for confirmation by qRT-PCR analyses, and expression patterns for most of these genes were confirmed, providing further insight into the roles of these genes in pig muscle development. This study revealed both developmental and breed type-specific patterns of gene expression in foetal pig skeletal muscle, including genes not previously associated with myogenesis. This information will contribute to future pig genetic improvement efforts.

Identification of differentially expressed genes in chickens differing in muscle glycogen content and meat quality

Background

The processing ability of poultry meat is highly related to its ultimate pH, the latter being mainly determined by the amount of glycogen in the muscle at death. The genetic determinism of glycogen and related meat quality traits has been established in the chicken but the molecular mechanisms involved in variations in these traits remain to be fully described. In this study, Chicken Genome Arrays (20 K) were used to compare muscle gene expression profiles of chickens from Fat (F) and Lean (L) lines that exhibited high and low muscle glycogen content, respectively, and of individuals exhibiting extremely high (G+) or low (G-) muscle glycogen content originating from the F₂ cross between the Fat and Lean lines. Real-time RT-PCR was subsequently performed to validate the differential expression of genes either selected from the microarray analysis or whose function in regulating glycogen metabolism was well known.

Results

Among the genes found to be expressed in chicken P. major muscle, 197 and 254 transcripts appeared to be differentially expressed on microarrays for the F vs. L and the G+ vs. G- comparisons, respectively. Some involved particularly in lipid and carbohydrate metabolism were selected for further validation studies by real-time RT-PCR. We confirmed that, as in mammals, the down-regulation of CEBPB and RGS2 coincides with a decrease in peripheral adiposity in the chicken, but these genes are also suggested to affect muscle glycogen turnover through their role in the cAMP-dependent signalling pathway. Several other genes were suggested to have roles in the regulation of glycogen storage in chicken muscle. PDK4 may act as a glycogen sensor in muscle, UGDH may compete for glycogen synthesis by using UDP-glucose for glucoronidation, and PRKAB1, PRKAG2, and PHKD may impact on glycogen turnover in muscle, through AMP-activated signalling pathways.

Conclusions

This study is the first stage in the understanding of molecular mechanisms underlying variations in poultry meat quality. Large scale analyses are now required to validate the role of the genes identified and ultimately to find molecular markers that can be used for selection or to optimize rearing practices.

Gene expression profiling of skeletal muscle of nursing piglets

To gain insight into the regulation mechanism associated with the rapid gain in skeletal muscle during neonatal period, gene expression profiles of skeletal muscle of nursing pigs was investigated using Affymetrix Porcine GeneChip. A total of 1094 transcripts were detected as differential expression over time course tested (p<0.01, q<0.05). With combinative use of partitioning around medoid and hierarchical clustering, three clusters of transcripts with distinct temporal expression were defined. Gene functional categories and pathways, particularly involved in cell signaling, cell cycle, cell adhesion, ECM-receptor interaction, glycolysis, protein synthesis and degradation, and intracellular transport, were identified. Moreover, we showed 49 of the differentially expressed genes within published QTL regions or with marked deletion effects. Our study demonstrates previously uncharacterized changes in transcription accompanying early postnatal growth of skeletal muscle of pigs. It has highlighted potential cascades and important candidates for further investigation on controlling of postnatal muscle growth.

Transcriptome analysis to identify differential gene expression affecting meat quality in heavy Italian pigs

Suppressive subtractive hybridization (SSH) was used to analyse the muscle transcriptome and identify genes affecting meat quality within an Italian pig population of Large White and Landrace purebred individuals. Seven phenotypes were recorded at slaughter: dorsal fat thickness, ham fat thickness, ham fat coverage, muscle compactness, marbling, meat colour and colour uniformity. Two subtractive libraries were created from longissimus dorsi tissue of selected pigs with extreme phenotypes for meat quality. Eighty-four differentially expressed ESTs were identified, which showed homology to expressed pig sequences and/or to genomic pig sequences produced within the pig genome project. Sixty-eight sequences were mapped on the pig genome, and most of these sequences co-localized with the same chromosomal positions as QTLs that have been previously identified for meat quality. Thirty sequences, including eight matching known genes previously related to muscle metabolic pathways, were selected to statistically validate their differential expression. Association analysis and t-test results indicated that 28 ESTs of the 30 analysed were associated with phenotypes investigated here and have significant differential expression levels (P≤ 0.05) between the two tails of the phenotypic distribution.

Identification of suitable reference genes for gene expression analysis of pork meat quality and analysis of candidate genes associated with the trait drip loss

The aim of this study was to identify a set of stably expressed endogenous control genes for quantitative PCR analysis of mRNA expression in the porcine LTL muscle and to subsequently perform expression analysis of potential candidate genes associated with drip loss. Expression stability of seven commonly used reference genes was examined in n=60 pigs from three independent populations of different genetic backgrounds. The genes examined were: ACTB, ATP5G1, B2M, GPX1, RPL4, TBP and YWHAZ. GeNorm analysis of expression stability identified B2M, RPL4 and TBP as consistently stable in each breed examined. Analysis of meat samples divergent for water holding capacity identified positive and negative associations between drip loss and gene expression using B2M, RPL4 and TBP as endogenous controls. Specifically, expression of COL1A1 increased significantly with increasing drip loss while expression of CAST decreased significantly with increasing drip loss. This study therefore indicates the use of B2M, RPL4 and TBP as suitable endogenous controls for gene expression analysis of the porcine LTL muscle. Further study is recommended to identify the detailed roles of COL1A1 and CAST with respect to the development of drip loss.

Additive effects of 19 porcine SNPs on growth rate, meat content and selection index

A total of 306 boars (108 Large White and 198 Landrace) were genotyped for 52 candidate SNPs to determine which of the polymorphisms influence growth rate, meat content and selection index. The effects of SNPs were estimated by a mixed linear model including a random additive polygenic animal effect, fixed effects of SNPs including additive, and pairwise additive-by-additive epistases, year*season of birth, breed and RYR1 genotype. In order to estimate all possible pairwise SNP combinations without overparameterising the model a stochastic approach was adopted. A total of 1 350 replications of the model were generated, each containing five randomly selected SNPs. The final estimates of the fixed effects of the model equaled an average out of the replications. The hypothesis of a nonzero effect of SNP was tested by the Wald test. Among 4 257 estimates calculated, many significant (P<0.01), but mostly minor effects (below 1 phenotypic standard deviation) were recorded. The selected SNPs will be further investigated to determine which may be used in MAS.

Using microarrays to identify positional candidate genes for QTL: the case study of ACTH response in pigs

BACKGROUND

Microarray studies can supplement QTL studies by suggesting potential candidate genes in the QTL regions, which by themselves are too large to provide a limited selection of candidate genes. Here we provide a case study where we explore ways to integrate QTL data and microarray data for the pig, which has only a partial genome sequence. We outline various procedures to localize differentially expressed genes on the pig genome and link this with information on published QTL. The starting point is a set of 237 differentially expressed cDNA clones in adrenal tissue from two pig breeds, before and after treatment with adrenocorticotropic hormone (ACTH).

RESULTS

Different approaches to localize the differentially expressed (DE) genes to the pig genome showed different levels of success and a clear lack of concordance for some genes between the various approaches. For a focused analysis on 12 genes, overlapping QTL from the public domain were presented. Also, differentially expressed genes underlying QTL for ACTH response were described. Using the latest version of the draft sequence, the differentially expressed genes were mapped to the pig genome. This enabled co-location of DE genes and previously studied QTL regions, but the draft genome sequence is still incomplete and will contain many errors. A further step to explore links between DE genes and QTL at the pathway level was largely unsuccessful due to the lack of annotation of the pig genome. This could be improved by further comparative mapping analyses but this would be time consuming.

CONCLUSION

This paper provides a case study for the integration of QTL data and microarray data for a species with limited genome sequence information and annotation. The results illustrate the challenges that must be addressed but also provide a roadmap for future work that is applicable to other non-model species.