Skip to main content

Advertisement

Non-invasive score identifies ultrasonography-diagnosed non-alcoholic fatty liver disease and predicts mortality in the USA

Article metrics

Abstract

Background

Several non-invasive prediction scores for non-alcoholic fatty liver disease (NAFLD) have been developed, but their performance has not been compared and validated in the same population, and whether these prediction scores can predict clinical outcomes remains unknown. In this study, we aimed to validate and compare the performance of four NAFLD prediction scores: fatty liver index, hepatic steatosis index, lipid accumulation product, and NAFLD liver fat score (LFS), and to evaluate the ability of the best NAFLD prediction score to predict mortality.

Methods

We analyzed data from the National Health and Nutrition Examination Survey conducted in 1988 to 1994, and subsequent follow-up data for mortality up to December 31, 2006. NAFLD was defined by ultrasonographic detection of hepatic steatosis in the absence of other known liver diseases.

Results

In a group of 5,184 participants, LFS consistently showed the highest area under the curve for predicting the presence of NAFLD. During a median follow-up of 14.7 years (range 0.1 to 18.2 years) and 83,830.5 person-years, participants in the high LFS group (LFS ≥1.257) had a higher cardiovascular and liver-related mortality than participants in the low (LFS ≤− 1.413; cardiovascular hazard ratio (HR) = 2.24, 95% CI 1.03 to 4.88; liver HR = 31.25, 95% CI 3.13 to 333.33) or intermediate (-1.413 < LFS < 1.257; cardiovascular HR = 2.3, 95% CI 1.19 to 4.48; liver HR = 30.3, 95% CI 4 to 250) LFS groups in the fully adjusted model. Similar results were obtained when LFS was treated as a continuous variable.

Conclusions

LFS is the best non-invasive prediction score for NAFLD, and people with a high LFS score have an increased risk for cardiovascular and liver-related mortality.

Background

Non-alcoholic fatty liver disease (NAFLD) represents a spectrum of progressive liver disease ranging from simple steatosis to non-alcoholic steatohepatitis (NASH), fibrosis, and cirrhosis, in the absence of excessive alcohol consumption. NAFLD is regarded as a hepatic manifestation of metabolic syndrome (MetS) [1], therefore the presence of NAFLD is not only strongly associated with liver-related mortality, but also with diseases related to the MetS, such as diabetes and cardiovascular diseases [2],[3]. As NAFLD is highly prevalent and affects up to 30% of the general adult population [4], screening for and diagnosing NAFLD has become an important issue in public health to prevent NAFLD-related complications and reduce healthcare costs.

Liver biopsy remains the “gold standard” for NAFLD diagnosis; however, this is an invasive technique making it impractical to be used widely. Ultrasonography is therefore the recommended first-line imaging technique in clinical practice, although it is known to have limited sensitivity [3]. Other non-invasive tools have been developed for diagnosing NAFLD, such as computed tomography and proton magnetic resonance spectroscopy (1H-MRS). However, these tools are expensive and time-consuming, and are not considered cost-effective for large-scale NAFLD screening. Recently, five biomarker-based non-invasive prediction scores of NAFLD have been developed: SteatoTest [5], fatty liver index (FLI) [6], NAFLD liver fat score (LFS) [7], lipid accumulation product (LAP) [8], and hepatic steatosis index (HSI) [9]. These scores are derived from simple clinical risk factors and biomarkers, and can therefore potentially be used for large-scale NAFLD screening. However, different definitions and techniques were used to define NAFLD in the original studies, and the performances of these scores have not been validated, evaluated, and compared directly in a large general population. In addition, whether these non-invasive scores of NAFLD can predict clinical outcome remains largely unknown.

In this study, we aimed to validate and evaluate the performance of these non-invasive prediction scores of NAFLD in predicting ultrasonography-diagnosed NAFLD in a representative general adult population in the USA (cross-sectional NAFLD prediction cohort), and to test if the marker can predict mortality in the general population (prospective mortality prediction cohort).

Methods

Participant recruitment

Data from the third National Health and Nutrition Examination Survey (NHANES III) were used [10]. NHANES III was conducted by the National Center for Health Statistics (NCHS) from 1988 to 1994, using a stratified multistage probability sample that represented the civilian non-institutionalized population in the USA. Participants gave written consent before participation, and ethics approval was obtained from the Human Subjects Committee of the US Department of Health and Human Services.

