The Association Between Congestive Heart Failure and One-Year Mortality After Surgery in Singaporean Adults: A Secondary Retrospective Cohort Study Using Propensity-Score Matching, Propensity Adjustment, and Propensity-Based Weighting

Background: Although congestive heart failure (CHF) is considered a risk factor for postoperative mortality, reliable quantification of the relationship between CHF and postoperative mortality risk is limited. We aimed to investigate the association between CHF and 1-year mortality after surgery in a large cohort of the Singaporean population.

Methods: In this retrospective cohort study, the study population included 69,032 adult patients who underwent surgery at Singapore General Hospital between 1 January 2012 and 31 October 2016. The target independent and dependent variables were CHF and 1-year mortality after surgery, respectively. Propensity score was estimated using a non-parsimonious multivariable logistic regression model. Multivariable adjustment, propensity score matching, propensity score adjustment, and propensity score-based weighting Cox proportional-hazards regression were performed to investigate the association between CHF and 1-year mortality after surgery.

Results: The multivariate-adjusted hazard ratio (HR) in the original cohort was 1.39 (95% confidence interval (CI): 1.20–1.61, P < 0.001). In additional propensity score adjustment, the HR between CHF and 1-year mortality after surgery was 1.34 (95% CI: 1.15–1.56, P < 0.001). In the propensity score-matched cohort, the multivariate-adjusted Cox proportional hazard regression model analysis showed participants with CHF had a 54% increased risk of 1-year mortality after surgery (HR 1.54, 95% CI: 1.19–1.98, P < 0.001). The multivariate-adjusted HR of the inverse probability of treatment-weighted and standardised mortality ratio-weighted cohorts was 1.34 (95% CI: 1.10–1.62, P = 0.004) and 1.24 (95% CI: 1.17–1.32, P < 0.001), respectively.

Conclusion: CHF is an independent risk factor for 1-year mortality after surgery in patients undergoing surgery. Depending on the statistical method, patients with CHF had a 24–54% increased risk of 1-year all-cause mortality after surgery. This provides a reference for optimising clinical decision-making, improving preoperative consultation, and promoting clinical communication.

Introduction

Congestive heart failure (CHF) is the leading cause of morbidity and mortality worldwide (1). Its prevalence is more than 22 million cases worldwide, with the incidence of 2 million new cases per year (2). The 5-year mortality rate is 62% for women and 75% for men after CHF (3). With advancements in medical care, an increasing number of patients with CHF may select surgery to treat certain diseases. Therefore, it is imperative to improve our understanding of the impact of HF on postoperative outcomes. Although CHF has been recognised as a risk factor for postoperative mortality and has been incorporated into several risk indicators (4–8), few studies have examined the relationship between CHF and postoperative mortality. The hazard ratio (HR)/odds ratio (OR) produced by previous studies on the relationship between CHF and the prognosis of surgical patients varied widely (9–11).

Previous studies have mainly applied traditional parsimonious regression models based on analytical adjustments to control for confounders. However, such models may still result in a bias owing to unmeasured or residual confounding of the model and overfitting (12, 13). Research methods based on propensity score (PS) are considered the core alternative for controlling the confounding of observational research. Both large and small sample theories show that adjustment for the scalar PS is sufficient to remove bias due to all observed covariates (14, 15). Several adjustment methods incorporating the estimated PS have been proposed, including matching, regression adjustment, and weighting (13, 15–18).

Therefore, based on the current status of research on the impact of CHF on the prognosis of surgical patients, no studies have applied PS-based methods to investigate this relationship. We used a large sample of patients discharged from a general hospital in Singapore to study the differential impact of CHF on surgical outcomes using various statistical models.

Materials and Methods

Data Source and Participants

This was a post-mortem analysis of a large vertical cohort established by Yilin Eileen Sim’s team in Singapore. The analysed data were stored in the Dryad database¹ by Yilin et al. The author of the original study waived all copyrights and related ownership of these data. Therefore, we can use these data for secondary analysis without infringing on the rights of the author. The data were obtained from a published paper: Sim et al. (19).

These clinical records were obtained from the hospital’s clinical information system [Sunrise Clinical Manager (SCM); Allscripts, IL, United States]. The mortality and follow-up data of the original study were paired with their national electronic health records. The variables included in the database file were as follows: sex, age, race, history of previous cerebrovascular accidents (CVA), history of ischaemic heart disease (IHD), history of diabetes mellitus on insulin (DMI), history of CHF, creatinine category, priority of surgery, surgical risk classification, American Society of Anaesthesiologists Physical Status (ASA-PS), red cell distribution width (RDW) category, stage of chronic kidney disease (CKD), degree of anaemia, type of anaesthesia (general or regional anaesthesia), follow-up days, and survival status (19, 20). The original study was approved by the Institutional Review Board (SingHealth CIRB 2014/651/D), which waived the requirement for individual informed consent (19). As reported elsewhere, the current study is a secondary analysis of the original study and did not receive ethical approval (21).

Study Sample

