Communities in English
Advertise on Engormix

Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria

Published: November 6, 2018
By: Byung-Whi Kong 1, Kentu Lassiter 1, Alissa Piekarski-Welsher 1, Sami Dridi 1, Antonio Reverter-Gomez 2, Nicholas James Hudson 2, Walter Gay Bottje 1. / 1 Department of Poultry Science, Center of Excellence for Poultry Science, University of Arkansas, Fayetteville, AR, USA; 2 CSIRO Livestock Industries, Queensland Bioscience Precinct, 306 Carmody Road, St.Lucia, QLD, Australia.

As feed represents 60 to 70% of the cost of raising an animal to market weight, feed efficiency (the amount of dry weight intake to amount of wet weight gain) remains an important genetic trait in animal agriculture. To gain greater understanding of cellular mechanisms of feed efficiency (FE), shotgun proteomics was conducted using in-gel trypsin digestion and tandem mass spectrometry on breast muscle samples obtained from pedigree male (PedM) broilers exhibiting high feed efficiency (FE) or low FE phenotypes (n = 4 per group). The high FE group had greater body weight gain (P = 0.004) but consumed the same amount of feed (P = 0.30) from 6 to 7 wk resulting in higher FE (P < 0.001). Over 1800 proteins were identified, of which 152 were different (P < 0.05) by at least 1.3 fold and ≤15 fold between the high and low FE phenotypes. Data were analyzed for a modified differential expression (DE) metric (Phenotypic Impact Factors or PIF) and interpretation of protein expression data facilitated using the Ingenuity Pathway Analysis (IPA) program. In the entire data set, 228 mitochondrial proteins were identified whose collective expression indicates a higher mitochondrial expression in the high FE phenotype (binomial probability P < 0.00001). Within the top up and down 5% PIF molecules in the dataset, there were 15 mitoproteome proteins up-regulated and only 5 down-regulated in the high FE phenotype. Pathway enrichment analysis also identified mitochondrial dysfunction and oxidative phosphorylation as the number 1 and 5 differentially expressed canonical pathways (upregulated in high FE) in the proteomic dataset. Upstream analysis (based on DE of downstream molecules) predicted that insulin receptor, insulin like growth receptor 1, nuclear factor, erythroid 2-like 2, AMP activated protein kinase (α subunit), progesterone and triiodothyronine would be activated in the high FE phenotype whereas rapamycin independent companion of target of rapamycin, mitogen activated protein kinase 4, and serum response factor would be inhibited in the high FE phenotype. The results provide additional insight into the fundamental molecular landscape of feed efficiency in breast muscle of broilers as well as further support for a role of mitochondria in the phenotypic expression of FE.

