Skip to main content

ORIGINAL RESEARCH article

Front. Med., 04 November 2021
Sec. Gastroenterology
Volume 8 - 2021 | https://doi.org/10.3389/fmed.2021.745931

CT-Based Radiomics Score Can Accurately Predict Esophageal Variceal Rebleeding in Cirrhotic Patients

Dongxiao Meng1,2 Yingnan Wei2 Xiao Feng2 Bing Kang3,4 Ximing Wang3,4 Jianni Qi5 Xinya Zhao3,4* Qiang Zhu1,2*
  • 1Department of Gastroenterology, Shandong Provincial Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
  • 2Department of Gastroenterology, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
  • 3Department of Radiology, Shandong Provincial Hospital, Cheeloo College of Medicine, Shandong University, Jinan, China
  • 4Department of Radiology, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
  • 5Department of Central Laboratory, Shandong Provincial Hospital Affiliated to Shandong University, Jinan, China

Purpose: This study aimed to develop a radiomics score (Rad-score) extracted from liver and spleen CT images in cirrhotic patients to predict the probability of esophageal variceal rebleeding.

Methods: In total, 173 cirrhotic patients were enrolled in this retrospective study. A total of 2,264 radiomics features of the liver and spleen were extracted from CT images. Least absolute shrinkage and selection operator (LASSO) Cox regression was used to select features and generate the Rad-score. Then, the Rad-score was evaluated by the concordance index (C-index), calibration curves, and decision curve analysis (DCA). Kaplan–Meier analysis was used to assess the risk stratification ability of the Rad-score.

Results: Rad-scoreLiver, Rad-scoreSpleen, and Rad-scoreLiver−Spleen were independent risk factors for EV rebleeding. The Rad-scoreLiver−Spleen, which consisted of ten features, showed good discriminative performance, with C-indexes of 0.853 [95% confidence interval (CI), 0.776–0.904] and 0.822 (95% CI, 0.749–0.875) in the training and validation cohorts, respectively. The calibration curve showed that the predicted probability of rebleeding was very close to the actual probability. DCA verified the usefulness of the Rad-scoreLiver−Spleen in clinical practice. The Rad-scoreLiver−Spleen showed good performance in stratifying patients into high-, intermediate- and low-risk groups in both the training and validation cohorts. The C-index of the Rad-scoreLiver−Spleen in the hepatitis B virus (HBV) cohort was higher than that in the non-HBV cohort.

Conclusion: The radiomics score extracted from liver and spleen CT images can predict the risk of esophageal variceal rebleeding and stratify cirrhotic patients accordingly.

Introduction

Esophageal variceal (EV) bleeding is one of the most serious complications in cirrhotic patients with portal hypertension (1). Although several recommended treatments are applied, patients who recover from the first episode of EV bleeding have a high risk of 1-year rebleeding (approximately 60%), with a mortality rate of up to 33% (2). EV rebleeding may lead to a series of complications, such as hepatic encephalopathy, spontaneous bacterial peritonitis, and liver failure, eventually making the patients lose opportunities for other remedial measures. Thus, prediction of rebleeding and the identification of patients at high risk of rebleeding after endoscopic therapy are urgent issues that could help improve the prognosis of cirrhotic patients (3).

In clinical practice, the most important predictor for variceal rebleeding is the size of the varices determined with endoscopy (4). However, the compliance of patients is affected by its expensive and invasive properties. Hepatic venous pressure gradient (HVPG) has been widely proven to be a strong predictive factor for EV bleeding in patients with cirrhosis (5), but it is available only in specialized hepatology units, which restricts its widespread use (6). Some researchers explored several non-invasive models, such as the portal vein diameter (7), Child-Pugh score (8) and model for end-stage liver disease (MELD) score (9), to predict esophageal variceal rebleeding in cirrhotic patients. However, the predictive performance of these non-invasive tool is still controversial.

Radiomics is an emerging field that extracts innumerable quantitative medical features from imaging into high-dimensional data using many image characterization algorithms (10, 11). It has great diagnostic and prognostic value in many fields of non-neoplastic liver lesions, such as liver fibrosis (12, 13), hypertension and EV bleeding (14, 15). Most radiomics-related studies on EV bleeding have focused mainly on the prediction of the severity of EV and the presence of EV bleeding. However, the prediction of esophageal variceal rebleeding based on radiomics has not yet been reported.

