- 1Harvard Medical School, Harvard University, Boston, MA, United States
- 2Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA, United States
- 3Department of Metabolism, Digestion and Reproduction, Faculty of Medicine, Imperial College London, London, United Kingdom
- 4Australian National Phenome Centre, Computational and Systems Medicine, Health Futures Institute, Perth, WA, Australia
- 5School of Biomedical Sciences, The University of Western Australia, Perth, WA, Australia
- 6Perron Institute for Neurological and Translational Science, Nedlands, WA, Australia
- 7Centre for Computational and Systems Medicine, Health Futures Institute, Murdoch University, Perth, WA, Australia
- 8Translational Genetics, Telethon Kids Institute, Perth, WA, Australia
- 9WA Department of Health, Burns Service of Western Australia, Perth, WA, Australia
- 10Faculty of Medicine, Institute of Global Health Innovation, London, United Kingdom
Introduction: Burn injury in children causes prolonged systemic effects on physiology and metabolism leading to increased morbidity and mortality, yet much remains undefined regarding the metabolic trajectory towards specific health outcomes.
Methods: A multi-platform strategy was implemented to evaluate the long-term immuno-metabolic consequences of burn injury combining metabolite, lipoprotein, and cytokine panels. Plasma samples from 36 children aged 4–8 years were collected 3 years after a burn injury together with 21 samples from non-injured age and sex matched controls. Three different 1H Nuclear Magnetic Resonance spectroscopic experiments were applied to capture information on plasma low molecular weight metabolites, lipoproteins, and α-1-acid glycoprotein.
Results: Burn injury was characterized by underlying signatures of hyperglycaemia, hypermetabolism and inflammation, suggesting disruption of multiple pathways relating to glycolysis, tricarboxylic acid cycle, amino acid metabolism and the urea cycle. In addition, very low-density lipoprotein sub-components were significantly reduced in participants with burn injury whereas small-dense low density lipoprotein particles were significantly elevated in the burn injured patient plasma compared to uninjured controls, potentially indicative of modified cardiometabolic risk after a burn. Weighted-node Metabolite Correlation Network Analysis was restricted to the significantly differential features (q <0.05) between the children with and without burn injury and demonstrated a striking disparity in the number of statistical correlations between cytokines, lipoproteins, and small molecular metabolites in the injured groups, with increased correlations between these groups.
Discussion: These findings suggest a ‘metabolic memory’ of burn defined by a signature of interlinked and perturbed immune and metabolic function. Burn injury is associated with a series of adverse metabolic changes that persist chronically and are independent of burn severity and this study demonstrates increased risk of cardiovascular disease in the long-term. These findings highlight a crucial need for improved longer term monitoring of cardiometabolic health in a vulnerable population of children that have undergone burn injury.
1. Introduction
Non-severe burn injury, defined as less than 20% total body surface area (TBSA), is a common childhood medical insult, accounting for over 85% of pediatric burn cases in Australia (1). Recent research has demonstrated that even non-severe childhood burns are associated with a significant increase in hospitalization for a range of morbidities, including cardiovascular diseases, suggestive of a sustained systemic response to this type of injury (2). However, non-severe burn injury has traditionally been less well-studied than severe burns (2). Severe burn injury induces a global multisystemic effect, far beyond the initial site of injury, and the subsequent hypermetabolic response has been shown to persist for several years post-injury, even when no scarring remains (3). However, the extent and nature of this response in non-severe burns is less well established.
The burn-induced hypermetabolic response is thought to result from a physiological drive toward healing, which causes elevated catabolism (4), inflammation and sustained stress hormone release (5). The hypermetabolic state has been found to be persistent in many patients with severe burns, with its resolution linked to the longevity and duration of exposure to elevated levels of catecholamines, cortisol and glucagons (6). Serum cortisol measurements, and a range of cytokines including IL-6 and TNF-α have been reported to remain significantly increased for up to 3 years post-burn injury in children (7). This heightened proinflammatory state is typically accompanied by hyperglycemia and increased insulin resistance, independent of burn severity (8). Post-burn dysregulation of metabolism and immune regulation has also been associated with increased production of reactive oxygen species (ROS) and reactive nitrogen species such as nitric oxide and peroxynitrite (9). While severe burns have been well documented to induce this chronic inflammatory and hypermetabolic response, the mechanisms behind these chronic aberrations and whether these changes persist after non-severe burns, is less clear (10, 11).
One of the hallmarks of burn injury is the high rate of hospital admissions for cardiovascular disease (CVD) compared to the general population within the first few years post-injury (12). However, some researchers report no such elevation in CVD risk (13). Our prior work, deploying targeted metabolic phenotyping strategies to map the hypermetabolic state present in children long after recovery from a non-severe burn injury, demonstrated increased plasma levels of multiple amino acids and the neurotoxic metabolite quinolinic acid (14) but did not assess whether there are sustained metabolic changes that might be associated with the observed increased risk of cardiovascular disease (15). In order to address this substantial research gap, here we utilize 1H Nuclear Magnetic Resonance (NMR) spectroscopy to measure a wide panel of plasma metabolites. We apply three different pulse sequences to capture different molecular panels including lipoproteins (quantified based on the standard 1D NMR pulse sequence); small molecules (spin echo pulse sequence to attenuate contributions of macromolecules) and the JEDI pulse sequence (combining three spectroscopic editing techniques: diffusion; relaxation; J-editing to suppress small molecules and large proteins/lipoproteins) to capture biomarkers of inflammation. Together these pulse sequences enabled assessment of the metabolic landscape associated with the risk and incidence of cardiovascular disease (16). We identify a systemic signature of dysregulated lipoproteins associated with prior burn injury that may be linked to altered CVD risk post burn-injury.
2. Results
2.1. Metabolic phenotypes indicate persistent metabolic effects of burn injury in children
Differentiation of children with prior burn injury versus no burn injury was apparent for all three data sets focusing on lipoproteins, small molecules, and the inflammatory NMR panel (17). Orthogonal Partial Least Squares-Discriminant Analysis (OPLS-DA) is a robust statistical strategy to interrogate features in high-dimensional data, utilizing supervised regression analysis to define variation for class discrimination (18). The OPLS-DA models of all three sets of spectral data indicated there were significant metabolic differences between the prior burn injury and the control groups. The strongest classifier of burn injury was based on the standard 1D NMR pulse sequence dataset suggesting that burn injury was defined by a broad metabolic response as the standard profile contains contributions from low molecular weight molecules, lipoproteins, and glycoproteins. The cross-validated OPLS-DA model yielded a predictive value (Q2Y) of 0.52 and was predominantly driven by higher relative concentrations of glucose, amino acids (with the exception of glutamate and glycine), pyruvate, creatine, N-acetyl signals from α-1-acid glycoproteins (GlycA and GlycB) representing the burn injury group, with lower concentrations of acetate, acetoacetate, glutamate, and 3-methyl histidine (Figures 1A,B). The contribution of these low molecular weight metabolites differentiating the prior burn injury group from controls was further assessed using metabolite concentrations quantified using B.I.-LISA™ software (19). Glutamine, alanine, creatine, phenylalanine, and glucose were found to be directly correlated with burn injury, based on Student’s t test, while controlling the false-discovery rate (p < 0.05) and the Cliff’s delta value (Supplementary Figure S2B), whereas glutamic acid and 3-D-hydroxybutyric acid were inversely associated (Table 1). 3-D-hydroxybutyric acid was higher in the control group plasma samples but did not retain significance after controlling the false discovery rate (Figures 1C,D). There was a trend toward higher levels of other ketone bodies (acetone and acetoacetate) in the control group, but these were not statistically significant.
Figure 1. OPLS-DA scores plots, with torth (orthogonal) representing the greatest amount of systematic variation in the data set uncorrelated to Y (burn injury or healthy control groups), and tpred (predictive) representing the greatest amount of correlated variation to Y. Corresponding model loadings coefficients are shown, differentiating the NMR plasma spectra from children with prior burn injury (red) from the control group (blue) based on the standard 1D NMR pulse sequence (A,B); the CPMG pulse sequence (C,D) and the JEDI-PGPE pulse sequence (E,F). Key: 1- Lipid: -CH2CH3, 2-Isoleucine, 3-Leucine, 4-Valine, 5- Lipid: -CH2CH2CH2–6-Alanine, 7- Lipid: -CH2CH2CO 8-Lysine 9-Glutamate, 10-Methioine, 11-Glutamine, 12-Acetone, 13-Pyruvate, 14-Tyrosine, 15-Creatine, 16-Glucose, 17-Glycine, 18-3-methylhistidine, 19-Phenylalanine, D1-Lipid: -CH2CH3 D2- Lipid: -CH2CH2CH2- D3- Lipid: -CH2CH2CO D4- Lipid: -CH2CH=CH- D5- Lipid: -CH2CO D6- Lipid: -CH=CH-CH2-CH=CH- D7- SPC: +NCH2CH2OP- D8-GlycA/GlycB N-acetyl sugar moieties, D9- SPC:-+NCH2CH2OP- D10- Lipid: -CH=CH-.
Whereas the standard 1-Dimensional (1D) NMR experiment reflects the global molecular signature and composition, the Carr-Purcell Meiboom-Gill (CPMG) spin-echo experiment, effects spectral editing based on the proton T2 relaxation properties of proton signals by attenuating peaks from large molecules with slow rotational correlation times such as lipoproteins (20). In contrast JEDI-PGPE pulse sequences use pulsed field gradients to enhance signals from molecules with high segmental motional freedom but slow translational diffusion molecules that are contained within certain lipoprotein sub-compartments (21). Since the OPLS-DA model constructed from the standard 1D NMR spectra (Q2Y = 0.52) was stronger than either the CPMG (Q2Y 0.46; Figures 1C,D) or JEDI-PGPE (Q2Y 0.26; Figures 1E,F) models, this suggests that a broad range of low molecular weight molecules (most prevalent through 1D NMR spectral resolution), lipoproteins and glycoproteins are disrupted in children who have suffered prior burn injury. The CPMG model loadings (Figure 1D) further highlighted the association of citrate and formate with thermal injury whereas GlycA and GlycB were the strongest discriminatory signals differentiating burn from control group in the JEDI-PGPE data set (Figure 1F). The signals arising from the supramolecular phospholipid composite (SPC) peak were reduced in the burn injury group, but were not significantly differential. The corresponding PCA scores plots are provided in Supplementary Figure S1.
Metabolite-metabolite correlations were explored for key metabolites characteristic of prior burn injury using STOCSY (22) (Supplementary Figure S3). To investigate the association between 3-D-hydroxybutyric acid and ketogenesis, a correlation analysis driven from the apex of the acetoacetate peak (Supplementary Figure S3A) demonstrated a strong association between the three ketone bodies and 3-D-hydroxybutyric acid, further suggesting that burn injury was associated with a systematic decrease in plasma ketone bodies. Glucose was directly associated with GlycA and HDL lipoproteins (Supplementary Figure S3B); pyruvate was weakly associated with alanine and lactate (Supplementary Figure S3C); creatine was positively associated with multiple amino acids (Supplementary Figure S3D) and glutamine and glutamate were strongly inversely correlated (Supplementary Figure S3E).
2.2. Prior burn injury affects the plasma lipoprotein composition
PCA modeling of lipoproteins, quantified according to the Bruker iVDR B.I.-LISA™ method (19), showed a cluster of burn injury patients in the lower left-hand quadrant of the scores plot (Figure 2A). The differential signature of the plasma lipoprotein parameters was further defined in the OPLS-DA model (R2X 0.437, R2Y 0.492, Q2Y 0.321) and the eruption plot based on the model loadings versus Cliff’s delta statistic (Figures 2B,D). Prior burn injury was characterized by relatively higher plasma concentrations of LDL subfractions 5 and 6, particularly particle number, total cholesterol, free cholesterol, triglycerides, and Apolipoprotein B. Total LDPN and LDAB were also elevated in the burn injury group although this association was not significant after adjusting for multiple testing. In contrast, multiple VLDL5 subfractions were present in lower concentrations in the prior burn injury group in comparison to the controls as were L1TG and L4TG (Figures 2C,D; Table 2). Although the group size was small, the prior burn injury group was subdivided into burns caused by scalds versus flame. LDL-1 and LDL-4 triglyceride fractions were associated with scalds whereas the L6PL subfraction was associated with flame burns (Supplementary Figure S4).
Figure 2. PCA and OPLS-DA models of quantified lipoprotein subfractions. (A) PCA scores plot of burn injury patients and healthy control profiles (R2Y 0.987, Q2Y 0.933); (B) Scores plot of the equivalent OPLS-DA model (R2X 0.437, R2Y 0.492, Q2Y 0.321); (C) Variable importance eruption plot of predictive lipoprotein parameters based on OPLS-DA component loadings (tpred) and univariate effect sizes (Cliff’s delta). Co-ordinates are colored by the FDR-adjusted value of p (<0.05) from the statistical group comparison; (D) Significant lipoprotein parameters ranked by Cliff’s delta. Significant lipoprotein parameters distinguishing burn injury source (scald vs. flame) from multivariate linear regression modeling can be found in Supplementary Figure S4. L6CH, Low-Density Lipoprotein- 6 Subclass Cholesterols; L6PL, Low-Density Lipoprotein- 6 Subclass Phospholipids; L6PN, Low-Density Lipoprotein- 6 Subclass Particle Number; L6AB, Low-Density Lipoprotein- 6 Subclass Apolipoprotein-B100; L6FC, Low-Density Lipoprotein- 6 Subclass Free Cholesterol; L6TG, Low-Density Lipoprotein- 6 Subclass Triglycerides; L5CH, Low-Density Lipoprotein- 5 Subclass Cholesterols; L5AB, Low-Density Lipoprotein- 5 Subclass Apolipoprotein-B100; L5PN, Low-Density Lipoprotein- 5 Subclass Particle Number; L5PL, Low-Density Lipoprotein- 5 Subclass Phospholipids; LDPN, Low-Density Lipoprotein- Particle Number; LDAB, Low-Density Lipoprotein- Apolipoprotein-B100; H2A2, High-Density Lipoprotein- 2 Subclass Apolipoprotein-A2; L5FC, Low-Density Lipoprotein- 5 Subclass Free Cholesterol; H4TG, High-Density Lipoprotein- 4 Subclass Triglycerides; L1TG, Low-Density Lipoprotein- 1 Subclass Triglycerides; V5PL, Very Low-Density Lipoprotein- 5 Subclass Phospholipids; V5CH, Very Low-Density Lipoprotein- 5 Subclass Cholesterols; L4TG, Low-Density Lipoprotein- 4 Subclass Triglycerides; V5FC, Very Low-Density Lipoprotein- 5 Subclass Free Cholesterols.
Table 2. Quantified plasma lipoproteins significantly differentiating children with prior burn injury from healthy controls (p < 0.05).
2.3. Correlation of and other inflammatory markers lipoproteins with cytokines
α-1-acid glycoprotein N-acetyl signals (GlycA and GlycB) have long been known to be associated with systemic inflammation (23), Spearman’s correlation of lipoprotein parameters and JEDI-PGPE integrals for GlycA, GlycB and SPC were calculated and presented two distinct profiles corresponding to non-burn controls (Figure 3A) and burn injury (Figure 3B). In the control group there were negative correlations between GlycA and IL-7, IL-10 and GM-CSF (Figure 4A), whereas GlycB was positively correlated with several cytokines including IFN-γ, IL-13, IL-17a, IL-13, IL-6 and TNF-α, signatures absent in burn injury children(Figure 4B). In contrast, in the burn injury group both GlycA and GlycB showed strong direct correlations for IL-2, IL-6, IL-10, IL-13, TNF-α and IFN-γ reflecting the hyperinflammatory consequences associated with previous burn injury (7). Of these cytokines, IL-2, TNF-α and IFN-γ showed independent association with burn injury.
Figure 3. Spearman’s correlation between lipoproteins and JEDI-PGPE quantified integrals for (A) non-burn controls and (B) burn injury patients; and between cytokines and differential lipoproteins (C) non-burn controls and (D) burn injury patients. GlycA, α-1-acid glycoprotein signal A; GlycB, α-1-acid glycoprotein signal B; SPC1, Supramolecular Phospholipid Composite Integral 1; SPC2, Supramolecular Phospholipid Composite Integral 2; SPC3, Supramolecular Phospholipid Composite Integral 3; TotalSPC, Total Supramolecular Phospholipid Composite; GM-CSF, Granulocyte-macrophage colony-stimulating factor; IFNγ, Interferon γ; IL-10, Interleukin 10, IL-12(p70), Interleukin 12 heterodimeric 70kDA; IL-13, Interleukin 13; IL-17A, Interleukin 17A; IL-1B, Interleukin 1B; IL-2, Interleukin 2; IL-5, Interleukin 5; IL-6, Interleukin 6; IL-7, Interleukin 7; IL-8, Interleukin 8; TNF-α, Tumour Necrosis Factor-α.
Figure 4. (A,B) Circular correlation map for healthy control (A) and prior burn injury (B) groups. (C–R) Weighted-node correlation networks showing inter-connectedness of significant parameters (p < 0.05) from the integrated three panels (cytokine, metabolites and lipoprotein) for control and burn injury groups highlighting key nodes. Additionally, IL-6 and IL-8 are included in the correlation maps based on their importance following burn injury. L6CH, Low-Density Lipoprotein- 6 Subclass Cholesterols; L6PL, Low-Density Lipoprotein- 6 Subclass Phospholipids; L6PN, Low-Density Lipoprotein- 6 Subclass Particle Number; L6AB, Low-Density Lipoprotein- 6 Subclass Apolipoprotein-B100; L6FC, Low-Density Lipoprotein- 6 Subclass Free Cholesterol; L6TG, Low-Density Lipoprotein- 6 Subclass Triglycerides; L5CH, Low-Density Lipoprotein- 5 Subclass Cholesterols; L5AB, Low-Density Lipoprotein- 5 Subclass Apolipoprotein-B100; L5PN, Low-Density Lipoprotein- 5 Subclass Particle Number; L5PL, Low-Density Lipoprotein- 5 Subclass Phospholipids; LDPN, Low-Density Lipoprotein- Particle Number; LDAB, Low-Density Lipoprotein- Apolipoprotein-B100; H2A2, High-Density Lipoprotein- 2 Subclass Apolipoprotein-A2; L5FC, Low-Density Lipoprotein- 5 Subclass Free Cholesterol; H4TG, High-Density Lipoprotein- 4 Subclass Triglycerides; L1TG, Low-Density Lipoprotein- 1 Subclass Triglycerides; V5PL, Very Low-Density Lipoprotein- 5 Subclass Phospholipids; V5CH, Very Low-Density Lipoprotein- 5 Subclass Cholesterols; L4TG, Low-Density Lipoprotein- 4 Subclass Triglycerides; V5FC, Very Low-Density Lipoprotein- 5 Subclass Free Cholesterols; GlycA, α-1-acid glycoprotein signal A; GlycB, α-1-acid glycoprotein signal B; IFNγ, Interferon γ; IL-2, Interleukin 2; IL-6, Interleukin 6; IL-7, Interleukin 7; IL-8, Interleukin 8; TNFα, Tumor Necrosis Factor α.
The same contrasting correlation profile for the control and prior burn (Figures 3C,D) groups was also evident in the correlation structure between the significant lipoprotein parameters and the cytokines. The cytokine correlation map for the full list of lipoproteins is provided in Supplementary Figure S5. The small dense LDL lipoprotein parameters from subclasses 5 and 6 were weakly anticorrelated with most measured cytokines in the control group, particularly GM-CSF, IL-7, IL-10, and IL-12, whereas the VLDL subfraction 5 parameters (V5PL, V5CH, V5FC) were weakly directly associated with most cytokines, particularly IL-2, IL-10 and IL-12. This pattern largely reversed in the prior burn injury group with a general positive correlation between cytokines and LDL-6 lipoprotein parameters with notable correlations with IL-7 and IL-12. VLDL5 sub-particles were inversely related with the global cytokine profile for the burn group. Plasma HDTG concentrations were lower in the burn group. Interestingly, H2TG, H3TG, H3A1 and H3A2 were positively correlated to GlycA in burn injury but showed no correlation in non-burn signatures.
The global correlation structure in these data can be visualized as circular correlation maps, illustrating Spearman’s correlations for the panels of metabolites, lipoproteins and cytokines that significantly differentiate the prior burn participants from the control group (Figure 4). Additionally, IL-6 and IL-8 were included in the network based on their previous reported sustained association with burn injury (7). A contrasting correlation pattern between cytokines and metabolites/lipoproteins was observed when the control group was compared to the prior burn injury group. Whereas the overall correlation between individual cytokines and metabolites/lipoproteins was slightly weaker in the prior burn injury group due to the presence of sub-phenotypes within this group (Supplementary Figure S6), the number of connections between nodes for the correlation networks was substantially greater for most metabolites and cytokines in the burn injury group, indicating a co-ordinated immune-driven systemic response to burn injury. The notable exceptions were H4TG and glucose, which demonstrated a greater number of connections to cytokines in the control group. Phenylalanine was the node with the greatest number of weighted correlations (n = 15), linked with LDL and VLDL parameters, TNF-α, IL-6, IL-8, glutamate, alanine, creatine and 3-D-hydroxybutyrate in the prior burn injury group, whereas for the control group, there was a lower number of significant associations (n = 6), losing the association with the LDL parameters, TNF-α and IL-8 correlation, despite a cluster of strong immuno-metabolic connections presented in the circular correlation map (Figure 4A). GlycB correlated with 9 parameters, including 3 cytokines in the burn injury group, while in the control group, GlycB only mapped to glutamine, alanine and creatine (all elevated following burn injury). Similarly, the small dense lipoprotein parameters (e.g., L6PN) had higher weighted node correlations with cytokines and alanine and phenylalanine in the burn group.
3. Discussion
3.1. Immuno-inflammatory signatures of burn injury
We identified an inflammatory signature in the plasma profile of children with prior burn injury based on elevated concentrations of GlycB. This is consistent with our previous observation of elevated cytokines TNF-α (1.31-fold), IL-2 (1.18-fold), IL-7 (1.63-fold), and IFN-γ indicating that children with prior burns manifested low grade chronic inflammation. GlycA, and GlycB, represent glycosylated amino sugars in sidechains of a composite of five main acute phase glycoproteins: α-1-acid glycoprotein, α-1-antichymotrypsin, α-1-antitrypsin, haptoglobin and transferrin. GlycA and GlycB are robust markers of inflammation, found to be superior in accuracy to traditional pathology measures such as C-Reactive Protein (CRP) (24). Elevations in these N-acetylated glycoproteins have been reported in association with non-ischemic heart failure (25), vascular aging (26), rheumatoid arthritis (27), SARS-CoV-2 infection (21), insulin resistance (28) and atherosclerosis (29). Glycosylation patterns of serum proteins are typically expected during acute phase injury. Connelly et al., found in children suffering acute Kawasaki disease, elevated GlycA, increased total and small LDL particle numbers and decreased HDL particle numbers (p < 0.001) (30). Similar signatures were observed in this cohort of children, with elevated GlycA and GlycB and elevated small dense LDL particle numbers (mainly L6PN). GlycB also demonstrated significant correlations with TNF-α, IFN-γ, IL-2, IL-6, IL-7 and IL-8 exclusively in burn injury children, suggesting that it is associated with the systemic response to inflammation post burn injury (due to overlap with lipoprotein peaks, the higher concentrations of GlycA in the prior burn injury group did not reach statistical significance). Burn-induced increase in TNF-α concentrations can increase the rate of lipolysis and hence production of free fatty acids (31), which can contribute to the hyperglycemia and increased insulin resistance associated with post burn injury (32).
Johnston et al., demonstrated elevated levels of TNF-α (1.31-fold), IL-2 (1.18-fold), IL-7 (1.63-fold), and IFN-γ (1.18-fold) 3 years post injury in the pediatric burn injury patients in comparison to healthy controls and described a significantly reduced antibody response to diphtheria, tetanus, and pertussis vaccine antigens (7). Numerous studies have reported high levels of serum cytokines and chemokines post burn injury including IL-7, IL-10, IL-12, macrophage inflammatory protein-1b with some studies noting the sustained elevation of these cytokines for at least 3 years post injury (33). Sustained elevation of cytokines such as TNF-α, IL-2 and IL-17 at 3 years post burn may suggest a transition from a Th1 response in acute burn injury toward a Th2 dominated immune response in the chronic phase (6).
3.2. The metabolic ‘memory’ of prior burn injury is associated with distinctive lipoprotein signatures
Multiple plasma lipoproteins were systematically altered in participants with prior burn injury. Specifically, we found higher concentrations of small dense LDL subfractions (LDL5 and LDL6) to be associated with burn injury. These small dense LDL particles are known to have greater atherogenic potential than larger LDL subfractions (LDL1-4) and in particular, small dense LDL cholesterol is more predictive of cardiovascular disease than that of total LDL-C (34). A large population-based study also found that increased long-term risk of ischemic heart disease in men was related to preferential accumulation of small dense LDL particles but not large LDL particles (35). This raises the question as to whether the association of LDL5 and LDL6 cholesterol and particle number could have implications with respect to increased risk of cardiovascular disease following burn injury. In one study, higher LDL particle number was consistently associated with increased risk for cardiovascular disease, independent of other lipid measurements, but other LDL parameters were not associated with cardiovascular disease after adjustment for cholesterol concentrations (36). There is substantial evidence demonstrating associations between burn injury and cardiometabolic diseases (37). Hypertriglyceridemia, hypocholesterolemia, and hypophospholipidemia have been reported to occur following thermal injury (38) with excessive accumulation of triacylglycerols in liver tissue (39). Concentrations of both low- and high-density lipoprotein have been reported to be reduced following burn injury, which is consistent with depleted blood cholesterol and phospholipid levels that occur following burn injury (40). While the altered lipid profiles may certainly contribute the increased risk of cardiometabolic disease, the systemic effect of burn injury is complex and it is likely that immune dysfunction and inflammatory mediators, also modulate cardiac damage (41).
3.3. Burn injury causes persistent changes in metabolite profiles reflecting hypermetabolism and hyperglycemia
Acute burn trauma is associated with long-term increased risk of developing multiple morbidities and with greater all-cause mortality (7). Burn injury exerts a profound physiological and metabolic impact resulting in raised resting energy expenditure, inflammation, altered cardiac and organ function and hypermetabolism and these effects persist for several years following the acute injury (6, 42). In the current study, we applied different NMR-based assays and integrated the metabolic response with cytokine data. All three NMR datasets indicated that prior burn injury, sustained 3 years prior to the study, induced a prolonged impact on the plasma phenotype with systematic differences in comparison with the non-injured group. These results are consistent with the known incidence of hypermetabolism and hyperglycemia, which are initiated a few days following injury and continue for several years after the burn event (43, 44), but which are not necessarily associated with the degree or severity of the burn.
We found a trend toward lower mean plasma concentrations of glucose for the control group in comparison with the burn injury group, which contributed to the differential weighting in the CPMG NMR data set (low molecular weight molecules) but was not significant for the univariate comparison using the quantified metabolite data. Burn injury stimulates insulin production and produces insulin resistance in liver, skeletal muscle, and adipose tissue, associated with post-receptor alterations such as phosphorylation of the insulin receptor substrate-1 (IRS-1) in the absence of changes in insulin receptor binding (45). Hyperglycemia during the acute phase of burn injury has been shown to be predictive of clinical outcome with high glucose levels being associated with higher incidence of infection, sepsis, and mortality (46). Similarly, abnormal insulin sensitivity indicative of peripheral and whole-body insulin resistance has been noted.
In this study, we find distinctly higher concentrations of branched chain amino acids (BCAAs) including valine, leucine, isoleucine and aromatic amino acids (AAAs: phenylalanine, tyrosine) in the prior burn injury group. Both branched chain amino acids (BCAAs) and aromatic amino acids (AAAs) have been associated with increased risk of developing type 2 diabetes with various mechanisms being proposed for their role in development of type 2 diabetes (47, 48). Leucine and other amino acids can induce pancreatic beta cells and activate the mammalian target of rapamycin (mTOR) mitogenic signaling pathway, which results in early beta cell dysfunction (49, 50). However, BCAAs can also enhance glycogen production from non-insulin dependent pathways via protein kinase C or phosphatidylinositol 3-kinase (51). The upregulation of BCAA and AAA production is consistent with the known 2.21-fold higher increase in developing type 2 diabetes within 5 years post burn injury (52). The subtle elevation of glucose and BCAAs along with the increase in serum phenylalanine levels may be indicative of an underlying impact on insulin sensitivity in the prior burn group.
One of the strongest changes in the burn cohort was the elevation of glutamine and depletion of glutamate, in comparison to healthy controls, with a 2.5-fold higher glutamine:glutamate ratio. In our previous targeted mass spectrometric study, we found similar evidence of hypermetabolism with significantly higher concentrations of 11 amino acids and quinolinic acid. Several of the amino acids overlapped with those measured by 1H NMR spectroscopy including alanine, phenylalanine, with several other amino acids demonstrating a similar trend, although not statistically significant (14). In contrast, with our observation of a higher glutamine: glutamate ratio in the prior burn injury group, most of the literature characterizing the metabolic phenotype of metabolic dysregulation [e.g., in relation to obesity (53), type-2 diabetes (54, 55) and recovery from viral infections such as long-COVID (56)] points to higher concentrations of glutamate and a lower glutamine: glutamate ratio driving the metabolic dysfunction. Similarly in incubated enterocytes obtained from burn-injury participants, increased utilization of glutamine, has been shown to lead to the increased formation of glutamate and alanine (57). The difference between the results found in the current study and other literature is difficult to reconcile and is complicated by the fact that most of the literature studies describe adult populations. Therefore the discrepancy in glutamine: glutamate ratio may be influenced by the fact that glutamine and glutamate production differs between adults and children (58).
The generalized hypermetabolic response associated with burn injury is not dependent on the size of burn and the responses elicited by a burn covering 30% or 50% of the total body surface area are the same (48). Jeschke et al. showed that unlike the hyperinflammatory response that persists for around 6 weeks post-injury, the hypermetabolic state and hyperglycemia can persist for several years (5). Here we show that prior burn injury is characterized by higher plasma glucose concentrations for at least 3 years after non-severe burns. Statistical correlation driven from the apex of the α-anomeric glucose signal at δ5.25 showed a weak correlation with GlycA/GlycB peaks (r < 0.25) indicating that the hyperglycemia may be associated with associated with inflammation since GlycA has been shown to reflect inflammation and is strongly correlated with high sensitivity C-reactive protein (hsCRP) (49). In fact, GlycA has been suggested to be a more robust marker of inflammation than C-reactive protein (CRP) (19). Following burn injury, it has been shown that metabolism shifts toward increased glycolysis, glycogenolysis, gluconeogenesis, lipolysis and proteolysis resulting from a severe energy deprivation at the cellular level promoting a hypermetabolic state (50). Altered hormone levels (increased circulating cortisol, catecholamines and glucagon) following burn injury also contribute to the hypermetabolic state (4) and serve to increase the cycling of glucose and free fatty acids (36).
Hypermetabolism induced by thermal injury has been associated with a switch from white to beige or brown fat metabolism, which may be related to increased mitochondrial mass and uncoupling protein 1 expression (59). Consistent with disrupted mitochondrial metabolism, we find prior burn injury to be associated upregulation of the tricarboxylic acid cycle as indicated by significantly higher concentrations of citrate in the multivariate models. Elevated pyruvic acid concentrations in the prior burn injury group may also impact on altered mitochondrial metabolism and has previously been reported to be increased 7 days post burn (60). It has been suggested that following burn injury, pyruvate is diverted from the TCA cycle toward lactate and alanine synthesis with Valbuena et al. finding that 60% of pyruvate, lactate and alanine carbons was glucose-derived following burn injury (61). Here we see increased plasma pyruvate and alanine without a concomitant increase in lactate concentrations. Another theory, proposed by Wolfe et al., is that there is an increased rate of substrate recycling in burn patients resulting in enhanced glycolytic-gluconeogenic and triglyceride-fatty acid cycling, which could serve to provide more metabolic flexibility to adapt efficiently to changes in burn-induced energy-substrate demands (62). Part of this energy recycling involves an upregulation of the Cahill cycle, which serves the purpose of transporting toxic nitrogenous waste from the muscles to the liver by converting L-glutamate and pyruvate into α-ketoglutarate and L-alanine (63), facilitating the transport of alanine to the liver where alanine is oxidatively deaminated to form pyruvate. The levels of alanine and aspartate aminotransferases have been shown to increase post-acute burn injury. However, Chiarelli et al. found that these transaminase levels do not always drop immediately after acute injury (64).
Phenylalanine concentrations were significantly higher in the plasma of the burn injury group than controls. Higher levels of phenylalanine have been noted in association with burn injury (primarily in the acute phase) and the elevation was associated with the hypermetabolic state (65). In severe cases of burn, even though the total free amino acids in plasma dropped, plasma phenylalanine and the phenylalanine to tyrosine ratio was consistently higher and was strongly associated with death and weight loss in both animals and patients (66). With regard to the current study, in the burn group, but not the control, phenylalanine correlated positively with multiple cytokines and small low density lipoprotein particles but negatively with 3-D-hydroxybutyrate. Collectively, there was a reduction of plasma ketone bodies in the burn injury group, although this was only significant for 3-D-hydroxybutyrate (Supplementary Figure S3A). Abbott et al. observed that in the initial phase of recovery from burn injury, patients failed to mount a ketonemic response to starvation monitored in non-burn participants following starvation (67).
3.4. Study strengths and limitations
This study utilizes novel metabolic profiling strategies to interrogate cardiometabolic risk in an overlooked vulnerable population. Here, we have demonstrated the elevated cardiometabolic risk for young children who have previously suffered from a burn injury. While prior research has established that these individuals are disproportionately affected by the development of comorbidities later in life, few follow up studies such as this exist. Although this current study was statistically underpowered for assessing if there are metabolic differences between different sources of burn, we did observe a difference between scald and flame burns scalds versus flame in the LDL lipoprotein signature with the L6PL subfraction associating with flame burns and lower density triglyceride subfractions associating with scalds. In a mouse model of thermal injury, IFN-γ was reported to be higher in scald-burned mice than in sham- and flame-burned mice suggesting that the type of burn injury may exert a differential effect on the immune-metabolic phenotypes (68). Given that there is a different chemical background to flame versus scald burns, with flame burns inducing greater formation of new protein degradation products, it is feasible that the systemic metabolic response may differ and warrants further investigation in a larger cohort. The authors recognize the limitations and constraints of capturing metabolic profiles at a single timepoint, and thus follow-up studies will be necessary.
Conclusion
The integrated multi-modal signature of plasma samples was indicative of chronic hypermetabolic and hyperinflammatory processes persisting at least 3 years post burn injury, driven by an immuno-metabolic memory. These processes were evident in the metabolic, glycoprotein, lipoprotein, and cytokine profiles. The weighted linkages between the various molecular panels underscored the systemic response to inflammation suggesting that the body mounts a prolonged response to the initial trauma with co-ordinated response between the small dense lipoprotein particles and a Th2-dominated cytokine profile. These findings support prior observation of a chronic hypermetabolic state following acute injury, driven by an immuno-metabolic memory. Further investigation of this systems response may uncover new understanding of the association of burn injury with increased long-term morbidity and mortality.
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 at: https://zenodo.org/, accession number: 7344673.
Ethics statement
The studies involving human participants were reviewed and approved by Child and Adolescent Health Service WA (approval numbers: 2015219EP; 1111EP; 768EP). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
SB: conceptualization, data curation, formal analysis, data visualization and manuscript preparation. SL and BJ: investigation, methodology, and manuscript preparation. DH, LW, and NG: investigation and methodology. SHB: investigation, methodology, and resources. VF: funding acquisition and manuscript preparation. MF and FW: conceptualization and manuscript preparation. EH and JN: conceptualization, funding acquisition, formal analysis, and manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
MF is supported by the Stan Perron Center of Excellence in Childhood Burns and the Perth Children’s Hospital Foundation. MRFF for funding the Australian National Phenome Center for development of methods used in this research. We thank the UK Medical Research Council (MRC), the Imperial College London MRC Doctoral Training Program and the Imperial College London Dean’s Award for funding SB, the Department of Jobs, Tourism, Science and Innovation, Government of Western Australian Premier’s Fellowship for funding EH, and the ARC Laureate Fellowship funding for EH.
Acknowledgments
The authors thank the Spinnaker Health Research Foundation, WA, The McCusker Foundation, WA, The Western Australian State Government. They also thank Susan L. Prescott and Debra J. Palmer for the use of control pediatric samples.
Conflict of interest
The 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.
Publisher’s note
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1105163/full#supplementary-material
Abbreviations
TBSA, Total Body Surface Area; NMR, Nuclear Magnetic Resonance; CPMG, Carr-Purcell Meiboom-Gill; JEDI-PGPE, J-Edited DiFFusional Pulsed Gradient Echo Experiment; PCA, Principal Component Analysis; OPLS-DA, Orthogonal Partial Least Squares-Discriminant Analysis; SD, Standard Deviation; GlycA, α-1-acid glycoprotein signal A; GlycB, α-1-acid glycoprotein signal B; SPC, Supramolecular Phospholipid Composite; STOCSY, Statistical TOtal Correlation SpectroscopY.
References
1. Tracy, L , Rosenblum, S , and Gabbe, B . Burns registry of Australia and New Zealand (BRANZ) 2018/19 annual report. (2020). Available at: https://www.monash.edu/__data/assets/pdf_file/0009/2191752/BRANZ-annual-report18_19_FINAL.pdf.
2. Krishnamoorthy, V , Ramaiah, R , and Bhananker, SM . Pediatric burn injuries. Int J Crit Illn Inj Sci. (2012) 2:128. doi: 10.4103/2229-5151.100889
3. Williams, FN , Herndon, DN , and Jeschke, MG . The hypermetabolic response to Burn injury and interventions to modify this response. Clin Plast Surg. (2009) 36:583–96. doi: 10.1016/j.cps.2009.05.001
4. Kraft, R , Herndon, DN , Finnerty, CC , Shahrokhi, S , and Jeschke, MG . Occurrence of multiorgan dysfunction in pediatric burn patients: incidence and clinical outcome. Ann Surg. (2014) 259:381–7. doi: 10.1097/SLA.0b013e31828c4d04
5. Herndon, DN , and Tompkins, RG . Support of the metabolic response to burn injury. Lancet. (2004) 363:1895–902. doi: 10.1016/S0140-6736(04)16360-5
6. Jeschke, MG , Gauglitz, GG , Kulp, GA , Finnerty, CC , Williams, FN , Kraft, R, et al. Long-term persistance of the pathophysiologic response to severe burn injury. PLoS One. (2011) 6:e21245. doi: 10.1371/journal.pone.0021245
7. Johnson, BZ , McAlister, S , McGuire, HM , Palanivelu, V , Stevenson, A , Richmond, P, et al. Pediatric Burn survivors have Long-term immune dysfunction with diminished vaccine response. Front Immunol. (2020) 11:1481. doi: 10.3389/fimmu.2020.01481
8. Rehou, S , Mason, S , Burnett, M , and Jeschke, MG . Burned adults develop profound glucose intolerance. Crit Care Med. (2016) 44:1059–66. doi: 10.1097/CCM.0000000000001605
9. Weidinger, A , and Kozlov, AV . Biological activities of reactive oxygen and nitrogen species: oxidative stress versus signal transduction. Biomol Ther. (2015) 5:472–84. doi: 10.3390/biom5020472
10. Porter, C , Tompkins, RG , Finnerty, CC , Sidossis, LS , Suman, OE , and Herndon, DN . The metabolic stress response to burn trauma: current understanding and therapies. Lancet. (2016) 388:1417–26. doi: 10.1016/S0140-6736(16)31469-6
11. Auger, C , Samadi, O , and Jeschke, MG . The biochemical alterations underlying post-burn hypermetabolism. Biochim Biophys Acta (BBA) - Mol Basis Dis. (2017) 1863:2633–44. doi: 10.1016/j.bbadis.2017.02.019
12. Duke, JM , Randall, SM , Fear, MW , O’Halloran, E , Boyd, JH , Rea, S, et al. Long term cardiovascular impacts after burn and non-burn trauma: a comparative population-based study. Burns. (2017) 43:1662–72. doi: 10.1016/j.burns.2017.08.001
13. Leung, B , Younger, JF , Stockton, K , Muller, M , and Paratz, J . Cardiovascular risk profile in burn survivors. Burns. (2017) 43:1411–7. doi: 10.1016/j.burns.2017.07.010
14. Begum, S , Johnson, BZ , Morillon, AC , Yang, R , Bong, SH , Whiley, L, et al. Systemic long-term metabolic effects of acute non-severe paediatric burn injury. Sci Rep. (2022) 12:1–10. doi: 10.1038/s41598-022-16886-w
15. O’Halloran, E , Shah, A , Dembo, L , Hool, L , Viola, H , Grey, C, et al. The impact of non-severe burn injury on cardiac function and long-term cardiovascular pathology. Sci Rep. (2016) 6:34650. doi: 10.1038/srep34650
16. Kostara, CE , Papathanasiou, A , Psychogios, N , Cung, MT , Elisaf, MS , Goudevenos, J, et al. NMR-based lipidomic analysis of blood lipoproteins differentiates the progression of coronary heart disease. J Proteome Res. (2014) 13:2585–98. doi: 10.1021/pr500061n
17. Nitschke, P , Lodge, S , Kimhofer, T , Masuda, R , Bong, SH , Hall, D, et al. J-edited DIffusional proton nuclear magnetic resonance spectroscopic measurement of glycoprotein and supramolecular phospholipid biomarkers of inflammation in human serum. Anal Chem. (2022) 94:1333–41. doi: 10.1021/acs.analchem.1c04576
18. Bylesjö, M , Rantalainen, M , Cloarec, O , Nicholson, JK , Holmes, E , and Trygg, J . OPLS discriminant analysis: combining the strengths of PLS-DA and SIMCA classification. J Chemom. (2006) 20:341–51. doi: 10.1002/cem.1006
19. Jiménez, B , Holmes, E , Heude, C , Tolson, RF , Harvey, N , Lodge, SL, et al. Quantitative lipoprotein subclass and low molecular weight metabolite analysis in human serum and plasma by 1H NMR spectroscopy in a multilaboratory trial. Anal Chem. (2018) 90:11962–71. doi: 10.1021/acs.analchem.8b02412
20. Nicholson, JK , Buckingham, MJ , and Sadler, PJ . High resolution 1H n.m.r. studies of vertebrate blood and plasma. Biochem J. (1983) 211:605–15. doi: 10.1042/bj2110605
21. Holmes, E , Lodge, S , Hall, D , Schaefer, H , Spraul, M , Embade, N, et al. Diffusion and relaxation edited proton NMR spectroscopy of plasma reveals a high-fidelity supramolecular biomarker signature of SARS-CoV-2 infection. Anal Chem. (2021) 93:3976–86. doi: 10.1021/acs.analchem.0c04952
22. Cloarec, O , Dumas, ME , Craig, A , Barton, RH , Trygg, J , Hudson, J, et al. Statistical Total correlation spectroscopy: An exploratory approach for latent biomarker identification from metabolic 1H NMR data sets. Anal Chem. (2005) 77:1282–9. doi: 10.1021/ac048630x
23. Bell, JD , Brown, JCC , Nicholson, JK , and Sadler, PJ . Assignment of resonances for ‘acute-phase’ glycoproteins in high resolution proton NMR spectra of human blood plasma. FEBS Lett. (1987) 215:311–5. doi: 10.1016/0014-5793(87)80168-0
24. Ballout, RA , and Remaley, AT . GlycA: a new biomarker for systemic inflammation and cardiovascular disease (CVD) risk assessment. J Lab Precis Med. (2020) 5:17–7. doi: 10.21037/jlpm.2020.03.03
25. Cediel, G , Teis, A , Codina, P , Julve, J , Domingo, M , Santiago-Vacas, E, et al. GlycA and GlycB as inflammatory markers in chronic heart failure. Am J Cardiol. (2022) 181:79–86. doi: 10.1016/j.amjcard.2022.07.019
26. Llauradó, G , Amigó, N , Fuertes-Martín, R , Romero, A , Cano, A , Albert, L, et al. Measurement of serum N-Glycans in the assessment of early vascular aging (arterial stiffness) in adults with type 1 diabetes. Diabetes Care. (2022) 45:2430–8. doi: 10.2337/dc22-0331
27. Fuertes-Martín, R , Taverner, D , Vallvé, JC , Paredes, S , Masana, L , Correig Blanchar, X, et al. Characterization of 1H NMR plasma glycoproteins as a new strategy to identify inflammatory patterns in rheumatoid arthritis. J Proteome Res. (2018) 17:3730–9. doi: 10.1021/acs.jproteome.8b00411
28. Serés-Noriega, T , Giménez, M , Perea, V , Blanco, J , Vinagre, I , Pané, A, et al. Quantification of glycoproteins by nuclear magnetic resonance associated with preclinical carotid atherosclerosis in patients with type 1 diabetes. Nutr Metab Cardiovasc Dis. (2021) 31:2099–108. doi: 10.1016/j.numecd.2021.03.021
29. An, L , Liu, Q , Feng, H , Bai, X , Dang, Y , Li, C, et al. Increased glycoprotein acetylation is associated with high cardiac event rates: analysis using coronary computed tomography angiography. Anatol J Cardiol. (2018) 20:152–8. doi: 10.14744/AnatolJCardiol.2018.01058
30. Connelly, MA , Shimizu, C , Winegar, DA , Shalaurova, I , Pourfarzib, R , Otvos, JD, et al. Differences in GlycA and lipoprotein particle parameters may help distinguish acute Kawasaki disease from other febrile illnesses in children. BMC Pediatr. (2016) 16:151. doi: 10.1186/s12887-016-0688-5
31. Rydén, M , Dicker, A , van Harmelen, V , Hauner, H , Brunnberg, M , Perbeck, L, et al. Mapping of early signaling events in tumor necrosis factor-alpha -mediated lipolysis in human fat cells. J Biol Chem. (2002) 277:1085–91. doi: 10.1074/jbc.M109498200
32. Mecott, GA , Al-Mousawi, AM , Gauglitz, GG , Herndon, DN , and Jeschke, MG . The role of hyperglycemia in burned patients: evidence-based studies. Shock. (2010) 33:5–13. doi: 10.1097/SHK.0b013e3181af0494
33. Sierawska, O , Małkowska, P , Taskin, C , Hrynkiewicz, R , Mertowska, P , Grywalska, E, et al. Innate immune system response to Burn damage—focus on cytokine alteration. Int J Mol Sci. (2022) 23:716. doi: 10.3390/ijms23020716
34. Ivanova, EA , Myasoedova, VA , Melnichenko, AA , Grechko, AV , and Orekhov, AN . Small dense low-density lipoprotein as biomarker for atherosclerotic diseases. Oxidative Med Cell Longev. (2017) 2017:1–10. doi: 10.1155/2017/1273042
35. St-Pierre, AC , Cantin, B , Dagenais, GR , Mauriège, P , Bernard, PM , Després, JP, et al. Low-density lipoprotein subfractions and the Long-term risk of ischemic heart disease in men. Arterioscler Thromb Vasc Biol. (2005) 25:553–9. doi: 10.1161/01.ATV.0000154144.73236.f4
36. Ip, S , Lichtenstein, AH , Chung, M , Lau, J , and Balk, EM . Systematic review: association of low-density lipoprotein subtractions with cardiovascular outcomes. Ann Intern Med. (2009) 150:474–84. doi: 10.7326/0003-4819-150-7-200904070-00007
37. Barrett, LW , Fear, VS , Waithman, JC , Wood, FM , and Fear, MW . Understanding acute burn injury as a chronic disease. Burns Trauma. (2019) 7:23. doi: 10.1186/s41038-019-0163-2
38. Coombes, EJ , Shakespeare, PG , and Batstone, GF . Lipoprotein changes after burn injury in man. J Trauma. (1980) 20:971–5. doi: 10.1097/00005373-198011000-00012
39. Mert, S , Bulutoglu, B , Chu, C , Dylewski, M , Lin, FM , Yu, YM, et al. Multiorgan metabolomics and Lipidomics provide new insights into fat infiltration in the liver, muscle wasting, and liver–muscle crosstalk following Burn injury. J Burn Care Res. (2021) 42:269–87. doi: 10.1093/jbcr/iraa145
40. Pratt, VC , Tredget, EE , Clandinin, MT , and Field, CJ . Fatty acid content of plasma lipids and erythrocyte phospholipids are altered following burn injury. Lipids. (2001) 36:675–82. doi: 10.1007/s11745-001-0772-y
41. Willis, MS , Carlson, DL , Dimaio, JM , White, MD , White, DJ , Adams, GA, et al. Macrophage migration inhibitory factor mediates late cardiac dysfunction after burn injury. Am J Physiol Heart Circ Physiol. (2005) 288:H795–804. doi: 10.1152/ajpheart.00189.2004
42. Elnaggar, RK , Osailan, AM , Mahmoud, WS , Alqahtani, BA , and Azab, AR . Beyond the acute phase: understanding relationships among cardiorespiratory response to exercises, physical activity levels, and quality of life in children after Burn injuries. J Burn Care Res. (2022) 43:827–33. doi: 10.1093/jbcr/irab203
43. Badoiu, SC , Miricescu, D , Stanescu-Spinu, II , Totan, AR , Badoiu, SE , Costagliola, M, et al. Glucose metabolism in Burns–what happens? Int J Mol Sci. (2021) 22:5159. doi: 10.3390/ijms22105159
44. Hart, DW , Wolf, SE , Mlcak, R , Chinkes, DL , Ramzy, PI , Obeng, MK, et al. Persistence of muscle catabolism after severe burn. Surgery. (2000) 128:312–9. doi: 10.1067/msy.2000.108059
45. Carter, EA , Burks, D , Fischman, AJ , White, M , and Tompkins, RG . Insulin resistance in thermally-injured rats is associated with post-receptor alterations in skeletal muscle, liver and adipose tissue. Int J Mol Med. (2004) 14:653–8.
46. Jeschke, MG , Kraft, R , Emdad, F , Kulp, GA , Williams, FN , and Herndon, DN . Glucose control in severely thermally injured pediatric patients: what glucose range should be the target? Ann Surg. (2010) 252:521–8. doi: 10.1097/SLA.0b013e3181f2774c
47. Gauglitz, GG , Herndon, DN , Kulp, GA , Meyer, WJ , and Jeschke, MG . Abnormal insulin sensitivity persists up to three years in pediatric patients post-burn. J Clin Endocrinol Metab. (2009) 94:1656–64. doi: 10.1210/jc.2008-1947
48. Zhang, L , Zhang, Y , Ma, Z , Zhu, Y , and Chen, Z . Altered amino acid metabolism between coronary heart disease patients with and without type 2 diabetes by quantitative 1H NMR based metabolomics. J Pharm Biomed Anal. (2021) 206:114381. doi: 10.1016/j.jpba.2021.114381
49. Krebs, M , Krssak, M , Bernroider, E , Anderwald, C , Brehm, A , Meyerspeer, M, et al. Mechanism of amino acid-induced skeletal muscle insulin resistance in humans. Diabetes. (2002) 51:599–605. doi: 10.2337/diabetes.51.3.599
50. Williams, FN , Herndon, DN , Suman, OE , Lee, JO , Norbury, WB , Branski, LK, et al. Changes in cardiac physiology AFTER severe BURN injury. J Burn Care Res. (2011) 32:269–74. doi: 10.1097/BCR.0b013e31820aafcf
51. Yoon, MS . The emerging role of branched-chain amino acids in insulin resistance and metabolism. Nutrients. (2016) 8:405. doi: 10.3390/nu8070405
52. Duke, JM , Randall, SM , Fear, MW , Boyd, JH , Rea, S , and Wood, FM . Diabetes mellitus after injury in burn and non-burned patients: a population based retrospective cohort study. Burns. (2018) 44:566–72. doi: 10.1016/j.burns.2017.10.019
53. Ren, W , Xia, Y , Chen, S , Wu, G , Bazer, FW , Zhou, B, et al. Glutamine metabolism in macrophages: a novel target for obesity/type 2 diabetes. Adv Nutr. (2019) 10:321–30. doi: 10.1093/advances/nmy084
54. Sikorski, C , Azab, S , de Souza, RJ , Shanmuganathan, M , Desai, D , Teo, K, et al. Serum metabolomic signatures of gestational diabetes in south Asian and white European women. BMJ Open Diabetes Res Care. (2022) 10:e002733. doi: 10.1136/bmjdrc-2021-002733
55. Elshorbagy, AK , Turner, C , Bastani, N , Refsum, H , and Kwok, T . The association of serum sulfur amino acids and related metabolites with incident diabetes: a prospective cohort study. Eur J Nutr. (2022) 61:3161–73. doi: 10.1007/s00394-022-02872-5
56. Holmes, E , Wist, J , Masuda, R , Lodge, S , Nitschke, P , Kimhofer, T, et al. Incomplete systemic recovery and metabolic Phenoreversion in post-acute-phase nonhospitalized COVID-19 patients: implications for assessment of post-acute COVID-19 syndrome. J Proteome Res. (2021) 20:3315–29. doi: 10.1021/acs.jproteome.1c00224
57. Ardawi, MSM , and Newsholme, EA . Maximal activities of glutaminase and some enzymes of glycolysis and ketone body utilization and rates of utilization of glutamine, glucose and ketone bodies by intestinal mucosa after burn injury. Burns. (1987) 13:438–44. doi: 10.1016/0305-4179(87)90220-8
58. Armstrong, MD , and Stave, U . A study of plasma free amino acid levels. II Normal values for children and adults. Metabolism. (1973) 22:561–9. doi: 10.1016/0026-0495(73)90069-3
59. Patsouris, D , Qi, P , Abdullahi, A , Stanojcic, M , Chen, P , Parousis, A, et al. Burn induces Browning of the subcutaneous White adipose tissue in mice and humans. Cell Rep. (2015) 13:1538–44. doi: 10.1016/j.celrep.2015.10.028
60. Aussel, C , Cynober, L , Lioret, N , Coudray-Lucas, C , Vaubourdolle, M , Saizy, R, et al. Plasma branched-chain keto acids in burn patients. Am J Clin Nutr. (1986) 44:825–31. doi: 10.1093/ajcn/44.6.825
61. Valbuena, GN , Rizzardini, M , Cimini, S , Siskos, AP , Bendotti, C , Cantoni, L, et al. Metabolomic analysis reveals increased aerobic glycolysis and amino acid deficit in a cellular model of amyotrophic lateral sclerosis. Mol Neurobiol. (2016) 53:2222–40. doi: 10.1007/s12035-015-9165-7
62. Wolfe, RR , Herndon, DN , Jahoor, F , Miyoshi, H , and Wolfe, M . Effect of severe Burn injury on substrate cycling by glucose and fatty acids. N Engl J Med. (2010) 317:403–8. doi: 10.1056/NEJM198708133170702
63. Karmen, A , Wroblewski, F , and Ladue, JS . Transaminase activity in human blood. J Clin Invest. (1955) 34:126–33. doi: 10.1172/JCI103055
64. Chiarelli, A , Casadei, A , Pornaro, E , Siliprandi, L , and Mazzoleni, F . Alanine and aspartate aminotransferase serum levels in burned patients: a long-term study. J Trauma. (1987) 27:790–4. doi: 10.1097/00005373-198707000-00017
65. Chang, XJ , Yang, CC , Hsu, WS , Xu, WZ , and Shih, TS . Serum and erythrocyte amino acid pattern: studies on major burn cases. Burns. (1983) 9:240–8. doi: 10.1016/0305-4179(83)90053-0
66. Stinnett, JD , Alexander, JW , Watanabe, C , BG, MM , Fischer, JE , Morris, MJ, et al. Plasma and skeletal muscle amino acids following severe burn injury in patients and experimental animals. Ann Surg. (1982) 195:75–89. doi: 10.1097/00000658-198201001-00012
67. Abbott, WC , Birkhahn, RH , Schiller, WR , Long, CL , and Blakemore, WS . The effect of major thermal injury on plasma ketone body levels. J Parenter Enter Nutr. (1985) 9:153–8. doi: 10.1177/0148607185009002153
68. Tschöp, J , Martignoni, A , Reid, MD , Adediran, SG , Gardner, J , Noel, GJ, et al. Differential immunological phenotypes are exhibited after scald and flame burns. Shock. (2009) 31:157–63. doi: 10.1097/SHK.0b013e31817fbf4d
69. Valvis, SM , Waithman, J , Wood, FM , Fear, MW , and Fear, VS . The immune response to skin trauma is dependent on the etiology of injury in a mouse model of Burn and excision. J Invest Dermatol. (2015) 135:2119–28. doi: 10.1038/jid.2015.123
70. Dona, AC , Jiménez, B , Schäfer, H , Humpfer, E , Spraul, M , Lewis, MR, et al. Precision high-throughput proton NMR spectroscopy of human urine, serum, and plasma for large-scale metabolic phenotyping. Anal Chem. (2014) 86:9887–94. doi: 10.1021/ac5025039
71. Nicholson, JK , Foxall, PJD , Spraul, M , Farrant, RD , and Lindon, JC . 750 MHz 1H and 1H-13C NMR spectroscopy of human blood plasma. Anal Chem. (1995) 67:793–811. doi: 10.1021/ac00101a004
72. Benjamini, Y , and Hochberg, Y . Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
73. Langfelder, P , and Horvath, S . WGCNA: An R package for weighted correlation network analysis. BMC Bioinformatics. (2008) 9:1–13. doi: 10.1186/1471-2105-9-559
Keywords: metabolomics, precision medicine, cardiometabolic disease risk, child health, burn injury, NMR, inflammation, cytokines
Citation: Begum S, Lodge S, Hall D, Johnson BZ, Bong SH, Whiley L, Gray N, Fear VS, Fear MW, Holmes E, Wood FM and Nicholson JK (2023) Cardiometabolic disease risk markers are increased following burn injury in children. Front. Public Health 11:1105163. doi: 10.3389/fpubh.2023.1105163
Edited by:
Shikai Yu, Tongji University, ChinaReviewed by:
Katrien Vandoorne, Technion Israel Institute of Technology, IsraelMathilde Sengoelge, Karolinska Institutet (KI), Sweden
Copyright © 2023 Begum, Lodge, Hall, Johnson, Bong, Whiley, Gray, Fear, Fear, Holmes, Wood and Nicholson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fiona M. Wood, fiona.wood@health.wa.gov.au; Jeremy K. Nicholson, j.nicholson@imperial.ac.uk