We studied people aged 20 to 74 years who participated in the NHANES III survey. Laboratory tests were carried out in a mobile examination center (n = 14,797) (see Additional file 1: Figure S1). Because all non-invasive score formulae require levels of biomarkers in fasting blood, we included participants with blood taken after fasting for at least 8 hours fasting (n = 9,268). Of those, participants with factors that can confound the diagnosis of NAFLD (including excessive alcohol consumption, defined as >21 drinks/week in men and >14 drinks/week in women [4]; viral hepatitis, defined as positive serum hepatitis B surface antigen and positive serum hepatitis C antibody; iron overload, defined as transferrin saturation ≥50%; or pregnancy) were excluded (n = 1,089). The LFS formula includes fasting insulin level, therefore we further excluded participants who were using insulin or other medications for diabetes (n = 315). This left 7,864 participants, and once the appropriate exclusion criteria were adopted for each aim (see Additional file 1: Figure S1), we had 5,184 and 5,892 participants included in the analysis of Aim 1 (evaluation of the performance of non-invasive prediction scores of NAFLD in predicting ultrasound-diagnosed NAFLD) and Aim 2 (evaluation of the relationship between non-invasive prediction scores of NAFLD and mortality), respectively.

Definition of NAFLD

In the original NHANES III between 1988 and 1994, gall bladder ultrasonography video images were recorded using a Toshiba Sonolayer SSA-90A and Toshiba video recorder. Between 2009 and 2010, hepatic steatosis (fatty liver) was assessed by archived video images being re-reviewed by three ultrasonography readers (trained by a board-certified radiologist specializing in hepatic imaging), who graded the presence of fat within the hepatic parenchyma.

The following information was recorded on a standard paper collection form: 1) presence of liver-to-kidney contrast; 2) degree of the brightness of the liver parenchyma; 3) presence of deep beam attenuation; 4) presence of echogenic walls in the small intrahepatic vessels, and (5) definition of the gallbladder walls. Finally, an overall primary finding was given based on the presence or absence of each of the five parameters. The liver was graded as having no, mild, moderate, or severe hepatic steatosis. Of the 13,983 participants with hepatic imaging records, 13,856 of them could be graded [7]. Detailed descriptions and procedures have been provided previously [7],[11]. In the absence of a standard definition, we defined NAFLD as moderate or severe hepatic steatosis, and non-NAFLD as no or mild hepatic steatosis [12]. The overall intra-rater and inter-rater - statistics for reliability of the dichotomized outcomes (“no or mild” and “moderate or severe”) were 0.77 (95% CI 0.73 to 0.82) and 0.70 (95% CI 0.64 to 0.76), respectively [11].

Non-invasive markers of NAFLD

The non-invasive markers of NAFLD were calculated based on the equations reported in the literature [6]-[9]. In brief, FLI includes body mass index (BMI), γ-glutamyltranspeptidase, triglycerides, and waist circumference; HSI includes aspartate aminotransferase (AST)/alanine aminotransferase (ALT) ratio, BMI, diabetes, and sex; LAP includes sex, triglycerides, and waist circumference; and LFS includes AST/ALT ratio, diabetes, fasting AST level, fasting insulin level, and MetS. The SteatoTest was not included in the current study, as this test is a commercially one, and the calculation formula is not disclosed. The threshold used in the current study also adopted the cutoff points suggested in the literature: the high/low cutoff points were ≥1.257/≤ -1.413, ≥30/<30, and ≥30/<30 for LFS, HSI, and FLI, respectively [6],[7],[9].

Mortality follow-up

In NHANES III, cause of death was coded using the International Classification of Diseases. 10th Revision (ICD-10). ICD codes I00 to I78 and E10 to E14 were used to assess cardiovascular and diabetes mortality, respectively, as in our previous studies [13],[14]. Malignancy and liver mortality were defined by the Underlying Cause of Death (UCOD)_113 20 to 23, 25 to 26, and 43, and UCOD_113 15, 24, and 93 to 95, respectively, as in the literature [15]. The length of follow-up was the time from the study examination date to death or to December 31, 2006, whichever was earlier.

Definition of diabetes, hypertension, and MetS