Consistent with the original study, patients who underwent surgery at Singapore General Hospital from 1 January 2012 to 31 October 2016 and were aged ≥18 years were included. The original research exclusion criteria were as follows: (1) patients who underwent transplantation and burn surgeries, (2) patients who had no surgery, and (3) patients who underwent minor surgeries (19). The final dataset consisted of 97,443 participants. In the present study, we excluded participants with missing CHF information (n = 28,411). Figure 1 depicts the participant selection process. Finally, our study included a total of 69,032 participants for secondary analysis.

FIGURE 1

Figure 1. Flowchart of study participants. It showed the process of screening participants.

Variable Definitions

Laboratory examinations were performed within 90 days of surgery. These included serum creatinine level, serum haemoglobin, and RDW. Anaemia was defined by the World Health Organization gender-based classification of anaemia severity (19, 22). The severity of anaemia is characterised as follows: mild anaemia (haemoglobin concentration: male 11–12.9 g/dL, female 11–11.9 g/dL); moderate anaemia (haemoglobin concentration: 8.0–10.9 g/dL), and severe anaemia (haemoglobin concentration: 8.0 g/dL). RDW was reported as a coefficient of variation (percentage) of red blood cell volume; levels above 15.7% were defined as high RDW, and the normal reference range for RDW was 10.9–15.7%. The target independent variable was the presence or absence of CHF history obtained at baseline. The dependent variable was 1-year mortality event during follow-up. Of the study patients, 88.7% were followed up for 1 year and mean duration of follow-up was 258 days (19).

Missing Data Processing

Missing data in observational studies are a frequently encountered problem that cannot be fully avoided. Missing data accounted for 2.71% of all variables in the dataset analysed in this study (Supplementary Table 1). To reduce the deviation caused by missing covariates, which cannot reflect the statistical efficiency of the target sample in the modelling process, missing data in this study adopted multiple imputations (23, 24). In the present study, the analysis steps in the original and imputed datasets were calculated and compared. Similar core results were obtained using the original and imputed datasets; thus, we report the results for the imputed dataset (Supplementary Table 2).

Statistical Analysis

Participants were stratified by CHF, continuous variables are expressed as median (quartile) (skewed distribution) or mean ± standard deviation (normal distribution), and categorical variables are expressed as a frequency or percentage. Two-sample t-tests were performed for normally distributed continuous variables, Wilcoxon rank-sum tests for non-normally distributed continuous variables and ordered categorical variables, and chi-square tests for categorical variables.

The study goal was to use different methods to control for confounding factors and evaluate the impact of CHF on the 1-year postoperative mortality of patients who underwent surgery. Specifically, we applied five methods: the original cohort multivariate-adjusted Cox proportional-hazards regression model and four PS methods. The four PS methods included PS matching multivariate-adjusted Cox proportional-hazards regression model; PS adjustment Cox proportional-hazards regression model; inverse probability of treatment-weighted (IPTW) model, and standardised mortality ratio (SMR)-weighted multivariate-adjusted Cox proportional-hazards regression model. (1) PS was estimated using a non-parsimonious multivariable logistic regression model wherein CHF was used as the independent variable, and all baseline characteristics listed in Table 1 were covariates (25). The variables used for matching included age, sex, race, CVA, IHD, DMI, stage of CKD, degree of anaemia, priority of surgery and surgical risk classification, ASA-PS, RDW category, type of anaesthesia, matching using a 1:1 matching protocol, no replacement (greedy matching algorithm), and calliper width of 0.0005. Standardised differences (SDs) were estimated for all baseline covariates before and after matching to assess pre-matched and post-matched balance. SD of <20.0% for a given covariate indicate a relatively small imbalance (26–28). We used a multivariate-adjusted Cox proportional hazard regression model to evaluate the association between CHF and 1-year mortality after surgery in the PS-matched cohort. (2) PS adjustment: a multivariable Cox proportional-hazards model with the same strata and covariates, with additional adjustment for PS. The analysis included all patients (29–31). (3) The IPTW estimator estimates the treatment effect in a population whose risk factor distribution is equal to that found in all study subjects. IPTW was calculated as the inverse of the PS of patients with combined CHF and the inverse of (1 - PS) for non-CHF patients. The IPTW model was applied to generate a weighted cohort (15, 31, 32). IPTW multivariate-adjusted Cox proportional-hazards regression model has the same strata and covariates with inverse probability weighting according to PS. The analysis included all patients (31). (4) The SMR-weighted estimator estimates the treatment effect in a population whose distribution of risk factors is equal to that found in the treated study subjects only. SMR-weighted analyses used 1 as the value of the CHF group and the probability of propensity PS/(1 - PS) for the non-CHF group as weights and estimated the standardised effect measure for the CHF group (exposed group) as the standard population (15, 18). SMR-weighted multivariate-adjusted Cox proportional-hazards regression model, with the same strata and covariates with inverse probability weighting according to PS. All patients were included in the analysis. Cumulative mortality was used to describe mortality (15, 33). When added to the model, the covariates that changed the hazard ratios by at least 10% were considered confounders and were adjusted for in the multivariate analysis (34). In addition, Kaplan–Meier analysis and log-rank tests were performed to compare 1-year mortality after surgery.