In this study, we constructed and validated a radiomics score (Rad-score) derived from radiomics features of the liver and spleen in cirrhotic patients to predict the risk of rebleeding. Moreover, the Rad-score was used to stratify patients into high-, intermediate- and low-risk groups.

Materials and Methods

This retrospective study was approved by the institutional review broad, and the requirement for written informed consent was waived.

Patients

In this study, data from 173 patients diagnosed with cirrhosis between January 2011 and December 2019 were retrospectively analyzed. The patient inclusion criteria were as follows: (1) patients who had recovered from a first episode of EV bleeding and there was no bleeding for at least 5 consecutive days; (2) abdominal computed tomography (CT) scan and HVPG measurement were performed before endoscopic variceal ligation within 2 weeks; and (3) at least 1 year of follow-up after endoscopic therapies. The patient exclusion criteria were as follows: (1) previous therapy including splenectomy, endoscopic variceal ligation, tissue adhesive injection, or usage of non-selective beta blocker to prevent rebleeding; (2) confirmed to have hepatocellular carcinoma based on a histologic examination of the liver; (3) non-sinusoidal portal hypertension (e.g., hepatic cavernoma, Budd-Chiari syndrome); and (4) no contrast-enhanced CT images and HVPG measurement. The recruitment process is shown in Supplementary Figure 1.

Definitions of Rebleeding and Therapy

The endpoint of the study was EV rebleeding during the 1-year follow-up. EV rebleeding is defined as the occurrence of new esophageal variceal bleeding after a period of 24 h or more from the 24-h point of stable vital signs and hematocrit/hemoglobin following the first episode of EV bleeding (4, 16).

The first episode of EV bleeding was controlled by measures including fluid resuscitation and medication administration (somatostatin and proton pump inhibitors). After recovering from a first episode of EV bleeding, all patients received secondary prevention of EV bleeding according to sixth Baveno Consensus (Baveno VI) (3), namely, the combination of endoscopic variceal ligation and non-selective beta blocker. Moreover, all hepatitis B-related cirrhotic patients received antiviral therapy.

Endoscopic examination was performed by experienced endoscopists. During the examination, the form, location and bleeding signs of varices were noted, and the size of the varices was classified as small, medium or large corresponding to <30, 30–60, or >60%, respectively, of the maximum theoretical size (17).

CT Image Acquisition and Analysis

Contrast-enhanced CT scans were performed using a 320-detector CT scanner (Aquilion ONE, TOSHIBA) and a 64-detector CT scanner (Discovery, GE Healthcare). Non-enhanced CT scans were first acquired, followed by three post-contrast CT scans in three phases: arterial, portal vein and delayed. Arterial phase scanning started ~20–30 s after injection, portal phase scanning was started 30–40 s after the beginning of the arterial phase, and delayed phase scanning was started 40–60 s after the beginning of the portal phase scanning. The following parameters were used: tube voltage, 120 kV; tube current, 150–600 mAs; 80 × 0.5 mm or 64 × 0.625 mm detector collimation, matrix, 512 × 512; slice thickness, 5 mm; and pitch, 1.388 or 0.984. All patients received an intravenous, non-ionic contrast medium (iodine concentration, 370 mg/mL; volume, 1.5–2.0 mL/kg of body weight; Omnipaque 350, GE Healthcare, Shanghai, China) at a rate of 3–5 mL/s. Two imaging-based indexes including diameters of portal vein and spleen vein were assessed.

Image Segmentation and Radiomics Feature Extraction

Regions of interest (ROIs) were drawn around the whole liver and spleen slice-by-slice using 3D-slicer software version 4.10.2 (Boston, USA) by a radiologist (Reader 1, Z.X.Y.) with 12 years of working experience in abdominal imaging. In ROIs of the liver and spleen, each ROI was as close as possible to the margin but excluded large vascular structures and artifacts to avoid adjacent organs such as the gallbladder, intestine, stomach, kidney and mesentery (Figure 1).

FIGURE 1
www.frontiersin.org

Figure 1. Selection and three-dimensional reconstruction of Regions of interest (ROIs) in the liver and spleen. Delineation of the liver (A) and spleen (B) as ROIs and three-dimensional reconstruction of ROIs by using 3D-slicer software for the extraction of radiomics features.