Diabetes was defined according to the latest American Diabetes Association (ADA) guideline, which includes fasting glucose ≥126 mg/dl, random plasma glucose ≥200 mg/dl, or A1C ≥ 6.5. Patients were considered to have hypertension if they had systolic blood pressure (SBP) ≥140 or diastolic blood pressure (DBP) ≥90 mmHg, or if they were receiving anti-hypertensive drug therapy. MetS was defined according to the joint scientific statement on harmonizing MetS [16], that is, having three or more of the following factors: 1) elevated blood pressure (SBP ≥ 130 mmHg and/or DBP ≥85 mmHg and/or being in receipt of anti-hypertension drug therapy); 2) elevated triglycerides (≥150 mg/dl (1.7 mmol/l) and/or being in receipt of drug treatment for elevated triglycerides); 3) reduced high-density liproprotein (HDL) cholesterol (<40 mg/dl (1.0 mmol/l) in men and <50 mg/dl (1.3 mmol/l) in women and/or being in receipt of drug treatment for elevated HDL cholesterol); 4) elevated fasting glucose (≥100 mg/dl (5.6 mmol/l) and/or being in receipt of treatment for elevated glucose); and 5) large waist circumference (>102 cm in men and >88 cm in women of European descent). Liver fat percentage was estimated using the equation reported in the same literature as the LFS [7]. The equation includes the same variables as the LFS, but with a different calculation.

Statistical analysis

To assess model discrimination, we calculated the area under curve (AUC) for the receiver operating characteristic (ROC) for each non-invasive score of NAFLD. The difference between two AUCs was compared using the maximum likelihood estimation method [17] and implemented using ROCKIT [18]. Sensitivity, specificity, positive likelihood ratio (+LR), negative likelihood ratio (-LR), and corresponding 95% CIs were also calculated. The non-invasive NAFLD measurement with the best performance (in terms of AUC for ROC) was selected and evaluated for its association with mortality.

In the Cox proportional hazard regression model, non-invasive score was modeled as threshold and continuous variables. Using the lower threshold as the reference, the hazard ratio (HR) and 95% CI for the highest threshold were calculated using the simple and fully adjusted Cox regression models. In the simple model, we adjusted for age and sex. In the full model, we adopted the adjustment model of a recent study related to NAFLD fibrosis and mortality [15], which includes age, sex, race/ethnicity, income, education, diabetes, hypertension, use of lipid-lowering medication, smoking, drinking, history of cardiovascular disease (CVD), waist circumference, dietary caffeine intake, HDL cholesterol, triglycerides, transferrin saturation and C-reactive protein. P ≤ 0.05 was considered significant. The proportional hazards assumption was evaluated by including time-dependent covariates in the full regression model; the overall test of proportional hazards was not significant (P > 0.05) suggesting that the proportional assumption was valid. To gain additional insight into the potential nonlinearity of the effect of LFS, we examined the Cox regression models using penalized spline. Two degrees of freedom (df) used in the spline because the model had the lowest Akaike's information criterion (AIC) (best fit) when df = 2. Sample weights that accounted for the unequal probabilities of selection, oversampling, and non-response were applied in all analyses using the complex sampling module in SPSS (V18.0; SPSS Inc, Chicago, IL, USA) or R software (V2.15.0) [19]. All values presented were weighted to represent the civilian population of the USA.

We also evaluated the ability in risk reclassification using integrated discrimination improvement (IDI) [20] and category-less net reclassification improvement (NRI) [20]. IDI was used to compare the difference in discrimination slopes [21], while category-less NRI was used to compare classifications from two models for changes by outcome for a net calculation of changes in the right direction. Estimated risk of death of different models was calculated using the equation of 1/(1 + exp(-1×XBeta)). Analyses were performed using R software (V2.15.0) [19].

Results

Cross-sectional NAFLD prediction cohort

Of the 5,184 participants included for the AUC study, 18.4% (16.5 to 20.4%) had NAFLD. The characteristics of this cohort are provided in Table 1. For NAFLD prediction, LFS was the best performer for predicting NAFLD, with an AUC of 0.771 (P < 0.001), whereas the lowest AUC (0.732) was observed for HSI (Table 2). Using maximum likelihood estimation, the difference between the AUC of LFS and other markers (FLI, LAP, and HSI) was statistically significant (all P < 0.01). Interestingly, the diagnostic accuracy of these markers differed by race/ethnicity (Table 2). The sensitivity and specificity, and the + LR and -LR of the suggested high and low cutoff points for excluding/including NAFLD are provided in Table 3 and Additional file 1: Table S1, respectively. The raw number used to calculate the diagnostic accuracy and the characteristics of true and false positive, and true and false negative (based on the LFS threshold) are provided in Additional file 1: Tables S2 and S3.

Table 1 Characteristics of participants in the cross-sectional NAFLD prediction cohort according to NAFLD status a
Table 2 Quality of prediction scores in predicting NAFLD
Table 3 Sensitivity and specificity of exclusion/inclusion cutoff points of LFS, FLI, and HSI in the cross-sectional NAFLD prediction cohort and literature

Prospective mortality prediction cohort

