Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 26 March 2021
Sec. Cancer Imaging and Image-directed Interventions
This article is part of the Research Topic Ultrasound in Oncology: Application of Big Data and Artificial Intelligence View all 15 articles

Machine Learning-Based Ultrasomics Improves the Diagnostic Performance in Differentiating Focal Nodular Hyperplasia and Atypical Hepatocellular Carcinoma

Wei Li&#x;Wei Li1†Xiao-Zhou Lv&#x;Xiao-Zhou Lv2†Xin ZhengXin Zheng1Si-Min RuanSi-Min Ruan1Hang-Tong HuHang-Tong Hu1Li-Da ChenLi-Da Chen1Yang HuangYang Huang1Xin LiXin Li3Chu-Qing ZhangChu-Qing Zhang4Xiao-Yan XieXiao-Yan Xie1Ming Kuang,Ming Kuang1,5Ming-De Lu,Ming-De Lu1,5Bo-Wen Zhuang*Bo-Wen Zhuang1*Wei Wang*Wei Wang1*
  • 1Department of Medical Ultrasonics, Institute of Diagnostic and Interventional Ultrasound, Ultrasomics Artificial Intelligence X-Lab, The First Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China
  • 2Department of Traditional Chinese Medicine, The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, China
  • 3Research Center, GE Healthcare, Shanghai, China
  • 4Zhongshan School of Medicine, Sun Yat-sen University, Guangzhou, China
  • 5Department of Hepatobiliary Surgery, The First Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China

Background: The typical enhancement patterns of hepatocellular carcinoma (HCC) on contrast-enhanced ultrasound (CEUS) are hyper-enhanced in the arterial phase and washed out during the portal venous and late phases. However, atypical variations make a differential diagnosis both challenging and crucial. We aimed to investigate whether machine learning-based ultrasonic signatures derived from CEUS images could improve the diagnostic performance in differentiating focal nodular hyperplasia (FNH) and atypical hepatocellular carcinoma (aHCC).

Patients and Methods: A total of 226 focal liver lesions, including 107 aHCC and 119 FNH lesions, examined by CEUS were reviewed retrospectively. For machine learning-based ultrasomics, 3,132 features were extracted from the images of the baseline, arterial, and portal phases. An ultrasomics signature was generated by a machine learning model. The predictive model was constructed using the support vector machine method trained with the following groups: ultrasomics features, radiologist’s score, and combination of ultrasomics features and radiologist’s score. The diagnostic performance was explored using the area under the receiver operating characteristic curve (AUC).

Results: A total of 14 ultrasomics features were chosen to build an ultrasomics model, and they presented good performance in differentiating FNH and aHCC with an AUC of 0.86 (95% confidence interval [CI]: 0.80, 0.89), a sensitivity of 76.6% (95% CI: 67.5%, 84.3%), and a specificity of 80.5% (95% CI: 70.6%, 85.9%). The model trained with a combination of ultrasomics features and the radiologist’s score achieved a significantly higher AUC (0.93, 95% CI: 0.89, 0.96) than that trained with the radiologist’s score (AUC: 0.84, 95% CI: 0.79, 0.89, P < 0.001). For the sub-group of HCC with normal AFP value, the model trained with a combination of ultrasomics features, and the radiologist’s score remain achieved the highest AUC of 0.92 (95% CI: 0.87, 0.96) compared to that with the ultrasomics features (AUC: 0.86, 95% CI: 0.74, 0.89, P < 0.001) and radiologist’s score (AUC: 0.86, 95% CI: 0.79, 0.91, P < 0.001).

Conclusions: Machine learning-based ultrasomics performs as well as the staff radiologist in predicting the differential diagnosis of FNH and aHCC. Incorporating an ultrasomics signature into the radiologist’s score improves the diagnostic performance in differentiating FNH and aHCC.

Introduction

The typical enhancement pattern of hepatocellular carcinoma (HCC) on contrast-enhanced ultrasound (CEUS) is characterized by hyper-enhancement in the arterial phase and wash out during the portal venous and late phases (1). However, atypical variations occur, especially in some well-differentiated tumors, accounting for 5–41% of HCC cases; such lesions may show sustained hyper-/iso-enhancement in the portal venous and late phases and are defined as atypical HCC (aHCC) (24). Meanwhile, most benign focal liver lesions show complete hyper- or iso-enhancement in the portal venous and late phases, making differential diagnosis both crucial and challenging (5, 6). This diagnostic difficulty could be resolved using CEUS techniques, such as micro-flow imaging to further characterize the enhancement features in the arterial phase, e.g., a spoke-wheel artery for focal nodular hyperplasia (FNH) and chaotic vessel for HCC (710). However, the interpretation of features involves the experience of radiologists, making inter-reader variability inevitable.

