Elevated MMP-8 levels, inversely associated with BMI, predict mortality in mechanically ventilated patients: an observational multicenter study
Critical Care volume 27, Article number: 290 (2023)
The present study aimed to investigate the correlation between weight status and mortality in mechanically ventilated patients and explore the potential mediators.
Three medical centers encompassing 3301 critically ill patients receiving mechanical ventilation were assembled for retrospective analysis to compare mortality across various weight categories of patients using machine learning algorithms. Bioinformatics analysis identified genes exhibiting differential expression among distinct weight categories. A prospective study was then conducted on a distinct cohort of 50 healthy individuals and 193 other mechanically ventilated patients. The expression levels of the genes identified through bioinformatics analysis were quantified through enzyme-linked immunosorbent assay (ELISA).
The retrospective analysis revealed that overweight individuals had a lower mortality rate than underweight individuals, and body mass index (BMI) was an independent protective factor. Bioinformatics analysis identified matrix metalloproteinase 8 (MMP-8) as a differentially expressed gene between overweight and underweight populations. The results of further prospective studies showed that overweight patients had significantly lower MMP-8 levels than underweight patients ((3.717 (2.628, 4.191) vs. 2.763 (1.923, 3.753), ng/ml, P = 0.002). High MMP-8 levels were associated with increased mortality risk (OR = 4.249, P = 0.005), indicating that elevated level of MMP-8 predicts the mortality risk of underweight patients receiving mechanical ventilation.
This study provides evidence for a protective effect of obesity in mechanically ventilated patients and highlights the potential role of MMP-8 level as a biomarker for predicting mortality risk in this population.
Obesity augments the susceptibility to an array of illnesses, such as hypertension, hyperlipidemia, type 2 diabetes, coronary heart disease, stroke, and certain cancers, and is linked to amplified mortality and morbidity across diverse populations . However, recent studies have suggested that mildly obese has a protective effect in specific clinical settings . In some populations with chronic diseases, overweight or mildly obese patients exhibit a lowered mortality rate, i.e., the phenomenon is termed the obesity paradox [3,4,5]. Whether an obesity paradox exists in critically ill patients undergoing mechanical ventilation is yet controversial . Furthermore, racial and ethnic disparities are noted in the obesity paradox . Hence, additional studies are required to comprehend the association between weight status and mortality in diverse mechanically ventilated populations of different ethnicities.
Mechanical ventilation is a widely used respiratory support technique in intensive care medicine to assist or replace breathing in patients with respiratory failure. The repetitive opening and closing of alveoli during mechanical ventilation can cause shear stress and consequent mechanical damage, leading to the collapse of alveolar septa or deflection of fluid-filled alveoli, thereby resulting in uneven inflation within adjacent alveoli . This uneven intraalveolar inflation can lead to uneven oxygen concentrations within the alveoli and locally high oxygen concentrations, which can lead to oxidative stress [7, 9,10,11]. Oxidative stress can cause damage to cellular components, impairing cell function and survival, and initiate pathological processes, such as lung damage and alveolar edema, which have an adverse impact on the effectiveness of mechanical ventilation [12, 13].
In the present study, we aimed to investigate the correlation between weight status and mortality in mechanically ventilated patients and explore the potential role of the expression of matrix metalloproteinase-8 (MMP-8), the reactive oxygen species (ROS)-related gene, in the obesity paradox among East Asian populations receiving mechanical ventilation. Herein, we conducted a retrospective analysis of clinical data from 3301 mechanically ventilated patients in three medical centres to examine the correlation between patient body mass index (BMI) and hospital mortality and validate the existence of the obesity paradox in the East Asian population using machine learning algorithms. Next, we obtained data from the gene expression omnibus (GEO) database for bioinformatic analysis and identified the oxidative stress-related gene matrix metalloproteinase-8 (MMP-8) with differential expression between different weight groups. Finally, we conducted a prospective study in a defined cohort of 243 subjects, using enzyme-linked immunosorbent assay (ELISA) to measure plasma MMP-8 levels to validate the differential expression of MMP-8 in patients with different body weights and investigate the correlation between MMP-8 and mortality. This study provided evidence for the protective effect of overweight or mildly obese in mechanically ventilated patients and highlighted the potential role of MMP-8 expression as a biomarker for predicting mortality risk in this population.
Study design and data collection
This study was designed using both retrospective and prospective approaches, and clinical data were collected from three medical centres located in Wuhan, China: Tongji Hospital Sino-French New City Branch in Caidian District, Tongji Hospital Qiaokou Branch in Hankou District, and Tongji Hospital Optics Valley Branch in Wuhan East Lake High-Tech Development Zone. We conducted a retrospective study by collecting data from a total of 3301 critically ill patients who underwent mechanical ventilation over the past 5 years and analysed the correlation between their BMI and in-hospital mortality. The validation cohort was obtained from a prospective study conducted between October 2014 and September 2018 . Within this study, patients were categorized based on their utilization of mechanical ventilation, and exclusion criteria were applied to obtain a final cohort comprising 50 healthy adults and 193 patients who had undergone mechanical ventilation. The study design is depicted in Fig. 1.
Inclusion, exclusion, and subgroup criteria
Patients who were mechanically ventilated for at least 10 h were included in this study. Our collection of mechanically ventilated patients included patients treated with invasive mechanical ventilation, non-invasive mechanical ventilation, and high-flow nasal cannula. The following groups were excluded: (1) individuals < 18-years-old; (2) patients with acute cardiovascular and cerebrovascular events that resulted in death within 24 h; (3) those with severe trauma that resulted in death within 24 h. Additionally, the prospective cohort did not include patients with solid tumours. In this study, we collected a predominantly Chinese population, and the grouping criteria for BMI referred to Chinese standards: underweight: BMI < 18.5 kg/m2; normal range: 18.5 ≤ BMI < 24 kg/m2; overweight: BMI ≥ 24 kg/m2. The age grouping criteria pertain to the definition of older people’s age in developing countries and are categorized into two groups: those < 60-years-old and those ≥ 60-years-old.
The data were obtained from electronic medical records and encompassed various parameters, such as age, gender, BMI, arterial partial pressure of oxygen (PaO2), partial pressure of carbon dioxide, (PCO2); liver function indicators, kidney function indicators, coagulation function, sequential organ failure assessment (SOFA) score, and age-adjusted Charlson comorbidity index (aCCI). In-hospital mortality was the primary clinical outcome. Plasma samples were collected from a prospective cohort within 4 h of receiving mechanical ventilation. The locally weighted regression scatter plot smoothing (LOWESS) was utilized for preliminary exploration of the relationship between BMI and in-hospital mortality rate.
Regression discontinuity design (RDD)
Herein, we employed restrictive cubic spline (RCS) analyses to evaluate the nonlinear correlation between BMI and mortality rate. We utilized RDD to evaluate the causal correlation between BMI and mortality rate. RDD analysis is similar to a randomized experiment . It uses the sudden change point of the independent variable to explore the changes in the dependent variable, thereby inferring the causal correlation between the independent and dependent variables . We utilized five machine learning models, including GaussianNB, LogisticRegression, DecisionTreeClassifier, RandomForestClassifier, and GradientBoostingClassifier, to obtain the optimal risk prediction model for predicting mortality risk in mechanically ventilated patients by comparing the area under the receiver operating characteristic (ROC) curve and model performance. The RDD with BMI was considered the grouping variable, and individual mortality risk was the outcome variable.
Bioinformatic analysis was performed to identify DEGs (differentially expressed genes) between overweight and underweight populations. To procure the dataset comprising individuals of varying weight statuses, we retrieved the GSE18897 dataset from the GEO database, consisting of 20 underweight and 20 overweight patients . Similarly, we retrieved the mechanical ventilation dataset to obtain the GSE21802 dataset, which comprised four healthy individuals and six mechanically ventilated groups . The source of the specimen was whole blood collected on the first day of admission. The gene expression profiles were extracted using the “GEOquery package” and the “limma package,” followed by the normalization of raw data [19, 20]. The intersection of the two gene sets was obtained to identify the DEGs during mechanical ventilation in patients with varying weight categories.
Identifying differentially expressed genes
DEGs were identified based on the following criteria: |Log2Fold-Change|> 1 and P < 0.05. The variance threshold method in machine learning was applied to screen the feature genes. This method is a feature selection technique that relies on the variance of a feature. It operates on the principle of eliminating features with a variance below a specified threshold, as these features might have a low impact on the prediction of the target variable. Typically, features with high variance exhibit strong predictive capabilities for the target variable. Diversified approaches can be employed to compute feature variance, and one popular method is the use of median absolute deviation (MAD), a robust variance estimation technique that remains unaffected by outliers. The MAD value was calculated for each gene, and the median MAD value × 1.4826 was employed as the threshold for variance screening. This process filtered out the genes with high variance and identified the feature genes. Finally, the top 20 genes based on their characteristic value were visualized using a heat map.
To obtain the co-expressed genes, we intersected the set of differential genes with the feature genes identified by machine learning algorithms. The protein–protein interactions (PPI) were established using the online database STRING (string-db.org) with medium confidence (0.4) . To confirm the hub genes using topological analysis methods of the degree algorithm, we utilized the CytoHubba plugin, which is a plugin of Cytoscape 3.8.0 (University of California, San Diego, CA, USA) . Next, we employed Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis to demonstrate the molecular function and critical pathways linked to the intersection genes using the package “org.Hs.eg.db” and “clusterProfiler” [24, 25].
Next, to identify hub genes associated with oxidative stress, we downloaded oxidative stress-related genes from the “regulation of reactive oxygen species metabolic process” gene set in the MSigDB database (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp). To investigate specific target genes related to the regulation of the ROS pathway, we used the intersection of hub genes and oxidative stress-related genes. Mechanically ventilated populations were frequently accompanied by various underlying diseases. We attempted to analyse the expression of MMP-8 in various diseases, as the origin of a critical illness often determines the risk of mortality . Therefore, we analysed DEGs expression levels in 15 non-neoplastic (renal failure, liver failure, coronary heart disease, chronic heart failure, infective endocarditis, stroke, idiopathic pulmonary fibrosis, asthma, pneumonia, sepsis, Crohn’s disease, ulcerative colitis, primary biliary cholangitis, primary sclerosing cholangitis, chronic obstructive pulmonary disease) and 20 neoplastic diseases to investigate whether differences in underlying diseases have an impact on DEGs expression levels (Additional file 1: Fig. S1). The non-neoplastic sample expression profile was retrieved from the Gene Expression Omnibus (GEO) database and analysed using the GEO2R online tool (http://www.ncbi.nlm.nih.gov/geo/geo2r). The tumor sample expression profile were sourced from The Cancer Genome Atlas (TCGA) database and analysed with Sangerbox .
Plasma sample collection
Within 4 h of admission, the whole blood samples of mechanically ventilated patients admitted to the intensive care unit (ICU) were collected for MMP-8 estimation. The blood samples were placed in tubes containing EDTA and a mixture of protease inhibitors (Sigma-Aldrich, MO, USA). The plasma was obtained by centrifugation of the blood at 3000 g, 4 °C for 15 min. The remaining blood samples were clarified by centrifugation and stored at − 80 °C.
Plasma MMP-8 levels were measured using a commercially available ELISA kit (HM11140, Bio-swamp, Wuhan, China) according to the manufacturer’s instructions (http://web.bio-swamp.com/productSpecification/29188).
We performed bioinformatics analysis using R software (version 4.1.0) and machine learning using anaconda3 software. Stata14.0 software was utilized for statistical analysis. To assess the associations between the variables under investigation, Spearman correlation and partial correlation analyses were performed. Continuous variables were expressed as mean ± standard deviation (SD) or median (lower quartile—upper quartile) based on their distribution. Categorical variables were presented as frequencies and percentages. For continuous variables, when the data met the normal distribution and satisfied the homogeneity of variance test, t-test was employed for two-group comparison, and one-way analysis of variance (ANOVA) was used for three-group comparison. When the data met the normal distribution but not the homogeneity of variance test, the Welch t-test was utilized for two-group comparison, and Welch one-way ANOVA was used for three-group comparison. When the data did not meet the normal distribution, Wilcoxon test was used for two-group comparison, and Kruskal–Wallis test was employed for three-group comparison. The chi-square test was utilized to compare the categorical variables between groups. Logistic regression analysis was utilized to evaluate the correlation between weight status and mortality, adjusted for age, sex, aCCI, and SOFA score. The significance level is set at two-tailed P < 0.05.
This study was approved by the ethics committee of the hospital where the study was conducted (Wuhan, China; Approval No. TJ-IRB20230313). This study was conducted in line with the Helsinki Declaration of 1975. As the study was observational study in nature, informed consent was exempted.
Demographic and clinical characteristics of the study population
As shown in Table 1, the retrospective cohort consisted of 3301 mechanically ventilated patients, with the majority being young individuals (51.4%) with a male predominance (67.5%). The median aCCI score was 2 (1–3), and the median SOFA score was 8 (7–10). The overall mortality rate was 17.15%. The collective median BMI was 21.36 (19.86–23.5) kg/m2, while the median BMI for survivors was 21.47 (19.98–23.72) kg/m2 and 20.73 (19.18–21.95) kg/m2 for non-survivors (P < 0.001). Additionally, no gender-based variations were observed between the group of surviving and deceased individuals. Moreover, no significant differences were detected in the inhaled oxygen concentration, white blood cell (WBC) count, and 24-h urine output (P > 0.05). Furthermore, we display the primary diagnosis for admission and the primary cause for mechanical ventilation in Additional file 2: Fig. S2.
BMI is an independent protective factor against in-hospital mortality
Logistic regression-based univariate and multivariate analyses, as well as, variance inflation factor (VIF) test determined BMI is an independent protective factor. In the univariate analysis, BMI showed a significant protective effect with a P-value < 0.001. In the multivariate analysis, after adjusting for other potential confounding factors, BMI remained a significant protective factor with an adjusted odds ratio of 0.898 [95% confidence interval (CI): 0.873–0.925) (Table 2). Additionally, the VIF values for all variables were < 10, indicating no presence of multicollinearity (Additional file 3: Table S1). Therefore, our findings suggested a protective role of BMI against in-hospital mortality.
Correlation between weight status and in-hospital mortality
The correlation between BMI and in-hospital mortality was initially investigated by utilizing the LOWESS curves. The analysis revealed a U-shaped association between BMI and in-hospital mortality, indicating that a slight increase in BMI led to a decrease in mortality (Fig. 2a). Furthermore, in order to gain a more in-depth understanding of the correlation between BMI and in-hospital mortality, we classified BMI into five groups based on quintiles. It was observed that the odds ratio (OR) of BMI decreased progressively with increasing quintile in any model, and the trend was statistically significant (all P for trend < 0.05, Fig. 2b–j). A non-linear correlation between BMI and in-hospital mortality was discovered in all subgroups of diseases after analysing various comorbidity subgroups (all P < 0.01, Fig. 2k–p).
As shown in Fig. 3a, a nonlinear correlation was further detected between BMI and in-hospital mortality using RCS analysis (all P < 0.001). This correlation remained significant even after adjusting for potential confounding factors (all P < 0.05). Additionally, the different BMI ranges were compared with in-hospital mortality using the normal range of BMI as a reference. The observed trend is similar to that of RCS analysis (Additional file 3: Table S2). Of the five machine learning algorithm models, the Gradient Boost machine learning algorithm showed the highest accuracy (Fig. 3b, c). The Gradient Boost machine learning algorithm was used to calculate the risk score for patients presenting in-hospital mortality. Interaction analysis showed the interaction of BMI with age (P for interaction < 0.001), Lac (P for interaction = 0.084), BUN (P for interaction = 0.004), AST (P for interaction = 0.014), PCO2 (P for interaction = 0.001) and SOFA (P for interaction < 0.001, Fig. 3d-i). Those covariates was controlled by conducting a RDD analysis.
RDD shows a significant reduction in mortality at BMI = 22.85 kg/m2 (P = 0.004, Fig. 3j), indicating a causal association between BMI and in-hospital mortality. In addition, sensitivity tests, continuity tests, placebo tests, and balance checks were conducted to demonstrate the robustness of the results (Fig. 3k–o, Additional file 3: Table S3). Therefore, these data indicate that the obesity paradox exist in East Asian populations receiving mechanical ventilation.
MMP-8 was differentially expressed in underweight and overweight populations
A total of 173 DEGs were identified in the GSE18897 dataset using bioinformatics analysis, of which 147 were upregulated and 26 were downregulated (Fig. 4a–f). In the GSE21802 dataset, 1054 DEGs were identified, including 487 upregulated and 567 downregulated (Fig. 4g–l). Enrichment analysis was conducted separately for the differential genes in each dataset, and GO molecular function enrichment analysis revealed enrichment of reactive oxygen metabolism-related processes (Fig. 4e, k). KEGG pathway enrichment analysis indicated the presence of oxidative stress-related signaling pathways, such as the Foxo (Forkhead box O) and HIF1 (Hypoxia-inducible factor 1) signaling pathways, in both datasets (Fig. 4f, l).
The variance threshold method identified 8926 and 15,490 signature genes in the GSE18897 and GSE21802 datasets, respectively (Fig. 5a, b). The intersection of the DEGs and signature gene datasets revealed 27 co-expressed genes (Fig. 5c). Further GO molecular function and KEGG pathway enrichment analysis of these 27 co-expressed genes suggested enrichment of oxidative stress-related signals (Fig. 5d, e). The PPI network and Cytoscape software identified 17 hub genes (Fig. 5f). To further identify the oxidative stress-related hub genes, an analysis was conducted between the hub genes and oxidative stress-related genes using a Venn diagram; the oxidative stress-related hub gene MMP-8 was obtained (Fig. 5g). Significant differences were noted in MMP-8 RNA levels in both datasets (GSE18897: 6.144 (4.034–11.662) vs. 2.657 (2.197–3.069), P < 0.001, Fig. 5h; GSE21802: 119.926 (116.764–126.732) vs. 508.967 (332.292–1570.505), P = 0.011, Fig. 5i), indicating that MMP-8 may have differential expression in patients with varying weights while receiving mechanical ventilation.
After examining the levels of MMP-8 expression in non-oncological illnesses, noteworthy distinctions were discovered in severe asthma, pneumonia, and sepsis (|Log2Fold-Change|> 1 and P < 0.05). Subsequently, upon further scrutinizing the MMP-8 expression levels in both tumor and normal tissues, noteworthy variations were discovered across various tumor types, with particularly heightened expression evident in gastrointestinal tumor and lung cancer (Additional file 1: Fig. S1, Additional file 3: Table S4).
Elevated MMP-8 levels predict the mortality risk of underweight patients receiving mechanical ventilation
To authenticate the findings of the bioinformatics analysis, we executed a prospective study comprising 50 healthy individuals and 193 mechanically ventilated patients. Demographics and clinical characteristics of the study population are shown in Table 3. In comparison to the healthy population, patients receiving mechanical ventilation exhibited a considerably higher concentration of plasma MMP-8 (2.702 (2.305–3.313) vs. 3.027(2.336–3.959) ng/ml, P = 0.030). Subsequently, we categorized all mechanically ventilated patients into three subgroups based upon their BMI levels: underweight (BMI < 18.5 kg/m2, n = 55), normal weight (BMI: 18.5–24 kg/m2, n = 90), and overweight (BMI > 24 kg/m2, n = 48). The demographic characteristics of the three groups, including age and gender, were not significantly different. Notably, plasma MMP-8 levels, as well as in-hospital mortality rates, varied significantly among the three groups of patients (all P < 0.05). Further analysis demonstrated that overweight patients had considerably lower MMP-8 levels than underweight patients ((3.717 (2.628- 4.191) vs. 2.763 (1.923- 3.753), ng/ml, P = 0.002). The mortality rates of patients in underweight, normal weight, and overweight groups were 13.6%, 15.6%, and 7.4%, respectively (Table 4). In contrast to the normal weight group, the obese group showed a prolonged mechanical ventilation duration (30 (24, 37) VS. 44.5 (26, 70), P < 0.001), albeit with a statistically significant reduction in ICU hospitalization time (123.5 (64–225) vs. 83 (56.5–149), P = 0.036). These findings together suggest that higher levels of MMP-8 indicate a higher risk of mortality for underweight patients undergoing mechanical ventilation.
Associations of MMP-8 with different parameters
The heat map demonstrated that MMP-8 expression levels correlated with numerous laboratory indicators (Additional file 4: Fig. 3). In a comparison of subgroups by weight status, MMP-8 levels were significantly higher in the underweight group than in the normal weight and overweight groups (all P < 0.05). Spearman's correlation analysis further showed a significant negative correlation between MMP-8 and BMI (r = − 0.143, P < 0.05). The correlations retained their significance even after patriating out the covariation estimates (all P < 0.05).
The dose–response correlation between patient MMP-8 expression levels and in-hospital mortality was examined using the LOWESS curve, which revealed an inverse S-shaped association between MMP-8 expression levels and ICU mortality. As illustrated in Fig. 6, patients were categorized into six groups based on MMP-8 concentration, and the mortality rates were computed for each group. The results of the LOWESS curve analysis were supported by the histograms. Furthermore, the S-shaped association was confirmed using a cubic term regression model (Additional file 3: Table S5). An S-shaped correlation was established between MMP-8 expression levels and in-hospital mortality, with two extreme points of (2.322, -0.929) and (3.867,0.903), respectively, having a trinomial coefficient (P < 0.001).
Due to an inverse “S” correlation between MMP-8 and in-hospital mortality, we conducted a segmented logistic regression analysis using the two vertices as cutoff points. The results revealed a significant association between high levels of MMP-8 and increased mortality risk. Specifically, MMP-8 levels (2.322–3.867 ng/mL) were associated with an odds ratio of 4.249 (P = 0.005) for mortality risk (Table 5). However, the results were not significant when the MMP was > 3.867 ng/mL or < 2.322 ng/mL (all P > 0.05), indicate that there is an increased risk of mortality in underweight patients undergoing mechanical ventilation with an increase in MMP-8 levels (OR (95% CI) = 4.249(1.546–11.675), P = 0.005). Therefore, MMP-8 level may serve as a potential predictive indicator of mortality risk in this population.
Finally, we explore the value of MMP-8 as a molecular marker to predict in-hospital mortality. The study put forward a composite index that combined MMP-8 levels and relevant clinical parameters, including BMI, Lac, PAO2, SAO2, FIO2, and BUN. The composite index showed promising potential for diagnosis, outperforming the reference SOFA score in predicting in-hospital mortality (0.869 (0.816–0.922) vs. 0.600 (0.519–0.681), P < 0.001, Additional file 4: Fig. S3, Additional file 3: Table S6).
The obesity paradox is observed in East Asian populations receiving mechanical ventilation, with obesity emerging as an independent protective factor for mechanically ventilated patients based on the findings of our retrospective cohort study. Subsequently, we conducted bioinformatic analysis to identify potential mediators and eventually identified MMP-8 as a pivotal element in the emergence of this paradox in mechanically ventilated patients. We subsequently assessed the level of MMP-8 expression in various disease gene sets and found a significant increase in its expression level among the sepsis, severe asthma, and pneumonia groups as compared to healthy individuals. This finding suggests that hypoxia and infection are the principal factors underlying the rise in MMP-8. Given that critically ill patients requiring mechanical ventilation often exhibit a high incidence of hypoxia and infection, the elevated levels of MMP-8 observed in this cohort can likely be attributed to these common medical conditions.
To corroborate our findings, we carried out a prospective cohort study, which once again demonstrated a lower mortality rate in the overweight group when compared to the underweight group, confirming previous research highlighting the existence of the obesity paradox in critically ill patients . Furthermore, we opted not to include tumor patients in our prospective cohort since MMP-8 expression levels vary considerably across different types of tumours. We also observed differential expression of MMP-8 in patients with varying weight statuses, indicating its potential as a biomarker for predicting clinical outcomes in this group. Currently, the underlying mechanisms of this paradox are inconclusive and can be attributed to various factors, including enhanced immune function, increased energy reserves, and reduced catabolic responses. Our study contributes new insights into the paradox from the perspective of oxidative stress.
Our study also found an inverted S-shaped association between MMP-8 and hospital mortality. Due to the lack of dose–response studies on MMP-8 levels in mechanically ventilated populations, we propose a potential mechanism for the inverse S-shaped association between MMP-8 and in-hospital mortality. High-serum MMP-8 levels have been reported to be associated with white blood cell counts in the literature. Furthermore, existing studies suggest that elevated levels of MMP-8 in the serum are linked to systemic inflammation and adverse outcomes , indicating that the expression of MMP-8 in patients is related to the body's inflammatory response to external stimuli. Prior to the first nadir, the inflammatory response is inadequate, and as MMP-8 gradually rises, the inflammatory response strengthens, leading to reduced mortality (P = 0.083). When MMP-8 exceeds the threshold, a threshold effect occurs, indicating an excessive inflammatory response in the organism and causing mortality to rise rapidly with increasing MMP-8 levels (P = 0.005). Once the maximum point is surpassed, a saturation effect takes place, and the association between in-hospital mortality and MMP-8 is no longer statistically significant (P = 0.138).
MMP-8, known as collagenase-2 or neutrophil collagenase, belongs to the MMP family and expressed extensively within various cell types, including endothelial cells, vascular smooth muscle cells, neutrophils, and macrophages [30, 31]. Neutrophils predominantly store MMP-8 in gelatinase granules, which are released at the site of inflammation during inflammatory responses [32, 33]. Under normal physiological conditions, MMP-8 is stored as an inactive proenzyme and is biologically non-functional. However, in the presence of inflammation and ROS, MMP-8 undergoes activation through a cysteine conversion mechanism, ultimately converting from an inactive to an active form. Additionally, hypoxia enhances the process of neutrophil degranulation, leading to the transfer of multiple proteases to the extracellular compartment . These phenomena might account for the rise in MMP-8 levels during mechanical ventilation.
The present study also discovered elevated levels of MMP-8 in lean patients, consistent with previous reports . In a mouse model of blunt chest trauma, lean mice had significantly higher levels of MMP-8 in the first 6 h compared to obesity . Additionally, a population-based study showed that plasma MMP-8 levels were lower in obese women than in lean women , consistent with our conclusion that MMP-8 expression differs across weight categories. Our study also demonstrated a correlation between MMP-8 and disease severity in mechanically ventilated patients, with higher MMP-8 levels associated with increased in-hospital mortality. It has been demonstrated in the colon cancer population that heightened serum MMP-8 levels are correlated with inflammatory responses and unfavorable outcomes, thus substantiating the notion . Therefore, MMP-8 may serve as a potential biomarker for predicting in-hospital mortality in this patient population. Two potential reasons contribute to these findings: Firstly, mechanical ventilation is known to cause some degree of lung damage, and the level of MMP-8 is closely linked to lung tissue damage . Secondly, increased levels of MMP-8 have been linked to poor prognosis, and clinical studies conducted on patients with severe COVID-19 lung disease revealed significantly higher levels of MMP-8 and MMP19 in the non-survival group compared to the survival group . The underlying mechanism may be that the overexpression of MMP-8 increases lipid peroxidation and exacerbates lung injury .
This study has the following strengths. To begin with, it provides evidence for the protective effect of obesity in mechanically ventilated patients. This finding has clinical implications for the management of critically ill patients that may guide decision-making regarding nutritional support and weight management in this population. Moreover, this study highlights the potential role of MMP-8 expression as a biomarker for predicting mortality risk in mechanically ventilated patients. The identification of a biomarker accurately predicts that mortality risk can have significant implications for patient care, allowing early interventions and potentially improving outcomes. Finally, the study takes a multidisciplinary approach, integrating retrospective and prospective designs, bioinformatic analysis, and ELISA to gain a nuanced understanding of the underlying mechanisms and potential clinical implications of these findings. The multidisciplinary approach allows for comprehensive comprehension of the mechanisms underlying these findings and their potential clinical implications. Nevertheless, the present study has several limitations. Firstly, the population of critically ill patients receiving mechanical ventilation is highly heterogeneous, encompassing a range of diseases that cause severe respiratory failure, such as severe pneumonia, interstitial pneumonia, and chronic obstructive pulmonary disease. Moreover, different comorbidities, such as age, genetic susceptibility, and sources of infection, can complicate the level of MMP-8 expression. Although the analysis incorporated SOFA and aCCI scores to assess the states of various diseases, it is still challenging to eliminate infections as confounding factors. Secondly, due to the clinical nature of this study and the absence of in vitro cellular experiments, differences in MMP-8 mRNA and protein expression levels under cellular experiments were not analyzed. Thirdly, there are limitations in the BMI indicator, which only takes into account height and weight and cannot distinguish between muscle and fat. Also, BMI does not reflect fat distribution. Finally, our study, conducted from a clinical perspective, has revealed a correlation between MMP-8 and BMI, suggesting that MMP-8 may be a potential underlying cause of the obesity paradox. However, to further explore the causal relationship between MMP-8 and low BMI, we require the support of fundamental scientific methods. Thus, this will also be a key direction of our future research.
This study provides new insights into the correlation between weight status and mortality rates in mechanically ventilated patients. Overweight and mild obesity is an independent protective factor, and MMP-8 is identified as a potential biomarker for predicting the outcomes in non-neoplastic condition. Further studies are required to validate the findings and explore the underlying mechanisms, which could have significant clinical implications and improve patient outcomes.
Availability of data and materials
The data and code used in this study could be obtained from the corresponding author upon reasonable request.
Activated partial thromboplastin time
Acute respiratory distress syndrome
Bladder urothelial carcinoma
Body mass index
Breast invasive carcinoma
Blood urea nitrogen
Colon adenocarcinoma/rectum adenocarcinoma esophageal carcinoma
Differentially expressed genes
Enzyme-linked immunosorbent assay
Fraction of inspired oxygen
Forkhead box O
Gene expression omnibus
Hypoxia-inducible factor 1
Head and Neck squamous cell carcinoma
Kyoto encyclopedia of genes and genomes
Pan-kidney cohort (KICH + KIRC + KIRP)
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Liver hepatocellular carcinoma
Lung squamous cell carcinoma
Median absolute deviation
Mean arterial pressure
Arterial partial pressure of oxygen
Regression discontinuity design
- ROC curve:
Receiver operating characteristic curve
Restrictive cubic spline
Sequential organ failure assessment
Stomach and esophageal carcinoma
The cancer genome atlas
Uterine corpus endometrial carcinoma
White blood cell count
Bray GA, Kim KK, Wilding JPH, World Obesity F. Obesity: a chronic relapsing progressive disease process. A position statement of the World Obesity Federation. Obes Rev. 2017;18(7):715–723.
Zhou D, Wang C, Lin Q, Li T: The obesity paradox for survivors of critically ill patients. Critical Care 2022, 26(1).
Horwich TB, Fonarow GC, Clark AL. Obesity and the obesity paradox in heart failure. Prog Cardiovasc Dis. 2018;61(2):151–6.
Persaud SR, Lieber AC, Donath E, Stingone JA, Dangayach NS, Zhang X, Mocco J, Kellner CP. Obesity paradox in intracerebral hemorrhage national inpatient sample analysis. Stroke. 2019;50(4):999–1002.
Han SJ, Boyko EJ. The evidence for an obesity paradox in type 2 diabetes mellitus. Diabetes Metab J. 2018;42(3):179–87.
De Jong A, Wrigge H, Hedenstierna G, Gattinoni L, Chiumello D, Frat J-P, Ball L, Schetz M, Pickkers P, Jaber S. How to ventilate obese patients in the ICU. Intensive Care Med. 2020;46(12):2423–35.
Kleine C-E, Moradi H, Streja E, Kalantar-Zadeh K. Racial and ethnic disparities in the obesity paradox. Am J Kidney Dis. 2018;72(5):S26–32.
Dresse C, Joris JL, Hans GA. Mechanical ventilation during anaesthesia: Pathophysiology and clinical implications. Trends in Anaesthesia and Critical Care. 2012;2(2):71–5.
Hafner C, Pramhas S, Schaubmayr W, Assinger A, Gleiss A, Tretter EV, Klein KU, Scharbert G. Brief high oxygen concentration induces oxidative stress in leukocytes and platelets: a randomized cross-over pilot study in healthy male volunteers. Shock. 2021;56(3):384–95.
Arora S, Guleria R. Oxidative stress and inflammation in critically ill patients on mechanical ventilation. Chest. 2009;136(4):31.
Santus P, Corsico A, Solidoro P, Braido F, Di Marco F, Scichilone N. Oxidative stress and respiratory system: pharmacological and clinical reappraisal of N-acetylcysteine. Copd-J Chronic Obstruct Pulmonary Dis. 2014;11(6):705–17.
Joelsson JP, Asbjarnarson A, Sigurdsson S, Kricker J, Valdimarsdottir B, Thorarinsdottir H, Starradottir E, Gudjonsson T, Ingthorsson S, Karason S. Ventilator-induced lung injury results in oxidative stress response and mitochondrial swelling in a mouse model. Lab Animal Res. 2022. https://doi.org/10.1186/s42826-022-00133-.
Ward PA. Oxidative stress: acute and progressive lung injury. In: Oxidative/Nitrosative Stress and Disease. Edited by Laskin DL, vol. 1203; 2010: 53–59.
Ruan H, Li YZ, Zhang Q, Wang BR, Wu R, Li SS, Ran X. Identification and clinical validation of hypoxia-inducible factor 1a protein as the potential biomarker in patients with sepsis. Shock. 2023;59(6):855–63.
Venkataramani AS, Bor J, Jena AB. Regression discontinuity designs in healthcare research. BMJ-British Med J. 2016. https://doi.org/10.1136/bmj.i1216.
Guo Z, Liu L, Yu F, Cai Y, Wang J, Gao Y, Ping Z. The causal association between body mass index and type 2 diabetes mellitus-evidence based on regression discontinuity design. Diabetes Metab Res Rev. 2021;37(8):3455.
Ghosh S, Dent R, Harper ME, Stuart J, McPherson R. Blood gene expression reveal pathway differences between diet-sensitive and resistant obese subjects prior to caloric restriction. Obesity. 2011;19(2):457–63.
Bermejo-Martin JF, Martin-Loeches I, Rello J, Anton A, Almansa R, Xu L, Lopez-Campos G, Pumarola T, Ran L, Ramirez P, et al. Host adaptive immunity deficiency in severe pandemic influenza. Critical Care. 2010;14(5):1–12.
Ritchie ME, Phipson B, Wu D, Hu YF, Law CW, Shi W, Smyth GK. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):13.
Sean D, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846–7.
Daszykowski M, Kaczmarek K, Heyden YV, Walczak B. Robust statistics in data analysis—A review basic concepts. Chemom Intell Lab Syst. 2007;85(2):203–19.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Yu GC, Wang LG, Han YY, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Carbon S, Douglass E, Good BM, Unni DR, Harris NL, Mungall CJ, Basu S, Chisholm RL, Dodson RJ, Hartline E, et al. The gene ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021;49(D1):D325–34.
Esteban A, Anzueto A, Frutos F, Alia I, Brochard L, Stewart TE, Benito S, Epstein SK, Apezteguia C, Nightingale P, et al. Characteristics and outcomes in adult patients receiving mechanical ventilation - A 28-day international study. Jama J Am Med Assoc. 2002;287(3):345–55.
Shen W, Song Z, Zhong X, Huang M, Shen D, Gao P, Qian X, Wang M, He X, Wang T, et al. Sangerbox: A comprehensive, interaction-friendly clinical bioinformatics analysis platform. Imeta. 2022;1(3):e36.
Ong CWM, Elkington PT, Brilha S, Ugarte-Gil C, Tome-Esteban MT, Tezera LB, Pabisiak PJ, Moores RC, Sathyamoorthy T, Patel V, et al. Neutrophil-Derived MMP-8 drives AMPK-dependent matrix destruction in human pulmonary tuberculosis. Plos Pathogens. 2015. https://doi.org/10.1371/journal.ppat.1004917.
Sirnio P, Tuomisto A, Tervahartiala T, Sorsa T, Klintrup K, Karhu T, Herzig K-H, Makela J, Karttunen TJ, Salo T, et al. High-serum MMP-8 levels are associated with decreased survival and systemic inflammation in colorectal cancer. Br J Cancer. 2018;119(2):213–9.
Berg G, Miksztowicz V, Schreier L. Metalloproteinases in metabolic syndrome. Clin Chim Acta. 2011;412(19–20):1731–9.
Juurikka K, Butler GS, Salo T, Nyberg P, Astrom P. The Role of MMP8 in Cancer: A Systematic Review. Int J Mol Sci. 2019;20(18):4506.
Van Lint P, Libert C. Matrix metalloproteinase-8: cleavage can be decisive. Cytokine Growth Factor Rev. 2006;17(4):217–23.
Ye S. Putative targeting of matrix metalloproteinase-8 in atherosclerosis. Pharmacol Ther. 2015;147:111–22.
Hoenderdos K, Lodge KM, Hirst RA, Chen C, Palazzo SGC, Emerenciana A, Summers C, Angyal A, Porter L, Juss JK, et al. Hypoxia upregulates neutrophil degranulation and potential for tissue injury. Thorax. 2016;71(11):1030–8.
Andrade VL, Petruceli E, Belo VA, Andrade-Fernandes CM, Caetano Russi CV, Bosco AA, Tanus-Santos JE, Sandrim VC. Evaluation of plasmatic MMP-8, MMP-9, TIMP-1 and MPO levels in obese and lean women. Clin Biochem. 2012;45(6):412–5.
Xu P, Gaertner F, Gihring A, Liu C, Burster T, Wabitsch M, Knippschild U, Paschke S. Influence of obesity on remodeling of lung tissue and organization of extracellular matrix after blunt thorax trauma. Respir Res. 2020;21(1):1–5.
Sepper R, Prikk K. Proteomics: Is it an approach to understand the progression of chronic lung disorders? J Proteome Res. 2004;3(2):277–81.
da Silva-Neto PV, do Valle VB, Fuzo CA, Fernandes TM, Toro DM, Fraga-Silva TFC, Basile PA, de Carvalho JCS, Pimentel VE, Perez MM et al: Matrix Metalloproteinases on Severe COVID-19 Lung Disease Pathogenesis: Cooperative Actions of MMP-8/MMP-2 Axis on Immune Response through HLA-G Shedding and Oxidative Stress. Biomolecules 2022. https://doi.org/10.3390/biom12050604.
This study was supported by funding from the National Natural Science Foundation of China (Grant No. 82271358), the Scientific Research Foundation for Returned Overseas Chinese Scholars of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, and the Talent Project of Public Health in Hubei Province (Grant No. 2022SCZ048).
Ethics approval and consent to participate
This study was approved by the ethics committee (TJ-IRB20230313).
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Differential expression of MMP-8 in diffrent diseases.
. The distribution of diseases among hospitalized patients.
. The relationship between MMP-8 expression levels and various clinical parameters.
. Appendix of Supplementary Table.
About this article
Cite this article
Ruan, H., Li, Ss., Zhang, Q. et al. Elevated MMP-8 levels, inversely associated with BMI, predict mortality in mechanically ventilated patients: an observational multicenter study. Crit Care 27, 290 (2023). https://doi.org/10.1186/s13054-023-04579-3