Funding provided by USDA-NIFA (#2013–01953), Arkansas Biosciences Institute (Little Rock, AR), McMaster Fellowship (AUS to WB) and the Agricultural Experiment Station (Univ. of Arkansas, Fayetteville).

The projected doubling of the human population by 2060 will require a 100% increase in food production coming from plants and animals with most of this increase coming from new technology and greater efficiency [1]. In animal production agriculture, feed is the highest input cost (50 to 70% of total) in raising an animal to market weight and feed costs can spike as observed during the severe drought in the US in 2012 [2]. Great strides have been made in animal agriculture production efficiency through the selection of animals for feed efficiency (FE, weight gained to feed consumed), feed conversion ratio (FCR, feed to gain), or residual feed intake (RFI, the actual amount of feed intake that is above or below predicated feed intake in a group of animals). These are all labor intensive but effective methods that require measuring feed intake and body weight gain on individual animals. Development of biomarker selection tools to be used in commercial breeding programs can contribute significantly towards increasing animal production efficiency.
Global gene expression studies have been conducted to gain a greater understanding of the cellular basis of FE in three different broiler models. Using individually phenotyped birds, gene expression has been investigated in; 1) duodenal tissue in broilers selected for low and high RFI [3], 2) breast muscle tissue in commercial broilers [4], and 3) breast muscle within a pedigree male (PedM) broiler line [5]. Based on the cDNA microarray results [5, 6], we hypothesized that gene expression in the low FE phenotype resulted from inherent gene expression that was modulated by mitochondrial ROS [7]. Higher oxidative stress was observed in the low FE PedM broiler phenotype in the form of higher mitochondrial reactive oxygen species (ROS) production and pervasive protein oxidation in isolated mitochondrial and whole tissue homogenates [8–10]. In contrast, Zhou et al. [4] hypothesized that increased oxidative stress in high FE commercial broilers was due to higher inflammatory gene expression (cytokines or myokines) associated with muscle remodeling and development.
Despite the potential of transcriptome analyses for discovering important genes and pathways associated with economically important production traits, gene transcript profiles provide only part of the picture as they may be subject to feedback control by their protein products. Potential drawbacks of transcriptome analyses include inconsistencies of corresponding protein expression levels [11–15] and lack of information on post-translational modifications (e.g. phosphorylation). Proteomic approaches are promising as they can circumvent some of the issues associated with transcriptomics. Various proteomic studies in poultry have recently emerged including reports on muscle [16, 17], eggs [18–20], and liver [21] as well as host-virus interactions [22–26]. Muscle proteomics in livestock has been extensively studied for physiological relevance associated with specific diseases, metabolic conditions, and production traits (See review by Picard et al. [27]). Most studies have relied on two dimensional gel electrophoresis (2DGE)- based protein identification methods. The 2DGE of crude tissue extracts usually underestimates the presence of low abundance elements, integral membrane proteins, regulatory proteins (e.g. transcription factors), and components with very high molecular masses in complex tissues, such as muscle [28] that contain very high levels of muscle structural subunits.
Shotgun proteomics based on high-resolution mass spectrometry (MS) allows quantification of thousands of proteins along with their modifications, localization, turnover, and interaction partners [29, 30]. Shotgun proteomics does not focus on specific proteins and thus offers a hypothesis-free and systems-wide analysis that complements antibody-based approaches. Advances in all areas of the MS-based proteomics workflow, including sample preparation, liquid chromatography (LC)-MS and computational analysis, have made it a preferred analytic tool to characterize protein dynamics in multiple functional dimensions; and proteomics based on LC-MS has essentially replaced previous tools such as 2DGE. The shotgun proteomics method has been used in poultry [22].
The major goal of this study was to conduct experiments in order to gain new insight into the cellular basis of FE by analysis of proteomic data obtained in PedM broilers exhibiting high or low FE phenotypes. We screened the data for patterns of differential expression (DE) and assessed the mitochondrial proteome in particular detail. We used DE and cluster analysis to identify groups of proteins expressed at different levels between the two groups of birds. Furthermore, a commercial software program (Ingenuity Pathway Analysis, IPA; Qiagen, Valencia, CA) was used to develop predictions of activation or inhibition of upstream regulators, molecules, and pathways of DE proteins. The analysis was conducted on the same set of samples that were used for microarray examination of global gene expression associated with FE [5], thus starting to develop a proteogenomic profile associated within this particular FE model.
Materials and Methods
Ethics Statement
The present study was conducted in accordance with the recommendations in the guide for the care and use of laboratory animals of the National Institutes of Health. All procedures for animal care complied with the University of Arkansas Institutional Animal Care and Use Committee (IACUC): Protocol #14012.
Animals and tissue (breast muscle) collection
Proteomic analysis was conducted on a subset of breast muscle samples previously used in global gene expression analysis [5]. Phenotyping of animals in this study followed procedures outlined previously [8]. Briefly, feed efficiency (FE, amount of body weight gain/amount of feed consumed) was determined between 6 and 7 wk of age on a group of 100 pedigree male (PedM) broilers housed in individual cages. Birds were provided access to feed and water ad libitum. All birds received the same corn-soybean based diet (20.5% protein, 3,280 kcal/kg) during the feed efficiency trial. From this group of 100 birds with the highest and lowest FE were selected. Feed efficiencies for the low- and high-FE groups (n = 4 per group) in the present study were 0.46 ± 0.01 and 0.65 ± 0.01, respectively. Higher efficiency was obtained by greater weight gain without a difference in feed intake during the week of FE phenotyping. The difference in mean FE (0.19) is identical to that reported in our initial investigation [8]. Birds were humanely killed, tissues were quickly removed, and flash frozen in liquid nitrogen, and stored at -80°C.
Protein extraction
Protein extraction methodology followed procedures by Iqbal et al. [31] with modifications. Briefly, total proteins were extracted from breast muscles of 4 birds per group (high and low FE), which were used in cDNA microarray [5]. Approximately 100 mg of frozen muscle tissue was minced and homogenized in 700 μl of ice-cold tissue lysis buffer containing 150 mM NaCl, 0.1% Triton X-100, 0.5% sodium deoxycholate, 0.1% SDS, 50 mM Tris-HCl pH 8.0, and protease and phosphatase inhibitors (Life Technologies, Grand Island, NY) with 0.9–2.0 mm stainless steel bullet blender beads (Next Advance, Averill Park, NY). The samples were centrifuged at 4°C for 20 min at 14,000xg and the supernatant was carefully collected, aliquoted and stored at -80°C until use. Protein concentration was determined by Bradford assay using a Synergy HT multi-mode microplate reader (BioTek, Winooski, VT). Extracted proteins were analyzed by SDS-PAGE gel electrophoresis (Life Technologies, Grand Island, NY) to assess protein quality by the presence of distinct protein bands (Fig 1).
Shotgun proteomics
Extracted individual proteins were subjected to shotgun proteomics analysis by in-gel trypsin digestion and tandem mass spectrometry (MS/MS) at the University of Arkansas Medical Science (UAMS) Proteomics Core Lab (Little Rock, AR). In total, 25 gel slices per sample were analyzed. Raw mass spectrometric data were analyzed by database searching using the Mascot (Matrix Science, Boston, MA) search engine and the UniProtKB (http://www.uniprot.org/ help/uniprotkb) database. Search results were compiled using Scaffold program (Proteome Software, Portland, OR).
Statistics and Bioinformatic Analyses
Data were normalized based on total spectral counts for each individual sample. Evidence of success in normalization of samples is provided in S1 Fig showing no differences in actin, myosin or tubulin expression that serve as standard internal house-keeping proteins. Expression of GAPDH, often used as an internal reference check, was down-regulated (-1.45 fold, P < 0.05) in the high compared to low FE group (See Table 1 below). To generate MA plots, we added a nominal 0.05 to all protein abundance values in order to handle the presence of 0’s, log2 transformed the data to stabilize the variance, and computed average (A) and minus (M) at the group level (high FE minus low FE) for all 1817 proteins (See Fig 2A). The A values were 0 justified by adding 4 to all values.
Phenotypic impact factor (PIF) and Hierarchical cluster analysis. To identify and interpret DE proteins, we employed two complementary approaches. Firstly, we identified the edges of the MA distribution of all proteins using a modified DE called Phenotypic Impact Factor (PIF) [32], taking the liberal nominal (P < 0.10) approach of identifying the 5% most up regulated and 5% most down regulated proteins. Each PIF is computed by multiplying the average abundance by DE for all proteins. This has a number of appealing numerical characteristics such as tracking the edge of the MA plot distribution and de-emphasizing lowly abundant proteins which are inherently noisier as they are close to the detection limit of the technology. These 182 extreme DE-PIF proteins were assessed for functional enrichment using the two list import function in GOrilla [33] which compares the functional contents of protein lists using hypergeometric statistics.
Next, we identified the top 40 proteins by PIF (20 up- and 20 down-regulated) using hierarchical cluster analysis. We imported the matrix with as many columns as birds (8) and as many rows as proteins (40), where each cell contains the log2 transformed abundance value for that protein and individual into PermutMatrix software [34], normalizing on rows. We applied hierarchical clustering on both rows and columns and exported the resulting dendrogram as a JPEG (Fig 3).
Bayes Methods and Ingenuity Pathway Analysis. A moderated t-statistic and its corresponding P-value was used based on empirical Bayes methods [35] for each peptide. Ingenuity Pathway Analysis (IPA; Qiagen, Valencia, CA; http://www.ingenuity.com) software was used for functional annotation of proteins, as well as canonical pathway analysis, upstream analysis, and network discovery. Upstream regulator analysis by IPA is based on prior knowledge of expected effects between transcriptional regulators and their target genes from published literature citations stored in the IPA program. The analysis determines how many known targets or regulators are within the user’s dataset, and compares the DE to the expected relationship in the literature. If the observed direction of change is mostly consistent with either activation or inhibition of the transcriptional regulator, then a prediction is made and an activation z score is generated that is also based on literature-derived regulation direction (i.e. “activating” or “inhibiting”). Activation z scores > 2.0 indicate that a molecule is activated whereas activation z scores of < -2.0 indicate that a target molecule is inhibited. Qualified predictions can also be made of high, moderate, or low absolute z scores of 1.7 or more, 1.5 to 1.7, or less than 1.5 and greater than 1.0, respectively. The p-value of overlap measures whether there is a statistically significant overlap between the dataset molecules and those regulated by an upstream regulator is calculated using Fisher’s Exact Test, and significance is attributed to p-values < 0.05.
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 1
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 2
Quantifying the mitoproteome
To estimate mitochondrial content between the groups, we downloaded the human mitoproteome from the Mitominer website (http://mitominer.mrc-mbu.cam.ac.uk/release-3.1/begin. do). This gave us a list of 1046 nuclear- and mitochondrial-encoded mitochondrial proteins of which we identified 228 matches in our proteomics data set. We generated a proxy for tissue mitochondria content by overlaying the MA data of the mitoproteome onto the MA plot and assessing its deviation from 50:50 equilibrium by binomial statistics (Fig 2A).
Results and Discussion
I. Global distribution of proteins and the mitochondrial proteome
The global distribution of DE for the 1817 proteins is illustrated in Fig 2A and 2B. The plot is symmetrical and centered on 0 (for DE and protein abundance). Names of proteins that are abbreviated can be found in Tables 1 or 2 or in the figure legends. As with all MA plots, the noise is greater for low abundant proteins on the left hand side because the data approaches the detection limit of the technology, and the spread of DE data is greater. For abundance (A) values less than 4.0, there are often missing values in a subset of samples, so these DE estimates can arguably be considered less reliable. Some prominent outlier proteins have been labeled (Fig 2A) and include the mitochondrial energy transfer proteins (SLC25A4 and VDAC2) upregulated in high FE, and the slow muscle isoforms (TPM2 and ACTC) down-regulated in high FE. The muscle fiber composition observation is consistent with an RNA sequencing analysis performed on the same samples (unpublished observations, manuscript in review) and observations of muscle whitening made in other FE selected agricultural species, such as pigs (e.g. [36]). Proteins expected to be very abundant in fast type IIB muscle fibers such as the glycolysis proteins GAPDH, PKM, and LDHA; and the fast fiber subunits MYH2 and MYOM2 are among the most abundant proteins on average across the entire data set. The CKM protein which transfers phosphate between ATP and creatine phosphate, a core energetic process in muscle cells, is also highly abundant in all samples. These results reflect the fact that broiler breast muscle is 100% type IIB fibers at a histological level [37]. The presence of some slow subunits in our data (and associated RNAseq data, not shown) reaffirms that fiber composition is more of a mosaic at the molecular level than it is at the whole fiber level
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 3
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 4
The extreme 10% from the entire dataset (comprising the 5% up- and 5% down-regulated PIF proteins analysis; i.e. 91 up- and 91 down-regulated, 182 PIF proteins total, See Fig 2B) in high FE were subject to functional enrichment analysis in GOrilla. The up-regulated 91 proteins were enriched for intermediate filament (FDR Q-value < 0.0001) based on the presence of DES, KRT3, KRT4, KRT5, KRT15, KRT75, LMNB2, DSP, SYNM and JUP). The downregulated list received no significant formal enrichment by this approach
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 5
The mitoproteome of 228 proteins was significantly skewed towards the High FE birds, indicating a higher tissue mitochondria content (binomial probability P < 0.00001; Fig 2A). These findings are consistent with the results of global gene expression analysis (RNAseq) that was carried out on the same set of breast muscle samples as the present proteomic dataset (unpublished observations, manuscript in review). The mitochondrial proteome expression elevation was likely not the result of increased numbers of mitochondria as mitochondrial DNA (a proxy for mitochondrial numbers) was the same in both HFE and LFE samples (S2 Fig). Additionally, Kim et al. [38] provided clear evidence that mtDNA copy number or mitochondrial mass was not associated with increased mitochondrial numbers and/or increased mitochondrial protein expression in chicken embryo fibroblasts. Within the top up or down 5% of phenotype impact factors (PIF), there were 15 up-regulated mitoproteome proteins (SLC25A4, KRT5, IDE, NDUFV2, TUBA8, VDAC2, ACSL1, NDUFB5, COX4I1, UQCRFS1, NARS, PTRH2, PPA2, C1QBP, PRDX3) and only 5 down-regulated proteins (COQ9, MRPL44, ACSF2, ALDH9A1, SOD1) (Fig 2B).
In the cluster dendrogram (Fig 3), red and green represent up- and down-regulated expression in the high FE phenotype, respectively. Cluster analysis for the extreme 40 DE proteins (determined by PIF analysis) clearly discriminates the 8 birds into the FE groups they originated from. Two birds, B1013 and G1004, have a slight tendency to be outliers but are still more similar to members of their own group than to the other FE group. Among the proteins present in this nominal top 40 are mitochondrial proteins (NDUFV2, SLC25A4, IDE and KRT5), non-mitochondrial metabolism (PHKG1), and muscle sarcomeric structure (TPM1, TPM2, ACTA1 and CAPZB). As shown in Fig 3, clustering proteins by patterns of co-expression across the 8 samples produces some functional groupings of cytoskeletal keratins (KRT3, KRT5, and KRT14), muscle structural isoforms (TPM1 and TPM2), and 2 different protein isoforms of ARF (ARF1 and ARF4). Not surprisingly, different protein fragments of the same protein are always strongly co-expressed. This latter observation can be taken as a form of technical validation for our approach.
II. Pathway Analysis
The 152 DE proteins ( 1.3 fold, P < 0.05) in breast muscle of high and low FE PedM broilers are presented in Table 1. The proteins in Table 1 are listed in order of the greatest to the least differentially expressed in the high FE compared to low FE phenotype. This DE ranking is related to the PIF metric output described above, but not identical because only fold difference (without multiplying by protein abundance) is reported. The DE proteins in Table 1 were organized by Ingenuity Pathway Analysis (IPA) into pathways and functional groupings, as well as predicted activation or inhibition of upstream regulators. Because of the amount of information on protein expression available in this dataset, we will focus on a) canonical pathways associated with mitochondrial function and oxidative phosphorylation, and b) selected upstream regulators.
Mitochondrial function—oxidative phosphorylation. Mitochondrial function and oxidative phosphorylation are listed as the first and fifth top canonical pathways, respectively, in the proteomic dataset (Table 3). This is consistent with the skewing of the mitoproteome towards the high FE birds as indicated above. Projection of the up-regulated mitochondrial proteins in the high FE phenotype (from Table 1) onto the canonical pathway of oxidative phosphorylation is presented in Fig 4A. From these data, the IPA program predicted that Complex I, III, IV and V activities would be greater (shown in orange in Fig 4B) in breast muscle of the high FE compared to the low FE phenotype in this study. The IPA program also predicted that Complex II would be activated in the high FE phenotype using a 1.2 fold expression cutoff (P <0.05, data not shown). These predicted increases in respiratory chain activities are in agreement with the general increase in mitoproteome expression discussed previously. The predicted increases in complex activities also concurs with results from enzymatic analysis conducted previously on muscle, liver and duodenal tissue obtained from the same pedigree male broiler line indicated that there was a general increase in complex activities in the high compared to low FE phenotype [10]. Similarly, complex activities were higher in the lumbar multifidus (LM) muscle of Ghezel male lambs with low RFI (more feed efficient) compared to those with high RFI (less efficient) [39]. We hypothesized that the lower complex activities in the low FE phenotype were due to a pervasive protein oxidation from mitochondrial ROS production observed in tissue homogenate and mitochondrial fractions [9, 10].
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 6
Two up-regulated proteins in the high FE breast muscle integral to mitochondrial function are adenine nucleotide translocase (ANT or SLC25A4) and voltage dependent anion channel 1 and 2 (VDAC1 and VDAC2) (Table 1). ANT and VDAC2 were also identified as PIF molecules as described above. With the exception of a few lipophilic compounds (e.g. acetaldehydes, short chain fatty acids and molecular oxygen) all metabolites that enter or leave mitochondria must pass through the VDAC channel in the outer mitochondrial membrane [40]. The channels play an important role in transmitting redox information between the mitochondria and nucleus [41] in part by the ability of superoxide to move out of the mitochondria through VDAC channels [42]. The combined activities of ANT and VDAC are responsible for shuttling ADP and ATP between the mitochondrial matrix and the cytosol thus facilitating cellular ‘work’ while providing substrate for ATP synthase that uses proton motive force generated by the respiratory chain to produce ATP from ADP and Pi.
The ANT-VDAC coupling in conjunction with mitochondrial creatine kinase (mtCK) facilitates the transfer of high energy phosphate bonds from ATP to creatine to form phosphocreatine (see review [43]). In the model presented by Schlattner et al. [43], mtCK, located in the intramembranous space, has contact sites with both VDAC on the outer membrane and ANT on the inner mitochondrial membrane. This arrangement facilitates direct channeling of ATP to mtCK that catalyzes the transfer of phosphate groups to creatine to form phosphocreatine that is released to the cytosol. A large cytosolic pool of phosphocreatine can ‘....then be used as a temporal buffer to maintain constant global and local ATP/ADP ratios over a wide range of workloads...and spatial energy buffering; i.e. for an energy shuttle between ATP-providing or consuming processes.” [43].
There are at least four subunit isoforms of CK expressed in vertebrate tissues: two cytosolic forms, CKM (muscle) and CKB (brain), and two mitochondrial (mt) CK isoforms. In vivo, CKM and CKB isoforms combine to provide three typical dimeric cytosolic isozymes: MM-, MB- and BB-CK [44]. In the proteomics dataset, CKB was up regulated 8.9 fold whereas CKM was down regulated -1.47 fold (Table 1) compared to the low FE phenotype; the mtCK isoform was not detected. Quest et al [45] suggested that tissue-specific dimerization of B-CK monomers in chicken tissues could represent a mechanism for specifying the intracellular distribution of this cytosolic isoform that in turn could influence energy transfer by phosphocreatine. How this differential expression of CKM and CKB relates to energy transduction in the high and low FE phenotypes is not apparent at this time, but further studies are warranted. We have also observed that mitochondrial CK (CKMT1A) was up-regulated in high FE in the RNA sequencing data being analyzed concurrently (manuscript in preparation).
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 7
Another up-regulated mitochondrial protein in the high FE phenotype was peroxiredoxin III (PRDX3, 2.85 fold, Table 1). In mammals, PRDX3 is located exclusively in the mitochondria and functions like a GSH peroxidase in converting hydrogen peroxide to water [46–48] and maintains mitochondrial mass and membrane potential [49]. Increased PRDX3 expression in high FE mitochondria could, therefore, contribute to lower mitochondrial ROS production and lower protein oxidation [10] and higher expression of electron transport chain proteins in the present study (Fig 4).
The increased expression of mitochondrial electron transport chain proteins shown in Fig 4 contrasts with recent reports in 1) female Large White pigs (115 kg, from lines divergently selected for high or low RFI [36]), and 2) in the top 5 and bottom 5 of a group of 238 castrated Yorkshire purebred boars individually phenotyped for RFI [50]. Using a cDNA microarray, Vincent et al. [36] reported that mRNA expression of genes encoding for protein subunits for Complex I (ND1, ND2, ND4, NDUFB9), and isocitrate dehydrogenase, were down-regulated in the low RFI pig. It is also contrary to the reduction in slow subunits (TPM2 and ACTC) we observed in the high FE birds [5, 6], as it is well established across all species that slow, aerobic fibers have a higher mitochondrial content than fast, glycolytic fibers. Vincent et al. [36] also reported lower ATP5A1 (a subunit of Complex V, ATP synthase) and, in contrast to what was mentioned above, lower mitochondrial CK expression in the low RFI females compared to the high RFI females. There are at least three possible explanations for differences between the present study and Vincent et al. [36]; 1) species difference (swine vs broiler), 2) sex (female vs. male), and 3) protein expression does not always follow gene expression differences. With respect to point 3 above, projection of the DE genes in the cDNA microarray dataset from Kong et al. [5] onto the oxidative phosphorylation canonical pathway using the overlay function in IPA, revealed up-regulation of genes and prediction of activation of Complex I but down-regulation of genes and prediction of inhibition of Complex III and IV activities (data not shown).
Using qPCR, Jing et al. [50] reported decreased mRNA expression of 16 of 17 mitochondrial DNA encoded proteins in low RFI castrated boars including four from Complex I and one from Complex III along with decreased expression of the 5’AMP activated protein kinase gamma 2 (AMPKγ2) and peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PPARPGC1A, or PGC1-α. PGC1-α is a master regulator of mitochondrial biogenesis [51–53] and AMPK is an energy sensing molecule that responds to increased AMP to ATP levels by increasing energy production (e.g. fatty acid β oxidation, glycolysis, oxidative phosphorylation) and decreasing energy consumption reactions (e.g. gluconeogenesis, fatty acid synthesis) [54–57]. Thus, the decreased gene expression of mitochondrial DNA-encoded proteins reported by Jing et al. [50] is congruent with the decreased expression of AMPK and PGC1-α. Similarly, increased expression of respiratory chain proteins in the current study is congruent with increased AMPK mRNA expression in the high FE broiler [6] as well as the predicted activation of PRKAA2 (AMPKα2, catalytic subunit of AMPK), and moderate activation of PPARGC1α (PGC1-α) derived from the proteomic dataset discussed below. However, a fundamental difference (besides species) between Jing et al [50] and the current proteomics study is the use of castrated males vs. intact males in the present study. Broilers typically begin to show sexual dimorphism (increased comb and wattle size and heavier body weights) around 4 to 5 wk of age; two weeks prior to when FE phenotyping occurred in our FE model. Shifts in steroid hormone levels between castrated and intact males could have a number of effects on nuclear and mitochondrial gene and protein expression; especially since thyroid and steroid hormone receptors are found in mitochondria as well as the nucleus [58].
Upstream Regulators. The IPA program was used to predict activity of upstream regulators based on DE proteins in the high and low FE phenotypes. The predictions are based on activation z-score and P-value of overlap generated from literature based citations for each differentially expressed molecule/protein. A P overlap value of < 0.05 and an activation z-score of > 2.0 or < -2.0 indicates that a specific molecule is predicted to be activated or inhibited, respectively. It is also possible to provide qualified predictions of activation or inhibition classified as strong, moderate, or weak that is based on a combination of the activation z-score and P-value of overlap. A list of differentially expressed proteins upon which the prediction was based (shown in red and green letter for up-regulated and down-regulated proteins in the high FE phenotype, respectively) is provided in Table 4.
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 8
A. Activated Upstream Regulators. 1. Insulin receptor (INSR), and insulin-like growth factor 1 receptor (IGF1R). Based on expression of downstream targets, INSR and IGF1R were predicted to be activated in the breast muscle from the high FE phenotype (Table 4). A thorough review of the insulin receptor signaling pathway in birds has presented differences that exist between mammals and birds [59]. After binding to its receptor, insulin activates INSR tyrosine kinase activity and phosphorylates downstream signaling molecules that includes the insulin receptor substrate (IRS) family leading to glucose homeostasis and energy storage [60– 64]. Insulin signaling is dependent upon endocytosis of the receptor that occurs in conjunction with caveolae and caveolin proteins (particularly caveolin 1, CAV-1) on the surface membrane of target cells [65]. In this regard, the 9.4 fold up-regulation of CAV-1[6] (Table 1), which is specific for insulin in mammals [65] would facilitate insulin receptor signaling in the high FE phenotype. Caveolae-mediated endocytosis signaling was the third ranked canonical pathway in the proteomics dataset (Table 3).
As in mammals, insulin increases glucose uptake in chicken muscle [66]. However, this is mediated by glucose transporter 8 (Glut-8) which is the avian analog to mammalian Glut-4 [67, 68]. A well characterized component of insulin receptor signaling is increased protein synthesis mediated by PI3K-Akt and mTOR (mammalian target of rapamycin) pathways (see Fig 27.3; in Dupont et al. [59]). Insulin-like growth factors (IGFs) along with PI3K/Akt signaling regulate muscle growth in various organisms (e.g. [69, 70]) including chickens [71]. Insulinlike growth factor 2 (IGF2) was up-regulated in high FE commercial broiler breast muscle [4] and in the latissimus dorsi (LD) muscle of castrated pigs [50] which supports a role of insulin-like growth factor in feed efficiency. Chen et al. [72] reported that mRNA expression of IGF binding protein 3 (IGFBP3) was higher in liver of more efficient low RFI angus bull calves (progeny of third generation of selection) and indicated that this lower expression likely contributed to lower fat deposition compared to high RFI bull calves. Expression of several genes of the PI3K signaling pathway, including PI3K, p85, pIK3R1, PP2A, MAP2K4, Ras, and Jnk, were up-regulated in the high FE breast muscle tissue [6] and led us to hypothesize that mTOR played a role in the phenotypic expression of FE [73]. Additionally, Lee et al. [3] reported that several genes in the mTOR pathway were up-regulated in liver and duodenum of broilers with low RFI (i.e. higher efficiency). The results of the present study along with Zhou et al. [4] and Lee et al. [3] indicate that insulin signaling and mTOR pathways may be biomarker candidates for FE selection. There is precedent for involvement of insulin signaling in FE in cattle and swine, particularly in young animals as recently reviewed by Davis et al. [74].
Downstream target proteins of INSR and IGF1R in Table 4 include five electron transport chain proteins; ATP50, COX4I1, NDUV2, UQCRC1 and UQCRC2. Insulin has been recognized to affect mitochondrial function early on in mitochondrial physiology investigations (e.g. [75–77]. Insulin signaling is integral to electron transport chain activity by maintaining the NAD+/NADH ratio that sustains signaling for mitochondrial biogenesis mediated by PGC1-α (see review, [78]). The PGC1-α protein had a moderate (1.531) activation z-score in high FE muscle (Table 4, PPARGC1α). AMP activated protein kinase (AMPK), which acts as an energetic sensor by detecting changes in the AMP/ATP ratio in the cell, stimulates mitochondrial energy production and mitochondrial biogenesis and is also tied to the insulin signaling pathway [78]. In this regard, the catalytic subunit of AMPK (PRKAA2, or AMPKα2) was predicted to be activated (1.951 activation z-score, Table 4) and mRNA expression of AMPK and two subunits (AMPKα and AMPKγ) were elevated in the high FE phenotype [6]. In addition, mTOR mRNA expression was elevated in PedM broiler males and in male quail selected for high FE [79].
In the present study, insulin degrading enzyme (IDE) is the third highest up-regulated protein in high FE breast muscle (11.7 fold, Table 1) and was also present in the 91 proteins deemed upregulated by the PIF metric described previously. This data concurs with up-regulated IDE mRNA expression (1.4 fold) in the high FE phenotype in the cDNA microarray data (See [5]) deposited in Gene Expression Omnibus (GEO; accession number: GSE24963; http:// www.ncbi.nlm.nih.gov/geo). IDE is involved in insulin signal transduction and IDE-mediated degradation of insulin has been linked to inhibition of protein degradation via a proteosomal mechanism [80–84]. Tundo et al. [85] reported that IDE also possesses heat-shock protein activity; thus IDE may also have protein chaperoning and scaffolding activities. Interestingly, Leissring et al. [86] reported a mitochondrial IDE isoform that could be instrumental in optimizing mitochondrial function. Mitochondrial proteins that are encoded by nuclear DNA typically contain a targeting sequence that must be cleaved upon entry into the mitochondria in order to be functional. Leissring et al. [86] provided evidence that IDE was capable of cleaving mitochondrial leader peptides that would play an important role in promoting optimal mitochondrial protein structure and subsequently support optimal mitochondrial function.
2. Nuclear factor, erythroid 2-like 2 (NFE2L2). Expression of the transcription factor NFE2L2 is increased in response to oxidative stress and enhances antioxidant protection by modifying genes with the antioxidant response element in their promotor regions (e.g. [87, 88]). From previous findings of higher mitochondrial ROS production and pervasive protein oxidation [8, 10], we anticipated that NFE2L2 would be activated in the low FE PedM. However, NFE2L2 was activated in the high FE PedM (Table 1), which concurs with findings in commercial broilers reported by Zhou et al. [4].
There are differences between the molecules that were used by the IPA program in the predicted activation of NFE2L2 between the present study using PedM broilers and that of Zhou et al. [4] conducted with commercial broilers. In the present study, of the downstream molecules used in predicting NFE2L2 activation (Table 4), only two antioxidant proteins were differentially expressed; glutathione peroxidase (GPX1, cytosolic) and superoxide dismutase (SOD1) that were up- and down-regulated in the high FE phenotype, respectively. As indicated above, PRDX3 was up-regulated in high FE muscle, but was not included in the downstream molecules used in predicting NFE2L2 activation. In the commercial male broiler FE model [4], up-regulated antioxidants in the high FE included 3 glutathione transferase enzymes, heme oxygenase 1, extracellular superoxide dismutase (SOD3), and thioredoxin (TXN). Whereas the ryanodine receptor (RYR3) was down-regulated in the high FE commercial broiler [4], RYR3 protein expression was up-regulated in high FE pedigree males in the present study. In addition, IDE, three proteosomal proteins (PSMA5, PSMB1, and PSMD1), ADP ribosylation factor 1 (ARF1), and NARS (asparaginyl-tRNA synthetase), (target molecules used in the prediction of NFE2L2 activation [Table 4]), did not appear in the predicted activation of NFE2L2 in high FE commercial broilers [4]. Higher expression of IDE could serve to regulate the insulin signaling and play a role in importing proteins into the mitochondria as discussed previously. Enhanced proteosomal activity in high FE seems counterintuitive for promoting greater efficiency, yet we observed higher mRNA expression of 19S proteasome in high FE breast muscle [73]. In addition, UBQLN1, another proteosomal protein, is down regulated in high FE (Fig 3) so the evidence for enhanced proteosomal activity in high FE requires further investigation.
Increased ARF1 expression in the high FE PedM broiler is of particular interest to us as this molecule was recently reported to play a critical role in mitochondrial homeostasis [89, 90]. Using inhibitory RNA to deplete ARF1, it was shown that muscle cells in C. elegans displayed altered mitochondrial morphology and diminished function [89]. The interaction of ARF1 with mitochondria enables the ‘...recruitment of a degradation machinery to mitochondria to remove toxic mitofusin/Fzo1 clusters; and it allows the extension of autophagy sequestration membranes needed for mitophagy to clear damaged mitochondria.” [90]. Thus, increased ARF1 expression could help in maintaining overall mitochondrial function in cells by removing mitochondria, or parts of mitochondria, that are under performing. Elevated ARF1 expression in high FE PedM broiler muscle represents an additional link between mitochondria and FE. Interestingly, cluster analysis showed ARF4 to be highly co-expressed with ARF1 across the 8 samples (Fig 3), implying the possibility of shared function between the two proteins. Indeed, Nakai et al. [91] found ARF1 and ARF4 are collectively required for the integrity of recycling endosomes, with endosomes playing a role in lysosomal degradation.
Zhou et al. [4] indicated that the oxidative stress in the high FE commercial broiler was due to muscle development; i.e. that higher NFE2L2 activity was in response to inflammatory-mediated oxidative stress as there was up-regulation of a number of myokines and cytokines in their RNAseq dataset. This contrasts with the evidence of lower oxidative stress in high FE PedM broilers [10]. Additional evidence of higher stress in the low FE pedigree male broiler included the up-regulation of corticotropin releasing hormone, actin polymerization and actin-myosin stress fiber formation, heat shock proteins, and suppressor of cytokine signaling (SOC3) [5–7]. As body weights were similar between the report by Zhou et al. [4] and the present study, differences in oxidative stress cannot be explained by differences in growth rate between the studies. This leads us to hypothesize that the predicted activation of NFE2L2 in the high FE PedM broiler may play a role in directing cellular processes to occur, rather than responding to oxidative stress as indicated by Zhou et al. [4] in commercial broilers.
3. Progesterone and Triiodothyronine. Progesterone and triiodothyronine (T3) were predicted to be activated (z-activation scores of 1.954 and 1.951, respectively) in breast muscle of the high FE PedM broiler phenotype (Table 4). Implants containing progesterone and estrogen have been used for years to improve growth and FE in the beef industry (see review by Preston, [92]). These implants mediate their effects by stimulating growth hormone activity and the insulin/insulin- like growth factor signaling pathway (e.g. [93–95]). Thus, activation of progesterone could contribute to insulin receptor activation described above. In a separate but concurrent RNA sequencing study on the same breast muscle tissue samples, we have additional indication that progesterone is likely to be centrally involved in activation of several transcription factors in the high FE phenotype (unpublished observations).
The predicted activation of triiodothyronine (T3, Table 4) follows the increase in the mitoproteome in high FE but is surprising since it has long been associated with increased energy expenditure through increased basal metabolic rate. However, a report by Clement et al. [96] investigating the impact of thyroid hormone on human skeletal muscle in vivo has several remarkable parallels to the present proteomic investigation. Clement et al. [96] reported that expression of 15 genes encoding mitochondrial electron transport chain proteins were up-regulated by thyroid hormone (75 μg/day for 14 d) similar to the up-regulation of electron transport proteins shown in high FE (Fig 1). From both in vivo and in vitro studies, Enriquez et al. (1999) reported that mitochondrial RNA synthesis was directly stimulated by thyroid hormone. Based on findings of Enriquez et al. [97] and others, Psarra and Sekeris [58] indicated that the mitochondrial thyroid hormone receptor “...is a mitochondrial, thyroid hormone activated transcription factor which in analogy to the action of the nuclear transcription factors, would exert its modulatory effect by interaction with the mitochondrial transcriptome.” Clement et al. [96] also reported that several proteasome and ubiquitin-related genes were up-regulated, which concurs with up-regulated 19S proteosome mRNA expression in the high FE phenotype [73] and proteosomal proteins in the present study (see target molecules for NFE2L2, Table 4). Recently, Piekarski et al. [79] reported that gene expression of key members of the autophagy pathway were higher in the high FE phenotype, which agrees with the up-regulation of Beclin associated with thyroid hormone reported by [96]. These anecdotal similarities do not prove that thyroid hormone is involved in the phenotypic expression of FE, but do warrant further investigation.
Cystathionine gamma lyase (CTH) is also a downstream target molecule of T3 (Table 4). As this enzyme is involved in the transsulfuration pathway in the conversion of methionine to cysteine, it conceivably could play a role in maintaining glutathione, which is the major intracellular antioxidant, as cysteine is the rate limiting amino acid of glutathione synthesis [98].
B. Inhibited Upstream Regulators. 1. Rapamycin independent companion of target of rapamycin (RICTOR). We previously reported that there was decreased mRNA expression of numerous actin-myosin cytoskeletal fibers in the high FE phenotype [5, 6]. RICTOR plays an important role in actin-cytoskeletal formation as impaired formation was observed in RICTOR knockout mice [99, 100]. Thus, inhibition of RICTOR in the high FE PedM broiler phenotype could play a role in the lower cytoskeletal fiber gene expression that we reported previously [5, 6]. Inhibition of serum response factor (SRF) (see below), could also contribute to lower cytoskeletal organization observed in the high FE PedM broiler. [101] reported that the RICTOR component of mTORC2 is integral in actin-stress fiber development that stimulates bone marrow derived mesenchymal stem cells to become osteoblasts. Thus, the observation that actin stress fiber formation was down-regulated in the high FE phenotype [6] could be directly related to lower RICTOR activity. Unlike ablation of regulatory associated protein of mTOR, Complex 1 (RAPTOR) activity of mTORC1, which is critical for muscle function, decreased RICTOR activity has no effect on muscle function [102].
The predicted inhibition of RICTOR, however, is incongruent with insulin signaling described above. Insulin stimulation of glucose uptake in muscle of mammals is carried out by mTORC2-mediated facilitation of glucose transport 4 (Glut-4) receptors [103]. Kumar et al. [104] reported that there was a decrease in insulin-mediated glucose uptake and an increase in basal glycogen synthase activity in male adult mTORC2 knockout mice. To our knowledge, effects of mTORC2 on the presumed glucose transporter protein (Glut-8) and glycogen synthase activity has not been reported in avian species.
2. Mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4). In the current study, MAP4K4 was predicted to be inhibited in the high FE PedM broiler phenotype (Table 4). Guntur et al. [105] reported that MAP4K4 acts as a negative regulator of PPARγ-mediated stimulation of translation in adipocytes by suppressing the mTOR signaling pathway. More recently, Wang et al. [106] demonstrated that myotube formation was enhanced in MAP4K4 mutant C2C12 cells via a myogenic factor 5 (Myf5)-dependent mechanism and by other mitogen-activated protein kinase (MAPK)-signaling pathways leading the researchers to conclude that MAP4K4 was a novel suppressor of skeletal muscle differentiation [106]. Thus, inhibition of MAP4K4 activity could contribute to muscle development in the high FE PedM broiler muscle in the present study.
3. Serum Response Factor (SRF). Serum response factor (SRF) was predicted to be inhibited in the high FE phenotype muscle (Table 4). SRF signaling is stimulated by MAPK or RhoA pathways that enhance nuclear transcription primarily of cytoarchitecture-related genes [107]. SRF was reported to bind to the promoter region of several cytoskeletal-contractile genes in fully differentiated muscle tissue [108–110]. Miano et al. [111] reviewed the literature on SRF and presented a list of ~120 cytoskeletal-contractile genes stimulated by SRF-mediated transcription (see Table 2 on p. C73 [111]). We hypothesized that the higher degree of cytoskeletal organization would contribute to the phenotypic expression of low FE due to the increased energy expenditure needed to maintain this cellular architecture [5, 6]. From the GEO database that was deposited with this microarray set, expression of 27 genes (approximately 25% of the SRF target genes listed in Table 1 of [111]) were DE in the cDNA microarray dataset; of these 27 genes, 25 were down-regulated in the high FE phenotype (from [5]) which concurs with the predicted inhibition of SRF (and its effect on cytoarchitecture gene transcription) in the current proteomic dataset. Also, SRF mRNA expression was down-regulated 1.25 fold in the high FE phenotype (unreported observation, [5]).
III. Summary and Conclusions
The diagram shown in Fig 5 presents a summary of the results that are presented in this shotgun proteomic study conducted in a pedigree male (PedM) broiler line exhibiting high or low FE phenotypes. The results of previous studies conducted using this PedM broiler line that indicated that the high FE phenotype had better respiratory chain coupling and increased respiratory chain complex activity [8, 10], appears to be recapitulated in the proteomic analysis indicating that 4 of 5 electron transport chain complexes were predicted to be activated in the high FE pedigree broiler male (Fig 4B) and the entire mitoproteome is skewed towards the high FE birds (Fig 2A). Of considerable interest are the predictions that both progesterone and thyroid hormone would be activated in the high FE phenotype. It is well known that exogenous treatment with progesterone (as one part of a combination hormone growth promotant treatment) plays a role in improving FE in cattle that is mediated by the insulin signaling pathway. With the exception of NFE2L2 and SRF, all other upstream regulators had at least one mitochondrial electron transport protein as a target molecule that was used in their respective predictions (see Table 4). In particular, the increase in mitochondrial protein expression would be supported by moderate activation predicted for two transcription factors, PPARGC1α and PPARα, associated with mitochondrial biogenesis. With the lower oxidative stress that we have consistently observed in the high FE phenotype (e.g. [10]), we hypothesize that the predicted activation of NFE2L2 (which is typically increased in response to oxidative stress), reflects an inherent difference that helps orchestrate gene expression in the high FE phenotype. The predicted activation of SRF might be important in enhanced cytoarchitectural gene expression reported in the low FE PedM broiler phenotype [5, 8]. We believe that further analysis of the proteomic and genomic datasets that we have obtained in the high and low FE PedM broilers will generate hypotheses that can then be mechanistically tested both in vivo and in vitro to yield further insights into fundamental understanding of the cellular basis of feed efficiency
Proteomics of Breast Muscle Tissue Associated with the Phenotypic Expression of Feed Efficiency within a Pedigree Male Broiler Line: I. Highlight on Mitochondria - Image 9
The broad findings of the current study indicate that upstream analysis implicates involvement of progesterone and triiodothyronine that would increase expression of mitochondrial proteins as indicated in Fig 2. This would be complemented by PGC1α and PPARA that would stimulate mitochondrial biogenesis. The predicted activation of NFE2L2 would enhance proteosomal activity (in agreement with Bottje et al. [2014]) and could enhance mechanisms associated with protein refolding. NFE2L2 activation was also attributed to increased expression of several proteins (PPIB, EIF4G2, and PPIB) as well as decreased inhibition exerted by RICTOR that collectively could result in increased protein synthesis. NFE2L2 predicted activation was also based on increased expression of ARF1 (ADP ribosylation factor 1) that may help in maintaining mitochondrial homeostasis [89, 90]. Increased expression of ANT and VDAC in the inner and outer mitochondrial membranes would be instrumental in exchange of ATP and ADP to and from the cytosol.
Supporting Information
S1 Fig. Fold expression values of house-keeping proteins were similar for Low (green bar) and high FE (red bar) resulting from normalization procedures (see text for details).Proteins include: AGTN2 (alpha-actin 2), MYH1G (myosin heavy chain 1 gamma), MYH7B (myosin heavy chain 7 beta), TBB5 (tubulin beta 5 chain), TBB7 (tubulin beta 7 chain), TBA1 (tubulin alpha 1 chain). Values represent the mean ± SE (n = 4). (TIF)
S2 Fig. Polymerase chain reaction (PCR) results for determining mitochondrial (mt) DNA content in breast muscle obtained from Pedigree Male Broilers exhibiting for low or high feed efficiency phenotypes (LFE and HFE). Values represent mean ± SE (n = 4). Primer information: mtDNA-D loop Forward—ACACCTGCGTTGCGTCCTA; Reverse—ACGCAAACCGTCTCA TCGA; 18S rRNA gene Forward—TCCCCTCCCGTTACTTGGAT; Reverse—GCGCTCGTCGG CATGTA. (TIF)
The authors wish to thank Cobb-Vantress, Inc. (Siloam Springs, AR) for providing access to the pedigree male broiler breeder line that was investigated in this study and for prior research support in the area of mitochondrial function and feed efficiency.
Author Contributions
Conceived and designed the experiments: WB BK NH. Performed the experiments: WB BK NH. Analyzed the data: WB BK NH AR-G. Contributed reagents/materials/analysis tools: BK WB NH AR-G. Wrote the paper: WB BK NH AR-G KL AP SD. Critical review of manuscript: AR-G.
This article was originally published in PLoS ONE 11(5): e0155679. doi:10.1371/journal.pone.0155679. This is an Open Access article distributed under the terms of the Creative Commons Attribution License.

1. Simmons J. (2010) Technology’s role in the 21st century—Elanco Animal Health. Why agriculture needs technology to help meet a growing demand for safe; nutritious and affordable food.

2. Nixon R, Lowry A. (2012) Drought forces reductions in U.S. crop forecasts. New York Times Business. Available: http://www.nytimes.com/2012/08/11/business/projections-for-corn-yield-falls-to-17-yearlow.html.

3. Lee J, Karnuah AB, Rekaya R, Anthony NA, Aggrey SE (2015) Transcriptomic analysis to elucidate the molecular mechanisms that underlie feed efficiency in meat-type chickens. Molecular Genetics and Genomics doi: 10.1007/s00438-015-1025-7

4. Zhou N, Lee WR, Abasht B (2015) Messenger RNA sequencing and pathway analysis provide novel insights into the biological basis of chickens' feed efficiency. BMC Genomics 16: 195–215. doi: 10. 1186/s12864-015-1364-0 PMID: 25886891

5. Kong B, Song J, Lee JY, Hargis BM, Wing T, Lassiter K et al. (2011) Gene expression in breast muscle associated feed efficiency in a single male broiler line using a chicken 44k microarray. I. top differentially expressed genes. Poultry Science 90: 2535–2547. doi: 10.3382/ps.2011-01435 PMID: 22010239

6. Bottje WG, Kong BW, Song JJ, Lee JY, Hargis BM, Lassiter K et al. (2012) Gene expression in breast muscle associated feed efficiency in a single male broiler line using a chicken 44k microarray II. Differentially expressed focus genes. Poultry Science 91: 2576–25887. PMID: 22991544

7. Bottje W, Kong BW (2013) Feed efficiency: Mitochondrial function to global gene expression. Journal of Animal Science 91: 1582–1593. doi: 10.2527/jas.2012-5787 PMID: 23148240

8. Bottje WG, Iqbal M, Tang Z, Cawthon DC, Okimoto R, Wing T, et al. (2002) Association of mitochondrial function with feed efficiency within a single genetic line of male broilers. Poultry Science 81: 546–555. PMID: 11989755

9. Ojano-Dirain C, Bottje W, Wing T, Cooper M (2005) Glutathione and respiratory chain complex activities in duodenal mitochondria from broilers with low and high feed efficiency. Poultry Science 84: 782–788. PMID: 15913191

10. Bottje W, Carstens GE (2009) Association of mitochondrial function and feed efficiency in poultry and livestock species. Journal of Animal Science 87: E48–E63. doi: 10.2527/jas.2008-1379 PMID: 19028862

11. Taniguchi Y, Choi PJ, Li GW, Chen H, Babu M, Hearm J, et al. (2010) Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science 329: 533–538. doi: 10. 1126/science.1188308 PMID: 20671182

12. Greenbaum D, Colangelo C, Williams K, Gerstein M (2003) Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biology 4: 117. PMID: 12952525

13. Pascal LE, True LD, Campbell DS, Deutsch EW, Risk M, Colman IM, et al. (2008) Correlation of mRNA and protein levels: Cell type-specific gene expression of cluster designation antigens in the prostate. BMC Genomics 9: 246. doi: 10.1186/1471-2164-9-246 PMID: 18501003

14. Guo Y, Xiao P, Lei S, Deng F, Xiao GG, Liu Y, et al. (2008) How is mRNA expression predictive for protein expression? A correlation study on human circulating monocytes. Acta Biochimica Biophysica Sin (Shanghai) 40: 426–436.

15. Nie L, Wu G, Zhang W (2006) Correlation of mRNA expression and protein abundance affected by multiple sequence features related to translational efficiency in desulfovibrio vulgaris: A quantitative analysis. Genetics 174: 2229–2243. PMID: 17028312

16. Zanetti E, Molette C, Chambon C, Pinguet J, Remignon H, Cassandro M (2011) Using 2-DE for the differentiation of local chicken breeds. Proteomics 2613: 2613–2619.

17. Sokale A, Peebles ED, Zhai W, Pendarvis K, Burgess S, Perchan T (2011) Proteome profile of the pipping muscle in broiler embryos. Proteomics 11: 4262–4265. doi: 10.1002/pmic.201000795 PMID: 21834139

18. Kaweewong K, Garnjanagoonchorn W, Jirapakkul W, Roytrakul S (2013) Solubilization and identification of hen eggshell membrane proteins during different times of chicken embryo development using the proteomic approach. Protein Journal 32: 297–308. doi: 10.1007/s10930-013-9487-0 PMID: 23636516

19. Liu Y, Qiu N, Ma M (2013) Comparative proteomic analysis of hen egg white proteins during early phase of embryonic development by combinatorial peptide ligand library and matrix-assisted laser desorption ionization-time of flight. Poultry Science 92: 1897–1904. doi: 10.3382/ps.2012-02986 PMID: 23776278

20. Sun C, Xu G, Yang N (2013) Differential label-free quantitative proteomic analysis of avian eggshell matrix and uterine fluid proteins associated with eggshell mechanical property. Proteomics 13: 3523– 3536. doi: 10.1002/pmic.201300286 PMID: 24151251

21. Mu P, Xu M, Zhang L, Wu K, Wu J, Chen O (2013) Proteomic changes in chicken primary hepatocytes exposed to T-2 toxin are associated with oxidative stress and mitochondrial enhancement. Proteomics 13: 3175–3188. doi: 10.1002/pmic.201300015 PMID: 24030954

22. Chien KY, Blackburn K, Liu HC, Goshe MB (2012) Proteomic and phosphoproteomic analysis of chicken embryo fibroblasts infected with cell culture-attenuated and vaccine strains of Marek's disease virus. Journal of Proteome Research 11: 5663–5677. doi: 10.1021/pr300471y PMID: 23106611

23. Fan Z, Hu X, Zhang Y, Yu C, Qian K, et al. (2012) Proteomics of DF-1 cells infected with avian leukosis virus subgroup. Journal of Virus Research 167: 314–321. doi: 10.1016/j.virusres.2012.05.016 PMID: 22664373

24. Haq K, Brisbin JT, Thanthrige-Don N, Heidari M, Sharif S (2010) Transcriptome and proteome profiling of host responses to Marek's disease virus in chickens. Veterinary Immunology and Immunopathology 138: 292–302. doi: 10.1016/j.vetimm.2010.10.007 PMID: 21067815

25. Hu X, Qin A, Qian K, Shao H, Yu C, Mio J (2012) Analysis of protein expression profiles in the thymus of chickens infected with Marek's disease virus. Virology Journal 9: 256. doi: 10.1186/1743-422X-9- 256 PMID: 23116199

26. Wu Y, Peng C, Xu L, Zheng X, Liao M, Jin Y (2012) Proteome dynamics in primary target organ of infectious bursal disease virus. Proteomics 12: 1844–1859. doi: 10.1002/pmic.201100479 PMID: 22623289

27. Picard B, Berri C, Lefaucheur L, Molette C, Sayd T, Terlouw C (2010) Skeletal muscle proteomics in livestock production. Brief Functional Genomics 9: 259–278.

28. Rabilloud T, Chevallet M, Luche S, Lelong C. (2010) Two-dimensional gel electrophoresis in proteomics: Past, present and future. Journal of Proteomics 73: 2064–2077. doi: 10.1016/j.jprot.2010.05. 016 PMID: 20685252

29. Mann M, Kulak NA, Nagaraj N, Cox J. (2013) The coming age of complete, accurate, and ubiquitous proteomes. Molecular Cell 49: 583–590. doi: 10.1016/j.molcel.2013.01.029 PMID: 23438854

30. Altelaar AF, Munoz J, Heck AJ. (2013) Next-generation proteomics: Towards an integrative view of proteome dynamics. Nat Rev Genet 14: 35–48. doi: 10.1038/nrg3356 PMID: 23207911

31. Iqbal M, Pumford NR, Tang Z, Lassiter K, Wing T, Cooper M, et al. (2004) Low feed efficient broilers within a single genetic line exhibit higher oxidative stress and protein expression in breast muscle with lower mitochondrial complex activity. Poultry Science 83: 474–484. PMID: 15049502

32. Hudson NJ, Reverter A, Dalrymple BP (2009) A differential wiring analysis of expression data correctly identifies the gene containing the causal mutation. PLoS Computational Biology 5: 1–15.

33. Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z (2009) GOrilla: A tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics 10: 1471–2105.

34. Caraux G, Pinloche S (2005) PermutMatrix: A graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics 21: 1280–1281. PMID: 15546938

35. Smyth GK. (2004) Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistics of Applied Genetics and Molecular Biology 3: 1–24.

36. Vincent A, Louveau I, Gondret G, Trefeu C, Gilbert H, Lefaucheur L (2015) Divergent selection for residual feed intake affects the transcriptomic and proteomic profiles of pig skeletal muscle. Journal of Animal Science 93: 2745–2758. doi: 10.2527/jas.2015-8928 PMID: 26115262

37. Kiessling KH (1977) Muscle structure and function in the goose, quail, pheasant, guinea hen, and chicken. Comparative Biochemistry and Physiology B 57: 287–292.

38. Kim H, You S, Kim I-J, Farris J, Foster LK, Foster DN (2001) Increased mitochondrial gene transcription in immortal DF-1 cells. Experimental Cellular Research 265: 339–347.

39. Sharifabadi HR, Zamir MJ, Rowghani E, Bottje WG (2012) Relationship between the activity of mitochondrial respiratory chain complexes and feed efficiency in fat-tailed ghezel lambs. Journal of Animal Science 90: 1807–1815. doi: 10.2527/jas.2011-4791 PMID: 22147474

40. Lemasters JJ, Holmuhamedov E (2006) Voltage-dependent anion channel (VDAC) as mitochondrial governator—thinking outside the box. Biochemica Biophysica Acta 1762: 181–190.

41. Galganska J, Karachitos A, Wojtkowska M, Stoienia O, Budzinska M, Kimta H (2010) Communication between mitochondria and nucleus: Putative role for VDAC in reduction/oxidation mechanism. Biochemica Biophysica Acta 1797: 1276–1280.

42. Han D, Antunes F, Canali R, Retton D, Cadenas E (2003) Voltage-dependent anion channels control the release of the superoxide anion from mitochondria to cytosol. Journal of Biological Chemistry 278: 5557–5563. PMID: 12482755

43. Schlattner U, Tokarska-Schlattner M, Wallimann T (2006) Mitochondrial creatinine kinase in health and disease. Biochemica Biophysica Acta 1762: 164–180.

44. Walliman T, Wyss SM, Brdiczka D, Nicolay K, Eppenberger HM (1992) Intracellular compartmentation, structure and function of creatine kinase isoenzymes in tissues with high and fluctuating energy demands: The 'phosphocreatine circuit' for cellular energy homeostasis. Biochemistry Journal 281: 21–40.

45. Quest AFG, Eppenberger HM, Walliman T (1990) Volume 262, number 2, 299–304 FEBS 08264 march 1990 two different B-type creatine kinase subunits dimerize in a tissue-specific. FEBS Letters 262: 299–204.

46. Watabe S, Hiro T, Yamamoto Y, Fugioka Y, Hasegawa H, Yago N (1997) SP-22 is a thioredoxindependent peroxide reductase in mitochondria. European Journal of Biochemistry 249: 52–60. PMID: 9363753

47. Chang T, Cho C, Park C, Yu S, Kang SW, Rhee SG (2004) Peroxiredoxin III, a mitochondrion-specific peroxidase, regulates apoptotic signaling by mitochondria. Journal of Biological Chemistry 279: 41975–41984. PMID: 15280382

48. Noh YH, Baek JY, Jeong W, Rhee SG, Chang T (2009) Sulfiredoxin translocation into mitochondria plays a crucial role in reducing peroxidized peroxiredoxin III. Journal of Biological Chemistry 284: 8470–8477. doi: 10.1074/jbc.M808981200 PMID: 19176523

49. Wonsey DR, Zeller KI, Dang CV (2002) The c-myc target gene PRDX3 is required for mitochondrial homeostasis and neoplastic transformation. Proceedings of the National Academy of Science (USA) 99: 6649–6654.

50. Jing L, Hou Y, Miao Y, Li X, Cao J, Brameld JM et al. (2015) Transcriptome analysis of mRNA and miRNA in skeletal muscle indicates an important network for differential residual feed intake in pigs. Scientific Reports 5:11953: doi: 10.1038/srep11953 PMID: 26150313

51. Nisoli E, Clementi E, Paolucci C, Cozzi V, Tonello C, Sciorati C, et al. (2003) Mitochondrial biogenesis in mammals: The role of endogenous nitric oxide. Science 299: 896–899. PMID: 12574632

52. Nisoli E, Clementi E, Moncada S, Carruba MO (2004) Mitochondrial biogenesis as a cellular signaling framework. Biochemical Pharmacology 67: 1–15. PMID: 14667924

53. Irrcher I, Adhihetty PJ, Sheehan T, Joseph A, Hood DA (2003) PPARγ coactivator-1α expression during thyroid hormone- and contractile activity-induced mitochondrial adaptations. American Journal of Physiology Cell Physiology 284: C1669–C1677. PMID: 12734114

54. Zhou G, Myers R., Li Y., Chen Y., Shen X., Fenyk-Melody J, et al. (2001) Role of AMP-activated protein kinase in mechanism of metformin action. Journal of Clinical Investigation 108: 1167–1174. PMID: 11602624

55. Hardie DG, Scott JW, Pan DA, Hudson ER (2003) Management of cellular energy by the AMP-activated protein kinase system. FEBS Letters 546: 113–120. PMID: 12829246

56. Hardie DG (2005) New roles for the LKB1-AMPK pathway. Current Opinions in Cell Biology 17: 167– 173.

57. Carling D (2004) The AMP-activated protein kinase cascade—A unifying system for energy control. Trends in Biochemical Science 29: 18–24.

58. Psarra AG, Sekeris CE (2008) Steroid and thyroid hormone receptors in mitochondria. Life 60: 210– 223. doi: 10.1002/iub.37 PMID: 18344181

59. Dupont J, Rideau N, Simon J (2015) Endocrine pancreas. In: Scanes CG, editor. Sturkie's Avian Physiology ( 6 edition). New York: Elsevier. pp. 613–631.

60. Eldar-Finkelman H, Schreyer SA, Shinohara MM, Leboeuf RC, Krebs EG (1999) Increased glycogen synthase kinase-3 activity in diabetes- and obesity-prone C57BL/6J mice. Diabetes 48: 1662–1666. PMID: 10426388

61. White MF (2002) IRS proteins and the common path to diabetes. American Journal of Physiology Endocrinology and Metabolism 283: E413–E422. PMID: 12169433

62. Baumann CA, Brady MJ, Saltiel AR (2001) Activation of glycogen synthase by insulin in 3T3-L1 adipocytes involves c-cbl-associating protein (CAP)-dependent and CAP-independent signaling pathways. Journal of Biological Chemistry 276: 6065–6068. PMID: 11222622

63. Baumann CA, Ribon C, Kanzaki M, Thurmond DC, Mora S, Bickle PE, et al. (2000) CAP defines a second signaling pathway required for insulin-stimulated glucose transport. Nature 407: 202–207. PMID: 11001060

64. Khan AH, Pessin JE (2002) Insulin regulation of glucose uptake: A complex interplay of intracellular signaling pathways. Diabetologia 45: 1475–1483. PMID: 12436329

65. Cohen AW, Combs TP, Scherer PE, Lisanti MP (2003) Role of caveolin and caveolae in insulin signaling and diabetes. American Journal of Physiology Endocrinology and Metabolism 285: E1151– E1160. PMID: 14607781

66. Tokushima R, Takahaski K, Sato K, Akiba Y (2005) Glucose uptake in vivo in skeletal muscles of insulin-injected chicks. Comparative Biochemistry and Physiology B 141: 43–48.

67. Kono T, Nishida M, Nishiki Y, Seki Y, Sato K, Akiba Y (2005) Characterization of glucose transporter (GLUT) gene expression in broiler chickens. British Poultry Science 46: 510–516. PMID: 16268111

68. Zhao JP, Bao J, Wang XJ, Jiano HC, Song ZG, Lin J (2012) Altered gene and protein expression of glucose transporter 1 underlies dexamethazone inhibition of insulin-stimulated glucose uptake in chicken muscles. Journal of Animal Science 90: 4337–4345. doi: 10.2527/jas.2012-5100 PMID: 22859751

69. Schiaffino S, Mammucari C (2011) Regulation of skeletal muscle growth by the IGF1-Akt/PKB pathway: Insights from genetic models. Skeletal Muscle 1: 1–4.

70. Fuentes EN, Björnsson BT, Valdés JA, Einarsdottir IE, Lorca B, Alvarez M (2011) IGF-I/PI3K/Akt and IGF-I/MAPK/ERK pathways in vivo in skeletal muscle are regulated by nutrition and contribute to somatic growth in the fine flounder. (2011) American Journal of Physiology Regulation and Integration of Comparative Physiology 300: R1532–R1542.

71. Duclos MJ (2005) Insulin-like growth factor-I (IGF-1) mRNA levels and chicken muscle growth. J of Physiology and Pharmacology 56 (Suppl 3): 25–35.

72. Chen Y, Gondro C, Quinn K, Herd RM, Parnell PF, Vaneslow B (2011) Global gene expression profiling reveals genes expressed differentially in cattle with high and low residual feed intake. Animal Genetics 42: 475–490. doi: 10.1111/j.1365-2052.2011.02182.x PMID: 21906099

73. Bottje WG, Kong B, Lee JY, Washington T, Baum JI, Dridi S, et al. (2014) Potential roles of mTOR and protein degradation pathways in the phenotypic expression of feed efficiency in broilers. Journal of Biochemical Physiology 3: 1–8.

74. Davis ME, Wick MP, Maquivar MG (2012) Hormonal regulation of feed efficiency. In: Hill R, editor. Feed efficiency in the beef industry. Ames, Iowa: Wiley-Blackwell. pp. 175–198.

75. Krebs EG, Eggleston LV (1938) The effect of insulin on oxidations in isolated muscle tissue. Biochemical Journal 32: 913–925. PMID: 16746703

76. Shorr E, Barker SB (1939) In vitro action of insulin on minced avian and mammalian muscle. Biochemical Journal 33: 1798–1809. PMID: 16747098

77. Hall JC, Sordahl LA, Stepko PL (1960) The effect of insulin on oxidative phosphorylation in normal and diabetic mitochondria. Journal of Biological Chemistry 235: 1536–1539. PMID: 14399044

78. Cheng Z, Tseng Y, White MF (2010) Insulin signaling meets mitochondria in metabolism. Trends in Endocrinology Metabolism 21: 589–598. doi: 10.1016/j.tem.2010.06.005 PMID: 20638297

79. Piekarski A (2015) Autophagy and its potential role in stress and feed efficiency in avian tissues. Ph. D. Dissertation, University of Arkansas, Fayetteville.

80. Fawcett J (2014) Biologically active insulin-derived peptides. Protein and Peptide Letters 21: 584– 588. PMID: 24559166

81. Smith OL, Wong CY, Gelfand RA (1989) Skeletal muscle proteolysis in rats with acute streptozocininduced diabetes. Diabetes 38: 1117–1122. PMID: 2670642

82. Price SR, Bailey JL, Wang X, Jurkovitz C, England BK, Ding X, et al. (1996) Muscle wasting in insulinopenic rats results from activation of the ATP-dependent, ubiquitin-proteasome proteolytic pathway by a mechanism including gene transcription. Clin Invest 98: 1703–1708.

83. Merforth S, Osmers A, Dahlmann B (1999) Alterations of proteasome activities in skeletal muscle tissue of diabetic rats. Molecular Biology of Reproduction 26 83–87.

84. Galban VD, Evangelista EA, Migliorini RH, do Carmo Kettelhut I (2001) Role of ubiquitin–proteasomedependent proteolytic process in degradation of muscle protein from diabetic rabbits. Molecular and Cellular Biochemistry 225: 35–41. PMID: 11716362

85. Tundo GR, Sbardella D, Ciacio C, Bianculli A, Orlandi A, Desimio MG, et al. (2013) Insulin-degrading enzyme (IDE): A novel heat shock-like protein. Journal of Biological Chemistry 288: 2281–2289. doi: 10.1074/jbc.M112.393108 PMID: 23188819

86. Leissring MA, Farris W, Wu X, Christodoulou DC, Haigis MC, Guarente L, et al. (2004) Alternative translation initiation generates a novel isoform of insulin-degrading enzyme targeted to mitochondria. Biochemical Journal 383: 439–446. PMID: 15285718

87. Kobayashi A, Kang M, Okawa H, Ohtsuji M, Zenke Y, Chiba T, et al. (2004) Oxidative stress sensor Keap1 functions as an adaptor for Cul3-based E3 ligase to regulate proteosomal degradation of Nrf2. Molecular and Cellular Biology 24(16): 7130–7139. PMID: 15282312

88. Kobayashi A, Kang M, Watai Y, Tong KI, Shibata T, Koli U, et al. (2006) Oxidative and electrophilic stresses activate Nrf2 through inhibition of ubiquitination activity of Keap1. Molecular and Cellular Biology 26: 221–229. PMID: 16354693

89. Ackema KB, Hench J, Böckler S, Wang SC, Sauder S, Mergantaler H, et al. (2014) The small GTPase Arf1 modulates mitochondrial morphology and function. EMBO Journal 33: 2659–2675. doi: 10. 15252/embj.201489039 PMID: 25190516

90. Rabouille C (2014) Old dog, new tricks: Arf1 required for mitochondria homeostasis. EMBO Journal 33: 2604–2605. doi: 10.15252/embj.201489899 PMID: 25230932

91. Nakai W, Kondo Y, Saitoh A, Naito T, Nakayama K, Shin HW (2013) ARF1 and ARF4 regulate recycling endosomal morphology and retrograde transport from endosomes to the golgi apparatus. Molecular and Cellular Biology 16: 2570–2581.

92. Preston RL (1999) Hormone containing growth promoting implants in farmed livestock. Advances in Drug Delivery Review 38: 123–138.

93. Johnson BJ, Hathaway MR, Anderson PT, Meiske JC, Dayton WR (1996) Stimulation of circulating insulin-like growth of feedlot steers due to administration of a combined trenbolone acetate and estradiol implant in feedlot cattle. Journal of Animal Science 74: 372–379. PMID: 8690673

94. Johnson BJ, White ME, Hathaway MR, Christians CJ (1998) Effect of a combined trenbolone acetate and estradiol (E) and recombinant bovine somatotropin (bST) 2 estradiol implant on steady-state IGFI mRNA concerning the anabolic response of feedlot steers. Journal of Animal Science 76: 491–497. PMID: 9498357

95. Thompson SH, Boxhorn LK, Allen RE (1989) Trenbolone alters the responsiveness of skeletal muscle trenbolone acetate to fibroblast growth factor and insulin-like growth factor 1. Endocrinology 124: 2110–2117. PMID: 2707149

96. Clement K, Viguerie N, Diehn M, Alizadeh PB, Thalamas C, Storey JD, et al. (2002) In vivo regulation of human skeletal muscle gene expression by thyroid hormone. Genome Research 12: 281–291. PMID: 11827947

97. Enriquez JA, Fernandez P, Garrido-Perez N, Lopez-Perez MJ, Perez-Martos A, Montova J (1999) Direct regulation of mitochondrial RNA synthesis by thyroid hormone. Molecular and Cellular Biology 19: 657–670. PMID: 9858589

98. Stipanuk MH (2004) Sulphur amino acid metabolism: Pathways for production and removal of homocysteine and cysteine. Annual Review of Nutrition 24: 539–577. PMID: 15189131

99. Jacinto E, Loewith A, Schmidt A, Lin S, Ruegg MA, Hall A, et al. (2004) Mammalian TOR complex 2 controls the actin cytoskeleton and is rapamycin insensitive. Nature Cell Biology 6: 1122–1128. PMID: 15467718

100. Sarbassov DD, Ali SM, Kim DH, Guertin DA, Latek RR, Erdoiment-Bromage H, et al. (2004) Rictor, a novel binding partner of mTOR, defines a rapamycin-insensitive and raptor-independent pathway that regulates the cytoskeleton. Current Biology 14: 1296–1302. PMID: 15268862

101. Sen B, Xie Z, Case N, Thompson WR, Uzer G, Styner M, et al. (2014) mTORC2 regulates mechanically induced cytoskeletal reorganization and lineage selection in marrow derived mesenchymal stem cells. Journal of Bone and Mineral Research 29: 78–89. doi: 10.1002/jbmr.2031 PMID: 23821483

102. Bentzinger CF, Romanino K, Cloetta D, Lin S, Mascarenhjas JB, Oliveri F, et al. (2008) Skeletal muscle-specific ablation of raptor, but not rictor, causes metabolic changes and results in muscle dystrophy. Cell Metabolism 8: 411–424. doi: 10.1016/j.cmet.2008.10.002 PMID: 19046572

103. Kim YB, Peroni OD, Minokoshi Y, Aschenbach WG, Minokoshi Y, Kotani A, et al. (2005) Muscle-specific deletion of Glut4 glucose transporter alters multiple regulatory steps in glycogen metabolism. Molecular and Cellular Biology 25: 9713–9723. PMID: 16227617

104. Kumar A, Harris TE, Keller SR, Choi KM, Magnuson MA, Lawrence JC (2008) Muscle-specific deletion of rictor impairs insulin-stimulated glucose transport and enhances basal glycogen synthase activity. Molecular and Cellular Biology 28: 61–70. PMID: 17967879

105. Guntur KVP, Guilherme A, Chawla A, Czech MP (2010) Map4k4 negatively regulates peroxisome proliferator-activated receptor (PPAR) gamma protein translation by suppressing the mammalian target of rapamycin (mTOR) signaling pathway in cultured adipocytes. Journal of Biological Chemistry 285: 6595–6603. doi: 10.1074/jbc.M109.068502 PMID: 20038583

106. Wang M, Amano SU, Flach RJR, Chawla A, Aouadi M, Czech MP (2013) Identification of Map4K4 as a novel suppressor of skeletal muscle differentiation. Molecular and Cellular Biology 33: 678–687. doi: 10.1128/MCB.00618-12 PMID: 23207904

107. Gineitis D, Treisman R (2001) Differential usage of signal transduction pathways define two types of serum response factor target gene. Journal of Biological Chemistry 276: 24531–24539. PMID: 11342553

108. Boxer LM, Prywes R, Roeder RG, Kedes L (1989) The sarcomeric actin CArG-binding factor is indistinguishable from the c-fos serum response factor. Molecular and Cellular Biology 9: 515–522. PMID: 2710114

109. Lee TC, Chow KL, Fang P, Schwartz RJ (1991) Activation of skeletal actin gene transcription: The cooperative formation of serum response factor-binding complexes over positive cis-acting promoter serum response elements displaces a negative-acting nuclear factor enriched in replicating myoblasts and nonmyogenic cells. Molecular Biology of the Cell 11: 5090–5100.

110. Minty A, Kedes L (1986) Upstream regions of the human cardiac actin gene that modulate its transcription in muscle cells: Presence of an evolutionarily conserved repeated motif. Molecular and Cellular Biology 6: 2125–2136. PMID: 3785189

111. Miano JM, Long X, Fujiwara K (2007) Serum response factor: Master regulator of the actin cytoskeleton and contractile apparatus. American Journal of Physiology Cell Physiology 292: C70–C81. PMID: 16928770

Related topics
Byungwhi Kong
University of Arkansas (USA)
Kentu R. Lassiter
University of Arkansas (USA)
Alissa Welsher
Sami Dridi
University of Arkansas (USA)
Walter Bottje
University of Arkansas (USA)
Show more
Join to be able to comment.
Once you join Engormix, you will be able to participate in all content and forums.
* Required information
Would you like to discuss another topic? Create a new post to engage with experts in the community.
Create a post
Join Engormix and be part of the largest agribusiness social network in the world.