The Hyperbilirubinemia and Potential Predictors Influence on Long-Term Outcomes in Sepsis: A Population-Based Propensity Score-Matched Study

Objective: Although hyperbilirubinemia has been associated with mortality in patients who are critically ill, yet no clinical studies dissect the effect of dynamic change of hyperbilirubinemia on long-term septic prognosis. The study aims to investigate the specific stages of hyperbilirubinemia and potential risk factors on long-term outcomes in patients with sepsis. Methods: In this retrospective observational cohort study, patients with sepsis, without previous chronic liver diseases, were identified from the Medical Information Mart for the Intensive Care III MIMIC-III database. We used propensity scores (PS) to adjust the baseline differences in septic patients with hyperbilirubinemia or not. The multivariate Cox was employed to investigate the predictors that influence a clinical outcome in sepsis. Results: Of 2,784 patients with sepsis, hyperbilirubinemia occurred in 544 patients (19.5%). After PS matching, a survival curve demonstrated that patients with sepsis with the new onset of total bilirubin (TBIL) levels more than or equal to 5 mg/dl survived at significantly lower rates than those with TBIL levels <5 mg/dl. Multivariate Cox hazard analysis showed that patients with TBIL at more than or equal to 5 mg/dl during sepsis exhibit 1.608 times (95% CI: 1.228–2.106) higher risk of 1-year mortality than those with TBIL levels <5 mg/dl. Also, age above 65 years old, preexisting malignancy, a respiratory rate above 30 beats/min at admission, serum parameters levels within 24-h admission, containing international normalized ratio (INR) above 1.5, platelet <50*10∧9/L, lactate above 4 mmol/L, and bicarbonate <22 or above 29 mmol/L are the independent risk factors for long-term mortality of patients with sepsis. Conclusions: After PS matching, serum TBIL levels at more than or equal to 5 mg/dl during hospitality are associated with increased long-term mortality for patients with sepsis. This study may provide clinicians with some cutoff values for early intervention, which may improve the prognosis of patients with sepsis.


INTRODUCTION
Sepsis is defined as a life-threatening acute organ dysfunction secondary to infection. A high incidence rate and high mortality of sepsis make it one of the leading causes of death as a global health priority (1,2). Notably, sepsis-induced organ dysfunction is an important predictor for poor prognosis (3)(4)(5). The liver plays a central role in homeostasis, immune surveillance, inflammation, and bacterial clearance (6,7). A large body of evidence has suggested that the liver is the main target of sepsis and decompensation of liver function can trigger overwhelming inflammation, immune response, and organ damage in sepsis (8,9). However, despite extensively studied lung, kidney, and heart injury in the course of sepsis, the question of whether dysfunction of the liver is associated with mortality or a poor outcome in sepsis remains unresolved. Hepatic dysfunction and hyperbilirubinemia commonly occur in patients who are critically ill with an incidence rate of 40% and up to 20% in patients with bacterial infection (10).
Hyperbilirubinemia may result from bacterial products or as a consequence of the response of the host to infection. The etiology of hyperbilirubinemia in patients who are critically ill is multifactorial, probably cholestasis or sclerosing cholangitis caused by circulating endotoxins, inflammation, hypoxia hepatitis, lower liver perfusion and ischemia, genetic and metabolic variations, and so on (11)(12)(13)(14)(15)(16). Hyperbilirubinemia has been shown to represent an important marker of mortality and poor outcomes for patients who are critically ill (17,18). Patients with a preexisting liver deficit, like cirrhosis, have a worse outcome of sepsis than ones without liver dysfunction due to impaired immunity (19). The underlying mechanism may relate to the reversal of bilirubin transport from intrahepatic toward the circulation, which is beneficial for relieving the highenergy burden for hepatocytes and serves as a metabolic and inflammatory stress response as well (20). Hence, elevated serum bilirubin might indicate the impairment of energy consumption due to liver injury, and total serum bilirubin level has been widely recognized as a powerful maker for assessing hepatic function compared with other serum activities on laboratory tests.
Thus far, several studies have demonstrated the potential effect of liver dysfunction on short-term mortality of patients with sepsis by using different definitions of hepatic dysfunction (21)(22)(23). Whereas, it is still unclear of the long-term impact of dynamic change of hyperbilirubinemia in patients with sepsis. The appendant result derived from PROWESS-SHOCK (Prospective Recombinant Human Activated Protein C Worldwide Evaluation in Severe Sepsis and Septic Shock) trial showed that liver dysfunction is associated with 180-day Abbreviations: MIMIC, medical information mart for intensive care; STROBE, STrengthening the Reporting of OBservational studies in Epidemiology; ICD-9, International Classification of Diseases, 9th revision; IRB, institutional review board; SOFA, sequential organ failure assessment; TBIL, total bilirubin; CHF, congestive heart failure; AFIB, atrial fibrillation; CAD, coronary artery disease; COPD, chronic obstructive pulmonary disease; ALT, alanine aminotransaminase; AST, aspartate aminotransferase; CI, confidence interval; HR, hazard ratio; OR, odds ratio; pCO 2 , partial pressure of CO 2 ; INR, international normalized ratio; BAIBF, bile acid-independent bile flow; USA, United States. mortality in patients with septic shock by using criteria of serum bilirubin ≥20 mmol/L (24). However, as hyperbilirubinemia is multifactorial, there is no large study dissecting the effect of liver dysfunction on long-term mortality of more than 180 days in patients with general sepsis. Therefore, we performed a large cohort study to analyze the dynamic change of hyperbilirubinemia and its influence on outcomes in patients with sepsis by using an openly available US-based critical care database named Medical Information Mart for Intensive Care (MIMIC)-III v 1.4, which includes 52,963 ICU admissions. Our results revealed a significant association between total bilirubin (TBIL) levels with long-term mortality in patients with sepsis.

