Liver injury on admission linked to worse outcomes in COVID-19: an analysis of 14,138 patients
Introduction
Spread of the SARS-CoV-2 virus occurred rapidly across the globe, resulting in one of the most challenging pandemics of the modern era (1). Despite the acquisition of a significant amount of knowledge about this peculiar virus, we continue to uncover new information. It targets not only the respiratory system but other organs as well including the liver (2). Given the potential for involvement of multiple systems, it is important for providers to recognize extra-pulmonary manifestations. The impact of coronavirus disease 2019 (COVID-19) on the medical and non-medical realms of humanity has been profound, and the efforts to overcome the pandemic continue to be underway with hope emerging with the arrival of vaccines (3). An abundance of information regarding the impact of the virus on the respiratory system has become available. The implications for other organs such as the liver are less well studied in comparison to the lungs though a plethora of studies has emerged over the past year with attempts to further solidify the effects on organs such as the liver. The impact on the liver can be variable and the degree of liver injury can vary, though patients may present with only mild abnormalities in liver chemistries. In this nation-wide multicenter study, we investigated liver chemistry abnormalities present on admission and the impact of those abnormalities on mortality, intubation, and hospital length of stay (LOS). Baseline characteristics of patients with abnormalities in liver chemistry are also evaluated with certain factors linked to worse outcomes. We present the following article in accordance with the STROBE reporting checklist for cohort studies (available at https://tgh.amegroups.com/article/view/10.21037/tgh-21-94/rc).
Methods
Objectives
Our objectives were to (I) assess baseline clinical and demographic characteristics in patients with COVID-19 infection who had and did not have abnormalities in liver chemistries [alkaline phosphatase (ALP), aspartate aminotransferase (AST), alanine aminotransferase (ALT), and total bilirubin (Tbili)] and (II) evaluate associations between abnormalities in liver chemistries and the primary outcomes of in-hospital death, intubation, and hospital LOS.
We conducted a retrospective cohort study using the Cerner Real-World DataTM (CRWD) de-identified COVID-19 patient cohort (4). CRWD was extracted from the electronic medical record (EMR) of hospitals with which Cerner had a data use agreement. All authors had access to the study data and reviewed and approved the final manuscript.
Ethical approval statement
The University of Missouri Institutional Review Board (Federalwide Assurance ID: FWA00002876) deemed this study of de-identified patient data to not constitute human subjects research; therefore, direct informed consent from individual patients was not required (IRB #2034442). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
Data source and participants
The December 2020 version (the 2020Q3 CRWD COVID database release) of Cerner’s COVID-19 patient cohort contains EMR data from 87 U.S.-based health systems from December 2019 through part of September 2020. Information about the more than 100 demographic, clinical, and non-clinical variables included in the dataset has been previously published (5).
We limited our analysis to the first encounter at which each patient had a lab test with a positive result consistent with COVID-19 infection. Logical Observation Identifiers Names and Codes (LOINC) terms used to define SARS-CoV-2 infection are in Table S1. Within this patient subset, we identified patients who had a documented laboratory result for ALP, AST, ALT, and Tbili during the 24-hour before or after the documented start time of the COVID encounter. LOINC codes used to identify liver chemistry test results are in Table S2. Patients who did not have documented laboratory results for all four liver chemistries within 24 hours of the start time of the COVID encounter were excluded. Figure 1 contains a flow diagram depicting patient selection and inclusion criteria.
We categorized patients into two groups, based on whether they did or did not have at least one abnormal liver chemistry result within 24 hours of the start time of the COVID encounter. Thresholds used to define liver chemistry abnormalities were ALP >105 U/L, AST >40 U/L, ALT >40 U/L, and Tbili >1.2 mg/dL (6,7). International Classification of Diseases, Ninth and Tenth Revision, Clinical Modification (ICD-9-CM and ICD-10-CM) diagnosis codes were used to identify patients who were associated with at least one ICD code indicative of any of the following conditions: diabetes mellitus (DM), chronic obstructive pulmonary disease (COPD), coronary artery disease (CAD), liver disease, history of nicotine dependence or exposure to environmental tobacco smoke, malignancy, heart failure (HF), and chronic kidney disease or end-stage renal failure (CKD-ESRD) (see Table S3). Case definitions and codes used to denote the primary outcomes of in-hospital death, mechanical ventilation, and hospital LOS are in Table S4. Hospital LOS was grouped into two categories: LOS ≤3 days and LOS >3 days. We used raw or calculated body mass index (BMI) measurements documented within 6 months of the COVID encounter to categorize patients as obese (BMI ≥30) or not obese (BMI <30). Medication data were used to identify patients who had been prescribed an oral proton pump inhibitor (PPI) or H2 blocker during the 60 days prior to the COVID encounter. Medication names and search terms for PPIs and H2 blockers are in Tables S5,S6. Serum albumin laboratory results were used to identify patients who did and did not have decreased serum albumin (<3.5 g/dL) (see Table S2). Patients with missing data for any of the covariates or outcomes of interest were excluded from analysis.
Statistical analysis
Analyses were conducted on Cerner’s HealtheDataLab platform, a data science ecosystem hosted by Amazon Web Services (AWS) (8). The analytics engine powering the HealtheDataLab platform is Apache Spark version 2.4.4 (Apache Software Foundation, Wilmington, DE). Data extraction and analyses were performed using Python version 3.7.6 (Python Software Foundation) and R version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria).
We used Pearson’s chi-squared tests to detect significant differences in categorical variables for patients with and without abnormal liver chemistries. Welch’s two-sample t-tests were used to make comparisons of liver chemistry (ALT, AST, ALP, Tbili) and serum albumin results. All other continuous variables were analyzed using independent samples t-tests. A P value of <0.05 was considered significant.
Retrospective analysis of two patient groups (i.e., patients with and without abnormal liver chemistry results) can be adversely impacted by dissimilarities existing between the two groups at baseline. Propensity score (PS) analysis is used to reduce this potential confounding. In PS analysis, a propensity score is calculated for each patient via logistic regression. In this study, the propensity score reflects the conditional probability that a given patient has at least one abnormal liver chemistry result (ALT, AST, ALP, and/or Tbili) given the patient’s values for the other covariates used to calculate the PS (9). The distribution of the covariates is equal for patients in each group who have the same propensity score (10). By using propensity scores to match patients across each group, we control for baseline differences in the covariates and can more accurately estimate the effect of the variable of interest (i.e., abnormalities in liver chemistries) on the outcomes of interest (9).
To perform PS analysis, we first categorized each patient into one of two groups: (I) patients with one or more abnormal liver chemistry results (n=8,883), or (II) patients with no abnormal liver chemistry results (n=5,255). Covariates used to calculate a propensity score for each patient included the following: age (as a quantitative variable); gender; race/ethnicity [grouped into the following categories: (I) non-Hispanic White; (II) non-Hispanic Black; (III) Hispanic; and (IV) Asian American/Pacific Islander, American Indian/Alaska Native, or non-Hispanic Other]; blood type; BMI (as a quantitative variable); history of smoking or exposure to tobacco smoke; presence of each of seven different comorbidities, including DM, COPD, CAD, malignancy, HF, CKD-ESRD, and pre-existing liver disease; serum albumin; recent use of an H2 blocker; and recent use of a PPI. We used MatchIt to calculate a PS for each patient via logistic regression. We then used the full matching algorithm in MatchIt to match each patient with at least one abnormal liver chemistry result with one or more individuals who had no abnormal liver chemistry results. Patients who did not have liver chemistry abnormalities were assigned weights proportional to those patients who did have liver chemistry abnormalities. To optimize the matching process and reduce bias, we specified a caliper of 0.25 of the standard deviation of the logit (9). After using propensity scores to match patients across each group, we used stepwise (bidirectional elimination) logistic regression to estimate the effect of abnormalities in liver chemistries on the risk of all three outcomes (death, intubation, and longer hospital LOS). We then used rbounds to conduct Rosenbaum’s sensitivity tests to quantify the sensitivity of our findings to hidden bias resulting from unmeasured covariates.
Lastly, we conducted separate analyses that used stepwise (bidirectional elimination) logistic regression, adjusting only for elevations in each of the four liver chemistries, to analyze associations between elevated liver chemistries and risk of experiencing each outcome.
Results
Baseline demographics and characteristics of liver injury
Among 102,438 patients who had an encounter associated with laboratory-confirmed COVID-19 infection, a total of 14,138 patients met inclusion criteria, had no missing data for the covariates of interest, and had no missing data for all three outcomes (see Figure 1). Demographic and clinical characteristics of these patients are outlined in Table 1.
Table 1
Characteristic | Overall | Patients with normal liver chemistries | Patients with ≥1 abnormal liver chemistries | P value |
---|---|---|---|---|
Total No. | 14,138 | 5,255 | 8,883 | |
Age, years [median (IQR)] | 62.0 (47.0, 73.0) | 64.0 (48.0, 76.0) | 61.0 (47.0, 72.0) | <0.001 |
Age categories, years, n (%) | ||||
0–24 | 686 (4.9) | 221 (4.2) | 465 (5.2) | 0.007 |
25–49 | 3,479 (24.6) | 1,278 (24.3) | 2,201 (24.8) | 0.555 |
50–74 | 6,964 (49.3) | 2,409 (45.8) | 4,555 (51.3) | <0.001 |
75+ | 3,009 (21.2) | 1,347 (25.7) | 1,662 (18.7) | <0.001 |
Male sex, n (%) | 6,939 (49.1) | 2,221 (42.3) | 4,718 (53.1) | <0.001 |
Race/ethnicity, n (%) | ||||
Non-Hispanic White | 4,708 (33.3) | 2,026 (38.6) | 2,682 (30.2) | <0.001 |
Hispanic White | 4,437 (31.4) | 1,388 (26.4) | 3,049 (34.3) | <0.001 |
Non-Hispanic Black | 2,436 (17.2) | 986 (18.8) | 1,450 (16.3) | <0.001 |
Hispanic Other | 1,486 (10.5) | 523 (10.0) | 963 (10.8) | 0.102 |
American Indian or Alaska Native | 447 (3.2) | 116 (2.2) | 331 (3.7) | <0.001 |
Asian or Pacific Islander | 334 (2.4) | 92 (1.8) | 242 (2.7) | <0.001 |
Non-Hispanic Other | 239 (1.7) | 107 (2.0) | 132 (1.5) | 0.017 |
Hispanic Black | 51 (0.4) | 17 (0.3) | 34 (0.4) | 0.672 |
Blood type, n (%) | ||||
O+ | 6,616 (46.8) | 2,361 (44.9) | 4,255 (47.9) | 0.001 |
A+ | 4,380 (31.0) | 1,693 (32.2) | 2,687 (30.2) | 0.015 |
B+ | 1,581 (11.2) | 584 (11.1) | 997 (11.2) | 0.862 |
O- | 578 (4.1) | 220 (4.2) | 358 (4.0) | 0.682 |
A- | 517 (3.7) | 209 (4.0) | 308 (3.5) | 0.130 |
AB+ | 353 (2.5) | 150 (2.9) | 203 (2.3) | 0.041 |
B- | 106 (0.7) | 36 (0.7) | 70 (0.8) | 0.559 |
AB- | 7 (0.0) | 2 (0.0) | 5 (0.1) | 0.936 |
BMI, kg/m2 [mean (SD)] | 30.7 (8.1) | 30.3 (7.9) | 31.0 (8.2) | <0.001 |
BMI categories, kg/m2, n (%) | ||||
<18.5 | 422 (3.0) | 163 (3.1) | 259 (2.9) | 0.564 |
18.5–24.9 | 2,564 (18.1) | 1,075 (20.5) | 1,489 (16.8) | <0.001 |
25.0–29.9 | 4,104 (29.0) | 1,507 (28.7) | 2,597 (29.2) | 0.492 |
30.0–34.9 | 3,286 (23.2) | 1,207 (23.0) | 2,079 (23.4) | 0.567 |
35.0–39.9 | 1,876 (13.3) | 649 (12.4) | 1,227 (13.8) | 0.014 |
≥40 | 1,886 (13.3) | 654 (12.4) | 1,232 (13.9) | 0.017 |
History of nicotine dependence or exposure to tobacco smoke, n (%) | 4,962 (35.1) | 1,924 (36.6) | 3,038 (34.2) | 0.004 |
Diabetes mellitus, n (%) | 6,620 (46.8) | 2,450 (46.6) | 4,170 (46.9) | 0.724 |
Chronic obstructive pulmonary disease, n (%) | 2,060 (14.6) | 921 (17.5) | 1,139 (12.8) | <0.001 |
Coronary artery disease, n (%) | 4,086 (28.9) | 1,606 (30.6) | 2,480 (27.9) | 0.001 |
Malignancy, n (%) | 2,092 (14.8) | 858 (16.3) | 1,234 (13.9) | <0.001 |
Heart failure, n (%) | 3,023 (21.4) | 1190 (22.6) | 1,833 (20.6) | 0.005 |
Chronic kidney disease/end-stage renal disease, n (%) | 3,859 (27.3) | 1510 (28.7) | 2,349 (26.4) | 0.003 |
Liver disease, n (%) | 2,485 (17.6) | 668 (12.7) | 1,817 (20.5) | <0.001 |
Alcohol-induced liver disease, n (%) | 209 (1.5) | 18 (0.3) | 191 (2.2) | <0.001 |
Viral liver disease, n (%) | 465 (3.3) | 144 (2.7) | 321 (3.6) | 0.006 |
NASH/fatty liver disease, n (%) | 1,085 (7.7) | 288 (5.5) | 797 (9.0) | <0.001 |
Other liver disease, n (%) | 1,468 (10.4) | 349 (6.6) | 1,119 (12.6) | <0.001 |
Albumin, g/dL [mean (SD)] | 3.3 (0.7) | 3.4 (0.6) | 3.2 (0.7) | <0.001 |
Albumin <3.5 g/dL, n (%) | 8,407 (59.5) | 2,687 (51.1) | 5,720 (64.4) | <0.001 |
ALT, U/L [median (IQR)] | 28.0 (18.0, 47.0) | 19.0 (13.0, 25.0) | 41.0 (25.0, 64.0) | <0.001 |
AST, U/L [median (IQR)] | 38.0 (25.0, 60.0) | 25.0 (19.0, 32.0) | 52.0 (39.0, 78.0) | <0.001 |
ALP, U/L [median (IQR)] | 81.0 (64.0, 109.0) | 71.0 (58.0, 85.0) | 95.0 (68.0, 132.0) | <0.001 |
Total bilirubin, mg/dL [median (IQR)] | 0.5 (0.4, 0.8) | 0.4 (0.3, 0.6) | 0.6 (0.4, 0.9) | <0.001 |
H2 blocker within 60 days, n (%) | 98 (0.7) | 34 (0.6) | 64 (0.7) | 0.686 |
PPI within 60 days, n (%) | 346 (2.4) | 147 (2.8) | 199 (2.2) | 0.044 |
Length of stay, days [mean (SD)] | 11.0 (12.9) | 9.2 (11.9) | 12.0 (13.4) | <0.001 |
Deceased, n (%) | 2,363 (16.7) | 557 (10.6) | 1,806 (20.3) | <0.001 |
Intubated, n (%) | 2,306 (16.3) | 557 (10.6) | 1,749 (19.7) | <0.001 |
Length of stay >3 days, n (%) | 10,267 (72.6) | 3,472 (66.1) | 6,795 (76.5) | <0.001 |
ALP, alkaline phosphatase; ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; g/dL, grams per deciliter; H2 blocker, histamine 2 receptor agonist; IQR, interquartile range; kg/m2, kilograms per meter squared; mg/dL, milligrams per deciliter; NASH, nonalcoholic steatohepatitis; PPI, proton pump inhibitor; SD, standard deviation; U/L, units per liter.
We compared the patient characteristics in the group with liver abnormalities to those in the group lacking the abnormalities. Median age of patients with derangements in liver chemistries was somewhat lower [61.0 (IQR, 47.0–72.0) vs. 64.0 (IQR, 48.0–76.0) years; P<0.001]. Male gender was more frequently associated with abnormalities in liver chemistries (53.1% vs. 42.3%, P<0.001). Non-Hispanic White individuals accounted for the greatest number of total patients (n=4,708, 33.3%) in the study, followed by Hispanic White (n=4,437, 31.4%), Non-Hispanic Black (n=2,436, 17.2%), Hispanic Other (n=1,486, 10.5%), American Indian or Alaska Native (n=447, 3.2%), Asian or Pacific Islander (n=334, 2.4%), Non-Hispanic Other (n=239, 1.7%), and Hispanic Black (n=51, 0.4%). Figure 2 shows the proportion of patients in each racial/ethnic group who had liver abnormalities. Blood type O+ was the most common blood group type in these patients (n=6,616, 46.8%), followed by blood types A+ (n=4,380, 31.0%) and B+ (n=1,581, 11.2%).
We also evaluated the presence of comorbidities such as obesity (defined as BMI ≥30 kg/m2), history of nicotine dependence or exposure to environmental tobacco smoke, DM, COPD, CAD, malignancy, HF, CKD-ESRD, and liver disease; other factors included recent prescriptions for a PPI or H2 blocker (see Table 1). Figure 3 shows the proportion of patients with and without liver abnormalities who were diagnosed with each comorbid condition. The comorbidity most commonly seen in this patient population was obesity (n=7,048, 49.8%), followed by diabetes (n=6,620, 46.8%), history of nicotine dependence or exposure to tobacco smoke (n=4,962, 35.1%), CAD (n=4,086, 28.9%), CKD/ESRD (n=3,859, 27.3%), HF (n=3,023, 21.4%), liver disease (n=2,485, 17.6%), malignancy (n=2,092, 14.8%) and COPD (n=2,060, 14.6%). The only co-morbidities that were more prevalent in the group with the abnormal liver chemistries were obesity (51.1% vs. 47.8%, P<0.001), DM (46.9% vs. 46.6%, P=0.724), and liver disease (20.5% vs. 12.7%, P<0.001); however, the difference in DM prevalence between patients with and without abnormal liver chemistries was not significant. Mean albumin was slightly lower in the group with abnormal labs (3.2 vs. 3.4 mg/dL, P<0.001). Prescriptions for H2 blockers were documented slightly more often in those with abnormal liver chemistries (0.7% vs. 0.6%, P=0.686), though the difference was not significant. PPIs were prescribed less frequently in those with abnormal liver chemistries (2.2% vs. 2.8%, P=0.044).
The liver chemistries that were included in the analysis were ALT, AST, ALP and Tbili. A total of 37.2% (5,255/14,138) of patients had normal liver chemistries, and 62.8% (8,883/14,138) of patients had one or more abnormal chemistries, during the 24-hour before or after the documented start time of their COVID encounter. Among patients with abnormal liver chemistries (n=8,883), we calculated the R ratio and observed the cholestatic pattern in 74.1% of patients (n=6,580), mixed in 21.5% of patients (n=1,909) and hepatocellular in 4.4% of patients (n=394) (see Figure 4).
Variables with known associations with prognosis
Prior to conducting PS analysis, we used stepwise (bidirectional elimination) multiple logistic regression to analyze associations between gender, obesity, history of smoking or exposure to tobacco smoke, pre-existing liver disease, decreased albumin, and recent PPI use, and risk for death, intubation, and increased hospital LOS (see Table 2). These factors have been identified in previous studies as having significant associations with adverse patient outcomes in COVID-19 infection.
Table 2
Outcome | Characteristic | Unadjusted | Adjusteda | Adjustedb | |||||
---|---|---|---|---|---|---|---|---|---|
Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | ||||
Mortality | Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | |||||
Gender, male | 1.49 (1.39–1.59) | <0.001 | 1.45 (1.35–1.55) | <0.001 | 1.45 (1.35–1.55) | <0.001 | |||
Healthy BMI (BMI <30 kg/m2) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Obese (BMI ≥30 kg/m2) | 0.87 (0.81–0.94) | <0.001 | 0.92 (0.85–0.99) | 0.033 | 0.92 (0.85–0.99) | 0.033 | |||
No history of nicotine dependence or exposure | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Nicotine dependence or exposure† | 1.17 (1.09–1.26) | <0.001 | 1.07 (0.99–1.15) | 0.107 | 1.07 (0.99–1.15) | 0.107 | |||
No liver disease | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Liver disease, any type | 1.47 (1.35–1.59) | <0.001 | 1.41 (1.29–1.53) | <0.001 | 1.41 (1.29–1.53) | <0.001 | |||
Normoalbuminemia (≥3.5 mg/dL) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Hypoalbuminemia (<3.5 mg/dL) | 2.00 (1.87–2.14) | <0.001 | 1.96 (1.83–2.10) | <0.001 | 1.96 (1.83–2.10) | <0.001 | |||
No PPI prescription documented in the last 60 days | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
PPI prescription documented in the last 60 days | 1.50 (1.24–1.80) | <0.001 | 1.41 (1.16–1.71) | 0.001 | 1.41 (1.16–1.71) | 0.001 | |||
Intubation | Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | |||||
Gender, male | 1.55 (1.45–1.66) | <0.001 | 1.57 (1.47–1.68) | <0.001 | 1.57 (1.47–1.68) | <0.001 | |||
Healthy BMI (BMI <30 kg/m2) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Obese (BMI ≥30 kg/m2) | 1.17 (1.09–1.26) | <0.001 | 1.24 (1.15–1.33) | <0.001 | 1.24 (1.15–1.33) | <0.001 | |||
No history of nicotine dependence or exposure | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Nicotine dependence or exposure† | 0.98 (0.91–1.06) | 0.656 | 0.90 (0.83–0.98) | 0.013 | 0.90 (0.83–0.98) | 0.010 | |||
No liver disease | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Liver disease, any type | 1.41 (1.29–1.53) | <0.001 | 1.39 (1.27–1.51) | <0.001 | 1.39 (1.27–1.51) | <0.001 | |||
Normoalbuminemia (≥3.5 mg/dL) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Hypoalbuminemia (<3.5 mg/dL) | 1.74 (1.62–1.87) | <0.001 | 1.73 (1.61–1.85) | <0.001 | 1.73 (1.60–1.85) | <0.001 | |||
No PPI prescription documented in the last 60 days | 1 (reference) | 1 (reference) | (dropped)‡ | – | |||||
PPI prescription documented in the last 60 days | 0.87 (0.66–1.12) | 0.274 | 0.83 (0.63–1.08) | 0.169 | (dropped)‡ | – | |||
LOS >3 days | Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | |||||
Gender, male | 1.13 (1.12–1.15) | <0.001 | 1.13 (1.12–1.15) | <0.001 | 1.13 (1.12–1.15) | <0.001 | |||
Healthy BMI (BMI <30 kg/m2) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Obese (BMI ≥30 kg/m2) | 1.01 (0.99–1.03) | 0.558 | 1.03 (1.01–1.05) | 0.003 | 1.03 (1.01–1.05) | 0.003 | |||
No history of nicotine dependence or exposure | 1 (reference) | 1 (reference) | (dropped)‡ | – | |||||
Nicotine dependence or exposure† | 1.02 (1.00–1.04) | 0.028 | 1.00 (0.98–1.03) | 0.729 | (dropped)‡ | – | |||
No liver disease | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Liver disease, any type | 1.06 (1.03–1.08) | <0.001 | 1.04 (1.01–1.07) | 0.003 | 1.04 (1.01–1.07) | 0.003 | |||
Normoalbuminemia (≥3.5 mg/dL) | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Hypoalbuminemia (<3.5 mg/dL) | 1.21 (1.20–1.22) | <0.001 | 1.21 (1.20–1.22) | <0.001 | 1.21 (1.20–1.22) | <0.001 | |||
No PPI prescription documented in the last 60 days | 1 (reference) | 1 (reference) | (dropped)‡ | – | |||||
PPI prescription documented in the last 60 days | 0.98 (0.91–1.05) | 0.603 | 0.96 (0.88–1.03) | 0.242 | (dropped)‡ | – |
a, Logistic regression without variable selection; b, Stepwise (bidirectional elimination) logistic regression; †, refers to current use or past history of dependency on nicotine products and/or documented exposure to environmental tobacco smoke; ‡, term was dropped from stepwise logistic regression model results. BMI, body mass index; CI, confidence interval; kg/m2, kilograms per meter squared; LOS, length of stay; mg/dL, milligrams per deciliter; PPI, proton pump inhibitor.
Mortality was significantly associated with male gender [adjusted relative risk (ARR), 1.45; 95% CI: 1.35–1.55; P<0.001], obesity (ARR, 0.92; 95% CI: 0.85–0.99; P=0.033), pre-existing liver disease (ARR, 1.41; 95% CI: 1.29–1.53; P<0.001), decreased serum albumin (ARR, 1.96; 95% CI 1.83–2.10; P<0.001) and PPI use (ARR, 1.41; 95% CI: 1.16–1.71; P=0.001).
Intubation was significantly associated with male gender (ARR, 1.57; 95% CI: 1.47–1.68; P<0.001), obesity (ARR, 1.24; 95% CI: 1.15–1.33; P<0.001), history of smoking or exposure to tobacco smoke (ARR, 0.90; 95% CI: 0.83–0.98; P=0.010), pre-existing liver disease (ARR 1.39; 95% CI: 1.27–1.51; P<0.001), and decreased serum albumin (ARR 1.73; 95% CI: 1.60–1.85; P<0.001). Hospital LOS >3 days was significantly associated with male gender (ARR, 1.13; 95% CI: 1.12–1.15; P<0.001), obesity (ARR, 1.03; 95% CI: 1.01–1.05; P=0.003), pre-existing liver disease (ARR, 1.04; 95% CI: 1.01–1.07; P=0.003), and decreased serum albumin (ARR, 1.21; 95% CI: 1.20–1.22; P<0.001).
Propensity score analysis
The full matching algorithm matched 8,879 (out of 8,883) individuals with abnormal liver chemistries with all 5,255 individuals who had only normal liver chemistries. In the matched cohort, the standardized differences in the covariates were all <0.06, indicating robust matching (9) (see Table S7). After matching, we used weighted logistic regression to analyze the impact of abnormalities in liver chemistries on the risk of death, intubation, and prolonged hospital LOS in patients with and without abnormalities in liver chemistries. Holding all other predictors fixed, the estimated risk of death was 70% higher for patients with at least one abnormal liver chemistry (ARR, 1.70; 95% CI: 1.58–1.82; P<0.001; see Table 3); and estimated risk of intubation was 44% higher for these patients (ARR, 1.44; 95% CI: 1.33–1.54; P<0.001; see Table 3). Estimated risk of hospital LOS >3 days was also slightly increased for patients with at least one abnormal liver chemistry (ARR, 1.08; 95% CI: 1.06–1.10; P<0.001; see Table 3).
Table 3
Outcome | Characteristic | Unadjusted | Adjusteda | Adjustedb | |||||
---|---|---|---|---|---|---|---|---|---|
Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | ||||
Mortality | Age† | 1.04 (1.04–1.04) | <0.001 | 1.04 (1.04–1.04) | <0.001 | 1.04 (1.04–1.04) | <0.001 | ||
Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Gender, male | 1.45 (1.35–1.55) | <0.001 | 1.44 (1.34–1.55) | <0.001 | 1.44 (1.34–1.55) | <0.001 | |||
Race/ethnicity, non–Hispanic White | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Race/ethnicity, non–Hispanic Black | 0.69 (0.61–0.79) | <0.001 | 0.88 (0.77–1.00) | 0.047 | 0.89 (0.78–1.01) | 0.065 | |||
Race/ethnicity, AAPI/AIAN/non–Hispanic Other | 0.89 (0.77–1.03) | 0.112 | 1.13 (0.97–1.31) | 0.128 | 1.14 (0.98–1.32) | 0.099 | |||
Race/ethnicity, Hispanic | 0.91 (0.84–0.99) | 0.033 | 1.12 (1.02–1.23) | 0.016 | 1.13 (1.03–1.23) | 0.009 | |||
Diabetes mellitus | 1.63 (1.53–1.73) | <0.001 | 1.23 (1.13–1.33) | <0.001 | 1.23 (1.14–1.33) | <0.001 | |||
Chronic obstructive pulmonary disease | 1.43 (1.31–1.56) | <0.001 | 0.98 (0.88–1.10) | 0.781 | (dropped)§ | – | |||
Coronary artery disease | 1.79 (1.68–1.91) | <0.001 | 1.12 (1.03–1.23) | 0.011 | 1.14 (1.04–1.24) | 0.004 | |||
Liver disease | 1.41 (1.31–1.52) | <0.001 | 1.42 (1.31–1.55) | <0.001 | 1.43 (1.32–1.55) | <0.001 | |||
BMI† | 0.99 (0.98–0.99) | <0.001 | 1.01 (1.01–1.02) | <0.001 | 1.01 (1.01–1.02) | <0.001 | |||
PPI prescription in the last 60 days | 1.41 (1.15–1.70) | 0.001 | 1.03 (0.81–1.29) | 0.812 | (dropped)§ | – | |||
H2 blocker prescription in the last 60 days | 1.45 (0.98–2.03) | 0.059 | 1.09 (0.69–1.64) | 0.699 | (dropped)§ | – | |||
Serum albumin† | 0.53 (0.50–0.56) | <0.001 | 0.61 (0.57–0.65) | <0.001 | 0.61 (0.57–0.65) | <0.001 | |||
Blood type, O+ | 1 (reference) | 1 (reference) | (dropped)§ | – | |||||
Blood type, A+ | 0.96 (0.88–1.05) | 0.362 | 0.93 (0.85–1.02) | 0.107 | (dropped)§ | – | |||
Blood type, B+ | 0.94 (0.83–1.06) | 0.326 | 1.00 (0.88–1.14) | 0.978 | (dropped)§ | – | |||
Blood type, O– | 0.88 (0.72–1.06) | 0.184 | 0.88 (0.71–1.08) | 0.234 | (dropped)§ | – | |||
Blood type, A– | 1.04 (0.85–1.27) | 0.681 | 0.92 (0.73–1.14) | 0.458 | (dropped)§ | – | |||
Blood type, AB+ | 1.04 (0.82–1.30) | 0.769 | 1.04 (0.80–1.32) | 0.781 | (dropped)§ | – | |||
Blood type, B– | 0.79 (0.47–1.28) | 0.345 | 0.63 (0.36–1.07) | 0.092 | (dropped)§ | – | |||
Blood type, AB– | 1.81 (0.45–4.20) | 0.368 | 0.87 (0.16–3.07) | 0.863 | (dropped)§ | – | |||
Malignancy | 1.43 (1.31–1.56) | <0.001 | 1.06 (0.95–1.17) | 0.299 | (dropped)§ | – | |||
Heart failure | 1.67 (1.55–1.79) | <0.001 | 1.06 (0.96–1.16) | 0.263 | (dropped)§ | – | |||
CKD/ESRD | 1.90 (1.79–2.02) | <0.001 | 1.31 (1.20–1.42) | <0.001 | 1.32 (1.22–1.43) | <0.001 | |||
Nicotine dependence or exposure‡ | 1.11 (1.03–1.19) | 0.008 | 0.80 (0.73–0.88) | <0.001 | 0.80 (0.74–0.88) | <0.001 | |||
Abnormal liver chemistry result present | 1.51 (1.40–1.62) | <0.001 | 1.70 (1.58–1.83) | <0.001 | 1.70 (1.58–1.82) | <0.001 | |||
Intubation | Age† | 1.01 (1.01–1.01) | <0.001 | 1.01 (1.01–1.01) | <0.001 | 1.01 (1.01–1.01) | <0.001 | ||
Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Gender, male | 1.51 (1.41–1.61) | <0.001 | 1.57 (1.47–1.69) | <0.001 | 1.57 (1.47–1.69) | <0.001 | |||
Race/ethnicity, non–Hispanic White | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Race/ethnicity, non–Hispanic Black | 0.89 (0.79–1.00) | 0.050 | 0.91 (0.80–1.03) | 0.128 | 0.91 (0.80–1.03) | 0.126 | |||
Race/ethnicity, AAPI/AIAN/non–Hispanic Other | 1.73 (1.55–1.93) | <0.001 | 1.59 (1.41–1.79) | <0.001 | 1.59 (1.41–1.79) | <0.001 | |||
Race/ethnicity, Hispanic | 1.04 (0.96–1.13) | 0.360 | 0.94 (0.85–1.03) | 0.185 | 0.94 (0.85–1.03) | 0.177 | |||
Diabetes mellitus | 1.61 (1.51–1.72) | <0.001 | 1.44 (1.33–1.55) | <0.001 | 1.44 (1.33–1.55) | <0.001 | |||
Chronic obstructive pulmonary disease | 1.03 (0.92–1.14) | 0.605 | 1.00 (0.89–1.13) | 0.954 | (dropped)§ | – | |||
Coronary artery disease | 1.15 (1.06–1.24) | <0.001 | 0.92 (0.83–1.01) | 0.090 | 0.92 (0.83–1.01) | 0.090 | |||
Liver disease | 1.36 (1.25–1.47) | <0.001 | 1.34 (1.24–1.46) | <0.001 | 1.34 (1.24–1.46) | <0.001 | |||
BMI† | 1.01 (1.01–1.02) | <0.001 | 1.02 (1.02–1.03) | <0.001 | 1.02 (1.02–1.03) | <0.001 | |||
PPI prescription in the last 60 days | 0.74 (0.56–0.98) | 0.034 | 0.65 (0.48–0.87) | 0.004 | 0.65 (0.48–0.87) | 0.004 | |||
H2 blocker prescription in the last 60 days | 1.16 (0.75–1.72) | 0.505 | 1.09 (0.68–1.67) | 0.705 | (dropped)§ | – | |||
Serum albumin† | 0.60 (0.57–0.64) | <0.001 | 0.63 (0.59–0.67) | <0.001 | 0.63 (0.59–0.67) | <0.001 | |||
Blood type, O+ | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Blood type, A+ | 0.91 (0.84–0.99) | 0.030 | 0.93 (0.85–1.02) | 0.107 | 0.93 (0.85–1.02) | 0.108 | |||
Blood type, B+ | 0.89 (0.79–1.01) | 0.066 | 0.93 (0.82–1.05) | 0.250 | 0.93 (0.82–1.05) | 0.252 | |||
Blood type, O– | 0.73 (0.59–0.90) | 0.003 | 0.82 (0.66–1.02) | 0.070 | 0.82 (0.66–1.02) | 0.070 | |||
Blood type, A– | 0.79 (0.62–0.99) | 0.036 | 0.80 (0.63–1.01) | 0.064 | 0.80 (0.63–1.01) | 0.065 | |||
Blood type, AB+ | 0.84 (0.65–1.09) | 0.193 | 0.83 (0.63–1.08) | 0.176 | 0.83 (0.63–1.08) | 0.175 | |||
Blood type, B– | 0.52 (0.28–0.95) | 0.033 | 0.48 (0.25–0.90) | 0.020 | 0.48 (0.25–0.90) | 0.020 | |||
Blood type, AB– | 0.85 (0.11–3.55) | 0.865 | 0.66 (0.08–3.18) | 0.667 | 0.66 (0.08–3.18) | 0.666 | |||
Malignancy | 0.76 (0.68–0.86) | <0.001 | 0.72 (0.63–0.81) | 0.000 | 0.72 (0.63–0.81) | <0.001 | |||
Heart failure | 1.30 (1.20–1.41) | <0.001 | 1.13 (1.02–1.24) | 0.019 | 1.13 (1.02–1.24) | 0.017 | |||
CKD/ESRD | 1.33 (1.23–1.43) | <0.001 | 1.09 (1.00–1.20) | 0.055 | 1.10 (1.00–1.20) | 0.054 | |||
Nicotine dependence or exposure‡ | 0.86 (0.79–0.93) | <0.001 | 0.75 (0.69–0.83) | <0.001 | 0.75 (0.69–0.82) | <0.001 | |||
Abnormal liver chemistry result present | 1.36 (1.27–1.47) | <0.001 | 1.43 (1.33–1.54) | <0.001 | 1.44 (1.33–1.54) | <0.001 | |||
LOS >3 days | Age† | 1.01 (1.01–1.01) | <0.001 | 1.01 (1.01–1.01) | <0.001 | 1.01 (1.01–1.01) | <0.001 | ||
Gender, female | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Gender, male | 1.13 (1.11–1.14) | <0.001 | 1.11 (1.10–1.13) | <0.001 | 1.11 (1.10–1.13) | <0.001 | |||
Race/ethnicity, non–Hispanic White | 1 (reference) | 1 (reference) | 1 (reference) | ||||||
Race/ethnicity, non–Hispanic Black | 0.90 (0.86–0.94) | <0.001 | 0.96 (0.92–0.99) | 0.018 | 0.96 (0.93–1.00) | 0.046 | |||
Race/ethnicity, AAPI/AIAN/non–Hispanic Other | 1.04 (1.00–1.08) | 0.066 | 1.05 (1.01–1.10) | 0.020 | 1.06 (1.02–1.10) | 0.009 | |||
Race/ethnicity, Hispanic | 0.97 (0.95–1.00) | 0.039 | 1.00 (0.97–1.03) | 0.901 | 1.00 (0.97–1.03) | 0.879 | |||
Diabetes mellitus | 1.17 (1.16–1.18) | <0.001 | 1.10 (1.08–1.12) | <0.001 | 1.10 (1.08–1.12) | <0.001 | |||
Chronic obstructive pulmonary disease | 1.07 (1.04–1.10) | <0.001 | 0.97 (0.93–1.01) | 0.208 | 0.97 (0.93–1.01) | 0.192 | |||
Coronary artery disease | 1.12 (1.10–1.14) | <0.001 | 0.98 (0.95–1.01) | 0.195 | 0.98 (0.95–1.01) | 0.213 | |||
Liver disease | 1.06 (1.04–1.08) | <0.001 | 1.04 (1.02–1.07) | 0.002 | 1.04 (1.02–1.07) | 0.002 | |||
BMI† | 1.00 (1.00–1.00) | 0.025 | 1.01 (1.00–1.01) | <0.001 | 1.01 (1.00–1.01) | <0.001 | |||
PPI prescription in the last 60 days | 0.93 (0.85–0.99) | 0.034 | 0.81 (0.73–0.90) | <0.001 | 0.81 (0.73–0.90) | <0.001 | |||
H2 blocker prescription in the last 60 days | 0.97 (0.82–1.09) | 0.631 | 0.92 (0.75–1.06) | 0.312 | (dropped)§ | – | |||
Serum albumin † | 0.70 (0.68–0.72) | <0.001 | 0.76 (0.74–0.79) | <0.001 | 0.76 (0.74–0.78) | <0.001 | |||
Blood type, O+ | 1 (reference) | 1 (reference) | (dropped)§ | – | |||||
Blood type, A+ | 0.99 (0.97–1.02) | 0.616 | 0.99 (0.97–1.02) | 0.672 | (dropped)§ | – | |||
Blood type, B+ | 0.99 (0.95–1.02) | 0.458 | 1.02 (0.98–1.06) | 0.279 | (dropped)§ | – | |||
Blood type, O– | 0.91 (0.85–0.97) | 0.001 | 0.92 (0.86–0.98) | 0.009 | (dropped)§ | – | |||
Blood type, A– | 1.00 (0.94–1.06) | 0.920 | 0.98 (0.91–1.04) | 0.579 | (dropped)§ | – | |||
Blood type, AB+ | 1.05 (0.98–1.11) | 0.193 | 1.04 (0.97–1.11) | 0.256 | (dropped)§ | – | |||
Blood type, B– | 0.97 (0.82–1.09) | 0.616 | 0.94 (0.78–1.07) | 0.391 | (dropped) § | – | |||
Blood type, AB– | 1.38 (0.00–1.38) | 0.929 | 1.38 (0.00–1.38) | 0.931 | (dropped)§ | – | |||
Malignancy | 1.06 (1.03–1.08) | <0.001 | 0.97 (0.94–1.01) | 0.120 | 0.97 (0.94–1.01) | 0.125 | |||
Heart failure | 1.13 (1.11–1.15) | <0.001 | 1.02 (0.99–1.05) | 0.207 | 1.02 (0.99–1.05) | 0.228 | |||
CKD/ESRD | 1.15 (1.13–1.17) | <0.001 | 1.05 (1.02–1.07) | 0.003 | 1.05 (1.02–1.07) | 0.003 | |||
Nicotine dependence or exposure‡ | 1.00 (0.98–1.03) | 0.660 | 0.93 (0.90–0.96) | <0.001 | 0.93 (0.90–0.96) | <0.001 | |||
Abnormal liver chemistry result present | 1.07 (1.05–1.09) | <0.001 | 1.08 (1.06–1.10) | <0.001 | 1.08 (1.06–1.10) | <0.001 |
a, Logistic regression without variable selection; b, Stepwise (bidirectional elimination) logistic regression; †, analyzed as a continuous variable. Odds ratio interpretations: (I) a one–unit change in age is equivalent to a change of 1 year, (II) a one–unit change in BMI is equivalent to a change of 1 kg/m2, and (III) a one–unit change in serum albumin is equivalent to a change of 1 mg/dL; ‡, refers to current use or past history of dependency on nicotine products and/or documented exposure to environmental tobacco smoke; §, term was dropped from stepwise logistic regression model results. AAPI, Asian American or Pacific Islander; AIAN, American Indian or Alaska Native; BMI, body mass index; CI, confidence interval; CKD/ESRD, Chronic kidney disease/End–stage renal disease; H2 blocker, histamine 2 receptor agonist; LOS, length of stay; PPI, proton pump inhibitor.
Treatment effect models were analyzed for sensitivity to hidden, or unmeasured, covariates via Rosenbaum’s sensitivity tests. Sensitivity analyses demonstrated that the findings pertaining to the outcomes of mortality and intubation were highly robust to hidden bias from unmeasured covariates. The treatment effect estimates for the outcomes of death and intubation were both robust to hidden bias up to Γ=11. For hospital LOS, the treatment effect became insignificant at Γ=1, indicating that the finding pertaining to hospital LOS was not robust to hidden bias from unmeasured covariates.
Effect of presence of liver injury
To further examine the impact of specific liver chemistry derangements on patient outcomes and to identify the abnormalities associated with the highest risk for poor outcomes, we used three stepwise (bidirectional elimination) multiple logistic regression models to analyze whether abnormalities in ALT, AST, ALP, and total bilirubin were associated with increased risk of death, intubation, and increased hospital LOS (see Table 3). These models, which adjusted only for elevations in the four liver chemistries, used the following thresholds to define abnormalities for these parameters: ALP >105 U/L, AST >40 U/L, ALT >40 U/L, and Tbili >1.2 mg/dL (6,7). Mortality was significantly associated with derangements in ALT (ARR, 0.69; 95% CI: 0.63–0.76; P<0.001), AST (ARR, 2.27; 95% CI: 2.12–2.42; P<0.001), and Tbili (ARR, 1.50; 95% CI: 1.34–1.67; P<0.001). Risk for intubation was significantly associated with derangements in ALT (ARR, 0.90; 95% CI: 0.82–0.98; P=0.021), AST (ARR, 2.12; 95% CI: 1.97–2.27; P<0.001), and Tbili (ARR, 1.17; 95% CI: 1.03–1.33; P=0.014). Hospital LOS >3 days was significantly associated with abnormalities in ALT (ARR, 0.93; 95% CI: 0.90–0.96; P<0.001), AST (ARR, 1.19; 95% CI: 1.18–1.21; P<0.001), and ALP (ARR, 0.96; 95% CI: 0.94–0.99; P=0.003). These findings are outlined in Table 4.
Table 4
Outcome | Characteristic | Unadjusted | Adjusteda | Adjustedb | |||||
---|---|---|---|---|---|---|---|---|---|
Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | Relative risk (95% CI) |
P value | ||||
Mortality | ALT >40 U/L | 1.20 (1.12–1.30) | <0.001 | 0.69 (0.63–0.76) | <0.001 | 0.69 (0.63–0.76) | <0.001 | ||
AST >40 U/L | 2.05 (1.92–2.17) | <0.001 | 2.27 (2.12–2.42) | <0.001 | 2.27 (2.12–2.42) | <0.001 | |||
ALP >105 U/L | 1.19 (1.10–1.29) | <0.001 | 1.08 (0.99–1.17) | 0.073 | (dropped)† | – | |||
Total bilirubin >1.2 mg/dL | 1.76 (1.59–1.94) | <0.001 | 1.50 (1.34–1.67) | <0.001 | 1.50 (1.34–1.67) | <0.001 | |||
Intubation | ALT >40 U/L | 1.42 (1.32–1.52) | <0.001 | 0.90 (0.82–0.98) | 0.022 | 0.90 (0.82–0.98) | 0.021 | ||
AST >40 U/L | 2.05 (1.93–2.18) | <0.001 | 2.12 (1.97–2.27) | <0.001 | 2.12 (1.97–2.27) | <0.001 | |||
ALP >105 U/L | 1.10 (1.02–1.19) | 0.020 | 1.00 (0.92–1.09) | 0.984 | (dropped)† | – | |||
Total bilirubin >1.2 mg/dL | 1.40 (1.24–1.57) | <0.001 | 1.17 (1.03–1.33) | 0.016 | 1.17 (1.03–1.33) | 0.014 | |||
Length of stay >3 days | ALT >40 U/L | 1.07 (1.05–1.09) | <0.001 | 0.93 (0.90–0.96) | <0.001 | 0.93 (0.90–0.96) | <0.001 | ||
AST >40 U/L | 1.17 (1.16–1.18) | <0.001 | 1.19 (1.18–1.21) | <0.001 | 1.19 (1.18–1.21) | <0.001 | |||
ALP >105 U/L | 0.99 (0.96–1.01) | 0.258 | 0.96 (0.94–0.99) | 0.003 | 0.96 (0.94–0.99) | 0.003 | |||
Total bilirubin >1.2 mg/dL | 1.06 (1.02–1.09) | 0.004 | 1.01 (0.97–1.05) | 0.595 | (dropped)† | – |
a, Logistic regression without variable selection; b, Stepwise (bidirectional elimination) logistic regression; †, term was dropped from stepwise logistic regression model results. ALP, alkaline phosphatase; ALT, alanine aminotransferase; AST, aspartate aminotransferase; CI, confidence interval; mg/dL, milligrams per deciliter; U/L, units per liter.
Discussion
Coronavirus-19 infection can affect numerous organs including the liver. The hepatotoxic effect is thought to be predominantly mediated by its interaction with the angiotensin-converting enzyme 2 (ACE2) receptor in the liver, especially on the cholangiocytes manifesting as intrahepatic cholestasis (11). Microvesicular steatosis and portal inflammation has been seen on autopsies of patients with the infection (12). Other possible mechanisms of hepatotoxicity include direct damage to the hepatocytes causing viral hepatitis and indirect damage may occur secondary to systemic immune mediated response or from medications that may be used for treatment of the infection itself (13). The degree of hepatic inflammation and injury can vary from being mild to severe and potentially acute fulminant liver injury can occur.
Male gender has been associated with an increased risk of hepatic damage (14,15). In our study, 53.1% of those with abnormalities in liver chemistries were male. Male gender was also associated with increased risk of mortality (RR 1.45), intubation (RR 1.57), and prolonged hospitalization (RR 1.13). This suggests that males may be at a higher risk for COVID-19 related liver injury; however, further evaluation is necessary to define this association.
Obesity has been linked with adverse outcomes, including increased mortality (16). It impairs the cytotoxic cellular response to viral infections by modifying the reaction of cytokines and through its effects on hormones like leptin (17). Our study indicates that obesity increases risk for intubation (RR 1.24) and prolonged hospitalization (RR 1.03). Interestingly, obesity was associated with a slightly decreased risk for mortality (RR 0.92). While the reason for this decreased mortality risk is not fully clear, it is possible that, in this patient population, obesity may impact mortality risk through interactions with other patient factors.
Smoking is another factor to consider in patients with SARS-CoV-2 infection and is associated with worse outcomes in some studies and an unexpectedly beneficial effect in others. An increased expression of the ACE2 gene in the lungs and a higher number of ACE2-producing goblet cells was found in those who have ever smoked as compared to those who have never smoked in a study performed by Cai et al. (18). A large study based in Iran saw that cigarette smoking was associated with an increased risk of mortality related to COVID-19 though a negative correlation with hospitalization (19). Another study that was carried out in Wuhan found that history of smoking was linked to progression of disease (20). Conflicting data was noted by Lippi et al. in their meta-analysis on Chinese patients that revealed active smoking was not linked to a higher risk for progression of severe disease in cases of SARS-CoV-2 infection (21). Several studies have revealed a potentially beneficial effect of smoking and this may be due to possible confounding variables (22). In our model that adjusted for gender, obesity, smoking history, presence of chronic liver disease, decreased albumin and PPI use, the association between smoking and mortality was not significant. However, smoking was linked to a reduced risk for intubation (RR 0.90). Hence further studies are needed to analyze the effects of smoking on the course of SARS-CoV-2 infection.
The presence of pre-existing liver disease was also associated with an increased relative risk of mortality (RR 1.41), intubation (RR 1.39) and prolonged hospital LOS (RR 1.04). This is an important comorbid condition to recognize due to its negative impact on infection outcomes.
Hypoalbuminemia was also associated with increased relative risk for mortality (RR 1.96), intubation (1.73), and prolonged hospitalization (RR 1.21). Serum albumin levels were lower in patients with abnormal chemistries. Low albumin levels have been shown in several studies to be linked to a more severe phenotypic presentation of COVID-19 infection (23-25). Hypoalbuminemia has also been shown to correlate with a poor prognosis including a higher probability of respiratory decompensation with need for respiratory support (20).
Recent use of PPIs (within 60 days of admission) was linked to increased relative risk of mortality (RR 1.41). This association has been reported in many other studies, with PPIs being linked to worse outcomes and increased mortality (26-28). PPI therapy causes hypochlorhydria with suppression of the effect of the gastric acid, hence promoting the survival of the virus and the resultant virulent effects (27,28). Hypochlorhydria also increases the bacterial load in the stomach and can lead to small intestinal bacterial overgrowth (27,28). In cases of micro-aspiration, the bacteria from the stomach can cause a secondary infection in the lungs and augment the chances of developing Acute Respiratory Distress Syndrome and worse outcomes including increased mortality (27,28). Indirect inhibition of nitric oxide synthase leads to thrombosis and inflammation with negative cardiovascular effects (27,28). They also exert a negative impact on neutrophil function and hence impair the immune response to the virus (27,29). In our stepwise models, PPI use was not associated with higher risk for intubation or prolonged hospitalization.
Small studies have revealed that there is an association between liver injury and worse outcomes and prognosis (7,11,13,30). To the best of our knowledge, this is the largest study to evaluate associations between elevated liver chemistries on admission and the outcomes of mortality, intubation, and duration of hospitalization. Propensity score analysis indicated that the relative risk of mortality was 1.70 times higher, and the relative risk of intubation was 1.44 times higher, in those with abnormal liver chemistries. Risk for prolonged hospitalization was also 1.08 times higher in this patient cohort, though this finding was highly sensitive to unmeasured covariates.
The study conducted by Medetalibeyoglu et al., with 554 subjects, they found that patients with a high ALT-AST ratio during their hospitalization had increased mortality and requirement for ICU level of care (7). They noted that AST/ALT ratio >1 led to worse outcomes including higher mortality (7). Lei et al. reviewed 5,771 patients in a multicenter retrospective study and evaluated levels of AST, ALP, ALT and total bilirubin during the course of their hospitalization (31). It showed that initially AST increased then the ALT in severe cases and the ALP increase was noted to be modest with mild fluctuations in total bilirubin (31). Hence they concluded that the AST elevations were linked to the highest mortality (31). A prospective study that consisted of 217 patients showed that increase in ALT, GGT (gamma-glutamyl transferase) and low albumin levels on hospital admission were all linked to severe infection and unfavorable outcomes such ICU stay and mechanical ventilation (32).
Similarly, in our study, we saw that derangements in liver chemistries on admission were associated with patient outcomes. In models that adjusted for elevations in all four chemistries, elevation in AST was associated with the highest risks for mortality, intubation, and prolonged hospitalization. In the same adjusted models, ALT elevation was associated with a decreased relative risk for mortality, intubation, and prolonged hospitalization. Elevated ALP was linked to decreased risk for prolonged hospitalization. Elevation in total bilirubin was associated with increased risks for mortality and intubation. Based on the calculation of the R ratio, the predominant pattern of liver enzyme elevation seen is this study was the cholestatic pattern.
A meta-analysis that included sixty-four studies showed that in severe cases of COVID-19 infection, the prevalence of an increased AST level was significantly higher (33). The worldwide prevalence of an increased AST was also higher as compared to ALT, ALP, total bilirubin and gamma-glutamyltransferase (33). Another large study conducted in China showed that 39.4% of patients had an AST value >40 U/L and 28.1% had ALT >40 U/L and these elevations were found more commonly in those with severe disease (30). Slightly different findings were portrayed in another meta-analysis with ALT and bilirubin levels being significantly increased and serum albumin levels being significantly decreased in those with severe infection (23). Hence though the data thus far shows that elevation in AST tends to occur more frequently and is more profound in those with severe disease, this finding needs further investigation.
As further data emerge and our understanding of the virus evolves, there will be additional clarity with regards to hepatic involvement secondary to COVID-19 infection and outcomes. Understanding the patient profile and baseline comorbid conditions is important as they impact the clinical course and are linked to severity of infection. Of all the liver chemistries, elevated AST was associated with the greatest risk for negative outcomes. However, the source of AST in the body is not solely the liver and the AST elevation can be driven potentially from muscle injury as well. The liver injury can occur from other hepatotoxic agents such as medications that are used in the treatment for the virus or other medications used during hospitalization. These factors need to be considered when evaluating patients with liver injury. A routine liver panel is obtained during admission for most hospitalized patients. If abnormalities are detected, patients should be monitored closely to follow the trend.
Limitations and strengths
The retrospective design of our study is a significant limitation. We also did not account for elevations in liver chemistries that occurred more than 24 hours after the start time of the COVID encounter. Further study is needed to analyze the impact of derangements in liver chemistries that occur more than 24 hours after admission. We also chose to exclude patients who had missing data for any of the covariates or outcomes analyzed in this study. A major strength of our study is its large, geographically diverse sample of patients with laboratory-confirmed COVID-19 infection. To address limitations inherent in retrospective cohort analysis, we used propensity score analysis to reduce confounding due to imbalance in observed covariates.
Conclusions
Using a large, nationwide cohort that includes more than 14,000 patients with laboratory-confirmed COVID-19 infection, we observed that liver chemistry abnormalities at the time of admission impact risk for mortality, intubation, and prolonged hospital LOS. The greatest risk for all three adverse outcomes is associated with elevated AST. A liver chemistry panel should be obtained with routine laboratory tests when patients test positive for COVID-19 infection. These results should be used for risk stratification. Clinicians should be cautious when derangements in liver chemistries are noted on presentation as they could be harbingers of eventual decompensation. Further studies are needed to develop more precise prediction models and scoring systems.
Acknowledgments
Funding: EMT is supported by the National Institutes of Health (No. 5T32LM012410). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://tgh.amegroups.com/article/view/10.21037/tgh-21-94/rc
Data Sharing Statement: Available at https://tgh.amegroups.com/article/view/10.21037/tgh-21-94/dss
Peer Review File: Available at https://tgh.amegroups.com/article/view/10.21037/tgh-21-94/prf
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tgh.amegroups.com/article/view/10.21037/tgh-21-94/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The University of Missouri Institutional Review Board (Federalwide Assurance ID: FWA00002876) deemed this study of de-identified patient data to not constitute human subjects research; therefore, direct informed consent from individual patients was not required (IRB #2034442).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Baloch S, Baloch MA, Zheng T, et al. The Coronavirus Disease 2019 (COVID-19) Pandemic. Tohoku J Exp Med 2020;250:271-8. [Crossref] [PubMed]
- Gavriatopoulou M, Korompoki E, Fotiou D, et al. Organ-specific manifestations of COVID-19 infection. Clin Exp Med 2020;20:493-506. [Crossref] [PubMed]
- Polack FP, Thomas SJ, Kitchin N, et al. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine. N Engl J Med 2020;383:2603-15. [Crossref] [PubMed]
- Cerner Corporation. Cerner provides access to de-identified patient data for COVID-19 research and vaccine development (news release). 2020 Apr 24 (cited 2021 Sept 16). Available online: https://www.cerner.com/newsroom/cerner-provides-access-to-de-identified-patient-data-for-covid-19-research-and-vaccine-development
- Qureshi AI, Baskett WI, Huang W, et al. Facilitating the study of relationships between COVID-19 and cardiovascular health outcomes using Cerner Real-World COVID-19 deidentified dataset. HealthCare Res J 2020;1:17-28.
- European Association for the Study of the Liver. EASL Clinical Practice Guidelines: management of cholestatic liver diseases. J Hepatol 2009;51:237-67. [Crossref] [PubMed]
- Medetalibeyoglu A, Catma Y, Senkal N, et al. The effect of liver test abnormalities on the prognosis of COVID-19. Ann Hepatol 2020;19:614-21. [Crossref] [PubMed]
- Ehwerhemuepha L, Gasperino G, Bischoff N, et al. HealtheDataLab - a cloud computing solution for data science and advanced analytics in healthcare with application to predicting multi-center pediatric readmissions. BMC Med Inform Decis Mak 2020;20:115. [Crossref] [PubMed]
- Lee J, Little TD. A practical guide to propensity score analysis for applied clinical research. Behav Res Ther 2017;98:76-90. [Crossref] [PubMed]
- Austin PC. An Introduction to Propensity Score Methods for Reducing the Effects of Confounding in Observational Studies. Multivariate Behav Res 2011;46:399-424. [Crossref] [PubMed]
- Zhang C, Shi L, Wang FS. Liver injury in COVID-19: management and challenges. Lancet Gastroenterol Hepatol 2020;5:428-30. [Crossref] [PubMed]
- Xu Z, Shi L, Wang Y, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med 2020;8:420-2. [Crossref] [PubMed]
- Gu J, Han B, Wang J. COVID-19: Gastrointestinal Manifestations and Potential Fecal-Oral Transmission. Gastroenterology 2020;158:1518-9. [Crossref] [PubMed]
- Zhang H, Liao YS, Gong J, et al. Clinical characteristics and risk factors for liver injury in COVID-19 patients in Wuhan. World J Gastroenterol 2020;26:4694-702. [Crossref] [PubMed]
- Xie H, Zhao J, Lian N, et al. Clinical characteristics of non-ICU hospitalized patients with coronavirus disease 2019 and liver injury: A retrospective study. Liver Int 2020;40:1321-6. [Crossref] [PubMed]
- Tartof SY, Qian L, Hong V, et al. Obesity and Mortality Among Patients Diagnosed With COVID-19: Results From an Integrated Health Care Organization. Ann Intern Med 2020;173:773-81. [Crossref] [PubMed]
- Rojas-Osornio SA, Cruz-Hernández TR, Drago-Serrano ME, et al. Immunity to influenza: Impact of obesity. Obes Res Clin Pract 2019;13:419-29. [Crossref] [PubMed]
- Cai G, Bossé Y, Xiao F, et al. Tobacco Smoking Increases the Lung Gene Expression of ACE2, the Receptor of SARS-CoV-2. Am J Respir Crit Care Med 2020;201:1557-9. [Crossref] [PubMed]
- Khorrami Z, Nili S, Sharifi H, et al. Association of cigarette smoking, obesity, and underlying medical conditions with COVID-19 hospitalization and mortality in Iran: A nationwide retrospective ecological study. Med J Islam Repub Iran 2020;34:133. [Crossref] [PubMed]
- Liu W, Tao ZW, Wang L, et al. Analysis of factors associated with disease outcomes in hospitalized patients with 2019 novel coronavirus disease. Chin Med J (Engl) 2020;133:1032-8. [Crossref] [PubMed]
- Lippi G, Sanchis-Gomar F, Henry BM. Active smoking and COVID-19: a double-edged sword. Eur J Intern Med 2020;77:123-4. [Crossref] [PubMed]
- Wenzl T. Smoking and COVID-19 - A review of studies suggesting a protective effect of smoking against COVID-19, EUR 30373 EN, Publications Office of the European Union, Luxembourg, 2020, ISBN 978-92-76-22062-6 (online), doi:
10.2760/564217 (online), JRC121837.10.2760/564217 - Ahmed J, Rizwan T, Malik F, et al. COVID-19 and Liver Injury: A Systematic Review and Meta-Analysis. Cureus 2020;12:e9424. [Crossref] [PubMed]
- Aziz M, Fatima R, Lee-Smith W, et al. The association of low serum albumin level with severe COVID-19: a systematic review and meta-analysis. Crit Care 2020;24:255. [Crossref] [PubMed]
- Chen G, Wu D, Guo W, et al. Clinical and immunological features of severe and moderate coronavirus disease 2019. J Clin Invest 2020;130:2620-9. [Crossref] [PubMed]
- Lee SW, Ha EK, Yeniova AÖ, et al. Severe clinical outcomes of COVID-19 associated with proton pump inhibitors: a nationwide cohort study with propensity score matching. Gut 2021;70:76-84. [Crossref] [PubMed]
- Hariyanto TI, Prasetya IB, Kurniawan A. Proton pump inhibitor use is associated with increased risk of severity and mortality from coronavirus disease 2019 (COVID-19) infection. Dig Liver Dis 2020;52:1410-2. [Crossref] [PubMed]
- Corsonello A, Lattanzio F, Bustacchini S, et al. Adverse Events of Proton Pump Inhibitors: Potential Mechanisms. Curr Drug Metab 2018;19:142-54. [Crossref] [PubMed]
- Namazi MR, Jowkar F. A succinct review of the general and immunological pharmacologic effects of proton pump inhibitors. J Clin Pharm Ther 2008;33:215-7. [Crossref] [PubMed]
- Guan WJ, Ni ZY, Hu Y, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med 2020;382:1708-20. [Crossref] [PubMed]
- Lei F, Liu YM, Zhou F, et al. Longitudinal Association Between Markers of Liver Injury and Mortality in COVID-19 in China. Hepatology 2020;72:389-98. [Crossref] [PubMed]
- Weber S, Hellmuth JC, Scherer C, et al. Liver function test abnormalities at hospital admission are associated with severe course of SARS-CoV-2 infection: a prospective cohort study. Gut 2021;70:1925-32. [Crossref] [PubMed]
- Wijarnpreecha K, Ungprasert P, Panjawatanan P, et al. COVID-19 and liver injury: a meta-analysis. Eur J Gastroenterol Hepatol 2021;33:990-5. [Crossref] [PubMed]
Cite this article as: Chela HK, Tallon EM, Baskett W, Gangu K, Tahan V, Shyu CR, Daglilar E. Liver injury on admission linked to worse outcomes in COVID-19: an analysis of 14,138 patients. Transl Gastroenterol Hepatol 2023;8:4.