Of the several NAFLD non-invasive prediction scores tested, the LFS gave the best performance. As our second aim was to investigate the relationship of the best non-invasive score with outcome, we tested if LFS was associated with mortality. Table 4 shows the characteristics of the participants. During a median follow-up of 14.7 years (range 0.1 to 18.2 years) and 83,830.5 person-years, 793, 311, 209, 58, and 17 participants died from all, cardiovascular-, malignancy-, diabetes-, and liver-related causes, respectively. Higher LFS were associated with all causes of mortality tested, except malignancy-related causes.

Table 4 Characteristics of participants in the prospective mortality prediction cohort according to different LFS thresholds

The results of Cox regression analysis are shown in Table 5. Participants in the high LFS group had a 60% higher risk (HR = 1.6; 95% CI 1.01 to 2.54; P = 0.048 in full model) of all-cause mortality than the intermediate LFS group. For cardiovascular mortality, participants in the high LFS group was associated with 2.24-fold (95% CI 1.03 to 4.88; P = 0.042 in full model) and 2.3-fold (95% CI 1.19 to 4.48; P = 0.015 in full model) increase in risk of death compared with the low and intermediate LFS groups. For liver mortality, participants in the high LFS group had a 31.25-fold (95% CI 3.13 to 333.33; P = 0.004 in full model) and 30.3-fold (95% CI 4 to 250; P = 0.001 in full model) increase in risk of death compared with the low and intermediate LFS groups. The Kaplan-Meier survival curves for cardiovascular and liver-related mortality are provided in Figure 1. When LFS was treated as continuous variable, a one-unit increase of LFS was associated with increased mortality of all-cause and cardiovascular-, liver-, and diabetes-related mortality, with HRs of 1.09 (95% CI 1.01 to 1.19; P = 0.039), 1.11 (95% CI 1.03 to 1.19; P = 0.006), 1.32 (95% CI 1.12 to 1.55; P = 0.001), and 1.21 (95% CI 1.02 to 1.44; P = 0.034), respectively (Table 5), after full adjustment. Similar results were obtained after further adjustment of the NAFLD fibrosis score (NFS). The relationship between LFS and cardiovascular mortality as examined by penalized regression spline is shown in Figure 2.

Table 5 Association between LFS and mortality
Figure 1
figure1

Kaplan-Meier curves of mortality curves according to different liver fat score (LFS) thresholds. (a) Cardiovascular-related and (b) liver-related mortality.

Figure 2
figure2

Association between liver fat score (LFS) and cardiovascular mortality via penalized regression splines.

Age, sex, hypertension, and diabetes are commonly used in assessing risk of mortality. We therefore evaluated whether addition of LFS (categorical: low/intermediate/high risk of NAFLD) in a basic clinical model composed of these traditional risk factors could improve the risk prediction. Risk reclassification with IDI showed a modest positive shift to improvement when LFS was added in the basic clinical model (IDI: 0.0131; 95% CI 0.009 to 0.017; P < 0.001). Similar result was observed using category-less NRI (NRI: 0.133; 95% CI 0.054 to 0.211; P < 0.001).

Discussion

Using a large, nationally representative cohort with more than 10 years of follow-up and ultrasonographic data, we have demonstrated that LFS is the best prediction score for ultrasonography-diagnosed NAFLD, and can predict mortality, including cardiovascular- and liver-related mortality.

It is important to find an easy and cost-effective way to screen for NAFLD. Of the several non-invasive scores we tested, LFS showed the best performance in identifying ultrasonography-diagnosed NAFLD. Notably, LFS was derived using 1H-MRS-diagnosed NAFLD, whereas FLI, HSI, and LAP were derived from ultrasonography-diagnosed NAFLD. For those scores derived from ultrasonography, the definitions of the NAFLD were also somewhat different (see Additional file 1: Table S4). Ultrasonography is a semi-quantitative imaging technique, and the definitions of NAFLD differed between studies. By contrast, 1H-MRS is by far the most sensitive and quantitative imaging tool in identifying hepatic steatosis. This could be the reason why the non-invasive score (LFS) derived from 1H-MRS performed better and more robustly in identifying cases in the current study. As there is no standard definition of ultrasonography-diagnosed NAFLD for good measurement, we used three additional definitions to test the performance of different non-invasive indices (see Additional file 1: Table S5), and the LFS still came out best. We evaluated whether combining all prediction scores (combined score) could improve the NAFLD prediction. The AUC of the combined score increased to 0.782 (95% CI 0.766 to 0.798), suggesting that there are unique NAFLD predicting components being captured in different prediction scores.