Study Population
We conducted the cohort study according to the STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) statement by using the MIMIC-III database, a large, integrated, de-identified, open-free, comprehensive clinical dataset, comprised of all the patients admitted to the ICUs Beth Israel Deaconess Medical Center in Boston, MA, from June 2001 to October 2012. For all of the data that are deidentified, patient consent or ethics approval will not be needed. Demographics records, laboratory results, radiology examinations, diagnosis, clinical treatment parameters, and dates of death were also concluded. The diagnosed diseases by the physician were according to the International Classification of Diseases, 9th revision (ICD-9) on patient discharge. Since the study was an analysis of a third-party anonymized publicly available database with preexisting institutional review board (IRB) approval, approval from our institution was exempted.
Sepsis was defined according to Sepsis-3 criteria: the suspected infection and Sequential Organ Failure Assessment (SOFA) score was of 2 points or more (25,26). Hyperbilirubinemia was diagnosed in patients with the new onset of serum TBIL at more than or equal to 2 mg/dl during the hospitalization. The included criteria were as follows: age ≥18 years old; patients without previous chronic liver diseases according to the recorded ICD-9 codes, including liver cirrhosis; at ICU admission more than 24 h; missing data <50%. The excluded criteria were preexisting bilirubin at more than or equal to 2 mg/dl before admission, and previous chronic liver diseases like chronic hepatitis, acute-onchronic liver failure, cirrhosis, liver cancer, hepatobiliary ductrelated tumors, or acute liver conditions like drug or toxininduced hepatitis. After being included, all the patients with sepsis were divided into a hyperbilirubinemia group (serum TBIL during hospitalization at more than or equal to 2 mg/dl) and a non-hyperbilirubinemia group (serum TBIL <2 mg/dl).

