Abstract

Objectives. The purpose of this study was to determine whether there is a dose-response relationship between body mass index (BMI) and all-cause mortality in patients after coronary revascularization. Methods. The MIMIC-III database (version 1.4) was used as the sample population. For variables with less than 10% of values missing, we used the mice package of R software for multiple imputations. Cox regression was used to determine the risk factors of all-cause mortality in patients. RCSs were used to observe the relationship between BMI and all-cause mortality. Additional subgroup and sensitivity analyses were also performed to explore whether the conclusion can be applied to specific groups. Results. Both univariate and multivariate Cox models indicated that the mortality risk was lower for overweight patients than for normal-weight patients (). In RCS models, BMI had a U-shaped relationship with all-cause mortality of patients after coronary artery bypass grafting (CABG) ( for nonlinearity = 0.0028). There was a weak U-shaped relationship between BMI and all-cause mortality after percutaneous coronary intervention (PCI), but the nonlinear relationship between these two parameters was not significant ( for nonlinearity = 0.1756). Conclusions. The obesity paradox does exist in patients treated with CABG and PCI. RCS analysis indicated that there was a U-shaped relationship between BMI and all-cause mortality in patients after CABG. After sex stratification, the relationship between BMI and all-cause mortality in male patients who received PCI was L-shaped, while the nonlinear relationship among females was not significant.

1. Introduction

Coronary artery revascularization is an important method for treating coronary heart disease. Coronary artery bypass grafting (CABG) and percutaneous coronary intervention (PCI) are the two most common operations in coronary revascularization. The body mass index (BMI) is a prognostic factor for all-cause mortality in many diseases, with the prognosis being better in people with a normal BMI than in those who are overweight or obese [1, 2]. However, studies have indicated that some overweight or even obese patients have a better prognosis than those with a normal BMI, which is referred to as the obesity paradox [3, 4].

Gruberg et al. were the first to mention the obesity paradox in patients after revascularization [5]. That study found that patients with a normal BMI were found to have higher hospitalization and 1-year mortality rates than those who were overweight or obese. This indicates that being overweight or obese may protect patients and reduce mortality. However, some studies have attributed the obesity paradox to factors other than BMI [6, 7].

The relationship between BMI and all-cause mortality in patients after coronary revascularization remains unclear. Most studies have included BMI as a continuous or categorical variable when performing logistic or Cox regression analysis, which does illustrate changes in BMI effectively [810]. However, basic research is still needed to explore this relationship.

The purpose of this study was to determine whether there is a dose-response relationship between BMI and all-cause mortality in coronary heart disease patients undergoing coronary artery revascularization, thereby determining whether there is an obesity paradox in the prognosis of patients who underwent coronary artery revascularization. We hypothesized that it is a nonlinear relationship and used restricted cubic splines (RCSs) to determine the dose-response relationship.

2. Methods

2.1. Study Design and Population

The MIMIC-III database contains data on 53,423 adult patients (16 years or older) admitted to an intensive care unit between 2001 and 2012. The MIMIC-III database (version 1.4) was sampled, and the project was approved by the Institutional Review Boards of the Beth Israel Deaconess Medical Center (Boston, Massachusetts) and the Massachusetts Institute of Technology (Cambridge, Massachusetts). Consent was not required from individual patients since all of the protected health information of the project has been deidentified [11, 12]. We completed recognized courses for protecting human research participants, including the requirements of the HIPAA (Health Insurance Portability and Accountability Act), and signed a data usage agreement. The Institutional Review Boards of the Beth Israel Deaconess Medical Center and Massachusetts Institute of Technology have approved the use of the MIMIC-III database by any researcher who meets the data user requirements, and the requirement for patient’s informed consent is waived.

We used the ninth edition of the International Classification of Diseases (Clinical Modification) codes to identify and analyze all patients in the MIMIC-III database whose primary operation was revascularization (including CABG and PCI). BMI was calculated using the heights and weights recorded in the database. As the main study indicator, all missing values for BMI were deleted. The MIMIC-III database is connected to the social security database to record the follow-up times and outcomes of patients. Our study outcome was all-cause mortality as registered by the social security bureau of the patient. The follow-up times were reported in days.

We also extracted patient demographic, laboratory, and vital-sign indicators from the database. These factors were possible confounding factors in the relationship between BMI and all-cause mortality. All indicators were extracted from the diagnoses_icd, admissions, patients, icustays, labevents, and procedures_icd parameters in the database. The flowchart for data inclusion and exclusion is displayed in Figure 1.