TABLE 1

Table 1. Baseline characteristics before and after propensity score matching.

We performed a series of sensitivity analyses to test the robustness of the results; first, because of the extreme difference in the mortality OR among the low- and high-propensity strata. Mortality was lower in patients with low PS. Therefore, we excluded participants with a PS of <0.05 and performed a sensitivity analysis using five models. Second, bias may arise because substantial CHF information was missing. We applied multiple imputations to estimate missing CHF values (n = 28,411). The imputation model included age, sex, race, CHF, CVA, IHD, DMI, stage of CKD, degree of anaemia, priority of surgery and surgical risk classification, ASA-PS, RDW category, and type of anaesthesia. Multiple statistical models were used to analyse the association between CHF and postoperative mortality in patients to verify the reliability of our analysis after excluding participants with missing CHF information. In addition, we explored the possibility of unmeasurable confounding between CHF and 1-year mortality after surgery by calculating E-values (35).

Differences in age and CKD stage remained between the CHF and non-CHF groups after PS matching. This may be associated with increased mortality in surgical patients (36). This may lead to overestimation of the relationship between CHF and 1-year mortality after surgery. Therefore, we stratified the participants according to CKD (stages 1–2 vs. 3–5) (37, 38). Pre-specified subgroup analyses were performed based on these two characteristics. The subgroups were based on age and CKD stage. Each stratification adjusted for all factors except the stratification factor itself. Only the corresponding matched pairs in the same subgroup were selected to maintain the balance of baseline characteristics between the CHF and non-CHF groups in subgroup analyses.

All results are reported according to the STROBE statement (39). All analyses were performed using the statistical software packages R² (The R Foundation) and Empower-Stats³ (X&Y Solutions, Inc., Boston, MA, United States). All tests were two-tailed, and P-values of <0.05 were regarded as statistically significant.

Results

A total of 69,032 participants (48.32% men and 51.68% women) were included in the analysis. Among them, 1,413 (2.05%) had CHF and 67,619 (97.95%) did not have CHF. The number of participants aged 18–29, 30–49, 50–69, and ≥70 was 7,288 (10.562%), 19,174 (27.78%), 30,269 (43.85%) and 12,301 (17.82%), respectively. PS was estimated using a non-parsimonious multivariable logistic regression model, and CHF was used as the independent variable. The variables used for matching included age, sex, race, CVA, IHD, DMI, stage of CKD, degree of anaemia, priority of surgery and surgical risk classification, ASA-PS, RDW category, type of anaesthesia, matching using a 1:1 matching protocol, no replacement (greedy matching algorithm), and calliper width of 0.0005. Before PS matching, there were differences in almost all baseline characteristics between the CHF and non-CHF groups (Table 1). After one-to-one matching using PS analysis, 1,315 patients with CHF were successfully matched with 1,315 non-CHF participants. After matching, except race and age, the SD for almost all variables were <20.0%, indicating that the PSs were well matched. In other words, there was only a small difference in baseline characteristics between the CHF and non-CHF groups (Table 1). The logistic model used to estimate PS yielded a c-statistic of 0.811 (Figure 2).

FIGURE 2

Figure 2. The ROC curve of propensity score to predict one-year mortality after surgery. It showed that the logistic model used to estimate the propensity score yielded a c-statistic of 0.811.

Supplementary Figure 1 showed that in the original cohort, the mean PS was 0.017 ± 0.049 for the CHF group and 0.163 ± 0.159 for the non-CHF group (P < 0.01). However, the mean PS in the matched population was 0.166 ± 0.125 for both groups (P = 0.995). (Supplementary Figure 2). The probability density functions of the PS for the CHF and non-CHF groups are summarised in Figure 3. As expected, the distribution of PS for the non-CHF group shifted toward 0 and for the CHF group shifted somewhat toward 1. The figure also illustrates that the overlap of PS for the CHF and non-CHF groups is limited to a narrow range (Figure 3A and Supplementary Figure 3A). The probability density functions of PS for the CHF and non-CHF groups after matching are summarised in Figure 3B. As expected, the distribution of PS for the CHF and non-CHF groups remained basically the same (Supplementary Figure 3B).

FIGURE 3

Figure 3. The probability density functions of the propensity score for CHF and non-CHF participants (A) before matching and (B) after matching. It showed the probability density functions of the propensity score for congestive heart failure and non-congestive heart failure participants before matching and after matching.

Propensity Stratum-Specific Effects

For 1-year mortality after surgery, the gradients across levels of PS for the CHF and non-CHF groups were strong and unexpectedly different. Supplementary Table 3 summarises the proportions of patient mortality during the follow-up period in the CHF and non-CHF groups according to PS percentiles. We noted some important considerations. First, the PS of very few patients in the CHF group was below the 70th percentile of the overall PS. Second, the mortality rate of both groups increased as PS increased. The associated empirical OR for 1-year mortality after surgery increased from 0.821 in the 70th–80th percentiles of PS to 1.247 in the 99th percentile (Supplementary Table 3).

