Using a multiomics approach to unravel a septic shock specific signature in skeletal muscle

Baptiste Duceau, Michael Blatzer, Jean Bardon, Thibault Chaze, Quentin Giai Gianetto, Florence Castelli, François Fenaille, Lucie Duarte, Thomas Lescot, Christophe Tresallet, Bruno Riou, Mariette Matondo, Olivier Langeron, Pierre Rocheteau, Fabrice Chrétien, Adrien Bouglé, Baptiste Duceau, Michael Blatzer, Jean Bardon, Thibault Chaze, Quentin Giai Gianetto, Florence Castelli, François Fenaille, Lucie Duarte, Thomas Lescot, Christophe Tresallet, Bruno Riou, Mariette Matondo, Olivier Langeron, Pierre Rocheteau, Fabrice Chrétien, Adrien Bouglé

Abstract

Sepsis is defined as a dysregulated host response to infection leading to organs failure. Among them, sepsis induces skeletal muscle (SM) alterations that contribute to acquired-weakness in critically ill patients. Proteomics and metabolomics could unravel biological mechanisms in sepsis-related organ dysfunction. Our objective was to characterize a distinctive signature of septic shock in human SM by using an integrative multi-omics approach. Muscle biopsies were obtained as part of a multicenter non-interventional prospective study. Study population included patients in septic shock (S group, with intra-abdominal source of sepsis) and two critically ill control populations: cardiogenic shock (C group) and brain dead (BD group). The proteins and metabolites were extracted and analyzed by High-Performance Liquid Chromatography-coupled to tandem Mass Spectrometry, respectively. Fifty patients were included, 19 for the S group (53% male, 64 ± 17 years, SAPS II 45 ± 14), 12 for the C group (75% male, 63 ± 4 years, SAPS II 43 ± 15), 19 for the BD group (63% male, 58 ± 10 years, SAPS II 58 ± 9). Biopsies were performed in median 3 days [interquartile range 1-4]) after intensive care unit admission. Respectively 31 patients and 40 patients were included in the proteomics and metabolomics analyses of 2264 proteins and 259 annotated metabolites. Enrichment analysis revealed that mitochondrial pathways were significantly decreased in the S group at protein level: oxidative phosphorylation (adjusted p = 0.008); branched chained amino acids degradation (adjusted p = 0.005); citrate cycle (adjusted p = 0.005); ketone body metabolism (adjusted p = 0.003) or fatty acid degradation (adjusted p = 0.008). Metabolic reprogramming was also suggested (i) by the differential abundance of the peroxisome proliferator-activated receptors signaling pathway (adjusted p = 0.007), and (ii) by the accumulation of fatty acids like octanedioic acid dimethyl or hydroxydecanoic. Increased polyamines and depletion of mitochondrial thioredoxin or mitochondrial peroxiredoxin indicated a high level of oxidative stress in the S group. Coordinated alterations in the proteomic and metabolomic profiles reveal a septic shock signature in SM, highlighting a global impairment of mitochondria-related metabolic pathways, the depletion of antioxidant capacities, and a metabolic shift towards lipid accumulation.ClinicalTrial registration: NCT02789995. Date of first registration 03/06/2016.

Conflict of interest statement

The authors declare no competing interests.

© 2022. The Author(s).

Figures

