ABSTRACT
-
Background/Aims
We aimed to determine the association between blood urea level and incident cirrhosis, hepatic decompensation, and hepatocellular carcinoma in chronic liver disease (CLD) patients.
-
Methods
The association between blood urea level and liver fibrosis/liver-related events were evaluated on continuous scale with restricted cubic spline curves based on generalized additive model or Cox proportional hazards models. Then, the above associations were evaluated by urea level within intervals.
-
Results
Among 4,282 patients who had undergone liver stiffness measurement (LSM) by transient elastography, baseline urea level had a U-shaped association with LSM and hepatic decompensation development after a median follow-up of 5.5 years. Compared to patients with urea of 3.6–9.9 mmol/L, those with urea ≤3.5 mmol/L (adjusted hazard ratio [aHR], 4.15; 95% confidence interval [CI], 1.68–10.24) and ≥10 mmol/L (aHR, 5.22; 95% CI, 1.86–14.67) had higher risk of hepatic decompensation. Patients with urea ≤3.5 mmol/L also had higher risk of incident cirrhosis (aHR, 3.24; 95% CI, 1.50–6.98). The association between low urea level and incident cirrhosis and hepatic decompensation was consistently observed in subgroups by age, gender, albumin level, and comorbidities. The U-shaped relationship between urea level and LSM was validated in another population screening study (n=917). Likewise, urea ≤3.5 mmol/L was associated with a higher risk of incident cirrhosis in a territory-wide cohort of 12,476 patients with nonalcoholic fatty liver disease at a median follow-up of 9.9 years (aHR, 1.27; 95% CI, 1.03–1.57).
-
Conclusions
We identified a U-shaped relationship between the urea level and liver fibrosis/incident cirrhosis/hepatic decompensation in patients with CLD.
-
Keywords: Urea; Liver cirrhosis; Fibrosis; Non-alcoholic fatty liver disease; Hepatitis B
Study Highlights
Increased blood urea level indicates impaired renal function, and decreased urea level can be due to poor liver function and protein-energy malnutrition. The association between low blood urea level and liver fibrosis/liver-related events has not been evaluated in previous studies. There was a U-shaped relationship between the urea level and the severity of liver disease. Patients with low urea level had increased risk of liver fibrosis, incident cirrhosis and hepatic decompensation. Clinicians should watch out for advanced liver disease in patients with low urea level.
Graphical Abstract
INTRODUCTION
Chronic liver diseases (CLDs) represent a global health problem and account for approximately 2 million deaths per year worldwide [
1]. Liver fibrosis is the result of liver injury from various CLDs and is the common pathway towards cirrhosis and its complications. The level of fibrosis inversely correlates with liver function and the quality of life [
2].
Cirrhosis is the terminal stage of progressive liver fibrosis, and is associated with increased risk of hepatocellular carcinoma (HCC) and complications of portal hypertension [
3]. Cirrhosis is classified as compensated or decompensated disease based on the absence or presence of portal hypertension-related complications (ascites, spontaneous bacterial peritonitis, variceal bleeding), hepatic encephalopathy, and/or hepatorenal syndrome [
4]. Most morbidity and mortality in patients with liver cirrhosis can be attributed to the transition to the decompensated state. Patients with compensated cirrhosis are expected to survive more than 10 years, compared with a median survival of 2 years in patients with hepatic decompensation [
5,
6]. Therefore, it is important to identify CLD patients at high risk of developing cirrhosis and decompensation.
Urea is a non-toxic, nitrogenous organic end product of protein metabolism, facilitating 80–90% of nitrogen elimination from human body [
7]. Increased blood urea level indicates impaired renal function, and decreased urea level can be due to poor liver function and protein-energy malnutrition [
8]. Serum creatinine, another marker of renal function, is a component of the Model for Endstage Liver Disease score, which is used for prognostication and liver transplant prioritization in patients with various liver conditions [
4]. Accordingly, high urea level also correlates with poor liver outcomes. Recently, several experimental studies revealed that urea synthesis is impaired even in early CLD and this correlates with hepatic fibrosis [
9-
11]. Additionally, urea cycle dysregulation contributes to carcinogenesis [
12]. Therefore, CLD patients with low urea level may have more advanced disease and be at risk of disease progression.
With this background, we hypothesize that there is a U-shaped relationship between blood urea level and the severity of liver disease, with both low and high urea level indicative of at-risk patients. To test this hypothesis, we studied three cohorts and demonstrated consistent results in the association with liver stiffness measurement (LSM) and various liver-related events (LREs).
MATERIALS AND METHODS
Study populations
LSM cohort
This was a
post-hoc analysis of the LSM cohort study [
13]. Between August 2012 and March 2016 at the Prince of Wales Hospital, patients who underwent transient elastography (TE) examination for suspected or known CLDs were included in this study. Patients who had history of cancer, hepatic decompensation, unreliable LSM (defined below), or with missing data of blood urea level were excluded. In addition, patients with LSM >15 kPa (suspected of cirrhosis) or clinical diagnosis of cirrhosis before baseline were further excluded for the analysis of the association between urea level and incident cirrhosis.
Magnetic resonance spectroscopy (MRS) cohort
The association between urea level and LSM was validated in an independent cohort (MRS cohort). The MRS cohort was a population screening study including 917 subjects randomly recruited from the census database of the Hong Kong Government (May 2008 to September 2010) [
14]. All the subjects underwent TE and proton-magnetic resonance spectroscopy and did not have cancer or hepatic decompensation at the time of assessment.
Nonalcoholic fatty liver disease (NAFLD) cohort
The association between urea level and hepatic decompensation/incident cirrhosis/HCC was validated in another independent cohort (NAFLD cohort) using data from the Clinical Data Analysis and Reporting System (CDARS) of the Hospital Authority, Hong Kong. CDARS facilitates the retrieval of clinical data captured from different operational systems for analysis and reporting and provides good quality information to support retrospective clinical and management decisions by integrating the clinical data residing in data warehouse. It covers the electronic health records and laboratory results from all public hospitals and clinics in Hong Kong and represents data of approximately 80% of the local population. The NAFLD cohort included 12,476 patients with a diagnosis code of NAFLD (International Classification of Diseases, ninth revision, clinical modification [ICD-9-CM] code 571.8) between January 2000 and June 2020. Patients who had cancer or hepatic decompensation before baseline were excluded.
The study protocol was approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee. All subjects in the MRS cohort provided informed written consent, whereas the Ethics Committee exempted the need for consenting for the LSM and NAFLD cohorts because of their retrospective nature.
Clinical data
The baseline date was defined as the date when TE was performed (LSM and MRS cohorts) or the date when NAFLD was first coded (NAFLD cohort). Data collected at baseline included the complete blood count, prothrombin time and international normalized ratio, liver and renal biochemistries, fasting lipids, glucose, and glycated hemoglobin. Use of lipid lowering medications, antiviral medications, and anti-platelet agents were also recorded. All diagnoses and events were extracted from the CDARS based on ICD-9-CM codes, which were subsequently confirmed manually by two independent investigators using patient medical records (
Supplementary Table 1).
Chronic hepatitis B or C were diagnosed by positive serology tests for serum hepatitis B surface antigen or anti-hepatitis C virus antibodies, respectively. In LSM and MRS cohorts, NAFLD was diagnosed by the presence of hepatic steatosis seen on ultrasonography and/or histology or as evidenced by an elevated controlled attenuation parameter (CAP) ≥248 dB/m after exclusion of alternative liver diseases and secondary causes of fatty liver. Autoimmune and cholestatic liver diseases (primary biliary cholangitis and primary sclerosing cholangitis) were diagnosed based on standard serological, cholangiographic, and histologic criteria.
The primary outcome was new onset hepatic decompensation. Secondary outcomes include incident cirrhosis, and HCC. Hepatic decompensation was defined as new-onset ascites, spontaneous bacterial peritonitis, variceal bleeding, hepatic encephalopathy, and/or hepatorenal syndrome. Liver cirrhosis, which included compensated and decompensated cirrhosis, was identified by the ICD-9-CM diagnosis codes for cirrhosis and its related complications. The diagnosis of HCC was established based on liver tumor features seen on dynamic imaging techniques, tumor histology, and/or high alpha-fetoprotein level (above 400 ng/mL) in LSM and MRS cohorts, and by the ICD-9-CM diagnosis codes in NAFLD cohort. Only the first occurrence of these events was studied. Patients who did not develop study endpoints due to death or loss to follow-up were censored at the date of their last recorded blood test or clinical encounter.
TE
LSM and CAP measurements were performed by experienced operators using the FibroScan502 Touch machine (Echosens, Paris, France) as described previously. All patients were fasted for at least 2 hours before the procedure. The final LSM (in kPa) and CAP values (in dB/m) were represented by the median of ≥10 measurements. LSM was considered reliable only if at least 10 successful acquisitions were obtained and the interquartile range (IQR)-to-median ratio of the 10 acquisitions was ≤0.3 for LSM values of ≥7.1 kPa [
15].
Statistical analysis was performed by Statistical Package for Social Science (SPSS version 25.0; IBM, Armonk, NY, USA) and R version 4.0.5; RStudio, Boston, MA, USA). Continuous variables were expressed in mean±standard deviation or median (IQR) as appropriate. Differences between subgroups were analyzed using chi-squared test or Fisher’s exact test for categorical parameters and one-way analysis of variance and Kruskal Wallis test for continuous parameters as appropriate.
Non-linear association
The non-linear relationship between urea level and LSM was evaluated on a spline curve based on generalized additive model in the LSM cohort. The associations between urea level and hepatic decompensations/incident cirrhosis/HCC were evaluated on a continuous scale with restricted cubic spline based on Cox proportional hazards models with four knots at the 5th, 35th, 65th and 95th percentiles [
16].
Stratification of LSM cohort by urea level
Rather than dichotomizing our population at a single cutoff of urea level, it would be more clinically relevant to define abnormally low or high urea levels because both conditions may be associated with increased risk of LREs [
17]. We first divided the LSM cohort into the derivation and validation cohorts with the ratio of 1:1 by the simple randomization method. In the derivation cohort, we identified a dividing point with the smallest hazard ratio (HR) for hepatic decompensation. In other words, the risk of hepatic decompensation began to increase in patients with a urea level below or above that dividing point. We then selected two cutoffs among patients below or above the dividing point. The optimal cutoffs for urea to predict hepatic decompensation were determined by the ‘survminer’ package from R by the minimum
P-value-based approach in the derivation cohort. Then, the cutoffs were validated by bootstrap with 5,000 replications in the validation cohort [
18]. Subsequently, patients in the LSM cohort were divided into three groups by the cutoffs, i.e., the low, moderate, and high urea level groups.
Survival analysis
We conducted univariate analyses to determine clinical parameters associated with the study endpoints (hepatic decompensation/incident cirrhosis/HCC). To test the predictive value of urea level based on the derived cutoffs, we categorized urea into low, moderate, or high in these models. The multivariable Cox regression model using a backward stepwise approach was performed among covariates that were associated with study endpoints on univariate analysis (P<0.05) to determine independent factors. Adjusted HR and 95% confidence interval (CI) of the risk factors were computed. The Kaplan-Meier method was used to estimate the cumulative probability of study endpoints. All statistical tests were two-sided. Statistical significance was taken as P<0.05.
The non-linear association between urea level and LSM was validated in the MRS cohort, and the associations between different levels of urea and hepatic decompensation/incident cirrhosis/HCC were validated in the NAFLD cohort.
RESULTS
The LSM cohort
During the study period, a total of 5,078 patients with known or suspected CLD underwent TE (
Fig. 1). We excluded 292 patients with cancer (5.8%), 503 (9.9%) with missing urea value, and 19 (0.3%) who had hepatic decompensation before baseline. A further 239 patients (4.7%) were excluded from the analysis for incident cirrhosis, because they had LSM >15 kPa at baseline. In the final cohort of 4,282 subjects (
Table 1), 55.3% were males and the mean age was 55.6 years.
The median follow-up was 5.5 years (IQR, 3.0–6.4). During the follow-up period, 46 patients (1.0%) developed hepatic decompensation, with ascites being the most common initial decompensating event. Twenty-nine patients developed new-onset ascites, including 20 in the moderate urea level group, four in the low urea level group, and five in the high urea level group. Cirrhosis was the main contributor of ascites in those patients. One hundred patients (2.3%) developed HCC.
Supplementary Table 2 lists the number of patients with different LREs.
A U-shaped association between LSM and urea level was identified in the LSM cohort (
Supplementary Fig. 1A). As shown in the
Figure 2, urea showed strong U-shaped associations with decompensation and incident cirrhosis after adjustment with confounding variables. For both outcomes, the risk was the lowest when the urea level was around 6 mmol/L and then increased in both directions (
P for non-linearity <0.05). There was no significant non-linear association between urea level and HCC (data not shown).
Additionally, we performed the analysis to study for the association between creatinine level/urea-to-creatinine ratio (UCR) and hepatic decompensation. As showed in the
Supplementary Figure 2, there was no significant non-linear association between creatinine level and hepatic decompensation. There was a U-shaped association between the UCR and risk of decompensation, however, the non-linear association between UCR and risk decompensation was mainly contributed by the urea level.
Next, we determined the best cutoffs of urea to identify patients at risk of hepatic decompensation. As indicated in the previous section, a urea level of 6 mmol/L was the dividing point with the risk of incident cirrhosis and hepatic decompensation starting to increase in both directions. We divided the LSM cohort into the derivation and validation cohorts (
Supplementary Table 3). Based on the minimum
P-value-based approach and bootstrap validation with 5,000 replications, the urea cutoffs of 3.5 mmol/L (63 mg/dL) and 10 mmol/L (180 mg/dL) were selected as optimal to predict hepatic decompensation, and the patients were divided into low (≤3.5 mmol/L), moderate (3.6–9.9 mmol/L), and high (≥10 mmol/L) urea groups (
Supplementary Fig. 3).
At the derived cutoffs, 344 (8.0%), 3,760 (87.8%), and 178 patients (4.2%) were in the low, moderate, and high urea level groups respectively. The median urea level in the LSM cohort was 5.1 mmol/L (IQR, 4.3–6.2), and 3.2 mmol/L (IQR, 2.9–3.4), 5.2 mmol/L (IQR, 4.4-6.2), 12.9 mmol/L (IQR, 10.9–16.9) in the low, moderated, and high urea level groups, respectively. Patients in the low urea group were younger and had fewer comorbidities compared with the moderate and high urea groups (
Table 1).
Compared with the moderate urea group, both low and high urea levels predicted a significantly higher risk of hepatic decompensation (
Fig. 3A, log-rank
P<0.001). Significant variables associated with hepatic decompensation on univariate analysis were low or high urea levels, older age, diagnosis of viral hepatitis or NAFLD, hypercholesterolemia, elevated alanine aminotransferase, elevated bilirubin, low serum albumin, low platelet count, and increased LSM (
Table 2). On multivariate analysis, urea level remained one of independent predictors of hepatic decompensation. Patients with low urea level (HR, 4.15; 95% CI, 1.68–10.24) had a comparable risk of decompensation as those with high urea level (HR, 5.22; 95% CI, 1.86–14.67) when using the moderate urea level as the reference. Other independent predictors of hepatic decompensation were elevated alanine aminotransferase, low serum albumin, low platelet count, and increased LSM.
High and low urea levels were also associated with incident cirrhosis (
Fig. 3B). After adjusting for potential confounding variables, low urea level was still significantly associated with incident cirrhosis (HR, 3.24; 95% CI, 1.50–6.98), but not high urea level (HR, 1.45; 95% CI, 0.49–4.29). As shown in the
Supplementary Figure 4 and
Supplementary Table 4, there was no significant association between the urea level and HCC.
As shown in the
Figure 4, both low and high urea levels were associated with significantly higher risk of hepatic decompensation compared with the moderate urea group in almost all subgroups. In patients with hypercholesterolemia, the HR for hepatic decompensation in patients with low urea level was of similar magnitude as the other subgroups, but fell short of statistical significance.
Likewise, both low and high urea levels were associated with incident cirrhosis in most subgroups, with the exception of a few subgroups with small patient numbers (
Fig. 5). Using a more conservative LSM cutoff of ≥12 kPa to exclude patients with pre-existing cirrhosis, a low urea level at baseline remained significantly associated with incident cirrhosis.
The MRS cohort was collected from the general population, and those subjects had very low risk of hepatic decompensation (
Supplementary Table 5). Therefore, we used this cohort to demonstrate a U-shaped association between urea level and LSM (
Supplementary Fig. 1B).
The baseline characteristics of the NAFLD cohort are shown in
Supplementary Table 6. One thousand eight hundred thirty-eight (14.7%), 10,091 (80.9%), and 547 patients (4.4%) were in the low, moderate, and high urea level groups, respectively. The median follow-up was 9.9 years (IQR, 8.2–12.1). Compared with the moderate urea group, high urea level at baseline was associated with increased incidence of hepatic decompensation and incident cirrhosis (
Supplementary Fig. 5, log-rank
P<0.001). Low urea level also increased the risk of hepatic decompensation and incident cirrhosis (
Table 3). By multivariate analysis, high urea level remained as an independent factor associated with hepatic decompensation (adjusted HR, 3.71), but low urea level was no longer significant after adjustment. In contrast, both high and low urea levels were independent factors for incident cirrhosis. Similar to what was observed in the LSM cohort, low urea level was not associated with HCC in the NAFLD cohort (
Supplementary Fig. 6,
Supplementary Table 7).
DISCUSSION
To our knowledge, this is the first study to describe a U-shaped relationship between urea level and the severity and outcome of CLD. We showed that patients with urea level of 3.5–10 mmol/L had the lowest risk of advanced liver disease and hepatic decompensation, while the risk was increased when urea was below or above this range. This is close to the normal range of urea used in different laboratories. Importantly, while clinicians focused on the use of high urea level to identify patients with renal dysfunction previously, we provided clinical meaning to low urea level.
Abnormal urea levels are usually caused by abnormal production or excretion of urea. Ureagenesis, which covers the final and irreversible transformation of amino nitrogen to urea nitrogen, is an essential function of the liver involved in whole body nitrogen homeostasis [
19]. The capacity for urea synthesis is decreased in patients with CLD due to several etiologies. Glavind and colleagues [
20] showed that patients with alcoholic hepatitis had marked reduction in urea synthesis compared with healthy controls using functional hepatic nitrogen clearance (FHNC). Reduced FHNC was also observed in patients with viral hepatitis and NAFLD [
9,
21]. The expression of enzymes involved in the urea cycle is also reduced in pre-cirrhotic CLD [
10]. Furthermore, reduced ureagenesis is associated with the accumulation of ammonia, which in turn activates hepatic stellate cells and promotes fibrosis progression [
11]. Nonetheless, previous studies lacked clinical outcome data, and it is impractical to apply FHNC and enzyme expression assays in routine clinical settings. Therefore, the present study offers a simple, yet important message on the use of low urea level to identify and predict advanced liver disease.
Unlike a mechanistic study, it is difficult to determine if low urea was the cause or result of advanced liver disease, though previous studies, as discussed above, suggest it might be both. In very advanced liver disease, low albumin levels reduce ureagenesis by the downregulation of protein metabolism [
22]. However, we excluded patients with pre-existing hepatic decompensation in all three cohorts, and the association between low urea level and clinical outcomes remained robust after adjustment for the albumin level.
In terms of potential clinical application, identification of these patients with low urea may prompt consideration for earlier intervention. In viral hepatitis patients, low urea level may herald the need to commence antiviral therapy if not already done so. Indeed, in our current study, patients with moderate urea levels had the highest rates of antiviral therapy (
Table 1). Furthermore, Laursen et al. [
23] showed early normalization of reduced urea synthesis capacity after commencing antiviral therapy. This is not to say clinicians should initiate antiviral therapy based on the urea level. We should follow the treatment guidelines instead. A low urea level may be a marker of advanced liver disease and prompt clinicians to perform further tests to assess the disease severity and guide treatment. In patients with alcohol-related liver disease and NAFLD, those with low urea levels could be prioritized for referral to more intensive alcohol cessation counselling or NAFLD clinical trials, respectively. As discussed above, patients with low urea levels should be assessed for malnutrition and sarcopenia, and commenced on aggressive nutritional supplementation if present. Clearly these hypotheses require future study to evaluate if urea levels can be improved with the above interventions and if they correspond with improvements in patient prognosis.
With regard to patients with high urea levels, it is well known that urea is primarily excreted in the urine and renal dysfunction results in its accumulation in blood. Renal dysfunction is common in patients with acute and CLDs [
24]. In the present study, more than 90% of patients with high urea level had chronic kidney disease (CKD), and CKD was a risk factor for hepatic decompensation. However, high urea level continued to be associated with hepatic decompensation after we excluded patients with CKD. Previous studies demonstrated that stressful situations like pain, recent surgery, uncontrolled diabetes, and inflammation located outside of the liver can also lead to upregulated ureagenesis [
20].
Our study has several limitations. First, although the LSM and MRS cohorts were prospectively recruited, the NAFLD cohort was retrospective in nature and was based on a real-world registry. However, the latter provided a huge sample size with a long follow-up duration. Second, although we validated the findings in multiple independent cohorts, our study only included ethnic Chinese, and therefore our findings may not be generalizable to other populations. Third, we assessed our patients using TE, which is a good but imperfect noninvasive test for liver fibrosis. However, it would be impractical and unethical to conduct liver biopsy on all subjects, particularly those in the general population. Fourth, the studied population was heterogeneous and the number of primary outcomes was small. To demonstrate consistent results, we have performed various subgroup analysis and adjusted for important confounding factors. Finally, treatments of CLD might alter the natural history and in turn affect the urea level. However, the present studies included CLD patients of different etiologies, and it would be difficult to adjust for the varying types and efficacies of treatments.
In conclusion, we identified a U-shaped relationship between the urea level and liver fibrosis/incident cirrhosis/hepatic decompensation in patients with CLD. Clinicians should recognize low urea level as a potential marker and predictor of advanced liver disease.
ACKNOWLEDGMENTS
This study was supported in part by a direct grant from The Chinese University of Hong Kong (2019.043).
FOOTNOTES
-
Authors’ contributions
Guarantor of the article: Vincent Wai-Sun Wong
Specific author contributions: Vincent Wong has full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Huapeng Lin was responsible for the study concept and design. Huapeng Lin, Grace Wong, Xinrong Zhang, Terry Yip, Ken Liu, Vicki Hui and Vincent Wong were responsible for the acquisition of data. Huapeng Lin, Terry Yip, Yee Kit Tse and Vincent Wong were responsible for data analysis. Huapeng Lin and Vincent Wong drafted the manuscript. All authors had significant contribution to the intellectual content of the manuscript and approved its final version.
-
Conflicts of Interest
Grace Wong has served as an advisory committee member for Gilead Sciences and Janssen, and as a speaker for Abbott, Abbvie, Bristol-Myers Squibb, Echosens, Gilead Sciences, Janssen and Roche.
Terry Yip has served as an advisory committee member and a speaker for Gilead Sciences.
Henry Chan has served as an advisor for AbbVie, Aligos, Arbutus, Hepion, GSK, Gilead Sciences, Janssen, Merck, Roche, Vaccitech, VenatoRx, Vir Biotechnology; and as a speaker for Gilead Sciences, Mylan, and Roche.
Vincent Wong has served as a consultant or advisory board member for 3V-BIO, AbbVie, Allergan, Boehringer Ingelheim, Center for Outcomes Research in Liver Diseases, Echosens, Gilead Sciences, Hanmi Pharmaceutical, Intercept, Inventiva, Merck, Novartis, Novo Nordisk, Perspectum Diagnostics, Pfizer, ProSciento, Sagimet Biosciences, TARGET PharmaSolutions, and Terns; and a speaker for Abbott, AbbVie, Bristol-Myers Squibb, Echosens, and Gilead Sciences. He has received a grant from Gilead Sciences for fatty liver research.
SUPPLEMENTAL MATERIAL
Supplementary material is available at Clinical and Molecular Hepatology website (
http://www.e-cmh.org).
Supplementary Table 1.
ICD-9-CM diagnosis codes for hepatic decompensation, liver cirrhosis, and HCC
Supplementary Table 3.
The baseline characteristics of the patients from the derivation and validation cohorts (LSM cohort)
Supplementary Table 5.
Clinical characteristics of patients with different urea levels in the MRS cohort
Supplementary Table 6.
Clinical characteristics of patients with different urea levels in the NAFLD cohort.
Supplementary Figure 1.
A non-linear relationship between urea level (mmol/L) and LSM (kPa) in the (A) LSM and (B) MRS cohorts. The LSM was plotted against urea level with a curve indicating the U-shaped relationship between the urea level and LSM after adjustment with patients age, gender, viral hepatitis infection status, NAFLD, alanine transaminase, total bilirubin, and albumin levels. LSM, liver stiffness measurement; Ref., reference; CI, confidence interval; MRS, magnetic resonance spectroscopy; NAFLD, nonalcoholic fatty liver disease.
Supplementary Figure 2.
The restricted cubic spline on the association between (A) serum creatinine level, (B) creatinine level after log transformation, (C) urea-to-creatinine ratio and hepatic decompensation in the liver stiffness measurement cohort. The solid red lines represented the adjusted hazard ratios (natural logarithm value), with dashed lines representing the 95% confidence intervals (CIs) derived from restricted cubic spline with four knots. Reference lines for no association are indicated by the dashed lines at a log HR of zero. The purple regions showed the fraction of the population. HR, hazard ratio; Cr, creatinine; UCR, urea-to-creatinine ratio.
Supplementary Figure 3.
Kaplan-Meier analysis of the cumulative probability of hepatic decompensation in the (A) derivation and (B) validation cohorts from liver stiffness measurement cohort. As the cutoff points (3.5 and 10 mmol/L) were selected out, the probability of hepatic decompensation was compared in patients with urea level ≤3.5 mmol/L (low level), 3.6–9.9 mmol/L (moderate level), and >10 mmol/L (high level) in the derivation and validation cohorts.
Supplementary Figure 4.
Kaplan-Meier analysis of the cumulative probability of hepatocellular carcinoma among patients with low (≤3.5 mmol/L), moderate (3.6–9.9 mmol/L), or high (≥10 mmol/L) urea levels in liver stiffness measurement cohort.
Supplementary Figure 5.
Kaplan-Meier analysis of the cumulative probability of (A) hepatic decompensation and (B) incident cirrhosis among patients with low (≤3.5 mmol/L), moderate (3.6–9.9 mmol/L), or high (≥10 mmol/L) urea levels in the nonalcoholic fatty liver disease cohort.
Supplementary Figure 6.
Kaplan-Meier analysis of the cumulative probability of hepatocellular carcinoma among patients with low (≤3.5 mmol/L), moderate (3.6–9.9 mmol/L), or high (≥10 mmol/L) urea levels in the nonalcoholic fatty liver disease cohort.
Figure 1.Study flow diagram. LSM, liver stiffness measurement; HCC, hepatocellular carcinoma.
Figure 2.A U-shaped relationship between urea level and (A) hepatic decompensation and (B) incident cirrhosis after adjustment for the patients’ age, gender, viral hepatitis infection status, NAFLD, alanine transaminase, total bilirubin, and albumin levels in the LSM cohort. The solid red lines represented the adjusted hazard ratios (natural logarithm value), with dashed lines representing the 95% confidence intervals derived from restricted cubic spline with four knots. Reference lines for no association are indicated by the dashed lines at a log HR of zero. The purple regions showed the fraction of the population with different levels of urea. The black points indicated the concentration of urea level with the lowest risk of hepatic decompensation. HR, hazard ratio; Ref., reference; CI, confidence interval; NAFLD, nonalcoholic fatty liver disease; LSM, liver stiffness measurement.
Figure 3.Kaplan-Meier analysis of the cumulative probability of (A) hepatic decompensation and (B) incident cirrhosis among patients with low (≤3.5 mmol/L), moderate (3.6-9.9 mmol/L), or high (≥10 mmol/L) urea levels in the liver stiffness measurement cohort.
Figure 4.Subgroup analysis of the association between urea level and hepatic decompensations risk in the liver stiffness measurement cohort. Subgroups of age less than 60 years, female, NAFLD, and chronic kidney disease (CKD) were not included due to a small number of patients with high urea level. CI, confidence interval; NAFLD, nonalcoholic fatty liver disease.
Figure 5.Subgroup analysis of the association between urea level and incident cirrhosis in the liver stiffness measurement (LSM) cohort. Subgroups of age less than 60 years, NAFLD, and chronic kidney disease (CKD) were not included due to a small number of patients with high urea level. CI, confidence interval; NAFLD, nonalcoholic fatty liver disease.
Table 1.Clinical characteristics of patients with different urea levels in the LSM cohort
Table 1.
|
Characteristic |
All |
Low urea level |
Moderate urea level |
High urea level |
P-value |
|
Total |
4,282 |
344 |
3,760 |
178 |
|
|
Age (years) |
55.6±12.3 |
47.7±11.5 |
55.9±12.0 |
65.6±9.9 |
<0.001 |
|
Male sex |
2,367 (55.3) |
125 (36.3) |
2,132 (56.7) |
110 (61.8) |
<0.001 |
|
HBV infection |
1,375 (32.1) |
125 (36.3) |
1,218 (32.4) |
32 (18.0) |
<0.001 |
|
HCV infection |
131 (3.1) |
6 (1.7) |
124 (3.3) |
1 (0.6) |
0.039 |
|
NAFLD |
1,581 (36.9) |
106 (30.8) |
1,410 (37.5) |
65 (36.5) |
0.048 |
|
Type 2 diabetes mellitus |
994 (23.2) |
66 (19.2) |
808 (21.5) |
120 (67.4) |
<0.001 |
|
Hypertension |
2,290 (53.5) |
109 (31.7) |
2007 (53.4) |
151 (84.8) |
<0.001 |
|
Hypercholesterolemia |
2364 (55.2) |
147 (42.7) |
2066 (54.9) |
151 (84.8) |
<0.001 |
|
Body mass index (kg/m2) |
26 (23–28) |
25 (22–28) |
26 (23–28) |
27 (24–29) |
0.021 |
|
ALT (IU/L) |
26 (18–39) |
26 (16–44) |
26 (19–39) |
21 (16–29) |
<0.001 |
|
Albumin (g/L) |
43.6±3.1 |
43.4±3.4 |
43.8±3.0 |
40.2±3.8 |
<0.001 |
|
Total bilirubin (μmol/L) |
10.5 (7.9–14.0) |
10.3 (7.5–14.2) |
10.6 (8.1–14.1) |
7.1 (5.6–10.0) |
<0.001 |
|
ALP (IU/L) |
67 (56–82) |
66 (53–79) |
67 (56–81) |
76 (62–91) |
<0.001 |
|
AFP (μg/L) |
2.8 (1.9–4.1) |
2.7 (1.8–4.3) |
2.8 (2.0–4.1) |
2.0 (1.5–2.8) |
0.034 |
|
WCC (×109/L) |
6.6±1.9 |
6.4±2.1 |
6.6±1.9 |
7.4±1.9 |
<0.001 |
|
Platelet (×109/L) |
214 (176–255) |
229 (188–276) |
213 (174–253) |
214 (183–252) |
<0.001 |
|
Prothrombin time (seconds) |
10.7 (10.1–11.3) |
10.7 (10.3–11.3) |
10.7 (10.1–11.3) |
10.4 (9.8–11.1) |
0.054 |
|
INR |
1.0 (0.9–1.1) |
1.0 (1.0–1.0) |
1.0 (0.9–1.0) |
1.0 (0.9–1.0) |
0.079 |
|
Triglycerides (mmol/L) |
1.2 (0.8–1.7) |
1.1 (0.8–1.5) |
1.2 (0.8–1.7) |
1.5 (1.0–2.4) |
<0.001 |
|
Cholesterol (mmol/L) |
4.5 (3.9–5.2) |
4.7 (4.0–5.3) |
4.5 (3.9–5.2) |
4.3 (3.7–4.8) |
<0.001 |
|
FBS (mmol/L) |
6.0 (5.1–7.7) |
5.3 (4.8–7.0) |
6.0 (5.1–7.7) |
6.9 (5.6–8.5) |
<0.001 |
|
Haemoglobin A1c (%) |
6.7 (5.8–7.8) |
6.1 (5.6–7.5) |
6.7 (5.8–7.7) |
7.6 (6.8–8.7) |
<0.001 |
|
Positive HBeAg*
|
213 (19.1) |
27 (24.5) |
184 (18.6) |
2 (12.5) |
0.253 |
|
HBV DNA* (log10 IU/mL) |
4 (3–6) |
4 (3–6) |
4 (3–6) |
2 (1-5) |
0.554 |
|
Creatinine (µmol/L) |
75 (62–89) |
62 (54–72) |
75 (63–88) |
184 (132–256) |
<0.001 |
|
eGFR (mL/min/1.73 m2) |
92.7 (76.9–104.2) |
107.2 (97.8–115.5) |
92.3 (78.1–103.2) |
29.0 (19.2–42.5) |
<0.001 |
|
Urea (mmol/L) |
5.1 (4.3–6.2) |
3.2 (2.9–3.4) |
5.2 (4.4–6.2) |
12.9 (10.9–16.9) |
<0.001 |
|
LSM (kPa) |
5.9 (4.5–7.9) |
5.4 (4.3–7.2) |
5.9 (4.5–7.9) |
6.3 (4.9–9.0) |
<0.001 |
|
CAP (dB/m) |
262 (217–310) |
255.5 (204.0–303.5) |
263.0 (218.0–311.0) |
254.0 (209.0–309.5) |
0.028 |
|
Medication use |
|
|
|
|
|
|
Statins |
1,807 (42.2) |
88 (25.6) |
1575 (41.9) |
144 (80.9) |
<0.001 |
|
Anti-platelet drugs |
834 (19.5) |
36 (10.5) |
704 (18.7) |
94 (52.8) |
<0.001 |
|
Antiviral treatment |
625 (14.6) |
48 (14.0) |
562 (14.9) |
15 (8.4) |
0.052 |
Table 2.Factors associated with hepatic decompensation and incident cirrhosis in the LSM cohort
Table 2.
|
Variable |
Univariable analysis
|
Multivariable analysis
|
|
HR (95% CI) |
P-value |
HR (95% CI) |
P-value |
|
Hepatic decompensation |
|
|
|
|
|
|
Urea level, low vs. moderate level |
3.52 (1.68–7.41) |
<0.001 |
4.15 (1.68–10.24) |
0.002 |
|
Urea level, high vs. moderate level |
5.56 (2.32–13.37) |
<0.001 |
5.22 (1.86–14.67) |
0.002 |
|
Age, per year increase |
1.05 (1.02–1.08) |
0.001 |
1.04 (1.00–1.08) |
0.027 |
|
Male sex, vs. female |
1.77 (0.95–3.27) |
0.071 |
|
|
|
Viral hepatitis, yes vs. no |
4.06 (2.21–7.46) |
<0.001 |
1.17 (0.47–2.91) |
0.737 |
|
NAFLD, yes vs. no |
0.26 (0.12–0.59) |
0.001 |
0.51 (0.17–1.48) |
0.212 |
|
Diabetes, yes vs. no |
1.82 (1.00–3.32) |
0.050 |
|
|
|
Hypertension, yes vs. no |
0.65 (0.37–1.17) |
0.156 |
|
|
|
Hypercholesterolemia, yes vs. no |
0.55 (0.31–0.99) |
0.047 |
0.81 (0.38–1.75) |
0.591 |
|
ALT, per IU/L increase |
1.03 (1.02-–1.04) |
<0.001 |
1.00 (1.00–1.01) |
0.027 |
|
Albumin, per g/L increase |
0.74 (0.69–0.78) |
<0.001 |
0.87 (0.81–0.93) |
<0.001 |
|
Total bilirubin, per μmol/L increase |
1.07 (1.05–1.09) |
<0.001 |
0.99 (0.97–1.02) |
0.547 |
|
Platelet, per 1×109/L increase |
0.98 (0.97–0.99) |
<0.001 |
0.99 (0.98–0.99) |
<0.001 |
|
LSM, per kPa increase |
1.09 (1.08–1.10) |
<0.001 |
1.07 (1.05–1.10) |
<0.001 |
|
Incident cirrhosis |
|
|
|
|
|
Urea level, low vs. moderate level |
3.26 (1.67–6.34) |
0.001 |
3.24 (1.50–6.98) |
0.003 |
|
Urea level, high vs. moderate level |
3.60 (1.42–9.12) |
0.007 |
1.45 (0.49–4.29) |
0.508 |
|
Age, per year increase |
1.03 (1.01–1.06) |
0.010 |
1.03 (0.99–1.05) |
0.079 |
|
Male sex, vs. female |
1.09 (0.65–1.85) |
0.726 |
|
|
|
Viral hepatitis, yes vs. no |
3.48 (2.04–5.92) |
<0.001 |
0.72 (0.36–1.41) |
0.334 |
|
NAFLD, yes vs. no |
0.14 (0.06–0.34) |
<0.001 |
0.19 (0.07–0.51) |
0.001 |
|
Diabetes, yes vs. no |
0.54 (0.25–1.14) |
0.103 |
|
|
|
Hypertension, yes vs. no |
0.66 (0.39–1.11) |
0.119 |
|
|
|
Hypercholesterolemia, yes vs. no |
0.35 (0.20–0.61) |
<0.001 |
0.57 (0.30–1.09) |
0.087 |
|
ALT, per IU/L increase |
1.00 (1.00–1.01) |
0.059 |
|
|
|
Albumin, per g/L increase |
0.86 (0.79–0.93) |
<0.001 |
0.90 (0.82–0.99) |
0.031 |
|
Total bilirubin, per μmol/L increase |
1.05 (1.02–1.07) |
<0.001 |
1.03 (0.99–1.06) |
0.126 |
|
Platelet, per 1×109/L increase |
0.98 (0.98–0.99) |
<0.001 |
0.99 (0.98–0.99) |
0.004 |
|
LSM, per kPa increase |
1.35 (1.25–1.45) |
<0.001 |
1.32 (1.21–1.44) |
<0.001 |
Table 3.Factors associated with hepatic decompensation and incident cirrhosis in the NAFLD cohort
Table 3.
|
Variable |
Univariable analysis
|
Multivariable analysis
|
|
HR (95% CI) |
P-value |
HR (95% CI) |
P-value |
|
Hepatic decompensation |
|
|
|
|
|
|
Urea level, low vs. moderate level |
1.63 (1.09–2.43) |
0.016 |
1.19 (0.78–1.80) |
0.421 |
|
Urea level, high vs. moderate level |
10.84 (6.72–17.46) |
<0.001 |
3.71 (2.21–6.23) |
<0.001 |
|
Age, per year increase |
1.06 (1.04–1.07) |
<0.001 |
1.03 (1.01–1.04) |
<0.001 |
|
Male sex, vs. female |
1.98 (1.43–2.73) |
<0.001 |
1.71 (1.22–2.39) |
0.002 |
|
Diabetes, yes vs. no |
2.57 (1.79–3.70) |
<0.001 |
2.20 (1.50–3.22) |
<0.001 |
|
Hypertension, yes vs. no |
1.78 (1.25–2.55) |
0.002 |
1.14 (0.77–1.67) |
0.520 |
|
ALT, per IU/L increase |
0.99 (0.99–1.00) |
0.038 |
0.99 (0.99–1.00) |
0.012 |
|
Albumin, per g/L increase |
0.85 (0.83–0.87) |
<0.001 |
0.87 (0.85–0.89) |
<0.001 |
|
Total bilirubin, per μmol/L increase |
1.01 (1.00–1.01) |
<0.001 |
1.01 (1.00–1.01) |
<0.001 |
|
Platelet, per 1×109/L increase |
0.99 (0.98–0.99) |
<0.001 |
0.99 (0.99–1.00) |
<0.001 |
|
Incident cirrhosis |
|
|
|
|
|
Urea level, low vs. moderate level |
1.58 (1.30–1.93) |
<0.001 |
1.27 (1.03–1.57) |
0.024 |
|
Urea level, high vs. moderate level |
4.28 (2.98–6.16) |
<0.001 |
2.26 (1.56–3.28) |
<0.001 |
|
Age, per year increase |
1.05 (1.04–1.06) |
<0.001 |
1.03 (1.03–1.04) |
<0.001 |
|
Male sex, vs. female |
1.63 (1.39–1.90) |
<0.001 |
1.41 (1.19–1.67) |
<0.001 |
|
Diabetes, yes vs. no |
1.26 (1.07–1.48) |
0.006 |
1.10 (0.93–1.30) |
0.263 |
|
Hypertension, yes vs. no |
1.27 (1.08–1.51) |
0.005 |
0.97 (0.81–1.16) |
0.718 |
|
ALT, per IU/L increase |
1.00 (1.00–1.01) |
<0.001 |
1.00 (1.00–1.01) |
<0.001 |
|
Albumin, per g/L increase |
0.88 (0.87–0.89) |
<0.001 |
0.91 (0.89–0.93) |
<0.001 |
|
Total bilirubin, per μmol/L increase |
1.01 (1.00–1.01) |
<0.001 |
1.00 (1.00–1.01) |
0.002 |
|
Platelet, per 1×109/L increase |
0.99 (0.98–0.99) |
<0.001 |
0.99 (0.98–0.99) |
<0.001 |
Abbreviations
controlled attenuation parameter
Clinical Data Analysis and Reporting System
functional hepatic nitrogen clearance
International Classification of Diseases
liver stiffness measurement
magnetic resonance spectroscopy
nonalcoholic fatty liver disease
REFERENCES
- 1. Younossi ZM, Stepanova M, Younossi Y, Golabi P, Mishra A, Rafiq N, et al. Epidemiology of chronic liver diseases in the USA in the past three decades. Gut 2020;69:564-568.
- 2. Wong VW, Adams LA, de Lédinghen V, Wong GL, Sookoian S. Noninvasive biomarkers in NAFLD and NASH - current progress and future promise. Nat Rev Gastroenterol Hepatol 2018;15:461-478.
- 3. Tsochatzis EA, Bosch J, Burroughs AK. Liver cirrhosis. Lancet 2014;383:1749-1761.
- 4. Yip TC, Chan HL, Tse YK, Lam KL, Lui GC, Wong VW, et al. On-treatment improvement of MELD score reduces death and hepatic events in patients with hepatitis B-related cirrhosis. Am J Gastroenterol 2018;113:1629-1638.
- 5. Calzadilla-Bertot L, Vilar-Gomez E, Wong VW, Romero-Gomez M, Aller-de la Fuente R, Wong GL, et al. ABIDE: an accurate predictive model of liver decompensation in patients with nonalcoholic fatty liver-related cirrhosis. Hepatology 2021;73:2238-2250.
- 6. Lee HW, Yip TC, Tse YK, Wong GL, Kim BK, Kim SU, et al. Hepatic decompensation in cirrhotic patients receiving antiviral therapy for chronic hepatitis B. Clin Gastroenterol Hepatol 2021;19:1950-1958.e7.
- 7. Pundir CS, Jakhar S, Narwal V. Determination of urea with special emphasis on biosensors: a review. Biosens Bioelectron 2019;123:36-50.
- 8. Vilstrup H. Synthesis of urea after stimulation with amino acids: relation to liver function. Gut 1980;21:990-995.
- 9. Thomsen KL, Grønbæk H, Glavind E, Hebbard L, Jessen N, Clouston A, et al. Experimental nonalcoholic steatohepatitis compromises ureagenesis, an essential hepatic metabolic function. Am J Physiol Gastrointest Liver Physiol 2014;307:G295-G301.
- 10. De Chiara F, Heebøll S, Marrone G, Montoliu C, Hamilton-Dutoit S, Ferrandez A, et al. Urea cycle dysregulation in non-alcoholic fatty liver disease. J Hepatol 2018;69:905-915.
- 11. De Chiara F, Thomsen KL, Habtesion A, Jones H, Davies N, GraciaSancho J, et al. Ammonia scavenging prevents progression of fibrosis in experimental nonalcoholic fatty liver disease. Hepatology 2020;71:874-892.
- 12. Lee JS, Adler L, Karathia H, Carmel N, Rabinovich S, Auslander N, et al. Urea cycle dysregulation generates clinically relevant genomic and biochemical signatures. Cell 2018;174:1559-1570.e22.
- 13. Liu K, Wong VW, Lau K, Liu SD, Tse YK, Yip TC, et al. Prognostic value of controlled attenuation parameter by transient elastography. Am J Gastroenterol 2017;112:1812-1823.
- 14. Wong VW, Chu WC, Wong GL, Chan RS, Chim AM, Ong A, et al. Prevalence of non-alcoholic fatty liver disease and advanced fibrosis in Hong Kong Chinese: a population study using protonmagnetic resonance spectroscopy and transient elastography. Gut 2012;61:409-415.
- 15. Boursier J, Zarski JP, de Ledinghen V, Rousselet MC, Sturm N, Lebail B, et al. Determination of reliability criteria for liver stiffness evaluation by transient elastography. Hepatology 2013;57:1182-1191.
- 16. Lee DH, Keum N, Hu FB, Orav EJ, Rimm EB, Willett WC, et al. Predicted lean body mass, fat mass, and all cause and cause specific mortality in men: prospective US cohort study. BMJ 2018;362:k2575.
- 17. Salazar MC, Rosen JE, Wang Z, Arnold BN, Thomas DC, Herbst RS, et al. Association of delayed adjuvant chemotherapy with survival after lung cancer surgery. JAMA Oncol 2017;3:610-619.
- 18. Wu W, Xu W, Sun W, Zhang D, Zhao J, Luo Q, et al. Forced vital capacity predicts the survival of interstitial lung disease in anti-MDA5 positive dermatomyositis: a multi-centre cohort study. Rheumatology (Oxford) 2021 Mar 25;doi: 10.1093/rheumatology/keab305.
- 19. Eriksen PL, Vilstrup H, Rigbolt K, Suppli MP, Sørensen M, Heebøll S, et al. Non-alcoholic fatty liver disease alters expression of genes governing hepatic nitrogen conversion. Liver Int 2019;39:2094-2101.
- 20. Glavind E, Aagaard NK, Grønbæk H, Møller HJ, Orntoft NW, Vilstrup H, et al. Alcoholic hepatitis markedly decreases the capacity for urea synthesis. PLoS One 2016;11:e0158388.
- 21. Foster GR, Irving WL, Cheung MC, Walker AJ, Hudson BE, Verma S, et al. Impact of direct acting antiviral therapy in patients with chronic hepatitis C and decompensated cirrhosis. J Hepatol 2016;64:1224-1231.
- 22. Bernardi M, Angeli P, Claria J, Moreau R, Gines P, Jalan R, et al. Albumin in decompensated cirrhosis: new concepts and perspectives. Gut 2020;69:1127-1138.
- 23. Laursen TL, Sandahl TD, Kazankov K, Eriksen PL, Kristensen LH, Holmboe CH, et al. Early normalization of reduced urea synthesis capacity after direct-acting antiviral therapy in hepatitis C cirrhosis. Am J Physiol Gastrointest Liver Physiol 2020;319:G151-G156.
- 24. Ginès P, Schrier RW. Renal failure in cirrhosis. N Engl J Med 2009;361:1279-1290.
