- Research article
- Open Access
Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties
BMC Genomics volume 14, Article number: 533 (2013)
Physiological processes aiding the conversion of muscle to meat involve many genes associated with muscle structure and metabolic processes. MicroRNAs regulate networks of genes to orchestrate cellular functions, in turn regulating phenotypes.
We applied weighted gene co-expression network analysis to identify co-expression modules that correlated to meat quality phenotypes and were highly enriched for genes involved in glucose metabolism, response to wounding, mitochondrial ribosome, mitochondrion, and extracellular matrix. Negative correlation of miRNA with mRNA and target prediction were used to select transcripts out of the modules of trait-associated mRNAs to further identify those genes that are correlated with post mortem traits.
Porcine muscle co-expression transcript networks that correlated to post mortem traits were identified. The integration of miRNA and mRNA expression analyses, as well as network analysis, enabled us to interpret the differentially-regulated genes from a systems perspective. Linking co-expression networks of transcripts and hierarchically organized pairs of miRNAs and mRNAs to meat properties yields new insight into several biological pathways underlying phenotype differences. These pathways may also be diagnostic for many myopathies, which are accompanied by deficient nutrient and oxygen supply of muscle fibers.
Muscle is the major energy consuming and storage organ. An imbalance of nutrients, energy, and oxygen supply-and-demand in muscle cells is evident following cardiac muscle or skeletal muscle attack, injury, or damage. The consequences of these imbalances depend on muscle structure and metabolism and, thus, the muscle’s entire complement of proteins and their expression patterns. Similar changes, i.e., termination of nutrient and energy supplies and anoxia, also occur in muscle cells post mortem. Indeed, these changes underlie the conversion of muscle to meat in food production. The physiological processes accompanying the change of muscle to meat involve expression of many genes associated with muscle structure and metabolic features [1, 2]. Genes active in the muscle could therefore potentially have pathogenetic effects by disturbing muscle energy and oxygen homeostasis in vivo, as well as conferring traits related to meat quality post mortem.
Meat quality is complex and is affected by genetic and environmental factors as well as slaughtering procedures . The conversion of muscle to meat is important not only as an economic factor in pork production, but also because these events mimic pathological processes associated with muscle injury or damage in humans. Post mortem traits for meat quality and carcass are influenced by a complex network of gene interactions in muscle; therefore, elucidating the relationships between genes and how these genes, in turn, influence meat quality and carcass traits is critical for developing a comprehensive understanding of the muscle to meat conversion as well as muscle pathologic processes including muscle atrophy, dystrophy, and hypoxia. Additionally, pigs share many genomic and physiological similarities with humans and, therefore, provide a good model to study the genetic determination of complex traits and as a biomedical model [4, 5].
Recent advances in functional genomic screening, which can help determine molecular processes underlying phenotypic differences [6–8], have identified roles for microRNAs (miRNAs) in regulation of myogenesis [9–11] and adipogenesis [12–14]. miRNAs are small, non-coding RNA molecules of approximately 22 nucleotides. The primary miRNA transcript has a stem-loop structure that is recognized and cleaved via RNA processing enzymes to produce a double-stranded duplex. miRNAs target mRNA transcripts via base-pair complementarity, typically in the 3′ untranslated region [15, 16], but also in the coding sequence . This targeting can induce transcript cleavage, degradation, destabilization, or repression of translation, thereby modulating protein levels. It has been recently shown that reduction of transcript level account for most of the regulatory, repressive effects of miRNAs . Target genes that are regulated by miRNAs through degradation of their respective transcripts consequently show negative correlation of their mRNA with the miRNA regulator. Moreover, one miRNA can target several — even hundreds of — genes. Therefore, a unique approach for identifying miRNA-mRNA regulatory modules was recently introduced, whereby paired miRNA-mRNA expression profiles were constructed to predict putative target genes of miRNAs [19, 20].
Many studies used the network analysis for dissecting the complex traits [21, 22]. Weighted gene co-expression network analysis (WGCNA)  has been successfully applied in a variety of different settings [24–28]. WGCNA groups genes into modules based on their co-expression across a set of samples and finally relates these modules to the traits of interest in order to elucidate relevant modules or genes.
In order to identify groups of co-expressed genes (mRNAs) and on the hierarchically superior level of miRNAs that are correlated with organismal traits related to carcass and meat quality, we first applied weighted gene co-expression network analysis (WGCNA) and subsequently we adapted the paired expression profile approach. We identified co-expression networks regulated by miRNA after filtering of negatively-correlated miRNA-mRNA pairs and predicting target genes. The integration of miRNA and mRNA expression analyses as well as network analysis enabled us to interpret the differentially-regulated genes from a systems perspective, yielding new insight into several biological pathways underlying phenotypic differences.
Meat quality and carcass traits phenotypes
Elucidating the relationships between genes and how these genes, in turn, influence muscle metabolic and structural properties is critical for developing a comprehensive understanding of the muscle to meat process as well as muscle pathologic and regenerative processes related to muscle atrophy, dystrophy and hypoxia. In total, 207 performance-tested crossbred pig [PI × (DL × DE)] samples were used to investigate meat quality and carcass traits. Descriptions of 7 carcass traits and 13 meat quality traits, as well as means and standard deviations, analysed in this study are shown in Table 1. High correlation coefficients were found between the same biochemical and biophysical parameters measured at different positions and at different time points post mortem, like fatness traits or pH. The other cluster covers the traits drip loss, protein content, and conductivity. This cluster was negatively correlated to pH or fatness (Additional file 1: Figure S1).
Gene co-expression network construction for mRNA
To investigate the role of transcriptional networks in muscle, we performed a weighted gene co-expression network analysis (WGCNA) using expression data from M. longissimus dorsi necropsies of the 207 performance-tested crossbred pigs [PI × (DL × DE)]. Expression analysis using GeneChip Porcine Genome Arrays (Affymetrix) containing 24,123 probe sets identified 11,191 probe sets with consistent expression according to MAS5 analysis; these were used for further analysis. Using WGCNA, residuals derived from the mixed-model analysis of expression levels of 11,191 probe sets were used for constructing the muscle transcriptional network. WGCNA grouped genes into 22 modules based on patterns of co-expression. Each module was labelled with a unique color identifier and was characterized for enrichment of genes of specific gene ontology (GO) categories (Table 2). To represent the gene expression profiles of the highly correlated genes inside a given module, we used the first principal component, which is referred to as the module eigengene (ME). We tested each ME for correlation with meat and carcass traits.
Module-trait associations of mRNA
Sets of genes (modules) with common expression patterns that were associated with particular traits were identified based on the correlation between ME and organismal phenotype. We identified five modules that significantly associated with meat quality. Modules dark-turquoise and orange were correlated positively to pH traits and negatively to drip loss (ME[dark-turquoise]: pH24MLD r = 0.34, p = 5.3 × 10-7, DL r = -0.19, p = 5.6 × 10-3; ME[orange]: pH24MLD r = 0.32, p = 3.7 × 10-6, DL r = -0.31, p = 5.8 ×10-6) (Figure 1). Module dark-turquoise (31 annotated genes) was highly enriched for genes belonging to the cluster “glucose metabolic process” (GO: 0006006) and the KEGG-pathway “insulin signaling” with an enrichment score (ES) of 2.65. Module orange (26 annotated genes) was enriched for transcripts of the functional annotation clusters “response to wounding”, “defense response” and “inflammatory response” (ES = 2.42). Modules red, black, and tan were correlated negatively to pH traits and positively to drip loss (ME[red]: pH45MLD r = -0.22, p = 1.8 × 10-3, DL r = 0.20, p = 3.9 × 10-3; ME[black]: pH45MLD r = -.23, p = 8.8 × 10-4, DL r = 0.16, p = 1.8 × 10-2; ME[tan]: pH45MLD r = -0.22, p = 1.4 × 10-3, DL r = 0.19, p = 6.9 × 10-3) (Figure 1). Further, modules red (315 annotated genes), black (436 annotated genes), and tan (154 annotated genes) were enriched for genes of the top functional annotation clusters of “mitochondrial ribosome”, “mitochondrion”, and “extracellular matrix” with ES of 10.23, 15.15, and 27.05, respectively. Only one module (ME[dark-orange]) showed association with traits related to fatness (Figure 1).
Co-expression networks and module-trait associations for miRNA
Transcriptional networks of muscle microRNAs were studied with a WGCNA using miRNA expression data on M. longissimus dorsi from the same animals as above. Residuals, derived from the analysis of expression levels after correction for systematic effects according to the mixed model, were used for constructing the muscle miRNA transcriptional networks, i.e., modules. Nine modules were identified (Figure 2). Only 2 modules were associated with meat quality at a significance level of p < 0.05. Module purple was correlated positively to LF24MLD at p = 0.03 and negatively to pH24MLD and IMP24MLD at p = 0.04 and p = 0.03, respectively. Module purple consisted of 8 miRNA families (miR-17, miR-30. miR221, miR-185, miR-324, miR362, miR-500, and miR-542). Module blue was positively correlated to pH45MLD and IMP24MLD at p = 0.02 and p = 0.01, respectively. Module blue comprised 29 miRNA families (let-7, miR-15, miR-17, miR-31, miR-95, miR-103, miR-105, miR-122, miR-124, miR-130, miR-138, miR-154, miR-184, miR-185, miR-197, miR-202, miR-204, miR-212, miR-214, miR-320, miR-326, miR-335, miR-346, miR-383, miR-467, miR-491, miR-744, miR-1224, and miR-1296).
Individual miRNA expression profiles and correlated traits
In addition to miRNA modules, the expression of 675 individual miRNAs probe sets, which corresponded to 513 unique sequences belonging to 159 miRNA families, were profiled and examined for correlation with meat and carcass traits. In total, 225 miRNA-phenotype pairs revealed correlations at p < 0.01. Table 3 lists the top 20 miRNAs with highly significant correlations to phenotypes (p < 0.005).
Endogenous correlation of expression profiles between miRNA and mRNA
We performed pairwise correlation coefficient analysis to evaluate association of expression levels between 675 miRNA probe sets and 11,191 mRNA probe sets. Among the 7,553,925 Pearson correlation coefficients, we detected significant correlation in 5,933 miRNA-mRNA pairs at p-values ≤ 8.97 × 10-5 (FDR = 0.1). The 5,933 pairs comprised 408 miRNA probe sets belonging to 128 miRNA families that were correlated with 2,296 mRNA probe sets. Of these 5,933 pairs, 4,005 and 1,928 pairs showed positive and negative correlations, respectively. Positive correlations tended to be more dramatic than negative correlations; the correlation between miR-122 and VTN was the most significant (FDR = 4.33x10-11). The most significant negative correlation was between miR-154 and LOC387820 (FDR = 1.1x10-5). The most frequently involved miRNA family was miRNA-221, which was correlated with 616 mRNAs. In total 96 miRNA families showed significantly negative correlations with groups of up to 253 genes. We evaluated GO classification for each miRNA-correlated gene set of more than 50 genes (Table 4). The most striking findings were from gene sets that were negatively correlated with miR-23, miR-30, miR17, miR154, and miR-132. For miR-23 and miR-17 the set of negatively correlated genes was highly enriched for genes belonging to the clusters “translation” (GO:0006412) and “translational elongation” (GO:0006414). The set of genes negatively correlated with miR-30 was enriched for “cytoskeletal protein binding” (GO:0008092). The set of genes negatively correlated with miR-154 was enriched for “threonine-type peptidase activity” (GO:0070003), while miR132-correlated genes were enriched for “proteasome complex” (GO:0000502).
Integration of negative correlation of miRNA and mRNA with module-trait association and target prediction
A total of 1,928 pairs of miRNAs and mRNAs that showed negative correlations at p-values ≤ 8.97 × 10-5 (FDR = 0.1) belonged to 1,073 mRNA probe sets (929 gene) and 264 miRNA probe sets (96 miRNA families). Of these, 286 pairs were assigned to modules dark-turquoise, red, black, and tan, which showed correlation with traits related to meat quality (Figure 3). However, no genes in module orange were negatively correlated with miRNA at FDR < 0.1. Only one gene (CREM) in the module dark-turquoise was negatively correlated with miR-153 (r = -0.34 p = 2.11x10-06 FDR = 0.02). In module black, 69 genes were negatively correlated with 21 miRNA families, totaling 140 miRNA-mRNA pairs. In total, 101 and 43 pairs of miRNAs and mRNAs were identified in modules red (22 miRNAs and 52 genes) and tan (17 miRNAs and 26 genes). Known genes belonging to modules dark-turquoise, red, black, and tan and their negatively-correlated miRNAs are shown in Table 5. Out of 1,928 pairs of miRNAs and mRNAs that showed negative correlations at p-values ≤ 8.97 × 10-5 (FDR = 0.1), 62 pairs were assigned to modules blue and purple, which correlated with meat quality (Figure 4). These 62 pairs of miRNAs and mRNAs belonged to 14 miRNA families and 40 genes (Table 6).
Further, TargetScan and RNAhybrid were used to scan miRNA and mRNA sequences (porcine RefSeq), to obtain additional evidence for their functional link; these sequences corresponded to 1,928 pairs of negatively-correlated miRNAs and mRNAs. In total, 474 pairs of miRNA and mRNA were confirmed by either of the two in silico prediction methods: 331, 195, and 32 miRNA-mRNA pairs were predicted by RNAhybrid, TargetScan, or both methods, respectively. The 474 miRNA-mRNA pairs covered 121 probe sets of miRNAs (65 miRNA families) and 331 targets probe sets (297 genes). When focusing on genes of the trait-correlated modules dark-turquoise, orange, red, black, and tan, 73 out 286 pairs of miRNAs and mRNAs were confirmed with at least one in silico method. These 73 pairs comprised 26 miRNA families and 51 genes (Tables 5 and 6 in bold).
Here, we present an integrative approach to identify transcriptomic differences that may contribute to variation of the kinetics of metabolic processes under diminished oxygen and nutrition supply that is evident during muscle conversion to meat. The speed and extend of the switch from aerobic to anaerobic ATP production, until final total failure of energy production, and of protein degradation processes largely affect meat quality . In order to identify functional networks of genes contributing to these processes an approach was used based on a multi-level integration of weighted gene co-expression network analysis (WGCNA) of mRNA and miRNA with mRNA-miRNA pair correlation and miRNA target prediction.
mRNA abundance and co-expression networks linked to muscle and meat properties
We used transcriptional network analysis to identify co-expression modules (dark-turquoise, orange, red, black, and tan) that correlated to meat quality phenotype. These modules were highly enriched for genes involved in “extracellular matrix”, “glucose metabolic process”, and “mitochondrion” (“oxidative phosphorylation” KEGG_PATHWAY); i.e. processes affecting structural and metabolic properties.
A dominant role of mitochondria is the production of ATP by oxidative phosphorylation that depends on oxygen supply. When oxygen is limited (post mortem or during prolonged vigorous exercise) the glucose metabolism occurs by anaerobic respiration, a process that is independent of the mitochondria. A shift from aerobic to anaerobic metabolism - favouring the production of lactic acid - results in a pH decline post mortem and thereby influence the meat quality . So the biological process of mitochondria as well as the way of glucose metabolism play a significant role in the muscle cell and finally impact on meat quality. Indeed, mitochondrial dysfunction resulting in decreased cellular energy production is also responsible for a variety of human myopathies and cardiomyopathies [30–33].
Genes assigned to the GO category of “extracellular matrix” encode proteins belonging to the myofibrillar scaffold. The characteristics of the myofibrillar scaffold and the kinetics of their fragmentation were associated with tenderness and water-holding capacity of meat. In particular, the proteolysis of muscle proteins affects the shrinkage of myofibrils, the development of pores in the cell membranes, so called drip channels, and the non-covalent binding of water molecules [34, 35]. Collagens are major constituents of the extracellular matrix (ECM). In our study there were many collagen genes that are reported to be correlated with various muscle disorders [36–39]. For example, collagen type VI (COLVI), an important component of skeletal muscle ECM, is involved in maintaining tissue integrity . Col6a1-/- mice show a complete absence of collagen VI chains and display a myopathic phenotype, abnormal mitochondria, and increased apoptosis of muscle fibres [41, 42].
miRNA abundance and co-expression networks linked to muscle and meat properties
Our finding of the relevance of mitochondrial metabolic pathways, including oxidative phosphorylation, and muscle structural protein composition to post mortem processes affecting meat quality is in line with our previous results obtained in other populations [1, 6, 7]. Additionally, while the previous studies focused on trait-associated mRNA expression, here another hierarchical level in the regulatory network relevant to processes occurring under conditions of insufficient oxygen, energy, and nutrient supplies is provided. In fact, miRNA was integrated into this study as a regulator molecule of muscle transcripts. miRNAs with identical seed sequences (the same family)  or that are closely located on the same chromosome (the same cluster) [44, 45] have similar expression trends. This was confirmed in our study, where most of the modules consisted of the same families of miRNAs or miRNAs located on the same chromosome. In this study, marginal association of miRNA co-expression modules to organismal traits was found compared to mRNA co-expression modules. On the one hand, co-expressed miRNAs of the same family or cluster might not regulate the same trait. On the other hand, this may be caused by indirect regulation of organismal traits by miRNAs via their effect on mRNA transcripts. Accordingly, individual miRNA correlation to phenotypes was also considered. Recent studies have revealed key roles for miRNAs in the regulation of skeletal muscle differentiation, and changes in miRNA expression are associated with various skeletal muscle disorders [46–48]. In this study, several miRNAs were correlated with carcass and meat quality traits. This includes miR-221, previously identified in studies of myotube maturation and in the maintenance of the myofibrillar organization  and found to contribute to muscle pathogenetic mechanisms . Interestingly, miR-133, which showed highest correlation with pH24MLD and pH24MSM, has been widely studied for roles in the regulation of skeletal muscle development, including in proliferation and myogenesis [10, 51] as well as muscle disorders . Recently, a study reported that mice with genetic deletions of miR-133a-1 and miR-133a-2 developed adult-onset centronuclear myopathy in type II (fast-twitch) myofibres, which was accompanied by impaired mitochondrial function, fast-to-slow myofibre conversion, and disarray of muscle triads . These are changes of muscle structure and metabolism that also impact meat quality. In addition to its well established role in translation (Table 4), miR-23a was also recently identified as a key regulator of skeletal muscle differentiation and is predicted to target multiple adult fast myosin heavy chain (Myh) genes, including Myh1, 2, and 4. For fat traits, miR-103 and miR-107 were highly correlated (Table 3). This is consistent with previous reports of miR-103 being involved in adipogenesis, lipid metabolism, and adipocyte differentiation [54, 55] and of miR-103/107 being involved in glucose homeostasis and insulin sensitivity .
In this study we showed for the first time more complex correlations among miRNAs and between miRNAs and post mortem organismal phenotypes in swine, while also confirming previous studies in human and mouse muscle as well as C2C12 myoblasts.
Links between miRNA and mRNA that relate to muscle and meat properties
This study also sought to evaluate to what extent the co-expression modules of trait-associated mRNAs are themselves regulated by miRNAs. A regulatory link between miRNA and mRNA and a functional link to the organismal phenotype was suspected if (1) the mRNA belonged to either one of the co-expression modules associated with the traits (i.e., MEdark-turquoise, MEorange, MEred, and MEtan), (2) mRNA abundance was significantly negatively correlated with its miRNA regulator, and (3) the mRNA was predicted to be a target gene of the respective miRNA. Therefore, RefSeqs of the genes with present calls from the 3′-IVT-Affymetrix arrays were explored to predict the targets of miRNAs by either seed sequence complementarity  or by thermodynamics-based modeling of RNA:RNA duplex interactions . Currently, no publicly-accessible database covers porcine miRNAs and their predicted target genes. Moreover, annotation of porcine genes is not yet finalized. Accordingly, the target predictions should be interpreted with caution.
Interestingly, no miRNA regulator was identified by negative correlation or target prediction in the module orange, which was enriched for genes related to “response to wounding”, “inflammatory response”, and “defense response”. Genes assigned to biofunctions related to response to exogenes stimuli, change their transcription rate immediately due to many factors. As suggested by finding no correlated miRNAs – their regulation of transcription may occur without major involvement of miRNA. However, many genes in the module orange were previously confirmed as transcriptional regulators in myogenesis or were located in QTL regions for muscle fiber traits like BTG2, EGR1, ANKRDS1 and FOS [59, 60]. Interestingly, the genes in module orange like Egr1, FOS and JUN that are associated with oxidative stress response were found upregulated in muscle in response to mechanical ventilation and immobilization in a porcine model for critical illness myopathy (CIM) [61, 62].
Co-expressed genes in module dark-turquoise were significantly associated with meat quality and based on the current knowledge of gene functions some links among them are suggested to be relevant For example, one member of this module, CREM, is a transcription factor binding to cAMP-responsive elements (CREs) in the promoters of various genes. This transcription factor plays important roles in various organismal functions [63–67]. Crem inactivation or knockout has been shown to increase the rate of apoptosis in testis tissue [68, 69]. The main cellular change associated with apoptosis processes also occur during post-mortem . Post mortem biochemical processes in muscle lead to pH decline. A high expression of CREM being positively correlated with pH at 24 hours, may indicate a slowdown of apoptosis related post mortem processes paralleling anaerobic metabolic processes that led to a decrease of pH. This indicates that the abundance of CREM transcripts in muscle plays a significant role in meat quality. Further, CREM and miR-153 were highly negatively correlated (FDR < 0.01) which is known to induce apoptosis in a glioblastoma cell line DBTRG-05MG . Thus a consistent link of effects on apoptosis and mRNA-miRNA interaction is obvious. miR-153 also inhibits the protein kinase B (PKB/Akt) pathway by reducing the protein level of insulin receptor substrate-2 (IRS2) . As recently shown, miR-135a targets IRS2 levels by binding to its 3′UTR and this interaction regulates skeletal muscle insulin signaling . Insulin signaling plays a pivotal role in the regulation of glucose uptake by skeletal muscle . The glucose uptake in skeletal muscle has large effects on meat characteristics [75, 76]. In our study, Irs-2 also belonged to module dark-turquoise and was negatively correlated with waterholding capacity related traits like DL (p = 0.004) and positively with pH (pH24MLD, p = 7.05E-07; pH24MSM, p = 1.30E-03; pH45MLD, p = 3.39E-03). Thus another plausible functional link of members of the module dark-turquoise and miRNAs can be shown.
Mitochondria supply energy for physiological functions and play a significant role in the regulation of other cellular events including apoptosis, calcium homeostasis, and production of reactive oxygen species. Mitochondrial metabolism is affected by miRNA regulation . Here, we found many miRNAs being negatively correlated to target genes of modules red, black, which were enriched for genes related to mitochondrial pathways. Indeed, miR-338 modulates energy metabolism, oxidative phosphorylation, and mitochondrial functions [78, 79], and miR-15b, -16, -195, and -424 decrease cellular ATP levels in cardiomyocytes . Additionally, miR-181c can enter and target the mitochondrial genome, ultimately causing electron transport chain complex IV remodeling and mitochondrial dysfunction . Here we found miR-181 was correlated with DL and pH24MLD as well as S100 calcium binding protein A6 (S100A6). Further, mitochondrial genes like Kif1b, Atp6v1g1, Atp5c1, Park7, Chchd2, Ruvbl1, Mrps7, and Mrpl48 were highly negatively correlated with, and some of them were predicted as targets of, miR-15, -154, -17, -503, -214, -330, -342, and -375.
Many genes of modul tan were assigned in the GO category “extracellular matrix” including Col1a1, Col1a2, Col3a1, Col5a1, Col5a2, Col6a1, Col6a3, Col12a1, Col14a1, Crispld2, Ddn, Fbn1, Nid1, Nid2, Ogn, and Postn. These genes were negatively correlated with, and some of them were predicted as targets of, miRNA. Interestingly, Col6a1 was found as a target for many miRNAs including miR-132, miR-205, miR-320, miR-330, miR-375, miR-423, and miR-425. Moreover, miR-205 was identified as a regulator of Col1a1, Maged1, and Dpysl3. Gandellini et al. (2012)  reported that miR-205 controls the deposition of laminin-332 and its receptor integrin-β4 as well as participates in a network involving ΔNp63α, which is essential for maintenance of the basement membrane in prostate epithelium. Similarly, other miRNAs, including miR-29, miR-133, and miR-30, are involved in the regulation of development and maintenance of extracellular matrix of bone and muscle [83, 84].
Much evidence suggests that a group of miRNAs (cluster and family) may contribute to the regulation of a set of common targets [46, 85, 86], and are, therefore, associated with phenotypes. WGCNA was used here to group miRNA products and revealed 2 modules associated with meat quality (purple and blue). The miRNA from these modules were negatively correlated to mRNAs, and some of these were predicted as targets. Most miRNAs in module purple were related to genes in the categories “ubiquitin” or “protein catabolic process” (Ube2e1, Ube2b1, Ubap1, Igfbp5, Herc1, Psmb3, Flj25076). Differential expression of genes of the ubiquitin system depending on muscle and meat quality was previously shown; only recently also association of genes of the ubiquitin system with meat quality was reported [1, 87]. In particular, miR-324 was highly negatively correlated and predicted to target Ube2e1, Ube2v1, Bap1, Lpf and Cic.
Most previous muscle mRNA and miRNA expression studies focused on cardiac muscle or skeletal muscle injury [50, 88]. In these injuries, the degree of damage results from an imbalance of energy, nutrients, and oxygen supply-and-demand in muscle cells. Similarly, nutrient, energy, and oxygen depletion occurs post mortem. Many changes in expression associated with muscle injury would therefore overlap with post mortem processes in conversion of muscle to meat, and vice versa. In this regard, functional annotation of mRNA co-expression and trait-correlated expression identified key post mortem pathways and functions, including glucose metabolic process, mitochondrial metabolic pathways, and muscle structural components, involved in muscle-to-meat conversion that will be relevant to muscle injury as well.
In this study, for the first time, expression and co-expression of miRNAs—functioning as a fine-tuning of mRNA transcription and translation—was integrated with mRNA transcript abundance measures and phenotypic data on meat quality. By this an additional hierarchical level, i.e. miRNA affecting mRNA, was considered in the molecular regulation of muscle-to-meat conversion. miRNAs are necessary for proper skeletal and cardiac muscle development and function, and have a profound influence on multiple myopathies, such as hypertrophy, dystrophy, and conduction defects. Consequently, an expression biomarker panel (whether from mRNA or miRNA) derived from this study may not only be predictive for quality of meat post mortem, but also for many muscle pathologic processes including muscle atrophy, dystrophy, and hypoxia [89, 90]. The abundance of mRNAs and their fine-tuning by corresponding miRNAs in molecular pathways related to mitochondrial metabolic balance and oxidative stress, cell proliferation and differentiation, as well as muscle structural protein composition play an important role in these myopathies and meat maturation. Functional studies of the interactions among and between mRNAs and μRNAs will provide additional experimental data for validation of the relationships on the level of mRNAs, miRNAs and organismal phenotype that were stressed in this study.
Animals, tissue collection, and phenotyping
Animal care and tissue collection procedures followed the guidelines of the German Law of Animal Protection, and the experimental protocol was approved by the Animal Care Committee of the FBN. This study was based on trait measurement and expression profile association analyses done with 207 performance-tested pigs from commercial herds of the crossbreed Pietrain × (German Large White × German Landrace). Animals were raised and slaughtered under standardized conditions in the experimental facilities of the Leibniz Institute for Farm Animal Biology (FBN). Sample collection was performed thoroughly after exsanguination, tissue samples were rapidly dissected, snap-frozen in liquid nitrogen and stored at -80°C. The average age of the pigs at sampling was ~180 days. Technological parameters of meat quality, i.e., pH-value, conductivity, and colour, were measured by using Star-series equipment (Rudolf Matthaeus Company, Germany). Measures of pH and conductivity were at 45 min post mortem (pH45) and 24 h post mortem (pH24), in both M. longissimus dorsi between 13th/14th rib (pH45MLD, pH24MLD, LF45MLD, LF24MLD) and the ham (M. semimembranosus) (symbol: pH24MSM, LF24MSM). Muscle colour was measured at 24 h post mortem by Opto-Star (Matthaeus, Klausa, Germany). Drip loss was scored based on a bag method with a size-standardized sample from the M. longissimus dorsi collected at 24 h post mortem and weighed, suspended in a plastic bag, held at 4°C for 48 h, and re-weighed . To determine cooking loss, a loin cube was taken from the M. longissimus dorsi, weighed, placed in a polyethylene bag, and incubated in water at 75°C for 50 minutes. The bag was then immersed in flowing water at room temperature for 30 minutes, and the solid portion was re-weighed. Thawing loss was determined similarly after at least 24 h freezing at -20°C. Drip loss, cooking loss, and thawing loss were calculated as a percentage of weight loss based on the start weight of a sample. Shear force was measured using Instron-4310 equipment, and average values of four replicates were used for analyses.
Customized miRNA microarrays design
Our custom porcine miRNA array was designed from 284 pig miRNAs obtained from the miRBASE (miRBase 14.0). Because miRNAs are highly conserved between closely-related species , we could predict novel porcine miRNA candidates by inter-species alignments requiring 100% mature miRNA similarity. Accordingly, we used previously known miRNA sequences from humans and mice, as well as other species, to perform BLAST searches against the porcine genome database porcine; 391 miRNA candidates were identified. In total, 675 miRNAs probe sets, corresponding to 513 unique sequences belonging to 159 miRNA families, were used for hybridisation with the target samples described above. Microarray data related to all samples were deposited in the Gene Expression Omnibus public repository (GEO accession number: GSE41294: GSM1013731-GSM1013920).
Customized microarrays, pre-processing, and normalization of miRNA
Total miRNA was isolated with Qiagen miReasy Mini kit and RNeasy MinElute Clean up kit (Qiagen, Hilden, Germany) according to manufacturer’s protocol for small RNA. Quality and quantity of isolated total RNA and miRNA were determined using an Agilent 2100 Bioanalyzer for RNA (Agilent Technologies, Santa Clara, CA). Affymetrix customized microarrays from our porcine miRNA candidate dataset were used. Targets for hybridisation were prepared from miRNA with the FlashTag™ Biotin RNA Labeling Kit for Affymetrix GeneChip miRNA arrays (Genisphere, Hatfield, PA, USA) according to manufacturer’s recommendations. Briefly, 250 ng of miRNA of each individual were poly(A)-tailed using ATP–poly-A-Polymerase, then FlashTag Biotin end-labelled. After hybridisation of biotin-labelled complementary RNA, chips were washed and processed to detect biotin-containing transcripts by Streptavidin-PE (Phycoerythrin) conjugate, then were scanned on GeneChip scanner 3000 7G (Affymetrix, Santa Clara, US). Data were extracted from the images, and spots were quantified and processed by quality filtering. Hybridisation quality was assessed in all samples by using JMP Genomics 5 utilising Robust Multi-array Average (RMA) background correction and log2 transformations. To acquire the expression value, data were normalized between chips using the quantile normalization method.
Whole-genome expression profiling (mRNA)
Gene expression profiling of M. longissimus dorsi samples of pigs was conducted with the same animals (207) as for miRNA. In brief, total RNA was isolated using TRI Reagent (Sigma, Taufkirchen, Germany) and used for target preparation for microarray hybridisation. According to Affymetrix protocols, 500 ng of total RNA were reverse-transcribed into cDNA, transcribed into cRNA, and labelled using Affymetrix One cycle synthesis and labelling kit (Affymetrix, UK) to prepare antisense biotinylated RNA targets. Quality of hybridisation was assessed in all samples following manufacturer’s recommendations. Data were analysed with the Affymetrix GCOS 1.1.1 software, using global scaling to a target signal of 500. Data were processed with MAS5.0 to generate cell intensity files (present or absent). Quantitative expression levels of present transcripts were estimated using the PLIER algorithm (Probe Logarithmic Intensity Error) for normalization that was implemented in Expression Console (Affymetrix). Based on BLAST comparison of the Affymetrix porcine target sequences with the porcine genome sequence (Ensembl_Sscrofa_10), 20,689 of the 24,123 probe sets on the Affymetrix Porcine GeneChip were localized and annotated [92, 93]. Microarray data related to all samples were deposited in the Gene Expression Omnibus public repository [GEO accession number: GSE32112: GSM796045-GSM796251].
Pre-processing of phenotype and expression data
Phenotypes and expression levels were adjusted for systematic effects by analysis of variance performed with the procedure “Mixed” of the SAS software package (SAS version 9.1 SAS Institute, Cary, NC) before analysing their correlation and by using co-expression network. Sex and RYR genotype was used as a fixed effect, “sire” and “slaughter day” as random effects, and “carcass weight” as a covariate. Subsequently, the residuals of log2-transformed expression intensities (miRNA and mRNA) and muscle phenotype were used for further analysis.
Pearson correlation of miRNA expression level and gene expression level was calculated using 190 individuals; correction for multiple testing was done by controlling the FDR level (q-value according to Storey and Tibshirani, 2003, ) at 10%.
Weighted gene Co-expression network analysis (WGCNA)
A weighted gene co-expression network was constructed for 207 muscle biopsies using the blockwise Modules function from the WGCNA package in R . Residuals of gene expression, after correcting the effect, were used for WGCNA. The blockwise Modules function allows the entire dataset of 11,191 probe sets by mRNA and 675 miRNA to be utilised in the construction of the weighted gene co-expression network.
Extremely outlying individuals were removed from the following analysis based on hierarchically clustered using the average linkage function, and common Euclidean distance. Pearson correlation matrix of all gene-gene comparisons were calculated across all microarrays. Adjacency matrix was then calculated using the correlation matrix of the expression sets. Finally, the Topological Overlap Matrix (TOM) was converted from the adjacency matrix and used to derive a TOM-based distance matrix for the hierarchical clustering of expressions. In the next step, modules of expression profiles, (i.e. sets of genes with high topological overlap) were formed based on hierarchical clustering, with empirically specified minimal module size (30 for gene expression, 10 for miRNA).
According to the WGCNA methodology, rather than traditional distance or correlation based similarity measures, it utilizes the topological overlap matrix Ω = [ω ij ] (TOM),
where x i and x j are the gene expression profiles of the i-th and j-th gene and α ij is the adjacency. TOM based distance matrix is a robust and powerful measurement in building co-expression network. Selection of appropriate value for the power β were derived according to the pickSoftThresholding function of the WGCNA package [23, 95]. Accordingly, by manually inspecting the fit of the scale free topology model with the candidate β-values for each set of expression profile, minimal β-values giving a coefficient of determination R2 higher than 90% were adopted.
Modules were further merged based on the dissimilarity between their “eigengenes”, which were defined as the first principal component of each module. Genes that were not assigned to another module were assigned to the grey module. A threshold of 0.2 for the dissimilarity as recommended by the WGCNA author was used. Module–trait associations were estimated using the correlation between the module eigengene and the phenotype, which allows easy identification of expression set (module) highly correlated to the phenotype. For each expression profile, Gene Significance (GS) was calculated as the absolute value of the correlation between expression profile and each trait; module membership (MM) was defined as the correlation of expression profile and each module eigengene, enabling further identification of key players in the regulation network.
Gene ontology and pathway enrichment analysis
We performed a gene ontology (GO) enrichment analysis for network modules using the Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/[96, 97]. Each analysis was performed using the functional annotation clustering option. Functional annotation clustering combines single categories with a significant overlap in gene content and assigns an enrichment score (ES, defined as the -log10 of the geometric mean of unadjusted p-value s for each single term in the cluster) to each cluster, making interpretation of the results more straightforward. To assess the significance of functional clusters, we created 22 sets of 11,191 probe sets corresponding to 8,036 genes (size of the average module identified in this study).
Predicting porcine targets of miRNAs by RNA hybrid and TargetScan
We used two methods to predict the targets of porcine miRNA. First, we predicted targets using the computational software RNAhybrid (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid), which detects the most energetically favourable hybridisation sites of a small RNA within a large RNA . Here, we tested the miRNA probe sets with the following parameters: number of hits per target = 1, energy cutoff = -25 kcal/mol, and maximal internal or bulge loop size per side = 4. Most targets found were located on the 3′-UTR of genes. Second, TargetScan (http://www.targetscan.org) was used to detect target gene candidates based on seed complementarity on UTR database 6.0 and our porcine RefSeq transcript with our miRNA seed sequence . TargetScan was applied considering both conserved and non-conserved targets. The porcine RefSeq transcripts, which derived from 11,191 probe sets that showed consistent expression in porcine muscle, were used as input targets for RNAhybrid and TargetScan.
Ponsuksili S, Murani E, Phatsara C, Schwerin M, Schellander K, Wimmers K: Porcine muscle sensory attributes associate with major changes in gene networks involving CAPZB, ANKRD1, and CTBP2. Funct Integr Genomics. 2009, 9: 455-471. 10.1007/s10142-009-0131-1.
Ponsuksili S, Murani E, Phatsara C, Schwerin M, Schellander K, Wimmers K: Expression quantitative trait loci analysis of genes in porcine muscle by quantitative real time RT-PCR compared to microarray data. Heredity. 2010, 105 (3): 309-317. 10.1038/hdy.2010.5.
Sellier P: The future role of molecular genetics in the control of meat production and meat quality. Meat Sci. 1994, 36 (1–2): 29-44.
Bode G, Clausing P, Gervais F, Loegsted J, Luft J, Nogues V, Sims J: The utility of the minipig as an animal model in regulatory toxicology. J Pharmacol Toxicol Methods. 2010, 62 (3): 196-220. 10.1016/j.vascn.2010.05.009.
Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, Rogel-Gaillard C, Park C, Milan D, Megens HJ: Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012, 491 (7424): 393-398. 10.1038/nature11622.
Ponsuksili S, Jonas E, Murani E, Phatsara C, Srikanchai T, Walz C, Schwerin M, Schellander K, Wimmers K: Trait correlated expression combined with expression QTL analysis reveals biological pathways and candidate genes affecting water holding capacity of muscle. BMC Genomics. 2008, 9: 367-10.1186/1471-2164-9-367.
Ponsuksili S, Murani E, Phatsara C, Jonas E, Walz C, Schwerin M, Schellander K, Wimmers K: Expression Profiling of Muscle Reveals Transcripts Differentially Expressed in Muscle That Affect Water-Holding Capacity of Pork. J Agric Food Chem. 2008, 56: 10311-10317. 10.1021/jf800881y.
Ponsuksili S, Murani E, Brand B, Schwerin M, Wimmers K: Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model. J Lipid Res. 2011, 52 (4): 668-678. 10.1194/jlr.M013342.
Callis TE, Deng Z, Chen JF, Wang DZ: Muscling through the microRNA world. Exp Biol Med. 2008, 233 (2): 131-138. 10.3181/0709-MR-237.
Huang TH, Zhu MJ, Li XY, Zhao SH: Discovery of Porcine microRNAs and Profiling from Skeletal Muscle Tissues during Development. PLoS One. 2008, 3 (9): e3225-10.1371/journal.pone.0003225.
van Rooij E, Liu N, Olson EN: MicroRNAs flex their muscles. Trends in Genet. 2008, 24 (4): 159-166. 10.1016/j.tig.2008.01.007.
Walden TB, Timmons JA, Keller P, Nedergaard J, Cannon B: Distinct expression of muscle-specific microRNAs (myomirs) in brown adipocytes. J Cell Physiol. 2009, 218 (2): 444-449. 10.1002/jcp.21621.
Xie H, Lim B, Lodish HF: MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes. 2009, 58 (5): 1050-1057. 10.2337/db08-1299.
Ortega FJ, Moreno-Navarrete JM, Pardo G, Sabater M, Hummel M, Ferrer A, Rodriguez-Hermosa JI, Ruiz B, Ricart W, Peral B, Fernández-Real JM: MiRNA expression profile of human subcutaneous adipose and during adipocyte differentiation. PLoS One. 2010, 5 (2): e9022-10.1371/journal.pone.0009022.
Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol. 2009, 10 (2): 126-139. 10.1038/nrm2632.
Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11 (3): 228-234. 10.1038/ncb0309-228.
Fang Z, Rajewsky N: The impact of miRNA target sites in coding sequences and in 3′UTRs. PLoS ONE. 2011, 6 (3): e18067-10.1371/journal.pone.0018067.
Guo H, Ingolia NT, Weissman JS, Bartel DP: Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. 2010, 466 (7308): 835-840. 10.1038/nature09267.
Peng X, Li Y, Walters KA, Rosenzweig ER, Lederer SL, Aicher LD, Proll S, Katze MG: Computational identification of hepatitis C virus associated microRNA-mRNA regulatory modules in human livers. BMC Genomics. 2009, 10: 373-10.1186/1471-2164-10-373.
Liu H, Brannon AR, Reddy AR, Alexe G, Seiler MW, Arreola A, Oza JH, Yao M, Juan D, Liou LS, Ganesan S, Levine AJ, Rathmell WK, Bhanot GV: Identifying mRNA targets of microRNA dysregulated in cancer: with application to clear cell Renal Cell Carcinoma. BMC Syst Biol. 2010, 4: 51-10.1186/1752-0509-4-51.
Chen Y, Zhu J, Lum PY, Yang X, Pinto S, MacNeil DJ, Zhang C, Lamb J, Edwards S, Sieberts SK: Variations in DNA elucidate molecular networks that cause disease. Nature. 2008, 452: 429-435. 10.1038/nature06757.
Gargalovic PS, Imura M, Zhang B, Gharavi NM, Clark MJ, Pagnon J, Yang WP, He A, Truong A, Patel S: Identification of inflammatory gene modules based on variations of human endothelial cell responses to oxidized lipids. Proc Natl Acad Sci USA. 2006, 103: 12741-12746. 10.1073/pnas.0605457103.
Langfelder P, Horvath S: WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 2008, 9 (1): 559-10.1186/1471-2105-9-559.
Saris CG, Horvath S, van Vught PW, van Es MA, Blauw HM, Fuller TF, Langfelder P, DeYoung J, Wokke JH, Veldink JH, van den Berg LH, Ophoff RA: Weighted gene co-expression network analysis of the peripheral blood from Amyotrophic Lateral Sclerosis patients. BMC Genomics. 2009, 10: 405-10.1186/1471-2164-10-405.
Plaisier CL, Horvath S, Huertas-Vazquez A, Cruz-Bautista I, Herrera MF, Tusie-Luna T, Aguilar-Salinas C, Pajukanta P: A systems genetics approach implicates USF1, FADS3, and other causal candidate genes for familial combined hyperlipidemia. PLoS Genet. 2009, 5 (9): e1000642-10.1371/journal.pgen.1000642.
Farber CR: Identification of a gene module associated with BMD through the integration of network analysis and genome-wide association data. J Bone Miner Res. 2010, 25 (11): 2359-2367. 10.1002/jbmr.138.
DiLeo MV, Strahan GD, den Bakker M, Hoekenga OA: Weighted correlation network analysis (WGCNA) applied to the tomato fruit metabolome. PLoS One. 2011, 6 (10): e26683-10.1371/journal.pone.0026683.
Lancu OD, Darakjian P, Malmanger B, Walter NA, McWeeney S, Hitzemann R: Gene networks and haloperidol-induced catalepsy. Genes Brain Behav. 2012, 11 (1): 29-37. 10.1111/j.1601-183X.2011.00736.x.
Rosenvold K, Andersen HJ: Factors of significance, for pork quality - a review. Meat science. 2003, 64: 219-237. 10.1016/S0309-1740(02)00186-9.
Figueiredo PA, Mota MP, Appell HJ, Duarte JA: The role of mitochondria in aging of skeletal muscle. Biogerontology. 2008, 9 (2): 67-84. 10.1007/s10522-007-9121-7.
Romanello V, Guadagnin E, Gomes L, Roder I, Sandri C, Petersen Y, Milan G, Masiero E, Del Piccolo P, Foretz M, Scorrano L, Rudolf R, Sandri M: Mitochondrial fission and remodelling contributes to muscle atrophy. EMBO J. 2010, 29 (10): 1774-1785. 10.1038/emboj.2010.60.
Dorn GW: Mitochondrial dynamics in heart disease. Biochim Biophys Acta. 2013, 1833 (1): 233-241. 10.1016/j.bbamcr.2012.03.008.
Irwin MH, Parameshwaran K, Pinkert CA: Mouse models of mitochondrial complex I dysfunction. Int J Biochem Cell Biol. 2013, 45 (1): 34-40. 10.1016/j.biocel.2012.08.009.
Bee G, Anderson AL, Lonergan SM, Huff-Lonergan E: Rate and extent of pH decline affect proteolysis of cytoskeletal proteins and water-holding capacity in pork. Meat Sci. 2007, 76 (2): 359-365. 10.1016/j.meatsci.2006.12.004.
Huff-Lonergan E, Lonergan SM: Mechanisms of water-holding capacity of meat: The role of postmortem biochemical and structural changes. Meat Sci. 2005, 71 (1): 194-204. 10.1016/j.meatsci.2005.04.022.
Camacho Vanegas O, Bertini E, Zhang RZ, Petrini S, Minosse C, Sabatelli P, Giusti B, Chu ML, Pepe G: Ullrich scleroatonic muscular dystrophy is caused by recessive mutations in collagen type VI. Proc Natl Acad Sci USA. 2001, 98 (13): 7516-7521. 10.1073/pnas.121027598.
Lampe AK, Bushby KM: Collagen VI related muscle disorders. J Med Genet. 2005, 42 (9): 673-685. 10.1136/jmg.2002.002311.
Merlini L, Martoni E, Grumati P, Sabatelli P, Squarzoni S, Urciuolo A, Ferlini A, Gualandi F, Bonaldo P: Autosomal recessive myosclerosis myopathy is a collagen VI disorder. Neurology. 2008, 71 (16): 1245-1253. 10.1212/01.wnl.0000327611.01687.5e.
Gualandi F, Urciuolo A, Martoni E, Sabatelli P, Squarzoni S, Bovolenta M, Messina S, Mercuri E, Franchella A, Ferlini A, Bonaldo P, Merlini L: Autosomal recessive Bethlem myopathy. Neurology. 2009, 73 (22): 1883-1891. 10.1212/WNL.0b013e3181c3fd2a.
Allamand V, Briñas L, Richard P, Stojkovic T, Quijano-Roy S, Bonne G: ColVI myopathies: where do we stand, where do we go?. Skelet Muscle. 2011, 23 (1): 30-
Irwin WA, Bergamin N, Sabatelli P, Reggiani C, Megighian A, Merlini L, Braghetta P, Columbaro M, Volpin D, Bressan GM, Bernardi P, Bonaldo P: Mitochondrial dysfunction and apoptosis in myopathic mice with collagen VI deficiency. Nat Genet. 2003, 35 (4): 367-371. 10.1038/ng1270.
Grumati P, Coletto L, Sabatelli P, Cescon M, Angelin A, Bertaggia E, Blaauw B, Urciuolo A, Tiepolo T, Merlini L, Maraldi NM, Bernardi P, Sandri M, Bonaldo P: Autophagy is defective in collagen VI muscular dystrophies, and its reactivation rescues myofiber degeneration. Nat Med. 2010, 16 (11): 1313-1320. 10.1038/nm.2247.
Friedman RC, Farh KKH, Burge CB, Bartel DP: Most mammalian mrnas are conserved targets of micrornas. Genome Res. 2009, 19 (1): 92-105.
Altuvia Y, Landgraf P, Lithwick G, Elefant N, Pfeffer S, Aravin A, Brownstein MJ, Tuschl T, Margalit H: Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33 (8): 2697-2706. 10.1093/nar/gki567.
Tanzer A, Stadler PF: Molecular evolution of a microrna cluster. J Mol Biol. 2004, 339 (2): 327-335. 10.1016/j.jmb.2004.03.065.
Williams AH, Liu N, van Rooij E, Olson EN: MicroRNA control of muscle development and disease. Curr Opin Cell Biol. 2009, 21 (3): 461-469. 10.1016/j.ceb.2009.01.029.
Eisenberg I, Eran A, Nishino I, Moggio M, Lamperti C, Amato AA, Lidov HG, Kang PB, North KN, Mitrani-Rosenbaum S, Flanigan KM, Neely LA, Whitney D, Beggs AH, Kohane IS, Kunkel LM: Distinctive patterns of microRNA expression in primary muscular disorders. Proc Natl Acad Sci USA. 2007, 104 (43): 17016-17021. 10.1073/pnas.0708115104.
Eisenberg I, Alexander MS, Kunkel LM: miRNAS in normal and diseased skeletal muscle. J Cell Mol Med. 2009, 13 (1): 2-11.
Cardinali B, Castellani L, Fasanaro P, Basso A, Alema S, Martelli F, Falcone G: Microrna-221 and microrna-222 modulate differentiation and maturation of skeletal muscle cells. PLoS One. 2009, 4 (10): e7607-10.1371/journal.pone.0007607.
Greco S, Perfetti A, Fasanaro P, Cardani R, Capogrossi MC, Meola G, Martelli F: Deregulated MicroRNAs in Myotonic Dystrophy Type 2. PLoS One. 2012, 7 (6): e39732-10.1371/journal.pone.0039732.
Xie SS, Huang TH, Shen Y, Li XY, Zhang XX, Zhu MJ, Qin HY, Zhao SH: Identification and characterization of microRNAs from porcine skeletal muscle. Anim Genet. 2010, 41 (2): 179-190. 10.1111/j.1365-2052.2009.01991.x.
Liu N, Bezprozvannaya S, Shelton JM, Frisard MI, Hulver MW, McMillan RP, Wu Y, Voelker KA, Grange RW, Richardson JA, Bassel-Duby R, Olson EN: Mice lacking microRNA 133a develop dynamin 2–dependent centronuclear myopathy. J Clin Invest. 2011, 121 (8): 3258-3268. 10.1172/JCI46267.
Wang L, Chen X, Zheng Y, Li F, Lu Z, Chen C, Liu J, Wang Y, Peng Y, Shen Z, Gao J, Zhu M, Chen H: MiR-23a inhibits myogenic differentiation through down regulation of fast myosin heavy chain isoforms. Exp Cell Res. 2012, 318 (18): 2324-2334. 10.1016/j.yexcr.2012.06.018.
Li G, Wu Z, Li X, Ning X, Li Y, Yang G: Biological role of microRNA-103 based on expression profile and target genes analysis in pigs. Mol Biol Rep. 2011, 38 (7): 4777-4786. 10.1007/s11033-010-0615-z.
McGregor RA, Choi MS: microRNAs in the regulation of adipogenesis and obesity. Curr Mol Med. 2011, 11 (4): 304-316. 10.2174/156652411795677990.
Trajkovski M, Hausser J, Soutschek J, Bhat B, Akin A, Zavolan M, Heim MH, Stoffel M: MicroRNAs 103 and 107 regulate insulin sensitivity. Nature. 2011, 474 (7353): 649-653. 10.1038/nature10112.
Grimson A, Farh KKH, Johnston KW, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA Targeting Specificity in Mammals: Determinants beyond Seed Pairing. Mol Cell. 2007, 27 (1): 91-105. 10.1016/j.molcel.2007.06.017.
Mazière P, Enright AJ: Prediction of microRNA targets. Drug Discov Today. 2007, 12 (11–12): 452-458.
Rajan S, Chu Pham Dang H, Djambazian H, Zuzan H, Fedyshyn Y, Ketela T, Moffat J, Hudson TJ, Sladek R: Analysis of early C2C12 myogenesis identifies stably and differentially expressed transcriptional regulators whose knock-down inhibits myoblast differentiation. Physiol Genomics. 2012, 44 (2): 183-197. 10.1152/physiolgenomics.00093.2011.
Reiner G, Heinricy L, Müller E, Geldermann H, Dzapo V: Indications of associations of the porcine FOS proto-oncogene with skeletal muscle fibre traits. Anim Genet. 2002, 33 (1): 49-55. 10.1046/j.1365-2052.2002.00805.x.
Banduseela VC, Ochala J, Chen YW, Göransson H, Norman H, Radell P, Eriksson LI, Hoffman EP, Larsson L: Gene expression and muscle fiber function in a porcine ICU model. Physiol Genomics. 2009, 39 (3): 141-159. 10.1152/physiolgenomics.00026.2009.
Banduseela VC, Chen YW, Göransson Kultima H, Norman HS, Aare S, Radell P, Eriksson LI, Hoffman EP, Larsson L: Impaired autophagy, chaperone expression and protein synthesis in response to critical illness interventions in porcine skeletal muscle. Physiol Genomics. 2013, 45 (12): 477-486. 10.1152/physiolgenomics.00141.2012.
Don J, Stelzer G: The expanding family of CREB/CREM transcription factors that are involved with spermatogenesis. Mol Cell Endocrinol. 2002, 187 (1–2): 115-124.
Müller FU, Lewin G, Baba HA, Bokník P, Fabritz L, Kirchhefer U, Kirchhof P, Loser K, Matus M, Neumann J, Riemann B, Schmitz W: Heart-directed expression of a human cardiac isoform of cAMP-response element modulator in transgenic mice. J Biol Chem. 2005, 280 (8): 6906-6914. 10.1074/jbc.M407864200.
Wu X, Jin W, Liu X, Fu H, Gong P, Xu J, Cui G, Ni Y, Ke K, Gao Z, Gao Y: Cyclic AMP response element modulator-1 (CREM-1) involves in neuronal apoptosis after traumatic brain injury. J Mol Neurosci. 2012, 47 (2): 357-367. 10.1007/s12031-012-9761-1.
Juang YT, Wang Y, Solomou EE, Li Y, Mawrin C, Tenbrock K, Kyttaris VC, Tsokos GC: Systemic lupus erythematosus serum IgG increases CREM binding to the IL-2 promoter and suppresses IL-2 production through CaMKIV. J Clin Invest. 2005, 115 (4): 996-1005.
Xu WD, Zhang YJ, Wang W, Li R, Pan HF, Ye DQ: Role of CREM in systemic lupus erythematosus. Cell Immunol. 2012, 276 (1–2): 10-15.
Nantel F, Monaco L, Foulkes NS, Masquilier D, LeMeur M, Henriksén K, Dierich A, Parvinen M, Sassone-Corsi P: Spermiogenesis deficiency and germ-cell apoptosis in CREM-mutant mice. Nature. 1996, 380 (6570): 159-162. 10.1038/380159a0.
Kosir R, Juvan P, Perse M, Budefeld T, Majdic G, Fink M, Sassone-Corsi P, Rozman D: Novel insights into the downstream pathways and targets controlled by transcription factors CREM in the testis. PLoS One. 2012, 7 (2): e31798-10.1371/journal.pone.0031798.
Lonergan EH, Zhang W, Lonergan SM: Biochemistry of postmortem muscle - Lessons on mechanisms of meat tenderization. Meat Sci. 2010, 86 (1): 184-195. 10.1016/j.meatsci.2010.05.004.
Xu J, Liao X, Wong C: Downregulations of B-cell lymphoma 2 and myeloid cell leukemia sequence 1 by microRNA 153 induce apoptosis in a glioblastoma cell line DBTRG-05MG. Int J Cancer. 2010, 126 (4): 1029-1035.
Xu J, Liao X, Lu N, Liu W, Wong CW: Chromatin-modifying drugs induce miRNA-153 expression to suppress Irs-2 in glioblastoma cell lines. Int J Cancer. 2011, 129 (10): 2527-2531. 10.1002/ijc.25917.
Agarwal P, Srivastava R, Srivastava AK, Ali S, Datta M: miR-135a targets IRS2 and regulates insulin signaling and glucose uptake in the diabetic gastrocnemius skeletal muscle. Biochim Biophys Acta. 2013, 1832 (8): 1294-1303. 10.1016/j.bbadis.2013.03.021.
Kubota T, Kubota N, Kumagai H, Yamaguchi S, Kozono H, Takahashi T, Inoue M, Itoh S, Takamoto I, Sasako T, Kumagai K, Kawai T, Hashimoto S, Kobayashi T, Sato M, Tokuyama K, Nishimura S, Tsunoda M, Ide T, Murakami K, Yamazaki T, Ezaki O, Kawamura K, Masuda H, Moroi M, Sugi K, Oike Y, Shimokawa H, Yanagihara N, Tsutsui M, Terauchi Y, Tobe K, Nagai R, Kamata K, Inoue K, Kodama T, Ueki K, Kadowaki T: Impaired insulin signaling in endothelial cells reduces insulin-induced glucose uptake by skeletal muscle. Cell Metab. 2011, 13 (3): 294-307. 10.1016/j.cmet.2011.01.018.
Milan D, Jeon JT, Looft C, Amarger V, Robic A, Thelander M, Rogel-Gaillard C, Paul S, Iannuccelli N, Rask L, Ronne H, Lundström K, Reinsch N, Gellin J, Kalm E, Roy PL, Chardon P, Andersson L: A mutation in PRKAG3 associated with excess glycogen content in pig skeletal muscle. Science. 2000, 288 (5469): 1248-1251. 10.1126/science.288.5469.1248.
Andersson L: Identification and characterization of AMPK gamma 3 mutations in the pig. Biochem Soc Trans. 2003, 31 (Pt 1): 232-235.
Li P, Jiao J, Gao G: Prabhakar BS. Control of mitochondrial activity by miRNAs. J Cell Biochem. 2012, 113 (4): 1104-1110. 10.1002/jcb.24004.
Aschrafi A, Schwechter AD, Mameza MG, Natera-Naranjo O, Gioio AE, Kaplan BB: MicroRNA-338 regulates local cytochrome c oxidase IV mRNA levels and oxidative phosphorylation in the axons of sympathetic neurons. J Neurosci. 2008, 28 (47): 12581-12590. 10.1523/JNEUROSCI.3338-08.2008.
Aschrafi A, Kar AN, Natera-Naranjo O, Macgibeny MA, Gioio AE, Kaplan BB: MicroRNA-338 regulates the axonal expression of multiple nuclear-encoded mitochondrial mRNAs encoding subunits of the oxidative phosphorylation machinery. Cell Mol Life Sci. 2012, Jul 8. [Epub ahead of print]: PMID:22773120
Nishi H, Ono K, Iwanaga Y, Horie T, Nagao K, Takemura G, Kinoshita M, Kuwabara Y, Mori RT, Hasegawa K, Kita T, Kimura T: MicroRNA-15b modulates cellular ATP levels and degenerates mitochondria via Arl2 in neonatal rat cardiac myocytes. J Biol Chem. 2010, 285 (7): 4920-4930. 10.1074/jbc.M109.082610.
Das S, Ferlito M, Kent OA, Fox-Talbot K, Wang R, Liu D, Raghavachari N, Yang Y, Wheelan SJ, Murphy E, Steenbergen C: Nuclear miRNA regulates the mitochondrial genome in the heart. Circ Res. 2012, 110 (12): 1596-1603. 10.1161/CIRCRESAHA.112.267732.
Gandellini P, Profumo V, Casamichele A, Fenderico N, Borrelli S, Petrovich G, Santilli G, Callari M, Colecchia M, Pozzi S, De Cesare M, Folini M, Valdagni R, Mantovani R, Zaffaroni N: miR-205 regulates basement membrane deposition in human prostate: implications for cancer development. Cell Death Differ. 2012, 19 (11): 1750-1760. 10.1038/cdd.2012.56. PMID: 22555458
Li Z, Hassan MQ, Jafferji M, Aqeilan RI, Garzon R, Croce CM, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Biological functions of miR-29b contribute to positive regulation of osteoblast differentiation. J Biol Chem. 2009, 284 (23): 15676-15684. 10.1074/jbc.M809787200.
Duisters RF, Tijsen AJ, Schroen B, Leenders JJ, Lentink V, van der Made I, Herias V, van Leeuwen RE, Schellings MW, Barenbrug P, Maessen JG, Heymans S, Pinto YM, Creemers EE: miR-133 and miR-30 regulate connective tissue growth factor: implications for a role of microRNAs in myocardial matrix remodeling. Circ Res. 2009, 104 (2): 170-178. 10.1161/CIRCRESAHA.108.182535.
Carraro G, El-Hashash A, Guidolin D, Tiozzo C, Turcatel G, Young BM, De Langhe SP, Bellusci S, Shi W, Parnigotto PP, Warburton D: Mir-17 family of micrornas controls fgf10-mediated embryonic lung epithelial branching morphogenesis through mapk14 and stat3 regulation of e-cadherin distribution. Dev Biol. 2009, 333 (2): 238-250. 10.1016/j.ydbio.2009.06.020.
Lee YB, Bantounas I, Lee DY, Phylactou L, Caldwell MA, Uney JB: Twist-1 regulates the mir-199a/214 cluster during development. Nucleic Acids Res. 2009, 37 (1): 123-128. 10.1093/nar/gkn920.
Huynh TP, Muráni E, Maak S, Ponsuksili S, Wimmers K: UBE3B and ZRANB1 polymorphisms and transcript abundance are associated with water holding capacity of porcine M. longissimus dorsi. Meat Sci. 2013, 95 (2): 166-172. 10.1016/j.meatsci.2013.04.033.
Nielsen S, Scheele C, Yfanti C, Akerström T, Nielsen AR, Pedersen BK, Laye MJ: Muscle specific microRNAs are regulated by endurance exercise in human skeletal muscle. J Physiol. 2010, 588 (Pt 20): 4029-4037.
Zhang P, Chen X, Fan M: Signaling mechanisms involved in disuse muscle atrophy. Med Hypotheses. 2007, 69 (2): 310-321. 10.1016/j.mehy.2006.11.043.
Jackman RW, Kandarian SC: The molecular basis of skeletal muscle atrophy. Am J Physiol Cell Physiol. 2004, 287 (4): C834-C843. 10.1152/ajpcell.00579.2003.
Honikel KO: Wasserbindungsvermögen von Fleisch. Mitteilungsblatt der BAFF. 1986, 6: 7150-7154.
Naraballobh W, Chomdej S, Murani E, Wimmers K, Ponsuksili S: Annotation and in silico localization of the Affymetrix GeneChip Porcine Genome Array. Arch Tierz. 2010, 53: 230-238.
Ponsuksili S, Du Y, Murani E, Schwerin M, Wimmers K: Elucidating molecular networks that either affect or respond to plasma cortisol concentration in target tissues of liver and muscle. Genetics. 2012, 192 (3): 1109-1122. 10.1534/genetics.112.143081.
Storey JD, Tibshirani R: Statistical significance for genome-wide experiments. PNAS. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.
Zhang B, Horvath S: A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005, 4: 17-
Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.
Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4 (1): 44-57.
Rehmsmeier M, Steffen P, Höchsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA. 2004, 10: 1507-1517. 10.1261/rna.5248604.
Lewis BP, Burge CB, Bartel DP: Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets. Cell. 2005, 120: 15-20. 10.1016/j.cell.2004.12.035.
The authors thank Annette Jugert and Joana Bittner for excellent technical help. This research was supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG; PO753/6-1 and WI1754/14-1).
The authors declare that they have no competing interests.
SP and KW conceived the study. SP analysed the microarray data and drafted the manuscript; YD made the WGCNA; FH helped in target prediction analysis. PS, EM, and MS helped in sampling and data collection and drafting the manuscript; KW discussed and contributed to data interpretation and helped in drafting the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Ponsuksili, S., Du, Y., Hadlich, F. et al. Correlated mRNAs and miRNAs from co-expression and regulatory networks affect porcine muscle and finally meat properties. BMC Genomics 14, 533 (2013). https://doi.org/10.1186/1471-2164-14-533
- Gene Ontology
- Meat Quality
- miRNA Family
- Carcass Trait
- Drip Loss