In contrast to the traditional practice of treating medical images as pictures intended solely for visual interpretation, radiomics features could reflect not only the macroscopic manifestation but also the cellular and molecular nature of tissues (1113). Radiomics offers a vast scale of imaging biomarkers that could potentially assist in detecting and diagnosing, evaluating the prognosis and predicting the therapeutic response, and monitoring the disease status of cancer (11, 12, 1416). Machine learning-based ultrasomics approaches, derived from radiomics, involve the analysis and transformation of ultrasound images into large sets of quantitative data and have been identified as potential alternatives to detect and classify lesions (17, 18).

Recently, few applications of machine learning in HCC diagnosis have been reported (19, 20). Most machine learning systems have demonstrated excellent diagnostic performance, with the area under the receiver operating characteristic curve (AUC) of 0.89-0.97 for HCC characterization (19, 20). Gatos et al. applied radiomics to segment and classify focal liver lesions on non-enhanced T2-weighted images, providing a non-invasive method for assessing liver lesions (21). Some studies have shown that multi-modal ultrasound images also perform well for the detection and classification of focal liver lesions (19, 22, 23). However, most studies have only compared the diagnostic performance between machine learning systems and radiologists. The influence of the performance of these systems on radiologists when used in clinical practice has not been evaluated. Thus, the added clinical value of machine learning systems to observers is necessary to determine and validate.

The purpose of our research was to develop a machine learning-based ultrasomics approach to assess ultrasomics features for improving the diagnostic performance in differentiating FNH and aHCC.

Patients and Methods

Patients

This retrospective analysis obtained ethical approval and waived the informed consent requirement. From December 2013 to January 2018, 119 patients with FNH and 107 patients with aHCC lesions were included in the study based on the inclusion and exclusion criteria. The inclusion criteria were as follows: (a) CEUS was performed; (b) lesions were visually hyper-enhanced during the arterial phase and sustained hyper- or iso-enhanced during the portal venous and late phases; (c) HCC was diagnosed by pathological examinations and FNH was confirmed by pathological examinations or supported by CT or MRI findings with a minimum 1 year follow-up; and (d) no treatment was conducted before CEUS. Patients were excluded if they had multiple tumors. Baseline clinical trial data, including age, gender, and some blood test, such as hepatitis background and alpha-fetoprotein (AFP), were performed no more than 7 days before or after the CEUS examination.

Image Acquisition

US examinations were performed using an Aplio 500 scanner (Canon Medical Systems, Tokyo, Japan), equipped with a 375BT convex transducer (frequency, 3.5 MHz) and an Aixplorer scanner (Supersonic, Paris, France) with an SC6-1 curvilinear transducer (frequency, 1–6 MHz). Contrast harmonic imaging (CHI) and contrast pulse sequencing (CPS) were used with a mechanical index of 0.06–0.10. Baseline ultrasonography was performed to scan the liver thoroughly before CEUS. Additionally, the target lesions were identified and observed carefully during the baseline observation in B mode. The imaging settings, such as the gain, depth, and focus, were optimized for each examination. After the CHI or CPS mode was activated, a bolus intravenous injection of 2.4 mL of SonoVue (Bracco, Milan, Italy) was administered, followed by flushing with 5 mL of saline. The targeted lesion was observed continuously for 5 minutes. The arterial, portal venous, and late phases were defined as 10–30 seconds, 31–120 seconds, and 121–300 seconds after injection, respectively. CEUS examinations were performed by one of two radiologists (WW and X-YX) with at least 10 years’ experience of performing CEUS. Three images from the same section, which showed the maximum observation of the target lesion, were taken from each patient, (a) a baseline ultrasound; (b) an arterial phase image at the enhancement time of 25–30 seconds; and (c) a portal phase image at the enhancement time of 60–70 seconds.

Radiologist’s Scoring

Two staff radiologists (B-WZ and L-DC) reviewed ultrasound images and videos retrospectively, and they had more than 5 years of experience in assessing liver CEUS data. The radiologists were not involved in the feature extraction process below. All patient identification information from the images was removed, and the researchers were unaware of all the clinicopathological information. The diagnostic criteria for HCC and FNH were based on the 2012 guidelines issued by the European Federation of Societies for Ultrasound in Medicine and Biology (EFSUMB) (2). The diagnostic criteria for HCC were the manifestation of basket pattern and/or chaotic vessels (Figure 1A) and non-enhanced areas (Figure 1B). The diagnostic criteria of CEUS feature for FNH were centrifugal enhancement (Video 1), spoke-wheel artery (Figure 1C), unenhanced central scar (Figure 1D), and feeding artery.