Measures and Variable Definition
For the patients in the study, we retrieved demographic and admission information from the database during the first 24 h of ICU admission, including age, gender, ethnicity (White, Hispanic, Black, or other), weight, time of admission or discharge, the severity of illness parameters: the SOFA score (clarified into four different strata: 2-4, 5-9, more than or equal to 10 scores) and the Elixhauser comorbidity score, and vital signs: heart rate (clarified into two different strata: <100 beats/min, more than or equal to 100 beats/min), and respiratory rate (clarified into two different strata: <30 beats/min, more than or equal to 30 beats/min). In addition, we routinely collected laboratory parameters within the first 24 h of ICU admission, including maximum levels of white blood cell count (WBC); levels of hemoglobin (g/dl); hematocrit; platelet levels (stratified into five different scales: more than or equal to 250 * 10 ∧ 9/L, at 150-249 * 10 ∧ 9/L, at 100-149 * 10 ∧ 9/L, at 50-99 * 10 ∧ 9/L, and <50 * 10 ∧ 9/L); serum potassium/sodium/chloride levels; serum bicarbonate levels (stratified into three different scales: <22 mmol/L, at 22-29 mmol/L, and more than 29 mmol/L); serum blood urea nitrogen (BUN)/creatinine; serum lactate levels (stratified into four different scales: at 0-2 mmol/L, at 2.1-4 mmol/L, at 4.1-10 mmol/L, and more than 10 mmol/L); international normalized ratio (INR, stratified into two different scales: ≤1.5, more than 1.5); arterial blood gas with PH, partial pressure of oxygen (pO2), partial pressure of carbon dioxide (pCO2); alanine transaminase (ALT)/aspartate transaminase (AST). We also collected parameters in the course of sepsis: use of mechanical ventilation, use of vasopressor agents, and use of sedative drugs. By collecting in 24, 48, and 72 h, 7 days, and the day of discharge after ICU admission, serum levels of TBIL were classified into four different scales: at 0-1.9, 2-4.9, 5-10 mg/dl, and more than 10 mg/dl.
We also included preexisting medical comorbidities according to the recorded ICD-9 codes, including congestive heart failure (CHF), renal disease, atrial fibrillation (AFIB), coronary artery disease (CAD), chronic obstructive pulmonary disease (COPD), stroke, and malignant tumor. The Kidney Disease: Improving Global Outcomes (KDIGO) clinical practice guidelines were used to define acute kidney injury (AKI). We conducted the followup at day 30 (30 days), 90 days, 180 days, and 1 year from the database.

Statistical Analysis
Statistical analysis was performed using SPSS 23 (SPSS, Inc., Chicago, IL). The parameters with missing data of more than 50% were excluded from our study. The mean-value imputation algorithm was selected to substitute missing values. Baseline characteristics and clinical parameters after ICU admission between the hyperbilirubinemia group and the control group were compared. All continuous variables were expressed as means (SD) or medians (interquartile range, IQR) by using either Student t-test or Mann-Whitney U test as appropriate. Categorical variables were compared by the chi-square test or Fisher's exact test. The survival curves for patients with sepsis were plotted using the Kaplan-Meier method, and differences between the curves were assessed using the log-rank test. To identify the association between the dynamic change and the levels of hyperbilirubinemia and a long-term outcome in sepsis, Cox regression analysis was used. Univariate Cox proportional hazards regression of clinical parameters was performed to identify potential predictors. Hazard ratios (HRs) were calculated with 95% CIs as an estimate of the risk associated with a particular variable. To determine independent predictors of the composite end points, variables in univariate Cox analysis with p < 0.1 were entered into multivariate Cox proportional hazards regression, and the predictors were performed with a likelihood ratio-forward selection.
To account for selection bias and potential confounding factors between groups in comparison of an outcome, we used the propensity score (PS) matching (1:1) to balances covariates for those who had hyperbilirubinemia and those who had not (532 pairs). A multivariable logistic regression model with confounding baseline characteristics was used to calculate the PS for each patient as the predicted probability of the hyperbilirubinemia group. The following variables were adjusted ( Table 1): age; sex; race/ethnicity; preexisting medical conditions: CHF, AFIB, CAD, COPD, stroke, malignant tumor, and chronic renal disease; Elixhauser comorbidity score; respiratory rate; and biochemical parameters. One-to-one nearest-neighbor matching without replacement with a caliper width of.1 was conducted. We evaluated the balance test after matching, with no significant difference with chi-square (X 2 = 16.468, P = 0.870), which demonstrated a good balance following PS matching between the hyperbilirubinemia and non-hyperbilirubinemia groups.
After matching, survival curves were computed and plotted using the Kaplan-Meier method. Univariate and multivariate Cox proportional hazards regression models were constructed to explore the independent risk factors of influence on long-term outcomes in sepsis. A P < 0.05 (two-sided) was considered significant.