Interestingly, there was a difference in the diagnostic accuracy of the different non-invasive scores, with the lowest diagnostic accuracy being observed in black patients for all tested scores, suggesting that the clinical risk factors of NAFLD could be ethnicity-specific and particularly different in black populations. Like other disease predictions [22], deriving an ethnicity- or population-specific prediction model may be required to achieve a high accuracy of NAFLD prediction. Notably, although there was an observed difference in the diagnostic accuracy of LFS for NAFLD, no significant interaction (P > 0.05) between LFS and race/ethnicity on mortality was observed, therefore, no subgroup analysis was performed in the subsequent analyses.

LFS was calculated based on the AST/ALT ratio, diabetes, fasting AST level, fasting insulin level, and MetS. Given that diabetes and MetS are known to be associated with mortality, the association between LFS and mortality could be attributable to these factors. However, the components of MetS were adjusted for in the full model, suggesting that the association of LFS with mortality may be independent of these factors. A number of different organizations use different recommended waist circumference thresholds for abdominal obesity in defining MetS. In addition to the threshold suggested by ATPIII [16], we also used the population-specific threshold suggested by the International Diabetes Federation, and the findings remained unchanged (data not shown).

In the literature on FLI, LFS, and HSI, various high and low cutoff points have been proposed to include and exclude NAFLD [6],[7],[9]. The high cutoff point should have a high specificity and + LR, while the low cutoff point should have a high sensitivity and low −LR. In general, the diagnostic performance of the defined cutoff points of the NAFLD prediction scores was not the same as that originally reported in the literature because the sample populations were different (Table 3). The high cutoff point of LFS had a slightly higher specificity (96.4%) in the current study than the figure (95%) reported in the literature, meaning that study participants with a high LFS were very likely to have a higher risk of mortality.

Two previous studies validated the non-invasive prediction scores in adults [23],[24]. Koehler et al. validated FLI and LAP in 2,652 participants in the Rotterdam Study. FLI and LAP had an AUC of 0.813 and 0.786, respectively, in predicting ultrasonography-diagnosed NAFLD. FLI had a higher AUC in the Rotterdam Study than in the current study. Interestingly, the Rotterdam Study used the scoring protocol of Hamaguchi et al. [25], and the one used by NHANES was an algorithm derived based on that same publication. The Rotterdam Study was a population-based cohort study of elderly inhabitants of a district of Rotterdam, whereas NHANES was a nationally representative population-based study with participants of different races/ethnicities and age. LFS was previously validated in a study of 40 non-diabetic patients with biopsy-proven NAFLD and 85 healthy controls [23], which showed that LFS had an AUC of 0.86. Although the AUC from different validation studies cannot be compared directly, our study is in agreement with previous validation studies showing that LFS had the highest AUC, followed by FLI and LAP. Although no validation study has been performed for HSI, our study suggested that HSI is better than LAP as a predictor of NAFLD.

Although identifying people with NAFLD is important, identifying people with adverse clinical outcome is even more important, as NAFLD consists of a wide spectrum of conditions, ranging from simple steatosis to cirrhosis with varying prognosis. In concordance with previous NHANES reports [12],[15], our study did not reveal any significant association between ultrasonography-diagnosed NAFLD and mortality (data not shown). This finding is intriguing. Ultrasonography-diagnosed NAFLD is not associated with mortality, whereas LFS, a marker of NAFLD, is associated with mortality. This could be due to the reason mentioned earlier, namely, that LFS was derived using the sensitive and quantitative imaging tool 1H-MRS, whereas other NAFLD prediction scores were derived using a less sensitive semi-quantitative ultrasonography. In fact, we found no association of other markers of NAFLD with mortality (se Additional file 1: Table S5), further suggesting that 1H-MRS-derived LFS may be more superior in identifying NAFLD and predicting clinical outcome.

We then investigated which individual component of the LFS was associated with CVD mortality in the multivariable model, and found the only significant association observed Qa with fasting serum insulin (estimate of 1.02; 95% CI 1.01 to 1.03; P < 0.001), suggesting that fasting serum insulin may be the main driver for the observed association.

Notably, high LFS is also associated with low transferrin saturation (Table 4). In our previous study, we showed that low transferrin saturation was robustly associated with pre-diabetes [26]. These observations suggested that elevated insulin resistance might be the key factor leading to mortality in people with NAFLD, which may also explain why LFS can predict mortality whereas ultrasonography-diagnosed NAFLD cannot. Another possibility is that high LFS indicates the presence of other NAFLD-related conditions, such as NASH and fibrosis. In participants with NAFLD, high LFS is associated with high NFS (data not shown), which is a prediction score of NAFLD fibrosis [27], although NFS does not predict NAFLD (data not shown) nor is it significantly associated with mortality in the general population (see Additional file 1: Table S5). However, further adjustment of NFS revealed that the effect of LFS is independent of NFS (Table 5), suggesting that the association between LFS and mortality may be independent of NAFLD fibrosis. Future study is required to confirm our observations and to examine the underlying mechanisms.