One-Year Mortality After Surgery

Supplementary Table 4 showed the 1-year mortality after surgery of the CHF and non-CHF groups before and after PS matching. Before PS matching, 2,307 participants died during the follow-up period. The corresponding mortality rates in the CHF and non-CHF groups were 15.9% (95% confidence interval (CI): 14.0–17.8) and 3.1% (95% CI: 2.9–3.2), respectively. After PS matching, the difference in mortality between the two groups changed significantly; the corresponding mortality rates in the CHF and non-CHF groups were 14.8% (95% CI: 12.8–16.7) and 10.0% (95% CI: 8.4–11.7). Kaplan–Meier analysis showed that the CHF group had a higher 1-year mortality after surgery than the non-CHF group in the original cohort (P < 0.001). After PS matching, the difference in mortality between the two groups significantly reduced (Figure 4).

FIGURE 4

Figure 4. (A) Kaplan–Meier survival curve based on congestive heart failure in the original cohort. Kaplan–Meier analysis of one-year mortality after surgery based on congestive heart failure(CHF) and non-congestive heart failure (non-CHF) in the original cohort(log-rank, P < 0.0001). (B) Kaplan– Meier survival curve based on congestive heart failure in the propensity-score matching cohort. Kaplan–Meier analysis of one-year mortality after surgery based on congestive heart failure (CHF) and non-congestive heart failure (non-CHF) in the propensity-score matching cohort (log-rank, P = 0.0002).

Analysis Results Through Different Confounding Control Methods

We used the Cox proportional-hazards regression model to assess the association between CHF and 1-year mortality after surgery in the original, PS matching, and weighted cohorts. In the original cohort, CHF was significantly associated with 1-year mortality after surgery (HR, 5.65; 95% CI: 4.92–6.48, P < 0.001). In other words, compared with participants without CHF, those with CHF had a 4.65-fold increased risk of 1-year mortality after surgery. After multivariate adjustment (adjusted for sex, age, race, CVA, IHD, DMI, priority of surgery, surgical risk classification, ASA-PS, stage of CKD, degree of anaemia, and type of anaesthesia), the association still existed (HR: 1.39, 95% CI: 1.20–1.61, P < 0.001). In PS adjustment (adjusted variables include PSs and the same strata and covariates), the HR between CHF and 1-year mortality after surgery was 1.34 (95% CI: 1.15–1.56, P < 0.001). Second, in the PS-matched cohort, the multivariate-adjusted Cox proportional hazard regression model analysis showed that the HR between CHF and 1-year mortality after surgery was 1.54 (95% CI: 1.19–1.98, P < 0.001). The adjusted variables were the same as those used in the original cohort. Finally, in the weighted cohort, after multivariate adjustment, the SMR-weighted analysis yielded an HR of 1.34 (95% CI: 1.10–1.62, P < 0.001) and the result of the IPTW multivariate-adjusted Cox proportional-hazards regression model analysis showed that the HR between CHF and 1-year mortality after surgery was 1.24 (95% CI: 1.17–1.32, P < 0.001). It should be emphasised that in the original cohort multivariate adjustment, PS matching multivariate adjustment, IPTW, and SMR-weighted multivariate adjustment, the adjusted variables were the same, including sex, age, race, CVA, IHD, DMI, priority of surgery, surgical risk classification, ASA-PS, stage of CKD, degree of anaemia, and type of anaesthesia. In PS adjustment, the adjusted variables included other model-adjusted variables and PS (Table 2).

TABLE 2

Table 2. Associations between CHF and one-year postoperative mortality of surgical patients in the crude analysis, multivariable analysis, and four propensity-score methods analyses.

Sensitivity Analysis

We considered a significant difference in the associated empirical OR for 1-year mortality after surgery between the participants with low and high PSs (Supplementary Table 3). We further analysed patients with PS of ≥0.05. The crude HR for the restricted population was 1.81 (95% CI: 1.55–2.12, P < 0.001); the HRs for the five different methods were similar. The multivariate-adjusted HR in the original cohort was 1.59 (95% CI: 1.35–1.89, P < 0.001); PS adjustment revealed that the HR between CHF and 1-year mortality after surgery was 1.54 (95% CI: 1.30–1.85, P < 0.001). In the PS-matched cohort, after multivariate adjustment, HR was 1.55 (95% CI: 1.19–2.02, P < 0.001); the multivariate-adjusted HR of the IPTW and SMR-weighted cohorts was 1.49 (95% CI: 1.36–1.67, P < 0.001) and 1.45 (95% CI: 1.18–1.78, P < 0.001), respectively (Supplementary Table 5).