FIGURE 1
www.frontiersin.org

Figure 1 Typical features for HCC and FNH lesions. (A) the basket pattern and/or chaotic vessels; (B) non-enhancing areas (arrow); (C) spoke-wheel arteries; and (D) unenhanced central scar (arrow). Annotations of the ROI generated by the radiologists around the tumor outline are delineated in red.

For each group, the diagnostic confidence was scored using a subjective three-point scale (grade 1, definitely or most likely FNH; grade 2, indeterminate; and grade 3, most likely or definitely HCC). If there was inconsistency, we performed a consensus reading, and the consensus data were used for subsequent analysis.

Ultrasomics Feature Extraction

Digital imaging and communications in medicine (DICOM) images were used to extract ultrasomics features using the in-house designed Ultrasomics-Platform software (Version 1.0; Ultrasomics Artificial Intelligence X-lab, Guangzhou, China). After an image was imported, the radiologist drew a region of interest (ROI) on the largest cross-section along the tumor contour. Next, the software automatically extracted the features from the ROIs. In total, 1,044 features could be extracted from a single image. These 1,044 features extracted from a single image consisted of five categories of features: histogram parameters, textural parameters, form factor parameters, grey-level co-occurrence matrix (GLCM) parameters, and run length matrix (RLM) parameters. Detailed information on the features is provided in Supplementary Material S1. Finally, 3,132 features were extracted from the baseline US, arterial phase and portal phase of CEUS images of each patient. Initially, two radiologists (WL and YH, with at least 5 years of experience in performing US examinations) were required to trace out ROIs on the selected images. The inter-observer and intra-observer reproducibility in feature extraction were assessed and are described in Supplementary Material S2. The remaining images were delineated by the first radiologist.

Feature Selection and Model Development for Prediction

Of the 3,132 features from each patient, many were highly redundant, which could degrade the classification. We eliminated redundant features by using a two-step feature selection method. First, if two features were highly-correlated with a correlation coefficient higher than 0.95, one of the features was removed. Second, we eliminated features with an AUC less than 0.6. According to the Harrell guidelines for multivariate analysis, the number of events should be at least 10 times greater than that of the included covariates (24). The least absolute shrinkage and selection operator (LASSO) regression was used to perform the ultrasomics features selection in the training dataset. All ultrasomics feature values were normalized by using the mean and variance of the feature values to be within similar dynamic ranges.

A support vector machine (SVM) based on the radial basis function (RBF) kernel was trained from the selected feature subset produced by the preceding steps. The entire data set was randomly divided into a training dataset (comprising 80% of subjects) and a validation dataset (comprising the remaining 20% of subjects). The training dataset was used to construct a model, which was then evaluated using the validation dataset. A 10-fold cross-validation method was adopted to ensure the robustness of the classifiers to training and testing data. All processes were repeated 10 times with random seeds, generating 10 different training and validation datasets. We built the model using the training dataset and then evaluated it using the validation dataset repeatedly. Subsequently, the model with the best classification performance was selected as the best model.

Statistical Analysis

Descriptive statistics are summarized as the mean ± standard deviation (SD) or median and interquartile range. Comparisons between groups were tested using Student’s t test or the Mann-Whitney test for quantitative variables and the chi-squared test or Fisher’s test for qualitative variables.

A weighted kappa statistics test was used to assess the two radiologists’ scores. We evaluated the reproducibility of the ultrasomics feature extraction using the “irr” package in R. The LASSO regression was performed using the “glmnet” package.

