ORIGINAL RESEARCH article
The Systemic Alterations of Lipids, Alanine-Glucose Cycle and Inter-Organ Amino Acid Metabolism in Swine Model Confirms the Role of Liver in Early Phase of Septic Shock
- 1Department of Electronics, Information and Bioengineering, Politecnico di Milano, Milan, Italy
- 2Department of Environmental Health Sciences, Istituto di Ricerche Farmacologiche Mario Negri IRCCS, Milan, Italy
- 3Experimental Laboratory of Intensive Care, Erasme Hospital, Université Libre de Bruxelles, Brussels, Belgium
- 4Department of Intensive Care, Erasme Hospital, Université Libre de Bruxelles, Brussels, Belgium
Septic shock is a medical emergency and is one of the main causes of mortality in critically ill patients. Given the pathophysiological complexity of sepsis spectrum and progression in clinical settings, animal models become essential tools to improve patient care, and to understand key mechanisms that may remain masked from the heterogeneity of clinical practice. Our aim was to verify whether the metabolic constellations we previously reported for septic shock patients appear also in our septic shock swine model as systemic markers of early disturbances in energy metabolism and hepatic homeostasis. Septic shock was induced in anesthetized, instrumented, and ventilated adult swines by polymicrobial peritonitis. Hemodynamic and serial measurements of arterial and mixed venous blood gasses were made. Laboratory measurements and mass spectrometry-based targeted quantitative plasma metabolomics were performed in blood samples collected at baseline, at shock and at fully resuscitation after fluids and vasopressors administration. Data elaboration was performed by multilevel and multivariate analysis. Changes in hemodynamic, blood chemistry, and inflammatory markers were in line with a septic shock phenotype. Time course alteration of systemic metabolites were characterized by marked decreased in phosphatidylcholines and lysophosphatidylcholines species, altered alanine-glucose cycle and inter-organ amino acid metabolism, pointing toward an early hepatic impairment similarly to what we previously reported for septic shock. This is the first study in which an experimental swine model of septic shock recapitulates the main metabolic derangements reported in a clinical setting of shock. These events occur within hours from infections and may act as early metabolic features to assist in evaluating subclinical hepatic alterations and pave the way to improve the management of septic shock.
Septic shock is a life threatening condition that can develop subsequent to infection. It is a subset of sepsis with underlying circulatory and cellular/metabolic abnormalities associated with higher mortality rates (Angus et al., 2001; Singer et al., 2016). Septic shock causes progressive failure of vital homeostatic mechanisms culminating in immunosuppression, coagulopathy and microvascular dysfunction, which can lead to refractory hypotension, organ failure and death.
Early supportive therapy with fluid resuscitation and vasopressors to restore hemodynamics and reduce tissue hypoperfusion is decisive for the patient’s outcome and has figured in treatment guidelines for decades (Dellinger et al., 2013). However, mortality rate for septic shock may reach 40% even in the era of early recognition and treatment with today’s poor prognosis, mainly related to multi-organ dysfunction syndrome (MODS) (Rhee et al., 2017). The multiple organ involvement, which can lead to death within a few days, suggests that this concomitant process is systemic in nature and initiated in parallel with inflammatory response.
Given the pathophysiological complexity of sepsis spectrum and progression in clinical settings, animal models become essential tools to improve patient care, and to understand key mechanisms that may remain masked from the heterogeneity of clinical practice.
Several swine models of sepsis and septic shock have been developed to mimic the human disease progression and clinical settings by utilizing standard clinical treatments such as fluid resuscitation, antibiotics, and mechanical ventilation (Kubiak et al., 2011; Li et al., 2013; Revie et al., 2013; Corrêa et al., 2014; Chalkias et al., 2015; Vassal et al., 2015; Zhao et al., 2015; Duburcq et al., 2017). However, critiques on the relevance of experimental data regarding the translation to clinical sepsis/septic shock have been raised in the scientific community (Piper et al., 1996; Deitch, 1998; Esmon, 2004; Dyson and Singer, 2009; Opal and Patrozou, 2009). Aware of such criticism, our main goal in this study was to investigate whether our previous metabolic features characterizing septic shock patients (Ferrario et al., 2016; Cambiaghi et al., 2017) would have been observed also in a swine model of septic shock.
Previously, we found that lipids species alteration played an important role in individual patient responses to infection, where reduction in circulating lysophosphatidylcholines species, less change in plasmalogens concentration in combination with larger increment of alanine were associated with non-responsiveness to therapy (Cambiaghi et al., 2017). We speculated that alanine modulation indicated a possible alteration in the glucose-alanine cycle in the liver, providing a different picture of early liver functionality disturbances in the non-responder patients. Therefore, we aim to verify whether such metabolic constellations occurring in septic shock patients appear also in the septic shock swine model as systemic markers of early disturbances in energy metabolism and hepatic homeostasis, where liver acts not only as inflammatory barrier, but as organ actively involved in energy pathways.
In order to allow a robust comparison of the metabolomics data, we applied the same mass spectrometry-based target metabolomics strategy used for our previous clinical settings. Metabolomics analyses and biochemical/physiological measurements were conducted at three time points in septic shock and sham animals, i.e., at baseline, after shock development and after fully resuscitation, in order to capture changes in metabolic and physiologic trajectories.
Materials and Methods
We performed a controlled experimental study on a large animal model of septic shock induced by a polymicrobial peritonitis on adult swine in our Experimental Laboratory of Intensive Care (LA1230336), in the Université Libre de Bruxelles.
The local animal ethics committee (Comité Ethique du Bien-Être Animal) approved the present study (protocol 641 N) and we followed the EU Directive 2010/63/EU for animal experiments and the ARRIVE guidelines for animal research.
Nine pigs (Sus scrofa domesticus) of both sex, aged 4 to 6 months, with a mean (± SD) weight of 42,3 ± 3,9 kg, were obtained from a local farm (BE 400108–48). Animals were fasted for 18 h prior to the start of the experiment with free access to water. This is a closed chest and low invasive instrumentation animal model with all the catheters and introducers inserted percutaneously under echographic guidance.
The animals were sedated in their cage with an intramuscular injection of 1.5 mg/kg midazolam (Dormicum; Roche, Belgium) and 5 mg/kg azaperone (Stresnil, Eli Lilly Benelux, Belgium) in the neck. After transportation to the operating room, a 14 G peripheral venous line (14G, Surflo IV Catheter, Terumo Medical Company, Belgium) was placed into a vein of the ear to provide vascular access and a 4,5 Fr arterial catheter (Leader-Cath, Vygon, France) was placed in the left common femoral artery for invasive arterial pressure monitoring and blood samples collection. The animals underwent endotracheal intubation following induction of anesthesia with the intravenous injection of 3 μg/kg sufentanil (Sufenta Forte, Janssen-Cilag, Belgium), 1 mg/kg propofol (Propovet, Zoetis, Belgium) and 0.5 mg/kg of rocuronium (Esmeron, Organon, The Netherlands).
A central venous access for drugs infusion was obtained via a four lumen central venous line inserted inside the left femoral vein (Edwards Life Sciences, California, United States). General anesthesia and analgesia were achieved using continuous inhalation of sevoflurane 1,8 to 2,5 % MAC (Sevorane, Abbott, Belgium) and continuous infusion of sufentanil 1 to 4 μg/kg/h (adapted according to the response to painful stimulations, as a nasal septum pinching), in association with rocuronium continuous infusion to avoid shivering.
Mechanical ventilation was performed in volume-controlled mode (Primus, Draëger, Germany) with tidal volumes of 8 mL/kg and a PEEP set at 5 cm H2O. Arterial and venous blood gasses were analyzed live in the operating room (Cobas b-123, Roche, Switzerland) and the respiratory rate was adjusted to maintain PaCO2 values between 35 and 49 mmHg – aiming an arterial pH between 7,35 and 7,45 – while fraction of inspired oxygen (FiO2) was set at a value of 0,30 at baseline and adjusted to keep PaO2 > 90 mmHg. Fluid maintenance was provided by a balanced crystalloid (Plasmalyte, Baxter, Belgium) perfusion at 3 to 5 ml/kg/h, aiming a normal volemia at baseline, according to an arterial pulse pressure variation (PPV) ≤ 12% (and a negative fluid challenge in case of doubts, i.e., PPV remaining in the gray zone). Central temperature was maintained between 38,5 and 39,5°C with a warming blanket and hypoglycemia avoided by a continuous 10% glucose solution infusion (1 to 2 ml/kg/h).
After a supra-pubic mini-laparotomy, a 14 Fr Foley catheter (Beiersdorf AG, Germany) was introduced into the bladder to monitor urine output (UO) and two large bore abdominal drains were placed in each side of the abdominal cavity for later introduction of the feces. A 7 Fr introducer was inserted into the right external jugular vein and a pulmonary artery catheter (CCO; Edwards LifeSciences, California, United States) was advanced in a pulmonary artery for continuous cardiac output (CO), right heart pressures and mixed venous oxygen saturation monitoring. Both femoral artery pressure and pulmonary artery pressure signals were continuously displayed (SC9000, Siemens, Germany) and exported to an A/D recording station (Notocord Hem, Notocord, France), similarly to EKG and cardiac output signals.
After instrumentation, the animals were allowed to rest for approximately 2 h (premedication wash-out) after which the first baseline measurements and blood samples were taken (baseline, T1). The animals were then randomized for sepsis induction (shock group, n = 6) or a sham procedure (sham group, n = 3). Sepsis was induced by the intraperitoneal instillation, via the two abdominal drains, of 3 g/kg of autologous feces collected in the cage, filtered and diluted in 300 ml of glucose 10%, after which the drains were clamped. During septic shock onset, fluid maintenance was reduced to 1 ml/kg/h until the mean arterial pressure (MAP) has decreased below 50 mmHg owing to vasodilation and capillary leakage. Thereafter, fluid maintenance perfusion was moderately increased to avoid a further lethal decrease in blood pressure and to keep the animal alive for one more hour of severe hypotension (MAP goal between 45 and 50 mmHg), in order to consolidate the peripheral hypoperfusion and the multiple organ failure, after which the second time point T2, i.e., septic shock condition, was defined. Immediately after a series of hemodynamic measurements and blood samples, a full fluid resuscitation was initiated with both the same rate of the balanced crystalloid perfusion and an additional colloid perfusion (Geloplasma, Fresenius Kabi, France), aiming the same target as at baseline (PPV < 12% ± negative fluid challenge). The abdominal drains clamping was removed, to avoid fluid accumulation in the third space that otherwise would have led to an intra-abdominal hypertension. After 120 min of hemodynamic stabilization, defined by a stable MAP and no further increase in CO, additional hemodynamic measurements and blood samples were taken again (T3, first resuscitation). Finally, a vasopressor therapy was introduced, on top of the fluid resuscitation, with a continuous infusion of norepinephrine (Noradrenaline tartrate, Aguettant, Belgium) at a fixed dose of 0,3 μg/kg/min for one hour, after which the last series of hemodynamic measurements and blood samples were taken again (T4, full resuscitation). Animal were then euthanized with a potassium chloride injection and an overdose of thiopental.
In the septic shock group, two out of six animals were immediately fully resuscitated already with vasopressor because of the severity of the distributive shock, therefore they skipped the measurements at time point T3. For this reason, in the following we refer only to baseline (T1), shock (T2) and to fully resuscitation (T4) period.
At each time point, arterial blood samples were collected and EDTA-plasma isolated for metabolomics and laboratory analyses after immediate refrigerated centrifugation. In addition, hemoglobin concentration, blood lactate and electrolyte concentrations were measured (Cobas b-123, Roche, Switzerland) on top of the arterial and mixed-venous blood gasses already mentioned.
Targeted quantitative approach using a combined direct flow-injection and liquid chromatography (LC) tandem mass spectrometry (MS/MS) assay (AbsoluteIDQ 180 kit, Biocrates, Innsbruck, Austria) was applied for the metabolomics analysis of plasma swine samples. The assay quantifies 186 metabolites from six analyte groups: acylcarnitines, amino acids, biogenic amines, hexoses (sum of hexoses), glycerophospholipids and sphingomyelins. The method combines derivatization and extraction of analytes with the selective mass-spectrometric detection using multiple reaction monitoring (MRM) pairs. Methodological details and data preprocessing have been extensively reported in our previous articles (Ferrario et al., 2016; Cambiaghi et al., 2017). Measurable metabolites and their extended names and abbreviations are listed in Supplementary Table S1.
We used the Wilcoxon signed-rank test to assess significant changes from baseline to shock and resuscitation time points within the same group. We adopted the Wilcoxon rank-sum test to verify significant differences in the parameters values between the two groups (sham and shock animals) separately at each time point. Significance was considered with a p-value <0.05.
Multivariate and Multilevel Analyses of Metabolomics Data
We adopted the MultiLevel Simultaneous Component Analysis (MLSCA) and the multilevel partial least square discriminant analysis (MLPLSDA) model as proposed in De Noord and Theobald (2005) and Westerhuis et al. (2010). In the multilevel models the variation between subjects and the variation within subjects are separated.
The between subject variation in multilevel data analysis is performed on the average of the observations, whereas the within subject variation is performed on the net differences between paired observations. The initial step in multilevel model is therefore to separate the between subject variation from the within subject variation.
Total data matrix X (N × J) contains measurements for I subjects observed at Ki time points on J variables, where total number of observations is . An element xijki of matrix X, which contains a measurement of subject i on variable j at time point ki, can be decomposed as
where x.j. is the overall mean for variable j and xij. is the mean of subject i on variable j. The first term in Equation (1) is an offset that is constant across subjects and time points, the second term is a between-subject deviation, and the third term describes the within-subject deviation. Similarly, the sum of squares per variable can be separated in three parts, analogous to Analysis of Variance. The objective of these two-level model is to approximate the data, explaining the offset, between-subject, and within-subject variation as good as possible.
In the multilevel MLSCA model we assumed that the within –loadings matrices are the same for all the subjects. The MLSCA model for a subject i with data matrix Xi (Ki × J) containing the part of X with data from subject i (where i = 1,…, I) is the following:
where 1Ki is a Ki × 1 vector of ones, m (J × 1) contains the offsets of the J variables, tb,i (Rb × 1) contains the between-subject scores of subject i for Rb retained between-components, Pb (J × Rb) is the between-subject loading matrix, Tw,i (Ki × Rw) contains the within-subject scores for subject i for Rw retained within-components, Pw (J × Rw) is the within-subject loadings matrix and Ei (Ki × J) contains the residual.
In the multilevel MLSCA the offset, the between-subject loadings, and the within-subject loadings are supposed equal for all subject, and the within-subject loadings are supposed to be time invariant as well. This means that the interpretation of the components is equal for all subject. The between-subject and within-subject variations in the data are described in separate terms, and are thus not confounded as in ordinary PCA. Given these hypotheses, the three parts of the model can be solved separately.
The between-subject part can be obtained from the matrices Xb (I × J), which contain the mean vectors of the I subjects, by means of the singular value decomposition (SVD).
The within-subject part of the model can be estimated from the matrix Xw (N × J), which consists of the concatenated centered matrices Xi, by means of the singular value decomposition.
In the multilevel partial least square discriminant analysis (PLS-DA), instead of using the simple SVD as in MLSCA, the projections are guided by the outcomes (Y matrix).
In these study, we performed a MLSCA model on I = 6 septic shock animals at the three time points (Ki = 3) and we performed also two separate MLPLSDA: a model to analyze the within-subject variations occurring from baseline to shock and a second model to analyze the within-subject variations occurring from baseline to resuscitation. In such a way after the decomposition of the within- and between-subject variability a two class PLSDA is applied as the outcome consists in two distinct time point.
The models were analyzed on J = 17 parameters which are the metabolites under investigation (Glucogenic AA, Tot PC, LPC16+LPC18, Tot SM, Fisher ratio, PC ae C40:2, Alanine, Arginine, Citrulline, Glutamine, Glutamate, Ornithine, Tyrosine, ADMA, Creatinine, Glucose and Albumin).
In such a way, it is possible to understand in a multivariate fashion which species are most involved in determining the within-subject variation from baseline to shock condition and after resuscitation.
Finally, we applied a linear discriminant analysis (LDA) to two couples of specific molecular species (LPC and albumin, the ratio of alanine and glutamate and the ratio of glucose and lactate). Similarly to the MLSCA model, we considered I = 6 septic shock animals at the three time points (Ki = 3) and we performed a three class LDA to the within-subject variability estimated from the matrix Xw (N × J) where J = 2, i.e., the couple of species under investigation.
Changes in Hemodynamic, Blood Chemistry Parameters and Inflammatory Markers in Septic Shock Swine
Table 1 shows the changes induced by septic shock at each time point as regard the key hemodynamic parameters, arterial and mixed-venous blood gasses, hemoglobin concentration, blood lactate, and electrolyte concentrations. The data are in line with the characteristic phenotype of septic shock. During the initial phases of sepsis, vasodilation and capillary leak caused filling pressures to be low and successively the fluid resuscitation produced the typical hyperdynamic state with increased cardiac output and decreased systemic vascular resistance. Serum bicarbonate decrease and blood lactate levels increased, although not reaching extremely high levels. As shock progressed, metabolic acidosis worsened, and blood pH decreased in spite of resuscitation.
Values of cytokines at each time point are shown in Figure 1. In particular, we measured the tumor necrosis factor alpha (TNF-alpha), which is involved in systemic inflammation particularly during the acute phase reaction, and the interleukin 10 (IL-10), which is an anti-inflammatory cytokine capable to inhibit the synthesis of pro-inflammatory cytokines, such as TNF-alpha. As expected, both cytokine concentrations significantly increased and remained high after the shock was induced.
Figure 1. Cytokine plasma concentration. Boxplot of cytokine plasma concentration (pg/ml). Each circle represents the value of an animal. In the corner of each figure, boxplots of delta values are shown, i.e., the differences between shock and baseline values (BL) and the differences between the values after resuscitation (Resusc.) and baseline values. All the variations with respect to baseline were significant (Wilcoxon sign-rank test p < 0.05).
Time Course Alteration of Systemic Metabolic Features in Septic Shock Swine
The selected metabolite concentrations (μM) at each time point are reported in Table 2, while Figure 2 shows the delta (Δ) values of metabolites concentrations, as differences between shock vs baseline and between resuscitation vs baseline. In particular, we observed in shock a significant marked decrease of the overall amount of phosphatidylcholines species (total PC), of the sum of palmitoyl (C16) and stearoyl (C18) lysophosphatidylcholines species (LPC) at different saturation (C16:0+16:1+18:0+18:1+18:2) expressed as LPC16+LPC18, the plasmenylcholine PCaeC40:2, the overall amount of sphingomyelins (Total SM). A decrease trend in ornithine and the citrulline/glutamine ratio was also observed. As illustrated also in Figure 2, these metabolic features remained reduced even after resuscitation. From baseline to shock, the global arginine bioavailability ratio, i.e., arginine/(ornithine+citrulline) and the Fisher ratio decreased. On the contrary, the citrulline/arginine ratio increased from baseline to shock. Glutamine significantly increased from baseline to shock and remained elevated after resuscitation.
Figure 2. Variations of metabolites concentrations. Boxplot of delta values, i.e., the differences between shock and BL and the differences between the values after the resuscitation (Resusc.) and baseline values. Each circle refers to an animal. The symbols highlight the significant differences with respect to baseline values: ∗Wilcoxon sign-rank test p < 0.05, #Wilcoxon sign-rank test p = 0.063. The green shaded line marks Δ = 0. LPC16+LPC18 represents the sum of LPC species C16:0+16:1+18:0+18:1+18:2.
As regard the glucose-alanine cycle, we observed that alanine was significantly higher at resuscitation with respect the baseline values, glutamate significantly decreased after the resuscitation with respect to baseline and shock values and hexoses were significantly lower in shock with respect to baseline.
Liver Functionality Markers in Septic Shock Swine
Figure 3 shows the values distribution of blood analyses related to liver functionality. The concentration of albumin and gamma-glutamyltransferase (GGT) significantly decreased at either shock and resuscitation time points, whereas the concentration of aspartate aminotransferase (ASAT) increased significantly at the last time point. Alanine aminotransferase (ALAT) showed decreased concentrations from baseline, although the downward trend did not reach statistical significance. We computed also the ASAT/ALAT ratio. We found a significant increase from baseline to shock and in particular after the resuscitation [BL: 2.58 (1.56, 2.81), Shock: 2.89 (2.17, 2.94), Resusc.: 6.75 (6.21, 7.61)].
Figure 3. Laboratory values for each time point in shock group. Boxplot of laboratory values: each circle refers to an animal. In the corner of each figure, boxplots of delta values are shown, i.e., the differences between shock and BL and the differences between the values after resuscitation (Resusc.) and baseline values. The symbols highlight the significant difference with respect to baseline values: ∗Wilcoxon sign-rank test p < 0.05, and with respect to shock values: §Wilcoxon sign-rank test p < 0.05, °Wilcoxon sign-rank test p = 0.063. GGT, gamma-glutamyltransferase; ALAT, alanine aminotransaminase; ASAT, aspartate aminotransferase.
Metabolic Features Association in Septic Shock Swine
All the quantified metabolites are reported in Supplementary Table S2. Figure 4 highlights some associations between molecular species in septic shock model. We reasoned that since hypoalbuminemia was observed in our septic shock model and that albumin is a fatty acid transporter, it would have been plausible to find an association between these molecular species. Indeed, the upper panels show that both albumin and the sum of LPC16+LPC18 species decreased after shock. Similarly, when we focused on energy metabolic pathways, the glucose/lactate ratio decreased after shock passing from values higher than 1 at baseline to lower values than 1. This trend was paralleled by increased alanine/glutamate ratio particularly after the resuscitation. To note that the results of linear discriminant analysis (LDA) on the within-subject variability (lower panels) allowed the separation of each experimental phase from the others.
Figure 4. Relationships between lipids and albumin and between alanine and glucose. Upper panels: representation of concentrations of LPC16+LPC18 as sum of species C16:0+16:1+18:0+18:1+18:2 and albumin; the ratio between alanine and glutamate and the ratio between glucose and lactate. Each mark refers to an animal at the three time points (baseline, BL; shock; resuscitation, Resusc.). Lower panels represent the results of linear discriminant analysis (LDA) on the within-subject variability on the couples of species illustrated in the upper panel. x1 and x2 refers to the variation of each variable with respect to the subject mean and the lines mark the boundary obtained from LDA.
Multilevel and Multivariate Analysis of Metabolomics Data in Septic Shock Swines
Figure 5 shows the within-subject variation obtained from the MLSCA analysis. The baseline condition appears clearly separated from the other two time points.
Figure 5. Multilevel simultaneous component analysis of metabolites concentrations in shock animals. The first three principal components of the within-subject part of the MLSCA which are related to intra-subject changes among the different time points.
Table 3 illustrates the VIP scores of the MLPLSDA models. We can observe that the sum of LPC16+LPC18 species, phosphatidylcholines (PC) and sphingomyelins (SM) are ranked in the first positions of VIP score in both the two models, so that they can be considered important drivers of the within-subject variations occurring from baseline to shock and to resuscitation conditions. The bioavailability of arginine seems to have a greater role in the first model related to within-subject variations occurring from baseline to shock, as already reported in the univariate analysis. The VIP score of glutamate, which is a substrate in the gluconeogenesis from alanine, is higher than 1 in the second model only, whereas creatinine has less importance in the second model with respect to the first one.
Metabolic Features Comparison Between Septic Shock and Sham Swine
Figure 6 shows the distribution of the metabolite concentrations in shock and sham groups. In both shock and resuscitation period the overall amount of phosphatidylcholines species (total PC), the sum of LPC16+LPC18 species, the plasmenylcholine PCaeC40:2, the ornithine/citrulline ratio, the global arginine bioavailability ratio (GABR: arginine/ornithine+citrulline), were significantly lower in shock group with respect to sham animals, whereas alanine and citrulline/arginine were significantly higher. Hexoses and citrulline/glutamine were significantly enhanced during shock period only, whereas the Fisher ratio was significantly reduced after resuscitation in shock group with respect to sham animals.
Figure 6. Comparisons between shock and sham animals. Boxplot of metabolites concentration at each time point: baseline (BL), shock and resuscitation (Resusc.). The symbol †highlights the significant differences between shock and sham group (Wilcoxon rank sum test p < 0.05).
Since the etiology and pathophysiology of septic shock are complex, it is difficult to develop a clinically accurate animal model of septic shock with relevance as translational model of sepsis-induced shock. Nevertheless, our swine model replicates the hemodynamic dysfunctions and the systemic inflammation that are typically found in clinical setting, although it does not reproduce the comorbidities that may complicate the clinical course of human septic shock.
Importantly, our model displays an altered metabolic landscape pointing toward an early hepatic impairment similarly to what we previously reported for septic shock patients (Cambiaghi et al., 2017). We confirmed here that the time course of metabolome changes in septic shock swines are mainly characterized by reduced plasma phospholipids species trajectories, altered alanine-glucose cycle and deranged inter-organ amino acid metabolism, all metabolic traits reported in septic shock patients (Bruins et al., 2003; Murch et al., 2006; Kao et al., 2009; Ferrario et al., 2016; Ahn et al., 2017; Cambiaghi et al., 2017).
To our knowledge, this is the first study in which an experimental model of septic shock recapitulates the main metabolic derangements reported in a clinical setting of septic shock.
Shock swines are characterized by reduced levels of PC, SM and of a cluster of LPC species composed by C16:0-, C16:1-, C18:0-, C18:1-, and C18:2-LPC. Further, this metabolites pool is lower in shock animals with respect to sham animals in the corresponding period of shock and after resuscitation.
LPCs have a very complex metabolism and play key roles in inflammatory responses (Lin et al., 2005; Schmitz and Ruebsaamen, 2010). Interestingly LPC species enhance the suppressive function of CD4+ CD25+ regulatory T cells (Tregs) involved in anti-inflammatory signaling (Hasegawa et al., 2011; Osborn and Olefsky, 2012). Defects in Treg function are considered to play a role in hepatic inflammation (Speletas et al., 2011; Peiseler et al., 2012) and low levels of LPC16 and 18 species were reported in inflammatory liver disease (Tanaka et al., 2012; Lehmann et al., 2013; Maricic et al., 2014). These observations suggest that decreased circulating LPC might hamper their anti-inflammatory effects and conceivably, may promote an excessive immune response.
Because the liver is known to be a major source of circulating PC and LPC (Sekas et al., 1985; Baisted et al., 1988), it is plausible their plasma levels might be influenced by the functionality of hepatic lipid-modifying enzymes, whose expression can be regulated by pro-inflammatory cytokines such as TNF-α (Tanaka et al., 2012). It is likely that the production of pro-inflammatory cytokines devoted to regulate the hepatic response in sepsis may further injure hepatocytes.
Importantly, we found the same imbalance of LPC/PCs cycling observed in non-responsive septic shock patients (Cambiaghi et al., 2017), therefore strongly supporting that hepatic homeostasis is very early compromised. Such lipid levels might act as metabolic messengers in signaling early liver suffering (e.g., energy status) with consequent implication in the host response (Liu et al., 2014).
To note that both albumin and the observed LPC cluster decreased after shock in our model. Hypoalbuminemia is frequently observed during acute disease state, as albumin is a negative acute-phase protein. Hypoalbuminemia may be caused by the combination of increased transcapillary escape ratio, accelerated catabolism and impaired hepatic synthesis (Nicholson et al., 2000; Gatta et al., 2012). Albumin acts as fatty acid transporter (van der Vusse, 2009) and its reduction may play a critical role in the lipid economy of the body (Novak et al., 2017).
Therefore, our findings might suggest that the concomitant drop in albumin and the reduction in LPC species may represent a vicious cycle with detrimental effect on the intensity of inflammatory responses in the liver.
The plausible impairment in the early functionality of the liver in septic shock swine is also supported by the derangement of many interconnected energy-related metabolic pathways. The enhanced circulating alanine indicates an alteration in the glucose-alanine cycle as we previously observed in non-responsive septic shock patients (Cambiaghi et al., 2017). The higher plasma alanine in shock animals may be a sign of lower hepatic capacity for conversion of alanine to glucose, whose enhanced elaboration by the liver (hepatic gluconeogenesis) is a prominent feature of the solid organ response to injury and provides fuel to the cellular elements of the inflammatory response. The glucose/lactate ratio decreased after shock, indicating an imbalance between glycolysis and gluconeogenesis. This trend paralleled the increased alanine/glutamate ratio, particularly after the resuscitation. The alanine/glutamate ratio represents the alanine bioavailability from the impaired hepatic glucose-alanine cycle as part of the gluconeogenesis pathway. In fact, ALAT catalyzes the transfer of an amino group from alanine to α-ketoglutarate, to produce glutamate and pyruvate, which enters the gluconeogenesis and is used for glucose synthesis. Interestingly, we found a downward trend of ALAT in shocked and resuscitated animals, although not statistically significant.
As regards the transaminases trends, we may refer to subclinical hepatic alterations and not an acute liver dysfunction. The observational window after the insult is too short to have an hypoxic hepatitis, which is related to hypoperfusion in the presence of hypovolemia as well as inadequate cardiac output. Moreover, the increased levels of both ALAT and ASAT are observed in animal model of sepsis for longer time windows such as 24 h after the insult (Idrovo et al., 2015). However, in McGill (2016) the author reviews the limitation of using ALAT and ASAT for monitoring liver function. It is generally thought that aminotransferase elevations are due to cell damage with plasma membrane disruption. Apoptotic cell death can be thought of as cell implosion; a controlled process with minimal release of proteins into the extracellular space. However, oncotic necrosis, which can be thought of as cell explosion, can occur in some cells secondarily after induction of apoptosis. The hallmarks of apoptosis (e.g., cell shrinkage, chromatin condensation, apoptotic body formation) are clearly observable in the liver in hepatocyte apoptosis models during the early phase of injury when there is no significant elevation in serum ALAT (McGill, 2016).
Disturbance in metabolism is one of the most critical feature in sepsis/septic shock as response of the body to mobilize nutritional resources for organs and immune cells (Bruins et al., 2000, 2003; Whelan et al., 2014). In this complex and dynamic scenario, the time-course increase of plasma alanine and glutamine, major sources of nitrogen and carbon in inter-organ amino acid metabolism, may reflect a reduced rate in their utilization by organs in septic shock swine. The rise of circulating glutamine may support dysregulation in glutamine transport and extraction system, either suggesting reduced utilization by the gut and/or impaired hepatic transformation into urea cycle. Indeed, intermediates of the urea cycle were also disturbed in shock animals, where the global arginine bioavailability ratio (GABR: arginine/ornithine+citrulline) was significantly reduced compared to sham animals. GABR reduction may reduce NO, vasodilator essential to vascular homeostasis, and may lead to endothelial dysfunction (Kao et al., 2009).
The liver receives compounds absorbed or released by the gut through the portal vein before they gain access to the systemic circulation. The raised plasma level of glutamate, citrulline, arginine, seen in our septic shock swine model, are likely due to dysregulation of the gut-liver axis conversion of glutamine to citrulline and arginine. The gut accounts for the major part of this conversion, although splanchnic citrulline release to the systemic circulation may be limited by hepatic citrulline extraction (Van De Poll et al., 2007; Marini, 2016). To note that both arginine and citrulline plays an essential role in the immune response to sepsis and a disturbed arginine-nitric oxide (NO) metabolism is associated with deteriorating organ functions (Wijnands et al., 2012, 2015). Taken together our data are consistent with recent studies revealing liver dysfunction as an early event in sepsis (Nesseler et al., 2012; Recknagel et al., 2012; Wang et al., 2014).
The animal experimental data presented herein provide evidence that reduced plasma LPC/PC species trajectories, altered alanine-glucose cycle and deranged inter-organ amino acid exchange metabolism are early events in septic shock swine. These events occur within hours from infection and point to alterations in energy substrate utilization and hepatic biochemical perturbations. They may act as early metabolic features to assist in evaluating subclinical hepatic alterations. Being consistent with the metabolomics patterns differentiating non-responsive from responsive septic shock patients (Cambiaghi et al., 2017), such metabolic features may even represent therapeutic targets and therefore pave the way to improve the management of septic shock.
Metabolomics data are in the Supplementary Material; other experimental data are available upon request.
MF conceived and designed the study, did the statistical analysis, contributed to data interpretation, drafted, revised, and approved the manuscript. LB did the metabolomics analysis, contributed to data interpretation, revised, and approved the manuscript. FS and AH designed and developed the animal model, clinical animal data acquisition and interpretation, revised, and approved the manuscript. RP conceived the study, coordinated metabolomics analysis and data interpretation, drafted, and approved the manuscript.
This work was supported by EU FP7 Health Program, ShockOmics project, Grant #602706.
Conflict of Interest Statement
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.
We acknowledge Prof. Giuseppe Baselli, Politecnico di Milano, for his support and his leading role as Principal Investigator of the EU FP7 Health Program, ShockOmics project. We are grateful to all investigators participating in the ShockOmics project that contributed to discuss and design animal protocols: Dr. Giuseppe Ristagno, Dr. Bernardo Bollen Pinto, Dr. Karim Bendjelid, Dr. Fabio S. Taccone, Dr. Emiel H. Post, and Dr. Daniel De Backer.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.00011/full#supplementary-material
ALAT, alanine aminotransferase; ASAT, aspartate aminotransferase; FiO2, fraction of inspired oxygen; GGT, gamma-glutamyltransferase; LC, liquid chromatography; LDA, linear discriminant analysis; LPC, lysophospatidylcholine; LPC16+LPC18, the sum of palmitoyl (C16) and stearoyl (C18) lysophosphatidylcholines species at different saturation (C16:0+16:1+18:0+18:1+18:2); MAP, mean arterial pressure; MLPLSDA, MultiLevel partial least square discriminant analysis; MLSCA, MultiLevel simultaneous component analysis; MODS, multi-organ dysfunction syndrome; MRM, multiple reaction monitoring; MS/MS, tandem mass spectrometry; PC, phosphatidylcholine; PPV, arterial pulse pressure variation; SM, sphingomyelin.
Ahn, W. G., Jung, J. S., Kwon, H. Y., and Song, D. K. (2017). Alteration of lysophosphatidylcholine-related metabolic parameters in the plasma of mice with experimental sepsis. Inflammation 40, 537–545. doi: 10.1007/s10753-016-0500-6
Angus, D. C., Linde-Zwirble, W. T., Lidicker, J., Clermont, G., Carcillo, J., and Pinsky, M. R. (2001). Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Crit. Care Med. 29, 1303–1310. doi: 10.1097/00003246-200107000-00002
Bruins, M. J., Deutz, N. E. P., and And Soeters, P. B. (2003). Aspects of organ protein, amino acid and glucose metabolism in a porcine model of hypermetabolic sepsis. Clin. Sci. 104, 127–141. doi: 10.1042/CS20020275
Bruins, M. J., Soeters, P. B., and Deutz, N. E. P. (2000). Endotoxemia affects organ protein metabolism differently during prolonged feeding in pigs. J. Nutr. 130, 3003–3013. doi: 10.1093/jn/130.12.3003
Cambiaghi, A., Pinto, B. B., Brunelli, L., Falcetta, F., Aletti, F., Bendjelid, K., et al. (2017). Characterization of a metabolomic profile associated with responsiveness to therapy in the acute phase of septic shock. Sci. Rep. 7:9748. doi: 10.1038/s41598-017-09619-x
Chalkias, A., Spyropoulos, V., Koutsovasilis, A., Papalois, A., Kouskouni, E., and Xanthos, T. (2015). Cardiopulmonary arrest and resuscitation in severe sepsis and septic shock: a research model. Shock 43, 285–291. doi: 10.1097/SHK.0000000000000285
Corrêa, T. D., Jeger, V., Pereira, A. J., Takala, J., Djafarzadeh, S., and Jakob, S. M. (2014). Angiotensin II in septic shock: effects on tissue perfusion, organ function, and mitochondrial respiration in a porcine model of fecal peritonitis. Crit. Care Med. 42, e550–e559. doi: 10.1097/CCM.0000000000000397
Dellinger, R. P., Levy, M. M., Rhodes, A., Annane, D., Gerlach, H., Opal, S. M., et al. (2013). Surviving sepsis campaign: international guidelines for management of severe sepsis and septic shock: 2012. Crit. Care Med. 41, 580–637. doi: 10.1097/CCM.0b013e31827e83af
Duburcq, T., Durand, A., Dessein, A.-F., Vamecq, J., Vienne, J.-C., Dobbelaere, D., et al. (2017). Comparison of fluid balance and hemodynamic and metabolic effects of sodium lactate versus sodium bicarbonate versus 0.9% NaCl in porcine endotoxic shock: a randomized, open-label, controlled study. Crit. Care 21:113. doi: 10.1186/s13054-017-1694-1
Ferrario, M., Cambiaghi, A., Brunelli, L., Giordano, S., Caironi, P., Guatteri, L., et al. (2016). Mortality prediction in patients with severe septic shock: a pilot study using a target metabolomics approach. Sci. Rep. 6:20391. doi: 10.1038/srep20391
Hasegawa, H., Lei, J., Matsumoto, T., Onishi, S., Suemori, K., and Yasukawa, M. (2011). Lysophosphatidylcholine enhances the suppressive function of human naturally occurring regulatory T cells through TGF-β production. Biochem. Biophys. Res. Commun. 415, 526–531. doi: 10.1016/j.bbrc.2011.10.119
Idrovo, J. P., Yang, W. L., Jacob, A., Corbo, L., Nicastro, J., Coppa, G. F., et al. (2015). Inhibition of lipogenesis reduces inflammation and organ injury in sepsis. J. Surg. Res. 200, 242–249. doi: 10.1016/j.jss.2015.06.059
Kubiak, B. D., Albert, S. P., Gatto, L. A., Vieau, C. J., Roy, S. K., Snyder, K. P., et al. (2011). A clinically applicable porcine model of septic and ischemia/reperfusion-induced shock and multiple organ injury. J. Surg. Res. 166, e59–e69. doi: 10.1016/j.jss.2010.10.014
Lehmann, R., Franken, H., Dammeier, S., Rosenbaum, L., Kantartzis, K., Peter, A., et al. (2013). Circulating lysophosphatidylcholines are markers of a metabolically benign nonalcoholic fatty liver. Diabetes Care 36, 2331–2338. doi: 10.2337/dc12-1760
Li, C., Zhang, P., Cheng, X., and Chen, J. (2013). High-volume hemofiltration reduces the expression of myocardial tumor necrosis factor-alpha in septic shock pigs. Artif. Organs 37, 196–202. doi: 10.1111/j.1525-1594.2012.01536.x
Lin, P., Welch, E. J., Gao, X.-P., Malik, A. B., and Richard, D. Y. (2005). Lysophosphatidylcholine modulates neutrophil oxidant production through elevation of cyclic AMP. J. Immunol. 174, 2981–2989. doi: 10.4049/jimmunol.174.5.2981
Maricic, I., Girardi, E., Zajonc, D. M., and Kumar, V. (2014). Recognition of lysophosphatidylcholine by type II NKT cells and protection from an inflammatory liver disease. J. Immunol. 193, 4580–4589. doi: 10.4049/jimmunol.1400699
Murch, O., Collin, M., Sepodes, B., Foster, S. J., Mota-Filipe, H., and Thiemermann, C. (2006). Lysophosphatidylcholine reduces the organ injury and dysfunction in rodent models of Gram-negative and Gram-positive shock. Br. J. Pharmacol. 148, 769–777. doi: 10.1038/sj.bjp.0706788
Novak, F., Borovska, J., Vecka, M., Rychlikova, J., Vavrova, L., Petraskova, H., et al. (2017). Plasma phospholipid fatty acid profile is altered in both septic and non-septic critically ill: a correlation with inflammatory markers and albumin. Lipids 52, 245–254. doi: 10.1007/s11745-016-4226-x
Opal, S. M., and Patrozou, E. (2009). Translational research in the development of novel sepsis therapeutics: logical deductive reasoning or mission impossible? Crit. Care Med. 37, S10–S15. doi: 10.1097/CCM.0b013e3181921497
Peiseler, M., Sebode, M., Franke, B., Wortmann, F., Schwinge, D., Quaas, A., et al. (2012). FOXP3+ regulatory T cells in autoimmune hepatitis are fully functional and not reduced in frequency. J. Hepatol. 57, 125–132. doi: 10.1016/j.jhep.2012.02.029
Piper, R. D., Cook, D. J., Bone, R. C., and Sibbald, W. J. (1996). Introducing critical appraisal to studies of animal models investigating novel therapies in sepsis. Crit. Care Med. 24, 2059–2070. doi: 10.1097/00003246-199612000-00021
Recknagel, P., Gonnert, F. A., Westermann, M., Lambeck, S., Lupp, A., Rudiger, A., et al. (2012). Liver dysfunction and phosphatidylinositol-3-kinase signalling in early sepsis: experimental studies in rodent models of peritonitis. PLoS Med. 9:e1001338. doi: 10.1371/journal.pmed.1001338
Revie, J. A., Stevenson, D., Chase, J. G., Pretty, C. J., Lambermont, B. C., Ghuysen, A., et al. (2013). Evaluation of a model-based hemodynamic monitoring method in a porcine study of septic shock. Comput. Math. Methods Med. 2013:505417. doi: 10.1155/2013/505417
Rhee, C., Dantes, R., Epstein, L., Murphy, D. J., Seymour, C. W., Iwashyna, T. J., et al. (2017). Incidence and trends of sepsis in US hospitals using clinical vs claims data, 2009-2014. JAMA 318, 1241–1249. doi: 10.1001/jama.2017.13836
Singer, M., Deutschman, C. S., Seymour, C., Shankar-Hari, M., Annane, D., Bauer, M., et al. (2016). The third international consensus definitions for sepsis and septic shock (sepsis-3). JAMA 315, 801–810. doi: 10.1001/jama.2016.0287
Speletas, M., Argentou, N., Germanidis, G., Vasiliadis, T., Mantzoukis, K., Patsiaoura, K., et al. (2011). Foxp3 expression in liver correlates with the degree but not the cause of inflammation. Mediators Inflamm. 2011:827565. doi: 10.1155/2011/827565
Tanaka, N., Matsubara, T., Krausz, K. W., Patterson, A. D., and Gonzalez, F. J. (2012). Disruption of phospholipid and bile acid homeostasis in mice with nonalcoholic steatohepatitis. Hepatology 56, 118–129. doi: 10.1002/hep.25630
Van De Poll, M. C. G., Siroen, M. P. C., Van Leeuwen, P. A. M., Soeters, P. B., Melis, G. C., Boelens, P. G., et al. (2007). Interorgan amino acid exchange in humans: consequences for arginine and citrulline metabolism. Am. J. Clin. Nutr. 85, 167–172. doi: 10.1093/ajcn/85.1.167
Vassal, O., Bonnet, J., Barthelemy, A., Allaouchiche, B., Goy-Thollot, I., Louzier, V., et al. (2015). Renal haemodynamic response to amino acids infusion in an experimental porcine model of septic shock. Acta Anaesthesiol. Scand. 59, 598–608. doi: 10.1111/aas.12507
Westerhuis, J. A., van Velzen, E. J. J., Hoefsloot, H. C. J., and Smilde, A. K. (2010). Multivariate paired data analysis: multilevel PLSDA versus OPLSDA. Metabolomics 6, 119–128. doi: 10.1007/s11306-009-0185-z
Whelan, S. P., Carchman, E. H., Kautza, B., Nassour, I., Mollen, K., Escobar, D., et al. (2014). Polymicrobial sepsis is associated with decreased hepatic oxidative phosphorylation and an altered metabolic profile. J. Surg. Res. 186, 297–303. doi: 10.1016/j.jss.2013.08.007
Wijnands, K. A. P., Castermans, T. M. R., Hommen, M. P. J., Meesters, D. M., and Poeze, M. (2015). Arginine and citrulline and the immune response in sepsis. Nutrients 7, 1426–1463. doi: 10.3390/nu7031426
Wijnands, K. A. P., Vink, H., Briedé, J. J., van Faassen, E. E., Lamers, W. H., Buurman, W. A., et al. (2012). Citrulline a more suitable substrate than arginine to restore no production and the microcirculation during endotoxemia. PLoS One 7:e37439. doi: 10.1371/journal.pone.0037439
Keywords: swine, septic shock, metabolomics, liver functionality, energy metabolism, lactate
Citation: Ferrario M, Brunelli L, Su F, Herpain A and Pastorelli R (2019) The Systemic Alterations of Lipids, Alanine-Glucose Cycle and Inter-Organ Amino Acid Metabolism in Swine Model Confirms the Role of Liver in Early Phase of Septic Shock. Front. Physiol. 10:11. doi: 10.3389/fphys.2019.00011
Received: 20 August 2018; Accepted: 08 January 2019;
Published: 28 January 2019.
Edited by:Gaetano Santulli, Columbia University, United States
Reviewed by:Hai Li, Shanghai Jiao Tong University, China
Raymond J. Langley, University of South Alabama, United States
Copyright © 2019 Ferrario, Brunelli, Su, Herpain and Pastorelli. 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: Manuela Ferrario, email@example.com