If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Nonalcoholic fatty liver disease (NAFLD) is the most common form of liver disease worldwide. There are limited biomarkers that can detect progression from simple steatosis to nonalcoholic steatohepatitis (NASH). The purpose of our study was to utilize CT texture analysis to distinguish steatosis from NASH.
Methods
16 patients with NAFLD (38% male, median (interquartile range): age 57 (48–64) years, BMI 37.5 (35.0–46.8) kg/m2) underwent liver biopsy and abdominal non-contrast CT. CT texture analysis was performed to quantify gray-level tissue summaries (e.g., entropy, kurtosis, skewness, and attenuation) using commercially available software (TexRad, Cambridge England). Logistic regression analyses were performed to quantify the association between steatosis/NASH status and CT texture. ROC curve analysis was performed to determine sensitivity, specificity, AUC, 95% CIs, and cutoff values of texture parameters to differentiate steatosis from NASH.
Results
By histology, 6/16 (37%) of patients had simple steatosis and 10/16 (63%) had NASH. Patients with NASH had lower entropy (median, interquartile range (IQR): 4.3 (4.1, 4.8) vs. 5.0 (4.9, 5.2), P = 0.013) and lower mean value of positive pixels (MPP) (34.4 (21.8, 52.2) vs. 66.5 (57.0, 70.7), P = 0.009) than those with simple steatosis. Entropy values below 4.73 predict NASH with 100% (95%CI: 67–100%) specificity and 80% (50–100%) sensitivity, AUC: 0.88. MPP values below 54.0 predict NASH with 100% (67–100%) specificity and 100% (50–100%) sensitivity, AUC 0.90.
Conclusion
Our study provides preliminary evidence that CT texture analysis may serve as a novel imaging biomarker for disease activity in NAFLD and the discrimination of steatosis and NASH.
NAFLD encompasses a spectrum ranging from simple steatosis to nonalcoholic steatohepatitis (NASH). NASH can lead to cirrhosis, decompensated liver disease and hepatocellular carcinoma
and is the second leading indication for liver transplantation in the United States and predicted to become the leading indication within the next decade.
Clinical research in NAFLD and NASH has been limited by the lack of validated non-invasive techniques to assess liver steatosis, inflammation and fibrosis. Liver biopsy remains the gold standard for diagnosis of NAFLD, fibrosis and NASH and is unique in its ability to quantify lobular inflammation and hepatocyte ballooning and exclude secondary causes of liver disease. However, biopsy is invasive and associated with complications such as bleeding and infection,
Thus alternative, non-invasive approaches to liver disease staging are of high priority. Magnetic resonance imaging (MRI) and proton spectroscopy (1H-MRS)
Specificity of unenhanced CT for non-invasive diagnosis of hepatic steatosis: implications for the investigation of the natural history of incidental steatosis.
can predict fibrosis in NAFLD but require access to liver MRI and specialized sequences and post-processing. Transient elastography (FibroScan) is a technique that has good predictive accuracy for higher levels of fibrosis and cirrhosis, but it is operator dependent, requires specialized equipment and is limited in patients with obesity.
CT texture analysis is an objective approach to quantify gray-level patterns of tissues by using computer-assisted measurements that are independent of subjective visual interpretation.
An advantage of CT texture analysis over advanced MR techniques is the potential for opportunistic hepatic assessment on routine abdominal CTs performed for other purposes, which could be prospectively or retrospectively analyzed to detect assess for NAFLD/NASH.
Studies have used CT texture analysis in patients with neoplasms to differentiate benign from malignant lesions, assess tumor grade, or predict survival in patients with cancer.
Texture analysis of non-small cell lung cancer on unenhanced computed tomography: initial evidence for a relationship with tumour glucose metabolism and stage.
Recently, CT texture analysis, using in-house developed programs, has been used to assess hepatic fibrosis in chronic liver disease on contrast-enhanced and non-enhanced CTs; however, these programs are not widely available.
performed CT texture analysis using a commercially available texture analysis program but were only able to predict NASH in patients without suspicion of hepatic fibrosis based on hyaluronic acid levels. Therefore, additional data is needed that utilize commercially available CT texture programs.
The purpose of our study was to determine CT texture features of the liver using a commercially available texture analysis program to distinguish steatosis from NASH in a group of individuals with concurrent liver histology.
Methods
Cohort Selection
This is a retrospective cross-sectional study of patients enrolled in an institutional review board-approved and Health Insurance Portability and Accountability Act-compliant natural history study of patients with NAFLD.
All subjects provided written informed consent for enrollment into the parent natural history study prior to data collection. The study has been carried out in accordance with The Code of Ethics of the World Medical Association (Declaration of Helsinki). Patients in the repository who had undergone an abdominal CT were assessed for eligibility. Inclusion criteria were diagnosis of NAFLD and the availability of an abdominal non-contrast CT scan and liver biopsy within ±6 months of each other. A time interval of 6 months was chosen as liver histology is not expected to change if subjects maintain a stable weight and do not undergo procedures or therapy that affects hepatic steatosis. Therefore, weight change between the date of biopsy and date of CT scan was documented. Liver enzymes, aspartate aminotransferase (AST) and alanine aminotransferase (ALT), obtained within 4 weeks of the liver biopsy, were recorded. Exclusion criteria were weight loss surgery or loss of more than 1.5 kg/m2 BMI between the CT and biopsy. Additional exclusion criteria included excess alcohol intake (more than two alcoholic drinks per day in women or three in men) or possible drug induced NAFLD with history of chronic oral steroid, methotrexate or tamoxifen use (Figure 1).
Liver biopsies were reviewed by a single blinded hepatopathologist and scored for steatosis grade (0 = < 5%; 1 = 5%–33%; 2 = 33%–66%; 3 > 66%), lobular inflammation per 200 × (0 = no foci; 1 = < 2 foci; 2 = 2–4 foci; 3 = > 4 foci) and hepatocyte ballooning (0 = no ballooning; 1 = few; 2 = many) per Kleiner et al. NAFLD activity score (NAS) was then assigned from 0 to 8 as a sum of the steatosis grade, lobular inflammation and hepatocyte ballooning values.
Subjects were then classified as having simple steatosis (n = 6), defined by grade 1 or higher steatosis without ballooning, inflammation or fibrosis, or NASH (n = 10), defined by grade ≥1 in steatosis, lobular inflammation and hepatocyte ballooning per the American Association for the Study of Liver Disease (AASLD).
Non-contrast CTs of the abdomen were performed using 64-slice multidetector CT scanners (Light-Speed VCT; GE Medical Systems, Milwaukee, Wisconsin or Somatom Definition; Siemens Medical Solutions, Ann Arbor, Michigan). The following parameters were used: variable tube current (150 mA maximum depending on patient size), 120 kVp, pitch of 1–1.5, rotation time of 0.5–0.8 s, section thickness of 2.5 mm, and standard reconstruction algorithm.
CT Texture Analysis
CT texture analysis was performed by a research instructor with 6 years of experience (A.T.), who was supervised by a fellowship-trained radiologist with 15 years of experience (M.A.B.) using commercially available software (TexRad, Cambridge England). A single circular region of interest in the peripheral portion of the right hepatic lobe was drawn avoiding hepatic vessels. The texture analysis approach includes an initial filtration step to extract and enhance features of different sizes (i.e. fine, medium and coarse texture scales corresponding to 1–4 mm in size) and to reduce noise. The following texture parameters were obtained on 12-bit images (4096 gray-levels); skewness, the asymmetry of gray-level pixel distribution; kurtosis, the pointedness or peakedness of pixel distribution; entropy, the inhomogeneity of pixel distribution; and mean value of positive pixels (MPP), the average attenuation value of pixels or greater than zero without filter application (Supplementary Figure).
CT texture analysis is not currently an FDA-approved technique.
Statistical Analysis
Descriptive summaries were computed by disease group (NASH, steatosis). Categorical variables were summarized as frequencies and percentages, and group differences were assessed using Fisher's exact test. Continuous variables were summarized as median (interquartile range (IQR) = 25th and 75th percentiles), and group differences were assessed using Wilcoxon rank sum test. Box-plots were created to summarize the distribution of each texture variable by disease group.
Separate univariate logistic regression models were constructed to quantify the association between NASH/steatosis status and each texture value (entropy, kurtosis, MPP, and skewness). Both linear and non-linear representations of each texture value was considered, but there was insufficient evidence to justify the more flexible coding of the texture value (via a chunk test). Exponentiated parameter estimates, and their 95% confidence intervals, were calculated to summarize the odds ratio (OR) of having the diagnosis of NASH for each standard-deviation increase in texture score. A receiver operator characteristic (ROC) curve analysis was performed, and the associated area under the ROC curve (AUC) calculated, to quantify the discriminative ability of each texture variable to diagnosis NASH. Estimates, and bootstrapped confidence intervals of sensitivity, specificity, and accuracy were computed for cut-points (or thresholds) along the ROC curve. The “best” cut-point was defined using Youden's Index (the value that maximizes the average sensitivity and specificity). All analyses were performed using R 4.1.0.
We included 16 subjects, 10 women, 6 men, median (IQR) age 57 (48–64) years and BMI 38 (35–47) kg/m2. By histology, 6/16 (37%) of subjects had simple steatosis and 10/16 (63%) patients had NASH. Of the subjects with NASH, 8/10 (80%) had fibrosis ranging from stage 1a to 4. Differences in the distributions of age were not detected by disease group [median (IQR): 58 (51–64) vs. 50 (27–61) years, P = 0.28] but patients with NASH had a significantly lower BMI than those with steatosis only [median (IQR): BMI 36 (34–38) vs. 47 (42–50), P = 0.03]. The median time between biopsy and CT was 8 days (IQR 1–57 days). There were no significant differences in liver enzymes between the groups (Table 1).
Table 1Patient Demographics and CT Texture Variables by Histology Status.
Steatosis n = 6
NASH n = 10
P
Demographics
Age (in years)
50 (27–61)
58 (51–64)
0.28
BMI (in kg/m2)
47 (42–50)
36 (34–38)
0.03
Female, % (number of subjects)
50% (3)
70% (7)
0.79
AST (in U/L)
35 (8–63)
55 (36–73)
0.23
ALT (in U/L)
40 (17–77)
72 (47–97)
0.14
CT texture features
Entropy
5.01 (4.92–5.17)
4.31 (4.08–4.82)
0.013
Kurtosis
−0.18 (−0.21–−0.12)
−0.08 (−0.15–−0.04)
0.08
Mean of positive pixels (MPP)
66.5 (57.0–70.7)
34.4 (21.8–52.2)
0.009
Skewness
0.15 (0.04–0.24)
0.09 (0.03–0.14)
0.38
Categorical variables are presented as frequencies and percentages and continuous variables are presented as median and interquartile range (IQR).
Abbreviations: BMI, body mass index; AST, Aspartate transaminase; ALT, alanine transaminase; MPP, Mean of positive pixels. P-values are derived from the Wilcoxon rank sum and the Fisher's Exact test for continuous and categorical variables, respectively.
Patients with NASH had lower mean entropy [median (IQR): 4.31 (4.08–4.82) vs. 5.01 (4.92–5.17), P = 0.013] and lower MPP [median (IQR): 34.4 (21.8–52.2) vs. 66.5 (57.0–70.7), P = 0.009] compared to those with simple steatosis. Kurtosis and mean skewness were not significantly different between simple steatosis and NASH groups (P = 0.08 and P = 0.38, respectively). There was insufficient evidence to conclude the mean skewness distributions differed between disease groups [median (IQR): 0.09 (0.03–0.14) vs. 0.15 (0.04–0.24), P = 0.38] (Table 1).
Discriminative Ability of Texture Metrics to Differentiate Steatosis from NASH
ROC curve analysis was used to determine threshold values of texture parameters to distinguish steatosis from NASH. Using a threshold of 4.73 (95% CI 4.54–4.99) for entropy resulted in a sensitivity, specificity, and accuracy of 100% (95% CI 0.67–1.00), 80% (95% CI 0.50, 1.00), and 88% (95% CI 0.69, 1.00), respectively, and an AUC of 0.88 (95% CI: 0.72–1.00) (Figure 1 and Table 2). For MPP, using a threshold of 54, the sensitivity, specificity, and accuracy were 100% (95% CI 0.50–1.00), 100% (95% CI 0.50–1.00), and 88% (95% CI 0.69–1.00), and an AUC of 0.90 (95% CI 0.74–1.00) (Figure 2 and Table 2). Accuracy and AUC for the other texture parameters are summarized in Figure 2 and Table 2.
Table 2Preliminary Diagnostic Summaries by Texture Metric and Threshold Value.
Texture Parameter
Threshold Estimate (95% CI)
Specificity Estimate (95% CI)
Sensitivity Estimate (95% CI)
Accuracy Estimate (95% CI)
Entropy
4.73 (4.54, 4.99)
1.00 (0.67, 1.00)
0.80 (0.50, 1.00)
0.88 (0.69, 1.00)
Kurtosis
−0.10 (−0.17, −0.08)
1.00 (0.67, 1.00)
0.70 (0.40, 1.00)
0.81 (0.62, 0.94)
MPP
54.03 (39.69, 63.66)
1.00 (0.50, 1.00)
1.00 (0.50, 1.00)
0.88 (0.69, 1.00)
Skewness
0.17 (0.01, 0.35)
0.67 (0.17, 1.00)
0.90 (0.50, 1.00)
0.75 (0.56, 0.94)
Abbreviations: Mean of Positive Pixels (MPP). Confidence intervals were calculated using a stratified bootstrap resampling procedure. Bolded entries reflect the “best” threshold using the Youden Index.
Our study shows that patients with NASH have lower mean entropy and MPP values compared to patients with simple steatosis, and that entropy and MPP cutoffs could potentially diagnose NASH with high accuracy. These findings suggest that CT texture analysis might serve as a novel imaging biomarker for disease activity that could noninvasively distinguish individuals with NASH from those with steatosis only. The importance of this findings is that the ability to identify high-risk individuals from the large number of patients with simple steatosis using non-invasive techniques, could impact management including close monitoring or intervention of patients with NASH.
NAFLD is a significant and growing public health threat as the progression to NASH can lead to cirrhosis, decompensated liver disease and hepatocellular carcinoma.
There are currently no FDA-approved therapies for NAFLD and NASH, partially due to the lack of validated non-invasive techniques to assess hepatic steatosis, inflammation and fibrosis. In addition, easily accessible, noninvasive methods to assess hepatic inflammation and fibrosis would be useful in clinical care of patients with NAFLD to predict prognosis and to assess the response to weight loss.
Specificity of unenhanced CT for non-invasive diagnosis of hepatic steatosis: implications for the investigation of the natural history of incidental steatosis.
can quantify hepatic steatosis, non-invasive techniques that can detect the progression from hepatic steatosis to NASH are needed for identification of high-risk individuals. Transient elastography (FibroScan) is an ultrasound-based technique that has been utilized to assess for risk of severe fibrosis.
However, FibroScan requires specialized equipment and can be operator dependent and is less accurate in patients with obesity. Advanced liver MR techniques, such as MR elastography or LiverMultiScan can assess fibrosis or combined inflammation and fibrosis in NAFLD, but these techniques are often more costly and require specialized equipment and processing.
Abdominal CTs are frequently performed in clinical practice and our goal was to utilize commercially available CT texture software to analyze non-contrast CT scans that are easily accessible and often obtained in the workup of other conditions. In addition, non-contrast CT has been found to be accurate in assessing hepatic steatosis across the weight spectrum, including in individuals with severe obesity.
CT texture analysis is a technique that determines the spatial arrangement of pixel intensities using statistical models of image intensities. CT texture analysis of the liver has been used in patients with neoplasms to predict outcome.
Texture features of colorectal liver metastases on pretreatment contrast-enhanced CT may predict response and prognosis in patients treated with bevacizumab-containing chemotherapy: a pilot study including comparison with standard chemotherapy.
utilized an in-house CT texture analysis program to evaluate and demonstrated a fair ability to differentiate between lower and higher levels of fibrosis as compared to gold standard histopathology in subjects with all types of liver disease. However, the analysis program used in the study is not commercially available, limiting its use.
A second study compared in-house developed CT texture analysis to gold standard histopathology in 149 subjects with all types of liver disease. Fifteen texture features were analyzed averaged over 10 regions of interest per case and 4 to 7 texture features were deemed optimal to classify fibrosis with a reported accuracy rate of 67% for determination of cirrhosis versus less severe stages of hepatic fibrosis.
Naganawa et al. examined a group of 88 patients with NAFLD specifically and found that liver texture parameters were better at predicting NASH in patients with low suspicion for fibrosis, based on normal hyaluronic acid levels, as compared to patients with a high suspicion for fibrosis based on elevated hyaluronic acid levels.
Standard texture parameters alone were poor predictors of NASH as defined by gold standard histopathology in this study.
Our study is unique in that we have an extremely well phenotyped NAFLD population with strict limitations for weight loss or weight loss surgery between the time of liver biopsy and CT. Liver biopsies were read by a single, blinded hepatopathologist. We only analyzed non-enhanced CT scans and utilized a commercially available program, TexRad, that could be uniformly utilized across different institutions. We also uniformly analyzed one filter for analysis, which is a more valid approach than utilizing data from multiple filters, as performed in prior studies.
Thus, the suggestion that entropy and MPP might serve as a novel imaging biomarker for NASH versus simple steatosis from this smaller cohort is promising but requires validation in larger longitudinal studies.
The major limitation of our study includes the small sample size, which was limited due to the exclusion of subjects with weight loss surgery or more than an absolute BMI change of more than 1.5 kg/m2 between the time of the CT and biopsy. Another limitation includes the timeframe of 6 months between CT and liver biopsy, however, the median time between CT and biopsy was 8 days. We have chosen the time interval of 6 months as liver histology is not expected to change if subjects maintain a stable weight, and in addition excluded subjects who demonstrated a BMI change of more than 1.5 kg/m2 between the CT and biopsy.
In conclusion, unenhanced CT texture analysis might be used as a novel imaging biomarker for disease progression in NAFLD, particularly to identify individuals with NASH versus steatosis only or those with fibrosis versus no fibrosis. Additional research should be conducted on larger cohorts explore the potential of this technology for diagnosis and monitoring of NASH and/or fibrosis.
CREDIT AUTHORSHIP CONTRIBUTION STATEMENT
Made a substantial contribution to the concept or design of the work; or acquisition, analysis or interpretation of data: L.D., A. T., N.M., K. C., J.H., S. O., M. C., E. R., K. M., M. B.
Drafted the article or revised it critically for important intellectual content: L.D., A. T., J.H., E. R., K. M., M. B.
Approved the version to be published: L.D., A. T., N.M., K. C., J.H., S. O., M. C., E. R., K. M., M. B.
All authors have participated sufficiently in the work to take public responsibility for appropriate portions of the content.
Conflicts of interest
LED has received research support from Pfizer (drug donation), Perspectum Ltd. (scan analysis donation) and Lumos Pharma per investigator-initiated requests and has research support from Recordati. Dr. Dichtel is a Mass General Brigham Innovation Fellow hosted by Third Rock Ventures, a venture capital firm. She remains full-time at MGH during the period of this educational program (anticipated 10/1/2022–9/30/2024). Dr. Dichtel's financial interests were reviewed and are managed by MGH and Mass General Brigham in accordance with their conflict-of-interest policies. KEC serves on the scientific advisory board for Novo Nordisk, Theratechnologies and Bristol Myers Squibb (BMS) and has received grant funding from Boehringer-Ingelheim, BMS and Novartis. KKM has received investigator-initiated grants and study medication from Amgen and study drug from Pfizer and has or recently had equity in Amgen, Bristol-Myers Squibb, General Electric, Boston Scientific, and Becton Dickinson.
Acknowledgments
This work was supported by NIH K23 DK113220 (Dichtel), K24 HL092902 (Miller), and K24 DK109940 (Bredella).
Funding
This work was supported by NIH K23 DK113220 (Dichtel), K24 HL092902 (Miller), and K24 DK109940 (Bredella).
Supplementary data
The following is the Supplementary data to this article.
Supplementary FigureLiver CT texture analysis of patient with NASH (dark blue) and patient with steatosis and no histologic evidence of NASH (turquoise). The patient with NASH has higher kurtosis with higher peak of histogram, (−0.08 vs. −0.65) and lower skewness, indicated by right shift of histogram toward higher CT attenuation (−0.09 vs. 0.03), compared to the patient without NASH.
Specificity of unenhanced CT for non-invasive diagnosis of hepatic steatosis: implications for the investigation of the natural history of incidental steatosis.
Texture analysis of non-small cell lung cancer on unenhanced computed tomography: initial evidence for a relationship with tumour glucose metabolism and stage.
Texture features of colorectal liver metastases on pretreatment contrast-enhanced CT may predict response and prognosis in patients treated with bevacizumab-containing chemotherapy: a pilot study including comparison with standard chemotherapy.