In addition, bias may arise due to the excessive missing CHF information. We applied multiple imputations to estimate missing CHF values (n = 28,411). After imputation, there were 2,101 participants in the CHF group and 95,342 in the non-CHF group. We applied five statistical models to the analysis, which yielded similar results (Supplementary Table 6). The HRs of the original cohort-adjusted model, PS adjustment model, PS matching adjusted, IPTW model, and SMR-weighted model were 1.24, 1.26, 1.31, 1.16, and 1.31, respectively. Therefore, excluding participants with missing CHF information did not affect the core findings of this study and suggests that our results are robust.

Furthermore, we generated an E-value to assess the sensitivity to unmeasured confounding factors. The E-value (2.13) was lower than the relative risk (3.34) of unmeasured confounders and 1-year mortality after surgery, suggesting that unmeasured or unknown confounders had little effect on the relationship between CHF and 1-year mortality after surgery. Our primary findings were robust.

Subgroup Analysis

We performed subgroup analysis to examine the impact of potential confounders that may influence the association between CHF and 1-year mortality after surgery. Age and CKD stage were used as stratified variables to assess the trend of effect size. We observed interactions among the subgroups according to our specification. We found that age and CKD stage did not affect the relationship between 1-year mortality after surgery (all P-values for interaction < 0.05), respectively (Supplementary Table 7).

Discussion

This one-to-one PS-matched cohort study showed that CHF was associated with a higher risk of mortality 1 year after surgery. After PS matching, CHF had a significant association with 1-year mortality after surgery, and the risk of mortality increased by 54% in the population with CHF (HR = 1.54, 95% CI: 1.1–1.98, P < 0.001). In PS adjustment, the HR between CHF and 1-year mortality after surgery was 1.34 (95% CI: 1.15–1.58, P < 0.001). We applied a Cox proportional-hazards model based on PS-based weighting to further verify the association between CHF and 1-year mortality after surgery. In the IPTW and SMR-weighted cohorts, compared with participants without CHF, those with CHF had a 24% and 34% increase in the risk of mortality 1 year after surgery, respectively.

We found that different methods of controlling confounding factors resulted in different HRs. The results of the PS matching multivariate-adjusted Cox proportional-hazards regression model were slightly higher, whereas the original cohort multivariate-adjusted, PS-adjusted, and SMR-weighted multivariate-adjusted Cox proportional-hazards regression models had similar hazard ratios. In comparison, the results of the IPTW multivariate-adjusted Cox proportional-hazards regression model were lower (Table 2). The number of participants without CHF was many times greater than the number of those with CHF. PS matching usually results in a successful match for almost all patients with CHF; many unmatched patients without CHF were excluded from the analysis. As a result, the distribution of covariates in the (successfully) matched subpopulation would be close to that in the treated study population. Most patients in the CHF group were in the propensity strata with a high risk of associated mortality, and the SMR-weighted method estimated the average effect of CHF in a population whose distribution of risk factors was equal to that found in the CHF group. Thus, it was not surprising that the SMR-weighted estimate was closer to the PS-matched estimate.

In contrast, the IPTW model estimated the average CHF effect for the entire study population. Given that 90% of the study population was in the propensity strata associated with low empirical odds ratios, the results of the IPTW multivariate-adjusted Cox proportional-hazards regression model were lower. The IPTW estimate increased to 1.49 when the patients in these three strata were excluded by restricting the analysis to the subpopulation of the CHF and non-CHF groups whose PS was ≥0.05. Indeed, once we restricted the analysis to subjects with PS of ≥0.05, all adjustment methods provided approximate results, with HR fluctuating around 1.50. The results of all these methods showed that the four PS methods could control for confounding factors well.

CHF is an established risk factor for poor prognosis after surgery across a wide range of specialties. A study conducted in Sweden showed that the crude and adjusted HRs for 30-day mortality after elective surgery in patients with CHF were 5.36 and 1.79, respectively (9). An analysis of 21,560,996 surgical hospitalisations revealed that the adjusted OR for in-hospital perioperative mortality in patients was 2.15 (11). Another study found similar results; compared with patients without CHF, patients with CHF had a 96% increased risk of 30-day mortality (10). Our study complements the existing literature supporting the hypothesis that CHF increases the risk of postoperative mortality in patients undergoing surgery. However, the estimated HR for the relationship between CHF and postoperative mortality was lower than that reported in previous studies.

We analysed these inconsistent findings, which may be justified by the following possible explanations: (1) The research population was different, including age, sex and race. As the original data failed to define the surgical category, our study population included all cardiac and non-cardiac surgery patients. Other studies have generally analysed these two populations separately. (2) Sample sizes in these studies varied widely. (3) The studies were adjusted for different covariates. (4) Previous studies mainly used variable adjustments to control confounding factors; this traditional parsimonious regression model could result in a bias because of unmeasured or residual confounding or overfitting of the model (40). We used PS methods to control for confounding factors and verify the association between CHF and 1-year mortality after surgery. (5) The outcome variables were different. Previous studies have focused on in-hospital or 30-day mortality rates. However, the dependent variable in our study was 1-year mortality after surgery.

