Clin Mol Hepatol > Volume 28(1); 2022 > Article
Lin, Wong, Zhang, Yip, Liu, Tse, Hui, Lai, Chan, and Wong: U-shaped relationship between urea level and hepatic decompensation in chronic liver diseases

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.

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

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.

U-shaped relationship of urea with LSM and 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.

Cutoffs of urea level to predict hepatic decompensation

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).

Factors associated with LREs in the LSM cohort

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.

Subgroup analysis

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.

Validation in the MRS and NAFLD cohorts

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
cmh-2021-0188-suppl1.pdf
Supplementary Table 2.
Liver-related events during follow-up in the LSM cohort
cmh-2021-0188-suppl2.pdf
Supplementary Table 3.
The baseline characteristics of the patients from the derivation and validation cohorts (LSM cohort)
cmh-2021-0188-suppl3.pdf
Supplementary Table 4.
Factors associated with HCC in the LSM cohort
cmh-2021-0188-suppl4.pdf
Supplementary Table 5.
Clinical characteristics of patients with different urea levels in the MRS cohort
cmh-2021-0188-suppl5.pdf
Supplementary Table 6.
Clinical characteristics of patients with different urea levels in the NAFLD cohort.
cmh-2021-0188-suppl6.pdf
Supplementary Table 7.
Factors associated with HCC in the NAFLD cohorts
cmh-2021-0188-suppl7.pdf
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.
cmh-2021-0188-suppl8.pdf
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.
cmh-2021-0188-suppl9.pdf
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.
cmh-2021-0188-suppl10.pdf
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.
cmh-2021-0188-suppl11.pdf
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.
cmh-2021-0188-suppl12.pdf
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.
cmh-2021-0188-suppl13.pdf

Figure 1.
Study flow diagram. LSM, liver stiffness measurement; HCC, hepatocellular carcinoma.

cmh-2021-0188f1.jpg
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.

cmh-2021-0188f2.jpg
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.

cmh-2021-0188f3.jpg
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.

cmh-2021-0188f4.jpg
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.

cmh-2021-0188f5.jpg

cmh-2021-0188f6.jpg
Table 1.
Clinical characteristics of patients with different urea levels in the LSM cohort
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

Values are presented as mean±standard deviation, median (interquartile range), or number (%).

LSM, liver stiffness measurement; HBV, hepatitis B virus; HCV, hepatitis C virus; NAFLD, nonalcoholic fatty liver disease; ALT, alanine aminotransferase; ALP, alkaline phosphatase; AFP, alpha-fetoprotein; WCC, white cell count; INR, international normalized ratio; FBS, fasting blood sugar level; HBeAg, hepatitis B e antigen; eGFR, estimated glomerular filtration rate; CAP, controlled attenuation parameter.

* HBeAg and HBV DNA data are only shown among chronic hepatitis B patients. HCV RNA data are not reported as very few patients who had HCV RNA checked.

Table 2.
Factors associated with hepatic decompensation and incident cirrhosis in the LSM cohort
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

LSM, liver stiffness measurement; HR, hazard ratio; CI, confidence interval; NAFLD, nonalcoholic fatty liver disease; ALT, alanine aminotransferase.

Table 3.
Factors associated with hepatic decompensation and incident cirrhosis in the NAFLD cohort
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

NAFLD, nonalcoholic fatty liver disease; HR, hazard ratio; CI, confidence interval; ALT, alanine aminotransferase.

Abbreviations

CAP
controlled attenuation parameter
CDARS
Clinical Data Analysis and Reporting System
CI
confidence interval
CKD
chronic kidney disease
CLD
chronic liver disease
FHNC
functional hepatic nitrogen clearance
HCC
hepatocellular carcinoma
HR
hazard ratio
ICD-9-CM
International Classification of Diseases
IQR
interquartile range
KM
Kaplan-Meier
LRE
liver-related event
LSM
liver stiffness measurement
MRS
magnetic resonance spectroscopy
NAFLD
nonalcoholic fatty liver disease
TE
transient elastography
UCR
urea-to-creatinine ratio

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.
crossref pmid
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.
crossref pmid
3. Tsochatzis EA, Bosch J, Burroughs AK. Liver cirrhosis. Lancet 2014;383:1749-1761.
crossref pmid
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.
crossref pmid
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.
crossref pmid
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.
crossref pmid
7. Pundir CS, Jakhar S, Narwal V. Determination of urea with special emphasis on biosensors: a review. Biosens Bioelectron 2019;123:36-50.
crossref pmid
8. Vilstrup H. Synthesis of urea after stimulation with amino acids: relation to liver function. Gut 1980;21:990-995.
crossref pmid pmc
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.
crossref pmid
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.
crossref pmid
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.
crossref pmid
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.
crossref pmid pmc
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.
crossref pmid
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.
crossref pmid
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.
crossref pmid
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.
crossref pmid pmc
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.
crossref pmid pmc
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.
crossref
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.
crossref pmid
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.
crossref pmid pmc
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.
crossref pmid
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.
crossref pmid pmc
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.
crossref pmid
24. Ginès P, Schrier RW. Renal failure in cirrhosis. N Engl J Med 2009;361:1279-1290.
crossref pmid

Editorial Office
The Korean Association for the Study of the Liver
Room A1210, 53 Mapo-daero(MapoTrapalace, Dowha-dong), Mapo-gu, Seoul, 04158, Korea
TEL: +82-2-703-0051   FAX: +82-2-703-0071    E-mail: kasl@kams.or.kr
Copyright © The Korean Association for the Study of the Liver.         
COUNTER
TODAY : 2825
TOTAL : 1738833
Close layer