All observations of patients with known outcomes were classified into three datasets (1): radiologist’s score, (2) ultrasomics features, and (3) a combination of ultrasomics features and radiologist’s score. The diagnostic performance of the radiologist’s score was evaluated by plotting receiver operating characteristic (ROC) curves. The diagnostic performance in discriminating between FNH and aHCC is expressed as the AUC. The ultrasomics features and the combination of the ultrasomics features and radiologist’s score were further compared through an SVM classifier using the “rattle” package in R. The performance of the SVM model was tested using the AUC. Paired comparisons of AUC values were performed by a two-sided Wilcoxon signed-rank test at a significance level of 5%. The predictive sensitivity (SEN), specificity (SPE), positive predictive value (PPV), negative predictive value (NPV), positive likelihood ratio (+LR), and negative likelihood ratio (-LR) were calculated at a cut-off point that maximized the value of the Youden index. Comparisons among the three datasets were performed using the Delong test. Decision curve analysis (DCA) was performed with the “dca.R” function. All statistical tests were two-sided tests, and P< 0.05 indicated statistical significance. All Statistical analyses were performed using R version 3.3.3 (http://www.r-project.org/).

Results

Clinical Characteristics

The clinical characteristics are listed in Table 1. The study included 226 patients; 107 (47.3%) patients (mean age, 54.0 ± 11.9 years old) had a final diagnosis of HCC; and the remaining 119 (52.7%) patients (mean age, 34.5 ± 11.7 years old) had a final diagnosis of FNH. 20 FNH lesion were confirmed by pathological examinations (11 by biopsy, 9 by surgery), while 99 cases were supported by CT or MRI findings with a minimum one-year follow-up. No significant difference was found in the tumor number between the two groups (P=0.118). The average lesion size of FNH and HCC was 3.3 ± 1.8 cm (range: 0.8-10.2 cm) and 4.8 ± 3.4 cm (range: 0.8-18.6 cm), respectively.

TABLE 1
www.frontiersin.org

Table 1 Clinical Characteristics and Laboratory Information of the Patients.

Ultrasomics Signature Construction and Validation

After the feature selection and dimensional reduction process, 14 selected features were taken as the input of the SVM to train a prediction model, including 6 features derived from baseline US images, 3 from arterial phase images, and 4 from portal phase images (Figure 2, Supplementary Material S3). All feature values were normalized to achieve similar dynamic ranges. The parameter C which is used to control the error-margin trade-off was set at 1, and the kernel width sigma was 0.012. Next, the training and validation procedures for tumor classification were employed with 10-fold cross-validation.

FIGURE 2
www.frontiersin.org

Figure 2 Radiomics feature selection using the least absolute shrinkage and selection operator (LASSO) regression model. The 10-fold cross-validation (CV) process was repeated 50 times to generate the optimal penalization coefficient lambda (λ) in the LASSO model. The value of λ that produced the minimum average binomial deviance was used to select features. Dotted vertical lines were drawn at the optimal values using the minimum criteria and the 1 standard error of the minimum criteria (the 1-SE criteria). A λ value of 0.043 was chosen (the 1-SE criteria) according to 10-fold CV, where optimal λ resulted in 14 nonzero coefficients.

Diagnostic Performance of Ultrasomics Features and Radiologist’s Score Models

A total of 83 FNHs were correctly identified while 36 lesions were incorrectly identified as HCC, leading to a specificity of 69.8% by the radiologists. For combined model, 100 FNHs were correctly classified whereas 19 lesions were incorrectly assigned to HCC, resulting in a specificity of 84.0%. Comparing the performance of the radiologists' score and the combined model, twenty-four cases have a different result; consequently, the combined model leads to an additional 15 FNHs and 5 HCCs being correctly classified.

The model trained with the combination of the ultrasomics features and radiologist’s score performed significantly better (AUC: 0.93, 95% CI: 0.89, 0.96) than that trained with the ultrasomics features (AUC: 0.86, 95% CI: 0.80, 0.89, P < 0.001) and radiologist’s score (AUC: 0.84, 95% CI: 0.79, 0.89, P < 0.001). Overall, the model based on the radiologist’s score had the highest diagnostic SEN of 94.4% (95% CI: 88.2%, 97.9%) but the lowest SPE of 69.8% (95% CI: 60.7%, 77.8%) compared with the combined model (SEN: 93.5% [95% CI: 87.0%, 97.3%], SPE: 84.9% [95% CI: 77.1%, 90.8%]), and the ultrasomics features model (SEN: 76.6% [95% CI: 67.5%, 84.3%], SPE: 80.5% [95% CI: 70.6%, 85.9%]). Furthermore, when the ultrasomics features were combined with the radiologist’s score, the diagnostic performance was significantly improved in terms of the AUC, SPE, and PPV and + LR (AUC: 0.93, SPE: 84.9%, and PPV: 84.7%, +LR 6.2) compared with the performance of the other two models. The performance measurements of each dataset are reported in Table 2 based on each ROC curves to distinguish between FNH and aHCC (Figure 3).

TABLE 2
www.frontiersin.org

Table 2 Diagnostic Performance of the Three Models in Differentiating Focal Nodular Hyperplasia and Atypical Hepatocellular Carcinoma.

FIGURE 3
www.frontiersin.org

Figure 3 Receiver operating characteristic curves of the combination of ultrasomics features and radiologist’s score (blue curve), ultrasomics features (green curve), and radiologist’s score (orange curve). The areas under the curves are 0.93, 0.86, 0.84, respectively.

The DCA shows that within most reasonable threshold probability ranges, the combined model showed the highest overall net benefit than the radiologist’s score or ultrasound feature model. The DCA results for the three models are presented in Figure 4.

FIGURE 4
www.frontiersin.org

Figure 4 Decision curve analysis for each model. The y-axis measures the net benefit. The net benefit was calculated by summing the benefits (true positive results) and subtracting the harms (false-positive results), weighting the latter by a factor related to the relative harm of undetected cancer compared with the harm of unnecessary treatment. The combined model (yellow line) had the highest net benefit compared with the other two models (blue line and red line) and simple strategies, such as the follow-up of all patients (grey line) or no patients (horizontal black line), across the full range of threshold probabilities at which a patient would choose to undergo a follow-up imaging examination.

Validation in the Sub-group of HCC With Normal AFP Level

Forty patients were confirmed to develop HCC with normal AFP level. The model trained above performed significantly better (AUC: 0.92, 95% CI: 0.87, 0.96) with the combination of the ultrasomics features and radiologist’s score than that with the ultrasomics features (AUC: 0.86, 95% CI: 0.74, 0.89, P < 0.001) and radiologist’s score (AUC: 0.86, 95% CI: 0.79, 0.91, P < 0.001) (Table 3).

TABLE 3
www.frontiersin.org

Table 3 Validation in the sub-group of HCC with normal AFP value.

Discussion

In this study, we derived and validated an ultrasomics-based machine learning approach to analyze ultrasound images for the preoperative individualized diagnosis of FNH and aHCC. Our analysis reveals that the diagnostic performance of ultrasomics is comparable to that of a staff radiologist in differentiating between FNH and aHCC. Furthermore, when adding ultrasomics to the radiologist’s classification, the diagnostic performance was improved significantly with an AUC ranging from 0.84 to 0.93 (P < 0.001). Our study shows that ultrasomics may increase the diagnostic confidence of radiologists in CEUS examinations and potentially improve their accuracy when facing atypical features. Thus, clinicians would benefit from this decision-making process in the diagnosis of HCC.

In clinical practice, when radiologists face a lesion that shows hyper-enhancement in the arterial phase and sustained enhancement in the portal vein and late phases, it is difficult to diagnose HCC. However, the high sensitivity and negative predictive values of the radiologist would be useful in clinical practice for excluding disease; thus, HCC would be excluded by the radiologist if the result was considered to be FNH by ultrasomics. Therefore, this system would help to reduce unnecessary biopsies or active clinical treatment requested by experienced radiologists.

In contrast, ultrasomics, referred to as high-throughput computing, extracts innumerable quantitative features from US images (18). By transforming digital medical images into mineable high-dimensional data, ultrasomics yields features, such as textural features, that could objectively reflect the homogeneity or heterogeneity of an image. These patterns could represent enhancement features just as heterogeneity might represent chaotic vessels and necrosis. Focal liver lesions can be featured by typical features in the arterial phase and wash-out during the portal and late phases. In this study, the features displayed in the arterial phase could provide a major benefit for the diagnosis of liver tumors. In previous studies, we utilized a maximum intensity projection technique of micro-flow imaging and achieved higher spatial resolution and higher temporal resolution when detecting vessel contours. Compared to conventional CEUS features (AUC: 0.84), micro-flow imaging technology provided significant improvements over the detection rates achieved for the staff radiologists (AUC: 0.89) (10). In this study, ultrasomics features alone can achieve a similar diagnostic performance (AUC: 0.86) as micro-flow imaging (AUC: 0.868–0.873). Ultrasomics could reach such achievement because it analyzes textural features objectively and quantitatively to describe the intrinsic characteristics of tumors, in particular heterogeneous tumors. Ultrasomics analysis has already been applied to various types of disease, such as HCC, liver fibrosis, and breast cancer (18, 2527). The potential of ultrasomics has already been demonstrated for liver imaging in some studies (17, 18).

In this study, we additionally evaluated the benefit of ultrasomics in assisting doctors with the interpretation of medical images. As an interdisciplinary technology, it combines elements of imaging generation, digital image processing, statistical imaging, and knowledge engineering to manage the volume of information related to the diagnostic process and outcome prediction (17, 18, 25). In our study, the additional information in the combined model led to improved diagnostic performance (AUC: 0.93) and higher specificity of 84.9% compared with the ultrasomics (AUC: 0.86, specificity: 80.5%) and radiologist’s score (AUC: 0.84, specificity: 69.8%) models. The combined model was also comparable to the ML model based on multi-modal ultrasound images (AUC: 0.94, sensitivity: 91.0%, specificity: 86.0%) (19). The results are also comparable to and even better than those of MRI (AUC: 0.89, sensitivity: 82.2%, specificity: 71.4%), as previously reported (20). However, the use of artificial intelligence is not intended to replace expert diagnosticians because no solution is guaranteed and knowledge-based maintenance is required. Artificial intelligence is also affected by several elements, such as the source of images and the cognition of disease. Presently, most domains of large data have not tapped the full potential of artificial intelligence technology. However, rapid developments in the area will add more potential to the advantages. Therefore, the most important role of artificial intelligence is to help improve diagnostic accuracy and assist rather than replace clinicians in making treatment decisions. It is worth noting that the combined model greatly improved the diagnostic ability of radiologists. A similar conclusion was obtained in another study (28).

Our research has some limitations. First, this study was retrospective and conducted in one center. This may cause potential variations and selection bias in the patient population and imaging methods, which is difficult to generalize the outcomes to other agencies. Second, due to the relative rarity of aHCC, the sample size is relatively small, which may cause over-fitting to this particular population. Hence, large-scale multicenter studies are necessary for the future to validate the results. Third, only two radiologists were involved in the assessment of the basic imaging features and feature extraction. All outcomes were based on the features extracted by one radiologist, which may not be generalizable to all radiologists. Fourth, the machine and imaging settings in this study were inconsistent, which may affect the ultrasomics features (29, 30).

In conclusion, an ultrasomics approach was developed to investigate the association between the quantitative ultrasound features and pathological characteristics of tumors effectively and objectively. We evaluated the added value of ultrasomics to the radiologist, and this approach improved the performance of CEUS by providing quantitative and standardized criteria to radiologists, thereby enabling the more confident application of CEUS in detecting HCC to achieve better treatment planning. Our findings can assist clinicians in the differential diagnosis between FNH and aHCC accurately using CEUS images, and this allows for early and precise medical management and treatment.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of the First Affiliated Hospital of Sun Yat-Sen University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

WW, M-DL, MK, and X-YX contributed conception and design of the study. XZ collected data from clinical trials. WL and B-WZ drafted and revised the manuscript. X-ZL, WW, L-DC, and S-MR revised the manuscript for important intellectual content. X-ZL, XL, H-TH, YH, and C-QZ gave technical support. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by grants from the National Natural Science Foundation of China (No. 81601500 and 81701701).

Conflict of Interest

XL was employed by GE Healthcare.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.544979/full#supplementary-material

References

1. Martie A, Sporea I, Popescu A, Sirli R, Danila M, Serban C, et al. Contrast enhanced ultrasound for the characterization of hepatocellular carcinoma. Med Ultrason (2011) 13(2):108–13.

PubMed Abstract | Google Scholar

2. Claudon M, Dietrich CF, Choi BI, Cosgrove DO, Kudo M, Nolsoe CP, et al. Guidelines and good clinical practice recommendations for contrast enhanced ultrasound (CEUS) in the liver–update 2012: a WFUMB-EFSUMB initiative in cooperation with representatives of AFSUMB, AIUM, ASUM, FLAUS and ICUS. Ultraschall Med (2013) 34(1):11–29. doi: 10.1055/s-0032-1325499

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Nicolau C, Catala V, Vilana R, Gilabert R, Bianchi L, Sole M, et al. Evaluation of hepatocellular carcinoma using SonoVue, a second generation ultrasound contrast agent: correlation with cellular differentiation. Eur Radiol (2004) 14(6):1092–9. doi: 10.1007/s00330-004-2298-0

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Jang HJ, Kim TK, Burns PN, Wilson SR. Enhancement patterns of hepatocellular carcinoma at contrast-enhanced US: comparison with histologic differentiation. Radiology (2007) 244(3):898–906. doi: 10.1148/radiol.2443061520

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Dietrich CF. Contrast-Enhanced Ultrasound of Benign Focal Liver Lesions. Ultraschall der Med (2018) 40(1):12–29. doi: 10.1055/a-0668-5746

CrossRef Full Text | Google Scholar

6. Friedrich-Rust M, Klopffleisch T, Nierhoff J, Herrmann E, Vermehren J, Schneider MD, et al. Contrast-Enhanced Ultrasound for the differentiation of benign and malignant focal liver lesions: a meta-analysis. Liver Int (2013) 33(5):739–55. doi: 10.1111/liv.12115

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Wang W, Chen LD, Lu MD, Liu GJ, Shen SL, Xu ZF, et al. Contrast-enhanced ultrasound features of histologically proven focal nodular hyperplasia: diagnostic performance compared with contrast-enhanced CT. Eur Radiol (2013) 23(9):2546–54. doi: 10.1007/s00330-013-2849-3

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Yen YH, Wang JH, Lu SN, Chen TY, Changchien CS, Chen CH, et al. Contrast-enhanced ultrasonographic spoke-wheel sign in hepatic focal nodular hyperplasia. Eur J Radiol (2006) 60(3):439–44. doi: 10.1016/j.ejrad.2006.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Xu HX, Liu GJ, Lu MD, Xie XY, Xu ZF, Zheng YL, et al. Characterization of focal liver lesions using contrast-enhanced sonography with a low mechanical index mode and a sulfur hexafluoride-filled microbubble contrast agent. J Clin Ultrasound (2006) 34(6):261–72. doi: 10.1002/jcu.20234

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Li W, Wang W, Liu GJ, Chen LD, Wang Z, Huang Y, et al. Differentiation of Atypical Hepatocellular Carcinoma from Focal Nodular Hyperplasia: Diagnostic Performance of Contrast-enhanced US and Microflow Imaging. Radiology (2015) 275(3):870–9. doi: 10.1148/radiol.14140911

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Gillies RJ, Kinahan PE, Hricak H. Radiomics: Images Are More than Pictures, They Are Data. Radiology (2016) 278(2):563–77. doi: 10.1148/radiol.2015151169

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sun R, Limkin EJ, Vakalopoulou M, Dercle L, Champiat S, Han SR, et al. A radiomics approach to assess tumor-infiltrating CD8 cells and response to anti-PD-1 or anti-PD-L1 immunotherapy: an imaging biomarker, retrospective multicohort study. Lancet Oncol (2018) 19(9):1180–91. doi: 10.1016/S1470-2045(18)30413-3

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Banerjee S, Wang DS, Kim HJ, Sirlin CB, Chan MG, Korn RL, et al. A computed tomography radiogenomic biomarker predicts microvascular invasion and clinical outcomes in hepatocellular carcinoma. Hepatology (Baltimore Md) (2015) 62(3):792–800. doi: 10.1002/hep.27877

CrossRef Full Text | Google Scholar

14. Limkin EJ, Sun R, Dercle L, Zacharaki EI, Robert C, Reuze S, et al. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology. Ann Oncol (2017) 28(6):1191–206. doi: 10.1093/annonc/mdx034

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Huang Y, Liu Z, He L, Chen X, Pan D, Ma Z, et al. Radiomics Signature: A Potential Biomarker for the Prediction of Disease-Free Survival in Early-Stage (I or II) Non-Small Cell Lung Cancer. Radiology (2016) 281(3):947–57. doi: 10.1148/radiol.2016152234

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Prasanna P, Patel J, Partovi S, Madabhushi A, Tiwari P. Radiomic features from the peritumoral brain parenchyma on treatment-naive multi-parametric MR imaging predict long versus short-term survival in glioblastoma multiforme: Preliminary findings. Eur Radiol (2017) 27(10):4188–97. doi: 10.1007/s00330-016-4637-3

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang K, Lu X, Zhou H, Gao Y, Zheng J, Tong M, et al. Deep learning Radiomics of shear wave elastography significantly improved diagnostic performance for assessing liver fibrosis in chronic hepatitis B: a prospective multicentre study. Gut (2018) 68(4):729–41. doi: 10.1136/gutjnl-2018-316204

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Li W, Huang Y, Zhuang BW, Liu GJ, Hu HT, Li X, et al. Multiparametric ultrasomics of significant liver fibrosis: A machine learning-based analysis. Eur Radiol (2018) 29(3):1496–506. doi: 10.1007/s00330-018-5680-z

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yao Z, Dong Y, Wu G, Zhang Q, Yang D, Yu JH, et al. Preoperative diagnosis and prediction of hepatocellular carcinoma: Radiomics analysis based on multi-modal ultrasound images. BMC Cancer (2018) 18(1):1089. doi: 10.1186/s12885-018-5003-4

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wu J, Liu A, Cui J, Chen A, Song Q, Xie L. Radiomics-based classification of hepatocellular carcinoma and hepatic haemangioma on precontrast magnetic resonance images. BMC Med Imaging (2019) 19(1):23. doi: 10.1186/s12880-019-0321-9

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gatos I, Tsantis S, Karamesini M, Spiliopoulos S, Karnabatidis D, Hazle JD, et al. Focal liver lesions segmentation and classification in nonenhanced T2-weighted MRI. Med Phys (2017) 44(7):3695–705. doi: 10.1002/mp.12291

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Gatos I, Tsantis S, Spiliopoulos S, Skouroliakou A, Theotokas I, Zoumpoulis P, et al. A new automated quantification algorithm for the detection and evaluation of focal liver lesions with contrast-enhanced ultrasound. Med Phys (2015) 42(7):3948–59. doi: 10.1118/1.4921753

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Acharya UR, Koh JEW, Hagiwara Y, Tan JH, Gertych A, Vijayananthan A, et al. Automated diagnosis of focal liver lesions using bidirectional empirical mode decomposition features. Comput Biol Med (2018) 94:11–8. doi: 10.1016/j.compbiomed.2017.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Helmreich JE. Regression modeling strategies: with applications to linear models, logistic and ordinal regression, and survival analysis. 2nd edn. Switzerland: Springer International Publishing (2015). doi: 10.1007/978-3-319-19425-7

CrossRef Full Text | Google Scholar

25. Hu HT, Wang Z, Huang XW, Chen SL, Zheng X, Ruan SM, et al. Ultrasound-based radiomics score: a potential biomarker for the prediction of microvascular invasion in hepatocellular carcinoma. Eur Radiol (2019) 29(6):2890–901. doi: 10.1007/s00330-018-5797-0

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chen LD, Ruan SM, Liang JY, Yang Z, Shen SL, Huang Y, et al. Differentiation of intrahepatic cholangiocarcinoma from hepatocellular carcinoma in high-risk patients: A predictive model using contrast-enhanced ultrasound. World J Gastroenterol (2018) 24(33):3786–98. doi: 10.3748/wjg.v24.i33.3786

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Guo B, Ouyang F, Ouyang L, Liu Z, Lin S, Meng W, et al. Development and validation of an ultrasound-based nomogram to improve the diagnostic accuracy for malignant thyroid nodules. Eur Radiol (2018) 29(3):1518–26. doi: 10.1007/s00330-018-5715-5

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wang J, Wu CJ, Bao ML, Zhang J, Wang XN, Zhang YD. Machine learning-based analysis of MR radiomics can help to improve the diagnostic performance of PI-RADS v2 in clinically relevant prostate cancer. Eur Radiol (2017) 27(10):4082–90. doi: 10.1007/s00330-017-4800-5

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Yu X, Guo Y, Huang SM, Li ML, Lee WN. Beamforming effects on generalized Nakagami imaging. Phys Med Biol (2015) 60(19):7513–31. doi: 10.1088/0031-9155/60/19/7513

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Byra M, Wan L, Wong JH, Du J, Shah SB, Andre MP, et al. Quantitative Ultrasound and B-Mode Image Texture Features Correlate with Collagen and Myelin Content in Human Ulnar Nerve Fascicles. Ultrasound Med Biol (2019) 45(7):1830–40. doi: 10.1016/j.ultrasmedbio.2019.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: ultrasonography, machine learning, focal nodular hyperplasia, hepatocellular carcinoma, ultrasomics

Citation: Li W, Lv X-Z, Zheng X, Ruan S-M, Hu H-T, Chen L-D, Huang Y, Li X, Zhang C-Q, Xie X-Y, Kuang M, Lu M-D, Zhuang B-W and Wang W (2021) Machine Learning-Based Ultrasomics Improves the Diagnostic Performance in Differentiating Focal Nodular Hyperplasia and Atypical Hepatocellular Carcinoma. Front. Oncol. 11:544979. doi: 10.3389/fonc.2021.544979

Received: 23 March 2020; Accepted: 03 March 2021;
Published: 26 March 2021.

Edited by:

Youyong Kong, Southeast University, China

Reviewed by:

Yuming Jiang, Stanford University, United States
Jinghao Duan, Shandong University, China
Yan Wu, Second Affiliated Hospital of Harbin Medical University, China

Copyright © 2021 Li, Lv, Zheng, Ruan, Hu, Chen, Huang, Li, Zhang, Xie, Kuang, Lu, Zhuang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wei Wang, wangw73@mail.sysu.edu.cn; Bo-Wen Zhuang, zhuangbw3@mail.sysu.edu.cn

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.