Although the exact mechanism underlying the increased risk of postoperative death associated with CHF is unclear, the incidence of complications in this population may provide clues. A study found that, compared with patients without CHF, the incidence of crude complications in patients with CHF doubled (41). There is convincing evidence that this may be because patients with CHF have poor recovery, and even minor postoperative complications can significantly affect postoperative mortality (42).

Our study has two other strengths: (1) Our sample size was relatively larger than that of previous similar studies. (2) To the best of our knowledge, few cohort studies have used PS matching to explore the relationship between CHF and 1-year mortality after surgery. Research methods based on PSs are considered a core alternative for controlling the confounding effects of observational research.

The potential limitations of this study are as follows. First, the population included in this study was Singaporean and the race was mainly Chinese. Therefore, the universality of these results in other races requires further verification. Second, this study is a secondary analysis based on published data, so variables not included in the dataset cannot be adjusted. However, we used the E-value to evaluate the unmeasured confounding factors and found that our study was stable and reliable. Third, this study was based on a secondary analysis of published data and lacked some relevant information, such as recent surgical observations (perioperative and 30-day mortality), minor surgeries, and types of transplants. Fourth, PS methods balance the known confounding variables as much as possible. However, it could not ensure that all measured baseline characteristics were matched, and the influence of unknown variables was considered. However, we reduced the calliper width to 0.0005 to minimise the interference of some variables in the results. Fifth, differences in CKD stage and age remained between the CHF and non-CHF groups after PS matching. However, multivariate adjustment and subgroup analyses were performed. These analyses suggest that our results are robust. In addition, this observational study provides inferences about the association between CHF and 1-year mortality after surgery but cannot establish a causal relationship. Therefore, our findings need to be further validated by future prospective studies.

Conclusion

CHF is an independent risk factor for 1-year mortality after surgery in patients undergoing surgery. This study quantified the relationship between CHF and surgical patient outcomes by applying various statistical models and presented a range of HR (1.24–1.54). This provides a reference for optimising clinical decision-making, improving preoperative consultation, and promoting clinical communication.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board (Singhealth CIRB 2014/651/D). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author Contributions

YH, HH, and YL conceived the research, drafted the manuscript, and performed the statistical analysis. DL and LW revised the manuscript and designed this study. All authors have read and approved the final manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.858068/full#supplementary-material

Abbreviations

CHF, congestive heart failure; IPTW, inverse probability of treatment-weighted; SMR, standardised mortality ratio; HF, heart failure; PSM, propensity score matching; RCT, randomised controlled trials; CVA, previous cerebrovascular accidents; IHD, ischaemic heart disease; DMI, diabetes mellitus on insulin; ASA-PS, American Society of Anaesthesiologists Physical Status; RDW, red cell distribution width; HR, hazard ratios; SD, standardised differences; CAD, coronary artery disease.

Footnotes

1. ^ www.Datadryad.org
2. ^ http://www.R-project.org
3. ^ http://www.empowerstats.com

References

1. Yancy CW, Jessup M, Bozkurt B, Butler J, Casey DJ, Colvin MM, et al. 2017 ACC/AHA/HFSA focused update of the 2013 ACCF/AHA guideline for the management of heart failure: a report of the American college of cardiology/American heart association task force on clinical practice guidelines and the heart failure society of America. J Card Fail. (2017) 23:628–51. doi: 10.1016/j.cardfail.2017.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Roger VL, Go AS, Lloyd-Jones DM, Benjamin EJ, Berry JD, Borden WB, et al. Heart disease and stroke statistics–2012 update: a report from the American heart association. Circulation. (2012) 125:e2–220. doi: 10.1161/CIR.0b013e31823ac046

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ho KK, Anderson KM, Kannel WB, Grossman W, Levy D. Survival after the onset of congestive heart failure in Framingham heart study subjects. Circulation. (1993) 88:107–15. doi: 10.1161/01.cir.88.1.107

CrossRef Full Text | Google Scholar

4. Lee TH, Marcantonio ER, Mangione CM, Thomas EJ, Polanczyk CA, Cook EF, et al. Derivation and prospective validation of a simple index for prediction of cardiac risk of major noncardiac surgery. Circulation. (1999) 100:1043–9. doi: 10.1161/01.cir.100.10.1043

CrossRef Full Text | Google Scholar