After image segmentation, we used Python 3.8.3 based on pyradiomics (version 3.0; https://pyradiomics.readthedocs.io/en/latest/index.html) for feature extraction. A total of 2,264 features were extracted from the liver and spleen ROIs (1,132 from each organ). Furthermore, the images of 173 patients were segmented by another radiologist (Reader 2, W.X.M.) who specialized in abdominal imaging and had 26 years of working experience to evaluate reproducibility. Reader 1 outlined the ROIs again after 1 month to minimize recall bias. The interobserver reproducibility and intraobserver reproducibility of all extracted features were evaluated by intra/interclass correlation coefficients (ICCs). Features with ICCs > 0.75 were considered to have good reproducibility. In the reproducibility analysis, a total of 1,882 features (953 from liver images and 929 features from spleen images) were found to be sufficiently reproducible and stable (ICCs > 0.75).

Radiomics Feature Selection and Rad-Score Calculation

The radiomics workflow is presented in Figure 2. The training cohort was used for feature selection and model building, while the validation cohort was used to test model performance. To select the best features and avoid overfitting from the training cohort, we used the least absolute shrinkage and selection operator (LASSO) method and conducted 100 iterations of 10-fold cross-validations to develop a Lasso Cox regression model (18, 19). The coefficients of some features were decreased to zero by adding penalty terms through Lasso Cox regression, and the features with non-zero coefficients were then selected. Moreover, the optimal tuning parameter (λ) is the value for which the partial likelihood deviance is the minimum criterion. The significant features were weighted with their coefficients and summed to form the Rad-score (Rad-score = coefficient 1 × feature 1+ coefficient 2 × feature 2…) (20). The Rad-scoreLiver, Rad-scoreSpleen, and Rad-scoreLiver−Spleen were calculated by a linear combination of the selected features from the liver, spleen and a combination of both organs that were weighted by their own coefficients in the LASSO Cox regression model.

FIGURE 2
www.frontiersin.org

Figure 2. Flowchart for the radiomics analysis. (A) Regions of interest (ROIs) of the liver and spleen were segmented manually on all axial slices. (B) Three-dimensional reconstruction, texture analysis and feature extraction of ROIs. (C) For feature selection, the least absolute shrinkage and selection operator (LASSO) Cox method was used. A radiomics score was generated by a linear combination of selected features. (D) Calibration curves and decision curve analysis (DCA) were utilized to evaluate the Rad-scoreLiver−Spleen.

Assessment and Performance of the Rad-Score

Harrell's concordance index (C-index) and the hazard ratio (HR) were calculated to evaluate the predictive accuracy of the Rad-score. Kaplan–Meier survival analysis and the log-rank test were used to evaluate the stratification ability of each model. In addition, calibration curves were generated to assess the calibration of the Rad-score. Decision curve analysis (DCA) was performed to analyze the clinical usefulness of the Rad-score by measuring the net benefit at different threshold probabilities.

Statistical Analysis

Statistical analysis was conducted with R software (version 3.6.1; http://www.r-project.org). The following R packages were used: glmnet, for running LASSO Cox; psych, for calculating ICCs; survival, for building the Cox proportional risk model and drawing Kaplan–Meier curves; hmisc, for calculating the C-index; rms, for generating calibration curves; stdca, for plotting DCA results; stats, for Mann–Whitney U and chi-square-tests; survcomp, for comparison of different C-indexes; and SurvProb, for predicting EV rebleeding probabilities. All statistical tests were two-sided, and p-values < 0.05 were considered significant.

Results

Study Patients

A total of 173 patients were divided into a training set and a validation set at a ratio of 7:3 with a random sampling method; 121 patients constituted the training cohort, and the other 52 constituted the validation cohort. There was no significant difference in clinical characteristics between the two cohorts (p = 0.212–0.868; Table 1). During the follow-up periods, rebleeding occurred in 39 of 173 patients (22.5%) within 1 year.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of patients in the training and validation cohorts.

Radiomics Feature Extraction, Selection, and Rad-Score Calculation

After extracting features from ROIs, we obtained 7, 6, and 10 features with non-zero coefficients as the predictive radiomics features for the liver, spleen and both organs, respectively. The Rad-score formulas are follows:

1. Rad-scoreLiver = −0.690 (Wavelet-LHH _Glrlm _RunEntropy)

−0.395 (original_glcm_JointAverage)

−0.250 (wavelet.HL_firstorder_Skewness)

+0.036(wavelet.LH_glcm_ClusterProminence)

+0.218 (log-sigma-4-0-mm-3D_glszm_LargeAreaEmphasis)

+0.535 (wavelet.HL_glszm_GrayLevelVariance)

+0.707 (wavelet-HHL_glcm_InverseVariance)

2. Rad-scoreSpleen = −0.612 (log-sigma-2-0-Mm-3D_Firstorder_TotalEnergy)

−0.396 (wavelet-LHH_glszm_LargeAreaEmphasis)

+0.045 (wavelet-HLL _glszm_ZoneEntropy)

+0.162 (log-sigma-5-0-mm-3D_gldm_DependenceVariance)

+0.443 (wavelet-HL_firstorder_RootMeanSquared)

+0.540 (wavelet-HHL_glrlm_HighGrayLevelRunEmphasis)

3. Rad-scoreLiver−Spleen = −0.529 (Wavelet-LHH_Glrlm_RunEntropy)

−0.318 (wavelet-LHH_glcm_JointEnergy)

−0.214 (log-sigma-2-0-mm-3D_firstorder_TotalEnergy)

−0.151 (wavelet.HL_firstorder_Skewness)

−0.093 (wavelet-LHH_glszm_LargeAreaEmphasis)

+0.102 (log-sigma-5-0-mm-3D_glrlm_RunVariance)

+0.139 (log-sigma-4-0-mm-3D_glszm_LargeAreaEmphasis)

+0.158 (wavelet-HHL_glrlm_HighGrayLevelRunEmphasis)

+0.432 (wavelet.HL_glszm_GrayLevelVariance)

+0.548 (wavelet-HHL_glcm_InverseVariance)

Figure 3 represents the process of features selected with non-zero coefficients in the LASSO Cox regression model in the training cohort. The ICCs of the selected features are also described in Supplementary Table 1.

FIGURE 3
www.frontiersin.org

Figure 3. Feature selection using the least absolute shrinkage and selection operator (LASSO) Cox regression model. The partial likelihood deviance was plotted vs. log (λ). The tuning parameter (λ) was chosen in the LASSO Cox model via the minimum criteria. Dotted vertical lines were drawn at both the optimal and minimum values from left to right by using the minimum criteria and 1 standard error of the minimum criteria. (A) In the liver group, we examined the coefficients of the 953 radiomics features to identify 7 potential predictors. A λ-value of 0.0714, with log (λ), −2.6391, was chosen using 10-fold cross-validation. (B) In the spleen group, we examined the coefficients of the 929 radiomics features to identify 6 potential predictors. A λ-value of 0.0554, with log (λ) −2.8915, was chosen. (C) In the combined liver and spleen group, we examined the coefficients of the 1882 radiomics features to identify 10 potential predictors. A λ-value of 0.0675, with log (λ) −2.6956, was chosen.

Univariate and Multivariate Cox Regression Analysis of Rad-Score and Clinical Characteristics

Predictive factors for EV rebleeding are summarized in Table 2. In the univariate analysis, portal vein diameter [HR = 1.114; 95% confidence interval (CI) = 1.010–1.220; p = 0.002], Rad-scoreLiver (HR = 1.692; 95% CI = 1.132–2.262; p < 0.001), Rad-scoreSpleen (HR = 1.368; 95% CI = 1.019–1.741; p < 0.001) and Rad-scoreLiver−Spleen (HR = 3.025; 95% CI = 2.029–3.961; p < 0.001) showed a significant association with EV rebleeding. In the multivariate analysis, Rad-scoreLiver (HR = 1.355; 95% CI = 1.101–1.503; p =0.008), Rad-scoreSpleen (HR = 1.148; 95% CI = 1.007–1.396; p =0.034) and Rad-scoreLiver−Spleen (HR = 2.682; 95% CI = 1.793–3.512; p < 0.001) were independent risk factors for EV rebleeding.

TABLE 2
www.frontiersin.org

Table 2. HR analysis of clinical characteristics and Rad-scores for predicting EV rebleeding.

Performance of the Rad-Score for EV Rebleeding

To compare the predictive performance of Rad-scoreLiver, Rad-scoreSpleen and Rad-scoreLiver−Spleen for EV rebleeding, the C-index was calculated. Rad-scoreLiver−Spleen showed significantly better performance than Rad-scoreLiver and Rad-scoreSpleen, yielding a C-index of 0.853 (95% CI = 0.776–0.904) in the training cohort and 0.822 (95% CI = 0.749–0.875) in the validation cohort (Table 3). The calibration curves of the Rad-scoreLiver−Spleen at 3, 6, 9, and 12 months showed that the predicted probability was very close to the actual probability (Figures 4A,B). DCA showed that the Rad-scoreLiver−Spleen yielded more clinical net benefit under almost all threshold probabilities, indicating that the Rad-scoreLiver−Spleen is more practical than the Rad-scoreLiver and Rad-scoreSpleen for predicting esophageal variceal rebleeding in cirrhotic patients (Figures 4C,D).

TABLE 3
www.frontiersin.org

Table 3. C-indexes of different models.

FIGURE 4
www.frontiersin.org

Figure 4. Calibration curves and decision curve analysis of the Rad-scoreLiver−Spleen. Calibration curves of the Rad-scoreLiver−Spleen demonstrate its predictive performance for rebleeding at 3, 6, 9, and 12 months in the training cohort (A) and validation cohort (B). Decision curve analysis was performed to compare the performance of the Rad-scoreLiver−Spleen, Rad-scoreLiver and Rad-scoreSpleen in the training cohort (C) and validation cohort (D).

We also used clinical indexes, including HVPG, Child-Pugh score, MELD score and EV size, to predict the probability of rebleeding. Compared with the clinical indexes, the Rad-scoreLiver−Spleen exhibited a higher C-index in the training and validation cohorts (Table 3).

Risk Stratification for Predicting EV Rebleeding According to Rad-ScoreLiver-Spleen

Based on the cutoff values of the Rad-scoreLiver−Spleen determined by X-tile software (21), all patients were divided into 3 risk groups to predict rebleeding in the training (low-risk, −0.03–0.30; intermediate-risk, 0.31–0.61; high-risk, 0.62–0.9; Figures 5A,B) and validation (low-risk, −0.03–0.30; intermediate-risk, 0.31–0.61; high-risk, 0.62–0.9; Figures 5D,E) cohorts. The 12-month rebleeding probabilities among the 3 risk groups in the training cohort were 0.090, 0.202, and 0.407. Likewise, significant differences were observed in the validation cohort (12-month rebleeding probability: 0.097 for the low-risk group, 0.218 for the intermediate-risk group, and 0.436 for the high-risk group; Table 4). Kaplan–Meier curves showed that the cumulative incidences of rebleeding in the training and validation cohorts were accurately differentiated by the risk stratification system (Figures 5C,F).

FIGURE 5
www.frontiersin.org

Figure 5. X-tile analysis of the total Rad-scoreLiver−Spleen and survival curves stratified by the score calculated by the Rad-scoreLiver−Spleen in the training (A–C) and validation (D–F) cohorts.

TABLE 4
www.frontiersin.org

Table 4. Rebleeding rates and probabilities according to the risk stratification.

Subgroup Analysis for Predicting EV Rebleeding in Hepatitis B Virus and Non-HBV Cohorts

For the subgroup analysis in the training cohort, the Rad-scoreLiver−Spleen in the HBV group had significantly better performance than that in the non-HBV group (C-index, 0.903 vs. C-index, 0.791; P < 0.001). Significant differences were also observed in the validation cohort (C-index, 0.884 vs. C-index, 0.781; P < 0.001, Table 5).

TABLE 5
www.frontiersin.org

Table 5. Subgroup analysis of the C-indexes of the HBV and non-HBV cohorts.

Discussion

Non-invasive tools for predicting EV rebleeding and risk stratification have been highlighted in recent years. The present study developed a Rad-score extracted from features of both the liver and spleen to predict EV rebleeding. Our results showed that Rad-scoreLiver−Spleen was an independent significant predictive factor and achieved great predictive performance. In addition, Rad-scoreLiver−Spleen could stratify patients into low-, intermediate- and high-risk groups for predicting rebleeding probability. Thus, the Rad-scoreLiver−Spleen might be a promising tool to predict EV rebleeding in cirrhotic patients.

Based on the results of the LASSO Cox regression analysis, a total of 10 potential radiomics features were selected to calculate the Rad-scoreLiver−Spleen. Among these, run entropy, run variance and high gray level run emphasis measured the randomness and variance in the distribution of run lengths or higher gray-level values. Consistent with previous studies (14), a higher absolute value of high gray level run emphasis increased the possibility of EV bleeding. Joint energy and inverse variance were measures of homogeneous patterns in the image; if the image texture was relatively uniform and changed slowly between different regions, the inverse variance was increased. These features had a proper ratio for calculation of the Rad-score that could avoid overfitting and mainly reflected the texture complexity of the liver and spleen (15).

Our study revealed that Rad-scoreLiver, Rad-scoreSpleen, and Rad-scoreLiver−Spleen were independent risk factors for EV rebleeding, suggesting that radiomics features of the liver and spleen were closely related to variceal bleeding. It was consistent with previous studies reporting that radiomics has a potential role in diagnosing portal hypertension and EV bleeding (14, 15, 22, 23). This finding could be explained by the hepatic-related factors and splenomegaly contributed to the rise of portal pressure in cirrhotic patients (24, 25). EV size and HVPG were not independent predictors in our study, which might be explained by the fact that non-selective beta blocker treatment can decrease portal blood flow and variceal pressure, leading to a change in hemodynamics.

Endoscopy and HVPG measurement which were reported to be predictors for EV rebleeding (4, 5, 26), are highly limited by their invasiveness and are therefore not suitable for dynamic monitoring. In contrast, Rad-scoreLiver−Spleen is non-invasive and reproducible, it can extract quantitative features that reflect information related to all directions of the complex spatial structure of organs that are invisible to the human eye. Clinical physicians need only to upload CT images and select the ROI of the liver and spleen to perform the radiomics analysis and help to assess the risk of rebleeding in cirrhotic patients.

Cirrhotic patients usually undergo endoscopy every 3–6 months (16) after successful eradication of the varices. In order to reduce or avoid endoscopy examinations, it is of great significance for physicians to determine appropriate candidates for endoscopy according to risk stratification. In this study, Rad-scoreLiver−Spleen divided all patients into low-, intermediate- and high-risk groups (3). Patients in the low-risk group could avoid endoscopy, while for patients in the high-risk group, endoscopy was performed as soon as possible to prevent rebleeding. For patients in the intermediate-risk group, regular follow-up should be carried out every 3–6 months until the Rad-scoreLiver−Spleen reached the standard of the high-risk group.

In our study, 60.1% of patients had been infected with HBV, which remains the primary cause of cirrhosis in most Asian nations (27). Our results showed that Rad-scoreLiver−Spleen had a significantly better performance in the HBV group than that in the non-HBV group, indicating that Rad-scoreLiver−Spleen was particularly more suitable for the HBV population than for the non-HBV population.

There are several limitations to this study. First, this study was a single-center, retrospective analysis and subjected to the inherent limitations of such investigations. A multicenter, prospective study with a larger data set is needed. Second, we lacked hemodynamic data of the left gastric vein, portal vein, spleen vein, liver stiffness and spleen stiffness by ultrasound and transient elastography, which have proven to be good predictors of the degree of cirrhosis and the development of EV bleeding (28, 29). A future study comparing radiomics with other radiologic methods is needed.

Conclusions

Our findings demonstrated that the Rad-scoreLiver−Spleen could be used to predict the probability of EV rebleeding and stratify cirrhotic patients accordingly. The Rad-scoreLiver−Spleen might serve as a useful tool for clinicians involved in therapeutic decision-making and individualized patient counseling.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by Biomedical Research Ethic Committee of Shandong Provincial Hospital. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements. Written informed consent was not obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

DM and XZ devised the experiment. DM and YW developed and organized the paper. XF and BK designed the tables and figures. XW and JQ performed the data analysis. XZ and QZ participated in the revision of the manuscript. DM and XZ wrote the original draft. All authors read and approved the final manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (81770607).

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

Conflict of Interest

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

Publisher's Note

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.

Supplementary Material

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

Supplementary Figure 1. Flowchart of the study population.

Supplementary Table 1. Intra/interclass correlation coefficients (ICCs) of the radiomics model.

References

1. Garcia-Tsao G, Abraldes JG, Berzigotti A, Bosch J. Portal hypertensive bleeding in cirrhosis: risk stratification, diagnosis, and management: 2016 practice guidance by the American Association for the study of liver diseases. Hepatology. (2017) 65:310–35. doi: 10.1002/hep.28906

PubMed Abstract | CrossRef Full Text

2. Garcia-Tsao G, Bosch J. Management of varices and variceal hemorrhage in cirrhosis. N Engl J Med. (2010) 362:823–32. doi: 10.1056/NEJMra0901512

PubMed Abstract | CrossRef Full Text | Google Scholar

3. de Franchis R, Baveno VIF. Expanding consensus in portal hypertension: report of the Baveno VI consensus workshop: stratifying risk and individualizing care for portal hypertension. J Hepatol. (2015) 63:743–52. doi: 10.1016/j.jhep.2015.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Garcia-Tsao G, Sanyal AJ, Grace ND, Carey W. Practice Guidelines Committee of the American Association for the Study of Liver D, Practice Parameters Committee of the American College of G. Prevention and management of gastroesophageal varices and variceal hemorrhage in cirrhosis Hepatology. (2007) 46:922–38. doi: 10.1002/hep.21907

PubMed Abstract | CrossRef Full Text

5. Ripoll C, Groszmann R, Garcia-Tsao G, Grace N, Burroughs A, Planas R, et al. Hepatic venous pressure gradient predicts clinical decompensation in patients with compensated cirrhosis. Gastroenterology. (2007) 133:481–8. doi: 10.1053/j.gastro.2007.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Bosch J, Abraldes JG, Berzigotti A, Garcia-Pagan JC. The clinical use of HVPG measurements in chronic liver disease. Nat Rev Gastroenterol Hepatol. (2009) 6:573–82. doi: 10.1038/nrgastro.2009.149

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Moubarak E, Bouvier A, Boursier J, Lebigot J, Ridereau-Zins C, Thouveny F, et al. Portosystemic collateral vessels in liver cirrhosis: a three-dimensional MDCT pictorial review. Abdom Imaging. (2012) 37:746–66. doi: 10.1007/s00261-011-9811-0

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zardi EM, Di Matteo FM, Pacella CM, Sanyal AJ. Invasive and non-invasive techniques for detecting portal hypertension and predicting variceal bleeding in cirrhosis: a review. Ann Med. (2014) 46:8–17. doi: 10.3109/07853890.2013.857831

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Xu XD, Dai JJ, Qian JQ, Pin X, Wang WJ. New index to predict esophageal variceal bleeding in cirrhotic patients. World J Gastroenterol. (2014) 20:6989–94. doi: 10.3748/wjg.v20.i22.6989

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. Nat Rev Clin Oncol. (2017) 14:749–62. doi: 10.1038/nrclinonc.2017.141

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Aerts HJ, Velazquez ER, Leijenaar RT, Parmar C, Grossmann P, Carvalho S, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun. (2014) 5:4006. doi: 10.1038/ncomms5644

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Park HJ, Lee SS, Park B, Yun J, Sung YS, Shim WH, et al. Radiomics analysis of gadoxetic acid-enhanced MRI for staging liver fibrosis. Radiology. (2019) 292:269. doi: 10.1148/radiol.2019194012

PubMed Abstract | CrossRef Full Text | Google Scholar

13. 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. (2019) 68:729–41. doi: 10.1136/gutjnl-2018-316204

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Yang JQ, Zeng R, Cao JM, Wu CQ, Chen TW, Li R, et al. Predicting gastro-oesophageal variceal bleeding in hepatitis B-related cirrhosis by CT radiomics signature. Clin Radiol. (2019) 74:976.e1–e9. doi: 10.1016/j.crad.2019.08.028

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tseng Y, Ma L, Li S, Luo T, Luo J, Zhang W, et al. Application of CT-based radiomics in predicting portal pressure and patient outcome in portal hypertension. Eur J Radiol. (2020) 126:108927. doi: 10.1016/j.ejrad.2020.108927

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tripathi D, Stanley AJ, Hayes PC, Patch D, Millson C, Mehrzad H., et al. UK guidelines on the management of variceal haemorrhage in cirrhotic patients. Gut. (2015) 64:1680–704. doi: 10.1136/gutjnl-2015-309262

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Pagliaro L., Paolo Spina G., D'Amico G., Brocchi E., Caletti G., Cosentino F. Reliability of endoscopy in the assessment of variceal features. The Italian Liver Cirrhosis Project. J Hepatol. (1987) 4:93–8.

PubMed Abstract | Google Scholar

18. Yang L, Yang J, Zhou X, Huang L, Zhao W, Wang T, et al. Development of a radiomics nomogram based on the 2D and 3D CT features to predict the survival of non-small cell lung cancer patients. Eur Radiol. (2019) 29:2196–206. doi: 10.1007/s00330-018-5770-y

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Li W, Zhang L, Tian C, Song H, Fang M, Hu C, et al. Prognostic value of computed tomography radiomics features in patients with gastric cancer following curative resection. Eur Radiol. (2019) 29:3079–89. doi: 10.1007/s00330-018-5861-9

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Royston P, Altman DG. External validation of a cox prognostic model: principles and methods. BMC Med Res Methodol. (2013) 13:33. doi: 10.1186/1471-2288-13-33

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. (2004) 10:7252–9. doi: 10.1158/1078-0432.CCR-04-0713

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Liu F, Ning Z, Liu Y, Liu D, Tian J, Luo H, et al. Development and validation of a radiomics signature for clinically significant portal hypertension in cirrhosis (CHESS1701): a prospective multicenter study. EBioMedicine. (2018) 36:151–8. doi: 10.1016/j.ebiom.2018.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Wan S, Wei Y, Zhang X, Liu X, Zhang W, He Y, et al. Multiparametric radiomics nomogram may be used for predicting the severity of esophageal varices in cirrhotic patients. Ann Transl Med. (2020) 8:186. doi: 10.21037/atm.2020.01.122

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kim BH, Chung JW, Lee CS, Jang ES, Jeong SH, Kim N, et al. Liver volume index predicts the risk of esophageal variceal hemorrhage in cirrhotic patients on propranolol prophylaxis. Korean J Intern Med. (2019) 34:1233–43. doi: 10.3904/kjim.2018.120

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhou YJ, Zheng JN, Zhou YF, Han YJ, Zou TT, Liu WY, et al. Development of a prognostic nomogram for cirrhotic patients with upper gastrointestinal bleeding. Eur J Gastroenterol Hepatol. (2017) 29:1166–73. doi: 10.1097/MEG.0000000000000943

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Li GQ, Yang B, Liu J, Wang GC, Yuan HP, Zhao JR, et al. Hepatic venous pressure gradient is a useful predictor in guiding treatment on prevention of variceal rebleeding in cirrhosis. Int J Clin Exp Med. (2015) 8:19709–16.

PubMed Abstract | Google Scholar

27. Wong MCS, Huang JLW, George J, Huang J, Leung C, Eslam M, et al. The changing epidemiology of liver diseases in the Asia-Pacific region. Nat Rev Gastroenterol Hepatol. (2019) 16:57–73. doi: 10.1038/s41575-018-0055-0

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hino S, Kakutani H, Ikeda K, Uchiyama Y, Sumiyama K, Kuramochi A, et al. Hemodynamic assessment of the left gastric vein in patients with esophageal varices with color Doppler EUS: factors affecting development of esophageal varices. Gastrointest Endosc. (2002) 55:512–7. doi: 10.1067/mge.2002.122333

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Sandrin L, Fourquet B, Hasquenoph JM, Yon S, Fournier C, Mal F, et al. Transient elastography: a new noninvasive method for assessment of hepatic fibrosis. Ultrasound Med Biol. (2003) 29:1705–13. doi: 10.1016/j.ultrasmedbio.2003.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: portal hypertension, non-invasive, computed tomography, radiomics, esophageal variceal rebleeding

Citation: Meng D, Wei Y, Feng X, Kang B, Wang X, Qi J, Zhao X and Zhu Q (2021) CT-Based Radiomics Score Can Accurately Predict Esophageal Variceal Rebleeding in Cirrhotic Patients. Front. Med. 8:745931. doi: 10.3389/fmed.2021.745931

Received: 22 July 2021; Accepted: 15 October 2021;
Published: 04 November 2021.

Edited by:

Pradeep Kumar Shukla, University of Tennessee Health Science Center (UTHSC), United States

Reviewed by:

Shivendra Vikram Singh, St. Jude Children's Research Hospital, United States
Pankaj Taneja, Sharda University, India

Copyright © 2021 Meng, Wei, Feng, Kang, Wang, Qi, Zhao and Zhu. 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: Xinya Zhao, zhaoxinya2000@126.com; Qiang Zhu, zhuqiang@sdu.edu.cn

Download