Patient Characteristics Before Matching
There was also a significantly higher severity of illness in the hyperbilirubinemia group than the non-hyperbilirubinemia group, with higher rates of the SOFA score at more than 4 (83.9 vs. 56.1%; p < 0.001), higher levels of Elixhauser comorbidity index [5. (10) vs. 4. (9)], and higher frequency of the respiratory rate at more than or equal to 30 beats/min (45.6 vs. 38.4%; P = 0.002; Table 1). At 24 h after admission to ICU, the   Table 1).
The patients from the non-hyperbilirubinemia group all had TBIL levels below 1.9 mg/dl. By contrast, 62.1% of patients from the hyperbilirubinemia group had serum TBIL levels at 2-4.9 mg/dl, 25.9% had TBIL levels at 5-10 mg/dl, and 11.9% had TBIL levels at more than 10 mg/dl ( Table 2). In addition, the hyperbilirubinemia group had higher rates of vasopressor usage when compared with the non-hyperbilirubinemia group. The overall in-hospital mortality of patients with sepsis was 23.5%. The patients in the hyperbilirubinemia group had significantly higher rates of in-hospital mortality (16.2 vs. 10.7%, P = 0.002), 30-day mortality (21.7 vs. 16.3%, P = 0.011), 90-day mortality (25.4 vs. 19.5%, P = 0.010), 180-day mortality (27.2 vs. 22.1%, P = 0.039), and 1-year mortality after discharge (29.8 vs. 24.5%, P = 0.037) than in the non-hyperbilirubinemia group, respectively ( Table 2). Kaplan-Meier's analysis also showed that the 1-year survival rate was significantly lower in the hyperbilirubinemia group than in the non-hyperbilirubinemia group (P = 0.006) before matching (Figure 2A). Moreover, it is of note that TBIL levels more than or equal to 5 mg/dl at the hospital significantly increases the risk of mortality in sepsis (p < 0.001, Figure 2B).