Age, sex, and presence of diabetes or hypertension are simple risk factors that are commonly used by clinicians to evaluate mortality risk. We showed that LFS has an independent role in predicting mortality and improved risk reclassification. Similarly, although the Framingham Risk Score (FRS) was not intended for use in mortality prediction, we found that the associations of LFS with cardiometabolic disease related mortality were independent of FRS (see Additional file 1: Table S6). These findings suggest that abnormal liver function may play a role in mortality determination, independently of traditional risk factors.

Our study has several strengths. The study population is large, multiethnic, nationally representative, and well-characterized, with data on ultrasonography-diagnosed NAFLD, multiple risk factors, and potential confounders. The long follow-up and the large number of events provided ample statistical power. The wide range of collected data from NHANES III allowed construction of four different non-invasive prediction scores simultaneously, so that they could be compared in parallel and with different definitions of NAFLD.

Nevertheless, there are limitations. The major limitation of the current study is the use of ultrasonography-diagnosed NAFLD, which can lead to misclassification error. In the absence of a standard definition, we defined NAFLD as presence of moderate or severe hepatic steatosis, while non-NAFLD was defined as presence of no or mild hepatic steatosis, as in previous study. The case definition, especially when mild hepatic steatosis was defined as non-NAFLD, could be a potential source of bias. This classification could have led to underestimation of the NAFLD prevalence in the current study, which was reported to be 20% to 33% in the general population [28], although the lower prevalence observed could also be due to the lower prevalence of obesity in the current study [29]. It is acknowledged that ultrasonography has limited sensitivity and specificity in diagnosing NAFLD, especially when less than 33% of the liver parenchyma is infiltrated by fat [30],[31] and in the presence of liver cirrhosis that may lead to decreased hepatic steatosis. To confirm our findings, we defined NAFLD in different ways and still found that LFS was the best marker of NAFLD (Table 4; see Additional file 1: Table S7), and participants with mild hepatic steatosis also did not have increased mortality compared with those without hepatic steatosis (see Additional file 1: Table S8).

Although liver biopsy is considered the gold standard in diagnosing NAFLD, it is not justifiable to perform liver biopsy in large numbers of asymptomatic individuals, therefore ultrasonography is still considered an acceptable first-line screening procedure for NAFLD in clinical practice [32]. However, ultrasonography cannot distinguish between NASH, fibrosis or cirrhosis.

The prediction score named the SteatoTest, was not included in the current study; whether it is superior to LFS or otherwise requires further study.

In the Cox regression analysis, there were too few cases of liver-related mortality, which led to unreliable estimates, and could also be a potential source of bias, therefore cautious interpretation is required.

Conclusion

In conclusion, we found that 1H-MRS derived NAFLD prediction score LFS was the most robust non-invasive score identifying NAFLD in this US population and predicted mortality. NAFLD is highly prevalent, and can be associated with morbidity and mortality if left unidentified. Our findings suggest that LFS may be a promising tool for large-scale NAFLD screening. If confirmed in future studies, LFS may be a useful marker for large-scale NAFLD screening and prediction of long-term clinical outcomes.

Authors' contributions

CLC was responsible for study conception and design; data generation and analysis; and manuscript drafting. CLC, KSLL, ICKW, BMYC were responsible for data interpretation. All authors reviewed, commented upon, and approved the final submission.

Authors' information

CLC is the article guarantor.

Additional file