The purpose of our research was to determine the dose-response relationship between BMI and all-cause mortality in patients after revascularization and to determine the existence of the obesity paradox while adjusting for possible confounding factors.

2.2. Statistical Analysis

All data were extracted using the Structured Query Language (SQL). Since BMI was the main study indicator, all missing values were deleted, and other possible confounding factors with >10% of values missing were also deleted. We used the mice package of R software for multiple imputation for variables with <10% of values missing [13]. This package uses relevant random samples based on the distribution of predicted values to replace the missing values with estimated values. We constructed a dataset by merging five imputation datasets.

The World Health Organization (WHO) standards were used to divide BMI into underweight, normal weight, overweight, and obese categories. Age was divided into three categories: youth, middle age, and elderly. All categorical variables were expressed as numbers and percentages, and the chi-squared and Fisher’s exact tests were used to determine the differences between the two groups. All continuous variables were expressed as medians and interquartile-range values, and the Mann–Whitney U test was used to identify differences between the two groups.

The Kaplan–Meier curve and log-rank test were used to evaluate whether there are differences in survival rates between different operations. Cox regression was used to determine the risk factors of patient all-cause mortality. We constructed univariate models, a demographic adjustment model (model I), and a multivariate adjustment model (model II). A trend test was used to determine linear trends between BMI and all-cause mortality.

RCSs were used to assess the nonlinear relationship between BMI and all-cause mortality. We divided the distribution of BMI into quartiles as the knots of the nonlinear models. We analyzed the nonlinear relationship between BMI and patient all-cause mortality and determined the value of the nonlinear test.

We also performed a subgroup analysis of sex and age to determine whether a relationship existed within a specific subgroup. We performed the following sensitivity analyses to determine the stability of our results: (1) we removed all patients who died within 30 days of their operation and performed the RCS analysis again, (2) abnormal BMI values (more than three times the standard deviation) were deleted before the analysis, and (3) we deleted the missing values from the original data before performing multiple imputations and reanalyzed the data to explore the difference resulting from multiple imputations.

All statistical analyses were conducted using Navicat Premium and R software (version 3.6.2, The R Foundation for Statistical Computing, Vienna, Austria). All cited values were two-sided, and was considered statistically significant.

2.3. Patient and Public Involvement

This research was done without patient involvement. Patients were not invited to comment on the study design and were not consulted to develop patient-relevant outcomes or interpret the results. Patients were not invited to contribute to the writing or editing of this document for readability or accuracy.

3. Results

3.1. Population Characteristics

Our final sample comprised 4133 patients treated with coronary artery revascularization in the MIMIC-III database: CABG patients (86.9%, n = 3593) and PCI patients (13.1%, n = 540). Their baseline data are listed in Table 1, which indicates that most of the variables differed significantly between the two types of surgery (). Table 1 indicates that there was a very small number of underweight patients, with overweight patients comprising the largest proportion (approximately 39%). Most patients were male, married, and white. Figure 2 displays the survival curves after surgery for both types of patients. The figure suggests that there are significant statistical differences between the two, and so the two types of patients were analyzed separately.

3.2. Univariate and Multivariate Cox Regression Analyses

We first divided BMI into four categories based on the WHO standards and then performed univariate and multivariate Cox regression analyses with normal weight as the reference. After adjusting all of the variables, the Cox regression model showed that length of stay, Elixhauser Comorbidity Index (ECI), urine output, pH, potassium, prothrombin time, red blood cell distribution width (RDW), age, marital status, and BMI were the significant prognostic factors for CAGB patients, while ventilator support, ECI, urine output, heart rate, temperature, hematocrit, RDW, sex, and BMI were the significant prognostic factors for PCI patients. The Cox regression results are presented in detail in Supplementary Tables S1 and S2.

Tables 2 and 3 present the results for model I (with the demographic characteristics of age, sex, race, and marital status adjusted) and model II (with the risk factors derived from the multivariate Cox analysis adjusted). Table 2 lists the hazard ratio (HR) and 95% confidence interval (CI) values for CABG patients according to their BMI group. Both the univariate model and model II showed that the risk of mortality was lower for overweight patients than for normal-weight patients (). Table 3 lists the HRs and 95% CIs for PCI patients. All three models indicated that the mortality risk was lower for overweight PCI patients than for normal-weight patients (). The trend tests of the six models revealed that a significant linear trend was only present for the BMI univariate model of PCI patients (), which also suggested the presence of a nonlinear relationship between BMI and all-cause mortality.

3.3. RCS Analyses of Nonlinear Relationships