5. Virani SS, Alonso A, Benjamin EJ, Bittencourt MS, Callaway CW, Carson AP, et al. Heart disease and stroke statistics-2020 update: a report from the American heart association. Circulation. (2020) 141:e139–596. doi: 10.1161/CIR.0000000000000757

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bohsali F, Klimpl D, Baumgartner R, Sieber F, Eid SM. Effect of heart failure with preserved ejection fraction on perioperative outcomes in patients undergoing hip fracture surgery. J Am Acad Orthop Surg. (2020) 28:e131–8. doi: 10.5435/JAAOS-D-18-00731

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Hernandez AF, Whellan DJ, Stroud S, Sun JL, O’Connor CM, Jollis JG. Outcomes in heart failure patients after major noncardiac surgery. J Am Coll Cardiol. (2004) 44:1446–53. doi: 10.1016/j.jacc.2004.06.059

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lerman BJ, Popat RA, Assimes TL, Heidenreich PA, Wren SM. Association of left ventricular ejection fraction and symptoms with mortality after elective noncardiac surgery among patients with heart failure. JAMA. (2019) 321:572–9. doi: 10.1001/jama.2019.0156

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Faxén UL, Hallqvist L, Benson L, Schrage B, Lund LH, Bell M. Heart failure in patients undergoing elective and emergency noncardiac surgery: still a poorly addressed risk factor. J Card Fail. (2020) 26:1034–42. doi: 10.1016/j.cardfail.2020.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Turrentine FE, Sohn MW, Jones RS. Congestive heart failure and noncardiac operations: risk of serious morbidity, readmission, reoperation, and mortality. J Am Coll Surg. (2016) 222:1220–9. doi: 10.1016/j.jamcollsurg.2016.02.025

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Smilowitz NR, Banco D, Katz SD, Beckman JA, Berger JS. Association between heart failure and perioperative outcomes in patients undergoing non-cardiac surgery. Eur Heart J Qual Care Clin Outcomes. (2021) 7:68–75. doi: 10.1093/ehjqcco/qcz066

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Robins JM, Greenland S. The role of model selection in causal inference from nonexperimental data. Am J Epidemiol. (1986) 123:392–402. doi: 10.1093/oxfordjournals.aje.a114254

PubMed Abstract | CrossRef Full Text | Google Scholar

13. D’Agostino RJ. Propensity score methods for bias reduction in the comparison of a treatment to a non-randomized control group. Stat Med. (1998) 17:2265–81. doi: 10.1002/(sici)1097-0258(19981015)17:193.0.co;2-b

CrossRef Full Text | Google Scholar

14. Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivariate Behav Res. (2011) 46:399–424. doi: 10.1080/00273171.2011.568786

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kurth T, Walker AM, Glynn RJ, Chan KA, Gaziano JM, Berger K, et al. Results of multivariable logistic regression, propensity matching, propensity adjustment, and propensity-based weighting under conditions of nonuniform effect. Am J Epidemiol. (2006) 163:262–70. doi: 10.1093/aje/kwj047

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Rubin DB, Thomas N. Matching using estimated propensity scores: relating theory to practice. Biometrics. (1996) 52:249–64.

PubMed Abstract | Google Scholar

17. Cheng D, Chakrabortty A, Ananthakrishnan AN, Cai T. Estimating average treatment effects with a double-index propensity score. Biometrics. (2020) 76:767–77. doi: 10.1111/biom.13195

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Sato T, Matsuyama Y. Marginal structural models as a tool for standardization. Epidemiology. (2003) 14:680–6. doi: 10.1097/01.EDE.0000081989.82616.7d

CrossRef Full Text | Google Scholar

19. Sim YE, Wee HE, Ang AL, Ranjakunalan N, Ong BC, Abdullah HR. Prevalence of preoperative anemia, abnormal mean corpuscular volume and red cell distribution width among surgical patients in Singapore, and their influence on one year mortality. PLoS One. (2017) 12:e182543. doi: 10.1371/journal.pone.0182543

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Glance LG, Lustik SJ, Hannan EL, Osler TM, Mukamel DB, Qian F, et al. The surgical mortality probability model: derivation and validation of a simple risk prediction rule for noncardiac surgery. Ann Surg. (2012) 255:696–702. doi: 10.1097/SLA.0b013e31824b45af

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Mo Z, Hu H, Du X, Huang Q, Chen P, Lai L, et al. Association of evaluated glomerular filtration rate and incident diabetes mellitus: a secondary retrospective analysis based on a Chinese cohort study. Front Med (Lausanne). (2021) 8:724582. doi: 10.3389/fmed.2021.724582

PubMed Abstract | CrossRef Full Text | Google Scholar

22. World Health Organization Technical Report Series. Nutritional anaemias. Report of a WHO scientific group. World Health Organ Tech Rep Ser. (1968) 405:5–37.

PubMed Abstract | Google Scholar

23. Groenwold RH, White IR, Donders AR, Carpenter JR, Altman DG, Moons KG. Missing covariate data in clinical research: when and when not to use the missing-indicator method for analysis. CMAJ. (2012) 184:1265–9. doi: 10.1503/cmaj.110977

PubMed Abstract | CrossRef Full Text | Google Scholar

24. White IR, Royston P, Wood AM. Multiple imputation using chained equations: issues and guidance for practice. Stat Med. (2011) 30:377–99. doi: 10.1002/sim.4067

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ahmed A, Husain A, Love TE, Gambassi G, Dell’Italia LJ, Francis GS, et al. Heart failure, chronic diuretic use, and increase in mortality and hospitalization: an observational study using propensity score methods. Eur Heart J. (2006) 27:1431–9. doi: 10.1093/eurheartj/ehi890

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Normand ST, Landrum MB, Guadagnoli E, Ayanian JZ, Ryan TJ, Cleary PD, et al. Validating recommendations for coronary angiography following acute myocardial infarction in the elderly: a matched analysis using propensity scores. J Clin Epidemiol. (2001) 54:387–98. doi: 10.1016/s0895-4356(00)00321-8

