ORIGINAL RESEARCH article

Front. Vet. Sci., 04 August 2020

Sec. Animal Nutrition and Metabolism

Volume 7 - 2020 | https://doi.org/10.3389/fvets.2020.00458

Muscle Metabolome Profiles in Woody Breast-(un)Affected Broilers: Effects of Quantum Blue Phytase-Enriched Diet

  • 1. Center of Excellence for Poultry Science, University of Arkansas, Fayetteville, AR, United States

  • 2. Department of Animal Sciences, University of Arkansas, Fayetteville, AR, United States

  • 3. Biological and Small Molecule Mass Spectrometry Core, Department of Chemistry, University of Tennessee, Knoxville, Knoxville, TN, United States

  • 4. AB Vista, Marlborough, United Kingdom

Article metrics

View details

28

Citations

5,3k

Views

2,4k

Downloads

Abstract

Woody breast (WB) myopathy is significantly impacting modern broilers and is imposing a huge economic burden on the poultry industry worldwide. Yet, its etiology is not fully defined. In a previous study, we have shown that hypoxia and the activation of its upstream mediators (AKT/PI3K/mTOR) played a key role in WB myopathy, and supplementation of quantum blue (QB) can help to reduce WB severity via modulation of hypoxia-related pathways. To gain further insights, we undertook here a metabolomics approach to identify key metabolite signatures and outline their most enriched biological functions. Ultra performance liquid chromatography coupled with high resolution mass spectrometry (UPLC–HRMS) identified a total of 108 known metabolites. Of these, mean intensity differences at P < 0.05 were found in 60 metabolites with 42 higher and 18 lower in WB-affected compared to unaffected muscles. Multivariate analysis and Partial Least Squares Discriminant analysis (PLS-DA) scores plot displayed different clusters when comparing metabolites profile from affected and unaffected tissues and from moderate (MOD) and severe (SEV) WB muscles indicating that unique metabolite profiles are present for the WB-affected and unaffected muscles. To gain biologically related molecule networks, a stringent pathway analyses was conducted using IPA knowledge-base. The top 10 canonical pathways generated, using a fold-change −1.5 and 1.5 cutoff, with the 50 differentially abundant-metabolites were purine nucleotide degradation and de novo biosynthesis, sirtuin signaling pathway, citrulline-nitric oxide cycle, salvage pathways of pyrimidine DNA, IL-1 signaling, iNOS, Angiogenesis, PI3K/AKT signaling, and oxidative phosphorylation. The top altered bio-functions in term of molecular and cellular functions in WB-affected tissues included cellular development, cellular growth and proliferation, cellular death and survival, small molecular biochemistry, inflammatory response, free radical scavenging, cell signaling and cell-to-cell interaction, cell cycles, and lipid, carbohydrate, amino acid, and nucleic acid metabolisms. The top disorder functions identified were organismal injury and abnormalities, cancer, skeletal and muscular disorders, connective tissue disorders, and inflammatory diseases. Breast tissues from birds fed with high dose (2,000 FTU) of QB phytase exhibited 22 metabolites with significantly different levels compared to the control group with a clear cluster using PLS-DA analysis. Of these 22 metabolites, 9 were differentially abundant between WB-affected and unaffected muscles. Taken together, this study determined many metabolic signatures and disordered pathways, which could be regarded as new routes for discovering potential mechanisms of WB myopathy.

Introduction

Broiler (meat-type) chickens play a key role in worldwide meat production and support the livelihoods and food security of billions of people (1). In fact, poultry meat is highly regarded worldwide as one of the most efficient food sources with high nutrient and organoleptic quality, inexpensive, and without religious taboos (2). However, the emerging woody breast (WB) myopathy is significantly impacting modern broilers and is imposing a huge economic burden on the poultry industry worldwide due to on-farm culling and mortality, down-grading and condemnation at processing, as well as rejection from human consumption (3–7).

WB incidence has increased drastically from an average of 5% in 2012 to 29% in 2015 and has been reported to affect up to 50% in other flocks (8). WB myopathy is emerging on a global scale, already present in Finland (9), France (10), Italy (5), Spain (11), Brazil (12, 13), United Kingdom (14), Japan (15), and in several other countries, and it is negatively impacting global chicken meat production and quality.

Although the etiology of WB is still not well-defined, evidence indicated multifocal degeneration and necrosis of muscle tissue with infiltration of inflammatory cells (9, 16). Lesions associated with the myopathy appear to be aseptic, superficially-located, and include muscle fiber fragmentation, hyalinization, and swelling with replacement by fibrous connective tissue, as well as an influx of macrophages and lipid infiltration (9, 17). These clinical and microscopic changes result in palpable severe hardness of the breast muscle. At molecular and cellular levels, several recent omics studies including transcriptomics (18, 19), metabolomics (8, 20), and proteomics (21) have been conducted and identified numerous potential contributing factors. Of particular interest, hypoxia, oxidative stress, fiber-type switching, cellular damage, and altered intracellular calcium were predicted to play a key role in WB myopathy. In 2016, Abasht's group identified several key metabolites in 48d-old WB-affected purebred lines and commercial broilers (8). In a more recent study using 38d-old Ross 708 broilers and NMR technique, Wang et al. reported alteration of several metabolites including amino acids in WB-affected breasts (20).

In attempt to identify a mechanism-based strategy to reduce the severity of WB, we recently demonstrated that systemic and local hypoxia along with an activation of HIF-1α and its upstream mediators PI3K/AKT/mTOR pathways are responsible for the development of WB myopathy in 56d-old Cobb500 broilers (22). Phytase (quantum blue, QB)-enriched diets at 2,000 FTU/kg reduced the severity of WB via modulation of hypoxia- and oxygen homeostasis-related pathway (22). In the present study, we aimed to gain further insights in the WB etiology and the mode of QB action by using a mass spectrometric metabolomics approach.

Materials and Methods