Figure 1
Figure 1
Study design and general overview.
Figure 2
Figure 2
Principal component analyses of proteomic and metabolomic datasets. The two first dimensions are shown, defining the subspace maximizing the variance of the dataset. Every point represents an individual muscle proteome (A) or metabolome (B). The individuals are color-coded according to their group. Points that are close together tend to have similar proteome/metabolome. Large points and ellipses represent respectively the barycenter of each group and its 95% confidence interval. Dim indicates Dimension. Some samples were discarded due to quality checks. N = 31 patients analyzed for the proteomics (Septic shock n = 12, cardiogenic shock n = 9, brain dead n = 10) and N = 40 for the metabolomics (Septic shock n = 17, cardiogenic shock n = 6, brain dead n = 17).
Figure 3
Figure 3
Heat map of the differentially abundant metabolites. Each row represents a metabolite, each column a patient color-coded according to its group. The overlying dendrogram is a graphical representation of patient similarity assessed by the euclidean distance: patients in the same cluster are more similar than patients in two separate clusters. The relative intensities were scaled by rows. Some samples were discarded due to quality checks. N = 40 patients analyzed for the metabolomics (Septic shock n = 17, cardiogenic shock n = 6, brain dead n = 17).
Figure 4
Figure 4
Mitochondrial complexes were decreased in the S group. The subunits of each mitochondrial complexes are represented. The subunits that were not identified in the proteomic dataset are not shown (Complex I : MT-ND2, MT-ND3, MT-ND6; Complex III : MY-CYB, UQCRHL; Complex IV : COX4l2, COX6B2, COX7B, COX7B2, COX8C, COX10, COX15; Complex V : ATP5MF). The bar represents the mean difference (Log 2 Fold Change) between the S group and the two control groups, the error bar represents the 95% confidence interval. N = 31 patients analyzed for the proteomics (Septic shock n = 12, cardiogenic shock n = 9, brain dead n = 10). *p < 0.05; **p < 0.01; ***p < 0.001 (Tukey’s post-hoc test, Group S versus Reference).
Figure 5
Figure 5
Alterations of the antioxidant and ROS detoxifying proteins and metabolites. The bar represents the mean difference (Log 2 Fold Change) between the S group and the two control groups, the error bar represents the 95% confidence interval. N = 31 patients analyzed for the proteomics (Septic shock n = 12, cardiogenic shock n = 9, brain dead n = 10) and N = 40 for the metabolomics (Septic shock n = 17, cardiogenic shock n = 6, brain dead n = 17). *p < 0.05; **p < 0.01; ***p < 0.001 (Tukey’s post-hoc test, Group S versus Reference). Glrx, Glutaredoxin; GPx, Glutathione peroxidase; GR, Glutathione disulfide reductase; GSH-S, Glutathione synthetase; Prx, Peroxiredoxin; Sod, Superoxide dismutase; Trx, Thioredoxin; TrxR, Thioredoxin reductase. The mitochondrial location of an enzyme is indicated by (mt).