Supplementary Tables S3 and S4 list the HR and 95% CIs of the BMI quartiles for all-cause mortality. The first quartile (Q1) was taken as a reference group for comparison with the other groups to obtain the corresponding HRs. In CABG patients, compared with the Q1 of BMI, the HR for Q2 of BMI was 0.722 (95% CI = 0.586–0.888) for the univariate model, 0.759 (95% CI = 0.616–0.936) for model I, and 0.740 (95% CI = 0.599–0.913) for model II (all ). Unlike CABG patients, the Q2 group of PCI patients did not show significant differences, whereas the Q3 group did show significant differences. Compared with the Q1 of BMI, the HR for Q3 of BMI was 0.506 (95% CI = 0.332–0.773) for the univariate model, 0.572 (95% CI = 0.370–0.884) for model I, and 0.601 (95% CI = 0.392–0.924) for model II (all ).

We constructed RCS models to further analyze the relationship between BMI and patient mortality. Figure 3 shows the dose-response curves of BMI and all-cause mortality after adjusting for the factors that were significant in the Cox analysis. The dose-response analysis revealed a U-shaped curve between BMI and the risk of all-cause mortality in CABG patients ( for nonlinearity = 0.0028) (Figure 3(a)). In contrast, there was a weak U-shaped relationship between BMI and all-cause mortality after PCI, but no significant nonlinear relationship between these two variables ( for nonlinearity = 0.1756) (Figure 3(b)).

3.4. Subgroup Analyses

After adjusting for the corresponding confounding factors, the subgroup analysis of BMI and all-cause mortality of patients after CABG (Supplementary Table S5) showed that there was no significant interaction effect of sex ( for interaction = 0.330) or age ( for interaction = 0.883). This indicates that the relationship between BMI and CABG did not differ significantly with sex or age. However, Supplementary Table S6 showed that significant interaction was found in the PCI group for stratification according to sex ( for interaction = 0.006). This suggests that the impact of BMI on the mortality of PCI patients is affected by sex.

We, therefore, constructed RCS models to perform a sex-stratified analysis. As shown in Figure 4, the dose-response relationship appeared as U-shaped curves for all-cause mortality of BMI and CABG. However, the nonlinear trend was significant for males ( for nonlinearity = 0.0163) but not for females ( for nonlinearity = 0.1367). The sex-related difference appeared to be greater in PCI patients. There was an approximate L-shaped relationship between BMI and all-cause mortality in male patients ( for nonlinearity = 0.0085) but not in female patients ( for nonlinearity = 0.8574).

3.5. Sensitivity Analysis

After excluding patients who died within 30 days after the operation or had abnormal BMI values, the RCS analysis showed that the nonlinear relationship between BMI and all-cause mortality remained consistent with the previous results. That is, there was a significant U-shaped nonlinear relationship in CABG patients, but no significant nonlinear relationship in PCI patients. However, analyzing the complete data before performing multiple imputations revealed a U-shaped relationship between BMI and all-cause mortality for both CABG and PCI patients ( for nonlinearity <0.05) (Figure 5). It should be noted that there was no significant nonlinear relationship for PCI patients after performing multiple imputation.

4. Discussion

Higher BMI is usually a risk factor for adverse outcomes of cardiovascular disease and various complications [14]. However, many studies have found that overweight or obese people have a survival advantage compared to people who are overly thin or have a BMI within the normal range, which is called the obesity paradox [15, 16]. There have been some reports on the obesity paradox in patients after revascularization. However, most of these studies simply divided BMI into four categories according to international standards and segmented the relationship between BMI and outcome [17, 18]. In contrast, the present study used RCSs to explore the nonlinear relationship between BMI and all-cause mortality after coronary revascularization surgery and, thereby, expressed the data using continuous and smooth graphs, which is more intuitive for explaining the overall trend for ORs of BMI than using the traditional segmented analysis.

The present results show that the survival rate is higher for CABG patients than for PCI patients. In univariate and multivariate Cox regression analyses, BMI was a prognostic factor for all-cause mortality in patients regardless of whether they received CABG or PCI. In the RCS analysis, the models before and after adjustment all showed a U-shaped relationship between the BMI and all-cause mortality for CABG patients. For PCI patients, there was a weak U-shaped relationship between these variables, but no significant nonlinear relationship between them after performing multivariate adjustment.

It is especially noteworthy that our subgroup analysis of PCI patients revealed a significant interaction between sex and BMI, which means that sex may affect the relationship between BMI and mortality. There was no significant difference in the Cox regression results for the relationship between BMI and all-cause mortality among females. Correspondingly, in the RCS analysis, we only observed significant nonlinear relationships in males, which indicates that the U-shaped relationship may only apply to males and that further investigations are needed into the relationship for females.

Previous studies have found BMI to be a poor discriminator of females with different risks of coronary heart disease. Weight and BMI are now considered less important than previously thought [19]. Materko et al. showed that the existing BMI classification standards are not suitable for females, [20] which might explain why the relationship between BMI and all-cause mortality in females was not significant in the present study.

In addition, after removing all missing data without performing multiple imputations in the sensitivity analysis, the nonlinear relationship became significant in PCI patients. This result might have been due to there being too many missing values for PCI patients since only 97 patients remained after removing the missing values. Compared with the 3345 CABG patients, the much smaller sample for PCI patients might have introduced bias into the study. Therefore, the relationship between BMI and all-cause mortality in PCI patients needs to be investigated further with large samples in prospective cohorts. The results of the other sensitivity analysis are consistent with the results obtained when analyzing the full sample, which further confirms the present research conclusions.

The present research has yielded further evidence for the existence of the obesity paradox. However, unlike previous studies [2123], we used RCSs to obtain HR values, and hence, more-precise risk ranges. In CABG patients, taking a BMI of 25 kg/m2 as the reference value, patients with a BMI within the range of 25–32 kg/m2 had a lower HR, which corresponds to a higher survival rate. The reasons for this obesity paradox remain to be elucidated, but some possible explanations are given below. One hypothesis is that obese patients receive more aggressive drug treatment than normal-weight patients, resulting in better control of their various indicators [24, 25]. In addition, it has been reported that compared with normal-weight patients, obese hypertensives have a higher cardiac output, enlarged blood volume, and lower systemic vascular resistance. This means that, for a given level of arterial pressure, obese patients will have lower total peripheral resistance than lean patients, which could improve the survival rate of obese patients [26, 27]. It has also been shown that BMI is a poor indicator for categorizing obesity. Indicators such as waist circumference and waist-to-height ratio may better reflect the true degree of obesity. [28, 29].

This study was subject to some limitations. Firstly, the MIMIC-III database is from a single center, and so the present research results need to be verified in a multicenter population. Secondly, although we used Cox regression for analysis, cross-sectional studies cannot be used to determine causality, and so future prospective cohort studies are needed. Thirdly, the large number of missing values for PCI patients may have made the results unstable. Larger samples are needed to confirm the conclusions and sex differences for PCI patients.

5. Conclusion

The analysis performed in this study of the MIMIC-III database revealed that the obesity paradox does exist in patients treated with CABG and PCI. The RCS analysis showed that BMI has a U-shaped relationship with all-cause mortality in CABG patients, while this relationship is not significant in PCI patients. However, sex stratification revealed an L-shaped relationship between BMI and all-cause mortality in male PCI patients, while there was no significant nonlinear relationship in females. Multicenter and large-sample prospective studies are needed to further explore sex differences in the obesity paradox.

Abbreviations

BMI:Body mass index
CABG:Coronary artery bypass grafting
PCI:Percutaneous coronary intervention
RCS:Restricted cubic splines
HIPAA:Health Insurance Portability and Accountability Act
SQL:Structured Query Language
WHO:World Health Organization
ECI:Elixhauser Comorbidity Index
RDW:Red blood cell distribution width
HR:Hazard ratio
CI:Confidence interval.

Data Availability

Data were fully available at https://physionet.org/content/mimiciii-demo/1.4/.

Ethical Approval

The MIMIC-III database has received ethical approval from the Institutional Review Boards of both Beth Israel Deaconess Medical Center (Boston, MA, USA) and the Massachusetts Institute of Technology (Cambridge, MA, USA).

Informed consent was unnecessary because the MIMIC-III data are anonymous and publicly available.

Conflicts of Interest

The authors declare that they have no conflicts of interests.

Authors’ Contributions

Chengzhuo Li and Didi Han contributed equally to this work.

Acknowledgments

This study was supported by The National Social Science Foundation of China (grant/award no. 16BGL183).

Supplementary Materials

Table S1: selected variables by multivariate Cox regression analysis for CAGB patients. Table S2: selected variables by multivariate Cox regression analysis for PCI patients. Table S3: Cox regression analyses of the relationship between BMI quartiles and all-cause mortality for patients after CABG. Table S4: Cox regression analyses of the relationship between BMI quartiles and all-cause mortality for patients after PCI. Tables S5: subgroup analysis of the associations between BMI and all-cause mortality for patients after CABG. Tables S6: subgroup analysis of the associations between BMI and all-cause mortality of patients after PCI. Supplementary material related to this article can be found in Supplementary Tables S1-S6. (Supplementary Materials)