CrossRef Full Text | Google Scholar

27. Wu Y, Hu H, Cai J, Chen R, Zuo X, Cheng H, et al. Association of hypertension and incident diabetes in Chinese adults: a retrospective cohort study using propensity-score matching. BMC Endocr Disord. (2021) 21:87. doi: 10.1186/s12902-021-00747-0

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Brown JB, Gestring ML, Stassen NA, Forsythe RM, Billiar TR, Peitzman AB, et al. Geographic variation in outcome benefits of helicopter transport for trauma in the United States: a retrospective cohort study. Ann Surg. (2016) 263:406–12. doi: 10.1097/SLA.0000000000001047

PubMed Abstract | CrossRef Full Text | Google Scholar

29. McCaffrey DF, Griffin BA, Almirall D, Slaughter ME, Ramchand R, Burgette LF. A tutorial on propensity score estimation for multiple treatments using generalized boosted models. Stat Med. (2013) 32:3388–414. doi: 10.1002/sim.5753

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Koch B, Vock DM, Wolfson J. Covariate selection with group lasso and doubly robust estimation of causal effects. Biometrics. (2018) 74:8–17. doi: 10.1111/biom.12736

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Geleris J, Sun Y, Platt J, Zucker J, Baldwin M, Hripcsak G, et al. Observational study of hydroxychloroquine in hospitalized patients with covid-19. N Engl J Med. (2020) 382:2411–8. doi: 10.1056/NEJMoa2012410

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Robins JM, Hernán MA, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology. (2000) 11:550–60. doi: 10.1097/00001648-200009000-00011

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Qin H, Chen Z, Zhang Y, Wang L, Ouyang P, Cheng L, et al. Triglyceride to high-density lipoprotein cholesterol ratio is associated with incident diabetes in men: a retrospective study of Chinese individuals. J Diabetes Investig. (2020) 11:192–8. doi: 10.1111/jdi.13087

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Raynes-Greenow CH, Hadfield RM, Cistulli PA, Bowen J, Allen H, Roberts CL. Sleep apnea in early childhood associated with preterm birth but not small for gestational age: a population-based record linkage study. Sleep. (2012) 35:1475–80. doi: 10.5665/sleep.2192

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Haneuse S, VanderWeele TJ, Arterburn D. Using the e-value to assess the potential effect of unmeasured confounding in observational studies. JAMA. (2019) 321:602–3. doi: 10.1001/jama.2018.21554

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Mathew A, Devereaux PJ, O’Hare A, Tonelli M, Thiessen-Philbrook H, Nevis IF, et al. Chronic kidney disease and postoperative mortality: a systematic review and meta-analysis. Kidney Int. (2008) 73:1069–81. doi: 10.1038/ki.2008.29

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhang M, Lin S, Wang MF, Huang JF, Liu SY, Wu SM, et al. Association between NAFLD and risk of prevalent chronic kidney disease: why there is a difference between east and west? BMC Gastroenterol. (2020) 20:139. doi: 10.1186/s12876-020-01278-z

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Harris BN, Pipkorn P, Nguyen K, Jackson RS, Rao S, Moore MG, et al. Association of adjuvant radiation therapy with survival in patients with advanced cutaneous squamous cell carcinoma of the head and neck. JAMA Otolaryngol Head Neck Surg. (2019) 145:153–8. doi: 10.1001/jamaoto.2018.3650

PubMed Abstract | CrossRef Full Text | Google Scholar

39. von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP. The strengthening the reporting of observational studies in epidemiology (strobe) statement: guidelines for reporting observational studies. Int J Surg. (2014) 12:1495–9. doi: 10.1016/j.ijsu.2014.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Siu JM, McCarty JC, Gadkaree S, Caterson EJ, Randolph G, Witterick IJ, et al. Association of vessel-sealant devices vs conventional hemostasis with postoperative neck hematoma after thyroid operations. JAMA Surg. (2019) 154:e193146. doi: 10.1001/jamasurg.2019.3146

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Lerman BJ, Popat RA, Assimes TL, Heidenreich PA, Wren SM. Association between heart failure and postoperative mortality among patients undergoing ambulatory noncardiac surgery. JAMA Surg. (2019) 154:907–14. doi: 10.1001/jamasurg.2019.2110

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Khuri SF, Henderson WG, DePalma RG, Mosca C, Healey NA, Kumbhani DJ. Determinants of long-term survival after major surgery and the adverse effect of postoperative complications. Ann Surg. (2005) 242:326–41. doi: 10.1097/01.sla.0000179621.33268.83

CrossRef Full Text | Google Scholar