
==== Front
Clin Appl Thromb Hemost
Clin Appl Thromb Hemost
CAT
spcat
Clinical and Applied Thrombosis/Hemostasis
1076-0296
1938-2723
SAGE Publications Sage CA: Los Angeles, CA

39215507
10.1177/10760296241275138
10.1177_10760296241275138
Original Research Article
Serum Bilirubin Levels and Risk of Venous Thromboembolism among Influenza Patients: A Cohort Study
Sutton S. Scott PharmD 1
https://orcid.org/0000-0002-8568-0097
Magagnoli Joseph PhD 1
Cummings Tammy PhD 1
Hardin James W. PhD 2
1 Dorn Research Institute, Columbia VA Health Care System, Columbia, South Carolina; Department of Clinical Pharmacy and Outcomes Sciences, 15525 College of Pharmacy, University of South Carolina , Columbia, SC, USA
2 Dorn Research Institute, Columbia VA Health Care System, Columbia, South Carolina; Department of Epidemiology & Biostatistics, University of South Carolina, Columbia, SC, USA
Joseph Magagnoli, Department of Clinical Pharmacy and Outcomes Sciences, University of South Carolina College of Pharmacy, 715 Sumter Street 311B, Columbia, South Carolina 29208. Email: magagnol@mailbox.sc.edu
30 8 2024
Jan-Dec 2024
30 1076029624127513820 5 2024
17 7 2024
25 7 2024
© The Author(s) 2024
2024
SAGE Publications Inc. unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by-nc/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License (https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access page (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Objective

This study aimed to investigate the associations between total serum bilirubin levels and the incidence of venous thromboembolism (VTE) among patients with influenza infection.

Methods

A retrospective cohort study was conducted among outpatients with laboratory-confirmed influenza using data from the Veterans Affairs Informatics and Computing Infrastructure (VINCI). Propensity score weighting was applied to balance study groups across baseline covariates. Cox proportional hazards models assessed VTE risk by total bilirubin levels, adjusting for important covariates including age, sex, race, comorbidity index, BMI, and smoking status.

Results

A total of 487 patients with total bilirubin levels <0.3 mg/dL, 8608 patients with levels between 0.3–1 mg/dL, and 1148 patients with levels >1 mg/dL were included. Patients with bilirubin <0.3 mg/dL exhibited a 6-fold higher risk of VTE compared to those with levels 0.3–1 mg/dL within 30 days of infection (HR = 6.2, 95% CI = 1.46–26.42). Elevated risks were noted through 90 days post infection (HR = 4.71, 95% CI = (1.42–15.67))

Conclusions

Serum bilirubin levels, particularly below 0.3 mg/dL, were significantly associated with an increased risk of VTE among individuals with influenza. These findings suggest that lower bilirubin levels may contribute to heightened inflammatory responses and subsequent thromboembolic events in patients with influenza. The underlying mechanisms and potential therapeutic implications for VTE prevention among patients with acute respiratory infection warrants further consideration.

bilirubin
inflammation
influenza
viral infection
vte
venous
thrombosis
cohort
typesetterts19
cover-dateJanuary-December 2024
==== Body
pmcIntroduction

Bilirubin is a byproduct of heme catabolism and found to be inversely correlated with the incidence of thromboembolic and cardiovascular events.1–3 This relationship suggests a potential protective role for bilirubin, attributed to its antioxidant properties. Conversely, influenza virus infection is associated with an increased risk of venous thromboembolism (VTE), encompassing both deep vein thrombosis (DVT) and pulmonary embolism (PE).4–8

The mechanisms underlying influenza-associated thrombosis are multifaceted, with pro-oxidant and inflammatory responses playing pivotal roles. Influenza infection induces excessive production of reactive oxygen species (ROS),9–12 creating a pro-oxidant and inflammatory environment favorable to coagulation.13–15 Influenza triggers a cascade of proinflammatory cytokines, including IL-1β, TNF-α, and IL-6.16,17 These inflammatory mediators contribute to endothelial dysfunction and activate the coagulation pathway through tissue factor (TF) induction, creating a prothrombotic state.

The interaction between bilirubin levels and VTE after influenza infection is currently unexplored. We hypothesize that lower serum bilirubin levels during influenza infection may exacerbate the pro-oxidant state, potentially amplifying the risk of subsequent VTE. Our study aims to investigate the association between total serum bilirubin levels and the incidence of thromboembolic events in a national cohort of patients with influenza.

Materials and Methods

Data Source

This retrospective study examined VTE among outpatients with laboratory confirmed influenza. Data on patient demographics, medical claims history, hospitalization, and outpatient medication dispensation was obtained using the Veterans Affairs Informatics and Computing Infrastructure (VINCI). The study received Dorn Research Institute Institutional Review Board (IRB No.1139248) and VA Research and Development approval and was conducted in compliance with the Department of Veterans Affairs requirements. This paper used STROBE guidelines for the reporting of observational studies.

Cohort Creation

Patients were included in the study if they had a positive influenza polymerase chain reaction (PCR) or antigen test along with an outpatient diagnosis code indicating fever, cough, acute respiratory infection, or influenza. Encounters for influenza with a prior positive influenza lab result within 180 days were excluded. Patients were required to have a bilirubin laboratory taken on the influenza diagnosis date and were excluded if they had no bilirubin lab result, were hospitalized on the day of diagnosis or had a hospital admission or discharge date within 30 days prior to diagnosis. Patients with any outpatient dispense of an anticoagulant medication including warfarin, apixaban, rivaroxaban, edoxaban, dabigatran, or enoxaparin were excluded. Patients with a missing body mass index (BMI) were also excluded.

Study Outcome

Venous thromboembolism was categorized similar to prior work, 18 utilizing a combination of diagnosis codes and anticoagulant medication dispenses. VTE was indicated if patients had an inpatient diagnosis code indicating VTE, codes listed in Supplementary Table S1, or at least 1 outpatient VTE diagnosis code and an initial anticoagulant medication within 10 days after the outpatient diagnosis. Anticoagulant medications include an outpatient dispense record for warfarin, apixaban, rivaroxaban, edoxaban, dabigatran, or enoxaparin, a common procedure terminology (CPT) code indicating administration for an anticoagulant (J1642, J1643, J1644, J1650, J1652) or an IV infusion of heparin or fondaparinux. VTE incidence was evaluated until a maximum of 180 days post influenza diagnosis. Patients with records indicating a VTE within 182 days prior to influenza were excluded.

Statistical Analysis

Patients were categorized into mutually exclusive cohorts based on their total bilirubin levels on the diagnosis date. Cohorts were broken down into three mutually exclusive groups based on their bilirubin levels including less than 0.3 mg/dL, between 0.3 and 1 mg/dL and greater than 1 mg/dL. Baseline characteristics were summarized using means, proportions, and standard deviations. Statistical tests, including T-tests and Chi-square tests, assessed cohort differences. Standardized differences (Std.Diff) were calculated to quantify the degree of similarity between the groups, with values less than 0.2 or ideally, 0.1 indicating good balance between cohorts. To analyze the risk of VTE among the cohort groups the cumulative hazard of VTE was assessed through 180 days stratified by the bilirubin level. Cumulative event counts and the number of patients at risk are also presented.

To minimize the possibility of confounding we used propensity score weighting. The propensity score model was specified using a multinomial logistic regression model that includes all baseline covariates. Weights we created as 1 divided by the estimated propensity score for respective treatment received. Weights were truncated at the 99th percentile. Means and percentages were compared with standardized differences across weighted samples. Subsequently, weighted Cox proportional hazards models evaluated hazards of VTE among influenza patients with varying bilirubin levels, adjusting for covariates including age, sex, race, comorbidity index, BMI, smoking status, and year of influenza diagnosis. Models were fit at different follow-times including withing 30, 60,90, and 180 days. A weighted Cox model was fit to examine differences in mortality within 90 days of diagnosis.

2020 and Later Subgroup Analysis

Variation among strains of influenza over time along with other respiratory viruses such as COVID-19 could play a role in the development of VTE. We conducted a sub-analysis among patients with influenza diagnosed in 2020 and later. Among this subgroup, propensity score weighting was conducted by first estimating a multinomial logistic regression to create the propensity score specific for the subgroup. Inverse propensity score weights were then calculated by one divided by the probability of treatment received. Weights were then truncated at the 99th percentile. Weighted Cox models were fit to estimate hazards of VTE.

Results

Sample Characteristics

Table 1 displays the the baseline characterstics for the influenza patients included in the sample. Among these patients a total of 487 patients had bilirubin levels <0.3 mg/dL, 8608 patients had bilirubin levels between 0.3–1 mg/dL, and 1148 patients had bilirubin levels >1 mg/dL were included. Mean bilirubin levels for each cohort were as follows: 0.2 in the <0.3 mg/dL cohort, 0.58 in the 0.3–1 mg/dL cohort, and 1.47 in the >1 mg/dL cohort. Antiviral treatment usage did not exhibit a significant difference across bilirubin groups (p = 0.45, Std.Diff =0.038). Age was positively correlated across the groups (p < 0.001, Std.Diff = 0.348), with older age being associated with higher bilirubin levels. Race distribution did not significantly differ among the groups (p = 0.1, Std.Diff = 0.124), nor was there significant variation among BMI (p = 0.263, Std.Diff = 0.119), drug or alcohol abuse (p = 0.07, Std.Diff = 0.113). Higher proportions of smokers were reported in the lower bilirubin cohorts (p = 0.047, Std.Diff = 0.11). The Charlson comorbidity index exhibited significant differences (p = 0.011, Std.Diff = 0.122), reflecting higher average comorbidity burdens as bilirubin increased. Table 2 illustrates propensity-weighted characteristics of the cohort groups. Propensity weighting achieved a balanced distribution for covariates across the study cohorts, as reflected by low (all less than 0.1) standardized differences.

Table 1. Sample Characteristics.

	Variable	<0.3 mg/dL	0.3–1 mg/dL	>1 mg/dL	P.value	Std.Diff*	
N = 478	N = 8608	N = 1148	
						
Bilirubin (mg/dL)	0.2(0.04)	0.58(0.2)	1.47(0.64)	<0.001	2.782	
Antiviral treatment		310(64.8%)	5433(63.12%)	743(64.72%)	0.45	0.038	
Age		53.76(17.5)	57.65(17.69)	59.98(18.24)	<0.001	0.348	
Race	BLACK	122(25.523%)	1989(23.11%)	233(20.3%)	0.1	0.124	
	OTHER/UNKNOWN	37(7.741%)	816(9.48%)	116(10.11%)	0.1	0.124	
	WHITE	319(66.736%)	5803(67.41%)	799(69.6%)	0.1	0.124	
Sex	FEMALE	163(34.1%)	1027(11.93%)	47(4.09%)	<0.001	0.826	
	MALE	315(65.9%)	7581(88.07%)	1101(95.91%)	<0.001	0.826	
Charlson comorbidity index		1.24(2.2)	1.34(1.93)	1.51(2.14)	0.011	0.122	
BMI	<18.5	7(1.46%)	89(1.03%)	11(0.96%)	0.263	0.119	
	18.5–24.9	104(21.76%)	1497(17.39%)	196(17.07%)	0.263	0.119	
	25–29.9	147(30.75%)	2826(32.83%)	388(33.8%)	0.263	0.119	
	30+	220(46.03%)	4196(48.75%)	553(48.17%)	0.263	0.119	
Drug or alcohol abuse		33(6.90%)	411(4.78%)	49(4.27%)	0.07	0.113	
Current smoker		32(6.7%)	502(5.83%)	48(4.18%)	0.047	0.11	
Year of influenza diagnosis		2019.38(2.92)	2019.34(3.0)	2019.24(3.2)	0.572	0.046	
.*Maximum of pairwise standardized differences

Table 2. Propensity Weighted Characteristics.

		<0.3 mg/dL	0.3–1 mg/dL	>1 mg/dL	Std. Diff*	
Antiviral Treatment	64.10%	63.40%	63.90%	0.016	
Age	58.147	57.726	57.574	0.031	
Race	BLACK	24.40%	22.90%	22.80%	0.038	
	OTHER/UNKNOWN	9.10%	9.50%	9.20%	0.014	
	WHITE	66.50%	67.60%	68.00%	0.031	
Sex	FEMALE	12.10%	12.10%	11.90%	0.007	
	MALE	87.90%	87.90%	88.10%	0.007	
Charlson comorbidity index	1.436	1.35	1.375	0.044	
BMI	<18.5	1.30%	1.00%	1.10%	0.024	
	18.5–24.9	21.10%	17.40%	17.00%	0.104	
	25–29.9	32.50%	32.80%	32.80%	0.008	
	30+	45.20%	48.70%	49.00%	0.077	
Drug or alcohol abuse	5.20%	4.80%	4.60%	0.03	
Current smoker	6.20%	5.70%	5.50%	0.031	
Year of influenza diagnosis	2019.419	2019.323	2019.391	0.032	
.*Maximum of pairwise standardized differences

Bilirubin and the Risk of VTE

Figure 1 presents cumulative incidence of VTE across the study period among the original study data. Across the study period, patients with low bilirubin levels (<0.3 mg/dL) and those with bilirubin greater than 1 mg/dL have higher cumulative VTE hazards compared to the 0.3–1 mg/dL cohort. The cohort with the lowest bilirubin levels (<0.3 mg/dL) had the highest hazard of VTE over the study period.

Figure 1. Kaplan-Meier cumulative incidence curves.

Table 3 presents the results of propensity weighted Cox proportional hazards models evaluating the hazard of VTE associated with bilirubin levels across time points. Hazard ratios (HR) with 95% confidence intervals are reported. Patients with bilirubin levels lower than 0.3 mg/dL exhibited a six-fold increased hazard of VTE at 30 days (HR = 6.2, 95% CI = 1.46–26.42). At subsequent follow-up time points the magnitude of the HR reduced: 60 days (HR = 5.34, 95% CI = 1.59–18.01) and 90 days (HR = 4.71, 95% CI = 1.42 −15.67), 180 days (HR = 2.8, 95% CI = 0.87–8.99) compared to those with levels between 0.3 and 1.0 mg/dL. Patients with bilirubin levels greater than 1 mg/dL did not have a statistically significantly greater risk of VTE in any time period compared to patients with bilirubin levels between 0.3–1 mg/dL. However, while not statistically significant, point estimates for the patients with >1 mg/dL all indicated higher risk compared to patients with bilirubin levels between 0.3–1 mg/dL. Moreover, the HRs were more consistent over the time frame ranging from a low of 1.81 (180 days) to 2.6 (60 days). This is a noticeable difference compared to the hazard ratio estimated for the patients with <0.3 mg/dL which peaked at the earliest time point (30 days) and descended with all subsequent time points.

Table 3. Propensity Score Weighted Cox Models.

	30 days	60 days	90 days	180 days	90 day mortality	
Variable	HR(95% CI)	HR(95% CI)	HR(95% CI)	HR(95% CI)	HR(95% CI)	
Antiviral Treatment	0.98(0.34–2.85)	0.93(0.39–2.2)	0.87(0.38–1.99)	0.71(0.37–1.37)	0.5(0.28–0.88)	
Bilirubin (ref 0.3–1 mg/dL)						
 < 0.3 mg/dL	6.2(1.46–26.42)	5.34(1.59–18.01)	4.71(1.42–15.67)	2.8(0.87–8.99)	1.19(0.26–5.42)	
 > 1 mg/dL	2.31(0.63–8.49)	2.6(0.93–7.26)	2.25(0.82–6.2)	1.81(0.78–4.18)	2.19(1.06–4.52)	
Age	1.02(1.01–1.03)	1.02(1.01–1.03)	1.02(1.01–1.03)	1.01(1–1.02)	1.03(1.02–1.04)	
Race Other/unknown versus Black	0(0–0)	0(0–0)	0(0–0)	0(0–0)	1.59(0.5–5.05)	
White versus Black	0.86(0.28–2.64)	0.69(0.29–1.68)	0.76(0.32–1.79)	0.85(0.41–1.74)	1.62(0.7–3.76)	
Male	5.08(0.69–37.47)	7.06(0.96–51.7)	7.36(1.01–53.64)	2.86(0.59–13.85)	1.23(0.17–9.19)	
Charlson comorbidity index	1.01(0.86–1.19)	1.1(0.96–1.27)	1.13(0.99–1.3)	1.18(1.07–1.31)	1.32(1.23–1.41)	
Current smoker	1.57(0.21–11.75)	1.07(0.14–8.02)	1(0.14–7.48)	1.78(0.54–5.89)	0(0–0)	
Drug or alcohol abuse	0(0–0)	0.69(0.09–5.53)	0.53(0.06–4.54)	0.32(0.04–2.75)	0.84(0.2–3.64)	
Year of influenza diagnosis	0.97(0.86–1.11)	0.98(0.87–1.1)	0.92(0.81–1.04)	0.92(0.83–1.02)	0.93(0.86–1)	

Bilirubin and the Risk of VTE: Years 2020 Onward

Table 4 presents the results of Propensity Score Weighted Cox models for the period January 2020 and later, evaluating the impact of bilirubin on risk of VTE across the study period (30 days, 60 days, and 180 days). Bilirubin levels lower than 0.3 mg/dL demonstrate a significantly increased risk of VTE compared to bilirubin levels of 0.3 to 1.0 mg/dl. Moreover, similar to the previous result, hazards ratios were highest at the 30-day time point and decline at subsequent time points. Interestingly, patients with bilirubin greater than 1 mg/dL also have a higher risk compared to those with 0.3–1 mg/dL at 30-days and hazard ratio magnitudes taper by 180 days.

Table 4. Propensity Score Weighted Cox Models: January 2020 and Later.

	30 days	60 days	180 days	
Variable	HR(95% CI)	HR(95% CI)	HR(95% CI)	
Antiviral Treatment	0.7(0.17–2.85)	0.73(0.24–2.22)	0.86(0.35–2.12)	
Bilirubin (ref 0.3–1 mg/dL)				
 < 0.3 mg/dL	6.6(0.89–48.95)	5.24(1.08–25.54)	2.65(0.60–11.55)	
 > 1 mg/dL	5.39(1.21–24.05)	4.08(1.19–13.99)	2.48(0.88–7.01)	
Age	1.03(1.01–1.04)	1.02(1.01–1.04)	1.01(1–1.03)	
Race Other/unknown versus Black	0(0–0)	0(0–0)	0(0–0)	
White versus Black	0.63(0.15–2.64)	0.64(0.2–2.08)	0.63(0.26–1.55)	
Male	2.9(0.38–21.82)	3.83(0.51–28.65)	6.35(0.87–46.47)	
Charlson comorbidity index	1.05(0.88–1.25)	1.13(0.95–1.36)	1.24(1.09–1.4)	
Current smoker	4.05(0.5–32.75)	2.77(0.38–20.21)	4.1(1.2–14.06)	
Drug or alcohol abuse	0(0–0)	0.98(0.12–7.68)	0.6(0.08–4.66)	
Year of influenza diagnosis	0.79(0.57–1.11)	0.79(0.61–1.03)	0.77(0.6–0.99)	

Discussion

Influenza virus infection imposes a considerable public health burden, affecting healthcare systems and causing individual morbidity and mortality. Recent evidence suggests a potential connection between influenza and an elevated risk of venous thromboembolism (VTE), including DVT and pulmonary embolism (PE). Influenza virus infection is associated with a significant overproduction of ROS,9–12 contributing to viral pathogenesis. This increase in ROS contributes to influenza pathogenicity through multiple mechanisms, including activation of innate immunity and promotion of inflammation.10,11 Recent studies have demonstrated that the influenza viral matrix 1 protein exacerbates viral pathogenicity by inducing TLR4-mediated ROS production and apoptotic cell death. 10 Furthermore, mitochondrial ROS have been implicated in the pathological inflammation observed during influenza A virus infection in mice. 12

The elevated levels of ROS during influenza infection can promote a procoagulant state, which is conducive to thrombus formation. ROS has been shown to increase the expression of tissue factor (TF), a key initiator of the coagulation cascade.13,14 Additionally, ROS-mediated NADPH oxidase activation in vascular smooth muscle cells can lead to tissue factor-dependent procoagulant activity. 15

Beyond promoting a procoagulant state, ROS can also impair anticoagulant mechanisms. Oxidative stress has been demonstrated to inactivate activated protein C, an important anticoagulant protein. 19 Moreover, oxidation of thrombomodulin, another crucial anticoagulant factor, can inhibit its cofactor activity, potentially providing a rapid mechanism for modulating coagulation. 20

While the link between excessive ROS and procoagulant states during influenza infection is well-established, recent studies have also investigated endogenous factors, such as bilirubin, that could also have a role in VTE formation. Previous studies reported an inverse relationship between bilirubin levels and thromboembolic and cardiovascular events.1–3

Bilirubin is a metabolite of heme that scavenges ROS preventing oxidation. Bilirubin is oxidized to biliverdin by taking on oxygen free radicals and then subsequent conversion back to bilirubin through NADPH/biliverdin reductase.21–25 Bilirubin's ability to prevent oxidative stress is thought to be linked to its structure, with the ability of its exposed hydrogen molecule to combine with oxygen radicals.21,26,27 Moreover, bilirubin could act as a regulator of inflammation through its action on vascular cell adhesion molecule-1(VCAM-1). 28 Therefore, it is possible that, in the context of viral infections, patients with low bilirubin levels may be more at risk for complications such as VTE formation that are in part due to oxidative stress and inflammation.

Using a national cohort of patients with influenza and bilirubin lab results, we found that patients with total serum bilirubin levels below 0.3 mg/dL exhibited a significantly higher hazard of VTE following influenza infection compared to those with higher bilirubin levels. Moreover, we found that lower serum bilirubin levels were more significant predictors of VTE during follow-up periods closer to the influenza infection point suggesting that low bilirubin levels are particularly important in the context of acute viral or inflammatory illness.

Given the possibility of acute inflammation exacerbating VTE risk, we investigated if two medications, montelukast or metformin could alter the subsequent risk of VTE among our sample of influenza patients. Montelukast, a commonly used leukotriene receptor antagonist, has been previously identified as being protective against VTE 29 and preclinical evidence purports its benefit in preventing thrombus formation among broad cardiovascular benefits. 30 Prior work found metformin associated with lower mortality among obese patients with influenza. 31 Metformin has a growing body of research examining its anti-inflammatory effects,32,33 which could be beneficial in VTE prevention. For each medication patients were classified as either exposed or unexposed on their influenza infection date and Cox models were fit for each medication separately, estimating the hazard ratio (HR) at each follow up period. Neither montelukast nor metformin demonstrated any association with subsequent VTE among influenza patients (Supplementary Table S2).

Subgroup analysis of more recent influenza infections (year 2020+) revealed a consistent association between lower bilirubin levels and VTE risk. Different strains of influenza could result in different risks for thromboembolic events. For example, H1N1 influenza during the 2009 influenza season was particularly associated with thrombotic events. 6 The results of the current study, while not parsing out specific influenza virus variants, indicate that lower bilirubin levels remained an important predictor of VTE development among influenza patients in the post-2020 era.

Point estimates for bilirubin levels >1 mg/dL consistently pointed to an increased risk compared to 0.3–1 mg/dL levels, however in the main analysis the finding was not statistically significant. Interestingly in the post-2020 era, bilirubin levels greater than 1 mg/dL were found to increase VTE risk. This result, along with the consistent point estimates in the main analysis, suggests a possible U-shape risk profile with patients having low or high bilirubin levels conferring higher VTE risk. Our results suggesting a U-shape risk profile between bilirubin values and VTE is consistent with other results examining bilirubin and cardiovascular disease risk.34–38 While low bilirubin levels may allow an environment favorable to coagulation, excessively high bilirubin levels might be an indicator of liver damage, which can increase risk of VTE.39–43 Moreover, preclinical data suggest that increased levels of bilirubin could result in a switch from an antioxidant to a pro-oxidant effect. 44

Limitations of this study warrant consideration. Being a non-randomized study the potential for confounding exists. To minimize the confounding bias, we used propensity score weighting to make the cohorts as similar as possible. However, despite inclusion of important covariates, there may be other factors such as genetics, lifestyle factors and differences in the influenza strains not captured in the medical records that could impact the risk of VTE. Moreover, the incidence of other respiratory illnesses such as COVID-19 could impact the risk of VTE and this was not accounted for in our analysis.

The results of this study demonstrated a significant association between low total serum bilirubin levels and increased VTE risk in patients with influenza infection. Importantly, sub-analysis also suggests that increased total bilirubin levels may be associated with increased VTE risk. These results are consistent with previously observed U-shape risk profile found for bilirubin levels in cardiovascular disease. Further investigations to examine bilirubin and VTE risk among broader populations are warranted.

Supplemental Material

sj-docx-1-cat-10.1177_10760296241275138 - Supplemental material for Serum Bilirubin Levels and Risk of Venous Thromboembolism among Influenza Patients: A Cohort Study

Supplemental material, sj-docx-1-cat-10.1177_10760296241275138 for Serum Bilirubin Levels and Risk of Venous Thromboembolism among Influenza Patients: A Cohort Study by S. Scott Sutton, PharmD, Joseph Magagnoli, PhD, Tammy Cummings, PhD, and James W. Hardin, PhD in Clinical and Applied Thrombosis/Hemostasis

Additional Contributions: We thank Charles L. Bennett, MD, PhD, MPP, Professor, College of Pharmacy, University of South Carolina, for discussions and editing of the manuscript. Dr Bennett did not receive financial compensation for this contribution.

Author Disclosures: Sutton has received research grants from Boehringer Ingelheim, Gilead Sciences, Alexion Pharmaceuticals, and EMD Serono, all for projects unrelated to study. The other authors declare no competing interests. The content of this article is solely the responsibility of the authors and does not necessarily represent the official views of the US Department of Veterans Affairs, nor does mention of trade names, commercial products or organizations imply endorsement by the US government. This paper represents, in part, original research conducted using data from the Department of Veterans Affairs and is the result of work supported with resources and the use of facilities at the Dorn Research Institute, Columbia VA Health Care System, Columbia, South Carolina.

Analyses of the Veterans Health Administration Database were performed using data within the US Department of Veterans Affairs secure research environment, the VA Informatics and Computing Infrastructure (VINCI). The completeness, utility, accuracy, validity, and access methods are described on the VA website, https://www.virec.research.va.gov.

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding: Sutton, Magagnoli, and Cummings are supported by NIH grant R01DA054992 and the South Carolina Center for Rural and Primary Healthcare unrelated to this study.

ORCID iD: Joseph Magagnoli https://orcid.org/0000-0002-8568-0097

Supplemental Material: Supplemental material for this article is available online.
==== Refs
References

1 Duman H Özyurt S Erdoğan T Kara BY Durakoğlugil ME . The role of serum bilirubin levels in determining venous thromboembolism. J Vasc Surg Venous Lymphat Disord. 2019;7 (5 ):635-639. doi:10.1016/j.jvsv.2019.02.002 30922986
2 Marconi VC Duncan MS So-Armah K , et al. Bilirubin is inversely associated with cardiovascular disease among HIV-positive and HIV-negative individuals in VACS (veterans aging cohort study). J Am Heart Assoc. 2018;7 (10 ):e007792. Published 2018 May 2. doi:10.1161/JAHA.117.007792
3 Suh S Cho YR Park MK Kim DK Cho NH Lee MK . Relationship between serum bilirubin levels and cardiovascular disease. PLoS One. 2018;13 (2 ):e0193041. Published 2018 Feb 15. doi:10.1371/journal.pone.0193041.
4 Silverstein MD Heit JA Mohr DN Petterson TM O'Fallon WM Melton LJ . Trends in the incidence of deep vein thrombosis and pulmonary embolism: A 25-year population-based study. Arch Intern Med. 1998;158 (6 ):585-593. doi:10.1001/archinte.158.6.585 9521222
5 van Wissen M Keller TT Ronkes B , et al. Influenza infection and risk of acute pulmonary embolism. Thromb J. 2007;5 (1 ):16. Published 2007 Oct 16. doi:10.1186/1477-9560-5-16 17939867
6 Bunce PE High SM Nadjafi M Stanley K Liles WC Christian MD . Pandemic H1N1 influenza infection and vascular thrombosis. Clin Infect Dis. 2011;52 (2 ):e14-e17. doi:10.1093/cid/ciq125
7 Rubino R Imburgia C Bonura S Trizzino M Iaria C Cascio A . Thromboembolic events in patients with influenza: A scoping review. Viruses. 2022;14 (12 ):2817. Published 2022 Dec 17. doi:10.3390/v14122817 36560821
8 Smeeth L Cook C Thomas S Hall AJ Hubbard R Vallance P . Risk of deep vein thrombosis and pulmonary embolism after acute infection in a community setting. Lancet. 2006;367 (9516 ):1075-1079. doi:10.1016/S0140-6736(06)68474-2 16581406
9 Amatore D Sgarbanti R Aquilano K , et al. Influenza virus replication in lung epithelial cells depends on redox-sensitive pathways activated by NOX4-derived ROS. Cell Microbiol. 2015;17 (1 ):131-145. doi:10.1111/cmi.12343 25154738
10 Kim CU Lim D Kim YS , et al. Influenza viral matrix 1 protein aggravates viral pathogenicity by inducing TLR4-mediated reactive oxygen species production and apoptotic cell death. Cell Death Dis. 2023;14 (1 ):228. doi:10.1038/s41419-023-05749 36990977
11 Wang L Cao Z Wang Z Guo J Wen J . Reactive oxygen species associated immunoregulation post influenza virus infection. Front Immunol. 2022;13 :927593. Published 2022 Jul 29. doi:10.3389/fimmu.2022.927593 35967412
12 To EE Erlich JR Liong F , et al. Mitochondrial reactive oxygen species contribute to pathological inflammation during influenza A virus infection in mice. Antioxid Redox Signal. 2020;32 (13 ):929-942. doi:10.1089/ars.2019.7727 31190565
13 Golino P Ragni M Cirillo P , et al. Effects of tissue factor induced by oxygen free radicals on coronary flow during reperfusion. Nat Med. 1996;2 (1 ):35-40. doi:10.1038/nm0196-35 8564835
14 Cadroy Y Dupouy D Boneu B Plaisancié H . Polymorphonuclear leukocytes modulate tissue factor production by mononuclear cells: Role of reactive oxygen species. J Immunol. 2000;164 (7 ):3822-3828. doi:10.4049/jimmunol.164.7.3822 10725743
15 Herkert O Diebold I Brandes RP Hess J Busse R Görlach A . NADPH Oxidase mediates tissue factor-dependent surface procoagulant activity by thrombin in human vascular smooth muscle cells. Circulation. 2002;105 (17 ):2030-2036. doi:10.1161/01.cir.0000014611.28864.1e 11980681
16 Yang Y Tang H . Aberrant coagulation causes a hyper-inflammatory response in severe influenza pneumonia. Cell Mol Immunol. 2016;13 (4 ):432-442. doi:10.1038/cmi.2016.1 27041635
17 Branchford BR Carpenter SL . The role of inflammation in venous thromboembolism. Front Pediatr. 2018;6 :142. Published 2018 May 23. doi:10.3389/fped.2018.00142 29876337
18 Li A da Costa WL Jr Guffey D , et al. Developing and optimizing a computable phenotype for incident venous thromboembolism in a longitudinal cohort of patients with cancer. Res Pract Thromb Haemost. 2022;6 (4 ):e12733. Published 2022 May 25. doi:10.1002/rth2.12733
19 Nalian A Iakhiaev AV . Possible mechanisms contributing to oxidative inactivation of activated protein C: Molecular dynamics study. Thromb Haemost. 2008;100 (1 ):18-25. doi:10.1160/TH08-02-0124 18612533
20 Glaser CB Morser J Clarke JH , et al. Oxidation of a specific methionine in thrombomodulin by activated neutrophil products blocks cofactor activity. A potential rapid mechanism for modulation of coagulation. J Clin Invest. 1992;90 (6 ):2565-2573. doi:10.1172/JCI116151 1334978
21 Chen Z Vong CT Gao C , et al. Bilirubin nanomedicines for the treatment of reactive oxygen species (ROS)-mediated diseases. Mol Pharm. 2020;17 (7 ):2260-2274. doi:10.1021/acs.molpharmaceut.0c00337 32433886
22 Kundur AR Singh I Bulmer AC . Bilirubin, platelet activation and heart disease: A missing link to cardiovascular protection in Gilbert's syndrome? Atherosclerosis. 2015;239 (1 ):73-84. doi:10.1016/j.atherosclerosis.2014.12.042 25576848
23 Minetti M Mallozzi C Di Stasi AM Pietraforte D . Bilirubin is an effective antioxidant of peroxynitrite-mediated protein oxidation in human blood plasma. Arch Biochem Biophys. 1998;352 (2 ):165-174. doi:10.1006/abbi.1998.0584 9587403
24 McDonagh AF . The biliverdin-bilirubin antioxidant cycle of cellular protection: Missing a wheel? Free Radic Biol Med. 2010;49 (5 ):814-820. doi:10.1016/j.freeradbiomed.2010.06.001 20547221
25 Sedlak TW Saleh M Higginson DS Paul BD Juluri KR Snyder SH . Bilirubin and glutathione have complementary antioxidant and cytoprotective roles. Proc Natl Acad Sci U S A. 2009;106 (13 ):5171-5176. doi:10.1073/pnas.0813132106 19286972
26 Stocker R Yamamoto Y McDonagh AF Glazer AN Ames BN . Bilirubin is an antioxidant of possible physiological importance. Science. 1987;235 (4792 ):1043-1046. doi:10.1126/science.3029864 3029864
27 Asad SF Singh S Ahmad A Khan NU Hadi SM . Prooxidant and antioxidant activities of bilirubin and its metabolic precursor biliverdin: A structure–activity study. Chem Biol Interact. 2001;137 (1 ):59-74. doi:10.1016/S0009-2797(01)00209-5 11518564
28 Vogel ME Zucker SD . Bilirubin acts as an endogenous regulator of inflammation by disrupting adhesion molecule-mediated leukocyte migration. Inflamm Cell Signal. 2016;3 (1 ):e1178. doi:10.14800/ics.1178
29 Sutton SS Magagnoli J Cummings TH Hardin JW . Leukotriene receptor antagonism with montelukast as a possible therapeutic for venous thromboembolism prophylaxis: An observational study. Prostaglandins Other Lipid Mediat. 2022 Aug;161 :106649. doi:10.1016/j.prostaglandins.2022.106649. Epub 2022 May 18. PMID: 35595009.35595009
30 Alomair BM Al-Kuraishy HM Al-Gareeb AI , et al. Montelukast and acute coronary syndrome: The endowed drug. Pharmaceuticals (Basel). 2022;15 (9 ):1147. Published 2022 Sep 14. doi:10.3390/ph15091147 36145367
31 Cummings TH Magagnoli J Hardin JW Sutton SS . Patients with obesity and a history of metformin treatment have lower influenza mortality: A retrospective cohort study. Pathogens. 2022 Feb 19;11 (2 ):270. doi:10.3390/pathogens11020270. PMID: 35215211; PMCID: PMC8876732.35215211
32 Saisho Y . Metformin and inflammation: Its potential beyond glucose-lowering effect. Endocr Metab Immune Disord Drug Targets. 2015;15 (3 ):196-205. doi:10.2174/1871530315666150316124019 25772174
33 Wang Z Wang M Lin M Wei P . The immunomodulatory effects of metformin in LPS-induced macrophages: An in vitro study. Inflamm Res. 2024;73 (2 ):175-181. doi:10.1007/s00011-023-01827-8 38091014
34 Lai X Fang Q Yang L , et al. Direct, indirect and total bilirubin and risk of incident coronary heart disease in the Dongfeng-Tongji cohort. Ann Med. 2018;50 (1 ):16-25. doi:10.1080/07853890.2017.1377846 28879775
35 Troughton JA Woodside JV Young IS , et al. Bilirubin and coronary heart disease risk in the Prospective Epidemiological Study of Myocardial Infarction (PRIME). Eur J Cardiovasc Prev Rehabil. 2007;14 (1 ):79-84. doi:10.1097/01.hjr.0000230097.81202.9f 17301631
36 Breimer LH Wannamethee G Ebrahim S Shaper AG . Serum bilirubin and risk of ischemic heart disease in middle-aged British men. Clin Chem. 1995;41 (10 ):1504-1508. doi:10.1093/clinchem/41.10.1504 7586525
37 Gao Y Zhao LB Li K , et al. The J-shape association between total bilirubin and stroke in older patients with obstructive sleep apnea syndrome: A multicenter study. J Nutr Health Aging. 2023;27 (9 ):692-700. doi:10.1007/s12603-023-1965-2 37754208
38 Li C Wu W Song Y Xu S Wu X . The nonlinear relationship between total bilirubin and coronary heart disease: A dose-response meta-analysis. Front Cardiovasc Med. 2022;8 :761520. Published 2022 Jan 5. doi:10.3389/fcvm.2021.761520 35071344
39 Northup PG McMahon MM Ruhl AP , et al. Coagulopathy does not fully protect hospitalized cirrhosis patients from peripheral venous thromboembolism. Am J Gastroenterol. 2006;101 (7 ):1524-1528. doi:10.1111/j.1572-0241.2006.00588.x 16863556
40 Søgaard KK Horváth-Puhó E Grønbaek H Jepsen P Vilstrup H Sørensen HT . Risk of venous thromboembolism in patients with liver disease: A nationwide population-based case-control study. Am J Gastroenterol. 2009;104 (1 ):96-101. doi:10.1038/ajg.2008.34 19098856
41 Gulley D Teal E Suvannasankha A Chalasani N Liangpunsakul S . Deep vein thrombosis and pulmonary embolism in cirrhosis patients. Dig Dis Sci. 2008;53 (11 ):3012-3017. doi:10.1007/s10620-008-0265-3 18443906
42 Tsochatzis EA Senzolo M Germani G Gatt A Burroughs AK . Systematic review: Portal vein thrombosis in cirrhosis. Aliment Pharmacol Ther. 2010;31 (3 ):366-374. doi:10.1111/j.1365-2036.2009.04182.x 19863496
43 Ambrosino P Tarantino L Di Minno G , et al. The risk of venous thromboembolism in patients with cirrhosis. A systematic review and meta-analysis. Thromb Haemost. 2017;117 (1 ):139-148. doi:10.1160/TH16-06-0450 27761574
44 Bianco A Dvořák A Capková N , et al. The extent of intracellular accumulation of bilirubin determines its anti- or pro-oxidant effect. Int J Mol Sci. 2020;21 (21 ):8101. 10.3390/ijms21218101 33143041