Ethics Statements

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, and the protocol was approved by the University of Arkansas Animal Care and Use Committee (#16084).

Birds and Diets

The experimental design was previously described (22). Briefly, male broilers (Cobb 500, n = 576) were weighed at day one of hatch and randomly assigned to 6 body-weight matched groups in 48 floor pens (12 birds/pen, 8 pens/group) in an environmentally controlled house with ad libitum access to feed and water. Birds were fed, for 56 days, one of six dietary treatments in a complete randomized design. The diets were a nutrient adequate positive control (PC) diet formulated to meet Cobb 500 nutrition requirements. Myo-inositol (MI, Sigma-Aldrich, St. Louis, MO) was added to the PC diet at 0.30% to create a second diet (PC + MI). The third one was considered the negative control (NC) diet with a reduction of available phosphorus (avP), calcium and sodium by 0.15, 0.16, or 0.03%, respectively. The NC diet was then supplemented with 500, 1,000, or 2,000 phytase units (FTU)/kg to create diets four (NC+500FTU), five (NC+1,000 FTU) and six (NC+2,000 FTU), respectively. The phytase was QB (AB Vista, Marlborough, UK) with an expected activity of 5,000 FTU/g. The diet composition was reported previously (22).

WB Scoring

Birds (n = 85–90/group) were processed at d56 using a commercial inline system at the University of Arkansas Pilot Processing Plant (UAPPP). Breast filets were blind analyzed and macroscopically scored by one trained person (who did not know the different treatments), and classified to WB categories to the degree: 0, normal (NORM); 0.5–1.5, moderate (MOD) with mild hardening in the cranial area; and 2–3, severe (SEV) with severe hardening and hemorrhagic lesions in the cranial region.

Sample Collection and Preparation

Breast muscle samples were collected from cranial surface (S1) of unaffected birds and from cranial S1 (woody woody, WW) and caudal S2 (woody normal WN, apparent healthy area) area of WB-affected birds for both MOD and SEV category and from each treatment (n = 8, 96 WB-affected with 48 MOD and 48 SEV, and 48 non-affected birds) (Figure 1). Tissues were ground snap frozen in liquid nitrogen, and stored at −80°C until further analysis. Metabolites were extracted with 1.5 mL of extraction solvent (40:40:20 HPLC grade methanol: acetonitrile: water with formic acid at a final concentration of 0.1 M) pre-chilled at 4°C, and incubated at −20°C for 20 min. Samples were centrifuged (13,300 g, 5 min, 4°C) and supernatants were collected. Solvent was evaporated under a stream of nitrogen and metabolites were suspended with 300 μL of ultrapure water.

Figure 1

Figure 1

Woody breast categories and sampling locations. Birds were processed at d56, and breast filets were macroscopically scored. NORM, normal breast without harden area; MOD, moderate woody breast with a score between 0.5 and 1.5; SEV, severe woody breast with a score between 2 and 3. S1, cranial surface; S2, caudal surface; WN, woody normal; WW, woody woody.

Ultra-High Performance Liquid Chromatography—High Resolution Mass Spectrometry (UHPLC–HRMS) Metabolomic Analysis

UHPLC-HRMS analysis has been described previously (23). Briefly, metabolites were separated on a Dionex UltiMate 3000RS (Sunnyvale, CA) by injecting a 10 μL sample on a Synergy reverse phase Hydro-RP 100 Å, 100 mm × 2.00 mm, 2.5 μm pore size LC column (Phenomenex, Torrance, CA) kept at 25°C. The untargeted metabolomics method, adapted from (24), ran for 26 min with the application of multistep gradient. Two HPLC grade solvents were used in gradient steps to separate the analytes. Solvent A contained 97:3 H2O: MeOH with 11 mM tributylamine and 15 mM acetic acid and solvent B was 100% MeOH. The gradient was ran as follows: 0 min, 0% B; 5 min, 20% B; 13 min, 55% B; 15.5 min, 95% B; 19 min, 0% B; 25 min, 0% B with a flow rate of 200 μL/min. The eluent was introduced into the mass spectrometer via an electrospray ionization (ESI) source conjoined to an Exactive™ Plus Orbitrap Mass Spectrometer (Thermo Scientific, Waltham, MA) under established parameters of aux gas: 8; sheath gas: 25; sweep gas: 3; spray voltage: 3.00 kV; and capillary temperature: 300°C. The mass spectrometer parameters were set as follows: resolution: 140,000; automatic gain control (AGC): 3 × 106; maximum IT time: 100; scan range: 72–1,200 m/z. Raw data were acquired from the Xcalibur MS software (Thermo Electron Corp, Waltham, MA) and converted to mzML format by ProteoWizard tool MSConverter (25, 26). The converted data were analyzed using MAVEN (27) and peaks were annotated with a maximum allowed error of 5 ppm. Area under the chromatographic curve was integrated based upon an in-house verified list of metabolites using exact m/z and known retention times (28). All metabolite values were normalized based on the exact mass of the tissue extracted prior to all statistical calculations.

Data Analysis

Metabolites showing >1.5 fold differences and <0.05 p-value in the comparison between WB-affected and unaffected birds were considered differentially abundant. Ingenuity Pathway Analysis (IPA; Qiagen, Valencia, CA; http://www.ingenuity.com) software was used for functional annotation, canonical pathways analysis, upstream analysis, and network discovery based on human metabolome database (29). Heat maps were made by Cluster 3.0 (30) and Javatreeview 1.1 (31) using log2 transformed data. P-values were calculated using Student's T-test. Partial Least Squares Discriminant Analysis (PLS-DA) and Variable Importance in Projection (VIP) scores were constructed using the statistical package DiscriMiner in R version 3.6.1 (https://cran.r-project.org). These VIP scores provide a metric for determining how much influence a metabolite has on the group separation seen in the PLS-DA plots. A value > 1 indicates that the metabolite contributes to the group differentiation, and this was considered as a significant VIP score.

Results

Multivariate Analysis and Comparative Metabolomics Profile in WB-Affected and Unaffected Birds

The untargeted metabolomics profiling analyses identified 108 known metabolites and has been submitted to MetabolLights database (https://www.ebi.ac.uk/metabolights). Formal testing, pre-metabolite by two-sample t-tests, showed that there were 60 metabolites with differences in mean intensity at P < 0.05 (Table S1). Of these, 42 were at higher concentrations and 18 were at lower levels in WB-affected birds compared to their unaffected counterparts (Table S1). Most altered metabolites belonged to the nucleosides, nucleotides and analogs (35.06%), amino acids, peptides and analogs (22.95%), carbohydrates and carbohydrate conjugates (13.11%), dicarboxylic acids (4.92%), pyridines and derivatives (4.92%), and 1.64% of each of phosphate ester, furanones, organosulfonic acid, keto acids, fatty acids, hydroxyl acids, flavonoids, peptidomimetics, alcohols, amines, and lactams. When the metabolite profiles for each WB-affected and unaffected group were represented using a PLS-DA scores plot, two clusters were observed (Figure 2A). Similarly, PLS-DA analyses and score plots displayed different clusters when comparing metabolites profile from MOD and SEV WB-affected muscles (Figure 2B). Within apparent healthy muscles obtained from control S1 region, and MOD and SEV S2 regions, PLS-DA scores plots showed 3 separate clusters (Figure 2C). When all the samples were plotted together, 5 clusters were identified (Figure 2D). This indicates that a clear separation of affected- and unaffected-tissues on one hand, and WB category (MOD vs. SEV) on the other hand can be achieved and that unique metabolite profiles are present for the WB-affected and unaffected muscles. The PLS-DA combines features from principle component analysis and multiple regression and transforms a large number of potentially correlated variables into a smaller number of orthogonal variables (i.e., PLS1, PLS2) that discriminates between classes. In order to determine which metabolites drove the separation among the groups, each metabolite was assigned a VIP score representing its importance to the PLS-DA model with a score >1 indicating a significant contribution to the separated clusters observed. As shown in Table 1, 43 metabolites were found to be driving the separation between WB-affected and unaffected muscles, and 38 between MOD and SEV affected tissues in the PLS-DA plot. Of these, 22 common metabolites are involved in all the segregation (Table 1).

Figure 2

Figure 2

Partial Least Squares Discriminate Analysis (PLS-DA) scores plots. PLS-DA was constructed using the statistical package DiscriMiner in R version 3.6.1 and displayed different clusters when comparing metabolite profiles between WB-affected and unaffected birds (A), MOD and SEV WB (B), WN-MOD, and WN-SEV (C), and all the groups together (D).

Table 1

MetabolitesC vs. WBWN1.5 vs. WN3WW1.5 vs. WW3
Flavone+++
Glucose6phosphate+++
Tricarballylicacid+++
Guanidoaceticacid+++
Sedoheptulose1phosphate+++
UTP+++
3−Phosphoglycerate+++
GDP+−+
Histamine+++
NAD++−+
Deoxycytidine+++
Trehalose/Sucrose++−
2, 3−Biphosphoglycerate++−
Homocamosine+−+
Vanillin+−−
Phosphorylethanolamine+++
AMP/dGMP++−
Salicylate+++
Cytidine++−
UDP−N−acetylglucosamine+++
Myo−inositol+++
Methylsuccinicacid++−
Cysteine+−+
Glutaricacid++−
CDP−ethanolamine+++
Cystathionine+++
2−Oxoisovalerate+++
4−Pyridoxoate+++
2−Isopropylmatate+++
S−Adenosyl−L−homocysteine++−
Dephospho−CoA+++
2−Dehydro−D−gluconate+−−
Glucosaminephosphate+−+
Deoxyuridine+−+
Alpha−Ketoglutarate+−+
Citrulline+++
Hydroxypherrylpyruvate+++
Homocysteine+++
Shikimate+−+
Pyroglutamicacid+−−
Dimethylglycine+−+
NADH+−+
Ophthalmate+−−
Taurine−++
ADP−++
AICAR−+−
CDP−++
Proline−+−
NADH−+−
Glycerate−+−
Dopamine−+−
1−Methyladenosine−++
Deoxyribosephosphate−++
N−Carbamoyl−L−aspartate−−+
Flavone+++
Glucose 6 phosphate+++
Tricarballylic acid+++
Guanidoacetic acid+++
Sedoheptulose 1 phosphate+++
UTP+++
3-Phosphoglycerate+++
GDP+−+
Histamine+++
NAD++−+
Deoxycytidine+++
Trehalose/Sucrose++−
2,3-Biphosphoglycerate++−
Homocamosine+−+
Vanillin+−−
Phosphorylethanolamine+++
AMP/dGMP++−
Salicylate+++
Cytidine++−
UDP-N-acetylglucosamine+++
Myo-inositol+++
Methyl succinic acid++−
Cysteine+−+
Glutaric acid++−
CDP-ethanolamine+++
Cystathionine+++
2-Oxoisovalerate+++
4-Pyridoxoate+++
2-Isopropylmatate+++
S-Adenosyl-L-homocysteine++−
Dephospho-CoA+++
2-Dehydro-D-gluconate+−−
Glucosamine phosphate+−+
Deoxyuridine+−+
Alpha-Ketoglutarate+−+
Citrulline+++
Hydroxypherrylpyruvate+++
Homocysteine+++
Shikimate+−+
Pyroglutamic acid+−−
Dimethylglycine+−+
NADH+−+
Ophthalmate+−−
Taurine−++
ADP−++
AICAR−+−
CDP−++
Proline−+−
NADH−+−
Glycerate−+−
Dopamine−+−
1-Methyladenosine−++
Deoxyribose phosphate−++ +
N-Carbamoyl-L-aspartate−−+

Metabolites driving cluster separations based on VIP scorea.

a

VIP, Variable Importance in Projections scores were obtained from PLSDA using the DiscriMiner software package. +, indicates a VIP score >1 and − indicates a VIP score <1.

Metabolic and Canonical Pathways Analysis

For more stringent pathway analyses and to gain biologically related molecule networks, the above identified metabolites (108) were mapped into the IPA knowledge-base and analyzed to outline the most enriched biological functions. Using a cut-off of FDR adjusted p-value < 0.05 and a fold-change between −1.5 and 1.5, IPA analysis identified 50 differentially abundant (DA) metabolites between WB-affected and unaffected birds (Table 2). Among the 50 DA metabolites, 37 were significantly higher and 13 were significantly lower in the breast muscle of WB-affected birds compared to their unaffected counterparts (Table 2). The top 10 canonical pathways generated with DA-metabolites were purine nucleotide degradation and de novo biosynthesis, sirtuin signaling pathway, citrulline-nitric oxide cycle, salvage pathways of pyrimidine DNA, IL-1 signaling, iNOS, Angiogenesis, PI3K/AKT signaling, and oxidative phosphorylation (Table 3). The top altered bio-functions in term of molecular and cellular functions in breast muscle comparison between WB-affected and unaffected birds included cellular development, cellular growth and proliferation, cellular death and survival, small molecular biochemistry, inflammatory response, free radical scavenging, cell signaling and cell-to-cell interaction, cell cycles, and lipid, carbohydrate, amino acid, and nucleic acid metabolisms (Table 4). The top disease and disorder functions identified by IPA analysis were organismal injury and abnormalities, cancer, skeletal and muscular disorders, connective tissue disorders, and inflammatory diseases (Table 4). According to the canonical pathway analysis, 13 prediction networks were identified. The first top predicted network used 15 metabolites and focuses on amino acid metabolism, carbohydrate metabolism, molecular transport, and small molecules biochemistry (Figure 3). Metabolites like S-adenosylhomocysteine, UDP-N-acetylglucosamine, L-cysteine, trehalose, and deoxycytidine that characterized by 11, 5.34, 3.84, 2.83, and −3.32-fold changes in WB-affected compared to unaffected tissues, had hubs at the signaling molecule epidermal growth factor receptor (EGFR) (Figure 3A). EGFR is reported to be involved in cancer, apoptosis, and skeletal and muscular disorders (32–34). Network 2 used 14 molecules and centers on energy production, amino acid metabolism, and small molecule biochemistry and has a nexus to the signaling molecule extracellular-signal-related kinase 1/2 (Figure 3B) which has been shown to play pivotal roles in cancer, cell death, and muscle disorders (35, 36).

Table 2

HMDB IDMetabolite nameFold changeP-value
HMDB0001517AICAR3.2550.000153
HMDB0000044Ascorbic acid3.1010.000789
HMDB0001564CDP ethanolamine3.2390.0000712
HMDB0000058Cyclic AMP1.8980.0442
HMDB0000089Cytidine2.7820.00021
HMDB0011732D-2-ketogluconic acid3.490.000973
HMDB0000625D-gluconic acid3.0120.00353
HMDB0061388Dimethyl-alpha-ketoglutarate2.7490.0218
HMDB0060475DL-glutamic acid1.6940.0245
HMDB0001227dTMP3.4560.000453
HMDB0001409dUMP3.4430.000214
HMDB0000134Fumaric acid2.0860.0389
HMDB0000661Glutaric acid2.0030.0426
HMDB0001397GMP2.1490.00243
HMDB0000133Guanosine2.9370.0000876
HMDB0000157Hypoxanthine2.2160.00439
HMDB0000517L-arginine1.690.0159
HMDB0000191L-aspartic acid1.8750.0189
HMDB0000099L-cystathionine4.5760.000028
HMDB0000574L-cysteine3.840.000237
HMDB0000744Malic acid1.9690.0341
HMDB0001844Methylsuccinic acid1.980.0432
HMDB0001138N-acetyl-L-glutamate3.8380.000725
HMDB0006029N-acetyl-l-glutamine2.6660.012
HMDB0000828N-carbamoyl-L-aspartic acid3.0590.0363
HMDB0005765Ophthalmic acid2.1930.0062
HMDB0000224Phosphorylethanolamine6.2650.0000338
HMDB0000267Pyroglutamic acid2.9240.00764
HMDB0000939S-adenosylhomocysteine2.4390.000472
HMDB0060509Sedoheptulose 1-phosphate110.0000145
HMDB0000251Taurine3.0770.0000885
HMDB0000975Trehalose2.8310.0116
HMDB0000290UDP-N-acetylglucosamine5.3450.0000258
HMDB0000300Uracil2.1170.0148
HMDB0000296Uridine2.1430.00181
HMDB0000299Xanthosine2.2310.00221
HMDB0001554Xanthosine monophosphate2.4430.00699
HMDB0031193Tricarballylic acid−5.0570.0243
HMDB0000904Citrulline−2.2550.00496
HMDB0060465AMP−1.9830.00846
HMDB0000014Deoxycytidine−3.3250.0022
HMDB0000092Dimethylglycine−1.9320.0161
HMDB0003075Flavone−18.3840.0106
HMDB0001201GDP−2.2830.0248
HMDB0001254Glucosamine phosphate−2.4370.0171
HMDB0001401Glucose-6-phosphate−9.3830.000114
HMDB0000128Guanidoacetic acid−5.8520.0229
HMDB0000870Histamine−2.1650.0104
HMDB0000745Homocarnosine−2.4140.00113
HMDB0000902NAD+−2.8090.00122

Differentially abundant metabolites in woody breast-affected and unaffected birds.

Table 3

Canonical pathwaysMolecules-log (p-value)Ratio
Purine nucleotides degradation IIGMP, guanosine, hypoxanthine, NAD, xanthosine, xanthosine monophosphate6.000.353
Purine nucleotides de novo biosynthesis IIAICAR, fumaric acid, GDP, GMP, L-aspartic acid, NAD, xanthosine monophosphate5.580.233
Sirtuin signaling pathwayCitrulline, fumaric acid, glucose-6-phosphate, L-aspartic acid, malic acid, NAD4.160.182
Citrulline-nitric oxide cycleCitrulline, fumaric acid, L-arginine, L-aspartic acid4.160.364
Salvage pathways of pyrimidine deoxyribonucleotidesDeoxycytidine, dTMP, dUMP, uracil3.990.333
IL-1 signalingCyclic AMP, GDP2.820.667
iNOS signalingCitrulline, L-arginine2.820.667
Inhibition of angiogenesisCitrulline, L-arginine2.310.400
PI3K/AKT signalingCitrulline, L-arginine2.310.400
Oxidative phosphorylationFumaric acid, NAD1.690.200

Top canonical pathways enriched by observed metabolite alteration in WB myopathy.

Table 4

Diseases and functionsP-value# Molecules*
Organismal injury and abnormalities2.55E−06-4.54E−0226
Cancer2.55E−06-3.93E−0220
Skeletal and muscular disorders6.09E−03-3.41E−0210
Connective tissue disorders6.09E−03-3.41E−026
Inflammatory diseases6.09E−03-4.54E−026
Molecular and Cellular FunctionsP-value# Molecules*
Small molecule biochemistry1.46E−04-4.54E−0223
Cellular growth and proliferation3.38E−05-4.47E−0221
Cellular death and survival3.68E−05-4.54E−0219
Molecular transport1.46E−04-4.54E−0219
Cellular development3.38E−05-4.47E−0217
Lipid metabolism1.46E−04-4.54E−0217
Inflammatory response1.52E−03-4.54E−0217
Cellular function and maintenance3.68E−04-4.22E−0215
Free radical scavenging1.46E−04-3.03E−0214
Carbohydrate metabolism1.20E−03-4.54E−0214
Nucleic acid metabolism3.37E−03-4.54E−0213
Cellular assembly and organization3.68E−04-3.79E−0212
Connective tissue development and function3.90E−05-4.54E−0212
Cell morphology3.68E−04-4.47E−0211
Amino acid metabolism6.67E−04-4.47E−0211
Cell signaling6.80E−03-3.79E−0210
Cell-to-cell signaling and interaction3.68E−04-4.54E−0210
Protein synthesis4.82E−03-3.41E−029
Vitamin and mineral metabolism6.80E−03-3.41E−026
Cell cycle2.99E−03-4.36E−027
Immune cell trafficking3.37E−03-4.54E−027
Skeletal and muscular system development and function7.27E−06-4.47E−027
Energy production2.99E−03-4.36E−026

Top diseases, molecular, cellular, and physiological functions for DE-metabolites in WB-affected compared to unaffected birds.

*

Number of molecules with each network. The molecule networks were generated from association with IPA database and our input data. The P-value was calculated using Fisher's exact Test by IPA.

Figure 3

Figure 3

The top two networks (A,B) built with IPA program from metabolomics data that was determined on WB-affected muscles. Metabolites in red were upregulated and metabolites in green were down-regulated in WB-affected tissues. Networks constructed from the IPA knowledgebase by connecting the altered molecules are not limited by canonical pathway boundaries. EGFR, epidermal growth factor receptor; ERK, extracellular signal-regulated kinase.

Effect of QB on WB Incidence and Muscle Metabolome Profile

QB supplementation at high doses (1,000 and 2,000 FTU) reduced the severity of WB by ~5% compared to the positive control [Table 5, (22)]. Multivariate data analysis and PLS-DA score plots displayed 6 different clusters when comparing all diet supplementations (Figure 4A). Samples from birds fed with high dose of QB (NC+2,000 FTU) form a distinct and segregated cluster from the rest of the groups (Figure 4A). Further analysis using VIP scores revealed that 30 metabolites drove the separation among the groups. As illustrated in Figure 4B and Figure S1, heat maps showed the relative levels of metabolites in each group compared to the PC control. Breast tissues from QB (2,000 FTU)-fed birds exhibited 22 metabolites with significantly different levels compared to the PC group (Figure 4B). These metabolites belong to amino acids (47.80%), nucleosides (26.13%), carboxylic acids (8.65), carbohydrates (4.35%), phenylpyruvic acids (4.37%), amines (4.35), and indoles (4.35%). Of these 22 metabolites, only two (S-adenosyl-L-homocysteine and arginine) were found to be significantly higher and seven (tricarballylic acid, NAD+, glucosamine phosphate, citrulline, histamine, AMP/dGMP, and 4-pyridoxate) were found to be significantly lower in WB-affected compared to unaffected birds.

Table 5

ParametersaDiets
PCPC+MIONCNC+500 FTUNC+1,000 FTUNC+2,000 FTU
WB categoryb
NORM10.348.1924.245.085.6810.54
MOD68.9668.8569.6971.1878.4273.68
SEV20.6822.956.0723.7415.9015.78
Total incidence89.6491.875.7694.9294.3189.46

Effects of quantum blue supplementation on WB incidence (%) in broilers.

a

The scoring and incidence of WB were determined using 8 replicate pens/treatment and 8 birds/pen.

b

MOD, moderate; NORM, normal; SEV, severe.

Figure 4

Figure 4

Partial Least Squares Discriminate Analysis (PLS-DA) scores plots and heat maps of relative levels of most significant metabolites modulated by high dose quantum blue-enriched diets. (A) PLS-DA was constructed using the statistical package DiscriMiner in R version 3.6.1 and displayed different clusters when comparing metabolite profiles between all diets, (B) heat maps was build using log2 ratio and VIP score and comparing each diet to the PC group.

Discussion

WB is a stumbling block and significant economic burden to the poultry industry worldwide and for which there is no effective preventive strategy due to its unknown etiology. Several elegant high-throughput omics studies predicted the involvement of various pathways (18, 21, 37) and several seminal managerial and nutritional strategies including our own were conducted to reduce the severity of WB myopathy (22, 38, 39). Because metabolites play a key role in the maintenance and growth of organisms and are considerably useful to illustrate the underlying molecular disease-causing mechanisms (40, 41), we used here a metabolomics approach to identify key metabolite markers and outline the most enriched biological functions in WB myopathy.

Our high-throughput metabolomics analysis and PLS-DA scores plots demonstrated that distinct metabolic profiles can be detected for WB-affected and unaffected tissues. The altered metabolites encompassed all major molecular categories and critical metabolic pathways from amino sugar and nucleotide sugar metabolism, methionine and choline metabolism, and purine/pyrimidine salvage to Krebs and urea cycles. Arranging these DA metabolites using both IPA and KEGG metabolic pathways into a single map revealed a high level of connectivity and inter-dependence (summarized in Figure 5).

Figure 5

Figure 5

Interconnected altered metabolic pathways in WB-affected muscles. Altered metabolites from the top canonical pathway analysis along with their cross-independence in KEGG pathways. Red text metabolites represent elevated mean intensity in WB-affected, whereas green text represent reduced intensity. The figure is representative and not for scale.

Of particular interest and in agreement with previous study (8), the 10-fold decrease in glucose-6-phosphate indicates a potential dysregulation of both the glycolysis and glycogenesis in WB-affected muscle. It has been shown that a defect in glucose-6-phosphate concentration was a major contributor to reduce whole-body insulin-mediated glucose disposal rates and to reduce insulin action on glycogen synthase (42). Although it was not measured in the present study, dysregulation of glucose-6-phosphate levels indicated an alteration of glycogen content in muscle with WB (8, 43). This is supported here by the increased levels of AICAR, which has been reported to alter glycogen content via alteration of glycogen synthase and glycogen phosphorylase activities (44). AICAR is an analog of adenosine monophosphate and is capable of stimulating the master energy sensor, AMP-dependent protein kinase (AMPK), which is higher in WB-affected compared to unaffected muscles (data not shown) indicating an intracellular ATP depletion status. AICAR plays a crucial role in the recruitment of ATP-sensitive K channels to the sarcolemma and the regulation of calcium content (45), and calcium overloaded was found to be a hallmark of muscle dystrophies (46). AICAR has also been shown to alter mTOR signaling pathway in skeletal muscle (47) and our previous study demonstrated that WB myopathy is associated with mTOR dysmetabolism-induced hypoxia (22).

The dysregulation of glucose metabolism (reduction of glucose-6-phosphate) was accompanied by an accumulation of trehalose, ascorbate, and gluconate in WB-affected tissues. The presence of trehalose is surprising as it has not been found in higher vertebrate to date even though the enzyme trehalase required for trehalose hydrolysis has been reported (48). Trehalose is synthesized as a stress-responsive factor and has been implicated in cyto-protection against hypoxia, neurodegenerative disease, and oculopharyngeal muscular dystrophy (49–51). The glucose-derived molecule ascorbate functions as an enzyme cofactor and antioxidant, and ~40% of the body's ascorbate is stored in skeletal muscle (52). However, its significant increase in our experimental conditions suggested that this might due to excess of oxidative stress in WB-affected muscles (8). As high doses of antioxidants interfere with muscle signaling pathways and can be detrimental (53), it is possible that high concentration of ascorbate along with low concentration of flavone (−18 fold changes) might exacerbate the WB myopathy. The alteration of gluconate metabolism (increase of gluconate and 2-keto-L-gluconate levels) in WB suggested an activation of amino acid catabolism in order to provide carbon skeleton for anaplerosis which is evident here by increased concentrations of glutamate (converted to α-ketoglutarate) and aspartate (converted to oxaloacetate), main components in TCA cycle. The increased levels of α-ketoglutarate, fumarate, and malate in this study suggested a dysregulated cataplerosis pathway resulting in accumulation of these metabolites in the mitochondrial matrix. Interestingly, one of the 2-keto-L-gluconate binding target was the β-subunit of the mitochondrial ATP synthase resulting in reduced mitochondrial oxygen consumption and oxidative phosphorylation (54), which has been predicted by IPA analysis in this study. Furthermore, 2-keto-L-gluconate has been shown to regulate acid-base balance (55), which could be the reason for higher ultimate pH seen in WB-affected muscles (37). The accumulation of the above mentioned metabolites in association with mitochondrial dysfunction can lead to high levels of glucaric (4-fold changes) and glutaric acids (2-fold changes). Glucaric and glucoronic acids can act as acidogen and metabotoxin and their accumulation has been reported to cause damage to several organs via up-regulation of the production and release of endogenous inflammatory signaling molecules (56). Glutaric acid administration has been shown to alter mitochondrial complexes I-III and II-III and impair energy metabolism in rodent skeletal muscle (57).

The pentose phosphate pathway (PPP), paralleling an approximate 11-fold increase in sedoheptulose 1/7 phosphate, was altered in WB-affected tissues. Previous studies showed that hydrogen peroxide treatments elevated sedoheptulose 1/7 phosphate levels in hepatoma HepG2 cells (58). Because the PPP activation was reported to be involved in cyto-protection against oxidative damage (59), our data suggested that oxidative stress might increase the flux of glucose into the PPP pathway. Consistent with an altered redox homeostasis, the levels of antioxidant-associated metabolites was affected. For instance, the major low molecular weight antioxidant glutathione as well as its trans-sulfuration sources were increased in WB-affected tissues suggesting high metabolic stress (peroxides and free radicals) and intracellular detoxification of xenobiotics and cell damaging compounds (60). Our data suggested also that chicken muscle cells might serve as a continuous extrahepatic source of glutathione during local trauma where the release of stress hormone may greatly enhance the release of the antioxidant (61). Furthermore, the aminoethane sulfonic acid, taurine, was higher in WB-affected compared to unaffected muscles. In addition to its antioxidant activity, taurine was found to be involved in mitochondrial t-RNA and mitochondrial protein synthesis (62). Taurine also plays a crucial role in calcium homeostasis and muscle contractile function (63, 64), both of which were defective in WB myopathy.

The alteration of methionine and choline metabolism was demonstrated by higher levels of s-adenosyl-homocysteine (SAH), phosphorylethanolamine and CDP-ethanolamine in WB-affected tissues. S-adenosyl-L-homocysteine is the major donor of methyl groups in the synthesis of phospholipids, nucleotides, epinephrine, carnitine, and creatine. It has been reported that, when present at high levels, SAH can act as an immunotoxin (disrupts, limits, or destroys the function of immune cells) and it has been linked to metabolic alteration, cardiovascular disease, and end-stage renal disease (65, 66). Phosphorylethanolamine was found to be linked to oxidative stress in cultured rat lenses (67). The diacylglycerol-consuming CDP-ethanolamine pathway is a major route for phosphatidylethanolamine synthesis, which constitutes the second major abundant phospholipids (68). In addition to controlling muscle diacylglycerol levels, emerging evidence indicated that CDP-ethanolamine pathway plays a significant role in regulating muscle function and mitochondrial biology (69).

A noteworthy additional point is the alteration of the nucleotide (purine and pyrimidine) de novo synthesis, catabolism, and salvage. These processes are energy-expensive, but they are crucial for every living cell, and their alteration may lead to significant functional consequences. For example, guanosine, which plays a key role in muscle contraction and signaling pathways (70), has been reported to be toxic at high concentration (71). High levels of serum uracil concentration has been used as a predictor of severe fluoropyrimidine-associated toxicity (72). Similarly high concentrations of xanthine and hypoxanthine have been found to be toxic (73, 74). Because hypoxanthine is very low under physiological conditions, our data suggested that hypoxanthine was build up with the progression of WB myopathy with concomitant production of reactive oxygen species (75, 76). Its high levels have been shown to be associated with muscle damage (77). The alteration of the aforementioned metabolic pathways in WB myopathy lead probably to excessive consumption of ATP and NAD+ which is supported by the activation of AMPK (data not shown) and a decreased levels of NAD+ (our metabolomics data), respectively. It has been shown that oxidative stress induced a deficit in NAD+ which in turn promoted generation of methylglyoxal that glycated nucleic acid and protein and conduced to advanced glycation end product (AGE) formation (78, 79).

While some of our altered metabolites have been reported by Abasht and co-workers (for fumarate, malate, sedoheptulose 1/7 phosphate, glutathione, taurine, and nucleotides) (8), by Wang and colleagues (for taurine, hypoxanthine, and NAD+) (20), and by Soglia et al. (80) (for taurine, uracil, and hypoxanthine), several metabolites such as AICAR, trehalose, flavones, etc. constitute new molecular signatures, and they have not been reported previously. The common altered metabolites observed in all three studies might constitute conserved metabolic pathways in WB myopathy across all broiler strains. The distinct metabolites observed in our study might be strain- or age-specific as we used 56d-old Cobb500 compared to 38d-old Ross708 by Wang et al. (20), and 47-48d-old undisclosed purebred lines and commercial broiler by Abasht et al. (8). Furthermore, these studies used different metabolomics-based methods [for example NMR in Wang's (20) vs. UPLC–HRMS in this study]. The effects of intrinsic sensitivity of detection and the capacity of separation of complex metabolite mixture by these different methods are not ruled out.

Ameliorating WB incidence by reducing its severity is of uppermost interest to the poultry industry because it would help in improving bird welfare, meat quality, and overall poultry production sustainability. As we previously showed that QB supplementation reduced the severity of WB myopathy (22), we sought to determine here its effect on breast muscle metabolomics profiling. Our data suggested that QB might improve WB via modulation of s-adenosyl-homocystein, arginine, tricarballylic acid, NAD+, glucosamine phosphate, citrulline, histamine, AMP, and pyridoxate. Although the exact underlying mechanisms are not known at this time, it is possible that QB modulates the above mentioned metabolites via improved solubility and digestibility of dietary nutrients and thereby enhanced release and bioavailability of calcium, phosphorus, minerals, and metal co-factors for endogenous enzymes and inositol liberation (81). For instance, it has been shown that QB superdosing improved liberation of zinc (82) and zinc has been reported to interact with ATP/ADP/AMP-hydrolyzing enzymes and have an antioxidant role via protection of sulfhydryl group against oxidation and inhibition of the ROS production by transition metals (83). Similarly, Tang and co-workers showed that phosphorus had antioxidant properties via enhancing superoxide dismutase and catalase activities in yellow cat fish (84). In chickens, it has been demonstrated that QB superdosing increased the levels of coenzyme Q(10), retinol, and alpha-tocopherol in liver (85). As coenzyme Q(10), retinol, and alpha-tocopherol all have been shown to improve mitochondrial respiration and function (86–88), our data suggested that QB might reduce the severity of WB via an improvement in the muscle oxidative status and metabolic profile.

In summary, we determined many metabolic biomarkers and disordered pathways, which could be regarded as new routes for discovering potential mechanisms of WB myopathy. Further in-depth investigations are warranted to define the mechanisms by which QB phytase modulates the metabolomics profile-ameliorating WB incidence.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics statement

The animal study was reviewed and approved by University of Arkansas.

Author contributions

SD conceived and designed the study. EG, RC, AD, and SD conducted the experiments and analyzed the data. SH, HC, and SC performed the UPLC-HRMS and generated the metabolomics data. MB provided the QB. SD wrote the paper with a critical review by MB, BK, and MK. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by a grant from ABV (to SD). ABV had no role in conducting the research, generating the data, interpreting the results or writing the manuscript.

Conflict of interest

MB is employed by company AB Vista. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2020.00458/full#supplementary-material

References

  • 1.

    SmithJSonesKGraceDMcMillanSTarawaliSHerreroM. Beyond milk, meat, and eggs: Role of livestock in food and nutrition security. Anim Front. (2013) 3:6–13. 10.2527/af.2013-0002

  • 2.

    MarangoniFCorselloGCricelliCFerraraNGhiselliALucchinLet al. Role of poultry meat in a balanced diet aimed at maintaining health and wellbeing: an Italian consensus document. Food Nutr Res. (2015) 59:27606. 10.3402/fnr.v59.27606

  • 3.

    ChatterjeeDZhuangHBowkerBCRinconAMSanchez-BrambilaG. Instrumental texture characteristics of broiler pectoralis major with the wooden breast condition. Poultry Sci. (2016) 95:2449–54. 10.3382/ps/pew204

  • 4.

    TallentireCWLeinonenIKyriazakisI. Artificial selection for improved energy efficiency is reaching its limits in broiler chickens. Sci Rep. (2018) 8:1168. 10.1038/s41598-018-19231-2

  • 5.

    MudalalSLorenziMSogliaFCavaniCPetracciM. Implications of white striping and wooden breast abnormalities on quality traits of raw and marinated chicken meat. Anim Int J Anim Biosci. (2015) 9:728–34. 10.1017/S175173111400295X

  • 6.

    TijareVVYangFLKuttappanVAAlvaradoCZCoonCNOwensCM. Meat quality of broiler breast fillets with white striping and woody breast muscle myopathies. Poultry Sci. (2016) 95:2167–73. 10.3382/ps/pew129

  • 7.

    MadrugaMSda RochaTCde CarvalhoLMSousaAde Sousa NetoACCoutinhoDGet al. The impaired quality of chicken affected by the wooden breast myopathy is counteracted in emulsion-type sausages. J Food Sci Tech. (2019) 56:1380–8. 10.1007/s13197-019-03612-0

  • 8.

    AbashtBMutrynMFMichalekRDLeeWR. Oxidative stress and metabolic perturbations in wooden breast disorder in chickens. PLoS ONE. (2016) 11:e0153750. 10.1371/journal.pone.0153750

  • 9.

    SihvoHKImmonenKPuolanneE. Myodegeneration with fibrosis and regeneration in the pectoralis major muscle of broilers. Vet Pathol. (2014) 51:61–23. 10.1177/0300985813497488

  • 10.

    AlnahhasNBerriCChabaultMChartrinPBoulayMBourinMCet al. Genetic parameters of white striping in relation to body weight, carcass composition, and meat quality traits in two broiler lines divergently selected for the ultimate pH of the pectoralis major muscle. BMC Genet. (2016) 17:61. 10.1186/s12863-016-0369-2

  • 11.

    RussoEDrigoMLongoniCPezzottiRFasoliPRecordatiC. Evaluation of white striping prevalence and predisposing factors in broilers at slaughter. Poultry Sci. (2015) 94:1843–8. 10.3382/ps/pev172

  • 12.

    HuangXAhnDU. The incidence of muscle abnormalities in broiler breast meat - a review. Korean J Food Sci Anim Resou. (2018) 38:835–50. 10.5851/kosfa.2018.e2

  • 13.

    CeminHSVieiraSLStefanelloCKindleinLFerreiraTZFiremanAK. Broiler responses to increasing selenium supplementation using Zn-L-selenomethionine with special attention to breast myopathies. Poultry Sci. (2018) 97:1832–40. 10.3382/ps/pey001

  • 14.

    de BrotSPerezSShivaprasadHLBaikerKPolledoLClarkMet al. Wooden breast lesions in broiler chickens in the UK.Vet Record. (2016) 178:141. 10.1136/vr.103561

  • 15.

    KawasakiTIwasakiTYamadaMYoshidaTWatanabeT. Rapid growth rate results in remarkably hardened breast in broilers during the middle stage of rearing: a biochemical and histopathological study. PLoS ONE. (2018) 13:e0193307. 10.1371/journal.pone.0193307

  • 16.

    KuttappanVAHargisBMOwensCM. White striping and woody breast myopathies in the modern poultry industry: a review. Poultry Sci. (2016) 95:2724–33. 10.3382/ps/pew216

  • 17.

    PapahMBBrannickEMSchmidtCJAbashtB. Evidence and role of phlebitis and lipid infiltration in the onset and pathogenesis of wooden breast disease in modern broiler chickens. Avian Pathol. (2017) 46:623–43. 10.1080/03079457.2017.1339346

  • 18.

    MutrynMFBrannickEMFuWLeeWRAbashtB. Characterization of a novel chicken muscle disorder through differential gene expression and pathway analysis using RNA-sequencing. BMC Genomics. (2015) 16:399. 10.1186/s12864-015-1623-0

  • 19.

    AbashtBZhouNLeeWRZhuoZPeripolliE. The metabolic characteristics of susceptibility to wooden breast disease in chickens with high feed efficiency. Poultry Sci. (2019) 98:3246–56. 10.3382/ps/pez183

  • 20.

    WangYYangYPanDHeJCaoJWangHet al. Metabolite profile based on (1)H NMR of broiler chicken breasts affected by wooden breast myodegeneration. Food Chem. (2020) 310:125852. 10.1016/j.foodchem.2019.125852

  • 21.

    KuttappanVABottjeWRamnathanRHartsonSDCoonCNKongBWet al. Proteomic analysis reveals changes in carbohydrate and protein metabolism associated with broiler breast myopathy. Poultry Sci. (2017) 96:2992–9. 10.3382/ps/pex069

  • 22.

    GreeneEFleesJDadgarSMallmannBOrlowskiSDhamadAet al. Quantum blue reduces the severity of woody breast myopathy via modulation of oxygen homeostasis-related genes in broiler chickens. Front Physiol. (2019) 10:1251. 10.3389/fphys.2019.01251

  • 23.

    ClemmonsBAMartinoCPowersJBCampagnaSRVoyBHDonohoeDRet al. Rumen bacteria and serum metabolites predictive of feed efficiency phenotypes in beef cattle. Sci Rep. (2019) 9:19265. 10.1038/s41598-019-55978-y

  • 24.

    LuWClasquinMFMelamudEAmador-NoguezDCaudyAARabinowitzJD. Metabolomic analysis via reversed-phase ion-pairing liquid chromatography coupled to a stand alone orbitrap mass spectrometer. Anal Chem. (2010) 82:3212–21. 10.1021/ac902837x

  • 25.

    KessnerDChambersMBurkeRAgusDMallickP. ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics. (2008) 24:2534–6. 10.1093/bioinformatics/btn323

  • 26.

    ChambersMCMacleanBBurkeRAmodeiDRudermanDLNeumannSet al. A cross-platform toolkit for mass spectrometry and proteomics. Nat Biotechnol. (2012) 30:918–20. 10.1038/nbt.2377

  • 27.

    ClasquinMFMelamudERabinowitzJD. LC-MS data processing with MAVEN: a metabolomic analysis and visualization engine. Curr Protocol Bioinform. (2012) 14:Unit14.11. 10.1002/0471250953.bi1411s37

  • 28.

    BazurtoJVDearthSPTagueEDCampagnaSRDownsDM. Untargeted metabolomics confirms and extends the understanding of the impact of Aminoimidazole Carboxamide Ribotide (AICAR) in the metabolic network of Salmonella enterica. Microb Cell. (2017) 5:74–87. 10.15698/mic2018.02.613

  • 29.

    WishartDSFeunangYDMarcuAGuoACLiangKVazquez-FresnoRet al. HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. (2018) 46:D608–17. 10.1093/nar/gkx1089

  • 30.

    de HoonMJImotoSNolanJMiyanoS. Open source clustering software. Bioinformatics. (2004) 20:1453–4. 10.1093/bioinformatics/bth078

  • 31.

    SaldanhaAJ. Java treeview–extensible visualization of microarray data. Bioinformatics. (2004) 20:3246–8. 10.1093/bioinformatics/bth349

  • 32.

    TredaCPopedaMKsiazkiewiczMGrzelaDPWalczakMPBanaszczykMet al. EGFR activation leads to cell death independent of PI3K/AKT/mTOR in an AD293 cell line. PLoS ONE. (2016) 11:e0155230. 10.1371/journal.pone.0155230

  • 33.

    JacksonNMCeresaBP. EGFR-mediated apoptosis via STAT3. Exp Cell Res. (2017) 356:93–103. 10.1016/j.yexcr.2017.04.016

  • 34.

    MansourAMAbdelrahimMLaymonMElsherbeenyMSultanMShokeirAet al. Epidermal growth factor expression as a predictor of chemotherapeutic resistance in muscle-invasive bladder cancer. BMC Urol. (2018) 18:100. 10.1186/s12894-018-0413-9

  • 35.

    BoyerJGPrasadVSongTLeeDFuXGrimesKMet al. ERK1/2 signaling induces skeletal muscle slow fiber-type switching and reduces muscular dystrophy disease severity. JCI Insight. (2019) 5:e127356. 10.1172/jci.insight.127356

  • 36.

    DridiSHiranoYTaralloVKimYFowlerBJAmbatiBKet al. ERK1/2 activation is a therapeutic target in age-related macular degeneration. Proc Natl Acad Sci USA. (2012) 109:13781–6. 10.1073/pnas.1206494109

  • 37.

    CaiKShaoWChenXCampbellYLNairMNSumanSPet al. Meat quality traits and proteome profile of woody broiler breast (pectoralis major) meat. Poultry Sci. (2018) 97:337–46. 10.3382/ps/pex284

  • 38.

    LivingstonMLFerketPRBrakeJLivingstonKA. Dietary amino acids under hypoxic conditions exacerbates muscle myopathies including wooden breast and white stripping. Poultry Sci. (2019) 98:1517–27. 10.3382/ps/pey463

  • 39.

    BodleBCAlvaradoCShirleyRBMercierYLeeJT. Evaluation of different dietary alterations in their ability to mitigate the incidence and severity of woody breast and white striping in commercial male broilers. Poultry Sci. (2018) 97:3298–310. 10.3382/ps/pey166

  • 40.

    WangTJLarsonMGVasanRSChengSRheeEPMcCabeEet al. Metabolite profiles and the risk of developing diabetes. Nat Med. (2011) 17:448–53. 10.1038/nm.2307

  • 41.

    LeeDSParkJKayKAChristakisNAOltvaiZNBarabasiAL. The implications of human metabolic network topology for disease comorbidity. Pro Natl Acad Sci USA. (2008) 105:9880–5. 10.1073/pnas.0802208105

  • 42.

    OrtmeyerHKBodkinNLHansenBC. Relationship of skeletal muscle glucose 6-phosphate to glucose disposal rate and glycogen synthase activity in insulin-resistant and non-insulin-dependent diabetic rhesus monkeys. Diabetologia. (1994) 37:127–33. 10.1007/s001250050082

  • 43.

    BaldiGYenCNDaughtryMRBodmerJBowkerBCZhuangHet al. Exploring the factors contributing to the high ultimate ph of broiler pectoralis major muscles affected by wooden breast condition. Front Physiol. (2020) 11:343. 10.3389/fphys.2020.00343

  • 44.

    AschenbachWGHirshmanMFFujiiNSakamotoKHowlettKFGoodyearLJ. Effect of AICAR treatment on glycogen metabolism in skeletal muscle. Diabetes. (2002) 51:567–73. 10.2337/diabetes.51.3.567

  • 45.

    SukhodubAJovanovicSDuQBudasGClellandAKShenMet al. AMP-activated protein kinase mediates preconditioning in cardiomyocytes by regulating activity and trafficking of sarcolemmal ATP-sensitive K(+) channels. J Cel Physiol. (2007) 210:224–36. 10.1002/jcp.20862

  • 46.

    AltamiranoFLopezJRHenriquezCMolinskiTAllenPDJaimovichE. Increased resting intracellular calcium modulates NF-kappaB-dependent inducible nitric-oxide synthase gene expression in dystrophic mdx skeletal myotubes. J Biol Chem. (2012) 287:20876–87. 10.1074/jbc.M112.344929

  • 47.

    DrakeJCAlwaySEHollanderJMWilliamsonDL. AICAR treatment for 14 days normalizes obesity-induced dysregulation of TORC1 signaling and translational capacity in fasted skeletal muscle. Am J Physiol Regul Integr Comp Physiol. (2010) 299:R1546–54. 10.1152/ajpregu.00337.2010

  • 48.

    Sasai-TakedatsuMTaketaniSNagataNFurukawaTTokunagaRKojimaTet al. Human trehalase: characterization, localization, and its increase in urine by renal proximal tubular damage. Nephron. (1996) 73:179–85. 10.1159/000189037

  • 49.

    HaddadGGMaE. Neuronal tolerance to O2 deprivation in drosophila: novel approaches using genetic models. Neuroscientist. (2001) 7:538–50. 10.1177/107385840100700610

  • 50.

    TanakaMMachidaYNiuSIkedaTJanaNRDoiHet al. Trehalose alleviates polyglutamine-mediated pathology in a mouse model of Huntington disease. Nat Med. (2004) 10:148–54. 10.1038/nm985

  • 51.

    DaviesJESarkarSRubinszteinDC. Trehalose reduces aggregate formation and delays pathology in a transgenic mouse model of oculopharyngeal muscular dystrophy. Hum Mol Genet. (2006) 15:23–31. 10.1093/hmg/ddi422

  • 52.

    ToutainPLBechuDHidiroglouM. Ascorbic acid disposition kinetics in the plasma and tissues of calves. Am J Physiol. (1997) 273:R1585–97. 10.1152/ajpregu.1997.273.5.R1585

  • 53.

    PeterneljTTCoombesJS. Antioxidant supplementation during exercise training: beneficial or detrimental?Sports Med. (2011) 41:1043–69. 10.2165/11594400-000000000-00000

  • 54.

    ChinRMFuXPaiMYVergnesLHwangHDengGet al. The metabolite alpha-ketoglutarate extends lifespan by inhibiting ATP synthase and TOR. Nature. (2014) 510:397–401. 10.1038/nature13264

  • 55.

    TokonamiNMorlaLCentenoGMordasiniDRamakrishnanSKNikolaevaSet al. alpha-Ketoglutarate regulates acid-base balance through an intrarenal paracrine mechanism. J Clin Invest. (2013) 123:3166–71. 10.1172/JCI67562

  • 56.

    LewisSSHutchinsonMRZhangYHundDKMaierSFRiceKCet al. Glucuronic acid and the ethanol metabolite ethyl-glucuronide cause toll-like receptor 4 activation and enhanced pain. Brain Behav Immun. (2013) 30:24–32. 10.1016/j.bbi.2013.01.005

  • 57.

    Ferreira GdaCViegasCMSchuckPFToninARibeiroCACoelho DdeMet al. Glutaric acid administration impairs energy metabolism in midbrain and skeletal muscle of young rats. Neurochem Res. (2005) 30:1123–31. 10.1007/s11064-005-7711-9

  • 58.

    ChengMLLinJFHuangCYLiGJShihLMChiuDTet al. Sedoheptulose-1,7-bisphospate accumulation and metabolic anomalies in hepatoma cells exposed to oxidative stress. Oxid Med Cell Long. (2019) 2019:5913635. 10.1155/2019/5913635

  • 59.

    KuehneAEmmertHSoehleJWinnefeldMFischerFWenckHet al. Acute activation of oxidative pentose phosphate pathway as first-line response to oxidative stress in human skin cells. Mol Cell. (2015) 59:359–71. 10.1016/j.molcel.2015.06.017

  • 60.

    TraversoNRicciarelliRNittiMMarengoBFurfaroALPronzatoMAet al. Role of glutathione in cancer progression and chemoresistance. Oxid Med Cell Long. (2013) 2013:972913. 10.1155/2013/972913

  • 61.

    CotgreaveIAGoldschmidtLTonkonogiMSvenssonM. Differentiation-specific alterations to glutathione synthesis in and hormonally stimulated release from human skeletal muscle cells. FASEB J. (2002) 16:435–7. 10.1096/fj.01-0685fje

  • 62.

    SuzukiTSuzukiTWadaTSaigoKWatanabeK. Taurine as a constituent of mitochondrial tRNAs: new insights into the functions of taurine and human mitochondrial diseases. EMBO J. (2002) 21:6581–9. 10.1093/emboj/cdf656

  • 63.

    PionPDKittlesonMDRogersQRMorrisJG. Myocardial failure in cats associated with low plasma taurine: a reversible cardiomyopathy. Science. (1987) 237:764–8. 10.1126/science.3616607

  • 64.

    BkailyGJaaloukDSaderSShbakloHPothierPJacquesDet al. Taurine indirectly increases [Ca]i by inducing Ca2+ influx through the Na(+)-Ca2+ exchanger. Mol Cell Biochem. (1998) 188:187–97. 10.1007/978-1-4615-5763-0_20

  • 65.

    LoehrerFMAngstCPBrunnerFPHaefeliWEFowlerB. Evidence for disturbed S-adenosylmethionine : S-adenosylhomocysteine ratio in patients with end-stage renal failure: a cause for disturbed methylation reactions?Nephrol Dial. (1998) 13:656–61. 10.1093/ndt/13.3.656

  • 66.

    KerinsDMKouryMJCapdevilaARanaSWagnerC. Plasma S-adenosylhomocysteine is a more sensitive indicator of cardiovascular disease than plasma homocysteine. Am J Clin Nutr. (2001) 74:723–9. 10.1093/ajcn/74.6.723

  • 67.

    JerniganHMJrDesoukyMAGellerAMBlumPSEkambaramMC. Efflux and hydrolysis of phosphorylethanolamine and phosphorylcholine in stressed cultured rat lenses. Exp Eye Res. (1993) 56:25–33. 10.1006/exer.1993.1005

  • 68.

    BleijerveldOBBrouwersJFVaandragerABHelmsJBHouwelingM. The CDP-ethanolamine pathway and phosphatidylserine decarboxylation generate different phosphatidylethanolamine molecular species. J Biol Chem. (2007) 282:28362–72. 10.1074/jbc.M703786200

  • 69.

    SelathuraiAKowalskiGMBurchMLSepulvedaPRisisSLee-YoungRSet al. The CDP-ethanolamine pathway regulates skeletal muscle diacylglycerol content and mitochondrial biogenesis without altering insulin sensitivity. Cell Metab. (2015) 21:718–30. 10.1016/j.cmet.2015.04.001

  • 70.

    FujiwaraTItohTKubotaYKuriyamaH. Effects of guanosine nucleotides on skinned smooth muscle tissue of the rabbit mesenteric artery. J Physiol. (1989) 408:535–47. 10.1113/jphysiol.1989.sp017474

  • 71.

    BatiukTDSchnizlein-BickCPlotkinZDagherPC. Guanine nucleosides and Jurkat cell death: roles of ATP depletion and accumulation of deoxyribonucleotides. Am J Physiol Cell Physiol. (2001) 281:C1776–84. 10.1152/ajpcell.2001.281.6.C1776

  • 72.

    MeulendijksDHenricksLMJacobsBAWAlievADeenenMJde VriesNet al. Pretreatment serum uracil concentration as a predictor of severe and fatal fluoropyrimidine-associated toxicity. Br J Cancer. (2017) 116:1415–24. 10.1038/bjc.2017.94

  • 73.

    DahlbackMHeintzLRyrfeldtAStenbergK. Toxic effects of some xanthine derivatives with special emphasis on adverse effects on rat testes. Toxicology. (1984) 32:23–35. 10.1016/0300-483X(84)90031-3

  • 74.

    CasaliEBerniPSpisniABaricchiRPertinhezTA. Hypoxanthine: a new paradigm to interpret the origin of transfusion toxicity. Blood Transfus. (2016) 14:555–6.11. 10.2450/2015.0177-15

  • 75.

    Cantu-MedellinNKelleyEE. Xanthine oxidoreductase-catalyzed reactive species generation: a process in critical need of reevaluation. Redox Biol. (2013) 1:353–8. 10.1016/j.redox.2013.05.002

  • 76.

    BattelliMGBolognesiAPolitoL. Pathophysiology of circulating xanthine oxidoreductase: new emerging roles for a multi-tasking enzyme. Biochim Biophys Acta. (2014) 1842:1502–17. 10.1016/j.bbadis.2014.05.022

  • 77.

    Sanchis-GomarFPareja-GaleanoHPerez-QuilisCSantos-LozanoAFiuza-LucesCGaratacheaNet al. Effects of allopurinol on exercise-induced muscle damage: new therapeutic approaches?Cell Stress Chaperon. (2015) 20:3–13. 10.1007/s12192-014-0543-2

  • 78.

    HipkissAR. NAD+ availability and proteotoxicity. Neuromol Med. (2009) 11:97–100. 10.1007/s12017-009-8069-y

  • 79.

    ThornalleyPJ. Protein and nucleotide damage by glyoxal and methylglyoxal in physiological systems–role in ageing and disease. Drug metabolism and drug interactions (2008) 23:125–50. 10.1515/DMDI.2008.23.1-2.125

  • 80.

    SogliaFSilvaAKLiaoLMLaghiLPetracciM. Effect of broiler breast abnormality and freezing on meat quality and metabolites assessed by 1H-NMR spectroscopy. Poultry Sci. (2019) 98:7139–50. 10.3382/ps/pez514

  • 81.

    TamimNMAngelRChristmanM. Influence of dietary calcium and phytase on phytate phosphorus hydrolysis in broiler chickens. Poultry Sci. (2004) 83:1358–67. 10.1093/ps/83.8.1358

  • 82.

    MohannaCNysY. Changes in zinc and manganese availability in broiler chicks induced by vegetal and microbial phytases. Anim Feed Sci Tech. (1999) 77:241–53. 10.1016/S0377-8401(98)00254-5

  • 83.

    BrayTMBettgerWJ. The physiological role of zinc as an antioxidant. Free Radic Biol Med. (1990) 8:281–91. 10.1016/0891-5849(90)90076-U

  • 84.

    TangQWangCXieCJinJHuangY. Dietary available phosphorus affected growth performance, body composition, and hepatic antioxidant property of juvenile yellow catfish Pelteobagrus fulvidraco. Sci World J. (2012) 2012:987570. 10.1100/2012/987570

  • 85.

    KaradasFPirgozlievVPappasACAcamovicTBedfordMR. Effects of different dietary phytase activities on the concentration of antioxidants in the liver of growing broilers. J Anim Physiol Anim Nutr. (2010) 94:519–26. 10.1111/j.1439-0396.2009.00938.x

  • 86.

    BarbiroliBFrassinetiCMartinelliPIottiSLodiRCortelliPet al. Coenzyme Q10 improves mitochondrial respiration in patients with mitochondrial cytopathies. An in vivo study on brain and skeletal muscle by phosphorous magnetic resonance spectroscopy. Cell Mol Biol. (1997) 43:741–9.

  • 87.

    TripathySChapmanJDHanCYHogarthCAArnoldSLOnkenJet al. All-trans-retinoic acid enhances mitochondrial function in models of human liver. Mol Pharmacol. (2016) 89:560–74. 10.1124/mol.116.103697

  • 88.

    FengZLiuZLiXJiaHSunLTianCet al. alpha-Tocopherol is an effective phase II enzyme inducer: protective effects on acrolein-induced oxidative stress and mitochondrial dysfunction in human retinal pigment epithelial cells. J Nutr Biochem. (2010) 21:1222–31. 10.1016/j.jnutbio.2009.10.010

Summary

Keywords

woody breast, metabolomics, broilers, quantum blue, IPA

Citation

Greene E, Cauble R, Dhamad AE, Kidd MT, Kong B, Howard SM, Castro HF, Campagna SR, Bedford M and Dridi S (2020) Muscle Metabolome Profiles in Woody Breast-(un)Affected Broilers: Effects of Quantum Blue Phytase-Enriched Diet. Front. Vet. Sci. 7:458. doi: 10.3389/fvets.2020.00458

Received

30 April 2020

Accepted

22 June 2020

Published

04 August 2020

Volume

7 - 2020

Edited by

Vincenzo Tufarelli, University of Bari Aldo Moro, Italy

Reviewed by

Massimiliano Petracci, University of Bologna, Italy; Michael Babak Papah, University of Pennsylvania, United States

Updates

Copyright

*Correspondence: Sami Dridi

This article was submitted to Animal Nutrition and Metabolism, a section of the journal Frontiers in Veterinary Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics