Gene expression patterns associated with peroxisome proliferator-activated receptor (PPAR) signaling in the longissimus dorsi of Hanwoo (Korean cattle).
Intramuscular fat (IMF or marbling) is deposited adipose tissue within muscle fibers in skeletal muscle. Marbling plays an important role in determining the quality of meat and therefore its economic value. Identifying the genes that regulate marbling may therefore help to optimize meat production in the livestock industry. Over the last few decades, studies have attempted to identify the genes that cause phenotypic differences in marbling using quantitative trait loci mapping, genome-wide association studies, microarrays, and proteome analyses. These results have demonstrated that marbling is controlled by multiple genes; i.e., it is a quantitative trait. Therefore, an analysis of multiple genes, especially genes related to adipogenesis or lipid metabolism, may lead to a better understanding of the effects of these processes on the phenotypes of marbled skeletal muscle. Recently, He et al. (2013) reported that 13 tag single nucleotide polymorphisms in the peroxisome proliferator-activated receptor (PPAR) signaling pathway are associated with porcine meat quality traits. On investigating fat deposition and muscle growth, seven single nucleotide polymorphisms in the reverse cholesterol transport pathway were found to be significantly associated with meat quality phenotypes in a Wagyu x Limousin reference population (Daniels et al., 2010). We recently revealed that three mitogen-activated protein kinase pathway-related genes participated in regulating the IMF content of Hanwoo, Korean native cattle (Kim et al., 2011). The skeletal muscle genes within the glycolytic pathway, which include protein kinase adenosine monophosphate-activated gamma 3 non-catalytic subunit (PRKAG3), phosphoglycerate mutase 2 (PGAM2), and pyruvate kinase muscle 2 (PKM2), have been described as markers of back fat thickness and average daily gain in Italian large white pigs (Fontanesi et al., 2008). These genes are therefore of interest with regard to understanding the factors that affect the meat quality of livestock. The PPAR signaling pathway, which is involved in lipid metabolism, is known to be an important biological pathway in determining meat quality in mammals. Peroxisome proliferator-activated receptor gamma, the primary transcription regulator in PPAR signaling, is reported to be a key factor in controlling the transcription of many genes that are involved in biological pathways based on adipogenesis (Fernyhough et al., 2007). It remains to be determined which genes of the PPAR signaling pathway play a role in livestock meat quality and more specifically, phenotypic differences in marbling.
The objective of this study was to identify genes that are differentially expressed in low- and high-marbled muscle of Korean cattle (Hanwoo) and to investigate their potential as biomarkers pertaining to the IMF content of beef.
MATERIALS AND METHODS
Animals and sample preparation
Twenty samples of low- and high-marbled m. longissimus from Hanwoo (n = 10 animals in each group) steers were obtained from the junction between the 11th and 12th lumbar vertebrae within 30 min of slaughter. These animals were bred under the same feeding conditions at the Hanwoo Experiment Station of the National Institute of Animal Science (NIAS) in Korea. The selected tissues were placed in liquid nitrogen, ground to a fine powder using a mortar, and stored at -80[degrees]C. All experimental procedures and the care of the animals were conducted in accordance with the guidelines of the Animal Care and Use Committee of the NIAS in Korea. The fat content of the m. longissimus muscle samples was analyzed using the methods established by the Association of Official Analytical Chemists (Table 1).
Quantitative real-time polymerase chain reaction and statistical analysis
The mRNA levels were analyzed by quantitative realtime polymerase chain reaction (qRT-PCR) using genespecific primer sets (Table 2). Total RNA was prepared from each tissue sample (100 mg) using TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA). For qRT-PCR analysis, 2 pg RNA was reverse transcribed in a 20 pL reaction volume using random primers (Promega, Madison, WI, USA) and reverse transcriptase (SuperScript II Reverse Transcriptase, Invitrogen Life Technologies). The reactions were incubated at 65[degrees]C for 5 min, 42[degrees]C for 50 min, and then at 70[degrees]C for 15 min to inactivate the reverse transcriptase.
The qRT-PCR was performed using the 2x Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) with the 7500 Real-Time PCR System (Applied Biosystems, USA) using 10 [micro]M of each primer (Table 2). The PCR was performed for 2 min at 50[degrees]C and 10 min at 95[degrees]C, followed by 40 cycles of 95[degrees]C for 10 s, and then 60[degrees]C for 1 min. Following amplification, melting curve analysis was conducted to verify the specificity of the reactions. The end point used in the qRT-PCR quantification, Ct, was defined as the PCR threshold cycle number. The [DELTA]Ct value was determined by subtracting the [beta]-actin Ct value for each sample from the target Ct value. The gene expression stability value of the [beta]-actin gene was less than 0.05, which met the stability requirement to be a housekeeping gene, confirming its suitability as an internal housekeeping gene in this experiment. The qRT-PCR data were used to calculate the normalized expression values ([2.sup.-[DELTA]Ct]) for the statistical analysis.
To examine the association between IMF content and gene expression levels, a statistical analysis was conducted using the analysis of variance (ANOVA) model. This resulted in the following equation: Fat = [mu] + expression + age + residual. In this equation, Fat denotes the IMF content (%) of each animal, [mu] is the overall mean, Expression is a normalized gene expression value ([2.sup.-[DELTA]Ct]), and Age is the slaughtering age (months) as a covariate. To determine the major patterns and relationships in the gene expression data, we performed a principal component analysis (PCA) of the genes.
Construction of the gene co-expression network and direct interaction relationships
In gene co-expression networks, we refer to nodes as genes whose degrees indicate the number of links connected to the node. We used expression values for 48 genes and evaluated pairwise correlations between the gene expression profiles of each pair of genes using Pearson's correlation coefficients (denoted as r). The network-encoded gene coexpression data were recorded as binary information (connected = 1, unconnected = 0) using a "hard" threshold. The degree of a node is the number of connections or edges the node has to other nodes. The degree distribution of a network has a generalized power-law form of p(k) ~ [k.sup.-r], which is the defining property of a scale-free network. We investigated hard thresholding with the power adjacency function and selected a power of tau ([tau]) = 0.65 using the weighted gene co-expression analysis (WGCNA) package in R (Langfelder and Horvath, 2008).
To characterize the gene co-expression network and the direct interaction relationships among the differentially expressed genes, we used the Pathway Studio software package (version 6.0, Stratagene, La Jolla, CA, USA) program (Nikitin et al., 2003). To identify the gene co-expression network, we used a p-value cut-off of 0.05, a Pearson's correlation (r) cut-off of 0.65 according to the result of WGCNA, with the removal of 5% of the genes with the most stable expression. To characterize the overall gene co-expression network topology, we investigated the node degree (or connectivity), betweenness centrality (BC), and closeness centrality (CC) using Cytoscape (Shannon et al., 2003) The node degree is the number of connections or edges the node has to other nodes. A node with high BC has considerable influence over what flows in the network. The BC can serve as a global property, as it is a useful indicator for detecting bottlenecks in a network. After experimental validation, we also examined the relationships and the effects of regulations (positive, negative, and unknown) among the significantly differentially expressed genes based on data in the literature.
RESULTS AND DISCUSSION
Gene expression and statistical analysis using principal component analysis and analysis of variance
Ten steers with low IMF content levels (low-marbled group) and 10 steers with high IMF content levels (highmarbled group) were used in an ANOVA model to identify an association between the expression of mRNA and the IMF content. The average IMF content in m. longissimus tissue for the low-marbled group was 7.44 [+ or -] 0.46%; this value for the high-marbled group was 23.54 [+ or -] 0.88% (Table 1). As shown in Table 1, the IMF content (p < 0.01) was significantly different between groups.
To identify the genes related to marbling, we investigated 48 genes within the PPAR signaling pathway from the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (bta03320). Table 2 shows the gene expression levels of 48 genes. The qRT-PCR results showed that 30 genes had higher expression levels in the high-marbled group than in the low-marbled group. The effect of 48 genes on the IMF was estimated using an ANOVA model; the least-square means were also estimated using the Student's t-test to determine if there was a significant difference between the marbled groups. The expression levels of 10 genes: retinoid X receptor alpha (RXRA), peroxisome proliferator-activated receptor gamma (PPARG), phospholipid transfer protein (PLTP), stearoylCoA desaturase (SCD), nuclear receptor subfamily 1 group H member 3 (NR1H3), fatty acid binding protein 3 (FABP3), carnitine palmitoyltransferase II (CPT2), acyl-Coenzyme A dehydrogenase long chain (ACADL), acyl-Coenzyme A oxidase 2 branched chain (ACOX2), and fatty acid binding protein 4 (FABP4), showed significant differences in gene expression between the low- and high-marbled groups (p < 0.05) (Table 2). Of these, nine genes showed higher expression in the high-marbled group, and one gene, FABP3, showed a low gene expression pattern in the high-marbled group.
FABP3 is involved in the uptake, transport, and metabolism of fatty acids (FAs) and is mainly expressed in cardiac and skeletal muscle. Ovilo et al. (2002) reported that the FABP3 gene had important effects on the IMF content. Other studies have also revealed that FABP3 is positively associated with the IMF content (Li et al., 2007; Zhao et al., 2009). However, Cui et al. (2012) reported that the expression level of the FABP3 gene in two chicken breeds was lower in later growth stages. Higher IMF and triglyceride (TG) levels are deposited later in growth compared with earlier growth stages in the pectoralis major tissue. Recently, Ye et al. (2014) also reported that the FABP3 gene was down-regulated in higher IMF content muscle. In the present study, the FABP3 gene expression pattern was lower in the high-marbled group than in the low-marbled group of Hanwoo cattle. We therefore suggest that lower expression of the FABP3 gene in the longissimus muscle is not a result of fat deposition and is most likely regulated by other biological factors. However, further investigations are needed to confirm the relationship between the FABP3 gene and marbling traits in skeletal muscle and adipocytes.
The PCA is a useful tool for data simplification and for the visualization of relationships. Therefore, we applied PCA to the gene expression dataset. The first two principal components (PCs) explained approximately 73.3% of the total variance, allowing most of the information to be visualized in two dimensions. The PCA indicated that the most important pattern of gene expression (PC1, accounting for 60% of the variance in the data) was associated with differences in the IMF levels. Individual samples were clearly partitioned into two separate groups: high- and lowmarbled groups based on PC1. In this analysis, the first PC illustrated a link between the nine differentially expressed genes, which had a positive relationship according to PC1 (Figure 1). Whereas, the FABP3 gene had a negative relationship according to PC1. Our experimental results suggested that these genes warrant further investigation as metabolic indicators of marbling in the bovine model.
Construction of a gene co-expression network based on gene expression patterns
The nodes represent 48 genes within the PPAR signaling pathway, and the links between the nodes represent the associations between expression profiles across all samples. The absolute value of Pearson's correlation coefficient was calculated for all pairwise comparisons. This correlation matrix was transformed into a matrix of adjacency using a "hard" threshold ([tau]). When the hard threshold is 0.65, the plot between [log.sub.10] (p(k)) and [log.sub.10] (k) shows that the network mostly follows a scale-free topology ([R.sup.2] = 0.54, slope = -0.45). To find clusters of highly correlated genes, we used the Pathway Studio program. We found 18 correlated genes and identified four distinct gene co-expression sub-clusters (Figure 2). The overall trend observed in the large co-expression pattern was also in agreement with the PCA result showing 10 significantly differentially expressed genes. From the PCA result, the significantly up-regulated genes according to PC1 could be divided into two clusters. The first cluster contained RXRA, PLTP, FABP4, and PPARG. The second contained ACOX2, CPT2, SCD, NR1H3, and ACADL. The gene co-expression network also showed a similar pattern. The genes in the PPARG-FABP4-PLTP group were interrelated. The ACOX2-CPT2 group also showed connections. This suggested that GK plays a role in connecting the PPARG, FABP4, PLTP, ACOX2, and CPT2 groups. The main cluster was large with a mean degree of 2.8, and it consisted of 10 genes, all of which were upregulated in the high-marbled group. This cluster was connected by nine target genes, including one transcription factor, PPARG. The nine target genes were acyl-CoA synthetase long-chain family member 6 (ACSL6), apolipoprotein A-I (APOA1), PLTP, GK, CPT2, aquaporin 7 (AQP7), acyl-Coenzyme A oxidase 1 palmitoyl (ACOX1), ACOX2, and FABP4. From the gene correlations, the FABP4-PPARG pair was highly correlated in the cluster (r = 0.92). PPARG and FABP4 are reportedly key regulators in the PPAR signaling pathway associated with adipogenesis and marbling traits in beef cattle. Therefore, the expression levels of these genes in muscle are reported to be positively correlated with marbling (Michal et al., 2006; Ahn et al., 2014). This is consistent with the fact that PPARG is directly linked to an increase in FABP4 among its targets or coactivators.
From the network topology analysis, the GK gene has also been shown to mediate the expression levels of genes involved in the network with the largest BC and CC values (BC = 0.35, CC = 0.69). The BC value is an indicator of a global central node. The effect of removing nodes with large BC values is similar to that of removing hub nodes because large BC nodes are correlated with hub nodes. The BC also implies that a site is relatively central with regard to all other sites. Thus, such sites are advantageously located to act as intermediaries. Therefore, we investigated the communication between genes and confirmed that hub and large-BC genes are located in the topological center of the network by calculating the CC in the network. The GK is directly connected to the following genes: ACSL6, ACOX1, ACOX2, AQP7, and PPARG. It was found that GK is functional in muscle cells and contributes to intramuscular TG synthesis for fat storage. Rahib et al. (2009) suggested that GK affects transcription factors that are important in muscle differentiation from a transcriptomic analysis of GK in mouse skeletal muscle. This gene showed no significant differences in expression between the high- and low-marbled groups. However, it may operate at the interface between fat and carbohydrate metabolism. Moreover, its deletion leads to alterations in the expression levels of genes involved in lipid and carbohydrate metabolism, PPARG, AQP7, or fatty acid synthase (Rahib et al., 2007). Therefore, GK may contribute to the cooperation between various genes and PPARG as a master regulator of fat and glycerol metabolism within the PPAR pathway.
The functional relationships between the significantly differentially expressed genes
As a candidate pathway involved in fat metabolism and development, KEGG analysis indicated that the PPARs of the PPAR signaling pathway in cattle trigger the induction of several downstream biological pathways, such as the pathways for lipid metabolism, adipocyte differentiation, and gluconeogenesis. Our results showed that eight genes were involved in the downstream pathways as transcription regulators, with the exception of PPARG and RXRA. Therefore, we categorized the significantly differentially expressed genes into the downstream pathways and investigated the relationships between these genes.
Lipid metabolism: The following genes were categorized into the lipid metabolism pathway: PLTP (lipid transport); SCD (lipogenesis); liver X receptor a, nuclear receptor subfamily 1, group H, member (LR1H3) (cholesterol metabolism); FABP3 (fatty acid transport); CPT2, ACADL, and ACOX2 (fatty acid oxidation); and FABP4 (adipocyte differentiation). These genes may be responsible for the high IMF content of Hanwoo.
The liver X receptor [alpha] (LR1H3 or LXR[alpha]) is a member of the nuclear receptor family of transcription factors and is involved in controlling cholesterol and fatty acid metabolism. This process is governed by sterol regulatory element-binding protein-1c gene (SREBP-1c), which regulates the acetyl-CoA carboxylase, fatty acid synthase, and SCD genes involved in the pathway of fatty acid biosynthesis (lipogenesis). The synergy between LR1H3 and SREBP-1c increases the synthesis of TG, which accumulates as cytoplasmic fat droplets and thus increases their secretion by PLTP into plasma in very-low-density lipoproteins (VLDL) (Repa et al., 2000). VLDL-TG production is directly related to free fatty acids, which are taken up by hepatocytes or myocytes for re-esterification or metabolism through [beta]-oxidation occurring in the mitochondria or in peroxisomes. Recent reports have also shown that activated LR1H3 stimulates adipocyte differentiation through the induction of PPARG expression (Seo et al., 2004). These reports demonstrated that the pathway of the liver X receptor a (LXR[alpha]) gene through the induction of PPARG was related to fat metabolism processes such as TG formation and fatty acid oxidation. Earlier reports have also shown that the LR1H3 (Yu et al., 2006), SCD (Cho et al., 2011), and PLTP (Seong et al., 2013) genes are related to fat metabolism and IMF. Interestingly, these findings are also in agreement with our results, in which three genes were found to be up-regulated in the high-marbled group as opposed to the low-marbled group, and expression levels were positively related to the IMF content in Hanwoo (Table 2). Moreover, we found direct regulatory relationships among these genes (Figure 3). The LR1H3-SCD and LR1H3-PLTP pairs showed a direct interactive relationship with the significantly differentially expressed genes identified in our study. Liver X receptor (LXR) ligands induce the expression of the PLTP gene and then regulate plasma lipoprotein metabolism (Laffitte et al., 2003). The lipogenic gene SCD is regulated partly by the insulin-sensitive transcription factor LXR[alpha] (Hebbachi et al., 2008). Other relationships, such as those between the PPARG-RXRA, PPARG-SCD, and RXRA-FABP4 pairs, have been thoroughly investigated in many studies (Chan et al., 2009).
In mammals, the oxidation of fatty acids occurs in two major systems: the mitochondrial and the peroxisomal [beta]-oxidation systems. The first step of the [beta]-oxidation cycle is catalyzed by acyl-coenzyme A dehydrogenases and acyl-CoA oxidases in the two respective systems. Prior to oxidation, fatty acids are transported into the mitochondria by carnitine palmitoyl transferases, they then enter the tricarboxylic acid cycle, where they are used by the electron transport chain to produce adenosine triphosphate (ATP). TG is the major component of marbled fat in skeletal muscle and is found mostly within intramuscular adipocytes. During the TG formation process, the first activation reaction is the consumption of two molar equivalents of ATP in a process that demands a high level of energy. Interestingly, the chemical inhibition of the mitochondrial respiratory chain, which produces ATP from products of the citric acid cycle and fatty acid oxidation, leads to reduced fat accumulation in the adipocytes of Caenorhabditis elegans, indicating that the mitochondria play a key role in fat synthesis (McKay et al., 2003). Acyl-coenzyme A dehydrogenase encodes one of the adipocyte fatty acidbinding proteins that belongs to a superfamily of lipidbinding proteins that are involved in the development of fatness traits. Recently, Jiang et al. (2014) found that ACADL genes were positively correlated with the adipocyte volume. We also found that the expression levels of the CPT2, ACADL, and ACOX2 genes were increased in the high-marbled group in the present study. Moreover, our previous reports found that gene expression within the mitochondrial oxidative phosphorylation pathway, involving NADH dehydrogenase 1, NADH dehydrogenase 2, cytochrome oxidase III, and ATP synthase F0 subunit 6 (ATP6), strongly correlated with IMF content (Kim et al., 2008). Taken together, these findings indicate that fatty acid oxidative metabolism in skeletal muscle is positively correlated with IMF content and fatty acid use for TG formation by providing the ATP.
Adipocyte differentiation: FABP4 is involved in adipocyte differentiation in the PPAR signaling pathway. FABP4 has been examined in many studies with regard to its function in fatness traits. Jurie et al. (2007) reported that TG content is also strongly correlated with FABP4 activity and oxidative enzyme activities at different IMF levels in cattle. Intracellular fatty acids in intramuscular adipocytes bind to FABP4, facilitating the transport of fatty acids from the plasma membrane to sites of fatty acid oxidation or esterification into the TG. In the present study, we found that the FABP4 gene showed a higher level of gene expression in the high-marbled sample, with the level differing significantly among the groups (Table 2). According to these results, we suggest that more marbled tissue exhibits more FABP4 content and uses more fatty acid for TG formation.
In conclusion, the 48 gene expression profiles examined in this study revealed 10 genes in the PPAR signaling pathway that were differently expressed between the high- and low-marbled groups, indicating a significant association with marbling. From the gene co-expression relationships and a network topology analysis, the main network was identified as 10 genes with a gene connectivity value of 2.8, all of which were up-regulated in the high-marbled group. Finally, we investigated the functional relationships according to the downstream pathways. We found that the major downstream pathway associated with the differently expressed genes was fatty acid oxidation. Activation of PPARG promotes terminal differentiation through the induction of a range of genes important for TG formation through fatty acid oxidation among downstream pathways involving CPT2, ACADL, and ACOX2.
Ahn, J., X. Li, Y M. Choi, S. Shin, S.-A. Oh, Y Suh, T H. Nguyen, M. Baik, S. Hwang, and K. Lee. 2014. Differential expressions of G0/G1 Switch gene 2 and comparative gene identification-58 are associated with fat content in bovine muscle. Lipids 49:1-14.
Chan, S. S. Y, L. J. Schedlich, S. M. Twigg, and R. C. Baxter. 2009. Inhibition of adipocyte differentiation by insulin-like growth factor-binding protein-3. Am. J. Physiol. Endocrinol. Metab. 296:E654-E663.
Cho, E. S., S. G. Kwon, J. H. Kim, D. H. Park, T. W. Kim, J. Nam, I. S. Jang, J.-S. Choi, W. Y. Bang, and C. W. Kim. 2011. Study for the expression of Adiponectin, fatty acid binding protein (FABP) 4, stearoyl-CoA desaturase (SCD) genes and the methylation of SCD promoter in porcine muscle and fat tissues. African J. Agric. Res. 6:6425-6431.
Cui, H. X., R. R. Liu, G. P. Zhao, M. Q. Zheng, J. L. Chen, and J. Wen. 2012. Identification of differentially expressed genes and pathways for intramuscular fat deposition in pectoralis major tissues of fast-and slow-growing chickens. BMC Genomics 13:213.
Daniels, T. F., X. L. Wu, Z. Pan, J. J. Michal, R. W. Wright Jr, K. M. Killinger, M. D. MacNeil, and Z. Jiang. 2010. The reverse cholesterol transport pathway improves understanding of genetic networks for fat deposition and muscle growth in beef cattle. PloS One 5(12):e15203.
Fernyhough, M., E. Okine, G. Hausman, J. L. Vierck, and M. V Dodson. 2007. PPARy and GLUT-4 expression as developmental regulators/markers for preadipocyte differentiation into an adipocyte. Domest. Anim. Endocrinol. 33:367-378.
Fontanesi, L., R. Davoli, L. Nanni Costa, F. Beretti, E. Scotti, M. Tazzoli, F. Tassone, M. Colombo, L. Buttazzoni, and V. Russo. 2008. Investigation of candidate genes for glycolytic potential of porcine skeletal muscle: Association with meat quality and production traits in Italian Large White pigs. Meat Sci. 80:780-787.
He, K., Q. Wang, Z. Wang, and Y. Pan. 2013. Association study between gene polymorphisms in PPAR signaling pathway and porcine meat quality traits. Mamm. Genome 24:322-331.
Hebbachi, A. M., B. L. Knight, D. Wiggins, D. D. Patel, and G. F. Gibbons. 2008. Peroxisome proliferator-activated receptor a deficiency abolishes the response of lipogenic gene expression to re-feeding. Restoration of the normal response by activation of liver X receptor a. J. Biol. Chem. 283:4866-4876.
Jiang, A.-A., M. Z. Li, H. F. Liu, L. Bai, J. Xiao, and X. W. Li. 2014. Higher expression of acyl-CoA dehydrogenase genes in adipose tissues of obese compared to lean pig breeds. Genet. Mol. Res. 13:1684-1689.
Jurie, C., I. Cassar-Malek, M. Bonnet, C. Leroux, D. Bauchart, P. Boulesteix, D. W. Pethick, and J. F. Hocquette. 2007. Adipocyte fatty acid-binding protein and mitochondrial enzyme activities in muscles as relevant indicators of marbling in cattle. J. Anim. Sci. 85:2660-2669.
Kim, N. K., D. Lim, S. H. Lee, Y M. Cho, E. W. Park, C. S. Lee, B. S. Shin, T. H. Kim, and D. Yoon. 2011. Heat shock protein B1 and its regulator genes are negatively correlated with intramuscular fat content in the Longissimus Thoracis muscle of Hanwoo (Korean cattle) steers. J. Agric. Food Chem. 59:5657-5664.
Kim, N. K., J. H. Lim, M. J. Song, O. H. Kim, B. Y Park, M. J. Kim, I. H. Hwang, and C. S. Lee. 2008. Comparisons of longissimus muscle metabolic enzymes and muscle fiber types in Korean and western pig breeds. Meat Sci. 78:455-460.
Laffitte, B. A., S. B. Joseph, M. Chen, A. Castrillo, J. Repa, D. Wilpitz, D. Mangelsdorf, and P. Tontonoz. 2003. The phospholipid transfer protein gene is a liver X receptor target expressed by macrophages in atherosclerotic lesions. Mol. Cell. Biol. 23:2182-2191.
Langfelder, P. and S. Horvath. 2008. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559.
Li, B., H. N. Zerby, and K. Lee. 2007. Heart fatty acid binding protein is upregulated during porcine adipocyte development. J. Anim. Sci. 85:1651-1659.
McKay, R. M., J. P. McKay, L. Avery, and J. M. Graff. 2003. C. elegans: A model for exploring the genetics of fat storage. Dev. Cell 4:131-142.
Michal, J., Z. W. Zhang, C. T. Gaskins, and Z. Jiang. 2006. The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu x Limousin [F.sub.2] crosses. Anim. Genet. 37:400-402.
Nikitin, A., S. Egorov, N. Daraselia, and I. Mazo. 2003. Pathway studio-the analysis and navigation of molecular networks. Bioinformatics 19:2155-2157.
Ovilo, C., A. Clop, J. L. Noguera, M. A. Oliver, C. Barragan, C. Rodriguez, L. Silio, M. A. Toro, A. Coll, J. M. Folch, A. Sanchez, D. Babot, L. Varona, and M. Perez-Enciso. 2002. Quantitative trait locus mapping for meat quality traits in an IberianxLandrace F2 pig population. J. Anim. Sci. 80:2801-2808.
Rahib, L., N. K. MacLennan, S. Horvath, L. C. Liao, and K. M. Dipple. 2007. Glycerol kinase deficiency alters expression of genes involved in lipid metabolism, carbohydrate metabolism, and insulin signaling. Eur. J. Hum. Genet. 15:646-657.
Rahib, L., G. Sriram, M. K. Harada, J. C. Liao, and K. M. Dipple. 2009. Transcriptomic and network component analysis of glycerol kinase in skeletal muscle using a mouse model of glycerol kinase deficiency. Mol. Genet. Metab. 96:106-112.
Repa, J. J., G. Liang, J. Ou, Y. Bashmakov, J.-M. A. Lobaccaro, I. Shimomura, B. Shan, M. S. Brown, J. L. Goldstein, and D. J. Mangelsdorf. 2000. Regulation of mouse sterol regulatory element-binding protein-1c gene (SREBP-1c) by oxysterol receptors, LXR[alpha] and LXRp. Genes Dev. 14:2819-2830.
Seo, J. B., H. M. Moon, W. S. Kim, Y S. Lee, H. W. Jeong, E. J. Yoo, J. Ham, H. Kang, M.-G. Park, K. R. Steffensen, T. M. Stulnig, J. Gustafsson, S. D. Park, and J. B. Kim. 2004. Activated liver X receptors stimulate adipocyte differentiation through induction of peroxisome proliferator-activated receptor y expression. Mol. Cell. Biol. 24:3430-3444.
Seong, J., J. D. Oh, H. K. Lee, G. J. Jeon, and H. S. Kong. 2013. Identification and analysis of phospholipid transfer protein polymorphisms and their association with marbling score in Hanwoo (Korean cattle). Genet. Mol. Res. 12:731-737.
Shannon, P., A. Markiel, O. Ozier, N. S. Baliga, J. T. Wang, D. Ramage, N. Amin, B. Schwikowski, and T. Ideker. 2003. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498-2504.
Ye, Y, S. Lin, H. Mu, X. Tang, Y Ou, J. Chen, Y. Ma, and Y Li. 2014. Analysis of differentially expressed genes and signaling pathways related to intramuscular fat deposition in skeletal muscle of sex-linked dwarf chickens. Biomed Res. Int. 2014, Article ID 724274.
Yu, M., B. Geiger, N. Deeb, and M. F. Rothschild. 2006. Liver X receptor alpha and beta genes have the potential role on loin lean and fat content in pigs. J. Anim. Breed. Genet. 123:81-88.
Zhao, S., L. J. Ren, L. Chen, X. Zhang, M. L. Cheng, W. Z. Li, Y. Y Zhang, and S. Z. Gao. 2009. Differential expression of lipid metabolism related genes in porcine muscle tissue leading to different intramuscular fat deposition. Lipids 44:1029-1037.
Dajeong Lim, Han-Ha Chai, Seung-Hwan Lee, Yong-Min Cho, Jung-Woo Choi, and Nam-Kuk Kim (1) *
National Institute of Animal Science, Rural Development Administration, Jeonju 565-851, Korea
* Corresponding Author: Nam-Kuk Kim. Tel: +82-54-429-7870, Fax: +82-54-429-7879, E-mail: firstname.lastname@example.org
(1) Experiment Research Institute, National Agricultural Products Quality Management Service, Gimcheon 740-870, Korea.
Submitted Oct. 17, 2014; Revised Dec. 29, 2014; Accepted Jan. 30, 2015
Table 1. Summary statistics of tissue samples for gene expression analysis Group Age (months) Carcass weight (kg) Low 27.1 [+ or -] 0.276 409.2 [+ or -] 10.588 High 29.5 [+ or -] 0.477 478.17 [+ or -] 13.618 ** Group Intramuscular fat content (%) Low 7.44 [+ or -] 0.463 High 23.54 [+ or -] 0.881 ** *, ** Significant differences (p < 0.05 and 0.01, respectively) between low-and high-marbled groups determined using an analysis of variance model. Table 2. Quantitative real-time polymerase chain reaction results for the 48 genes identified by analysis of variance Gene Full name PPARA Peroxisome proliferator-activated receptor alpha RXRA Retinoid X receptor, alpha RXRB Retinoid X receptor, beta RXRG Retinoid X receptor, gamma APOA1 Apolipoprotein A-I APOA2 Apolipoprotein A-II PPARD Peroxisome proliferator-activated receptor delta PPARG Peroxisome proliferator-activated receptor gamma APOA5 Apolipoprotein A-V APOC3 Apolipoprotein C-III PLTP Phospholipid transfer protein SCD Stearoyl-coa desaturase CYP7A1 Similar to cholesterol 7alpha-hydroxylase CW8B1 Cytochrome P450, family 8, subfamily B, polypeptide FADS2 Fatty acid desaturase 2 LR1H3 Nuclear receptor subfamily 1, group H, member FABP1 Fatty acid binding protein 1 SLC27A4 Solute carrier family 27 (fatty acid transporter) member 1 FABP3 Fatty acid binding protein 3 SLC27A1 Solute carrier family 27 (fatty acid transporter), member 4 CD36 CD36 molecule (thrombospondin receptor) OLR1 Oxidized low density lipoprotein (lectin-like) receptor 1 LPL Lipoprotein lipase ACSL6 Acyl-coa synthetase long-chain family member 6 ACSL5 Acyl-coa synthetase long-chain family member 5 ACSL4 Acyl-coa synthetase long-chain family member 4 ACSL1 Acyl-coa synthetase long-chain family member 1 CBT1B Carnitine palmitoyltransferase 1B CPT1C Carnitine palmitoyltransferase 1C CPT2 Carnitine palmitoyltransferase II ACAA1 Acetyl-Coenzyme A acyltransferase 1 ACADL Acyl-Coenzyme A dehydrogenase, long chain SCPX Sterol carrier protein 2 ACOX1 Acyl-Coenzyme A oxidase 1, palmitoyl ACOX2 Acyl-Coenzyme A oxidase 2, branched chain PGAR Angiopoietin-like 4 FABP4 Fatty acid binding protein 4, adipocyte CAP Sorbin and SH3 domain containing 1 PLIN1 Perilipin ADIPOQ Adiponectin, C1Q and collagen domain containing MMP1 Matrix metallopeptidase 1 (interstitial collagenase) UCP1 Uncoupling protein 1 (mitochondrial, proton carrier) ILK Integrin-linked kinase UBC Ubiquitin C PCK1 Phosphoenolpyruvate carboxykinase 1 PCK2 Phosphoenolpyruvate carboxykinase 2 GK Glycerol kinase AQP7 Aquaporin 7 [beta]-Actin Actin, beta Gene Primer sequence Forward Reverse PPARA ATTCCGAAACCGCGGATCTCAAGT AAAGGCGGGTTGTTGTTGGTCTTC RXRA TTCACCAAGCACATCTGTGCCATC GCTTGTCAATCAGGCAGTCCTTGT RXRB TGGCAAACGGCTATGTGCAATCTG CTTACGGATGGTGCGCTTGAAGAA RXRG AAGAGGGAAGCTGTGCAGGAAGAA TATCTTCATGGCCACTGTTGGCAC APOA1 AGCTCAACCTGAAACTCCTGGACA TTTCCAGGTTGTCCCAGAACTCCT APOA2 GAAGCTGCTTGCATTGACTGTGCT AAAGTAGGCCTTGGCTTGGGTTTG PPARD ATGTCGCACAATGCCATCCGTTTC TGAAGTTCTTCAGGTAGGCGCTGT PPARG CGATATCGACCAACTGAACC ACCGGTGATTTGTCTGTCGT APOA5 AGCGCATCCAGCAGAACCTAGA TTGAGCTGCTTCCTCCGCAT APOC3 TGGACAAGATGCAGGGCTATGTGA AGCTGAAGCTCTCAGTCATCCAGT PLTP ACGTGGCCTTCTCTGAGTTCTTCT TCATCAGGACGATGCTCCCAAAGT SCD CCAGAGGAGGTACTACAAACCTG AGCCAGGTGACGTTGAGC CYP7A1 GCAATGAAAGCGGCCACTGAAGAA TACTGTCTAGCACGGGCATGTTGT CW8B1 ATGCGCAGAGGAAGCTAGACTTTG AAGAGATGCTGGGATTGAGCTCCT FADS2 ACGATTACGGCCACCTCTCTGTTT AGGCACCCTTTAAGTGACCGATGA LR1H3 CTGAGTTTGCCTTGCTCATTGCCA ATGTGTGTTGCAGCCTCTCTACCT FABP1 GGAAATCGTGCAGAATGGGAAGCA TTCCTGCTGAACCACTGCCTTGAT SLC27A4 AGCCTGGTCAAGTTCTGTTCTGGA AGAAGAGTCGATCATCCATGCCCT FABP3 ACATGAAGTCACTCGGTGTCGGTT TTCCTGTCATCTGCTGTGGTCTCA SLC27A1 CAGGGCACCAACGACAAGAAGATT AGCTCGTCCATCACTAGCACATCA CD36 TCCAGACAACCACTGTTTCTGCAC GCATTTGCCAATGTCTAGCACACC OLR1 ACTTCACCGCCAGAACCTGAATCT GTGGGCATCCAAAGACAAGCAGTT LPL TACCCTGCCTGAAGTTTCCAC CCCAGTTTCAGCCAGACTTTC ACSL6 ATTGTGTGTTTCACGAGTGGCACG CAAACATGTGAGCCAAAGGCAGGA ACSL5 TGACAGTCTGTGGGACAAGCTCAT AGCCAGCAGTGCATTCTGTTTGAC ACSL4 GCCTTTGGCTCATGTGCTTGAGTT TTTGCTGGACTGGTCAGAGAGTGT ACSL1 CCCAAAGGAGCAATGATCACGCAT AGCTCCATGACAGAGCATCACACA CBT1B AGGCTGCTGAGAAGCACCAGAATA ACTTGGAGACCACGTAAAGGCAGA CPT1C TGCACTGGGACATTCCAGACAAGA TGTCTGAAGAGAGATGGCAGCGTT CPT2 AACTCTTCACCGATGACAAGGCCA ATGAGCCTGGATTTCAGAGGCACT ACAA1 TATGACATTGGCATGGCTTGTGGG TGCTCAGCCACATTCTCTGAGGTT ACADL TCCCTGAGATGGTTGCAGGGAAAT TGAGCCTCTCGATTTGTGATCGCT SCPX AAGCTGCCAGAAGGTGCTATGAGA ACACAGTCCCAGTGCTTCGTAAGT ACOX1 ATGCTGATGAAACATGCCCAGGTG TTCAGACTGATGCCTCACAGCACT ACOX2 ATGGTGGTGGTACGTGTGGACATT AGAGTTTCTGCTGCTGAGTCTGGT PGAR AAATGACCTCAGATGGAGGCTGGA TCTCCGAAGCCATCTTTGTAGGCT FABP4 CGTGGGCTTTGCTACCAG TGGTTGATTTTCCATCCCAG CAP ACAAGAATCTGCACTGGGAAGGGA CCATTGGCATCGGCATTGGAAGAA PLIN1 TTGGGTGGTGTCGAGAAGATGGTA AAGGTGTGTCTAGAGAGAGTGTTGGC ADIPOQ TGTGTATCGCTCAGCGTTCAGTGT AGTAGAGTCCCGGAATGTTGCAGT MMP1 TCTGGAGCAATGTCACACCCTTGA TCCGATACCTGCACCTGGTTGAAA UCP1 TTCTGCACAACGGTTCTGTCCTCT TTGGGCACACTTGTGTACTGTCCT ILK ATCACACACTGGATGCCATACGGA TCCAATGCAAACTTCACAGCCTGG UBC ACGTCAAGGCCAAGATCCAAGACA TTCACGAAGATCTGCATCCCACCT PCK1 AACCTGGCCATGATGAACCCTACT TCTTCTGGATGGTCTTGATGGCGT PCK2 AAGAAGTGCTTTGCCCTTCGCATC AGGTTTGTCTTGCCACAAGCACTG GK TGGAAGGTGTGCCAATATCTGGGT TCTCTGCCGAGTTTGTAAGCCACT AQP7 TGTCCTGGAAGAAGTTTCCCGTGT ACAGTGTCATGTGGTCAGGAAGGT [beta]-Actin ACCTGCGAAGATGTGTCTGAACCA TGTACACCGTCAGCTTGTTGAGGA Gene Expression (1) p value (2) Low High PPARA 0.0997 0.1261 0.308 RXRA 0.0883 0.1685 0.042 * RXRB 0.0376 0.0325 0.784 RXRG 0.0150 0.0251 0.153 APOA1 0.0042 0.0054 0.248 APOA2 0.0003 0.0004 0.158 PPARD 0.0762 0.0631 0.465 PPARG 0.0098 0.0511 0.021 * APOA5 0.0001 0.0001 0.647 APOC3 0.0004 0.0007 0.113 PLTP 0.0205 0.0658 0.001 ** SCD 0.8597 3.0974 0.002 ** CYP7A1 0.0003 0.0004 0.439 CW8B1 0.0013 0.0015 0.640 FADS2 0.0177 0.0126 0.450 LR1H3 0.0098 0.0183 0.034 FABP1 0.0025 0.0024 0.740 SLC27A4 0.0417 0.0473 0.304 FABP3 1.7208 1.0294 0.023 * SLC27A1 0.0002 0.0003 0.471 CD36 0.2479 0.2026 0.544 OLR1 0.0002 0.0003 0.268 LPL 0.0820 0.1676 0.226 ACSL6 0.0001 0.0002 0.317 ACSL5 0.0008 0.0006 0.320 ACSL4 0.0024 0.0023 0.663 ACSL1 0.0230 0.0178 0.436 CBT1B 0.0055 0.0050 0.932 CPT1C 0.0002 0.0003 0.467 CPT2 0.0020 0.0076 0.001 * ACAA1 0.0101 0.0075 0.184 ACADL 0.0011 0.0036 0.002 * SCPX 0.0113 0.0085 0.439 ACOX1 0.0164 0.0185 0.399 ACOX2 0.0044 0.0090 0.008 ** PGAR 0.0192 0.0172 0.905 FABP4 0.0669 0.4280 0.002 ** CAP 0.0753 0.0810 0.488 PLIN1 0.0572 0.0740 0.401 ADIPOQ 0.0001 0.0001 0.088 MMP1 0.0001 0.0001 0.107 UCP1 0.0003 0.0002 0.070 ILK 0.2481 0.1484 0.053 UBC 0.1505 0.2632 0.450 PCK1 0.0233 0.0148 0.127 PCK2 0.0246 0.0372 0.406 GK 0.0145 0.0177 0.349 AQP7 0.3712 0.4460 0.363 [beta]-Actin (1) Expression shown as the mean of the normalized expression value of each gene within the low- and high-marbled groups. (2) p value was calculated using analysis of variance. *, ** Significant effects on intramuscular fat contents (p < 0.05, 0.01, respectively).
|Printer friendly Cite/link Email Feedback|
|Author:||Lim, Dajeong; Chai, Han-Ha; Lee, Seung-Hwan; Cho, Yong-Min; Choi, Jung-Woo; Kim, Nam-Kuk|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Aug 1, 2015|
|Previous Article:||Y-single nucleotide polymorphisms diversity in Chinese indigenous horse.|
|Next Article:||Single nucleotide polymorphism in the coding region of bovine chemerin gene and their associations with carcass traits in Japanese Black cattle.|