References

  1. 1.

    Targher G, Byrne CD: Clinical Review: Nonalcoholic fatty liver disease: a novel cardiometabolic risk factor for type 2 diabetes and its complications. J Clin Endocrinol Metab. 2013, 98: 483-495. 10.1210/jc.2012-3093.

  2. 2.

    Targher G, Day CP, Bonora E: Risk of cardiovascular disease in patients with nonalcoholic fatty liver disease. N Engl J Med. 2010, 363: 1341-1350. 10.1056/NEJMra0912063.

  3. 3.

    Anstee QM, Targher G, Day CP: Progression of NAFLD to diabetes mellitus, cardiovascular disease or cirrhosis. Nat Rev Gastroenterol Hepatol. 2013, 10: 330-344. 10.1038/nrgastro.2013.41.

  4. 4.

    Chalasani N, Younossi Z, Lavine JE, Diehl AM, Brunt EM, Cusi K, Charlton M, Sanyal AJ: The diagnosis and management of non-alcoholic fatty liver disease: practice Guideline by the American Association for the Study of Liver Diseases, American College of Gastroenterology, and the American Gastroenterological Association. Hepatology. 2012, 55: 2005-2023. 10.1002/hep.25762.

  5. 5.

    Poynard T, Ratziu V, Naveau S, Thabut D, Charlotte F, Messous D, Capron D, Abella A, Massard J, Ngo Y, Munteanu M, Mercadier A, Manns M, Albrecht J: The diagnostic value of biomarkers (SteatoTest) for the prediction of liver steatosis. Comp Hepatol. 2005, 4: 10-10.1186/1476-5926-4-10.

  6. 6.

    Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, Tiribelli C: The Fatty Liver Index: a simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006, 6: 33-10.1186/1471-230X-6-33.

  7. 7.

    Kotronen A, Peltonen M, Hakkarainen A, Sevastianova K, Bergholm R, Johansson LM, Lundbom N, Rissanen A, Ridderstråle M, Groop L, Orho-Melander M, Yki-Järvinen H: Prediction of non-alcoholic fatty liver disease and liver fat using metabolic and genetic factors. Gastroenterology. 2009, 137: 865-872. 10.1053/j.gastro.2009.06.005.

  8. 8.

    Bedogni G, Kahn HS, Bellentani S, Tiribelli C: A simple index of lipid overaccumulation is a good marker of liver steatosis. BMC Gastroenterol. 2010, 10: 98-10.1186/1471-230X-10-98.

  9. 9.

    Lee JH, Kim D, Kim HJ, Lee CH, Yang JI, Kim W, Kim YJ, Yoon JH, Cho SH, Sung MW, Lee HS: Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. 2010, 42: 503-508. 10.1016/j.dld.2009.08.002.

  10. 10.

    National Center for Health Statistics (NCHS). National Health and Nutrition Examination Survey Data. 1996, U.S. Department of Health and Human Services, Centers for Disease Control and Prevention [1988-1994], Hyattsville, MD

  11. 11.

    Hepatic Steatosis Ultrasound Images Assessment Procedures Manual: Hepatic Steatosis Ultrasound Images Assessment Procedures Manual. Available from http://www.cdc.gov/nchs/nhanes/nhanes3/HGUHS.htm, accessed 25/01/2014., [http://www.cdc.gov/nchs/data/nhanes/nhanes3/hepatic_steatosis_ultrasound_procedures_manual.pdf]

  12. 12.

    Lazo M, Hernaez R, Bonekamp S, Kamel IR, Brancati FL, Guallar E, Clark JM: Non-alcoholic fatty liver disease and mortality among US adults: prospective cohort study. BMJ. 2011, 343: d6891-10.1136/bmj.d6891.

  13. 13.

    Cheung CL, Lam KS, Cheung BM: Serum beta-2 microglobulin concentration predicts cardiovascular and all-cause mortality. Int J Cardiol. 2013, 168: 4811-4813. 10.1016/j.ijcard.2013.07.014.

  14. 14.

    Cheung CL, Lam KS, Cheung BM: Serum beta-2 microglobulin predicts mortality in people with diabetes. Eur J Endocrinol. 2013, 169: 1-7. 10.1530/EJE-13-0003.

  15. 15.

    Kim D, Kim WR, Kim HJ, Therneau TM: Association between noninvasive fibrosis markers and mortality among adults with nonalcoholic fatty liver disease in the United States. Hepatology. 2013, 57: 1357-1365. 10.1002/hep.26156.

  16. 16.

    Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, Fruchart JC, James WP, Loria CM, Smith SC: Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation. 2009, 120: 1640-1645. 10.1161/CIRCULATIONAHA.109.192644.

  17. 17.

    Metz CE, Herman BA, Shen JH: Maximum likelihood estimation of receiver operating characteristic (ROC) curves from continuously-distributed data. Stat Med. 1998, 17: 1033-1053. 10.1002/(SICI)1097-0258(19980515)17:9<1033::AID-SIM784>3.0.CO;2-Z.

  18. 18.

    Metz CE. ROCKIT. [], [http://metz-roc.uchicago.edu/MetzROC/software]

  19. 19.

    R project. [], [www.R-project.org]

  20. 20.

    Pencina MJ, D'Agostino RB, D'Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27: 157-172. 10.1002/sim.2929. discussion 207-112

  21. 21.

    Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, Pencina MJ, Kattan MW: Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. 2010, 21: 128-138. 10.1097/EDE.0b013e3181c30fb2.

  22. 22.

    Cheung EY, Bow CH, Cheung CL, Soong C, Yeung S, Loong C, Kung A: Discriminative value of FRAX for fracture prediction in a cohort of Chinese postmenopausal women. Osteoporos Int. 2012, 23: 871-878. 10.1007/s00198-011-1647-5.

  23. 23.

    Musso G, Gambino R, Durazzo M, Cassader M: Noninvasive assessment of liver disease severity with liver fat score and CK-18 in NAFLD: prognostic value of liver fat equation goes beyond hepatic fat estimation. Hepatology. 2010, 51: 715-717. 10.1002/hep.23255.

  24. 24.

    Koehler EM, Schouten JN, Hansen BE, Hofman A, Stricker BH, Janssen HL: External validation of the fatty liver index for identifying nonalcoholic fatty liver disease in a population-based study. Clin Gastroenterol Hepatol. 2013, 11: 1201-1204. 10.1016/j.cgh.2012.12.031.

  25. 25.

    Hamaguchi M, Kojima T, Itoh Y, Harano Y, Fujii K, Nakajima T, Kato T, Takeda N, Okuda J, Ida K, Kawahito Y, Yoshikawa T, Okanoue T: The severity of ultrasonographic findings in nonalcoholic fatty liver disease reflects the metabolic syndrome and visceral fat accumulation. Am J Gastroenterol. 2007, 102: 2708-2715. 10.1111/j.1572-0241.2007.01526.x.

  26. 26.

    Cheung CL, Cheung TT, Lam KS, Cheung BM: High ferritin and low transferrin saturation are associated with pre-diabetes among a national representative sample of U.S. adults. Clin Nutr. 2013, 32: 1055-1060. 10.1016/j.clnu.2012.11.024.

  27. 27.

    Angulo P, Hui JM, Marchesini G, Bugianesi E, George J, Farrell GC, Enders F, Saksena S, Burt AD, Bida JP, Lindor K, Sanderson SO, Lenzi M, Adams LA, Kench J, Therneau TM, Day CP: The NAFLD fibrosis score: a noninvasive system that identifies liver fibrosis in patients with NAFLD. Hepatology. 2007, 45: 846-854. 10.1002/hep.21496.

  28. 28.

    Browning JD, Szczepaniak LS, Dobbins R, Nuremberg P, Horton JD, Cohen JC, Grundy SM, Hobbs HH: Prevalence of hepatic steatosis in an urban population in the United States: impact of ethnicity. Hepatology. 2004, 40: 1387-1395. 10.1002/hep.20466.

  29. 29.

    Lazo M, Hernaez R, Eberhardt MS, Bonekamp S, Kamel I, Guallar E, Koteish A, Brancati FL, Clark JM: Prevalence of nonalcoholic fatty liver disease in the United States: the Third National Health and Nutrition Examination Survey, 1988-1994. Am J Epidemiol. 2013, 178: 38-45. 10.1093/aje/kws448.

  30. 30.

    Ryan CK, Johnson LA, Germin BI, Marcos A: One hundred consecutive hepatic biopsies in the workup of living donors for right lobe liver transplantation. Liver Transpl. 2002, 8: 1114-1122. 10.1053/jlts.2002.36740.

  31. 31.

    Saadeh S, Younossi ZM, Remer EM, Gramlich T, Ong JP, Hurley M, Mullen KD, Cooper JN, Sheridan MJ: The utility of radiological imaging in nonalcoholic fatty liver disease. Gastroenterology. 2002, 123: 745-750. 10.1053/gast.2002.35354.

  32. 32.

    Ratziu V, Bellentani S, Cortez-Pinto H, Day C, Marchesini G: A position statement on NAFLD/NASH based on the EASL 2009 special conference. J Hepatol. 2010, 53: 372-384. 10.1016/j.jhep.2010.04.008.

Download references

Acknowledgements

BMYC received support from the Sun Chieh Yeh Heart Foundation and CLC was supported by the Faculty Development Fund (Li Ka Shing Faculty of Medicine, University of Hong Kong), and the PDF/RAP Scheme (University of Hong Kong).

Author information

Correspondence to Ching-Lung Cheung or Bernard MY Cheung.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cheung, C., Lam, K.S., Wong, I.C. et al. Non-invasive score identifies ultrasonography-diagnosed non-alcoholic fatty liver disease and predicts mortality in the USA. BMC Med 12, 154 (2014) doi:10.1186/s12916-014-0154-x

Download citation

Keywords

  • NAFLD
  • NAFLD liver fat score
  • Mortality