Univariate and Multivariate Cox Hazard Analysis of Risk Factors for Mortality in Sepsis Before Matching
Having identified that several laboratory serum activities, such as serum platelet levels, bicarbonate levels, lactate levels, and INR levels that are significantly associated with hyperbilirubinemia in ICU, we tested whether these serum activities, as well as other clinical outcomes, could be predictors of sepsis prognosis. We first deployed the univariate Cox proportional hazard model to analyze the predictors of mortality in sepsis. We employed a variety of categorical variables with a reference variable: serum platelet levels with platelet levels above or equal to 250 * 10 ∧ 9/L being a reference variable; serum lactate levels with lactate levels ≤2 mmol/L being a reference variable; serum bicarbonate levels with bicarbonate levels at 22-29 mmol/L being a reference variable; serum TBIL levels with TBIL levels <2 mg/dl being a reference variable; the AKI KDIGO stage without AKI as a reference variable. Univariate analysis indicated that age at more than 65 years, weight, race in Black, Hispanic, or other, preexisting medical conditions (Malignancy, Stroke, AFIB, and renal disease), Elixhauser comorbidity index, respiratory rate at more than 30 m beats/min, serum INR levels at more than 1.5, serum platelet levels at <50 * 10 ∧ 9/L, serum lactate levels at more than 2 mmol/L, serum bicarbonate levels at more than 29 mmol/L, or <22 mmol/L, AKI at two or three stages, mechanical ventilation usage, vasoactive drug usage, and serum TBIL at more than or equal to 5 mg/dl were the significant risk factors in mortality in sepsis.
To test the prognostic predictors of mortality in sepsis, we performed multivariate analysis using the Cox proportional hazard model for all variables identified as significant by univariate analysis. Following control of confounders and Likelihood ratio (LR) forward elimination, results indicated (  year after admission of the last patients before matching, with a significant difference between liver injury and non-liver injury groups (Log-rank P = 0.006). (B) shows the survival curves for different stages of serum total bilirubin (TBIL) levels before propensity matching. Kaplan-Meier analysis showed survival time did not differ between the hyperbilirubinemia group and the normal group. Log-rank p < 0.001. (C) shows the survival curves between liver injury and non-liver injury groups after matching (log-rank P = 0.084); (D) shows the survival curves for different stages of serum TBIL levels after propensity matching (Log-rank p < 0.001).

PS Analysis
From the above results, we found that the baseline clinical information significantly influenced the hyperbilirubinemia group and long-term clinical outcomes. Thus, we applied propensity-score matching to minimize confounding biases. One-to-one propensity-score matching yielded a cohort of 532 patients in the hyperbilirubinemia group and 532 in the control group. Baseline characteristics for patients with and without hyperbilirubinemia were well-balanced after matching, as shown in Table 4. Standardized biases for all variables were.05 or less. The SOFA score at admission was excluded as a matching variable for one of its evaluated parameters, including TBIL. In the hyperbilirubinemia group after matching, the rates of the SOFA score at admission above 4 was significantly higher in the hyperbilirubinemia group when compared with the non-hyperbilirubinemia group (83.4 vs. 59.8%; p < 0.001).
The matched results also showed that the patients with sepsis with the new onset of hyperbilirubinemia were associated with significantly increased risks of mortality at the hospital (16. vs. 11.8%, P = 0.036), 30 days (21.4 vs. 16.7%, P = 0.038), 90 days (25 vs. 19.5%, P = 0.027), and 1 year (29.5 vs. 24.4%, P = 0.047), but not statistically increased in 90 days, 180 days, and 1-year mortality in the overall population ( Table 5). The duration of mechanical ventilation, vasoactive drug or sedative drug usage, and AKI stage did not significantly differ between the two groups. After matching, Kaplan-Meier's analysis showed that the patients with new onset hyperbilirubinemia had a lower long-term survival rate compared with the non-hyperbilirubinemia group but did not reach statistical significance (P = 0.084; Figure 2C). However, TBIL levels more than or equal to 5 mg/dl at the hospital also   significantly increases the risk of mortality in sepsis (p < 0.001, Figure 2D).
To further verify the levels of TBIL at the hospital on clinical outcomes, we then classified the patients into TBIL levels at the <5 mg/dl group and TBIL levels at the ≥5 mg/dl group. The results demonstrated that patients with sepsis with serum TBIL levels at ≥5 mg/dl during the hospital stage could bring significant poor outcomes on 30-day (28.1 vs. 17%), 180-day (34.7 vs. 19.4%), and 1-year mortality (37.7 vs. 21.7%), but with noninfluence on hospital mortality (18.6 vs. 12.8%) when compared with those with TBIL levels at the <5 mg/dl group (Supplementary Table 1).

DISCUSSION
In this large cohort study, we find the incidence rate of the new onset of hyperbilirubinemia during sepsis is 19.5%, which is considerably lower than the rate of hyperbilirubinemia in other patients who are critically ill (14,27). The difference is mainly due to the hyperbilirubinemia in our study that is defined as the new onset excluding past liver diseases, while previous studies included patients with preexisting chronic liver conditions (22,28). Our data show that, after PS matching, the new onset of hyperbilirubinemia during sepsis cannot significantly increase the risk of long-term mortality, but the patients with TBIL levels are more than or equal to 5 mg/dl do. Further multivariate Cox hazard analysis shows that age above 65 years old, other race, preexisting malignancy, respiratory rate above 30 beats/min at admission, serum parameters levels within 24-h admission containing INR above 1.5, platelet <50 * 10 ∧ 9/L, lactate above 4 mmol/L, bicarbonate <22 or above 29 mmol/L, serum TBIL during sepsis at more than or equal to 5 mg/dl are the independent risk factors in 1-year mortality of patients with sepsis. Furthermore, the low level of platelet, the INR level more than 1.5, the respiratory rate more than or equal to 30 (beats/min), the higher creatinine level, the higher arterial pCO2, and the lower bicarbonate level in the hyperbilirubinemia group compared with the non-hyperbilirubinemia group indicate that there is a close correlation between hyperbilirubinemia and the deterioration of hematologic, coagulation, respiratory, urinary, and acid-base balance systems function in sepsis. In stepwise multivariate Cox hazard analysis, Elixhauser comorbidity index, the respiratory rate above 30 beats/min at admission, serum INR above 1.5, serum platelet <50 * 10 ∧ 9/L, serum lactate more than 4 mmol/L, serum bicarbonate <22 mmol/L, or above 29 mmol/L at 24 h of admission, and mechanical ventilation as independent risk factors in 1-year mortality of patients with sepsis hint that hematologic, coagulation, respiratory, and acid-base systems that are imbalance at the early stage of sepsis predict poor prognosis of long-term mortality of sepsis. In addition, when serum TBIL at more than or equal to 5 mg/dl as a predictive factor in 1-year mortality of patients with sepsis consolidates our speculation: the new onset of hyperbilirubinemia reach at some extent aggravates sepsis prognosis.
Intriguingly, we observe that the hyperbilirubinemia group has a significantly higher rate of preexisting AFIB and malignancy but a lower rate of past COPD, stroke, and renal disease than the non-hyperbilirubinemia group. The mechanism of the potential relationship between the new onset of hyperbilirubinemia and past medical history in sepsis is not clear and needs further study. Multivariate Cox analysis results suggest that preexisting diseases, such as malignancy, affect the long-term survival of patients with sepsis, which is consistent with recent findings (29). Given that great bias was frequently generated when analyzing the sole effect of hyperbilirubinemia on sepsis, we used PS matching analysis to balance baseline characteristics to minimize confounding bias. After matching, baseline variants are well adjusted to parallel, thus reducing confounding bias to the maximum extent. Importantly, the new onset of hyperbilirubinemia when TBIL levels achieve at more than or equal to 5 mg/dl is associated with significantly increased long-term mortality of patients with sepsis after matching. Compared with the previous study, which did not use PS matching to strictly control possible confounders (24), our study provides solid evidence and comprehensive dissection of the relationship between the dynamic change of hyperbilirubinemia and sepsis prognosis, and, at the same time, we seek a cutoff value of TBIL levels associated with a poor outcome in sepsis, which may guide for early intervention for ICU physicians. Our results are consistent with previous basic research, which shows that the liver plays an important role in endotoxin-induced acute lung injury (30). Moreover, Zhang et al. (31) have proposed a Cox regression model with time-varying covariates, which is useful for evaluating the dynamic change of hyperbilirubinemia on clinical outcomes and can be used in further analysis. The liver is the key detoxification organ and the critical site to clear invasive pathogens and alleviate inflammation reaction by exerting innate immune system function (6,7). Meanwhile, liver dysfunction and associated hyperbilirubinemia lead to inflammation and immune response out of control and cascading organ damage in sepsis (8,9). Furthermore, we explore the independent risk factors affecting 1-year mortality of patients with sepsis after PS matching, and we find that old age, preexisting malignancy, respiratory rate above 30 beats/min at admission, serum parameters within 24h admission with INR above 1.5, serum platelet <50 * 10 ∧ 9/L, serum lactate more than 4 mmol/L, and bicarbonate <22 or above 29 mmol/L are the independent risk factors in 1year mortality of patients with sepsis. These parameters may provide early warning of the prognosis of patients with sepsis. Previous basic pieces of research show that sepsis-induced acute kidney and myocardial injury are age-dependent (32,33), and a recent clinical study has shown that the odds for mortality of patients with sepsis in ICU increase with age (34), which are consistent with our results. In qSOFA definition, expert consensus shows the respiratory rate of 22/min or greater, which is the predictor for a poor outcome of patients with sepsis (25); our results further demonstrate the respiratory rate of 30/min or greater at admission is the independent risk factor in long-term mortality. Abundant studies have focused on the correlation between the serum lactate level and shortterm mortality of patients with sepsis; the lactate level ≥2 mmol/L around is demonstrated having the predictive value (35,36). However, the cutoff of the lactate level for predicting long-term mortality of sepsis remains unclear. Our results show lactate above 4 mmol/L within 24-h admission is the independent predictor for 1-year mortality of sepsis, and patients with lactate above 10 mmol/L will bring a poor prognosis in sepsis. Of note, recent studies have shown that platelets play a vital role in immunological surveillance against pathogens invaders and contribute to innate immune system function (37,38). Wong et al. demonstrated that platelets collaborate with macrophages to fight against certain blood-borne infections. In addition, the absence of the platelet resulted in the platelet being unable to localize to the sites of infection, leading to rapid death of the Kupffer cells and endothelium, followed by more leakage of plasma out of blood vessels, and even host mortality (39). This may explain why patients with sepsis with lower-level platelets have poor outcomes in our study. Besides, we find that bipolar serum bicarbonate levels are another risk factor in the long-term mortality of patients with sepsis. Previous experimental studies indicate that endotoxemia significantly decreases bile acid-independent bile flow (BAIBF) and associated biliary HCO − 3 output (40), thus reducing the serum bicarbonate level. As there are no studies that have reported the correlation between the bicarbonate level and sepsis mortality, our study first shows that bicarbonate <22 or above 29 mmol/L is an independent predictor for 1-year mortality of sepsis.
To our knowledge, this is the first and largest study so far to dissect the correlation between the dynamic change of hyperbilirubinemia and outcomes of sepsis. However, there are several limitations to our study. First, we used the database from a single academic medical center in the USA; therefore, some of the cases from almost 10 years ago, as diagnosis or treatment strategies at that time, would be inconsistent with current guidelines, which brings great bias. There is residual confounding by variables not collected into the MIMIC-III database. However, we include all the patients with sepsis according to the uniform standard of Sepsis 3.0 and apply PSs analysis and match the baseline characteristics of the patients to eliminate confounding factors and decrease bias to the maximum. Second, there are a few missing data that bring bias, yet we delete the data whose missing percentage is larger than 50% to decrease bias. The meanvalue imputation algorithm is selected to substitute missing values. Third, the single-centered design restricts generalizability to apply our conclusion to other regions, while we use a large sample size and PS matching analysis to guarantee the quality of our study.

CONCLUSION
In conclusion, patients with serum TBIL at more than or equal to 5mg/dl during sepsis decrease survival rates after PS matching. In addition, we conclude that age above 65 years old, preexisting malignancy, respiratory rate above 30 beats/min at admission, serum parameter levels within 24-h admission, containing INR above 1.5, platelet <50 * 10 ∧ 9/L, lactate above 4 mmol/L, and bicarbonate <22 or above 29 mmol/L are independent risk predictors for long-term mortality of sepsis. Our study provides solid evidence and will rekindle the awareness of the risk factors leading to poor prognosis in sepsis.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Files, further inquiries can be directed to the corresponding authors.

ETHICS STATEMENT
The Ethics approval statement was approved by the Institutional Review Boards of Beth Israel Deaconess Medical Center (Boston, MA), and all the data were under a deidentification process to protect individual privacy and waived the need for a patient consent statement. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
MP, DQ, LZ, and ZH contributed to the conception and design of the research. MP and FD contributed to the acquisition and analysis of the data. MP, ZH, and LZ contributed to the interpretation of the results. MP and DQ drafted the manuscript. DQ, ZH, and LZ revised the manuscript. All authors have agreed to be fully accountable for ensuring the accuracy of the work, and approved the final manuscript.

FUNDING
This work was supported by Changsha Municipal Natural Science Foundation (grant number: kq2007078).