References

    1. Singer M, Deutschman CS, Seymour C, et al. The third international consensus definitions for sepsis and septic shock (sepsis-3) JAMA. 2016;315:801–810. doi: 10.1001/jama.2016.0287.
    1. Hotchkiss RS, Monneret G, Payen D. Sepsis-induced immunosuppression: from cellular dysfunctions to immunotherapy. Nat. Rev. Immunol. 2013;13:862–874. doi: 10.1038/nri3552.
    1. Angus DC, van der Poll T. Severe sepsis and septic shock. N. Engl. J. Med. 2013;369:840–851. doi: 10.1056/NEJMra1208623.
    1. Puthucheary ZA, Rawal J, McPhail M, et al. Acute skeletal muscle wasting in critical illness. JAMA. 2013;310:1591–1600. doi: 10.1001/jama.2013.278481.
    1. Batt J, Herridge M, Dos Santos C. From skeletal muscle weakness to functional outcomes following critical illness: a translational biology perspective. Thorax. 2019;74:1091–1098. doi: 10.1136/thoraxjnl-2016-208312.
    1. Brealey D, Brand M, Hargreaves I, et al. Association between mitochondrial dysfunction and severity and outcome of septic shock. Lancet. 2002;360:219–223. doi: 10.1016/S0140-6736(02)09459-X.
    1. Carré JE, Orban J-C, Re L, et al. Survival in critical illness is associated with early activation of mitochondrial biogenesis. Am. J. Respir. Crit. Care Med. 2010;182:745–751. doi: 10.1164/rccm.201003-0326OC.
    1. Vanhorebeek I, Latronico N, Van den Berghe G. ICU-acquired weakness. Intensive Care Med. 2020;46:637–653. doi: 10.1007/s00134-020-05944-4.
    1. Friedrich O, Reid MB, Van den Berghe G, et al. The sick and the weak: neuropathies/myopathies in the critically ill. Physiol. Rev. 2015;95:1025–1109. doi: 10.1152/physrev.00028.2014.
    1. Fernández-García M, Rojo D, Rey-Stolle F, et al. Metabolomic-based methods in diagnosis and monitoring infection progression. Metab. Interact. Infect. 2018;109:283–315. doi: 10.1007/978-3-319-74932-7_7.
    1. Lee J, Banerjee D. Metabolomics and the microbiome as biomarkers in sepsis. Crit. Care Clin. 2020;36:105–113. doi: 10.1016/j.ccc.2019.08.008.
    1. Neugebauer S, Giamarellos-Bourboulis EJ, Pelekanou A, et al. Metabolite profiles in sepsis: developing prognostic tools based on the type of infection. Crit. Care Med. 2016;44:1649–1662. doi: 10.1097/CCM.0000000000001740.
    1. Ferrario M, Cambiaghi A, Brunelli L, et al. Mortality prediction in patients with severe septic shock: A pilot study using a target metabolomics approach. Sci. Rep. 2016;6:20391. doi: 10.1038/srep20391.
    1. Parenica J, Jarkovsky J, Malaska J, et al. Infectious complications and immune/inflammatory response in cardiogenic shock patients: A prospective observational study. Shock Augusta Ga. 2017;47:165–174. doi: 10.1097/SHK.0000000000000756.
    1. Cuinet J, Garbagnati A, Rusca M, et al. Cardiogenic shock elicits acute inflammation, delayed eosinophilia, and depletion of immune cells in most severe cases. Sci. Rep. 2020;10:7639. doi: 10.1038/s41598-020-64702-0.
    1. Schwarz P, Custódio G, Rheinheimer J, et al. Brain death-induced inflammatory activity is similar to sepsis-induced cytokine release. Cell Transplant. 2018;27:1417–1424. doi: 10.1177/0963689718785629.
    1. Vandenbroucke JP, von Elm E, Altman DG, et al. Strengthening the reporting of observational studies in epidemiology (STROBE): explanation and elaboration. Int. J. Surg. Lond. Engl. 2014;12:1500–1524. doi: 10.1016/j.ijsu.2014.07.014.
    1. Giacomoni F, Le Corguillé G, Monsoor M, et al. Workflow4Metabolomics: A collaborative research infrastructure for computational metabolomics. Bioinformatics. 2015;31:1493–1495. doi: 10.1093/bioinformatics/btu813.
    1. Giacomoni F, Le Corguillé G, Monsoor M, et al. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47:D506–D515. doi: 10.1093/nar/gky1049.
    1. Boudah S, Olivier M-F, Aros-Calt S, et al. Annotation of the human serum metabolome by coupling three liquid chromatography methods to high-resolution mass spectrometry. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2014;966:34–47. doi: 10.1016/j.jchromb.2014.04.025.
    1. Lazar C, Gatto L, Ferro M, et al. Accounting for the multiple natures of missing values in label-free quantitative proteomics data sets to compare imputation strategies. J. Proteome Res. 2016;15:1116–1125. doi: 10.1021/acs.jproteome.5b00981.
    1. Schafer JL. Analysis of incomplete multivariate data. 1. Chapman and Hall/CRC; 1997.
    1. Bø TH, Dysvik B, Jonassen I. LSimpute: Accurate estimation of missing values in microarray data with least squares methods. Nucleic Acids Res. 2004;32:e34. doi: 10.1093/nar/gnh026.
    1. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30. doi: 10.1093/nar/28.1.27.
    1. Kanehisa M, Sato Y, Furumichi M, et al. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019;47:D590–D595. doi: 10.1093/nar/gky962.
    1. Mansmann U, Meister R. Testing differential gene expression in functional groups. Goeman’s global test versus an ANCOVA approach. Methods Inf. Med. 2005;44:449–453. doi: 10.1055/s-0038-1633992.
    1. Hummel M, Meister R, Mansmann U. GlobalANCOVA: exploration and assessment of gene group effects. Bioinformatics. 2008;24:78–85. doi: 10.1093/bioinformatics/btm531.
    1. Luo W, Friedman MS, Shedden K, et al. GAGE: Generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10:161. doi: 10.1186/1471-2105-10-161.
    1. Giai Gianetto Q, Wieczorek S, Couté Y, Burger T. A peptide-level multiple imputation strategy accounting for the different natures of missing values in proteomics data. Bioinformatics. 2020;16:1190.
    1. Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29:1830–1831. doi: 10.1093/bioinformatics/btt285.
    1. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32:2847–2849. doi: 10.1093/bioinformatics/btw313.
    1. Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 1995;57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x.
    1. Jiroutková K, Krajčová A, Ziak J, et al. Mitochondrial function in skeletal muscle of patients with protracted critical illness and ICU-acquired weakness. Crit. Care. 2015;19:448. doi: 10.1186/s13054-015-1160-x.
    1. Puthucheary ZA, Astin R, Mcphail MJW, et al. Metabolic phenotype of skeletal muscle in early critical illness. Thorax. 2018;73:926–935. doi: 10.1136/thoraxjnl-2017-211073.
    1. Murray-Stewart T, Dunston TT, Woster PM, Casero RA. Polyamine catabolism and oxidative damage. J. Biol. Chem. 2018;293:8736–18745. doi: 10.1074/jbc.TM118.003337.
    1. Sonee M, Martens JR, Mukherjee SK. Nicotinamide protects HCN2 cells from the free radical generating toxin, tertiary butylhydroperoxide (t-BuOOH) Neurotox. Res. 2002;4:595–599. doi: 10.1080/1029842021000045480.
    1. Picard M, Jung B, Liang F, et al. Mitochondrial dysfunction and lipid accumulation in the human diaphragm during mechanical ventilation. Am. J. Respir. Crit Care Med. 2012;186:1140–1149. doi: 10.1164/rccm.201206-0982OC.
    1. Langley RJ, Tsalik EL, van Velkinburgh JC, et al. An integrated clinico-metabolomic model improves prediction of death in sepsis. Sci .Transl. Med. 2013;5:19595. doi: 10.1126/scitranslmed.3005893.
    1. Khaliq W, Großmann P, Neugebauer S, et al. Lipid metabolic signatures deviate in sepsis survivors compared to non-survivors. Comput. Struct. Biotechnol. J. 2020;18:3678–3691. doi: 10.1016/j.csbj.2020.11.009.
    1. Gómez H, Kellum JA, Ronco C. Metabolic reprogramming and tolerance during sepsis-induced AKI. Nat. Rev. Nephrol. 2017;13:143–151. doi: 10.1038/nrneph.2016.186.
    1. Derde S, Hermans G, Derese I, et al. Muscle atrophy and preferential loss of myosin in prolonged critically ill patients. Crit. Care Med. 2012;40:79–89. doi: 10.1097/CCM.0b013e31822d7c18.
    1. Crossland H, Constantin-Teodosiu D, GardinerGreenhaff SMPL. Peroxisome proliferator-activated receptor γ agonism attenuates endotoxaemia-induced muscle protein loss and lactate accumulation in rats. Clin. Sci. Lond. Engl. 2017;131:1437–1447. doi: 10.1042/CS20170958.
    1. Menconi MJ, Arany ZP, Alamdari N, et al. Sepsis and glucocorticoids downregulate the expression of the nuclear cofactor PGC-1beta in skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 2010;299:E533–543. doi: 10.1152/ajpendo.00596.2009.
    1. Kaplan JM, Zingarelli B, Krallman K, et al. Phase 1 safety and pharmacokinetic study on the use of pioglitazone in critically ill patients with sepsis: a randomized clinical trial. Intensive Care Med. 2018;44:2006–2008. doi: 10.1007/s00134-018-5374-7.
    1. Lexell J, Henriksson-Larsén K, Sjöström M. Distribution of different fibre types in human skeletal muscles. 2. A study of cross-sections of whole m. vastus lateralis. Acta Physiol. Scand. 1983;117:115–122. doi: 10.1111/j.1748-1716.1983.tb07185.x.
    1. Johnson MA, Polgar J, Weightman D, Appleton D. Data on the distribution of fibre types in thirty-six human muscles. An autopsy study. J. Neurol. Sci. 1973;18:111–129. doi: 10.1016/0022-510x(73)90023-3.

Source: PubMed

3
Sottoscrivi